Crevasses pose severe risks for mountaineers and field glaciologists. Smaller cracks between 0.5 and 2 m are still dangerous, but often not visible in medium resolution satellite imagery. If they are snow covered, they are completely undetectable by optical sensors. We set out to develop an approach to detect potentially crevassed areas by a minimum of geometric data, and to make the method generally applicable to glacier regions. On Vernagtferner, we compared a reference dataset of crevasses observed in high-resolution optical imagery with the curvature of the ice surface and the spatial gradients in driving stress. Both parameters can be derived from a digital surface model and a bedrock model, derived from ice thickness measurements. The correlation patterns show that crevasses preferably form in convex areas and in areas where the driving stress rapidly increases. This corresponds with the theory of crevasse formation. Although the method still misclassifies larger parts, the approach has the potential to define probable non-crevassed areas as well as to aid the planning of safe routes.

To understand the mechanism of simultaneous drainage event related to supraglacial lakes on a debris-covered glacier, we investigated water-level variations of supraglacial lakes on the southern Inylchek Glacier in Kyrgyzstan. To examine these variations, we used daily aerial images for 2017–2019 from an uncrewed aerial vehicle that were converted to 15 cm-digital surface models and ortho-images. Our main results are as follows: (1) When one lake drained, the water levels of other lakes simultaneously increased, indicating that drainage water is shared with several lakes through a main englacial conduit. In one drainage event, a branched off englacial conduit clearly connected to a main englacial conduit. (2) Sometimes several lakes discharged simultaneously, indicating that several lakes had connected to a main englacial conduit that had opened. Such cases can cause larger-scale drainage than that from the opening of a branched off englacial conduit. (3) Simultaneous drainage occurred twice in the same year, each time through a different conduit, indicating that the main englacial conduit can be abandoned and reused. (4) In some lakes, the water level on the hydraulic gradient line increased gradually with nearly the same increase rate just before drainage. Such an increase may be an indicator of a possible simultaneous drainage event.

In the afternoon of 30 January 2017, a catastrophic outburst flood occurred in the Larsemann Hills (Princess Elizabeth Land, East Antarctica). The rapid drainage of both a thin supraglacial layer of water (near Boulder Lake) and Lake Ledyanoe into the englacial Lake Dålk provoked its overfill and outburst. As a result, a depression of 183 m × 220 m was formed in the place where Lake Dålk was located. This study summarises and clarifies the current state of knowledge on the flood that occurred in 2017. We present a phenomenological model of depression formation. We specify the reasons for the outburst of the system of lakes Boulder, Ledyanoe and Dålk. In addition, we carry out mathematical modelling of the outburst of each of the three lakes and estimate the flood severity. Outburst hydrographs, channel diameters, volume and duration of floods were calculated. Particular simulation results were validated with field data. In conclusion, we give an overview of the new outburst cycle of the lake system, which began in 2020 with the drainage of the lakes Boulder and Ledyanoe, and the new formation of Lake Dålk. Further research is required to improve our understanding of the lake system responses to changing external factors.

Glacier surges are periodic episodes of mass redistribution characterized by dramatic increases in ice flow velocity and, sometimes, terminus advance. We use optical satellite imagery to document five previously unexamined surge events of Sít’ Kusá (Turner Glacier) in the St. Elias Mountains of Alaska from 1983 to 2013. Surge events had an average recurrence interval of ~5 years, making it the shortest known regular recurrence interval in the world. Surge events appear to initiate in the winter, with speeds reaching up to ~25 m d−1. The surges propagate down-glacier over ~2 years, resulting in maximum thinning of ~100 m in the reservoir zone and comparable thickening at the terminus. Collectively, the rapid recurrence interval, winter initiation and down-glacier propagation suggest Sít’ Kusá's surges are driven by periodic changes in subglacial hydrology and glacier sliding. Elevation change observations from the northern tributary show a kinematic disconnect above and below an icefall located 23 km from the terminus. We suggest the kinematic disconnect inhibits drawdown from the accumulation zone above the icefall, which leads to a steady flux of ice into the reservoir zone, and contributes to the glacier's exceptionally short recurrence interval.

