Plasma drifts inferred from thermospheric neutral winds and temperature gradients observed at low latitudes

Low-latitude plasma drifts (zonal and meridional) in the F-region are inferred from the observed night-time thermospheric neutral wind velocities and temperature gradients, together with models for the neutral density (MSIS-86 model) and the electron density (IRI model). The thermospheric neutral winds and temperatures are derived from measurements of Doppler shifts and widths of the Oi 630.0 nm airglow emission line, respectively, using a Fabry-Perot interferometer at Cachoeira Paulista (23°S, 45°W), Brazil. The equations considered are the ideal gas law and the momentum equation for the thermosphere, which includes the time variation of the neutral wind, the pressure gradient which is related to the temperature and density gradients and the ion drag force. The present method to infer the night-time plasma drift using observed neutral parameters (time variation of neutral wind velocities and temperature gradients) showed results that are in reasonable agreement with our calculated plasma drifts and those observed in other low-latitude locations. On the other hand, it is surprising that sometimes the winds flow from the observed coldest sector to the hottest part of the thermosphere during many hours, suggesting that plasma drift can drive the neutral winds at low latitudes for a period of time.     

The isobaric F2-layer

Cyclic diagrams, obtained by plotting the daily variation of the ionospheric electron density NmF2 against the height hmF2, are drawn for typical conditions at Slough (52°N) and Watheroo (30°S). Using the MSIS86 thermospheric model to relate the heights hmF2 to values of atmospheric pressure, the F2-peak is found to lie at nearly the same pressure-level at any given local time, over a wide range of geophysical conditions (season, solar cycle, magnetic disturbance). As local time varies, the pressure level corresponding to hmF2 varies in a way that is mainly determined by the local time variation of the thermospheric winds. This is verified for noon and midnight, using the MSIS86 model to compute the winds. The noon values of peak electron density (NmF2) are fairly consistent with theory (using values of solar ionizing flux as quoted in the literature), but with some discrepancies—particularly at sunspot maximum—that are probably due to uncertainties in chemical composition, or to the effects of vibrational excitation of molecular nitrogen. Overall, the analysis shows a remarkable consistency between ionospheric theory, the data and the MSIS model.     

Modelling studies of the conjugate-hemisphere differences in ionospheric ionization at equatorial anomaly latitudes

The relative importance of the equatorial plasma fountain (caused by vertical E x B drift at the equator) and neutral winds in leading to the ionospheric variations at equatorial-anomaly latitudes, with particular emphasis on conjugate-hemisphere differences, is investigated using a plasmasphere model. Values of ionospherec electron content (IEC) and peak electron density (Nmax) computed at conjugate points in the magnetic latitude range 10–30° at longitude 158°W reproduce the observed seasonal, solar activity, and latitudinal variations of IEC and Nmax, including the conjugate-hemisphere differences. The model results show that the plasma fountain, in the absence of neutral winds, produces almost identical effects at conjugate points in all seasons; neutral winds cause conjugate-hemisphere differences by modulating the fountain and moving the ionospheres at the conjugate hemispheres to different altitudes.At equinox., the neutral winds, mainly the zonal wind, modulate the fountain to supply more ionization to the northern hemisphere during evening and night-time hours and, at the same time, cause smaller chemical loss in the southern hemisphere by raising the ionosphere. The gain of ionization through the reduction in chemical loss is greater than that supplied by the fountain and causes stronger premidnight enhancements. in IEC and Nmax (with delayed peaks) in the southern hemisphere at all latitudes (10–30°). The same mechanism, but with the hemispheres of more flux and less chemical loss interchanged, causes stronger daytime IEC in the northern hemisphere at all latitudes. At solstice, the neutral winds, mainly the meridional wind, modulate the fountain differently at different altitudes and latitudes with a general interhemispheric flow from the summer to the winter hemisphere at altitudes above the F-region peaks. The interhemispheric flow causes stronger premidnight enhancements in IEC and Nmax and stronger daytime Nmax in the winter hemisphere, especially at latitudes equatorward of the anomaly crest. The altitude and latitude distributions of the daytime plasma flows combined with the longer daytime period can cause stronger daytime IEC in the summer hemisphere at all latitudes.     

