An eikonal equation-based earthquake location method by inversion of multiple phase arrivals

Abstract

The precise determination of earthquake location is the fundamental basis in seismological community, and is crucial for analyzing seismic activity and performing seismic tomography. First arrivals are generally used to practically determine earthquake locations. However, first-arrival traveltimes are not sensitive to focal depths. Moreover, they cannot accurately constrain focal depths. To improve the accuracy, researchers have analyzed the depth phases of earthquake locations. The traveltimes of depth phases are sensitive to focal depths, and the joint inversion of depth phases and direct phases can be implemented to potentially obtain accurate earthquake locations. Generally, researchers can determine earthquake locations in layered models. Because layered models can only represent the first-order feature of subsurface structures, the advantages of joint inversion are not fully explored if layered models are used. To resolve the issue of current joint inversions, we use the traveltimes of three seismic phases to determine earthquake locations in heterogeneous models. The three seismic phases used in this study are the first P-, sPg- and PmP-waves. We calculate the traveltimes of the three seismic phases by solving an eikonal equation with an upwind difference scheme and use the traveltimes to determine earthquake locations. To verify the accuracy of the earthquake location method by the inversion of three seismic phases, we take the 2021 M_S6.4 Yangbi, Yunnan earthquake as an example and locate this earthquake using synthetic and real seismic data. Numerical tests demonstrate that the eikonal equation-based earthquake location method, which involves the inversion of multiple phase arrivals, can effectively improve earthquake location accuracy.