Quasi-Geostrophic Equation Overview

• Quasi-Geostrophic Equation is a PDE framework modeling large-scale, nearly-horizontal motions in geophysical fluids, emphasizing geostrophic balance and potential vorticity dynamics.
• It is derived from the rotating shallow-water and Boussinesq equations under asymptotic expansions, revealing a Hamiltonian structure with conserved invariants like energy and enstrophy.
• Numerical methods such as Arakawa’s scheme and spectral truncation preserve its Lie–Poisson structure, enabling accurate simulation of jet formation, cascades, and turbulent processes.

The quasi-geostrophic (QG) equation is a fundamental PDE framework modeling the large-scale, slow, nearly horizontal motions in geophysical fluids (atmosphere, oceans) in regimes of small Rossby and Froude numbers. It emerges as the formal, and often rigorous, limit of the rotating, stratified Boussinesq or shallow-water systems under asymptotic expansion around geostrophic and hydrostatic balance. Mathematically, the QG equation is an active scalar equation for the advection of potential vorticity (PV), exhibiting rich Hamiltonian and Lie–Poisson geometric structures, infinite hierarchies of invariants (Casimirs), and universal cascade phenomena. Both deterministic and stochastic QG systems, on planar bounded domains and the sphere, underpin rigorous theories of quasi-Lagrangian turbulence, baroclinic instability, and jet formation.

1. Formulation and Geometric Structure

The canonical barotropic QG model on a planar β-plane domain $M\subset\mathbb R^2$ governs the evolution of the potential vorticity $q(x,t)$,

$$\partial_t q + \mathbf{u}\cdot\nabla q = \text{forcing} - \text{dissipation},$$

where $\mathbf{u} = \nabla^\perp \psi$ is non-divergent ($\nabla\cdot\mathbf{u}=0$), and the streamfunction–PV relation is

$$q = \Delta\psi + \beta y - \frac{f_0}{H}(\psi - b),$$

with $\beta$ the meridional gradient of the Coriolis parameter, $f_0$ reference Coriolis, $H$ mean depth, and $b$ topography. The two-dimensional Poisson bracket,

$$\{\psi,q\} = \partial_x\psi\,\partial_y q - \partial_y\psi\,\partial_x q,$$

encodes the skew-symmetric, area-preserving advection. On the sphere, this generalizes to

$\partial_t q + \{\psi,q\} = 0,$

with $\{\cdot,\cdot\}$ the canonical Lie–Poisson bracket adapted to the spherical geometry and the appropriate potential vorticity definition incorporating variable Coriolis parameter $f(\phi)$ and curvature terms (Franken et al., 2023, Luesink et al., 2024).

Conserved Quantities and Hamiltonian Structure

For inviscid, unforced QG, there are two primary quadratic invariants:

• Kinetic energy: $E = \tfrac12\int |\nabla\psi|^2 \, dA$,
• Enstrophy: $Z = \tfrac12\int q^2 \, dA$, as well as infinitely many Casimir invariants: $\mathcal{C}_\xi(q) = \int_M \xi(q(x)) dA$ for any smooth $\xi$ (Franken et al., 2023, Luesink et al., 2024). The system can be written in Lie–Poisson Hamiltonian form: $$\partial_t q = \{q,\mathcal{H}\}_{LP}, \quad \mathcal{H}[q] = -\tfrac12\int (q-2\mu)\psi\,dA.$$ In particular, the QG equation is a paradigmatic Lie–Poisson system on the dual of an infinite-dimensional Lie algebra of streamfunctions, with the commutator structure of the bracket justified by the geometry of area-preserving diffeomorphisms (Modin et al., 14 Apr 2025).

2. Derivation from Geophysical Fluid Dynamics

Starting from the rotating shallow-water or 3D Boussinesq equations, QG arises under asymptotic expansions with small Rossby number ($\mathrm{Ro}\ll1$), Froude number, and aspect ratio. The leading balance enforces geostrophy ($f_0\mathbf{u}^\perp+\nabla p=0$) and hydrostasy. The amplitude expansion yields a closed advection equation for PV that includes planetary curvature and topography, connecting the QG regime closely to vertical columnar and stratified layering in geophysical flows (Luesink et al., 2024, Franken et al., 2023).

The spherical QG equation is derived with a potential vorticity

