Water and related chemistry in the solar system. A guaranteed time key programme for Herschel

“Water and related chemistry in the Solar System” is a Herschel Space Observatory Guaranteed-Time Key Programme. This project, approved by the European Space Agency, aims at determining the distribution, the evolution and the origin of water in Mars, the outer planets, Titan, Enceladus and the comets. It addresses the broad topic of water and its isotopologues in planetary and cometary atmospheres. The nature of cometary activity and the thermodynamics of cometary comae will be investigated by studying water excitation in a sample of comets. The D/H ratio, the key parameter for constraining the origin and evolution of Solar System species, will be measured for the first time in a Jupiter-family comet. A comparison with existing and new measurements of D/H in Oort-cloud comets will constrain the composition of pre-solar cometary grains and possibly the dynamics of the protosolar nebula. New measurements of D/H in giant planets, similarly constraining the composition of proto-planetary ices, will be obtained. The D/H and other isotopic ratios, diagnostic of Mars’ atmosphere evolution, will be accurately measured in H₂O and CO. The role of water vapor in Mars’ atmospheric chemistry will be studied by monitoring vertical profiles of H₂O and HDO and by searching for several other species (and CO and H₂O isotopes). A detailed study of the source of water in the upper atmosphere of the Giant Planets and Titan will be performed. By monitoring the water abundance, vertical profile, and input fluxes in the various objects, and when possible with the help of mapping observations, we will discriminate between the possible sources of water in the outer planets (interplanetary dust particles, cometary impacts, and local sources). In addition to these inter-connected objectives, serendipitous searches will enhance our knowledge of the composition of planetary and cometary atmospheres.     

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.     

The scattered disk as a source of Halley-type comets

We investigate a new dynamical mechanism for producing Halley-type comets from the scattered disk of comets. Levison and Duncan [Levison, H., Duncan, M., 1997. Icarus 127, 13-32] and Duncan and Levison [Duncan, M., Levison, H., 1997. Science 276, 1670-1672] showed that a significant number of objects leave the scattered disk by evolving to semi-major axes greater than 1000 AU. We find that once these objects reach semi-major axes on the order of ¹⁰⁴ AU, a significant fraction immediately have their perihelia driven inward by the galactic tides. Approximately 0.01% of the objects that reach ¹⁰⁴ AU then evolve onto orbits similar to the observed Halley-like comets due to gravitational interactions with the giant planets. The orbital element distribution resulting from this process is statistically consistent with observations. We discuss the implications of this model for the number of objects in the scattered disk in the text. The model predicts a temporal variation in the influx of HTCs with a period of ∼118 Myr. At the peak, the model predicts that there should be roughly 10 times as many HTCs as currently observed (i.e., there should be weak HTC showers). However, the model may inflate the importance of these showers because it does not include the effects of passing stars and giant molecular clouds.     

Extra-solar Oort cloud encounters and planetary impact rates

If other stars possess Oort-like comet clouds, then some such clouds will pass sufficiently close to our Sun to induce an influx of Extra-Solar comets through the planetary region. We investigate this source, and find that while the expected number of planetary impacts due to Extra-Solar Oort comets will dominate impacts caused by free interstellar comets, only a few such comets have impacted the terrestial planets over time.     

Ultraviolet absorption studies of H2O and other species in comet comae with satellite telescope-spectrometers

Ultraviolet absorption by H₂O and other species in the comae of comets could be detected by studying, with satellite telescope-spectrometers, the occultation of hot stars by comets. Such observations could produce the first direct detection of H₂O, the fundamental parent molecule in comet comae, and give measures of molecular level populations. The first instrument suitable for such observations will be the High Resolution Spectrograph on Space Telescope and, therefore, we consider its capabilities. We have used a Haser model to estimate the molecular column densities and to predict equivalent widths for lines of H₂O, OH, CO, and O as functions of time and angular distance from a comet with a high H₂O production rate. We have determined the minimum detectable equivalent widths, and therefore, the maximum angular separation from such a comet at which H₂O, OH, and CO could be studied. A conservative, statistical estimate shows that comets with high water production rates should pass near enough to about 10 to 100 stars suitable for absorption studies of the $\text{C}\text{}\text{X}$ band of H₂O (1240 Å). Estimated equivalent widths for CO, OH, and the resonance lines of C and O indicate that these species may also be detected.     

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.     

