arxiv: 2604.14123 · v1 · submitted 2026-04-15 · ⚛️ physics.ao-ph

Recognition: unknown

Revisiting the Dynamical Properties of Pedlosky's Two-Layer Model for Finite Amplitude Baroclinic Waves

Nicolas De Ro , Jonathan Demaeyer , St\'ephane Vannitsem

Authors on Pith no claims yet

Pith reviewed 2026-05-10 11:23 UTC · model grok-4.3

classification ⚛️ physics.ao-ph

keywords baroclinic wavestwo-layer modeldeterministic chaosbifurcation analysisEkman dissipationPedlosky modelLorenz-like truncation

0 comments

The pith

Pedlosky's two-layer baroclinic wave model becomes chaotic through period-doubling as Ekman dissipation rises from an integrable inviscid state.

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

The paper examines the dynamical regimes of finite-amplitude baroclinic waves in Pedlosky's two-layer model, focusing on the role of dissipation from Ekman boundary layers. It establishes that the model is integrable without dissipation. With increasing dissipation, a sequence of bifurcations occurs, leading to periodic oscillations, homoclinic orbits, and deterministic chaos identified via Lyapunov exponents in a low-order truncation equivalent to the Lorenz model. In equilibrating regimes, multiple attractors coexist depending on initial conditions. The study also assesses the robustness of these behaviors in higher-dimensional versions of the model.

Core claim

The geophysical state of the baroclinic wave in Pedlosky's model shows a rich diversity of dynamical regimes controlled by Ekman dissipation. The inviscid limit is integrable, while increasing dissipation causes bifurcations to deterministic chaos and coexistence of multiple attractors. The minimal low-order truncation, structurally equivalent to the Lorenz model, reveals the transition from stable wave equilibration to periodic oscillations terminating in homoclinic orbits and then to chaos via period-doubling.

What carries the argument

The bifurcation diagram of the low-order truncated system of ODEs, which is structurally equivalent to the Lorenz model and controls the wave amplitude transitions.

If this is right

Deterministic chaos induced by dissipation provides a genuine mechanism for destabilization of the baroclinic wave.
Initial-condition dependence allows multiple attractors to coexist in regimes where the wave equilibrates.
The sequence of bifurcations and the period-doubling route remain robust when the model is extended to higher dimensions.

Where Pith is reading between the lines

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

Atmospheric wave prediction may require explicit accounting for initial-condition sensitivity even at moderate dissipation levels.
Similar bifurcation structures could govern baroclinic instability in three-layer or continuously stratified models.
Climate simulations that parameterize Ekman friction might exhibit abrupt shifts between wave states as friction coefficients vary.

Load-bearing premise

The low-order truncation of the governing equations is sufficient to capture the essential nonlinear dynamics of the full two-layer system.

What would settle it

Direct numerical integration of the full two-layer partial differential equations showing no period-doubling cascade to chaos or absence of multiple coexisting attractors at the same dissipation parameter values.

Figures

Figures reproduced from arXiv: 2604.14123 by Jonathan Demaeyer, Nicolas De Ro, St\'ephane Vannitsem.

