Equatorial spread-F and neutral atmospheric turbulence : a review and a comparative anatomy

A number of satellite and rocket plasma density spectra obtained during equatorial spread-F conditions are presented and discussed in the light of similar measurements in the neutral atmosphere. We discuss this comparison in some detail and find both distinct similarities and subtle differences. The horizontal spectral measurements show a peak at an outer scale quite similar to the scale of the undulations caused by gravity wave interactions with the ionosphere. This feature is similar to a buoyancy subrange but it is easy to show that the amplitudes of the plasma fluctuations are too large to be directly driven by the neutral atmosphere. At intermediate scales the plasma fluctuations have a one-dimensional horizontal spectrum with a power law well described by a ($\text{−5}\text{3}$) slope. Once again it can be shown that the neutral fluid cannot be similarly structured at 400 km altitude due to the high viscosity coefficient. A plasma cascade process seems to be operating but it is not at all clear how the spectrum is formed. Furthermore, vertical power spectra seem to run the gamut in spectral form with slopes (n) varying in the range from −1 to −3. So the horizontal spectra are near universal in form, while the vertical spectra are quite variable.     

Equatorial spread-F: recent results and outstanding problems

Recent results and outstanding problems in the field of equatorial spread-F are reviewed. The discussion is organized about four distinct wavelength regimes: long (≳ 20 km). intermediate (20 km–100 m). transitional (100 m–10 m) and short (≲ 10 m). The intermediate waves are best understood, although the enhanced turbulent power near 1 km wavelength needs explanation, as does the saturation mechanism itself. The role of shear flow and irregularity seeding by ‘geophysical noise’ in the neutral atmosphere form important future research topics at the longest scales. Identification of drift waves in the transitional wavelength range has been a major new step in explaining the full spread-F phenomenon. This drift wave turbulence in the transition regime will be actively studied in the near future, particularly with regard to its role in determining the saturated amplitude of the intermediate waves, as well as in the anomalous diffusion of plasma at high altitudes. Although waves as short as 11 cm have been unambiguously detected and linear theories exist, the origin and amplitude spectra of these short wavelength waves inside topside plumes remain in some doubt and should also be examined in more detail.     

Large-scale TIDs and upper atmospheric gravity waves

Excitation of the guided acoustic-gravity waves in the upper thermosphere in response to enhanced auroral electrojets is calculated in the absence of dissipation under a fully ducted condition. It is shown that a model atmosphere terminated with an isothermal half-space supports a long-period, high-speed mode, which is the interface mode guided along the half-space termination of the atmosphere. The dispersion properties and the vertical distributions of the kinetic energy density of this mode are similar to those of the so called ‘gravity pseudomode’. The excitation of this mode is computed to show how the wave generation depends on the source mechanism (the Lorentz force and joule heating) and on the source altitude. Joule heating can generate the waves with appreciable amplitudes. On the other hand, the Lorentz force prevailing in the lower region cannot excite the waves with any observable amplitudes. The waves are intensified with increasing the heat source altitude. The gross features of the calculated waves indicate that the ducted thermospheric gravity waves are capable of producing observable thermospheric waves. It is therefore suggested that further examination of the excitation of the ducted acoustic-gravity waves undergoing partial reflections due to viscosity and thermal conduction should be useful for the theory of large-scale travelling ionospheric disturbances.     

The role of acoustic gravity waves in the generation of spread-F and ionospheric scintillation

In the aggregate, acoustic gravity waves in the F-region constitute a spectrum of geophysical noise extending from the frequencies involved in diurnal variations up to the Brunt-Väisälä buoyancy frequency. They drive a roughly uniform power spectrum of travelling ionospheric disturbances (TIDs) with vertical scales of the order of the atmospheric scale height H and with horizontal scales extending from the radius of the Earth down to H. It has been known since the 1950s that this permits multiple normals onto the F-region from an ionosonde, thereby creating the multiple-trace type of spread F on ionograms. At shorter scales the spectrum of TIDs decreases in strength and, below the mean free path of the neutral atmosphere, creates a spectrum of plasma turbulence aligned along the Earth's magnetic field. Progressively shorter scales are responsible for phase scintillation, for amplitude scintillation and for blur-type spread F on ionograms. A weak extension of the spectrum to scales less than the ion gyroradius is responsible for spread F and transequatorial propagation in the VHF band. Under evening conditions in equatorial regions a band of TIDs with wavelengths of the order of 600 km can, at times, have a phase velocity that matches the drift velocity of the plasma (Röttger 1978). This band of TIDs is then amplified until it breaks (Klostermeyer 1978). The associated explosive increase in plasma turbulence creates the plume phenomenon discovered by Woodmn and La Hoz (1976).     

