Sub-debris Melt Research Articles

Using thermal UAV imagery to model distributed debris thicknesses and sub-debris melt rates on debris-covered glaciers

AbstractSupraglacial debris cover regulates the melt rates of many glaciers in mountainous regions around the world, thereby modifying the availability and quality of downstream water resources. However, the influence of supraglacial debris is often poorly represented within glaciological models, due to the absence of a technique to provide high-precision, spatially continuous measurements of debris thickness. Here, we use high-resolution UAV-derived thermal imagery, in conjunction with local meteorological data, visible UAV imagery and vertically profiled debris temperature time series, to model the spatially distributed debris thickness across a portion of Llaca Glacier in the Cordillera Blanca of Peru. Based on our results, we simulate daily sub-debris melt rates over a 3-month period during 2019. We demonstrate that, by effectively calibrating the radiometric thermal imagery and accounting for temporal and spatial variations in meteorological variables during UAV surveys, thermal UAV data can be used to more precisely represent the highly heterogeneous patterns of debris thickness and sub-debris melt on debris-covered glaciers. Additionally, our results indicate a mean sub-debris melt rate nearly three times greater than the mean melt rate simulated from satellite-derived debris thicknesses, emphasising the importance of acquiring further high-precision debris thickness data for the purposes of investigating glacier-scale melt processes, calibrating regional melt models and improving the accuracy of runoff predictions.

Using ground-based thermal imagery to estimate debris thickness over glacial ice: fieldwork considerations to improve the effectiveness

AbstractDebris-covered glaciers are an important component of the mountain cryosphere and influence the hydrological contribution of glacierized basins to downstream rivers. This study examines the potential to make estimates of debris thickness, a critical variable to calculate the sub-debris melt, using ground-based thermal infrared radiometry (TIR) images. Over four days in August 2019, a ground-based, time-lapse TIR digital imaging radiometer recorded sequential thermal imagery of a debris-covered region of Peyto Glacier, Canadian Rockies, in conjunction with 44 manual excavations of debris thickness ranging from 10 to 110 cm, and concurrent meteorological observations. Inferring the correlation between measured debris thickness and TIR surface temperature as a base, the effectiveness of linear and exponential regression models for debris thickness estimation from surface temperature was explored. Optimal model performance (R2 of 0.7, RMSE of 10.3 cm) was obtained with a linear model applied to measurements taken on clear nights just before sunrise, but strong model performances were also obtained under complete cloud cover during daytime or nighttime with an exponential model. This work presents insights into the use of surface temperature and TIR observations to estimate debris thickness and gain knowledge of the state of debris-covered glacial ice and its potential hydrological contribution.

Controls on the relative melt rates of debris-covered glacier surfaces

Supraglacial debris covers 7% of mountain glacier area globally and generally reduces glacier surface melt. Enhanced energy absorption at ice cliffs and supraglacial ponds scattered across the debris surface leads these features to contribute disproportionately to glacier-wide ablation. However, the degree to which cliffs and ponds actually increase melt rates remains unclear, as these features have only been studied in a detailed manner for selected locations, almost exclusively in High Mountain Asia. In this study we model the surface energy balance for debris-covered ice, ice cliffs, and supraglacial ponds with a set of automatic weather station records representing the global prevalence of debris-covered glacier ice. We generate 5000 random sets of values for physical parameters using probability distributions derived from literature, which we use to investigate relative melt rates and to isolate the melt responses of debris, cliffs and ponds to the site-specific meteorological forcing. Modelled sub-debris melt rates are primarily controlled by debris thickness and thermal conductivity. At a reference thickness of 0.1 m, sub-debris melt rates vary considerably, differing by up to a factor of four between sites, mainly attributable to air temperature differences. We find that melt rates for ice cliffs are consistently 2–3× the melt rate for clean glacier ice, but this melt enhancement decays with increasing clean ice melt rates. Energy absorption at supraglacial ponds is dominated by latent heat exchange and is therefore highly sensitive to wind speed and relative humidity, but is generally less than for clean ice. Our results provide reference melt enhancement factors for melt modelling of debris-covered glacier sites, globally, while highlighting the need for direct measurement of debris-covered glacier surface characteristics, physical parameters, and local meteorological conditions at a variety of sites around the world.

