Skip to Content
📄 Fibonacci Laws — Read the paper
The ModelEccentricity

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).

MeasurementValue
1 AU (mean Earth-Sun distance)149,597,870.698828 km
Eccentricity (J2000)0.01671022
Offset distance2,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

Earth and PERIHELION-OF-EARTH orbiting in opposite directions

Two motions work in opposite directions:

MotionDirectionPeriod
Earth around EARTH-WOBBLE-CENTERClockwise~25,794 years
PERIHELION-OF-EARTH around SunCounter-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

ParameterValueNotes
Current eccentricity (J2000)0.01671022Measured, NASA Planetary Fact Sheet
Base eccentricity0.015386Model-derived arithmetic midpoint (derived mean = 0.0154454).
Maximum eccentricity~0.0167At winter solstice alignment
Minimum eccentricity~0.0140At summer solstice alignment
Variation amplitude±0.001356Half the range
Cycle period20,957 years335,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:

  1. Minimum eccentricity occurred in ~9,233 BC when perihelion aligned with the June solstice
  2. Maximum eccentricity occurred in 1246 AD when perihelion aligned with the December solstice
  3. 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:

AlignmentEccentricityLast OccurrenceNext Occurrence
Perihelion at December solsticeMaximum (~0.0167)1246 AD~22,203.344 AD
Perihelion at June solsticeMinimum (~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.

Graph showing Earth's eccentricity oscillating between minimum and maximum over one eccentricity cycle, with maximum at December solstice alignment

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:

QuestionDetails
1. The “~100k” simplificationMilankovitch’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 problemEccentricity 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. observedThe ~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 precessionEarth’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.

Eccentricity predictions compared: this model (purple) versus Newcomb (1898, blue), Harkness (1891, red), and Meeus (1998, green). All converge at J2000.

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.

Graph comparing the model's bounded eccentricity cycle (blue) against the standard Milankovitch eccentricity variation (grey) over several hundred thousand years

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 (H/8-H/8 = 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.

EccentricityEffect 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 TimingNorthern 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

AspectValue
Current eccentricity0.01671022 (decreasing)
Cycle period20,957 years
Range~0.0140 to ~0.0167
Maximum alignmentPerihelion at December solstice
Minimum alignmentPerihelion at June solstice
Last maximum1246 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:

PlanetWobble period (eccentricity cycle)H expression
Mercury31,935 yr2H/21
Venus14,045 yr8H/191
Earth20,957 yrH/16
Mars50,298 yr3H/20
Jupiter62,385 yr8H/43
Saturn16,559 yr4H/81
Uranus33,532 yr~H/10
Neptune26,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

  1. Eccentricity = orbital elongation - Currently 0.01671022, meaning ~5 million km difference between perihelion and aphelion
  2. 20,957-year cycle - From the meeting frequency of two counter-rotating motions
  3. Maximum at winter solstice alignment - When Earth and PERIHELION-OF-EARTH offsets add together
  4. Minimum at summer solstice alignment - When offsets partially cancel
  5. Currently decreasing - We passed maximum around 1246 AD
  6. 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.

Last updated on: