Skip to Content
📄 Fibonacci Laws — Read the paper
The ModelPredictions

Predictions

The Holistic model produces specific testable predictions that differ from conventional theory and can be checked against future observations. Predictions are grouped by timeline.

All predictions are measured from the 3D model. Every value on this page is read directly from the 3D Simulation  using objective measurement functions — not derived from theoretical assumptions. See Analysis & Export Tools for how measurements are taken.


Near-term predictions (decades)

1. Mercury’s “missing” perihelion precession will decrease

Mercury’s perihelion precession includes a ~43″/cy contribution attributed to General Relativity. The model attributes this to Earth’s wobble on its Axial Precession Orbit rather than space-time curvature. The geocentric rate should therefore decrease — by ~5.0″/century:

YearGeocentric (moving equinox)Heliocentric (ICRF)“Anomaly”
2000 AD~5,598.27~569.47~38.03
2100 AD~5,592.87~564.07~32.63

Falsifiability: if geocentric precession stays constant at ~5,604″/cy → standard GR supported; if it decreases from ~5,598.27 toward ~5,592.87″/cy → the model is supported.

Near-term test — BepiColombo: Mercury orbit insertion 21 November 2026; routine science from April 2027. The MORE instrument will measure Mercury’s orbit at 1–2 orders of magnitude better precision than MESSENGER. The model predicts ~574.61″/cy or lower vs MESSENGER’s 575.31″/cy — a 0.70″/cy difference, ~500× larger than MESSENGER’s uncertainty. The test is decisive if BepiColombo’s pipeline reports the raw measured perihelion advance rather than a GR-inclusive ephemeris fit total — methodology canonical at Mercury Precession §BepiColombo Test.

2. RA at maximum declination will shift from 6h

Coordinate-system note: in standard precessing equatorial coordinates the June solstice is at RA 6h by definition. This prediction refers to the Sun’s position in a fixed reference frame (ICRF) at maximum declination.

The Sun’s ICRF position at June solstice oscillates about a mean of 5h47m22s with amplitude ±3.20° (±12.8 min), a compound signature from two co-equal H-lattice cycles. It has sat near the top of that range throughout recorded history, reaching its maximum of 6h00m04s in 1675 AD — and it is only now beginning to decline, at ~5.5″/century.

The turnover is extremely flat: 6h00m03.1s at 1246.03125 AD, 6h00m04.0s at the 1675 peak, 6h00m03.4s at J2000 — about one second of RA across 750 years. The 1246.03125 AD perihelion–solstice alignment is the anchor epoch for H, not a feature of this curve: it sits at 99.0 % of maximum swing, indistinguishable from the peak at this resolution.

PropertyValue
Maximum reached1675 AD at 6h00m04s
Value at 1246.03125 AD6h00m03.1s (99.0 % of max swing)
Current rate (2000→2100)−5.5″/century
Oscillation amplitude±3.20° (±12.8 min RA)
Component 1 (inclination frame shift)~111,772 years (H/3)
Component 2 (obliquity)~41,915 years (H/8)
Period ratio8:3 (H/3 : H/8)
Component amplitudesequal, A/sin(ε̄) ≈ 1.60° each (A = 0.63604°)
Mean RA at max declination (fixed frame)5h47m22s / 17h47m22s
RA at maximum obliquity decreasing from peak at solstice–perihelion alignment epoch RA fluctuation pattern from obliquity and inclination interference

Falsifiable signature — two co-equal spectral peaks, not one. The prediction combines two co-equal H-lattice mechanisms with matched amplitudes:

  • H/3 inclination-frame shift (co-equal, ~96 arcmin peak): Earth’s orbital-plane inclination in ICRF oscillates ±0.63604° at the H/3 = ~111,772-yr period. This tilts the ecliptic frame relative to ICRF, shifting the RA at which the max-declination event occurs. Contribution A/sin(ε̄) ≈ 1.60°.
  • H/8 obliquity variation (co-equal, ~96 arcmin peak): Earth’s axial tilt oscillates over the H/8 = ~41,915-yr obliquity cycle with the same 0.63604° driver (the shared amplitude appears in the two-component obliquity formula ε(t) = ε̄ − A·cos(2π·t/(H/3)) + A·cos(2π·t/(H/8)) — see Obliquity & Inclination), contributing an equal A/sin(ε̄) ≈ 1.60° term with opposite sign. The two terms superpose as −sin(2π·t/(H/3)) + sin(2π·t/(H/8)) to give an exactly H-periodic composite signal with inner beat envelope at H/5 = ~67,063 yr, reaching ±12.5-min RA peak swing.
  • H/16 Equation of Center — residual sub-arcmin signal: Δλ_EoC = 2e·cos(ϖ) from perihelion precession would naively give ±8 min, but the max-declination event location absorbs the EoC displacement (SS coincides with true ecliptic longitude 90° regardless of eccentricity), leaving only ~5 arcsec residual at H/16 from the small e(t) variation.
  • Secondary mechanisms (deep-time H(t) drift, general precession H/13, tropical-calendar drift) contribute the remaining ~2 arcmin RMSE beyond the two-term closed form.

