Any model of snow avalanches must be able to reproduce velocity profiles. This is a key problem in avalanche science because the profiles are the result of a multitude of snow/ice particle interactions that, in the fend, define the rheology of flowing snow. Recent measurements on real-scale avalanches show that the velocity profiles change from a highly sheared profile at the avalanche front to a plug-like profile at the avalanche tail, preventing the application of a single, simple rheology to the avalanche problem. In this paper, we model not only the velocity profiles but also the evolution of the velocity profiles, by taking into account the production and decay of the kinetic energy of the random motion of the snow granules. We find that the generation of this random energy depends on the distribution of viscous shearing within the avalanche. Conversely, the viscous shearing depends on the magnitude of the random energy and therefore its collisional dissipation. Thus, there is a self–consistency problem that must be resolved in order to predict the amount of random energy and therefore the velocity profiles. We solve this problem by stating equations that describe the production and decay of random energy in avalanches. An important guide to the form of these equations is that the generation of random energy is irreversible. We show that our approach successfully accounts for measured profiles in natural avalanches.

Seismometers deployed over a 3 year period on icebergs in the Ross Sea and on the Ross Ice Shelf, Antarctica, reveal that impulsive sources of ocean surface waves are frequent (e.g. ∼200 events per year in the Ross Sea) in the ice-shelf and iceberg-covered environment of coastal Antarctica. The 368 events recorded by our field deployment suggest that these impulsive events are generated by glaciological mechanisms, such as (1) small-scale calving and edge wasting of icebergs and ice-shelf fronts, (2) edge-on-edge closing and opening associated with iceberg collisions and (3) possibly the impulsive opening of void space associated with ice-shelf rifting and basal crevasse formation. The observations described here provide a background of glaciogenic ocean-wave phenomena relevant to the Ross Sea and suggest that these phenomena may be exploited in the future (using more purposefully designed observation schemes) to understand iceberg calving and ice-shelf disintegration processes.

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.

The Himalaya has some of the largest glacier concentrations outside the polar regions. Despite this, long-term measurements detecting the impact of global warming and changing precipitation patterns on glaciers are rare. The Nanga Parbat massif in northern Pakistan is an exception. The cartographer and glaciologist R. Finsterwalder investigated glacier dynamics of this mountain massif in the 1930s, and several other studies document changes since then. The aim of this study is to detect and analyse the changes of Raikot Glacier over the past seven decades. We use a multitemporal and multiscale approach, based on repeat terrestrial images, additional historical data and remotely sensed imagery (Corona, ASTER, Landsat, QuickBird). The multitemporal approach covers the period 1934–2007. While the analyses show a total glacier retreat of ~200 m in 73 years, this general trend was interrupted by a significant glacier advance between the 1950s and 1980s. Although down-wasting processes can be inferred from an increase in debris-covered area, a general trend of reduced glacier thickness does not appear significant over the whole observation period.

In the past, numerical prediction of regional avalanche danger using statistical methods with meteorological input variables has shown insufficiently accurate results, possibly due to the lack of snowstratigraphy data. Detailed snow-cover data were rarely used because they were not readily available (manual observations). With the development and increasing use of snow-cover models this deficiency can now be rectified and model output can be used as input for forecasting models. We used the output of the physically based snow-cover model SNOWPACK combined with meteorological variables to investigate and establish a link to regional avalanche danger. Snow stratigraphy was simulated for the location of an automatic weather station near Davos, Switzerland, over nine winters. Only dry-snow situations were considered. A variety of selection algorithms was used to identify the most important simulated snow variables. Data mining and statistical methods, including classification trees, artificial neural networks, support vector machines, hidden Markov models and nearest-neighbour methods were trained on the forecasted regional avalanche danger (European avalanche danger scale). The best results were achieved with a nearest-neighbour method which used the avalanche danger level of the previous day as additional input. A cross-validated accuracy (hit rate) of 73% was obtained. This study suggests that modelled snow-stratigraphy variables, as provided by SNOWPACK, are able to improve numerical avalanche forecasting.

