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.

A Late Paleocene age for Greenland's Hiawatha impact structure.

The ~31-km-wide Hiawatha structure, located beneath Hiawatha Glacier in northwestern Greenland, has been proposed as an impact structure that may have formed after the Pleistocene inception of the Greenland Ice Sheet. To date the structure, we conducted 40Ar/39Ar analyses on glaciofluvial sand and U-Pb analyses on zircon separated from glaciofluvial pebbles of impact melt rock, all sampled immediately downstream of Hiawatha Glacier. Unshocked zircon in the impact melt rocks dates to ~1915 million years (Ma), consistent with felsic intrusions found in local bedrock. The 40Ar/39Ar data indicate Late Paleocene resetting and shocked zircon dates to 57.99 ± 0.54 Ma, which we interpret as the impact age. Consequently, the Hiawatha impact structure far predates Pleistocene glaciation and is unrelated to either the Paleocene-Eocene Thermal Maximum or flood basalt volcanism in east Greenland. However, it was contemporaneous with the Paleocene Carbon Isotope Maximum, although the impact's exact paleoenvironmental and climatic significance awaits further investigation.

Retreat of Humboldt Gletscher, North Greenland, Driven by Undercutting From a Warmer Ocean.

Humboldt Gletscher is a 100-km wide, slow-moving glacier in north Greenland which holds a 19-cm global sea level equivalent. Humboldt has been the fourth largest contributor to sea level rise since 1972 but the cause of its mass loss has not been elucidated. Multi-beam echo sounding data collected in 2019 indicate a seabed 200 m deeper than previously known. Conductivity temperature depth data reveal the presence of warm water of Atlantic origin at 0°C at the glacier front and a warming of the ocean waters by 0.9 ± 0.1°C since 1962. Using an ocean model, we reconstruct grounded ice undercutting by the ocean, combine it with calculated retreat caused by ice thinning to floatation, and are able to fully explain the observed retreat. Two thirds of the retreat are caused by undercutting of grounded ice, which is a physical process not included in most ice sheet models.

Ocean forcing drives glacier retreat in Greenland.

The retreat and acceleration of Greenland glaciers since the mid-1990s have been attributed to the enhanced intrusion of warm Atlantic Waters (AW) into fjords, but this assertion has not been quantitatively tested on a Greenland-wide basis or included in models. Here, we investigate how AW influenced retreat at 226 marine-terminating glaciers using ocean modeling, remote sensing, and in situ observations. We identify 74 glaciers in deep fjords with AW controlling 49% of the mass loss that retreated when warming increased undercutting by 48%. Conversely, 27 glaciers calving on shallow ridges and 24 in cold, shallow waters retreated little, contributing 15% of the loss, while 10 glaciers retreated substantially following the collapse of several ice shelves. The retreat mechanisms remain undiagnosed at 87 glaciers without ocean and bathymetry data, which controlled 19% of the loss. Ice sheet projections that exclude ocean-induced undercutting may underestimate mass loss by at least a factor of 2.

Ocean melting of the Zachariae Isstrøm and Nioghalvfjerdsfjorden glaciers, northeast Greenland.

Zachariae Isstrøm (ZI) and Nioghalvfjerdsfjorden (79N) are marine-terminating glaciers in northeast Greenland that hold an ice volume equivalent to a 1.1-m global sea level rise. ZI lost its floating ice shelf, sped up, retreated at 650 m/y, and experienced a 5-gigaton/y mass loss. Glacier 79N has been more stable despite its exposure to the same climate forcing. We analyze the impact of ocean thermal forcing on the glaciers. A three-dimensional inversion of airborne gravity data reveals an 800-m-deep, broad channel that allows subsurface, warm, Atlantic Intermediate Water (AIW) (+1.[Formula: see text]C) to reach the front of ZI via two sills at 350-m depth. Subsurface ocean temperature in that channel has warmed by 1.3[Formula: see text]C since 1979. Using an ocean model, we calculate a rate of ice removal at the grounding line by the ocean that increased from 108 m/y to 185 m/y in 1979-2019. Observed ice thinning caused a retreat of its flotation line to increase from 105 m/y to 217 m/y, for a combined grounding line retreat of 13 km in 41 y that matches independent observations within 14%. In contrast, the limited access of AIW to 79N via a narrower passage yields lower grounded ice removal (53 m/y to 99 m/y) and thinning-induced retreat (27 m/y to 50 m/y) for a combined retreat of 4.4 km, also within 12% of observations. Ocean-induced removal of ice at the grounding line, modulated by bathymetric barriers, is therefore a main driver of ice sheet retreat, but it is not incorporated in most ice sheet models.

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.

