Why we study the geological record for evidence of the solar oscillation about the galactic midplane

The Solar System oscillates about the plane defined by the disk of matter in our Galaxy. This oscillatory motion gives rise to a substantial modulation in the tidally induced flux of Oort cloud comets. An observational determination of the quasi-periodicity of this motion carries with it significant information about the population, distributions, dynamics and origins of short-period and long-period comets. An additional incentive for emphasizing such a study is the information about dark disk matter that a period can yield. If dark disk matter is completely negligible, the amplitude of the solar motion will be sufficiently large that the peak-to-trough flux ratio will be ≈ 2.5 and the plane-crossing period will exceed 40 Myr. Dark disk matter comparable in mass to bright disk matter and distributed in any manner is inconsistent with K-dwarf distributions and can be rejected as a working hypothesis. But if a modest fraction of the disk matter is dark and distributed like the interstellar medium, as is consistent with limits deduced from K-giant and K-dwarf velocity distributions, the peak-to-trough flux ratio can increase to a factor of 4 even though the solar z amplitude is decreased. In that case the period can be as little as 30 Myr and the implied Oort population is smaller by a factor of 3. We should carefully reconsider the geological record as a potential discriminator of these options.     

The Galactic Theory of Mass Extinctions: an Update

Astronomical and geological evidence is consistent with the hypothesis that mass extinctions of life on Earth are related to impacts of comets whose flux is partly modulated by the dynamics of the Milky Way Galaxy. Geologic evidence for impact (ejecta and large impact craters) has been found at times of mass extinction events, and the record of large dated craters shows a significant correlation with extinctions. Statistical analyses suggest that mass extinction events exhibit a periodic component of about 30 Myr, and periodicities of 30 ± 0.5 Myr and 35 ± 2 Myr have been extracted from sets of well-dated large impact craters. These results suggest periodic or quasi-periodic showers of impactors, probably Oort Cloud comets, with an approximately 30 or 36 Myr cycle. The best explanation for these proposed quasi-periodic comet showers involves the Sun's vertical oscillation through the galactic disk, which may have a similar cycle time between crossings of the galactic plane. Further refinement of the model will depend on the identification and quantification of the dark matter component in the galactic disk, and discovery and accurate dating of additional impact craters. This revised version was published online in July 2006 with corrections to the Cover Date.     

NEOs and Impacts: The Galactic Connection

The orbit of the solar system within the Galaxy is subject both to adiabatic tidal forces and to irregular buffeting caused by encounters with molecular clouds and spiral arms, while the long-period comet flux varies with the differential of these forces. Quasi-periodic episodes of Earth bombardment by Oort cloud comets are expected, with the half-period of the solar vertical oscillations being prominent. Various terrestrial datasets are examined and show evidence of a periodicity ∼27±1 Myr at a confidence level ≳ 99.5%, consistently with the expectations. Bombardment episodes thus appear to exert a controlling influence on global terrestrial processes. Conversely, the long-term terrestrial record provides an empirical constraint on the source regions of NEOs as well as on the dark matter content of the Galactic disc. During such episodes, climatic effects due to stratospheric dusting by very large comets are likely to be a major factor in geological and biological change. This revised version was published online in July 2006 with corrections to the Cover Date.     

A note on the total mass of comets in the solar system

The assumption that the very low albedo determined for Halley's comet is typical of all short period comets, taken together with the assumption that the average sizes of long and short period comets are approximately equal, leads to an increase in the total mass of comets in the solar system by almost two orders of magnitude. If gravitational ejection from the Uranus - Neptune zone during the later phases of planet formation is indeed responsible for the classical Oort cloud between 10⁴–10¹⁵ AU, then the mass of comets in this transplanetary region during cosmogonie times has to exceed the combined masses of Uranus and Neptune by over an order of magnitude. Furthermore, if the recent arguments for as many as 10¹⁴ comets in an inner Oort cloud between ~40– 10⁴AU are valid, then the total mass of comets in the solar system approaches 2% of a solar mass.     

Long-term effects of the Galactic tide on cometary dynamics

