A suite of subglacial water-pressure records from the 1996 summer field season at Trapridge Glacier, Yukon Territory, Canada, discloses a hydraulic event that cannot readily be explained by known forcings. We suggest that these records indicate covert failure of the pressure sensors caused by at least one large water-pressure pulse. The sign and magnitude of the pulse appears to have varied spatially and the pulse duration was less than the 2 min sampling interval of our data loggers. Laboratory experiments support this interpretation and indicate that the pulse magnitude exceeded 900 m of hydraulic head, roughly 15 times the ice-overburden pressure. Within glaciers, large water-pressure pulses can be generated when abrupt ice motion changes the volume of the subglacial hydraulic system.

The ability of synthetic-aperture radar (SAR) to record change in glaciers and icefields on seasonal to interannual time-scales is useful in maritime mountain regions where visible data are often obscured by clouds. A time series of RADARSAT and second European Remote-sensing Satellite (ERS-2) SAR images shows dramatic changes related to the onset and progression of glacier melting on the Juneau Icefield, southeast Alaska, U.S.A. Four “radar glacier zones” are interpreted from SAR images as snow that is dry during winter, an early snowmelt (M), a second phase of snowmelt (M2) and bare ice. These zones develop every year on the Juneau Icefield and other mid-latitude glaciers. Summer field observations on the Juneau Icefield during 1997 and 1998 constrain interpretations of the C-band radar glacier zones. Of the two zones that occur in melting snow (M, M2), M has low radar backscatter coefficients (σ° < −12), in contrast to the higher backscatter coefficients (σ° > −12) of the subsequent M2 zone. Snow moisture and surface roughness at the scale of the radar wavelength (5.6 cm) were measured to characterize the melt zones. Correlation length, wetness and grain-size in the two zones are not distinguishable in the late ablation season. Mean surface roughness, due to the presence of suncups, is higher in the M2 zone than in the M zone and probably causes the higher backscatter due to greater scattering.

The characteristics and sensitivities of a cold-based glacier on the Tibetan Plateau, where the summer monsoon provides most of the mass input to glaciers, are discussed using an energy-balance model incorporating the process of water refreezing. The model accurately represents the observational results related to the mass balance of Xiao Dongkemadi glacier on the central plateau during 1992/93. Our data revealed that the mass balance of cold glaciers cannot simply be described by the surface mass/heat balances, because about 20% of infiltrated water is refrozen and thus does not run off from the glacier. Model calculations demonstrate that glaciers in an arid environment can maintain their mass since the monsoon provides precipitation during the melting season. Snowfall in summer keeps surface albedo high and largely restrains ablation. Nevertheless, the calculations also make clear that glaciers on the plateau are more vulnerable than those of other regions because of summer accumulation. In the monsoon climate, warming would cause not only a decrease in accumulation, but also a drastic increase in ablation in combination with surface-albedo lowering. Therefore, although glaciers on and around the plateau can be sustained by summer accumulation, they are more vulnerable to warming than winter-accumulation-type glaciers.

During summer 1997, a 15 m firn core was recovered from Dasuopu glacier (28°23′ N, 85°44′ E; 7000 m a.s.l.) on the northwest margin of Xixabarngma Feng in the central Himalaya. Oxygen isotope values and concentrations of Ca2+, Mg2+, NH4+, SO42− and NO3− were measured over the 10 years of snow accumulation captured in the firn core. The seasonal variations of δ18O values and major-ion concentrations in the Dasuopu core indicate that summer monsoon and dust signals are clearly recorded in Dasuopu glacier. Annual variations in δ18O values are controlled by the amount effect, with more negative (i.e. lighter) δ18O values representing summer monsoon precipitation characteristic of tropical regions. Higher concentrations of Ca2+, Mg2+ and SO42− reflect the influx of mineral aerosols from the vast arid and semi-arid desert regions to the north and west during the spring dust-storm period. High spring concentrations of NH4+ and NO3− appear to reflect changes in regional biogenic-source strength.

General equations of ice dynamics are re-examined, using scale analysis, in order to derive a simplified thermomechanically coupled model for ice flow and heat transfer in ice caps filling volcanic craters. Relatively large aspect ratios between crater depths and diameters, low surface temperatures and intense volcanic heating are the principal characteristics of such craters. The conventional boundary-layer (shallow-ice) approximation is revised to account for these conditions and, in addition, the variable density of the snow, firn and bubbly ice. Large crater depths and intense bottom melting result in low longitudinal balance velocities, controlled by both shear and longitudinal stresses, and hence small surface slopes. In such situations ice can be assumed to be linearly viscous. A flowline model of the glacier dynamics is developed using this assumption. Explicit predictive formulas for ice-particle trajectories and age–depth relations, thus obtained, suggest that the age of ice at the bottom of glaciers in volcanic craters on Kamchatka Peninsula, Russia, may reach hundreds or thousands of years. Ice cores from these glaciers represent unique climatic and volcanic archives.

