This study provides a detailed structural glaciological analysis of changes in surface structures on the Larsen B ice shelf on the Antarctic Peninsula prior to its collapse in February–March 2002. Mapped features include the ice-shelf front, rifts, crevasses, longitudinal linear surface structures and meltwater features. We define domains on the ice shelf related to glacier source areas and demonstrate that, prior to collapse, the central Larsen B ice shelf consisted of four sutured flow units fed by Crane, Jorum, Punchbowl and Hektoria/Green/Evans glaciers. Between these flow units were ‘suture zones’ of thinner ice where the feeder glaciers merged. Prior to collapse, large open-rift systems were present offshore of Foyn Point and Cape Disappointment. These rifts became more pronounced in the years preceding break-up, and ice blocks in the rifts rotated because of the strong lateral shear in this zone. Velocity mapping of the suture zones indicates that the major rifts were not present more than about 20 years ago. We suggest that the ice shelf was preconditioned to collapse by partial rupturing of the sutures between flow units. This, we believe, was the result of ice-shelf front retreat during 1998–2000, reducing the lateral resistive stress on the upstream parts of the shelf and glacier flow units, ice-shelf thinning and pre-shelf-break-up glacier acceleration.

For three field seasons (2002/03, 2004/05, 2005/06) we have deployed a network of GPS receivers and seismometers around the tip of a propagating rift on the Amery Ice Shelf, East Antarctica. During these campaigns we detected seven bursts of episodic rift propagation. To determine whether these rift propagation events were triggered by short-term environmental forcings, we analyzed simultaneous ancillary data such as wind speeds, tidal amplitudes and sea-ice fraction (a proxy variable for ocean swell). We find that none of these environmental forcings, separately or together, correlated with rift propagation. This apparent insensitivity of ice-shelf rift propagation to short-term environmental forcings leads us to suggest that the rifting process is primarily driven by the internal glaciological stress. Our hypothesis is supported by order-of-magnitude calculations that the glaciological stress is the dominant term in the force balance. However, our calculations also indicate that as the ice shelf thins or the rift system matures and iceberg detachment becomes imminent, short-term stresses due to winds and ocean swell may become more important.

An annual net accumulation history of the high-elevation region of Devon Ice Cap, Nunavut, Canada, was reconstructed for the period 1963–2003 using five shallow firn cores. Annual net accumulation decreased significantly after 1989. To explain variability in the reconstructed annual net accumulation record, monthly and seasonal moisture-source probabilities were calculated for gridcells throughout the Arctic during 1979–2003. Seasonally, moisture-source probabilities reach a maximum in northern Baffin Bay in late summer/early fall and approach zero throughout the Arctic in winter. Late-summer/early-fall moisture-source probabilities were significantly higher around the North Open Water (NOW) Polynya during the 4 year period of highest annual net accumulation during the 1979–2003 period (1984–87), than during the 4 year period with the lowest annual net accumulation (1994–97). This is due to both a significant decrease in the sea-ice fraction and a significant increase in low-elevation atmospheric transport over the NOW area during the high net accumulation period. Anomalously low net accumulation and anomalously high firnification rates during the 1989–2003 period suggest that a change in ice dynamics, rather than a change in surface mass balance, may explain recent ice-cap thickening observed by laser altimetry.

Rutford Ice Stream and Carlson Inlet are neighbouring glaciers in West Antarctica. Rutford Ice Stream flows at speeds greater than 350 m a−1, whereas Carlson Inlet, which has some similar dimensions and supports a similar driving stress, flows 10–50 times slower. We discuss a range of observations concerning Carlson Inlet, and conclude that there is good indirect evidence that it is a relict ice stream, which ceased streaming more than 240 years BP, but sufficiently recently that its surface morphology, basal water content and basal morphology still retain characteristics produced by streaming. An analysis of expected subglacial drainage pathways indicates that Carlson Inlet is not streaming because it is currently starved of subglacial water, which is currently directed beneath Rutford Ice Stream. This current state of water piracy by Rutford Ice Stream is, however, sensitive to minor thickness changes on the ice streams; a ∼120 m (<4%) thickening of Rutford Ice Stream would divert almost all the subglacial water in the system towards Carlson Inlet and could reactivate its flow. The result highlights the importance of subglacial drainage in controlling ice-stream evolution and the requirement for ice-sheet models to couple ice flow with subglacial drainage.

