Surface height changes above three previously undetected subglacial lakes in northeastern Greenland are documented using CryoSat, DEMs and ICESat-2. Between 7 February and 6 March 2012, the central ice region (22.6 km2) above the largest lake dropped by ~37 m followed by a further drop of 12 m in the following 29 days. This implies a subglacial water outflow, or jökulhlaup, of at least 1 km3 at rates of hundreds of cubic meters per second. A comparable outflow occurred again between 23 July and 15 September 2019, with smaller outflows in the fall of 2014 and 2016. In contrast, a second smaller subglacial lake at a higher elevation had two subglacial outbursts of ~0.3 km3 in 2012 and 2019 but the lake filling was gradual and not strongly seasonal or episodic. Water remained in both lakes after the outflows but this may not be the case for the third smallest and lowest subglacial lake. While there appears to be some hydrological link between the three lakes, the flux of water moving under the ice in this area appears to be larger than would be expected from local summer melt. However, the source of the excess water remains uncertain.

The number of glacial lakes has grown globally concurrently with the retreat of glaciers in the last few decades, increasing the risk of potentially hazardous glacial lake outburst floods (GLOFs) and posing a threat to downstream communities. Norway has several known ice-dammed lakes that produce repeated GLOFs but as we show here, the existing GLOF database is incomplete and needs to be improved through continuous monitoring of glaciers and glacial lakes. This study examines the case of an ice cap in central Norway hosting at least four drainage-prone lakes. We reconstruct the sequence of lake drainage events through a combination of remote sensing, ground-truthing and citizen science while scrutinising the applicability of the PlanetScope imagery vs Sentinel-2 and Landsat-8 OLI products. As opposed to the Landsat imagery that often fails to resolve even the largest glacial lakes of Tystigbreen, both PlanetScope and Sentinel-2 are helpful in identifying previously unrecognised glacial lakes and undocumented drainage events. Our analysis suggests that a fusion of the two satellite products may be beneficial for automated tracking of glacial lake changes. We also demonstrate that local knowledge and systematic involvement of citizens in data collection have a potential to enrich GLOF databases.

The rapid formation of glacial lakes is one of the most conspicuous landscape changes caused by atmospheric warming in glacierised regions. However, relatively little is known about the history and current state of glacial lakes in the High Arctic. This study aims to address this issue by providing the first inventory of glacial lakes in Svalbard, focusing in particular on the post-Little Ice Age evolution of glacial lakes and their typology. To do so, we used aerial photographs and satellite imagery together with archival topographic data from 1936 to 2020. The inventory comprises the development of 566 glacial lakes (146 km2) that were still in direct contact with glaciers during the period 2008–2012. The results show a consistent increase in the total area of glacial lakes from the 1930s to 2020 and suggest an apparent link between climatic and geological factors, and the formation of specific lake dam types: moraine, ice, or bedrock. We also detected 134 glacial lake drainage events that have occurred since the 1930s. This study shows that Svalbard has one of the highest rates of glacial lake development in the world, which is an indicator of the overall dynamics of landscape change in the archipelago in response to climate change.

Proglacial braided river systems discharge large volumes of meltwater from ice sheets and transport coarse-grained sediments from the glaciated areas to the oceans. Here, we test the hypothesis if high-energy hydrological events can leave distinctive signatures in the sedimentary record of braided river systems. We characterize the morphology and infer a mode of formation of a 25 km long and 1–3 km wide Early Pleistocene incised valley recently imaged in 3-D seismic data in the Hoop area, SW Barents Sea. The fluvial system, named Bjørnelva River Valley, carved 20 m deep channels into Lower Cretaceous bedrock at a glacial paleo-surface and deposited 28 channel bars along a paleo-slope gradient of ~0.64 m km−1. The landform morphologies and position relative to the paleo-surface support that Bjørnelva River Valley was formed in the proglacial domain of the Barents Sea Ice Sheet. Based on valley width and valley depth, we suggest that Bjørnelva River Valley represents a braided river system fed by violent outburst floods from a glacial lake, with estimated outburst discharges of ~160 000 m3 s−1. The morphological configuration of Bjørnelva River Valley can inform geohazard assessments in areas at risk of outburst flooding today and is an analogue for landscapes evolving in areas currently covered by the Greenland and Antarctic ice sheets.

