Holistic view of our universe — Expanding Solar System Resonance Theory

The orbits of the eight planets, the wobble of Earth’s axis, and the rhythm of ice ages are all governed by the same physics. This holistic view puts them in one place: standard Newtonian dynamics carrying every planet from a single cited J2000 state, and Earth’s long periods read as beats against its lunisolar precession clock — one integrated account of solar system dynamics and Earth’s climate record, in a form you can run.

Two Counter-Rotating Motions

The model starts from a single observation: two of Earth’s precession motions rotate in opposite directions.

Earth’s two counter-rotating precessions. The perihelion direction advances counter-clockwise around the Sun at T ≈ ~111,582 years (apsidal precession); Earth’s Wobble Center circles Earth clockwise at T ≈ ~25,771 years (axial precession, the clock). Their ratio reads 4.330 today — a reading that wanders, not a constant.

Because the two motions run in opposite directions their rates add, and their beat is the ~20,936-year perihelion-of-date cycle. From this starting point the model derives quantities normally calculated separately: precession of the equinoxes, obliquity oscillation, eccentricity cycles and Milankovitch beat frequencies, length of days and years, and climate patterns (ice ages) — while the eight planets themselves ride the model’s own N-body engine, derived from one cited J2000 state with zero fitted constants.

Everything is read against one clock: Earth’s other long periods are beats of the precession clock with its orbit’s own motions, and the planets’ periods are the N-body engine’s eigenfrequencies. A 3D simulation visualises it all in one interactive view.

Why This Matters

Modern celestial mechanics — JPL Horizons, Laskar’s secular theory, high-precision ephemerides — achieves sub-arcsecond accuracy on Myr timescales and remains the right instrument for high-precision position prediction. This framework claims no better accuracy and no new physics. What it adds is integration:

1. Each phenomenon normally has its own model. Precession, obliquity, eccentricity, day and year lengths, eclipse timings and the climate spectrum are all driven by the same gravity, but are conventionally modelled independently, with separate parameterizations that never have to agree with one another.
2. Here they share one causal chain. One cited J2000 state and one derived clock feed every published quantity. Nothing is tuned per phenomenon, so a change anywhere has to stay consistent everywhere — which is what makes the whole thing falsifiable at once rather than one claim at a time.
3. The chain can be perturbed end to end. Multiply Jupiter’s GM by 1.01 and re-integrate: every planet’s apsidal rate moves, Mercury’s by +1.5″/cy — exactly one percent of Jupiter’s first-order Laplace–Lagrange share of Mercury’s precession — and Saturn’s window rate deepens its retrograde phase. The answer is not looked up; it falls out of the same run.
4. No comparable unified 3D view exists. Existing tools handle orbital mechanics or Milankovitch cycles in isolation. None bring axial precession, obliquity oscillation, perihelion precession, eccentricity cycles, day and year length variation, eclipses and all eight planets’ orbits into a single interactive scene.

Observation First — A Bottom-Up Approach

Standard paleoclimate theory takes the climate record as its starting point. The 100-kyr, 41-kyr, and 23-kyr cycles observed in ice cores and ocean sediments are interpreted by selecting which combinations of planetary motions best fit each peak — eigenmode beats from classical Laplace–Lagrange theory. The math is rigorous, but the climate-to-planet attribution is fitted to the data. Several major puzzles remain unexplained: the 100,000-year problem (eccentricity forcing is too weak to drive ice ages directly), the absence of the theoretically dominant 400-kyr cycle in post-MPT climate records, and the Mid-Pleistocene Transition ~1 million years ago.

This holistic view reverses the direction. It starts from directly observed orbital motion — Earth’s measured precession periods, every planet’s measured orbital cycles — anchored against JPL Horizons / WebGeoCalc ephemeris data (1800–2100 AD, every target inside 0.012°) and historical astronomical records. From these observations the model’s own N-body engine yields the planets’ secular modes, and Earth’s spin clock is derived from its year laws. Only then are the beats of those quantities — the climate formula’s 28 orbital lines — overlaid onto the climate record. The line periods are fixed by the model’s dynamics, not fitted to climate — only the per-line amplitudes are fitted to LR04 / CENOGRID / EPICA / CenCO2PIP.