What Can Thermal Imagery Tell Us About Glacier Melt Below Rock Debris?

Rock debris on the surface of a glacier can dramatically reduce the local melt rate, where the primary factor governing melt reduction is debris layer thickness. Relating surface temperature to debris thickness is a recurring approach in the literature, yet demonstrations of reproducibility have been limited. Here, I present the results of a field experiment conducted on the Canwell Glacier, Alaska, United States to constrain how thermal data can be used in glaciology. These datasets include, 1) a measured sub-daily “Østrem curve” time-series; 2) a time-series of high resolution thermal images capturing several segments of different debris thicknesses including the measurements from 1); 3) a thermal profile through a 38 cm debris cover; and 4) two Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite thermal images acquired within 2 and 3 min of a field-based thermal camera image. I show that, while clear sky conditions are when space-borne thermal sensors can image a glacier, this is an unfavorable time, limiting the likelihood that different thicknesses of debris will have a unique thermal signature. I then propose an empirical approach to estimate debris thickness and compare it to two recently published methods. I demonstrate that instantaneous calibration is essential in the previously published methods, where model parameters calibrated only 1 h prior to a repeat thermal image return diminished debris thickness estimates, while the method proposed here remains robust through time and does not appear to require re-calibration. I then propose a method that uses a time-series of surface temperature at one location and debris thickness to estimate bare-ice and sub-debris melt. Results show comparable cumulative melt estimates to a recently published method that requires an explicit/external estimate of bare ice melt. Finally, I show that sub-pixel corrections to ASTER thermal imagery can enable a close resemblance to high resolution, field-based thermal imagery. These results offer a deeper insight into what thermal data can and cannot tell us about surface debris properties and glacier melt.

Quantifying heterogeneous monsoonal melt on a debris-covered glacier in Nepal Himalaya using repeat uncrewed aerial system (UAS) photogrammetry

AbstractThe ablation zones of debris-covered glaciers in Himalaya exhibit heterogeneous processes and melt patterns. Although sub-debris melt is measured at ablation stakes, the high variability of debris thickness necessitates distributed melt measurements at the glacier scale. Focusing on Annapurna III Glacier, we used uncrewed aerial system (UAS) photogrammetry to estimate total volume loss and slope-perpendicular glacier melt between May and November 2019 using flow-corrected point clouds. Results indicated the average elevation change was −1.10 ± 0.19 m, while the mean melt was −0.87 m w.e., equating to a mean melt rate of −0.47 cm w.e. d−1. However, the spatial pattern was highly variable due to complex local processes necessitating future study over short intervals. The evaluation of specific areas showed the interplay of debris thickness variability, subseasonal debris redistribution, supraglacial channel reconfiguration and the imprint of relict ice cliffs in leading to contemporary melt rates. Ice cliffs had higher melt distances (mean −3.9 ± 0.19 m) compared to non-cliff areas (mean −0.75 ± 0.19 m) and were the predominant control on the spatial patterns of seasonal melt rates. Crucially, the definition of ice cliff areas from thinning data has a profound impact on derived melt rates and melt enhancement. Our study demonstrates the possibility and utility of deriving fully-distributed slope-perpendicular melt measurements.

Distributed Global Debris Thickness Estimates Reveal Debris Significantly Impacts Glacier Mass Balance.

Supraglacial debris affects glacier mass balance as a thin layer enhances surface melting, while a thick layer reduces it. While many glaciers are debris‐covered, global glacier models do not account for debris because its thickness is unknown. We provide the first globally distributed debris thickness estimates using a novel approach combining sub‐debris melt and surface temperature inversion methods. Results are evaluated against observations from 22 glaciers. We find the median global debris thickness is ∼0.15 ± 0.06 m. In all regions, the net effect of accounting for debris is a reduction in sub‐debris melt, on average, by 37%, which can impact regional mass balance by up to 0.40 m water equivalent (w.e.) yr‐1. We also find recent observations of similar thinning rates over debris‐covered and clean ice glacier tongues is primarily due to differences in ice dynamics. Our results demonstrate the importance of accounting for debris in glacier modeling efforts.

