A framework for the evaluation of turbulence closures used in mesoscale ocean large-eddy simulations

TL;DR

This paper introduces a methodology to evaluate and compare turbulence closures in mesoscale ocean models by analyzing their error landscapes in spectral transfer predictions, aiding in selecting optimal closures.

Contribution

The study develops an objective framework for assessing turbulence closures based on spectral transfer errors, and applies it to six different closures in 2D Navier-Stokes simulations.

Findings

Hyper-viscous closure best reproduces enstrophy cascade.

Viscous and Leith closures perform nearly as well at small scales.

Smagorinsky closure misrepresents enstrophy dissipation scales.

Abstract

We present a methodology to determine the best turbulence closure for an eddy-permitting ocean model through measurement of the error-landscape of the closure's subgrid spectral transfers and flux. We apply this method to 6 different closures for forced-dissipative simulations of the barotropic vorticity equation on a f-plane (2D Navier-Stokes equation). Using a high-resolution benchmark, we compare each closure's model of energy and enstrophy transfer to the actual transfer observed in the benchmark run. The error-landscape norms enable us to both make objective comparisons between the closures and to optimize each closure's free parameter for a fair comparison. The hyper-viscous closure most closely reproduces the enstrophy cascade, especially at larger scales due to the concentration of its dissipative effects to the very smallest scales. The viscous and Leith closures perform nearly as well, especially at smaller scales where all three models were dissipative. The Smagorinsky closure dissipates enstrophy at the wrong scales. The anticipated potential vorticity closure was the only model to reproduce the upscale transfer of kinetic energy from the unresolved scales, but would require high-order Laplacian corrections in order to concentrate dissipation at the smallest scales. The Lagrangian-averaged alpha-model closure did not perform successfully for forced 2D isotropic Navier-Stokes: small-scale filamentation is only slightly reduced by the model while small-scale roll-up is prevented. Together, this reduces the effects of diffusion.