File:Milankovitch Variations.png
According to Milankovitch Theory, the precession of the equinoxes and the apsides, variations in the tilt of the Earth's axis (obliquity) and changes in the eccentricity of the Earth's orbit are responsible for causing the observed 100 kyr cycle in ice ages by varying the amount of sunlight received by the Earth at different times and locations, particularly high northern latitude summer. These changes in the Earth's orbit are the predictable consequence of interactions between the Earth, its moon, and the other planets.
The orbital data shown are from Quinn et al. (1991). Principal frequencies for each of the three kinds of variations are labeled. The solar forcing curve (aka "insolation") is derived from July 1st sunlight at 65 °N latitude according to Jonathan Levine's insolation calculator [1]. The glacial data are from Lisiecki and Raymo (2005) and gray bars indicate interglacial periods, defined here as deviations in the 5 kyr average of at least 0.8 standard deviations above the mean.