Automatic recognition of T and teleseismic P waves by statistical analysis of their spectra: An application to continuous records of moored hydrophones

A network of moored hydrophones is an effective way of monitoring seismicity of oceanic ridges since it allows detection and localization of underwater events by recording generated T waves. The high cost of ship time necessitates long periods (normally a year) of autonomous functioning of the hydrophones, which results in very large data sets. The preliminary but indispensable part of the data analysis consists of identifying all T wave signals. This process is extremely time consuming if it is done by a human operator who visually examines the entire database. We propose a new method for automatic signal discrimination based on the Gradient Boosted Decision Trees technique that uses the distribution of signal spectral power among different frequency bands as the discriminating characteristic. We have applied this method to automatically identify the types of acoustic signals in data collected by two moored hydrophones in the North Atlantic. We show that the method is capable of efficiently resolving the signals of seismic origin with a small percentage of wrong identifications and missed events: 1.2% and 0.5% for T waves and 14.5% and 2.8% for teleseismic P waves, respectively. In addition, good identification rates for signals of other types (iceberg and ship generated) are obtained. Our results indicate that the method can be successfully applied to automate the analysis of other (not necessarily acoustic) databases provided that enough information is available to describe statistical properties of the signals to be identified.

[1]  Greg Ridgeway,et al.  Generalized Boosted Models: A guide to the gbm package , 2006 .

[2]  F. Waldhauser,et al.  Frequency‐magnitude distribution of microearthquakes beneath the 9°50′N region of the East Pacific Rise, October 2003 through April 2004 , 2008 .

[3]  W. Munk,et al.  Ocean Acoustic Tomograpthy: Notation , 1995 .

[4]  I. Daubechies,et al.  Biorthogonal bases of compactly supported wavelets , 1992 .

[5]  Guust Nolet,et al.  Modern mermaids: New floats image the deep Earth , 2011 .

[6]  F. Simons,et al.  Analysis of seafloor seismograms of the 2003 Tokachi‐Oki earthquake sequence for earthquake early warning , 2008 .

[7]  Colin MacBeth,et al.  The application of back‐propagation neural network to automatic picking seismic arrivals from single‐component recordings , 1997 .

[8]  S. Solomon,et al.  Microearthquake Characteristics of a Mid‐Ocean Ridge along‐axis high , 1992 .

[9]  R. Dziak,et al.  The June‐July 1993 seismo‐acoustic event at CoAxial segment, Juan de Fuca Ridge: Evidence for a lateral dike injection , 1995 .

[10]  Alberto Michelini,et al.  Automatic picking of P and S phases using a neural tree , 2006 .

[11]  Won Sang Lee,et al.  Tectonomagmatic activity and ice dynamics in the Bransfield Strait back-arc basin, Antarctica , 2010 .

[12]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

[13]  Anders la Cour-Harbo,et al.  Ripples in Mathematics , 2003 .

[14]  Mathilde Cannat,et al.  Regional seismicity of the Mid-Atlantic Ridge: observations from autonomous hydrophone arrays , 2010 .

[15]  Deborah K. Smith,et al.  Time-clustering behavior of spreading-center seismicity between 15 and 35 ◦ N on the Mid-Atlantic Ridge: observations from hydroacoustic monitoring , 2003 .

[16]  Alois Knoll,et al.  Gradient boosting machines, a tutorial , 2013, Front. Neurorobot..

[17]  E. Okal,et al.  Seismic and hydroacoustic tremor generated by colliding icebergs , 2008 .

[18]  R. Dziak,et al.  P wave detection thresholds, Pn velocity estimates, and T wave location uncertainty from oceanic hydrophones , 1999 .

[19]  R. Dziak.  Empirical relationship of T‐wave energy and fault parameters of northeast Pacific Ocean earthquakes , 2001 .

[20]  Deborah K. Smith,et al.  P- and T-Wave Detection Thresholds, Pn Velocity Estimate, and Detection of Lower Mantle and Core P-Waves on Ocean Sound-Channel Hydrophones at the Mid-Atlantic Ridge , 2004 .

[21]  J. Royer,et al.  Spatiotemporal distribution of the seismicity along the Mid‐Atlantic Ridge north of the Azores from hydroacoustic data: Insights into seismogenic processes in a ridge–hot spot context , 2012 .

[22]  E. Okal.  The generation of T waves by earthquakes , 2008 .

[23]  J. Friedman.  Greedy function approximation: A gradient boosting machine. , 2001 .

[24]  Yue Zhao,et al.  An artificial neural network approach for broadband seismic phase picking , 1999 .

[25]  Guust Nolet,et al.  On the potential of recording earthquakes for global seismic tomography by low‐cost autonomous instruments in the oceans , 2009 .

[26]  R. V. Allen,et al.  Automatic earthquake recognition and timing from single traces , 1978, Bulletin of the Seismological Society of America.

[27]  Y. Hello,et al.  Automatic discrimination of underwater acoustic signals generated by teleseismic P‐waves: A probabilistic approach , 2011 .

[28]  Deborah K. Smith,et al.  Hydroacoustic monitoring of seismicity at the slow‐spreading Mid‐Atlantic Ridge , 2002 .

[29]  E. Okal,et al.  Long-range detection of hydroacoustic signals from large icebergs in the Ross Sea, Antarctica , 2002 .

[30]  W. Sweldens.  The Lifting Scheme: A Custom - Design Construction of Biorthogonal Wavelets "Industrial Mathematics , 1996 .

[31]  Hadley Wickham,et al.  The Split-Apply-Combine Strategy for Data Analysis , 2011 .

[32]  J. Friedman.  Stochastic gradient boosting , 2002 .

[33]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[34]  G. Masters,et al.  Global P and PP traveltime tomography: rays versus waves , 2004 .

[35]  M. Tolstoy,et al.  Sound‐channel observations of ice‐generated tremor in the Indian Ocean , 2005 .

[36]  Hadley Wickham,et al.  Reshaping Data with the reshape Package , 2007 .

[37]  S. Solomon,et al.  Microearthquake characteristics and crustal velocity structure at 29°N on the Mid‐Atlantic Ridge: The architecture of a slow spreading segment , 1995 .

[38]  Deborah K. Smith,et al.  Aftershock sequences in the mid-ocean ridge environment: an analysis using hydroacoustic data , 2002 .

[39]  Christopher G. Fox,et al.  Monitoring Pacific Ocean seismicity from an autonomous hydrophone array , 2001 .

[40]  M. Ewing,et al.  The T phase of shallow-focus earthquakes , 1950 .

[41]  B. Cornuelle,et al.  Ocean acoustic tomography , 1988, Signal Recovery and Synthesis.