Expanding Solar System Resonance Theory (ESSRT) — The Lattice Through Deep Time

This page is the synthesis of the time dimension. For the mechanism that holds the lattice stable, see Physical Origin; for one application of the lattice to Earth’s climate, see Climate Summary; for the paleontological day-count tables that validate the time-evolution, see Supporting Evidence; for falsifiable predictions about deep past and far future, see Predictions.

1. What expands, and what stays the same

The Solar System Resonance Cycle 8H and its integer divisors form a lattice. The lattice has two parts:

• Integer labels — the integers N in 8H/N and H/N. These are structural constants: n = 65 for the obliquity main beat, n = 39 for Jupiter’s ecliptic perihelion, n = 16 for the perihelion harmonic, and so on across all the predicted periods in the model.
• Absolute periods — the duration of each cycle in years. These scale with H(t).

When H expands, every H/N and 8H/N expands by the same proportional factor. The lattice’s shape — the integer ratios between cycles — is preserved exactly. What changes is the time unit in which the lattice is measured.

Every row of this table moves by the same fractional amount per epoch. The lattice is rigid in its proportions; flexible in its scale. A useful intuition for the rate: one billion years ago, H was about ~80 % of its current value. The growth is slow on geological timescales but cumulative over the planet’s lifetime.

2. Two drivers, two physical processes

The expansion is driven by two independent physical processes that act on different parts of the system. They produce complementary effects on the lattice.

Driver 1 — Earth-Moon tidal evolution

The Moon raises a tidal bulge on Earth. Earth’s rotation drags that bulge slightly ahead of the Moon-Earth line. The misalignment creates a torque that slows Earth’s rotation and transfers angular momentum to the Moon’s orbit, pushing the Moon farther out. Modern lunar laser ranging measures the recession at 3.83 cm/yr; the long-term Phanerozoic average is ~3.43 cm/yr.

As Earth’s rotation slows, the day grows longer (LOD increases). Through the model’s structural relation H = 13 × axial precession period, a longer day means a longer H. So Driver 1 expands the temporal lattice:

Moon recedes → LOD grows → axial precession period grows → H grows → 8H grows → every H/N grows

This is the dominant driver in fractional terms. At Devonian, H is ~7.8 % smaller than today; the change comes entirely from Driver 1.

Driver 2 — Solar mass loss

The Sun loses mass through electromagnetic radiation (mass-energy via L_☉ / c²) and the solar wind, at a combined rate of about 9.3 × 10⁻¹⁴ of its mass per year. For each planet in orbit around a slowly-shrinking central mass, the adiabatic invariant a × M_☉ ≈ constant predicts the orbit expands as M_☉ decreases. Going to past, the Sun was more massive, so every planet’s orbit was smaller.

Every planet drifts by the same fractional amount. At Devonian, the entire solar system was ~35 ppm tighter — Mercury was at 57,907,130 km (now 57,909,176), Jupiter at 778,519,688 km (now 778,547,200), Neptune at 4,503,284,523 km (now 4,503,443,661); every distance moved by the same proportional amount. At Hadean, ~423 ppm tighter. Driver 2 rescales the spatial lattice uniformly. Planet orbit ratios and the L1 lattice’s integer structure are completely preserved across all epochs.

How the two drivers interact

The drivers act on different physical channels — tidal coupling for Earth-Moon angular momentum; gravitational binding for Sun-planet orbits. But the model’s structural near-invariant ties them at the observational level: the product H × (days/year) stays close to 122,471,920 days across geological time. This identity is exact at J2000 (the anchor). At deep time it drifts smoothly — about −71 ppm at Devonian, −845 ppm at Hadean — because Driver 2 shortens the year in seconds (the Sun was more massive in the past, so by Kepler’s third law the year was shorter). Driver 1 is ~100× larger in fractional terms than Driver 2, so for most paleoclimate purposes Driver 1 dominates the lattice expansion and Driver 2 contributes a small Gyr-scale correction.

3. The Hadean origin — and the model’s strongest self-validation

If H expanded from past to present, it must have been smaller in the past — perhaps much smaller. How small? Run the proper-physics formula backwards in time.

At Patterson’s Pb-Pb Earth age of 4.54 Gyr, the model places:

• Moon distance: 20,532 km = 3.22 Earth radii — just outside the Roche limit (~2.9 Earth radii, where a fluid Moon would tidally disintegrate)
• Length of day: ~5.0 hours
• H = ~69,837 yr (about 21 % of today)
• 8H ≈ 0.56 Myr (the smallest the cycle has ever been)

