Waveguide Parameters of 19.8 kHz Signal Propagating over a Long Path

The amplitude and phase of 19.8 kHz signal from navigational transmitter located in North West Cape, Australia, recorded at Suva, Fiji, have been utilized to determine the waveguide mode parameters. The propagation path is mixed over land and sea having Transmitter-Receiver Great Circle Path distance 6.7 Mm. The experimental values of the parameters were found to be consistent with the theoretical values calculated using the mode theory of VLF wave propagation in the waveguide.


Introduction
Extremely Low Frequency (ELF: 3-3000 Hz) and Very Low Frequency (VLF: 3-30 kHz) signals propagate to great distances in the waveguide bounded by the Earth's surface (ground or sea) and the lower ionosphere (∼60-140 km), consisting of the D and E-regions. This atmospheric waveguide is known as the Earth-Ionosphere Waveguide (EIWG). The D and E-regions are too high for balloons and too low for the satellite measurements. Radio sounding does not work, particularly at night, since electron densities in this region are low to reflect high frequency radio waves with ionosondes and incoherent radars. Therefore, the lower ionosphere remains the least studied region of the Earth's atmosphere [1]. The conductivity of the D-region increases exponentially with height which may also vary with latitude and longitude along the propagation path. At ELF/VLF, the Earth's surface and the lower ionosphere act as good electrical conductors. Electromagnetic waves reflect when incident upon conducting boundaries and are guided along the conducting structures.
The EIWG offers very little attenuation hence the guided ELF-VLF signals can be received literally around the world [2] and have long been used for long-distance communication, positioning and timings. Communications with submarines immersed in the conducting sea necessitated the use of VLF waves and later ELF waves, with their comparatively large skin depths in sea water [3]. Yokoyama and Tanimura [4] first observed diurnal variation of the amplitude of VLF (17.7 kHz and 22.9 kHz) signals propagated over long distances (>5000 km). Diurnal variations of VLF transmitter signals show amplitude minima associated with phase steps during sunrise and sunset transitions between transmitter and receiver [4][5][6] which could not be explained by single mode propagation theory. Wait [2] suggested that multiple modes are needed to explain VLF propagation in the EIWG. Crombie [5] proposed a model based on two modes being present in the nighttime portion of the propagation path and only one mode in the daytime portion of the propagation path in the EIWG. Clilverd et al. [6] have presented some interesting long-term (1990-1995) studies of VLF propagation over a long North-South path from Cutler (USA), Marine (NAA transmitter, 24 kHz) to Faraday, Antarctica. They found the timings of the minima in the received signal strength to be remarkably consistent from year to year. They also found that the timings of minima were consistent with modal conversion taking place as the daynight boundary (terminator) crossed the propagation path, at specific, consistent locations. The stability and reproducibility of the received amplitude and phase of VLF navigational transmitter signals make VLF propagation a useful tool for long-range communication and navigation system.
Purpose of this paper is to estimate the waveguide parameters from the observations of 19. 8

Theoretical Background
When a VLF wave is transmitted from a vertical antenna, a number of modes are excited in the EIWG. The electric field strength of the signal propagating in the waveguide can be represented by a superposition of discrete waveguide modes [2,7] expressed as: where E o is a constant depending on radiated power and frequency of operation (ω), h is the height of EIWG, A m is the excitation factor of the m th mode, α m is the attenuation rate of the m th mode, d is the distance between transmitter and receiver, t is the time, and v m is the phase velocity of the m th mode. The EIWG attenuation during day is higher than during the night and higher order waveguide modes have a larger spatial attenuation rate and are less efficiently radiated. The waveguide parameters can be estimated involving two modes in the night and one mode in the day. At large distances from transmitter (>5000 km), the mode 1 with lowest attenuation is dominant in the intensity. The diurnal change of the signal amplitude (E DN ) is defined as the ratio of nighttime to daytime signal strengths due to first mode. Using the expression of modal superposition (1), the value of E DN in terms of dB can be obtained as: where N and D refer to night and day. A N (dB) and A D (dB) are the excitation factors in dB given by A N (dB) = 20 Log 10 (A N ) and A D (dB) = 20 Log 10 (A D ) respectively. The diurnal phase shift normalized over 1 Mm (= 1000 km) represented by φ DN (μs Mm −1 ) is expressed as [7]: where c is the speed of light in free space and v is the phase velocity of VLF waves with mode 1 in the EIWG given by The modal interference spacing (D MS ) between the successive signal minima (fades) at receiver can be simply defined by considering the interference of two modes [5]: where λ 1 and λ 2 are the waveguide wavelengths of the first and second order modes, respectively. The Equation (4) can be simplified to as given by Crombie [5]:

Observational Data and Results
Experimental set-up consists of a short (1.5 m) whip antenna, pre-amplifier, VLF service unit coupled with preamplifier, and Software based Phase and Amplitude Logger termed "SoftPAL". The whip antenna receives the vertical electric field component of transverse magnetic (TM) mode of the VLF propagation. SoftPAL is a state of art data acquisition system recently (2006) developed by AD Instruments, New Zealand, which can log amplitudes (in dB above 1 μV/m) and phases (in degrees) of seven MSK (minimum shift key) VLF transmitter signals continuously with the time resolutions ranging from 10 milliseconds to 10 seconds, using a GPS based timing. Transmitter Receiver Great Circle Path (TRGCP) for NWC-Suva is 6.696 Mm. Typical diurnal variation of amplitude and phase of the signal on 10 December 2006, as recorded by SoftPAL, is presented in Figure 1. The daytime signal strength is larger than the nighttime which could be mostly due to the higher attenuation of EIWG for the daytime VLF propagation [8]. However, the received intensity of the signal depends on a number of factors including the excitation of modes, path attenuation, surface conductivity, ionosphere height, and the direction of propagation. The attenuation of VLF signals for the east-to west propagation is greater than for the west-east propagation and the difference in the attenuation decreases as the frequency increases which probably disappears above about 20 kHz [9]. Three amplitude minima labeled as SR 1 , SR 2 , SR 3 and SS 1 , SS 2 , SS 3 , respectively, are observed during sunrise and sunset transition along the TRGCP. As seen from Figure 1, the phase retards during the night and the rapid change of phase takes place at the time of the amplitude minima. Since D-region ionization is controlled by solar radiation, the reflection height of VLF waves, that is, height of the EIWG, at night, increases by about 10-15 km from its daytime value. As a consequence of the increase in the reflection heights, the phase velocity of the wave decreases and the signal is delayed (retarded) at the receiver. The amplitude data for selected 15 days in December 2006 are overplotted in Figure 2 to indicate the reproducibility of the amplitude values over a 24 hour period. The phase values are not plotted as the phase builds up over the day and deviates from reproducing similar values over the days. However, the stepwise phase advances coincides very well with the amplitude minima over any day and form of variation is almost the same. The signal variability is larger during the night than during the day indicating that propagation path is more stable in the day. The generation of successive minima at the receiver during sunrise and sunset transitions is generally accepted due to the destructive interference in the superposition of propagating modes with each other which may involve more than two modes in the night and dominant day mode at the terminator converted into a series of nighttime modes [6]. It can be seen from Figure 2 that the depth of fading varies from 20-22 dB during sunrise and 3-4 dB during sunset transitions along the TRGCP and fading is better shaped during sunrise transition. This could be due the less attenuation of waves propagating in the nighttime portion of the EIWG from transmitter to receiver, the shape of the lower ionospheric boundary and resulting scattering patterns, and the destructive interference of different number of modes at receiver during sunrise and sunset transitions. Also over the course of the propagation path, energy couples between modes which may result in the variation of depths of minima with the movement of the terminator between receiver and transmitter. Considering three different models for conductivity of the lower ionosphere, Devi et al. [10] carried out a waveguide mode analysis of 16 kHz VLF wave travelling great distances, to study its propagation parameters. They showed that for waveguide height less than 70 km, the change in phase due to changes in the phase velocity is smaller when ionosphere is assumed to have finite conductivity than when it is infinitely conducting. Using the data from Figure  theoretical values of the E DN and φ DN might be mainly due to the existence of higher modes (>1) particularly in the nighttime which have not been considered in the calculations for simplicity. Bainbridge and Inan [12] using the Long Wave Propagation Capability code for a long propagation path (NAA to Cheyenne, USA) found that the lowest order quasi-TM mode dominates during the daytime, but during the nighttime there are 4 to 5 significant modes with phase velocities ranging from 0.99-1.05 c and attenuation constants ranging from 0-10 dB/Mm. Using the amplitude and phase measurements of 16 kHz signal from Rugby (UK) VLF transmitter in India having TRCGP of 6.681 Mm, Joshi and Iyer [13] estimated the values E DN and φ DN , 14.0 dB and 7.85 μs Mm −1 , respectively. The value of φ DN obtained here compares well with results of Joshi and Iyer [13], whereas the value of E DN is lower by about 50% when compared with their value of E DN . The day-to-day variability in the times of minima is less during sunrise as compared to that during sunset ( Figure 2). Therefore, sunrise minima have been used to estimate the value of D MS using equation (5). From Figure 1, the separation between sunrise minima is found to be 77.5 minutes which corresponds to D MS of 2,150 km. The theoretical value of D MS from equation (5) for h = 90 km comes out 2,140 km which matches very well with the experimental value. The values of D MS estimated here both theoretically and experimentally are consistent with the results of Crombie [5] and Lynn [14] for the east-west VLF propagation in the frequency range of 13-22 kHz.

Conclusions
Preliminary results presented in this letter on the waveguide mode analysis of VLF waves at 19.8 kHz over a long propagation path, indicate a good consistency between experimental and theoretical values of waveguide model parameters. Detailed study of these parameters during different seasons by considering more modes in the nighttime will further enhance the understanding of this subject and appropriate the waveguide parameters at 19.8 kHz for NWC-Suva path.