Use of VLF transmissions in the location and mapping of lightning-induced ionisation enhancements (LIEs)

Lightning-induced ionisation enhancements (LIEs) are usually produced by short (~ 1 s) bursts of energetic electrons precipitated from the radiation belts in the process of amplifying whistlers. During their short life (~ 30 s) LIEs diffract or otherwise modify stable transmissions of VLF waves propagating in the two-dimensional Earth-ionosphere waveguide. This causes perturbations (‘Trimpis’) on the same time scale in the phase and amplitude of these VLF waves. Unless the LIEs are large (> 100 km) and smoothly varying (e.g., Gaussian distribution of ionisation enhancement) in the horizontal directions, the LIEs need not be on the great circle path (GCP) from VLF transmitter to receiver to produce Trimpis. Large and smooth LIEs produce ‘GCP Trimpis’, while small or structured LIEs produce ‘echo Trimpis’. The two can usually be distinguished, if Trimpi phase and amplitude are monitored and if the Trimpis are observed at several frequencies or on two or more spaced receivers simultaneously.If only GCP Trimpis are considered, the causative LIEs can be located and mapped by geometric optics using a network of receivers of sufficient density (spacing ~ 100 km) and a few transmitters. Provided all Trimpis are identified as GCP, their mere detection is sufficient for location. This is equivalent to locating the LIE ‘shadows’ cast onto arrays of spaced receivers by two or more transmitters. If this GCP identification is not made, or is just assumed, location and mapping (size estimation) errors can be quite large. At VLF (λ ~ 15 km) this geometric optics approach cannot be used to study the horizontal fine structure of LIEs since LIEs producing GCP Trimpis have no fine structure.Small or structured LIEs cast a diffraction pattern onto an array of spaced receivers. If both the phase and amplitude perturbation of echo Trimpis are measured at each receiver of the array, holographic techniques can be used to reconstruct the two-dimensional map or image of the causative LIEs. It is shown that, for a single system of one transmitter and a receiver array, this allows high resolution (~ 10 km) in the azimuthal dimension only. Equally high resolution in both horizontal dimensions can be achieved with two orthogonal systems. This technique works equally well on GCP Trimpis to map the causative LIEs (which are large and structureless) without incurring location errors thereby.     

VLF scattering from a column of ionization in the Earth-ionosphere waveguide

Choosing a highly idealized model, we analyze the scattering of VLF radio waves from a thin vertical column of ionization within the Earth-ionosphere waveguide. It is shown that mode conversion is produced if the column is limited in height and/or if the ionization is a function of height. However, the relative higher-mode content would be small in some situations such as in the case of the ionization produced by a cloud-to-ionosphere lightning discharge.     

Size and location of lightning-induced ionisation enhancements from measurement of VLF phase and amplitude perturbations on multiple antennas

A 600-km array of five Trimpi receivers (“elements”) has been set up in New Zealand broadside to the VLF (22.3 kHz) transmitter, NWC, some 6000 km west, with element separations varying from 8 km to 550 km. Although such a five-element array is inadequate for imaging of lightning-induced ionisation enhancements (LIEs) by VLF holography, or inverse scattering, estimates of LIE size and location can be made if the shape and form of the LIE can be guessed or assumed, with even fewer elements. With five elements, tests of the assumed model can be made as well.Owing to its transform properties, the simplest model to use for scattering inversion is the Gaussian LIE distribution. For this model, and for single mode propagation, an inversion process is derived here for the full range of LIE and path dimensions, ranging from those for which the receiver is in the diffraction far field to those in which “geometric optics” dominate. This inversion process has some validity for small LIEs of other shapes of simple form. For more extreme models, the dominance of geometry or diffraction can usually be established in individual cases which then allows simple scaling procedures to be used in scattering inversion.Some 70 Trimpi events were observed on all five elements during a single night in July. 1991 (late winter). These were used to determine LIE location and size, and to test the applicability of various LIE models. It was found that most LIEs that night occurred over the Tasman Sea near the great circle from the VLF transmitter, NWC, to Wellington, generally some 500 to 2000 km from Wellington, and with north-south dimensions of 100–250 km. Much longer east-west dimensions (oriented towards NWC) are suggested to account for the very strong Trimpis observed. While about half of these LIEs that night could have had a smooth lateral spread (e.g., Gaussian), the remainder required varying degrees of fine structure, from “flat” or Butterworth LIEs to multiple LIEs as might be expected from multiduct whistlers, to explain the observed diffraction pattern exhibiting maxima and minima as well as the wide angular range over which simultaneous Trimpis were observed.     

Very-low frequency signals (2–10 kHz) received at Palmer Station,Antarctica, from the 21.4 km dipole antenna at Siple Station, 1400 km distant

Radio signals transmitted from the unique experimental VLF transmitter at Siple Station (76°S, 84°W), Antarctica, as well as VLF signals from communication and navigation systems and waves that propagate in the ionosphere and magnetosphere in the whistler mode, are regularly received and analysed at Palmer Station (65°S, 64°W), Antarctica. The amplitude and polarization properties of the Siple signals are predicted using a ray optics analysis. The amplitude of the signal received from Siple varies with frequency; observed nulls in the signal spectrum, where thesignal amplitude/alls 5–10 dB below what might be expected, are explained by the ray analysis. The amplitude spectrum is observed to be very sensitive to ionospheric conditions. Whereas the arrival bearings of signals from VLF transmitters other than Siple are found to be within 5° of their expected values, which is consistent with their expected vertical polarization and the operation of the DF system, an approximately 90° anomaly in the apparent arrival bearing of the signals from Siple is attributed to the essentially horizontal polarization of the received signal. The anomaly is found to be consistent with the theory of operation of the DF system. Occasional anomalies greater than 90° are explained in terms of a combination of polarization error and a smaller multi-path error. Siple two-hop signals and whistlers propagating on a common magnetospheric path showed arrival bearings and other properties consistent with a path end point within 200km of Siple. This suggests that these signals were received at Palmer with essentially vertical polarization.     

Experimental studies of daytime LF radio propagation and their implications for D-region electron densities

Measurements are presented of interference phenomena in amplitude and phase of VLF and LF signals along propagation paths from central England to the Norwegian Sea. The data are interpreted by means of the ‘wave-hop’ propagation theory, incorporating full wave evaluation of ionospheric reflection coefficients with realistic D-region models. No published electron density profiles are found which completely satisfy the experimental data, but modified profiles are presented which provide a better fit to the observations.     

VLF propagation parameters derived from sferics observations at high southern latitudes

Sferics are electromagnetic pulses generated by lightning events. Their maximum spectral energy is in the frequency range below 15 kHz. These powerful natural VLF transmitters can be used to determine the propagation characteristics of the atmospheric wave guide between earth and ionospheric D layer along virtually every propagation path. A VLF-sferics-analyzer was operating at the German Antarctic von Neumayer Station from January to June 1983. This analyzer recorded sferics from distant lightning events in the frequency range between 5 and 9 kHz. The method of measurement is described. The data are evaluated, and the propagation characteristics of the atmospheric wave guide are determined as a function of azimuth and season. The result is compared with theoretical calculations. It is shown that the difference between west-to-east and east-to-west propagation is much smaller than theory predicts, indicating that the ionospheric D layer at high southern latitudes behaves less anisotropic with respect to VLF propagation than at mid-latitudes.     

The use of the phasor display in studying ionospheric radio echoes

The phase and amplitude of a radio pulse reflected from the ionosphere usually vary during the pulse. It is convenient to observe these variations using the X-Y mode of an oscilloscope to display the phasor of the echo. The variations are then seen as an oval or spiral shape traced out by the end point of the phasor. These shapes provide a sensitive method of detecting the presence of more than one echo, and are useful as a measure of dispersion.     

Rapid onset,rapid decay (RORD), phase and amplitude perturbations of VLF subionospheric transmissions

Rapid onset (few ms), rapid decay (~ls) perturbations or RORDs occur frequently on the west-to-east signal from NWC to Dunedin, more often than not with classic Trimpis. They do not appear on an NWC mimic signal directly injected into the antenna and so cannot be broadband bursts. There is no delay between the initiating sferic and RORD start, implying that they are produced not by whistler-induced electron precipitation but directly by lightning. Observations on a multi element array show that classic Trimpis and RORDs initiated by the same sferic usually come from measurably different directions, so the lightning-induced ionisation enhancements (LIEs) which cause them must be laterally displaced. They may also be vertically displaced to explain the differing decay rates (30s versus 1 s). We conclude that RORDs are VLF echoes from vertical columns of ionisation at around 40km altitude and having vertical dimensions of some tens of km and horizontal dimensions of 1–2km, since such a column would scatter sufficient signal to fit observed amplitudes. Cloud-to-ionosphere (CID) lightning discharges (also called “cloud-to-space” and “cloud-to-stratosphere” discharges) of these visible dimensions have been observed on mountain observatories and on board the Space Shuttle.     

Trimpi events on low latitude paths: An investigation of gyroresonance interactions at low L-values

We report on Trimpi events observed at Durban (L = 1.69, 29°53′S, 31°00′E) and investigate the efficacy of gyroresonance scattering in precipitating electrons into the atmosphere at low L (<2). The rate of occurrence of Trimpis at Durban is less than one per day. Our observations include a number of daytime events on OMEGA signals from La Reunion. Using the full relativistic equations of motion, a test particle simulation is employed to find the region in parameter space where large pitch angle scattering occurs. We find that at low L the conditions for pitch angle scattering are less favourable than at higher L (L ∼ 4). Resonant electrons have high (relativistic) energies, interaction times are of the order of milliseconds (T_i ∼ 5 ms) and large wave amplitudes (B_w ∼ 200 pT) are required at whistler frequencies to produce pitch angle changes of greater than 1°. Large pitch angle scattering is needed near Durban since particles near the loss cone will have been lost in the South Atlantic Geomagnetic Anomaly. We note that the radio frequencies transmitted into the magnetosphere from lightning are too low to give effective electron scattering at low L. We suggest an explanation for the low rate of occurrence of Trimpis at Durban.     

Lower ionosphere effect observed during the 30 June 1992 total solar eclipse

VLF radio signals (12.9 kHz) transmitted from Ω-Argentina (43°12′S, 65°24′W) were received in Atibaia, Brazil (23°11 'S, 46°33'W) during the total solar eclipse of 30 June 1992. The surface path of the totality crossed the VLF propagation path in the sunrise transition period causing a phase delay of 6.4 μs and an amplitude change of 1.3 dB. The ionospheric response to the Sun's obscuration was compared with the phase delays reported for several solar eclipses that occurred from 1966 to 1979. The results are mainly discussed in terms of the length of VLF propagation path affected. Some similarities between a sudden phase anomaly and a reversed eclipse effect are also raised.     

Satellite observations of signals from a Soviet mid-latitude VLF transmitter in the magnetic-conjugate region

The results of the Intercosmos-19 satellite experiment which receives VLF signals arriving via a magnetospheric path have been examined. The reception zone in the magnetically conjugate region (MCR) has been shown to be centred near the L-value (2.6) of the transmitter of 15 kHZ radio waves. The received signals arrive at the MCR with wave normal angles to the geomagnetic field, ψ, far from the resonance cone. These results indicate an effective amplification of the signal in the magnetosphere by 10–15 dB and effective ducting of VLF waves across the equatorial plane of the magnetosphere.     

Early Trimpi events from lightning-induced electric fields in the ionosphere: An alternative explanation

Two classes of ‘Trimpi’ modulation of VLF signals in the Earth-ionosphere waveguide have been identified in the literature. The more common type occurs l s or more after causative lightning strokes, the second in less than 100 ms. We explore the possibility that these early Trimpi events result from lighting-generated, electric field impulses lowering the mirror altitudes of trapped electrons. To overcome the mirror force on energetic electrons, upward-directed electric fields with strengths of a few tens of mV/m are required. This is well within the range of electric fields observed on sounding rockets above thunderstorms.     

On VLF radio wave reflection from distant mountain ranges—theory

We present an analysis for the reflection or scattering of VLF radio waves by gently sloping parallel ridges that may be well displaced from the great circle path between the transmitter and receiver antennae. It is shown that, under certain conditions, a spatial harmonic of the transverse profile of the terrain will produce a secondary signal that could be not far below the direct or primary signal. Such a resonant scattering mechanism may explain some of Thomson's [(1989), J. atmos. terr. Phys. 51, 339] observations in Dunedin, New ZeaLond.     

Notes on the diversity of the properties of radio bursts observed on the nightside of Venus

