<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Earth Sci.</journal-id>
<journal-title>Frontiers in Earth Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Earth Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-6463</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/feart.2021.652279</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Earth Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Understanding Complex Debris-Covered Glaciers: Concepts, Issues, and Research Directions</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Huo</surname> <given-names>Da</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1086495/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Bishop</surname> <given-names>Michael P.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Bush</surname> <given-names>Andrew B. G.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1253294/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Geography, Texas A&#x00026;M University</institution>, <addr-line>College Station, TX</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Earth and Atmospheric Sciences, University of Alberta</institution>, <addr-line>Edmonton, AB</addr-line>, <country>Canada</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Duncan Joseph Quincey, University of Leeds, United Kingdom</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Leif S. Anderson, University of Lausanne, Switzerland; Christoph Mayer, Bavarian Academy of Sciences and Humanities, Germany</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Da Huo <email>dh21&#x00040;tamu.edu</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Cryospheric Sciences, a section of the journal Frontiers in Earth Science</p></fn></author-notes>
<pub-date pub-type="epub">
<day>24</day>
<month>05</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>9</volume>
<elocation-id>652279</elocation-id>
<history>
<date date-type="received">
<day>12</day>
<month>01</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>04</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2021 Huo, Bishop and Bush.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Huo, Bishop and Bush</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>Understanding the climate-glacier dynamics of debris-covered glaciers is notoriously difficult given a multitude of controlling factors and feedback mechanisms involving climate forcing, debris-load properties, supraglacial water bodies, and multi-scale topographic effects. Recent studies have provided insights into controlling factors, and have reported the presence of anomalies that contradict the general consensus of the protective influence of debris loads on ablation dynamics. Nevertheless, numerous processes that regulate glacier dynamics at various spatial and temporal scales have not been adequately accounted for in statistical and numerical modeling studies. Furthermore, important feedbacks involving ablation, topography, irradiance, gravitational debris flux, and supraglacial ponding are often neglected or oversimplified in existing models, which poses a challenge to our understanding of conflicting field observations such as the accelerated mass loss of many Himalayan glaciers, and glacier-subsystem responses (ice-flow, debris flux, surface morphology, and supraglacial water bodies) to climate forcing. This paper provides insights into the complexity of debris-covered glacier systems by addressing concepts and issues associated with forcing factors and glacial processes, and highlights the importance of understanding system couplings and feedbacks. Specifically, we review recent studies on debris-covered glaciers and utilize simulation results based on the Baltoro Glacier in the central Karakoram to discuss important concepts and issues. Our results demonstrate that climate forcing, the properties and transport of debris, topography and supraglacial water bodies are the key controlling factors in a debris-covered glacier system, and that their coupled effects and positive feedbacks may increase the ice loss of a debris-covered glacier. We also recommend new research directions for future studies.</p></abstract>
<kwd-group>
<kwd>debris-covered glaciers</kwd>
<kwd>climate-glacier dynamics</kwd>
<kwd>debris load</kwd>
<kwd>supraglacial water bodies</kwd>
<kwd>glacier model</kwd>
<kwd>Baltoro glacier</kwd>
</kwd-group>
<counts>
<fig-count count="10"/>
<table-count count="1"/>
<equation-count count="22"/>
<ref-count count="115"/>
<page-count count="21"/>
<word-count count="16617"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1. Introduction</title>
<p>Debris-covered glaciers (DCGs) present a series of unique and complex challenges to those investigating climate-glacier dynamics and their role in downstream impacts such as freshwater availability, hydroelectric power generation, natural hazards and landscape evolution (Bishop et al., <xref ref-type="bibr" rid="B20">2002</xref>; Seong et al., <xref ref-type="bibr" rid="B100">2009</xref>; Immerzeel et al., <xref ref-type="bibr" rid="B58">2010</xref>; Reid and Brock, <xref ref-type="bibr" rid="B88">2010</xref>; Bush et al., <xref ref-type="bibr" rid="B27">2020</xref>). While debris-free glaciers are relatively well-understood and modeled, DCGs introduce a number of additional complex dynamical and thermodynamical feedbacks between the underlying glacier ice, the debris cover, the adjacent-terrain geomorphological conditions, and the overlying atmosphere. The thermal and radiative characteristics of the debris cover itself depend on geographic location and the surrounding lithology. The thickness of the debris cover, which so critically determines whether the debris acts as an insulator for the ice or accelerates ice ablation, depends on the steepness of the surrounding topography, regional erosion characteristics, and glacier dynamics, which in turn depend on atmosphere-glacier interactions that have not yet been fully quantified (Bolch et al., <xref ref-type="bibr" rid="B22">2012</xref>; Dobreva et al., <xref ref-type="bibr" rid="B35">2017</xref>; Farinotti et al., <xref ref-type="bibr" rid="B39">2020</xref>).</p>
<p>Currently, we have not achieved an agreement on the state and fate of many DCGs in High-Mountain Asia due to uncertainties related to climate, topographic evolution, debris-load forcings, and the development of supraglacial water bodies and ice-cliffs. Most researchers agree that debris insulation makes a glacier less sensitive to climate change (e.g., Reid and Brock, <xref ref-type="bibr" rid="B86">2014</xref>; Pratap et al., <xref ref-type="bibr" rid="B80">2015</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Vincent et al., <xref ref-type="bibr" rid="B108">2016</xref>). Recent observations and simulations, however, show that many DCGs exhibit comparable ice-loss to debris-free glaciers despite the presence of supraglacial debris (K&#x000E4;&#x000E4;b et al., <xref ref-type="bibr" rid="B61">2012</xref>; Immerzeel et al., <xref ref-type="bibr" rid="B57">2013</xref>; Fujita et al., <xref ref-type="bibr" rid="B42">2014</xref>; Salerno et al., <xref ref-type="bibr" rid="B97">2017</xref>). Therefore, it appears that DCGs may be more sensitive to climate change than previously thought, as a multitude of processes and feedbacks have not been accurately characterized. Existing DCG models do not provide a satisfactory explanation to these conflicting interpretations of climate-glacier dynamics, which is partially due to the over-simplified assumptions of important processes and feedbacks that involve ablation dynamics, debris flux and property variation, surface ponding, and topographic evolution.</p>
<p>Mountains that exhibit DCGs typically have substantial vertical relief and erosional forces that produce significant mass flux from the sidewalls onto the glacier. In the central Karakoram in Pakistan, extreme relief and complex climate systems generate some of the largest DCGs in the world (Bishop et al., <xref ref-type="bibr" rid="B16">2001</xref>; Hewitt, <xref ref-type="bibr" rid="B52">2005</xref>; Copland et al., <xref ref-type="bibr" rid="B33">2009</xref>). Precipitation rates can be highly variable in space and time caused by the dominance of regional climate systems and topography. A particular system such as the Westerlies or the Monsoon could dominate over an area over a specific time frame, and it is also possible that an area could receive precipitation from both systems for a period. Furthermore, topographic variation partially controls the degree of orographic precipitation and this would also vary depending upon which system dominants over an area (i.e., moisture conditions, wind velocity and direction). This very high degree of coupled climate-topographic interactions presents a challenge to the numerical modeling of climate-glacier dynamics.</p>
<p>The surface-energy balance determines the ablation rates on DCGs, and net long-wave and short-wave fluxes play a major role in that balance, but the atmospheric conditions through which that radiation must pass play an equally important role: if the atmosphere is cloudy then the incoming short- and long wave radiations are fundamentally altered. Additionally, atmospheric temperature above the surface is critical to determining precipitation type (rain, snow, or some combination of frozen/liquid water). For example, debris cover can provide a surface heat source that drives convection and precipitation in the overlying atmosphere (e.g., Collier et al., <xref ref-type="bibr" rid="B30">2015</xref>); depending on the temperature, that precipitation may fall as snow which would accumulate on the glacier ice in the absence of the debris, but given the presence of the warm debris the snow may melt and percolate down to the glacier surface, thereby enhancing ablation.</p>
<p>Supraglacial ponding is another critical process on DCGs. Supraglacial ponds and lakes play an important role in glacier ice-loss as well as glacial hydrology (Fountain and Walder, <xref ref-type="bibr" rid="B41">1998</xref>; Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>; Wessels et al., <xref ref-type="bibr" rid="B109">2002</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>, <xref ref-type="bibr" rid="B72">2018</xref>). Researchers have identified an increasing number of ponds on some large DCGs (Gibson et al., <xref ref-type="bibr" rid="B46">2017</xref>). Studies have also indicated that the increase in spatial density of supraglacial ponds can be facilitated by localized thinning, insufficient drainage, collapse of water channels, and overall lowering that decreases the slope of the glacier altitude profile in the ablation zone (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>). Unfortunately, many of these important mechanisms have not been adequately studied and are often neglected in existing glacier models and simulations.</p>
<p>Such complexities are only now beginning to be understood through the use of numerical models (e.g., Collier et al., <xref ref-type="bibr" rid="B31">2013</xref>; Rowan et al., <xref ref-type="bibr" rid="B94">2015</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Miles et al., <xref ref-type="bibr" rid="B72">2018</xref>). Field observations of en-debris and sub-debris processes can be difficult to obtain given time, economic, and logistic constraints. Nevertheless, field observations are essential for understanding the magnitude of specific properties and processes, providing constraints for input into numerical models, and for validation and accuracy assessment of numerical modeling efforts. However, many mechanisms that govern glacier dynamics can not be fully understood through field observations and remote sensing data. Consequently, numerical modeling that evaluates new parameterization schemes to account for processes and important feedback is sorely needed to better understand complex climate-glacier dynamics of DCGs.</p>
<p>The objectives of this paper are therefore to discuss a multitude of complex interactions and uncertainties that limit our understanding of how DCGs may be responding to climate change. We review our current understanding of debris-covered glaciers and identify gaps in knowledge about some processes and feedbacks that appear to play important roles with respect to ablation dynamics. We also utilize numerical simulations to demonstrate various concepts and processes. We specifically discuss various issues and use our simulation results to identify future research pathways that can help eliminate uncertainty and promote new empirical and numerical-modeling research. We acknowledge that we cannot account for all processes, nor utilize fully comprehensive parameterization schemes that accurately predict process rates (an objective of numerical modeling research). Rather, we attempt to provide for a first-order characterization and understanding of processes and feedbacks based upon the simulations that we have incorporated into our discussion. Furthermore, it should be noted that this work should not be viewed as a formal modeling paper, complete with parameter sensitivity and uncertainty analyses, as this is beyond the scope of our objectives.</p>
<p>We synthesize our treatment of this complex topic by organizing the paper based upon concepts related to climate-glacier system components. We first address climate forcing and specifically address radiative and precipitation forcing, where we identify some processes that should be considered regarding the surface energy balance in high relief mountains. We then discuss debris load properties and debris transport dynamics. We finally address various feedback mechanisms and system couplings on the glacier surface with respect to ablation, supraglacial water bodies, ice cliffs, and surface topography.</p>
</sec>
<sec id="s2">
<title>2. Climate Forcing</title>
<p>Radiative and precipitation forcing govern glacier dynamics, although the partitioning of these forcing components in any particular geographic area is not well-known with any degree of certainty, given a multitude of multi-scale topographic effects and our inability to accurately characterize and account for radiative and precipitation processes (Bishop et al., <xref ref-type="bibr" rid="B21">2019</xref>; Bush et al., <xref ref-type="bibr" rid="B27">2020</xref>). For example, it is thought that radiative forcing dominates in the Eastern Himalaya, while precipitation forcing has the biggest influence in the western Himalaya (Bolch et al., <xref ref-type="bibr" rid="B22">2012</xref>). Nevertheless, we require a much more detailed understanding of climate-glacier dynamics from an internal and external scale-dependent forcing perspective, as various aspects of climate, glacial, and geomorphological systems need to be accounted for. Unfortunately, we lack basic data regarding radiation transfer and precipitation components that govern glacier states, and we do not know the responses and sensitivity to climate change for many regions of the world (Bolch et al., <xref ref-type="bibr" rid="B22">2012</xref>; Dobreva et al., <xref ref-type="bibr" rid="B35">2017</xref>; Bishop et al., <xref ref-type="bibr" rid="B21">2019</xref>; Bush et al., <xref ref-type="bibr" rid="B27">2020</xref>; Farinotti et al., <xref ref-type="bibr" rid="B39">2020</xref>). Furthermore, the data sets that we tend to rely upon for such information (e.g., climate simulations and re-analysis data) do not accurately characterize local-scale radiative and precipitation fluxes due to inadequate representation of topography and/or processes (Bush et al., <xref ref-type="bibr" rid="B27">2020</xref>).</p>
<p>Energy balance calculation at the glacier surface (with or without debris) is the first step toward determining the low-level air temperature above the glacier, but it neglects potential temperature advection by atmospheric circulations (e.g., anabatic, katabatic, or frontal winds, which all advect temperature). Furthermore, during the daytime, heating of the debris cover itself can generate atmospheric convection above the debris leading to strong converging surface winds and potential precipitation, which can further alter the surface energy balance. Consequently, numerous radiative-transfer (RT) parameters, circulation dynamics and cloud properties, and convective, cyclonic and orographic precipitation mechanisms must be considered to understand glacier response and sensitivity to forcings.</p>
<sec>
<title>2.1. Irradiant Fluxes</title>
<p>The short-wave and long-wave radiant energy components of the surface-energy balance strongly regulate ablation dynamics. The short-wave surface irradiance (<italic>E</italic>) components include the direct irradiance from the sun (<italic>E</italic><sub><italic>b</italic></sub>), the diffuse-skylight irradiance due to atmospheric scattering (<italic>E</italic><sub><italic>d</italic></sub>) and the adjacent-terrain irradiance (<italic>E</italic><sub><italic>t</italic></sub>). All of these components are space-time-topography-wavelength dependent, and numerical modeling is the only way to account for the inherent scale-dependent operation of orbital, atmospheric and lithospheric processes that govern the variability of the magnitude of these RT parameters over glaciers (Bishop et al., <xref ref-type="bibr" rid="B21">2019</xref>). These radiation fluxes significantly contribute to ablation at a particular location depending upon the surrounding land cover and topographic conditions. Studies have estimated that the shortwave radiation usually accounts for 75% or more of total irradiance on glacier surfaces (Oerlemans and Klok, <xref ref-type="bibr" rid="B76">2002</xref>).</p>
<sec>
<title>2.1.1. Short-Wave Direct Irradiance</title>
<p>The direct irradiance is the dominant short-wave irradiance component (Proy et al., <xref ref-type="bibr" rid="B82">1989</xref>; Olson and Rupper, <xref ref-type="bibr" rid="B77">2019</xref>). It is strongly controlled by local and meso-scale topographic effects. It can be represented as:</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>E</mml:mi></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:msubsup></mml:mstyle><mml:msup><mml:mrow><mml:mi>E</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02193;</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo class="qopname">cos</mml:mo><mml:mi>i</mml:mi><mml:mi>S</mml:mi><mml:mi>d</mml:mi><mml:mo>&#x003BB;</mml:mo><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003BB;<sub><italic>l</italic></sub> is wavelength, <italic>E</italic><sup>0</sup> is the exo-atmospheric irradiance adjusted for Earth-Sun distance, <italic>T</italic><sup>&#x02193;</sup> is the downward atmospheric transmittance of radiation governed by atmospheric absorption and scattering, and &#x003B8;<sub><italic>s</italic></sub> is the solar zenith angle. Transmittance also accounts for the atmospheric optical depth that is partially governed by the relief structure of the landscape (i.e., altitude).</p>
<p>The cosine of the incidence angle (cos<italic>i</italic>) characterizes the local topographic effects on irradiance given solar and terrain geometry relationships. Specifically, the incidence angle of illumination (<italic>i</italic>) between the sun and normal to the ground surface is defined as:</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mo class="qopname">cos</mml:mo><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mo class="qopname">cos</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo class="qopname">cos</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mo class="qopname">sin</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo class="qopname">sin</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo class="qopname">cos</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003B8;<sub><italic>s</italic></sub> is the apparent solar zenith angle that accounts for altitude variations and atmospheric refraction given the atmospheric temperature and pressure profiles, &#x003B8;<sub><italic>t</italic></sub> represents the terrain slope angle, &#x003D5;<sub><italic>s</italic></sub> is the solar azimuth angle, and &#x003D5;<sub><italic>t</italic></sub> is the terrain slope-azimuth angle. Calculation of cos<italic>i</italic> is possible with the use of a digital elevation model, and negative values must be corrected to 0.0. Local topographic properties can have a significant effect on the magnitude of the direct irradiance that reaches the surface of a glacier. Ablation and meltwater transport, and gravitational sediment fluxes all alter the altitude, slope and slope azimuth of the ice and sediment surfaces of a glacier. A feedback between ablation-topography and surface irradiance controls the topographic evolution of a glacier surface (Huo et al., <xref ref-type="bibr" rid="B54">2020</xref>). The importance of the local topographic conditions on the direct irradiance cannot be overlooked. The cosine of the incidence angle can range from 0 to 1 based upon slope angle and slope-azimuth angle variations caused by sediment transport and ablation processes. Based upon our simulations in the Karakoram (see <xref ref-type="fig" rid="F1">Figure 1</xref>), direct irradiance values can be greater than 1,000 <italic>Wm</italic><sup>&#x02212;2</sup>. Consequently, rapidly evolving topography should be accounted for as this sub-parameter strongly regulates ablation, all other factors such as debris properties being equal.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Simulated short-wave direct-beam irradiance (<italic>E</italic><sub><italic>b</italic></sub>) component (wavelength range: 0.3 &#x003BC;m to 3 &#x003BC;m) for the Baltoro Glacier, compared to field measurements of short-wave irradiance at the same location on July 24th 2004. The discrepancy in magnitude is mainly caused by the diffuse-skylight irradiance and the adjacent-terrain irradiance.</p></caption>
<graphic xlink:href="feart-09-652279-g0001.tif"/>
</fig>
<p>For many valley glaciers, the topographic shading on the direct solar radiation can be significant due to complex topography (Olson and Rupper, <xref ref-type="bibr" rid="B77">2019</xref>). The <italic>S</italic> parameter in Equation (1) accounts for the meso-scale relief structure of the topography that governs the presence of cast shadows on the landscape. Cast shadows on the landscape are governed by the size of the solar disk and the Earth-Sun distance. More specifically, the solar angular width (&#x003B1;<sub><italic>s</italic></sub>(<italic>t</italic>),[degrees]) in relation to topographic relief, will govern the overall length of a cast shadow in the direction of &#x003D5;<sub><italic>s</italic></sub>. This angle and other geometry parameters can be used to compute the planimetric length of the umbra and penumbra subregions that collectively represent the total cast-shadow region. The umbra is the region over the landscape where <italic>E</italic><sub><italic>b</italic></sub> is totally obstructed by the topography. Conversely, the penumbra region receives a fraction of <italic>E</italic><sub><italic>b</italic></sub> because a portion of the solar disk is obstructed by the topography. In the penumbra region, the fraction of <italic>E</italic><sub><italic>b</italic></sub> is governed by the pixel location. Consequently, a location in the umbra region would exhibit a shadow coefficient of 0.0, while a location in the penumbra region would exhibit a coefficient value ranging from 0.0 &#x0003C; <italic>S</italic> &#x0003C; 1.0. Therefore, cast shadow locations can also include irradiance from <italic>E</italic><sub><italic>b</italic></sub>. Furthermore, because of the eccentricity of the Earth&#x00027;s orbit around the sun, the length of the cast shadow region varies non-linearly with time. Given the extreme relief in many mountain environments, cast shadows can significantly modulate direct irradiance over the course of a day, let alone over the course of an ablation season. The importance of cast shadows should also not be underestimated, as the <italic>S</italic> parameter strongly regulates the direct irradiance during the early morning and late afternoon, as governed by basin relief and solar geometry. Large areas of a glacier&#x00027;s surface can be in cast shadow (umbra) for hours during the day, thereby blocking the direct irradiance and decreasing the surface irradiance and ablation potential. It is more difficult to evaluate the significance of the penumbra on ablation, other than to say that the length of the penumbra in the direction of the solar azimuth angle causes a direct irradiance gradient that influences ablation depending upon its length. Numerical modeling and sensitivity analysis are required to have a better understanding of the entire <italic>S</italic> parameter. Consequently, the <italic>S</italic> parameter can have a significant to insignificant influence on glacier ablation depending upon environmental conditions and the total length of the umbra and penumbra and time that cast shadows are found over the glacier surface. Its modulating effects on ablation will increase during the ablation season as solar zenith angles increase.</p>
<p>Simulating the magnitude and distribution of <italic>E</italic><sub><italic>b</italic></sub> on a glacier surface is relatively straight-forward, given that topographic evolution resulting from ablation dynamics is accounted for with respect to local topographic properties and relief. We demonstrate relatively accurate simulated estimates of diurnal direct irradiance (1 day; July 24, 2004) using a SRTM 30 m DEM, and compared our simulation to field-based surface irradiance measurements collected on the same day by Mihalcea et al. (<xref ref-type="bibr" rid="B69">2008</xref>) for the same location on the Baltoro Glacier (<xref ref-type="fig" rid="F1">Figure 1</xref>). We accounted for orbital parameters based upon Berger (<xref ref-type="bibr" rid="B14">1978</xref>) and used the apparent solar zenith angle accounting for parallax and atmospheric refraction, assuming the temperature and pressure profiles of the atmosphere. In this simulation, however, we did not utilize ray-tracing to account for cast shadows. The simulation results are reasonable, the modeled <italic>E</italic><sub><italic>b</italic></sub> is less than measured <italic>E</italic> values in the field, as we did not account for <italic>E</italic><sub><italic>d</italic></sub> and <italic>E</italic><sub><italic>t</italic></sub> that is represented in the <italic>E</italic> values. Perhaps it is more important to notice the discrepancy in magnitude in the early morning and late afternoon, caused by cast shadows (i.e., shape at the base of the distribution is different). Rapidly changing magnitudes at these times clearly reveals the influence of cast shadows in <italic>E</italic>, vs. simulated <italic>E</italic><sub><italic>b</italic></sub> values. We would expect these differences to become more pronounced during the later stages of the ablation season as the solar zenith angles increases.</p>
<p>Challenges in estimating this dominant component more accurately involve better characterization of atmospheric conditions, use of higher-resolution DEMs for establishing initial conditions, and the development and evaluation of different parameterization schemes to account for those controlling parameters that govern the energy-balance at the glacier surface.</p>
</sec>
<sec>
<title>2.1.2. Short-Wave Diffuse-Skylight Irradiance</title>
<p>Atmospheric scattering will produce a hemispherical source of irradiance that contains isotropic and anisotropic components (Perez et al., <xref ref-type="bibr" rid="B79">1986</xref>; Proy et al., <xref ref-type="bibr" rid="B82">1989</xref>). The anisotropic nature of <italic>E</italic><sub><italic>d</italic></sub> has been the largest source of error for estimating this component, as it is governed by circumsolar brightening due to forward scattering by aerosols, and horizontal brightening due to Rayleigh scattering (Perez et al., <xref ref-type="bibr" rid="B79">1986</xref>). A Two-stream approximation model assumes that there is a Rayleigh-scattering component (<italic>E</italic><sub><italic>r</italic></sub>), an aerosol-scattering component (<italic>E</italic><sub><italic>a</italic></sub>), and a secondary ground/sky-backscattering component (<italic>E</italic><sub><italic>g</italic></sub>) caused by multiple interactions between the ground surface and the atmosphere (Proy et al., <xref ref-type="bibr" rid="B82">1989</xref>; Gueymard, <xref ref-type="bibr" rid="B47">1995</xref>; Zhang et al., <xref ref-type="bibr" rid="B115">2015</xref>). One can account for single-scattering or a multiple-scattering Rayleigh atmosphere. The single-scattering albedo is required to account for aerosol scattering which is a function of both wavelength and humidity. Aerosol type and models such as rural, urban, maritime or tropospheric can also be accounted for Gueymard (<xref ref-type="bibr" rid="B47">1995</xref>). The ground backscatter component can be modeled by accounting for the zonal ground reflectance from the overall ground and sky reflectance (Gueymard, <xref ref-type="bibr" rid="B47">1995</xref>). Numerous parameterization schemes are available that enable a reasonable characterization of this irradiance component for a horizontal surface (e.g., Bird and Riordan, <xref ref-type="bibr" rid="B15">1986</xref>; Gueymard, <xref ref-type="bibr" rid="B47">1995</xref>). Multi-scale topographic effects, however, need to be considered at each location on the landscape (Proy et al., <xref ref-type="bibr" rid="B82">1989</xref>; Bishop et al., <xref ref-type="bibr" rid="B21">2019</xref>).</p>
<p>Both local and meso-scale topographic properties govern <italic>E</italic><sub><italic>d</italic></sub> and Proy et al. (<xref ref-type="bibr" rid="B82">1989</xref>) provides a computation solution such that:</p>
<disp-formula id="E3"><label>(3)</label><mml:math id="M3"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>E</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x003C0;</mml:mi><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mstyle><mml:msup><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02193;</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003B1;</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo class="qopname">cos</mml:mo><mml:mi>I</mml:mi><mml:mo class="qopname">sin</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>L</italic><sup>&#x02193;</sup> is the downward radiance from sky-hemispherical incident directions &#x003B8;<sub><italic>i</italic></sub> and &#x003D5;<sub><italic>i</italic></sub>, and &#x003B8;<sub><italic>i</italic></sub> and &#x003D5;<sub><italic>i</italic></sub> are the zenith and azimuth angles from the hemisphere, &#x003B1;<sub><italic>s</italic></sub> is the solar-elevation angle (&#x003B1;<sub><italic>s</italic></sub> &#x0003D; &#x003C0;/2 &#x02212; &#x003B8;<sub><italic>s</italic></sub>), and <italic>I</italic> is the incidence angle of the direction defined by sky-hemisphere and terrain geometry, similar to Equation (2) using incident geometry. Such a computation is computationally expensive given that <italic>E</italic><sub><italic>d</italic></sub> is wavelength dependent, and hemispherical topographic shielding must be accounted for by limiting the integration based upon the maximum relief angle for all hemispherical azimuth directions.</p>
<p>Basin relief production is primarily governed by paleo-glacier erosion and uplift history. There is a time disconnect such that the majority of the relief production occurred in the past, however the resulting relief (i.e., topographic shielding) partly influences the modern-day diffuse-skylight irradiance. Its altitudinal spatial variability can therefore govern the spatial variability in ablation. We demonstrate this topographic control using radiative and ablation simulations over the Baltoro Glacier in the central Karakoram (<xref ref-type="fig" rid="F2">Figure 2</xref>). For this simulation, we did not perform the full integration of <italic>E</italic><sub><italic>d</italic></sub> because of the computational intensity of doing so over all wavelengths. Rather we used the algorithms from Bird and Riordan (<xref ref-type="bibr" rid="B15">1986</xref>) to estimate <italic>E</italic><sub><italic>d</italic></sub> on a horizontal surface, and then multiplied the results by the so-called skyview-factor coefficient that represents the relative magnitude of topographic shielding. The computation of this topographic parameter is highly scale dependent, and Bishop and Dobreva (<xref ref-type="bibr" rid="B18">2017</xref>) demonstrated that it is highly variable across mountain environments given the polygenetic nature of erosion dynamics and topographic evolution (for more parameter details see Dozier et al., <xref ref-type="bibr" rid="B36">1981</xref>; Wilson and Bishop, <xref ref-type="bibr" rid="B110">2013</xref>). Regarding the ablation simulations, we account for debris thickness distribution over the Baltoro Glacier from Mihalcea et al. (<xref ref-type="bibr" rid="B69">2008</xref>) and meltwater ponds to generate a temporally-averaged ablation rate over the glacier. We then compared the skyview-factor coefficient to simulated ablation over the glacier (<xref ref-type="fig" rid="F2">Figure 2</xref>). Clearly there is a decrease in the ablation rate with altitude that accounts for the collective influence of cooler atmospheric temperatures and an increase in topographic shielding (lower skyview coefficient). At about 4, 400 &#x02212; 4, 500 m (below Concordia; Concordia is at about 4,691 m) topographic shielding begins to decrease and is associated with the confluence of glacier tributaries. It then abruptly increases with maximum topographic shielding being spatially coincident with the minimum ablation rate for the profile. These simulation results suggest that meso-scale topographic effects also play an important role governing the magnitude of <italic>E</italic><sub><italic>d</italic></sub> and ablation depending upon paleo-glacier erosion dynamics and basin relief production (i.e., erosion-uplift dynamics).</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Topographic shading effect on the simulated surface ablation rates for the Baltoro Glacier in 2004. The ablation rates are averaged over 30 m altitude bins. The skyview factor represents the topographic shielding that governs the diffuse-skylight irradiance. The abrupt drop in the skyview is associated with the minimum ablation along the profile, which highlights this meso-scale topographic effect on glacier surface ablation.</p></caption>
<graphic xlink:href="feart-09-652279-g0002.tif"/>
</fig>
<p>The significance of this secondary surface-irradiance parameter on ablation is difficult to ascertain, as numerical simulation sensitivity analysis focused on irradiance partitioning is sorely needed. We have demonstrated, however, that mesoscale topographic effects on the diffuse irradiance potentially regulate ablation rates based upon glacier erosion histories (i.e., relief production). Simulation results indicate that the magnitude of this parameter can range from relatively low or moderate energy levels in the visible portion of the spectrum (generally 5&#x02013;200 <italic>Wm</italic><sup>&#x02212;2</sup>) depending upon topographic conditions. Spatial distribution patterns of higher diffuse energy are most likely found at higher altitudes where basin relief significantly decreases. Furthermore, the influence of landcover conditions can also increase the magnitude of this component given high albedo surfaces (i.e., snow) and secondary scattering back to the surface.</p>
</sec>
<sec>
<title>2.1.3. Short-Wave Adjacent-Terrain Irradiance</title>
<p>The aforementioned irradiance components and surface Bi-Directional Reflectance Distribution Function (BRDF) interact with the surrounding terrain geometry to produce the adjacent-terrain irradiance (<italic>E</italic><sub><italic>t</italic></sub>). This irradiance component is extremely complicated but important to characterize, because variations in land cover and topographic complexity strongly regulate its magnitude. Highly reflective materials such as felsic minerals, vegetation, and ice and snow coupled with steep slopes can cause significant irradiance that impacts ablation dynamics, and especially the retreat of ice-cliffs (Sakai et al., <xref ref-type="bibr" rid="B95">2002</xref>; Buri et al., <xref ref-type="bibr" rid="B26">2016</xref>).</p>
<p>The characteristics of surface reflectance properties can range between diffuse and specular (Iqbal, <xref ref-type="bibr" rid="B59">1983</xref>; Zhang et al., <xref ref-type="bibr" rid="B115">2015</xref>). We know, however, that the isotropic assumption of reflectance is not valid in mountains, and the surface BRDF (&#x003C1;<sub><italic>brdf</italic></sub>) must be accounted for at each location on the landscape, as reflectance from the surrounding terrain and other glacier surface locations can contribute to <italic>E</italic><sub><italic>t</italic></sub>. The BRDF is used to describe the anisotropic nature of surface reflectance characteristics. It is a scattering function that describes anisotropic reflectance given all input-output angles. It is characterized as:</p>
<disp-formula id="E4"><label>(4)</label><mml:math id="M4"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003C1;</mml:mi></mml:mrow><mml:mrow><mml:mi>b</mml:mi><mml:mi>r</mml:mi><mml:mi>d</mml:mi><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msubsup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:msubsup><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:msubsup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:msubsup><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>L</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>E</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where, <italic>L</italic> is the surface radiance and <italic>E</italic> is the surface irradiance (<italic>E</italic> &#x0003D; <italic>E</italic><sub><italic>b</italic></sub>&#x0002B;<italic>E</italic><sub><italic>d</italic></sub>&#x0002B;<italic>E</italic><sub><italic>t</italic></sub>). The effective illumination (<inline-formula><mml:math id="M5"><mml:msubsup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>) and viewing directions (<inline-formula><mml:math id="M6"><mml:msubsup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>) are a function of solar geometry and terrain geometry. Consequently, the topography has a significant affect on the BRDF (Zhang et al., <xref ref-type="bibr" rid="B115">2015</xref>), and the intimate mixture of materials and structure of the surface also influence scattering direction. The BRDF, however, is difficult to measure accurately, and is also the basis for the computation of surface albedo, another important RT parameter that is difficult to estimate.</p>
<p>The irradiance component <italic>E</italic><sub><italic>t</italic></sub> can be characterized in the following fashion:</p>
<disp-formula id="E5"><label>(5)</label><mml:math id="M7"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>E</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msubsup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msubsup><mml:mo>,</mml:mo><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mtext>&#x000A0;</mml:mtext><mml:msubsup><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02193;</mml:mi><mml:mi>&#x02191;</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo class="qopname">cos</mml:mo><mml:msub><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>L</italic> is the adjacent surface reflected radiance coming from the effective incident direction, &#x003D5;<sub><italic>i</italic></sub> is the hemispherical incident azimuth angle, &#x003B8;<sub><italic>i</italic></sub> is the incident vertical hemispherical zenith angle to account for terrain radiance above and below a pixel location, <italic>T</italic><sup>&#x02193;&#x02191;</sup> is the atmospheric transmittance given the optical depth of the atmosphere due to relief and propagation zenith angle through the atmosphere (&#x003B8;<sub><italic>v</italic></sub>) between two locations, cos<italic>I</italic><sub><italic>t</italic></sub> represents the terrain incidence angle given the influence of the local terrain geometry in relation to the incident directional geometry, and <italic>S</italic><sub><italic>t</italic></sub> representing terrain blocking of the adjacent surface radiance between any two points of the surrounding terrain. As a general rule-of-thumb, the computation of <italic>E</italic><sub><italic>t</italic></sub> should account for the terrain conditions extending out to some distance (e.g., 5 km) from each pixel location (Proy et al., <xref ref-type="bibr" rid="B82">1989</xref>). However, the appropriate distance in any direction can vary as a function of the anisotropic nature of the topography, given the orientation of the structural fabric of the topography related to differential erosion and tectonics. As such, the adjacent hemispherical distribution of locations that can contribute to <italic>E</italic><sub><italic>t</italic></sub> for any particular location can be extremely scale dependent, given that it is also a function of altitude, such that the contributing zones of adjacent-terrain irradiance may be spatially contiguous or spatially fragmented. To our knowledge, this level of complexity has not been evaluated with respect to determining its impact on ablation dynamics, although it most likely plays a role in ice-cliff retreat when considering the coupled influence of the long-wave adjacent-terrain irradiance term.</p>
<p>The short-wave adjacent-terrain irradiance is thought to be a tertiary surface irradiance component that does not significantly affect ablation. However, given our parameterization scheme and the nature of highly directional and diffuse reflectance components associated with the BRDF, coupled with changing glacier surface topography, solar zenith and solar azimuth angles, the <italic>E</italic><sub><italic>t</italic></sub> component is likely to be highly variable in space and time, and facilitate ablation in general. The significance of this component could be greater than the diffuse irradiance given specific topographic and adjacent land cover conditions, such as at depression areas that are surrounded by highly reflective felsic minerals and snow.</p>
</sec>
</sec>
<sec>
<title>2.2. Precipitation</title>
<p>Precipitation in mountainous regions is partly modified by avalanche redistribution of snow, and is determined by uplift in traveling frontal zones or by orographic uplift. The former will move spatially with the frontal system while the latter will remain stationary given the topography and the wind direction. Determining where near-constant orographically-forced precipitation occurs may lead to an explanation of why a particular glacier exists at one location but not in another nearby location.</p>
<p>Separating orographic precipitation from frontal precipitation is difficult from observations (in which all precipitation is combined). For orographic precipitation one needs to know the wind direction, the atmospheric moisture content, and the topography itself but one also has to track the properties of the air parcels along their trajectories to keep track of the amounts of condensation, evaporation, and precipitation as the air travels up and perhaps over multiple peaks, potentially losing all of its water content along the way. Some attempts at quantifying this process have been made using a combination of a high resolution DEM and downscaled reanalysis climate data (e.g., Dobreva et al., <xref ref-type="bibr" rid="B35">2017</xref>).</p>
<p>Regardless of the driving mechanism behind mountain precipitation, the precipitate-formation process is mostly a cold-based one given typical surface temperatures at high elevations and the adiabatic lapse rate, i.e., the precipitate-formation process occurs at temperatures below freezing and so is governed by the Bergeron-Findeisen process. According to this process, in a cloud that has a mix of ice crystals and supercooled water droplets, water vapor in the cloud preferentially deposits on the ice crystals because of the lower saturation vapor pressure over ice than over water. As the water vapor decreases through deposition on the ice, the supercooled droplets evaporate and lose their water to the growing ice crystals, which may then grow large enough to overcome any updrafts and fall to the surface. This process dominates mid-latitude precipitation (even during the summer) and at high altitudes.</p>
<p>Most precipitation in mountainous regions therefore begins frozen; whether it remains frozen on its way down to the surface is dependent on the atmospheric temperatures through which it falls. At high altitudes, it is more likely to remain frozen than, for example, at sea level in mid-latitudes, but if conditions are warm enough then it will melt completely and fall as rain (or, if the precipitate does not have enough time to melt completely it will fall as graupel, sleet, or freezing rain).</p>
<p>Snowfall above the equilibrium line altitude (ELA) is paramount for a glacier&#x00027;s survival. Climate change can alter the ELA up or down, depending on local conditions, but with a general warming trend ELAs on average are rising. If the ELA rises above the peak of the mountain then no glacier can exist. If atmospheric conditions above the ELA become conducive to rain then glacier ablation will be accelerated.</p>
<p>When debris cover is present additional factors must be taken into account. If the debris temperature is below freezing, snowfall would accumulate as on a debris-free glacier. If the debris cover is above freezing, however, then it will melt the snow leading to percolation of liquid water onto the glacier surface that can enhance ablation. Refreezing of the percolated water at the glacier surface is possible, with the release of latent heat into the interstitial spaces of the debris cover (Collier et al., <xref ref-type="bibr" rid="B32">2014</xref>, <xref ref-type="bibr" rid="B30">2015</xref>). The debris cover fundamentally alters the fluxes of sensible and latent heat to the overlying atmosphere (Collier et al., <xref ref-type="bibr" rid="B32">2014</xref>). In addition, a relatively warm debris cover (through daytime solar heating) can trigger convection in the atmosphere above and produce precipitation that melts on contact and increases surface ablation (Collier et al., <xref ref-type="bibr" rid="B30">2015</xref>). So DCGs have the potential to increase their own ablation rates through dynamic interactions with the overlying atmosphere as well as the underlying glacier.</p>
</sec>
</sec>
<sec id="s3">
<title>3. Supraglacial Debris Load</title>
<sec>
<title>3.1. Composition and Properties</title>
<p>Investigators are increasingly recognizing the significance of the supraglacial debris load as a major controlling factor in regulating ablation and the morphological evolution of DCGs (Adhikary et al., <xref ref-type="bibr" rid="B1">2000</xref>; Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>; Reid et al., <xref ref-type="bibr" rid="B87">2012</xref>; Rowan et al., <xref ref-type="bibr" rid="B94">2015</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Zhang et al., <xref ref-type="bibr" rid="B114">2016</xref>; Gibson et al., <xref ref-type="bibr" rid="B46">2017</xref>). Field data and numerical-modeling efforts (given various assumptions) demonstrate the general relationship between debris thickness and ablation rate, as thicker debris loads decrease the thermal energy available for ablation at the debris/ice interface, while thin and moisture laden debris can significantly accelerate ablation beyond that of debris-free conditions (Mattson, <xref ref-type="bibr" rid="B66">1993</xref>; Kayastha et al., <xref ref-type="bibr" rid="B62">2000</xref>; Reznichenko et al., <xref ref-type="bibr" rid="B90">2010</xref>). This general relationship, however, does not account for a variety of debris-load properties, processes, and feedbacks regarding ablation dynamics, as supported by findings related to the &#x0201C;debris-covered glacier anomaly&#x0201D; (Salerno et al., <xref ref-type="bibr" rid="B97">2017</xref>) and numerical modeling efforts (Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Zhang et al., <xref ref-type="bibr" rid="B114">2016</xref>; Gibson et al., <xref ref-type="bibr" rid="B46">2017</xref>; Huo et al., <xref ref-type="bibr" rid="B54">2020</xref>). Consequently, the multi-faceted nature of this controlling factor is not adequately understood, as there are many unknowns with regard to properties and processes that govern ablation dynamics, and research clearly reveals the relatively high spatio-temporal variability in debris properties and processes that should be accounted for Scherler et al. (<xref ref-type="bibr" rid="B99">2011b</xref>); Reid et al. (<xref ref-type="bibr" rid="B87">2012</xref>); Rowan et al. (<xref ref-type="bibr" rid="B94">2015</xref>); Yue et al. (<xref ref-type="bibr" rid="B112">2017</xref>). The complexity of the problem is highlighted by the relatively large number of debris, near-surface and scale-dependent topographic parameters that partially govern ablation. <xref ref-type="table" rid="T1">Table 1</xref> lists various physical and radiative-transfer properties that all contribute toward regulating the magnitude of ablation in space and time. It is also important to realize that location/scale specific parameters related to the topography and lithology beneath and surrounding a glacier also contribute to spatial and property variations in the debris load, as debris sources and mineralogical composition of debris and rocks is governed by the nature of the geological setting, and the geomorphological system within the glacial basin.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Debris-cover physical and radiative-transfer properties that regulate ablation in space and time.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>Physical properties</bold></th>
<th valign="top" align="left"><bold>Symbol [SI units]</bold></th>
<th valign="top" align="left"><bold>Relevance</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Debris thickness</td>
<td valign="top" align="left"><italic>h</italic><sub><italic>d</italic></sub> [m]</td>
<td valign="top" align="left">Governs distribution and magnitude of ablation and thermal gradient</td>
</tr>
<tr>
<td valign="top" align="left">Mineralogical/rock composition</td>
<td/>
<td valign="top" align="left">Governs radiation transfer properties and thermal gradient</td>
</tr>
<tr>
<td valign="top" align="left">Porosity/particle size</td>
<td valign="top" align="left">&#x003C6; [dimensionless]</td>
<td valign="top" align="left">Governs thermal gradient and bulk thermal conductivity</td>
</tr>
<tr>
<td valign="top" align="left">Moisture content</td>
<td valign="top" align="left"><italic>w</italic> [dimensionless]</td>
<td valign="top" align="left">Governs thermal properties and temperature profiles</td>
</tr>
<tr>
<td valign="top" align="left">Slope</td>
<td valign="top" align="left">&#x003B8;<sub><italic>t</italic></sub> [rad]</td>
<td valign="top" align="left">Governs surface irradiance and sediment transport</td>
</tr>
<tr>
<td valign="top" align="left">Slope Azimuth</td>
<td valign="top" align="left">&#x003D5;<sub><italic>t</italic></sub> [rad]</td>
<td valign="top" align="left">Governs surface irradiance and sediment transport</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left">Temperature</td>
<td valign="top" align="left"><italic>T</italic><sub><italic>d</italic></sub> [<italic>K</italic>]</td>
<td valign="top" align="left">Governs long-wave irradiant flux and convective atmospheric motion</td>
</tr>
<tr style="border-bottom: thin solid #000000;">
<td valign="top" align="left" colspan="3"><bold>Radiative-Transfer Properties</bold></td>
</tr> <tr>
<td valign="top" align="left">Albedo</td>
<td valign="top" align="left">&#x003B1; [dimensionless]</td>
<td valign="top" align="left">Governs irradiance</td>
</tr>
<tr>
<td valign="top" align="left">Emissivity</td>
<td valign="top" align="left">&#x003B5; [dimensionless]</td>
<td valign="top" align="left">Governs irradiance</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>An important first-order property that requires consideration is the particle size/distribution of the debris load. This property can be highly spatio-temporally variable over a glacier, and is regulated by source and sediment-transport dynamics. It can have a significant impact on ablation dynamics, as &#x0201C;particles&#x0201D; can include dust, fines, gravel, cobbles, rocks, isolated boulders, and boulder fields (<xref ref-type="fig" rid="F3">Figure 3</xref>; the full range of debris particle sizes). Such heterogeneous conditions are common on many glaciers around the world, and &#x0201C;particle&#x0201D; distributions can be significantly different on glaciers within a region, depending upon paleo-climate-glacier dynamics (Bush et al., <xref ref-type="bibr" rid="B27">2020</xref>). Such variations are the result of polygenetic sediment-transport dynamics involving mass movements, aeolean, suprafluvial, englacial, gravitational, and ice-flow transport mechanisms. Spatio-temporal variations of this aspect of the debris load are usually not considered, and homogeneous unconsolidated debris conditions are assumed, which causes uncertainties in ablation estimates.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Supraglacial debris cover conditions on the Baltoro Glacier during the summer of 2005. The nature of the supraglacial debris load can vary significantly with respect to particle size, composition, and debris thicknesses. <bold>(A)</bold> Cryoconite holes at Concordia at 4,500 m. These features form as wind blown dust, soot and particles that absorb radiation generate holes that contain water that contributes further to increase ablation and the depth and width of the hole. <bold>(B)</bold> Wind blown silt and sand deposits at the terminus of the glacier. Wind turbulence caused by up valley winds and altitude variations between the pro-glacial valley floor and the glacier surface results in wind-sorted supraglacial debris. <bold>(C)</bold> Supraglacial debris load consisting of gravel, pebbles, and cobble-sized particles. Notice the felsic and mafic compositional variations of particles. <bold>(D)</bold> Supraglacial debris dominated by larger rocks over the glacier surface. Rock spatial density greatly governs the spatial variation in ablation dues to shadowing and rock mass influence on the vertical debris load thermal gradient. <bold>(E)</bold> Debris load consisting of sediment, rocks, and boulders. Debris thickness is about 5 m. Debris thickness can be highly variable over a glacier surface depending upon a multitude of controlling factors. <bold>(F)</bold> Debris loads in the Himalaya frequently consist of boulder fields, that make it difficult to accurately assess or model thermal gradients, given shadowing and alternation of the wind field. Photo credit: Michael P. Bishop.</p></caption>
<graphic xlink:href="feart-09-652279-g0003.tif"/>
</fig>
<p>The composition of the debris load is another important property. It is rarely accounted for, although it can be taken into consideration by using estimates of the debris bulk conductivity. It not only governs the bulk thermal properties of the debris, but also the surface reflectance properties and the surface energy balance of the short-wave irradiance (i.e., albedo). Consequently, remote sensing analysis and absorption-feature modeling can be used for mineral detection and determining the composition of the debris using hyperspectral data. This approach can provide more explicit estimates of the surface composition and distribution of moisture-laden debris. Such data regarding the debris load are too logistically difficult to acquire in the field, but high-resolution hyperspectral remote-sensing capabilities are now available using uncrewed aerial systems (UAS) to better assess &#x0201C;particle&#x0201D; distributions, surface mineralogy and other parameters that are not usually accounted for (e.g., Fyffe et al., <xref ref-type="bibr" rid="B44">2020</xref>).</p>
<p>Debris load thickness has been shown to be highly variable (Mihalcea et al., <xref ref-type="bibr" rid="B69">2008</xref>; Fyffe et al., <xref ref-type="bibr" rid="B44">2020</xref>). Many glaciers exhibit low-frequency debris thickness variation (the generally thicker debris toward the terminus) controlled by surface ice-flow transport dynamics and englacial transport of sediment and rocks into the terminus region (Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Wirbel et al., <xref ref-type="bibr" rid="B111">2018</xref>). Higher-frequency debris thickness variations must also be considered, as the debris load is modulated by supra-fluvial transport, topography-driven gravitational transport, and other ablation dynamics including supraglacial lake development and evolution. Furthermore, a high percentage of debris cover lowers the glacier surface albedo that partially counters the debris insulation effect (Fyffe et al., <xref ref-type="bibr" rid="B44">2020</xref>). Such complexities have not been adequately investigated, and it is essential to note that complex feedbacks exist between sediment fluxes, ablation, topography and the surface-energy balance (Huo et al., <xref ref-type="bibr" rid="B54">2020</xref>). Although several approaches have been developed to estimate the distribution of debris thickness (e.g., Mihalcea et al., <xref ref-type="bibr" rid="B69">2008</xref>; Zhang et al., <xref ref-type="bibr" rid="B113">2011</xref>; Reid et al., <xref ref-type="bibr" rid="B87">2012</xref>; Rounce et al., <xref ref-type="bibr" rid="B92">2018</xref>), higher-frequency variations in debris thickness and the spatial variability of ablation is not known with any degree of certainty on most DCGs.</p>
<p>Field measurements suggest that debris temperature gradients are approximately stable and linear with depth during the ablation season (Rowan et al., <xref ref-type="bibr" rid="B93">2021</xref>), although the linearity and stability depend on the time scale, and linearity is also dependent upon the debris thickness. Assuming linearity, debris surface temperature can be used to estimate sub-debris ablation, which is also controlled by debris thickness, and thermal conductivity (Nakawo and Young, <xref ref-type="bibr" rid="B74">1981</xref>; Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>). Given shallow debris thicknesses, relative homogeneity in particle distribution, and lack of moisture-laden sediment, this first-order approximation assumes isotropic conditions with respect to various properties. With greater debris load thicknesses, however, porosity variation has the potential to alter the thermal gradient given the translocation of the fine size-fraction of particles and presence of moisture-laden sediment. Therefore, the thermal conductivity of debris can be difficult to determine given the complex composition and pore-space conditions. Nicholson and Benn (<xref ref-type="bibr" rid="B75">2006</xref>) presented an approach to address this issue using measured temperature profiles. Based on this model, the bulk thermal conductivity of the debris layer (<italic>k</italic><sub><italic>d</italic></sub>) can be written as:</p>
<disp-formula id="E6"><label>(6)</label><mml:math id="M8"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x003C1;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x003C1;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow></mml:msub><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mtext>&#x000A0;</mml:mtext><mml:mi>&#x003BA;</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>c</italic><sub><italic>d</italic></sub>, &#x003C1;<sub><italic>d</italic></sub>, <italic>c</italic><sub><italic>v</italic></sub>, &#x003C1;<sub><italic>v</italic></sub> are the specific heat capacity and density of debris, and specific heat capacity and density of the pore-space filler, respectively. &#x003C6; is the porosity, and &#x003BA; is the apparent thermal diffusivity of the debris, which can be determined using the thermal diffusion equation:</p>
<disp-formula id="E7"><label>(7)</label><mml:math id="M9"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mi>&#x003BA;</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x02202;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:msup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>T</italic><sub><italic>d</italic></sub> is the debris temperature, which is a function of depth (<italic>z</italic>) and time (<italic>t</italic>). &#x003BA; can then by determined from the measured temperatures-time and temperatures-depth variations. Another method was presented by Reid and Brock (<xref ref-type="bibr" rid="B88">2010</xref>), based on which, the debris was broken down to multiple thin layers to estimate the thermal conductivity and internal temperature profiles:</p>
<disp-formula id="E8"><label>(8)</label><mml:math id="M10"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003C1;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>&#x02202;</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:mi>z</mml:mi></mml:mrow></mml:mfrac><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:mi>z</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>Debris moisture conditions, mineral mixing, debris surface structure, and surface topography can also significantly influence the surface albedo that regulates debris temperature and ablation (Pelto, <xref ref-type="bibr" rid="B78">2000</xref>; Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>). Accurate estimation of the albedo based upon the surface BRDF is difficult to compute, and rapidly changing glacier-surface dynamics should cause the surface albedo to exhibit relatively high spatio-temporal variability (Takeuchi and Li, <xref ref-type="bibr" rid="B105">2008</xref>; Yue et al., <xref ref-type="bibr" rid="B112">2017</xref>). Therefore, better estimates of debris surface reflective parameters are required using more reliable parameterization schemes.</p>
<p>We demonstrate this using numerical simulations and accounting for the spectral reflectance of minerals that make of the debris load of the Baltoro Glacier. Specifically, we utilize a radiation-transfer approach (Coakley, <xref ref-type="bibr" rid="B29">2003</xref>) to compute the broadband surface albedo (&#x003B1;) based on the mixing of sediments and water at each grid cell, such that:</p>
<disp-formula id="E9"><label>(9)</label><mml:math id="M11"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mi>&#x003B1;</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:msubsup></mml:mstyle><mml:mi>L</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:msubsup></mml:mstyle><mml:mi>E</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>E</italic> is the total short-wave irradiance and <italic>L</italic> is the reflected surface radiance, which, based on an isotropic (Lambertian) assumption, can be written as (Bishop and Colby, <xref ref-type="bibr" rid="B17">2011</xref>):</p>
<disp-formula id="E10"><label>(10)</label><mml:math id="M12"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mi>L</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mfrac><mml:mrow><mml:mi>E</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>r</italic><sub><italic>mc</italic></sub> is the spectral reflectance of the mineralogical composite representing the debris, which is considered as a composite reflectance of different types of minerals. The debris lithological distribution over the Baltoro Glacier was based on Gibson et al. (<xref ref-type="bibr" rid="B46">2017</xref>), which indicates that the surface is mostly covered by a composite of gneiss (about 53%), granite (about 30%), and schist metasediment (about 17%). Simulated results in <xref ref-type="fig" rid="F4">Figure 4</xref> reveal significant spatial variations in the surface albedo over the lower ablation zone of the Baltoro Glacier caused by mineralogical and moisture variations. We performed the narrow-band to broad-band conversion of surface albedo using ASTER imagery based on the work of Liang (<xref ref-type="bibr" rid="B64">2001</xref>), and confirmed that our simulated values are very close to the remote-sensing approach. In the future, the use of semi-empirical BRDF models should be used with multi-angle reflectance imagery to generate higher quality albedo estimates for glacier surfaces.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Modeled surface albedo over the lower ablation zone of the Baltoro Glacier. The geomorphological conditions were estimated from an ASTER scene acquired on August 14th, 2004 (AST_08_00308142004054614).</p></caption>
<graphic xlink:href="feart-09-652279-g0004.tif"/>
</fig>
</sec>
<sec>
<title>3.2. Debris Transport</title>
<p>Debris loads are mobile due to ice-flow, gravitational movement, and fluvial processes (Anderson, <xref ref-type="bibr" rid="B6">2000</xref>; Benn and Evans, <xref ref-type="bibr" rid="B9">2014</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>, <xref ref-type="bibr" rid="B4">2018</xref>; Zhang et al., <xref ref-type="bibr" rid="B114">2016</xref>; van Woerkom et al., <xref ref-type="bibr" rid="B107">2019</xref>; Fyffe et al., <xref ref-type="bibr" rid="B44">2020</xref>). Although the exact glacial debris production mechanisms have not been fully understood, most investigators consider the main sources are from headwalls, sidewalls, basal erosion and moraines (Anderson, <xref ref-type="bibr" rid="B6">2000</xref>; Benn and Owen, <xref ref-type="bibr" rid="B12">2002</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; van Woerkom et al., <xref ref-type="bibr" rid="B107">2019</xref>). Studies have observed mass transport from landslides (Benn and Evans, <xref ref-type="bibr" rid="B9">2014</xref>), avalanches (Benn and Lehmkuhl, <xref ref-type="bibr" rid="B11">2000</xref>), and rockfalls from valley walls (Benn and Evans, <xref ref-type="bibr" rid="B9">2014</xref>). Englacial debris loads and rock encapsulated by basal erosion also contribute to the total debris mass balance, and several models have been developed to investigate the transport of englacial debris loads (Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Wirbel et al., <xref ref-type="bibr" rid="B111">2018</xref>).</p>
<p>One of the first models of supraglacial debris flux was developed by Anderson (<xref ref-type="bibr" rid="B6">2000</xref>), which describes the diffusive debris flux originated from medial moraines in valley glaciers. Anderson and Anderson (<xref ref-type="bibr" rid="B3">2016</xref>) and Anderson and Anderson (<xref ref-type="bibr" rid="B4">2018</xref>) developed a similar model that couples with ice dynamics to investigate the longer-term evolution of debris cover. Moore (<xref ref-type="bibr" rid="B73">2018</xref>) presented a comprehensive study on the gravitational debris transport across the ice surface, as the topography evolves due to melt. A theoretical framework was developed for assessing slope stability and gravitational mass transport accounting for meltwater balance on ablating ice. Modeled results provided insights into the geometry of stable slopes as a function of debris thickness and texture. Another recent study by Fyffe et al. (<xref ref-type="bibr" rid="B44">2020</xref>) provides valuable field observations on the rates and forms of debris transport. van Woerkom et al. (<xref ref-type="bibr" rid="B107">2019</xref>) studied the debris transport from lateral moraines using high resolution DEMs and found that the lateral moraines contribute to debris thickening along the margin of the glacier surface.</p>
<p>Here we discuss two dominant processes that govern debris movement on a glacier at different spatio-temporal scales: the gravitational process that operates at smaller spatio-temporal scales given complex glacier topography, and the advective-diffusive processes that control debris flux over larger spatio-temporal scales governed by the ice-flow.</p>
<sec>
<title>3.2.1. Gravitational Debris Fluxes</title>
<p>Gravitational debris movement (such as sliding and slumping) occurs on hillslopes and on glacier surfaces (Benn and Evans, <xref ref-type="bibr" rid="B9">2014</xref>; Fyffe et al., <xref ref-type="bibr" rid="B44">2020</xref>). Gravity-driven debris flux regulates the local thickness distribution of supraglacial debris during the ablation season when the surface topography is constantly changing under rapid melting, especially around melt hotspots such as supraglacial ponds and ice-cliffs. Field observations have identified sediment sliding or slumping off steepening ice-cliffs due to ice-cliff retreat and supraglacial lake expansion (Buri et al., <xref ref-type="bibr" rid="B26">2016</xref>; Miles et al., <xref ref-type="bibr" rid="B71">2017</xref>).</p>
<p>Gravity, internal friction and basal resistance force govern the gravitational movement of debris. Based on the debris flow model by Chen and Lee (<xref ref-type="bibr" rid="B28">2000</xref>), the unit net force acting on a debris column, <italic>F</italic>, can be written as:</p>
<disp-formula id="E12"><label>(11)</label><mml:math id="M14"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C1;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mi>g</mml:mi><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mfrac><mml:mo>-</mml:mo><mml:mi>k</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msqrt><mml:mrow><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msqrt></mml:mrow></mml:mfrac><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msqrt><mml:mrow><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msqrt></mml:mrow></mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E14"><label>(12)</label><mml:math id="M16"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>F</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C1;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mi>g</mml:mi><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mfrac><mml:mo>-</mml:mo><mml:mi>k</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>y</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msqrt><mml:mrow><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msqrt></mml:mrow></mml:mfrac><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msqrt><mml:mrow><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msqrt></mml:mrow></mml:mfrac><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>u</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>&#x003D5;</mml:mi></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003C1;<sub><italic>d</italic></sub> is debris density, <italic>g</italic> is the gravitational acceleration, <italic>z</italic><sub><italic>x</italic></sub> and <italic>z</italic><sub><italic>y</italic></sub> are the first derivatives of ice-surface elevation in the horizontal and vertical directions, respectively, <italic>k</italic> is the pressure ratio as defined by Chen and Lee (<xref ref-type="bibr" rid="B28">2000</xref>), <italic>u</italic><sub><italic>x</italic></sub> and <italic>u</italic><sub><italic>y</italic></sub> are debris velocity components, <italic>r</italic><sub><italic>u</italic></sub> is the constant pore-pressure ratio, and &#x003D5; is the dynamic internal friction angle of the debris. The first term on the right represents the gravitational force, the second term describes the internal friction of debris particles, and the third term represents the resistance force at the base of the debris column. Gravitational debris flux governs the local redistribution of debris thickness during the ablation season when the ice surfaces topography changes rapidly due to melt. The transport of a thin debris cover, however, is strongly controlled by more complex processes such as fluvial transport and rain wash (Fyffe et al., <xref ref-type="bibr" rid="B44">2020</xref>).</p>
</sec>
<sec>
<title>3.2.2. Ice-Flow and Debris Advection</title>
<p>Ice-flow describes the glacier motion that governs the advective debris transport from the production zone to the terminus (Benn et al., <xref ref-type="bibr" rid="B8">2012</xref>; Rowan et al., <xref ref-type="bibr" rid="B94">2015</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Wirbel et al., <xref ref-type="bibr" rid="B111">2018</xref>). Therefore, ice-flow dynamics play a fundamental role in debris transport as demonstrated by recent numerical simulations (e.g., Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Wirbel et al., <xref ref-type="bibr" rid="B111">2018</xref>).</p>
<p>Ice-flow velocity is governed by multiple factors such as the ice thickness and basal water pressure, and many DCGs exhibit high flow velocities even though they have stable termini (Copland et al., <xref ref-type="bibr" rid="B33">2009</xref>; Quincey et al., <xref ref-type="bibr" rid="B84">2009</xref>, <xref ref-type="bibr" rid="B83">2011</xref>). The surface velocity field can be determined from remote sensing analysis, estimating englacial ice-flow velocities can be challenging, although estimates are usually sufficient for computing particle trajectories if far enough away from the bed (about 80% of the ice depth). Therefore, most studies use ice-flow models (such as the shallow ice approximation) to solve the englacial velocity field (Herman and Braun, <xref ref-type="bibr" rid="B51">2008</xref>; Bueler and Brown, <xref ref-type="bibr" rid="B25">2009</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>), which usually requires ice thickness measurements from geophysical surveys using seismic or radio sounding (McNabb et al., <xref ref-type="bibr" rid="B67">2012</xref>). Modeled ice-flow velocities often suffer from high uncertainty (Farinotti et al., <xref ref-type="bibr" rid="B38">2017</xref>), however, more accurate models have been developed to account for basal melt, temperature-adjusted ice properties and rugged bed topography (e.g., Bueler and Brown, <xref ref-type="bibr" rid="B25">2009</xref>; Egholm et al., <xref ref-type="bibr" rid="B37">2011</xref>).</p>
<p>Debris advection is governed by the ice-flow velocity field (Benn and Evans, <xref ref-type="bibr" rid="B9">2014</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>; Wirbel et al., <xref ref-type="bibr" rid="B111">2018</xref>). A recent model was developed by Anderson and Anderson (<xref ref-type="bibr" rid="B3">2016</xref>), in which the englacial and supraglacial debris advection were modeled under a steady debris input to understand the mechanisms in the debris-glacier-climate system. Simulations indicated that debris has significant control on glacier length and gradients of ice discharge, ice thickness, and surface velocities. Their model demonstrated that high debris flux slows down the glacier and contributes to extending its length. Based on this, the rate of change of supraglacial debris thickness in the ablation zone can be written as (Anderson, <xref ref-type="bibr" rid="B6">2000</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>, <xref ref-type="bibr" rid="B4">2018</xref>):</p>
<disp-formula id="E15"><label>(13)</label><mml:math id="M17"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:mi>&#x003C6;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003C1;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>-</mml:mo><mml:mo>&#x02207;</mml:mo><mml:mo>&#x000B7;</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>u</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>s</mml:mtext></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>h</italic><sub><italic>d</italic></sub> is the surface debris thickness, <italic>t</italic> is time, <italic>c</italic><sub><italic>s</italic></sub> is the near-surface debris concentration, <italic>M</italic><sub><italic>s</italic></sub> is surface ablation rate, and <bold>u</bold><sub><bold>s</bold></sub> is surface ice velocity. The use of advection is valid for glacial debris transport because advection is defined as the transport of materials due to the bulk motion of a fluid, and glacier ice is a form of a viscoelastic fluid that is the basis for all modern ice-flow models. <xref ref-type="fig" rid="F5">Figure 5D</xref> shows the simulated supraglacial debris advection based on Equation (13) given the initial debris thickness (<xref ref-type="fig" rid="F5">Figure 5A</xref>) and surface velocity (<xref ref-type="fig" rid="F5">Figure 5C</xref>).</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Supraglacial debris conditions for the Baltoro Glacier. <bold>(A)</bold> The debris-thickness distribution over the Baltoro glacier in summer 2004 computed from ASTER surface kinetic temperature data following the approach by Mihalcea et al. (<xref ref-type="bibr" rid="B69">2008</xref>). <bold>(B)</bold> Debris rock type distribution over the Baltoro Glacier in 2004 based on the remote-sensing analysis by Gibson et al. (<xref ref-type="bibr" rid="B46">2017</xref>). <bold>(C)</bold> Surface ice-flow velocity of the Baltoro Glacier estimated from a Landsat-8 OLI panchromatic image pair acquired on September 15, 2004 and September 2, 2005. <bold>(D)</bold> Simulated debris transport due to ice-flow (advection process) over a 50-year period. Result represents the debris thickness distribution for year 50.</p></caption>
<graphic xlink:href="feart-09-652279-g0005.tif"/>
</fig>
<p>The transport of englacial debris and its contribution to the surface debris load is controlled by ice-flow dynamics, given that the vertical component of ice velocity usually points upward in the lower glacier where ice thickness increases toward the terminus. These processes have been described by Anderson (<xref ref-type="bibr" rid="B6">2000</xref>), and demonstrated in the simulations by Anderson and Anderson (<xref ref-type="bibr" rid="B3">2016</xref>) and Wirbel et al. (<xref ref-type="bibr" rid="B111">2018</xref>). The contribution of englacial debris load may explain the debris mass balance issue discussed by van Woerkom et al. (<xref ref-type="bibr" rid="B107">2019</xref>). The most recent englacial debris transport model was presented by Wirbel et al. (<xref ref-type="bibr" rid="B111">2018</xref>), in which both the advection and diffusion of englacial debris are accounted for:</p>
<disp-formula id="E16"><label>(14)</label><mml:math id="M18"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>&#x02202;</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mo>&#x02207;</mml:mo><mml:mo>&#x000B7;</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>D</mml:mi><mml:mo>&#x02207;</mml:mo><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:mo>&#x02207;</mml:mo><mml:mo>&#x000B7;</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>u</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>e</mml:mtext></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mi>r</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>c</italic><sub><italic>e</italic></sub> is the concentration of englacial debris, <italic>D</italic> is the diffusion coefficient, <bold>u</bold><sub><bold>e</bold></sub> is englacial ice velocity, and <italic>r</italic> represents englacial debris sources or sinks.</p>
<p>Remote-sensing approaches have been used to map debris thickness distributions (Mihalcea et al., <xref ref-type="bibr" rid="B69">2008</xref>; Juen et al., <xref ref-type="bibr" rid="B60">2014</xref>; Gibson et al., <xref ref-type="bibr" rid="B46">2017</xref>; Rounce et al., <xref ref-type="bibr" rid="B92">2018</xref>); the results, however, suffer from high uncertainties due to the instantaneous nature of satellite imagery and the complex properties of debris. Further investigations on debris production and transport modeling are required to better address the spatial patterns of surface melt and supraglacial morphological conditions.</p>
</sec>
</sec>
</sec>
<sec id="s4">
<title>4. Ablation Dynamics</title>
<p>Surface melt dominates the ablation of a DCG (Hambrey et al., <xref ref-type="bibr" rid="B49">2008</xref>; Reid and Brock, <xref ref-type="bibr" rid="B86">2014</xref>; Rowan et al., <xref ref-type="bibr" rid="B94">2015</xref>). The melt rate is strongly affected by the properties of debris loads and the presence of supraglacial water bodies due to their strong impact on the energy transfer dynamics (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>; Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>; Reid and Brock, <xref ref-type="bibr" rid="B88">2010</xref>; Fyffe et al., <xref ref-type="bibr" rid="B43">2014</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>).</p>
<p>Field measurement of ablation is usually conducted through monitoring ablation stakes at multiple locations on the glacier surface. Ablation rates at locations with different debris thickness are measured by several studies, and field data showed that debris thickness is a predominant factor that governs the sub-debris ablation rate (Mattson, <xref ref-type="bibr" rid="B66">1993</xref>; Kayastha et al., <xref ref-type="bibr" rid="B62">2000</xref>; Rounce et al., <xref ref-type="bibr" rid="B91">2015</xref>). Given the complex debris effects, empirical approaches such as the degree-day method for estimating ablation may be highly unreliable for DCGs (Braithwaite and Olesen, <xref ref-type="bibr" rid="B23">1990</xref>; Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>). Therefore, physics-based surface energy balance models are widely used for estimating sub-debris ablation (e.g., Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>; Reid and Brock, <xref ref-type="bibr" rid="B88">2010</xref>, <xref ref-type="bibr" rid="B86">2014</xref>; Fyffe et al., <xref ref-type="bibr" rid="B43">2014</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>), and several remote-sensing approaches have been developed to estimate the key parameters, such as the thickness and thermal properties of debris (e.g., Mihalcea et al., <xref ref-type="bibr" rid="B69">2008</xref>; Zhang et al., <xref ref-type="bibr" rid="B113">2011</xref>; Foster et al., <xref ref-type="bibr" rid="B40">2012</xref>; Juen et al., <xref ref-type="bibr" rid="B60">2014</xref>).</p>
<p>Several recent studies have successfully used the surface energy balance model to estimate ablation on DCGs. For example, Fyffe et al. (<xref ref-type="bibr" rid="B43">2014</xref>) used a distributed model to compute melt rate on the Miage glacier in Alps. Rounce et al. (<xref ref-type="bibr" rid="B91">2015</xref>) compared the influences of debris thermal conductivity, albedo, and surface roughness on the surface energy balance of the Imja-Lhotse Shar glacier in Nepal. These results highlighted that: (1) Unlike the thick debris layer, a thin debris is sensitive to air temperature variation and water content, therefore areas with wet thin debris must be identified. (2) Subdebris ablation is more sensitive to the changes in thermal conductivity than surface roughness and albedo. (3) It is important to account for the realistic topographic shielding and wind speed over the glacier surface in the energy balance model. (4) Many model parameters need to be calibrated using field measurements or remote sensing data as they can be highly site-specific. The energy balance model at the debris/air interface can be written as (Nakawo and Young, <xref ref-type="bibr" rid="B74">1981</xref>; Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>; Zhang et al., <xref ref-type="bibr" rid="B113">2011</xref>; Rounce et al., <xref ref-type="bibr" rid="B91">2015</xref>):</p>
<disp-formula id="E17"><label>(15)</label><mml:math id="M19"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>Q</italic><sub><italic>s</italic></sub> is the net shortwave radiation flux, <italic>Q</italic><sub><italic>l</italic></sub> is the net long-wave radiation flux, <italic>Q</italic><sub><italic>h</italic></sub> is the net sensible-heat flux, <italic>Q</italic><sub><italic>e</italic></sub> is the net latent-heat flux, and <italic>Q</italic><sub><italic>c</italic></sub> is the conductive heat flux into the debris which governs the ablation rate. Net radiation fluxes dominate glacier surface energy balance and significantly control surface ablation rates (Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>; Cuffey and Paterson, <xref ref-type="bibr" rid="B34">2010</xref>). The net shortwave radiation flux and the net long-wave radiation flux can be computed as:</p>
<disp-formula id="E18"><label>(16)</label><mml:math id="M20"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>E</mml:mi></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>E</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>E</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:mi>&#x003B1;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E19"><label>(17)</label><mml:math id="M21"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B5;</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mi>&#x003C3;</mml:mi><mml:msubsup><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B5;</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mi>&#x003C3;</mml:mi><mml:msubsup><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B5;</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mi>&#x003C3;</mml:mi><mml:msubsup><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msubsup><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>E</italic><sub><italic>b</italic></sub> is the direct-beam irradiance from the sun, <italic>E</italic><sub><italic>d</italic></sub> is the diffuse-skylight irradiance, <italic>E</italic><sub><italic>t</italic></sub> is the adjacent-terrain irradiance, and &#x003B1; is the surface albedo. For the long-wave radiation fluxes, &#x003B5;<sub><italic>a</italic></sub>, &#x003B5;<sub><italic>s</italic></sub> and &#x003B5;<sub><italic>t</italic></sub> are the emissivity for the air, glacier surface and adjacent terrain, respectively, &#x003C3; is the Stefan-Boltzmann constant, <italic>T</italic><sub>0</sub> is the air temperature which is a function of altitude, <italic>T</italic><sub><italic>s</italic></sub> is the glacier surface temperature, and <italic>T</italic><sub><italic>t</italic></sub> is the surface temperature of the surrounding terrain, which is responsible for the long-wave adjacent-terrain irradiance.</p>
<p>The energy balance at the debris/ice interface can be written as (Nakawo and Young, <xref ref-type="bibr" rid="B74">1981</xref>):</p>
<disp-formula id="E20"><label>(18)</label><mml:math id="M22"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msubsup><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02193;</mml:mi></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msubsup><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>Q</italic><sub><italic>m</italic></sub> is the heat flux used for sub-debris ice ablation, <inline-formula><mml:math id="M23"><mml:msubsup><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02193;</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the conductive heat flux from the debris, and <inline-formula><mml:math id="M24"><mml:msubsup><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the heat flux toward the ice that is not used for ablation, which is often negligible under a temperate ice assumption (Cuffey and Paterson, <xref ref-type="bibr" rid="B34">2010</xref>).</p>
<p>Most models assume constant heat storage and a linear debris-temperature gradient, then <italic>Q</italic><sub><italic>m</italic></sub> during the ablation season can be computed using the simplified heat-flux approach based on daily mean values at a minimum 24 h timescale (Nakawo and Young, <xref ref-type="bibr" rid="B74">1981</xref>; Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>):</p>
<disp-formula id="E21"><label>(19)</label><mml:math id="M25"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>k</italic><sub><italic>d</italic></sub> is the bulk thermal conductivity of the debris layer accounting for water in pore-spaces, <italic>T</italic><sub><italic>i</italic></sub> is the ice temperature. The sub-debris melt rate (<italic>M</italic><sub><italic>s</italic></sub>) can then be computed as (Nakawo and Young, <xref ref-type="bibr" rid="B74">1981</xref>; Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>):</p>
<disp-formula id="E22"><label>(20)</label><mml:math id="M26"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003C1;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>L</italic><sub><italic>f</italic></sub> is the latent heat of fusion for ice, and &#x003C1;<sub><italic>i</italic></sub> is the density of ice.</p>
<p>Precipitation can also contribute to surface melt (Fujita et al., <xref ref-type="bibr" rid="B42">2014</xref>; Fyffe et al., <xref ref-type="bibr" rid="B43">2014</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>), and several methods have been used to account for this process (e.g., Reid and Brock, <xref ref-type="bibr" rid="B88">2010</xref>; Fujita et al., <xref ref-type="bibr" rid="B42">2014</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>). Although more complicated processes have not been accounted for, such as the energy exchange between rainfall and heated debris that could eventually lead to additional ice melt. The latent and sensible heat fluxes have been found to be less significant with respect to melt variability compared to the shortwave and longwave fluxes in different climates (Sicart et al., <xref ref-type="bibr" rid="B101">2008</xref>).</p>
<p>Using the aforementioned equations, <xref ref-type="fig" rid="F6">Figure 6A</xref> shows the simulated ablation rates on the Baltoro Glacier in the Karakoram, using surface conditions presented by Mihalcea et al. (<xref ref-type="bibr" rid="B69">2008</xref>). Note that the modeled values are reasonably similar to field measurements on several Himalayan glaciers (Nicholson and Benn, <xref ref-type="bibr" rid="B75">2006</xref>), including Barpu Glacier (Khan, <xref ref-type="bibr" rid="B63">1989</xref>), Khumbu Glacier (Kayastha et al., <xref ref-type="bibr" rid="B62">2000</xref>), and Rakhiot Glacier (Mattson, <xref ref-type="bibr" rid="B66">1993</xref>). <xref ref-type="fig" rid="F6">Figure 6B</xref> is an example of the temporal variation in ablation rate as a function of debris thickness over an ablation season, which highlights the important role of the dynamics involving debris redistribution in regulating local ablation rates on the glacier surface, especially for areas with water bodies and ice-cliffs. <xref ref-type="fig" rid="F7">Figure 7</xref> depicts simulated surface ablation rates compared to the remote-sensing-based estimates (Mihalcea et al., <xref ref-type="bibr" rid="B69">2008</xref>). Both results show suppressed ablation in the terminus region and higher ablation around inter-moraine valleys corresponding to the difference in debris thickness (<xref ref-type="fig" rid="F5">Figure 5A</xref>).</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>The debris-thickness control on surface ablation rates. <bold>(A)</bold> June-averaged ablation rate as a function of debris thickness on the Baltoro Glacier (red) compared to <italic>in-situ</italic> measurements in the ablation season on other Himalayan glaciers (replotted), including Barpu Glacier (Khan, <xref ref-type="bibr" rid="B63">1989</xref>), Khumbu Glacier (Kayastha et al., <xref ref-type="bibr" rid="B62">2000</xref>), and Rakhiot Glacier (Mattson, <xref ref-type="bibr" rid="B66">1993</xref>). Note that the melt-enhancing effect of thin-debris is observed on all glaciers. <bold>(B)</bold> An example of the temporal variation in ablation rate vs. debris thickness around a supraglacial lake on the Baltoro Glacier. The debris layer is getting thinner as sediments are gradually sliding into the lake, which leads to an increase in ablation rate (the decrease of ablation rate in late ablation season is due to the decreasing irradiance).</p></caption>
<graphic xlink:href="feart-09-652279-g0006.tif"/>
</fig>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Surface ablation rates for the Baltoro Glacier modeled with two different methods: <bold>(A)</bold> The remote-sensing based approach (Mihalcea et al., <xref ref-type="bibr" rid="B69">2008</xref>), which uses surface temperature data acquired by satellite. <bold>(B)</bold> The surface energy balance based approach, which uses modeled surface irradiance. Both results are based on the same debris thickness distribution, assuming no water bodies and were produced to represent rates on August 14 2004. Note the high spatial variability in ablation rate caused by the debris thickness distribution (<xref ref-type="fig" rid="F5">Figure 5A</xref>).</p></caption>
<graphic xlink:href="feart-09-652279-g0007.tif"/>
</fig>
</sec>
<sec id="s5">
<title>5. Supraglacial Water Bodies and Ice Cliffs</title>
<p>Supraglacial water bodies and ice-cliffs elevate the ice loss of DCGs (Anderson, <xref ref-type="bibr" rid="B2">2014</xref>; Reid and Brock, <xref ref-type="bibr" rid="B86">2014</xref>; Steiner et al., <xref ref-type="bibr" rid="B103">2015</xref>; Buri et al., <xref ref-type="bibr" rid="B26">2016</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>, <xref ref-type="bibr" rid="B72">2018</xref>; Thompson et al., <xref ref-type="bibr" rid="B106">2016</xref>; Dobreva et al., <xref ref-type="bibr" rid="B35">2017</xref>; Mertes et al., <xref ref-type="bibr" rid="B68">2017</xref>; Huang et al., <xref ref-type="bibr" rid="B53">2018</xref>). The number of supraglacial ponds and lakes have been found to be increasing on many DCGs, which is related to the glacier morphological changes such as differential surface lowering and collapse of englacial channel roofs (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>, <xref ref-type="bibr" rid="B95">2002</xref>; Benn et al., <xref ref-type="bibr" rid="B10">2001</xref>; Quincey and Glasser, <xref ref-type="bibr" rid="B85">2009</xref>; Gibson et al., <xref ref-type="bibr" rid="B46">2017</xref>). These water bodies also play an important role in regulating water storage and drainage in a glacier system (Cuffey and Paterson, <xref ref-type="bibr" rid="B34">2010</xref>; Benn and Evans, <xref ref-type="bibr" rid="B9">2014</xref>), and cause englacial ablation via efficient heat transfer (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>; Gulley and Benn, <xref ref-type="bibr" rid="B48">2007</xref>; Benn et al., <xref ref-type="bibr" rid="B8">2012</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>; Mertes et al., <xref ref-type="bibr" rid="B68">2017</xref>).</p>
<p>Supraglacial lakes can be relatively large to moderate-sized water bodies (lake azimuthal distances have been observed to be up to 0.6 km in length on the Baltoro Glacier) that can be relatively short lived or exist over multiple years, and supraglacial ponds are smaller but more common water bodies that usually only appear during the ablation season. These water bodies are typically abundant in the ablation zone where the surface slope is gentle and heterogeneous surface lowering occurs (Reynolds, <xref ref-type="bibr" rid="B89">2000</xref>; Sakai et al., <xref ref-type="bibr" rid="B95">2002</xref>; Reid and Brock, <xref ref-type="bibr" rid="B86">2014</xref>). The evolution of many supraglacial lakes and ponds are governed by filling and draining cycles, and hydrological connectivity that are controlled by glacier surface and internal structures (Benn et al., <xref ref-type="bibr" rid="B10">2001</xref>, <xref ref-type="bibr" rid="B13">2017</xref>; Wessels et al., <xref ref-type="bibr" rid="B109">2002</xref>; Miles et al., <xref ref-type="bibr" rid="B71">2017</xref>). Studies also found that proglacial lakes also have a significant control on the glacier dynamics (Sutherland et al., <xref ref-type="bibr" rid="B104">2020</xref>). These processes have only been studied on a handful of glaciers in the field, and there are still many processes and feedbacks that need to be investigated.</p>
<p>Supraglacial ponds and lakes are efficient in transferring heat into glacier ice due to their low surface albedo and active convection (L&#x000FC;thje et al., <xref ref-type="bibr" rid="B65">2006</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>). Studies have estimated that the ablation rates around the lakes can be much higher than that of most debris-covered areas (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>; Thompson et al., <xref ref-type="bibr" rid="B106">2016</xref>). The simulations by Miles et al. (<xref ref-type="bibr" rid="B72">2018</xref>) revealed that ponds may be responsible for 1/8 of total ice loss in the Langtang valley, Nepal. Furthermore, supraglacial water bodies will most likely exhibit accelerated growth on Himalayan glaciers given current atmospheric temperature trends (Benn et al., <xref ref-type="bibr" rid="B10">2001</xref>). For example, Gibson et al. (<xref ref-type="bibr" rid="B46">2017</xref>) estimated that the number of water bodies on the Baltoro Glacier has increased from 234 in 2001 to 570 in 2012, and the area has tripled, which suggests that supraglacial water bodies can expand quickly on DCGs.</p>
<p>One of the first energy-balance models for supraglacial ponds on DCGs was developed by Sakai et al. (<xref ref-type="bibr" rid="B96">2000</xref>), which was used to address field observations on the Lirung Glacier in the Langtang Valley, Nepal. The core of this model is the energy-balance of a supraglacial lake that accounts for heat flux input and output due to meltwater, heat storage of the lake, and the bare-ice vs. debris covered areas beneath the lake surface:</p>
<disp-formula id="E23"><label>(21)</label><mml:math id="M27"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mi>Q</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>o</mml:mi><mml:mi>u</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mo>&#x00394;</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>Q</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>Q</italic> is the net radiation heat flux on the lake surface, <italic>Q</italic><sub><italic>in</italic></sub> is the input heat flux from meltwater inflow into the lake, which in most cases, can be neglected (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>), <italic>Q</italic><sub><italic>out</italic></sub> is the output heat flux from the water outflow from the lake, which dominates the englacial ablation such as the expansion of conduits (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>), &#x00394;<italic>Q</italic><sub><italic>t</italic></sub> is the change in heat storage of the lake, <italic>Q</italic><sub><italic>i</italic></sub> is the heat flux for ice melt at subaqueous ice-cliff, and <italic>Q</italic><sub><italic>d</italic></sub> is the heat flux for ice melt beneath the subaqueous debris layer.</p>
<p>A more recent model was presented by Miles et al. (<xref ref-type="bibr" rid="B70">2016</xref>), which improves the Sakai et al. (<xref ref-type="bibr" rid="B96">2000</xref>)&#x00027;s model in several aspects, such as the energy flux due to rainfall, more accurate turbulent fluxes, and the free convection of water that may circulate within the pond. In addition to the energy-balance, Miles et al. (<xref ref-type="bibr" rid="B70">2016</xref>) used a mass balance model to account for the changes in pond volume (&#x00394;<italic>V</italic>):</p>
<disp-formula id="E24"><label>(22)</label><mml:math id="M28"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mo>&#x00394;</mml:mo><mml:mi>V</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>l</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mi>o</mml:mi><mml:mi>u</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>V</italic><sub><italic>in</italic></sub> and <italic>V</italic><sub><italic>out</italic></sub> are the inflows and outflows, <italic>V</italic><sub><italic>i</italic></sub> is the subaqueous bare ice melt volume, <italic>V</italic><sub><italic>d</italic></sub> is the subaqueous subdebris melt volume, <italic>V</italic><sub><italic>le</italic></sub> is the volume of vapor exchanges, and <italic>V</italic><sub><italic>r</italic></sub> is the volume gained from rain. All terms in the energy-balance model and the mass-balance model have been discussed. Although some parameters are based on empirical relationships, this work provides the most comprehensive physics-based pond model to date. <xref ref-type="fig" rid="F8">Figure 8</xref> depicts the large number of supraglacial water bodies identified in the terminus region of the Baltoro Glacier. Our simulation results (<xref ref-type="fig" rid="F8">Figures 8C,D</xref>) based on the above model suggest that supraglacial water bodies make a significant contribution to the total ice loss and surface lowering in the ablation zone, and exhibit a non-linear seasonal trend toward the end of the ablation season, as the extent and amount of lakes and ponds gradually expand to yearly maximum.</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p><bold>(A)</bold> ASTER false-color composite imagery (VNIR bands) of the lower ablation zone of the Baltoro Glacier acquired on August 14, 2004. <bold>(B)</bold> A large number of supraglacial ponds and lakes mapped from the satellite data using the normalized difference water index. <bold>(C)</bold> Temporal variation in the simulated cumulative ice-volume loss for this area over the ablation season of 2004. Two scenarios (with and without lakes/ponds) are compared. <bold>(D)</bold> Comparison of the mean surface altitude variations for the two simulations. Note the increased non-linearity in ice loss and surface lowering when supraglacial lakes/ponds are present.</p></caption>
<graphic xlink:href="feart-09-652279-g0008.tif"/>
</fig>
<p>The dynamics of supraglacial lakes and ponds is often coupled with the evolution of ice-cliffs around them. Studies have suggested that this cliff/lake system plays an important role in governing the mass balance of DCGs, and most field observations and modeling results confirmed that ice-cliffs can make a significant contribution to the total ice mass loss on a DCG (Sakai et al., <xref ref-type="bibr" rid="B95">2002</xref>; Anderson, <xref ref-type="bibr" rid="B2">2014</xref>; Steiner et al., <xref ref-type="bibr" rid="B103">2015</xref>; Buri et al., <xref ref-type="bibr" rid="B26">2016</xref>), although many questions remain unresolved (Sakai et al., <xref ref-type="bibr" rid="B95">2002</xref>; Buri et al., <xref ref-type="bibr" rid="B26">2016</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>). Several forms of interactions exist between ice-cliffs and ponds/lakes. For example, lake undercut can cause the steepening of some ice-cliffs (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>, <xref ref-type="bibr" rid="B95">2002</xref>), and the melt of ice-cliff below the water surface is strongly affected by water fluctuations (Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>).</p>
<p>Sakai et al. (<xref ref-type="bibr" rid="B95">2002</xref>) found that the orientation and inclination of ice-cliffs significantly control the evolution of ice cliffs given the solar geometry and cast shadows. The geometry of the ice-cliff also controls the thickness distribution of debris cover which regulates meltwater production. Ice-cliffs are also controlled by englacial processes, and studies have attributed the formation of some ice-cliffs to the collapses of englacial conduits (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>, <xref ref-type="bibr" rid="B95">2002</xref>). Some recent simulations (Steiner et al., <xref ref-type="bibr" rid="B103">2015</xref>; Buri et al., <xref ref-type="bibr" rid="B26">2016</xref>) confirmed that the complex terrain around the cliff has a non-negligible effect on ice-cliff evolution due to the local shading and the adjacent terrain irradiance. Furthermore, ice-cliff evolution can also be governed by pond-water ablation that undercuts adjacent lake slopes. In addition, the ice topographic load coupled with the debris load causes significant variations in the ice stress fields, that partially control the process of calving. We have observed this process on a very large supraglacial lake near the terminus of the Liligo Glacier during the summer of 2005. Significant water depths coupled with wind and wave action can cause significant undercutting, thereby altering the ice stress fields which enables calving and the production of icebergs in glacial lakes. The ice cliffs that are generated from these processes typically have very steep slopes, and the debris loads above them are quickly transported. Collectively, various feedbacks promote ice-cliff retreat and rapid supraglacial lake expansion, that most likely represents a non-linear ablation response to climate forcing.</p>
<p>Understanding supraglacial water bodies and ice-cliff dynamics is critical to assessing DCG sensitivity to climate change, and valuable knowledge has been gained from the aforementioned studies. However, there are still many issues that need to be investigated including the filling and draining cycles that control the evolution of many supraglacial ponds, the connectivity between surface water bodies and englacial conduits, englacial ablation due to warm water outflow, rockfall from ice-cliffs, and modeling the edge effects and temperature gradients in ponds (Benn et al., <xref ref-type="bibr" rid="B10">2001</xref>, <xref ref-type="bibr" rid="B13">2017</xref>; Miles et al., <xref ref-type="bibr" rid="B70">2016</xref>, <xref ref-type="bibr" rid="B71">2017</xref>).</p>
</sec>
<sec id="s6">
<title>6. Feedback and System Couplings</title>
<p>The dynamics of a DCG system is controlled by the complex interactions between climate, debris load, ice dynamics, hydrologic conditions and topography (Scherler et al., <xref ref-type="bibr" rid="B98">2011a</xref>,<xref ref-type="bibr" rid="B99">b</xref>; Dobreva et al., <xref ref-type="bibr" rid="B35">2017</xref>; Gibson et al., <xref ref-type="bibr" rid="B46">2017</xref>). Unfortunately, few studies have investigated the couplings and feedback mechanisms in a DCG system, which is a major issue in characterizing DCGs, because isolating a small number of processes out of the full system increases uncertainty in results, as system couplings and feedbacks can significantly affect glacier dynamics (Rowan et al., <xref ref-type="bibr" rid="B94">2015</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>). Studies are often limited in spatio-temporal scale (Anderson, <xref ref-type="bibr" rid="B2">2014</xref>; Anderson et al., <xref ref-type="bibr" rid="B5">2021</xref>), because they rely on instantaneous remote-sensing data or field measurements at limited sites that do not describe multi-year or longer-term variations in glacier conditions (such as ice thickness, debris distribution and topography). Furthermore, many system couplings and feedbacks operate at very different spatio-temporal scales, making them difficult to investigate in the field. Therefore, numerical modeling is required to address the feedbacks and non-linear relationships between major glacier components. <xref ref-type="fig" rid="F9">Figure 9</xref> is a conceptual diagram showing the major systems couplings we identified between climate, mass balance, supraglacial water bodies, debris load, ice flow, and topography that govern the dynamics of a DCG.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p>Diagram illustrating the interactions between key controlling factors discussed in this paper that govern the dynamics of a debris-covered glacier. The feedbacks (dashed lines) are important to the system as they may increase a glacier&#x00027;s sensitivity to climate change over time.</p></caption>
<graphic xlink:href="feart-09-652279-g0009.tif"/>
</fig>
<sec>
<title>6.1. Feedbacks on Glacier Surface</title>
<p><xref ref-type="fig" rid="F10">Figure 10</xref> is a conceptual diagram that illustrates the feedbacks we identified between ablation, surface ponding and debris thickness. The most obvious positive feedback describes the relationship between meltwater production and supraglacial pond expansion: as a pond grows, its surface area and surrounding ice-cliff area get larger, they absorb more solar energy due to the lower albedo, causing enhanced melt and calving, which in turn, leads to further expansion of the pond. Another feedback involves the differential lowering of the glacier surface, the distribution of debris and the increased ponding of surface water: the spatially heterogeneous distribution of supraglacial water bodies and ice-cliffs governs the differential elevation change, which in turn, facilities the expansion of existing ponds and the formation of new ponds as suggested by multiple studies (Benn and Owen, <xref ref-type="bibr" rid="B12">2002</xref>; Benn et al., <xref ref-type="bibr" rid="B8">2012</xref>; Reid and Brock, <xref ref-type="bibr" rid="B86">2014</xref>; Rowan et al., <xref ref-type="bibr" rid="B94">2015</xref>; Salerno et al., <xref ref-type="bibr" rid="B97">2017</xref>; Miles et al., <xref ref-type="bibr" rid="B72">2018</xref>). The presence of supraglacial ponds and lakes also increases the spatial heterogeneity of debris thickness. Rounce et al. (<xref ref-type="bibr" rid="B92">2018</xref>) found that debris thinning is usually associated with the formation of supraglacial lakes and ice-cliffs, then the heterogeneous distribution of debris thickness further enhances the differential surface lowering which creates more depression zones for new ponds.</p>
<fig id="F10" position="float">
<label>Figure 10</label>
<caption><p>Conceptual diagram that illustrates the positive feedbacks between ablation, surface ponding and the spatial distribution of debris thickness.</p></caption>
<graphic xlink:href="feart-09-652279-g0010.tif"/>
</fig>
<p>Positive feedbacks also exist between the lowing of surface albedo and enhanced melt, as a wet surface often exhibits lower albedo. Pritchard et al. (<xref ref-type="bibr" rid="B81">2008</xref>) showed that this mechanism has significant influence on millennial-scale mass-loss of ice sheets. For certain locations on the glacier surface, such as on ice-cliffs, the decrease in albedo is more significant due to the thinner debris and high moisture content, i.e., the &#x0201C;dirty ice&#x0201D; that causes more ablation (Fyffe et al., <xref ref-type="bibr" rid="B44">2020</xref>), which is also coupled with the topographic factors (such as the orientation of the cliff wall and the high adjacent terrain irradiance on the cliffs) to form positive feedback that accelerates melting, causing faster ice-cliff backwasting than previously thought (Sakai et al., <xref ref-type="bibr" rid="B96">2000</xref>; Benn and Owen, <xref ref-type="bibr" rid="B12">2002</xref>; Reid and Brock, <xref ref-type="bibr" rid="B86">2014</xref>).</p>
</sec>
<sec>
<title>6.2. The Role of Glacier Surface Topography</title>
<p>Surface topography is an important factor that acts as a bridge between multiple components of a DCG, and studies have identified the following topographic effects on DCGs:</p>
<p>1. Surface topography has a direct impact on glacier surface ablation because slope, slope azimuth, and basin topographic shielding directly control the total amount of shortwave radiation received by the glacier surface (Arnold et al., <xref ref-type="bibr" rid="B7">2006</xref>; Garg et al., <xref ref-type="bibr" rid="B45">2017</xref>; Olson and Rupper, <xref ref-type="bibr" rid="B77">2019</xref>).</p>
<p>2. Surface topographic conditions control supraglacial hydrology, especially the development of ponds and lakes (Huo et al., <xref ref-type="bibr" rid="B55">2021</xref>). Water transport is enhanced with steeper gradients, while ponding is enhanced with lower gradients, therefore supraglacial water bodies are more abundant where the glacier surface is flatter (Reynolds, <xref ref-type="bibr" rid="B89">2000</xref>; Sakai et al., <xref ref-type="bibr" rid="B95">2002</xref>; Immerzeel et al., <xref ref-type="bibr" rid="B56">2014</xref>; Reid and Brock, <xref ref-type="bibr" rid="B86">2014</xref>). Topographic depressions caused by differential ablation and debris loading are also important, as they provide topographic sinks for meltwater to accumulate.</p>
<p>3. Topography also affects the distribution and transport of supraglacial debris. Observations and simulations have shown that the morphometric conditions of the sidewalls and headwalls of the glacial valley control sediment input to the glacier (Benn and Evans, <xref ref-type="bibr" rid="B9">2014</xref>; Anderson and Anderson, <xref ref-type="bibr" rid="B3">2016</xref>).</p>
<p>4. The adjacent-terrain irradiance due to complex topography coupled with debris cover contributes to the surface melt, and especially the retreat of ice-cliffs (Sakai et al., <xref ref-type="bibr" rid="B95">2002</xref>; Wessels et al., <xref ref-type="bibr" rid="B109">2002</xref>).</p>
</sec>
<sec>
<title>6.3. Basal Processes</title>
<p>The glacier bed is slowly evolving due to differential glacier erosion caused by variations in ice-flow dynamics, such that erosion is most likely at a maximum near the ELA due to relatively thick ice and high flow velocities (Herman et al., <xref ref-type="bibr" rid="B50">2011</xref>; Steer et al., <xref ref-type="bibr" rid="B102">2012</xref>). Another complexity is the basal hydrological conditions that govern sliding velocity, as high basal water pressure may lubricate the ice-bedrock interface (Cuffey and Paterson, <xref ref-type="bibr" rid="B34">2010</xref>; Quincey et al., <xref ref-type="bibr" rid="B83">2011</xref>). The higher sliding velocity leads to more frictional heating that further increases basal water pressure, which forms a positive feedback (Benn and Evans, <xref ref-type="bibr" rid="B9">2014</xref>). Many details about the interactions between basal sliding, erosion, subglacial ablation, basal hydrology, and subglacial debris transport need to be investigated in future studies.</p>
<p>Collectively, the DCG system is complicated by interactions and feedbacks between debris load, ablation, topography, ice-flow, basal processes, and supraglacial water bodies. An integrated modeling of these processes may help explain field observations over High-Mountain Asia, such as the rapid lowering on Himalayan DCGs (K&#x000E4;&#x000E4;b et al., <xref ref-type="bibr" rid="B61">2012</xref>; Immerzeel et al., <xref ref-type="bibr" rid="B57">2013</xref>; Fujita et al., <xref ref-type="bibr" rid="B42">2014</xref>) and the advancement of many glaciers in the Karakoram (Bolch et al., <xref ref-type="bibr" rid="B22">2012</xref>; Bishop et al., <xref ref-type="bibr" rid="B19">2014</xref>; Dobreva et al., <xref ref-type="bibr" rid="B35">2017</xref>; Farinotti et al., <xref ref-type="bibr" rid="B39">2020</xref>). Furthermore, it is also important to account for ice-flow dynamics being a potential cause of glacier thinning, as the re-distribution of ice mass, especially the declining ice discharge plays a dominant role in setting the thinning patterns on debris-free glaciers (e.g., Cuffey and Paterson, <xref ref-type="bibr" rid="B34">2010</xref>), and has also been shown to be important for DCGs (Vincent et al., <xref ref-type="bibr" rid="B108">2016</xref>; Brun et al., <xref ref-type="bibr" rid="B24">2018</xref>; Anderson et al., <xref ref-type="bibr" rid="B5">2021</xref>).</p>
<p>Improved parameterization schemes (more parameters and processes) suggest that climate glacier dynamics for DCGs are more complex than debris-free glaciers and that DCGs may be more sensitive to climate change than previously thought.</p>
<p>We recommend that more energy input components (e.g., the adjacent-terrain irradiance) and multi-scale topographic effects be accounted for in improved parameterization schemes. Our simulations indicate that scale dependent processes such as the gravitational debris transport and supraglacial ponding must also be accounted for given the coupling of climate and surface process and the non-linear responses in ablation and ice-mass loss. In all likelihood, our simulations probably underestimate the magnitude of non-linear responses due to these coupled processes/systems, as other processes such as calving were not accounted for. Furthermore, it is essential that future research focus on addressing questions about the significance of missing and existing processes so that we have an improved understanding of the degree of dominance and partitioning of key parameters and processes/systems that need to be incorporated into DCG models. This will require important comparative analysis of parameterization schemes and parameter sensitivity analysis to account for uncertainty in so many parameters and processes that govern climate-glacier dynamics.</p>
<p>More remote sensing and field-based studies are required to validate numerical modeling results. Further research is also required to determine if DCGs in different geographic locations exhibit similar or different sensitivity to climate change given variations in internal and external forcing factors. We speculate that such complex non-linear systems have the potential to reach tipping points based upon changing climate and geomorphological conditions within glacierized basins.</p>
</sec>
</sec>
<sec sec-type="conclusions" id="s7">
<title>7. Conclusions</title>
<p>Debris-covered glaciers represent complex non-linear systems that operate given the coupling of climate, geomorphological and glaciological processes and feedbacks. Many important processes and systems have not been adequately characterized. Consequently, our current understanding of DCG sensitivity to climate change is most likely oversimplified, as numerous processes and scale dependencies have not been accounted for. This has resulted in conflicting views about these glaciers, such as the Karakoram anomaly and the &#x0201C;debris-covered glacier anomaly,&#x0201D; where the lack of data and information about climate forcing and ablation dynamics force us to rely on empirical measurements that represent a static snapshot in space and time. This paper and our simulations have attempted to provide a more integrated synthesis into the complexity of DCG systems by addressing important concepts associated with the key external and internal forcing factors and glacial processes, highlighting the importance of feedback mechanisms, system couplings, and recognition of the high degree of uncertainty related to numerous properties and parameters that need to be taken into consideration. Our simulation results based on the Baltoro Glacier in the central Karakoram combined with a review of literature reveal that climate-DCG dynamics and responses can be significantly attributed to the following factors:</p>
<p>1. Climate forcing: solar radiation and precipitation are the main driving forces for ice loss or gain. An accurate surface energy balance model for debris-covered glaciers needs to account for short-wave and long-wave irradiance. Our simulations showed that multi-scale topographic effects on surface ablation can be significant for debris-covered glaciers, given glacier surface topographic evolution and differential sky-view conditions in the glacier valley. Other meso-scale topographic effects such as cast shadows and adjacent-terrain influences should also be considered.</p>
<p>2. Supraglacial debris: physical properties of debris, such as thickness, lithology, grain size, moisture content and thermal conductivity strongly govern the surface ablation rates on debris-covered glaciers. Many of these factors, however, have not been adequately characterized or studied. We showed that there are significant spatial variability in debris thickness, albedo, particle size and composition, which are also coupled with the dynamic movement of debris due to gravity and ice-flow. Such complexity needs to be accounted for in models and field studies to better assess debris-related effects on glacier mass balance.</p>
<p>3. Supraglacial water bodies and ice-cliffs: supraglacial ponds, lakes, and surrounding ice-cliffs are melt hotspots on debris-covered glaciers and are expanding in size on many glaciers. Our simulations based on the Baltoro Glacier suggest that supraglacial water bodies make a significant contribution to the total ice-mass loss over the ablation season. Further investigations are needed to address the evolution of supraglacial water bodies and ice-cliffs, and their interactions with debris load, topography and englacial processes.</p>
<p>4. System coupling and feedback: positive feedbacks on debris-covered glaciers have been identified in numerical simulations, such as the surface melt-lowering-ponding feedback, and the atmospheric convection feedback. These couplings and especially the positive feedbacks may govern the non-linearity of the glacier&#x00027;s response to climate forcing, which are represented as an acceleration in ice-mass loss and heterogeneous surface lowering. The combined effect of these processes may lead to the beginning of a critical transition of the glacier system that signifies an increasing level of glacier sensitivity to climate change.</p>
</sec>
<sec sec-type="data-availability-statement" id="s8">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.</p>
</sec>
<sec id="s9">
<title>Author Contributions</title>
<p>DH conducted the model simulations, made the figures, and wrote the sections 3.2&#x02013;6. MB proposed and designed the research, developed the parameterization schemes, wrote the sections 2.1 and 3.1, and edited the manuscript. AB wrote the section 2.2 and edited the manuscript. All authors contributed to the writing of the introduction and conclusions and the final editing.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</body>
<back>
<ack><p>We would like to thank Dr. Zhaohui Chi and Brennan Young at Texas A&#x00026;M University for discussion. AB recognizes support from the Canadian Natural Sciences and Engineering Research Council and the computational resources of Compute Canada.</p></ack>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Adhikary</surname> <given-names>S.</given-names></name> <name><surname>Masayoshi</surname> <given-names>N.</given-names></name> <name><surname>Seko</surname> <given-names>K.</given-names></name> <name><surname>Shakya</surname> <given-names>B.</given-names></name></person-group> (<year>2000</year>). <article-title>&#x0201C;Dust influence on the melting process of glacier ice: experimental results from Lirung glacier, Nepal Himalayas</article-title>, in <source>Debris-Covered Glaciers: Proceedings of an International Workshop</source> (<publisher-loc>Seattle, WA</publisher-loc>: <publisher-name>IAHS</publisher-name>), <fpage>43</fpage>.</citation></ref>
<ref id="B2">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>L. S.</given-names></name></person-group> (<year>2014</year>). <source>Glacier response to climate change: modeling the effects of weather and debris-cover</source> (Doctoral dissertation), University of Colorado at Boulder, Boulder, CO.</citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>L. S.</given-names></name> <name><surname>Anderson</surname> <given-names>R. S.</given-names></name></person-group> (<year>2016</year>). <article-title>Modeling debris-covered glaciers: response to steady debris deposition</article-title>. <source>Cryosphere</source> <volume>10</volume>:<fpage>1105</fpage>. <pub-id pub-id-type="doi">10.5194/tc-10-1105-2016</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>L. S.</given-names></name> <name><surname>Anderson</surname> <given-names>R. S.</given-names></name></person-group> (<year>2018</year>). <article-title>Debris thickness patterns on debris-covered glaciers</article-title>. <source>Geomorphology</source> <volume>311</volume>, <fpage>1</fpage>&#x02013;<lpage>12</lpage>. <pub-id pub-id-type="doi">10.1016/j.geomorph.2018.03.014</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>L. S.</given-names></name> <name><surname>Armstrong</surname> <given-names>W. H.</given-names></name> <name><surname>Anderson</surname> <given-names>R. S.</given-names></name> <name><surname>Buri</surname> <given-names>P.</given-names></name></person-group> (<year>2021</year>). <article-title>Debris cover and the thinning of Kennicott glacier, Alaska: <italic>in situ</italic> measurements, automated ice cliff delineation and distributed melt estimates</article-title>. <source>Cryosphere</source> <volume>15</volume>, <fpage>265</fpage>&#x02013;<lpage>282</lpage>. <pub-id pub-id-type="doi">10.5194/tc-15-265-2021</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname> <given-names>R. S.</given-names></name></person-group> (<year>2000</year>). <article-title>A model of ablation-dominated medial moraines and the generation of debris-mantled glacier snouts</article-title>. <source>J. Glaciol</source>. <volume>46</volume>, <fpage>459</fpage>&#x02013;<lpage>469</lpage>. <pub-id pub-id-type="doi">10.3189/172756500781833025</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Arnold</surname> <given-names>N. S.</given-names></name> <name><surname>Rees</surname> <given-names>W. G.</given-names></name> <name><surname>Hodson</surname> <given-names>A. J.</given-names></name> <name><surname>Kohler</surname> <given-names>J.</given-names></name></person-group> (<year>2006</year>). <article-title>Topographic controls on the surface energy balance of a high arctic valley glacier</article-title>. <source>J. Geophys. Res. Earth Surface</source> <fpage>111</fpage>. <pub-id pub-id-type="doi">10.1029/2005JF000426</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benn</surname> <given-names>D.</given-names></name> <name><surname>Bolch</surname> <given-names>T.</given-names></name> <name><surname>Hands</surname> <given-names>K.</given-names></name> <name><surname>Gulley</surname> <given-names>J.</given-names></name> <name><surname>Luckman</surname> <given-names>A.</given-names></name> <name><surname>Nicholson</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Response of debris-covered glaciers in the Mount Everest region to recent warming, and implications for outburst flood hazards</article-title>. <source>Earth Sci. Rev</source>. <volume>114</volume>, <fpage>156</fpage>&#x02013;<lpage>174</lpage>. <pub-id pub-id-type="doi">10.1016/j.earscirev.2012.03.008</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Benn</surname> <given-names>D.</given-names></name> <name><surname>Evans</surname> <given-names>D. J.</given-names></name></person-group> (<year>2014</year>). <source>Glaciers and Glaciation</source>. <publisher-loc>London</publisher-loc>: <publisher-name>Routledge</publisher-name>. <pub-id pub-id-type="doi">10.4324/9780203785010</pub-id><pub-id pub-id-type="pmid">29134193</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benn</surname> <given-names>D.</given-names></name> <name><surname>Wiseman</surname> <given-names>S.</given-names></name> <name><surname>Hands</surname> <given-names>K.</given-names></name></person-group> (<year>2001</year>). <article-title>Growth and drainage of supraglacial lakes on debris-mantled Ngozumpa glacier, Khumbu Himal, Nepal</article-title>. <source>J. Glaciol</source>. <volume>47</volume>, <fpage>626</fpage>&#x02013;<lpage>638</lpage>. <pub-id pub-id-type="doi">10.3189/172756501781831729</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benn</surname> <given-names>D. I.</given-names></name> <name><surname>Lehmkuhl</surname> <given-names>F.</given-names></name></person-group> (<year>2000</year>). <article-title>Mass balance and equilibrium-line altitudes of glaciers in high-mountain environments</article-title>. <source>Quat. Int</source>. <volume>65</volume>, <fpage>15</fpage>&#x02013;<lpage>29</lpage>. <pub-id pub-id-type="doi">10.1016/S1040-6182(99)00034-8</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benn</surname> <given-names>D. I.</given-names></name> <name><surname>Owen</surname> <given-names>L. A.</given-names></name></person-group> (<year>2002</year>). <article-title>Himalayan glacial sedimentary environments: a framework for reconstructing and dating the former extent of glaciers in high mountains</article-title>. <source>Quat. Int</source>. <volume>97</volume>, <fpage>3</fpage>&#x02013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1016/S1040-6182(02)00048-4</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benn</surname> <given-names>D. I.</given-names></name> <name><surname>Thompson</surname> <given-names>S.</given-names></name> <name><surname>Gulley</surname> <given-names>J.</given-names></name> <name><surname>Mertes</surname> <given-names>J.</given-names></name> <name><surname>Luckman</surname> <given-names>A.</given-names></name> <name><surname>Nicholson</surname> <given-names>L.</given-names></name></person-group> (<year>2017</year>). <article-title>Structure and evolution of the drainage system of a Himalayan debris-covered glacier, and its relationship with patterns of mass loss</article-title>. <source>Cryosphere</source> <volume>11</volume>, <fpage>2247</fpage>&#x02013;<lpage>2264</lpage>. <pub-id pub-id-type="doi">10.5194/tc-11-2247-2017</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Berger</surname> <given-names>A.</given-names></name></person-group> (<year>1978</year>). <article-title>Long-term variations of daily insolation and quaternary climatic changes</article-title>. <source>J. Atmos. Sci</source>. <volume>35</volume>, <fpage>2362</fpage>&#x02013;<lpage>2367</lpage>. <pub-id pub-id-type="doi">10.1175/1520-0469(1978)035&#x0003C;2362:LTVODI&#x0003E;2.0.CO;2</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bird</surname> <given-names>R. E.</given-names></name> <name><surname>Riordan</surname> <given-names>C.</given-names></name></person-group> (<year>1986</year>). <article-title>Simple solar spectral model for direct and diffuse irradiance on horizontal and tilted planes at the Earth&#x00027;s surface for cloudless atmospheres</article-title>. <source>J. Climate Appl. Meteorol</source>. <volume>25</volume>, <fpage>87</fpage>&#x02013;<lpage>97</lpage>. <pub-id pub-id-type="doi">10.1175/1520-0450(1986)025&#x0003C;0087:SSSMFD&#x0003E;2.0.CO;2</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Bonk</surname> <given-names>R.</given-names></name> <name><surname>Kamp</surname> <given-names>U.</given-names> <suffix>Jr.</suffix></name> <name><surname>Shroder</surname> <given-names>J. F.</given-names> <suffix>Jr.</suffix></name></person-group> (<year>2001</year>). <article-title>Terrain analysis and data modeling for alpine glacier mapping</article-title>. <source>Polar Geogr</source>. <volume>25</volume>, <fpage>182</fpage>&#x02013;<lpage>201</lpage>. <pub-id pub-id-type="doi">10.1080/10889370109377712</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Colby</surname> <given-names>J. D.</given-names></name></person-group> (<year>2011</year>). <article-title>Topographic normalization of multispectral satellite imagery</article-title>. <source>J. Glaciol</source>. <volume>55</volume>, <fpage>131</fpage>&#x02013;<lpage>146</lpage>. Available online at: <ext-link ext-link-type="uri" xlink:href="http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.727.9345&#x00026;rep=rep1&#x00026;type=pdf">http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.727.9345&#x00026;rep=rep1&#x00026;type=pdf</ext-link></citation></ref>
<ref id="B18">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Dobreva</surname> <given-names>I. D.</given-names></name></person-group> (<year>2017</year>). <article-title>Geomorphometry and mountain geodynamics: issues of scale and complexity</article-title>, in <source>Integrating Scale in Remote Sensing and GIS</source>, eds <person-group person-group-type="editor"><name><surname>Quattrochi</surname> <given-names>D. A.</given-names></name> <name><surname>Wentz</surname> <given-names>E.</given-names></name> <name><surname>Lam</surname> <given-names>N. S. N.</given-names></name> <name><surname>Emerson</surname> <given-names>C. W.</given-names></name></person-group> (<publisher-loc>Boca Raton, FL</publisher-loc>: <publisher-name>Taylor and Francis</publisher-name>), <fpage>189</fpage>&#x02013;<lpage>228</lpage>. <pub-id pub-id-type="doi">10.1201/9781315373720-8</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Shroder</surname> <given-names>J. F.</given-names></name> <name><surname>Ali</surname> <given-names>G.</given-names></name> <name><surname>Bush</surname> <given-names>A. B. G.</given-names></name> <name><surname>Haritashya</surname> <given-names>U. K.</given-names></name> <name><surname>Roohi</surname> <given-names>R.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Remote sensing of glaciers in Afghanistan and Pakistan</article-title>, in <source>Global Land Ice Measurements from Space</source>, eds <person-group person-group-type="editor"><name><surname>Kargel</surname> <given-names>J. S.</given-names></name> <name><surname>Leonard</surname> <given-names>G. J.</given-names></name> <name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>K&#x000E4;&#x000E4;b</surname> <given-names>A.</given-names></name> <name><surname>Raup</surname> <given-names>B.</given-names></name></person-group> (<publisher-loc>Berlin; Heidelberg</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>509</fpage>&#x02013;<lpage>548</lpage>. <pub-id pub-id-type="doi">10.1007/978-3-540-79818-7_23</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Shroder</surname> <given-names>J. F.</given-names> <suffix>Jr.</suffix></name> <name><surname>Bonk</surname> <given-names>R.</given-names></name> <name><surname>Olsenholler</surname> <given-names>J.</given-names></name></person-group> (<year>2002</year>). <article-title>Geomorphic change in high mountains: a western himalayan perspective</article-title>. <source>Glob. Planet. Change</source> <volume>32</volume>, <fpage>311</fpage>&#x02013;<lpage>329</lpage>. <pub-id pub-id-type="doi">10.1016/S0921-8181(02)00073-5</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Young</surname> <given-names>B. W.</given-names></name> <name><surname>Colby</surname> <given-names>J. D.</given-names></name> <name><surname>Furfaro</surname> <given-names>R.</given-names></name> <name><surname>Schiassi</surname> <given-names>E.</given-names></name> <name><surname>Chi</surname> <given-names>Z.</given-names></name></person-group> (<year>2019</year>). <article-title>Theoretical evaluation of anisotropic reflectance correction approaches for addressing multi-scale topographic effects on the radiation-transfer cascade in mountain environments</article-title>. <source>Remote Sens</source>. <volume>11</volume>:<fpage>2728</fpage>. <pub-id pub-id-type="doi">10.3390/rs11232728</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bolch</surname> <given-names>T.</given-names></name> <name><surname>Kulkarni</surname> <given-names>A.</given-names></name> <name><surname>K&#x000E4;&#x000E4;b</surname> <given-names>A.</given-names></name> <name><surname>Huggel</surname> <given-names>C.</given-names></name> <name><surname>Paul</surname> <given-names>F.</given-names></name> <name><surname>Cogley</surname> <given-names>J. G.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>The state and fate of Himalayan glaciers</article-title>. <source>Science</source> <volume>336</volume>, <fpage>310</fpage>&#x02013;<lpage>314</lpage>. <pub-id pub-id-type="doi">10.1126/science.1215828</pub-id><pub-id pub-id-type="pmid">22517852</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Braithwaite</surname> <given-names>R. J.</given-names></name> <name><surname>Olesen</surname> <given-names>O. B.</given-names></name></person-group> (<year>1990</year>). <article-title>Response of the energy balance on the margin of the greenland ice sheet to temperature changes</article-title>. <source>J. Glaciol</source>. <volume>36</volume>, <fpage>217</fpage>&#x02013;<lpage>221</lpage>. <pub-id pub-id-type="doi">10.3189/S0022143000009461</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brun</surname> <given-names>F.</given-names></name> <name><surname>Wagnon</surname> <given-names>P.</given-names></name> <name><surname>Berthier</surname> <given-names>E.</given-names></name> <name><surname>Shea</surname> <given-names>J. M.</given-names></name> <name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name> <name><surname>Kraaijenbrink</surname> <given-names>P. D. A.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Ice cliff contribution to the tongue-wide ablation of Changri Nup glacier, Nepal, central Himalaya</article-title>. <source>Cryosphere</source> <volume>12</volume>, <fpage>3439</fpage>&#x02013;<lpage>3457</lpage>. <pub-id pub-id-type="doi">10.5194/tc-12-3439-2018</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bueler</surname> <given-names>E.</given-names></name> <name><surname>Brown</surname> <given-names>J.</given-names></name></person-group> (<year>2009</year>). <article-title>Shallow shelf approximation as a &#x0201C;sliding law&#x0201D; in a thermomechanically coupled ice sheet model</article-title>. <source>J. Geophys. Res. Earth Surface</source> <fpage>114</fpage>. <pub-id pub-id-type="doi">10.1029/2008JF001179</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Buri</surname> <given-names>P.</given-names></name> <name><surname>Miles</surname> <given-names>E. S.</given-names></name> <name><surname>Steiner</surname> <given-names>J. F.</given-names></name> <name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name> <name><surname>Wagnon</surname> <given-names>P.</given-names></name> <name><surname>Pellicciotti</surname> <given-names>F.</given-names></name></person-group> (<year>2016</year>). <article-title>A physically based 3-D model of ice cliff evolution over debris-covered glaciers</article-title>. <source>J. Geophys. Res. Earth Surface</source> <volume>121</volume>, <fpage>2471</fpage>&#x02013;<lpage>2493</lpage>. <pub-id pub-id-type="doi">10.1002/2016JF004039</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Bush</surname> <given-names>A. B. G.</given-names></name> <name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Huo</surname> <given-names>D.</given-names></name> <name><surname>Chi</surname> <given-names>Z.</given-names></name> <name><surname>Tiwari</surname> <given-names>U.</given-names></name></person-group> (<year>2020</year>). <article-title>Issues in climate analysis and modeling for understanding mountain erosion dynamics</article-title>, in <source>Reference Module in Earth Systems and Environmental Sciences</source> (<publisher-loc>Amsterdam</publisher-loc>: <publisher-name>Elsevier</publisher-name>). <pub-id pub-id-type="doi">10.1016/B978-0-12-818234-5.00022-5</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>H.</given-names></name> <name><surname>Lee</surname> <given-names>C.</given-names></name></person-group> (<year>2000</year>). <article-title>Numerical simulation of debris flows</article-title>. <source>Can. Geotech. J</source>. <volume>37</volume>, <fpage>146</fpage>&#x02013;<lpage>160</lpage>. <pub-id pub-id-type="doi">10.1139/t99-089</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Coakley</surname> <given-names>J.</given-names></name></person-group> (<year>2003</year>). <article-title>Reflectance and albedo, surface</article-title>, in <source>Encyclopedia of the Atmosphere Sciences</source> eds. <person-group person-group-type="editor"><name><surname>Holton</surname> <given-names>J.</given-names></name> <name><surname>Curry</surname> <given-names>J.</given-names></name></person-group> (<publisher-loc>Amsterdam</publisher-loc>: <publisher-name>Elsevier</publisher-name>), <fpage>1914</fpage>&#x02013;<lpage>1923</lpage>. <pub-id pub-id-type="doi">10.1016/B0-12-227090-8/00069-5</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Collier</surname> <given-names>E.</given-names></name> <name><surname>Maussion</surname> <given-names>F.</given-names></name> <name><surname>Nicholson</surname> <given-names>L. I.</given-names></name> <name><surname>M&#x000F6;lg</surname> <given-names>T.</given-names></name> <name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name> <name><surname>Bush</surname> <given-names>A. B. G.</given-names></name></person-group> (<year>2015</year>). <article-title>Impact of debris cover on glacier ablation and atmosphere-glacier feedbacks in the Karakoram</article-title>. <source>Cryosphere</source> <volume>9</volume>, <fpage>1617</fpage>&#x02013;<lpage>1632</lpage>. <pub-id pub-id-type="doi">10.5194/tc-9-1617-2015</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Collier</surname> <given-names>E.</given-names></name> <name><surname>M&#x000F6;lg</surname> <given-names>T.</given-names></name> <name><surname>Maussion</surname> <given-names>F.</given-names></name> <name><surname>Scherer</surname> <given-names>D.</given-names></name> <name><surname>Mayer</surname> <given-names>C.</given-names></name> <name><surname>Bush</surname> <given-names>A. B. G.</given-names></name></person-group> (<year>2013</year>). <article-title>High-resolution interactive modelling of the mountain glacier-atmosphere interface: an application over the Karakoram</article-title>. <source>Cryosphere</source> <volume>7</volume>, <fpage>779</fpage>&#x02013;<lpage>795</lpage>. <pub-id pub-id-type="doi">10.5194/tc-7-779-2013</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Collier</surname> <given-names>E.</given-names></name> <name><surname>Nicholson</surname> <given-names>L. I.</given-names></name> <name><surname>Brock</surname> <given-names>B. W.</given-names></name> <name><surname>Maussion</surname> <given-names>F.</given-names></name> <name><surname>Essery</surname> <given-names>R.</given-names></name> <name><surname>Bush</surname> <given-names>A. B. G.</given-names></name></person-group> (<year>2014</year>). <article-title>Representing moisture fluxes and phase changes in glacier debris cover using a reservoir approach</article-title>. <source>Cryosphere</source> <volume>8</volume>, <fpage>1429</fpage>&#x02013;<lpage>1444</lpage>. <pub-id pub-id-type="doi">10.5194/tc-8-1429-2014</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Copland</surname> <given-names>L.</given-names></name> <name><surname>Pope</surname> <given-names>S.</given-names></name> <name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Shroder</surname> <given-names>J. F.</given-names></name> <name><surname>Clendon</surname> <given-names>P.</given-names></name> <name><surname>Bush</surname> <given-names>A. B. G.</given-names></name> <etal/></person-group>. (<year>2009</year>). <article-title>Glacier velocities across the central Karakoram</article-title>. <source>Ann. Glaciol</source>. <volume>50</volume>, <fpage>41</fpage>&#x02013;<lpage>49</lpage>. <pub-id pub-id-type="doi">10.3189/172756409789624229</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Cuffey</surname> <given-names>K. M.</given-names></name> <name><surname>Paterson</surname> <given-names>W. S. B.</given-names></name></person-group> (<year>2010</year>). <source>The Physics of Glaciers</source>. <publisher-loc>Oxford</publisher-loc>: <publisher-name>Academic Press</publisher-name>.</citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dobreva</surname> <given-names>I. D.</given-names></name> <name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Bush</surname> <given-names>A. B. G.</given-names></name></person-group> (<year>2017</year>). <article-title>Climate-glacier dynamics and topographic forcing in the Karakoram Himalaya: concepts, issues and research directions</article-title>. <source>Water</source> <volume>9</volume>:<fpage>405</fpage>. <pub-id pub-id-type="doi">10.3390/w9060405</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dozier</surname> <given-names>J.</given-names></name> <name><surname>Bruno</surname> <given-names>J.</given-names></name> <name><surname>Downey</surname> <given-names>P.</given-names></name></person-group> (<year>1981</year>). <article-title>A faster solution to the horizon problem</article-title>. <source>Comput. Geosci</source>. <volume>7</volume>, <fpage>145</fpage>&#x02013;<lpage>151</lpage>. <pub-id pub-id-type="doi">10.1016/0098-3004(81)90026-1</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Egholm</surname> <given-names>D. L.</given-names></name> <name><surname>Knudsen</surname> <given-names>M. F.</given-names></name> <name><surname>Clark</surname> <given-names>C. D.</given-names></name> <name><surname>Lesemann</surname> <given-names>J. E.</given-names></name></person-group> (<year>2011</year>). <article-title>Modeling the flow of glaciers in steep terrains: the integrated second-order shallow ice approximation (isosia)</article-title>. <source>J. Geophys. Res. Earth Surface</source> <fpage>116</fpage>. <pub-id pub-id-type="doi">10.1029/2010JF001900</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Farinotti</surname> <given-names>D.</given-names></name> <name><surname>Brinkerhoff</surname> <given-names>D. J.</given-names></name> <name><surname>Clarke</surname> <given-names>G. K. C.</given-names></name> <name><surname>F&#x000FC;rst</surname> <given-names>J. J.</given-names></name> <name><surname>Frey</surname> <given-names>H.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>How accurate are estimates of glacier ice thickness? Results from ITMIX, the Ice Thickness Models Intercomparison experiment</article-title>. <source>Cryosphere</source> <volume>11</volume>, <fpage>949</fpage>&#x02013;<lpage>970</lpage>. <pub-id pub-id-type="doi">10.5194/tc-11-949-2017</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Farinotti</surname> <given-names>D.</given-names></name> <name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name> <name><surname>de Kok</surname> <given-names>R. J.</given-names></name> <name><surname>Quincey</surname> <given-names>D. J.</given-names></name> <name><surname>Dehecq</surname> <given-names>A.</given-names></name></person-group> (<year>2020</year>). <article-title>Manifestations and mechanisms of the Karakoram glacier anomaly</article-title>. <source>Nat. Geosci</source>. <volume>13</volume>, <fpage>8</fpage>&#x02013;<lpage>16</lpage>. <pub-id pub-id-type="doi">10.1038/s41561-019-0513-5</pub-id><pub-id pub-id-type="pmid">31915463</pub-id></citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Foster</surname> <given-names>L.</given-names></name> <name><surname>Brock</surname> <given-names>B.</given-names></name> <name><surname>Cutler</surname> <given-names>M.</given-names></name> <name><surname>Diotri</surname> <given-names>F.</given-names></name></person-group> (<year>2012</year>). <article-title>A physically based method for estimating supraglacial debris thickness from thermal band remote-sensing data</article-title>. <source>J. Glaciol</source>. <volume>58</volume>, <fpage>677</fpage>&#x02013;<lpage>691</lpage>. <pub-id pub-id-type="doi">10.3189/2012JoG11J194</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fountain</surname> <given-names>A. G.</given-names></name> <name><surname>Walder</surname> <given-names>J. S.</given-names></name></person-group> (<year>1998</year>). <article-title>Water flow through temperate glaciers</article-title>. <source>Rev. Geophys</source>. <volume>36</volume>, <fpage>299</fpage>&#x02013;<lpage>328</lpage>. <pub-id pub-id-type="doi">10.1029/97RG03579</pub-id><pub-id pub-id-type="pmid">15703743</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujita</surname> <given-names>K.</given-names></name> <name><surname>Sakai</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Modelling runoff from a himalayan debris-covered glacier</article-title>. <source>Hydrol. Earth Syst. Sci</source>. <volume>18</volume>, <fpage>2679</fpage>&#x02013;<lpage>2694</lpage>. <pub-id pub-id-type="doi">10.5194/hess-18-2679-2014</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fyffe</surname> <given-names>C. L.</given-names></name> <name><surname>Reid</surname> <given-names>T. D.</given-names></name> <name><surname>Brock</surname> <given-names>B. W.</given-names></name> <name><surname>Kirkbride</surname> <given-names>M. P.</given-names></name> <name><surname>Diolaiuti</surname> <given-names>G.</given-names></name> <name><surname>Smiraglia</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>A distributed energy-balance melt model of an alpine debris-covered glacier</article-title>. <source>J. Glaciol</source>. <volume>60</volume>, <fpage>587</fpage>&#x02013;<lpage>602</lpage>. <pub-id pub-id-type="doi">10.3189/2014JoG13J148</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fyffe</surname> <given-names>C. L.</given-names></name> <name><surname>Woodget</surname> <given-names>A. S.</given-names></name> <name><surname>Kirkbride</surname> <given-names>M. P.</given-names></name> <name><surname>Deline</surname> <given-names>P.</given-names></name> <name><surname>Westoby</surname> <given-names>M. J.</given-names></name> <name><surname>Brock</surname> <given-names>B. W.</given-names></name></person-group> (<year>2020</year>). <article-title>Processes at the margins of supraglacial debris cover: quantifying dirty ice ablation and debris redistribution</article-title>. <source>Earth Surface Process. Landforms</source> <volume>45</volume>, <fpage>2272</fpage>&#x02013;<lpage>2290</lpage>. <pub-id pub-id-type="doi">10.1002/esp.4879</pub-id></citation></ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Garg</surname> <given-names>P. K.</given-names></name> <name><surname>Shukla</surname> <given-names>A.</given-names></name> <name><surname>Jasrotia</surname> <given-names>A. S.</given-names></name></person-group> (<year>2017</year>). <article-title>Influence of topography on glacier changes in the central Himalaya, India</article-title>. <source>Glob. Planet. Change</source> <volume>155</volume>, <fpage>196</fpage>&#x02013;<lpage>212</lpage>. <pub-id pub-id-type="doi">10.1016/j.gloplacha.2017.07.007</pub-id></citation></ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gibson</surname> <given-names>M. J.</given-names></name> <name><surname>Glasser</surname> <given-names>N. F.</given-names></name> <name><surname>Quincey</surname> <given-names>D. J.</given-names></name> <name><surname>Mayer</surname> <given-names>C.</given-names></name> <name><surname>Rowan</surname> <given-names>A. V.</given-names></name> <name><surname>Irvine-Fynn</surname> <given-names>T. D.</given-names></name></person-group> (<year>2017</year>). <article-title>Temporal variations in supraglacial debris distribution on Baltoro glacier, Karakoram between 2001 and 2012</article-title>. <source>Geomorphology</source> <volume>295</volume>, <fpage>572</fpage>&#x02013;<lpage>585</lpage>. <pub-id pub-id-type="doi">10.1016/j.geomorph.2017.08.012</pub-id><pub-id pub-id-type="pmid">22951191</pub-id></citation></ref>
<ref id="B47">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Gueymard</surname> <given-names>C.</given-names></name></person-group> (<year>1995</year>). <source>SMARTS2, A Simple Model of the Atmospheric Radiative Transfer of Sunshine: Algorithms and Performance Assessment</source>. Florida Solar Energy Center, Cocoa, FL.</citation></ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gulley</surname> <given-names>J.</given-names></name> <name><surname>Benn</surname> <given-names>D.</given-names></name></person-group> (<year>2007</year>). <article-title>Structural control of englacial drainage systems in himalayan debris-covered glaciers</article-title>. <source>J. Glaciol</source>. <volume>53</volume>, <fpage>399</fpage>&#x02013;<lpage>412</lpage>. <pub-id pub-id-type="doi">10.3189/002214307783258378</pub-id></citation></ref>
<ref id="B49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hambrey</surname> <given-names>M. J.</given-names></name> <name><surname>Quincey</surname> <given-names>D. J.</given-names></name> <name><surname>Glasser</surname> <given-names>N. F.</given-names></name> <name><surname>Reynolds</surname> <given-names>J. M.</given-names></name> <name><surname>Richardson</surname> <given-names>S. J.</given-names></name> <name><surname>Clemmens</surname> <given-names>S.</given-names></name></person-group> (<year>2008</year>). <article-title>Sedimentological, geomorphological and dynamic context of debris-mantled glaciers, Mount Everest (Sagarmatha) region, Nepal</article-title>. <source>Quat. Sci. Rev</source>. <volume>27</volume>, <fpage>2361</fpage>&#x02013;<lpage>2389</lpage>. <pub-id pub-id-type="doi">10.1016/j.quascirev.2008.08.010</pub-id></citation></ref>
<ref id="B50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Herman</surname> <given-names>F.</given-names></name> <name><surname>Beaud</surname> <given-names>F.</given-names></name> <name><surname>Champagnac</surname> <given-names>J.-D.</given-names></name> <name><surname>Lemieux</surname> <given-names>J.-M.</given-names></name> <name><surname>Sternai</surname> <given-names>P.</given-names></name></person-group> (<year>2011</year>). <article-title>Glacial hydrology and erosion patterns: a mechanism for carving glacial valleys</article-title>. <source>Earth Planet. Sci. Lett</source>. <volume>310</volume>, <fpage>498</fpage>&#x02013;<lpage>508</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsl.2011.08.022</pub-id></citation></ref>
<ref id="B51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Herman</surname> <given-names>F.</given-names></name> <name><surname>Braun</surname> <given-names>J.</given-names></name></person-group> (<year>2008</year>). <article-title>Evolution of the glacial landscape of the southern Alps of New Zealand: insights from a glacial erosion model</article-title>. <source>J. Geophys. Res. Earth Surface</source> <fpage>113</fpage>. <pub-id pub-id-type="doi">10.1029/2007JF000807</pub-id></citation></ref>
<ref id="B52">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hewitt</surname> <given-names>K.</given-names></name></person-group> (<year>2005</year>). <article-title>The Karakoram anomaly? Glacier expansion and the &#x0201C;elevation effect,&#x0201D; Karakoram Himalaya</article-title>. <source>Mountain Res. Dev</source>. <volume>25</volume>, <fpage>332</fpage>&#x02013;<lpage>340</lpage>. <pub-id pub-id-type="doi">10.1659/0276-4741(2005)025[0332:TKAGEA]2.0.CO;2</pub-id></citation></ref>
<ref id="B53">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>L.</given-names></name> <name><surname>Li</surname> <given-names>Z.</given-names></name> <name><surname>Han</surname> <given-names>H.</given-names></name> <name><surname>Tian</surname> <given-names>B.</given-names></name> <name><surname>Zhou</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>Analysis of thickness changes and the associated driving factors on a debris-covered glacier in the tienshan mountain</article-title>. <source>Remote Sens. Environ</source>. <volume>206</volume>, <fpage>63</fpage>&#x02013;<lpage>71</lpage>. <pub-id pub-id-type="doi">10.1016/j.rse.2017.12.028</pub-id></citation></ref>
<ref id="B54">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Huo</surname> <given-names>D.</given-names></name> <name><surname>Bishop</surname> <given-names>M. P.</given-names></name> <name><surname>Young</surname> <given-names>B. W.</given-names></name> <name><surname>Chi</surname> <given-names>Z.</given-names></name> <name><surname>Haritashya</surname> <given-names>U. K.</given-names></name></person-group> (<year>2020</year>). <article-title>Numerical modeling issues for understanding complex debris-covered glaciers</article-title>, in <source>Reference Module in Earth Systems and Environmental Sciences</source> (<publisher-loc>Elsevier</publisher-loc>). <pub-id pub-id-type="doi">10.1016/B978-0-12-818234-5.00019-5</pub-id></citation></ref>
<ref id="B55">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huo</surname> <given-names>D.</given-names></name> <name><surname>Chi</surname> <given-names>Z.</given-names></name> <name><surname>Ma</surname> <given-names>A.</given-names></name></person-group> (<year>2021</year>). <article-title>Modeling surface processes on debris-covered glaciers: a review with reference to the high mountain asia</article-title>. <source>Water</source> Amsterdam, <volume>13</volume>:<fpage>101</fpage>. <pub-id pub-id-type="doi">10.3390/w13010101</pub-id></citation></ref>
<ref id="B56">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Immerzeel</surname> <given-names>W.</given-names></name> <name><surname>Kraaijenbrink</surname> <given-names>P.</given-names></name> <name><surname>Shea</surname> <given-names>J.</given-names></name> <name><surname>Shrestha</surname> <given-names>A.</given-names></name> <name><surname>Pellicciotti</surname> <given-names>F.</given-names></name> <name><surname>Bierkens</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>High-resolution monitoring of Himalayan glacier dynamics using unmanned aerial vehicles</article-title>. <source>Remote Sens. Environ</source>. <volume>150</volume>, <fpage>93</fpage>&#x02013;<lpage>103</lpage>. <pub-id pub-id-type="doi">10.1016/j.rse.2014.04.025</pub-id></citation></ref>
<ref id="B57">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Immerzeel</surname> <given-names>W.</given-names></name> <name><surname>Pellicciotti</surname> <given-names>F.</given-names></name> <name><surname>Bierkens</surname> <given-names>M.</given-names></name></person-group> (<year>2013</year>). <article-title>Rising river flows throughout the twenty-first century in two Himalayan glacierized watersheds</article-title>. <source>Nat. Geosci</source>. <volume>6</volume>, <fpage>742</fpage>&#x02013;<lpage>745</lpage>. <pub-id pub-id-type="doi">10.1038/ngeo1896</pub-id></citation></ref>
<ref id="B58">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name> <name><surname>Van Beek</surname> <given-names>L. P.</given-names></name> <name><surname>Bierkens</surname> <given-names>M. F.</given-names></name></person-group> (<year>2010</year>). <article-title>Climate change will affect the asian water towers</article-title>. <source>Science</source> <volume>328</volume>, <fpage>1382</fpage>&#x02013;<lpage>1385</lpage>. <pub-id pub-id-type="doi">10.1126/science.1183188</pub-id><pub-id pub-id-type="pmid">20538947</pub-id></citation></ref>
<ref id="B59">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Iqbal</surname> <given-names>M.</given-names></name></person-group> (<year>1983</year>). <source>An Introduction to Solar Radiation</source>. <publisher-loc>Toronto, ON</publisher-loc>: <publisher-name>Academic Press</publisher-name>.</citation></ref>
<ref id="B60">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Juen</surname> <given-names>M.</given-names></name> <name><surname>Mayer</surname> <given-names>C.</given-names></name> <name><surname>Lambrecht</surname> <given-names>A.</given-names></name> <name><surname>Han</surname> <given-names>H.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Impact of varying debris cover thickness on ablation: a case study for KOXKAR glacier in the Tien Shan</article-title>. <source>Cryosphere</source> <volume>8</volume>:<fpage>377</fpage>. <pub-id pub-id-type="doi">10.5194/tc-8-377-2014</pub-id></citation></ref>
<ref id="B61">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>K&#x000E4;&#x000E4;b</surname> <given-names>A.</given-names></name> <name><surname>Berthier</surname> <given-names>E.</given-names></name> <name><surname>Nuth</surname> <given-names>C.</given-names></name> <name><surname>Gardelle</surname> <given-names>J.</given-names></name> <name><surname>Arnaud</surname> <given-names>Y.</given-names></name></person-group> (<year>2012</year>). <article-title>Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas</article-title>. <source>Nature</source> <volume>488</volume>, <fpage>495</fpage>&#x02013;<lpage>498</lpage>. <pub-id pub-id-type="doi">10.1038/nature11324</pub-id><pub-id pub-id-type="pmid">22914167</pub-id></citation></ref>
<ref id="B62">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kayastha</surname> <given-names>R. B.</given-names></name> <name><surname>Takeuchi</surname> <given-names>Y.</given-names></name> <name><surname>Nakawo</surname> <given-names>M.</given-names></name> <name><surname>Ageta</surname> <given-names>Y.</given-names></name></person-group> (<year>2000</year>). <article-title>Practical prediction of ice melting beneath various thickness of debris cover on Khumbu glacier, Nepal, using a positive degree-day factor</article-title>. <source>IAHS Publ</source>. <volume>264</volume>, <fpage>71</fpage>&#x02013;<lpage>81</lpage>.</citation></ref>
<ref id="B63">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Khan</surname> <given-names>M. I.</given-names></name></person-group> (<year>1989</year>). <source>Ablation on Barpu glacier, Karakoram Himalaya, Pakistan a study of melt processes on a faceted, debris-covered ice surface</source> (<publisher-loc>Master thesis</publisher-loc>), Wilfrid Laurier University, Waterloo, ON.</citation></ref>
<ref id="B64">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liang</surname> <given-names>S.</given-names></name></person-group> (<year>2001</year>). <article-title>Narrowband to broadband conversions of land surface albedo I: algorithms</article-title>. <source>Remote Sens. Environ</source>. <volume>76</volume>, <fpage>213</fpage>&#x02013;<lpage>238</lpage>. <pub-id pub-id-type="doi">10.1016/S0034-4257(00)00205-4</pub-id></citation></ref>
<ref id="B65">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>L&#x000FC;thje</surname> <given-names>M.</given-names></name> <name><surname>Pedersen</surname> <given-names>L.</given-names></name> <name><surname>Reeh</surname> <given-names>N.</given-names></name> <name><surname>Greuell</surname> <given-names>W.</given-names></name></person-group> (<year>2006</year>). <article-title>Modelling the evolution of supraglacial lakes on the west Greenland ice-sheet margin</article-title>. <source>J. Glaciol</source>. <volume>52</volume>, <fpage>608</fpage>&#x02013;<lpage>618</lpage>. <pub-id pub-id-type="doi">10.3189/172756506781828386</pub-id></citation></ref>
<ref id="B66">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mattson</surname> <given-names>L.</given-names></name></person-group> (<year>1993</year>). <article-title>Ablation on debris covered glaciers: an example from the Rakhiot glacier, Punjab, Himalaya</article-title>. <source>Intern. Assoc. Hydrol. Sci</source>. <volume>218</volume>, <fpage>289</fpage>&#x02013;<lpage>296</lpage>.</citation></ref>
<ref id="B67">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>McNabb</surname> <given-names>R.</given-names></name> <name><surname>Hock</surname> <given-names>R.</given-names></name> <name><surname>O&#x00027;Neel</surname> <given-names>S.</given-names></name> <name><surname>Rasmussen</surname> <given-names>L. A.</given-names></name> <name><surname>Ahn</surname> <given-names>Y.</given-names></name> <name><surname>Braun</surname> <given-names>M.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Using surface velocities to calculate ice thickness and bed topography: a case study at Columbia glacier, Alaska, USA</article-title>. <source>J. Glaciol</source>. <volume>58</volume>, <fpage>1151</fpage>&#x02013;<lpage>1164</lpage>. <pub-id pub-id-type="doi">10.3189/2012JoG11J249</pub-id></citation></ref>
<ref id="B68">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mertes</surname> <given-names>J. R.</given-names></name> <name><surname>Thompson</surname> <given-names>S. S.</given-names></name> <name><surname>Booth</surname> <given-names>A. D.</given-names></name> <name><surname>Gulley</surname> <given-names>J. D.</given-names></name> <name><surname>Benn</surname> <given-names>D. I.</given-names></name></person-group> (<year>2017</year>). <article-title>A conceptual model of supra-glacial lake formation on debris-covered glaciers based on GPR facies analysis</article-title>. <source>Earth Surface Process. Landforms</source> <volume>42</volume>, <fpage>903</fpage>&#x02013;<lpage>914</lpage>. <pub-id pub-id-type="doi">10.1002/esp.4068</pub-id></citation></ref>
<ref id="B69">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mihalcea</surname> <given-names>C.</given-names></name> <name><surname>Mayer</surname> <given-names>C.</given-names></name> <name><surname>Diolaiuti</surname> <given-names>G.</given-names></name> <name><surname>D&#x00027;agata</surname> <given-names>C.</given-names></name> <name><surname>Smiraglia</surname> <given-names>C.</given-names></name> <name><surname>Lambrecht</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>Spatial distribution of debris thickness and melting from remote-sensing and meteorological data, at debris-covered Baltoro glacier, Karakoram, Pakistan</article-title>. <source>Ann. Glaciol</source>. <volume>48</volume>, <fpage>49</fpage>&#x02013;<lpage>57</lpage>. <pub-id pub-id-type="doi">10.3189/172756408784700680</pub-id></citation></ref>
<ref id="B70">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Miles</surname> <given-names>E. S.</given-names></name> <name><surname>Pellicciotti</surname> <given-names>F.</given-names></name> <name><surname>Willis</surname> <given-names>I. C.</given-names></name> <name><surname>Steiner</surname> <given-names>J. F.</given-names></name> <name><surname>Buri</surname> <given-names>P.</given-names></name> <name><surname>Arnold</surname> <given-names>N. S.</given-names></name></person-group> (<year>2016</year>). <article-title>Refined energy-balance modelling of a supraglacial pond, Langtang Khola, Nepal</article-title>. <source>Ann. Glaciol</source>. <volume>57</volume>, <fpage>29</fpage>&#x02013;<lpage>40</lpage>. <pub-id pub-id-type="doi">10.3189/2016AoG71A421</pub-id></citation></ref>
<ref id="B71">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Miles</surname> <given-names>E. S.</given-names></name> <name><surname>Steiner</surname> <given-names>J.</given-names></name> <name><surname>Willis</surname> <given-names>I.</given-names></name> <name><surname>Buri</surname> <given-names>P.</given-names></name> <name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name> <name><surname>Chesnokova</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Pond dynamics and supraglacial-englacial connectivity on debris-covered Lirung glacier, Nepal</article-title>. <source>Front. Earth Sci</source>. <volume>5</volume>:<fpage>69</fpage>. <pub-id pub-id-type="doi">10.3389/feart.2017.00069</pub-id></citation></ref>
<ref id="B72">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Miles</surname> <given-names>E. S.</given-names></name> <name><surname>Willis</surname> <given-names>I.</given-names></name> <name><surname>Buri</surname> <given-names>P.</given-names></name> <name><surname>Steiner</surname> <given-names>J. F.</given-names></name> <name><surname>Arnold</surname> <given-names>N. S.</given-names></name> <name><surname>Pellicciotti</surname> <given-names>F.</given-names></name></person-group> (<year>2018</year>). <article-title>Surface pond energy absorption across four Himalayan glaciers accounts for 1/8 of total catchment ice loss</article-title>. <source>Geophys. Res. Lett</source>. <volume>45</volume>, <fpage>10</fpage>&#x02013;<lpage>464</lpage>. <pub-id pub-id-type="doi">10.1029/2018GL079678</pub-id><pub-id pub-id-type="pmid">31031450</pub-id></citation></ref>
<ref id="B73">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moore</surname> <given-names>P. L.</given-names></name></person-group> (<year>2018</year>). <article-title>Stability of supraglacial debris</article-title>. <source>Earth Surface Process. Landforms</source> <volume>43</volume>, <fpage>285</fpage>&#x02013;<lpage>297</lpage>. <pub-id pub-id-type="doi">10.1002/esp.4244</pub-id></citation></ref>
<ref id="B74">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nakawo</surname> <given-names>M.</given-names></name> <name><surname>Young</surname> <given-names>G. J.</given-names></name></person-group> (<year>1981</year>). <article-title>Field experiments to determine the effect of a debris layer on ablation of glacier ice</article-title>. <source>Ann. Glaciol</source>. <volume>2</volume>, <fpage>85</fpage>&#x02013;<lpage>91</lpage>. <pub-id pub-id-type="doi">10.3189/172756481794352432</pub-id></citation></ref>
<ref id="B75">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nicholson</surname> <given-names>L.</given-names></name> <name><surname>Benn</surname> <given-names>D. I.</given-names></name></person-group> (<year>2006</year>). <article-title>Calculating ice melt beneath a debris layer using meteorological data</article-title>. <source>J. Glaciol</source>. <volume>52</volume>, <fpage>463</fpage>&#x02013;<lpage>470</lpage>. <pub-id pub-id-type="doi">10.3189/172756506781828584</pub-id></citation></ref>
<ref id="B76">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Oerlemans</surname> <given-names>J.</given-names></name> <name><surname>Klok</surname> <given-names>E. J.</given-names></name></person-group> (<year>2002</year>). <article-title>Energy balance of a glacier surface: analysis of automatic weather station data from the morteratschgletscher, switzerland</article-title>. <source>Arctic Antarct. Alpine Res</source>. <volume>34</volume>, <fpage>477</fpage>&#x02013;<lpage>485</lpage>. <pub-id pub-id-type="doi">10.1080/15230430.2002.12003519</pub-id></citation></ref>
<ref id="B77">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Olson</surname> <given-names>M.</given-names></name> <name><surname>Rupper</surname> <given-names>S.</given-names></name></person-group> (<year>2019</year>). <article-title>Impacts of topographic shading on direct solar radiation for valley glaciers in complex topography</article-title>. <source>Cryosphere</source> <volume>13</volume>, <fpage>29</fpage>&#x02013;<lpage>40</lpage>. <pub-id pub-id-type="doi">10.5194/tc-13-29-2019</pub-id></citation></ref>
<ref id="B78">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pelto</surname> <given-names>M. S.</given-names></name></person-group> (<year>2000</year>). <article-title>Mass balance of adjacent debris-covered and clean glacier ice in the north Cascades, Washington</article-title>. <source>IAHS Publ</source>. <fpage>35</fpage>&#x02013;<lpage>42</lpage>.</citation></ref>
<ref id="B79">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Perez</surname> <given-names>R.</given-names></name> <name><surname>Stewart</surname> <given-names>R.</given-names></name> <name><surname>Arbogast</surname> <given-names>C.</given-names></name> <name><surname>Seals</surname> <given-names>R.</given-names></name> <name><surname>Scott</surname> <given-names>J.</given-names></name></person-group> (<year>1986</year>). <article-title>An anisotropic hourly diffuse radiation model for sloping surfaces: description, performance validation, site dependency evaluation</article-title>. <source>Solar Energy</source> <volume>36</volume>, <fpage>481</fpage>&#x02013;<lpage>497</lpage>. <pub-id pub-id-type="doi">10.1016/0038-092X(86)90013-7</pub-id></citation></ref>
<ref id="B80">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pratap</surname> <given-names>B.</given-names></name> <name><surname>Dobhal</surname> <given-names>D.</given-names></name> <name><surname>Mehta</surname> <given-names>M.</given-names></name> <name><surname>Bhambri</surname> <given-names>R.</given-names></name></person-group> (<year>2015</year>). <article-title>Influence of debris cover and altitude on glacier surface melting: a case study on Dokriani glacier, central Himalaya, India</article-title>. <source>Ann. Glaciol</source>. <volume>56</volume>, <fpage>9</fpage>&#x02013;<lpage>16</lpage>. <pub-id pub-id-type="doi">10.3189/2015AoG70A971</pub-id></citation></ref>
<ref id="B81">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pritchard</surname> <given-names>M. S.</given-names></name> <name><surname>Bush</surname> <given-names>A. B. G.</given-names></name> <name><surname>Marshall</surname> <given-names>S. J.</given-names></name></person-group> (<year>2008</year>). <article-title>Neglecting ice-atmosphere interactions underestimates ice sheet melt in millennial-scale deglaciation simulations</article-title>. <source>Geophys. Res. Lett</source>. <fpage>35</fpage>. <pub-id pub-id-type="doi">10.1029/2007GL031738</pub-id></citation></ref>
<ref id="B82">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Proy</surname> <given-names>C.</given-names></name> <name><surname>Tanr&#x000E9;</surname> <given-names>D.</given-names></name> <name><surname>Deschamps</surname> <given-names>P. Y.</given-names></name></person-group> (<year>1989</year>). <article-title>Evaluation of topographic effects in remotely sensed data</article-title>. <source>Remote Sens. Environ</source>. <volume>30</volume>, <fpage>21</fpage>&#x02013;<lpage>32</lpage>. <pub-id pub-id-type="doi">10.1016/0034-4257(89)90044-8</pub-id></citation></ref>
<ref id="B83">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quincey</surname> <given-names>D.</given-names></name> <name><surname>Braun</surname> <given-names>M.</given-names></name> <name><surname>Glasser</surname> <given-names>N. F.</given-names></name> <name><surname>Bishop</surname> <given-names>M.</given-names></name> <name><surname>Hewitt</surname> <given-names>K.</given-names></name> <name><surname>Luckman</surname> <given-names>A.</given-names></name></person-group> (<year>2011</year>). <article-title>Karakoram glacier surge dynamics</article-title>. <source>Geophys. Res. Lett</source>. <fpage>38</fpage>. <pub-id pub-id-type="doi">10.1029/2011GL049004</pub-id></citation></ref>
<ref id="B84">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quincey</surname> <given-names>D.</given-names></name> <name><surname>Copland</surname> <given-names>L.</given-names></name> <name><surname>Mayer</surname> <given-names>C.</given-names></name> <name><surname>Bishop</surname> <given-names>M.</given-names></name> <name><surname>Luckman</surname> <given-names>A.</given-names></name> <name><surname>Belo</surname> <given-names>M.</given-names></name></person-group> (<year>2009</year>). <article-title>Ice velocity and climate variations for Baltoro glacier, Pakistan</article-title>. <source>J. Glaciol</source>. <volume>55</volume>, <fpage>1061</fpage>&#x02013;<lpage>1071</lpage>. <pub-id pub-id-type="doi">10.3189/002214309790794913</pub-id></citation></ref>
<ref id="B85">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quincey</surname> <given-names>D. J.</given-names></name> <name><surname>Glasser</surname> <given-names>N. F.</given-names></name></person-group> (<year>2009</year>). <article-title>Morphological and ice-dynamical changes on the Tasman glacier, New Zealand, 1990-2007</article-title>. <source>Glob. Planet. Change</source> <volume>68</volume>, <fpage>185</fpage>&#x02013;<lpage>197</lpage>. <pub-id pub-id-type="doi">10.1016/j.gloplacha.2009.05.003</pub-id></citation></ref>
<ref id="B86">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reid</surname> <given-names>T.</given-names></name> <name><surname>Brock</surname> <given-names>B.</given-names></name></person-group> (<year>2014</year>). <article-title>Assessing ice-cliff backwasting and its contribution to total ablation of debris-covered miage glacier, Mont Blanc Massif, Italy</article-title>. <source>J. Glaciol</source>. <volume>60</volume>, <fpage>3</fpage>&#x02013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.3189/2014JoG13J045</pub-id></citation></ref>
<ref id="B87">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reid</surname> <given-names>T.</given-names></name> <name><surname>Carenzo</surname> <given-names>M.</given-names></name> <name><surname>Pellicciotti</surname> <given-names>F.</given-names></name> <name><surname>Brock</surname> <given-names>B.</given-names></name></person-group> (<year>2012</year>). <article-title>Including debris cover effects in a distributed model of glacier ablation</article-title>. <source>J. Geophys. Res. Atmos</source>. <fpage>117</fpage>. <pub-id pub-id-type="doi">10.1029/2012JD017795</pub-id></citation></ref>
<ref id="B88">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reid</surname> <given-names>T. D.</given-names></name> <name><surname>Brock</surname> <given-names>B. W.</given-names></name></person-group> (<year>2010</year>). <article-title>An energy-balance model for debris-covered glaciers including heat conduction through the debris layer</article-title>. <source>J. Glaciol</source>. <volume>56</volume>, <fpage>903</fpage>&#x02013;<lpage>916</lpage>. <pub-id pub-id-type="doi">10.3189/002214310794457218</pub-id></citation></ref>
<ref id="B89">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reynolds</surname> <given-names>J. M.</given-names></name></person-group> (<year>2000</year>). <article-title>On the formation of supraglacial lakes on debris-covered glaciers</article-title>. <source>IAHS Publ</source>. <fpage>153</fpage>&#x02013;<lpage>164</lpage>.</citation></ref>
<ref id="B90">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reznichenko</surname> <given-names>N.</given-names></name> <name><surname>Davies</surname> <given-names>T.</given-names></name> <name><surname>Shulmeister</surname> <given-names>J.</given-names></name> <name><surname>McSaveney</surname> <given-names>M.</given-names></name></person-group> (<year>2010</year>). <article-title>Effects of debris on ice-surface melting rates: an experimental study</article-title>. <source>J. Glaciol</source>. <volume>56</volume>, <fpage>384</fpage>&#x02013;<lpage>394</lpage>. <pub-id pub-id-type="doi">10.3189/002214310792447725</pub-id></citation></ref>
<ref id="B91">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rounce</surname> <given-names>D.</given-names></name> <name><surname>Quincey</surname> <given-names>D.</given-names></name> <name><surname>McKinney</surname> <given-names>D.</given-names></name></person-group> (<year>2015</year>). <article-title>Debris-covered energy balance model for Imja-Lhotse Shar glacier in the Everest region of Nepal</article-title>. <source>Cryosphere Discuss</source> <volume>9</volume>, <fpage>3503</fpage>&#x02013;<lpage>3540</lpage>. <pub-id pub-id-type="doi">10.5194/tcd-9-3503-2015</pub-id></citation></ref>
<ref id="B92">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rounce</surname> <given-names>D. R.</given-names></name> <name><surname>King</surname> <given-names>O.</given-names></name> <name><surname>McCarthy</surname> <given-names>M.</given-names></name> <name><surname>Shean</surname> <given-names>D. E.</given-names></name> <name><surname>Salerno</surname> <given-names>F.</given-names></name></person-group> (<year>2018</year>). <article-title>Quantifying debris thickness of debris-covered glaciers in the Everest region of Nepal through inversion of a subdebris melt model</article-title>. <source>J. Geophys. Res. Earth Surface</source> <volume>123</volume>, <fpage>1094</fpage>&#x02013;<lpage>1115</lpage>. <pub-id pub-id-type="doi">10.1029/2017JF004395</pub-id></citation></ref>
<ref id="B93">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rowan</surname> <given-names>A.</given-names></name> <name><surname>Nicholson</surname> <given-names>L.</given-names></name> <name><surname>Quincey</surname> <given-names>D.</given-names></name> <name><surname>Gibson</surname> <given-names>M.</given-names></name> <name><surname>Irvine-Fynn</surname> <given-names>T.</given-names></name> <name><surname>Watson</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Seasonally stable temperature gradients through supraglacial debris in the Everest region of Nepal, Central Himalaya</article-title>. <source>J. Glaciol</source>. <volume>67</volume>, <fpage>170</fpage>&#x02013;<lpage>181</lpage>. <pub-id pub-id-type="doi">10.1017/jog.2020.100</pub-id></citation></ref>
<ref id="B94">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rowan</surname> <given-names>A. V.</given-names></name> <name><surname>Egholm</surname> <given-names>D. L.</given-names></name> <name><surname>Quincey</surname> <given-names>D. J.</given-names></name> <name><surname>Glasser</surname> <given-names>N. F.</given-names></name></person-group> (<year>2015</year>). <article-title>Modelling the feedbacks between mass balance, ice flow and debris transport to predict the response to climate change of debris-covered glaciers in the Himalaya</article-title>. <source>Earth Planet. Sci. Lett</source>. <volume>430</volume>, <fpage>427</fpage>&#x02013;<lpage>438</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsl.2015.09.004</pub-id></citation></ref>
<ref id="B95">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sakai</surname> <given-names>A.</given-names></name> <name><surname>Nakawo</surname> <given-names>M.</given-names></name> <name><surname>Fujita</surname> <given-names>K.</given-names></name></person-group> (<year>2002</year>). <article-title>Distribution characteristics and energy balance of ice cliffs on debris-covered glaciers, Nepal Himalaya</article-title>. <source>Arctic Antarct. Alpine Res</source>. <volume>34</volume>, <fpage>12</fpage>&#x02013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.1080/15230430.2002.12003463</pub-id></citation></ref>
<ref id="B96">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sakai</surname> <given-names>A.</given-names></name> <name><surname>Takeuchi</surname> <given-names>N.</given-names></name> <name><surname>Fujita</surname> <given-names>K.</given-names></name> <name><surname>Nakawo</surname> <given-names>M.</given-names></name></person-group> (<year>2000</year>). <article-title>Role of supraglacial ponds in the ablation process of a debris-covered glacier in the nepal himalayas</article-title>. <source>IAHS Publ</source>. <fpage>119</fpage>&#x02013;<lpage>132</lpage>.</citation></ref>
<ref id="B97">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Salerno</surname> <given-names>F.</given-names></name> <name><surname>Thakuri</surname> <given-names>S.</given-names></name> <name><surname>Tartari</surname> <given-names>G.</given-names></name> <name><surname>Nuimura</surname> <given-names>T.</given-names></name> <name><surname>Sunako</surname> <given-names>S.</given-names></name> <name><surname>Sakai</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Debris-covered glacier anomaly? Morphological factors controlling changes in the mass balance, surface area, terminus position, and snow line altitude of Himalayan glaciers</article-title>. <source>Earth Planet. Sci. Lett</source>. <volume>471</volume>, <fpage>19</fpage>&#x02013;<lpage>31</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsl.2017.04.039</pub-id></citation></ref>
<ref id="B98">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Scherler</surname> <given-names>D.</given-names></name> <name><surname>Bookhagen</surname> <given-names>B.</given-names></name> <name><surname>Strecker</surname> <given-names>M. R.</given-names></name></person-group> (<year>2011a</year>). <article-title>Hillslope-glacier coupling: the interplay of topography and glacial dynamics in high Asia</article-title>. <source>J. Geophys. Res. Earth Surface</source> <fpage>116</fpage>. <pub-id pub-id-type="doi">10.1029/2010JF001751</pub-id></citation></ref>
<ref id="B99">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Scherler</surname> <given-names>D.</given-names></name> <name><surname>Bookhagen</surname> <given-names>B.</given-names></name> <name><surname>Strecker</surname> <given-names>M. R.</given-names></name></person-group> (<year>2011b</year>). <article-title>Spatially variable response of Himalayan glaciers to climate change affected by debris cover</article-title>. <source>Nat. Geosci</source>. <volume>4</volume>, <fpage>156</fpage>&#x02013;<lpage>159</lpage>. <pub-id pub-id-type="doi">10.1038/ngeo1068</pub-id></citation></ref>
<ref id="B100">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Seong</surname> <given-names>Y. B.</given-names></name> <name><surname>Owen</surname> <given-names>L. A.</given-names></name> <name><surname>Yi</surname> <given-names>C.</given-names></name> <name><surname>Finkel</surname> <given-names>R. C.</given-names></name></person-group> (<year>2009</year>). <article-title>Quaternary glaciation of Muztag Ata and Kongur Shan: evidence for glacier response to rapid climate changes throughout the late glacial and holocene in westernmost tibet</article-title>. <source>Geol. Soc. Am. Bull</source>. <volume>121</volume>, <fpage>348</fpage>&#x02013;<lpage>365</lpage>. <pub-id pub-id-type="doi">10.1130/B26339.1</pub-id></citation></ref>
<ref id="B101">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sicart</surname> <given-names>J.-E.</given-names></name> <name><surname>Hock</surname> <given-names>R.</given-names></name> <name><surname>Six</surname> <given-names>D.</given-names></name></person-group> (<year>2008</year>). <article-title>Glacier melt, air temperature, and energy balance in different climates: the Bolivian Tropics, the French Alps, and northern Sweden</article-title>. <source>J. Geophys. Res. Atmos</source>. <fpage>113</fpage>. <pub-id pub-id-type="doi">10.1029/2008JD010406</pub-id></citation></ref>
<ref id="B102">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Steer</surname> <given-names>P.</given-names></name> <name><surname>Huismans</surname> <given-names>R. S.</given-names></name> <name><surname>Valla</surname> <given-names>P. G.</given-names></name> <name><surname>Gac</surname> <given-names>S.</given-names></name> <name><surname>Herman</surname> <given-names>F.</given-names></name></person-group> (<year>2012</year>). <article-title>Bimodal plio-quaternary glacial erosion of fjords and low-relief surfaces in scandinavia</article-title>. <source>Nat. Geosci</source>. <volume>5</volume>, <fpage>635</fpage>&#x02013;<lpage>639</lpage>. <pub-id pub-id-type="doi">10.1038/ngeo1549</pub-id></citation></ref>
<ref id="B103">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Steiner</surname> <given-names>J. F.</given-names></name> <name><surname>Pellicciotti</surname> <given-names>F.</given-names></name> <name><surname>Buri</surname> <given-names>P.</given-names></name> <name><surname>Miles</surname> <given-names>E. S.</given-names></name> <name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name> <name><surname>Reid</surname> <given-names>T. D.</given-names></name></person-group> (<year>2015</year>). <article-title>Modelling ice-cliff backwasting on a debris-covered glacier in the Nepalese Himalaya</article-title>. <source>J. Glaciol</source>. <volume>61</volume>, <fpage>889</fpage>&#x02013;<lpage>907</lpage>. <pub-id pub-id-type="doi">10.3189/2015JoG14J194</pub-id></citation></ref>
<ref id="B104">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sutherland</surname> <given-names>J. L.</given-names></name> <name><surname>Carrivick</surname> <given-names>J. L.</given-names></name> <name><surname>Gandy</surname> <given-names>N.</given-names></name> <name><surname>Shulmeister</surname> <given-names>J.</given-names></name> <name><surname>Quincey</surname> <given-names>D. J.</given-names></name> <name><surname>Cornford</surname> <given-names>S. L.</given-names></name></person-group> (<year>2020</year>). <article-title>Proglacial lakes control glacier geometry and behavior during recession</article-title>. <source>Geophys. Res. Lett</source>. <volume>47</volume>:<fpage>e2020G</fpage>L088865. <pub-id pub-id-type="doi">10.1029/2020GL088865</pub-id></citation></ref>
<ref id="B105">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Takeuchi</surname> <given-names>N.</given-names></name> <name><surname>Li</surname> <given-names>Z.</given-names></name></person-group> (<year>2008</year>). <article-title>Characteristics of surface dust on &#x000DC;r&#x000FC;mqi glacier no. 1 in the tien Shan mountains, China</article-title>. <source>Arctic Antarct. Alpine Res</source>. <volume>40</volume>, <fpage>744</fpage>&#x02013;<lpage>750</lpage>. <pub-id pub-id-type="doi">10.1657/1523-0430(07-094)[TAKEUCHI]2.0.CO;2</pub-id></citation></ref>
<ref id="B106">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thompson</surname> <given-names>S.</given-names></name> <name><surname>Benn</surname> <given-names>D. I.</given-names></name> <name><surname>Mertes</surname> <given-names>J.</given-names></name> <name><surname>Luckman</surname> <given-names>A.</given-names></name></person-group> (<year>2016</year>). <article-title>Stagnation and mass loss on a himalayan debris-covered glacier: processes, patterns and rates</article-title>. <source>J. Glaciol</source>. <volume>62</volume>, <fpage>467</fpage>&#x02013;<lpage>485</lpage>. <pub-id pub-id-type="doi">10.1017/jog.2016.37</pub-id></citation></ref>
<ref id="B107">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>van Woerkom</surname> <given-names>T.</given-names></name> <name><surname>Steiner</surname> <given-names>J. F.</given-names></name> <name><surname>Kraaijenbrink</surname> <given-names>P. D.</given-names></name> <name><surname>Miles</surname> <given-names>E.</given-names></name> <name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name></person-group> (<year>2019</year>). <article-title>Sediment supply from lateral moraines to a debris-covered glacier in the Himalaya</article-title>. <source>Earth Surface Dyn</source>. <volume>7</volume>, <fpage>411</fpage>&#x02013;<lpage>427</lpage>. <pub-id pub-id-type="doi">10.5194/esurf-7-411-2019</pub-id></citation></ref>
<ref id="B108">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vincent</surname> <given-names>C.</given-names></name> <name><surname>Wagnon</surname> <given-names>P. W.</given-names></name> <name><surname>Shea</surname> <given-names>J. M.</given-names></name> <name><surname>Immerzeel</surname> <given-names>W. W.</given-names></name> <name><surname>Kraaijenbrink</surname> <given-names>P.</given-names></name> <name><surname>Shrestha</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Reduced melt on debris-covered glaciers: investigations from Changri Nup Glacier, Nepal</article-title>. <source>Cryosphere</source> <volume>10</volume>, <fpage>1845</fpage>&#x02013;<lpage>1858</lpage>. <pub-id pub-id-type="doi">10.5194/tc-10-1845-2016</pub-id></citation></ref>
<ref id="B109">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wessels</surname> <given-names>R. L.</given-names></name> <name><surname>Kargel</surname> <given-names>J. S.</given-names></name> <name><surname>Kieffer</surname> <given-names>H. H.</given-names></name></person-group> (<year>2002</year>). <article-title>Aster measurement of supraglacial lakes in the Mount Everest region of the Himalaya</article-title>. <source>Ann. Glaciol</source>. <volume>34</volume>, <fpage>399</fpage>&#x02013;<lpage>408</lpage>. <pub-id pub-id-type="doi">10.3189/172756402781817545</pub-id></citation></ref>
<ref id="B110">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Wilson</surname> <given-names>J. P.</given-names></name> <name><surname>Bishop</surname> <given-names>M. P.</given-names></name></person-group> (<year>2013</year>). <article-title>Geomorphometry</article-title>, in <source>Remote Sensing and GIScience in Geomorphology</source>, ed <person-group person-group-type="editor"><name><surname>Bishop</surname> <given-names>J. S.</given-names></name></person-group> (<publisher-loc>Waltham, MA</publisher-loc>: <publisher-name>Academic Press</publisher-name>), <fpage>162</fpage>&#x02013;<lpage>186</lpage>. <pub-id pub-id-type="doi">10.1016/B978-0-12-374739-6.00049-X</pub-id></citation></ref>
<ref id="B111">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wirbel</surname> <given-names>A.</given-names></name> <name><surname>Jarosch</surname> <given-names>A. H.</given-names></name> <name><surname>Nicholson</surname> <given-names>L.</given-names></name></person-group> (<year>2018</year>). <article-title>Modelling debris transport within glaciers by advection in a full-stokes ice flow model</article-title>. <source>Cryosphere</source> <volume>12</volume>, <fpage>189</fpage>&#x02013;<lpage>204</lpage>. <pub-id pub-id-type="doi">10.5194/tc-12-189-2018</pub-id></citation></ref>
<ref id="B112">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yue</surname> <given-names>X.</given-names></name> <name><surname>Zhao</surname> <given-names>J.</given-names></name> <name><surname>Li</surname> <given-names>Z.</given-names></name> <name><surname>Zhang</surname> <given-names>M.</given-names></name> <name><surname>Fan</surname> <given-names>J.</given-names></name> <name><surname>Wang</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Spatial and temporal variations of the surface albedo and other factors influencing Urumqi glacier no. 1 in Tien Shan, China</article-title>. <source>J. Glaciol</source>. <volume>63</volume>, <fpage>899</fpage>&#x02013;<lpage>911</lpage>. <pub-id pub-id-type="doi">10.1017/jog.2017.57</pub-id></citation></ref>
<ref id="B113">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Fujita</surname> <given-names>K.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name> <name><surname>Liu</surname> <given-names>Q.</given-names></name> <name><surname>Nuimura</surname> <given-names>T.</given-names></name></person-group> (<year>2011</year>). <article-title>Distribution of debris thickness and its effect on ice melt at Hailuogou glacier, southeastern Tibetan plateau, using <italic>in situ</italic> surveys and aster imagery</article-title>. <source>J. Glaciol</source>. <volume>57</volume>, <fpage>1147</fpage>&#x02013;<lpage>1157</lpage>. <pub-id pub-id-type="doi">10.3189/002214311798843331</pub-id></citation></ref>
<ref id="B114">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Hirabayashi</surname> <given-names>Y.</given-names></name> <name><surname>Fujita</surname> <given-names>K.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name> <name><surname>Liu</surname> <given-names>Q.</given-names></name></person-group> (<year>2016</year>). <article-title>Heterogeneity in supraglacial debris thickness and its role in glacier mass changes of the mount Gongga</article-title>. <source>Sci. China Earth Sci</source>. <volume>59</volume>, <fpage>170</fpage>&#x02013;<lpage>184</lpage>. <pub-id pub-id-type="doi">10.1007/s11430-015-5118-2</pub-id></citation></ref>
<ref id="B115">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>X.</given-names></name> <name><surname>Wen</surname> <given-names>J.</given-names></name> <name><surname>Liu</surname> <given-names>Q.</given-names></name> <name><surname>Yan</surname> <given-names>G.</given-names></name></person-group> (<year>2015</year>). <article-title>Improved topographic normalization for landsat tm images by introducing the MODIS surface BRDF</article-title>. <source>Remote Sens</source>. <volume>7</volume>, <fpage>6558</fpage>&#x02013;<lpage>6575</lpage>. <pub-id pub-id-type="doi">10.3390/rs70606558</pub-id></citation></ref>
</ref-list> 
</back>
</article>