Review article: The hydrology of debris-covered glaciers – state of the science and future research directions

authorsversionNon peer reviewe

Sub-regional variability in the influence of ice-contact lakes on Himalayan glaciers

Ice-contact lakes modify glacier geometry and dynamics by shifting the majority of mass loss from the ice surface to the terminus. Lake-terminating glaciers are known to experience greater thinning rates and higher velocities than land-terminating glaciers, but the controls on variability in surface elevation change and ice flow between lake-terminating glaciers in different regions remain poorly explored. We combined existing datasets of glacier velocity, surface elevation change and glacial lake area to characterise the evolution of 352 lake-terminating and land-terminating glaciers within three Himalayan sub-regions between 2000 and 2019. These analyses show that the influence of ice-contact lakes propagates up-glacier across only the lowermost 30% of the hypsometric distribution, even where lakes are well established. We find that ice-contact lakes only affect glacier behaviour when the lakes reach an advanced evolutionary stage; most clearly manifested in the Eastern Himalaya by statistically robust differences in glacier-wide surface elevation change between lake-terminating (–0.68 ± 0.05 m a–1) and land-terminating (–0.54 ± 0.04 m a–1) glaciers. These differences are driven by the presence of a greater number of well-developed ice-contact lakes in the Eastern Himalaya compared to in the Western and Central Himalaya, resulting from greater mass loss rates to date

Continuous borehole optical televiewing reveals variable englacial debris concentrations at Khumbu Glacier, Nepal

Surface melting of High Mountain Asian debris-covered glaciers shapes the seasonal water supply to millions of people. This melt is strongly influenced by the spatially variable thickness of the supraglacial debris layer, which is itself partially controlled by englacial debris concentration and melt-out. Here, we present measurements of deep englacial debris concentrations from debris-covered Khumbu Glacier, Nepal, based on four borehole optical televiewer logs, each up to 150 m long. The mean borehole englacial debris content is ≤ 0.7% by volume in the glacier’s mid-to-upper ablation area, and increases to 6.4% by volume near the terminus. These concentrations are higher than those reported for other valley glaciers, although those measurements relate to discrete samples while our approach yields a continuous depth profile. The vertical distribution of englacial debris increases with depth, but is also highly variable, which will complicate predictions of future rates of surface melt and debris exhumation at such glaciers

Hydrology of debris-covered glaciers in High Mountain Asia

The hydrological characteristics of debris-covered glaciers are known to be fundamentally different from those of clean-ice glaciers, even within the same climatological, geological, and geomorphological setting. Understanding how these characteristics influence the timing and magnitude of meltwater discharge is particularly important for regions where downstream communities rely on this resource for sanitation, irrigation, and hydropower, such as High Mountain Asia. The hydrology of debris-covered glaciers is complex: rugged surface topographies typically route meltwater through compound supraglacial-englacial systems involving both channels and ponds, as well as pathways that remain unknown. Low-gradient tongues that extend several kilometres retard water conveyance and promote englacial storage. Englacial conduits are frequently abandoned and reactivated as water supply changes, new lines of permeability are exploited, and drainage is captured due to high rates of surface and subsurface change. Seasonal influences, such as the monsoon, are superimposed on these distinctive characteristics, reorganising surface and subsurface drainage rapidly from one season to the next. Recent advances in understanding have mostly come from studies aimed at quantifying and describing supraglacial processes; little is known about the subsurface hydrology, particularly the nature (or even existence) of subglacial drainage. In this review, we consider in turn the supraglacial, englacial, subglacial, and proglacial hydrological domains of debris-covered glaciers in High Mountain Asia. We summarise different lines of evidence to establish the current state of knowledge and, in doing so, identify major knowledge gaps. Finally, we use this information to suggest six themes for future hydrological research at High Mountain Asian debris-covered glaciers in order to make timely long-term predictions of changes in the water they supply

Multiannual observations and modelling of seasonal thermal profiles through supraglacial debris in the Central Himalaya

Many glaciers in the Central Himalaya are covered with rock debris that modifies the transfer of heat from the atmosphere to the underlying ice. These debris-covered glaciers are experiencing rapid mass loss at rates that have accelerated during the last two decades. Quantifying recent and future glacier mass change requires understanding the relationship between debris thickness and ablation particularly through the summer monsoon season. We present air, near-surface and debris temperatures measured during three monsoon seasons at five sites on Khumbu Glacier in Nepal, and compare these results to similar measurements from two other debris-covered glaciers in this region. Seasonal debris temperature profiles are approximately linear and consistent between sites for thick (>?0.5?m) and thin (<?0.5?m) debris across thicknesses ranging from 0.26 to 2.0?m. The similarities between these multiannual data imply that they are representative of supraglacial debris layers in the monsoon-influenced Himalaya more generally. We compare three methods to calculate sub-debris ablation, including using our temperature measurements with a thermal diffusion model that incorporates a simplified treatment of debris moisture. Estimated ablation between 3 June and 11 October at around 5000?m above sea level ranged from 0.10?m water equivalent beneath 1.5?m of debris to 0.47?m water equivalent beneath 0.3?m debris. However, these values are small when compared to remotely observed rates of surface lowering, suggesting that mass loss from these debris-covered glaciers is greatly enhanced by supraglacial and englacial processes that locally amplify ablationauthorsversionPeer reviewe