Validation of Glacier Topographic Acquisitions from an Airborne Single-Pass Interferometer.

The airborne glacier and ice surface topography interferometer (GLISTIN-A) is a single-pass radar interferometer developed for accurate high-resolution swath mapping of dynamic ice surfaces. We present the first validation results of the operational sensor, collected in 2013 over glaciers in Alaska and followed by more exhaustive collections from Greenland in 2016 and 2017. In Alaska, overlapping flight-tracks were mosaicked to mitigate potential residual trends across-track and the resultant maps are validated with lidar. Furthermore, repeat acquisitions of Columbia Glacier collected with a three day separation indicate excellent stability and repeatability. Commencing 2016, GLISTIN-A has circumnavigated Greenland for 4 consecutive years. Due to flight hour limitations, overlapping swaths were not flown. In 2016, comparison with airborne lidar data finds that residual systematic errors exhibit evenly distributed small slopes (all less than 10 millidegrees) and nadir biases were typically less than 1 m. Similarly 2017 data exhibited up to meter-scale nadir biases and evenly distributed residual slopes with a standard deviation of ~10 millidegrees). All satisfied the science accuracy requirements of the Greenland campaigns (3 m accuracy across an 8 km swath).

Forty-six years of Greenland Ice Sheet mass balance from 1972 to 2018.

We reconstruct the mass balance of the Greenland Ice Sheet using a comprehensive survey of thickness, surface elevation, velocity, and surface mass balance (SMB) of 260 glaciers from 1972 to 2018. We calculate mass discharge, D, into the ocean directly for 107 glaciers (85% of D) and indirectly for 110 glaciers (15%) using velocity-scaled reference fluxes. The decadal mass balance switched from a mass gain of +47 ± 21 Gt/y in 1972-1980 to a loss of 51 ± 17 Gt/y in 1980-1990. The mass loss increased from 41 ± 17 Gt/y in 1990-2000, to 187 ± 17 Gt/y in 2000-2010, to 286 ± 20 Gt/y in 2010-2018, or sixfold since the 1980s, or 80 ± 6 Gt/y per decade, on average. The acceleration in mass loss switched from positive in 2000-2010 to negative in 2010-2018 due to a series of cold summers, which illustrates the difficulty of extrapolating short records into longer-term trends. Cumulated since 1972, the largest contributions to global sea level rise are from northwest (4.4 ± 0.2 mm), southeast (3.0 ± 0.3 mm), and central west (2.0 ± 0.2 mm) Greenland, with a total 13.7 ± 1.1 mm for the ice sheet. The mass loss is controlled at 66 ± 8% by glacier dynamics (9.1 mm) and 34 ± 8% by SMB (4.6 mm). Even in years of high SMB, enhanced glacier discharge has remained sufficiently high above equilibrium to maintain an annual mass loss every year since 1998.

Four decades of Antarctic Ice Sheet mass balance from 1979-2017.

We use updated drainage inventory, ice thickness, and ice velocity data to calculate the grounding line ice discharge of 176 basins draining the Antarctic Ice Sheet from 1979 to 2017. We compare the results with a surface mass balance model to deduce the ice sheet mass balance. The total mass loss increased from 40 ± 9 Gt/y in 1979-1990 to 50 ± 14 Gt/y in 1989-2000, 166 ± 18 Gt/y in 1999-2009, and 252 ± 26 Gt/y in 2009-2017. In 2009-2017, the mass loss was dominated by the Amundsen/Bellingshausen Sea sectors, in West Antarctica (159 ± 8 Gt/y), Wilkes Land, in East Antarctica (51 ± 13 Gt/y), and West and Northeast Peninsula (42 ± 5 Gt/y). The contribution to sea-level rise from Antarctica averaged 3.6 ± 0.5 mm per decade with a cumulative 14.0 ± 2.0 mm since 1979, including 6.9 ± 0.6 mm from West Antarctica, 4.4 ± 0.9 mm from East Antarctica, and 2.5 ± 0.4 mm from the Peninsula (i.e., East Antarctica is a major participant in the mass loss). During the entire period, the mass loss concentrated in areas closest to warm, salty, subsurface, circumpolar deep water (CDW), that is, consistent with enhanced polar westerlies pushing CDW toward Antarctica to melt its floating ice shelves, destabilize the glaciers, and raise sea level.

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.

Rapid submarine ice melting in the grounding zones of ice shelves in West Antarctica.

