pith. sign in

arxiv: 2604.25600 · v1 · submitted 2026-04-28 · ⚛️ physics.flu-dyn · cs.NA· math-ph· math.MP· math.NA

Minimum-enstrophy solutions in topographic quasi-geostrophic flow on the rotating sphere

Sagy Ephrati , Erik Jansson This is my paper

Pith reviewed 2026-05-07 15:01 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn cs.NAmath-phmath.MPmath.NA
keywords minimum enstrophyquasi-geostrophic flowrotating spheretopographynonlinear stabilityplanetary atmosphereszonal flow
0
0 comments X

The pith

Minimum-enstrophy states exist, are nonlinearly stable, and display latitude-dependent patterns in rotating spherical quasi-geostrophic flow with topography.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends the minimum-enstrophy principle from planar two-dimensional turbulence to the quasi-geostrophic equations on a rotating sphere that include bottom topography and the full nonlinear Coriolis term. It proves existence of the minimizing states and establishes their nonlinear stability at fixed energy. Analytical descriptions are given for limiting regimes set by rotation rate, topography amplitude, and energy level. Numerical solutions computed via a structure-preserving discretization reveal flows that trap around topographic features near the poles while forming zonal bands near the equator, with parameter values chosen to represent Jupiter's atmosphere. Time integration of perturbed states confirms that the solutions remain stable under the dynamics.

Core claim

Minimum-enstrophy solutions exist and are nonlinearly stable for the topographic quasi-geostrophic system on the rotating sphere; they admit closed-form asymptotic expressions for selected ranges of rotation, topography scale, and energy, and they produce latitude-dependent flows with polar topographic trapping and equatorial zonal structure.

What carries the argument

The variational minimization of potential enstrophy subject to fixed kinetic energy, performed in the presence of a latitude-dependent Coriolis parameter and bottom topography.

If this is right

  • The minimum-enstrophy states furnish explicit equilibrium profiles whose latitude dependence can be read off from the rotation rate and topography scale.
  • For Jupiter-like parameters the equilibria exhibit topographic trapping near the poles together with equatorial zonal flow.
  • Nonlinear stability guarantees that small perturbations remain bounded and do not grow under the ideal dynamics.
  • The same variational construction supplies candidate steady states for any prescribed energy and topography once the rotation rate is fixed.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the relaxation hypothesis holds, these states could serve as diagnostic templates for interpreting zonal-jet observations on gas-giant planets.
  • Adding weak dissipation or stochastic forcing would provide a direct test of whether the minimum-enstrophy attractor survives outside the ideal Hamiltonian setting.
  • The spherical geometry and latitude-dependent Coriolis term introduce qualitative differences from planar theory that can be probed by comparing solutions across different planetary rotation rates.

Load-bearing premise

The flow governed by the quasi-geostrophic equations on the sphere with topography will relax toward the minimum-enstrophy configuration at fixed energy.

What would settle it

A long-time integration of the structure-preserving discretization starting from a small perturbation of a computed minimum-enstrophy state that settles to a visibly different, higher-enstrophy configuration at the same energy.

Figures

Figures reproduced from arXiv: 2604.25600 by Erik Jansson, Sagy Ephrati.

Figure 1
Figure 1. Figure 1: Graphical summary of the paper. A necessary condition for enstrophy minimisation yields an Euler– view at source ↗
Figure 2
Figure 2. Figure 2: Topography used throughout the reported test cases. view at source ↗
Figure 3
Figure 3. Figure 3: Qualitative comparison of minimum-enstrophy solutions at different Rossby numbers. The colored back view at source ↗
Figure 4
Figure 4. Figure 4: Qualitative comparison of minimum-enstrophy solutions for different values of Lamb’s parameter view at source ↗
Figure 5
Figure 5. Figure 5: Qualitative comparison of minimum-enstrophy solutions for different energy values, at fixed Rossby view at source ↗
Figure 6
Figure 6. Figure 6: Minimum-enstrophy solution using the parameters of the Jovian atmosphere (see Table 1). In the top view at source ↗
Figure 7
Figure 7. Figure 7: Evolution of the relative deviation from the equilibrium view at source ↗
Figure 9
Figure 9. Figure 9: Log-linear plot of minimum-enstrophy Lagrange multiplier view at source ↗
read the original abstract

