The deadliest outburst flood from an englacial cavity occurred on Glacier de Tête Rousse in the Mont Blanc area, French Alps, in 1892. A subglacial reservoir was discovered in the same glacier in 2010 and drained artificially in 2010, 2011 and 2012 to protect the 3000 inhabitants downstream. The mechanism leading to the spontaneous refilling of the cavity following these pumping operations has been analyzed. For this purpose, the subglacial water volume changes between 2010 and 2013 were reconstructed. The size of the cavity following the pumping was found to have decreased from 53500 m3 in 2010 to 12 750 m3 in 2013. Creep and the partial collapse of the cavity roof explain a large part of the volume loss. Analysis of cavity filling showed a strong relationship between measured surface melting and the filling rate, with a time delay of 4–6 hours. A permanent input of 15 m3 d−1, not depending on surface melt, was also found. The meltwater and rain from the surface is conveyed to bedrock through crevasses and probably through a permeable layer of rock debris at the glacier bed. The drainage pathway permeability was estimated at 0.054 ms−1 from water discharge measurements and dye-tracing experiments.

New high-resolution velocity maps of the eastern St Elias Mountains, North America, are obtained from speckle tracking of winter 2011 and 2012 RADARSAT-2 image pairs. This includes the most complete velocity mapping to date of Hubbard Glacier, allowing for an upward revision of the Hubbard Glacier calving flux to 5.48 ± 1.16 km3 a−1. Combined with historical velocities from feature tracking of Landsat image pairs (1980s−2000s), and previously published results, these new velocity measurements allow for an evaluation of the interannual variability of motion at eight glaciers in this region, due to both long-term force-balance effects and surge dynamics. Multi-decadal velocities at the non-surge-type Kaskawulsh Glacier indicate little change along most of its length, except for the lowermost 10 km where deceleration has been pronounced since the late 1980s in a region that has undergone rapid recent thinning. Interannual variability of surge-type glaciers was high, with year-to-year velocity variations of up to several hundred m a−1. These glaciers were also characterized by distinct patterns of deceleration and/or acceleration along their length.

We develop a formal thin-plate treatment of the viscoelastic flexure of floating ice shelves as an initial step in treating various problems relevant to ice-shelf response to sudden changes of surface loads and applied bending moments (e.g. draining supraglacial lakes, iceberg calving, surface and basal crevassing). Our analysis is based on the assumption that total deformation is the sum of elastic and viscous (or power-law creep) deformations (i.e. akin to a Maxwell model of viscoelasticity, having a spring and dashpot in series). The treatment follows the assumptions of well-known thin-plate approximation, but is presented in a manner familiar to glaciologists and with Glen’s flow law. We present an analysis of the viscoelastic evolution of an ice shelf subject to a filling and draining supraglacial lake. This demonstration is motivated by the proposition that flexure in response to the filling/drainage of meltwater features on the Larsen B ice shelf, Antarctica, contributed to the fragmentation process that accompanied its collapse in 2002.

Controls on the global distribution of surge-type glaciers hold the keys to a better understanding of surge mechanisms. We investigate correlations between the distribution of surge-type glaciers and climatic and glacier geometry variables, using a new global geodatabase of 2317 surge-type glaciers. The highest densities of surge-type glaciers occur within an optimal climatic envelope bounded by temperature and precipitation thresholds. Across all regions with both surge-type and normal glaciers, the former are larger, especially at the cold, dry end of the climatic spectrum. A species distribution model, Maxent, accurately predicts the major clusters of surge-type glaciers using a series of climatic and glacier geometry variables, but under-predicts clusters found outside the climatically optimal surge zone. We interpret the results in terms of a new enthalpy cycle model. Steady states require a balance between enthalpy gains generated by the balance flux and losses via heat conduction and meltwater discharge. This condition can be most easily satisfied in cold, dry environments (thin, low-flux glaciers, efficient conductive heat losses) and warm, humid environments (high meltwater discharges). Intermediate conditions correspond to the optimal surge zone, where neither heat conduction nor runoff can effectively discharge enthalpy gains, and dynamic cycling can result.

