Planet formation: Compositional mixing and lunar compositional anomalies

Significant fractions of each planet's late-accreted mass originated not at its own distance from the Sun, but from a neighboring planet's orbit, according to results that follow from calculations by Wetherill (1975). “Late-accreted” refers to a loosely defined period after planets acquired most of their present mass. In an idealized model, Mercury, Venus, Earth, and Mars received 47, 45, 37, and 52% of their late-accreted mass from planetesimals formed closer to other planets. Resulting compositional anomalies in outer parts of early planets could be significant; atmospheric tests of Lewis's predicted S deficiency on Venus may be inconclusive.The Moon's orbit around Earth puts it in a special category: sorting occurs between Moon-impacting and Earth-impacting material according to approach velocity. In the above model, the moon receives 60% of its late-accreted mass from planetesimals formed near Venus' orbit. Distant planetesimals could be perturbed into the Earth-Moon system and cause major changes in the Moon's composition with only minor effect on Earth. The entire lunar bulk composition anomaly could be explained by plausible reservoirs of distant low-density material.     

Asteroids were born big

How big were the first planetesimals? We attempt to answer this question by conducting coagulation simulations in which the planetesimals grow by mutual collisions and form larger bodies and planetary embryos. The size frequency distribution (SFD) of the initial planetesimals is considered a free parameter in these simulations, and we search for the one that produces at the end objects with a SFD that is consistent with Asteroid belt constraints. We find that, if the initial planetesimals were small (e.g. km-sized), the final SFD fails to fulfill these constraints. In particular, reproducing the bump observed at diameter in the current SFD of the asteroids requires that the minimal size of the initial planetesimals was also ∼100 km. This supports the idea that planetesimals formed big, namely that the size of solids in the proto-planetary disk “jumped” from sub-meter scale to multi-kilometer scale, without passing through intermediate values. Moreover, we find evidence that the initial planetesimals had to have sizes ranging from 100 to several 100 km, probably even 1000 km, and that their SFD had to have a slope over this interval that was similar to the one characterizing the current asteroids in the same size range. This result sets a new constraint on planetesimal formation models and opens new perspectives for the investigation of the collisional evolution in the Asteroid and Kuiper belts as well as of the accretion of the cores of the giant planets.     

The fossilized size distribution of the main asteroid belt

Planet formation models suggest the primordial main belt experienced a short but intense period of collisional evolution shortly after the formation of planetary embryos. This period is believed to have lasted until Jupiter reached its full size, when dynamical processes (e.g., sweeping resonances, excitation via planetary embryos) ejected most planetesimals from the main belt zone. The few planetesimals left behind continued to undergo comminution at a reduced rate until the present day. We investigated how this scenario affects the main belt size distribution over Solar System history using a collisional evolution model (CoEM) that accounts for these events. CoEM does not explicitly include results from dynamical models, but instead treats the unknown size of the primordial main belt and the nature/timing of its dynamical depletion using innovative but approximate methods. Model constraints were provided by the observed size frequency distribution of the asteroid belt, the observed population of asteroid families, the cratered surface of differentiated Asteroid (4) Vesta, and the relatively constant crater production rate of the Earth and Moon over the last 3 Gyr. Using CoEM, we solved for both the shape of the initial main belt size distribution after accretion and the asteroid disruption scaling law . In contrast to previous efforts, we find our derived function is very similar to results produced by numerical hydrocode simulations of asteroid impacts. Our best fit results suggest the asteroid belt experienced as much comminution over its early history as it has since it reached its low-mass state approximately 3.9-4.5 Ga. These results suggest the main belt's wavy-shaped size-frequency distribution is a “fossil” from this violent early epoch. We find that most diameter D?120 km asteroids are primordial, with their physical properties likely determined during the accretion epoch. Conversely, most smaller asteroids are byproducts of fragmentation events. The observed changes in the asteroid spin rate and lightcurve distributions near D∼100-120 km are likely to be a byproduct of this difference. Estimates based on our results imply the primordial main belt population (in the form of D<1000 km bodies) was 150-250 times larger than it is today, in agreement with recent dynamical simulations.     

Origin and evolution of the earth-moon system

The origin and evolution of the Earth-Moon system is studied by comparing it to the satellite systems of other planets. The normal structure of a system of secondary bodies orbiting around a central body depends essentially on the mass of the central body. The Earth with a mass intermediate between Uranus and Mars should have a normal satellite system that consists of about half a dozen satellites each with a mass of a fraction of a percent of the lunar mass. Hence, the Moon is not likely to have been generated in the environment of the Earth by a normal accretion process as is claimed by some authors.Capture of satellites is quite a common process as shown by the fact that there are six satellites in the solar system which, because they are retrograde, must have been captured. There is little doubt that the Moon is also a captured satellite, but its capture orbit and tidal evolution are still incompletely understood.The Earth and the Moon are likely to have been formed from planetesimals accreting in particle swarms in Kepler orbits (jet streams). This process leads to the formation of a cool lunar interior with an outer layer accreted at increasingly higher temperatures. The primeval Earth should similarly have formed with a cool inner core surrounded in this case by a very strongly heated outer core and with a mantle accreted slowly and with a low average temperature but with intense transient heating at each individual impact site.     

Micrometeorite impact annealing of ice in the outer Solar System

The spectra of water ice on the surfaces of icy satellites and Kuiper Belt Objects (KBOs) indicate that the surface ice on these bodies is in a crystalline state. This conflicts with theoretical models, which predict that radiation (galactic cosmic rays and solar ultraviolet) should damage the crystalline structure of ice on geologically short timescales. Temperatures are too low in the outer Solar System for the ice to anneal, and reflectance spectra of these bodies should match those of amorphous solid water (ASW). We assess whether the kinetic energy deposited as heat by micrometeorite impacts on outer Solar System bodies is sufficient to anneal their surface ice down to a near-infrared optical depth . We calculate the kinetic energy flux from interplanetary micrometeorite impacts, including gravitational focusing. We also calculate the thermal diffusion of impact heat in various surfaces and the rate of annealing of ice. We conclude that the rate of annealing from micrometeorite impacts is sufficient to explain the crystallinity of ice on nearly all the surfaces of the saturnian, uranian and neptunian satellites. We discuss how the model can be used in conjunction with spectra of KBOs to probe dust fluxes in the Kuiper Belt.     

High-resolution simulations of the final assembly of Earth-like planets I. Terrestrial accretion and dynamics

The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.     

Experimental study on the impact fragmentation of core-mantle bodies: Implications for collisional disruption of rocky planetesimals with sintered core covered with porous mantle

Impact experiments of inhomogeneous targets such as layered bodies consisting of a dense core and porous mantle were conducted to clarify the effect of the layered structure on impact strength. The layered structure of small bodies could be the result of the thermal evolution of planetesimals in the solar nebula. So, the impact disruption of thermally evolved bodies with core-mantle structure is important for the origin of small bodies such as asteroids. We investigated the impact strength of rocky-layered bodies with porous mantle-sintered cores, which could be formed at an initial stage of thermal evolution. Spherical targets composed of soda-lime glass or quartz core and porous gypsum mantle were prepared as an analog of small bodies with a core-mantle structure, and the internal structure was changed. A nylon projectile was impacted at the impact velocity from 1 to 5 km/s. The impact strength of the core-mantle targets decreases with the increase of the core/target mass ratio (RCM) in the specific energy range from 1×¹⁰³ to 4×¹⁰⁴ J/kg. We observed two distinct destruction modes characterized by the damage to the core: one shows a damaged core and fractured mantle, and the other shows an intact core and broken mantle. The former mode was usually observed with increasing RCM, and the boundary condition of the core destruction () was experimentally found to be , where is the specific energy required to disrupt a glass core. From this empirical equation, it might be possible to discuss the destruction conditions of a thermally evolved body with a porous mantle-sintered core structure. We speculate that the impact strength of the body could be significantly reduced with the progress of internal evolution at the initial stage of thermal evolution.     

Bombardment as a cause of the lunar asymmetry

A comparison of the lunar frontside gravity field with topography indicates that low-density ( 2.9 g cm^–3) types of rock form a surface layer or crust of variable thickness: 40-60 km beneath terrae; 20-40 km beneath non-mascon maria; 0-20 km beneath mascon maria. The observed offset between lunar centers of mass and figure is consistent with farside crustal thicknesses of 40-50 km, similar to frontside terra thicknesses.The Moon is asymmetric in crustal thickness, and also in the distribution of maria and gamma radioactivity. Early bombardment of the Moon by planetesimals, in both heliocentric and geocentric orbits, is examined as a possible cause of the asymmetries. The presence of a massive companion (Earth) causes a spin-orbit coupled Moon to be bombarded non-uniformly. The most pronounced local concentration of impacts would have occurred on the west limb of the Moon, when it orbited close to the Earth, if low-eccentricity heliocentric planetesimals were still abundant in the solar system at that time.A very intense bombardment of this type could have redistributed crustal material on the Moon, thinning the west limb crust appreciably. This would have caused a change in position of the principal axes of inertia, and a reorientation of the spin-orbit coupled Moon such that the thinnest portion of its crust turned toward one of the poles. Erupting lavas would have preferentially flooded such a thin-crusted, low-lying area. This would have caused another readjustment of principal moments, and a reorientation of the Moon such that the mare areas tipped toward the equator. The north-south and nearside-farside asymmetries of mare distribution on the present Moon can be understood in terms of such a history.Paper dedicated to Prof. Harold C. Urey on the occasion of his 80th birthday on 29 April 1973.     

A note on the limits of stability for the restricted problem of three bodies as applied to the Sun-Earth-Moon system

The Sun-Earth-Moon system is modeled by the restricted problem of three bodies, and the curves of zero velocity are used to define the limits of stability of the Moon's orbit about the Earth. By holding the relative distances fixed, and maintaining the circular velocities of the Earth and Moon while their masses are varied by a common factor (=m _E/m _E=m _M/m _M), it is found that the possibility of the Moon leaving Earth orbit and orbiting the Sun exists for the range of values 0.005<<0.4.     

Formation of the terrestrial planets

The early phases of formation in the inner solar system were dominated by collisions and short-range dynamical interactions among planetesimals. But the later phases, which account for most of the differences among planets, are unsure because the dynamics are more subtle. Jupiter's influence became more important, leading to drastic clearing out of the asteroid belt and the stunting of Mars's growth. Further in, the effect of Jupiter-- both directly and indirectly, through ejection of mass in the outer solar system-- was probably to speed up the process without greatly affecting the outcome. The great variety in bulk properties of the terrestrial bodies indicate a terminal phase of great collisions, so that the outcome is the result of small-N statistics. Mercury, 65 percent iron, appears to be a residual core from a high-velocity collision. All planets appear to require a late phase of high energy impacts to erode their atmospheres: including the Earth, to remove CO₂ so that its ocean could form by condensation of water.Consistent with this model is that the largest collision, about 0.2 Earth masses, was into the proto-Earth, although the only property that appears to require it is the great lack of iron in the Moon. The other large differences between the Earth and Venus, angular momentum (spin plus satellite) and inert gas abundances, must arise from origin circumstances, but neither require nor forbid the giant impact. Venus's higher ratio of light to heavy inert gases argues for it receiving a large icy impactor, about 10^–6 Earth masses from far out, requiring some improbable dynamics to get a low enough approach velocity. Core formation in both planets probably started rather early during accretion.Some geochemical evidences argue for the Moon coming from the Earth's mantle, but are inconclusive. Large scale melting of the mantle by the giant impact would plausibly have led to stratification. But the "lock-up" at the end of turbulent mantle convection is a trade-off between rates: crystallization of constituents of small density difference versus overall freezing. Also, factors such as differences in melting temperatures and densities, melt compressibilities, and phase transitions may have had homogenizing effects in the subsequent mantle convection.     

