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.     

A modelling study of the effects of neutral air winds on electron content at mid-latitudes in winter

A modelling study of the effects of neutral air winds on the electron content of the mid-latitude ionosphere and protonosphere in winter has been made. The theoretical models are based on solutions of time dependent momentum and continuity equations for oxygen and hydrogen ions. The computations are compared with results from slant path observations of the ATS-6 radio beacon made at Lancaster (U.K.) and Boulder, Colorado (U.S.A.).It is found that the magnitude of the poleward neutral air wind velocity has a strong effect on the general magnitude of the electron content, but that the daily pattern of electron content variation is relatively insensitive to changes in the magnitude and phase of the wind pattern. These results are in contrast with the behaviour reported previously (Sethia et al., 1983) for summer conditions. However, the night-time electron content is increased by advancing the phase of the neutral air wind and decreased by retarding it. It appears that day-to-day variations in the electron content pattern in winter cannot be explained as effects of changing neutral air winds, which again contrasts with the findings for summer. As in summer, the wind has a major effect on the filling of the protonosphere, but in opposite sense.It is argued that the effect of the neutral air wind on the ionospheric and the protonospheric electron contents depends on the duration of the poleward wind in relation to daylight and on whether or not the wind reverses direction whilst the ionosphere is sunlit.     

Ion transition heights from topside electron density profiles

Theoretical electron density profiles are calculated for the topside ionosphere to determine the major factors controlling the profile shape. Only the mean temperature, the vertical temperature gradient and the $\text{O}{}^{\mathrm{+}}\text{H}{}^{\mathrm{+}}$ ion transition height are important. Vertical proton fluxes alter the ion transition height but have no other effect on the profile shape. Diffusive equilibrium profiles including only these three effects fit observed profiles, at all latitudes, to within experimental accuracy.Values of plasma temperature, temperature gradient and ion transition height h_tT were determined by fitting theoretical models to 60,000 experimental profiles obtained from Alouette l ionograms, at latitudes of 75°S–85°N near solar minimum. Inside the plasmasphere h_T varies from about 500 km on winter nights to 850 km on summer days. Diurnal variations are caused primarily by the production and loss of O⁺ in the ionosphere. The approximately constant winter night value of h_T is close to the level for chemical equilibrium. In summer h_T is always above the equilibrium level, giving a continual production of protons which travel along lines of force to aid in maintaining the conjugate winter night ionosphere. Outside the plasmasphere h_T is 300–600 km above the equilibrium level at all times. This implies a continual near-limiting upwards flux of protons which persists down to latitudes of about 60° at night and 50° during the day.     

The effects of neutral air winds on the electron content of the mid-latitude ionosphere and protonosphere in summer

The effects of neutral air winds on the electron content (N_T) and other parameters of the mid-latitude ionosphere have been modelled by means of mathematical solutions of the time-dependent continuity and momentum equations for oxygen and hydrogen ions. The geometry is chosen to represent a propagation path between a geosynchronous satellite and a ground station, and the computations are compared with results from slant path observations of the ATS-6 radio beacon made at Lancaster (U.K.) and Boulder, Colorado (U.S.A.).It is demonstrated that the electron content responds markedly to the magnitude and phase of the neutral air winds and that the effect induced by the wind on the electron content shows a consistent quantitative relationship with the wind velocity, especially during daytime. Reasonable variations in the phase and magnitude of the wind produce a range of daily electron content patterns which encompass the range of daily variations observed.The computations show that the wind gives rise to enhanced filling of the protonosphere. This shows as a depressed value of the shape factor (F), which by definition means that a greater fraction of the ionization is at higher altitudes. The depression of F is enhanced by a poleward wind and is suppressed or even superseded by an equatorward wind through changes of the electron density distribution with altitude.     

Nighttime proton fluxes at Millstone Hill

Incoherent scatter measurements of electron density and vertical O⁺ fluxes over Millstone Hill (42.6°N, 71.5°W) previously have been used to study the exchange of plasma between the ionosphere and the magnetosphere. During the daytime there is usually an upward flux of O⁺ ions above about 450 km that can be measured readily and equated to the escaping proton flux. At night the O⁺ fluxes usually are downwards everywhere owing to the decay of the F-layer, and it becomes difficult to detect effects due an arriving proton flux. In a new study of the nighttime fluxes, appeal was made to the estimated abundance of the H⁺ ions in the upper F-region which can be extracted from the observations. From a study of the behavior on 25 days over the interval 1969–1973, we conclude that in the daytime the flux always is upwards and close to its limiting value. This situation persists throughout the night in summer at times of high sunspot activity (e.g., 1969). There is a period of downward flux prior to ionospheric sunrise on winter nights whose duration increases with decreasing sunspot number. As sunspot minimum is approached (e.g., in 1973) downward fluxes are encountered for a brief period prior to ionospheric sunrise in summer also. Thus, over most parts of sunspot cycle, it appears that the protonosphere supplies ionization to the winter night ionosphere, while being maintained from the summer hemisphere. This helps explain the smallness of the day-to-night variations reported for the electron content of magnetospheric flux tubes near L = 4 in the American sector.     

