Origin of water in the terrestrial planets

Abstract— I examine the origin of water in the terrestrial planets. Late‐stage delivery of water from asteroidal and cometary sources appears to be ruled out by isotopic and molecular ratio considerations, unless either comets and asteroids currently sampled spectroscopically and by meteorites are unlike those falling to Earth 4.5 Ga ago, or our measurements are not representative of those bodies. However, the terrestrial planets were bathed in a gas of H, He, and O. The dominant gas phase species were H₂, He, H₂ O, and CO. Thus, grains in the accretion disk must have been exposed to and adsorbed H2 and water. Here I conduct a preliminary analysis of the efficacy of nebular gas adsorption as a mechanism by which the terrestrial planets accreted “wet.” A simple model suggests that grains accreted to Earth could have adsorbed 1‐3 Earth oceans of water. The fraction of this water retained during accretion is unknown, but these results suggest that examining the role of adsorption of water vapor onto grains in the accretion disk bears further study.     

Computer simulation of planet formation in a binary star system: Terrestrial planets

A model of planetary formation in a binary system with a small relative mass of primary is computed on the assumption of a mass transfer from the less massive component to the more massive one with no mass and angular momentum carried away from the system under consideration. At the last stage of mass transfer the condensed Moon-like objects (planetoids) are ejected through the inner Lagrange point of the primary Roche lobe with the outflow of gaseous matter.The whole system is considered in the plane of binary star rotation. Newtonian equations of motion are integrated with the initial conditions for the planetoids referred to as the coordinates and velocity of the inner Lagrangian point at the moments of planetoid ejections, all the pairwise gravitational interactions being included in computations but without a gas-drag. The mass transfer ceases at the primary relative mass 10^–3 which corresponds to the present Sun-Jupiter system. The total mass of planetoids approximates that of the terrestrial planets. Those are formed through coagulation of the planetoids with the effective radius of capture cross-section as an input parameter in the computer simulation. When the minimum separation between the pair of bodies becomes less than this radius they coalesce into a single body with their masses and momenta summed. If the effective radius value is under a certain limit the computer simulation yields the planetary system like that of terrestrial planets of the present Sun system.Numerical computations reveal the division of the planetoids into 4 groups along their distances from the Sun. Further, each group forms a single planet or a planet and a less massive body at the nearest orbits. The parameters of simulated planet orbits are close to the present ones and the interplanetary spacings are in accord with the Titius-Bode law.     

Meteoritical and dynamical constraints on the growth mechanisms and formation times of asteroids and Jupiter

Thermal models and radiometric ages for meteorites show that the peak temperatures inside their parent bodies were closely linked to their accretion times. Most iron meteorites come from bodies that accreted <0.5 Myr after CAIs formed and were melted by ²⁶Al and ⁶⁰Fe, probably inside 2 AU. Rare carbon-rich differentiated meteorites like ureilites probably also come from bodies that formed <1 Myr after CAIs, but in the outer part of the asteroid belt. Chondrite groups accreted intermittently from diverse batches of chondrules and other materials over a 4 Myr period starting 1 Myr after CAI formation when planetary embryos may already have formed at ∼1 AU. Meteorite evidence precludes accretion of late-forming chondrites on the surface of early-formed bodies; instead chondritic and non-chondritic meteorites probably formed in separate planetesimals. Maximum metamorphic temperatures in chondrite groups are correlated with mean chondrule age, as expected if ²⁶Al and ⁶⁰Fe were the predominant heat sources. Because late-forming bodies could not accrete close to large, early-formed bodies, planetesimal formation may have spread across the nebula from regions where the differentiated bodies formed. Dynamical models suggest that the asteroids could not have accreted in the main belt if Jupiter formed before the asteroids. Therefore Jupiter probably reached its current mass >3-5 Myr after CAIs formed. This precludes formation of Jupiter via a gravitational instability <1 Myr after the solar nebula formed, and strongly favors core accretion. Jupiter probably formed too late to make chondrules by generating shocks directly, or indirectly by scattering Ceres-sized bodies across the belt. Nevertheless, shocks formed by gravitational instabilities or Ceres-sized bodies scattered by planetary embryos may have produced some chondrules. The minimum lifetime for the solar nebula of 3-5 Myr inferred from the total spread of CAI and chondrule ages may exceed the median lifetime of 3 Myr for protoplanetary disks, but is well within the 1-10 Myr observed range. Shorter formation times for extrasolar planets may help to explain their unusual orbits compared to those of solar giant planets.     

