The steep increase in Greenland’s glacial earthquake activity detected by the Global Seismographic Network since the late 1990s suggests that a close inspection of these events might provide clues to the nature and origin of such seismic activity. Here we discuss the detection of large, unexpected seismic events of extraordinarily long duration (10–40 min) occurring about once every 2 days, and localized in the ice stream that feeds the Earth’s fastest-moving glacier (Jakobshavn Isbræ) from the east. These ‘glacial rumblings’ represent an ice-mass wasting process that is greater and more frequent than glacial earthquakes have suggested. Probably triggered by calving, the rumblings are all very similar regardless of duration, and all end with a sharp, earthquake-like event in which the largest seismic amplitude is in the rumbling and that might signal the collapse of large ice masses upstream. By calculating the total amount of seismic energy released as rumblings, we estimate that the maximum seasonal amount of ice moved seismogenically down the ice stream is up to 12 km3, or ∼30% of the average annual iceberg discharge in Jakobshavn.

A structural glaciological description and analysis of surface morphological features of the Larsen C ice shelf, Antarctic Peninsula, is derived from satellite images spanning the period 1963–2007. The data are evaluated in two time ranges: a comparison of a 1963 satellite image photomosaic with a modern digital mosaic compiled using 2003/04 austral summer data; and an image series since 2003 showing recent evolution of the shelf. We map the ice-shelf edge, rift swarms, crevasses and crevasse traces, and linear longitudinal structures (called ‘flow stripes’ or ‘streak lines’). The latter are observed to be continuous over distances of up to 200 km from the grounding line to the ice-shelf edge, with little evidence of changes in pattern over that distance. Integrated velocity measurements along a flowline indicate that the shelf has been stable for ∼560 years in the mid-shelf area. Linear longitudinal features may be grouped into 12 units, each related to one or a small group of outlet feeder glaciers to the shelf. We observe that the boundaries between these flow units often mark rift terminations. The boundary zones originate upstream at capes, islands or other suture areas between outlet glaciers. In agreement with previous work, our findings imply that rift terminations within such suture zones indicate that they contain anomalously soft ice. We thus suggest that suture zones within the Larsen C ice shelf, and perhaps within ice shelves more generally, may act to stabilize them by reducing regional stress intensities and thus rates of rift lengthening.

Horizontal ice-core sites, where ancient ice is exposed at the glacier surface, offer unique opportunities for paleo-studies of trace components requiring large sample volumes. Following previous work at the Pâkitsoq ice margin in West Greenland, we use a combination of geochemical parameters measured in the ice matrix (δ18Oice) and air occlusions (δ18Oatm, δ15N of N2 and methane concentration) to date ice layers from specific climatic intervals. The data presented here expand our understanding of the stratigraphy and three-dimensional structure of ice layers outcropping at Pâkitsoq. Sections containing ice from every distinct climatic interval during Termination I, including Last Glacial Maximum, Bølling/Allerød, Younger Dryas and the early Holocene, are identified. In the early Holocene, we find evidence for climatic fluctuations similar to signals found in deep ice cores from Greenland. A second glacial–interglacial transition exposed at the extreme margin of the ice is identified as another outcrop of Termination I (rather than the onset of the Eemian interglacial as postulated in earlier work). Consequently, the main structural feature at Pâkitsoq is a large-scale anticline with accordion-type folding in both exposed sequences of the glacial–Holocene transition, leading to multiple layer duplications and age reversals.

Sound knowledge of the ice volume and ice-thickness distribution of a glacier is essential for many glaciological applications. However, direct measurements of ice thickness are laborious, not feasible everywhere and necessarily restricted to a small number of glaciers. In this paper, we present a method to estimate the ice-thickness distribution and the total ice volume of alpine glaciers. This method is based on glacier mass turnover and principles of ice-flow mechanics. The required input data are the glacier surface topography, the glacier outline and a set of borders delineating different ‘ice-flow catchments’. Three parameters describe the distribution of the ‘apparent mass balance’, which is defined as the difference between the glacier surface mass balance and the rate of ice-thickness change, and two parameters define the ice-flow dynamics. The method was developed and validated on four alpine glaciers located in Switzerland, for which the bedrock topography is partially known from radio-echo soundings. The ice thickness along 82 cross-profiles can be reproduced with an average deviation of about 25% between the calculated and the measured ice thickness. The cross-sectional areas differ by less than 20% on average. This shows the potential of the method for estimating the ice-thickness distribution of alpine glaciers without the use of direct measurements.

The resistance to sliding and the extent of till deformation beneath soft-bedded glaciers depend on the spatially averaged level of effective stress , which is controlled by the distribution of water pressure at the bed. Major subglacial conduits that facilitate large-scale water transport are expected to be predominantly aligned with the direction of maximum hydraulic gradient, which is normally parallel to the slope of the glacier surface. When the basal heat flow promotes net melting or freezing, seepage transport can enable water exchange between these conduits and the rest of the basal surface area. For a simple glacier geometry with subglacial conduits that are aligned parallel to a uniform slope, the seepage transport is driven primarily by gradients in effective stress. Balance equations determine how varies with conduit spacing and the heat-flow regime. Considerations of thermodynamic equilibrium require that ice penetrates the pore space at high effective stress. Even when the glacier base experiences net melting, for a given heat-flow regime there are limits on the conduit spacing that can be attained before a finite till layer becomes partially frozen throughout. During net freezing, the resistance to flow through partially frozen sediments limits the steady-state conduit spacing. The partially frozen zone can actually be restricted to smaller thicknesses when the freezing rate is greater.

Debris cover over glaciers greatly affects their rate of ablation and is a sensitive indicator of glacier health. This study focuses on estimation of debris cover over Samudratapu glacier, Chenab basin, Himalaya, using optical remote-sensing data. Remote-sensing image data of IRS-1C LISS-III (September 2001), IRS-P6 AWiFS (September 2004) and Terra ASTER (September 2004) along with Survey of India topographical maps (1963) were used in the study. Supervised classification of topographically corrected reflectance image data was systematically conducted to map six land-cover classes in the glacier terrain: snow, ice, mixed ice and debris, debris, valley rock, and water. An accuracy assessment of the classification was conducted using the ASTER visible/near-infrared data as the reference. The overall accuracies of the glacier-cover maps were found to range from 83.7% to 89.1%, whereas the individual class accuracy of debris-cover mapping was found to range from 82% to 95%. This shows that supervised classification of topographically corrected reflectance data is effective for the extraction of debris cover. In addition, a comparative study of glacier-cover maps generated from remote-sensing data (supervised classification) of September 2001 and September 2004 and Survey of India topographical maps (1963) has highlighted the trends of glacier depletion and recession. The glacier snout receded by about 756 m from 1963 to 2004, and the total glacier area was reduced by 13.7 km2 (from 110 km2 in 1963). Further, glacier retreat is found to be accompanied by a decrease in mixed ice and debris and a marked increase in debris-cover area. The area covered by valley rock is found to increase, confirming an overall decrease in the glacier area. The results from this study demonstrate the applicability of optical remote-sensing data in monitoring glacier terrain, and particularly mapping debris-cover area.

In order to develop and evaluate a method for the determination of glacier volume from ice-thickness data, the volume of Schaufelferner, Austria, is calculated (1) by manual interpolation of ground-penetrating radar (GPR) data based on measurements at 36 locations in 1995, (2) by manual interpolation of 144 GPR measurements acquired for a higher-resolution estimate in 2003 and 2006, (3) by multiplying the mean of the measured ice-thickness data by the glacier area, (4) by automatic kriging of the 1995 GPR data and (5) by application of area/volume scaling algorithms to the Austrian glacier inventory data of 1969, 1997 and 2006. The so determined glacier volumes are compared with the ice-volume changes calculated from digital elevation models (DEMs) of the Austrian glacier inventories. The manually interpolated volumes based on the 1995 and 2003/06 GPR data yielded a volume loss only slightly different from volume loss calculated from the glacier inventories of 1997 and 2007. Other methods were not able to reproduce the volume losses of the glacier inventory DEMs. To assess the accuracy of deriving ice-thickness changes with GPR, repeated ice-thickness measurements at the same locations were carried out between 2005 and 2008.

Subgrain boundaries revealed as shallow sublimation grooves on ice sample surfaces are a direct and easily observable feature of intracrystalline deformation and recrystallization. Statistical data obtained from the EPICA Dronning Maud Land (EDML) deep ice core drilled in East Antarctica cannot detect a depth region of increased subgrain-boundary formation. Grain-boundary morphologies show a strong influence of internal strain energy on the microstructure at all depths. The data do not support the classical view of a change of dominating recrystallization regimes with depth. Three major types of subgrain boundaries, reflecting high mechanical anisotropy, are specified in combination with crystal-orientation analysis.

