We observed several swarms of repeating low-frequency (1–5 Hz) seismic events during a 3 week period in May–June 2010, near the summit of Mount Rainier, Washington, USA, that likely were a result of stick–slip motion at the base of alpine glaciers. The dominant set of repeating events (‘multiplets’) featured >4000 individual events and did not exhibit daytime variations in recurrence interval or amplitude. Volcanoes and glaciers around the world are known to produce seismic signals with great variability in both frequency content and size. The low-frequency character and periodic recurrence of the Mount Rainier multiplets mimic long-period seismicity often seen at volcanoes, particularly during periods of unrest. However, their near-surface location, lack of common spectral peaks across the recording network, rapid attenuation of amplitudes with distance, and temporal correlation with weather systems all indicate that ice-related source mechanisms are the most likely explanation. We interpret the low-frequency character of these multiplets to be the result of trapping of seismic energy under glacial ice as it propagates through the highly heterogeneous and attenuating volcanic material. The Mount Rainier multiplet sequences underscore the difficulties in differentiating low-frequency signals due to glacial processes from those caused by volcanic processes on glacier-clad volcanoes.

Many of Greenland’s marine-terminating outlet glaciers have undergone rapid retreat in the past decade, accompanied by accelerated flow and dynamic thinning. Superimposed on this pattern of retreat, these glaciers undergo seasonal variations in terminus position, corresponding roughly to wintertime advance and summertime retreat. We compiled near-daily time series of terminus position for five of Greenland’s largest outlet glaciers (Daugaard Jensen, Kangerdlugssuaq and Helheim glaciers in East Greenland, and Jakobshavn Isbræ and Rink Isbræ in West Greenland) using Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery. There are spatial differences in the timing of the onset of seasonal retreat among all the glaciers in our study, as well as variability in terminus behavior for individual glaciers from year to year. We examine whether this spatial and temporal variability is linked to above-freezing air temperatures or high sea surface temperatures, but find no simple relationship. Instead, we hypothesize that terminus geometry (ice thickness, subglacial topography, fjord bathymetry) exerts an important control on the response of marine-terminating glaciers to climate perturbations. Models for predicting outlet glacier response to climate change need to include this complex interaction between geometry and environmental forcing.

A new 10 year surface mass balance (SMB) record of Hurd and Johnsons Glaciers, Livingston Island, Antarctica, is presented and compared with earlier estimates on the basis of local and regional meteorological conditions and trends. Since Johnsons is a tidewater glacier, we also include a calving flux calculation to estimate its total mass balance. The average annual SMB over the 10 year observation period 2002–11 is −0.15 ± 0.10 m w.e. for Hurd Glacier and 0.05 ± 0.10 m w.e. for Johnsons Glacier. Adding the calving losses to the latter results in a total mass balance of −0.09 ± 0.10 m w.e. There has been a deceleration of the mass losses of these glaciers from 1957–2000 to 2002–11, which have nearly halved for both glaciers. We attribute this decrease in the mass losses to a combination of increased accumulation in the region and decreased melt. The increased accumulation is attributed to larger precipitation associated with the recent deepening of the circumpolar pressure trough, while the melt decrease is associated with lower summer surface temperatures during the past decade.

The subglacial water system of lower Whillans Ice Stream on the Siple Coast, West Antarctica, contains numerous connected subglacial lakes in three hydrological basins (northern, central and southern). We use Ice, Cloud and land Elevation Satellite (ICESat) data to derive estimates of lake volume change and regional thickness changes. By combining these results with a water budget model, we show that a uniform, localized thickness increase perturbed the hydropotential, resulting in a change in course of a major flow path within the system in 2005. Water originating from upper Whillans and Kamb Ice Streams that previously supplied the southern basin became diverted toward Subglacial Lake Whillans (SLW). This diversion led to a tenfold filling rate increase of SLW. Our observation suggests that water piracy may be common in the Siple Coast region, where the gentle basal relief makes the basal hydropotential particularly sensitive to small changes in ice thickness. Given the previously inferred connections between water piracy and ice-stream slowdown elsewhere in the region, the subtle and complex nature of this system presents new challenges for numerical models.