A new method of ice-core drilling uses an annulus of hot-water jets to melt out a cylindrical ice core. This lightweight device used in combination with a fast hot-water drill can quickly obtain ice cores from any depth.

This paper discusses results from the second phase of the European Ice Sheet Modelling Initiative (EISMINT). It reports the intercomparison of ten operational ice-sheet models and uses a series of experiments to examine the implications of thermomechanical coupling for model behaviour. A schematic, circular ice sheet is used in the work which investigates both steady states and the response to stepped changes in climate. The major finding is that the radial symmetry implied in the experimental design can, under certain circumstances, break down with the formation of distinct, regularly spaced spokes of cold ice which extended from the interior of the ice sheet outward to the surrounding zone of basal melt. These features also manifest themselves in the thickness and velocity distributions predicted by the models. They appear to be a common feature to all of the models which took part in the intercomparison, and may stem from interactions between ice temperature, flow and surface form. The exact nature of these features varies between models, and their existence appears to be controlled by the overall thermal regime of the ice sheet. A second result is that there is considerable agreement between the models in their predictions of global-scale response to imposed climate change.

Black Rapids Glacier is a 40 km long surge-type glacier in the central Alaska Range. In spring 1997 a wireline drill rig was set up at a location where the measured surface velocities are high and seasonal and annual velocity variations are large. The drilling revealed a layer of subglacial “till”, up to 7 m thick, that is believed to be water-saturated. At one location a string of instruments, containing three dual-axis tiltmeters and one piezometer, was successfully introduced into the till. The tiltmeters monitored the inclination of the borehole at the ice–till interface and at 1 and 2 m into the till, for 410 days. They showed that no significant deformation occurred in the upper 2 m of the till layer, and no significant amount of the basal motion was due to sliding of the ice over the till. The measured surface velocity at the drill site is about 60 m a−1, of which 20–30 m a can be accounted for by ice deformation. Almost the entire amount of basal motion, 30–40 m a−1, was taken up at a depth of > 2 m in the till, possibly in discrete shear layers, or as sliding of till over the underlying bedrock. We propose that the large-scale mobilization of such till layers is a key factor in initiating glacier surges.

Ice single crystals of various orientations containing various concentrations of H2SO4 up to 11.5 ppm were cut from large pucks of laboratory-grown ice. Constant-strain-rate compression tests were performed on the doped ice crystals both at −20°C at an axial strain rate of 1 × 10−5 s−1 and at −10°C at 1 × 106 s−1. The stress–strain curves 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. The tests clearly showed, for the first time, that this naturally occurring impurity dramatically decreases both the peak stress and the subsequent flow stress of ice single crystals. The decrease in the peak strength was related to the square root of the concentration of H2SO4 up to 11.5 ppm, suggesting that the solubility limit of H2SO4 in ice is at least 11.5 ppm. The sulfuric acid also appeared to increase the ductility of the ice. Preliminary examination of a doped ice single crystal by synchrotron X-ray topography suggested that sulfuric acid dramatically increases the grown-in dislocation density.

A numerical model of the ice-sheet/ice-shelf transition was used to investigate ice-sheet dynamics across the large subglacial lake beneath Vostok station, central East Antarctica. European Remote-sensing Satellite (ERS-1) altimetry of the ice surface and 60 MHz radio-echo sounding (RES) of the ice-sheet base and internal ice-sheet layering were used to develop a conceptual flowline across the ice sheet, which the model used as input. The model calculates horizontal and vertical velocities and stresses, from which particle flow paths can be obtained, and the ice-sheet temperature distribution. An inverse approach to modelling was adopted, where particle flow paths were forced to match those identified from internal RES layering. Results show that ice dynamics across the inflow grounding line are similar to an ice-sheet/ice-shelf transition. Model particle flow paths match internal RES layering when ice is (a) taken away from the ice base across the first 2 km of the flowline over the lake and (b) added to the base across the remainder of the lake. We contend that the process causing this transfer of ice is likely to be melting of ice and freezing of water at the ice–water interface. Other explanations, such as enhanced rates of accumulation over the grounding line, or three-dimensional convergent/divergent flow of ice are inconsistent with available measurements. Such melting and refreezing would be responsible for circulation and mixing of at least the surface layers of the lake water. Our model suggests that several tens of metres of refrozen “basal ice” would accrete from lake water to the ice sheet before the ice regrounds.

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.

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.

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.

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.

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.