Gaussian Process emulation of multi-model ice sheet and glacier projections

Author(s):  
Tamsin Edwards ◽  

<p><strong>The land ice contribution to global mean sea level rise has not yet been predicted with ice sheet and glacier models for the latest set of socio-economic scenarios (SSPs), nor with coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects (ISMIP6 and GlacierMIP) generated a large suite of projections using multiple models, but mostly used previous generation scenarios and climate models, and could not fully explore known uncertainties. </strong></p><p><strong>Here we estimate probability distributions for these projections for the SSPs using Gaussian Process emulation of the ice sheet and glacier model ensembles. We model the sea level contribution as a function of global mean surface air temperature forcing and (for the ice sheets) model parameters, with the 'nugget' allowing for multi-model structural uncertainty. Approximate independence of ice sheet and glacier models is assumed, because a given model responds very differently under different setups (such as initialisation). </strong></p><p><strong>We find that limiting global warming to 1.5</strong>°<strong>C </strong><strong>would halve the land ice contribution to 21<sup>st</sup> century </strong><strong>sea level rise</strong><strong>, relative to current emissions pledges: t</strong><strong>he median decreases from 25 to 13 cm sea level equivalent (SLE) by 2100. However, the Antarctic contribution does not show a clear response to emissions scenario, due to competing processes of increasing ice loss and snowfall accumulation in a warming climate. </strong></p><p><strong>However, under risk-averse (pessimistic) assumptions for climate and Antarctic ice sheet model selection and ice sheet model parameter values, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 cm SLE under current policies and pledges, with the 95<sup>th</sup> percentile exceeding half a metre even under 1.5</strong>°<strong>C warming. </strong></p><p><strong>Gaussian Process emulation can therefore be a powerful tool for estimating probability density functions from multi-model ensembles and testing the sensitivity of the results to assumptions.</strong></p>

2013 ◽  
Vol 26 (8) ◽  
pp. 2502-2513 ◽  
Author(s):  
N. Bouttes ◽  
J. M. Gregory ◽  
J. A. Lowe

Abstract During the last century, global climate has been warming, and projections indicate that such a warming is likely to continue over coming decades. Most of the extra heat is stored in the ocean, resulting in thermal expansion of seawater and global mean sea level rise. Previous studies have shown that after CO2 emissions cease or CO2 concentration is stabilized, global mean surface air temperature stabilizes or decreases slowly, but sea level continues to rise. Using idealized CO2 scenario simulations with a hierarchy of models including an AOGCM and a step-response model, the authors show how the evolution of thermal expansion can be interpreted in terms of the climate energy balance and the vertical profile of ocean warming. Whereas surface temperature depends on cumulative CO2 emissions, sea level rise due to thermal expansion depends on the time profile of emissions. Sea level rise is smaller for later emissions, implying that targets to limit sea level rise would need to refer to the rate of emissions, not only to the time integral. Thermal expansion is in principle reversible, but to halt or reverse it quickly requires the radiative forcing to be reduced substantially, which is possible on centennial time scales only by geoengineering. If it could be done, the results indicate that heat would leave the ocean more readily than it entered, but even if thermal expansion were returned to zero, the geographical pattern of sea level would be altered. Therefore, despite any aggressive CO2 mitigation, regional sea level change is inevitable.


2020 ◽  
Author(s):  
Thomas Kleiner ◽  
Jeremie Schmiedel ◽  
Angelika Humbert

<p>Ice sheets constitute the largest and most uncertain potential source of future sea-level rise. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) brings together a consortium of international ice sheet and climate models to explore the contribution from the Greenland and Antarctic ice sheets to future sea-level rise.</p> <p>We use the Parallel Ice Sheet Model (PISM, pism-docs.org) to carry out spinup and projection simulations for the Antarctic Ice Sheet. Our treatment of the ice-ocean boundary condition previously based on 3D ocean temperatures (initMIP-Antarctica) has been adopted to use the ISMIP6 parameterisation and 3D ocean forcing fields (temperature and salinity) according to the ISMIP6 protocol.</p> <p>In this study, we analyse the impact of the choices made during the model initialisation procedure on the initial state. We present the AWI PISM results of the ISMIP6 projection simulations and investigate the ice sheet response for individual basins. In the analysis, we distinguish between the local and non-local ice shelf basal melt parameterisation.</p>


2020 ◽  
Author(s):  
Jim Jordan ◽  
Hilmar Gudmundsson ◽  
Adrian Jenkins ◽  
Chris Stokes ◽  
Stewart Jamieson ◽  
...  

