Interference of tidal and gravity waves in the ionosphere and an associated sporadic E-layer

Observations made on 10 July 1987 with the EISCAT UHF radar are presented. The F-region measurements of both electron density and field-aligned ion velocity show that an upward propagating gravity wave with a period of about 1 h is present. The origin of the gravity wave is probably auroral. The E-region ion velocities show a tidal wave and both upward and downward propagating gravity waves. The gravity waves have three dominant periods with a possible harmonic relationship and similar vertical wavelengths. These waves are either reflected at a single reflection level, ducted between two levels, or they are generated in a non-linear interaction between gravity and tidal waves. The E-region electron density is dominated by particle precipitation. After a short burst of more intense precipitation, a sporadic E-layer forms at 105km and then disappears 40min later. Within this time, the layer rises and falls by a few kilometres, following closely the motion of a convergent null in the velocity profile. We suggest that the formation and destruction of this layer is controlled by both the precipitation, which indirectly provides a source of metal ions through charge exchange, and the superposition of gravity waves and the tidal wave.     

Density profiles of sporadic E-layers containing two metal ion species

Ionospheric plasma containing two types of metal ions is investigated under the action of the wind shear mechanism or, alternatively, an electric field causing convergent vertical plasma flow. It is shown that the different ion species are separately collected into thin sheets with a height difference ranging from some hundreds of meters to several kilometers. Theoretical density profiles for Mg⁺ and Fe⁺ ions are calculated assuming a screw-like wind structure or a strong auroral electric field. It is found that the two ion layers usually partially merge forming a single E_s-layer. If the height difference of the ion sheets is not too great as compared to their thicknesses, the E_s-profile is single peaked and approximately symmetric. With increasing layer separation the two sheets will gradually be discerned, until finally a double peaked profile is created. It is suggested that some of the observed complexities in E_s-profiles are caused by the presence of more than one monoatomic ion species.     

Incoherent scatter observations of thin ionization layers at Sondrestrom

High resolution incoherent-scatter observations of E-region thin (1–3 km) metallic ion layers are presented. Data were collected during three different periods from August 1990 to August 1991, in three different experimental modes. First, the antenna was directed vertically and the entire duty cycle was devoted to Barker coded multi-pulse [Zamlutti (1980) J. atmos. terr. Phys.42, 975–982] measurements to determine the densities and temperatures in the E-region with 300 m resolution. The second experiment measured the F-region electric field as well as the high resolution E-region densities. For the third experiment the antenna was scanned magnetic north-south while only the E-region densities were measured. The experiments were carried out on 16 different nights for a period of 4 h each night at a time near magnetic midnight. Thin ionization layers were observed on 12 of the 16 nights. The first experiment demonstrated that the thin layers are composed of a significant fraction of heavy metallic ions; assuming the layers are composed of a mixture of Fe⁺ and Mg⁺ a composition estimate of 63% Fe⁺ was obtained in one example. The second experiment investigated the relationship between the direction of the electric field and the presence of the thin layers. In these observations thin layers were only present when the electric field was pointed in the magnetic north-west or south-west quadrants, most frequently when the field was near magnetic west. Correlation between layer altitude and field direction was also observed, layers occurring at higher altitudes for fields directed in the north-west, and lower altitudes for fields directed to the south-west. The observations are compatible with the electric field mechanism for thin ionization layer formation. The scanning experiment showed that the layers were of a limited latitudinal extent, typically about 100 km up to a maximum of about 200 km.     

Weak nonlinear theory of the ionospheric response to atmospheric gravity waves in the F-region

