Earthquake Phase Association using a Bayesian Gaussian Mixture Model

TL;DR

This paper introduces GaMMA, a Bayesian Gaussian Mixture Model-based method for earthquake phase association that effectively clusters seismic phases and estimates earthquake parameters without requiring supervised training, especially useful for dense seismic data.

Contribution

The paper presents a novel unsupervised clustering approach using Gaussian mixture models for phase association and earthquake parameter estimation, eliminating the need for grid-search or supervised training.

Findings

GaMMA accurately associates phases in dense seismic sequences.

It provides reliable estimates of earthquake location and magnitude.

Effective on both synthetic and real earthquake data.

Abstract

Earthquake phase association algorithms aggregate picked seismic phases from a network of seismometers into individual earthquakes and play an important role in earthquake monitoring. Dense seismic networks and improved phase picking methods produce massive earthquake phase data sets, particularly for earthquake swarms and aftershocks occurring closely in time and space, making phase association a challenging problem. We present a new association method, the Gaussian Mixture Model Association (GaMMA), that combines the Gaussian mixture model for phase measurements (both time and amplitude), with earthquake location, origin time, and magnitude estimation. We treat earthquake phase association as an unsupervised clustering problem in a probabilistic framework, where each earthquake corresponds to a cluster of P and S phases with hyperbolic moveout of arrival times and a decay of amplitude with distance. We use a multivariate Gaussian distribution to model the collection of phase picks for an event, the mean of which is given by the predicted arrival time and amplitude from the causative event. We carry out the pick assignment for each earthquake and determine earthquake parameters (i.e., earthquake location, origin time, and magnitude) under the maximum likelihood criterion using the Expectation-Maximization (EM) algorithm. The GaMMA method does not require the typical association steps of other algorithms, such as grid-search or supervised training. The results on both synthetic test and the 2019 Ridgecrest earthquake sequence show that GaMMA effectively associates phases from a temporally and spatially dense earthquake sequence while producing useful estimates of earthquake location and magnitude.