Extensive inland thinning and speed-up of Northeast Greenland Ice Stream.

Over the past two decades, ice loss from the Greenland ice sheet (GrIS) has increased owing to enhanced surface melting and ice discharge to the ocean1-5. Whether continuing increased ice loss will accelerate further, and by how much, remains contentious6-9. A main contributor to future ice loss is the Northeast Greenland Ice Stream (NEGIS), Greenland's largest basin and a prominent feature of fast-flowing ice that reaches the interior of the GrIS10-12. Owing to its topographic setting, this sector is vulnerable to rapid retreat, leading to unstable conditions similar to those in the marine-based setting of ice streams in Antarctica13-20. Here we show that extensive speed-up and thinning triggered by frontal changes in 2012 have already propagated more than 200 km inland. We use unique global navigation satellite system (GNSS) observations, combined with surface elevation changes and surface speeds obtained from satellite data, to select the correct basal conditions to be used in ice flow numerical models, which we then use for future simulations. Our model results indicate that this marine-based sector alone will contribute 13.5-15.5 mm sea-level rise by 2100 (equivalent to the contribution of the entire ice sheet over the past 50 years) and will cause precipitous changes in the coming century. This study shows that measurements of subtle changes in the ice speed and elevation inland help to constrain numerical models of the future mass balance and higher-end projections show better agreement with observations.

Three-dimensional topology dataset of folded radar stratigraphy in northern Greenland.

We present a dataset of reconstructed three-dimensional (3D) englacial stratigraphic horizons in northern Greenland. The data cover four different regions representing key ice-dynamic settings in Greenland: (i) the onset of Petermann Glacier, (ii) a region upstream of the 79° North Glacier (Nioghalvfjerdsbræ), near the northern Greenland ice divide, (iii) the onset of the Northeast Greenland Ice Stream (NEGIS) and (iv) a 700 km wide region extending across the central ice divide over the entire northern part of central Greenland. In this paper, we promote the advantages of a 3D perspective of deformed englacial stratigraphy and explain how 3D horizons provide an improved basis for interpreting and reconstructing the ice-dynamic history. The 3D horizons are provided in various formats to allow a wide range of applications and reproducibility of results.

Crystal orientation fabric anisotropy causes directional hardening of the Northeast Greenland Ice Stream.

The dynamic mass loss of ice sheets constitutes one of the biggest uncertainties in projections of ice-sheet evolution. One central, understudied aspect of ice flow is how the bulk orientation of the crystal orientation fabric translates to the mechanical anisotropy of ice. Here we show the spatial distribution of the depth-averaged horizontal anisotropy and corresponding directional flow-enhancement factors covering a large area of the Northeast Greenland Ice Stream onset. Our results are based on airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal anisotropy and a rapid crystal reorganisation on the order of hundreds of years coinciding with the ice-stream geometry. Compared to isotropic ice, parts of the ice stream are found to be more than one order of magnitude harder for along-flow extension/compression while the shear margins are potentially softened by a factor of two for horizontal-shear deformation.

Greenland Mass Trends From Airborne and Satellite Altimetry During 2011-2020.

We use satellite and airborne altimetry to estimate annual mass changes of the Greenland Ice Sheet. We estimate ice loss corresponding to a sea-level rise of 6.9 ± 0.4 mm from April 2011 to April 2020, with a highest annual ice loss rate of 1.4 mm/yr sea-level equivalent from April 2019 to April 2020. On a regional scale, our annual mass loss timeseries reveals 10-15 m/yr dynamic thickening at the terminus of Jakobshavn Isbræ from April 2016 to April 2018, followed by a return to dynamic thinning. We observe contrasting patterns of mass loss acceleration in different basins across the ice sheet and suggest that these spatiotemporal trends could be useful for calibrating and validating prognostic ice sheet models. In addition to resolving the spatial and temporal fingerprint of Greenland's recent ice loss, these mass loss grids are key for partitioning contemporary elastic vertical land motion from longer-term glacial isostatic adjustment (GIA) trends at GPS stations around the ice sheet. Our ice-loss product results in a significantly different GIA interpretation from a previous ice-loss product.

Centennial response of Greenland's three largest outlet glaciers.

