skip to main content


Title: Acceleration of ice loss across the Himalayas over the past 40 years
Himalayan glaciers supply meltwater to densely populated catchments in South Asia, and regional observations of glacier change over multiple decades are needed to understand climate drivers and assess resulting impacts on glacier-fed rivers. Here, we quantify changes in ice thickness during the intervals 1975–2000 and 2000–2016 across the Himalayas, using a set of digital elevation models derived from cold war–era spy satellite film and modern stereo satellite imagery. We observe consistent ice loss along the entire 2000-km transect for both intervals and find a doubling of the average loss rate during 2000–2016 [−0.43 ± 0.14 m w.e. year −1 (meters of water equivalent per year)] compared to 1975–2000 (−0.22 ± 0.13 m w.e. year −1 ). The similar magnitude and acceleration of ice loss across the Himalayas suggests a regionally coherent climate forcing, consistent with atmospheric warming and associated energy fluxes as the dominant drivers of glacier change.  more » « less
Award ID(s):
1853881
NSF-PAR ID:
10145909
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Science Advances
Volume:
5
Issue:
6
ISSN:
2375-2548
Page Range / eLocation ID:
eaav7266
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Southern Andean glaciers contribute substantially to global sea-level rise. Unfortunately, mass balance estimates prior to 2000 are limited, hindering our understanding of the evolution of glacier mass changes over time. Elevation changes over 1976/1979 to 2000 derived from historical KH-9 Hexagon imagery and NASADEM provide the basis for geodetic mass balance estimates for subsets of the Northern Patagonian Icefield (NPI) and the Southern Patagonian Icefield (SPI), extending current mass balance observations by ∼20 years. Geodetic mass balances were −0.63 ± 0.03 m w.e. yr −1 for 63% of the NPI and −0.33 ± 0.05 m w.e. yr −1 for 52% of the SPI glacierized areas for this historical period. We also extend previous estimates temporally by 25% using NASADEM and ASTER elevation trends for the period 2000 to 2020, and find geodetic mass balances of −0.86 ± 0.03 m w.e. yr −1 for 100% of the NPI and −1.23 ± 0.04 m w.e. yr −1 for 97% of the SPI glacierized areas. 2000–2020 aggregations for the same areas represented in the 1976/1979 to 2000 estimates are −0.78 ± 0.03 m w.e. yr −1 in the NPI and −0.80 ± 0.04 m w.e. yr −1 on the SPI. The significant difference in SPI geodetic mass balance in the modern period for 100% vs. 52% of the glacierized area suggests subsampling leads to significant biases in regional mass balance estimates. When we compare the same areas in each time period, the results highlight an acceleration of ice loss by a factor of 1.2 on the NPI and 2.4 on the SPI in the 21st century as compared to the 1976/1979 to 2000 period. While lake-terminating glaciers show the most significant increase in mass loss rate from 1976/1979–2000 to 2000–2020, mass balance trends are highly variable within glaciers of all terminus environments, which suggests that individual glacier sensitivity to climate change is dependent on a multitude of morphological and climatological factors. 
    more » « less
  2. Abstract. Over the last century, northwestern Canada experienced some of the highest rates of tropospheric warming globally, which caused glaciers in the region to rapidly retreat. Our study seeks to extend the record of glacier fluctuations and assess climate drivers prior to the instrumental record in the Mackenzie and Selwyn mountains of northwestern Canada. We collected 27 10Be surface exposure ages across nine cirque and valley glacier moraines to constrain the timing of their emplacement. Cirque and valley glaciers in this region reached their greatest Holocene extents in the latter half of the Little Ice Age (1600–1850 CE). Four erratic boulders, 10–250 m distal from late Holocene moraines, yielded 10Be exposure ages of 10.9–11.6 ka, demonstrating that by ca. 11 ka, alpine glaciers were no more extensive than during the last several hundred years. Estimated temperature change obtained through reconstruction of equilibrium line altitudes shows that since ca. 1850 CE, mean annual temperatures have risen 0.2–2.3 ∘C. We use our glacier chronology and the Open Global Glacier Model (OGGM) to estimate that from 1000 CE, glaciers in this region reached a maximum total volume of 34–38 km3 between 1765 and 1855 CE and had lost nearly half their ice volume by 2019 CE. OGGM was unable to produce modeled glacier lengths that match the timing or magnitude of the maximum glacier extent indicated by the 10Be chronology. However, when applied to the entire Mackenzie and Selwyn mountain region, past millennium OGGM simulations using the Max Planck Institute Earth System Model (MPI-ESM) and the Community Climate System Model 4 (CCSM4) yield late Holocene glacier volume change temporally consistent with our moraine and remote sensing record, while the Meteorological Research Institute Earth System Model 2 (MRI-ESM2) and the Model for Interdisciplinary Research on Climate (MIROC) fail to produce modeled glacier change consistent with our glacier chronology. Finally, OGGM forced by future climate projections under varying greenhouse gas emission scenarios predicts 85 % to over 97 % glacier volume loss by the end of the 21st century. The loss of glaciers from this region will have profound impacts on local ecosystems and communities that rely on meltwater from glacierized catchments.

     
    more » « less
  3. Abstract

    The expansion of shrubs across the Arctic tundra may fundamentally modify land–atmosphere interactions. However, it remains unclear how shrub expansion pattern is linked with key environmental drivers, such as climate change and fire disturbance. Here we used 40+ years of high‐resolution (~1.0 m) aerial and satellite imagery to estimate shrub‐cover change in 114 study sites across four burned and unburned upland (ice‐poor) and lowland (ice‐rich) tundra ecosystems in northern Alaska. Validated with data from four additional upland and lowland tundra fires, our results reveal that summer precipitation was the most important climatic driver (r = 0.67,p < 0.001), responsible for 30.8% of shrub expansion in the upland tundra between 1971 and 2016. Shrub expansion in the uplands was largely enhanced by wildfire (p < 0.001) and it exhibited positive correlation with fire severity (r = 0.83,p < 0.001). Three decades after fire disturbance, the upland shrub cover increased by 1077.2 ± 83.6 m2 ha−1, ~7 times the amount identified in adjacent unburned upland tundra (155.1 ± 55.4 m2 ha−1). In contrast, shrub cover markedly decreased in lowland tundra after fire disturbance, which triggered thermokarst‐associated water impounding and resulted in 52.4% loss of shrub cover over three decades. No correlation was found between lowland shrub cover with fire severity (r = 0.01). Mean summer air temperature (MSAT) was the principal factor driving lowland shrub‐cover dynamics between 1951 and 2007. Warmer MSAT facilitated shrub expansion in unburned lowlands (r = 0.78,p < 0.001), but accelerated shrub‐cover losses in burned lowlands (r = −0.82,p < 0.001). These results highlight divergent pathways of shrub‐cover responses to fire disturbance and climate change, depending on near‐surface permafrost and drainage conditions. Our study offers new insights into the land–atmosphere interactions as climate warming and burning intensify in high latitudes.

     
    more » « less
  4. Abstract A mass-balance model using upper-air meteorological data for input was calibrated with surface mass balance measured mainly during 1977–78 at 67 sites on Columbia Glacier, Alaska, between 135 and 2645 m a.s.l. Root-mean-square error, model vs measured, is 1.0 m w.e. a −1 , with r 2 = 0.88. A remarkable result of the analysis was that both precipitation and the factor in the positive degree-day model used to estimate surface ablation were constant with altitude. The model was applied to reconstruct glacier-wide components of surface mass balance over 1948–2007. Surface ablation, 4 km 3 ice eq. a −1 (ice equivalent), has changed little throughout the period. From 1948 until about 1981, when drastic retreat began, the surface mass balance was positive but changes in glacier geometry were small, so the positive balance was offset by calving, ∼0.9 km 3 ice eq. a −1 . During retreat, volume loss of the glacier accounted for 92% of the iceberg production. Calving increased to ∼4.3 km 3 ice eq. a −1 from 1982 to 1995, and after that until 2007 to ∼8.0 km 3 ice eq. a −1 , which was about twice the loss by surface ablation, whereas prior to retreat it was only about a quarter as much. Calving is calculated as the difference between glacier-wide surface mass balance and geodetically determined volume change. 
    more » « less
  5. Abstract

    Globally, glaciers are shrinking in response to climate change, with implications for global sea level rise as well as downstream ecosystems and water resources. Sliding at the ice‐bed interface (basal motion) provides a mechanism for glaciers to respond rapidly to climate change. While the short‐term dynamics of glacier basal motion (<10 years) have received substantial attention, little is known about how basal motion and its sensitivity to subglacial hydrology changes over long (>50 year) timescales—this knowledge is required for accurate prediction of future glacier change. We compare historical data with modern estimates from field and satellite data at Athabasca Glacier and show that the glacier thinned by 60 m (−21%) over 1961–2020. However, a concurrent increase in surface slope results in minimal change in the average driving stress (−6 kPa and −4%). These geometric changes coincide with relatively uniform slowing (−15 m a−1and −45%). Simplified ice modeling suggests that declining basal motion accounts for most of this slow down (91% on average and 46% at minimum). A decline in basal motion can be explained by increasing basal friction resulting from geometric change in addition to increasing meltwater flux through a more efficient subglacial hydrologic system. These results highlight the need to include time‐varying dynamics of basal motion in glacier models and analyses. If these findings are generalizable, they suggest that declining basal motion reduces the flux of ice to lower elevations, helping to mitigate glacier mass loss in a warming climate.

     
    more » « less