The use of equivalent current systems in the interpretation of Geomagnetic Deep Sounding data

Summary. Many geomagnetic variation anomalies are probably caused by the channelling, through small-scale bodies, of electric currents induced in much larger conductors elsewhere. Consequently, the direct interpretation of anomalous magnetic fields by modelling the electromagnetic response of conductive structures may give misleading results. It is suggested that, rather than attempting to proceed directly from the electromagnetic fields to conductivity models, we should instead take the intermediate step of determining the distribution of anomalous current flow.
Maps of the anomalous fields over a conductive structure can be generated from inter-station transfer functions. If it is assumed that the internal currents are concentrated in a thin sheet at a specified depth, the equivalent current system in the sheet can be computed directly from the vertical magnetic field. The most straightforward method of performing this calculation is to compute the Fast Fourier Transform of the magnetic field data, and then to apply a wavenumber filter.
The presence of any vertical currents invalidates the thin sheet model. However, if the spatial distribution of a horizontal component of the anomalous magnetic field is also known, the presence of any vertical currents can be detected directly, and their position determined. The value of the methods is illustrated by applying them to the interpretation of a Geomagnetic Deep Sounding survey of the Kenya rift valley.     

The mapping of induced currents around the Kenya Rift: a comparison of techniques

Summary. Single-station and inter-station transfer functions are derived from data recorded by a magnetometer array in and around the Kenya Rift Valley. Different methods of presentation of the transfer function estimates are compared in terms of their ability to define lateral variations in electrical conductivity. When, as in East Africa, the conductivity structure is complex and three-dimensional, by far the most useful method of presentation is to use the transfer functions to simulate the anomalous internal fields associated with regional current flow at a particular azimuth. The use of inter-station rather than single-station transfer functions is almost essential where the amplitudes of anomalous horizontal fields are of the same size or greater than the normal fields. Horizontal field transfer functions prove to be just as useful as those usually calculated for the vertical component, and provide strong additional constraints on the internal current system.
Maps of simulated vertical and horizontal fields of internal origin indicate the importance of current channelling around the Rift Valley. A conductor just to the east of Nairobi apparently funnels regionally induced currents into two conductors associated with the Rift and dome; one at shallow depths beneath the Rift floor, and a deeper body to the east.     

The frequency response of the horizontal magnetic field for a conductive channel

Summary. Babour & Mosnier's results showing no frequency dependence between the anomalous horizontal magnetic field above a conductor and the difference between the horizontal magnetic fields above and below the conductor over a wide range of frequencies led them to conclude that this effect is due to current channelling. A two-dimensional numerical model of a conductive channel, with a uniform horizontal source field, shows the same effect over a wide range of frequencies. Thus local induction can show the same effect.     

E-polarization induction in two thin half-sheets

Summary. An analytical solution is obtained for the E-polarization problem of electromagnetic induction in two adjacent half-sheets underlain by a uniform conducting half-space. In this mode the inducing magnetic field is assumed horizontal, uniform and perpendicular to the discontinuity. The same model was previously solved under B-polarization by Dawson & Weaver. The present solution then completes the study of two-dimensional induction in the described model. Further, it extends both the analytic E-polarization solution of Weidelt by the inclusion of an underlying conductor and that of Raval, Weaver & Dawson by the inclusion of arbitrary conductance values for the two surface sheets. The solution may be used as an idealized model of the coast effect and allows detailed study of the field behaviour near the discontinuity. The horizontal magnetic field on each side of the surface layer has a finite jump discontinuity at the interface and the vertical magnetic field exhibits a logarithmic singularity there. If the right-hand conductance (say) becomes infinite, the horizontal magnetic field exhibits an algebraic singularity as the coastline is approached from the right, while the vertical magnetic field does likewise from the left. Calculations are presented for the same two models as discussed in B-polarization by Dawson & Weaver and the results are compared to values obtained from a more general numerical scheme. The electric current distribution inside the conducting half-space is depicted for the second model.     

The use of the vertical component in geomagnetic depth sounding

In previous papers, we showed how we found it possible to define the anomalous field by means of its horizontal components defined as the difference between the field recorded at a reference station and that recorded at stations on the anomaly. In the present short note, we extend this definition to the vertical component. We compare the previous results with those obtained by characterizing the anomaly by the total vertical component and by the differential vertical component. This comparison leads us to conclude that the differential components are more effective than the total vertical component for delineating an anomaly.     

Three-dimensional model for tides and surges with vertical eddy viscosity prescribed in two layers—I. Mathematical formulation

