We have obtained common offset, common midpoint (CMP) and borehole vertical (VRP) ground-penetrating radar profiles close to the margin of Falljökull, a small, steep temperate valley glacier situated in southeast Iceland. Velocity analysis of CMP and VRP surveys provided a four-layered velocity model. This model was verified by comparison between the depths of englacial reflectors and water channels seen in borehole video, and from the depths of boreholes drilled to the bed. In the absence of sediment within the glacier ice, radar velocity is inversely proportional to water content. Using mixture models developed by Paren and Looyenga, the variation of water content with depth was determined from the radar velocity profile. At the glacier surface the calculated water content is 0.23–0.34% (velocity 0.166 m ns−1), which rises sharply to 3.0–4.1% (velocity 0.149 m ns−1) at 28 m depth, interpreted to be the level of the piezometric surface. Below the piezometric surface the water content drops slowly to 2.4–3.3% (velocity 0.152 m ns−1) until ∼102 m depth where it falls to 0.09–0.14% (velocity 0.167 m ns−1). The water content of the ice then remains low to the glacier bed at about 112 m. These results suggest storage of a substantial volume of water within the glacier ice, which has significant implications for glacier hydrology, ice rheology and interpretations of both radar and seismic surveys.

New maps are presented of three internal ice-sheet radio-echo sounding (RES) layers in the region 73.5–75.75° S, 120–127° E (56 000 km2) around Dome C, central East Antarctica. These layers represent horizons of enhanced acidity resulting from volcanic aerosol deposition, identified from analogue RES data. They are continuous over the entire mapped area, and constitute deformation markers in the ice column. Internal RES layers were initially identified from discrete radar power reflection coefficient profiles and subsequently digitized directly from prints of ice-sheet cross-sections, acquired by continuous RES profiling. Georeferenced vector data are used to generate a 5 km gridcell raster of depth for each internal RES layer, as a basis for contour mapping. Ice deformation in the Dome C region is significant because this is the location of the European Project for Ice Coring in Antarctica. Since internal layers are isochronous, the one-dimensional ice-core data at Dome C can be correlated over the survey area to produce a three-dimensional context.

Precipitation at 2050 m on Blue Glacier, U.S.A., was measured daily from August 1957 through July 1958. Its correlation with a nearby lowland station with a good, long-term (1914–96) record is used to estimate precipitation on the glacier over that entire period. Average annual precipitation on Blue Glacier is 4500 mm w.e. Snowfall depends on the joint distribution of precipitation and temperature. Over the period 1948–96, for which twice-daily radiosonde observations are available, temperature at any elevation on the glacier is interpolated in the radiosonde profile to partition the precipitation as either rain or snow. Daily partitioning is preferred, especially during spring and autumn storms when averaging over longer periods may substantially under- or overestimate snowfall on the glacier. Prior to 1948, snowfall is estimated from the mean over 1948–96, in a particular month and elevation, of the fraction of the precipitation falling as snow. The standard error in the October–May snowfall at 2100 m is estimated to be 250 mm w.e. during the radiosonde era (1948–96) and 350 mm prior to that. For the first 10 years or so after mass-balance measurements began at Blue Glacier (1957), precipitation increased and winter temperature at 850 mbar (about 1450 m) decreased, but since then the trends have reversed. The combined effect, increasing snowfall until 1965 and decreasing since, closely parallels measured mass changes of Blue Glacier. When the average vertical profile of total annual snowfall is subjected to a hypothetical 1 K warming, the resulting reduction in snowfall is greatest at the glacier terminus and decreases up-glacier; the average over the entire glacier is 300 mm w.e.

