Interpretation of the carbon abundance in Saturn measured by Cassini

Spectral observations of Saturn from the far infrared spectrometer aboard the Cassini spacecraft [Flasar, F.M., et al., 2005. Temperatures, winds, and composition in the Saturnian system. Science 307, 1247-1251] have revealed that the C/H ratio in the planet is in fact about twice higher than previously derived from ground based observations and in agreement with the C/H value derived from Voyager IRIS by Courtin et al. [1984. The composition of Saturn's atmosphere at northern temperate latitudes from Voyager IRIS spectra - NH₃, PH₃, C₂H₂, C₂H₆, CH₃D, CH₄, and the Saturnian D/H isotopic ratio. Astrophys. J. 287, 899-916]. The implications of this measurement are reanalyzed in the present report on the basis that volatiles observed in cometary atmospheres, namely CO₂, CH₄, NH₃ and H₂S may have been trapped as solids in the feeding zone of the planet. CH₄ and H₂S may have been in the form of clathrate hydrates while CO₂ presumably condensed in the cooling solar nebula. Carbon may also have been incorporated in organics. Conditions of temperature and pressure ease the hydratation of NH₃. Such icy grains were included in planetesimals which subsequently collapsed into the hydrogen envelope of the planet, then resulting in C, N and S enrichments with respect to the solar abundance. Our calculations are consistent, within error bars, with observed elemental abundances on Saturn provided that the carbon trapped in planetesimals was mainly in the form of CH₄ clathrate and CO₂ ice (and maybe as organics) while nitrogen was in the form of NH₃ hydrate. Our approach has implications on the possible pattern of noble gases in Saturn, since we predict that contrary to what is observed in Jupiter, Ar and Kr should be in solar abundance while Xe might be strongly oversolar. The only way to verify this scenario is to send a probe making in situ mass spectrometer measurements. Our scenario also predicts that the ¹⁴N/¹⁵N ratio should be somewhat smaller in Saturn than measured in Jupiter by Galileo.     

Numerical simulations of type III planetary migration – II. Inward migration of massive planets

We present a numerical study of rapid, so-called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of inward migration, and study in detail its evolution and physics, concentrating on the structure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using the results of global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. The initial conditions are chosen to satisfy the condition for rapid inward migration. We find that type III migration can be divided into two regimes, fast and slow. The structure of the co-orbital region, mass accumulation rate and migration behaviour differ between these two regimes. All our simulations show a transition from the fast to the slow regime, ending type III migration well before reaching the star. The stopping radius is found to be larger for more massive planets and less massive discs. A sharp density drop is also found to be an efficient stopping mechanism. In the fast migration regime the migration rate and induced eccentricity are lower for less massive discs, but almost do not depend on planet mass. Eccentricity is damped on the migration time-scale.     

Adiabatic chaos in the Prometheus–Pandora system

The chaotic orbital motion of Prometheus and Pandora, the 16th and 17th satellites of Saturn, is studied. Chaos in their orbital motion, as found by Goldreich & Rappaport and Renner & Sicardy, is due to interaction of resonances in the resonance multiplet corresponding to the 121:118 commensurability of the mean motions of the satellites. It is shown rigorously that the system moves in adiabatic regime. The Lyapunov time (the 'time horizon of predictability' of the motion) is calculated analytically and compared to the available numerical–experimental estimates. For this purpose, a method of analytical estimation of the maximum Lyapunov exponent in the perturbed pendulum model of non-linear resonance is applied. The method is based on the separatrix map theory. An analytical estimate of the width of the chaotic layer is made as well, based on the same theory. The ranges of chaotic diffusion in the mean motion are shown to be almost twice as big compared to previous estimates for both satellites.     

Secular frequencies of 3-D exoplanetary systems

Using a 12th order expansion of the perturbative potential in powers of the eccentricities and the inclinations, we study the secular effects of two non-coplanar planets which are not in mean–motion resonance. By means of Lie transformations (which introduce an action–angle formulation of the Hamiltonian), we find the four fundamental frequencies of the 3-D secular three-body problem and compute the long-term time evolutions of the Keplerian elements. To find the relations between these elements, the main combinations of the fundamental frequencies common to these evolutions are identified by frequency analysis. This study is performed for two different reference frames: a general one and the Laplace plane. We underline the known limitations of the linear Laplace–Lagrange theory and point out the great sensitivity of the 3-D secular three-body problem to its initial values. This analytical approach is applied to the exoplanetary system Andromedae in order to search whether the eccentricities evolutions and the apsidal configuration (libration of ) observed in the coplanar case are maintained for increasing initial values of the mutual inclination of the two orbital planes. Anne-Sophie Libert is FNRS Research Fellow.     

The coating hypothesis for ammonia ice particles in Jupiter: Laboratory experiments and optical modeling

We report laboratory experiments and modeling calculations investigating the effect of a hydrocarbon coating on ammonia ice spectral signatures. Observational evidence and thermochemical models indicate an abundance of ammonia ice clouds in Jupiter's atmosphere. However, spectrally identifiable ammonia ice clouds are found covering less than 1% of Jupiter's atmosphere, notably in areas of strong vertical transport, indicating a short lifetime for the signature of ammonia absorption on condensed ammonia particles [Baines, K.H., Carlson, R.W., Kamp, L.W., 2002. Icarus 159, 74-94]. Current literature has suggested coating of ammonia ice particles by a hydrocarbon haze as a possible explanation for this paradox. The work presented here supports the inference of a coating effect that can alter or suppress ammonia absorption features. In the experiments, thin films of ammonia ices are deposited in a cryogenic apparatus, coated with hydrocarbons, and characterized by reflection-absorption infrared spectroscopy. We have observed the effects on the ammonia ice absorption features near 3 and 9 μm with coverage by thin layers of hydrocarbons. Modeling calculations of these multilayer thin films assist in the interpretation of the experimental results and reveal the important role of optical interference in altering the aforementioned ammonia spectral features. Mie and T-matrix scattering calculations demonstrate analogous effects for ammonia ice particles and investigate the relative effects of ammonia ice particle size, shape, and coating layer thickness on the ice particle spectral signatures.     

The formation of the Oort cloud in open cluster environments

We study the influence of an open cluster environment on the formation and current structure of the Oort cloud. To do this, we have run 19 different simulations of the formation of the Oort cloud for 4.5 Gyrs. In each simulation, the Solar System spends its first 100 Myrs in a different open cluster environment before transitioning to its current field environment. We find that, compared to forming in the field environment, the inner Oort cloud is preferentially loaded with comets while the Sun resides in the open cluster and that most of this material remains locked in the interior of the cloud for the next 4.4 Gyrs. In addition, the outer Oort cloud trapping efficiencies we observe in our simulations are lower than previous formation models by about a factor of 2, possibly implying an even more massive early planetesimal disk. Furthermore, some of our simulations reproduce the orbits of observed extended scattered disk objects, which may serve as an observational constraint on the Sun's early environment. Depending on the particular open cluster environment, the properties of the inner Oort cloud and extended scattered disk can vary widely. On the other hand, the outer portions of the Oort cloud in each of our simulations are all similar.     

Core formation in giant gaseous protoplanets

Sedimentation rates of silicate grains in gas giant protoplanets formed by disk instability are calculated for protoplanetary masses between 1 MSaturn to 10 MJupiter. Giant protoplanets with masses of 5 MJupiter or larger are found to be too hot for grain sedimentation to form a silicate core. Smaller protoplanets are cold enough to allow grain settling and core formation. Grain sedimentation and core formation occur in the low mass protoplanets because of their slow contraction rate and low internal temperature. It is predicted that massive giant planets will not have cores, while smaller planets will have small rocky cores whose masses depend on the planetary mass, the amount of solids within the body, and the disk environment. The protoplanets are found to be too hot to allow the existence of icy grains, and therefore the cores are predicted not to contain any ices. It is suggested that the atmospheres of low mass giant planets are depleted in refractory elements compared with the atmospheres of more massive planets. These predictions provide a test of the disk instability model of gas giant planet formation. The core masses of Jupiter and Saturn were found to be ∼0.25 M_⊕ and ∼0.5 M_⊕, respectively. The core masses of Jupiter and Saturn can be substantially larger if planetesimal accretion is included. The final core mass will depend on planetesimal size, the time at which planetesimals are formed, and the size distribution of the material added to the protoplanet. Jupiter's core mass can vary from 2 to 12 M_⊕. Saturn's core mass is found to be ∼8 M_⊕.