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 3 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

Named alongside it

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

EccentricityObliquityAngular momentum deficitAsteroid familyCaloric half-yearCassini stateChaosDisturbing functionDynamical ellipticityEigenmodeForced eccentricityInsolation

All concepts