Comparison of plasma flow and thermospheric circulation over northern Scandinavia using EISCAT and a Fabry-Perot interferometer

The EISCAT incoherent scatter radar, operating in a full tristatic mode, provided data on the ionospheric plasma drift above northern Scandinavia, during the 24 h period, 11 UT 25 November to 11 UT 26 November 1982. For the hours of darkness, 14 UT until 05 UT, observations of thermospheric winds were made by means of a ground-based Fabry-Perot interferometer (FPI) operated at Kiruna Geophysical Institute (21° E, 68° N). During this period, the radar observations describe well the ebbing and flowing of regions of strong convective ion flow associated with the auroral oval. As individual geomagnetic disturbances occur, the overall ion flow pattern intensifies and moves equatorward. The zonal thermospheric wind observed by the FPI responds rapidly to surges of the local ionospheric convection, while the meridional wind response is slower and apparently to much larger-scale features of the geomagnetic input to the high latitude thermosphere. From the data base, periods of strong heating of the ionospheric ions and of the thermospheric gas can be identified, which can be compared with Joule and particle heating rates deduced from the observations of ionospheric drifts, neutral winds, electron densities and auroral emission rates. A three-dimensional, time-dependent global thermospheric model is used to distinguish local and global features of the thermospheric wind field. Meridional and zonal wind components at 312 km may be theoretically derived from the EISCAT data using an appropriate model (MSIS) for neutral temperature. The EISCAT-derived meridional wind is within about 50 m s⁻¹ of the FPI observations throughout the period of joint observations. The EISCAT-derived zonal wind is systematically larger (by about 50%) than the FPI measurement, but the two independent measurements follow closely the same fluctuations in response to geophysical events until 03 UT, when the EISCAT solution is driven away from the FPI measurement by a sharp increase in both neutral and ion temperatures. Between 03 and 05 UT the EISCAT-derived zonal wind is 200–400 m s⁻¹ westward. Allowance for the neutral temperature rise would reduce the EISCAT values towards the very small zonal winds shown by the FPI during this period. We describe the relatively straightforward analysis required to derive the meridional wind from the radar data and the limitations inherent in the derivation of zonal wind, using the ion energy equation, due to the lack of precise knowledge of the background neutral temperature from the EISCAT data alone. For analysis of EISCAT ion drift observations at 312 km, the ground-based FPI temperature measurements do not improve the accuracy of the analysis, since the median altitude of the FPI measurement is probably in the range 180–240 km throughout the observation period. This median altitude and the temperature gradient both fluctuate in response to local geomagnetic events, while the temperature gradient may be considerably greater than that predicted by standard atmospheric models. When the neutral temperature is well known, or when there is a large enhancement of the ion temperature, the EISCAT-derived zonal wind exceeds the FPI measurement, but the consistency with which they correlate and follow ion-drag accelerations suggests that the differences are purely due to the considerable altitude gradients which are predicted by theoretical models.     

The trouble with thermospheric vertical winds: geomagnetic,seasonal and solar cycle dependence at high latitudes

The vertical wind component is frequently used to determine the zero-velocity baseline for measurements of thermospheric winds by Fabry-Perot and other interferometers. For many of the upper atmospheric emission lines from which Doppler shifts are determined, for example for the OI 630 nm emission, available laboratory sources are not convenient for long-term use at remote automatic observatories. Therefore, the assumption that the long-term average vertical wind is zero is frequently used to create a baseline from which the Doppler shifts corresponding with the line-of-sight wind from other observing directions can then be calculated. A data base consisting of 1242 nights of thermospheric wind measurements from Kiruna (68°N, 20°E), a high-latitude site, has been analysed. There are many interesting short-term fluctuations of the vertical wind which will be discussed in future papers. However, the mean vertical wind at Kiruna also has a systematic variation dependent on geomagnetic activity, season and solar cycle. This means that the assumption that the average value of the vertical wind is zero over the observing period cannot be used in isolation to determine the instrument reference or baseline. Despite this note of caution, even within the auroral oval, the assumption of a zero mean vertical wind can be used to derive a baseline which is probably valid within 5 ms⁻¹ during periods of quiet geomagnetic activity (K_p < 2), near winter solstice. During other seasons, and during periods of elevated geomagnetic activity, a systematic error in excess of 10 ms⁻¹ may occur.     

