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.     

Thermal proton flow in the plasmasphere: The morning sector

Vertical profiles of electron density obtained in the vicinity of the plasmapause using the Alouette-II topside sounder have been analyzed to assess the presence of H⁺ flow in the topside ionosphere. The observations in the midnight sector show clearly the presence of the plasmapause; i.e. there is a sharp boundary separating the poleward regions of polar wind H⁺ flow and the more gentle conditions of the plasmasphere where light ions are present in abundance. In contrast, in the sunlit morning sector upwards H⁺ flow is deduced to be present to invariant latitudes as low as 48° (L = 2·2) in the regions normally known to be well inside the plasmasphere. The upwards H⁺ flux is sufficiently large (3 × 10⁸ ions cm^?2 sec^?1) that the plasmapause cannot be seen in the latitudinal electron density contours of the topside ionosphere. The cause for this flow remains unknown but it may be a result of a diurnal refilling process.     

Dynamical behavior of thermal protons in the mid-latitude ionosphere and magnetosphere

Satellite and other observations have shown that H⁺ densities in the mid-latitude topside ionosphere are greatly reduced during magnetic storms when the plasmapause and magnetic field convection move to relatively low L-values. In the recovery phase of the magnetic storm the convection region moves to higher L-values and replenishment of H⁺ in the empty magnetospheric field tubes begins. The upwards flow of H⁺, which arises from O⁺—H charge exchange, is initially supersonic. However, as the field tubes fill with plasma, a shock front moves downwards towards the ionosphere, eventually converting the upwards flow to subsonic speeds. The duration of this supersonic recovery depends strongly on the volume of the field tube; for example calculations indicate that for L = 5 the time is approximately 22 hours. The subsonic flow continues until diffusive equilibrium is reached or a new magnetic storm begins. Calculations of the density and flux profiles expected during the subsonic phase of the recovery show that diffusive equilibrium is still not reached after an elapsed time of 10 days and correspondingly there is still a net loss of plasma from the ionosphere to the magnetosphere at that time. This slow recovery of the H⁺ density and flux patterns, following magnetic storms, indicates that the mid-latitude topside ionosphere may be in a continual dynamic state if the storms occur sufficiently often.     

The topside ionosphere above Arecibo at equinox during sunspot maximum

The coupled time-dependent O⁺ and H⁺ continuity and momentum equations and O⁺, H⁺ and electron heat balance equations are solved simultaneously within the L = 1.4 (Arecibo) magnetic flux tube between an altitude of 120 km and the equatorial plane. The results of the calculations are used in a study of the topside ionosphere above Arecibo at equinox during sunspot maximum. Magnetically quiet conditions are assumed.The results of the calculations show that the L = 1.4 magnetic flux tube becomes saturated from an arbitrary state within 2–3 days. During the day the ion content of the magnetic flux tube consists mainly of O⁺ whereas O⁺ and H⁺ are both important during the night. There is an altitude region in the topside ionosphere during the day where ion-counterstreaming occurs with H⁺ flowing downward and O⁺ flowing upward. The conditions causing this ion-counterstreaming are discussed. There is a net chemical gain of H⁺ at the higher altitudes. This H⁺ diffuses both upwards and downwards whilst O⁺ diffuses upwards from its solar e.u.v. production source which is most important at the lower altitudes. During the night the calculated O⁺ and H⁺ temperatures are very nearly equal whereas during the day there are occasions when the H⁺ temperature exceeds the O⁺ temperature by about 300 K.     

A comparison of the temperature and density structure in high and low speed thermal proton flows