This is the model’s strongest self-validation. No Hadean constraint was used to calibrate the formula — the only anchors are modern LOD (24 hr at J2000), Wells’s Phanerozoic tidal rate (0.00526 hr/Ma), and Farhat et al. 2022’s deep-time tidal-evolution curve calibrated on independent geological data. The formula puts the Moon at the Roche limit at exactly Patterson’s radiometric Earth age. Two independent measurement chains — paleoclimate / paleontology on one side, Pb-Pb radiometry on the other — converge on the same answer.

The non-coincidence

A second, independent self-validation falls out of the same numbers. Wells’s tidal rate of 0.00526 hr/Ma was calibrated entirely on Phanerozoic fossil data (0–500 Ma — corals, bivalves, rhythmites). Linearly extrapolated, it predicts LOD = 0 at 24 hr ÷ 0.00526 hr/Ma ≈ 4.56 Gyr ago. Patterson’s Pb-Pb radiometric Earth age is 4.54 Gyr. The two numbers agree to within 0.5 % — the precision of the underlying measurements. This is not a coincidence: it is the structural signature of a single, internally consistent picture of Earth-Moon evolution. The 8H cycle has a beginning, and the beginning is when the Moon’s orbit was just stable enough to exist as a satellite.

4. The proper-physics formula

The model’s LOD(t) and H(t) are computed by a two-layer formula:

Layer 1 — Moon-distance polynomial, calibrated against Farhat et al. 2022:

a_Moon(t) = a_Moon_now × (1 + α₁·t_Ma + α₃·t_Ma³ + α₄·t_Ma⁴)

The polynomial captures the time evolution of the Moon’s semi-major axis from t = 0 (today) backwards through the entire 4.54-Gyr deep-time record. The cubic and quartic terms model the non-linear tidal-dissipation history without requiring an explicit Q(t) integration.

Layer 2 — angular-momentum conservation (Earth-Moon system):

LOD(t) = 2π · I_E / (L_total − M_M · √(GM_EM · a_Moon(t)) · √(1 − e²))

Given Moon distance, LOD follows from conservation of the Earth-Moon total angular momentum L_total. From LOD, every other deep-time quantity follows: H(t), the lattice periods H/N and 8H/N, the Earth-Moon distance, the tidal-lock asymptote.

Verification at modern, Devonian, and Hadean is reproducible via scripts/devonian_cross_check.py in the source repository. The formula matches:

• Modern J2000 LOD = 24.000 hr (anchor, exact)
• Devonian (380 Ma) LOD = 22.12 hr → days/year = 396.21 vs Wells 1963’s coral count of ~400 (1 % match)
• Hadean (4.54 Gyr) Moon at 20,532 km = 3.22 Earth radii (Roche-limit self-validation, ~10 % outside Roche)

5. The structural near-invariant

A central identity ties the two drivers together: the product H × (days/year) is approximately 122,471,920 days. This value is the J2000 anchor — exact today by construction — and drifts smoothly with geological time (−71 ppm at Devonian, −845 ppm at Hadean). The physical interpretation: Earth rotates approximately the same number of times — about 122 million rotations — during each H cycle, across the planet’s entire history. The near-invariance arises from a cancellation: Driver 1 makes H longer in years while shortening days/year (LOD grows); Driver 2 shortens the tropical year in seconds (the Sun was more massive in the past). The two scalings almost cancel, leaving the rotation count per H cycle near-invariant but not exactly constant.

At J2000 the anchor value computes as 335,317 years × 365.2422036 days/year ≈ 122,471,920 days. Deep-time drift comes from Driver 2.

The Phanerozoic drift (< 100 ppm) is well within the precision of paleontological day-count measurements — Wells’s coral rings have ±1–2 % uncertainty per epoch. The Gyr-scale drift is real and explicitly modelled.

6. The Lunar Precession Invariant

A second deep-time invariant falls out of the framework — this one governing the Moon’s apsidal and nodal precession. Where the day-count near-invariant (§5) ties Earth’s rotation count to one H cycle, the Lunar Precession Invariant ties the Moon’s perigee and node advance to the same lattice:

T_apsidal × H = constant   (perigee advance, ICRF frame, in years)
T_nodal   × H = constant   (node regression, ICRF frame, in years)

Equivalently, the count of lunar apsidal/nodal cycles per H scales as (H(t)/H₀)². As the Moon recedes and H grows, the Moon’s precession periods grow in lock-step, so their product with H stays fixed at every epoch.

J2000 anchors and derived invariant value:

The invariant value is then derived from these anchors:

T_apsidal × H = H₀² / N_apsidal = 335,317² / 37,900 = 2,966,688 yr²
T_nodal   × H = H₀² / N_nodal   = 335,317² / 18,015 = 6,241,326 yr²

These numbers are empirically anchored (one structural H₀, one observed lunar period) — they are not Fibonacci-derived integers like 8H = 2,682,536 yr. What is structural is the claim that this value stays preserved across every epoch.

Where it comes from

Brown’s lunar perturbation theory gives the apsidal and nodal rates as proportional to m² at leading order, where m = n_Sun / n_Moon ≈ 1/13.37 is the ratio of solar to lunar mean motion. In period form (Brouwer-Clemence):

T_apsidal ∝ T_year² / T_sidereal_month     (Brown m² scaling, leading order)

This sets the scaling form. The leading m² term alone gives T_apsidal ≈ 17.8 yr — roughly double the observed 8.85 yr. This is the historical Newton-Clairaut problem: Newton’s leading m² calculation in the Principia (1687) gave half the observed perigee rate, and Clairaut showed in 1749 that the m³ and higher terms approximately double it to recover the observed value. The model therefore anchors the J2000 magnitude from observation and adopts the H² scaling form (which matches Brown m² leading order) to propagate to deep time. No polynomial corrections — the scaling is purely structural on the H-lattice.

Aligning, not replacing. The Lunar Precession Invariant does not compete with Brown’s lunar theory but sharpens it: the m²-leading-order scaling is adopted as a structural law, with the J2000 magnitude anchored from observation and propagated to deep time without polynomial corrections. Where Brown’s expansion derives the absolute period from the underlying m-series, ESSRT treats the period as one anchored input and the (H/H₀)² evolution as the structural claim — the two views agree on the scaling form by construction.

Under Driver 1, the Moon recedes → its sidereal month grows; angular-momentum conservation simultaneously slows Earth’s spin → LOD grows → H grows (via the structural relation H = 13 × axial precession). The two periods evolve together, and the framework’s choice N_apsidal(t) = N₀ × (H(t)/H₀)² makes T_apsidal × H an exact constant at every epoch.

A stricter invariant than the day-count

The Lunar Precession Invariant is the structurally cleanest of the model’s deep-time invariants. Where the day-count H × days/year (§5) drifts at deep time because the year in seconds slowly shortens under Driver 2 (~−850 ppm at Hadean), the Lunar Precession Invariant — expressed in year-units rather than seconds — does not involve the year-length and is therefore preserved by Driver 1 alone, exact at every epoch:

The T_apsidal × H column is the J2000-anchored value 2,966,688 yr² (exact: H₀² / N_apsidal,J2000 = 2,966,688.40 yr²), held constant by the framework’s (H/H₀)² scaling — not a structurally-derived integer. Both H (integer) and T_apsidal (6 decimals) are display-rounded, so the hand-reproduction T × H ≈ 2,966,688 is accurate to within ±10 yr². Same pattern for nodal: T_nodal × H = 6,241,326 yr² at every age.

A modern observer sees the Moon’s perigee advance once every ~8.85 yr; a Devonian observer would have seen it advance once every ~9.60 yr (slower, because H was smaller and the Moon was closer). But both observers count N₀ × (H/H₀)² precession cycles per H cycle — the same structural ratio at every epoch.

The deep-time invariant family

The model’s deep-time invariants form a small family — relations preserved across all epochs under the drivers, distinct from the J2000-static Fibonacci Laws. The Lunar Precession Invariant is the structurally cleanest member:

The Fibonacci Laws describe how the lattice is organized at J2000; the deep-time invariants describe how that organization survives as H(t) expands.

7. Validation against the geological record

The model’s deep-time predictions match the empirical paleo-LOD record across multiple independent measurement techniques.

The Phanerozoic record (0–380 Ma) is matched to within 1.2 %. The 620 Ma point (Williams 2000) sits in a known transition interval — the late Cryogenian Snowball Earth boundary — where the calibrated tidal-evolution curve passes between Williams’s direct rhythmite count and the modern Phanerozoic rate. See Supporting Evidence for the full validation discussion.

Wu et al. 2024 — independent 650-Myr cyclostratigraphy reconstruction

