A hanging glacier at the east face of Weisshorn, Switzerland, broke off in 2005. We were able to monitor and measure surface motion and icequake activity for 25 days up to 3 days prior to the break-off. The analysis of seismic waves generated by the glacier during the rupture maturation process revealed four types of precursory signals of the imminent catastrophic rupture: (1) an increase in seismic activity within the glacier; (2) a change in the size–frequency distribution of icequake energy; (3) a modification in the structure of the waiting-time distributions between two successive icequakes; and (4) a correlation between the seismic activity and the log-periodic oscillations of the surface velocities superimposed on the global acceleration of the glacier during the rupture maturation. Analysis of the seismic activity led us to identify two regimes: a stable phase with diffuse damage and an unstable and dangerous phase characterized by a hierarchical cascade of rupture instabilities where large icequakes are triggered.

There is growing recognition of the significance of biologically active supraglacial dust (cryoconite) for glacial mass balance and ecology. Nonetheless, the processes controlling the distribution, transport and fate of cryoconite particles in the glacial system remain somewhat poorly understood. Here, using a 216 hour time series of plot-scale (0.04 m2) images, we quantify the small-scale dynamics of cryoconite on Longyearbreen, Svalbard. We show significant fluctuations in the apparent cryoconite area and dispersion of cryoconite over the plot, within the 9 day period of observations. However, the net movement of cryoconite across the ice surface averaged only 5.3 mm d−1. High-resolution measurements of cryoconite granule motion showed constant, random motion but weak correlation with meteorological forcing factors and no directional trends for individual particle movement. The high-resolution time-series data suggest that there is no significant net transport of dispersed cryoconite material across glacier surfaces. The areal coverage and motion of particles within and between cryoconite holes appears to be a product of differential melting leading to changes in plot-scale microtopography, local meltwater flow dynamics and weather-dependent events. These subtle processes of cryoconite redistribution may be significant for supraglacial albedo and have bearing on the surface energy balance at the glacier scale.

The cold-based termini of polythermal glaciers are usually assumed to adhere strongly to an immobile substrate and thereby supply significant resistance to the flow of warm-based ice up-glacier. This compressive environment is commonly thought to uplift basal sediment to the surface of the glacier by folding and thrust faulting. We present model and field evidence from the terminus of Storglaciären, Sweden, showing that the cold margin provides limited resistance to flow from up-glacier. Ice temperatures indicate that basal freezing occurs in this zone at 10−1 −10−2 m a−1, but model results indicate that basal motion at rates greater than 1 m a−1 must, nevertheless, persist there for surface and basal velocities to be consistent with measurements. Estimated longitudinal compressive stresses of 20–25 kPa within the terminus further indicate that basal resistance offered by the cold-based terminus is small. These results indicate that where polythermal glaciers are underlain by unlithified sediments, ice-flow trajectories and sediment transport pathways may be affected by subglacial topography and hydrology more than by the basal thermal regime.

Stakes at 2 km intervals were installed in January 1997 and remeasured in February 1998, January 1999, January 2005 and during the 2007/08 austral summer along a 1248 km traverse route from Zhongshan station to Dome A, East Antarctica. Based on topographical parameters, meteorological features and the records of ∼650 stakes and six stake arrays, the route is divided into five zones. We find that the snow accumulation rate decreases with increasing altitude as one progresses inland, except in the zone 800–1128 km from the coast, where the average annual accumulation rate is higher than in the zone 524–800 km from the coast. The Dome A zone (1128–1248 km) has the lowest accumulation rate (35 kg m−2 a−1, 2005–08) due to having the highest elevation and being furthest from the coast. The surface mass balance in the region 202–1128 km from the coast exhibits no temporal change from 1999–2005 to 2005–08, but there is a change in the accumulation distribution. The zone from 202 to 524 km shows a decrease in surface mass balance from 84 kg m−2 a−1 in 1999–2005 to 67 kg m−2 a−1 in 2005–08, while the zone between 800 and 1128 km shows an increase from 67 kg m−2 a−1 in 1999–2005 to 75 kg m−2 a−1 in 2005–08.

