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.     

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.     

Modelling of transequatorial propagation of HF radio waves

In the geometrical optics approximation, a synthesis oblique ionogram of ionospheric and magnetospheric HF radio wave signals propagating between magnetic conjugate points has been carried out. The magnetospheric HF propagation is considered for a model of the waveguide formed by field-aligned irregularities with depleted electron density. The characteristic peculiarities of the magnetospheric mode have been determined: (i) strong disperion of the group delay with a frequency at 14–18 MHz, from − 1.4 to 0.6 ms/MHz for magnetically conjugate points at geomagnetic latitudes φ = 30°, 40° and 50°, respectively, (ii) spreading ∼ 1–2 ms, and (iii) a possibility of propagation between magnetic conjugates points at moderately low geomagnetic latitudes φ₀ ∼ 30–40° at frequencies exceeding 1.5 times the maximum usable frequency (MUF) of multi-hop ionospheric propagation.     

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.     

Dispersion of chirp pulses by the ionosphere

Nonlinearity of the phase time delay vs frequency of the ionospheric channel results in frequency dispersion. This distorts wideband signals and leads to amplitude reduction and ‘elongation’ of narrow pulses. Its effect on frequency modulated continuous wave signals (FMCW or chirp) is to broaden the width of the compressed pulse and to produce a chirp signal at the output of the detector instead of the ideal sine function. Several workers have studied this distortion and derived expressions which related the width of the compressed pulse to the first order derivative of the group time delay vs frequency. These expressions are limited to the case when the bandwidth of the channel is greater than the bandwidth of the transmitted chirp signal. In this paper a generalized expression for the time-bandwidth product of the chirp signal at the output of the detector is derived. Numerical calculations for the width of the compressed pulse vs its time-bandwidth product for various window functions are presented. The results are then applied to experimental data obtained over a short skywave radio link.     

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.     

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.     

Ionospheric ELF radio signal generation due to LF and/or MF radio transmissions—I. Experimental results

A new non-linear ionospheric effect has been discovered from the analysis of ground based FLF/VLF radio data obtained in Scandinavia during 1979. Apart from signals of natural origin, timing signals (six pips which occurred on the hour) were received. The pips of frequency 1000±0.5 Hz, duration 105 ± 8 ms, and field strength ~ 0.1 pT at Sodankyla, Finland, exhibit a favoured source location ~ 150 km south-south-east of Sodankyla. A close association between the reception of these ELF pips and the auroral electrojet is demonstrated by the positive temporal correlation between ELF pip generation and periods of enhanced local magnetic activity, and also by the spatial correlation between source location and the latitude over F'inland at which riometer absorption is a maximum. In the evening and midnight sectors the latter is interpreted as indicating the electrojet position. The originating signals are shown to emanate from one or more Soviet LF/MF broadcast transmitters, all of which are several hundred kilometres or more from the favoured generation region.     

Anomalous interference in Omega VLF wave propagation on east-to-west equatorial paths

The phase of Omega Haiku (Hawaii U.S.A.) signals has been measured at 13.6 and 10.2 kHz on-board ship from Tokyo to Fremantle, Australia, during the early parts of the winters of 1979 and 1980. Short-term (∼- 1 h) fluctuations are observed on the phase of the Haiku signals as received around the geomagnetic equator at a distance of 8000 km west of the transmitter. Phase cycle slippings take place frequently in association with the phase fluctuations, the occurrence frequency of which is a maximum at 6 S geomagnetic latitude. These propagation anomalies are a consequence of the passage of the ship through an interference pattern the spacing distance of which is about 11 km at 13.6 kHz —one-half of the wave-length of the transmitted wave. It is concluded that the interference is caused by the signal propagated directly from the transmitter and the long-path signal propagating a distance of 32.000 km in the west-to-east direction. This result implies that the attenruation of the east-to-west propagating Omega wave is anomalously great at the equator, in pood agreement with calculated values based on the Galejs' anisotropic waveguide model.     

LHR whistlers and lhr spherics in the outer ionosphere

Special types of VLF signals, which follow whistlers and spherics and have an anomalous dispersion near the lower hybrid resonance (LHR) frequency, have been observed on the low-altitude Intercosmos satellites. These signals have been named LHR whistlers and LHR spherics, respectively. A mechanism is suggested for the formation of their spectra, based on the peculiarities of quasi-resonance wave propagation at frequencies near the LHR frequencies. It is shown that the large dispersion observed may be accounted for by a significant increase in the propagation time of the wave as its frequency approaches the maximum in the LHR frequency profile.     

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.     

The effect of the mid-latitude ionospheric trough on whistler mode ducting during magnetic storms

Whistler-mode signals observed at Faraday, Antarctica (65° S, 64° W, Λ=50.8°) show anomalous changes in group delay and Doppler shift with time during the main phase of intense geomagnetic activity. These changes are interpreted as the effect of refracting signals into and out of ducts near L=2.5 by electron concentration gradients associated with edges of the mid-latitude ionospheric trough. The refraction region is observed to propagate equatorwards at velocities in the range 20–85 ms⁻¹ during periods of high geomagnetic activity (K_p ≥ 5), which is in good agreement with typical trough velocities. Model estimates of the time that the trough edges come into view from Faraday show a good correlation with the observed start times of the anomalous features. Whistler-mode signals observed at Dunedin, New Zealand (46° S, 171° E, Λ=52.5°) that have propagated at an average L-shell of 2.2 (Λ=47.6°) do not show such trough-related changes in group delay. These observations are consistent with a lower occurrence of the trough at lower invariant latitudes.     

