Climate intervention on a high-emissions pathway could delay but not prevent West Antarctic Ice Sheet demise

<jats:title>Abstract</jats:title><jats:p>Solar radiation modification (SRM) is increasingly discussed as a tool to reduce or avert global warming and concomitantly the risk of ice-sheet collapse, as is considered possible for the West Antarctic Ice Sheet (WAIS). Here we analyse the impact of stratospheric aerosol injections on the centennial-to-millennial Antarctic sea-level contribution using an ice-sheet model. We find that mid-twenty-first-century large-scale SRM could delay but ultimately not prevent WAIS collapse in a high-emissions scenario. On intermediate-emissions pathways, SRM could be an effective tool to delay or even prevent an instability of WAIS if deployed by mid-century. However, SRM interventions may be associated with substantial risks, commitments and unintended side effects; therefore, emissions reductions to prevent WAIS collapse seem to be the more practical and sensible approach at the current stage.</jats:p>

section_content: Solar radiation modification (SRM) is increasingly discussed as a tool to reduce or avert global warming and concomitantly the risk of ice-sheet collapse, as is considered possible for the West Antarctic Ice Sheet (WAIS). 

Here we analyse the impact of stratospheric aerosol injections on the centennial-to-millennial Antarctic sea-level contribution using an ice-sheet model. We find that mid-twenty-first-century large-scale SRM could delay but ultimately not prevent WAIS collapse in a high-emissions scenario. On intermediate-emissions pathways, SRM could be an effective tool to delay or even prevent an instability of WAIS if deployed by mid-century. However, SRM interventions may be associated with substantial risks, commitments and unintended side effects; therefore, emissions reductions to prevent WAIS collapse seem to be the more practical and sensible approach at the current stage. 

The notion of offsetting global warming by means of solar climate intervention, or solar radiation modification (SRM), is at the centre of an intense debate on the intrinsic governance dilemma [1] [2] [3]. Stratospheric aerosol deposition, the injection of particles reflecting incoming solar radiation into the stratosphere as a means to limit global warming, is one form of SRM. Substantial SRM interventions promise, in theory, to stop global warming quasi-instantaneously 4, 5 and at relatively limited cost 6, 7. The door to limit climate warming to well below 2 °C via greenhouse gas mitigation within the twenty-first century is rapidly closing 8, 9. As every decade of delayed emissions reduction increases the committed warming 10, SRM interventions might be seriously considered in the future. 

While the effectiveness of SRM to halt anthropogenic climate change has been quantified by modelling studies [11] [12] [13], the regional impacts of planetary SRM remain highly uncertain 14. This is due to an insufficient understanding of the responses of the hydrological cycle 14 and atmospheric chemistry 4 to such interventions. Potential side effects of SRM range from shifting monsoon patterns 15 to changes in the meridional temperature gradients 16, in atmospheric and oceanic circulations 17, 18 and in the statistics of the modes of climate variability 18, all of which are insufficiently studied 3. In summary, SRM could have impacts on regional weather patterns detrimental to society and the biosphere [19] [20] [21], and other still-unknown effects, while not addressing the direct adverse effects of rising atmospheric CO 2 such as ocean acidification 22. Furthermore, the effect of SRM on the polar ice sheets remains largely unexplored. 

SRM research to date has been confined to model simulations and laboratory studies. Field experiments that mimic natural analogues such as volcanic eruptions (for example, the Stratospheric Controlled Perturbation Experiment 23 ) have encountered societal and political barriers to deployment. The physical, governance and ethical challenges of SRM interventions to combat anthropogenic climate change continue to be publicly debated 3, 24. Issues concern the risk of delaying (or even inhibiting) mitigation efforts, unilateral deployment of SRM, termination shock 25 and intergenerational injustice 2. One premise of SRM interventions would be to prevent the triggering of tipping points 26, 27 in the climate system, such as a potential marine ice-sheet instability of the West Antarctic Ice Sheet (WAIS) through self-sustained grounding-line retreat 28 assuming that it is not under way already 29. 

Numerous studies are dedicated to the question whether the observed grounding-line retreat in West Antarctica is due to an ongoing marine ice-sheet instability 30 or simply a response to current https://doi. org/10. 1038/s41558-023-01738-w and four SRM climate scenarios (Fig. 1 ) and extended to the year for this study. We consider four different starting dates of stratospheric aerosol injections under the RCP 8. 5 scenario, either in the year or later in 2040, 2060 and 2080, henceforth called SRM85-20, -40, -60, -80, with equivalent nomenclature used for SRM under the RCP 4. 5 scenario. After the year 2300 (end of HadGEM2-ES simulations), we repeat the last 30 years of each simulation output until the year ce 3000 (an in-depth description of climate forcing and ice-sheet model set-up is provided in Methods).

section_content: In the high-emissions scenario RCP 8. 5, Antarctic surface air temperature anomalies (with respect to the period 1995-2014) peak around 13 °C (Fig. 1a ) in the year 2300. A southward shift of the westerly winds associated with increased ocean upwelling (Extended Data Figs. 2 and 3 ) combined with warming of the subsurface circumantarctic Southern Ocean (adjacent to calving or ice front, averaged over 400-700 m depth) leads to 1 °C and 3. 5 °C warmer ocean waters around the Antarctic continent by the years 2100 and 2300, respectively (Fig. 1b ). Widespread surface melt above the Ross and Filchner-Ronne ice shelves occurs within the second half of the twenty-first century under RCP 8. 5, dominating the overall Antarctic negative surface mass balance (Fig. 1c ). Locally, surface melt occurs on the grounded parts of the ice sheet after 2100 as well. By contrast, neither RCP 4. 5 nor RCP 2. 6 shows sustained surface melt by the year 2300. However, subsurface Southern Ocean warming in RCP 4. 5 is twice as large as in RCP 2. 6 (1. 5 °C versus 0. 7 °C) in 2300, and the global surface air temperature anomaly reaches 3 °C (year 2300). In all scenarios, the subsurface Southern Ocean continues to warm in the year 2300, albeit at a reduced pace (Fig. 1b ).
section_name: Antarctic climate forcing under RCP scenarios

section_content: Instantaneous deployment of a planetary-scale SRM in the year (SRM85-20) leads to an almost immediate (after a few years) levelling of the mean Antarctic surface air temperature anomaly below 2 °C, but the subsurface ocean warming by about 0. 5 °C (similarly to RCP 2. 6) continues until 2050, after which ocean temperatures only slowly increase until the year 2300 (Fig. 1b ). Likewise, in RCP 2. 6, the surface air temperature anomaly does not exceed 2 °C global warming, with a long-term stabilization around 1. 5 °C until the year 2300. Overall, the ocean variability and the resulting ice-shelf mass-balance anomalies 31. Answering this question is challenging due to the intrinsic uncertainties of ice-sheet models 32 and the applied climate forcing 33. It is definite that offsetting sea level would require serious emissions-mitigation efforts 34 ; however, it is unclear whether grounding lines in West Antarctica might continue retreating even under current climate conditions 35, 36. A lack of progress on decarbonizing the economy and society could potentially lead to substantial ice-sheet loss and multi-metre sea-level rise on centennial timescales due to melting of glaciers, the Greenland Ice Sheet and the Antarctic Ice Sheet (AIS) and ocean thermal expansion 37, 38. Considering the severe physical, ecological and societal disruptions caused by sea-level rise, pressure could be mounting to employ some form of SRM on a global, regional or local scale sometime during the twenty-first century. 

However, the literature so far available on the impact of SRM on the polar ice sheets is very limited 39 and restricted mostly to the analysis of atmospheric and oceanic conditions in and around Greenland and Antarctica 17, [40] [41] [42]. To our knowledge, explicit ice dynamic modelling studies under SRM scenarios for Antarctica do not yet exist, while more unconventional climate intervention ideas, such as artificial mass deposition in West Antarctica 43 or building subsea walls blocking warm ocean currents 44, have been investigated. 

In this Article, we investigate the long-term response of the AIS in a geoengineered future compared with a low-emissions strong-mitigation (representative concentration pathway (RCP) 2. 6), an intermediate-emissions stabilization (RCP 4. 5) and a high-emissions no-mitigation (RCP 8. 5) scenario. We are aware that the high-emissions no-mitigation scenario paints a rather pessimistic picture of future GHG emissions as current climate policies are more congruent with the RCP 4. 5 scenario, and emissions-mitigation pledges even leave a chance to stay within 2 °C warming 8. However, it is sensible to assume that RCP 8. 5 is a scenario under which SRM interventions would be seriously considered. 

