A rapidly converging initialisation method to simulate the present-day Greenland ice sheet using the GRISLI ice sheet model (version 1.3)

• Published in: Geoscientific Model Development
• Main Authors: Le Clec'H, Sébastien, Quiquet, Aurélien, Charbit, Sylvie, Dumas, Christophe, Kageyama, Masa, Ritz, Catherine
• Other Authors: Laboratoire des Sciences du Climat et de l'Environnement Gif-sur-Yvette (LSCE), Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS)-Direction de Recherche Fondamentale (CEA) (DRF (CEA)), Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA), Modélisation du climat (CLIM), Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS)-Direction de Recherche Fondamentale (CEA) (DRF (CEA)), Institut des Géosciences de l’Environnement (IGE), Institut de Recherche pour le Développement (IRD)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National de la Recherche Scientifique (CNRS)-Université Grenoble Alpes 2016-2019 (UGA 2016-2019 )-Observatoire des Sciences de l'Univers de Grenoble (Fédération OSUG), ANR-15-CE01-0015,AC-AHC2,Circulation atmosphérique et changement de cycle hydrologique pour l'Arctique(2015)
• Format: Article in Journal/Newspaper
• Language: English
• Published: CCSD 2019
• Subjects: [SDU.STU.CL]Sciences of the Universe [physics]/Earth Sciences/Climatology; [SDU.STU.GL]Sciences of the Universe [physics]/Earth Sciences/Glaciology; Greenland; Ice Sheet; Ice Shelf
• Online Access: https://hal.science/hal-02328234; https://hal.science/hal-02328234v1/document; https://hal.science/hal-02328234v1/file/gmd-12-2481-2019.pdf; https://doi.org/10.5194/gmd-12-2481-2019

International audience Providing reliable projections of the ice sheet contribution to future sea-level rise has become one of the main challenges of the ice sheet modelling community. To increase confidence in future projections, a good knowledge of the present-day state of ice flow dynamics, which is critically dependent on basal conditions, is strongly needed. The main difficulty is tied to the scarcity of observations at the ice-bed interface at the scale of the whole ice sheet, resulting in poorly constrained parameterisations in ice sheet models. To circumvent this drawback, inverse modelling approaches can be developed to infer initial conditions for ice sheet models that best reproduce available data. Most often such approaches allow for a good representation of the mean present-day state of the ice sheet but are accompanied with unphys-ical trends. Here, we present an initialisation method for the Greenland ice sheet using the thermo-mechanical hybrid GRISLI (GRenoble Ice Shelf and Land Ice) ice sheet model. Our approach is based on the adjustment of the basal drag coefficient that relates the sliding velocities at the ice-bed interface to basal shear stress in unfrozen bed areas. This method relies on an iterative process in which the basal drag is periodically adjusted in such a way that the simulated ice thickness matches the observed one. The quality of the method is assessed by computing the root mean square errors in ice thickness changes. Because the method is based on an adjustment of the sliding velocities only, the results are discussed in terms of varying ice flow enhancement factors that control the deformation rates. We show that this factor has a strong impact on the minimisation of ice thickness errors and has to be chosen as a function of the internal thermal state of the ice sheet (e.g. a low enhancement factor for a warm ice sheet). While the method performance slightly increases with the duration of the minimisation procedure, an ice thickness root mean square error (RMSE) of 50.3 m ...