A model of the terrestrial ionosphere in the altitude interval 50–4000 km I. Atomic ions (H+, He+, N+, O+)

An empirical model of atomic ion densities (H⁺, He⁺, N⁺, O⁺) is presented up to 4000 km altitude as a function of time (diurnal, annual), space (position, altitude) and solar flux (F10.7) — using observations of satellites (AE-B, AE-C, AE-D, AE-E, ISIS-2, OGO-6) and rockets during quiet geophysical conditions (K _p 3). The numerical treatment is based upon harmonic functions for the horizontal pattern and cubic splines for the vertical structure.The ion densities increase with increasing height up to a maximum (depending roughly on the ion mass) and decrease beyond that with increasing altitude. Above 200 km, O⁺ is the main ionic constituent being replaced at approximately 800 km (depending on latitude, local time, etc.) by H⁺. Around polar regions the light ions, H⁺ and He⁺, are depleted (polar wind) and the heavier ones enhanced. During local summer conditions the ion densities increase around polar latitudes and correspondingly decrease during local winter, except He⁺ which reflects the opposite pattern. Diurnal variations are intrinsically coupled to the individual plasma layers: N⁺ and O⁺ peak, in general, during daytime, while the amplitudes and phases of H⁺ and He⁺ change strongly with altitude and latitude. Earth, Moon and Planets Review article.     

The topside equatorial ionosphere during daytime: Elevated plasma temperatures,the transequatorial O+ breeze and the dynamics of minor ions

Steady-state calculations are performed for the daytime equatorial F2-region and topside ionosphere. Values are calculated of the electron and ion temperatures and the concentrations and field-aligned velocities of the ions O⁺, H⁺ and He⁺. Account is taken of upward $\text{E × B}$ drift, a summer-winter horizontal neutral air wind and heating of the electron gas by thermalization of fast photoelectrons.The calculated plasma temperatures are in accord with experiment: at the equator there is an isothermal region from about 400–550 km altitude, with temperatures of about 2400 K around 800 km altitude. The transequatorial O⁺ breeze flux from summer to winter in the topside ionosphere is not greatly affected by the elevated plasma temperatures. The field-aligned velocities of H⁺ and He⁺ depend strongly on the O⁺ field-aligned velocity and on the presence of large temperature gradients. For the minor ions, ion-ion drag with O⁺ cannot be neglected for the topside ionosphere.     

Helium ion outflow from the terrestrial ionosphere

Extensive calculations have been made of the behaviour of He⁺ for situations where ion outflow occurs from the topside ionosphere. For these circumstances, steady state solutions for the He⁺ continuity, momentum and energy equations have been obtained self-consistently, yielding density, velocity and temperature profiles of He⁺ from 200 to 2000 km altitude. To model the high latitude topside ionosphere, a range of background H⁺O⁺ ionospheres was considered with variations in the H⁺ outflow velocity, the presence of a perpendicular electric field and different peak O⁺ densities. In addition, the atmospheric density of neutral helium was chosen to model typical observed winter and summer densities. From our studies we have found that: (a) The outflowing He⁺ has density profiles of similar shape to those of H⁺, for basically different reasons; (b) The effect of the perpendicular electric field differs considerably for H⁺ and He⁺. This difference stems from the fact that an electric field acts to alter significantly the O⁺ density at high altitudes and this, in turn, changes the H⁺ escape flux through the O⁺+H charge exchange reaction. A similar situation does not occur for He⁺ and therefore the He⁺ escape flux exhibits a negligibly small change with electric field; (c) The fractional heating of He⁺ due to the He⁺O⁺ relative flow is not as effective in heating He⁺ as the H⁺O⁺ relative flow is in heating H⁺; (d) During magnetospheric disturbances when the N₂ density at the altitude of the He⁺ peak (600 km) can increase by a factor as large as 50, the He⁺ peak density decreases only by approximately a factor of 2; and (e) The He⁺ escape flux over the winter pole is approximately a factor of 20 greater than the He⁺ escape flux over the summer pole. Consequently, on high latitude closed field lines there could be an interhemispheric He⁺ flux from winter to summer.     

Modeling of global variations and ring shadowing in Saturn's ionosphere

