Lightning-induced perturbations on VLF subionospheric transmissions

Phase and amplitude perturbations on VLF subionospheric transmissions from transmitter NWC to Dunedin have been studied on both MSK frequencies and at spaced receivers, 9 km apart. In any one event (a ‘Trimpi’) the phase and amplitude perturbation can be expressed in terms of a perturbation phasor. This is generally believed to be the result of lightning-induced electron precipitation (LEP) producing a localized increase in ionization near the normal reflection height for subionospheric (waveguide) VLF waves. Most of the Trimpis received on the NWC-Dunedin path can be best explained if the LEP ionization is sufficiently localized so that it acts as a scattering centre for the subionospheric VLF wave from the transmitter. It is then this scattered wave or echo at the receiver which makes the perturbation phasor. We call these ‘echo Trimpis’. The phase of the echo relative to the direct signal will differ on spaced antennae if the angle of arrival of the two signals differ. Similarly, this relative phase will vary with frequency if the group delay of the signals differ. Thus measurement of these differences allows location of the scattering centres, and so too the LEP. Locations made show a significant grouping in a region where the lightning intensity is high. This and other features strongly suggest that these echo Trimpis originate from local (southern hemisphere) lightning. This and other reasons are suggested to explain the high proportion of echo Trimpis on this path.     

Standing waves in VLF transmissions produced by geographical discontinuities

It has generally been accepted that the fields constituting the waveguide modes reflected from natural discontinuities in the lower boundary of the earth-ionosphere waveguide have a negligible amplitude compared to the fields constituting the incident modes at VLF (Kirchoff's approximation). Recent measurements, in the region of the South Island of New Zealand, show cases where such reflections play a significant role in the determination of VLF field strength and phase.     

VLF transmissions observed at anomalously high latitude above the ionosphere in the conjugate hemisphere

Ariel 3 and 4 satellite observations of the GBR 16 kHz and NAA 17.8 kHz transmissions above the ionosphere in the conjugate hemisphere show that their wave-fields generally show a rapid reduction in signal strength for geomagnetic latitudes greater than 55°–60°. Sometimes, however, the signal strength has been observed to be high in the invariant latitude range > 60°. At certain times during these observations, the signal showed clear evidence of amplification, whilst at other times the pattern of signal strength was displaced to higher latitude with the signal strength integrated over latitude being unchanged from that normally observed.It is shown that the plasmapause can guide both the NAA and GBR signals but that the efficiency of this guiding depends on the plasmapause position. The important condition is found that the plasmapause must be situated sufficiently equatorwards that half the equatorial electron gyrofrequency at the plasmapause position is greater than (or approximately equal to) the transmitter signal frequency. Ray-tracing calculations in a realistic magnetosphere model indicate that for the 16 kHz GBR signal, the efficiency of guiding falls off for L_pp, (the L-value of the plasmapause) > 3.0 and guiding effectively ceases for L_pp > 3.5.Guidance by the plasmapause results in a wave-field at higher latitude than for non-guided propagation. This will only occur when, following geomagnetic storms, the plasmapause position is at a sufficiently low L-value. This is in agreement with the experimental observations of anomalously high latitude signal reception following strong magnetic storms (K_p ≥ 4+).     

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.     

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.     

In-situ observations of magnetospheric electron scattering by a VLF transmitter

High resolution pitch angle measurements of outer zone electrons in the energy range 12 keV−1.6 MeV were obtained at high altitude in the region of the high power VLF transmitter UMS [300 kW radiated at 17.1 kHz (Watt A. D., 1967, VLF Radio Engineering, Pergamon Press, Oxford)] while a resonant wave-particle interaction was in progress. Additional complementary electron measurements in the range of 36–316 keV were obtained in the drift loss cone by another satellite at low altitude along the drift path 75° east of the interaction region. The data from the low-altitude satellite confirm that UMS was precipitating particles in the inner zone, in the slot, and in the outer zone at the time that the high-altitude satellite was obtaining its data. The high-altitude pitch angle distributions indicate that, for this event, two types of scattering interactions were in progress. Particles with small pitch angles, up to 17.2° at the Equator, were being removed, resulting in an enhanced loss cone. Particles which were mirroring between 6500 km and the altitude of the spacecraft (7200) km were also being strongly scattered, resulting in a relative minimum in the pitch angle distribution around 90°. The data are interpreted as indicating that a cyclotron mode interaction with UMS waves was precipitating electrons with equatorial pitch angles up to 17.2° and that another process, perhaps electrostatic (ES) waves arising from the UMS radiations through a mode-conversion process, was present in the region above 6500 km and was efficiently scattering those particles which mirrored in that region     

Periodic and quasiperiodic VLF emissions

A review of periodic and quasiperiodic VLF emissions observed at ground-based stations and in the Earth's magnetosphere is presented. Emissions with periods below 10 s are divided into three main groups: periodic emissions, hisslers and pulsing hiss, while the quasiperiodic emissions with periods above 10 s are divided into two main groups: QP1, which are definitely related to geomagnetic pulsations, and QP2, which are not obviously related to them, based on ground observations. However, it is pointed out that all types of quasiperiodic emissions with periods over 10 s observed onboard satellites show some association with geomagnetic pulsations which suggests that both QP1 and QP2 can have similar mechanisms for their generation. Different approaches to the theoretical modelling of periodic and quasiperiodic emissions are discussed.     

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.     

Experimental daytime VLF ionospheric parameters

Measured field strengths from VLF transmitters are used to determine improved daytime values of ionospheric parameters to enable improved VLF propagation predictions. These parameters are the traditional H′ (height in km) and β (sharpness in km⁻¹) as used by Wait and by NOSC in their Earthionosphere waveguide computer program. They are found by comparing the predictions of the NOSC program with the observed VLF field strengths over both long and short paths.Experimental observations from two nearly north-south paths are used to determine the solar zenith angle dependence of both H′ and β for low latitude (or summer mid-latitude) conditions. These results are then used to predict the daytime variations in VLF field strengths with solar zenith angle (and hence time) on other suitable paths and good agreement is found with measurements made on these paths.The absolute value of β for overhead Sun is found to be 0.45 km⁻¹ and is principally determined by the attenuation on the very long, west to east, fully sunlit, 14.4 Mm path from NWC (Australia, 22°S) to San Francisco (37°N), after applying small corrections for the solar zenith angle variations along the path at midday. Further support is obtained from results from the 8.6 Mm path NDT (Japan) to San Francisco, an 8.2 Mm path NPM (Hawaii) to New Zealand, and an east to west 7.5 Mm path from NPM to Townsville, Australia. The conditions studied are solar maximum. The frequencies studied are 15–30 kHz.     

Recent findings on VLF/ELF sferics

Our recent activity on VLF/ELF sferics is reviewed, and this paper is composed of two parts. The first is a new method of estimating the propagation distance and ionospheric reflection height by using sophisticated signal processing for the dispersive tails of tweek sferics near the cutoff frequencies. In the second part, we have carried out the first attempt to apply our field-analysis direction finding (developed for whistlers) to tweek sferics; the corresponding experimental results and their interpretation are presented.     

Multi-station VLF direction-finding measurements in eastern Canada

The directions of propagation, in the earth-ionosphere waveguide, of multi-component two-hop whistlers recorded on 10 July 1972 by four VLF goniometer receivers in eastern Canada have been determined. Using the bearings of these great circle paths, triangulation of several whistler exit-points has been accomplished. The L-values of the whistler exit-points determined by this method are systematically lower than those expected from their nose frequencies, by ~ 0.6. Various explanations are discussed for this effect. The most satisfactory is that the whistler waves leave through the side of the ducts (in which they had propagated for most of their path through the magnetosphere) at an altitude of a few thousand kilometres, and then are refracted to lower L-values before exiting from the lower ionosphere. The results are consistent with both the duct termination altitude predicted by Bernhardt and Park (1977) for the appropriate conditions and also with the observed upper cut-off frequency of the whistlers.     

On the frequency shift in modulated VLF emissions

The nature of the shift between the frequency of geomagnetic pulsations registered on the Earth and the modulation frequency of VLF emissions registered on board satellites is analyzed. The results obtained are compatible with the mechanism of VLF emissions modulations founded on the interaction between magnetosonic waves propagating perpendicular to the geomagnetic field and the energetic electrons exciting VLF emissions. It is shown that the frequency shift considered can be used for pulsations source localization in the magnetosphere.     

Temperatures in the plasmasphere determined from VLF observations

Satellite and ground-based VLF recordings were made at SANAE, Antarctica from 1976 to 1979. In this paper we combine ground and satellite observations to determine temperatures in the plasmasphere. Scale heights in the plasmasphere are determined at high altitudes using a diffusive equilibrium model and measurements of equatorial electron densities and densities at about 3000 km. The temperatures corresponding to these scale heights show a gradual increase with increasing L-value and sharp increases of about 2000 K just inside the plasmapause.     

ELF and VLF wave generation by modulated HF heating of the current carrying lower ionosphere

This paper presents further experimental results on ionospheric current modulation, using powerful amplitude modulated HF waves produced by the new heating facility at Ramfjordmoen near Tromsø, Norway. As a result of the current modulation, waves in the ULF, ELF and VLF range can be efficiently generated. The experiments discussed here cover the range from low ELF up to 7 kHz. The observed signal strengths are of the order 1 pT. Decomposition of the received ELF/VLF waves into R- and L-mode shows that both modes are usually of comparable strength. The signal strength as a function of modulation frequency shows pronounced maxima at multiples of approximately 2 kHz. The paper also presents a brief theoretical discussion of the processes involved in the generation of ELF/VLF waves by HF induced current modulation.     

Re-radiation of VLF radio waves from mountain ranges

VLF signals transmitted from Hawaii, Japan and Australia (NPM, NDT, NWC) and scattered from large mountain ranges such as the Andes or Rockies are monitored in New Zealand in the presence of the much larger direct signals. Measurements of the amplitudes of these indirect signals are reported in the range 0.l-2.0μVm⁻¹.Sudden changes in the refractive index along the path, such as those caused by abrupt changes in ground conductivity, are shown to give rise to reflected amplitudes much smaller than those measured. However, the measured amplitudes are found to be comparable with reasonable estimates of the re-radiation from typical surface currents induced by the incident wave in mountains of suitable height, ground conductivity and steepness.     

A simple theoretical model of the diffraction of VLF electromagnetic waves by the Antarctic icecap

The Antarctic continent has been modelled as a spherical cap whose pole is coincident with that of the South Pole, which totally absorbs VLF radio waves attempting to propagate over it. The propagation of Omega navigation signals around this model icecap has then been computed using Kirchhoff diffraction theory. Spherical caps extending to 66.5 and 75.5°S have been found to accurately model the signals from Omega La Reunion and Argentina, respectively, received on flights between Christchurch, New Zealand and Scott Base in Antarctica, up to the boundary of the theoretical icecap. These model icecaps were found to be good fits to the boundary of the Antarctic continent, when measured at the 1–1.5 km contour of ice thickness, in the region where the VLF waves diffracted around the icecap. The good agreement obtained between the experimental field strength data and those computed theoretically, using only simple diffraction theory, suggests that coastal refraction plays at most only a secondary role in circumpolar propagation.     

Reflection of VLF radio waves from distant mountain ranges

A computer cross-correlation technique is being used to determine the group delays and directions of arrival of man-made subionospheric VLF signals which have reached the receiver by paths other than the direct great circle path. The 200 baud MSK signals transmitted by NWC, NPM and NLK allow time resolution to at least 5 ms and, with 15 min of integration, the sensitivity can be as low as about 0.1 μV m⁻¹ in quiet conditions. Reflections from the Andes, the Rockies and the mountains of S.E. Asia have now been identified at Dunedin, New Zealand. Round-the-world and round-the-world-the-other-way signals have also been observed.     

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.     

Waveguide model analysis of VLF wave propagation at 16 kHz

Preliminary results of measurement of phase and amplitude of 16 kHz VLF signals propagated in the Earth-ionosphere waveguide are presented. From the measured values of diurnal phase, amplitude changes and sunrise fadings the relevant waveguide parameters are estimated.     

Doppler shift pulsations on whistler mode signals from a VLF transmitter

Whistler mode signals from the NAA transmitter (24 kHz) received at Faraday, Antarctica are processed to obtain the Doppler shift at a much higher time resolution than has previously been possible. This has allowed the observation of pulsations of about 13 mHz frequency which are believed to be associated with hydromagnetic waves in the magnetosphere. The pulsations are observed separately on signals with a number of discrete group delay features that can be interpreted as individual whistler ducts. Using the measured pulsation phase over the array of ducts the phase velocity and wave normal direction of the hydromagnetic wave in the equatorial plane are estimated. The direction of propagation is consistent with a source on the dayside magnetopause.The association between whistler mode Doppler shifts and hydromagnetic waves has been reported before but not, as far as we are aware, using an experimental technique that allows measurements on individual ducts in order to determine the direction of propagation of the hydromagnetic wave.