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

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

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.     

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.     

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

Energetic Neutral Atoms (ENA) at Mars: Properties of the hydrogen atoms produced upstream of the martian bow shock and implications for ENA sounding technique around non-magnetized planets

We have studied the interaction of fast solar wind hydrogen atoms with the martian atmosphere by a three-dimensional Monte Carlo simulation. These energetic neutral hydrogen atoms, H-ENAs, are formed upstream of the martian bow shock. Both H-ENAs scattered and non-scattered from the martian atmosphere/exosphere were studied. The colliding H-ENAs were found to scatter both to the dayside and nightside. On the dayside they contribute to the so-called H-ENA albedo. On the nightside the heated and scattered hydrogen atoms were found also in the martian wake. The density, the energy distribution function and the direction of the velocity of H-ENAs on the nightside are presented. The present study describes a novel “ENA sounding” technique in which energetic neutral atoms are used to derive information of the properties of planetary exosphere and atmosphere in a similar manner as the solar wind photons are used to derive atmospheric densities by measuring the scattered UV light. A detailed study of the direction and energy of the scattered and non-scattered H-ENAs suggest that the ENA sounding is a method to study the interaction between the planetary atmosphere and the solar wind and to monitor the density, and likely also the magnetization, of the planetary upper atmosphere. Already present-day ENA instrument should be capable to detect the analyzed particle fluxes.     

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.     

Ulf Waves at the Proton Gyrofrequency Upstream of the Martian bow Shock, Based on Mars Global Surveyor's Measurements

Small amplitude MHD waves at the proton gyrofrequency were detected and analyzed using the magnetic field measurements obtained by the Mars Global Surveyor, in the region upstream of the Martian bow shock, between days 87 and 255 of 1998. The origin of the observed ULF waves was found to be protons back-streaming from the bow shock and pick-up protons from the planetary hydrogen exosphere, interacting with the incoming solar wind. A small contribution of oxygen ions has also been detected. Furthermore, the spatial distribution of the waves was investigated, using the spacecraft's position during the time the data were obtained.     

Neutral and ion-exospheres in the solar wind with applications to mercury

A model of planetary neutral and ion-exospheres in the solar wind is formulated for weak or lunar like solar-wind interaction with a planet. The neutral exosphere model allows for density and temperature variations and for rotation at the exobase. The ion-exosphere is produced by ionization of the neutral exosphere in the solar wind and its density distribution is obtained by solving the continuity equation in the drift approximation. Applying to Mercury a surface temperature distribution inferred from infra-red data and a vanishing bound neutral flux at the base, He and He⁺ density distributions are found. When the He atmosphere of Mercury is due entirely to the surface bombardment by solar wind He⁺⁺, the resulting He⁺ density is found to vary from 1.5 × 10⁻¹ to 10⁻³ cm⁻³ over the range 1.5–5 planetocentric radii on the dayside. These densities are found to be detectable by typical solar-wind plasma instruments. The possible effects of cyclotron-resonant scattering by interplanetary magnetic field fluctuations are examined and shown to be negligible. An electromagnetic plasma instability, triggered by the birth of ions in the exosphere, is shown to be important for the thermalization of the energy mode transverse to the interplanetary magnetic field, allowing more ions to be detected by solar-wind ion probes.     

Upstream proton cyclotron waves at Venus

Data from the magnetometer MAG aboard the Venus Express S/C are investigated for the occurrence of cyclotron wave phenomena upstream of the Venus bow shock. For an unmagnetized planet such as Venus and Mars the neutral exosphere extends into the on-flowing solar wind and pick-up processes can play an important role in the removal of particles from the atmosphere. At Mars upstream proton cyclotron waves were observed but at Venus they were not yet detected. From the MAG data of the first 4 months in orbit we report the occurrence of proton cyclotron waves well upstream from the planet, both outside and inside of the planetary foreshock region; pick-up protons generate specific cyclotron waves already far from the bow shock. This provides direct evidence that the solar wind is removing hydrogen from the Venus exosphere. Determining the role the solar wind plays in the escape of particles from the total planetary atmosphere is an important step towards understanding the evolution of the environmental conditions on Venus. The continual observations of the Venus Express mission will allow mapping the volume of escape more accurately, and determine better the present rate of hydrogen loss.     

History of Martian volatiles: Implications for organic synthesis

