Sea level rise (SLR) may impose substantial economic costs to coastal communities worldwide, but characterizing its global impact remains challenging because SLR costs depend heavily on natural characteristics and human investments at each location – including topography, the spatial distribution of assets, and local adaptation decisions. To date, several impact models have been developed to estimate the global costs of SLR. Yet, the limited availability of open-source and modular platforms that easily ingest up-to-date socioeconomic and physical data sources restricts the ability of existing systems to incorporate new insights transparently. In this paper, we present a modular, open-source platform designed to address this need, providing end-to-end transparency from global input data to a scalable least-cost optimization framework that estimates adaptation and net SLR costs for nearly 10 000 global coastline segments and administrative regions. Our approach accounts both for uncertainty in the magnitude of global mean sea level (g.m.s.l.) rise and spatial variability in local relative sea level rise. Using this platform, we evaluate costs across 230 possible socioeconomic and SLR trajectories in the 21st century. According to the latest Intergovernmental Panel on Climate Change Assessment Report (AR6), g.m.s.l. is likely to rise during the 21st century by 0.40–0.69 m if late-century warming reaches 2 ∘C and by 0.58–0.91 m with 4 ∘C of warming (Fox-Kemper et al., 2021). With no forward-looking adaptation, we estimate that annual costs of sea level rise associated with a 2 ∘C scenario will likely fall between USD 1.2 and 4.0 trillion (0.1 % and 1.2 % of GDP, respectively) by 2100, depending on socioeconomic and sea level rise trajectories. Cost-effective, proactive adaptation would provide substantial benefits, lowering these values to between USD 110 and USD 530 billion (0.02 and 0.06 %) under an optimal adaptation scenario. For the likely SLR trajectories associated with 4 ∘C warming, these costs range from USD 3.1 to 6.9 trillion (0.3 % and 2.0 %) with no forward-looking adaptation and USD 200 billion to USD 750 billion (0.04 % to 0.09 %) under optimal adaptation. The Intergovernmental Panel on Climate Change (IPCC) notes that deeply uncertain physical processes like marine ice cliff instability could drive substantially higher global sea level rise, potentially approaching 2.0 m by 2100 in very high emission scenarios. Accordingly, we also model the impacts of 1.5 and 2.0 m g.m.s.l. rises by 2100; the associated annual cost estimates range from USD 11.2 to 30.6 trillion (1.2 % and 7.6 %) under no forward-looking adaptation and USD 420 billion to 1.5 trillion (0.08 % to 0.20 %) under optimal adaptation. Our modeling platform used to generate these estimates is publicly available in an effort to spur research collaboration and support decision-making, with segment-level physical and socioeconomic input characteristics provided at https://doi.org/10.5281/zenodo.7693868 (Bolliger et al., 2023a) and model results at https://doi.org/10.5281/zenodo.7693869 (Bolliger et al., 2023b).
more »
« less
Ice sheet contributions to future sea-level rise from structured expert judgment
Despite considerable advances in process understanding, numerical modeling, and the observational record of ice sheet contributions to global mean sea-level rise (SLR) since the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change, severe limitations remain in the predictive capability of ice sheet models. As a consequence, the potential contributions of ice sheets remain the largest source of uncertainty in projecting future SLR. Here, we report the findings of a structured expert judgement study, using unique techniques for modeling correlations between inter- and intra-ice sheet processes and their tail dependences. We find that since the AR5, expert uncertainty has grown, in particular because of uncertain ice dynamic effects. For a +2 °C temperature scenario consistent with the Paris Agreement, we obtain a median estimate of a 26 cm SLR contribution by 2100, with a 95th percentile value of 81 cm. For a +5 °C temperature scenario more consistent with unchecked emissions growth, the corresponding values are 51 and 178 cm, respectively. Inclusion of thermal expansion and glacier contributions results in a global total SLR estimate that exceeds 2 m at the 95th percentile. Our findings support the use of scenarios of 21st century global total SLR exceeding 2 m for planning purposes. Beyond 2100, uncertainty and projected SLR increase rapidly. The 95th percentile ice sheet contribution by 2200, for the +5 °C scenario, is 7.5 m as a result of instabilities coming into play in both West and East Antarctica. Introducing process correlations and tail dependences increases estimates by roughly 15%.
more »
« less
- Award ID(s):
- 1663807
- NSF-PAR ID:
- 10105828
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 116
- Issue:
- 23
- ISSN:
- 0027-8424
- Page Range / eLocation ID:
- 11195 to 11200
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract The ice sheets covering Antarctica and Greenland present the greatest uncertainty in, and largest potential contribution to, future sea level rise. The uncertainty arises from a paucity of suitable observations covering the full range of ice sheet behaviors, incomplete understanding of the influences of diverse processes, and limitations in defining key boundary conditions for the numerical models. To investigate the impact of these uncertainties on ice sheet projections we undertook a structured expert judgement study. Here, we interrogate the findings of that study to identify the dominant drivers of uncertainty in projections and their relative importance as a function of ice sheet and time. We find that for the 21st century, Greenland surface melting, in particular the role of surface albedo effects, and West Antarctic ice dynamics, specifically the role of ice shelf buttressing, dominate the uncertainty. The importance of these effects holds under both a high‐end 5°C global warming scenario and another that limits global warming to 2°C. During the 22nd century the dominant drivers of uncertainty shift. Under the 5°C scenario, East Antarctic ice dynamics dominate the uncertainty in projections, driven by the possible role of ice flow instabilities. These dynamic effects only become dominant, however, for a temperature scenario above the Paris Agreement 2°C target and beyond 2100. Our findings identify key processes and factors that need to be addressed in future modeling and observational studies in order to reduce uncertainties in ice sheet projections.
-
Observational estimates of Antarctic ice loss have accelerated in recent decades, and worst-case scenarios of modeling studies have suggested potentially catastrophic sea level rise (~2 meters) by the end of the century. However, modeled contributions to global mean sea level from the Antarctic ice-sheet (AIS) in the 21st century are highly uncertain, in part because ice-sheet model parameters are poorly constrained. Individual ice-sheet model runs are also deterministic and not computationally efficient enough to generate the continuous probability distributions required for incorporation into a holistic framework of probabilistic sea-level projections. To address these shortfalls, we statistically emulate an ice-sheet model using Gaussian Process (GP) regression. GP modeling is a non-parametric machine-learning technique which maps inputs (e.g. forcing or model parameters) to target outputs (e.g. sea-level contributions from the Antarctic ice-sheet) and has the inherent and important advantage that emulator uncertainty is explicitly quantified. We construct emulators for the last interglacial period and an RCP8.5 scenario, and separately for the western, eastern, and total AIS. Separate emulation of western and eastern AIS is important because their evolutions and physical responses to climate forcing are distinct. The emulators are trained on 196 ensemble members for each scenario, composed by varying the parameters of maximum rate of ice-cliff wastage and the coefficient of hydrofracturing. We condition the emulators on last interglacial proxy sea-level records and modern GRACE measurements and exclude poor-fitting ensemble members. The resulting emulators are sampled to produce probability distributions that fill intermediate gaps between discrete ice-sheet model outcomes. We invert emulated high and low probability sea-level contributions in 2100 to explore 21st century evolution pathways; results highlight the deep uncertainty of ice-sheet model physics and the importance of using observations to narrow the range of parameters. Our approach is designed to be flexible such that other ice-sheet models or parameter spaces may be substituted and explored with the emulator.more » « less
-
more » « less
Abstract Sea level rise (SLR) is a long‐lasting consequence of climate change because global anthropogenic warming takes centuries to millennia to equilibrate for the deep ocean and ice sheets. SLR projections based on climate models support policy analysis, risk assessment and adaptation planning today, despite their large uncertainties. The central range of the SLR distribution is estimated by process‐based models. However, risk‐averse practitioners often require information about plausible future conditions that lie in the tails of the SLR distribution, which are poorly defined by existing models. Here, a community effort combining scientists and practitioners builds on a framework of discussing physical evidence to quantify high‐end global SLR for practitioners. The approach is complementary to the IPCC AR6 report and provides further physically plausible high‐end scenarios. High‐end estimates for the different SLR components are developed for two climate scenarios at two timescales. For global warming of +2°C in 2100 (RCP2.6/SSP1‐2.6) relative to pre‐industrial values our high‐end global SLR estimates are up to 0.9 m in 2100 and 2.5 m in 2300. Similarly, for a (RCP8.5/SSP5‐8.5), we estimate up to 1.6 m in 2100 and up to 10.4 m in 2300. The large and growing differences between the scenarios beyond 2100 emphasize the long‐term benefits of mitigation. However, even a modest 2°C warming may cause multi‐meter SLR on centennial time scales with profound consequences for coastal areas. Earlier high‐end assessments focused on instability mechanisms in Antarctica, while here we emphasize the importance of the timing of ice shelf collapse around Antarctica. This is highly uncertain due to low understanding of the driving processes. Hence both process understanding and emission scenario control high‐end SLR.
-
Abstract. Multi-meter sea level rise (SLR) is thought to be possible within the next few centuries, with most of the uncertainty originating from the Antarctic land ice contribution. One source of uncertainty relates to the ice sheet model initialization. Since ice sheets have a long response time (compared to other Earth system components such as the atmosphere), ice sheet model initialization methods can have significant impacts on how the ice sheet responds to future forcings. To assess this, we generated 25 different ice sheet spin-ups, using the Community Ice Sheet Model (CISM) at a 4 km resolution. During each spin-up, we varied two key parameters known to impact the sensitivity of the ice sheet to future forcing: one related to the sensitivity of the ice shelf melt rate to ocean thermal forcing (TF) and the other related to the basal friction. The spin-ups all nudge toward observed thickness and enforce a no-advance calving criterion, such that all final spin-up states resemble observations but differ in their melt and friction parameter settings. Each spin-up was then forced with future ocean thermal forcings from 13 different CMIP6 models under the Shared Socioeconomic Pathway (SSP)5-8.5 emissions scenarioand modern climatological surface mass balance data. Our results show that the effects of the ice sheet and ocean parameter settings used during the spin-up are capable of impacting multi-century future SLR predictions by as much as 2 m. By the end of this century, the effects of these choices are more modest, but still significant, with differences of up to 0.2 m of SLR. We have identified a combined ocean and ice parameter space that leads to widespread mass loss within 500 years (low friction and high melt rate sensitivity). To explore temperature thresholds, we also ran a synthetically forced CISM ensemble that is focused on the Amundsen region only. Given certain ocean and ice parameter choices, Amundsen mass loss can be triggered with thermal forcing anomalies between 1.5 and 2 ∘C relative to the spin-up.Our results emphasize the critical importance of considering ice sheet and ocean parameter choices during spin-up for SLR predictions and suggest the importance of including glacial isostatic adjustment in ice sheet simulations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.1817205116