Surface ponds on Antarctic fast ice were examined by measuring temperature, salinity and concentrations of chlorophyll a (Chl-a), dissolved inorganic carbon (DIC) and nutrients (NO3 + NO2, PO4 and SiO2) in the surface pond water and under-ice water. Sea-ice cores were also collected from the bottom of a surface pond (pond-ice core) and from a site away from the pond (bare-ice core). Time-series measurements of surface pond water temperature showed that it varied with solar radiation rather than with air temperature. Comparison of water properties between surface pond water and under-ice water suggested that DIC and nutrients were consumed by biological productivity during pond formation. Depth profiles of nutrient concentrations in the pond-ice core suggested the remineralization of organic matter at the bottom of the surface pond. The Chl-a concentration was lower at the bottom of the pond-ice core than in the bare-ice core, suggesting that surface pond formation reduces ice algae abundance in sea ice because meltwater flushes algae from the porous sea ice into the under-ice water.

To investigate the subsurface hydrological characteristics of an overdeepened cirque glacier, nine boreholes were drilled to the bed of West Washmawapta Glacier, British Columbia, Canada, in summer 2007. All holes were surveyed with a video camera, and four were subsequently instrumented with a combination of pressure transducers, thermistors and conductivity sensors. Diurnal pressure and temperature records indicate the presence of a hydraulically connected subglacial drainage system towards the northern glacier margin. Hydraulic jacking in the overdeepening, controlled by changing water volume in the marginal zone, potentially impacts basal ice flow and erosion. The presence of a sediment layer underlying the glacier also likely impacts hydrology and ice dynamics. Influx of warm groundwater into the basal system raises subglacial water temperatures above the pressure-melting point (pmp) and induces diurnal water temperature fluctuations of as much as 0.8°C; water temperatures above the pmp could affect basal melt rates and the development of subglacial drainage systems. These observations suggest that the characteristics of the subglacial drainage system substantially affect patterns of flow and erosion by this small cirque glacier.

There are many advantages to formulating an ice-sheet model in terms of a variational principle. In particular, this applies to the specification of boundary conditions, which might otherwise be problematic to implement. Here we focus primarily on the frictional basal sliding boundary condition in a non-Newtonian Stokes model. This type of boundary condition is particularly difficult because it is heterogeneous, requiring both a Dirichlet (no-penetration) condition normal to the bed, and a Neumann (frictional sliding) condition tangential to the bed. In general, Neumann conditions correspond to natural boundary conditions in a variational principle; that is, they arise naturally in the variational formulation and thus need not be explicitly specified. While the same is not necessarily true of Dirichlet conditions, it is possible to enforce a no-penetration condition using Lagrange multipliers within the variational principle so that the Dirichlet condition becomes a natural boundary condition. Thus, in the case of ice sheets, all relevant boundary conditions may be incorporated in the variational functional, making them particularly easy to discretize. For the Stokes model, the resulting basal boundary condition is valid for arbitrary topographic slopes. Here we apply the same methodology to the Blatter– Pattyn higher-order approximate model, which is ordinarily limited to small basal slopes by the smallaspect-ratio approximation. We introduce a modification that improves on the accuracy of the standard Blatter–Pattyn model for all values of the basal slope, as we demonstrate in the slow sliding regime for which analytical results are available. The remaining error is due to the effects of the small-aspect-ratio approximation in the Blatter–Pattyn model.

In snow, acoustic emissions originate from the breaking of bonds between snow crystals and the formation of cracks. Previous research has shown that acoustic signals emanate from a natural snowpack. The relation between these signals and the stability of the snowpack has thus far remained elusive. Studies on other hazardous gravitational processes suggest that damage accumulation precedes major failure. If increased cracking activity could be detected in snow this might be used for avalanche prediction. We report on the development of a seismic sensor array to continuously monitor acoustic emissions in an avalanche start zone. During three winters, over 1400 sensor days of continuous acoustic data were collected. With the aid of automatic cameras and a microphone the main types of background noise were identified. Seismic signals generated by avalanches were also identified. Spectrograms from seismic signals generated by avalanches exhibit a unique triangular shape unlike any source of background noise, suggesting that automatic detection and classification of events is possible. Furthermore, discriminating between loose-snow and snow-slab avalanches is possible. Thus far we have not identified precursor events for natural dry-snow slab avalanche release. Detailed investigation of one dry-snow slab avalanche showed that signals observed prior to the release originated from background noise or small loose-snow avalanches.

