Through the detection of surface deformation in response to water movement, recent satellite studies have demonstrated the existence of subglacial lakes in Antarctica that fill and drain on timescales of months to years. These studies, however, were confined to specific regions of the ice sheet. Here we present the first comprehensive study of these ‘active’ lakes for the Antarctic ice sheet north of 86° S, based on 4.5 years (2003–08) of NASA’s Ice, Cloud and land Elevation Satellite (ICESat) laser altimeter data. Our analysis has detected 124 lakes that were active during this period, and we estimate volume changes for each lake. The ICESat-detected lakes are prevalent in coastal Antarctica, and are present under most of the largest ice-stream catchments. Lakes sometimes appear to transfer water from one to another, but also often exchange water with distributed sources undetectable by ICESat, suggesting that the lakes may provide water to or withdraw water from the hydrologic systems that lubricate glacier flow. Thus, these reservoirs may contribute pulses of water to produce rapid temporal changes in glacier speeds, but also may withdraw water at other times to slow flow.

Many glacier snouts in the Himalaya are known to be stagnant and exhibiting low surface gradients, conditions that are conducive to the formation of glacial lakes impounded either by the terminal moraine or by the remnant glacier snout. In this study, we use interferometry and feature-tracking techniques to quantify the extent of stagnation in 20 glaciers across the Everest (Qomolangma; Sagarmatha) region, and subsequently we examine the relationship between local catchment topography and ice dynamics. The results show that only one of the studied glaciers, Kangshung Glacier, is dynamic across its entire surface, with flow rates greater than 40 m a−1 being recorded in high-elevation areas. Twelve other glaciers show some evidence of flow, but are generally characterized by long, stagnant tongues, indicating widespread recession and in situ decay. The remaining seven glaciers show no evidence of flow in any of the available datasets. Hypsometric data suggest that catchment topography plays an important role in controlling glacier flow regimes, with those fed by wide, high-altitude accumulation areas showing the most extensive active ice, and those originating at low elevations exhibiting large areas of stagnant ice. Surface profiles extracted from a SRTM digital elevation model indicate that stagnant snouts are characterized by very low (<2°) surface angles and that down-wasting is the prevalent ablation pattern in the study area.

Glaciers are widely recognized as key indicators of climate change, and their meltwater plays an important role in hydropower production in Norway. Since the last glacier inventory was compiled in northern Norway in the 1970s, marked fluctuations in glacier length and mass balance have been reported for individual glaciers, and the current overall glacier state is thus not well known. Within the framework of the Global Land Ice Measurements from Space (GLIMS) initiative, we have created a new inventory for 489 glaciers in the Svartisen region, northern Norway, using a Landsat Enhanced Thematic Mapper Plus (ETM+) satellite scene from 7 September 1999 and automated multispectral glacier mapping (thresholded band ratios). In addition, visual inspection and correction of the generated glacier outlines has been applied. Adverse snow conditions and uncertain drainage divides made glacier mapping challenging in some regions of the study site. Glacier outlines from 1968, as digitized from a topographic map, were used for a quantitative change assessment for a selection of 300 glaciers. The overall area change of this sample from 1968 to 1999 was close to zero, but with a strongly increasing scatter towards smaller glaciers, large area gains (>50%) for small glaciers (<1 km2), and an unexpected stronger relative area loss towards the wetter coast. The overall size changes are small (<1%) for the three largest ice masses in the study region (Vestisen, Østisen and Blåmannsisen).

In 2003–07 summer balance was measured at altitudes between 1700 and 3382 m a.s.l. on two glaciers on Mount Rainier, Washington State, USA (46.85° N,121.72° W; 4400 m a.s.l.): south-facing Nisqually Glacier and east-northeast-facing Emmons Glacier. Upper-air temperatures at the nearest gridpoint in the NCEP/NCAR reanalysis database are used in a distributed (over altitude) positive-degreeday (PDD) model. For each glacier the model used the same coefficients at all altitudes, for all years. The rms model error was 0.65 (r2 = 0.87) and 0.78 m a−1 w.e. (r2 = 0.93) for Nisqually and Emmons Glaciers, respectively. Although PDD work generally uses different coefficients for snow and ice surfaces, and the duration of exposure of those surfaces varies with altitude, error in this single-coefficient model is nearly uncorrelated with altitude. Values of coefficients obtained are within the range of those found in other PDD work. The degree-day coefficient, however, differs markedly between the two glaciers, and is shown to be controlled by the difference between them in vertical gradient of measured summer balance. It is smaller for Nisqually Glacier, where solar radiation is a stronger contributor to melt; and larger for Emmons Glacier, where it is a weaker contributor. Over 1948–2007, when the model calibrated over 2003–07 was applied to the upper-air temperatures, estimated summer balance was ∼0.4 m a−1 less negative over 1962–83 than before and ∼0.6 m a−1 less negative than after, corresponding roughly with changes of the northeast Pacific sea-surface temperatures.