The nonlinear ionospheric response to atmospheric gravity waves is studied in an approximate fashion using a new approach. The concept of nonlinear travelling ionospheric disturbances (TIDs) is outlined, and the nonlinear behaviour of atmospheric gravity waves is calculated. A principal result is that harmonics are generated which cause the wave velocity perturbation to deform. The ionospheric response is investigated by solving the continuity equation for ionization in the F-region. The distortion of the TIDs waveform produced by the nonlinear interactions is depicted. The nonlinear TIDs depart seriously from a cosinusoidal wave described by previous linear TID theory. The distorted TIDs appear as ‘sharp peak’ and ‘sawtooth’ waveform shapes. The ‘peaks’ can be upward or downward, and the ‘sawteeth’ forward or backward, depending on the wave parameters. The nonlinearly distorted TIDs show a good agreement with various observed ionospheric irregularities produced by atmospheric gravity waves.     

Effect of the large-scale convection electric field structure on the formation of thin ionization layers at high latitudes

A three-dimensional simulation of the high-latitude ionosphere was applied to investigate the geographical distribution of E-region thin ionization layers which may be formed by the action of the convection electric field. The simulation model computes the ion densities (O⁺, O⁺₂, N⁺, N⁺₂, NO⁺, Fe⁺), and temperatures as a function of altitude, latitude, and longitude. The stationary state momentum and continuity equations are solved for each ion species, then the energy equation is solved for electrons, neutrals, and a generic ion having the mean ion mass and velocity. The various electric field patterns of the Heppner and Maynard [(1987) J. geophys. Res.92, 4467–4489] convection electric field model were applied and the ionization density pattern was examined after a time sufficient for the formation of thin layers (≈2000 s). It was found that large areas of thin ionization layers were formed for each of the electric field patterns examined. Southward IMF B_z conditions resulted in thin layers forming in the pre-midnight sector in the latitude range north of about 70° to about 80°, and after midnight between 60 and 70°. For northward B_z conditions, the layers were mainly in the pre-midnight sector and covered a latitude range from about 60 to 80°.     

Nocturnal high-latitude E-region in winter during extremely quiet conditions

The quiet night-time E-region at high latitudes has been studied using the EISCAT UHF radar. Data from three subsequent nights during a long period of low magnetic activity are shown and typical features of electron density are described. The background electron density is observed to be 5·10⁹ m⁻³ or smaller. Two types of enhancements above this level are observed ; one is due to charged particle precipitation associated with the F-region trough and the other is composed of sporadic-E layers due to waves in the neutral atmosphere. The sporadic-E is observed to exist almost continuously and to exhibit a regular diurnal behaviour. In addition to the typical afternoon and morning sequential layers, a third major descending layer is formed at night after the passage of the F-region trough The afternoon layer disappears simultaneously with the enhancement of the northward trough-associated electric field and the night-time layer appears at high altitudes after the field has again been reduced to a small value. It is suggested that metal ions from low altitudes are swept by the electric field to the upper E-region where they are again compressed to the night-time layer. A set of steeply descending weaker layers, merging to the main night-time layer are also observed. These layers are most probably caused by atmospheric gravity waves. Theoretical profiles for molecular ions indicate that the strongest layers are necessarily composed of metal ions but, during times when the layers are at their weakest, they may be mainly composed of molecular ions.     

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.     

Formation of ionisation layers responsible for blanketing Es during daytime counter-electrojet over magnetic equator

A model using photochemistry and transport due to electric fields and gravity wave winds has been used to explain the formation of ionisation layers observed over an equatorial station Thumba (dip 0°47′S) with rocket-borne Langmuir probes during two daytime counter-electrojet periods. These layers were seen as blanketing E_s-layers with an ionosonde at Thumba. Convergence of the metallic ions due to three-dimensional gravity wave winds and a westward electric field appears to be mainly responsible for the observed ionisation layer over the equator.     

High-resolution measurements of magnetospheric electric fields

Observations made at EISCAT suggest that the plasma velocity measured in the F-region above Tromsø can vary substantially on a timescale of a minute or so. The high-resolution measurements made using alternating codes during the ERRRIS experiment have confirmed this result by showing that the rapid variations of plasma velocity measured directly correspond exactly to the variations of ion temperature in the rmupper-E and lower-F region caused by frictional heating, and the variations of electron temperature in the E-region, caused by wave turbulence heating.     

Thermospheric wind measurements with EISCAT

