[1] 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.

Abstract. In response to climate change, most glaciers are losing mass and hence contribute to sea-level rise. Repeated and accurate mapping of their surface topography is required to estimate their mass balance and to extrapolate/calibrate sparse field glaciological measurements. In this study we evaluate the potential of sub-meter stereo imagery from the recently launched Pléiades satellites to derive digital elevation models (DEMs) of glaciers and their elevation changes. Our five evaluation sites, where nearly simultaneous field measurements were collected, are located in Iceland, the European Alps, the central Andes, Nepal and Antarctica. For Iceland, the Pléiades DEM is also compared to a lidar DEM. The vertical biases of the Pléiades DEMs are less than 1 m if ground control points (GCPs) are used, but reach up to 7 m without GCPs. Even without GCPs, vertical biases can be reduced to a few decimetres by horizontal and vertical co-registration of the DEMs to reference altimetric data on ice-free terrain. Around these biases, the vertical precision of the Pléiades DEMs is ±1 m and even ±0.5 m on the flat glacier tongues (1σ confidence level). Similar precision levels are obtained in the accumulation areas of glaciers and in Antarctica. We also demonstrate the high potential of Pléiades DEMs for measuring seasonal, annual and multiannual elevation changes with an accuracy of 1 m or bet-