Search

In-situ observations from the Filchner-Ronne Ice Shelf (FRIS) have uncovered dominant time scales of variability in basal melting and circulation beneath this extensive ice shelf. In particular, the data characterize mechanisms of seasonal and inter-annual variability in sub-ice shelf properties, and show that the amplitude of the variability over the past thirty years is very modest.Because accurate representation of variability under present-day climate is an obvious prerequisite for earth system models that aim to project climate under a future change, this new observational understanding presents an opportunity for critical evaluation and improvement of existing models. We focus on the Energy Exascale Earth System Model (E3SM) and through a series of simulations we investigate the impact of ocean mixing parameterizations on the variability in the FRIS cavity.

Evolution of ice sheets contribute to sea-level change globally by exchanging mass with the ocean, and regionally by causing the solid Earth deformation and perturbation of the Earth’s rotation and gravitational field, so-called “gravitational, rotational and deformational (GRD) effects”. In the last decade, much work has been done to establish the importance of coupling GRD effects particularly in modeling of marine-based ice sheets (e.g., West Antarctic Ice Sheet; WAIS) to capture the interactions between ice sheets, sea level and the solid Earth at the grounding lines. However, coupling of GRD effects has not yet been done widely within the ice-sheet modeling community; for example, GRD effects were not included in any of the ice sheet models that contributed to the most recent recent ice-sheet model intercomparison through 2100 (Ice Sheet Model Intercomparison Project for CMIP6: ISMIP6-2100; Serrousi et al., 2020) cited by the latest IPCC AR6 report.In this work, we couple the US Department of Energy’s MPAS-Albany Land Ice model (which was one of the models that participated in the ISMIP6-2100 project) to a 1D sea-level model and perform coupled simulations of Antarctica under the new ISMIP6-2300 protocol in which climate forcing is extended beyond 2100 to 2300. Comparing to the standalone ice-sheet simulations with fix bed topography without GRD effects, the results from our coupled simulations show multi-decadal to centennial-scale delays in the retreat of the Thwaites glacier in the West Antarctica. Our results further suggest that the strength of the negative feedback of sea-level changes on the WAIS retreat becomes weaker as the strength of the applied forcing increases, implying the pertinence of our commitment to limiting greenhouse gas emissions. In addition, within our coupled ice sheet-sea level modeling frame, we introduce a new workflow work in which the ISMIP6 protocol-provided ocean thermal forcing is re-extrapolated based on the updated ocean bathymetry. Our preliminary results indicate that bedrock uplift due to ice mass loss can block the bottom warm ocean, providing additional negative feedback, but also can block cold water when/if the vertical ocean temperature profile gets inverted due to climate change (e.g., as represented in the UKESM model - SSP585 scenario results).

We present our approach to coupling an ocean component (MPAS-Ocean, Model for Prediction Across Scales-Ocean) to an ice sheet component (MALI, MPAS-Albany Land Ice) within an Earth System Model (E3SM, the Energy Exascale Earth System Model) developed by the US Department of Energy. First, we present an extrapolation technique, similar to the ISMIP6 (Ice Sheet Modeling for CMIP6) protocol, that can be used in the absence of evolving grounding lines in the ocean component. This technique, while crude, can be used in both Greenland fjords and ice-shelf cavities as a stop-gap in situations where the ocean component cannot capture the topographic evolution (e.g. because the ocean grid is too coarse or full coupling has not yet been completed). Second, we demonstrate progress on a fully conservative wetting-and-drying technique using the idealized MISOMIP1 (Marine Ice Sheet-Ocean Intercomparison Project, phase 1) experiments within E3SM.

High-resolution simulation of global climate physics enables us to model how the climate may change under a variety of future scenarios. Such simulations produce vast amounts of information and dense datasets. If interrogated in tandem, these datasets can provide holistic, vital information on Earth's many integrated systems by revealing the manifold interrelated properties of the atmosphere, ocean, and polar ice, framed by real-world terrain in three-dimensional space as they vary over time. To accomplish this, climate scientists have joined with computer scientists and an artist to develop techniques enabling scientists to see these relationships. The impact of ocean water properties on Antarctic ice shelves illustrates the benefit of this analysis in understanding land ice melt rates and thus sea-level rise.

