Tidal dissipation,orbital evolution,and the nature of Saturn's inner satellites

Estimates of tidal damping times of the orbital eccentricities of Saturn's inner satellites place constraints on some satellite rigidities and dissipation functions Q. These constraints favor rock-like rather than ice-like properties for Mimas and probably Dione. Photometric and other observational data are consistent with relatively higher densities for these two satellites, but require lower densities for Tethys, Enceladus, and Rhea. This leads to a nonmonotonic density distribution for Saturn's inner satellites, apparently determined by different mass fractions of rocky materials. In spite of the consequences of tidal dissipation for the orbital eccentricity decay and implications for satellite compositions, tidal heating is not an important contributor to the thermal history of any Saturnian satellite.     

Orbital evolution of Uranus’s new outer satellites

Data on three recently discovered satellites of Uranus are used to determine basic evolutional parameters of their orbits: the extreme eccentricities and inclinations, as well as the circulation periods of the pericenter arguments and of the longitudes of the ascending nodes. The evolution is mainly investigated by analytically solving Hill’s double-averaged problem for the Uranus-Sun-satellite system, in which Uranus’s orbital eccentricity e _U and inclination i _U to the ecliptic are assumed to be zero. For the real model of Uranus’s evolving orbit with e _U≠0 and i _U≠0, we refine the evolutional parameters of the satellite orbits by numerically integrating the averaged system. Having analyzed the configuration and dynamics of the orbits of Uranus’s five outer satellites, we have revealed the possibility of their mutual crossings and obtained approximate temporal estimates.     

Galilean satellites: Evolutionary paths in deep resonance

The Laplace resonance among the inner three Galilean satellites (mean motions n₁ − 3n₂ + 2n₃ = 0) has stable configurations in “deep resonance,” i.e., where mean motions taken by pairs are in ratios very close to 2:1. The present satellite configuration, with the resonance variable φ ≡ λ₁ − 3λ₂ + 2λ₃ stable at 180°, is unstable near this exact commensurability. But there is a continuous path of stable conditions branching from φ = 180° to higher and lower values of φ and toward very deep resonance, according to a theory extended to third order in orbital eccentricity. This path provides a track for tidal evolution of the system. Thus, scenarios involving evolution (probably episodic) from deep resonance are viable, and eliminate the requirement by the alternative equilibrium hypothesis for rapid tidal dissipation in Jupiter. Evolution out from deep resonance is consistent with the free eccentricity of Ganymede, the free libration of φ, and observational constraints on Io's secular acceleration. Also, the relatively large forced eccentricities in deep resonance may have controlled geophysical processes in the satellites by much greater tidal heating and global stress than at present.     

Two-micron spectrophotometry of Uranus and its rings

Multiaperture K photometry and 2.0- to 2.5-μm spectrophotometry of Uranus and its ring system are presented. The photometric results are used, together with a previously published measurement, to set limits on the geometric albedos of Uranus and the rings at ～2.2 μm: (0.74 ± 0.02) × 10(su?4) ≤ p_K (Uranus) ≤ (1.5 ± 0.3) × 10^?4, and (2.7 ± 0.6) × 10^?2 ≤ p_K (rings) ≤ (3.4 ± 0.1) × 10^?2. Reflectance spectra of Uranus and Uranus plus rings show features in the planet's spectrum which are attributed to gaseous CH₄ absorption, and a 2.20-μm feature in the combined spectrum which may be due to the rings. This feature is tentatively identified with either the 2.26-μm absorption feature of NH₃ frost, or the 2.2-μm OH band exhibited by certain silicate minerals. The results of JHK photometry of Uranus' satellite, Ariel (U1), indicate that the infrared colors of this object are very similar to those of the satellites U2, U3, and U4.     

Tidally Induced Volcanism

The dissipation of tidal energy causes the ongoing silicate volcanism on Jupiter's satellite, Io, and cryovolcanism almost certainly has resurfaced parts of Saturn's satellite, Enceladus, at various epochs distributed over the latter's history. The maintenance of tidal dissipation in Io and the occurrence of the same on Enceladus depends crucially on the maintenance of the respective orbital eccentricities by the existence of mean motion resonances with nearby satellites. A formation of the resonances among the Galilean satellites by differential expansion of the satellite orbits from tides raised on Jupiter by the satellites means the onset of the volcanism on Io could be relatively recent. If, on the other hand, the resonances formed by differential migration from resonant interactions of the satellites with the disk of gas and particles from which they formed, Io would have been at least intermittently volcanically active throughout its history. Either means of assembling the Galilean satellite resonances lead to the same constraint on the dissipation function of Jupiter Q _J 10⁶, where the currently high heat flux from Io seems to favor episodic heating as Io's eccentricity periodically increases and decreases. Either of the two models might account for sufficient tidal dissipation in the icy satellite Enceladus to cause at least occasional cryovolcanism over much of its history. However, both models are assumption-dependent and not secure, so uncertainty remains on how tidal dissipation resurfaced Enceladus.     

