An efficient decoupled 3-D axial anisotropic resistivity inversion for magnetotelluric data with OpenMP parallelization
Abstract We have developed two novel axial anisotropic inversion codes for magnetotelluric (MT) data: a full axial inversion and a decoupled axial inversion. Both codes are based on the data space Gauss–Newton inversion method used in the WSINV3DMT code, enhanced with OpenMP parallelization and Inte...
Saved in:
| Main Authors: | , , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
SpringerOpen
2025-01-01
|
| Series: | Earth, Planets and Space |
| Subjects: | |
| Online Access: | https://doi.org/10.1186/s40623-024-02113-5 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1832585967309946880 |
|---|---|
| author | Wisart Thongyoy Weerachai Siripunvaraporn Puwis Amatayakul Tawat Rung-Arunwan Kiyoshi Baba |
| author_facet | Wisart Thongyoy Weerachai Siripunvaraporn Puwis Amatayakul Tawat Rung-Arunwan Kiyoshi Baba |
| author_sort | Wisart Thongyoy |
| collection | DOAJ |
| description | Abstract We have developed two novel axial anisotropic inversion codes for magnetotelluric (MT) data: a full axial inversion and a decoupled axial inversion. Both codes are based on the data space Gauss–Newton inversion method used in the WSINV3DMT code, enhanced with OpenMP parallelization and Intel MKL PARDISO direct solver for improved computational efficiency. The full axial inversion searches for three axial anisotropic resistivity elements $${\rho }_{xx}$$ ρ xx , $${\rho }_{yy}$$ ρ yy , and $${\rho }_{zz}$$ ρ zz . Based on our prior influence studies, we identified contrasting influence patterns: $${\rho }_{xx}$$ ρ xx has a strong influence on Z xy and Z yy , whereas $${\rho }_{yy}$$ ρ yy has a strong influence on Z yx and Z xx , while $${\rho }_{zz}$$ ρ zz has negligible impact on all responses. We, therefore, proposed a decoupled axial inversion. This novel technique splits the full axial inversion into two independent modes: $${\rho }_{xx}$$ ρ xx - and $${\rho }_{yy}$$ ρ yy -modes. The $${\rho }_{xx}$$ ρ xx -mode requires just the Z xy and Z yy responses to recover $${\rho }_{xx}$$ ρ xx , while the $${\rho }_{yy}$$ ρ yy -mode requires the other two responses. By splitting the inversion, each mode requires half the data parameters and reduces model parameters to one-third compared to the full axial inversion. The substantial parameter reduction is the key factor leading to significantly faster processing times and lower memory requirements. Both codes were validated with the synthetic model. Utilizing OpenMP and direct solvers, the full axial inversion processed 256 MT stations in 65 min per iteration or a total of 195 min for three iterations to converge. The decoupled inversion achieved significantly faster processing, requiring just 32 min per iteration or a total of 96 min to finish the inversion. In addition, the decoupled inversion requires just a fraction of the memory used in the full axial inversion. This makes it practical to operate even on the standard personal computers of current technology. We further applied the decoupled inversion to real MT data acquired in the Northwestern Pacific Ocean. The resulting 3-D inverted $${\rho }_{xx}$$ ρ xx and $${\rho }_{yy}$$ ρ yy models support the existence of the anisotropy occurring at 60–120 km depth beneath the ocean floor, in agreement with previous studies. Graphical abstract |
| format | Article |
| id | doaj-art-14fad40d62f14c49852b29301ff877c1 |
| institution | Kabale University |
| issn | 1880-5981 |
| language | English |
| publishDate | 2025-01-01 |
| publisher | SpringerOpen |
| record_format | Article |
| series | Earth, Planets and Space |
| spelling | doaj-art-14fad40d62f14c49852b29301ff877c12025-01-26T12:19:33ZengSpringerOpenEarth, Planets and Space1880-59812025-01-0177112110.1186/s40623-024-02113-5An efficient decoupled 3-D axial anisotropic resistivity inversion for magnetotelluric data with OpenMP parallelizationWisart Thongyoy0Weerachai Siripunvaraporn1Puwis Amatayakul2Tawat Rung-Arunwan3Kiyoshi Baba4Department of Physics, Faculty of Science, Mahidol UniversityDepartment of Physics, Faculty of Science, Mahidol UniversityDepartment of Physics, Faculty of Science, Mahidol UniversityCurl-E Geophysics Co., Ltd.Earthquake Research Institute, The University of TokyoAbstract We have developed two novel axial anisotropic inversion codes for magnetotelluric (MT) data: a full axial inversion and a decoupled axial inversion. Both codes are based on the data space Gauss–Newton inversion method used in the WSINV3DMT code, enhanced with OpenMP parallelization and Intel MKL PARDISO direct solver for improved computational efficiency. The full axial inversion searches for three axial anisotropic resistivity elements $${\rho }_{xx}$$ ρ xx , $${\rho }_{yy}$$ ρ yy , and $${\rho }_{zz}$$ ρ zz . Based on our prior influence studies, we identified contrasting influence patterns: $${\rho }_{xx}$$ ρ xx has a strong influence on Z xy and Z yy , whereas $${\rho }_{yy}$$ ρ yy has a strong influence on Z yx and Z xx , while $${\rho }_{zz}$$ ρ zz has negligible impact on all responses. We, therefore, proposed a decoupled axial inversion. This novel technique splits the full axial inversion into two independent modes: $${\rho }_{xx}$$ ρ xx - and $${\rho }_{yy}$$ ρ yy -modes. The $${\rho }_{xx}$$ ρ xx -mode requires just the Z xy and Z yy responses to recover $${\rho }_{xx}$$ ρ xx , while the $${\rho }_{yy}$$ ρ yy -mode requires the other two responses. By splitting the inversion, each mode requires half the data parameters and reduces model parameters to one-third compared to the full axial inversion. The substantial parameter reduction is the key factor leading to significantly faster processing times and lower memory requirements. Both codes were validated with the synthetic model. Utilizing OpenMP and direct solvers, the full axial inversion processed 256 MT stations in 65 min per iteration or a total of 195 min for three iterations to converge. The decoupled inversion achieved significantly faster processing, requiring just 32 min per iteration or a total of 96 min to finish the inversion. In addition, the decoupled inversion requires just a fraction of the memory used in the full axial inversion. This makes it practical to operate even on the standard personal computers of current technology. We further applied the decoupled inversion to real MT data acquired in the Northwestern Pacific Ocean. The resulting 3-D inverted $${\rho }_{xx}$$ ρ xx and $${\rho }_{yy}$$ ρ yy models support the existence of the anisotropy occurring at 60–120 km depth beneath the ocean floor, in agreement with previous studies. Graphical abstracthttps://doi.org/10.1186/s40623-024-02113-5Electrical anisotropyMagnetotelluricInversionData spaceMarine electromagneticNorthwestern Pacific |
| spellingShingle | Wisart Thongyoy Weerachai Siripunvaraporn Puwis Amatayakul Tawat Rung-Arunwan Kiyoshi Baba An efficient decoupled 3-D axial anisotropic resistivity inversion for magnetotelluric data with OpenMP parallelization Earth, Planets and Space Electrical anisotropy Magnetotelluric Inversion Data space Marine electromagnetic Northwestern Pacific |
| title | An efficient decoupled 3-D axial anisotropic resistivity inversion for magnetotelluric data with OpenMP parallelization |
| title_full | An efficient decoupled 3-D axial anisotropic resistivity inversion for magnetotelluric data with OpenMP parallelization |
| title_fullStr | An efficient decoupled 3-D axial anisotropic resistivity inversion for magnetotelluric data with OpenMP parallelization |
| title_full_unstemmed | An efficient decoupled 3-D axial anisotropic resistivity inversion for magnetotelluric data with OpenMP parallelization |
| title_short | An efficient decoupled 3-D axial anisotropic resistivity inversion for magnetotelluric data with OpenMP parallelization |
| title_sort | efficient decoupled 3 d axial anisotropic resistivity inversion for magnetotelluric data with openmp parallelization |
| topic | Electrical anisotropy Magnetotelluric Inversion Data space Marine electromagnetic Northwestern Pacific |
| url | https://doi.org/10.1186/s40623-024-02113-5 |
| work_keys_str_mv | AT wisartthongyoy anefficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT weerachaisiripunvaraporn anefficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT puwisamatayakul anefficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT tawatrungarunwan anefficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT kiyoshibaba anefficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT wisartthongyoy efficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT weerachaisiripunvaraporn efficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT puwisamatayakul efficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT tawatrungarunwan efficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization AT kiyoshibaba efficientdecoupled3daxialanisotropicresistivityinversionformagnetotelluricdatawithopenmpparallelization |