The analysis of historical avalanche data is important when developing accurate hazard maps. The record of snow-avalanche disasters on Sakhalin and the Kuril Islands is incomplete, due to the historical division into periods of Japanese and Russian rule. Here we combine and analyze data from Japanese and Russian sources to reconstruct a continuous record of avalanche catastrophes in the region from 1910 to 2010. Despite the relatively small scale of the majority of catastrophic avalanches, with a total vertical drop < 200 m, we document evidence that places the region among the most avalanche-affected areas in the world. In total, 756 fatalities and > 238 injuries have occurred in 275 incidents over a 100 year period (two-thirds of those killed were Japanese). This death toll is higher than that in Canada, New Zealand or Iceland, or non-recreational fatalities in France. A wave of avalanche disasters (1930s–60s) following intense colonization of Sakhalin and the Kuril Islands is evident. Although this ‘wave’ could be considered a local issue of the past, many presently developing countries may face similar situations. The fatality rate has decreased over time, due to social factors, and differs from that of any other region, in its absence of deaths through recreational activities. Although in recent years the fatality rate is lower than that of Iceland or the USA, the per capita avalanche casualty rate on Sakhalin and the Kuril Islands remains among the highest in the world.

Shadowgraphy was employed to study snow saltation in boundary-layer wind tunnel experiments with fresh, naturally deposited snow. The shadowgraphy method allowed for a temporally and spatially high-resolution investigation of snow particle characteristics within a measurement area of up to 50 mm × 50 mm. Snow particle size and number characteristics, and their variation with height in the saltation layer, were analysed. The following observations and findings were made for the saltation layer: (1) the particle number decreases exponentially with height, (2) the mean particle diameter is fairly constant, with a very slight tendency to decrease with height, (3) the maximum particle diameter decreases linearly with height, and (4) the snow particle size distribution can be adequately described by gamma probability density functions. The shape and scale parameters of the gamma distribution were found to vary systematically, though only slightly, with height over ground and between experiments with different snowpack characteristics.

This paper provides an introduction to adjoint methods, which are used to find the gradient of an objective function, as required by optimization algorithms. Examples are included, culminating in a data-assimilation problem from glaciology.

Compression experiments on laboratory-grown columnar S2 saltwater ice loaded triaxially through proportional loading at T = –20°C at applied strain rates of demonstrate that plastic (P) faulting is a mode of failure in saltwater ice when rapidly loaded under a high degree of confinement. In terms of microstructure, mechanical behavior and strength, saltwater ice that fails via P-faulting is almost indistinguishable from columnar S2 freshwater ice that fails via P-faulting loaded under the same conditions. The results also demonstrate that saltwater ice loaded rapidly may exhibit yet another mode of failure, in addition to P-faulting, through what appears to be a mechanism of pore collapse.

Rock glaciers are periglacial alpine landforms that are found in many locations worldwide. Whereas well-developed models of deformation are established for traditional alpine glaciers, rock glacier deformation is poorly understood. Geophysical data from Lone Peak Rock Glacier (LPRG), southwest Montana, USA, are paired with lidar bare-earth 1 m digital elevation model (DEM) analysis to explore potential genetic relationships between internal composition, structure and regularly spaced arcuate transverse ridges expressed at the rock glacier surface. The internal composition of LPRG is heterogeneous, with frozen debris and clean ice overlain by an unconsolidated talus mantle. Upslope-dipping, clearly distinguished reflectors in the ground-penetrating radar (GPR) longitudinal survey at LPRG correspond to transverse ridges. The spacing and slope of individual features at the surface and in the subsurface were measured and compared and are found to be similar. The structures observed at LPRG and other rock glaciers are similar to structures detected in glaciotectonically altered sediment, ice-cored moraines and other rock glacier settings. This finding suggests that transverse ridges on rock glaciers may be used as geomorphic indicators of internal deformation. This study contributes to the body of research on the application of GPR to rock glaciers, and is the first to directly pair and analyze individual surface topographic features with internal structures.

In the deep ice cores drilled at the GRIP, NGRIP and GISP2 sites in Greenland and at Byrd Station and the summit of Law Dome in Antarctica, the mean crystal size increases with depth in the shallow subsurface and reaches steady values at intermediate depth. This behaviour has been attributed to the competition between grain-boundary migration driven crystal growth and crystal polygonization, but the effects of changing crystal dislocation density and non-equiaxed crystal shape in this competition are uncertain. We study these effects with a simple model. It describes how the mean height and width of crystals evolve as they flatten under vertical compression, and as crystal growth and polygonization compete. The polygonization rate is assumed to be proportional to the mean dislocation density across crystals. Migration recrystallization, which can affect crystal growth via strain-induced grain boundary migration but whose impact on the mean crystal size is difficult to quantify for ice at present, is not accounted for. When applied to the five ice-core sites, the model simulates the observed crystal-size profiles well down to the bottom of their steady regions, although the match for Law Dome is less satisfactory. Polygonization rate factors retrieved for the sites range from 10–5 to 10–2 a–1. We conclude that since crystal size and dislocation density evolve in a strongly coupled manner, consistent modelling requires multiple differential equations to track both of these variables. Future ice-core analysis should also determine crystal size in all three principal directions.

The Wind River Range in Wyoming, USA, contains the largest concentration of glacial mass in the Rocky Mountains of the contiguous USA. Despite this distinction, only a few field or remotely sensed studies providing glacier volume changes have been published. The current study focuses on Continental Glacier located on the northern end of the range and uses two field datasets (high-accuracy GPS surface elevation points and ice-penetrating radar transects of the glacier bed) to create a three-dimensional model of glacier volume. Current surface elevations are compared with historical elevation data to calculate surface elevation change over time. An average thinning rate of 13.8 ± 7.8 m (0.30 ± 0.17 m a–1) between 1966 and 2012 was found. Surface elevation change rates varied across the glacier, ranging from +0.30 to –0.98 m a–1. Taking into account variable melt rates across the glacier, along with a glacial volume of 72.1 × 106 ± 10.8 × 106 m3, we estimate that Continental Glacier will be reduced in volume by 43% over the next 100 years and will disappear completely over the next 300–400 years, if current climatic conditions persist.

Petermann Gletscher, northwest Greenland, drains 4% of the Greenland ice sheet into Nares Strait. Its floating ice shelf retreated from 81 to 48 km in length during two large calving events in 2010 and 2012. We document changes in the three-dimensional ice-shelf structure from 2000 to 2012, using repeated tracks of airborne laser altimetry and ice radio-echo sounding, ICESat laser altimetry and MODIS visible imagery. The recent ice-shelf velocity, measured by tracking surface features between flights in 2010 and 2011, is ~1.25 km a−1, ~15–30% faster than estimates made before 2010. The steady- state along-flow ice divergence represents 6.3 Gta−1 mass loss through basal melting (~5Gta−1) and surface melting and sublimation (~1.0Gta−1). Airborne laser altimeter data reveal thinning, both along a thin central channel and on the thicker ambient ice shelf. From 2007 to 2010 the ice shelf thinned by ~5 m a−1, which represents a non-steady mass loss of ~4.1 Gta−1. We suggest that thinning in the basal channels structurally weakened the ice shelf and may have played a role in the recent calving events.

Three stations near the calving front of the Ross Ice Shelf, Antarctica, recorded GPS data through a full spring–neap tidal cycle in November 2005. The data revealed a diurnal horizontal motion that varied both along and transverse to the long-term average velocity direction, similar to tidal signals observed in other ice shelves and ice streams. Based on its periodicity, it was hypothesized that the signal represents a flow response of the Ross Ice Shelf to the diurnal tides of the Ross Sea. To assess the influence of the tide on the ice-shelf motion, two hypotheses were developed. The first addressed the direct response of the ice shelf to tidal forcing, such as forces due to sea-surface slopes or forces due to sub-ice-shelf currents. The second involved the indirect response of ice-shelf flow to the tidal signals observed in the ice streams that source the ice shelf. A finite-element model, based on viscous creep flow, was developed to test these hypotheses, but succeeded only in falsifying both hypotheses, i.e. showing that direct tidal effects produce too small a response, and indirect tidal effects produce a response that is not smooth in time. This nullification suggests that a combination of viscous and elastic deformation is required to explain the observations.

Snow stratigraphy and water percolation are key contributing factors to avalanche formation. So far, only destructive methods can provide this kind of information. Radar technology allows continuous, non-destructive scanning of the snowpack so that the temporal evolution of internal properties can be followed. We installed an upward-looking ground-penetrating radar system (upGPR) at the Weissfluhjoch study site (Davos, Switzerland). During two winter seasons (2010/11 and 2011/12) we recorded data with the aim of quantitatively determining snowpack properties and their temporal evolution. We automatically derived the snow height with an accuracy of about ± 5 cm, tracked the settlement of internal layers (± 7 cm) and measured the amount of new snow (± 10 cm). Using external snow height measurements, we determined the bulk density with a mean error of 4.3% compared to manual measurements. Radar-derived snow water equivalent deviated from manual measurements by 5%. Furthermore, we tracked the location of the dry-to-wet transition in the snowpack until water percolated to the ground. Based on the transition and an independent snow height measurement it was possible to estimate the volumetric liquid water content and its temporal evolution. Even though we need additional information to derive some of the snow properties, our results show that it is possible to quantitatively derive snow properties with upGPR.

Ice-thickness measurements in Antarctic ice-shelf grounding zones are necessary for calculating the mass balance of individual catchments, but remain poorly constrained for most of the continent. We describe a new inverse modelling optimization approach to estimate ice thickness in the grounding zone of Antarctic outlet glaciers and ice shelves using spatial patterns of tide-induced flexure derived from differential interferometric synthetic aperture radar (InSAR). We demonstrate that the illposedness of the inverse formulation of the elastic-plate equations for bending can be controlled by regularization. In one dimension, the model recreates smooth, synthesized profiles of ice thickness from flexure information to within 1–2%. We test the method in two dimensions and validate it in the grounding zone of Beardmore Glacier, a major outlet glacier in the Transantarctic Mountains, using interferograms created from TerraSAR-X satellite imagery acquired in 2012. We compare our results with historic and modern ice-thickness data (radio-echo sounding from 1967 and ground-penetrating radar from 2010). We match both longitudinal and transverse thickness transects to within 50 m root-mean-square error using an effective Young’s modulus of 1.4 GPa. The highest accuracy is achieved close to the grounded ice boundary, where current estimates of thickness based on surface elevation measurements contain a systematic bias towards thicker ice.

The Randolph Glacier Inventory (RGI) is a globally complete collection of digital outlines of glaciers, excluding the ice sheets, developed to meet the needs of the Fifth Assessment of the Intergovernmental Panel on Climate Change for estimates of past and future mass balance. The RGI was created with limited resources in a short period. Priority was given to completeness of coverage, but a limited, uniform set of attributes is attached to each of the ~198 000 glaciers in its latest version, 3.2. Satellite imagery from 1999–2010 provided most of the outlines. Their total extent is estimated as 726 800 ± 34 000 km2. The uncertainty, about ±5%, is derived from careful single-glacier and basin-scale uncertainty estimates and comparisons with inventories that were not sources for the RGI. The main contributors to uncertainty are probably misinterpretation of seasonal snow cover and debris cover. These errors appear not to be normally distributed, and quantifying them reliably is an unsolved problem. Combined with digital elevation models, the RGI glacier outlines yield hypsometries that can be combined with atmospheric data or model outputs for analysis of the impacts of climatic change on glaciers. The RGI has already proved its value in the generation of significantly improved aggregate estimates of glacier mass changes and total volume, and thus actual and potential contributions to sea-level rise.

The potential for future dynamical instability of Pine Island Glacier, West Antarctica, has been addressed in a number of studies, but information on its past remains limited. In this study we use airborne radio-echo sounding (RES) data acquired over Pine Island Glacier to investigate past variations in accumulation pattern. In the dataset a distinctive pattern of layers was identified in the central part of the glacier basin. We use these layers as chronological identifiers in order to construct elevation maps of the internal stratigraphy. The observed internal layer stratigraphy is then compared to calculated stratigraphy from a three-dimensional ice-flow model that has been forced with different accumulation scenarios. The model results indicate that the accumulation pattern is likely to have changed at least twice since the deposition of the deepest identified layer. Additional RES data linked to the Byrd ice core provide an approximate timescale. This timescale suggests that the layers were deposited at the beginning of or during the Holocene period. Thus the widespread changes occurring in the coastal extent of the West Antarctic ice sheet at the end of the last glacial period could have been accompanied by changes in accumulation pattern.

During summer 2011, seismic activity in the ablation zone of the western Greenland ice sheet (GrIS) was monitored using a network of three-component seismometers. The seismic record includes a large variety of icequakes and seismic tremors that demonstrate a clear correlation with subglacial water flow. We verified the existence of well-known shallow icequakes (related to surface crevasse formation), deep icequakes (located at 100–160 m depth) and narrow-banded short-term seismic tremors (tens of seconds in duration). In addition, we present previously unreported long-term tremors lasting several hours. Using attenuation of the measured tremor amplitude, we locate the epicentre of this long-term tremor to a large moulin within our study area. Between 3 and 11 Hz, our continuous seismic record is dominated by this ‘moulin tremor’ and shows strong correlation with the water level of the generating moulin. We argue that monitoring of icequake and glacial tremor sources bears high potential for investigating glacier hydraulics and dynamics, and is thus an ideal supplement to traditional glaciological measurements.

Marine ice is an important component of ice shelves in Antarctica. It accretes in substantial amounts at weak points and below ice shelves. It is likely to exhibit peculiar rheological properties, which are crucial to understanding its potential role in stabilizing ice-shelf flow. Due to its location and consolidation processes, marine ice can present a variety of textures which are likely to influence its rheological properties. We present a new dataset of unconfined uniaxial compression experiments on folded marine ice samples that have been cut at various angles to the folds. Texture and fabric analyses are described ‘before’ and ‘after’ the deformation experiment. It is shown that, in the given stress configuration, the geometry of the anisotropy controls the rheological behaviour of the marine ice. During secondary creep, folded marine ice is harder to deform than weakly textured ice when compressed parallel or perpendicular to the folds' hinges, while the reverse is true for ice compressed at 45°. The observed range of values for the n exponent in Glen's flow law is between 2.1 and 4.1. Surprisingly, we see that tertiary creep tends to develop at a higher total strain than for randomly oriented impurity-free meteoric ice.

Distributed energy-balance melt models have rarely been applied to glaciers with extensive supraglacial debris cover. This paper describes the development of a distributed melt model and its application to the debris-covered Miage glacier, western Italian Alps, over two summer seasons. Sub-debris melt rates are calculated using an existing debris energy-balance model (DEB-Model), and melt rates for clean ice, snow and partially debris-covered ice are calculated using standard energy-balance equations. Simulated sub-debris melt rates compare well to ablation stake observations. Melt rates are highest, and most sensitive to air temperature, on areas of dirty, crevassed ice on the middle glacier. Here melt rates are highly spatially variable because the debris thickness and surface type varies markedly. Melt rates are lowest, and least sensitive to air temperature, beneath the thickest debris on the lower glacier. Debris delays and attenuates the melt signal compared to clean ice, with peak melt occurring later in the day with increasing debris thickness. The continuously debris-covered zone consistently provides 30% of total melt throughout the ablation season, with the proportion increasing during cold weather. Sensitivity experiments show that an increase in debris thickness of 0.035 m would offset 18C of atmospheric warming.