A rigid mast, equipped with five multi-turn potentiometers, was bolted to the bedrock floor of a frontal cavity at Glacier de Tsanfleuron, Switzerland. Each potentiometer was linked by a thin cord to an anchor emplaced in the adjacent ice wall, enabling basal ice motion to be measured continually over a 6 day period in the late summer. Results indicate that basal ice velocity increased quasi-linearly away from the glacier bed, rising from 10.66 to 11.82 mm d−1 between 25 and ∼ 265 mm above the ice– bed interface. Extrapolation of this gradient indicates that ∼ 10.55 mm d−1 may have been accommodated by pure slip between the ice and the glacier bed. Basal ice motion is temporally variable, generally being characterized by tens of minutes to hours of little or no motion interspersed with rapid-motion “events” lasting for between < 120 s (a single measurement interval) and 360 s. These motion events were responsible for speeds of up to > 400 mm d−1 over individual measurement intervals. The magnitude–frequency distribution of these events is consistent with a pattern of infrequent and large “slips” initiated at the ice–bed interface, that are manifested as more frequent and smaller motion events some tens of cm above that interface.

We investigate how temperature gradient and initial density influence depth-hoar growth in snow and seek to better define the range of conditions under which cohesive, or hard, depth hoar forms. Samples of 400 kg m−3 sieved snow were exposed to temperature gradients of 20–80°C m−1, and samples of four different densities were exposed to a 40°C m temperature gradient. Following exposure to temperature gradients, penetrometer tests were made on samples to determine the presence of solid and/or hard depth hoar. Grain bond orientation was analyzed in section planes by two-dimensional stereological techniques where hard depth hoar developed. Results indicate that hard cohesive depth hoar forms from rounded-grain snow having a density of 400 kg m−3 or greater, following exposure to a temperature gradient of 20°C m−1 or greater. Hard depth hoar appears to consist of solid-type depth-hoar grains connected by necks, with vertically preferred directions of grain elongation and organization of grain-to-grain chains. This work corroborates Atikaya’s (1974) results, but extends his observation of formation of hard depth hoar to weaker temperature gradients for high-density snow. Our results also indicate that hard depth hoar is composed of faceted solid-type (anhedral) grains.

During an ice-tank experiment, samples were taken to study the processes of acquisition and alteration of the gas properties in young first-year sea ice during a complete growth–warming–cooling cycle. The goal was to obtain reference levels for total gas content and concentrations of atmospheric gases (O2, N2, CO2) in the absence of significant biological activity. The range of total gas-content values obtained (3.5–18 mL STP kg−1) was similar to previous measurements or estimates. However, major differences occurred between current and quiet basins, showing the role of the water dynamics at the ice–water interface in controlling bubble nucleation processes. Extremely high CO2 concentrations were observed in all the experiments (up to 57% in volume parts). It is argued that these could have resulted from two unexpected biases in the experimental settings. Concentrations in bubbles nucleated at the interface are controlled by diffusion both from the ice–water interface towards the well-mixed reservoir and between the interface water and the bubble itself. This double kinetic effect results in a transition of the gas composition in the bubbles from values close to solubility in sea water toward values close to atmospheric, as the ice cover builds up.

We investigate limitations of the one-dimensional elastic-beam model to detect grounding line and thickness of an ice shelf from a differential interferogram. Spatial limitations due to grounding-line curvature and variable ice thickness are analyzed by comparison with two-dimensional plate flexure. Temporal limitations from the tide-dependent shift of the grounding line are analyzed by superpositions of four tidal flexure profiles representing differential interferograms. (i) At scales greater than one ice thickness, seaward protrusions of the grounding line are well represented by the elastic-beam model, while landward embayments of the same scale produce significant misplacements >10% of the ice thickness. (ii) For reasonable spatial variations of shelf thickness, the elastic-beam model gives reliable estimates of grounding-line position and unfractured mean ice thickness near the grounding line. (iii) For about 20% of superpositions of four tidal flexure profiles, the resulting grounding-line misplacements exceed the physical tidal shift of the grounding line by factors >2. For differential tide levels <10% of a 1 m tide dynamics, a physical shift of the grounding line of 0.3 km per metre of tide can lever misplacements of >2 km. Examples of real interferometric profiles from West Antarctic ice shelves corroborate our results.

