Debris-covered ice is widespread in mountain regions with debris an important control on surface ice melt and glacier retreat. Quantifying debris cover extent and its evolution through time over large regions remains a challenge. This study develops two Normalized Difference Supraglacial Debris Indices for mapping debris-covered ice based on thermal and near Infrared Landsat 8 bands. They were calibrated with field data. Validation suggests that they have a high level of accuracy. They are then applied to Landsat data for 2016 to produce the first detailed glacier inventory of the Afghanistan Hindu Kush Himalaya that includes debris cover. 3408 glaciers were identified which, for those ⩾0.01 km2 in area, gives an ice cover of 2,222 ± 11 km2 and a debris cover of 619 ± 40 km2. Principal components analysis was used to identify the most influential drivers of debris-covered ice extent. Lower proportions of debris cover were associated with glaciers with a higher elevation range, that were larger, longer and wider. These relations were statistically clearer when the dataset was broken down into climate and geological zones. A glaciers continue to shrink, the proportion of debris cover will become higher, making it more important to map debris cover reliably.

We assess the composition and geometry of four individual rock glaciers in Alaska, Wyoming and Colorado by measuring their radio wave speed and applying these results to ground-penetrating radar depth corrections and dielectric mixing models. Our method includes a correction for subsurface reflector dip angle, which we show can lead to an incorrect determination of wave speeds using common midpoint configurations. By observing the radar properties of the rock glaciers and their supraglacial debris, we find that some of the sites exhibit nearly pure ice cores, and all of the sites indicate volumetric ice fractions >50%. These results have implications for terrestrial glaciology and hydrology because the present ice volume is connected to past ice accumulation and subsurface ice preservation, which may affect the future availability of alpine water resources. An understanding of the processes that govern rock glacier evolution over a wide range of latitudes and elevations will also contribute to the exploration of planetary surfaces such as Mars, which hosts a significant population of debris-covered glaciers. Our subsurface composition and geometry estimates will inform simulations of rock glacier formation and evolution to test hypothesized ice origin mechanisms along with the preservation of climate signals.

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

A supraglacial debris layer controls energy transfer to the ice surface and moderates ice ablation on debris-covered glaciers. Measurements of vertical temperature profiles within the debris enables the estimation of thermal diffusivities and sub-debris ablation rates. We have measured the debris-layer temperature profiles at 16 locations on Satopanth Glacier (central Himalaya) during the ablation seasons of 2016 and 2017. Debris temperature profile data are typically analysed using a finite-difference method, assuming that the debris layer is a homogeneous one-dimensional thermal conductor. We introduce three more methods for analysing such data that approximate the debris layer as either a single or a two-layered conductor. We analyse the performance of all four methods using synthetic experiments and by comparing the estimated ablation rates with in situ glaciological observations. Our analysis shows that the temperature measurements obtained at equispaced sensors and analysed with a two-layered model improve the accuracy of the estimated thermal diffusivity and sub-debris ablation rate. The accuracy of the ablation rate estimates is comparable to that of the in situ observations. We argue that measuring the temperature profile is a convenient and reliable method to estimate seasonal to sub-seasonal variations of ablation rates in the thickly debris-covered parts of glaciers.

Surface energy-balance models are commonly used in conjunction with satellite thermal imagery to estimate supraglacial debris thickness. Removing the need for local meteorological data in the debris thickness estimation workflow could improve the versatility and spatiotemporal application of debris thickness estimation. We evaluate the use of regional reanalysis data to derive debris thickness for two mountain glaciers using a surface energy-balance model. Results forced using ERA-5 agree with AWS-derived estimates to within 0.01 ± 0.05 m for Miage Glacier, Italy, and 0.01 ± 0.02 m for Khumbu Glacier, Nepal. ERA-5 data were then used to estimate spatiotemporal changes in debris thickness over a ~20-year period for Miage Glacier, Khumbu Glacier and Haut Glacier d'Arolla, Switzerland. We observe significant increases in debris thickness at the terminus for Haut Glacier d'Arolla and at the margins of the expanding debris cover at all glaciers. While simulated debris thickness was underestimated compared to point measurements in areas of thick debris, our approach can reconstruct glacier-scale debris thickness distribution and its temporal evolution over multiple decades. We find significant changes in debris thickness over areas of thin debris, areas susceptible to high ablation rates, where current knowledge of debris evolution is limited.