The continuity, momentum and energy hydrodynamic equations for an H⁺-O⁺ topside ionosphere have been solved self-consistently for steady state conditions similar to those found outside the plasmasphere. Results are given for undisturbed and trough conditions with a range of H⁺ outflow velocities yielding subsonic and supersonic flow. In the formulation of the equations, account was taken of the velocity dependence of ion-neutral, ion-ion and ion-electron collision frequencies. In addition, parallel stress and the nonlinear acceleration term were retained in the H⁺ momentum equation. Results computed from this model show that, as a result of Joule (frictional) heating, the H⁺ temperature rises with increasing outflow velocity in the subsonic flow regime, reaching a maximum value of about 4000 K. For supersonic flow other terms in the H⁺ momentum equation become important and alter the H⁺ velocity profile such that convection becomes a heat sink in the 1000–1500 km altitude range. This, together with the reduced Joule heating resulting from the high-speed velocity dependence of the H⁺ collision frequencies, results in a decrease in the H⁺ temperature as the outflow velocity increases. However, for all outward flows the H⁺ temperature remains substantially greater than the O⁺ temperature. With identical upper boundary velocities, the H⁺ flow velocity is higher at low altitudes for trough conditions compared with non-trough conditions, but the H⁺ temperature in the trough is lower. The form of the H⁺ density profiles for supersonic flow does not in general differ greatly from those obtained with wholly subsonic flow conditions.     

Thermal response properties of the Earth's ionospheric plasma

The thermal response of the Earth's ionospheric plasma is calculated for various suddenly applied electron and ion heat sources. The time-dependent coupled electron and ion energy equations are solved by a semi-automatic computational scheme that employs Newton's method for coupled vector systems of non-linear parabolic (second order) partial differential equations in one spatial dimension. First, the electron and composite ion energy equations along a geomagnetic field line are solved with respect to a variety of ionospheric heat sources that include: thermal conduction in the daytime ionosphere; heating by electric fields acting perpendicular to the geomagnetic field line; and heating within a stable auroral red are (SAR-arc). The energy equations are then extended to resolve differential temperature profiles, first for two separate ion species (H⁺, O⁺) and then for four separate ion species (H⁺, He⁺, N⁺, O⁺) in addition to the electron temperature. The electron and individual ion temperatures are calculated for conditions within a night-time SAR-arc excited by heat flowing from the magnetosphere into the ionosphere, and also for typical midlatitude daytime ionospheric conditions. It is shown that in the lower ionosphere all ion species have the same temperature; however, in the topside ionosphere above about 400 km, ion species can display differential temperatures depending upon the balance between thermal conduction, heating by collision with electrons, cooling by collisions with the neutrals, and energy transfer by inter-ion collisions. Both the time evolution and steady-state distribution of such ion temperature differentials are discussed.The results show that below 300km both the electrons and ions respond rapidly (<30s) to variations in direct thermal forcing. Above 600 km the electrons and ions display quite different times to reach steady state, depending on the electron density: when the electron density is low the electrons reach steady state temperatures in 30 s, but typically require 700 s when the density is high; the ions, on the other hand, reach steady state in 700 s when the density is high, and 1500–2500 s when the density is low. Between 300 and 600 km, a variety of thermal structures can exist, depending upon the electron density and the type of thermal forcing; however steady state is generally reached in 200–1000 s.     

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.     

Plasma transport in the topside venus ionosphere

We have studied the extent to which certain transport processes affect ion composition and heat flow in the daytime, topside Venus ionosphere. Particular attention is given to the conditions that prevailed during the Mariner 5 measurements, at which time the topside Venus ionosphere appeared to be in a state of diffusive equilibrium. We have found that the ion composition is sensitive to the ion temperature, the ion temperature gradient, and to relative drifts between the ion species of a few $\text{m}\text{sec}$. The electron density, on the other hand, is very insensitive to these parameters. As a consequence, ionospheric models of the topside Venus ionosphere are not likely to yield definitive information about the ion composition, the thermal structure or the flow conditions, since at present only electron density profiles are available for testing model predictions. We have also found that a relative drift between the ion species of a few $\text{m}\text{sec}$ induces an ion heat flow that is equivalent to a $\text{1}\text{K}\text{km}$ temperature gradient. This induced heat flow could influence the energy balance in the topside Venus ionosphere.     

Effects of convection electric field on the thermal plasma flow between the ionosphere and the protonosphere

