Predictions
The Holistic Universe Model produces specific, testable predictions about astronomical values that differ from conventional theory. These predictions can be tested against future observations.
All predictions are based on measurements from the 3D model. Every prediction on this page originates from values measured directly in the 3D Simulation — not from theoretical assumptions. The simulation produces raw data (year lengths, day lengths, precession periods, orbital parameters) using objective measurement functions; the predictions here follow from that data. See Analysis & Export Tools for how measurements are taken.
Predictions are grouped by timeline - how soon they could potentially be tested. Some require decades of observation; others require centuries.
Near-Term Predictions (Decades)
These predictions could show measurable differences within the next 50-200 years.
1. Mercury’s “Missing” Perihelion Precession Will Decrease
Current theory: Mercury’s perihelion precession includes a ~43 arcsecond/century contribution from relativistic effects (General Relativity). This value has been stable at ~43″ since Newcomb’s 1882 measurement.
Model prediction: The “missing” precession may be due to Earth’s wobble on its Axial Precession Orbit rather than space-time curvature.
What actually changes: Mercury’s precession can be measured against two reference directions (both in the ecliptic plane):
- ~575″/century — relative to fixed stars (ICRF), the inertial frame defined by distant quasars
- ~5,604″/century — relative to the moving vernal equinox (~575 + ~5,028.8 equinox drift from Earth’s axial precession)
The equinox-based value (~5,604″) is what was historically measured on Earth. The model predicts this value will decrease by ~5.0″/century:
| Year | Relative to moving equinox (geocentric) | Relative to fixed stars / ICRF (heliocentric) | “Anomaly” (derived) |
|---|---|---|---|
| 2000 AD | ~5,598.26″ (~569.46 + 5,028.8) | ~569.46″ | ~38.02″ |
| 2100 AD | ~5,592.85″ (~564.05 + 5,028.8) | ~564.05″ | ~32.61″ |
The ICRF value (~575″) is not directly measured — it is derived from the geocentric measurement (~5,604″) by subtracting the equinox drift (~5,028.8″). The “anomaly” (~43″) is then obtained by subtracting the Newtonian contribution (~532″). Since the Newtonian contribution is constant, the “anomaly” decreases if the geocentric rate decreases.
Falsifiability: This is a clear test between the two interpretations:
- If geocentric precession remains constant at ~575 + 5,028.8 = ~5,604″/century: Standard GR explanation is supported, model’s alternative is refuted
- If geocentric precession decreases from ~5,598.26″ toward ~5,592.85″/century: Model’s explanation is supported
Near-term test — BepiColombo: ESA’s BepiColombo mission arrives at Mercury on 21 November 2026 (delayed from Dec 2025 due to thruster issues), with routine science operations starting April 2027. The MORE (Mercury Orbiter Radio science Experiment) instrument will measure Mercury’s orbit with 1–2 orders of magnitude better precision than MESSENGER. The model predicts BepiColombo will measure ~574.61″/century or lower versus MESSENGER’s 575.31″/century — a difference of 0.70″/century, which is ~500× larger than MESSENGER’s measurement uncertainty (±0.0015″/century). This makes it the quickest available test of the model’s Mercury prediction, provided BepiColombo’s analysis pipeline reports the raw measured perihelion advance rather than a GR-inclusive ephemeris fit total — see the BepiColombo test methodology.
See Mercury Precession: The BepiColombo Test for a detailed breakdown of the two possible outcomes and what each would mean for the model. See also Scientific Background for the full scientific discussion including academic critiques.
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.
Current observation: The Sun’s maximum declination (June solstice), when expressed in ICRF coordinates, occurs near Right Ascension 6h.
Model prediction: This ICRF position peaked at exactly 6h in 1246.03125 AD and is now slowly shifting. By ~6,000 AD, the June solstice RA (in ICRF) will be ~5h58m22s.
| Property | Value |
|---|---|
| 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 |
The underlying pattern follows a ~41,915-year cycle driven by obliquity and inclination interference:
What this predicts: The model claims that the Sun’s celestial position at the moment of maximum declination, when expressed in ICRF coordinates (not precessing equatorial coordinates), varies slightly over millennia. This would manifest as a very small timing offset in when the Sun reaches its northernmost point relative to the fixed stars.
Verification method: Testing requires:
- Expressing the Sun’s position in ICRF at each solstice
- Tracking whether the ICRF position varies systematically over centuries
- Precision of ~1 arcsecond over centuries (currently achievable)
The current shift rate is ~17 arcseconds per century (accelerating as we move away from the 1246.03125 AD peak), well within modern astrometric precision.
Independent confirmation from standard theory: The IAU 2006 precession framework (Capitaine et al. 2003 ) gives the general precession in right ascension as m_A = p_A × cos(ε) − χ_A. Because this depends on cos(ε), and obliquity oscillates with a ~41,040-year period (Laskar et al. 1993 ), the precession rate in RA itself oscillates with the same period. A back-of-envelope calculation confirms the amplitude at ~±10 minutes of RA — matching this prediction. See Supporting Evidence for the full analysis.
3. Jupiter and Saturn Perihelion Trends Will Continue
Current prediction (WebGeocalc): Jupiter and Saturn’s perihelion precession rates are predicted to change pattern.
Model prediction: The current trends will simply continue as-is without pattern change.
The model predicts perihelion precession periods of H/5 = 67,063 years for Jupiter (prograde) and H/8 = 41,915 years for Saturn (retrograde) in the ecliptic frame, connected through the Fibonacci identity 3 + 5 = 8 (see Fibonacci Laws). Saturn’s observed ecliptic-retrograde rate is ~-3,400″/cy (WebGeoCalc), consistent with the model’s prediction.
The key disagreement: Standard theory attributes Saturn’s ecliptic-retrograde motion to a transient phase of the Great Inequality — an ~883-year oscillation from the Jupiter-Saturn 5:2 resonance — which should reverse within ~450 years. The model treats it as permanent. For Jupiter, the disagreement is in period: secular theory implies ~305,000 years vs the model’s 67,063 years. Jupiter’s observed inclination trend favours the shorter period (~3″/cy error vs ~8.5″/cy). See Supporting Evidence §13 and §9 for the full analysis.
4. No Major Planet Nine (Key Differentiator)
Conventional hypothesis (Batygin & Brown 2016-2025): An undiscovered 9th planet of mass 4-10 M_Earth at 290-700 AU shepherds extreme trans-Neptunian objects (ETNOs). This hypothesis is itself contested within the astronomy community (OSSOS null result, no detection after a decade of searching).
Model prediction: A two-tier test rejects all proposed Planet Nine candidates. Primary: Law-4 compliance — the observed eccentricity must be consistent with the framework’s amplitude prediction e_amp = K · sin(tilt) · √d / (√m · a^(3/2)). All Batygin-Brown / Siraj candidates fail by 4–7 orders of magnitude. Secondary: canonical v-balance search — adding a 9th body to the 8-planet balance crashes the closure for any mass above ~10⁻⁴ M_Earth. The combined Law-4 threshold for compatibility at high-e ETNO orbits sits near ~2 × 10⁻¹⁰ M_Earth (~2-km rocky asteroid).
Quantitative result (secondary v-balance test): A canonical 7,558,272-configuration search (matching the model’s own balance-search implementation that produced the 767 surviving 8-planet configurations) found that the best achievable 9-planet balance with a 5 M_Earth body at 450 AU is min(Law 3, Law 5) ≈ 8.96% — vs. the current 8-planet baseline of 99.9975% on Law 3 and 99.8632% on Law 5. Even with full freedom to regroup every planet’s in-phase/anti-phase assignment and reassign all Fibonacci d-values, no configuration accommodates a multi-Earth-mass 9th body.
| Candidate | Mass (M_E) | Distance (AU) | Best 9-planet balance | v-balance verdict |
|---|---|---|---|---|
| Batygin & Brown 2016 | 10.0 | 700 | 1.25% | REJECT |
| Batygin & Brown 2021 | 6.2 | 380 | 10.07% | REJECT |
| Siraj et al. 2025 | 4.4 | 290 | 13.19% | REJECT |
| Pluto-mass test | 0.0022 | 460 | 94.71% | MARGINAL |
| Ceres-mass test | 0.00016 | 460 | 99.99% | ACCEPT under v-balance only |
Under the primary Law-4 compliance test, even the Ceres-mass candidate fails (the body’s orbital eccentricity 0.25 is 4,270× larger than its maximum predicted Law-4 amplitude). The two-tier verdict is that no candidate at proposed parameters is framework-compatible.
Near-term test — Vera Rubin Observatory (LSST): The 10-year LSST survey (2025-2035) is expected to either find Planet Nine or rule it out conclusively. Either outcome teaches us something definitive:
| LSST outcome | Model prediction | Conventional Batygin/Brown |
|---|---|---|
| M ≥ 1 M_Earth body found at 300-700 AU with e ≈ 0.2-0.6 | FALSIFIED (Law-4 violation by 5+ orders) | confirmed |
| M ~ Ceres-mass body found at 300-500 AU with high e | FALSIFIED (Law-4 violation by 3+ orders) | weakened |
| Tiny body (≲ 2 × 10⁻¹⁰ M_Earth) found at similar orbit | consistent | n/a |
| No detection above ~10⁻¹⁰ M_Earth by 2035 | consistent | FALSIFIED |
| ETNO clustering disappears with more discoveries | consistent | FALSIFIED |
A single detection of a ≳ Ceres-mass body in a high-e Batygin-Brown-style orbit falsifies the model — either via Law-4 violation (primary, very strong) or via v-balance disruption (secondary, still strong). The conventional hypothesis, by contrast, can be re-parameterized 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.
See Planet Nine: A Falsifiable Prediction for the full canonical search, ETNO data, alternative explanations (observational bias + ancient stellar flyby), and references.
5. Small Classical-Belt KBO Obliquity Clustering (Bidirectional Law 4 — Testable by LSST)
Conventional theory: Orbital eccentricity (a property of the orbit) and axial obliquity (a property of rotation) are treated as independent state variables. Their long-term evolutions are weakly coupled through tidal and spin-orbit-resonance effects, but no first-principles theory links them as a single algebraic identity.
Model prediction: Law 4 — e_amp = K · sin(tilt) · √d / (√m · a^(3/2)) — is bidirectional. Solving for sin(tilt) produces a quantitative obliquity prediction from a body’s mass, semi-major axis, secular eccentricity amplitude, and Fibonacci d-slot:
sin(tilt) = e_amp · √m · a^(3/2) / (K · √d)
For the 8 IAU planets this identity closes to <1% (already validated). For non-planet bodies the framework partitions sharply: bodies whose actual e_amp ≈ Law-4 intrinsic prediction are “Law-4-compliant” (their dynamics evolve under their own (m, a, d) algebra); bodies whose actual e_amp is dominated by external forcing (resonance pumping, scattering, galactic tides) are “non-planet” in the framework’s quantitative sense — see Fibonacci Laws Derivation §Law-5 residual decomposition. Pluto’s actual e_amp ≈ 0.025 decomposes into Law-4 intrinsic (~0.001) + Neptune-resonance external (~0.024) — a 1:24 split, with Law 4 correctly capturing the intrinsic baseline.
The testable claim: 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 is that small cold-classical-belt KBOs (the regime where external forcing on e_amp is weakest) should have axial obliquities clustering near arcsin(K · √d · a^(-3/2) · √m / e) per body’s parameters — i.e., a predicted clustering near ~36° rather than a uniform-random distribution.
| Test outcome (LSST, ~10³ small-TNO rotation-pole measurements by 2030–2035) | Verdict |
|---|---|
| Small classical-belt KBO obliquities cluster near the Law-4 prediction (~30°–50°) | Bidirectional Law 4 confirmed beyond the 8 planets |
| Obliquity distribution is uniform-random (~57° mean, Lambert’s law on a sphere) | Bidirectional reading restricted to the 8-planet domain |
The random expectation (⟨sin(tilt)⟩ = π/4 ≈ 0.785, mean obliquity ≈ 57.3°) differs from the model prediction (sin(tilt) ≈ 0.596, ~36.6°) by ~30% — a sample of ~10² obliquity measurements at ±10° precision would resolve it.
Scope caveat — no individually named TNO is per-body testable. Every catalogued/named TNO (Pluto, Eris, Haumea, Makemake, Quaoar, Sedna, …) is too large (>700 km) at its observed eccentricity to be Law-4-admissible — they all sit above the admissibility curve. Comets (67P) and main-belt asteroids (Ceres) also fail: Law 4 admits them numerically but predicts ~0° and ~17° respectively while measurements give 52° and 4°. These bodies are externally forced (Jupiter scattering, non-gravitational outgassing, Yarkovsky/YORP) — the prediction does not apply at the per-body level. The framework’s positive predictive zone is exclusively 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. Bodies that have cleared their neighborhood evolve under their own dynamics — their e_amp is set by Law-4 intrinsic algebra. Bodies that have not are externally forced. Soter (2006) proposed a mass-discriminant Λ to formalise the IAU criterion quantitatively, but Λ never became canonical; Law-4 closure provides a different quantitative formalisation that returns a clean pass/fail and matches the same 8-body partition. The 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.
Near-term test — Vera Rubin Observatory (LSST): The same survey that delivers the Planet Nine verdict (prediction #4) will also constrain small-TNO rotation poles via lightcurve inversion. By ~2035, the bidirectional reading of Law 4 will either gain empirical support beyond the 8 planets (positive clustering near 36°) or be restricted to the 8-planet domain (uniform-random distribution). Either outcome teaches us something definitive — and the same instrument disambiguates two distinct framework predictions in one survey.
Medium-Term Predictions (Centuries)
These predictions would become measurable over several centuries.
6. Obliquity
Current theory: Obliquity (23.4393° in 2000 AD) will decrease until year ~13,900 AD, reaching a minimum of ~22.6° (per Chapront et al. & Laskar formulas).
Model prediction: The model agrees with current theory through year ~13,700 AD. The next obliquity minimum lands at ~22.51° around year 13,665 AD, after which obliquity rises again on its ~41,915-year cycle. Both model and standard theory bottom out near the same epoch (~13,665–13,900 AD); they differ on what happens next — the model predicts a clean reversal back toward the mean, while polynomial extrapolations diverge.
The full long-term oscillation envelope over the Earth Fundamental Cycle is ~22.21° – ~24.72°, which closely tracks the simple mean ± 2A theoretical range of the two-cosine model — the H/3 (inclination) and H/8 (obliquity) components align tightly enough each cycle for the envelope to essentially reach its extrema. The high end (~24.72°) sits slightly above Laskar’s standard maximum (~24.5°), which is a distinct model prediction worth testing against paleoclimate proxies.
7. Longitude of Perihelion
Current theory: J. Meeus’s formula calculates longitude of perihelion, with solstice-perihelion alignment in 1246.03125 AD.
Model prediction: The model matches Meeus’s formula closely until ~3000 AD. After that, values diverge significantly.
8. Gregorian Calendar Drift
Current situation: The Gregorian calendar year (365.2425 days) doesn’t match the actual solar year (~365.2422036 days).
Model prediction:
- By year 6486 AD: June solstice will be on June 17, 20:00 UTC
- By year 11,725 AD: June solstice will be on June 17, 10:00 UTC (4 days earlier than today)
9. Analemma Shape Changes
Current understanding: The analemma (figure-8 pattern of the Sun’s position) depends on perihelion position, eccentricity, and obliquity.
Model prediction: The analemma will shift forward in time with the perihelion precession cycle. Its width changes with eccentricity (~20,957-year cycle), while its length changes with obliquity (fluctuating between 22.21° and 24.72°).
10. Axial Precession Period Will Reach a Minimum, Then Increase
Current theory: The axial precession period is decreasing (Capitaine et al. 2003 polynomial formula).
Model prediction: The axial precession period is currently below the mean (~25,771 years vs mean ~25,794 years) and still decreasing. Both the model and the Capitaine formula agree on this current trend. The key difference: the model predicts the period will reach a minimum of ~25,312 years around year 12,431 AD, then rise again to a maximum of ~26,051 years around year 32,340 AD, then continue oscillating on the perihelion precession cycle. The Capitaine formula is a polynomial extrapolation that does not model this oscillation.
The key difference: Both agree the period is currently decreasing. The model predicts an eventual reversal (oscillation) that the Capitaine polynomial does not capture. The underlying cause is the interaction between lengthening day and shortening solar year.
Long-Term Predictions (Millennia)
These predictions distinguish the model from conventional theory over thousands of years.
11. Eccentricity (Key Differentiator)
Current theory: Earth’s orbital eccentricity (0.01671 in 2000 AD) will decrease toward ~0 by year ~27,000 AD.
Model prediction: Eccentricity will reach minimum of ~0.0140 much earlier - around year 11,725 AD - and then increase again.
Why this matters: This is a sharp divergence from current theory. If eccentricity begins rising around 11,725 AD, it would be strong evidence for the model.
12. Inclination
Current theory: Earth’s inclination to the invariable plane (1.57869° in 2000 AD) will decrease, but no minimum is predicted.
Model prediction: Inclination will decrease to minimum ~0.845° in year 32,682 AD, then increase again in a ~111,772-year cycle.
13. Earth Rotation / Length of Day / Delta-T
Current theory: Earth’s rotation is generally slowing due to tidal friction, causing the mean solar day to grow longer than 86,400 SI seconds.
Model prediction: LOD will slightly increase until ~5,823 AD, then decrease until ~22,718 AD, before increasing again. Short-term fluctuations (like the 2020–present speedup) occur around this long-term trend.
Delta-T implications: The model predicts LOD varies cyclically over millennia, which would affect the long-term accumulation of Delta-T.
Sidereal day and stellar day implications: Both sidereal day and stellar day follow the same cyclical pattern as solar day.
14. Solar Year Length in Days
Current theory (Laskar): Solar year in days is slowly decreasing until year 10,900 AD (based on fixed 86,400-second day).
Model prediction: Solar year in days will decrease until 12,394 AD.
15. Sidereal Year in Seconds
Current theory (Chapront et al.): Sidereal year in seconds is slowly increasing until year 15,600 AD.
Model prediction: The sidereal year in seconds is fixed at 31,558,149.76 seconds - it’s the anchor point of the model.
16. All Precession Movements Are Related
Current theory: Different precession movements (ecliptic, axial, etc.) are largely unrelated to each other.
Model prediction: All precession movements are related and follow a clear pattern repeating every ~20,957-year perihelion cycle.
Structural Predictions
These predictions relate to the model’s framework rather than time-varying parameters.
17. Invariable Plane Tilt
Model prediction: Earth’s path relative to the invariable plane has a mean tilt of ~1.48113° with amplitude ~0.63603°. This specific value should be observable and shared by all planetary movements.
Additionally, Jupiter and Saturn precession will be found to be directly connected:
- Jupiter’s precession = Determines Earth’s Ecliptic precession
- Saturn’s precession = Determines Earth’s Obliquity cycle
Saturn is unique among the planets: it is the only planet whose perihelion precesses obviously retrograde in the ecliptic frame (confirmed by JPL WebGeoCalc at ~-3,400 arcsec/century; see Supporting Evidence §13), and it is the only planet that is anti-phase — its cosine sign is flipped relative to all other planets. This places Saturn alone in the anti-phase group, while all seven other planets belong to the in-phase group. The full anti-phase alignment — Saturn at maximum while all others at minimum — occurs once per Solar System Resonance Cycle (8H). These two groups are not arbitrary: the angular-momentum-weighted inclination amplitudes of the in-phase group and the anti-phase group cancel to 99.9975%, keeping the invariable plane balanced. Saturn single-handedly carries the entire anti-phase side of this balance. The same two groups also satisfy an independent eccentricity balance condition (99.8632%). See Fibonacci Laws: Derivation for the full derivation.
In the 3D Simulation , you can explore this balance visually via Tools > Invariable Plane Inspector.
Climate Prediction
18. Temperature Trends
Model prediction: As the inclination tilt effect decreases and the axial tilt effect approaches its midpoint, the warmer period will end, transitioning toward a longer ice age period.
Climate is influenced by many factors (solar cycles, volcanic activity, human impact). This prediction concerns the long-term orbital contribution to climate patterns.
19. Next Natural Glaciation Peak at ~38,000 yr from now
Model prediction: The 8H Orbital Forcing Formula (25 integer-divisor components fit on LR04, R² = 0.232) extrapolates forward from t = 0 (≈ 2000 AD) to identify both the peak interglacial warmth ahead and the next predicted natural glacial maxima:
| Years from now | AD date | Orbital signal C(t) | Phase |
|---|---|---|---|
| 0 (today) | 2000 AD | −0.44 | Holocene interglacial |
| ~5,700 | ~7,700 AD | −0.56 | Peak orbital interglacial warmth — cooling begins after this |
| ~32,300 | ~34,300 AD | 0.00 | Orbital signal crosses from interglacial-favoring to glacial-favoring |
| ~38,000 | ~40,000 AD | +0.25 | Next natural glaciation onset |
| ~81,500 | ~83,500 AD | +0.45 | moderate glacial |
| ~110,500 | ~112,500 AD | +0.32 | mild glacial |
| ~154,500 | ~156,500 AD | +0.67 | strong glacial |
| ~194,500 | ~196,500 AD | +0.92 | Strongest in next quarter-million years |
The orbital signal C(t) is the normalised δ¹⁸O proxy (negative = warmer/interglacial; positive = colder/glacial). The model places peak orbital warmth around ~7,700 AD, after which sustained orbital cooling begins. The signal becomes glacial-favoring (crosses zero) around ~34,300 AD, and the first major glacial peak follows ~6,000 years later at ~40,000 AD.
Validation against past glacials: the formula correctly identifies the Holocene as interglacial (C(0) negative), places MIS 6 within 2 kyr (predicted 138 kyr BP vs. observed ~140 kyr BP), and predicts the Last Glacial Maximum at 29 kyr BP — 9 kyr earlier than the observed LGM at ~20 kyr BP. The 9 kyr lag is expected: orbital forcing predicts the insolation drive, not the ice-sheet response; peak ice volume lags peak orbital cooling by several kyr (standard climate physics).
Comparison with established literature: Consistent with the classical Berger & Loutre (2002) prediction (~50 kyr to next glaciation) within model uncertainty.
Orbital forcing is not climate: The formula captures the orbital-forcing component only — not ice-sheet hysteresis, carbon-cycle feedbacks, or anthropogenic CO₂. The ~5,700 yr peak-warmth and ~38,000 yr glacial-onset predictions are when the orbital signal shifts; the actual surface-temperature response can lag by thousands of years (ice sheets carry thermal memory) and can be displaced or overridden by non-orbital forcings. Ganopolski et al. (2016) found that moderate-emission anthropogenic CO₂ scenarios may delay the next natural glaciation by 50+ kyr. These predictions describe when the orbital clock makes a phase transition possible, not when surface climate necessarily follows.
Pacing shift toward obliquity-band intervals. The intervals between the predicted glacial peaks above are roughly 43, 29, 44, 40, 49 kyr — clustered near the obliquity band (~40 kyr), not the ~100-kyr pacing that dominated the past 700 kyr of LR04. The forward projection therefore looks more like a return to the pre-MPT “41-kyr world” of the Early Pleistocene than a continuation of the late-Pleistocene 100-kyr regime. The same orbital prediction was reached from a different angle by Berger & Loutre (2002) . See Orbital Forcing Formula §“Pacing shift” for the physical interpretation and three climate-response scenarios.
20. Planet Obliquity Cycles (Two-Component Structure)
Model prediction: Every planet with an obliquity cycle follows the same two-component structure as Earth. The obliquity is the sum of two cosines with equal amplitude — one at the ICRF perihelion period (inclination component, negative sign) and one at the obliquity cycle period (positive sign):
obliquity(t) = mean − A × cos(ICRF period) + A × cos(obliquity cycle)
Three obliquity cycles are already confirmed by observation:
| Planet | Predicted cycle | Observed | Error | Source |
|---|---|---|---|---|
| Mercury | 894,179 yr (8H/3) | ~895,000 yr | 0.2% | Bills 2005 |
| Earth | ~41,915 yr (H/8) | ~41,000 yr | 2% | Laskar+ 1993 |
| Mars | 127,740 yr (8H/21) | ~124,800 yr | 2.4% | Ward 1973; Laskar+ 2004 |
Three remain as testable predictions:
| Planet | Predicted cycle | Current literature | How to test |
|---|---|---|---|
| Jupiter | 167,659 yr (H/2) | “No regular cycle” (Saillenfest+ 2020) | Long-term spin-axis integration |
| Saturn | 111,772 yr (H/3) | “No regular cycle” (Saillenfest+ 2021) | Long-term spin-axis integration |
| Uranus | 167,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 tidally damped at 177°, Neptune frozen at ~28°).
Mean obliquity: The formula midpoint (around which the two-cosine oscillation is centered) typically shifts slightly from the J2000 snapshot. Most planets have midpoints within ±0.2° of their J2000 values. See Obliquity: A Universal Pattern for the full analysis.
Summary: Verification Pathways
| Prediction | Timeframe | Type |
|---|---|---|
| Mercury geocentric precession decrease | Decades | Differs from GR prediction |
| RA shift from 6h | Centuries | New observable |
| Jupiter/Saturn perihelion trend | Decades | Differs from WebGeocalc |
| No major Planet Nine (≳ 2 × 10⁻¹⁰ M_Earth at high-e ETNO orbits, Law-4 compliance) | Decades (LSST 2030-2035) | Key differentiator |
| Small classical-belt KBO obliquity clustering (Law 4 bidirectional, ~36° predicted for 100 km / 45 AU / e≈0.05) | Decades (LSST 2030-2035) | Bidirectional Law 4 |
| Axial precession reversal | Centuries | Differs from Capitaine |
| Eccentricity minimum at 11,725 AD | Millennia | Key differentiator |
| LOD variation (~5,823→~22,718 AD) | Millennia | Differs from current theory |
| Next natural glaciation ~38,000 yr ahead (~40,000 AD) | Millennia | 8H formula forward projection |
| Invariable plane tilt 1.48113° | Structural | New observable |
| Jupiter/Saturn/Uranus obliquity cycles | Long-term | Testable by N-body integration |
The quickest ways to test this model would be:
- BepiColombo data (~2027): The model predicts ~574.61″/century or lower versus MESSENGER’s 575.31″/century — a 0.70″/century difference, ~500× larger than measurement uncertainty (if the analysis pipeline reports the raw measurement; see the BepiColombo test methodology)
- Vera Rubin Observatory / LSST (2030-2035) — delivers two distinct framework verdicts. (a) Planet Nine: a single detection of any ≳ Ceres-mass body at 300-700 AU with e ≈ 0.2-0.6 falsifies the model via Law-4 violation; conversely, no detection above ~10⁻¹⁰ M_Earth in high-e orbits falsifies the conventional Planet Nine hypothesis (see Planet Nine). (b) Small-KBO obliquity clustering: ~10³ small-TNO rotation-pole measurements will test whether sub-200 km cold-classical-belt KBOs cluster near the Law-4 bidirectional prediction (~36° for 100 km / 45 AU / e≈0.05) rather than the uniform-random ~57° expectation.
- Noticing the RA at max obliquity beginning to shift from 6h
- Monitoring Mercury’s geocentric precession — the model predicts decrease; GR predicts it stays constant
Only time will tell if these predictions prove correct.
Many of these predictions can be verified using the formulas on the Formulas page.