Orbital resonances in the inner neptunian system: I. The 2:1 Proteus-Larissa mean-motion resonance

We investigate the orbital resonant history of Proteus and Larissa, the two largest inner neptunian satellites discovered by Voyager 2. Due to tidal migration, these two satellites probably passed through their 2:1 mean-motion resonance a few hundred million years ago. We explore this resonance passage as a method to excite orbital eccentricities and inclinations, and find interesting constraints on the satellites' mean density () and their tidal dissipation parameters (Qs>10). Through numerical study of this mean-motion resonance passage, we identify a new type of three-body resonance between the satellite pair and Triton. These new resonances occur near the traditional two-body resonances between the small satellites and, surprisingly, are much stronger than their two-body counterparts due to Triton's large mass and orbital inclination. We determine the relevant resonant arguments and derive a mathematical framework for analyzing resonances in this special system.     

Turbulent viscosity and Jupiter's tidal Q

A recent estimate of tidal dissipation by turbulent viscosity in Jupiter's convective interior predicts that the current value of the planet's tidal Q ～ 5 × 10⁶. We point out a fundamental error in this calculation, and show that turbulent dissipation alone implies that at present Q ～ 5 × 10¹³. Our reduced estimat for the rate of tidal dissipation shows conclusively that tidal torques have produced only negligible modifications of the orbits of the Galilean satellites over the age of the solar system.     

The Effect of Tidal Interaction with a Gas Disk on Formation of Terrestrial Planets

We have performed N-body simulation on final accretion stage of terrestrial planets, including the effect of damping of eccentricity and inclination caused by tidal interaction with a remnant gas disk. As a result of runway and oligarchic accretion, about 20 Mars-sized protoplanets would be formed in nearly circular orbits with orbital separation of several to ten Hill radius. The orbits of the protoplanets would be eventually destabilized by long-term mutual gravity and/or secular resonance of giant gaseous planets. The protoplanets would coalesce with each other to form terrestrial planets through the orbital crossing. Previous N-body simulations, however, showed that the final eccentricities of planets are around 0.1, which are about 10 times higher than the present eccentricities of Earth and Venus. The obtained high eccentricities are the remnant of orbital crossing. We included the effect of eccentricity damping caused by gravitational interaction with disk gas as a drag force (“gravitational drag”) and carried out N-body simulation of accretion of protoplanets. We start with 15 protoplanets with 0.2M⊕ and integrate the orbits for 10⁷ years, which is consistent with the observationally inferred disk lifetime (in some runs, we start with 30 protoplanets with 0.1M⊕). In most runs, the damping time scale, which is equivalent to the strength of the drag force, is kept constant throughout each run in order to clarify the effects of the damping. We found that the planets' final mass, spatial distribution, and eccentricities depend on the damping time scale. If the damping time scale for a 0.2M⊕ mass planet at 1 AU is longer than 10⁸ years, planets grow to Earth's size, but the final eccentricities are too high as in gas-free cases. If it is shorter than 10⁶ years, the eccentricities of the protoplanets cannot be pumped up, resulting in not enough orbital crossing to make Earth-sized planets. Small planets with low eccentricities are formed with small orbital separation. On the other hand, if it is between 10⁶ and 10⁸ years, which may correspond to a mostly depleted disk (0.01-0.1% of surface density of the minimum mass model), some protoplanets can grow to about the size of Earth and Venus, and the eccentricities of such surviving planets can be diminished within the disk lifetime. Furthermore, in innermost and outermost regions in the same system, we often find planets with smaller size and larger eccentricities too, which could be analogous to Mars and Mercury. This is partly because the gravitational drag is less effective for smaller mass planets, and partly due to the “edge effect,” which means the innermost and outermost planets tend to remain without collision. We also carried out several runs with time-dependent drag force according to depletion of a gas disk. In these runs, we used exponential decay model with e-folding time of 3×10⁶ years. The orbits of protoplanets are stablized by the eccentricity damping in the early time. When disk surface density decays to ?1% of the minimum mass disk model, the damping force is no longer strong enough to inhibit the increase of the eccentricity by distant perturbations among protoplanets so that the orbital crossing starts. In this disk decay model, a gas disk with 10⁻⁴-10⁻³ times the minimum mass model still remains after the orbital crossing and accretional events, which is enough to damp the eccentricities of the Earth-sized planets to the order of 0.01. Using these results, we discuss a possible scenario for the last stage of terrestrial planet formation.     

