Analytical description of physical librations of saturnian coorbital satellites Janus and Epimetheus

Janus and Epimetheus are famously known for their distinctive horseshoe-shaped orbits resulting from a 1:1 orbital resonance. Every 4 years these two satellites swap their orbits by a few tens of kilometers as a result of their close encounter. Recently Tiscareno et al. (Tiscareno, M.S., Thomas, P.C., Burns, J.A. [2009]. Icarus 204, 254-261) have proposed a model of rotation based on images from the Cassini orbiter. These authors inferred the amplitude of rotational librational motion in longitude at the orbital period by fitting a shape model to Cassini ISS images. By a quasi-periodic approximation of the orbital motion, we describe how the orbital swap impacts the rotation of the satellites. To that purpose, we have developed a formalism based on quasi-periodic series with long- and short-period librations. In this framework, the amplitude of the libration at the orbital period is found proportional to a term accounting for the orbital swap. We checked the analytical quasi-periodic development by performing a numerical simulation and find both results in good agreement. To complete this study, the results obtained for the short-period librations are studied with the help of an adiabatic-like approach.     

Dynamical classification of trans-neptunian objects: Probing their origin, evolution, and interrelation

The orbital structure of trans-neptunian objects (TNOs) in the trans-neptunian belt (Edgeworth-Kuiper belt) and scattered disk provides important clues to understand the origin and evolution of the Solar System. To better characterize these populations, we performed computer simulations of currently observed objects using long-arc orbits and several thousands of clones. Our preliminary analysis identified 622 TNOs, and 65 non-resonant objects whose orbits penetrate that of at least one of the giant planets within 1 Myr (the centaurs). In addition, we identified 196 TNOs locked in resonances with Neptune, which, sorted by distance from the Sun, are 1:1 (Neptune trojans), 5:4, 4:3, 11:8, 3:2, 18:11, 5:3, 12:7, 19:11, 7:4, 9:5, 11:6, 2:1, 9:4, 16:7, 7:3, 12:5, 5:2, 8:3, 3:1, 4:1, 11:2, and 27:4. Kozai resonant TNOs are found inside the 3:2, 5:3, 7:4, and 2:1 resonances. We present detailed general features for the resonant populations (i.e., libration amplitude angles, libration centers, Kozai libration amplitudes, etc.). Taking together the simulations of Lykawka and Mukai [Lykawka, P.S., Mukai, T., 2007. Icarus 186, 331-341], an improved classification scheme is presented revealing five main classes: centaurs, resonant, scattered, detached and classical TNOs. Scattered and detached TNOs (non-resonant) have q (perihelion distance) <37 AU and q>40 AU, respectively. TNOs with 37 AU<q<40 AU occupy an intermediate region where both classes coexist. Thus, there are no clear boundaries between the scattered and detached regions. We also securely identified a total of 9 detached TNOs by using 4-5 Gyr orbital integrations. Classical objects are non-resonant TNOs usually divided into cold and hot populations. Their boundaries are as follows: cold classical TNOs (i?5°) are located at 37 AU<a<40 AU (q>37 AU) and 42 AU<a<47.5 AU (q>38 AU), and hot classical TNOs (i>5°) occupy orbits with 37 AU<a<47.5 AU (q>37 AU). However, a more firm classification is found with i>10° for hot classical TNOs. Lastly, we discuss some implications of our classification scheme comparing all TNOs with our model and other past models.     

Mercury's capture into the 3/2 spin-orbit resonance including the effect of core-mantle friction

The rotation of Mercury is presently captured in a 3/2 spin-orbit resonance with the orbital mean motion. The capture mechanism is well understood as the result of tidal interactions with the Sun combined with planetary perturbations [Goldreich, P., Peale, S., 1966. Astron. J. 71, 425-438; Correia, A.C.M., Laskar, J., 2004. Nature 429, 848-850]. However, it is now almost certain that Mercury has a liquid core [Margot, J.L., Peale, S.J., Jurgens, R.F., Slade, M.A., Holin, I.V., 2007. Science 316, 710-714] which should induce a contribution of viscous friction at the core-mantle boundary to the spin evolution. According to Peale and Boss [Peale, S.J., Boss, A.P., 1977. J. Geophys. Res. 82, 743-749] this last effect greatly increases the chances of capture in all spin-orbit resonances, being 100% for the 2/1 resonance, and thus preventing the planet from evolving to the presently observed configuration. Here we show that for a given resonance, as the chaotic evolution of Mercury's orbit can drive its eccentricity to very low values during the planet's history, any previous capture can be destabilized whenever the eccentricity becomes lower than a critical value. In our numerical integrations of 1000 orbits of Mercury over 4 Gyr, the spin ends 99.8% of the time captured in a spin-orbit resonance, in particular in one of the following three configurations: 5/2 (22%), 2/1 (32%) and 3/2 (26%). Although the present 3/2 spin-orbit resonance is not the most probable outcome, we also show that the capture probability in this resonance can be increased up to 55% or 73%, if the eccentricity of Mercury in the past has descended below the critical values 0.025 or 0.005, respectively.     