Habitable zones around main sequence stars

A one-dimensional climate model is used to estimate the width of the habitable zone (HZ) around our Sun and around other main sequence stars. Our basic premise is that we are dealing with Earth-like planets with CO2/H2O/N2 atmospheres and that habitability requires the presence of liquid water on the planet's surface. The inner edge of the HZ is determined in our model by loss of water via photolysis and hydrogen escape. The outer edge of the HZ is determined by the formation of CO2 clouds, which cool a planet's surface by increasing its albedo and by lowering the convective lapse rate. Conservative estimates for these distances in our own Solar System are 0.95 and 1.37 AU, respectively; the actual width of the present HZ could be much greater. Between these two limits, climate stability is ensured by a feedback mechanism in which atmospheric CO2 concentrations vary inversely with planetary surface temperature. The width of the HZ is slightly greater for planets that are larger than Earth and for planets which have higher N2 partial pressures. The HZ evolves outward in time because the Sun increases in luminosity as it ages. A conservative estimate for the width of the 4.6-Gyr continuously habitable zone (CHZ) is 0.95 to 1.15 AU. Stars later than F0 have main sequence lifetimes exceeding 2 Gyr and, so, are also potential candidates for harboring habitable planets. The HZ around an F star is larger and occurs farther out than for our Sun; the HZ around K and M stars is smaller and occurs farther in. Nevertheless, the widths of all of these HZs are approximately the same if distance is expressed on a logarithmic scale. A log distance scale is probably the appropriate scale for this problem because the planets in our own Solar System are spaced logarithmically and because the distance at which another star would be expected to form planets should be related to the star's mass. The width of the CHZ around other stars depends on the time that a planet is required to remain habitable and on whether a planet that is initially frozen can be thawed by modest increases in stellar luminosity. For a specified period of habitability, CHZs around K and M stars are wider (in log distance) than for our Sun because these stars evolve more slowly. Planets orbiting late K stars and M stars may not be habitable, however, b ecause they can become trapped in synchronous rotation as a consequence of tidal damping. F stars have narrower (log distance) CHZ's than our Sun because they evolve more rapidly. Our results suggest that mid-to-early K stars should be considered along with G stars as optimal candidates in the search for extraterrestrial life.     

A Postulated Planetary Collision, the Terrestrial Planets, the Moon and Smaller Solar-System Bodies

In a scenario produced by the Capture Theory of planetary formation, a collision between erstwhile solar-system giant planets, of masses 798.75 and 598.37 M _⊕, is simulated using smoothed-particle hydrodynamics. Due to grain-surface chemistry that takes place in star-forming clouds, molecular species containing hydrogen, with a high D/H ratio taken as 0.01, form a layer around each planetary core. Temperatures generated by the collision initiate D–D reactions in these layers that, in their turn, trigger a reaction chain involving heavier elements. The nuclear explosion shatters and disperses both planets, leaving iron-plus-silicate stable residues identified as a proto-Venus and proto-Earth. A satellite of one of the colliding planets, captured or retained by the proto-Earth core, gave the Moon; two massive satellites released into heliocentric orbits became Mercury and Mars. For the Moon and Mars, abrasion of their surfaces exposed to collision debris results in hemispherical asymmetry. Mercury, having lost a large part of its mantle due to massive abrasion, reformed to give the present high-density body. Debris from the collision gave rise to asteroids and comets, much of the latter forming an inner reservoir stretching outwards from the inner Kuiper Belt that replenishes the Oort Cloud when it is depleted by a severe perturbation. Other features resulting from the outcome of the planetary collision are the relationship of Pluto and Triton to Neptune, the presence of dwarf planets and light-atom isotopic anomalies in meteorites.     

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.     

Relationship of comets and planets

