This paper introduces a new way to investigate in situ processes, the wireless multi-sensor probe, as part of an environmental sensor network. Instruments are housed within a ‘probe’ which can move freely and so behave like a clast. These were deployed in the ice and till at Briksdalsbreen, Norway. The sensors measure temperature, resistivity, case stress, tilt angle and water pressure and send their data to a base station on the glacier surface via radio links. These data are then forwarded by radio to a reference station with mains power 2.5 km away, from where they are sent to a web server in the UK. The system deployed during 2004/05 was very successful and a total of 859 probe days worth of data from the ice and till were collected, along with GPS, weather and diagnostic data about the system.

Past accumulation rate can be estimated from the measured number-density of bubbles in an ice core and the reconstructed paleotemperature, using a new technique. Density increase and grain growth in polar firn are both controlled by temperature and accumulation rate, and the integrated effects are recorded in the number-density of bubbles as the firn changes to ice. An empirical model of these processes, optimized to fit published data on recently formed bubbles, reconstructs accumulation rates using recent temperatures with an uncertainty of 41% (P < 0.05). For modern sites considered here, no statistically significant trend exists between mean annual temperature and the ratio of bubble number-density to grain number-density at the time of pore close-off; optimum modeled accumulation-rate estimates require an eventual ~2.02 ± 0.08 (P < 0.05) bubbles per close-off grain. Bubble number-density in the GRIP (Greenland) ice core is qualitatively consistent with independent estimates for a combined temperature decrease and accumulation-rate increase there during the last 5 kyr.

This work presents a method of mapping deformation-related sublimation patterns, formed on the surface of ice specimens, at microscopic resolution (3–4 μm pixel−1). The method is based on the systematic sublimation of a microtomed piece of ice, prepared either as a thick or a thin section. The mapping system consists of an optical microscope, a CCD video camera and a computer-controlled xy-stage. About 1500 images are needed to build a high-resolution mosaic map of a 4.5 × 9 cm section. Mosaics and single images are used to derive a variety of statistical data about air inclusions (air bubbles and air clathrate hydrates), texture (grain size, shape and orientation) and deformation-related features (subgrain boundaries, slip bands, subgrain islands and loops, pinned and bulged grain boundaries). The most common sublimation patterns are described, and their relevance for the deformation of polar ice is briefly discussed.

Terminus fluctuation and area change in mountain glaciers in the middle Chinese Tien Shan over the past four decades are examined, based on Landsat multispectral scanner (MSS), Enhanced Thematic Mapper Plus (ETM+), Système Probatoire pour l’Observation de la Terre (SPOT) highresolution visible (HRV) imagery, topographical maps and China Glacier Inventory data using Geographical Information System (GIS) and remote-sensing techniques. This study shows that the glaciers in the study area have undergone continuous recession between 1963 and 2000. Over the past 37 years the total glacier area decreased from 55 × 106 m2 to 48 × 106 m2 (a 13% decrease). The glaciers at low elevations decreased more rapidly than those at high elevations. Seven of the ten large valley glaciers in the study area retreated by at least 150 ± 49 m during this period. The average rate of recession for these glaciers was 4–8 m a−1. The recession rate between 1986 and 2000 was larger than between 1963 and 1986. The glacier changes in the study area were similar to those in the Ürümqi river basin in the northern Chinese Tien Shan. It is apparent that temperature rise over the past four decades has had a profound impact on these glaciers.

Ice-marginal Hidden Creek Lake, Alaska, USA, outbursts annually over the course of 2–3 days. As the lake fills, survey targets on the surface of the ‘ice dam’ (the glacier adjacent to the lake) move obliquely to the ice margin and rise substantially. As the lake drains, ice motion speeds up, becomes nearly perpendicular to the face of the ice dam, and the ice surface drops. Vertical movement of the ice dam probably reflects growth and decay of a wedge of water beneath the ice dam, in line with established ideas about jökulhlaup mechanics. However, the distribution of vertical ice movement, with a narrow (50–100 m wide) zone where the uplift rate decreases by 90%, cannot be explained by invoking flexure of the ice dam in a fashion analogous to tidal flexure of a floating glacier tongue or ice shelf. Rather, the zone of large uplift-rate gradient is a fault zone: ice-dam deformation is dominated by movement along high-angle faults that cut the ice dam through its entire thickness, with the sense of fault slip reversing as the lake drains. Survey targets spanning the zone of steep uplift gradient move relative to one another in a nearly reversible fashion as the lake fills and drains. The horizontal strain rate also undergoes a reversal across this zone, being compressional as the lake fills, but extensional as the lake drains. Frictional resistance to fault-block motion probably accounts for the fact that lake level falls measurably before the onset of accelerated horizontal motion and vertical downdrop. As the overall fault pattern is the same from year to year, even though ice is lost by calving, the faults must be regularly regenerated, probably by linkage of surface and bottom crevasses as ice is advected toward the lake basin.