Summary. The mathematical basis of a spectral method for the numerical solution of the three-dimensional hydrodynamic equations for tides and surges is described. Vertical eddy viscosity is prescribed in two layers, upper and lower. In each layer, horizontal components of current are expanded through the vertical in terms of a set of eigenfunctions. Coefficients of these expansions are evaluated in the horizontal and through time using a two-dimensional numerical time-stepping procedure. Thence the three-dimensional current structure is determined.     

Transient electromagnetic fields about an infinitesimally long grounded horizontal electric dipole on the surface of a uniform half-space

Summary. Laplace and Bessel Transforms are used to solve for the transient behaviour of the electromagnetic fields after switching off a steady current in a grounded infinitesimal horizontal dipole on the surface of a uniformly conducting half-space. Simple analytic expressions, which are valid for times sufficiently long after the switch that displacement terms can be ignored, are obtained on the surface of the half-space for the electric field and the time derivative of the magnetic field. At the instant of switching an infinitesimally long image becomes established directly under the source dipole. It is the diffusion of this image which gives the vertical magnetic field and horizontal electric fields their transient behaviour. During the transient, there is also a decaying charge distribution on the surface.     

Telluric measurements and differential geomagnetic soundings in the Rhinegraben (period range: 1–125 s)

Summary. In 1976, seven stations measuring the variations of the telluric and geomagnetic fields in the period range 1–125 s were operated in the southern part of the Rhinegraben. The study of the recordings shows that the telluric field is linearly polarized according to a direction perpendicular to that of the horizontal anomalous magnetic fields and that telluric and anomalous magnetic fields have the same time dependence. The conducted currents responsible for the anomaly flow probably into the superficial conductive layer.     

Surface-wave propagation in an ocean basin with an anisotropic upper mantle: numerical modelling

Summary. The characteristics of surface-wave propagation in ocean basins are examined numerically for models with two types of anisotropic alignment in the upper mantle: one resulting from glide-plane slip in olivine with horizontal or vertical slip-planes, and the other from syntectonic recrystallization of olivine in a zone of horizontal shear. Glide-plane slip can cause highly anomalous inclined-Rayleigh particle-motion in the third-generalized mode (corresponding to the isotropic second-Rayleigh mode). The amplitude of this anomaly is rather insensitive to details of the structure. Syntectonic recrystallization can cause an anomalous combination of inclined-and tilted-Rayleigh motion in all modes. The variation with period of the amplitude of the anomaly in the fundamental mode can indicate the approximate depth to the anisotropic layer. In both types of alignment, the sense of tilt and the inclination varies with direction of propagation in a manner characteristic of the structural symmetry.     

Bimodal electromagnetic induction in non-uniform thin sheets with an application to the northern Pyrenean induction anomaly

In order to handle the distortion of large-scale induced electric currents by local conductivity anomalies, the problem of electromagnetic induction in non-uniform thin sheets has been reformulated in terms of an integral equation over the anomalous domain. This formulation considers in the layered substratum in addition to toroidal currents also the poloidal current mode (vertical current loops), at the expense that two scalar functions have to be determined. Simple formulas for the required kernels are derived. The algorithm is applied to model the gross features of the northern Pyrenean induction anomaly. It is suggested that this pronounced anomaly results from a conductive channel between the Atlantic Ocean and the Mediterranean Sea.     

Electromagnetic response of an ocean-coast model to E-polarization induction

summary . An ocean-coast model which consists of a uniformly conducting half-space screened by a perfectly conducting half-plane (the model ocean) is studied. On the land the electric field decreases continuously to zero as the coast is approached. The horizontal magnetic field component is found to vary rapidly, but remains finite; the vertical component on the other hand, increases to infinity at the coast. On the surface of the model ocean as well as on the sea floor, electric field and vertical magnetic field are both nil, but the horizontal magnetic field becomes singular as the seashore is approached. This horizontal magnetic field however, is different on the sea floor and at the ocean surface, because the integrated ocean current is finite, even growing to infinity as the shore is approached. The very large ocean currents near the shore act as an extremely long line antenna, which radiates far afield. This antenna feature explains the very long range of the ocean-coast effects observed under E -polarization induction, compared to the corresponding H -polarization effects where no such antenna-like feature occurs. A similarly large difference of ranges can be expected for all shallow structures with large lateral conductivity contrasts. The present study may therefore be of some interest in relation to geomagnetic depth soundings by the inductive and magnetotelluric methods, as well as in understanding the ocean-coast effect known for some time from records of coastal observatories.     

Significance of the sign changing of the imaginary arrows in geomagnetic induction investigation

