Heterogeneity in supraglacial debris thickness and its role in glacier mass changes of the Mount Gongga

In the Tibetan Plateau, many glaciers have extensive covers of supraglacial debris in their ablation zones, which affects glacier response to climate change by altering ice melting and spatial patterns of mass loss. Insufficient debris thickness data make it difficult to analyze regional debris-cover effects. Maritime glaciers of the Mount Gongga have been characterized by a substantial reduction in glacier area and ice mass in recent decades. The thermal property of the debris layer estimated from remotely sensed data reveals that debris-covered glaciers are dominant in this region, on which the proportion of debris cover to total glacier area varies from 1.74% to 53.0%. Using a physically-based debris-cover effect assessment model, we found that although the presence of supraglacial debris has a significant insulating effect on heavily debris-covered glaciers, it accelerates ice melting on ~10.2% of total ablation zone and produces rapid wastage of ~25% of the debris-covered glaciers, leading to the similar mass losses between the debris-covered and debris-free glaciers. Widespread debris cover also facilitates the development of active terminus regions. Regional differences in debris-cover effects are apparent, highlighting the importance of debris cover for understanding glacier mass changes in the Tibetan Plateau and other mountain ranges around the world.     

Evolution of debris cover on glaciers of the Eastern Alps,Austria, between 1996 and 2015

Debris cover on glaciers is an important component of glacial systems as it influences climate–glacier dynamics and thus the lifespan of glaciers. Increasing air temperatures, permafrost thaw and rock faces freshly exposed by glacier downwasting in accumulation zones result in increased rockfall activity and debris input. In the ablation zone, negative mass balances result in an enhanced melt-out of englacial debris. Glacier debris cover thus represents a clear signal of climate warming in mountain areas. To assess the temporal development of debris on glaciers of the Eastern Alps, Austria, we mapped debris cover on 255 glaciers using Landsat data at three time steps. We applied a ratio-based threshold classification technique and analysed glacier catchment characteristics to understand debris sources better. Across the Austrian Alps, debris cover increased by more than 10% between 1996 and 2015 while glaciers retreated in response to climate warming. Debris cover distribution shows significant regional variability, with some mountain ranges being characterised by mean debris cover on glaciers of up to 75%. We also observed a general rise of the mean elevation of debris cover on glaciers in Austria. The debris cover distribution and dynamics are highly variable due to topographic, lithological and structural settings that determine the amount of debris delivered to and stored in the glacier system. Despite strong variation in debris cover, all glaciers investigated melted at increasing rates. We conclude that the retarding effects of debris cover on the mass balance and melt rate of Austrian glaciers is strongly subdued compared with other mountain areas. The study indicates that, if this trend continues, many glaciers in Austria may become fully debris covered. However, since debris cover seems to have little impact on melt rates, this would not lead to prolonged existence of debris-covered ice compared with clean ice glaciers.     

The formation of supraglacial debris covers by primary dispersal from transverse englacial debris bands

Glaciological controls on debris cover formation are investigated from the perspective of primary dispersal of supraglacial debris across a melting ice surface. This involves the migration of angled debris septa outcrops across a melting, thinning glacier ablation zone. Three measures of a glacier's ability to evacuate supraglacial debris are outlined: (1) a concentration factor describing the focusing of englacial debris into specific supraglacial mass loads; (2) the rate of migration of a septum outcrop relative to the local ice surface; and (3) a downstream velocity differential between a slower septum outcrop and the faster ice surface velocity. Measures (1) and (2) are inversely related, while measure (3) increases down‐glacier to explain why slow‐moving, thinning ice rapidly becomes debris covered. Data from Glacier d'Estelette (Italian Alps) are used to illustrate these processes, and to explore the potential for debris cover formation and growth in different glaciological environments. The transition from a ‘clean’, transport‐dominated to a debris‐covered ablation‐dominated glacier is explained by the melting out of more closely‐spaced debris septa, in combination with the geometric interactions of angled septa and ice surface in a field of reducing flow and increasing ablation. The growth and shrinkage of debris covers are most sensitive to glaciological changes at glaciers with gently‐dipping debris‐bearing foliation, but less sensitive at high‐compression glaciers whose termini are constrained by moraine dams and other forms of obstruction. These findings show that a variety of debris‐covered glacier types will show a spectrum of response characteristics to negative mass balance. Copyright © 2013 John Wiley & Sons, Ltd.     