Mass balances were measured on four glaciers in Taylor Valley, Antarctica, from 1993 to 2001. We used a piecewise linear regression, which provided an objective assessment of error, to estimate the mass balance with elevation. Missing measurements were estimated from linear regressions between points and showed a significant improvement over other methods. Unlike temperate glaciers the accumulation zone of these polar glaciers accumulates mass in summer and winter and the ablation zone loses mass in both seasons. A strong spatial trend of smaller mass-balance values with distance inland (r2 = 0.80) reflects a climatic gradient to warmer air temperatures, faster wind speeds and less precipitation. Annual and seasonal mass-balance values range only several tens of millimeters in magnitude and no temporal trend is evident. The glaciers of Taylor Valley, and probably the entire McMurdo Dry Valleys, are in equilibrium with the current climate, and contrast with glacier trends elsewhere on the Antarctic Peninsula and in temperate latitudes.

Representation of till rheology in glaciological models of ice motion over deformable sediments has, until now, focused largely on two end-member cases: (1) linear, or mildly non-linear, viscous rheology and (2) (nearly) plastic rheology. Most laboratory and in situ experiments support the latter model. Hindmarsh (1997) and Fowler (2002, 2003) proposed that experimental results represent the behavior of small till samples (characteristic length scales of ~0.1 to ~1 m) but that till behaves viscously over length scales that are relevant to determination of ice-flow rates in glaciers and ice sheets (~1 km or more). Observations of short speed-up events on the ice plain of Whillans Ice Stream, West Antarctica, provide an opportunity to compare the in situ rheology of this till, integrated over ~10–100 km, with the rheology of till from beneath the same ice stream determined on small laboratory samples and in local borehole experiments. This comparison indicates that the rheology of the subglacial till beneath Whillans Ice Stream is independent of scale.

A newly designed eight-arm caliper–inclinometer probe for borehole logging is presented. The caliper sounding gives eight points of cross-sections of the borehole at 1 cm intervals along its axis. Regression circles and ellipses are calculated and used as obvious hypotheses to be tested in order to obtain information on strain rates in the horizontal planes and isotropic or anisotropic borehole closure. The scratching of identifiable marks into the wall of the borehole enables material displacement to be tracked, in particular in the vertical direction, thus providing information for determining vertical strain rates along the borehole. Measurements in two boreholes drilled on high-altitude glaciers in the Swiss Alps are used to demonstrate both the potential and the limitations of the probe and of the mathematical methods presented.

Using a global climate model coupled with an ocean and a sea-ice model, we compare the effects of doubling CO2 and halving CO2 on sea-ice cover and connections with the atmosphere and ocean. An overall warming in the 2 × CO2 experiment causes reduction of sea-ice extent by 15%, with maximum decrease in summer and autumn, consistent with observed seasonal sea-ice changes. The intensification of the Northern Hemisphere circulation is reflected in the positive phase of the Arctic Oscillation (AO), associated with higher-than-normal surface pressure south of about 50° N and lower-than-normal surface pressure over the high northern latitudes. Strengthening the polar cell causes enhancement of westerlies around the Arctic perimeter during winter. Cooling, in the 0.5 × CO2 experiment, leads to thicker and more extensive sea ice. In the Southern Hemisphere, the increase in ice-covered area (28%) dominates the ice-thickness increase (5%) due to open ocean to the north. In the Northern Hemisphere, sea-ice cover increases by only 8% due to the enclosed land/sea configuration, but sea ice becomes much thicker (108%). Substantial weakening of the polar cell due to increase in sea-level pressure over polar latitudes leads to a negative trend of the winter AO index. The model reproduces large year-to-year variability under both cooling and warming conditions.

Crystal-orientation fabric (COF) has a large influence on ice-sheet flow. Earlier radar studies have shown that COF-based birefringence occurs within ice sheets. Radio-wave scattering in polar ice results from changing physical properties of permittivity and conductivity that arise from differing values of density, acidity and COF. We present an improved mathematical model that can handle all these phenomena together. We use this matrix-based model to study the two-way propagation of depolarized radio waves that scatter at both isotropic and anisotropic boundaries. Based on numerical simulations, we demonstrate how COF affects the radar signals in terms of radar polarization and frequency. We then compare the simulated features with VHF radar data obtained at two contrasting inland sites in East Antarctica, where COF is known from ice-core studies. These two sites are Dome Fuji, located near a dome summit, and Mizuho, located in a converging ice-flow region. Data at Dome Fuji are dominated by typical features resulting from birefringence. In contrast, both birefringence and anisotropic scattering affect the radar data at Mizuho. We argue that radar methods can be used to determine principal axes and strength of birefringence in the ice sheets.

