Automated seismic event location by waveform coherence analysis

SUMMARY Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterized bylowsignal-to-noiseratio,automated locationmethods arerequested to be noise robust and sufficiently accurate. Most of the standard automated location routines arebasedontheautomatedpicking,identificationandassociationofthefirstarrivalsofPandS waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this paper, we propose a picking free earthquake location method based on the use of the short-term-average/long-term-average (STA/LTA) traces at different stations as observed data. For the P phases, we use the STA/LTA traces of the vertical energy function, whereas for the S phases, we use the STA/LTA traces of a second characteristic function, which is obtained using the principal component analysis technique. In order to locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a 3-D grid, we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial coordinates of the seismic event. A pilot application was performed in the Campania-Lucania region (southern Italy) using a seismic network (Irpinia Seismic Network) with an aperture of about 150km. We located 196 crustal earthquakes (depth < 20km) with magnitude range 1.1 < ML < 2.7. A subset of these locations were compared with accurate manual locations refined by using a double-difference technique. Our results indicate a good agreement with manual locations. Moreover, our method is noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals.

[1]  George A. McMechan,et al.  Determination of source parameters by wavefield extrapolation , 1982 .

[2]  Francesca Pacor,et al.  Ground-Motion Simulations for the 1980 M 6.9 Irpinia Earthquake (Southern Italy) and Scenario Events , 2011 .

[3]  E. Tessmer,et al.  Reverse modelling for seismic event characterization , 2005 .

[4]  Volker Oye,et al.  Orientation of three-component geophones in the San Andreas Fault observatory at depth Pilot Hole, Parkfield, California , 2005 .

[5]  Christopher John Young,et al.  An automated local and regional seismic event detection and location system using waveform correlation , 1999 .

[6]  Robert Granat,et al.  Real-time Earthquake Location Using Kirchhoff Reconstruction , 2005 .

[7]  R. V. Allen,et al.  Automatic phase pickers: Their present use and future prospects , 1982 .

[8]  F. Waldhauser,et al.  A Double-Difference Earthquake Location Algorithm: Method and Application to the Northern Hayward Fault, California , 2000 .

[9]  T. A. Stabile,et al.  Anatomy of a microearthquake sequence on an active normal fault , 2012, Scientific Reports.

[10]  A Cichowicz,et al.  An automatic S-phase picker , 1993, Bulletin of the Seismological Society of America.

[11]  Claudio Satriano,et al.  Automatic Picker Developments and Optimization: A Strategy for Improving the Performances of Automatic Phase Pickers , 2012 .

[12]  M. Weber,et al.  The effect of low-velocity sediments on the mislocation vectors of the GRF array , 1992 .

[13]  Lion Krischer,et al.  ObsPy: A Python Toolbox for Seismology , 2010 .

[14]  Göran Ekström,et al.  Global Detection and Location of Seismic Sources by Using Surface Waves , 2006 .

[15]  Robert E. Sheriff,et al.  Encyclopedic dictionary of applied geophysics , 2002 .

[16]  Christopher John Young,et al.  A comparison of select trigger algorithms for automated global seismic phase and event detection , 1998, Bulletin of the Seismological Society of America.

[17]  Vincenzo Convertito,et al.  A prototype system for earthquake early-warning and alert management in southern Italy , 2010 .

[18]  Manfred Baer,et al.  An automatic phase picker for local and teleseismic events , 1987 .

[19]  Paul Michaels,et al.  Use of Principal Component Analysis to Determine Down-Hole Tool Orientation and Enhance SH-Waves , 2001 .

[20]  Christopher John Young,et al.  An Automatic, Adaptive Algorithm for Refining Phase Picks in Large Seismic Data Sets , 2002 .

[21]  Francesco Grigoli,et al.  Automated Seismic Event Location by Travel‐Time Stacking: An Application to Mining Induced Seismicity , 2013 .

[22]  H. Kao,et al.  Delineating complex spatiotemporal distribution of earthquake aftershocks: an improved Source-Scanning Algorithm , 2012 .

[23]  Aldo Zollo,et al.  Short Note S-Wave Identification by Polarization Filtering and Waveform Coherence Analyses , 2012 .