Atmospheric winds calculated from diurnal changes in the mid-latitude ionosphere

Diurnal variations in the electron content (N_t) and peak density (N_m) of the ionosphere are calculated using a full time-varying model which includes the effects of electric fields, interhemispheric fluxes and neutral winds. The calculation is iterated, adjusting the assumed hourly values of neutral wind until a good match is obtained with mean experimental values of N_t and N_m. Using accurate ionospheric data for quiet conditions at 35°S and 43°S, winds are derived for summer, equinox and winter conditions near solar maximum and solar minimum. Solar maximum results are also obtained at 35°N. Changes in the neutral wind are found to be the major cause of seasonal changes in the ionosphere, and of differences between the two hemispheres. Calculated winds show little variation with latitude, but the winds increase by about 30% at solar minimum (in equinox and winter). The HWM90 wind model gives daytime winds which are nearly twice too large near solar maximum. The theoretical VSH model agrees better with observed daytime variations, and both models fit the observed winds reasonably well at night. Results indicate that modelling of the quiet, mid-latitude ionosphere should be adequate for many purposes when improved wind models are available. Model values for the peak height of the ionosphere are also provided; these show that wind calculations using servo theory are unreliable from sunrise to noon and for several hours after sunset.     

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.     

The dynamics of the poleward edge of the trough deduced using data from a meridional chain of vertical-incidence sounding stations

Data from four ionospheric stations located along the 902E meridian in the range 55–702 of corrected geomagnetic latitude, were used to construct latitude-time electron density distributions in the F2-layer peak for 17 winter nights of 1982–1983. It is concluded that under stationary convection conditions the poleward edge of the trough during the nighttime displaces only 0.5-l°/h, that is, significantly less than obtained from existing analytical models of the ionospheric trough. When the stationarity is upset (due to the development of a substorm or abrupt changes of the north-south component of the IMF), the poleward edge of the trough is observed to displace abruptly equatorward. In the substorm expansion phase these displacements can amount to 4–5° in less than an hour. Such displacements of the poleward edge in the evening hours can characterize the dynamics of the inner edge of the plasma sheet.     

Effect of the electric field on the equatorial ionospheric plasma distribution

It is known that on a counter electrojet day the noontime electron density at the equator shows enhanced values with no bite-out. The consequences of the absence of the normal equatorial electrojet on the electron density distribution at the equatorial station Kodaikanal (dip latitude 1.4°N, long. 77.5°E) and at an anomaly crest location Ahmedabad (dip latitude 18°N, long. 73°E) are discussed for a strong electrojet (SEJ) day and a counter electrojet (CEJ) day. The electron density distribution with height for a pair of SEJ and CEJ days at the two equatorial stations Kodaikanal and Huancayo (dip latitude 1°N, long. 75°W) are studied. The F-region peak height, hm and the semi-thickness parameter ym on the SEJ day followed a similar variation pattern. On the CEJ days ym exhibited a substantially low and mostly flattened daytime variation compared to the peaked values on the SEJ day. An attempt is made to interpret these differences in terms of the changes in the vertical drift pattern resulting from the E × B drift of plasma at the equator and the varying recombination rate β, which is also a height dependent and a local time dependent parameter.     

Can the high latitude ionosphere support large field-aligned ion drifts?

Ion velocities perpendicular and parallel to the geomagnetic field have recently been deduced by Smith et al. from bistatic measurements at 71° geomagnetic latitude in the afternoon sector. The results of this experiment include large (>400 m s⁻¹) downward ion velocities parallel to the magnetic field that persist for hours, small (100 m s⁻¹) ion velocities perpendicular to the magnetic field and electron density profiles with extremely narrow full-width at half-maximum. The explanation of these results was that the ionospheric flux tubes observed were near the terminator, and thus, sunlit at the top and in darkness at the bottom. The difference in production between the top and bottom of the flux tube creates an excess of ions at the top, which rapidly diffuse downwards. A three-dimensional, time-dependent model of the ionosphere has been used to test this explanation. Numerical experiments were performed to determine upper limits for the downward ion velocity. Assuming reasonable vertically-induced ion drifts due to either neutral winds or plasma convection, these upper limits were substantially smaller than the measurements. The location of the terminator was found to contribute a maximum of about 60 m s⁻¹ to the vertical ion velocity due to diffusion in a partially illuminated flux tube. In an attempt to explain the narrow density profiles without invoking an additional ionization source, the downward force in the model was arbitrarily increased, as would occur due to parallel electric fields in the ionosphere. Since the interpretation of these measurements as large field-aligned flows seems untenable by a model thought to be consistent with the currently accepted physics of the atmosphere, an alternate hypothesis is presented. If the common volume measurement is made in a region of O⁺ precipitation, then the line profile would not be Doppler shifted when viewed off-zenith. Therefore, the field-aligned velocities would be small, and the narrow width of the profiles would be due to enhanced electron densities in an O⁺ arc.     

