1887
  1. Home
  2. A-Z Publications
  3. Geophysical Prospecting
  4. Volume 66, Issue 9
  5. Article
Volume 66, Issue 9
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Electrical anisotropy in earth media increases the complexity of magnetotelluric responses. Magnetotelluric models based on anisotropic media must be developed to fully understand observed data. This paper presents a three‐dimensional algorithm for calculating magnetotelluric responses of arbitrary anisotropic media in the frequency domain. Using a staggered‐grid finite difference method, the model space is discretized into rectangular blocks with electric fields on the edges of each block and the magnetic fields normal to the faces of each block. The electric field Helmholtz vector equation that considers a full 3 × 3 conductivity tensor is calculated numerically under two orthogonal polarizations. In calculating the boundary values on the four sides of the three‐dimensional anisotropic model, we adopt different procedures for calculating the two‐dimensional responses on the sides in the and directions. The responses for a layered anisotropic model and a three‐dimensional isotropic model calculated with this algorithm are compared with the corresponding analytical and numerical solutions, respectively. The comparisons show that the algorithm's approximations are highly precise for a wide frequency band. A typical two‐dimensional anisotropic model and a general three‐dimensional anisotropic model were also constructed, and their responses were calculated. These anisotropic models have ordinary structures but can produce phase rolling out of quadrant magnetotelluric responses, which indicates that considering electrical anisotropy may improve our interpretation of observed data. Using this algorithm, we can model the observed data from the northern Qaidam Basin in northern Tibet, where ultrahigh‐pressure metamorphic rocks are exposed along an old suture, and seismic anisotropy was indicated in neighbouring areas. The phase tensors of the magnetotelluric sites at this location show large skew angles, and the corresponding phase splits are distinct in the off‐diagonal impedance elements. Although the isotropic three‐dimensional electrical structure can model the profile data well, the structure shows a sequence of conductive and resistive bodies in the mid‐lower crust of the north Qaidam Basin, which is very spatially inhomogeneous, and a simple intrinsic anisotropic body can also produce similar surficial responses. Using the three‐dimensional anisotropic algorithm, we found an equivalent anisotropic replacement for this area. The results of the three‐dimensional anisotropy modelling of the magnetotelluric data from the northern Tibetan Plateau show the valuable applicability of the three‐dimensional anisotropic algorithm in testing the qualitative presumption of electrical anisotropy.

© 2018 European Association of Geoscientists & Engineers © 2018 European Association of Geoscientists & Engineers
Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12690
2018-10-23
2024-04-23

  • From This Site

    /content/journals/10.1111/1365-2478.12690
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. AbramoviciF.1974. The forward magnetotelluric problem for an inhomogeneous and anisotropic structure. Geophysics39, 56–68.
     [Google Scholar]
  2. AbramoviciF. and ShohamY.1977. Inversion of anisotropic magnetotelluric data. Geophysical Journal International of the Royal Astronomical Society50, 55–74.
     [Google Scholar]
  3. AlumbaughD.L., NewmannG.A., PrevostL. and ShadidJ.1996. Three‐dimensional wide‐band electromagnetic modeling on massively parallel computers. Radio Science31, 1–23.
     [Google Scholar]
  4. AvdeevD.B.2005. Three‐dimensional electromagnetic modeling and inversion from theory to application. Surveys in Geophysics. 26, 767–799.
     [Google Scholar]
  5. BabaK., ChaveA.D., EvansR.L., HirthG. and MackieR.L.2006. Mantle dynamics beneath the East Pacific Rise at 17°S: Insights from the mantle electromagnetic and tomography (MELT) experiment. Journal of Geophysical Research Solid Earth111, B02101.
     [Google Scholar]
  6. BahrK., BantinM., JantosC., SchneiderE. and StorzW.2000. Electrical anisotropy from electromagnetic array data: Implications for the conduction mechanism and for distortion at long periods. Physics of the Earth and Planetary Interiors119, 237–257.
     [Google Scholar]
  7. BahrK. and SimpsonF.2002. Electrical anisotropy below slow and fast‐moving plates: Paleoflow in the upper mantle? Science295, 1270–1272.
     [Google Scholar]
  8. BrasseH., KapinosG., LiY., MütschardL., SoyerW. and EydamD.2009. Structural electrical anisotropy in the crust at the South‐Central Chilean continental margin as inferred from geomagnetic transfer functions. Physics of the Earth and Planet Interiors173, 7–16.
     [Google Scholar]
  9. CaldwellT.G., BibbyH.M. and BrownC.2004. The magnetotelluric phase tensor. Geophysical Journal International158, 457–469.
     [Google Scholar]
  10. DavydychevaS., DruskinV. and HabashyT.2003. An efficient finite‐difference scheme for electromagnetic logging in 3D anisotropic inhomogeneous media. Geophysics68, 1525–1536.
     [Google Scholar]
  11. DavydychevaS. and FrenkelM. A.2010. Review of 3D EM modeling and interpretation methods for triaxial induction and propagation resistivity well logging tools. PIERS Proceedings, Cambridge, USA, July 2010. Curan.
  12. de la Kethulle de RyhoveS. and MittetR.2014. 3D marine magnetotelluric modeling and inversion with the finite‐difference time‐domain method. Geophysics79, 269–286.
     [Google Scholar]
  13. DekkerD.L. and HastieL.M.1980. Magneto‐telluric impedances of an anisotropic layered Earth model. Geophysical Journal International of the Royal Astronomical Society61, 11–20.
     [Google Scholar]
  14. EatonD.W., JonesA.G. and FergusonI.J.2004. Lithospheric anisotropy structure from collocated teleseismic and magnetotelluric observations: Great Slave Lake shear zone, Northern Canada. Geophysical Research Letters31, L19614.
     [Google Scholar]
  15. EgbertG.D.1990. Comments on ‘concerning dispersion relations for the magnetotelluric impedence tensor’ by Yee E. and Paulson K.V. Geophysical Journal International of the Royal Astronomical Society102, 1–8.
     [Google Scholar]
  16. EgbertG.D. and KelbertA.2012. Computational recipes for electromagnetic inverse problems. Geophysical Journal International189, 251–267.
     [Google Scholar]
  17. EiselM. and HaakV.1999. Macro‐anisotropy of the electrical conductivity of the crust: A magnetotelluric study of the German continental deep drilling site (KTB). Geophysical Journal International136, 109–122.
     [Google Scholar]
  18. EvansR.L., HirthG., BabaK., ForsythD., ChaveA., and MackieR.2005. Geophysical evidence from the melt area for compositional controls on oceanic plates. Nature437, 249–252.
     [Google Scholar]
  19. FomenkoE.Y. and MogiT.2002. A new computation method for a staggered grid of 3D EM fields conservative modeling. Earth Planets and Space54, 499−509.
     [Google Scholar]
  20. GatzemeierA. and MoorkampM.2005. 3D modelling of electrical anisotropy from electromagnetic array data: Hypothesis testing for different upper mantle conduction mechanisms. Physics of the Earth and Planetary Interiors149, 225–242.
     [Google Scholar]
  21. HeiseW., CaldwellT.G., BibbyH.M. and BrownC.2006. Anisotropy and phase splits in magnetotellurics. Physics of the Earth and Planetary Interiors158, 107–121.
     [Google Scholar]
  22. HeiseW. and PousJ.2001. Effects of anisotropy on the two‐dimensional inversion procedure. Geophysical Journal International147, 610–621.
     [Google Scholar]
  23. HeiseW. and PousJ.2003. Anomalous phases exceeding 90° in magnetotellurics: Anisotropic model studies and a field example. Geophysical Journal International155, 308–318.
     [Google Scholar]
  24. HerquelG., WittlingerG. and GuilbertJ.1995. Anisotropy and crustal thickness of Northern‐Tibet. New constraints for tectonic modeling. Geophysical Research Letters22, 1925–1928.
     [Google Scholar]
  25. IchiharaH. and MogiT.2009. A realistic 3‐D resistivity model explaining anomalous large magnetotelluric phases: the L‐shaped conductor model. Geophysical Journal International179, 14–17.
     [Google Scholar]
  26. JaysavalP., ShantsevD.V., de la Kethulle de RyhoveS. and BrattelandT.2016. Fully anisotropic 3‐D EM modelling on a Lebedev grid with a multigrid pre‐conditioner. Geophysical Journal International207, 1554−1572.
     [Google Scholar]
  27. JonesA.G., EvansR.L., MullerM.R., HamiltonM.P., MiensopustM.P., GarciaX., et al. 2009. Area selection for diamonds using magnetotellurics: Examples from southern Africa. Lithos112, 83–92.
     [Google Scholar]
  28. JungeA.2016. 3D anisotropic modelling and EM transfer functions. The 23rd Electromagnetic Induction Workshop, Chiang Mai, Thailand.
  29. KelbertA., MeqbelN., EgbertG.D. and TandonK.2014. ModEM: A modular system for inversion of electromagnetic geophysical data. Computers & Geoscience66, 40–53.
     [Google Scholar]
  30. KellettR.L., MareschalM. and KurtzR.D.1992. A model of lower crustal electrical anisotropy for the Pontiac subprovince of the Canadian shield. Geophysical Journal International111, 141–150.
     [Google Scholar]
  31. KumarG.P. and ManglikA.2012. Electrical anisotropy in the main central thrust zone of the Sikkim Himalaya: Inference from anomalous MT phase. Journal of Asian Earth Sciences57, 120–127.
     [Google Scholar]
  32. KurtzR.D., CravenJ.A., NiblettE.R. and StevensR.A.1993. The conductivity of the crust and mantle beneath the Kapuskasing uplift: Electrical anisotropy in the upper mantle. Geophysical Journal International113, 483–498.
     [Google Scholar]
  33. KuzminA., LuisierM. and SchenkO.2013. Fast methods for computing selected elements of the Greens function in massively parallel nanoelectronic device simulations. In: Euro‐Par 2013 Parallel Processing (eds F.Wolf , B.Mohr and D.Mey ), pp. 533–544. Berlin: Springer.
     [Google Scholar]
  34. LebedevV.I.1964. Difference analogies of orthogonal decompositions of basic differential operators and some boundary value problems. Soviet Computational Mathematics and Mathematical Physics4, 449–465.
     [Google Scholar]
  35. LeibeckerJ., GatzemeierA., HönigM., KurasO. and SoyerW.2002. Evidence of electrical anisotropic structures in the lower crust and the upper mantle beneath the Rhenish shield. Earth and Planetary Science Letters202, 289–302.
     [Google Scholar]
  36. Le PapeF., JonesA.G., VozarJ. and WeiW.2012. Penetration of crustal melt beyond the Kunlun Fault into northern Tibet. Nature Geoscience5, 330–335.
     [Google Scholar]
  37. LezaetaP. and HaakV.2003. Beyond magnetotelluric decomposition: Induction, current channeling, and magnetotelluric phases over 90°. Journal of Geophysical Research Solid Earth108, 2305.
     [Google Scholar]
  38. LiY.2002. A finite‐element algorithm for electromagnetic induction in two‐dimensional anisotropic conductivity structures. Geophysical Journal International148, 389–401.
     [Google Scholar]
  39. LiddellM., UnsworthM. and PekJ.2016. Magnetotelluric imaging of anisotropic crust near Fort McMurray, Alberta: Implications for engineered geothermal system development. Geophysical Journal International205, 1365–1381.
     [Google Scholar]
  40. LiuY., XuZ. and LiY.2018. Adaptive finite element modelling of three‐dimensional magnetotelluric fields in general anisotropic media. Journal of Applied Geophysics151, 113–124.
     [Google Scholar]
  41. LivelybrooksD., MareschalM., BlaisE. and SmithJ.T.1996. Magnetotelluric delineation of the trillabelle massive sulfide body in Sudbury, Ontario. Geophysics61, 971–986.
     [Google Scholar]
  42. MackieR.L., MaddenT.R. and WannamakerP.E.1993. Three‐dimensional magnetotelluric modeling using difference equations: Theory and comparisons to integral equation solutions. Geophysics58, 215–226.
     [Google Scholar]
  43. MaddenT.R. and MackieR.1989. Three‐dimensional magnetotelluric modeling and inversion. Proceedings of the IEEE77, 318–333.
     [Google Scholar]
  44. MannJ.E.1965. The importance of anisotropic conductivity in magnetotelluric interpretation. Journal of Geophysical Research70, 2940–2942.
     [Google Scholar]
  45. MartinelliP. and OsellaA.1997. MT forward modeling of 3‐D anisotropic electrical conductivity structures using the Rayleigh‐Fourier method. Journal of Geomagnetism and Geoelectricity49, 1499–1518.
     [Google Scholar]
  46. MartíA.2014. The role of electrical anisotropy in magnetotelluric responses: From modelling and dimensionality analysis to inversion and interpretation. Surveys in Geophysics35, 179–218.
     [Google Scholar]
  47. MiensopustM.P., QueraltP. and JonesA.G.2013. Magnetotelluric 3‐D inversion–a review of two successful workshops on forward and inversion code testing and comparison. Geophysical Journal International193, 1216–1238.
     [Google Scholar]
  48. O'BrienD.P. and MorrisonH.F.1967. Electromagnetic fields in an n‐layer anisotropic half‐space. Geophysics32, 668–677.
     [Google Scholar]
  49. PadilhaA.L., VitorelloI., PáduaM.B. and BolognaM.S.2006. Lithospheric and sublithospheric anisotropy beneath central‐southeastern Brazil constrained by long period magnetotelluric data. Physics of the Earth and Planetary Interiors158, 190–209.
     [Google Scholar]
  50. PekJ. and SantosFAM.2002. Magnetotelluric impedances and parametric sensitivities for 1‐D anisotropic layered media. Computers & Geosciences28, 939–950.
     [Google Scholar]
  51. PekJ. and VernerT.1997. Finite‐difference modelling of magnetotelluric fields in two‐dimensional anisotropic media. Geophysical Journal International128, 505–521.
     [Google Scholar]
  52. ReddyI.K. and RankinD.1971. Magnetotelluric effect of dipping anisotropies. Geophysical Prospecting19, 84–97.
     [Google Scholar]
  53. ReddyI.K. and RankinD.1975. Magnetotelluric response of laterally inhomogeneous and anisotropic media. Geophysics40, 1035–1045.
     [Google Scholar]
  54. SchenkO. and GärtnerK.2004. Solving unsymmetric sparse systems of linear equations with PARDISO. Journal of Future Generation Computer Systems20, 475–487.
     [Google Scholar]
  55. SimpsonF.2001. Resistance to mantle flow inferred from the electromagnetic strike of the Australian upper mantle. Nature412, 632–635.
     [Google Scholar]
  56. SiripunvarapornW., EgbertG. and LenburyY.2002. Numerical accuracy of magnetotelluric modeling: A comparison of finite difference approximations. Earth Planets and Space54, 721–725.
     [Google Scholar]
  57. SmithJ.T.1996. Conservative modelling of 3‐D electromagnetic fields. Geophysics61, 1308–1324.
     [Google Scholar]
  58. TapponnierP., XuZ., RogerF., MeyerB., ArnaudN., WittlingerG., et al. 2001. Oblique stepwise rise and growth of the Tibet plateau. Science294, 1671−1677.
     [Google Scholar]
  59. WangT. and FangS.2001. 3‐D electromagnetic anisotropy modeling using finite differences. Geophysics66, 1386–1398.
     [Google Scholar]
  60. WannamakerP.E.2005. Anisotropy versus heterogeneity in continental solid earth electromagnetic studies: Fundamental response characteristics and implications for physicochemical state. Surveys in Geophysics26, 733–756.
     [Google Scholar]
  61. WeckmannU., RitterO. and HaakV.2003. A magnetotelluric study of the Damara Belt in Namibia. Physics of the Earth and Planetary Interiors138, 91–112.
     [Google Scholar]
  62. WeideltP.1999. 3D conductivity models: Implications of electrical anisotropy. In: Three‐Dimensional Electromagnetics, (eds M.Oristaglio and B.Spies ), pp. 119–137. Society of Exploration Geophysicists Publishers.
     [Google Scholar]
  63. WeissC.J. and NewmanG.A.2002. Electromagnetic induction in a fully 3‐D anisotropic earth. Geophysics67, 1104–1114.
     [Google Scholar]
  64. WeissC.J. and NewmanG.A.2003. Electromagnetic induction in a generalized 3D anisotropic earth, Part 2: The LIN preconditioner. Geophysics68, 922–930.
     [Google Scholar]
  65. WuC., XuT., BadalJ., WuZ. and TengJ.2015. Seismic anisotropy across the Kunlun fault and their implications for northward transforming lithospheric deformation in northeastern Tibet. Tectonophysics659, 91–101.
     [Google Scholar]
  66. XiaoQ., ShaoG., YuG., CaiJ. and WangJ.2016. Electrical resistivity structures of the Kunlun–Qaidam–Qilian system at the northern Tibet and their tectonic implications. Physics of the Earth and Planetary Interiors255, 1–17.
     [Google Scholar]
  67. YangJ., XuZ., SongS., ZhangJ., WuC., ShiR., et al. 2001. Discovery of coesite in the north Qaidam early Palaeozoic ultrahigh pressure (UHP) metamorphic belt, NW China. Sciences de la Terre et des planètes333, 719–724.
     [Google Scholar]
  68. YeeK.S.1966. Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media. IEEE Transactions on Antennas & Propagation14, 302–307.
     [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12690
Loading
Three‐dimensional magnetotelluric responses for arbitrary electrically anisotropic media and a practical application
Geophysical Prospecting 66, 1764 (2018); https://doi.org/10.1111/1365-2478.12690
/content/journals/10.1111/1365-2478.12690
/content/journals/10.1111/1365-2478.12690
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Arbitrary electrical anisotropy; Finite difference; Magnetotelluric sounding; Northern Tibet; Three‐dimensional modelling

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error