Next-order balanced model captures submesoscale physics and statistics

Abstract

Using nonlinear simulations in two settings, we demonstrate that QG$^\mathrm{+1}$, a potential-vorticity (PV) based next-order-in-Rossby balanced model, captures several aspects of ocean submesoscale physics. In forced-dissipative 3D simulations under baroclinically unstable Eady-type background states, the statistical equilibrium turbulence exhibits long cyclonic tails and a plethora of rapidly-intensifying ageostrophic fronts. Despite that the model requires setting an explicit, small value for the fixed scaling Rossby number, the emergent flows are nevertheless characterized by $O(f)$ vorticity and convergence, as observed in upper-ocean submesoscale flows. Simulations of QG$^\mathrm{+1}$ under the classic strain-induced frontogenesis set-up show realistic frontal asymmetry and a finite time blow-up, quantitatively comparable to simulations of the semigeostrophic equations. The inversions in the QG$^\mathrm{+1}$ model are straightforward elliptic problems, allowing for the reconstruction of all flow fields from the PV and surface buoyancy, while avoiding the semigeostrophic coordinate transformation. Together, these results suggest QG$^\mathrm{+1}$ as a useful tool for studying upper-ocean submesoscale dynamics.