Origin and early evolution of the planetary system

It is shown by linear stability analysis that a preplanetary (presatellite) disk of dust and gas with Keplerian velocity field can become unstable due to the collective self-gravity of the disk. The radial distribution of rings, which may result from this instability, is derived. These rings later on can be the formation sites for planets around the Sun and for satellites around the planets. The derived orbits are shown to be in good agreement with that of the planets and the satellites (of Jupiter, Saturn, and Uranus). Predictions and conclusions seem to be possible for the existence of three yet unknown Uranian satellites, the origin of the early Moon and the possible radial extension of the planetary system.     

Models of Jupiter and Saturn after Galileo mission

New models of Jupiter are based on observational data provided by the Galileo spaceprobe, which considerably improved previously existing estimates of the helium abundance in the atmosphere of Jupiter. These data yield for Jupiter’s atmosphere 20% of the solar oxygen abundance and do not agree with the results of the analysis of the collision of comet Shoemaker-Levy 9 with Jupiter (10 times the solar value). Therefore, both the models of Jupiter with water-depleted and water-enriched atmosphere are considered. By analogy with Jupiter, trial models of Saturn with a water-depleted external envelope are also developed. The molecular-metallic phase transition pressure of hydrogen P_m was taken to be 1.5, 2 and 3 Mbar. Since Saturn’s internal molecular envelope is noticeably enriched in the IR-component (its weight concentration, 0.25–0.30, being by a factor of 3–4 higher than in Jupiter), the phase transition pressure in Saturn can be lower than in Jupiter. In the constructed models, the IR-core masses are 3–3.5 M_⊕ for Jupiter and 3–5.5 M_⊕ for Saturn. Jupiter’s and Saturn’s IR-cores can be considered embryos onto which the accretion of the gas occurred during the formation of the planets. The mass of the hydrogen–helium component dispersed in the zone of planetary formation constitutes ≈2–5 planetary masses for Jupiter and ≈11–14 planetary masses for Saturn.     

Formation of terrestrial planets in a dissipating gas disk with Jupiter and Saturn

We have performed N-body simulations on final accretion stage of terrestrial planets, including the eccentricity and inclination damping effect due to tidal interaction with a gas disk. We investigated the dependence on a depletion time scale of the disk, and the effect of secular perturbations by Jupiter and Saturn. In the final stage, terrestrial planets are formed through coagulation of protoplanets of about the size of Mars. They would collide and grow in a decaying gas disk. Kominami and Ida [Icarus 157 (2002) 43-56] showed that it is plausible that Earth-sized, low-eccentricity planets are formed in a mostly depleted gas disk. In this paper, we investigate the formation of planets in a decaying gas disk with various depletion time scales, assuming disk surface density of gas component decays exponentially with time scale of τ_gas. Fifteen protoplanets with are initially distributed in the terrestrial planet regions. We found that Earth-sized planets with low eccentricities are formed, independent of initial gas surface density, when the condition (τ_cross+τ_growth)/2?τ_gas?τ_cross is satisfied, where τ_cross is the time scale for initial protoplanets to start orbit crossing in a gas-free case and τ_growth is the time scale for Earth-sized planets to accrete during the orbit crossing stage. In the cases satisfying the above condition, the final masses and eccentricities of the largest planets are consistent with those of Earth and Venus. However, four or five protoplanets with the initial mass remain. In the final stage of terrestrial planetary formation, it is likely that Jupiter and Saturn have already been formed. When Jupiter and Saturn are included, their secular perturbations pump up eccentricities of protoplanets and tend to reduce the number of final planets in the terrestrial planet regions. However, we found that the reduction is not significant. The perturbations also shorten τ_cross. If the eccentricities of Jupiter and Saturn are comparable to or larger than present values (∼0.05), τ_cross become too short to satisfy the above condition. As a result, eccentricities of the planets cannot be damped to the observed value of Earth and Venus. Hence, for the formation of terrestrial planets, it is preferable that the secular perturbations from Jupiter and Saturn do not have significant effect upon the evolution. Such situation may be reproduced by Jupiter and Saturn not being fully grown, or their eccentricities being smaller than the present values during the terrestrial planets' formation. However, in such cases, we need some other mechanism to eliminate the problem that numerous Mars-sized planets remain uncollided.     

Orbital Evolution of the Sun–Jupiter–Saturn–Uranus–Neptune Four-Planet System on Long-Time Scales

