Characterizing Lithospheric Transverse Isotropy via Non-double Couple Components of Moment Tensors

The seismic moment tensor, which represents the equivalent body-force system of the seismic source (Backus and Mulcahy, 1973), may exhibit non-double couple components (NDCs) when the earthquake occurs on a planer fault if the source medium is anisotropic (Aki and Richards, 1981; Kawasaki and Tanimoto, 1981). Kawakatsu (1991, GRL) reported that the NDCs of the moment tensors for shallow earthquakes from the Harvard CMT catalog (Dziewonski et al.,1981; predecessor of GCMT) exhibit a systematic characteristic dependent on faulting types. Specifically, the sign of NDC on average systematically switches between normal-faulting and reverse-faulting. The average NDC parameter ε (Giardini, 1983) is negative for thrust faulting and positive for normal faulting. This behavior can be explained if the source region is transversely isotropic with a vertical symmetry axis (radially anisotropic). In fact, the transverse isotropy model of PREM at a depth of 24.4 km predicts the observed systematic NDC pattern, although the magnitude is slightly underestimated, indicating the potential to enhance our understanding of the lithospheric transverse isotropy using the NDC of the moment tensors.

To investigate the lithospheric transverse isotropy structure utilizing the NDCs of the moment tensors, we propose a novel inversion scheme, building upon the approaches employed by Vavrycuk (2004) and Li, Zheng, et al. (2018) for deep and intermediate-depth earthquakes, but with necessary modifications to address shallow sources (Kawakatsu, 1996, GJI). Synthetic tests conducted under conditions of random faulting indicate the potential to constrain the S-wave anisotropy (ξ) and the fifth parameter (η_κ; Kawakatsu, 2016, GJI). However, in realistic scenarios where a predominant stress regime influences earthquake occurrence to limit the diversity of faulting types, a significant correlation between these two parameters is anticipated, especially in regional-scale cases. Preliminary application of this method to real data sourced from the GCMT catalog suggests that the lithospheric transverse isotropy of PREM serves as a suitable initial model. However, some adjustments may be necessary, particularly regarding the fifth parameter, to enhance the model's fidelity in representing observed NDCs of the moment tensors.