Abstract

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.

Highlights

  • Time reversal (TR) is a universal concept that can be found in numerous physical sciences, including meteorology (e.g. Talagrand and Courtier, 2007), geodynamics (e.g. Bunge et al, 2003), ground water modelling (e.g. Sun, 1994) and seismology

  • 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

  • Kremers 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

Read more

Summary

Introduction

Time reversal (TR) is a universal concept that can be found in numerous physical sciences, including meteorology (e.g. Talagrand and Courtier, 2007), geodynamics (e.g. Bunge et al, 2003), ground water modelling (e.g. Sun, 1994) and seismology. The time-invariance of perfectly elastic wave propagation provides the intuitive justification for the TR imaging of seismic sources: Seismograms u0(xr ,t) recorded at positions xr (r = 1,...,n) are reversed in time, re-injected as sources at their respective receiver locations and the resulting wave field u(x,t) is . Tarantola, 1988; Tromp et al, 2004; Fichtner et al, 2006; Fichtner, 2010) leads to a similar result: The gradient of the misfit χ with respect to the source parameters is given in terms of the time-reversed wave field generated by adjoint sources that radiate the misfit from the receiver positions back into the Earth model. We provide an application to the strong-motion data recorded during the 2000 Tottori (Japan) earthquake

Numerical method
Single point source
Multiple point sources
Quantitative assessment of focussing for point sources
Synthetic finite source simulations
Dense regular grid of stations
Station arrays
Weighting of adjoint sources
Time-reverse imaging of the Tottori earthquake source
Discussion
Large-amplitude surface waves
The missing sink
Lack of prior information
Incomplete knowledge of the 3-D Earth structure
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call