We employ the Parallel Ice Sheet Model (PISM 45 ) in a range of parameterizations (Extended Data Table 1 ) and at 16 (default configuration), 8 and 4 km grid resolution (Extended Data Fig. 1 ). PISM is forced with previously produced 11 1a, c ). However, intermittent regional ocean temperature difference between RCP 2. 6 and SRM85-20 can be as high as 1. 0-1. 5 °C (Supplementary Fig. 4 ). Importantly, we do not observe increased ocean upwelling and associated changes in temperature in the Southern Ocean under the SRM scenarios, in contrast to a previous study 40 (Supplementary Note 1). While the warming of the subsurface Southern Ocean until the year 2300 in SRM85-60/80 is very similar to that in RCP 4. 5, the surface air temperature in SRM85-60/80 is higher than under RCP 4. 5, resulting in a more negative surface mass balance of the Antarctic ice shelves (Fig. 1c ). It is important to note that HadGEM2-ES ranks among the more-sensitive models in the Coupled Model Intercomparison Project Phase 5 (CMIP5) ensemble with respect to both ocean warming and Antarctic surface mass balance 46. However, the negative surface mass-balance anomalies computed in HadGEM2-ES for the end of the twenty-first century are comparable to recent estimates from a regional climate model 47 in both magnitude and spatial patterns (Extended Data Fig. 4 ; see also Supplementary Note 2 for a discussion of future surface mass-balance changes). We consider this a feature of our forcing as in reality SRM interventions would most likely be considered on a high climate sensitivity trajectory, that is, in a world in which severe climate impacts materialize sooner.
section_name: Persistent Southern Ocean warming despite climate stabilization

section_content: Under RCP 8. 5, the collapse of the WAIS is well under way by the year 2300 (Fig. 2 ). The individual contribution of the WAIS to sea-level rise amounts to 0. 3-0. 8 m (see Fig. 2a ), and that of the entire AIS is 0. 6-1. 1 m (Fig. 2b and Extended Data Table 2 ). This is due to the surface-melt-induced loss of both the Ross and Filchner-Ronne ice shelves and the rapid retreat of Thwaites glacier in response to the ocean thermal forcing in the Amundsen Sea sector. It is important to note that, while ice-sheet retreat and sea-level contribution by the year ce 3000 (WAIS: 2. 1-2. 8; AIS: 2. 9-4. 0 m; Fig. 2a, b and Extended Data Table 2 ) are largely independent of the grid resolutions applied here, the timing of Thwaites collapse is not (Extended Data Fig. 5 ). 

At 16 km horizontal resolution, the inroad to WAIS collapse in RCP 8. 5 begins with the collapse of the Ross Ice Shelf, only then followed by retreat of the Thwaites grounding line. This is in response to the widespread surface melt above the Ross Ice Shelf in the second half of the twenty-first century, leading to the loss of the whole ice shelf and subsequent acceleration of the major Ross-sector ice streams (Fig. 3e, f ). Thus, the collapse of Thwaites is rather a secondary response to the rapid retreat opposite of the WAIS ice divide and ongoing drawdown of the WAIS surface elevation. This suggests a drainage-basin-scale interaction of a marine ice-sheet instability (MISI in one basin triggering MISI in a neighbouring basin) as discussed earlier 48. 

By contrast, at 4 and 8 km resolutions, the collapse of Thwaites Glacier is initiated much earlier as the ice-sheet model response to elevated ocean temperatures is more pronounced. The retreat is initially confined to the Amundsen Sea Embayment and subsequently engulfs the whole WAIS. Here we consider an ice-sheet model ensemble that is consistent with the model sensitivity highlighted in Ice Sheet Model Intercomparison Project for CMIP6 (Supplementary Fig. 5 ). 

In East Antarctica, the Wilkes Basin grounding line along the George V coast remains stable; that is, there is no ongoing forced or self-sustained retreat until the year ce 3000 independent of the emissions scenario (Fig. 3e, f ). This is despite mass loss from this drainage sector. The stability of the George V coast grounding line is robust in all tested resolutions (16, 8 and 4 km) and in agreement with a recent high-resolution study determining the ice sheet's future response after complete loss of ice-shelf buttressing force 49, as well as another study constraining its dynamics during the Last Interglacial 50. However, we note that in the pre-industrial control simulations, ice cover along the George V coast is elevated compared with observations (Supplementary Fig. 8 ). This means that the modelled response might underestimate ice loss in this region. 

Totten Glacier is retreating, and the glaciers of the Recovery Basin are losing mass rapidly (Fig. 3e, f ). However, under RCP 8. 5, mass gains (~0. 2-0. 3 m sea-level-equivalent ice volume) in Dronning Maud Land and central East Antarctica (~0. 4-0. 5 m, respectively) partly compensate for grounding-line retreat in East Antarctica. Interestingly, the substantial grounding-line retreat of the Recovery, Bailey and Slessor glaciers is balanced by mass gains in the hinterland of these glaciers. 