$$q = \Delta\psi + 2\Omega\sin\phi - \gamma\sin^2\phi\,\psi,$$

where the Lamb parameter $\gamma=4\Omega^2R^2/gH$ quantifies the influence of stratification and geometry (Franken et al., 2023, Modin et al., 14 Apr 2025). The advection law $\partial_t q + \{\psi,q\}=0$ is structurally identical to the planar case, but with spherical differential geometry.

The rigorous link between the full stratified–rotating Boussinesq system and QG is subtle: convergence to the QG model holds provided the ratio of rotation to stratification avoids resonance ($\lambda\neq1$), but breaks down at exact resonance or if $\lambda\to1$ slowly as $\mathrm{Ro},\mathrm{Fr}\to0$ (Jo et al., 2022).

3. Numerical Discretization and Structure-Preserving Integrators

Structure-preserving numerical methods for QG exploit its Lie–Poisson and Casimir structure to guarantee long-time numerical fidelity. Major discretization strategies include:

• Arakawa's Energy–Enstrophy Scheme: This finite-difference method on a structured grid exactly conserves discrete energy and enstrophy, preserving phase-space geometric structure (Zanna, 2012).
• Spectral Zeitlin Truncation: The infinite-dimensional Lie–Poisson algebra is projected onto a sequence of finite-dimensional matrix algebras ($\mathfrak u(N)$), with the PV projected onto finite bands of spherical harmonics. This replaces the Poisson bracket with matrix commutators, exactly preserving discrete Casimirs and the Hamiltonian up to high accuracy (Franken et al., 2023).
• Finite Volume/Primal–Dual Meshes: Conservative finite-volume schemes constructed on unstructured, coastal-conforming primal–dual meshes utilize cell-to-edge and cell-to-vertex operators to approximate differential and advection terms, with advection schemes that conserve PV and potential enstrophy up to time discretization errors (Chen et al., 2017).
• Isospectral Time Integration: Midpoint and Runge–Kutta-type integrators built for coadjoint orbits precisely maintain Casimirs and nearly conserve the Hamiltonian, ensuring stable numerical statistics even over tens of thousands of time units (Franken et al., 2023, Luesink et al., 2024).

In practice, these discretizations reliably reproduce fundamental QG dynamics such as jet formation, inverse and forward cascades, and statistical equilibria, with numerical errors in invariants tightly controlling accuracy.

4. Geometric and Analytical Properties

Geometric research on the QG equation has established that the flow of QG (on the sphere, for instance) is the geodesic equation corresponding to a right-invariant weak Riemannian metric on an infinite-dimensional group (the quantomorphism group or contactomorphisms, or for 2D flows, the Lie algebra of area-preserving vector fields): $\llangle u,v\rrangle = \int_M f\,(\rho^2-\Delta)\,g\,dA,$ for suitable streamfunctions $f,g$ (Modin et al., 14 Apr 2025). The parameter $\gamma$ enters both the underlying metric and the curvature, controlling the stability and appearance of conjugate points along geodesics (trade winds).

For the surface QG (SQG) and generalized SQG families, the equation becomes a geodesic flow on the group of volume-preserving diffeomorphisms with a fractional Sobolev $\dot H^{-(1-\alpha/2)}$-metric. The Fredholm property of the exponential map and distribution of conjugate points varies sharply between the SQG and 2D Euler endpoints, implying geometric phase-transition phenomena (Bauer et al., 2023, Washabaugh, 2015).

5. Analytical Results: Well-posedness, Stochastic Forcing, and Stationary Problems

Well-posedness theory for the QG equation distinguishes between:

• Barotropic QG: For the inviscid scenario, global weak solutions exist under suitable boundary conditions (free-surface, no-flux, or no-slip), with uniform a priori $L^\infty$ bounds on PV and classical uniqueness for $L^\infty$ initial data (Chen, 2017). For forced-dissipative systems, strong solutions and global attractors exist in the subcritical dissipative regime $\alpha>1/2$ (Dlotko et al., 2014).
• SQG and gSQG: The active scalar equation for surface potential temperature $$\theta_t + u\cdot\nabla\theta + \kappa(-\Delta)^\alpha\theta=0$$, $u = \nabla^\perp(-\Delta)^{-\frac\alpha2}\theta$ splits into subcritical ($\alpha>1/2$), critical ($\alpha=1/2$), and supercritical ($0\leq\alpha<1/2$) regimes, with global regularity in the former, delicate balance and only small data global existence in the latter (Chen, 2021).
• Stochastic QG: Both martingale and strong (pathwise-unique) solutions are available for the stochastic QG equation under multiplicative or additive noise. For subcritical $\alpha>1/2$, existence and uniqueness of strong Markovian solutions (invariant measures, exponential mixing) is established; in the critical/supercritical regime, existence and Markov selection hold—uniqueness remains elusive (Röckner et al., 2011, Inahama et al., 2018).
• Stationary QG: The stationary equation $-\Delta\theta + u\cdot\nabla\theta = f$ admits sharp well-posedness in the critical homogeneous Besov spaces $\dot B^{2/p-1}_{p,q}(\mathbb{R}^2)$ for $(p,q)\in([1,4)\times[1,\infty])\cup\{(4,2)\}$, but is ill-posed outside this range; the improvement over the stationary Navier–Stokes case is due to the enhanced bilinear structure of QG (Fujii et al., 14 Mar 2025).

6. Nonlinear Dynamics: Jets, Cascades, and Fronts

High-resolution numerical and theoretical investigations reveal a diversity of QG flow phenomena:

• Zonal Jets: Numerical simulations of spherical QG show robust formation of east–west jets; for large Lamb parameter $\gamma$, jets intensify near the equator and decay toward the poles, with the jet amplitude’s latitudinal attenuation consistent with atmospheric observations (Franken et al., 2023, Luesink et al., 2024).
• Inverse/Forward Cascades: Kinetic energy spectra exhibit a $-5/3$ slope at large scales (inverse energy cascade) and $-3$ at small scales (direct enstrophy cascade). The spectrum’s anisotropy (the Rhines barrier, preferential zonalization) emerges in the QG regime (Franken et al., 2023).
• Sharp Fronts and Filamentation: The singular limit of QG with discontinuous initial data leads to contour dynamics and front models, with nonlocal, singular equations governing the motion of sharp PV fronts. Cubic weakly nonlinear approximations for gSQG fronts yield finite-time singularity formation and modulational instability (Hunter et al., 2017, Hunter et al., 2018).
• Travelling-wave Solutions: There exist explicit families of smooth, localized, steady and travelling-wave solutions (including vortex–antivortex pairs) for the inviscid SQG and generalized equations, obtained via variational methods for the associated nonlocal elliptic problems (Gravejat et al., 2017).

7. Data Assimilation, Statistical Mechanics, and Numerical Monitoring

Precise conservation of energy, enstrophy, and phase-space volume undergirds the physical fidelity of QG simulations and underlies ergodic statistical parameter estimation. Volume-preserving integrators, combined with energy–enstrophy preserving Arakawa discretizations, yield robust ensemble averages and microcanonical statistical predictions provided that discretization errors in invariants remain small over long runs (Zanna, 2012, Jolly et al., 2016). Data assimilation with coarse observations is rigorously justified for subcritical dissipation; exponential synchronization to the true state is achieved provided nudging parameters satisfy explicit balance–resolution constraints (Jolly et al., 2016).

---

References:

• Zeitlin truncation and Lie–Poisson methods: (Franken et al., 2023)
• Geometric derivation and structure-preserving simulation: (Luesink et al., 2024)
• Geodesic/Euler–Arnold structure on the quantomorphism group: (Modin et al., 14 Apr 2025)
• Energy–enstrophy conserving numerical schemes: (Zanna, 2012)
• Finite-volume coastal mesh numerics: (Chen et al., 2017)
• Stochastic QG theory: (Röckner et al., 2011, Inahama et al., 2018)
• Stationary QG well-/ill-posedness: (Fujii et al., 14 Mar 2025)
• Spherical barotropic QG and free surface: (Chen, 2017)
• Geometric analysis, Fredholm property, and criticality: (Bauer et al., 2023, Washabaugh, 2015)
• Travelling waves/variational solutions in SQG: (Gravejat et al., 2017)
• Sharp front dynamics and weakly nonlinear approximations: (Hunter et al., 2017, Hunter et al., 2018)
• Non-convergence regimes for Boussinesq to QG: (Jo et al., 2022)