Abstract:At present, the 3-D magnetotelluric (MT) unstructured finite element modeling approaches often adopt the block-wise uniform conductivity parameterization method, in which the Earth conductivity model is discretized into plenty of tetrahedrons with uniform conductivity. However, that block-wise uniform parameterization method is not only difficult to accurately simulate actual continuous subsurface conductivity distribution, but also would substantially increase the unknowns in the subsequent inversion, which largely aggravates the non-uniqueness of inversion. In this study, we propose a 3-D unstructured finite element modeling strategy of MT problems allowing for continuous variation of conductivity in each tetrahedron. First, the curl–curl governing differential equation of electrical field for the 3-D MT problems and the corresponding integral weak form are derived successively. Next, the open-source tetrahedral mesh generator Tetgen is utilized to discretize the computational domain, both the conductivity and the electric field within each tetrahedral element are approximated by using linear interpolation with node and edge vector basis functions, respectively. Then, the closed-form expression of the resulting elemental integral is derived and the Newman-type boundary condition is loaded. The resulting system of linear equations is solved by using the direct LU decomposition method. Further, to enhance the accuracy of the solution, a goal-oriented adaptive technology based on the continuity of the normal component of the current density is adopted. Finally, the correctness of the proposed algorithm is verified by using a layered model with continuous conductivity distribution and a 3-D topographic model. While the necessity of considering the continuous variation of conductivity in blocks is proved by using two 3-D continuum models. The results show that the proposed algorithm can obviously improve the accuracy of the responses at the sites above the region with continuous conductivity distribution and will greatly reduce the unknowns of the subsequent inversion, compared with the traditional algorithm with block-wise uniform parameterization.

A Remote-Reference-Free Parametric Method for Robust Magnetotelluric Impedance Tensor Estimation Evaluated at Several USArray Sites

Parametric Magnetotelluric Impedance Tensor Estimation

An Evaluation of Robust Remote Reference and Parametric Magnetotelluric Transfer Function Estimation.

A Parametric Method for Robust Magnetotelluric Transfer Function Estimation Evaluated with Data at Different Sampling Frequencies

Parametric Electromagnetic Transfer Function Estimation at USArray Site MNF34

Two-dimensional Magnetotelluric Modelling by the Boundary Element Method.

Data-Driven Confidence Intervals for Parametric Magnetotelluric Impedance Tensor Estimates

A FEASIBILITY STUDY OF RAPID TWO-DIMENSIONAL MT CURVE COMPARISON INVERSION

Modeling the 3D Terrain Effect on MT by the Boundary Element Method

Analysis Of Integrating Ert Data Of Different Array Types: Simulation Study And Observations On An Archaeological Site

Comparison of Stable Maximum Likelihood Estimator with Traditional Robust Estimator in Magnetotelluric Impedance Estimation

Comparisons of the Magnetotelluric experimental results by Multi‐site remote reference and Multi‐site superposition

Reducing Reliance on Observation Duration of Magnetotelluric Impedance Estimation With an Improved Instantaneous Spectrum-Based Method

Robust inversion of 1D magnetotelluric data using the Huber loss function

3D magnetotelluric inversion using a limited-memory quasi-Newton optimization

Three-Dimensional Unstructured Finite Element Modeling of Magnetotelluric Problems Allowing for Continuous Variation of Conductivity in Each Block

An Alternating Direction Method of Multipliers Algorithm for 1-D Magnetotelluric Anisotropic Inversion Using Fourier Series Expansion

Appraisal of the Magnetotelluric and Magnetovariational Transfer Functions' Selection in a 3-D Inversion

Unstructured Grid Finite Element Modeling of the Three-Dimensional Magnetotelluric Responses in a Model With Arbitrary Conductivity and Magnetic Susceptibility Anisotropies

A limited-memory quasi-Newton inversion for 1D magnetotellurics

Estimation of the efficiency of unbiased predictive risk estimator in the inversion of 2D magnetotelluric data