We present an analysis of temporal trends in ∼55 000 avalanches recorded between 1946 and 2010 in the French Alps and two north/south subregions. First, Bayesian hierarchical modelling is used to isolate low-, intermediate- and high-frequency trends in the mean avalanche occurrence and runout altitude per year/winter. Variables are then combined to investigate their correlation and the recent evolution of large avalanches. Comparisons are also made to climatic and flow regime covariates. The results are important for risk assessment, and the development of new high-altitude climate proxies. At the entire French Alps scale, a major change-point exists in ∼1978 at the heart of a 10 year period of high occurrences and low runout altitudes corresponding to colder and snowier winters. The differences between this change-point and the beginning/end of the study period are 0.1 avalanche occurrences per winter and per path and 55 m in runout altitude. Trends before/after are well correlated, leading to enhanced minimal altitudes for large avalanches at this time. A marked upslope retreat (80 m for the 10 year return period runout altitude) accompanied by a 12% decrease in the proportion of powder snow avalanches has occurred since then, interrupted from about 2000. The snow-depth and temperature control on these patterns seems significant (R = 0.4–0.6), but is stronger at high frequencies for occurrences, and at lower frequencies for runout altitudes. Occurrences between the northern and southern French Alps are partially coupled (R∼0.4, higher at low frequencies). In the north, the main change-point was an earlier shift in ∼1977, and winter snow depth seems to be the main control parameter. In the south, the main change-point occurred later, ∼1979–84, was more gradual, and trends are more strongly correlated with winter temperature.

A new densification model, which simulates the effect of impurities on the densification of polar firn, is presented. The classical densification models of Herron and Langway (1980) and Pimienta and Barnola (Barnola and others, 1991) are modified by assuming that the activation energy for deformation is reduced by the impurities. Motivated by recent observations, the impurity effect is formulated on an empirical basis using the seasonally varying Ca2+ ion concentration. Excellent agreement between simulated and measured high-resolution density profiles confirms the new approach. The same parameterization applies for Greenland and Antarctica despite the one order of magnitude difference in impurity concentration. The new models allow us, for the first time, to simulate the density layering in firn down to the firn–ice transition. Our results emphasize the importance of impurities and density layering for the air entrapment and for dating gas records of deep ice cores, in particular for glacial climate conditions where the impurity concentrations are 10–100-fold higher than in modern firn.

We investigated the large-scale (10–1000 m) and small-scale (mm–cm) variations in size, number and arrangement of air bubbles in the EPICA Dronning Maud Land (EDML) (Antarctica) ice core, down to the end of the bubble/hydrate transition (BHT) zone. On the large scale, the bubble number density shows a general correlation with the palaeo-temperature proxy, δ18O, and the dust concentration, which means that in Holocene ice there are fewer bubbles than in glacial ice. Small-scale variations in bubble number and size were identified and compared. Above the BHT zone there exists a strong anticorrelation between bubble number density and mean bubble size. In glacial ice, layers of high number density and small bubble size are linked with layers with high impurity content, identified as cloudy bands. Therefore, we regard impurities as a controlling factor for the formation and distribution of bubbles in glacial ice. The anticorrelation inverts in the middle of the BHT zone. In the lower part of the BHT zone, bubble-free layers exist that are also associated with cloudy bands. The high contrast in bubble number density in glacial ice, induced by the impurities, indicates a much more pronounced layering in glacial firn than in modern firn.

