Ocean Kinetic Energy Backscatter Parametrization on Unstructured Grids : Impact on Global Eddy-Permitting Simulations

• Published in: Journal of Advances in Modeling Earth Systems
• Main Authors: Juricke, Stephan, Danilov, Sergey, Koldunov, Nikolay, Oliver, Marcel, Sidorenko, Dmitry
• Format: Article in Journal/Newspaper
• Language: English
• Published: Wiley-Blackwell 2020
• Subjects: Antarctic; The Antarctic; Antarc*; Sea ice
• Online Access: https://edoc.ku.de/id/eprint/30012/; https://doi.org/10.1029/2019MS001855

In this study we demonstrate the potential of a kinetic energy backscatter scheme for use in global ocean simulations. Ocean models commonly employ (bi)harmonic eddy viscosities causing excessive dissipation of kinetic energy in eddy-permitting simulations. Overdissipation not only affects the smallest resolved scales but also the generation of eddies through baroclinic instabilities, impacting the entire wave number spectrum. The backscatter scheme returns part of this overdissipated energy back into the resolved flow. We employ backscatter in the FESOM2 multiresolution ocean model with a quasi-uniform 1/4° urn:x-wiley:jame:media:jame21051:jame21051-math-0001 mesh. In multidecadal ocean simulations, backscatter increases eddy activity by a factor 2 or more, moving the simulation closer to observational estimates of sea surface height variability. Moreover, mean sea surface height, temperature, and salinity biases are reduced. This amounts to a globally averaged bias reduction of around 10% for each field, which is even larger in the Antarctic Circumpolar Current. However, in some regions such as the coastal Kuroshio, backscatter leads to a slight overenergizing of the flow and, in the Antarctic, to an unrealistic reduction of sea ice. Some of the bias increases can be reduced by a retuning of the model, and we suggest related adjustments to the backscatter scheme. The backscatter simulation is about 2.5 times as expensive as a simulation without backscatter. Most of the increased cost is due to a halving of the time step to accommodate higher simulated velocities.