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.

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.

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.

A recently developed dating method for glacier ice, based on the analysis of radiocarbon in carbonaceous aerosol particles, is thoroughly investigated. We discuss the potential of this method to achieve a reliable dating using examples from a mid- and a low-latitude ice core. Two series of samples from Colle Gnifetti (4450 m a.s.l., Swiss Alps) and Nevado Illimani (6300 m a.s.l., Bolivian Andes) demonstrate that the 14C ages deduced from the water-insoluble organic carbon fraction represent the age of the ice. Sample sizes ranged between 7 and 100 μg carbon. For validation we compare our results with those from independent dating. This new method is thought to have major implications for dating non-polar ice cores in the future, as it provides complementary age information for time periods not accessible with common dating techniques.

Dry-snow slab avalanches initiate from a failure in a weak snow layer below a cohesive slab. Snow is considered as a porous ice structure, and the strength distribution of the single elements of this structure, i.e. grains and bonds between grains, shows a high degree of disorder. On the bond or microstructural level, the failure process is believed to start if the fracturing of bonds between snow grains is not balanced by the formation of new bonds. We use a statistical fracture model – a fibre bundle model – to study the failure process in a weak snow layer. The model consists of fibres of various strengths representing single snow grains between two rigid plates which represent the slab above and the substratum below the weak layer. The fibres deform in a linear elastic manner and break instantly at their rupture strength. Broken fibres may sinter (re-bond) and regain strength after a finite sintering time. We show that the different characteristic times for breaking and sintering lead to the rate dependence of snow strength. This is, to our knowledge, the first statistical model to reproduce the ductile-to-brittle transition which snow exhibits with increasing strain rate. When the model is applied to simulate experimental stress–strain curves for different strain rates, the model and experimental results are in fair agreement.

Historically, dry-snow metamorphism has been classified by the thermal environment and thermodynamic processes in a snowpack. Snow experiencing predominantly macroscopic isothermal conditions develops different microstructure than snow subjected to large temperature gradients. As such, much previous research has been categorized by and limited to specific thermal conditions. The current research expands a generalized approach for the movement of heat and mass to include a snow crystal kinetic growth model. An existing spiral defect propagation theory for kinetic growth on simple faceted geometry is utilized. Primary crystal habit as a function of temperature is incorporated. A model of heat and mass transfer through an ice and pore structure is coupled with phase-change thermodynamics during kinetic growth. A kinetic growth microstructure model is developed and integrated into heat and mass transfer representations, which are solved using finite-difference techniques. The kinetic morphology model approximates frequently observed hopper-type crystals. The snow microstructure is allowed to change at every step, resulting in a transient description of kinetic growth metamorphism. Variable kinetic growth rates are demonstrated based on temperature and on crystallographic orientation relative to a temperature gradient. Crystals preferentially aligned with the temperature gradient have significantly higher growth rates, supporting previous observations of predominant crystal habits developing under temperature gradient conditions. Grain-size dispersion increase with time is demonstrated and supported experimentally in the literature. A dominant grain growth theory based on crystallographic orientation that has been previously postulated is supported. A broad range of metamorphic geometric parameters and thermal conditions may now be simulated with a single model.

We investigate the evolution of snow temperature, water content, density and stable water isotopes of δ18O at four Arctic snow-pit sites during early-season melt, in order to understand the effects of melt on snowpack stratigraphies and seasonal isotopic signals. We relate isotopic changes observed at these sites to temperature reconstructions derived from a 33 year firn-core record drilled on the same icefield. Decreases in seasonal isotopic amplitudes observed at all but one snow-pit site coincide with the percolation of more enriched meltwater into the snowpack, suggesting that meltwater percolation is the dominant process causing isotopic redistribution in Arctic snowpacks during the melt season. The decrease in isotopic range was accompanied by increases in mean δ18O values at all snow-pit sites. Positive degree-day (PDD) calculations are used to relate the amount of melt observed at the low-elevation snow-pit sites to the firn-core site. Results based on PDD values suggest an average overestimation of 1.1°C in average annual temperature reconstructions from the firn-core site from 1967 to 2006, with the possibility of errors in excess of 3°C during high-melt years.

Basal radar reflectivity is the most important measurement for the detection of subglacial water. However, dielectric loss in the overlying ice column complicates the determination of basal reflectivity. Dielectric attenuation is a function of ice temperature and impurity concentration. Temperature distribution is a function of climate history, basal heat flow and vertical strain rate, all of which can be partially inferred from the structure of dated internal layers. Using 11 dated layers, isotope records from Dome C, East Antarctica, and a model of the spatial variation of geothermal flux, we calculate the vertical strain rate and accumulation-rate history, allowing identification of areas where the basal melt rate exceeds 1.5 mm a−1. The accumulation-rate history and vertical strain rates are then used as inputs for a transient temperature model. The model outputs for the present-day temperature distribution are then combined with depth-dependent ionic concentrations to model dielectric loss and infer basal reflectivity. The resulting reflection coefficients are consistent (∼−5 dB) across a variety of subglacial water bodies. We also identify a high reflectivity >−15 dB in Concordia Trench and along suspected subglacial water-flow routes in Vincennes Basin. Highland areas tend to have highly variable reflection coefficients near −30 dB, consistent with an ice–bedrock interface. This combined model also identifies three areas of enhanced basal melting along Concordia Ridge, Concordia Subglacial Lake and Vincennes Basin. Melt at Concordia Subglacial Lake exceeds 5 mm a−1. The inferred basal melt at these locations is not possible without enhanced geothermal flux. We demonstrate how radar-sounding data can provide both input and verification for a self-consistent model of vertical strain, vertical temperature distribution and meltwater distribution.

In January 2000 nine snow sensors were deployed in Dronning Maud Land, Antarctica, along a 355 km transect from Kvitkuven near the shelf edge via the Finnish research station, Aboa, to the Amundsenisen plateau. The purpose was to test a sensor system for spatio-temporal variations in temperature across the snow–air interface and snow accumulation/ablation, which includes atmospheric net balance and migrating snow dunes. In the dry snow conditions, environmental static electricity interfered with data transfer; several sensors were disabled early, while the longest record reached 6 months. Along the main transect, the year 2000 mass balance ranged from 52 to 221 mm w.e., largely following spatial patterns seen by other researchers. The level increased toward the edge of the ice sheet; unloading of drifting snow as the slope flattens, rather than increased snowfall as previously thought, may be responsible for this. At the Högisen ice dome site the mass balance was 897 mm w.e., possibly due to unloading of wind-blown snow, as katabatic winds are locally forced uphill. This mechanism is important to maintain such topographic features along the Antarctic ice sheet margin. Major precipitation events occurred at 3–5 week intervals and much of the precipitation fell before mid-June. The daily signal in temperature disappeared after the autumn equinox, then the spectrum displayed a broad peak at synoptic frequencies. Potential temperature decreased towards the shelf edge, displaying a pattern consistent with strong inversions and suggesting that strong evaporative cooling is associated with the katabatic winds.

Taku Glacier, Alaska, USA, is currently in the advance stage of the tidewater glacier cycle. We investigated the near-terminus dynamics by measuring surface velocities, surface elevation changes, ice thickness and ablation. Velocities vary on sub-daily, diurnal, seasonal and interannual timescales. Flowline modeling shows that the modeled surface velocities are sensitive to changes in till yield strength and thus effective basal pressures. The glacier bed deepens in the up-glacier direction and this imposes a minimum subglacial water pressure necessary for water to drain along the bed. In a simple model we impose water-pressure gradients based on phreatic surfaces of constant slopes to simulate the winter–summer transitions. This proves sufficient to explain an observed early-season switch from compressional to block flow. Velocities also vary between years. Changing basal conditions can result in lower horizontal velocities, which decrease the ice supply to the terminus and result in temporary surface lowering. But a decrease in ice flux to the terminus must lead to ice storage further upstream, and that ice mass will eventually reach the terminus. This can explain the observed episodic nature of terminus advance.