The framework makes a structural commitment the standard approach does not: the orbital layer carries only the engine’s own beat lines, a sharply falsifiable claim. A climate peak with no counterpart among those beats would refute the orbital layer.

The canonical Climate Formula on LR04 δ¹⁸O over the past 700,000 years. R² = 0.779 post-MPT — the 28 engine-derived orbital lines (with the 405-kyr carbon-thermostat family) and the step transitions capture most of the post-MPT climate signal.

From Software to Theory

This framework began as a software project — a 3D simulation that reproduces the positions of the Sun, Moon, and all seven planets, verified against JPL Horizons ephemeris data and historical astronomical observations. The same run simultaneously produces obliquity, perihelion precession, and the planets’ inclinations and eccentricities — quantities normally computed by separate models. Its inputs are kept in three ledgers rather than a single parameter count (Mathematical Foundations): the planetary side has no free parameters at all, and the fitted constants live only in the time-domain correction stack.

From this, a theoretical framework emerged:

• Observable patterns in obliquity, eccentricity, orbital longitudes, and planetary motions could be captured in closed-form formulas
• Precise measurements of day lengths and year lengths in the 3D model also reduced to formulas
• A clock emerged: Earth’s axial precession period, derived from the year laws, against which every other Earth period is a beat — and which lengthens through geological time on the measured tidal history
• The planets moved onto dynamics: an N-body engine seeded by one cited J2000 state supplies the secular modes; an earlier integer-label reading of the planets’ periods (the “Fibonacci Laws”, a Jupiter–Saturn “lock”, Saturn’s “anti-phase” balance role, a formation-epoch KAM story) was tested against the engine’s own eigenfrequencies, found to carry no information, and retired with its record (Earth’s Clock)

The simulation also serves as a test bed: existing results such as Souami & Souchay’s (2012) invariable plane definitions can be verified, and six testable predictions — each naming the observation that would falsify it — can be checked against future observations.

Headline Findings — Climate Cluster

The model’s climate work is grounded in four independent proxy records (LR04 / CENOGRID / EPICA / CenCO2PIP). Three results stand out:

• R² = 0.779 — the 28 engine-derived orbital lines (L1, with the 405-kyr carbon-thermostat family inside) plus the boundary-condition step transitions (L3) on post-MPT LR04; 0.730 on EPICA CO₂, 0.692 on CenCO2PIP 0–66 Ma. (The comb-era formula with 33 integer-labelled lines fitted 0.87; the shipped physical-line formula pays the honest price of carrying only the engine’s own beats.) These are in-window descriptions: fitted on one half of a regime window, the lines have no skill on the other half (T7, pre-registered); the model’s own orbital histories predict 0.357 of post-MPT LR04 out of sample — the prediction, published beside the description.
• No fourth layer — adding classical Berger 1978 insolation features (ε, e, e·sin ϖ, e·cos ϖ), with the model’s own orbit or with Laskar 2004 / 2010, on top of the formula yields no cross-window-stable gain in any regime. Classical insolation features alone explain R² = 0.2362 of post-MPT LR04 — 3× less than the orbital lines alone. The beat lines carry the insolation-relevant variance.
• Two tiers at deep time — the formula’s precession-band and obliquity lines ride Earth’s lengthening clock; its eccentricity-band lines are planetary beats and do not. Every newly dated Precambrian section tests the split (Expanding Resonance).

Read the Climate Summary → for the full synthesis statement and empirical case in one table.

An open framework, not a closed theory

Explore the Model

Start reading: How It Works →

Quick Facts

• The engine: Newtonian N-body from one cited J2000 state, zero fitted constants — it carries every planet’s orbit and Earth’s own eccentricity, inclination, node and perihelion of date (Moon & Planets)
• The clock — axial precession: 25,771.4 years at J2000 (derived; lengthens through geological time)
• Apsidal precession: ~111,582 years (4.330 clock periods — a reading that wanders, between 1.07 and 9.85)
• Perihelion precession (of date): ~20,936 years (0.8124 clock periods)
• Obliquity cycle: ~41,224 years (the beat of the clock against the orbit plane’s nodal turn)
• Obliquity range: ~22.21° – ~24.72°
• The fitted anchor: 335,317 years — the unit of the correction bases, 13.011 clock periods; bookkeeping, not a cycle