I have copied some of the main graphs to my own website to assure display here. Here's the crux of the paper:
Origin of the 100 kyr Glacial Cycle
Figure 2. Spectral fingerprints in the vicinity of the 100 kyr peak: (a) for data from Site 607; (b) for data of the SPECMAP stack; (c) for a model with linear response to eccentricity, calculated from the results of Quinn et al. (ref 6); (d) for the nonlinear ice-sheet model of Imbrie and Imbrie (ref 22); and (e) for a model with linear response to the inclination of the Earth's orbit (measured with respect to the invariable plane). All calculations are for the period 0-600 ka. The 100 kyr peak in the data in (a) and (b) do not fit the fingerprints from the theories (c) and (d), but are a good match to the prediction from inclination in (e). return to beginning
Far more important to our present analysis, however, is the fact that the predicted 100 kyr "eccentricity line" is actually split into 95 and 125 kyr components, in serious conflict with the single narrow line seen in the climate data. The splitting of this peak into a doublet is well known theoretically (see, e.g., ref 5), but in comparisons with data the two peaks in the eccentricity were merged into a single broad peak by the poor resolution of the Blackman-Tukey algorithm (as was done, for example, in ref 8). The single narrow peak in the climate data was likewise broadened, and it appeared to match the broad eccentricity feature.***
Figure 3. Variations of the inclination vector of the Earth's orbit. The inclination i is the angle between this vector and the vector of the reference frame; Omega is the azimuthal angle = the angle of the ascending node (in astronomical jargon).. In (A), (B), and (C) the measurements are made with respect to the zodiacal (or ecliptic) frame, i.e. the frame of the current orbit of the Earth. In (D), (E), and (F) the motion has been trasformed to the invariable frame, i.e. the frame of the total angular momentum of the solar system. Note that the primary period of oscillation in the zodiacal frame (A) is 70 kyr, but in the invariable plane (D) it is 100 kyr.