<p>The East Antarctic Ice Sheet (EAIS) is the single largest potential contributor to future global mean sea level rise, containing a water mass equivalent of 53 m. Recent work has found the overall mass balance of the EAIS to be approximately in equilibrium, albeit with large uncertainties. However, changes in oceanic conditions have the potential to upset this balance. This could happen by both a general warming of the ocean and also by shifts in oceanic conditions allowing warmer water masses to intrude into ice shelf cavities.</p><p>We use the Úa numerical ice-flow model, combined with ocean-melt rates parameterized by the PICO box mode, to predict the future contribution to global-mean sea level of the EAIS. Results are shown for the next 100 years under a range of emission scenarios and oceanic conditions on a region by region basis, as well as for the whole of the EAIS. </p>


2020 ◽  
Author(s):  
Samuel Helsen ◽  
Sam Vanden Broucke ◽  
Alexandra Gossart ◽  
Niels Souverijns ◽  
Nicole van Lipzig

<p>The Totten glacier is a highly dynamic outlet glacier, situated in E-Antarctica, that contains a potential sea level rise of about 3.5 meters. During recent years, this area has been influenced by sub-shelf intrusion of warm ocean currents, contributing to higher basal melt rates. Moreover, most of the ice over this area is grounded below sea level, which makes the ice shelf potentially vulnerable to the marine ice sheet instability mechanism. It is expected that, as a result of climate change, the latter mechanisms may contribute to significant ice losses in this region within the next decades, thereby contributing to future sea level rise. Up to now, most studies have been focusing on sub-shelf melt rates and the influence of the ocean, with much less attention for atmospheric processes (often ignored), which also play a key-role in determining the climatic conditions over this region. For example: surface melt is important because it contributes to hydrofracturing, a process that may lead to ice cliff instabilities. Also precipitation is an important atmospheric process, since it determines the input of mass to the ice sheet and contributes directly to the surface mass balance. In order to perform detailed studies on these processes, we need a well-evaluated climate model that represents all these processes well. Recently, the COSMO-CLM<sup>2</sup> (CCLM<sup>2</sup>) model was adapted to the climatological conditions over Antarctica. The model was evaluated by comparing a 30 year Antarctic-wide hindcast run (1986-2016) at 25 km resolution with meteorological observational products (Souverijns et al., 2019). It was shown that the model performance is comparable to other state-of-the-art regional climate models over the Antarctic region. We now applied the CCLM<sup>2</sup> model in a regional configuration over the Totten glacier area (E-Antarctica) at 5 km resolution and evaluated its performance over this region by comparing it to climatological observations from different stations. We show that the performance for temperature in the high resolution run is comparable to the performance of the Antarctic-wide run. Precipitation is, however, overestimated in the high-resolution run, especially over dome structures (Law-Dome). Therefore, we applied an orographic smoothening, which clearly improves the precipitation pattern with respect to observations. Wind speed is overestimated in some places, which is solved by increasing the surface roughness. This research frames in the context of the PARAMOUR project. Within PARAMOUR, CCLM<sup>2 </sup>is currently being coupled to an ocean model (NEMO) and an ice sheet model (f.ETISh/BISICLES) in order to understand decadal predictability over this region.</p>


2014 ◽  
Vol 7 (5) ◽  
pp. 1933-1943 ◽  
Author(s):  
W. Chang ◽  
P. J. Applegate ◽  
M. Haran ◽  
K. Keller

Abstract. Computer models of ice sheet behavior are important tools for projecting future sea level rise. The simulated modern ice sheets generated by these models differ markedly as input parameters are varied. To ensure accurate ice sheet mass loss projections, these parameters must be constrained using observational data. Which model parameter combinations make sense, given observations? Our method assigns probabilities to parameter combinations based on how well the model reproduces the Greenland Ice Sheet profile. We improve on the previous state of the art by accounting for spatial information and by carefully sampling the full range of realistic parameter combinations, using statistically rigorous methods. Specifically, we estimate the joint posterior probability density function of model parameters using Gaussian process-based emulation and calibration. This method is an important step toward calibrated probabilistic projections of ice sheet contributions to sea level rise, in that it uses data–model fusion to learn about parameter values. This information can, in turn, be used to make projections while taking into account various sources of uncertainty, including parametric uncertainty, data–model discrepancy, and spatial correlation in the error structure. We demonstrate the utility of our method using a perfect model experiment, which shows that many different parameter combinations can generate similar modern ice sheet profiles. This result suggests that the large divergence of projections from different ice sheet models is partly due to parametric uncertainty. Moreover, our method enables insight into ice sheet processes represented by parameter interactions in the model.


