U.S. Department of Energy Office of Biological and Environmental Research

BER Research Highlights

Adaptive Mesh Refinement Versus Subgrid Friction Interpolation in Simulations of Antarctic Ice Dynamics
Published: May 13, 2016
Posted: May 10, 2017

Dynamic-adaptive mesh ice sheet shows importance of mesh spacing to capture instability of ice-sheets.

The Science
The 2007 & 2013 Intergovernmental Panel on Climate Change (IPCC) reports highlighted the need for better projections of the Antarctic contribution to sea-level rise (SLR). The West Antarctic Ice Sheet (WAIS) is vulnerable to marine forcing from warm-water incursion into subshelf cavities, causing melting and thinning of ice shelves and weakening or eliminating the buttressing exerted on feeder ice streams.  The result is accelerated ice loss and increased SLR. Theory and modeling suggest that ice-sheet models require very fine (sub-km) resolution to correctly model these effects. This study uses extreme melt forcing scenario to examine resolution requirements for millennial-scale whole-continent simulations.

The Impact
This is the first fully-resolved whole-continent millennial-scale Antarctic simulations to evaluate resolution requirements and it is broadly useful as a method to ensure confidence in model results. The Berkeley Adaptive Mesh Refinement (AMR) ice sheet model (BISICLES) fully resolves dynamically important regions like grounding lines. The use of a basic subgrid-scale friction interpolation scheme results in factor-of-two reduction in resolution requirements. Under-resolved-mesh simulations consistently underestimate SLR contributions. Even moderate under-resolution results in erroneous predictions of timing and mechanisms of WAIS collapse, while severe under-resolution results in missing WAIS collapse entirely. This is the first application of BISICLES, a state-of-the-art AMR ice sheet model, to millennial-scale whole-continent projections of Antarctic response to extreme climate forcing. AMR allows for full resolution of dynamically changing grounding lines in regional- and continental-scale simulations.

At least in conventional hydrostatic ice-sheet models, the numerical error associated with grounding line dynamics can be reduced by modifications to the discretization scheme. These involve altering the integration formulae for the basal traction and/or driving stress close to the grounding line and exhibit lower — if still first-order — error in the MISMIP3d experiments. MISMIP3d may not represent the variety of real ice streams, in that it lacks strong lateral stresses, and imposes a large basal traction at the grounding line. We study resolution sensitivity in the context of extreme forcing simulations of the entire Antarctic ice sheet, using the BISICLES adaptive mesh ice-sheet model with two schemes: the original treatment, and a scheme, which modifies the discretization of the basal traction. The second scheme does indeed improve accuracy — by around a factor of two — for a given mesh spacing, but < 1 km resolution is still necessary. For example, in coarser resolution simulations Thwaites Glacier retreats so slowly that other ice streams divert its trunk. In contrast, with < 1 km meshes, the same glacier retreats far more quickly and triggers the final phase of West Antarctic collapse a century before any such diversion can take place.

Contacts (BER PM)
Dorothy Koch
Dorothy.Koch@doe.science.gov (301-903-0105)

(PI Contact)
Dan Martin
Lawrence Berkeley National Lab

Earth System Modeling, SciDAC PISCEES project


 “Adaptive mesh refinement versus Subgridfriction interpolation in simulations of Antarctic Ice Dynamics”, Cornford, Martin, Lee, Payne, Ng, Annals of Glaciology, 57 (73), DOI: 10.1017/aog.2016.13, 2016 (Reference link)

Topic Areas:

  • Research Area: Earth and Environmental Systems Modeling

Division: SC-23.1 Climate and Environmental Sciences Division, BER


BER supports basic research and scientific user facilities to advance DOE missions in energy and environment. More about BER

Recent Highlights

Aug 24, 2019
New Approach for Studying How Microbes Influence Their Environment
A diverse group of scientists suggests a common framework and targeting of known microbial processes [more...]

Aug 08, 2019
Nutrient-Hungry Peatland Microbes Reduce Carbon Loss Under Warmer Conditions
Enzyme production in peatlands reduces carbon lost to respiration under future high temperatures. [more...]

Aug 05, 2019
Amazon Forest Response to CO2 Fertilization Dependent on Plant Phosphorus Acquisition
AmazonFACE Model Intercomparison. The Science Plant growth is dependent on the availabi [more...]

Jul 29, 2019
A Slippery Slope: Soil Carbon Destabilization
Carbon gain or loss depends on the balance between competing biological, chemical, and physical reac [more...]

Jul 15, 2019
Field Evaluation of Gas Analyzers for Measuring Ecosystem Fluxes
How gas analyzer type and correction method impact measured fluxes. The Science A side- [more...]

List all highlights (possible long download time)