We tested whether spaceborne interferometric synthetic aperture radar (InSAR) could be used to reveal patterns of redistribution of wind-drifted snow in arctic Alaska. Based on a simple geometric model, we found that lateral variations in new-snow (assuming a density of 0.3 g cm−3) accumulation of > 11 cm, or redistribution of the existing snow into dunes of half this height, could produce decorrelation of C-band interferograms. Comparison of interferograms with field observations for two periods from winter 1993/94, one with wind but little new snow, and the second with wind and new snow, indicates the interferograms delineated areas where the snow depth had changed due to drifting. Striking patterns of windward scouring and leeward deposition were revealed. The interferograms also showed that during one high-wind event, conspicuous interferometric bands a few kilometers wide and 30 km long were formed downwind of a mountain ridge. We speculate that these bands were caused by large-scale alterations in the concentrations of moving snow particles, a finding consistent with ground observations of alternating bands of clear air and blizzard for the same area and similar to phenomena observed with the Advanced Very High Resolution Radiometer and the Geodetic Earth Observing Satellite during blizzards in North Dakota and Iowa, U.S.A.

Measurements of the electrical conductivity of subglacial water provide a useful complement to measurements of pressure and turbidity. In the summer season, fluctuations of conductivity can be attributed to changes in water transport, water provenance and subglacial residence time. These explanations are unlikely to apply during the winter season because surface melt sources are not active and the subglacial water system is predominantly unconnected. Thus, fluctuations in water conductivity during the winter months seem paradoxical. To introduce a quantitative basis for comprehending such phenomena, we develop an interpretative model of the hydrochemical interaction between a water-filled borehole and a subglacial aquifer. The electrical conductivity of water near the borehole–aquifer contact is affected not only by diffusion but also by advective transport of solute between the two reservoirs in response to pressure forcing of the system. Using records of ice strain, water pressure and electrical conductivity from unconnected boreholes in Trapridge Glacier, we demonstrate that changes in borehole geometry caused by ice-strain events provide a plausible mechanism for at least some of the observed fluctuations of electrical conductivity. Conductivity records provide information regarding advective coupling of the borehole–aquifer system that is not available from pressure records alone.

Sediment production at a terrestrial section of the ice-sheet margin in West Greenland is dominated by debris released through the basal ice layer. The debris flux through the basal ice at the margin is estimated to be 12–45 m3 m−1 a−1. This is three orders of magnitude higher than that previously reported for East Antarctica, an order of magnitude higher than sites reported from in Norway, Iceland and Switzerland, but an order of magnitude lower than values previously reported from tidewater glaciers in Alaska and other high-rate environments such as surging glaciers. At our site, only negligible amounts of debris are released through englacial, supraglacial or subglacial sediment transfer. Glaciofluvial sediment production is highly localized, and long sections of the ice-sheet margin receive no sediment from glaciofluvial sources. These findings differ from those of studies at more temperate glacial settings where glaciofluvial routes are dominant and basal ice contributes only a minor percentage of the debris released at the margin. These data on debris flux through the terrestrial margin of an outlet glacier contribute to our limited knowledge of debris production from the Greenland ice sheet.

The freezing of sea water to the base of an ice shelf can give rise to large patches of accumulated ice, a phenomenon known as marine ice. In this study a numerical method is presented for calculating the thickness of the marine-ice layer using an ice- shelf-ocean model. The present-day modeling paradigm of ice-shelf–ocean interaction usually involves the fixed specification of the ice-shelf geometry while the ocean circulation in the cavity beneath the ice shelf evolves freely. This approach relies on several assumptions, such as steady-state ice-shelf thickness and ice-shelf flow fields, in order to make reasonable quantitative estimates of the thermodynamic exchange processes occurring at the ice-shelf base. This paper discusses the impact of these and other assumptions on the estimation of the thickness of the marine-ice layer. Model simulation results are presented for an idealized ice-shelf–ocean configuration as a demonstration of the feasibility of the numerical method. A sensitivity analysis is given so as to quantify the relative uncertainty in the marine-ice thickness that arises from uncertainties in the model input parameters, these being principally the ice-shelf flow field, the basal accumulation rate and the ice-shelf thickness field.

