Supporting Evidence from Current Science

The Holistic Universe Model makes claims that challenge several established theories. This page collects external scientific evidence — 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 ~100,000-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,669 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 ~95,000 and ~125,000 years. But the climate record shows a single narrow peak near ~100,000 years. These don’t match.

2. The 400,000-year absence: Eccentricity’s theoretically strongest component (~400,000 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 41,000 years (obliquity-dominated) to ~100,000 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 ~100,000-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): 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 “~100,000-year” signal is actually ~111,669 years — the inclination precession period. The ~10% discrepancy between 100k and 111k 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 111k

Why the ~111k cycle only emerged at the MPT

The inclination precession cycle (H/3 = 111,669 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 ~111k 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 41,000-year obliquity pacemaker, allowing the longer ~111,000-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,770 yr) and inclination precession (111,669 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): 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 corrugated patterns at rational resonances with Saturn’s moons — another dramatic demonstration.

What this means for the model

The model’s 13:3 ratio is not numerology — it reflects the maximally stable orbital configuration 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 Holistic-Year 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–2022)

The model predicts that Earth’s Length of Day (LOD) follows a 20,938-year cycle linked to perihelion precession, with LOD currently decreasing (Earth speeding up) since the 1246 AD maximum.

What was observed

Starting in 2020, Earth unexpectedly began rotating faster:

• 2020: The 28 shortest days since atomic clocks were invented
• June 29, 2022: The shortest day ever recorded — 1.59 milliseconds under 24 hours
• 2019–2022: Average LOD shifted from +0.39 ms to −0.25 ms 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 that we are past the LOD maximum (1246 AD) and that Earth’s rotation should be gradually speeding up over the coming millennia, reaching a minimum LOD around 11,715 AD. The 2020–2022 speedup is qualitatively consistent with this prediction.

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 20,938-year LOD cycle 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 1.08 × 10⁻⁵ to 1.46 × 10⁻⁷ 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 Upcoming Decisive 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.

This is a binary test: either the rate decreased or it didn’t.

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 ~100,000-year glacial signal is actually ~111,669 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 111k?) 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 111k 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,876-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,000-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 between ~22.1° and ~24.5°, 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,002 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

10. Neptune’s Ecliptic Inclination: Observational Uncertainties

The model predicts that Neptune’s ecliptic inclination is decreasing, while JPL’s current trend figures show it increasing. This is the only planet where the model and JPL disagree on the direction of inclination change. However, several independent lines of evidence suggest that Neptune’s observed inclination trend carries significantly more uncertainty than is commonly appreciated.

Neptune’s observational arc is uniquely short

Neptune was discovered in 1846 and completed its first full orbit only in 2011 — just 165 years of observation covering barely one orbital period. For comparison, Jupiter has been tracked for over 400 years (multiple orbits), and even Uranus (discovered 1781) has completed nearly three orbits. Neptune’s orbital elements are therefore constrained by the least observational data of any planet.

The practical consequences are significant. Folkner et al. (2014) — the authors of JPL’s DE430/DE431 planetary ephemeris — note that Neptune’s position uncertainty is “several thousand kilometers.” This is because:

• No continuous spacecraft ranging: The only spacecraft to visit Neptune was Voyager 2 during a single flyby in August 1989. Unlike Mars, Jupiter, and Saturn (which benefit from ongoing orbiter missions providing continuous ranging data), Neptune’s orbit is constrained primarily by ground-based optical astrometry — inherently less precise than radar or spacecraft ranging.
• Limited stellar occultation data: While stellar occultations provide high-precision position fixes, they are rare for Neptune and cover only a few decades.
• Pre-discovery observations: Galileo observed Neptune in 1612–1613 without recognizing it as a planet. These observations extend the arc to ~400 years but with low precision.

Reference: Folkner et al. 2014, Interplanetary Network Progress Report 42-196

JPL’s inclination rate is sensitive to the fit interval

JPL publishes two sets of approximate orbital elements with different time ranges:

The rate changes by 60% between the two intervals. For well-constrained planets like Jupiter, the corresponding rates differ by less than 10%. This sensitivity suggests that Neptune’s inclination trend has not yet converged to a reliable long-term value. JPL’s own documentation warns: “The elements are not valid outside the given time-interval over which they were fit.”

The s₈ eigenmode vs the model’s 2H period

Laskar’s secular perturbation theory shows that Neptune’s inclination is dominated by the s₈ eigenfrequency (~−0.692″/yr), which drives an oscillation with a period of approximately 1.87 million years. The model instead uses a period of 2H = 670,016 years — roughly 3× shorter than the s₈ eigenmode period.