The processes that govern freshwater flux from the Antarctic Ice Sheet (AIS)—ice-shelf basal melting and iceberg calving—are generally poorly represented in current Earth System Models (ESMs). The processes governing ocean flows onto the Antarctic continental and into ice-shelf cavities can only be captured accurately at resolutions significantly higher than those in typical CMIP-class ESMs. The Energy Exascale Earth System Model (E3SM) from the US Department of Energy supports regional refinement in all components, allowing global modeling with high resolution in regions of interest. Here, we present fully coupled results from an ocean/sea-ice mesh that has high resolution (12 km) on the Antarctic continental shelf and much of the Southern Ocean and low resolution (~30 to 60 km) over the rest of the globe. E3SM includes Antarctic ice-shelf cavities with fixed geometry and calculates ice-shelf basal melt rates from the heat and freshwater fluxes computed by the ocean component. In addition, E3SM permits prescribed forcing from a climatology of iceberg melt, providing a more realistic representation of these freshwater fluxes than found in many ESMs. With these new capabilities, E3SM version 2 produces realistic and stable ice-shelf basal melt rates across the continent. We show preliminary results of modeled ice-shelf basal melt rates across a range of Antarctic ice-shelves under pre-industrial and historical climate forcing, as well as the impacts of these added capabilities on the region’s climate. We show that the use of a mesoscale eddy parameterization, tapered with the mesh resolution, reduces biases even in the 12-km region where some eddies are resolved. The accurate representation of these processes within a coupled ESM is an important step towards reducing uncertainties in projections of the Antarctic response to climate change and Antarctica's contribution to global sea-level rise.

Using a numerical ice flow model, we study changes in ice shelf buttressing and grounding-line flux due to localized ice thickness perturbations, a proxy for localized changes in sub-ice-shelf melting. From our experiments, applied to idealized (MISMIP+) and realistic (Larsen C) ice shelf domains, we identify a correlation between a locally derived buttressing number on the ice shelf, based on the first principal stress, and changes in the integrated grounding-line flux. The origin of this correlation, however, remains elusive from the perspective of a theoretical or physically based understanding. This and the fact that the correlation is generally much poorer when applied to realistic ice shelf domains motivate us to seek an alternative approach for predicting changes in grounding-line flux. We therefore propose an adjoint-based method for calculating the sensitivity of the integrated grounding-line flux to local changes in ice shelf geometry. We show that the adjoint-based sensitivity is identical to that deduced from pointwise, diagnostic model perturbation experiments. Based on its much wider applicability and the significant computational savings, we propose that the adjoint-based method is ideally suited for assessing grounding-line flux sensitivity to changes in sub-ice-shelf melting.

The Southern Ocean overturning circulation is driven by winds, heat fluxes, and freshwater sources. Among these sources of freshwater, Antarctic sea ice formation and melting play the dominant role. Even though ice-shelf melt is relatively small in magnitude, it is located close to regions of convection, where it may influence dense water formation. Here, we explore the impacts of ice-shelf melting on Southern Ocean water-mass transformation (WMT) using simulations from the Energy Exascale Earth System Model (E3SM) both with and without the explicit representation of melt fluxes from beneath Antarctic ice shelves. We find that ice-shelf melting enhances transformation of Upper Circumpolar Deep Water, converting it to lower density values. While the overall differences in Southern Ocean WMT between the two simulations are moderate, freshwater fluxes produced by ice-shelf melting have a further, indirect impact on the Southern Ocean overturning circulation through their interaction with sea ice formation and melting, which also cause considerable upwelling. We further find that surface freshening and cooling by ice-shelf melting cause increased Antarctic sea ice production and stronger density stratification near the Antarctic coast. In addition, ice-shelf melting causes decreasing air temperature, which may be directly related to sea ice expansion. The increased stratification reduces vertical heat transport from the deeper ocean. Although the addition of ice-shelf melting processes leads to no significant changes in Southern Ocean WMT, the simulations and analysis conducted here point to a relationship between increased Antarctic ice-shelf melting and the increased role of sea ice in Southern Ocean overturning.