Snow avalanches possess two types of kinetic energy: the kinetic energy associated with the mean velocity in the downhill direction and the kinetic energy associated with individual particle velocities that vary from the mean. The mean kinetic energy is directional; the kinetic energy associated with the velocity fluctuations is non-directional in the sense that it is connected to random particle movements. However, the rigid, basal boundary directs the random fluctuation energy into the avalanche. Thus, the random energy flux is converted to free mechanical energy which lifts and dilates the avalanche flow mass, changing the flow density and increasing the normal (dispersive) pressure and, as a consequence, changing the flow resistance. In this paper we derive macroscopic relations that link the production of the random kinetic energy to the perpendicular acceleration of the avalanche’s center of mass. We show that a single burst of fluctuation energy will produce pressures that oscillate around the hydrostatic pressure. Because we do not include a damping process, the oscillations of the center of mass remain, even if the production of random kinetic energy stops. We formulate relationships that can be used within the framework of depth-averaged mass and momentum equations that are often used to simulate snow avalanches in realistic terrain.

We have measured the glacier area changes in the central Southern Alps, New Zealand, between 1978 and 2002 and have compiled the 2002 glacier outlines using an image scene from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER). Three automated classification methods were tested: (1) band ratio, (2) normalized-difference snow index and (3) supervised classification. The results were compared with the glacier outlines photo-interpreted from the ASTER data, and were further validated using GPS-aided field mapping of selected test glaciers. The ASTER 3/4 band ratio provided the best results. However, all the classification methods failed to extract extensive debris-covered parts of the glaciers. Therefore, the photo-interpreted 2002 outlines were used when comparing with the existing 1978 glacier inventory derived from aerial photographs. Our results show a ∼17% reduction of glacier area, mainly driven by the retreat of the large valley glaciers. Despite the large climatic gradient from west to east, glaciers on both sides of the Main Divide lost similar percentages of area, except Franz Josef and Fox Glaciers which advanced. Smaller glaciers were found to have changed very little in the study period.

We derive mass changes of the Greenland ice sheet (GIS) for 2003–07 from ICESat laser altimetry and compare them with results for 1992–2002 from ERS radar and airborne laser altimetry. The GIS continued to grow inland and thin at the margins during 2003–07, but surface melting and accelerated flow significantly increased the marginal thinning compared with the 1990s. The net balance changed from a small loss of 7 ± 3 Gt a−1 in the 1990s to 171 ± 4 Gt a−1 for 2003–07, contributing 0.5 mm a−1 to recent global sea-level rise. We divide the derived mass changes into two components: (1) from changes in melting and ice dynamics and (2) from changes in precipitation and accumulation rate. We use our firn compaction model to calculate the elevation changes driven by changes in both temperature and accumulation rate and to calculate the appropriate density to convert the accumulation-driven changes to mass changes. Increased losses from melting and ice dynamics (17–206 Gt a−1) are over seven times larger than increased gains from precipitation (10–35 Gt a−1) during a warming period of ∼2 K (10 a)−1 over the GIS. Above 2000 m elevation, the rate of gain decreased from 44 to 28 Gt a−1, while below 2000 m the rate of loss increased from 51 to 198 Gt a−1. Enhanced thinning below the equilibrium line on outlet glaciers indicates that increased melting has a significant impact on outlet glaciers, as well as accelerating ice flow. Increased thinning at higher elevations appears to be induced by dynamic coupling to thinning at the margins on decadal timescales.

Rain-on-snow events trigger immediate and delayed avalanches as liquid water penetrates the snowpack. We present results from an extreme rain-on-snow event that triggered a glide avalanche near Snoqualmie Pass, Washington, USA. Snoqualmie Pass recorded 463 cm of snowfall from 13 December 2008 to 6 January 2009. This period of snowfall was followed by a strong southwesterly tropical flow that resulted in an extreme rain-on-snow event. Sensors at Snoqualmie Pass recorded 285 mm of precipitation over a 52 hour period. Flooding, slush flows, landslides and avalanches resulted from the influx of precipitation. Snow heights decreased rapidly over the period, with settlement rates approaching 80 mm h−1. Liquid water infiltrated and flowed through the snowpack within a few hours of the arrival of rain, yet many of the major avalanches occurred 12–30 or more hours after the onset of rain and water outflow. A glide avalanche occurred ∼30 hours after the onset of rain and the establishment of drainage through the snowpack. Increasing glide rates correlate with periods of rapid snow settlement. Here glide rates approached 670 mm h−1. Although glide and settlement rates increased during periods of intense precipitation, glide failure occurred 8 hours after peak precipitation and outflow.