The polarisation of whistlers received on the ground near L = 4

The observed polarisation of the horizontal magnetic components of whistler mode signals received at Halley, Antarctica (L≈ 4.3), is in many cases that expected from a simple model of the transionospheric and sub-ionospheric propagation in the southern hemisphere; i.e. right-hand elliptical (field vectors rotate clockwise, looking towards the source) for ionospheric exit points close to the receiver, tending towards linear for more distant exit points. This suggests it may be possible to use the observed polarisation to estimate the propagation distance. However, in other cases, in certain frequency ranges, left-hand elliptically polarised signals have been observed. More realistic models do predict polarisation reversals at certain frequencies and exit point to receiver distances, but not over such a wide frequency range as has sometimes been observed. Also, in some cases, signals with nearly right-hand circular polarisation have been observed for exit points at large distances where linear polarisation would be expected.     

Experimental estimates of electron density variations at the reflection height of VLF signals

To study the behaviour of the electron concentration at the reflection level of very low frequency (VLF) waves, two years of phase and amplitude records of the 12.9 kHz signals emitted from Omega-Argentina (43.20°S; 294.60°E) and received at Tucumán (26.90°S; 294.70°E) have been used. The experimental results are compared with values derived from the International Reference Ionosphere model (IRI-79). The experimental data show a seasonal variation not predicted by the model. Differences are explained in terms of changes of night-time atomic oxygen concentration, which control the electron density profile at the base of the night-time D-region, not taken into account in the IRI model. Values of atomic oxygen necessary to explain VLF data are comparable with published data.     

Mid-latitude range spread and travelling ionospheric disturbances

This paper first discusses some early results (most hitherto unpublished) on off-vertical reflections which result from tilted isoionic contours associated with the passage of travelling ionospheric disturbances (TIDs). Azimuth of arrival and zenith angle information on these returning signals is discussed together with the role of these signals in producing both resolved and unresolved range spread on ionograms. Some ray tracing through a model ionosphere which incorporates wavelike structures is shown to predict fixed-frequency patterns (on recordings of virtual height vs. time) of converging and diverging satellite traces similar to those observed in practice. New experimental evidence is then presented to suggest (from the N(h) analyses performed) that TID wavetrains of several cycles are effective in producing spread traces on ionograms by specular reflections from the tilted surfaces of each cycle of these wavetrains. Ionograms from a modern ionosonde show that range spread consists primarily of discrete satellite traces even though at times these traces overlap and vary in intensity as a function of frequency thus creating a diffuse range spread ionogram.     

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.     

Spectral characteristics of high frequency waves backscattered by small scale F-region irregularities: evidence of strong sub-auroral ion flow

The spectra of high frequency waves backscattered at night by small scale (10–20 m) sub-auroral F-region irregularities often exhibit large Doppler shifts and widths in the local time sector 2000–2400. After local midnight the Doppler shifts and the widths of the spectra decrease rapidly. We present examples of experimental data, obtained with the two coherent backscatter radars of the EDIA¹ experiment, showing the spectral characteristics just mentioned. From the Doppler shift measured at the two sites we deduced the perpendicular velocity of the irregularities, which can reach values as high as 2000 ms ⁻¹. These observations are interpreted using results of theoretical models which predict strong sub-auroral ion flow in the trough region.     

Plasmaspheric zonal electric fields and coupling fluxes from the Dunedin VLF Doppler experiment,for 180°E,L ≈ 2.3, at solstice and equinox

Group delays and Doppler shifts from ducted whistler-mode signals are measured using the VLF Doppler experiment at Dunedin, New Zealand (45.8°S, 170.5°E). Equatorial zonal electric field and plasmasphere-ionosphere coupling fluxes are determined for L ≈ 2.3 at June solstice and equinox during magnetically quiet periods. The general features of the electric field measured at Dunedin agree with those predicted from ionospheric dynamo theory with a (1,−2) tidal component. Some seasonal variations are observed, with the electric field measured during equinox being smaller and predominantly westward during the night. The electric field at June solstice is also westward during the evening and for part of the night, but turns sharply eastward during the pre-dawn and dawn period at the duct entry site. The June electric field appears to follow a diurnal variation whereas the equinox electric field shows a possible 4-hourly periodic variation. Seasonal variations in the neutral wind pattern, altering the configuration of the ionospheric dynamo field, are the probable cause of the seasonal differences in the electric field. The seasonal variation of the coupling fluxes can be explained by the alteration of the E x B drift pattern, caused by the changes in the electric field.     

Regimes of ionospheric turbulence from fractal analysis of satellite radio signal scintillations

Samples of amplitude scintillations of the radio signal from a geostationary satellite obtained at a midlatitude station near Irkutsk were processed. For calculating the fractal dimensionalities the Grassberger and Procaccia [(1983) Physica D9, 189] algorithm was used. Results of the data processing tend to divide into two groups. One group includes those realizations for which it was possible to obtain reliable estimates of dimensionality. Three of the seven realizations considered were in this group, and the fractal dimensionalities were found to be low (3.12 4.5). The other data fall within the second group; a reliable estimate of dimensionality for them is unobtainable in terms of the method used. We suppose that this is attributable to the high dimensionality of the process. Power spectra of the signals of this group are close to those with an exponent of −2. The spectra of the signals of the first group are markedly steeper. On the basis of the data analyzed it is supposed that there exist two modes of ionospheric turbulence in midlatitudes, namely the mode with low dimensionality typical of localized turbulent processes, and the mode with high dimensionality typical of homogeneous turbulence that covers an extensive region of the ionosphere.     

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.