This document describes the limitations of the current ascending node calculations in the Holistic Universe Model. The model calculates ascending node precession from a geocentric perspective, which produces values that differ significantly from JPL's published rates.
Scope note (ESSRT). The geometric formula
dΩ/dε = −sin(Ω)/tan(i)and the rate comparison framework are scale-invariant — they apply at any epoch. JPL's published rates use a 1800-2050 fitting interval and the J2000-fixed ecliptic frame; the model's rates use the H-year obliquity cycle as driver. Under ESSRT, H(t) evolves at deep time via Drivers 1 (LOD growth) and 2 (Kepler) — sub-ppm over the modern-era comparison window used here, but conceptually relevant for understanding which side of the geocentric-vs-heliocentric distinction inherits epoch dependence.
The ascending node calculation in calculateDynamicAscendingNodeFromTilts() models how the ascending node appears to shift when observed from Earth's tilting reference frame.
Core formula:
dΩ/dε = -sin(Ω) / tan(i)
Where:
Ω= ascending node longitudeε= Earth's obliquity (driver of the calculation)i= planet's orbital inclination
Direction determination:
- Compares planet's ecliptic inclination against Earth's invariable plane inclination
- If Earth incl > planet incl: ascending node increases when obliquity decreases
- If Earth incl < planet incl: ascending node decreases when obliquity decreases
This calculation captures the geometric transformation effect: when Earth's reference frame (the ecliptic) tilts over time, the line where other planetary orbits cross this plane shifts accordingly. This is a coordinate transformation effect, not a physical movement of the orbital planes.
JPL Approximate Positions of the Planets
JPL's rates represent actual gravitational precession of each planet's orbital plane in inertial space:
- Reference frame: Fixed J2000 mean ecliptic and equinox
- Method: Best-fit Keplerian elements to DE200 ephemeris
- Time span: 1800-2050 AD (250 years)
- Cause: Gravitational perturbations from other planets
| Planet | Rate (°/century) | Rate (arcsec/century) |
|---|---|---|
| Mercury | -0.12534081 | -451 |
| Venus | -0.27769418 | -1000 |
| Earth | 0.0 | 0 |
| Mars | -0.29257343 | -1053 |
| Jupiter | +0.20469106 | +737 |
| Saturn | -0.28867794 | -1039 |
| Uranus | +0.04240589 | +153 |
| Neptune | -0.00508664 | -18 |
| Planet | Our Model (arcsec/cy) | JPL (arcsec/cy) | Ratio | Notes |
|---|---|---|---|---|
| Mercury | -285 | -451 | 0.63 | Same direction |
| Venus | -770 | -1000 | 0.77 | Same direction |
| Mars | -1100 | -1053 | 1.04 | Close match |
| Jupiter | +2000 | +737 | 2.7 | Same direction, large difference |
| Saturn | -1000 | -1039 | 0.96 | Close match |
| Uranus | +3300 | +153 | 21.6 | Large discrepancy |
| Neptune | -1130 | -18 | 62.8 | Large discrepancy |
- Inner planets and Saturn: Reasonable agreement (within ~40%)
- Jupiter: Same direction but 2.7× larger in model
- Uranus and Neptune: Extremely large discrepancies (21× and 63×)
| Aspect | Our Model | JPL Rates |
|---|---|---|
| Physical cause | Reference frame rotation | Gravitational torques |
| What moves | The coordinate system | The orbital plane itself |
| Reference frame | Moving (current) ecliptic | Fixed J2000 ecliptic |
| Timescale | H-year obliquity cycle | 10^4 to 10^7 year precession periods |
Conclusion: The model and JPL measure fundamentally different effects. They are not directly comparable and should not be expected to match.
The formula dΩ/dε = -sin(Ω) / tan(i) has tan(i) in the denominator:
| Planet | Inclination | tan(i) | Amplification Factor |
|---|---|---|---|
| Mercury | 7.0° | 0.123 | 8× |
| Saturn | 2.5° | 0.044 | 23× |
| Mars | 1.85° | 0.032 | 31× |
| Neptune | 1.77° | 0.031 | 32× |
| Jupiter | 1.31° | 0.023 | 44× |
| Uranus | 0.77° | 0.013 | 77× |
Uranus and Neptune have very low inclinations, causing extreme amplification. This may indicate the formula is not appropriate for near-coplanar orbits.
JPL's ascending node rates have inherent limitations:
- Short fitting interval: 250 years vs. precession periods of 10,000+ years
- Linear approximation: True nodal motion is a superposition of multiple oscillating eigenmodes
- Purpose-driven accuracy: Optimized for position accuracy, not rate accuracy
- Small rates are noisy: Neptune's -18 arcsec/cy rate is difficult to measure precisely over 250 years
From JPL's documentation:
"The elements are not intended to represent any sort of mean; they are simply the result of being adjusted for a best fit."
For perihelion precession, the model includes a translation from geocentric to heliocentric rates via the "extra fluctuation" factor. A similar translation may be needed for ascending nodes but has not been implemented.
Perihelion example (Mercury at J2000):
- Heliocentric baseline rate: 531.4 arcsec/cy
- Geocentric rate (predictive formula): 569.2 arcsec/cy
- Extra fluctuation: +37.7 arcsec/cy
If applied to ascending nodes:
- Geocentric rate: -285 arcsec/cy
- Plus extra fluctuation: +37.7 arcsec/cy
- Result: -247 arcsec/cy
This still does not match JPL's -451 arcsec/cy, suggesting the translation is more complex for ascending nodes than for perihelion.
| Effect | Cause | Rate Example | Our Model? |
|---|---|---|---|
| Gravitational precession | Planet-planet perturbations | Mercury: -451 arcsec/cy | No |
| Ecliptic plane precession | Perturbations on Earth's orbit | ~47 arcsec/cy | No |
| Reference frame tilt | Earth's obliquity/inclination cycles | Mercury: -285 arcsec/cy | Yes |
JPL's rates can theoretically be derived from Lagrange-Laplace secular perturbation theory, which involves:
- B Matrix: Encodes gravitational interactions between all planets
- Eigenfrequencies (s₁, s₂, ... s₈): Fundamental nodal precession frequencies
- Superposition: Each planet's nodal motion is a sum of eigenmodes
However, there is no simple closed-form formula to directly calculate rates like Mercury's -451 arcsec/cy. The full calculation requires:
- Planetary masses
- Semi-major axes
- Laplace coefficients (special functions)
- Matrix eigenvalue decomposition
The current implementation correctly models the geocentric reference frame effect. For visualizing how ascending nodes appear to shift from Earth's perspective over long timescales, the model is valid.
Two options:
Option A: Use JPL rates directly
const ascendingNodeRates = {
mercury: -451, // arcsec/century
venus: -1000,
mars: -1053,
jupiter: +737,
saturn: -1039,
uranus: +153,
neptune: -18
};Option B: Derive from perihelion precession pattern
This would require understanding how the geocentric-to-heliocentric translation differs between perihelion and ascending node precession. Current investigation suggests this is not a simple additive relationship.
The current model should be understood as calculating a geocentric coordinate transformation effect, not matching JPL's heliocentric gravitational precession rates. Both are physically meaningful but measure different phenomena.
- JPL Approximate Positions of the Planets
- Secular Evolution of Planetary Orbits (Laplace-Lagrange theory)
- La2010: Orbital Solution for Earth (A&A)
- Dynamic Ascending Node Calculation - Main implementation documentation
- Dynamic Ascending Node Calculation - Implementation details
- Dynamic Elements Overview - Master overview of all dynamic systems
- Constants Reference - All constants and source values
- Expanding Solar System Resonance Theory - Deep-time scaling of H(t)