Latitudinal librations of Mercury with a fluid core

In the framework of the space missions to Mercury, an accurate model of rotation is needed. Librations around the 3:2 spin-orbit resonance as well as latitudinal librations have to be predicted with the best possible accuracy. In this paper, we use a Hamiltonian analysis and numerical integrations to study the librations of Mercury, both in longitude and latitude. Due to the proximity of the period of the free libration in longitude to the orbital period of Jupiter, the 88-day and 11.86-year contributions dominate Mercury’s libration in longitude (with the Hermean parameters chosen). The amplitude of the libration in latitude is much smaller (under 1 arcsec) and should not be detected by the space missions. Nevertheless, we point out that this amplitude could be much larger (up to several tens of arcsec) if the free period related to the libration in latitude approaches the period of the Jupiter-Saturn Great Inequality (883 years). Given the large uncertainties on the planetary parameters, this new resonant forcing on Mercury’s libration in latitude should be borne in mind.     

Did tidal deformation power the core dynamo of Mars?

We first show that 7 out of the 20 giant impact basins of Mars recently reported by Frey [Frey, H., 2008. Geophys. Res. Lett. 35. L13203] trace a great circle on Mars. The other five basins trace another great circle and still the other three basins trace yet another great circle. The latter great circle is in good agreement with the pre-Tharsis equator of Mars that is estimated from modeling crustal magnetic anomalies [Arkani-Hamed, J., 2001. Geophys. Res. Lett. 28, 3409-3412] and diagonalizing the moment of inertia of Mars after removing the loading effects of Tharsis bulge [Sprenke, K.F., Baker, L.L., Williams, A.F., 2005. Icarus 174, 486-489]. It is shown in this paper that the three great circles were likely the equatorial plane of Mars at certain times and Mars experienced appreciable polar wander. The great circles also indicate that the asteroids that created the basins were satellites of Mars whose orbits decayed in time through spin-orbit coupling with tidally deforming Mars, and eventually impacted on the planet creating the giant basins at around 4 Ga. The orbital dynamics of four largest asteroids show that they could have orbited Mars for several hundred million years if they were retrograde satellites. Continual elliptical straining of otherwise circular fluid streamlines of the liquid core of Mars by tidal deformation could have exerted a strong strain that was large enough to overcome dissipation and excite the elliptical instability inside the core. We investigate the physical properties of the martian core that are required to allow the tidal deformation to power the core dynamo, i.e., the growth time of the elliptical instability to become shorter than the dissipation time. The tidal energy dissipation rate inside Mars caused by even only one of the 4 largest asteroids is found to be over two orders of magnitude greater than the magnetic energy dissipation rate in the core, indicating that if only one of the 4 largest asteroids were orbiting in retrograde sense, it would have likely powered the core dynamo of Mars for several hundred million years.     

Chaos induced by De Haerdtl inequality in the Galilean system

This paper aims at studying the long-term orbital consequences of the perturbations related to De Haerdtl inequality, a current quasi-commensurability between the Galilean satellites of Jupiter Ganymede and Callisto. We used the method of Frequency Map Analysis to detect a chaotic behavior in a 5-bodies system where every inequality has been dropped, except of De Haerdtl one. We also used Frequency Analysis to draw the behavior of the arguments likely to become resonant, in several numerical integrations. We show that De Haerdtl inequality might have induced chaos in the past if Ganymede's and Callisto's eccentricities have been higher than 4×¹⁰⁻³. Moreover, we enlight the influence of Jupiter's obliquity on this chaos. We also enlight some aspects of this chaotic behavior, showing for instance stable chaos and single resonances. The main result of this study is that De Haerdtl inequality should be taken into account in every study of the long term orbital evolution of the Galilean satellites.     

Saturn A ring surface mass densities from spiral density wave dispersion behavior