The first three-dimensional properties of polar firn obtained by X-ray microtomography are used to study the microstructure of snow on a 15 m deep firn core from West Antarctica. The snow is found to undergo coarsening down to approximately 2.5 m depth before grain growth and densification become the prevalent mechanisms of microstructure change. In contrast to previous assumptions, distinct anisotropy of the ice and pore geometry is observed throughout the profile, with a maximum at 2.5 m depth. The air permeability and the degree of anisotropy vary with depth and can be linked to short-term changes in accumulation rate via the residence time for which a certain snow layer stays in the uppermost 2.5 m. Patterns of the degree of anisotropy and air permeability of buried polar firn are relative indicators of past accumulation rates.

Visualization and quantification of snow structures at a scale of a few millimetres is important in understanding the mechanical, thermal and electromagnetic properties of snow. Surface sections and, to an even greater degree, three-dimensional (3-D) reconstructions of cast snow samples are difficult to prepare, and automatic image processing is notoriously difficult and often requires manual evaluation. Here, we present a new method to measure the 3-D structure of cast snow samples. Snow samples cast with diethyl phthalate (DEP) and frozen are cut to a sample size a few centimetres in diameter and up to 10 cm in height. The ice of these samples is then sublimated in high vacuum and the remaining negative structure (replica) is imaged using X-ray microtomography (micro-CT). The accuracy of the method is demonstrated by comparing micro-CT scans of the original snow structure and the replica. The method described here allows easy transportation of samples, requires little manual interaction, has a very high spatial resolution of up to 10 μm and is environmentally friendly.

We use a simple energy-conservation model and a model based on Röthlisberger’s theory for steady-state water flow in a subglacial conduit to model water movement between lakes in the Adventure subglacial trench region of East Antarctica during a 1996–98 jökulhlaup. Using available field evidence to constrain the models suggests that water flow would likely be accommodated in a tunnel with a cross-sectional area of 36 m2 and a value for k (the reciprocal of Manning’s roughness parameter) larger than the 12.5 m1/3 s−1 previously calculated. We also use Nye’s theory for time-dependent conduit water flow to model the temporal evolution of conduit discharge, cross-sectional area, water pressure and lake draining and filling during the flood. We initially assume one source and one sink lake. We perform sensitivity tests on the input parameter set, matching modeled source- and sink-lake depth changes with measured surface elevation data. Using a simple function for vertical ice deformation in which surface deformation scales linearly to the lake depth change, we find the scaling factor is of the order 4 × 10−3 of the ice thickness. The most likely value of k lies in the range 55–68 m1/3 s−1, and the ratio of source to sink-lake radii is approximately 1 : 1.4. Finally, we experiment using Nye’s theory to model water movement between one source and three sink lakes. The model fails to produce the observed patterns of water movement as indicated by the surface deformation data.

Recurring jökulhlaups from ice-dammed lakes often form irregular time sequences that are seemingly unpredictable. Using the flood dates of Merzbacher Lake, Kyrgyzstan, as an example, we study these sequences through a model of lake filling and drainage where flood events initiate at a threshold water depth. Even with a constant threshold, model simulation can explain key aspects of the Merzbacher flood sequence. General analysis of model dynamics reveals a pacing mechanism that links one flood to the next, and which may be represented mathematically as an iterative map. This theory clarifies how environmental factors govern the long-term pattern of flood timings and their frequency distribution in the year. A reconstruction of the past level of Merzbacher Lake also suggests that its flood-initiation threshold decreases with the rate of lake-level rise. These results may help us understand how to forecast future outbursts from jökulhlaup lakes.

Estimates of glacier mass balance using geodetic methods can differ significantly from estimates using direct glaciological field-based measurements. To determine if such differences are real or methodological, there is a need to improve uncertainty estimates in both methods. In this paper, we focus on the uncertainty of geodetic methods and describe a geostatistical technique that takes into account the spatial correlation of the elevation differences when calculating spatially averaged elevation changes. We apply this method to the western Svartisen ice cap, Norway, using elevation differences from the surrounding bedrock derived from stereophotogrammetry. We show that the uncertainty is not only dependent on the standard error of the individual elevation differences but is also dependent on the size of the averaging area and the scale of the spatial correlation. To assess if the geostatistical analysis made over bedrock is applicable to glacier surfaces, we use concurrent photogrammetrical and laser scanning data from bedrock and a range of glacier surfaces to evaluate the dependency of the geostatistical analysis on the surface type. The estimated geodetic mass balance, and its uncertainty, is −2.6 ± 0.9 m w.e. for the period 1968–85, and −2.0 ± 2.2 m w.e. for 1985–2002.

