Theoretical ion densities in the lower ionosphere

We have solved the coupled momentum and continuity equations for NO⁺, O₂⁺, and O⁺ions in the E- and F-regions of the ionosphere. This theoretical model has enabled us to examine the relative importance of various processes that affect molecular ion densities. We find that transport processes are not important during the day; the molecular ions are in chemical equilibrium at all altitudes. At night, however, both diffusion and vertical drifts induced by winds or electric fields are important in determining molecular ion densities below about 200 km. Molecular ion densities are insensitive to the O⁺ density distribution and so are little affected by decay of the nocturnal F-region or by processes, such as a protonospheric flux, that retard this decay. The O⁺ density profile, on the other hand, is insensitive to molecular ion densities, although the O⁺ diffusion equation is formally coupled to molecular ion densities by the polarization electrostatic field. Nitric oxide plays an important role in determining the NO⁺ to O₂⁺ ratio in the E-region, particularly at night. Nocturnal sources of ionization are required to maintain the E-region through the night. Vertical velocities induced by expansion and contraction of the neutral atmosphere are too small to affect ion densities at any altitude.     

The source of midlatitude metallic ions at F-region altitudes

A survey of metallic ions detected by the Bennett Ion Mass Spectrometers flown on the Atmosphere Explorer satellites, including both circular and eccentric orbital configurations, shows that patches of these ions of meteoric origin are frequently present during magnetically active periods on the bottomside of the F-layer at middle and high latitudes. In particular the F-region metals statistically tend to appear at night in the vicinity of the main ionospheric trough (in a band of invariant latitudes approx. 10 degrees wide) and on the day side of the polar cap. These distributions were previously associated with the expected dynamics of ions in the F-region above 140 km where meridional neutral wind drag and convection electric fields are the dominant ion transport mechanisms. However, the main meteor deposition layer—the presumed source region of the metals—is located below 100 km where these transport mechanisms do not prevail. It is demonstrated that the Pedersen ion drifts driven by intense electric fields such as those associated with sub-auroral ion drifts (SAID) are sufficient to transport the long-lived metallic ions upward from the main meteor layer to altitudes where the drag of equatorial directed neutral winds and electric field convection can support them against the downward pull of gravity and transport them to other locations. The spatial and temporal distribution of the middle and high latitude F-region metals are consistent with the known characteristics of the electric fields and with the expected F-region ion dynamics.     

Atmospheric expansion from Joule heating

An expression for the vertical velocity of the neutral atmosphere in the F-region is derived for Joule heating by the electric field that drives the auroral electrojet. When only vertical expansion is allowed, it is found that the vertical wind must always increase monotonically with altitude. The heating rate is proportional to the F-region ion density, so that appreciable heating, even during high electric fields, requires some production mechanism of ionization such as auroral secondary ionization or solar photoionization, in the lower F-region. Once started at night, when an ionizing source is present in the lower F-region, the expansion of the atmosphere transports ionization upward, thereby increasing the heating rate, and hence the expansion rate, i.e. positive feedback. Electric field strengths and F-region ion densities of 50 mV/m and 2 × 10¹¹e/m³, respectively, will produce vertal neutral wind speeds of several tens of m/sec in the 300–500 km altitude range. During periods of high magnetic activity, i.e. high electric field, Joule heating can produce large increases in the relative N₂ concentration in the upper F-region; computations made with a simple model suggest that tenfold increases can occur at 400 km altitude $\text{1}\text{2?1}\text{hr}$ after the onset of magnetic activity, a result in agreement with satellite observations. When the Joule heating theory is applied to incoherent scatter data taken during one period of high heating, the horizontal electric field in the F-region is found to decrease markedly, possibly approaching zero as the field penetrates a weak, discrete auroral arc; the decrease began 10–20 km from the arc.     

Time-dependent behaviour of a stable auroral red arc excited by an electric field

We examine the electric field hypothesis as a possible explanation of a stable auroral red arc. An electric field perpendicular to the geomagnetic field in the ionosphere heats the ambient F-region electrons and ions. Given large enough electric fields, the electrons can be heated sufficiently to excite the OI (¹D) term of atomic oxygen by electron impact, giving rise to the λ6300 emission characteristic of the red arc. The electron and ion heating rates are determined by the relative drift between the plasma and neutral gas.     

Determination of electric fields and neutral wind velocities in the ionospheric E- and F-regions from observations of artificial ion clouds