A time-dependent one-dimensional model of Saturn's ionosphere has been developed as an intermediate step towards a fully coupled Saturn Thermosphere-Ionosphere Model (STIM). A global circulation model (GCM) of the thermosphere provides the latitude and local time dependent neutral atmosphere, from which a globally varying ionosphere is calculated. Four ion species are used (H⁺, H⁺₂, H⁺₃, and He⁺) with current cross-sections and reaction rates, and the SOLAR2000 model for the Sun's irradiance. Occultation data from the Voyager photopolarimeter system (PPS) are adapted to model the radial profile of the ultraviolet (UV) optical depth of the rings. Diurnal electron density peak values and heights are generated for all latitudes and two seasons under solar minimum and solar maximum conditions, both with and without shadowing from the rings. Saturn's lower ionosphere is shown to be in photochemical equilibrium, whereas diffusive processes are important in the topside. In agreement with previous 1-D models, the ionosphere is dominated by H⁺ and H⁺₃, with a peak electron density of ∼¹⁰⁴ electrons cm⁻³. At low- and mid-latitudes, H⁺ is the dominant ion, and the electron density exhibits a diurnal maximum during the mid-afternoon. At higher latitudes and shadowed latitudes (smaller ionizing fluxes), the diurnal maximum retreats towards noon, and the ratio of [H⁺]/[H⁺₃] decreases, with H⁺₃ becoming the dominant ion at altitudes near the peak (∼1200-1600 km) for noon-time hours. Shadowing from the rings leads to attenuation of solar flux, the magnitude and latitudinal structure of which is seasonal. During solstice, the season for the Cassini spacecraft's encounter with Saturn, attenuation has a maximum of two orders of magnitude, causing a reduction in modeled peak electron densities and total electron column contents by as much as a factor of three. Calculations are performed that explore the parameter space for charge-exchange reactions of H⁺ with vibrationally excited H₂, and for different influxes of H₂O, resulting in a maximum diurnal variation in electron density much weaker than the diurnal variations inferred from Voyager's Saturn Electrostatic Discharge (SED) measurements. Peak values of height-integrated Pedersen conductivities at high latitudes during solar maximum are modeled to be ∼42 mho in the summer hemisphere during solstice and ∼18 mho during equinox, indicating that even without ionization produced by auroral processes, magnetosphere-ionosphere coupling can be highly variable.     

Regions of He+ dominance in the high-latitude topside ionosphere

Observations of the occurrence of He⁺ dominance in the topside ionosphere are discussed. An earlier model of the behaviour of high-latitude H⁺ and O⁺ thermal plasma (Quegan et al., 1982) is extended to include He⁺ as a major ion. Calculations using the extended model show that plasma convection is likely to play a key rôle in producing regions of He⁺ dominance. Suggested conditions for He⁺ dominance are listed and their applicability to observed He⁺ behaviour is discussed.     

Effects of solar XUV variations on the thermosphere and the variability of the photoionization efficiency parameter during the solar cycle 21

Making use of the latest available semi-empirical atmospheric models, solar XUV radiations rates of photoionization and absorbed energy profiles have been graphically presented showing the latitudinal, seasonal and solar cycle variations. The photoionization limits of the major neutral constitutents of the terrestrial atmosphere O₂, O, and N₂ that occur at wavelengths 102.7, 91.2, and 79.6 nm, respectively have been quantified by showing the photoionization rates of O ₂ ⁺ , O⁺, and N ₂ ⁺ for different spectral groups both under quiet and different solar flare conditions. The variability of the photoionization efficiency parameter which is height-dependent, from winter to summer, for solar minimum to solar maximum for four significantly different latitudes under local noon conditions have been investigated during the solar cycle 21. More energy is required to produce an electron-ion pair in a denser atmosphere than in a thinner atmosphere and hence more energy is being deposited in the height range between 100–120 km which itself manifests in raising the electron gas temperatures higher than the neutral gas temperatures.     

Diffusive equilibrium distributions of He+ in the plasmasphere

Recent satellite observations of thermal ion composition in the near-equatorial plasmasphere have shown that He⁺ comprises 5–10% typically and occasionally 25% or more of the total thermal ion density. A steady state diffusive equilibrium model for the distributions of H⁺, He⁺ and O⁺ along a plasmaspheric flux tube is used to elicit effects that may help explain these observed high He⁺ fractional concentrations. The model indicates that both the ionospheric composition and the temperature distribution along the flux tubes are important factors controlling the equatorial He⁺ composition, through the plasma scale height and thermal diffusion effects. Direct comparison of the model results with thermal ion observations by ISEE-1 indicates that the effects incorporated into the model may explain some of the elevated He⁺ concentrations. In some instances, however, effects not included in the model may also be of importance.     

Relative flow of H+ and O+ ions in the topside ionosphere at mid-latitudes