An accurate three-dimensional reduced model (shallow-ice approximation) flow with velocity depending on all three spatial coordinates is constructed for the commonly adopted isotropic viscous law with temperature-dependent rate factor. The solution is for steady flow with a prescribed temperature distribution, but can be extended to flow with a coupled energy balance, and to unsteady flow. The accuracy hinges on the reduction to a two-point ordinary differential equation problem for the surface profile, on an unknown span, for which established accurate numerical methods are available. This is achieved by setting one horizontal velocity component in elliptic cylindrical coordinates to zero, but the other two components depend on all three spatial variables. While not of direct physical interest, such an ‘exact’ solution is valuable as a test solution for the large-scale numerical codes commonly used in ice-sheet modelling, which have not yet been subjected to such a comparison.

Previously, it was shown that the reduction in peak stress of single crystal ice at −20°C due to H2SO4, Δσ, could be described by Δσ = kC1/2, where C is the sulfuric acid concentration and k is a constant. Here we show that the strength reduction due to H2SO4 is more general and that the reduction is greater as the temperature decreases. Δσ was again found to be proportional to C1/2 at −10°C, but k at −10°C was significantly less than at −20°C.

The effects of temperature and seasonal air-mass trajectories on stable water isotopes in alpine snowpacks are investigated using meteorological and snow-pit data at two alpine field sites in the Canadian Rocky Mountains: Haig Glacier, Alberta, and Opabin Glacier, British Columbia. Snow pits were sampled through three accumulation seasons (October–June, 2004/05, 2005/06 and 2006/07) for δ18O, δD, temperature and density. The isotopic characteristics of precipitation over these time periods, including the local meteoric waterline and average δ18O, δD and deuterium excess, were defined using this dataset. Individual snowfall events over the three seasons were identified in the accumulation records from both sites and then fit to snow-pit stratigraphies to determine their mean isotopic characteristics. A trajectory classification was produced for all events, and the key meteorological characteristics of each trajectory class were investigated using data from alpine field sites and a suite of meteorological records from the region. An analysis of the relative influences of temperature and air-mass trajectory on snow isotope ratios reveals some separation in mean δ18O between storm classes. However, the separation appears to be driven primarily by the mean temperature of each class rather then being a direct effect of vapour pathway.

We describe a high-precision method, now in use in our laboratory, for measuring the CO2 mixing ratio of ancient air trapped in polar ice cores. Occluded air in ice samples weighing ∼8–15 g is liberated by crushing with steel pins at −35°C and trapped at −263°C in a cryogenic cold trap. CO2 in the extracted air is analyzed using gas chromatography. Replicate measurements for several samples of high-quality ice from the Siple Dome and Taylor Dome Antarctic ice cores have pooled standard deviations of <0.9 ppm. This high-precision technique is directly applicable to high-temporal-resolution studies for detection of small CO2 variations, for example CO2 variations of a few parts per million on millennial to decadal scales.

We describe the development of a low-frequency airborne radar specifically designed for the sounding of temperate ice. The system operates at a central frequency of 1 MHz and consists of an impulse transmitter with an output voltage up to 5000 V and a digital receiver with a maximum gain of 80 dB. The radar was deployed on board a CASA 212 aircraft, which also carries a laser altimeter, an inertial navigation system, a digital camera and a GPS receiver. A description of the radar system is provided, as well as preliminary results obtained at Glaciar Tyndall, Campo de Hielo Sur (Southern Patagonia Icefield), where an ice depth of 670 m was reached.

Recent work has shown that surface-to-bed drainage systems re-form annually on parts of the Greenland ice sheet and some High Arctic glaciers, leading to speed-up events soon after the onset of summer melt. Surface observations and geophysical data indicate that such systems form by hydrologically driven fracture propagation (herein referred to as ‘hydrofracturing’), although little is known about their characteristics. Using speleological techniques, we have explored and surveyed englacial drainage systems formed by hydrofracturing in glaciers in Svalbard, Nepal and Alaska. In Hansbreen, Svalbard, vertical shafts were followed through ∼60 m of cold ice and ∼10 m of temperate basal ice to a subglacial conduit. Deep hydrofracturing occurred at this site due to a combination of extensional ice flow and abundant surface meltwater at a glacier confluence. The englacial drainage systems in Khumbu Glacier, Nepal, and Matanuska Glacier, Alaska, USA, formed in areas of longitudinal compression and transverse extension and consist of vertical slots that plunge down-glacier at angles of 55° or less. The occurrence of englacial drainages initiated by hydrofracturing in diverse glaciological regimes suggests that it is a very widespread process, and that surface-to-bed drainage can occur wherever high meltwater supply coincides with ice subjected to sufficiently large tensile stresses.

