Wave form and spectral distribution of the electromagnetic field of lightning currents

In a resonant wave guide model of lightning currents, two impulse type standing waves can exist: aperiodic waves of the Bruce-Golde form (type 1), and damped oscillations (type 2). The electromagnetic waves generated by these two types of lightning currents are calculated for various distances and compared with observations. It is shown that the measured wave forms of sferics at distances smaller than about 300 km of return strokes (R strokes) as well as of intracloud strokes (K strokes) are generated mainly by type 2 lightning currents. The channel parameters like channel length and channel diameter derived from the observed sferics are 19 km and 4.6 cm, respectively, for the average R stroke, and 4 km and 1.6 cm for typical K strokes. The large values of the lengths probably correspond to the real lengths of oblique and tortuous channels, rather than to their vertical elevation. The finite electric conductivity of the earth modifies the high frequency component of the wave forms. With decreasing conductivity and/or increasing distance, the rise times of the radiation component to its first maximum increases and the maximum amplitude decreases.Typical rise times for R-strokes are about 3 μs consistent with the observations if the electric conductivity of the earth is of the order of 10⁻³ S/m. The spectral functions of the wave forms are also calculated. The spectral amplitude of the average type 2 R-stroke has its maximum near 4 kHz, and that of the type 2 K-stroke maximizes near 35 kHz. Within the high frequency region at frequencies greater than about 300 kHz. the spectral amplitudes decay proportional to the reciprocal third power of the frequency. The radiation component in the far field contains 7% (16%) of the total electromagnetic energy of type 2 R (K) strokes.     

Temporal behaviour of lightning HF radiation at 3 MHz near the time of first return strokes

Simultaneous recordings of the electric field and high frequency emissions (HF) at 3 MHz at times close to the first return stroke are presented. It is shown that there is no time delay between the beginning of the return stroke electric field and the build up of the HF. The delay observed in previous studies may be due to propagation effects. The results in turn suggest that the first return stroke from the very beginning is a strong source of HF at 3 MHz.     

Phase velocity studies of 34-cm E-region irregularities observed at Millstone Hill

Phase velocity observations at E-region heights made with the Millstone Hill 440 MHz radar find no evidence of an ion acoustic limiting speed for phase speeds observed near 0° magnetic aspect angle. Under most circumstances the phase speed increases steadily with increasing backscattered power amplitude. For a 34cm volume backscatter cross-section, σ_v, less than ∼5 × 10⁻¹³ m⁻¹, the phase speed is at or below the usual ion acoustic speed in the E-region (350m/s), and increases only slowly with the observed backscattered power amplitude (∼50 m/s per 10dB). At higher power levels, the phase speed exceeds 350 m/s, reaching values in excess of 750 m/s at times, and increases more rapidly with backscattered power (∼200 m/s per 10dB). Phase velocity/time maps observed over a 3° span of latitude suggest that many features of the phase speeds observed are directly related to changes in the ambient convection electric field in the E-region due to changing activity conditions or the effects of superimposed magnetospheric pulsations.     

Auroral radar observations at Siple Station,Antarctica

A dual coherent 50 MHz auroral radar installed at Siple Station, Antarctica, continuously monitored echo intensity and mean-Doppler velocity from 29 December 1976 until late March 1978. We describe here the experimental technique and present some statistical results including yearly averaged echooccurrence patterns and irregularity drift velocity characteristics. Our results show that:

• 1.(1) the irregularity drifts are approximately westward (poleward electric fields) in the evening, and eastward (equatorward electric fields) in the morning following the electric field reversal in the region of the Harang discontinuity;
• 2.(2) the Harang discontinuity under disturbed conditions (average K_p = 5₀) as seen by both radars is located around 2100–2300 MLT;
• 3.(3) the relationship between the irregularity drift velocity and the actual electron drift velocity is strongly dependent on the angle between the radar beam and the earth's magnetic field, as predicted by linear theory.

     

Broad-band auroral VLF hiss and inverted-V electron precipitation in the polar magnetosphere

