Suprathermal oxygen atoms in the Martian upper atmosphere: Contribution of the proton and hydrogen atom precipitation

This is a study of the kinetics and transport of hot oxygen atoms in the transition region (from the thermosphere to the exosphere) of the Martian upper atmosphere. It is assumed that the source of the hot oxygen atoms is the transfer of momentum and energy in elastic collisions between thermal atmospheric oxygen atoms and the high-energy protons and hydrogen atoms precipitating onto the Martian upper atmosphere from the solar-wind plasma. The distribution functions of suprathermal oxygen atoms by the kinetic energy are calculated. It is shown that the exosphere is populated by a large number of suprathermal oxygen atoms with kinetic energies up to the escape energy 2 eV; i.e., a hot oxygen corona is formed around Mars. The transfer of energy from the precipitating solar-wind plasma protons and hydrogen atoms to the thermal oxygen atoms leads to the formation of an additional nonthermal escape flux of atomic oxygen from the Martian atmosphere. The precipitation-induced escape flux of hot oxygen atoms may become dominant under the conditions of extreme solar events, such as solar flares and coronal mass ejections, as shown by recent observations onboard NASA’s MAVEN spacecraft (Jakosky et al., 2015).     

Stochastic models of hot planetary and satellite coronas: A hot oxygen corona of Mars

The processes of the kinetics and transport of hot oxygen atoms in the transition (between thermosphere and exosphere) region of the upper atmosphere of Mars are studied. The reaction of dissociative recombination of the main ionospheric ion O ₂ ⁺ with thermal electrons in the ionosphere of Mars is considered as a source of hot oxygen atoms. The distribution of suprathermal oxygen atoms by kinetic energy is calculated. It is shown that the exosphere is populated by a considerable number of suprathermal oxygen atoms with kinetic energies just below the escape energy of 2 eV; that is, a hot oxygen corona of Mars is formed.     

Stochastic models of hot planetary and satellite coronas: A photochemical source of hot oxygen in the upper atmosphere of Mars

The processes of the kinetics and transport of hot oxygen atoms in the upper atmosphere of Mars are studied. A reaction of dissociative recombination of the main ionospheric ion O ₂ ⁺ with thermal electrons is considered as a photochemical source of suprathermal oxygen atoms. Oxygen atoms are formed in the dissociative recombination reaction with an excess of kinetic energy of about 0.4–4 eV and lose that energy in elastic and inelastic collisions with the ambient thermal atmospheric gas. The altitude distributions of the concentrations of neutral and ionized components, as well as their temperatures, were taken from Krasnopolsky (2002). Unlike the models published earlier, detailed calculations of the formation, collisional kinetics, and transport of suprathermal oxygen atoms in the thermosphere-exosphere transition region of the upper atmosphere of Mars have been made for the first time. For this, we used a stochastic model of the formation of a hot planetary corona (Shematovich, 2004). It has been shown that the considered photochemical source of suprathermal oxygen leads to the formation of the hot corona and to higher nonthermal losses of oxygen from the upper atmosphere of Mars due to escape fluxes. The detailed energy spectra of the fluxes of suprathermal atomic oxygen were calculated for the thermosphere-exosphere transition region of the Martian atmosphere.Translated from Astronomicheskii Vestnik, Vol. 39, No. 1, 2005, pp. 26–37.Original Russian Text Copyright © 2005 by Krestyanikova, Shematovich.     

Suprathermal hydrogen produced by the dissociation of molecular hydrogen in the extended atmosphere of exoplanet HD 209458b

The ionization and dissociation of molecular hydrogen by the ultraviolet (UV) radiation of the parent star lead to the formation of hydrogen atoms with an excess of kinetic energy and, thus, are an important source of suprathermal hydrogen atoms in the upper atmosphere of exoplanet HD 209458b. Contemporary aeronomical models did not investigate these processes because they assumed the fast local thermalization of the hot atoms of hydrogen by elastic collisions. However, the kinetics and transfer of these atoms were not calculated in detail, because they require the solving of the Boltzmann equation for a nonthermal atom population. This work estimates the effect of the UV radiation of the parent star and the accompanying photocleacton flux on the production of the suprathermal fraction of atomic hydrogen in the H₂ → H transition region. We also consider the formation of the escaping flux of Hatoms created by this effect in the upper atmosphere of HD 209458b. We calculate the production rate and energy spectrum of the hydrogen atoms with excess kinetic energy during the dissociation of H₂. Using the numerical stochastic model created by Shematovich (2004) for a hot planetary corona, we investigate the molecular-scale kinetics and transfer of suprathermal hydrogen atoms in the upper atmosphere and the emergent flux of atoms evaporating from the atmosphere. The latter is estimated as 3.4 × 10¹² cm⁻² s⁻¹ for a moderate stellar activity level of UV radiation, which leads to a planetary atmosphere evaporation rate of 3.4 × 10⁹ g s⁻¹ due to the process of the dissociation of H₂. This estimate is close to the observational value of ∼10¹⁰ g s⁻¹ for the rate of atmospheric loss of HD 209458b.     

