Satellite observations of zonal electric fields near sunrise in the equatorial ionosphere

We report here on a number of examples of anomalous enhancements of eastward electric fields near sunrise in the equatorial ionospheric F-region. These examples were selected from the data base of the equatorial satellite, San Marco D (1988), which measured ionospheric electric fields during a period of solar minimum. The eastward electric fields reported correspond to vertical plasma drifts. The examples studied here are similar in signature and polarity to the pre-reversal electric field enhancements seen near sunset from ground-based radar systems. The morphology of these sunrise events, which are observed on about 14% of the morning-side satellite passes, are studied as a function of local zonal velocity, magnetic activity, geographic longitude and altitude. The nine events studied occur at locations where the zonal plasma flow is generally measured to be eastward, but reducing as a function of local time and at satellite longitudes where the magnetic declination has the opposite polarity as the declination of the sunrise terminator.     

The equatorial electrojet and associated currents as seen in Magsat data

Magsat data are re-examined with regard to the presence and character of fields due to the equatorial electrojet and meridional currents at dawn and dusk local times. Dip-latitude organized field variations at dawn are:

• 1.(1) extremely weak,
• 2.(2) extremely variable with longitude,
• 3.(3) inconsistent with the pattern expected from a line or narrow sheet current.

It is shown that the use of Magsat dusk data can ‘contaminate’ a main field model, introducing apparent equatorial electrojet effects into the dawn data.Fields due to the equatorial electrojet and (presumably) associated meridional currents are clearly present in the dusk data. They show a variation with longitude which is apparently associated with the longitudinal variation of the strength, or square of the strength, of the main field in the E-region. Also evident is a variation with time of the year, although data are available for only a six month period. The meridional currents are generally minimum during January and February and maximum either during November and December or March and April, depending upon longitude. The E-region horizontal currents are minimum in November and December and maximum in March and April, except for − 30° to −90° longitude when the maximum occurs in January and February.Assuming that field gradients in local time are considerably smaller than field gradients in dip-latitude, current densities are estimated to be 1–3.6μA/m² for the horizontal current at 110km and about 10–20 × 10⁻⁹ A/m² for the vertical currents at 400km altitude. These results confirm and extend earlier results of Takeda and Maeda.Most models of the electrojet system in the literature disagree severely with these measurements either because their scope is inadequate or because of the wind system they assume. Those models which best describe the data invoke an eastward wind and/or an eastward electric field at dusk local time.     

Worldwide atmospheric gravity-wave study in the European sector 1985–1990

Four campaigns of the Worldwide Atmospheric Gravity-wave Study (WAGS) have taken place in the European sector. On many occasions the onset of auroral activity in the evening and midnight sector, as indicated by EISCAT measurements of the electric field, was associated after a suitable delay with the detection of periodic ionospheric disturbance travelling southward over the U.K. at speeds between 500 and 1000 m s⁻¹. The velocity and wavelength of the TIDs corresponded to large-scale atmospheric gravity-waves. The characteristic periods of the travelling disturbances were similar to the intrinsic time scales of the auroral activity for periods of 40 min or more, but variations on a time scale of 20 min or less were strongly attenuated. The r.m.s. amplitude of the auroral electric field was proportional to the r.m.s. amplitude of the HF Doppler-shift associated with the gravity-wave. The time-lag between the onset of strong auroral activity and the arrival of the travelling disturbance over the U.K. was usually about an hour, suggesting a source region about 2000 km north. Similar levels of activity in the afternoon did not appear to produce strong waves in the far field. This is possibly due to ion-drag in the daytime ionosphere although the effects of the lower sensitivity of the HF Doppler-network during daytime must also be considered.     

Electric field generation at the magnetospheric boundary for northward IMF