Theoretical results on the daily variation of O⁺ and H⁺ field-aligned velocities in the topside ionosphere are presented. The results are for an L = 3 magnetic field tube under sunspot minimum conditions at equinox. They come from calculations of time-dependent O⁺ and H⁺ continuity and momentum balance in a magnetic field tube which extends from the lower F2 region to the equatorial plane (Murphy et al., 1976).There are occasions when ion counterstreaming occurs, with the O⁺ velocity upward and H⁺ velocity downward. The conditions causing this counterstreaming are described: the H⁺ layer is descending whilst O⁺ is supplied from below either to increase the O⁺ concentration at fixed heights or to replace O⁺ ions lost by charge exchange with neutral H. It is suggested that the results of observations at Arecibo by Vickrey et al. (1976) of O⁺ and H⁺ concentrations and counterstreaming velocities are significantly affected by $\text{E}\text{×}\text{B}$ drift.     

Constraints on bulk composition,seasonal variation,and global dynamics of Pluto's atmosphere

Cosmic abundance, vapor pressure, and molecular weight considerations restrict the likely gas candidates for Pluto's atmosphere to Ne, N₂, CO, O₂, and Ar, in addition to the already detected CH₄. The vapor pressures and cosmic abundances of these gases indicate that all except Ne should be saturated in Pluto's atmosphere. The vapor pressure of Ne is so high that the existence of solid or liquid Ne on Pluto's surface is very unlikely; cosmic abundance arguments imply that Ne cannot attain saturation in Pluto's atmosphere. At both perihelion, N₂ should dominate the saturated gases. CO₂ should have the next highest mixing ratio, followed by O₂ and Ar. CH₄ should have the smallest mixing ratio. Because vapor pressures of these gases vary with temperature at diverse rates, the bulk and constituent mixing ratios of Pluto's atmosphere should vary with season. Between perihelion and aphelion, the column abundance of CH₄ may change by a factor of 260 while that of N₂ changes by only a factor of 52. The potential seasonal variation of Pluto's atmosphere was investigated by considering the behavior of these gases when individually mixed with CH₄. The effects of diurnal and latitudinal variation of insolation and eclipses on the atmosphere also were investigated. Seasonal effects are shown to dominate. It was shown that the atmospheric bulk may not be a minimum near aphelion but rather at intermediate distances from the Sun during summer/winter inadequate ice deposits may allow the atmosphere to collapse by freezing out over winter latitudes. If the atmosphere does not collapse, its weight is sufficient to keep it distributed uniformly around Pluto's surface. In this case, the atmosphere tends to regulate the surface temperature to a seasonally dependent value which is uniform over the globe.Finally, the likely global circulation regimes for each model atmosphere as a function of temperature were investigated and it was concluded that if CH₄, O₂, or CO dominates the atmosphere, Pluto will exhibit cyclic variations between an axially symmetric circulation system at perihelion and a baroclinic wave regime at aphelion. However, if N₂ dominates, as is likely, the wave regime should hold continuously. If the atmosphere collapses to a thin halo during summer/winter seasons, only a weak, symmetric circulation should occur.     

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

The role of ion-molecule reactions in the conversion of N2O5 to HNO3 in the stratosphere

We have investigated the role of several ion-molecule reactions in the conversion of N₂O₅ to HNO₃. In the proposed conversion, an N₂O₅ molecule would react with an H₂O molecule clustered to an inert ion to produce two HNO₃ molecules. Subsequent clustering of an H₂O molecule to the inert ion would make the reaction catalytic. If such an ion-catalysed conversion of N₂O₅ to HNO₃ occurs, it would probably play a role in the stratospheric chemistry at high latitudes in winter. In this paper we present reaction rate constant measurements made in a flowing afterglow apparatus for hydrated H₃O⁺, H⁺(CH₃CN)_m (m = 1, 2, 3), and several negative ions reacting with N₂O₅. Slow rate constants were found for these ions for hydration levels that are predominant in the stratosphere. With the known stratospheric ion density, these slow rate constants preclude significant N₂O₅ conversion by ion-molecule reactions.     

Changes in thermospheric composition inferred from twilight O+(2P) emission

Measurements of the twilight enhancement of airglow emission from O⁺(²P) near 7325 Å reveal major changes which accompany geomagnetic activity, no significant distance between evening and morning and an increase in brightness paralleling the approach to solar maximum. The principal source for O⁺(²P) is direct photoionization from O(³P) but at low solar activity there appears to be a contribution from another source in early twilight which may be local photoelectron ionization into O⁺(²P). The geomagnetic and solar effects appear to reflect changes in the O and N₂ density in the thermosphere; ground based twilight measurements of O⁺ emissions thus provide a simple means for monitoring thermospheric structure from 300 km to ～ 500 km at solar minimum and to ～600 km at solar maximum.     