Thermospheric wind measurements with the EISCAT UHF radar around the evening Harang discontinuity are presented both in the E- and F-layers. Within the E-layer auroral oval the Lorentz and Coriolis force are shown to be more or less in balance. The neutral velocity is a factor of the order of two smaller than the ion velocity and is on average advanced 90° in a clockwise direction compared to the ion velocity. In the low electron density region just before the Harang discontinuity and outside the auroral oval a large (~250 m s⁻¹), thermally dominated neutral wind is closely followed by the ion wind in the antisolar direction. There is also a large downward flow present just before the Harang discontinuity. In the F-layer the neutral wind approximately follows the ion convection pattern, except for a couple of hours after the sudden change in the ion convection just after the passage of the evening Harang discontinuity. The close resemblance between the equilibrium ion and neutral flow when the neutral-ion collision frequency is close to twice the Earth's angular velocity may be connected to back pressures created by Joule heating in the case of an appreciable ion-neutral velocity difference.     

On the role of auroral electric fields in the formation of low altitude sporadic-E and sudden sodium layers

The characteristics of metallic and molecular ion sporadic-E (E_s) layers, formed by the action of strong electric fields at auroral latitudes, are examined using computer simulations. It is found that, for electric fields directed between northward and westward (northern hemisphere), thin metallic ion layers (<2 km thick) can be formed above about 105 km altitude. For electric fields directed from westward, through southward, to south-eastward, slightly thicker (4–6 km thick) metallic ion layers can form between 90 and 105 km altitudes. Thin layers of molecular ions can be formed by electric fields directed between north and west if the ion density is low. Examples of E_s layers observed by the EISCAT radar, together with simultaneous observations of electric fields and ion drifts are presented which show good agreement with the simulations. The relationship between the lower-altitude E_s layers and sudden sodium layers (SSLs) is discussed leading to an explanation of some of the characteristics of SSLs at high latitude. A possible involvement of smoke particles in the formation of both E_s layers and SSLs is proposed.     

A comparison of three methods of measuring tidal oscillations in the lower thermosphere using EISCAT common programmes

The EISCAT Common Programme can be used in three ways to monitor tidal oscillations in the lower thermosphere. In Common Programme One (CPI) tristatic observations provide measurements of the ion-velocity vector at several heights in the E-region and one height in the F-region. In Common Programme Two (CP2) monostatic measurements give profiles of ion velocity in the E-region while tristatic measurements give continuous measurements of ion velocity in the F-region. From the ion velocities and the ion-neutral collision frequency, the vector of the E-region neutral wind can be determined and both east-west and north-south components of the diurnal, semi-diurnal and ter-diurnal oscillations can be identified. CP1 and CP2 also provide profiles of the field-aligned ion velocity, and these can be used to calculate the north-south component of the neutral wind without knowing the ion-neutral collision frequency, but the result is affected by any vertical component of neutral velocity. The three methods are compared and the advantages of CP2 demonstrated.     

Phase velocity studies of 34-cm E-region irregularities observed at Millstone Hill

Phase velocity observations at E-region heights made with the Millstone Hill 440 MHz radar find no evidence of an ion acoustic limiting speed for phase speeds observed near 0° magnetic aspect angle. Under most circumstances the phase speed increases steadily with increasing backscattered power amplitude. For a 34cm volume backscatter cross-section, σ_v, less than ∼5 × 10⁻¹³ m⁻¹, the phase speed is at or below the usual ion acoustic speed in the E-region (350m/s), and increases only slowly with the observed backscattered power amplitude (∼50 m/s per 10dB). At higher power levels, the phase speed exceeds 350 m/s, reaching values in excess of 750 m/s at times, and increases more rapidly with backscattered power (∼200 m/s per 10dB). Phase velocity/time maps observed over a 3° span of latitude suggest that many features of the phase speeds observed are directly related to changes in the ambient convection electric field in the E-region due to changing activity conditions or the effects of superimposed magnetospheric pulsations.     

Observations of ion layer motions during the AIDA campaign

