The longest continuous glaciological mass-balance time-series in the intertropical zone of South America goes back to 1991 on Glaciar Zongo, Bolivia. Photogrammetric and hydrological data have been used to (1) check the specific net balance over long periods and (2) extend the mass-balance time series over the last 50 years. These data reveal a bias in the glaciological mass balance which can be explained by the field-measurement sampling network. Our study shows a large temporal variability of the surface mass balances in the ablation area and reveals strong relationships between independent surface mass-balance data coming from selected ablation areas with numerous data. It demonstrates the very large contribution (80%) of low-elevation ranges (one-third of the surface) to the specific mass balance and, consequently, the importance of the reduction of the area of the tongue. With these new results, Glaciar Zongo offers the longest and most accurate mass-balance series in any Andean country. The dataset shows that Glaciar Zongo experienced a relatively steady state over the period 1956–75, with even a slight mass gain over 1963–75, and a rapid and continuous decrease since then.

Where polythermal glaciers have frozen margins that buttress otherwise temperate-based sliding ice, longitudinal compression can strongly influence ice-flow trajectory, and consequently sediment transport paths. Past efforts to model flow in the vicinity of a basal thermal transition (BTT) have generally relied on simplified boundary conditions or rheological idealizations, making these model results difficult to apply to real glacier termini. Herein, we present results of numerical simulations using a power-law rheology and with boundary conditions that better represent the frozen margin. Model results indicate that a transition to a non-sliding frozen margin causes a decline in surface velocity made possible by upward ice flow, implying either enhanced ablation for steady-state simulations or the formation of a surface bulge. Permitting ice loss by ablation combined with numerical smoothing of the basal slip transition subdues basal stress concentrations and thereby inhibits development of structural discontinuities such as thrust faults. Upward ice flow is accommodated by vertical extension up-glacier of the BTT. This strain regime can potentially account for key structural features in polythermal glacier termini without appealing to thrusting.

Following recommendations from the International Commission on Snow and Ice for a world glacier inventory, an inventory of glaciers in China was carried out by Chinese glaciologists from 1978 to 2002. Each glacier was measured from aerial photographs and topographical maps and 34 parameters recorded. These parameters were then analyzed statistically for the various river systems in China. Twelve volumes of the Glacier Inventory of China (GIC) have been published, consisting of 22 parts in 21 books. The data were subsequently abridged into a Concise GIC, published in Chinese (2005) and in English (2008), to make the glacier inventory more accessible and better adapted for assessing glacier response to climate change. After the GIC was completed, new aerial photographs became available and remote-sensing techniques became more common. To investigate glacier changes since completion of the first GIC, a second Glacier Inventory of China was initiated in 2007. This 5 year project, supported by the Ministry of Science and Technology, will be undertaken mainly using remote-sensing techniques.

We conducted two-dimensional continuous multi-offset georadar surveys on Bench Glacier, south-central Alaska, USA, to measure the distribution of englacial water. We acquired data with a multichannel 25 MHz radar system using transmitter–receiver offsets ranging from 5 to 150 m. We towed the radar system at 5–10 kmh–1 with a snow machine with transmitter/receiver positions established by geodetic-grade kinematic differentially corrected GPS (nominal 0.5 m trace spacing). For radar velocity analyses, we employed reflection tomography in the pre-stack depth-migrated domain to attain an estimated 2% velocity uncertainty when averaged over three to five wavelengths. We estimated water content from the velocity structure using the complex refractive index method equation and use a three-phase model (ice, water, air) that accounts for compression of air bubbles as a function of depth. Our analysis produced laterally continuous profiles of glacier water content over several kilometers. These profiles show a laterally variable, stratified velocity structure with a low-water-content (~0–0.5%) shallow layer (~20–30 m) underlain by high-water-content (1–2.5%) ice.

