A theoretical study of the ionospheric F region equatorial anomaly—II. results in the American and Asian sectors

When observed noontime values of the maximum electron density, NMAX(F2), in the ionospheric F2 region are plotted as a function of magnetic latitude, a curve is produced which has two peaks, one on either side of the dip equator at ±16° dip latitude. This paper theoretically investigates the daily variation of this latitudinal distribution in NMAX(F2) (the so-called Appleton or equatorial anomaly) and specifically attempts to account for the longitudinal differences observed between the American and Asian sectors.In Part II, models of the neutral atmosphere, production, loss and diffusion rates, neutral wind, and electric field are described and the electron densities obtained by solving the continuity equation utilizing these models are presented. In each sector, the extent to which the equatorial anomaly's daily variation is affected by changes in the geomagnetic field configuration, neutral wind, and $\text{E}\text{×}\text{B}$ drift is examined. It is found that development of the anomaly is most sensitive to the electric field model assumed, and that the observed differences at the magnetic equator between the American and Asian sectors could be accounted for by an upward $\text{E}\text{×}\text{B}$ drift which commences an hour or two earlier in the Asian sector.     

A theoretical study of the ionospheric F region equatorial anomaly—I. Theory

When observed noontime values of the maximum electron density, NMAX(F2), in the ionospheric F2 region are plotted as a function of magnetic latitude, a curve is produced which has two peaks, one on either side of the dip equator at ± 16° dip latitude. This paper theoretically investigates the daily variation of this latitudinal distribution of NMAX(F2) (the so-called Appleton or equatorial anomaly) and specifically attempts to account for the longitudinal differences observed between the American and Asian sectors.Part I outlines the theory involved in solving the time-dependent plasma continuity equation in which production, loss, and transport of ionization are taken into account, where the effects of neutral wind, ambipolar diffusion and $\text{E}\text{×}\text{B}$ drift are included in the transport term. By describing the geomagnetic field in two equivalent ways, $\text{B}\text{= ? ▽γ}$ and $\text{B}\text{= ▽α × ▽β}$, where α, β and γ are known magnetic scalar potentials, the spherical r, θ and φ space coordinates of the continuity equation are transformed to coordinates which define directions parallel and perpendicular to the magnetic field thus putting the equation in a form suitable for numerical integration.     

Plasma bubbles near the dawn terminator in the topside ionosphere

The physical properties of plasma bubbles in the topside ionosphere near the dawn terminator are investigated. It is assumed that the bubbles result from either a Rayleigh-Taylor or an $\text{E}\text{×}\text{B}$ instability on the bottom side of the F-layer. While the E-region is in darkness, the top and bottomsides of the ionosphere are electrically decoupled and the motion of bubbles can be described by non-linear, two-dimensional theory. After sunrise, electric fields within the bubbles discharge through the conducting lower ionosphere. The upward drift of the bubbles is effectively halted. To achieve a dayside state of diffusive equilibrium the bubbles slowly begin to collapse from the bottom.     

Lunar tidal oscillations in the horizontal ionospheric drift at the equator

The effect of lunar tides on the apparent ionospheric drift velocity (V) for an equatorial station Thumba (0.6°S dip) is computed by using nearly six years of data at fixed solar hours. Significant tides are observed in the E-region drifts, particularly around 12.00 hr and in the F-region drifts around 15.00 hr.A good correlation in the phases of the lunar tides in H and V is found to exist, suggesting a strong electrojet control of the horizontal ionospheric drift around these hours.     

A study of the relationship between geo-magnetic storms and ionospheric disturbances at mid-latitudes

The problem of the ionospheric disturbances associated with geomagnetic storms is examined with the goal of searching for a relationship between the time-developments of the two phenomena. Faraday rotation measurements of total electron content (N_T) are used to monitor the ionospheric F-region at a mid-latitude site, while a variety of geomagnetic parameters are examined as possible ways of following the geomagnetic variations. The ionospheric and geomagnetic data taken during 28 individual storms from 1967 to 1969 are used to search for a predictive scheme which can be tested using data from 17 storms in 1970. The specific aim is to find the geomagnetic parameter whose time-development can best forecast whether or not the ionospheric response will include an initial positive phase prior to the normally extended period of F-region depletions. Correlations between N_T and the geomagnetic indices K_p, and equatorial Dst(H) prove to be wholly inadequate. The local times of main-phase-onset (MPO) determined from the equatorial Dst(H) indices as well as from local horizontal component data, also prove to be unsatisfactory. The best correlations are obtained using local measurements of the total geomagnetic field (F). These results show that a storm commencement (SC) will produce an enhancement in n_t during the afternoon period following the SC unless there is an intervening post-midnight period with a strong depression of the geomagnetic field. Operationally this is taken to be a depression in F of at least 100γ near 03:00 LT     

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

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.     

The structure of the polar ionosphere during exceptionally quiet periods

