Seismometers deployed over a 3 year period on icebergs in the Ross Sea and on the Ross Ice Shelf, Antarctica, reveal that impulsive sources of ocean surface waves are frequent (e.g. ∼200 events per year in the Ross Sea) in the ice-shelf and iceberg-covered environment of coastal Antarctica. The 368 events recorded by our field deployment suggest that these impulsive events are generated by glaciological mechanisms, such as (1) small-scale calving and edge wasting of icebergs and ice-shelf fronts, (2) edge-on-edge closing and opening associated with iceberg collisions and (3) possibly the impulsive opening of void space associated with ice-shelf rifting and basal crevasse formation. The observations described here provide a background of glaciogenic ocean-wave phenomena relevant to the Ross Sea and suggest that these phenomena may be exploited in the future (using more purposefully designed observation schemes) to understand iceberg calving and ice-shelf disintegration processes.

To contribute to the understanding of the influence of climate on glacial erosion and on orogenic processes, we report contemporary glacial erosion rates from a calving glacier in the Southern Andes and elucidate the influence of ice dynamics on erosion. Using seismic profiles of sediments collected in a proglacial fjord and a documented history of retreat, we determine the time-varying sediment flux of Glaciar Marinelli as a measure of basin-wide erosion rates, and compare these rates with the annual ice budget reconstructed using NCEP–NCAR reanalysis climate data from 1950 to 2005. The rate of erosion of the largest tidewater glacier in Tierra del Fuego averaged 39 ± 16 mm a−1 during the latter half of the 20th century, with an annual maximum approaching 130 mm a−1 following a decade of rapid retreat. A strong correlation emerges between the variable rate of ice delivery to the terminus and the erosion rate, providing quantitative insight into the relationship between ice fluxes and glacial erosion rates. For Glaciar Marinelli, as for other calving glaciers for which suitable data exist, the marked retreat and thinning over the past 50 years have resulted in a period of accelerated basal sliding and unusually rapid erosion.

The glaciers on Yulong mountain (5596 m a.s.l.), China, are the southernmost in mainland Eurasia. The largest is Baishui glacier No. 1. Up to 90% of the annual precipitation there falls in summer, when warm moisture-rich air masses associated with the southwest monsoon reach the area from the Indian Ocean. The winter climate is influenced by air masses with a continental origin and by the southern branch of the westerlies. The snow that accumulates on Baishui glacier No. 1 includes marine aerosols associated with the summer monsoon, and dust brought in winter from central/west Asia, Africa or the Thar Desert area. Studies in May 2006 at two sites, one in the accumulation area (4900 m) and one in the ablation area (4750 m), revealed differences between the ionic composition of the snow that had accumulated in the 2005/06 winter and that of the snow which had been deposited during the preceding summer monsoon. Differences in the chemistry of the summer-accumulated snow at the two sites probably reflected local differences in ablation and elution rates, rather than differences in ion supplies. Differences in the chemistry of the winter-accumulated snow may reflect the influence of up-valley winds, which bring more crustal material to the lower site.

No single flow model can simulate all possible glaciers and ice sheets without violating fundamental tenets of computational science. The root cause is not one of numerical sophistication, precision or accurate initial conditions. Instead, using flow and transport as data transmission, glaciers inadvertently function as information processors. This computational capability confers a level of complexity that inherently limits our ability to accurately and efficiently predict glacier flow and therefore, for example, to forecast those aspects of climate systems that depend on glaciers. In particular, even with considerable future advancements in glacier physics, computational theory shows that no dramatic improvements in numerical speed are likely when compared to today’s glacier models. Therefore, to increase speed and resolution, the next generation of climate and sea-level models must rely on simulations tailored to specific ice-sheet geometries rather than general-purpose glacier flow models. However, because glaciers process information, entirely new computation-theoretic advances in glaciology are possible, and concepts from information entropy may help to define new glacier scaling relationships and identify which geometries will be most problematic for modeling.

This paper presents a newly derived velocity field of the Brunt Ice Shelf and Stancomb-Wills Ice Tongue in East Antarctica. Velocity is derived using automated feature-tracking techniques on RADARSAT-1 imagery from the 1997 Antarctic Mapping Mission-1 (AMM-1) and 2000 Modified Antarctic Mapping Mission (MAMM). We measure velocities up to 1350 m a−1 at the ice front. Closer investigation of the velocity field and radar images reveals an obvious divergence of contemporary flowlines from relict flow stripes along a portion of Stancomb-Wills Ice Tongue. Analysis of relict flow stripes can reveal glacial events in the past. Assuming relatively steady flow since the event occurred, we utilize the current velocity field to reconstruct former stripe configuration. We believe the flowline deviation points to a change in dynamics that occurred approximately 140 years ago, perhaps due to thinning and an associated shift in grounding-line position. This seems to be supported by historic observations and bottom topography data.

