We compare rates of surface-elevation change on the Greenland ice sheet derived from European Remote-sensing Satellite-2 (ERS-2) radar-altimeter data with those obtained from laser-altimeter data collected over nearly the same time periods. Radar-altimeter data show more rapid thickening (9 ± 1 cm a−1 above 1500 m elevation in the north, and 3 ± 1 cm a−1 above 2000 m in the south) than the laser estimates, possibly caused by a lifting of the radar-reflection horizon associated with changes in the snowpack, such as those caused by progressively increased surface melting, as summer temperatures rise. Over all the ice sheet above 2000 m, this results in an ERS-derived volume balance ∼75 ± 15 km3 a−1 more positive than that from laser data. This bias between laser and radar estimates of elevation change varies spatially and temporally, so cannot at present be corrected without independent surveys such as those presented here. At lower elevations, comparison of detailed repeat laser surveys over Jakobshavn Isbræ with ERS results over the same time interval shows substantial ERS underestimation of ice-thinning rates. This results partly from missing data because of ‘bad’ radar waveforms over the very rough surface topography, and partly from the tendency for large radar footprints to sample preferentially local high points in the topography, thus missing regions of most rapid thinning along glacier depressions.

We use an ice-flow model to demonstrate how flow variations initiated in the marginal zone of an ice sheet affect flow farther inland through longitudinal (along-flow) coupling. Our findings allow for an alternate interpretation of seasonal accelerations observed near the equilibrium line of the Greenland ice sheet (Zwally and others, 2002). We demonstrate that these observations can be explained by accelerations initiated up to 12 km closer to the margin where the ice is ∼40% thinner, is heavily crevassed, experiences a seasonal doubling of velocity, and where the ablation rate, surface meltwater flux and ice temperature are likely higher. Our modeling and observations suggest that conditions and processes normally found near ice-sheet margins are adequate for explaining the observations of Zwally and others (2002). This and considerations of the likely subglacial hydrology in the marginal zone lead us to suggest that seasonal accelerations may have limited impact on ice-sheet mass balance even in the face of climate warming.

To assess the potential volume of a glacial lake outburst flood (GLOF) more precisely than in previous studies, we analyze ground survey data and remote-sensing digital elevation models (DEMs) around glacial lakes in the Lunana region, Bhutan. Based on a DEM generated from differential GPS ground surveys, we first evaluate the relative accuracies of DEMs produced by the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and the Shuttle Radar Topography Mission (SRTM). Root-mean-square errors of the altitudinal difference between these DEMs and ground survey data were 11.0 m for ASTER and 11.3 m for SRTM. These errors are similar to those of previous studies. We show that a topographical classification allows a better estimate of elevation on lakes/ponds, riverbeds and glaciers due to their flat surfaces, while the relative accuracy is worse over moraines and hill slopes due to their narrow ridges and steep slopes. Using the optical satellite images and the ground survey data, we re-evaluate the GLOF volume in 1994 as (17.2 ± 5.3) × 106 m3. We show GLOF-related information such as distance, altitudinal difference and gradient at possible outburst points where the lake level is higher than the neighboring riverbed and/or glacial lake.

The age–depth relationship of the Vostok (Antarctica) ice core has been reconstructed in the depth interval 3300–3347 m, by comparing three gas properties in ice (CO2, CH4 and δ18Oatm) with those in the EPICA Dome C (Antarctica) core. Fourteen Vostok depths were examined in this interval, and it was found that nine samples are uniquely dated if candidate ages are restricted to the interval between 400 and 650 kyr. One of these samples is uniquely dated without restriction. The analysis supports previous reports that this section contains ice from Termination V, but that the stratigraphic order of ice is reversed. The top of the overturned layer lies between 3316 and 3319 m. At least one other stratigraphic disturbance was found between 3340 and 3343 m, as indicated by another reversal of the age–depth relationship. Finally, the oldest ice in this section is dated at ≥440 kyr, confirming the existence of ice from the cold marine isotope stage (MIS) 12 interval.

