On the discretization of the ice thickness distribution in the NEMO3.6-LIM3 global ocean–sea ice model

• Published in: Geoscientific Model Development
• Main Authors: Massonnet, François, Barthélemy, Antoine, Worou, Koffi, Fichefet, Thierry, Vancoppenolle, Martin, rousset, clement, Moreno-Chamarro, Eduardo
• Other Authors: Centre Georges Lemaître for Earth and Climate Research Louvain (TECLIM), Earth and Life Institute Louvain-La-Neuve (ELI), Université Catholique de Louvain = Catholic University of Louvain (UCL)-Université Catholique de Louvain = Catholic University of Louvain (UCL), Nucleus for European Modeling of the Ocean (NEMO R&D ), Laboratoire d'Océanographie et du Climat : Expérimentations et Approches Numériques (LOCEAN), Muséum national d'Histoire naturelle (MNHN)-Institut de Recherche pour le Développement (IRD)-Institut national des sciences de l'Univers (INSU - CNRS)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS)-Institut Pierre-Simon-Laplace (IPSL (FR_636)), École normale supérieure - Paris (ENS-PSL), Université Paris Sciences et Lettres (PSL)-Université Paris Sciences et Lettres (PSL)-Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Institut national des sciences de l'Univers (INSU - CNRS)-École polytechnique (X), Institut Polytechnique de Paris (IP Paris)-Institut Polytechnique de Paris (IP Paris)-Centre National d'Études Spatiales Toulouse (CNES)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS)-Université Paris Cité (UPCité)-École normale supérieure - Paris (ENS-PSL), Université Paris Sciences et Lettres (PSL)-Université Paris Sciences et Lettres (PSL)-Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-École polytechnique (X), Institut Polytechnique de Paris (IP Paris)-Institut Polytechnique de Paris (IP Paris)-Centre National d'Études Spatiales Toulouse (CNES)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS)-Université Paris Cité (UPCité)-Muséum national d'Histoire naturelle (MNHN)-Institut de Recherche pour le Développement (IRD)-Institut national des sciences de l'Univers (INSU - CNRS)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS)-Institut Pierre-Simon-Laplace (IPSL (FR_636)), Institut Polytechnique de Paris (IP Paris)-Institut Polytechnique de Paris (IP Paris)-Centre National d'Études Spatiales Toulouse (CNES)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS)-Université Paris Cité (UPCité), Barcelona Supercomputing Center - Centro Nacional de Supercomputacion (BSC-CNS), European Commission (APPLICATE (grant no. 727862), European Project: 641727,H2020,H2020-SC5-2014-two-stage,PRIMAVERA(2015)
• Format: Article in Journal/Newspaper
• Language: English
• Published: HAL CCSD 2019
• Subjects: [SDU.STU.OC]Sciences of the Universe [physics]/Earth Sciences/Oceanography; [SDU.STU.GL]Sciences of the Universe [physics]/Earth Sciences/Glaciology; Antarc*; Antarctic; Sea ice
• Online Access: https://hal.sorbonne-universite.fr/hal-02286889; https://hal.sorbonne-universite.fr/hal-02286889/document; https://hal.sorbonne-universite.fr/hal-02286889/file/gmd-12-3745-2019.pdf; https://doi.org/10.5194/gmd-12-3745-2019

International audience The ice thickness distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to discretize it remains an open question. Here, we use the ocean-sea ice general circulation model, Nucleus for European Modelling of the Ocean (NEMO) version 3.6 and Louvain-la-Neuve sea Ice Model (LIM) version 3 (NEMO3.6-LIM3), forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ∼ 30 % and ∼ 10 % of the reference values (run with five categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ∼ 4 and ∼ 2 m for the Arctic and Antarctic , respectively, is a prerequisite for allowing the storage of deformed ice and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean sea ice state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories , has a primary ...