Morphometric evolution of Everest region debris-covered glaciers

Debris-covered glaciers in the central Himalaya have now experienced several decades of sustained ice loss, manifested predominantly in glacier surface lowering. In particular, glacier surfaces of low longitudinal gradient and low ice surface velocity have developed locally complex surface topographies and undergone profound changes in supraglacial hydrology. In this study we examine the development of complex ice surface topography across six debris-covered glaciers in the Everest region over the last four decades via a new metric of glacier surface relief applied to Digital Elevation Models (DEMs). We focus in on Khumbu Glacier, and use fine spatial and temporal resolution DEMs covering a period of 28 months to quantify the contemporary contribution of ice cliff and supraglacial pond expansion to overall mass loss from stagnant areas of ice. On the broader scale, we find three common long-term changes in glacier surface topography, (1) glacier-wide expansion of high relief topography in response to ice cliff and supraglacial pond network evolution, (2) up-glacier expansion of high local relief zones that may be caused by differential sub-debris melt beneath thin debris, and (3) increase in glacier surface relief proximal to glacier termini caused by supraglacial stream incision where linked proglacial-supraglacial hydrological networks exist. Overall, we contend that these topographic measurements will be important for understanding glacier surface water storage and also the energy balance of a debris-covered glacier surface, both of which could exacerbate future ice loss and downstream meltwater supply.

Incorporating moisture content in surface energy balance modeling of a debris-covered glacier

Abstract. Few surface energy balance models for debris-covered glaciers account for the presence of moisture in the debris, which invariably affects the debris layer's thermal properties and, in turn, the surface energy balance and sub-debris melt of a debris-covered glacier. We adapted the interactions between soil, biosphere, and atmosphere (ISBA) land surface model within the SURFace EXternalisée (SURFEX) platform to represent glacier debris rather than soil (referred to hereafter as ISBA-DEB). The new ISBA-DEB model includes the varying content, transport, and state of moisture in debris with depth and through time. It robustly simulates not only the thermal evolution of the glacier–debris–snow column but also moisture transport and phase changes within the debris – and how these, in turn, affect conductive and latent heat fluxes. We discuss the key developments in the adapted ISBA-DEB and demonstrate the capabilities of the model, including how the time- and depth-varying thermal conductivity and specific heat capacity depend on evolving temperature and moisture. Sensitivity tests emphasize the importance of accurately constraining the roughness lengths and surface slope. Emissivity, in comparison to other tested parameters, has less of an effect on melt. ISBA-DEB builds on existing work to represent the energy balance of a supraglacial debris layer through time in its novel application of a land surface model to debris-covered glaciers. Comparison of measured and simulated debris temperatures suggests that ISBA-DEB includes some – but not all – processes relevant to melt under highly permeable debris. Future work, informed by further observations, should explore the importance of advection and vapor transfer in the energy balance.

Quantifying Debris Thickness of Debris‐Covered Glaciers in the Everest Region of Nepal Through Inversion of a Subdebris Melt Model

AbstractDebris‐covered glaciers are ubiquitous in the Himalaya, and supraglacial debris significantly alters how glaciers respond to climate forcing. Estimating debris thickness at the glacier scale, however, remains a challenge. This study inverts a subdebris melt model to estimate debris thickness for three glaciers in the Everest region from digital elevation model difference‐derived elevation change. Flux divergences are estimated from ice thickness and surface velocity data. Monte Carlo simulations are used to incorporate the uncertainties associated with debris properties, flux divergence, and elevation change. On Ngozumpa Glacier, surface lowering data from 2010 to 2012 and 2012 to 2014 are used to calibrate and validate the method, respectively. The debris thickness estimates are consistent with existing in situ measurements. The method performs well over both actively flowing and stagnant parts of the glacier and is able to accurately estimate thicker debris (>0.5 m). Uncertainties associated with the thermal conductivity and elevation change contribute the most to uncertainties of the debris thickness estimates. The surface lowering associated with ice cliffs and supraglacial ponds was found to significantly reduce debris thickness, especially for thicker debris. The method is also applied to Khumbu and Imja‐Lhotse Shar Glaciers to highlight its potential for regional application.