Electron content measurements of equatorial anomalous ionosphere using the geostationary satellite signals

Measurements of electron content (N_T) near the crest of the equatorial ionosphere anomaly in South America have been made and analysed to investigate N_T variations with solar hour, solar rotation and geomagnetic storms. The annual mean of diurnal ratio, defined as the ratio of the maximum to the minimum electron content of the day is found to be 5.0. Anomalous increases in night time electron content are observed with maxima around 2100 LMT and 2300 LMT during summer and equinoctial months. These increases are found to be linked with vertical motion of the F-layer. Spatial resonance in equatorial F-layer plasma appears to be the possible cause of these increases.     

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

Electron content modelling: The significance of protonospheric contents

The principal advance of the ATS-6 satellite beacon experiment was the ability to deduce continuously the electron content along the entire slant path from ground-based measurements of the signal group delay. This feature has been exploited in conjunction with the more usual Faraday rotation technique to separate the total electron content into ionospheric and protonospheric components. The physical validity of the deduced quantities is investigated using a mathematical model of the plasmasphere in which integration of the time-dependent continuity and momentum equations for oxygen and hydrogen ions along selected L shells yields the ion concentrations and field-aligned fluxes. The ion concentrations are then integrated along the propagation path to various ground stations from ATS-6 to give computed values for comparison with observations. The mathematical model is used with different sets of atmospheric parameters to investigate the significance of ionospheric and protonospheric contents as derived from beacon data.The calculated electron concentrations are able to reproduce mid-latitude equinoctial electron content observations. The shape parameters τ and F can also be simulated by day, but night-time values do not match the observations well, a greater protonospheric content being required. The calculations show that the quantity N_p, which is readily derived from ATS-6 observations, may be interpreted as the slant H⁺ content above some fixed height in the case of some stations (but not others) if the plasmasphere is reasonably full. The total slant content of H⁺ is approx. twice the value of N_p, though it appears that for the Lancaster raypath a closer relationship exists between N_p and the H⁺ tube content at L = 1.8. In general,N_p is most closely related to the tube content for an L value slightly greater than the minimum L intersected along the raypath.     

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.     

Observed heating effects of conjugate photoelectrons

A gridded spherical electrostatic analyzer aboard Injun 5 has been used to measure fluxes of thermal and hyperthermal electrons at subauroral latitudes in the midnight sector of the northern ionosphere between altitudes of 2500 and 850 km. Due to the offset between the geomagnetic and geographic poles hyperthermal fluxes, consisting of energetic photoelectrons that have escaped from the sunlit southern hemisphere are observed along orbits over the Atlantic Ocean and North America but not over Asia. The ambient electron temperatures (T_e) near 2500 km have their highest values at trough latitudes for all longitudes. At altitudes near 1000 km elevated electron temperatures in the trough were not a consistent feature of the data. Equatorward of the trough, in the longitude sector to which conjugate photoelectrons have access, T_e ～ 4000 K at 2500 km and ～ 3000 K at 1000 km. For regions with the conjugate point in darkness T_e ? 2300 K over the 1000–2500 km altitude range. The effective thermal characteristics of conjugate photoelectrons are studied as functions of altitude and latitude. The observations indicate that (1) at trough latitudes elevated electron temperatures in the topside ionosphere are mostly produced by sources other than conjugate photoelectrons, and (2) at subtrough latitudes, in the Alantic Ocean-North American longitude sector, conjugate photoelectrons contribute significantly to the heating of topside electrons. Much of the conjugate photoelectron energy is deposited at altitudes >2500 km then conducted along magnetic field lines into the 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.     

Coupling between the F-region and protonosphere: Numerical solution of the time-dependent equations

This paper describes a new method of solution of the time-dependent continuity and momentum equations for H⁺ and O⁺ in mid-latitude magnetic field tubes from the F-region to the equator. For each ion the equations are expressed as an integro-differential equation. This equation is treated as an ordinary differential equation and solved by a searching method. By means of this method, the distribution of H⁺ in the O⁺?H⁺ transition region and the protonosphere can be investigated and the influence of H⁺ fluxes on the F layer examined.As an example of application of the method a suggestion by Park (1971) about observed night-time enhancements of N_mF2 is examined. He suggested that lowering of the F layer some hours after a magnetic substorm may cause N_mF2 to increase because of increased ion influx from the protonosphere. In the present calculations the Flayer is maintained around a constant height for some time and then abruptly lowered. Under the conditions adopted the resulting increase in downward H⁺ flux is sufficient to maintain N_mF2 against the increased recombination but not to increase N_mF2 significantly. It is emphasised that these results are not conclusive.     

