对火星轨道变化问题的最后解释

cuses particularly on the inner four terrestrial planets for which the orbital time-scales are much shorter than those of the outer five planets. As we can see clearly from the planar orbital configurations shown in Figs 2 and 3， orbital positions of the terrestrial planets differ little between the initial and final part of each numerical integration， which spans several Gyr. The solid lines denoting the present orbits of the planets lie almost within the swarm of dots even in the final part of integrations (b) and (d). This indicates that throughout the entire integration period the almost regular variations of planetary orbital motion remain nearly the same as they are at present.

Vertical view of the four inner planetary orbits (from the z -axis direction) at the initial and final parts of the integrationsN±1. The axes units are au. The xy -plane is set to the invariant plane of Solar system total angular momentum.(a) The initial part ofN+1 ( t = 0 to 0.0547 × 10 9 yr).(b) The final part ofN+1 ( t = 4.9339 × 10 8 to 4.9886 × 10 9 yr).(c) The initial part of N?1 (t= 0 to ?0.0547 × 109 yr).(d) The final part ofN?1 ( t =?3.9180 × 10 9 to ?3.9727 × 10 9 yr). In each panel， a total of 23 684 points are plotted with an interval of about 2190 yr over 5.47 × 107 yr . Solid lines in each panel denote the present orbits of the four terrestrial planets (taken from DE245).

The variation of eccentricities and orbital inclinations for the inner four planets in the initial and final part of the integration N+1 is shown in Fig. 4. As expected， the character of the variation of planetary orbital elements does not differ significantly between the initial and final part of each integration， at least for Venus， Earth and Mars. The elements of Mercury， especially its eccentricity， seem to change to a significant extent. This is partly because the orbital time-scale of the planet is the shortest of all the planets， which leads to a more rapid orbital evolution than other planets; t