Journal metrics

Journal metrics

  • IF value: 4.252 IF 4.252
  • IF 5-year value: 4.890 IF 5-year 4.890
  • CiteScore value: 4.49 CiteScore 4.49
  • SNIP value: 1.539 SNIP 1.539
  • SJR value: 2.404 SJR 2.404
  • IPP value: 4.28 IPP 4.28
  • h5-index value: 40 h5-index 40
  • Scimago H index value: 51 Scimago H index 51
Volume 11, issue 9 | Copyright
Geosci. Model Dev., 11, 3747-3780, 2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Model description paper 18 Sep 2018

Model description paper | 18 Sep 2018

MPAS-Albany Land Ice (MALI): a variable-resolution ice sheet model for Earth system modeling using Voronoi grids

Matthew J. Hoffman1, Mauro Perego2, Stephen F. Price1, William H. Lipscomb1,3, Tong Zhang1, Douglas Jacobsen1, Irina Tezaur4, Andrew G. Salinger2, Raymond Tuminaro2, and Luca Bertagna2 Matthew J. Hoffman et al.
  • 1Fluid Dynamics and Solid Mechanics Group, Los Alamos National Laboratory, P.O. Box 1663, MS B216, Los Alamos, NM, USA
  • 2Center for Computing Research, Sandia National Laboratories, P.O. Box 5800, MS 1320, Albuquerque, NM, USA
  • 3National Center for Atmospheric Research, Boulder, CO, USA
  • 4Quantitative Modeling and Analysis, Sandia National Laboratories, P.O. Box 969, MS 9159, Livermore, CA, USA

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.

Download & links
Publications Copernicus
Short summary
MPAS-Albany Land Ice (MALI) is a new variable-resolution land ice model that uses unstructured grids on a plane or sphere. MALI is built for Earth system modeling on high-performance computing platforms using existing software libraries. MALI simulates the evolution of ice thickness, velocity, and temperature, and it includes schemes for simulating iceberg calving and the flow of water beneath ice sheets and its effect on ice sliding. The model is demonstrated for the Antarctic ice sheet.
MPAS-Albany Land Ice (MALI) is a new variable-resolution land ice model that uses unstructured...