Small-scale turbulent flow below ice shelves is regionally isolated and difficult to measure and simulate. Yet these small-scale processes, which regulate heat and salt transfer between the ocean and ice shelves, can affect sea-level rise by altering the ability of Antarctic ice shelves to “buttress” ice flux to the ocean. In this study, we improve our understanding of turbulence below ice shelves by means of large-eddy simulations at sub-meter resolution, capturing boundary layer mixing at scales intermediate between laboratory experiments or direct numerical simulations and regional or global ocean circulation models. Our simulations feature the development of an ice-shelf ocean boundary layer through dynamic ice melting in a regime with low thermal driving, low ice-shelf basal slope, and strong shear driven by the geostrophic flow. We present a preliminary assessment of existing ice-shelf basal melt parameterizations adopted in single component or coupled ice-sheet and ocean models on the basis of a small parameter study. While the parameterized linear relationship between ice-shelf melt rate and far-field ocean temperature appears to be robust, we point out a little-considered relationship between ice-shelf basal slope and melting worthy of further study.

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60km in the midlatitudes and 30km at the equator and poles), and river transport (55km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model's strong aerosol-related effective radiative forcing (ERFari+aci=−1.65W/m2) and high equilibrium climate sensitivity (ECS=5.3K).

The Energy Exascale Earth System Model (E3SM) is a new coupled Earth system model sponsored by the U.S Department of Energy. Here we present E3SM global simulations using active ocean and sea ice that are driven by the Coordinated Ocean‐ice Reference Experiments II (CORE‐II) interannual atmospheric forcing data set. The E3SM ocean and sea ice components are MPAS‐Ocean and MPAS‐Seaice, which use the Model for Prediction Across Scales (MPAS) framework and run on unstructured horizontal meshes. For this study, grid cells vary from 30 to 60 km for the low‐resolution mesh and 6 to 18 km at high resolution. The vertical grid is a structured z‐star coordinate and uses 60 and 80 layers for low and high resolution, respectively. The lower‐resolution simulation was run for five CORE cycles (310 years) with little drift in sea surface temperature (SST) or heat content. The meridional heat transport (MHT) is within observational range, while the meridional overturning circulation at 26.5°N is low compared to observations. The largest temperature biases occur in the Labrador Sea and western boundary currents (WBCs), and the mixed layer is deeper than observations at northern high latitudes in the winter months. In the Antarctic, maximum mixed layer depths (MLD) compare well with observations, but the spatial MLD pattern is shifted relative to observations. Sea ice extent, volume, and concentration agree well with observations. At high resolution, the sea surface height compares well with satellite observations in mean and variability.

The land ice contribution to global mean sea level rise has not yet been predicted 1 using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models 2–8, but primarily used previous-generation scenarios 9 and climate models 10, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios 11,12 using statistical emulation of the ice sheet and glacier models. We find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.

The Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP) is a community effort sponsored by the Climate and Cryosphere (CliC) project. MISOMIP aims to design and coordinate a series of MIPs—some idealized and realistic—for model evaluation, verification with observations, and future projections for key regions of the West Antarctic Ice Sheet (WAIS). The first phase of the project, MISOMIP1, was an idealized, coupled set of experiments that combined elements from the MISMIP+ and ISOMIP+ standalone experiments for ice-sheet and ocean models, respectively. These MIPs had 3 main goals: 1) to provide simplified experiments that allow model developers to compare their results with those from other models; 2) to suggest a path for testing components in the process of developing a coupled ice sheet-ocean model; and 3) to enable a large variety of parameter and process studies that branch off from these basic experiments.Here, we describe preliminary analysis of the MISOMIP1 results. Eight models in 14 configurations participated in the MIP. In keeping with analysis of the MISMIP+ experiment, we find that the choice of basal friction parameterizations in the ice-sheet component (Weertman vs. Coulomb limited) has a particularly significant impact on the rate of ice-sheet retreat but the choice of stress approximation (SSA, SSA* or L1Lx) seems to have little impact. Models with Coulomb-limited basal friction also tend to be those with the highest melt rates, confirming a positive feedback between melt and retreat in the MISOMIP1 configuration seen in previous work. The ocean component’s treatment of the boundary layer below the ice shelf also has a significant impact on melt rates and resulting retreat, consistent with findings based on ISOMIP+. Feedbacks between the components lead to localized features in the melt rates and the ice geometry not seen in standalone simulations, though the ~2-km horizontal and ~20-m vertical resolution of these simulations appears to be too coarse to produce long-lived, sub-ice-shelf channels seen at higher resolution.

