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

May 10, 2019
Quantifying Decision Uncertainty in Water Management via a Coupled Agent-Based Model
Considering risk perception can improve the representation of human decision-making processes in age [more...]

May 09, 2019
Projecting Global Urban Area Growth Through 2100 Based on Historical Time Series Data and Future Scenarios
Study provides country-specific urban area growth models and the first dataset on country-level urba [more...]

May 05, 2019
Calibrating Building Energy Demand Models to Refine Long-Term Energy Planning
A new, flexible calibration approach improved model accuracy in capturing year-to-year changes in bu [more...]

May 03, 2019
Calibration and Uncertainty Analysis of Demeter for Better Downscaling of Global Land Use and Land Cover Projections
Researchers improved the Demeter model’s performance by calibrating key parameters and establi [more...]

Apr 22, 2019
Representation of U.S. Warm Temperature Extremes in Global Climate Model Ensembles
Representation of warm temperature events varies considerably among global climate models, which has [more...]

List all highlights (possible long download time)