2020 ◽  
Author(s):  
Michiel van den Broeke ◽  
Brice Noël ◽  
Leo van Kampenhout ◽  
Willem-Jan van de Berg

<p>The mass balance of the Greenland ice sheet (GrIS, units Gt per year) equals the surface mass balance (SMB) minus solid ice discharge across the grounding line. As the latter is definite positive, an important threshold for irreversible GrIS mass loss occurs when long-term average SMB becomes negative. For this to happen, runoff (mainly meltwater, some rain) must exceed mass accumulation (mainly snowfall minus sublimation). Even for a single year, this threshold has not been passed since at least 1958, the first year with reliable estimates of SMB components, although recent years with warm summers (e.g. 2012 and 2019) came close. Simply extrapolating the recent (1991-present) negative SMB trend into the future suggests that the SMB = 0 threshold could be reached before ~2040, but such predictions are extremely uncertain given the very large interannual SMB variability, the relative brevity of the time series and the uncertainty in future warming. In this study we use a cascade of models, extensively evaluated with in-situ and remotely sensed (GRACE) SMB observations, to better constrain the future regional warming threshold for the 5-year average GrIS SMB to become negative. To this end, a 1950-2100 climate change run with the global model CESM2 (app. 100 km resolution) was dynamically downscaled using the regional climate model RACMO2 (app. 11 km), which in turn was statistically downscaled to 1 km resolution. The result is a threshold regional Greenland warming of close to 4 degrees. We then use a range of CMIP5 and CMIP6 global climate models to translate the regional value into a global warming threshold for various warming scenarios, including its timing this century. We find substantial differences, ranging from stabilization before the threshold is reached in the RCP/SSP2.6 scenarios with a limited but still significant sea-level rise contribution (< 5 cm by 2100) to an imminent crossing of the warming threshold for the RCP/SSP8.5 scenarios with substantial and ever-growing contributions to sea level rise (> 10 cm by 2100). These results stress the need for strong mitigation to avoid irreversible GrIS mass loss. We finish by discussing the caveats and uncertainties of our approach.</p>


2020 ◽  
Author(s):  
Emily Hill ◽  
Sebastian Rosier ◽  
Hilmar Gudmundsson ◽  
Matthew Collins

<p>Mass loss from the Antarctic Ice Sheet is the main source of uncertainty in projections of future sea-level rise, with important implications for coastal regions worldwide. Enhanced melt beneath ice shelves could destabilise large parts of the ice sheet, and further increase ice loss. Despite advances in our understanding of feedbacks in the ice sheet-ice shelf-ocean system, future projections of ice loss remain poorly constrained in many parts of Antarctica. In particular, there is ongoing debate surrounding the future of the Filchner-Ronne Ice Shelf (FRIS) region. The FRIS has remained relatively unchanged in recent decades, but an increase in air and ocean temperatures in the neighbouring Weddell Sea, could force rapid retreat in the near future. Indeed, previous modelling work has suggested the potential for widespread infiltration of warm water beneath the ice shelf in the second half of the twenty-first century, leading to a drastic increase in basal melting.</p><p>Here, we use the ice flow model Úa alongside the ocean box model PICO (Potsdam Ice-shelf Cavity mOdel) to understand the key physical processes and model variability in future projections of sea level rise from the FRIS region. We investigate uncertain model parameters associated with ice dynamics, surface melting and precipitation, ocean temperature forcing, and parameters relating to the strength of basal melt generated by PICO. We optimise the prior distributions of parameters in PICO using observations and a Bayesian approach, leading to improved posterior distributions for use in the following stages of uncertainty quantification. We then run our forward model through the 21<sup>st</sup> century for various RCP scenarios and extensive random sampling of uncertain parameters to train an emulator. From this, we present probabilistic projections of potential sea level rise from the FRIS region for different future climate change scenarios, together with a sensitivity analysis to identify the most important parameters that contribute to uncertainty in these projections.</p>


2020 ◽  
Author(s):  
Stefan Hofer ◽  
Charlotte Lang ◽  
Charles Amory ◽  
Christoph Kittel ◽  
Alison Delhasse ◽  
...  