In RCP 8. 5, Antarctica contributes 3. 2-4. 0 m to global sea level by the year ce 3000 (Fig. 2b and Extended Data Table 2 ) and still loses mass at a rate 5-10 times faster than observed today (Fig. 2d and Supplementary Fig. 6 ). According to our simulations, the impact of runaway climate change on Antarctic ice-mass change will be discernible from mitigation or SRM scenarios only by the second half of the twenty-first century. Its imprint on global sea level diverges even later (Fig. 2a, b ), in agreement with recent model estimates 51, 52 but contrasting studies where hydrofracturing and cliff failure are considered 53. We find that for both RCP 8. 5 and RCP 4. 5, WAIS collapse is irreversible; that is, there is no recovery on millennial timescales (Supplementary Note 3 and Supplementary Figs. 2 and 3 ).
section_name: Antarctic ice loss under RCP 8.5

section_content: The only scenarios mostly ensuring WAIS stability (absence of MISI or ongoing grounding-line retreat) are RCP 2. 6 (76% of ensemble members) and early SRM interventions, meaning planetary SRM intervention starting in either the year 2020 (RCP 8. 5; Fig. 2a, e ) or the year 2040 (RCP 4. 5; Fig. 4b ). This is not surprising as the climatic drivers in RCP 2. 6 and SRM85-20 are very similar. However, despite the strong emissions reductions in RCP 2. 6 or the stabilization of global temperatures in SRM85-20, the subsurface circumantarctic ocean has a committed warming of an additional 0. 5-0. 7 °C by 2300 (regionally, this warming can reach 1. 0-1. 5 °C; Supplementary Fig. 4 ). 

In the few model runs in which Thwaites Glacier is already close to or beyond a tipping point today, that is, a configuration in which a self-sustained unforced retreat of the grounding line ensues, this warming can be sufficient to trigger WAIS collapse under several parameter sets considered here. This explains the high-end tail of AIS sea-level contributions in RCP 2. 6 (1. 38 m) and SRM85-20 (1. 08 m) by the year ce 3000 (Fig. 2a, e ). Ocean temperatures in the Amundsen, Weddel and Ross seas are up to 0. 5, 1. 0 and 1. 5 °C warmer during the middle twenty-first century (Supplementary Fig. 4 ), which could be one factor behind the higher sea-level range in RCP 2. 6. 

WAIS collapse in lower-emissions/early SRM scenarios unfolds via the disintegration of Thwaites Glacier associated with accelerating grounding-line retreat in the coming centuries. Both the Filchner-Ronne and Ross ice shelves remain intact during the next centuries in the SRM85-20, RCP 2. 6 and RCP 4. 5 scenarios; thus, the only path towards WAIS collapse is due to elevated basal melt underneath the ice shelves of the Amundsen Sea sector. Importantly, the chances of WAIS collapse in SRM85-20 and RCP 2. 6 are nearly identical, with very similar median long-term sea-level contributions (~0. 50 and 0. 75 m by the year ce 3000, respectively). Notably, while some simulations in the ensemble do project negative sea-level contributions for RCP 2. 6 and SRM85-20, those matching the historical record do not.
section_name: Ice-sheet stability under rapid decarbonization or early SRM