We present techniques for obtaining large (∼100 L STP) samples of ancient air for analysis of 14C of methane (14CH4)and other trace constituents. Paleoatmospheric 14CH4 measurements should constrain the fossil fraction of past methane budgets, as well as provide a definitive test of methane clathrate involvement in large and rapid methane concentration ([CH4]) increases that accompanied rapid warming events during the last deglaciation. Air dating to the Younger Dryas–Preboreal and Oldest Dryas–Bølling abrupt climatic transitions was obtained by melt extraction from old glacial ice outcropping at an ablation margin in West Greenland. The outcropping ice and occluded air were dated using a combination of δ15N of N2, δ18O of O2, δ18Oice and [CH4] measurements. The [CH4] blank of the melt extractions was <4 ppb. Measurements of δ18O and δ15N indicated no significant gas isotopic fractionation from handling. Measured Ar/N2, CFC-11 and CFC-12 in the samples indicated no significant contamination from ambient air. Ar/N2, Kr/Ar and Xe/Ar ratios in the samples were used to quantify effects of gas dissolution during the melt extractions and correct the sample [CH4]. Corrected [CH4] is elevated over expected values by up to 132 ppb for most samples, suggesting some in situ CH4 production in ice at this site.

A 5 year record of data from an automatic weather station (AWS) operating in the ablation zone of Storbreen, Norway, has been used to calculate the local surface energy and mass balance. The AWS observations cover five mass-balance years with an unusually strong mass deficit on Storbreen. The average energy flux (Q) contributing to melt for the period 2001–06 is 113 W m−2. Of this, the net shortwave radiation flux is the dominant contributor (92 W m−2), followed by the sensible heat flux (20 W m−2) and the latent heat flux (9 W m−2). The net longwave radiation (–6 W m−2) and the subsurface heat flux (–2 W m−2) contribute negatively to the budget. Net radiation thus produces 76% of the melt, while the turbulent fluxes and the subsurface heat flux produce 24% of the total melt. The seasonal mean incoming shortwave radiation is remarkably constant between the years, whereas variations in temperature and reflected shortwave radiation (albedo) explain most of the interannual variation in melt. The modelled ablation compares well with the measured ablation from stake readings. The sensitivity of the energy-balance model was examined by varying the surface roughness length of momentum and the sensitivity of the calculated melt by perturbations of temperature, wind speed and relative humidity.

The ionic and isotopic characteristics of bulk waters emanating from the cold-based Longyearbreen, central Svalbard, in 2004 are examined to determine lithological, hydrological and glaciological controls on water composition, solute provenance and chemical denudation. The geology consisted of reactive coal seams and associated sedimentary rocks. Acidity caused by microbial-mediated oxidation of sulfides and, to a lesser extent, nitrogen-bearing minerals was neutralized by congruent dissolution of dolomite and incongruent weathering of silicates in open-system subglacial drainage channels. The ablation season was divided into an early melt season, a peak-flow period and a late melt season. The runoff distribution during these periods was 1.7%, 89.7% and 8.6%, respectively, whereas the solute flux distribution was 1.9%, 82.1% and 16.0%, respectively. Comparisons between different annual solute flux estimation methods indicated that extrapolation of peak-flow period data significantly underestimated both the early- and late-melt-season solute fluxes. About 3.8% of the solutes derived from sea-salt spray, 0.7% from acid aerosol deposition and 95.5% from crustal/organic sources. The physical and chemical conditions resulted in diffusion of CO2 rather than atmospheric drawdown. The cation-equivalent weathering rate and the crustal solute yield were 322 ΣmEq+m−2 a−1 and 22 t km−2 a−1, respectively, which are within the regional range of Svalbard. However, the chemical weathering intensity was as high as 940 ΣmEq+ m−3 owing to the relatively low specific discharge of 0.34 m a−1.

