Scientific Background

This page provides in-depth scientific context for readers who want to understand how the framework 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 Integer Lattice
2. Physical Mechanisms
3. Comparison with Standard Precession Theory
4. The Mercury Perihelion Question — Resolved
5. Eccentricity Cycles and Milankovitch Theory
6. Mathematical Framework
7. Open Questions
8. References

1. The Integer Lattice

This section’s title is kept for its links; its content is the restatement. The Earth-side layer of the framework is one derived clock and a set of beats: Earth’s mean lunisolar precession period (25,771.4 years at J2000, the beat of the model’s own sidereal and tropical year laws) is the clock, and the other long periods follow from it and two orbital motions read from the model’s N-body chain — the apsidal period (~111,582 yr) and the nodal period (~68,751 yr): 1/T_peri = 1/T_p + 1/T_aps gives the perihelion-of-date period, 1/T_obl = 1/T_p − 1/T_s₃ the obliquity cycle. The model’s small periodic correction terms are Fourier series on a fixed grid whose unit is a fitted anchor of 335,317 years — a bounded harmonic basis under pre-registered falsification tests, not a parameter count. The planets themselves ride the model’s own N-body engine, derived from one cited J2000 state with zero fitted constants (see Mathematical Foundations).

The retired integer lattice and “six Fibonacci Laws”. Earlier versions of the model presented Earth’s periods as small-integer divisions of the anchor (H/3, H/5, H/8, H/13, H/16, with the closures as integer identities), and a set of planet-side laws — amplitude constants ψ and K predicting all eight planets’ inclination and eccentricity amplitudes from divisor and mass, two seven-against-Saturn balance identities, and a mirror-symmetric configuration uniqueness claim. Re-evaluated with the model’s own N-body engine, the Earth periods sit 0.09–2.5 % off the integer fractions with ratios that wander far from them across ±26 kyr, the planets’ fractions landed on the grid at chance level, the balances hold only approximately with dynamical long-term means (~98% for eccentricity), the Saturn e-prediction missed by 3.6%, and the configuration-uniqueness claim proved input-dependent. The framing is retired, including the “gas-giant lock” once read between Jupiter’s ICRF perihelion, Saturn’s ecliptic perihelion and Earth’s obliquity beat — that beat is the clock against Earth’s own nodal mode s₃. The published literature on Fibonacci-like period ratios in planetary systems (Molchanov 1968, Aschwanden 2018, Pletser 2019, the KAM theorem) is context for resonance-chain physics in general, not evidence for that reading. The full re-evaluation record is doc 109 in the 3d repository ; the retired framing’s account is on Earth’s Clock. What the model states today is the composed clock with its measured tidal history, the beats of that clock with the planets’ modes, and the two-tier deep-time prediction those beats make.

2. Physical Mechanisms

The framework’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 2006) and the model’s own J2000 reading 25,771.4 years — the beat of its certified sidereal and tropical year laws — agree to eight parts in a million. The period is a reading, not a constant: it wanders between 24,801 and 26,559 years across ±26 kyr as the orbit modulates the torques, and it lengthens at geological time (Earth’s Clock). Throughout this document, ~25,771 years refers to the IAU value; ~25,771 years to the model’s J2000 reading.

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: Perihelion Precession

When axial and apsidal precession combine, they produce perihelion precession — the mean meeting cycle that determines when Earth is closest to the Sun relative to the seasons:

Perihelion precession period = 1 / (1/T_axial + 1/T_apsidal)
                           = 1 / (1/~25,771 + 1/~111,582)
                           ≈ ~20,936 years

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

This ~21k-year cycle is the formula’s single-period result. The actual observed Milankovitch climate signal — what is conventionally called climatic precession — sits at a different period (~23.7 kyr dominant peak, with additional peaks at ~22.4 and ~19.0 kyr; Berger 1978). The two are distinct: the simple beat formula uses one apsidal precession value (~112 kyr) and yields a single ~21-kyr mean, but apsidal precession actually has internal eigenmode structure (planet-specific g_j sub-modes) that splits the observed signal into the multi-peak Berger spectrum centred near 23.7 kyr. See §5 for the spectral details.

The Model’s Representation

The framework represents these same physical phenomena using a different mathematical formulation:

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 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 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 the statements differ in kind rather than in the next few millennia:

Earth’s deep eccentricity minimum is not a divergence: the model’s own N-body engine forecasts the La2004-class deep minimum, the apsidal-period law being its epoch-local chart (see Eccentricity).

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 column below is the scene device’s harmonic obliquity law, kept here for the comparison it was built for — the published obliquity is the derived one, from a single precession equation on the engine’s own orbit plane, tracking La2004 to 50″ rms over 13 kyr (Obliquity). It agrees with La2004 to within ~0.001° from -1000 BC through 5,000 AD — essentially exact over this 6 kyr window — and within 0.1° from 7,000 AD through 12,000 AD. The largest near-term discrepancy (~0.37°) occurs at 10,000 BC, and by 20,000 AD the device and La2004 differ by ~0.29°, with La2004 oscillating back upward while the device’s bounded harmonic continues. The model’s published obliquity is the derived one — the spin axis integrated on its own N-body orbit plane — which tracks La2004 to 50″ rms over the last 13 kyr and 687″ rms over the last million years (Obliquity).

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)

The model’s eccentricity is its own N-body series, so this comparison is two integrations of the same secular system (theory-vs-theory):

Assessment: Both decline together through the coming minimum near 27,000 AD, and the two integrations differ only in the last digit shown: the engine’s Venus-driven mode g₂ reads 7.4524″/yr against La2004’s 7.452″/yr, with the Moon’s quadrupole term carried in the integration; the minimum’s depth is a consistency check, not a model-vs-standard discriminator (see Eccentricity).

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:

• Four named structural constants on the Earth spin-and-time side plus a gated fitted-correction stack; zero free parameters on the planetary orbits, with the per-planet long-period descriptors accounted in the device ledger (the three-ledger accounting of Mathematical Foundations)
• The Earth-side correction terms are a bounded harmonic basis, not a derivation from celestial mechanics; the planets are N-body dynamics

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 is narrower than it once was: both are N-body integrations with the first-order relativistic correction. DE441 carries the full minor-planet population and is fitted to the observational record; the model integrates nine bodies (the Moon’s quadrupole term and the three largest asteroids included) from one cited J2000 state with no fitted constants on that path.

4. The Mercury Perihelion Question — RESOLVED to General Relativity

This section examines what was one of the most debated aspects of the framework: the once-proposed 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 determinations measure the same ecliptic-plane motion — the difference is the reference direction. The ~575″ value is relative to fixed stars (ICRF), the inertial frame defined by distant quasars. The classical 5,599.74″ value is relative to the moving vernal equinox, which the equinox precession carries backward — it is what was historically measured before ICRF existed. Each is a complete determination inside its own constant system: the classical total contains Newcomb’s precession constant (5,025.645″/century), not the modern IAU 2006 one, so the two totals are not inter-convertible by adding the modern precession rate to the ICRF value — the anomaly is extracted by subtraction within one system, where the equinox precession cancels.

The Classical Geocentric Total

The commonly cited “~5,600″/century” geocentric total is Clemence’s (1947)  equinox-of-date determination, built on Newcomb’s equinox precession rate of 5,025.645″/century. The full Clemence breakdown (Berche & Medina, 2024 , Table 2):

Newcomb’s equinox precession (5,025.645″) has since been superseded — the IAU 2006 precession model (P03)  gives 5,028.796″/century (Lieske 1976: 5,029.097″) — but that does not license “updating” Clemence’s total by swapping constants: his 5,599.74″ is a joint solution of 1765–1937 longitudes in which the precession constant, the planetary masses and the orbit corrections were determined together. No equinox-referred determination exists in the modern system, and one cannot be manufactured by adding the modern precession rate to the ICRF ranging value — the anomaly is only ever formed by subtraction within a single system.

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:

Smulsky’s N-body integration (5,601.9″) and Berche & Medina’s analytical review (5,599.7″) agree on the geocentric total at J2000 while attributing the ~43″ surplus over the Newtonian planetary perturbations to different causes (a compound solar rotation, and General Relativity). Smulsky’s results also show the geocentric total decreasing slightly between epochs (5,602.9″ at 1950 → 5,601.9″ at 2000).

The framework’s contribution to this comparison is an identity rather than a third total: Mercury’s ecliptic perihelion advance in the model is the ecliptic rate 531.44″/century, and projected into right ascension at the J2000 perihelion longitude and the IAU 2006 obliquity it reads 574.14″/century — an excess of 42.71″ against the relativistic advance of 42.98″ derived from the same constants. The identity is stated for all planets, reproduces the anomaly for Mercury only, and is closed as a coincidence — see the settled statement at Mercury Precession.

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: Three Layers, Answered

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

