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.     

Suprathermal oxygen and hydrogen atoms in the upper Martian atmosphere

The processes of kinetics and transport of hot oxygen and hydrogen atoms in the transition (from the thermosphere to the exosphere) region of the upper Martian atmosphere are studied. The reaction of dissociative recombination of the principal ionospheric ion O ₂ ⁺ with thermal electrons in the ionosphere of Mars served as the source of hot oxygen atoms. The process of momentum and energy transfer in elastic collisions between hot oxygen atoms and atmospheric hydrogen atoms with thermal energies was regarded as the source of hot hydrogen atoms. The kinetic energy distribution functions are determined for suprathermal oxygen and hydrogen atoms. It is shown that the exosphere is populated with a significant number of suprathermal oxygen atoms with kinetic energies ranging up to the escape energy of 2 eV (i.e., the hot oxygen Martian corona is formed). The transfer of energy from hot oxygen atoms to thermal hydrogen atoms creates an additional nonthermal flux of atomic hydrogen escaping from 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.     

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: 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.     

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.     

Stochastic models of hot planetary and satellite coronas: Total water loss in the Martian atmosphere

Estimates of the total thermal and nonthermal losses of hydrogen and the total nonthermal loss of oxygen from the atmosphere of Mars are discussed, and their ratio is analyzed. It is shown that an H to O ratio of 2:1 has not been achieved in any of the current models of various authors. The closest ratio, H:O = 4:1, has been obtained by Krestyanikova and Shematovich (2006) in the model of formation of a hot oxygen corona.     

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.     

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.     

Rosetta-Alice observations of exospheric hydrogen and oxygen on Mars

The European Space Agency’s Rosetta spacecraft, en route to a 2014 encounter with comet 67P/Churyumov-Gerasimenko, made a gravity assist swing-by of Mars on 25 February 2007, closest approach being at 01:54 UT. The Alice instrument on board Rosetta, a lightweight far-ultraviolet imaging spectrograph optimized for in situ cometary spectroscopy in the 750-2000 Å spectral band, was used to study the daytime Mars upper atmosphere including emissions from exospheric hydrogen and oxygen. Offset pointing, obtained five hours before closest approach, enabled us to detect and map the H i Lyman-α and Lyman-β emissions from exospheric hydrogen out beyond 30,000 km from the planet’s center. These data are fit with a Chamberlain exospheric model from which we derive the hydrogen density at the 200 km exobase and the H escape flux. The results are comparable to those found from the Ultraviolet Spectrometer experiment on the Mariner 6 and 7 fly-bys of Mars in 1969. Atomic oxygen emission at 1304 Å is detected at altitudes of 400-1000 km above the limb during limb scans shortly after closest approach. However, the derived oxygen scale height is not consistent with recent models of oxygen escape based on the production of suprathermal oxygen atoms by the dissociative recombination of .     

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 .     

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.     

Comment on the papers “Production of hot nitrogen atoms in the martian thermosphere” by F. Bakalian and “Monte Carlo computations of the escape of atomic nitrogen from Mars” by F. Bakalian and R.E. Hartle

In recent articles published in Icarus, Bakalian [2006. Icarus 183, 69-78] discusses and computes the production rates of hot nitrogen atoms in the martian thermosphere due to N₂ photodissociation and N⁺₂ dissociative recombination, and Bakalian and Hartle [2006. Icarus 183, 55-68] use a Monte Carlo code to compute the escape rates of nitrogen atoms from Mars due to photodissociation of N₂, dissociative recombination of N⁺₂, and pickup ion escape due to photoionization of N atoms above the ionopause. Bakalian concludes that “photodissociation of N₂ is the dominant escape mechanism in the martian atmosphere.” We will show that this conclusion is not supportable. In addition, both papers contain scientific errors, misrepresentations, inaccurate referencing, lack of proper attribution, and they fail to place these investigations into the existing extensive body of work on this subject.     

Suprathermal electron fluxes on the nightside of Mars: ASPERA-3 observations

Recently aurora-type UV emissions were discovered on the nightside of Mars [Bertaux, J.-L., Leblanc, F., Witasse, O., et al., 2005. Discovery of an aurora on Mars. Nature 439, doi:10.1038/nature03603]. It was suggested that these emissions are produced by suprathermal electrons with energies of tens of eV, rather than by the electrons with spectra peaked above 100 eV [Leblanc, F., Witasse, O., Winningham J., et al., 2006. Origin of the martian aurora observed by spectroscopy for investigation of characteristics of the atmosphere of Mars (SPICAM) onboard Mars Express. J. Geophys. Res. 111, A09313, doi:10.1029/2006JA011763]. In this paper we present observations of fluxes of suprathermal electrons (E_e≈30-100 eV) on the Martian nightside by the ASPERA-3 experiment onboard the Mars Express spacecraft. Narrow spikes of suprathermal electrons are often observed in energy-time spectrograms of electron fluxes at altitudes between 250 and 600 km. These spikes are spatially organized and form narrow strips in regions with strong upward or downward crustal magnetic field. The values of electron fluxes in such events generally could explain the observed auroral UV emissions although a question of their origin (transport from the dayside or local precipitation) remains open.     

