Obliquity pacing of the late Pleistocene glacial terminations

The 100,000-year timescale in the glacial/interglacial cycles of the late Pleistocene epoch (the past approximately 700,000 years) is commonly attributed to control by variations in the Earth's orbit. This hypothesis has inspired models that depend on the Earth's obliquity (approximately 40,000 yr; approximately 40 kyr), orbital eccentricity (approximately 100 kyr) and precessional (approximately 20 kyr) fluctuations, with the emphasis usually on eccentricity and precessional forcing. According to a contrasting hypothesis, the glacial cycles arise primarily because of random internal climate variability. Taking these two perspectives together, there are currently more than thirty different models of the seven late-Pleistocene glacial cycles. Here we present a statistical test of the orbital forcing hypothesis, focusing on the rapid deglaciation events known as terminations. According to our analysis, the null hypothesis that glacial terminations are independent of obliquity can be rejected at the 5% significance level, whereas the corresponding null hypotheses for eccentricity and precession cannot be rejected. The simplest inference consistent with the test results is that the ice sheets terminated every second or third obliquity cycle at times of high obliquity, similar to the original proposal by Milankovitch. We also present simple stochastic and deterministic models that describe the timing of the late-Pleistocene glacial terminations purely in terms of obliquity forcing.