Three layers of processing separate the raw measurement from the reported 575″/century, and each was once an open question for the model. The model’s own bookkeeping has since answered all three (methodology at Mercury Precession; the measured decomposition is the simulator’s Cycles tab and the closure gate of doc 13 §1.8):

1. Reference frame transformation. The model does carry a frame term of the anomaly’s size. Its own Earth-frame rate for Mercury’s perihelion — the right ascension of the perihelion direction in the scene’s equator — is 578.79″/century, and it decomposes exactly: the ecliptic advance 531.44″ + the equatorial projection 42.71″ (the RA scale is stretched by dα/dλ − 1 where the perihelion sits) + the obliquity-rate term 4.31″ + the scene’s own marker conventions (~1.4″). The projection term is a coordinate quantity and nothing else: it is identically zero in ecliptic longitude, it does not depend on Earth’s orbital motion (measured: six different eccentricity histories leave the rate unchanged to 0.01″), and it swings between −44″ and +48″ over Mercury’s 244-kyr perihelion cycle, changing sign every 61 kyr.

2. Newtonian subtraction. The Newtonian ~532″ is not a long-term mean but the instantaneous rate computed from the planetary masses for the present configuration (Le Verrier 1859; Clemence 1947: Venus ≈ 278″, Jupiter ≈ 154″, Earth ≈ 90″, Saturn ≈ 7″, Mars ≈ 2.5″). The beating of the secular eigenmodes — the “sometimes faster, sometimes slower” of a perihelion that is the phase of a sum of rotating vectors — is already inside that number. It cannot absorb a further 43″ without changing Venus’s or Jupiter’s mass, and those are fixed to 10⁻⁸ by their satellites and spacecraft flybys.

3. Which coordinate the observations are in. The classical chains reduced meridian positions to ecliptic longitude before any orbit was fitted (Le Verrier’s 187 equations; Clemence’s 5,599.7″ against the equinox of date), and the ranging chains (Park et al. 2017 , Pireaux & Rozelot 2003: 42.980 ± 0.002″) fit the orbit in ICRF Cartesian coordinates. Neither uses the equatorial projection, and both contain the 43″. The transit record (1677–2019, alignments that no coordinate choice can move) runs at ~574″, not 531″.

What this means for the model’s argument: the ~42.98″ is present in every equator-free measurement, while the model’s 42.71″ is present only in the equatorial coordinate. The two agree to 0.6 % — and by 0.27″/century they disagree, at roughly fifty times the ranging uncertainty. The same projection rule gives −1.9″ for Venus, +493″ for Earth and −101″ for Mars where +8.62″, +3.84″ and +1.35″ are measured. The identity is a characterised coincidence (the exact match falls at a perihelion longitude of 77.9°, which Mercury reaches around the year 2250), not the cause. The frame hypothesis for Mercury is therefore closed by the model’s own measurements; what survives is the exact decomposition itself, which is now a gate.

Academic Critiques, Briefly

A short record of the critiques the literature carries, none of which changes the settled attribution: Křížek (2023 ) notes the anomaly is a small difference of two large numbers (~96 km/yr of perihelion motion at Mercury’s distance) — but MESSENGER’s ~0.8 m ranging precision resolves 96 km/yr about 120,000× over, so smallness is no refuge. Gerber (1898) published a formula returning the same ~43″/cy from finite gravity-propagation speed, 17 years before Einstein — the consensus (von Laue) is that his derivation was flawed and the numerical agreement coincidental. The historically interesting circularity concern (modern ephemerides carry GR inside the fit) does not touch the anomaly’s existence, which Le Verrier established in 1859, fifty-six years before General Relativity.

The Record of the Projection Reading

The once-proposed alternative held that the ~43″ arises from the equatorial projection of Mercury’s ecliptic advance — a coordinate identity using three inputs (Mercury’s device ecliptic period, the IAU J2000 perihelion longitude, the IAU 2006 obliquity): the slope dα/dλ at Mercury’s perihelion longitude is 1.08036, so the 531.44″ ecliptic advance reads 574.14″ in right ascension — an excess of 42.71″ against the relativistic 42.98″. The reading was closed by its own tests: the ranging precision (42.980 ± 0.002″ leaves the 0.27″ shortfall at >100σ), the all-planet failure of the identical projection, the absence of any equatorial step in either determination chain, and the model’s own N-body engine requiring the 1PN term to reproduce the observed window rate. What survives, as a gate rather than a claim: the model’s Earth-frame rates for all seven planets are exact projections of their ecliptic advances (the closure gate pins the decomposition to 1″/cy at 1900/2000/2100) — real coordinate bookkeeping, needed to read the simulator’s equatorial-frame exports, and unrelated to the anomaly.

