Numerical delicacies associated with the use of isoneutral mixing operators in ocean models

Florian Lemarié 1, * Laurent Debreu 1 Patrick Marchesiello 2
* Auteur correspondant
1 MOISE - Modelling, Observations, Identification for Environmental Sciences
Inria Grenoble - Rhône-Alpes, LJK - Laboratoire Jean Kuntzmann, INPG - Institut National Polytechnique de Grenoble
Abstract : Ocean models usually rely on a tracer mixing operator which diffuses along isoneutral directions. This requirement is imposed by the highly adiabatic nature of the oceanic interior, and a numerical simulation needs to respect these small levels of dianeutral mixing to maintain physically realistic results. This is a key issue nowadays in oceanic numerical models (e.g., Hansen et al., 2011 "Earth's energy imbalance and implications"). For non-isopycnic models, rotated mixing operators must therefore be used. In continuation of the pioneering work of Griffies et al. (1998) and Beckers et al. (2000) in this field, we, first, exhaustively present the performance of various space-time discretizations in terms of stability, accuracy, tracer variance dissipation and min-max violations for the harmonic and biharmonic rotated operators. From this study, we can anticipate the flaws of the different schemes in practical situations. Because global climate models are now targeting increasingly higher horizontal resolution, the question of the viability of an isoneutral biharmonic operator is not only relevant for the regional modeling community but also for the ocean climate community. A new way of handling the temporal discretization of this type of operator is thus introduced. This scheme requires only the resolution of a simple one-dimensional tridiagonal system in the vertical direction to provide the same stability limit of the non-rotated operator. The results are illustrated by idealized numerical experiments of the diffusion of a passive tracer along isoneutral directions as well as fully realistic eddy-resolving and eddy-permitting configurations. Those numerical results show that rotated operators must be used with care and can sometimes lead to an accumulation of dispersive errors which can be seen in T/S diagrams. Furthermore, we will offer a discussion on some possible alternatives to the use of rotated operators.
Type de document :
Communication dans un congrès
European Geosciences Union General Assembly 2013, Apr 2013, Vienna, Austria. 2013
Liste complète des métadonnées

https://hal.inria.fr/hal-00826006
Contributeur : Florian Lemarié <>
Soumis le : samedi 25 mai 2013 - 18:08:49
Dernière modification le : mercredi 11 avril 2018 - 01:58:56

Identifiants

  • HAL Id : hal-00826006, version 1

Citation

Florian Lemarié, Laurent Debreu, Patrick Marchesiello. Numerical delicacies associated with the use of isoneutral mixing operators in ocean models. European Geosciences Union General Assembly 2013, Apr 2013, Vienna, Austria. 2013. 〈hal-00826006〉

Partager

Métriques

Consultations de la notice

368