Scientific Background

This page provides in-depth scientific context for readers who want to understand how the Holistic Universe Model relates to established astronomical theory. It addresses physical mechanisms, compares predictions with standard models, and acknowledges limitations and open questions.

Quick Reference

Table of Contents

1. The Fibonacci Laws
2. Physical Mechanisms
3. Comparison with Standard Precession Theory
4. The Mercury Perihelion Question
5. Eccentricity Cycles and Milankovitch Theory
6. The 100,000-Year Glacial Cycle
7. Mathematical Framework
8. Open Questions
9. References

1. The Fibonacci Laws

The quantitative core of the Holistic Universe Model is a set of six Fibonacci Laws that connect planetary precession periods, inclinations, and eccentricities through Fibonacci numbers and mass-weighted quantities. Together they predict orbital properties for all eight planets from a single timescale: the Earth Fundamental Cycle H = 335,317 years. For the full accessible treatment, see Fibonacci Laws of Planetary Motion; for the mathematical derivation, see Fibonacci Laws Derivation.

The Six Laws at a Glance

Law 1 — Fibonacci Cycle Hierarchy: Dividing H by successive Fibonacci numbers (3, 5, 8, 13, …) produces the major precession periods of the solar system. The corresponding frequencies obey the Fibonacci addition rule: 1/T₃ + 1/T₅ = 1/T₈, at every level. The ratios between consecutive periods converge toward the golden ratio.

Law 2 — Inclination Constant: Each planet’s mass-weighted inclination amplitude η = amplitude × √m, multiplied by a Fibonacci divisor d specific to that planet, equals a single universal constant: d × η = ψ. The constant ψ is derived from Earth’s parameters: ψ = d_E × amplitude_E × √m_E, with zero free parameters beyond the master cycle. All eight planets satisfy this law exactly (by construction).

Law 3 — Inclination Balance: The angular-momentum-weighted inclination oscillations of seven planets (Mercury through Neptune, excluding Saturn) balance against Saturn’s alone. The balance reaches 99.9975% (with dual-balanced eccentricities, which enter via angular momentum) — a consequence of the invariable plane’s stability. Saturn’s unique role arises because it is the only planet whose inclination oscillation phase differs from the other seven.

Law 4 — Eccentricity Amplitude Constant: A single constant K predicts all eight eccentricity oscillation amplitudes from Fibonacci divisors, mass, distance, and axial tilt: e_amp = K × sin(tilt) × √d / (√m × a^1.5). K = 3.4149 × 10⁻⁶, derived from Earth. This is the eccentricity analog of ψ (Law 2). See Law 4.

Law 5 — Eccentricity Balance: Using the same Fibonacci divisors and Saturn-vs-seven grouping as Law 3, the eccentricity-weighted contributions balance to 99.8632%. This is independent of the inclination balance — different weight formulas, same Fibonacci structure. The balance predicts Saturn’s eccentricity to ~0.27% accuracy.

Law 6 — Saturn-Jupiter-Earth Resonance: Saturn’s ecliptic-retrograde precession (H/8) creates a closed resonance loop with Jupiter (H/5) and Earth (H/3). The Fibonacci relationship 3 + 5 = 8 is the frequency-domain expression of this three-body resonance: 1/(H/3) + 1/(H/5) = 1/(H/8). This resonance is the physical mechanism that links the Fibonacci timescale (Law 1) to the orbital structure (Laws 2–5).

Scientific Status

The laws formalize a pattern that has independent support in the peer-reviewed literature. Fibonacci-related frequency ratios in planetary orbits were first documented by Molchanov (1968)  in Icarus, confirmed by Aschwanden (2018)  in ~60% of 75 solar system period ratios, and extended to exoplanets by Aschwanden & Scholkmann (2017)  in 73% of 932 planet pairs. Pletser (2019)  confirmed that orbits near Fibonacci ratios are associated with more regular, less inclined, and more circular configurations.

The theoretical explanation comes from the KAM theorem (Kolmogorov 1954, Arnold 1963, Moser 1962): in perturbed dynamical systems, orbits whose frequency ratios are “most irrational” — closest to the golden ratio, toward which Fibonacci ratios converge — are maximally stable. Greene and Mackay (1979) confirmed computationally that the golden invariant torus is the last to break. Morbidelli and Giorgilli (1995) demonstrated super-exponential stability near golden-ratio frequency ratios in the asteroid belt.

What the Holistic Model adds is the quantitative framework: specific laws that predict numerical values for all eight planets’ inclinations and eccentricities. A significance analysis over the 4 empirical tests (direct joint permutation test — model-independent, no distributional assumptions, correlation baked into the joint null by construction) yields a combined p-value spanning 1.5 × 10⁻⁴ (permutation null, conservative) to 1.0 × 10⁻⁶ (log-uniform Monte Carlo over 9 tests, more powerful) — equivalently 3.62–4.75σ across the three null distributions, comfortably above the conventional 3σ “evidence” threshold but short of the particle-physics 5σ “discovery” threshold. Whether this reflects a deep physical principle or an elaborate numerical coincidence is the central question this document examines.

2. Physical Mechanisms

The Holistic Universe Model’s two counter-rotating motions correspond directly to two well-established astronomical phenomena: axial precession and apsidal precession. These are not invented by the model - they are standard astronomy with known physical causes.

The Two Precessions in Standard Astronomy

Key fact: These two precessions move in opposite directions. This is well-documented in the scientific literature:

“The apsidal precession direction is opposite from the axial precession, thus climatic precession cycles experienced by the planet are more rapid than the axial precession cycles.” — Global Climate Change Organization 

Axial Precession: The Physical Mechanism

Axial precession (also called “precession of the equinoxes”) is caused by gravitational torque from the Sun and Moon acting on Earth’s equatorial bulge:

The physics:

• Earth is not a perfect sphere - it bulges at the equator (oblateness J₂ ≈ 0.00108)
• The equatorial diameter is ~43 km larger than the polar diameter
• The Sun and Moon exert differential gravitational pull on this bulge
• This creates a torque perpendicular to Earth’s rotation axis
• The torque causes the rotation axis to precess (wobble like a spinning top)

Direction: The equinoxes drift westward along the ecliptic at ~50.3 arcseconds per year. When viewed from above the North Pole, the celestial pole traces a clockwise circle.

Period: ~25,771 years currently (varies slightly over time)

Note on values: The current measured value ~25,771 years (IAU) is below the model’s mean of ~25,794 years (335,317 / 13 = 25,793.62 years) and still decreasing. The model predicts this trend will eventually reverse, with the period increasing back toward the mean (see Predictions). Throughout this document, ~25,771 years refers to the current measured value; ~25,794 years refers to the model’s mean value.

Key references:

• Capitaine, N., Wallace, P.T., & Chapront, J. (2003). “Expressions for IAU 2000 precession quantities.” A&A, 412, 567-586.
• Britannica: Precession of the Equinoxes 

Apsidal Precession: The Physical Mechanism

Apsidal precession (also called “perihelion precession”) is caused by gravitational perturbations from other planets:

The physics:

• Each planet’s gravitational pull slightly deflects Earth’s orbit
• These perturbations accumulate over time
• The net effect rotates the entire orbital ellipse around the Sun
• Jupiter contributes the most (~60%), followed by Venus and Saturn

Direction: Earth’s perihelion advances in a prograde direction (same as orbital motion). When viewed from above the North Pole, this is counter-clockwise.

Period: ~112,000 years for Earth’s ellipse to complete one full rotation relative to the fixed stars.

Calculation method (Gauss): Treat other planets as uniform concentric rings centered on the Sun, with mass equal to planetary mass and radius equal to mean orbital distance. This averages the gravitational interactions over complete orbits.

Key references:

• University of Texas: Perihelion Precession of Planets 
• NASA: Earth’s Orbital Precession 

Why Opposite Directions?

The opposite directions arise from different physical causes:

This is not a coincidence or assumption - it’s a consequence of the underlying physics.

The Combined Effect: Climatic Precession

When axial and apsidal precession combine, they produce climatic precession - the cycle that determines when Earth is closest to the Sun relative to the seasons:

Climatic precession period = 1 / (1/T_axial + 1/T_apsidal)
                           = 1 / (1/25,771 + 1/111,717)
                           ≈ 20,940 years

The formula uses addition (not subtraction) because the precessions move in opposite directions, so they “meet” more frequently.

This ~21k-year cycle is one of the Milankovitch cycles  that influence Earth’s climate.

The Model’s Representation

The Holistic Universe Model represents these same physical phenomena using a different mathematical framework:

Important: The model does not invent new motions or claim different physics. It provides an alternative mathematical representation of the same observable phenomena - similar to how both geocentric and heliocentric coordinates can accurately describe planetary positions.

3. Comparison with Standard Precession Theory

What Standard Theory Predicts

The IAU 2006 precession model (Capitaine et al. 2003) provides high-precision predictions:

Vondrák, Capitaine & Wallace (2011) extended the IAU 2006 precession expressions from a few centuries to ±200,000 years — the same timescale over which the Holistic model operates. Their long-term expressions use Fourier-type series fitted to numerical integrations (Mercury 6 package with Laskar 1993 solutions) and achieve accuracy comparable to IAU 2006 near J2000, a few arcseconds over historical timescales, and a few tenths of a degree at the ±200,000-year endpoints. This makes Vondrák et al. (2011) the most direct comparison standard for evaluating the Holistic model’s long-term precession predictions.

Where They Agree

For periods of ±2,000 years around the present, the model closely matches established theory:

• Obliquity values: Within ±0.01° of Laskar (1993) and Chapront et al. (2002)
• Longitude of perihelion: Matches Meeus (1998) within ±0.1°
• Precession rate: Matches IAU within 0.01%

Where They Diverge

For longer timescales, predictions differ:

Comparison with JPL DE440/441 Ephemeris

The JPL Development Ephemeris (DE440/441) is the gold standard for solar system dynamics, achieving sub-arcsecond accuracy for inner planets over centuries. A fair evaluation of the model requires direct comparison.

About DE440/441:

• Published: Park et al. 2021 
• Time span: DE440 covers 1550-2650 AD; DE441 extends to ±13,000 years
• Accuracy: ~0.1 mas (milliarcseconds) for inner planets over centuries
• Method: Full numerical integration with GR corrections
• Data sources: Planetary radar, spacecraft ranging, VLBI, optical observations

Orbital Element Comparison (J2000 Epoch)

Assessment: The model matches DE440 at J2000 because J2000 values were used as inputs during calibration. This match is expected and does not validate the model.

Obliquity Predictions (Model vs La2004)

Assessment: The model agrees with La2004 to within ~0.001° from -1000 BC through 5,000 AD — essentially exact over this 6 kyr window. Agreement remains within 0.1° from 7,000 AD through 12,000 AD. The largest near-term discrepancy (~0.37°) occurs at 10,000 BC, where La2004’s long-term envelope is somewhat lower than the model’s. By 20,000 AD the model and La2004 diverge by ~0.29°, with La2004 oscillating back upward while the model continues its bounded H/8 cycle.

Longitude of Perihelion (Model vs Meeus/DE440)

Assessment: Good agreement across the range where Meeus’s polynomial formula is valid (~±1000 years from J2000). At J2000 the model uses the observed value (102.947°, matching DE440), while Meeus’s polynomial gives 102.937° — a +0.01° offset that propagates through the polynomial’s projections.

Eccentricity Predictions (Key Divergence)

This is where the model differs most significantly from standard theory:

Assessment: Major divergence beyond ~10,000 years. The model predicts a minimum eccentricity of ~0.0140 at 11,725 AD; Laskar predicts continued decrease to 0.00263 near 27,000 AD. This is the model’s primary differentiating prediction.

How to Verify These Comparisons

Anyone can verify the model’s predictions against JPL data:

1. JPL Horizons (ssd.jpl.nasa.gov/horizons ):

   • Query Earth’s orbital elements for any date within DE440/441 range
   • Compare eccentricity, obliquity, longitude of perihelion

2. Model Calculator:

   • Use the formulas at Formulas
   • Enter any year and calculate the model’s predictions

3. 3D Simulation:

   • The Interactive 3D Simulation displays all values in real-time

Limitations of This Comparison

DE440/441 limitations:

• Based on ~100 years of precise tracking data
• Long-term extrapolations (>centuries) are modeled, not measured
• Chaotic behavior limits predictability beyond ~50 Myr (Laskar et al. 2011)

Model limitations:

• 6 free parameters, all governing the Earth simulation; the planet configuration is uniquely determined by mirror-symmetry constraints (no additional degrees of freedom)
• No physical derivation from celestial mechanics
• Fibonacci ratios are assumed, not derived

Important: For timescales beyond a few hundred years (the high-precision observation era), neither the model nor DE440/441 can be directly verified against observations. Both extrapolate from the same modern precision data; DE441’s nominal ±13,000-year validity comes from numerical integration calibrated against that modern data, not from independent verification across its full span. The methodological difference: DE441 uses full N-body integration with GR; the model uses a parameterized formula.

4. The Mercury Perihelion Question

This section examines one of the most debated aspects of the Holistic Universe Model: the alternative explanation for Mercury’s ~43 arcsecond/century perihelion precession “anomaly.”

Historical Context

Mercury’s perihelion precession was a crucial test for gravitational theory:

Timeline:

• 1859: Urbain Le Verrier identifies a ~38″/century discrepancy between observed Mercury precession and Newtonian prediction (using telescope observations of Mercury transits)
• 1882: Simon Newcomb refines the value to ~43″/century
• 1915: Einstein’s General Relativity predicts ~43″/century from space-time curvature, derived from the standard formula Δϖ_GR = 6πGM/(ac²(1−e²)) per orbit
• 1960s onward: Radar ranging from Earth improves measurement precision
• 2011-2015: MESSENGER spacecraft orbits Mercury, enabling radio ranging measurements
• 2017: Park et al. publish MESSENGER analysis: total precession = 575.3100 ± 0.0015″/century

The Measurement Breakdown

Key point: Both the ~575″ and ~5,604″ values are measured in the ecliptic plane — the difference is the reference direction. The ~575″ value is relative to fixed stars (ICRF), the inertial frame defined by distant quasars. The ~5,604″ value is relative to the moving vernal equinox, which drifts backward at ~5,028.8″/century due to Earth’s axial precession. This equinox-based value is what was historically measured before ICRF corrections existed, and it is what is actually experienced from Earth’s reference frame.

The Geocentric Total: ~5,600 Is an Approximation

