Inversion of magnetotelluric data using integral equation approach with variable sensitivity domain: application to earthscope MT data

We have developed a multi-level parallel magnetotelluric (MT) integral equation based inversion program which uses variable sensitivity domain. The limited sensitivity of the data, which decreases with increasing frequency, is exploited by a receiver sensitivity domain, which also varies with frequency. We assess the effect of inverting principal, full impedance tensor, and full tensor with magnetovariational data (tipper). We first apply this method to several models and then invert the EarthScope MT data. We recover well the prominent features in the area including resistive structure associated with the Juan de Fuca slab subducting beneath the northwestern United States, the conductive zone of partially melted material above the subducting slab at the Cascade volcanic arc, conductive features in the Great Basin resulting from prior subduction and in the area of Yellowstone associated with the hot spot, and resistive areas to the east corresponding to the older and more stable cratons. Our results agree well with other publications which furthers the confidence in the inversion of the EarthScope MT project.

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

[2]  G. Keller,et al.  Advanced theory of deep geomagnetic sounding , 1984 .

[3]  M. Zhdanov,et al.  A multigrid integral equation method for large-scale models with inhomogeneous backgrounds , 2008 .

[4]  Brandon Schmandt,et al.  Complex subduction and small-scale convection revealed by body-wave tomography of the western United States upper mantle , 2010 .

[5]  Michael S. Zhdanov,et al.  Contraction integral equation method in three‐dimensional electromagnetic modeling , 2002 .

[6]  Zonghou Xiong.  Electromagnetic modeling of 3-D structures by the method of system iteration using integral equations , 1992 .

[7]  D. Sandwell,et al.  Global seafloor topography from dense satellite altimetry and sparse ship soundings , 1997 .

[8]  R. Carlson,et al.  Slab fragmentation, edge flow and the origin of the Yellowstone hotspot track , 2011 .

[9]  Mei Xue,et al.  Slab‐plume interaction beneath the Pacific Northwest , 2010 .

[10]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

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

[12]  Philip E. Wannamaker,et al.  Advances in three-dimensional magnetotelluric modeling using integral equations , 1991 .

[13]  K. Dueker,et al.  Temperature of the plume layer beneath the Yellowstone hotspot , 2008 .

[14]  Michael S. Zhdanov,et al.  Geophysical Electromagnetic Theory and Methods , 2009 .

[15]  Michael I. Jordan,et al.  Geodynamics of the Yellowstone hotspot and mantle plume: Seismic and GPS imaging, kinematics, and mantle flow , 2009 .

[16]  Matthew J. Fouch,et al.  The Yellowstone Hotspot: Plume or Not? , 2012 .

[17]  R. Carlson,et al.  Three‐dimensional seismic velocity structure of the northwestern United States , 2008 .

[18]  Robert B. Smith,et al.  Three‐dimensional inversion of large‐scale EarthScope magnetotelluric data based on the integral equation method: Geoelectrical imaging of the Yellowstone conductive mantle plume , 2011 .

[19]  R. D. Hibbs,et al.  Electromagnetic induction in three dimensional structures for various source fields. , 1978 .

[20]  P. Bedrosian,et al.  Structure and tectonics of the northwestern United States from EarthScope USArray magnetotelluric data , 2014 .

[21]  E. Humphreys,et al.  Seismically imaged relict slab from the 55 Ma Siletzia accretion to the northwest United States , 2011 .

[22]  Adam Schultz,et al.  Deep electrical resistivity structure of the northwestern U.S. derived from 3-D inversion of USArray magnetotelluric data , 2014 .

[23]  P. M. Berg,et al.  Iterative forward and inverse algorithms based on domain integral equations for three-dimensional electric and magnetic objects , 2004 .

[24]  K. Dueker,et al.  Crust and upper mantle velocity structure of the Yellowstone hot spot and surroundings , 2008 .

[25]  Robert B. Smith,et al.  Models of lithosphere and asthenosphere anisotropic structure of the Yellowstone hot spot from shear wave splitting , 2005 .

[26]  A. Pommier,et al.  Laboratory measurements of electrical conductivities of hydrous and dry Mount Vesuvius melts under pressure , 2008 .

[27]  Mario Martinez,et al.  Magnetotelluric observations across the Juan de Fuca Subduction System in the EMSLAB Project , 1989 .

[28]  G. W. Hohmann.  Three-Dimensional Induced Polarization and Electromagnetic Modeling , 1975 .

[29]  M. Zhdanov,et al.  Inverse Theory and Applications in Geophysics , 2015 .

[30]  R. Allen,et al.  Lithosphere-asthenosphere interaction beneath the western United States from the joint inversion of body-wave traveltimes and surface-wave phase velocities , 2011 .

[31]  Anna Avdeeva,et al.  uasi-Newton optimization , 2009 .

[32]  M. Zhdanov,et al.  Geoelectrical Structure of the Lithosphere and Asthenosphere beneath the Northwestern United States , 2012 .

[33]  A. Raiche.  An Integral Equation Approach to Three-Dimensional Modelling , 1974 .

[34]  Michael S. Zhdanov,et al.  3D inversion of airborne electromagnetic data using a moving footprint , 2010 .

[35]  F. X. Bostick,et al.  Geothermal significance of magnetotelluric sounding in the Eastern Snake River Plain-Yellowstone Region , 1977 .

[36]  Michael S. Zhdanov,et al.  Integral equation method for 3-D modeling of electromagnetic fields in complex structures with inhomogeneous background conductivity in marine CSEM applications , 2005 .