A scale‐aware subgrid model for quasi‐geostrophic turbulence

• Published in: Journal of geophysical research. Oceans 2017-02, Vol.122 (2), p.1529-1554
• Main Authors: Bachman, Scott D., Fox‐Kemper, Baylor, Pearson, Brodie
• Format: Article
• Language: English
• Subjects: Advection; baroclinic instability; Boundary layers; Cascades; Climate; Climate models; Closures; Computational fluid dynamics; Convection; Cyclones; Diapycnal mixing; diffusion; Dynamics; eddies; Energy; Energy consumption; Enstrophy; Equator; Filters; Geophysics; Geostrophic turbulence; High resolution; Large eddy simulation; LES; Marine; Mathematical models; mesoscale; Methods; Numerical stability; Ocean models; Oceanic turbulence; Oceanography; Oceans; Parameter sensitivity; parameterization; potential enstrophy; Primitive equation models; Primitive equations; quasi‐geostrophic; Resolution; Robustness (mathematics); Scale (ratio); scale‐aware; Simulation; Spectra; subgrid model; Test procedures; Turbulence; Turbulent flow; Two dimensional models; Uncertainty; Viscosity; Vortices
• ISSN: 2169-9275; 2169-9291
• DOI: 10.1002/2016JC012265

This paper introduces two methods for dynamically prescribing eddy‐induced diffusivity, advection, and viscosity appropriate for primitive equation models with resolutions permitting the forward potential enstrophy cascade of quasi‐geostrophic dynamics, such as operational ocean models and high‐resolution climate models with O(25) km horizontal resolution and finer. Where quasi‐geostrophic dynamics fail (e.g., the equator, boundary layers, and deep convection), the method reverts to scalings based on a matched two‐dimensional enstrophy cascade. A principle advantage is that these subgrid models are scale‐aware, meaning that the model is suitable over a range of grid resolutions: from mesoscale grids that just permit baroclinic instabilities to grids below the submesoscale where ageostrophic effects dominate. Two approaches are presented here using Large Eddy Simulation (LES) techniques adapted for three‐dimensional rotating, stratified turbulence. The simpler approach has one nondimensional parameter, Λ, which has an optimal value near 1. The second approach dynamically optimizes Λ during simulation using a test filter. The new methods are tested in an idealized scenario by varying the grid resolution, and their use improves the spectra of potential enstrophy and energy in comparison to extant schemes. The new methods keep the gridscale Reynolds and Péclet numbers near 1 throughout the domain, which confers robust numerical stability and minimal spurious diapycnal mixing. Although there are no explicit parameters in the dynamic approach, there is strong sensitivity to the choice of test filter. Designing test filters for heterogeneous ocean turbulence adds cost and uncertainty, and we find the dynamic method does not noticeably improve over setting Λ = 1. Key Points A subgrid closure for mesoscale ocean LES is introduced The new closure is flow and scale‐aware and specifies the Gent and McWilliams transport coefficient The new closure compares favorably to extant closures in truncating the cascades of energy and potential enstrophy