Summary. In this investigation, we carry out a two-dimensional study of the dependence of the imaginary Parkinson arrows on the frequency of the inducing geomagnetic field. Our results demonstrate that the imaginary arrows reverse direction as the inducing period varies. Therefore, we consider that there is no way to fix a consistent sign convention for the imaginary arrows even when the time factor is taken into account. We find that in the twodimensional case there exists a characteristic period T _c at which the phase difference between the vertical and horizontal magnetic components is zero. It is anticipated that T _c is related to the parameters of the conductivity anomaly and the status of the half-space host.     

2-D or 3-D interpretation of conductivity anomalies: example of the Rhine-Graben conductivity anomaly

Summary. Evidence of a conductivity anomaly in the Rhine-Graben was first given about 15 years ago and consequently led to the definition of various models of induction in the region for periods ranging from a few minutes to a few hours. These models reflect two antagonistic ways of explaining the observed anomalous variations of the magnetic field: direct induction in a two-dimensional (2-D) structure or static distortion of telluric currents by the resistive crystalline Vosges (France) and Schwarzwalde (Germany) massifs. We discuss the two approaches using a simple formalism. In particular, we show that the self-induction related to the anomalous currents flowing in the Rhine-Graben is negligible for periods larger than 1000 s, and that, even though the static distortion of telluric currents does account for the observed anomaly, 2-D models can explain some of its features. We also show how the channelled currents are induced in the large sedimentary basins surrounding the area under study.
An experimental verification of this result is given.     

Modelling of solar quiet magnetic field variations near a conductivity anomaly

Summary. A simplified model of the solar quiet-time ionospheric current system is used to calculate the induced currents in a model earth. The conductivity is assumed to be constant below a depth of about 400 km and zero above that depth. The current induced in the north—south conductivity anomaly under the Rocky Mountains is then estimated from the time-varying potential difference between points at 30 and 45° latitude at the surface of the conducting sphere. The purpose of these calculations is to investigate whether variations in the latitude of the northern hemisphere current system vortex will substantially alter the relationship between the observed magnetic field components at the Earth's surface and the local magnetic field gradient caused by the conductivity anomaly. We find that a 10° shift in the latitude of the ionospheric current focus causes a change of 6 per cent or less in the transfer function from the field components to the gradient in the total field. Thus such latitude shifts cannot explain much of the magnetic field gradient variation at periods near 24 hr that has been observed near Boulder, Colorado.     

Impulsive Flux Transfer Events and Solar Flares

Experimental and theoretical aspects of the fundamental process of magnetic flux transfer ('field line reconnection') are reviewed. Explosively rapid events are observed in a laboratory magnetoplasma device designed to test certain aspects of neutral point theories of the solar flare mechanism. In these events, ultra-fast release of magnetic energy through reconnection of field lines is effected by a conduction mode instability in which an abrupt upward transition to anomalous resistivity takes place in the neighbourhood of the X-point. The events are triggered when the X-point current density exceeds the conduction mode instability threshold. The energy that is released is that which is associated with the induced current system that supports field line reconnection under ordinary quiescent conditions. Many characteristics of the laboratory events and their associated phenomena, when appropriately scaled, show remarkable agreement with corresponding flare observational data. Some new directions for further development of flare theory are therefore recommended. These experiments indicate that a number of assumptions commonly made in theoretical analysis of the reconnection process are inappropriate. The flux transfer process can be well understood theoretically in terms of a simple circuit analogue which models the laboratory process.     

Two-dimensional induction in a thin sheet of variable integrated conductivity at the surface of a uniformly conducting earth

Summary. The forward solution of the general two-dimensional problem of induction in a model earth comprising a uniformly conducting half-space covered by a thin sheet of variable integrated conductivity is obtained. Unlike some previous treatments of similar problems, the method presented here does not require the field to be separated into its normal and anomalous parts. Both the E - and B -polarization modes of induction are considered and in each case the solution is expressed in terms of the horizontal component of the electric field satisfying, on the surface of the conductor, a singular integral equation whose kernel is a well-known analytic function. A recently published solution of the coast effect is included as a special case. The numerical procedure for solving the integral equations is described and some illustrative calculations are presented.     

Forward modelling of direct current and low-frequency electromagnetic fields using integral equations

