Energization of positive ions in the cometary foreshock region

The energization of positive ions in front of a cometary bow shock is investigated. Ions produced by ionization of the cometary neutrals interact with the solar wind protons to produce, among other waves, large amplitude oscillations of the ambient magnetic field. Such oscillations are convected towards the comet at the unperturbed solar wind speed far from the shock and at a lower speed closer to the shock (due to the solar wind mass loading) ; hence, they can energize the suprathermal ions by Fermi acceleration. The spatial extension of the acceleration region is of the order of 10⁶ km and the resulting ion energy spectrum is harder than in the Earth's bow shock case. The energization of cometary ions produces an additional deceleration of the solar wind, such that the cometary bow shock of Halley-type comet may be regarded as a “cosmic ray shock”.     

Cometary impact and the magnetization of the Moon

Collisions of comets with planetary bodies are capable of impressing patterns of magnetization onto them that match those observed for the Moon and possibly for Mercury. The ambient solar wind magnetic field is briefly but strongly enhanced as the large partially ionized cometary atmosphere is compressed against the planetary surface. Just at the time of peak field enhancement, the solid part of the comet collides with the surface and the compressed fields are permanently imprinted by shock magnetization.     

Spectral survey of the Comet Halley atmosphere during an occultation of a star

On November 15, 1985 a star was occulted by Comet Halley. Seven consequent spectra of the star has been obtained with 5-minute exposition with the TV scanner at 6-meter telescope. Simultaneously seven spectra of the cometary atmosphere were obtained at the projected distance 40” from the star. A method is developed to calculate the optical thickness of the cometary atmosphere as a function of wave length and the nucleocentric distance. This function is given in the spectral range 3925 Å−4850 Å.     

Currents in the cometary atmosphere

If the structure of the magnetic field and electric current in the cometary type I tail can be represented by an electric current circuit, disruption of the cross-tail current system may lead to a current discharging through the cometary ionosphere, and the dissipation of the magnetic energy stored in the tail. From the point of view of energy budget, a tail-aligned magnetic field on the order of 10γ will be sufficient to produce a strong ionization effect of the cometary atmosphere.     

MHD Simulation of Comets: The Plasma Environment of Comet Hale-Bopp

MHD simulation results of the interaction of the expanding atmosphere of comet Hale-Bopp with the magnetized solar wind are presented. At the upstream boundary a supersonic and superalfvénic solar wind enters into the simulation box 25 million km upstream of the nucleus. The solar wind is continuously mass loaded with cometary ions originating from the nucleus. The effects of photoionization, recombination and ion-neutral frictional drag are taken into account in the model. The governing equations are solved on an adaptively refined unstructured Cartesian grid using our MUSCL-type upwind numerical technique, MAUS-MHD (Multiscale Adaptive Upwind Scheme for MHD). The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the diamagnetic cavity of the comet. Detailed simulation results for the plasma environment of comet Hale-Bopp for slow and fast solar wind conditions are presented. We also calculate synthetic H₂O⁺, CO⁺ and soft x-ray images for observing conditions on April 11, 1997. This revised version was published online in July 2006 with corrections to the Cover Date.     

The multiple modes of interaction of the solar wind with a comet as it approaches the Sun

A quasi-steady 1-D hydrodynamic model, with mass addition, has been used to study the various modes of interaction of the solar wind with a medium-bright, H₂O-dominated comet (such as P/Halley) approaching the Sun.At large heliocentric distances (d 5 AU) the solar wind penetrates unimpeded on to the surface. As the comet moves further in, mass loading of the solar wind by heavy ions from the fledgling cometary atmosphere causes it to slow down, thereby causing a significant enhancement of the interplanetary field. Still further in at d 3.14 AU, as the mass loading reaches a critical value, a collision-less standing shock is formed in the solar wind upstream of the nucleus. As d decreases further, the distance of this shock from the nucleus increases. The cometary atmosphere becomes dense enough to stand off the solar wind ahead of the nuclear surface and form a well defined tangential discontinuity surface (or ionopause) only when d reaches the value 2.65 AU. When d 2.65 AU an inner shock could, in principle, also form within the cometary ionosphere, although its existence would depend on the detailed thermodynamics of the cometary ionosphere. Resolution of this question is beyond the scope of the present analysis.The conclusions of the present study would be qualitatively valid for other comets having sizes, surface optical properties and chemical compositions, different from those adopted here. The helio-centric distances at which the various transitions take place from the one mode of solar wind interaction to another, would, of course, be different, with all these distances being smaller for less active comets.     

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.     

