Perihelion Precession

Every planet’s perihelion (closest approach to the Sun) slowly rotates around the Sun — a phenomenon called 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.

The rate of precession can be estimated by averaging gravitational pulls from other planets over many orbits. The table below compares the WebGeocalc observed rate (1800–2000 trend), the Holistic Universe Model’s prediction, and two theoretical predictions — Fitzpatrick’s textbook secular value (Standish & Williams 1992) and a 19th-century analytical approximation (Lagrange–Laplace, “L-L Theory”). All values at epoch J2000, relative to fixed stars (ICRF — the inertial reference frame defined by distant quasars):

WebGeocalc = JPL ephemeris 1800–2000 trend. Model = Holistic Universe Model J2000 prediction (fit to WebGeocalc). Fitzpatrick = long-term secular average (Standish & Williams 1992). L-L Theory = analytical first-order approximation (Lagrange–Laplace, 19th century).

Heads-up before continuing: Mercury’s WebGeocalc trend above (~572″/cy) is ~3″/cy below the canonical ~575″/cy used in the textbook anomaly story below. That gap matters — the standard “anomaly = 575 − 532 = 43″/cy” equation needs the higher figure to recover General Relativity’s predicted contribution. 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 (see The BepiColombo Test below).

The model’s data source is WebGeocalc , NASA’s tool for calculating orbital elements. The graph below shows Mercury’s perihelion movement over 6 centuries:

The Mercury “Anomaly”

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

The Numbers

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) is derived by subtracting that drift, so it’s measured against fixed stars instead. The equinox drift cancels out when you subtract, so the ~43″ discrepancy is the same in both frames.

The equinox-based values (~5,600″) are what was historically measured on Earth before ICRF corrections existed. The commonly cited “~5,600″” is a rounded approximation from Clemence (1947), who used ~5,025″ for the equinox precession rate available at the time. Modern measurements give ~5,028.8″/cy — the figure 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 using N-body integration of Newton’s inverse-square law over time — typically using JPL’s DE-series planetary ephemeris, not a closed-form equation. The first specific Mercury calculation was by Le Verrier (1859, ~527″/century), who 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 L-L Theory column’s higher ~554″/century is the analytical first-order version of the same Newtonian calculation (see the callout earlier on this page for details).

The Standard Explanation

The standard explanation for this 43 arcsecond discrepancy is Einstein’s General Relativity (1915):

• Space-time is curved near massive objects
• Mercury, being closest to the Sun, experiences the strongest curvature
• This causes an additional precession of ~43 arcseconds per 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 Reference-Frame Effect

When we observe Mercury from Earth, we’re not observing from a fixed point. Earth itself is moving in multiple ways:

1. Axial precession — Earth’s spin axis wobbles over ~25,794 years, tracing a small circle around its mean orientation
2. Inclination precession — the direction toward Earth’s perihelion (its closest approach to the Sun) drifts the opposite way over ~111,772 years

These two motions change Earth’s orientation axis over time. When we measure Mercury’s perihelion precession from this moving reference frame, we get a different value than the “true” heliocentric rate.

Two Interpretations Compared

The standard view and the model agree on the observed numbers; they disagree on the cause and stability of the anomaly:

The key difference: GR’s +~43″ is a fixed physical effect, so the geocentric value stays constant at ~5,604″. The model’s Earth-frame offset is variable — it was ~43.01″ around 1900, essentially identical to Newcomb’s 1882 measurement (the value Einstein later derived from GR), but has decreased to ~38″ by 2000. Under the model, the historical 1900 match is exactly what time-varying frame effects predict for that specific epoch; under GR it should persist permanently. The model predicts the geocentric value will continue to decrease over time.

A Key Prediction

If this model is correct, the observed total precession should change over time as Earth’s precession cycles progress:

Key point: Both columns decrease by ~5.0″/century, so the “anomaly” shrinks too. In Einstein’s era (~1900) it was ~43.01″ — matching the famous value — but the model predicts it will keep decreasing.

The Oscillation Structure

The decline shown above is not smooth — the model predicts Mercury’s Earth-frame precession rate oscillates with a dominant period of ~7,451 years (1/45 of the Earth Fundamental Cycle). Over the full 335,317-year Earth Fundamental Cycle, the fluctuation ranges from -180″ to +202″/century around the ~531.4″ baseline. The current era (~+38.02″) happens to be near the historical ~43.01″ value because we are on the descending slope of one oscillation.

The BepiColombo Test

ESA’s BepiColombo mission arrives at Mercury on 21 November 2026 (delayed from the original December 2025 date 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.

Mismatch with MESSENGER’s 575.31″

The model sets Mercury’s base perihelion period using the Fibonacci-derived fraction:

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

At epoch ~2000, the Earth-frame offset adds ~38.0″, giving an ICRF rate of ~569.46″/century — not 575.31″. MESSENGER reported 575.31 ± 0.0015″/century in ICRF (heliocentric) coordinates, and BepiColombo will report in the same frame.

The model does not force a match to MESSENGER’s value for three reasons. First, we do not yet know whether 575.31″/century is a stable long-term rate or an epoch-specific measurement — BepiColombo will answer this question. Second, the ICRF value of 575.31″ is not a direct measurement: it is derived from the geocentric measurement (~5,604″/cy) by subtracting the equinox drift, then the Newtonian contribution. Whether those subtractions fully account for all reference-frame effects is precisely the question the model raises. Third, 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 (~532″/cy) and adding the assumed GR contribution (~43″/cy). If BepiColombo’s analysis pipeline applies the same GR-inclusive fit, any change in the underlying perihelion advance from frame effects may be absorbed into a slightly different best-fit β, into residuals, or into the orbital baseline — rather than showing up cleanly as a drift in the reported total.

For context, the full geocentric picture (including Earth’s axial precession of ~5,028.8″/century) is:

Model at epoch ~2000:      531.4 + 38.0 + 5,028.8 = ~5,598.26″/century (geocentric)
MESSENGER at epoch ~2013:  575.31 + 5,028.8        = ~5,604″/century (geocentric)

The model predicts this geocentric rate will decrease over time — from ~5,598.26″ (epoch 2000) toward ~5,592.85″ (epoch 2100) — because the Earth-frame offset (currently ~38″) is shrinking as Earth’s reference frame moves (axial precession, the slow drift of Earth’s perihelion direction). GR predicts it remains constant at ~5,604″/century.

Two Possible Outcomes

Scenario A — Geocentric rate has decreased (supports the model): If the geocentric precession (according to MESSENGER ~5,604″/century) has decreased, this would be evidence that the precession rate is changing over time — exactly what the model predicts. The corresponding ICRF value would be ~574.61″/century or lower (vs MESSENGER’s 575.31″). In this case:

• The “missing” ~43 arcseconds would be a variable quantity, not a fixed relativistic effect
• The model’s Earth-frame interpretation gains support
• We could then refine Mercury’s base rate to match the observed trend more precisely

This conclusion holds as long as BepiColombo’s analysis pipeline reports the raw measured perihelion advance rather than a GR-inclusive fit total.

Scenario B — Geocentric rate is unchanged (supports GR): If the geocentric precession remains constant at ~5,604″/century (corresponding to ~575.31″ ICRF), within measurement uncertainty, this would be evidence that the precession rate is constant — consistent with General Relativity’s prediction that the ~43″ is a fixed space-time curvature effect. In this case:

• The model’s alternative explanation for Mercury’s perihelion anomaly would be refuted
• The ~43 arcsecond advance is a permanent physical effect, not a reference frame artifact

Alternatively, this outcome could indicate that any model-predicted drift has been absorbed by the GR-inclusive analysis pipeline (different best-fit β, residuals, or orbital baseline) — see the third reason in §Mismatch with MESSENGER’s 575.31″ above.

MESSENGER vs BepiColombo

Values shown in ICRF (heliocentric) as reported by the missions. The geocentric equivalents (adding ~5,028.8″/cy equinox drift) are ~5,604″ for MESSENGER and ~5,598.26″ (model) or ~5,604″ (GR) for BepiColombo.

The predicted difference of 0.70″/century is ~500× larger than MESSENGER’s measurement uncertainty — large enough to be decisive if BepiColombo’s analysis pipeline reports the raw measured perihelion advance, not a GR-inclusive fit total (see the third reason above).

For the full scientific discussion including measurement uncertainties and academic critiques, see Scientific Background: The Mercury Perihelion Question.

Solar Oblateness Uncertainty

An often-overlooked systematic uncertainty in the standard Mercury GR test involves the Sun’s gravitational quadrupole moment (J₂), caused by its oblateness:

A 2022 study (MDPI Remote Sensing 14(17), 4139 ) found that if a periodic J₂ component exceeding 0.04% of J₂ exists and is not accounted for, it could falsely confirm or contradict General Relativity in BepiColombo’s measurements.

BepiColombo will improve J₂ determination by 1–2 orders of magnitude. However, the time-variable component remains a systematic uncertainty that must be carefully modeled when interpreting the precession results. This does not invalidate the GR test, but it highlights that the standard Mercury test has an unresolved systematic uncertainty rarely discussed in popular presentations.

Perihelion Precession Across the Solar System

The Holistic Universe Model calculates perihelion precession for all planets. Each planet has a perihelion point — the location of its closest approach to the Sun — that itself slowly drifts, creating the precession we observe:

† 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 seen for the other planets.

Notice Venus’s fluctuation range (~-1,353 to +1,231″/cy) is ~7× larger than Mercury’s (~-180 to +202″/cy) — even though Venus is much closer to circular than Mercury. This is exactly 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) rather than Venus’s own orbital geometry. See Scientific Background: Questions for This Interpretation (Q6) for the full discussion.

The Mean column shows long-term baselines; the J2000 column shows the current epoch-specific rate. They differ because Earth’s reference frame is moving — the same effect detailed in The Model’s Alternative Explanation above, applied to every planet, not just Mercury.

In the Interactive 3D Simulation, you can see all planetary perihelion points:

1. Open the Show / Hide folder and enable the perihelion objects for each planet (e.g., “PERIHELION Mercury”)
2. Set “1 second equals” to “1000 years”
3. Press “Run” to see the perihelion points precess around the Sun

Key Takeaways

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