Eccentricity
Eccentricity measures how elliptical Earth’s orbit is. A value of 0 would be a perfect circle; higher values mean a more elongated ellipse. Earth’s current eccentricity is 0.01671022 - a nearly circular orbit.
What Eccentricity Means in Practice
The eccentricity value (0.01671022) represents the offset distance between the center of Earth’s orbit and the Sun, expressed as a fraction of the orbital radius (1 AU).
| Measurement | Value |
|---|---|
| 1 AU (mean Earth-Sun distance) | 149,597,870.698828 km |
| Eccentricity (J2000) | 0.01671022 |
| Offset distance | 2,499,813 km |
| Perihelion distance | ~147.1 million km |
| Aphelion distance | ~152.1 million km |
| Difference | ~5 million km |
This means:
- At perihelion (closest, ~January 3): Earth is ~147.1 million km from the Sun
- At aphelion (farthest, ~July 4): Earth is ~152.1 million km from the Sun
- Earth receives about 7% more solar energy at perihelion than at aphelion
The 20,957-Year Cycle
In the Holistic Universe Model, eccentricity changes in a predictable 20,957-year cycle - not the ~100,000 and ~400,000-year cycles predicted by Milankovitch theory.
The Mechanism
Two motions work in opposite directions:
| Motion | Direction | Period |
|---|---|---|
| Earth around EARTH-WOBBLE-CENTER | Clockwise | ~25,794 years |
| PERIHELION-OF-EARTH around Sun | Counter-clockwise | ~111,772 years |
Because they move in opposite directions, they meet more frequently than either cycle alone:
Meeting frequency = 1/25,794 + 1/111,772 = 1/20,957
Therefore: They meet every 20,957 years
This is the perihelion precession cycle.
Why Alignment Affects Eccentricity
The PERIHELION-OF-EARTH defines where Earth’s closest approach to the Sun occurs. Earth orbits the EARTH-WOBBLE-CENTER at a small radius (~202,881 km).
When Earth and PERIHELION-OF-EARTH are on the same side of EARTH-WOBBLE-CENTER:
- Their distances add together
- Maximum eccentricity (~0.0167)
When Earth and PERIHELION-OF-EARTH are on opposite sides:
- Their distances partially cancel
- Minimum eccentricity (~0.0140)
Eccentricity Values
| Parameter | Value | Notes |
|---|---|---|
| Current eccentricity (J2000) | 0.01671022 | Measured, NASA Planetary Fact Sheet |
| Base eccentricity | 0.015386 | Model-derived arithmetic midpoint (derived mean = 0.0154454). |
| Maximum eccentricity | ~0.0167 | At winter solstice alignment |
| Minimum eccentricity | ~0.0140 | At summer solstice alignment |
| Variation amplitude | ±0.001356 | Half the range |
| Cycle period | 20,957 years | 335,317 ÷ 16 |
How Base Eccentricity Was Derived
The base value (0.015386) — the arithmetic midpoint of the eccentricity cycle — cannot be measured directly because we only have observations from recent centuries. It was derived using three constraints:
- Minimum eccentricity occurred in ~9,233 BC when perihelion aligned with the June solstice
- Maximum eccentricity occurred in 1246 AD when perihelion aligned with the December solstice
- Current eccentricity (0.01671022) is near the maximum and decreasing
The 3D Simulation was calibrated to satisfy all three constraints, yielding:
- Base (arithmetic midpoint) = 0.015386
- Amplitude = ±0.001356
The Solstice Connection
Eccentricity extremes correlate with solstice alignments:
| Alignment | Eccentricity | Last Occurrence | Next Occurrence |
|---|---|---|---|
| Perihelion at December solstice | Maximum (~0.0167) | 1246 AD | ~22,203.344 AD |
| Perihelion at June solstice | Minimum (~0.0140) | ~9,233 BC | ~11,725 AD |
Why this correlation?
When perihelion aligns with the December solstice (Northern Hemisphere winter), Earth and PERIHELION-OF-EARTH are positioned such that their orbital offsets add together. When aligned with the June solstice, they partially cancel.
Current Status
We passed maximum eccentricity around 1246 AD. The current value (0.01671022) is:
- Decreasing toward the mean
- Will reach minimum (~0.0140) around 11,725 AD
- Will return to maximum around 22,203.344 AD
Why Not Milankovitch’s 100k/400k Cycles?
The conventional Milankovitch theory proposes eccentricity cycles of ~100,000 and ~400,000 years. The model proposes a simpler 20,957-year cycle instead. Here’s why:
Open Questions in Conventional Eccentricity Theory
The Milankovitch eccentricity cycles are well-established in the literature, but several open questions remain:
| Question | Details |
|---|---|
| 1. The “~100k” simplification | Milankovitch’s actual calculations give ~95k and ~125k cycles. The commonly cited “~100k” is a rounded average that does not correspond to a single physical cycle. |
| 2. The 100,000-year problem | Geological temperature records show a dominant ~100k pattern, but no clear ~400k periodicity — despite the ~400k eccentricity cycle being the strongest in theory. This is a recognized unsolved problem in paleoclimatology. |
| 3. The energy problem | Eccentricity changes affect total annual insolation by only ~0.2%. How this small signal drives major glacial cycles remains debated — most proposals invoke amplification mechanisms (ice-albedo feedback, CO₂ feedbacks). |
| 4. Modeled vs. observed | The ~95k, ~125k, and ~400k cycles are derived from Jupiter-Saturn gravitational resonance models, not directly measured in the geological record. The match between theory and observation is approximate. |
| 5. Inclination precession | Earth’s orbital inclination precesses at ~67k years (vs ecliptic) or ~111k years (vs ICRF). This cycle was not part of Milankovitch’s original framework and is not included in standard eccentricity calculations, though it may contribute to the observed ~100k signal. |
Status as of 2025: The 100,000-year problem remains actively debated and unsolved. Barker et al. (2025, Science) investigated the distinct roles of precession, obliquity, and eccentricity in Pleistocene glacial cycles — still unable to resolve which parameter dominates. The Mid-Pleistocene Transition — the shift from 41-kyr to ~100-kyr glacial cycles around 1 million years ago — remains “one of paleoclimatology’s great unsolved puzzles.”
Notably, Muller & MacDonald’s 1997 PNAS paper showed that the spectral shape of the ~100k climate signal is incompatible with eccentricity’s split-peak spectrum (95k + 125k). This spectral evidence has never been refuted — only Muller’s proposed dust mechanism was rejected.
Independent dating methods (speleothem U-Th dating, O₂/N₂ ratio dating) exist that could test whether the true period is closer to 100k or 111k without circular orbital tuning. See Supporting Evidence for details.
See Scientific Background: Eccentricity Cycles and Milankovitch Theory for detailed analysis of these issues and the “100,000-year problem.”
The ~100k Pattern in Ice Cores
Ice core data does show a roughly ~100,000-year pattern in glacial cycles. The model proposes this actually reflects the inclination precession cycle (~111,772 years), not eccentricity:
~100,000 years in ice cores ≈ 111,772 years (inclination precession)
The ~10% discrepancy may be due to dating uncertainties in ice core chronology. See Scientific Background: Ice Core Chronology for detailed analysis.
Comparison with Standard Formulas
The model’s eccentricity predictions are compared with polynomial formulas from Newcomb (1898), Harkness (1891), and Meeus (1998). All four converge at J2000 (e = 0.0167). Standard polynomials predict continued decrease toward ~0.01 by 20,000 AD, while this model predicts bounded oscillation within ~0.0140–~0.0167 with minimum at ~11,725 AD followed by increase.
Long-term predictions differ because:
- The model uses a single 20,957-year cycle with bounded oscillation
- Standard theory uses ~100k/400k cycles predicting continued decrease
- Direct measurements only exist for recent centuries
Model vs. Milankovitch: Eccentricity Over 300,000 Years
The graph below shows both predictions over 300,000 years — from 100,000 years into the future to 200,000 years into the past. The model’s 20,957-year cycle (blue) oscillates within a narrow, bounded range of ~0.0140–~0.0167 around a base value of 0.015386 (red). The standard Milankovitch eccentricity (grey) varies over much larger amplitudes (up to ~0.06) on ~100,000 and ~400,000-year timescales.
The two predictions diverge significantly: the model predicts eccentricity never leaves its narrow band, while standard theory predicts it has varied by a factor of ~4 over the past 200,000 years. Since direct measurements only cover recent centuries, neither prediction can be verified for deep time — but they offer testable, fundamentally different forecasts.
Saturn Coupling — an Additional Effect
The eccentricity curve shown above reflects only Earth’s own 20,957-year perihelion cycle. Saturn’s eccentricity is independently predicted by Law 5 — the global eccentricity balance equation determines Saturn’s value from the other seven planets to 0.27%. The two predictions are not independent: Saturn participates in the same balance system that includes Earth, so changes in Earth’s eccentricity feed into Saturn’s via Law 5.
The physical Earth–Saturn coupling in the model comes from Saturn being the only planet whose precession formula requires Earth’s time-varying obliquity and eccentricity as inputs (GROUP 15 terms; see Formulas), and from Saturn’s perihelion cycle ( = 41,915 years, ecliptic-retrograde) decomposing via Fibonacci addition into Jupiter’s perihelion (H/5) + Earth’s inclination precession (H/3) — see Law 6.
A note on axial precession. The Earth-Saturn connection is purely orbital. Saturn’s axial (spin-axis) precession has a period of ~1.8 million years, driven by a spin-orbit resonance with Neptune’s nodal mode (Saillenfest et al. 2021). This is unrelated to the Fibonacci timescale hierarchy.
Important: Both the model’s predictions and standard Milankovitch predictions for ancient/future eccentricity are theoretical. Neither can be directly verified for times before ~1900 AD.
Climate Implications
Eccentricity affects Earth’s climate through two mechanisms:
1. Total Annual Energy
Higher eccentricity gives Earth slightly more total annual solar energy. The orbit-averaged flux scales as 1/√(1−e²) — perihelion’s intense, close-range flux more than compensates for the longer time spent near aphelion.
| Eccentricity | Effect on Annual Insolation |
|---|---|
| Maximum (~0.0167) | ~0.014% more than circular |
| Minimum (~0.0140) | ~0.010% more than circular |
| Difference | ~0.004% |
This effect is small - too small alone to cause ice ages.
2. Seasonal Contrast
The more important effect is when perihelion occurs relative to seasons:
| Perihelion Timing | Northern Hemisphere Effect |
|---|---|
| January (current) | Milder winters, cooler summers |
| July (~11,725 AD) | Hotter summers, colder winters |
When perihelion occurs during Northern Hemisphere winter (as now), winters are slightly milder. When it occurs during summer, seasonal contrasts increase.
Summary
| Aspect | Value |
|---|---|
| Current eccentricity | 0.01671022 (decreasing) |
| Cycle period | 20,957 years |
| Range | ~0.0140 to ~0.0167 |
| Maximum alignment | Perihelion at December solstice |
| Minimum alignment | Perihelion at June solstice |
| Last maximum | 1246 AD |
| Next minimum | ~11,725 AD |
Eccentricity Cycles for Other Planets
The same two-counter-rotating-motion principle applies to every planet. Each planet has its own wobble period — the meeting frequency of its axial precession and ICRF perihelion precession. This is the period over which the planet’s eccentricity completes one full oscillation:
| Planet | Wobble period (eccentricity cycle) | H expression |
|---|---|---|
| Mercury | 31,935 yr | 2H/21 |
| Venus | 14,045 yr | 8H/191 |
| Earth | 20,957 yr | H/16 |
| Mars | 50,298 yr | 3H/20 |
| Jupiter | 62,385 yr | 8H/43 |
| Saturn | 16,559 yr | 4H/81 |
| Uranus | 33,532 yr | ~H/10 |
| Neptune | 26,825 yr | ~2H/25 |
For Earth, the wobble period coincides with the perihelion precession period (H/16) because Earth’s axial precession (H/13) and ICRF perihelion precession (H/3) meet at this rate. For other planets, the two component periods are different, so the wobble period is a derived beat frequency.
Important distinction: The wobble period (eccentricity cycle) is NOT the same as the perihelion ecliptic period. For Earth they happen to coincide because the beat frequency of axial precession (H/13) and ICRF perihelion (H/3) equals H/16 (since 13 + 3 = 16); for other planets they differ. The 3D simulation’s Grand Holistic Octave panel shows all six cycle types per planet (axial, perihelion ecliptic, ICRF, ascending node, obliquity, eccentricity) — each as an integer divisor of 8H.
Calculate Eccentricity at Any Year
To calculate eccentricity values for any year, see the Formulas page which provides the complete formulas.
Key Takeaways
- Eccentricity = orbital elongation - Currently 0.01671022, meaning ~5 million km difference between perihelion and aphelion
- 20,957-year cycle - From the meeting frequency of two counter-rotating motions
- Maximum at winter solstice alignment - When Earth and PERIHELION-OF-EARTH offsets add together
- Minimum at summer solstice alignment - When offsets partially cancel
- Currently decreasing - We passed maximum around 1246 AD
- Simpler than Milankovitch - One cycle (20,957 years) instead of multiple overlapping cycles (100k/400k)
Continue to Days & Years to learn how these cycles affect the length of our days and years.