A lightweight, portable drilling system for coring up to 500 m depths has been developed and field-tested. The drilling system includes four major components: (1) an electromechanical (EM) dry-hole drill; (2) an ethanol thermal electric drill; (3) a drill set-up with a 500 m cable capacity; and (4) a controller unit. The system may be switched quickly from a dry-hole EM drill to an antifreeze thermal electric drill. This lightweight system makes ice-core drilling more cost-efficient, and creates a minimal environmental impact. The new EM drill, which recovers 100 mm diameter, 1 m long pieces of ice core, is 3.2 m long and weighs 35 kg. This drill and the drilling set-up were recently tested at the Raven (former Dye 2) site, southern Greenland, where a core was recovered to 122 m. The thermal drill is 2.9 m long and weighs 25 kg. It produces 100 mm diameter, 2.1 m long pieces of ice core, and was tested to 315 m in Franz Josef Land, Eurasian Arctic. The drilling set-up with a 250 m cable weighs about 100 kg (or 128 kg for 500 m of cable). After minor adjustments this drill system retrieved cores of better quality than those recovered by other drill systems under similar glaciological conditions. After adjustments to optimize its performance, the drill retrieved 5.25 m of core per hour over the depth range 0–21 m.

Five ice cores have been retrieved from a transect close to the terminus of Glacier de Tsanfleuron, Switzerland. The cores extend from the ice surface to the glacier bed, and are 3.5–44.8 m long. Stratigraphic logging based on bubble size and density reveals the presence of a highly metamorphosed basal ice layer, about 10 m thick, from which all traces of bubble-rich ice have been removed. This bubble-poor ice, which corresponds closely with clear-facies ice observed in cavities beneath numerous temperate-based glaciers, contrasts with the overlying bubble-rich or bubble-foliated englacial ice and the underlying debris-rich and bubble-free dispersed-facies basal ice.

Down-core patterns in major-ion composition, stable-isotope composition and total gas content and composition are generally consistent with formation of clear-facies ice by deformation-related metamorphism of bubbly, englacial ice. In addition, isotopic data suggest that storage of downward-percolating meltwaters occurs close to the upper surface of the clear-facies ice layer, perhaps reflecting a local variation in ice permeability across the transition from englacial to clear-facies ice. Enrichment in crustally derived ionic species is noted in the lowermost decimetres of the debris-free, clear-facies ice that immediately overlies debris-rich dispersed-facies basal ice. This ionic enrichment in debris-free ice is interpreted in terms of active inter-granular meltwater flow within some decimetres of the glacier bed.

Time-dependent ice-sheet modelling of a 176 000 km2 area in western Dronning Maud Land, East Antarctica, provided information on the ice sheet’s response to six climate-change scenarios. Another experiment was done to study changes in ice thickness, flow and basal temperature conditions between the present ice configuration and a simulated maximum palaeo-ice sheet. The input to the model included new datasets of bed and surface topography compiled for this study. The results of the six climate-change experiments, including a 0.5°C per century global-warming scenario, show that the ice sheet has a robust behaviour with respect to the different climate changes. The maximum change in ice volume was <5% of the initial volume in all climate runs. This is for only relatively short-term climate changes without major changes in global sea level, and also a simulated ice sheet without an ice shelf. The modelled long-term response time of the ice sheet, 20 kyr or more, indicates that the ice sheet may still be adjusting to the climate change that ended the Last Glacial Maximum. In the maximum palaeo-ice-sheet simulation, with a 5°C climate cooling and the grounding line located at the continental-shelf margin, ice thickness increased drastically downstream from the Heimefrontfjella mountain range but remained basically unaffected on the upstream polar plateau. Compared to present conditions, complex changes in basal temperatures were observed. The extent of areas with basal melting increased, for example in the deep trough of the Veststraumen ice stream. Areas at intermediate elevations in the landscape also experienced increased basal temperatures, with significant areas reaching the melting point. In contrast, high-altitude areas that today are clearly cold-based, such as around Heimefrontfjella and Vestfjella and the Högisen dome, experienced a 5–10°C decrease in basal temperatures in the palaeo-ice-sheet reconstruction. The results suggest that the alpine landscape within these mountain regions was formed by wet-based local glaciers and ice sheets prior to the late Cenozoic.

Surface flow in a 10 000 km2 expanse of the onset area of Ice Stream D, West Antarctica, was measured by repeat, precise global positioning system surveys over a 1 year interval. The pattern of velocity and strain rate shows the development of Ice Stream D, the major flow into which originates south of Byrd station and follows the course of a deep bed channel. Plotting of the driving stress vs the ratio of velocity and ice thickness identifies the onset of streaming flow (roughly 140 km downstream of Byrd station) as a transition between deformation flow and sliding flow. Along the kinematic center line of the developing ice stream, the ice rheology is linear at stresses below 0.6 bar, and appears temperate at the base well before the onset of streaming is reached. The onset corresponds to a maximum driving stress of 0.8 bar. It occurs downstream of a slight increase in longitudinal strain rate where stronger along-flow lineations are apparent in Landsat imagery, and after the ice has passed the center of an overdeepening in the bed channel. No current deviation from equilibrium is detected in this region, but a set of flow stripes misaligned with present flow indicates significant changes in flow have occurred in the past.

The spatial distribution of accumulation across Siple Dome, West Antarctica, is determined from analysis of the shapes of internal layers detected by radio-echo sounding (RES) measurements. A range of assumed accumulation patterns is used in an ice-flow model to calculate a set of internal layer patterns. Inverse techniques are used to determine which assumed accumulation pattern produces a calculated internal layer pattern that best matches the shape of internal layers from RES measurements. All of the observed internal layer shapes at Siple Dome can be matched using a spatially asymmetric accumulation pattern which has been steady over time. Relative to the divide, the best-fitting accumulation pattern predicts 40% less accumulation 30 km from the divide on the south flank of Siple Dome and 15–40% more accumulation 30 km from the divide on the north flank. The data also allow the possibility for a small time variation of the pattern north of the divide. The mismatch between the calculated and the observed layer shapes is slightly reduced when the accumulation rate north of the divide is higher in the past (> 5kyr BP) than at present. Sensitivity tests show that the predicted change in the spatial accumulation pattern required to cause the slight Siple Dome divide migration (inferred from previous studies) would be detectable in the internal layer pattern if it persisted for > 2 kyr. Our analysis reveals no evidence that such a change has occurred, and the possible change in accumulation distribution allowed by the data is in the opposite sense. Therefore, it is unlikely that the Siple Dome divide migration has been caused by a temporal change in the spatial pattern of accumulation. This conclusion suggests the migration may be caused by elevation changes in Ice Streams C and D at the boundaries of Siple Dome.

An interferogram reflecting the motion of Black Rapids Glacier, Alaska, U.S.A., was obtained from two European Remote-sensing Satellite (ERS-1) syntheticaperture radar (SAR) images, acquired on 22 and 25 January 1992. We investigate whether the interferometric data are quantitatively consistent with terrestrial velocity measurements along three transverse profiles. These terrestrial data are from different years (1987, 1990, 1996) and cover different periods (6–28 April, 23 May–7 July and a whole year) than the SAR interferogram. Terrestrial ice velocity at the date of the SAR imagery is obtained via seasonal and annual corrections that are calculated from other terrestrial velocity measurements available at higher time resolution for selected sites on the glacier. Interferometric and terrestrial velocity are in excellent agreement if a (terrestrially measured) surface-normal velocity component (v⊥) is properly accounted for. This suggests that both the interferometric velocities and the conversions of terrestrial data to the winter period are reliable. The terrestrial velocity measurements show that ice flow in the upper ablation area (14–16 km sites) changes from longitudinal compression in mid-winter (v⊥ = + 0.82 cm d−1) to moderate longitudinal extension during summer (v⊥ ≈ –0.25 cm d−1). In the lower ablation area, the seasonal variations of the longitudinal strain rate are much smaller: +1.0 and about +0.85 cm d−1 for the respective mid-winter and summer values of v⊥ at the 20 km site.

