Semi-analytical eddy-viscosity and backscattering closures for 2D geophysical turbulence

Abstract: Physics-based closures such as eddy-viscosity and backscattering models are widely used for large-eddy simulation (LES) of geophysical turbulence for applications including weather and climate prediction. However, these closures have parameters that are often chosen empirically. Here, for the first time, we semi-analytically derive the parameters of the Leith and Smagorinsky eddy-viscosity closures and the Jansen-Held backscattering closure for 2D geophysical turbulence. The semi-analytical derivation provides these parameters up to a constant that can be estimated from the turbulent kinetic energy spectrum of a few snapshots of direct numerical simulation (DNS) or other high-fidelity (eddy resolving) simulations, or even obtained from earlier analytical work based on renormalization group. The semi-analytically estimated closure parameters agree with those obtained from online (a-posteriori) learning in several setups of 2D geophysical turbulence in our earlier work. LES with closures that use these parameters can correctly reproduce the key statistics of DNS, including those of the extreme events and interscale energy and enstrophy transfers, and outperform the baselines (dynamic Leith and Smagorinsky and the latter with standard parameter).