Perihelion Precession

Every planet’s perihelion (closest approach to the Sun) slowly rotates around the Sun — perihelion precession. Standard celestial mechanics attributes this to gravitational perturbations from other planets. The Holistic Universe Model proposes that part of the apparent rate is also a reference-frame effect from Earth’s own motion.

A perihelion rate lives in one of two coordinates, and the two must not be compared with each other: (a) ecliptic longitude — what every observer publishes (relative to the stars, or, adding the equinox precession, of date); (b) right ascension in an equatorial frame — what the simulation’s Earth-frame export measures, which no observer has ever published. The table below keeps them apart: the WebGeocalc trend, Fitzpatrick’s secular value (Standish & Williams 1992) and the Lagrange–Laplace approximation are all (a); the model’s lattice rate is (a); the model’s Earth-frame value is (b), shown for completeness only.

WebGeocalc = JPL ephemeris 1800–2000 trend, ecliptic longitude (a). Model lattice rate = 360° per the planet’s 8H/N period, ecliptic (a). Model Earth-frame RA = the lattice rate projected into the simulation’s equatorial frame plus the obliquity-rate term (b) — the quantity the Earth-frame export and the predictive formula produce; it is not an observable and is not fit to WebGeocalc. Fitzpatrick = long-term secular average (Standish & Williams 1992). L-L Theory = analytical first-order approximation (Lagrange–Laplace, 19th century).

Heads-up. Mercury’s WebGeocalc 1800–2000 trend (~572″/cy) is ~3″/cy below the canonical ~575″/cy used in the textbook anomaly story below. The standard “anomaly = 575 − 532 = 43″/cy” equation needs the higher figure. The ~575 comes from MESSENGER’s 2013 spacecraft snapshot at J2000; whether that’s a stable long-term rate is exactly what BepiColombo (2027) will test.

The Mercury “Anomaly”

Mercury’s perihelion precession is historically significant because of a famous discrepancy debated for over a century.

Both columns describe the same physical motion in the ecliptic plane — the difference is the reference direction. The right column (geocentric) is what’s directly measured on Earth, against the moving vernal equinox which drifts backward at ~5,028.8″/century due to Earth’s axial precession. The left column (ICRF) subtracts that drift, measuring against fixed stars. The equinox drift cancels in subtraction, so the ~43″ discrepancy is the same in both frames.

The equinox-based “~5,600″” figure dates from Clemence (1947), who used the then-current 5,025″ equinox precession rate. Modern measurements give ~5,028.8″/cy — used throughout this page — so the corresponding total is ~5,604″/cy.

Origin of the ~532″ Newtonian prediction

The ~532 arcseconds/century Newtonian prediction comes from gravitational perturbations by all other planets. Since Mercury is the innermost planet, every other planet pulls its perihelion forward (prograde):

These contributions are calculated by N-body integration of Newton’s inverse-square law over time — typically using JPL’s DE-series planetary ephemeris, not a closed-form equation. Le Verrier (1859, ~527″/century) first flagged the discrepancy with observation; Newcomb (1882, ~532″/century) refined it to the canonical figure underlying the ~43″/century GR anomaly. Newton himself never computed this — Mercury’s precession was not measured precisely until well after his death (1727).

The standard explanation

The standard explanation for the 43 arcsecond discrepancy is Einstein’s General Relativity (1915): space-time is curved near massive objects; Mercury, closest to the Sun, experiences the strongest curvature, producing an additional precession of ~43 arcseconds/century given by Δϖ_GR = 6πGM / (ac²(1−e²)) per orbit. This was one of the first major confirmations of Einstein’s theory.

The Model’s Alternative Explanation

The projection

In the model Mercury’s perihelion advances in ecliptic longitude at exactly the lattice rate, 360° per 8H/11 = 531.44″/century — numerically the Newtonian planetary-perturbation value. An observer on Earth does not measure that longitude directly: the sky is measured in right ascension, against the equator, which is tilted to the ecliptic by the obliquity ε. A direction advancing uniformly in ecliptic longitude λ advances in right ascension α at the rate

dα/dλ = cos ε / (cos²λ + sin²λ cos²ε)