We discuss three different processes which generate electric fields at the magnetopause during northward interplanetary magnetic field (IMF) conditions. These are (1) Petschek-type magnetic field reconnection, (2) magnetic field diffusion, and (3) viscous-like interaction resulting from the Kelvin-Helmholtz instability. For northward IMF all three processes lead to the formation of a boundary layer on closed magnetic field lines adjacent to the magnetospheric boundary. The thickness of the boundary layer depend on Petschek's parameter in the first case, the magnetic Reynolds number in the second case, and an effective Reynolds number in the third case. In each case coupling between the boundary layer and the ionosphere occurs via field-aligned currents. These field-aligned currents result from the penetration into the polar ionosphere of the electric field generated at the magnetospheric boundary. These currents are closed by a transverse current in the boundary layer and the associated Lorentz force causes a decrease of the kinetic energy of the solar wind plasma inside the boundary layer. As a result of this velocity decrease the thickness of the boundary layer increases on both flanks of the magnetosphere near the equatorial plane. The convergence of the boundary layer on the dawn and dusk sides leads to antisunward plasma flow in the magnetospheric tail.     

Modeling equatorial ionospheric electric fields

This paper gives a brief overview of the processes responsible for the equatorial electric field, and reviews relevant modeling work of these processes, with emphases on basic aspects and recent progress. Modeling studies have been able to explain most of the observed features of equatorial electric fields, although some uncertainties remain. The strong anisotropy of the conductivity and the presence of an east-west electric field lead to a strong vertical polarization electric field in the lower ionosphere at the magnetic equator, whose magnitude can be limited by plasma irregularities. Local winds influence the structure of the equatorial polarization field in both the E and F regions. The evening pre-reversal enhancement of the eastward electric field has been modeled by considering a combination of effects due to the presence of a strong eastward wind in the F region and to east-west gradients of the conductivity, current, and wind. Models of coupled thermosphere-ionosphere dynamics and electrodynamics have demonstrated the importance of mutual-coupling effects. The low-latitude east-west electric field arises mainly from the global ionospheric wind dynamo and from the magnetospheric dynamo, but models of these dynamos and of their coupling have not yet attained accurate predictive capability.     

The morphology of dayside Birkeland currents

Intense (10⁵ A) electric currents flow into and from the Earth's two polar ionospheres near magnetic noon. These currents, called Birkeland or magnetic field-aligned currents, are the agent by which momentum couples from the flowing solar wind plasma to drive plasma motions in the high latitude ionosphere. Coupling is strongest when the interplanetary magnetic field (IMF) has a southward component and when this occurs there exist two principal regions of Birkeland current near magnetic noon called the region 1 and the cusp systems. We present a simple model bringing theoretical order to the many patterns proposed previously for the morphology of these dayside Birkeland currents as observed by orbiting satellites in the topside polar ionosphere. Specifically we show that the cusp Birkeland current system is not a latitudinally separate region but is instead the extension in longitude of the region 1 Birkeland current from either dawn or dusk; which particular one depends on the sign of the east-west (Y) component of the IMF. The presence of an IMF Y-component therefore leads to two region 1 current systems near magnetic noon, with the poleward one being that previously called the ‘cusp’ system.     

Vertical neutral wind in the equatorial F-region deduced from electric field and ion density measurements

DC electric field and ion density measurements near density depletion regions (that is, equatorial plasma bubbles) are used to estimate the vertical neutral wind speed. The measured zonal electric field in a series of density depletions crossed by the San Marco D satellite at 01.47-01.52 UT on 25 October 1988, can be explained if a downward neutral wind of 15–30 m s⁻¹ exists. Simultaneously, the F-region plasma was moving downward at a speed of 30–50 m s⁻¹ These events appear in the local time sector of 23.002̄23.15 in which strong downward neutral winds may occur. Indeed, airglow measurements suggest that downward neutral velocities of 25–50 m s⁻¹ are possible at times near midnight in the equatorial F-region.     

Plasmaspheric zonal electric fields and coupling fluxes from the Dunedin VLF Doppler experiment,for 180°E,L ≈ 2.3, at solstice and equinox

Group delays and Doppler shifts from ducted whistler-mode signals are measured using the VLF Doppler experiment at Dunedin, New Zealand (45.8°S, 170.5°E). Equatorial zonal electric field and plasmasphere-ionosphere coupling fluxes are determined for L ≈ 2.3 at June solstice and equinox during magnetically quiet periods. The general features of the electric field measured at Dunedin agree with those predicted from ionospheric dynamo theory with a (1,−2) tidal component. Some seasonal variations are observed, with the electric field measured during equinox being smaller and predominantly westward during the night. The electric field at June solstice is also westward during the evening and for part of the night, but turns sharply eastward during the pre-dawn and dawn period at the duct entry site. The June electric field appears to follow a diurnal variation whereas the equinox electric field shows a possible 4-hourly periodic variation. Seasonal variations in the neutral wind pattern, altering the configuration of the ionospheric dynamo field, are the probable cause of the seasonal differences in the electric field. The seasonal variation of the coupling fluxes can be explained by the alteration of the E x B drift pattern, caused by the changes in the electric field.     

Low latitude electrodynamic plasma drifts : a review

Radar and radio measurements have provided detailed information on the dependence of F-region electrodynamic drifts on height, season, solar cycle and magnetic activity. Recently, satellite ion drift and electric field probes have determined the variation of low latitude ionospheric drifts over a large range of altitudes and latitudes. The general characteristics of the quiet time plasma can be explained as resulting from E- and F-region dynamo and interhemispheric coupling processes. The low latitude and equatorial zonal and upward/poleward components of the plasma drift respond differently to geomagnetic activity. Disturbance dynamo effects are responsible for the drift perturbations following periods of enhanced magnetic activity. The prompt penetration of high latitude electric fields to lower latitudes produces large perturbations on the upward/poleward drifts, but has no significant effect on the low latitude and equatorial zonal drifts. A number of processes such as ‘overshielding’, ‘fossil wind’ and magnetic reconfiguration were suggested as being responsible for the direct penetration of high latitude electric fields to lower latitudes. Detailed low latitude and global numerical models were used to study the characteristics of low latitude and equatorial plasma drifts and their response to changes in the polar cap potential drop or in the high latitude field-aligned currents. These models can reproduce the latitudinal variation of the perturbation electric fields and their diurnal variations, but are still unable to account for several aspects of the experimental data as a result of the complexity of the high latitude and magnetospheric processes involved.     

Meridional electric currents in the equatorial F-region

A study of the boundary conditions for the equatorial thermospheric transport equations by the authors has led to the theoretical prediction of the vertical electric field at the base of the F-region. Earlier, this result was applied to the calculation of the zonal wind field in the equatorial F-region. In this work, the aforementioned model is applied to the calculation of the F-region electric current field in the meridional plane as a function of time and the east-west magnetic field generated by these currents. In particular, the field at sunset is compared with the observations made by Magsat.     

Electric field sources in the quiet plasmasphere from whistler observations

Existing evidence for the ionospheric dynamo being the source of quiet time electric fields in the plasmasphere is reviewed. Part of a 24 h set of whistler data recorded continuously at Sanae, Antarctica (L = 4), during quiet magnetic (average K_p = 1) is analysed to obtain westward electric fields in the equatorial plane. These electric fields are examined as a function of L-value in order to infer their source. It is found that for periods of outward flow of plasma during the noon-midnight local time period, the electric fields are consistent with the dominant source being the ionospheric dynamo. There is some evidence that during the evening period of inward flow the electric fields are magnetospheric in origin, although this could also be consistent with a refined dynamo model. The observed whistler duct convection patterns do not fit either of two theoretical models, which invoke a magnetospheric field but not a dynamo field.     

What determines the latitudinal extent of the auroral acceleration region?

Characteristic scales associated with auroral precipitation are investigated on the basis of quasistatic magnetotail models, resistive MHD simulations of magnetotail dynamics, and a general relation between parallel electric fields and velocity shear. Since the inverted-V precipitation region of discrete auroras (on the dusk side) is associated with upward flowing, region 1, currents, we investigate the distribution of these currents first. The overall distribution of region 1 type field-aligned currents and their dynamic changes can be explained by characteristic scales in the magnetotail and their mapping to the ionosphere. The quiet time region 1 currents are associated with the decrease of tail flaring. Their overall extent in the north-south direction is closely related to the scale height of the cross-tail current. Dynamic region 1 currents are related to the velocity shear of earthward flow, which can be generated by a tail instability. In that case the peaks of the enhanced region 1 currents are found to map closer to midnight and to lower latitudes than the quiet region 1 currents, consistent with average observations [Iijima and Potemra (1976a), J. geophys. Res.81, 2165]. On the basis of a general relation between parallel electric fields and ‘slippage’ in the plasma transport [Hesse and Schindler (1988), J. geophys. Res.93, 5559; Schindleret al. (1991), Astrophys. J.380, 293], we make estimates of the spatial extent of nonideal regions, where parallel electric fields may exist. For a plausible model of substorm reconfiguration, we find a latitudinal extent of about 7 km for a time scale of 1 min and a integrated parallel electric field of 5 kV. The length scale is proportional to the time scale. The sign of this parallel electric field is consistent with downward acceleration of electrons on the dusk side. The spatial extent of the parallel electric field region is independent of the microscopic generation mechanism if the time scale and the characteristic parallel potential difference (i.e. the integrated parallel electric field) are independent of this mechanism.     

Zonal irregularity drifts and neutral winds measured near the magnetic equator in Peru

Measurements of zonal irregularity drifts were made by the spaced receiver scintillation and radar interferometer techniques from Huancayo and Jicamarca, respectively. The Fabry-Perot Interferometer operated at Arequipa provided the zonal neutral winds. These simultaneous measurements were performed during evening hours in the presence of equatorial spread-F on three nights in October 1988. The zonal drift of 3-m irregularities obtained with the 50-MHz radar showed considerable variation as a function of altitude. The drift of hundreds of m-scale irregularities obtained by the scintillation technique agreed with the drift of 3-m irregularities when the latter were measured near the F-peak. The neutral winds, on the other hand, sometimes exceeded the irregularity drifts by a factor of two. This is a possible result of the partial reduction of the vertical polarization electric field in the F-region caused by the effects of integrated Pedersen conductivity of the off-equatorial night-time E-region coupled to the F-region at high altitudes above the magnetic equator.     

Simulations of the behaviour of the polar thermosphere and ionosphere for northward IMF

The dynamics and structure of the polar thermosphere and ionosphere within the polar regions are strongly influenced by the magnetospheric electric field. The convection of ionospheric plasma imposed by this electric field generates a large-scale thermospheric circulation which tends to follow the pattern of the ionospheric circulation itself. The magnetospheric electric field pattern is strongly influenced by the magnitude and direction of the interplanetary magnetic field (IMF), and by the dynamic pressure of the solar wind. Previous numerical simulations of the thermospheric response to magnetospheric activity have used available models of auroral precipitation and magnetospheric electric fields appropriate for a southward-directed IMF. In this study, the UCL/Sheffield coupled thermosphere/ionosphere model has been used, including convection electric field models for a northward IMF configuration. During periods of persistent strong northward IMF B_z, regions of sunward thermospheric winds (up to 200 m s⁻¹) may occur deep within the polar cap, reversing the generally anti-sunward polar cap winds driven by low-latitude solar EUV heating and enhanced by geomagnetic forcing under all conditions of southward IMF B_z. The development of sunward polar cap winds requires persistent northward IMF and enhanced solar wind dynamic pressure for at least 2–4 h, and the magnitude of the northward IMF component should exceed approximately 5 nT. Sunward winds will occur preferentially on the dawn (dusk) side of the polar cap for IMF B_y negative (positive) in the northern hemisphere (reverse in the southern hemisphere). The magnitude of sunward polar cap winds will be significantly modulated by UT and season, reflecting E-and F-region plasma densities. For example, in northern mid-winter, sunward polar cap winds will tend to be a factor of two stronger around 1800 UT, when the geomagnetic polar cusp is sunlit, then at 0600 UT, when the entire polar cap is in darkness.     

On the structure and circulation of the polar thermosphere under magnetically quiet conditions

A theoretical study is presented bearing on the thermospheric circulation and composition at polar latitudes. The observed motions and density perturbations in N₂, O and He have signatures which may be understood in terms of two different source mechanisms. We consider electric field momentum coupling and Joule heating as well as interactions between both processes. A spectral model in terms of vector spherical harmonics (with magnetic coordinates) is used, delineating the diurnal and mean (time independent) components. The important non-linearities are evaluated in configuration space. The electric field model of Volland and the global average density and temperature variations of Hedin (MSIS) are adopted as input. Our analysis leads to the following conclusions. (1) The vortex type double cell polar circulation (zonal wave number m = 1) is primarily driven by collisional momentum transfer from electric field induced ion convection. (2) Because of the thermospheric low pass filter, a large time independent component (zonal wave number m = 0) is produced by Joule heating; the heavier species (N₂) being concentrated where the lighter ones (O, He) are depleted, and vice versa due to wind induced diffusion (3) The electric field driven vortex circulation redistributes the mass and energy in the time independent density and temperature variations (from Joule heating), producing primarily diurnal variations (m = 1) in the temperature and composition near the pole and at auroral latitudes, again the heavier and lighter species varying out of phase. The above results are in substantial agreement with observations. It is worth noting that momentum rectification associated with the diurnally varying electric field and conductivity induces a weak zonally symmetric (single cell) prograde polar vortex. However, this motion is partially compensated by a retrograde vortex from geostrophic balance due to Joule heating, which dominates near the pole. These motions are small compared with the diurnally varying component in the polar circulation.     

Equatorial penetration of magnetic disturbance effects in the thermosphere and ionosphere

The neutral dynamic and electrodynamic coupling between high and low latitudes, and the mutual interactions between these two processes, are investigated. For 22 March 1979, when a sudden increase in magnetic activity occurred, we have analyzed the following experimental data: (a) neutral densities and cross-track neutral winds as a function of latitude (0°–80°) near 200 km from a satellite-borne accelerometer; (b) hourly mean H-component magnetic data from the Huancayo Observatory (0.72°S, 4.78°E; dipole geomagnetic coordinates) magnetometer; and (c) hourly mean foF2 measurements from the ionosonde at Huancayo. Comparisons are also made with a self-consistent thermosphere-ionosphere general circulation model and with observationally-based empirical models of winds and density.In concert with the increase in magnetic activity to K_p levels of 5–7, a nighttime (2230 LT) westward intensification of the neutral wind approaching 400 ± 100 ms⁻¹ occurred near the magnetic equator on 22 March 1979, accompanied by a 35% increase in neutral mass density. About 2 h after each of two substorm commencements associated with periods of southward IMF, ∼100γ and ∼200γ reductions in the daytime Huancayo H-component (corrected for ring current effects) are interpreted in terms of ∼0.5 and ∼1.0 mVm⁻¹ westward perturbation electric fields, respectively. An intervening 2-hour period of northward IMF preceded a positive equatorial magnetic perturbation of about 200γ. Time scales for field variations are a few hours, suggesting that processes other than Alfven shielding are involved. Variations in f₀F2 (∼ ± 1.0 MHz) over Huancayo are consistent with the inferred electric fields and magnetic variations. Similar equatorial perturbations are found through examination of other magnetic disturbances during 1979.     

Electrodynamic structural features of the auroral zone as deduced using data from a meridional chain of ionospheric and magnetic stations

This study has used ionospheric and magnetic observational data obtained at a meridional chain of stations during the high latitude geophysical experiment ‘Taimir-82’ in the winter of 1982–1983. Mean statistical latitude-time distributions of the occurrence probability of various types of E_s, their blanketing frequency and of the amplitude of geomagnetic field H-variations have been constructed. Based on these distributions and taking the E_s properties into account, an analysis is made of the mutual correspondence of large-scale structures of the auroral ionosphere and ionospheric currents.Ionospheric currents flow mainly in the region of high E-layer ionization. With increasing magnetic activity, the zone of currents and the zone of ionization expand simultaneously toward lower latitudes. The evening eastward electrojet and the morning westward electrojet are localized inside the zone of diffuse auroral precipitation which is responsible for the formation of E_s type r. The equatorial part of the midnight westward electrojet is also located in the zone of diffuse precipitation which coincides also with the region of maximum ionization of the E-layer. The polar part of this electrojet, which extends far into the dusk sector, is located in the zone of discrete auroral precipitation (a type E_s). Whereas there exists in the meridional cross-section quite a definite relationship between the Harang discontinuity and ionospheric parameters, such a relationship is not manifested in the zonal cross-section of the Harang discontinuity.     

Analytical modelling of the electric field distribution in the high-latitude ionosphere

An analytical approach is implemented for self-consistent modelling of the high-latitude convection electric field. Input parameters are determined as distributions of field-aligned currents and height-integrated conductivity. The high-latitude ionosphere is approximated with an arbitrary number (N) of concentric rings. The height-integrated conductivity (∑) is independent of co-latitude within any ring, but depends on the longitude ~ sin λ. The field-aligned currents flow only along the boundaries of each ring and are presented by Fourier series in longitude. The analytical solution for the potential φ as a function of longitude is also presented as a Fourier series. An analytical solution is obtained for the potential dependencies on co-latitude. For the extreme case, when the integrated conductivity does not depend on longitude, this solution coincides with the analytical results, obtained by other authors. Based on this solution, the potential distribution in the high-latitude ionosphere, an example with N = 5 is shown, the values of conductivity and field-aligned currents being similar to those values used by other authors.     

Diurnal propagating tides in the low-latitude middle atmosphere

Using an equivalent gravity wave f-plane model it is shown that longitude variations in diurnal insolation absorption by tropospheric H₂O can account for longitudinal variations of at least ± 12–15% about zonal mean values in the diurnal wind amplitude at low latitudes (0–20°) between 80 and 100 km, by virtue of the non-migrating propagating tidal modes which are excited. Phase variations of about ± 0.75 h also occur. These percentage variations are conservative estimates, since the background migrating (1,1,1) mode appears to be slightly (20–25%) overestimated in amplitude. In addition, the assumed eddy dissipation values, which appear necessary to model the breaking (1,1,1) mode, are larger than generally considered ‘reasonable’ by photochemical modellers. For a photochemically more reasonable eddy diffusion profile, estimates of longitude differences in diurnal wind amplitude are quite similar to the above values below 87 km, but increase to ± 17–25% near 100 km, with accompanying phase variations of ± 1–2 h about zonal mean values. In addition, it is shown that radiative damping by CO₂ parameterized by a scale-dependent Newtonian cooling coefficient accounts for no more than a 20% reduction in the amplitudes of diurnal propagating tides above 80 km.     

Longitudinal differences and inter-annual variations of zonal wind in the tropical stratosphere and troposphere

A quantitative assessment has been made of the longitude-dependent differences and the interannual variations of the zonal wind components in the equatorial stratosphere and troposphere, from the analysis of rocket and balloon data for 1979 and 1980 for three stations near ±8.5° latitude (Ascension Island at 14.4°W, Thumba at 76.9°E and Kwajalein at 67.7°E) and two stations near 21.5° latitude (Barking Sands at 159.6°W and Balasore at 86.9°E). The longitude-dependent differences are found to be about 10–20 m s⁻¹ (amounting to 50–200% in some cases) for the semi-annual oscillation (SAO) and the annual oscillation (AO) amplitudes, depending upon the altitude and latitude. Inter-annual variations of about 10 m s⁻¹ also exist in both oscillations. The phase of the SAO exhibits an almost 180° shift at Kwajalein compared to that at the other two stations near 8.5°, while the phase of the AO is independent of longitude, in the stratosphere.The amplitude and phase of the quasi-biennial oscillation (QBO) are found to be almost independent of longitude in the 18–38 km range, but above 40 km height the QBO amplitude and phase have different values in different longitude sectors for the three stations near ±8.5° latitude. The mean zonal wind shows no change from 1979 to 1980, but in the troposphere at 8.5° latitude strong easterlies prevail in the Indian zone, in contrast to the westerlies at the Atlantic and Pacific stations.