Upper limits for the Martian exospheric number density during the Planet B/Nozomi mission

Dissociative recombination (DR) of ionospheric O₂⁺ ions is an important source of suprathermal atomic oxygen in the exosphere as previous studies about the Martian upper atmosphere have shown. Because of the weaker gravitational attraction a hot oxygen corona on Mars should be denser than that observed on Venus. Since the most important mechanism for the production of the hot oxygen atoms in the Martian exosphere is DR, we investigated the variability of this production mechanism depending of solar activity. The Japanese Nozomi spacecraft will have the possibility to detect with the neutral mass spectrometer (NMS) for the first time in-situ the theoretically predicted hot oxygen corona on Mars, if the corona number density above the cold background atmosphere is of the order of 10,000 cm⁻³. Due to a problem in the propulsion system Nozomi failed its planned arrival rendevouzs with Mars in October 1999 and will, therefore, arrive at the red planet not before January 2004. Solar activity will reach its maximum in 2001, so the related production rate of hot oxygen atoms will be in the medium range during the new arrival date of Nozomi. We used the ionospheric profiles from the Viking mission for low solar activity conditions (F_10.7≈70) and the Mariner 9 mission with a solar activity of about 120 for medium solar wind activity. The latter is comparable to the level we expect for the Mars arrival of Nozomi. The resulting influence of the hot oxygen corona number density distribution was calculated with a Monte Carlo technique. This technique is used to compute a hot particle density distribution function. We studied the atomic diffusion process in the Martian atmosphere by simulating the collision probability, particle direction and energy loss after collisions by generating random numbers. Compared to previous studies we have improved the Monte Carlo model by using more and smaller altitude steps and more detailed treatment of particles with a temporary downward motion. This has resulted in an increased amount of collisions and a shift to lower energies in the energy spectrum. Our results show that the hot oxygen component should begin to dominate above the cold background atmosphere at an altitude of about 500 km above the Martian surface. The NMS instrument on board of Nozomi should detect the hot oxygen component after its arrival at Mars in January 2004, at an altitude of about 600 km above the Martian surface. Since the solar activity will decrease during the mission the measurements during the first orbits will be the most significant ones. The first in-situ measurements of the hot oxygen number density would be very important for adjusting atmospheric escape models by separating ballistic, satellite and escape trajectories of the hot oxygen atoms, which are significant for studies of the evolution and solar wind interaction of the Martian atmosphere.     

Stochastic models of hot planetary and satellite coronas: Atomic oxygen in Europa’s corona

This paper analyzes the formation, kinetics, and transport of hot oxygen atoms in the atmosphere of the Jovian satellite Europa. Atmospheric sources of suprathermal oxygen atoms are assumed to be represented by the processes of dissociation of molecular oxygen, which is the main component of the atmosphere, by solar UV radiation and electron fluxes from the inner magnetosphere of Jupiter, as well as by the reaction of dissociative recombination of the main ionospheric ion O ₂ ⁺ which thermal electrons. It is shown that dissociation in Europa’s near-surface atmosphere is balanced by the processes of the loss of atomic oxygen due to the effective escape of suprathermal oxygen atoms into the inner magnetosphere of Jupiter along the orbit of Europa and due to ionization by magnetospheric electrons and catalytic recombination of oxygen atoms on the icy surface of the satellite. It thus follows that atomic oxygen is only a small admixture to the main atmospheric component—molecular oxygen—in the near-surface part of the atmosphere. However, the outer exospheric layers of Europa’s atmosphere are populated mostly by suprathermal oxygen atoms. The near-surface molecular envelope of Europa is therefore surrounded by a tenuous extended corona of hot atomic oxygen.     

Helium in the Martian atmosphere: Thermal loss considerations

