Subgrid Parameterization of Eddy, Meanfield and Topographic Interactions in Simulations of an Idealized Antarctic Circumpolar Current.

The ocean circulation dynamics can be represented as a high‐dimensional multi‐scale nonlinear system with inhomogenous meanfields and topography. Numerical simulations of the ocean dynamics resolve the large scales of motion on a computational grid, with the unresolved subgrid interactions parameterized. These simulations are highly dependent upon the grid resolution, unless the subgrid terms are appropriately parameterized. There are five fundamental classes of subgrid interactions: eddy‐eddy; eddy‐topographic; eddy‐meanfield; meanfield‐meanfield; and meanfield‐topographic. Scale dependent parameterizations representing each of these interaction classes are presented here in oceanic flows for the first time. Subgrid parameterizations are calculated and validated in baroclinic quasi‐geostrophic simulations of idealized Antarctic Circumpolar Current flows with representative mean currents and ocean floor topography. The parameterization coefficients are derived from the coarse grained statistics of high resolution reference simulations in spectral space. Stochastic and deterministic parameterizations are developed for the eddy‐eddy interactions, and deterministic forms for the remaining classes. The kinetic energy spectra and meanfield resulting from large eddy simulations (LES) adopting these coefficients accurately replicate those of the reference simulation. The eddy‐eddy interactions are dominant, but all classes need to be parameterized for best results. For LES where baroclinic instability is explicitly resolved the stochastic variants out‐perform the deterministic ones across all scales. When baroclinic instability is not explicitly resolved, the stochastic variants out‐perform the deterministic cases at the large scales, but introduce some distortions at the smallest resolved scales. This study provides an assessment of the relevant strengths of the subgrid interaction classes in an idealized yet representative ocean. Plain Language Summary: Numerical simulation of oceanic circulations is incredibly challenging due in part to the large range of scales. There are oscillations on the size of entire ocean basins, through to fine scale centimeter sized turbulence, and everything in between. Global simulations of the ocean typically capture the large scale features on a discrete grid. To ensure these simulations are physically representative, the influence of the missing smaller scales needs to be captured accurately. Here we present a first calculation of parameterizations governing the interactions between the resolved dynamics and the missing small scale eddies, meanfield and topographic features. Reduced resolution simulations incorporating these parameterizations without tuning have statistics in very good agreement with the high resolution simulations from which they were derived. Key Points: Subgrid models of all interaction classes between eddies, meanfield and topography are calculated from reference global ocean simulationsInter‐eddy interactions are modeled both stochastically and deterministically, with deterministic models used for the other interactionsCoarse simulations using these subgrid coefficients without tuning are in very good agreement with those of the high‐resolution simulation [ABSTRACT FROM AUTHOR]