Intercomparison Project For CMIP6 Research Articles

Local forcing mechanisms challenge parameterizations of ocean thermal forcing for Greenland tidewater glaciers

Abstract. Frontal ablation has caused 32 %–66 % of Greenland Ice Sheet mass loss since 1972, and despite its importance in driving terminus change, ocean thermal forcing remains crudely incorporated into large-scale ice sheet models. In Greenland, local fjord-scale processes modify the magnitude of thermal forcing at the ice–ocean boundary but are too small scale to be resolved in current global climate models. For example, simulations used in the Ice Sheet Intercomparison Project for CMIP6 (ISMIP6) to predict future ice sheet change rely on the extrapolation of regional ocean water properties into fjords to drive terminus ablation. However, the accuracy of this approach has not previously been tested due to the scarcity of observations in Greenland fjords, as well as the inability of fjord-scale models to realistically incorporate icebergs. By employing the recently developed IceBerg package within the Massachusetts Institute of Technology general circulation model (MITgcm), we here evaluate the ability of ocean thermal forcing parameterizations to predict thermal forcing at tidewater glacier termini. This is accomplished through sensitivity experiments using a set of idealized Greenland fjords, each forced with equivalent ocean boundary conditions but with varying tidal amplitudes, subglacial discharge, iceberg coverage, and bathymetry. Our results indicate that the bathymetric obstruction of external water is the primary control on near-glacier thermal forcing, followed by iceberg submarine melting. Despite identical ocean boundary conditions, we find that the simulated fjord processes can modify grounding line thermal forcing by as much as 3 °C, the magnitude of which is largely controlled by the relative depth of bathymetric sills to the Polar Water–Atlantic Water thermocline. However, using a common adjustment for fjord bathymetry we can still predict grounding line thermal forcing within 0.2 °C in our simulations. Finally, we introduce new parameterizations that additionally account for iceberg-driven cooling that can accurately predict interior fjord thermal forcing profiles both in iceberg-laden simulations and in observations from Kangiata Sullua (Ilulissat Icefjord).

Mass loss of the Antarctic ice sheet until the year 3000 under a sustained late-21st-century climate

AbstractIce-sheet simulations of Antarctica extending to the year 3000 are analysed to investigate the long-term impacts of 21st-century warming. Climate projections are used as forcing until 2100 and afterwards no climate trend is applied. Fourteen experiments are for the ‘unabated warming’ pathway, and three are for the ‘reduced emissions’ pathway. For the unabated warming path simulations, West Antarctica suffers a much more severe ice loss than East Antarctica. In these cases, the mass loss amounts to an ensemble average of ~3.5 m sea-level equivalent (SLE) by the year 3000 and ~5.3 m for the most sensitive experiment. Four phases of mass loss occur during the collapse of the West Antarctic ice sheet. For the reduced emissions pathway, the mean mass loss is ~0.24 m SLE. By demonstrating that the consequences of the 21st century unabated warming path forcing are large and long term, the results present a different perspective to ISMIP6 (Ice Sheet Model Intercomparison Project for CMIP6). Extended ABUMIP (Antarctic BUttressing Model Intercomparison Project) simulations, assuming sudden and sustained ice-shelf collapse, with and without bedrock rebound, corroborate a negative feedback for ice loss found in previous studies, where bedrock rebound acts to slow the rate of ice loss.

ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model

Abstract. The future retreat rate for marine-based regions of the Antarctic Ice Sheet is one of the largest uncertainties in sea-level projections. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) aims to improve projections and quantify uncertainties by running an ensemble of ice sheet models with atmosphere and ocean forcing derived from global climate models. Here, the Community Ice Sheet Model (CISM) is used to run ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution. Using multiple combinations of sub-ice-shelf melt parameterizations and calibrations, CISM is spun up to steady state over many millennia. During the spin-up, basal friction parameters and basin-scale thermal forcing corrections are adjusted to optimize agreement with the observed ice thickness. The model is then run forward for 550 years, from 1950–2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, the ocean forcing triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner–Ronne and Ross sectors. Mass loss accelerates late in the 21st century and then rises steadily for several centuries without leveling off. The resulting ocean-forced sea-level rise at year 2500 varies from about 150 to 1300 mm, depending on the melt scheme and ocean forcing. Further experiments show relatively high sensitivity to the basal friction law, moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves, and low sensitivity to the stress-balance approximation. The Amundsen sector exhibits threshold behavior, with modest retreat under many parameter settings but complete collapse under some combinations of low basal friction and high thermal forcing anomalies. Large uncertainties remain, as a result of parameterized sub-shelf melt rates, simplified treatments of calving and basal friction, and the lack of ice–ocean coupling.

