Marine ice sheet dynamics: Hysteresis and neutral equilibrium

Abstract : The stability of marine ice sheets and outlet glaciers is mostly controlled by the dynamics of their grounding line, i.e., where the bottom contact of the ice changes from bedrock or till to ocean water. The last report of the Intergovernmental Panel on Climate Change has clearly underlined the poor ability of models to capture the dynamics of outlet glaciers. Here we present computations of grounding line dynamics on the basis of numerical solutions of the full Stokes equations for ice velocity, coupled with the evolution of the air ice– and sea ice–free interfaces. The grounding line position is determined by solving the contact problem between the ice and a rigid bedrock using the finite element code Elmer. Results of the simulations show that marine ice sheets are unstable on upsloping beds and undergo hysteresis under perturbation of ice viscosity, confirming conclusions from boundary layer theory. The present approach also indicates that a 2-D unconfined marine ice sheet sliding over a downsloping bedrock does not exhibit neutral equilibrium. It is shown that mesh resolution around the grounding line is a crucial issue. A very fine grid size (<100 m spacing) is needed in order to achieve consistent results.