Fourteen englacial conduits were mapped within 2 km of the terminus of the temperate Matanuska Glacier, Alaska, USA, to ice depths of 65 m using speleological techniques. Detailed three-dimensional maps of the conduits were made over 3 years to characterize conduit relationships with glacier structural features and to track conduit evolution through time. All conduits consisted of single unbranching passages that followed fractures in the ice. All conduits were either too constricted to continue or became water-filled at their deepest explored point and were not able to be followed to the glacier bed. Conduit morphology varied systematically with the orientation of the glacier principal stresses, allowing them to be categorized into two broad classes. The first class of conduits were formed by hydrostatic crevasse penetration where a large supraglacial stream intersected longitudinal crevasses. These conduits plunged toward the glacier bed at angles of 30–40°. The second class of conduits formed where smaller streams sank into the glacier on shear crevasses. Many of these conduits changed direction dramatically where they intersected transverse crevasses at depth. These results suggest that the conduits observed in this study formed along fractures and, over their surveyed length, were not affected by gradients in ice overburden pressure.

An analysis of the local topographic setting of very small (<0.4 km2) glaciers within a small region of the Monashee Mountains, British Columbia, was conducted to investigate its influence on recent changes in the extent of these glaciers. Net changes in glacier extent were determined from a detailed manual comparison of remotely sensed imagery acquired in 1951, 2001 and 2004. Most of the 86 glaciers included in the study showed no observable net change in area over the study period, while six glaciers retreated, four disappeared entirely and only one advanced. Indices derived to characterize elements of the local topographic setting that might affect the local mass balance suggest that most of the glaciers are situated in locations that favor ice preservation by enhancing mass input and/or reducing ablation rates. Glaciers situated in less favorable settings generally either decreased in area or disappeared. The results suggest that most of the glaciers studied have retreated as far as they are likely to under the climatic conditions of the late 20th century.

Spectral reflectances in the visible to near-infrared wavelength range (350–1050 nm) were measured on the surface of Gulkana Glacier, Alaska Range, Alaska, USA, at six different elevations from May to September 2001. The measurements showed variable spectral reflectances on the glacier surface, and the spectra changed over time and elevation. The spectra in May were those of a typical wet snow surface: the reflectance ranged approximately from 0.4 to 0.8, and generally decreased as the wavelength increased. As snowmelt proceeded, the reflectance gradually fell, and specific absorptions appeared in the spectrum at wavelength ranges of 400–600 and 670–680 nm, corresponding to a red snow surface due to snow algal bloom. As the ice surface was exposed, the spectrum finally shifted to a lower (0.1–0.4) and flat reflectance curve, typical of an ice surface contaminated with dust. Analyses of the dust on the glacier surface showed that its quantity and characteristics changed spatially and seasonally and were correlated with spectral albedo. The results showed that the temporal and spatial variations in spectral reflectance on this glacier are due not only to physical properties of the glacier surface, but also to biogenic materials on the surface (e.g red-colored snow algal cells and cryoconite).

Since the mid-19th century, most glaciers have been losing area and volume. This loss of area has not been homogeneous in time and space, and direct observations are sparse, making regional and global estimates of glacier change difficult. This paper focuses on developing a single index for monitoring glacier change, one that would be particularly useful for remote-sensing applications. We combine the results of direct glacier mass-balance observations B, total glacier area S and accumulation area Sc derived from maps or remotely sensed images. Using the accumulation–area ratio (AAR = Sc/S), we note the differences between observed AAR, time-averaged 〈AAR〉 and the equilibrium state AAR0, as determined by its value at B = 0 from a regression of B(AAR). We suggest that αd = (〈AAR〉 – AAR0)/AAR0 quantifies the difference between the currently observed state of glaciers and their equilibrium state and measures the delay in the dynamic response of S relative to the climatic response of Sc. Using all available observations for the period 1961–2004, αd ≈ −65% for tropical glaciers, which implies their rapid shrinkage as S continues to decrease and ‘catch up’ with Sc. During the same period, mid-latitude and polar glaciers show less negative values of αd. Of 86 glaciers from all latitudes and regions, only 11 show positive αd at any time between 1961 and 2004. Averaged over 1961–2004, αd is −15.1 ± 2.2%, and 〈B〉 is −360 ± 42 mm a−1 w.e. Values for AAR0 range between about 40% and 80%, but the bulk of the equilibrium values are between 50% and 60%. The average AAR0 is 57.9 ± 0.9% and has remained stable over time (the equilibrium AAR has not changed with climate). Overall, the observed negative αd suggests a committed retreat of glaciers and their continuing contribution to sea level even if global temperature is held constant.