Accretion of terrestrial planets from oligarchs in a turbulent disk

We have investigated the final accretion stage of terrestrial planets from Mars-mass protoplanets that formed through oligarchic growth in a disk comparable to the minimum mass solar nebula (MMSN), through N-body simulation including random torques exerted by disk turbulence due to Magneto-Rotational Instability. For the torques, we used the semi-analytical formula developed by Laughlin et al. [Laughlin, G., Steinacker, A., Adams, F.C., 2004. Astrophys. J. 608, 489-496]. The damping of orbital eccentricities (in all runs) and type-I migration (in some runs) due to the tidal interactions with disk gas is also included. Without any effect of disk gas, Earth-mass planets are formed in terrestrial planet regions in a disk comparable to MMSN but with too large orbital eccentricities to be consistent with the present eccentricities of Earth and Venus in our Solar System. With the eccentricity damping caused by the tidal interaction with a remnant gas disk, Earth-mass planets with eccentricities consistent with those of Earth and Venus are formed in a limited range of disk gas surface density (∼10⁻⁴ times MMSN). However, in this case, on average, too many (?6) planets remain in terrestrial planet regions, because the damping leads to isolation between the planets. We have carried out a series of N-body simulations including the random torques with different disk surface density and strength of turbulence. We found that the orbital eccentricities pumped up by the turbulent torques and associated random walks in semimajor axes tend to delay isolation of planets, resulting in more coagulation of planets. The eccentricities are still damped after planets become isolated. As a result, the number of final planets decreases with increase in strength of the turbulence, while Earth-mass planets with small eccentricities are still formed. In the case of relatively strong turbulence, the number of final planets are 4-5 at 0.5-2 AU, which is more consistent with Solar System, for relatively wide range of disk gas surface density (∼¹⁰⁻⁴-¹⁰⁻² times MMSN).     

Nonuniform cratering of the Moon and a revised crater chronology of the inner Solar System

We model the cratering of the Moon and terrestrial planets from the present knowledge of the orbital and size distribution of asteroids and comets in the inner Solar System, in order to refine the crater chronology method. Impact occurrences, locations, velocities and incidence angles are calculated semi-analytically, and scaling laws are used to convert impactor sizes into crater sizes. Our approach is generalizable to other moons or planets. The lunar cratering rate varies with both latitude and longitude: with respect to the global average, it is about 25% lower at (±65°N, 90°E) and larger by the same amount at the apex of motion (0°N, 90°W) for the present Earth-Moon separation. The measured size-frequency distributions of lunar craters are reconciled with the observed population of near-Earth objects under the assumption that craters smaller than a few kilometers in diameter form in a porous megaregolith. Varying depths of this megaregolith between the mare and highlands is a plausible partial explanation for differences in previously reported measured size-frequency distributions. We give a revised analytical relationship between the number of craters and the age of a lunar surface. For the inner planets, expected size-frequency crater distributions are calculated that account for differences in impact conditions, and the age of a few key geologic units is given. We estimate the Orientale and Caloris basins to be 3.73 Ga old, and the surface of Venus to be 240 Ma old. The terrestrial cratering record is consistent with the revised chronology and a constant impact rate over the last 400 Ma. Better knowledge of the orbital dynamics, crater scaling laws and megaregolith properties are needed to confidently assess the net uncertainty of the model ages that result from the combination of numerous steps, from the observation of asteroids to the formation of craters. Our model may be inaccurate for periods prior to 3.5 Ga because of a different impactor population, or for craters smaller than a few kilometers on Mars and Mercury, due to the presence of subsurface ice and to the abundance of large secondaries, respectively. Standard parameter values allow for the first time to naturally reproduce both the size distribution and absolute number of lunar craters up to 3.5 Ga ago, and give self-consistent estimates of the planetary cratering rates relative to the Moon.     

