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.     

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.     

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.     

The spatial coherence of travelling ionospheric disturbances at mid-latitudes

A multifrequency HF Doppler sounder and four spaced receivers were operated near Alma-Ata to form a three-dimensional array of reflection points of HF radio waves. The spacings of reflection points ranged from 5 to 80 km in the vertical and from 30 to 65 km in the horizontal. The purpose of the experiment was to estimate the spatial coherence of travelling ionospheric disturbances (TIDs). Estimation of the coherence length (the distance at which the coherence falls to e⁻¹) in both vertical and horizontal planes is carried out. The coherence often shows peaks at frequencies exceeding the Brunt-Väisälä frequency. Measurements of the slant coherencies have given the opportunity to study the coherence as a function of orientation.     

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.     

Elongation of midlatitude scintillation irregularities

This paper reports a study of the length of the midlatitude F-region irregularities which cause scintillations. The length of scintillation irregularities is usually measured by the Full Correlation Analysis (FCA) method. If observations from a typical 3-receiver array, with small spacing, are used the length is always found to he only a few kilometers. Basic ionospheric theory indicates that the irregularities should he much longer than this. In this experiment scintillations of beacon transmissions from polar-orbiting satellites were observed on a 3-receiver array with 310 m spacing, and a fourth receiver 4.6 km away was used to check the lengths measured. Analysis of the scintillation observations from the 310 m array used the FCA method. In all cases the FCA lengths, from the small array, were much shorter than those measured by the distant receiver. Measurements using the 4.6 km receiver gave an average ratio of 44.5 and a half-length (in ionosphere) of 9.4 km. It appears that random errors in the correlation functions cause the FCA to underestimate the length of the ground pattern of irregularities when the spacing of the receivers is much less than the pattern size.     

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.     

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.     

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.     

The height distribution of radio meteors: observations at 2 MHz

Conventional meteor radars, operating at wavelengths of around 5–15 m, are unable to detect high-altitude meteors due to the wavelength-dependent echo ceiling. It is suggested that the ‘missing mass’ in the 10⁻⁶–10⁻² g range of interplanetary material is in fact a high-velocity component which is normally undetected since it ablates at high altitude. This contention is supported by previous work. In this paper we describe measurements of the heights of radio meteors (limiting magnitude about +7) at a wavelength of 150 m (frequency 2 MHz), for which the echo ceiling is above 140 km. The resultant true height distribution is found to peak at ~ 104 km, about 10 km above the peak found by conventional meteor radars. The majority of meteors are detected at or above this peak, and substantial numbers are seen right up to 140 km. It is therefore concluded that the ‘missing mass’, comprising the vast majority of the meteoric input to the atmosphere, ablates well above 100 km.     

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.     

Heights of MF radar scatter (1986/87) and the wind field (55–95 km): Saskatoon,Canada

The Saskatoon MF radar (2.2 MHz) at 52°N, 107°W, has been used to measure the heights of occurrence of radar scatter during four seasons, and twelve months of 1986/87. Mean winds, and gravity waves are also available, by the spaced antenna method and from the same radar echoes. Certain heights, called elsewhere ‘preferred heights’, are identified near 60km, 70km, 75km in summer, and 80–86 km. Several layers have seasonal and diurnal variations. Associations with electron density gradients (rocket data), mean wind shear in summer, and gravity wave amplitude-minima in the equinoxes are effectively demonstrated. Case studies, involving 3 h data sets of radar scatter and wind elaborate the comparison: gravity waves of long period (τ > 6 h) are shown to modulate the scattering process.     

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.     

The Siple VLF transmitter as a multi-frequency probe of the Earth-ionosphere waveguide