Dynamic behavior of the coupled ionosphere-protonosphere system in the magnetospheric convection electric field has been theoretically studied for two plasmasphere models. In the first model, it is assumed that the whole plasmasphere is in equilibrium with the underlying ionosphere in a diurnal average sense. The result for this model shows that the plasma flow between the ionosphere and the protonosphere is strongly affected by the convection electric field as a result of changes in the volume of magnetic flux tubes associated with the convective cross-L motion. Since the convection electric field is assumed to be directed from dawn to dusk, magnetic flux tubes expand on the dusk side and contract on the dawn side when rotating around the earth. The expansion of magnetic flux tubes on the dusk side causes the enhancement of the upward H⁺ flow, whereas the contraction on the dawn side causes the enhancement of the downward H⁺ flow. Consequently, the H⁺ density decreases on the dusk side and increases on the dawn side. It is also found that significant latitudinal variations in the ionospheric structures result from the L-dependency of these effects. In particular, the H⁺ density at 1000 km level becomes very low in the region of the plasmasphere bulge on the dusk side. In the second model, it is assumed that the outer portion of the plasmasphere is in the recovery state after depletions during geomagnetically disturbed periods. The result for this model shows that the upward H⁺ flux increases with latitude and consequently the H⁺ density decreases with latitude in the region of the outer plasmasphere. In summary, the present theoretical study provides a basis for comparison between the equatorial plasmapause and the trough features in the topside ionosphere.     

Diffusion coefficients for three major ions in the topside ionosphere

Published experimental data on ion composition in the topside ionosphere are examined. For certain features (the light ion trough, the high-latitude trough, the high-latitude hole and the mid-latitude total ion concentration trough) it is pointed out that the number of major ions present may be 3 or more. Transport equations derived by Schunk and co-workers are extended to include the case of three major ions in the topside ionosphere. Specific calculations are made for the O⁺, H⁺ and He⁺ ions and the behaviour of the diffusion coefficients is discussed. From a model of the high-latitude ionization hole, described by Heelis et al., representative concentration and temperature profiles are obtained. These profiles are used to demonstrate further the behaviour of the ion diffusion coefficients.     

Interhemispheric flow of thermal plasma in a closed magnetic flux tube at mid-latitudes under sunspot minimum conditions

The coupled H⁺ and O⁺ time-dependent continuity and momentum equations are solved within a region of the L = 3 magnetic flux tube lying between (and including) the F2-layers of conjugate hemispheres. The method of solution is an extended and modified version of the Murphy et al. (1976) method. The model is used to study the coupling between the F2-layers of conjugate hemispheres during magnetically quiet periods.The results of the calculations strongly indicate that the protonosphere acts as a reservoir, with variable H⁺ content, which prevents direct coupling between the F2-layers of conjugate hemispheres. However there is generally a significant interhemispheric flow of plasma. This flow is caused by conditions in the summer and winter topside ionospheres and it maintains continuity in the plasma concentration within the protonosphere. There are times when the direction of flow is from the winter hemisphere to the summer hemisphere. It is suggested that maintenance of the winter F2-layer at night is not assisted directly by the F2-layer of the conjugate summer hemisphere.It is shown that during the first few days of protonosphere replenishment after a magnetic storm there is an upflow of H⁺ in the topside ionosphere at all times in the summer hemisphere. There is also an upflow of H⁺ during the daytime in both hemispheres. A comparison with the results obtained when the interhemispheric H⁺ flux is held permanently at zero shows that both F2-layers are little affected by the interhemispheric H⁺ flux. Nevertheless both F2-layers are affected by the H⁺ tube content of the protonosphere. When the H⁺ flux at 1000 km in one hemisphere is much greater than the H⁺ flux at 1000 km in the conjugate hemisphere, there is a corresponding signature in the interhemispheric H⁺ flux.The results suggest that there is insufficient time between magnetic storms for complete replenishment of the protonosphere to occur.     

Induced proton flow in a collapsing post-sunset ionosphere and protonosphere

The effect of the onset of post-sunset conditions on thermal proton flow is examined for mid-latitudes by numerical solution of the equations of continuity, momentum and energy balance for H⁺ and O⁺. Results are calculated for a dipole magnetic field tube situated at L = 4 and acceleration terms are included in the momentum equations. Proton flow into the ionosphere results from decay of the F2-layer. Changes in temperatures and temperature gradients following sunset may not enhance the H⁺ flow. Under extreme conditions the H⁺ flow remains subsonic. It seems unlikely that an interhemispheric flux of protons can directly maintain the nighttime F2-layer.     

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.     

Thermal ion flows in the topside auroral ionosphere and the effects of low-altitude,transverse acceleration