Solar flare track densities in interplanetary dust particles: The determination of an asteroidal versus cometary source of the zodiacal dust cloud

Dust particles that are larger than 1 μm, when injected into the Solar System from comets and asteroids, will spiral into the Sun due to the Poynting-Robertson effect. During the process of spiraling in, such dust particles accumulate solar flare tracks in their component minerals. The accumulated track density for a given dust grain is a function of the duration of its space exposure and its distance from the Sun. Using a computer model, it was determined that the expected track density distributions from grains produced by comets are very different from those produced by asteroids. Individual asteroids produce populations of particles that arrive at 1 AU with scaled track density distributions containing “spikes,” while comets supply particles with a flatter and wider distribution of track densities. Particles with track densities above 3 × 10⁷ (sϱA/v) tracks/cm² have probably been exposed to solar flare tracks prior to injection into the interplanetary medium and are therefore likely to be asteroidal. Particles with track densities below 0.7 × 10⁷(sϱA/v) tracks/cm² must be derived from comets or Earth-crossing asteroids. Earth-crossing asteroids are not responsible for all the dust collected at 1 AU since they cannot produce the large track densities observed in some of the interplanetary dust particles collected in the stratosphere. The track densities observed in the stratospheric dust fall within the predicted range, but there is at present an insufficient number of carefully determined densities to make strong statements about the sources of the present dust population.     

Asteroid (21) Lutetia as a remnant of Earth’s precursor planetesimals

Isotopic and chemical compositions of meteorites, coupled with dynamical simulations, suggest that the main belt of asteroids between Mars and Jupiter contains objects formed in situ as well as a population of interlopers. These interlopers are predicted to include the building blocks of the terrestrial planets as well as objects that formed beyond Neptune ( [Bottke et al., 2006] , [Levison et al., 2009] and [Walsh et al., 2011] ). Here we report that the main belt asteroid (21) Lutetia – encountered by the Rosetta spacecraft in July 2010 – has spectral (from 0.3 to 25 μm) and physical (albedo, density) properties quantitatively similar to the class of meteorites known as enstatite chondrites. The chemical and isotopic compositions of these chondrites indicate that they were an important component of the formation of Earth and other terrestrial planets. This meteoritic association implies that Lutetia is a member of a small population of planetesimals that formed in the terrestrial planet region and that has been scattered in the main belt by emerging protoplanets (Bottke et al. 2006) and/or by the migration of Jupiter (Walsh et al. 2011) early in its history. Lutetia, along with a few other main-belt asteroids, may contains part of the long-sought precursor material (or closely related materials) from which the terrestrial planets accreted.     

Planetary impact probabilities for long-period comets

Planetary impact probabilities for long-period (near-parabolic) comets are determined by averaging Öpik's equations over inclination and perihelion distance for each planet. These averaged values compare well with the results of more elaborate Monte Carlo calculations. The impact probabilities are proportional to the square of the normalized capture radius of each planet, which in turn is a function of the planet's radius and mass, so that the major planets have the highest impact probabilities. Encounter velocities have an average value of $\text{3}\text{1}\text{2}$ times the planetary orbital velocity but the most probable encounter velocities are slightly higher than this for the terrestrial planets and slightly lower for the major planets. Comparison of the impact probabilities with the cratering record, corrected for gravity and velocity effects, indicates that long-period comets may account for 3 to 9% of the observed large crattes (diameter > 10 km) on the terrestrial planets. The inclination and perihelion properties of the impact probabilities obtained from numerical averaging provide a simple method for determining the impact probabilities for nonuniform distributions. The perihelion distribution of long period comets from J. A. Fernandez ((1981) Astron. Astrophys.96, 26–35) results in a crater production rate quite similar throughout the solar system, unlike that of a uniform perihelion distribution.     

Refractory Metal Nuggets – Formation of the first condensates in the Solar Nebula

As gas flowed from the solar accretion disk or Solar Nebula onto the proto-Sun, magnetic pressure gradients in the solar magnetosphere and the inner Solar Nebula provided an environment where some of this infalling flow was diverted to produce a low pressure, high temperature, gaseous, “infall” atmosphere around the inner Solar Nebula. The pressure in this inner disk atmosphere was mainly dependant on the accretion flow rate onto the star. High flow rates implied relatively high pressures, which decreased over time as the accretion rate decreased.In the first hundred thousand years after the formation of the Solar Nebula, accretional flow gas pressures were high enough to create submicron-sized Refractory Metal Nuggets (RMNs) – the precursors to Calcium Aluminum Inclusions (CAIs). Optimal temperatures and pressures for RMN formation may have occurred between 20,000 and 100,000 years after the formation of the Solar Nebula. It is possible that conditions were conducive to RMN/CAI formation over an 80,000 year timescale. The “infall” atmosphere and the condensation of refractory particles within this atmosphere may be observable around the inner disks of other protostellar systems.The interaction of forces from magnetic fields with the radiation pressure from the proto-Sun and the inner solar accretion disk potentially produced an optical-magnetic trap above and below the inner Solar Nebula, which provided a relatively stable environment in which the RMNs/proto-CAIs could form and grow. These RMN formation sites only existed during accretion events from the proto-solar disk onto the proto-Sun. As such, the formation and growth time of a particular RMN was dependent on the timescale of its nascent accretion event.Observational evidence suggests that RMNs were the nucleation particles for CAIs. As a consequence, the observed bimodal distribution of ²⁶Al in CAIs, where some CAIs have ²⁶Al while others do not, is probably due to the injection ²⁶Al during the short CAI formation period, where ²⁶Al was not present when the first CAIs were formed.     

The “Shiva Hypothesis”: Impacts, mass extinctions, and the galaxy

The “Shiva Hypothesis”, in which recurrent, cyclical mass extinctions of life on Earth result from impacts of comets or asteroids, provides a possible unification of important processes in astrophysics, planetary geology, and the history of life. Collisions with Earth-crossing asteroids and comets ≥ a few km in diameter are calculated to produce widespread environmental disasters (dust clouds, wildfires), and occur with the proper frequency to account for the record of five major mass extinctions (from ≥ 10⁸ Mt TNT impacts) and ~ 20 minor mass extinctions (from 10⁷–10⁸ Mt impacts) recorded in the past 540 million years. Recent studies of a number of extinctions show evidence of severe environmental disturbances and mass mortality consistent with the expected after-effects (dust clouds, wildfires) of catastrophic impacts. At least six cases of features generally considered diagnostic of large impacts (e.g., large impact craters, layers with high platinum-group elements, shock-related minerals, and/or microtektites) are known at or close to extinction-event boundaries. Six additional cases of elevated iridium levels at or near extinction boundaries are of the amplitude that might be expected from collision of relatively low-Ir objects such as comets. The records of cratering and mass extinction show a correlation, and might be explained by a combination of periodic and stochastic impactors. The mass extinction record shows evidence for a periodic component of about 26 to 30 Myr, and an ~ 30 Myr periodic component has been detected in impact craters by some workers, with recent pulses of impacts in the last 2–3 million years, and at ~ 35, 65, and 95 million years ago. A cyclical astronomical pacemaker for such pulses of impacts may involve the motions of the Earth through the Milky Way Galaxy. As the Solar System revolves around the galactic center, it also oscillates up and down through the plane of the disk-shaped galaxy with a half-cycle ~ 30±3 Myr. This cycle should lead to quasi-periodic encounters with interstellar clouds, and periodic variations in the galactic tidal force with maxima at times of plane crossing. This “galactic carrousel” effect may provide a viable perturber of the Oort Cloud comets, producing periodic showers of comets in the inner Solar System. These impact pulses, along with stochastic impactors, may represent the major punctuations in earth history.     

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.     

Cratering of the terrestrial planets by Apollo objects