An enhanced temperature index model for debris-covered glaciers accounting for thickness effect

Debris-covered glaciers are increasingly studied because it is assumed that debris cover extent and thickness could increase in a warming climate, with more regular rockfalls from the surrounding slopes and more englacial melt-out material. Debris energy-balance models have been developed to account for the melt rate enhancement/reduction due to a thin/thick debris layer, respectively. However, such models require a large amount of input data that are not often available, especially in remote mountain areas such as the Himalaya, and can be difficult to extrapolate. Due to their lower data requirements, empirical models have been used extensively in clean glacier melt modelling. For debris-covered glaciers, however, they generally simplify the debris effect by using a single melt-reduction factor which does not account for the influence of varying debris thickness on melt and prescribe a constant reduction for the entire melt across a glacier.In this paper, we present a new temperature-index model that accounts for debris thickness in the computation of melt rates at the debris-ice interface. The model empirical parameters are optimized at the point scale for varying debris thicknesses against melt rates simulated by a physically-based debris energy balance model. The latter is validated against ablation stake readings and surface temperature measurements. Each parameter is then related to a plausible set of debris thickness values to provide a general and transferable parameterization. We develop the model on Miage Glacier, Italy, and then test its transferability on Haut Glacier d’Arolla, Switzerland.The performance of the new debris temperature-index (DETI) model in simulating the glacier melt rate at the point scale is comparable to the one of the physically based approach, and the definition of model parameters as a function of debris thickness allows the simulation of the nonlinear relationship of melt rate to debris thickness, summarised by the Østrem curve. Its large number of parameters might be a limitation, but we show that the model is transferable in time and space to a second glacier with little loss of performance. We thus suggest that the new DETI model can be included in continuous mass balance models of debris-covered glaciers, because of its limited data requirements. As such, we expect its application to lead to an improvement in simulations of the debris-covered glacier response to climate in comparison with models that simply recalibrate empirical parameters to prescribe a constant across glacier reduction in melt.     

Morphology and evolution of supraglacial hummocks on debris-covered Himalayan glaciers

Thick supraglacial debris layers often have an undulating, hummocky topography that influences the lateral transport of debris and meltwater and provides basins for supraglacial ponds. The role of ablation and other processes associated with supraglacial debris in giving rise to this hummocky topography is poorly understood. Characterizing hummocky topography is a first step towards understanding the feedbacks driving the evolution of debris-covered glacier surfaces and their potential impacts on mass balance, hydrology and glacier dynamics. Here we undertake a geomorphological assessment of the hummocky topography on five debris-covered glaciers in the Everest region of the central Himalaya. We characterize supraglacial hummocks through statistical analyses of their vertical relief and horizontal geometry. Our results establish supraglacial hummocks as a distinct landform. We find that a typical hummock has an elongation ratio of 1.1:1 in the direction of ice flow, length of 214 ± 109 m and width of 192 ± 88 m. Hummocky topography has a greater amplitude across-glacier (15.4 ± 10.9 m) compared to along the glacier flow line (12.6 ± 8.3 m). Consequently, hummock slopes are steeper in the across-glacier direction (8.7 ± 4.3°) than in the direction of ice flow (5.6 ± 4.0°). Longer, wider and higher-amplitude hummocks are found on larger glaciers. We postulate that directional anisotropy in the hummock topography arises because, while the pattern of differential ablation driving topography evolution is moderated by processes including the gravitational redistribution of debris across the glacier surface, it also inherits an orientation preference from the distribution of englacial debris in the underlying ice. Our morphometric data inform future efforts to model these interactions, which should account for additional factors such as the genesis of supraglacial ponds and ice cliffs and their impact on differential ablation.     