The Apparent-Fluctuation Machinery

Distinct from the anomaly and fully alive: each planet’s apparent Earth-frame precession rate oscillates while its ecliptic advance stays at the ecliptic rate. The oscillation is the equatorial projection of the ecliptic motion — one closed formula with no fitted constants, Earth’s obliquity of date the only Earth input — so it is a reference-frame effect calculable for every planet from its period and J2000 perihelion longitude alone, no observation of the planet required. The range scales with the ecliptic rate: Saturn, on the shortest perihelion period, swings widest; Neptune, on the longest, least. The formula: Formulas §7; the per-planet ecliptic/Earth-frame table: Mercury Precession.

Where the Model Stands

The model’s position, in one paragraph: General Relativity supplies the ~43″/century, and the model carries it as a derived, zero-fitted supplement — the 1PN term computed from the same constants (GM☉, c, a, e) and gated against the model’s own N-body engine, giving the stated total of 574.4″/century. The once-proposed projection alternative is closed (see the record above); what the model adds beyond the standard account is the apparent-fluctuation machinery — the reference-frame layer that any Earth-based reading of a planetary rate must pass through, predicted for all seven planets from time alone. The full settled statement lives on the Mercury Precession page.

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 Eccentricity: the N-body Series

The model’s eccentricity is the N-body series on every path — Earth’s eccentricity vector integrated from one cited J2000 state, which carries the 405-kyr metronome (405.6 kyr in the ±10-Myr run, rock record 405.6) with its 123.8- and 94.9-kyr companions and follows the La2004-class deep minimum (Supporting Evidence §14). The published value, the Sun’s equation of centre behind the eclipse and cardinal-point instants and the Moon’s eccentricity channel read this one source, anchored at J2000 to the observed eccentricity and its rate (−4.2037e-5/century). The perihelion direction is a separate quantity on the ~20,936-year perihelion-of-date cycle (December-solstice alignment ~1246.03125 AD). Full statement and epoch-by-epoch comparison: Eccentricity.

The Standing Challenges, Briefly

• Laskar’s computed range (~0.0005–0.058). The model’s series spans the same class of range: it is the same secular system integrated from the same kind of seed, reproducing the leading eigenfrequency g₅ to 4.2574 ″/yr against Laskar’s 4.2575. What the model rejects is the attribution of the geological ~100k climate signal to eccentricity forcing: adding Laskar’s own e(t) and ϖ(t) to the climate formula adds nothing (the insolation null test), and the tuning-independent speleothem centroid sits on the formula’s planetary beats, not on eccentricity’s 95k/125k/405k structure.
• No standard mechanism links spin-axis orientation to orbital shape. None is claimed: the 1246 AD perihelion–solstice alignment constrains the perihelion direction (the perihelion-of-date cycle); the orbit’s shape is the planets’ secular dynamics. Separate quantities, separate cycles.
• Berger’s multi-peak climatic-precession spectrum (~19–24 kyr) vs one perihelion-of-date period. The formula carries Berger’s peaks as its five climatic-precession lines p+gᵢ (19.0–23.8 kyr) — the clock rate plus each planetary apsidal mode, read from the model’s own integration; the perihelion-of-date period (~20,936 yr) is the beat with Earth’s own apsidal motion and sits near the centroid of the spread. (The comb-era integer reading of the same peaks is the record in Supporting Evidence §7.)

Model vs. Laskar at a Glance

Where it stands: the 100,000-year problem remains unsolved in the mainstream (Barker 2025, Mitsui 2025, Lisiecki 2023). At the current record length (T ≈ 1.2 Myr) the Rayleigh resolution near 110 kyr is ~10 kyr, so 95k, 100k and 112k are spectrally collinear; the attribution therefore rests on the discriminating tests rather than on the peak position. Two facts sit awkwardly for direct Earth-eccentricity forcing: the ~100k climate signal lacks eccentricity’s split-peak structure (Muller & MacDonald 1997 ), and eccentricity’s theoretically strongest ~405-kyr component is largely absent from the last 1.2 Myr of climate records; the U-Th-dated Cheng2016 speleothem record rules out a ~10% chronology offset as the explanation. The full empirical case, including the ice-core dating methods: Supporting Evidence §1; epoch-by-epoch numbers: Eccentricity: Numerical Comparison.

