Despite their extreme elevation, glaciers on the Tibetan Plateau are losing mass in response to atmospheric warming, the pattern of which purportedly reflects regional contrasts in climate. Here we examine the evolution of glaciers along ~500 km of the Tanggula Shan, Central-Eastern Tibetan Plateau. Using remotely sensed datasets, we quantified changes in glacier mass, area and surface velocity, and compared these results to time series of meteorological observations, in order to disentangle drivers of glacier mass loss since the 1960s. Glacier mass loss has increased (from −0.21 ± 0.12 m w.e. a−1 in 1960s–2000s, to −0.52 ± 0.18 m w.e. a−1 in 2000s–2015/18) in association with pervasive positive temperature anomalies (up to 1.85°C), which are pronounced at the end of the now lengthened ablation season. However, glacier mass budget perturbations do not mirror the magnitude of temperature anomalies in sub-regions, thus additional factors have heightened glacier recession. We show how proglacial lake expansion and glacier surging have compounded glacier recession over decadal/multi-decadal time periods, and exert similar influence on glacier mass budgets as temperature changes. Our results demonstrate the importance of ice loss mechanisms not often incorporated into broad-scale glacier projections, which need to be better considered to refine future glacier runoff estimates.

Himalayan glaciers have been shrinking and losing mass rapidly since 1970s with an enhanced rate after 2000. The shrinkage is, however, quite heterogeneous and it is important to document individual glacier characteristics and their changes at the basin scale. We present an updated glacier inventory of the Upper Alaknanda Basin (UAB), Central Himalaya for the year 2020 and report area, debris cover and length changes for the periods 1994–2006 and 2006–2020 based on remote-sensing data. We identified 198 glaciers, comprising an area of 354.6 ± 8.5 km2, and classified them according to their size and morphology. The glaciers of the basin lost 4.2 ± 2.9% (0.16 ± 0.11% a−1) of their frontal area (from 368.6 ± 9.2 to 353.0 ± 5.3 km2) from 1994 to 2020. The average retreat rate was higher in the period 2006–2020 (13.3 ± 1.8 m a−1) in comparison to 1994–2006 (9.3 ± 1.9 m a−1). However, the area change rate was similar for the two periods (0.14 ± 0.27% a−1 for 1994–2006 and 0.16 ± 0.19% a−1 for 2006–2020). The debris-covered area has increased by 13.4 ± 4.4% from 1994 to 2020. A comparison with previous studies in UAB indicates consistent area loss of ~0.15% a−1 since the 1960s.

Despite previous studies, glacier–lake interactions and future lake development in the Poiqu River basin, central Himalaya, are still not well understood. We mapped glacial lakes, glaciers, their frontal positions and ice flow from optical remote sensing data, and calculated glacier surface elevation change from digital terrain models. During 1964–2017, the total glacial-lake area increased by ~110%. Glaciers retreated with an average rate of ~1.4 km2 a−1 between 1975 and 2015. Based on rapid area expansion (>150%), and information from previous studies, eight lakes were considered to be potentially dangerous glacial lakes. Corresponding lake-terminating glaciers showed an overall retreat of 6.0 ± 1.4 to 26.6 ± 1.1 m a−1 and accompanying lake expansion. The regional mean glacier elevation change was −0.39 ± 0.13 m a−1 while the glaciers associated with the eight potentially dangerous lakes lowered by −0.71 ± 0.05 m a−1 from 1974 to 2017. The mean ice flow speed of these glaciers was ~10 m a−1 from 2013 to 2017; about double the mean for the entire study area. Analysis of these data along with climate observations suggests that ice melting and calving processes play the dominant role in driving lake enlargement. Modelling of future lake development shows where new lakes might emerge and existing lakes could expand with projected glacial recession.

