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.25~569.45~38.01
2100 AD~5,592.84~564.04~32.60

Falsifiability: if geocentric precession stays constant at ~5,604″/cy → standard GR supported; if it decreases from ~5,598.25 toward ~5,592.84″/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 peaked at exactly 6h in 1246.03125 AD and is now slowly shifting. By ~6,000 AD it will be ~5h58m22s.

PropertyValue
Mean RA at max declination (fixed frame)~5h48m50s / ~17h48m50s
Oscillation amplitude±11 minutes
Cycle period~41,915 years
Peak value~6h00m00s (reached in 1246.03125 AD)
Predicted value by 6,000 AD~5h58m22s
RA at maximum obliquity decreasing from peak at solstice–perihelion alignment epoch RA fluctuation pattern from obliquity and inclination interference

Current shift rate ~17 arcsec/century, well within modern astrometric precision. The mechanism and amplitude are independently confirmed by standard IAU 2006 precession theory: m_A = p_A × cos(ε) − χ_A inherits the obliquity period because it depends on cos(ε) (Capitaine et al. 2003 ; Laskar et al. 1993 ). A back-of-envelope calculation gives ~±10 minutes of RA — matching the model’s ±11 minutes. Canonical: 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.4393° 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,665 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.

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,312 yr around year 12,431 AD, then a rise to ~26,051 yr around year 32,341 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 ~5,372 yr at Hadean, ~23,776 yr at Devonian, ~~25,794 yr today, and rises to ~26,974 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 ~0.5 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 non-tidal Earth-rotation speedup component (the Munk-MacDonald hypothesis) on top, attributing the gap between pure-tidal prediction and the historical eclipse record to glacial isostatic adjustment + core-mantle coupling. The model predicts LOD will increase until ~30,000 AD, then decrease until ~2,000 AD, before increasing again — short-term fluctuations (like the 2020–present speedup) sit on top of this trend. Both sidereal day and stellar day follow the same cyclical pattern. The model’s pure-tidal ΔT formula has been empirically confirmed against 19 documented historical solar eclipses spanning -762 to 1654 CE (19/19 events visible at the documented site vs Stephenson empirical fit’s 17/19) — no non-tidal-speedup component is required by the historical record. Canonical: Timekeeping; empirical validation: Historical Eclipse Validation.

Solar day length cyclical pattern over Earth Fundamental Cycle

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: decrease until 30,000 AD.

Solar year in days decreasing over time

15. Sidereal year in seconds

Chapront et al.: slowly increasing until ~15,600 AD. The model: fixed at 31,558,149.76 seconds — the anchor point of the entire framework.

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.63603°. 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 Earth’s Hadean origin 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. Hadean Moon at the Roche limit

At Patterson’s Pb-Pb radiometric Earth age of 4.54 Gyr, the proper-physics two-layer LOD formula places the Moon at 20,532 km = 3.22 R_E — just outside the Roche limit (~2.9 R_E, where a fluid Moon would tidally disintegrate). No Hadean constraint was used in the fit; the formula is anchored at modern LOD = 24 hr and calibrated against Farhat 2022’s Phanerozoic deep-time tidal-evolution curve. Any future radiometric or paleoclimate technique that refines either Earth’s age or the Hadean Moon distance should remain consistent with this match. A measured Hadean Moon distance outside ~2.5–5 R_E at Patterson’s epoch falsifies the prediction.

21. Devonian H ≈ 309,083 yr matches Wells 1963 to within 1 %

The model’s Devonian (380 Ma) Earth Fundamental Cycle is H = 309,083 yr, producing 396.21 tropical days/year — matching Wells 1963’s directly-counted Devonian coral growth rings (~400 days/yr) to within 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).


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 variation (~30,000 → ~2,000 AD)MillenniaDiffers 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

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: