Numerical simulation of the penetration and reflection of a whistler beam incident on the lower ionosphere at very low latitude

By the full-wave algorithm with Fourier synthesis, 3-D propagation of a whistler beam incident on the pre-dawn lower ionosphere at very low latitude is numerically investigated. Processes of transmission, reflection, and coupling with the Earth-ionosphere waveguide are discussed via the wave energy and polarisation distributions and their dependence on the wave parameters and the ionospheric profile (such as the E_s-layer). It is shown that the dominant wave above 90 km altitude has the propagation characteristics of the magneto-ionic whistler mode, and absorption, spreading, reflection and mode conversion mainly occur at, and are greatly affected by, the bottom of the ionosphere. It is found that the transmitted energy density along the Earth's surface is reduced by 20 dB or more. Beam transmission loss varies asymmetrically with the incident angle, but changes little with the frequency. In the region 150 km (for 5 kHz) away from the ‘exit area’ where whistlers emerge, the bearing measurements using ground-based VLF direction-finders may be in error because direction-finding algorithms assume plane wave propagation. Only a small portion (about −25 dB at 5 kHz) of the incident energy is reflected up to an altitude of 150 km, and major reflection takes place in a small range of altitude at the bottom of the ionosphere with little spreading and lateral shift with respect to the incident beam. Reflection is enhanced considerably at lower frequency. Our results also suggest that an E_s-layer or an ionospheric gradient refracting waves to higher latitudes would be favorable factors for multi-hop echoes to be received on the ground.     

Penetration of hydromagnetic waves of cylindrical symmetry through the high latitude ionosphere

The problem is formulated, and boundary conditions are developed, in order to solve numerically the equations for the penetration of hydromagnetic waves of horizontal cylindrical symmetry through a stratified high latitude ionosphere and atmosphere due to sources above the ionosphere. There are two orthogonal polarization modes. The numerical results are shown for various cases. It is found that near the ground the horizontal components of both the electric and magnetic fields along the circumference direction are much smaller than the radial components.     

Comparisons between EISCAT observations and model calculations of the high latitude ionosphere

Calculations using a numerical model of the convection dominated high latitude ionosphere are compared with observations made by EISCAT as part of the UK-POLAR Special Programme. The data used were for 24–25 October 1984, which was characterized by an unusually steady IMF, with Bz < 0 and By > 0; in the calculations it was assumed that a steady IMF implies steady convection conditions. Using the electric field models of Heppner and Maynard (1983) appropriate to By > 0 and precipitation data taken from Spiroet al. (1982), we calculated the velocities and electron densities appropriate to the EISCAT observations. Many of the general features of the velocity data were reproduced by the model. In particular, the phasing of the change from eastward to westward flow in the vicinity of the Harang discontinuity, flows near the dayside throat and a region of slow flow at higher latitudes near dusk were well reproduced. In the afternoon sector modelled velocity values were significantly less than those observed. Electron density calculations showed good agreement with EISCAT observations near the F-peak, but compared poorly with observations near 211 km. In both cases, the greatest disagreement occurred in the early part of the observations, where the convection pattern was poorly known and showed some evidence of long term temporal change. Possible causes for the disagreement between observations and calculations are discussed and shown to raise interesting and, as yet, unresolved questions concerning the interpretation of the data. For the data set used, the late afternoon dip in electron density observed near the F-peak and interpreted as the signature of the mid-latitude trough is well reproduced by the calculations. Calculations indicate that it does not arise from long residence times of plasma on the nightside, but is the signature of a gap between two major ionization sources, viz. photoionization and particle precipitation.     

A new aspect of daily variations of the geomagnetic field at low latitude

Data from the unique network of low latitude geomagnetic observatories in India extending from the dip equator to the northern focus of the Sq current system have shown a new type of Sq current distribution different from those associated with the normal or the counter electrojet currents. On 3 December 1985 both the horizontal as well as the vertical components of the geomagnetic field at Annamalainagar showed maximum values around the midday hours. The abnormal feature described seems to be rather a rare phenomenon. The solar daily range of H field is found to be fairly constant from the dip equator up to about 12° dip latitude, suggesting the complete absence of the equatorial enhancement of ΔH, typical of the equatorial electrojet. The cancellation of the equatorial electrojet is suggested to be caused by a westward flowing current system much wider than the conventional equatorial electrojet. This additional current system could be due to the excitation of certain tidal modes at low latitudes on such abnormal days.     

Ray-tracing the paths of very low latitude whistler-mode signals

VLF whistler-mode signals with very low group delays (75–160 ms) received at night in Dunedin, N.Z., from the 23.4 kHz MSK transmissions of NPM, Hawaii (21.5°N, 158°W), are explained by ray-tracing along unducted paths. The typical vertical and horizontal electron density gradients of the night equatorial ionosphere are found to be sufficient to explain not only the typical group delays but also their decrease during the night and the typical frequency shifts observed on these signals. An important feature appears to be the decreasing starting and finishing latitudes (and the decreasing maximum height of the path) during the course of the night. The amplitude of the signals in relation to the expected collisional absorption in the ionosphere is discussed. A simple but quite accurate analytical expression suitable for ray-tracing is derived for the night electron density in the height range 170–1400 km, based on non-isothermal diffusive equilibrium and O⁺/O friction.     

Production and control of ion-cyclotron instabilities in the high latitude ionosphere by high power radio waves

The possible generation and suppression of ion-cyclotron waves in a collisional plasma by external high power electromagnetic (EM) waves with frequency close to the local upper-hybrid frequency is considered. It is shown that the ion cyclotron instability can be destabilized (stabilized) for ω₀ <ω_UH (ω₀ > ω_UH), where ω₀ is the pump frequency of the EM wave. The results are applied to naturally occurring ion-cyclotron instabilities in the high latitude ionosphere.     

A method of determining horizontal structure of the ionosphere from backscatter ionograms

We have developed a technique lor oblique backscatter sounding (OBS) ionogram inversion as a diagnostic tool for the horizontally inhomogeneous structure of the ionosphere. Input data for the method include the leading edge of a backscalter ionogram that is measured through soundings in a given direction, and the vertical electron density profile measured over the sounding station or over some other site lying in the sounding direction. The method may be useful for reconstructing the two-dimensional electron density distribution in a vertical plane aligned with the direction of sounding. The inverse problem has been solved using the Newton Konlorovich method and the Tikhonov regularization method. The algorithm we have developed was tested against model data, that is, OBS ionograms synthesized using geometrical optics calculations for different models of the inhomogeneous ionosphere. Test results demonstrate that our method converges reliably, is stable to measurement errors and provides a good accuracy of reconstruction of inhomogeneous structures with scales of 100 2000 km. This indicates that this method shows promise as an operative remote diagnostic tool for ionospheric irregularities of natural and artificial origin.     

Satellite observations of zonal electric fields near sunrise in the equatorial ionosphere

We report here on a number of examples of anomalous enhancements of eastward electric fields near sunrise in the equatorial ionospheric F-region. These examples were selected from the data base of the equatorial satellite, San Marco D (1988), which measured ionospheric electric fields during a period of solar minimum. The eastward electric fields reported correspond to vertical plasma drifts. The examples studied here are similar in signature and polarity to the pre-reversal electric field enhancements seen near sunset from ground-based radar systems. The morphology of these sunrise events, which are observed on about 14% of the morning-side satellite passes, are studied as a function of local zonal velocity, magnetic activity, geographic longitude and altitude. The nine events studied occur at locations where the zonal plasma flow is generally measured to be eastward, but reducing as a function of local time and at satellite longitudes where the magnetic declination has the opposite polarity as the declination of the sunrise terminator.     

F-region neutral winds from ionosonde measurements of hmF2 at low latitude magnetic conjugate regions

The behavior of the F2-peak height difference Δh_mF2, between low latitude magnetic conjugate points, is known to be governed by thermospheric winds blowing along the magnetic meridian. Ground based ionosonde measurements of h_mF2, at two pairs of magnetic conjugate stations, have been analysed in conjunction with the results of a realistic dynamic computer model of the tropical ionospheric F-region, to determine thermospheric wind velocities. The behavior of monthly average values of the sum, at conjugate points, of the thermospheric horizontal wind velocity component in the magnetic meridian, at low latitudes, has been inferred for months of solstice and equinox, as well as for periods of low and high solar activity.     

A study of wave packets and phase jumps in low latitude Pc3 geomagnetic pulsation

An array of four low latitude induction coil magnetometer stations has been used to study the spatial and temporal characteristics of Pc3 pulsations over a longitudinal range of 17° at L = 1.8 to 2.7 in southeast Australia. A preliminary study of individual Pc3 wave packet structure at the azimuthal stations has established the existence of phase jumps between wave packets at low latitudes, similar to those observed at synchronous orbit and at higher latitude ground stations. However, there did not appear to be any obvious pattern in phase jump occurrences between stations or signal components.     

Acceleration of electrons in the ionosphere under the action of intense radio-waves near electron cyclotron harmonics

The effect of electron acceleration by intense ionospheric plasma turbulence induced by a high-power radio wave is studied theoretically, under conditions when the turbulence frequency is nearly equal to harmonics of the electron-cyclotron frequency. The turbulence located in a thin layer at the region of the pump-wave upper-hybrid resonance gives rise to the formation of an intense electron distribution function tail at high energies. Given the resonantly absorbed fraction of the pump-wave energy and the typical plasma turbulence scale, one can estimate the turbulence energy density and calculate the modified electron distribution function.     

