Magnetotelluric 3-D inversion—a review of two successful workshops on forward and inversion code testing and comparison

Over the last half decade the need for, and importance of, three-dimensional (3-D) modelling of magnetotelluric (MT) data have increased dramatically and various 3-D forward and inversion codes are in use and some have become commonly available. Comparison of forward responses and inversion results is an important step for code testing and validation prior to ‘production’ use. The various codes use different mathematical approximations to the problem (finite differences, finite elements or integral equations), various orientations of the coordinate system, different sign conventions for the time dependence and various inversion strategies. Additionally, the obtained results are dependent on data analysis, selection and correction as well as on the chosen mesh, inversion parameters and regularization adopted, and therefore, a careful and knowledge-based use of the codes is essential. In 2008 and 2011, during two workshops at the Dublin Institute for Advanced Studies over 40 people from academia (scientists and students) and industry from around the world met to discuss 3-D MT inversion. These workshops brought together a mix of code writers as well as code users to assess the current status of 3-D modelling, to compare the results of different codes, and to discuss and think about future improvements and new aims in 3-D modelling. To test the numerical forward solutions, two 3-D models were designed to compare the responses obtained by different codes and/or users. Furthermore, inversion results of these two data sets and two additional data sets obtained from unknown models (secret models) were also compared. In this manuscript the test models and data sets are described (supplementary files are available) and comparisons of the results are shown. Details regarding the used data, forward and inversion parameters as well as computational power are summarized for each case, and the main discussion points of the workshops are reviewed. In general, the responses obtained from the various forward models are comfortingly very similar, and discrepancies are mainly related to the adopted mesh. For the inversions, the results show how the inversion outcome is affected by distortion and the choice of errors, as well as by the completeness of the data set. We hope that these compilations will become useful not only for those that were involved in the workshops, but for the entire MT community and also the broader geoscience community who may be interested in the resolution offered by MT.

[1]  Gary D. Egbert,et al.  WSINV3DMT: Vertical magnetic field transfer function inversion and parallel implementation , 2009 .

[2]  Yongwimon Lenbury,et al.  Numerical accuracy of magnetotelluric modeling: A comparison of finite difference approximations , 2002 .

[3]  D. Oldenburg,et al.  Inversion of electromagnetic data: An overview of new techniques , 1990 .

[4]  V. Červ,et al.  Modelling and analysis of electromagnetic fields in 3D inhomogeneous media , 1990 .

[5]  Alan G. Jones,et al.  Calculations of voltages for magnetotelluric modelling of a region with near-surface inhomogeneities , 1989 .

[6]  K. Spitzer,et al.  3D finite element simulation of magnetotelluric fields using unstructured grids , 2007 .

[7]  Myung Jin Nam,et al.  3D magnetotelluric modelling including surface topography , 2007 .

[8]  William Rodi,et al.  The inverse problem , 2012 .

[9]  Gary D. Egbert,et al.  Computational recipes for electromagnetic inverse problems , 2012 .

[10]  Dmitry B. Avdeev,et al.  3D magnetotelluric inversion using a limited-memory quasi-Newton optimization , 2009 .

[11]  Ralph-Uwe Börner,et al.  Numerical Modelling in Geo-Electromagnetics: Advances and Challenges , 2010 .

[12]  Alan G. Jones,et al.  Static shift of magnetotelluric data and its removal in a sedimentary basin environment , 1988 .

[13]  Alan G. Jones,et al.  Multisite, multifrequency tensor decomposition of magnetotelluric data , 2001 .

[14]  S. Hautot,et al.  Three‐dimensional magnetotelluric inversion of large data sets: Case study of Pasfield Lake (Saskatchewan) for mineral exploration , 2011 .

[15]  J. T. Smith,et al.  Three-dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example , 1994 .

[16]  Gregory A. Newman,et al.  Three‐dimensional induction logging problems, Part I: An integral equation solution and model comparisons , 2002 .

[17]  J. Ledo,et al.  2D inversion of 3D magnetotelluric data: The Kayabe dataset , 1999 .

[18]  S. Hautot,et al.  3‐D magnetotelluric inversion and model validation with gravity data for the investigation of flood basalts and associated volcanic rifted margins , 2007 .

[19]  G. Newman,et al.  Three-dimensional magnetotelluric inversion using non-linear conjugate gradients , 2000 .

[20]  D. Avdeev,et al.  Detecting a salt dome overhang with magnetotellurics: 3D inversion methodology and synthetic model studies , 2012 .

[21]  R. C. Bailey,et al.  Decomposition of magnetotelluric impedance tensors in the presence of local three-dimensional galvanic distortion , 1989 .

[22]  Dmitry B. Avdeev,et al.  Three-Dimensional Electromagnetic Modelling and Inversion from Theory to Application , 2005 .

[23]  William Rodi,et al.  3-D magnetotelluric inversion for resource exploration , 2001 .

[24]  Juanjo Ledo,et al.  2-D Versus 3-D Magnetotelluric Data Interpretation , 2005 .

[25]  Alan G. Jones,et al.  The magnetotelluric method : theory and practice , 2012 .

[26]  Colin Farquharson,et al.  Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction , 2011 .

[27]  Max A. Meju,et al.  Three-dimensional joint inversion for magnetotelluric resistivity and static shift distributions in complex media , 2006 .

[28]  Michael S. Zhdanov,et al.  Methods for modelling electromagnetic fields Results from COMMEMI—the international project on the comparison of modelling methods for electromagnetic induction , 1997 .

[29]  Douglas W. Oldenburg,et al.  An Algorithm For the Three-dimensional Inversion of Magnetotelluric Data , 2002 .

[30]  Alan G. Jones,et al.  Artefacts of isotropic inversion applied to magnetotelluric data from an anisotropic Earth , 2011 .

[31]  Adam Schultz,et al.  Introduction to MT-DIW2 Special Issue , 1997 .

[32]  Alan G. Jones.  Three‐dimensional galvanic distortion of three‐dimensional regional conductivity structures: Comment on “Three‐dimensional joint inversion for magnetotelluric resistivity and static shift distributions in complex media” by Yutaka Sasaki and Max A. Meju , 2011 .

[33]  Alan G. Jones,et al.  The COPROD2 Dataset: Tectonic Setting, Recorded MT Data, and Comparison of Models , 1993 .

[34]  Yongwimon Lenbury,et al.  Three-dimensional magnetotelluric inversion : data-space method , 2005 .

[35]  Alan G. Jones,et al.  The BC87 Dataset: Tectonic Setting, Previous EM Results, and Recorded MT Data , 1993 .

[36]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[37]  R. W. Groom,et al.  Analytic investigations of the effects of near‐surface three‐dimensional galvanic scatterers on MT tensor decompositions , 1991 .

[38]  Y. Sasaki.  Full 3-D inversion of electromagnetic data on PC , 2001 .

[39]  Weerachai Siripunvaraporn,et al.  Three-Dimensional Magnetotelluric Inversion: An Introductory Guide for Developers and Users , 2011, Surveys in Geophysics.