Traditional obliquity-only estimates (Laskar et al. 1993 ) predict a single spectral peak at the obliquity period. The Holistic model’s dual mechanism predicts two co-equal peaks in an 8:3 ratio at the H/3 (~111,772 yr) and H/8 (~41,915 yr) periods — with the H/16 EoC peak reduced to a sub-arcsec residual by the max-declination event definition. Testable via long-baseline VLBI or historical solstice records; current shift rate is ~5.5″/century, within modern astrometric precision but demanding a long baseline. Empirical validation: least-squares fit against model-computed SS events sampled at 1-year cadence across one full Earth Fundamental Cycle gives RMSE 2.1 arcmin for the two-term basis (verified via computeSolsticeRA in the reference implementation). Canonical derivation: Supporting Evidence §7.

WebGeocalc displays apparent pattern changes for Jupiter and Saturn perihelion precession over the next centuries. The model predicts the current trends continue as-is — Jupiter prograde at perihelion period 8H/39 = 68,783 yr and Saturn retrograde at −8H/65 = 41,270 yr (8H-lattice secular periods). Saturn’s observed ecliptic-retrograde rate ~-3,400″/cy (WebGeoCalc) is consistent with the model.

WebGeocalc prediction of Jupiter and Saturn perihelion precession changes Jupiter perihelion 1900–2500 AD continued trend Saturn perihelion 1900–2500 AD continued trend

The disagreement: standard theory attributes Saturn’s ecliptic-retrograde motion to a transient phase of the Great Inequality (~883-yr oscillation from the 5:2 resonance), reversing within ~450 years. The model treats it as permanent. For Jupiter the disagreement is in period — secular theory implies ~305,000 yr vs the model’s 68,783 yr; Jupiter’s observed inclination trend favours the shorter period (~3″/cy error vs ~8.5″/cy). Full analysis: Supporting Evidence §12 and §8.

4. No major Planet Nine (key differentiator)

The Batygin–Brown hypothesis posits an undiscovered 4–10 M_Earth planet at 290–700 AU shepherding extreme TNOs. A two-tier test against the model rejects all proposed candidates: (primary) Law-4 compliance — the observed eccentricity must match e_amp = K · sin(tilt) · √d / (√m · a^(3/2)); all proposed candidates fail by 4–7 orders of magnitude. (secondary) v-balance — adding a 9th body crashes the canonical closure for any mass above ~10⁻⁴ M_Earth. Combined threshold for compatibility at high-e ETNO orbits: ~2 × 10⁻¹⁰ M_Earth (~2-km rocky asteroid).

CandidateMass (M_E)Distance (AU)Best 9-planet balanceVerdict
Batygin & Brown 201610.07001.25%REJECT
Brown & Batygin 20216.238010.07%REJECT
Siraj et al. 20254.429013.19%REJECT
Pluto-mass test0.002246094.71%MARGINAL
Ceres-mass test0.0001646099.99%ACCEPT (v-balance only)

Under the primary Law-4 test, even the Ceres-mass candidate fails (its eccentricity 0.25 is 4,270× larger than the maximum Law-4 amplitude). Verdict: no candidate at proposed parameters is framework-compatible.

Near-term test — Vera Rubin Observatory (LSST 2025–2035):

LSST outcomeModel predictionConventional Batygin/Brown
M ≥ 1 M_Earth at 300–700 AU, e ≈ 0.2–0.6FALSIFIEDconfirmed
Ceres-mass body at 300–500 AU, high eFALSIFIEDweakened
Tiny body (≲ 2 × 10⁻¹⁰ M_Earth) at similar orbitconsistentn/a
No detection above ~10⁻¹⁰ M_Earth by 2035consistentFALSIFIED
ETNO clustering disappears with more discoveriesconsistentFALSIFIED

A single detection of any ≳ Ceres-mass body in a high-e Batygin–Brown-style orbit falsifies the model. The conventional hypothesis can be re-parameterised to fit nearly any detection — making it harder to falsify cleanly. In Popper’s terms, the model’s prediction is more vulnerable and therefore stronger. Full search, ETNO data, and alternative explanations: Planet Nine.

5. Small classical-belt KBO obliquity clustering (Law 4 bidirectional)

Law 4 is bidirectional. Solving the closure for sin(tilt) gives an obliquity prediction from a body’s mass, semi-major axis, eccentricity amplitude, and Fibonacci d-slot:

sin(tilt) = e_amp · √m · a^(3/2) / (K · √d)