Ice motion over Lake Vostok, Antarctica, is measured using repeat-pass synthetic-aperture radar (SAR) interferometry. The coverage of the lake and the components of the vector field are resolved using 10 overlapping data takes from ascending and descending look directions. Seventy-day temporal baselines provide the sensitivity required to observe the range of ice motion (0–6 m a−1) over the lake and the adjacent ice sheet. It is remarkable that the scattering field remained coherent over these time separations. This is critical for interferometric analysis and can be attributed to the low surface accumulation and low air temperature at this elevation. The regional flow of the ice sheet around Lake Vostok is from west to east, perpendicular to the surface elevation contours. As the ice flows past the grounding line, a southward component of motion develops that is correlated with the north–south surface slope along the length of the lake. The surface velocity increases slowly from the northern tip of the lake and then more rapidly south of 77° S. At Vostok station, the ice motion is 4.2 m a−1. Across the lake and away from boundary effects, the down-lake flow pattern takes on a parabolic profile with maximum velocity close to the center line of the lake. The overall influence of the subglacial lake is the addition of a down-lake motion component to the prevailing west–east motion of the ice sheet. As a result, we estimate 10% of the mass flowing onto the lake is diverted south. Reconstructions based on the Vostok ice core indicate that the ice was grounded up-glacier from the core site approximately 5000 years ago. This suggests a minimum freezing rate of 40 mm a−1 for the subglacial accretion ice, 10 times greater than that inferred from thermodynamic modeling of the upper 2 km of the ice core.

Permafrost existed around and under marginal parts of the southern Laurentide ice sheet during the Last Glacial Maximum. The presence of permafrost was important in determining the extent, form and dynamics of ice lobes and the landforms they produced because of influences on resistance to basal motion and subglacial hydrology. We develop a two-dimensional time-dependent model of permafrost and glacier-ice dynamics along a flowline to examine: (i) the extent to which permafrost survives under an advancing ice lobe and how it influences landform development and hydrology, and (ii) the influence of permafrost on ice motion and surface profile. The model is applied to the Green Bay lobe, which terminated near Madison, Wisconsin, during the Last Glacial Maximum. Simulations of ice advance over permafrost indicate that the bed upstream of the ice-sheet margin was frozen for 60–200 km at the glacial maximum. Permafrost remained for centuries to a few thousand years under advancing ice, and penetrated sufficiently deep (tens of meters) into the underlying aquifer that drainage of basal meltwater became inefficient, likely resulting in water storage beneath the glacier. Our results highlight the influence of permafrost on subglacial conditions, even though uncertainties in boundary conditions such as climate exist.

This paper documents the glaciological structures associated with the surge of Skeiðarárjökull, Iceland, in 1991. These structures are interpreted as units of stratified ice, low-angle fractures, vertical and sub-vertical fractures (crevasse traces) and thrusts. The inferred thrusts are debris-rich and, unusually, have both down-glacier and up-glacier dips close to the ice margin. Sediment infills consist of either massive sand or horizontally stratified sand units. The most significant debris-rich structures on the glacier surface, however, are supraglacial crevasse and conduit fills, which contain either massive or horizontally stratified silts, sands and granule-gravels. These sediments infill both vertical fractures (relict crevasses) and englacial conduits. At the stratigraphic base of these sediment fills there is evidence of syn-sedimentary deformation, suggesting that sedimentation occurred during crevasse closure and continued thereafter. We argue that these structures relate to an episode of supraglacial meltwater flow during the 1991 surge, caused by the build-up of subglacial water pressure in a linked-cavity system or some similar distributed drainage system beneath the glacier. The development of this high-level drainage route may have helped regulate basal water pressures and therefore the active phase of the surge. The idea that the supraglacial leakage of subglacial water may have played a role in terminating the surge is explored.

We combine European Remote-sensing Satellite (ERS-1) radar altimeter surface elevations (Fricker and others, 2000) with six different accumulation distributions to compute balance fluxes for the Lambert Glacier–Amery Ice Shelf drainage system. These interpolated balance fluxes are compared with fluxes derived from in situ measurements of ice thickness and velocity at 73 stations of the Lambert Glacier basin traverse and at 11 stations further downstream, to assess the system’s state of balance. For the upstream line we obtain a range of imbalance estimates, from −23.8% to +19.9% of the observed flux, reflecting the sensitivity to the accumulation distributions. For some of the accumulation distributions the imbalance estimates vary significantly between different parts of the line. Imbalance estimates for the downstream line range from −17.7% to +70.2%, with four of the estimates exceeding +30%, again reflecting the sensitivity of the result to input accumulation, and strongly suggesting that the mass balance of the region between the two lines is positive. Our results confirm the importance of accurate estimates of accumulation in ice-sheet mass-balance studies. Furthermore, they suggest that it is not possible to accurately determine the state of balance of large Antarctic drainage basins on the basis of currently available accumulation distributions.