Helium concentrations in the Martian atmosphere are estimated assuming that the helium production on Mars, comparable to its production on Earth, via the radioactive decay of uranium and thorium, is in steady state equilibrium with its thermal escape. Although non-thermal losses would tend to reduce the estimated concentrations, these concentrations are not necessarily an upper limit since higher production rates and/or a possibly lower effective exospheric temperature over the solar activity cycle could increase them to even higher values. The computed helium concentration at the Martian exobase (200 km) is 8 × 10⁶ atoms cm^?3. Through the lower exosphere, the computed helium concentrations are 30–200 times greater than the Mariner-measured atomic hydrogen concentrations. It follows that helium may be the predominant constituent in the Martian lower exosphere and may well control the orbital lifetime of Mars-orbiting spacecraft. The estimated helium mixing ratio is greater at the Martian turbopause than at the terrestrial turbopause, and the helium column density in the lower Martian atmosphere may be comparable to that on Earth.     

Influence of the hot oxygen corona on the satellite drag in the Earth’s upper atmosphere

Calculation results on the possible influence of the hot oxygen fraction on the satellite drag in the Earth’s upper atmosphere on the basis of the previously developed theoretical model of the hot oxygen geocorona are presented. Calculations have shown that for satellites with orbits above 500 km, the contribution from the corona is extremely important. Even for the energy flux Q ₀ = 1 erg cm⁻² s⁻¹, the contribution of the hot oxygen can reach tens of percent; and considering that real energy fluxes are usually higher, one can suggest that for extreme solar events, the contribution of hot oxygen to the atmospheric drag of the satellite will be dominant. For lower altitudes, the contribution of hot oxygen is, to a considerable degree, defined by the solar activity level. The calculations imply that for the daytime polar atmosphere, the change of the solar activity level from F _10.7 ∼ 200 to F _10.7 ∼ 70 leads to an increase in the ratio of the hot oxygen partial pressure to the thermal oxygen partial pressure by a factor of almost 30, from 0.85 to 25%. The transition from daytime conditions to nighttime conditions almost does not change the contribution from suprathermal particles. The decrease of the characteristic energy of precipitating particles, i.e., for the case of charged particles with a softer energy spectrum, leads to a noticeable increase of the contribution of the suprathermal fraction, by a factor of 1.5–2. It has been ascertained that electrons make the main contribution to the formation of the suprathermal fraction; and with the increase of the energy of precipitating electrons, the contribution of hot oxygen to the satellite drag also increases proportionally. Thus, for a typical burst, the contribution of the suprathermal fraction is 30% even at relatively high solar activity F _10.7 = 135.     

Stochastic Models of Hot Planetary and Satellite Coronas: Suprathermal Nitrogen in Titan's Upper Atmosphere

The processes of dissociation and dissociative ionization of molecular nitrogen by solar UV radiation and by the accompanying flux of photoelectrons, as well as sputtering of the atmosphere by fluxes of magnetospheric ions and pick-up ions, are the main sources of translationally excited (hot, or suprathermal) nitrogen atoms and molecules in Titan's upper atmosphere. Since Titan does not possess an intrinsic magnetic field, ions from Saturn's magnetosphere can penetrate into the outer layers of Titan's atmosphere and sputter atoms and molecules from the atmosphere in momentum-transfer and charge exchange collisions. Atmospheric sputtering by corotating nitrogen ions and carbon-containing pick-up ions, as well as photodissociation-related losses, was considered previously by Lammer and Bauer (1993) and Shematovich et al. (2001, 2003). In this paper we investigate the processes of the formation and evolution of the fraction of suprathermal nitrogen atoms and molecules in the transition region of Titan's upper atmosphere using the previously developed Monte Carlo model for hot satellite and planetary coronas (Shematovich, 1999, 2004). It is established that the suprathermal nitrogen fraction in the transition region of Titan's upper atmosphere includes nitrogen atoms and molecules but the suprathermal nitrogen concentration is relatively small owing to high rates of escape from the atmosphere and to the efficient thermalization of suprathermal nitrogen at the altitudes of the relatively dense lower thermosphere. However, the scale height for suprathermal nitrogen in the transition region is much higher than that for the ambient atmospheric gas. Therefore, suprathermal nitrogen becomes one of the dominant components in the outer exosphere.     

Energy transfer in oxygen-hydrogen collisions

A formula is obtained for the rate coefficient for the production of a particle with a specific energy in an elastic collision with the atoms of a thermal gas. Accurate calculations are carried out of the energy transfer in elastic collisions of oxygen and hydrogen atoms. The fractions of hydrogen atoms with sufficient energy to escape from Venus produced in collisions of oxygen atoms with energies of 2.6 eV in a gas of atomic hydrogen at 100, 200 and 300 K are respectively 5.1, 6.9 and 8.5%.The mutual diffusion coefficient D of oxygen and hydrogen atoms is obtained. It can be represented by Dn = 7.25 × 10¹⁷T^0.71cm^?1s^?1 where n is the total atom density and T is the temperature.     

Chemistry and evolution of Titan's atmosphere

The chemistry and evolution of Titan's atmosphere is reviewed in the light of the scientific findings from the Voyager mission. It is argued that the present N₂ atmosphere may be Titan's initial atmosphere rather than photochemically derived from an original NH₃ atmosphere. The escape rate of hydrogen from Titan is controlled by photochemical production from hydrocarbons. CH₄ is irreversibly converted to less hydrogen rich hydrocarbons, which over geologic time accumulate on the surface to a layer thickness of ～0.5 km. Magnetospheric electrons interacting with Titan's exosphere may dissociate enough N₂ into hot, escaping N atoms to remove ～0.2 of Titan's present atmosphere over geologic time. The energy dissipation of magnetospheric electrons exceeds solar e.u.v. energy deposition in Titan's atmosphere by an order of magnitude and is the principal driver of nitrogen photochemistry. The environmental conditions in Titan's upper atmosphere are favorable to building up complex molecules, particularly in the north polar cap region.     

Titan's corona: The contribution of exothermic chemistry

The contribution of exothermic ion and neutral chemistry to Titan's corona is studied. The production rates for fast neutrals N₂, CH₄, H, H₂, ³CH₂, CH₃, C₂H₄, C₂H₅, C₂H₆, N(⁴S), NH, and HCN are determined using a coupled ion and neutral model of Titan's upper atmosphere. After production, the formation of the suprathermal particles is modeled using a two-stream simulation, as they travel simultaneously through a thermal mixture of N₂, CH₄, and H₂. The resulting suprathermal fluxes, hot density profiles, and energy distributions are compared to the N₂ and CH₄ INMS exospheric data presented in [De La Haye, V., Waite Jr., J.H., Johnson, R.E., Yelle, R.V., Cravens, T.E., Luhmann, J.G., Kasprzak, W.T., Gell, D.A., Magee, B., Leblanc, F., Michael, M., Jurac, S., Robertson, I.P., 2007. J. Geophys. Res., doi:10.1029/2006JA012222, in press], and are found insufficient for producing the suprathermal populations measured. Global losses of nitrogen atoms and carbon atoms in all forms due to exothermic chemistry are estimated to be and .     

The calcium exosphere of Mercury

A tenuous calcium atmosphere at Mercury, principally seen in the polar regions, was first observed in July, 1998, using the High Resolution Echelle Spectrograph (HIRES) at the W.M. Keck I telescope (Bida et al., Nature 404, 159, 2000). We report four years of observations of the calcium exosphere of Mercury, confirming the initial findings of a very tenuous atmosphere. These observations show a persistent but spatially variable blue shift, indicating an excess velocity toward the observer of up to 3 km s⁻¹, with an average excess velocity of 2.2 km s⁻¹ above the south pole. In addition, the line profiles reveal a hot corona at the equivalent of 12,000-20,000 K in a thermalized atmosphere, indicating a large range of motion with respect to the observer. The calcium is not confined to the polar-regions: rare and low Ca abundance is seen in the equatorial regions. Strong emission was seen anti-sunward on 3 May 2002. Apparent weak emission on the sunward hemisphere may be due to scattered light from the surface, or may indicate a high latitude source. We show that the likely source of the calcium is either impact vaporization in the form of CaO and clusters, which are subsequently photo-dissociated, or ion-sputtering of atoms, molecules and ions. The column abundance is somewhat, but not strongly, correlated with solar activity. We predict a very hot (probably escaping) oxygen component to the hermean exosphere.     

Stochastic Models of Hot Planetary and Satellite Coronas

The hot planetary and satellite coronas are populated by the suprathermal particles produced in the transition region between the collision-dominated and free-molecule atmospheric layers under the external effects of electromagnetic and corpuscular solar radiation and magnetospheric plasma. We construct a numerical stochastic model to investigate both the local formation and kinetics of suprathermal particles and their transport to exospheric heights from underlying atmospheric layers. In contrast to other commonly used approaches, the suggested numerical model is suitable for studying the flows of atmospheric gas weakly and strongly perturbed by suprathermal particles, i.e., for studying the formation of hot planetary and satellite coronas proper. Highly efficient Monte-Carlo algorithms with weighted particles underlie the numerical implementation of the model. This numerical model is used to investigate the following: (i) the hot oxygen corona of Europa, a Jovian satellite, which is an example of a highly nonequilibrium near-surface atmosphere; and (ii) the nonthermal losses of nitrogen from Titan, a Saturnian satellite, when suprathermal atoms and molecules of nitrogen are only a small admixture to the surrounding thermal molecular nitrogen—the main atmospheric component of Titan.     