**Figure 2.** Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Trajectories of the Hamiltonian system ( [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Stability of the two fixed points [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. Basins of attraction of the toy model ( [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6. Solutions of the toy model ( [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗

**Figure 7.** Figure 7: FIG. 7. Maximal LCE of the toy model ( [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗

**Figure 8.** Figure 8: FIG. 8. Typical nature of the toy model’s asymptotic dynamics as [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗

**Figure 9.** Figure 9: FIG. 9. Some of the stable periodic orbits of the toy model listed in Tab. [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗

**Figure 10.** Figure 10: FIG. 10. Bifurcation diagram of the toy model ( [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗

**Figure 11.** Figure 11: FIG. 11. Continuation of the first stable periodic orbit found at [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗

**Figure 12.** Figure 12: FIG. 12. Poincaré sections on the plane [PITH_FULL_IMAGE:figures/full_fig_p014_12.png] view at source ↗

**Figure 13.** Figure 13: FIG. 13. Basins of attraction of the model ( [PITH_FULL_IMAGE:figures/full_fig_p017_13.png] view at source ↗

**Figure 14.** Figure 14: FIG. 14. Maximal LCE of the model ( [PITH_FULL_IMAGE:figures/full_fig_p017_14.png] view at source ↗

**Figure 15.** Figure 15: FIG. 15. Bifurcation diagram of the model ( [PITH_FULL_IMAGE:figures/full_fig_p018_15.png] view at source ↗

**Figure 16.** Figure 16: FIG. 16. Continuation of the first stable periodic orbit found at [PITH_FULL_IMAGE:figures/full_fig_p018_16.png] view at source ↗

**Figure 17.** Figure 17: FIG. 17. Absolute difference between maximal Lyapunov exponents [PITH_FULL_IMAGE:figures/full_fig_p019_17.png] view at source ↗

**Figure 18.** Figure 18: FIG. 18. Time evolution of a trajectory from initial condition [PITH_FULL_IMAGE:figures/full_fig_p023_18.png] view at source ↗

read the original abstract

Baroclinic instability is a fundamental mechanism driving atmospheric dynamics. In this work, we revisit Pedlosky's two-layer model for finite amplitude baroclinic waves - a seminal framework for studying the unstable growth of finite perturbations - leveraging modern nonlinear techniques and computational resources. We show that the geophysical state of the baroclinic wave exhibits a rich diversity of dynamical regimes governed by the level of dissipation induced by Ekman boundary layers. In the inviscid limit, we demonstrate that the model is integrable. Upon increasing dissipation, the system undergoes a complex sequence of bifurcations. On one hand, deterministic chaos, identified by means of the Lyapunov exponents, provides a genuine mechanism for destabilization of the wave. On the other hand, in regimes where the wave equilibrates, dependence on the initial condition is crucial, eventually leading to the coexistence of multiple attractors. We study the governing equations of the model and their truncation to a finite-dimensional system of ordinary differential equations, together with the minimal low-order truncated system which is structurally equivalent to the Lorenz model. Its bifurcation diagram allows for elucidating the transition of the wave amplitude from stable equilibration to periodic oscillations - terminating in homoclinic orbits - and, ultimately, deterministic chaos through a period-doubling route. We finally comment on the robustness of these features for higher-dimensional models.

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.

Desk Editor's Note private letter to a colleague

This paper maps out period-doubling to chaos and coexisting attractors in the low-order truncation of Pedlosky's model, with a nod to higher truncations, but the quantitative link between them stays thin.

read the letter

This paper takes Pedlosky's two-layer model for finite-amplitude baroclinic waves and works out its bifurcation structure under varying Ekman dissipation using modern tools. The key new pieces are the demonstration of integrability in the inviscid limit, the period-doubling route to chaos identified via Lyapunov exponents, the termination in homoclinic orbits, and the coexistence of multiple attractors in certain dissipation ranges. They also note that the minimal truncation is structurally equivalent to the Lorenz system and discuss how these behaviors appear in higher truncations. The work does a good job applying Lyapunov analysis and bifurcation diagrams to this specific model. It is a direct numerical study of the truncated ODEs, so the results are reproducible in principle from the equations given. The abstract indicates they examine both the full governing equations and the low-order system, which is the right approach. The main limitation is that the headline claims about the geophysical baroclinic wave rest on the truncation. While the paper comments on robustness for higher-dimensional versions, the stress-test concern holds some weight: there is no clear quantitative evidence in the provided summary that the attractor structure or exponent values converge as more modes are added. If higher modes alter the bifurcation points or introduce new instabilities, the reported diversity of regimes could be truncation-dependent. That said, for a low-order model study this is not unusual, and the equivalence to Lorenz is a known feature. Readers working on atmospheric predictability or low-order GFD models will find the explicit regimes and diagrams useful. It is incremental rather than foundational, but the details on multiple attractors and the dissipation dependence are worth having on record. I recommend sending this to peer review. The core calculations are straightforward and the topic is standard, so referees can assess the truncation validity and the numerical methods without much trouble.

Referee Report

2 major / 2 minor

Summary. The manuscript revisits Pedlosky's two-layer model for finite-amplitude baroclinic waves. It demonstrates integrability of the governing equations in the inviscid limit and, with increasing Ekman-layer dissipation, a sequence of bifurcations in a low-order truncation that is structurally equivalent to the Lorenz model. This truncation yields stable equilibration, periodic oscillations terminating in homoclinic orbits, a period-doubling route to deterministic chaos (via Lyapunov exponents), and coexistence of multiple attractors. The paper studies the full governing equations and their truncations, and comments on robustness in higher-dimensional models.

Significance. If the truncation analysis holds, the work supplies a concrete bridge between classical baroclinic-instability theory and low-dimensional chaotic dynamics, with explicit identification of integrability and a Lyapunov-based mechanism for wave destabilization. The numerical construction of bifurcation diagrams and exponent spectra constitutes a clear methodological advance over purely analytic treatments of the same model.

major comments (2)

[§3] §3 (low-order truncation): the structural equivalence to the Lorenz equations is asserted without an explicit change-of-variables or coefficient-matching derivation; because the subsequent bifurcation diagram and period-doubling route rest entirely on this equivalence, the mapping must be shown in detail.
[§5] §5 (robustness discussion): the claim that the reported sequence of regimes persists in higher-dimensional truncations is supported only by qualitative remarks; no quantitative diagnostics (convergence of Lyapunov spectra, basin volumes, or attractor statistics between the minimal and next-order truncations) are provided, leaving open the possibility that omitted modes alter the bifurcation structure.

minor comments (2)

[Methods] Numerical details for the Lyapunov-exponent calculation (integration method, step size, transient discard, and convergence criteria) are not stated; these should be supplied so that the reported exponents can be reproduced.
[Figures] Bifurcation diagrams lack explicit annotation of the critical dissipation values at which the homoclinic and period-doubling transitions occur.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major point below and indicate the changes planned for the revised version.

read point-by-point responses

Referee: [§3] §3 (low-order truncation): the structural equivalence to the Lorenz equations is asserted without an explicit change-of-variables or coefficient-matching derivation; because the subsequent bifurcation diagram and period-doubling route rest entirely on this equivalence, the mapping must be shown in detail.

Authors: We agree that the manuscript asserts structural equivalence without supplying the explicit transformation. In the revised version we will insert a dedicated derivation in §3 that performs the change of variables, lists the resulting coefficients, and verifies that they match those of the Lorenz system up to a linear rescaling. This will make the subsequent bifurcation analysis fully traceable. revision: yes
Referee: [§5] §5 (robustness discussion): the claim that the reported sequence of regimes persists in higher-dimensional truncations is supported only by qualitative remarks; no quantitative diagnostics (convergence of Lyapunov spectra, basin volumes, or attractor statistics between the minimal and next-order truncations) are provided, leaving open the possibility that omitted modes alter the bifurcation structure.

Authors: The referee correctly observes that the robustness section remains qualitative. We will augment §5 with quantitative comparisons: Lyapunov exponent spectra and estimates of basin volumes will be computed for both the minimal (3-mode) and next-order (5-mode) truncations at representative dissipation values. These diagnostics will be presented to demonstrate that the sequence of regimes and the period-doubling route remain intact. revision: yes

Circularity Check

0 steps flagged

No circularity: results follow from direct truncation and numerical analysis of the model equations

full rationale

The paper starts from Pedlosky's two-layer equations, performs an explicit truncation to a finite-dimensional ODE system, identifies a minimal truncation that is structurally equivalent to the Lorenz equations by direct reduction of the governing PDEs, and then computes integrability (inviscid limit), Lyapunov exponents, and bifurcation diagrams via standard numerical methods on those ODEs. No parameter is fitted to data and then relabeled as a prediction; no load-bearing uniqueness theorem is imported via self-citation; the truncation is presented as an approximation whose robustness is separately commented upon for higher-dimensional versions. All central claims are therefore obtained from the model's own equations rather than by construction from prior outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only abstract available; no explicit free parameters, axioms, or invented entities are stated beyond the standard assumptions of the Pedlosky two-layer quasigeostrophic setup.

pith-pipeline@v0.9.0 · 5550 in / 1055 out tokens · 43513 ms · 2026-05-10T11:23:51.710043+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

4 extracted references

[1]

Lyapunov vectors and error growth patterns in a t21l3 quasigeostrophic model,

1J. R. Holton and G. J. Hakim,An Introduction to Dynamic Meteorology, fifth edition ed. (Academic Press, Boston, 2013). 2G. K. Vallis,Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation, 2nd ed. (Cambridge University Press, 2017). 3J. Pedlosky,Geophysical Fluid Dynamics(Springer New York, New York, NY , 1987). 4H. Poincaré,Sci...

2013
[2]

Kolmogorov entropy and nu- merical experiments,

Kuznetsov, B. Oldeman, R. Paffenroth, B. Sandstede, X. Wang, and C. Zhang,AUTO-07P: Continuation and Bifurcation Software for Ordinary Differential Equations(2007), software and documentation. 12That was not available at the time of the first attempt to characterize the bifurcation diagram of the model10. 13G. Benettin, L. Galgani, and J.-M. Strelcyn, “Ko...

2007
[3]

York, New York, NY , 1989). 27S. H. Strogatz,Nonlinear Dynamics and Chaos: With Applications to

1989
[4]

auto-auto: A python layer for automati- cally running the auto-07p continuation software,

Physics, Biology, Chemistry, and Engineering, 2nd ed. (CRC Press, Boca Raton, 2015). 28J. Demaeyer and O. Hamilton, “auto-auto: A python layer for automati- cally running the auto-07p continuation software,” Journal of Open Source Software10, 8079 (2025). 29C. Sparrow,The Lorenz Equations: Bifurcations, Chaos, and Strange At- tractors(Springer New York, N...

2015