Modeling the low latitude ionospheric total electron content

The undisturbed ambient total electron content of the ionosphere in the equatorial region exhibits two characteristic features:

• 1.(i) a longitudinal behavior of the post-sunset variation of the ionization near the crests of the equatorial anomaly
• 2.(ii) an enhancement at lower latitudes following the post-sunset decay. During high solar activity periods the southern crest of the equatorial anomaly in the African longitude sector is characterized by a post-sunset maximum often exceeding the afternoon maximum. In the Indian and other longitude zones, the post-sunset peak is not so prominent. Instead, a ledge is obtained in the corresponding local time period. At lower magnetic latitudes, the ionization decays very rapidly around sunset, but an enhancement lasting 2–4 h is observed afterwards.

Numerical solution of the plasma continuity equation, including the effects of ionization production by solar ultraviolet radiation, loss through charge exchange and transport by diffusion, electrodynamic drift and neutral wind, has been used to investigate the above two features. It is found that the pre-reversal peak of the E × B drift at the magnetic equator around sunset is the dominant mechanism responsible for the post-sunset behavior near the crests of the equatorial anomaly. The zonal wind causes an asymmetry of the total content in the northern and southern hemispheres. In African longitudes, where the magnetic declination is about 20°W, the southern crest is more developed at the expense of the northern counterpart. The north-south asymmetry is practically absent in the Asian sector, with its low (< 5°) declination angle. In the Pacific area, an easterly declination (about 9°E) results in a higher post-sunset ionization at the northern crest, although the asymmetry is less pronounced than that in the African zone. The night-time enhancement at lower latitudes has been found to be controlled by the post-sunset increase in the vertical drift, possibly also modulated by the neutral wind.     

Tropospheric gravity waves and their possible association with medium-scale travelling ionospheric disturbances

Analysis of pressure fluctuations observed over a period of several days using an array of microbarographs has shown the existence of long trains of gravity waves with two or more waves often present simultaneously. Meteorological data from radiosonde ascents indicates that many of the waves have a velocity which matches that of the background wind at some level within the troposphere. Generally this height corresponds to that of a frontal zone marking the transition between air masses and it is suggested that the waves may have been generated by shear flow instability within the frontal layer. Theoretical considerations, based on a three-layer model troposphere, show that some of the observed waves could have been ducted in or near the frontal zone. Some evidence is found to indicate that a non-linear wave-wave interaction between pairs of waves occurring simultaneously in the frontal zone could yield secondary waves with the characteristics of the gravity waves which had been observed in the thermosphere at appropriate times and whose group paths were traced to source regions in the troposphere in the general vicinity of the microbarograph array.     

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).     

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 case study of a high latitude ionospheric electron density depletion

Incoherent scatter measurements inside and outside an ionospheric electron density depletion are described. The density depletion is probably caused by an enhancement of NO⁺-ions and subsequent dissociative recombination. The NO⁺-ions are increased because high electric fields present at the geographical location of the density depletion speed up the reaction O⁺ + N₂ → NO⁺ + N. The electron as well as the ion temperature within the density depletion are strongly enhanced, the latter due to Joule heating, also caused by the electric field.     

A study of ionospheric electron density deviations during two great storms

A large foF2 data base has been collected for the two great storms occurring in March and October 1989. Plots of foF2 deviations versus time and latitude show the large depressions cannot be caused by a composition change alone; the storm induced dynamo electric field must also play an important role, especially for low latitudes. The observed ionospheric response is hemispherically asymmetric, with the autumnal hemisphere suffering a longer lasting, more latitudinally extended, and deeper depression.     

GPS satellite measurements: ionospheric slab thickness and total electron content

A preliminary analysis was made of ionospheric slab thickness, τ, and total electron content, TEC, for southern Australia using GPS satellite measurements. It was found that at mid-latitudes τ has similar overall diurnal, seasonal and latitudinal variations in the southern hemisphere as in the northern hemisphere. However, there are appreciable differences between τ in the two hemispheres which would justify appropriate modifications to ionospheric models based on northern hemisphere data before being applied confidently to the southern hemisphere. The usefulness of GPS satellites together with ionosondes over a spread of latitudes was demonstrated in determining long-term variations of TEC and τ over a large area. It was concluded that as few as four GPS receivers could provide TEC for the whole of Australia in real-time, though approximately six receivers in convenient locations would be required in practice.     

Production of electromagnetic field disturbances due to the interaction between acoustic gravity waves and the ionospheric plasma

The problem of electromagnetic field disturbances produced by the interaction between winds of acoustic gravity waves (AGW) origin and the ionospheric plasma has been considered. It is shown that, when not allowing the electrostatic approach, electromagnetic field disturbances represent shear Alfvén and compressional modes modified by ionospheric Pedersen and Hall conductivities. It is further shown that the quasielectrostatic Alfvén type disturbances give the main contribution to electric field perturbations. Magnetic field perturbations due to Alfvén and compressional modes have the same order of magnitude. Two numerical models for simulation of the problem under consideration have been developed. The first model is intended for the simulation of Alfvén type disturbance production and transmission into the magnetosphere, taking into account the dipole geometry of the geomagnetic field, but a mutual transformation of Alfvén and compressional modes is ignored. The second model is constructed for the simulation of both electromagnetic field disturbance production and their mutual transformation in the ionosphere. The results of numerical simulations with these models show that there is an opportunity for AGW activity monitoring in the lower thermosphere by ground-and satellite-based recordings of magnetic and electric field variations.     

On medium scale travelling ionospheric disturbances as seen through total electron content power spectra