A theoretical reconstruction of the history of Martian volatiles indicates that Mars probably possessed a substantial reducing atmosphere at the outset of its history and that its present tenous and more oxidized atmosphere is the result of extensive chemical evolution. As a consequence, it is probable that Martian atmospheric chemical conditions, now hostile with respect to abiotic organic synthesis in the gas phase, were initially favorable. Evidence indicating the chronology and degradational history of Martian surface features, surface mineralogy, bulk volatile content, internal mass distribution, and thermal history suggests that Mars catastrophically developed a substantial reducing atmosphere as the result of rapid accretion. This atmosphere probably persisted—despite the direct and indirect effects of hydrogen escape—for a geologically short time interval during, and immediately following, Martian accretion. That was the only portion of Martian history when the atmospheric environment could have been chemically suited for organic synthesis in the gas phase. Subsequent gradual degrassing of the Martian interior throughout Martian history could not sustain a reducing atmosphere due to the low intensity of planet-wide orogenic activity and the short atmospheric mean residence time of hydrogen on Mars. During the post-accretion history of Mars, the combined effects of planetary hydrogen escape, solar-wind sweeping, and reincorporation of volatiles into the Martian surface produced and maintained the present atmosphere.     

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.     

Kinetic models for the exospheres of Jupiter and Saturn

The exospheric theory based on the Kappa velocity distribution function (VDF) is used to model the exosphere of the giant planets Jupiter and Saturn. Such Kappa velocity distribution functions with an enhanced population of suprathermal particles are indeed often observed in space plasmas and in the space environment of the planets. The suprathermal particles have significant effects on the escape flux, density and temperature profiles of the particles in the exosphere of the giant planets. The polar wind flux becomes several orders larger when suprathermal electrons are considered, so that the planetary ionosphere becomes then a significant source for their inner magnetosphere. Moreover, the number density of the particles decreases slower as a function of the altitude when a Kappa distribution is considered instead of a Maxwellian one. Two-dimensional maps of density are calculated for typical values of the temperatures. The exospheric formalism is also applied to study the escape flux from the exospheres of Io and Titan, respectively, moons of Jupiter and Saturn.     

Monte Carlo study of interaction between solar wind plasma and Venusian upper atmosphere

Isolated events of proton and alpha particle precipitation in the Venusian atmosphere were recorded with the use of the ASPERA-4 analyzer on board the ESA Venus Express spacecraft. Using a Monte Carlo simulation method for calculation of proton and alpha particle precipitations in the Venusian atmosphere, reflected and upward directed particle fluxes have been found. It has been found that only a vanishing percentage of protons and alpha particles are backscattered to the Venusian exosphere when neglecting the induced magnetic field and under conditions of low solar activity. Accounting for the induced field drastically changes the situation: the backscattered by the atmosphere energy fluxes increase up to 44% for the horizontal magnetic field B = 20 nT, measured for Venus, for the case of precipitating protons, and up to 64%, for alpha particles. The reflected energy fluxes increase to about 100% for both protons and alpha particles as the field grows to 40 nT, i.e., the atmosphere is protected against penetration of solar wind particles.     

Non-thermal water loss of the early Mars: 3D multi-ion hybrid simulations

In this study we analyze the non-thermal loss rates of O⁺, O₂⁺ and CO₂⁺ ions over the last 4.5 billion years (Gyr) in the Martian history by using a 3D hybrid model. For this reason we derived the past solar wind conditions in detail. We take into account the intensified particle flux of the early Sun as well as an Martian atmosphere, which was exposed to a sun's extreme ultraviolet (EUV) radiation flux 4.5 Gyr ago that was 100 times stronger than today. Furthermore, we model the evolution of the interplanetary magnetic field by a Weber & Davis solar wind model. The ‘external’ influences of the Sun's radiation flux and solar wind flux lead to the formation of an ionospheric obstacle by photoionization, charge exchange and electron impact. For the early Martian conditions we could show that charge exchange was the dominant ionization mechanism. Several hybrid simulations for different stages in the evolution of the Martian atmosphere, at 1, 2, 5, 10, 30 and 100 EUV, were performed to analyze the non-thermal escape processes by ion pick-up, momentum transfer from the solar wind to the ionosphere and detached ionospheric plasma clouds. Our results show a non-linear evolution of the loss rates. Using mean solar wind parameters the simulations result in an oxygen loss equivalent to the depth of a global Martian ocean of about 2.6 m over the last 4.5 Gyr. The induced magnetic field strength could be increased up to about 2000 nT. A simulation run with high solar wind density results in an oxygen loss of a Martian ocean up to 205 m depth during 150 million years after the sun reached the zero age mean sequence (ZAMS).