We analyze our earlier data on the numerical integration of the equations of motion for 274 short-period comets (with the period P<200 yr) on a time interval of 6000 yr. As many as 54 comets had no close approaches to planets, 13 comets passed through the Saturnian sphere of action, and one comet passed through the Uranian sphere of action. The orbital elements of these 68 comets changed by no more than ±3 percent in a space of 6000 yr. As many as 206 comets passed close to Jupiter. We confirm Everhart’s conclusion that Jupiter can capture long-period comets with q = 4–6 AU and i < 9° into short-period orbits. We show that nearly parabolic comets cross the solar system mainly in the zone of terrestrial planets. No relationship of nearly parabolic comets and terrestrial planets was found for the epoch of the latest apparition of comets. Guliev’s conjecture about two trans-Plutonian planets is based on the illusory excess of cometary nodes at large heliocentric distances. The existence of cometary nodes at the solar system periphery turns out to be a solely geometrical effect.     

江苏科技信息网之天文数据库

中国科学院紫金山天文台受江苏省科工委、江苏省科技信息中心委托研制的天文数据库,是省科技信息网的重要组成部分。该信息网包括工业、农业、土壤、地理等较为全面的科技信息,将于年底全面开通。天文数据库已于十月中旬上网。它包括了阴阳历对照、行星特殊天象、太阳、恒星、人造地球卫星等方面的内容。用ＦｏｘＰｒｏ５ ．０ 建立的天文数据库管理系统,包括了成百个数据库表,数据量５ 兆左右,记录约为５ 万余条。天体力学部和天体物理部等有关部门,在提供数据方面,给予了大力支持,在此一并表示感谢。     

Theoretical predictions of deuterium abundances in the Jovian planets

Using current concepts for the origin of the Jovian planets and current constraints on their interior structure, we argue that the presence of large amounts of “ice” (H₂O, CH₄, and NH₃) in Uranus and Neptune indicates temperatures low enough to condense these species at the time Uranus and Neptune formed. Yet such low temperatures imply orders-of-magnetude fractionation effects for deuterium into the “ice” component if isotopic equilibration can occur. Our models thus imply that Uranus and Neptune should have a D/H ratio at least four times primordial, contrary to observation for Uranus. We find that the Jovian and Saturnian D/H should be close to primordial regardless of formation scenario. The Uranus anomaly could indicate that there was a strong initial radial gradient in D/H in the primordial solar nebula, or that Uranus is so inactive that no significant mixing of its interior has occurred over the age of the solar system. Observation of Neptune's atmospheric D/H may help to resolve the problem.     

Reassessing the formation of the inner Oort cloud in an embedded star cluster

We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser, R., Duncan, M.J., Levison, H.F. [2006]. Icarus 184, 59–82) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar encounters with the outer Solar System as hyperbolic fly-bys. Second, we have included the expulsion of the cluster gas, a feature that was absent previously. Third, we have used several models for the initial conditions and density profile of the cluster – either a Hernquist or Plummer potential – and chose other parameters based on the latest observations of embedded clusters from the literature. These other parameters result in the stars being on radial orbits and the cluster collapses. Similar to previous studies, in our simulations the inner Oort cloud is formed from comets being scattered by Jupiter and Saturn and having their pericentres decoupled from the planets by perturbations from the cluster gas and other stars. We find that all inner Oort clouds formed in these clusters have an inner edge ranging from 100 AU to a few hundred AU, and an outer edge at over 100,000 AU, with little variation in these values for all clusters. All inner Oort clouds formed are consistent with the existence of (90377) Sedna, an inner Oort cloud dwarf planetoid, at the inner edge of the cloud: Sedna tends to be at the innermost 2% for Plummer models, while it is 5% for Hernquist models. We emphasise that the existence of Sedna is a generic outcome. We define a ‘concentration radius’ for the inner Oort cloud and find that its value increases with increasing number of stars in the cluster, ranging from 600 AU to 1500 AU for Hernquist clusters and from 1500 AU to 4000 AU for Plummer clusters. The increasing trend implies that small star clusters form more compact inner Oort clouds than large clusters. We are unable to constrain the number of stars that resided in the cluster since most clusters yield inner Oort clouds that could be compatible with the current structure of the outer Solar System. The typical formation efficiency of the inner Oort cloud is 1.5%, significantly lower than previous estimates. We attribute this to the more violent dynamics that the Sun experiences as it rushes through the centre of the cluster during the latter’s initial phase of violent relaxation.     

