Skip to Main content Skip to Navigation
New interface
Conference papers

Full-Stokes finite element modelling of the Greenland ice-sheet using inverse methods.

Fabien Gillet-Chaulet 1, * Olivier Gagliardini 1 Maëlle Nodet 2 Catherine Ritz 1 Gaël Durand 1 Thomas Zwinger 3 Hakime Seddik 4 Ralph Greve 4 
* Corresponding author
2 MOISE - Modelling, Observations, Identification for Environmental Sciences
Inria Grenoble - Rhône-Alpes, LJK - Laboratoire Jean Kuntzmann, Grenoble INP - Institut polytechnique de Grenoble - Grenoble Institute of Technology
Abstract : The contribution of the Greenland ice sheet to sea level rise has continuously increased during the last decade and is a concern for the next centuries. Within the two international efforts ice2sea and SeaRise, many groups are now running prognostic simulations of the Greenland ice sheet. The quality of the model predictions depend primary on the good description of the physical processes involved and on a good initial state reproducing the current observations (geometry and surface velocities). But, these two points remain weaknesses of the current ice flow models: (i) The mass loss is equally split between surface processes and ice dynamics (Van den Broeke et al., Science, 2009). For the dynamical part, the ice discharge is localized on small outlet glaciers and is impossible to capture in large scale ice flow models using shallow approximations. (ii) The quantity and quality of the observations has considerably increased during the last decades, but data assimilation techniques are still under-employed during the spin-up procedures. Here, we model the flow of ice of the whole Greenland ice sheet using the full-Stokes finite element code Elmer/Ice. The finite element mesh is generated using the anisotropic mesh adaptation tool YAMS. The resulting mesh size varies from few tens of kilometers in the central parts, down to a 1km in the main outlet glaciers. Two inverse methods are used to initialize the basal sliding coefficient field to reduce the mismatch between the present day surface velocities and the model surface velocities:(i) a Robin Inverse problem recently proposed by Arthern and Gudmundsson (J. Glaciol., 2010), and (ii) a control method using the fact that the Stokes equations are self-adjoint for a Newtonian rheology (Morlighem et al., GRL, 2010). We show that magnitude and repartition of the modelled present day ice discharge compares well with the observations. Finally, transient forecast simulations of the Greenland ice sheet for the 21st century are run to investigate the model sensitivity to the initialization procedure.
Complete list of metadata
Contributor : Maëlle Nodet Connect in order to contact the contributor
Submitted on : Tuesday, April 19, 2011 - 8:16:52 PM
Last modification on : Tuesday, July 5, 2022 - 8:38:21 AM


  • HAL Id : inria-00587251, version 1


Fabien Gillet-Chaulet, Olivier Gagliardini, Maëlle Nodet, Catherine Ritz, Gaël Durand, et al.. Full-Stokes finite element modelling of the Greenland ice-sheet using inverse methods.. EGU General Assembly, Apr 2011, Vienne, Austria. ⟨inria-00587251⟩



Record views