The seven Wu et al. anchors deserve special mention because they constitute an entirely independent measurement chain. Wu, Y., Malinverno, A., Meyers, S. R. & Hinnov, L. A. (2024) — “A 650-Myr history of Earth’s axial precession frequency and the evolution of the Earth-Moon system derived from cyclostratigraphy”, Science Advances, doi:10.1126/sciadv.ado2412  — apply a TimeOptB Bayesian inversion to cyclostratigraphic records (Milankovitch-driven sedimentary cycles) across 650 Myr of the Phanerozoic. Unlike paleontology (Wells 1963 coral growth bands, Pannella 1972 bivalves, Williams 2000 tidal rhythmites), Wu’s method does not depend on biological growth-band counting — it derives LOD, Earth-Moon distance, and axial precession rate directly from spectral analysis of climate proxies, with full ±2σ uncertainty propagation at 100-Myr intervals.

The model matches Wu’s 7-anchor reconstruction within −1.70 % to +3.11 % across the entire 650-Myr range, with a 7-point mean residual under 1.5 %. This is independent validation across an entirely separate measurement technique — paleontology and cyclostratigraphy converging on the same H(t) evolution. The match at 400 Ma is particularly significant: Wu’s cyclostratigraphy-derived 403.03 days/year sits within 0.7 % of Wells’s 1963 paleontological count of ~400 days/year, and the model reproduces both within 1.3 %.

Historical eclipses — confirmation in the past 2,400 years

The proper-physics formula also predicts a closed-form ΔT — the deviation between Terrestrial Time and Earth-rotation-tracking Universal Time. Two complementary tests against the documented historical eclipse record confirm the formula:

• Solar visibility test (19 events, -762 to 1654 CE): the model — pure-tidal Farhat + viscoelastic α(t) GIA correction — explains 19/19 events as visible at the documented observation site vs 17/19 for the empirical fit of Stephenson-Morrison. The underlying Moon polynomial agrees with NASA’s Five Millennium Catalog within ±15 minutes back to 2,500 years before J2000 (n = 11 canonical events, mean residual 6.9 min). See Solar Eclipse Validation.
• Lunar timing test (270 primary-source observations from Stephenson 2016, -720 BCE to 1280 CE): the higher-resolution test (Babylonian, Greek, Chinese, Arab traditions) reaches a mean residual of 26.7 min vs NASA’s 20.0 min under the L1-orbital-coupled α(t) refinement — a 6.7-min gap on top of the ~20-min per-observation noise floor. Four independent observation traditions agree on the model’s residual magnitude to within ±200 s after detrending. An independent solar cross-validation on 89 events (L-7) confirms the same physics. See Lunar Eclipse Validation.

The two tests together resolve a long-standing ambiguity in the eclipse-validation literature. The full Munk-MacDonald-scale (~5-6 ms/century) non-tidal Earth-rotation speedup postulate is rejected by the historical record: the model fits without it. A dominant GIA-scale (~0.6 ms/century) non-tidal channel is included via the α(t) viscoelastic correction — measured independently by satellite gravimetry (Cox & Chao 2002 + Peltier ICE-5G(VM2) multi-mode rheology), with zero parameters fitted to eclipse data. A smaller fractional non-tidal secular rate (~0.5 ms/century) is additionally detected in the lunar-timing residual but not currently modelled. The medieval residual (~1,000 s peak in the 840–1020 CE window; exact peak year is reference-conditional) decomposes structurally into a framework-native 8H/1851 lattice harmonic (= 73 × Jupiter-Saturn synodic; structural prediction, amplitude calibration deferred), the fractional non-tidal drift, and observation noise — see Lunar Eclipse Validation for the full decomposition and hypothesis-testing details.

8. The Expanding-Universe parallel

The two frameworks share a logical shape: a single scale parameter grows monotonically while the discrete structure it carries — integer ratios, dimensionless invariants — survives the change. Setting them side-by-side helps fix what does and does not vary across geological time. The parallel is pedagogical, not physical — cosmic expansion is metric and universal, while the expanding resonance is Newtonian and bounded to the solar system — but the shape of the claim is familiar:

Cosmic expansion is driven by a fundamental property of spacetime and dominates at the largest scales; the expanding resonance is driven by Newtonian mechanics plus the Sun’s mass-loss rate and acts on geological-to-astronomical timescales.

9. The asymptotic future — when does the cycle end?

If the cycle had a beginning at the Roche limit, it must have an end at the opposite extreme. The Earth-Moon system asymptotically approaches tidal lock: Earth’s spin period equals the Moon’s orbital period, both at about 47 days in current units. From angular-momentum conservation:

• Tidal-lock Moon distance: 555,623 km = 87.1 Earth radii (currently 60.3 Earth radii)
• Approach timescale: ~50 Gyr (well beyond the Sun’s red-giant phase at +5 Gyr)
• Days per year at tidal lock: ~7.8 (vs current 365.24)