Abstract— Absolute and relative cratering rates on the terrestrial planets have been calculated using the same asteroidal collision model and Monte Carlo program used for previous studies of the terrestrial meteorite flux, the steady-state number of Apollo-Amor objects, and the orbital distribution of both meteorites and Apollo-Amor objects. The most straightforward result is that projectiles from the asteroid belt appear to provide about one-third the observed present-day production of terrestrial craters larger than 10 km in diameter. When uncertainties in the calculations and observations are included, it cannot be excluded that the entire terrestrial cratering flux is asteroidal. On the other hand, assumption of an additional Apollo-Amor source of extinct comets, in the same quantity permitted by Apollo-Amor observations, provides better agreement with the observed cratering rate. In addition, a significant (e.g., ～30%) terrestrial contribution from active long and short period comets is acceptable within the uncertainties of the assumptions required. The ratios of the cratering rates on the different terrestrial planets are somewhat sensitive to the assumed source. A purely asteroidal source predicts a martian cratering rate per unit area about four times that on Earth, whereas the difference is reduced to about a factor of two for the mixed asteroid-extinct comet source. The opposite effect is found for Mercury. As discussed by previous authors, the predicted lunar cratering rate is significantly higher than that observed. It is not clear whether this is a result of scaling to impacts on a body considerably smaller than Earth, or if it indicates an increase in the cratering flux during the Phanerozoic.     

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.     

Hidden Mass in the Asteroid Belt

The total mass of the asteroid belt is estimated from an analysis of the motions of the major planets by processing high precision measurements of ranging to the landers Viking-1, Viking-2, and Pathfinder (1976-1997). Modeling of the perturbing accelerations of the major planets accounts for individual contributions of 300 minor planets; the total contribution of all remaining small asteroids is modeled as an acceleration caused by a solid ring in the ecliptic plane. Mass M_ring of the ring and its radius R are considered as solve-for parameters. Masses of the 300 perturbing asteroids have been derived from their published radii based mainly on measured fluxes of radiation, making use of the corresponding densities. This set of asteroids is grouped into three classes in accordance with physical properties and then corrections to the mean density for each class are estimated in the process of treating the observations. In this way an improved system of masses of the perturbing asteroids has been derived.The estimate M_ring≈(5±1)×10⁻¹⁰M_⊙ is obtained (M_⊙ is the solar mass) whose value is about one mass of Ceres. For the mean radius of the ring we have R≈2.80 AU with 3% uncertainty. Then the total mass M_belt of the main asteroid belt (including the 300 asteroids mentioned above) may be derived: M_belt≈(18±2)×10⁻¹⁰M_⊙. The value M_belt includes masses of the asteroids which are already discovered, and the total mass of a large number of small asteroids—most of which cannot be observed from the Earth. The second component M_ring is the hidden mass in the asteroid belt as evaluated from its dynamical impact onto the motion of the major planets.Two parameters of a theoretical distribution of the number of asteroids over their masses are evaluated by fitting to the improved set of masses of the 300 asteroids (assuming that there is no observational selection effect in this set). This distribution is extrapolated to the whole interval of asteroid masses and as a result the independent estimate M_belt≈18×10⁻¹⁰M_⊙ is obtained which is in excellent agreement with the dynamical finding given above.These results make it possible to predict the total number of minor planets in any unit interval of absolute magnitude H. Such predictions are compared with the observed distribution; the comparison shows that at present only about 10% of the asteroids with absolute magnitude H<14 have been discovered (according to the derived distribution, about 130,000 such asteroids are expected to exist).     

Complex of Meteoroid Orbits with Eccentricities Near 1 and Higher

In our work, the method that can help to predict the existence of distant objects in the Solar system is demonstrated. This method is connected with statistical properties of a heliocentric orbital complex of meteoroids with high eccentricities. Heliocentric meteoroid orbits with high eccentricities are escape routes for dust material from distant parental objects with near-circular orbits to Earth-crossing orbits. Ground-based meteor observations yield trajectory information from which we can derive their place of possible origin: comets, asteroids, and other objects (e.g. Kuiper Objects) in the Solar system or even interstellar space. Statistical distributions of radius vectors of nodes, and other parameters of orbits of meteoroids contain key information about position of greater bodies. We analyze meteor orbits with high eccentricities that were registered in 1975–1976 in Kharkiv (Ukraine). The orbital data of the Kharkiv electronic catalogue are received from observations of radiometeors with masses 10⁻⁶−10⁻³ g.     

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.     

Dynamical and collisional evolution of Halley-type comets

