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.     

Geoelectromagnetic induction in a heterogeneous sphere:a new three-dimensional forward solver using a conservative staggered-grid finite difference method

A conservative staggered-grid finite difference method is presented for computing the electromagnetic induction response of an arbitrary heterogeneous conducting sphere by external current excitation. This method is appropriate as the forward solution for the problem of determining the electrical conductivity of the Earth's deep interior. This solution in spherical geometry is derived from that originally presented by Mackie et al. (1994 ) for Cartesian geometry. The difference equations that we solve are second order in the magnetic field H , and are derived from the integral form of Maxwell's equations on a staggered grid in spherical coordinates. The resulting matrix system of equations is sparse, symmetric, real everywhere except along the diagonal and ill-conditioned. The system is solved using the minimum residual conjugate gradient method with preconditioning by incomplete Cholesky decomposition of the diagonal sub-blocks of the coefficient matrix. In order to ensure there is zero H divergence in the solution, corrections are made to the H field every few iterations. In order to validate the code, we compare our results against an integral equation solution for an azimuthally symmetric, buried thin spherical shell model ( Kuvshinov & Pankratov 1994 ), and against a quasi-analytic solution for an azimuthally asymmetric configuration of eccentrically nested spheres ( Martinec 1998 ).     

Geomagnetic induction in a heterogeneous sphere: fully three-dimensional test computations and the response of a realistic distribution of oceans and continents

Long-period geomagnetic data can resolve large-scale 3-D mantle electrical conductivity heterogeneities which are indicators of physiochemical variations found in the Earth's dynamic mantle. A prerequisite for mapping such heterogeneity is the ability to model accurately electromagnetic induction in a heterogeneous sphere. A previously developed finite element method solution to the geomagnetic induction problem is validated against an analytic solution for a fully 3-D geometry: an off-axis spherical inclusion embedded in a uniform sphere. Geomagnetic induction is then modelled in a uniform spherical mantle overlain by a realistic distribution of oceanic and continental conductances. Our results indicate that the contrast in electrical conductivity between oceans and continents is not primarily responsible for the observed geographic variability of long-period geomagnetic data. In the absence of persistent high-wavenumber magnetospheric disturbances, this argues strongly for the existence of large-scale, high-contrast electrical conductivity heterogeneities in the mid-mantle. Lastly, for several periods the geomagnetic anomaly associated with a mid-mantle spherical inclusion is calculated. A high-contrast inclusion can be readily detected beneath the outer shell of oceans and continents. A comparison between observed and computed c responses suggests that the mid-mantle contains more than one order of magnitude of lateral variability in electrical conductivity, while the upper mantle contains at least two orders of magnitude of lateral variability in electrical conductivity.