Supporting Evidence

The Holistic Universe Model makes claims that challenge several established theories. This page collects external scientific evidence from current science — published papers, unresolved problems, and recent observations — that independently support or align with the model’s framework.

1. The 100,000-Year Problem (Still Unsolved)

The dominant ~100k-year glacial cycle visible in ice core records is one of the most persistent unsolved problems in paleoclimatology. The Holistic Universe Model proposes this signal reflects the inclination precession cycle (~111,772 years), not eccentricity.

Why eccentricity is problematic

Eccentricity changes Earth’s annual insolation by only ~0.2% — far too small to drive major ice ages without invoking unverified amplification mechanisms. Three specific problems persist:

1. Spectral mismatch: Eccentricity’s spectrum shows a split peak at ~95k and ~125k years. But the climate record shows a single narrow peak near ~100k years. These don’t match.

2. The 400k-year absence: Eccentricity’s theoretically strongest component (~400k years) is largely absent from climate records of the past 1.2 million years. If eccentricity drives ice ages, its dominant cycle should appear in the data.

3. The Mid-Pleistocene Transition: Around 1 million years ago, glacial cycles shifted from ~41k years (obliquity-dominated) to ~100k years — with no change in orbital forcing. Multiple competing hypotheses have been proposed; as of 2025, none are certain. This remains “one of paleoclimatology’s great unsolved puzzles.”

Peer-reviewed support: Muller & MacDonald (1997)

The proposal that orbital inclination — not eccentricity — drives the ~100k-year cycle was published in the Proceedings of the National Academy of Sciences:

“The shape of the peak is incompatible with both linear and nonlinear models that attribute the cycle to eccentricity.” — Muller & MacDonald, PNAS 94(16), 8329–8334 

Their spectral and bispectral analyses showed that Earth’s orbital inclination relative to the invariable plane provides a better match to both the shape and phase of the climate signal.

Recent research (2024–2025)

The 100,000-year problem remains actively debated:

• Barker et al. (2025, Science, 387, eadp3491): Investigated the distinct roles of precession, obliquity, and eccentricity in Pleistocene glacial cycles — still unable to resolve which parameter dominates.
• Mitsui et al. (2025, Earth System Dynamics, 16, 1569–1584): Found that “the ~100 kyr spectral peak actually aligns with the 95 kyr eccentricity peak” — showing that even the peak identification itself is debated.
• Lisiecki (2023, Nature Geoscience): “Precession pacing of Late Pleistocene ice-sheet changes” — found precession plays a more important role than obliquity during Late Pleistocene ice-sheet changes, further complicating the standard picture.

How the model addresses this

The model proposes the “~100k-year” signal is actually ~111,772 years — the inclination precession period. The ~10% discrepancy between ~100k and ~112k may fall within ice core dating uncertainties, particularly since:

• The spectral peak in climate data spans 80–120 ka
• Many deep-time chronologies rely on orbital tuning (adjusting dates to match Milankovitch predictions), which is circular when testing Milankovitch theory
• Non-orbitally-tuned dating methods (O₂/N₂ ratio, U-Th speleothems) could independently test whether the true period is closer to ~100k or ~112k

Why the ~112k cycle only emerged at the MPT

The inclination precession cycle (H/3 = ~111,772 yr) is not a recent phenomenon — it is a formation-epoch feature of the solar system. Earth’s orbital inclination relative to the invariable plane was set when the protoplanetary disk dissipated ~4.5 billion years ago, frozen into place by the same KAM-optimal process that organized all Fibonacci relationships (see Physical Origin). No known mechanism — impacts, stellar encounters, chaotic diffusion — can change a planet’s orbital plane on million-year timescales. The inclination cycle has been operating continuously for billions of years.

The question is not “what started the ~112k cycle?” but “why did it become climatically visible only ~1 Myr ago?” Two mechanisms are supported by evidence:

1. Interplanetary dust concentration (Farley 1995, Nature 376, 153): Helium-3 measurements in deep-sea sediments show a real increase in interplanetary dust accretion beginning at ~1 Ma. The cause is likely an asteroid family breakup that created new dust concentrated near the invariable plane. As Earth’s orbit oscillates above and below this plane with the H/3 period, a denser dust concentration would make the inclination cycle climatically relevant for the first time — the orbital forcing didn’t change; the medium it acts through did. Muller & MacDonald (1997) originally proposed dust accretion as the climate mechanism for inclination forcing; while their specific model was questioned, Farley’s ³He evidence for a dust increase at the MPT remains unchallenged.

2. Ice sheet threshold (Willeit et al. 2019, Science Advances 5, eaav7337): Progressive CO₂ decline and removal of easily-erodible regolith allowed ice sheets to grow past a critical size where they could survive obliquity maxima. This “silenced” the ~41k-year obliquity pacemaker, allowing the longer ~112k-year inclination signal — always present in the orbital dynamics — to emerge as the dominant climate cycle. The orbital forcing didn’t change; the climate system’s sensitivity did.

Both mechanisms are consistent with the model’s framework: the Fibonacci orbital balance is a permanent feature of formation, while the MPT marks when the inclination cycle became detectable in climate records.

2. Fibonacci Ratios in Orbital Mechanics (KAM Theory)

The model’s 13:3 Fibonacci ratio between axial precession (~25,794 yr) and inclination precession (~111,772 yr) may appear to be a numerical coincidence. However, there is a rigorous theoretical reason for Fibonacci ratios to appear in stable orbital systems.

The KAM Theorem

The Kolmogorov–Arnold–Moser (KAM) theorem (1954–1963) proves that in perturbed dynamical systems, orbits with “most irrational” frequency ratios are maximally stable against perturbation. The key insight:

• Orbits with frequency ratios that are simple fractions (like 2:1 or 3:1) create resonances — repeated gravitational kicks that destabilize the orbit
• Orbits with “irrational” frequency ratios avoid these resonances
• The golden ratio (φ ≈ 1.618), to which successive Fibonacci ratios converge, is the most irrational number in a precise mathematical sense — it is hardest to approximate by ratios of small integers