Ice streams are a fact of ice-sheet dynamics, draining up to 90% of the ice. Thermal convection in ice below the density inversion is a speculation. An attempt is made to meld the two in such a way that the speculation becomes an explanation for the fact.

This paper aims to interpret the temporal variations of microwave brightness temperature (at 19 and 37GHz and at vertical and horizontal polarizations) in Antarctica using a physically based snow dynamic and emission model (SDEM). SDEM predicts time series of top-of-atmosphere brightness temperature from widely available surface meteorological data (ERA-40 re-analysis). To do so, it successively computes the heat flux incoming the snowpack, the snow temperature profile, the microwaves emitted by the snow and, finally, the propagation of the microwaves through the atmosphere up to the satellite. Since the model contains several parameters whose value is variable and uncertain across the continent, the parameter values are optimized for every 50 km × 50 km pixel. Simulation results show that the model is inadequate in the melt zone (where surface melting occurs on at least a few days a year) because the snowpack structure and its temporal variations are too complex. In contrast, the accuracy is reasonably good in the dry zone and varies between 2 and 4 K depending on the frequency and polarization of observations and on the location. At the Antarctic scale, the error is larger where wind is usually stronger, suggesting either that meteorological data are less accurate in windy regions or that some neglected processes (e.g. windpumping, surface scouring) are important. At Dome C, in calm conditions, a detailed analysis shows that most of the error is due to inaccuracy of the ERA-40 air temperature (∼2 K). Finally, the paper discusses the values of the optimized parameters and their spatial variations across the Antarctic.

In order to find environmental signals based on the dust and calcium-ion concentrations in ice cores, we determine the constituent elements of residue particles obtained after melting ice samples. We have designed a sublimating system that operates at −45°C, below the eutectic temperatures of major salts. This system permits us to obtain a great many non-volatile particles. After studying the non-volatile particles, we immersed them in water to remove soluble particles and compounds. We thereby analyzed a total of 1272 residue particles (from the melted sample), 2418 non-volatile particles (after sublimation) and 1463 insoluble particles taken from five sections of Last Glacial Maximum ice from the Dome Fuji (Antarctica) ice core. Their constituent elements were determined by scanning electron microscopy/energy-dispersive X-ray spectrometry (SEM-EDS) and compared to the dust, calcium-ion and sodium-ion concentrations measured by ion chromatography. Our results indicate that >99.9% of the insoluble particles contain silicon but no sulfur, nitrogen or chlorine. A significant number of the non-volatile particles, however, contain sulfur and chlorine. We conclude that insoluble dust consists mostly of silicate, that almost all calcium ions originate from calcium sulfate and that almost all sodium ions originate from sodium sulfate and sodium chloride.

The puzzling fact that Jakobshavn Isbræ, West Greenland, is flowing very fast but without any significant seasonal velocity changes, despite big amounts of surface-derived meltwater entering the ice stream (Echelmeyer and Harrison, 1990), has been explained by a combination of different types of measurements. It is now well established from seismic measurements and radio-echo sounding that Jakobshavn Isbræ flows through a deeply eroded subglacial trench that, even 50 km inland of the grounding line, extends as far as 1500 m below sea level (Clarke and Echelmeyer, 1996; Legarsky and Huang, 2006). Temperature measurements in boreholes down to 65% of the 2500 m thick ice stream at site B, some 50 km upstream of the calving front (Fig. 1b; Iken and others, 1993), were used to infer the presence of a substantial layer of temperate ice, the thickness of which was estimated to be at least 300 m by modeling and matching internal layering structures (Funk and others, 1994; Lüthi and others, 2002). The presence of a thick layer of temperate ice under very high driving stress allows for high ice-deformation rates, which contribute substantially to the observed fast flow velocities. Basal motion, while certainly important, seems to be barely influenced by the seasonal meltwater input (Echelmeyer and Harrison, 1990).

Clusters of surge-type glaciers occur in many areas peripheral to the Greenland ice sheet (Weidick, 1988). In East Greenland (68–72° N), 30–70% of glaciers are of surge type, but only five of these have been described (Jiskoot and others, 2003, and references therein). Regional surge characteristics and dynamics have been used to suggest a hydrologically controlled surge mechanism, with surge behaviour that is more Alaskan-type than Svalbard-type in this region (Jiskoot and others, 2001, 2003; Murray and others, 2002, 2003; Pritchard and others, 2005).