Historical records are valuable for assessing glacier change in the Canadian High Arctic. Ommanney’s (1969) detailed inventory of Axel Heiberg Island glaciers, based on photography from 1958–59, has been revisited, converted into digital format and compared to glacier extents mapped from 1999–2000 satellite imagery. Our results show that the island-wide ice coverage decreased by 15.92 km2 in the 42 year period, a loss of <1%. However, two trends are apparent: one of advance or minor retreat from basins hosting outlet glaciers from Müller and Steacie Ice Caps, and one of significant retreat, on the order of 50–80%, for independent ice masses, which include valley glaciers, mountain glaciers, glacierets, and ice caps smaller than 25 km2. If the contributions to ice advance of only three surging glaciers are removed, then the island-wide ice loss approaches 60 km2. Furthermore, it is notable that 90% of ice masses smaller than 0.2 km2 disappeared entirely during the 42 year study period, an observation confirmed by field studies. Successful predictions from the original inventory are highlighted and the likely mechanisms driving the observed advances and retreats are discussed.

Impact pressures of snow avalanches have been measured at the Swiss Vallée de la Sionne experimental test site using two kinds of sensor placed at different locations in the avalanche flow. Pressures measured in a fast dry-snow avalanche and a slow wet-snow avalanche are compared and discussed. The pressures recorded using the two types of sensor in the dense flow of a dry-snow avalanche agree well, showing negligible dependence on the measurement device. On the other hand, significantly different pressures are measured in the slow dense flow of a wet-snow avalanche. This is attributed to the slow drag and bulk flow of this type of avalanche, leading to the formation and collapse of force-chain structures against the different surfaces of the sensors. At a macroscopic scale, limit state analysis can be used to explain such a mechanism by a shear failure occurring between freely flowing snow and a confined snow volume against the sensor, according to a Mohr–Coulomb failure criterion. The proposed model explains (1) how impact pressure can be up to eight times higher than hydrostatic snow pressure in wet cohesive slow avalanches and (2) its dependence on sensor geometry.

We have developed a technique in which we use a borehole video camera and post-processing software to make a record of the optical brightness as a function of depth in polar firn. We call this method borehole optical stratigraphy. To measure firn compaction, we note the positions of optical features on the borehole wall detected by an initial ‘baseline’ log. We track the displacements of these features in subsequent logs. The result provides a measurement of the relative vertical motion and thus compaction of the firn over the survey period. We have successfully used this system at Summit, Greenland, to measure the depth distribution of firn column shortening experienced in a borehole over three 1 year periods. The uppermost 30 m of the firn at Summit is compacting as predicted by a simple steady-state model, implying that the firn density profile at Summit is at or close to steady state over the past ∼70 years.

We report on satellite and ground-based observations that link glacier motion with subglacial hydrology beneath Skeiðarárjökull, an outlet glacier of Vatnajökull, Iceland. We have developed a technique that uses interferometric synthetic aperture radar (InSAR) data, from the European Remote-sensing Satellite (ERS-1/-2) tandem mission (1995–2000), to detect localized anomalies in vertical ice motion. Applying this technique we identify an area of the glacier where these anomalies are frequent: above the subglacial course of the river Skeiðará, where we observed uplift of 0.15–0.20 m d−1 during a rainstorm and a jökulhlaup, and subsidence at a slower rate subsequent to rainstorms. A similar pattern of motion is apparent from continuous GPS measurements obtained at this location in 2006/07. We argue that transient uplift of the ice surface is caused by water accumulating at the glacier base upstream of an adverse bed slope where the overburden pressure decreases significantly over a short distance. Most of the frictional energy of the flowing water is therefore needed to maintain water temperature at the pressure-melting point. Hence, little energy is available to enlarge water channels sufficiently by melting to accommodate sudden influxes of water to the base. This causes water pressure to exceed the overburden pressure, enabling uplift to occur.

