In areas of high relief, many glaciers have extensive covers of supraglacial debris in their ablation zones, which alters both rates and spatial patterns of melting, with important consequences for glacier response to climate change. Wastage of debris-covered glaciers can be associated with the formation of large moraine-dammed lakes, posing risk of glacier lake outburst floods (GLOFs). In this paper, we use observations of glaciers in the Mount Everest region to present an integrated view of debris-covered glacier response to climate change, which helps provide a long-term perspective on evolving GLOF risks. In recent decades, debris-covered glaciers in the Everest region have been losing mass at a mean rate of 0.32 m yr¹, although in most cases there has been little or no change in terminus position. Mass loss occurs by 4 main processes: (1) melting of clean ice close to glacier ELAs; (2) melting beneath surface debris; (3) melting of ice cliffs and calving around the margins of supraglacial ponds; and (4) calving into deep proglacial lakes. Modelling of processes (1) and (2) shows that Everest-region glaciers typically have an inverted ablation gradient in their lower reaches, due to the effects of a down-glacier increase in debris thickness. Mass loss is therefore focused in the mid parts of glacier ablation zones, causing localised surface lowering and a reduction in downglacier surface gradient, which in turn reduce driving stress and glacier velocity, so the lower ablation zones of many glaciers are now stagnant. Model results also indicate that increased summer temperatures have raised the altitude of the rain-snow transition during the summer monsoon period, reducing snow accumulation and ice flux to lower elevations. As downwasting proceeds, formerly efficient supraglacial and englacial drainage networks are broken up, and supraglacial lakes form in hollows on the glacier surface. Ablation rates around supraglacial lakes are typically one or two orders of magnitude greater than sub-debris melt rates, so extensive lake formation accelerates overall rates of ice loss. Most supraglacial lakes are 'perched' above hydrological base level, and are susceptible to drainage if they become connected to the englacial drainage system. Speleological surveys of conduits show that large englacial voids can be created by drainage of warm lake waters along pre-existing weaknesses in the ice. Roof collapses can open these voids up to the surface, and commonly provide the nuclei of new lakes. Thus, by influencing both lake drainage and formation, englacial conduits exert a strong control on surface ablation rates. An important threshold is crossed when downwasting glacier surfaces intersect the hydrological base level of the glacier. Base-level lakes formed behind intact moraine dams can grow monotonically, and in some cases can pose serious GLOF hazards. Glacier termini can evolve in different ways in response to the same climatic forcing, so that potentially hazardous lakes will form in some situations but no...

Seasonal acceleration of the Greenland Ice Sheet is influenced by the dynamic response of the subglacial hydrologic system to variability in meltwater delivery to the bed via crevasses and moulins (vertical conduits connecting supraglacial water to the bed of the ice sheet). As the melt season progresses, the subglacial hydrologic system drains supraglacial meltwater more efficiently, decreasing basal water pressure and moderating the ice velocity response to surface melting. However, limited direct observations of subglacial water pressure mean that the spatiotemporal evolution of the subglacial hydrologic system remains poorly understood. Here we show that ice velocity is well correlated with moulin hydraulic head but is out of phase with that of nearby (0.3-2 kilometres away) boreholes, indicating that moulins connect to an efficient, channelized component of the subglacial hydrologic system, which exerts the primary control on diurnal and multi-day changes in ice velocity. Our simultaneous measurements of moulin and borehole hydraulic head and ice velocity in the Paakitsoq region of western Greenland show that decreasing trends in ice velocity during the latter part of the melt season cannot be explained by changes in the ability of moulin-connected channels to convey supraglacial melt. Instead, these observations suggest that decreasing late-season ice velocity may be caused by changes in connectivity in unchannelized regions of the subglacial hydrologic system. Understanding this spatiotemporal variability in subglacial pressures is increasingly important because melt-season dynamics affect ice velocity beyond the conclusion of the melt season.

Penetration of surface meltwater to the bed of the Greenland Ice Sheet each summer causes an initial increase in ice speed due to elevated basal water pressure, followed by slowdown in late summer that continues into fall and winter. While this seasonal pattern is commonly explained by an evolution of the subglacial drainage system from an inefficient distributed to efficient channelized configuration, mounting evidence indicates that subglacial channels are unable to explain important aspects of hydrodynamic coupling in late summer and fall. Here we use numerical models of subglacial drainage and ice flow to show that limited, gradual leakage of water and lowering of water pressure in weakly connected regions of the bed can explain the dominant features in late and post melt season ice dynamics. These results suggest that a third weakly connected drainage component should be included in the conceptual model of subglacial hydrology.

On uncrevassed regions of polythermal glaciers, englacial conduits can form by incision of supraglacial stream channels followed by roof closure. The origin and evolution of examples in Longyearbreen, Svalbard, and Khumbu Glacier, Nepal, were determined by speleological survey. The development of perennial incised channels requires that incision is significantly faster than glacier surface ablation, and thus will be favoured by high meltwater discharges in combination with cool climatic conditions or thick debris cover. Incised canyons can become blocked by drifted winter snow, refrozen meltwater, ice rafting from non-local sources (allochthonous breccias) and roof collapses (autochthonous breccias). Conduit closure can also occur in response to ice creep, particularly at depth. Following isolation from the surface, englacial conduits continue to evolve by vadose incision down to local base level. In the case of Longyearbreen, incision allowed the channel to reach the glacier bed, but on Khumbu Glacier deep incision is prevented because an effectively impermeable terminal moraine provides a high base level for the glacier drainage system. During our period of observations, deeper parts of the Longyearbreen conduit became blocked by a combination of ice accumulation and creep, causing the stream course to be re-routed to higher levels. On that glacier, incision, blockage and upward re-routing are cyclic. We conclude that 'cut and closure' is the dominant mechanism of englacial conduit formation on uncrevassed regions of polythermal glaciers.

No abstract