A review of in situ and remote-sensing data covering the ice shelves of the Antarctic Peninsula provides a series of characteristics closely associated with rapid shelf retreat: deeply embayed ice fronts; calving of myriad small elongate bergs in punctuated events; increasing flow speed; and the presence of melt ponds on the ice-shelf surface in the vicinity of the break-ups. As climate has warmed in the Antarctic Peninsula region, melt-season duration and the extent of ponding have increased. Most break-up events have occurred during longer melt seasons, suggesting that meltwater itself, not just warming, is responsible. Regions that show melting without pond formation are relatively unchanged. Melt ponds thus appear to be a robust harbinger of ice-shelf retreat. We use these observations to guide a model of ice-shelf flow and the effects of meltwater. Crevasses present in a region of surface ponding will likely fill to the brim with water. We hypothesize (building on Weertman (1973), Hughes (1983) and Van der Veen (1998)) that crevasse propagation by meltwater is the main mechanism by which ice shelves weaken and retreat. A thermodynamic finite-element model is used to evaluate ice flow and the strain field, and simple extensions of this model are used to investigate crack propagation by meltwater. The model results support the hypothesis.

Towards accounting for the dynamic response of glaciers and ice caps in the estimation of their contribution to sea-level rise due to global warming, a mass-balance degree-day model is coupled to a geometric glacier model. The ice dynamics are treated implicitly in the geometric model by using scaling parameters that have been extensively investigated in the literature. The model is tested by presenting a case-study of the glacier Hintereisferner, Austrian Alps. The results are compatible with geomorphological data and other modelling studies. An estimate is made of the volume decrease due to initial disequilibrium. An extensive sensitivity study using generalized glacier shapes and sizes allows a comparison of results with dynamic theory. According to the geometric model, glaciers with a narrowing channel change more with a change in mass balance than glaciers with a widening channel, due to their shape and the way in which that shape changes with a changing climate. Also their response time is longer. As time progresses after a mass-balance perturbation, the longer response time for continental glaciers compared to glaciers with larger mass turnover offsets the effect of their smaller static sensitivity. Thus, although for the next century we may expect greater changes in volume from alpine glaciers, the equilibrium or committed change is greater for the continental glaciers.

We have defined the complex bed topography for a section of a small temperate glacier using 50 MHz monostatic short-pulse radar data and a synthetic-aperture array-processing method. The data were collected on a 100 m by 340 m array grid in the upper stem of Gulkana Glacier, central Alaska, U.S.A. The array processing was based on a modified three-dimensional (3-D) Kirchhoff migration integral and implemented with a synthetic-aperture approach that uses sequences of overlapping sub-arrays to generate depth images in vertical planes. Typical sub-array beam patterns are generally <5° at the −6 dB level, giving a flashlight-like searching capability without distorting the wavelet shape. The bed topography was constructed using normal reflections picked from 3-D array depth images. In some instances reflections were imaged outside the data-cover-age area. The bed surface dips steeply, both parallel and transverse to the direction of ice flow. The maximum observed depth is roughly 140 m. The 3-D method resolved bed dips up to 45°. In regions of steepest dip, it improved depth accuracy by 36% compared with raw data, and by 15% compared with standard two-dimensional (2-D) migration. Over 12 dB of signal-to-noise improvement and improved spatial resolution was achieved compared to raw data and 2-D migration. False bottom layering seen in the raw data and in 2-D migrations is not observed in the 3-D array results. Furthermore, loss of bottom reflections is shown by the 3-D migration to be attributable to the dip and curvature of the reflector, and not scattering losses or signal clutter from englacial inclusions.