1983 receptions of subionospheric signals radiated from Siple, Antarctica (L = 4.3) to neighboring stations Palmer (L = 2.3), Halley (L = 4.3), and South Pole (Λ = 74°), each ~ 1500 km from the horizontal (magnetically east-west) VLF transmitting antenna at Siple, were found to be strongly dependent upon azimuth and upon signal frequency. At Palmer, located equatorward in the broadside direction with respect to the antenna, signals near 2.5 kHz were often well defined, while the third harmonic of the transmitted signal, near 7.5 kHz, was not detected. Meanwhile, at Halley, the third harmonic was regularly observed and directionally stable, while the fundamental was often weak or undetectable. The field strength of the third harmonic component at Halley exceeded by ~ 40dB the level of the fundamental, when both were normalized to the same antenna input power. The large size of these effects is attributed in part to antenna properties that favor the endfire direction (toward Halley) at the 3d harmonic of the antenna half wave resonance frequency, and in general provide greater efficiency at higher frequencies. Other factors are high waveguide attenuation in the 2–4 kHz range and azimuth dependent differences in the propagating modes. The observed effects represent a way of extending the effective frequency range of the narrowband Siple antenna system, and also, by using the new crossed dipole configuration at Siple, of selectively probing certain regions of the Earth-ionosphere waveguide.     

