Decadal-scale, high-resolution geodetic measurements of glacier thinning have transformed our understanding of glacier response to climate change. Annual glacier mass balance can be estimated using remote-sensing proxies like snow-line altitude. These methods require field data for calibration, which are not available for most glaciers. Here we propose a method that combines multiple remotely-sensed proxies to obtain robust estimates of the annual glacier-wide balance using only remotely-sensed decadal-scale geodetic mass balance for calibration. The method is tested on Chhota Shigri, Argentière and Saint-Sorlin glaciers in the Himalaya and the Alps between 2001 and 2020, using four remotely-sensed proxies – the snow-line altitude, the minimum summer albedo over the glacier and two statistics of normalised difference snow index over the off-glacier area around the ablation zone. The reconstructed mass balance compares favourably with the corresponding glaciological field data (correlation coefficient 0.81 − 0.90, p < 0.001; root mean squared error 0.38 − 0.43 m w.e. a−1). The method presented may be useful to study interannual variability in mass balance on glaciers where no field data are available.

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.

A volume-area scaling relation is commonly used to estimate glacier volume or its future changes on a global scale. The presence of an insulating supraglacial debris cover alters the mass-balance profile of a glacier, potentially modifying the scaling relation. Here, the nature of scaling relations for extensively debris-covered glaciers is investigated. Theoretical arguments suggest that the volume-area scaling exponent for these glaciers is ~7% smaller than that for clean glaciers. This is consistent with the results from flowline-model simulations of idealised glaciers, and the available data from the Himalaya. The best-fit scale factor for debris-covered Himalayan glaciers is ~60% larger compared to that for the clean ones, implying a significantly larger stored ice volume in a debris-covered glacier compared to a clean one having the same area. These results may help improve scaling-based estimates of glacier volume and future glacier changes in regions where debris-covered glaciers are abundant.

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.

Avalanching is a prominent source of accumulation on glaciers that have high and steep valley-walls surrounding their accumulation zones. These glaciers are typically characterised by an extensive supraglacial debris cover and a low accumulation area ratio. Despite an abundance of such glaciers in the rugged landscapes of the High Himalaya, attempts to quantify the net avalanche contribution to mass balance and its long-term variation are almost missing. We first discuss diagnostic criteria to identify strongly avalanche-fed glaciers. Second, we develop an approximate method to quantify the magnitude of the avalanche accumulation exploiting its expected control on the dynamics of these glaciers. The procedure is based on a simplified flowline model description of the glacier concerned and utilises the known glaciological mass-balance, velocity and surface-elevation profiles of the glacier. We apply the method to three Himalayan glaciers and show that the data on the recent dynamics of these glaciers are consistent with a dominant contribution of avalanches to the total accumulation. As a control experiment, we also simulate another Himalayan glacier where no significant avalanche contribution is expected, and reproduce the recent changes in that glacier without any additional avalanche contribution.

Modelling the response of Himalayan glaciers to rapid climate change is an important problem. The poorly understood effects of debris cover and the variable response of the glaciers have made it difficult to understand their dynamics. We propose a simple model for debris-covered glaciers and validate it against data from Dokriani Glacier, India. Numerical investigations of the model show that the response of debris-covered glaciers to a warming climate has two timescales. There is a period when the glacier loses ice by thinning but the front is almost stationary and it develops a long, slow-flowing tongue. This stationary period, which can be >100 years for glaciers with a large extent of debris cover, is negligible for bare glaciers. The quasi-stagnant tongue does not develop in response to cooling. An analysis of remote-sensing data in the light of these results indicates that the variable response of the glaciers in the Himalaya is consistent with a climate that is warming on average, but has considerable spatial variability in the warming rates. We estimate the average warming rate to be about the same as the global average.

A temperature reconstruction in the glacierized Himalaya over the past centuries using glacial length fluctuation records is challenging due to the abundance of debris-covered glaciers and a scarcity of glacial length fluctuation data. Using idealized flowline model simulations, we show that supraglacial debris cover significantly alters the length fluctuations only when the debris cover is very thick. An expanded database of length fluctuation records for 43 glaciers in the Himalaya and Karakoram is compiled and a standard linear inversion procedure is applied to a subset of 34 glaciers in this database. The reconstructed temperature anomaly during 1860-2010 indicates a continued warming of the region with a total temperature change of ~1.6 K. A close resemblance of the regional temperature anomaly to global trends is seen.

We have compiled and analysed available records and data on the shrinkage of Satopanth (SPG) and Bhagirath Kharak (BKG) Glaciers, Uttarakhand, India, during the period 1936–2013. We estimate the mean retreat rates of the snouts of SPG and BKG for this period at 9.7±0.8m a–1 and 7.0±0.6m a–1 respectively. We have also revised the estimates of the area vacated during the period 1956–2013 to be 0.27 ± 0.05 km2 and 0.17 ± 0.04 km2 for SPG and BKG respectively, corresponding to front-averaged retreat rates of 5.7±0.6m a–1 and 6.0±0.9m a–1. The study revealed an average thinning of glacial ice in the lower ablation zone of SPG of 9 ± 11 m in the past 51 years. We observed that while the fronts of SPG and BKG depicted in the Survey of India topographic map published in 1962 are inconsistent with other available records, the elevation contours are consistent with them.