Large-eddy simulations are used to investigate boundary layer turbulence and its control on ice-shelf basal melt rates in Antarctic settings. We present simulations at relatively low thermal driving and low ice-shelf basal slopes, resulting in simulated melt rates from 10s cm/yr to several m/yr. Our results are broadly consistent with the linear relationships between far-field thermal driving and melt rate and between ice-shelf slope and melt rate reported by previous studies. The simulated thermal exchange coefficient is lower than recommended values; thermal exchange becomes less efficient as stratification increases. In our simulations, shear production of turbulent kinetic energy outweighs buoyant production, as found below Larsen C Ice Shelf through recent microstructure measurements. We also find that turbulent intensity and melt rate vary significantly with the orientation between the ice-shelf slope and the far-field flow, even at low ice-shelf slopes. Our results suggest that numerical ocean models employing the standard ice-shelf melt parameterization will underestimate slope effects on ice-shelf melt rates even if they capture the mean buoyancy effects on boundary layer flow. The proposed slope effects would modify feedbacks between ocean circulation and ice-shelf geometry and tidal variability in ice-shelf melt rates.

The processes responsible for freshwater flux from the Antarctic Ice Sheet (AIS) -- ice-shelf basal melting and iceberg calving -- are generally poorly represented in current Earth System Models (ESMs). Here, we document the first effort to date at simulating the ocean circulation and exchanges of heat and freshwater within ice-shelf cavities in a coupled ESM, the Department of Energy's Energy Exascale Earth System Model (E3SM). As a step towards full ice-sheet coupling, we implemented static Antarctic ice-shelf cavities and the ability to calculate ice-shelf basal melt rates from the heat and freshwater fluxes computed by the ocean model. In addition, we added the capability to prescribe forcing from iceberg melt, allowing us to realistically represent the other dominant mass loss process from the AIS. In global, low resolution (i.e., non-eddying ocean) simulations, we find high sensitivity of modeled ocean/ice shelf interactions to the ocean state, which can result in a tipping point to high melt regimes under certain ice shelves, presenting a significant challenge to representing the ocean/ice shelf system in a coupled ESM. We show that inclusion of a spatially dependent parameterization of eddy-induced transport reduces biases in water mass properties on the Antarctic continental shelf. With these improvements, E3SM produces realistic and stable ice-shelf basal melt rates across the continent under pre-industrial climate forcing. We also show preliminary results using an ocean/sea-ice grid that makes use of E3SM’s regional-refinement capability, where increased resolution (down to 12km) is placed in the Southern Ocean around Antarctica, bypassing the need for parameterization of eddy-induced transport in this region. The accurate representation of these processes within a coupled ESM is an important step towards reducing uncertainties in projections of the Antarctic response to climate change and Antarctica's contribution to global sea-level rise.

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law, and to the treatment of stresses and melt rates near the grounding line. We present the effects of these numerical choices on marine ice-sheet dynamics in the Community Ice Sheet Model (CISM). In the experimental framework of the Marine Ice Sheet Model Intercomparison Project (MISMIP+), we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice-sheet geometry, forcing, or model numerics. In the MISMIP+ framework, the ice flow is also sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ∼2 km (arguably 4 km) when the connectivity is low or moderate, and ∼1 km (arguably 2 km) when the connectivity is strong.

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 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 several combinations of sub-ice-shelf melt schemes, CISM is spun up to steady state over many millennia. During the spin-up, basal-friction and thermal-forcing parameters are adjusted to optimize agreement with the observed ice thickness. The model is then run forward to year 2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, ocean warming triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner-Ronne and Ross sectors. The ocean-forced sea-level rise in 2500 varies from about 150 mm to 1300 mm, depending on the melt scheme and ocean forcing applied. Further experiments show relatively high sensitivity to the basal friction law, and moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves. 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.

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.

Projections of the sea level contribution from the Greenland and Antarctic ice sheets (GrIS and AIS) rely on atmospheric and oceanic drivers obtained from climate models. The Earth System Models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6) generally project greater future warming compared with the previous Coupled Model Intercomparison Project phase 5 (CMIP5) effort. Here we use four CMIP6 models and a selection of CMIP5 models to force multiple ice sheet models as part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We find that the projected sea level contribution at 2100 from the ice sheet model ensemble under the CMIP6 scenarios falls within the CMIP5 range for the Antarctic ice sheet but is significantly increased for Greenland. Warmer atmosphere in CMIP6 models results in higher Greenland mass loss due to surface melt. For Antarctica, CMIP6 forcing is similar to CMIP5 and mass gain from increased snowfall counteracts increased loss due to ocean warming.

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ∼2 km, as opposed to ∼0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ∼2 km (arguably 4 km) when the connectivity is low or moderate and ∼1 km (arguably 2 km) when the connectivity is strong.

Sub-ice-shelf melting plays a critical role in the dynamics of the Antarctic Ice Sheet and also feeds back on the regional climate, transforming ocean properties (e.g., affecting deep-water production and sea-ice formation). A full understanding of these processes, as well as the ability to project their response to a changing climate, requires Earth System Models (ESMs) that include coupling with ice-sheet processes. However, biases in regional Antarctic climate can be amplified through sub-ice-shelf melting, and biased melt rates can have significant adverse effects on ice-sheet model initialization and evolution. In preparation for inclusion of dynamic ice sheets in ESMs, this presentation discusses our recent experience in understanding the causes of biases in ocean properties on the Antarctic continental shelf and their relationship to ice-shelf melting. Differences in model behavior across configurations and simulations using the Energy Exascale Earth System Model (E3SM) demonstrates a sensitivity of melt rates to climate. We assess the sensitivity of those melt rates to changes in the region’s climate, including freshening on the continental shelf and shoaling of the thermocline. We also show that ice-shelf meltwater feeds back onto the climate, for example, by affecting melting under neighboring ice shelves, sometimes dramatically so. We demonstrate that significant reductions in melt-rate biases can be achieved through modifications to ocean model mixing parameterizations. This work charts a path forward for configuring ESMs to produce realistic Antarctic melt rates.

Predictions of ice shelf melting depend on dynamical insights into ocean boundary layers below ice shelves. Fundamental questions regarding the nature of stratified turbulence below the sloped and ablating ice shelf base remain. Laboratory experiments, direct numerical simulations, and observations have yielded important insights, but have yet to produce a robust relationship between ice shelf melt rates and shear- and buoyancy-driven mixing. This relationship is the target of our Large-Eddy Simulations (LES) of the ice-shelf ocean boundary layer. Several new developments were applied to the LES code PALM to produce dynamic melting as well as tides. In this presentation, we demonstrate these new model capabilities. We contrast profiles of vertical turbulent fluxes of heat, salt and momentum across different simulated ice shelf settings: cold, shear-dominated settings vs. warm, buoyancy-dominated settings. We also discuss our recent work toward a new ice-shelf melt parameterization for use in large-scale ocean models on the basis of these simulations. A new melt parameterization is a critical component of ongoing ice-ocean coupling efforts, both to place melt rate predictions on a more physical footing and to achieve convergence with vertical ocean model resolution, on which current parameterizations fail.

We seek to understand causal connections between changes in sub-ice shelf melting, ice shelf buttressing, and grounding-line flux. Using a numerical ice flow model, we study changes in ice shelf buttressing and grounding line flux due to localized ice thickness perturbations – a proxy for changes in sub-ice shelf melting – applied to idealized (MISMIP+) and realistic (Larsen C) domains. From our experiments, we identify a correlation between a locally derived buttressing number on the ice shelf, based on the first principal stress, and changes in the integrated grounding line flux. The origin of this correlation, however, remains elusive from a physical perspective; while local thickness perturbations on the ice shelf (thinning) generally correspond to local increases in buttressing, their integrated impact on changes at the grounding line are exactly the opposite (buttressing at the grounding line decreases and ice flux at the grounding line increases). This and additional complications encountered when examining realistic domains motivates us to seek an alternative approach, an adjoint-based method for calculating the sensitivity of the integrated grounding line flux to local changes in ice shelf geometry. We show that the adjoint-based sensitivity is identical to that deduced from pointwise, diagnostic model perturbation experiments. Based on its much wider applicability and the significant computational savings, we propose that the adjoint-based method is ideally suited for assessing grounding line flux sensitivity to changes in sub-ice shelf melting.

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 standalone 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 standalone 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 standalone ice sheet simulations, document the experimental framework and implementation, as well as present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups.

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, ISMIP6 projections of ocean-forced Antarctic Ice Sheet evolution are illustrated using the Community Ice Sheet Model (CISM). 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 nudge the ice thickness toward observed values. The model is then run forward for 500 years, 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, including the Amundsen, Filchner-Ronne, and Ross Basins. Mass loss accelerates late in the 21st century and rises steadily over the next several centuries without leveling off. The resulting ocean-forced SLR at year 2500 varies from about 10 cm to nearly 2 m, depending on the melt scheme and model forcing. Relatively little ice loss is simulated in East Antarctica. Large uncertainties remain, as a result of parameterized basal melt rates, missing ocean and ice sheet physics, and the lack of ice–ocean coupling.