section_content: Current decarbonization policies and actions correspond to a middle-of-the-road scenario resulting in warming of global temperatures by ~2. 8 °C by the end of this century 54. We investigate the impact of a continuous SRM intervention under the RCP 4. 5 scenario following the same approach as in the case of RCP 8. 5 (employment in 2020, 2040,. Black line depicts present-day observed grounding line 59, 60. Bathymetry/topography plotted for reference in the background 59, 60.
section_name: SRM in a middle-of-the-road scenario

section_content: https://doi. org/10. 1038/s41558-023-01738-w 2060 and 2080; in the case of SRM4. 5, SO 2 emissions are two-thirds lower compared with SRM8. 5 by the end of the twenty-first century as shown in Supplementary Fig. 7 ). Surprisingly, the committed ice loss until the year ce 3000 is very similar to the high-emissions pathway in case of a late intervention (2080; red dashed line in Fig. 4b, d, f ). However, most ensemble members (65%) suggest a stabilization of the grounding line under SRM45-40 (SRM started in 2040), whereas a long-term collapse is triggered for SRM85-40 in 70% of the simulations. This indicates that WAIS collapse might be prevented by SRM under an RCP 4. 5 scenario even for mid-century interventions, however with a substantial risk of failure (35%).
section_name: Article

section_content: Delay of SRM85 by 20, 40 or 60 years increases both the committed Southern Ocean warming (Fig. 1b ) and surface melt above the Ross Ice Shelf and therefore the probability of WAIS collapse (70%, 76%, 82%; Fig. 4 ). Maximum (median) sea-level contributions in these scenarios amount to 1. 6 (0. 9) m, 1. 9 (1. 2) m and 2. 2 (1. 4) m, respectively. Moreover, if SRM is started in the year 2060, widespread surface melt over the Ross Ice Shelf is already present (Fig. 5 ), resulting in collapse of the ice shelf and loss of buttressing. If WAIS collapse is already triggered, the SRM intervention would have to reverse, not merely halt, climate warming due to the nonlinear increase in ice loss and hysteresis effects of the ice sheet 55 (Supplementary Figs. 1 and 2 ). Already in RCP 4. 5, thus under current climate policies, the chance of long-term WAIS stability is reduced to a third of the ensemble members, and employment of SRM in the second half of this century mostly fails to prevent the long-term demise of the WAIS.
section_name: Committed ice-sheet collapse despite climate intervention

section_content: Our simulations suggest that the most effective path to prevent total long-term WAIS collapse is rapid decarbonization or, depending on GHG emissions, the early deployment of planetary SRM. The latter is politically challenging and associated with unintended and potentially damaging side effects. We believe that a robust understanding of the regional impacts of planetary SRM, or any kind of climate intervention, is imperative before considering large-scale employment. Thus, at the current stage, we suggest that the best chance to guarantee the long-term stability of the WAIS is to reduce GHG emissions to net zero with no delay, thereby keeping the goal of the Paris Agreement 56 to limit global warming to well below 2 °C in reach. While in the worst-case emissions scenario (RCP 8. 5) WAIS collapse and East Antarctic Ice Sheet retreat are unavoidable according to our model simulations, even RCP 4. 5 would lead to a high chance of a long-term demise of WAIS. This committed retreat could be slowed should SRM interventions as considered here materialize during the twenty-first century, providing a larger time window for adaptation measures. We use a stand-alone ice-sheet model; interactive coupling to a climate model or a different model initialization would affect the timing of ice-sheet tipping and pacing of retreat. Refer to Methods and Supplementary Information for an in-depth discussion of caveats associated with our modelling choices such as parameterization, model spin-up, climate forcing and mismatches with respect to present-day observations (for example, Supplementary Notes 4 and 5). 

It is important to note that we have not considered the possibility of reversing temperature conditions under high-emissions scenarios 4. Reducing Antarctic surface and subsequently ocean temperatures via SRM after surpassing a certain climate threshold (2 °C) during the twenty-first century could further delay the demise of the WAIS. Limited polar SRM interventions 57 combined with strong emissions-mitigation policies that aim at preventing the tipping of the WAIS would be an additional modelling scenario. However, this requires more elaborate model set-ups and preferably interactive coupling between the ice sheet and the climate system. However, as long as atmospheric GHG concentrations remain high, SRM would most likely have to be continued for centuries 58, perhaps even millennia, to avoid WAIS collapse. This would be an intergenerational commitment whose cost and risks, for example due to permanent changes in the lower stratosphere and the danger of termination shock, are potentially immense but difficult to quantify. Our simulations show, in agreement with previous studies, that the impact of the adopted emissions pathways on the future AIS mass imbalance will not be discernible within this century. However, the longer actions towards limiting global warming are delayed, be it via emissions mitigation and/or SRM, the more likely a collapse of WAIS will become.
section_name: Discussion

section_content: Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi. org/10. 1038/s41558-023-01738-w.
section_name: Online content

section_content: https://doi. org/10. 1038/s41558-023-01738-w Open Access This article is licensed under a Creative Commons Attribution 4. 0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/ licenses/by/4. 0/. © The Author(s) 2023 https://doi. org/10. 1038/s41558-023-01738-w
section_name: Article

section_name: Methods

section_content: To investigate the response of the AIS to a wide range of future emissions and SRM scenarios, we employ the PISM 45 with a range of parameterizations controlling ice flow and ice-shelf mass balance. Details of the spin-up of PISM to the year 1860 are given in the following. The spin-up is followed by a 145 yr common historical forcing period after which the individual ensembles are branched off in 2006 and integrated until the year ce 3000. We constrain the model ensemble by focusing on those runs that best match the recent historical record (1992-2017) of sea-level-equivalent ice volume, mass balance and annual sea-level contribution, as well as ice thickness and grounding-line positions. We employ PISM (v. 1. 2) in the so-called shallow shelf-shallow ice hybrid mode, which is a superposition of the shallow shelf and shallow ice approximations of the stress balance. This allows for model representation of both slow non-sliding (inland ice) and fast-flowing (glaciers) parts of the ice. To allow for 540 simulations over millennial timescales, we use a relatively coarse horizontal resolution of 16 km. The position of the grounding line is established by the flotation criterion employing a subgrid interpolation scheme 62. Basal friction and stresses are computed accordingly. This enables a more adequate representation of grounding-line dynamics at coarse resolution. However, we repeat selected simulations on 8 and 4 km grids to assess the robustness of the respective coarse-resolution results (Extended Data Fig. 5 ). We note that the computation of basal melt rate at the grounding line is not affected by the subgrid interpolation but happens on the native grid. Recently, it has been suggested that seawater intrusions at the grounding line can increase melt rates locally 63. A subgrid extrapolation of basal melt at the grounding line is one way to mimic this effect. However, it is not yet possible to quantify the increase in melt rates in situ, and a subgrid melt extrapolation almost certainly would overestimate local melt rates and bias the model to be more prone to grounding-line retreat. 

In the model ensemble, we vary two parameters relevant for ice flow (Extended Data Table 1 ). In PISM, the till friction angle controls the yield strength in the Mohr-Coulomb criterion with low values of the till friction angle resulting in low yield strength. In our simulations, the till friction angle is linearly dependent on depth, with a maximum value of 30 and a minimum value between 3 and 8. For sliding, we employ a power law ('pseudo-plastic') with an exponent q of 0. 60 or 0. 75. The preceding parameters are chosen as they create a near-equilibrium ice sheet with reasonable grounding-line positions and ice cover under pre-industrial climate forcing (Supplementary Figs. 8 and 9 ). For a full list of parameters, see Extended Data Table 1. 

Due to a lack of direct observations, basal ice-shelf melt, the major driver behind current grounding-line retreat, is highly uncertain. Therefore, we employ eight different basal ice-shelf-melt sensitivities using the Potsdam Ice-Shelf Cavity mOdel (PICO) basal shelf melt model 64, 65. 

Here we vary the heat exchange coefficient that controls the turbulent heat exchange from the ambient ocean across the boundary layer beneath the ice-shelf base (Extended Data Table 1 ). Accounting for different parameterizations of ice-sheet flow and basal friction, this amounts to 60 simulations per climate forcing and allows us to consider a wide range of potential ice-sheet responses. The sensitivity of ice-shelf melt to the heat exchange coefficient is shown in Supplementary Fig. 10. Recent observations for the Thwaites Eastern Ice Shelf 66 have shown that subshelf melt rates suggested by current regional ocean models or simplified parameterizations can considerably deviate from observations (as far as they are available). While we try to match currently estimated melt rate patterns adequately, substantial deviations from reality are to be expected. 

Calving amounts to about half of the currently estimated ice-shelf mass loss 67, 68. Due to lack of a general physical fracture and melt-induced calving law, we employ a heuristic calving relationship that depends on the horizontal strain-rate field ('eigencalving' 69 ). We remove all floating ice thinner than 10 m to prevent unrealistic ice-shelf geometries. Finally, we prescribe the current ice-shelf geometry as the maximum extent of floating ice to exclude unrealistic buttressing effects not in line with observations.
section_name: Ice-sheet model parameterizations, set-up and resolution

section_content: We carry out a separate control simulation under constant pre-industrial forcing for each simulation and apply additional ensemble constraints on the basis of observational records. We constrain the whole model ensemble by calibration against observations of the Antarctic sea-level contribution 70, the ice-shelf mass loss due to calving and basal melting 67, 68 and the present-day observed ice-mass change and grounding-line positions 61 (Supplementary Fig. 11 ).
section_name: Historical observational ensemble constraints

section_content: The model spin-up is divided into three steps with the goal to relax the ice sheet to present-day conditions. In a first step, the thermal state of the ice sheet is initiated via a 200,000 yr simulation under constant present-day forcing with a fixed ice-sheet geometry. The present-day surface climate forcing is taken from the regional climate model RACMO 71. Subsequently, we run a 2,000 yr pre-industrial spin-up with pre-industrial climate anomalies (100 yr repeated) from HadGEM2-ES in which the ice sheet is allowed to freely evolve. These anomalies are added onto the RACMO surface mass balance and surface air temperature data and onto the World Ocean Atlas 2018 temperature and salinity fields as used in the Ice Sheet Model Intercomparison Project for CMIP6 34. Ocean conditions are averaged between depths of 700 and 400 m, representative of the approximate depth of Antarctic ice-shelf cavities. The ensemble spread of the pre-industrial model spin-up is shown in Supplementary Fig. 9. In a last step, the model is run from 1860 until 2005 under common historical forcing from HadGEM2-ES, after which the different climate scenarios are branched off. The historical runs are restarted from a best-fit member of the pre-industrial spin-up ensemble (Supplementary Fig. 8 ). Note that the higher-resolution simulations (8 km and 4 km) are restarted and regridded from the 16 km pre-industrial spin-up.
section_name: Model spin-up

section_content: We consider four different starting dates for SRM under RCP 8. 5 (2020, 2040, 2060 and 2080). While SRM85-20 is modelled explicitly in HadGEM2-ES, we construct the forcing for later starting dates by adding the post-2020 modelled climate anomalies of SRM85-20 onto RCP 8. 5 at the years 2040, 2060 and 2080, respectively. This means we follow the RCP 8. 5 scenario until the respective starting dates, upon which we continue the forcing, adding the SRM85-20 anomalies to the RCP 8. 5 climate state at the year 2040, 2060 or 2080. This implies that after the start of SRM, the climate system responds similarly, regardless of whether SRM is started in the year 2020 or later during the twenty-first century. We consider this to be a reasonable assumption as the modelled twenty-first-century temperature response to RCP 8. 5 remains approximately linear throughout.
section_name: Start of SRM

section_content: We force PISM with anomalies (relative to pre-industrial) annual mean surface mass balance, surface air temperature, ocean temperature and salinity fields regridded from the native HadGEM2-ES-grid to the PISM-grid at 16, 8 or 4 km lateral resolution. The anomalies are added to the present-day reference field from RACMO 71 and World Ocean Data 2018 72. The HadGEM2-ES climate projections cover the period from 1860 to 2300. To continue the climate forcing until the year ce 3000, we repeat the last 30 years (2270-2300) of the forcing until the year ce 3000. This means that we ignore natural or anthropogenic changes in atmospheric GHG concentrations and their temperature response after the year 2300. https://doi. org/10. 1038/s41558-023-01738-w
section_name: Model forcing

section_content: HadGEM2-ES is a CMIP5-generation Earth system model (ESM) with reduced complexity compared with CMIP6-generation ESMs such as its successor, UKESM1 73. HadGEM2-ES is a fully coupled atmosphere-ocean ESM with an atmospheric horizontal resolution of N96 (1. 875° × 1. 25°) and 38 vertical levels extending to 40 km altitude and an oceanic horizontal resolution of 1° (reducing to 1/3° at the Equator) with 40 vertical levels 11. HadGEM2-ES employs the single-moment Coupled Large-scale Aerosol Simulator for Studies in Climate (CLASSIC) aerosol module, which includes a fully interactive sulfur cycle that represents sulfate in Aitken, accumulation and dissolved modes 74. HadGEM2-ES has been utilized for many SRM studies independently and as part of GeoMIP 11, 12, 74, and the CLASSIC aerosol module remains in operational use for daily air quality forecasts in the United Kingdom 75. A full description of HadGEM2-ES and its contribution to CMIP5 is provided in refs. 46,73. The design of the SRM simulations is to inject sulfur dioxide (SO 2 ) into the model stratosphere at such a rate as to maintain global-mean near-surface air temperature at 1. 5 °C above pre-industrial levels, where the reference temperature is derived from HadGEM2-ES's 140 yr pre-industrial control simulation 15. Three baseline GHG concentration scenarios from the CMIP5 protocol were chosen to span the range of plausible futures: the low-emissions strong-mitigation RCP 2. 6, the intermediate-emissions stabilization RCP 4. 5 and the high-emissions no-mitigation RCP 8. 5 76, 77. Ref. 15 provides detail on how the temporally evolving SO 2 emission rates were determined for each SRM scenario, which pre-empts the online feedback algorithm used in GLENS (Geoengineering Large Ensemble Project) 8. For the purpose of this study, the original GEO2. 6 (SRM26), GEO4. 5 (SRM45) and GEO8. 5 (SRM85) SRM simulations were extended to the year 2300 following the RCP extension protocol outlined in ref. 71 and using the same algorithm for calculating SO 2 injection rates as described in ref. 11.
section_name: Earth system model configuration and simulation design

section_content: Model sensitivity. Overall, we consider our ice-sheet model set-up to have a relatively low sensitivity to climate change compared with previous studies with PISM (for example, ref. 78) and to other models currently in use (Supplementary Fig. 5 ) 34. Consequently, we conclude that the modelled effectiveness (preventing ice-sheet instability) of SRM interventions shown here represents a conservative estimate with respect to ice-sheet stability (more-sensitive models would show less effectiveness of SRM). This, however, underlines our assessment that large-scale SRM interventions during the twenty-first century, if not employed early, are not sufficient to prevent WAIS collapse, assuming intermediate-or high-emissions scenarios. 

Stand-alone ice-sheet model. As we run PISM in stand-alone mode (not coupled to a climate model), the climate forcing does not include feedback from ice-sheet topographic changes and freshwater hosing from ice-shelf melt/calving and surface run-off. The WAIS undergoes the most drastic elevation changes (ice-sheet collapse) until the year ce 3000; therefore, both freshwater input into the Southern Ocean and surface temperature lapse rate corrections would be largest in this area. A recent study in which ice-ocean-atmosphere feedbacks were incorporated via offline coupling of PISM to a climate model 78 suggests that the combined feedbacks accelerate ice-sheet loss in Antarctica until the year 2100. However, also in East Antarctica, where grounding-line retreat is mostly contained, surface lowering and increased ice discharge take place that might affect the magnitude of ice-sheet retreat in a coupled model set-up. In summary, it is difficult to quantify the extent to which the missing ice-to-ocean and ice-to-atmosphere feedbacks would affect ice-sheet stability. 

Model resolution. We run the ice-sheet model on a relatively coarse resolution of 16 km to allow for a larger ensemble over millennial timescales. While the subgrid grounding-line interpolation is supposed to make the results more independent from model resolution, we still find different ice-sheet responses at higher resolutions of 8 and 4 km (Extended Data Figs. 4 and 5 ). The 8 km ice-sheet configuration at the year ce 3000 is very similar to the 16 km runs, but the pacing of retreat (onset of marine ice-sheet instability) is different. At 8 and 4 km grid resolutions, Thwaites Glacier shows a nonlinear retreat earlier compared with the 16 km simulations. However, the George V Coast grounding line bordering the Wilkes Subglacial Basin remains stable (until the year ce 3000) regardless of the resolution tested here. Taken together, this also suggests that SRM interventions would be less effective at higher resolutions, further strengthening our line of argument. 

Initialization. Another issue pertaining to all ice-sheet modelling efforts is the initial configuration of the ice sheet at the beginning of the simulations. As we do not employ inversion of, for example, basal friction but rely on more general heuristics, ice-sheet thickness and grounding-line positions in critical areas are not identical to present-day observations. For example, the grounding line of Thwaites Glacier is located slightly farther upstream and WAIS surface elevation is lower than presently observed, which could favour faster retreat. Furthermore, some East Antarctic coastal regions are thicker than observations and hence more resilient to ocean warming. 

Model bias in ice-shelf tributaries and surface velocities. See Supplementary Information.
section_name: Caveats of study design

section_name: Acknowledgements

section_content: Calculations were performed on UBELIX, the high-performance computing cluster at the University of Bern. We thank G. Jansen for excellent technical support on the HPC infrastructure. J. S. acknowledges funding from the Deutsche Forschungsgemeinschaft under grant no. SU 1166/1-1. T. F. S. acknowledges funding from the Swiss National Science Foundation (grant 200020_200492 ) and from the European Union 's Horizon 2020 research and innovation programme under grant agreement no. 820970 (project TiPES ; this is TiPES contribution no. 199 ). T. L. F. acknowledges funding from the Swiss National Science Foundation (grant PP00P2_198897 ) and the European Union 's Horizon 2020 research and innovation programme under grant agreement no. 820989 (project COMFORT ).
section_name: Acknowledgements

section_name: Data availability

section_content: The data generated and analysed during the current study are available from the corresponding author. The data underlying the main figures are uploaded to Zenodo (https://doi. org/10. 5281/zenodo. 7640508) (ref. 79).
section_name: Data availability

section_content: PISM is openly available from https://github. com/pism/pism.
section_name: Code availability

section_name: Author contributions

section_content: J. S., T. F. S. and T. L. F. devised the study. J. S. set up the ice-sheet model and ran the ensemble simulations. A. J. carried out the prolonged solar radiation modification experiments. C. W. carried out the reversibility tests. J. S. analysed the ice-sheet model results T. L. F. carried out analysis on the ocean component of HadGEM2-ES. J. S. led the writing of the paper with contributions from all authors.
section_name: Author contributions

section_content: The authors declare no competing interests.
section_name: Competing interests

section_content: Extended data is available for this paper at https://doi. org/10. 1038/s41558-023-01738-w.
section_name: Additional information

section_content: The online version contains supplementary material available at https://doi. org/10. 1038/s41558-023-01738-w. 

Correspondence and requests for materials should be addressed to J. Sutter.
section_name: Supplementary information

section_content: Feldmann and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. 

Reprints and permissions information is available at www. nature. com/reprints.
section_name: Peer review information Nature Climate Change thanks Johannes