Influence of the magnetic field on the density distribution of solar wind protons and cometary ions in the shock layer ahead of cometary ionospheres

The “mass loading” of the solar wind by cometary ions produced by the photoionization of neutral molecules outflowing from the cometary nucleus plays a major role in the interaction of the solar wind with cometary atmospheres. In particular, this process leads to a decrease in the solar wind velocity with a transition from supersonic velocities to subsonic ones through the bow shock. The so-called single-fluid approximation, in which the interacting plasma flows are considered as a single fluid, is commonly used in modeling such an interaction. However, it is occasionally necessary to know the distribution of parameters for the components of the interacting plasma flows. For example, when the flow of the cometary dust component in the interplanetary magnetic field is considered, the dust particle charge, which depends significantly on the composition of the surrounding plasma, needs to be known. In this paper, within the framework of a three-dimensional magnetohydrodynamic model of the solar wind flow around cometary ionospheres, we have managed to separately obtain the density distributions of solar wind protons and cometary ions between the bow shock and the cometary ionopause (in the shock layer). The influence of the interplanetary magnetic field on the position of the point of intersection between the densities with the formation of a region near the ionopause where the proton density is essentially negligible compared to the density of cometary ions is investigated. Such a region was experimentally detected by the Vega-2 spacecraft when investigating Comet Halley in March 1986. The results of the model considered below are compared with some experimental data obtained by the Giotto spacecraft under the conditions of flow around Comets Halley and Grigg–Skjellerup in 1986 and 1992, respectively. Unfortunately, our results of calculations on Comet Churyumov–Gerasimenko are only predictive in character, because the trajectory of the Rosetta spacecraft, which manoeuvred near its surface for several months, is complex.     

The effect of heavy ions on the formation and structure of cometary bow shocks

We investigate the interaction of heavy cometary ions with the solar wind and the formation of a bow shock in front of a comet by means of a hybrid (particle ion, fluid electron) simulation code that solves self-consistently for the electromagnetic fields and the motion of the charged particles. This kinetic treatment of the solar wind protons and the heavy cometary ions allows us to examine two important issues. One is the effect of the velocity distribution function of the heavy ions on the shock formation and structure, and the other is the degree of coupling between the two ion species. The result of this study indicate that at high Mach numbers the shock structure is highly dependent upon the velocity distribution of the heavy ions. For example, when the newly created ions comprise a ring distribution in the solar wind frame, most of them turn around downstream of the shock surface and reenter the upstream region to form a large foot that extends about a heavy ions gyroradius upstream of the shock. On the other hand, heavy ions which have been picked up by the solar wind and possesses a Maxwellian distribution can mostly penetrate the shock without returning upstream and affecting the shock structure as much. In either case, however, at high Mach numbers the shock strength is the same. At low Mach numbers, where the shock is weak, the velocity distribution of heavy ions has a smaller effect on the formation of shock and its structure. In this regime, the degree of coupling between the cometary ions and the solar wind protons and the corresponding critical Mach number (at which a shock should begin to form) are determined from a set of Rankine-Hugoniot relations. The results of the simulations suggest that some coupling does occur (evidently, through the electromagnetic fields, since there are no particle collisions in the calculations), but less than that expected from magnetohydrodynamics. For low Mach numbers, it is also shown that shocks have a transitory nature, where they are continuously formed by the protons and subsequently destroyed by the heavy ions.     

Direct statistical simulation of the near-surface layers of a cometary atmosphere. II: A nonspherical nucleus