Constraints on the Formation Regions of Comets from their D:H Ratios

Studies of the D:H ratio in H₂O within the Solar nebula provide a relationship between the degree of enrichment of deuterium and the distance from the young Sun. In the context of cometary formation, such models suggest that comets which formed in different regions of the Solar nebula should have measurably different D:H ratios. We aim to illustrate how the observed comets can give information about the formation regions of the reservoirs in which they originated. After a discussion of the current understanding of the regions in which comets formed, simple models of plausible formation regions for two different cometary reservoirs (the Edgeworth–Kuiper belt and the Oort Cloud) are convolved with a deuterium-enrichment profile for the pre-solar nebula. This allows us to illustrate how different formation regions for these objects can lead to great variations in the deuterium enrichment distributions that we would observe in comets today. We also provide an illustrative example of how variations in the population within a source region can modify the resulting observational profile. The convolution of a deuterium-enrichment profile with examples of proto-cometary populations gives a feel for how observations could be used to draw conclusions on the formation region of comets which are currently fed into the inner Solar system from at least two reservoirs. Such observations have, to date, been carried out on only three comets, but future work with instruments such as ALMA and Herschel should vastly improve the dataset, leading to a clearer consensus on the formation of the Oort cloud and Edgeworth–Kuiper belt.     

Encounter frequency of Halley-like comets with the planets

We present the results of a numerical study on encounter frequencies of fictitious Halley-like comets with the planets in a dynamical model of the solar system, in which we take into account the gravitational forces of the Sun and the planets Venus through Neptune. The change of the orbital elements during a close approach with a planet was carefully monitored with the aid of a thoroughly tested numerical integration method with automatic step size control. We computed the encounter frequencies of the comets' orbits using two different spheres of influence and compared the results. In both cases, it turned out that the encounter frequency of the fictitious Halley-like comets with Jupiter and Saturn is about a factor 10 to 100 higher than for the other planets. Concerning the changes of the semi-major axes and inclinations our results show that an increase and decrease of these elements is equally probable after an encounter.     

Report of the IAU/IAG Working Group on cartographic coordinates and rotational elements: 2006

Every three years the IAU/IAG Working Group on Cartographic Coordinates and Rotational Elements revises tables giving the directions of the poles of rotation and the prime meridians of the planets, satellites, minor planets, and comets. This report introduces improved values for the pole and rotation rate of Pluto, Charon, and Phoebe, the pole of Jupiter, the sizes and shapes of Saturn satellites and Charon, and the poles, rotation rates, and sizes of some minor planets and comets. A high precision realization for the pole and rotation rate of the Moon is provided. The expression for the Sun’s rotation has been changed to be consistent with the planets and to account for light travel time     

Extrasolar planets and the rotation and axisymmetric mass-loss of evolved stars

I examine the implications of the recently found extrasolar planets on the planet-induced axisymmetric mass-loss model for the formation of elliptical planetary nebulae (PNe). This model attributes the low departure from spherical mass-loss of upper asymptotic giant branch (AGB) stars to envelope rotation which results from deposition of orbital angular momentum of the planets. Since about half of all PNe are elliptical, i.e., have low equatorial to polar density contrast, it was predicted that about 50 per cent of all Sun-like stars have Jupiter-like planets around them, i.e., a mass about equal to that of Jupiter, M _J, or more massive. In the light of the new findings that only 5 per cent of Sun-like stars have such planets, and a newly proposed mechanism for axisymmetric mass-loss, the cool magnetic spots model, I revise this prediction. I predict that indeed ∼50 per cent of PN progenitors do have close planets around them, but the planets can have much lower masses, as low as ∼0.01 M _J, in order to spin-up the envelopes of AGB stars efficiently. To support this claim, I follow the angular momentum evolution of single stars with main-sequence mass in the range of 1.3–2.4 M_⊙ , as they evolve to the post-AGB phase. I find that single stars rotate much too slowly to possess any significant non-spherical mass-loss as they reach the upper AGB. It seems, therefore, that planets, in some cases even Earth-like planets, are sufficient to spin-up the envelope of these AGB stars for them to form elliptical PNe. The prediction that on average several such planets orbit each star, as in the Solar system, still holds.     