This means orbits with golden-ratio-related frequencies are the last to become unstable under perturbation.

Observational evidence

Fibonacci ratios appear throughout the solar system and beyond:

• Pletser (2019, Astrophysics and Space Science 364:158): Orbital period ratios in solar planetary and satellite systems preferentially cluster near Fibonacci fractions (~60% vs ~40% for non-Fibonacci). These orbits are associated with more regular, less inclined, and more circular configurations.

• Aschwanden & Scholkmann (2017, New Astronomy 58:107): Found that the most prevalent harmonic ratios in 73% of 932 exoplanet pairs are Fibonacci fractions (2:1, 3:2, 5:3).

• Kirkwood Gaps: The asteroid belt shows dramatic gaps at simple integer resonances with Jupiter (3:1, 5:2, 7:3, 2:1) — while the regions between these resonances are stable. This is KAM theory in visible action.

• Saturn’s rings: Show gaps and density waves at rational resonances with Saturn’s moons — another dramatic demonstration.

What this means for the model

The model’s 13:3 ratio is therefore consistent with the maximally stable orbital configurations predicted by dynamical systems theory. The solar system has naturally evolved toward these configurations over billions of years.

The fact that the same Fibonacci numbers (3, 5, 8, 13) that divide the Earth Fundamental Cycle also appear in exoplanetary systems strengthens the case that this is a fundamental feature of gravitational dynamics, not a fitting artifact.

The model extends beyond period ratios: six independent Fibonacci Laws connect planetary precession periods, eccentricities, and inclination amplitudes through Fibonacci numbers and the mass-weighted quantity $e \times \sqrt{m}$, predicting orbital elements for all 8 planets with zero free parameters. See the full technical derivation for details.

3. Earth’s Rotation Speedup (2020–present)

The model predicts that Earth’s Length of Day (LOD) varies cyclically over millennia. LOD will slightly increase until ~5,823 AD, then decrease until ~22,718 AD. Short-term fluctuations are superimposed on this long-term trend.

What was observed

Starting in 2020, Earth unexpectedly began rotating faster:

• 2020: The 28 shortest days since atomic clocks were invented
• July 5, 2024: The shortest day ever recorded — 1.66 milliseconds under 24 hours
• 2019–2024: Average LOD shifted from +0.39 ms to a multi-year negative trend relative to 86,400 seconds

Scientists are puzzled

This speedup was not predicted by standard models. Nick Stamatakos of the IERS Directing Board stated  they “run into trouble predicting more than six months or one year ahead.” The ENSO cycle, core-mantle coupling, and atmospheric angular momentum explain some short-term variation, but the multi-year trend remains poorly understood.

How the model aligns

The model predicts LOD will slightly increase until ~5,823 AD before entering a longer decreasing phase until ~22,718 AD. The 2020–present speedup is a short-term fluctuation — qualitatively consistent with the model’s prediction that LOD varies cyclically, but not evidence of the long-term trend reversal itself.

4. Day Length Stalled for 1 Billion Years

A landmark 2023 paper in Nature Geoscience fundamentally challenges the assumption that Earth’s rotation has slowed monotonically due to tidal friction.

The finding

Mitchell & Kirscher (2023) analyzed geological constraints on Precambrian day length and found that Earth’s day length stalled at approximately 19 hours for roughly 1 billion years during the mid-Proterozoic (2.0–1.0 Ga). They proposed that atmospheric thermal tides from solar heating balanced the decelerative torque of lunar oceanic tides, temporarily stabilizing Earth’s rotation.

Why this matters for the model

This finding proves two things:

1. LOD dynamics are more complex than simple tidal deceleration — additional mechanisms can influence or even reverse the tidal slowing
2. Cyclical LOD behavior is physically possible — if atmospheric tides could halt rotational slowing for a billion years, other mechanisms could create cyclical variations

The model proposes a millennial-scale LOD variation superimposed on the long-term tidal trend. The Mitchell & Kirscher finding establishes that such complex rotational dynamics are not unprecedented.

Reference: Mitchell & Kirscher, 2023, Nature Geoscience 16, 567 

5. Solar Oblateness Uncertainty and Mercury

The standard test of General Relativity through Mercury’s perihelion precession assumes that all non-GR contributions (planetary perturbations, solar oblateness) are precisely known. Recent research questions this assumption.

The problem

The Sun’s gravitational quadrupole moment (J₂) — caused by its oblateness — contributes to Mercury’s precession. However:

• J₂ is not constant: It varies with the solar magnetic activity cycle (~11 years)
• Historical measurements disagree: Published J₂ values have ranged from ~10⁻⁵ (older oblateness-based estimates) to ~10⁻⁷ (modern helioseismology) depending on the method
• J₂ mimics GR: The solar oblateness contribution has the same temporal signature as the relativistic precession, making them difficult to separate

The risk for BepiColombo

A 2022 study found that if a periodic J₂ component exceeding 0.04% of J₂ exists and is not accounted for, it could falsely confirm or contradict General Relativity in BepiColombo’s measurements. This creates a systematic uncertainty in the standard Mercury GR test that is rarely discussed.

The model does not claim GR is wrong — but it notes that the standard Mercury test has an unresolved systematic uncertainty that weakens its status as a clean confirmation.

Reference: The Influence of Dynamic Solar Oblateness on Tracking Data Analysis, MDPI Remote Sensing, 2022 

6. BepiColombo: The Near-Term Test

The ESA/JAXA BepiColombo mission provides the most immediate opportunity to test the model.

Updated timeline

The test

The Mercury Orbiter Radio science Experiment (MORE) will measure Mercury’s orbit with 1–2 orders of magnitude better precision than MESSENGER.