The recent dynamic behaviour of Karakoram glaciers is expected to differ from that shown by glaciers in the central and eastern Himalaya because of regional variations in precipitation and temperature trends. However, there are insufficient quantitative data to support or confute such hypotheses. We present velocity data covering the period 1993–2008 for Baltoro Glacier, one of the longest glaciers in the Karakoram. Velocity measurements were made using cross-correlation feature tracking applied to European Remote-sensing Satellite (ERS-1 and -2) and Envisat advanced synthetic aperture radar (ASAR) data, supplemented by differential global positioning system (DGPS) measurements. We find a gradual acceleration of the glacier during the early 2000s, in particular during winter months. Multi-seasonal data reveal a large difference between summer and winter flow characteristics, but only in the upper ablation zone. Summer 2005 was a particularly dynamic period following from the heavy winter snowfall of 2004, indicating the importance of basal meltwater availability for glacier flow. Transverse velocity profiles indicate that Baltoro Glacier undergoes ‘block’ flow across much of the upper ablation zone during the summer, which we interpret as evidence of widespread basal sliding. The DGPS data confirm the rapid increase in flow detected during 2005. Modelled climatic data reveal decreasing summer temperatures and increasing precipitation over the study period, helping to explain the observed dynamic variations and their differences from glaciers elsewhere in the Himalaya.

Increased mass losses from the Greenland ice sheet and inferred contributions to sea-level rise have heightened the need for hydrologic observations of meltwater exiting the ice sheet. We explore whether temporal variations in ice-sheet surface hydrology can be linked to the development of a downstream sediment plume in Kangerlussuaq Fjord by comparing: (1) plume area and suspended sediment concentration from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery and field data; (2) ice-sheet melt extent from Special Sensor Microwave/Imager (SSM/I) passive microwave data; and (3) supraglacial lake drainage events from MODIS. Results confirm that the origin of the sediment plume is meltwater release from the ice sheet. Interannual variations in plume area reflect interannual variations in surface melting. Plumes appear almost immediately with seasonal surface-melt onset, provided the estuary is free of landfast sea ice. A seasonal hysteresis between melt extent and plume area suggests late-season exhaustion in sediment supply. Analysis of plume sensitivity to supraglacial events is less conclusive, with 69% of melt pulses and 38% of lake drainage events triggering an increase in plume area. We conclude that remote sensing of sediment plume behavior offers a novel tool for detecting the presence, timing and interannual variability of meltwater release from the ice sheet.

Low-elevation blue-ice areas in Antarctica can contain significant amounts of subsurface liquid water during summer and may experience internal melt–freeze cycles due to absorption of shortwave radiation. An existing 2-D (x-z) model has been used to study the phenomenon and this is improved by changing the lower boundary condition from a no-flux condition to one that lets heat through, and by changing the bulk optical properties to spectral values. Both changes made the model more realistic. The optical and thermal boundary conditions, together with the optical attenuation coefficient, have a large effect on the amount of water produced in the ice. Our results show that if the lower boundary condition is changed from no flux to radiating, subsurface melting is reduced dramatically. Using a spectral albedo produces less melting than using a corresponding bulk albedo, the other variables left unchanged. If optical properties are changed to spectral values, the melting is different than using bulk values. If the linear extent of the blue-ice area is <0.6 m, subsurface melting is diminished. We found the spatial variability of snow has a significant effect on subsurface melting. Subsurface melting was found to be impossible for albedos >0.7 and the subsurface ponds can persist over the winter if the albedo is <0.4.

A gamma-ray density gauge can provide high-resolution and high-precision density measurements of firn and ice cores. This study describes the design, gamma-ray energy optimization and mass attenuation coefficient calibration of the Maine Automated Density Gauge Experiment (MADGE), a portable, field-operable gamma-ray density gauge used on overland traverses in East Antarctica. The MADGE instrument uses a 241Am gamma-ray source, a pulse-mode counting system and electronic core diameter calipers to collect high-precision (±0.004 g cm−3) density data from 3–8 cm diameter firn and ice cores. The data are collected at a 3.3 mm spatial resolution and an average throughput of 1.5 m h−1 for 5 cm diameter cores.

We investigate the dynamics of a small surge-type valley glacier as part of a study to characterize glacier response to climate in the Donjek Range, southwest Yukon, Canada. Pole displacements were measured using kinematic GPS techniques during three consecutive summer field seasons. Measured surface velocities range from <10 m a−1 over the lower 1500 m of the 5 km long glacier to a maximum of ∼25–35 m a−1 over the upper 3500 m. Basal velocities along an approximate flowline are reconstructed from the measured surface velocities using inverse methods. Control tests are used to validate the inversion scheme, and sensitivity tests are performed to evaluate the influence of the flow-law coefficient, shape factor and longitudinal averaging length. Inversion of the real data shows that basal motion accounts for 50–100% of the total surface motion along the flowline. Based on these results, and several other lines of evidence, we suggest this glacier may be undergoing a slow surge.