# Benefits of the Brinkman Volume Penalisation Method for the Ice-Shelf Melt Rates Produced by Z-coordinate Ocean Models

Antoine-Alexis Nasser, Nicolas C. Jourdain, Pierre Mathiot, Gurvan Madec

► To cite this version:

Antoine-Alexis Nasser, Nicolas C. Jourdain, Pierre Mathiot, Gurvan Madec. Benefits of the Brinkman Volume Penalisation Method for the Ice-Shelf Melt Rates Produced by Z-coordinate Ocean Models. EGU24 - European Geosciences Union General Assembly, EGU, Apr 2024, Vienne, Austria. 10.5194/egusphere-egu24-17676. insu-04725899

HAL Id: insu-04725899

https://insu.hal.science/insu-04725899v1

Submitted on 12 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Benefits of the Brinkman Volume Penalisation Method for the Ice-Shelf Melt Rates Produced by Z-coordinate Ocean Models

Antoine-Alexis Nasser<fnref n="1" />, Nicolas C. Jourdain<fnref n="1" />, Pierre Mathiot<fnref n="1" />, and Gurvan Madec<fnref n="2" />,<fnref n="3" />

Antarctic ice-shelf basal melting is a major source of uncertainty in sea level rise projections. A persistent challenge in simulating the ice-shelf-ocean interactions in z-coordinate ocean models is the introduction of artificial steps, leading to the generation of noise that impacts both melting and ocean currents. This study explores the potential of the Brinkman Volume Penalisation (BVP) method (Debreu et al. 2020, 2022) to address the recurrent issue of steps in ice-shelf-ocean models. While penalisation methods are typically applied to land topography, here, the method is generalised to ice-shelf interactions with oceans. This approach introduces porous cells that are half-ice, half-ocean, combined with a permeability parameter (friction within porous cells) to model the blocking effect of the ice draft. A unique aspect of this method is its ability to spread the penalisation region, thereby reducing model sensitivity to numerical level changes. We assess the potential benefits of the BVP approach within the idealised ice-shelf configuration ISOMIP+ as presented by Asay-Davis et al. (2016). First, a new calculation of the horizontal pressure gradient is formulated using the BVP approach, which eliminates residual biases in ocean currents down to zero machine precision. Second, the spreading of the penalised interface significantly reduces noise in the melt rates, enabling a smooth response of the ocean beneath the ice-shelf without the need for further mesh refinement. Other simulations are used to investigate the sensitivity of basal melting and freezing in the penalised configuration to changes in numerical parameters (e.g. spatial resolution). These results pave the way for a better numerical treatment of ice-shelves in earth system models.

<footnote n="1">¹Univ. Grenoble Alpes/CNRS/IRD/G-INP/INRAE, Institut des Geosciences de l'Environnement, Grenoble, France</footnote>

<footnote n="2">²Univ. Grenoble Alpes, Inria, CNRS, Grenoble INP, LJK, Grenoble, France</footnote>

<footnote n="3">³LOCEAN Laboratory, Sorbonne University-CNRS-IRD-MHNM, Paris F-75005, France</footnote>