Bakaninbreen, Svalbard, was in the late phase of a prolonged surge during 1995–96. In order to investigate the nature of the basal mechanical and hydraulic regime, sets of tilt sensors equipped with pore-pressure sensors were installed into the basal till layer. Water-pressure transducers were also installed in the boreholes above the ice–till interface. During 1995, the measured mean strain rate in the till was 0.44 a−1 up-glacier of the surge front, while down-glacier the rate was 0.04 a−1. Our interpretation indicates that deformation was the dominant component of basal motion. Strain rates fluctuated strongly, and periods of both positive (down-glacier) and negative (up-glacier) strain occurred. No consistent relationship was apparent between effective pressure and strain rate within the basal till layer. The viscosity of the basal sediments ranged from 1.3 × 1012 to 3.4 × 1012 Pa s. Such high values suggest that the till layer beneath Bakaninbreen may be partially frozen. Analysis of simultaneous measurements of pore-water and borehole water pressure suggests the hydraulic diffusivity of the till layer ranges from 1.2 to 5.0 × 10−6 m2 s−1, and the hydraulic conductivity from 3.2 × 10−7 to 7.7 × 10−8 m s−1.

This research investigates snow stability on the eastern side of a small mountain range in southwest Montana, U.S.A., on one mid-season day and one late-season day during the 1996/97 winter. Although previous research has addressed snow stability at smaller spatial scales, this is the first field-based study to investigate snow stability (as measured by stability tests) over a mountain range in order to better understand its spatial distribution and the implications for predicting dry-slab avalanches. Using helicopter access, six two-person sampling teams collected data from over 70 sites on each of the two sampling days. Variables for terrain, snowpack, snow strength and snow stability were generated from the field data, and analyzed using descriptive statistics, correlation analysis and multiple regression. Results from the first sampling day show stability is only weakly linked to terrain, snowpack and snow-strength variables due to consistently stormy weather conditions leading up to that day. The second field day’s results demonstrate a stronger relationship between stability and the other variables due to more variable weather conditions that ranged from periods of sunshine to storms. On both days stability decreased on high-elevation, northerly-facing slopes. The data-structure complexity provides insights into the difficulties faced by both scientists and conventional avalanche forecasters in predicting snow avalanches.

Glaciar Chacaltaya is an easily accessible glacier located close to La Paz, Bolivia. Since 1991, information has been collected about the evolution of this glacier since the Little Ice Age, with a focus on the last six decades. The data considered in this study are monthly mass-balance measurements, yearly mappings of the surface topography and a map of the glacier bed given by ground-penetrating radar survey. A drastic shrinkage of ice has been observed since the early 1980s, with a mean deficit about 1 m a−1 w.e. From 1992 to 1998, the glacier lost 40% of its average thickness and two-thirds of its total volume, and the surface area was reduced by >40%. With a mean estimated equilibrium-line altitude lying above its upper reach, the glacier has been continuously exposed to a dominant ablation on the whole surface area. If the recent climatic conditions continue, a complete extinction of this glacier in the next 15 years can be expected. Glaciar Chacaltaya is representative of the glaciers of the Bolivian eastern cordilleras, 80% of which are small glaciers (<0.5 km2). A probable extinction of these glaciers in the near future could seriously affect the hydrological regime and the water resources of the high-elevation basins.

In the low latitudes there is an absence of major thermal seasonality, yet there are three different climate regimes related to global circulation patterns and their seasonal oscillation: the humid inner tropics, the dry subtropics and, intermediate between these two, the outer tropics. For the respective glacier regimes the vertical profiles of specific mass balance (VBPs) are modeled considering vertical gradients of accumulation, air temperature and albedo, the duration of the ablation period and a factor for the ratio between melting and sublimation. The model is first calibrated with data from Hintereisferner, Austrian Alps, and is then applied to tropical conditions. The simulated VBP matches well the measured profiles from Irian jaya and Mount Kenya. Due to lack of field evidence, the subtropical VBP cannot be verified directly. However, application of the respective model versions separately to the humid and dry seasons of the outer-tropical Glaciar Uruashraju, Cordillera Blanca, Peru, provides reasonable results. Glaciers in the humid inner tropics are considered to be most sensitive to variations in air temperature, while dry subtropical glaciers are most sensitive to changes in air humidity. The two seasons of the outer tropics have to be viewed from these different perspectives.