The minimum-enstrophy theory of Bretherton and Haidvogel postulates that two-dimensional turbulent systems evolve to a state that minimises enstrophy at a fixed energy level. We extend this to the rotating spherical quasi-geostrophic setting, accounting for bottom topography and the fully nonlinear Coriolis effect, resulting in latitude-dependent effects not present in planar approximations. We prove existence and nonlinear stability of minimum-enstrophy solutions and describe analytically asymptotic regimes for certain rates of rotation, topography scales, and energy values. We compute the minimum-enstrophy solutions by a structure-preserving method for the quasi-geostrophic equations on the sphere. We apply the method to a range of parameter values, including those describing Jupiter's atmosphere. The results reveal a distinct latitude dependence of the flow, with a tendency for topographical trapping near the poles and zonal flow near the equator, depending on the chosen parameters. The predicted nonlinear stability is confirmed numerically by integrating perturbed solutions using a structure-preserving time discretisation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper extends the Bretherton-Haidvogel minimum-enstrophy postulate to the topographic quasi-geostrophic equations on the rotating sphere with a fully nonlinear Coriolis term. It proves existence and nonlinear stability of the minimum-enstrophy states at fixed energy, derives analytical asymptotic regimes for selected rotation rates, topography scales, and energy levels, computes the states via a structure-preserving spatial discretization, applies the method to Jupiter-relevant parameters (showing polar topographic trapping and equatorial zonal flows), and numerically confirms nonlinear stability by time-integrating small perturbations of the computed minimizers with a structure-preserving scheme.

Significance. If the central claims hold, the work supplies a rigorous variational and computational framework for equilibrium states in spherical geophysical flows that incorporates curvature and latitude dependence absent from planar models. The combination of existence proofs, parameter-specific asymptotics, and structure-preserving numerics strengthens the mathematical foundation for large-scale turbulence closures in rotating fluids and offers testable predictions for planetary atmospheres. The structure-preserving methods and analytical limits are clear strengths that could be adopted more broadly in GFD modeling.

major comments (2)
  1. [Abstract] Abstract and introduction: the physical relevance of the constructed minimum-enstrophy states rests on the assumption that the spherical topographic QG dynamics relax toward them from generic initial data, yet the manuscript only establishes existence, nonlinear stability, and stability under small perturbations; no long-time integrations from turbulent, random, or high-enstrophy initial conditions are reported to test whether enstrophy actually decreases to the variational minimum (as required to extend the Bretherton-Haidvogel selection mechanism).
  2. [Numerical experiments] Numerical experiments section: the structure-preserving time discretization is used solely to verify that perturbed minimizers remain close; without accompanying relaxation tests or comparison against direct enstrophy decay from non-equilibrium states, the numerical evidence does not address the core dynamical selection postulate being extended to the sphere.
minor comments (2)
  1. [Formulation] Notation for the nonlinear Coriolis term and the precise definition of the energy and enstrophy functionals on the sphere should be cross-referenced explicitly between the variational formulation and the numerical scheme to avoid ambiguity for readers.
  2. [Asymptotics] The asymptotic regimes are described analytically but would benefit from a short table summarizing the leading-order balances for each parameter regime (rotation rate, topography scale, energy) to improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. We address the major comments point by point below, clarifying the scope of our contributions while acknowledging the limitations noted.