Lyman-alpha observations of venera-9 and 10 I. The non-thermal hydrogen population in the exosphere of venus

The dayside hydrogen exosphere was observed in October–November 1975 with a Lymanalpha photometer carried on board Venera-9 and 10. In addition to intensity measurements, the use of a hydrogen cell allowed for the first time linewidth measurements. Both intensity and linewidth measurements below 1500 km of altitude are well fitted by a single exospheric component (T_c = 500 ± 100 K, n_c = 1.5 × 10⁴ atom cm^?3at 250 km). Above 3000 km, the measured linewidth increased sharply, to decrease again above 4500 km. This feature is interpreted as the signature of an additional population of “hot” atoms circulating on satellite orbits, created just behind the bow-shock by charge-exchange collisions (with an efficiency of 4%) between the neutral atoms and the solar wind protons, which became turbulent after bow-shock crossing. The density ratio of “hot” to standard population is of the order of 10% around 3500 km of altitude.     

The early thermal evolution of Mars

Hf‐W isotopic systematics of Martian meteorites have provided evidence for the early accretion and rapid core formation of Mars. We present the results of numerical simulations performed to study the early thermal evolution and planetary scale differentiation of Mars. The simulations are confined to the initial 50 Myr (Ma) of the formation of solar system. The accretion energy produced during the growth of Mars and the decay energy due to the short‐lived radio‐nuclides ²⁶Al, ⁶⁰Fe, and the long‐lived nuclides, ⁴⁰K, ²³⁵U, ²³⁸U, and ²³²Th are incorporated as the heat sources for the thermal evolution of Mars. During the core‐mantle differentiation of Mars, the molten metallic blobs were numerically moved using Stoke's law toward the center with descent velocity that depends on the local acceleration due to gravity. Apart from the accretion and the radioactive heat energies, the gravitational energy produced during the differentiation of Mars and the associated heat transfer is also parametrically incorporated in the present work to make an assessment of its contribution to the early thermal evolution of Mars. We conclude that the accretion energy alone cannot produce widespread melting and differentiation of Mars even with an efficient consumption of the accretion energy. This makes ²⁶Al the prime source for the heating and planetary scale differentiation of Mars. We demonstrate a rapid accretion and core‐mantle differentiation of Mars within the initial ~1.5 Myr. This is consistent with the chronological records of Martian meteorites.     

The sodium and hydrogen gas clouds of Io

Models are developed to describe the spatial distribution of gases emitted by Io and are applied to recent observations which indicate extensive gas clouds of hydrogen and sodium in orbit around Jupiter. Hydrogen and sodium atoms are emitted from Io with velocities in the range 2 to 3 km sec^?1, with fluxes of about 10¹⁰ and 10⁸cm^?2sec^?1 for hydrogen and sodium respectively. Hydrogen atoms may be formed by photodecomposition of gases such as NH₃ or H₂S released from the satellite surface and may escape thermally from an exosphere whose temperature is about 500 K. Sodium may be ejected from the surface by energetic particles or by ultraviolet radiation and it appears that a non-thermal mechanism drawing energy from Jupiter's magnetic field is required in order to account for its release to space.     

Europa's Oxygen Exosphere and Its Magnetospheric Interaction

The newly detected oxygen atmosphere of Europa is modeled by invoking charged particle sputtering with H₂O and O₂molecules as the main ejecta. The magnetospheric corotating ions could provide the required source strength (∼3 × 10²⁶sec⁻¹) of O₂molecules if a fraction (∼20%) of the exospheric ions were recycled to Europa's surface where they produce additional sputtering product. Two exospheric components are expected to form: an extended corona with a size of a few satellite radii which is composed of sputtered molecules in ballistic motion, and a thermal population with a surface density of 10⁸–10⁹cm⁻³and a scale height of about 20 km. The electron impact ionization of this exosphere would lead to an Io-like interaction with the jovian magnetosphere with a field-aligned Birkeland current of about 5 × 10⁵A.