Publication: WSINV3DMT: Vertical magnetic field transfer function inversion and parallel implementation
Issued Date
2009-04-01
Resource Type
ISSN
00319201
Other identifier(s)
2-s2.0-62949241218
Rights
Mahidol University
Rights Holder(s)
SCOPUS
Bibliographic Citation
Physics of the Earth and Planetary Interiors. Vol.173, No.3-4 (2009), 317-329
Suggested Citation
Weerachai Siripunvaraporn, Gary Egbert WSINV3DMT: Vertical magnetic field transfer function inversion and parallel implementation. Physics of the Earth and Planetary Interiors. Vol.173, No.3-4 (2009), 317-329. doi:10.1016/j.pepi.2009.01.013 Retrieved from: https://repository.li.mahidol.ac.th/handle/20.500.14594/27521
Research Projects
Organizational Units
Authors
Journal Issue
Thesis
Title
WSINV3DMT: Vertical magnetic field transfer function inversion and parallel implementation
Author(s)
Other Contributor(s)
Abstract
We describe two extensions to the three-dimensional magnetotelluric inversion program WSINV3DMT (Siripunvaraporn, W., Egbert, G., Lenbury, Y., Uyeshima, M., 2005, Three-dimensional magnetotelluric inversion: data-space method. Phys. Earth Planet. Interiors 150, 3-14), including modifications to allow inversion of the vertical magnetic transfer functions (VTFs), and parallelization of the code. The parallel implementation, which is most appropriate for small clusters, uses MPI to distribute forward solutions for different frequencies, as well as some linear algebraic computations, over multiple processors. In addition to reducing run times, the parallelization reduces memory requirements by distributing storage of the sensitivity matrix. Both new features are tested on synthetic and real datasets, revealing nearly linear speedup for a small number of processors (up to 8). Experiments on synthetic examples show that the horizontal position and lateral conductivity contrasts of anomalies can be recovered by inverting VTFs alone. However, vertical positions and absolute amplitudes are not well constrained unless an accurate host resistivity is imposed a priori. On very simple synthetic models including VTFs in a joint inversion had little impact on the inverse solution computed with impedances alone. However, in experiments with real data, inverse solutions obtained from joint inversion of VTF and impedances, and from impedances alone, differed in important ways, suggesting that for structures with more realistic levels of complexity the VTFs will in general provide useful additional constraints. © 2009 Elsevier B.V.