We use models constrained by remotely sensed data from Pine Island and Thwaites Glaciers, West Antarctica, to infer basal properties that are difficult to observe directly. The results indicate strong basal melting in areas upstream of the grounding lines of both glaciers, where the ice flow is fast and the basal shear stress is large. Farther inland, we find that both glaciers have ‘mixed’ bed conditions, with extensive areas of both bedrock and weak till. In particular, there are weak areas along much of Pine Island Glacier’s main trunk that could prove unstable if it retreats past the band of strong bed just above its current grounding line. In agreement with earlier studies, our forward ice-stream model shows a strong sensitivity to small perturbations in the grounding line position. These results also reveal a large sensitivity to the assumed bed (sliding or deforming) model, with non-linear sliding laws producing substantially greater dynamic response than earlier simulations that assume a linear-viscous till rheology. Finally, comparison indicates that our results using a plastic bed are compatible with the limited observational constraints and theoretical work that suggests an upper bound exists on maximum basal shear stress.

We investigate the transferability of an enhanced temperature-index melt model that was developed and tested on Haut Glacier d’Arolla, Switzerland, in the 2001 season. The model’s empirical parameters (temperature factor, TF, and shortwave radiation factor, SRF) are recalibrated for: (1) other locations on Haut Glacier d’Arolla; (2) subperiods of distinct meteorological conditions; (3) different years on Haut Glacier d’Arolla; and (4) other glaciers in different years. The model parameters are optimized against simulations of an energy-balance model validated against ablation observations. Results are compared with those obtained with the original parameters. The model works very well when applied to other sites, seasons and glaciers, with the exception of overcast conditions. Differences are due to underestimation of high melt rates. The parameter values are associated with the prevailing energy-balance conditions, showing that high SRF are obtained on clear-sky days, whereas higher TF are typical of locations where glacier winds prevail and turbulent fluxes are high. We also provide a range of parameters clearly associated with the site’s location and its meteorological characteristics that could help to assign parameter values to sites where few data are available.

In summer 2005, two pilot snow/firn cores were obtained at 5365 and 5206 m a.s.l. on Fedchenko glacier, Pamirs, Tajikistan, the world’s longest and deepest alpine glacier. The well-defined seasonal layering appearing in stable-isotope and trace element distribution identified the physical links controlling the climate and aerosol concentration signals. Air temperature and humidity/precipitation were the primary determinants of stable-isotope ratios. Most precipitation over the Pamirs originated in the Atlantic. In summer, water vapor was re-evaporated from semi-arid regions in central Eurasia. The semi-arid regions contribute to non-soluble aerosol loading in snow accumulated on Fedchenko glacier. In the Pamir core, concentrations of rare earth elements, major and other elements were less than those in the Tien Shan but greater than those in Antarctica, Greenland, the Alps and the Altai. The content of heavy metals in the Fedchenko cores is 2–14 times lower than in the Altai glaciers. Loess from Afghan–Tajik deposits is the predominant lithogenic material transported to the Pamirs. Trace elements generally showed that aerosol concentration tended to increase on the windward slopes during dust storms but tended to decrease with altitude under clear conditions. The trace element profile documented one of the most severe droughts in the 20th century.

Broadband radiation schemes (parameterizations) are commonly used tools in glacier mass-balance modelling, but their performance at high altitude in the tropics has not been evaluated in detail. Here we take advantage of a high-quality 2 year record of global radiation (G ) and incoming longwave radiation (L ↓) measured on Kersten Glacier, Kilimanjaro, East Africa, at 5873 m a.s.l., to optimize parameterizations of G and L ↓. We show that the two radiation terms can be related by an effective cloud-cover fraction neff , so G or L ↓ can be modelled based on neff derived from measured L ↓ or G, respectively. At neff = 1, G is reduced to 35% of clear-sky G, and L ↓ increases by 45–65% (depending on altitude) relative to clear-sky L ↓. Validation for a 1 year dataset of G and L ↓ obtained at 4850 m on Glaciar Artesonraju, Peruvian Andes, yields a satisfactory performance of the radiation scheme. Whether this performance is acceptable for mass-balance studies of tropical glaciers is explored by applying the data from Glaciar Artesonraju to a physically based mass-balance model, which requires, among others, G and L ↓ as forcing variables. Uncertainties in modelled mass balance introduced by the radiation parameterizations do not exceed those that can be caused by errors in the radiation measurements. Hence, this paper provides a tool for inclusion in spatially distributed mass-balance modelling of tropical glaciers and/or extension of radiation data when only G or L ↓ is measured.