Through the detection of surface deformation in response to water movement, recent satellite studies have demonstrated the existence of subglacial lakes in Antarctica that fill and drain on timescales of months to years. These studies, however, were confined to specific regions of the ice sheet. Here we present the first comprehensive study of these ‘active’ lakes for the Antarctic ice sheet north of 86° S, based on 4.5 years (2003–08) of NASA’s Ice, Cloud and land Elevation Satellite (ICESat) laser altimeter data. Our analysis has detected 124 lakes that were active during this period, and we estimate volume changes for each lake. The ICESat-detected lakes are prevalent in coastal Antarctica, and are present under most of the largest ice-stream catchments. Lakes sometimes appear to transfer water from one to another, but also often exchange water with distributed sources undetectable by ICESat, suggesting that the lakes may provide water to or withdraw water from the hydrologic systems that lubricate glacier flow. Thus, these reservoirs may contribute pulses of water to produce rapid temporal changes in glacier speeds, but also may withdraw water at other times to slow flow.

We present a comprehensive new inventory of surge-type glaciers on the Novaya Zemlya archipelago, using high-resolution (up to 4 m) satellite imagery from 1976/77 (Hexagon), 1989 (Landsat TM), 2001 (Landsat ETM+) and 2006 (ASTER). A total of 692 glaciers and their forelands were observed for glaciological and geomorphological criteria indicative of glacier surging (e.g. looped moraines, heavy surface crevassing, surface potholes, thrust-block moraines, concertina eskers). This enabled the identification of 32 potential surge-type glaciers (compared with four previously identified) representing 4.6% of the total but 18% by glacier area. We assess the characteristics of surge-type glaciers. Surge-type glaciers are statistically different from non-surge-type glaciers in terms of their area, length, surface slope, minimum elevation, mid-range elevation and terminus type. They are typically long (median length 18.5 km), large (median area 106.8 km2) outlet glaciers, with relatively low overall surface slopes (median slope 1.7°) and tend to terminate in water (marine or lacustrine). They are predominantly directed towards and located in the more maritime western region of the Russian Arctic, and we suggest that surge occurrence might be related to large and complex catchment areas that receive increased delivery of precipitation from the Barents Sea.

To contribute to the understanding of the influence of climate on glacial erosion and on orogenic processes, we report contemporary glacial erosion rates from a calving glacier in the Southern Andes and elucidate the influence of ice dynamics on erosion. Using seismic profiles of sediments collected in a proglacial fjord and a documented history of retreat, we determine the time-varying sediment flux of Glaciar Marinelli as a measure of basin-wide erosion rates, and compare these rates with the annual ice budget reconstructed using NCEP–NCAR reanalysis climate data from 1950 to 2005. The rate of erosion of the largest tidewater glacier in Tierra del Fuego averaged 39 ± 16 mm a−1 during the latter half of the 20th century, with an annual maximum approaching 130 mm a−1 following a decade of rapid retreat. A strong correlation emerges between the variable rate of ice delivery to the terminus and the erosion rate, providing quantitative insight into the relationship between ice fluxes and glacial erosion rates. For Glaciar Marinelli, as for other calving glaciers for which suitable data exist, the marked retreat and thinning over the past 50 years have resulted in a period of accelerated basal sliding and unusually rapid erosion.

A fresh-snow sampling campaign was conducted during the late austral summer of 2006 in the accumulation zone of Glaciar Marinelli, located in the Cordillera Darwin, Tierra del Fuego, Chile. Snow samples were analyzed for stable isotopes (δ18O), major soluble ions (Na+, K+, Ca2+, Mg2+, Cl−, NO3−, SO42−, MS−) and major and trace elements (Na, Mg, Al, Fe, Ca, Sr, Cd, Cs, Ba, La, Ce, Pr, Dy, Ho, Er, Bi, U, As, Ti, V, Cr, Mn). The dominance of marine chemistry resembles that in studies from Patagonian glaciers. Snow chemistry was dominantly loaded by marine species (Cl−, Na+ and ssSO42−), while contributions of crustal species (e.g. Al and Fe) were very low. Empirical orthogonal function analysis suggests two possible dust sources, one represented by Al and Fe and the other by La, Ce and Pr. Enrichment-factor calculations suggest the majority of elements are within average upper-crustal ratios, but major enrichments of Bi and Cd (hundreds of times) suggest possible anthropogenic sources. Linear correlation of δ18O and barometric pressure (r = 0.60, p < 0.007) suggests a potential ‘amount effect’ relationship between depleted δ18O ratios and stronger storm conditions (e.g. greater precipitation). The snow-chemistry records from Glaciar Marinelli are the first measured in Tierra del Fuego, the southernmost glaciated region outside Antarctica.

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.