Debris thickness patterns on debris-covered glaciers

Many debris-covered glaciers have broadly similar debris thickness patterns: surface debris thickens and tends to transition from convex- to concave-up-down glacier. We explain this pattern using theory (analytical and numerical models) paired with empirical observations. Down glacier debris thickening results from the conveyor-belt-like nature of the glacier surface in the ablation zone (debris can typically only be added but not removed) and from the inevitable decline in ice surface velocity toward the terminus. Down-glacier thickening of debris leads to the reduction of sub-debris melt and debris emergence toward the terminus. Convex-up debris thickness patterns occur near the up-glacier end of debris covers where debris emergence dominates (ablation controlled). Concave-up debris thickness patterns occur toward glacier termini where declining surface velocities dominate (velocity controlled). A convex-concave debris thickness profile inevitably results from the transition between ablation-control and velocity-control down-glacier. Debris thickness patterns deviating from this longitudinal shape are most likely caused by changes in hillslope debris supply through time. By establishing this expected debris thickness pattern, the effects of climate change on debris cover can be better identified.

Impact of debris cover on glacier ablation and atmosphere–glacier feedbacks in the Karakoram

Abstract. The Karakoram range of the Hindu-Kush Himalaya is characterized by both extensive glaciation and a widespread prevalence of surficial debris cover on the glaciers. Surface debris exerts a strong control on glacier surface-energy and mass fluxes and, by modifying surface boundary conditions, has the potential to alter atmosphere–glacier feedbacks. To date, the influence of debris on Karakoram glaciers has only been directly assessed by a small number of glaciological measurements over short periods. Here, we include supraglacial debris in a high-resolution, interactively coupled atmosphere–glacier modeling system. To investigate glaciological and meteorological changes that arise due to the presence of debris, we perform two simulations using the coupled model from 1 May to 1 October 2004: one that treats all glacier surfaces as debris-free and one that introduces a simplified specification for the debris thickness. The basin-averaged impact of debris is a reduction in ablation of ~ 14 %, although the difference exceeds 5 m w.e. on the lowest-altitude glacier tongues. The relatively modest reduction in basin-mean mass loss results in part from non-negligible sub-debris melt rates under thicker covers and from compensating increases in melt under thinner debris, and may help to explain the lack of distinct differences in recent elevation changes between clean and debris-covered ice. The presence of debris also strongly alters the surface boundary condition and thus heat exchanges with the atmosphere; near-surface meteorological fields at lower elevations and their vertical gradients; and the atmospheric boundary layer development. These findings are relevant for glacio-hydrological studies on debris-covered glaciers and contribute towards an improved understanding of glacier behavior in the Karakoram.

Modeling runoff from the inylchek glaciers and filling of ice‐dammed lake merzbacher, central tian shan

We developed a spatially distributed conceptual melt and runoff model able to cope with debris covered glacier parts. Subdebris melt was calculated based on empirical equations derived from field observations. The model was applied at the nylchek glaciers. While runoff from outhern nylchek flows directly into nylchek iver, runoff from orthern nylchek is blocked by outhern nylchek and forms ake erzbacher. This lake is drained frequently by outburst floods. To calibrate the model, a multi‐algorithm, genetically adaptive multi‐objective method () was used. Daily runoff observations from nylchek iver as well as glacier mass balance gradients from other glaciers in the ian han were used as objective functions in three split samples. Runoff simulations yielded ash–utcliffe coefficients mostly higher than 0.9, and the mass balance gradients were simulated in an acceptable range. The results of the different split samples vary considerably. Considering the used objective functions and additional data such as the equilibrium line altitude, we chose the most promising parameter set to calculate the filling of ake erzbacher. This analysis indicated answers to several questions concerning the hydrological system of the lake. Apart from calving at the ice dam formed by outhern nylchek glacier, the lake seemed so be filled predominantly by meltwater from orthern nylchek glacier. The lake has a permanent runoff, also during the filling process. After the outburst, the channels seemed to be closed quickly.