The influence of neutral winds on He+ distributions in the equatorial ionosphere

The effects of F-region neutral winds on the distribution of He⁺ in the equatorial ionosphere have been examined using a theoretical model and an observational data set. It is shown by the model that components of neutral wind in the magnetic meridian up to only 50 m s^? can produce He⁺ gradients in the northern and southern sectors of a flux tube that differ by more than 80%. This is associated with interhemisphere transport velocities of He⁺ as large as 15 m s^?1 at 800 km. A substantial latitude gradient in the He⁺ distribution across the dip equator also results from the redistribution of He⁺ The changes in the He⁺ concentration at the dip equator and the latitude distribution of He⁺ in response to different neutral wind components is determined from the model and used to construct longitude distributions of He⁺ to compare with observations made at equinox. Good agreement between the calculations and observations is obtained both at the geographic and geomagnetic equators using the relationship between neutral winds, interhemispheric transport velocity and He⁺ concentration derived from the model. If these relationships can be extrapolated to accommodate the different conditions expected during solstice, we can also discuss the He⁺ distributions expected during this season.     

A study of plasmaspheric density distributions for diffusive equilibrium conditions

We have modelled the plasmaspheric density distribution for a range of solar cycle, seasonal and diurnal conditions with a magnetic flux tube dependent diffusive equilibrium model by using experimentally determined values of ionospheric parameters at 675 km as boundary conditions.Data is presented in terms of plasmaspheric H⁺ and He⁺ density contours, total flux tube content and equatorial plasma density for a range of L-values from 1.15 to 3.0. The variation of equatorial density with L-value shows good agreement with the $\text{1}\text{L}{}^4$ dependence observed experimentally.The results show that the model predicts larger solar cycle and diurnal variation in equatorial plasma density than observed using whistler techniques. However, the whistler method requires a model to deduce the equatorial density and is therefore open to interpretation.Seasonal variations are rather artifical since in this general model we have not attempted to match equatorial densities for flux tubes emanating from the winter and summer hemispheres.     

Photochemistry of the martian atmosphere: Seasonal, latitudinal, and diurnal variations

There is a significant progress in the observational data relevant to Mars photochemistry in the current decade. These data are not covered by and sometimes disagree with the published models. Therefore we consider three types of models for Mars photochemistry. A steady-state model for global-mean conditions is currently the only way to calculate the abundances of long living species (H₂, O₂, and CO). However, our model does not fit the observed CO abundance using gas-phase chemistry and reasonable values of heterogeneous loss of odd hydrogen on the water ice aerosol. The second type of the calculated models is steady-state models for local conditions. The MGS/TES data on temperature profiles, H₂O, and dust are input parameters for these models. The calculations have been made for nine seasonal points spread over the martian year and for twelve latitudes with a step of 10° for each season. The only adopted heterogeneous reaction is a weak loss of H₂O₂ on water ice with probability of 5×¹⁰⁻⁴. The results are in good agreement with the recent observations of the O₂ dayglow at 1.27 μm and the O₃ and H₂O₂ abundances. Global maps of the seasonal and latitudinal behavior of these species have been made. The third type of models is a time-dependent model for local conditions. These models show that odd hydrogen quickly converts to H₂O₂ at the nighttime and the chemistry is switched off while the association of O, the heterogeneous loss of H₂O₂, and eddy diffusion continue. This requires significant changes in the global-mean and local steady-state models discussed above, and these changes have been properly done. The calculated diurnal variations of Mars photochemistry are discussed. The martian photochemistry at low and middle latitudes is significantly different in the aphelion period at LS=10°-130° from that in the remaining part of the year.     

Long-term spectroscopic observations of Mars using IRTF/CSHELL: Mapping of O2 dayglow, CO, and search for CH4

Long-term spectroscopic observations of the O₂ dayglow at 1.27 μm result in a map of the latitudinal and seasonal behavior of the dayglow intensity for the full martian year. The O₂ dayglow is a sensitive tracer of Mars' photochemistry, and this map reflects variations of Mars' photochemistry at low and middle latitudes. It may be used to test photochemical models. Long-term observations of the CO mixing ratio have been also combined into the seasonal-latitudinal map. Seasonal and latitudinal variations of the mixing ratios of CO and the other incondensable gases (N₂, Ar, O₂, and H₂) discovered in our previous work are caused by condensation and sublimation of CO₂ to and from the polar regions. They reflect dynamics of the atmosphere and polar processes. The observed map may be used to test global circulation models of the martian atmosphere. The observed global abundances of CO are in reasonable agreement with the predicted variations with the 11-year solar cycle. Despite the perfect observing conditions, methane has not been detected using the IRTF/CSHELL with a 3σ upper limit of 14 ppb. This upper limit does not rule out the value of 10 ppb observed using the Canada-France-Hawaii Telescope and the Mars Express Planetary Fourier Spectrometer.     