The model’s proper-physics formula remains valid for projections up to about +3 Gyr forward; beyond that the polynomial saturates at the tidal-lock distance and the formula returns no value. For longer projections, the angular-momentum boundary itself takes over.

The model’s effective predictive domain spans about 7.5 Gyr — from Earth-Moon genesis to the formula’s +3 Gyr horizon at the tidal-lock distance. The modern epoch sits at about 61 % through this domain. Beyond +3 Gyr the proper-physics formula returns no value, though the Earth-Moon system continues to exist (and the Sun continues to shine) until the red-giant phase at +5 Gyr.

Future-projection uncertainty

The Farhat 2022 polynomial is calibrated against past anchors and extrapolated forward beyond J2000. For sub-200-Myr forward projections the model is well within accepted tidal-modelling ranges. At +1 Gyr the model’s projected H of 435,488 yr corresponds to ~32-hour LOD and ~421,000 km Moon distance — these values sit ~10 % above the predictions of independent tidal models (Bills & Ray 1999, Touma & Wisdom 1994) which predict slower late-future recession as tidal force weakens (∝ 1/r⁶) with Moon distance. The polynomial’s smooth-extrapolation behavior does not capture this physical slowdown. For internal lattice ratios the difference is irrelevant (every period scales together), but for absolute long-term forward use the +500 Myr to +1 Gyr range should be read as an upper-bound estimate, not a precise prediction.

10. Falsifiable predictions

The Expanding Solar System Resonance Theory (ESSRT) makes specific testable claims about how the lattice evolved and how it will continue to evolve.

1. Hadean Moon at the Roche limit. The proper-physics formula places the Moon at 3.22 Earth radii at Patterson’s Pb-Pb Earth age of 4.54 Gyr — within ~10 % of the Roche limit. No Hadean constraint was used in the fit. Any new radiometric or paleoclimate technique that refines either Earth age or Hadean Moon distance should remain consistent with this match.
2. Devonian H ≈ 309,083 yr. The model’s Devonian H produces 396.21 days/year, matching Wells 1963’s coral count of ~400 to within 1 %. New high-precision paleontological day-count techniques (e.g., expanded Torreites bivalve sampling) should reproduce this match.
3. Integer-label invariance. L1 lattice fits to paleoclimate spectra at Devonian, Permian, and Cretaceous epochs should find the same set of integers (n = 65 for obliquity main, n = 22/25/28 for short-eccentricity, etc.) — only with rescaled absolute periods.
4. Future tidal-lock asymptote at 87.1 Earth radii. Lunar laser ranging extrapolated forward through full tidal Q-decay should converge on the angular-momentum boundary at 555,623 km Moon distance.
5. Cheng cross-proxy persistence. The L1 lattice should continue to fit Cheng 2016’s Asian-Monsoon δ¹⁸O record at R² ≈ 0.68 as new high-precision speleothem chronologies extend or refine the record.
6. The classical precession of the equinoxes scales with H(t). Earth’s axial precession period (H/13) was ~23,776 yr at Devonian and ~5,372 yr at Hadean. The classical observable known since Hipparchus is not a fixed astronomical constant but a now-snapshot of a slowly-evolving cycle. Modern-era rate of change (~0.5 yr per Myr) is small but in principle observable in high-precision IAU precession-rate measurements over decades. Detail and the deep-time table: Precession §Precession through deep time.
7. Lunar precession scales as (H/H₀)² across geological time. The Moon’s apsidal period was ~9.60 yr at the Devonian and ~16.43 yr at −2.5 Gyr, with T_apsidal × H ≈ 2,966,688 yr² at every epoch (analogously T_nodal × H ≈ 6,241,326 yr²). Independent deep-time reconstructions of the lunar apsidal period (via spectral analysis of tidal rhythmites or cyclostratigraphic records sensitive to perigean modulation) should reproduce the H₀²/H(t) track within ~1%; a Phanerozoic apsidal period substantially off this track would falsify the invariant. Full derivation: §6 The Lunar Precession Invariant.

The full list of testable predictions (25 total, including these deep-time claims and 19 spanning near-term through climate timeframes) lives in Predictions.

See also

• Physical Origin — the mechanism that holds the lattice stable (action-angle closure for the climate lattice; KAM for planetary spacing)
• Climate Summary — Earth’s climate as one application of the modern lattice
• Climate Formula — the canonical L1 + L2 + L3 climate formula in detail
• Supporting Evidence — paleo-day-count validation tables (Wells, Williams, Pannella, de Winter, Cheng)
• Predictions — falsifiable predictions including deep-time claims