A new method of interpreting the behaviour of artificial ion clouds released in the Earth's ionosphere is presented. It is shown that values for the ionospheric electric field, neutral wind velocities and, in some circumstances, ion collision frequency, can be deduced from a study of the motion and deformation of the ion clouds, including those released in the E-region.     

The calculated and observed diurnal variation of the ionosphere over Millstone Hill on 23–24 March 1970

A numerical model of current F-region theory is use to calculate the diurnal variation of the mid-latitude ionospheric F-region over Millstone Hill on 23–24 March 1970, during quiet geomagnetic conditions. From the solar EUV flux, the model calculates at each altitude and time step primary photoelectron spectra and ionization rates of various ion species. The photoelectron transport equation is solved for the secondary ionization rates, photoelectron spectra, and various airglow excitation rates. Five ion continuity equations that include the effects of transport by diffusion, magnetospheric-ionospheric plasma transport, electric fields, and neutral winds are solved for the ion composition and electron density. The electron and ion temperatures are also calculated using the heating rates determined from chemical reactions, photoelectron collisions, and magnetospheric-ionospheric energy transport. The calculations are performed for a diurnal cycle considering a stationary field tube co-rotating with the Earth; only the vertical plasma drift caused by electric fields perpendicular to the geomagnetic field line is allowed but not the horizontal drift. The boundary conditions used in the model are determined from the incoherent scatter radar measurements of T_e, T_i and O⁺ flux at 800km over Millstone Hill (Evans, 1971a). The component of the neutral thermospheric winds along the geomagnetic field has an important influence on the overall ionospheric structure. It is determined from a separate dynamic model of the neutral thermosphere, using incoherent scatter radar measurements.The calculated diurnal variation of the ionospheric structure agrees well with the values measured by the incoherent scatter radar when certain restrictions are placed on the solar EUV flux and model neutral atmospheric compositions. Namely, the solar EUV fluxes of Hinteregger (1970) are doubled and an atomic oxygen concentration of at least $\text{10}{}^{11}\text{cm}{}^3$ at 120 km is required for the neutral model atmosphere. Calculations also show that the topside thermal structure of the ionosphere is primarily maintained by a flow of heat from the magnetosphere and the night-time F2-region is maintained in part by neutral winds, diffusion, electric fields, and plasma flow from the magnetosphere. The problem of maintaining the calculated night-time ionosphere at the observed values is also discussed.     

Electrical coupling of the E- and F-regions and its effect on F-region drifts and winds

Electric currents, generated by thermospheric winds, flow along the geomagnetic field lines linking the E-and F-regions. Their effects on the electric field distribution are investigated by solving the electrical and dynamical equations. The input data include appropriate models of the F-region tidal winds, the thermospheric pressure distribution and the E-and F-layer concentrations. At the magnetic equator, the calculated neutral air wind at 240 km height has a prevailling eastward component of 55 m sec^-1 and the west-east and vertical ion drifts agree in their general form with incoherent scatter data from Jicamarca     

Magnetic disturbance effects in the E region of the ionosphere

Detailed studies of the daytime E-region critical frequency at Aberystwyth (geomagnetic latitude +56°) show clear evidence for changes associated with both the axially-symmetric (Dst) and asymmetric (DS) components of the disturbance magnetic field. Comparison of the sensitivity of the E-region peak density to these two influences shows that the changes cannot entirely (if at all) be ascribed to the influence of electric currents in the region. It is suggested that a major role is played by dynamical influences associated with the neutral air “storm circulation” which distributes the energy fed into the auroral region to lower latitudes.     

Auroral E-region: Ion composition and nitric oxide

Data from eight auroral ion composition measurements, seven of which have been reported in the literature, are analyzed and compared in terms of a single model format. We find, contrary to conclusions published previously for two of the experiments, that there is no discrepancy concerning O⁺ ions. In general, the mean CIRA 1972 neutral model is found to be quite suitable as a representative of the major gas composition required for auroral E-region calculations which agree with the data. Nitric oxide profiles inferred from analysis of the data range from about normal non-auroral E-region nitrix oxide distributions with peak concentrations near 10⁸ cm^?3 to profiles with peak populations near 10⁹ cm^?3. Although the higher concentrations are generally correlated with intense aurora, we acknowledge that the length and strength of auroral activity prior to the individual rocket flights can have an even greater bearing, at times, on the NO “snapshot” profile deduced from the auroral ion composition data.     