The values of the Voellmy friction parameters of 735 historical avalanches that have occurred along 26 paths in the Chamonix valley, France, since 1958 are back-analyzed with a depth-averaged hydraulic model, including sub-models for erosion, entrainment and deposition. For each path, the longitudinal and crosswise topographic profiles were derived from a high-resolution digital elevation model acquired by laser scanning. The initial snow depth and snow cohesion, as well as various physical properties of snow, were computed from numerical simulations of the detailed snowpack model Crocus fed by the SAFRAN meteorological analysis. For each event, the full ranges of the two friction parameters were scanned and the pairs of friction parameters for which the run-out altitude is found close enough to the observed one (with an uncertainty of ±5 m) were retained. Statistical class analysis was used to investigate the correlation between the obtained friction coefficients and the snow physical properties. No evident trend with the snow parameters was found for the inertial friction coefficient. For the static friction coefficient, an increasing trend with temperature and density was observed, as well as a decreasing trend with liquid water content and initial snow depth. Although modeling assumptions and limitations regarding data and the calibration procedure should be kept in mind, these trends are worth noting, allowing avalanche simulations to be refined to take into account prevailing weather and snow conditions.

Monitoring of ice-shelf and sub-ice-shelf ocean temperatures represents an important component in understanding ice-sheet stability. Continuous monitoring is challenging due to difficult surface access, difficulties in penetrating the ice shelf, and the need for long-term operation of non-recoverable sensors. We aim to develop rapid lightweight drilling and near-continuous fiber-optic temperature-monitoring methods to meet these challenges. During November 2011, two instrumented moorings were installed within and below the McMurdo Ice Shelf (a sub-region of the Ross Ice Shelf, Antarctica) at Windless Bight. We used a combination of ice coring for the upper portion of each shelf borehole and hot-point drilling for penetration into the ocean. The boreholes provided temporary access to the ice-shelf cavity, into which distributed temperature sensing (DTS) fiber-optic cables and conventional pressure/temperature transducers were installed. The DTS moorings provided near-continuous (in time and depth) observations of ice and ocean temperatures to a depth of almost 800 m beneath the ice-shelf surface. Data received document the presence of near-freezing water throughout the cavity from November through January, followed by an influx of warmer water reaching ∼150 m beneath the ice-shelf base during February and March. The observations demonstrate prospects for achieving much higher spatial sampling of temperature than more conventional oceanographic moorings.

Repeated Global Navigation Satellite Systems (GNSS) observations were carried out at 50 surface markers in the Vostok Subglacial Lake (East Antarctica) region between 2001 and 2011. The horizontal ice flow velocity vectors were derived with accuracies of 1 cm a−1 and 0.5°, representing the first reliable information on ice flow kinematics in the northern part of the lake. Within the lake area, ice flow velocities do not exceed 2 m a−1. The ice flow azimuth is southeast in the southern part of the lake and turns gradually to east-northeast in the northern part. In the northern part, as the ice flow enters the lake at the western shore, the velocity decreases towards the central lake axis, then increases slightly past the central axis. In the southern part, a continued acceleration is observed from the central lake axis across the downstream grounding line. Based on the observed flow velocity vectors and ice thickness data, mean surface accumulation rates are inferred for four surface segments between Ridge B and Vostok Subglacial Lake and show a steady increase towards the north.

In the simulation of skiing the force between ski and snow is a decisive factor. We decompose the reaction force into a penetration force normal to the snow surface, a shear force and friction. Two portable measurement devices were developed to study the penetration and shear forces for compacted snow on groomed ski slopes. The penetration force was assessed by measuring the penetration depth of a ski-tool loaded normal to the snow surface. For the shear force the tangential load was measured when the snow began to fail. Overall 236 penetration and 108 shear experiments were conducted on different types of snow. The penetration force was proportional to the volume of snow displaced by the ski-tool. The failure shear force was proportional to the penetration depth multiplied by the length of the tool. The constants of proportionality, HV and Sf, are material parameters of snow. The snow hardness, HV, varied between 0.04 and 90 N mm–3 and the failure shear stress, Sf, between 0.04 and 0.40 N mm–2. In another investigation, skiing turns were simulated using the presented snow reaction forces. Maximum deviations between computed and real trajectories were <1 % of the overall length of the runs.