Numerical modelling of heating in porous planetesimal collisions

Collisions between planetesimals at speeds of several kilometres per second were common during the early evolution of our Solar System. However, the collateral effects of these collisions are not well understood. In this paper, we quantify the efficiency of heating during high-velocity collisions between planetesimals using hydrocode modelling. We conducted a series of simulations to test the effect on shock heating of the initial porosity and temperature of the planetesimals, the relative velocity of the collision and the relative size of the two colliding bodies. Our results show that while heating is minor in collisions between non-porous planetesimals at impact velocities below 10 km s⁻¹, in agreement with previous work, much higher temperatures are reached in collisions between porous planetesimals. For example, collisions between nearly equal-sized, porous planetesimals can melt all, or nearly all, of the mass of the bodies at collision velocities below 7 km s⁻¹. For collisions of small bodies into larger ones, such as those with an impactor-to-target mass ratio below 0.1, significant localised heating occurs in the target body. At impact velocities as low as 5 km s⁻¹, the mass of melt will be nearly double the mass of the impactor, and the mass of material shock heated by 100 K will be nearly 10 times the mass of the impactor. We present a first-order estimate of the cumulative effects of impact heating on a porous planetesimal parent body by simulating the impact of a population of small bodies until a disruptive event occurs. Before disruption, impact heating is volumetrically minor and highly localised; in no case was more than about 3% of the parent body heated by more than 100 K. However, heating during the final disruptive collision can be significant; in about 10% of cases, almost all of the parent body is heated to 700 K (from an initial temperature of ∼300 K) and more than a tenth of the parent body mass is melted. Hence, energetic collisions between planetesimals could have had important effects on the thermal evolution of primitive materials in the early Solar System.     

Organic-chemical clues to the theory of impacts as a cause of mass extinctions

The reasons for the mass extinctions, which occur from time to time in Earth's history-as, e.g., the dinosaur extinction at the Cretaceous/Tertiary boundary 65 myr ago - are still not satisfactorily cleared up. A possible reason might be the impact of one or several comets of several kilometers in diameter. In this paper the astrophysical background of this hypothesis and organic-chemical processes during an impact will be discussed.Quantitative estimations are given, which show that the amount of organic substances brought to the Earth may be of the same order of magnitude as the normal biological production of organic material.Investigations are proposed to examine the organic-chemical composition of profiles of the Cretaceous/Tertiary boundary and other boundaries, at which mass extinction had occurred, in order to find anomalies as consequences of impacts.     

Size-frequency distributions of fragments from SPH/N-body simulations of asteroid impacts: Comparison with observed asteroid families

We investigate the morphology of size-frequency distributions (SFDs) resulting from impacts into 100-km-diameter parent asteroids, represented by a suite of 161 SPH/N-body simulations conducted to study asteroid satellite formation [Durda, D.D., Bottke, W.F., Enke, B.L., Merline, W.J., Asphaug, E., Richardson, D.C., Leinhardt, Z.M., 2004. Icarus 170, 243-257]. The spherical basalt projectiles range in diameter from 10 to 46 km (in equally spaced mass increments in logarithmic space, covering six discrete sizes), impact speeds range from 2.5 to 7 km/s (generally in 1 km/s increments), and impact angles range from 15° to 75° (nearly head-on to very oblique) in 15° increments. These modeled SFD morphologies match very well the observed SFDs of many known asteroid families. We use these modeled SFDs to scale to targets both larger and smaller than 100 km in order to gain insights into the circumstances of the impacts that formed these families. Some discrepancies occur for families with parent bodies smaller than a few tens of kilometers in diameter (e.g., 832 Karin), however, so due caution should be used in applying our results to such small families. We find that ∼20 observed main-belt asteroid families are produced by the catastrophic disruption of D>100 km parent bodies. Using these data as constraints, collisional modeling work [Bottke Jr., W.F., Durda, D.D., Nesvorný, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H.F., 2005b. Icarus 179, 63-94] suggests that the threshold specific energy, , needed to eject 50% of the target body's mass is very close to that predicted by Benz and Asphaug [Benz, W., Asphaug, E., 1999. Icarus 142, 5-20].     