The effects of neutral winds on the propagation of medium-scale atmospheric gravity waves at mid-latitudes

The influence of neutral winds on the propagation of medium-scale atmospheric gravity waves at mid-latitudes is investigated. A 3-dimensional neutral wind model is developed and used together with an atmospheric model in a gravity-wave ray-tracing analysis. It is demonstrated that the thermospheric wind can act as a filter for waves travelling at unfavorable angles to the mean flow, via the mechanisms of reflection and critical coupling. This wind filtering action rotates clockwise diurnally through 360° in the northern hemisphere. Observational evidence is presented which supports these predictions. Extensive modelling indicates that (a) faster and longer period waves are least affected by the neutral winds and (b) fixed-height (e.g. HF Doppler) observations of medium scale gravity waves is only likely to be possible for waves generated locally (within 500–1000 km). Waves generated at greater distances are probably dissipated before reaching the observation region.     

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.     

Neutral density cells in the high latitude thermosphere—1. Solar maximum cell morphology and data analysis

A new high latitude thermospheric neutral density structure has been revealed in NCAR-TIGCM simulations at 120–350 km altitude. The structure consists of density cells above 50° latitude with radii of approximately 1000 km. There are between two to four cells present depending on the altitude and magnetic activity. For example, at 200 km under magnetically active conditions, the density structure consists of four cells: low density cells are located near dawn and dusk and high density cells are located near noon and midnight. Density variations among the cells range from 5 to 50% for magnetically quiet and active conditions respectively. The cells are present at all seasons, for a wide range of magnetic activity levels, and at solar minimum and solar maximum. The density cell morphology is established for equinox solar maximum as a function of altitude and magnetic activity. Departures of the cell structure from this morphology due to seasonal and solar cycles are discussed. The cell morphology provides a new framework in which to interpret lower thermospheric density data. Data to test and confirm the model predictions were provided by the SETA-1 satellite.     

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.     

Seasonal variations of equatorial night-time thermospheric meridional winds

Using h'F data at two equatorial stations, night-time equatorial thermospheric meridional winds have been deduced for a period of two years to study their seasonal characteristics. It has been found that the thermospheric wind shows trans-equatorial flow from summer to winter hemisphere. During equinoxes the flow is mainly equatorward with a reversal to poleward direction around midnight hours. The abatement and reversal of equatorward wind which is weaker in summer compared to equinoxes is attributed to Midnight Temperature Maximum (MTM). The results of the present investigation are compared with those at other equatorial stations and also with the empirical model of Hedin et al. (1991).     

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.     

Scintillation at high latitudes during winter magnetic storms

Scintillation observations are described which were made at Kiruna in northern Sweden during three magnetic storm periods in the winter of 1984–1985. The results were obtained using transmissions from the multisatellite NNSS system, so that it has been possible to chart the development of scintillation activity over some 20° of geomagnetic latitude as a function of time for several days throughout each storm. A region of strong scintillation at the highest latitudes near magnetic noon is a common feature on all but the quietest days. This feature, probably associated with soft particle precipitation into the cusp, shows an abrupt boundary which moves equatorwards as the disturbance develops. In the magnetic midnight sector two latitudinally separate zones of scintillation are found, patchy at high latitudes although more sustained in the auroral zone. An absence of auroral scintillations around midnight UT can be followed by prolonged intense scintillation activity at auroral latitudes during the early morning hours on some disturbed days.     

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.     

Neutral density cells in the high latitude thermosphere—2. Mechanisms