High-precision stellar abundances of the elements: methods and applications

Efficient spectrographs at large telescopes have made it possible to obtain high-resolution spectra of stars with high signal-to-noise ratio and advances in model atmosphere analyses have enabled estimates of high-precision differential abundances of the elements from these spectra, i.e. with errors in the range 0.01–0.03 dex for F, G, and K stars. Methods to determine such high-precision abundances together with precise values of effective temperatures and surface gravities from equivalent widths of spectral lines or by spectrum synthesis techniques are outlined, and effects on abundance determinations from using a 3D non-LTE analysis instead of a classical 1D LTE analysis are considered. The determination of high-precision stellar abundances of the elements has led to the discovery of unexpected phenomena and relations with important bearings on the astrophysics of galaxies, stars, and planets, i.e. (i) Existence of discrete stellar populations within each of the main Galactic components (disk, halo, and bulge) providing new constraints on models for the formation of the Milky Way. (ii) Differences in the relation between abundances and elemental condensation temperature for the Sun and solar twins suggesting dust-cleansing effects in proto-planetary disks and/or engulfment of planets by stars; (iii) Differences in chemical composition between binary star components and between members of open or globular clusters showing that star- and cluster-formation processes are more complicated than previously thought; (iv) Tight relations between some abundance ratios and age for solar-like stars providing new constraints on nucleosynthesis and Galactic chemical evolution models as well as the composition of terrestrial exoplanets. We conclude that if stellar abundances with precisions of 0.01–0.03 dex can be achieved in studies of more distant stars and stars on the giant and supergiant branches, many more interesting future applications, of great relevance to stellar and galaxy evolution, are probable. Hence, in planning abundance surveys, it is important to carefully balance the need for large samples of stars against the spectral resolution and signal-to-noise ratio needed to obtain high-precision abundances. Furthermore, it is an advantage to work differentially on stars with similar atmospheric parameters, because then a simple 1D LTE analysis of stellar spectra may be sufficient. However, when determining high-precision absolute abundances or differential abundance between stars having more widely different parameters, e.g. metal-poor stars compared to the Sun or giants to dwarfs, then 3D non-LTE effects must be taken into account.     

Comets and the formation of planets

A morphological study of the physical and dynamical processes of planet formation is presented, with emphasis on the intermediary role of comet nuclei. Although guided by a particular model of the evolution of the pre-planetary solar nebula, implying the freezing-out of hydrogen in the region of the giant planets, the derivations and conclusions are of wider import, applicable to other cosmogonic models as well as to certain phases of star formation. The items evaluated physically, dynamically, or statistically comprise: (1) the total number mass of comets in Oort's cloud; (2) a re-evaluation of the diameters and masses of cometary nuclei; (3) the processes of nucleation from gravitational and Boltzmann instabilities of gaseous media to agglomerations of particulate matter as conditioned by inbuilt angular momentum; (4) the statistical-dynamical conditions and time scales of orbital interaction of comets with the planets and the consequences of disintegration.A consistent model proposes the formation of comets and planets in pre-planetary rings of the residual solar nebula, with subsequent ejection, chiefly by Jupiter, of the comets to Oort's sphere. Screening by absorbing matter is not only probable, but necessary to protect the comets from dis-integration during the process of ejection.Paper dedicated to Prof. H. C. Urey on the occasion of his 80th birthday on 29 April, 1973.This work has been currently supported by grants from the National Aeronautics and Space Administration, Goddard Space Flight Center, Greenbelt, Maryland.     

On the relationship of long-period comets to known and unknown planets

Statistical aspects of the question of the dynamical relationship of long-period comets to giant planets (Saturn, Uranus, Neptune), dwarf planets (Pluto, 2003 EL₆₁, 2003 UB₃₁₃), and five hypothetical planets are investigated. Data for 859 and 888 comets (2005, 2006) with periods P > 200 years are used in the statistics. No comets of the Kreutz, Marsden, Kracht, and Meyer groups are considered. The minimum interorbital distances of comets from the listed planetary bodies are mainly analyzed. Effects testifying to the kinematical relationship of some of the comets to planets have been established by testing the cometary data relative to 67 planes. The distribution of the perihelia of long-period comets is also discussed.