Linking the collisional history of the main asteroid belt to its dynamical excitation and depletion

The main belt is believed to have originally contained an Earth mass or more of material, enough to allow the asteroids to accrete on relatively short timescales. The present-day main belt, however, only contains ∼5×¹⁰⁻⁴ Earth masses. Numerical simulations suggest that this mass loss can be explained by the dynamical depletion of main belt material via gravitational perturbations from planetary embryos and a newly-formed Jupiter. To explore this scenario, we combined dynamical results from Petit et al. [Petit, J. Morbidelli, A., Chambers, J., 2001. The primordial excitation and clearing of the asteroid belt. Icarus 153, 338-347] with a collisional evolution code capable of tracking how the main belt undergoes comminution and dynamical depletion over 4.6 Gyr [Bottke, W.F., Durda, D., Nesvorny, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H., 2005. The fossilized size distribution of the main asteroid belt. Icarus 175, 111-140]. Our results were constrained by the main belt's size-frequency distribution, the number of asteroid families produced by disruption events from diameter D>100 km parent bodies over the last 3-4 Gyr, the presence of a single large impact crater on Vesta's intact basaltic crust, and the relatively constant lunar and terrestrial impactor flux over the last 3 Gyr. We used our model to set limits on the initial size of the main belt as well as Jupiter's formation time. We find the most likely formation time for Jupiter was 3.3±2.6 Myr after the onset of fragmentation in the main belt. These results are consistent with the estimated mean disk lifetime of 3 Myr predicted by Haisch et al. [Haisch, K.E., Lada, E.A., Lada, C.J., 2001. Disk frequencies and lifetimes in young clusters. Astrophys. J. 553, L153-L156]. The post-accretion main belt population, in the form of diameter D?1000 km planetesimals, was likely to have been 160±40 times the current main belt's mass. This corresponds to 0.06-0.1 Earth masses, only a small fraction of the total mass thought to have existed in the main belt zone during planet formation. The remaining mass was most likely taken up by planetary embryos formed in the same region. Our results suggest that numerous D>200 km planetesimals disrupted early in Solar System history, but only a small fraction of their fragments survived the dynamical depletion event described above. We believe this may explain the limited presence of iron-rich M-type, olivine-rich A-type, and non-Vesta V-type asteroids in the main belt today. The collisional lifetimes determined for main belt asteroids agree with the cosmic ray exposure ages of stony meteorites and are consistent with the limited collisional evolution detected among large Koronis family members. Using the same model, we investigated the near-Earth object (NEO) population. We show the shape of the NEO size distribution is a reflection of the main belt population, with main belt asteroids driven to resonances by Yarkovsky thermal forces. We used our model of the NEO population over the last 3 Gyr, which is consistent with the current population determined by telescopic and satellite data, to explore whether the majority of small craters (D<0.1-1 km) formed on Mercury, the Moon, and Mars were produced by primary impacts or by secondary impacts generated by ejecta from large craters. Our results suggest that most small craters formed on these worlds were a by-product of secondary rather than primary impacts.     

Radial migration of planetesimals