For the 8 IAU planets this closes to <1%. For a representative 100-km classical-belt KBO at a ≈ 45 AU, e ≈ 0.05, mass ~10⁻¹² M☉, Law 4 predicts axial obliquity ≈ 36.6° at d = 55. The population-statistical claim: small cold-classical-belt KBOs (the regime where external forcing on e_amp is weakest) should cluster near this prediction, rather than show the uniform-random ~57° mean expected from Lambert’s law on a sphere.

LSST outcome (~10³ small-TNO rotation-pole measurements by 2030–2035)Verdict
Small classical-belt KBO obliquities cluster near ~30°–50°Bidirectional Law 4 confirmed beyond the 8 planets
Obliquity distribution is uniform-random (~57° mean)Bidirectional reading restricted to the 8-planet domain

The ~30% gap between the model prediction (~36.6°) and the random expectation (~57°) is resolvable with ~10² obliquity measurements at ±10° precision.

Scope caveat: no individually named TNO is per-body testable — every catalogued/named TNO (Pluto, Eris, Haumea, Makemake, Quaoar, Sedna, …) is too large at its observed eccentricity to be Law-4-admissible. Comets and main-belt asteroids fail too (Jupiter scattering, outgassing, Yarkovsky/YORP). The framework’s positive predictive zone is the sub-200 km low-e cold-classical-belt population (10⁴–10⁵ bodies enumerated by Col-OSSOS), tested statistically across many rotation-pole measurements.

