Abstract
Himalayan glaciers are undergoing rapid mass loss but rates of contemporary change lack long-term (centennial-scale) context. Here, we reconstruct the extent and surfaces of 14,798 Himalayan glaciers during the Little Ice Age (LIA), 400 to 700 years ago. We show that they have lost at least 40 % of their LIA area and between 390 and 586 km3 of ice; 0.92 to 1.38 mm Sea Level Equivalent. The long-term rate of ice mass loss since the LIA has been between − 0.011 and − 0.020 m w.e./year, which is an order of magnitude lower than contemporary rates reported in the literature. Rates of mass loss depend on monsoon influence and orographic effects, with the fastest losses measured in East Nepal and in Bhutan north of the main divide. Locally, rates of loss were enhanced with the presence of surface debris cover (by 2 times vs clean-ice) and/or a proglacial lake (by 2.5 times vs land-terminating). The ten-fold acceleration in ice loss we have observed across the Himalaya far exceeds any centennial-scale rates of change that have been recorded elsewhere in the world.
Similar content being viewed by others
Introduction
Meltwater released by Himalayan glaciers forms the headwaters of the major river systems in Asia, supporting food and energy production downstream, as well as maintaining a range of ecosystems and ecosystem services1,2,3. Ongoing recession and thinning of Himalayan glaciers raises concerns about the sustainability of water supply in the region1,3,4,5,6. Recent rates of ice loss across the region exhibit marked spatial and temporal variability4,5 but there is consensus that recession has accelerated over recent decades7,8. However, there has been no consideration of whether the magnitude and rates of glacier mass loss presently reported across the Himalaya are unusual in a longer-term (centennial-scale) context.
Long-term fluctuations of mountain glaciers are known to be controlled primarily by regional climate and topography9,10. Variations in the South Asian monsoon affect the southerly High Mountain Asia (HMA) mountain ranges and the mid-latitude westerlies are particularly relevant to understanding the glaciology of the western mountain ranges11,12 and both these two atmospheric systems interact with the Tibetan anticyclone, but the exact relationships that describe glacier response to changes in precipitation and temperature are poorly constrained13. Local factors such as the presence or absence of a debris layer or proglacial lake can play an important role14,15,16, although the overall impact of surface debris on ablation remains unclear17,18. Internally-consistent methods that robustly quantify changes exhibited by hundreds to thousands of glaciers at a time are required if such knowledge gaps are to be filled7,8,19.
In this study, we made the first assessment of Himalayan glacier changes over a centennial time scale. Specifically, this assessment was achieved using a combination of geomorphological mapping and palaeo-ice surface reconstruction techniques. We inferred the former extent (Fig. 1A,B) and surface elevation (Fig. 1C) of 14,798 glaciers during the Little Ice Age (LIA) maximum. The LIA was a period of pronounced climate cooling that culminated between 400 and 700 years ago across the Himalaya13. It represents the last period of widespread glacier expansion in the Himalaya and is therefore the benchmark position from which modern glaciers are currently receding. Our regional-scale analysis spanned ~ 2300 km of the Himalaya, across 20 degrees of longitude from Spiti Lahaul, India in the west, to Bhutan in the east and is unprecedented in its spatio-temporal coverage and density. Mapped glacier extents were used to quantify spatially-distributed elevation change (Fig. 1D) and hence glacier-specific volume change and mass balance from the LIA to present.
Glacier area, volume and mass changes
Himalayan glaciers (n = 14,798) covered ~ 28,000 km2 at the LIA maximum. The contemporary area of these glaciers (n = 19,484) according to the Randolph Glacier Inventory20 is ~ 19,600 km2. LIA glacier coverage across the Himalaya was therefore at least 40% greater than at present. This is a conservative area loss estimate given our mapping protocol; see Supplementary Information (SI). Areal losses since the LIA were proportionally consistent between the major regions of Spiti Lahaul, West Nepal, East Nepal and Bhutan (Fig. 2).
Total Himalayan glacier volume loss since the LIA maximum was between 390 km3 and 586 km3. Relatively small losses can be seen around the terminus areas of most glaciers, with maximum surface lowering occurring several kilometers up-glacier, particularly on large valley-floor glacier tongues (1.7% of total glacier area) where lowering commonly exceeds 100 m (Fig. 1D). Rates of surface lowering ranged from − 0.03 m/year to − 0.18 m/year (Fig. 2A) and the greatest rates of change were in the eastern Himalaya (Fig. 2A).
Regional contributions to absolute volume loss were disproportionate given that the total number and density of glaciers are both highly spatially variable. For example, at the LIA maximum, Spiti Lahaul had 55% of the total number of glaciers and 40% of the glacierised area but contributed only 30% to the total measured volume change across the broader region. In contrast, East Nepal, with just 16% of the total number of glaciers at the LIA maximum, and 25% of the glacierised area, contributed 35% of the total ice volume loss. In comparison, West Nepal and Bhutan had 20% and 9% of glaciers at the LIA, and contributed 22% and 13% of volume loss, respectively.
Estimates of long-term mass balance do not differ markedly between sub-regions, with East Nepal and Bhutan showing the most negative median rates overall, at − 0.013 to − 0.023 m w.e./year and − 0.012 to − 0.022 m w.e./year, respectively (Fig. 2B). The same two regions also showed the greatest within-region heterogeneity, with distinctly wider inter-quartile ranges in mass balance when compared to the two regions of Spiti Lahaul and West Nepal further west (Fig. 2B).
Influence of debris cover and terminus type
To understand spatially-distributed differences in glacial mass balance rates since the LIA, glaciers were assigned to one of four mutually exclusive groups depending on their (contemporary) surface characteristics (debris-covered vs clean-ice) and their (contemporary) terminus environment (lake vs land). Combined groups were very heterogenous in number of glaciers: debris-lake (1%), debris-land (7%), clean-lake (4%) and clean-land (88%).
Using generalised linear models (GLM), we identified significant differences in mass balance of all glaciers north versus all glaciers south of the main divide, between regions, and between glacier-terminus types (Fig. 3). See SI for summary statistics and median rates of mass loss for each group. The difference in glacier changes between north and south (Supplementary Fig. SI 2) most likely reflects large-scale atmosphere–topography interactions, which serve to bring orographically-enhanced precipitation to the southern slopes for example, as well as more local effects such as shading, particularly for those glaciers flowing to the north. East Nepal experienced the most negative mass balances for all terminus types (Fig. 3). In all regions, clean-lake, debris-land, and debris-lake glaciers had more negative mass balances compared to clean-land glaciers (Fig. 3; Supplementary Fig. SI 2). Clean-lake glaciers account for only 4% of the total number of glaciers in the Himalaya but contributed a disproportionately large amount to the total volume loss of around 14% (Supplementary Fig. SI 2).
Lake-terminating glaciers (clean-lake and debris-lake) had a median rate of mass loss that was two and a half times more negative (Fig. 3; p-value = < 0.01; Supplementary Table SI2) than that of clean-land or debris-land glaciers. Examining these lake-terminating glaciers south and north of the main divide shows that they lost 40.6 km3 ice and 69.4 km3 ice, representing 13% and 40% of total ice loss on each side of the divide, respectively (Supplementary Fig. SI 2).
The presence/absence of a debris-cover exerts a strong control on rates of mass loss that is largely identical north and south of the topographic divide; debris-land and debris-lake glaciers combined contributed 47% (south) and 46% (north) to total volume loss (Supplementary Fig. SI 2), despite comprising only 6% and 9% of the number of glaciers, respectively. Separately, it is debris-covered glaciers without a lake that make the greatest contribution to mass loss because of their larger number, accounting for 38% of the loss overall. The median mass balance of all debris-land glaciers has been two times more negative than that of clean-land glaciers (Fig. 3; p-value = < 0.01; Supplementary Table SI2).
Comparing the median rates of mass loss between debris-land and debris-lake we find that lakes have enhanced mass loss by − 0.01 m w.e./year. Similarly, comparing the median of the clean-lake group with that of the debris-lake group suggests an enhancement in mass loss by debris cover of − 0.07 m w.e./year. The additive effect of these two glaciological attributes becomes apparent when reminded that clean-ice glaciers had a median mass loss of 0.01 m w.e./year and that debris-covered glaciers that terminate in lakes have experienced a median rate of mass loss of − 0.027 m w.e./year, which is the most negative glacier-specific mass balances of any group of glaciers in all regions (Fig. 3).
Controls on volume loss and mass balance
Our comprehensive dataset of LIA glacier extents allows us to make the first empirical investigation of the controls on long-term (centennial-scale) Himalayan glacier change. It is well known from quantification of recent (decadal) glacier changes during the satellite era that glaciers towards the east of the Himalaya are particularly sensitive to changes in air temperature21, amplified by variations in the seasonal precipitation regime22,23. These studies also show rates of mass loss to be largely uniform across Spiti Lahaul and East and West Nepal, with slightly elevated rates of mass loss recorded in Bhutan. The east–west pattern likely reflects spatially-variable temperature changes during the late twentieth/early twenty-first century24 and a recent weakening of the monsoon25. Our data, on the other hand, capture long-term and aggregated glacier changes, which are likely affected by changing temperature and precipitation patterns across different sub-regions and the subsequent timing of deglaciation26. The longitudinal gradient (p = 0.03; r2 = 0.06) in our results of both surface lowering and mass balance (Supplementary Table SI3), which shows a clear increase in mass loss in eastern regions compared to those in the west (Fig. 2), is much stronger than suggested by studies focusing on decadal-scale mass-balance and surface dynamical changes5,27.
Local-scale variations in terminus condition and surface cover are superimposed on these regional trends and play a crucial role in determining rates of change (Fig. 3). Debris-cover is a key factor in this local variability, with our results suggesting their long-term rate of mass loss outstrips that of their clean-ice counterparts by a factor of 3. This is perhaps surprising given that debris covers exceeding a critical thickness (~ 0.02 m) tend to reduce melt when compared to clean-ice glaciers at the same elevation28,29. However, several studies of contemporary mass loss have found comparable thinning rates over clean and debris-covered glaciers4,30. The physical processes underlying this anomaly are thought to include locally-enhanced ablation owing to ice cliff and pond development18,31,32 and decreased ice fluxes from the accumulation area, particularly where tongues have become dynamically detatched33. Over centennial timescales, a more plausible explanation is that these glaciers have been sustained at lower elevations than would otherwise be possible and therefore have large areas of ice available for ablation due to climatic forcing14,34,35.
Lakes at glacier termini exert a number of local thermo-mechanical effects causing ice surface drawdown and glacier acceleration, and hence enhanced ice mass loss7,36,37. Glaciers developing lakes at their termini therefore have an evolution of mass balance that is non-linear over decades to centuries15,38. Although glaciers with lakes are low in number (relative to those without) across the Himalaya, our results show, unequivocally, that they contribute a disproportionate amount to overall ice loss. It is likely therefore that lake development will exert an increasingly strong control on future Himalayan glacier mass balance as they grow in both number and volume39.
Comparison with contemporary rates of change across the Himalaya
Centennial-scale rates of change provide a longer-term context to reported changes of Himalayan glaciers over the last few decades (Supplementary Table SI3). Crucially, our analysis reveals that glacier changes (median − 0.011 to − 0.020 m w.e./year) since the LIA maximum (1300 CE to 1600 CE, respectively) are an order of magnitude lower than reported for recent decades; for example Gardelle et al.4 for all Himalayan sub-regions (mean − 0.35 m w.e./year), Kääb et al.40 for 2003 to 2008 (mean − 0.42 m w.e./year), Brun et al.5 for 2000 to 2016 (mean − 0.36 m w.e./year), Maurer et al.8 for 1975 to 2000 (mean − 0.21 m w.e./year) and for 2000 to 2016 (mean − 0.43 m w.e./year) and Shean et al.41 for 2000 to 2018 (mean − 0.19 m w.e./year), albeit with some slight variations in the exact sub-regions studied. This comparison demonstrates that the mass balance of glaciers across the Himalaya has become dramatically more negative in recent decades in response to climatic forcing.
Comparisons to other world regions
Studies from several other parts of the world have shown that recent decadal rates of glacier mass loss are greater than centennial-scale average losses since LIA, but not on the same scale as those detected for the Himalaya. For example, compared to the LIA-present average rate, there has been a doubling of mass loss for the Vatnajökull ice cap in Iceland between 2002 and 201042, a 23% acceleration of mass losses in NE Greenland between 1980 and 201443, a doubling of mass loss since 1986 for Patagonia44,45, and a doubling of losses for the Southern Alps of New Zealand between 2009 to 201946. The order of magnitude increase in Himalayan glacier mass loss in recent decades compared to the post-LIA average therefore represents the most dramatic glacier response of any world region.
In conclusion, this study provides the first self-consistent measurement of Himalayan glacier extent and surface elevation during the LIA maximum across the entire Himalaya. The area loss of Himalayan glaciers since the LIA has been at least 40% of the LIA extent. Ice volume loss since the LIA has been between 390 and 586 km3, which equates to between 332 and 498 Gt mass equivalent and to between 0.92 and 1.38 mm Sea Level Equivalent (SLE). For the first time, we have identified that the centennial rate of mass loss since the LIA is an order of magnitude lower than the rates reported for the last decades. We suggest that nested scales of effects influence these rates of mass loss over a centennial time scale, namely (i) regional climate via the differences north versus south of the main divide and between morpho-climatic regions, and (ii) local topography and glacier morphology via terminus environment and debris-cover. The greatest rates of mass loss are in the eastern Himalaya. Whilst clean-land glaciers produced the most mass loss overall due to their number and size, glaciers with debris-cover had more negative mass balances and those with lakes the most negative mass balances.
Our study highlights three important points. Firstly, by comparison to other world regions the magnitude of the acceleration in glacier mass loss across the Himalayan region is exceptional. Secondly, whilst centennial-scale changes in climate forcings alter spatial gradients in regional glacier mass balance, it is local topography and glaciological attributes such as debris cover and the presence of a proglacial lake that promote rapid mass loss7,14,47. Thirdly, those attributes need to be accounted for in numerical models15,38 because they important drivers of mass loss. Overall, quantification of past glacier variability over a centennial timescale should help validate climate-glacier numerical models and thus produce more reliable projections of future mass loss under a warming climate48.
Methods
Little Ice Age moraine mapping
LIA glacier extents were mapped by extending the RGI 6.0 inventory outlines20 down-valley to the crests of moraines interpreted to represent the LIA maximum. Specifically, we delimited the innermost prominent moraine ridge down-valley of the contemporary glacier termini, where the ridge was often sharp-crested and with unstable, steep sides devoid of vegetation. These innermost ridges have in some valleys been dated to the LIA and contrast with other parts of multi-ridged moraine complexes that are more subdued, vegetated and more likely represent late Holocene ice advances13. Furthermore, choice of the innermost moraines maintains consistency, permitting replicability, and provides a minimum estimate of LIA glacier extent and hence a minimum LIA volume and mass.
Ice surface reconstruction
For each of our LIA outlines, ablation zones were delineated by calculating areas below the equilibrium line altitude (ELA). Glacier specific ELAs were defined using the Area-Altitude Balance Ratio (AABR) method, automated using code developed by Pellitero et al.49. We specified a BR (balance ratio) of 1.75 as suggested by Rea50 for Himalayan glaciers. The area below each ELA was then extracted in an automated fashion, as described by Carrivick et al.43. The vertices of each ablation area were converted to points to enable extraction of elevations from a digital elevation model (DEM); i.e. of terminal and lateral moraine crests and of the ELA. A surface was interpolated between those points to represent the LIA glacier surface. Calculating the difference between our LIA surface and the contemporary DEM indicated the surface lowering that has occurred since the LIA to the date of the DEM. We then converted that to a glacier-specific volume change and a mean annual rate of change. The effect of DEM source (resolution or timing; Supplementary Fig. SI 1) on our calculations of volume was negligible; either using the ALOS 30 m DEM (www.eorc.jaxa.jp/)51 or the HMA 8 m DEM (www.nsidc.org)52. The effect of our choice of interpolation method on volume calculations was also small (Supplementary Table SI1).
Rates of change and uncertainty
We report non-parametric statistics because our mass loss data are non-normally distributed, as they pertain only to glacier ablation areas. Rather than volume uncertainty, by far the biggest effect on our rates of change is the date ascribed to the LIA maximum. Rowan’s13 analysis shows that the LIA across the Himalaya was between 1300 and 1600 CE. Therefore, recognizing that not all Himalayan glaciers will have reached their maximum extent at the same time, we provide reasonable lower and upper bounds by computing rates of glacier change for both 400- and 700-year time periods. Between the time of the LIA maximum and the present, we acknowledge that rates of glacier mass loss may have varied. Episodes of stagnation and minor advances may have occurred, but these cannot be accounted for in this study. The SI details our consideration of the timing of the LIA maximum across the Himalaya, while also providing a description of our datasets, spatial analysis methods and our uncertainty assessments, and a summary of our statistical assessment of differences between groups (north/south of main divide, regions, terminus type; Supplementary Tables SI2 and SI 3).
Data availability
Our study utilized publically available datasets as cited and referenced. Our study produced new datasets; glacier outlines for the LIA and an ice surface reconstruction for the LIA, which are available as shapefiles and geotiffs, respectively, from https://doi.org/10.5518/939.
References
Immerzeel, W. W., van Beek, L. P. H. & Bierkens, M. F. P. Climate change will affect the Asian water towers. Science 328(5984), 1382–1385 (2010).
Immerzeel, W. W. et al. Importance and vulnerability of the world’s water towers. Nature 577(7790), 364–369 (2020).
Bolch, T. et al. The state and fate of Himalayan glaciers. Science 336(6079), 310–314 (2012).
Gardelle, J., Berthier, E., Arnaud, Y. & Kääb, A. Region-wide glacier mass balances over the Pamir-Karakoram-Himalaya during 1999–2011. Cryosphere 7(4), 1263–1286 (2013).
Brun, F., Berthier, E., Wagnon, P., Kääb, A. & Treichler, D. A spatially resolved estimate of High Mountain Asia glacier mass balances from 2000 to 2016. Nat. Geosci. 10(9), 668–673 (2017).
Pritchard, H. D. Asia’s glaciers are a regionally important buffer against drought. Nature 545(7653), 169–174 (2017).
King, O., Bhattacharya, A., Bhambri, R. & Bolch, T. Glacial lakes exacerbate Himalayan glacier mass loss. Sci. Rep. 9(1), 18145 (2019).
Maurer, J. M., Schaefer, J. M., Rupper, S. & Corley, A. Acceleration of ice loss across the Himalayas over the past 40 years. Sci. Adv. 5(6), eaav7266 (2019).
Owen, L. A. Latest Pleistocene and Holocene glacier fluctuations in the Himalaya and Tibet. Quatern. Sci. Rev. 28(21), 2150–2164 (2009).
McGrath, D., Sass, L., O’Neel, S., Arendt, A. & Kienholz, C. Hypsometric control on glacier mass balance sensitivity in Alaska and northwest Canada. Earth’s Future 5(3), 324–336 (2017).
Benn, D. I. & Owen, L. A. The role of the Indian summer monsoon and the mid-latitude westerlies in Himalayan glaciation: Review and speculative discussion. J. Geol. Soc. 155(2), 353–363 (1998).
Mölg, T., Maussion, F. & Scherer, D. Mid-latitude westerlies as a driver of glacier variability in monsoonal High Asia. Nat. Clim. Change. 4(1), 68–73 (2014).
Rowan, A. V. The ‘Little Ice Age’ in the Himalaya: A review of glacier advance driven by Northern Hemisphere temperature change. Holocene 27(2), 292–308 (2017).
Rowan, A. V., Egholm, D. L., Quincey, D. J. & Glasser, N. F. Modelling the feedbacks between mass balance, ice flow and debris transport to predict the response to climate change of debris-covered glaciers in the Himalaya. Earth Planet. Sci. Lett. 430, 427–438 (2015).
Sutherland, J. L. et al. Proglacial lakes control glacier geometry and behavior during recession. Geophys. Res. Lett. 47(19), e2020GL088865 (2020).
King, O., Dehecq, A., Quincey, D. & Carrivick, J. Contrasting geometric and dynamic evolution of lake and land-terminating glaciers in the central Himalaya. Glob. Planet. Change 167, 46–60 (2018).
Pellicciotti, F. et al. Mass-balance changes of the debris-covered glaciers in the Langtang Himal, Nepal, from 1974 to 1999. J. Glaciol. 61(226), 373–386 (2015).
Brun, F. et al. Ice cliff contribution to the tongue-wide ablation of Changri Nup Glacier, Nepal, central Himalaya. Cryosphere 12(11), 3439–3457 (2018).
Farinotti, D. et al. Substantial glacier mass loss in the Tien Shan over the past 50 years. Nat. Geosci. 8(9), 716–722 (2015).
RGI Consortium. Randolph Glacier Inventory—A Dataset of Global Glacier Outlines: Version 6.0: Technical Report, Global Land Ice Measurements from Space. Colorado, USA (2017).
Peng, X. et al. Late Holocene glacier fluctuations in the Bhutanese Himalaya. Glob. Planet. Change 187, 103137 (2020).
Oerlemans, J. & Fortuin, J. P. F. Sensitivity of glaciers and small ice caps to greenhouse warming. Science 258(5079), 115–117 (1992).
Fujita, K. Influence of precipitation seasonality on glacier mass balance and its sensitivity to climate change. Ann. Glaciol. 48, 88–92 (2008).
Krusic, P. J. et al. Six hundred thirty-eight years of summer temperature variability over the Bhutanese Himalaya. Geophys. Res. Lett. 42(8), 2988–2994 (2015).
Salerno, F. et al. Weak precipitation, warm winters and springs impact glaciers of south slopes of Mt. Everest (central Himalaya) in the last 2 decades (1994–2013). Cryosphere 9(3), 1229–1247 (2015).
Shekhar, M. et al. Himalayan glaciers experienced significant mass loss during later phases of little ice age. Sci. Rep. 7(1), 10305 (2017).
Dehecq, A. et al. Twenty-first century glacier slowdown driven by mass loss in High Mountain Asia. Nat. Geosci. 12(1), 22–27 (2019).
Østrem, G. Ice melting under a thin layer of moraine, and the existence of ice cores in moraine ridges. Geogr. Ann. 41(4), 228–230 (1959).
Kayastha, R. B., Takeuchi, Y., Nakawo, M. & Ageta, Y. Practical prediction of ice melting beneath various thickness of debris cover on Khumbu Glacier, Nepal, using a positive degree-day factor. Int. Assoc. Hydrol. Sci. Publ. 264, 71–81 (2000).
Kääb, A., Berthier, E., Nuth, C., Gardelle, J. & Arnaud, Y. Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas. Nature 488(7412), 495–498 (2012).
Lovell, A. M., Carr, J. R. & Stokes, C. R. Spatially variable glacier changes in the annapurna conservation area, Nepal, 2000 to 2016. Remote Sens. 11(12), 1452 (2019).
Watson, C. S., Quincey, D. J., Carrivick, J. L. & Smith, M. W. Ice cliff dynamics in the Everest region of the Central Himalaya. Geomorphology 278, 238–251 (2017).
Rowan, A. V. et al. The role of differential ablation and dynamic detachment in driving accelerating mass loss from a debris-covered Himalayan glacier. J. Geophys. Res. Earth Surf. 126, e2020JF005761 (2021).
Scherler, D., Bookhagen, B. & Strecker, M. R. Spatially variable response of Himalayan glaciers to climate change affected by debris cover. Nat. Geosci. 4(3), 156–159 (2011).
Anderson, L. S. & Anderson, R. S. Modeling debris-covered glaciers: Response to steady debris deposition. Cryosphere 10(3), 1105–1124 (2016).
Tsutaki, S. et al. Contrasting thinning patterns between lake- and land-terminating glaciers in the Bhutanese Himalaya. Cryosphere 13(10), 2733–2750 (2019).
Carrivick, J. L. & Tweed, F. S. A global assessment of the societal impacts of glacier outburst floods. Glob. Planet. Change 144, 1–16 (2016).
Carrivick, J. L., Tweed, F. S., Sutherland, J. L. & Mallalieu, J. Toward numerical modeling of interactions between ice-marginal proglacial lakes and glaciers. Front. Earth Sci. 8, 500 (2020).
Linsbauer, A. et al. Modelling glacier-bed overdeepenings and possible future lakes for the glaciers in the Himalaya—Karakoram region. Ann. Glaciol. 51(71), 119–130 (2016).
Kääb, A., Treichler, D., Nuth, C. & Berthier, E. Brief Communication: Contending estimates of 2003–2008 glacier mass balance over the Pamir–Karakoram–Himalaya. Cryosphere 9(2), 557–564 (2015).
Shean, D. E. et al. A systematic, regional assessment of high mountain asia glacier mass balance. Front. Earth Sci. 7, 363 (2020).
Hannesdóttir, H., Björnsson, H., Pálsson, F., Aðalgeirsdóttir, G. & Guðmundsson, S. Changes in the southeast Vatnajökull ice cap, Iceland, between ~ 1890 and 2010. Cryosphere 9(2), 565–585 (2015).
Carrivick, J. L. et al. Accelerated volume loss in glacier ablation zones of NE Greenland, Little Ice Age to present. Geophys. Res. Lett. 46(3), 1476–1484 (2019).
Glasser, N. F., Harrison, S., Jansson, K. N., Anderson, K. & Cowley, A. Global sea-level contribution from the Patagonian Icefields since the Little Ice Age maximum. Nat. Geosci. 4(5), 303–307 (2011).
Davies, B. J. & Glasser, N. F. Accelerating shrinkage of Patagonian glaciers from the Little Ice Age (~AD 1870) to 2011. J. Glaciol. 58(212), 1063–1084 (2012).
Carrivick, J. L., James, W. H. M., Grimes, M., Sutherland, J. L. & Lorrey, A. M. Ice thickness and volume changes across the Southern Alps, New Zealand, from the little ice age to present. Sci. Rep. 10(1), 13392 (2020).
Brun, F. et al. Heterogeneous influence of glacier morphology on the mass balance variability in high mountain Asia. J. Geophys. Res. Earth Surf. 124(6), 1331–1345 (2019).
Khan, S. A. et al. Centennial response of Greenland’s three largest outlet glaciers. Nat. Commun. 11(1), 5718 (2020).
Pellitero, R. et al. A GIS tool for automatic calculation of glacier equilibrium-line altitudes. Comput. Geosci. 82, 55–62 (2015).
Rea, B. R. Defining modern day Area-Altitude Balance Ratios (AABRs) and their use in glacier-climate reconstructions. Quatern. Sci. Rev. 28(3–4), 237–248 (2009).
Tadono, T. et al. Precise global DEM generation by ALOS PRISM. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. II–4, 71–76 (2014).
Shean, D. High Mountain Asia 8-meter DEM Mosaics Derived from Optical Imagery, Version 1. (NASA National Snow and Ice Data Center Distributed Active Archive Center, 2017).
Acknowledgements
W. James was in receipt of NERC studentship WHMJ, NE/K500847/1 for a part of this work. A. Rowan shared the Little Ice Age moraine date database. O. King helped acquire and manage the ALOS DEM. F. Brun shared shapefiles of the morpho-climatic regions.
Author information
Authors and Affiliations
Contributions
This study was conceptualized by D.J.Q. and S.J.C and J.L.C. Glacier LIA extents were mapped by E.L. across the eastern Himalaya, by E.L. and S.J. across the central Himalaya, and by J.L.C. across the western Himalaya. D.J.Q. and S.J.C. checked the mapping. J.L.C. developed the analysis workflow and led the data analysis with assistance from E.L. for glacier attributes, from L.E.B. for the statistical analysis and from W.H.M.J. on derivation of glacier ablation areas. J.L.C. prepared all the figures. All authors interpreted the results. J.L.C. led the writing of the paper and all authors contributed to writing.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lee, E., Carrivick, J.L., Quincey, D.J. et al. Accelerated mass loss of Himalayan glaciers since the Little Ice Age. Sci Rep 11, 24284 (2021). https://doi.org/10.1038/s41598-021-03805-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-03805-8
This article is cited by
-
Characteristics and changes of glacial lakes and outburst floods
Nature Reviews Earth & Environment (2024)
-
Decadal evolution of the glacial lakes: a comprehensive study on spatial mapping and dynamics (2000–2020) in Baspa River Basin, India
Proceedings of the Indian National Science Academy (2024)
-
An ice-snow avalanche triggered small glacial lake outburst flood in Birendra Lake, Nepal Himalaya
Natural Hazards (2024)
-
The effects of climatic change and inter-annual variability on glacier retreat from ~ 1850s AD moraines in the Kuoqionggangri peak region, southern Tibetan Plateau
Climate Dynamics (2024)
-
Future emergence of new ecosystems caused by glacial retreat
Nature (2023)