Causes of ice-sheet layering at ice depths greater than about 900 mina transect between Dome C and the Transantarctic Mountains are examined using 60 MHz radar data, collected in the 1970s by the U.K.–U.S.–Danish collaboration. Normally, a dual-frequency technique is required for accurate determination of internal reflection mechanisms. However, by extracting the depth-related features of 60 MHz radar profiles and comparing them with the dual-frequency data collected by the Japanese Antarctic Research Expedition, we have identified a simple method to estimate internal reflection mechanisms. Two zones can be distinguished: (1) the CA zone, where change in electrical conductivity due to variation in acidity is the major cause of internal reflection, and (2) the PCOF zone, where change in dielectric permittivity due to crystal-orientation fabrics is the major cause of internal reflections. Our analysis shows that the radar data reveal the development of PCOF zones in regions where large amounts of ice shearing are expected. This analysis shows how a similar interpretation of the full radar-data archive may reveal information on internal reflection mechanisms across a large part of the East Antarctic ice sheet.

Thwaites Glacier, the second largest ice stream in West Antarctica, drains an area of 166 500 ± 2000 km2 which accumulates 55 ± 5 Gt a−1 (or 60 ± 6 km3 ice a−1) into the Amundsen Sea, unrestrained by an ice shelf. Using interferometric synthetic-aperture radar (InSAR) data collected by the European Remote-sensing Satellites (ERS-1 and -2) in 1996, an output flux of 71 ±7 Gt a−1 (or 77 ± 8 km3 ice a−1) is estimated at the grounding line, where ice thickness is deduced from hydrostatic equilibrium. A similar flux, 70 ± 7 Gt a−1 (or 76 ± 8 km3 ice a−1), is obtained at a gate located 20 km upstream, where ice thickness was measured in 1978 by ice-sounding radar. Total accumulation in between the two gates is 1.6 Gt a−1, or 1.8 km3 ice a−1. Ice discharge therefore exceeds mass accumulation by 30 ± 15%, and Thwaites Glacier must be thinning and retreating at present. The InSAR data show that the glacier floating ice tongue exerts no back pressure on the inland ice, calves into tabular icebergs along a significant fraction of its grounding line, and has a grounding-line thickness which exceeds a prior-calculated limit for stability. Glacier thinning is confirmed at the coast by the detection of a 1.4 ± 0.2 km retreat of its grounding line between 1992 and 1996 with InSAR, which implies 3.2 ± 0.6 m ice a−1 thinning at the glacier center and less near the sides. These results complement the decimeter-scale annual surface lowering observed with satellite radar altimetry several hundred km inland of the grounding line. The magnitude of ice thinning estimated at the coast, however, rules out temporal changes in accumulation as the explanation for surface lowering. Ice thinning must be due to changes in ice flow.

Debris-laden basal ice is exposed along an ice cliff near Hamna Glacier, Sôya Coast, East Antarctica. The basal ice is about 6.8 m thick and shows conspicuous stratigraphic features. The upper 5.5 m consists of alternating layers of bubble-free and bubbly ice. δ values of the bubble-free ice layers are enriched by 2.4 ±1.0‰ (standard deviation) for δ18O compared to values of neighboring bubbly ice layers above, and slopes of δ18O vs δD are close to 8. Such layers are suggested to have been formed by refreezing of meltwater in an open system. In contrast, part of the bubbly ice layers shows neutral profiles for stable isotopes, suggesting that these ice masses are undisturbed ice-sheet ice which was not affected by melting and freezing. The massive alternating layers are thus considered to have been formed by folding of refrozen and non-melted layers. The lower 1.3 m consists predominantly of bubble-free massive ice. The profile of co-isotopic values shows a change of about 3.0‰ for δ18O at the interface between bubble-free and bubbly ice. Since the isotopic change occurred over a wider thickness than the upper 5.5 m, the basal ice is suggested to have been formed by refreezing of meltwater on a larger scale than the upper 5.5 m.