A number of glaciers in the Canadian High Arctic are composed primarily of cold ice, but the ice at or near their beds reaches the pressure-melting point (PMP) in the ablation zone. Past modelling studies have suggested that the basal temperatures of some of these glaciers reach the PMP where they should not, indicating that they are not in thermal equilibrium with present-day surface air temperatures. To investigate the possible reasons for thermal disequilibria in such glaciers, a two-dimensional ice temperature model was used to simulate the inferred thermal characteristics of John Evans Glacier, Ellesmere Island. Results indicate that while surface refreezing and historical ice-thickness changes have had a warming effect upon basal ice temperatures, supraglacial meltwater reaching the glacier bed provides the single most critical heat source for explaining the apparent thermal disequilibrium between present-day inferred ice–bed temperatures and those modelled under present-day boundary conditions.

Sensible-heat flux is obviously important for glacier ablation but is difficult to measure routinely. Sensible-heat flux can be estimated from wind-speed and temperature data using a dimensionless heat-transfer coefficient. Values of the heat-transfer coefficient are evaluated for six sites by correlating measured melt energy with a wind–temperature variable (product of daily mean wind speed, temperature and mean atmospheric pressure for the altitude in question). Data are available for short periods from two sites in Arctic Canada and two sites in North Greenland, and for hundreds of days of record at Nordbogletscher and Qamanârssûp sermia in South and West Greenland, respectively. Average transfer coefficients for four out of the six sites are close to 0.003, which is in reasonable agreement with values reported elsewhere, while larger values of 0.0047 and 0.0057 are found at the other two sites. Heat-transfer coefficients are also estimated on a monthly basis for the two long records, and substantial variations are found, suggesting that the method should not be used for <20–30 days of data. The present study is based on manually observed ablation and climate data, but the approach could be updated to use data from automatic recording stations using modern sensors.

Radar profiles of bed echo intensity can survey conditions at the ice–bed interface and test for the presence or absence of water. However, extracting information about basal conditions from bed echo intensities requires an estimate of the attenuation loss through the ice. We used the relationship between bed echo intensities from constant-offset radar data and ice thickness to estimate depth-averaged attenuation rates at several locations on and near Kamb Ice Stream (KIS), West Antarctica. We found values varying from 29 dBkm–1 at Siple Dome to 15 dBkm–1 in the main trunk region of KIS, in agreement with a previous measurement and models. Using these attenuation-rate values, we calculated the relative bed reflectivity throughout our KIS surveys and found that most of the bed in the trunk has high basal reflectivities, similar to those obtained in the location of boreholes that found water at the bed. Areas of lower bed reflectivity are limited to the sticky spot, where a borehole found a dry bed, and along the margins of KIS. These results support previous hypotheses that the basal conditions at locations like the sticky spot on KIS control its stagnation and possible reactivation.

This paper reviews the origins of the World Glacier Inventory program and Canada’s participation. It tracks the program of the Canadian federal government from its initiation to its demise. Today, progress is geared towards the international GLIMS and GlobGlacier programs, to which scientists from several Canadian universities are making contributions. Achievements during the past decade, towards providing data on Canada’s glacierized areas in the Arctic and in western Canada, are documented.

We discuss the inclusion of the subsurface heat-conduction flux into the calculation of the energy balance and ablation at the glacier–atmosphere interface. Data from automatic weather stations are used to force an energy-balance model at several locations on alpine glaciers and at one site in the dry Andes of central Chile. The heat-conduction flux is computed using a two-layer scheme, assuming that 36% of the net shortwave radiation is absorbed by the surface layer and that the rest penetrates into the snowpack. We compare simulations conducted with and without subsurface heat flux. Results show that assuming a surface temperature of zero degrees leads to a larger overestimation of melt at the sites in the accumulation area (10.4–13.3%) than in the ablation area (0.5–2.8%), due to lower air temperatures and the presence of snow. The difference between simulations with and without heat conduction is also high at the beginning and end of the ablation season (up to 29% for the first 15 days of the season), when air temperatures are lower and snow covers the glacier surface, while they are of little importance during periods of sustained melt at all the locations investigated.