As a result of climate warming, glacial meltwater discharge has been increasing in Greenland. During the summers of 2015 and 2016, there were rapid increases in discharge from Qaanaaq Glacier in northwestern Greenland. These discharges resulted in floods that destroyed the road linking the settlement of Qaanaaq to Qaanaaq Airport. Field measurements were performed and a numerical model of glacier runoff was developed to quantify these discharges. The high discharge associated with the 2015 flood, estimated at 9.1 m3 s−1 (hourly mean), resulted from intensive glacier melting due to warm air temperature and strong winds, while the high discharge associated with the 2016 flood resulted from heavy rainfall (90 mm d−1) that led to a peak discharge estimated at 19.9 m3 s−1. The developed model, when used to investigate future glacier runoff under warming conditions, revealed a nonlinear increase in glacial melt with increasing temperature. Additionally, the model forecasted a threefold increase in total summer discharge, owing to a 4 °C rise in temperature. Thus, this study quantified the impact of a changing climate on glacier runoff, which gives insight into future risks of flood hazards along the coast of Greenland.

During the last few decades, the lake-terminating glaciers in the Himalaya have receded faster than the land-terminating glaciers as proglacial lakes have exacerbated the mass loss of their host glaciers. Monitoring the impacts of glacier recession and dynamics on lake extent and water volume provides an approach to assess the mass interplay between glaciers and proglacial lakes. We describe the recession of Longbasaba Glacier and estimate the mass wastage and its contribution to the water volume of its proglacial lake. The results show that the glacier area has decreased by 3% during 1988–2018, with a more variable recession prior to 2008 than in the last decade. Longbasaba Lake has expanded by 164% in area and 237% in water volume, primarily as a result of meltwater inflow produced from surface lowering of the glacier. Over the periods 1988–2000 and 2000–18, the mass loss contributed by glacier thinning has decreased from 81 to 61% of the total mass loss, accompanied by a nearly doubled contribution from terminus retreat. With the current rate of retreat, Longbasaba glacier is expected to terminate in its proglacial lake for another four decades. The hazard risk of this lake is expected to continue to increase in the near future because of the projected continued glacier mass loss and related lake expansion.

Using an ensemble of close- and long-range remote sensing, lake bathymetry and regional meteorological data, we present a detailed assessment of the geometric changes of El Morado Glacier in the Central Andes of Chile and its adjacent proglacial lake between 1932 and 2019. Overall, the results revealed a period of marked glacier down wasting, with a mean geodetic glacier mass balance of −0.39 ± 0.15 m w.e.a−1 observed for the entire glacier between 1955 and 2015 with an area loss of 40% between 1955 and 2019. We estimate an ice elevation change of −1.00 ± 0.17 m a−1 for the glacier tongue between 1932 and 2019. The increase in the ice thinning rates and area loss during the last decade is coincident with the severe drought in this region (2010–present), which our minimal surface mass-balance model is able to reproduce. As a result of the glacier changes observed, the proglacial lake increased in area substantially between 1955 and 2019, with bathymetry data suggesting a water volume of 3.6 million m3 in 2017. This study highlights the need for further monitoring of glacierised areas in the Central Andes. Such efforts would facilitate a better understanding of the downstream impacts of glacier downwasting.

On 30 September 2017, an Air France Airbus A380-800 suffered a failure of its fourth engine while over Greenland. This failure resulted in the loss of the engine fan hub, fan blades and surrounding structure. An initial search recovered 30 pieces of light debris, but the primary part of interest, a ~220 kg titanium fan hub, was not recovered because it had a different fall trajectory than the light debris, impacted into the ice-sheet's snow surface, and was quickly covered by drifting snow. Here we describe the methods used for the detection of the fan hub and details of the field campaigns. The search area included two crevasse fields of at least 50 snow-covered crevasses 1 to ~30 m wide with similar snow bridge thicknesses. After 21 months and six campaigns, using airborne synthetic aperture radar, ground-penetrating radar, transient electromagnetics and an autonomous vehicle to survey the crevasse fields, the fan hub was found within ~1 m of a crevasse at a depth of ~3.3 to 4 m and was excavated with shovels, chain saws, an electric winch, sleds and a gasoline heater, by workers using fall-arrest systems.

Since the mid-1970s radio-echo soundings have been conducted on Iceland's temperate glaciers. Since then, low-frequency radar technology has furthered the study of most of the island's ice caps. Their masses and volumes have been quantified and detailed subglacial topographic maps produced which demarcate glacial drainage basins and identify subglacial lakes and volcanoes. Even internal tephra layers have been charted. The resulting data have been used to force and validate models of past and future glacier evolution. Many practical applications in glacier hydrology have come into being, including hydropower management, road and bridge planning and the prediction of catastrophic flood paths from subglacial eruption sites. Finally, emerging landscapes can now be foreseen in places where glaciers may soon disappear. These achievements would not have been possible without the advances in RES technology.