NCAR-TIGCM simulations predict mesoscale cellular structures in the high latitude neutral density at altitudes from 120–350 km. During magnetically active conditions, the density structure at 200 km consists of low-density cells near dawn and dusk and high-density cells near noon and midnight. Mechanisms causing the structured density cells are a result of thermosphere-ionosphere coupling and can be explained in terms of dynamic meteorology. For example, at high latitudes ion drag causes the neutral circulation to flow cyclonically in the dawn sector and anticyclonically in the dusk sector. Low densities are contained within the cyclonic circulation at all altitudes. Below about 170 km, the densities inside the anticyclonic flow are high, while above that altitude densities within the anticyclonic flow are low. While typical dynamic meteorology explains low densities in the centre of cyclonic circulation and high densities inside anticyclonic circulation, the dusk low-density cell in the centre of anticyclonic flow is unexpected. The anticyclonic dusk low-density cell is explained by anomalous antibaric flow due to high-speed winds. 120 km and 200 km altitudes are used to demonstrate the relationship between the high latitude densities and winds as well as the effect of joule heating and auroral particle precipitation on the density structures.     

On the IRI model predictions for the low-latitude ionosphere

Measurements of ionospheric electron density vertical profiles, carried out at a magnetic equatorial station located at Fortaleza (4°S, 38°W; dip latitude 2°S) in Brazil, are analyzed and compared with low-latitude electron density profiles predicted by the International Reference Ionosphere (IRI) model. The analysis performed here covers periods of high (1979/1980) and low (1986) solar activities, considering data obtained under magnetically quiet conditions representative of the summer, winter and equinox seasons. Some discrepancies are found to exist between the observed and the IRI model-predicted ionospheric electron density profiles. For high solar activity conditions the most remarkable one is the observed fast upward motion of the F-layer just after sunset, not considered in the IRI model and which precedes the occurrence of nighttime ionospheric plasma irregularities. These discrepancies are attributed mainly to dynamical effects associated with the low latitude E × B electromagnetic plasma drifts and the thermospheric neutral winds, which are not satisfactorily reproduced either in the CCIR numerical maps or in the IRI profile shapes. In particular, the pre-reversal enhancement in the vertical E × B plasma drifts around sunset hours has a great influence on the nighttime spatial distribution of the low-latitude ionospheric plasma. Also, the dynamical control exerted by the electromagnetic plasma drifts and by the thermospheric neutral winds on the low-latitude ionospheric plasma is strongly dependent on the magnetic declination angle at a given longitude. These important longitudinal and latitudinal dependences must be considered for improvement of IRI model predictions at low latitudes.     

Some further results on long term mesospheric and lower thermospheric wind observations by the Arecibo UHF radar

A second series of long term mesospheric and lower thermospheric wind observations was conducted at Arecibo (18.4°N, 66.8°W) between 6 and 20 March 1981 using the UHF Doppler radar, following the first observations in August 1980 (Hirota et al., 1983). Zonal and meridional wind velocities were measured during the morning (8–10 LT) and afternoon (13–15 LT) periods. The mean wind profile averaged over the entire observational period shows the predominance of the diurnal tide. The fluctuating wind vector rotates clockwise relative to height with a characteristic vertical scale of about 10 km. The phase difference inferred by a cross correlation analysis between morning and afternoon profiles indicates that the dominant period is about 20–30 h. This oscillation is discussed in relation to internal inertia-gravity waves observed by the same radar in the lower stratosphere. On the other hand, wind fluctuation with a vertical scale larger than 20 km shows a substantial day-to-day variation with a period of 5–8 days. This long period oscillation shows a good correlation with the global scale geopotential height anomalies at 1 mb (46–48 km) observed by the Tiros-N satellite at 20°N. Our evidence suggests that westward travelling planetary-scale waves with zonal wavenumber one may propagate up to the lower thermosphere.     

Mean winds at 60–90 km observed with the MU radar (35°N)

Mean winds at 60–90 km altitudes observed with the MU radar (35°N, 136°E) in 1985–1989 are presented in this paper. The zonal wind at 70 km became westward and eastward in summer and winter, respectively, with a maximum amplitude of 45 m s⁻¹ westward in early July and 80 m s⁻¹ eastward at the end of November. The meridional wind below 85 km was generally northward with the amplitudes less than 10 m s⁻¹. In September to November, the meridional wind at 75–80 km becomes as large as 20–30 m s⁻¹. Those zonal wind profiles below 90 km show good coincidence with the CIRA 1986 model, except for the latter half of winter, from January to March, when the observational result showed a much weaker eastward wind than the CIRA model. The height of the reversal of the summer wind from westward to eastward was determined as being 83–84 km, which is close to the CIRA 1986 model of 85 km. The difference between the previous meteor radar results at 35–40°N, which showed the reversal height below 80 km, could be due to interannual variations or the difference in wind measurement technique. In order to clarify that point, careful comparative observations would be necessary. These mean winds were compared with Adelaide MF radar observations, and showed good symmetry between the hemispheres, including the summer reversal height, except for the short period of eastward winds above Kyoto and the long period over Adelaide.     

Spectral characteristics of high frequency waves backscattered by small scale F-region irregularities: evidence of strong sub-auroral ion flow

The spectra of high frequency waves backscattered at night by small scale (10–20 m) sub-auroral F-region irregularities often exhibit large Doppler shifts and widths in the local time sector 2000–2400. After local midnight the Doppler shifts and the widths of the spectra decrease rapidly. We present examples of experimental data, obtained with the two coherent backscatter radars of the EDIA¹ experiment, showing the spectral characteristics just mentioned. From the Doppler shift measured at the two sites we deduced the perpendicular velocity of the irregularities, which can reach values as high as 2000 ms ⁻¹. These observations are interpreted using results of theoretical models which predict strong sub-auroral ion flow in the trough region.     

Observations of winds in the Antarctic summer mesosphere using the spaced antenna technique

Observations of winds in the 60–100 km height range were made at Mawson (68°S, 63°E) during December 1981 and January 1982 with the MF spaced antenna technique. The prevailing winds are in accord with other recent observations made at high latitudes and show a peak in the zonal wind near 80 km with westward winds of 30 m s ⁻¹. The meridional winds maximize near 90 km with an equatorward flow of 10 m s⁻¹. The diurnal tidal components are in reasonable agreement with recent model predictions, especially in phase. The amplitudes tend to be larger than the model values. The semidiurnal tide is not as stable as the diurnal tide and shows evidence for interference effects between different modes.     

Observations of neutral winds in the polar cap during northward IMF

Long term remote observations of neutral winds at F-region altitudes have been performed at Thule Air Base (lat. 76.5°N, long. 69.0°W), Greenland, and Søndre Strømfjord (lat. 67.0°N, long. 50.9°W), Greenland. The former site is very close to the geomagnetic pole, while the latter site is within the polar cap for several hours each night on either side of geomagnetic midnight. Wind data corresponding to clear sky conditions and Kp ⩽ 4 were sorted according to the sign of the IMF B_z component. The averaged maximum poleward flow near midnight LST was reduced by approximately one third during B_z northward conditions. If the magnitude of B_y was less than the magnitude of the northward B_z component, then the averaged poleward flow was further reduced by one half. In addition, if B_z > 5 nT, then sunward directed horizontal neutral winds were observed at the very highest latitudes near noon LST.     

Summertime stratospheric wind measurements above the South Pole

Observations of the mean wind flow and wave motions in the stratosphere at the South Pole are presented. The atmospheric motions are determined from the tracking of a high altitude, zero-pressure balloon launched from Amundsen-Scott Station during the austral summer of 1985–1986. The balloon position was precisely monitored by an optical theodolite for a large portion of the flight so that small scale motions could be resolved. The mean flow above the pole was approximately 3ms⁻¹. Atmospheric motions characteristic of internal gravity waves were observed with an intrinsic period of approximately 4.5 h and vertical and horizontal wavelengths of approximately 2.5km and 125km, respectively. The horizontal perturbation velocity of the observed waves was large compared to the mean horizontal flow velocity. The implication is that wave motions play a dominant role in the transport of stratospheric constituents in regions where the mean winds are light, such as over the South Pole during austral summer.     

Neutral wind measurements over Waltair using meteor wind radar

A coherent pulse Doppler Radar System at Waltair (17°43′N, 83°18′E) has been used to measure neutral winds in the altitude range 80–110 km. Preliminary results obtained during July and August 1979 are presented. The diurnal component is found to be the dominant tidal mode for most of the time. The results are compared with meteor wind observations at other stations.     

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.