Milankovitch Cycles
**Milutin Milankovitch** highlighted how slow changes in **Earth's orbital eccentricity**, **obliquity** (axial tilt), and **precession** of the equinoxes modulate the **distribution and seasonality** of **solar heating**, especially at **high northern latitudes** where ice sheets respond. These orbital bands (~100 kyr, ~41 kyr, ~19–23 kyr) appear in **paleoclimate** records, though **full** glacial cycles require **carbon cycle** and **ice-albedo** feedbacks not modeled here. The graph integrates a **deliberately simplified** mean daily insolation at **65°N** over one year using a declination pattern tied to obliquity, a **Keplerian** eccentricity proxy on distance, and a **precession** phase slider—pedagogy only.
Who it's for: Earth science and astronomy survey courses linking orbital mechanics to climate context.
Key terms
- Milankovitch cycles
- Eccentricity
- Obliquity
- Precession
- Insolation
- Ice ages
- Paleoclimate
How it works
**Milankovitch** (orbital) climate forcing bundles **eccentricity** (modulating Earth–Sun distance over the year), **obliquity** (tilt of the spin axis, changing seasonal contrasts), and **precession** (where seasons fall on the ellipse). Together they shift **high-latitude summer insolation**—a pacemaker for **ice-age** cycles on ~10–100 kyr scales when combined with feedbacks. This page plots a **highly simplified** mean daily insolation at **65°N** through one year using a circular-orbit **declination** pattern with your obliquity, scaled by an **eccentricity** proxy from Kepler’s **true anomaly**, and shifted by a **precession** phase. It is **pedagogy**, not a coupled climate model.
Key equations
Frequently asked questions
- Does this simulator predict the next ice age?
- No. It illustrates orbital forcing magnitudes, not coupled ice-sheet or CO₂ dynamics.
- Why emphasize 65°N?
- Northern high-latitude summer insolation is a common index in Milankovitch literature because of Laurentide ice-sheet sensitivity; it is a convention, not a universal truth for all glaciers.
More from Astronomy & The Sky
Other simulators in this category — or see all 35.
Supernova Light Curves
Schematic Ia rise/decay vs II-P plateau; standard-candle note.
Apparent vs Absolute Magnitude
m = M + 5 log₁₀(d/10 pc); distance modulus; flux ratios.
Spin–Orbit Resonance
Moon 1:1 lock vs Mercury 3:2; schematic animations.
Pulsar Lighthouse
Rotating beam cone, pulse profile; timing / ms pulsars context.
Meteor Shower & Radiant
Earth crosses comet debris; radiant on a star field (schematic).
Cosmological Expansion (FLRW)
a(t), z, χ and c/H vs time; flat Ω_m + Λ (toy ΛCDM).