Radio-echo sounding data are used to investigate bed roughness beneath the three enhanced-flow tributaries of Slessor Glacier, East Antarctica. Slow-moving inter-tributary areas are found to have rough beds, while the bed of the northernmost tributary is relatively smooth. A reconstruction of potential subglacial drainage routing indicates that water would be routed down this tributary, and investigations of basal topography following isostatic recovery reveal that the bed would have been below sea level in preglacial times, so marine sediments may have accumulated here. Together, these factors are further support for the dominance of basal motion in this tributary, reported elsewhere. Conversely, although the other two Slessor tributaries may have water routed beneath them, they would not have been below sea level before the growth of the ice sheet, so cannot be underlain by marine sediments. They are also found to be rough, and, within the range of uncertainties, it is likely that basal motion does not play a major role in the flow of these tributaries. Perhaps the most interesting area, however, is a deep trough where flow rates are currently low but the bed is as smooth as the northern Slessor trough. It is proposed that, although ice deformation currently dominates in this trough, basal motion may have occurred in the past, when the ice was thicker.

We analyzed the profiles of ionic chemical species in three 500 mm sections of an ice core from Dome Fuji, Antarctica, dated 3.0, 8.9 and 13.3 kyr BP (before present), and compared the profiles to those in the surface snow. The 3.0 and 8.9 kyr sections are from the Holocene and the 13.3 kyr section slightly predates the Holocene. The analyses were done on 2 mm thick slices within each section. At each depth, the primary ionic species were Na+, H+, Cl– and SO42 A The SO42, Na+ and Mg2+ levels varied with depth in each section over distances ranging from several millimeters to several centimeters. Also, the correlation coefficients between Na+ and SO4 and between Mg2+ and SO42 for each depth were 0.90 or greater, in contrast to the value of 0.59 or less in the surface snow (defined here as 0–3.4 m from the surface). These results suggest that almost all Na+ and Mg2+ in the Holocene ice exists as Na2SO4 and MgSO4 salts, and the formation of these salts occurs not only in the atmosphere during transport, but also in the firn layer.

A newly developed hammer was used to insert two autonomous probes 0.8 m and 2.1 m into clast-rich subglacial till under Black Rapids Glacier, Alaska, USA. Both probes were instrumented with a dual-axis tilt sensor and a pore-water pressure transducer. The data are compared to a 75 day record of surface velocities. Till deformation at depth was found to be highly seasonal: it is significant during an early-season speed-up event, but during long periods thereafter measured till deformation rates are negligible. Both tilt records show rotation around the probe axis, which indicates a change in tilt direction of about 30°. The tilt records are very similar, suggesting spatial homogeneity on the scale of the probe separation (4 m horizontal and 3.3 m vertical). There is evidence that during much of the year sliding of ice over till or deformation of a thin till layer (<20 cm) accounts for at least two-thirds of total basal motion. Basal motion accounts for 50–70% of the total surface motion. The inferred amount of ice–till sliding is larger than that found at the same location in a previous study, when surface velocities were about 10% lower. We suggest that variations in ice–till coupling account for the observed variations in mean annual speed.

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.

The effect of spatial variations in ice thickness, accumulation rate and lateral flow divergence on radar-detected isochrone geometry in ice sheets is computed using an analytical method, under assumptions of a steady-state ice-sheet geometry, a steady-state accumulation pattern and a horizontally uniform velocity shape function. By using a new coordinate transform, we show that the slope of the isochrones (with a normalized vertical coordinate) depends on three terms: a principal term which determines the sign of the slope, and two scale factors which can modify only the amplitude of the slope. The principal term depends only on a local characteristic time (ice thickness divided by accumulation rate minus melting rate) between the initial and final positions of the ice particle. For plug flow, only the initial and final values have an influence. Further applications are a demonstration of how the vertical velocity profile can be deduced from sharp changes in isochrone slopes induced by abrupt steps in bedrock or mass balance along the ice flow. We also demonstrate ways the new coordinate system may be used to test the accuracy of numerical flow models.

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.

Wind tunnel measurements of snowdrift in a turbulent, logarithmic velocity boundary layer have been made in Davos, Switzerland, using natural snow. Regression analysis gives the drift threshold friction velocity (u*t), assuming an exponential drift profile and a simple drift to friction velocity relationship. Measurements over 15 snow covers show that u*t is influenced more by snow density and particle size than by ambient temperature and humidity, and varies from 0.27 to 0.69 ms–1. Schmidt’s threshold algorithm and a modified version used in SNOWPACK (a snow-cover model) agree well with observations if small bond sizes are assumed. Using particle hydraulic diameters, obtained from image processing, Bagnold’s threshold parameter is 0.18. Roughness lengths (z0) vary between snow covers but are constant until the start of drift. Threshold roughness lengths are proportional to . The influence of macroscopic objects on the roughness length is shown by the lower values measured over the smooth and flat snow surface of the wind tunnel (0.04 ≤ z0 ≤ 0.13 mm), compared to field measurements. Mean drifting-snow grain sizes for mainly new and partly decomposed snow are 100–175 μm, and independent of surface particle size.