Short-term glacier velocity variations typically occur when a water input is accommodated by an increase in the subglacial water pressure. Although these velocity variations have been well documented on many glaciers, few studies have considered them on glaciers where heavy rain and glacier melt occur year-round. This study investigates the relationship between water inputs and glacier velocity on Franz Josef Glacier, New Zealand. We installed six GNSS stations across the lower glacier during austral summer 2010/11 and one station during summer 2012/13. Glacier velocity remained elevated at all stations for ∼7 days following large rain events. During diurnal melt events, we find velocity variations in the early afternoon (12:00–16:00) at 600 m a.s.l. and in the late evening (20:00–01:00) at 400 m a.s.l. We hypothesize that the late-evening velocity variations occurred as an upstream region of high subglacial water pressures and accelerated ice motion propagated downstream. This mechanism may also explain the increased longitudinal compression and transverse extension across the lower glacier during speed-up events. Our results indicate that the subglacial drainage system likely decreases in efficiency upstream and that the water input variability can still cause short-term velocity variations despite the large year-round water inputs.

Despite the important role of supraglacial debris in ablation, knowledge of debris thickness on Himalayan glaciers is sparse. A recently developed method based on reanalysis data and thermal band satellite imagery has proved to be potentially suitable for debris thickness estimation without the need for detailed field data. In this study, we further develop the method and discuss possibilities and limitations arising from its application to a glacier in the Himalaya with scarce in situ data. Surface temperature patterns are consistent for 13 scenes of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and Landsat 7 imagery and correlate well with incoming shortwave radiation and air temperature. We use an energy-balance approach to subtract these radiation or air temperature effects, in order to estimate debris thickness patterns as a function of surface temperature. Both incoming shortwave and longwave radiation are estimated with reasonable accuracy when applying parameterizations and reanalysis data. However, the model likely underestimates debris thickness, probably due to incorrect representation of vertical debris temperature profiles, the rate of heat storage and turbulent sensible heat flux. Moreover, the uncertainty of the result was found to increase significantly with thicker debris, a promising result since ablation is enhanced by thin debris of 1–2 cm.

Marine ice accretes at the base of ice shelves, often infilling open structural weaknesses and increasing ice-shelf stability. However, the timing and location of marine ice formation remain poorly understood. This study determines marine ice source water composition and origin by examining marine ice crystal morphology, water isotope and solute chemistry in ice samples collected from the southern McMurdo Ice Shelf (SMIS), Antarctica. The measured co-isotopic record together with the output of a freezing model for frazil crystals indicate a spatio-temporally varying water source of sea water and relatively fresher water, such as melted meteoric or marine ice. This is in agreement with the occurrence of primarily banded and granular ice crystal facies typical for frazil ice crystals that nucleate in a supercooled mixture of water masses. We propose that marine ice exposed at the surface of SMIS, which experiences summer melt, is routed to the ice-shelf base via the tide crack. Here frazil crystals nucleate in a double diffusion mechanism of heat and salt between two water masses at their salinity-dependent freezing point. Recycling of previously formed marine ice facilitates ice-shelf self-sustenance in a warming climate.

For many regions, glacier inaccessibility results in sparse geometric datasets for use as model initial conditions (e.g. along the central flowline only). In these cases, two-dimensional (2-D) flowline models are often used to study glacier dynamics. Here we systematically investigate the applicability of a 2-D, first-order Stokes approximation flowline model (FLM), modified by shape factors, for the simulation of land-terminating glaciers by comparing it with a 3-D, ‘full’-Stokes ice-flow model (FSM). Based on steady-state and transient, thermomechanically uncoupled and coupled computational experiments, we explore the sensitivities of the FLM and FSM to ice geometry, temperature and forward model integration time. We find that, compared to the FSM, the FLM generally produces slower horizontal velocities, due to simplifications inherent to the FLM and to the underestimation of the shape factor. For polythermal glaciers, those with temperate ice zones, or when basal sliding is important, we find significant differences between simulation results when using the FLM versus the FSM. Over time, initially small differences between the FLM and FSM become much larger, particularly near cold/temperate ice transition surfaces. Long time integrations further increase small initial differences between the two models. We conclude that the FLM should be applied with caution when modelling glacier changes under a warming climate or over long periods of time.

We investigate past and future trends in glacier runoff and the associated hydrological impacts on river runoff in the Hailuogou catchment, a highly glacierized catchment with extensive debris cover in the southeastern Tibetan Plateau, using a catchment-scale glacio-hydrological model. Past trends in various runoff components of the catchment indicate that glacier runoff has been a large component of total runoff, contributing ∼53.4% of total runoff during the period 1952–2013. Future changes in runoff calculated using the outputs of ten global climate models for representative concentration pathway (RCP) 4.5 and RCP8.5 reveal that glacier runoff plays different roles in the water supply of the catchment for the two RCPs, and the discrepancies between the two RCPs increase in the second half of this century, leading to a considerable difference in the hydrological regime of the catchment. In particular, changes are more remarkable under RCP8.5, under which all glaciers are projected to retreat dramatically and total runoff to decrease slightly by the end of this century. An experimental analysis, in which no debris cover is assumed on glacier ablation zones, indicates that excess meltwater from the debris-covered area provides an 8.1% increase in total runoff relative to the no-debris assumption case.

We analyse the long-term (1961–2012) distributed surface mass balance and firn evolution of the Kongsvegen and Holtedahlfonna glacier systems in northwestern Svalbard. We couple a surface energy-balance model to a firn model, with forcing provided from regional climate model output. In situ observational data are used to calibrate model parameters and validate the output. The simulated area-averaged surface mass balance for 1961–2012 is slightly positive (0.08 mw.e.a−1), which only fractionally compensates for mass loss by calving. Refreezing of percolating water in spring/summer (0.13 m w.e. a−1) and stored water in fall/winter (0.18 m w.e. a−1) provides a buffer for runoff. Internal accumulation, i.e. refreezing below the previous year’s summer surface in the accumulation zone, peaks up to 0.22 m w.e. a−1, and is unaccounted for by stake observations. Superimposed ice formation in the lower accumulation zone ranges as high as 0.25 m w.e. a−1. A comparison of the periods 1961–99 and 2000–12 reveals 21% higher annual melt rates since 2000 and a 31% increase in runoff, which can only in part be ascribed to recent warmer and drier conditions. In response to firn line retreat, both albedo lowering (snow/ice–albedo feedback) and lower refreezing rates (refreezing feedback) further amplified runoff.

Observations show that glaciers around the world are in retreat and losing mass. Internationally coordinated for over a century, glacier monitoring activities provide an unprecedented dataset of glacier observations from ground, air and space. Glacier studies generally select specific parts of these datasets to obtain optimal assessments of the mass-balance data relating to the impact that glaciers exercise on global sea-level fluctuations or on regional runoff. In this study we provide an overview and analysis of the main observational datasets compiled by the World Glacier Monitoring Service (WGMS). The dataset on glacier front variations (∼42 000 since 1600) delivers clear evidence that centennial glacier retreat is a global phenomenon. Intermittent readvance periods at regional and decadal scale are normally restricted to a subsample of glaciers and have not come close to achieving the maximum positions of the Little Ice Age (or Holocene). Glaciological and geodetic observations (∼5200 since 1850) show that the rates of early 21st-century mass loss are without precedent on a global scale, at least for the time period observed and probably also for recorded history, as indicated also in reconstructions from written and illustrated documents. This strong imbalance implies that glaciers in many regions will very likely suffer further ice loss, even if climate remains stable.