We have undertaken an analysis of the Voyager photopolarimeter (PPS) stellar occultation data of Saturn's A ring. The Voyager PPS observed the bright star δ Scorpii as it was occulted by Saturn's main rings during the spacecraft flyby of the Saturn system in 1981. The occultation measurement produced a ring profile with radial resolution of approximately 100 m, and radial structure is evident in the profile down to the resolution limit. We have applied an autoregressive technique to the data for estimating the power spectrum as a function of radius, which has allowed us to identify 40 spiral density waves in Saturn's A ring, associated with the strongest torques due to forcing from the moons. The majority of the detected waves are observed to disperse linearly in regions beginning a few kilometers from the resonance location. We have used the dispersion behavior for those waves to calculate local surface mass densities in the vicinity of each wave. We find that the inner three-quarters of the A ring (up to the beginning of the Encke gap) has an average surface mass density of , while the outer region has an average surface mass density of . The two regions have different mean surface mass densities with a significance of approximately 0.999993, as estimated with a T-statistic, which corresponds to about 4.5σ. While the mean optical depth of the A ring increases slightly with increasing distance from Saturn, we find that it is not significantly correlated with the surface mass density; the two quantities having a linear Pearson's correlation coefficient of r_corr≈−0.03. The variation of mass density, independent of PPS optical depth, is consistent with previous conjectures that the particle size distribution and composition are not constant across the entire A ring, particularly in the very outer portion. We estimate the mass of Saturn's A ring from our surface mass density estimates as 4.9×10²¹ gm, or 8.61×10⁻⁹ of the mass of Saturn, roughly equivalent to the mass of a 110-km diameter icy satellite. This mass is almost 25% smaller than estimates from previous studies, but is well within the expected errors of the derived mass densities. We also identified three previously unstudied features which exhibit linear dispersion. The strongest of these features is tentatively identified as the Janus 13:11 density wave. The other two features do not fall near any known satellite resonances and may represent density waves created by previously undetected satellites.     

On fine orbit selection for particular geodetic and oceanographic missions involving passage through resonances

 The identification of mean semi-major axes (suitably defined) for satellite orbits to satisfy a variety of requirements for geodesy, geophysics and oceanography, in terms of repeat orbits (with orbital resonances), is investigated. Various options for the definition of semi-major axis, from the viewpoint of satellite dynamics, are described. Simple simulations of the expected resonant changes in inclination are presented, and tools for the analysis of orbit resonances to extract certain lumped harmonic coefficients of the geopotential (e.g. from the very precise CHAMP orbit) are resurrected. Finally, a preliminary example of the 46th-order resonance analysis possible for CHAMP, based on the mean orbital elements produced by GFZ (GeoForschungs Zentrum) for ephemeris prediction, is presented. Received: 10 July 2001 / Accepted: 17 July 2002 Correspondence to: J. Klokočník at Ondřejov Observatory Acknowledgements. We thank Prof. Dr. Ch. Reigber, Dr. P. Schwintzer, Dr. T. Gruber and Dr. R. K?nig from GFZ Potsdam for various consultations and discussions, and for the CHAMP two-line mean elements. This investigation was performed under the aegis of CEDR (Center for Earth's Dynamics Research, Prague-Ondřejov); it has been supported by project LN00A005 (provided by the Ministry of Education of the Czech Republic) and by grant A 3004 of the Grant Agency of the Academy of Sciences of the Czech Republic.     

An efficient, low-velocity, resonant mechanism for capture of satellites by a protoplanet

Numerical simulations of the gravitational scattering of planetesimals by a protoplanet reveal that a significant fraction of scattered planetesimals can become trapped as so-called quasi-satellites in heliocentric 1:1 co-orbital resonance with the protoplanet. While trapped, these resonant planetesimals can have deep low-velocity encounters with the protoplanet that result in temporary or permanent capture onto highly eccentric prograde or retrograde circumplanetary orbits. The simulations include solar nebula gas drag and use planetesimals with diameters ranging from ∼1 to ∼1000 km. Initial protoplanet eccentricities range from ep=0 to 0.15 and protoplanet masses range from 300 Earth-masses (M_⊕) down to 0.1M_⊕. This mass range effectively covers the final masses of all planets currently thought to be in possession of captured satellites—Jupiter, Saturn, Neptune, Uranus, and Mars. For protoplanets on moderately eccentric orbits (ep?0.1) most simulations show from 5-20% of all scattered planetesimals becoming temporarily trapped in the quasi-satellite co-orbital resonance. Typically, 20-30% of the temporarily trapped quasi-satellites of all sizes came within half the Hill radius of the protoplanet while trapped in the resonance. The efficiency of the resonance trapping combined with the subsequent low-velocity circumplanetary capture suggests that this trapped-to-captured transition may be important not only for the origin of captured satellites but also for continued growth of protoplanets.     

Origin of scattered disk resonant TNOs: Evidence for an ancient excited Kuiper belt of 50 AU radius

Resonance occupation of trans-neptunian objects (TNOs) in the scattered disk (>48 AU) was investigated by integrating the orbits of 85 observed members for 4 Gyr. Twenty seven TNOs were locked in the 9:4, 16:7, 7:3, 12:5, 5:2, 8:3, 3:1, 4:1, 11:2, and 27:4 resonances. We then explored mechanisms for the origin of the resonant structure in the scattered disk, in particular the long-term 9:4, 5:2, and 8:3 resonant TNOs (median 4 Gyr), by performing large scale simulations involving Neptune scattering and planetary migration over an initially excited planetesimals disk (wide range of eccentricities and inclinations). To explain the formation of Gyr-resident populations in such distant resonances, our results suggest the existence of a primordial planetesimal disk of at least 45-50 AU radius that suffered a dynamical perturbation leading to 0.1-0.3 or greater eccentricities and a range of inclinations up to ∼20° during early stages of the Solar System history, before planetary migration.