Seasonal sediment fluxes forcing supraglacial melting on the Wright Lower Glacier,McMurdo Dry Valleys,Antarctica

Cold‐based glaciers exist in low temperature and low humidity environments in which shortwave radiation is the largest source of energy to the glacier surface and the energy budget is very sensitive to the surface albedo. Consequently, the presence of relatively low volumes of debris on glacier surfaces has a significant impact on the timing, magnitude and rate of ablation at the surface. The aim of this study is to understand how the presence of sediment on the glacier surface at the start of the melt season can affect meltwater generation and delivery on a cold‐based glacier. A combination of field measurements, energy balance modelling and chemical mixing modelling were used on the Wright Lower Glacier, McMurdo Dry Valleys, Antarctica, between October 2005 and January 2006 to address this aim. In this system, sediment was transported onto the glacier surface during the winter months (March–October) by foehn winds, which reduced surface albedo at the start of the summer melt season. The areas of the glacier on which sediment accumulated began to melt earlier than other parts of the glacier and experienced a longer melt season. Over the study period, the total ablation on the dirty surfaces was nine times greater than for clean ice. Ablation on the dirty surfaces is dominated by melting, whereas sublimation dominates the clean ice. As the sediment was unevenly distributed over the glacier surface, the variation in melt amount and timing drove the development of a cryoconite hole system. Copyright © 2012 John Wiley & Sons, Ltd.     

Geomorphological evolution of a debris-covered glacier surface

There exists a need to advance our understanding of debris-covered glacier surfaces over relatively short timescales due to rapid, climatically induced areal expansion of debris cover at the global scale, and the impact debris has on mass balance. We applied unpiloted aerial vehicle structure-from-motion (UAV-SfM) and digital elevation model (DEM) differencing with debris thickness and debris stability modelling to unravel the evolution of a 0.15 km² region of the debris-covered Miage Glacier, Italy, between June 2015 and July 2018. DEM differencing revealed widespread surface lowering (mean 4.1 ± 1.0 m a^-1; maximum 13.3 m a^-1). We combined elevation change data with local meteorological data and a sub-debris melt model, and used these relationships to produce high resolution, spatially distributed maps of debris thickness. These maps were differenced to explore patterns and mechanisms of debris redistribution. Median debris thicknesses ranged from 0.12 to 0.17 m and were spatially variable. We observed localized debris thinning across ice cliff faces, except those which were decaying, where debris thickened. We observed pervasive debris thinning across larger, backwasting slopes, including those bordered by supraglacial streams, as well as ingestion of debris by a newly exposed englacial conduit. Debris stability mapping showed that 18.2–26.4% of the survey area was theoretically subject to debris remobilization. By linking changes in stability to changes in debris thickness, we observed that slopes that remain stable, stabilize, or remain unstable between periods almost exclusively show net debris thickening (mean 0.07 m a^-1) whilst those which become newly unstable exhibit both debris thinning and thickening. We observe a systematic downslope increase in the rate at which debris cover thickens which can be described as a function of the topographic position index and slope gradient. Our data provide quantifiable insights into mechanisms of debris remobilization on glacier surfaces over sub-decadal timescales, and open avenues for future research to explore glacier-scale spatiotemporal patterns of debris remobilization. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Properties of natural supraglacial debris in relation to modelling sub‐debris ice ablation

As debris‐covered glaciers become a more prominent feature of a shrinking mountain cryosphere, there is increasing need to successfully model the surface energy and mass balance of debris‐covered glaciers, yet measurements of the processes operating in natural supraglacial debris covers are sparse. We report measurements of vertical temperature profiles in debris on the Ngozumpa glacier in Nepal, that show: (i) conductive processes dominate during the ablation season in matrix‐supported diamict; (ii) ventilation may be possible in coarse surface layers; (iii) phase changes associated with seasonal change have a marked effect on the effective thermal diffusivity of the debris. Effective thermal conductivity determined from vertical temperature profiles in the debris is generally ~30% higher in summer than in winter, but values depend on the volume and phase of water in the debris. Surface albedo can vary widely over small spatial scales, as does the debris thickness. Measurements indicate that debris thickness is best represented as a probability density function with the peak debris thickness increasing down‐glacier. The findings from Ngozumpa glacier indicate that the probability distribution of debris thickness changes from positively skewed in the upper glacier towards a more normal distribution nearer the terminus. Although many of these effects remain to be quantified, our observations highlight aspects of spatial and temporal variability in supraglacial debris that may require consideration in annual or multi‐annual distributed modelling of debris‐covered glacier surface energy and mass balance. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling the hydrological response of debris‐free and debris‐covered glaciers to present climatic conditions in the semiarid Andes of central Chile

