Exploring the potentials and limitations of the time-reversal imaging of finite seismic sources

Abstract. The characterisation of seismic sources with time-reversed wave fields is developing into a standard technique that has already been successful in numerous applications. While the time-reversal imaging of effective point sources is now well-understood, little work has been done to extend this technique to the study of finite rupture processes. This is despite the pronounced non-uniqueness in classic finite source inversions. The need to better constrain the details of finite rupture processes motivates the series of synthetic and real-data time reversal experiments described in this paper. We address questions concerning the quality of focussing in the source area, the localisation of the fault plane, the estimation of the slip distribution and the source complexity up to which time-reversal imaging can be applied successfully. The frequency band for the synthetic experiments is chosen such that it is comparable to the band usually employed for finite source inversion. Contrary to our expectations, we find that time-reversal imaging is useful only for effective point sources, where it yields good estimates of both the source location and the origin time. In the case of finite sources, however, the time-reversed field does not provide meaningful characterisations of the fault location and the rupture process. This result cannot be improved sufficiently with the help of different imaging fields, realistic modifications of the receiver geometry or weights applied to the time-reversed sources. The reasons for this failure are manifold. They include the choice of the frequency band, the incomplete recording of wave field information at the surface, the excitation of large-amplitude surface waves that deteriorate the depth resolution, the absence of a sink that should absorb energy radiated during the later stages of the rupture process, the invisibility of small slip and the neglect of prior information concerning the fault geometry and the inherent smoothness of seismologically inferred Earth models that prevents the beneficial occurrence of strong multiple-scattering. The condensed conclusion of our study is that the limitations of time-reversal imaging – at least in the frequency band considered here – start where the seismic source stops being effectively point-localised.

[1]  G. Festa Earthquake Source Inversion Blindtest: Initial Results and Further Developments , 2007 .

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

[3]  Michel Campillo,et al.  Frequency domain inversion of strong motions: Application to the 1992 Landers earthquake , 1995 .

[4]  P. Courtier,et al.  Variational Assimilation of Meteorological Observations With the Adjoint Vorticity Equation. I: Theory , 2007 .

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

[6]  William S. Hodgkiss,et al.  An initial demonstration of underwater acoustic communication using time reversal , 2002 .

[7]  Andreas Fichtner,et al.  Full seismic waveform tomography for upper-mantle structure in the Australasian region using adjoint methods , 2009 .

[8]  C. S. Clay,et al.  Reproducibility of signal transmissions in the ocean , 1965 .

[9]  P. Spudich,et al.  A global search inversion for earthquake kinematic rupture history: Application to the 2000 western Tottori, Japan earthquake , 2007 .

[10]  Andreas Fichtner,et al.  Simulation and Inversion of Seismic Wave Propagation on Continental Scales Based on a Spectral - Element Method , 2009 .

[11]  Mathias Fink,et al.  Time‐reversal imaging of seismic sources and application to the great Sumatra earthquake , 2006 .

[12]  Michel Campillo,et al.  The 2000 Tottori earthquake: A shallow earthquake with no surface rupture and slip properties controlled by depth , 2005 .

[13]  Qinya Liu,et al.  Tomography, Adjoint Methods, Time-Reversal, and Banana-Doughnut Kernels , 2004 .

[14]  A. Fichtner,et al.  Efficient numerical surface wave propagation through the optimization of discrete crustal models—a technique based on non-linear dispersion curve matching (DCM) , 2008 .

[15]  P. Shearer,et al.  A High-Frequency Secondary Event During the 2004 Parkfield Earthquake , 2007, Science.

[16]  V. Hjorleifsdottir Earthquake source characterization using 3D numerical modeling , 2007 .

[17]  Alexander Sutin,et al.  Single-channel time reversal in elastic solids , 2004 .

[18]  M. Fink,et al.  Time reversal processing in ultrasonic nondestructive testing , 1995, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[19]  T. Dahm,et al.  Automated procedure for point and kinematic source inversion at regional distances , 2010 .

[20]  M. Fink,et al.  Time reversal acoustics , 1996, IEEE Ultrasonics Symposium, 2004.

[21]  George A. McMechan,et al.  Imaging of earthquake sources in Long Valley Caldera, California, 1983 , 1985 .

[22]  P. Drummond,et al.  Time reversed acoustics , 1997 .

[23]  Brad Artman,et al.  Source location using time‐reverse imaging , 2010 .

[24]  Mickael Tanter,et al.  Multiwave imaging and super resolution , 2010 .

[25]  P. Courtier,et al.  Variational Assimilation of Meteorological Observations With the Adjoint Vorticity Equation. Ii: Numerical Results , 2007 .

[26]  Andreas Fichtner,et al.  The adjoint method in seismology – I. Theory , 2006 .

[27]  B. Kennett Geophysics: A new way to estimate seismic source parameters , 1983, Nature.

[28]  Peter M. Shearer,et al.  Extent, duration and speed of the 2004 Sumatra–Andaman earthquake imaged by the Hi-Net array , 2005, Nature.

[29]  F. Cotton,et al.  Frequency Domain Inversion of Strong Motions: Application to the Blind Test for Kinematic Source Inversion (EC-project SPICE) , 2007 .

[30]  Albert Tarantola,et al.  Theoretical background for the inversion of seismic waveforms including elasticity and attenuation , 1988 .

[31]  F. Scherbaum,et al.  Acoustic imaging of earthquake sources from the Chalfant Valley, 1986, aftershock series , 1994 .

[32]  Gilberto Saccorotti,et al.  Time reverse location of seismic long-period events recorded on Mt Etna , 2011 .

[33]  Erik H. Saenger,et al.  Time reverse modeling of low‐frequency microtremors: Application to hydrocarbon reservoir localization , 2008 .

[34]  R. Guyer,et al.  Tremor source location using time reversal: Selecting the appropriate imaging field , 2009 .

[35]  A. Fichtner Full Seismic Waveform Modelling and Inversion , 2011 .

[36]  C. R. Hagelberg,et al.  Mantle circulation models with variational data assimilation: inferring past mantle flow and structure from plate motion histories and seismic tomography , 2001 .

[37]  Roux,et al.  Robust Acoustic Time Reversal with High-Order Multiple Scattering. , 1995, Physical review letters.

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

[39]  M. Dumbser,et al.  An arbitrary high-order Discontinuous Galerkin method for elastic waves on unstructured meshes — III. Viscoelastic attenuation , 2007 .

[40]  M. Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - I. The two-dimensional isotropic case with external source terms , 2006 .