Tidal winds at mesopause altitudes over Arecibo (18°N, 67°W), 5–11 April 1989 (AIDA '89)

An imaging Doppler interferometer (IDI) radar was operated during the three AIDA '89 campaigns in Puerto Rico over the period March–May of 1989. The output of the IDI analysis characterizes radar scattering in terms of a number of discrete ‘scattering points,’ also referred to as ‘multiple scattering centers,’ IDI/MSC for short. For each of these points the three-dimensional location, radial velocity and amplitude and phase are determined, similar to the output of meteor radars. We have applied the conventional Groves [(1959) J. atmos. terr. Phys. 16, 344–356] meteor wind radar analysis to the scattering points to produce the mean apparent motions over the height range from 70 to 110 km which are presented here. The mean apparent motion of the scattering centers is the quantity that would correspond to the neutral atmosphere wind or bulk motion if the scattering points are physical entities (such as turbulent eddies) whose motions are determined solely by advection. This is the quantity which is treated as the ‘wind’ in the analysis which follows and which should be compared to the wind measurements as deduced from the other methods employed during this campaign. There is, however, a caveat which supports the contention of Hineset al. [(1993) J. atmos. terr. Phys. 55, 241–287] that extreme care must be used in interpreting the velocities measured by partial reflection radars as winds. The current application of the Groves method of analysis has revealed motions from which one would infer a typical equatorial easterly circulation, with mean meridional circulation becoming significant only above 96 km. A periodogram analysis of the complete data interval (5–11 April) has shown the diurnal tide to be the most significant feature of the wind field at these altitudes, with zonal amplitudes up to some 50 m/s and meridional amplitudes approximately half this value. The 12 and 6 h tides become as significant as the diurnal above 100 km. The two day (48 ± 5 h) wave is the next most significant feature, with zonal amplitude increasing with height up to 30 m/s at 110km. The semidiurnal tide is not at all well developed below 100 km. However, analysis on a day by day basis reveals a significant semidiurnal component which is not phase coherent over the total interval. Mean vertical velocities are of the order of tens of centimeters per second and are considered to be more realistic than the meters per second velocities usually inferred from analyses of meteor trail drifts.     

Diurnal propagating tides in the low-latitude middle atmosphere

Using an equivalent gravity wave f-plane model it is shown that longitude variations in diurnal insolation absorption by tropospheric H₂O can account for longitudinal variations of at least ± 12–15% about zonal mean values in the diurnal wind amplitude at low latitudes (0–20°) between 80 and 100 km, by virtue of the non-migrating propagating tidal modes which are excited. Phase variations of about ± 0.75 h also occur. These percentage variations are conservative estimates, since the background migrating (1,1,1) mode appears to be slightly (20–25%) overestimated in amplitude. In addition, the assumed eddy dissipation values, which appear necessary to model the breaking (1,1,1) mode, are larger than generally considered ‘reasonable’ by photochemical modellers. For a photochemically more reasonable eddy diffusion profile, estimates of longitude differences in diurnal wind amplitude are quite similar to the above values below 87 km, but increase to ± 17–25% near 100 km, with accompanying phase variations of ± 1–2 h about zonal mean values. In addition, it is shown that radiative damping by CO₂ parameterized by a scale-dependent Newtonian cooling coefficient accounts for no more than a 20% reduction in the amplitudes of diurnal propagating tides above 80 km.     

The influence of geomagnetic activity on the upper mesosphere lower thermosphere in the auroral zone. II. Horizontal winds

A spaced antenna partial reflection radar located at Mawson, Antarctica (67°S, 63°E, invariant latitude 70°S), has been used to measure the horizontal wind field in the height range 70–110 km. Three years of data (1985–1987) from the radar have been analysed in order to investigate correlations between geomagnetic activity (determined from the local K-index) and the horizontal wind. Results are analysed using a randomization technique and show that larger winds are measured during geomagnetically active periods in both the raw (or unfiltered) wind values and in the medium-frequency (2–6 h period) and high-frequency (1–3 h period) components. The raw winds tend to be shifted towards the geographic NW to NE quadrant in the early morning hours during high K-times. The observed correlation is seen down to 86 km and shows a seasonal dependence. The mean r.m.s. velocity of the radar scatterers and the angular spread of the return echoes are also found to be correlated with geomagnetic activity. The medium- and high-frequency components of the wind are polarized in the magnetic zonal direction during all seasons of the year.     

Coherent-array HF Doppler sounding of traveling ionospheric disturbances: I. Basic technique

We present an introduction to the use of phase-coherent, multi-receiver HF Doppier sounding arrays for measuring the horizontal velocity of traveling ionospheric disturbances (TID's). The point of departure is the theorem of Pfister (1971, J. atmos. terr. Phys. 33, 999) relating ray Doppler to ray zenith angle for a monostatic full reflection sounder. Retaining the simple model of a specular, smooth ionospheric reflector which is deformed by a propagating undulation, we first generalize the theorem to bistatic sounding geometry and then include the effects of amplitude in addition to phase. Next, these results are cast into an algorithm for treating multi-receiver phase sounders containing many diverse baselines, in order to obtain an accurate and unambiguous solution in the plane of wave slowness (inverse of velocity). The point spread function of this solution is controlled by process bandwidth and by array geometry. We illustrate the coherent-array approach using data from an eight receiver array during passage of a TID.     

Studies of high latitude mesospheric turbulence by radar and rocket 1: energy deposition and wave structure

During early spring, 1985, the MAE-3 (Middle Atmospheric Electrodynamics) Program was conducted at Poker Flat Research Range, Alaska to study the origin of wintertime mesospheric echoes observed with the Poker Flat MST radar there, by probing the mesosphere with in situ rocket measurements when such echoes occurred. Pre-launch criteria required the appearance of echoes exhibiting some wave structure on the MST radar display; these could be met even under weak precipitation conditions with riometer absorption near or above 1.0 dB. Two morning rockets were launched under such conditions, the first (31.048) on 29 March 1985, at 1703 UT and the second (31.047) on 1 April 1985, at 1657 UT. Both payloads were deployed on a high altitude parachute near a 95 km apogee to provide a stable platform for data acquisition within the mesosphere (below 80 km). Each payload carried a solid state detector to measure energetic electrons between 0.1 and 1.0 MeV and an NaI crystal detector to measure x-rays from >5 to >80 keV. Payload 31.048 also carried a positive ion ‘turbulence’ probe which measured ion density changes (ΔN_i/N_i) during payload descent, whereas 31.047 carried a nose tip ‘turbulence’ probe designed to measure electron density changes (ΔN_e/N_e) during upleg ram conditions plus a Gerdien condenser for the measurement of bulk ion properties during downleg. The energy deposition curves for each event exhibited peak deposition rates between 75 and 80 km with a half width of 16–18 km, almost exclusively induced by precipitating relativistic electrons. They also showed a maximum bottomside gradient between 65 and 75 km. Radar echoes and atmospheric turbulence were observed in the same altitude domain, consistent with the anticipated need for adequate free thermal electron gradients to make such phenomena visible on the radar. The vertical wave structure from radar echoes was found to be consistent with that observed in horizontal wind and temperature profiles measured by Datasondes flown shortly after each large rocket. An analysis of the wave structure from radar data has shown that although large scale waves (λ_z ~ 7 km) were found to be present, a higher frequency shorter wavelength (∼ 1–3 km) component probably played a more significant role in modulating the signal-to-noise structure of the radar echoes.