The AIDA-89 campaign has yielded the most comprehensive set of low-latitude incoherent scatter radar power profiles and derived electron concentration results ever made. These results have been used to study the time-height trajectories of 80–150 km ion layers and serve to gauge both the periodicity and variability of ion layer structure throughout the campaign. Features of the AIDA ion layer trajectories point to a dynamics ‘zoo’ of processes ranging from multiday-period waves, tides and acoustic-gravity waves (AGWs) to geomagnetic storm effects and evidence of coupled neutral sodium and ion layer/plasma processes. The semidiurnal and diurnal tides are evidenced in the almost always present layers, the Tidal Ion Layers (TILs), which are identified by their regular and periodic trajectories that also display regions of variable mixing or confluence of the various tides. The TILs are contrasted with the truly sporadic layers that include sporadic E and sporadic intermediate layers. The sporadic layers may be formed due to interaction of the tidal wind system with AGWs. The formation process may involve horizontal as well as vertical ion convergence mechanisms and/or various non-linear effects. Limits to the study derive from volume undersampling due to use of the single radar beam.     

Sudden sodium layers in polar latitudes

We discuss in this paper sudden sodium layers (SSLs), which we observe with a sodium lidar instrument at Andenes, Norway (69°N). We speak of a SSL if, in a narrow altitude range (typically less than 2km), the Na density increases over the normal Na density by a factor of at least 2 within 5 min. Between December 1985 and November 1987, we have observed 42 such layers in 378 h of lidar measurements. This number increases to 75 if we only require an increase of a factor of 1.5 within 8 min. At our observation site, SSLs have the following properties: (a) they develop between 90 and 110 km altitude, (b) they develop between 20 and 02 LT, (c) their appearance shows a strong, positive correlation with that of ƒ-type E_s layers, and (d) their appearance does not show a strong correlation with either riometer absorption or meteor showers. We discuss a number of potential processes for SSL formation. SSLs above 100km can be formed in ƒ-type E_s layers by the conversion of Na⁺ ions into neutral Na. The development of SSLs below 95 km requires the presence of an additional reservoir of Na, such as Na-bearing molecules, ions, and/or ‘smoke’ particles. We also evaluate the proposal that SSLs are the outcome of single meteoroids entering the upper atmosphere, a proposition for which we find little observational support.     

Evidence for non-Maxwellian ion velocity distributions in the F-region

A study has been made of data taken with EISCAT using the Common Program CP-3-C (F-region meridian scan) which shows that regions of enhanced ion temperature (in excess of 3000K at all three EISCAT stations) are found on most days when K_p exceeds 2 or 3, usually accompanied by ion drift velocities of more than 1 km s⁻¹. These periods are often accompanied by anisotropy of the ion temperature and abnormally low apparent electron temperature, consistent with the presence of a non-Maxwellian ion velocity distribution such as would result from large but not exceptional ion drifts. Data for a selected period have been fitted using theoretical ion velocity distributions based on the relaxation collision model and assuming that the ion composition is 100% O⁺. The results confirm the presence of non-Maxwellian distributions, but a detailed comparison with theory reveals some discrepancies, indicating that the analysis may need to be extended to include effects due to, for example, molecular ions and instabilities.     

The solar eclipse of 26 February 1979: analysis and modelling of primary pair production,electron density and ionic composition

The solar eclipse of 26 February 1979 was observed from Red Lake, Canada, (52 °N, 91 °W) where totality occurred at about 1053 local time. Several research groups and government agencies participated in an extensive ground- and rocket-based observational program directed at the middle atmosphere. At the time of the eclipse, an extensive geomagnetic storm was in progress and the middle atmosphere was undergoing temperature and circulation changes associated with a stratospheric warming. Concurrent observations of atmospheric constituents, solar radiation, electron flux and other middle atmosphere parameters were obtained as inputs for a D-region predictive chemical computer code, DAIRCHEM, tailored to eclipse conditions. Ion pair production rates were computed by an E-region infrared radiance model and were used as necessary source function input values for DAIRCHEM computations. The computations yielded predictions of electron and total positive ion densities about totality. The positive ion measurements of a supersonic Gerdien condenser and a subsonic blunt probe during the eclipse were in agreement with the model computations and provided normalizing summations of total positive ions for the interpretation of mass spectrometer measurements. The chemical computer code identified principal routes for increase and removal of key species such as O₂⁺, NO⁺, hydrated clusters and negative ions. The dominant precursor ion for pair production hydrates was O₂⁺ and the chemistry was characteristic of the disturbed D-region.     

