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.

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.

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.

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.

The Menzies Range in the southern Prince Charles Mountains, Antarctica, records at least four intervals of Cenozoic terrestrial glacigene sedimentation, and two periods of glacial erosion. The oldest Cenozoic strata, here named the Pardoe Formation, are >240 m thick, and consist of variable diamicts with subordinate sandstones and minor laminated lacustrine siltstones. The Pardoe Formation overlies a rugged erosion surface cut into Precambrian basement. Two subsequent Cenozoic sequences are here named informally the Trail diamicts and the younger Amphitheatre diamicts. The latter infilled the lower regions of an extremely rugged erosion surface, many components of which still dominate the present topography. The palaeodrainage of this erosion surface is markedly discordant with that of the older erosion surface underlying the Pardoe Formation. These three depositional events and the two associated erosion surfaces record warmer climates and increased snow accumulation under conditions of temperate wet-based glaciation. During the excavation of the sub-Amphitheatre diamict erosion surface, the East Antarctic ice sheet was either absent, further inland or the height of its surface relative to the Menzies Range was considerably lower than at present. The fourth and youngest depositional episode, recorded by a veneer of boulder gravel distributed along the northern flank of the Menzies Range, is from dry-based glacier ice, and assumed to be <2.6 Myr.

We present a comparison of surface velocities in 1996, derived from interferometric synthetic aperture radar, with an estimate of the long-term, depth-averaged velocity based on the assumption of steady-state flow for both Pine Island Glacier and its neighbour, Thwaites Glacier, West Antarctica. The results show that the former was close to balance conditions in 1996 (within 9%). The ice flux and velocity of the latter is significantly different in magnitude and distribution from that required to maintain the basin in a state of mass balance over the long term. The balance flux was found to be 32 ± 19% less than the measured outgoing flux. We conclude that the mass imbalance and dramatic difference in the pattern of flow is most likely due to a substantial change in the flow regime of Thwaites Glacier in the recent past.

Aerial photographs from 1947 and 1966, satellite optical imagery from 1973 and 1980, and interferometric synthetic aperture radar (InSAR) data from 1992, 1996 and 2000 are employed to detect ice-shelf changes in Pine Island Bay, Antarctica. The front position of the fast-flowing central ice shelf did not migrate discernibly over the past 50 years. New cracks and rifts appeared in the 1990s, however, that reveal a major weakening of the ice shelf. At the grounding-line center, the ice shelf thinned 21 m in 8 years. The northern, slow-moving ice shelf also shows signs of decay: (1) its calving front is retreating at an accelerating rate; and (2) the ice shelf is slowly unpinning from its bedrock anchors. These changes are taking place in a region well beyond the temperature-dependent limit of viability of ice shelves, and hence differ from those observed along the Antarctica Peninsula. They are likely due to a change in oceanic forcing, not to a change in air temperature. One possibility is that the documented intrusion of warm circumpolar deep water on the continental shelf has increased basal melting compared to that required to maintain the ice shelf in a state of mass balance, and that this has triggered a general retreat of ice in this sector.

Nineteen former valley glaciers were reconstructed for their Last Glacial Maximum (LGM) extents in the northern Uinta Mountains, Utah, U.S.A. Mean equilibrium-line altitudes (ELAs) calculated by four methods (accumulation–area ratio, toe–headwall altitude ratio, lateral moraines and cirque floors) range from 3050 to 3300 m a.s.l. Modern mean summer temperatures (Ts) at the ELAs range from 8.7° to 11.2°C, while modern winter precipitation (P) ranges from 354 to 590 mm snow water equivalent (SWE). Based on the difference in elevation of mean ELAs across the range, LGM P values must have ranged from 940 to 3040 mm SWE, assuming the modern summer lapse rate was the same during the LGM. A Ts depression of 5.5°C is required for these precipitation values to plot in the range of modern ELA values. The reconstructed increase in P at the western end of the range is 10 times the modern increase, reflecting the influence of pluvial Lake Bonneville. Assuming ELA depression (ΔELA) resulted from this P increase and a uniform 5.5°C Ts decrease, the regional LGM ΔELA was approximately 900 m.

We compare satellite albedo images of Vatnajökull, Iceland, with mass-balance measurements for the years 1991–99. We find that the equilibrium line is mostly not visible when it is located above its position of the previous year(s). Equilibrium-line detection is further hindered by clouds and a gradual transition between ice and firn or snow. Consequently, firn-line elevation at the end of the melting season is not particularly useful for estimating the annual mass balance. Instead, we propose to study the mean albedo of the entire ice cap throughout the melting season so that all available information about the surface albedo is taken into account. The mean net potential global radiation, which can be estimated from the mean surface albedo alone, both depends on and influences summer melt. It also depends on winter precipitation and, integrated over the melting season, is found to relate linearly to the specific mass balance B (r = 0.87 and 0.94 for different outlets of Vatnajökull). B can be estimated quantitatively when this relation is known and qualitatively when it is not. The uncertainty in the satellite-derived value of B is 0.5–0.8 m w.e., which for Vatnajökull corresponds to about 27% of the interannual variability of B.

