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.

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.

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.

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.

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.

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.

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.

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.

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).

Undercutting of marine-terminating glaciers in West Greenland.

Marine-terminating glaciers control most of Greenland's ice discharge into the ocean, but little is known about the geometry of their frontal regions. Here we use side-looking, multibeam echo sounding observations to reveal that their frontal ice cliffs are grounded deeper below sea level than previously measured and their ice faces are neither vertical nor smooth but often undercut by the ocean and rough. Deep glacier grounding enables contact with subsurface, warm, salty Atlantic waters (AW) which melts ice at rates of meters per day. We detect cavities undercutting the base of the calving faces at the sites of subglacial water (SGW) discharge predicted by a hydrological model. The observed pattern of undercutting is consistent with numerical simulations of ice melt in which buoyant plumes of SGW transport warm AW to the ice faces. Glacier undercutting likely enhances iceberg calving, impacting ice front stability and, in turn, the glacier mass balance.

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.

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.

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.

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.

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.

Chapter 1. Impacts of the oceans on climate change.

The oceans play a key role in climate regulation especially in part buffering (neutralising) the effects of increasing levels of greenhouse gases in the atmosphere and rising global temperatures. This chapter examines how the regulatory processes performed by the oceans alter as a response to climate change and assesses the extent to which positive feedbacks from the ocean may exacerbate climate change. There is clear evidence for rapid change in the oceans. As the main heat store for the world there has been an accelerating change in sea temperatures over the last few decades, which has contributed to rising sea-level. The oceans are also the main store of carbon dioxide (CO2), and are estimated to have taken up approximately 40% of anthropogenic-sourced CO2 from the atmosphere since the beginning of the industrial revolution. A proportion of the carbon uptake is exported via the four ocean 'carbon pumps' (Solubility, Biological, Continental Shelf and Carbonate Counter) to the deep ocean reservoir. Increases in sea temperature and changing planktonic systems and ocean currents may lead to a reduction in the uptake of CO2 by the ocean; some evidence suggests a suppression of parts of the marine carbon sink is already underway. While the oceans have buffered climate change through the uptake of CO2 produced by fossil fuel burning this has already had an impact on ocean chemistry through ocean acidification and will continue to do so. Feedbacks to climate change from acidification may result from expected impacts on marine organisms (especially corals and calcareous plankton), ecosystems and biogeochemical cycles. The polar regions of the world are showing the most rapid responses to climate change. As a result of a strong ice-ocean influence, small changes in temperature, salinity and ice cover may trigger large and sudden changes in regional climate with potential downstream feedbacks to the climate of the rest of the world. A warming Arctic Ocean may lead to further releases of the potent greenhouse gas methane from hydrates and permafrost. The Southern Ocean plays a critical role in driving, modifying and regulating global climate change via the carbon cycle and through its impact on adjacent Antarctica. The Antarctic Peninsula has shown some of the most rapid rises in atmospheric and oceanic temperature in the world, with an associated retreat of the majority of glaciers. Parts of the West Antarctic ice sheet are deflating rapidly, very likely due to a change in the flux of oceanic heat to the undersides of the floating ice shelves. The final section on modelling feedbacks from the ocean to climate change identifies limitations and priorities for model development and associated observations. Considering the importance of the oceans to climate change and our limited understanding of climate-related ocean processes, our ability to measure the changes that are taking place are conspicuously inadequate. The chapter highlights the need for a comprehensive, adequately funded and globally extensive ocean observing system to be implemented and sustained as a high priority. Unless feedbacks from the oceans to climate change are adequately included in climate change models, it is possible that the mitigation actions needed to stabilise CO2 and limit temperature rise over the next century will be underestimated.

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.