Note: this comparison assumes BepiColombo’s analysis pipeline reports the raw measured perihelion advance. If a GR-inclusive ephemeris fit is used (as in Park et al. 2017 , where Mercury’s orbit is estimated jointly with the PPN parameter β), the model-predicted drift may be absorbed into the fit parameters rather than appearing as a clean change in the reported total.

This is a near-term observable test, provided BepiColombo’s analysis pipeline reports the raw measured perihelion advance rather than a GR-inclusive ephemeris fit total (see Mercury Precession Test methodology). Under that condition, either the rate has measurably decreased (supports the model) or it hasn’t (supports GR).

See Mercury Precession: The Model’s Alternative for the full analysis.

7. Independent Dating Methods

Several dating methods exist that are completely independent of orbital tuning — meaning they could test whether the ~100k-year glacial signal is actually ~111,772 years.

Speleothems (cave deposits)

• Dated by uranium-thorium (U-Th) decay — no orbital assumptions
• Cheng et al. (2016, Science) found ~100k patterns with consistent timing
• The exact spectral peak position (~100k vs ~112k?) deserves reanalysis

O₂/N₂ ratio dating

• Kawamura et al. (2007, Nature): trapped air O₂/N₂ ratio correlates with local summer insolation
• Constrains timing to precession cycles (~23 ka), not 100k cycles
• Provides an independent orbital constraint without assuming eccentricity drives climate

Tidal rhythmites

• Sedimentary records preserving ancient tidal cycles
• Provide constraints on ancient Length of Day
• Show discrepancies with simple tidal deceleration models
• Support the existence of complex rotational dynamics

The opportunity

A spectral reanalysis of non-orbitally-tuned climate records could distinguish between the ~100k and ~112k hypotheses. The data exists — it needs to be analyzed with this specific question in mind.

8. Solstice RA Oscillation (Confirmed by Standard Precession Theory)

The model predicts that the Sun’s Right Ascension at maximum declination (June solstice), expressed in ICRF coordinates, oscillates with a ~41,915-year period and ±11 minutes amplitude — peaking at exactly 6h in 1246 AD.

Standard theory confirms the mechanism

The IAU 2006 precession framework (Capitaine et al. 2003 ) gives the general precession in right ascension as:

m_A = p_A × cos(ε) − χ_A

Because m_A 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. The RA of the solstice in ICRF doesn’t advance linearly — it speeds up and slows down with the obliquity cycle.

Quantitative agreement

The envelope calculation: obliquity oscillates over ~22.21° – ~24.72°, giving a ~1.8% variation in cos(ε). At m ≈ 46.1”/yr, this produces ~0.83”/yr rate variation. Accumulated over a half-cycle (~20,500 yr), this yields ~10 minutes of RA — matching the model’s ±11 minutes.

Why this matters

This effect is implicit in standard precession equations but has never been separately highlighted as an observable prediction. The model identifies it as a measurable consequence of the obliquity cycle; standard theory independently confirms both the mechanism and the magnitude. The current shift rate of ~17 arcseconds/century (accelerating as we move away from the 1246 AD peak) is well within modern astrometric precision.

9. Jupiter and Saturn Secular Eigenfrequencies Are Stable (Laskar)

The model predicts that Jupiter’s and Saturn’s perihelion precession trends will simply continue as-is — without the pattern change that short-term ephemeris tools like WebGeocalc appear to show. Laskar’s secular theory independently confirms that the underlying precession rates are constant over millions of years.

The Great Inequality: why short-term tools show “pattern changes”

Jupiter and Saturn are locked in a near 5:2 mean-motion resonance. This produces the Great Inequality — a ~883-year quasi-periodic oscillation first identified by Kepler and explained by Laplace in 1786. The oscillation is large enough to dominate instantaneous perihelion rates over century timescales:

Saturn’s fitted perihelion rate changes sign between the two intervals. This is not a physical inconsistency — it demonstrates that the Great Inequality’s amplitude in Saturn’s perihelion longitude is comparable to or larger than the secular drift over century timescales. JPL’s own documentation warns: “The elements are not valid outside the given time-interval over which they were fit.”

Any tool that displays the full numerical ephemeris over centuries (including WebGeocalc) will show these oscillations as apparent “pattern changes” in the precession rate.

Secular theory: the long-term dynamics are stable

Laskar’s orbital solutions (La2004 , La2010 ) decompose planetary eccentricity evolution into eigenmodes with frequencies g₁ through g₈. Each planet’s actual perihelion motion is a superposition of all eigenmodes — the instantaneous precession rate is the angular velocity of the combined eccentricity vector, not any single eigenfrequency. For Jupiter, the dominant mode is g₅ (4.257″/yr) with a significant g₆ admixture (amplitude ratio ~2.8:1). For Saturn, g₆ (28.245″/yr) dominates but g₅ contributes nearly as strongly (ratio ~1.4:1), which is why Saturn’s instantaneous rate can even reverse sign.

Laskar found that g₅, g₆, and s₆ are “practically stable over time” — remaining essentially constant over at least 50 million years. The giant planet orbits are far less chaotic than the inner planets. Since the eigenfrequencies that govern Jupiter and Saturn are stable, the long-term secular dynamics — and the average precession trends built from them — do not change.

Laplace’s proof: the Great Inequality averages to zero

In 1786, Laplace proved that the Great Inequality is truly periodic — it does not produce a net secular drift in orbital elements. Over ~10 full cycles (~9,000 years), the oscillation averages out completely. Brouwer & van Woerkom (1950)  later showed that eliminating the Great Inequality from the Hamiltonian introduces small corrections to the secular eigenfrequencies (additional modes g₉, g₁₀), but these modify without destabilizing the fundamental rates.

Inclination cross-validation: H/5 outperforms ~305,000 years

An independent test comes from Jupiter’s inclination trend. The model’s perihelion precession period (H/5 = 67,063 years) is coupled to the inclination dynamics through the Fibonacci resonance structure. When this period is used, the predicted inclination trend matches JPL’s observed values with an error of ~3 arcseconds/century. When the secular theory’s g₅ eigenfrequency period (~305,000 years) is used instead, the inclination trend error jumps to ~8.5 arcseconds/century — nearly 3× worse.

This cross-validation against independent JPL inclination data suggests the model’s H/5 period may better describe Jupiter’s dynamics than the g₅ eigenfrequency alone. The tension between these two periods — and the fact that JPL’s own inclination data favors the shorter one — deserves further investigation.

What this means for the model

The model’s J2000 values are based on the trend from 1900 to 2000 from WebGeocalc , which better captures the actual direction of perihelion movement than single-epoch snapshots (see Mercury Precession for details). No observational data spans the thousands of years needed to directly observe a “pattern change.” Predictions about future trend reversals — from any model, including this one — are necessarily theoretical extrapolations. The model’s approach is to take the observed 100-year trend as the most reliable starting point.

The model’s prediction that perihelion trends “continue as-is” is further supported by two independent lines of evidence:

1. Laskar’s secular theory confirms the underlying eigenfrequencies are stable over millions of years — the Great Inequality is periodic, not a genuine trend change
2. JPL inclination data independently favors the model’s H/5 period over the theoretical ~305,000-year g₅ period

Ascending node periods: integer divisors of 8H

Each planet’s ascending node regression period takes the form 8H/N for an integer N, with Jupiter and Saturn locked to a shared N=36. Across all 7 fitted planets, the 8H/N integers reproduce JPL’s J2000-fixed-frame ascending-node trends with a cumulative residual error of ~5.8″/century (≈0.8″/century per planet).

Why this matters: The model derives all eight ascending-node periods from a single constant (H = 335,317), with each rate as an integer divisor of 8H. Laskar’s secular theory measures these as 8 independent eigenfrequencies with no known structural relationship to each other.

Observability limitation: These periods describe motion over 50,000–2,000,000 year timescales. With ~4,000 years of recorded astronomy, the periods cannot be verified by direct observation of a complete cycle.

10. Obliquity Amplitude: Berger (1978) Dominant Term

The model claims both the axial tilt effect and the inclination tilt effect oscillate by the same amplitude of ±0.63603°, combining to produce the full obliquity range of ~22.21° – ~24.72°. Independent support comes from the standard Fourier decomposition of obliquity.

Berger’s obliquity series

Berger (1978) decomposed Earth’s obliquity into 47 quasi-periodic terms. The dominant term (frequency s₃ + k) has:

This dominant amplitude of 0.684° is within 8% of the model’s 0.63603°. The five largest terms are:

The dominant term overwhelmingly controls the obliquity signal — roughly 3× larger than the next term. Berger & Loutre (2001) confirmed the value is essentially unchanged in updated solutions, and Laskar et al. (2004) numerical integrations produce the same spectral peak.

The physical mechanism matches

The frequency s₃ + k arises from exactly the two motions the model identifies:

• k (~50.47”/yr): the axial precession rate — Earth’s spin axis precessing prograde
• s₃ (~−18.85”/yr): the eigenmode most strongly associated with Earth’s orbital plane — precessing retrograde

Because these move in opposite directions, their rates add: 50.47 + 18.85 = 69.32”/yr for the meeting frequency, giving a half-cycle of 1,296,000/69.32 ≈ 18,700 years and a full cycle of ~37,400 years. The standard ~41k-year period includes contributions from the secondary terms (s₄ + k, s₆ + k) which shift the effective average period upward.

The model’s mechanism — two counter-rotating precessions producing a combined obliquity cycle — is the same physics that Berger’s frequency decomposition describes mathematically.

What is not explained

No published work explains from first principles why the dominant amplitude is specifically ~0.684° (or the model’s 0.63603°). It emerges from the full solution of the coupled spin-orbit-planetary perturbation equations. The amplitude depends on planetary masses, orbital configurations, and the Moon’s stabilizing torque (Laskar, Joutel & Robutel 1993). Without the Moon, obliquity could vary chaotically between 0° and 85°.

Climatic precession comparison

Berger’s climatic precession spectrum (e × sin ϖ̄) has a multi-peak structure distributed across ~19-24 kyr, with no single dominant term. The proper comparison frame is the Solar System Resonance Cycle (2,682,536 yr) — the period over which every planetary cycle returns to alignment in the model. Each Berger peak corresponds to an integer-fraction divisor of this cycle.

Every Berger climatic-precession peak matches an integer-fraction divisor of the Solar System Resonance Cycle within <0.4%, including Saturn (which falls outside the conventional ~19-24 kyr “climatic precession band” but is recovered by the same framework):

The structural decomposition: each integer n splits as n = 104 + δ, where:

• 104 = 8 × 13 is the integer corresponding to axial precession (k = H/13 of the Earth Fundamental Cycle gives 104 sub-divisions in one resonance cycle)
• δ is the planet’s eigenfrequency contribution: Jupiter +9, Mercury +12, Venus +16, Earth +36, Mars +37, Saturn +59

This mirrors the secular-theory decomposition climatic precession = k + g_j, but expressed entirely as integer fractions of the resonance cycle.

The model’s natural Fibonacci centroid sits at integer n = 128, equivalent to H/16 = ~20,957 yr in the Earth Fundamental Cycle (since 128 = 8 × 16 = 8 × (13+3), corresponding to k + the inclination-precession contribution H/3). The Berger spread is centered around n = 128 with planetary offsets that are partially Fibonacci-derivable:

What this establishes: the model’s H/n hierarchy correctly predicts the centroid of Berger’s climatic precession spectrum (H/16, equivalent to integer n = 128 of the resonance cycle), and the largest planetary offsets (Mars +13, Venus −8) are direct Fibonacci numbers. The remaining inner-planet offsets (±12, −15) are within 1-2 of Fibonacci numbers; Saturn’s larger +35 offset is non-Fibonacci, which may reflect that Saturn’s contribution falls outside the conventional climatic precession band entirely.

Open question: whether the planetary eigenfrequencies (δ = 9, 12, 16, 36, 37, 59 measured from the axial integer 104) can be derived from a Fibonacci-organized planetary mass distribution remains an open question. The fact that 3 of the 6 centroid-relative offsets are at or adjacent to Fibonacci numbers (±8, ±12, +13) suggests this is a tractable extension rather than a fundamental limitation for the inner planets and Mars; Saturn’s larger offset (+35 from centroid, +59 from axial) is the outlier.

Planetary perihelion periods cluster in the same band

A separate observation reinforces the structural alignment: each planet’s perihelion period in ICRF (model values from model-values.ts) sits in the same 16-42 kyr band that contains the Berger climatic-precession spectrum and the obliquity cycle:

Excluding Earth (the outlier at the inclination cycle), every planet’s ICRF perihelion period falls within 16-42 kyr — the same band that contains both Berger spectral families.

A triple identity at H/8 = ~41,915 yr is particularly striking:

• Jupiter’s ICRF perihelion period = H/8 (resonance cycle integer n = 64)
• Saturn’s ecliptic perihelion period = H/8 (retrograde, n = 64)
• Obliquity cycle period = H/8 (Berger’s dominant s₃ + k term)

The same divisor — H/8 — describes Jupiter’s perihelion motion in inertial space, Saturn’s perihelion motion in Earth’s ecliptic frame, and Earth’s obliquity oscillation. This is a structural identity, not a numerical coincidence: the obliquity cycle in the model is mathematically Saturn’s ecliptic perihelion period, with Jupiter’s ICRF perihelion sharing the same divisor on the inertial side. It supports the interpretation that the H/n hierarchy organizes the entire planetary system, not just Earth.

References:

• Berger, A. (1978). “Long-term variations of daily insolation and Quaternary climatic changes.” J. Atmos. Sci., 35, 2362–2367.
• Berger, A. & Loutre, M.F. (2001). “Amplitude and Frequency Modulations of the Earth’s Obliquity.” J. Climate, 14(6), 1043–1054.
• Laskar, J. et al. (2004). “A long-term numerical solution for the insolation quantities of the Earth.” A&A, 428, 261–285.
• Laskar, J., Joutel, F. & Robutel, P. (1993). “Stabilization of the Earth’s obliquity by the Moon.” Nature, 361, 615–617.

11. Why Equal Amplitudes? The Physics of Balanced Systems

The model’s claim that both obliquity components oscillate by the same amplitude (±0.63603°) is not just a simplification — it reflects a physical principle. Several independent branches of physics predict that equal amplitudes are the preferred state for coupled counter-rotating systems.

The coupled oscillator argument

In classical mechanics, two coupled oscillators produce normal modes. The key result: when two coupled oscillators are effectively identical (same restoring torque, same effective moment of inertia), their normal modes have exactly equal amplitudes. Unequal amplitudes arise only when the oscillators are fundamentally asymmetric.

The model’s obliquity formula has two cosine terms of equal amplitude: the inclination tilt effect at the ICRF perihelion period (~111,772 yr, H/3 — driven by planetary perturbations) and the axial tilt effect at the obliquity cycle period (~41,915 yr, H/8 — the Fibonacci beat of the H/5 ecliptic precession and H/3 inclination precession, corresponding to Berger’s s₃ + k beat). Equal amplitudes imply these two mechanisms are effectively equivalent oscillators — a physically meaningful claim about the dynamics, not just a parameter choice. The ±0.63603° amplitude is not observable on the H/13 (~25,794 yr) equinox precession cycle, which rotates the axis direction but does not change the tilt angle.

Conservation of angular momentum (Noether’s theorem)

Two counter-rotating motions with equal amplitudes produce zero net angular momentum transfer. This is required by Noether’s theorem in a system with rotational symmetry. The analogy is precise: a linearly polarized electromagnetic wave is the superposition of left- and right-circular polarization with exactly equal amplitudes. The equality is not a coincidence — it is mandated by the symmetry. If the amplitudes were unequal, the wave would be elliptically polarized, carrying net angular momentum, which would violate the symmetry.

Applied to obliquity: the two component motions that produce the observed ~41k-year obliquity cycle correspond in standard physics to Berger’s s₃ + k beat — the f_axial − f_nodal = f_obliquity identity (13 − 5 = 8 in the model’s H/n indices). Fibonacci closure within the H/n system means the equivalent identity 3 + 5 = 8 also holds (the model’s H/3 inclination precession and H/5 ecliptic precession sum to H/8 in frequency); both forms describe the same cycle. Having equal amplitudes, these motions conserve angular momentum about the relevant axis. An unbalanced system (e.g., 0.684° + some other value) would imply a net angular momentum flux — requiring an external source or sink.

Energy equipartition

The equipartition theorem states that in equilibrium, energy distributes equally among equivalent degrees of freedom. For two oscillatory modes with equivalent effective stiffness, this directly predicts equal amplitudes. While planetary systems are not thermal systems, the long-term secular evolution of orbital elements can exhibit equipartition-like behaviour over millions of years.

The virial theorem

In gravitationally bound systems at equilibrium, energies maintain fixed ratios dictated by the force law. For harmonic oscillators, the virial theorem requires that average kinetic energy equals average potential energy — a 1:1 balance. Equal amplitudes between two equivalent modes is consistent with this equilibrium condition. Crucially, balanced energy distributions are not fine-tuned states but attractors — systems that don’t satisfy the virial relation evolve toward states that do.

Precedents in nature

Equal-amplitude decompositions reflecting underlying symmetries appear throughout physics:

In every case, equal amplitudes reflect a symmetry of the underlying system. Unequal amplitudes indicate broken symmetry.

The parsimony argument (formalized)

Beyond Occam’s Razor as a philosophical principle, Bayesian model selection quantifies the advantage: a model with one amplitude parameter is inherently preferred over a model with two independent parameters, unless data demand the extra complexity. A 2022 paper (Jiang et al., J. Royal Society Interface) demonstrated that models with structural symmetries receive additional preference in model selection — beyond mere parameter counting — because symmetries encode conservation laws and fundamental physical properties.

The equal-amplitude model is not just simpler (one parameter instead of two). It implies a symmetry — the two oscillation modes are dynamically equivalent. This symmetry could reflect angular momentum conservation, energy equipartition, or the coupled-oscillator identity of the two precession mechanisms. An unbalanced model would require explaining why the symmetry is broken.

References:

• Goldstein, H., Poole, C. & Safko, J. (2002). Classical Mechanics, 3rd ed. — Normal modes of coupled oscillators.
• Noether, E. (1918). “Invariante Variationsprobleme.” — Symmetry and conservation laws.
• Laskar, J., Joutel, F. & Robutel, P. (1993). “Stabilization of the Earth’s obliquity by the Moon.” Nature, 361, 615–617.
• Jiang, B. et al. (2022). “Occam’s razor gets a new edge: the use of symmetries in model selection.” J. R. Soc. Interface, 19, 20220324.

12. Milankovitch Beat Frequency Structure

The five Milankovitch-type cycles — apsidal precession, nodal regression, obliquity, axial precession, and climatic precession — are not independent. Standard orbital mechanics derives two of them as beat frequencies of the others. The model’s identification of all five as H/n, where n ∈ 16, reveals that these beat frequency relationships are Fibonacci identities.

The standard physics

Vervoort et al. (2022) present the two fundamental beat frequency equations of Milankovitch theory:

Obliquity period:            1/P_axial  −  1/P_nodal   =  1/P_obliquity
Climatic precession period:  1/P_axial  +  1/P_apsidal  =  1/P_climatic

The first arises because obliquity measures the angle between Earth’s spin axis and orbital-plane normal — both precess in the same direction but at different rates, producing a beat. The second arises because axial precession moves the equinox backward while apsidal precession moves perihelion forward — opposite directions, so their frequencies add.

Fibonacci closure

With all five cycles expressed as H/n, the beat frequency of any two cycles H/a and H/b is:

1/(H/a) − 1/(H/b) = (a − b)/H = 1/(H/(a−b))

This means the result is again an H/n cycle whenever the index difference is meaningful. The Fibonacci subtraction property — each number equals the difference of the next two — guarantees exactly this. The full set of relationships:

Only two of these are physically independent (the first two); the rest follow algebraically. The Fibonacci structure makes the system self-consistent: any two of the five cycles determine the other three.

Comparison with standard values

Standard values: axial precession from IAU 2006; obliquity from Laskar et al. (2004); Inclination Precession from the total secular perturbation rate (~11.6″/yr including GR); Ecliptic Precession from the s₃ eigenfrequency (~18.85″/yr); climatic precession from the combined axial + apsidal rate.

Dual identification of H/3 and H/5

In the model’s Fibonacci Laws framework, H/3 = ~111,772 yr is the Inclination Precession period (Earth’s orbital plane completing one full oscillation cycle) and H/5 = 67,063 yr is the Ecliptic Precession period (the precession of Earth’s orbital plane around the invariable plane), which also equals Jupiter’s perihelion precession period. The beat frequency analysis reveals these same values also correspond to Earth-centric Milankovitch parameters:

• H/3 serves as Earth’s apsidal precession in the climatic precession equation
• H/5 serves as the nodal regression of Earth’s orbit in the obliquity equation

This dual identification is not contradictory — it suggests the solar system’s secular frequencies are organized such that Earth’s precession rates, Jupiter’s precession rate, and the inclination precession rate all fall on the same Fibonacci ladder defined by H.

Why this is significant

Economy of parameters: Standard astronomy requires 5 independent measurements to characterize the 5 Milankovitch cycles — each period emerges from different physics (tidal torques, planetary perturbations, spin-orbit coupling). The H framework needs just one number (H = 335,317) plus the Fibonacci index assignments 16. If H is chosen to match any single cycle perfectly, the other four become predictions — and they all land within 0.3–2.8%.

The closure is non-trivial: The beat frequency relationships are not optional — they are required by physics. Obliquity must equal the beat of axial and nodal precession; climatic precession must equal the sum of axial and apsidal rates. If the indices were arbitrary — say 17 — these physical equations would not close: 14 − 7 ≠ 10, and 14 + 3 ≠ 17. The Fibonacci subtraction property (each number = difference of the next two) is exactly the condition needed for the physical beat equations to close within the H/n system. Non-Fibonacci indices would fail.

Secular frequencies ≠ orbital periods: Fibonacci ratios in orbital periods are well-documented (Pletser 2019 , Aschwanden 2017 ) and understood via KAM theory — they arise from mean-motion resonances during the disk phase. But Milankovitch cycles are secular frequencies, not mean-motion. They depend on the full N-body gravitational coupling: all 8 planetary masses, all semi-major axes, the Moon’s stabilizing torque, even GR corrections. No known mechanism organizes secular frequencies into Fibonacci patterns. Finding them there extends the Fibonacci structure into a completely different dynamical regime.

Three orders of magnitude: The pattern spans from H/16 = ~20,957 yr (climatic precession) through H/3 = ~111,772 yr (apsidal precession) to 13H = 4,359,121 yr (resonant libration) — a factor of 208 in timescale — all organized by the same constant H and Fibonacci arithmetic.

Fibonacci multiples: deep-time cycles

The Fibonacci ladder extends beyond H itself. Fibonacci multiples of H match established deep-time geological cycles:

3H ≈ 1 Myr: The beat between Mercury’s and Jupiter’s apsidal eigenfrequencies (g₁ = 5.579″/yr and g₅ = 4.258″/yr in La2004 ) produces an eccentricity modulation with period g₁−g₅ ≈ 980,000 yr. This ~1 Myr cycle appears in geological records as a modulation of the short eccentricity signal. 3H = 1,005,951 yr matches it within 2.6%.

13H ≈ 4.4 Myr: The ~4.5 Myr cycle arises from the resonant argument θ = 2(g₄−g₃) − (s₄−s₃), a nonlinear coupling between the eccentricity and inclination systems of Earth and Mars. Unlike the simpler beat-frequency cycles, this is a libration period — the timescale on which the resonant angle oscillates rather than circulating. Boulila et al. (2020, Palaeogeography) measured it at ~4.5 Myr (range 3.7–4.8 Myr) in Mesozoic–Cenozoic sedimentary records. 13H = 4,359,121 yr falls within this range at ~3.1% from the central estimate.

Both 3 and 13 are Fibonacci numbers, extending the pattern from sub-divisions (H/3, H/5, H/8, H/13) to multiples. The 5H and 8H multiples do not match known present-epoch cycles — the pattern is selective, not universal.

Eigenfrequency convergence at H/3 and H/5

The Laskar (La2004) secular solution decomposes planetary eccentricity and inclination into eight eigenfrequencies each (g₁–g₈ for eccentricity, s₁–s₈ for inclination). Remarkably, multiple independent eigenfrequency combinations converge on the same H/n values:

H/3 = ~111,772 yr matches three independent combinations:

H/5 = 67,063 yr matches two independent combinations:

The ~111 kyr region is surprisingly crowded: three physically distinct cycles (apsidal precession, an eccentricity beat, and an inclination beat) all cluster within ±1.3% of H/3. Similarly, H/5 matches both a single-planet eigenfrequency (s₃) and an inter-planet beat (|s₂−s₆|). This convergence is not required by any known theory — the eigenfrequencies depend on all planetary masses and semi-major axes, and there is no reason a priori for combinations involving Mercury, Venus, Earth, and Saturn to produce the same period.

References:

• Vervoort et al. 2022, “System Architecture and Planetary Obliquity,” The Astronomical Journal, 164, 130
• Laskar et al. 2004, “A long-term numerical solution for the insolation quantities of the Earth,” A&A, 428, 261–285
• Boulila et al. 2020, “Toward a robust and consistent middle Eocene astronomical timescale,” Palaeogeography, 549, 109702

13. Saturn’s Ecliptic-Retrograde Perihelion Precession: A New Explanation

The Holistic Universe Model offers a new theoretical framework for a clearly observed but poorly explained phenomenon: Saturn’s longitude of perihelion moves retrograde in the ecliptic frame (opposite to orbital motion) at the current epoch, despite secular perturbation theory predicting prograde precession for all planets.

The observation

JPL’s WebGeoCalc  tool — which computes geometric quantities directly from the SPICE ephemeris kernels — shows that Saturn’s longitude of perihelion (ϖ = Ω + ω) is decreasing over the period 1900–2000 AD. The two components move in opposite directions:

This is confirmed by the JPL Keplerian elements table (Standish & Williams 1992 ), which gives Saturn dϖ/dt = -0.419 deg/century in the 1800–2050 fit — the only major planet (along with Neptune at -0.322 deg/century) with a negative rate.

Where the disagreement lies: ω, not Ω. Decomposing ϖ into its two components using the JPL Standish tables reveals that the ascending node is uncontested — both fit intervals agree it is retrograde:

The ascending node barely changes between the two intervals. The argument of periapse completely reverses — from retrograde in the short interval to strongly prograde in the long interval. This is because the Great Inequality affects the eccentricity eigenfrequencies (g-type, which govern ω) but not the inclination eigenfrequencies (s-type, which govern Ω). The entire debate about Saturn’s perihelion direction reduces to a single question: is ω permanently retrograde (this model) or only transiently retrograde due to the Great Inequality (standard theory)?

The standard explanation: the Great Inequality

Standard celestial mechanics attributes this retrograde observation to a transient phase of the Great Inequality — the ~900-year oscillation caused by the near-5:2 mean-motion resonance of Jupiter and Saturn. The history of this explanation:

Laplace’s Great Inequality theory predicts that Saturn’s longitude of perihelion has a long-term secular rate of approximately +19.5 arcsec/yr (prograde) — dominated by the g₆ eigenfrequency of +22.44 arcsec/yr (Fitzpatrick , based on Murray & Dermott). The JPL 6000-year fit (3000 BC – 3000 AD) confirms +19.50 arcsec/yr. Under this view, the current retrograde motion is a temporary phase of the oscillation — the sign should reverse at some point within the ~900-year cycle.

Why the theory requires prograde. This is not a choice but a structural constraint of the framework. Laplace-Lagrange secular theory decomposes long-term perihelion evolution into eight eigenfrequencies (g₁–g₈), computed from a coupling matrix whose elements depend on planetary masses and orbital distances. The mathematical structure of this matrix — positive diagonal elements (self-coupling), negative off-diagonal elements (planet-planet coupling) — produces eigenvalues that are all positive (prograde): g₁ = +5.59″/yr through g₈ = +0.67″/yr. There is no eigenmode in the theory that can produce permanent retrograde apsidal precession. If the observation shows retrograde, the framework must attribute it to a periodic (non-secular) perturbation — the Great Inequality is the only available mechanism. The theory drives the interpretation of the observation: Saturn is declared “really” prograde because the mathematics of the framework cannot produce any other answer.

Why the standard explanation is incomplete

However, the Great Inequality theory has limitations when applied to Saturn’s perihelion:

1. Never directly verified: No complete 900-year cycle of Saturn’s perihelion precession has ever been observed. The Great Inequality itself was derived from mean longitude (position along orbit), not from the longitude of perihelion. Its effect on ϖ operates through small correction modes (g₉, g₁₀ from Brouwer & van Woerkom 1950 ) that are described as “relatively small-amplitude.”

2. Magnitude problem: The observed retrograde rate (~-3,400 arcsec/century from WebGeoCalc) and the predicted secular rate (+1950 arcsec/century) differ in both magnitude and sign — the observation is retrograde while the theory predicts prograde, giving a total discrepancy of ~5350 arcsec/century. Attributing this entirely to the Great Inequality’s modulation of the perihelion requires an oscillation amplitude larger than what the small correction modes predict.