The Amery Ice Shelf, East Antarctica, undergoes high basal melt rates near the southern limit of its grounding line where 80% of the ice melts within 240 km of becoming afloat. A considerable portion of this later refreezes downstream as marine ice. This produces a marine ice layer up to 200 m thick in the northwest sector of the ice shelf concentrated in a pair of longitudinal bands that extend some 200 km all the way to the calving front. We drilled through the eastern marine ice band at two locations 70 km apart on the same flowline. We determine an average accretion rate of marine ice of 1.1 ± 0.2 m a−1, at a reference density of 920 kg m−3 between borehole sites, and infer a similar average rate of 1.3 ± 0.2 m a−1 upstream. The deeper marine ice was permeable enough that a hydraulic connection was made whilst the drill was still 70–100 m above the ice-shelf base. Below this marine close-off depth, borehole video imagery showed permeable ice with water-filled cavities and individual ice platelets fused together, while the upper marine ice was impermeable with small brine-cell inclusions. We infer that the uppermost portion of the permeable ice becomes impermeable with the passage of time and as more marine ice is accreted on the base of the shelf. We estimate an average closure rate of 0.3 m a−1 between the borehole sites; upstream the average closure rate is faster at 0.9 m a−1. We estimate an average porosity of the total marine ice layer of 14–20%, such that the deeper ice must have even higher values. High permeability implies that sea water can move relatively freely through the material, and we propose that where such marine ice exists this renders deep parts of the ice shelf particularly vulnerable to changes in ocean properties.

The automatic weather station (AWS) on the snout of the Vadret da Morteratsch, Switzerland, has delivered a unique 12 year meteorological dataset from the ablation zone of a temperate glacier. This dataset can be used to study multi-annual trends in the character of the surface energy budget. Since 2003 there has been a substantial darkening of the glacier tongue due to the accumulation of mineral and biogenic dust. The typical surface albedo in summer has dropped from 0.32 to 0.15. We have analysed the implications of the lowered albedo for the energy balance and the annual ablation. For the 4 year period 2003–06, the decreased albedo caused an additional removal of about 3.5 m of ice. Calculations with an energy-balance model show that the same increase in ablation is obtained by keeping the ice albedo fixed to 0.32 and increasing the air temperature by 1.7 K. Our analysis confirms that for retreating glaciers the deposition of dust from exposed side moraines on the glacier surface constitutes an important feedback mechanism. The mineral dust stimulates the growth of algae, lowers the surface albedo, enhances the melt rates, and thereby facilitates the further retreat of the glacier snout.

Snow cover in the Tibetan Plateau is highly variable in space and time and plays a key role in ecological processes of this cold-desert ecosystem. Resolution of passive microwave data is too low for regional-scale estimates of snow cover on the Tibetan Plateau, requiring an alternate data source. Optically derived snow indices allow for more accurate quantification of snow cover using higher-resolution datasets subject to the constraint of cloud cover. This paper introduces a new optical snow index and assesses four optically derived MODIS snow indices using Landsat-based validation scenes: MODIS Snow-Covered Area and Grain Size (MODSCAG), Relative Multiple Endmember Spectral Mixture Analysis (RMESMA), Relative Spectral Mixture Analysis (RSMA) and the normalized-difference snow index (NDSI). Pearson correlation coefficients were positively correlated with the validation datasets for all four optical snow indices, suggesting each provides a good measure of total snow extent. At the 95% confidence level, linear least-squares regression showed that MODSCAG and RMESMA had accuracy comparable to validation scenes. Fusion of optical snow indices with passive microwave products, which provide snow depth and snow water equivalent, has the potential to contribute to hydrologic and energy-balance modeling in the Tibetan Plateau.

Interest in crevasses and associated ice-fracture processes has recently increased due to recognition of the importance of calving glaciers to the mass balance of the cryosphere, as well as the importance of fractures in glacier hydrology. Recently developed calving criteria make use of models which predict crevasse depth from surface strain rates, but these models have rarely been tested against observations. In this study, we present data on crevasse depth and surface strain rates, and compare the measured values with results of two crevasse-depth models: a simple function proposed by Nye and a linear elastic fracture mechanics (LEFM) model developed by Van der Veen. Our results indicate that both models predict crevasse depths within the correct order of magnitude. The LEFM model, incorporating measured values of crevasse spacing and tuned for fracture toughness, performs better in predicting crevasse depths, but where lack of input data precludes such tuning, the results are similar to Nye’s model predictions. We conclude that both models may be used to calculate crevasse depths in calving models, although the Nye function is undoubtedly much simpler to implement within an ice-dynamics model.