Tidal dissipation in the solid cores of the major planets

If the orbital resonances in the Jovian and Saturnian satellite systems are the result of orbital evolution due to tidal dissipation then the present rates of energy dissipation ($\text{E}\text{dot}$) are >2 × 10²⁰ ergs sec^?1 (Jupiter) and ?2 × 10¹⁶ ergs sec^?1 (Saturn). These values of $\text{E}\text{dot}$ can be accounted for if the planets have rocky cores with volumes equal to those suggested by current models of the interiors and if the material of these cores is both solid and imperfectly elastic (Q_e ～ 34). The calculated values of Q_e are not strongly dependent on either the rigidity of the core or the densities of the core and the mantle. Thus, these quantities need not be known precisely. It may be significant that approximately the same value of Q_e is needed for all the major planets (Jupiter, Saturn, and Uranus) even though the values of $\text{E}\text{dot}$ for these planets differ by a factor greater than 10⁴.     

Formation of the regular satellites of giant planets in an extended gaseous nebula II: satellite migration and survival

For a satellite to survive in the disk the time scale of satellite migration must be longer than the time scale for gas dissipation. For large satellites (∼1000 km) migration is dominated by the gas tidal torque. We consider the possibility that the redistribution of gas in the disk due to the tidal torque of a satellite with mass larger than the inviscid critical mass causes the satellite to stall and open a gap (W.R. Ward, 1997, Icarus 26, 261-281). We adapt the inviscid critical mass criterion to include gas drag, and m-dependent nonlocal deposition of angular momentum. We find that such a model holds promise of explaining the survival of satellites in the subnebula, the mass versus distance relationship apparent in the saturnian and uranian satellite systems, the concentration of mass in Titan, and the observation that the satellites of Jupiter get rockier closer to the planet whereas those of Saturn become increasingly icy. It is also possible that either weak turbulence (close to the planet) or gap-opening satellite tidal torque removes gas on a similar time scale (10⁴-10⁵ years) as the orbital decay time of midsized (200-700 km) regular satellites forming in the inner disk (inside the centrifugal radius (I. Mosqueira and P.R. Estrada, 2003, Icarus, this issue)). We argue that Saturn’s satellite system bridges the gap between those of Jupiter and Uranus by combining the formation of a Galilean-sized satellite in a gas optically thick subnebula with a strong temperature gradient, and the formation of smaller satellites, closer to the planet, in a disk with gas optical depth ?1, and a weak temperature gradient.Using an optically thick inner disk (given gaseous opacity), and an extended, quiescent, optically thin outer disk, we show that there are regions of the disk of small net tidal torque (even zero) where satellites (Iapetus-sized or larger) may stall far from the planet. For our model these outer regions of small net tidal torque correspond roughly to the locations of Callisto and Iapetus. Though the precise location depends on the (unknown) size of the transition region between the inner and outer disks, the result that Saturn’s is found much farther out (at ∼3rcS, where rcS is Saturn’s centrifugal radius) than Jupiter’s (at ∼ 2rcJ, where rcJ is Jupiter’s centrifugal radius) is mostly due to Saturn’s less massive outer disk and larger Hill radius. However, despite the large separation between Ganymede and Callisto and Titan and Iapetus, the long formation and migration time scales for Callisto and Iapetus (I. Mosqueira and P.R. Estrada, 2003, Icarus, this issue) makes it possible (depending on the details of the damping of acoustic waves) that the tidal torque of Ganymede and Titan clears the gas disk out to their location, thus stranding Callisto and Iapetus far from the planet. Either way, our model provides an explanation for the presence of regular satellites outside the centrifugal radii of Jupiter and Saturn, and the absence of such a satellite for Uranus.     

Orbital history of the Martian satellites with inferences on their origin

