High-precision earthquake locations in Switzerland using regional secondary arrivals in a 3-D velocity model

We present a new approach to relocate earthquakes in the greater western Alpine region using main crustal phases ( Pg , Pn , PmP ) that takes advantage of recent developments in P -wave velocity models and modelling of the Moho topography in the region, as well as the ability to track reflected and refracted phases in three-dimensional (3-D) heterogeneous media. Our approach includes a new 3-D P -wave velocity model for Switzerland and surrounding regions that combines a first-order Moho discontinuity based on local earthquake tomography (LET) and controlled-source seismology (CSS) information and 3-D seismic velocity information based on LET. Traveltimes for the main crustal phases ( Pg , Pn , PmP ) are computed using a fast marching method. We use a non-linear, probabilistic approach to relocate earthquakes that has been extended to include the use of secondary phases. We validate our approach using synthetic data, which was computed for a real earthquake and different combinations of available phases ( Pg , Pn , PmP ). We also applied our approach to relocate four selected earthquakes, two shallow and two deep crustal events in the northern Alpine foreland, for which independent information (ground truth information) on their focal depths exist. Our results demonstrate that the precision and accuracy of focal depth estimates can be greatly improved if secondary phases are used. This gain is a combined effect of an improved range of take-off angles and the use of differential traveltimes between first and secondary arriving phases. Our results also show that reliable information on the Moho depth is crucial to obtain accurate focal depths, if Pn or PmP phases are used in the relocation process. Finally, our approach demonstrates that proper identification of the main crustal phases in combination with an appropriate model parametrization in the forward solver will significantly improve earthquake locations.

[1]  Chao-ying Bai,et al.  2-D/3-D irregular shortest-path ray tracing for multiple arrivals and its applications , 2010 .

[2]  Walter H. F. Smith,et al.  New, improved version of generic mapping tools released , 1998 .

[3]  Donat Fäh,et al.  Earthquakes in Switzerland and surrounding regions during 2007 , 2008 .

[4]  C. Thurber,et al.  Seismic velocity structure and event relocation in Kazakhstan from secondary P phases , 1992, Bulletin of the Seismological Society of America.

[5]  Donat Fäh,et al.  Earthquakes in Switzerland and surrounding regions during 2006 , 2007 .

[6]  N. Deichmann,et al.  High-resolution 3-D P-wave model of the Alpine crust , 2009 .

[7]  Anthony Lomax,et al.  A Reanalysis of the Hypocentral Location and Related Observations for the Great 1906 California Earthquake , 2005 .

[8]  Urs Kradolfer,et al.  Initial reference models in local earthquake tomography , 1994 .

[9]  Jean Virieux,et al.  Probabilistic Earthquake Location in 3D and Layered Models , 2000 .

[10]  Keith L. McLaughlin,et al.  Collection of a Reference Event Set for Regional and Teleseismic Location Calibration , 2004 .

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

[12]  E. Kissling.  Deep structure of the Alps—what do we really know? , 1993 .

[13]  D. Giardini,et al.  Earthquakes in Switzerland and surrounding regions during 2005 , 2006 .

[14]  D. Giardini,et al.  Earthquakes in Switzerland and Surrounding Regions During 1999 , 2000 .

[15]  Felix Waldhauser,et al.  Three dimensional interface modelling with two-dimensional seismic data: the Alpine crust-mantle boundary , 1998 .

[16]  Tobias Diehl,et al.  Consistent phase picking for regional tomography models: application to the greater Alpine region , 2009 .

[17]  Joan S. Gomberg,et al.  The effect of S-wave arrival times on the accuracy of hypocenter estimation , 1990, Bulletin of the Seismological Society of America.

[18]  P. Molnar,et al.  Focal depths of intracontinental and intraplate earthquakes and their implications for the thermal and mechanical properties of the lithosphere , 1983 .

[19]  S. Parolai,et al.  Anomalously deep earthquakes in northwestern Italy , 1999 .

[20]  M. Sambridge,et al.  Multiple reflection and transmission phases in complex layered media using a multistage fast marching method , 2004 .

[21]  K. Priestley,et al.  Earthquake focal depths, effective elastic thickness, and the strength of the continental lithosphere , 2000 .

[22]  Donat Fäh,et al.  Earthquakes in Switzerland and surrounding regions during 2004 , 2005 .

[23]  C. Thurber,et al.  Advances in seismic event location , 2000 .

[24]  E. Kissling,et al.  Model parametrization in seismic tomography: a choice of consequence for the solution quality , 2001 .

[25]  J. Luetgert.  MacRay; interactive two-dimensional seismic raytracing for the Macintosh , 1992 .

[26]  R. A. Bailey,et al.  The earth's crust , 2002 .

[27]  F. Waldhauser,et al.  High-resolution teleseismic tomography of upper-mantle structure using an a priori three-dimensional crustal model , 2002 .

[28]  S. Schmid,et al.  The arc of the western Alps in the light of geophysical data on deep crustal structure , 2000 .

[29]  Xiaoping Yang,et al.  Location Calibration Data for CTBT Monitoring at the Prototype International Data Center , 2001 .

[30]  E. Engdahl,et al.  Global teleseismic earthquake relocation with improved travel times and procedures for depth determination , 1998, Bulletin of the Seismological Society of America.

[31]  E. Kissling,et al.  Investigating effects of 3-D ray tracing methods in local earthquake tomography , 2001 .

[32]  D. Giardini,et al.  Evidence for an active fault below the northwestern Alpine foreland of Switzerland , 2007 .

[33]  Domenico Giardini,et al.  Earthquakes Induced By the Stimulation of an Enhanced Geothermal System Below Basel (Switzerland) , 2009 .

[34]  H. Kao,et al.  Hypocentre determination offshore of eastern Taiwan using the Maximum Intersection method , 2004 .

[35]  Malcolm Sambridge,et al.  A practical grid-based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous media , 2006 .

[36]  K. Priestley,et al.  Lithospheric structure and deep earthquakes beneath India, the Himalaya and southern Tibet , 2008 .

[37]  M. Sambridge,et al.  The Fast Marching Method: An Effective Tool for Tomographic Imaging and Tracking Multiple Phases in Complex Layered Media , 2005 .

[38]  S. Mueller.  A New Model of the Continental Crust , 2013 .

[39]  M. Sambridge,et al.  Wave front evolution in strongly heterogeneous layered media using the fast marching method , 2004 .

[40]  Peter Bormann,et al.  The IASPEI Standard Seismic Phase List , 2003 .

[41]  E. Kissling,et al.  Combining controlled-source seismology and local earthquake tomography to derive a 3-D crustal model of the western Alpine region , 2012 .

[42]  Shutian Ma.  Focal Depth Determination for Moderate and Small Earthquakes by Modeling Regional Depth Phases sPg, sPmP, and sPn , 2010 .

[43]  Nicholas Deichmann,et al.  Structural and rheological implications of lower-crustal earthquakes below northern Switzerland , 1992 .

[44]  Stefan Wiemer,et al.  Probabilistic earthquake location in complex three‐dimensional velocity models: Application to Switzerland , 2003 .

[45]  Constraining Event Depths and Crustal Velocities Using Regional Depth Phases , 2009 .

[46]  A. Michelini,et al.  High-resolution Rayleigh-wave velocity maps of central Europe from a dense ambient-noise data set , 2012 .