Perihelion Precession
Every planet’s perihelion (closest approach to the Sun) slowly rotates around the Sun — a phenomenon called perihelion precession. This happens because of gravitational perturbations from other planets.
The rate of precession can be estimated using Lagrange-Laplace secular perturbation theory — an 18th-19th century mathematical framework that treats gravitational interactions as averaged effects over complete orbits. The table below compares observed rates, L-L theoretical predictions, and the Holistic Universe Model’s predictions at epoch J2000 (all relative to fixed stars / ICRF):
| Planet | Fitzpatrick (″/cy) | L-L Theory (″/cy) | Model at J2000 (″/cy) |
|---|---|---|---|
| Mercury | ~575 | ~554 | ~573 |
| Venus | ~205 | ~1,207 | ~427 |
| Earth | ~1,145 | ~1,279 | ~1,164 |
| Mars | ~1,628 | ~1,775 | ~1,567 |
| Jupiter | ~655 | ~751 | ~1,797 |
| Saturn | ~1,950 | ~1,859 | ~-3,385 |
| Uranus | ~334 | ~275 | ~1,073 |
| Neptune | ~36 | ~67 | ~204 |
Why do these values differ so much? The “Fitzpatrick” column (from Table 5.1 ) represents a single-epoch snapshot at J2000. Perihelion precession rates fluctuate significantly over time — these fluctuations are clearly visible in WebGeocalc data. A single measurement is not representative of the long-term direction of movement. The model’s J2000 values are based on the trend from 1900 to 2000 from WebGeocalc, which better captures the actual direction of perihelion movement. For Jupiter and Saturn in particular, the model predicts their current perihelion trends will simply continue as-is, while WebGeocalc predicts a pattern change — see Predictions: Jupiter and Saturn Perihelion Trends. The model’s long-term mean rates are shown below.
Source: Fitzpatrick and L-L columns from Fitzpatrick, R. An Introduction to Celestial Mechanics, University of Texas at Austin (Table 5.1 ). L-L theory assumes nearly circular, coplanar orbits. Venus shows the poorest agreement because its unusually low eccentricity (~0.007) makes its perihelion point extremely sensitive to small perturbations. More precise values (like Mercury’s ~532″) come from detailed N-body calculations rather than L-L theory alone.
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
| Measurement | Relative to fixed stars (ICRF) | Relative to moving equinox |
|---|---|---|
| Observed precession | ~575″/century | ~5,604″/century (~575 + ~5,028.8) |
| Newtonian prediction | ~532″/century | ~5,561″/century (~532 + ~5,028.8) |
| Discrepancy | ~43″/century | ~43″/century |
Both columns are measured in the ecliptic plane — the difference is the reference direction. The left column measures against fixed stars (ICRF) — the inertial frame defined by distant quasars. The right column measures against the moving vernal equinox, which drifts backward at ~5,028.8″/century due to Earth’s axial precession. 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 the older ~5,025″ equinox precession value).
Where Does the ~532″ Come From?
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):
| Planet | Contribution | Percentage |
|---|---|---|
| Venus | ~278 arcsec/century | ~52% |
| Jupiter | ~154 arcsec/century | ~29% |
| Earth | ~90 arcsec/century | ~17% |
| Saturn | ~7 arcsec/century | ~1% |
| Mars, Uranus, Neptune | ~3 arcsec/century | less than 1% |
| Total (Newtonian) | ~532 arcsec/century | 100% |
These contributions are calculated using detailed N-body perturbation calculations — more precise than the first-order L-L estimates shown above, which overestimate Mercury’s total Newtonian rate at ~554″ due to their simplifying assumptions (circular, coplanar orbits).
Data Sources
The perihelion precession data used in the model comes from WebGeocalc , NASA’s tool for calculating orbital elements.
The graph shows Mercury’s perihelion movement over 6 centuries confirming the observed ~575″/century rate (ICRF). Subtracting the Newtonian prediction (~532″) from this observed rate gives the ~43″/century discrepancy — or equivalently, in equinox-based terms: ~5,604″ (observed) minus ~5,561″ (Newtonian) = ~43″.
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 exactly ~43 arcseconds per century
This was one of the first major confirmations of Einstein’s theory.
The Model’s Alternative Explanation
Alternative Proposal: The Holistic Universe Model proposes that the ~43 arcsecond discrepancy may not be caused by relativistic effects, but by Earth’s reference frame motion. This is a testable alternative interpretation, not a claim that General Relativity is wrong. See Scientific Background for detailed discussion including academic critiques and measurement uncertainties.
Why the Reference Frame Matters
When we observe Mercury from Earth, we’re not observing from a fixed point. Earth itself is moving in multiple ways:
- Earth wobbles around the EARTH-WOBBLE-CENTER over ~25,684 years (axial precession)
- PERIHELION-OF-EARTH orbits the Sun in the opposite direction over ~111,296 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 what is observed — they disagree on what causes the “extra” precession beyond Newtonian mechanics. The geocentric value (~5,600″) is what was historically measured on Earth; the ICRF value (~575″) and the “anomaly” (~43″) are derived from it by subtracting the equinox drift (~5,028.8″):
| Standard (GR) | Model (epoch 2000) | |
|---|---|---|
| Newtonian perturbations | ~532″/cy | ~533.7″/cy (model baseline) |
| Additional advance | +~43″ (space-time curvature, constant) | +~39″ (Earth-frame offset, variable) |
| = Observed (ICRF) | ~575″/cy | ~573″/cy |
| + Equinox drift | +~5,028.8″ | +~5,028.8″ |
| = Observed (geocentric) | ~5,604″/cy (constant) | ~5,601″/cy (decreasing) |
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″ around 1900 (matching Einstein’s era) but has decreased to ~39″ by 2000. The model predicts the geocentric value will continue to decrease over time.
The diagram below shows this in detail — Mercury’s perihelion positions at epoch 1900 and 2000, with the model’s baseline (~533.7″/century) and the observed rate (~572.5″/century ICRF / ~5,601.3″/century geocentric at epoch 2000):
A Key Prediction
If this model is correct, the observed total precession should change over time as Earth’s precession cycles progress:
| Year | Relative to moving equinox (geocentric) | Relative to fixed stars / ICRF (heliocentric) | “Anomaly” (observed − ~533.7″ model baseline, cf. ~532″ Newtonian) |
|---|---|---|---|
| 1800 | ~5,609.20″/century | ~580.35″/century | ~46.64″ |
| 1900 | ~5,605.44″/century | ~576.62″/century | ~42.91″ |
| 2000 | ~5,601.34″/century | ~572.54″/century | ~38.82″ |
| 2100 | ~5,596.89″/century | ~568.09″/century | ~34.38″ |
Key point: The equinox-based geocentric values (~5,600″) are what was historically measured on Earth — they include the drift of the vernal equinox (~5,028.8″/century). The ICRF heliocentric values (~575″) are obtained by subtracting this equinox drift. Both decrease by ~4–4.5″/century. The “anomaly” decreases accordingly. In Einstein’s era (~1900) it was ~43″ — matching the famous value — but the model predicts it will continue to decrease.
The standard General Relativity explanation predicts the geocentric value remains constant at ~575 + 5,028.8 = ~5,604″/century. The model predicts this value will decrease — from ~5,601″ (2000) toward ~5,597″ (2100) — as Earth’s precession cycles progress.
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,420 years (1/45 of the Holistic-Year). This oscillation arises from the beat frequency between Earth’s inclination precession (111,296 years) and ecliptic precession (66,778 years):
1/66,778 − 1/111,296 = 1/166,890 ≈ 1/7,420 × 1/22.5Over the full 333,888-year Holistic cycle, the fluctuation ranges from −157″ to +174″/century around the ~534″ baseline. The current era (~+39″) happens to be near the historical ~43″ value because we are on the descending slope of one oscillation. The two independent measurement methods — direct 3D simulation output and analytical formula (R² = 0.999) — produce matching results. See Formula Derivation: Mercury Key Combination Periods for the complete frequency mixing analysis.
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.
Why the Model Does Not Match MESSENGER’s 575.31″
The model sets Mercury’s base perihelion period using the Fibonacci-derived fraction:
Holistic Year / (1 + 3/8) = 333,888 / 1.375 = ~242,828 years → ~533.7″/century base rateAt epoch ~2000, the Earth-frame offset adds ~38.8″, giving an ICRF rate of ~572.54″/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 because 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.
For context, the full geocentric picture (including Earth’s axial precession of ~5,028.8″/century) is:
Model at epoch ~2000: 533.7 + 38.8 + 5,028.8 = ~5,601″/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,601″ (epoch 2000) toward ~5,597″ (epoch 2100) — because the Earth-frame offset (currently ~39″) is shrinking as Earth’s reference frame moves (axial precession, perihelion counter-movement). GR predicts it remains constant at ~5,604″/century.
Two Possible Outcomes
Scenario A — Rate has decreased (supports the model): If BepiColombo measures ~574.69″/century or lower (ICRF) versus MESSENGER’s 575.31″/century, this would be evidence that the precession rate is changing over time — exactly what the model predicts (~0.045″/century per year of decrease). 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.
Scenario B — Rate is unchanged (supports GR): If BepiColombo measures approximately 575.31″/century (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
The Numbers
| MESSENGER (~2013) | BepiColombo (~2027) | Difference | |
|---|---|---|---|
| Model predicts | 575.31″/cy (measured) | ~574.69″/cy or lower | −0.62″/cy or more |
| GR predicts | 575.31″/cy (measured) | ~575.31″/cy | ~0 (constant) |
| Measurement precision | ±0.0015″/cy |
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,603″ (model) or ~5,604″ (GR) for BepiColombo.
The predicted difference of 0.62″/century is ~400× larger than MESSENGER’s measurement uncertainty — making this a decisive test.
The model does not claim General Relativity is wrong as a theory - only that this particular phenomenon (Mercury’s perihelion anomaly) may have an alternative explanation based on reference frame motion. As Křížek and Somer note in their book Mathematical Aspects of Paradoxes in Cosmology (2023), “the definitive proof of these famous theories is still to be delivered.”
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:
| Issue | Detail |
|---|---|
| J₂ is not constant | It varies with the solar magnetic activity cycle (~11 years) |
| Measurements disagree | Published J₂ values have ranged from 1.08 × 10⁻⁵ to 1.46 × 10⁻⁷ depending on the method |
| J₂ mimics GR | 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(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.
Visualizing Planetary Perihelions
The Holistic Universe Model calculates perihelion precession for all planets. Each planet has a PERIHELION-POINT that orbits the Sun, creating the precession we observe:
| Planet | PERIHELION-POINT Orbital Period | Direction | Mean (″/cy) | At J2000 (″/cy) |
|---|---|---|---|---|
| Mercury | ~242,828 years | Prograde | ~533.7 | ~573 |
| Venus | ~667,776 years | Prograde | ~194 | ~427 |
| Earth | ~111,296 years | Prograde | ~1,164 | ~1,160 |
| Mars | ~77,051 years | Prograde | ~1,682 | ~1,567 |
| Jupiter | ~66,778 years | Prograde | ~1,941 | ~1,797 |
| Saturn | ~41,736 years | Ecliptic-retrograde | ~-3,106 | ~-3,385 |
| Uranus | ~111,296 years | Prograde | ~1,164 | ~1,073 |
| Neptune | ~667,776 years | Prograde | ~194 | ~204 |
Note: The Mean column shows the long-term average rate over each planet’s full PERIHELION-POINT cycle — derived directly from the orbital period (1,296,000″ / period in centuries). The At J2000 column shows the model’s epoch-specific rate based on the 1900–2000 trend from WebGeocalc . The difference between Mean and J2000 reflects fluctuations caused by Earth’s reference frame motion and gravitational perturbations. Prograde means the same direction as orbital motion (counter-clockwise when viewed from above the North Pole). Saturn’s perihelion precession is ecliptic-retrograde (clockwise in the ecliptic frame), moving opposite to all other planets as seen from Earth.
The mean values are the model’s baseline precession rates — the corrected long-term averages without fluctuations. The values observed from Earth differ because Earth’s own reference frame is moving (see The Model’s Alternative Explanation above). For Mercury, the mean baseline is ~534″/century, but the J2000 rate is ~573″/century — a difference of ~39″ that the model attributes to Earth’s wobble and perihelion counter-movement. This offset varies over time and affects all planets, not just Mercury.
In the Interactive 3D Simulation, you can see all planetary perihelion points:
- Open the Show / Hide folder and enable the perihelion objects for each planet (e.g., “PERIHELION Mercury”)
- Set “1 second equals” to “1000 years”
- Press “Run” to see the perihelion points precess around the Sun
The perihelion points form a spiral pattern that evolves over one Holistic-Year (333,888 years).
Predictive Formulas for All Planets: The model includes predictive formulas for all 8 planets (Mercury through Neptune) that require only a year as input—no observations needed. These achieve R² > 0.99 accuracy for all planets. See Formulas: Predictive Formulas for complete implementations.
Key Takeaways
| Question | Answer |
|---|---|
| What is perihelion precession? | The slow rotation of a planet’s closest approach point around the Sun |
| What’s Mercury’s cycle? | ~242,828 years (prograde) |
| What’s the “anomaly”? | ~43 arcsec/century difference between observed (~575″) and Newtonian (~532″) |
| Standard explanation? | General Relativity’s space-time curvature (~43 arcsec/century) |
| Model’s alternative? | Earth’s reference frame motion may explain the anomaly |
| Testable prediction? | Model: geocentric decreases ~5,601→~5,597″; GR: constant ~5,604″ |
| Full discussion? | Scientific Background |
- All planets have perihelion precession - Mercury has a ~242,828-year prograde cycle
- The ~43 arcsecond “anomaly” has been attributed to General Relativity since 1915
- The model proposes an alternative - Earth’s reference frame motion may create the apparent anomaly
- Falsifiable prediction - if geocentric precession decreases from ~5,601″ toward ~5,597″/century (~572.54→~568.09″ heliocentric), it supports the model; if it stays constant at ~5,604″, it refutes the alternative
- Reference direction matters — ~575″/century relative to fixed stars (ICRF) vs ~5,604″/century relative to the moving vernal equinox (~575 + ~5,028.8 equinox drift). Both are in the ecliptic plane; the difference is whether you measure against fixed stars or Earth’s moving equinox
Continue to Mathematical Foundation to see the formal mathematical framework behind the model.