Recent Viking results indicate the Martian satellites are composed of carbonaceous chondritic material, suggesting that Phobos and Deimos were once asteroids captured by Mars. On the other hand, the low eccentricities and inclinations of their orbits on the equator of Mars argue against that hypothesis. This paper presents detailed calculations of the tidal evolution of Phobos and Deimos, considering dissipation in both Mars and its satellites simultaneously and using a new method applicable for any value of the eccentricity. In particular, including precession of the satellites' orbits indicates that they have always remained close to their Laplacian plane, so that the orbital planes of Phobos and Deimos switched from near the Martian orbital plane to the Martian equator once the perturbations due to the planetary oblateness dominated the solar perturbations, as they do presently. The results show that Deimos has been little affected by tides, but several billion (10⁹) years ago, Phobos was in a highly eccentric orbit lying near the common plane of the solar system. This outcome is obtained for very reasonable values of dissipation inside Mars and inside Phobos. Implications for the origin of the Martian satellites are discussed.     

Peculiarities of Uranus’s satellite system

The absence of Uranus’s equatorial satellites in the region of approximately equal influence of its oblateness and solar perturbations is explained in terms of an improved physical model. This model is more complete than the previously studied case of an integrable averaged problem. The model improvement stems from the fact that the inclination of Uranus’s equator to the ecliptic differs by 90° and that the orbital evolution of Uranus due to secular planetary perturbations is taken into account. The lifetime of Uranus’s hypothetical satellites in orbits with semimajor axes 1.3–7 million km can be estimated by numerically integrating the evolution equations to be ～10⁴ yr. This is the time scale on which the evolution of the orbits leads to their intersection with the orbits of inner satellites.     

The discovery of faint irregular satellites of Uranus

We report the discovery of four new uranian irregular satellites in our deep, m_R∼25.4, optical search around that planet. The orbital properties of these satellites are diverse. There is some grouping of inclinations and one of the satellites appears to be inside the Kozai resonant zone of Uranus. Further, we find that the differential size distribution of satellites is rather shallow compared to objects in the asteroid and Kuiper belts, going as ∼r^−2.4. We also report a strong coupling between semi-major axis and orbital eccentricity. We comment on the apparent paradox between the inclination grouping, shallow size distribution, and orbital correlation as they relate to the likelihood of a collisional origin for the uranian irregulars. The currently observed irregulars appear to be consistent with a disruptive formation process and a collisional origin for Uranus' obliquity.     

Secular ring instability in the protoplanetary accretion disk

The basic parameters describing the angular momentum distribution within the Uranus system and of its tidal evolution have been estimated. The nine satellites orbiting under the synchronous zone of Uranus is the maximum number in the solar system and it makes the Uranus system different compared with any other in the Solar system, however the satellites in question are relatively small and their contribution of the tidal dynamics of the system is small compared with that due to UI and UV. The time for existence of the nine satellites as integrated bodies can be estimated as 1.4 × 10⁹ y (UVI) and more. The total tidal decrease in the Uranus angular velocity of rotation is estimated as 7 × 10^–9s^–1.     

Secular resonances and the origin of eccentricities of Mars and the asteroids

A mechanism is treated for the origin of the eccentricities of the asteroids and of Mars: secular resonances associated with the dissipation of a primitive solar nebula. The nebula is modeled as a two-dimensional disk; a closed-form, convergent integral is derived to represent its disturbing function. Dissipation of this nebula gives rise to “excitation waves”, produced by the variable location of the secular resonances, which can excite the eccentricity of Mars, and scatter asteroidal eccentricities through the observed ranges. By requiring that these ranges match the observed values as a functions of semimajor axis, one infers: (a) the primordial eccentricities of Jupiter and Saturn initially had amplitudes different from present-day values, but these amplitudes approached the present values toward the end of nebular dissipation; (b) the nebular dissipation time scale may have been of the order of (few) × 10⁴ years as the dissipation neared completion (but this depends on the validity of linear equations which model the inherently nonlinear asteroidal eccentricity pumping); (c) it is reasonable to propose a common origin for the eccentricies of Mars and the asteroids. A simple extension of the model also accounts for the quasi-Gaussian distribution of the number density of asteroidal eccentricities.     

Editorial

