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, plus, for Mercury above all, the relativistic advance — and so does the model: Mercury’s stated total is the lattice rate plus a derived relativistic supplement (see below). What the model adds is coordinate bookkeeping done exactly: the apparent rates an Earth observer’s equatorial frame produces are exact projections of the ecliptic advances (a gate verifies the decomposition for all seven planets), which matters for reading the simulation’s Earth-frame exports.

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 against the moving equinox 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 classical column is measured against the moving vernal equinox, which the equinox precession (~5,026″/century in the Newcomb constant system Clemence used) carries backward; the ICRF column is measured against the fixed stars. The equinox precession cancels in the subtraction, which is why the ~43″ discrepancy comes out the same in both frames. Each column is a complete determination inside its own constant system and epoch — the two totals are not inter-convertible by adding the modern precession rate to the ICRF value; no equinox-referred determination has been made in the modern system, and the subtraction must stay inside one column.

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 Settled Position — and a Coincidence Worth the Record

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.

The verdict

The question is settled, against the projection. Three independent tests close it: (1) Precision — ranging pins the inertial excess at 42.980 ± 0.002″/century; the projection gives 42.71 — a 0.27″ shortfall, over a hundred times the measurement uncertainty. (2) The all-planet test — the identical projection fails for every other planet (table above), while ranging confirms the other planets’ relativistic advances directly (Venus 8.62, Earth 3.84, Mars 1.35″/cy). (3) No mechanism — the classical chains form the perihelion correction in ecliptic longitude (Le Verrier’s Paris meridian series alone, re-solved as pure geocentric longitudes with no transit and no modern constant in the chain, yields the excess: δπ′ = 43.1 ± 16.5″/cy), and the ranging chains form no equatorial angle anywhere; there is no step on which the projection slope could act. Decisively, the model’s own N-body engine requires the relativistic term to match the observed window rate. What survives is the identity itself, kept as the record: a striking 0.6 % numerical coincidence at Mercury’s particular perihelion longitude — published together with the table that refutes it. The two effects were never the same kind of thing: the relativistic advance is a physical effect in the ecliptic, the same in every frame; the projection excess is a property of the measuring coordinate, drifting as λ and ε move (~+0.3″/cy per millennium) and swinging between −44″ and +48″ over Mercury’s full perihelion cycle.

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 makes no prediction of a measurable decline of the anomaly within the era of precise measurement — its forward prediction is General Relativity’s.

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). The model’s prediction for BepiColombo is the same as General Relativity’s — ~575.31″/century again — because the model carries the relativistic advance as its derived supplement. What the mission adds for the model is the same thing it adds for GR: a second high-precision epoch bounding any secular change of the advance, and a much sharper solar-J₂ separation (below).

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 agrees. The historical 43″ (Le Verrier, Newcomb) came mainly from transits of Mercury — timing events in heliocentric ecliptic geometry — with meridian observations reduced to ecliptic coordinates on positions, not rates. No equatorial angle is formed in any of these chains — one of the three findings that closed the projection reading (above).
3. What the fit actually determines. 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) — the advance is constrained through the joint fit rather than read off as an angle. The precision of that constraint (±0.0015″/cy, with the EPM residual to the relativistic rate at −0.004 ± 0.005″/cy) is exactly what made the projection candidate’s 0.27″ shortfall a decisive refutation.

What BepiColombo tests

Values in ICRF as reported by the missions. The model and General Relativity make the same prediction here — the model ships the relativistic supplement — so BepiColombo is a consistency check for both, and its real gain for this page’s subject is the solar-J₂ separation below. For the historical record of the projection candidate and its closure, 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 for the interpretation.

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.