The behaviour of the midlatitude F-region at the time of four magnetospheric substorms during the night of 29/30 October 1973

The data from observations of the geomagnetic field, ionospheric parameters and atmospheric emissions, carried out at four midlatitude station in Bulgaria are analysed. The observations refer to the geomagnetic disturbance on 28/30 October 1973 (K_p^max = 7) and also to a very quiet period before it. It is shown that all four geomagnetic substorms during the night of 29/30 October influenced the midlatitude F-region. This is indicated by a lowering of the height of the F-region by ca. 50–70 km. Owing to this downward drift of ionisation the dissociative recombination and the intensity of the red line is accordingly increased. As an explanation of this phenomenon we suggest the action of the electric fields, which can at the same time be transported from the magnetosphere to the ionosphere.     

The influence of convection electric fields on thermal proton outflow from the ionosphere

The continuity, momentum and energy hydrodynamic equations for an O⁺-H⁺ ionosphere have been solved self-consistently for steady state conditions when a perpendicular (convection) electric field is present. Comparison of the H⁺ temperature profiles obtained with and without the electric field show that the effect of the electric field is to enhance the H⁺ temperature at high altitudes from about 3600 to 6400 K. Due to ion heating by the electric field, there is a net reduction of O⁺ in the F2-region as compared with the case of a non-convecting ionosphere. When the reduction of O⁺ is neglected, the electric field acts to increase the H⁺ outward flux from 8.3 × 10⁷ to 2.7 × 10⁸ cm^?2 sec^?1 for average ionospheric conditions. However, when the reduction of O⁺ is included, there is a net reduction in the outward H⁺ flux. Nevertheless, the convection electric field still results in an increase in the rate of depletion of the F-re m^?1 electric field.     

Ionosphere-plasmasheet field-aligned currents and parallel electric fields

Two kinetic models for the auroral topside ionosphere are compared. The collisionless plasma distributed along an auroral magnetic field line behaves like a non-Ohmic conducting medium with highly non-linear characteristic curves relating the parallel current density to the potential difference between the cold ionosphere and the hot plasmasheet region. The (zero-electric current) potential difference, required to balance the current carried by the precipitating plasmasheet particles and the current transported by the outflowing ionospheric particles, depends on the ratio n_ps.e/n_th.e and T_ps.e/T_th.e of the plasmasheet and ionospheric electron densities and temperatures. When in the E-region the magnetic field lines are interconnected by a high conductivity plasma the resulting field-aligned currents driven by the magnetospheric potential distribution are limited by the integrated Pedersen conductivity of the ionospheric layers. These currents are not related to the parallel electric field intensity as they would be in Ohmic materials. The parallel electric field intensity is necessarily determined by the local quasi-neutrality of the plasma.     

The role of neutral winds and ionospheric electric field in forming stable sporadic E-layers

Using the combined measurements from a rocket flight through a stable intense sporadic E-layer, we examine the shortcomings of conventional wind shear theory of ion layer formation, principally in underestimating the role of the ambient ionospheric electric field. Our results imply that the ionospheric electric field may control the stability and precise location of such ionisation layers within a region of convergent ion flow.     

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.     

The solar M-region problem—An old problem now facing its solution?

The solar M-region problem is briefly reviewed. The Mustel and the Allen-Saemundsson M-region schools are discussed in the light of (a) recent results on coronal structures and solar wind variations, and (b) statistical analyses of coronal-geomagnetic correlations. From this discussion it is suggested that the M-regions should be identified with the central portion of magnetically open solar regions, or coronal holes. With reference to the papers by Gulbrandsen (1973b, 1974), it is concluded that such an identification is in fact favoured by the great majority of analyses published on the M-region problem during the last 40 yr.It is assumed that historically, the M-region conflict has evolved from the scientific pecularity that geomagnetic phenomena which are actually related to magnetically open structures, were generally not studied in relation to such structures, but to closed magnetic configurations.     

Modeling the midlatitude F-region ionospheric storm using east-west drift and a meridional wind

Under magnetically quiet conditions, ionospheric plasma in the midlatitude F-region corotates with the Earth and relative east-west drifts are small compared to the corotation velocity. During magnetic storms, however, the enhanced dawn-to-dusk magnetospheric convection electric field often penetrates into the midlatitude region, where it maps into the ionosphere as a poleward electric field in the 18:00 LT sector, producing a strong westward plasma drift. To evaluate the ionospheric response to this east-west drift, the time-dependent O⁺ continuity equation is solved numerically, including the effects of production by photoionization, loss by charge exchange and transport by diffusion, neutral wind and $\text{E}\text{×}\text{B}$ drift. In this investigation only the neutral wind's meridional component and east-west $\text{E}\text{×}\text{B}$ drift are included. It is found that an enhanced equatorward wind coupled with westward drift produces an enhancement in the peak electron density (NMAX(F2)) and in the electron content (up to 1000 km) in the afternoon sector and a subsequent greater-than-normal decay in ionization after 18:00 LT. These results agree in general with midlatitude F-region ionospheric storm observations of NMAX(F2) and electron content which show an afternoon enhancement over quiet-time values followed by an abrupt transition to lower-than-normal values. Westward drift appears to be a sufficient mechanism in bringing about this sharp transition.     

Ion momentum and energy transfer rates for charge exchange collisions

The rates of momentum and energy transfer have been obtained for charge exchange collisions between ion and neutral gases having arbitrary Maxwellian temperatures T_i and T_n and bulk transport velocities c_i and c_n. The results are directly applicable to the F-region of the ionosphere where O⁺ - O charge is the dominant mechanism affecting ion momentum and energy transfer.     

Photoionization rates in the night-time E- and F-region ionosphere1

Quantitative estimates of ionization sources that maintain the night-time E- and F-region ionosphere are given. Starlight (stellar continuum radiation in the spectral inverval 911–1026 Å) and resonance scattering of solar Ly-β into the night sector are the most important sources in the E-region and are capable of maintaining observable electron densities of order (1–4) × 10³ cm^?3. Starlight ionization rates have substantial variations (factors of 2–4) with latitude and time of year since the brightest stars in the night sky occur in the southern Milky Way and Orion regions. In the lower F-region the major O⁺ source in the equatorial ionosphere is 910 Å radiation from the O⁺ recombination in the F2-region, whereas in the extratropical ionosphere interplanetary 584 Å radiation only exceeds resonance scattering of solar 584 and 304 Å radiation as the dominant O⁺ source during the month of December.     

Antarctic substorm events observed by sounding rockets,ionization of the D- and E-regions by auroral electrons

The ionization structure of the auroral arc was measured on a sounding rocket which penetrated into a bright auroral arc. The E-region electron density becomes large (2 ～ 5 × 10⁵ el/cm³ only in the moving auroral arc, whose N₂⁺ 4278 Å brightness is 1 ～ 2·5 kR. The electron density in the D-region beneath the lower boundary of the arc (75 ～ 98 km in altitude) is also considerably enhanced to 2 ～ 5 × 10⁴ el/cm³.The observed E-region electron density can be interpreted theoretically as due to the direct ionization by precipitating electrons, whose energy spectrum is approximately represented by an exponential type having the characteristic energy of 2 keV. The correlation between the electron density and the N₂⁺ 4278 Å brightness can be reasonably explained by considering the simultaneous effects on the ionization and the optical excitation caused by the primary electrons having a flux of 9 × 10⁹ el/cm²/sec per 1 kR of the 4278 Å emission.Further analyses using the electron density data from four other sounding rockets have shown that the D-region ionization has good correlations to the cosmic noise absorption (CNA) and the magnetic substorm activities observed simultaneously at the ground station, whereas it has poor correlation to the same quantity of the E-region measured in the same experiment. It is found that the observed D-region ionization is much larger than that predicted by the theory which takes into account the Bremsstrahlung X-ray ionization along with the direct impact ionization when it is applied to the precipitating electron flux spectrum consistent to the E-region ionization and optical excitation.After all the present experimental results suggest a dual nature of the electron precipitation spectrum in the substorm, i.e. the softer part which is localized in the auroral arc and the harder part which is spatially wide-spread over the substorm area.     

Measurement of nitric oxide and related parameters in the equatorial mesosphere and lower thermosphere

A sounding rocket was launched in March 1982 from Thumba, India, shortly after sunrise. The measurements included the concentration of nitric oxide and ozone, the total ion density and the Lyman-α flux. Hence most parameters important for the formation of the D-region during daytime are available with the exception of solar radiation other than Lyman-α which only becomes important above 95 km. The mutual agreement is satisfactory which adds weight to the measurements.