A distributed energy-balance melt model of an alpine debris-covered glacier

AbstractDistributed energy-balance melt models have rarely been applied to glaciers with extensive supraglacial debris cover. This paper describes the development of a distributed melt model and its application to the debris-covered Miage glacier, western Italian Alps, over two summer seasons. Sub-debris melt rates are calculated using an existing debris energy-balance model (DEB-Model), and melt rates for clean ice, snow and partially debris-covered ice are calculated using standard energy-balance equations. Simulated sub-debris melt rates compare well to ablation stake observations. Melt rates are highest, and most sensitive to air temperature, on areas of dirty, crevassed ice on the middle glacier. Here melt rates are highly spatially variable because the debris thickness and surface type varies markedly. Melt rates are lowest, and least sensitive to air temperature, beneath the thickest debris on the lower glacier. Debris delays and attenuates the melt signal compared to clean ice, with peak melt occurring later in the day with increasing debris thickness. The continuously debris-covered zone consistently provides 30% of total melt throughout the ablation season, with the proportion increasing during cold weather. Sensitivity experiments show that an increase in debris thickness of 0.035 m would offset 18C of atmospheric warming.

Distribution of debris thickness and its effect on ice melt at Hailuogou glacier, southeastern Tibetan Plateau, using in situ surveys and ASTER imagery

AbstractDebris cover is widely present in glacier ablation areas of the Tibetan Plateau, and its spatial distribution greatly affects glacier melt rates. High-resolution in situ measurements of debris thickness on Hailuogou glacier, Mount Gongga, southeastern Tibetan Plateau, show pronounced inhomogeneous debris distribution. An analysis of transverse and longitudinal profiles indicates that the ground-surveyed debris thicknesses and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER)-derived thermal resistances of debris layers correlate strongly over the entire ablation area. Across- and along-glacier patterns of ASTER-derived thermal resistance correspond well with spatial patterns of debris thickness, which may reflect large-scale variations in the extent and thickness of the debris cover. The ice melt rate variability over the ablation area simulated by a surface energy-balance model that considered thermal resistance of the debris layer indicates clearly the crucial role of debris and its spatial continuity in modifying the spatial characteristics of melt rates. Because of the inhomogeneous distribution of debris thickness, about 67% of the ablation area on Hailuogou glacier has undergone accelerated melting, whereas about 19% of the ablation area has experienced inhibited melting, and the sub-debris melt rate equals the bare-ice melt rate in only 14% of the ablation area.

Glaciological results of the 2005 expedition to Inylchek Glacier, Central Tian Shan

Like many other glaciers in Central Asia, Southern Inylchek glacier in the Kyrgyz Tian Shan is covered by supraglacial moraine, which drastically influences melt rates and complicates the estimation of ablation. The quantification of sub-debris melt from simple parameters is still an unsolved problem, but also essential to predict future yield from high mountains. Snow cover and glacier ice are the main water storages for the surrounding arid lowlands and a better understanding of ablation processes is the prerequisite for a sustainable water resources management. Another interesting feature of Southern Inylchek glacier is the existence of an ice dammed lake in a tributary valley, which is drained regularly by outburst floods. Improvements in predicting these floods would lower the risk potential for the downstream population. The main objectives of a group of glaciologists which participated in an expedition to the glacier in 2005 were to investigate melt rates on debris covered glacier parts and to quantify the ice flux into the glacier lake. The results of their field experiments are reported in this paper.