Spatial patterns in residual bed reflection power (BRPr), derived from ground-based radio-echo sounding, were mapped and interpreted in terms of the thermal and hydrological conditions at the base of a high-Arctic polythermal glacier (John Evans Glacier, Ellesmere Island, Canada). BRPr is the residual from a statistical relationship between measured bed reflection power and ice thickness that describes the rate of dielectric loss with depth in the glacier. We identified three types of thermal structure: (a) Positive BRPr and an internal reflecting horizon occur over the glacier terminus. The reflecting horizon is interpreted as the boundary between warm and cold ice, and suggests the presence of a warm basal layer. (b) Positive BRPr occurs without an internal reflector in the upper part of the ablation zone. This suggests that ice is at the pressure-melting point only at the bed. (c) Negative BRPr without an internal reflector occurs in all other regions, suggesting cold ice at the bed. Where BRPr is positive, its pattern is similar to the pattern of subglacial water flow predicted from the form of the subglacial hydraulic equipotential surface. This suggests that hydrological conditions at the glacier bed are a major control on BRPr, probably because the dielectric contrast between ice and water is higher than that between ice and other subglacial materials.

A new model, the variable width/depth ratio (VWDR) model, is used to analyze longitudinal variations in cross-section morphology along glacial valleys. In the VWDR, the cross-sectional shape of a valley is expressed as a function of the width/depth ratio at various heights above the valley floor. Two parameters, m, a measure of the breadth of the valley floor, and n, a measure of the steepness of the valley sides, appear in the model. We have used the VWDR model to study morphological variations of cross-sections along glacial valleys in the middle Tien Shan mountains, China, and find that: (1) in valleys without tributaries, m increases (the valley floor becomes wider) and n becomes more negative (valley sides become steeper) from the head of the valley to the equilibrium line, and then m decreases and n becomes less negative to the end of the valley; (2) in valleys with tributaries, a similar pattern is observed, with an oscillating maximum in m and minimum in n in those sections where a tributary enters the main valley. These characteristics are believed to reflect a maximization of glacial erosion potential in the vicinity of the equilibrium line and in locations of confluence.

We examine the ability of synthetic aperture radar (SAR) to detect the equilibrium line on the glacier Austre Okstindbreen, Norway, using multipolarization SAR images in C- and L-band acquired with the Electromagnetic Institute of the Technical University of Denmark’s airborne EMISAR sensor during the European Multisensor Airborne Campaign EMAC ’95. The late-summer snowline, approximating the equilibrium line, cannot be seen on the SAR images. Instead, photographs from Austre Okstind-breen show that a distinct boundary visible on the C-band SAR images corresponds to the firn line created by old snow from many previous years. This is better seen on the cross-polarized SAR images (HV and VH), which in general reveal more detail than the HH- and VV-polarized images. We model the stratigraphy from net balance and glacier velocity data to calculate the firn-line altitude (FLA). Modelled FLA and the observed boundary are separated by 50 m in elevation, but considering errors during co-registration and modelling we conclude that the observed boundary on Austre Okstindbreen is the firn line. Monitoring FLA rather than equilibrium-line altitude (ELA) for mass-balance studies with remote sensing is therefore suggested.

We use a self-gravitating viscoelastic model of the Earth and a dynamically consistent marine ice-sheet model to study the relationships between marine ice-sheet dynamics, relative sea level, basal topography and bedrock dynamics. Our main conclusion is that sea-level change and lithospheric coupling are likely to have played limited roles in the postglacial retreat of marine ice sheets. The postglacial rise in sea level would only have caused at the most around 100 km of grounding-line retreat for an ice sheet of similar dimensions to the West Antarctic ice sheet, compared with the several hundred km of retreat which has occurred in the Ross Sea. There is no evidence that reverse slopes lead to instability. Incorporating coupling with lithospheric dynamics does not produce markedly different effects. The implication of these studies is that marine ice-sheet retreat is the result of physical mechanisms other than lithospheric coupling and sea-level rise.

Marine ice sheets with mechanics described by the shallow-ice approximation by definition do not couple mechanically with the shelf. Such ice sheets are known to have neutral equilibria. We consider the implications of this for their dynamics and in particular for mechanisms which promote marine ice-sheet retreat. The removal of ice-shelf buttressing leading to enhanced flow in grounded ice is discounted as a significant influence on mechanical grounds. Sea-level rise leading to reduced effective pressures under ice streams is shown to be a feasible mechanism for producing postglacial West Antarctic ice-sheet retreat but is inconsistent with borehole evidence. Warming thins the ice sheet by reducing the average viscosity but does not lead to grounding-line retreat. Internal oscillations either specified or generated via a MacAyeal–Payne thermal mechanism promote migration. This is a noise-induced drift phenomenon stemming from the neutral equilibrium property of marine ice sheets. This migration occurs at quite slow rates, but these are sufficiently large to have possibly played a role in the dynamics of the West Antarctic ice sheet after the glacial maximum. Numerical experiments suggest that it is generally true that while significant changes in thickness can be caused by spatially uniform changes, spatial variability coupled with dynamical variability is needed to cause margin movement.

The region where two active tributaries feed into the now stagnant Ice Stream C (ISC), West Antarctica, is thickening. In this region, we observe a correlation between faster ice flow (the tributaries) and elevated topography. We conclude that stagnation of ISC resulted in compression and thickening along the tributaries, eventually forming a “bulge” on the ice-sheet surface. Modern hydraulic potential gradients would divert basal meltwater from ISC to Ice Stream B (ISB). These gradients are primarily controlled by the bulge topography, and so likely formed subsequent to trunk stagnation. As such, we argue against “water piracy” as being the cause for ISC’s stagnation. Kinematic-wave theory suggests that thickness perturbations propagate downstream over time, but that kinematic-wave speed decreases near the stagnant trunk. This and modest diffusion rates combine to trap most of the tributary-fed ice in the bulge region. Using interfero-metric synthetic aperture radar velocity measurements, we observe that half of the ice within ISC’s southern tributary flows into ISB. That flow pattern and other observations of non-steady flow in the region likely result from stagnation-induced thickening along upper ISC combined with a longer period of thinning on upper ISB. If current trends in thickness change continue, more ice from upper ISC will be diverted to ISB.

The sediment in icebergs offers an opportunity to sample the otherwise inaccessible beds of modern tidewater glaciers. This paper presents a discussion of the benefits and problems of using iceberg debris to make inferences about the nature of the beds of modern tidewater glaciers. As an example, we present data obtained from icebergs calved from fast-flowing tidewater glaciers in Kongsfjorden, Svalbard. Sediment samples obtained from 12 icebergs show that clast-rich muddy diamicton is the dominant facies, although muddy gravel and mud are also present in some debris layers. Sediment texture and clast-shape analyses indicate that most debris in the icebergs is of basal origin. Micro-morphological analyses of thin sections of the diamicton also indicate pervasive deformation of the sediment. Although individual icebergs only represent small samples of the bed of a tidewater glacier, and problems remain in linking iceberg sediments to their precalving transport history beneath the glacier, this method provides a potentially fruitful way of linking sediments to glacier dynamics.

Measurements of the surface and internal layer geometry from ice-penetrating radar and global positioning system surveys on three inter-ice-stream ridges in West Antarctica (Siple Dome, ridge DE and ridge BC) are examined with ice-flow models to infer (1) the history of the divide position at each site and (2) the spatial pattern of accumulation across the ridges. We find that the divide position is most steady at Siple Dome, somewhat steady at ridge DE and highly variable at ridge BC. Data from Siple Dome and ridge DE show evidence for steady northward motion of the ice divide for the past few thousand years. The layers beneath ridge BC suggest a 5 km northward shift of the divide position within the past several hundred years. Assuming the divide shifts are all due to changing elevation of the bounding ice streams, we infer the relative elevation history for segments of Ice Streams B–E. The northward displacement of the divide for all ridges implies a progressive relative thinning of the ice streams from E to B, with most dramatic recent thinning (100 m in <103 years) of Ice Stream B relative to Ice Stream C. Analysis of the internal layer pattern across the ridges indicates a south–north accumulation gradient with higher accumulation rates on the northern flanks of the ridges in all three cases. The inferred accumulation distribution is nearly uniform on the northern flanks, decreases sharply within a few ice thicknesses across the divides, and then decreases gradually farther to the south. The north/south decrease is strongest for ridge DE and weakest for ridge BC. This spatial pattern and the reduction in gradient strength with distance from the Amundsen Sea is consistent with the hypothesis that storms from the Amundsen Sea carry moisture first south then west over West Antarctica and deposit more snow on the windward side of ridges due to orographic lifting. This pattern has been stable for at least the past several thousand years.