Sensitivity of Greenland ice sheet projections to spatial resolution in higher-order simulations: the Alfred Wegener Institute (AWI) contribution to ISMIP6 Greenland using the Ice-sheet and Sea-level System Model (ISSM)

Abstract. Projections of the contribution of the Greenland ice sheet to future sea-level rise include uncertainties primarily due to the imposed climate forcing and the initial state of the ice sheet model. Several state-of-the-art ice flow models are currently being employed on various grid resolutions to estimate future mass changes in the framework of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). Here we investigate the sensitivity to grid resolution of centennial sea-level contributions from the Greenland ice sheet and study the mechanism at play. We employ the finite-element higher-order Ice-sheet and Sea-level System Model (ISSM) and conduct experiments with four different horizontal resolutions, namely 4, 2, 1 and 0.75 km. We run the simulation based on the ISMIP6 core climate forcing from the MIROC5 global circulation model (GCM) under the high-emission Representative Concentration Pathway (RCP) 8.5 scenario and consider both atmospheric and oceanic forcing in full and separate scenarios. Under the full scenarios, finer simulations unveil up to approximately 5 % more sea-level rise compared to the coarser resolution. The sensitivity depends on the magnitude of outlet glacier retreat, which is implemented as a series of retreat masks following the ISMIP6 protocol. Without imposed retreat under atmosphere-only forcing, the resolution dependency exhibits an opposite behaviour with approximately 5 % more sea-level contribution in the coarser resolution. The sea-level contribution indicates a converging behaviour below a 1 km horizontal resolution. A driving mechanism for differences is the ability to resolve the bedrock topography, which highly controls ice discharge to the ocean. Additionally, thinning and acceleration emerge to propagate further inland in high resolution for many glaciers. A major response mechanism is sliding, with an enhanced feedback on the effective normal pressure at higher resolution leading to a larger increase in sliding speeds under scenarios with outlet glacier retreat.

The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6

Abstract. The Greenland ice sheet is one of the largest contributors to global mean sea-level rise today and is expected to continue to lose mass as the Arctic continues to warm. The two predominant mass loss mechanisms are increased surface meltwater run-off and mass loss associated with the retreat of marine-terminating outlet glaciers. In this paper we use a large ensemble of Greenland ice sheet models forced by output from a representative subset of the Coupled Model Intercomparison Project (CMIP5) global climate models to project ice sheet changes and sea-level rise contributions over the 21st century. The simulations are part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We estimate the sea-level contribution together with uncertainties due to future climate forcing, ice sheet model formulations and ocean forcing for the two greenhouse gas concentration scenarios RCP8.5 and RCP2.6. The results indicate that the Greenland ice sheet will continue to lose mass in both scenarios until 2100, with contributions of 90±50 and 32±17 mm to sea-level rise for RCP8.5 and RCP2.6, respectively. The largest mass loss is expected from the south-west of Greenland, which is governed by surface mass balance changes, continuing what is already observed today. Because the contributions are calculated against an unforced control experiment, these numbers do not include any committed mass loss, i.e. mass loss that would occur over the coming century if the climate forcing remained constant. Under RCP8.5 forcing, ice sheet model uncertainty explains an ensemble spread of 40 mm, while climate model uncertainty and ocean forcing uncertainty account for a spread of 36 and 19 mm, respectively. Apart from those formally derived uncertainty ranges, the largest gap in our knowledge is about the physical understanding and implementation of the calving process, i.e. the interaction of the ice sheet with the ocean.

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 of the Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the form of simulations from coupled ice sheet–climate models and stand-alone ice sheet models, overseen by the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). This paper describes the experimental setup for process-based sea level change projections to be performed with stand-alone Greenland and Antarctic ice sheet models in the context of ISMIP6. The ISMIP6 protocol relies on a suite of polar atmospheric and oceanic CMIP-based forcing for ice sheet models, in order to explore the uncertainty in projected sea level change due to future emissions scenarios, CMIP models, ice sheet models, and parameterizations for ice–ocean interactions. We describe here the approach taken for defining the suite of ISMIP6 stand-alone ice sheet simulations, document the experimental framework and implementation, and present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups.

Twenty-first century ocean forcing of the Greenland ice sheet for modelling of sea level contribution

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.

CMIP5 model selection for ISMIP6 ice sheet model forcing: Greenland and Antarctica

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.