The Greenland Ice Sheet is the largest land ice contributor to sea level rise. This will continue in the future but at an uncertain rate and observational estimates are limited to the last few decades. Understanding the long-term glacier response to external forcing is key to improving projections. Here we use historical photographs to calculate ice loss from 1880-2012 for Jakobshavn, Helheim, and Kangerlussuaq glacier. We estimate ice loss corresponding to a sea level rise of 8.1 ± 1.1 millimetres from these three glaciers. Projections of mass loss for these glaciers, using the worst-case scenario, Representative Concentration Pathways 8.5, suggest a sea level contribution of 9.1-14.9 mm by 2100. RCP8.5 implies an additional global temperature increase of 3.7 °C by 2100, approximately four times larger than that which has taken place since 1880. We infer that projections forced by RCP8.5 underestimate glacier mass loss which could exceed this worst-case scenario.

A large impact crater beneath Hiawatha Glacier in northwest Greenland.

We report the discovery of a large impact crater beneath Hiawatha Glacier in northwest Greenland. From airborne radar surveys, we identify a 31-kilometer-wide, circular bedrock depression beneath up to a kilometer of ice. This depression has an elevated rim that cross-cuts tributary subglacial channels and a subdued central uplift that appears to be actively eroding. From ground investigations of the deglaciated foreland, we identify overprinted structures within Precambrian bedrock along the ice margin that strike tangent to the subglacial rim. Glaciofluvial sediment from the largest river draining the crater contains shocked quartz and other impact-related grains. Geochemical analysis of this sediment indicates that the impactor was a fractionated iron asteroid, which must have been more than a kilometer wide to produce the identified crater. Radiostratigraphy of the ice in the crater shows that the Holocene ice is continuous and conformable, but all deeper and older ice appears to be debris rich or heavily disturbed. The age of this impact crater is presently unknown, but from our geological and geophysical evidence, we conclude that it is unlikely to predate the Pleistocene inception of the Greenland Ice Sheet.

Spatial and temporal Antarctic Ice Sheet mass trends, glacio-isostatic adjustment, and surface processes from a joint inversion of satellite altimeter, gravity, and GPS data.

We present spatiotemporal mass balance trends for the Antarctic Ice Sheet from a statistical inversion of satellite altimetry, gravimetry, and elastic-corrected GPS data for the period 2003-2013. Our method simultaneously determines annual trends in ice dynamics, surface mass balance anomalies, and a time-invariant solution for glacio-isostatic adjustment while remaining largely independent of forward models. We establish that over the period 2003-2013, Antarctica has been losing mass at a rate of -84 ± 22 Gt yr-1, with a sustained negative mean trend of dynamic imbalance of -111 ± 13 Gt yr-1. West Antarctica is the largest contributor with -112 ± 10 Gt yr-1, mainly triggered by high thinning rates of glaciers draining into the Amundsen Sea Embayment. The Antarctic Peninsula has experienced a dramatic increase in mass loss in the last decade, with a mean rate of -28 ± 7 Gt yr-1 and significantly higher values for the most recent years following the destabilization of the Southern Antarctic Peninsula around 2010. The total mass loss is partly compensated by a significant mass gain of 56 ± 18 Gt yr-1 in East Antarctica due to a positive trend of surface mass balance anomalies.

Geodetic measurements reveal similarities between post-Last Glacial Maximum and present-day mass loss from the Greenland ice sheet.

Accurate quantification of the millennial-scale mass balance of the Greenland ice sheet (GrIS) and its contribution to global sea-level rise remain challenging because of sparse in situ observations in key regions. Glacial isostatic adjustment (GIA) is the ongoing response of the solid Earth to ice and ocean load changes occurring since the Last Glacial Maximum (LGM; ~21 thousand years ago) and may be used to constrain the GrIS deglaciation history. We use data from the Greenland Global Positioning System network to directly measure GIA and estimate basin-wide mass changes since the LGM. Unpredicted, large GIA uplift rates of +12 mm/year are found in southeast Greenland. These rates are due to low upper mantle viscosity in the region, from when Greenland passed over the Iceland hot spot about 40 million years ago. This region of concentrated soft rheology has a profound influence on reconstructing the deglaciation history of Greenland. We reevaluate the evolution of the GrIS since LGM and obtain a loss of 1.5-m sea-level equivalent from the northwest and southeast. These same sectors are dominating modern mass loss. We suggest that the present destabilization of these marine-based sectors may increase sea level for centuries to come. Our new deglaciation history and GIA uplift estimates suggest that studies that use the Gravity Recovery and Climate Experiment satellite mission to infer present-day changes in the GrIS may have erroneously corrected for GIA and underestimated the mass loss by about 20 gigatons/year.