This parallels the situation with Jupiter (section 9), where the model’s H/5 = 67,002 years is also ~3× shorter than the theoretical g₅ eigenfrequency period (~305,000 years). For Jupiter, the shorter period produces better agreement with JPL’s observed inclination trend (3″/cy error vs 8.5″/cy). Whether the same holds for Neptune requires further investigation — the observational baseline is currently too short to distinguish between these two timescales.

What we observe over ~250 years is a tiny linear segment of either oscillation — equivalent to sampling less than 0.04% of even the shorter 2H cycle. Whether that segment happens to be on the “increasing” or “decreasing” portion depends critically on where Neptune currently sits in its cycle. A small change in the initial conditions — well within the observational uncertainty — could place Neptune on the other side of the oscillation’s turning point.

Ecliptic reference frame drift

All published inclination values are measured relative to the ecliptic — Earth’s orbital plane. But the ecliptic itself precesses relative to the invariable plane (the solar system’s fundamental reference) with a period of ~70,000 years. This adds an apparent secular drift to every planet’s ecliptic inclination that has nothing to do with the planet’s actual orbital evolution.

For inner planets with large inclination rates, this frame drift is a small correction. But for Neptune — whose intrinsic inclination rate is among the smallest of any planet — the ecliptic drift is a significant fraction of the total observed rate. The model works in the invariable-plane frame, which does not suffer from this artifact.

What this means for the model

The disagreement between the model and JPL on Neptune’s inclination direction is real — but it is the least well-constrained of all eight planetary inclination trends:

1. Shortest observational arc of any planet (barely one orbit)
2. No continuous spacecraft ranging (single Voyager 2 flyby)
3. 60% rate sensitivity to the choice of fit interval
4. Long-period oscillation (670,016 yr in the model, ~1.87 Myr in secular theory) sampled over just 250 years
5. Ecliptic frame drift contributes a significant fraction of the observed rate

None of this proves JPL’s trend is wrong — but it does mean that Neptune’s inclination direction is far less certain than for the other seven planets, where model and JPL agree. Rather than a “clean, decisive test,” this disagreement may require decades of additional observation to resolve definitively.

11. Obliquity Amplitude: Berger (1978) Dominant Term

The model claims both the axial tilt and inclination tilt oscillate by the same amplitude of ±0.634°, combining to produce the full obliquity range of ~22.21°–24.71°. 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.634°. 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 ~41,000-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.634°). 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°.

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.

12. Why Equal Amplitudes? The Physics of Balanced Systems

The model’s claim that both obliquity components oscillate by the same amplitude (±0.634°) 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 treats the axial tilt oscillation (driven by lunisolar torques, period ~25,770 yr) and the inclination tilt oscillation (driven by planetary perturbations, period ~111,669 yr) as two coupled counter-rotating motions. Equal amplitudes imply these two mechanisms are effectively equivalent oscillators — a physically meaningful claim about the dynamics, not just a parameter choice.

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 clockwise axial precession and counterclockwise inclination precession, having equal amplitudes, 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.

13. 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,669 yr is the Inclination Precession period (Earth’s orbital plane completing one full oscillation cycle) and H/5 = 67,002 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,008) 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,938 yr (climatic precession) through H/3 = 111,669 yr (apsidal precession) to 13H = 4,340,544 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,001,664 yr matches it within 2.2%.

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,340,544 yr falls within this range at ~1.4% 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,669 yr matches three independent combinations:

H/5 = 67,002 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

14. 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 (~-3400 arcsec/century from WebGeoCalc) differs from the predicted secular rate (+1950 arcsec/century) by a factor of ~2.5 with opposite sign — a 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 solar system’s Fibonacci cycle hierarchy (Law 1):

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

This is the only planet assigned ecliptic-retrograde apsidal precession in the model. The ecliptic-retrograde direction is what makes Saturn the unique “pivot” in the Fibonacci framework:

• Law 3 (Inclination Balance): Seven planets (phase 203°) balance against Saturn alone (phase 23°) — to 100%
• Law 5 (Eccentricity Balance): The same 7-vs-1 grouping independently balances eccentricities — to 100%
• 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,385″/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 Holistic-Year 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 -3100 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 -3100 to -3400 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 rate = −H/21 = −81.5”/yr): the general precession (H/13 ≈ 50.3”/yr) is subtracted from the ecliptic rate (−H/8 = −31.1”/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 →