Generalized Lorenz Model · redesigned note

A Shared Cubic Landscape: From the Generalized Lorenz Model to Lorenz 1963, NLS Amplitudes, and the Cubic EBM

This page reframes the Generalized Lorenz Model (GLM) as a bridge: from the classical Lorenz (1963) model, to simplified non-dissipative Lorenz dynamics, to the amplitude equation of the nonlinear Schrödinger (NLS) equation, and finally to first-order cubic energy-balance models (EBMs).

Main message

A hierarchy connected by reductions, not by analogy alone

The GLM extends Lorenz-type convection models by adding higher-wavenumber temperature modes. The classical Lorenz (1963) model appears as the three-mode core. When selected dissipative terms are removed or simplified, Lorenz-type dynamics reduce to a second-order cubic equation. That same cubic structure appears in the amplitude equation of the nonlinear Schrödinger equation and in the overdamped limit becomes a first-order cubic EBM.

1

GLM

Higher-dimensional Lorenz models with added Fourier mode pairs.

2

Lorenz 1963

The three-mode core recovered when added modes are omitted.

3

Non-dissipative Lorenz

A second-order cubic reduction with centers, oscillations, and solitary-wave solutions.

4

NLS amplitude

Traveling-wave envelopes satisfy a related second-order cubic ODE.

5

Cubic EBM

The overdamped first-order limit shares the same potential landscape.

GLM and the classical core

From the generalized Lorenz model to Lorenz (1963)

The generalized Lorenz model may be written as a core Lorenz subsystem coupled to additional mode pairs \((Y_j,Z_j)\). For \(j=1,\dots,N\), where \(N=(M-3)/2\), one representative form is:

GLM core and added modes
\[ \begin{aligned} \frac{dX}{d\tau} &= \sigma(Y-X),\\ \frac{dY}{d\tau} &= -XZ+rX-Y,\\ \frac{dZ}{d\tau} &= XY-XY_1-bZ,\\ \frac{dY_j}{d\tau} &= jXZ_{j-1}-(j+1)XZ_j-d_{j-1}Y_j,\\ \frac{dZ_j}{d\tau} &= (j+1)XY_j-(j+1)XY_{j+1}-\beta_j Z_j . \end{aligned} \]

Here \(X\) is the stream-function amplitude and \((Y,Z),(Y_1,Z_1),\ldots\) are temperature-mode amplitudes at increasing wavenumbers. The extra nonlinear term \(-XY_1\) in the \(Z\)-equation represents the influence of added Fourier modes.

Lorenz (1963) three-mode core
\[ \frac{dX}{d\tau}=\sigma(Y-X),\qquad \frac{dY}{d\tau}=-XZ+rX-Y, \qquad \frac{dZ}{d\tau}=XY-bZ. \]

The L63 model is therefore not separate from the GLM; rather, it is the three-dimensional truncation obtained by dropping the added modes and the feedback they introduce. When the linear dissipative terms are ignored, the L63 model becomes nondissipative, as shown below.

Second-order reduction

The cubic oscillator equation and the non-dissipative Lorenz model

Near the unstable critical point, the dynamics of a pendulum is governed by a damped cubic oscillator equation. This second-order equation serves as the organizing bridge between the nondissipative Lorenz model and the bulk Energy Balance Model (EBM):

Damped cubic oscillator equation
\[ \frac{d^2X}{dt^2}+\mu\frac{dX}{dt}=X-\frac{X^3}{6}. \]

When \(\mu=0\), the damping vanishes and the equation becomes the non-dissipative Lorenz model in reduced second-order form. It preserves inertia and permits oscillatory or solitary-wave behavior while retaining the same cubic force term.

For \(\mu\ne0\), if the second-derivative term is ignored, the same equation gives the overdamped first-order form:

Overdamped first-order limit
\[ \mu\frac{dX}{dt}=X-\frac{X^3}{6}, \]

which is the unforced prototype of the first-order cubic EBM. Thus the pathway is: damped cubic oscillator \(\rightarrow\) non-dissipative Lorenz model when \(\mu=0\), and damped cubic oscillator \(\rightarrow\) first-order EBM when \(\mu\ne0\) and the inertial term is neglected.

Interpretation. The non-dissipative Lorenz model is not a climate EBM yet; it is a conservative, second-order dynamical system. The first-order EBM is obtained by removing inertia, but both inherit the same cubic force term and therefore the same unforced potential landscape.
Wave-envelope connection

The NLS amplitude equation shares the same cubic skeleton

For the nonlinear Schrödinger equation, a traveling-wave ansatz \(\psi(x,t)=h(\xi)e^{i(kx-\omega t)}\), \(\xi=x-ct\), yields a real amplitude equation of the form:

NLS traveling-wave amplitude
\[ h''(\xi)+\delta h(\xi)+\gamma h^3(\xi)=0, \qquad \delta=2\omega-k^2,\quad \gamma=-2\kappa. \]

For the focusing case and appropriate signs of \(\delta\) and \(\gamma\), this amplitude equation is mathematically equivalent to the non-dissipative Lorenz cubic reduction after rescaling. This is the bridge from Lorenz-type amplitude dynamics to solitary-wave and envelope dynamics.

Overdamped limit

The first-order cubic EBM keeps the landscape and removes inertia

Starting from the same damped cubic oscillator equation, keep \(\mu\ne0\) but neglect the inertial term \(d^2X/dt^2\). The result is the overdamped first-order system. Adding external forcing \(F\) gives the cubic EBM:

Cubic energy-balance model
\[ \dot X = F+aX-bX^3, \qquad a>0,\quad b>0. \]

Here \(X\) may be interpreted as a temperature perturbation or regime-state variable, while \(F\) tilts the potential. For a range of \(F\), the cubic EBM admits two stable equilibria and one unstable equilibrium, capturing bistability, hysteresis, and tipping behavior.

Second-order non-dissipative Lorenz form

Has inertia. Stable points behave as centers in the conservative limit. Oscillation and solitary-wave interpretations are possible.

First-order cubic EBM

No inertia. Motion is gradient-like: trajectories move down the tilted potential toward stable climate states.

Shared potential landscape

The key common structure is the double-well potential

Both the second-order cubic Lorenz reduction and the first-order cubic EBM can be analyzed using the same potential function. Defining the force term as \(H(X;F)=F+aX-bX^3=-\partial V/\partial X\), the potential is:

Potential function
\[ V(X;F)=-FX-\frac{a}{2}X^2+\frac{b}{4}X^4. \]

For \(F=0\), the two minima represent stable regimes and the central maximum represents an unstable barrier. Increasing \(F\) tilts the landscape, shifting equilibrium locations and eventually producing a saddle-node bifurcation: the geometric onset of tipping.

F = 0.00This diagram is drawn locally with plain JavaScript; it does not use MathJax or external plotting libraries.

Same landscape

The second-order and first-order equations share \(V(X;F)\). Stability of equilibria is determined by the curvature \(V_{XX}\): minima are stable, maxima are unstable.

Different motion on the landscape

The second-order system can oscillate across the landscape; the first-order EBM descends it. The geometry is shared, but the time evolution differs.

Selected references and links
  1. Lorenz, E. N. (1963). Deterministic nonperiodic flow. Journal of the Atmospheric Sciences, 20, 130–141.
  2. Shen, B.-W. (2019). Aggregated Negative Feedback in a Generalized Lorenz Model. International Journal of Bifurcation and Chaos. https://doi.org/10.1142/S0218127419500378
  3. Shen, B.-W. (2025). From Airy’s Equation to the Non-Dissipative Lorenz Model: Turning Points, Quantum Tunneling, and Solitary Waves. Encyclopedia. MDPI page
  4. Shen, B.-W. (2026). Tipping, Nonlinear Feedback, and Bistability in a Cubic Energy Balance Model. International Journal of Bifurcation and Chaos. DOI: 10.1142/S0218127426501257.
  5. Shen et al. (2022). The Dual Nature of Chaos and Order in the Atmosphere. Atmosphere. https://doi.org/10.3390/atmos13111892