We apply the process‐based, distributed TOPKAPI‐ETH glacio‐hydrological model to a glacierized catchment (19% glacierized) in the semiarid Andes of central Chile. The semiarid Andes provides vital freshwater resources to valleys in Chile and Argentina, but only few glacio‐hydrological modelling studies have been conducted, and its dominant hydrological processes remain poorly understood. The catchment contains two debris‐free glaciers reaching down to 3900 m asl (Bello and Yeso glaciers) and one debris‐covered avalanche‐fed glacier reaching to 3200 m asl (Piramide Glacier). Our main objective is to compare the mass balance and runoff contributions of both glacier types under current climatic conditions. We use a unique dataset of field measurements collected over two ablation seasons combined with the distributed TOPKAPI‐ETH model that includes physically oriented parameterizations of snow and ice ablation, gravitational distribution of snow, snow albedo evolution and the ablation of debris‐covered ice. Model outputs indicate that while the mass balance of Bello and Yeso glaciers is mostly explained by temperature gradients, the Piramide Glacier mass balance is governed by debris thickness and avalanches and has a clear non‐linear profile with elevation as a result. Despite the thermal insulation effect of the debris cover, the mass balance and contribution to runoff from debris‐free and debris‐covered glaciers are similar in magnitude, mainly because of elevation differences. However, runoff contributions are distinct in time and seasonality with ice melt starting approximately four weeks earlier from the debris‐covered glacier, what is of relevance for water resources management. At the catchment scale, snowmelt is the dominant contributor to runoff during both years. However, during the driest year of our simulations, ice melt contributes 42 ± 8% and 67 ± 6% of the annual and summer runoff, respectively. Sensitivity analyses show that runoff is most sensitive to temperature and precipitation gradients, melt factors and debris cover thickness. Copyright © 2016 John Wiley & Sons, Ltd.     

Ice flow‐unit influence on glacier structure,debris entrainment and transport

The links between structural glaciology, glacial debris entrainment and transport have been established in a number of different glacier settings. Here we document the structural evolution of a temperate Alpine valley glacier (Vadrec del Forno, Switzerland) and demonstrate that individual flow units within the glacier have very different structural and debris characteristics. The glacier consists of a broad accumulation area with multiple basins feeding a relatively narrow tongue and is formed from six distinct flow units. Each flow unit has its own characteristic structural assemblage. Flow units that narrow rapidly down‐glacier are dominated by primary stratification that has evolved into longitudinal foliation. In contrast, wider flow units preferentially develop an axial planar foliation. Glacier structure plays a limited role in the entrainment of debris, which is more strongly influenced by ice‐marginal rockfall and avalanche inputs onto the glacier surface. However, once entrained, glacier structure controls the reorientation and redistribution of debris within the ice mass. By taking a whole‐glacier approach to describing glacier structure and debris transport, we conclude that individual flow units are unique with regard to structure and debris transfer. Copyright © 2013 John Wiley & Sons, Ltd.     

Implications of decadal to century scale glacio‐hydrological change for water resources of the Hood River basin,OR, USA