The number of observed Halley-type comets is hundreds of times less than predicted by models (Levison, H.F., Dones, L., Duncan, M.J. [2001]. Astron. J. 121, 2253–2267). In this paper we investigate the impact of collisions with planetesimals on the evolution of Halley-type comets. First we compute the dynamical evolution of a sub-set of 21 comets using the Mercury integrator package over 100 Myr. The dynamical lifetime is determined to be of the order of 10⁵–10⁶ years in agreement with previous work. The collisional probability of Halley-type comets colliding with known asteroids, a simulated population of Kuiper-belt objects, and planets, is calculated using a modified, Öpik-based collision code. Our results show that the catastrophic disruption of the cometary nucleus has a very low probability of occurring, and disruption through cumulative minor impacts is concluded to be negligible. The dust mantle formed from ejected material falling back to the comet’s surface is calculated to be less than a few centimeters thick, which is insignificant compared to the mantle formed by volatile depletion, while planetary encounters were found to be a negligible disruption mechanism.     

On the origin of the unusual orbit of Comet 2P/Encke

The orbit of Comet 2P/Encke is difficult to understand because it is decoupled from Jupiter—its aphelion distance is only 4.1 AU. We present a series of orbital integrations designed to determine whether the orbit of Comet 2P/Encke can simply be the result of gravitational interactions between Jupiter-family comets and the terrestrial planets. To accomplish this, we integrated the orbits of a large number of objects from the trans-neptunian region, through the realm of the giant planets, and into the inner Solar System. We find that at any one time, our model predicts that there should be roughly 12 objects in Encke-like orbits. However, it takes roughly 200 times longer to evolve onto an orbit like this than the typical cometary physical lifetime. Thus, we suggest that (i) 2P/Encke became dormant soon after it was kicked inward by Jupiter, (ii) it spent a significant amount of time inactive while rattling around the inner Solar System, and (iii) it only became active again as the ν₆ secular resonance drove down its perihelion distance.     

Model consideration of the bombardment event of the asteroidal belt by the planetesimals scattered from the Jupiter zone

The temporal evolutions of the planetesimals scattered from the Jupiter zone for different masses of the proto-Jupiter [(a) 0.1 and (b) 1.0 of the present mass] are investigated. Due to the combined effects of the orbital evolution of the planetesimals and the elimination of these projectiles either via impact capture or injection into escape velocity by the outer planets, the whole scattering process lasts about 10⁸ yr for case (a) and about 10⁷ yr for case (b). The longer time scale may be a good estimate for the accretion time interval of Jupiter while the shorter one (10⁷) gives the upper time limit of the late heavy-bombardment epoch of the terrestrial planets due to planetesimals scattered from the Jupiter zone. The limiting value of the encounter velocity U at the end of the scattering process is ≈0.6. Consideration of the collisional interaction of these projectiles with the asteroids indicates that the corresponding bombardment effect could be rather appreciable. Also, the asteroids on the inner edge of the main asteroid belt would have been bombarded more severely than those on the outer edge. From this point of view, the structure of the asteroidal belt could be affected significantly not only by Jupiter's gravitational perturbation effect but also by its early scattering process.     

Origin of water in the inner Solar System: A kinetic Monte Carlo study of water adsorption on forsterite

The origin of water in the inner Solar System is not well understood. It is believed that temperatures were too high in the accretion disk in the region of the terrestrial planets for hydrous phases to be thermodynamically stable. Suggested sources of water include direct adsorption of hydrogen from the nebula into magma oceans after the terrestrial planets formed, and delivery of asteroidal or cometary material from beyond the zone of the terrestrial planets. We explore a new idea, direct adsorption of water onto grains prior to planetary accretion. This hypothesis is motivated by the observation that the accretion disk from which our planetary system formed was composed of solid grains bathed in a gas dominated by hydrogen, helium, and oxygen. Some of that hydrogen and oxygen combined to make water vapor. We examine quantitatively adsorption of water onto grains in the inner Solar System accretion disk by exploring the adsorption dynamics of water molecules onto forsterite surfaces via kinetic Monte Carlo simulations. We conclude that many Earth oceans of water could be adsorbed.