EGU26-12617, updated on 14 Mar 2026
https://doi.org/10.5194/egusphere-egu26-12617
EGU General Assembly 2026
© Author(s) 2026. This work is distributed under
the Creative Commons Attribution 4.0 License.
Oral | Monday, 04 May, 17:20–17:30 (CEST)
 
Room L2
IPSL-CM-Elmer/Ice: a new coupled ice sheet – climate model
Lucas Bastien1, Pierre Mathiot1, Nicolas C. Jourdain1, Cécile Agosta2, Justine Caillet3, Arnaud Caubel2, Sylvie Charbit2, Mondher Chekki1, Julie Deshayes4, Christophe Dumas2, Gaël Durand1, Fabien Gillet-Chaulet1, Olivier Marti2, Cyrille Mosbeux1, and Etienne Vignon5
Lucas Bastien et al.
  • 1Univ. Grenoble Alpes/CNRS/INRAE/IRD/Grenoble INP, IGE, Grenoble, France
  • 2Laboratoire des Sciences du Climat et de l'Environnement/IPSL/CEA/CNRS/UPS, Gif-sur-Yvette/Paris, France
  • 3Thayer School of Engineering, Dartmouth College, Hanover, NH, USA
  • 4LOCEAN Laboratory/Sorbonne Université/CNRS/IRD/MNHN, Paris, France
  • 5Laboratoire de Météorologie Dynamique/IPSL/Sorbone Universités/CNRS, Paris, France

The contribution of ice sheets to future sea level rise remains highly uncertain, and complex positive feedback mechanisms can lead to accelerating melt in a warming climate. Yet, few climate models explicitly represent ice flow of Greenland and Antarctica, or their interactions with the rest of the climate system.

Here we present the coupling of the Elmer/Ice ice sheet model with the IPSL-CM7 climate model. Two-way coupling with the atmospheric and oceanic components of IPSL-CM7 (LMDZ and NEMO, respectively) occurs every simulated year. On the atmospheric side, the surface mass balance from LMDZ is used to force the ice sheet model. In this coupling step, a positive degree day scheme is used to re-calculate surface melt and runoff for Greenland to yield more realistic results. The elevation of the LMDZ domain’s bottom surface is in turn updated to account for the new ice sheet geometry provided by Elmer/Ice. On the ocean side, sub-ice shelf melting is explicitly represented where NEMO's resolution allows it and is extrapolated near the grounding line and under small ice shelves, where the cavity geometry is not resolved by the ocean model. NEMO’s computational domain is updated yearly to account for new icy or wet cells.

We then present the results of two 100-year simulations, which were conducted to test the robustness of the coupling and the behaviour of the model in a warming climate. The first  simulation has a constant pre-industrial atmospheric CO2 concentration, whereas in the second one  the CO2 concentration increases by 1% every year. We describe some interesting features that emerge due to increasing CO2 concentrations, such as the transition from cold to warm water on the continental shelf of the Amundsen Sea, and a retreat of the grounding line in this region.

While still in its early stage of development, this work is an important milestone in the addition of interactive ice sheets within the IPSL-CM7 climate model. Future developments include interactive ice fronts, which are currently fixed in the model, and the possibility of uncovering solid ground as ice sheets retreat.

How to cite: Bastien, L., Mathiot, P., Jourdain, N. C., Agosta, C., Caillet, J., Caubel, A., Charbit, S., Chekki, M., Deshayes, J., Dumas, C., Durand, G., Gillet-Chaulet, F., Marti, O., Mosbeux, C., and Vignon, E.: IPSL-CM-Elmer/Ice: a new coupled ice sheet – climate model, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-12617, https://doi.org/10.5194/egusphere-egu26-12617, 2026.