The Galilean satellites Io, Europa, and Ganymede interact through several stable orbital resonances where $\text{λ}{}_1\text{? 2λ}{}_2\text{+}\text{ω}{}_1\text{= 0, λ}{}_1\text{? 2λ}{}_2\text{+}\text{ω}{}_2\text{= 180°, λ}{}_2\text{? 2λ}{}_3\text{+}\text{ω}{}_2\text{= 0}$ and λ₁ ? 3λ₂ + 2λ₃ = 180°, with λ_i being the mean longitude of the ith satellite and $\text{ω}{}_{\mathrm{i}}$ the longitude of the pericenter. The last relation involving all three bodies is known as the Laplace relation. A theory of origin and subsequent evolution of these resonances outlined earlier (C. F. Yoder, 1979b, Nature279, 747–770) is described in detail. From an initially quasi-random distribution of the orbits the resonances are assembled through differential tidal expansion of the orbits. Io is driven out most rapidly and the first two resonance variables above are captured into libration about 0 and 180° respectively with unit probability. The orbits of Io and Europa expand together maintaining the 2:1 orbital commensurability and Europa's mean angular velocity approaches a value which is twice that of Ganymede. The third resonance variable and simultaneously the Laplace angle are captured into libration with probability ～0.9. The tidal dissipation in Io is vital for the rapid damping of the libration amplitudes and for the establishment of a quasi-stationary orbital configuration. Here the eccentricity of Io's orbit is determined by a balance between the effects of tidal dissipation in Io and that in Jupiter, and its measured value leads to the relation $\text{k}{}_1\text{?}{}_1\text{/Q}{}_1\text{≈ 900k}{}_{\mathrm{<mtext>J</mtext>}}\text{/Q}{}_{\mathrm{<mtext>J</mtext>}}$ with the k's being Love numbers, the Q's dissipation factors, and f a factor to account for a molten core in Io. This relation and an upper bound on Q₁ deduced from Io's observed thermal activity establishes the bounds 6 × 10⁴ < Q_J < 2 × 10⁶, where the lower bound follows from the limited expansion of the satellite orbits. The damping time for the Laplace libration and therefore a minimum lifetime of the resonance is 1600 Q_J years. Passage of the system through nearby three-body resonances excites free eccentricities. The remnant free eccentricity of Europa leads to the relation $\text{Q}{}_2\text{/?}{}_2\text{? 2 × 10}{}^{\mathrm{?4}}\text{Q}{}_{\mathrm{<mtext>J</mtext>}}$ for rigidity μ₂ = 5 × 10¹¹ dynes/cm². Probable capture into any of several stable 3:1 two-body resonances implies that the ratio of the orbital mean motions of any adjacent pair of satellites was never this large.A generalized Hamiltonian theory of the resonances in which third-order terms in eccentricity are retained is developed to evaluate the hypothesis that the resonances were of primordial origin. The Laplace relation is unstable for values of Io's eccentricity e₁ > 0.012 showing that the theory which retains only the linear terms in e₁ is not valid for values of e₁ larger than about twice the current value. Processes by which the resonances can be established at the time of satellite formation are undefined, but even if primordial formation is conjectured, the bounds established above for Q_J cannot be relaxed. Electromagnetic torques on Io are also not sufficient to relax the bounds on Q_J. Some ideas on processes for the dissipation of ideal energy in Jupiter yield values of Q_J within the dynamical bounds, but no theory has produced a Q_J small enough to be compatible with the measurements of heat flow from Io given the above relation between Q₁ and Q_J. Tentative observational bounds on the secular acceleration of Io's mean motion are also shown not to be consistent with such low values of Q_J. Io's heat flow may therefore be episodic. Q_J may actually be determined from improved analysis of 300 years of eclipse data.     

On the Inclination Distribution of the Jovian Irregular Satellites

Irregular satellites—moons that occupy large orbits of significant eccentricity e and/or inclination I—circle each of the giant planets. The irregulars often extend close to the orbital stability limit, about 1/3-1/2 of the way to the edge of their planet's Hill sphere. The distant, elongated, and inclined orbits suggest capture, which presumably would give a random distribution of inclinations. Yet, no known irregulars have inclinations (relative to the ecliptic) between 47 and 141°.This paper shows that many high-I orbits are unstable due to secular solar perturbations. High-inclination orbits suffer appreciable periodic changes in eccentricity; large eccentricities can either drive particles with ∼70°<I<110° deep into the realm of the regular satellites (where collisions and scatterings are likely to remove them from planetocentric orbits on a timescale of 10⁷-10⁹ years) or expel them from the Hill sphere of the planet.By carrying out long-term (10⁹ years) orbital integrations for a variety of hypothetical satellites, we demonstrate that solar and planetary perturbations, by causing particles to strike (or to escape) their planet, considerably broaden this zone of avoidance. It grows to at least 55°<I<130° for orbits whose pericenters freely oscillate from 0 to 360°, while particles whose pericenters are locked at ±90° (Kozai mechanism) can remain for longer times.We estimate that the stable phase space (over 10 Myr) for satellites trapped in the Kozai resonance contains ∼10% of all stable orbits, suggesting the possible existence of a family of undiscovered objects at higher inclinations than those currently known.     

