Effects of charge exchange involving H and H+ in the upper atmosphere

It is argued that there is a terrestrial loss of hydrogen as ions which includes the polar wind but extends effectively down to a latitude in the range 45–50° invariant. In daytime and for much of the night-time the flux is close to the limiting value for H⁺ flow through the topside ionosphere. It is argued that the flux decreases rapidly with increasing solar activity, following the decrease in neutral hydrogen concentration. It has been found that as solar activity increases the Jeans escape flux increases, and the charge exchange escape flux increases until moderate solar activity levels are reached. As solar activity increases from moderate to high levels, the charge exchange escape may decrease again. A new budget for terrestrial hydrogen loss over the solar cycle is given. The global flux of hydrogen ions outward from the ionosphere is comparable with estimates of the plasma sheet loss rates, and this flux, together with some solar wind plasma, is an attractive source for the plasma sheet.The energetic neutrals produced from the charge exchange of ring current ions with thermal-energy neutrals in the exosphere produce the optical emission of the equatorial aurora, which can be related to ion production rates near and above the E-region. The ionization production is adequate to explain the enhancements in ion production observed during magnetic storms at Arecibo.     

A model of the electron and ion temperatures in the ionosphere

A model (empirical) of the electron and ion temperatures (T_eT_i) is presented in the altitude interval 50–4000 km as a function of time (diurnal, annual), space (position, altitude) and solar flux (F_10.7). Using observations of six satellites (AE-C, AE-D, AE-E, ISIS-1, ISIS-2, OGO-6), five incoherent scatter stations (Arecibo, Chatanika, Jicamarca, Millstone Hill, St Santin) and rocket measurements, this model describes the global gross features of the ionosphere during quiet geophysical conditions (K_p⩽3). The numerical analysis is based on spectral decomposition; the horizontal structure is represented by spherical harmonies and Fourier series, and the vertical structure by spline functions. The electron temperature is, in general, very similar to the ion temperature below ∼90 km. Up to approx. 1500 km, the electron temperature is, on an average, distinctly higher than the ion temperature. Above ∼2000 km, however, the ion temperature is quickly catching up and attains somewhat below 4000 km the same magnitude as the electron temperature.     

Nonthermal escape of hydrogen and deuterium from Venus and implications for loss of water

Nonthermal escape processes responsible for the escape of hydrogen and deuterium from Venus are examined for present and past atmospheres. Three mechanisms are important for the escape of hydrogen from the present atmosphere: (a) charge exchange of plasmaspheric H⁺ with exospheric H, (b) impact of exospheric hot O atoms on H, and (c) ion molecule reactions involving O⁺ and H₂. However, in the past when the H abundance was higher, the charge-exchange mechanism would be the strongest. The H escape flux increases rapidly with increasing hydrogen abundance in the upper atmosphere and saturates at a value of 1 × 10¹⁰ cm^?2 sec^?1 emerging primarily from the day side when the H mixing ratio at the homopause is 2 × 10^?3. This corresponds to an H₂O mixing ratio of 1 × 10^?3 at the cold trap and ～15% at the surface. Deuterium would also escape by the charge-exchange mechanism and a D/H enrichment by a factor of ～1000 over the nonthermal escape regime is expected, which could have lasted over the last 3 billion years. Coincidentally, the onset of hydrodynamic flow leading to efficient H escape occurs just at the H₂O mixing ratio at which the charge-exchange escape flux saturates. Thus it is possible that Venus has lost an Earth-equivalent ocean of water over geologic time. If so, either the D/H enrichment has been kept low by modest outgassing of juvenile water or Venus started out with a D/H ratio of ～4.0 × 10^?6.     

Helium escape from the Earth''s atmosphere: The charge exchange mechanism revisited

We have studied the escape of neutral helium from the terrestrial atmosphere through exothermic charge exchange reactions between He⁺ ions and the major atmospheric constituents N₂, O₂, and O. Elastic collisions with the neutral background particles were treated quantitatively using a recently developed kinetic theory approach. An interhemispheric plasma transport model was employed to provide a global distribution of He⁺ ions as a function of altitude, latitude and local solar time and for different levels of solar ionization. Combining these ion densities with neutral densities from an MSIS model and best estimates for the reaction rate coefficients of the charge exchange reactions, we computed the global distribution of the neutral He escape flux. The escape rates show large diurnal and latitudinal variations, while the global average does not vary by more than a factor of three over a solar cycle. We find that this escape mechanism is potentially important for the overall balance of helium in the Earth's atmosphere. However, more accurate values for the reaction rate coefficients of the charge exchange reactions are required to make a definitive assessment of its importance.     