Three-dimensional ionospheric currents and fields generated by the atmospheric global circuit current

Three-dimensional ionospheric currents and fields generated by atmospheric global circuit currents, using the distribution of air-Earth currents as a lower boundary conditions of the ionosphere, have been studied. The air-Earth currents are obtained taking geomagnetic and orographic effects into account, under the assumption of an ionosphere with infinite conductivity. Three dimensional ionospheric currents due to thunderstorm sources are calculated, considering the conductivity distribution in the ionosphere and the configuration of the magnetic field. The calculated potential difference in the ionosphere is 55 V and according to our model the horizontal electric field is too weak to affect the ionosphere and magnetosphere significantly. Horizontal currents are not distributed uniformly, but preferably in the day-side hemisphere and especially in the equatorial region, and vertical currents and fields do not simply decrease with altitude near the equator because of anisotropy and nonuniformity in the conductivity.     

Measurement of vertical phase and group velocities of atmospheric gravity waves by VHF radar

A systematic method of deriving from MST radar data the group velocity and phase velocity of the atmospheric wave along the radar beam direction is proposed and verified by a series of numerical simulations. We apply the method to two data sets measured by Chung-Li radar under different background wind conditions. It is found that the vertical group velocity and phase velocity are mostly in the opposite direction when the background wind is weak. The energy source of downgoing wave packets was evidently related to the instability in the upper height range (10.5–11.7 km) where strong wind shear existed. When the background wind and wind shear are stronger, the vertical group and phase velocities may propagate in the same direction. We also found from numerical simulation and data analysis that the wave packet of gravity waves following power law spectrum are short-lived. A by-product of the group velocity measurement is that the horizontal wavelength may also be deduced from a vertical radar beam measurement from the dispersion relation if it is valid.     

Recent measurements of middle atmospheric electric fields and related parameters

In 1989, two series of rocket measurements were carried out to investigate middle atmosphere electric fields. The measurements were taken both in the Northern Hemisphere on Heiss Island (80°37′N and 58°03′E) and in the Southern Hemisphere in the Indian Ocean (40–60°S and ~45°E) on board the research vessel ‘Akademik Shirshov’. Along with the vertical electric fields, aerosol content and positive ion density were also measured. Some of the rocket launches were made during the extremely strong solar proton events (SPE) of October 1989. The experiments showed the strong variability of the electric fields in the middle atmosphere at polar and high middle latitudes. In all the measurements the maximum of the vertical electric field height profile in the lower mesosphere was observed to be more than ~ 1 V/m. The electric field strength and the field direction at maximum varied considerably among the launches. A maximum value of + 12 V/m was detected at a height of about 58 km at 58°30′S on 21 October 1989 during the SPE. The simultaneous measurements of the electric field strength, positive ion density and aerosols point out both an ion -aerosol interaction and a connection between the mesospheric electric fields and aerosol content.     

A test of the Hines dispersion equation for atmospheric gravity waves

Simultaneous observations of an ionospheric wave by two incoherent scatter facilities and three Faraday-rotation polarimeters have provided measurements of the frequency, vertical wavelength, horizontal wavelength and direction of propagation of the wave. These measured values confirm the Hines dispersion equation for atmospheric gravity waves.     

The generation and propagation of atmospheric gravity waves from activity in the auroral electrojet