The transition from inland- to streaming-style ice flow near to and upstream from the onset to Ice Stream D, West Antarctica, is investigated using the force-balance technique. Basal drag provides the majority of the flow resistance over the study area but is substantially modified by non-local stress gradients. Lateral drag increases with distance downstream, balancing ∼50–100% of the driving stress at the onset. Longitudinal stress gradients (LSG) are also found to be significant, an observation that distinguishes ice flow in this region from the inland- and streaming-flow regimes that bound it, in which LSG are usually negligible. LSG decrease the spatial variability in basal drag and sliding speed and increase the area of the bed over which frictional melting occurs. Overall, LSG decrease the resistive influence of basal stress concentrations and increase the spatial uniformity of basal sliding. These observations suggest that streaming flow develops as an integrated response to the physical interaction between the ice and its bed over an extended region upstream from the onset, rather than being solely due to changes in basal characteristics at the onset. An implication is that non-steady-flow behavior upstream from the onset may ultimately propagate downstream and result in non-steady behavior at the onset.

To investigate the spatial distribution of the energy- and mass-balance fluxes of a glacier, a two-dimensional mass-balance model was developed and applied to Morteratschgletscher, Switzerland. The model is driven by meteorological input from four synoptic weather stations located in the vicinity of Morteratschgletscher. The model results were compared to observations made on the glacier. The calculated mean specific mass balance is −0.47 m w.e. for 1999, and 0.23 m w.e. for 2000. Net shortwave radiation shows a minimum at around 3350 m a.s.l., due to the effects of shading, slope, aspect, reflection from the slopes, and obstruction of the sky. Ignoring these effects results in a 37% increase in the annual incoming shortwave radiation on the glacier, causing 0.34 m w.e. more ablation. A 1°C change in the air temperature results in a shift of 0.67 m w.e. in the mean specific mass balance, while altering the precipitation by 10% causes a change of 0.17 m w.e.

At a site on the ice sheet adjacent to the Jakobshavn ice stream in West Greenland, ice deformation rates and temperatures have been measured in boreholes to the bedrock at 830 m depth. Enhanced deformation rates were recorded just below the Holocene–Wisconsin transition at 680 m depth. A 31 m layer of temperate ice and the temperature minimum of −22°C at 520 m depth were detected. The good agreement of these data with results of a two-dimensional thermomechanically coupled flow model implies that the model input is adequate. Discrepancies between modelled and measured temperature profiles on a flowline at the ice-stream centre have been attributed to effects not accounted for by the model. We have suggested that the convergent three-dimensional flow leads to a vertical extension of the basal ice entering the stream. A thick basal layer of temperate and Wisconsin ice would explain the fast flow of this ice stream. As a test of this hypothesis, the new core-borehole conductivity (CBC) method has been used to compare conductivity sequences from the ice stream to those of the adjacent ice sheet. The correlation thus inferred suggests that the lowest 270 m of the ice sheet correspond to the lowermost 1700 m of the stream, and, consequently, that the lower part of the ice stream has experienced a very large vertical extension.

Seasonal and interannual variations in surface elevation at the Greenland summit are modeled using a new temperature-dependent formulation of firn densification and are compared with elevations from European Remote-sensing Satellite (ERS-1/-2) radar altimetry. The rate constant and activation energy, usually set as constants in the Arrhenius-type relation, are strongly temperature-dependent, based on measurements of crystal-growth rates. A multiplicative factor in the densification rate accounts for differences between densification and grain-growth rates and is chosen to match the modeled and measured density profiles from 0 to 40 m. The stronger temperature dependence produces a significant seasonal cycle in the densification rate in the upper firn. Much of the densification and consequent surface lowering occur within 3 months in late spring/early summer, followed by a build-up from accumulation. Modeled elevation changes, using automatic weather station measurements of temperature and accumulation and modeled precipitation, agree well with observations. The respective seasonal amplitudes are 18 and 25 cm peak-to-peak with minima in mid-summer. The minimum elevation in 1995 is driven mainly by a temporary accumulation decrease and secondarily by warmer temperatures. Increased compaction driven by a summer warming trend decreases the modeled elevation (1992–99) by 20 cm, but accumulation increases in latter years dominate the overall 4.2 cm a−1 trend.