Supraglacial debris is significant in many regions and complicates modeling of glacier melt, which is required for predicting glacier change and its influences on hydrology and sea-level rise. Temperature-index models are a popular alternative to energy-balance models when forcing data are limited, but their transferability among glaciers and inherent uncertainty have not been documented in application to debris-covered glaciers. Here, melt factors were compiled directly from published studies or computed from reported melt and MERRA-2 air temperature for 27 debris-covered glaciers around the world. Linear mixed-effects models were fit to predict melt factors from debris thickness and variables including debris lithology and MERRA-2 radiative exchange. The models were tested by leave-one-site-out cross-validation based on predicted melt rates. The best model included debris thickness (fixed effect) and glacier and year (random effects). Predictions were more accurate using MERRA-2 than on-site air temperature data, and pooling MERRA-2-derived and reported melt factors improved cross-validation accuracy more than including additional predictors such as shortwave or longwave radiation. At one glacier where monthly ablation was measured over 4 years, seasonal variation of melt factors suggested that heat storage significantly affected the relation between melt and energy exchange at the debris surface.

The ice-cored Galena Creek Rock Glacier, Wyoming, USA, has been the subject of a number of studies that sought to determine the origin of its ice. We present new observations of the rock glacier's internal structure from ground-penetrating radar to constrain ice and debris distribution and accumulation. We imaged dipping reflectors in the center of the glacier that are weak and discontinuous, in contrast to strong reflectors toward the edge of the cirque beneath large debris-avalanche chutes. These reflectors form a network of concave-up, up-glacier dipping layers. We interpret these as englacial debris bands formed by large debris falls buried by subsequent ice and snow accumulation. They are discontinuous where ice outpaces debris accumulation, but with sufficient debris accumulation an interleaved pattern of ice and debris layers can form. We propose a model in which the ice in these interleaved layers is snowfall preserved by debris-facilitated accumulation. Large debris falls that occur in early spring bury sections of the snowpack, which are then preserved through summer and incorporated into the rock glacier body over time. This study highlights the importance of sequential accumulation of ice and debris for understanding the dynamics of rock glaciers and debris-covered glaciers.

Glaciological ablation is computed from point-scale data at a few ablation stakes that are usually regressed as a function of elevation and averaged over the area-elevation distribution of a glacier. This method is contingent on a tight control of elevation on local ablation. However, in debris-covered glaciers, systematic and random spatial variations of debris thickness modify the ablation rates. We propose and test a method to compute sub-debris ablation where stake data are interpolated as a function of debris-thickness alone and averaged over the debris-thickness distribution at different parts of the glacier. We apply this method on Satopanth Glacier located in Central Himalaya utilising ~1000 ablation measurements obtained from a network of up to 56 stakes during 2015–2017. The estimated mean sub-debris ablation ranges between 1.5±0.2 to 1.7±0.3 cm d−1. We show that the debris-thickness-dependent regression describes the spatial variability of the sub-debris ablation better than the elevation dependent regression. The uncertainties in ablation estimates due to the corresponding uncertainties in the measurement of ablation and debris-thickness distribution, and those due to interpolation procedures are estimated using Monte Carlo methods. Possible biases due to a finite number of stakes used are also investigated.

Surface debris covers much of the western portion of the McMurdo Ice Shelf and has a strong influence on the local surface albedo and energy balance. Differential ablation between debris-covered and debris-free areas creates an unusual heterogeneous surface of topographically low, high-ablation, and topographically raised (‘pedestalled’), low-ablation areas. Analysis of Landsat and MODIS satellite imagery from 1999 to 2018, alongside field observations from the 2016/2017 austral summer, shows that pedestalled relict lakes (‘pedestals’) form when an active surface meltwater lake that develops in the summer, freezes-over in winter, resulting in the lake-bottom debris being masked by a high-albedo, superimposed, ice surface. If this ice surface fails to melt during a subsequent melt season, it experiences reduced surface ablation relative to the surrounding debris-covered areas of the ice shelf. We propose that this differential ablation, and resultant hydrostatic and flexural readjustments of the ice shelf, causes the former supraglacial lake surface to become increasingly pedestalled above the lower topography of the surrounding ice shelf. Consequently, meltwater streams cannot flow onto these pedestalled features, and instead divert around them. We suggest that the development of pedestals has a significant influence on the surface-energy balance, hydrology and flexure of the ice shelf.

