Plane Calibration

The model uses two distinct ascending node calculation systems. This page documents the refinement of planetary ascending nodes to the invariable plane — a calibration that achieves sub-arcsecond accuracy for all planets.

Background: Souami & Souchay (2012)

Souami & Souchay (2012) published the definitive modern determination of the solar system’s invariable plane, including ascending node longitudes for all planets. Their work provided:

• The invariable plane orientation at J2000: inclination 1°34’43.3” (~1.57869°) to the ecliptic, ascending node at 107°34’56” (~107.582°)
• Planetary inclinations to the invariable plane
• Ascending node positions on the invariable plane for each planet

These values form the foundation of the model’s invariable plane calculations. However, a systematic issue arises when combining S&S values with the model’s dynamic inclination framework.

The Problem: Epoch-Specific vs Mean Inclinations

The S&S ascending node values are calibrated to the J2000 epoch. The Holistic Universe Model uses mean inclinations (time-averaged over each planet’s full oscillation cycle) as the baseline for its dynamic inclination formula:

i(t) = i_mean + A × cos(ω̃_ICRF(t) - φ)

Earth’s inclination to the invariable plane illustrates the issue:

When using Earth’s mean inclination (1.48113°) in the geometric relationship between orbital planes, the original S&S ascending nodes produce ecliptic inclinations that deviate from observed JPL values. The ascending nodes need recalibration to work correctly within the mean-inclination framework.

The Solution: J2000-Verified Ascending Nodes

Geometric Principle

Four parameters suffice: any planet’s ecliptic inclination is fully determined by its invariable-plane elements (i_inv, Ω_inv) and Earth’s invariable-plane elements (i_Earth, Ω_Earth).

The angle between any two orbital planes (the ecliptic inclination) is determined by three quantities:

1. Each plane’s inclination to the invariable plane
2. The difference in their ascending node positions

This relationship follows directly from spherical trigonometry:

cos(i_ecl) = cos(i_planet) · cos(i_earth) + sin(i_planet) · sin(i_earth) · cos(ΔΩ)

Where:

• i_ecl = ecliptic inclination (the angle between the planet’s orbit and Earth’s orbit)
• i_planet = planet’s inclination to the invariable plane (from S&S 2012)
• i_earth = Earth’s inclination to the invariable plane (1.57869° at J2000)
• ΔΩ = difference in ascending node longitudes (planet minus Earth)

Since i_ecl is known from JPL observations, and i_planet and i_earth are known from S&S, the ascending node difference ΔΩ is the only unknown. Solving:

cos(ΔΩ) = [cos(i_ecl) - cos(i_planet) · cos(i_earth)] / [sin(i_planet) · sin(i_earth)]

Two Independent Verification Methods

The calibrated ascending nodes were computed using two completely independent methods:

Both methods produce identical results, confirming the geometric validity of the approach. The numerical method minimizes the error between calculated and target ecliptic inclinations using dot products of orbital plane normal vectors. The analytical method solves the spherical triangle directly.

Results: Calibrated Ascending Nodes

Comparison Table

*Pluto is a dwarf planet (IAU 2006), included here for completeness; it does not participate in the model’s eight-planet inclination/eccentricity balance laws.

Key Observations

• Adjustments range from -6.00° (Pluto) to +5.97° (Jupiter)
• Mercury and Uranus require the smallest corrections (< 1°)
• Jupiter requires the largest positive adjustment (+5.97°)
• All verified values reproduce JPL ecliptic inclinations to better than 0.0001°

Why This Is Not Circular Reasoning

A potential concern: are we simply fitting free parameters to match observations?

The answer is no, because the three inputs come from independent sources:

The ascending node is not a free parameter — it is geometrically determined by the other two independently measured quantities. Given the planet’s tilt to the invariable plane and its tilt to the ecliptic, only one (or two, with a sign ambiguity) ascending node position is mathematically possible.

The sign ambiguity (two arccos solutions) is resolved by choosing the solution closest to the S&S original value — a physically motivated choice since the S&S values are already close to correct.

Dynamic Inclination: Coupling ICRF Perihelion to Tilt

Each planet’s inclination on the invariable plane is directly coupled to its ICRF perihelion longitude (not its ascending node). The formula is:

i(t) = i_mean + A × cos(ω̃_ICRF(t) - φ)

Where:

• ω̃_ICRF(t) = ω̃_J2000 + (360° / T_ICRF) × (t − 2000) — the precessing ICRF perihelion longitude
• T_ICRF = ICRF perihelion period (ecliptic rate minus general precession H/13)
• A = inclination amplitude ψ / (d × √m) from Law 2
• φ = per-planet inclination cycle anchor (Earth ~21.77°, Saturn ~120.38°, varying per planet)

The ascending node also precesses on the invariable plane, but at a different rate than the perihelion. Its role in the model is to set the direction of each planet’s angular momentum perturbation vector — the quantity that determines the vector balance of the invariable plane. The ICRF perihelion longitude drives the inclination phase (the magnitude of the oscillation), while the ascending node sets the direction in which that magnitude acts. Both are needed: the inclination phase is what oscillates over time, and the ascending node direction is what makes the vectors cancel.

Per-Planet Cycle Anchors

Each planet has its own cycle anchor φ — the ICRF perihelion longitude at which the planet reaches its inclination extremum (MAX for in-phase planets, MIN for Saturn). These are set by the System Reset epoch (2,649,854 BC, n=7). For example, Earth’s cycle anchor is ~21.77° and Saturn’s is ~120.38°.

The model has two balance groups: seven planets are in-phase (cosine sign positive), while Saturn alone is anti-phase (cosine sign flipped). This 7+1 split is what enables the inclination balance (Law 3) — Saturn’s contribution must equal the sum of all seven other planets’ contributions. The full anti-phase alignment — Saturn at maximum while all others at minimum — occurs once per Solar System Resonance Cycle (8H).

The cycle anchor for each planet is derived directly from its own ICRF perihelion longitude at the System Reset, ensuring full consistency between the inclination oscillation and the underlying perihelion dynamics.

Complete Parameters

The Period column is the inclination oscillation period — equal to the planet’s ICRF perihelion period (T_ICRF), which is what the formula at the top of this section uses. Earth is the sole prograde-in-ICRF planet; all other planets are retrograde in ICRF (their ecliptic rates are slower than the general precession H/13).

Period Ratios

The inclination oscillation periods (ICRF frame) are derived from each planet’s ecliptic perihelion period via the subtractive Fibonacci relation (ICRF rate = ecliptic rate − general precession H/13). The result is that all ICRF periods are simple fractions of the 335,317-year Earth Fundamental Cycle (H):

Dual Coordinate Systems

The simulation tracks ascending nodes in two coordinate systems because height calculations require ecliptic-frame coordinates while long-term dynamics use inertial coordinates:

*Periods shown are Earth’s (H/3 ICRF, H/16 ecliptic). Other planets have their own per-planet ICRF periods (see Complete Parameters table above). Earth itself is a special case for the ecliptic ascending node: it precesses at −H/5 (~67,063 years, ecliptic precession) rather than H/16 — see the per-planet residuals section below.

The relationship between the two rates follows from the precession meeting frequency:

1/P_ecliptic = 1/P_ICRF + 1/P_axial
1/~20,957 ≈ 1/~111,772 + 1/~25,794

This is the same relationship that produces the perihelion precession period from the axial and inclination precession periods — another manifestation of the model’s unified framework.

Ecliptic Ascending Node System

Separate from the invariable plane system, the model also tracks how ascending nodes shift on the ecliptic as Earth’s axial tilt (obliquity) changes. This uses a different formula:

dΩ/dε = -sin(Ω) / tan(i)

Where:

• Ω = ascending node longitude on the ecliptic
• ε = obliquity (axial tilt)
• i = orbital inclination to the ecliptic

As the obliquity oscillates over the ~41,915-year cycle, the ecliptic itself tilts slightly, shifting where planetary orbits cross it. This is an ecliptic-frame effect that operates independently of the invariable plane node precession.

Observability of Ascending Node Periods

The model assigns each fitted planet’s ascending-node precession period as 8H/N — a simple integer divisor (N) of the Solar System Resonance Cycle (8H = ~2.68 million years). These periods describe motion over 50,000–2,000,000 year timescales (Laskar’s measured eigenfrequencies fall in the same range). With only ~4,000 years of recorded astronomical observations at most, humanity cannot directly observe a complete cycle for any planet.

This creates an important observability limitation: short-term measurements (≤200 years) cannot distinguish between:

Across all 7 fitted planets, the 8H/N integers reproduce JPL’s J2000-fixed-frame ascending-node trends with a cumulative residual error of ~5.8″/century (≈0.8″/century per planet). Over observable timescales, different integer choices produce indistinguishable motion.

Per-planet residuals

The cumulative ~5.8″/century error is the sum of absolute per-planet residuals between the model’s 8H/N trend and the JPL-reported trend. Earth’s ascending node is set by formula at −H/5 (the ecliptic precession period of ~67,063 years, not a fitted parameter) and is therefore excluded from the comparison.

Most residuals are well below 1″/century. Venus (3.2″/cy) and Saturn (1.7″/cy) dominate the cumulative error; together they account for ~5 of the ~5.8″/century total. The remaining five planets each contribute < 0.5″/century. All 7 fitted planets match the JPL-reported direction of motion (Dir. ✓).

The scalar balance is what’s verifiable: Laws 3 and 5 (inclination + eccentricity balance) operate on the J2000 snapshot — these are testable with current data. The ascending node dynamics, like Earth’s eccentricity cycle predictions, are theoretical extrapolations that will only be observationally testable on much longer timescales.

Verification and Reproducibility

All calibration work is fully reproducible through publicly available scripts:

Summary

Return to Invariable Plane or explore Formulas for the complete formula set.