Milankovitch cycles describes the collective effects of changes in the Earth's movements upon its climate, named after Serbian geophysicist and astronomer Milutin Milanković, who in the 1920s had theorized that variations in eccentricity, axial tilt, and precession of the Earth's orbit strongly influenced climatic patterns on Earth through orbital forcing.
The Earth's axis completes one full cycle of precession approximately every 26,000 years. At the same time, the elliptical orbit rotates more slowly. The combined effect of the two precessions leads to a 21,000-year period between the astronomical seasons and the orbit. In addition, the angle between Earth's rotational axis and the normal to the plane of its orbit (obliquity) oscillates between 22.1 and 24.5 degrees on a 41,000-year cycle. It is currently 23.44 degrees and decreasing at a rate of about 0.013° per century.
Similar astronomical theories had been advanced in the 19th century by Joseph Adhemar, James Croll and others, but verification was difficult due to the absence of reliably dated evidence and doubts as to exactly which periods were important. Not until the advent of deep-ocean cores and a seminal paper by Hays, Imbrie, and Shackleton, "Variations in the Earth's Orbit: Pacemaker of the Ice Ages", in Science (1976) did the theory attain its present state.
As the Earth rotates around its axis and revolves around the Sun, several quasi-periodic variations occur due to gravitational interactions. Although the curves have a large number of sinusoidal components, a few components are dominant. Milankovitch studied changes in the orbital eccentricity, obliquity, and precession of Earth's movements. Such changes in movement and orientation alter the amount and location of solar radiation reaching the Earth. This is known as solar forcing (an example of radiative forcing). Changes near the north polar area, about 65 degrees North, are considered important due to the great amount of land. Land masses respond to temperature change more quickly than oceans, which have a higher effective heat capacity, because of the mixing of surface and deep water and the fact that the specific heat of solids is generally lower than that of water.
Orbital shape (eccentricity)
The Earth's orbit is an ellipse. The eccentricity is a measure of the departure of this ellipse from circularity. The shape of the Earth's orbit varies in time between nearly circular (with the lowest eccentricity of 0.000055) and mildly elliptical (highest eccentricity of 0.0679) with the mean eccentricity of 0.0019 as geometric or logarithmic mean. The major component of these variations occurs on a period of 413,000 years (eccentricity variation of ±0.012). A number of other terms vary between components 95,000 and 125,000 years (with a beat period 400,000 years), and loosely combine into a 100,000-year cycle (variation of −0.03 to +0.02). The present eccentricity is 0.017 and decreasing.
If the Earth were the only planet orbiting our Sun, the eccentricity of its orbit would not perceptibly vary even over a period of a million years. The Earth's eccentricity varies primarily due to interactions with the gravitational fields of Jupiter and Saturn. As the eccentricity of the orbit evolves, the semi-major axis of the orbital ellipse remains unchanged. From the perspective of perturbation theory used in celestial mechanics to compute the evolution of the orbit, the semi-major axis is an adiabatic invariant. According to Kepler's third law the period of the orbit is determined by the semi-major axis. It follows that the Earth's orbital period, the length of a sidereal year, also remains unchanged as the orbit evolves.
Orbital shape and Temperature
As the semi-minor axis is decreased with the eccentricity increase, the seasonal changes increase. But the mean solar irradiation for the planet changes only slightly for small eccentricity, due to Kepler's second law. (The intensity of seasons is not primarily governed by Earth's distance from the sun. See Season.)
The relative increase in solar irradiation at closest approach to the Sun (perihelion) compared to the irradiation at the furthest distance (aphelion) is slightly larger than four times the eccentricity. For the current orbital eccentricity this amounts to a variation in incoming solar radiation of about 6.8%, while the current difference between perihelion and aphelion is only 3.4% (5.1 million km). Perihelion presently occurs around January 3, while aphelion is around July 4. When the orbit is at its most elliptical, the amount of solar radiation at perihelion will be about 23% more than at aphelion.
The higher eccentricity also causes extra behavior due to precession and axial tilt. The true global summer does not appear to be the warmer part of the year in the North (until an eon from now when balanced). Seasons always arrive early (see Mean tropical year), but not the same for differing seasons, since elliptical motion (speed, angle) with respect to the sun differs per season arrival.
|Year||Northern Hemisphere||Southern Hemisphere||Date: GMT||Season duration|
|2005||Winter solstice||Summer solstice||21 December 2005 18:35||88.99 days|
|2006||Spring equinox||Autumn equinox||20 March 2006 18:26||92.75 days|
|2006||Summer solstice||Winter solstice||21 June 2006 12:26||93.65 days|
|2006||Autumn equinox||Spring equinox||23 September 2006 4:03||89.85 days|
|2006||Winter solstice||Summer solstice||22 December 2006 0:22||88.99 days|
|2007||Spring equinox||Autumn equinox||21 March 2007 0:07||92.75 days|
|2007||Summer solstice||Winter solstice||21 June 2007 18:06||93.66 days|
|2007||Autumn equinox||Spring equinox||23 September 2007 9:51||89.85 days|
|2007||Winter solstice||Summer solstice||22 December 2007 06:08|
Orbital mechanics requires that the length of the seasons be proportional to the areas of the seasonal quadrants, so when the eccentricity is extreme, the Earth's orbital motion becomes more nonuniform and the lengths of the seasons change. When autumn and winter occur at closest approach, as is the case currently in the northern hemisphere, the earth is moving at its maximum velocity and therefore autumn and winter are slightly shorter than spring and summer. Thus, summer in the northern hemisphere is 4.66 days longer than winter and spring is 2.9 days longer than autumn. But as the orientation of Earth's orbit changes relative to the Vernal Equinox due to apsidal precession, the way the length of the seasons are altered by the nonuniform motion changes, since different sections of the orbit are involved. When the Earth's apsides are aligned with the equinoxes the length of Spring and Summer (together) equals that of Autumn and Winter. When they are aligned with the solstices either Spring and Summer or Autumn and Winter will be at its longest. Increasing the eccentricity lengthens the time spent near aphelion and shortens the time near perihelion.
Changes to the eccentricity do not by themselves change the length of the anomalistic year or the Earth's mean motion along its orbit since they are both functions of the semi-major axis.
Axial tilt (obliquity)
The angle of the Earth's axial tilt (obliquity of the ecliptic) varies with respect to the plane of the Earth's orbit. These slow 2.4° obliquity variations are roughly periodic, taking approximately 41,000 years to shift between a tilt of 22.1° and 24.5° and back again. When the obliquity increases, the amplitude of the seasonal cycle in insolation increases, with summers in both hemispheres receiving more radiative flux from the Sun, and winters less. Conversely, when the obliquity decreases, summers receive less insolation and winters more.
But these changes of opposite sign in summer and winter are not of the same magnitude everywhere on the Earth's surface. At high latitude the annual mean insolation increases with increasing obliquity, while lower latitudes experience a reduction in insolation. Cooler summers are suspected of encouraging the onset of an ice age by melting less of the previous winter's precipitation. Because most of the planet's snow and ice lies at high latitude, it can be argued that lower obliquity favors ice ages for two reasons: the reduction in overall summer insolation and the additional reduction in mean insolation at high latitude.
Scientists using computer models to study more extreme tilts than those that actually occur have concluded that climate extremes at high obliquity would be particularly threatening to advanced forms of life that presently exist on Earth. They noted that high obliquity would not likely sterilize a planet completely, but would make it harder for fragile, warm-blooded land-based life to thrive as it does today.
Currently the Earth is tilted at 23.44 degrees from its orbital plane, roughly halfway between its extreme values. The tilt is in the decreasing phase of its cycle, and will reach its minimum value around the year 11,800 CE ; the last maximum was reached in 8,700 BCE. This trend in forcing, by itself, tends to make winters warmer and summers colder (i.e. milder seasons), as well as cause an overall cooling trend.
Precession is the trend in the direction of the Earth's axis of rotation relative to the fixed stars, with a period of roughly 26,000 years. This gyroscopic motion is due to the tidal forces exerted by the Sun and the Moon on the solid Earth; both contribute roughly equally to this effect.
When the axis points toward the Sun in perihelion (i.e. the north pole is pointed towards the Sun), the northern hemisphere has a greater difference between the seasons while the southern hemisphere has milder seasons. When the axis points away from the Sun in perihelion (i.e. the south pole is pointed towards the Sun), the southern hemisphere has a greater difference between the seasons while the northern hemisphere has milder seasons. The hemisphere that is in summer at perihelion receives much of the corresponding increase in solar radiation, but that same hemisphere in winter at aphelion has a colder winter. The other hemisphere will have a relatively warmer winter and cooler summer.
When the Earth's axis is aligned such that aphelion and perihelion occur near the equinoxes, the northern and southern hemispheres will have similar contrasts in the seasons.
At present, perihelion occurs during the southern hemisphere's summer, and aphelion is reached during the southern winter. Thus the southern hemisphere seasons are somewhat more extreme than the northern hemisphere seasons, when other factors are equal.
In addition, the orbital ellipse itself precesses in space, in an irregular fashion, completing a full cycle every 112,000 years relative to the fixed stars. This happens primarily as a result of interactions with Jupiter and Saturn. Smaller contributions are also made by the sun's oblateness and by the effects of general relativity that are well known for Mercury.
The apsidal precession combined with the axial precession means that it takes between 20,800 and 29,000 years for the perihelion to return to the position in the year. On average this shortens the period of the precession of the equinoxes with respect to the perihelion from 25,771.5 to 23,000 years. Apsidal precession occurs in the plane of the Ecliptic and alters the orientation of the Earth's orbit relative to the Ecliptic. In combination with changes to the eccentricity it alters the length of the seasons.
The inclination of Earth's orbit drifts up and down relative to its present orbit. Milankovitch did not study this three-dimensional movement. This movement is known as "precession of the ecliptic" or "planetary precession".
More recent researchers noted this drift and that the orbit also moves relative to the orbits of the other planets. The invariable plane, the plane that represents the angular momentum of the Solar System, is approximately the orbital plane of Jupiter. The inclination of Earth's orbit drifts up and down relative to its present orbit with a cycle having a period of about 70,000 years. The inclination of the Earth's orbit has a 100,000-year cycle relative to the invariable plane. This is very similar to the 100,000-year eccentricity period. This 100,000-year cycle closely matches the 100,000-year pattern of glacial events.
It has been proposed that a disk of dust and other debris exists in the invariable plane, and this affects the Earth's climate through several possible means. The Earth presently moves through this plane around January 9 and July 9, when there is an increase in radar-detected meteors and meteor-related noctilucent clouds.
A study of the chronology of Antarctic ice cores using oxygen-nitrogen ratios in air bubbles trapped in the ice, which appear to respond directly to the local insolation, concluded that the climatic response documented in the ice cores was driven by northern hemisphere insolation as proposed by the Milankovitch hypothesis (Kawamura et al., Nature, 23 August 2007, vol 448, pp 912–917). This is an additional validation of the Milankovitch hypothesis by a relatively novel method.
Because the observed periodicities of climate fit so well with the orbital periods, the orbital theory has overwhelming support. Nonetheless, there are several difficulties in reconciling theory with observations.
The 100,000-year problem is that the eccentricity variations have a significantly smaller impact on solar forcing than precession or obliquity – according to theory – and hence might be expected to produce the weakest effects. However, the greatest observed response in regard to the ice ages is at the 100,000-year timescale, even though the theoretical forcing is smaller at this scale. During the last 1 million years, the strongest climate signal is the 100,000-year cycle. In addition, despite the relatively great 100,000-year cycle, some have argued that the length of the climate record is insufficient to establish a statistically significant relationship between climate and eccentricity variations. Various explanations for this discrepancy have been proposed, including frequency modulation or various feedbacks (from carbon dioxide, cosmic rays, or from ice sheet dynamics). Some models can reproduce the 100,000-year cycles as a result of non-linear interactions between small changes in the Earth's orbit and internal oscillations of the climate system.
Stage 5 problem
The stage 5 problem refers to the timing of the penultimate interglacial (in marine isotopic stage 5) that appears to have begun ten thousand years in advance of the solar forcing hypothesized to have caused it (also known as the causality problem, because the putative effect precedes the cause).
Effect exceeds cause
The effects of these variations are primarily believed to be due to variations in the intensity of solar radiation upon various parts of the globe. Observations show climate behavior is much more intense than the calculated variations. Various internal characteristics of climate systems are believed to be sensitive to the insolation changes, causing amplification (positive feedback) and damping responses (negative feedback).
The unsplit peak problem
The unsplit peak problem refers to the fact that eccentricity has cleanly resolved variations at both the 95 and 125 ka periods. A sufficiently long, well-dated record of climate change should be able to resolve both frequencies. However, some researchers interpret climate records of the last million years as showing only a single spectral peak at 100 ka periodicity
The transition problem
The transition problem refers to the switch in the frequency of climate variations 1 million years ago. From 1–3 million years, climate had a dominant mode matching the 41 ka cycle in obliquity. After 1 million years ago, this switched to a 100 ka variation matching eccentricity, for which no reason has been established.
Identifying dominant factor
Milankovitch believed that decreased summer insolation in northern high latitudes was the dominant factor leading to glaciation, which led him to (incorrectly) deduce an approximate 41 ka period for ice ages. Subsequent research has shown that ice age cycles of the Quaternary glaciation over the last million years have been at a 100,000-year period, leading to identification of the 100 ka eccentricity cycle as more important, although the exact mechanism remains obscure.
Present and future conditions
As mentioned above, at present, perihelion occurs during the southern hemisphere's summer and aphelion during the southern winter. Thus the southern hemisphere seasons should tend to be somewhat more extreme than the northern hemisphere seasons. The relatively low eccentricity of the present orbit results in a 6.8% difference in the amount of solar radiation during summer in the two hemispheres.
Since orbital variations are predictable, if one has a model that relates orbital variations to climate, it is possible to run such a model forward to "predict" future climate. Two caveats, however, are necessary: that anthropogenic effects may modify or even overwhelm orbital effects; and that the mechanism by which orbital forcing influences climate is not well understood. In the most prominent anthropogenic example, orbital forcing from the Milankovitch cycles has been in a cooling phase for millennia, but that cooling trend was reversed in the 20th and 21st centuries due to warming caused by increased anthropogenic greenhouse gas emissions.
The amount of solar radiation (insolation) in the Northern Hemisphere at 65° N seems to be related to occurrence of an ice age. Astronomical calculations show that 65° N summer insolation should increase gradually over the next 25,000 years. A regime of eccentricity lower than the current value will last for about the next 100,000 years. Changes in northern hemisphere summer insolation will be dominated by changes in obliquity ε. No declines in 65° N summer insolation, sufficient to cause a glacial period, are expected in the next 50,000 years.
An often-cited 1980 study by Imbrie and Imbrie determined that, "Ignoring anthropogenic and other possible sources of variation acting at frequencies higher than one cycle per 19,000 years, this model predicts that the long-term cooling trend that began some 6,000 years ago will continue for the next 23,000 years."
Observed on other planets
Other planets in the Solar System have been discovered to have Milankovitch cycles. Mostly these cycles are not as intense or complex as the Earth's cycles, but they do have a global geological impact with respect to the movement of mobile solids such as water or nitrogen ices or hydrocarbon lakes. The known affected planets are:
- Karney, Kevin. "Variation in the Equation of Time" (PDF).
- Hays, J. D.; Imbrie, J.; Shackleton, N. J. (1976). "Variations in the Earth's Orbit: Pacemaker of the Ice Ages". Science. 194 (4270): 1121–1132. doi:10.1126/science.194.4270.1121. PMID 17790893.
- Girkin, Amy Negich (2005). A Computational Study on the Evolution of the Dynamics of the Obliquity of the Earth (PDF) (Master of Science thesis). Miami University.
- Laskar, J; Fienga, A.; Gastineau, M.; Manche, H (2011). "La2010: A New Orbital Solution for the Long-term Motion of the Earth." (PDF). Astronomy & Astrophysics. 532 (A889): A89. Bibcode:2011A&A...532A..89L. doi:10.1051/0004-6361/201116836.
- Berger A.; Loutre M.F.; Mélice J.L. (2006). "Equatorial insolation: from precession harmonics to eccentricity frequencies" (PDF). Clim. Past Discuss. 2 (4): 519–533. doi:10.5194/cpd-2-519-2006.
- Data from United States Naval Observatory
- Williams, D.M., Pollard, P. (2002). "Earth-like worlds on eccentric orbits: excursions beyond the habitable zone" (PDF). Inter. J. Astrobio. 1: 21–9. doi:10.1017/s1473550402001064.
- van den Hewel, E. P. J. (1966). "On the Precession as a Cause of Pleistocene Variations of the Atlantic Ocean Water Temperatures". Geophysical Journal International. 11: 323–336.
- Richard A Muller, Gordon J. F. MacDonald (1997). "Glacial Cycles and Astronomical Forcing". Science. 277 (5323): 215–8. Bibcode:1997Sci...277..215M. doi:10.1126/science.277.5323.215.
- "Origin of the 100 kyr Glacial Cycle: eccentricity or orbital inclination?". Richard A Muller. Retrieved March 2, 2005.
- Milankovitch, Milutin (1998) . Canon of Insolation and the Ice Age Problem. Belgrade: Zavod za Udz̆benike i Nastavna Sredstva. ISBN 86-17-06619-9.; see also "Astronomical Theory of Climate Change".
- Wunsch, Carl (2004). "Quantitative estimate of the Milankovitch-forced contribution to observed Quaternary climate change". Quaternary Science Reviews. 23 (9–10): 1001–12. Bibcode:2004QSRv...23.1001W. doi:10.1016/j.quascirev.2004.02.014.
- Rial, J.A. (October 2003), "Earth's orbital Eccentricity and the rhythm of the Pleistocene ice ages: the concealed pacemaker" (PDF), Global and Planetary Change, 41 (2): 81–93, doi:10.1016/j.gloplacha.2003.10.003, archived from the original (PDF) on 2011-07-20
- Ghil, Michael (1994). "Cryothermodynamics: the chaotic dynamics of paleoclimate". Physica D. 77 (1–3): 130–159. Bibcode:1994PhyD...77..130G. doi:10.1016/0167-2789(94)90131-7.
- Gildor H, Tziperman E (2000). "Sea ice as the glacial cycles' climate switch: Role of seasonal and orbital forcing". Paleoceanography. 15 (6): 605–615. Bibcode:2000PalOc..15..605G. doi:10.1029/1999PA000461.
- Zachos JC, Shackleton NJ, Revenaugh JS, Pälike H, Flower BP (April 2001). "Climate response to orbital forcing across the Oligocene-Miocene boundary". Science. 292 (5515): 27–48. Bibcode:2001Sci...292..274Z. doi:10.1126/science.1058288. PMID 11303100.
- "Nonlinear coupling between 100 ka periodicity of the paleoclimate records in loess and periodicities of precession and semi-precession" (PDF). ProQuest.
- Imbrie and Imbrie; Ice Ages, solving the mystery, p 158
- Imbrie, Hays, Shackleton Science 1976
- Shackleton, Berger, Peltier An alternative astronomical calibration of the lower Pleistocene timescale based on ODP Site 677 Trans Roy Soc Edinb 1990
- Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume Ayako Abe-Ouchi et al Nature 500 2013
- F. Varadi; B. Runnegar; M. Ghil (2003). "Successive Refinements in Long-Term Integrations of Planetary Orbits" (PDF). The Astrophysical Journal. 592: 620–630. Bibcode:2003ApJ...592..620V. doi:10.1086/375560.
- Kaufman, D. S.; Schneider, D. P.; McKay, N. P.; Ammann, C. M.; Bradley, R. S.; Briffa, K. R.; Miller, G. H.; Otto-Bliesner, B. L.; Overpeck, J. T.; Vinther, B. M.; Abbott, M.; Axford, M.; Bird, Y.; Birks, B.; Bjune, H. J. B.; Briner, A. E.; Cook, J.; Chipman, T.; Francus, M.; Gajewski, P.; Geirsdottir, K.; Hu, A.; Kutchko, F. S.; Lamoureux, B.; Loso, S.; MacDonald, M.; Peros, G.; Porinchu, M.; Schiff, D.; Seppa, C.; Seppa, H.; Arctic Lakes 2k Project Members (2009). "Recent Warming Reverses Long-Term Arctic Cooling". Science. 325 (5945): 1236–1239. doi:10.1126/science.1173983. PMID 19729653.
"Arctic Warming Overtakes 2,000 Years of Natural Cooling". UCAR. September 3, 2009. Retrieved 19 May 2011.
Bello, David (September 4, 2009). "Global Warming Reverses Long-Term Arctic Cooling". Scientific American. Retrieved 19 May 2011.
- "NOAA Paleoclimatology Program – Orbital Variations and Milankovitch Theory".
- J Imbrie; J Z Imbrie (1980). "Modeling the Climatic Response to Orbital Variations". Science. 207 (4434): 943–953. Bibcode:1980Sci...207..943I. doi:10.1126/science.207.4434.943. PMID 17830447.
- Berger A, Loutre MF (2002). "Climate: An exceptionally long interglacial ahead?". Science. 297 (5585): 1287–8. doi:10.1126/science.1076120. PMID 12193773.
- A. Ganopolski, R. Winkelmann & H. J. Schellnhuber (2016). "Critical insolation–CO2 relation for diagnosing past and future glacial inception". Nature. 529: 200–203. doi:10.1038/nature16494.
- Schorghofer, Norbert (23 September 2008). "Temperature response of Mars to Milankovitch cycles" (PDF). Geophysical Research Letters. Ifa.hawaii.edu. 35.
- Schorghofer, Norbert (2008). "Temperature response of Mars to Milankovitch cycles". Geophysical Research Letters. 35 (18): L18201. Bibcode:2008GeoRL..3518201S. doi:10.1029/2008GL034954.
- "3.5 Modeling Milankovitch cycles on Mars (2010 – 90; Annual Symp Planet Atmos)". Confex.
- "Hydrocarbon lakes on Titan – Alex Hayes (SETI Talks)". YouTube.
- Nicholos Wethington (30 November 2009). "Lake Asymmetry on Titan Explained".
- "Sun Blamed for Warming of Earth and Other Worlds". LiveScience.com.
- "Matt McIrvin's Steam-Operated World of Yesteryear – Warming on other planets".
- The oldest reference for Milankovitch cycles is: M. Milankovitch, Mathematische Klimalehre und Astronomische Theorie der Klimaschwankungen, Handbuch der Klimatologie, Band I, Teil A,Berlin, Verlag von Gebrüder Borntraeger, 1030.
- Roe G (2006). "In defense of Milankovitch". Geophysical Research Letters. 33 (24): L24703. Bibcode:2006GeoRL..3324703R. doi:10.1029/2006GL027817. This shows that Milankovitch theory fits the data extremely well, over the past million years, provided that we consider derivatives.
- Kaufmann R. K.; Juselius K. (2016), "Testing competing forms of the Milankovitch hypothesis", Paleoceanography, 31, doi:10.1002/2014PA002767.
- Edvardsson S, Karlsson KG, Engholm M (2002). "Accurate spin axes and solar system dynamics: Climatic variations for the Earth and Mars". Astronomy and Astrophysics. 384 (2): 689–701. Bibcode:2002A&A...384..689E. doi:10.1051/0004-6361:20020029. This is the first work that investigated the derivative of the ice volume in relation to insolation (page 698).
- Zachos J, Pagani M, Sloan L, Thomas E, Billups K (2001). "Trends, Rhythms, and Aberrations in Global Climate 65 Ma to Present". Science. 292 (5517): 686–693. Bibcode:2001Sci...292..686Z. doi:10.1126/science.1059412. PMID 11326091.
This review article discusses cycles and great-scale changes in the global climate during the Cenozoic Era.
- Pälike, H., R.D. Norris, J.O. Herrle, P.A. Wilson, H.K. Coxall, C.H. Lear, N.J. Shackleton, A.K. Tripati, and B.S. Wade (2006). "The Heartbeat of the Oligocene Climate System". Science. 314 (5807): 1894–1898. doi:10.1126/science.1133822.
A 13-million-year continuous record of Oligocene climate from the equatorial Pacific reveals a pronounced “heartbeat” in the global carbon cycle and periodicity of glaciations.
|Wikimedia Commons has media related to Milankovitch cycles.|
|The Wikibook Historical Geology has a page on the topic of: Milankovitch cycles|
- Ice Age – Milankovitch Cycles – National Geographic Channel
- The Milankovitch band, Internet Archive of American Geophysical Union lecture
- Some history of the adoption of the Milankovitch hypothesis (and an alternative)
- More detail on orbital obliquity also matching climate patterns
- "Milutin Milankovitch". On the Shoulders of Giants. Retrieved January 15, 2010.
- The Seasons
- "Climate Forcing Data". NOAA. National Centers for Environmental Information.
Text: includes (calculated) data on orbital variations over the last 50 million years and for the coming 20 million years.
- The Orbital Simulations; by Varadi, Ghil and Runnegar (2003) provide another, slightly different series for orbital eccentricity, and also a series for orbital inclination
- ABC: Earth wobbles linked to extinctions
- Milankovitch Cycles & Ice Ages