We report on further studies of radio wave bursts detected by the Orbiting Electric Field Detector (OEFD) on PVO in the nightside ionosphere of Venus. We have tested a total of 25 cases of wave burst activity for evidence of whistler-mode propagation to the spacecraft from impulsive subionospheric sources. As in a previous study of 11 of these cases (Sonwalkar et al., 1991), we find at least two distinct classes of events, one, mostly involving bursts at 100 Hz only, that passes certain tests for whistler-mode propagation, and another, mostly involving bursts in two or more of the four PVO narrowband channels (at 100 Hz, 730 Hz, 5.4 kHz, and 30 kHz), that fails to pass the tests. The subionospheric lightning hypothesis continues to be tenable as a candidate explanation for many of the 100 Hz-only events, but its plausibility could be better evaluated if mechanisms could be found to explain the existence of a significant number of 100 Hz-only cases that do not pass all the applicable whistler-mode tests, as well as the existence at a wide range of altitudes of multichannel cases that are clearly not propagating whistler-mode waves. The wideband bursts are often observed at altitudes above 1000 km and frequently occur in regions of locally reduced electron density. Those observed at high altitude (and possibly at low altitude as well) are believed to be generated near the spacecraft, possibly by an as yet unknown mechanism responsible for similar burst observations made near Earth and other planets.     

ELF/VLF propagation measurements in the Atlantic during 1989

The vertical electric field component was measured by a group of the Ukrainian Institute of Radio Astronomy on board the Professor Zubov scientific vessel during April 1989 at latitudes from 30°S to 50°N. Results of the amplitude measurements in the Atlantic of natural ELF radio signals and those from the VLF navigation system “Omega” at its lowest frequency of 10.2 kHz are given. Characteristics were obtained of the moving ship as the field-site for the ELF observations. Variations in the ELF radio noise amplitude recorded at tropical latitudes agree with the computed data for the model of three continental centres of lightning activity. The VLF results were obtained by the “beat” technique providing the simplest narrow-band amplitude registration. Range dependencies of the field amplitudes from A (Norway), B (Liberia) and F (Argentina) stations have been analysed. The VLF attenuation factor was estimated for the ambient day conditions along the four cardinal directions. This allowed the detection of a statistically significant attenuation difference between the east-west and west-east propagation paths. The VLF radio signal was also used as a probe to evaluate the effective height of the vertical electric antenna and to calibrate the ELF noise amplitudes.     

The diffraction of VLF radio waves by the Antarctic ice cap

Measurements of the amplitude and phase of VLF radio signals from the Omega transmitters on La Reunion Island and in Argentina have been made on routine Antarctic re-supply nights from Christchurch, New Zealand. It has been found that when the propagation paths to the transmitters cross the Antarctic ice cap, the direct path signals are very rapidly attenuated below the receiver noise level, the dominant signal source then being provided by the radio waves diffracting around the edge of the ice cap. These results have been made possible by the simultaneous use of the phase and amplitude data in a synthetic aperture antenna type analysis.     

Characteristics and frequency of occurrence of Trimpi events recorded during 1982 at Sanae,Antarctica

‘Trimpi’ amplitude perturbations on VLF signals received at Sanae, Antarctica, have been identified using a new computerised technique. Our survey of 1982 data, taken during magnetically disturbed times, shows that events of short duration (<25 s) constitute 60% of all events detected and that all events found are amplitude attenuations with deviations from quiescent levels ranging up to 90%. It is unusual, at Sanae, to observe the causative whistler with a Trimpi event. This, together with further evidence from Trimpi occurrence statistics, may suggest that the gyroresonant interactions responsible for some of the events occur with non-ducted whistler mode waves. A method for estimating the extent of the precipitation region is presented.     

The effects of ionospheric horizontal electron density gradients on whistler mode signals

Whistler mode group delays observed at Faraday, Antarctica (65° S, 64° W) and Dunedin, New Zealand (46° S, 171° E) show sudden increases of the order of hundreds of milliseconds within 15 minutes. These events (‘discontinuities’) are observed during sunrise or sunset at the duct entry regions, close to the receiver's conjugate point. The sudden increase in group delay can be explained as a tilting of the up-going wave towards the sun by horizontal electron density gradients associated with the passage of the dawn/dusk terminator. The waves become trapped into higher L-shell ducts. The majority of the events are seen during June-August and can be understood in terms of the orientation of the terminator with respect to the field aligned ducts. The position of the source VLF transmitter relative to the duct entry region is found to be important in determining the contribution of ionospheric electron density gradients to the L-shell distribution of the whistler mode signals.     

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.