Planetesimals orbiting a protostar in a circumstellar disk are affected by gravitational interaction among themselves and by gas drag force due to disk gas. Within the Kyoto model of planetesimal accretion, the migration rate is interpreted as the inverse of the planetary formation time scale. Here, we study time scales of gravitational interaction and gas drag force and their influence on planetesimal migration in detail. Evaluating observations of 86 T Tauri stars (Beckwithet al., 1990), we find the mean radial temperature profile of circumstellar disks. The disk mass is taken to be 0.01M in accordance with minimum mass models and observed T Tauri disks. The time scale of gravitational interaction between planetesimals is studied analogously to Chandrasekhar's stellar dynamics. Hence, Chandrasekhar's coefficient , defined as the fraction between the mean separation of planetesimals and the impact parameter, plays an important role in determining the migration rate. We find ln to lie between 5 and 10 within the protosolar disk. Our result is that, at the stage of disk evolution considered here, gas drag force affects the radial migration of planetesimals by a few orders of magnitude more than gravitational interaction.Paper presented at the Conference on Planetary Systems: Formation, Evolution, and Detection held 7–10 December, 1992 at CalTech, Pasadena, California, U.S.A.     

On the accretion of the earth and moon

The bodies which produced the premare impact craters on the moon contained a much higher proportion of smaller bodies in the earliest observable times than subsequently. This suggests that the earth and moon accreted from small objects with only an occasional large planetoid.If the earliest observable lunar craters are 4.3 × 10⁹ yr old, the half-life of the primitive planetesimals which produced the giant lunar craters larger than 161 km in diameter, was 143 × 10⁶ yr, while the half-life of the primitive planetesimals which produced lunar craters larger than 1 km in diameter was only 88 × 10⁶ yr. The half-life of the bodies which produced 1 km craters was still shorter, about 75 × 10⁶ yr.     

Against all odds? Forming the planet of the HD 196885 binary

HD 196885 Ab is the most ??extreme?? planet-in-a-binary discovered to date, whose orbit places it at the limit for orbital stability. The presence of a planet in such a highly perturbed region poses a clear challenge to planet-formation scenarios. We investigate this issue by focusing on the planet-formation stage that is arguably the most sensitive to binary perturbations: the mutual accretion of kilometre-sized planetesimals. To this effect we numerically estimate the impact velocities dv amongst a population of circumprimary planetesimals. We find that most of the circumprimary disc is strongly hostile to planetesimal accretion, especially the region around 2.6 AU (the planet??s location) where binary perturbations induce planetesimal-shattering dv of more than 1 kms^?1. Possible solutions to the paradox of having a planet in such accretion-hostile regions are (1) that initial planetesimals were very big, at least 250 km (2) that the binary had an initial orbit at least twice the present one, and was later compacted due to early stellar encounters (3) that planetesimals did not grow by mutual impacts but by sweeping of dust (the ??snowball?? growth mode identified by Xie et al., in Astrophys J 724:1153, 2010b), or (4) that HD 196885 Ab was formed not by core-accretion but by the concurrent disc instability mechanism. All of these 4 scenarios remain however highly conjectural.     

Protoatmospheres and Surface Environment of Protoplanets

Protoatmospheres and surface environment of terrestrial protoplanets during the oligarchic accretion phase and the giant impacts phase are discussed from theoretical points of view. Mars-sized protoplanets form during the stage of the oligarchic growth. Since protoplanets are formed from more or less ‘local’ planetesimals, the surface environment of the accreting protoplanets depends on availability of volatile material in planetesimals. Even if no volatile-bearing planetesimals are available, a gravitationary captured solar composition atmosphere is formed during accretion. In such cases the surface temperature is always kept under the melting temperature of mantle silicate and only a subsurface magma ocean is formed. Core formation proceeds under dry conditions, and volatile elements are not partitioned into metallic iron. Accretion of water-bearing planetesimals results in impact degassing. A surface hydrous magma ocean forms in response to the thermal blanketing effect of the proto-atmosphere. Then, some volatile materials dissolve into the magma ocean. If we consider reaction with metallic iron, the proto-atmosphere is likely to be rich in hydrogen. In addition, a large amount of hydrogen may be partitioned into metallic iron under high pressure, and delivered to the core. In the stage of giant impacts, both dry and water-bearing protoplanets collide on the proto-Earth. Substantial amount of proto-atmosphere (including water vapor) survives giant impacts. Moreover, giant impacts on protoplanets with oceans result in relative concentration of water against other gases.