High rates of surface uplift and horizontal velocities were measured during a hydrologically induced spring speed-up event. Spatial patterns of surface uplift are analyzed to estimate components of vertical motion due to flow along an inclined bed and vertical strain. Areas are identified where surface uplift was most likely due in part to the opening or enlargement of subglacial cavities by bed separation. Results suggest a widespread enlargement of subglacial cavities during the event, and survival of residual cavities after the event. The spatial pattern of cavity enlargement closely matches previously identified axes of preferential subglacial drainage. It is suggested that localized cavity opening along axes of preferential drainage may constitute the initial stage in the seasonal development of channelized subglacial drainage. It is concluded that spatial and temporal variations in glacier motion may play an active role in determining the structure and rate of development of subglacial drainage during the summer melt season.

Deglaciated bedrock surfaces in limestone areas often exhibit extensive patterning by solutional furrows and carbonate deposits that occur in close association with undulations in the bed topography. These features clearly result from subglacial dissolution and precipitation of calcite on the bed — induced, for instance, by melting and freezing in a regelation water film — but little is known about the observed morphology. In particular, it is intriguing that (i) the solutional furrows, whose formation requires explanation, are collectively organized into arcuate patterns, with characteristic spacing, and (ii) a fluted or “spiculed” surface texture is ubiquitous on the calcite deposits. Herein, we propose specific mechanisms for such patterning based on a theory where chemical processes in the water film are coupled to regelation physics. Solutional furrows reflect locally enhanced dissolution along stoss surfaces, where CO2-rich bubbles advected in the ice from up-glacier come into contact with the bed. The bubbles form as CO2 is exsolved from freezing film water at the lee of bed bumps. The flutings on the deposit are inherently the manifestation of a spatial instability at the interface where calcite precipitation occurs. Complex interactions underlie some of the striking glacier-bed features shaped by subglacial chemical processes.

A major discontinuity in the variation of δ18O (δD) with altitude in high mountains was first seen in data from Mount Logan, YukonTerritory, Canada (Holdsworth and others, 1991). The profile of δ vs altitude revealed three well-defined regions: (1) a lower, monotonic, fractionation sequence below ∼3 km; (2) a middle layer, typically 1–2 km thick, within which δ values are nearly constant or stepped with altitude, and (3) part of another fractionation sequence in the “quasi-geostrophic flow region” above ∼5.3 km. The middle region was inferred to be a “mixed layer”, combining moisture from regions (1) and (3). This type of structure is now seen to occur on other high-altitude mountains, including Cerro Aconcagua, Argentina, where observations reach almost 7 km. The new observations confirm the general occurrence of a multi-layered atmosphere during precipitation at high-altitude glacier sites. This structure is linked to synoptic-scale polar cyclones, where the middle layer is identified as being the signature of the warm-front zone. These results have implications for the common practice of using a specific, spatially derived, isotopic thermometer in the time domain for the paleoclimatic interpretation of high-altitude ice-core δ records.

Stable-water-isotope data (δD and δ18O) from three groups of samples (fresh-snow and snow-pit samples collected on Qomolangma (Mount Everest) and Xixa-bangma during field seasons 1997, 1998 and 2001, and precipitation samples collected at Tingri station during summer 2000) are presented and used to survey the isotopic composition of precipitation over the northern slope of the central Himalaya. Multi-year snow-pit samples on Qomolangma have a local meteoric water-line (slope = 8) close to the global value. Deuterium excess (d = δD – 8δ18O) values at Tingri are much lower than those in fresh snow from Qomolangma, probably due to differences in moisture source and air-mass trajectories as well as local weather conditions. There is no obvious seasonal trend for d values in the Qomolangma region. A negative relationship exists between δ18O and d values in both fresh snow on Qomolangma and precipitation at Tingri. Fresh-snow samples collected from different altitudes on Xixabangma allow us to investigate the altitude effect on δ18O values in snow. Of four storm events, only one has an obvious altitude effect on δ18O variation and a very low gradient of −0.1 % per 100 m elevation.