Knowledge of the structure of the polar ionosphere during exceptionally quiet periods is basic for studying complicated ionospheric behaviors during disturbances. On the basis of data from an airborne ionosonde as well as a meridian chain of ground-basedionosondes, the circumpolar structure of the E,-and F-regions is elucidated. There are two circumpolar zones of E-region ionization with differing characteristics. The first is an auroral E,-layer and/or retarded type sporadic E-band that has previously (Whalen et al., 1971) been found to be identical with the continuous aurora. The second is a zone of non-retarded type spora die E located poleward of the former band. In general, discrete auroras are co-located with the latter. The main trough, a prominent feature of the night sector F-region, is most pronounced in the early morning. The main trough is bounded on the poleward side by a well defined ‘wall’ of F-region ionization. The night sector poleward trough wall is located approximately three degrees of latitude equatorward of the auroral oval. A ‘plateau’ of F-region ionization extends from the poleward trough wall to the auroral oval.     

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.     

Development of the auroral absorption substorm: Studies of pre-onset phase and sharp onset using an extensive riometer network

The development of an auroral absorption substorm has been studied using riometer measurements in the northern hemisphere. In the events studied, the onset is preceded by an absorption bay which begins to develop $\text{1?1}\text{1}\text{2}\text{h}$ before the onset. The bay may occur between L-values 3–19 and can cover as much as 150° of geomagnetic longitude, generally in the same longitudinal sector where the substorm breaks up and to the west of it. Whereas the substorm breaks up at or near the midnight meridian, the preceding bay may, in some geophysical conditions, appear in the afternoon sector. The preceding bay moves southward with a velocity between 60 and 600 ms^?1, intensifying during the movement. This equatorward movement is consistent with an E × B drift in a cross-magnetotail electric field of between 0.5 and 1 mV m^?1. The absorption at the onset exceeds that in the bay, and in the sector of break up the absorption shows a minimum just before the onset; to the west-of the break up the preceding bay continues its southward movement. In 14 cases studied, the sharp onset moved to the west with a velocity of 1–31 km s^?1, median 6 km s^?1. The onset was seen at higher L-values to the west than in the break-up sector. This applied also to the preceding bay. Whereas most onsets showed westward movement, in only about half of the cases studied was there movement towards the east. The injection area affected during the first minute of the onset was typically 1–2 L-value units, but as much as 30° of geomagnetic longitude. The onset later spread to cover 1–10 L-value units, and up to 130° of longitude. The contouring method used in the analysis of the data from the riometer is described in the Appendix.     

Empirical model for general circulation of the atmosphere at ionospheric levels above 100 km

This paper presents an empirical model for space-time distribution of the basic parameters of the general circulation of the atmosphere at ionospheric levels (E-and F-regions). The model is based on the results of a physico-statistical analysis of experimental data on the measurement of horizontal ionospheric drifts by close spaced receivers, carried out by the world network of stations in 1958–1970. This model allows an evaluation of the motion parameters at a given latitude, local time, season and the level of solar activity to be made. The limitations and shortcomings of the model are discussed, the results are compared with theoretical and semi-empirical schemes of the atmospheric general circulation, as well as with data of both rocket measurements of wind and drift measurements of plasma by the method of incoherent scatter of radio waves. The physics of the results obtained are stressed. The characteristics of the model are tested and defined using the materials af the coordinated program of drift measurements in the E-region from 8 stations of the northern hemisphere in 1971–1974. The characteristics of motions at higher latitudes and the longitudinal effect are discussed.     

Field-aligned and ionospheric currents

Zmuda and Armstrong (1974) showed that the field-aligned currents consist of two pairs; one is located in the morning sector and the other in the evening sector. Our analysis of magnetic records from the TRIAD satellite suggests that in each pair the poleward field-aligned current is more intense than the equatorward current, a typical ratio being 2:1. This difference has a fundamental importance in understanding the coupling between the magnetosphere and the ionosphere. We demonstrate this importance by computing the ionospheric current distribution by solving the continuity equation $\text{▽ .}\text{I}\text{= j}{}_{\mathrm{\parallel }}$ using the “observed” distribution of j_∥ for several models of the ionosphere with a high conductive annular ring (simulating the auroral oval).It is shown that the actual field-aligned and ionospheric current system is neither a simple Birkeland type, Boström type nor Zmuda-Armstrong type, but is a complicated combination of them. The relative importance among them varies considerably, depending on the conductivity distribution, the location of the peak of the field-aligned currents, etc. Further, it is found that the north-south segment of ionospheric current which connects the pair of the field-aligned currents in the morning sector does not close in the same meridian and has a large westward deflection. Thus, it has an appreciable contribution to the westward electrojet. One of the model calculations shows that the entire north-south closure current contributes to the westward electrojet.     

Empirical model of the E-region electron temperature around noon and at low magnetic activity

Electron temperature measurements made with Langmuir probes at E-region heights together with deviative absorption data show that the electrons are not in thermal equilibrium with the neutrals. Moreover, for very quiet days (A_p ? 7, K_p ? 1⁺ throughout the whole day) and hours close to noon the quotients between the electron and neutral gas temperature profiles have a similar behaviour. In this paper T_e profiles measured in situ with Langmuir probes and T_n, profiles given by neutral atmosphere models, both in the specified ionospheric conditions, are used to compute $\text{T}{}_{\mathrm{e}}\text{T}{}_{\mathrm{n}}$ profiles. Each of the profiles thus obtained is fitted by a Lorentzian curve and the variation with F_10.7 of its parameters is also fitted by simple mathematical expressions.     