We examine patterns of localized surface elevation change in lower Mercer and Whillans Ice Streams, West Antarctica, which we interpret as subglacial water movement through a system of lakes and channels. We detect and measure the lake activity using repeat-track laser altimetry from ICESat and image differencing from MODIS image pairs. A hydrostatic-potential map for the region shows that the lakes are distributed across three distinct hydrologic regimes. Our analysis shows that, within these regimes, some of the subglacial lakes appear to be linked, with drainage events in one reservoir causing filling and follow-on drainage in adjacent lakes. We also observe changes near ice raft ‘a’ in lower Whillans Ice Stream, and interpret them as evidence of subglacial water and other changes at the bed. The study provides quantitative information about the properties of this complex subglacial hydrologic system, and a relatively unstudied component of ice-sheet mass balance: subglacial drainage across the grounding line.

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.

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.

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.

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.

Glacier runoff from the Yangtze River source region (YRSR), China, is estimated for the period 1961–2000 using a degree-day approach. In the investigation area, glacier runoff accounts for 11.0% of the total river runoff during the period 1961–2000. In the 1990s its contribution to river runoff rises to 17.0%. Due to the current rate of glacier decline, the impact of glacier runoff on river runoff has recently increased in the source region. Based on two different climate-change scenarios derived from ECHAM5/MPI-OM, future glacier runoff is assessed for the period 2001–50. In all climate-change scenarios, annual glacier runoff shows a significant increase due to intensified ice melting. There is an increase in glacier runoff during spring and early summer, yet a significant decrease in late summer. This study highlights the current and future impact of glacier runoff on river runoff in the YRSR.

This paper examines the basal ice sequence of Rhone Glacier, a cold-based glacier in the McMurdo Dry Valleys, Antarctica, using isotopic and solute chemistry data. Three different ice facies are identified: englacial, amber and stratified. The englacial facies is clean, bubbly ice of meteoric origin and is underlain by an amber facies. Amber ice is a characteristic of cold alpine glaciers in the McMurdo Dry Valleys and is distinctive for its high solute concentrations and much higher strain rates compared with the overlying englacial ice and the underlying stratified ice. Analysis of the stratified facies reveals an isotopic signature indicative of melt then refreeze processes and it is most likely associated with apron entrainment at the margin. By contrast, the amber facies has a co-isotopic slope of 8 and plots on a meteoric waterline. The inclusion of impurities in the amber ice reveals prolonged contact with the bed, and its depleted isotopic signature is consistent with ice formed during a cooler period. Comparison of the basal sequence of Rhone Glacier with other cold-based glaciers in the McMurdo Dry Valleys reveals strong similarities between valley-side glaciers (e.g. Meserve and Rhone Glaciers), whereas valley-floor glacier basal sequences (e.g. Suess Glacier) are characterized by structurally complex amalgamations of ice and debris.

We investigate the major sources of backscatter at 5.3 GHz, within the superimposed ice and firn areas of a polythermal glacier. Two ground-penetrating radar systems, an 800 MHz impulse system and a polarimetric 5.3 GHz frequency-modulated continuous-wave system, are used to acquire along-glacier profiles in the accumulation area of Kongsvegen, Svalbard. The 800 MHz response is used to map reflection horizons in the glacier. Using cores from the superimposed ice and firn areas, the causes of these reflection horizons, in terms of snow, firn and ice layers, are investigated. Superimposing the reflection horizons on the co-polarized and cross-polarized 5.3 GHz profile, we are able to determine how the 5.3 GHz frequency responds to the different media. Scattering at rough interfaces and volume scattering occur in the superimposed ice area and are apparently caused by air-bubble number, size and distribution. In the firn the strongest return originates from below the previous summer surface, consistent with previous findings. At approximately the same depth, strong incoherent scattering begins. The rapid decrease in coherent reflections indicates the significance of scattering in the firn.