which is larger than 1 where the perihelion happens to sit. At Mercury’s J2000 perihelion longitude (77.457°) and the IAU 2006 obliquity the slope is 1.08036, so the 531.44″ ecliptic advance reads 574.14″/century in right ascension — an excess of 42.71″/century over the ecliptic rate. The general-relativistic advance, derived from the same constants the model uses (GM☉, c, Mercury’s semi-major axis and eccentricity: 6π GM/(c²a(1−e²)) per orbit), is 42.98″/century.

Nothing here is fitted. Three inputs — the 8H/11 divisor, the IAU perihelion longitude, the IAU 2006 obliquity — and one coordinate identity. The agreement is 0.6 %, not exact: ranging determinations pin the inertial excess at 42.980 ± 0.002″/century (Pireaux & Rozelot 2003; Pitjeva’s EPM2008 residual to the relativistic rate is −0.004 ± 0.005″/century), so the projected 42.71″ falls 0.27″/century short of the measured value — many times its uncertainty.

What the model actually measures

The Earth-frame rate the simulation exports for Mercury is ~579.83″/century, flat over the last thousand years and oscillating around the ecliptic value over a full Earth Fundamental Cycle. It decomposes exactly (verified by a gate for all seven planets at 1900, 2000 and 2100, to better than 1″/century):

So the “+48″” the simulation shows above the lattice rate is two coordinate effects: the projection of the advance (+42.71) and the slow decrease of the obliquity (+4.31). The first is the number that coincides with the relativistic anomaly.

The same projection for every planet

The statement is not made for Mercury alone. The identical projection, with each planet’s own lattice rate and J2000 perihelion longitude:

The projection reproduces the relativistic advance for Mercury and for no other planet: where the other planets’ GR advances are small, their projection excesses are large and of either sign. This table is published with the claim as its own test.

Two interpretations compared

Whether the projection is the cause of the anomaly is a question the model states rather than settles: the modern 42.98″ is obtained from ranging fits in an inertial frame with GR inside the dynamical model, and the classical values from transit timings — chains in which no equatorial step is evident. What is established is the identity itself: the model’s ecliptic advance, projected into the equatorial frame, equals the anomaly to within 1 %.

What changes over time

The projection excess is not a constant. It depends on where the perihelion sits relative to the equinox (λ) and on the obliquity (ε), both of which move: the excess is 42.71″ at J2000 and drifts by about +0.3″/century per millennium. General Relativity’s advance is constant. Over a century of ranging the difference is a few hundredths of an arcsecond, at the edge of what the missions resolve; over the full 335,317-year cycle the Earth-frame rate ranges from -47″ to +48″/century around the 531.4″ ecliptic rate.

The model’s Earth-frame RA rate (coordinate b) at 1800, 1900, 2000 and 2100 is 579.84, 579.84, 579.83 and 579.82″/century: essentially flat over the historical record. The model therefore makes no prediction of a measurable decline of the anomaly within the era of precise measurement.

The BepiColombo Test

ESA’s BepiColombo mission arrives at Mercury on 21 November 2026 (delayed from December 2025 due to thruster issues), with orbital commissioning completing around March 2027 and routine science operations starting April 2027. The Mercury Orbiter Radio science Experiment (MORE) will measure Mercury’s orbit with 1–2 orders of magnitude better precision than MESSENGER. This provides the first opportunity to compare two high-precision measurement epochs — MESSENGER (~2013) and BepiColombo (~2027) — separated by ~14 years.

What the ranging missions measure

MESSENGER reported 575.31 ± 0.0015″/century in ICRF coordinates, and BepiColombo will report in the same frame. Both are ranging determinations: Earth–spacecraft distances fitted with a dynamical model that has General Relativity inside it (the PPN parameter β fitted jointly, β ≈ 1). No angle is measured in an Earth-based frame anywhere in that chain, so an equatorial projection cannot appear in it. Under General Relativity BepiColombo returns ~575.31″/century again; under the projection statement the ranging chain is blind to the projection and returns the same number. This test does not discriminate between GR and the projection statement, and the model does not claim that it does. What would discriminate them is the slow drift of the projection excess with λ and ε (about +0.3″/century per millennium, against a constant relativistic advance) — a few hundredths of an arcsecond per century between the two missions, below their stated uncertainties.

