Widespread seawater intrusions beneath the grounded ice of Thwaites Glacier, West Antarctica.

Warm water from the Southern Ocean has a dominant impact on the evolution of Antarctic glaciers and in turn on their contribution to sea level rise. Using a continuous time series of daily-repeat satellite synthetic-aperture radar interferometry data from the ICEYE constellation collected in March-June 2023, we document an ice grounding zone, or region of tidally controlled migration of the transition boundary between grounded ice and ice afloat in the ocean, at the main trunk of Thwaites Glacier, West Antarctica, a strong contributor to sea level rise with an ice volume equivalent to a 0.6-m global sea level rise. The ice grounding zone is 6 km wide in the central part of Thwaites with shallow bed slopes, and 2 km wide along its flanks with steep basal slopes. We additionally detect irregular seawater intrusions, 5 to 10 cm in thickness, extending another 6 km upstream, at high tide, in a bed depression located beyond a bedrock ridge that impedes the glacier retreat. Seawater intrusions align well with regions predicted by the GlaDS subglacial water model to host a high-pressure distributed subglacial hydrology system in between lower-pressure subglacial channels. Pressurized seawater intrusions will induce vigorous melt of grounded ice over kilometers, making the glacier more vulnerable to ocean warming, and increasing the projections of ice mass loss. Kilometer-wide, widespread seawater intrusion beneath grounded ice may be the missing link between the rapid, past, and present changes in ice sheet mass and the slower changes replicated by ice sheet models.

Melt rates in the kilometer-size grounding zone of Petermann Glacier, Greenland, before and during a retreat.

Warming of the ocean waters surrounding Greenland plays a major role in driving glacier retreat and the contribution of glaciers to sea level rise. The melt rate at the junction of the ocean with grounded ice-or grounding line-is, however, not well known. Here, we employ a time series of satellite radar interferometry data from the German TanDEM-X mission, the Italian COSMO-SkyMed constellation, and the Finnish ICEYE constellation to document the grounding line migration and basal melt rates of Petermann Glacier, a major marine-based glacier of Northwest Greenland. We find that the grounding line migrates at tidal frequencies over a kilometer-wide (2 to 6 km) grounding zone, which is one order of magnitude larger than expected for grounding lines on a rigid bed. The highest ice shelf melt rates are recorded within the grounding zone with values from 60 ± 13 to 80 ± 15 m/y along laterally confined channels. As the grounding line retreated by 3.8 km in 2016 to 2022, it carved a cavity about 204 m in height where melt rates increased from 40 ± 11 m/y in 2016 to 2019 to 60 ± 15 m/y in 2020 to 2021. In 2022, the cavity remained open during the entire tidal cycle. Such high melt rates concentrated in kilometer-wide grounding zones contrast with the traditional plume model of grounding line melt which predicts zero melt. High rates of simulated basal melting in grounded glacier ice in numerical models will increase the glacier sensitivity to ocean warming and potentially double projections of sea level rise.

Automatic delineation of glacier grounding lines in differential interferometric synthetic-aperture radar data using deep learning.

Delineating the grounding line of marine-terminating glaciers-where ice starts to become afloat in ocean waters-is crucial for measuring and understanding ice sheet mass balance, glacier dynamics, and their contributions to sea level rise. This task has been previously done using time-consuming, mostly-manual digitizations of differential interferometric synthetic-aperture radar interferograms by human experts. This approach is no longer viable with a fast-growing set of satellite observations and the need to establish time series over entire continents with quantified uncertainties. We present a fully-convolutional neural network with parallel atrous convolutional layers and asymmetric encoder/decoder components that automatically delineates grounding lines at a large scale, efficiently, and accompanied by uncertainty estimates. Our procedure detects grounding lines within 232 m in 100-m posting interferograms, which is comparable to the performance achieved by human experts. We also find value in the machine learning approach in situations that even challenge human experts. We use this approach to map the tidal-induced variability in grounding line position around Antarctica in 22,935 interferograms from year 2018. Along the Getz Ice Shelf, in West Antarctica, we demonstrate that grounding zones are one order magnitude (13.3 ± 3.9) wider than expected from hydrostatic equilibrium, which justifies the need to map grounding lines repeatedly and comprehensively to inform numerical 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.

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.

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.

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.