We introduce a model for integrating the effects of Galactic tides on Oort cloud comets, which involves two procedures, according to the values of the osculating semi-major axis a and eccentricity e. Ten simulations of the dynamics of 10⁶ comets over 5 Gyr are performed using this model. We thus investigate the long-term effects of the Galactic tide with and without a radial component, the effects of the local density of the Galactic disk, and those of the Oort constants. Most of the results may be understood in terms of the integrability or non-integrability of the system. For an integrable system, which occurs for moderate semi-major axes with or without radial component, the dynamics is explained by periodic variation of the cometary perihelion, inducing the depletion of the outer region of the Oort cloud, a constant flux from the inner region after 500 Myr, and the quick formation of a reservoir of comets with argument of perihelion near 26.6°. When the system is non-integrable, the efficiency of the tide in reducing the cometary perihelion distance is enhanced both by replenishing the Oort cloud domain from which comets are sent toward the planetary system, and by reducing the minimal value that the perihelion distance may reach. No effects of varying the Oort constants were observed, showing that the flat rotation curve is a satisfactory approximation in Oort cloud dynamics.     

Long-term evolution of Oort Cloud comets: capture of comets

We test different possibilities for the origin of short-period comets captured from the Oort Cloud. We use an efficient Monte Carlo simulation method that takes into account non-gravitational forces, Galactic perturbations, observational selection effects, physical evolution and tidal splittings of comets. We confirm previous results and conclude that the Jupiter family comets cannot originate in the spherically distributed Oort Cloud, since there is no physically possible model of how these comets can be captured from the Oort Cloud flux and produce the observed inclination and Tisserand constant distributions. The extended model of the Oort Cloud predicted by the planetesimal theory consisting of a non-randomly distributed inner core and a classical Oort Cloud also cannot explain the observed distributions of Jupiter family comets. The number of comets captured from the outer region of the Solar system are too high compared with the observations if the inclination distribution of Jupiter family comets is matched with the observed distribution. It is very likely that the Halley-type comets are captured mainly from the classical Oort Cloud, since the distributions in inclination and Tisserand value can be fitted to the observed distributions with very high confidence. Also the expected number of comets is in agreement with the observations when physical evolution of the comets is included. However, the solution is not unique, and other more complicated models can also explain the observed properties of Halley-type comets. The existence of Jupiter family comets can be explained only if they are captured from the extended disc of comets with semimajor axes of the comets   a <5000 au  . The original flattened distribution of comets is conserved as the cometary orbits evolve from the outer Solar system era to the observed region.     

Impact cratering and the Oort Cloud

We calculate the expected flux profile of comets into the planetary system from the Oort Cloud arising from Galactic tides and encounters with molecular clouds. We find that both periodic and sporadic bombardment episodes, with amplitudes an order of magnitude above background, occur on characteristic time-scales ∼25–35 Myr. Bombardment episodes occurring preferentially during spiral arm crossings may be responsible both for mass extinctions of life and the transfer of viable microorganisms from the bombarded Earth into the disturbing nebulae. Good agreement is found between the theoretical expectations and the age distribution of large, well-dated terrestrial impact craters of the past 250 Myr. A weak periodicity of ∼36 Myr in the cratering record is consistent with the Sun's recent passage through the Galactic plane, and implies a central plane density  ∼0.15 M_⊙ pc⁻³  . This leaves little room for a significant dark matter component in the disc.     

Near-parabolic cometary flux in the outer solar system

The flux of near-parabolic comets in the outer-planetary region is estimated on the presumption that the major planets and the galactic tide control the dynamics of comets. It is found that the flux of the Oort cloud comets (semi-major axis > 20000 AU) is similar to the case of a strong comet shower derived on the presumption that the galactic tidal force were not operative. On the other hand, the flux of comets with semi-major axes <- 20000 AU is found to be an increasing function of q (perihelion distance) until q reaches 20 AU, while for a 45000 AU it is a rapidly increasing function for q 12 AU. In other words, for comets of the inner extension of the Oort cloud the planetary perturbation acts as a strong barrier for them to penetrate into the inner planetary region.     

Origin of oort cloud comets in the interstellar space