We recorded the burial times of temperature sensors mounted on a specially constructed tower to determine snow accumulation during individual storms in the summit caldera of Mount Wrangell, Alaska, USA, (62° N, 144° W; 4100 m a.s.l.) during the accumulation year June 2005 to June 2006. The experiment showed most of the accumulation occurred in episodic large storms, and half of the total accumulation was delivered in late summer. The timing of individual events correlated well with storms recorded upwind, at Cordova, the closest Pacific coastal weather station (200 km south-southeast), although the magnitude of events showed only poor correlation. Hence, snow accumulation at Mount Wrangell appears to be a reflection of synoptic-scale regional weather systems. The accumulation at Mount Wrangell’s summit (>2.5 m w.e.) exceeded the precipitation at Cordova. Although the direct relationship between accumulation of individual storms at the summit of Mount Wrangell and precipitation events at Cordova may be unique in the region, it is useful for interpreting ice cores obtained on Mount Wrangell. This is especially the case here because the high rate of accumulation allows high time resolution within the core.

Sea ice surviving the summer melt season to become multi-year ice in the Arctic Ocean is of interest because multi-year ice significantly affects the ice-thickness distribution and the dynamics and thermodynamics of the ice pack in subsequent seasons. However, the amount of ice surviving summer melting has not been well determined because the time of the minimum ice area varies from region to region. A concept of local temporal minimum (LTM) accounts for non-simultaneity of the melt–freeze transition by determining the minima ice concentrations (CLTM) on local spatial scales. CLTM are calculated for 25 km gridcells using 24 years (1979–2002) of satellite passive-microwave data. The total area of ice surviving the summer melt (ALTM) is given by spatial integration of CLTM. Over 24 years, the average ALTM is 2.6 × 106 km2 (excluding ∼0.7 × 105 km2 above 84° N). In contrast, the average area (3.8 × 106 km2) of all ice types (ASM), measured when the total (simultaneous) ice cover is a minimum in daily maps in mid-September, is an often-used estimate of ice surviving the summer melting that is ∼45% too large. Over 24 years, the ALTM decreased by 9.5 ± 2.2% (10 a)−1 (0.27 ± 0.06 × 106 km2 (10 a)−1), which is similar to the rate of decline of ASM and about three times the rate of the annual average. The time-of-occurrence of the LTM averaged over the perennial ice pack increased by 8 days from around 11 to 19 August, indicating a later ending of the melt season by about 3 days (10 a)−1 as the summer pack declines. Estimates of multi-year ice in midwinter from passive microwave observations are ∼17% smaller than ALTM, suggesting that the microwave algorithm does not measure all the multi-year ice.

Water levels were measured in boreholes spaced along the entire length of Bench Glacier, Alaska, USA, for a period in excess of 2 years. Instrumented boreholes were arranged as nine pairs along the center line of the glacier and an orthogonal grid of 16 boreholes in a 3600 m2 region at the center of the ablation area. Diurnal fluctuations of the water levels were found to be restricted to the late melt season. Pairs of boreholes spaced along the length of the ablation area often exhibited similar fluctuations and diurnal changes in water levels. Three distinct and independent types of diurnal fluctuations in water level were observed in clusters of boreholes within the grid of boreholes. Head gradients suggest water did not flow between clusters, and a single tunnel connecting the boreholes could not explain the observed pattern of diurnal water-level fluctuations. Inter-borehole and borehole-cluster connectivity suggests the cross-glacier width of influence of a segment of the drainage system connected to a borehole was limited to tens of meters. A drainage configuration whereby boreholes are connected to a somewhat distant tunnel by drainage pipes of differing lengths, often hundreds of meters, is shown with a numerical test to be a plausible explanation for the observed borehole behavior.

Annual equilibrium-line altitude (ELA) and surface mass balance of Glacier Blanc, Ecrins region, French Alps, were reconstructed from a 25 year time series of satellite images (1981–2005). The remote-sensing method used was based on identification of the snowline, which is easy to discern on optical satellite images taken at the end of the ablation season. In addition, surface mass balances at the ELA were reconstructed for the same period using meteorological data from three nearby weather stations. A comparison of the two types of series reveals a correlation of r > 0.67 at the 0.01 level of significance. Furthermore, the surface mass balances obtained from remote-sensing data are consistent with those obtained from field measurements on five other French glaciers (r = 0.76, p < 0.01). Also consistent for Glacier Blanc is the total mass loss (10.8 m w.e.) over the studied period. However, the surface mass balances obtained with the remote-sensing method show lower interannual variability. Given that the remote-sensing method is based on changes in the ELA, this difference probably results from the lower sensitivity of the surface mass balance to climate parameters at the ELA.

