Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
The Cryosphere, 11, 179-190, 2017
https://doi.org/10.5194/tc-11-179-2017
© Author(s) 2017. This work is distributed under
the Creative Commons Attribution 3.0 License.
Research article
25 Jan 2017
A comparison of two Stokes ice sheet models applied to the Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d)
Tong Zhang1,3,4, Stephen Price2, Lili Ju3, Wei Leng5, Julien Brondex6, Gaël Durand6, and Olivier Gagliardini6 1State Key Laboratory of Severe Weather (LASW), Chinese Academy of Meteorological Sciences, Beijing, China
2Fluid Dynamics and Solid Mechanics Group, Los Alamos National Laboratory, Los Alamos, NM, USA
3Department of Mathematics and Interdisciplinary Mathematics Institute, University of South Carolina, Columbia, SC, USA
4State Key Laboratory of Cryospheric Sciences, Chinese Academy of Sciences, Lanzhou, China
5State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences, Beijing, China
6Université Grenoble Alpes, CNRS, IRD, IGE, 38000 Grenoble, France
Abstract. We present a comparison of the numerics and simulation results for two "full" Stokes ice sheet models, FELIX-S (Leng et al. 2012) and Elmer/Ice (Gagliardini et al. 2013). The models are applied to the Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d). For the diagnostic experiment (P75D) the two models give similar results ( <  2 % difference with respect to along-flow velocities) when using identical geometries and computational meshes, which we interpret as an indication of inherent consistencies and similarities between the two models. For the standard (Stnd), P75S, and P75R prognostic experiments, we find that FELIX-S (Elmer/Ice) grounding lines are relatively more retreated (advanced), results that are consistent with minor differences observed in the diagnostic experiment results and that we show to be due to different choices in the implementation of basal boundary conditions in the two models. While we are not able to argue for the relative favorability of either implementation, we do show that these differences decrease with increasing horizontal (i.e., both along- and across-flow) grid resolution and that grounding-line positions for FELIX-S and Elmer/Ice converge to within the estimated truncation error for Elmer/Ice. Stokes model solutions are often treated as an accuracy metric in model intercomparison experiments, but computational cost may not always allow for the use of model resolution within the regime of asymptotic convergence. In this case, we propose that an alternative estimate for the uncertainty in the grounding-line position is the span of grounding-line positions predicted by multiple Stokes models.

Citation: Zhang, T., Price, S., Ju, L., Leng, W., Brondex, J., Durand, G., and Gagliardini, O.: A comparison of two Stokes ice sheet models applied to the Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d), The Cryosphere, 11, 179-190, https://doi.org/10.5194/tc-11-179-2017, 2017.
Publications Copernicus
Download
Short summary
Stokes-flow models are the highest-fidelity representation of the equations governing ice sheet flow and they are often treated as the standard against which other models are compared in model benchmark activities. We compare two different Stokes models applied to a canonical set of idealized marine ice sheet experiments and demonstrate that the solutions converge with increasing grid resolution. This provides confidence in the use of Stokes models for generating test case solution metrics.
Stokes-flow models are the highest-fidelity representation of the equations governing ice sheet...
Share