Comets must form a major part of the interstellar medium. The solar system provides a flux of comets into the interstellar space and there is no reason to suspect that many other stars and their surrounding cometary systems would not make a similar contribution. Occasionally interstellar comets must pass through the inner solar system, but Whipple (1975) considers it unlikely that such a comet is among the known cases of apparently hyperbolic comets. Even so the upper limit for the density of unobserved interstellar comets is relatively high.In addition, we must consider the possibility that comets are a genuine component of interstellar medium, and that the Oort Cloud is merely a captured part of it (McCrea, 1975). Here we review various dynamical possibilities of two-way exchange of comet populations between the Solar System and the interstellar medium. We describe ways in which a traditional Oort Cloud (Oort, 1950) could be captured from the interstellar medium. However, we note that the so called Kuiper belt (Kuiper, 1951) of comets cannot arise through this process. Therefore we have to ask how necessary the concept of the yet unobserved Kuiper belt is for the theory of short period comets.There has been considerable debate about the question whether short period comets can be understood as a captured population of the Oort Cloud of comets or whether an additional source has to be postulated. The problem is made difficult by the long integration times of comet orbits through the age of the Solar System. It would be better to have an accurate treatment of comet-planet encounters in a statistical sense, in the form of cross sections, and to carry out Monte Carlo studies. Here we describe the plan of action and initial results of the work to derive cross sections by carrying out large numbers of comet — planet encounters and by deriving approximate analytic expressions for them. Initially comets follow parabolic orbits of arbitrary inclination and perihelion distance; cross sections are derived for obtaining orbits of given energy and inclination after the encounter. The results are used in subsequent work to make evolutionary models of the comet population.     

Long-Period Comets and the Oort Cloud

Long-period (LP) comets, Halley-type (HT) comets, and even some comets of the Jupiter family, probably come from the Oort cloud, a huge reservoir of icy bodies that surrounds the solar system. Therefore, these comets become important probes to learn about the distant Oort cloud population. We review the fundamental dynamical properties of LP comets, and what is our current understanding of the dynamical mechanisms that bring these bodies from the distant Oort cloud region to the inner planetary region. Most new comets have original reciprocal semimajor axes in the range2 × 10^-5 < 1/a_orig < 5 × 10^-5AU^-1. Yet, this cannot be taken to represent the actual space distribution of Oort cloud comets, but only the region in the energy space in which external perturbers have the greatest efficiency in bringing comets to the inner planetary region. The flux of Oort cloud comets in the outer planetary region is found to be at least several tens times greater than the flux in the inner planetary region. The sharp decrease closer to the Sun is due to the powerful gravitational fields of Jupiter and Saturn that prevent most Oort cloud comets from reaching the Earth’s neighborhood (they act as a dynamical barrier). A small fraction of ～10^-2 Oort cloud comets become Halley type (orbital periods P < 200 yr), and some of them can reach short-period orbits with P < 20 yr. We analyze whether we can distinguish the latter, very ‘old” LP comets, from comets of the Jupier family coming from the Edgeworth-Kuiper belt.     

Embedded star clusters and the formation of the Oort cloud: II. The effect of the primordial solar nebula

This paper deals with Oort cloud formation while the Sun was in an embedded cluster and surrounded by its primordial nebula. This work is a continuation of Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82], building on the model presented therein, and adding the aerodynamic drag and gravitational potential of the primordial solar nebula. Results are presented of numerical simulations of comets subject to the gravitational influence of the Sun, Jupiter, Saturn, star cluster and primordial solar nebula; some of the simulations included the gravitational influence of Uranus and Neptune as well. The primordial solar nebula was approximated by the minimum-mass Hayashi model [Hayashi, C., Nakozawa, K., Nakagawa, Y., 1985. In: Black, D.C., Matthews, M.S. (Eds.). Protostars and Planets II. Univ. of Arizona Press, Tucson, AZ] whose inner and outer radii have been truncated at various distances from the Sun. A comet size of 1.7 km was used for most of our simulations. In all of our simulations, the density of the primordial solar nebula decayed exponentially with an e-folding time of 2 Myr. It turns out that when the primordial solar nebula extends much beyond Saturn or Neptune, virtually no material will end up in the Oort cloud (OC) during this phase. Instead, the majority of the material will be on circular orbits inside of Jupiter if the inner edge of the disk is well inside Jupiter's orbit. If the disk's inner edge is beyond Jupiter's orbit, most comets end up on orbits in exterior mean-motion resonances with Saturn when Uranus and Neptune are not present. In those cases where the outer edge of the disk is close to Saturn or Neptune, the fraction of material that ends up in the subsequently formed OC is much less than that found in Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82] for the same cluster densities. This implies that for comets of roughly 2 km in size, the presence of the primordial solar nebula hinders OC formation. A byproduct of some of our simulations are endresults with a substantial fraction of the comets in the Uranus-Neptune scattered disk. A subsequent followup of this material is planned for the near future. In order to determine the effect of the size of the comets on OC formation efficiency, a set of runs with the same initial conditions but different cometary radii have been performed as well, from which it is determined that the threshold comet size to begin producing significant Oort clouds is roughly 20 km. This implies that the presence of the primordial solar nebula acts as a size-sorting mechanism, with large bodies unaffected by the gas drag and ending up in the OC while small bodies remain trapped in the planetary region, in the models studied.     