Solar cycle and seasonal variations of the low latitude OI 630 nm nightglow

Regular zenith measurements of the OI 630 nm nightglow emission have been carried out at Cachoeira Paulista (22.7°S, 45.0°W; geomag. 11.9°S), Brazil, since 1975. The long series of observations during the period 1975–1982, including the ascending phase of the last solar cycle, permitted studies of solar cycle effects and seasonal variations. A large intensity increase, about seven times, from low solar activity to high solar activity has been observed. Also, the seasonal-nocturnal intensity variations show large changes between years of low and high solar activity. The characteristics of the variations observed are closely related to the equatorial electric field variations, since the observation site is under the southern equatorial ionospheric anomaly crest.     

Determination of the turbulent spectrum in the ionosphere by a tomographic method

A theory of tomographic reconstruction of the statistical properties of the random turbulent ionospheric plasma is presented. Derived integral equations for the coherence functions of the measured fields allow the determination of inhomogeneous layer coordinates and the reconstruction of cross-sections of the electron density correlation functions. For statistically homogeneous layers and a transmitter on board a moving satellite with a linear receiving array on the ground, we have the possibility of determining the three-dimensional correlation function structure or its spectrum using a set of two-dimensional cross-sections. One receiver allows the reconstruction of the spectrum of the two-dimensional cross-section of the correlation function. We also consider the solution of the inverse problem for non-homogeneous fluctuations. In this case the distribution of the electron density fluctuations, its variance and the correlation coefficient, characterizing the spatial structure of fluctuations may be reconstructed by a tomographic technique. Experimental results on the identification of the layer height of the irregularities and on the spectrum of the two-dimensional cross-sections of the correlation function measurements are presented.     

Very unusual low latitude whistlers with additional traces of the Earth-ionosphere waveguide propagation effect

The present paper reports very unusual whistlers strongly influenced by the Earth-ionosphere waveguide propagation after emerging from the ionosphere, as observed simultaneosly at our two stations, Sakushima (geomag. lat. 24°) and Kagoshima (20°). These unsual whistlers are characterized by clearly exhibiting additional dispersion effects near the cut-off frequencies of the 1st and 2nd order modes of waveguide propagation and, to our knowledge, they are a new finding. All the subionospheric dispersion is deduced to occur between the ionospheric exit point and the receiver. Detailed spectral analysis, after extracting the small waveguide dispersion effect from the overall spectrum by taking the beat with the appropriate pseudo-whistler, has enabled us to determine the propagation distance of the ionospheric exit region from each station. These distances have then been used to locate the ionospheric exit region, which is found to be about 3000 km east of the stations and in the local sunrise time sector. The generation mechanism of such unusual whistlers is discussed in terms of the joint influences of the ionospheric transmission mechanism (longitudinal gradient of the ionosphere, wave scattering by density irreglarities) and magnetospheric propagation and characteristics of ducts.     

Variations of the lower ionosphere parameters measured by the resonance scattering method

This paper presents the results of diagnostics of the lower ionospheric parameters by using the resonance scattering method (RSM) of radio waves on periodic artificial irregularities (PAI). The following ionospheric parameters have been measured by the heating facilities ‘Zimenki’ and ‘Sura’: electron density including the E-F interlayer valley and the lowest height of the ionosphere; the velocity of the vertical wind; the relaxation times of the PAI characterizing coefficients of ambipolar diffusion, and coefficients of attachment and detachment of electrons from negative ions. Variations of these parameters have been considered as a function of height, local time, season, solar activity as well as being induced by passing acoustic-gravity waves and by ionospheric disturbances.     

Measurement of the eddy diffusion coefficient of the middle atmosphere from a balloon at low latitude

A great deal of uncertainty exists concerning the distribution of the vertical eddy diffusion coefficient of the middle atmosphere. A new technique has been developed in this laboratory in which chemical clouds are released in the middle atmosphere from a balloon platform. The expansion of the cloud is monitored by ground photography, from which the value of the eddy diffusion coefficient is calculated. The experiment was successfully tested on 9 March 1985 at Hyderabad (17.5°N, 78.6°E), India. The value of the coefficient was found to be of the order of 10⁴ cm² s⁻¹ in the altitude range 10–20 km.     

Diagnostics of the lower ionosphere by the method of resonance scattering of radio waves

Results are presented on measurements of lower ionosphere parameters (electron concentration, coefficient of ambipolar diffusion, etc.) by the method of resonance scattering of radio waves by periodic artificial irregularities. The method of resonance scattering is based on the generation of periodic irregularities in the ionospheric plasma by powerful radio emission and investigation of the characteristics of the back scattering of diagnostic radio waves by these irregularities.