Modulation phase of 140 MHz with respect to 360 MHz signals from ATS 6 satellite recorded at Slough (51.5°N, 0.6°W geographic latitude and longitude; 54.3° geomagnetic latitude) at one-minute interval are power spectral analysed to derive dominant periodicities corresponding to Travelling Ionospheric Disturbances. From the significant peaks in the spectra, an occurrence peak at periods between 10 and 15 min and a secondary peak at 60–65 min are seen. From cross-spectral analysis of the same records from three stations separated by a few hundred kilometers for a short period, the speed and azimuth of the propagating disturbances are determined. During the day-time, most of the waves in the period range 30–100 min are seen to propagate at azimuths of 90–160°. At night-time they propagate poleward. Theoretical computations of the azimuth response of TEC to typical gravity waves, including the effect of neutral winds, show that the observed azimuths of propagation are in reasonably good agreement with theoretical predictions.     

Temporal variations of ionospheric waves in the D- and F-regions

It is shown that mesoscale ionospheric wave disturbances in the D- and F-regions regularly occur at all times of day, night and season (characteristic periods ~100, 24, 12.6 min) and are a characteristic property of dynamic processes in the ionosphere.     

A quartic dispersion equation for internal gravity waves in the thermosphere

A new quartic dispersion equation in the square of the complex vertical wave number is derived by employing the ‘shallow atmosphere’ approximation and an ion drag approximation. These approximations allow the coefficients of the quartic equation to be given in terms of the corresponding cubic equation, which neglects the Coriolis force and the zonal ion drag component, but modified to take into account these neglected effects. Coupling between the extraordinary viscosity wave mode and the other three wave modes is highlighted and numerical solutions are compared for this quartic equation, an exact eighth order equation and the cubic equation. For the first time the validity of using the ‘shallow atmosphere’ approximation to describe internal gravity wave motions is demonstrated.     

Ground and in situ excitation of waves in the ionospheric plasma

A review is presented of seven papers the contents of which range from ULF to VLF wave excitation in the ionospheric plasma by ground-based radio wave and acoustic wave sources, to in situ plasma wave excitation by satellite- and rocket-borne radio transmitters.     

Climatology of gravity waves in the middle atmosphere

An attempt is made at the statistical analysis of small-scale disturbances in the stratosphere and mesosphere with the aid of meteorological rocket observations at many stations from 77°N to 8°S for several years.By applying a high-pass filter to daily rocket data in the height range 20–65 km, wind and temperature fluctuations with characteristic vertical scales close to or less than 10 km are obtained, which are considered to be due to internal gravity waves. Results are expressed in terms of parameters which tend to emphasize smallscale vertical fluctuations and which should provide qualitative measures of gravity wave activity.It is found that the gravity wave activity shows a notable annual cycle in higher latitudes with the maximum in wintertime, while it shows a semiannual cycle in lower latitudes with the maxima around equinoxes. It is also found from the standard deviation around the monthly mean that the temporal variability of gravity waves is very large.     

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.     

Optical observations of gravity waves in the auroral zone

Night-time observations of O(¹D) λ630 nm and O(¹S) λ558 nm thermospheric emissions were made at Mawson, Antarctica (67.6°S, 62.9°E) from 1982 to 1989, using a three-field photometer. Crossspectral analysis of the data was used to extract frequencies and horizontal trace velocities of periodic structures. Structures in the λ630 nm emission were characteristic of large-scale waves, and those in the λ558 nm emission were characteristic of medium-scale waves. The results showed distinct polarisation of the propagation azimuths; waves in the λ630 nm emission propagated approximately northwestward throughout the 8 yr period, whilst propagation azimuths of waves in the λ558 nm emission appeared to be solar-cycle-dependent. It is suggested that waves observed in the λ630 nm emission were of predominantly auroral electrojet origin, whilst those observed in the λ558 nm emission were of both auroral and tropospheric origin.     

Solar cycle variation in the total electron content at Sagamore Hill

The TEC data obtained at Sagamore Hill observatory by using ATS-3 beacon signal during the period from November 1967 to December 1976 have been used to analyze the solar cycle variations of total electron content at invariant latitude 54°. The correlation coefficients between TEC and sunspot number were found to be largest if 12 month running mean values were used, and to be smallest if monthly mean values were used. By the method of linear regression analysis, the contour charts for real diurnal and seasonal variations of TEC at given sunspot number were constructed and described. The diurnal variation of TEC was represented by the sum of its diurnal mean and first three harmonic components. The solar cycle variations of these components were also given.     

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.     

Mechanisms for observed total electron content pulsations at mid latitudes

Total electron content variations in the Pc3–Pc4 range of frequencies of the order of 4 parts in 10⁴ have been reported in apparent correlation with simultaneous ground based magnetic pulsation observations. By means of a term-by-term analysis of the continuity equation for electrons, the plausibility of various mechanisms is investigated. The most likely explanation is in terms of localized increases in the electron density at F-region heights caused by the field-aligned (compressional) component of the pulsation magnetic field. The analysis predicts a tendency for the amplitude of the TEC pulsations to vary in antiphase with ground-based measurements of the north-south component of the pulsation field.     

Nonlinear disturbances in the ionosphere due to acoustic gravity waves

The influence of the higher harmonics of an internal gravity wave on the formation of nonlinear quasi-periodic disturbances in the F-region of the Earth's ionosphere is considered. It is shown that the Boussinesq approximation cannot be used in describing a plane nonlinear gravity wave as nonlinearities associated with the compressibility of the atmosphere have to be taken into account.