Enhanced submarine ice-shelf melting strongly controls ice loss in the Amundsen Sea embayment (ASE) of West Antarctica, but its magnitude is not well known in the critical grounding zones of the ASE's major glaciers. Here we directly quantify bottom ice losses along tens of kilometres with airborne radar sounding of the Dotson and Crosson ice shelves, which buttress the rapidly changing Smith, Pope and Kohler glaciers. Melting in the grounding zones is found to be much higher than steady-state levels, removing 300-490 m of solid ice between 2002 and 2009 beneath the retreating Smith Glacier. The vigorous, unbalanced melting supports the hypothesis that a significant increase in ocean heat influx into ASE sub-ice-shelf cavities took place in the mid-2000s. The synchronous but diverse evolutions of these glaciers illustrate how combinations of oceanography and topography modulate rapid submarine melting to hasten mass loss and glacier retreat from West Antarctica.

Observed latitudinal variations in erosion as a function of glacier dynamics.

Glacial erosion is fundamental to our understanding of the role of Cenozoic-era climate change in the development of topography worldwide, yet the factors that control the rate of erosion by ice remain poorly understood. In many tectonically active mountain ranges, glaciers have been inferred to be highly erosive, and conditions of glaciation are used to explain both the marked relief typical of alpine settings and the limit on mountain heights above the snowline, that is, the glacial buzzsaw. In other high-latitude regions, glacial erosion is presumed to be minimal, where a mantle of cold ice effectively protects landscapes from erosion. Glacial erosion rates are expected to increase with decreasing latitude, owing to the climatic control on basal temperature and the production of meltwater, which promotes glacial sliding, erosion and sediment transfer. This relationship between climate, glacier dynamics and erosion rate is the focus of recent numerical modelling, yet it is qualitative and lacks an empirical database. Here we present a comprehensive data set that permits explicit examination of the factors controlling glacier erosion across climatic regimes. We report contemporary ice fluxes, sliding speeds and erosion rates inferred from sediment yields from 15 outlet glaciers spanning 19 degrees of latitude from Patagonia to the Antarctic Peninsula. Although this broad region has a relatively uniform tectonic and geologic history, the thermal regimes of its glaciers range from temperate to polar. We find that basin-averaged erosion rates vary by three orders of magnitude over this latitudinal transect. Our findings imply that climate and the glacier thermal regime control erosion rates more than do extent of ice cover, ice flux or sliding speeds.

A reconciled estimate of ice-sheet mass balance.

We combined an ensemble of satellite altimetry, interferometry, and gravimetry data sets using common geographical regions, time intervals, and models of surface mass balance and glacial isostatic adjustment to estimate the mass balance of Earth's polar ice sheets. We find that there is good agreement between different satellite methods--especially in Greenland and West Antarctica--and that combining satellite data sets leads to greater certainty. Between 1992 and 2011, the ice sheets of Greenland, East Antarctica, West Antarctica, and the Antarctic Peninsula changed in mass by -142 ± 49, +14 ± 43, -65 ± 26, and -20 ± 14 gigatonnes year(-1), respectively. Since 1992, the polar ice sheets have contributed, on average, 0.59 ± 0.20 millimeter year(-1) to the rate of global sea-level rise.

Accumulation and erosion of Mars' south polar layered deposits.

Mars' polar regions are covered with ice-rich layered deposits that potentially contain a record of climate variations. The sounding radar SHARAD on the Mars Reconnaissance Orbiter mapped detailed subsurface stratigraphy in the Promethei Lingula region of the south polar plateau, Planum Australe. Radar reflections interpreted as layers are correlated across adjacent orbits and are continuous for up to 150 kilometers along spacecraft orbital tracks. The reflectors are often separated into discrete reflector sequences, and strong echoes are seen as deep as 1 kilometer. In some cases, the sequences are dipping with respect to each other, suggesting an interdepositional period of erosion. In Australe Sulci, layers are exhumed, indicating recent erosion.

Subsurface radar sounding of the south polar layered deposits of Mars.

The ice-rich south polar layered deposits of Mars were probed with the Mars Advanced Radar for Subsurface and Ionospheric Sounding on the Mars Express orbiter. The radar signals penetrate deep into the deposits (more than 3.7 kilometers). For most of the area, a reflection is detected at a time delay that is consistent with an interface between the deposits and the substrate. The reflected power from this interface indicates minimal attenuation of the signal, suggesting a composition of nearly pure water ice. Maps were generated of the topography of the basal interface and the thickness of the layered deposits. A set of buried depressions is seen within 300 kilometers of the pole. The thickness map shows an asymmetric distribution of the deposits and regions of anomalous thickness. The total volume is estimated to be 1.6 x 10(6) cubic kilometers, which is equivalent to a global water layer approximately 11 meters thick.