We have computed the velocity of Langjökull ice cap (LIC), Iceland, using InSAR data acquired by the European Remote-sensing Satellite in 1994 and an ice-surface elevation model derived from a GPS survey conducted in 1997. The velocity map covers 68% of the 925 km2 of LIC. Eight principal outlet glaciers were identified, along which ice velocities increased from zero at the ice divide to a maximum of 75 m a−1. When the InSAR velocity dataset is compared with surface GPS measurements of ice flow determined in summer 2001, there is generally good agreement, suggesting interannual and seasonal stability for much of the ice cap during this period. However, non-zero ice velocities at the termini of four of the eight glaciers surveyed suggest that these glaciers were advancing in 1994. Comparison of our GPS and InSAR velocities suggests that Sudurjökull experienced elevated velocity in late February 1994, which, when other studies are considered, appears to be consistent with a glacier surge.

The recent dynamic behaviour of Karakoram glaciers is expected to differ from that shown by glaciers in the central and eastern Himalaya because of regional variations in precipitation and temperature trends. However, there are insufficient quantitative data to support or confute such hypotheses. We present velocity data covering the period 1993–2008 for Baltoro Glacier, one of the longest glaciers in the Karakoram. Velocity measurements were made using cross-correlation feature tracking applied to European Remote-sensing Satellite (ERS-1 and -2) and Envisat advanced synthetic aperture radar (ASAR) data, supplemented by differential global positioning system (DGPS) measurements. We find a gradual acceleration of the glacier during the early 2000s, in particular during winter months. Multi-seasonal data reveal a large difference between summer and winter flow characteristics, but only in the upper ablation zone. Summer 2005 was a particularly dynamic period following from the heavy winter snowfall of 2004, indicating the importance of basal meltwater availability for glacier flow. Transverse velocity profiles indicate that Baltoro Glacier undergoes ‘block’ flow across much of the upper ablation zone during the summer, which we interpret as evidence of widespread basal sliding. The DGPS data confirm the rapid increase in flow detected during 2005. Modelled climatic data reveal decreasing summer temperatures and increasing precipitation over the study period, helping to explain the observed dynamic variations and their differences from glaciers elsewhere in the Himalaya.

Inclusions affect the behavior of ice, and their characteristics help us understand the formation history of the ice. Recently, a low-temperature magnesium sulfate salt was discovered. This paper describes this naturally occurring MgSO4·11H2O mineral, meridianiite, derived from salt inclusions in sea ice of Lake Saroma, Japan and in Antarctic continental core ice. Its occurrence is confirmed by using micro-Raman spectroscopy to compare Raman spectra of synthetic MgSO4·11H2O with those of the inclusions.

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.

Using airborne and spaceborne high-resolution digital elevation models and laser altimetry, we present estimates of interannual and multi-decadal surface elevation changes on the Bering Glacier system, Alaska, USA, and Yukon, Canada, from 1972 to 2006. We find: (1) the rate of lowering during 1972–95 was 0.9 ± 0.1 m a−1; (2) this rate accelerated to 3.0 ± 0.7 m a−1 during 1995–2000; and (3) during 2000–03 the lowering rate was 1.5 ± 0.4 m a−1. From 1972 to 2003, 70% of the area of the system experienced a volume loss of 191 ± 17 km3, which was an area-average surface elevation lowering of 1.7 ± 0.2 m a−1. From November 2004 to November 2006, surface elevations across Bering Glacier, from McIntosh Peak on the south to Waxell Ridge on the north, rose as much as 53 m. Up-glacier on Bagley Ice Valley about 10 km east of Juniper Island nunatak, surface elevations lowered as much as 28 m from October 2003 to October 2006. NASA Terra/MODIS observations from May to September 2006 indicated muddy outburst floods from the Bering terminus into Vitus Lake. This suggests basal–englacial hydrologic storage changes were a contributing factor in the surface elevation changes in the fall of 2006.