To understand the geological origins of minerals in cryoconite and the nutrients sources for microbes on glaciers, we analyzed the Sr–Nd isotopic ratios of the four mineral fractions in cryoconites including saline, carbonate, phosphate, silicate and the organic fraction obtained from Gulkana Glacier in Alaska. The isotopes in the silicate mineral fraction exhibited spatial variation within the glacier (87Sr/86Sr: 0.704533–0.709563, εNd (0): −16.0 to 0.5), which can be explained by the different mixing ratios of the two distinct sources: one of the sources is lateral and terminal moraines or soil, and the other is the medial moraine of the glacier. The minerals in the cryoconite at the lower sites in the glacier are likely derived from the former source, whereas those at the upper sites are from latter sources. The mineralogical and elemental compositions also support mixing of the silicate minerals from the two local sources. The Sr isotopic ratios of the organic fraction also showed spatial variation on the glacier in the middle sites – a trend similar to those of the phosphate fraction. The results suggest that the organic matter is mostly the byproducts of microbes using the phosphate minerals as a nutrient source.

Light-absorbing impurities (LAI) can darken snow and ice surfaces, reduce snow/ice albedo and accelerate melt. Efforts to allocate the relative contribution of different LAI to snow/ice albedo reductions have been limited by uncertainties in the optical properties of LAI. We developed a new method to measure LAI spectral reflectance at the submicron scale by modifying a Hyperspectral Imaging Microscope Spectrometer (HIMS). We present the instrument's internal calibration, and the overall small influence of a particle's orientation on its measured reflectance spectrum. We validated this new method through the comparison with a field spectroradiometer by measuring different standard materials. Measurements with HIMS at the submicron scale and the bulk measurements of the same standard materials with the field spectroradiometer are in good agreement with an average deviation between the spectra of 3.2% for the 400–1000 nm wavelength range. The new method was used (1) to identify BC (black carbon), mineral dust including hematite and the humic substances present in an environmental sample from Plaine Morte glacier and (2) to collect the individual reflectance spectra of each of these types of impurity. The results indicate that this method is applicable to heterogeneous samples such as the LAI found in snow and ice.

The thickness of supraglacial debris cover controls how it impacts the ablation rate of underlying glacier ice, yet this quantity remains challenging to measure, particularly at glacier scales. We present a relatively straightforward, and cost-effective method to estimate debris thickness exposed above ice cliffs using simplified geometrical measurements from a high-resolution digital surface model (DSM), derived from a terrestrial photographic survey and a Structure from Motion with Multi-View Stereo workflow (SfM-MVS). As the ice surface relief beneath the debris cover is unknown, we assume it to be horizontal and provide error bounds based on characteristic ice-surface slope at the visible debris/ice interface. Debris thickness around the three sampled ice cliffs was highly variable (interquartile range of 0.80–2.85 m) and negatively skewed with a mean thickness of 2.08 ± 0.68 m. Manual, and high-frequency radar, determinations of debris thickness in the same area show similar thickness distributions, but statistically different mean debris thickness, due to local heterogeneity. Debris thickness values derived in this study all exceed estimates from satellite surface temperature inversions. Wider application of the method presented here would provide useful data for improving debris thickness approximations from satellite imagery.

