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.

Suppressed basal melting in the eastern Thwaites Glacier grounding zone.

Thwaites Glacier is one of the fastest-changing ice-ocean systems in Antarctica1-3. Much of the ice sheet within the catchment of Thwaites Glacier is grounded below sea level on bedrock that deepens inland4, making it susceptible to rapid and irreversible ice loss that could raise the global sea level by more than half a metre2,3,5. The rate and extent of ice loss, and whether it proceeds irreversibly, are set by the ocean conditions and basal melting within the grounding-zone region where Thwaites Glacier first goes afloat3,6, both of which are largely unknown. Here we show-using observations from a hot-water-drilled access hole-that the grounding zone of Thwaites Eastern Ice Shelf (TEIS) is characterized by a warm and highly stable water column with temperatures substantially higher than the in situ freezing point. Despite these warm conditions, low current speeds and strong density stratification in the ice-ocean boundary layer actively restrict the vertical mixing of heat towards the ice base7,8, resulting in strongly suppressed basal melting. Our results demonstrate that the canonical model of ice-shelf basal melting used to generate sea-level projections cannot reproduce observed melt rates beneath this critically important glacier, and that rapid and possibly unstable grounding-line retreat may be associated with relatively modest basal melt rates.

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.

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.

Storstrømmen and L. Bistrup Bræ, North Greenland, Protected From Warm Atlantic Ocean Waters.

Storstrømmen and L. Bistrup Bræ are 20- and 10-km wide, surge type glaciers in North Greenland in quiescent phase that terminate in the southernmost floating ice tongue in East Greenland. Novel multi-beam echo sounding data collected in August 2020 indicate a seabed at 350-400 m depth along a relatively uniform ice shelf front, 100 m deeper than expected, but surrounded by shallower terrain (<100 m) over a 30-km wide region that blocks the access of warm, salty, subsurface Atlantic Intermediate Water (AIW) at +1.6°C. Conductivity temperature depth data reveal waters in front of the glaciers at -1.8°C not connected to AIW in the outer fjord, Dove Bugt. The recent grounding line retreat of the glaciers is attributed to glacier thinning at its ablation rate, with little influence of ocean waters, which illustrates the fundamental importance of knowing the bathymetry of glacial fjords.

The International Bathymetric Chart of the Southern Ocean Version 2.

The Southern Ocean surrounding Antarctica is a region that is key to a range of climatic and oceanographic processes with worldwide effects, and is characterised by high biological productivity and biodiversity. Since 2013, the International Bathymetric Chart of the Southern Ocean (IBCSO) has represented the most comprehensive compilation of bathymetry for the Southern Ocean south of 60°S. Recently, the IBCSO Project has combined its efforts with the Nippon Foundation - GEBCO Seabed 2030 Project supporting the goal of mapping the world's oceans by 2030. New datasets initiated a second version of IBCSO (IBCSO v2). This version extends to 50°S (covering approximately 2.4 times the area of seafloor of the previous version) including the gateways of the Antarctic Circumpolar Current and the Antarctic circumpolar frontal systems. Due to increased (multibeam) data coverage, IBCSO v2 significantly improves the overall representation of the Southern Ocean seafloor and resolves many submarine landforms in more detail. This makes IBCSO v2 the most authoritative seafloor map of the area south of 50°S.

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.

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.

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.

The International Bathymetric Chart of the Arctic Ocean Version 4.0.

Bathymetry (seafloor depth), is a critical parameter providing the geospatial context for a multitude of marine scientific studies. Since 1997, the International Bathymetric Chart of the Arctic Ocean (IBCAO) has been the authoritative source of bathymetry for the Arctic Ocean. IBCAO has merged its efforts with the Nippon Foundation-GEBCO-Seabed 2030 Project, with the goal of mapping all of the oceans by 2030. Here we present the latest version (IBCAO Ver. 4.0), with more than twice the resolution (200 × 200 m versus 500 × 500 m) and with individual depth soundings constraining three times more area of the Arctic Ocean (∼19.8% versus 6.7%), than the previous IBCAO Ver. 3.0 released in 2012. Modern multibeam bathymetry comprises ∼14.3% in Ver. 4.0 compared to ∼5.4% in Ver. 3.0. Thus, the new IBCAO Ver. 4.0 has substantially more seafloor morphological information that offers new insights into a range of submarine features and processes; for example, the improved portrayal of Greenland fjords better serves predictive modelling of the fate of the Greenland Ice Sheet.

Earth's water reservoirs in a changing climate.

Progress towards achieving a quantitative understanding of the exchanges of water between Earth's main water reservoirs is reviewed with emphasis on advances accrued from the latest advances in Earth Observation from space. These exchanges of water between the reservoirs are a result of processes that are at the core of important physical Earth-system feedbacks, which fundamentally control the response of Earth's climate to the greenhouse gas forcing it is now experiencing, and are therefore vital to understanding the future evolution of Earth's climate. The changing nature of global mean sea level (GMSL) is the context for discussion of these exchanges. Different sources of satellite observations that are used to quantify ice mass loss and water storage over continents, how water can be tracked to its source using water isotope information and how the waters in different reservoirs influence the fluxes of water between reservoirs are described. The profound influence of Earth's hydrological cycle, including human influences on it, on the rate of GMSL rise is emphasized. The many intricate ways water cycle processes influence water exchanges between reservoirs and thus sea-level rise, including disproportionate influences by the tiniest water reservoirs, are emphasized.

Pathways of ocean heat towards Pine Island and Thwaites grounding lines.

In the Amundsen Sea, modified Circumpolar Deep Water (mCDW) intrudes into ice shelf cavities, causing high ice shelf melting near the ice sheet grounding lines, accelerating ice flow, and controlling the pace of future Antarctic contributions to global sea level. The pathways of mCDW towards grounding lines are crucial as they directly control the heat reaching the ice. A realistic representation of mCDW circulation, however, remains challenging due to the sparsity of in-situ observations and the difficulty of ocean models to reproduce the available observations. In this study, we use an unprecedentedly high-resolution (200 m horizontal and 10 m vertical grid spacing) ocean model that resolves shelf-sea and sub-ice-shelf environments in qualitative agreement with existing observations during austral summer conditions. We demonstrate that the waters reaching the Pine Island and Thwaites grounding lines follow specific, topographically-constrained routes, all passing through a relatively small area located around 104°W and 74.3°S. The temporal and spatial variabilities of ice shelf melt rates are dominantly controlled by the sub-ice shelf ocean current. Our findings highlight the importance of accurate and high-resolution ocean bathymetry and subglacial topography for determining mCDW pathways and ice shelf melt rates.

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.

Origin of Circumpolar Deep Water intruding onto the Amundsen and Bellingshausen Sea continental shelves.

Melting of West Antarctic ice shelves is enhanced by Circumpolar Deep Water (CDW) intruding onto the Amundsen and Bellingshausen Seas (ABS) continental shelves. Despite existing studies of cross-shelf and on-shelf CDW transports, CDW pathways onto the ABS originating from further offshore have never been investigated. Here, we investigate CDW pathways onto the ABS using a regional ocean model. Simulated CDW tracers from a zonal section across 67°S (S04P) circulate along the Antarctic Circumpolar Current (ACC) and Ross Gyre (RG) and travel into ABS continental shelf after 3-5 years, but source locations are shifted westward by ~900 km along S04P in 2001-2006 compared to 2009-2014. We find that simulated on- and off-shelf CDW is ~0.1-0.2 °C warmer in the 2009-2014 case than in the 2001-2006 case together with changes in simulated ocean circulation. These differences are primarily caused by lateral, rather than surface, boundary conditions, implying that large-scale atmospheric and ocean circulations are able to control CDW pathways and thus off- and on-shelf CDW properties.

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.