The quantitative assessment of glacier changes as well as improved modeling of climate-change impacts on glaciers requires digital vector outlines of individual glacier entities. Unfortunately, such a glacier inventory is still lacking in many remote but extensively glacierized gions such as the Canadian Arctic. Multispectral satellite data in combination with digital elevation models (DEMs) a particularly useful for creating detailed glacier inventory data including topographic information for each entity. In this study, we extracted glacier outlines and a DEM using two adjacent Terra ASTER scenes acquired in August 2000 for a remote region on southern Baffin Island, Canada. Additionally, Little Ice Age (LIA) extents we digitized from trimlines and moraines visible on the ASTER scenes, and Landsat MSS and TM scenes from the years 1975 and 1990 we used to assess changes in glacier length and area. Because automated delineation of glaciers is based on a band in the shortwave infrared, we have developed a new semi-automated glacier-mapping approach for the MSS sensor. Wrongly classified debris-coved glaciers, water bodies and attached snowfields we corrected manually for both ASTER and MSS. Glacier drainage divides we manually digitized by combining visual interptation with DEM information. In this first paper, we describe the applied methods for glacier mapping and the glaciological challenges encounted (e.g. data voids, snow cover, ice caps, tributaries), while the second paper ports the data analyses and the derived changes.

We present the results of several radio-echo sounding surveys carried out on Johnsons and Hurd Glaciers, Livingston Island, Antarctica, between the 1999/2000 and 2004/05 austral summer campaigns, which included both radar profiling and common-midpoint measurements with low (20– 25 MHz)- and high (200 MHz)-frequency radars. The latter have allowed us to estimate the radio-wave velocity in ice and firn and the corresponding water contents in temperate ice, which vary between 0 and 1.6% depending on the zone. Maximum ice thickness is ~200 m, with a mean value of 93.6±2.5 m. Total ice volume is 0.968±0.026 km3, for an area of 10.34±0.03km2. The subglacial relief of Johnsons Glacier is quite smooth, while that of Hurd Glacier shows numerous overdeepenings and peaks. The radar records suggest that Hurd Glacier has a polythermal structure, contrary to the usual assumption that glaciers in Livingston Island are temperate. This is also supported by other dynamical and geomorphological evidence.

The geomorphology of palaeo-ice-stream beds and the internal structure of underlying tills can provide important information about the subglacial conditions during periods of fast flow and quiescence. This paper presents observations from three-dimensional seismic data, revealing the geomorphology of buried beds of the Bjørnøyrenna (Bear Island Trough) ice stream, the main drainage outlet of the former Barents Sea ice sheet. Repeated changes in ice dynamics are inferred from the observed successions of geomorphic features. Megablocks, aligned in long chains parallel to inferred ice-stream flowlines, and forming dipping plates that are thrust one on top of another, are taken as evidence for conditions of compressive ice flow. Mega-scale glacial lineations (MSGL) and pull-apart of underlying sediment blocks suggest extensional flow. The observed pattern of megablocks and rafts overprinted by MSGL indicates a change in ice dynamics from a compressional to an extensional flow regime. Till stiffening, due to subglacial freezing, is the favoured mechanism for creating switches in sub-ice-stream conditions. The observed pattern of geomorphic features indicates that periods of slowdown or quiescence were commonly followed by reactivation and fast flow during several glaciations, suggesting that this may be a common behaviour of marine ice streams.

