WSJointInv2D-MT-DCR: An efficient joint two-dimensional magnetotelluric and direct current resistivity inversion

An efficient joint two-dimensional direct current resistivity (DCR) and magnetotelluric (MT) inversion, referred to as WSJointInv2D-MT-DCR, was developed with FORTRAN 95 based on the data space Occam's inversion algorithm. Our joint inversion software can be used to invert just the MT data or the DCR data, or invert both data sets simultaneously to get the electrical resistivity structures. Since both MT and DCR surveys yield the same resistivity structures, the two data types enhance each other leading to a better interpretation. Two synthetic and a real field survey are used here to demonstrate that the joint DCR and MT surveys can help constrain each other to reduce the ambiguities occurring when inverting the DCR or MT alone. The DCR data increases the lateral resolution of the near surface structures while the MT data reveals the deeper structures. When the MT apparent resistivity suffers from the static shift, the DCR apparent resistivity can serve as a replacement for the estimation of the static shift factor using the joint inversion. In addition, we also used these examples to show the efficiency of our joint inversion code. With the availability of our new joint inversion software, we expect the number of joint DCR and MT surveys to increase in the future. Joint inversion software for 2D magnetotelluric and direct current resistivity.Both data help enhance each other reducing the ambiguities of the interpretation.The software is free for academic uses.

[1]  Chatchai Vachiratienchai,et al.  An efficient inversion for two-dimensional direct current resistivity surveys based on the hybrid finite difference-finite element method , 2013 .

[2]  C. Farquharson,et al.  Three-dimensional inversion of magnetotelluric data for mineral exploration: An example from the McArthur River uranium deposit, Saskatchewan, Canada , 2009 .

[3]  Nong Wu,et al.  Rapid Two-Dimensional Inversion of COPROD2 Data , 1993 .

[4]  Bülent Tezkan,et al.  A radiomagnetotelluric survey on an oil‐contaminated area near the Brazi Refinery, Romania , 2005 .

[5]  Klaus Spitzer,et al.  Magnetotelluric static shift and direct current sensitivity , 2001 .

[6]  B. B. Bhattacharya,et al.  Audiomagnetotelluric studies to trace the hydrological system of thermal fluid flow of Bakreswar Hot Spring, Eastern India: A case history , 2010 .

[7]  R. Barker,et al.  Rapid least-squared inversion of apparent resisitivity pseudosections by a quasi-Newton method , 1996 .

[8]  M. Loddo,et al.  Shallow to intermediate resistivity features of the Colfiorito Fault System inferred by DC and MT survey , 2008 .

[9]  M. Uyeshima,et al.  Interpretation of two-dimensional magnetotelluric profile data with three-dimensional inversion: synthetic examples , 2005 .

[10]  Fernando A. Monteiro Santos,et al.  2D joint inversion of dc and scalar audio-magnetotelluric data in the evaluation of low enthalpy geothermal fields , 2007 .

[11]  K. Spitzer,et al.  A method to determine the magnetotelluric static shift from DC resistivity measurements in practice , 2010 .

[12]  Yasuo Ogawa,et al.  On Two-Dimensional Modeling Of Magnetotelluric Field Data , 2002 .

[13]  A. Hördt,et al.  A joint application of radiomagnetotellurics and transient electromagnetics to the investigation of a waste deposit in Cologne (Germany) , 1996 .

[14]  Thomas Kalscheuer,et al.  Non-linear model error and resolution properties from two-dimensional single and joint inversions of direct current resistivity and radiomagnetotelluric data , 2010 .

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

[16]  B. Tezkan,et al.  A new broadband radiomagnetotelluric instrument: applications to near surface investigations , 2008 .

[17]  Bülent Tezkan,et al.  Two‐dimensional joint inversion of radiomagnetotelluric and direct current resistivity data , 2008 .

[18]  P. Turberg,et al.  Radiomagnetotelluric mapping, groundwater numerical modelling and 18-Oxygen isotopic data as combined tools to determine the hydrogeological system of a landslide prone area , 2006 .

[19]  Yutaka Sasaki,et al.  Two‐dimensional joint inversion of magnetotelluric and dipole‐dipole resistivity data , 1989 .

[20]  Robert L. Parker,et al.  Optimal one-dimensional inversion and bounding of magnetotelluric apparent resistivity and phase measurements , 1996 .

[21]  Gary D. Egbert,et al.  Data space conjugate gradient inversion for 2-D magnetotelluric data , 2007 .

[22]  Lena Persson,et al.  Delineation of a quick clay zone at Smørgrav, Norway, with electromagnetic methods under geotechnical constraints , 2013 .

[23]  Y. Ogawa,et al.  Evidence for middle Triassic to Miocene dual subduction zones beneath the Shan-Thai terrane, western Thailand from magnetotelluric data , 2013 .

[24]  F. Santos,et al.  1D joint inversion of AMT and resistivity data acquired over a graben , 1997 .

[25]  B. Tezkan,et al.  Radiomagnetotelluric and Direct Current Resistivity measurements for the characterization of conducting soils , 2007 .

[27]  S. Wood,et al.  Exploring the shallow geothermal fluid reservoir of Fang geothermal system, Thailand via a 3-D magnetotelluric survey , 2016 .

[28]  W. Rodi.  A Technique for Improving the Accuracy of Finite Element Solutions for Magnetotelluric Data , 1976 .

[29]  W. Siripunvaraporn,et al.  An efficient data space conjugate gradient Occam's method for three‐dimensional magnetotelluric inversion , 2011 .

[30]  A. Malehmir,et al.  Integrated 2D modeling and interpretation of geophysical and geotechnical data to delineate quick clays at a landslide site in southwest Sweden , 2014 .

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

[32]  S. L. Fontes,et al.  Regional aquifer mapping using combined VES-TEM-AMT/EMAP methods in the semiarid eastern margin of Parnaiba Basin, Brazil , 1999 .

[33]  Thomas Kalscheuer,et al.  2D joint inversion of RMT and ERT data versus individual 3D inversion of full tensor RMT data. An example from Trecate site in Italy , 2012 .

[34]  Niklas Linde,et al.  Characterization of a fractured granite using radio magnetotelluric (RMT) data , 2004 .

[35]  Chatchai Vachiratienchai,et al.  Two-dimensional direct current (DC) resistivity inversion: Data space Occam's approach , 2008 .

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

[37]  Y. Ogawa,et al.  A 3-D conductivity model of the Australian continent using observatory and magnetometer array data , 2014 .

[38]  S. Constable,et al.  Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data , 1990 .

[39]  On the 3-D inversion of vertical electrical soundings: Application to the South Ismailia area—Cairo desert road, Cairo, Egypt , 2008 .

[40]  Í. Vitorello,et al.  Magnetotelluric soundings across the Taubaté Basin, Southeast Brazil , 2002 .

[41]  G. Buselli,et al.  Groundwater contamination monitoring with multichannel electrical and electromagnetic methods , 2001 .

[42]  Bülent Tezkan,et al.  Groundwater contamination in the Roorkee area, India: 2D joint inversion of radiomagnetotelluric and direct current resistivity data , 2012 .

[43]  William Rodi,et al.  Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion , 2001 .

[44]  W. Siripunvaraporn,et al.  A pilot magnetotelluric survey for geothermal exploration in Mae Chan region, northern Thailand , 2015 .

[45]  J. Craven,et al.  Exploration for unconformity-type uranium deposits with audiomagnetotelluric data : A case study from the McArthur River mine, Saskatchewan, Canada , 2006 .

[46]  Gary D. Egbert,et al.  An efficient data-subspace inversion method for 2-D magnetotelluric data , 2000 .

[47]  S. Wilson,et al.  The applicability of earth resistivity methods for saline interface definition , 2006 .