Iterative tomographic analysis based on automatic refined picking

The ever-growing size of data sets for active and passive seismic imaging makes the availability of automatic procedures for rapid analysis more and more valuable. Such procedures are especially important for time-critical applications like emergency decisions or re-orienting of ongoing seismic surveys. In this paper a new, iterative scheme for 3D traveltime tomography is presented. The technique, based on a tool originally developed for earthquake data, uses cross-correlation to examine waveform similarity and to adjust arrival times on seismic sections. A preliminary set of reference arrival times is first corrected by the cross-correlation lag and then used to build an initial 3D tomographic velocity model through a standard inversion code; traveltimes calculated from this model are then taken as new reference arrivals and the process of pick adjustment is repeated. The result is a tomographic image, upgraded and refined at each iteration of the procedure. The test performed on the waveform data set recorded during the 2001 SERAPIS active seismic survey in the gulfs of Naples and Pozzuoli (Southern Italy) shows that the 3D iterative tomography scheme produces a velocity image of the structure of the Campi Flegrei caldera which is consistent with the results from previous studies, employing just a fraction of the time needed by a human analyst to identify first breaks. We believe that this technique can be effectively employed for rapid analysis of large data-sets within time-critical or time-dependent tasks and for automatic 4D tomographic investigations.

[1]  Aldo Zollo,et al.  Seismic images and rock properties of the very shallow structure of Campi Flegrei caldera (southern Italy) , 2009 .

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

[3]  A. Zollo,et al.  The Bay of Naples (southern Italy): Constraints on the volcanic structures inferred from a dense seismic survey , 2004 .

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

[5]  Clifford H. Thurber,et al.  Hypocenter-velocity structure coupling in local earthquake tomography , 1992 .

[6]  Elijah Polak,et al.  Computational methods in optimization , 1971 .

[7]  P. Shearer.  Application to the Whittier Narrows California aftershock sequence , 1997 .

[8]  D. Schmitt,et al.  First-break timing: Arrival onset times by direct correlation , 1999 .

[9]  Gregory C. Beroza,et al.  Foreshock sequence of the 1992 Landers, California, earthquake and its implications for earthquake nucleation , 1995 .

[10]  P. Gori,et al.  Time-Resolved Seismic Tomography Detects Magma Intrusions at Mount Etna , 2006, Science.

[11]  Fred W. Klein,et al.  Deep fault plane geometry inferred from multiplet relative relocation beneath the south flank of Kilauea , 1994 .

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

[13]  T. Lay,et al.  Modern Global Seismology , 1995 .

[14]  B. Gelchinsky,et al.  AUTOMATIC PICKING OF FIRST ARRIVALS AND PARAMETERIZATION OF TRAVELTIME CURVES , 1983 .

[15]  Reinoud Sleeman,et al.  Robust automatic P-phase picking: an on-line implementation in the analysis of broadband seismogram recordings , 1999 .

[16]  Clifford H. Thurber,et al.  Dome growth behavior at Soufriere Hills Volcano, Montserrat, revealed by relocation of volcanic event swarms, 1995-1996 , 2004 .

[17]  Stephen D. Malone,et al.  High precision relative locations of earthquakes at Mount St. Helens, Washington , 1987 .

[18]  Luca D'Auria,et al.  Evidence for the buried rim of Campi Flegrei caldera from 3‐d active seismic imaging , 2003 .

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

[20]  F. Boschetti,et al.  A fractal-based algorithm for detecting first arrivals on seismic traces , 1996 .

[21]  John C. Lahr,et al.  Three‐dimensional P and S wave velocity structure of Redoubt Volcano, Alaska , 1996 .

[22]  Clifford H. Thurber,et al.  Double-Difference Tomography: The Method and Its Application to the Hayward Fault, California , 2003 .

[23]  Michael E. Murat,et al.  AUTOMATED FIRST ARRIVAL PICKING: A NEURAL NETWORK APPROACH1 , 1992 .

[24]  Clifford H. Thurber,et al.  Automatic P-Wave Arrival Detection and Picking with Multiscale Wavelet Analysis for Single-Component Recordings , 2003 .

[25]  Jean Virieux,et al.  An efficient algorithm for double‐difference tomography and location in heterogeneous media, with an application to the Kilauea volcano , 2005 .

[26]  Hiroaki Niitsuma,et al.  Multiplet-Clustering Analysis Reveals Structural Details within the Seismic Cloud at the Soultz Geothermal Field, France , 2003 .