Modelling of electromagnetic fields in thin heterogeneous layers with application to field generation by volcanoes—theory and example |
| |
| Authors: | B Sh Singer & E B Fainberg |
| |
| Institution: | Cooperative Research Centre for Australian Mineral Exploration Technology, Building E5A, Macquarie University, North Ryde, NSW 2109, Australia. E-mail:;, Institute of Geoelectromagnetic Investigations, Russian Academy of Sciences, 142092 Moscow Region, Russia |
| |
| Abstract: | When interpreting electromagnetic fields observed at the Earth's surface in a realistic geophysical environment it is often necessary to pay special attention to the effects caused by inhomogeneities of the subsurface sedimentary and/or water layer and by inhomogeneities of the Earth's crust. The inhomogeneities of the Earth's crust are expected to be especially important when the electromagnetic field is generated by a source located in a magma chamber of a volcano. The simulation of such effects can be carried out using generalized thin-sheet models, which were independently introduced by Dmitriev (1969 ) and Ranganayaki & Madden (1980 ). In the first part of the paper, a system of integral equations is derived for the horizontal current that flows in the subsurface inhomogeneous conductive layer and for the vertical current crossing the inhomogeneous resistive layer representing the Earth's mantle. The terms relating to the finite thickness of the laterally inhomogeneous part of the model are retained in the equations. This only marginally complicates the equations, whilst allowing for a significant expansion of the approximation limits.
The system of integral equations is solved using the iterative dissipative method developed by the authors in the period from 1978 to 1988. The method can be applied to the simulation of the electromagnetic field in an arbitrary inhomogeneous medium that dissipates the electromagnetic energy. When considered on a finite numerical grid, the integral equations are reduced to a system of linear equations that possess the same contraction properties as the original equations. As a result, the rate at which the iterative-perturbation sequence converges to the solution remains independent of the numerical grid used for the calculations. In contrast to previous publications on the method, aspects of the algorithm implementation that guarantee its effectiveness and robustness are discussed here. |
| |
| Keywords: | crust generalized thin sheet iterative dissipative method subsurface inhomogeneities volcano |
|
|