A new method for ice-fabric analysis which is easily automated is described. The technique relies on a set of digital images of an ice thin section viewed through crossed polarizers and rotated to various angles. From the acquired set of images, the c-axis orientation of all the grains in the image can be determined. The theoretical basis for the technique is described in detail, including corrections for refraction, as well as other corrections for off-axis grains and polarization/reflection effects. An experiment was performed comparing results of the new technique with those of the standard technique. Excellent agreement was obtained. It is expected that an automated system using this technique (currently almost complete) will have a large impact on the amount of information and physical knowledge which can be extracted from ice fabrics in the future.

A degree-day model is used to assess the sensitivity of the mass balance of five Swiss glaciers to temperature changes. The model uses temperature data extrapolated from nearby climate stations, and is tuned by varying precipitation to make the model fit the observed distribution of mass balance with altitude. Once the model is tuned, the effect of temperature change is simulated by recalculating the mass balance with the same parameters as before, but with a temperature increase of 1°C throughout the year. The largest mass-balance changes, involving increased ablation of > 1 m w.e. a−1 °C−1, occur at the snout, with a progressively smaller increase with altitude. The area-averaged sensitivities for the five glaciers are −0.7 to −0.9 m w.e. a−1 °C−1. If annual precipitation also increased by 20% it would partly offset the effect of the 1°C higher temperatures but could not compensate for it.

Single crystals with a wide variety of orientations were cut from large pucks of laboratory-grown ice. Constant-strain-rate compression tests were performed on the crystals either at an axial strain rate of 1 × 10−5 s−1at –20°C or at axial strain rates from 1 × 10−6 s−1 to 1 × 10−4 s−1 at –10°C. In agreement with previous studies of ice flow, the compression tests showed a linearly rising stress with increasing strain, followed by a sharply declining stress after reaching a peak. With further strain, the sharp decline in stress slowed and the flow stress approached a plateau that was only weakly dependent on strain. For all crystallographic orientations, it was found that Schmid’s (critical resolved shear stress) law was obeyed by the peak stress. Slip lines clearly showed that basal slip was the deformation mode.

Modeling the reflectance and transmittance of strong volume scatterers is a delicate task. Slightly different approaches can lead to different results, making comparisons difficult. Here a simple, analytic multiple-scattering model is presented as a possible reference for comparisons and also for better understanding of the physics involved. The model quantifies the transmittance and reflectance of homogeneously distributed scatterers within slabs of any thickness. The simplicity of the model is given by the one-dimensional geometry, a system consisting of freely arranged ice lamellae in air. Although direct application of the model will be limited, it gives a spectral description of ice clouds and snowpacks over a very broad spectral range from microwave to ultraviolet. As well as the transmittance and reflectance, the model gives the emittance through Kirchhoff’s law. Comparison with other models shows, on the one hand, agreement with current snow models in the spectral description, and on the other, some quantitative inconsistencies between all of them. It appears that the lamella pack produces the same optical spectra as an average snow model, with spherical ice grains whose radius corresponds to about the lamella thickness, whereas microwave spectra appear to be slightly different.

The evolution of an inter-ice-stream ridge flanked by stagnated ice streams is simulated using a finite-difference, continuity ice-flow model. The model tracks the elevation of small-scale topographic undulations on the ice surface (“scars”) which form at ice-stream margins, and shows that after ice-stream stagnation these surface features are lifted onto the flanks of the evolving ridge before they are carried downslope by ice flow. The model is applied to the stagnant ice streams bounding Siple Dome, West Antarctica: “Siple Ice Stream” (SIS) on the northeast flank near Ice Stream D, and the “Duckfoot” area (DF) on the south flank near Ice Stream C. The volume-adjustment time-scale corresponding to the evolution of Siple Dome and these stagnant ice-stream areas is 1500–2000 years. The present geometry and elevation of the scar features, in addition to measurements of the present mass flux across the ridge, are used to estimate stagnation ages for SIS and DF. These measurements suggest that both SIS and DF stagnated 200–500 years ago.

