A conservative low-order model for Boussinesq baroclinic fronts

TL;DR

This paper develops a low-order, nonlinear model of Boussinesq baroclinic fronts that captures the finite-time adjustment mechanisms between turbulent eddies and flow balance, extending quasi-balanced models.

Contribution

It introduces a five-dimensional nonlinear ODE system derived from the Boussinesq equations, explicitly resolving the inertial lag of secondary circulation.

Findings

The model conserves total energy and cross-frontal density gradient magnitude.

It describes the adiabatic adjustment as a rotation of the density gradient slope.

The dynamics are non-Hamiltonian due to turbulent closure effects.

Abstract

The internal dynamics of baroclinic fronts are governed by a fundamental interplay: turbulent eddies systematically act to disrupt thermal wind balance, with baroclinic eddies flattening isopycnals and barotropic momentum fluxes intensifying the primary jet, while the ageostrophic overturning circulation acts to restore it. In quasi-balanced models, this restorative adjustment is assumed instantaneous, locking the flow onto a balanced manifold. To conceptually track this mechanism when the adjustment takes a finite time, we construct a low-order model that spans from $\mathcal{O}(1)$ Rossby numbers down to the quasi-balanced limit. Formulated from the continuous Boussinesq equations under a $Ro^2 Ri \sim 1$ scaling, which constrains the horizontal length scale to the Rossby deformation radius, the derivation yields a closed, nonlinear five-dimensional ODE system. The degrees of freedom consist of the domain-averaged along-front vertical shear, the cross-frontal overturning vorticity, the horizontal and vertical buoyancy gradients, and the total eddy energy. We identify two constants of motion that constrain the evolution of the mean flow: the total energy (kinetic energy of the along- and cross-frontal flows, mean potential energy, and eddy energy) and the magnitude of the domain-averaged cross-frontal density gradient. Notably, while the system is energetically conservative, the parameterized turbulent closure renders the dynamics strictly non-Hamiltonian. Bounded by these invariants, the adiabatic adjustment of the front physically reduces to a continuous rotation of the density gradient's slope. By explicitly resolving the inertial lag of the secondary circulation, this framework isolates the individual mechanisms governing frontal adjustment and tracks their continuous dynamic interplay.