Three-dimensional magnetotelluric inversion considering electrical anisotropy with synthetic and real data

Electrical anisotropy increases the non-uniqueness of magnetotelluric (MT) inversions. However, it is still useful to develop MT anisotropic inversion to provide new knowledge for this field. In this study, we developed a three-dimensional (3-D) MT inversion code considering electrical anisotropy. The forward procedure is based on a staggered-grid finite-difference method and solves the governing equation for the electric fields with a direct method. The values of a traditional objective function and its gradient could then be calculated. The inversion procedure was completed with a well-developed subroutine of non-linear conjugate gradient method. Synthetic tests for different isotropic models indicate the validity of this inversion code. Furthermore, we performed 3-D anisotropic inversions for synthetic and real data. The results indicate that vertical resistivity is still indistin-guishable in the anisotropic inversion results. In our synthetic tests, it becomes more difficult to recover low-resistivity components of anisotropic anomalies as the depth increases, and the vertical sizes of anisotropic anomalies are reduced, as characterized by a much shallower bottom boundary. The anisotropic strike angle is not as stable as the principal resistivities in the inversion process. Under the same changes in model parameters, the strike angle may cause larger deviations of the responses than the principal horizontal resistivities. With this code, we revisited the field data collected at the northern margin of the Tibetan Plateau. The inversion results indicate a possible anisotropic anomaly exists below the boundary between the plateau and its northern neighbor. This anisotropic anomaly could represent an area of concentration of shear deformation.