The study presents the results of numerical simulations of mass-transfer processes in the near-surface layer of the cometary nucleus and in the inner part of the cometary atmosphere, which is formed under the action of solar radiation. The gas-kinetic model of the inner part of the cometary atmosphere surrounding a spherical nucleus (Skorov et al., 2004) is extended to the case of a nonspherical nucleus with axial symmetry. After high-resolution images of comets 19P/Borrelly and Wild 2 have been obtained by Deep Space 1 and Stardust spacecraft, such an extension seems to be vital and important. The nucleus and the inner part of the coma are closely related to each other because of the permanent exchange of energy and mass; therefore, they are modeled consistently. As in the first part of our study, the boundary conditions at the inner boundary of the simulation domain, which are necessary for gas-kinetic simulations, were determined from the self-consistent model of heat and mass transfer in a porous cometary nucleus that was developed earlier by the authors. The model took into account the volumetric character of the radiation absorption in a porous sublimating medium, the kinetic regime of the transport of sublimation products in the pores, and the backward gas fluxes from the coma due to intermolecular collisions. We considered different models of the nucleus structure that determined the effective gas production. Using the direct simulation Monte Carlo method, we computed the two-dimensional gas flow from a heterogeneous nonspherical cometary nucleus. The simulations were performed using the SMILE software. The parallel computer implementation of the software made it possible to calculate the spatial structure of the gas flow for the entire circumnucleus zone.     

Alfvén wave plasma turbulence during solar wind-comet interaction

The large differences in drift velocities between the solar wind protons and the picked-up ions of cometary origin cause the Alfvén waves (among others) to become unstable and generate turbulence. A self-consistent treatment of such instabilities has to take into account that these cometary ions affect the solar wind plasma in a decisive way. With the help of a previously developed formalism one finds the correct Alfvén instability criterion, which is here nondispersive, in contrast to recent calculations where the cometary ions are treated as a low-density, high-speed, and non-neutral beam through an otherwise undisturbed solar wind. The true bulk speed of the combined solar wind plus cometary ion plasma clearly shows the mass-loading and deceleration of the solar wind near the cometary nucleus, indicating a bow shock. The instability criterion is also used to determine the region upstream where the Alfvén waves can be unstable, based upon recent observations near comet Halley.     

An experimental study of the isotopic enrichment in Ar, Kr, and Xe when trapped in water ice

The isotopic enrichment of argon, krypton, and xenon, when trapped in water ice, was studied experimentally. The isotopes were found to be enriched according to their (m1/m2)1/2 ratio. These enrichment factors could be useful for comparison among the uncertain cosmic or solar isotopic ratios, the hopeful in situ cometary ratio, and those in Earth's atmosphere, in the context of cometary delivery of volatiles to Earth.     

Determination of the height of the “meteoric explosion”

When cosmic bodies of asteroidal and cometary origin, with a size from 20 to approximately 100 m, enter dense atmospheric layers, they are destroyed with a large probability under the action of aerodynamic forces and decelerated with the transfer of their energy to the air at heights from 20–30 to several kilometers. The forming shock wave reaches the Earth’s surface and can cause considerable damage at great distances from the entry path similar to the action of a high-altitude explosion. We have performed a numerical simulation of the disruption (with allowance for evaporation of fragments) and deceleration of meteoroids having the aforesaid dimensions and entering the Earth’s atmosphere at different angles and determined the height of the equivalent explosion point generating the same shock wave as the fall of a cosmic body with the given parameters. It turns out that this height does not depend on the velocity of the body and is approximately equal to the height at which this velocity is reduced by half. The obtained results were successfully approximated by a simple analytical formula allowing one to easily determine the height of an equivalent explosion depending on the dimensions of the body, its density, and angle of entry into the atmosphere. A comparison of the obtained results with well-known approximate analytical (pancake) models is presented and an application of the obtained formula to specific events, in particular, to the fall of the Chelyabinsk meteorite on February 15, 2013, and Tunguska event of 1908, is discussed.     

The volatile fraction of the cometary nucleus

In order to prepare a flyby mission to Comet Encke, six different sources of information on the possible chemical composition of the cometary nucleus are compared. These are: the neutral and charged radicals and molecules observed in cometary spectra; the chemical composition of type I carbonaceous chondrites; the meteor spectra; the metallic ions collected in the upper atmosphere and correlated with the meteor shower associated with Comet Encke; and finally the volatile molecules observed in a volatile-rich sample of lunar soil, that were interpreted as a possible cometary impact. Possible molecular abundances for the volatile fraction of Comet Encke are tentatively proposed.     

The dual origin of the terrestrial atmosphere