Spatio-temporal variations of the recently determined accumulation rate are investigated using ground-penetrating radar (GPR) measurements and firn-core studies. The study area is located on Ritscherflya in western Dronning Maud Land, Antarctica, at an elevation range 1400–1560 m. Accumulation rates are derived from internal reflection horizons (IRHs), tracked with GPR, which are connected to a dated firn core. GPR-derived internal layer depths show small relief along a 22 km profile on an ice flowline. Average accumulation rates are about 190 kg m−2 a−1 (1980–2005) with spatial variability (1σ) of 5% along the GPR profile. The interannual variability obtained from four dated firn cores is one order of magnitude higher, showing 1σ standard deviations around 30%. Mean temporal variations of GPRderived accumulation rates are of the same magnitude or even higher than spatial variations. Temporal differences between 1980–90 and 1990–2005, obtained from two dated IRHs along the GPR profile, indicate temporally non-stationary processes, linked to spatial variations. Comparison with similarly obtained accumulation data from another coastal area in central Dronning Maud Land confirms this observation. Our results contribute to understanding spatio-temporal variations of the accumulation processes, necessary for the validation of satellite data (e.g. altimetry studies and gravity missions such as Gravity Recovery and Climate Experiment (GRACE)).

Experiments have been conducted on the French full-scale experimental site at Lautaret pass to improve our understanding of the action of snow avalanches on obstacles. The ultimate objective is to provide realistic pressure distribution models suitable for use in civil engineering design and to eliminate the restrictive assumptions currently used in this field. We focus on the feasibility of using the inverse method to quantify the action of the avalanche from its effects on realistic structures rather than from sensors placed directly in the flow. This approach takes into account the interactions between the flow and the obstacle and ensures that the result is effectively the action experienced by the obstacle. The inverse analysis procedure is developed and validated using both numerical and laboratory tests. In situ tests carried out at the Lautaret site to determine the avalanche action at different scales confirm the reliability of this original approach. Its intrinsic characteristics make it especially suitable for application to different structures to provide new knowledge in this complex field.

We investigate snowpack properties at a site in west-central Greenland with ground-penetrating radar (GPR), supplemented by stratigraphic records from snow pits and shallow firn cores. GPR data were collected at a validation test site for CryoSat (T05 on the Expéditions Glaciologiques Internationales au Groenland (EGIG) line) over a 100 m × 100 m grid and along 1 km sections at frequencies of 500 and 800 MHz. Several internal reflection horizons (IRHs) down to a depth of 10 m were tracked. IRHs are usually related to ice-layer clusters in vertically bounded sequences that obtain their initial characteristics near the surface during the melt season. Warm conditions in the following melt season can change these characteristics by percolating meltwater. In cold conditions, smaller melt volumes at the surface can lead to faint IRHs. The absence of simple mechanisms for internal layer origin emphasizes the need for independent dating to reliably interpret remotely sensed radar data. Our GPR-derived depth of the 2003 summer surface of 1.48 m (measured in 2004) is confirmed by snow-pit observations. The distribution of IRH depths on a 1 km scale reveals a gradient of increasing accumulation to the northeast of about 5 cm w.e. km−1. We find that point measurements of accumulation in this area are representative only over several hundred metres, with uncertainties of about 15% of the spatial mean.

