Accelerated global glacier mass loss in the early twenty-first century.

Glaciers distinct from the Greenland and Antarctic ice sheets are shrinking rapidly, altering regional hydrology1, raising global sea level2 and elevating natural hazards3. Yet, owing to the scarcity of constrained mass loss observations, glacier evolution during the satellite era is known only partially, as a geographic and temporal patchwork4,5. Here we reveal the accelerated, albeit contrasting, patterns of glacier mass loss during the early twenty-first century. Using largely untapped satellite archives, we chart surface elevation changes at a high spatiotemporal resolution over all of Earth's glaciers. We extensively validate our estimates against independent, high-precision measurements and present a globally complete and consistent estimate of glacier mass change. We show that during 2000-2019, glaciers lost a mass of 267 ± 16 gigatonnes per year, equivalent to 21 ± 3 per cent of the observed sea-level rise6. We identify a mass loss acceleration of 48 ± 16 gigatonnes per year per decade, explaining 6 to 19 per cent of the observed acceleration of sea-level rise. Particularly, thinning rates of glaciers outside ice sheet peripheries doubled over the past two decades. Glaciers currently lose more mass, and at similar or larger acceleration rates, than the Greenland or Antarctic ice sheets taken separately7-9. By uncovering the patterns of mass change in many regions, we find contrasting glacier fluctuations that agree with the decadal variability in precipitation and temperature. These include a North Atlantic anomaly of decelerated mass loss, a strongly accelerated loss from northwestern American glaciers, and the apparent end of the Karakoram anomaly of mass gain10. We anticipate our highly resolved estimates to advance the understanding of drivers that govern the distribution of glacier change, and to extend our capabilities of predicting these changes at all scales. Predictions robustly benchmarked against observations are critically needed to design adaptive policies for the local- and regional-scale management of water resources and cryospheric risks, as well as for the global-scale mitigation of sea-level rise.

Dynamic vulnerability revealed in the collapse of an Arctic tidewater glacier.

Glacier flow instabilities can rapidly increase sea level through enhanced ice discharge. Surge-type glacier accelerations often occur with a decadal to centennial cyclicity suggesting internal mechanisms responsible. Recently, many surging tidewater glaciers around the Arctic Barents Sea region question whether external forces such as climate can trigger dynamic instabilities. Here, we identify a mechanism in which climate change can instigate surges of Arctic tidewater glaciers. Using satellite and seismic remote sensing observations combined with three-dimensional thermo-mechanical modeling of the January 2009 collapse of the Nathorst Glacier System (NGS) in Svalbard, we show that an underlying condition for instability was basal freezing and associated friction increase under the glacier tongue. In contrast, continued basal sliding further upstream increased driving stresses until eventual and sudden till failure under the tongue. The instability propagated rapidly up-glacier, mobilizing the entire 450 km2 glacier basin over a few days as the till entered an unstable friction regime. Enhanced mass loss during and after the collapse (5-7 fold compared to pre-collapse mass losses) combined with regionally rising equilibrium line altitudes strongly limit mass replenishment of the glacier, suggesting irreversible consequences. Climate plays a paradoxical role as cold glacier thinning and retreat promote basal freezing which increases friction at the tongue by stabilizing an efficient basal drainage system. However, with some of the most intense atmospheric warming on Earth occurring in the Arctic, increased melt water can reduce till strength under tidewater glacier tongues to orchestrate a temporal clustering of surges at decadal timescales, such as those observed in Svalbard at the end of the Little Ice Age. Consequently, basal terminus freezing promotes a dynamic vulnerability to climate change that may be present in many Arctic tidewater glaciers.

Digital elevation model and orthophotographs of Greenland based on aerial photographs from 1978-1987.

Digital Elevation Models (DEMs) play a prominent role in glaciological studies for the mass balance of glaciers and ice sheets. By providing a time snapshot of glacier geometry, DEMs are crucial for most glacier evolution modelling studies, but are also important for cryospheric modelling in general. We present a historical medium-resolution DEM and orthophotographs that consistently cover the entire surroundings and margins of the Greenland Ice Sheet 1978-1987. About 3,500 aerial photographs of Greenland are combined with field surveyed geodetic ground control to produce a 25 m gridded DEM and a 2 m black-and-white digital orthophotograph. Supporting data consist of a reliability mask and a photo footprint coverage with recording dates. Through one internal and two external validation tests, this DEM shows an accuracy better than 10 m horizontally and 6 m vertically while the precision is better than 4 m. This dataset proved successful for topographical mapping and geodetic mass balance. Other uses include control and calibration of remotely sensed data such as imagery or InSAR velocity maps.

Spatial and temporal distribution of mass loss from the Greenland Ice Sheet since AD 1900.

The response of the Greenland Ice Sheet (GIS) to changes in temperature during the twentieth century remains contentious, largely owing to difficulties in estimating the spatial and temporal distribution of ice mass changes before 1992, when Greenland-wide observations first became available. The only previous estimates of change during the twentieth century are based on empirical modelling and energy balance modelling. Consequently, no observation-based estimates of the contribution from the GIS to the global-mean sea level budget before 1990 are included in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Here we calculate spatial ice mass loss around the entire GIS from 1900 to the present using aerial imagery from the 1980s. This allows accurate high-resolution mapping of geomorphic features related to the maximum extent of the GIS during the Little Ice Age at the end of the nineteenth century. We estimate the total ice mass loss and its spatial distribution for three periods: 1900-1983 (75.1 ± 29.4 gigatonnes per year), 1983-2003 (73.8 ± 40.5 gigatonnes per year), and 2003-2010 (186.4 ± 18.9 gigatonnes per year). Furthermore, using two surface mass balance models we partition the mass balance into a term for surface mass balance (that is, total precipitation minus total sublimation minus runoff) and a dynamic term. We find that many areas currently undergoing change are identical to those that experienced considerable thinning throughout the twentieth century. We also reveal that the surface mass balance term shows a considerable decrease since 2003, whereas the dynamic term is constant over the past 110 years. Overall, our observation-based findings show that during the twentieth century the GIS contributed at least 25.0 ± 9.4 millimetres of global-mean sea level rise. Our result will help to close the twentieth-century sea level budget, which remains crucial for evaluating the reliability of models used to predict global sea level rise.

Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas.

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.