The Invariable Plane

The invariable plane is the solar system’s true reference plane — the plane perpendicular to the total angular momentum vector, fixed in space while everything else moves around it. Imagine a spinning top: no matter how it wobbles, the spin axis stays constant.

The model’s ~111,772-year inclination precession is Earth’s orbital plane tilting relative to this plane (see How It Works). Without a fixed reference that motion cannot be defined. Beyond providing a reference, the invariable plane reveals structure: planetary inclinations measured against the (moving) ecliptic look irregular, but against the invariable plane the same data shows smooth, predictable oscillations and a clean Fibonacci structure.

Planetary inclinations

Every planet’s orbit is tilted relative to the invariable plane. The tilts oscillate in predictable patterns around mean values.

*Venus and Saturn have retrograde apsidal precession (opposite direction to the other planets).

†Periods shown are ecliptic-frame perihelion periods, which is what direct observation reports. Earth is the special case: its row uses the ICRF perihelion period (H/3 = ~111,772 yr) because Earth is the only planet whose ICRF perihelion is prograde — see Why Earth Is Unique. The inclination-oscillation formula below uses ICRF perihelion for all planets.

Notable patterns

• Earth and Uranus both have an H/3 period (~111,772 yr) — Earth’s is the ICRF perihelion / inclination-oscillation period; Uranus’s is the ecliptic perihelion period. A numerical coincidence, not a shared physical period.
• Jupiter’s ecliptic perihelion period (68,783 yr = 8H/39) sits close to Earth’s ecliptic precession period — suggesting Jupiter’s gravitational torque drives the precession of Earth’s orbital plane.
• Saturn’s ecliptic perihelion period (41,270 yr = 8H/65) coincides exactly with Jupiter’s ICRF perihelion period at 8H/65 — the gas-giant lock that drives Earth’s obliquity (k+s₃ climate beat). Earth’s own obliquity Fibonacci anchor is H/8 = 8H/64, one lattice integer off. Canonical: Fibonacci Laws §Law 6; see also Predictions §17 — Invariable Plane Tilt.
• Saturn is the only planet whose longitude of perihelion moves obviously retrograde in the ecliptic frame at the current epoch (~-3,400 arcsec/cy from JPL WebGeoCalc). The model interprets this as a permanent feature (8H/65 = 41,270 yr, ecliptic-retrograde), not a transient oscillation — full analysis at Supporting Evidence §12.
• All periods are simple fractions or multiples of the 335,317-year Earth Fundamental Cycle (J2000 anchor; the H/N and 8H/N integer-divisor labels are invariant, but the literal year counts rescale at geological time — see Expanding Resonance).

Inclination oscillation

Each planet’s inclination to the invariable plane oscillates around its mean value over its precession period.

Gravitational pull from the other planets causes each orbit to precess around the invariable plane, producing two coupled effects: nodal precession (the ascending node rotates around the plane) and inclination oscillation (the tilt angle oscillates toward and away from the plane). The general formula:

inclination(t) = mean + amplitude × cos(ϖ_ICRF(t) − phase_angle)

ϖ_ICRF(t) is the planet’s ICRF perihelion longitude (ecliptic rate minus general precession H/13). The phase_angle is the planet’s own ICRF perihelion longitude at the Balanced Year (302,635 BC), which determines whether the planet is in-phase (cosine sign positive) or anti-phase (sign flipped). The ascending node precesses on the invariable plane at a different rate; its role is to set the direction of the angular momentum perturbation vector that determines the vector balance below. The J2000-verified ascending node values — refined from Souami & Souchay (2012) using spherical trigonometry — match JPL ecliptic inclinations to <0.0001°.

Earth’s oscillation

For how Earth’s inclination drives obliquity, see Obliquity & Inclination.

Why not use the ecliptic?

The ecliptic (Earth’s orbital plane) is commonly used as a reference, but it moves:

Measured against the ecliptic, planetary inclinations look chaotic because the reference plane itself oscillates. The clean structure visible against the invariable plane disappears:

A planet’s ecliptic inclination depends on both the planet’s tilt to the invariable plane and Earth’s tilt to the invariable plane — both oscillating at different rates, hence the complex appearance. The J2000 values below are the JPL-observed targets used in the ascending node calibration:

Visualising the invariable plane

In the Interactive 3D Simulation:

1. The invariable plane is shown as a reference grid (under celestial tools).
2. Earth’s orbit is tilted ~1.57869° to this plane.
3. Earth moves above and below the plane during its yearly orbit.

Earth crosses the plane twice per year. The crossing dates shift slowly over millennia as inclination and node precess. Around J2000:

Maximum distance above or below is small — about 4 million km at the extremes.

Why the plane is stable: vector balance

The invariable plane is stable because all planetary angular momentum perturbations cancel as 2D vectors. Each planet’s orbital tilt creates a “force” on the plane with a direction (ascending node Ω) and magnitude (L × sin i). For the plane to be stable, all 8 forces must sum to zero. The force is highly concentrated in the outer planets:

The 42/58 split of L × sin(i_J2000) magnitudes is not a balance by itself — the actual Law 3 scalar balance uses a different weight formula (w_j = √(m·a(1-e²))/d × amp) and produces the 99.9974% balance. The model achieves 99.9974% scalar balance — the structural weights of the in-phase group closely match those of the anti-phase group. This is a genuine physical constraint that selects the Fibonacci d-values. The vector balance (whether the 2D forces cancel at all times) depends on how the ascending nodes precess: in the model’s single-mode approximation, vector balance degrades to ~72% over time as nodes go out of sync; in a multi-mode representation (7 eigenfrequencies per planet, matching secular perturbation theory), vector balance is maintained at 100% — but this property holds for any frequencies, not just specific ones.

Ascending node periods: a structural prediction

Each planet’s ascending node regression period takes the form 8H/N for integer N, with Jupiter and Saturn locked at a shared N=36. Across all 7 fitted planets, the 8H/N integers reproduce JPL’s J2000-fixed-frame ascending-node trends with cumulative residual ~5.8″/century (~0.8″/century per planet).

The model derives all eight from a single constant H. Laskar’s secular theory measures them as eight independent eigenfrequencies with no known structural relationship to each other.

Compute inclination at any year

The full closed-form expressions for planetary inclinations to the invariable plane at arbitrary year are in Formulas. For the underlying ascending node values and how they were derived: Plane Calibration.

Continue to Moon & Planets for how the model handles lunar and planetary motion, or Plane Calibration for the technical methodology behind the refined ascending node values used throughout this page.