Himalayan glaciers are a focus of public and scientific debate. Prevailing uncertainties are of major concern because some projections of their future have serious implications for water resources. Most Himalayan glaciers are losing mass at rates similar to glaciers elsewhere, except for emerging indications of stability or mass gain in the Karakoram. A poor understanding of the processes affecting them, combined with the diversity of climatic conditions and the extremes of topographical relief within the region, makes projections speculative. Nevertheless, it is unlikely that dramatic changes in total runoff will occur soon, although continuing shrinkage outside the Karakoram will increase the seasonality of runoff, affect irrigation and hydropower, and alter hazards.

Glaciers are among the best indicators of terrestrial climate variability, contribute importantly to water resources in many mountainous regions and are a major contributor to global sea level rise. In the Hindu Kush-Karakoram-Himalaya region (HKKH), a paucity of appropriate glacier data has prevented a comprehensive assessment of current regional mass balance. There is, however, indirect evidence of a complex pattern of glacial responses in reaction to heterogeneous climate change signals. Here we use satellite laser altimetry and a global elevation model to show widespread glacier wastage in the eastern, central and south-western parts of the HKKH during 2003-08. Maximal regional thinning rates were 0.66 ± 0.09 metres per year in the Jammu-Kashmir region. Conversely, in the Karakoram, glaciers thinned only slightly by a few centimetres per year. Contrary to expectations, regionally averaged thinning rates under debris-mantled ice were similar to those of clean ice despite insulation by debris covers. The 2003-08 specific mass balance for our entire HKKH study region was -0.21 ± 0.05 m yr(-1) water equivalent, significantly less negative than the estimated global average for glaciers and ice caps. This difference is mainly an effect of the balanced glacier mass budget in the Karakoram. The HKKH sea level contribution amounts to one per cent of the present-day sea level rise. Our 2003-08 mass budget of -12.8 ± 3.5 gigatonnes (Gt) per year is more negative than recent satellite-gravimetry-based estimates of -5 ± 3 Gt yr(-1) over 2003-10 (ref. 12). For the mountain catchments of the Indus and Ganges basins, the glacier imbalance contributed about 3.5% and about 2.0%, respectively, to the annual average river discharge, and up to 10% for the Upper Indus basin.

High Mountain Asia hosts the largest glacier concentration outside the polar regions. These glaciers are important contributors to streamflow in one of the most populated areas of the world. Past studies have used methods that can only provide regionally-averaged glacier mass balances to assess the High Mountain Asia glacier contribution to rivers and sea level rise. Here we compute the mass balance for about 92 % of the glacierized area of High Mountain Asia using time series of digital elevation models derived from satellite stereo-imagery. We calculate an average region-wide mass balance of -16.3 ± 3.5 Gt yr-1 (-0.18 ± 0.04 m w.e. yr-1) between 2000 and 2016, which is less negative than most previous estimates. Region-wide mass balances vary from -4.0 ± 1.5 Gt yr-1 (-0.62 ± 0.23 m w.e. yr-1) in Nyainqentanglha to +1.4 ± 0.8 Gt yr-1 (+0.14 ± 0.08 m w.e. yr-1) in Kunlun, with large intra-regional variability of individual glacier mass balances (standard deviation within a region ˜0.20 m w.e. yr-1). Specifically, our results shed light on the Nyainqentanglha and Pamir glacier mass changes, for which contradictory estimates exist in the literature. They provide crucial information for the calibration of the models used for projections of future glacier response to climatic changes, models that presently do not capture the pattern, magnitude and intra-regional variability of glacier changes in High Mountain Asia.

Abstract. The recent evolution of Pamir-Karakoram-Himalaya (PKH) glaciers, widely acknowledged as valuable high-altitude as well as mid-latitude climatic indicators, remains poorly known. To estimate the region-wide glacier mass balance for 9 study sites spread from the Pamir to the Hengduan Shan (eastern Himalaya), we compared the 2000 Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) to recent (2008–2011) DEMs derived from SPOT5 stereo imagery. During the last decade, the region-wide glacier mass balances were contrasted with moderate mass losses in the eastern and central Himalaya (−0.22 ± 0.12 m w.e. yr−1 to −0.33 ± 0.14 m w.e. yr−1) and larger losses in the western Himalaya (−0.45 ± 0.13 m w.e. yr−1). Recently reported slight mass gain or balanced mass budget of glaciers in the central Karakoram is confirmed for a larger area (+0.10 ± 0.16 m w.e. yr−1) and also observed for glaciers in the western Pamir (+0.14 ± 0.13 m w.e. yr−1). Thus, the "Karakoram anomaly" should be renamed the "Pamir-Karakoram anomaly", at least for the last decade. The overall mass balance of PKH glaciers, −0.14 ± 0.08 m w.e. yr−1, is two to three times less negative than the global average for glaciers distinct from the Greenland and Antarctic ice sheets. Together with recent studies using ICESat and GRACE data, DEM differencing confirms a contrasted pattern of glacier mass change in the PKH during the first decade of the 21st century.

Abstract.There are an increasing number of digital elevation models (DEMs) available worldwide for deriving elevation differences over time, including vertical changes on glaciers. Most of these DEMs are heavily post-processed or merged, so that physical error modelling becomes difficult and statistical error modelling is required instead. We propose a three-step methodological framework for assessing and correcting DEMs to quantify glacier elevation changes: (i) remove DEM shifts, (ii) check for elevation-dependent biases, and (iii) check for higher-order, sensor-specific biases. A simple, analytic and robust method to co-register elevation data is presented in regions where stable terrain is either plentiful (case study New Zealand) or limited (case study Svalbard). The method is demonstrated using the three global elevation data sets available to date, SRTM, ICESat and the ASTER GDEM, and with automatically generated DEMs from satellite stereo instruments of ASTER and SPOT5-HRS. After 3-D co-registration, significant biases related to elevation were found in some of the stereoscopic DEMs. Biases related to the satellite acquisition geometry (along/cross track) were detected at two frequencies in the automatically generated ASTER DEMs. The higher frequency bias seems to be related to satellite jitter, most apparent in the backlooking pass of the satellite. The origins of the more significant lower frequency bias is uncertain. ICESat-derived elevations are found to be the most consistent globally available elevation data set available so far. Before performing regional-scale glacier elevation change studies or mosaicking DEMs from multiple individual tiles (e.g. ASTER GDEM), we recommend to co-register all elevation data to ICESat as a global vertical reference system.