A Kinematic Kinetic Energy Backscatter Parametrization: From Implementation to Global Ocean Simulations

• Published in: Journal of Advances in Modeling Earth Systems
• Main Authors: S. Juricke, S. Danilov, N. Koldunov, M. Oliver, D. V. Sein, D. Sidorenko, Q. Wang
• Format: Article in Journal/Newspaper
• Language: English
• Published: American Geophysical Union (AGU) 2020
• Subjects: kinetic energy backscatter; viscosity closures; mesoscale eddies; ocean modeling; eddy‐permitting resolution; overdissipation; Physical geography; GB3-5030; Oceanography; GC1-1581; Sea ice
• Online Access: https://doi.org/10.1029/2020MS002175; https://doaj.org/article/cf765c74416b485482fcceeb2f3717e0

Abstract Ocean models at eddy‐permitting resolution are generally overdissipative, damping the intensity of the mesoscale eddy field. To reduce overdissipation, we propose a simplified, kinematic energy backscatter parametrization built into the viscosity operator in conjunction with a new flow‐dependent coefficient of viscosity based on nearest neighbor velocity differences. The new scheme mitigates excessive dissipation of energy and improves global ocean simulations at eddy‐permitting resolution. We find that kinematic backscatter substantially raises simulated eddy kinetic energy, similar to an alternative, previously proposed dynamic backscatter parametrization. While dynamic backscatter is scale‐aware and energetically more consistent, its implementation is more complex. Furthermore, it turns out to be computationally more expensive, as it applies, among other things, an additional prognostic subgrid energy equation. The kinematic backscatter proposed here, by contrast, comes at no additional computational cost, following the principle of simplicity. Our primary focus is the discretization on triangular unstructured meshes with cell placement of velocities (an analog of B‐grids), as employed by the Finite‐volumE Sea ice‐Ocean Model (FESOM2). The kinematic backscatter scheme with the new viscosity coefficient is implemented in FESOM2 and tested in the simplified geometry of a zonally reentrant channel as well as in a global ocean simulation on a 1/4° mesh. This first version of the new kinematic backscatter needs to be tuned to the specific resolution regime of the simulation. However, the tuning relies on a single parameter, emphasizing the overall practicality of the approach.