Estimating regional-scale rates of glacier mass change from repeat airborne and satellite altimetry requires assumptions about the density of the material being added to, or removed from, the glacier surface. This may vary in both space and time, especially in a period of warming climate such as has occurred in Arctic Canada since 2005. Here we compare firn properties measured on 20 shallow cores recovered from elevations of 1400–1900 m a.s.l. on Devon Ice Cap in 2012 with equivalent measurements from 34 cores recovered from the same locations between 2004 and 2011. The average density of the top 2.5 m of the firn column has increased by 13–80%. As a result, measured rates of density increase in the top 1.25 mw.e. of the firn layer have likely thinned this layer at rates of 0.021–0.168 m a−1 over the period 2004–12. These thinning rates are significant relative to local rates of surface height change measured by repeat airborne laser altimetry between 2004 and 2012 (−0.077 to −0.252 m a−1) and indicate that, during periods of climate warming or cooling, firn densification rate changes must be taken into account when estimating mass change rates from measured height changes.

An attempt is made to obtain and quantify the mechanical properties of two common types of seasonal snow on the ground. Different samples of natural snow whose metamorphism had stabilized (such as would remain on a road throughout winter in a cold, snowy area) were gathered and tested using mesoscale indentation tests (metrics on the order of mm to cm). Results from the stress vs displacement curves from indentation indicated that (1) first peak strength decreased, according to a power law, with increasing indenter size and was not affected by snow average grain size, (2) plateau strength decreased with increasing indenter size, and snow compaction strength might be calculated from these data, and (3) mean energy absorption density during indentation was independent of indenter size in some size ranges, and decreased with increasing indenter size in other size ranges.

El Niño and La Niña affect global climate and atmospheric circulation to determine winter temperature and precipitation patterns. Both winter temperatures and the associated precipitation patterns have effects on mountain snow deposition and snow avalanche occurrences. Approximately 25 000 slab avalanches from 30 winters were analyzed in relation to snowfall patterns contrasted for El Niño and La Niña winters for two avalanche areas with different snow climates in British Columbia (BC), Canada. La Niña winters were shown to produce more snow, more avalanches and a higher percentage of dry avalanches than wet avalanches. The data and analysis show that the avalanche patterns depend on the altitude and snow climate. Analysis of snowfall and accident data from the Andes of Chile suggests behavior opposite to BC. El Niño winters in central Chile produced the most snow and, by inference, the most avalanches. This paper is the first to show the links of El Niño and La Niña to snow avalanche activity.

The evaluation of avalanche release depths constitutes a great challenge for risk assessment in mountainous areas. This study focuses on slab avalanches, which generally result from the rupture of a weak layer underlying a cohesive slab. We use the finite-element code Cast3M to build a mechanical model of the slab/weak-layer system, taking into account two key ingredients for the description of avalanche release: weak-layer heterogeneity and stress redistribution via slab elasticity. The system is loaded by increasing the slope angle until rupture. We first examine the cases of one single and two interacting weak spots in the weak layer, in order to validate the model. We then study the case of heterogeneous weak layers represented through Gaussian distributions of the cohesion with a spherical spatial covariance. Several simulations for different realizations of weak-layer heterogeneity are carried out and the influence of slab depth and heterogeneity correlation length on avalanche release angle distributions is analyzed. We show, in particular, a heterogeneity smoothing effect caused by slab elasticity. Finally, this mechanically based probabilistic model is coupled with extreme snowfall distributions. A sensitivity analysis of the predicted distributions enables us to determine the values of mechanical parameters that provide the best fit to field data.