Herein, we present results obtained from time-lapse imagery acquired by a digital single-lens reflex camera during 2014–18 used to monitor the Planpincieux Glacier on the Italian side of the Grandes Jorasses (Mont Blanc massif). We processed the images using image cross-correlation to measure the surface kinematics of the most active lobe of the glacier that presents a high velocity and frequent ice detachments. During the monitoring, we observed two or three periods of sharp acceleration per year that culminated with large break-offs followed by analogous decelerations. Overall, we registered more than 350 failures with a volume >100 m3, of which, 14 events had volumes larger than 10 000 m3. The study identified a monotonic positive relationship between the velocity and failure volume that may be used to estimate the volume of the collapses before an event. We identified the thresholds of velocity and acceleration that characterise the activation of the speed-up periods. The study allowed the characterisation of three different instability processes that lead to the break-off of ice chunks from the glacier terminus: (i) disaggregation, (ii) slab fracture and (iii) water tunnelling failure which can be differentiated based on the rheology, the volume involved and the trigger process.

Donjek Glacier has an unusually short and regular surge cycle, with eight surges identified since 1935 from aerial photographs and satellite imagery with a ~12 year repeat interval and ~2 year active phase. Recent surges occurred during a period of long-term negative mass balance and cumulative terminus retreat of 2.5 km since 1874. In contrast to previous work, we find that the constriction where the valley narrows and bedrock lithology changes, 21 km from the terminus, represents the upper limit of surging, with negligible surface speed or elevation change up-glacier from this location. This positions the entire surge-type portion of the glacier in the ablation zone. The constriction geometry does not act as the dynamic balance line, which we consistently find at 8 km from the glacier terminus. During the 2012–2014 surge event, the average lowering rate in the lowest 21 km of the glacier was 9.6 m a−1, while during quiescence it was 1.0 m a−1. Due to reservoir zone refilling, the ablation zone has a positive geodetic balance in years immediately following a surge event. An active surge phase can result in a strongly negative geodetic mass balance over the surge-type portion of the glacier.

Twin glaciers collapsed in 2016 near Aru Co, western Tibet and caused extreme loss to human beings. In this study, we attempted to track the dynamics of glaciers in the region, for example the glacier area and mass changes in Aru Co for the period 1971–2016, which were determined using topographic maps and Landsat images and ASTER-derived DEMs (2011–16), the Shuttle Radar Terrain Mission DEM (2000) and topographic maps (1971). Our results showed that the glacier area of Aru Co decreased by −0.4 ± 4.1% during 1971–2016. The geodetic mass-balance results showed that the glaciers in Aru Co lost mass at a rate of −0.15 ± 0.30 m w.e. a−1 during 1971–99, while they gained mass at a rate of 0.33 ± 0.61 m w.e. a−1 for the period 1999–2016. The twin glaciers experienced a larger negative mass budget than the others in the region before 1999. This process produced large amounts of meltwater, followed by a sustained increase in the meltwater on the pressure melting point, possibly in response to a period of positive mass balance (1999–2016) and then, transferred to the glacier bed until the glaciers collapsed.

Kyagar glacier is located in the Chinese Karakoram mountains. The glacier tongue entirely blocks the riverbed in the upper Shaksgam valley and impounds a glacial lake, which was the source of several violent and disastrous glacial lake outburst floods (GLOFs). A GLOF early warning system was implemented between 2011 and 2013. We present an integrative analysis of the hazard potential of Kyagar lake, taking into account the ice flow dynamics of Kyagar glacier as well as the recent surface mass-balance response to climate change. Comparison of two high-resolution digital elevation models (DEMs) for the ice dam shows surface lowering rates of >5ma– 1 between 2002 and 2011, leading to a significant reduction in the maximum potential lake volume. However, two DEMs covering the entire glacier for the period 2000–10 indicate mass gains in its central part, and flow speed measurements show an acceleration in this region. This pattern of local ice-thickness changes combined with varying ice flow velocities is typical for surge-type glaciers. The velocity of the glacier surface and of the ice dam between 2011 and 2012 are analyzed at high temporal and spatial resolution, based on feature tracking of synthetic aperture radar (SAR) images.

Very large volumes of ice break off regularly from Taconnaz hanging glacier, French Alps. During winter, when the snow mantle is unstable, these collapses can trigger very large avalanches that represent a serious threat to inhabited areas below. Photogrammetric measurements have been performed over 1 year to assess the volume and frequency of the largest collapses. Major collapses occur when the glacier reaches a critical geometry. After a major ice collapse, the glacier is in a minimal position and subsequently recharges over 6 months to reach the maximum position again. This critical geometry is a necessary but not sufficient condition for further large collapses. Large collapses do not systematically occur in the maximum position, as ice is often removed by disintegration into small ice blocks. For two major collapses, the volume of ice breaking off has been assessed at ~275 000 m3. Photogrammetric measurements were used to determine an ice flux of 820 000 m3 a–1 through the studied ice stream, in agreement with an assessment based on ice-flow modeling. This ice flux estimation was used to determine the average ice volumes breaking off over surveyed periods.