A linear stability analysis of a two-dimensional flow of an isothermal ice sheet interacting with the ocean is considered. The set of boundary conditions determining motion of the grounding line is adopted to describe hydrostatic equilibrium of ice in water and a cubic dependence of the mass flow rate on ice thickness. The numerical analysis shows that the zero-growth (zero-eigenvalue) mode found for linear bed slopes and constant accumulation rates indeed determines neutral equilibrium and separates stable and unstable solutions. It is also argued that, provided some conditions of regularity of the solutions are satisfied, finding only one stable and one unstable solution would be enough to ascertain that the condition determining a zero eigenvalue also determines neutral equilibrium. This supports the intuitive understanding of ice-sheet stability: ice sheets are stable on bed slopes that ensure that the mass flow rate at the grounding line increases faster than the cumulative ice accumulation rate at the surface when the grounding line is perturbed; and ice sheets are unstable otherwise.

Many glacier snouts in the Himalaya are known to be stagnant and exhibiting low surface gradients, conditions that are conducive to the formation of glacial lakes impounded either by the terminal moraine or by the remnant glacier snout. In this study, we use interferometry and feature-tracking techniques to quantify the extent of stagnation in 20 glaciers across the Everest (Qomolangma; Sagarmatha) region, and subsequently we examine the relationship between local catchment topography and ice dynamics. The results show that only one of the studied glaciers, Kangshung Glacier, is dynamic across its entire surface, with flow rates greater than 40 m a−1 being recorded in high-elevation areas. Twelve other glaciers show some evidence of flow, but are generally characterized by long, stagnant tongues, indicating widespread recession and in situ decay. The remaining seven glaciers show no evidence of flow in any of the available datasets. Hypsometric data suggest that catchment topography plays an important role in controlling glacier flow regimes, with those fed by wide, high-altitude accumulation areas showing the most extensive active ice, and those originating at low elevations exhibiting large areas of stagnant ice. Surface profiles extracted from a SRTM digital elevation model indicate that stagnant snouts are characterized by very low (<2°) surface angles and that down-wasting is the prevalent ablation pattern in the study area.

High-resolution airborne lidar data are used to produce digital elevation models (DEMs) of an arctic valley glacier (midre Lovénbreen, Svalbard) at resolutions of 2.5–2000 m, using three different interpolation schemes. These data are used in a distributed model of solar radiation loading for glaciers. When the mean of all lidar measurements within a DEM cell is used to calculate cell height, the differences between the finest- (2.5 m) and coarsest-resolution (2000 m) DEMs for the calculated annual whole-glacier spatial means of total potential direct-beam solar radiation, potential duration of direct-beam solar radiation, and intensity of potential direct-beam solar radiation are 20%, 56% and −23% of the 2.5 m DEM values respectively. A resolution change from 2.5 m to 200 m affects the whole-glacier spatial mean summer net solar radiative flux by an average of 5%, and the summer melt production from the glacier by an average of 3% compared with the 2.5 m DEM values, for the years 2001–03. These changes are largely driven by underestimation of shading by surrounding topography at coarser DEM resolutions. This dependency is reduced in the second and third interpolation schemes, especially at resolutions finer than 50 m, which use the maximum lidar height measurement in some or all DEM cells. These results suggest that resolutions of ∼50 m are the coarsest that should be adopted in high-resolution glacier surface energy-balance models for glaciers of similar size and in similar topographic situations to midre Lovénbreen, and that the impact of DEM resolution on calculated solar radiation receipts can be reduced by an appropriate choice of DEM interpolation scheme.

We provide the first direct evidence that a number of water-soluble compounds, in particular calcium sulfate (CaSO4·2H2O) and calcium carbonate (CaCO3), are present as solid, micron-sized inclusions within the Greenland GRIP ice core. The compounds are detected by two independent methods: micro-Raman spectroscopy of a solid ice sample, and energy-dispersive X-ray spectroscopy of individual inclusions remaining after sublimation. CaSO4·2H2O is found in abundance throughout the Holocene and the last glacial period, while CaCO3 exists mainly in the glacial period ice. We also present size and spatial distributions of the micro-inclusions. These results suggest that water-soluble aerosols in the GRIP ice core are dependable proxies for past atmospheric conditions.

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.