In 1972 and 1995, shallow ice temperatures were measured at identical locations and depths on polythermal McCall Glacier, Brooks Range, Alaska, U.S.A. Mean annual ice temperatures at 10 m depth have systematically increased by > 1 K for the ablation area (1400–1900 m), while closer to the firn area, where meltwater percolation and refreezing play a role, they remained approximately unchanged. Interpreting these findings in terms of climate change requires careful consideration of the observed thinning of the glacier, which causes lowering of the surface through an existing vertical temperature gradient. Such temperature gradients can be particularly large in the ablation areas of polythermal glaciers; on McCall Glacier they are on the order of 0.2 K m−1. We also study the evolution of a 75 m deep temperature profile measured at one location in 1972 using a one-dimensional heat-diffusion model. We find that this profile was in approximate equilibrium with the mean surface temperature extrapolated from the 1972 data. Using the observed rate of surface lowering and measured rates of vertical advection, we find that both the measured temperature change and the vertical temperature gradient at 10 m depth can be reproduced only if the mean annual surface temperature on McCall Glacier has increased by 1.1 ± 0.3 K between 1972 and 1995. This result is consistent with the observed trend toward more negative mass balances on the glacier in the 1990s.

Structural glaciological, sedimentological and geophysical techniques are used to provide new insight concerning the formation of band ogives and associated structures at Bas Glacier d’Arolla, Switzerland. Sedimentary stratification, crevasse traces and transverse foliation are identified as planar structures in the lower icefall and glacier tongue. Stratification and crevasse traces are progressively deformed into, and enhance, the transverse foliation found in the glacier tongue. Three-dimensional geometry has been defined using ground-penetrating radar, which portrays four main characteristics: (i) deep reflectors interpreted as the ice/bed interface, (ii) alternating reflection-rich and reflection-poor zones interpreted as ogives, (iii) up-glacier-dipping reflectors, interpreted as planar structures, and (iv) down-glacier-dipping reflectors of uncertain origin. At the glacier surface, each band ogive consists of a light and dark band. The dark bands contain more intense foliation which, on differential weathering, traps fine debris. Clasts and clear ice of basal character within dark ogive bands suggest that basal ice has been raised to the glacier surface. The most applicable model for the formation of band ogives at Bas Glacier d’Arolla is a refinement of Posamentier’s (1978) “reverse faulting” hypothesis. In this context, multiple shear zones are formed, through which basal ice is uplifted to the glacier surface to produce the dark, foliated ogive bands. This model fits observations reported from other glaciers with band ogives.

The numerical integration of the heat-conduction equation is one of the main components in a thermodynamic sea-ice model. The spatial resolution in the ice normally varies from a minimum of three layers up to a few tens of layers. The temporal resolution varies from a few minutes up to hours. In this paper the impact of numerical resolution on the prediction of a one-dimensional thermodynamic ice model is studied. Analytical solutions for idealized cases were derived and compared with the numerical results. For the full ice model, groups of simulations were made, applying average climatic weather-forcing data corresponding to the ice-freezing, ice-thermal equilibrium and ice warm-up seasons. Special attention was paid to the effect of model spatial resolution. Early in the freezing season, the influence of resolution on model predictions is not significant. When the shortwave radiation becomes large, its absorption within the ice or snow cover was found to modulate the effect of numerical resolution on predictions of ice temperature and surface heat fluxes (e.g. the model run with a coarse spatial resolution yielded large daily variations in surface temperature). Resolution also affects the in-ice temperature profile. For process studies, a two-layer scheme for the solar radiation penetrating into the ice is suitable for a fine-spatial-resolution ice model.

In the past, two-dimensional images of internal snow structure have been obtained through plane surface sections or thin sections. These techniques are time-consuming and necessarily destroy the snow specimen. Computed tomography (CT) allows similar images to be obtained, but in a more efficient and non-destructive manner. To demonstrate the methodology, a CT scanner was used to obtain cross-sectional images over time of a snow sample undergoing kinetic-growth metamorphism. Densities calculated from the CT images correlated well to density measured using a traditional method. A procedure was developed that allows the CT image to be converted to binary in an objective manner. Employing innovative stereological software, the microstructural properties (grain diameter, bond diameter, neck length and intercept length) of the snow were then measured from the two-dimensional CT images. The presented methodology provides significant improvements over previous techniques, requiring less time and labor to obtain high-quality microstructural data.

Helicopter-borne 135 MHz short-pulse radar profiles of the Bagley Ice Field, southern Alaska, U.S.A., have been analyzed to determine the nature of radio-wave propagation through crevassed, stratified, wet and draining firn. Diffractions characterize the firn horizons along the trunk, and reflections characterize those of overlying snow-fields. Dynamic stacking was used to form and determine firn depths, and unstacked diffraction analysis to determine firn-layer properties. Refractive indices range from 4.1 to 4.5 for the near-surface, from 1.7 to 2.1 at about 13–17 m depth along the main trunk, and to 2.6 to 58 m depth within a snowfield. Average trunk values correspond with volumetric water contents of about 0.09, which agree with values measured for other glaciers. The analysis of the airborne-recorded diffractions is confirmed by deriving the approximate refractive index of water from a supraglacial lake-bottom diffraction. The lack of snowfield diffractions suggests that crevassing and not firn structure caused the trunk diffractions. The reasonable values of the indices imply that the diffractions originated from single points or edges orthogonal to the profile transect, and they predict low interlayer transmission losses. The snowfield penetration suggests that several hundred meters might be penetrated in uncrevassed accumulation zones with improved system design.