A fast computation of the diurnal secondary ion production in the ionosphere of Titan

We propose an analytic model that allows rapid computation of the secondary ion production due to electron impact from the primary photo-production in the ionosphere of Titan. The model parameters are given for each of the 5 major ion productions (N⁺₂, CH⁺₄, N⁺, CH⁺₃, N⁺⁺₂) as well as for the electron production.     

Structure and time variations of the Jovian ionosphere

The problem of the ionospheric formation in the Jovian upper atmosphere is examined. By adopting two plausible atmospheric models, we solve coupled time-dependent continuity equations for ions H₂⁺, H₅⁺, H⁺, H₃⁺ and H_eH⁺ simultaneously. It is shown that both radiative and three body association of H⁺ to H₂ are important for the determination of the structure of the Jovian ionosphere. The maximum electron density in the daytime is found to be about 10⁵ cm^?3. It is also shown that diurnal variation with large-amplitude can exist in the Jovian ionosphere.     

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.     

The slab thickness of the mid-latitude ionosphere

The thickness of the peak of the ionosphere depends primarily on the temperature T _n of the neutral gas, and corresponds approximately to an α-Chapman layer at a temperature of 0.87T _n. The overall slab thickness, as given by Faraday rotation measurements, is then τ =0.22 _n + 7km. Expansion of the topside ionosphere, and changes in the E-andFl-regions increase τ by about 20 km during the day in summer. Near solar minimum τ is increased by a lowering of the O ⁺/H ⁺ transition height; if the neutral temperature T _n is estimated, this height can be obtained from observed values of τ.Hourly values of slab thickness were determined over a period of 6 yr at 34°S and 42°S. Near solar maximum the night-time values were about 260 km in all seasons. The corresponding neutral temperatures agree with satellite drag values; they show a semiannual variation of 14 per cent and a seasonal change of 5 per cent. Daytime values of τ were about 230 km in winter and 320 km in summer, implying a seasonal change of 30 per cent in T _n. Temperatures increase steadily throughout the day in all seasons, with a rapid post-sunset cooling in summer. Downwards movements produce a large peak in τ at 0600 hr in winter. A large upwards flux, equal to about 40 per cent of the maximum (limiting) value, reduces τ for several hours after sunrise in winter. The slab thickness increases near solar minimum showing a reduction of the O ⁺/H ⁺ transition height to about 700 km in summer and 500 km in winter.     

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.     

Simultaneous storm-time increases of the ionospheric total electron content and the geomagnetic field in the dusk sector

Using the Faraday rotation technique with the ATS-3 satellite, it has been possible to monitor changes in the total electron content (N_T) of the mid-latitude ionosphere during the first day of 20 geomagnetic storms. Our analysis has shown that during the positive phase (ΔN_T > 0) of ionospheric storms the absolute magnitude of the increase in N_T exhibits a very pronounced maximum near sunset. The mean value of ΔN_T at 17:00 LT is more than five times the average ΔN_T value at local noon. This effect is basically independent of the storm commencement time and is usually associated with substantial local enhancements of the total geomagnetic field. The N_T enhancements are discussed in terms of a contraction and draining of the plasmasphere. A model is presented in which the dawn-dusk electric field responsible for the magnetospheric convection slows down the corotational motion of the plasmaspheric ionization in the dusk sector. This braking action causes a ‘pile up’ of the plasma and the magnetic field along the entire dusk sector.     

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.     

Ionospheric composition in SAR-arcs

The high electron temperatures existing within SAR-arcs can result in enhanced vibrational excitation of atmospheric N₂ molecules and, as a consequence, increase the rate coefficient of the reaction, O⁺ + N₂ → NO⁺ + N. This results in a change in the relative abundance of O⁺ and NO⁺+ in the SAR-arc region compared with that in the undisturbed ionosphere. Theoretical ion density profiles were computed by a triple ion analysis solving the mass, momentum and energy equations for O⁺, NO⁺ and O⁺₂ ions self-consistently. Although the electron temperature dependence of the recombination rate of NO⁺ is not well known, the results show that for a range of expected recombination rates NO⁺ still remains the dominant ion up to ca. 320 km at night within a bright SAR-arc. Studies were also made of the relative importance of a downward O⁺ flux and an upward ion drift in maintaining the F-region under SAR-arc conditions. It was found that the upward drift caused a marked increase in the NO⁺/O⁺ transition altitude as high as 460 km at night. However, for typical drift speeds up to 50 m sec^?1 the peak electron density was lower than experimental observations. The effect of a large, short-duration perpendicular electric field on the SAR-arc ion and electron density profiles was found to be small. In all cases considered the magnitude of the enhanced NO⁺ density as a result of vibrationally excited N₂ molecules was sufficient to prevent the electron density within the night-time SAR-arc from becoming vanishingly small.