read point-by-point responses

  1. Referee: [Abstract] Abstract and introduction: the physical relevance of the constructed minimum-enstrophy states rests on the assumption that the spherical topographic QG dynamics relax toward them from generic initial data, yet the manuscript only establishes existence, nonlinear stability, and stability under small perturbations; no long-time integrations from turbulent, random, or high-enstrophy initial conditions are reported to test whether enstrophy actually decreases to the variational minimum (as required to extend the Bretherton-Haidvogel selection mechanism).

    Authors: We agree that the manuscript does not include long-time integrations from generic turbulent or high-enstrophy initial conditions to directly demonstrate enstrophy decay to the variational minimum. The work extends the Bretherton-Haidvogel postulate by establishing rigorous existence and nonlinear stability of the minimum-enstrophy states for the spherical topographic QG system, including latitude-dependent effects. The nonlinear stability result implies local attractivity for small perturbations, providing mathematical support for the selection mechanism. Full relaxation tests from arbitrary initial data are computationally demanding and lie outside the present scope, which centers on the variational characterization, analytical asymptotics, and structure-preserving computation of the equilibria. We will add a clarifying paragraph in the introduction and conclusions to explicitly state the scope and note that dynamical relaxation remains a conjecture supported by the stability analysis, analogous to the original planar theory. revision: partial

  2. Referee: [Numerical experiments] Numerical experiments section: the structure-preserving time discretization is used solely to verify that perturbed minimizers remain close; without accompanying relaxation tests or comparison against direct enstrophy decay from non-equilibrium states, the numerical evidence does not address the core dynamical selection postulate being extended to the sphere.

    Authors: The numerical section employs the structure-preserving discretization to compute the minimum-enstrophy states and to confirm their nonlinear stability via small-perturbation time integrations, as required to validate the analytical results. We acknowledge that this does not include relaxation experiments tracking enstrophy decay from non-equilibrium states. Such tests would provide additional numerical support for the dynamical selection but require substantial extra resources for long-time spherical simulations across parameter regimes. The current evidence substantiates the stability of the predicted states, which is a necessary ingredient for the extended postulate. We will revise the numerical experiments discussion to include an explicit statement of this limitation and the rationale for focusing on stability verification. revision: partial

Circularity Check

0 steps flagged

No significant circularity; variational proofs and direct numerical verification are independent of the imported postulate

full rationale

The paper imports the minimum-enstrophy selection principle from Bretherton-Haidvogel (prior literature, no author overlap) as a modeling assumption and then proves existence plus nonlinear stability of the resulting states via variational arguments on the sphere. Numerical work consists of structure-preserving time integration applied only to small perturbations of the computed minimizers to confirm the stability prediction; no parameters are fitted to data and then relabeled as predictions, no self-citation chain supports a uniqueness theorem, and no ansatz is smuggled. The derivation chain therefore remains self-contained against external benchmarks and does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claim rests on the imported minimum-enstrophy variational principle together with the quasi-geostrophic approximation on the sphere; numerical applications introduce free parameters for rotation rate, topography amplitude, and energy level that are varied but not derived from first principles.

free parameters (3)
  • rotation rate
    Varied across a range of values including those for Jupiter; controls the strength of the Coriolis term.
  • energy level
    Held fixed while enstrophy is minimized; appears as a constraint in the variational problem.
  • topography scale
    Amplitude and spatial scale of bottom topography treated as input parameters.
axioms (2)
  • domain assumption The system evolves toward a state of minimum enstrophy at fixed energy
    Postulated by Bretherton and Haidvogel for planar flows and extended without additional derivation to the spherical case with nonlinear Coriolis.
  • domain assumption Quasi-geostrophic balance holds on the sphere
    Standard approximation invoked to obtain the governing equations.

pith-pipeline@v0.9.0 · 5486 in / 1441 out tokens · 141339 ms · 2026-05-07T15:01:02.329287+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

