The outer solar system for 200 million years

A special-purpose computer is used to integrate the orbits of the outer five planets for more than 100 Myr into the future and more than 100 Myr into the past. The strongest features in the Fourier transforms of the orbital elements of the Jovian planets can be identified with the frequencies predicted by linear secular theory. Many of the weaker features in the Fourier spectra are identified as linear combinations of the basic frequencies. Serious differences are noted between the present measurements and the predictions of Bretagnon (1974). The amplitude of the 3.796 Myr period libration of Pluto's longitude of perihelion is modulated with a period of 34 Myr. Very long periods, on the order of 137 Myr, are also seen. The orbit of Pluto is stable for the duration of the integration; the maximum Liapunov characteristic exponent is less than 10 to the -6.8 power/yr.