A model subglacial drainage system, coupled to an ice-dammed reservoir that receives a time-varying meltwater input, is described and analysed. In numerical experiments an ice-marginal lake drains through a subglacial channel, producing periodic floods, and fills with meltwater at a rate dependent on air temperature, which varies seasonally with a peak value of Tm. The analysis reveals regions of Tm parameter space corresponding to ‘mode locking’, where flood repeat time is independent of Tm; resonance, where decreasing Tm counter-intuitively increases flood size; and chaotic dynamics, where flood cycles are sensitive to initial conditions, never repeat and exhibit phase-space mixing. Bifurcations associated with abrupt changes in flood size and timing within the year separate these regions. This is the first time these complex dynamics have been displayed by a glaciohydraulic model and these findings have implications for our understanding of ice-marginal lakes, moulins and subglacial lakes.

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.

To better understand the risk of disasters due to glacial lake outburst floods (GLOFs), we synthetically analyze the spatial distribution and evolution of moraine-dammed lakes and potentially dangerous glacial lakes (PDGLs) in the Chinese Himalaya. Our county-based assessment of GLOF disaster risk combines PDGL outburst hazard, regional exposure, vulnerability of exposed elements and adaptation capability (risk management) using the analytic hierarchy process. We synthetically analyze the disaster risk using the weighted comprehensive method. Remote-sensing data show there are 329 moraine-dammed lakes (>0.02 km2; total area 125.43 km2) in the Chinese Himalaya, of which 116 (total area 49.49 km2) are identified as PDGLs. The zones at highest risk of GLOF disaster are mainly located in Nyalam, Tingri, Dinggyê, Lhozhag, Kangmar and Zhongba, in the mid-eastern Himalaya. Lowest-risk zones are located in the eastern Himalaya. On the county scale, Lhozhag and Lhunze have the highest hazard degrees and exposure, while Zhongba and Zando have the highest degree of vulnerability and lowest adaptation capacity. Our regionalization results for GLOF disaster risk are consistent with the distribution of historical disaster sites across the Chinese Himalaya.

In snow avalanche long-term forecasting, existing risk-based methods remain difficult to use in a real engineering context. In this work, we expand a quasi analytical decisional model to obtain simple formulae to quantify risk and to perform the optimal design of an avalanche dam in a quick and efficient way. Specifically, the exponential runout model is replaced by the Generalized Pareto distribution (GPD), which has theoretical justifications that promote its use for modelling the different possible runout tail behaviours. Regarding the defence structure/flow interaction, a simple law based on kinetic energy dissipation is compared with a law based on the volume stored upstream of the dam, whose flexibility allows us to cope with various types of snow. We show how a detailed sensitivity study can be conducted, leading to intervals and bounds for risk estimates and optimal design values. Application to a typical case study from the French Alps, highlights potential operational difficulties and how they can be tackled. For instance, the highest sensitivity to the runout tail type and interaction law is found at abscissas of legal importance for hazard zoning (return periods of 10–1000 a), a crucial result for practical purposes.

The ablation areas of debris-covered glaciers typically consist of a complex mosaic of surface features with contrasting processes and rates of mass loss. This greatly complicates glacier response to climate change, and increases the uncertainty of predictive models. In this paper we present a series of high-resolution DEMs and repeat lake bathymetric surveys on Ngozumpa Glacier, Nepal, to study processes and patterns of mass loss on a Himalayan debris-covered glacier in unprecedented detail. Most mass loss occurs by melt below supraglacial debris, and melt and calving of ice cliffs (backwasting). Although ice cliffs cover only ~5% of the area of the lower tongue, they account for 40% of the ablation. The surface debris layer is subject to frequent re-distribution by slope processes, resulting in large spatial and temporal differences in debris-layer thickness, enhancing or inhibiting local ablation rates and encouraging continuous topographic inversion. A moraine-dammed lake on the lower glacier tongue (Spillway Lake) underwent a period of rapid expansion from 2001 to 2009, but later experienced a reduction of area and volume as a result of lake level lowering and sediment redistribution. Rapid lake growth will likely resume in the near future, and may eventually become up to 7 km long.