The origin of the terrestrial atmosphere is one of the most puzzling enigmas in the planetary sciences. It is suggested here that two sources contributed to its formation, fractionated nebular gases and accreted cometary volatiles. During terrestrial growth, a transient gas envelope was fractionated from nebular composition. This transient atmosphere was mixed with cometary material. The fractionation stage resulted in a high Xe/Kr ratio, with xenon being more isotopically fractionated than krypton. Comets delivered volatiles having low Xe/Kr ratios and solar isotopic compositions. The resulting atmosphere had a near-solar Xe/Kr ratio, almost unfractionated krypton delivered by comets, and fractionated xenon inherited from the fractionation episode. The dual origin therefore provides an elegant solution to the long-standing “missing xenon” paradox. It is demonstrated that such a model could explain the isotopic and elemental abundances of Ne, Ar, Kr, and Xe in the terrestrial atmosphere.     

The neutral atmospheres of comets

In this paper we have endeavored to critically evaluate our present understanding of cometary atmospheres. Following a brief introduction of the significance of the study of cometary atmospheres (Section 1), the relevant photometric and spectroscopic observations are summarized in Section 2.The interaction with the solar radiation, with regard to both the excitation of the observed species as well as the dissociation of stable molecules evaporating from the nucleus, is considered in Sections 3 and 4. The gas phase chemistry likely to take place in the dense inner coma is next considered in Section 5.The exospheric and hydrodynamic models of the expanding cometary atmosphere are considered in detail in Section 6, and both their limitations as well as possible improvements are discussed.The observed chemical composition of the neutral atmosphere and the inferred chemical composition of the volatile component of the nucleus, together with possible variations between different classes of comets is next considered in Section 7, and their possible cosmogonic significance is discussed.In conclusion, some of the important directions in which future research should progress, in order to provide more complete and secure knowledge of cometary atmospheres, are stressed (Section 8).Astrophysics and Space Science Review Paper.     

Interrelationships among meteoric metals,meteors, interplanetary dust,micrometeorites, and meteorites

Abstract— Meteor science, aeronomy, and meteoritics are different disciplines with natural interfaces. This paper is an effort to integrate the chemistry and mineralogy of collected interplanetary dust particles (IDPs), micrometeorites, and meteorites with meteoric data and with atmospheric metal abundances. Evaporation, ablation, and melting of decelerating materials in the Earth's atmosphere are the sources of the observed metal abundances in the upper atmosphere. Many variables ultimately produce the materials and phenomena we can analyze, such as different accretion and parent‐body histories of incoming extraterrestrial materials, different interactions of meteors with the Earth's middle atmosphere, meteor data reduction, and complex chemical interactions of the metals and ions with the ambient atmosphere. The IDP‐like and unequilibrated ordinary chondrite matrix materials are reasonable sources for observed meteoric and atmospheric metals. The hypothesis of hierarchical dust accretion predicts that low, correlated refractory element abundances in cometary meteors may be real. It implies that the CI or cosmic standard is not useful to appreciate the chemistry of incoming petrologically heterogeneous cometary matter. The quasi steady‐state metal abundances in the lower thermosphere and upper mesosphere are derived predominantly from materials with cometary orbital characteristics and velocities such as comets proper and near‐Earth asteroids. The exact influence of atmospheric chemistry on these abundances still needs further evaluation. Metal abundances in the lower mesosphere and upper stratosphere region are mostly from materials from the asteroidal belt and the Kuiper belt.     

Amino acid survival in large cometary impacts

Abstract— A significant fraction of the Earth's prebiotic volatile inventory may have been delivered by asteroidal and cometary impacts during the period of heavy bombardment. The realization that comets are particularly rich in organic material seemed to strengthen this suggestion. Previous modeling studies, however, indicated that most organics would be entirely destroyed in large comet and asteroid impacts. The availability of new kinetic parameters for the thermal degradation of amino acids in the solid phase made it possible to readdress this question. We present the results of new high-resolution hydrocode simulations of asteroid and comet impact coupled with recent experimental data for amino acid pyrolysis in the solid phase. Differences due to impact velocity as well as projectile material have been investigated. Effects of angle of impacts were also addressed. The results suggest that some amino acids would survive the shock heating of large (kilometer-radius) cometary impacts. At the time of the origins of life on Earth, the steady-state oceanic concentration of certain amino acids (like aspartic and glutamic acid) delivered by comets could have equaled or substantially exceeded concentrations due to Miller-Urey synthesis in a CO₂-rich atmosphere. Furthermore, in the unlikely case of a grazing impact (impact angle ～5° from the horizontal), an amount of some amino acids comparable to that due to the background steady-state production or delivery would be delivered to the early Earth.