Chaos & Predictability Notation How to run

Notation

The symbols and sign conventions here are mandatory across every chapter and every chaoslib docstring. A chapter that needs a new symbol adds it here in the same commit. Where the two course texts disagree, this book follows Kalnay (2003) for data assimilation and Palmer & Hagedorn (2006) for predictability.

Units are SI unless a standard domain alternative is noted.

Dynamical systems

SymbolMeaningNotes
xxstate vectorxRnx \in \mathbb{R}^n; nn is the state dimension
x˙=f(x)\dot x = f(x)the dynamicsautonomous unless stated
M\mathcal{M}the nonlinear model mapx(t+τ)=M(x(t))x(t+\tau) = \mathcal{M}(x(t))
tttimeMTU (model time units) for L63; time units for L96
τ\taulead time / optimisation windowalways positive
J\mathbf{J}Jacobian fi/xj\partial f_i/\partial x_jof the continuous dynamics
M\mathbf{M}tangent linear propagatorof the discrete map over τ\tau
M ⁣\mathbf{M}^{\!\top}the adjointtranspose under the Euclidean inner product
δx\delta xperturbation / errorδ0=δx(0)\delta_0 = \|\delta x(0)\|
λi\lambda_iLyapunov exponentsordered λ1λ2\lambda_1 \ge \lambda_2 \ge \dots, units of inverse time
λ(x,τ)\lambda(x,\tau)finite-time (local) exponenta property of a state, not the attractor
σi\sigma_isingular values of M\mathbf{M}ordered descending; σ1\sigma_1 is optimal growth
hKSh_{KS}Kolmogorov–Sinai entropyλi>0λi\sum_{\lambda_i>0}\lambda_i, nats per unit time
DKYD_{KY}Kaplan–Yorke dimensionfrom the spectrum
D2D_2correlation dimensionfrom a sampled trajectory
TdT_derror-doubling timeln2/λ1\ln 2/\lambda_1
Kτ\mathcal{K}^\tauKoopman operator(Kτg)(x)=g(Mτ(x))(\mathcal{K}^\tau g)(x) = g(\mathcal{M}^\tau(x)); linear on observables
ggan observable, or a dictionary of thema function of the state, not a state
K\mathbf{K}the finite-dictionary compression of Kτ\mathcal{K}^\tauwhat EDMD fits
ϕ\phia Koopman eigenfunctionKτϕ=eλτϕ\mathcal{K}^\tau\phi = e^{\lambda\tau}\phi
ν\nuinvariant (natural) measure“climate” as a distribution on state space
Aν\langle A\rangle_\nuensemble averageAdν\int A\,d\nu; the object a climatology estimates
AˉT\bar A_Ttime average over a window TTthe estimator a control run provides
τint\tau_{\rm int}integrated autocorrelation timeρ(s)ds\int\rho(s)\,ds; a time, not a lag count
τerg\tau_{\rm erg}ergodic timehow long to sample the whole attractor, not just its fast part
W1W_11-Wasserstein distance between measurescarries the units of the variable
DND_N^*star discrepancy of a samplesampling error of the worst indicator observable

Time-unit convention. For Lorenz 63, 1 MTU is read as 5\approx 5 atmospheric days; for Lorenz 96 with F=8F=8, 1 time unit is read as 5 days. These are interpretive conventions, not derivations — state them as such.

Specific systems

SymbolMeaning
σ,ρ,β\sigma, \rho, \betaLorenz 63 parameters (Prandtl, Rayleigh, geometry)
ρH\rho_HHopf threshold, σ(σ+β+3)/(σβ1)\sigma(\sigma+\beta+3)/(\sigma-\beta-1)
C±C^\pmthe non-trivial Lorenz 63 fixed points
N,FN, FLorenz 96 site count and forcing
rrlogistic-map parameter
μ\mudouble-well tilt; the fold is at μc=2/33\mu_c = 2/3\sqrt3
θ,ω\theta, \omegapendulum angle and angular velocity
L,g,mL, g, mpendulum length, gravity, mass

Note the collisions, all of them deliberate rather than tidied away, because each symbol is standard in its own literature. K\mathbf{K} is the Kalman gain of chapters 18 to 20 and the finite-dictionary Koopman operator of chapter 31; the two never appear in the same chapter, and chapter 31 says so where it introduces it. gg is gravity for the pendulum and an observable for the Koopman operator. μ\mu is the double-well tilt of chapters 27 and 30 and is the conventional symbol for an invariant measure; this book writes the measure ν\nu throughout and reserves μ\mu for the tilt. And σ\sigma is the Lorenz 63 Prandtl number, the noise amplitude of a stochastic run, and the conventional symbol for a singular value. Chapters using both write singular values as σi\sigma_i with an explicit subscript, and say so in the text.

Predictability and error growth

SymbolMeaning
E(t)E(t)forecast error (RMS unless stated)
EE_\inftysaturation error — the climatological level
Δt\Delta tgain in forecast horizon
D(pq)D(p\|q)relative entropy (Kullback–Leibler), nats
I(X;Y)I(X;Y)mutual information, nats
H(p)H(p)Shannon entropy, nats

Information units are nats (natural logarithm) everywhere. Divide by ln2\ln 2 for bits; never mix the two in one figure.

Data assimilation (following Kalnay 2003)

SymbolMeaning
xbx^bbackground (prior) state
xax^aanalysis (posterior) state
xtx^ttrue state
yyobservations
H\mathbf{H}observation operator
B\mathbf{B}background error covariance (fixed, as in 3D-Var)
Pb,Pa\mathbf{P}^b, \mathbf{P}^abackground / analysis error covariance (evolving)
R\mathbf{R}observation error covariance
K\mathbf{K}Kalman gain
J(x)J(x)variational cost function
NNensemble size

Sign convention for innovations: d=yHxbd = y - \mathbf{H}x^b, observation minus background. The analysis increment is Kd\mathbf{K}d, so a positive innovation moves the analysis towards the observation.

Ensembles and verification

SymbolMeaning
xˉ\bar xensemble mean
ssensemble spread (RMS about the mean, N1N-1 normalisation)
CRPScontinuous ranked probability score
BSBrier score
ACCanomaly correlation coefficient

Calibration identity: a reliable ensemble satisfies RMS spread == RMS error of the ensemble mean. State any departure from it as under- or over-dispersion, not as “error”.