<p>Future climate projections show a marked increase in Greenland Ice Sheet (GrIS) runoff<br>during the 21st century, a direct consequence of the Polar Amplification signal. Regional<br>climate models (RCMs) are a widely used tool to downscale ensembles of projections from<br>global climate models (GCMs) to assess the impact of global warming on GrIS melt and<br>sea level rise contribution. Initial results of the CMIP6 GCM model intercomparison<br>project have revealed a greater 21st century temperature rise than in CMIP5 models.<br>However, so far very little is known about the subsequent impacts on the future GrIS<br>surface melt and therefore sea level rise contribution. Here, we show that the total GrIS<br>melt during the 21st century almost doubles when using CMIP6 forcing compared to the<br>previous CMIP5 model ensemble, despite an equal global radiative forcing of +8.5 W/m2<br>in 2100 in both RCP8.5 and SSP58.5 scenarios. The total GrIS sea level rise contribution<br>from surface melt in our high-resolution (15 km) projections is 17.8 cm in SSP58.5, 7.9 cm<br>more than in our RCP8.5 simulations, despite the same radiative forcing. We identify a<br>+1.7°C greater Arctic amplification in the CMIP6 ensemble as the main driver behind the<br>presented doubling of future GrIS sea level rise contribution</p>


2012 ◽  
Vol 4 (3) ◽  
pp. 212-229 ◽  
Author(s):  
K. Kvale ◽  
K. Zickfeld ◽  
T. Bruckner ◽  
K. J. Meissner ◽  
K. Tanaka ◽  
...  

Abstract Anthropogenic emissions of greenhouse gases could lead to undesirable effects on oceans in coming centuries. Drawing on recommendations published by the German Advisory Council on Global Change, levels of unacceptable global marine change (so-called guardrails) are defined in terms of global mean temperature, sea level rise, and ocean acidification. A global-mean climate model [the Aggregated Carbon Cycle, Atmospheric Chemistry and Climate Model (ACC2)] is coupled with an economic module [taken from the Dynamic Integrated Climate–Economy Model (DICE)] to conduct a cost-effectiveness analysis to derive CO2 emission pathways that both minimize abatement costs and are compatible with these guardrails. Additionally, the “tolerable windows approach” is used to calculate a range of CO2 emissions paths that obey the guardrails as well as a restriction on mitigation rate. Prospects of meeting the global mean temperature change guardrail (2° and 0.2°C decade−1 relative to preindustrial) depend strongly on assumed values for climate sensitivity: at climate sensitivities >3°C the guardrail cannot be attained under any CO2 emissions reduction strategy without mitigation of non-CO2 greenhouse gases. The ocean acidification guardrail (0.2 unit pH decline relative to preindustrial) is less restrictive than the absolute temperature guardrail at climate sensitivities >2.5°C but becomes more constraining at lower climate sensitivities. The sea level rise and rate of rise guardrails (1 m and 5 cm decade−1) are substantially less stringent for ice sheet sensitivities derived in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report, but they may already be committed to violation if ice sheet sensitivities consistent with semiempirical sea level rise projections are assumed.


Ocean Science ◽  
2021 ◽  
Vol 17 (1) ◽  
pp. 181-186
Author(s):  
Aslak Grinsted ◽  
Jens Hesselbjerg Christensen

Abstract. Recent assessments from the Intergovernmental Panel on Climate Change (IPCC) imply that global mean sea level is unlikely to rise more than about 1.1 m within this century but will increase further beyond 2100. Even within the most intensive future anthropogenic greenhouse gas emission scenarios, higher levels are assessed to be unlikely. However, some studies conclude that considerably greater sea level rise could be realized, and a number of experts assign a substantially higher likelihood of such a future. To understand this discrepancy, it would be useful to have scenario-independent metrics that can be compared between different approaches. The concept of a transient climate sensitivity has proven to be useful to compare the global mean temperature response of climate models to specific radiative forcing scenarios. Here, we introduce a similar metric for sea level response. By analyzing the mean rate of change in sea level (not sea level itself), we identify a nearly linear relationship with global mean surface temperature (and therefore accumulated carbon dioxide emissions) both in model projections and in observations on a century scale. This motivates us to define the “transient sea level sensitivity” as the increase in the sea level rate associated with a given warming in units of meters per century per kelvin. We find that future projections estimated on climate model responses fall below extrapolation based on recent observational records. This comparison suggests that the likely upper level of sea level projections in recent IPCC reports would be too low.


Sign in / Sign up

Export Citation Format

Share Document