EISCAT measurements of the electric field in the auroral electrojet are compared with the signature of TIDs propagating equatorward as observed by an HF-Doppler network. At night-time the onset of auroral activity is usually followed by the arrival of a TID at lower latitude. Cross-correlation of the time variations of the electric field measured by EISCAT and the frequency offset recorded by the HF-Doppler system confirms a relationship between the auroral activity and the gravity wave, indicating both the travel time and the periodicity of the wave. The relationship is especially close under quiet conditions when the cross-correlation coefficient is typically 60%, significant at 0.1%. When the observed electric field is used as input to a thermosphere-ionosphere coupled global model it predicts the time signature of the observed HF-Doppler variation reasonably well but seriously underestimates the amplitude of the disturbance. Examination of this discrepancy may lead to a better understanding of the mechanisms involved in the generation and propagation of atmospheric gravity waves.     

Wind and temperature effects on F-region medium-scale gravity waves estimated using a multi-layer atmospheric model

Diurnal variations in the propagation direction of atmospheric gravity waves, and the travelling ionospheric disturbances to which they give rise, have been observed in many experimental observations and several modelling studies have demonstrated that this is primarily due to the corresponding diurnal rotation in the direction of the thermospheric wind. Other variations have been attributed to seasonal or other effects, but the effects of variations in the thermospheric temperature have not previously been analysed in detail. We present results from a study of the propagation of gravity waves through a layered atmosphere in which the thermospheric wind and temperature are derived from a three-dimensional time-dependent model. The analysis has been carried out for a range of wave speeds and periods, and for a range of times, seasons and propagation azimuths. Results suggest that a significant diurnal variation in the transmission coefficient for waves propagating through the thermosphere exists with seasonally dependent maxima. Transmission increases for increasing wave period up to about 50 min, after which it remains approximately constant. Maximum transmission occurs for wave phase speeds around 200–250 m/s and falls to zero for speeds less than about 100 m/s. An exception to this rule occurs for waves with periods less than 40 min and speeds less than 50 m/s for which significant transmission appears to be theoretically possible.     

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.     

A Van de Graaf source mechanism for middle atmospheric vertical electric fields

It is proposed that meteoric and other debris descending through the mesosphere constitute a natural Van de Graaf generator for vertical electric fields within the mesosphere. Dust and aerosol particles falling from above 85 km are charged negatively in the upper D-region. Charge is lost in the region below 70 km. This net charge transport creates a vertical polarization electric field. Calculated fields are in the range of 10 mV/m for the average input of meteoric debris. Observed vertical electric fields are confined to a few occasions when large fields of the order of 4 V/m are observed to maximize at 65 km. Calculated fields from this model also maximize at this altitude, but a special event with increased dust density or another mechanism to increase relative vertical velocity is required to explain the large fields. Such large values are the exception rather than the rule for D-region vertical electric fields.     

Atmospheric gravity waves and sporadic-E

Observations of atmospheric gravity waves have identified several occasions when a wave-group propagating upwards can be associated first with a sporadic-E layer whose total electron content varies periodically, and later with a disturbance in the F-region varying with the same period. This paper reports four examples of such behaviour.     

Equatorial ionospheric electric fields during magnetospheric disturbances: local time/longitude dependences from recent EITS campaigns

Data sets collected during a few coordinated Equatorial Ionosphere-Thermosphere System (EITS) observational campaign periods, mainly from the Brazilian and Asian longitude sectors are analysed in this paper. Ionosonde magnetometer and Ionospheric Electron Content (IEC) data from the EITS-1 and -2 campaigns (during March and December 1991) are complemented by interplanetary magnetic field and some ground based data sets from other campaigns. The analysis focuses on the response of the equatorial ionospheric heights and ionization anomaly to disturbance electric fields, identified as a direct penetration electric field associated with IMF B_z changes and development of the ring current (especially the asymmetric component), and that produced by a disturbance zonal neutral wind. New evidence on the local time and longitudinal dependences of these electric fields constitute the main results of this paper. Especially, a large eastward electric field (associated with the asymmetric ring current) in the dusk-dawn sector causes significant expansion of the EIA in this sector, and amplification of the evening prereversal uplift of the F-layer over Brazil. Significant inhibition of the evening prereversal electric field enhancement seems to be produced by the disturbance zonal wind associated with the magnetic disturbances prevailing several hours earlier. Some tentative evidence on the Brazilian dusk sector disturbance field being larger than that of the Asian dusk sector support the existence of a longitude asymmetry in the intensity of the disturbance electric field.     

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.     

