Supraglacial Networks Research Articles

Seasonal evolution of the supraglacial drainage network at Humboldt Glacier, northern Greenland, between 2016 and 2020

Abstract. Supraglacial rivers and lakes are important for the routing and storage of surface meltwater during the summer melt season across the Greenland Ice Sheet (GrIS) but remain poorly mapped and quantified across the northern part of the ice sheet, which is rapidly losing mass. Here we produce, for the first time, a high-resolution record of the supraglacial drainage network (including both rivers and lakes) and its seasonal behaviour at Humboldt Glacier, a wide-outlet glacier draining a large melt-prone hydrologic catchment (13 488 km2), spanning the period 2016 to 2020 using 10 m spatial resolution Sentinel-2 imagery. Our results reveal a perennially extensive yet interannually variable supraglacial network extending from an elevation of 200 m a.s.l. to a maximum of ∼ 1440 m a.s.l. recorded in 2020, with limited development of the network observed in the low-melt years of 2017 and 2018. The supraglacial drainage network is shown to cover an area ranging between 966 km2 (2018) and 1566 km2 (2019) at its maximum seasonal extent, with spatial coverage of up to 2685 km2 recorded during the early phases of the melt season when a slush zone is most prominent. Up-glacier expansion and the development of an efficient supraglacial drainage network as surface runoff increases and the snowline retreats is clearly visible. Preconditioning of the ice surface following a high-melt year is also observed, with an extreme and long-lasting 2019 melt season and over-winter persistence of liquid lakes, followed by low snow accumulation the following spring, culminating in earlier widespread exposure of the supraglacial drainage network in 2020 compared to other years. This preconditioning is predicted to become more common with persistent warmer years into the future. Overall, this study provides evidence of a persistent, yet dynamic, supraglacial drainage network at this prominent northern GrIS outlet glacier and advances our understanding of such hydrologic processes, particularly under ongoing climatic warming and enhanced runoff.

Extraordinary runoff from the Greenland ice sheet in 2012 amplified by hypsometry and depleted firn retention

Abstract. It has been argued that the infiltration and retention of meltwater within firn across the percolation zone of the Greenland ice sheet has the potential to buffer up to ∼ 3.6 mm of global sea-level rise (Harper et al., 2012). Despite evidence confirming active refreezing processes above the equilibrium line, their impact on runoff and proglacial discharge has yet to be assessed. Here, we compare meteorological, melt, firn stratigraphy and discharge data from the extreme 2010 and 2012 summers to determine the relationship between atmospheric forcing and melt runoff at the land-terminating Kangerlussuaq sector of the Greenland ice sheet, which drains into the Watson River. The 6.8 km3 bulk discharge in 2012 exceeded that in 2010 by 28 %, despite only a 3 % difference in net incoming melt energy between the two years. This large disparity can be explained by a 10 % contribution of runoff originating from above the long-term equilibrium line in 2012 caused by diminished firn retention. The amplified 2012 response was compounded by catchment hypsometry; the disproportionate increase in area contributing to runoff as the melt-level rose high into the accumulation area.Satellite imagery and aerial photographs reveal an extensive supraglacial network extending 140 km from the ice margin that confirms active meltwater runoff originating well above the equilibrium line. This runoff culminated in three days with record discharge of 3100 m3 s−1 (0.27 Gt d−1) that peaked on 11 July and washed out the Watson River Bridge. Our findings corroborate melt infiltration processes in the percolation zone, though the resulting patterns of refreezing are complex and can lead to spatially extensive, perched superimposed ice layers within the firn. In 2012, such layers extended to an elevation of at least 1840 m and provided a semi-impermeable barrier to further meltwater storage, thereby promoting widespread runoff from the accumulation area of the Greenland ice sheet that contributed directly to proglacial discharge and global sea-level rise.

Supraglacial channel inception: Modeling and processes

Supraglacial drainage systems play a key role in glacial hydrology. Nevertheless, physical processes leading to spatial organization in supraglacial networks are still an open issue. In the present work we thus address from a quantitative point of view the question of what is the physics leading to widely observed patterns made up of evenly spaced channels. To this aim, we set up a novel mathematical model describing a condition antecedent channel formation, i.e., the down-glacier flow of a distributed meltwater film. We then perform a linear stability analysis to assess whether the ice-water interface undergoes a morphological instability compatible with observed patterns. The instability is detected, its features depending on glacier surface slope, ice friction factor, and water as well as ice thermal conditions. By contrast, in our model channel spacing is solely hydrodynamically driven and relies on the interplay between pressure perturbations, flow depth response, and Reynolds stresses. Geometrical features of the predicted pattern are quantitatively consistent with available field data. The hydrodynamic origin of supraglacial channel morphogenesis suggests that alluvial patterns might share the same physical controls.