Early Pleistocene Glacial Cycles and the Integrated Summer Insolation Forcing
Abstract
Long-term variations in Northern Hemisphere summer insolation are generally thought to control glaciation. But the intensity of summer insolation is primarily controlled by 20,000-year cycles in the precession of the equinoxes, whereas early Pleistocene glacial cycles occur at 40,000-year intervals, matching the period of changes in Earth's obliquity. The resolution of this 40,000-year problem is that glaciers are sensitive to insolation integrated over the duration of the summer. The integrated summer insolation is primarily controlled by obliquity and not precession because, by Kepler's second law, the duration of the summer is inversely proportional to Earth's distance from the Sun.
Get full access to this article
View all available purchase options and get full access to this article.
Already a Subscriber?Sign In
Supplementary Material
File (huybers.som.pdf)
References and Notes
1
J. Adhémar, Révolutions de la Mer (Carilian-Goeury et V. Dalmont, Paris, 1842).
2
E. Bard, C. R. Geosci.336, 603 (2004).
3
J. Croll, Philos. Mag.28, 121 (1864).
4
M. Milankovitch, Kanon der Erdbestrahlung und seine Andwendung auf das Eiszeitenproblem (Royal Serbian Academy, Belgrade, 1941).
5
J. D. Hays, J. Imbrie, N. J. Shackleton, Science194, 1121 (1976).
6
J. Imbrie, J. Z. Imbrie, Science207, 943 (1980).
7
J. Imbrie et al., Paleoceanography8, 699 (1993).
8
J. Imbrie et al., Paleoceanography7, 701 (1992).
9
R. Tiedemann, M. Sarnthein, N. J. Shackleton, Paleoceanography9, 619 (1994).
10
L. Lisiecki, M. Raymo, Paleoceanography20,2005).
11
M. Raymo, K. Nisancioglu, Paleoceanography18,2003).
12
K. Nisancioglu, thesis, Massachusetts Institute of Technology (2004).
13
P. Clark, R. Alley, D. Pollard, Science286, 1104 (1999).
14
E. Tziperman, H. Gildor, Paleoceanography18,2003).
15
D. Paillard, Nature391, 378 (1998).
16
P. Huybers, C. Wunsch, Nature434, 491 (2005).
17
A. Berger, X. Li, M. Loutre, Quat. Sci. Rev.18, 1 (1999).
18
R. DeConto, D. Pollard, Palaeogeogr. Palaeoclimatol. Palaeoecol.198, 39 (2003).
19
A. Berger, M. F. Loutre, Earth Planet. Sci. Lett.111, 369 (1992).
20
Daily average surface temperatures are estimated by using the network of 8892 World Meteorological Organization (WMO) stations above 30°N for the years 1994 to 1999. All stations that have greater than 80% data coverage are used. Data gaps are filled by linear interpolation. Stations are standardized to 1 km of elevation assuming a lapse rate of 6.5°C/km, were binned according to 1° or 10° latitude bands (as indicated in the text), and are then averaged together. Lastly, each of the six consecutive seasonal cycles are averaged together, yielding average annual cycles as a function of latitude.
21
C. Wunsch, J. Clim.18, 4374 (2005).
22
W. Paterson, Physics of Glaciers (Pergamon Press, Oxford, ed. 3, 1994).
23
R. Braithwaite, Y. Zhang, J. Glaciol.152, 7 (2000).
24
The use of a constant value for τ illustrates the concept of summer energy. A more detailed description would take into account that τ is expected to be spatially and temporally variable, depending on factors such as elevation, albedo, clouds, heat transport, and greenhouse gas concentrations. Note, however, that results are not sensitive to plausible choices of τ and that values less than 325 W/m2 yield broadly consistent summer energies (fig. S1). Summer energy values at 65°N are given in table S1.
25
The relationship between insolation intensity and insolation energy is more precisely illustrated by noting that I ∝ 1/r2, where I is insolation intensity and r is the distance from the Earth to the Sun. Conservation of angular momentum (or, equivalently, Kepler's second law) dictates that dt ∝ r2dλ, where dt is an infinitesimal change in time and dλ the corresponding change in solar longitude. The energy received by the Earth is then J = Idt ∝ dλ. In contrast with I, the J between any two solar longitudes is independent of r and, thus, independent of the precession of the equinoxes.
26
Are past changes in summer energy sufficient to cause the waxing and waning of ice sheets? Although a full answer requires a realistic model of Pleistocene climate, some indication is provided by modern glacial changes: A 2.4 W/m2 global average increase in radiative forcing caused by greenhouse gases (34) has apparently led to a general decrease in glacial mass (35), suggesting that glaciers are sensitive to relatively small changes in the radiation budget.
27
C. Wunsch, Clim. Dyn.20, 353 (2003).
28
Materials and Methods are available as supporting material on Science Online.
29
Amplitude cross correlation was computed by pairing local maxima in insolation with the nearest (in time) maximum in the rate of change of ice volume. Before identifying maxima, both the δ18O record and the summer energy were smoothed by using an 11-ky tapered window. There are 34 local maxima in summer energy at 65°N between 2 and 1 My ago and another 34 between 1 My ago and the present. Squared cross correlations of 0.4 and higher have P values of less than 0.01. Spectral and coherence analysis is performed by using Thomson's multitaper method (36).
30
J. C. Zachos, N. J. Shackleton, J. S. Revenaugh, H. Pälike, B. P. Flower, Science292, 274 (2001).
31
Similar with the early Pleistocene, late Pleistocene ice volume change has significant variability concentrated at the obliquity band, which is in phase and coherent with summer energy (P = 0.01). That the obliquity component of summer energy varies symmetrically between the hemispheres helps explain the symmetry of glacial variations between the hemispheres. Also, the increase in summer energy near 420 ky ago, absent in measures of summer insolation forcing relying on intensity, helps explain the corresponding glacial termination.
32
M. E. Raymo, Paleoceanography12, 577 (1997).
33
M. Raymo, Annu. Rev. Earth Planet. Sci.22, 353 (1994).
34
J. Houghtonet al., Eds., Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge Univ. Press, New York, 2001).
35
J. Oerlemans, Science308, 675 (2005); published online 3March 2005 (
36
D. Percival, A. Walden, Spectral Analysis for Physical Applications (Cambridge Univ. Press, Cambridge, 1993).
37
This paper benefited from discussion with E. Boyle, B. Curry, M. Raymo, P. Stone, E. Tziperman, and C. Wunsch. J. Levine provided valuable assistance in calculating the insolation. The NSF paleoclimate program supported this work under grant no. ATM-0455470.
Information & Authors
Information
Published In

Science
Volume 313 | Issue 5786
28 July 2006
28 July 2006
Copyright
American Association for the Advancement of Science.
Article versions
You are viewing the most recent version of this article.
Submission history
Received: 23 January 2006
Accepted: 9 June 2006
Published in print: 28 July 2006
Authors
Metrics & Citations
Metrics
Article Usage
Altmetrics
Citations
Export citation
Select the format you want to export the citation of this publication.
Cited by
- Organic carbon burial is paced by a ~173-ka obliquity cycle in the middle to high latitudes, Science Advances, 7, 28, (2021)./doi/10.1126/sciadv.abf9489
- Chronology and Eccentricity Phasing for the Early Turonian Greenhouse (∼93–94 Ma): Constraints on Astronomical Control of the Carbon Cycle, Paleoceanography and Paleoclimatology, 36, 4, (2021).https://doi.org/10.1029/2020PA004188
- Observationally Constrained Cloud Phase Unmasks Orbitally Driven Climate Feedbacks, Geophysical Research Letters, 48, 6, (2021).https://doi.org/10.1029/2020GL091873
- Climate and Astronomical Cycles, Paleoclimatology, (385-404), (2021).https://doi.org/10.1007/978-3-030-24982-3_28
- Spatiotemporal patterns of sediment deposition on the northern slope of the South China Sea in the last 150,000 years, Journal of Palaeogeography, 10, 1, (2021).https://doi.org/10.1186/s42501-021-00102-3
- Oscillatory reactions, Chemical Kinetics, (579-602), (2021).https://doi.org/10.1016/B978-0-444-64039-0.00008-X
- Duelling narratives of chironomids and pollen explain climate enigmas during The Last Glacial-Interglacial transition in North Island New Zealand, Quaternary Science Reviews, 263, (106997), (2021).https://doi.org/10.1016/j.quascirev.2021.106997
- Orbital forcing in southern Africa: Towards a conceptual model for predicting deep time environmental change from an incomplete proxy record, Quaternary Science Reviews, 265, (107050), (2021).https://doi.org/10.1016/j.quascirev.2021.107050
- Evolution of paleo-climate and seawater pH from the late Permian to postindustrial periods recorded by boron isotopes and B/Ca in biogenic carbonates, Earth-Science Reviews, 215, (103546), (2021).https://doi.org/10.1016/j.earscirev.2021.103546
- Late Pleistocene Emergence of Crystalline Canadian Shield Sources in Sediments of the Northern Gulf of Mexico, Paleoceanography and Paleoclimatology, 36, 4, (2021).https://doi.org/10.1029/2020PA004082
- See more
Loading...
View Options
Get Access
Log in to view the full text
AAAS login provides access to Science for AAAS Members, and access to other journals in the Science family to users who have purchased individual subscriptions.
- Become a AAAS Member
- Activate your AAAS ID
- Purchase Access to Other Journals in the Science Family
- Account Help
Log in via OpenAthens.
Log in via Shibboleth.
More options
Purchase digital access to this article
Download and print this article for your personal scholarly, research, and educational use.
Buy a single issue of Science for just $15 USD.
View options
PDF format
Download this article as a PDF file
Download PDF