The effects of temperature and seasonal air-mass trajectories on stable water isotopes in alpine snowpacks are investigated using meteorological and snow-pit data at two alpine field sites in the Canadian Rocky Mountains: Haig Glacier, Alberta, and Opabin Glacier, British Columbia. Snow pits were sampled through three accumulation seasons (October–June, 2004/05, 2005/06 and 2006/07) for δ18O, δD, temperature and density. The isotopic characteristics of precipitation over these time periods, including the local meteoric waterline and average δ18O, δD and deuterium excess, were defined using this dataset. Individual snowfall events over the three seasons were identified in the accumulation records from both sites and then fit to snow-pit stratigraphies to determine their mean isotopic characteristics. A trajectory classification was produced for all events, and the key meteorological characteristics of each trajectory class were investigated using data from alpine field sites and a suite of meteorological records from the region. An analysis of the relative influences of temperature and air-mass trajectory on snow isotope ratios reveals some separation in mean δ18O between storm classes. However, the separation appears to be driven primarily by the mean temperature of each class rather then being a direct effect of vapour pathway.

We present estimates of internal shear rates of real-scale avalanches that are based on velocity measurements. Optical velocity sensors installed on the instrument pylon at the Swiss Vallée de la Sionne test site are used to measure flow velocities at different flow heights of three large dry and wet snow avalanches. Possible sources of error in the correlation analysis of the time-lagged reflectivity signals measured by optical sensors are identified for real-size avalanches. These include spurious velocities due to noise and elongated peaks. An appropriate choice of the correlation length is essential for obtaining good velocity estimates. Placing restrictions on the maximum possible accelerations in the flow improves the analysis of the measured data. Coherent signals are found only in the dense flowing cores. We observe the evolution of shear rates at different depths between the front and tail of the flowing avalanche. At the front, large shear rates are found throughout the depth; at the tail, plug flows overriding highly sheared layers near the bottom of the flow are observed. The measured velocities change strongly with height above the ground and fluctuations around the measured mean velocity can be identified. We find that the dense flows are laminar, undergoing a transition from supercritical to subcritical flow behaviour from the head to the tail. Furthermore, we provide real-scale experimental evidence that the mean shear rate and the magnitude of velocity fluctuations increase with the mean discharge.

Melt factors for snow (ks) and ice (ki) were derived from specific mass-balance data and regionally interpolated daily air-temperature series at nine glaciers in the western Cordillera of Canada. Fitted ks and ki were relatively consistent across the region, with mean values (standard deviations) of 3.04 (0.38) and 4.59 (0.59) mm d−1 °C−1, respectively. The interannual variability of melt factors was investigated for two long-term datasets. Calculated annually, snow- and ice-melt factors were relatively stable from year to year; standard deviations for snowmelt factors were 0.48 (17%) and 0.42 (18%) at Peyto and Place Glaciers, respectively, while standard deviations of ice-melt factors were 1.17 (25%) and 0.81 (14%). While fitted values of ks are comparable to those presented in previous observational and modeling studies, fitted ki are substantially and consistently lower across the region. Fitted melt factors were sensitive to the choice of lapse rate used in the air-temperature interpolation. Melt factors fitted to mass-balance data from a single site (Place Glacier) provided reasonable summer balance predictions at most other sites representing both maritime and continental climates, although there was a tendency for under-prediction at several sites. The combination of regionally interpolated air temperatures and a degree-day model appears capable of generating first-order estimates of regional summer balance, which can provide a benchmark against which to judge the predictive ability of more complex (e.g. energy balance) models applied at a regional scale. Mass-balance sensitivity analyses indicate that a temperature increase of 1 K will increase summer ablation in the region by 0.51 m w.e. a−1 on average.

Discriminating between glacier variations due to natural climate variability and those due to true climate change is crucial for the interpretation and attribution of past glacier changes, and for the expectations of future changes. We explore this issue for the well-documented glaciers of Mount Baker in the Cascades Mountains of Washington State, USA, using glacier histories, glacier modeling, weather data and numerical weather model output. We find that natural variability alone is capable of producing kilometer-scale excursions in glacier length on multi-decadal and centennial timescales. Such changes are similar in magnitude to those attributed to a global Little Ice Age. The null hypothesis, that no climate change is required to explain the glacier fluctuations in this setting, cannot be rejected. These results for Mount Baker glaciers are also consistent with an earlier study analyzing individual glaciers in Scandinavia and the Alps. The principle that long-timescale fluctuations of glacier length can be driven by short-timescale fluctuations in climate reflects a robust and fundamental property of stochastically forced physical systems with memory. It is very likely that this principle also applies to other Alpine glaciers and that it therefore complicates interpretations of the relationship between glacier and climate history. However, the amplitude and timescale of the length fluctuations depends on the details of the particular glacier geometry and climatic setting, and this remains largely unevaluated for most glaciers.

We describe a high-precision method, now in use in our laboratory, for measuring the CO2 mixing ratio of ancient air trapped in polar ice cores. Occluded air in ice samples weighing ∼8–15 g is liberated by crushing with steel pins at −35°C and trapped at −263°C in a cryogenic cold trap. CO2 in the extracted air is analyzed using gas chromatography. Replicate measurements for several samples of high-quality ice from the Siple Dome and Taylor Dome Antarctic ice cores have pooled standard deviations of <0.9 ppm. This high-precision technique is directly applicable to high-temporal-resolution studies for detection of small CO2 variations, for example CO2 variations of a few parts per million on millennial to decadal scales.

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).

Increased mass losses from the Greenland ice sheet and inferred contributions to sea-level rise have heightened the need for hydrologic observations of meltwater exiting the ice sheet. We explore whether temporal variations in ice-sheet surface hydrology can be linked to the development of a downstream sediment plume in Kangerlussuaq Fjord by comparing: (1) plume area and suspended sediment concentration from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery and field data; (2) ice-sheet melt extent from Special Sensor Microwave/Imager (SSM/I) passive microwave data; and (3) supraglacial lake drainage events from MODIS. Results confirm that the origin of the sediment plume is meltwater release from the ice sheet. Interannual variations in plume area reflect interannual variations in surface melting. Plumes appear almost immediately with seasonal surface-melt onset, provided the estuary is free of landfast sea ice. A seasonal hysteresis between melt extent and plume area suggests late-season exhaustion in sediment supply. Analysis of plume sensitivity to supraglacial events is less conclusive, with 69% of melt pulses and 38% of lake drainage events triggering an increase in plume area. We conclude that remote sensing of sediment plume behavior offers a novel tool for detecting the presence, timing and interannual variability of meltwater release from the ice sheet.

Low-elevation blue-ice areas in Antarctica can contain significant amounts of subsurface liquid water during summer and may experience internal melt–freeze cycles due to absorption of shortwave radiation. An existing 2-D (x-z) model has been used to study the phenomenon and this is improved by changing the lower boundary condition from a no-flux condition to one that lets heat through, and by changing the bulk optical properties to spectral values. Both changes made the model more realistic. The optical and thermal boundary conditions, together with the optical attenuation coefficient, have a large effect on the amount of water produced in the ice. Our results show that if the lower boundary condition is changed from no flux to radiating, subsurface melting is reduced dramatically. Using a spectral albedo produces less melting than using a corresponding bulk albedo, the other variables left unchanged. If optical properties are changed to spectral values, the melting is different than using bulk values. If the linear extent of the blue-ice area is <0.6 m, subsurface melting is diminished. We found the spatial variability of snow has a significant effect on subsurface melting. Subsurface melting was found to be impossible for albedos >0.7 and the subsurface ponds can persist over the winter if the albedo is <0.4.

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.

Deep, pervasive shear deformation of the bed to high strains (>100) may have been primarily responsible for flow and sediment transport of the Lake Michigan lobe of the Laurentide ice sheet. To test this hypothesis, we sampled at 0.2 m increments a basal till from one advance of the lobe (Batestown till) along vertical profiles and measured fabrics due to both anisotropy of magnetic susceptibility and sand-grain preferred orientation. Unlike past fabric studies, interpretations were guided by results of laboratory experiments in which this till was deformed in simple shear to high strains. Fabric strengths indicate that more than half of the till sampled has a <5% probability of having been sheared to moderate strains (7–30). Secular changes in fabric azimuth over the thickness of the till, probably due to changing ice-flow direction as the lobe receded, indicate that the bed accreted with time and that the depth of deformation of the bed did not exceed a few decimeters. Orientations of principal magnetic susceptibilities show that the state of strain was commonly complex, deviating from bed-parallel simple shear. Deformation is inferred to have been focused in shallow, temporally variable patches during till deposition from ice.

We describe the development of a low-frequency airborne radar specifically designed for the sounding of temperate ice. The system operates at a central frequency of 1 MHz and consists of an impulse transmitter with an output voltage up to 5000 V and a digital receiver with a maximum gain of 80 dB. The radar was deployed on board a CASA 212 aircraft, which also carries a laser altimeter, an inertial navigation system, a digital camera and a GPS receiver. A description of the radar system is provided, as well as preliminary results obtained at Glaciar Tyndall, Campo de Hielo Sur (Southern Patagonia Icefield), where an ice depth of 670 m was reached.

Surface ice-flow velocities measured at stakes on Rutford Ice Stream, West Antarctica, covering a period of ∼25 years are analysed for evidence of temporal variations in flow. No indications of significant long-term changes in flow are found. Earlier observations have shown significant tidally related variations in flow speed. We conclude that temporal variability on Rutford Ice Stream, West Antarctica, appears limited to tidal periods of days and weeks, and weaker interannual variation, possibly related to long-period tides or seasonal effects, while long-term (decadal) changes in flow speed are either absent or smaller than ∼0.1% a−1.

This study fills the gap in hydrologic measurements of runoff exiting a part of the Greenland ice sheet (GrIS), the Kangerlussuaq drainage area, West Greenland. The observations are of value for obtaining knowledge about the terrestrial freshwater and sediment output from part of the GrIS and the strip of land between the ice sheet and the ocean, in the context of varying ice-sheet surface melt and influx entering the ocean. High-resolution stage, discharge and suspended sediment load show a decrease in runoff of ∼25% and in sediment load of ∼40% from 2007 to 2008 in response to a decrease in the summer accumulated number of positive degree-days. During the 2007 and 2008 runoff seasons, jökulhlaups were observed at Kangerlussuaq, drained from an ice-dammed lake at the margin of the GrIS.

Aircraft laser-altimeter surveys during the 1990s showed near-coastal parts of the Greenland ice sheet to be thinning; despite slow thickening at higher elevations, the ice sheet lost mass to the ocean. Many outlet glaciers thinned more rapidly than could be explained by increased melting during the recent warmer summers, indicating dynamic imbalance between glacier velocity and upstream snow accumulation. Results from more recent surveys, presented here, show that thinning rates have increased in most coastal regions. For almost half of the surveys, these increases might have resulted from increases in summer melting, but rapid thinning on others is indicative of dynamic changes that increased with time. In particular, thinning rates on the three fastest glaciers increased to tens of m a−1 after 2000, and other observations show an approximate doubling in their velocities. The deep beds of these glaciers appear to have a strong influence on rates of grounding-line retreat and thickness change, with periods of glacier acceleration and rapid thinning initiated by flotation and break-up of lightly grounded glacier snouts or break-up of floating ice tongues. Near-simultaneous thinning of these widely separated glaciers suggests that warming of deeper ocean waters might be a common cause. Nearby glaciers without deep beds are thinning far more slowly, suggesting that basal lubrication as a result of increased surface melting has only a marginal impact on Greenland outlet-glacier acceleration

Two previous theoretical analyses of data from the GRIP, Vostok and Byrd ice cores, presenting age–depth correlations, grain growth and dislocation-density evolution, are re-examined. It is found that the age–depth correlations are inconsistent with the idealized flow with unchanging history adopted, but that good correlations can be obtained by relaxing those restrictions. A modified grain-growth relation is proposed, consistent with the distinct growth profiles of the Vostok and other two cores, which can be solved simultaneously with the given dislocation-density evolution equation. These are solved for all three cores with the given parameters, and the depth profiles of grain diameter and dislocation density at the present time are determined with the new age–depth correlation and with that shown empirically in the papers. The varying flow history influences the age–depth correlation, and hence the depth profiles, which is important both for the interpretation of core data, and for the determination of constitutive variables at each depth at the present time.