Recent advances in both ocean modeling and melt parameterization in ice-sheet models point the way toward coupled ice sheet–ocean modeling, which is needed to quantify Antarctic mass loss and the resulting sea-level rise. The latest Antarctic ocean modeling shows that complex interactions between the atmosphere, sea ice, icebergs, bathymetric features, and ocean circulation on many scales determine which water masses reach ice-shelf cavities and how much heat is available to melt ice. Meanwhile, parameterizations of basal melting in standalone ice-sheet models have evolved from simplified, depth-dependent functions to more sophisticated models, accounting for ice-shelf basal topography, and the evolution of the sub-ice-shelf buoyant flow. The focus of recent work has been on better understanding processes or adding new model capabilities, but a broader community effort is needed in validating models against observations and producing melt-rate projections. Given time, community efforts in coupled ice sheet–ocean modeling, already underway, will tackle the considerable challenges involved in building, initializing, constraining, and performing projections with coupled models, leading to reduced uncertainties in Antarctica’s contribution to future sea-level rise.

Understanding ice shelf–ocean interaction is fundamental to projecting the Antarctic ice sheet response to a warming climate. Numerical ice shelf–ocean models are a powerful tool for simulating this interaction, yet are limited by inherent model weaknesses and scarce observations, leading to parameterisations that are unverified and unvalidated below ice shelves. We explore how different models simulate ice shelf–ocean interaction using the 2nd Ice Shelf–Ocean Model Intercomparison Project (ISOMIP+) framework. Vertical discretisation and resolution of the ocean model are shown to have a significant effect on ice shelf basal melt rate, through differences in the distribution of meltwater fluxes and the calculation of thermal driving. Z-coordinate models, which generally have coarser vertical resolution in ice shelf cavities, may simulate higher melt rates compared to terrain-following coordinate models. This is due to the typically higher resolution of the ice–ocean boundary layer region in terrain following models, which allows better representation of a thin meltwater layer, increased stratification, and as a result, better insulation of the ice from water below. We show that a terrain-following model, a z-level coordinate model and a hybrid approach give similar results when the effective vertical resolution adjacent to the ice shelf base is similar, despite each model employing different paradigms for distributing meltwater fluxes and sampling tracers for melting. We provide a benchmark for thermodynamic ice shelf–ocean interaction with different model vertical coordinates and vertical resolutions, and suggest a framework for any future ice shelf–ocean thermodynamic parameterisations.

Ocean-induced melting below ice shelves is one of the dominant drivers for mass loss from the Antarctic Ice Sheet at present. An appropriate representation of sub-shelf melt rates is therefore essential for model simulations of marine-based ice sheet evolution. Continental-scale ice sheet models often rely on simple melt-parameterizations, in particular for long-term simulations, when fully coupled ice–ocean interaction becomes computationally too expensive. Such parameterizations can account for the influence of the local depth of the ice-shelf draft or its slope on melting. However, they do not capture the effect of ocean circulation underneath the ice shelf. Here we present the Potsdam Ice-shelf Cavity mOdel (PICO), which simulates the vertical overturning circulation in ice-shelf cavities and thus enables the computation of sub-shelf melt rates consistent with this circulation. PICO is based on an ocean box model that coarsely resolves ice shelf cavities and uses a boundary layer melt formulation. We implement it as a module of the Parallel Ice Sheet Model (PISM) and evaluate its performance under present-day conditions of the Southern Ocean. We identify a set of parameters that yield two-dimensional melt rate fields that qualitatively reproduce the typical pattern of comparably high melting near the grounding line and lower melting or refreezing towards the calving front. PICO captures the wide range of melt rates observed for Antarctic ice shelves, with an average of about 0.1 m a−1 for cold sub-shelf cavities, for example, underneath Ross or Ronne ice shelves, to 16 m a−1 for warm cavities such as in the Amundsen Sea region. This makes PICO a computationally feasible and more physical alternative to melt parameterizations purely based on ice draft geometry.