Supraglacial debris thickness is a key control on the surface energy balance of debris-covered glaciers, yet debris thickness measurements are sparse due to difficulties of data collection. Here we use ground-penetrating radar (GPR) to measure debris thickness on the ablation zone of Lirung Glacier, Nepal. We observe a strong, continuous reflection, which we interpret as the ice surface, through debris layers of 0.1 to at least 2.3 m thick, provided that appropriate acquisition parameters were used while surveying. GPR measurements of debris thickness correlate well with pit measurements of debris thickness (r = 0.91, RMSE = 0.04 m) and two-way travel times are consistent at tie points (r = 0.97). 33% of measurements are <0.5 m, so sub-debris melting is likely important in terms of mass loss on the debris-covered tongue and debris thickness is highly variable over small spatial scales (of order 10 m), likely due to local slope processes. GPR can be used to make debris thickness measurements more quickly, over a wider debris thickness range, and at higher spatial resolution than any other means and is therefore a valuable tool with which to map debris thickness distribution on Himalayan glaciers.

A history and conception of glacial events for the central Karakoram Himalaya, proposed some 80 years ago by Giotto Dainelli, are largely accepted today. However, certain deposits identified as terminal moraine complexes marking glacial episodes were actually emplaced by rock avalanches. In the Skardu and Shigar intermontane basins of Baltistan, at least 15 rock avalanche events were previously mapped as moraine or till. Criteria used for distinguishing these catastrophic landslide deposits emphasize homogeneous lithology of rubble and matrix, clast shape, facies characteristics, the large scale unity of emplacement, and morphological relations to valley topography. The deposits of three events, at Katzarah, Satpura, and the north end of Shigar Valley, have been reconstructed in detail. Thick supraglacial debris does not result in similar deposits. Extensive valley fills, river terraces, large sediment fans, and lacustrine sediments formerly attributed to late-glacial conditions are reinterpreted as postglacial events involving rock avalanches that interrupted fluvial development. Existing reconstructions of glaciations are left in doubt, especially late-glacial events in the central Karakoram, as are the roles assigned to Karakoram, main Indus Valley, and western Himalayan ice.

A glacierized terrain comprises different land covers, and their mapping using satellite data is challenged by their spectral similarity. We propose a hierarchical knowledge-based classification (HKBC) approach for differentiation of glacier terrain classes and mapping of glacier boundaries, using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery and Global Digital Elevation Model (GDEM). The methodology was tested over Kolahoi Glacier, Kashmir Himalaya. For the sequential extraction of various glacier terrain classes, several input layers were generated from the primary datasets by applying image-processing techniques. Noticeable differences in temperature and spectral response between supraglacial debris and periglacial debris facilitated the development of a thermal glacier mask and normalized-difference debris index, which together with slope enabled their differentiation. These and the other layers were then used in several discrete tests in HKBC, to map various glacier terrain classes. An ASTER visible near-infrared image and 42 field points were used to validate results. The proposed approach satisfactorily classified all the glacier terrain classes with an overall accuracy of 89%. The Z-test reveals that results obtained from HKBC are significantly (at 95% confidence level) better than those from a maximum likelihood classifier (MLC). Glacier boundaries obtained from HKBC were found to be plausibly better than those obtained from MLC and visual interpretation.

The ice surface of the glacier in the 7 km long valley Scharffenbergbotnen in Heimefrontfjella, Dronning Maud Land, is characterized by blue ice fields and is partly covered by supraglacial moraine complexes. The debris cover of the supraglacial moraines is generally less than 50 cm thick and rests on glacial ice, and most of the surface morphology of ridges and numerous sink holes reflects irregularities in the underlying ice. The debris consists of angular to sub-angular clasts of local bedrock. We interpret the moraines to have an entirely supraglacial history, where colluvium and pre-existing moraines were brought into Scharffenbergbotnen primarily from the south/south-west by an advancing glacier during the last ice maximum in the region. The debris cover of the supraglacial moraines extends up to 200–250 m above the present ice surface on surrounding slopes in Scharffenbergbotnen, and generally to less than 100 m above the present ice surface on slopes outside the valley. We interpret this as marking the ice surface elevation in the area during the last ice maximum (probably LGM). The deposition of supraglacial moraines in the valley at that time, and their survival in the area until the present day, indicate that a local ablation centre and probably blue ice fields were present in Scharffenbergbotnen at LGM, and has been a persistent feature since. In a wider context, we argue that supraglacial moraines constitute a hitherto unexploited resource when reconstructing the former extent of blue ice areas.