Topside ionospheric profiles are used to study the upward field-aligned flow of thermal O⁺ at high latitudes. On the majority of the field lines outside the plasmasphere, the mean flux is approximately equal to the mean polar wind measured by spacecraft at greater altitudes. This is consistent with the theory of thermal light ion escape supported, via charge exchange, by upward O⁺ flow at lower heights. Events of larger O⁺ flow are detected at auroral latitudes and their occurrence is found to agree with that of transversely accelerated ions within the topside ionosphere and the magnetosphere. The effects of low altitude heating of O⁺ by oxygen cyclotron waves, driven by downward field-aligned currents, are considered as a possible common cause of these two types of event.     

Ion temperature troughs induced by a meridional neutral air wind in the night-time equatorial topside ionosphere

A mathematical model has been developed to calculate consistent values for the O⁺ and H⁺ concentrations and field-aligned velocities and for the O⁺, H⁺ and electron temperatures in the night-time equatorial topside ionosphere. Using the results of the model calculations a study is made to establish the ability of F-region neutral air winds to produce observed ion temperature distributions and to investigate the characteristics of ion temperature troughs as functions of altitude, latitude and ionospheric composition. Solar activity conditions that give exospheric neutral gas temperatures 600 K, 800 K and 1000 K are considered.It is shown that the O⁺-H⁺ transition height represents an altitude limit above which ion cooling due to adiabatic expansion of the plasma is extremely small. The neutral atmosphere imposes a lower altitude limit since the neutral atmosphere quenches any ion cooling which field-aligned transport tends to produce. The northern and southern edges of the ion temperature troughs are shown to be restricted to a range of dip latitudes, the limiting dip latitudes being determined by the magnetic field line geometry and by the functional form of the F-region neutral air wind velocity. Both these parameters considerably influence the interaction between the neutral air and the plasma within magnetic flux tubes.     

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.     

The transient response of the topside ionosphere to precipitation

A numerical time-dependent model of the topside and F-layer ionosphere is used to describe how the density of O⁺ ions and the plasma temperatures change as a result of transient electron precipitation with a soft energy spectrum (ca. 100 eV per electron). The response time for electron gas heating is about 2 min; for changes in topside scale height it is from 5 to 15 min, depending on altitude; and for changes in F-layer peak density, it is more than an hour. The low-density topside ion gas is thermally isolated on a short time scale; consequently the ion temperature responds almost adiabatically to volume changes. A transient precipitation event (of, say, 10 min duration) initiates a disturbance that propagates upward at approximately the sonic upeed in the plasma (ca. 2km/s), growing in amplitude with height. Such an event has little effect on the density at the peak of the F layer. An element of ionosphere that drifts horizontally in an antisunward direction through the magnetospheric cleft and into the polar cap recieves some ionization from the cleft, but not enough to be decisive in its survival. The collapse of the topside when heating is removed increases temporarily the density of the F layer.     

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.     

The atmosphere of Io from Pioneer 10 radio occultation measurements

The occultation of the Pioneer 10 spacecraft by Io (JI) provided an opportunity to obtain two S-band radio occultation measurements of its atmosphere. The dayside entry measurements revealed an ionosphere having a peak density of about 6 × 10⁴ elcm^?3 at an altitude of about 100 km. The topside scale height indicates a plasma temperature of about 406 K if it is composed of Na⁺ and 495 K if N₂⁺ is principal ion. A thinner and less dense ionosphere was observed on the exit (night side), having a peak density of 9 × 10³ elcm^?3 at an altitude of 50 km. The topside plasma temperature is 160 K for N₂^? and 131 K for Na⁺. If the ionosphere is produced by photoionization in a manner analogous to the ionospheres of the terrestrial planets, the density of neutral particles at the surface of Io is less than 10¹¹?10¹² cm³, corresponding to a surface pressure of less than 10^?8 to 10^?9 bars. Two measurements of its radius were also obtained yielding a value of 1830 km for the entry and 192 km for the exit. The discrepancy between these values may indicate an ephemeris uncertainty of about 45 km. The two measurements yield an average radius of 1875 km, which is not in agreement with the results of the Beta Scorpii stellar occultation.