In 2005 the ongoing retreat of Rhonegletscher, Switzerland, led to the formation of a proglacial lake. To investigate the influence of proglacial lake formation on the dynamics and evolution of glaciers, we measured horizontal flow velocity, vertical ice motion and water levels in boreholes with high spatial resolutions during the summer seasons of 2007–09. Annual flow speeds near the terminus increased by a factor of 2.7 from 2005/06 to 2007/08, and exceeded 20 m a−1 in 2009. The velocity increased towards the glacier front, indicating that the ice was thinning under a longitudinally stretching flow regime. Our observations show that the increase in flow speed near the terminus was due to increases in basal motion as a result of ice thinning. During summer 2009, the ice surface at the terminus moved vertically upwards by up to 4.69 m as the ice began to float on the lake. The observed ice motion can be explained by the upward bending of marginal ice and the formation of faults that cut through the entire ice thickness. We predict that if the current ice thinning continues, the basal water pressure will exceed the pressure exerted by the ice overburden, and the glacier will progressively disintegrate over an expanding area.

Studies addressing the response of glaciers to climate change have so far analyzed the effect of long-term trends in a particular set of meteorological variables only, implicitly assuming an unaltered climatic variability. Here a framework for distinguishing between year-to-year, month-to-month and day-to-day variability is proposed. Synthetically generated temperature and precipitation time series following the same long-term trend but with altered variability are then used to force an ice-dynamics model set up for Rhonegletscher, Swiss Alps. In the case of temperature, variations in the day-to-day variability are shown to have a larger effect than changes at the yearly scale, while in the case of precipitation, variability changes are assessed as having negligible impact. A first set of scenarios is used to show that compared to reference, doubling the temperature variability can reduce glacier ice volume by up to 64% within half a decade. A second set derived from the results of the European ENSEMBLES project, however, shows that such changes are expected to remain below 8% even for extreme scenarios. Although the latter results relativize the importance of the effect in the near future, the analyses indicate that at least caution is required when assuming ‘unchanged variability’.

Supraglacial lakes (SGLs) affect the dynamics of the Greenland ice sheet by storing runoff and draining episodically. We investigate the evolution of SGLs as reported in three datasets, each based on automated classification of satellite imagery. Although the datasets span the period 2001–10, there are differences in temporal sampling, and only the years 2005–07 are common. By subsampling the most populous dataset, we recommend a sampling frequency of one image per 6.5 days in order to minimize uncertainty associated with poor temporal sampling. When compared with manual classification of satellite imagery, all three datasets are found to omit a sizeable (29, 48 and 41 %) fraction of lakes and are estimated to document the average size of SGLs to within 0.78, 0.48 and 0.95 km2. We combine the datasets using a hierarchical scheme, producing a single, optimized, dataset. This combined record reports up to 67% more lakes than a single dataset. During 2005–07, the rate of SGL growth tends to follow the rate at which runoff increases in each year. In 2007, lakes drain earlier than in 2005 and 2006 and remain absent despite continued runoff. This suggests that lakes continue to act as open surface–bed conduits following drainage.

With mass loss from the Greenland ice sheet accelerating and spreading to higher latitudes, the quantification of mass discharge in the form of icebergs has recently received much scientific attention. Here we make use of very low-frequency (0.001–0.01 Hz) seismic data from three permanent broadband stations installed in the summers of 2009–10 in northwest Greenland in order to monitor local calving activity. At these frequencies, calving seismograms are dominated by a tilt signal produced by local ground flexure in response to fjord seiching generated by major iceberg calving events. A simple triggering algorithm is proposed to detect calving events from large calving fronts with potentially no user interaction. Our calving catalogue identifies spatial and temporal differences in calving activity between Jakobshavn Isbræ and glaciers in the Uummannaq district ∼200 km to the north. The Uummannaq glaciers show clear seasonal fluctuations in seiche-based calving detections as well as seiche amplitudes. In contrast, the detections at Jakobshavn Isbræ show little seasonal variation, which may be evidence for an ongoing transition to winter calving activity. The results offer further evidence that seismometers can provide efficient and inexpensive monitoring of calving fronts.