Estimation of the traditional transverse electric (TE) and transverse magnetic (TM) impedances of the magnetotelluric tensor for two-dimensional structures can be decoupled from the estimation of the strike direction with significant implications when dealing with galvanic distortions. Distortion-free data are obtainable by combining a quadratic equation with the phase tensor. In the terminology of Groom–Bailey, the quadratic equation provides amplitudes and phases that are immune to twist, and the phase tensor provides phases immune to both, twist and shear. On the other hand, distortion-free strike directions can be obtained using Bahr's approach or the phase tensor. In principle, this is all that is needed to proceed to a two-dimensional (2D) interpretation. However, the resulting impedances are strike ignorant because they are invariant under coordinate system rotation, and if they are to be related to a geological strike, they must be linked to a particular direction. This is an additional ambiguity to the one of 90° arising in classic strike-determination methods, which must be resolved independently. In this work, we use the distortion model of Groom–Bailey to resolve the ambiguity by bringing back the coupling between impedances and strike in the presence of galvanic distortions. Our approach is a hybrid between existing numerical and analytical methods that reduces the problem to a binary decision, which involves associating the invariant impedances with the correct TE and TM modes. To determine the appropriate association, we present three algorithms. Two of them require optimizing the fit to the data, and the third one requires a comparison of phases. All three keep track of possible crossings of the phase curves providing a clear-cut solution. Synthetic and field data illustrate the performance of the three schemes.Graphical
Read full abstract