Law-4 closure ≡ “cleared the neighborhood”: the framework’s planet/non-planet partition coincides with the IAU’s 2006 third criterion. Soter (2006) proposed a mass-discriminant Λ but it never became canonical; Law-4 closure provides a different quantitative formalisation that returns a clean pass/fail and matches the same 8-body partition. Agreement is not tautological: Planet Nine candidates fail Law-4 closure by 4–7 orders of magnitude (prediction #4), and the small-KBO obliquity-clustering signature here is the second non-circular extension.


Medium-term predictions (centuries)

6. Obliquity

Standard theory: obliquity (23.4394° in 2000 AD) decreases until ~13,900 AD, reaching ~22.6° (Chapront et al. & Laskar). The model agrees through ~13,700 AD; the next obliquity minimum lands at ~22.51° around year 13,664 AD, after which obliquity rises again on the ~41,915-year cycle. Both bottom out near the same epoch; they differ on what happens next — the model predicts a clean reversal back toward the mean; polynomial extrapolations diverge.

The full long-term envelope is ~22.21° – ~24.72°, with the high end (~24.72°) sitting slightly above Laskar’s standard maximum (~24.5°) — a distinct prediction worth testing against paleoclimate proxies.

The ~41,915-year cycle period quoted is the J2000-anchored value (H/8); at geological timescales the cycle itself drifts as H(t) evolves — at Devonian (380 Ma) H was ~306,189 yr (so H/8 obliquity cycle was correspondingly shorter); in 200 Myr H will be ~352,600 yr (and the cycle proportionally longer). See Expanding Resonance for the deep-time layer; the within-H Fourier pattern this prediction describes is essentially independent of that drift at the millennial scale.

Obliquity prediction showing decrease until minimum then rising

7. Longitude of perihelion

The model matches Meeus’s formula closely until ~3000 AD; afterward the two diverge. The model completes 360° in a mean period of ~20,957 yr; Meeus’s polynomial extrapolation deviates increasingly outside its fit window.

Longitude of perihelion matching Meeus for several millennia Long-term longitude of perihelion showing divergence from Meeus after several millennia

8. Gregorian calendar drift

The Gregorian year (365.2425 days) doesn’t match the actual solar year (~365.2422036 days):

  • By 6486 AD: June solstice on June 17, 20:00 UTC
  • By 11,725 AD: June solstice on June 17, 10:00 UTC (4 days earlier than today)

9. Analemma shape changes

The analemma shifts forward in time with the perihelion precession cycle. Its width changes with eccentricity (~20,957-year cycle); its length changes with obliquity (fluctuating between 22.21° and 24.72°).

Predicted analemma shape changes over perihelion precession cycles

10. Axial precession period reaches a minimum, then increases

Current Capitaine polynomial: axial precession period is decreasing. The model agrees on the current trend but predicts a reversal: minimum of ~25,314 yr around year 12,411 AD, then a rise to ~26,047 yr around year 32,343 AD, then continued oscillation on the perihelion precession cycle. The Capitaine polynomial does not model this oscillation. The underlying cause is the interaction between lengthening day and shortening solar year.

On geological timescales, the mean itself drifts. The model’s H/13 axial precession period was ~4,991 yr at Hadean, ~23,553 yr at Devonian, ~~25,794 yr today, and rises to ~27,123 yr at +200 Myr. The classical “precession of the equinoxes” — known since Hipparchus and traditionally treated as a fixed astronomical constant — is a now-snapshot of an H(t)-evolving cycle. The modern rate of change is ~6 yr per Myr: small but in principle observable in high-precision IAU precession-rate measurements over decades. Framework: Expanding Resonance.

Axial precession period decreasing to a minimum then increasing — an oscillation the Capitaine polynomial does not model

Long-term predictions (millennia)

11. Eccentricity (key differentiator)

Standard theory: Earth’s orbital eccentricity (0.01671 in 2000 AD) decreases toward ~0 by ~27,000 AD. The model predicts it reaches a minimum of ~0.0140 much earlier — around year 11,725 AD — and then increases again. Canonical: Eccentricity.

Eccentricity reaching minimum then rising, contradicting conventional theory

12. Inclination

Standard theory: Earth’s inclination to the invariable plane (1.57869° in 2000 AD) decreases, with no minimum predicted. The model predicts inclination decreases to ~0.845° in year 32,682 AD, then increases again on the ~111,772-year cycle.

Inclination decreasing to predicted minimum

13. Earth rotation / LOD / Delta-T

Standard theory: Earth’s rotation is monotonically slowing due to tidal friction. The empirical Stephenson polynomial encodes a Munk-MacDonald-scale (~5-6 ms/cy) non-tidal Earth-rotation speedup component on top, attributing the apparent gap to glacial isostatic adjustment + core-mantle coupling. The model predicts LOD monotonically increases under Driver 1 (Earth-Moon tidal recession) at a modern rate of ~1.7 ms/century, with a ~250 ms peak-to-peak Milankovitch oscillation superimposed from the refined α(t) coupling to the Climate Formula’s L1 orbital layer (see below). Short-term fluctuations (like the 2020–present speedup) sit on top of this trend. Both sidereal day and stellar day follow the same tidal-recession + Milankovitch pattern.

Two timescales coexist. At geological timescales (Myr+) the model agrees with standard tidal evolution — LOD monotonically increases per ESSRT Driver 1 (Earth-Moon tidal evolution), the same mechanism standard theory invokes. See Expanding Resonance for the deep-time layer. On top of that monotonic tidal drift, the refined α(t) coupling to the Climate Formula’s L1 orbital layer adds a ~250 ms peak-to-peak Milankovitch oscillation at ~100 kyr timescales — visible as wiggles on the tidal trend, matching known glacial-interglacial cycles. Standard theory’s polynomial extrapolation doesn’t model this Milankovitch oscillation; the model adds it as a specific physical prediction. The disagreement is therefore about the superimposed climate-cycle structure — not about the monotonic deep-time trend.

The model’s ΔT formula has been empirically confirmed against the historical record on two independent eclipse tracks, with a third, weaker comparison against a paleoclimate proxy: (a) a 26-event eclipse alignment audit on documented solar eclipses spanning -762 BCE to 2026 CE (18/26 events with framework umbra reaching the observation site within the ±4-hour scan window; 8/26 pure geographic misses); (b) the higher-resolution lunar timing test on 267 primary-source observations (20.2-min mean |residual|, matching NASA within 13 s and beating NASA on 117/267 events); and (c) an out-of-sample comparison against the Bond 2001 IRD drift-ice record at Pearson r = +0.36 in the 4,000 BC – 1,800 AD window — an open correspondence that fails its null tests, and which the two eclipse tracks do not depend on. The full Munk-MacDonald postulate (~5-6 ms/cy) is rejected by the historical record; a dominant GIA-scale channel (dLOD/dt = -0.35 ms/cy at J2000, Cox & Chao 2002 satellite anchor with Peltier ICE-6G factor-2.0) is included via the L1-orbital-coupled α(t) correction, and a smaller fractional non-tidal secular rate (~0.5 ms/cy) is detected in the residual but not currently modelled — see Lunar Eclipse Validation §1 for the quantitative statement and §6 for the three-component decomposition of the medieval residual. Canonical: Timekeeping; empirical validation: Solar Eclipse Validation and Lunar Eclipse Validation.

LOD is coupled to the same orbital forcing that drives climate. The L1-orbital-coupled α(t) correction above uses the same L1 orbital layer that fits the LR04 δ¹⁸O record in the Climate Formula — one mechanism, two observables. Ice sheets grow and retreat under Milankovitch orbital forcing on ~100 kyr timescales; that ice-mass redistribution modulates Earth’s polar moment α on the same timescale; and α modulates LOD. The model predicts that length-of-day carries a ~250 ms peak-to-peak Milankovitch oscillation in phase with the glacial cycle, superimposed on the smooth tidal-recession trend. LOD peaks (α max) fall at glacial extrema — MIS 6 (~140 ka BP), MIS 2 / LGM (~22 ka BP), and the projected next-glacial (~60,500 AD, matching Prediction 18) — and LOD troughs (α min) fall at interglacial extrema — MIS 7e (~215 ka BP), MIS 5e Eemian (~125 ka BP), and today’s Holocene. The underlying “ice → J₂ → α → LOD” mechanism is confirmed on decadal timescales by Cheng, Tapley & Ries 2011  (LAGEOS satellite J₂ shifting from linear GIA-driven decrease to acceleration around 1998, attributed to polar ice-mass loss). The 100-kyr extrapolation follows from the same physics but is not yet directly observable — modern satellite records span only ~50 years. The unifying claim: planetary gravity, climate, and Earth’s rotation are one system, connected by ice mass as the mediator. Predictions 13 and 18 are the same physical mechanism viewed through two observables (LOD vs δ¹⁸O). See Lunar Eclipse Validation §4 for the derivation.

Solar day length under refined α(t) coupling over −248 kyr to +102 kyr — model curve (blue) oscillates around the long-term climate mean (dashed purple), with MIS peak-age labels aligned to LOD extrema

The same coupling in the rate domain: over the last two glacial cycles the Tidal + GIA baseline rides above the tidal line through ice-sheet growth and dips sharply at the two deglaciations — with the LR04 record itself (inverted δ¹⁸O, the very dataset the L1 layer is fitted against) overlaid for direct comparison. The deepest negative-rate excursions land at Termination II (→ Eemian) and Termination I (→ Holocene).

Framework dLOD/dt from 200,000 BC to 2050 with the LR04 benthic stack (inverted δ¹⁸O, green) on the secondary axis — the Tidal + GIA baseline carries the ~41/100-kyr α(t) modulation, the 4-flag stack the millennial envelope, and the deepest spin-up excursions align with the two deglaciations

14. Solar year length in days

Laskar: solar year in days decreases until ~10,900 AD (assuming a fixed 86,400-second day). The model: the solar year in days decreases secularly under Driver 1 (LOD growth → fewer, longer days per year), with the H/8-band Fourier harmonics superimposed at the millennial scale — no near-term reversal.

At geological timescales the days-per-year ratio drifts on the same monotonic trend driven by ESSRT — both numerator (sidereal year in seconds, Driver 2) and denominator (LOD, Driver 1) evolve. At Devonian (380 Ma) the model gives ~399.96 tropical days/year (matching Wells 1963 coral growth bands to within 1 %; see Prediction 21). See Expanding Resonance for the full deep-time layer.

Solar year in days decreasing over time

15. Sidereal year in seconds

Chapront et al.: slowly increasing until ~15,600 AD. The model: held at the J2000 anchor value of 31,558,149.76 seconds within the modern-era scope, where it serves as the framework’s reference orbital period. At geological timescales the orbital period itself drifts via solar mass loss (Kepler’s third law) — see Expanding Resonance.

Standard theory treats ecliptic, axial, and other precession motions as largely unrelated. The model predicts they all follow a clear pattern repeating every ~20,957-year perihelion cycle.

All precession movements related in eccentricity-cycle pattern

Structural predictions

17. Invariable plane tilt

Earth’s path relative to the invariable plane has mean tilt ~1.48113° with amplitude ~0.63604°. The same structure governs all planetary motion: Jupiter’s precession determines Earth’s ecliptic precession, Saturn’s precession determines Earth’s obliquity cycle.

Earth's invariable plane tilt at mean inclination with characteristic oscillation amplitude Jupiter's invariable plane tilt Saturn's invariable plane tilt

Saturn is the only planet whose perihelion precesses retrograde in the ecliptic frame (confirmed by JPL WebGeoCalc at ~-3,400 arcsec/cy — canonical at Supporting Evidence §12) and the only anti-phase planet (cosine sign flipped). Saturn alone carries the entire anti-phase side; the seven other planets form the in-phase group. The angular-momentum-weighted inclination amplitudes cancel between groups to 99.9974%, keeping the invariable plane balanced; the same partition independently satisfies the eccentricity balance (99.8636%). Full derivation: Fibonacci Laws Derivation §Law 3. Visualise via Tools > Invariable Plane Inspector in the 3D Simulation .


Climate prediction

18. Next natural glaciation peak at ~60,500 AD

The canonical 3-layer Climate Formula extrapolated forward from t ≈ 2000 AD identifies the next predicted glacial maxima and interglacial peaks:

Climate Formula forward projection from −250 to +250 kyr. Past glacial maxima (LGM, MIS 6, MIS 8) correctly identified within ≤9 kyr; future markers: next glacial onset ~60,500 AD, warmest interglacial ~166,000 AD, strongest glacial ~198,500 AD.
Years from nowAD dateC(t) normalizedNote
~58,500~60,500 AD+2.27Next natural glaciation onset
~106,000~108,000 AD+0.24mild
~153,000~155,000 AD−0.99(local max, interglacial-range)
~164,000~166,000 AD−2.31Warmest interglacial in window
~196,500~198,500 AD+2.48Strongest glaciation in next 250 kyr

The signal C(t) is the normalised δ¹⁸O proxy from the post-MPT regime fit (negative = warmer/interglacial; positive = colder/glacial). The Holocene is correctly identified as interglacial; MIS 6 is placed within ~2 kyr; the LGM is predicted within ~9 kyr (the expected ice-sheet response lag).

Comparison with established forecasts:

FrameworkNext glacial onset (kyr from now)Mechanism
Berger & Loutre (2002) ~50Astronomical insolation + LLN-2D climate model
Loutre & Berger (2003) ~50–100Same + low-CO₂ scenarios
Tzedakis et al. (2012) ~50 (analog-bound)MIS 19c past-interglacial analog
Ganopolski et al. (2016) ~50 (no CO₂) / ~100 (moderate) / ≥500 (high)CLIMBER-2 with explicit CO₂ feedback
Holistic Climate Formula~5832-integer 8H lattice + L2 carbon thermostat + L3 step transitions

The Holistic ~58.5 kyr (~60,500 AD) prediction sits in the consensus range and matches Berger & Loutre 2002 within ~16%. What distinguishes it is the mechanism: standard frameworks attribute the 100-kyr cycle to direct eccentricity forcing; the model attributes it to the s₁ − s₄ nodal eigenmode beat at n = 25 = 107.3 kyr — a planet-pair orbital-plane coupling. Full comparison and mechanistic distinction: Climate Formula.

Orbital forcing is not climate: the formula captures the orbital component (L1) + silicate-weathering thermostat (L2) + discrete Cenozoic step transitions (L3). It does not model ice-sheet hysteresis, CO₂ amplification feedbacks beyond L2, regional asymmetries, or anthropogenic CO₂. The ~58,500 yr glacial-onset prediction is when the orbital clock makes a phase transition possible, not when surface climate necessarily follows — ice sheets carry thermal memory. Ganopolski et al. (2016)  found moderate-emission anthropogenic CO₂ may delay the next natural glaciation by 50+ kyr (high-emission ≥100+ kyr).

The pacing departs from the post-MPT 100-kyr regime: two strong glaciations at ~58.5 kyr and ~196.5 kyr from now (138 kyr apart, not the regular ~100-kyr drumbeat), with mostly weak intermediate wiggles. The L1 orbital component is exactly 8H-periodic and matches the late-Pliocene “41-kyr world” 8H ago. Physical interpretation and three climate-response scenarios: Climate Formula §Pacing shift.

19. Planet obliquity cycles (two-component structure)

Every planet with an obliquity cycle follows the same two-component formula as Earth — the sum of two cosines with equal amplitude, one at the ICRF perihelion period and one at the obliquity cycle period:

obliquity(t) = mean − A × cos(ICRF period) + A × cos(obliquity cycle)

Three obliquity cycles are already confirmed:

PlanetPredicted cycleObservedErrorSource
Mercury894,179 yr (8H/3)~895,000 yr0.2%Bills & Comstock 2005
Earth~41,915 yr (H/8)~41,000 yr2%Laskar+ 1993
Mars127,740 yr (8H/21)~124,800 yr2.4%Ward 1973; Laskar+ 2004

Three remain testable:

PlanetPredicted cycleCurrent literatureHow to test
Jupiter167,659 yr (H/2)“No regular cycle” (Saillenfest+ 2020)Long-term spin-axis integration
Saturn111,772 yr (H/3)“No regular cycle” (Saillenfest+ 2021)Long-term spin-axis integration
Uranus167,659 yr (H/2)“Frozen at ~98°” (Saillenfest+ 2022)Extremely long timescale simulation

Venus and Neptune have obliquity cycle = |ICRF perihelion period| (auto-derived from their ecliptic periods, tidally damped). The two-component formula cancels exactly, producing constant obliquity — consistent with observations (Venus 177°, Neptune ~28°). Two-component formula canonical at Obliquity §A Universal Pattern.


Deep-time predictions

These predictions extend the model’s testable claims across geological time, from the Earth-Moon genesis through the modern epoch to the far-future tidal-lock asymptote. They follow from the proper-physics two-layer LOD formula (Driver 1 — Earth-Moon tidal evolution) combined with the adiabatic invariant a × M_☉ = const (Driver 2 — solar mass loss). Framework: Expanding Resonance.

20. Earth-Moon genesis at the rigid Roche limit at the giant-impact age

The proper-physics two-layer LOD formula, run backwards, crosses the rigid Roche limit (9,471 km = 1.49 R_E) at ~4.498 Ga — the canonical giant-impact Moon-formation age (~4.5 Ga, dated independently by isotope geochemistry), ~40 Myr after Patterson’s Pb-Pb Earth age of 4.54 Gyr. No Hadean constraint was used in the fit; the formula is anchored at modern LOD = 24 hr with the LLR-observed recession rate and calibrated against Farhat 2022’s deep-time tidal-evolution curve. Any future refinement of the Moon-formation age or the genesis distance should remain consistent with this match. A genesis epoch outside ~4.4–4.55 Ga, or a genesis distance far from the Roche zone, falsifies the prediction.

21. Devonian H ≈ 306,189 yr matches Wells 1963 to within 0.1 %

The model’s Devonian (380 Ma) Earth Fundamental Cycle is H = 306,189 yr, producing 399.96 tropical days/year — matching Wells 1963’s directly-counted Devonian coral growth rings (~400 days/yr) to within 0.1 %. New high-precision paleontological day-count techniques (expanded Torreites bivalve sampling, recalibrated coral growth-band analysis, tidal-rhythmite inter-annual cycles) should reproduce this match. A measured Devonian count outside ~390–405 days/yr falsifies the prediction. Canonical validation table: Supporting Evidence §14.

22. Integer-label invariance across geological time

The L1 lattice integers (n = 65 for the climate-recorded obliquity main beat, n = 39 for Jupiter’s ecliptic perihelion, n = 16 for Earth’s perihelion harmonic, and the full 32-integer set) are structural constants of the solar system’s secular dynamics. The absolute periods rescale as H(t) expands across geological time, but the integer labels themselves do not change. Independent cyclostratigraphic L1 lattice fits to Devonian, Permian, and Cretaceous proxy spectra should find the same set of integers identified in the post-MPT LR04 fit — only with rescaled absolute periods. A deep-time spectrum that fits a substantively different integer set (e.g., a Devonian lattice missing n = 65 or with a new structurally-dominant integer absent post-MPT) falsifies the framework’s invariance claim.

23. Future tidal-lock asymptote at 87.1 R_E

The Earth-Moon system asymptotically approaches tidal lock at Moon distance 555,623 km = 87.1 R_E (currently 60.3 R_E), reached at roughly 50 Gyr from now — well beyond the proper-physics formula’s +3 Gyr predictive horizon and beyond the Sun’s red-giant phase at +5 Gyr. Lunar laser ranging extrapolated forward through full tidal-Q decay should converge on this angular-momentum boundary. A converged asymptote outside ~80–95 R_E falsifies the angular-momentum calculation.

24. Cheng cross-proxy persistence

The L1 lattice currently fits Cheng 2016’s Asian-Monsoon δ¹⁸O record at R² = 0.68 on its independent U-Th radiometric chronology (see Supporting Evidence §13). As new high-precision speleothem chronologies extend the record or refine its dating, the same 32-integer lattice should continue to fit with comparable R². A substantially degraded fit (R² dropping below ~0.5 on a refined or extended Cheng record) would falsify both the lattice’s cross-proxy persistence and the integer-label invariance claim above (§22).

25. Lunar precession scales as (H/H₀)² across geological time

The Lunar Precession Invariant T_apsidal × H = 2,966,728 yr² (year-units, equinox-of-date frame) is held exact at every epoch by the framework’s (H/H₀)² scaling, with the analogous T_nodal × H = 6,241,369 yr² for the nodal mode. The Moon’s apsidal period was ~9.69 yr at the Devonian and ~17.15 yr at −2.5 Gyr (J2000 anchor: 8.848 yr); the nodal period was ~20.38 yr at Devonian (J2000 anchor: 18.613 yr). The H² scaling form matches the leading m² order of Brown’s lunar perturbation theory (the historical Newton-Clairaut problem); the framework sharpens Brown’s theory by anchoring the J2000 magnitude from observation and adopting the structural (H/H₀)² scaling for deep-time evolution, with no polynomial corrections.

Independent deep-time reconstructions of the lunar apsidal period — via spectral analysis of tidal rhythmites, cyclostratigraphic records sensitive to perigean modulation, or future high-precision paleo-tidal sediment studies — should reproduce the H₀²/H(t) track within ~1%. A Phanerozoic apsidal period substantially off this track would falsify the invariant. Full derivation: Expanding Resonance §6: The Lunar Precession Invariant.

26. Deglacial spin-up leads the interglacial optimum

Sustained negative dLOD/dt occurs nowhere in the framework’s 200-kyr record except during the two major deglaciations — GIA-channel rate minima of ~−1.2 ms/cy at ~10,600 BC (Termination I) and ~127,800 BC (Termination II). The mechanism makes this a leading indicator by construction: the GIA rate term tracks the melting rate — α follows ice volume (more ice ⇒ mantle displaced equatorward ⇒ larger α), so dα/dt ∝ d(ice)/dt and melting drives α down — so the rate minimum marks maximum melting speed, and peak interglacial warmth follows ~5–10 kyr later, when melting completes (LR04 warm peaks at ~121,000 BC and the Holocene plateau — the derivative leads the integral). See the 200-kyr chart in Prediction 13 and Timekeeping §When the day got shorter.

Falsifiable test: any eclipse-independent paleo-rotation record should show the spin-up episode preceding the Holocene Climatic Optimum and the Eemian peak by several kyr. A rotation record showing the spin-up concurrent with or after peak warmth would falsify the α(t) coupling’s phase structure.

Scope caveat: causality runs climate → rotation (ice → J₂ → α → LOD); the rotational signature indicates warming already in progress — a leading indicator of peak warmth, not a driver of warming. At millennial scale the stack channel follows the concurrent sign rule on the accumulated offset — warm ↔ Σstack > 0, i.e. a longer day — not on its rate; the lagged relationship is specific to the deglacial GIA channel, whose sign is the opposite (see Timekeeping §Two channels).


27. The cardinal-point braid law: seasonal asymmetries are one rotating vector

The four cardinal points (solstices and equinoxes) carry no independent physics: their event-time offsets from the mean tropical year are one complex function — the perihelion vector of the H/16 cycle — read through four fixed rotations, δ_X = Im[e^(i·n·λ_X)·W_n(t)] with λ = 0°/90°/180°/270° (SS→AE→WS→VE). The law is verified at four levels (Days & Years §The braid law): 90° carrier quadrature with equal amplitudes to 0.08°/0.29% (unenforced by the fit, so unfakeable), the order-2 pair pattern, the full fine-structure spectrum captured by coefficients shared across all four points, and statistical equality of the four points’ residuals. The decisive sign test: the geometry requires the spectrum to be quadrature-locked counter-rotating; both rotation senses were fitted with identical freedom, and only the required sense captures anything — the same discrimination logic as heterodyne sideband detection.

Falsifiable, forward-looking content: the law fixes the differential drift rates of the four per-point tropical years with zero per-point freedom. At J2000 the model predicts (relative to the common mean): SS +0.4, WS −0.4, VE +1.5, AE −1.6 s/century — two near-opposite pairs summing to zero, the quadrature structure expressed in rates. Century-scale precision timing of equinoxes and solstices (IMCCE/Meeus-class ephemeride series) must reproduce this pattern; the per-point year lengths themselves already match the IAU values at J2000 to 0.2–3.2 s. The seasonal pattern rotates with perihelion: by ~11,725 AD the winter/summer solstice roles swap. And any future high-precision decomposition of cardinal-point timing residuals must find the fine structure in the counter-rotating sideband only — a coherent co-rotating component at measurable amplitude would falsify the one-vector geometry outright.


Verification pathways

PredictionTimeframeType
Mercury geocentric precession decreaseDecadesDiffers from GR
RA shift from 6hCenturiesNew observable
Jupiter/Saturn perihelion trendDecadesDiffers from WebGeocalc
No major Planet Nine (Law-4 compliance, ≳ 2 × 10⁻¹⁰ M_Earth at high-e ETNO orbits)Decades (LSST 2030–2035)Key differentiator
Small classical-belt KBO obliquity clustering (~36° at 100 km / 45 AU / e≈0.05)Decades (LSST 2030–2035)Bidirectional Law 4
Axial precession reversalCenturiesDiffers from Capitaine
Eccentricity minimum at 11,725 ADMillenniaKey differentiator
LOD growth-rate modulation (stack trough ~2,206 AD → peak ~2,740 AD)CenturiesDiffers from standard theory
Next natural glaciation ~60,500 AD (~58,500 yr ahead)MillenniaClimate Formula forward projection
Invariable plane tilt 1.48113°StructuralNew observable
Jupiter/Saturn/Uranus obliquity cyclesLong-termTestable by N-body integration
Deglacial spin-up leads interglacial optimum by ~5–10 kyrDeep time (paleo-rotation records)New observable

The four quickest tests:

  1. BepiColombo data (~2027) — model predicts ~574.61″/cy or lower vs MESSENGER’s 575.31″/cy (if the pipeline reports raw measurement; see methodology)
  2. Vera Rubin Observatory / LSST (2030–2035) — delivers two verdicts simultaneously: (a) Planet Nine — any ≳ Ceres-mass detection at 300–700 AU with e ≈ 0.2–0.6 falsifies the model; no detection above ~10⁻¹⁰ M_Earth falsifies the conventional hypothesis; (b) small-KBO obliquity clustering — ~10³ rotation-pole measurements test whether sub-200 km cold-classical-belt KBOs cluster near ~36° rather than uniform-random ~57°.
  3. The RA at max obliquity shifting from 6h
  4. Mercury’s geocentric precession — model predicts decrease; GR predicts constant.

Most predictions on this page can be verified directly using the formulas at Formulas.


← Supporting Evidence

Last updated on: