Statistical Approaches to Detecting Transient Signals in GPS: Results from the 2009–2011 Transient Detection Exercise

We present the results of our participation in four phases of the Southern California Earthquake Center (SCEC) transient detection exercise (Lohman and Murray, 2013). In each phase, a blind test was conducted in which sets of synthetic Global Positioning Systems (GPS) data were released and a deadline set for submission of detection results. For each data set, the presence or absence of transient events was to be determined, and the location and time of each specified. After all submissions were received, the ground‐truth information about any transient signals in the data was released. The synthetic data sets were generated by FAKENET, a software package that produces realistic GPS time series that include secular motion and seasonal signals as well as realistic noise and distributions of missing data (Agnew, 2013). Station locations in the synthetic data set were a subset of GPS installations in the western United States. In this work, we pursue a purely data‐driven approach to transient detection, rather than one based on an assumption of an underlying physical model. We view this approach as having two important advantages. First, it facilitates the detection of events, which might happen through a previously unknown, but geophysically interesting, physical process or on a previously unknown fault or structure. Second, it enables the detection of events, which have nothing to do with the solid earth but, although not the subject of the SCEC exercise, have scientific or practical merit. These include signals such as those due to atmospheric phenomena as well as signals, which result from hardware faults or failures and software processing glitches. The four phases of the exercise were conducted over approximately two years. As a result, our approach evolved as we learned from experiences in the initial phases as well as from other ongoing work during that period. The …

[1]  Frank H. Webb,et al.  Slow But Not Quite Silent , 2003, Science.

[2]  Duncan Carr Agnew Realistic Simulations of Geodetic Network Data: The Fakenet Package , 2013 .

[3]  Michael B. Heflin,et al.  The coseismic geodetic signature of the 1999 Hector Mine earthquake , 2000 .

[4]  Biing-Hwang Juang,et al.  Hidden Markov Models for Speech Recognition , 1991 .

[5]  Brian Kan-Wing Mak,et al.  Subspace distribution clustering hidden Markov model , 2001, IEEE Trans. Speech Audio Process..

[6]  A. Donnellan,et al.  Combined GPS and InSAR Models of Postseismic Deformation from the Northridge Earthquake , 2002 .

[7]  Marlon E. Pierce,et al.  Modeling and On-the-Fly Solutions for Solid Earth Sciences: Web Services and Data Portal for Earthquake Early Warning System , 2008, IGARSS 2008 - 2008 IEEE International Geoscience and Remote Sensing Symposium.

[8]  Gráinne McGuire,et al.  A Bayesian Model for Detecting Past Recombination Events in DNA Multiple Alignments , 2000, J. Comput. Biol..

[9]  Timothy Ian Melbourne,et al.  Precursory transient slip during the 2001 Mw = 8.4 Peru earthquake sequence from continuous GPS , 2002 .

[10]  Steve J. Young,et al.  State clustering in hidden Markov model-based continuous speech recognition , 1994, Comput. Speech Lang..

[11]  R. Lohman,et al.  The SCEC Geodetic Transient‐Detection Validation Exercise , 2013 .

[12]  Tim Melbourne,et al.  Periodic Slow Earthquakes from the Cascadia Subduction Zone , 2002, Science.

[13]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[14]  Galip Aydin,et al.  Analysis of streaming GPS measurements of surface displacement through a web services environment , 2007, 2007 IEEE Symposium on Computational Intelligence and Data Mining.

[15]  Yehuda Bock,et al.  Spatiotemporal filtering using principal component analysis and Karhunen-Loeve expansion approaches for regional GPS network analysis , 2006 .

[16]  Biing-Hwang Juang,et al.  Mixture autoregressive hidden Markov models for speech signals , 1985, IEEE Trans. Acoust. Speech Signal Process..

[17]  Lawrence R. Rabiner,et al.  A minimum discrimination information approach for hidden Markov modeling , 1989, IEEE Trans. Inf. Theory.

[18]  Andrea Donnellan,et al.  GPS observations of fault afterslip and upper crustal deformation following the Northridge earthquake , 1998 .

[19]  Y. Okada Surface deformation due to shear and tensile faults in a half-space , 1985 .

[20]  Hiromichi Tsuji,et al.  Silent fault slip following an interplate thrust earthquake at the Japan Trench , 1997, Nature.

[21]  Xuelong Li,et al.  Gabor-Based Region Covariance Matrices for Face Recognition , 2008, IEEE Transactions on Circuits and Systems for Video Technology.

[22]  Kai-Fu Lee,et al.  Context-independent phonetic hidden Markov models for speaker-independent continuous speech recognition , 1990 .

[23]  S. Miyazaki,et al.  A slow thrust slip event following the two 1996 Hyuganada Earthquakes beneath the Bungo Channel, southwest Japan , 1999 .

[24]  L. Rivera,et al.  Coseismic Deformation from the 1999 Mw 7.1 Hector Mine, California, Earthquake as Inferred from InSAR and GPS Observations , 2002 .

[25]  J. Parker,et al.  Analysis and modeling of Southern California deformation , 2000 .

[26]  Dirk P. Kroese,et al.  Kernel density estimation via diffusion , 2010, 1011.2602.

[27]  Kenneth W. Hudnut,et al.  THE SOUTHERN CALIFORNIA INTEGRATED GPS NETWORK (SCIGN) , 2001 .

[28]  Yehuda Koren,et al.  Factor in the neighbors: Scalable and accurate collaborative filtering , 2010, TKDD.

[29]  Claire Monteleoni,et al.  Tracking climate models , 2011, CIDU.

[30]  Fatih Murat Porikli,et al.  Human Detection via Classification on Riemannian Manifolds , 2007, 2007 IEEE Conference on Computer Vision and Pattern Recognition.

[31]  Jerome R. Bellegarda,et al.  Tied mixture continuous parameter modeling for speech recognition , 1990, IEEE Trans. Acoust. Speech Signal Process..

[32]  Hsiao-Wuen Hon,et al.  Speaker-independent phone recognition using hidden Markov models , 1989, IEEE Trans. Acoust. Speech Signal Process..

[33]  A. Farago,et al.  Algorithm to find the global optimum of left-to-right hidden Markov model parameters , 1989 .

[34]  D. Haussler,et al.  Hidden Markov models in computational biology. Applications to protein modeling. , 1993, Journal of molecular biology.

[35]  Andrea Donnellan,et al.  Influence of anelastic surface layers on postseismic thrust fault deformation , 2000 .

[36]  R. Polikar,et al.  Ensemble based systems in decision making , 2006, IEEE Circuits and Systems Magazine.

[37]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[38]  Timothy Ian Melbourne,et al.  Rapid postseismic transients in subduction zones from continuous GPS , 2002 .

[39]  Kenneth W. Hudnut,et al.  Detection of crustal deformation from the Landers earthquake sequence using continuous geodetic measurements , 1993, Nature.

[40]  Fatih Murat Porikli,et al.  Pedestrian Detection via Classification on Riemannian Manifolds , 2008, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[41]  H. Dragert,et al.  Episodic Tremor and Slip on the Cascadia Subduction Zone: The Chatter of Silent Slip , 2003, Science.

[42]  Lawrence R. Rabiner,et al.  A tutorial on hidden Markov models and selected applications in speech recognition , 1989, Proc. IEEE.