Despite its large area covered by glaciers and ice caps, detailed glacier inventory data are not yet available for most parts of Baffin Island, Canada. Automated classification of satellite data could help to overcome the data gaps. Along-track stereo sensors allow the derivation of a digital elevation model (DEM) and glacier outlines from the same point in time, and are particularly useful for this task. While part I of this study describes the remote-sensing methods, in part II we present an analysis of the derived glacier inventory data for 662 glaciers and an application to glacier volume and volume-change calculations. Among other things, the analysis reveals a mean glacier elevation of 990 m, with a weak dependence on aspect and a close agreement of the arithmetic mean with the statistical mean elevation as derived from the DEM. A strong scatter of mean slope is observed for glaciers <1 km2, and the derived glacier thickness differs by a factor of two for glaciers of the same size. For the period from about 1920 to 2000 the relative area change is –12.5% (264 glaciers), with a strong dependence on glacier size. Mean mass loss as derived from volume changes is about –0.15 mw.e. a–1.

Accumulation rates and their spatio-temporal variability are important boundary conditions for ice-flow models. The depths of radar-detected internal layers can be used to infer the spatial variability of accumulation rates. Here we infer accumulation rates from three radar layers (26, 35 and 41 ka old) in the Vostok Subglacial Lake region using two methods: (1) the local-layer approximation (LLA) and (2) a combination of steady-state flowband modeling and formal inverse methods. The LLA assumes that the strain-rate history of a particle traveling through the ice sheet can be approximated by the vertical strain-rate profile at the current position of the particle, which we further assume is uniform. The flowband model, however, can account for upstream strain-rate gradients. We use the LLA to map accumulation rates over a 150 km × 350 km area, and we apply the flowband model along four flowbands. The LLA accumulation-rate map shows higher values in the northwestern corner of our study area and lower values near the downstream shoreline of the lake. These features are also present but less distinct in the flowband accumulation-rate profiles. The LLA-inferred accumulation-rate patterns over the three time periods are all similar, suggesting that the regional pattern did not change significantly between the start of the Holocene and the last ~20 ka of the last Glacial Period. However, the accumulation-rate profiles inferred from the flowband model suggest changes during that period of up to 1 cma–1 or ~50% of the inferred values.

Due to the lack of measurements of ice velocity, mass balance, glacier geometry and other baseline data, model-based studies of glacial systems in the Nepal Himalaya are very limited. Here a numerical ice-flow model has been developed for glacier AX010 in order to study its relation to local climate and investigate the possible causes of its general retreat since the end of the Little Ice Age. First, an attempt is made to simulate the historical front variations, considering each climatic parameter separately. Good agreement between the observations and model projections can be obtained under the assumption that variations in glacier front position are a response to changes in temperature alone. The same assumption is made about future changes to explore the 21st-century evolution of the glacier. Under a no-change scenario, the glacier will retreat by another ∽600m by AD 2100, whereas it is projected to vanish completely during this century for all trends with a temperature rise larger than +2.5˚C by AD 2100 with respect to the 1980–99 mean. With constant precipitation at the 1980–99 mean, the model predicts that the glacier will cease to exist at AD 2083, 2056 or 2049 if the temperature rises linearly by 3˚C, 4.5˚C or 6˚C respectively by the end of this century. With an additional range of precipitation changes between –30% and +30%, the life expectancy of glacier AX010 varies by 18, 6 and 2 years for the respective temperature rises. Thus the role of precipitation becomes minimal for the higher trends of temperature rise.

Conventional and reference-surface mass-balance data from Gulkana and Wolverine Glaciers, Alaska, USA, are used to address the questions of how rapidly these glaciers are adjusting (or ‘responding’) to climate, whether their responses are stable, and whether the glaciers are likely to survive in today’s climate. Instability means that a glacier will eventually vanish, or at least become greatly reduced in volume, if the climate stabilizes at its present state. A simple non-linear theory of response is presented for the analysis. The response of Gulkana Glacier is characterized by a timescale of several decades, but its stability and therefore its survival in today’s climate are uncertain. Wolverine seems to be responding to climate more slowly, on the timescale of one to several centuries. Its stability is also uncertain, but a slower response time would make it more susceptible to climate changes.