MPAS-Albany Land Ice (MALI): A Variable-Resolution Ice Sheet Model for Earth System Modeling Using Voronoi Grids

TitleMPAS-Albany Land Ice (MALI): A Variable-Resolution Ice Sheet Model for Earth System Modeling Using Voronoi Grids
Publication TypeJournal Article
Year of Publication2018
AuthorsHoffman, Matthew J., Perego Mauro, Price Stephen F., Lipscomb William H., Zhang Tong, Jacobsen Douglas, Tezaur Irina, Salinger Andrew G., Tuminaro Raymond, and Bertagna Luca
JournalGeoscientific Model Development
Volume11
Number9
Pages3747-3780
Date Published09/2018
Abstract

We introduce MPAS-Albany Land Ice (MALI) v6.0, a new variable-resolution land ice model that uses unstructured Voronoi grids on a plane or sphere. MALI is built using the Model for Prediction Across Scales (MPAS) framework for developing variable-resolution Earth system model components and the Albany multi-physics code base for the solution of coupled systems of partial differential equations, which itself makes use of Trilinos solver libraries. MALI includes a three-dimensional first-order momentum balance solver (Blatter–Pattyn) by linking to the Albany-LI ice sheet velocity solver and an explicit shallow ice velocity solver. The evolution of ice geometry and tracers is handled through an explicit first-order horizontal advection scheme with vertical remapping. The evolution of ice temperature is treated using operator splitting of vertical diffusion and horizontal advection and can be configured to use either a temperature or enthalpy formulation. MALI includes a mass-conserving subglacial hydrology model that supports distributed and/or channelized drainage and can optionally be coupled to ice dynamics. Options for calving include eigencalving, which assumes that the calving rate is proportional to extensional strain rates. MALI is evaluated against commonly used exact solutions and community benchmark experiments and shows the expected accuracy. Results for the MISMIP3d benchmark experiments with MALI's Blatter–Pattyn solver fall between published results from Stokes and L1L2 models as expected. We use the model to simulate a semi-realistic Antarctic ice sheet problem following the initMIP protocol and using 2km resolution in marine ice sheet regions. MALI is the glacier component of the Energy Exascale Earth System Model (E3SM) version 1, and we describe current and planned coupling to other E3SM components.

URLhttp://dx.doi.org/10.5194/gmd-11-3747-2018
DOI10.5194/gmd-11-3747-2018
Funding Program: 
Journal: Geoscientific Model Development
Number: 9
Volume: 11

We introduce MPAS-Albany Land Ice (MALI) v6.0, a new variable-resolution land ice model that uses unstructured Voronoi grids on a plane or sphere. MALI is built using the Model for Prediction Across Scales (MPAS) framework for developing variable-resolution Earth system model components and the Albany multi-physics code base for the solution of coupled systems of partial differential equations, which itself makes use of Trilinos solver libraries. MALI includes a three-dimensional first-order momentum balance solver (Blatter–Pattyn) by linking to the Albany-LI ice sheet velocity solver and an explicit shallow ice velocity solver. The evolution of ice geometry and tracers is handled through an explicit first-order horizontal advection scheme with vertical remapping. The evolution of ice temperature is treated using operator splitting of vertical diffusion and horizontal advection and can be configured to use either a temperature or enthalpy formulation. MALI includes a mass-conserving subglacial hydrology model that supports distributed and/or channelized drainage and can optionally be coupled to ice dynamics. Options for calving include eigencalving, which assumes that the calving rate is proportional to extensional strain rates. MALI is evaluated against commonly used exact solutions and community benchmark experiments and shows the expected accuracy. Results for the MISMIP3d benchmark experiments with MALI's Blatter–Pattyn solver fall between published results from Stokes and L1L2 models as expected. We use the model to simulate a semi-realistic Antarctic ice sheet problem following the initMIP protocol and using 2km resolution in marine ice sheet regions. MALI is the glacier component of the Energy Exascale Earth System Model (E3SM) version 1, and we describe current and planned coupling to other E3SM components.

DOI: 10.5194/gmd-11-3747-2018
Year of Publication: 2018
Citation:
Hoffman, M, M Perego, S Price, W Lipscomb, T Zhang, D Jacobsen, I Tezaur, A Salinger, R Tuminaro, and L Bertagna.  2018.  "MPAS-Albany Land Ice (MALI): A Variable-Resolution Ice Sheet Model for Earth System Modeling Using Voronoi Grids."  Geoscientific Model Development 11(9): 3747-3780, doi:10.5194/gmd-11-3747-2018.