Taylor Glacier, an outlet of the East Antarctic ice sheet, flows through the Transantarctic Mountains and terminates in the Dry Valleys. Understanding how this glacier fluctuates is important for studies of glacial geology, paleoclimate, ice dynamics and ecology. Sublimation is the primary mass-loss process for most of the glacier. Four years of specific balance measurements from the ablation zone show sublimation rates up to 40 cm a−1. We used data from an array of weather stations as inputs to a model for latent heat flux and hence sublimation rate. Calculated and measured ablation rates agree to within uncertainties, indicating that wind speed and vapour pressure gradient (a function of temperature and humidity) are the governing variables, as expected from theory. Measurements and model results together allowed us to examine the spatial and temporal variations of sublimation on the glacier. On average, sublimation is about two times faster in summer than winter. Rapid sublimation occurs during storms and katabatic wind events, but such periods contribute less to the annual total than do slow, persistent losses. Spatially, sublimation reaches a maximum midway along the glacier, where descending surface air currents are focused by the topography of the aptly named tributary, Windy Gully.

Ice crystal orientation fabric data from ice cores contain important information concerning the internal structure and the flow behaviour of ice sheets. When ice cores are recovered from the Antarctic and Greenland ice sheets, crystal orientation measurements are performed immediately to obtain fundamental physical property information. In the past, we have measured the c-axis orientation of ice crystals by a manual optical method using a universal stage. This method is very time-consuming, involving tedious work in a cold laboratory. Recently, automated systems have been developed that enable measurement of c-axis orientation, grain size and other microstructures. However, in order to detect the full crystal orientation of an ice crystal, we also need information on its a-axis orientation. A variety of other crystal orientation measurement methods have previously been discussed, but some shortcomings for ice-core studies cannot be neglected. We have developed a crystal-orientation analysing device using the Laue X-ray diffraction method. As this device can measure the orientations of all crystal axes with high accuracy, it is possible to obtain new microstructure information on natural ice crystals. For the first time, we are able to quantify very low subgrain misorientation angles in polar ice-core samples, allowing us to investigate micro-deformation features of individual crystals. Here we discuss the analysis process, which is customized to measure standard ice thin sections, and show preliminary results.

Spectral measurements of solar radiation in the band 400–900 nm were performed above and inside the snowpack in two locations in Finland, using a spectroradiometer. The transmittance and extinction coefficient were estimated for different snow layers. Four small candle-shaped photosynthetically active radiation (PAR) sensors were also used to measure the transfer of PAR inside the snowpack. In addition to the light measurements, physical characterization of snow stratigraphy was done, including the thickness, density, hardness (hand test), salinity, and grain size and shape (photographs of crystals). The transmittance varied from <1% (0–12 cm layer) to 80% (0–4 cm layer), and the extinction coefficient was between 0.03 cm−1 (4–8 cm layer) and 0.8 cm −1 (0–4 cm layer). The physical properties of the snow varied considerably between locations and days. The density of the snow varied between 140 and 480 kg m−3.

Ice cores from small ice caps provide valuable climatic information, additional to that of Greenland and Antarctica. However, their integrity is usually compromised by summer meltwater percolation. To determine to what extent this can affect such ice cores, we performed high-resolution tritium measurements on samples from two ice cores from Spitsbergen covering the period AD1955–75. The very sharp and distinct peaks in the tritium precipitation record are subject to several post-depositional processes. We developed a model that uses the precipitation record as input and incorporates the three most important processes (radioactive decay, isotope diffusion and meltwater percolation). Results are compared with measured tritium and density profiles. Both ice-core records contain sharp bomb peaks in the pre-1963 period. It is shown that these peaks would be much smoother in the absence of melt. In this case the main effect of melt and the refreezing of percolation water is the formation of ice layers that form barriers for firn diffusion; thus melt paradoxically results in better preservation of the annual isotope signals. Conversely, for the period after 1963 the main effect of melt is a stronger smoothing of the tritium profiles.