We present a semi-analytical, unifying approach for modelling the electromagnetic response of 3-D bodies excited by low-frequency electric and magnetic sources. We write the electric and magnetic fields in terms of power series of angular frequency, and show that to obey Maxwell's equations, the fields must be real when the exponent is even, and imaginary when it is odd. This leads to the result that the scattering equations for direct current fields and for fields proportional to frequency can both be explicitly formulated using a single, real dyadic Green's function. Although the underground current flow in each case is due to different physical phenomena, the interaction of the scattering currents is of the same type in both cases. This implies that direct current resistivity, magnetometric resistivity and electric and magnetic measurements at low induction numbers can all be modelled in parallel using basically the same algorithm. We make a systematic derivation of the quantities required and show that for these cases they can all be expressed analytically. The problem is finally formulated as the solution of a system of linear equations. The matrix of the system is real and does not depend on the type of source or receiver. We present modelling results for different arrays and apply the algorithm to the interpretation of field data. We assume the standard dipoledipole resistivity array for the direct current case, and vertical and horizontal magnetic dipoles for induction measurements. In the case of magnetometric resistivity we introduce a moving array composed of an electric dipole and a directional magnetometer. The array has multiple separations for depth discrimination and can operate in two modes. The mode where the predominant current flow runs along the profile is called MMR-TM. This mode is more sensitive to lateral variations in resistivity than its counterpart, MMR-TE, where the mode of conduction is predominantly perpendicular to the profile.     

Electromagnetic induction fields in the deep ocean north-east of Hawaii: implications for mantle conductivity and source fields

Summary. We have analysed a thirty-six day recording of the natural electric and magnetic field variations obtained on the deep ocean floor north-east of Hawaii. The electromagnetic fields are dominated by tides which have an appreciable oceanic component, especially in the east electric and north magnetic components. The techniques of data analysis included singular value decomposition (SVD) to remove uncorrelated noise. There are three degrees of freedom in the data set for periods longer than five hours, indicating a correlation of the vertical magnetic field and the horizontal components, suggesting source field inhomogeneity. Tensor response functions were calculated using spectral band averaging with both SVD and least squares techniques and rotated to the principal direction. One diagonal component, determined mainly by the north electric and east magnetic fields, is not interpretable as a one-dimensional induction phenomenon. The other diagonal term of the response function indicates a rapid rise in conductivity to 0.05 mho m⁻¹ near 160 km. No decrease in conductivity below this depth is resolvable. Polarization analysis of the magnetic field indicates moving source fields with a wavelength near 5000 km. Model studies suggest that the two dimensionality in the response function may be caused by motion in the ionospheric current system.     

Three-dimensional induction in a non-uniform thin sheet at the surface of a uniformly conducting earth

Summary. A new method for solving problems in three-dimensional electromagnetic induction in which the Earth is represented by a uniformly conducting half-space overlain by a surface layer of variable conductance is presented. Unlike previous treatments of this type of problem the method does not require the fields to be separated into their normal and anomalous parts, nor is it necessary to assume that the anomalous region is surrounded by a uniform structure; the model may approach either an E- or a B -polarization configuration at infinity. The solution is expressed as a vector integral equation in the horizontal electric field at the surface. The kernel of the integral is a Green's tensor which is expressed in terms of elementary functions that are independent of the conductance. The method is applied to an illustrative model representing an island near a bent coastline which extends to infinity in perpendicular directions.     

A study of two polar magnetic substorms with a two-dimensional magnetometer array

Summary. The paper reports studies of the three-dimensional magnetospheric—ionospheric current systems which produced polar magnetic substorms on 1974 September 7 and September 18. The data were magnetic perturbation fields observed with a two-dimensional array of 23 three-component magnetometers located in western Canada beneath the auroral oval. In an earlier study of a substorm of September 11 (Bannister & Gough) the fields fitted calculated field for a Boström Type 1 current loop with field-aligned currents at east and west ends of the ionospheric segment, and with uniform current density across the width. The substorms here reported could not be modelled with uniform current density. An inverse method due to Oldenburg was therefore used to estimate current density distributions, and satisfactory fits of calculated to observed field resulted. Each substorm was modelled at six representative epochs. In general the principal ionospheric current seem by the array was westward. At four epochs of the September 7 substorm and throughout the September 18 substorm, significant eastward ionospheric current (or its equivalent in terms of the fields produced) was observed north of the westward electrojet. Northwestward bends in the ionospheric current segments were found at four epochs on September 7 and at three epochs on September 18. As in the September 11 substorm (Paper 1), these bends were either west of or close to magnetic midnight. In some cases the bends may follow the auroral oval, but in others they are sharper and may be associated with the Harang discontinuity. East of geomagnetic the ionospheric currents tend to run in a constant geomagnetic midnight latitude range. The developments of the three substorms, of September 7, 11 (Paper 1) and 18, are compared. They showed a variety of shifts in longitude, though all moved eastward relative to magnetic midnight.