Puncak Jaya, Irian Jaya, Indonesia, contains the only remaining tropical glaciers in East Asia. The extent of the ice masses on Puncak Jaya has been mapped from high-resolution IKONOS satellite images acquired on 8 June 2000 and 11 June 2002. Exclusive of Southwall Hanging Glacier, the ice extent on Puncak Jaya was 2.326 km2 and 2.152 km2 in 2000 and 2002, respectively. From 2000 to 2002, the Puncak Jaya glaciers lost a surface area of 0.174 km2 or 7.48% of their 2000 ice extent. Comparison of the IKONOS-based glacier extents with previous glacier extents demonstrates a continuing reduction of ice area on Puncak Jaya. By 2000, ice extent on Puncak Jaya had reduced by 88% of its maximum neoglacial extent. Between 1992 and 2000 Meren Glacier disappeared entirely. All remaining ice masses on Puncak Jaya continue their retreat from their neoglacial maxima. Comparison of 2000/2002 ice extents with previous extents suggests that these glaciers have not experienced accelerating rates of retreat during the last half of the 20th century. If the recession rates observed from 2000 to 2002 continue, the remaining ice masses on Puncak Jaya will melt within 50 years.

Complex glacier motion in the Grove Mountains region, Antarctica, is measured using four-pass differential synthetic aperture radar (SAR) interferometry (InSAR). The components of the motion vector field are resolved using a 44 day-separation Japanese Earth Resources Satellite-1 (JERS-1) InSAR pair and a European Remote-sensing Satellite-1/-2 (ERS-1/-2) tandem InSAR pair. The 44 day temporal baseline provides the sensitivity required to observe the range of ice motion (around 8–10ma–1), and the 1 day short baseline provides the best choice for glacier digital elevation model reconstruction. It is remarkable that the scattering field of the JERS-1 pair remained coherent over the long time interval and the interferometric fringes are clear. The overall ice flow is from east to west, downslope and towards Lambert Glacier. The regional flow is obstructed by nunataks extending north–south, with two wide gaps. Two narrow glaciers flow past the nunataks and rejoin each other at the downstream end. Regional morphology, and the resolved flow in the Grove Mountains area, suggests that ice flow is channeled throughout this eastern flank of Lambert Glacier.

We have made side-by-side measurements in several boreholes at Summit, Greenland, using borehole optical stratigraphy (BOS) and neutron-scattering density logging techniques. The BOS logs show strong positive correlation at shallow depths with neutron-scattering logs taken in the same borehole. This supports the hypothesis that BOS detects changes in density. The positive correlation between returned brightness and density decreases with depth and finally becomes negative. We conclude this inversion of correlation is related to changes in densification regime from grain-boundary sliding to pressure sintering.

Azimuth dependence of a normalized radar cross-section (σº) over the Greenland ice sheet is modeled with a simple surface scattering model. The model assumes that azimuth anisotropy in surface roughness at scales of 3–300 m is the primary mechanism driving the modulation. To evaluate the contribution of azimuth anisotropy in surface roughness to the radar backscatter, the model is compared to models based on isotropic surface roughness. The models are inverted to estimate snow surface properties using σº measurements from the C-band European Remote-sensing Satellite advanced microwave instrument in scatterometer mode. Results indicate that the largest mesoscale rms surface slopes are found in the lower portions of the dry snow zone. Estimates of the preferential direction in surface roughness are highly correlated with katabatic wind fields over Greenland, which is consistent with wind-formed sastrugi as the dominant mechanism causing azimuth modulation of σº. The maximum improvement of the azimuth modulation surface model compared to its isotropic counterparts occurs in the lower regions of the dry snow zone where the azimuth variability of σº is the largest. In regions with azimuth modulation over 1 dB, the mean root-mean-square error estimate of the azimuth-dependent surface scattering model is 0.46 dB compared with 0.70 dB for similar models using isotropic roughness.

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 sedimentary sequence of till overlying a gravel aquifer was instrumented with water-pressure transducers prior to a small, anticipated surge of the margin of the glacier Breiðamerkurjökull in Iceland. The records of water pressure at each transducer site show a well-defined temporal sequence of hydraulic regimes that reflect the changing recharge of surface-derived meltwater, the pressure drop along the drainage pathway and the pattern of ice loading. The poroelastic and water-pressure response of glacially overridden sediments to the recharge rate is determined in the frequency domain through an analytic solution. This permits the in situ conductivity, compressibility and consolidation states of subglacial sediments to be derived, and reveals aquifer-scale compressibility that produces an important water-pressure wave associated with the advancing glacier. The model is then used to explore how varying conductivity/compressibility, largely determined by granulometry, can determine drainage states and instabilities that may have a large impact on glacier/ice-sheet dynamics, and how the drainage time of surface water to the bed can determine the frequency response of subglacial groundwater regimes and their influence on subglacial sediment stability. Mismatches between model predictions and specific events in water-pressure records are used to infer processes that are not incorporated in the model: hydrofracturing that changes the hydraulic properties of subglacial sediments; the impact on groundwater pressure of subglacial channel formation; upwelling beyond the glacier margin; and rapid variations in the state of consolidation. The poroelastic model also suggests how seismic methods can be developed further to monitor hydraulic conditions at the base of an ice sheet or glacier.