44 extracted references · 44 canonical work pages

  1. [1]

    V. I. Arnold. An a priori estimate in the theory of hydrodynamic stability.Izv. Vyssh. Uchebn. Zaved. Mat.[Sov. Math. J.], 5:3, 1966

    work page 1966

  2. [2]

    V. I. Arnold and B. A. Khesin.Topological methods in hydrodynamics. Springer, 1998

    work page 1998

  3. [3]

    Biskamp and E

    D. Biskamp and E. Schwarz. On two-dimensional magnetohydrodynamic turbulence.Physics of Plasmas, 8(7):3282–3292, 2001

    work page 2001

  4. [4]

    Boffetta and R

    G. Boffetta and R. E. Ecke. Two-dimensional turbulence.Annual review of fluid mechanics, 44(1):427–451, 2012

    work page 2012

  5. [5]

    Bordemann, J

    M. Bordemann, J. Hoppe, P. Schaller, and M. Schlichenmaier.gl(∞) and geometric quantization.Communi- cations in Mathematical Physics, 138(2):209–244, 1991

    work page 1991

  6. [6]

    Bordemann, E

    M. Bordemann, E. Meinrenken, and M. Schlichenmaier. Toeplitz quantization of K¨ ahler manifolds and gl(N), N→ ∞limits.Communications in Mathematical Physics, 165(2):281–296, 1994

    work page 1994

  7. [7]

    Brands, P

    H. Brands, P. Chavanis, R. Pasmanter, and J. Sommeria. Maximum entropy versus minimum enstrophy vortices.Physics of Fluids, 11(11):3465–3477, 1999. 15

    work page 1999

  8. [8]

    F. P. Bretherton and D. B. Haidvogel. Two-dimensional turbulence above topography.Journal of Fluid Mechanics, 78(1):129–154, 1976

    work page 1976

  9. [9]

    G. F. Carnevale and J. S. Frederiksen. Nonlinear stability and statistical mechanics of flow over topography. Journal of Fluid Mechanics, 175:157–181, 1987

    work page 1987

  10. [10]

    J. G. Charney and M. E. Stern. On the stability of internal baroclinic jets in a rotating atmosphere.Journal of Atmospheric Sciences, 19(2):159–172, 1962

    work page 1962

  11. [11]

    Cifani, M

    P. Cifani, M. Viviani, E. Luesink, K. Modin, and B. J. Geurts. Casimir preserving spectrum of two-dimensional turbulence.Physical Review Fluids, 7(8):L082601, 2022

    work page 2022

  12. [12]

    Cifani, M

    P. Cifani, M. Viviani, and K. Modin. An efficient geometric method for incompressible hydrodynamics on the sphere.Journal of Computational Physics, 473:111772, 2023

    work page 2023

  13. [13]

    G. P. Cressman. Barotropic divergence and very long atmospheric waves.Monthly weather review, 86(8):293– 297, 1958

    work page 1958

  14. [14]

    R. Daley. Linear non-divergent mass-wind laws on the sphere.Tellus A: Dynamic Meteorology and Oceanog- raphy, 35(1):17–27, 1983

    work page 1983

  15. [15]

    Ephrati, E

    S. Ephrati, E. Jansson, A. Lang, and E. Luesink. An exponential map free implicit midpoint method for stochastic Lie–Poisson systems.arXiv preprint arXiv:2408.16701, 2024

    work page arXiv 2024

  16. [16]

    Ephrati, E

    S. Ephrati, E. Jansson, and K. Modin. On spectral scaling laws for averaged turbulence on the sphere.Physica D: Nonlinear Phenomena, page 134808, 2025

    work page 2025

  17. [17]

    Franken, E

    A. Franken, E. Luesink, S. Ephrati, and B. Geurts. Critical latitude in global quasi-geostrophic flow on a rotating sphere.arXiv preprint arXiv:2409.05432, 2024

    work page arXiv 2024

  18. [18]

    A. D. Franken, M. Caliaro, P. Cifani, and B. J. Geurts. Zeitlin truncation of a shallow water quasi-geostrophic model for planetary flow.Journal of Advances in Modeling Earth Systems, 16(6):e2023MS003901, 2024

    work page 2024

  19. [19]

    A. D. Franken, E. Luesink, S. R. Ephrati, and B. J. Geurts. Casimir preserving numerical method for global multi-layer quasi-geostrophic turbulence.Journal of Computational Physics, 538:114155, 2025

    work page 2025

  20. [20]

    B. Gallet. Two-dimensional turbulence above topography: condensation transition and selection of minimum enstrophy solutions.Journal of Fluid Mechanics, 988:A13, 2024

    work page 2024

  21. [21]

    D. D. Holm, J. E. Marsden, T. Ratiu, and A. Weinstein. Nonlinear stability of fluid and plasma equilibria. Physics reports, 123(1-2):1–116, 1985

    work page 1985

  22. [22]

    J. Hoppe. Diffeomorphism groups, quantization, andSU(∞).International Journal of Modern Physics A, 4(19):5235–5248, 1989

    work page 1989

  23. [23]

    W. Kahan. Accurate eigenvalues of a symmetric tridiagonal matrix. Technical Report CS41, Computer Science Department, Stanford University, Stanford, CA, July 1966

    work page 1966

  24. [24]

    H. Kuo. Finite-amplitude three-dimensional harmonic waves on the spherical earth.Journal of Atmospheric Sciences, 16(5):524–534, 1959

    work page 1959

  25. [25]

    J. H. Lacasce, A. Pal´ oczy, and M. Trodahl. Vortices over bathymetry.Journal of Fluid Mechanics, 979:A32, 2024

    work page 2024

  26. [26]

    Luesink, A

    E. Luesink, A. Franken, S. Ephrati, and B. Geurts. Geometric derivation and structure-preserving simulation of quasi-geostrophy on the sphere.arXiv preprint arXiv:2402.13707, 2024

    work page arXiv 2024

  27. [27]

    Modin and M

    K. Modin and M. Roop. Spatio-temporal Lie–Poisson discretization for incompressible magnetohydrodynamics on the sphere.IMA Journal of Numerical Analysis, page draf024, 2025

    work page 2025

  28. [28]

    Modin and M

    K. Modin and M. Viviani. A Casimir preserving scheme for long-time simulation of spherical ideal hydrody- namics.Journal of Fluid Mechanics, 884:A22, 2020. 16

    work page 2020

  29. [29]

    Modin and M

    K. Modin and M. Viviani. A brief introduction to matrix hydrodynamics.Journal of Computational Dynamics, 14(0):17–35, 2026

    work page 2026

  30. [30]

    P. J. Morrison. A paradigm for joined Hamiltonian and dissipative systems.Physica D: Nonlinear Phenomena, 18(1-3):410–419, 1986

    work page 1986

  31. [31]

    V. K. Priya, S. S. Patil, K. Seshasayanan, and R. Lakkaraju. Two-dimensional turbulence over topography of varying roughness.Journal of Fluid Mechanics, 1033:A33, 2026

    work page 2026

  32. [32]

    P. B. Rhines. Waves and turbulence on a beta-plane.Journal of Fluid Mechanics, 69(3):417–443, 1975

    work page 1975

  33. [33]

    Roop and S

    M. Roop and S. Ephrati. Thermal quasi-geostrophic model on the sphere: Derivation and structure-preserving simulation.Physics of Fluids, 37(9), 2025

    work page 2025

  34. [34]

    R. Salmon. Geostrophic turbulence.Topics in ocean physics, 30:78, 1982

    work page 1982

  35. [35]

    Salmon, G

    R. Salmon, G. Holloway, and M. C. Hendershott. The equilibrium statistical mechanics of simple quasi- geostrophic models.Journal of Fluid Mechanics, 75(4):691–703, 1976

    work page 1976

  36. [36]

    L. Z. Sans´ on, A. Gonz´ alez-Villanueva, and L. Flores. Evolution and decay of a rotating flow over random topography.Journal of fluid mechanics, 642:159–180, 2010

    work page 2010

  37. [37]

    W. H. Schubert, R. K. Taft, and L. G. Silvers. Shallow water quasi-geostrophic theory on the sphere.Journal of Advances in Modeling Earth Systems, 1(2), 2009

    work page 2009

  38. [38]

    Siegelman and W

    L. Siegelman and W. R. Young. Two-dimensional turbulence above topography: Vortices and potential vorticity homogenization.Proceedings of the National Academy of Sciences, 120(44):e2308018120, 2023

    work page 2023

  39. [39]

    J. Theiss. Equatorward energy cascade, critical latitude, and the predominance of cyclonic vortices in geostrophic turbulence.Journal of physical oceanography, 34(7):1663–1678, 2004

    work page 2004

  40. [40]

    W. T. Verkley. A balanced approximation of the one-layer shallow-water equations on a sphere.Journal of the atmospheric sciences, 66(6):1735–1748, 2009

    work page 2009

  41. [41]

    E. W. Weisstein. Wigner 3j-symbol. FromMathWorld– A Wolfram Web Resource. Accessed 2026-03-02

    work page 2026

  42. [42]

    V. Zeitlin. Finite-mode analogs of 2D ideal hydrodynamics: Coadjoint orbits and local canonical structure. Physica D: Nonlinear Phenomena, 49(3):353–362, 1991

    work page 1991

  43. [43]

    V. Zeitlin. Self-consistent finite-mode approximations for the hydrodynamics of an incompressible fluid on nonrotating and rotating spheres.Physical review letters, 93(26):264501, 2004

    work page 2004

  44. [44]

    Zeitlin.Geophysical fluid dynamics: understanding (almost) everything with rotating shallow water models

    V. Zeitlin.Geophysical fluid dynamics: understanding (almost) everything with rotating shallow water models. Oxford University Press, 2018. A Gaunt coefficients Below, we briefly elaborate on the triple product of spherical harmonics appearing in Section 3.2. We letY l,m(φ, θ) denote the complex spherical harmonic of degreeland orderm, and will omit the c...

    work page 2018