Auroral zone effects on hydrogen geocorona structure and variability

The instantaneous structure of planetary exospheres is determined by the time history of energy dissipation, chemical, and transport processes operative during a prior time interval set by intrinsic atmospheric time scales. The complex combination of diurnal and magnetospheric activity modulations imposed on the Earth's upper atmosphere no doubt produce an equally complex response, especially in hydrogen, which escapes continuously at exospheric temperatures. Vidal-Madjar and Thomas (1978) have discussed some of the persistent large scale structure which is evident in satellite ultraviolet observations of hydrogen, noting in particular a depletion at high latitudes which is further discussed by Thomas and Vidal-Madjar (1978). The latter authors discussed various causes of the H density depletion, including local neutral temperature enhancements and enhanced escape rates due to polar wind H⁺ plasma flow or high latitude ion heating followed by charge exchange. We have reexamined the enhancement of neutral escape by plasma effects including the recently observed phenomenon of low altitude transverse ion acceleration. We find that, while significant fluxes of neutral H should be produced by this phenomenon in the auroral zone, this process is probably insufficient to account for the observed polar depletion. Instead, the recent exospheric temperature measurements from the Dynamics Explorer-2 spacecraft suggest that neutral heating in and near the high latitude cusp may be the major contributor to depleted atomic hydrogen densities at high latitudes.     

Fractionation of hydrogen and deuterium on Venus due to collisional ejection

The fractionation factor f is important for interpreting the current escape fluxes of H and D on Venus and how the D/H ratio has evolved. The escape flux is currently governed by the two processes of charge exchange and collisional ejection by fast oxygen atoms. Using a best-fit parameterized equation for the O-H scattering angle phase function, more accurate branching ratios for the oxygen ion dissociation and including the effects of the initial energy and momentum of the ions and electrons, as well as for the hydrogen and deuterium gas, we have reanalyzed the collisional ejection process. Our analysis produces improved values for the efficiency of H and D escape as a function of the ionospheric temperature. From our results we propose the reduction of the hydrogen flux for collisional ejection from 8 to 3.5 x 10(6) cm-2 s-1. Assuming that collisions leading to escape occur mostly in the region between 200 and 400 km, the revised D/H fractionation factor due to collisional ejection is 0.47, where previously the process had been considered completely discriminating against deuterium escape (or f approximately 0.) The resulting deuterium flux is 3.1 x 10(4) cm-2 s-1, roughly 6 times the flux due to charge exchange, making collisional ejection the dominant escape mechanism for deuterium on Venus.     

Global atomic hydrogen density derived from OGO-6 Lyman α measurements

Measurements of the Lyman α airglow intensity were made between June 1969 and June 1970 by a u.v. photometer experiment on the OGO-6 satellite. The data for the zenith intensity at altitudes between 400 and 1100 km were fitted to theoretical airglow models to derive atomic hydrogen density n_c at a reference altitude, taken to be 500 km. n_c was determined for each of 286 orbits throughout the year. The mean exospheric temperature T∞(J) during this period varied from 900 to 1300 K according to the Jacchia model. The solar Lyman α flux at line center F₀ was also determined over each 90-min orbit in the model-fitting procedure. F₀ was found to be correlated with sunspot number, in agreement with previous results. A nearly-exact linear relationship was found for F₀, when averaged over ‘bins’ which are 20 sunspot numbers in width. n_c was found to be inversely correlated with T∞(J); however the dependence is not that predicted by steady-state models whose only escape mechanism is Jeans evaporative escape. Unless the total atmospheric loss rate depends upon 27-day changes in the solar EUV, which is unlikely, an additional upper atmospheric loss is required in order that the total loss remain constant with T∞(J). This extra loss may be largely due to charge-exchange reactions in the exosphere, wherein energetic protons are converted to fast hydrogen atoms, as suggested previously by a number of authors. An additional result is suggested by the apparent spherical symmetry of the inferred density, namely that the familiar diurnal variation of hydrogen is absent at the high latitudes preferentially sampled by the OGO-6 data.     

The diurnal variation of atomic hydrogen

This note critically examines the relative importance of several effects which influence the diurnal variation of atomic hydrogen abundance near the critical level.It is pointed out that the neglect of exospheric hydrogen in a recent theoretical treatment causes an overestimation of the diurnal variation at high exospheric temperatures, and an underestimation at low exospheric temperatures. The fluxes due to lateral flow are large compared to other fluxes only to the extent that the actual diurnal variation is very different from the diurnal variation corresponding to zero net lateral flow, which does not seem to be the case in the real atmosphere. Two effects which are probably important are charge exchange reactions with thermal oxygen ions, resulting in a diurnal exchange with the plasmasphere; and charge exchange reactions with high velocity protons, resulting in enhanced escape and diurnal variation.     

