---
title: Lorenz-84 Attractor Analysis
url: https://www.emergentmind.com/topics/lorenz-84-attractor
type: topic
---

# Lorenz-84 Attractor Analysis

Searching arXiv for recent papers on the Lorenz-84 attractor and related context.
The Lorenz–84 attractor is the chaotic invariant set associated with the Lorenz–84 system, a low-order model of large-scale atmospheric dynamics introduced by Lorenz. In the parameter regime \(a=0.25\), \(b=4.0\), \(F=8.0\), and \(G=1.0\), it is described as a weakly dissipative chaotic attractor whose geometry is substantially thicker than canonical strongly dissipative examples such as Lorenz–63 or Rössler. Recent structural analysis identifies it as a nontrivial case of toroidal chaos organized around a period-2 cavity and generated by a multidirectional stretching mechanism, while later work on seasonally forced and noise-perturbed variants recasts the attractor as a non-autonomous, stochastic regime structure rather than a single static invariant set [2507.04921].

## 1. Governing equations and phase-space conventions

The autonomous Lorenz–84 model studied in the structural analysis is
\[
\left\{
\begin{array}{l}
\dot{x} = - y^2 - z^2 -a x + aF \\
\dot{y} = xy - bxz - y + G \\
\dot{z} = bxy +xz -z \;.
\end{array}
\right.
\]
Here \(a\) is the damping or relaxation parameter in the \(x\)-equation, \(b\) is the coupling parameter between \(x,y,z\), \(F\) is the external forcing in the \(x\)-equation, and \(G\) is the forcing in the \(y\)-equation. In the later seasonally forced formulation, \(x\) is the strength of the zonal jet stream, while \(y\) and \(z\) are the amplitudes of the cosine and sine phases of planetary eddies; terms such as \(xy\) and \(xz\) represent exchange of energy between the jet and eddies, \(-y^2-z^2\) represents jet energy loss to eddies, and \(-bxz\) and \(bxy\) represent advection or displacement of eddies by the mean flow [2510.26005].

For topological analysis of the autonomous chaotic regime, the flow is rotated. First, a sign change \((x',y',z')=(-x,y,-z)\) is applied, then a rotation by \(\theta=1.5\) rad in the \((y',z')\)-plane:
\[
\left\{
\begin{aligned}
x_r &= -x \\
y_r &= y \cos \theta + z \sin \theta \\
z_r &= y \sin \theta - z \cos \theta .
\end{aligned}
\right.
\]
In these rotated coordinates, the flow evolves clockwise in the chosen projection. The Poincaré section is
\[
\mathcal{P} = \{ (x_{r,n}, z_{r,n}) ~ | ~ y_{r,n} = 0, ~\dot{y}_{r} < 0\}.
\]
A central fact for the Lorenz–84 attractor is that this section is two-dimensional and thick rather than effectively one-dimensional [2507.04921].

## 2. Weak dissipation and the limits of classical structure analysis

The Lorenz–84 attractor at \(a=0.25\), \(b=4.0\), \(F=8.0\), \(G=1.0\) is already known to be weakly dissipative, with Kaplan–Yorke dimension around \(D_{KY}\sim 2.39\). This places it in a different geometric class from nearly two-dimensional, strongly dissipative attractors. In the classical setting, strong transverse contraction makes the Poincaré section thin, permits parameterization by a single scalar coordinate, and yields a first return map that is a thin one-dimensional curve. The standard route is then thin section \(\to\) scalar return map \(\to\) 2D template.

For Lorenz–84, that route becomes ambiguous because the Poincaré section is genuinely bidimensional, scalar projection loses the one-to-one property, return-map branches are thick rather than sharp, and several stretching directions seem to coexist. The paper illustrates this with the scalar coordinate
\[
\rho_n = \frac{-z_r-\min(z_r)}{\max(z_r)-\min(z_r)},
\]
for which the first return map is not bijective enough to support standard analysis. A more informative but still artificial partition is introduced through
\[
\gamma_n =
\left\{
\begin{aligned}
&\frac{-z_r-\min(z_r)}{1.5-\min(z_r)} && \text{if } z_r \geq -1.5 \text{ and } x_r \leq -0.67 \\
&\frac{-z_r-\min(z_r)}{\max(z_r)-\min(z_r)} + 1 && \text{else}.
\end{aligned}
\right.
\]
Even this construction does not yield a direct classical template derivation because thickness remains essential rather than negligible [2507.04921].

A frequent misconception is to treat all low-dimensional chaotic attractors as if they admitted the same thin-template reduction. The Lorenz–84 case shows that weak dissipation can make the section itself part of the problem: the section is not merely noisy around a one-dimensional skeleton, but structurally bidimensional. This suggests that the failure of classical methods is not accidental but tied to the attractor’s volumetric organization.

## 3. Three-dimensional organization: toroidal chaos around a period-2 cavity

The principal structural result is that the Lorenz–84 attractor is a nontrivial toroidal chaotic attractor organized around an internal cavity associated with a period-2 orbit [2507.04921]. The authors observed an apparently empty region in the Poincaré section, with periodic points arranged around it, then ran simulations from that region, found trajectories remaining there for a very long transient, and extracted the organizing period-2 orbit by optimization. Escape from this region to the chaotic attractor required 361 returns, and once trajectories reached the chaotic attractor they never returned to the hole region.

This establishes an internal cavity that lies inside the attractor’s bounding volume and belongs to the basin of attraction, but is not part of the attractor itself. Topologically, the attractor is argued to be externally bounded by a genus-1 torus and internally pierced by a genus-1 cavity of period 2. The resulting organization is therefore toroidal, but not of the simplest type centered on a period-1 core. The authors explicitly state that this appears to be the first revealed case of toroidal chaos organized around a cavity of period 2 [2507.04921].

The reconstructed geometry shows an alternating left-right and right-left circulation around this cavity. The flow stretches and folds around the excluded region, with branches on one side mapping to the other and back again. The period-1 orbit also exists, but the internal organization is not described as a genus-1 torus of period 1 because that region is too densely occupied. A plausible implication is that the topological role of “the hole” in Lorenz–84 cannot be inferred from low-period periodic orbits alone; it depends on the occupancy structure of the surrounding chaotic layers.

## 4. Color tracer mapping and multidirectional stretching

To recover the attractor’s three-dimensional organization, the structural study introduces color tracer mapping. Instead of collapsing the section dynamics to a scalar return map, the method assigns colors to a region of the section and tracks those colors under the return dynamics. Formally, the paper defines a discrete map \(g\),
\[
g \colon \mathbb{R}^n ~~ \to ~~ \mathbb{R}^n,
\]
a color palette map \(C^\alpha\),
\[
C^{\alpha} \colon \mathbb{R}^n ~~ \to ~~ \mathbb{R}^3,
\]
and a color map
\[
\Gamma_p^{\alpha} \colon \mathbb{R}^{n+3} ~~ \to ~~ \mathbb{R}^{n+3},
\]
with
\[
\{\vec{X}_j, C^{\alpha}(\vec{X}_j) \} ~~ \mapsto ~~ \{\vec{X}_{j+p}, C^{\alpha}(\vec{X}_j) \},
\]
where \(p\) is the return number. In practice, a localized region of the Poincaré section is colored with a vertical gradient, propagated one step backward and one step forward, and compared across three consecutive iterations \(\Gamma^{\alpha}_{n-1},\Gamma^{\alpha}_n,\Gamma^{\alpha}_{n+1}\), after which a continuous deformation is manually reconstructed [2507.04921].

The method reveals how local strips are stretched, where they fold, whether continuity is preserved, how branches split and recombine, and how the attractor wraps in three dimensions. The reported Lorenz–84 outcome is continuous deformation with no tearing, and a structure reconstructible modulo \(2\pi\), with the remaining angular ambiguity resolved using the period-2 cavity and its positive twist.

The same reconstruction supports the paper’s second major claim: the chaos-generation mechanism is multidirectional stretching rather than simple one-directional stretching and folding. A single branch in the upper-left part of the section is stretched and folded so as to generate three separate branches around the central cavity. In the \(\gamma_n\)-based return map, the dynamics can be partitioned into six branches \(A,B,D,C,E,F\), with \(A,B\) mapping part I to part II, \(D,C\) mapping part I to part I, and \(E,F\) mapping part II to part I. Part I therefore stretches and folds into two destinations at once: one back to itself and one to the other part. The authors emphasize this as a visible signature of multidirectional stretching and suggest that it may be a general feature, perhaps even a necessary condition, for weakly dissipative chaos [2507.04921].

## 5. Branched-manifold representation and topological validation

Despite the attractor’s thickness, the paper constructs a reduced branched-manifold representation. The reconstruction begins with a 3D skeleton containing initially nine branches, with branches \(c\)-\(c'\) and \(d\)-\(d'\) kept separate because they mix with different branches. A simplified skeleton is then obtained by removing short branches not structurally necessary, emphasizing the period-6 orbit embedded in the skeleton, the period-2 cavity in the middle, and the branch architecture responsible for stretching and folding. The authors finally present an artificial flattened representation: a two-dimensional branched manifold embedded through the flow, not a standard thin-attractor template but a reduced topological model [2507.04921].

Two related representations are distinguished. The more geometric squeezed template has seven branches \(a,b,c,d,e,f,g\), including split branches \(c\!-\!c'\) and \(d\!-\!d'\). The more algebraic reduced template has six branches \(A,B,D,C,E,F\), designed for symbolic and topological validation. It retains branch ordering, branch parity, symbolic coding of periodic orbits, torsions and permutations via a linking matrix, and the special role of the period-2 cavity. The linking matrix is
\[
\begin{array}{cc}
& \begin{matrix} \;\;\;A & B & D & C & E & F \end{matrix} \\
\begin{matrix} A \\ B \\ D \\ C \\ E \\ F \end{matrix}
&
\left.
\begin{bmatrix}
\;0 & \;0  & 1  & 1  & 1 & 1 \\
\;0 & 1 & 1 & 1 & 1 & 1 \\
1 & 1 & \;0 & \;0 & \;0 & \;0 \\
1 & 1 & \;0 & 1 & 1 & 1 \\
1 & 1 & \;0 & 1 & \;0 & \;0 \\
1 & 1 & \;0 & 1 & \;0 & 1
\end{bmatrix}
\!\!\!\!
\right|
\end{array}.
\]

Validation uses unstable periodic orbits extracted directly from the Poincaré section because the \(\rho_n\)-based scalar return map is not sufficiently bijective. The extracted set consists of one period-1 orbit, three period-5 orbits, one period-6 orbit, and additionally a period-2 orbit inside the cavity. Their symbolic codes in the six-branch partition are
\[
p1: C,\qquad p2: BE,
\]
\[
p5\text{-A}: AFD_3CD_5,\qquad
p5\text{-B}: EB_2FDB_5,\qquad
p5\text{-C}: DBFC_4C_5,
\]
\[
p6: B_1E_2B_3E_4CD.
\]
The reported numerical linking numbers include
\[
lk(p1,p5\text{-A})=2,\quad
lk(p1,p5\text{-B})=2,\quad
lk(p1,p5\text{-C})=2,\quad
lk(p1,p6)=3,\quad
lk(p1,p2)=1,
\]
\[
lk(p5\text{-A},p5\text{-B})=10,\quad
lk(p5\text{-A},p5\text{-C})=10,\quad
lk(p5\text{-B},p5\text{-C})=10,
\]
\[
lk(p5\text{-A},p6)=12,\quad
lk(p5\text{-B},p6)=12,\quad
lk(p5\text{-C},p6)=12,
\]
\[
lk(p5\text{-A},p2)=4,\quad
lk(p5\text{-B},p2)=4,\quad
lk(p5\text{-C},p2)=4,\quad
lk(p6,p2)=5.
\]
Theoretical linking numbers are computed from crossings induced by initial branch separation \(N_{\text{sep}}\), torsions and permutations from the linking matrix \(M_{ij}\), and insertion or squeezing crossings \(N_{\text{ins}}\), with the linking number equal to half the total crossing count. The agreement between theoretical and numerical linking numbers is presented as the main evidence that the proposed branched-manifold description captures the attractor’s topology correctly [2507.04921].

The authors also record several limitations. Thickness-induced ambiguity makes branch boundaries approximate; the template is validated but may not be unique; non-conventional insertion crossings probably arise from flattening a thick structure as if it were dissipative; and placement of longer-period orbits may require refinement. The resulting template is therefore a validated first structural and topological model rather than a final exact flattening.

## 6. Non-autonomous, seasonal, and stochastic Lorenz–84 attractors

A later study examines a seasonally forced and noise-perturbed Lorenz–84 model in which
\[
\frac{dx}{dt} = -y^{2} - z^{2} - ax + aF(t), \qquad
\frac{dy}{dt} = xy - bxz - y + G, \qquad
\frac{dz}{dt} = bxy + xz - z,
\]
with
\[
F(t) = F_{0} + F_{1} \cos(\omega t),
\]
where \(F_0=7\) and \(F_1=2\). The stochastic perturbation is introduced in discrete-time form as
\[
\mathbf{x}_{t} = f(\mathbf{x}_{t-1}) + \varepsilon_{t}, \qquad
\varepsilon_t \sim \mathcal{N}\!\left(0,\, m \cdot \lvert \mathbf{x}_{t-1} \rvert \right).
\]
In this setting, the object under study is no longer the simple invariant chaotic attractor familiar from autonomous Lorenz systems, but a more complex non-autonomous, noise-perturbed attractor-like structure whose geometry and visitation statistics vary with season [2510.26005].

The paper describes this regime structure in terms of broad exploration of phase space, frequent irregular excursions, high-eddy-amplitude regimes, latent hidden states occupying different regions of state space, and one particularly unstable dangerous regime. The principal observable is the eddy amplitude
\[
|Y|+|Z|,
\]
with a threshold of \(2.4\) shown in the main figures as a high-eddy threshold; the appendix also uses the 97th percentile of \(|Y|+|Z|\), while the local Lyapunov exponent setup discusses the 90th percentile as a threshold for high-eddy regimes in long simulations. For fixed forcing values \(F=5,6,7,8\), larger \(F\), especially \(F=8\), produces the most chaotic behavior, with stronger and more frequent high-eddy excursions. Under continuous \(F(t)\), trajectories move through a family of seasonally varying regimes rather than a fixed attractor [2510.26005].

The paper analyzes these regimes through local Lyapunov exponents and a non-homogeneous Hidden Markov Model. The NHMM transition probabilities are
\[
P(S_{t+1} = j \mid S_t = i,\, C(t)) = P_{ij}(t)
\]
with
\[
P_{ij}(t) = \frac{\exp(\beta_{ij,0} + \beta_{ij,1} C(t))}
{\sum_{k=1}^K \exp(\beta_{ik,0} + \beta_{ik,1} C(t))},
\]
where the covariate \(C(t)\) is the seasonal forcing. Model selection by BIC chose a 9-state specification. State-conditioned distributions in the \((y,z)\) plane occupy identifiable regions with distinct density patterns, and one latent state, state 4, is singled out as the most unstable part of the attractor-like structure: its local Lyapunov exponent density is sharply peaked above \(3.5\), most other states are centered between \(1.5\) and \(2.5\), and state 4 accounts for more than \(90.5\%\) of all exceedances above the 97th percentile threshold [2510.26005].

This reframes the Lorenz–84 attractor in the non-autonomous stochastic setting as a seasonally evolving stochastic regime landscape. The paper is explicit that one should be cautious about speaking as if there were a single static invariant attractor in the classical autonomous sense. A common misunderstanding is therefore to transfer the autonomous notion of “the attractor” unchanged to the seasonally forced noisy model; the paper instead treats the relevant object through trajectories, thresholds, latent states, local instability, and seasonally modulated transition probabilities.

## 7. Relation to Lorenz-type attractors and rigorous dynamical context

The Lorenz–84 attractor is frequently compared with Lorenz-type chaotic sets, but the available material distinguishes sharply between direct analysis of Lorenz–84 and rigorous theory developed for other systems. A computer-assisted proof for the Shimizu–Morioka system establishes the existence of a Lorenz attractor for an open set of parameter values and clarifies what mathematicians mean by a Lorenz-type attractor: a geometric singular hyperbolic or pseudo-hyperbolic object organized by a saddle equilibrium, homoclinic bifurcation structure, and area expansion along separatrices [1711.10404].

That paper does not discuss the Lorenz–84 model specifically, yet it provides transferable context. In particular, it emphasizes pseudo-hyperbolicity, the Shilnikov criterion based on a symmetric homoclinic butterfly with zero saddle value and separatrix value \(0<|A|<2\), and the distinction between robust chaos and full structural stability. This does not establish that the Lorenz–84 attractor itself is a Lorenz attractor in the geometric singular-hyperbolic sense. Instead, it supplies a rigorous framework for understanding what would count as a Lorenz-type attractor and how such behavior can be certified on an open parameter region [1711.10404].

Within that broader context, the Lorenz–84 attractor occupies a distinctive position. The structural analysis portrays it neither as a standard strongly dissipative template-type attractor nor as a simple toroidal chaos organized around a period-1 hole. It is instead a thick toroidal chaos with a period-2 cavity and multidirectional stretching, and in its seasonally forced noisy extension it becomes a regime-partitioned, seasonally modulated stochastic structure. A plausible implication is that Lorenz–84 sits at an intersection of several traditions in dynamical-systems research: low-order atmospheric modeling, topological analysis of three-dimensional chaos, and probabilistic coarse-graining of non-autonomous stochastic dynamics.

Source: https://www.emergentmind.com/topics/lorenz-84-attractor