Astronomy & Celestial Mechanics
the sky, computed
Concept

Milankovitch cycles — where it appears

The periodic changes in a planet's obliquity, orbital eccentricity and orientation, and the climate variations they pace. On Earth the obliquity term has a period of about 41,000 years and dominates the older part of the ice-age record.

Named by 8 essays across 2 fields — each of them below, with the objects they name alongside it.

Earth's eccentricity is a sum of 8 sinusoids. The eccentricity of Earth over 800 thousand years, from the Laplace–Lagrange solution for all eight planets — the secular matrix built from the JPL masses and semi-major axes, symmetrised, and diagonalised by Jacobi rotations. It runs between 0.0035 and 0.0436, and it has no period, because it is a sum of 8 incommensurable frequencies. The two largest contributions to this planet are the modes at 3.73 and 7.33 arcseconds per year, drawn as the flat lines: those are constants, and everything moving in the figure is their beat. The check is the quadratic form ½ΣΛe², which a symmetric secular matrix conserves exactly and which drifts by 6.7e-16 across the whole interval — computed from the same curves the figure draws and from nothing the eigenvalues were fitted to. The true angular momentum deficit, Σ Λ(1 − √(1−e²)), drifts by 7.1e-4, and that difference is not an error either: the two agree only to fourth order in e, and Mercury at 0.206 supplies almost all of the gap.

No planet has an eccentricity of its own

Strip the short-period terms out of the planetary equations and what is left is a linear system. Its eigenvectors are modes of the whole solar system, and the number a catalogue quotes for a planet's eccentricity turns out to be a reading of a clock rather than a property of the planet.

orbits · Secular theory Two equilibria become four, and three worlds sit near the join. The Cassini equilibria of a spin axis, drawn against the ratio of its own precession rate to the rate at which its orbit plane turns, for an orbit inclination of 1.5 degrees. Each column of dots is the full set of obliquities at which the two precessions keep step at that ratio, found by root-finding rather than by tracing a remembered curve. Below α cos ε/|g| = 1.135 there are two such obliquities and above it there are four, and the figure checks both counts on either side of the join. The three marked bodies are placed by their own measured precession constants: the Earth with the Moon at 2.67, safely on the four-state side; the Earth without it at 0.86; and Mars at 1.06. Two of the three sit within a few tenths of the bifurcation, which is the whole reason their obliquities are not constants: near the join the equilibria are close together, the libration around them is wide, and a body pushed between neighbouring resonances wanders. The Moon's contribution to the Earth's precession constant is what moves the first mark away from that region, and the second mark is the same planet with that contribution removed. This is a two-frequency model of a many-frequency system, and the real chaos comes from the overlap of resonances it does not contain.

A tilt that is not a constant

The Earth's axis leans by 23.4 degrees, and that lean is what makes the seasons. It is also a dynamical variable with its own equilibria, its own resonances and its own chaos — and on Mars the same variable has swung between nearly zero and sixty degrees without anything having to happen.

sky · Obliquity June sunlight at 65°N, against where perihelion sits. Daily-mean insolation at latitude 65 degrees north on the June solstice, at an obliquity of 23.44 degrees, against the longitude of perihelion measured from the March equinox — the angle that precesses right round in about twenty-one thousand years. Three eccentricities are drawn. The swing is ±10.0 per cent at e = 0.05 and ±1.0 per cent at e = 0.005, in proportion to the eccentricity, because the Sun–Earth distance on a fixed date carries e cos of the precession angle and that is first order. Over the same range of eccentricity the annual mean at this latitude moves by 0.12 per cent, because the annual mean carries 1/√(1−e²) and that is second order. The two together are the whole of the precession term in Milankovitch's theory: the eccentricity does almost nothing to how much sunlight the Earth receives and a great deal to when it arrives, and the ice sheets of the northern hemisphere respond to the summer they might melt in rather than to the year's total. It also explains why the precession signal disappears when the orbit is nearly circular: multiply a large angular swing by a vanishing eccentricity and there is nothing left, which is what the innermost curve here is.

An average that precession cannot move

Sunlight arrives as the inverse square of the distance and time passes as its square, so the two cancel exactly in a year's integral. The longitude of perihelion therefore changes the annual mean insolation at every latitude by precisely nothing — and changes June at 65°N by ten per cent.

orbits · Orbital averages A hotter summer that is also a shorter one. Daily-mean sunlight at 65°N through one year beginning at the March equinox, for an orbit of eccentricity 0.0167 and obliquity 23.44°, drawn twice: with perihelion at the June solstice, when the northern summer is closest to the Sun, and with perihelion at the December solstice, half a precession cycle later. The first summer peaks 6.9% higher, and it is also shorter, because the Earth moves fastest near perihelion: days above 275 W m⁻² number 139 against 144. The energy delivered on those days — the shaded area above the line — is 5.055 against 4.989 GJ m⁻², a difference of 1.3%. The second law makes the intensity and the duration trade against each other almost exactly, so a summer measured by how much energy it delivers barely feels precession at all.

A hotter summer that is also a shorter one

When the Earth passes closest to the Sun in northern summer, the northern summer is brighter — and, by Kepler's second law, shorter by exactly the right amount to deliver almost the same energy. Measured by the energy it delivers above a melting threshold rather than by its brightest day, a summer barely feels the precession cycle at all and feels the tilt of the axis several times more strongly.

orbits · Orbital averages The same tilt, and ice ages of two different lengths. Top: the obliquity anomaly in degrees from the five largest terms of the standard series, over 2 million years of model time; it is dominated by the 41,000-year cycle and has no power at all near 100,000 years. Bottom: ice volume in a threshold model driven by it — the ice grows steadily, and it collapses whenever volume plus a multiple of the obliquity anomaly crosses a threshold, marked by dots. For the first 1 million years the threshold is low and the ice collapses at nearly every obliquity maximum, on average every 43 thousand years. At the vertical line the threshold is raised, and nothing else changes: the ice now survives one or two obliquity maxima and collapses at the second or third, every 108 thousand years on average. The forcing is identical on both sides; the ice sheet counts it differently.

An ice sheet that skips a beat of the tilt

The last eight hundred thousand years of ice ages come roughly every hundred thousand years, and the Earth's orbit has no forcing of any strength at that period. A model in which an ice sheet collapses only when it is both large and warmed by a high-obliquity summer reproduces the rhythm with nothing but the forty-one-thousand-year tilt — the ice sheet grows through one or two maxima and falls on the next — and a slow rise in how large the ice must grow turns the older forty-one-thousand-year world into the newer one without any change in the orbit.

orbits · Orbital averages The summer sunshine at Mars's north pole over the last million years. The sunshine reaching Mars's north pole at the summer solstice, relative to its average, over the last 1,000 thousand years (blue), from a forcing built from the main periods of Mars's recent orbit: the tilt varying by about 3.5 degrees in 120,000 years and 1 degree in 99,000, the eccentricity by 0.03 in 95,000 years and 0.02 in 2.4 million, and the perihelion going round the seasons once in 51,000 years. The part due to the tilt alone is dashed. The sunshine swings by a factor of 2.06 between its extremes: the tilt sets a slow swell with a period of 120,000 years, and the precession, whose strength the eccentricity sets, cuts a faster ripple of 51,000 years into it, strongest when the orbit is most eccentric. These two periods, and their ratio of 2.35, are what a record of the sunshine should carry.

Layers that keep the rhythm and lose the dates

Mars's polar caps are stacks of ice and dust laid down while the planet's tilt and orbit changed the summer sunshine, and the question is whether they can be read as a record of that sunshine. They can — for the rhythm. Even when the ice accumulates faster under weak sunshine than under strong, the cycles of the tilt and the precession survive in the depth of the stack in their correct proportion. What does not survive is the date of any one layer.

sky · Obliquity The rhythms in the rock, set by the Moon's distance. The periods of the climatic precession cycles — the precession constant plus the perihelion frequency of each of the planets that dominate the Earth's eccentricity — and of the obliquity cycle, the precession constant plus the node frequency s₃, against the Earth–Moon distance. Today they are 23.8, 22.4, 19.0 thousand years for precession and 41 for obliquity. At 341,000 km they were 14.7, 14.2, 12.8 and 20.0. The planetary frequencies are held at today's values, which the chaos of the inner solar system makes uncertain beyond a few tens of millions of years; the precession constant is the term that moves.

A rhythm in the rock that measures the Moon

The Earth's axis precesses because the Sun and Moon pull on its equatorial bulge, and two-thirds of the pull is the Moon's. When the Moon was nearer and the day shorter, the axis turned faster, and the climate cycles that precession drives were shorter. Sediments laid down 1.4 billion years ago record those cycles — and with nothing in the rock but layer thicknesses, the ratio of two rhythms gives the Moon's distance at the time.

sky · Precession Two old points, and the history they draw back to the Moon's birth. The Earth–Moon distance against time before the present, with the estimates from tidal rhythmites (620 million years), cyclic sediments (1.4 billion) and banded iron formation (2.46 billion, 321.8 thousand km, ±6.2). Today's tidal efficiency, scaled as the inverse sixth power of the distance, brings the Moon to the Earth 1.54 billion years ago. The efficiency fitted between the two oldest points — 21% of today's — held further back brings it out from the Earth 4.78 billion years ago, close to the Moon's age of about 4.5 billion. The constant law is an assumption; what the points measure is the average between them.

A Moon measured in banded iron

Layers of iron-rich sediment laid down 2.46 billion years ago in what is now Western Australia keep time with the Earth's precession, and the precession's period says how far away the Moon was: about 321,800 kilometres, with a day of seventeen hours. Together with an older-is-closer point from 1.4 billion years ago, it measures how hard the tides worked between them — about a fifth as hard as today. Held further back, that efficiency brings the Moon out from the Earth a few hundred million years too early, which says the young tides worked harder than the middle-aged ones.

sky · Precession

Named alongside it

The objects these essays reach for when they reach for this one.

ObliquityEccentricityInsolationPrecessionPrecession of the equinoxesAngular momentumCaloric half-yearCyclostratigraphyDynamical ellipticityIce agesLength of dayLongitude of perihelion

All concepts