The Role of Differential Ablation and Dynamic Detachment in Driving Accelerating Mass Loss From a Debris-Covered Himalayan Glacier

Sustained mass loss from Himalayan glaciers is causing supraglacial debris to expand and thicken, with the expectation that thicker debris will suppress ablation and extend glacier longevity. However, debris-covered glaciers are losing mass at similar rates to clean-ice glaciers in High Mountain Asia. This rapid mass loss is attributed to the combined effects of; (a) low or reversed mass balance gradients across debris-covered glacier tongues, (b) differential ablation processes that locally enhance ablation within the debris-covered section of the glacier, for example, at ice cliffs and supraglacial ponds, and (c) a decrease in ice flux from the accumulation area in response to climatic warming. Adding meter-scale spatial variations in supraglacial debris thickness to an ice-flow model of Khumbu Glacier, Nepal, increased mass loss by 47% relative to simulations assuming a continuous debris layer over a 31-year period (1984–2015 CE) but overestimated the reduction in ice flux. Therefore, we investigated if simulating the effects of dynamic detachment of the upper active glacier from the debris-covered tongue would give a better representation of glacier behavior, as suggested by observations of change in glacier dynamics and structure indicating that this process occurred during the last 100 years. Observed glacier change was reproduced more reliably in simulations of the active, rather than entire, glacier extent, indicating that Khumbu Glacier has passed a dynamic tipping point by dynamically detaching from the heavily debris-covered tongue that contains 20% of the former ice volume

Be‐10 dating of ice‐marginal moraines in the Khumbu Valley, Nepal, Central Himalaya, reveals the response of monsoon‐influenced glaciers to Holocene climate change

The dynamic response of large mountain glaciers to climatic forcing operates over timescales of several centuries and therefore understanding how these glaciers change requires observations of their behavior through the Holocene. We used Be-10 exposure-age dating and geomorphological mapping to constrain the evolution of glaciers in the Khumbu Valley in the Everest region of Nepal. Khumbu and Lobuche Glaciers are surrounded by high-relief lateral and terminal moraines from which seven glacial stages were identified and dated to 7.4 ± 0.2, 5.0 ± 0.3, 3.9 ± 0.1, 2.8 ± 0.2, 1.3 ± 0.1, 0.9 ± 0.02, and 0.6 ± 0.16 ka. These stages correlate to each of the seven latest Holocene regional glacial stages identified across the monsoon-influenced Himalaya, demonstrating that a coherent record of high elevation terrestrial palaeoclimate change can be extracted from dynamic mountain landscapes. The time-constrained moraine complex represents a catchment-wide denudation rate of 0.8–1.4 mm a−1 over the last 8 kyr. The geometry of the ablation area of Khumbu Glacier changed around 4 ka from a broad, shallow ice tongue to become narrower and thicker as restricted by the topographic barrier of the terminal moraine complex

Field-based research directions for investigating the interior of high-elevation debris-covered glaciers

The debris that covers the ablation areas of high-elevation debris-covered glaciers contributes to the distinctive features and processes occurring both on and within such glaciers. Despite recent advances, knowledge of the subsurface environments of high-elevation debris-covered glaciers is still extremely limited. In particular, targeted field-based data are needed to parameterise and refine the projections of these glaciers in numerical models. Here, we outline the current understanding of the internal properties of high-elevation debris-covered glaciers based on direct field-based methods and suggest future research directions for field-based studies

Quantifying ice cliff evolution with multi-temporal point clouds on the debris-covered Khumbu Glacier, Nepal

Measurements of glacier ice cliff evolution are sparse, but where they do exist, they indicate that such areas of exposed ice contribute a disproportionate amount of melt to the glacier ablation budget. We used Structure from Motion photogrammetry with Multi-View Stereo to derive 3-D point clouds for nine ice cliffs on Khumbu Glacier, Nepal (in November 2015, May 2016 and October 2016). By differencing these clouds, we could quantify the magnitude, seasonality and spatial variability of ice cliff retreat. Mean retreat rates of 0.30–1.49 cm d−1 were observed during the winter interval (November 2015–May 2016) and 0.74–5.18 cm d−1 were observed during the summer (May 2016–October 2016). Four ice cliffs, which all featured supraglacial ponds, persisted over the full study period. In contrast, ice cliffs without a pond or with a steep back-slope degraded over the same period. The rate of thermo-erosional undercutting was over double that of subaerial retreat. Overall, 3-D topographic differencing allowed an improved process-based understanding of cliff evolution and cliff-pond coupling, which will become increasingly important for monitoring and modelling the evolution of thinning debris-covered glaciers