LOC-FLOW: An End-to-End Machine Learning-Based High-Precision Earthquake Location Workflow

Abstract

Abstract The ever-increasing networks and quantity of seismic data drive the need for seamless and automatic workflows for rapid and accurate earthquake detection and location. In recent years, machine learning (ML)-based pickers have achieved remarkable accuracy and efficiency with generalization, and thus can significantly improve the earthquake location accuracy of previously developed sequential location methods. However, the inconsistent input or output (I/O) formats between multiple packages often limit their cross application. To reduce format barriers, we incorporated a widely used ML phase picker—PhaseNet—with several popular earthquake location methods and developed a “hands-free” end-to-end ML-based location workflow (named LOC-FLOW), which can be applied directly to continuous waveforms and build high-precision earthquake catalogs at local and regional scales. The renovated open-source package assembles several sequential algorithms including seismic first-arrival picking (PhaseNet and STA/LTA), phase association (REAL), absolute location (VELEST and HYPOINVERSE), and double-difference relative location (hypoDD and GrowClust). We provided different location strategies and I/O interfaces for format conversion to form a seamless earthquake location workflow. Different algorithms can be flexibly selected and/or combined. As an example, we apply LOC-FLOW to the 28 September 2004 Mw 6.0 Parkfield, California, earthquake sequence. LOC-FLOW accomplished seismic phase picking, association, velocity model updating, station correction, absolute location, and double-difference relocation for 16-day continuous seismic data. We detected and located 3.7 times (i.e., 4357) as many as earthquakes with cross-correlation double-difference locations from the Northern California Earthquake Data Center. Our study demonstrates that LOC-FLOW is capable of building high-precision earthquake catalogs efficiently and seamlessly from continuous seismic data.