The problem of describing ice dynamics has been faced by many researchers; in this paper a fully three-dimensional model for ice dynamics is presented and tested. Using an approach followed by other researchers, ice is considered a non-linear incompressible viscous fluid so that a fluid-dynamic approach can be used. The model is based on the full three-dimensional Stokes equations for the description of pressure and velocity fields, on the Saint-Venant equation for the description of the free-surface time evolution and on a constitutive law derived from Glen’s law for the description of ice viscosity. The model computes the complete pressure field by considering both the hydrostatic and hydrodynamic pressure components; it is time-evolutive and uses high-order numerical approximation for equations and boundary conditions. Moreover it can deal with both constant and variable viscosity. Three theoretical tests and two applications to Priestley Glacier, Antarctica, are presented in order to evaluate the performance of the model and to investigate important phenomena of ice dynamics such as the influence of viscosity on pressure and velocity fields, basal sliding and flow over perturbed bedrocks. All these applications demonstrate the importance of treating the complete pressure and stress fields.

Observations of surface motion and ice deformation from 2002–03 were used to infer mean stress fields in a cross-section of Black Rapids Glacier, Alaska, USA, over seasonal timescales. Basal shear stresses in a well-defined zone north of the center line (orographic left) were approximately 7% and 16% lower in spring and summer, respectively, than in winter. Correspondingly higher stresses were found near the margins. These changes in the basal shear stress distribution were sufficiently large to cause mean surface velocities to be 1.2 and 1.5 times larger in spring and summer than in winter. These results were inferred with a simple inverse finite-element flow model that can successfully reproduce bulk surface velocities and tiltmeter data. Stress redistribution between the well-defined zone and the margins may also occur over much shorter time periods as a result of rapidly changing basal conditions (ice–bed decoupling or enhanced till deformation), thereby causing large variations in surface velocity and strongly influencing the glacier’s net motion during summer.

Precise measurements of basal melting have been made at a series of 14 sites lying within a few kilometres of the grounding line of the Ronne Ice Shelf, Antarctica, where the ice thickness ranges from 1570 to 1940 m. The study was conducted over the course of 1 year and included a detailed survey of the horizontal deformation, as well as phase-sensitive radar measurements of the vertical displacement of both internal reflecting horizons and the ice-shelf base. Results from the surface survey show that the long-term viscous strain rate is modulated at tidal frequencies by (probably) elastic strains of order 10−5 per metre of tidal elevation. The radar measurements show a similar modulation of the long-term thinning/thickening of the ice shelf, with thickness oscillations up to a few centimetres in range. The long-term trends in ice thickness determined at points moving with the ice-shelf flow are consistent with a steady-state thickness profile. Vertical strain rates within the ice shelf were determined from the relative motion of internal reflectors. At two sites the observations were sufficient to discern the effect of tidal bending about a neutral surface 60% of the way down the ice column. Coincident measurements of horizontal and vertical strain imply a Poisson’s ratio of 0.5, and this combined with the asymmetric bending gives rise to the observed oscillations in thickness. At a number of sites the longterm viscous strain rates were found to be a linear function of depth. For an ice shelf this is an unexpected result. It can be attributed to the presence of significant vertical shear stresses set up close to the grounding line where the ice is still adjusting to flotation. Additional vertical motion arising from firn compaction was observed within the upper layers of the ice shelf. The additional motion was consistent with the assumption that firn density is a function only of the time since burial by steady surface accumulation. With both spatial and temporal fluctuations in the vertical strain rate accurately quantified it was possible to estimate the vertical motion of the ice-shelf base in response. Differences between the calculated and observed motion of the basal reflector arise because of basal melting. Derived melt rates at the 14 sites ranged from −0.11 ±0.31 to 2.51 ±0.10 m a−1, with a mean of 0. 85 m a−1 and a standard deviation of 0.69 m a−1, and showed no signs of significant sub-annual temporal variability. There was no obvious global correlation with either ice thickness or distance from the grounding line, although melt rates tended to decrease downstream along each of the flowlines studied. Previous estimates of basal melting in this region have been obtained indirectly from an assumption that the ice shelf is locally in equilibrium and have included a broad range of values. Only those at the lower end of the published range are consistent with the directly measured melt rates reported here.

Generalized numerical models of sub-debris ice ablation are preferable to empirical approaches for predicting runoff and glacier response to climate change, as empirical methods are site-specific and strongly dependent upon the conditions prevailing during the measurement period. We present a modified surface energy-balance model to calculate melt beneath a surface debris layer from daily mean meteorological variables. Despite numerous simplifications, the model performs well and modelled melt rates give a good match to observed melt rates, suggesting that this model can produce reliable estimates of ablation rate beneath debris layers several decimetres thick. This is a useful improvement on previous models which are inappropriate for thick debris cover.