A polar map of the occurrence rate of broad-band auroral VLF hiss in the topside ionosphere was made by a criterion of simultaneous intensity increases more than 5 dB above the quiet level at 5, 8, 16 and 20 kHz bands, using narrow-band intensity data processed from VLF electric field (50 Hz–30 kHz) tapes of 347 ISIS passes received at Syowa Station, Antarctica, between June 1976 and January 1983.The low-latitude contour of occurrence rate of 0.3 is approximately symmetric with respect to the 10–22 MLT (geomagnetic local time) meridian. It lies at 74° around 10 MLT, and extends down to 67° around 22 MLT. The high-latitude contour of 0.3 lies at invariant latitude of about 82° for all geomagnetic local times. The polar occurrence map of broad-band auroral VLF hiss is qualitatively similar to that of inverted-V electron precipitation observed by Atmospheric Explorer.(AE-D) (Huffman and Lin, 1981, American Geophys. Union, Geophysics Monograph, No. 25, p. 80), especially concerning the low-latitude boundary and axial symmetry of the 10–22 h MLT meridian.The frequency range of the broad-band auroral VLF hiss is discussed in terms of whistler Aode Cerenkov radiation by inverted-V electrons (1–30 keV) precipitated from the boundary plasma sheet. High-frequency components, above 12 kHz of whistler mode Cerenkov radiation from inverted-V electrons with energy below 40 keV, may be generated at altitudes below 3200 km along geomagnetic field lines at invariant latitudes between 70 and 77°. Low-frequency components below 2 kHz may be generated over a wide region at altitudes below 6400 km along the same field lines. Thus, the frequency range of the downgoing broad-band auroral hiss seems to be explained by the whistler mode Cerenkov radiation generated from inverted-V electrons at geocentric distances below about 2 R_E (Earth's radius) along polar geomagnetic field lines of invariant latitude from 70 to 77°, since the whistler mode condition for all frequencies above 1 kHz of the downgoing hiss is not satisfied at geocentric distance of 3 r_e on the same field lines.     

Electric fields and currents in the equatorial electrojet deduced from VHF radar observations—III. Comparison of observed ΔH values with those estimated from measured electric fields

From VHF backscatter radar measurements at Thumba (dip: 56′S) of the phase velocities of type II irregularities in the equatorial electrojet (EEJ), electric field (E_y) values are estimated for different times of the day. Using the electric field values thus deduced and the Pedersen and Hall conductivities calculated using model values of electron densities and the collision frequencies of ions and electrons, the height integrated current intensity in the EEJ is estimated. The surface level geomagnetic field perturbation ΔH produced by this ionospheric current is then calculated. The calculated values of ΔH are compared with observed values of ΔH (after subtracting the magnetospheric contribution of D_st) for a number of days. The comparisons show good agreement between observed ΔH values and those calculated from measured electric fields. The agreement is found to be good even when type I irregularities are present at higher altitudes in the EEJ. This comparative study demonstrates the validity of estimating electric field values from VHF radar measurements and it indicates the possibility of deducing electric field values from ground level ΔH values, at least for statistical studies.     

Comparison of stratospheric ozone profiles retrieved from microwave-radiometer and Dobson-spectrometer data

Passive Microwave Remote Sensors (MRS) can provide information on the composition of the atmosphere by measuring the thermal radiation emitted from rotational transitions of atmospheric molecules Ozone profiles simultaneously obtained from an MRS and a Dobson instrument (‘Umkehr’ method) are compared over a time period of approximately four months. The microwave measurements yield ozone concentrations which are 20–30% higher than the ‘Umkehr’ values. A critical, though not well known, parameter for the microwave inversion procedure is the foreign gas pressure broadening parameter (C) for the observed 142GHz ozone resonance. Throughout the intercomparison we used a value of 3 MHz mb⁻¹. There is recent theoretical and experimental indication that C is more likel y to be as low as 2.5 MHz mb⁻¹. If we use this new value for C all microwave retrieved profiles decrease by 20–25%, thus leading to a far better agreement with the ‘Umkehr’ results. Our measurements therefore strongly support the proposed value of 2.5 MHz mb⁻¹. A final answer on MRS ‘Umkehr’ correlation accuracy cannot be given. We feel that comparison on a day-to-day basis may be rather meaningless and monthly mean values should be used. On the other hand, there was relatively little change in these mean values during the intercomparison period The MRS showed its potential to retrieve ozone profiles also under adverse meteorological conditions, such as cloud cover or fog.     

Incoherent scatter observations of thin ionization layers at Sondrestrom

High resolution incoherent-scatter observations of E-region thin (1–3 km) metallic ion layers are presented. Data were collected during three different periods from August 1990 to August 1991, in three different experimental modes. First, the antenna was directed vertically and the entire duty cycle was devoted to Barker coded multi-pulse [Zamlutti (1980) J. atmos. terr. Phys.42, 975–982] measurements to determine the densities and temperatures in the E-region with 300 m resolution. The second experiment measured the F-region electric field as well as the high resolution E-region densities. For the third experiment the antenna was scanned magnetic north-south while only the E-region densities were measured. The experiments were carried out on 16 different nights for a period of 4 h each night at a time near magnetic midnight. Thin ionization layers were observed on 12 of the 16 nights. The first experiment demonstrated that the thin layers are composed of a significant fraction of heavy metallic ions; assuming the layers are composed of a mixture of Fe⁺ and Mg⁺ a composition estimate of 63% Fe⁺ was obtained in one example. The second experiment investigated the relationship between the direction of the electric field and the presence of the thin layers. In these observations thin layers were only present when the electric field was pointed in the magnetic north-west or south-west quadrants, most frequently when the field was near magnetic west. Correlation between layer altitude and field direction was also observed, layers occurring at higher altitudes for fields directed in the north-west, and lower altitudes for fields directed to the south-west. The observations are compatible with the electric field mechanism for thin ionization layer formation. The scanning experiment showed that the layers were of a limited latitudinal extent, typically about 100 km up to a maximum of about 200 km.     

Meridional electric currents in the equatorial F-region

A study of the boundary conditions for the equatorial thermospheric transport equations by the authors has led to the theoretical prediction of the vertical electric field at the base of the F-region. Earlier, this result was applied to the calculation of the zonal wind field in the equatorial F-region. In this work, the aforementioned model is applied to the calculation of the F-region electric current field in the meridional plane as a function of time and the east-west magnetic field generated by these currents. In particular, the field at sunset is compared with the observations made by Magsat.     

Lightning transient fields in the atmosphere—low ionosphere

The complexity of the electric properties of the atmosphere, which influences the excitation and propagation of electromagnetic signals generated by lightning discharges, gives rise to a diversity of the typical signal wave forms and amplitudes in the various space-time domains. This causes difficulties and controversies in the interpretation of the experimental data. In geophysical problems of such a kind, considerations accounting for a broadband spectrum are particularly important.In the present paper, lightning induced ULF-VLF fields in the daytime inhomogeneous atmosphere and D-region of the ionosphere above the thundercloud are analyzed for a standard model of the vertical return stroke. The developed analytical approach makes it possible to unite the canonical model of the electromagnetic radiation into a free-half-space, the discharge (ELF-VLF spectrum) and ‘post-discharge’ (ULF spectrum) low-frequency models and to clarify the temporal scales in the mechanism of lightning coupling with the atmosphere and low ionosphere with respect to the local relaxation time, discharge characteristic time scales and observer location.     

Midlatitude measurements of the ionospheric electric field during the ALADDIN programme

Three measurements of ionospheric electric field were made during the 24 h ALADDIN rocket programme at Wallops Island (37°50′N, 75°29′W) on June 29–30, 1974. The first of these used a double probe instrument, flown at 1500 Local Solar Time, and the second and third measurements were made by barium cloud releases at evening and morning twilight. These three electric field vectors have been compared with the predictions of a number of models of electric field due to the dynamo effects of various atmospheric tides, and also of a possible magnetospheric origin. On the assumption that the measurements were made at a location equatorward of the afternoon convergence and poleward of the morning divergence in the electric field patterns related to the Sq current cystem, Stening's model of the diurnal variation of the electric field induced by the (1, −2) tidal model at the time of the Summer solstice correctly predicts the directions of the observed electric field. Forbes and Lindzen's model, incorporating the three major propagating tidal modes as well as the evanescent (1, −2) mode, also bears an acceptable relationship to the ALADDIN electric field directions. The ALADDIN E-field magnitudes are comparable with those obtained by ground-based observations (incoherent scatter) from Millstone Hill and from Saint Santin but are about half of Stening's model values, and three times those of Forbes and Lindzen.While the Millstone Hill E-field directions are compatible with the ALADDIN observations, Saint Santin E-field directions, at the same latitude but 75° difference in longitude, are distinctly different from ALADDIN, implying that longitudinal differences are significant.     

Effect of the large-scale convection electric field structure on the formation of thin ionization layers at high latitudes

A three-dimensional simulation of the high-latitude ionosphere was applied to investigate the geographical distribution of E-region thin ionization layers which may be formed by the action of the convection electric field. The simulation model computes the ion densities (O⁺, O⁺₂, N⁺, N⁺₂, NO⁺, Fe⁺), and temperatures as a function of altitude, latitude, and longitude. The stationary state momentum and continuity equations are solved for each ion species, then the energy equation is solved for electrons, neutrals, and a generic ion having the mean ion mass and velocity. The various electric field patterns of the Heppner and Maynard [(1987) J. geophys. Res.92, 4467–4489] convection electric field model were applied and the ionization density pattern was examined after a time sufficient for the formation of thin layers (≈2000 s). It was found that large areas of thin ionization layers were formed for each of the electric field patterns examined. Southward IMF B_z conditions resulted in thin layers forming in the pre-midnight sector in the latitude range north of about 70° to about 80°, and after midnight between 60 and 70°. For northward B_z conditions, the layers were mainly in the pre-midnight sector and covered a latitude range from about 60 to 80°.     

On the structure and circulation of the polar thermosphere under magnetically quiet conditions

A theoretical study is presented bearing on the thermospheric circulation and composition at polar latitudes. The observed motions and density perturbations in N₂, O and He have signatures which may be understood in terms of two different source mechanisms. We consider electric field momentum coupling and Joule heating as well as interactions between both processes. A spectral model in terms of vector spherical harmonics (with magnetic coordinates) is used, delineating the diurnal and mean (time independent) components. The important non-linearities are evaluated in configuration space. The electric field model of Volland and the global average density and temperature variations of Hedin (MSIS) are adopted as input. Our analysis leads to the following conclusions. (1) The vortex type double cell polar circulation (zonal wave number m = 1) is primarily driven by collisional momentum transfer from electric field induced ion convection. (2) Because of the thermospheric low pass filter, a large time independent component (zonal wave number m = 0) is produced by Joule heating; the heavier species (N₂) being concentrated where the lighter ones (O, He) are depleted, and vice versa due to wind induced diffusion (3) The electric field driven vortex circulation redistributes the mass and energy in the time independent density and temperature variations (from Joule heating), producing primarily diurnal variations (m = 1) in the temperature and composition near the pole and at auroral latitudes, again the heavier and lighter species varying out of phase. The above results are in substantial agreement with observations. It is worth noting that momentum rectification associated with the diurnally varying electric field and conductivity induces a weak zonally symmetric (single cell) prograde polar vortex. However, this motion is partially compensated by a retrograde vortex from geostrophic balance due to Joule heating, which dominates near the pole. These motions are small compared with the diurnally varying component in the polar circulation.     

Accelerated drying of water-wetted materials in electric fields

The electric currents, I, associated with drying rates of discs of paper toweling and cotton cloth wetted with water have been measured as functions of several parameters: E, the strength of the electric field in which the discs were placed; T, the temperature; H, the percent relative humidity. Drying rates increased monotonically with E, and typical drying times decreased by factors of up to about 10 as E rose from 0 to above 7 kV/cm. Measured electric currents over this range of E increased 3–4 orders of magnitude. A solid dielectric object placed in the electric field caused I to be reduced proportionately to the area of the disc that was blocked by the object.     

The dynamic F2-layer over Arecibo

The idealized ‘servo’ model of the ionospheric F2-layer, developed by Rishbeth, ganguly and Walker (1978), is used to simulate the observed behaviour of the daytime F2-peak at Arecibo for sunspot minimum. Taking the east-west electric field to be given by the observed plasma drift velocity perpendicular to the magnetic field, the theoretical equations are integrated using a trial-and-error approach to match the observed values of field-parallel plasma velocity, and the height and electron density of the F2-peak. From the calculation is determined empirically the meridional pressure-gradient force associated with the meridional neutral-air wind. The local time variation during the day is found to be consistent with the semidiurnal variation given by the MSIS atmospheric model of Hedinet al. (1977a, b), though with a phase shift that varies with season; on some days in the fall the pressure-gradient force displays a strong equatorward ‘surge’ in the evening. The values of F2-layer loss and diffusion coefficients needed to match the data are broadly consistent with the MSIS model. The analysis thus validates the MSIS model by way of ionospheric parameters quite independent of the data from which MSIS was originally derived.     

The E-region Rocket/Radar Instability Study (ERRRIS): scientific objectives and campaign overview

The E-region Rocket/Radar Instability Study (Project ERRRIS) investigated in detail the plasma instabilities in the low altitude (E-region) auroral ionosphere and the sources of free energy that drive these waves. Three independent sets of experiments were launched on NASA sounding rockets from Esrange, Sweden, in 1988 and 1989, attaining apogees of 124, 129 and 176km. The lower apogee rockets were flown into the unstable auroral electrojet and encountered intense two-stream waves driven by d.c. electric fields that ranged from 35 to 115 mV/m. The higher apogee rocket returned fields and particle data from an active auroral arc, yet observed a remarkably quiescent electrojet region as the weak d.c. electric fields (~ 10–15 mV/m) there were below the threshold required to excite two-stream waves. The rocket instrumentation included electric field instruments (d.c. and wave), plasma density fluctuation (δn/n) receivers, d.c. fluxgate magnetometers, energetic particle detectors (ions and electrons), ion drift meters, and swept Langmuir probes to determine absolute plasma density and temperature. The wave experiments included spatially separated sensors to provide wave vector and phase velocity information. All three rockets were flown in conjunction with radar backscatter measurements taken by the 50MHz CUPRI system, which was the primary tool used to determine the launch conditions. Two of the rockets were flown in conjunction with plasma drift, density, and temperature measurements taken by the EISCAT incoherent scattar radar. The STARE radar also made measurements during this campaign. This paper describes the scientific objectives of these rocket/radar experiments, provides a summary of the geophysical conditions during each launch, and gives an overview of the principal rocket and radar observations.     

Observation of low energy particle precipitation at low altitude in the equatorial zone

Precipitation of protons (~ 1 MeV) in the equatorial zone was investigated by the Phoenix-1 experiment on board the S81-1 mission from May–November, 1982. The protons show a precipitation peak along the line of minimum magnetic field strength with a full width at half maximum (FWHM) of 13°. The index of equatorial pitch angle distribution is q ~ 19. The peak proton flux shows a fifth-power altitude dependence, and the proton flux shows approximately a factor of 3 times increase in 1982 compared to that in 1969 due, possibly, to the stronger (~ 1.2 times) solar maximum conditions of 10.7cm radiation in 1982.     

Power spectra of VHF intensity scintillations from F2- and E-region ionospheric irregularities

An experiment is described for the routine study of scintillations and ionospheric irregularities at high-latitudes using NNSS satellites with additional coordinated observations by means of the EISCAT ionospheric radar facility. Early results, obtained during the development phase of the experiment, are presented of the power spectra of intensity fluctuations at 150 MHz observed at the equatorwards edge of the high-latitude irregularity zone. The spectra of 165 samples of night-time scintillation recorded during October 1982 to May 1983 show a spectral index with a mean value of −3.58 and a steepening of the spectral slope with increasing S₄. Some examples of scintillation arising from irregularities at E-layer height show spectral indices of magnitude generally smaller than for F-region cases. A few spectra have been found with a clear break in spectral slope at around 10 Hz, suggesting two regimes for irregularities of different scale sizes.