The four-planet problem is solved by constructing an averaged semi-analytical theory of secondorder motion by planetary masses. A discussion is given of the results obtained by numerical integration of the averaged equations of motion for the Sun–Jupiter–Saturn–Uranus–Neptune system over a time interval of 10 Gyr. The integration is based on high-order Runge–Kutta and Everhart methods. The motion of the planets is almost periodic in nature. The eccentricities and inclinations of the planetary orbits remain small. Short-period perturbations remain small over the entire interval of integration. Conclusions are drawn about the resonant properties of the motion. Estimates are given for the accuracy of the numerical integration.     

Building the terrestrial planets: Constrained accretion in the inner Solar System

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include: (1) the orbits, in particular the small eccentricities, and (2) the masses of the terrestrial planets - Mars’ relatively small mass in particular has not been adequately reproduced in previous simulations; (3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; (4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and (5) Earth’s relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N = 1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the “Nice model”. We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars’ mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. Our simulations leave us with a problem: we can reasonably satisfy the observed constraints (except for Earth’s water) with a configuration of Jupiter and Saturn that is at best marginally consistent with models of the outer Solar System, as it does not allow for any outer planet migration after a few Myr. Alternately, giant planet configurations which are consistent with the Nice model fail to reproduce Mars’ small size.     

Generation of equatorial jets by large-scale latent heating on the giant planets

Three-dimensional numerical simulations show that large-scale latent heating resulting from condensation of water vapor can produce multiple zonal jets similar to those on the gas giants (Jupiter and Saturn) and ice giants (Uranus and Neptune). For plausible water abundances (3-5 times solar on Jupiter/Saturn and 30 times solar on Uranus/Neptune), our simulations produce ∼20 zonal jets for Jupiter and Saturn and 3 zonal jets on Uranus and Neptune, similar to the number of jets observed on these planets. Moreover, these Jupiter/Saturn cases produce equatorial superrotation whereas the Uranus/Neptune cases produce equatorial subrotation, consistent with the observed equatorial-jet direction on these planets. Sensitivity tests show that water abundance, planetary rotation rate, and planetary radius are all controlling factors, with water playing the most important role; modest water abundances, large planetary radii, and fast rotation rates favor equatorial superrotation, whereas large water abundances favor equatorial subrotation regardless of the planetary radius and rotation rate. Given the larger radii, faster rotation rates, and probable lower water abundances of Jupiter and Saturn relative to Uranus and Neptune, our simulations therefore provide a possible mechanism for the existence of equatorial superrotation on Jupiter and Saturn and the lack of superrotation on Uranus and Neptune. Nevertheless, Saturn poses a possible difficulty, as our simulations were unable to explain the unusually high speed (∼) of that planet’s superrotating jet. The zonal jets in our simulations exhibit modest violations of the barotropic and Charney-Stern stability criteria. Overall, our simulations, while idealized, support the idea that latent heating plays an important role in generating the jets on the giant planets.     

Numerical periodic orbits of charged grains around magnetic planets

We apply a numerical searching method to investigate three-dimensional periodic orbits of charged dust particles in planetary magnetospheres. A classic generalized Stormer model of magnetic planets along with the parameters of Saturn is employed. More periodic orbits are found, besides the already known circular periodic orbits in or parallel to the equatorial plane. We divide all these orbits into six categories based on their appearances. By calculating the characteristic multipliers of the orbits, we investigate the stabilities of these periodic orbits.     

Planetary theories with the aid of the expansions of elliptic functions

For coplanar circular orbits, the mutual perturbations between two bodies can be expressed in term of the argument of Jacobian elliptic functions instead of the difference of the mean longitudes. For a given pair of planets, such a change of time variable improves the convergence of the developments. At the first order of planetary masses an integration of Lagrange's equations for the osculating elements is performed. When compared to classical developments the results are reduced by an important factor. The method is then extended to the mutual perturbations of Jupiter and Saturn, at any order of planetary masses, either with Fourier series with two arguments, or with one argument solely, taking advantage of the close commensurability of the mean motions.     

Terrestrial planet formation with strong dynamical friction

We have performed 8 numerical simulations of the final stages of accretion of the terrestrial planets, each starting with over 5× more gravitationally interacting bodies than in any previous simulations. We use a bimodal initial population spanning the region from 0.3 to 4 AU with 25 roughly Mars-mass embryos and an equal mass of material in a population of ∼1000 smaller planetesimals, consistent with models of the oligarchic growth of protoplanetary embryos. Given the large number of small planetesimals in our simulations, we are able to more accurately treat the effects of dynamical friction during the accretion process. We find that dynamical friction can significantly lower the timescales for accretion of the terrestrial planets and leads to systems of terrestrial planets that are much less dynamically excited than in previous simulations with fewer initial bodies. In addition, we study the effects of the orbits of Jupiter and Saturn on the final planetary systems by running 4 of our simulations with the present, eccentric orbits of Jupiter and Saturn (the EJS simulations) and the other 4 using a nearly circular and co-planar Jupiter and Saturn as predicted in the Nice Model of the evolution of the outer Solar System [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465] (the CJS simulations). Our EJS simulations provide a better match to our Solar System in terms of the number and average mass of the final planets and the mass-weighted mean semi-major axis of the final planetary systems, although increased dynamical friction can potentially improve the fit of the CJS simulations as well. However, we find that in our EJS simulations, essentially no water-bearing material from the outer asteroid belt ends up in the final terrestrial planets, while a large amount is delivered in the CJS simulations. In addition, the terrestrial planets in the EJS simulations receive a late veneer of material after the last giant impact event that is likely too massive to reconcile with the siderophile abundances in the Earth's mantle, while the late veneer in the CJS simulations is much more consistent with geochemical evidence.     

A possible long-lived asteroid population at the equilateral Lagrangian points of Saturn

The Lagrangian equilateral points of a planetary orbit are points of equilibrium that trail at 60°, ahead (L4) or behind (L5), the trajectory of a planet. Jupiter is the only major planet in our Solar system harbouring a known population of asteroids at those locations. Here we report the existence of orbits close to the Lagrangian points of Saturn, stable at time-scales comparable to the age of the Solar system. By scaling with respect to the Trojan population we have estimated the number of objects that would populate the regions, which gives a significant figure. Moreover, mutual physical collisions over the age of the Solar system would be very rare, so the evaporation rate of this swarm arising from mutual interactions would be very low. A population of asteroids not self-collisionally evolved after their formation stage would be the first to be observed in our planetary system. Our present estimations are based on the assumption that the capture efficiency at Saturn's equilateral points is comparable with the one corresponding to Jupiter, thus our figures may be taken as upper limits. In any case, observational constraints on their number would provide fundamental clues to our understanding of the history of the outer Solar system. If they existed, the surface properties and size distribution of those objects would represent unusually valuable fossil records of our early planetary system.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation when the two planets are in migration

In a previous paper (Hou et al. in Celest Mech Dyn Astron 119:119–142, 2014a), the problem of dynamical symmetry between two Jupiter triangular libration points (TLPs) with Saturn’s perturbation in the present configuration of the two planets was studied. A small short-time scale spatial asymmetry exists but gradually disappears with the time going, so the planar stable regions around the two Jupiter TLPs should be dynamically symmetric from a longtime perspective. In this paper, the symmetry problem is studied when the two planets are in migration. Several mechanisms that can cause asymmetries are discussed. Studies show that three important ones are the large short-time scale spatial asymmetry when Jupiter and Saturn are in resonance, the changing orbits of Jupiter and Saturn in the planet migration process, and the chaotic nature of Trojan orbits during the planet migration process. Their joint effects can cause an observable difference to the two Jupiter Trojan swarms. The thermal Yarkovsky effect is also found to be able to cause dynamical differences to the two TLPs, but generally they are too small to be practically observed.     

Planets in the asteroid belt

Abstract— The main asteroid belt has lost >99.9% of its solid mass since the time at which the planets were forming, according to models for the protoplanetary nebula. Here we show that the primordial asteroid belt could have been cleared efficiently if much of the original mass accreted to form planetsized bodies, which were capable of perturbing one another into unstable orbits. We provide results from 25 N‐body integrations of up to 200 planets in the asteroid belt, with individual masses in the range 0.017–0.33 Earth masses. In the simulations, these bodies undergo repeated close encounters which scatter one another into unstable resonances with the giant planets, leading to collision with the Sun or ejection from the solar system. In response, the giant planets' orbits migrate radially and become more circular. This reduces the size of the main‐belt resonances and the clearing rate, although clearing continues. If ～3 Earth masses of material was removed from the belt this way, Jupiter and Saturn would initially have had orbital eccentricities almost twice their current values. Such orbits would have made Jupiter and Saturn 10–100x more effective at clearing material from the belt than they are on their current orbits. The time required to remove 90% of the initial mass from the belt depends sensitively on the giant planets' orbits, and weakly on the masses of the asteroidal planets. 18 of the 25 simulations end with no planets left in the belt, and the clearing takes up to several hundred million years. Typically, the last one or two asteroidal planets are removed by interactions with planets in the terrestrial region     

2-Gyr Simulation of the Oort-Cloud Formation. I. Introduction on a New Model of the Outer Oort-Cloud Formation

The study of the reservoirs of small bodies in the Solar System can help us to refine our theory of the origin and evolution of the whole planetary system we live in. In this contribution, we introduce a simulation of the evolution of an initial proto-planetary disc (PPD) for 2 Gyr period, in which 10,038 studied test particles (TPs), representing the disc, are perturbed by four giant planets in their current orbits and having their current masses. In addition, Galactic-tide and stellar perturbations are considered. The outer Oort cloud (OC) appears to be formed within about 0.5 Gyr. At 2 Gyr, only 0.14% of all TPs reside in the outer OC, according to our simulation. This is the largest discrepancy with the results of previous authors. The TPs in the outer OC originate from almost all regions of the PPD. Specifically, 14, 7, 29, 21, and 29% of all TPs, that are the members of the outer OC at 2 Gyr, originate in the Jupiter, Saturn, Uranus, Neptune, and Edgeworth-Kuiper-belt regions, respectively. The largest number of the TPs, 6,669 (66.4%), was ejected into the interstellar space. Besides other results, we found a dominance of high galactic inclinations of outer-OC orbits. Obviously, this is a consequence of the action of Galactic tide.     

Dynamical evolution of a cometary swarm in the outer planetary region

The dynamical evolution of bodies under the gravitational influence of the accreting proto-Uranus and proto-Neptune is investigated. The main aim of this study is to analyze the interrelations between the accretion of Uranus and Neptune with other processes of cosmological importance as, for example, the formation of a cometary reservoir from bodies placed into near-parabolic orbits by planetary perturbations and the scattering of bodies to the region of the terrestrial planets. Starting with a mass ratio (initial mass/present mass) of 0.1, Uranus and Neptune acquire masses close to their present ones in a time scale of 10⁸ years. Neptune is found to be the most important contributor of comets to the cometary reservoir. The time scale of bodies scattered by Neptune to reach near-parabolic orbits (semimajor axes a > 10⁴ AU)is about 10⁹ years. The contribution of Uranus was partially inhibited because a large part of the residual bodies of its accretion zone fell under the strong gravitational influence of Jupiter and Saturn. A significant fraction of the bodies dispersed by Uranus and Neptune reached the region of the terrestrial planets in a time scale of some 10⁸ years.     

Restricted problem stability and the solar system

Reference periodic orbits are determined accurately for the planets of the solar system, using a restricted problem model with the Sun and Jupiter as the two primaries. The prediction is verified that stability of the planetary orbits should imply stability of those reference orbits that are simple-periodic.     

Kinetic models for the exospheres of Jupiter and Saturn

The exospheric theory based on the Kappa velocity distribution function (VDF) is used to model the exosphere of the giant planets Jupiter and Saturn. Such Kappa velocity distribution functions with an enhanced population of suprathermal particles are indeed often observed in space plasmas and in the space environment of the planets. The suprathermal particles have significant effects on the escape flux, density and temperature profiles of the particles in the exosphere of the giant planets. The polar wind flux becomes several orders larger when suprathermal electrons are considered, so that the planetary ionosphere becomes then a significant source for their inner magnetosphere. Moreover, the number density of the particles decreases slower as a function of the altitude when a Kappa distribution is considered instead of a Maxwellian one. Two-dimensional maps of density are calculated for typical values of the temperatures. The exospheric formalism is also applied to study the escape flux from the exospheres of Io and Titan, respectively, moons of Jupiter and Saturn.     

The transmission of mass and angular momentum from a satellite or planetary system to its primary

In this paper main implication of basic properties detected in the satellite systems of Jupiter, Saturn and Uranus, and presented by the author in an earlier contribution (Barricelli, 1971b) are investigated. The similarity between the primary periods in the three systems, their apparent relation to the axial rotation periods of the three planets and other features suggesting that collisions with the planetary surfaces may have played a role in the evolution of the three satellite systems are interpreted by assuming that in each case a satellite of unusually large size was originally disintegrated at the Roche limit of its primary. The disintegration of large satellites and their fusion with the respective planets is assumed to be a normal feature in the latest stage of planetary growth and the main cause of axial rotation in the respective planets.These assumptions make it possible to give a selfconsistent interpretation of the similarity between the axial rotation periods of the three planets and their relation to the primary periods (as defined by Barricelli, 1971b) in the three systems.Similar assumptions when applied to the Earth-Moon system make it possible to understand why the Moon, in its closest approach to the Earth is found to have been almost exactly at the Roche limit (Gerstenkorn, 1955; MacDonald, 1964), a coincidence which is too good to be accidental. According to this interpretation our Moon is a portion (representing about one third) of our original satellite, which survived its approach to the Roche limit and the ensuing fusion process with the Earth. It can be shown (see text) that under certain conditions this could leave a residual satellite with a stationary distance from the Earth (which in retrospect would be identified as its lowest distance from the Earth) at the Roche limit.The only other case in which we have observational evidence of parts of a satellite surviving its fusion process at the Roche limit is represented by the rings of Saturn and possibly the small innermost satellite Janus which seems to have been feeding on the rings.     

Turbulent convection in rapidly rotating spherical shells: A model for equatorial and high latitude jets on Jupiter and Saturn

The origin of zonal jets on the jovian planets has long been a topic of scientific debate. In this paper we show that deep convection in a spherical shell can generate zonal flow comparable to that observed on Jupiter and Saturn, including a broad prograde equatorial jet and multiple alternating jets at higher latitudes. We present fully turbulent, 3D spherical numerical simulations of rapidly rotating convection with different spherical shell geometries. The resulting global flow fields tend to be segregated into three regions (north, equatorial, and south), bounded by the tangent cylinder that circumscribes the inner boundary equator. In all of our simulations a strong prograde equatorial jet forms outside the tangent cylinder, whereas multiple jets form in the northern and southern hemispheres, inside the tangent cylinder. The jet scaling of our numerical models and of Jupiter and Saturn is consistent with the theory of geostrophic turbulence, which we extend to include the effect of spherical shell geometry. Zonal flow in a spherical shell is distinguished from that in a full sphere or a shallow layer by the effect of the tangent cylinder, which marks a reversal in the sign of the planetary β-parameter and a jump in the Rhines length. This jump is manifest in the numerical simulations as a sharp equatorward increase in jet widths—a transition that is also observed on Jupiter and Saturn. The location of this transition gives an estimate of the depth of zonal flow, which seems to be consistent with current models of the jovian and saturnian interiors.     

Solar activity and planetary luminosity

The Chree superposition analysis of the luminosities of the planets Jupiter, Saturn, Uranus and Neptune indicates a correlation between solar activity and planetary luminosity. The variations of the solar constant in the visible range are considered to be too small to explain the observed changes in brightness. The interaction of solar extreme ultraviolet or solar wind particles with the atmospheres of these planets is probably responsible for the increased albedo during periods of high solar activity.     

Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2009

Every three years the IAU Working Group on Cartographic Coordinates and Rotational Elements revises tables giving the directions of the poles of rotation and the prime meridians of the planets, satellites, minor planets, and comets. This report takes into account the IAU Working Group for Planetary System Nomenclature (WGPSN) and the IAU Committee on Small Body Nomenclature (CSBN) definition of dwarf planets, introduces improved values for the pole and rotation rate of Mercury, returns the rotation rate of Jupiter to a previous value, introduces improved values for the rotation of five satellites of Saturn, and adds the equatorial radius of the Sun for comparison. It also adds or updates size and shape information for the Earth, Mars?? satellites Deimos and Phobos, the four Galilean satellites of Jupiter, and 22 satellites of Saturn. Pole, rotation, and size information has been added for the asteroids (21) Lutetia, (511) Davida, and (2867) ?teins. Pole and rotation information has been added for (2) Pallas and (21) Lutetia. Pole and rotation and mean radius information has been added for (1) Ceres. Pole information has been updated for (4) Vesta. The high precision realization for the pole and rotation rate of the Moon is updated. Alternative orientation models for Mars, Jupiter, and Saturn are noted. The Working Group also reaffirms that once an observable feature at a defined longitude is chosen, a longitude definition origin should not change except under unusual circumstances. It is also noted that alternative coordinate systems may exist for various (e.g. dynamical) purposes, but specific cartographic coordinate system information continues to be recommended for each body. The Working Group elaborates on its purpose, and also announces its plans to occasionally provide limited updates to its recommendations via its website, in order to address community needs for some updates more often than every 3 years. Brief recommendations are also made to the general planetary community regarding the need for controlled products, and improved or consensus rotation models for Mars, Jupiter, and Saturn.