skip to main content


Title: A protocol for calculating basal melt rates in the ISMIP6 Antarctic ice sheet projections
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
Award ID(s):
1744792
NSF-PAR ID:
10206399
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
14
Issue:
9
ISSN:
1994-0424
Page Range / eLocation ID:
3111 to 3134
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Mass loss from the Amundsen Sea sector of the West Antarctic Ice Sheet has increased in recent decades, suggestive of sus- tained ocean forcing or an ongoing, possibly unstable, response to a past climate anomaly. Lengthening satellite records appear to be incompatible with either process, however, revealing both periodic hiatuses in acceleration and intermittent episodes of thinning. Here we use ocean temperature, salinity, dissolved-oxygen and current measurements taken from 2000 to 2016 near the Dotson Ice Shelf to determine temporal changes in net basal melting. A decadal cycle dominates the ocean record, with melt changing by a factor of about four between cool and warm extremes via a nonlinear relationship with ocean temperature. A warm phase that peaked around 2009 coincided with ice-shelf thinning and retreat of the grounding line, which re-advanced during a post-2011 cool phase. These observations demonstrate how discontinuous ice retreat is linked with ocean variability, and that the strength and timing of decadal extremes is more influential than changes in the longer-term mean state. The non- linear response of melting to temperature change heightens the sensitivity of Amundsen Sea ice shelves to such variability, possibly explaining the vulnerability of the ice sheet in that sector, where subsurface ocean temperatures are relatively high. 
    more » « less
  4. The loss or thinning of buttressing ice shelves and accompanying changes in grounding-zone stress balance are commonly implicated as the primary trigger for grounding-line retreat, such as that observed in Amundsen Sea outlet glaciers today. Ice-shelf thinning is mostly attributed to the presence of warm ocean waters beneath the shelves. However, climate model projections show that summer air temperatures could soon exceed the threshold for widespread meltwater production on ice-shelf surfaces. This has serious implications for their future stability, because they are vulnerable to water-induced flexural stresses and water-aided crevasse penetration, termed ‘hydrofracturing’. Once initiated, the rate of shelf loss through hydrofracturing can far exceed that caused by sub-surface melting, and could result in the complete loss of some buttressing ice shelves, with marine grounding lines suddenly becoming calving ice fronts. In places where those exposed ice fronts are thick (>900m) and crevassed, deviatoric stresses can exceed the strength of the ice and the cliff face will fail mechanically, leading to rapid calving like that seen in analogous settings such as Jakobshavn on Greenland. Here we explore the implications of hydrofacturing and subsequent ice-cliff collapse in a warming climate, by parameterizing these processes in a hybrid ice sheet-shelf model. Model sensitivities to meltwater production and to ice-cliff calving rate (a function of cliff height above the stress balance threshold triggering brittle failure) are calibrated to match modern observations of calving and thinning. We find the potential for major ice-sheet retreat if global mean temperature rises more than ~2ºC above preindustrial. In the model, Antarctic calving rates at thick ice fronts are not allowed to exceed those observed in Greenland today. This may be a conservative assumption, considering the very different spatial scales of Antarctic outlets, such as Thwaites. Nonetheless, simulations following a ‘worst case’ RCP8.5 scenario produce rates of sea-level rise measured in cm per year by the end of this century. Clearly, the potential for brittle processes to deliver ice to the ocean, in addition to viscous and basal processes, needs to be better constrained through more complete, physically based representations of calving. 
    more » « less
  5. null (Ed.)
    This study presents near future (2020–2044) temperature and precipitation changes over the Antarctic Peninsula under the high-emission scenario (RCP8.5). We make use of historical and projected simulations from 19 global climate models (GCMs) participating in Coupled Model Intercomparison Project phase 5 (CMIP5). We compare and contrast GCMs projections with two groups of regional climate model simulations (RCMs): (1) high resolution (15-km) simulations performed with Polar-WRF model forced with bias-corrected NCAR-CESM1 (NC-CORR) over the Antarctic Peninsula, (2) medium resolution (50-km) simulations of KNMI-RACMO21P forced with EC-EARTH (EC) obtained from the CORDEX-Antarctica. A further comparison of historical simulations (1981–2005) with respect to ERA5 reanalysis is also included for circulation patterns and near-surface temperature climatology. In general, both RCM boundary conditions represent well the main circulation patterns of the historical period. Nonetheless, there are important differences in projections such as a notable deepening and weakening of the Amundsen Sea Low in EC and NC-CORR, respectively. Mean annual near-surface temperatures are projected to increase by about 0.5–1.5 ◦ C across the entire peninsula. Temperature increase is more substantial in autumn and winter ( ∼ 2 ◦C). Following opposite circulation pattern changes, both EC and NC-CORR exhibit different warming rates, indicating a possible continuation of natural decadal variability. Although generally showing similar temperature changes, RCM projections show less warming and a smaller increase in melt days in the Larsen Ice Shelf compared to their respective driving fields. Regarding precipitation, there is a broad agreement among the simulations, indicating an increase in mean annual precipitation ( ∼ 5 to 10%). However, RCMs show some notable differences over the Larsen Ice Shelf where total precipitation decreases (for RACMO) and shows a small increase in rain frequency. We conclude that it seems still difficult to get consistent projections from GCMs for the Antarctic Peninsula as depicted in both RCM boundary conditions. In addition, dominant and common changes from the boundary conditions are largely evident in the RCM simulations. We argue that added value of RCM projections is driven by processes shaped by finer local details and different physics schemes that are introduced by RCMs, particularly over the Larsen Ice Shelf. 
    more » « less