This paper assesses a two-dimensional, vertically integrated ice model for its numerical properties in the calculation of ice-sheet evolution on a sloping bed using the shallow-ice approximation. We discuss the influence of initial conditions and individual model parameters on the model’s numerical behaviour, with emphasis on varying spatial discretizations. The modelling results suffer badly from numerical problems. They show a strong dependence on gridcell size and we conclude that the widely used gridcell spacing of 20 km is too coarse. The numerical errors are small in each single time-step, but increase non-linearly over time and with volume change, as a result of feedback of the mass balance with height. We propose a new method for the calculation of the surface gradient near the margin, which improves the results significantly. Furthermore, we show that we may use dimension analysis as a tool to explain in which situations numerical problems are to be expected.

We present a technique that modifies and extends down-hole target methods to provide absolute measures of uncertainty in radar-reflector depth of origin. We use ice-core profiles to model wave propagation and reflection, and then cross-correlate the model results with radio-echo sounding (RES) data to identify the depth of reflector events. Stacked traces recorded with RES near the EPICA drill site in Dronning Maud Land, Antarctica, provide reference radargrams, and dielectric properties along the deep ice core form the input data to a forward model of wave propagation that produces synthetic radargrams. Cross-correlations between synthetic and RES radargrams identify differences in propagation wave speed. They are attributed to uncertainties in pure-ice permittivity and are used for calibration. Removing conductivity peaks results in the disappearance of related synthetic reflections and enables the unambiguous relation of electric signatures to RES features. We find that (i) density measurements with g-attenuation or dielectric profiling are too noisy below the firn–ice transition to allow clear identification of reflections, (ii) single conductivity peaks less than 0.5 m wide cause the majority of prominent reflections beyond a travel time of about 10 µs (~900m depth) and (iii) some closely spaced conductivity peaks within a range of 1–2m cannot be resolved within the RES or synthetic data. Our results provide a depth accuracy to allow synchronization of age–depth profiles of ice cores by RES, modeling of isochronous internal structures, and determination of wave speed and of pure-ice properties. The technique successfully operates with dielectric profiling and electrical conductivity measurements, suggesting that it can be applied at other ice cores and drill sites.

Proportional straining experiments have been performed on columnar-grained S2 freshwater ice biaxially compressed across the columns at –10°C at a strain rate of (4.5 ± 1.5) × 10–3 s–1. The results are compared with those obtained earlier (Iliescu and Schulson, 2004) from the same kind of material deformed to terminal failure under the same conditions, but through proportional loading. The exercise shows that the biaxial strength is practically independent of the path taken, at least under low confinement where Coulombic shear faulting limits terminal failure. First-year sea ice is expected to exhibit the same behavior.

Glaciers occupy an area of ~1600km2 in the Caucasus Mountains. There is widespread evidence of retreat since the Little Ice Age, but an up-to-date regional assessment of glacier change is lacking. In this paper, satellite imagery (Landsat Thematic Mapper and Enhanced Thematic Mapper Plus) is used to obtain the terminus position of 113 glaciers in the central Caucasus in 1985 and 2000, using a manual delineation process based on a false-colour composite (bands 5, 4, 3). Measurements reveal that 94% of the glaciers have retreated, 4% exhibited no overall change and 2% advanced. The mean retreat rate equates to ~8ma–1, and maximum retreat rates approach ~38 m a–1. The largest (>10 km2) glaciers retreated twice as much (~12ma–1) as the smallest (<1 km2) glaciers (~6ma–1), and glaciers at lower elevations generally retreated greater distances. Supraglacial debris cover has increased in association with glacier retreat, and the surface area of bare ice has reduced by ~10% between 1985 and 2000. Results are compared to declassified Corona imagery from the 1960s and 1970s and detailed field measurements and mass-balance data for Djankuat glacier, central Caucasus. It is concluded that the decrease in glacier area appears to be primarily driven by increasing temperatures since the 1970s and especially since the mid-1990s. Continued retreat could lead to considerable changes in glacier runoff, with implications for regional water resources.

Automatic c-axes analyzers have been developed over the past few years, leading to a large improvement in the data available for analysis of ice crystal texture. Such an increase in the quality and quantity of data allows for stricter statistical estimates. The current textural parameters, i.e. fabric (crystallographic orientations) and microstructure (grain-boundary networks), are presented. These parameters define the state of the polycrystal and give information about the deformation undergone by the ice. To reflect the findings from automatic measurements, some parameter definitions are updated and new parameters are proposed. Moreover, a MATLAB® toolbox has been developed to extract all the textural parameters. This toolbox, which can be downloaded online, is briefly described.

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.