On crosswell diffusive time-domain electromagnetic tomography

We investigate the reconstruction of a conductive target using crosswell time-domain electromagnetic tomography in the diffusive limit. The work is a natural extension of our ongoing research in the modification of time-domain methods for the rugged marine mid-ocean-ridge environment, an environment characterized by extreme topography and pronounced variations in crustal conductivity on all scales. We have proved both in theory and in practice that 'traveltime', the time taken for an electromagnetic signal to be identified at a receiver following a change of current in the transmitter, is an excellent, robust estimator of average conductivity on a path between transmitter and receiver. A simple estimate of the traveltime for a parallel electric dipole-dipole system is the time at which the derivative of the electric field with respect to logarithmic time at the receiver reaches its maximum. We have derived the fundamental relationship between the traveltime and the conductivity of the medium for a uniform whole-space. We have applied the concept of the traveltime inversion to the related crosswell problem and demonstrated reconstructions of finite targets based on tomographic analyses. Results show that the crosswell time-domain electromagnetic tomography can supply useful information, such as the location and shape of a conductive target.     

A variational principle for the subseismic wave equation

Summary. A variational principle is developed for the subseismic wave equation governing the normal modes of the outer liquid core with frequencies below seismic frequencies (>300 μHz). The calculation of these modes is important both in determining the core contribution to the Earth's dynamical response to tidal and other forces and because their detection at the surface could provide valuable new insight into the density structure of the core, critical to theories of the geomagnetic dynamo. Included as a special case is a variational principle for the Poincaré equation governing inertial oscillations studied in the laboratory by Aldridge and others. This opens up the possibility of 'tracking' laboratory results over to the real Earth.     

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.     

The forward problem of electromagnetic induction: accurate finite-difference approximations for two-dimensional discrete boundaries with arbitrary geometry

We describe finite-difference approximations to the equations of 2-D electromagnetic induction that permit discrete boundaries to have arbitrary geometrical relationships to the nodes. This allows finite-difference modelling with the flexibility normally ascribed to finite-element modelling. Accuracy is demonstrated by comparison with finite-element computations. We also show that related approximations lead to substantially improved accuracy in regions of steep, but not discontinuous, conductivity gradient.     

A modified Lax-Wendroff correction for wave propagation in media described by Zener elements

A modified Lax-Wendroff correction for wave propagation in attenuating and dispersive media described by Zener elements is presented. As opposed to the full correction, this new technique is explicit and offers large computational savings. The technique may be applied to a wide variety of hyperbolic problems. Here, the concept is illustrated for wave propagation in visco-acoustic media.     

Adaptive unstructured grid finite element simulation of two-dimensional magnetotelluric fields for arbitrary surface and seafloor topography

We present an adaptive unstructured triangular grid finite element approach for effectively simulating plane-wave diffusive electromagnetic fields in 2-D conductivity structures.
The most striking advantage of irregular grids is their potential to incorporate arbitrary geometries including surface and seafloor topography. Adaptive mesh refinement strategies using an a posteriori error estimator yield most efficient numerical solutions since meshes are only refined where required.
We demonstrate the robustness of this approach by comparison with analytical solutions and previously published numerical simulations. Maximum errors may systematically be reduced to, for example, 0.8 per cent for the apparent resistivity and 0.2° in the phase.
An additional accuracy study of the thickness of the air layer in E-polarization suggests to keep a minimum thickness depending on lateral conductivity contrasts within the earth.
Furthermore, we point out the new quality and flexibility of our simulation technique by addressing two marine magnetotelluric applications. In the first case, we discuss topographic effects associated with a synthetic sinusoidal sea bottom model and in the second case, we show a close-to-reality scenario using real bathymetry data from the East Pacific Rise at 17°S.