Final states of two-dimensional turbulence above large-scale topography: stationary vortex solutions and barotropic stability
Pith reviewed 2026-05-19 06:41 UTC · model grok-4.3
The pith
An empirical model superposing Gaussian vortices on linear background flows produces locally stationary solutions for final states of two-dimensional topographic turbulence.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The final states of 2D topographic turbulence consist of a background flow satisfying a linear PV-streamfunction relation together with localized vortices that, after the background is subtracted, follow a sinh-like PV-streamfunction relation. An empirical model that combines this linear background flow with Gaussian vortices centered at the topographic extrema accurately reproduces the quasi-stationary states obtained in simulations and yields locally stationary solutions to the inviscid governing equation. Linear stability analyses of these stationary vortex solutions reveal background-flow-dependent stability: cyclone-elevation and anticyclone-depression configurations are stable at low背景
What carries the argument
Empirical model formed by adding Gaussian vortex profiles centered at topographic extrema to a linear background flow that satisfies the stationary PV-streamfunction relation.
If this is right
- The model supplies explicit vortex solutions for quasi-stationary final states of 2D topographic turbulence.
- Stability of each vortex-topography configuration depends on the background energy level.
- The observed vortex-topography correlations in simulations are explained by the energy-dependent stability thresholds.
- Both the sinh-like relation and the Gaussian profiles persist under complex topography and high-energy conditions.
Where Pith is reading between the lines
- The same construction might be applied to other topographic shapes to test whether Gaussian centering remains a good approximation.
- The stability criteria could be used to forecast which vortex signs are favored over seamounts or trenches in oceanic flows at different energy levels.
- Weak viscosity could be added to the model to examine how the stationary states slowly erode or reorganize over long times.
Load-bearing premise
Localized vortices in the final states can be accurately represented by Gaussian profiles centered exactly at the topographic extrema once the linear background flow is subtracted.
What would settle it
High-resolution simulations in which the potential-vorticity profiles of the locked vortices deviate measurably from the Gaussian shape after background subtraction, or in which the vortices fail to remain centered on the extrema.
Figures
read the original abstract
The final states of freely decaying two-dimensional (2D) topographic turbulence consist of a background flow and localized vortices. While the background flow satisfies a linear potential vorticity (PV)-streamfunction relation, the vortex structures remain poorly understood. To address this gap and ensure oceanic relevance, we examine quasi-stationary final states of 2D turbulence over a sinusoidal topography featuring a bump and a dip, where two oppositely signed vortices are locked to the topographic extrema. After subtracting the background flow, the vortices exhibit a "sinh"-like PV-streamfunction relation, as observed in flat-bottom turbulence. Motivated by Gaussian vortex profiles in flat-bottom turbulence, we propose an empirical model combining the background flow with Gaussian vortices centered at the topographic extrema. This model accurately reproduces quasi-stationary states and yields locally stationary solutions to the inviscid governing equation. We further test the model under complex topography and high-energy conditions, confirming that the "sinh"-like trend and Gaussian profiles are robust features of localized vortices. Linear stability analyses of these stationary vortex solutions reveal background flow-dependent stability: cyclone/elevation and anticyclone/depression configurations are stable at low background energy, while anticyclone/elevation and cyclone/depression configurations are stable at high background energy. These findings align with vortex-topography correlations observed in simulations across energy regimes. Our results provide explicit vortex solutions for quasi-stationary final states of 2D topographic turbulence and elucidate the mechanism underlying vortex-topography correlations through stability analyses of vortices embedded in topographic background flows.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the final states of freely decaying two-dimensional turbulence over sinusoidal topography containing a bump and a dip. It reports that these states comprise a linear background flow satisfying a PV-streamfunction relation together with oppositely signed localized vortices locked to the topographic extrema. After subtracting the background, the vortices display a sinh-like PV-streamfunction relation. The authors introduce an empirical model that superposes the linear background with Gaussian vortex profiles centered at the extrema; they claim this construction reproduces the quasi-stationary states observed in simulations, furnishes locally stationary solutions of the inviscid equation, and remains robust under more complex topography and higher energies. Linear stability analysis of the resulting vortex solutions is shown to be background-energy dependent, with cyclone-over-elevation and anticyclone-over-depression configurations stable at low background energy and the opposite pairings stable at high background energy; these stability properties are said to explain the vortex-topography correlations seen across the simulations.
Significance. If the empirical construction is shown to satisfy the stationary condition and the stability results are confirmed, the work supplies explicit, reproducible vortex solutions for topographic turbulence final states and a mechanistic account of preferred vortex-topography alignments. Such solutions would be useful for initializing or parameterizing oceanic models over rough bathymetry and for interpreting satellite or drifter observations of coherent vortices.
major comments (2)
- [Empirical model section] Empirical model section: the assertion that the superposition of the linear background flow, sinusoidal topography, and Gaussian vortex streamfunctions produces locally stationary solutions requires an explicit verification that the composite potential vorticity q remains a single-valued function of the composite streamfunction ψ near the extrema. The background satisfies a linear relation by construction and the isolated vortices exhibit a sinh-like trend, but the Laplacian of the added Gaussian profile plus the topographic term h does not automatically preserve a functional relation q = F(ψ); the manuscript should report a quantitative check (e.g., scatter-plot collapse or maximum deviation from a fitted F) for the composite fields.
- [Stability analysis section] Stability analysis section: the reported background-energy dependence of stability (cyclone/elevation stable at low energy, anticyclone/elevation stable at high energy) is load-bearing for the explanation of observed correlations, yet the precise nondimensional energy thresholds separating the two regimes are not stated. The manuscript should define these thresholds in terms of the nondimensional parameters used in the linear stability eigenvalue problem and show that the transition occurs within the range of energies realized in the turbulence simulations.
minor comments (2)
- Figure captions should explicitly state the nondimensional parameters (e.g., topographic amplitude, Reynolds number, initial energy) corresponding to each panel so that the reader can map the displayed fields to the stability regimes discussed in the text.
- The phrase “locally stationary solutions” is used repeatedly; a brief parenthetical definition or reference to the precise diagnostic (e.g., residual of the stationary PV equation) would remove ambiguity.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments, which help clarify the presentation of our empirical model and stability results. We address each major comment below and will incorporate the suggested clarifications into the revised manuscript.
read point-by-point responses
-
Referee: [Empirical model section] Empirical model section: the assertion that the superposition of the linear background flow, sinusoidal topography, and Gaussian vortex streamfunctions produces locally stationary solutions requires an explicit verification that the composite potential vorticity q remains a single-valued function of the composite streamfunction ψ near the extrema. The background satisfies a linear relation by construction and the isolated vortices exhibit a sinh-like trend, but the Laplacian of the added Gaussian profile plus the topographic term h does not automatically preserve a functional relation q = F(ψ); the manuscript should report a quantitative check (e.g., scatter-plot collapse or maximum deviation from a fitted F) for the composite fields.
Authors: We agree that an explicit quantitative verification strengthens the claim that the composite fields yield locally stationary solutions. Although the background satisfies a linear relation by construction and the isolated vortices follow a sinh-like relation, the superposition with the Gaussian profile and topography does not automatically guarantee a single-valued q(ψ) everywhere. In the revised manuscript we will add a scatter plot of composite q versus ψ in the vicinity of each topographic extremum, together with the maximum deviation from a fitted functional form (linear for the background component and sinh-like for the vortex component). This will quantify the accuracy of the empirical construction as a locally stationary solution. revision: yes
-
Referee: [Stability analysis section] Stability analysis section: the reported background-energy dependence of stability (cyclone/elevation stable at low energy, anticyclone/elevation stable at high energy) is load-bearing for the explanation of observed correlations, yet the precise nondimensional energy thresholds separating the two regimes are not stated. The manuscript should define these thresholds in terms of the nondimensional parameters used in the linear stability eigenvalue problem and show that the transition occurs within the range of energies realized in the turbulence simulations.
Authors: We acknowledge that explicit thresholds improve the connection between the linear stability analysis and the simulation results. In the revised manuscript we will define the critical nondimensional background energies (expressed in terms of the parameters of the eigenvalue problem, such as the nondimensional topographic amplitude and the background flow strength) that separate the two stability regimes. We will also show that these thresholds fall within the range of background energies attained in the freely decaying turbulence simulations, thereby confirming that the observed vortex-topography correlations are consistent with the stability properties. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper constructs an empirical model by superposing a linear background flow (which satisfies a linear PV-streamfunction relation by design) with Gaussian vortex profiles motivated by flat-bottom observations and simulation data. It then verifies that this composite approximately satisfies the inviscid stationary condition locally and performs independent linear stability analysis on the resulting solutions, comparing outcomes to separate numerical simulations. No load-bearing step reduces by the paper's own equations to a fitted parameter renamed as a prediction, nor does any central claim rely on a self-citation chain or imported ansatz that is itself unverified. The derivation remains self-contained against external simulation benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- Gaussian vortex amplitude and width parameters
axioms (2)
- standard math The governing equations are the inviscid 2D Euler (or quasi-geostrophic) equations with topography entering the potential vorticity.
- domain assumption Final states of freely decaying topographic turbulence consist of a linear background flow plus localized vortices locked to topographic extrema.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
empirical model combining the background flow with Gaussian vortices centered at the topographic extrema... yields locally stationary solutions to the inviscid governing equation
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leancostAlphaLog_high_calibrated_iff echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
vortices... exhibit a 'sinh'-like PV-streamfunction relation
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year eprint label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence ...
-
[2]
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...
-
[3]
Agoua, W. , Yin, X.-Y. , Wu, T. & Bos, W. J. T. 2025 Coexistence of two equilibrium configurations in two-dimensional turbulence . Phys. Rev. Fluids 10 (3), 034604
work page 2025
-
[4]
Benzi, R. , Paladin, G. , Patarnello, S. , Santangelo, P. & Vulpiani, A. 1986 Intermittency and coherent structures in two-dimensional turbulence . J. Phys. A: Math. Gen. 19 (18), 3771
work page 1986
-
[5]
Benzi, R. , Patarnello, S. & Santangelo, P. 1987 On the Statistical Properties of Two-Dimensional Decaying Turbulence . EPL 3 (7), 811
work page 1987
-
[6]
Benzi, R. , Patarnello, S. & Santangelo, P. 1988 Self-similar coherent structures in two-dimensional decaying turbulence . J. Phys. A: Math. Gen. 21 (5), 1221
work page 1988
-
[7]
Bouchet, F. & Venaille, A. 2012 Statistical mechanics of two-dimensional and geophysical flows . Phys. Rep. 515 (5), 227--295
work page 2012
-
[8]
Bretherton, F. P. & Haidvogel, D. B. 1976 Two-dimensional turbulence above topography . J. Fluid Mech. 78 (1), 129--154
work page 1976
-
[9]
Carnevale, G. F. & Frederiksen, J. S. 1987 Nonlinear stability and statistical mechanics of flow over topography . J. Fluid Mech. 175 , 157--181
work page 1987
-
[10]
Carnevale, G. F. , Kloosterziel, R. C. & Van Heijst, G. J. F. 1991 Propagation of barotropic vortices over topography in a rotating tank . J. Fluid Mech. 233 , 119--139
work page 1991
-
[11]
Charney, J. G. & DeVore, J. G. 1979 Multiple Flow Equilibria in the Atmosphere and Blocking . J. Atmos. Sci. 36 (7), 1205--1216
work page 1979
-
[12]
Constantinou, N. , Wagner, G. , Siegelman, L. , Pearson, B. & Pal \'o czy, A. 2021 Geophysicalflows. jl: Solvers for geophysical fluid dynamics problems in periodic domains on cpus gpus . J. Open Source Softw. 6 (60)
work page 2021
-
[13]
Constantinou, N. C. 2018 A Barotropic Model of Eddy Saturation . J. Phys. Oceanogr. 48 (2), 397--411
work page 2018
-
[14]
Constantinou, N. C. & Hogg, A. M. 2019 Eddy Saturation of the Southern Ocean : A Baroclinic Versus Barotropic Perspective . Geophys. Res. Lett. 46 (21), 12202--12212
work page 2019
-
[15]
Constantinou, N. C. & Young, W. R. 2017 Beta-plane turbulence above monoscale topography . J. Fluid Mech. 827 , 415--447
work page 2017
-
[16]
DiBattista, M. T. & Majda, A. J. 2001 Equilibrium Statistical Predictions for Baroclinic Vortices : The Role of Angular Momentum . Theoret. Comput. Fluid Dynamics 14 (5), 293--322
work page 2001
-
[17]
Dritschel, D. G. , Qi, W. & Marston, J. B. 2015 On the late-time behaviour of a bounded, inviscid two-dimensional flow . J. Fluid Mech. 783 , 1--22
work page 2015
-
[18]
Eyink, G. L. 2014 Turbulence theory iii. https://www.ams.jhu.edu/ eyink/TurbulenceIII/, course notes, The Johns Hopkins University
work page 2014
-
[19]
Eyink, G. L. & Sreenivasan, K. R. 2006 Onsager and the theory of hydrodynamic turbulence . Rev. Mod. Phys. 78 (1), 87--135
work page 2006
-
[20]
Funakoshi, S.i , Sato, T. & Miyazaki, T. 2012 Maximum entropy state of the quasi-geostrophic bi-disperse point vortex system: Bifurcation phenomena under periodic boundary conditions . Fluid Dyn. Res. 44 (3), 031407
work page 2012
-
[21]
Gallet, B. 2024 Two-dimensional turbulence above topography: Condensation transition and selection of minimum enstrophy solutions . J. Fluid Mech. 988 , A13
work page 2024
- [22]
-
[23]
Hill, J. A. 2017 Jupiter's string of pearls . Caltech Mag. 80 (1), 38--38
work page 2017
-
[24]
Jim \'e nez, J. , Moffatt, H. K. & Vasco, C. 1996 The structure of the vortices in freely decaying two-dimensional turbulence . J. Fluid Mech. 313 , 209--222
work page 1996
-
[25]
Joyce, G. & Montgomery, D. 1973 Negative temperature states for the two-dimensional guiding-centre plasma . J. Plasma Phys. 10 (1), 107--121
work page 1973
-
[26]
2007 Generation and Stability of a Quasi-Permanent Vortex in the Lofoten Basin
K \"o hl, A. 2007 Generation and Stability of a Quasi-Permanent Vortex in the Lofoten Basin . J. Phys. Oceanogr. 37 (11), 2637--2651
work page 2007
-
[27]
LaCasce, J. H. , Pal \'o czy, A. & Trodahl, M. 2024 Vortices over bathymetry . J. Fluid Mech. 979 , A32
work page 2024
-
[28]
Lucarini, V. , Blender, R. , Herbert, C. , Ragone, F. , Pascale, S. & Wouters, J. 2014 Mathematical and physical ideas for climate science . Rev. Geophys 52 (4), 809--859
work page 2014
-
[29]
Matthaeus, W.H. , Stribling, W.T. , Martinez, D. , Oughton, S. & Montgomery, D. 1991 a\/ Decaying, two-dimensional, Navier-Stokes turbulence at very long times . Physica D 51 (1-3), 531--538
work page 1991
-
[30]
Matthaeus, W. H. , Stribling, W. T. , Martinez, D. , Oughton, S. & Montgomery, D. 1991 b\/ Selective decay and coherent vortices in two-dimensional incompressible turbulence . Phys. Rev. Lett. 66 (21), 2731--2734
work page 1991
-
[31]
Mcwilliams, J. C. 1984 The emergence of isolated coherent vortices in turbulent flow . J. Fluid Mech. 146 , 21--43
work page 1984
-
[32]
Mcwilliams, J. C. 1990 The vortices of two-dimensional turbulence . J. Fluid Mech. 219 (-1), 361
work page 1990
-
[33]
1990 Statistical mechanics of Euler equations in two dimensions
Miller, J. 1990 Statistical mechanics of Euler equations in two dimensions . Phys. Rev. Lett. 65 (17), 2137--2140
work page 1990
-
[34]
Modin, K. & Viviani, M. 2020 A Casimir preserving scheme for long-time simulation of spherical ideal hydrodynamics . J. Fluid Mech. 884 , A22
work page 2020
-
[35]
Montgomery, D. & Joyce, G. 1974 Statistical mechanics of “negative temperature” states . Phys. Fluids 17 (6), 1139--1145
work page 1974
-
[36]
Montgomery, D. , Matthaeus, W. H. , Stribling, W. T. , Martinez, D. & Oughton, S. 1992 Relaxation in two dimensions and the ``sinh- Poisson '' equation . Phys. Fluids A 4 (1), 3--6
work page 1992
-
[37]
1949 Statistical hydrodynamics
Onsager, L. 1949 Statistical hydrodynamics . Nuovo Cim 6 (2), 279--287
work page 1949
-
[38]
1991 A maximum-entropy principle for two-dimensional perfect fluid dynamics
Robert, R. 1991 A maximum-entropy principle for two-dimensional perfect fluid dynamics . J Stat. Phys. 65 (3), 531--553
work page 1991
-
[39]
Robert, R. & Sommeria, J. 1991 Statistical equilibrium states for two-dimensional flows . J. Fluid Mech. 229 (-1), 291
work page 1991
- [40]
-
[41]
Salmon, R. , Holloway, G. & Hendershott, M. C. 1976 The equilibrium statistical mechanics of simple quasi-geostrophic models . J. Fluid Mech. 75 (4), 691--703
work page 1976
-
[42]
Santangelo, P. , Benzi, R. & Legras, B. 1989 The generation of vortices in high-resolution, two-dimensional decaying turbulence and the influence of initial conditions on the breaking of self-similarity . Phys. Fluids A 1 (6), 1027--1034
work page 1989
- [43]
-
[44]
Siegelman, L. & Young, W. R. 2023 Two-dimensional turbulence above topography: Vortices and potential vorticity homogenization . PNAS 120 (44), e2308018120
work page 2023
-
[45]
Siegelman, L. , Young, W. R. & Ingersoll, A. P. 2022 Polar vortex crystals: Emergence and structure . Proc. Natl. Acad. Sci. USA. 119 (17), e2120486119
work page 2022
-
[46]
Solodoch, A. , Stewart, A. L. & McWilliams, J. C. 2021 Formation of Anticyclones above Topographic Depressions . J. Phys. Oceanogr. 51 (1), 207--228
work page 2021
-
[47]
Sommeria, J. , Staquet, C. & Robert, R. 1991 Final equilibrium state of a two-dimensional shear layer . J. Fluid Mech. 233 , 661--689
work page 1991
-
[48]
Thompson, A. F. 2010 Jet formation and evolution in baroclinic turbulence with simple topography . J. Phys. Oceanogr. 40 (2), 257--278
work page 2010
-
[49]
Vallis, G. K. & Maltrud, M. E. 1993 Generation of mean flows and jets on a beta plane and over topography . J. Phys. Oceanogr. 23 (7), 1346--1362
work page 1993
-
[50]
Zhang, L.-F. & Xie, J.-H. 2024 Spectral condensation and bidirectional energy transfer in quasi-geostrophic turbulence above small-scale topography . Phys. Fluids 36 (8), 086601
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.