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

International audience 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 Greenlan...

Full description

Bibliographic Details
Main Authors: Gillet-Chaulet, Fabien, Gagliardini, Olivier, Nodet, Maëlle, Ritz, Catherine, Durand, Gaël, Zwinger, Thomas, Seddik, Hakime, Greve, Ralph
Other Authors: Laboratoire de glaciologie et géophysique de l'environnement (LGGE), Observatoire des Sciences de l'Univers de Grenoble (OSUG), Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB Université de Savoie Université de Chambéry )-Centre National de la Recherche Scientifique (CNRS)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB Université de Savoie Université de Chambéry )-Centre National de la Recherche Scientifique (CNRS)-Centre National de la Recherche Scientifique (CNRS), Modelling, Observations, Identification for Environmental Sciences (MOISE), Inria Grenoble - Rhône-Alpes, Institut National de Recherche en Informatique et en Automatique (Inria)-Institut National de Recherche en Informatique et en Automatique (Inria)-Laboratoire Jean Kuntzmann (LJK), Université Pierre Mendès France - Grenoble 2 (UPMF)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Centre National de la Recherche Scientifique (CNRS)-Université Pierre Mendès France - Grenoble 2 (UPMF)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Centre National de la Recherche Scientifique (CNRS), Scientific Computing Ltd (CSC), Institute of Low Temperature Science Sapporo, Hokkaido University Sapporo, Japan
Format: Conference Object
Language:English
Published: HAL CCSD 2011
Subjects:
Online Access:https://inria.hal.science/inria-00587251
Description
Summary:International audience 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 ...