Peripheral Structures of Planetary Systems

We analyze the conditions for the formation and time evolution of peripheral comet structures of solar-type planetary systems. In the Solar system, these include the Kuiper belt, the Oort cloud, the comet spear, and the Galactic comet ring that marks the Galactic orbit of the Sun. We consider the role of the viscosity of a protoplanetary gas–dust disk, major planets, field stars, globular clusters, giant molecular clouds, and the Galactic gravitational field in the formation of these peripheral structures marked by comets and asteroids. We give a list of the closest past and future passages of neighboring stars through the solar Oort cloud that perturb the motion of its comets and, thus, contribute to the enhancement of its cometary activity, on the one hand, and to the replenishment of the solar comet spear with new members, on the other hand.     

Periodic variation of Oort Cloud flux and cometary impacts on the Earth and Jupiter

Time variation in impact probability is studied by assuming that the periodic flux of the Oort Cloud comets within 15 au arises from the motion of the Sun with respect to the Galactic mid-plane. The periodic flux clearly shows up in the impact rate of the captured Oort Cloud cometary population, with a phase shift caused by the orbital evolution. Depending on the assumed flux of comets and the size distribution of comets, the impact rate of the Oort Cloud comets of 1 km in diameter or greater is from 5 to 700 impacts Myr⁻¹ on the Earth and from 0.5 to 70 impacts per 1000 yr on Jupiter. The relative fractions of impacts are 0.09, 0.11, 0.26 and 0.54 for long-period comets, Halley type comets, Jupiter family comets and near-Earth objects, respectively. For Jupiter, the corresponding fractions in the first three categories are 0.18, 0.31 and 0.51. If we consider physical fading of comet activity that is compatible with the observations, then the impact rates of active comets are two orders of magnitude smaller than the total impact rates by all kinds of comets and cometary asteroids of size 1 km or greater.     

The key role of massive stars in Oort cloud comet dynamics

The effects of a sample of 1300 individual stellar encounters spanning a wide range of parameter values (mass, velocity and encounter distance) are investigated. Power law fits for the number of injected comets demonstrate the long range effect of massive stars, whereas light stars affect comets mainly along their tracks. Similarly, we show that the efficiency of a star to fill the phase space region of the Oort cloud where the Galactic tides are able to inject comets into the observable region - the so-called “tidally active zone” (TAZ) - is also strongly dependent on the stellar mass. Power laws similar to those for direct injection are obtained for the efficiency of stars to fill the TAZ. This filling of the tidally active zone is crucial for the long term flux of comets from the Oort cloud. Based on long-term Monte Carlo simulations using a constant Galactic tide and a constant flux of stellar encounters, but neglecting the detailed effects of planetary perturbations, we show that this flux essentially results from a two step mechanism: (i) the stellar injection of comets into the TAZ; and (ii) the tidal injection of TAZ comets into the loss cone. We find that single massive stars are able to induce “comet drizzles” - corresponding to an increase of the cometary flux of about 40% - which may last for more than 100 Myr by filling the TAZ to a higher degree than normal. It appears that the stars involved in this process are the same that cause comet showers.     

Long-period comet flux in the planetary region: Dynamical evolution from the Oort cloud