Exospheric density and escape fluxes of atomic isotopes on Venus and Mars

Energetic neutrals in dissociative recombinations near or above the exobase provide an important component of exospheric density and escape fluxes. Plasma thermal velocities provide the main contribution to the velocity spread and an exact integral for the escape flux applicable in marginal cases is found for a simple atmosphere and collisional cut-off. Atomic fragments from recombination of diatomic oxygen and nitrogen ions in the Venus and Mars atmospheres are examined and density integrals derived. The oxygen escape flux on Mars is half that previously estimated and there is very little isotope preference supplementing diffusive separation. However, escape of the heavier ¹⁵N isotope is low by a factor two. Reinterpretation of its 75% enrichment as detected by Viking leads to a range 0.4–1.4 mbar for the primeval nitrogen content on Mars.     

Changes in the Martian atmosphere induced by auroral electron precipitation

Typical auroral events in the Martian atmosphere, such as discrete and diffuse auroral emissions detected by UV spectrometers onboard ESA Mars Express and NASA MAVEN, are investigated. Auroral electron kinetic energy distribution functions and energy spectra of the upward and downward electron fluxes are obtained by electron transport calculations using the kinetic Monte Carlo model. These characteristics of auroral electron fluxes make it possible to calculate both the precipitation-induced changes in the atmosphere and the observed manifestations of auroral events on Mars. In particular, intensities of discrete and diffuse auroral emissions in the UV and visible wavelength ranges (Soret et al., 2016; Bisikalo et al., 2017; Gérard et al., 2017). For these conditions of auroral events, the analysis is carried out, and the contribution of the fluxes of precipitating electrons to the heating and ionization of the Martian atmosphere is estimated. Numerical calculations show that in the case of discrete auroral events the effect of the residual crustal magnetic field leads to a significant increase in the upward fluxes of electrons, which causes a decrease in the rates of heating and ionization of the atmospheric gas in comparison with the calculations without taking into account the residual magnetic field. It is shown that all the above-mentioned impact factors of auroral electron precipitation processes should be taken into account both in the photochemical models of the Martian atmosphere and in the interpretation of observations of the chemical composition and its variations using the ACS instrument onboard ExoMars.     

Intensity variation of atomic oxygen red line in the day airglow with solar activity

Rates of production of O(¹ D) atoms in the upper atmosphere by photodissociation of O₂, dissociative recombination of O₂ ⁺, NO⁺ and electron impact excitation of O(³ P) have been calculated for low, medium and high levels of solar activity. Variations with solar activity, of neutral and ionic composition, electron and neutral temperatures of the upper atmosphere and solar extreme ultraviolet fluxes incident on it have been taken into consideration.Emission rates ofOi red line (6300Å) have been computed taking into account the deactivation both by molecular oxygen and nitrogen. It has been shown that the integrated intensity from low to high activity period varies by approximately an order of magnitude in agreement with the results of experimental observations.     

A sensitive search for nitric oxide in the lower atmospheres of Venus and Mars: Detection on Venus and upper limit for Mars

High-resolution spectra of Venus and Mars at the NO fundamental band at 5.3 μm with resolving power ν/δν=76,000 were acquired using the TEXES spectrograph at NASA IRTF on Mauna Kea, Hawaii. The observed spectrum of Venus covered three NO lines of the P-branch. One of the lines is strongly contaminated, and the other two lines reveal NO in the lower atmosphere at a detection level of 9 sigma. A simple photochemical model for NO and N at 50-112 km was coupled with a radiative transfer code to simulate the observed equivalent widths of the NO and some CO₂ lines. The derived NO mixing ratio is 5.5±1.5 ppb below 60 km and its flux is . Predissociation of NO at the (0-0) 191 nm and (1-0) 183 nm bands of the δ-system and the reaction with N are the only important loss processes for NO in the lower atmosphere of Venus. The photochemical impact of the measured NO abundance is significant and should be taken into account in photochemical modeling of the Venus atmosphere. Lightning is the only known source of NO in the lower atmosphere of Venus, and the detection of NO is a convincing and independent proof of lightning on Venus. The required flux of NO is corrected for the production of NO and N by the cosmic ray ionization and corresponds to the lightning energy deposition of . For a flash energy on Venus similar to that on the Earth (∼¹⁰⁹ J), the global flashing rate is ∼90 s⁻¹ and ∼6 km⁻² y⁻¹ which is in reasonable agreement with the existing optical observations. The observed spectrum of Mars covered three NO lines of the R-branch. Two of these lines are contaminated by CO₂ lines, and the line at 1900.076 cm⁻¹ is clean and shows some excess over the continuum. Some photochemical reactions may result in a significant excitation of NO (v=1) in the lowest 20 km on Mars. However, quenching of NO (v=1) by CO₂ is very effective below 40 km. Excitation of NO (v=1) in the collisions with atomic oxygen is weak because of the low temperature in the martian atmosphere, and we do not see any explanation of a possible emission of NO at 5.3 μm. Therefore the data are treated as the lack of absorption with a 2 sigma upper limit of 1.7 ppb to the NO abundance in the lower atmosphere of Mars. This limit is above the predictions of photochemical models by a factor of 3.