Information about the ice volume stored in glaciers is of high importance for sustainable water management in many arid regions of Central Asia. Several methods to estimate the ice volume exist. However, none of them take the specific characteristics of flat terminus debris-covered glaciers into account. We present a method for deriving spatially-distributed ice thickness for debris-covered dendritic glaciers, which are common not only in Central Tien Shan but also in several other mountain ranges in High Asia. The method relies on automatically generated branch lines, observed surface velocities and surface topographic parameters as basic input. Branch lines were generated using Thiessen polygons and Dijkstra's path algorithm. Ice thicknesses for four debris-covered glaciers – South Inylchek, Kaindy, Tomur and Koxkar glaciers – have been estimated along the branch line network solving the equation of laminar flow. For Koxkar and South Inylchek glaciers, respectively, maximum thicknesses of ~250 and 380 m were estimated. These results differ by ~50 m compared with GPR measurements with an uncertainty for the debris-covered parts of ~40%. Based on geodetic mass balances, we estimate that the investigated glaciers lost between 6 and 28% of their volume from 1975 to the early 2000s.

Despite renewed efforts to better understand glacier change and recognize glacier change trends in the Andes, relatively large areas in the Andes of Argentina and Chile are still not investigated. In this study, we report on glacier elevation and mass changes in the outer region of the Northern and Southern Patagonian Icefields in the Southern Patagonian Andes. A newly-compiled Landsat ETM+ derived glacier inventory (consisting of 2253 glaciers and ~1314 ± 66 km2 of ice area) and differencing of the SRTM and SPOT5 DEMs were used to derive glacier-specific elevation changes over the 2000–12 period. The investigated glaciers showed a volume change of −0.71 ± 0.55 km3 a−1, yielding a surface lowering of 0.52 ± 0.35 m a−1 on average and an overall mass loss of 0.46 ± 0.37 m w.e. a−1. Highly variable individual glacier responses were observed and interestingly, they were less negative than previously reported for the neighboring Patagonian Icefields.

Here, we present a comprehensive assessment of Siachen Glacier (East Karakoram), in terms of its area and elevation change, velocity variations and mass budget, utilizing different satellite datasets including Landsat, Hexagon, Cartosat-I, Shuttle Radar Topography Mission, Envisat Advanced Synthetic Aperture Radar and Japanese Advanced Land Observing Satellite Phased Array-type L-band SAR. The total areal extent of Siachen Glacier did not change significantly between 1980 and 2014; however the exposed-ice area decreased during that period. The terminus of the glacier has experienced substantial downwasting (on average 30 m) over the period of 1999–2007, followed by a retreat of the transition between exposed and debris-covered ice by a distance of 1.3 km during the short span 2007–14. The spatial patterns of the elevation difference and velocity are heterogeneous over the large areal extent of Siachen Glacier. The average velocity of the entire glacier, as computed between 11 December 2008 and 26 January 2009, was 12.3 ± 0.4 cm d−1, while those estimated separately for the accumulation and ablation regions were 9.7 ± 0.4 cm d−1 and 20.4 ± 0.4 cm d−1, respectively. The mass budget of Siachen Glacier is estimated to be –0.03 ± 0.21 m w.e. a−1 for the period of 1999–2007.

Thinning rates for the debris-covered Gangotri Glacier and its tributary glaciers during the period 1968–2014, length variation and area vacated at the snout from 1965 to 2015, and seasonal variation of ice-surface velocity for the last two decades have been investigated in this study. It was found that the mass loss of Gangotri and its tributary glaciers was slightly less than those reported for other debris-covered glaciers in the Himalayan regions. The average velocity during 2006–14 decreased by ~6.7% as compared with that during 1993–2006. The debris-covered area of the main trunk of Gangotri Glacier increased significantly from 1965 until 2015 with the maximum rate of increase (0.8 ± 0.2 km2 a−1) during 2006–15. The retreat (~9.0 ± 3.5 m a−1) was less in recent years (2006–2015) but the down-wasting (0.34 ± 0.2 m a−1) in the same period (2006–2014) was higher than that (0.20 ± 0.1 m a−1) during 1968–2006. The study reinforced the established fact that the glacier length change is a delayed response to climate change and, in addition, is affected by debris cover, whereas glacier mass balance is a more direct and immediate response. Therefore, it is recommended to study the glacier mass balance and not only the glacier extent, to conclude about a glacier's response to climate change.