Annual-layer thickness data, spanning AD 1534–2001, from an ice core from East Rongbuk Col on Qomolangma (Mount Everest, Himalaya) yield an age–depth profile that deviates systematically from a constant accumulation-rate analytical model. The profile clearly shows that the mean accumulation rate has changed every 50–100 years. A numerical model was developed to determine the magnitude of these multi-decadal-scale rates. The model was used to obtain a time series of annual accumulation. The mean annual accumulation rate decreased from ∼0.8 m ice equivalent in the 1500s to ∼0.3 m in the mid-1800s. From ∼1880 to ∼1970 the rate increased. However, it has decreased since ∼1970. Comparison with six other records from the Himalaya and the Tibetan Plateau shows that the changes in accumulation in East Rongbuk Col are broadly consistent with a regional pattern over much of the Plateau. This suggests that there may be an overarching mechanism controlling precipitation and mass balance over this area. However, a record from Dasuopu, only 125 km northwest of Qomolangma and 700 m higher than East Rongbuk Col, shows a maximum in accumulation during the 1800s, a time during which the East Rongbuk Col and Tibetan Plateau ice-core and tree-ring records show a minimum. This asynchroneity may be due to altitudinal or seasonal differences in monsoon versus westerly moisture sources or complex mountain meteorology.

Despite the large amount of subglacial lakes present underneath the East Antarctic ice sheet and the melt processes involved, the hydrology beneath the ice sheet is poorly understood. Changes in subglacial potential gradients may lead to subglacial lake outbursts, discharging excess water through a subglacial drainage system underneath the ice sheet. Such processes can eventually lead to an increase in ice flow. In this paper, a full Stokes numerical ice-sheet model was employed which takes into account the ice flow over subglacial water bodies in hydrostatic equilibrium with the overlying ice. Sensitivity experiments were carried out for small perturbations in ice flow and basal melt rate as a function of ice thickness, general surface slope, ice viscosity and lake size, in order to investigate their influence on the subglacial potential gradient and the impact on subglacial lake drainage. Experiments clearly demonstrate that small changes in surface slope are sufficient to start and sustain episodic subglacial drainage events. Lake drainage can therefore be regarded as a common feature of the subglacial hydrological system and may influence, to a large extent, the present and future behavior of large ice sheets.

In order to examine the effects of solutes on recrystallization and subsequent grain growth in ice, both doped and undoped ice single crystals were extruded through a 120° equal-channel angular extrusion jig, in order to impart a large shear strain (∼1.15). Upon subsequent annealing at −3°C, the original single crystals recrystallized, in most cases to a new single crystal with a different orientation. Increasing the solute concentration (for H2SO4 to ∼200–300 ppb, and for NaCl, KCl and MgSO4 to >5 ppm) was found to significantly retard the growth and possibly, for H2SO4-doped ice, the nucleation of new grains in the strained ice single crystals. This is indicative of how soluble impurities can retard grain growth in ice cores. It was also found that the migrating grain boundaries surrounding the newly formed grains contained large concentrations of impurities, often observed as filaments. These could have formed by the grain boundaries sweeping up impurities from the lattice into the boundary or by their diffusion to the boundary – mechanisms whereby impurities could be concentrated into the grain boundaries in ice cores – although the latter mechanism seems unlikely since it would require very high diffusion rates.

During 2000–07, five giant icebergs (B15A, B15J, B15K, C16 and C25) adrift in the southwestern Ross Sea, Antarctica, were instrumented with global positioning system (GPS) receivers and other instruments to monitor their behavior in the near-coastal environment. The measurements show that collision processes can strongly influence iceberg behavior and delay their progress in drifting to the open ocean. Collisions appear to have been a dominant control on the movement of B15A, the largest of the icebergs, during the 4-year period it gyrated within the limited confines of Ross Island, the fixed Ross Ice Shelf and grounded C16. Iceberg interactions in the near-coastal regime are largely driven by ocean tidal effects which determine the magnitude of forces generated during collision and break-up events. Estimates of forces derived from the observed drift trajectories during the iceberg-collisioninduced calving of iceberg C19 from the Ross Ice Shelf, during the iceberg-induced break-off of the tip of the Drygalski Ice Tongue and the break-up of B15A provide a crude estimate of the stress scale involved in iceberg calving. Considering the total area the vertical face of new rifts created in the calving or break-up process, and not accounting for local stress amplification near rift tips, this estimated stress scale is 104 Pa.