MPE StudioMath of Planet Earth
Lab · Exploration 1

The Butterfly and the Forecast

Deterministic does not mean predictable.

Explore how rounded initial conditions, unstable dynamics, and the prediction target determine how long a detailed forecast remains useful.

This Exploration develops the complete argument while letting you change the assumptions and watch what survives.
full-precision restartrounded restartnearly the same initial statedetails separate
A tiny initial discrepancy can remain hidden for a time, then grow into a completely different detailed forecast.
Playing
Thumbnail for Why Prediction Has Limits (Even with Perfect Models)

Watch first · short video

Why Prediction Has Limits (Even with Perfect Models)

Charlotte Moser introduces why deterministic dynamics can still lose predictability.

Watch video ↗Applied Mathematics in Geosciences · Episode 4
Related: From Chaos to Probability ↗
Then explore it yourself ↓
01

The restart that did not repeat

An exact rule is not the same as an exactly known starting state. Round the state, restart the same equations, and follow what happens.

State held by the computer0.506127
printed to 3 decimal places
State printed on the page0.506
restart from the printed value
Illustrative restart vector(0.506, -5.800, 24.200)

The first number echoes the reported rounding example. The three-component Lorenz-63 state is a modern illustration, not Lorenz's historical state.

0510152025time after restartstate variable x(t)trajectories shadow each othertrajectory details divergeoriginalrestart
The two x trajectories begin close together. At the displayed time their full-state separation is 4.62e-4.
Initial separation
5.380e-4
Current separation
4.619e-4
First tolerance crossing
4.30 time units

Same equations. Nearly the same initial state. Detailed trajectory information has a finite lifetime.

02

The direction of an error matters

A small initial error does not only have a size—it also has a direction. In an unstable system, perturbations pointing in different directions can grow at very different rates.

Uncertainty has a geometry

δx(t + Δt) ≈ Aδx(t)

The matrix A is a local approximation of the dynamics. It tells us what happens to a very small perturbation over a short interval: the perturbation may rotate, stretch, or compress.

Both perturbations begin with exactly the same size. One points along a direction strongly amplified by the dynamics; the other points along a direction that is temporarily compressed.

1 · initial perturbation δx(t)2 · local dynamics A3 · evolved perturbation δx(t + Δt)same initial sizedifferent amplificationrotated, stretched, and compressed
Both initial vectors have length 1. Vector A is amplified by 2.39, while vector B is amplified by 0.60.
Warm orange vector
2.39× amplification
Teal vector
0.60× amplification
Show the mathematics

δx(t + Δt) ≈ Aδx(t)

A=2.090.690.690.86

The matrix A is a local approximation of the dynamics. It tells us what happens to a very small perturbation over a short interval: the perturbation may rotate, stretch, or compress.

Knowing only the size of the initial error is therefore not enough. Its direction in state space matters too. Forecast uncertainty can grow much faster in some directions than in others.

03

Trajectory versus structure

Nearby paths can lose detailed agreement while remaining organized by the same dynamical structure.

state variable xstate variable znearby starting statesdifferent pathssame bounded attractor
The trajectories now differ by 0.000 in full-state distance. Both are produced by the same deterministic equations.

Detailed trajectory agreement can be lost while the structure organizing the motion persists.

04

How much forecast time does precision buy?

A forecast horizon depends on the initial error, how errors grow, the prediction target, and the tolerance chosen for that target.

error(t) ≈ initial error × exp(λt)

Tp ≈ (1/λ) log(tolerance / initial error)

10-810-610-410-2100lead timeforecast error (log scale)chosen toleranceTₚ = 13.8
The selected forecast horizon is 13.8 model time units.
Predictability horizon
13.8
Horizon gained vs. default
0.0 time units
Initial-error reduction
1×

A 1-fold reduction in initial error changes forecast time logarithmically. For 100×, ΔTₚ = (1/λ) log(100), not 100× the time.

A predictability horizon is chosen relative to a target and tolerance. It is not a universal moment when all prediction becomes impossible.

05

From one forecast to a family of futures

An ensemble begins from a family of plausible states. What it represents depends on which uncertainties the experiment includes.

lead timeforecast variablecompact initial uncertaintygrowing family of futures
18 members sample initial-state uncertainty. The darker line, when shown, is the ensemble mean, not truth.

More members sample the uncertainties already included here. They do not automatically repair a missing source of uncertainty.

Ensemble spread becomes a probability statement only after the forecasting procedure has been evaluated against outcomes.

An ensemble is a scientific object, not a decorative fan of lines.

06

Chaos is not randomness

These are two different kinds of perfectly valid models. The difference is not that one model is correct and the other is noisy or incorrect. The difference is where uncertainty enters the dynamics.

Deterministic model

Same model + same initial state → same trajectory

evolution time

Run 1: exact overlap with the reference trajectory.

A deterministic model contains no random input. Once the equations and the initial condition are fixed, the future trajectory is fixed. Even a chaotic deterministic model behaves this way. Chaos means that slightly different initial conditions can eventually produce very different trajectories—not that repeated runs from the identical initial condition are random.

Stochastic model

Same model + same initial state → different realizations

evolution time

previous realizationcurrent realization

A stochastic model contains a random term as part of the model itself. Therefore, even when every run begins from exactly the same initial state, different realizations can evolve differently because they experience different random inputs.

Deterministicdx/dt = f(x)

Stochasticdx = f(x)dt + σdWt

The term dWt represents new random forcing entering during the evolution.

Deterministic chaos

  • same initial state → same trajectory
  • nearby initial states → trajectories may diverge
  • uncertainty is amplified by the dynamics

Stochastic evolution

  • same initial state → different realizations are possible
  • random forcing enters during the evolution
  • uncertainty can continually be introduced

Chaos is sensitivity, not randomness. A deterministic chaotic system amplifies differences that are already present. A stochastic system can introduce new differences as it evolves.

Stochastic terms can represent unresolved processes statistically. Their use does not imply that the underlying physical world is literally random.

Chapter takeaways
01

Determinism is not the same as predictability.

An exact rule determines a unique trajectory only after a starting point is specified; a real forecast begins from an estimate rather than the exact state.

02

Uncertainty has its own dynamics and geometry.

It can be stretched, rotated, and compressed, and its most rapidly growing direction can change with the state.

03

Losing one trajectory does not mean losing all predictability.

Distributions, regime probabilities, averages, and forced responses may remain informative, while slowly evolving variables can retain useful memory.

04

An ensemble is a scientific object, not a decorative fan of lines.

It is meaningful only when it represents the uncertainties that matter and is evaluated against outcomes.

What the Lorenz model leaves out

Lorenz-63 omits spatial structure, moisture, external forcing, ocean coupling, unresolved scales, observational networks, and model error.

  • Exponential growth is a local approximation.
  • Real errors may contract temporarily, rotate, interact across scales, and saturate.
  • Chaos is only one source of forecast failure.
  • Higher numerical resolution does not automatically make a conclusion more credible.
Precision of calculation is not the same as certainty of conclusion.
Sources and further reading
  • Edward N. Lorenz, Deterministic Nonperiodic Flow
  • Edward N. Lorenz, The Essence of Chaos
  • Operational ensemble-forecast documentation from major weather centers