This study analyzes the evolution of 2 × 10⁵ orbits with initial parameters corresponding to the orbits of comets of the Oort cloud under the action of planetary, galactic, and stellar perturbations over 2 × 10⁹ years. The dynamical evolution of comets of the outer (orbital semimajor axes a > 10⁴ AU) and inner (5 × 10³ < a (AU) < 10⁴) parts of the comet cloud is analyzed separately. The estimates of the flux of “new” and long-period comets for all perihelion distances q in the planetary region are reported. The flux of comets with a > 10⁴ AU in the interval 15 AU < q < 31 AU is several times higher than the flux of comets in the region q < 15 AU. We point out the increased concentration of the perihelia of orbits of comets from the outer cloud, which have passed several times through the planetary system, in the Saturn-Uranus region. The maxima in the distribution of the perihelia of the orbits of comets of the inner Oort cloud are located in the Uranus-Neptune region. “New” comets moving in orbits with a < 2 × 10⁴ AU and arriving at the outside of the planetary system (q > 25 AU) subsequently have a greater number of returns to the region q < 35 AU. The perihelia of the orbits of these comets gradually drift toward the interior of the Solar System and accumulate beyond the orbit of Saturn. The distribution of the perihelia of long-period comets beyond the orbit of Saturn exhibits a peak. We discuss the problem of replenishing the outer Oort cloud by comets from the inner part and their subsequent dynamical evolution. The annual rate of passages of comets of the inner cloud, which replenish the outer cloud, in the region q < 1 AU in orbits with a > 10⁴ AU (～ 5.0 × 10^?14 yr^?1) is one order of magnitude lower than the rate of passage of comets from the outer Oort cloud (～ 9.1 × 10^?13 yr^?1).     

The Scattered Disk Population and the Oort Cloud

The trans-Neptunian belt has been subject to a strong depletion that has reduced its primordial population by a factor of one hundred over the solar system's age. One by-product of such a depletion process is the existence of a scattered disk population in transit from the belt to other places, such as the Jupiter zone, the Oort cloud or interstellar space. We have integrated the orbits of the scattered disk objects (SDOs) so far discovered by 2500 Myr to study their dynamical time scales and the probability of falling in each of the end states mentioned above, paying special attention to their contribution to the Oort cloud. We found that their dynamical half-time is close to 2.5 Gyr and that about one third of the SDOs end up in the Oort cloud.     

The Role of Giant Molecular Clouds in the Evolution of the Oort Comet Cloud

We estimated the gravitational influence of giant molecular clouds passing near the Solar system on the orbital evolution of Oort cloud comets. We performed a comparative analysis of the accuracies of the following two methods of allowance for the perturbations from giant molecular clouds: the impulse approximation and numerical integration. The impulse approximation yields fairly accurate estimates of the change in the energy of Oort cloud comets and the probability of their ejection under the influence of a molecular cloud if the path of the Solar system does not cross its boundary and if the molecular cloud may be treated as a point perturbing mass. The comet survival probability in the Oort cloud depends significantly on the internal structure of the perturbing molecular cloud and the impact parameter of the encounter. The most massive injection of comets into the planetary region and their ejection from the Oort cloud take place if the Solar system passes through a giant molecular cloud composed of several high-mass condensations. In this case, most of the comets injected into the planetary region were initially comets of the inner Oort cloud (a 10^–4 AU) with high orbital eccentricities.     

Tidal gravitational forces: The infall of “new” comets and comet showers

The tidal gravitational field of the Galaxy directed into the galactic plane changes the angular momentum of comets in the Oort cloud. For comet orbits with semimajor axis greater than 2 × 10⁴ AU, the change of angular momentum in one orbit is sufficient to bring comets from the Oort cloud into the visible region, causing the infall of “new” comets. The limiting size orbit is weakly dependent on the angle between the major axis of the comet orbit and the galactic plane. The flux of comets into the inner Solar System caused by the galactic tidal field will be continuous and nearly isotropic. This effect appears to exclude any determination of the trajectories of passing stars by analysis of the angular distribution of new comets. The production of intense comet showers by the tidal field of a solar companion or of an interstellar cloud is considered. We show that the direction of a solar companion cannot be found from the present distribution of observable comets. The frequency of comet showers induced by encounters with interstellar clouds is found to be much lower than that from passing stars, and the tidal fields of interstellar clouds are not strong enough to cause comet showers of sufficient intensity to result in Earth impacts.     

Orbital element distributions in the Oort cloud

In this paper we consider orbital element distributions for comets moving on admissible orbits in the Oort cloud and distributions for some functions that depend on the orbital elements. Moreover, we find the probability of an event that an arbitrarily chosen admissible orbit belongs to the set (r) of orbital elements and the distribution of circular velocities in the cloud.     

Implications of the observed distributions of very long period comet orbits

With allowance for galactic perturbations and observational error, the observed distributions of sizes an orientations of very long period comets are consistent with a uniform distribution of comets within the Oort Cloud.