skip to main content


Title: Experimental protocol for sea level projections from ISMIP6 stand-alone ice sheet models
Abstract. Projection of the contribution of ice sheets to sea level change as part ofthe Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the formof simulations from coupled ice sheet–climate models and stand-alone icesheet models, overseen by the Ice Sheet Model Intercomparison Project forCMIP6 (ISMIP6). This paper describes the experimental setup forprocess-based sea level change projections to be performed with stand-aloneGreenland and Antarctic ice sheet models in the context of ISMIP6. TheISMIP6 protocol relies on a suite of polar atmospheric and oceanicCMIP-based forcing for ice sheet models, in order to explore the uncertaintyin projected sea level change due to future emissions scenarios, CMIPmodels, ice sheet models, and parameterizations for ice–ocean interactions.We describe here the approach taken for defining the suite of ISMIP6stand-alone ice sheet simulations, document the experimental framework andimplementation, and present an overview of the ISMIP6 forcing to beused by participating ice sheet modeling groups.  more » « less
Award ID(s):
2021699 1643733 1744792
NSF-PAR ID:
10191305
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
The Cryosphere
Volume:
14
Issue:
7
ISSN:
1994-0424
Page Range / eLocation ID:
2331 to 2368
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. The ice sheet model intercomparison project for CMIP6 (ISMIP6) effort brings together the ice sheet and climate modeling communities to gain understanding of the ice sheet contribution to sea level rise. ISMIP6 conducts stand-alone ice sheet experiments that use space- and time-varying forcing derived from atmosphere–ocean coupled global climate models (AOGCMs) to reflect plausible trajectories for climate projections. The goal of this study is to recommend a subset of CMIP5 AOGCMs (three core and three targeted) to produce forcing for ISMIP6 stand-alone ice sheet simulations, based on (i) their representation of current climate near Antarctica and Greenland relative to observations and (ii) their ability to sample a diversity of projected atmosphere and ocean changes over the 21st century. The selection is performed separately for Greenland and Antarctica. Model evaluation over the historical period focuses on variables used to generate ice sheet forcing. For stage (i), we combine metrics of atmosphere and surface ocean state (annual- and seasonal-mean variables over large spatial domains) with metrics of time-mean subsurface ocean temperature biases averaged over sectors of the continental shelf. For stage (ii), we maximize the diversity of climate projections among the best-performing models. Model selection is also constrained by technical limitations, such as availability of required data from RCP2.6 and RCP8.5 projections. The selected top three CMIP5 climate models are CCSM4, MIROC-ESM-CHEM, and NorESM1-M for Antarctica and HadGEM2-ES, MIROC5, and NorESM1-M for Greenland. This model selection was designed specifically for ISMIP6 but can be adapted for other applications. 
    more » « less
  2. null (Ed.)
    Abstract. Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution inresponse to different climate scenarios and assess the mass loss that would contribute tofuture sea level rise. However, there is currently no consensus on estimates of the future massbalance of the ice sheet, primarily because of differences in the representation of physicalprocesses, forcings employed and initial states of ice sheet models. This study presentsresults from ice flow model simulations from 13 international groups focusing on the evolutionof the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet ModelIntercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from theCoupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climatemodel results. Simulations of the Antarctic ice sheet contribution to sea level rise in responseto increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent(SLE) under Representative ConcentrationPathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment withconstant climate conditions and should therefore be added to the mass loss contribution underclimate conditions similar to present-day conditions over the same period. The simulated evolution of theWest Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighingthe increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelfcollapse, here assumed to be caused by large amounts of liquid water ponding at the surface ofice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without iceshelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, thecalibration of these melt rates based on oceanic conditions taken outside of ice shelf cavitiesand the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario basedon two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared tosimulations done under present-day conditions for the two CMIP5 forcings used and displaylimited mass gain in East Antarctica. 
    more » « less
  3. Abstract. Changes in ocean temperature and salinity are expected to be an important determinant of the Greenland ice sheet's future sea level contribution. Yet, simulating the impact of these changes in continental-scale ice sheet models remains challenging due to the small scale of key physics, such as fjord circulation and plume dynamics, and poor understanding of critical processes, such as calving and submarine melting. Here we present the ocean forcing strategy for Greenland ice sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change Sixth Assessment Report. Beginning from global atmosphere–ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland ice sheet models, termed the “retreat” and “submarine melt” implementations. The retreat implementation parameterises glacier retreat as a function of projected subglacial discharge and ocean thermal forcing, is designed to be implementable by all ice sheet models and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios, respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly constrained parameterisations and, as such, further research on submarine melting, calving and fjord–shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland ice sheet models to be systematically and consistently forced by the ocean for the first time and should result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change. 
    more » « less
  4. null (Ed.)
    Abstract. Climate model projections have previously been used to compute ice shelf basal melt rates in ice sheet models, but the strategies employed – e.g., ocean input, parameterization, calibration technique, and corrections – have varied widely and are often ad hoc. Here, a methodology is proposed for the calculation of circum-Antarctic basal melt rates for floating ice, based on climate models, that is suitable for ISMIP6, the Ice Sheet Model Intercomparison Project for CMIP6 (6th Coupled Model Intercomparison Project). The past and future evolution of ocean temperature and salinity is derived from a climate model by estimating anomalies with respect to the modern day, which are added to a present-day climatology constructed from existing observational datasets. Temperature and salinity are extrapolated to any position potentially occupied by a simulated ice shelf. A simple formulation is proposed for a basal melt parameterization in ISMIP6, constrained by the observed temperature climatology, with a quadratic dependency on either the nonlocal or local thermal forcing. Two calibration methods are proposed: (1) based on the mean Antarctic melt rate (MeanAnt) and (2) based on melt rates near Pine Island's deep grounding line (PIGL). Future Antarctic mean melt rates are an order of magnitude greater in PIGL than in MeanAnt. The PIGL calibration and the local parameterization result in more realistic melt rates near grounding lines. PIGL is also more consistent with observations of interannual melt rate variability underneath Pine Island and Dotson ice shelves. This work stresses the need for more physics and less calibration in the parameterizations and for more observations of hydrographic properties and melt rates at interannual and decadal timescales. 
    more » « less
  5. Abstract

    The Antarctic ice sheet (AIS) will be a dominant contributor to global mean sea level rise in the 21st century but remains a major source of uncertainty. The Ice Sheet Model Intercomparison for CMIP6 (ISMIP6) is an ensemble of continental‐scale models for studying the evolution of the AIS and projecting its future contribution to sea level. Due to their complexity and computational cost, ISMIP6 simulations are sparse and generated infrequently. Emulators are smaller‐scale models that approximate ISMs and enable experimentation and exploration into the drivers of sea level change. We introduce a neural network (NN) emulator to approximate the ISMIP6 ensemble, using a variational Long Short‐Term Memory (LSTM) with Monte Carlo dropout to quantify single‐projection uncertainty. The proposed NN emulator is compared to a Gaussian Process (GP) emulator on four criteria: accuracy of point estimates and predictive distributions of individual model projections, approximation of the ensemble projections, and model training time. The NN predicts more accurately on single projections, with a mean absolute error of 0.46 mm Sea Level Equivalent (SLE) versus 0.73 mm SLE for the GP, and has more accurate uncertainty estimates. The NN emulator also better approximates the ensemble distribution of ISMIP6 model projections, with a Kullback‐Leibler divergence of 18.26 versus 199.14 for GP at the projection year 2100. The NN enables more accurate experimentation with a reduced runtime, offering a new tool for understanding the important role of regional precipitation, ice sheet drainage systems, and interannual and longer timescale dynamics.

     
    more » « less