A simple rule to determine which insolation cycles lead to interglacials

The pacing of glacial–interglacial cycles during the Quaternary period (the past 2.6 million years) is attributed to astronomically driven changes in high-latitude insolation. However, it has not been clear how astronomical forcing translates into the observed sequence of interglacials. Here we show that before one million years ago interglacials occurred when the energy related to summer insolation exceeded a simple threshold, about every 41,000 years. Over the past one million years, fewer of these insolation peaks resulted in deglaciation (that is, more insolation peaks were ‘skipped’), implying that the energy threshold for deglaciation had risen, which led to longer glacials. However, as a glacial lengthens, the energy needed for deglaciation decreases. A statistical model that combines these observations correctly predicts every complete deglaciation of the past million years and shows that the sequence of interglacials that has occurred is one of a small set of possibilities. The model accounts for the dominance of obliquity-paced glacial–interglacial cycles early in the Quaternary and for the change in their frequency about one million years ago. We propose that the appearance of larger ice sheets over the past million years was a consequence of an increase in the deglaciation threshold and in the number of skipped insolation peaks.

[1]  M. Maslin,et al.  The role of orbital forcing in the Early Middle Pleistocene Transition , 2015 .

[2]  Maureen E. Raymo,et al.  Matuyama 41,000-year cycles: North Atlantic Ocean and northern hemisphere ice sheets , 1986 .

[3]  N. Opdyke,et al.  Oxygen-Isotope and Paleomagnetic Stratigraphy of Pacific Core V28-239 Late Pliocene to Latest Pleistocene , 1976 .

[4]  André Berger,et al.  On the Structure and Origin of Major Glaciation Cycles .2. the 100,000-year Cycle , 1993 .

[5]  Barry Saltzman,et al.  A first-order global model of late Cenozoic climatic change II. Further analysis based on a simplification of CO2 dynamics , 1991, Transactions of the Royal Society of Edinburgh: Earth Sciences.

[6]  K. Maasch.  Statistical detection of the mid-Pleistocene transition , 1988 .

[7]  A. Ganopolski,et al.  The role of orbital forcing, carbon dioxide and regolith in 100 kyr glacial cycles , 2011 .

[8]  A. Abe‐Ouchi,et al.  Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume , 2013, Nature.

[9]  R. Edwards,et al.  The Asian monsoon over the past 640,000 years and ice age terminations , 2016, Nature.

[10]  P. Gibbard,et al.  Early–Middle Pleistocene transitions: Linking terrestrial and marine realms , 2015 .

[11]  H. Elderfield,et al.  Evolution of Ocean Temperature and Ice Volume Through the Mid-Pleistocene Climate Transition , 2012, Science.

[12]  V. Brovkin,et al.  Can we predict the duration of an interglacial , 2012 .

[13]  Albert Lunde,et al.  A Paleoclimate Model of Northern Hemisphere Ice Sheets , 1981, Quaternary Research.

[14]  D. V. Husen.  Chapter 2 – Quaternary Glaciations in Austria , 2011 .

[15]  Didier Paillard,et al.  The timing of Pleistocene glaciations from a simple multiple-state climate model , 1998, Nature.

[16]  M. Raymo,et al.  The 41 kyr world: Milankovitch's other unsolved mystery , 2003 .

[17]  Milutin Milankovictch,et al.  Canon of insolation and the ice-age problem : (Kanon der Erdbestrahlung und seine Anwendung auf das Eiszeitenproblem) Belgrade, 1941. , 1998 .

[18]  L. Lisiecki.  Atlantic overturning responses to obliquity and precession over the last 3 Myr , 2014 .

[19]  Annabel M. Imbrie-Moore,et al.  A phase-space model for Pleistocene ice volume , 2010, 1104.3610.

[20]  P. Huybers.  Combined obliquity and precession pacing of late Pleistocene deglaciations , 2011, Nature.

[21]  N. Shackleton,et al.  Distribution, stratigraphic position and age of ash layer , 1975 .

[22]  M. Loutre,et al.  Modelling northern hemisphere ice volume over the last 3 Ma , 1999 .

[23]  G. Haug,et al.  Strengthening of North American dust sources during the late Pliocene (2.7 Ma) , 2012 .

[24]  Peter Huybers,et al.  Early Pleistocene Glacial Cycles and the Integrated Summer Insolation Forcing , 2006, Science.

[25]  F. Parrenin,et al.  Terminations VI and VIII (∼ 530 and ∼ 720 kyr BP) tell us the importance of obliquity and precession in the triggering of deglaciations , 2012 .

[26]  Joseph G Ibrahim,et al.  Properties and Implementation of Jeffreys’s Prior in Binomial Regression Models , 2008, Journal of the American Statistical Association.

[27]  N. Shackleton.  Oxygen Isotope Analyses and Pleistocene Temperatures Re-assessed , 1967, Nature.

[28]  N. Shackleton.  Pliocene stable isotope stratigraphy of site 846 , 1995 .

[29]  P. Ashwin,et al.  The middle Pleistocene transition as a generic bifurcation on a slow manifold , 2015, Climate Dynamics.

[30]  Wallace S. Broecker,et al.  The Carbon cycle and atmospheric CO[2] : natural variations Archean to present , 1985 .

[31]  Eli Tziperman,et al.  On the mid-Pleistocene transition to 100-kyr glacial cycles and the asymmetry between glaciation and deglaciation times , 2003 .

[32]  K. Lawrence,et al.  North Atlantic climate evolution through the Plio-Pleistocene climate transitions , 2010 .

[33]  Jacques Laskar,et al.  A long-term numerical solution for the insolation quantities of the Earth , 2004 .

[34]  R. Wilkinson,et al.  Global sensitivity analysis of the climate–vegetation system to astronomical forcing: an emulator-based approach , 2014 .

[35]  Peter John Huybers,et al.  Glacial variability over the last two million years: an extended depth-derived agemodel, continuous obliquity pacing, and the Pleistocene progression , 2007 .

[36]  D. Hodell,et al.  Late Neogene history of deepwater ventilation in the Southern Ocean , 2006 .

[37]  Gerald H. Haug,et al.  Southern Ocean dust–climate coupling over the past four million years , 2011, Nature.

[38]  V. Brovkin,et al.  Interglacials of the last 800,000 years , 2016 .

[39]  André Berger,et al.  An alternative astronomical calibration of the lower Pleistocene timescale based on ODP Site 677 , 1990, Transactions of the Royal Society of Edinburgh: Earth Sciences.

[40]  J. D. Hays,et al.  Variations in the Earth ' s Orbit : Pacemaker of the Ice Ages Author ( s ) : , 2022 .

[41]  Yair Rosenthal,et al.  Pleistocene sea-surface temperature evolution: Early cooling, delayed glacial intensification, and implications for the mid-Pleistocene climate transition , 2013 .

[42]  Nicholas J Shackleton,et al.  Oxygen Isotope and Palaeomagnetic Stratigraphy of Equatorial Pacific Core V28-238: Oxygen Isotope Temperatures and Ice Volumes on a 105 Year and 106 Year Scale , 1973, Quaternary Research.

[43]  Fridtjof Nansen,et al.  The polar oceans and their role in shaping the global environment : the Nansen centennial volume , 1994 .

[44]  Carl Wunsch,et al.  Consequences of pacing the Pleistocene 100 kyr ice ages by nonlinear phase locking to Milankovitch forcing , 2006 .

[45]  P. Gibbard,et al.  Quaternary Glaciations - Extent and Chronology: A closer look , 2011 .

[46]  D. Pollard.  A coupled climate-ice sheet model applied to the Quaternary Ice Ages , 1983 .

[47]  Maureen E. Raymo,et al.  The timing of major climate terminations , 1997 .

[48]  D. Hodell,et al.  Mode transitions in Northern Hemisphere glaciation: co-evolution ofmillennial and orbital variability in Quaternary climate , 2016 .

[49]  C. Robert,et al.  Computational methods for Bayesian model choice , 2009, 0907.5123.

[50]  N. Pisias,et al.  THE EVOLUTION OF PLEISTOCENE CLIMATE: A TIME SERIES APPROACH , 1981 .

[51]  Didier Paillard,et al.  The Antarctic ice sheet and the triggering of deglaciations , 2004 .

[52]  D. Macayeal.  A Catastrophe Model of the Paleoglimate , 1979, Journal of Glaciology.

[53]  David Archer,et al.  The middle Pleistocene transition: characteristics, mechanisms, and implications for long-term changes in atmospheric pCO2 , 2006 .

[54]  L. Bopp,et al.  Impact of oceanic processes on the carbon cycle during the last termination , 2011 .

[55]  E. Wolff,et al.  Interglacial and glacial variability from the last 800 ka in marine, ice and terrestrial archives , 2010 .

[56]  M. Raymo,et al.  A Pliocene‐Pleistocene stack of 57 globally distributed benthic δ18O records , 2005 .