Optical interferometer measurements of thermospheric temperature at Kavalur (12.5°N, 78.5°E), India

First results on the behaviour of thermospheric temperature over Kavalur (12.5°N, 78.5°E geographic; 2.8°N geomagnetic latitude) located close to the geomagnetic equator in the Indian zone are presented. The results are based on measurements of the Doppler width of O(¹D) night airglow emission at 630 nm made with a pressure-scanned Fabry-Perot interferometer (FPI) on 16 nights during March April 1992. The average nighttime (2130-0430 IST) thermospheric temperature is found to be consistently higher than the MSIS-86 predictions on all but one of the nights. The mean difference between the observed nightly temperatures and model values is 269 K with a standard error of 91 K. On one of the nights (9/10 April 1992, Ap = 6) the temperature is found to increase by ~250 K around 2330 IST and is accompanied by a ‘midnight collapse’ of the F-region over Ahmedabad (23°N, 72°E, dip 26.3°N). This relationship between the temperature increase at Kavalur and F-region height decrease at Ahmedabad is also seen in the average behaviour of the two parameters. The temperature enhancement at Kavalur is interpreted as the signature of the equatorial midnight temperature maximum (MTM) and the descent of the F-region over Ahmedabad as the effect of the poleward neutral winds associated with the MTM.     

Variability of winds and temperatures in the lower thermosphere

High-resolution daytime incoherent scatter radar measurements of plasma temperatures and drifts in the ionospheric E-region above Millstone Hill (42.6°N, 71.5°W) have been used to derive horizontal neutral winds and temperatures in the lower thermosphere (105–130 km) during five multi-day campaigns in 1987–1991. The underlying semi-diurnal tidal component has been determined from the observations, with characteristic average amplitudes of 50 ± 15 m/s and 30 ± 10 K. Phase propagation with altitude follows the expected structure of semi-diurnal tidal modes, but reveals complex coupling of tidal modes, particularly above 115 km. Day-to-day variability in the winds and temperatures is large, and the deviations from the semi-diurnal harmonic can exceed 40 m/s and 50 K. No strong correlations have so far been found with geophysical parameters to explain the observed variability.     

The effect of the electric field and neutral winds on E-region ion densities and conductivities at low latitudes

A model to calculate electron densities and electrical conductivities in the ionospheric E-region at low latitudes has been developed. Calculations have been performed under photochemical equilibrium and including plasma transport due to the electric field and neutral winds. Results have been compared with observations at Arecibo (18.15°N, 66.20°W), Thumba (8°32′N, 76°51′E) and SHAR (14.0°N, 80.0° E). Good agreement is obtained for Arecibo. For Thumba and SHAR agreement is satisfactory for altitudes above 110 km. Below 100 km, model predictions are too low in comparison with the observed data. The effect of plasma transport on electron densities and Hall and Pedersen conductivities is investigated in detail. A combination of neutral winds and a downward (or westward) electric field can compress the plasma into a thin layer. An upward electric field along with the neutral winds gives rise to a broad, multilayered profile. The ratio of height-integrated Hall to Pedersen conductivities changes from 1.2 to 2 in some cases.     

Thermosphere dynamics

Some recent investigations of thermosphere dynamics, carried out in the U.S.S.R., are reviewed briefly. The global empirical models of thermospheric motions are obtained on the basis of ground-based HF and meteor radar measurements of ionospheric irregularities drifts. Numerical modelling of large scale thermospheric electrodynamics for the low and mid-latitudes for quiet geomagnetic conditions, is presented. Disturbances of thermospheric wind systems from high latitude heat sources are considered. The response of lower thermosphere dynamics due to variations of solar and geomagnetic activity are discussed.     

North-south asymmetry in the ionospheric equatorial anomaly in the African and the West Asian regions produced by asymmetrical thermospheric winds

The effect of asymmetrical thermospheric winds on NmF2 at the dip I = 30° and its magnetic conjugate point have been computed for equinox conditions to study asymmetry in the ionospheric equatorial anomaly in the African and West Asian regions. The wind models of I11 et al. and Chan and Walker have been used in our computations. During the daytime, due to the winds NmF2 in the northern crest becomes greater than NmF2 in the southern crest; at night the reverse is true in both regions. It is shown that the observed asymmetry in NmF2 at the equatorial crest in the African sector can be well explained by considering the effects of asymmetrical winds with respect to those in the West Asian sector.     

Simulations of the seasonal variations of the thermosphere and ionosphere using a coupled,three-dimensional,global model,including variations of the interplanetary magnetic field

The University College London Thermospheric Model and the Sheffield University Ionospheric Convection Model have been integrated and improved to produce a self-consistent coupled global thermospheric/high latitude ionospheric model. The neutral thermospheric equations for wind velocity, composition, density and energy are solved, including their full interactions with the evolution of high latitude ion drift and plasma density, as these respond to convection, precipitation, solar photoionisation and changes of the thermosphere, particularly composition and wind velocity. Four 24 h Universal Time (UT) simulations have been performed. These correspond to positive and negative values of the IMF BY component at high solar activity, for a level of moderate geomagnetic activity, for each of the June and December solstices. In this paper we will describe the seasonal and IMF reponses of the coupled ionosphere/thermosphere system, as depicted by these simulations. In the winter polar region the diurnal migration of the polar convection pattern into and out of sunlight, together with ion transport, plays a major role in the plasma density structure at F-region altitudes. In the summer polar region an increase in the proportion of molecular to atomic species, created by the global seasonal thermospheric circulation and augmented by the geomagnetic forcing, controls the plasma densities at all Universal Times. The increased destruction of F-region ions in the summer polar region reduces the mean level of ionization to similar mean levels seen in winter, despite the increased level of solar insolation. In the upper thermosphere in winter for BY negative, a tongue of plasma is transported anti-sunward over the dusk side of the polar cap. To effect this transport, co-rotation and plasma convection work in the same sense. For IMF BY positive, plasma convection and co-rotation tend to oppose so that, despite similar cross-polar cap electric fields, a smaller polar cap plasma tongue is produced, distributed more centrally across the polar cap. In the summer polar cap, the enhanced plasma destruction due to enhancement of neutral molecular species and thus a changed ionospheric composition, causes F-region plasma minima at the same locations where the polar cap plasma maxima are produced in winter.     

A new role for the plasmaphere in “quiet” ionospheric electrodynamics

Global scale longitudinal gradients of pressure in the plasmasphere may be formed naturally by ionospheric processes, or caused by electrostatic fields of ionospheric dynamo origin. It is shown that plasmaspheric gradients of pressure, orthogonal both to the magnetic field (B) and to grad B, generate geophysically significant field-aligned currents. Considering the ionosphere and plasmasphere as a coupled electrodynamic system, these currents alter non-negligibly the self-consistent ionospheric electric field and current. Criteria are established for this coupling mechanism (a kind of plasmaspheric impedance) to be significant. This has implications for the relationships of ionospheric electric fields and currents, F-region drifts, and magnetic variations, due to upper atmosphere tides and winds.     

Solar activity effects on equatorial plasma bubble zonal velocity and its latitude gradient as measured by airglow scanning photometers

Using a new mode of scanning 630-nm photometer operation the zonal velocities of ionospheric plasma depletions were measured over Cachoeira Paulista in Brasil in two east-west planes tilted 30°N and 30° S with respect to zenith. The measurements cover a time period of approximately 2 years, from January 1988 to January 1990, a period marked by significant increase in solar activity of the ongoing cycle. The results have permitted a rather detailed evaluation of the local time and latitude variations in the zonal plasma bubble velocity as a function of solar activity. Although the mean trend in the velocity local time variation is a decrease from early evening to post-midnight hours, a strong tendency for velocity peaks is observed near 21 LT and midnight. The velocities as well as their height (latitude) gradients show perceivable increases with solar activity represented as sunspot numbers. The present results are compared with the ambient plasma velocities measured using the Jicamarca radar by Fejer el al. (1985), J. Geophys. Res. 90, 12249, with that measured on board the DE 2 satellite on the equatorial latitudes by Coley and Heelis (1989), J. geophys. Res. 94, 6751, and with various theoretical calculations, in an attempt to bring out the salient features of the plasma dynamics of the equatorial ionosphere.     

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.     

Computer simulations of the seasonal variations of the ionospheric equatorial anomaly in East Asia under solar minimum conditions

The previous dynamical, computer simulation model of the ionosphere at low latitudes of Chan H. F. and Walker G. O. (1984a, J. atmos. terr. Phys. 46, 1103; 1984b, J. atmos. terr. Phys. 46, 1113) has been modified to (1) include photoionization of molecular species NO⁺, N₂⁺ and O₂⁺ below 300km, (2) decouple the ionization and wind calculations below 180 km and (3) expand the geographical coverage to 46°N-30°S latitude. The first two modifications improved the model stability and the latter reduced the effect of the lateral boundaries on the equatorial anomaly. Results are presented for the representative seasonal months of January, April and July for East Asia, during solar minimum, comprising latitudinal-local standard time (120°E) contour plots of (1) the atmospheric pressure, (2) the computed meridional wind at 300 km, (3) the foF2 and (4) hmF2, together with latitudinal profiles of foF2 and N_T (electron content) showing the daytime development and nighttime decay of the equatorial anomaly.Comparisons have been made between the computer simulations and various experimental measurements of foF2, M(3000) F2 and N_T obtained in East Asia during periods of low solar activity. Most of the gross features of the development and decay of the equatorial anomaly at the various seasons were reproducible by the model simulations, the best agreement occurring for the equinoctial month of April.     

The effects of ionospheric horizontal electron density gradients on whistler mode signals

Whistler mode group delays observed at Faraday, Antarctica (65° S, 64° W) and Dunedin, New Zealand (46° S, 171° E) show sudden increases of the order of hundreds of milliseconds within 15 minutes. These events (‘discontinuities’) are observed during sunrise or sunset at the duct entry regions, close to the receiver's conjugate point. The sudden increase in group delay can be explained as a tilting of the up-going wave towards the sun by horizontal electron density gradients associated with the passage of the dawn/dusk terminator. The waves become trapped into higher L-shell ducts. The majority of the events are seen during June-August and can be understood in terms of the orientation of the terminator with respect to the field aligned ducts. The position of the source VLF transmitter relative to the duct entry region is found to be important in determining the contribution of ionospheric electron density gradients to the L-shell distribution of the whistler mode signals.     

OI 630 nm imaging observations of equatorial plasma depletions at 16° S dip latitude

Equatorial ionospheric irregularities in the F-layer have been the subject of intensive experimental and theoretical investigations during recent years. The class or irregularities which continues to receive much attention is characterized by large scale plasma depletions, generally referred to as ionospheric plumes and bubbles. The OI 630.0 nm F-region night-glow emissions arising from recombination processes can be used to observe the dynamics of transequatorial ionospheric plasma bubbles and smaller scale plasma irregularities. In a collaborative project between the Center for Space Physics of Boston University and Brazil's National Institute for Space Research (INPE), an all-sky imaging system was operated at Cachoeira Paulista (22.7° S, 45.0° W, dip latitude 15.8° S), between March 1987 and October 1991. In addition to the imager, photometer and VHP polarimeter observations were conducted at Cachoeira Paulista, with ionospheric soundings carried out at both C. Paulista and Fortaleza, the latter at 3.9° S, 38.4° W, 3.7° S dip latitude. For this longitude, the observed seasonal variation of the airglow depletions shows a maximum from October through March and a very low occurrence of airglow depletions from April through September. This long series of OI 630.0 nm imaging observations has permitted us to determine that when there are extended plumes, the altitudes affected over the magnetic equator often exceed 1500 km and probably exceed 2500 km at times, the maximum projection that can be seen from Cachoeira Paulista. This holds true even during years of low solar flux.