6. Mathematical Framework

Model Geometry and Key Periods, in Brief

The simulation works in ICRS-aligned coordinates (barycentric origin, J2000 equinox and pole) and represents the two precessions with two reference points. EARTH-WOBBLE-CENTER: Earth orbits it clockwise in ~25,771 years; it marks the precession direction and does not enter e(t). PERIHELION-OF-EARTH: orbits the Sun counter-clockwise in ~111,582 years; it marks a direction, not a measured distance — Earth’s eccentricity is the N-body series (Eccentricity). Its angular position is Earth’s longitude of perihelion.

One clock governs the Earth chain — the mean lunisolar precession period, 25,771.4 years at J2000, derived from the year laws and lengthening at geological time (see Expanding Resonance). The correction bases’ unit is the fitted anchor of 335,317 years, determined empirically (see Calibration Transparency below). Every Earth cycle is a dynamical value read against the clock:

The perihelion-of-date period is the frame-arithmetic addition identity: the perihelion direction turns at the meeting frequency of the two counter-rotating motions, 1/T_peri = 1/T_p + 1/T_aps. The scene device’s equations — its harmonic obliquity law, the inclination cosine, the perihelion-longitude progression — are stated with all constants in Formulas; the published obliquity, inclination and eccentricity are the derived quantities of Ledger 1 (Mathematical Foundations), and the retired integer reading of the periods is recorded on Earth’s Clock.

Comparison with Standard Formulas

Key observation: The model matches current observed values well; at deep time its N-body series follows the same secular dynamics Laskar integrates, while its spin clock lengthens on the measured tidal history. Both approaches are theoretical extrapolations that cannot be directly verified for ancient/future periods.

Error Analysis, in Brief

The scene device’s uncertainty budget has four sources. The anchor year (t₀ = -302,635) is fitted to J2000 values, not independently constrained — it shifts the phase of all cycles but no period. The device’s amplitude (A = 0.63607°, ±0.01°) derives from the observed obliquity range (~22.1° to ~24.5°); its combined range is 22.21°–24.72°. The mean values (mean obliquity 23.41353°, base eccentricity 0.015386) carry sub-0.001-class uncertainty. The periods are the device’s series means, which differ from the standard instantaneous values by construction (obliquity ~41,224 vs Berger’s ~41k yr). Propagated: at 3000 AD, obliquity ±0.02° and eccentricity ±0.0001; at 12,000 AD, ±0.2° and ±0.001 (full year-by-year comparison: §3 — Obliquity Predictions). One caveat governs all of it: agreement with Laskar over the next ~10,000 years doesn’t validate the model — both are extrapolations from the same J2000 starting conditions, and significant divergence only appears beyond ~50,000 years.

Short-Term Perturbations

The model is secular by design: it carries the five long cycles (axial, apsidal, ecliptic, obliquity, perihelion-direction) and deliberately averages out everything faster — the 18.6-year lunar nutation (±9″), the Chandler wobble, Jupiter’s and Saturn’s ~12/~29.5-year orbital perturbations. Over millennial predictions these contribute <0.01° to obliquity, <0.0001 to eccentricity, <0.1° to longitude — within the stated uncertainties. The model is not designed to compete with JPL DE440/441 for short-term positions (use Horizons for those); it becomes testable where the secular trends dominate — averaged over decades to centuries, compared against other long-term solutions (Laskar, Berger), or at phenomena like the timing of the eccentricity minimum.

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’s accounting is three ledgers, not a single parameter count (Mathematical Foundations): zero free parameters on the planetary side (one cited J2000 state), four named structural constants on the Earth spin-and-time side, and a gated fitted-correction stack. (The comb era stated “6 free parameters” with a planet configuration fixed by a mirror-symmetry search; the search’s uniqueness claim proved input-dependent and is retired.)

For comparison:

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

