A comparison of two Stokes ice sheet models applied to the Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d)
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 planview models (MISMIP3D). For the diagnostic experiment (P75D) the two models give very 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 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 largely due to slightly different choices in the implementation of basal boundary conditions used by the two models. Based on current understanding, neither set of implementations is more or less favorable, in which case we propose that the span of different grounding line positions from these two models provides a measure of uncertainty when treating the results from full-Stokes models as a metric for accuracy in model intercomparison experiments. More importantly, we show that as grid resolution increases the grounding line positions for FELIX-S and Elmer/Ice appear to converge. We conclude that future model intercomparisons using full-Stokes models as a metric should include more than one model, to provide both additional confidence in the results from full-Stokes models and a measure of their uncertainty.