Strain Green's Tensors, Reciprocity, and Their Applications to Seismic Source and Structure Studies

Green's function approach is widely used in modeling seismic wave- form. The representation theorem expresses the wave field as the inner product of the moment tensor and the spatial gradients of the Green's tensor. Standard practice in waveform calculations has been to compute the Green's tensors first and then obtain their gradients by numerical differentiation. The reciprocity of the Green's tensor enables us to express the wave field explicitly in terms of the strain Green's tensor, a third-order tensor composed of the spatial gradients of the Green's tensor elements. We propose here to use the strain Green's tensors rather than the Green's tensors themselves in computing the waveforms. By bypassing the need for Green's tensors and directly using the strain Green's tensors, we can improve the computa- tional efficiency in waveform modeling while eliminating the possible errors from numerical differentiation. The strain Green's tensor elements are also directly related to the partial derivatives of the waveforms with respect to moment tensor elements and structural parameters. Through the inversion of the focal mechanisms of 27 small events in the Los Angeles region, we demonstrate the effectiveness of the strain Green's tensor database approach in quickly recovering source parameters based on realistic 3D models. We show that the same database can also be used to improve the efficiency and accuracy in computing the Frechet kernels for tomography inver- sions.

[1]  Luis Rivera,et al.  A note on the dynamic and static displacements from a point source in multilayered media , 2002 .

[2]  D. Komatitsch,et al.  Simulations of Ground Motion in the Los Angeles Basin Based upon the Spectral-Element Method , 2004 .

[3]  Robert W. Graves,et al.  The SCEC Southern California Reference Three-Dimensional Seismic Velocity Model Version 2 , 2000 .

[4]  Kim B. Olsen,et al.  Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models , 2005 .

[5]  Masayuki Kikuchi,et al.  Inversion of complex body waves , 1982 .

[6]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[7]  Toshiro Tanimoto,et al.  Waveforms of long-period body waves in a slightly aspherical earth model , 1993 .

[8]  Jeroen Tromp,et al.  Spectral-element moment tensor inversions for earthquakes in Southern California , 2004 .

[9]  Barbara Romanowicz,et al.  Toward real-time estimation of regional moment tensors , 1996, Bulletin of the Seismological Society of America.

[10]  T. Jordan,et al.  Seismic structure of the upper mantle in a central Pacific corridor , 1996 .

[11]  Michel Bouchon,et al.  Teleseismic body wave radiation from a seismic source in a layered medium , 1976 .

[12]  J. Shaw P-wave seismic velocity structure derived from sonic logs and industry reflection data in the Los An , 2003 .

[13]  Peter M. Shearer,et al.  A New Method for Determining First-Motion Focal Mechanisms , 2002 .

[14]  Thomas H. Jordan,et al.  Three‐dimensional Fréchet differential kernels for seismicdelay times , 2000 .

[15]  Monica D. Kohler,et al.  Mantle Heterogeneities and the SCEC Reference Three-Dimensional Seismic Velocity Model Version 3 , 2003 .

[16]  Lupei Zliu,et al.  Moho depth variation in southern California from teleseismic receiver functions , 2000 .

[17]  Robert W. Clayton,et al.  A Reciprocity Method for Multiple-Source Simulations , 2001 .

[18]  T. Jordan,et al.  Structural sensitivities of finite-frequency seismic waves: a full-wave approach , 2006 .

[19]  Thomas H. Jordan,et al.  Generalized seismological data functionals , 1992 .

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

[21]  Guust Nolet,et al.  Wavefront healing: a banana–doughnut perspective , 2001 .

[22]  Guoqing Lin,et al.  Southern California Hypocenter Relocation with Waveform Cross- Correlation, Part 2: Results Using Source-Specific Station Terms and Cluster Analysis , 2005 .

[23]  Thomas H. Jordan,et al.  Teleseismic inversion for the second-degree moments of earthquake space-time distributions , 2001 .

[24]  Thomas H. Jordan,et al.  Finite-Moment Tensor of the 3 September 2002 Yorba Linda Earthquake , 2005 .

[25]  Donald V. Helmberger,et al.  Advancement in source estimation techniques using broadband regional seismograms , 1996, Bulletin of the Seismological Society of America.

[26]  Hiroo Kanamori,et al.  Moho depth variation in southern California from teleseismic receiver functions , 2000 .

[27]  Robert W. Graves,et al.  Modeling path effects in three-dimensional basin structures , 1992, Bulletin of the Seismological Society of America.

[28]  Robert W. Graves,et al.  Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions. 1. Strong motions , 2001 .

[29]  3-D Sensitivity Kernels for Surface Wave Observables , 2002 .

[30]  Masayuki Kikuchi,et al.  Inversion of complex body waves—III , 1991, Bulletin of the Seismological Society of America.

[31]  Guust Nolet,et al.  Three‐dimensional sensitivity kernels for surface wave observables , 2004 .

[32]  Martin J. Siegert,et al.  EOS Trans. AGU , 2003 .

[33]  T. Jordan,et al.  High‐resolution, two‐dimensional vertical tomography of the central Pacific mantle using ScS reverberations and frequency‐dependent travel times , 1998 .

[34]  Thorne Lay,et al.  Quantitative Seismology, Second Edition , 2003 .

[35]  Egill Hauksson,et al.  Crustal structure and seismicity distribution adjacent to the Pacific and North America plate boundary in southern California , 2000 .

[36]  Guust Nolet,et al.  Fréchet kernels for finite-frequency traveltimes—I. Theory , 2000 .