EISCAT observations of ion composition and temperature anisotropy in the high-latitude F-region

The papers by Winseret al. [(1990) J. atmos. terr. Phys.52, 501] and Häggström and Collis [(1990) J. atmos. terr. Phys.52, 519] used plasma flows and ion temperatures, as measured by the EISCAT tristatic incoherent scatter radar, to investigate changes in the ion composition of the ionospheric F-layer at high latitudes, in response to increases in the speed of plasma convection. These studies reported that the ion composition rapidly changed from mainly O⁺ to almost completely (>90%) molecular ions, following rapid increases in ion drift speed by >1 km s⁻¹. These changes appeared inconsisent with theoretical considerations of the ion chemistry, which could not account for the large fractions of molecular ions inferred from the obsevations. In this paper, we discuss two causes of this discrepancy. First, we reevaluate the theoretical calculations for chemical equilibrium and show that, if we correct the derived temperatures for the effect of the molecular ions, and if we employ more realistic dependences of the reaction rates on the ion temperature, the composition changes derived for the faster convection speeds can be explained. For the Winser et al. observations with the radar beam at an aspect angle of ϕ = 54.7° to the geomagnetic field, we now compute a change to 89% molecular ions in < 2 min, in response to the 3 km s⁻¹ drift. This is broadly consistent with the observations. But for the two cases considered by Häggström and Collis, looking along the field line (ϕ = 0°), we compute the proportion of molecular ions to be only 4 and 16% for the observed plasma drifts of 1.2 and 1.6 km s⁻¹, respectively. These computed proportions are much smaller than those derived experimentally (70 and 90%). We attribute the differences to the effects of non-Maxwellian, anisotropic ion velocity distribution functions. We also discuss the effect of ion composition changes on the various radar observations that report anisotropies of ion temperature.     

The E-region Rocket/Radar Instability Study (ERRRIS): scientific objectives and campaign overview

The E-region Rocket/Radar Instability Study (Project ERRRIS) investigated in detail the plasma instabilities in the low altitude (E-region) auroral ionosphere and the sources of free energy that drive these waves. Three independent sets of experiments were launched on NASA sounding rockets from Esrange, Sweden, in 1988 and 1989, attaining apogees of 124, 129 and 176km. The lower apogee rockets were flown into the unstable auroral electrojet and encountered intense two-stream waves driven by d.c. electric fields that ranged from 35 to 115 mV/m. The higher apogee rocket returned fields and particle data from an active auroral arc, yet observed a remarkably quiescent electrojet region as the weak d.c. electric fields (~ 10–15 mV/m) there were below the threshold required to excite two-stream waves. The rocket instrumentation included electric field instruments (d.c. and wave), plasma density fluctuation (δn/n) receivers, d.c. fluxgate magnetometers, energetic particle detectors (ions and electrons), ion drift meters, and swept Langmuir probes to determine absolute plasma density and temperature. The wave experiments included spatially separated sensors to provide wave vector and phase velocity information. All three rockets were flown in conjunction with radar backscatter measurements taken by the 50MHz CUPRI system, which was the primary tool used to determine the launch conditions. Two of the rockets were flown in conjunction with plasma drift, density, and temperature measurements taken by the EISCAT incoherent scattar radar. The STARE radar also made measurements during this campaign. This paper describes the scientific objectives of these rocket/radar experiments, provides a summary of the geophysical conditions during each launch, and gives an overview of the principal rocket and radar observations.