Full Stokes modeling of marine ice sheets: influence of the grid size

Using the finite-element code Elmer, we show that the full Stokes modeling of the ice-sheet/ice-shelf transition we propose can give consistent predictions of grounding-line migration. Like other marine ice-sheet models our approach is highly sensitive to the chosen mesh resolution. However, with a grid size down to <5 km in the vicinity of the grounding line, predictions start to be robust because: (1) whatever the grid size (<5 km) the steady-state grounding-line position is sensibly the same (6 km standard deviation), and (2) with a grid-size refinement in the vicinity of the grounding line (200m), the steady-state solution is independent of the applied perturbation in fluidity, provided this perturbation remains monotonic.