Three-dimensional inversion of magnetotelluric phase tensor data

SUMMARY Recent increase in the application of 3-D inversion of magnetotelluric (MT) data is facilitated by the availability of several 3-D inversion codes, which led to improved interpretation of MT data. However, still the galvanic effects continue to pose problems in interpretation of the MT data. We have addressed this problem using the phase tensor (PT) concept and developed a scheme based on a modification of the sensitivity matrix in the 3-D inversion code of Siripunvaraporn et al., which enables us to directly invert the phase tensor elements. We have used this modified code for PT inversion of MT data and evaluated its efficacy in reducing the galvanic effects through a few examples of inversion of synthetic data and its application on real data. The synthetic model study suggests that the prior model (mo) setting is important in retrieving the true model. Comparison of results obtained from conventional impedance inversion and the proposed PT inversion method suggests that, even in the presence of the galvanicdistortion(duetonear-surfacecheckerboardanomaliesinourcase),thenewinversion algorithmreliablyretrievestheregionalconductivitystructurewhenthepriormodelorregional resistivity value level can be estimated with sufficient accuracy. Data errors were propagated to the PT elements using delta method while inverting the real data from USArray. The PT results compare very well with those from full tensor inversion published earlier, signifying the efficacy of this new inversion scheme.

[1]  G. Egbert,et al.  Regional conductivity structure of Cascadia: Preliminary results from 3D inversion of USArray transportable array magnetotelluric data , 2008 .

[2]  H. Bibby,et al.  Three‐dimensional modelling of magnetotelluric data from the Rotokawa geothermal field, Taupo Volcanic Zone, New Zealand , 2008 .

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

[4]  Gary D. Egbert,et al.  Robust multiple‐station magnetotelluric data processing , 1997 .

[5]  K. Bahr,et al.  Interpretation of the magnetotelluric impedance tensor; regional induction and local telluric distortion , 1988 .

[6]  B. Efron The jackknife, the bootstrap, and other resampling plans , 1987 .

[7]  A. Bennett,et al.  TOPEX/POSEIDON tides estimated using a global inverse model , 1994 .

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

[9]  Gylfi Páll Hersir,et al.  Joint 1D inversion of TEM and MT data and 3D inversion of MT data in the Hengill area, SW Iceland , 2010 .

[10]  M. Meju Geoelectromagnetic Exploration For Natural Resources: Models, Case Studies And Challenges , 2002 .

[11]  H. Bibby,et al.  Three‐dimensional electrical resistivity image of magma beneath an active continental rift, Taupo Volcanic Zone, New Zealand , 2010 .

[12]  M. Zhdanov,et al.  Large-scale three-dimensional inversion of EarthScope MT data using the integral equation method , 2010 .

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

[14]  Vladimir I. Dmitriev,et al.  Magnetotellurics in the Context of the Theory of Ill-Posed Problems , 2002 .

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

[16]  Alan D. Chave,et al.  On electric and magnetic galvanic distortion tensor decompositions , 1994 .

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

[18]  John R. Booker,et al.  The Magnetotelluric Phase Tensor: A Critical Review , 2013, Surveys in Geophysics.

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

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

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

[22]  M. Uyeshima EM Monitoring of Crustal Processes Including the Use of the Network-MT Observations , 2007 .