The commonly cited “~5,600″/century” geocentric total originates from Clemence (1947) , who used Newcomb’s 19th-century equinox precession rate of 5,025.645″/century. The full Clemence breakdown (Berche & Medina, 2024 , Table 2):

However, Newcomb’s equinox precession (5,025.645″) has since been updated. The IAU 2006 precession model (P03)  gives a rate of 5,028.796″/century (Lieske 1976: 5,029.097″). With the modern rate, the geocentric total should be ~5,604″ rather than ~5,600″ — the literature simply hasn’t updated this rounded figure.

Independent N-body computations confirm this: Smulsky (2011) , working at the Institute of Earth’s Cryosphere (Siberian Branch, Russian Academy of Sciences), computed Mercury’s geocentric perihelion rotation using a fundamentally different approach from classical perturbation theory. His Galactica program — a Fortran-based N-body numerical integrator — simultaneously solves the gravitational equations for all solar system bodies treated as point masses, with integration spans covering up to 100 million years. Rather than using analytical approximations, Galactica performs direct numerical integration of the full equations of motion.

Smulsky’s analysis also introduces a compound model of the Sun’s rotation, distributing solar mass symmetrically across bodies in the equatorial plane to simulate solar oblateness and rotational effects. He argues this compound solar rotation accounts for the ~53″/century surplus over Newtonian planetary perturbations (~530″) — offering an alternative to the general relativistic explanation (~43″). His computed geocentric values are epoch-dependent:

The convergence is notable: three independent approaches — Smulsky’s N-body integration (5,601.9″), Berche & Medina’s analytical review (5,599.7″), and the Holistic Universe Model (5,598.26″) — all arrive at geocentric totals near 5,598–5,602″ at epoch J2000, despite using different methodologies, different software, and different theoretical frameworks for explaining the anomalous component. The model’s value is closest to Berche & Medina’s analytical result.

Furthermore, Smulsky’s results show the geocentric total decreasing between epochs (5,602.9″ at 1950 → 5,601.9″ at 2000, a drop of ~1″ over 50 years). This epoch-dependence aligns qualitatively with the model’s prediction of a systematic decrease over time due to Earth’s precession cycles:

The geocentric values are what is actually measured on Earth. The standard theory predicts these values remain constant (~5,604″). The model predicts a decrease of ~5.0″/century (averaged over 1800–2100) — a testable difference.

The Standard Explanation (General Relativity)

General Relativity predicts additional perihelion precession due to space-time curvature near the Sun:

Δφ = 6πGM / (c²a(1-e²)) per orbit

For Mercury: ~0.1036″ per orbit × 415.2 orbits/century ≈ 43.0″/century

This is not a free parameter - it’s calculated directly from:

• G (gravitational constant)
• M (solar mass)
• c (speed of light)
• a (Mercury’s semi-major axis)
• e (Mercury’s eccentricity)

Modern verification (Park et al. 2017 ):

• MESSENGER spacecraft orbited Mercury from March 2011 to April 2015
• Radio ranging between Earth tracking stations and MESSENGER provided precise distance measurements
• Combined with Earth’s known position, this yields Mercury’s position in ICRF coordinates
• Result: 575.3100 ± 0.0015″/century total precession
• PPN parameters: (β-1) = (-2.7 ± 3.9) × 10⁻⁵
• Range measurement precision: ~0.8 meter RMS

Historical measurement methods:

• 1859-1882 (Le Verrier, Newcomb): Telescope observations of Mercury transits across the Sun
• 1960s-2000s: Radar ranging from Earth to Mercury’s surface
• 1974-75 (Mariner 10): Two flybys provided limited gravity field data
• 2011-2015 (MESSENGER): First spacecraft to orbit Mercury, enabling unprecedented precision

The Measurement Chain: An Open Question

The 575″/century value is reported as Mercury’s precession “relative to ICRF.” But how is this actually measured?

The measurement chain:

1. Earth tracking stations ←→ Radio signals ←→ MESSENGER (orbiting Mercury)
2. Round-trip time → Distance from Earth to MESSENGER
3. Earth's position in ICRF (calculated from Earth orientation models)
4. Mercury's position = Earth's position + measured distance vector
5. Track Mercury's longitude of perihelion over years → precession rate

The critical dependency: Step 3 requires knowing Earth’s position in ICRF. This comes from Earth orientation models that account for:

• Earth’s rotation (UT1)
• Polar motion
• Precession and nutation
• Length of day variations

The open questions: Three layers of processing separate the raw measurement from the reported 575″/century:

1. Reference frame transformation: Is the reported value truly “Mercury in ICRF” or is it actually “Mercury relative to Earth, then transformed to ICRF”? If Earth’s long-period precession motions (axial ~26k years, apsidal ~~112k years) have any systematic modeling errors, these would propagate into Mercury’s calculated position.

2. Newtonian subtraction: The heliocentric 575″ is not measured directly — what is actually measured is the geocentric precession (~5,604″/century). The ~5,604″ value is relative to the moving vernal equinox, which drifts backward at ~5,028.8″/century due to Earth’s axial precession. The Newtonian contribution (~532″) is then subtracted to obtain the heliocentric rate, and the ~43″ residual is attributed to GR. Whether this subtraction fully accounts for all reference frame effects is precisely the question.

3. GR-inclusive ephemeris fit: A further methodological caveat applies to modern measurements specifically. The reported “575.31″/cy” comes from fitting a GR-inclusive ephemeris to spacecraft ranging data (Park et al. 2017  fit Mercury’s orbit jointly with all planets, 343 asteroids, and the PPN parameter β, finding β ≈ 1). This is conceptually equivalent to measuring the Newtonian baseline (~532″/cy) and adding the assumed GR contribution (~43″/cy). If BepiColombo applies the same GR-inclusive analysis, any change in the underlying perihelion advance from frame effects may be absorbed into a slightly different best-fit β, into residuals, or into the orbital baseline — rather than appearing as a clean drift in the reported total.

Why this matters for the model’s argument: The model proposes that the ~43″ “anomaly” may arise from how Earth’s reference frame motion affects the measurement. If the ICRF transformation, the Newtonian subtraction, or the GR-inclusive fit doesn’t perfectly account for Earth’s precession, a residual would appear in Mercury’s calculated precession — and would look like an “anomaly.”

This is a technical question that requires detailed analysis of the IERS (International Earth Rotation and Reference Systems Service) Earth orientation parameters and their uncertainties over long timescales.

A Broader Precedent: Reference Frame Assumptions in Cosmology

The question of whether reference frame assumptions can produce measurement artifacts is not unique to Mercury’s perihelion. A strikingly parallel debate is playing out in cosmology around the Hubble tension — the persistent ~8% discrepancy between the expansion rate measured from the early universe (H₀ ≈ 67.4 km/s/Mpc from the CMB) and the local universe (H₀ ≈ 73.2 km/s/Mpc from the distance ladder).

Local H₀ measurements convert observed redshifts to the CMB rest frame, assuming the CMB dipole (~370 km/s) is purely kinematic. However, multiple independent datasets now challenge this assumption:

• The cosmic dipole anomaly: The dipole in distant quasar and radio source counts is 2–5× larger than predicted from the CMB kinematic dipole, rejected at >5σ by multiple groups (Secrest et al. 2021 , Dam et al. 2023 , Wagenveld et al. 2025 ). This suggests the CMB rest frame may not be the correct rest frame for matter.
• H₀ anisotropy: The measured Hubble constant varies with direction on the sky at 3–4σ significance (Boubel et al. 2024 , Hu et al. 2024 ), with a dipolar pattern consistent with bulk flow contamination.
• Rest frame choice matters: Wiltshire et al. (2013)  found with decisive Bayesian evidence that the Hubble flow is more uniform in the Local Group rest frame than in the CMB frame — implying the standard CMB-frame correction may itself introduce a systematic bias.
• Tilted cosmology: Tsagas (2021–2024)  showed that in General Relativity (but not Newtonian gravity), observers moving relative to the CMB frame can measure a different deceleration parameter — meaning bulk motion can create the illusion of cosmic acceleration.

The structural parallel to Mercury’s perihelion question is direct:

This does not claim that the same physical mechanism explains both anomalies. The Mercury question involves Earth orientation models at solar-system scale; the Hubble tension involves cosmological rest frame choice. The parallel is in reasoning structure: both cases illustrate how reference frame assumptions embedded in measurement pipelines can produce apparent anomalies that are interpreted as requiring new physics — when the underlying issue may be the reference frame itself.

The cosmic dipole anomaly, at >5σ, demonstrates that reference frame questions in precision measurement science are not merely theoretical concerns. They are active, unresolved problems at the frontier of observational cosmology.

Planetary Contributions to Newtonian Precession

The ~532″/century Newtonian prediction comes from gravitational perturbations. The precise Clemence (1947) values are shown in the geocentric breakdown above. Rounded summary:

Source: These values derive from Lagrange-Laplace secular perturbation theory, originally calculated by Le Verrier and Newcomb, refined by Clemence (1947), and updated with modern ephemeris data.

Analytical vs. Numerical Methods: A Known Discrepancy

When calculating planetary contributions using first-order Laplace-Lagrange secular perturbation theory, the analytical results typically overestimate precession by ~3-4% compared to full numerical integration (as used in JPL DE440/441 ephemerides).

Why the analytical method overestimates:

1. First-order approximation only: The secular theory uses only first-order terms in planetary masses. Higher-order terms (mass², mass³, etc.) are neglected, which accumulates errors especially for Jupiter’s large mass influence.

2. Periodic terms assumed to cancel: Secular theory assumes that periodic perturbations (short-term oscillations) perfectly average to zero over complete orbits. In reality, they don’t fully cancel—some “residual” effects remain that only numerical integration captures.

3. Limited eccentricity/inclination corrections: The classical formulas assume nearly circular, coplanar orbits. Mercury has the highest eccentricity (0.206) and inclination (7°) of the inner planets, making these corrections more significant.

4. No indirect (cascading) effects: When Venus perturbs Mercury, it also slightly shifts Earth’s position, which then affects Mercury differently. These second-order cascading effects require full N-body integration to model correctly.

Individual planet accuracy: For individual planetary contributions, the analytical method shows 5-50% deviations from numerical integration, though these errors partially compensate in the total sum. The full per-planet comparison is documented in the project’s technical notes  (project documentation, not peer-reviewed).

Important note on circularity: The canonical ~532″ value has a complex history:

• Historically (Le Verrier through Clemence): Calculated independently using pure Newtonian mechanics
• Modern ephemerides (JPL DE440/441): Include GR effects in numerical integration, so the “Newtonian contribution” is often derived by subtracting the theoretical GR value (~43″) from the total

This creates a potential circularity: if ~532″ = 575″ (observed) - 43″ (GR prediction), then using ~532″ to “confirm” GR involves assuming GR is correct. This is one of Křížek’s critiques of the Mercury perihelion test.

Historical context addressing circularity: The original discovery of the Mercury anomaly by Le Verrier (1859) was entirely pre-relativistic. Le Verrier calculated the planetary perturbations using Newtonian mechanics only and found a discrepancy of ~38″/century. This calculation predated Einstein’s General Relativity by 56 years. The subsequent refinement to ~43″ by Newcomb (1882) was also purely Newtonian. Thus, the existence of an anomaly was established independently of GR. The circularity concern applies primarily to modern precision values where GR-based ephemerides are used.

Uncertainties and Academic Critiques

Several researchers have questioned aspects of the Mercury perihelion test:

Křížek’s critique (Křížek & Somer, Mathematical Aspects of Paradoxes in Cosmology, 2023):

• Notes the lack of explicit error bars on the Newtonian ~531″ calculation
• Points out that the ~43″ result comes from subtracting two large, uncertain numbers
• Argues this is mathematically “ill-conditioned” (small errors in inputs produce large errors in output)
• Calculates that the “missing” precession corresponds to only ~96 km/year movement

The 96 km/year calculation (Křížek 2015 , Křížek 2019 ):

The ~43″/century discrepancy translates to a surprisingly small physical distance:

Mercury's perihelion distance: 46,000,000 km (0.307 AU)
Circumference at perihelion: 2π × 46,000,000 = 289,026,524 km
Full circle: 360° = 1,296,000 arcseconds

Arc length per arcsecond: 289,026,524 / 1,296,000 = 223.04 km

43 arcseconds/century = 223.04 × 43 = 9,591 km/century
                      = 95.9 km/year ≈ 96 km/year

This means the entire GR “correction” amounts to Mercury’s perihelion shifting by 96 km per year.

Note on measurement precision: While 96 km/year may seem small, modern astrometry easily achieves this precision. MESSENGER achieved ~0.8 meter RMS range precision, meaning 96 km is approximately 120,000× larger than the measurement uncertainty. The smallness of 96 km relative to astronomical distances does not make it difficult to measure.

For comparison (Corda 2023 ): The Solar System barycenter (center of mass) shifts by approximately 1,000 km per day due to planetary motions - much larger than Mercury’s 96 km/year perihelion shift.

Caveat on this comparison: The barycenter motion is well-characterized in modern ephemerides (JPL DE series) and is explicitly corrected for in coordinate transformations. The comparison illustrates the scale of motions that must be accounted for, but does not directly demonstrate that the 96 km/year is a measurement artifact - the standard position is that barycentric corrections are accurately handled. The comparison shows that precision astrometry involves accounting for motions of this magnitude, making a 96 km/year residual significant if it exists.

The model’s explanation: The 96 km/year can be explained by Earth’s reference frame motion rather than relativistic effects.

The model’s two motions have these physical distances per year:

Earth around its wobble center:
Location: At Earth (~1 AU from Sun)
Radius: ~202,880.73 km (from Earth's center)
Circumference: 2π × 202,880.73 = 1,274,737.21 km
Period: ~~25,794 years
Movement: 1,274,737.21 / ~25,794 = 49 km/year (clockwise)

Earth's perihelion point around Sun:
Location: Around the Sun (at radius 0.015386 AU from Sun)
Radius: ~2,301,680.63 km
Circumference: 2π × 2,301,680.63 = 14,461,885.94 km
Period: ~~111,772 years
Movement: 14,461,885.94 / ~111,772 = 129 km/year (counter-clockwise)

Key distinction: The ~49 km/year is Earth’s actual physical motion at its location. The ~129 km/year is the motion of Earth’s perihelion point around the Sun, which affects observations made from Earth at 1 AU distance.

Why it’s not simple arithmetic: The 96 km/year is not simply 49 + 129 or 129 − 49. The relationship involves:

1. Angular projection: Earth’s motion must be projected onto the direction of Mercury’s perihelion
2. Distance ratio: The effect at Mercury’s orbit (0.307 AU) differs from the effect at Earth’s orbit (1 AU)
3. Phase relationship: The counter-rotating motions create interference patterns over time

The Interactive 3D Simulation computes these geometric relationships directly. The apparentRaFromPdA function transforms Mercury’s true perihelion position to its apparent position as seen from Earth’s moving reference frame, producing the predicted decrease — from ~5,598.26″ (2000) toward ~5,592.85″ (2100) in the geocentric frame that is actually observed on Earth.

Analytical Formulas: The model derives the ~43″ fluctuation both empirically from the 3D simulation and analytically with closed-form formulas. These formulas reproduce the simulation results, confirming:

• The vector geometry of Earth’s two precession motions
• The projection onto Mercury’s orbital plane
• The time-varying phase relationship between the cycles

These formulas demonstrate that the ~43″/century fluctuation at year 1900 emerges mathematically from the configured cycle periods — it is not an empirical fit but a consequence of the model’s fundamental structure.

Verification: The unified ~2,400-term predictive formula matches the 3D simulation output with R² = 0.999999 for Mercury across the full 335,317-year cycle, confirming that the geometric transformations in apparentRaFromPdA are correctly computing the combined effect of Earth’s two precession motions. For the complete derivation and coefficient breakdown, see Formula Derivation.

Simulation verification over the full cycle:

The 3D simulation calculates Mercury’s apparent precession across the complete 335,317-year Earth Fundamental Cycle:

The fluctuation ranges from -180″ to +202″/century over the full cycle. At year 1900 — the era when Newcomb (1882) established the canonical ~43″ figure later cited as evidence for Einstein’s General Relativity (1915) — the model’s anomaly is ~43.01″, essentially matching it. By year 2000 the model’s anomaly has decreased to ~38.02″ — the model’s distinctive prediction of a variable (rather than constant) anomaly. The pattern is non-sinusoidal because the fluctuation results from the interference of multiple periodic components (see Formula Derivation for the 106-term breakdown).

Baseline comparison: The model’s baseline is close to the standard literature value (~532″):

Technical note on the baseline: The model’s baseline is determined by Mercury’s perihelion precession period:

Mercury perihelion period: 243,867 years (H × 8/11 Fibonacci fraction)
Baseline precession: 1,296,000″ / 243,867 × 100 = 531.4″/century

The Mercury period (243,867 years = 335,317 × 8/11) follows from the Fibonacci-fraction pattern discovered in the solar system’s orbital periods.

Impact on the model’s claim: This discrepancy does NOT invalidate the testable prediction. The prediction concerns whether the observed geocentric precession changes over time, not the absolute baseline. If observations show:

• Constant ~575 + 5,028.8 = ~5,604″/century (geocentric) → GR is supported regardless of baseline
• Decreasing from ~5,598.26″ (2000) toward ~5,592.85″/century (2100) → model’s interpretation gains support regardless of baseline

The heliocentric values (~569.46″ at 2000, ~564.05″ at 2100) are derivatives — what is actually measured on Earth is the geocentric total (~5,598.26″).

However, even this small ~0.6″ baseline discrepancy warrants explanation. A proper reconciliation with standard ephemerides would strengthen the argument.

The mathematical relationship:

The ~96 km/year and the ~43″/century represent the same physical quantity at different scales:

At Mercury's perihelion distance (0.307 AU = 45.93 million km):

43 arcsec/century = 43 × (π/180) × (1/3600) radians/century
                  = 2.085 × 10⁻⁴ radians/century

Arc length = radius × angle
           = 45,930,000 km × 2.085 × 10⁻⁴
           = 9,574 km/century
           = 95.7 km/year ≈ 96 km/year

This confirms: The ~96 km/year is simply the physical distance corresponding to the angular anomaly (~43″/century) when measured at Mercury’s perihelion distance from the Sun.

Alternative GR formula (velocity-based, from Vankov 2010 ):

The standard GR precession can also be calculated using orbital velocity instead of GM/a:

Δφ = 6π(v/c)² / (1-e²) × (orbits per century)

Where:
  v = 47.87 km/s (Mercury's mean orbital velocity)
  c = 299,792.458 km/s (speed of light)
  e = 0.2056 (Mercury's eccentricity)
  orbits/century = (days/year × 100) / Mercury's orbital period
                 = 36525 / 87.969 = 415.2 (using Julian century = 36525 days)

Calculation:
  (v/c)² = (47.87/299792.458)² = 2.55 × 10⁻⁸
  6π(v/c)² = 4.81 × 10⁻⁷ radians per orbit
  ÷ (1-e²) = 4.81 × 10⁻⁷ / 0.9577 = 5.02 × 10⁻⁷ radians per orbit
  × (180/π) × 3600 = 0.1036 arcseconds per orbit
  × 415.2 orbits = 43.0 arcseconds/century ✓

This formula produces the same ~43″/century result.

Historical note: This formula was first published by Paul Gerber in 1898 (Gerber 1898 ) - 17 years before Einstein’s General Relativity. Gerber assumed that gravity propagates at the speed of light, arriving at the identical mathematical result.

Important context: The mainstream physics community considers Gerber’s derivation flawed - his assumptions lacked proper physical justification, and Max von Laue argued that Gerber’s potential does not produce the correct equations of motion when consistently applied. The consensus view is that Gerber’s correct result was a mathematical coincidence rather than a genuine theoretical insight. Einstein’s derivation, based on the geometric structure of spacetime, is considered the physically sound explanation.

Why this is mentioned: Despite these criticisms, Gerber’s work demonstrates that the same numerical formula (~43″/century for Mercury) can emerge from different theoretical frameworks. This historical fact is relevant when evaluating whether the ~43″ value uniquely confirms GR, or whether other approaches might also produce this result.

Note on alternative critiques: A paper by Nguyen (vixra:2402.0138 ) questions whether using instantaneous velocity in the GR formula can logically demonstrate spacetime curvature over an entire orbit. However, viXra is an open-access repository with no peer review, and the standard physics response is that using instantaneous rates to compute cumulative effects is precisely how calculus and differential equations work — the precession per orbit is the integral of the instantaneous precession rate over the orbital path. This critique is mentioned for completeness but is not considered a serious challenge to GR by the physics community.

The Model’s Alternative Interpretation

The Holistic Universe Model proposes the ~43″ discrepancy arises from Earth’s reference frame motion, not relativistic space-time curvature:

The argument:

1. All observations of Mercury are made from Earth
2. Earth undergoes two precession motions:
   • Axial precession: ~25,794 year clockwise cycle (Earth around its wobble center)
   • Apsidal precession: ~111,772 year counter-clockwise cycle (Earth’s perihelion point around Sun)
3. These combined motions shift Earth’s orientation relative to fixed references
4. The ~43″/century “anomaly” may reflect this observer-frame effect

Geometric mechanism (conceptual):

Imagine standing on a slowly rotating platform while trying to measure the position of a distant object:

1. The true position: Mercury’s perihelion precesses at ~531-532″/century (Newtonian) in the heliocentric frame
2. Your observation platform rotates: Earth’s observation direction shifts due to:
   • The wobble center: Earth orbits a point ~202,881 km away (clockwise, ~25,794 years)
   • The perihelion point: The reference direction shifts as Earth’s perihelion orbits the Sun (counter-clockwise, ~111,772 years)
3. The apparent position differs from true position: Because you’re rotating relative to the “fixed” background stars (quasars), your measurement of Mercury’s perihelion includes your own motion

The key insight: Your observation axis is not fixed in space. As Earth’s wobble and perihelion motions progress, your “straight ahead” direction changes relative to the distant quasars. This makes Mercury’s perihelion appear to be in a slightly different position than its true heliocentric location.

Analogy: Standing on a merry-go-round measuring the angle to a distant building. Your measurement changes not because the building moves, but because your reference frame rotates.

Important caveat on this analogy: The ICRF (International Celestial Reference Frame) is specifically designed to eliminate reference frame rotation effects by defining positions relative to distant quasars. Standard astrometry corrects for Earth’s known rotations. The model’s argument is more subtle: it proposes that the long-period components of Earth’s motion (~25,794 and ~111,772 year cycles) may not be fully captured in the standard IAU precession models used to transform coordinates to ICRF. This is a technical claim that would require detailed analysis of the IERS Earth orientation parameters to verify.

The calculation: The model’s prediction is not theoretical - it comes directly from the Interactive 3D Simulation. The apparentRaFromPdA function in the simulation calculates Mercury’s apparent perihelion position as observed from Earth by:

1. Computing the geometric angle between Mercury’s perihelion direction and Earth’s perihelion direction
2. Accounting for Earth’s position on its wobble cycle around the wobble center
3. Transforming to the apparent position as seen from Earth’s moving reference frame

This produces two outputs for each planet, visible in the 3D simulation as:

• <planet> (heliocentric): The true precession rate measured against fixed stars (ICRF) — the geometric angle between the planet’s perihelion point and the Sun. For Mercury at J2000: ~569.46″/century.
• <planet> (geocentric): The apparent precession rate as seen from Earth’s moving reference frame, computed by apparentRaFromPdA. This adds the equinox drift (~5,028.8″/century) to produce the value actually measured on Earth. For Mercury at J2000: ~5,598.26″/century.

The model predicts the geocentric value decreases from ~5,598.26″ (2000) toward ~5,592.85″ (2100), a ~5″/century drift. The calculation emerges from the configured movements in the 3D model, not from fitting parameters to match GR.

Analytical Formula for Planetary Precession Fluctuation

The fluctuation formula was derived from analysis of the 3D simulation data spanning the complete 335,317-year Earth Fundamental Cycle. The key insight is that Mercury’s fluctuation arises from three interacting movements that create frequency mixing through amplitude modulation.

For the complete formulas, planetary parameters, amplitude scaling relationships, and combination periods, see the Planetary Precession Fluctuation section in the Formulas reference.

Mercury Fluctuation Formula

A simple two-term formula only achieves R² ≈ 0.20 for Mercury across the full Earth Fundamental Cycle. The full formula uses actual observed angles (Mercury Perihelion, Earth Perihelion, Obliquity, Eccentricity, Earth Rate Deviation) from the model data and achieves R² = 1.0000. For the complete formula and coefficients, see the Mercury Fluctuation Formula section; the year-only unified ~2,400-term predictive system achieves R² = 0.999999, RMSE = 0.0974″/century.

Predictive Formulas (Time-Only Input)

A key validation of the Holistic Model is that planetary precession can be predicted without observing the planet’s perihelion position. Since planetary perihelions precess at known rates, we can calculate their positions from time alone:

$\theta_M(t) = \theta_{M,0} + \Delta\theta_M + \frac{360°}{\text{<V k="mercuryPeriPeriod" />}} \times t$

$\theta_V(t) = \theta_{V,0} + \Delta\theta_V + \frac{360°}{\text{<V k="venusPeriPeriod" />}} \times t$

Where $\Delta\theta$ are the angle corrections that align theoretical precession to the model:

Predictive Formula Accuracy (Unified ~2,400-Term System):

All planets can be predicted with >99.998% accuracy (R² > 0.99998 for all 7 planets) using only:

• Time (t)
• Earth formulas: θ_E, obliquity, eccentricity, ERD
• Precession periods and angle corrections

No observation of any planet’s perihelion position is required.

This strongly validates the model’s core claim: planetary precession “anomalies” are reference frame effects calculable entirely from Earth’s perspective. The formulas use analytical ERD (true derivative of the 25-harmonic Earth perihelion formula, RMSE ~0.0006° vs actual orbital data).

Physical Interpretation

The fluctuation represents how Earth’s changing observation direction affects the apparent position of any planet’s perihelion. For Mercury, three movements interact:

1. Earth’s effective perihelion (~20,957 years)

• Created by the combination of axial precession (~25,794 years) and true perihelion precession (~111,772 years)
• This is the COMMON component affecting observations of ALL planets from Earth

2. Earth’s true perihelion (~111,772 years)

• The orbital period of Earth’s perihelion point
• This was MISSING from simpler formulas and explains the failure at distant epochs

3. Mercury’s perihelion (243,867 years)

• Mercury’s perihelion point orbits the Sun
• Creates a double-angle (2×) pattern due to orbital symmetry

4. Frequency mixing (amplitude modulation)

• When Earth’s and Mercury’s angular rates combine, they create NEW frequencies
• The |sin(θ_E - θ_M)| term acts as an amplitude modulator
• This produces sidebands at 7,192 years and 28,237 years
• These mixing products dominate the fluctuation pattern

┌─────────────────────────────────────────────────────────────────────────┐
│                    FREQUENCY MIXING VISUALIZATION                       │
│                                                                         │
│    Three input frequencies mix to create output spectrum:               │
│                                                                         │
│    Input:                      Output (after mixing):                   │
│    ├─ ~20,957 yr (Earth eff)    ├─ 7,192 yr  (2×diff + sum) ← STRONG    │
│    ├─ ~111,772 yr (Earth true)  ├─ 19,299 yr (sum)                      │
│    └─ 243,867 yr (Mercury)     ├─ 22,928 yr (diff)                      │
│                                ├─ 28,237 yr (Mercury internal) ← STRONG │
│                                ├─ ~111,772 yr (true perihelion)         │
│                                └─ 121,933 yr (Mercury / 2)              │
│                                                                         │
│    The |sin(θ_E - θ_M)| modulator creates the sidebands                 │
│    just like AM radio signal mixing creates upper/lower sidebands       │
└─────────────────────────────────────────────────────────────────────────┘

This frequency mixing is why the simple two-term formula fails across the full Earth Fundamental Cycle — it misses the dominant 7,192-year and 28,237-year mixing products.

Venus: A Contrasting Case

Venus presents a different case: its low eccentricity (0.00678 vs Mercury’s 0.20564) makes the perihelion poorly defined, so the fluctuation is dominated by Earth Rate Deviation (ERD² × periodic) rather than geometric modulation.

This supports the reference-frame interpretation:

• Planets with well-defined perihelia (high eccentricity, e.g. Mercury) → geometric modulation dominates
• Planets with poorly-defined perihelia (low eccentricity, e.g. Venus) → Earth’s rate variations dominate

GR predicts a constant relativistic correction of ~8.6″/century for Venus using the same spacetime-curvature formula as Mercury. The model’s fluctuation is a different quantity — it measures deviation from baseline due to reference-frame effects, and Venus’s fluctuation varies by hundreds of arcseconds over the Earth Fundamental Cycle precisely because the poorly-defined perihelion primarily reflects Earth’s rate variation, not a constant relativistic effect.

For Venus formula details (R² = 0.999997, RMSE, Python implementation, group structure), see the Venus subsection in the Formulas reference.

Other Planets

The Holistic Universe Model extends beyond Mercury and Venus to predict precession fluctuations for all seven planets. Each planet’s perihelion precession period follows a Fibonacci-fraction ratio of the Earth Fundamental Cycle (H = 335,317 years), with a stable baseline plus a fluctuation that varies over the full cycle:

Key observations:

• All planets achieve R² ≈ 1.0000 — the model explains essentially all variance
• Saturn shows ecliptic-retrograde precession (negative baseline), correctly captured by the model. This ecliptic-retrograde motion is directly confirmed by JPL WebGeoCalc (~-3,400 arcsec/century); the model’s prediction is -3,372 arcsec/century (H/8), matching the observed rate. Standard theory attributes this to a transient phase of the Great Inequality (~900-year Jupiter-Saturn oscillation, Laplace 1784); the model treats it as a permanent feature. See Supporting Evidence §13 for the full comparison

See Formula Derivation for implementation details and accuracy metrics.

Questions for This Interpretation

Q1: Why hasn’t the anomaly changed since 1882?

• Le Verrier (1859): 38″/century
• Newcomb (1882): 43″/century
• Modern: 42.98″/century
• The value has been remarkably stable for 140+ years

Model response: This is a significant challenge to the model’s interpretation. The model predicts ~5.0″/century change in observed precession.

However, several factors may explain this:

1. Modern methods don’t measure the raw perihelion advance: Le Verrier (1859) and Newcomb (1882) measured Mercury’s perihelion advance directly via telescope observations of solar transits, recording the angle in the geocentric/ecliptic frame relative to the moving vernal equinox. Modern measurements (radar from the 1960s, MESSENGER 2011-2015 — Park et al. 2017 ) measure spacecraft ranges, then compute the perihelion advance via a global GR-inclusive ephemeris fit (DE440/441 ) that estimates Mercury’s orbit jointly with all planets, 343 asteroids, and the PPN parameter β. The ~575″/cy figure is reported in ICRF, derived (not directly measured) by subtracting equinox drift. Modern values therefore can’t directly falsify the geocentric drift the model predicts — the methodology assumes GR holds in the fit and reports a residual relative to that assumption.
2. Method standardization: After Newcomb’s work became the standard, subsequent measurements used similar methods and reference frames, potentially stabilizing around ~43″ regardless of actual slight variations.
3. The stability itself is the test: If the model is correct, the geocentric precession (currently ~569.46 + 5,028.8 = ~5,598.26″/century) should decrease toward ~564.05 + 5,028.8 = ~5,592.85″/century over the coming century. It is this geocentric value — what is actually measured on Earth — that the model predicts will drift. The “anomaly” (observed minus Newtonian) would decrease accordingly. This prediction is falsifiable.

Honest assessment: The remarkable stability of ~43″ since 1882 is partially confounded by this methodological shift. Modern measurements that explicitly track the raw geocentric perihelion advance — rather than computing GR residuals from a global fit — would provide a more definitive test of the model’s prediction.

Q2: If ICRF measurements already correct for Earth’s motion, how can the model claim a residual Earth-frame bias?

The challenge:

• ICRF is defined by distant quasars — essentially fixed in space
• Standard reductions (precession, nutation, polar motion) translate observations from the moving Earth frame into ICRF, so Earth’s motion should already be removed
• Therefore the model’s “reference-frame effect” appears redundant

Model response: The model does NOT claim there’s an unknown Earth motion. It claims the combination of two known motions creates a time-varying observational bias:

1. Earth’s axial precession (~25,771 years, ~5,028.8″/century) - this IS corrected in standard reductions
2. Earth’s apsidal precession (~111,772 years) - the slow rotation of Earth’s orbital ellipse

The model proposes that standard corrections apply these as constant rates, but the actual effect on Mercury observations varies over the 335,317-year Earth Fundamental Cycle (the LCM of both periods). This variation manifests as the ~5.0″/century change in observed precession that the model predicts.

The key claim: It’s not that a motion is missing from corrections, but that the interference pattern between two known long-period cycles produces time-dependent residuals. Whether this is physically valid or the corrections already account for this requires detailed analysis of the IERS coordinate transformation pipeline.

See The Measurement Chain above for details on how Earth’s position is used in the measurement process.

Q3: If General Relativity has been independently confirmed by other tests, isn’t it likely correct for Mercury too?

The challenge:

• GPS requires GR corrections to function (time dilation)
• LIGO has directly detected gravitational waves
• Light bending was measured during the 1919 eclipse (Eddington)
• Shapiro delay (radar signal time delay) has been confirmed
• If GR is correct in all these contexts, the simplest hypothesis is it’s also correct for Mercury

Model response: The model does NOT claim GR is wrong as a theory. It proposes that this specific test (Mercury perihelion) may have an alternative explanation based on reference frame effects, while other GR effects remain valid.

Why this is not self-contradictory:

• GPS time dilation: Measures local clock rates, not angular positions relative to distant objects. Reference frame rotation doesn’t affect local time measurement.
• Gravitational waves (LIGO): Detects local spacetime strain using laser interferometry. No dependency on ICRF or Earth’s orbital motion.
• Light bending (Eddington): Measured angular deflection during a single event (eclipse). Short timescale (~hours) means Earth’s precession motions are negligible.
• Shapiro delay: Measures radar signal travel time. Again, local measurement not dependent on long-period reference frame effects.

The key distinction: Mercury’s perihelion precession is uniquely sensitive to reference frame effects because it requires comparing angular positions over decades to centuries. The model proposes that the long-period components of Earth’s motion (~25,794 and ~111,772 year cycles) may not be fully corrected in these measurements. Other GR tests either work on shorter timescales or measure local physical quantities independent of angular reference frames.

Caveat: This distinction does not prove the model is correct. It explains why an alternative interpretation for Mercury’s precession could be consistent with other confirmed GR effects.

Q4: Measurement precision is now ~±0.0015″/century. How can the anomaly change at all?

The challenge:

• MESSENGER (Park et al. 2017) measured 575.31 ± 0.0015″/century in ICRF — unprecedented precision
• No drift has been observed in recent decades
• Any change in the rate at this scale should already have been detected

Model response: The model predicts the observed precession RATE (measured in ″/century) will change over time. Let’s clarify what this means:

Model prediction (heliocentric rate of change):
Year 2000: ~569.46″/century (geocentric: ~5,598.26″)
Year 2100: ~564.05″/century (geocentric: ~5,592.85″)
Change in rate over 100 years: ~5.0″/century

MESSENGER mission (2011-2015):
Measured: 575.31 ± 0.0015″/century at epoch ~2013

BepiColombo (orbit November 2026, science operations April 2027):
Gap from MESSENGER: ~14 years
Model's expected rate change over 14 years: ~0.70″/century
Model predicts: ~574.61″/century or lower
  (vs MESSENGER's 575.31″/century)
Predicted gap is ~500× larger than MESSENGER precision (±0.0015″/century)

Key distinction: MESSENGER measured the precession rate at one epoch with high precision. The question is whether this rate will be the same when measured again years later.

Methodological caveat: The reported “575.31″/cy” comes from fitting a GR-inclusive ephemeris to spacecraft data (Park et al. 2017  fit the PPN parameter β jointly, finding β ≈ 1). Conceptually, this is equivalent to measuring the Newtonian baseline (~532″/cy) and adding the assumed GR contribution (~43″/cy). If BepiColombo’s analysis pipeline applies the same GR-inclusive fit, any change in the underlying perihelion advance from frame effects may be absorbed into a slightly different best-fit β, into residuals, or into the orbital baseline — rather than showing up cleanly as a drift in the reported total. A definitive test requires reporting the raw measured perihelion advance independent of an assumed GR baseline, or explicitly tracking the Newtonian and GR contributions separately between epochs.

What the model predicts:

• Over 14 years (MESSENGER → BepiColombo), the rate should decrease by ~0.70″/century
• This is ~500× larger than MESSENGER’s measurement uncertainty — easily detectable if real and if the analysis pipeline reports it transparently (see methodological caveat above)

What GR predicts: The rate should be constant at ~575.31 + 5,028.8 = ~5,604″/century geocentric (within measurement uncertainty)

Current status: No drift can be detected with only one high-precision measurement epoch (MESSENGER). BepiColombo will provide the second epoch needed for this test — provided the analysis pipeline reports the raw measured perihelion advance, not a GR-inclusive fit total. See Mercury Precession: The BepiColombo Test for a detailed breakdown of the two possible outcomes and what each would mean for the model.

Q5: Why does the model’s “artifact” match the GR prediction at all?

The challenge:

• GR predicts ~42.98″/century from fundamental constants (G, M, c, a, e)
• The observed anomaly closely matches this value
• If the anomaly is a reference frame artifact, why does it equal the value GR predicts from first principles?
• This match seems to leave little room for an alternative explanation

Model response: This is the strongest single argument against the model’s interpretation. The match is real and striking. Three responses:

1. The historical 43″ value was established at the exact epoch where the model also predicts ~43″: The “anomaly” entered the literature when Le Verrier (1859) flagged the discrepancy and Newcomb (1882) refined it to ~43″/century — late-19th-century methodology, effective epoch ~1900. The model’s predicted reference-frame anomaly at ~1900 is ~43.01″ — essentially identical to Newcomb’s number. This is not a tuning artifact: the model is deterministic, and the time-evolution of Earth’s perihelion precession produces a value that happens to be ~43″ precisely when Newcomb measured it. By J2000 the model has decreased to ~38.02″, and the model predicts continued decrease. Under GR, 43″ is a permanent fundamental constant; under the model, the 1900 match is exactly what time-varying frame effects predict for that specific epoch, and the rate should drift. The decisive test is whether the value still equals 43″ a century later (GR) or has drifted (model) — addressed in Q4 and the BepiColombo Test.

2. GR-inclusive measurement methodology: As discussed in §3.4 The Measurement Chain, modern measurements fit a GR-inclusive ephemeris. The reported ~43″ residual is what comes out of fitting Mercury’s orbit jointly with the PPN parameter β (which converges to ≈1). The β fit itself is not circular — β could have come out different — but the methodology constrains the residual to behave as GR predicts.

3. Alternative derivation by Gerber (1898): Paul Gerber derived the same ~43″ formula 17 years before Einstein, from the assumption that gravity propagates at the speed of light. Mainstream physics considers Gerber’s derivation flawed, but it shows that the same numerical result can emerge from non-GR assumptions.

Honest assessment: At a single epoch, the match is strong evidence for GR. The decisive test is whether the rate stays constant (GR) or drifts (model) — see Q4 above and the BepiColombo Test.

Q6: If GR’s universal formula matches perihelion precession for all 8 planets, how can the model claim Mercury’s anomaly is a frame effect but not the others?

The challenge:

• GR predicts perihelion precession for all planets using the same formula: Δϖ_GR = 6πGM/(ac²(1−e²)) per orbit
• Venus: ~8.6″/century relativistic contribution
• Earth: ~3.8″/century
• Mars: ~1.4″/century
• The same universal formula matches across the inner solar system — strong evidence for GR’s general validity

(These are theoretical predictions calculated from the GR formula using each planet’s orbital elements. See Clemence 1947  for historical derivations.)

Model response: The model is not claiming Mercury is uniquely a frame effect while the others are GR — rather, the model proposes that all planet “anomalies” arise from how Earth’s reference-frame motion projects onto each planet’s geometry. The model predicts all 7 non-Earth planets with R² ≈ 1.0000 across the full Earth Fundamental Cycle (see Other Planets above). The mechanism varies per planet:

1. Well-defined vs poorly-defined perihelia produce different signatures: As discussed in Venus: A Contrasting Case, high-eccentricity planets (Mercury, e ≈ 0.20564) show geometric modulation — Earth’s frame motion projects onto the well-defined perihelion direction. Low-eccentricity planets (Venus, e ≈ 0.00678) show ERD² × periodic effects — the poorly-defined perihelion primarily reflects Earth’s own rate variations. The Fluctuation Range column in the Other Planets table makes this visible: Mercury’s range is ~-180 to +202″/cy; Venus’s is ~-1,353 to +1,231″/cy (~7× larger), exactly because of the eccentricity contrast.

2. The same GR-inclusive measurement caveat applies to all planets: As discussed in Q4 and §3.4 The Measurement Chain, modern measurements of all planets fit a GR-inclusive ephemeris (Park et al. 2017 , DE440/441 ) with PPN β as a free parameter. The “match with GR’s universal formula” emerges from a fit that constrains β ≈ 1 across the dataset — strong evidence at one epoch, but doesn’t directly test whether the underlying perihelion advances drift over time.

3. The decisive test is the same: rate stability vs drift: Under GR, all planet perihelion advances should be permanent constants. Under the model, all should drift on Earth Fundamental Cycle timescales — the magnitude of drift varies per planet (see the Fluctuation Range column in Other Planets). BepiColombo provides Mercury’s near-term test; longer baselines for outer planets would be needed.

Scientific Position Summary

The model’s position: This is not a claim that General Relativity is wrong — GR has been confirmed independently across many other tests. For Mercury’s perihelion specifically, the model offers an alternative interpretation that makes a different time-evolution prediction, allowing BepiColombo and longer-baseline future measurements to distinguish the two.

5. Eccentricity Cycles and Milankovitch Theory

This section provides a fair presentation of Milankovitch theory and modern orbital solutions, then examines the model’s alternative proposal.

Milankovitch Theory: A Fair Presentation

Milutin Milankovitch (1879-1958) was a Serbian mathematician and astronomer who developed the astronomical theory of climate change. His work, culminating in Canon of Insolation and the Ice-Age Problem (1941), proposed that Earth’s ice ages are driven by variations in solar radiation received at high northern latitudes during summer.

The Milankovitch cycles and their constituents:

† Axial and apsidal precession are not Milankovitch climate drivers in their own right — they are the two physical motions that combine to produce the climatic precession (the third Milankovitch cycle). Listed here for completeness.

*The “climatic precession” is what determines insolation timing — where the equinoxes fall relative to perihelion. Berger (1978)  identified dominant periods near ~23.7, ~22.4, and ~19.0 kyr, jointly summarized as ~23,000 years in popular accounts. The math: 1/T_climatic = 1/T_axial + 1/T_apsidal ≈ 1/25,800 + 1/112,000 ≈ 1/21,000 yr (mean); the multiple spectral peaks arise because the apsidal precession itself has internal structure from different planetary perturbation modes. Milankovitch’s original 1941 work used different numbers based on then-current ephemerides.

Key insight: Milankovitch identified that summer insolation at 65°N is the critical parameter for ice sheet growth/decay. When northern summers are cool (low insolation), snow survives year-round and ice sheets can grow.

Historical validation: The theory was largely ignored until Hays, Imbrie & Shackleton (1976)  demonstrated that deep-sea sediment records show spectral peaks at the predicted Milankovitch frequencies. This landmark paper, “Variations in the Earth’s Orbit: Pacemaker of the Ice Ages,” established Milankovitch theory as the foundation of paleoclimatology.

The Eccentricity Spectrum: What Milankovitch Actually Calculated

Modern long-term integrations (Laskar et al. 2004 , the La2004 solution — full N-body integration of all 8 planets, Moon, solar oblateness, and GR corrections, valid for ~50 Myr beyond which chaos limits predictability) decompose Earth’s eccentricity variation into spectral components driven by interactions between the inner planets’ orbital precession frequencies (g₂ Venus, g₃ Earth, g₄ Mars, g₅ Jupiter):

Eccentricity variations are quasi-periodic, not strictly periodic — the dominant terms involve interactions between planetary orbital frequencies, and the ~100k-year “cycle” cited in paleoclimate literature is actually the combined effect of the ~95k and ~125k components, producing a quasi-periodic signal with average period near 100k years. There is no single ~100k spectral peak in eccentricity itself, but the combination produces something that looks like one in time-domain data — which is what climate records typically capture.

Modern Orbital Solutions (Laskar et al.)

Beyond the spectral decomposition above, the Laskar group has produced refined long-term orbital solutions including La2010  (Laskar et al. 2011 — updated planetary masses; provides eccentricity, obliquity, and precession for Earth).

Laskar’s eccentricity predictions (from La2004):

Physical basis: These predictions derive from Lagrange-Laplace secular perturbation theory, which calculates how planetary gravitational interactions cause slow orbital changes. The mathematics involves:

• Fourier decomposition of orbital elements
• Secular (long-term averaged) perturbation equations
• Numerical integration over millions of years

The “100,000-Year Problem”

Despite Milankovitch theory’s success, a major puzzle remains:

The paradox:

• Eccentricity causes only ~0.2% variation in total annual solar energy
• Obliquity causes ~10% variation in polar summer insolation
• Yet for the past ~1 million years, ice ages follow a ~100k pattern, not the stronger ~41k obliquity signal

This is genuinely puzzling: If orbital forcing drives ice ages, why does the weakest forcing (eccentricity) produce the strongest climate signal?

Mainstream proposed solutions:

1. Ice sheet nonlinear dynamics (Imbrie et al. 1993 ):

   • Ice sheets have internal dynamics with ~100k timescales
   • Small eccentricity forcing triggers large ice sheet responses
   • Threshold effects and hysteresis create apparent ~100k cycles

2. Eccentricity modulates precession (Raymo 1997 ):

   • Precession’s climate effect depends on eccentricity
   • High eccentricity amplifies precession’s seasonal contrast
   • The ~100k signal is precession amplitude modulation, not direct eccentricity forcing

3. Carbon cycle feedbacks (Paillard 1998 ):

   • Ocean-atmosphere CO₂ exchange has long time constants
   • Eccentricity cycles modulate carbon storage in oceans
   • The ~100k climate response is amplified by carbon feedbacks

4. Antarctic ice sheet control (Raymo et al. 2006 ):

   • Southern Hemisphere ice sheets may be more sensitive to eccentricity
   • The ~100k signal originates from Antarctic, not Greenland

The 100,000-year problem remains “one of the most significant unresolved questions in climate science” (Imbrie et al. 1993). No single explanation has achieved consensus. Recent work continues to debate the question:

• 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 peak identification 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
• The Mid-Pleistocene Transition — the shift from 41-kyr to ~100-kyr glacial cycles around 1 million years ago — remains one of paleoclimatology’s great unsolved puzzles

The Model’s Alternative: A ~20,957-Year Eccentricity Cycle

The Holistic Universe Model proposes a fundamentally different view of Earth’s eccentricity:

The claim: Earth’s eccentricity varies in a ~20,957-year cycle with range ~0.0140 to ~0.0167, not the ~100k/~400k cycles with range 0.0005-0.058 predicted by Laskar.

The mechanism (from Eccentricity):

Two counter-rotating motions:
1. Earth around its wobble center: ~25,794 years (clockwise)
2. Earth's perihelion point around Sun: ~111,772 years (counter-clockwise)

Meeting frequency = 1/~25,794 + 1/~111,772 = 1/~20,957

They meet every ~20,957 years → eccentricity cycle

The solstice-eccentricity correlation (claimed by the model):

Model’s key dates:

• Maximum eccentricity: ~1246.03125 AD (perihelion at December solstice)
• Next minimum: ~11,725 AD (perihelion at June solstice)

Questions for This Eccentricity Claim

This claim requires careful examination because it contradicts well-established celestial mechanics:

Q1: How can the model claim eccentricity varies on a geometric meeting frequency when standard physics ties eccentricity to orbital energy?

The challenge:

• In standard physics, eccentricity is determined by orbital energy and angular momentum — conserved quantities that change only through gravitational perturbations (the basis of Laskar’s secular theory)
• Standard axial precession (gyroscopic wobble of Earth’s spin axis) is an angular reorientation; it does not change the orbit’s energy or angular momentum
• The model’s claim that two precession cycles’ meeting frequency produces an eccentricity oscillation has no direct analog in standard celestial mechanics

Model response: The model proposes a fundamentally different mechanism for eccentricity variation than Laskar’s secular theory: instead of gravitational perturbations modifying orbital energy and angular momentum, the model proposes a geometric construction in which Earth orbits a moving wobble center, producing an apparent eccentricity oscillation as that wobble center precesses around the Sun. This is a substantive alternative that, if correct, replaces Laskar’s secular theory for Earth’s eccentricity rather than supplementing it. The model owes a first-principles physical justification for the wobble-center construction itself — this is an open theoretical question.

Q2: How does the model reconcile its narrow ~0.003 range with Laskar’s computed ~0.057 range over millions of years?

The challenge:

• Modern orbital integrations (Laskar et al. 2004, La2004) compute that Earth’s eccentricity has varied between ~0.0005 and ~0.058 over the past several million years — a range of ~0.057
• The model’s ~20,957-year cycle would predict a range of only ~0.003 (~0.0167 − ~0.0140), about 20× narrower
• This range contradiction holds independently of the climate-proxy debate: Laskar’s range is computed from secular perturbation theory using well-verified planetary masses and Newton’s laws, not from proxies

Model response: The model proposes that Laskar’s 0.0005–0.058 range is incorrect. Laskar’s range is computed from secular perturbation theory; it has not been directly observed because the relevant timescales (millions of years) are far longer than direct precision-measurement records (~centuries). Geological proxies typically cited as evidence (ice cores, marine sediments) involve orbital tuning and have circularity concerns — see §6 (The 100,000-Year Glacial Cycle). The model rejects the 95k/125k/405k cycle structure for eccentricity entirely; the geological ~100k climate signal usually attributed to Laskar’s eccentricity forcing may actually reflect inclination precession (~112k), not eccentricity. The two predictions are mutually exclusive: future direct observations of Earth’s eccentricity over the next centuries to millennia will discriminate.

Q3: Why would the model’s eccentricity correlate with solstice alignment when no mainstream mechanism connects spin-axis orientation to orbital shape?

The challenge:

• The model claims eccentricity peaks when perihelion aligns with the December solstice (and minimizes at June solstice)
• In standard secular theory, eccentricity and longitude of perihelion evolve via coupled Lagrange-Laplace equations driven by gravitational perturbations — neither is determined by Earth’s axial-tilt orientation
• Solstice dates are determined by Earth’s spin axis and orbit geometry independently of eccentricity magnitude
• No mainstream mechanism connects spin-axis–perihelion alignment to orbital shape

Model response: Under the model, the solstice-eccentricity correlation is a direct geometric consequence of the wobble-center construction: when perihelion aligns with the December solstice, the wobble center’s offset adds constructively to the orbital eccentricity (maximum e); when it aligns with the June solstice, the offset partially cancels (minimum e). This is internally consistent with the model’s geometric framework but has no analog in standard celestial mechanics — and the wobble-center construction itself still requires first-principles physical justification.

Q4: How does the model’s single H/16 climatic precession period account for Berger’s multi-peak spectrum (~19-24 kyr)?

The challenge:

• Berger’s classical Fourier expansion of e × sin ϖ̄ produces a multi-peak climatic-precession spectrum, with terms at 23,716 / 22,428 / 23,159 / 19,155 / 18,976 / 16,469 yr (each corresponding to a planet-specific eigenmode g_j + k for Jupiter, Venus, Mercury, Earth, Mars, and Saturn respectively)
• The model presents climatic precession as a single H/16 = ~20,957-year cycle
• A single-peak description would seem incompatible with the observed multi-peak structure of secular theory

Model response: At the level of the Solar System Resonance Cycle (~2,682,536 yr — the period over which every planetary cycle returns to alignment), every Berger climatic-precession peak matches an integer-fraction divisor of this cycle within <0.4% — including Saturn (16,469 yr → integer n = 163, giving 16,457 yr at 0.07%). Each integer decomposes as n = 104 + δ, where 104 corresponds to axial precession (k = H/13) and δ is the planet’s eigenfrequency contribution. The largest planetary offsets from the centroid (n = 128, equal to H/16) — Mars +13 and Venus −8 — are direct Fibonacci numbers. The model’s H/16 is the centroid of the Berger spread; the spread itself is recoverable from integer fractions of the resonance cycle. The same framework also exposes a triple identity at H/8 = ~41,915 yr: Jupiter’s ICRF perihelion period, Saturn’s ecliptic perihelion period, and the obliquity cycle all share this divisor. See Supporting Evidence §10: Climatic precession comparison for the full table.

What would further validate this claim:

• A first-principles derivation of the wobble-center geometry from observed Earth-Sun dynamics
• Or independent observational confirmation that apparent eccentricity measured from Earth peaks at the predicted solstice alignments

Detailed Comparison: Model vs. Laskar

Where they agree:

• Current eccentricity value (~0.0167)
• Current trend (decreasing)

Where they fundamentally differ:

• Cycle period (~~20,957 yr vs ~95k–400k yr)
• Amplitude range (~0.003 vs ~0.057)
• Physical mechanism

A subtle observation: La2004 predicts e(J2000) = 0.01670, while the modern DE440 ephemeris (observation-fitted) gives e(J2000) = 0.01671022 — a difference of ~10⁻⁵, with La2004 slightly lower than observed. This means La2004’s secular decrease has already overshot the actual trajectory by J2000. The magnitude is small (within typical secular-theory precision), but if the trend continues, the actual eccentricity decrease may proceed slower than La2004 predicts — closer to the model’s slower decrease. The model is calibrated to DE440 by construction, so it agrees with the observed J2000 value. See the eccentricity comparison image for the curves visualized.

Evidence Assessment

Evidence supporting Laskar/Milankovitch:

1. Fundamental physics: Laskar’s calculations derive from Newton’s laws and observed planetary masses. The mathematics is well-established secular perturbation theory.

2. Geological proxy agreement: Ice cores, marine sediments, and speleothems all show ~100k cycles over the past million years (though with circularity concerns - see §6).

3. Spectral analysis: Climate records show power at Milankovitch frequencies (~100k, ~41k, ~23k climatic precession), with the ~23k attributed to climatic precession (not eccentricity). There is no clearly identified ~~20,957-year peak attributed specifically to eccentricity in the records.

Evidence supporting the model’s alternative:

For the full discussion of model-supporting evidence (100k problem, Muller & MacDonald 1997, MPT mechanisms, recent research), see Supporting Evidence §1: The 100,000-Year Problem. Headline points:

• The 100,000-year problem remains unsolved (Barker 2025, Mitsui 2025, Lisiecki 2023)
• Geological proxies have circularity concerns from orbital tuning (LR04, ice core chronologies)
• Laskar’s predictions for deep time (>~10,000 years) are unverified extrapolations
• Spectral evidence (Muller & MacDonald 1997 ) shows the ~100k climate signal is incompatible with eccentricity’s split-peak structure (~95k + ~125k)
• La2004’s J2000 value (0.01670) is slightly lower than DE440’s observed value (0.01671022) — a small overshoot consistent with the model’s slower-decrease prediction (see “subtle observation” above)
• The ~400k-year absence: eccentricity’s theoretically strongest component (~405k from g₂−g₅) is largely absent from climate records of the past 1.2 Myr (Muller’s work on ice ages ) — difficult to explain if eccentricity drives climate

Testable Predictions

The numerical comparison of model vs La2004 at specific epochs (5,000 / ~11,725 / 27,000 AD) is in Eccentricity: Numerical Comparison. Direct verification of either prediction requires geological timescales.

Indirect tests (potentially shorter timeframe):

• If the model is correct, the ~100k glacial cycles should show evidence of ~112k periodicity when dated without orbital tuning
• If Laskar is correct, improved dating methods should continue to confirm ~100k timing

Connection to the 100,000-Year Glacial Cycle (§6)

The model’s eccentricity claim is linked to its ice core argument:

The model proposes:

1. The ~100k glacial signal is actually ~112k (inclination precession)
2. Ice core dating has ~10% systematic error from orbital tuning
3. The model’s ~20,957-year eccentricity cycle is separate from the glacial ~100k pattern

Implication: If the model is correct about ice core chronology, it would support reinterpreting the ~100k climate signal as inclination rather than eccentricity. This doesn’t directly validate the ~20,957-year eccentricity cycle, but it removes one objection (that ~100k proxies contradict it).

An open tension: The dominant climate signal in proxy records sits at ~100k years, while the model identifies the inclination precession period as ~112k — a ~10% gap. The model addresses this by arguing that ice core dating has ~10% systematic error from orbital tuning, so the actual period could be ~112k once corrected. Independent dating methods (U-Th speleothems, O₂/N₂ ratio) could test whether the true climate period is closer to ~100k or ~112k. Until then, this discrepancy remains unresolved.

See §6 (The 100,000-Year Glacial Cycle) for detailed analysis of dating methods and the circularity problem.

Scientific Position Summary

The model’s position: This is not a claim that planetary gravitational interactions are absent — N-body perturbations are well-established physics. For Earth’s eccentricity specifically, the model offers an alternative interpretation: that the dominant variation reflects a single ~20,957-year cycle tied to the meeting frequency of Earth’s two precession motions, rather than the quasi-periodic Lagrange-Laplace solution. The two predictions diverge enough over centuries-to-millennia baselines that future high-precision ephemerides could distinguish them.

6. The 100,000-Year Glacial Cycle

This section examines the dominant ~100,000-year cycle in glacial-interglacial climate records and asks whether it reflects eccentricity (Milankovitch) or inclination precession (the model’s alternative). Ice core chronology is the primary evidence base, so the section first presents the standard view of how ice cores are dated, then examines the model’s reinterpretation and the evidence for and against it.

How Ice Cores Are Dated

Modern ice core chronology uses multiple independent methods:

1. Annual layer counting (most precise for recent ice):

• Visual stratigraphy (summer/winter layers)
• Chemical signatures (seasonal dust, sea salt, isotopes)
• Electrical conductivity measurements (ECM)
• Precision: ±1% for Holocene (~11,700 years), ±2-3% for glacial periods
• Limitation: Layers become too thin to resolve beyond ~60,000-100,000 years

2. Volcanic markers (absolute tie points):

• Sulfate spikes from known eruptions
• Tephra (volcanic ash) layers with unique chemical signatures
• Examples: Toba (74 ka), Laacher See (12.9 ka), Campanian Ignimbrite (39 ka)
• Strength: Provides independent absolute dates where identified

3. Gas synchronization (global correlation):

• Methane (CH₄) variations are globally synchronous (within ~50 years)
• Links Greenland and Antarctic records precisely
• Allows transfer of well-dated Greenland chronology to Antarctic cores
• Precision: ±50-200 years for the synchronization itself
• Limitation: Still requires one record to have an independent chronology

4. O₂/N₂ ratio dating (Kawamura et al. 2007 ):

• Trapped air’s O₂/N₂ ratio correlates with local summer insolation
• Provides an independent orbital constraint without assuming eccentricity cycles
• Used to construct and validate the DFO-2006 chronology for the Dome Fuji core (and cross-checked against AICC2012 in the overlap interval)
• Key feature: This method constrains timing to precession cycles (~23 ka), not ~100 ka cycles

5. Radiometric dating (uranium-series, ¹⁴C):

• Radiocarbon (¹⁴C) useful to ~50 ka
• U-series dating of speleothems synchronized to ice cores
• Strength: Completely independent of orbital assumptions

6. Orbital tuning (matches isotopes to insolation):

• Adjusts chronology so δ¹⁸O variations match calculated insolation changes
• Assumes Milankovitch theory is correct
• Critical issue: This method is circular if used to test Milankovitch theory!

Major Ice Core Chronologies

AICC2012 (Antarctic Ice Core Chronology 2012) is the current community standard for long Antarctic records. It uses:

• Layer counting where possible
• Gas synchronization between cores
• Glaciological ice flow modeling
• Minimal orbital tuning (only for transitions where other constraints are weak)

Chronology Uncertainties

The Circularity Problem

Standard orbital tuning works by adjusting the ice core timescale so that δ¹⁸O (ice isotope) variations match calculated summer insolation at high latitudes. This assumes:

1. Climate responds predictably to insolation
2. Milankovitch’s insolation calculations are correct
3. The response time (phase lag) is known

This is circular when testing Milankovitch theory. If you tune your chronology to match Milankovitch predictions, you cannot then use that chronology to validate Milankovitch theory.

The model’s argument: If chronologies are orbitally tuned assuming ~100k eccentricity cycles, they would artificially compress a true ~112k signal to appear as ~100k.

Counter-argument: Modern chronologies (especially AICC2012) minimize orbital tuning:

• Use multiple independent constraints
• Only apply tuning where other methods fail
• The O₂/N₂ method constrains precession timing independently
• Volcanic markers provide absolute tie points unaffected by tuning

The Model’s Claim

The Holistic Universe Model proposes that the dominant ~100k-year pattern in glacial-interglacial cycles reflects the inclination precession period (~111,772 years), not Milankovitch’s eccentricity cycles (~95k/~125k). This would require that ice core chronologies have a systematic error of ~10% for older ice.

This is a significant claim that deserves careful examination.

The Inclination Hypothesis (Muller & MacDonald 1997)

The model’s claim that the ~100k glacial signal reflects inclination precession is not unique. A peer-reviewed hypothesis in the Proceedings of the National Academy of Sciences made a similar argument:

Muller & MacDonald (1997) : “Spectrum of 100-kyr glacial cycle: Orbital inclination, not eccentricity”

Key findings:

1. Spectral mismatch with eccentricity: The observed ~100k peak in climate data is narrow and well-defined. However, the eccentricity spectrum shows peaks at ~95k and ~125k years - not a single ~100k peak. The climate data doesn’t match eccentricity’s spectral shape.

2. Inclination provides better match: Earth’s orbital inclination relative to the invariable plane (perpendicular to the solar system’s angular momentum vector) has a dominant ~100k period that closely matches the climate spectrum.

3. Bispectral analysis: The authors found that the bispectrum (which reveals phase relationships between frequency components) of the climate data matches the inclination signal, but not eccentricity.

4. Proposed mechanism: The authors suggested extraterrestrial accretion (meteoroids or interplanetary dust) as the link between inclination and climate - though this mechanism remains speculative.

Relevance to the model:

Important distinction: Muller & MacDonald found that inclination changes match the ~100k climate signal as currently dated. The Holistic Universe Model goes further, proposing that the true period is ~112k and the ~100k appearance is a dating artifact. These are related but distinct claims.

Status of the Muller-MacDonald hypothesis: This hypothesis remains controversial. The proposed mechanism (extraterrestrial accretion) was subsequently rejected by the community. However, the spectral evidence — the mismatch between eccentricity’s split spectral peak (95k + 125k) and the climate record’s single narrow peak near ~100k — has not been definitively refuted and remains a recurring concern in spectral analyses of glacial cycles. The Holistic Universe Model provides an alternative mechanism (inclination precession at ~111,772 years from two counter-rotating reference points) that does not rely on dust. The “100,000-year problem” remains unresolved, and the Muller-MacDonald work demonstrates that serious scientists have questioned the eccentricity explanation.

The Mechanism: Inclination → Obliquity → Climate

Muller & MacDonald (1997) proposed extraterrestrial dust accretion as the link between inclination and climate; that mechanism was rejected. The Holistic Universe Model proposes a different link that uses only standard Milankovitch physics.

The model’s obliquity formula (see Supporting Evidence §10–§11) contains two equal-amplitude cosine terms:

1. Axial tilt effect at H/8 = ~41,915 yr — matches mainstream’s ~41k obliquity cycle (Berger’s s₃ + k beat)
2. Inclination tilt effect at H/3 = ~111,772 yr — Earth’s orbital plane oscillation vs. the invariable plane

Both contribute ±0.63603° to Earth’s tilt. The prediction is therefore that Earth’s obliquity oscillates at two periods, not one: ~41k AND ~112k.

The implied climate chain is:

Inclination precession (~112k yr) → modulates obliquity at ~112k yr → modulates seasonal insolation contrast → drives glacial cycles at ~112k yr

Every step after “modulates obliquity” is mainstream Milankovitch physics. The mechanism needs no dust, no exotic forcing — only that mainstream secular theory has distributed the H/3 inclination tilt component across small spectral terms rather than recognizing it as one ~112k peak.

Visual evidence — Vostok ice core (Petit et al. 1999):

Both cycles are overlaid on Vostok temperature and CO₂ data:

• The blue ~112k curve tracks the major glacial-interglacial envelope; ~3.5 cycles span the past 400 ka, aligning with the four interglacial peaks (Holocene, Eemian, MIS 7, MIS 9, MIS 11)
• The orange ~41k curve tracks the secondary oscillations within each glacial period
• The Last Glacial Maximum coincides with the model’s most recent maximum-inclination year (23,204 BC) — within the conventional LGM phase (~24,500 to ~17,000 BC), with the standard ~4 kyr lag to peak ice (~19,000 BC)
• Three inclination cycles span one Earth Fundamental Cycle (H = 335,317 yr), corresponding to the three major glacial-interglacial transitions in the visible record

The temperature signal is visually decomposable into the sum of these two oscillations — direct empirical support for the model’s two-component obliquity claim.

Testability: A spectral analysis of Earth’s obliquity over 0-800 ka (Berger 1978, La2004) for power at ~112k would confirm or refute the prediction analytically. If obliquity has measurable ~112k power, the inclination → climate mechanism is established without invoking new physics.

Questions for This Climate Reinterpretation

For the model’s claim to hold (that ~100k is actually ~112k), several conditions must be met. The following questions probe the strongest objections.

Q1: How can a ~10% systematic dating bias affect every independent method uniformly?

The challenge: Layer counting, volcanic markers, gas synchronization, U-series speleothem correlations, and radiometric methods are all independent of orbital tuning. Yet all major Antarctic (EPICA, Vostok, Dome Fuji), Greenland (NGRIP, GRIP, GISP2), and marine chronologies converge on ~100k for the dominant glacial pattern. For the model’s claim to hold, all of these methods would need to carry the same ~10% bias — and no mechanism has been proposed for why they would all err in the same direction. Furthermore, the required ~10% offset sits at the very upper limit of established uncertainty estimates: the Chronology Uncertainties table above shows ±5-10% only for ice older than 400 ka, with much tighter bounds (±1-2% to ±4-6%) for younger ice that still shows the ~100k pattern.

Model response: AICC2012 still applies orbital tuning for ice older than ~400 ka, where other methods are weakest. Marine δ¹⁸O stacks (e.g., Lisiecki & Raymo 2005 ) explicitly use orbital tuning, so they cannot independently confirm ice core chronology. Recent cycles (under ~400 ka) have stronger non-tuning constraints and would be the place to look for systematic offsets — a reanalysis targeting this question is not yet published.

Q2: What about non-orbitally-tuned dating like the O₂/N₂ ratio?

The challenge: Kawamura et al. (2007)  demonstrated that O₂/N₂ ratio dating uses precession cycles (~23 ka), not ~100k cycles, to constrain timing. It validates the AICC2012 chronology within ~±2%. If this method is independent of any ~100k assumption and still confirms conventional timing, the ~10% systematic-bias hypothesis is hard to defend.

Model response: The O₂/N₂ method constrains precession-band timing (~23 ka), not the ~100k peak directly. A reanalysis specifically targeting whether the ~100k climate peak shifts toward ~112k when O₂/N₂-anchored chronologies are used would test this directly. To date, this question has not been the explicit focus of any published spectral analysis.

Q3: Is the spectral peak really broad enough for ~~112k to fit?

The challenge: The dominant glacial signal in spectral analyses appears as a peak near ~100k. If this peak is sharp enough to exclude ~112k, the model’s reinterpretation fails; if it’s broad enough to span 80-120 ka, both ~100k and ~112k are in scope.

Model response: Published spectral analyses (beginning with Hays, Imbrie & Shackleton 1976  and subsequent refinements) typically report the peak spanning 80-120 ka — wide enough to accommodate ~112k without contradiction. More importantly, Muller & MacDonald (1997)  showed that the climate peak’s shape (a single narrow peak rather than the 95k + 125k split predicted by eccentricity) fits inclination better than eccentricity regardless of which exact period is chosen. The “100,000-year problem” — the persistent failure of eccentricity to mechanistically explain the dominance of ~100k cycles over the expected ~41k obliquity signal — adds independent reason to reconsider the eccentricity attribution.

Evidence Assessment

This subsection catalogs evidence not already covered in Q1-Q3 above and “The Inclination Hypothesis” subsection.

Additional evidence supporting conventional ~100k chronology:

1. Speleothem correlations (Cheng et al. 2016 ): Cave deposits dated by U-Th (completely independent of ice) show the same ~100k pattern with consistent timing.

2. Marine sediment record details (Lisiecki & Raymo 2005 ): The LR04 stack of 57 benthic δ¹⁸O records shows ~100k cycles using biostratigraphy and magnetostratigraphy in addition to orbital tuning.

Evidence supporting the model’s alternative:

The core case for the model’s reinterpretation is laid out in Q1-Q3 above and in “The Inclination Hypothesis (Muller & MacDonald 1997)” subsection. For broader context — including MPT mechanisms, recent literature (Barker 2025, Mitsui 2025), and the ~400k-year absence problem — see Supporting Evidence §1: The 100,000-Year Problem.

Testable Predictions

If the model is correct:

1. Future improvements in ice core dating (better volcanic correlations, new radiometric methods) should reveal systematic offsets
2. The ~100k pattern should show slight stretching when dated without orbital tuning
3. Marine sediment records dated without tuning should show ~112k rather than ~100k

If the conventional chronology is correct:

1. New dating methods should continue to validate existing chronologies
2. The ~100k spectral peak should tighten with better data, not shift toward ~112k
3. Independent speleothem and marine records should maintain ~100k timing

Connection to Eccentricity Cycles (§5)

The model’s ice core claim is the mirror of its eccentricity claim. §5 argues that Earth’s eccentricity follows a single ~~20,957-year cycle rather than Laskar’s quasi-periodic beat; §6 argues that the ~100k climate signal in proxies reflects inclination precession (~~111,772 yr), not eccentricity. Together they form a self-consistent reinterpretation: if either fails empirically, the other is harder to defend.

See §5 (Eccentricity Cycles and Milankovitch Theory) above and Eccentricity for the full model explanation in the main documentation.

Scientific Position Summary

The model’s position: This is not a claim that Milankovitch theory is wholly wrong — orbital obliquity and precession remain well-supported climate drivers. For the dominant ~100k glacial cycle specifically, the model offers an alternative interpretation: the climate signal reflects inclination precession (~112k) acting through a second obliquity oscillation at H/3, alongside the well-known H/8 (~41k) obliquity cycle. The mechanism is therefore standard Milankovitch insolation forcing applied to an obliquity component that mainstream secular theory has distributed across small spectral terms rather than recognizing as one ~112k peak. Independent (non-orbitally-tuned) chronologies and direct spectral analysis of obliquity solutions could test this prediction.

7. Mathematical Framework

Coordinate System Definitions

ICRS (International Celestial Reference System):

• Origin: Solar system barycenter
• X-axis: Toward vernal equinox (J2000)
• Z-axis: Toward celestial north pole (J2000)
• Fixed relative to distant quasars

Model Reference Points:

EARTH-WOBBLE-CENTER:

• Definition: A reference point that Earth orbits around, used to simulate the axial precession “wobble”

• Location: At the center of Earth’s Axial Precession Orbit (APO), approximately 202,881 km from Earth’s center (0.001356 AU)

• Motion: Earth orbits this point clockwise (as viewed from above the North Pole) in ~25,794 years

• Physical interpretation: Represents the combined gravitational torque from Moon and Sun that causes Earth’s axis to precess

• Mathematical representation: In the 3D simulation, Earth is placed at orbital radius 0.001356 AU from this center point

• Derivation of the 0.001356 AU value: This is a calibrated model parameter, not a directly measured distance. It was determined by:

  1. Setting the orbital period to ~25,794 years (mean axial precession)
  2. Requiring the simulation to produce correct precession rate (~50.29″/year)
  3. Iterating the radius until the obliquity variation amplitude matches observations (~0.63603°)

  The value has no independent physical meaning - it’s the radius that makes the mathematical model reproduce observed precession behavior. Think of it as similar to how Ptolemy’s epicycle radii were calibrated to match planetary positions without corresponding to physical objects.

PERIHELION-OF-EARTH:

• Definition: A reference point near the Sun that represents the center of Earth’s elliptical orbit and the direction of Earth’s perihelion
• Location: Orbits the Sun at radius ~2,301,681 km (0.015386 AU) - this is the base eccentricity (arithmetic midpoint of the eccentricity cycle)
• Motion: Orbits the Sun counter-clockwise (prograde) in ~111,772 years
• Physical interpretation: Represents Earth’s apsidal precession - the slow rotation of Earth’s orbital ellipse due to planetary perturbations (mainly Jupiter)
• Mathematical representation: In the 3D simulation, this point’s angular position determines Earth’s longitude of perihelion

Derivation of Key Periods

The Holistic Universe Model derives all precession cycles from a single master period: 335,317 years (the “Earth Fundamental Cycle”). This value is determined empirically by fitting to the observed axial precession period.

Step 1: Starting from observed axial precession

The IAU 2006 precession constant gives:

Observed precession rate = 50.2879"/year (with small corrections)

This yields an instantaneous axial precession period of approximately ~25,771 years. However, this rate varies over time due to changing obliquity. The model derives its base value from the mean rate:

Axial precession (mean) ≈ 25,793.62 years

Step 2: Deriving the Earth Fundamental Cycle

The model proposes that precession cycles follow Fibonacci ratios. Given axial precession corresponds to the Fibonacci number 13:

Earth Fundamental Cycle = 25,793.62 × 13 = 335,317 years

Step 3: Deriving other periods

All other cycles follow from dividing by Fibonacci numbers:

Note on Fibonacci sequence: 1, 1, 2, 3, 5, 8, 13, 21, 34… The divisors 3, 8, 13, 16 are from this sequence (16 = 13 + 3, a Fibonacci-like combination). See Formula Derivation: Fibonacci Hierarchy for how this pattern extends to planetary cycles.

Core Simulation Equations

The 3D simulation calculates orbital parameters using these relationships:

1. Obliquity (axial tilt)

ε(t) = ε₀ + A × [-cos(2π(t-t₀)/P₃) + cos(2π(t-t₀)/P₈)]

Where:
ε₀ = 23.41354° (mean obliquity)
A  = 0.63603° (amplitude)
P₃ = 335,317/3 = ~111,772 years (inclination precession)
P₈ = 335,317/8 = ~41,915 years (obliquity cycle)
t₀ = -302,635 (anchor year)

2. Orbital inclination (to invariable plane)

i(t) = i₀ - A × cos(2π(t-t₀)/P₃)

Where:
i₀ = 1.48113° (mean inclination)
A  = 0.63603° (same amplitude as obliquity)
P₃ = ~111,772 years

3. Eccentricity

The eccentricity formula is more complex, involving the meeting frequency of two counter-rotating motions:

Meeting period = 1 / (1/P₁₃ + 1/P₃) = P₁₆

Where:
P₁₃ = 25,793.62 years (axial precession, clockwise)
P₃  = ~111,772 years (inclination precession, counter-clockwise)
P₁₆ = ~20,957 years (resulting eccentricity cycle)

The full eccentricity equation (see Formulas) produces values oscillating between ~0.0140 and ~0.0167.

4. Longitude of perihelion

λ(t) = 270° + 360°(t-t₀)/P₁₆ + corrections

With sinusoidal corrections for the apsidal motion.

Comparison with Standard Formulas

Key observation: The model matches current observed values well, but diverges significantly from Laskar/Berger predictions for deep time. Both approaches are theoretical extrapolations that cannot be directly verified for ancient/future periods.

Error Analysis

Sources of uncertainty in model calculations:

1. Anchor year (t₀ = -302,635)

• Derived by fitting to J2000 observed values
• Uncertainty: The anchor year is not independently constrained
• Effect: Shifts phase of all cycles but not periods

2. Amplitude (A = 0.63603°)

• Derived from observed obliquity range
• Current observed range: 22.1° to 24.5° over ~41k years
• Model predicts: 22.21° to 24.72° (combined range from two ±0.63603° components)
• Uncertainty: ±0.01° based on current measurements

3. Mean values

• Mean obliquity (23.41354°): Based on IAU J2000 value, uncertainty <0.001°
• Base eccentricity (0.015386): Model-derived arithmetic midpoint of eccentricity cycle (time-averaged mean = 0.015386)

4. Period values

Propagated uncertainty in predictions:

For year 3000 AD (1,000 years from J2000):

• Obliquity: Model predicts 23.3103° ± 0.02°
• Eccentricity: Model predicts 0.01657 ± 0.0001

For year 12,000 AD (10,000 years from J2000):

• Obliquity: Model predicts 22.5358° ± 0.2°
• Eccentricity: Model predicts 0.01403 ± 0.001

Comparison with La2004 predictions: See the master table in §3.2 — Obliquity Predictions for the full year-by-year comparison.

Important caveat: Agreement with Laskar over the next ~10,000 years doesn’t validate the model, as both are extrapolations from similar J2000 starting conditions. Significant divergence only appears beyond ~50,000 years.

Short-Term Perturbations

The model describes long-term cycles (thousands to hundreds of thousands of years). A valid question is: how does it handle short-term variations (years to decades)?

What the Model Includes

The model explicitly includes:

What the Model Does NOT Include

Why These Are Omitted

1. The model describes secular (long-term) trends only

The model’s equations produce smooth curves representing mean values over the precession cycles. Short-term oscillations are considered “noise” around these secular trends.

Analogy: Climate models distinguish between:

• Weather (short-term variations, days to weeks)
• Climate (long-term averages, decades to millennia)

Similarly, the Holistic Universe Model describes “orbital climate” (secular trends) not “orbital weather” (short-term perturbations).

2. Short-term effects average to zero

Over timescales longer than the perturbation period, these effects average out:

• The 18.6-year nutation averages to zero over ~40 years
• Jupiter’s ~12-year perturbations average out over ~100 years
• The Chandler wobble averages to zero over ~10 years

3. Practical impact is small

For the model’s primary predictions (obliquity, eccentricity, longitude of perihelion over millennia), short-term perturbations contribute:

• < 0.01° to obliquity predictions
• < 0.0001 to eccentricity predictions
• < 0.1° to longitude predictions

These are within the model’s stated uncertainties.

Comparison with Standard Astronomy

Standard astronomical ephemerides (JPL DE440/441) include all these perturbations through numerical integration. This is why DE440/441 achieves sub-arcsecond precision for short timescales.

The Holistic Universe Model is not designed to compete with ephemerides for short-term predictions. Its purpose is to:

1. Provide a unified framework for understanding long-term precession cycles
2. Make predictions about millennial-scale trends
3. Offer an alternative interpretation of the Fibonacci ratios in these cycles

For precise positions on any given day, use JPL Horizons. For understanding the ~26k-year and ~112k-year cycles, the model provides a conceptual framework.

What This Means for Testability

Short-term tests are difficult because short-term perturbations mask the secular trends. The model’s predictions become testable when:

1. Averaging over sufficient time to remove short-term noise (decades to centuries)
2. Comparing with other long-term models (Laskar, Berger) rather than daily ephemerides
3. Focusing on phenomena where short-term perturbations don’t dominate (e.g., eccentricity minimum timing)

Calibration Transparency

A common criticism of phenomenological models is circular reasoning: if you tune parameters to match data, then cite that match as evidence, you’ve proven nothing. This section explicitly addresses this concern.

The Circularity Problem

The concern (valid):

“If 335,317 was found by fitting to observations, then claiming the model ‘matches observations’ is meaningless. You’ve just done curve-fitting.”

This is a legitimate scientific concern. Any model with adjustable parameters can be made to fit data. The question is: what can the model predict that it wasn’t trained on?

Degrees of Freedom Analysis

The model has 6 free parameters, all governing the Earth simulation. The planetary Fibonacci configuration is uniquely determined by mirror-symmetry constraints (an exhaustive search yields a single solution) and adds no additional degrees of freedom.

For comparison:

• Laskar’s (1993) obliquity formula has ~6 free parameters
• Standard precession theory uses multiple fitted constants

The model has a similar number of free parameters to standard approaches.

What Was Used to Find 335,317?

Direct inputs (the model was explicitly fitted to these):

1. J2000 year lengths - Only H ≈ 335,317 produces solar year (365.242190 days) and sidereal year (365.256363 days) matching observations. See Days & Years for how these derive from obliquity and eccentricity.
2. 1246 AD perihelion-solstice alignment - From Meeus’s formula
3. J2000 longitude of perihelion (102.947°) - The progression from 90° to 102.947° over 754 years
4. J2000 obliquity (23.439°) - For setting mean value
5. Observed obliquity range (22.1° to 24.5°) - For setting amplitude

What Was NOT Used?

Genuine predictions (these values were NOT used in calibration):

Important: The eccentricity and inclination values were checked AFTER 335,317 was determined. They were not used in the fitting process.

The Climate Cycle Question

Problematic: The “eight constraints” in Mathematical Foundations include “Climate Cycles (3 × ~100k pattern)”. This is potentially circular:

1. If ~100k climate cycles were used to find 335,317 → Cannot use them as validation
2. If they were checked afterward → Valid validation

Honest answer: The climate cycle pattern was known to the author when searching for 335,317. It was part of the motivation for the search. However, the model doesn’t match the conventional ~100k; it proposes ~112k. The model’s claim is that conventional chronology has a ~10% systematic error.

This is neither purely input nor purely prediction - it’s a reinterpretation of existing data.

What Would Constitute Independent Validation?

The model can be tested by observations that:

1. Were not used in calibration
2. Cannot be adjusted after the fact
3. Differ meaningfully from standard theory

Strong tests:

Weak tests (similar predictions):

Honest Assessment

What the model CANNOT claim:

• That matching J2000 values validates the model (they were inputs)
• That matching the ~100k climate cycle validates the model (it was known during construction)
• That any parameter fitted to data constitutes evidence

What the model CAN claim:

• Obliquity predictions at dates other than J2000 agree with Laskar
• Perihelion longitude predictions at dates other than 1246 AD/J2000 agree with Meeus
• Eccentricity and inclination emerge from the structure without being used as inputs

Scientific standard: The model should be judged by its testable predictions that differ from standard theory, not by how well it reproduces data used in its construction.

8. Open Questions

The model acknowledges several unresolved questions:

Fundamental Questions

1. Why Do Fibonacci Ratios Appear in Precession Cycles?

The model proposes that Earth’s precession cycles follow Fibonacci ratios (3, 8, 13, 16). While the model fits observed data well, why these ratios appear is not explained by known physics.

The observation:

Inclination precession : Axial precession = ~111,772 : ~25,794 ≈ 4.33 : 1 = 13 : 3
Obliquity cycle : Perihelion precession = ~41,915 : ~20,957 = 2 : 1 = 8 : 4 (Fibonacci-related)

Fibonacci patterns in nature - supporting context:

The model’s author notes that Fibonacci ratios appear throughout nature and the solar system. This is not unique to the Holistic Universe Model:

Fibonacci in the solar system:

A 1968 study by A.M. Molchanov in Icarus analyzed orbital period ratios across the solar system and found that many planetary orbital resonances approximate Fibonacci fractions. More recent work has confirmed this pattern:

Aschwanden (2018)  analyzed 75 orbital period ratios in the solar system (planets and moons) and found that approximately 60% match Fibonacci fractions within measurement uncertainty.

Possible explanations:

1. KAM theorem (strongest explanation): The Kolmogorov–Arnold–Moser theorem (1954–1963) rigorously proves that in perturbed dynamical systems, orbits with “most irrational” frequency ratios are maximally stable against perturbation. 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. Fibonacci ratios (3/2, 5/3, 8/5, 13/8…) converge to φ, so they represent near-maximally stable configurations. The Kirkwood Gaps in the asteroid belt — dramatic depletions at simple integer resonances with Jupiter — are KAM theory in visible action.

2. Observational confirmation: Pletser (2019, Astrophysics and Space Science 364:158) analyzed orbital period ratios in solar planetary and satellite systems and found that ~60% preferentially cluster near Fibonacci fractions (vs ~40% for non-Fibonacci), with these orbits associated with more regular, less inclined, and more circular configurations. Aschwanden & Scholkmann (2017) found Fibonacci harmonic ratios in 73% of 932 exoplanet pairs — extending the pattern well beyond our solar system.

3. Coincidence: The ratios might be approximate coincidences. The human tendency to find patterns (apophenia) may overstate the significance.

4. Selection effect: We observe the current solar system because it’s stable. Unstable configurations would have been disrupted long ago. This doesn’t explain why Fibonacci specifically, but explains why we see stable ratios.

The model’s position: The model notes that if the solar system is a “balanced system” - gravitationally relaxed over 4.5 billion years - Fibonacci ratios may emerge naturally. However, the model does not derive these ratios from first principles; they are empirically fitted.

Scientific status: Fibonacci patterns in orbital mechanics are documented in peer-reviewed literature. The specific claim that Earth’s precession cycles follow exact Fibonacci divisors of 335,317 years is not supported by mainstream astronomy. This remains an open question requiring either:

• A physical derivation from gravitational dynamics
• High-precision measurements confirming the exact periods
• Or demonstration that the apparent pattern is coincidental

2. Why Are the Amplitudes Equal (~0.63603°)?

The model uses the same amplitude (0.63603°) for both:

• Axial tilt variation (obliquity oscillation)
• Orbital inclination variation (relative to the invariable plane)

The observation:

Obliquity range: 22.21° to 24.72° (amplitude 0.63603°)

Inclination range: 0.845° to 2.117° (amplitude 0.63603°)

This equality produces the observed obliquity behavior when both components combine in the model’s formula:

ε(t) = 23.41354° + 0.63603° × [-cos(inclination phase) + cos(obliquity phase)]

Why might this be significant?

1. Conservation principle: Equal amplitudes could indicate energy or angular momentum being exchanged between the two oscillation modes. In coupled oscillator systems, equal amplitudes sometimes emerge from conservation laws.

2. Coincidence: The equality might be approximate and not exact. Current measurements may not be precise enough to detect small differences.

3. Calibration artifact: The model derives both amplitudes by fitting to the same observed obliquity range (22.1° to 24.5°). The equality might be an artifact of this fitting procedure rather than a physical constraint.

Comparison with standard theory:

Standard orbital mechanics calculates obliquity and inclination variations independently:

• Obliquity: 22.1° to 24.5° over ~41,040 years (Laskar)
• Inclination (to ecliptic): ~0° to ~3° over ~100k years

Standard theory does not predict equal amplitudes; the similarity in the model is a feature of its mathematical construction.

Scientific status: The equal amplitude assumption is not derived from physical principles. It is a simplifying assumption that fits current data well but may not hold over longer timescales.

3. Why 335,317 Years Specifically?

The Earth Fundamental Cycle (335,317 years) is an empirically fitted value, not a derived constant.

The fitting process:

The model establishes that year lengths depend on orbital parameters (see Days & Years for details):

• The sidereal year in seconds is fixed (Earth’s orbital period relative to fixed stars)
• Obliquity drives the solar year length
• Eccentricity drives the sidereal year in days
• Day length = sidereal year (seconds) / sidereal year (days)

The 3D simulation calculates year lengths for any Earth Fundamental Cycle value. When testing different values, only H ≈ 335,317 produces year lengths matching J2000 observations:

• Solar year: 365.242190 days
• Sidereal year: 365.256363 days

The Fibonacci connection:

The value 335,317 ≈ ~25,794 × 13, where ~25,794 is the mean axial precession period and 13 is a Fibonacci number. This relationship fits the model’s Fibonacci-fraction pattern for orbital periods, but remains unexplained — it is an observation, not a derivation.

Testable Questions

4. Will Mercury’s “Anomaly” Change?

• Model predicts decrease from ~40 to ~34″/century over ~5,000 years
• Current measurements show no drift (uncertainty ~0.003″/century/decade)
• Requires continued precision observation
• See Mercury Precession for full analysis

5. Will Eccentricity Follow the ~20,957-Year Cycle?

• Model predicts minimum eccentricity (~0.0140) around 11,725 AD
• Laskar predicts continued slow decrease toward ~0.005 over millions of years
• Testable over geological timescales through proxy records
• Near-term (next few millennia): both predictions similar; divergence grows with time

9. References

Primary Sources Used in the Model

1. Precession theory:

   • Capitaine, N., Wallace, P.T., & Chapront, J. (2003). “Expressions for IAU 2000 precession quantities.” Astronomy & Astrophysics, 412, 567-586.

2. Obliquity calculations:

   • Laskar, J., Robutel, P., Joutel, F., et al. (2004). “A long-term numerical solution for the insolation quantities of the Earth.” Astronomy & Astrophysics, 428, 261-285.
   • Laskar, J. (1993). “Orbital, precessional and insolation quantities for the Earth from -20 Myr to +10 Myr.” Astronomy & Astrophysics, 270, 522-533.
   • Vondrák, J., Capitaine, N., & Wallace, P. (2011). “New precession expressions, valid for long time intervals.” Astronomy & Astrophysics, 534, A22. Link 

3. Perihelion calculations:

   • Meeus, J. (1998). Astronomical Algorithms (2nd ed.). Willmann-Bell.

4. Invariable plane:

   • Souami, D., & Souchay, J. (2012). “The solar system’s invariable plane.” Astronomy & Astrophysics, 543, A133.

5. Planetary ephemerides:

   • Park, R.S., et al. (2021). “The JPL Planetary and Lunar Ephemerides DE440 and DE441.” The Astronomical Journal, 161, 105.
   • Fienga, A., Laskar, J., Kuchynka, P., et al. (2011). “The INPOP10a planetary ephemeris and its applications in fundamental physics.” Celestial Mechanics and Dynamical Astronomy, 111, 363-385.
   • Pitjeva, E.V. (2010). “EPM ephemerides and relativity.” Proceedings of the IAU Symposium, 261, 170-178.

Climate and Ice Core References

6. Milankovitch theory:

   • Hays, J.D., Imbrie, J., & Shackleton, N.J. (1976). “Variations in the Earth’s orbit: Pacemaker of the ice ages.” Science, 194, 1121-1132.
   • Berger, A. (1978). “Long-term variations of daily insolation and Quaternary climatic changes.” Journal of the Atmospheric Sciences, 35(12), 2362-2367. Link  — Identified the dominant climatic precession periods at ~23.7, ~22.4, and ~19.0 kyr.
   • Berger, A. (1988). “Milankovitch theory and climate.” Reviews of Geophysics, 26(4), 624-657.

7. 100,000-year problem:

   • Imbrie, J., et al. (1993). “On the structure and origin of major glaciation cycles.” Paleoceanography, 8(6), 699-735.

8. Ice core data and chronology:

   • Petit, J.R., Jouzel, J., Raynaud, D., et al. (1999). “Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica.” Nature, 399(6735), 429-436.
   • Veres, D., et al. (2013). “The Antarctic ice core chronology (AICC2012): an optimized multi-parameter and multi-site dating approach for the last 120 thousand years.” Climate of the Past, 9, 1733-1748. Link 
   • Parrenin, F., et al. (2007). “The EDC3 chronology for the EPICA Dome C ice core.” Climate of the Past, 3, 485-497. Link 
   • Rasmussen, S.O., et al. (2014). “A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records.” Quaternary Science Reviews, 106, 14-28.
   • Kawamura, K., et al. (2007). “Northern Hemisphere forcing of climatic cycles in Antarctica over the past 360,000 years.” Nature, 448, 912-916. Link 

9. Marine sediment chronology:

   • Lisiecki, L.E., & Raymo, M.E. (2005). “A Pliocene-Pleistocene stack of 57 globally distributed benthic δ¹⁸O records.” Paleoceanography, 20, PA1003. Link 

10. Speleothem chronology:

    • Cheng, H., et al. (2016). “The Asian monsoon over the past 640,000 years and ice age terminations.” Science, 352, 343-347.

11. Inclination hypothesis:

    • Muller, R.A., & MacDonald, G.J. (1997). “Spectrum of 100-kyr glacial cycle: Orbital inclination, not eccentricity.” Proc. Natl. Acad. Sci. U.S.A., 94(16), 8329-8334. Link 

12. Milankovitch original work:

    • Milankovitch, M. (1941). Canon of Insolation and the Ice-Age Problem. Royal Serbian Academy Special Publication 132. (English translation: Israel Program for Scientific Translations, 1969)

13. Modern orbital solutions:

    • Laskar, J., Fienga, A., Gastineau, M., & Manche, H. (2011). “La2010: A new orbital solution for the long-term motion of the Earth.” Astronomy & Astrophysics, 532, A89. Link 

14. 100,000-year problem mechanisms:

    • Ridgwell, A.J., Watson, A.J., & Raymo, M.E. (1999). “Is the spectral signature of the 100 kyr glacial cycle consistent with a Milankovitch origin?” Paleoceanography, 14(4), 437-440. Link 
    • Raymo, M.E. (1997). “The timing of major climate terminations.” Paleoceanography, 12(4), 577-585.
    • Paillard, D. (1998). “The timing of Pleistocene glaciations from a simple multiple-state climate model.” Nature, 391, 378-381. Link 
    • Raymo, M.E., Lisiecki, L.E., & Nisancioglu, K.H. (2006). “Plio-Pleistocene Ice Volume, Antarctic Climate, and the Global δ¹⁸O Record.” Science, 313, 492-495.

15. Recent 100-kyr and climate cycle research (2023–2025):

    • Barker, S., Lisiecki, L.E., Knorr, G., Nuber, S., & Tzedakis, P.C. (2025). “Distinct roles for precession, obliquity, and eccentricity in Pleistocene 100-kyr glacial cycles.” Science, 387(6737), eadp3491. DOI 
    • Mitsui, T., Ditlevsen, P., Boers, N., & Crucifix, M. (2025). “100 kyr ice age cycles as a timescale-matching problem.” Earth System Dynamics, 16, 1569–1584. DOI 
    • Lisiecki, L.E. (2023). “Precession pacing of Late Pleistocene ice-sheet changes.” Nature Geoscience.

16. Day length dynamics:

    • Mitchell, R.N., & Kirscher, U. (2023). “Mid-Proterozoic day length stalled by tidal resonance.” Nature Geoscience, 16, 567. Link 

Mercury Perihelion References

17. Historical:

    • Le Verrier, U.J. (1859). “Lettre de M. Le Verrier à M. Faye sur la théorie de Mercure.” Comptes Rendus, 49, 379-383.
    • Newcomb, S. (1882). Astronomical Papers of the American Ephemeris, Vol. 1.
    • Clemence, G.M. (1947). “The Relativity Effect in Planetary Motions.” Reviews of Modern Physics, 19, 361.
    • Berche, B. & Medina, E. (2024). “The advance of Mercury’s perihelion: a historical review.” arXiv:2402.04643. arXiv Link  — Comprehensive historical review reproducing Clemence’s full breakdown table.

18. Modern measurements and precession standards:

    • Park, R.S., et al. (2017). “Precession of Mercury’s Perihelion from Ranging to the MESSENGER Spacecraft.” The Astronomical Journal, 153, 121. Link  — Full PDF (MIT) . Note: estimates Mercury’s orbit jointly with all planets, 343 asteroids, and the PPN parameter β; reports the result in ICRF, not as a directly-measured geocentric perihelion advance.
    • Park, R.S., Folkner, W.M., Williams, J.G., & Boggs, D.H. (2021). “The JPL Planetary and Lunar Ephemerides DE440 and DE441.” The Astronomical Journal, 161, 105. Link  — the GR-inclusive global fit underlying modern Mercury-perihelion analyses.
    • Pitjeva, E.V., & Pitjev, N.P. (2013). “Relativistic effects and dark matter in the Solar system from observations of planets and spacecraft.” Monthly Notices of the Royal Astronomical Society, 432, 3431-3437.
    • Fienga, A., Laskar, J., Kuchynka, P., et al. (2011). “The INPOP10a planetary ephemeris and its applications in fundamental physics.” Celestial Mechanics and Dynamical Astronomy, 111, 363-385. — Independent French ephemeris confirming planetary orbital parameters.
    • Smulsky, J.J. (2011). “New Components of the Mercury’s Perihelion Precession.” Natural Science, 3(4), 268-274. doi:10.4236/ns.2011.34034  — Independent N-body integration via Galactica program yielding geocentric total of 5,601.9″/century.
    • Hilton, J.L., et al. (2006). “Report of the International Astronomical Union Division I Working Group on Precession and the Ecliptic.” Celestial Mechanics and Dynamical Astronomy, 94, 351-367. doi:10.1007/s10569-006-0001-2  — Defines the IAU 2006 general precession rate of 5,028.796″/century.

19. Academic critiques and alternative derivations:

    • Křížek, M., & Somer, L. (2023). Mathematical Aspects of Paradoxes in Cosmology. Springer. Link 
    • Křížek, M. (2015). “On the Perihelion Precession.” PDF  — Contains the 96 km/year calculation.
    • Křížek, M. (2019). “Numerical Modeling of the Anomalous Perihelion Precession of Mercury.” Astronomical Journal of Bulgaria, 27. PDF 
    • Vankov, A.A. (2010). “General Relativity Problem of Mercury’s Perihelion Advance Revisited.” arXiv:1008.1811. arXiv Link  — Contains alternative velocity-based GR formula.
    • Nguyen, A.K. (2024). “Einstein’s Spacetime Curvature Claim Belied By One Second Loophole Of His Own Perihelion Precession Equation.” viXra:2402.0138. PDF  — Analysis of one-second sampling inconsistency in GR perihelion equations.
    • Gerber, P. (1898). “Die räumliche und zeitliche Ausbreitung der Gravitation” (The Spatial and Temporal Propagation of Gravity). Zeitschrift für Mathematik und Physik, 43, 93-104. Wikipedia  — Published the same perihelion precession formula 17 years before Einstein.
    • Corda, C. (2023). “On the existence of precession of planets’ orbits in Newtonian gravity.” Qeios. Link  — Notes that the Solar System barycenter shifts ~1000 km/day, much larger than Mercury’s 96 km/year perihelion shift.

20. Solar gravitational quadrupole:

    • Mecheri, R., & Abdelatif, T. (2022). “Secular Variations of the Sun’s Gravitational Quadrupole Moment and Their Impact on GR Tests.” Remote Sensing, 14(19), 4798.

21. Solar oblateness and BepiColombo:

    • (2022). “The Influence of Dynamic Solar Oblateness on Tracking Data Analysis from Planetary Missions.” Remote Sensing, 14(17), 4139. Link 

Fibonacci and Orbital Resonance References

22. DNA and Golden Ratio:

    • Yamagishi, M.E.B., & Shimabukuro, A.I. (2008). “Nucleotide frequencies in human genome and Fibonacci numbers.” Bull. Math. Biol., 70(3), 643-653. PubMed 

23. Phyllotaxis (plant Fibonacci patterns):

    • Prusinkiewicz, P., & Lindenmayer, A. (1990). The Algorithmic Beauty of Plants. Springer. Link 

24. Planetary orbital resonances:

    • Molchanov, A.M. (1968). “The resonant structure of the Solar System.” Icarus, 8(1-3), 203-215. ScienceDirect 

25. Solar system Fibonacci analysis:

    • Aschwanden, M.J. (2018). “Self-organizing systems in planetary physics: Harmonic resonances of planet and moon orbits.” New Astronomy, 58, 107-123. arXiv 

26. Fibonacci in planetary period ratios:

    • Pletser, V. (2019). “Fibonacci Numbers and the Golden Ratio in Biology, Physics, Astrophysics, Chemistry and Technology: A Non-Exhaustive Review.” Astrophysics and Space Science, 364, 158.

Last updated: May 2026

Return to Mathematical Foundations | Predictions | Formulas