Scientific Background

This page places the framework beside established astronomical theory. It names the physical mechanisms the model’s two counter-rotating motions stand for, compares the model’s values with IAU 2006, Laskar’s secular solutions and the JPL ephemeris epoch by epoch, records how the Mercury perihelion question was posed and settled, and gives Milankovitch theory a fair hearing before stating where the model’s climate attribution differs. The model’s own construction — how the clock is derived, how the anchor was fitted, and the parameter accounting — is on Mathematical Foundations; the concepts are on the Model pages.

This page establishes:

• The two motions are standard axial and apsidal precession, with known causes and opposite directions — Physical Mechanisms.
• Near J2000 the model matches IAU 2006, La2004 and DE440; the divergences are at deep time and in kind — Comparison with Standard Precession Theory.
• Mercury’s ~43″/century anomaly is General Relativity’s; the projection reading was tested and closed — The Mercury Perihelion Question.
• The model’s eccentricity is the same secular dynamics Laskar integrates; what it disputes is the 100-kyr climate attribution — Eccentricity Cycles and Milankovitch Theory.
• The model’s geometry in brief, and the questions the model’s own engine has settled — Model Geometry in Brief, Resolved and Open Questions.

Terms used on this page

1. 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,635)
                           ≈ ~20,938 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 §4 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.

2. 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 and Diverge

For periods of ±2,000 years around the present, the model closely matches established theory: obliquity within ±0.01° of Laskar (1993) and Chapront et al. (2002), longitude of perihelion within ±0.1° of Meeus (1998), the precession rate within 0.01 % of IAU. The quantity-by-quantity ledger — what the model reproduces at J2000 against the IAU, JPL and Souami & Souchay references, and the two topics on which it proposes a different explanation — has one home: Supporting Evidence: Scope. For longer timescales the statements differ in kind rather than in the next few millennia. The axial-precession period is, in the model, a reading that wanders between 24,801 and 26,559 years across ±26 kyr, and lengthens at geological time on the composed lunisolar rate (Expanding Resonance). The standard long-term expressions are instead a polynomial in time (Capitaine) or Vondrák 2011’s Fourier series to ±200 kyr. The obliquity band is the beat 2π/(ψ̇(t) − |s₃|) of the lengthening clock against the nodal mode, tracking La2004 to 13.5″ rms over the last million years (Obliquity), where polynomial expressions diverge beyond a few tens of kyr. Earth’s deep eccentricity minimum is not a divergence: the model’s own N-body engine forecasts the La2004-class deep minimum (Eccentricity). Mercury’s ~43″ anomaly is General Relativity’s, carried as a derived supplement (§3 below); Saturn’s window-era retrograde perihelion is the Great-Inequality phase, confirmed by the model’s own engine (Supporting Evidence §9). Both the model and the secular solutions are theoretical extrapolations that cannot be directly verified for ancient or future epochs.

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 (Mathematical Foundations: calibration vs prediction).

Obliquity Predictions (Model vs La2004)

Assessment: The Model column 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 1.2″ rms over the last 100 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 1.2″ rms over the last 100 kyr and 13.5″ rms over the last million years (Obliquity).

The same device law against a second reference, Chapront et al. (2002):

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).

Anyone can verify these comparisons against JPL data; the sources and the steps are listed under Mathematical Foundations: reproducibility.

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
• 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.

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.

3. The Mercury Perihelion Question — RESOLVED to General Relativity

This section places Mercury’s perihelion precession in its literature: the history of the anomaly, the classical equinox-of-date determination against the modern ICRF one, the standard explanation, and the critiques the literature carries.

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 Classical Geocentric Total

The modern ranging determination (575.31 ± 0.0015″/century, Park et al. 2017) is relative to the fixed stars (ICRF); the classical determination is relative to the moving equinox of date. Both measure the same ecliptic-plane motion, each is a complete determination inside its own constant system, and the anomaly is extracted by subtraction within one system, where the equinox precession cancels (Mercury Precession sets the two columns side by side). 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 adds no third total: its Mercury rate is the Newtonian ecliptic advance plus the relativistic supplement derived from the same constants (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

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.

4. 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,938-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,938 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.

5. Model Geometry 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,635 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). Every Earth cycle is a dynamical value read against it. 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 five periods with their derivations are tabulated on Mathematical Foundations; the correction bases’ unit, the fitted anchor of 335,317 years, and how it was found are there too (the clock and the anchor). 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.

6. Resolved and Open Questions

1. Why 335,317 Years Specifically?

The fitted anchor is an empirically fitted value, not a derived constant: the eight constraints it satisfies, the fitting process and its relation to the clock (13.011 clock periods — a fit convention, not structure) are on Mathematical Foundations. The model does not explain why that number rather than another; it is the fit’s result on the 1246 AD alignment.

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

• Resolved with the frame reading’s closure (§3 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

3. 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 §4) — 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

7. 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