3. No independent confirmation: The claim that Saturn’s perihelion precession is “really” prograde at +19.5 arcsec/yr rests on the assumption that the Great Inequality’s mean-longitude theory applies equally to the longitude of perihelion. This has not been independently tested for ϖ.

The Holistic Universe Model’s explanation

The model proposes that Saturn’s ecliptic-retrograde perihelion precession is not a transient oscillation but a permanent feature of the model’s Fibonacci structure — specifically, of the H/8 triple identity in Law 6 where Saturn’s ecliptic perihelion period coincides with Earth’s obliquity cycle and Jupiter’s ICRF perihelion period:

Saturn’s perihelion precesses retrograde in the ecliptic frame with a period of H/8 = ~41,915 years.

Saturn’s anti-phase role (its cosine sign flipped relative to the other seven planets) makes it the unique “pivot” in the Fibonacci framework:

• Law 3 (Inclination Balance): Seven in-phase planets balance against Saturn alone (anti-phase) — to 99.9975%
• Law 5 (Eccentricity Balance): The same 7-vs-1 grouping independently balances eccentricities — to 99.8632%
• Law 6 (Resonance Loop): Saturn’s ecliptic-retrograde H/8 creates a closed beat-frequency triangle with Jupiter (H/5) and Earth (H/3): 3 + 5 = 8

The model’s predicted rate, including the missing advance of perihelion:

The geocentric prediction of -3,372″/cy closely matches the observed WebGeoCalc value of ~-3,400″/cy.

The 3D simulation  implements Saturn’s ecliptic-retrograde perihelion precession directly, producing a visual representation of how the seven-against-one balance operates over the full Earth Fundamental Cycle.

Comparison of the two theories

How to distinguish the two theories

The two theories make a decisive, testable prediction: the Great Inequality theory requires Saturn’s ecliptic-frame perihelion precession rate to oscillate between prograde and retrograde over the ~900-year cycle. The Holistic Universe Model predicts it to remain retrograde at approximately -3,092.0 arcsec/century indefinitely.

Over decades-to-centuries of continued high-precision ephemeris tracking (Cassini legacy data, future Saturn missions), the trend in dϖ/dt should either:

• Show curvature toward zero and eventually become prograde → supports Great Inequality theory
• Remain steady near -3,092.0 to -3,400 arcsec/century → supports the Holistic Universe Model

Implications if the ecliptic-retrograde motion is permanent. The model and standard theory disagree on ICRF direction. The model predicts Saturn is retrograde in any fixed frame (ICRF period H/21, rate ≈ -81.2″/yr): the general precession (period H/13, rate ≈ 50.245″/yr) is subtracted from the ecliptic rate (period H/8, rate ≈ -30.9″/yr), deepening the retrograde. Standard long-term secular theory predicts prograde (+19.5”/yr). Notably, Standish Table 1 (1800–2050, J2000 ecliptic) shows Saturn currently retrograde at −0.419°/cy even in the fixed J2000 frame — consistent with the model’s prediction. Standard theory attributes this to the transient Great Inequality; the model predicts the retrograde is permanent. If future observations confirm that the ecliptic rate remains retrograde indefinitely, it would indicate that the standard secular coupling matrix is incomplete: either a fundamentally different organizational principle (such as Fibonacci-structured KAM stability) governs Saturn’s long-term perihelion evolution, or the long-term secular average is itself incorrect.

High-precision ephemeris analyses

Saturn’s perihelion precession has been the subject of intensive study using the most precise planetary ephemerides available. In 2008, Pitjeva detected a small anomalous retrograde residual in Saturn’s precession — the amount left over after subtracting all known Newtonian and general-relativistic effects:

The EPM2008 anomaly (-6 mas/cy) was detected using early Cassini spacecraft ranging data (2004–2006). Iorio (2009)  showed that no standard Newtonian or Einsteinian effect could explain this retrograde residual — not planetary perturbations, solar oblateness, asteroid belt mass, trans-Neptunian objects, general relativity, or modified gravity theories (MOND, DGP braneworld). The anomaly was purely retrograde, with no known mechanism to produce it.

Later ephemerides (INPOP10a in 2011, EPM2011 in 2013) found the residual consistent with zero, suggesting it may have been an artifact of the limited early Cassini data span. Modern ephemerides (DE440, EPM2017, INPOP19a) do not report a significant anomalous residual for Saturn.

References:

• Standish, E.M. & Williams, J.G. (1992), “Keplerian Elements for Approximate Positions of the Major Planets,” JPL 
• Fitzpatrick, R., “Secular Evolution of Planetary Orbits” (Table 10.1) 
• Brouwer, D. & van Woerkom, A.J.J. (1950), “The Secular Variation of the Orbital Elements of the Principal Planets,” Astronomical Papers, 13, 81–107
• Wilson, C. (1985), “The great inequality of Jupiter and Saturn: from Kepler to Laplace,” Archive for History of Exact Sciences, 33, 15–290
• Iorio, L. (2009), “The recently determined anomalous perihelion precession of Saturn,” The Astronomical Journal, 137, 3615–3618
• Pitjeva, E.V. (2010), “EPM ephemerides and relativity,” Proc. IAU Symp. 261, 170–178
• Pitjeva, E.V. & Pitjev, N.P. (2013), “Relativistic effects and dark matter in the Solar system from observations of planets and spacecraft,” MNRAS, 432, 3431–3437
• Fienga, A. et al. (2011), “The INPOP10a planetary ephemeris and its applications in fundamental physics,” CeMDA, 111, 363–385
• Iorio, L. (2010), “The perihelion precession of Saturn, planet X/Nemesis and MOND,” The Open Astronomy Journal, 3, 1–6

Summary

← Mathematical Foundation | Predictions →