Glacier volume and ice thickness distribution are important variables for water resource management in Norway and the assessment of future glacier changes. We present a detailed assessment of thickness distribution and total glacier volume for mainland Norway based on data and modelling. Glacier outlines from a Landsat-derived inventory from 1999 to 2006 covering an area of 2692 ± 81 km2 were used as input. We compiled a rich set of ice thickness observations collected over the past 30 years. Altogether, interpolated ice thickness measurements were available for 870 km2 (32%) of the current glacier area of Norway, with a total ice volume of 134 ± 23 km3. Results indicate that mean ice thickness is similar for all larger ice caps, and weakly correlates with their total area. Ice thickness data were used to calibrate a physically based distributed model for estimating the ice thickness of unmeasured glaciers. The results were also used to calibrate volume–area scaling relations. The calibrated total volume estimates for all Norwegian glaciers ranged from 257 to 300 km3.

The Greenland ice sheet is experiencing dramatic melt that is likely to continue with rapid Arctic warming. However, the proportion of meltwater stored before reaching the global ocean remains difficult to quantify. We use NASA MODIS surface reflectance data to estimate river discharge from two West Greenland rivers – the Watson River near Kangerlussuaq and the Naujat Kuat River near Nuuk – over the summers of 2000–12. By comparison with in situ river discharge observations, ‘inundation–discharge’ relations were constructed for both rivers. MODIS-based total annual discharges agree well with total discharge estimated from in situ observations (86% of summer discharge in 2009 to 96% in 2011 at the Watson River, and 106% of total discharge in 2011 to 104% in 2012 at the Naujat Kuat River). We find, however, that a time-lapse camera, deployed at the Watson River in summer 2012, better captures the variations in observed discharge, benefiting from fewer data gaps due to clouds. The MODIS-derived estimates indicate that summer discharge has not significantly increased over the last decade, despite a strong warming trend. Also, meltwater runoff estimates derived from the regional climate model RACMO2/GR for the drainage basins are higher than our reconstructions of river discharge. These results provide indirect evidence for a considerable component of water storage within the glacio-hydrological system.

Phase velocities and attenuation in snow cannot be explained by the widely used elastic or viscoelastic models for acoustic wave propagation. Instead, Biot’s model of wave propagation in porous materials should be used. However, the application of Biot’s model is complicated by the large property space of the underlying porous material. Here constant properties for ice and air, and empirical relationships are used to estimate unknown porous properties from snow porosity. Using this set of equations, phase velocities and plane wave attenuation of shear- and compressional waves are predicted as functions of porosity or density. For light snow the peculiarity was found that the velocity of the first compressional wave is lower than that of the second compressional wave that is commonly referred to as the ‘slow’ wave. The reversal of the velocities comes with an increase of attenuation for the first compressional wave. This is in line with the common observation that sound is strongly absorbed in light snow. The results have important implications for the use of acoustic waves to evaluate snow properties and to numerically simulate wave propagation in snow.

Simulations of the post-Little Ice Age evolution of three outlet glaciers of southeast Vatnajökull, Iceland – Skálafellsjökull, Heinabergsjökull and Fláajökull – are presented. A coupled shallow-ice-approximation ice-flow and degree-day mass-balance model is applied that is calibrated with a 14 year record of in situ mass-balance measurements. The measured mass balance cannot be realistically represented by constant horizontal and vertical precipitation gradients. High-resolution (1 km) precipitation fields, derived from downscaled orographic atmospheric circulation models of precipitation, are required to capture the spatial variation of the winter mass balance. The observed ice volume around 1890 (15–30% larger than in 2000) can be simulated with 1°C lower temperatures and a 20% reduction in the annual precipitation, relative to the reference climate period, 1980–2000. The sensitivity of each glacier’s annual balance to a change in temperature is −1.51 to −0.97 m w.e. a−1 °C−1 and +0.16 to +0.65 m w.e. a−1 for a 10% increase in precipitation. A steady-state experiment applying a step increase in temperature of 2°C (3°C), and precipitation increase of 10%, results in a >50% (80–90%) decrease in ice volume.