Some recent studies revealed that Himalayan glaciers were shrinking at an accelerated rate since the beginning of the 21st century. However, the climatic causes for this shrinkage remain unclear given that surface energy balance studies are almost nonexistent in this region. In this study, a point-scale surface energy balance analysis was performed using in situ meteorological data from the ablation zone of Chhota Shigri Glacier over two separate periods (August 2012 to February 2013 and July to October 2013) in order to understand the response of mass balance to climatic variables. Energy balance numerical modelling provides quantification of the surface energy fluxes and identification of the factors affecting glacier mass balance. The model was validated by comparing the computed and observed ablation and surface temperature data. During the summer-monsoon period, net radiation was the primary component of the surface energy balance accounting for 80 % of the total heat flux followed by turbulent sensible (13%), latent (5%) and conductive (2%) heat fluxes. A striking feature of the energy balance is the positive turbulent latent heat flux, suggesting re-sublimation of moist air at the glacier surface, during the summer-monsoon characterized by relatively high air temperature, high relative humidity and a continual melting surface. The impact of the Indian Summer Monsoon on Chhota Shigri Glacier mass balance has also been assessed. This analysis demonstrates that the intensity of snowfall events during the summer-monsoon plays a key role on surface albedo (melting is reduced in the case of strong snowfalls covering the glacier area), and thus is among the most important drivers controlling the annual mass balance of the glacier. The summer-monsoon air temperature, controlling the precipitation phase (rain versus snow and thus albedo), counts, indirectly, also among the most important drivers.

We present the updated glaciological mass balance (MB) of Chhota Shigri Glacier, the longest continuous annual MB record in the Hindu-Kush Karakoram Himalaya (HKH) region. Additionally, 4 years of seasonal MBs are presented and analyzed using the data acquired at an automatic weather station (AWS-M) installed in 2009 on a lateral moraine (4863ma.s.l.). The glaciological MB series since 2002 is first recalculated using an updated glacier hypsometry and then validated against geodetic MB derived from satellite stereo-imagery between 2005 (SPOT5) and 2014 (Pléiades). Chhota Shigri Glacier lost mass between 2002 and 2014 with a cumulative glaciological MB of –6.72mw.e. corresponding to a mean annual glacier-wide MB (B_a) of –0.56mw.e. a^–1. Equilibrium-line altitude (ELA₀) for the steady-state condition is calculated as ~4950ma.s.l., corresponding to an accumulation–area ratio (AAR₀) of ~61%. Analysis of seasonal MBs between 2009 and 2013 with air temperature from AWS-M and precipitation from the nearest meteorological station at Bhuntar (1050ma.s.l.) suggests that the summer monsoon is the key season driving the interannual variability of B_a for this glacier. The intensity of summer snowfall events controls the B_a evolution via controlling summer glacier-wide MB (B_s).

The Himalayan glaciers contribute significantly to regional water resources. However, limited field observations restrict our understanding of glacier dynamics and behaviour. Here, we investigated the long-term in situ mass balance, meteorology, ice velocity and discharge of the Chhota Shigri Glacier. The mean annual glacier-wide mass balance was negative, −0.46 ± 0.40 m w.e. a⁻¹ for the period 2002–2019 corresponding to a cumulative wastage of −7.87 m w.e. Winter mass balance was 1.15 m w.e. a⁻¹ and summer mass balance was −1.35 m w.e. a⁻¹ over 2009–2019. Surface ice velocity has decreased on average by 25–42% in the lower and middle ablation zone (below 4700 m a.s.l.) since 2003; however, no substantial change was observed at higher altitudes. The decrease in velocity suggests that the glacier is adjusting its flow in response to negative mass balance. The summer discharge begins to rise from May and peaks in July, with a contribution of 43%, followed by 38% and 19% in August and September, respectively. The discharge pattern closely follows the air temperature. The long-term observation on the ‘Chhota Shigri – a benchmark glacier’, shows a mass wastage which corresponds to the slowdown of the glacier in the past two decades.

We present the first-ever mass-balance (MB) observation (2014–19), reconstruction (between 1978 and 2019) and sensitivity of debris-free Stok glacier (33.98°N, 77.45°E), Ladakh Region, India. In-situ MB was negative throughout the study period except in 2018/19 when the glacier witnessed a balanced condition. For MB modelling, three periods were considered based on the available data. Period I (1978–87, 1988/89) witnessed a near balance condition (−0.03 ± 0.35 m w.e. a⁻¹) with five positive MB years. Whereas Period II (1998–2002, 2003–09) and III (2011–19) experienced high (−0.9 ± 0.35 m w.e. a⁻¹) and moderate (−0.46 ± 0.35 m w.e. a⁻¹) negative MBs, respectively. Glacier area for these periods was derived from the Corona, Landsat and PlanetScope imageries using a semi-automatic approach. The in-situ and modelled MBs were in good agreement with RMSE of 0.23 m w.e. a⁻¹,R²= 0.92,P< 0.05. The average mass loss was moderate (−0.47 ± 0.35 m w.e. a⁻¹) over 28 hydrological years between 1978 and 2019. Sensitivity analysis showed that the glacier was more sensitive to summer temperature (−0.32 m w.e. a⁻¹°C⁻¹) and winter precipitation (0.12 m w.e. a⁻¹for ± 10%). It was estimated that ~27% increase in precipitation is required on Stok glacier to compensate for the mass loss due to 1°C rise in temperature.

Understanding changes in glacier mass balances is essential for investigating climate changes. However, glacier‐wide mass balances determined from geodetic observations do not provide a relevant climatic signal as they depend on the dynamic response of the glaciers. In situ point mass balance measurements provide a direct signal but show a strong spatial variability that is difficult to assess from heterogeneous in situ measurements over several decades. To address this issue, we propose a nonlinear statistical model that takes into account the spatial and temporal changes in point mass balances. To test this model, we selected four glaciers in different climatic regimes (France, Bolivia, India, and Norway) for which detailed point annual mass balance measurements were available over a large elevation range. The model extracted a robust and consistent signal for each glacier. We obtained explained variances of 87.5, 90.2, 91.3, and 75.5% on Argentière, Zongo, Chhota Shigri, and Nigardsbreen glaciers, respectively. The standard deviations of the model residuals are close to measurement uncertainties. The model can also be used to detect measurement errors. Combined with geodetic data, this method can provide a consistent glacier‐wide annual mass balance series from a heterogeneous network. This model, available to the whole community, can be used to assess the impact of climate change in different regions of the world from long‐term mass balance series.