Recombination of electrons with H3+ and H5+ ions

The dissociative recombination coefficients α for capture of electrons by H₃⁺ and H₅⁺ ions have been determined as a function of electron temperature T_e using a microwave afterglow-mass spectrometer apparatus. At ion and neutral temperatures T_u+ = T_n = 240 K, the coefficient α (H₃⁺) is found to vary slowly with T_e at first, decreasing from 1.6 × 10^?7 cm³/s at T_e = 240 K to 1.2 × 10^?7 cm³/s at T_e = 500 K, thereafter falling as T_e^?1 over the range 500 K ? T_e, ? 3000 K. These results, which have a ± 20% uncertainty, agree satisfactorily over the common energy range (0.03–0.36 eV) with the recombination cross sections determined in merged beam measurements by Auerbach et al. At T₊ = T_n = 128 K, the coefficient α(H₅⁺) is found to be (1.8 ± 0.3) × 10^?6 [T_e(K)/300]^?0.69 cm³/s over the range 128 K ? T_e ? 3000 K, with a more rapid decrease, as T_e^?1, between 3000 K and 5500 K. The implications of these results for modelling planetary atmospheres and interstellar clouds are briefly touched on.     

Effects of photoelectron heating and interhemisphere transport on day-time plasma temperatures at low latitudes

The thermal balance of the plasma in the day-time equatorial F region is examined. Steady-state solutions of electron and ion temperatures are obtained, assuming the ions are O⁺ and H⁺. The theoretical concentrations of O⁺ and H⁺ and the field-aligned velocity were obtained following Moffett and Hanson (1973), while theoretical photoelectron heating rates of the electron gas were taken from Swartz et al. (1975).The results demonstrate the gross features in the electron and ion temperatures as observed at the Jicamarca Observatory and in the ion temperatures observed on the OGO-6 satellite. The rapid increase in electron temperature above 500 km at the magnetic equator is due to heating by photoelectrons created at higher latitudes and travelling up along the field lines. The rapid increase in ion temperature is due to good thermal contact with the electrons rather than the neutrals. It is shown that field-aligned interhemispheric thermal plasma flows appreciably affect these temperatures, and that, with a net plasma flow from the summer hemisphere to the winter hemisphere, the temperatures are higher in the winter hemisphere. These effects are related to the character of the ion temperature minimum observed by OGO-6 near the magnetic equator.     

Frank-Condon factors for H2O+ molecular bands

Frank-Condon factors for H₂O⁺ bands have been calculated. They are used to estimate the photon scattering coefficient g_8.0 for the (8,0) band at 6158 Å.     

Quenching of O+(2D) by electrons in the thermosphere

A major loss process for the metastable species, O⁺(²D), in the thermosphere is quenching by electrons

$\text{O}{}^{\mathrm{+}}\text{(}{}^2\text{D) + e → O}{}^{\mathrm{+}}\text{(}{}^4\text{S) + e}$

.To date no laboratory measurement exists for the rate coefficient of this reaction. Thermospheric models involving this process have thus depended on a theoretically calculated value for the rate coefficient and its variation with electron temperature. Earlier studies of the O⁺(²D) ion based on the Atmosphere Explorer data gathered near solar minimum, could not quantify this process. However, Atmosphere Explorer measurements made during 1978 exhibit electron densities that are significantly enhanced over those occurring in 1974, due to the large increases that have occurred in the solar extreme ultraviolet flux. Under such conditions, for altitudes ? 280 km, the electron quenching process becomes the major loss mechanism for O⁺(²D), and the chemistry of the N⁺₂ ion, from which the O⁺(²D) density is deduced, simplifies to well determined processes. We are thus able to use the in situ satellite measurements made during 1978 to derive the electron quenching rate coefficient. The results confirm the absolute magnitude of the theoretical calculation of the rate coefficient, given by the analytical expression k(T_e) = 7.8 × 10^?8 (T_e/300)^?0.5cm³s^?1. There is an indication of a stronger temperature dependence, but the agreement is within the error of measurement.