Analysis of air extracted from bubbles of polar ice showed the anthropogenic increase of atmospheric CO2 concentration during the past few hundred years (preindustrial concentration 280 ppmv) and, unexpectedly, that the concentration also increased due to natural causes from 200 to 280 ppmv during the transition from the last glacial epoch to the Holocene. However, more detailed reconstructions based on ice-core analyses suggested that some of the additional observed variations were actually due to the modification of CO2 concentration in the bubbles by chemical reactions between impurities in the ice. Detailed analyses of acidity, the carbonate concentration, the concentration of oxidation agents like hydrogen peroxide (H2O2), and organic compounds (e.g. formaldehyde (HCHO)) along short intervals of ice cores, representing only a few annual layers, allowed us to investigate these chemical reactions. The records are not conclusive but they indicate that the oxidation of organic compounds is at least as important as acid–carbonate reactions. Ice containing a low mean carbonate and a low H2O2 concentration, as well as a small scattering of CO2 results obtained on several adjacent samples representing a few annual layers at most, is at present the best guarantee for reliable results of the atmospheric CO2 concentration.

Balance velocities for the Greenland ice sheet have been calculated from a new digital elevation model (DEM), accumulation rates and an existing ice-thickness grid, using a fully two-dimensional finite-difference scheme. The pattern of velocities is compared with velocities derived from synthetic-aperture radar (SAR) interferometry for three different regions of the ice sheet. Differences between the two estimates of velocity highlight the respective strengths and weaknesses of the datasets and techniques used. A comparison with ten global positioning system-derived velocities indicates that the balance-velocity scheme and input datasets used here provide a remarkably good representation of the velocity distribution inland from the margins. These balance-velocity data, therefore, could help constrain numerical ice-sheet models. The balance velocities were found to be unreliable close to the ice-sheet margins due to larger errors in ice thickness, surface slope and ablation rate in this region. Comparison of the balance velocities with SAR interferometry in the region of the “Northeast Greenland Ice Stream” indicates the importance of the smoothing distance that must be applied to the DEM before calculating balance velocities. A smoothing distance of 20 times the ice thickness gave good agreement between the two measures of velocity.

Snow algae in shallow ice cores (7 m long) from Yala Glacier in the Langtang region of Nepal were examined for potential use in ice-core dating. Ice-core samples taken at 5350 m a.s.l. in 1994 contained more than seven species of snow algae. In a vertical profile of the algal biomass, 11 distinct algal layers were observed. Seasonal observation in 1996 at the coring site indicated most algal growth occurred from late spring to late summer. Pit observation in 1991, 1992 and 1994 indicated that algal layer formation takes place annually. δ18O, chemical ions (Na+, Cl−, SO42− and NO3−) and microparticles failed to show any clear seasonal variation, particularly at depths exceeding 2 m, possibly due to heavy meltwater percolation. Snow algae in ice cores would thus appear to be accurate boundary markers of annual layers and should prove useful for ice-core dating in Himalayan-type glaciers.

A detailed study of a proglacial bedrock site and a subglacial cavity of an outlet of Øksfjordjøkelen, Norway, is presented together with observations from the foreland of Konowbreen, Spitsbergen. Striation directions and subglacial observations indicate that local ice-flow paths were highly variable, deviating at angles of approximately 90° from the main ice-flow direction. Stepped bedrock topography appears conducive to the production of highly variable ice-flow paths, because the high bed roughness creates a locally variable stress regime within the ice, including low-pressure, lee-side areas into which ice can flow. If ice flow is sustained along a specific path and the ice contains debris, then abrasion should produce an erosional bedform. Models are proposed whereby locally variable ice-flow patterns could produce erosional bedforms, which would be described as p-forms, purely through mechanical abrasion.

We present results from satellite imagery, ice-motion surveys and ice-penetrating radar studies of part of the north margin of Ice Stream C, one of the ice streams draining the West Antarctic ice sheet to the “Ross Embayment”. Our studies suggest that the shutdown of Ice Stream C about 150 years ago was not a single event, but a sequence involving stagnation of ice and migration of the ice-stream boundary. Ground-based studies confirm the inference from imagery that a series of former shear zones exist, decreasing in age towards the ice-stream center. A region of ice-stream trunk, including a former margin, lies sheared and folded between the (recent) inner and (older) outer margins of the area. ice-motion and topographic surveys give some constraint on the time of shutdown of the outer margin. The results provide a forum for discussing shutdown mechanisms. Possible causes for the stepwise migration of the north margin of Ice Stream C include a gradual decrease in ice flux, a reduction in the available water or hydrostatic pressure in the basal till, or a freezing of the till layers on the northern side.