A co-accretional model of satellite formation

Numerical calculation of a simple accretion model including the effects of tidal friction indicate that coformation is tenable only if the planet's Q is less than about 10³. The parameter which most strongly affects the final mass ratio of the pair is the time at which the secondary embryo is introduced. Our model yields the proper Moon-Earth mass ratio if the Moon embryo is introduced when the Earth is only about $\text{1}\text{10}$ of its final mass. The lunar orbit remains at about 10 Earth radii throughout most of the growth.This model of satellite formation overcomes two difficulties of the “circumterrestrial cloud” model of Ruskol (1960, 1963, 1972): (1) The difficulty of accumulating a mass as great as the entire Moon before gravitational instability reduces the cloud to a small number of moonlets is removed. (2) The differences between terrestrial and outer planet satellite systems is easily understood in terms of the differences in Q between these planets. The high Q of the outer planets does not allow a satellite embryo to survive a significant portion of the accretion process, thus only small bodies which formed very late in the accumulation of the planet remain as satellites. The low Q of the terrestrial planets allows satellite embryos of these planets to survive during accretion, thus massive satellites such as the Earth's Moon are expected. The present lack of such satellites of the other terrestrial planets may be the result of tidal evolution, either infall following primary despinning (Burns, 1973) or escape due to increase in orbit eccentricity.     

Occurrence of central peak craters on the Uranian satellites: Implications for surface structure and comparison with the Jovian and Saturnian systems

Craters with central peaks occur on the Uranian satellites Ariel, Umbriel, Titania, and Oberon; but do not occur on Miranda. The inelastic surface of Miranda is apparently due to the heavy tectonic reworking of its surface. A theory of expansion/contraction is proposed to explain the tectonic history of Miranda. The existence of central peak craters on the four largest satellites of Uranus implies that they have surface strengths similar to those of the Saturnian satellites and silicate bodies of the inner solar system which all have central peak craters. The absence of central peak craters on Miranda implies that it has an inelastic surface similar to those of the Jovian ice satellites Ganymede and Callisto whose surfaces do not contain central peak craters.     

On the stability of Earth-like planets in multi-planet systems

We present a continuation of our numerical study on planetary systems with similar characteristics to the Solar System. This time we examine the influence of three giant planets on the motion of terrestrial-like planets in the habitable zone (HZ). Using the Jupiter–Saturn–Uranus configuration we create similar fictitious systems by varying Saturn’s semi-major axis from 8 to 11 AU and increasing its mass by factors of 2–30. The analysis of the different systems shows the following interesting results: (i) Using the masses of the Solar System for the three giant planets, our study indicates a maximum eccentricity (max-e) of nearly 0.3 for a test-planet placed at the position of Venus. Such a high eccentricity was already found in our previous study of Jupiter–Saturn systems. Perturbations associated with the secular frequency g ₅ are again responsible for this high eccentricity. (ii) An increase of the Saturn-mass causes stronger perturbations around the position of the Earth and in the outer HZ. The latter is certainly due to gravitational interaction between Saturn and Uranus. (iii) The Saturn-mass increased by a factor 5 or higher indicates high eccentricities for a test-planet placed at the position of Mars. So that a crossing of the Earth’ orbit might occur in some cases. Furthermore, we present the maximum eccentricity of a test-planet placed in the Earth’ orbit for all positions (from 8 to 11 AU) and masses (increased up to a factor of 30) of Saturn. It can be seen that already a double-mass Saturn moving in its actual orbit causes an increase of the eccentricity up to 0.2 of a test-planet placed at Earth’s position. A more massive Saturn orbiting the Sun outside the 5:2 mean motion resonance (a _S  ≥9.7 AU) increases the eccentricity of a test-planet up to 0.4.