Glaciers are widely recognized as key indicators of climate change, and their meltwater plays an important role in hydropower production in Norway. Since the last glacier inventory was compiled in northern Norway in the 1970s, marked fluctuations in glacier length and mass balance have been reported for individual glaciers, and the current overall glacier state is thus not well known. Within the framework of the Global Land Ice Measurements from Space (GLIMS) initiative, we have created a new inventory for 489 glaciers in the Svartisen region, northern Norway, using a Landsat Enhanced Thematic Mapper Plus (ETM+) satellite scene from 7 September 1999 and automated multispectral glacier mapping (thresholded band ratios). In addition, visual inspection and correction of the generated glacier outlines has been applied. Adverse snow conditions and uncertain drainage divides made glacier mapping challenging in some regions of the study site. Glacier outlines from 1968, as digitized from a topographic map, were used for a quantitative change assessment for a selection of 300 glaciers. The overall area change of this sample from 1968 to 1999 was close to zero, but with a strongly increasing scatter towards smaller glaciers, large area gains (>50%) for small glaciers (<1 km2), and an unexpected stronger relative area loss towards the wetter coast. The overall size changes are small (<1%) for the three largest ice masses in the study region (Vestisen, Østisen and Blåmannsisen).

The glaciers on Yulong mountain (5596 m a.s.l.), China, are the southernmost in mainland Eurasia. The largest is Baishui glacier No. 1. Up to 90% of the annual precipitation there falls in summer, when warm moisture-rich air masses associated with the southwest monsoon reach the area from the Indian Ocean. The winter climate is influenced by air masses with a continental origin and by the southern branch of the westerlies. The snow that accumulates on Baishui glacier No. 1 includes marine aerosols associated with the summer monsoon, and dust brought in winter from central/west Asia, Africa or the Thar Desert area. Studies in May 2006 at two sites, one in the accumulation area (4900 m) and one in the ablation area (4750 m), revealed differences between the ionic composition of the snow that had accumulated in the 2005/06 winter and that of the snow which had been deposited during the preceding summer monsoon. Differences in the chemistry of the summer-accumulated snow at the two sites probably reflected local differences in ablation and elution rates, rather than differences in ion supplies. Differences in the chemistry of the winter-accumulated snow may reflect the influence of up-valley winds, which bring more crustal material to the lower site.

From 1999 to 2001 a 724 m deep ice core was drilled on Akademii Nauk ice cap, Severnaya Zemlya, to gain high-resolution proxy data from the central Russian Arctic. Despite strong summertime meltwater percolation, this ice core provides valuable information on the regional climate and environmental history. We present data of stable water isotopes, melt-layer content and major ions from the uppermost 57 m of this core, covering the period 1883–1998. Dating was achieved by counting seasonal isotopic cycles and using reference horizons. Multi-annual δ18O values reflect Eurasian sub-Arctic and Arctic surface air-temperature variations. We found strong correlations to instrumental temperature data from some stations (e.g. r = 0.62 for Vardø, northern Norway). The δ18O values show pronounced 20th-century temperature changes, with a strong rise about 1920 and the absolute temperature maximum in the 1930s. A recent decrease in the deuterium-excess time series indicates an increasing role of the Kara Sea as a regional moisture source. From the multi-annual ion variations we deduced decreasing sea-salt aerosol trends in the 20th century, as reflected by sodium and chloride, whereas sulphate and nitrate are strongly affected by anthropogenic pollution.

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.

We use an empirically validated high-resolution three-dimensional ice-sheet model to investigate the mass-balance regime, flow mechanisms and subglacial characteristics of a simulated Younger Dryas Stadial ice cap in Scotland, and compare the resulting model forecasts with geological evidence. Input data for the model are basal topography, a temperature forcing derived from GRIP δ18O fluctuations and a precipitation distribution interpolated from modern data. The model employs a positive-degree-day scheme to calculate net mass balance within a domain of 112500 km2, which, under the imposed climate, gives rise to an elongate ice cap along the axis of the western Scottish Highlands. At its maximum, the ice cap is dynamically and thermally zoned, reflecting topographic and climatic controls, respectively. In order to link these palaeoglaciological conditions to geological interpretations, we calculate the relative balance between sliding and creep within the simulated ice cap, forecast areas of the ice cap with the greatest capacity for basal erosion and predict the likely pattern of subglacial drainage. We conclude that ice flow in central areas of the ice cap is largely due to internal deformation, and is associated with geological evidence of landscape preservation. Conversely, the distribution of streamlined landforms is linked to faster-flowing ice whose velocity is predominantly the result of basal sliding. The geometry of the main ice mass focuses subglacial erosion in the mid-sections of topographic troughs, and produces glaciohydraulic gradients that favour subglacial drainage through low-order arterial routes.