Modelling the heat flow in small, cold high-altitude glaciers is important for the interpretation of paleoclimatic data from ice cores. Coupled glacier-flow and heat-flow models are presented that incorporate the densification, heat advection and possible phase transitions at the permafrost boundaries within the bedrock. Marked bends observed in the temperature profiles from two recent boreholes on Colle Gnifetti, Swiss Alps, are interpreted with the help of a transient heat-flow model, driven with a temperature history. The conclusion is that substantial warming of the mean firn temperature at shallow depths has taken place over the last few decades. This has not been observed before in cold-firn regions of the Alps. Modelled heat fluxes in the Monte Rosa massif are strongly influenced by the mountain topography. This leads to a spatial variability of the temperature gradient near the glacier base which has been observed in boreholes to the bedrock. In order to match the measured temperature profiles in the glacier, the vertical heat flux at great depth must be set to an extremely low value. It is shown with the help of the transient heat-flow model that this is a paleoclimatic effect, possibly enhanced by a degrading permafrost base.

We have mapped Antarctic blue-ice areas using the U.S. National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR) Antarctica cloud-free image mosaic established by the United States Geological Survey. The mosaic consists of 38 scenes acquired from 1980 to 1994. Our results show that approximately 60 000 km2 of blue ice exist for each of the two main types of blue ice: “melt-induced” and “wind-induced”. Normally, the former type is located on slopes in coastal areas where climate conditions (i.e. persistent winds and temperature), together with favourable surface orientation, sustain conditions for surface and near surface melt. The latter blue-ice category occurs near mountains or on outlet glaciers, often at higher elevations, where persistent winds erode snow away year-round, and combined with sublimation creates areas of net ablation. Furthermore, we have identified an additional area of 121 000 km2 as having potential for blue ice. However, in these areas features such as mixed pixels, glazed snow surfaces, crevasses and/or shadows make interpretation more uncertain. In conclusion, a conservative estimate of Antarctic blue-ice area coverage by this method is found to be 120 000 km2 (∼0.8% of the Antarctic continent), with a potential maximum of 241 000 km2 (∼1.6% of the Antarctic continent).

In 1997 a 121 m ice core was retrieved from Lomonosovfonna, the highest ice field in Spitsbergen, Svalbard (1250 m a.s.l.). Radar measurements indicate an ice depth of 126.5 m, and borehole temperature measurements show that the ice is below the melting point. High-resolution sampling of major ions, oxygen isotopes and deuterium has been performed on the core, and the results from the uppermost 36 m suggest that quasi-annual signals are preserved. The 1963 radioactive layer is situated at 18.5–18.95 m, giving a mean annual accumulation of 0.36 m w.e. for the period 1963–96. The upper 36 m of the ice core was dated back to 1920 by counting layers provided by the seasonal variations of the ions in addition to using a constant accumulation rate, with thinning by pure shear according to Nye (1963). The stratigraphy does not seem to have been obliterated by meltwater percolation, in contrast to most previous core sites on Svalbard. The anthropogenic influence on the Svalbard environment is illustrated by increased levels of sulphate, nitrate and acidity. Both nitrate and sulphate levels started to increase in the late 1940s, remained high until the late 1980s and have decreased during the last 15 years. The records of δ18O, MSA (methane-sulphonic acid), and melt features along the core agree with the temperature record from Longyearbyen and the sea-ice record from the Barents Sea at a multi-year resolution, suggesting that this ice core reflects local climatic conditions.