Recent observations of increased discharge through fast-flowing outlet glaciers and ice streams motivate questions concerning the inland migration of regions of fast flow, which could increase drawdown of the ice-sheet interior. To investigate one process that could lead to inland migration we conduct experiments with a two-dimensional, full-stress, transient ice-flow model. An initial steady state is perturbed by initiating a jump in sliding speed over a fraction of the model domain. As a result, longitudinal-stress gradients increase frictional melting upstream from the slow-to-fast sliding transition, and a positive feedback between longitudinal-stress gradients, basal meltwater production and basal sliding causes the sliding transition to migrate upstream over time. The distance and speed of migration depend on the magnitude of the perturbation and on the degree of non-linearity assumed in the link between basal stress and basal sliding: larger perturbations and/or higher degrees of non-linearity lead to farther and faster upstream migration. Migration of the sliding transition causes the ice sheet to thin over time and this change in geometry limits the effects of the positive feedback, ultimately serving to impede continued upstream migration.

To investigate the internal accumulation on Storglaciären, Sweden, we couple a multi-layer snow model to a distributed energy- and mass-balance model. The snow model describes the temperature, density and water-content evolution of the snow/ice pack and includes the processes of percolation and refreezing of water. The model is run for the period 9 May–2 September 1999 and validated against weather station and mass-balance observations on the glacier. The model performs reasonably well, with an average summer mass balance for the location of stake observations of –1.56 m w.e. compared to –1.59 m w.e. observed. However, the amount of melt is overestimated in the higher parts of the accumulation area and underestimated in the lower parts of the ablation area. The modelled mass balance is most sensitive to the albedo parameterization, the chosen momentum and scalar roughness lengths and all parameters related to snowfall. The modelled internal accumulation is +0.25 m w.e., which amounts to about 20% of the winter accumulation and results in a positive net balance for 1999 of +0.23 m w.e. The modelled internal accumulation is most sensitive to the initial subsurface temperature profile and the irreducible water content.

Ground-penetrating radar (GPR) and satellite ERS-2 synthetic aperture radar (SAR) are used to map the thickness and extent of the superimposed ice (SI) zone on the surge-type glacier Kongsvegen, Svalbard. GPR imagery shows sub-horizontal SI layers lying unconformably above a discrete boundary. Below this boundary, the ice has a GPR signature similar to that of ice further down-glacier in the ablation zone. This boundary is posited to represent the closing of crevasses that were created during the last surge of Kongsvegen in ∼1948. Open crevasses would have interrupted the formation of sheet layers of SI due to efficient vertical drainage of the snowpack. Aerial photographs suggest that the crevasses closed sometime in the period 1956–66. A classified SAR image from 2003 is used to delineate the extent of the SI zone. The SI extent in the SAR image agrees well with the SI zone mapped by GPR. Using the SI spatial depth distribution, we estimate the mean annual accumulation of superimposed ice to be 0.16 ± 0.06 m w.e. a−1 (locally up to 0.43 m a−1 w.e.). This corresponds to ∼15–33% of the local winter balance and ∼5–10% of the total winter balance measured since 1987.

A daily time series of ‘clear-sky’ surface temperature has been compiled of the Greenland ice sheet (GIS) using 1 km resolution moderate-resolution imaging spectroradiometer (MODIS) land-surface temperature (LST) maps from 2000 to 2006. We also used mass-concentration data from the Gravity Recovery and Climate Experiment (GRACE) to study mass change in relationship to surface melt from 2003 to 2006. The mean LST of the GIS increased during the study period by ∼0.27°C a−1. The increase was especially notable in the northern half of the ice sheet during the winter months. Melt-season length and timing were also studied in each of the six major drainage basins. Rapid (<15 days) and sustained mass loss below 2000 m elevation was triggered in 2004 and 2005 as recorded by GRACE when surface melt begins. Initiation of large-scale surface melt was followed rapidly by mass loss. This indicates that surface meltwater is flowing rapidly to the base of the ice sheet, causing acceleration of outlet glaciers, thus highlighting the metastability of parts of the GIS and the vulnerability of the ice sheet to air-temperature increases. If air temperatures continue to rise over Greenland, increased surface melt will play a large role in ice-sheet mass loss.

Subglacial water threatens the stability of a deep ice sheet. The Greenland ice sheet has been the subject of extensive radar surveys under NASA’s PARCA program and others, yet the extent of basal water has not previously been explored in detail. This paper provides a method for recovering basal radar reflection characteristics that yield robust discrimination of subglacial water, using PARCA data, by narrowing the spread of measured echo intensities. Sources of validation are provided, and they support the procedure and its results. The method can be extended to the entire Greenland radar dataset. It is not dependent on details of radar characteristics but relies on sufficient penetration performance, high-fidelity data acquisition and appropriate treatment of signals and interface statistics.

The surface mass balance (SMB) at Dome Fuji, East Antarctica, was estimated using 36 bamboo stakes (grid of 6 × 6, placed at 20 m intervals) from 1995 to 2006. The heights of the stake tops from the snow surface were measured at 0.5 cm resolution twice monthly in 1995, 1996, 1997 and 2003, and once a year for the rest of the study period. To account for snow settling, the average snow density at the stake base during the measurements was used for converting the stake-height data to SMB. The annual SMB from 1995 to 2006 at Dome Fuji was 27.3 ± 1.5 kg m−2 a−1. This result agrees well with the annual SMB from AD 1260 to 1993 (26.4 kg m−2 a−1) estimated from volcanic signals in the Dome Fuji ice core. Over the period 1995–2006, there were 37 (8.6% of the measurements) negative or zero annual SMB results. Variation in the multi-year averages of annual SMB decreased with the square root of the number of observation years, and 10 years of observations of a single stake allowed the estimation of annual SMB at ±10% accuracy. The frequency distributions of annual and monthly SMB were examined. The findings clarify the complex behavior of the annual and monthly SMB at Dome Fuji, which will be common phenomena in areas of low snow accumulation of the interior of the Antarctic ice sheet.

The seasonal deuterium excess signal of fresh snow samples from Neumayer station, coastal Dronning Maud Land, Antarctica, was studied to investigate the relationship between deuterium excess and precipitation origin. An isotope model was combined with a trajectory model to determine the relative influence of different moisture sources on the mean annual course of the deuterium excess, focusing on the phase lag between δ18O and excess d. Whereas the annual course of δ18O always shows an austral summer maximum, which clearly depends on local temperature and the annual course of moisture source-area parameters, the deuterium excess of the fresh snow samples shows maximum values already in spring. There can be many different reasons for the time lag between δ18O and deuterium excess in an ice core, including post-depositional processes and changes in the moisture source of precipitation. The use of fresh snow samples enabled us to exclude post-depositional processes and study solely the influence of precipitation origin. Changes in the moisture source connected to systematic changes in the general atmospheric circulation can have a strong influence on the phase lag between deuterium excess and δ18O, which has to be taken into account for climatic interpretation of stable-isotope profiles from ice cores.

Based on the field data at Koryto glacier, Kamchatka Peninsula, Russia, we constructed a one-dimensional numerical glacier model which fits the behaviour of the glacier. The analysis of meteorological data from the nearby station suggests that the recent rapid retreat of the glacier since the mid-20th century is likely to be due to a decrease in winter precipitation. Using the geographical data of the glacier terminus variations from 1711 to 1930, we reconstructed the fluctuation in the equilibrium-line altitude by means of the glacier model. With summer temperatures inferred from tree-ring data, the model suggests that the winter precipitation from the mid-19th to the early 20th century was about 10% less than that at present. This trend is close to consistent with ice-core results from the nearby ice cap in the central Kamchatka Peninsula.

Rapid thinning and velocity increase on major Greenland outlet glaciers during the last two decades may indicate that these glaciers became unstable as a consequence of the Jakobshavn effect (Hughes, 1986), with terminus retreat leading to increased discharge from the interior and consequent further thinning and retreat. To assess whether recent trends deviate from longer-term behavior, we measured glacier surface elevations and terminus positions for Jakobshavn Isbræ, West Greenland, using historical photographs acquired in 1944, 1953, 1959, 1964 and 1985. These results were combined with data from historical records, aerial photographs, ground surveys, airborne laser altimetry and field mapping of lateral moraines and trimlines, to reconstruct the history of changes since the Little Ice Age (LIA). We identified three periods of rapid thinning since the LIA: 1902–13, 1930–59 and 1999–present. During the first half of the 20th century, the calving front appears to have been grounded and it started to float during the late 1940s. The south and north tributaries exhibit different behavior. For example, the north tributary was thinning between 1959 and 1985 during a period when the calving front was stationary and the south tributary was in balance. The record of intermittent thinning, combined with changes in ice-marginal extent and position of the calving front, together with changes in velocity, imply that the behavior of the lower parts of this glacier represents a complex ice-dynamical response to local climate forcings and interactions with drainage from the interior.

We conducted a seismic and radar survey of the central part of midtre Lovénbreen, a small, polythermal valley glacier in Svalbard. We determined the physical properties of the material beneath the glacier by measuring the reflection coefficient of the bed by comparing the energy of the primary and multiple reflections, and deriving the acoustic impedance. By making reasonable assumptions about the properties of the basal ice, we determined the acoustic impedance of the bed material as (6.78 ± 1.53) × 106 kg m−2 s−1. We interpret the material beneath the glacier to be permafrost with up to 50% ice, and we speculate that the material may be frozen talus similar to a deposit observed directly by others beneath another Svalbard glacier. The implication for midtre Lovénbreen is that the basal material beneath the present glacier is not able to support fast flow. We conclude that midtre Lovénbreen has most likely had limited capability for faster flow in the past, with motion dominated by internal deformation. Midtre Lovénbreen is used as a ‘study glacier’ for the scientific community in Svalbard, and a large number of studies have been based there. Our results show that it cannot be used as an analogue for larger glaciers in Svalbard, having distinct basal boundary conditions.

A passive seismic study was carried out underneath Glacier d’Argentière, Mont Blanc, France, where an array of seismometers was installed in a subglacial access tunnel. The data show a very high emissivity from the glacier. Fracturing can be discriminated from serac falls using the signal characteristics. We apply seismic array methods to locate the sources of these signals, using a two-step grid search in the parameter space. Four clusters of activity are found close to the network, showing that this fracturing does not take place uniformly over the glacier, but rather in isolated small zones. We compute a local magnitude using regional earthquakes for calibration. The magnitudes follow a classical Gutenberg–Richter law in the range ML = −3 to 0.15, showing that no characteristic size events dominate the process. We suggest that those spatial clusters of icequakes could reveal the heterogeneous nature of the friction at the base of the glacier, with patches of high frictional stresses locally generating intense fracturing within the ice mass.

If basal-water discharge and pressure are sufficiently high, a soft-bedded glacier will slip over its bed by ploughing, the process in which particles that span the ice–bed interface are dragged across the bed surface. Results of laboratory experiments indicate that resistance to ploughing can decrease with increasing ploughing velocity (velocity weakening). During ploughing at various velocities (15–400 m a−1), till was compacted in front of idealized particles, causing pore pressures there that were orders of magnitude higher than the ambient value. This excess pore pressure locally weakened the till in shear, thereby decreasing ploughing resistance by a factor of 3.0–6.6 with a six-fold increase in ploughing velocity. Characteristic timescales of pore-pressure diffusion and compaction down-glacier from ploughing particles depend on till diffusivity, ploughing velocity and sizes of ploughing particles. These timescales accurately predict the ranges of these variables over which excess pore pressure and velocity weakening occurred. Existing ploughing models do not account for velocity weakening. A new ploughing model with no adjustable parameters predicts ploughing resistance to no worse than 38% but requires that excess pore pressures be measured. Velocity weakening by this mechanism may affect fast glacier flow, sediment transport by bed deformation and basal seismicity.

We use the invariant Γ integral of fracture mechanics to calculate the frontal pressure and resistance to the downward motion a of snow avalanche. A basic characteristic property of the snowpack, termed the entrainment toughness, is introduced. From an analysis of the non-entrainment frictional mechanisms of avalanches, we find the necessary condition for a fracture-entrainment regime, and from an analysis of limiting equilibrium of gravitational force and frontal resistance, the necessary condition equation for the start of avalanches. We then derive the governing equations for the dynamics of avalanches, using a point-mass approach with entrainment taken into account. The governing equations are used to numerically simulate the Vallée de la Sionne (Switzerland) avalanche of 7 February 2003.