Synthetic-aperture radar interferometry data and airborne ice-sounding radar (ISR) data are employed to obtain modern estimates of the inland ice production from Nioghalvfjerdsbræ (NB) and Zachariae Isstrøm (ZI), the two largest glaciers draining the northeast sector of the Greenland ice sheet. Ice fluxes are measured at the grounding line (14.2 ±1 km3 ice a−1 for NB and 10.8 ±1 km3 ice a−1 for ZI) with an ice thickness deduced from ice-shelf hydrostatic equilibrium, and along an ISR profile collected upstream of the grounding line (14.3 ± 0.7 km3 ice a−1 for NB and 11.6 ± 0.6 km3 ice a−1 for ZI). Balance fluxes calculated from a map of snow accumulation and model predictions of surface melt are 11.9 ± 2 km3 ice a−1 for NB and 10.0 ± 2 km3 ice a−1 for ZI at the grounding line, and 12.2 and 10.3 km3 ice a−1, respectively, at the ISR line. The two glaciers therefore exhibit a negative mass balance equivalent to 14% of their balance flux, with a ±12% uncertainty. Independently, we detect a retreat of the grounding line of NB between 1992 and 1996 which is larger at the glacier center (920 ± 250 m) than on the sides (240 ± 50 m). The corresponding ice-thinning rates (2 ± 1 m a−1 at the glacier center and 0.6 ± 0.3 m a−1 on the sides) are too large to be accommodated by temporal changes in ablation or accumulation, and must be due to dynamic thinning.

The ice-sheet margin at Russell Glacier, West Greenland, advanced ∼7 m a−1 between 1968 and 1999. As the ice advanced over moraine ridges, small changes in position caused major changes in the routing of proglacial water and sediment. These included changes in the distribution of ice-marginal lakes, in the periodic drainage of ice-dammed lakes, in the routing and sediment content of meltwater draining into the proglacial zone, and in the release of sediment from the moraines by erosion and mass movements. Proglacial hydrology and sediment flux appear to be controlled not simply by glacier mass balance, but by evolving ice-marginal geomorphology, which must be accounted for in palaeoenvironmental interpretation of proglacial sediments.

A climate–glacier model was used to reconstruct Late-glacial climate conditions from two case-study glaciers at 18° and 22° S in the arid (sub)tropical western Andes of northern Chile. The model uses (i) the geometry of the Late-glacial maximum glaciation, (ii) modern diurnal and annual cycles, amplitudes and lapse rates of the climate, (iii) empirical–statistical sublimation, melt and accumulation models developed for this area, and (iv) dynamic ice flow through two known cross-sections for steady-state conditions. The model is validated with modern conditions and compares favorably with the glaciological features of today. The mass-balance model calculates the modern equilibrium-line altitude at 18° S as high as 5850 m (field data 5800 m), whereas no glaciers exist in the fully arid southern area at 22° S despite altitudes above 6000 m and continuous permafrost. For Late-glacial times, the model results suggest a substantial increase in tropical summer precipitation (ΔP = +840 (− 50/+ 10) mm a−1 for the northern test area; +1000 (− 10/+ 30) mm a−1 for the southern test area) and a moderate temperature depression (ΔT = −4.4 (− 0.1/+ 0.2) °C at 18° S; −3.2 (±0.1) °C at 22° S). Extratropical frontal winter precipitation (June–September) was <15% of the total annual precipitation. A scenario with higher winter precipitation from the westerlies circulation belt does not yield a numerical solution which matches the observed geometry of the glaciers. Therefore, we conclude that an equatorward displacement of the westerlies must be discarded as a possible explanation for the late Pleistocene glaciation in the Andes of northern Chile.