In glacier‐fed rivers, melting of glacier ice sustains streamflow during the driest times of the year, especially during drought years. Anthropogenic and ecologic systems that rely on this glacial buffering of low flows are vulnerable to glacier recession as temperatures rise. We demonstrate the evolution of glacier melt contribution in watershed hydrology over the course of a 184‐year period from 1916 to 2099 through the application of a coupled hydrological and glacier dynamics model to the Hood River basin in Northwest Oregon, USA. We performed continuous simulations of glaciological processes (mass accumulation and ablation, lateral flow of ice and heat conduction through supra‐glacial debris), which are directly linked with seasonal snow dynamics as well as other key hydrologic processes (e.g. evapotranspiration and subsurface flow). Our simulations show that historically, the contribution of glacier melt to basin water supply was up to 79% at upland water management locations. We also show that supraglacial debris cover on the Hood River glaciers modulates the rate of glacier recession and progression of dry season flow at upland stream locations with debris‐covered glaciers. Our model results indicate that dry season (July to September) discharge sourced from glacier melt started to decline early in the 21st century following glacier recession that started early in the 20th century. Changes in climate over the course of the current century will lead to 14–63% (18–78%) reductions in dry season discharge across the basin for IPCC emission pathway RCP4.5 (RCP8.5). The largest losses will be at upland drainage locations of water diversions that were dominated historically by glacier melt and seasonal snowmelt. The contribution of glacier melt varies greatly not only in space but also in time. It displays a strong decadal scale fluctuations that are super‐imposed on the effects of a long‐term climatic warming trend. This decadal variability results in reversals in trends in glacier melt, which underscore the importance of long‐time series of glacio‐hydrologic analyses for evaluating the hydrological response to glacier recession. Copyright © 2016 John Wiley & Sons, Ltd.     

Till deposition by glacier submarginal,incremental thickening

Macro‐ and micro‐scale sedimentological analyses of recently deposited tills and complex push/squeeze moraines on the forelands of Icelandic glaciers and in a stacked till sequence at the former Younger Dryas margin of the Loch Lomond glacier lobe in Scotland are used to assess the depositional processes involved in glacier submarginal emplacement of sediment. Where subglacial meltwater is unable to flush out subglacial sediment or construct thick debris‐rich basal ice by cumulative freeze‐on processes, glacier submarginal processes are dictated by seasonal cycles of refreezing and melt‐out of tills advected from up‐ice by a combination of lodgement, deformation and ice keel and clast ploughing. Although individual till layers may display typical A and B horizon deformation characteristics, the spatially and temporally variable mosaic of subglacial processes will overprint sedimentary and structural signatures on till sequences to the extent that they would be almost impossible to classify genetically in the ancient sediment record. At the macro‐scale, Icelandic tills display moderately strong clast fabrics that conform to the ice flow directions documented by surface flutings; very strong fabrics typify unequivocally lodged clasts. Despite previous interpretations of these tills as subglacial deforming layers, micro‐morphological analysis reveals that shearing played only a partial role in the emplacement of till matrixes, and water escape and sediment flowage features are widespread. A model of submarginal incremental thickening is presented as an explanation of these data, involving till slab emplacement over several seasonal cycles. Each cycle involves: (1) late summer subglacial lodgement, bedrock and sediment plucking, subglacial deformation and ice keel ploughing; (2) early winter freeze‐on of subglacial sediment to the thin outer snout; (3) late winter readvance and failure along a decollement plane within the till, resulting in the carriage of till onto the proximal side of the previous year's push moraine; (4) early summer melt‐out of the till slab, initiating porewater migration, water escape and sediment flow and extrusion. Repeated reworking of the thin end of submarginal till wedges produces overprinted strain signatures and clast pavements. Copyright © 2005 John Wiley & Sons, Ltd.     

Depositional features and transportation mechanism of valley-filling Iwasegawa and Kaida debris avalanches, Japan

 The depositional features of two valley-filling debris avalanche deposits were studied to reveal their transportation and depositional mechanisms. The valley-filling Iwasegawa debris avalanche deposit (ca. 0.1 km³) is distributed along the valleys at the southeastern foot of Tashirodake Volcano, northern Honshu, Japan. Debris-avalanche blocks range in size from <35 m proximally to <10 m in the distal zone and consist dominantly of fragile materials. Debris-avalanche matrix percentages increase from 35–60% in the proximal zone to 95% in the distal zone. The debris-avalanche matrix is greater in volume (80–90%) at the bottom and margins of the deposit. Normal grading of large clasts and reverse grading of wood logs and branches occur within the debris-avalanche matrix. Preferred orientation of 311 wood logs and branches within the deposit coincide with the interpreted local flow direction. The basal part of the deposit is characterized by (1) erosional features and incorporated clasts of underlying material; (2) a higher proportion (30–50%) of incorporated clasts than the upper part; and (3) reverse grading of clasts. The valley-filling Kaida debris avalanche deposit (50 000 y B.P., >0.3 km³) is distributed along the valleys at the eastern-southeastern foot of Ontake Volcano, central Japan. Debris-avalanche blocks range in size from <25 m proximally to <7 m in the medial zone. Debris-avalanche matrix percentages increase from 50–70% in the proximal zone to 80% in the distal zone. The debris-avalanche matrix is more abundant (80–90%) at the bottom part of the deposit. Deformation structures observed in the debris-avalanche blocks include elongation, folding, conjugate reverse faults, and numerous minor faults in unconsolidated materials. Lithic components within the debris-avalanche matrix tend to have a higher percentage of plucked clasts from the adjacent underlying formations. A Bingham "plug flow" model is consistent with the transportation and depositional mechanisms of the valley-filling debris avalanches. In the plug of the debris avalanche, fragile blocks were transported without major rupturing due to relatively small shear stresses in regions of small strain rate. The debris-avalanche matrix was mainly produced by shearing at the bottom and margins of the avalanche. Valley-filling debris avalanches tend to have smaller debris-avalanche blocks and larger amounts of debris-avalanche matrix than do unconfined debris avalanches. These differences may be due to disaggregation of debris-avalanche blocks by shearing against valley walls and interaction between debris-avalanche blocks and valley walls. Oriented wood logs and branches, reverse grading of clasts at the base, and a higher proportion of incorporated clasts at the base are interpreted to result from shearing along the bottom and valley walls. Received: 25 March 1998 / Accepted: 10 October 1998     

Heterogeneous water storage and thermal regime of supraglacial ponds on debris‐covered glaciers

The water storage and energy transfer roles of supraglacial ponds are poorly constrained, yet they are thought to be important components of debris‐covered glacier ablation budgets. We used an unmanned surface vessel (USV) to collect sonar depth measurements for 24 ponds to derive the first empirical relationship between their area and volume applicable to the size distribution of ponds commonly encountered on debris‐covered glaciers. Additionally, we instrumented nine ponds with thermistors and three with pressure transducers, characterizing their thermal regime and capturing three pond drainage events. The deepest and most irregularly‐shaped ponds were those associated with ice cliffs, which were connected to the surface or englacial hydrology network (maximum depth = 45.6 m), whereas hydrologically‐isolated ponds without ice cliffs were both more circular and shallower (maximum depth = 9.9 m). The englacial drainage of three ponds had the potential to melt ~100 ± 20 × 10³ kg to ~470 ± 90 × 10³ kg of glacier ice owing to the large volumes of stored water. Our observations of seasonal pond growth and drainage with their associated calculations of stored thermal energy have implications for glacier ice flow, the progressive enlargement and sudden collapse of englacial conduits, and the location of glacier ablation hot‐spots where ponds and ice cliffs interact. Additionally, the evolutionary trajectory of these ponds controls large proglacial lake formation in deglaciating environments. Copyright © 2017 John Wiley & Sons, Ltd.     

Subglacial observations from øksfjordjøkelen,North Norway

Stepped bedrock topography at the snout of a small outlet glacier from Øksfjordjøkelen, North Norway, produces an extensive subglacial cavity system which stretches some 70m across and 100m up-glacier, giving access beneath ice ≤50 m thick. Inside the cavity, regelation ice, clean glacier ice and deforming basal ice have been observed. Samples were taken and basal debris concentrations at the glacier sole were found to vary between 0.005 and 15.38 per cent by volume. The basal ice velocity has been determined using a linear variable differential transformer attached to an analogue recorder, and also by means of measured displacements of ice crews and clasts embedded in the basal ice. Velocities were found to differ both spatially and temporally from a maximum of 2.55 mm h¹ to a minimum of 0.3 mm h^?1. The measurements and observations, which have been related to present theory, show how spatially averaged values for a number of variables could lead to inaccuracies in predicted erosion values, certainly at a local scale. On the exposed foreland, jointcontrolled lee-side faces provide evidence for extensive subglacial plucking (here taken to mean the removal of preloosened bed material and/or material resulting from bed failure). Indeed, in the cavity the early stages of removal of joint-controlled blocks by ice deformation along joints have been observed. The importance of debris-rich basal ice is shown in the formation of large striations (up to 500cm × 16cm × 2cm) present on the foreland.     

Debris cover on thaw slumps and its insulative role in a warming climate

Thaw slumps in ice-rich permafrost can retreat tens of metres per summer, driven by the melt of subaerially exposed ground ice. However, some slumps retain an ice-veneering debris cover as they retreat. A quantitative understanding of the thermal regime and geomorphic evolution of debris-covered slumps in a warming climate is largely lacking. To characterize the thermal regime, we instrumented four debris-covered slumps in the Canadian Low Arctic and developed a numerical conduction-based model. The observed surface temperatures >20° C and steep thermal gradients indicate that debris insulates the ice by shifting the energy balance towards radiative and turbulent losses. After the model was calibrated and validated with field observations, it predicted sub-debris ice melt to decrease four-fold from 1.9 to 0.5 mas the thickness of the fine-grained debris quadruples from 0.1 to 0.4 m. With warming temperatures, melt is predicted to increase most rapidly, in relative terms, for thick (∼0.5–1.0 m) debris covers. The morphology and evolution of the debris-covered slumps were characterized using field and remote sensing observations, which revealed differences in association with morphology and debris composition. Two low-angle slumps retreated continually despite their persistent fine-grained debris covers. The observed elevation losses decreased from ∼1.0 m/yr where debris thickness ∼0.2 mto 0.1 m/yr where thickness ∼1.0 m. Conversely, a steep slump with a coarse-grained debris veneer underwent short-lived bursts of retreat, hinting at a complex interplay of positive and negative feedback processes. The insulative protection and behaviour of debris vary significantly with factors such as thickness, grain size and climate: debris thus exerts a fundamental, spatially variable influence on slump trajectories in a warming climate. © 2020 John Wiley & Sons, Ltd.     

Flow prediction in high altitude ungauged catchments: A case study in the Italian Alps (Pantano Basin,Adamello Group)

We develop a simple model to evaluate the daily flow discharges in the ablation season for the 11 km² Pantano basin in the Retiche Italian Alps, based upon the data gathered during a three years field campaign. The Pantano basin embeds the Venerocolo debris covered and the Avio debris free glaciers, covering 2.14 km² in the Adamello Group, where the widest Italian glacier Adamello is located. First, degree-day models based upon air temperature are tuned to calculate snow and ice melt at daily scale. Glaciers’ meteorological data are collected from an automatic weather station (AWS), operating on the glacier during summer 2007. The melt factors in the debris covered areas of the glacier are estimated against debris thickness, using a data driven parameterization. The flow discharge from the catchment is estimated using semi distributed flow routing for the ablation seasons of four years, from 2006 to 2009. The predicted discharges are compared to those derived from inverse reservoir's routing at the Benedetto lake, catching the basin outflow. The proposed approach is valuable as a tool to investigate the hydrology of poorly gauged glacierized areas, including those with debris covered ice, widely diffused and yet poorly understood. Pending accurate parameterization the approach is usable for water resources evaluation and for long term assessment of the climate change impact on the glacierized areas within the Alps.     

Interrogating the time and processes of development of the Las Liebres rock glacier,central Chilean Andes,using a numerical flow model

The Las Liebres rock glacier is a large (~2.2 km long) Andean rock glacier whose internal composition and kinematics are known from previous studies. We investigate its development by posing and testing the following null hypothesis: the rock glacier has developed from a constant supply of debris and ground ice in periglacial conditions and resulting creep of the ice‐rock mixture. A rheological model was formulated based on recent advances in the study of ice‐rock mixture rheology, and calibrated on the known surface velocities and internal composition of the rock glacier. We show that the rock glacier viscosity is inversely related to both water and debris fractions, in agreement with recent field and theoretical studies of ice‐rock mixture rheology. Taking into account the possible variations in water fraction, the model was used to estimate the time spans of development (0.91–7.11 ka), rates of rock wall retreat (0.44–4.18 mm/a), and rates of ground ice formation (0.004–0.026 m/a) for the rock glacier. These results support the null hypothesis of a periglacial origin of the Las Liebres rock glacier. Copyright © 2016 John Wiley & Sons, Ltd.     

Debris entrainment by basal freeze‐on and thrusting during the 1995–1998 surge of Kuannersuit Glacier on Disko Island,west Greenland

Kuannersuit Glacier, a valley glacier on Disko Island in west Greenland, experienced a major surge from 1995 to 1998 where the glacier advanced 10·5 km and produced a ～65 m thick stacked sequence of debris‐rich basal ice and meteoric glacier ice. The aim of this study is to describe the tectonic evolution of large englacial thrusts and the processes of basal ice formation using a multiproxy approach including structural glaciology, stable isotope composition (δ¹⁸O and δD), sedimentology and ground‐penetrating radar. We argue that the major debris layers that can be traced in the terminal zone represent englacial thrusts that were formed early during the surge. Thrust overthrow was at least 200–300 m and this lead to a 30 m thick repetition of basal ice at the ice margin. It is assumed that the englacial thrusting was initiated at the transition between warm ice from the interior and the cold snout. The basal debris‐rich ice was mainly formed after the thrusting phase. Two sub‐facies of stratified basal ice have been identified; a lower massive ice facies (S_M) composed of frozen diamict enriched with heavy stable isotopes overlain by laminated ice facies (S_L) consisting of millimetre thick lamina of alternating debris‐poor and debris‐rich ice. We interpret the stratified basal ice as a continuum formed mainly by freeze‐on processes and localized regelation. First laminated basal ice is formed and as meltwater is depleted more sediment is entrained and finally the glacier freezes to the base and massive diamict is frozen‐on. The increased ability to entrain sediments may partly be associated with higher basal freezing rates enhanced by loss of frictional heat from cessation of fast flow and conductive cooling through a thin heavily crevassed ice during the final phase of the glacier surge. The dispersed basal ice facies (D) was mainly formed by secondary processes where fine‐grained sediment is mobilized in the vein system of ice. Our results have important implications for understanding the significance of basal ice formation and englacial thrusting beneath fast‐flowing glaciers and it provides new information about the development of landforms during a glacier surge. Copyright © 2010 John Wiley & Sons, Ltd.     

Stability of supraglacial debris

Rock debris on the surface of ablating glaciers is not static, and is often transported across the ice surface as relief evolves during melt. This supraglacial debris transport has a strong influence on the spatial distribution of melt, and is implicated in the formation of hummocky glacial topography in deglaciated terrain. Furthermore, as ice‐dammed lakes and ice‐cored slopes become increasingly common in deglaciating watersheds, there is rising concern about hazards to humans and infrastructure posed by mass‐wasting of ice‐cored debris. The existing quantitative framework for describing these debris transport processes is limited, making it difficult to account for transport in mass balance, hazard assessment, and landscape development models. This paper develops a theoretical framework for assessing slope stability and gravitational mass transport in a debris‐covered ice setting. Excess water pressure at the interface between ablating ice and lowering debris is computed by combining Darcy's law with a meltwater balance. A limit‐equilibrium slope stability analysis is then applied to hypothetical debris layers with end‐member moisture conditions derived from a downslope meltwater balance that includes production and seepage. The resulting model system constrains maximum stable slope angles and lengths that vary with debris texture, thickness, and the rate of meltwater production. Model predictions are compared with field observations and with digital elevation model (DEM)‐derived terrain metrics from two modern debris‐covered glaciers on Mount Rainier, USA. Copyright © 2017 John Wiley & Sons, Ltd.