A theoretical investigation of sources of large and medium scale atmospheric gravity waves in the auroral oval

Many papers have been published to devise models to describe the sources of large and medium scale atmospheric gravity waves (AGWs) in the auroral oval ionosphere. One of the models proposed by Chimonas and Hines [(1970) Planet. Space Sci.18, 565] calculated the relative importance of Lorentz force and Joule heating as sources of AGWs in the auroral regions based on certain assumptions. In this paper, we develop a general theory to describe the behavior of the AGW source terms. It has been found that the source terms which generate AGWs are closely related to the velocities and frequencies of AGWs, and that the Lorentz force is dominant in generating the vertical velocity perturbation of large scale AGWs. The formulas which determine the source term contributions are derived. This relationship gives us the possibility to predict what kind of AGW will be generated by observing the source terms, or conversely perhaps to deduce some of the source characteristics by measuring properties of the traveling ionospheric disturbances (TIDs).     

Evaluation of the vertical flux of energy into the thermosphere from medium scale gravity waves generated by the jet stream

Since the paper of Hines (1960), the heating of the thermosphere by gravity waves is well known. Until recently there were no statistical evaluations of their energy. The present paper is a study of this subject. From experimental data of a Faraday rotation experiment and with the help of a measurement theory and a propagation model in a dissipative atmosphere, we evaluated the mean energy flux for the medium scale gravity waves detected during July 1974. We found 0.11 erg cm⁻² s⁻¹ for waves with horizontal phase velocity between 100 and 180 m s⁻¹. We also evaluated the heights of energy deposition which are from 120 to 180 km. This study shows that the heating by gravity waves may be important for the thermospheric equilibrium.     

The generation and propagation of atmospheric gravity waves observed during the Worldwide Atmospheric Gravity-wave Study (WAGS)

During the Worldwide Atmospheric Gravity-wave Study (WAGS) in October 1985, the EISCAT incoherent scatter radar was used to observe the generation of atmospheric gravity waves in the auroral zone in conjunction with a network of magnetometers and riometers. At the same time a chain of five ionosondes, an HF-Doppler system, a meteor radar and a radio telescope array were used to monitor any waves propagating southwards over the U.K.The EISCAT measurements indicated that in the evening sector both Joule heating and Lorentz forcing were sufficiently strong to generate waves, and both frequently showed an intrinsic periodicity caused by periodic variation in the magnetospheric electric field.Two occasions have been examined in detail where the onset of a source with intrinsic periodicity was followed by a propagating wave of the same period which was detected about an hour later, travelling southwards at speeds of over 300 m s⁻¹, by the ionosondes and the HF-Doppler radar. In both cases the delay in arrival was consistent with the observed velocity, which suggests a direct relationship between a source in the auroral zone and a wave observed at mid-latitude.     

Generation of gravity waves by inhomogeneous heating of the atmosphere

Local variation of atmospheric heating which might occur in inhomogeneities of various constituents such as ozone or molecular oxygen may generate gravity waves. These perturbations are induced by the terminator crossing constituent inhomogeneities of short lifetime. The quasi-point heating model developed here shows that the largest amplitude must appear vertically above the source, where the perturbation frequency is close to the Vaisala-Brunt frequency. Numerical calculations not band limited in the frequency suggest several characteristics of the perturbation.     

Dynamic coupling between the lower and upper atmosphere by tides and gravity waves

Results of a General Circulation Model simulation of the dynamics of the middle atmosphere are shown focusing our attention to the tidal wave mean flow interaction and propagation of migrating diurnal and semidiurnal tides in the model. It is shown that migrating tidal waves are well simulated and the amplitude growth with height is effectively suppressed by the convective adjustment in the model. It is also shown that the dissipating solar diurnal tide plays an important role in inducing mean zonal winds in the low latitude region of the lower thermosphere. The behavior of non-migrating diurnal tides is also analyzed to show that non-migrating diurnal tides have significant amplitudes in the lower thermosphere. It is suggested that the non-migrating diurnal tide, which propagates against background mean zonal winds, has the possibility to propagate into the middle to high latitude region due to the Doppler effect.