[24]  C. Satriano,et al.  Automatic Picker Developments and Optimization: FilterPicker—a Robust, Broadband Picker for Real-Time Seismic Monitoring and Earthquake Early Warning , 2012 .

[25]  V. Oye,et al.  Automated microearthquake location using envelope stacking and robust global optimization , 2010 .

[26]  T. A. Stabile,et al.  Fault Delineation and Regional Stress Direction from the Analysis of Background Microseismicity in the southern Apennines, Italy , 2012 .

[27]  George A. McMechan,et al.  Two methods for determining geophone orientations from VSP data , 2006 .

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

[29]  Rongjiang Wang,et al.  A simple orthonormalization method for stable and efficient computation of Green's functions , 1999, Bulletin of the Seismological Society of America.

[30]  Frederik Tilmann,et al.  Coalescence microseismic mapping , 2013 .

[31]  Matthias Ohrnberger,et al.  Tracking the rupture of the Mw = 9.3 Sumatra earthquake over 1,150 km at teleseismic distance , 2005, Nature.

[32]  John E. Vidale,et al.  Complex polarization analysis of particle motion , 1986 .

[33]  G. Beroza,et al.  Full waveform earthquake location : Application to seismic streaks on the Calaveras Fault, California , 2007 .

[34]  M. Becquey,et al.  Three-component sonde orientation in a deviated well , 1990 .

[35]  H. Kao,et al.  The Source‐Scanning Algorithm: mapping the distribution of seismic sources in time and space , 2004 .

[36]  A. Zollo,et al.  An improved 1-D seismic velocity model for seismological studies in the Campania–Lucania region (Southern Italy) , 2013 .

[37]  P. Favali,et al.  Ocean bottom seismometers deployed in Tyrrhenian Sea , 2002 .

[38]  Tobias Diehl,et al.  Automatic S-Wave Picker for Local Earthquake Tomography , 2009 .

[39]  Matthias Ohrnberger,et al.  Spatio‐temporal source characteristics of the 26 December 2004 Sumatra earthquake as imaged by teleseismic broadband arrays , 2005 .

[40]  Gaetano Festa,et al.  Twin ruptures grew to build up the giant 2011 Tohoku, Japan, earthquake , 2012, Scientific Reports.

[41]  Shawn Larsen,et al.  Elastic modeling initiative, part III: 3-D computational modeling , 1998 .

[42]  M. Bischoff,et al.  Coal Mining Induced Seismicity in the Ruhr Area, Germany , 2010 .

[43]  Stewart Greenhalgh,et al.  Orientation of a downhole triaxial geophone , 1995 .

[44]  Martin Hensch.  On the interrelation of fluid-induced seismicity and crustal deformation at the Columbo submarine volcano (Aegean Sea, Greece) , 2010 .

[45]  L. Geiger.  Herdbestimmung bei Erdbeben aus den Ankunftszeiten , 1910 .

[46]  Kenneth S. Miller,et al.  Complex Linear Least Squares , 1973 .

[47]  Jianxin Yuan,et al.  Geophone orientation and coupling in three‐component sea‐floor data: a case study , 1999 .

[48]  Y. Nakamura,et al.  Large‐offset seismic surveying using ocean‐bottom seismographs and air guns: Instrumentation and field technique , 1987 .

[49]  P. Podvin,et al.  Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools , 1991 .

[50]  Joachim Wassermann,et al.  Locating the sources of volcanic explosions and volcanic tremor at Stromboli volcano (Italy) using beam-forming on diffraction hyperboloids , 1997 .

[51]  H. Kao,et al.  FAST TRACK PAPER: Rapid identification of earthquake rupture plane using Source-Scanning Algorithm , 2007 .

[52]  Aldo Zollo,et al.  Real-time evolutionary earthquake location for seismic early warning , 2008 .

[53]  Vincenzo Convertito,et al.  An Advanced Seismic Network in the Southern Apennines (Italy) for Seismicity Investigations and Experimentation with Earthquake Early Warning , 2007 .

[54]  Kelly B. Knowlton,et al.  Polarization measurement uncertainty on three-component VSP , 1996 .

[55]  H. D. Leslie,et al.  Automated Microseismic Event Detection and Location by Continuous Spatial Mapping , 2005 .