Numerous studies of the ice caps in Greenland and Antarctica have observed accumulations of transuranic radionuclides and fission products from nuclear weapons testing, particularly during the period 1945–75. Recently, the concentrations of radionuclides in the annually deposited surface layers of Agassiz Ice Cap, Ellesmere Island, Canadian Arctic, from 1945 to the present have been measured and have demonstrated a continuous record of deposition of 137Cs and 239,240Pu in ice and snow. In this study, 3He-ingrowth mass spectrometry has been used to measure the low levels of tritium (3H) in some of these samples. Pre-nuclear-bomb tritium levels in ice-core samples were approximately 12 TU in high-latitude meteoric waters and 3–9 TU in mid-latitude meteoric waters. Comparisons of 3H levels and 3H/137Cs + 239,240Pu ratios, which were quite low during the earliest fission-bomb detonations (1946–51) and substantially higher during thermonuclear hydrogen-fusion bomb testing (1952–64), provide a clear indication of the type of nuclear device detonated. This finding accords with the results from other ice-core studies of the distribution of anthropogenic radionuclides from bomb fallout.

We propose to quantify the climate sensitivity of the mean specific balance B of a glacier by a seasonal sensitivity characteristic (SSC). The SSC gives the dependence of B on monthly anomalies in temperature and precipitation. It is calculated from a mass-balance model. We show and discuss examples for Franz-Josef Glacier (New Zealand), Nigardsbreen (Norway), Hintereisferner (Austria), Peyto Glacier (Canadian Rockies), Abramov Glacier (Kirghizstan) and White Glacier (Canadian Arctic). With regard to the climate sensitivity of B, the SSCs clearly show that summer temperature is the most important factor for glaciers in a dry climate. For glaciers in a wetter climate, spring and fall temperatures also make a significant contribution to the overall sensitivity. The SSC is a 2 × 12 matrix. Multiplying it with monthly perturbations of temperature and precipitation for a particular year yields an estimate of the balance for that year. We show that, with this technique, mass-balance series can be (re)constructed from long meteorological records or from output of atmospheric models.

Buried layers of surface hoar are the failure plane for many slab avalanches, including fatal human-triggered avalanches in various mountain regions. These layers may persist as weak layers in the snow cover for weeks or months. It is therefore essential for operational avalanche forecasters to monitor the evolution of persistent weak layers, such as buried surface hoar. Traditional grain-shape observations of isolated grains with a magnifier and crystal screen do not show bonding that is decisive for strength. In this study we used in situ microphotography and observations of texture to complement strength measurements from shear frame tests. Buried layers of surface hoar consist of crystals most of which extend from the layer below to the layer above, and may exhibit a columnar or truss-like structure. Observations and measurements show that texture and crystal size change little over periods of up to several months during which the snowpack remains dry. Under these conditions, layer thickness decreases while density and strength increase. Based on field measurements, we argue that the increase in strength is primarily due to penetration of the surface-hoar crystals into the adjacent layers, especially at the bottom of the buried surface-hoar layer, where bonding is critical. The weak bonding at the bottom implies that shear failure occurs at the lower interface rather than within the weak layer. On slopes, we find that surface-hoar crystals that were initially surface-normal are tilted downslope faster than predicted by published shear strain rates for settled snow, indicating that shear strain is concentrated in these layers. The characteristic texture of buried surface hoar (columnar or truss-like) permits collapsing at the time of fracture. The gravitational energy released by the displacement of the slab may contribute to the extensive fracture propagation associated with buried surface-hoar layers.

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.

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.

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.

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.