Rapidly-rising jökulhlaups, or glacial outburst floods, are a phenomenon with a high potential for damage. The initiation and propagation processes of a rapidly-rising jökulhlaup are still not fully understood. Seismic monitoring can contribute to an improved process understanding, but comprehensive long-term seismic monitoring campaigns capturing the dynamics of a rapidly-rising jökulhlaup have not been reported so far. To fill this gap, we installed a seismic network at the marginal, ice-dammed lake of the A.P. Olsen Ice Cap (APO) in NE-Greenland. Episodic outbursts from the lake cause flood waves in the Zackenberg river, characterized by a rapid discharge increase within a few hours. Our 6 months long seismic dataset comprises the whole fill-and-drain cycle of the ice-dammed lake in 2012 and includes one of the most destructive floods recorded so far for the Zackenberg river. Seismic event detection and localization reveals abundant surface crevassing and correlates with changes of the river discharge. Seismic interferometry suggests the existence of a thin basal sedimentary layer. We show that the ballistic part of the first surface waves can potentially be used to infer medium changes in both the ice body and the basal layer. Interpretation of time-lapse interferograms is challenged by a varying ambient noise source distribution.

How landscapes respond to, and evolve from, large jökulhlaups (glacial outburst floods) is poorly constrained due to limited observations and detailed monitoring. We investigate how melt of glacier ice transported and deposited by multiple jökulhlaups during the 2010 eruption of Eyjafjallajökull, Iceland, modified the volume and surface elevation of jökulhlaup deposits. Jökulhlaups generated by the eruption deposited large volumes of sediment and ice, causing significant geomorphic change in the Gígjökull proglacial basin over a 4-week period. Observation of these events enabled robust constraints on the physical properties of the floods which informs our understanding of the deposits. Using ground-based LiDAR, GPS observations and the satellite-image-derived ArcticDEMs, we quantify the post-depositional response of the 60 m-thick Gígjökull sediment package to the meltout of buried ice and other geomorphic processes. Between 2010 and 2016, total deposit volume reduced by −0.95 × 106 m3 a−1, with significant surface lowering of up to 1.88 m a−1. Surface lowering and volumetric loss of the deposits is attributed to three factors: (i) meltout of ice deposited by the jökulhlaups; (ii) rapid melting of the buried Gígjökull glacier snout; and (iii) incision of the proglacial meltwater system into the jökulhlaup deposits.

Mountain glaciers and ice sheets often host marginal and subglacial lakes that are hydraulically connected through subglacial drainage systems. These lakes exhibit complex dynamics that have been the subject of models for decades. Here we introduce and analyze a model for the evolution of glacial lakes connected by subglacial channels. Subglacial channel equations are supplied with effective pressure boundary conditions that are determined by a simple lake model. While the model can describe an arbitrary number of lakes, we solve it numerically with a finite element method for the case of two connected lakes. We examine the effect of relative lake size and spacing on the oscillations. Complex oscillations in the downstream lake are driven by discharge out of the upstream lake. These include multi-peaked and anti-phase filling–draining events. Similar filling–draining cycles have been observed on the Kennicott Glacier in Alaska and at the confluence of the Whillans and Mercer ice streams in West Antarctica. We further construct a simplified ordinary differential equation model that displays the same qualitative behavior as the full, spatially-dependent model. We analyze this model using dynamical systems theory to explain the appearance of filling–draining cycles as the meltwater supply varies.

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.

Discharge and water temperature measurements in the Skaftá river and measurements of the lowering of the ice over the subglacial lake at the western Skaftá cauldron, Vatnajökull, Iceland, were made during a rapidly rising glacial outburst flood (jökulhlaup) in September 2006. Outflow from the lake, flood discharge at the glacier terminus and the transient subglacial volume of floodwater during the jökulhlaup are derived from these data. The 40 km long initial subglacial path of the jökulhlaup was mainly formed by lifting and deformation of the overlying ice, induced by water pressure in excess of the ice overburden pressure. Melting of ice due to the heat of the floodwater from the subglacial lake and frictional heat generated by the dissipation of potential energy in the flow played a smaller role. Therefore this event, like other rapidly rising jökulhlaups, cannot be explained by the jökulhlaup theory of Nye (1976). Instead, our observations indicate that they can be explained by a coupled subglacial-sheet–conduit mechanism where essentially all of the initial flood path is formed as a sheet by the propagation of a subglacial pressure wave.

The subglacial hydrologic system exerts strong controls on the dynamics of the overlying ice, yet the parameters that govern the evolution of this system are not widely known or observable. To gain a better understanding of these parameters, we invert a spatially averaged model of subglacial hydrology from observations of ice surface velocity and outlet stream discharge at Kennicott Glacier, Wrangell Mountains, AK, USA. To identify independent parameters, we formally non-dimensionalize the forward model. After specifying suitable prior distributions, we use a Markov-chain Monte Carlo algorithm to sample from the distribution of parameter values conditioned on the available data. This procedure gives us not only the most probable parameter values, but also a rigorous estimate of their covariance structure. We find that the opening of cavities due to sliding over basal topography and turbulent melting are of a similar magnitude during periods of large input flux, though turbulent melting also exhibits the greatest uncertainty. We also find that both the storage of water in the englacial system and the exchange of water between englacial and subglacial systems are necessary in order to explain both surface velocity observations and the relative attenuation in the amplitude of diurnal signals between input and output flux observations.

GPS campaigns on glaciers during jökulhlaups show how subglacial floods affect glacier motion and shed light on the dynamics of such floods. Three such campaigns have been carried out on southern and western Vatnajökull, southeast Iceland, over known jökulhlaup paths. Two slowly rising jökulhlaups from Grímsvötn and two rapidly rising jökulhlaups from the western and eastern Skaftá cauldrons were captured in these campaigns, with maximum discharge ranging from 240 to 3300 m3 s−1. Glacier surface movements measured in these campaigns are presented along with the corresponding discharge curves. The measurements are interpreted as indicating: (1) initiation of rapidly rising jökulhlaups with a propagating subglacial pressure wave, (2) decreased glacier basal friction during jökulhlaups, (3) subglacial accumulation of water in slowly rising jökulhlaups and (4) lifting of the glacier caused by subglacial water pressure exceeding overburden in both rapidly and slowly rising jökulhlaups. The latter two observations are inconsistent with assumptions that are typically made in theoretical and numerical modelling of jökulhlaups. Both viscous and elastic deformation of the glacier as well as turbulent hydraulic fracture at the ice/bedrock interface are important in the dynamics of the subglacial pressure wave at the front of rapidly rising jökulhlaups.

To simulate a glacial lake outburst flood, we employ a comprehensive physically based general two-phase mass flow model (Pudasaini, 2012). This model accounts for a strong interaction between the solid and fluid phases and incorporates buoyancy and other dominant physical aspects of the mass flows such as enhanced non-Newtonian viscous stress, virtual mass force and generalized drag. Our real two-phase mass flow simulation describes explicit evolution of the solid and fluid phases and the debris bulk as a whole, akin to torrential debris flows or debris floods during glacial lake outburst floods (GLOFs). The emptying of a lake following rapid collapse of a restraining dam, the consequent downslope motion of a mixed solid–fluid mass, and the tendency of the mass to form extruding plumes are analyzed in detail for different flow configurations, volumes, conduit geometries and boundary conditions. The solid and fluid phases evolve completely differently and reveal fundamentally different dynamical behaviours. During the flow, the relatively long fluid tail follows the solid-rich dense frontal surge head. The bulk debris develops into a frontal and side levee as derived from the initial frontal moraine dam. Results show that our high-resolution, unified simulation strategies and the advanced model equations can be applied to study the flow dynamics of a wide range of geophysical mass flows such as snow and rock–ice avalanches, debris flows and flash floods as well as GLOFs. This may help substantially in forming a basis for appropriate mitigation measures against potential natural hazards in high mountain slopes and valleys.

Glaciers in the Bugyai Kangri are located in a transition zone from southeast Tibet, where monsoonal temperate glaciers dominate, to inner Tibet, where continental glaciers dominate. Here we analyze glacier and glacial lake changes in this region using multi-year inventories based on Landsat images from 1981–2013. Results show that the total area of 141 glaciers in the region decreased by 30.44 ±0.89 km2 from 198.35 ±9.54 km2 (1980s) to 167.93 ±4.52 km2 (2010s). The annual area shrinkage rate (–0.48% a–1) is lower than that reported for southeastern Tibet but higher than that of inner Tibet. Both the number and total area of glacial lakes increased between 1981 and 2013. Among all lakes, proglacial lakes contribute most (~81 %) to the expansion. The total area of ten proglacial lakes increased by 150.3 ± 13.17% and of these ten lakes the four that expanded most sharply showed increased calving at their upper margins, resulting in more rapid retreat of lake-terminating glaciers than land-terminating glaciers. Owing to rapid calving, several lakes may undergo further growth in the near future, increasing the potential risk of glacial lake outburst floods.