Outlying plasmasphere structure detected by whistlers

Whistlers recorded at Eights ($\text{L ? 4}$) and Byrd ($\text{f ? 7}$), Antarctica have been used to study large-scale structure in equatorial plasma density at geocentric distances $\text{?3–6}\text{R}{}_{\mathrm{E}}$. The observations were made during conditions of magnetic quieting following moderate disturbance. The structures were detected by a “scanning” process involving relative motion, at about one tenth of the Earth's angular velocity or greater, between the observed density features and the observing whistler station or stations. Three case studies are described, from 26 March 1965, 11 May 1965 and 29 August 1966. The cases support satellite results by showing outlying high density regions at $\text{?4–6}\text{R}{}_{\mathrm{E}}$ that are separated from the main plasmasphere by trough-like depressions ranging in width from $\text{?0.2}$ to 1 R_E. The structures evidently endured for periods of 12 hr or more. In the cases of deepest quieting their slow east-west motions with respect to the Earth are probably of dynamo origin. The cases observed during deep quieting (11 May 1965 and 29 August 1966) suggest the approximate rotation with the Earth of structure formed during previous moderate disturbance activity in the dusk sector. The third case, from 26 March 1965, may represent a structure formed near local midnight. The reported structures appear to be closely related to the bulge phenomenon. The present work supports other experimental and theoretical evidence that the dusk sector is one of major importance in the generation of outlying density structure. It is inferred that irregularities of the type reported here regularly develop near 4–5 R_E during moderate substorm activity. This research suggests that at least a major class of the density structures that develop near 4 R_E are tail-like in nature, joined to the main body of the plasmasphere. The apparent disagreement with Chappell's results from OGO 5, which are interpreted as showing regions of “detached” plasma beyond 5 R_E, may be related to the pronounced spatial structure of electric fields observed in high-latitude ionospheric regions that are conjugate to the magnetospheric regions in which the OGO-5 observations were made.     

Interpretation of single-site scintillation measurements to estimate F-region drift velocities

A method of estimating ionospheric drift velocities using single-site scintillation measurements is applied to determine a correlation coefficient of 0.55 between magnetic activity and F-region drift velocity near the auroral ionosphere. This method is based on the relationship between the drift velocity and the scintillation spectral breakpoint.     

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.     

Electric fields in diffuse aurora

The ionospheric electric field has been measured in the E region above the Churchill auroral research range under quiet and under disturbed conditions. Results were obtained $\text{1}\text{1}\text{2}$ and $\text{2}\text{1}\text{2}$ hr before local midnight over an altitude range of 115–165 km. The instruments and analysis differ from those used by other workers. An unusually advantageous vehicle motion resulted in dipole measurements along the magnetic field being modulated by the vehicle motion. Under quiet conditions and in the presence of a diffuse, east-west 2 kR auroral arc, the predominant vector component of the electric field was also quiet and between 35 and 40 $\text{mV}\text{m}$ perpendicular to the magnetic field, southward. Parallel to the magnetic field, the vector component increased from ?17 mV/m at 130 km, reversed direction at 160 km during the latter third of the flight and fluctuated around + 6 mV/m between 155 and 135 km on the descent. Under disturbed conditions during the recovery phase of a large magnetic storm, the electric field was also more disturbed; however, there was no significant electric field along B. Analysis of effects caused when parts of the measurement system are connected by a common magnetic field line, and when one of the probes lies in the wake of the vehicle, shows that measurement perturbations produced by those effects are dominated by the magnetic field line connections and that wake effects are relatively unimportant.     

Changes in the ionospheric profile and the faraday factor M with Kp

Using incoherent scatter data from Millstone Hill, we investigated the variations in the shape of the daytime, mid-latitude ionospheric electron density profile associated with changes in geomagnetic activity. The analysis performed was to deduce the dependence upon the 3-hr geomagnetic index Kp of h(N_m), h(0·7 N_m) above and below N_m, the plasma scale height H_T in the range 500–1000 km, and the ratio $\text{N(1000)}\text{N(h}{}_{\mathrm{m}}\text{)}$. The electron density data used spanned the solar maximum years 1968–1971. Daytime data from the period 1000 to 1600 LT were averaged separately for summer, winter and spring-fall. It is shown that the mean value $\text{M}$ of the factor M = B cos θ sec χ used by Titheridge (1972) to relate the Faraday rotation $\text{Ω}$ from a geostationary satellite to the total electron content N_N up to 2000 km is practically the same (to within 1–2 per cent) as the $\text{M}$ value used to relate the N_T and $\text{Ω}$ values both computed up to 1000 km. Taking advantage of this identity, we have used the linear relationship obtained between the ionospheric parameters and Kp to deduce the height at which $\text{M}$ should be evaluated as a function of Kp.     

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.