The model’s lattice rate for Mercury,

H / (1 + 3/8) = 335,317 / 1.375 = ~243,867 years → ~531.4″/century

is the ecliptic advance; it is 42.98″ short of the ranging value, which contains the relativistic term. Three caveats on the ranging value remain:

1. Epoch. 575.31″/century is a determination at the MESSENGER epoch; BepiColombo gives a second epoch ~14 years later, which bounds any drift at the level of the missions’ uncertainties.
2. The classical chain is different. The historical 43″ (Le Verrier, Newcomb) came mainly from transits of Mercury — timing events in heliocentric ecliptic geometry — with meridian observations reduced from equatorial to ecliptic coordinates on positions, not rates. Whether an equatorial projection can enter any of these chains is the open question the projection reading has to answer; the model states the identity, not the mechanism.
3. GR-inclusive fit caveat. The reported ”575.31″/cy” comes from fitting a GR-inclusive ephemeris to spacecraft ranging data (Park et al. 2017  fit the PPN parameter β jointly, finding β ≈ 1) — conceptually equivalent to measuring the Newtonian baseline and adding the assumed GR contribution. If BepiColombo’s analysis pipeline applies the same GR-inclusive fit, any change in the underlying perihelion advance from frame effects may be absorbed into a slightly different best-fit β, into residuals, or into the orbital baseline — rather than showing up cleanly as a drift in the reported total.

What BepiColombo will and will not test

Values in ICRF as reported by the missions. This test does not probe the projection. For the full scientific discussion including measurement uncertainties and academic critiques, see Scientific Background §4.

Solar Oblateness Uncertainty

The standard Mercury GR test has a rarely-discussed systematic uncertainty: the Sun’s gravitational quadrupole moment (J₂), caused by its oblateness, is not constant — it varies with the solar magnetic activity cycle (~11 years), and published J₂ values have ranged from ~10⁻⁵ to ~10⁻⁷ depending on the method. The solar oblateness contribution has the same temporal signature as the relativistic precession, making them difficult to separate. A 2022 study (MDPI Remote Sensing 14:4139 ) found that an unaccounted-for periodic J₂ component exceeding 0.04% of J₂ could falsely confirm or contradict GR in BepiColombo’s measurements. BepiColombo will improve J₂ determination by 1–2 orders of magnitude, but the time-variable component remains a systematic uncertainty. Detail: Supporting Evidence §5.

Perihelion Precession Across the Solar System

The model calculates perihelion precession for all planets. Each planet has a perihelion point (location of closest approach to the Sun) that slowly drifts:

† Earth’s range comes from its own Earth Rate Deviation (ERD, the deviation of Earth’s perihelion rate from its mean), not the unified 7-planet fluctuation formula. ERD is the underlying cause of the apparent fluctuations for the other planets.

The Mean column is the long-term average over each planet’s full perihelion cycle. The At J2000 column is the model’s epoch-specific rate. The Fluctuation Range is the deviation from the mean over the full Earth Fundamental Cycle. Prograde means counter-clockwise from above the North Pole. Saturn’s perihelion precesses ecliptic-retrograde (clockwise in the ecliptic frame) — see Supporting Evidence §12.

Venus’s fluctuation range (~-29 to +41″/cy) is ~7× larger than Mercury’s despite Venus being much closer to circular. This is what the model predicts: Venus’s poorly-defined perihelion (eccentricity ~0.00678, vs Mercury’s ~0.20564) primarily reflects variations in Earth’s own perihelion rate (ERD), not Venus’s own orbital geometry. See Scientific Background §4 (Q6) for the full discussion.

The Mean and J2000 columns differ because Earth’s reference frame is moving — the same effect detailed above, applied to every planet.

In the Interactive 3D Simulation: open the Show / Hide folder, enable each planet’s perihelion object (e.g., “PERIHELION Mercury”), set “1 second equals” to “1000 years”, press Run.

Key Takeaways

Continue to Mathematical Foundation to see the formal framework behind the model.