A formulation is described for finding the surface vertical E field throughout the lateral extent of the spherical Earth-ionosphere waveguide at extremely low frequencies (ELF) below 100 kHz for a completely general model of values of the sine of the waveguide eigenangle. These values are specified at elements of a mesh point array that extends throughout the lateral extent of the spherical shell guide. The solution is based on the condition that the vertical E field at the surface must satisfy the reduced two-dimensional wave equation at every localized group of mesh points. It has been found practical to use a mesh size of 312.5 km halfway to the antipode, which implies a mesh size of $0 km at a distance 1 Mm from the transmitter. oeigensolutions are found for any far range from the transmitter for which there is no transverse variation of properties, and these are combined with the more general solution which applies to the region of transverse inhomogeneity. The latter may include the entire waveguide. An example is given for illustration. It is expected that use of the formulation will make it possible to describe the combined effects of inhomogeneities such as land-sea boundaries around continents and large seas and islands, the day-night terminator, and particle precipitation on the amplitude and phase of the vertical E field.
INTRODUCTION
For a number of years there has been interest in long-range extremely low frequency (ELF) communications [Wait, 1974]. The ELF waves travel within the spherical shell waveguide formed by the Earth as a lower boundary and by a diffuse upper boundary formed principally by the ionosphere E layer. The propagation of ELF waves in a homogeneous spherical shell waveguide has been described by various workers [e.g., Galejs, 1972]. Inhomogeneities in the real Earth-ionosphere waveguide, however, have a pronounced effect on propagation. The difference in the properties of the day and night hemispheres represents a large-scale inhomogeneity in the spherical shell waveguide [Large and Wait, 1968]. Inhomogeneities in the Earth-ionosphere waveguide also include differences in ground conductivity over continents, seawater, and ice caps. Anomalies in ELF propagation are also caused by particle precipitation in the polar regions [lmhof et al., 1978], and can be caused by resonance between a layer of sporadic E This paper is not subject to U.S. copyright. Published in 1989 by the American Geophysical Union. Paper number 88RS03707. and the lower edge of the E layer [Pappert and Moler Pappert, 1980Pappert, , 1985aPappert, , 1986.Wait [1964] introduces the use of the reduced two-dimensional wave equation for describing the effects of inhomogeneities in the Earth-ionosphere waveguide on single-mode VLF propagation. He also points out the optical analogy of the approach with propagation of light through a lens of refractive index equal to the sine of the waveguide eigenangle. This approach is used in the theory that follows. It is assumed that the propagation of ELF waves may be described by ...