{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T01:46:46Z","timestamp":1760147206218,"version":"build-2065373602"},"reference-count":60,"publisher":"MDPI AG","issue":"2","license":[{"start":{"date-parts":[[2023,1,16]],"date-time":"2023-01-16T00:00:00Z","timestamp":1673827200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"National Key Research and Development Program of China","award":["2021YFB3202104","42030806","42074120","42274093","42174167","20204BCJL23058"],"award-info":[{"award-number":["2021YFB3202104","42030806","42074120","42274093","42174167","20204BCJL23058"]}]},{"name":"National Natural Science Foundation of China","award":["2021YFB3202104","42030806","42074120","42274093","42174167","20204BCJL23058"],"award-info":[{"award-number":["2021YFB3202104","42030806","42074120","42274093","42174167","20204BCJL23058"]}]},{"name":"Jiangxi Provincial Academic Leaders (Youth) Training Program","award":["2021YFB3202104","42030806","42074120","42274093","42174167","20204BCJL23058"],"award-info":[{"award-number":["2021YFB3202104","42030806","42074120","42274093","42174167","20204BCJL23058"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Remote Sensing"],"abstract":"<jats:p>The traditional three-dimensional (3D) magnetotelluric (MT) forward modeling using Krylov subspace algorithms has the problem of low modeling efficiency. To improve the computational efficiency of 3D MT forward modeling, we present a novel geometric multigrid algorithm for the finite element method. We use the vector finite element to discretize Maxwell\u2019s equations in the frequency domain and apply the Dirichlet boundary conditions to obtain large sparse complex linear equations for the solution of EM responses. To improve the convergence of the solution at low frequencies we use the divergence correction to correct the electric field. Then, we develop a V-cycle geometric multigrid algorithm to solve the linear equations system. To demonstrate the efficiency and effectiveness of our geometric multigrid method, we take three synthetic models (COMMEMI 3D-2 model, Dublin test model 1, modified SEG\/EAEG salt dome model) and compare our results with the published ones. Numerical results show that the geometric multigrid algorithm proposed in this paper is much better than the commonly used Krylov subspace algorithms (such as SOR-GMRES, ILU-BICGSTAB, SOR-BICGSTAB) in terms of the iteration number, the solution time, and the stability, and thus is more suitable for large-scale 3D MT forward modeling.<\/jats:p>","DOI":"10.3390\/rs15020537","type":"journal-article","created":{"date-parts":[[2023,1,17]],"date-time":"2023-01-17T02:58:16Z","timestamp":1673924296000},"page":"537","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":5,"title":["A Geometric Multigrid Method for 3D Magnetotelluric Forward Modeling Using Finite-Element Method"],"prefix":"10.3390","volume":"15","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-9436-6190","authenticated-orcid":false,"given":"Xianyang","family":"Huang","sequence":"first","affiliation":[{"name":"College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China"}]},{"given":"Changchun","family":"Yin","sequence":"additional","affiliation":[{"name":"College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China"}]},{"given":"Luyuan","family":"Wang","sequence":"additional","affiliation":[{"name":"College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-3634-9832","authenticated-orcid":false,"given":"Yunhe","family":"Liu","sequence":"additional","affiliation":[{"name":"College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China"}]},{"given":"Bo","family":"Zhang","sequence":"additional","affiliation":[{"name":"College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China"}]},{"given":"Xiuyan","family":"Ren","sequence":"additional","affiliation":[{"name":"College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China"}]},{"given":"Yang","family":"Su","sequence":"additional","affiliation":[{"name":"College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China"}]},{"given":"Jun","family":"Li","sequence":"additional","affiliation":[{"name":"College of Geo-Exploration Sciences and Technology, Jilin University, Changchun 130026, China"}]},{"given":"Hui","family":"Chen","sequence":"additional","affiliation":[{"name":"State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang 330013, China"}]}],"member":"1968","published-online":{"date-parts":[[2023,1,16]]},"reference":[{"key":"ref_1","first-page":"295","article-title":"On determining electric characteristics of the deep layers of the Earth\u2019s crust","volume":"73","author":"Tikhonov","year":"1950","journal-title":"Doklady"},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"605","DOI":"10.1190\/1.1437915","article-title":"Basic theory of the magneto-telluric method of geophysical prospecting","volume":"18","author":"Cagniard","year":"1953","journal-title":"Geophysics"},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"450","DOI":"10.1016\/j.jappgeo.2008.02.002","article-title":"Three-dimensional inversion of magnetotelluric data for mineral exploration: An example from the McArthur River uranium deposit, Saskatchewan, Canada","volume":"68","author":"Farquharson","year":"2009","journal-title":"J. Appl. Geophys."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"1628","DOI":"10.1093\/gji\/ggac029","article-title":"Application of multiscale magnetotelluric data to mineral exploration: An example from the east Tennant region, Northern Australia","volume":"229","author":"Jiang","year":"2022","journal-title":"Geophys. J. Int."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"23","DOI":"10.1016\/j.jappgeo.2014.04.003","article-title":"Theoretical assessment of 3-D magnetotelluric method for oil and gas exploration: Synthetic examples","volume":"106","author":"Zhang","year":"2014","journal-title":"J. Appl. Geophys."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"1005","DOI":"10.1007\/s10712-017-9439-x","article-title":"Magnetotelluric studies for hydrocarbon and geothermal resources: Examples from the Asian region","volume":"38","author":"Patro","year":"2017","journal-title":"Surv. Geophys."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"269","DOI":"10.1016\/j.jappgeo.2008.12.001","article-title":"Magnetotelluric and aeromagnetic investigations for assessment of groundwater resources in Parnaiba basin in Piaui State of North-East Brazil","volume":"68","author":"Chandrasekhar","year":"2009","journal-title":"J. Appl. Geophys."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"763","DOI":"10.1111\/j.1365-246X.2006.02948.x","article-title":"Error propagation in electromagnetic transfer functions: What role for the magnetotelluric method in detecting earthquake precursors?","volume":"165","author":"Hanekop","year":"2006","journal-title":"Geophys. J. Int."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"290","DOI":"10.1016\/j.epsl.2013.12.026","article-title":"Deep electrical resistivity structure of the northwestern US derived from 3-D inversion of USArray magnetotelluric data","volume":"402","author":"Meqbel","year":"2014","journal-title":"Earth Planet. Sci. Lett."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"1716","DOI":"10.1190\/1.1442984","article-title":"Advances in three-dimensional magnetotelluric modeling using integral equations","volume":"56","author":"Wannamaker","year":"1991","journal-title":"Geophysics"},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"L08307","DOI":"10.1029\/2011GL046953","article-title":"Three-dimensional inversion of large-scale EarthScope magnetotelluric data based on the integral equation method: Geoelectrical imaging of the Yellowstone conductive mantle plume","volume":"38","author":"Zhdanov","year":"2011","journal-title":"Geophys. Res. Lett."},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"215","DOI":"10.1190\/1.1443407","article-title":"Three-dimensional magnetotelluric modeling using difference equations\u2014Theory and comparisons to integral equation solutions","volume":"58","author":"Mackie","year":"1993","journal-title":"Geophysics"},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"923","DOI":"10.1029\/94RS00326","article-title":"Three-dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example","volume":"29","author":"Mackie","year":"1994","journal-title":"Radio Sci."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"40","DOI":"10.1016\/j.cageo.2014.01.010","article-title":"ModEM: A modular system for inversion of electromagnetic geophysical data","volume":"66","author":"Kelbert","year":"2014","journal-title":"Comput. Geosci."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"1943","DOI":"10.1137\/S1064827599360741","article-title":"Fast finite volume simulation of 3D electromagnetic problems with highly discontinuous coefficients","volume":"22","author":"Haber","year":"2001","journal-title":"SIAM J. Sci. Comput."},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"E287","DOI":"10.1190\/geo2013-0312.1","article-title":"A finite-volume solution to the geophysical electromagnetic forward problem using unstructured grids","volume":"79","author":"Jahandari","year":"2014","journal-title":"Geophysics"},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"700","DOI":"10.1093\/gji\/ggt154","article-title":"A goal-oriented adaptive finite-element approach for plane wave 3-D electromagnetic modelling","volume":"194","author":"Ren","year":"2013","journal-title":"Geophys. J. Int."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"556","DOI":"10.1007\/s11770-018-0703-8","article-title":"3D magnetotelluric inversions with unstructured finite-element and limited-memory quasi-Newton methods","volume":"15","author":"Cao","year":"2018","journal-title":"Appl. Geophys."},{"key":"ref_19","unstructured":"Heagy, L.J., Capriotii, J., Kuttai, J., Cowan, D., Perez, F., Hamman, J., Banihirwe, A., and Paul, K. (2020, January 7\u201311). Advances in Magnetotelluric modelling and inversion with SimPEG. Proceedings of the AGU Fall Meeting Abstracts, Online."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"105030","DOI":"10.1016\/j.cageo.2021.105030","article-title":"3D magnetotelluric modeling using high-order tetrahedral N\u00e9d\u00e9lec elements on massively parallel computing platforms","volume":"160","author":"Modesto","year":"2022","journal-title":"Comput. Geosci."},{"key":"ref_21","unstructured":"Jin, J.-M. (2015). The Finite Element Method in Electromagnetics, John Wiley & Sons."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"E101","DOI":"10.1190\/geo2010-0398.1","article-title":"A finite element multifrontal method for 3D CSEM modeling in the frequency domain","volume":"77","author":"Morgan","year":"2012","journal-title":"Geophysics"},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"615","DOI":"10.1093\/gji\/ggt006","article-title":"Finite element based inversion for time-harmonic electromagnetic problems","volume":"193","author":"Schwarzbach","year":"2013","journal-title":"Geophys. J. Int."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"74","DOI":"10.1093\/gji\/ggv410","article-title":"3-D magnetotelluric inversion including topography using deformed hexahedral edge finite elements and direct solvers parallelized on SMP computers\u2013Part I: Forward problem and parameter Jacobians","volume":"204","author":"Kordy","year":"2016","journal-title":"Geophys. J. Int."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"E337","DOI":"10.1190\/geo2015-0580.1","article-title":"A goal-oriented adaptive finite-element method for 3D scattered airborne electromagnetic method modeling","volume":"81","author":"Yin","year":"2016","journal-title":"Geophysics"},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"WB71","DOI":"10.1190\/geo2017-0544.1","article-title":"Adaptive finite element for 3D time-domain airborne electromagnetic modeling based on hybrid posterior error estimation","volume":"83","author":"Zhang","year":"2018","journal-title":"Geophysics"},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"5912813","DOI":"10.1109\/TGRS.2022.3153026","article-title":"3D Finite-Element Forward Modeling of Airborne EM Systems in Frequency-Domain Using Octree Meshes","volume":"60","author":"Han","year":"2022","journal-title":"IEEE Trans. Geosci. Remote Sens."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"1189","DOI":"10.1093\/gji\/ggx358","article-title":"3-D minimum-structure inversion of magnetotelluric data using the finite-element method and tetrahedral grids","volume":"211","author":"Jahandari","year":"2017","journal-title":"Geophys. J. Int."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"204","DOI":"10.1016\/j.jappgeo.2018.08.013","article-title":"Direct solutions of 3-D magnetotelluric fields using edge-based finite element","volume":"159","author":"Xiong","year":"2018","journal-title":"J. Appl. Geophys."},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"170","DOI":"10.1007\/s11770-008-0024-4","article-title":"Three-dimensional magnetotellurics modeling using edgebased finite-element unstructured meshes","volume":"5","author":"Liu","year":"2008","journal-title":"Appl. Geophys."},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"1460","DOI":"10.1093\/gji\/ggt071","article-title":"Efficient pre-conditioned iterative solution strategies for the electromagnetic diffusion in the Earth: Finite-element frequency-domain approach","volume":"193","author":"Um","year":"2013","journal-title":"Geophys. J. Int."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"678","DOI":"10.1093\/gji\/ggt027","article-title":"A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling","volume":"193","author":"Puzyrev","year":"2013","journal-title":"Geophys. J. Int."},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"164","DOI":"10.1016\/j.cageo.2014.09.008","article-title":"3D controlled-source electromagnetic modeling in anisotropic medium using edge-based finite element method","volume":"73","author":"Cai","year":"2014","journal-title":"Comput. Geosci."},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1016\/j.jappgeo.2017.12.009","article-title":"Three-dimensional magnetotelluric modeling in anisotropic media using edge-based finite element method","volume":"149","author":"Xiao","year":"2018","journal-title":"J. Appl. Geophys."},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"699","DOI":"10.1016\/j.jappgeo.2011.09.025","article-title":"Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction","volume":"75","author":"Farquharson","year":"2011","journal-title":"J. Appl. Geophys."},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"104324","DOI":"10.1016\/j.jappgeo.2021.104324","article-title":"Three-dimensional edge-based finite element modeling of magnetotelluric data in anisotropic media with a divergence correction","volume":"189","author":"Zhou","year":"2021","journal-title":"J. Appl. Geophys."},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"E367","DOI":"10.1190\/geo2020-0520.1","article-title":"Three-dimensional controlled-source electromagnetic forward modeling by edge-based finite element with a divergence correction","volume":"86","author":"Tang","year":"2021","journal-title":"Geophysics"},{"key":"ref_38","unstructured":"Trottenberg, U., Oosterlee, C.W., and Schuller, A. (2000). Multigrid, Elsevier."},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"633","DOI":"10.1111\/j.1365-2478.2006.00558.x","article-title":"A multigrid solver for 3D electromagnetic diffusion","volume":"54","author":"Mulder","year":"2006","journal-title":"Geophys. Prospect."},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"129","DOI":"10.1007\/s00791-007-0064-y","article-title":"Geophysical modelling of 3D electromagnetic diffusion with multigrid","volume":"11","author":"Mulder","year":"2008","journal-title":"Comput. Vis. Sci."},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"1554","DOI":"10.1093\/gji\/ggw352","article-title":"Fully anisotropic 3-D EM modelling on a Lebedev grid with a multigrid pre-conditioner","volume":"207","author":"Jaysaval","year":"2016","journal-title":"Geophys. J. Int."},{"key":"ref_42","first-page":"334","article-title":"Stable electromagnetic modeling using a multigrid solver on stretching grids: The magnetotelluric example","volume":"13","author":"Li","year":"2016","journal-title":"IEEE Geosci. Remote Sens. Lett."},{"key":"ref_43","doi-asserted-by":"crossref","first-page":"E121","DOI":"10.1190\/geo2021-0275.1","article-title":"An efficient multigrid solver based on a four-color cell-block Gauss-Seidel smoother for 3D magnetotelluric forward modeling","volume":"87","author":"Guo","year":"2022","journal-title":"Geophysics"},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"1834","DOI":"10.1093\/gji\/ggac152","article-title":"An efficient cascadic multigrid solver for 3-D magnetotelluric forward modelling problems using potentials","volume":"230","author":"Pan","year":"2022","journal-title":"Geophys. J. Int."},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"1442","DOI":"10.1093\/gji\/ggu086","article-title":"Algebraic multigrid preconditioning within parallel finite-element solvers for 3-D electromagnetic modelling problems in geophysics","volume":"197","author":"Koldan","year":"2014","journal-title":"Geophys. J. Int."},{"key":"ref_46","doi-asserted-by":"crossref","first-page":"154","DOI":"10.1007\/s11770-017-0605-1","article-title":"Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method","volume":"14","author":"Chen","year":"2017","journal-title":"Appl. Geophys."},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"2073","DOI":"10.1093\/gji\/ggx511","article-title":"Electrical Resistivity Tomography using a finite element based BFGS algorithm with algebraic multigrid preconditioning","volume":"212","author":"Codd","year":"2018","journal-title":"Geophys. J. Int."},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"1798","DOI":"10.1007\/s11430-020-9786-9","article-title":"3D finite-element modeling of Earth induced electromagnetic field and its potential applications for geomagnetic satellites","volume":"64","author":"Yao","year":"2021","journal-title":"Sci. China Earth Sci."},{"key":"ref_49","doi-asserted-by":"crossref","unstructured":"Briggs, W.L., Henson, V.E., and McCormick, S.F. (2000). A Multigrid Tutorial, SIAM.","DOI":"10.1137\/1.9780898719505"},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"57","DOI":"10.1007\/BF01389668","article-title":"A new family of mixed finite elements in \u211d3","volume":"50","year":"1986","journal-title":"Numer. Math."},{"key":"ref_51","doi-asserted-by":"crossref","first-page":"1319","DOI":"10.1190\/1.1444055","article-title":"Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator","volume":"61","author":"Smith","year":"1996","journal-title":"Geophysics"},{"key":"ref_52","doi-asserted-by":"crossref","first-page":"1646","DOI":"10.1137\/S1064827596305829","article-title":"Flexible multiple semicoarsening for three-dimensional singularly perturbed problems","volume":"19","author":"Washio","year":"1998","journal-title":"SIAM J. Sci. Comput."},{"key":"ref_53","doi-asserted-by":"crossref","first-page":"431","DOI":"10.1016\/j.jcp.2007.03.033","article-title":"A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation","volume":"224","author":"Riyanti","year":"2007","journal-title":"J. Comput. Phys."},{"key":"ref_54","doi-asserted-by":"crossref","first-page":"1088","DOI":"10.1109\/TMAG.2009.2012627","article-title":"Robustness of nested multigrid method for edge-based finite element analysis","volume":"45","author":"Watanabe","year":"2009","journal-title":"IEEE Trans. Magn."},{"key":"ref_55","doi-asserted-by":"crossref","first-page":"133","DOI":"10.1016\/S0926-9851(97)00013-X","article-title":"Methods for modelling electromagnetic fields results from COMMEMI\u2014The international project on the comparison of modelling methods for electromagnetic induction","volume":"37","author":"Zhdanov","year":"1997","journal-title":"J. Appl. Geophys."},{"key":"ref_56","doi-asserted-by":"crossref","first-page":"1216","DOI":"10.1093\/gji\/ggt066","article-title":"Magnetotelluric 3-D inversion\u2014A review of two successful workshops on forward and inversion code testing and comparison","volume":"193","author":"Miensopust","year":"2013","journal-title":"Geophys. J. Int."},{"key":"ref_57","doi-asserted-by":"crossref","first-page":"721","DOI":"10.1186\/BF03351724","article-title":"Numerical accuracy of magnetotelluric modeling: A comparison of finite difference approximations","volume":"54","author":"Siripunvaraporn","year":"2002","journal-title":"Earth Planets Space"},{"key":"ref_58","doi-asserted-by":"crossref","unstructured":"Farquharson, C.G., Oldenburg, D.W., Haber, E., and Shekhtman, R. (2002). An algorithm for the three-dimensional inversion of magnetotelluric data. SEG Technical Program Expanded Abstracts 2002, Society of Exploration Geophysicists.","DOI":"10.1190\/1.1817336"},{"key":"ref_59","doi-asserted-by":"crossref","first-page":"277","DOI":"10.1111\/j.1365-2478.2007.00614.x","article-title":"3D magnetotelluric modelling including surface topography","volume":"55","author":"Nam","year":"2007","journal-title":"Geophys. Prospect."},{"key":"ref_60","doi-asserted-by":"crossref","first-page":"949","DOI":"10.1190\/1.1437054","article-title":"SEG\/EAEG 3-D modeling project: 2nd update","volume":"13","author":"Aminzadeh","year":"1994","journal-title":"Lead. Edge"}],"container-title":["Remote Sensing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2072-4292\/15\/2\/537\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,10]],"date-time":"2025-10-10T18:07:41Z","timestamp":1760119661000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2072-4292\/15\/2\/537"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,1,16]]},"references-count":60,"journal-issue":{"issue":"2","published-online":{"date-parts":[[2023,1]]}},"alternative-id":["rs15020537"],"URL":"https:\/\/doi.org\/10.3390\/rs15020537","relation":{},"ISSN":["2072-4292"],"issn-type":[{"type":"electronic","value":"2072-4292"}],"subject":[],"published":{"date-parts":[[2023,1,16]]}}}