A model of thermospheric temperature and composition

An empirical model of thermospheric temperature (T_∞T₁₂₀, and s) and composition (H, He, N, O, N₂, O₂, and Ar) was derived from measurements of 8 satellites (AE-C, AE-E, AEROS-A, AEROS-B, ARIEL-3, ESRO-4, OGO-6, and SAN MARCO-3) and 4 incoherent scatter stations (Arecibo, Jicamarca, Millstone Hill, and St Santin). The altitude covered extends from 120 km up to about 600 km over the time period 1967 to 1976. The analytical framework used in the model resembles closely the MSIS setup: time independent terms, solar flux terms, geomagnetic activity (K_p) effect, annual (semiannual) and diurnal (semidiurnal, terdiurnal) variations, longitudinal terms, the U.T. effect, and corrections compensating for deviations from diffusive equilibrium at altitudes below 200 km. The model describes quiet to medium disturbed geomagnetic conditions (K_p ? 4) at solar fluxes (10.7cm) ranging from 60 to 180 × 10^?22 Wm^?2Hz^?1. To get an impression of the accuracy presently obtained, the model is compared with MSIS, Jacchia (1977), and the models of Thuillier (T_∞ and Engebretson (N). The best agreement is found for the temperature and the constituents He, O, and N₂ with increasing deviations in the order of H, N, Ar, and O₂.     

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

The regulation of hydrogen and oxygen escape from Mars

It is shown that under present conditions the Jeans escape flux of hydrogen from Mars in the form of H and H₂ is constrained to be the same as twice the non-thermal (McElroy, 1972) escape of O atoms. The mediation of the chemical chain that recombines CO₂ plays an essential role in regulating the escape of hydrogen to match that of oxygen, confirming a mechanism postulated by McElroy and Donahue (1972). It is also shown that if the oxygen flux changes, a change in the O₂ mixing ratio results and the consequence is to induce a large change in the odd hydrogen concentration, and consequently in H₂ production and hydrogen escape. The effect is to stabilize the hydrogen escape flux at twice the O flux. It is shown that surface chemistry should not change the operation of this mechanism but has consequences for the eddy coefficient variation at low altitudes. There is a strong correlation between low humidity, large solar zenith angles and large O₃ abundances. The effect of argon in a mixing ratio as large as 0.3 on these results is also investigated.     

N2 escape rates from Pluto's atmosphere

Hydrodynamic escape of N₂ molecules from Pluto's atmosphere is calculated under the assumption of a high density, slow outflow expansion driven by solar EUV heating by N₂ absorption, near-IR and UV heating by CH₄ absorption, and CO cooling by rotational line emission as a function of solar activity. At 30 AU, the N₂ escape rate varies from in the absence of heating, but driven by an upward thermal heat conduction flux from the stratosphere, for lower boundary temperatures varying from 70-100 K. With solar heating varying from solar minimum to solar maximum conditions and a calculated lower boundary temperature, 88.2 K, the N₂ escape rate range is , respectively. LTE rotational line emission by CO reduces the net solar heat input by at most 35% and plays a minor role in lowering the calculated escape rates, but ensures that the lower boundary temperature can be calculated by radiative equilibrium with near-IR CH₄ heating. While an upward thermal conduction heat flux at the lower boundary plays a fundamental role in the absence of heating, with solar heating it is downward at solar minimum, and is, at most, 13% of the integrated net heating rate over the range of solar activity. For the arrival of the New Horizons spacecraft at Pluto in July 2015, predictions are lower boundary temperature, T₀∼81 K, and N₂ escape rate , and peak thermospheric temperature ∼103 K at 1890 km, based on expected solar medium conditions.     

Transterminator ion flow in the Martian ionosphere

The upper ionospheres of Mars and Venus are permeated by the magnetic fields induced by the solar wind. It is a long-standing question whether these fields can put the dense ionospheric plasma into motion. If so, the transterminator flow of the upper ionosphere could explain a significant part of the ion escape from the planets atmospheres. But it has been technically very challenging to measure the ion flow at energies below 20 eV. The only such measurements have been made by the ORPA instrument of the Pioneer Venus Orbiter reporting speeds of 1-5 km/s for O⁺ ions at Venus above 300 km altitude at the terminator ( [Knudsen et al., 1980] and [Knudsen et al., 1982]). At Venus the transterminator flow is sufficient to sustain a permanent nightside ionosphere, at Mars a nightside ionosphere is observed only sporadically. We here report on new measurements of the transterminator ion flow at Mars by the ASPERA-3 experiment on board Mars Express with support from the MARSIS radar experiment for some orbits with fortunate observation geometry. We observe a transterminator flow of O⁺ and O₂⁺ ions with a super-sonic velocity of around 5 km/s and fluxes of 0.8×10⁹/cm² s. If we assume a symmetric flux around the terminator this corresponds to an ion flow of 3.1±0.5×10²⁵/s half of which is expected to escape from the planet. This escape flux is significantly higher than previously observed on the tailside of Mars. A possible mechanism to generate this flux can be the ionospheric pressure gradient between dayside and nightside or momentum transfer from the solar wind via the induced magnetic field since the flow velocity is in the Alfvénic regime. We discuss the implication of these new observations for ion escape and possible extensions of the analysis to dayside observations which may allow us to infer the flow structure imposed by the induced magnetic field.     

Monte Carlo computations of the escape of atomic nitrogen from Mars

Monte Carlo simulations were carried out to compute the escape flux of atomic nitrogen for the low and high solar activity martian thermospheres. The total escape of atomic nitrogen at low and high solar activities was found to be 3.03×¹⁰⁵ and , respectively. The escape flux of atomic nitrogen at low and high solar activities from photodissociation of N₂ was found to be 2.75×¹⁰⁵ and , respectively. The remainder of the contribution is from dissociative recombination, which is only important at high solar activity were it comprises about 25% of the total escape. The relative contributions to the total N escape flux from thermal motion of the background atmosphere, winds and co-rotation, and photoionization and subsequent solar wind pickup are also considered here. We find that the total predicted escape fluxes are observed to increase by 20 and 25% at low and high solar activities owing to thermal motion of the background atmosphere. At low and high solar activities, we find that the co-rotation and wind velocities combined translate to a maximum transferable energy of ∼0.0103 and 0.0181 eV, respectively, and that the total escape flux contribution from winds and co-rotation is negligible. Photoionization was found to be a minor process only impacting those source atoms produced with energies close to the escape energy, between 1.5 and 2 eV. The contributions to the total escape fluxes at low and high solar activities from photoionization and subsequent solar wind pickup are found to be about 8 and 13%, respectively.     

Venus O+ pickup ions: Collected PVO results and expectations for Venus Express

Observations of oxygen pickup ions by the plasma analyzer on the Pioneer Venus Orbiter (PVO) Mission arguably launched broad interest in solar wind erosion of unmagnetized planet atmospheres, and its potential evolutionary effects. Oxygen pickup ions may play key roles in the removal of the oxygen excess left behind from the photodissociation of water vapor by enabling direct escape, additional sputtering of oxygen when they impact the exobase, and escape as energetic neutrals produced in charge exchange reactions with the ambient exospheric oxygen and hydrogen. Although the PVO observations were compromised by an ∼8 keV energy limit for O⁺ detection, a lack of ion composition capability, and the limited sampling and data rate of the plasma analyzer which was designed for solar wind monitoring, these measurements provide our best information about the extended O⁺ exosphere and wake at Venus. Here we show the full picture of the spatial distribution and energies of the O⁺ ion observations collected by the plasma analyzer during PVO's ∼5000 orbit tour. A model of O⁺ test particles launched in the circum-Venus fields described by an MHD simulation of the solar wind interaction is used to help interpret the PVO observations and to anticipate the expanded view of Venus O⁺ escape that will be provided by the ASPERA-4 experiment on Venus Express.     

The aeronomy of the upper atmosphere of Venus

Density profiles for CO, O, and O₂ in the Cytherean atmosphere above 90 km are plotted with eddy diffusion coefficient (K) as a parameter, subject to the constraint that the mixing ratios of CO and O₂ approach their observed value or values under the observed upper limit at the lower boundary. It is then shown that the value of K puts upper limits on the amount of hydrogen (in the form of H₂O, HCl, and H₂) the atmosphere near 90km can contain. This value is a function of the density and temperature of hydrogen at the critical level and the magnitude of the total escape flux, where unspecified flux mechanisms other than thermal are postulated ad hoc. In general these constraints call for large values of K to accomodate the atomic hydrogen produced by measured mixing ratios of HCl and H₂O. Hence they constrain thee amount of O in the upper atmosphere to values well under 1% at 130 km unless there are very large hydrogen escape fluxes, 10⁷ cm^?2sec^?1 or larger. The freedom to assume arbitrary amounts of H₂ in the atmosphere is also restricted. We suggest either very effective escape mechanisms—despite low exospheric hydrogen densities—or novel excitation mechanisms for O(3³S) and O(3⁵S) in the upper atmosphere.     

Backscatter of solar resonance radiation—I

A simple model which essentially neglects temperature effects, but which includes gravity, radiation pressure, photoionization and charge exchange, is used to calculate the angular dependence of the intensity of solar Lyman α resonantly scattered from neutral interstellar hydrogen which has penetrated the solar system; the results are then compared with the observations. The resonant scattering of HeI λ584 is also treated.     

A model of the terrestrial ionosphere in the altitude interval 50–4000 km II. Molecular ions (N2 +, NO+, O2 +) and electron density

Empirical models of molecular ion densities (N₂ ⁺, NO⁺, O₂ ⁺) and the electron density (N _e ) are presented in the altitude interval 50–4000 km as functions of time (diurnal, annual), space (position, altitude) and solar flux (F _10.7). Using observations of 6 satellites (AE-C, AE-D, AE-E, ALOUETTE-2, ISIS-1, ISIS-2), 4 incoherent scatter stations (Arecibo, Jicamarca, Millstone Hill, St Santin) and more than 700 D-region profiles, this model describes the global gross features of the ionosphere for quiet geophysical conditions (K _p 3).The molecular ion densities and the electron density increase with increasing altitude up to a maximum (or several maxima) - and decrease from thereon with increasing height. Between ~80 and 200 km, the main ionic constituents are NO⁺ and O₂ ⁺; below ~80 km cluster ions are predominating. During local summer conditions the molecular ions and N _e increase around polar latitudes and decrease correspondingly during local winter. The diurnal variations are intrinsically coupled to the individual plasma layers; in general, the molecular ion and electron densities are enhanced during daytime and depleted during nighttime (for details and exceptions, see text).     

Loss of water from Venus. I. Hydrodynamic escape of hydrogen

A one-dimensional photochemical-dynamic model is used to study hydrodynamic loss of hydrogen from a primitive, water-rich atmosphere on Venus. The escape flux is calculated as a function of the H₂O mixing ratio at the atmospheric cold trap. The cold-trap mixing ratio is then related in an approximate fashion to the H₂O concentration in the lower atmosphere. Hydrodynamic escape should have been the dominant loss process for hydrogen when the H₂O mass mixing ratio in the lower atmosphere exceeded ～0.1. The escape rate would have depended upon the magnitude of the solar ultraviolet flux and the atmospheric euv heating efficiency and, to a lesser extent, on the O₂ content of the atmosphere. The time required for Venus to have lost the bulk of a terrestrial ocean of water is on the order of a billion years. Deutrium would have been swept away along with hydrogen if the escape rate was high enough, but some D/H enrichment should have occurred as the escape rate slowed down.     

The intensity variation of the atomic oxygen red line during morning and evening twilight on 9–10 April 1969

During the evening of 9 April and the morning of 10 April 1969, the twilight zenith intensity of the atomic oxygen red line OI(³P-¹D) at 6300 Å was measured at the Blue Hill Observatory (42°N, 17°W). At the same time incoherent scatter radar data were being obtained at the Millstone Hill radar site 50 km distant. We have used a diurnal model of the mid-latitude F-region to calculate the ionospheric structure over Millstone Hill conditions similar to 9–10 April 1969. The measured electron temperature, ion temperature, and electron density at 800 km are used as boundary conditions for the model calculations. The diurnal variation of neutral composition and temperature were obtained from the OGO-6 empirical model and the neutral winds were derived from a semiempirical three-dimensional dynamic model of the neutral thermosphere. The solar EUV flux was adjusted to yield reasonable agreement between the calculated and observed ionospheric properties.This paper presents the results of these model computations and calculations of the red line intensity. The 6300 Å emission includes contributions from photoelectron excitation, dissociative recombination, Schumann-Runge photodissociation and thermal electron impact. The variations of these four components for morning and evening twilight between 90–120° solar zenith angles, and their relative contributions to the total 6300 Å emission line intensity, are presented and the total is compared to the observations. For this particular day the Schumann-Runge photodissociation component, calculated using the solar fluxes tabulated by Ackermann (1970), is the dominant component of the morning twilight 6300 Å emission. During evening twilight it is necessary to utilize a lower O₂ density than for the morning twilight in order to bring the calculated and observed 6300 Å emission rates into agreement. The implication that there may be a diurnal variation in the O₂ density at the base of the thermosphere is discussed in the light of available experimental data and current theoretical ideas.