The model publishes no single free-parameter count: its inputs sit in three ledgers, and the planetary ledger has none at all — the orbits come from one cited J2000 state (Mathematical Foundations). The fitted constants live in the time-domain correction stack.

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 a climate-cycle compatibility requirement — the post-MPT ~100-kyr band. 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 the anchor, and it was part of the motivation for the search. The model no longer attributes the ~100-kyr band to Earth’s apsidal period: the shipped climate formula carries the band as the planets’ eccentricity beats read from the model’s own modes, and Earth’s own apsidal period is not a line of the formula (Climate Formula). The anchor’s coincidence with the ~100-kyr band is therefore neither an input to the present formula nor a prediction of it; the formula’s own status as description versus prediction is measured separately (T7, on the Climate Formula page).

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.

7. Open Questions

The model acknowledges several unresolved questions:

Fundamental Questions

1. Were the Obliquity Amplitudes Equal? — retired

The model’s former obliquity device used one amplitude (0.63607°) for two cosine components — an “axial tilt” term on the obliquity beat and an “inclination tilt” term on the apsidal period — and this section asked whether that equality was physics, coincidence or a calibration artifact. The question is closed: the equality was an artifact of the construction (both cosines were fitted to the same observed obliquity range, ~22.1° to ~24.5°), and the model no longer decomposes the obliquity at all.

The present state: the published obliquity is derived by integrating the spin axis’s precession about Earth’s orbit normal, with the orbit normal read from the model’s own N-body nodal modes and one anchor, the J2000 precession rate. It reproduces La2004 to 50″ rms over the last 13,000 years with no fitted constants, and its dominant period is the beat of the precession clock against the nodal mode s₃ (~41,224 years). Earth’s orbital inclination to the invariable plane is a separate quantity on the nodal modes, not a component of the obliquity.

Comparison with standard theory: this is the same construction standard secular theory uses — obliquity from the spin axis against a perturbed orbit plane (Laskar: ~22.1° to ~24.5° over ~41,040 years). The model’s distinctive statement is not a decomposition but the clock the beat rides on at deep time (Expanding Resonance).

2. Why 335,317 Years Specifically?

The fitted anchor (335,317 years, the unit of the correction bases) 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 the J2000 anchor (Earth’s orbital period relative to fixed stars; held constant within the modern-era scope, drifts at deep time per Expanding Resonance)
• 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 anchor 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 relation to the clock:

The anchor is 13.011 × ~25,771 years, the J2000 axial precession period — not 13 exactly: the anchor was fitted on the perihelion-of-date alignment, not on the axial rate, so the 0.086 % gap is the fit’s convention. The two scale together at deep time and their ratio is a fit constant, not structure; the former reading “H = 13 × axial precession, 13 a Fibonacci number” is retired.

Testable Questions

3. Will Mercury’s “Anomaly” Change? — RESOLVED: no

• Resolved with the frame reading’s closure (§4 above): the model now attributes the excess to General Relativity and predicts it constant, as GR does
• Current measurements show no drift (uncertainty ~0.003″/century/decade) — consistent with the settled position
• See Mercury Precession for the settled statement and the record

4. Does the Model’s Eccentricity Follow the Secular Solutions?

• The model’s own N-body series and La2004 decline together into the same minimum region near ~27,000 AD (0.00265 vs 0.00263, see §5) — two integrations of one secular system
• The discriminator at deep time is not the secular solutions but the rock record: the long-eccentricity metronome (405.6 kyr), which the model’s own gravity reproduces at 405.6 kyr
• The Moon’s part is measured: with its quadrupole term in the integration the engine’s g₂ reads 7.4524″/yr against La2004’s 7.452″/yr; with the Earth–Moon pair merged into a single point mass the mode reads low and the beat long

8. 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.
   • Chapront, J., Chapront-Touzé, M., & Francou, G. (2002). “A new determination of lunar orbital parameters, precession constant and tidal acceleration from LLR measurements.” Astronomy & Astrophysics, 387, 700-709. Link 

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., Joutel, F., & Boudin, F. (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.” Nature, 534(7609), 640-646.

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. — Mercury precession refinement to ~43″/century.
    • Newcomb, S. (1898). Tables of the Motion of the Earth on its Axis and Around the Sun. Astronomical Papers Prepared for the Use of the American Ephemeris and Nautical Almanac, Vol. VI, Part I (Washington: Bureau of Equipment, Navy Department). — Source of the classical eccentricity polynomial.
    • 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.” European Journal of Physics, 45, 055601. DOI: 10.1088/1361-6404/ad48ca  · arXiv:2402.04643  — 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 

Return to Mathematical Foundations | Predictions | Formulas