This study explores the link between area increase of Imja Tsho (Lake) and changes of Imja Glacier (area ~25km2) under the influence of climate change using multitemporal satellite imagery and local climate data. Between 1962 and 2013, Imja Lake expanded from 0.03±0.01 to 1.35±0.05 km2 at a rate of 0.026±0.001 km2 a-1. The mean glacier-wide flow velocity was 37±30ma-1 during 1992–93 and 23±15ma-1 during 2013–14, indicating a decreasing velocity. A mean elevation change of –1.29±0.71ma-1 was observed over the lower part of the glacier in the period 2001–14, with a rate of –1.06±0.63ma-1 in 2001–08 and –1.56±0.80ma-1 in 2008–14. We conclude that the decrease in flow velocity is mainly associated with reduced accumulation due to a decrease in precipitation during the last few decades. Furthermore, glacier ablation has increased due to increasing maximum temperatures during the post-monsoon months. Decreased glacier flow velocities and increased mass losses induce the formation and subsequent expansion of glacial lakes under favourable topographic conditions.

Thick debris cover on glaciers can significantly reduce ice melt. However, several studies have suggested that debris-covered glaciers in the Himalaya might have lost mass at a rate similar to debris-free glaciers. We reconstruct elevation and mass changes for the debris-covered glaciers of the upper Langtang valley, Nepalese Himalaya, using a digital elevation model (DEM) from 1974 stereo Hexagon satellite data and the 2000 SRTM (Shuttle Radar Topography Mission) DEM. Uncertainties are high in the accumulation areas, due to data gaps in the SRTM and difficulties with delineation of the glacier borders. Even with these uncertainties, we obtain thinning rates comparable to those of several other studies in the Himalaya. In particular, we obtain a total mass balance for the investigated debris-covered glaciers of the basin of –0.32 ± 0.18 m w.e. a−1. However, there are major spatial differences both between glaciers and within any single glacier, exhibiting a very distinct nonlinear mass-balance profile with elevation. Through analysis of surface velocities derived from Landsat ETM+ imagery, we show that thinning occurs in areas of low velocity and low slope. These areas are prone to a general, dynamic decay of surface features and to the development of supraglacial lakes and ice cliffs, which may be responsible for a considerable increase in overall glacier ablation.

The Randolph Glacier Inventory (RGI) is a globally complete collection of digital outlines of glaciers, excluding the ice sheets, developed to meet the needs of the Fifth Assessment of the Intergovernmental Panel on Climate Change for estimates of past and future mass balance. The RGI was created with limited resources in a short period. Priority was given to completeness of coverage, but a limited, uniform set of attributes is attached to each of the ~198 000 glaciers in its latest version, 3.2. Satellite imagery from 1999–2010 provided most of the outlines. Their total extent is estimated as 726 800 ± 34 000 km2. The uncertainty, about ±5%, is derived from careful single-glacier and basin-scale uncertainty estimates and comparisons with inventories that were not sources for the RGI. The main contributors to uncertainty are probably misinterpretation of seasonal snow cover and debris cover. These errors appear not to be normally distributed, and quantifying them reliably is an unsolved problem. Combined with digital elevation models, the RGI glacier outlines yield hypsometries that can be combined with atmospheric data or model outputs for analysis of the impacts of climatic change on glaciers. The RGI has already proved its value in the generation of significantly improved aggregate estimates of glacier mass changes and total volume, and thus actual and potential contributions to sea-level rise.

We investigated area changes in glaciers covering an area of ∼200 km2 in the Tista basin, Sikkim, Eastern Indian Himalaya, between ∼1990 and 2010 using Landsat Thematic Mapper (TM) and Indian Remote-sensing Satellite (IRS) images and related the changes to debris cover, supraglacial lakes and moraine-dam lakes. The glaciers lost an area of 3.3 ± 0.8% between 1989/90 and 2010. More detailed analysis revealed an area loss of 2.00 ± 0.82, 2.56 ± 0.61 and 2.28 ± 2.01 km2 for the periods 1989–97, 1997–2004/05 and 2004–2009/10, respectively. This indicates an accelerated retreat of glaciers after 1997. On further analysis, we observed (1) the formation and expansion of supraglacial lakes on many debris-covered glaciers and (2) the merging of these lakes over time, leading to the development of large moraine-dam lakes. We also observed that debris-covered glaciers with lakes lose a greater area than debris-covered glaciers without lakes and debris-free glaciers. The climatic data for 24 years (1987–2011), measured at the Gangtok meteorological station (1812 m a.s.l.), showed that the region experienced a 1.0°C rise in the summer minimum temperature and a 2.0°C rise in the winter minimum temperature, indicating hotter summers and warmer winters. There was no significant trend in the total annual precipitation. We find that glacier retreat is caused mainly by a temperature increase and that debris-covered glaciers can retreat at a faster rate than debris-free glaciers, if associated with lakes.

Glacier outlines are mapped for the upper Bhagirathi and Saraswati/Alaknanda basins of the Garhwal Himalaya using Corona and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite images acquired in 1968 and 2006, respectively. A subset of glaciers was also mapped using Landsat TM images acquired in 1990. Glacier area decreased from 599.9 ± 15.6 km2 (1968) to 572.5 ± 18.0 km2 (2006), a loss of 4.6 ± 2.8%. Glaciers in the Saraswati/Alaknanda basin and upper Bhagirathi basin lost 18.4 ± 9.0 km2 (5.7 ± 2.7%) and 9.0 ± 7.7 km2 (3.3 ± 2.8%), respectively, from 1968 to 2006. Garhwal Himalayan glacier retreat rates are lower than previously reported. More recently (1990–2006), recession rates have increased. The number of glaciers in the study region increased from 82 in 1968 to 88 in 2006 due to fragmentation of glaciers. Smaller glaciers (<1 km2) lost 19.4 ± 2.5% (0.51 ± 0.07% a−1) of their ice, significantly more than for larger glaciers (>50 km2) which lost 2.8 ± 2.7% (0.074 ± 0.071 % a−1). From 1968 to 2006, the debris-covered glacier area increased by 17.8 ± 3.1% (0.46 ± 0.08% a−1) in the Saraswati/Alaknanda basin and 11.8 ± 3.0% (0.31 ± 0.08% a−1) in the upper Bhagirathi basin. Climate records from Mukhim (∼1900 m a.s.l.) and Bhojbasa (∼3780 m a.s.l.) meteorological stations were used to analyze climate conditions and trends, but the data are too limited to make firm conclusions regarding glacier–climate interactions.

Multitemporal space imagery from 1962 (Corona KH-4), 1992 (Landsat TM), 2001 and 2005 (Terra ASTER) was used to investigate the glacier changes in the Khumbu Himal, Nepal. The ice coverage in the investigation area decreased by about 5% between 1962 and 2005, with the highest retreat rates occurring between 1992 and 2001. The debris coverage increased concomitantly with the decrease in total glacier area. The clean-ice area decreased by >10%. Digital terrain model (DTM) generation from the early Corona KH-4 stereo data in this high-relief terrain is time-consuming, and the results still contain some elevation errors. However, these are minor in the snow-free areas with gentle slopes. Thus comparison of the surfaces of the debris-covered glacier tongues based on the Corona DTM and an ASTER DTM is feasible and shows the downwasting of the debris-covered glaciers. The highest downwasting rates, more than 20 m (>0.5 m a−1), can be found near the transition zone between the active and the stagnant glacier parts of the debris-covered glacier tongues. The downwasting is lower, but still evident, in the active ice areas and at the snout with thick debris cover. All investigated debris-covered glaciers in the study area show similar behaviour. The estimated volume loss for the investigated debris-covered glacier tongues is 0.19 km3.