We have developed a system for measuring a vertical strain-rate profile in the firn on polar ice sheets using a readily available video camera to detect metal bands inserted in an air-filled hole. We used this system in 1995 and 1996 at Taylor Dome, Antarctica. We use density measurements combined with our strain rates to infer vertical velocities. From our velocities we calculate a steady-state depth–age scale for the firn at Taylor Dome. The age of a visible ash layer from 79.1 m is 675 ± 25 years; this ash can be correlated with ash found at 97.2 m in a recent ice core at Siple Dome, West Antarctica.

Observations of the motion and basal conditions of Worthington Glacier, Alaska, U.S.A., during late-winter and spring melt seasons revealed no evidence of a relationship between water pressure and sliding velocity. Measurements included borehole water levels (used as a proxy for basal water pressure), surface velocity, englacial deformation, sliding velocity, and time-lapse videography of subglacial water flow and bed characteristics. The boreholes were spaced 10–15 m apart; six were instrumented in 1997, and five in 1998. In late winter, the water-pressure field showed spatially synchronous fluctuations with a diurnal cycle. The glacier’s motion was relatively slow and non-cyclic. In spring, the motion was characterized by rapid, diurnally varying sliding. The basal water pressure displayed no diurnal signal, but showed high-magnitude fluctuations and often strong gradients between holes. This transition in character of the basal water-pressure field may represent a seasonal evolution of the drainage system from linked cavities to a network of isolated patches and conduits. These changes occurred as the glacier was undergoing a seasonal-velocity peak. The apparent lack of correlation between sliding velocity and water pressure suggests that local-scale water pressure does not directly control sliding during late winter or early in the melt season.

As part of a larger program to measure and model vertical strain around Siple Dome on the West Antarctic ice sheet, we developed a new sensor to accurately and stably record displacements. The sensors consist of optical fibers, encased in thin-wall stainless-steel tubes, frozen into holes drilled with hot water, and stretched from the surface to various depths (up to 985 m) in the ice sheet. An optical system, connected annually to the fibers, reads out their absolute lengths with a precision of about 2 mm. Two sets of five sensors were installed in the 1997/98 field season: one set is near the Siple Dome core hole (an ice divide), and a second set is on the flank 7 km to the north (the ice thickness at both sites is approximately 1000 m). The optical-fiber length observations taken in four field seasons spanning a 3 year interval reveal vertical strain rates ranging from −229 ± 4 ppm a−1 to −7 ± 9 ppm a−1. In addition to confirming a non-linear constitutive relationship for deep ice, our analysis of the strain rates indicates the ice sheet is thinning at the flank and is in steady state at the divide.

Grain-size is an important but not well-known characteristic of snow at the surface of Antarctica. In the past, grain-size has been reported using various methods, the reliability, reproducibility and intercomparability of which is not warranted. In this paper, we present and recommend, depending on available logistical support, three techniques of snow-grain sampling and/or imaging in the field as well as an original digital image-processing method, which we have proved provides reproducible and intercomparable measures of a snow grain-size parameter, the mean convex radius. Results from more than 500 samples and 3000 images of snow grains are presented, which yield a still spatially limited yet unprecedentedly wide picture of near-surface snow grain-size distribution from fieldwork in Antarctica. In particular, except at sites affected by a very particular meteorology, surface grains in the interior of the ice sheet are uniformly small (0.1–0.2 mm). The climate-related increase of grain-size with depth through metamorphism is, as expected, not spatially uniform. Our Antarctic snow grain-size database will continue to grow as field investigations bring new samples, images and measures of snow grain.