Orbital variations (eccentricity, obliquity, and precession) alter the distribution and timing of solar radiation on Earth's surface over timescales of 20–400 ka. These variations are caused by gravitational interactions among Earth, Moon, and other planets and are predictable for millions of years using orbital mechanics. Milankovitch theory proposes that orbital forcing is the primary driver of glacial-interglacial cycles, a hypothesis strongly supported by paleoclimate records showing clear orbital periodicities in ice volume, temperature, and atmospheric CO₂.
Calculate how orbital parameters vary with time using published orbital solutions (Laskar et al.). Compute the seasonal and latitudinal distribution of insolation and see how glacial sensitivity varies by latitude and season.
The total solar energy received by Earth barely changes with orbital variations (<0.1%); what matters is the DISTRIBUTION of insolation. Also, orbital forcing is necessary but not sufficient for glaciation; internal climate feedbacks amplify the small insolation signal.
From Kepler's laws you know that Earth's orbit is not a perfect circle — it is a slightly elliptical path around the Sun, and the planet moves faster near perihelion (closest approach) and slower near aphelion. From Earth's energy balance you know that the amount and distribution of solar radiation drives climate. Milankovitch theory connects these two ideas: the shape of Earth's orbit, the tilt of its axis, and the wobble of that tilt all change slowly and predictably, and these changes alter how solar radiation is distributed across latitudes and seasons.
There are three orbital parameters. Eccentricity describes how elliptical the orbit is, cycling from nearly circular to slightly more elliptical and back with a dominant period of about 100,000 years (and a secondary ~400,000-year cycle). When eccentricity is higher, the difference in Earth-Sun distance between perihelion and aphelion is greater, creating stronger insolation asymmetry between the two halves of the year. Obliquity is the tilt of Earth's rotational axis relative to its orbital plane, cycling between about 22.1° and 24.5° with a period of roughly 41,000 years. Higher obliquity means more extreme seasons — hotter summers and colder winters at high latitudes. Precession refers to the slow wobble of Earth's rotation axis, completing a cycle in about 23,000 years. Precession determines which hemisphere is tilted toward the Sun at perihelion, affecting whether Northern or Southern Hemisphere summers coincide with Earth being closest to the Sun.
Here is the crucial misconception to avoid: orbital variations change total annual solar input by less than 0.1%. If you added up all the sunlight Earth receives in a year across all latitudes, it barely changes. What changes is the *distribution* — specifically, how much sunlight reaches high latitudes in summer. The leading theory for why this matters is that glaciations grow when high northern latitudes (where most of Earth's landmass lies) receive insufficient summer sunlight to melt the previous winter's snow. Ice that survives summer persists, reflects more sunlight, and the ice sheet grows. The relevant signal is summer insolation at ~65°N, not global annual total.
Even with the right orbital forcing, the observed magnitude of glacial-interglacial cycles — roughly 5–8°C in global mean temperature, with polar regions changing by far more — requires amplification. This is where feedbacks come in. As initial cooling allows ice to spread, Earth's albedo increases (ice-albedo feedback), reducing absorbed solar radiation further. As the ocean cools, it absorbs more CO₂ from the atmosphere, reducing the greenhouse effect (CO₂ feedback). Both are positive feedbacks that amplify the initial orbital nudge into a full glaciation. Orbital forcing provides the *pacing* and *timing* of ice ages; feedbacks provide the *amplitude*.
The evidence for Milankovitch theory is compelling and comes primarily from paleoclimate proxies: deep-sea sediment cores and Antarctic ice cores record oxygen isotope ratios (reflecting ice volume and temperature), CO₂, and dust over hundreds of thousands of years. Spectral analysis of these records reveals strong periodicities at exactly the orbital frequencies — ~100 ka, ~41 ka, and ~23 ka — matching what orbital mechanics predicts. This correspondence between astronomical theory and geological data, assembled decades after Milankovitch's original calculations, is one of the great confirmations in Earth science.