Rapid Full-Wave Centroid Moment Tensor (CMT) Inversion in a Three-Dimensional Earth Structure Model for Earthquakes in Southern California

SUMMARY A central problem of seismology is the inversion of regional waveform data for models of earthquake sources. In regions such as Southern California, preliminary 3-D earth structure models are already available, and efficient numerical methods have been developed for 3-D anelastic wave-propagation simulations. We describe an automated procedure that utilizes these capabilities to derive centroid moment tensors (CMTs). The procedure relies on the use of receiver-side Green’s tensors (RGTs), which comprise the spatial-temporal displacements produced by the three orthogonal unit impulsive point forces acting at the receivers. We have constructed a RGT database for 219 broad-band stations in Southern California using a tomographically improved version of the 3-D SCEC Community Velocity Model Version 4.0 (CVM4) and a staggered-grid finite-difference code. Finite-difference synthetic seismograms for any earthquake in our modelling volume can be simply calculated by extracting a small, source-centred volume from the RGT database and applying the reciprocity principle. The partial derivatives needed for the CMT inversion can be generated in the same way. We have developed an automated algorithm that combines a grid-search for suitable focal mechanisms and hypocentre locations with a Gauss–Newton optimization that further refines the grid-search results. Using this algorithm, we have determined CMT solutions for 165 small to medium-sized earthquakes in Southern California. Preliminary comparison with the CMT solutions provided by the Southern California Seismic Network (SCSN) shows that our solutions generally provide better fit to the observed waveforms. When applied to a large number of earthquakes, our algorithm may provide a more robust CMT catalogue for earthquakes in Southern California.

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

[2]  G. McMechan,et al.  Source-Parameter Estimation by Full Waveform Inversion in 3D Heterogeneous, Viscoelastic, Anisotropic Media , 2001 .

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

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

[5]  Tiankai Tu,et al.  High Resolution Forward And Inverse Earthquake Modeling on Terascale Computers , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[6]  T. Okamoto.  Full waveform moment tensor inversion by reciprocal finite difference Green’s function , 2002 .

[7]  En-Jui Lee,et al.  Classification of Seismic Windows Using Artificial Neural Networks , 2011, ICCS.

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

[9]  Kim B. Olsen,et al.  Estimation of Q for Long-Period (>2 sec) Waves in the Los Angeles Basin , 2003 .

[10]  En-Jui Lee,et al.  Rapid Processing of Synthetic Seismograms Using Windows Azure Cloud , 2010, 2010 IEEE Second International Conference on Cloud Computing Technology and Science.

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

[12]  H. Kanamori,et al.  Seismic structure of the Transverse Ranges, California , 1977 .

[13]  Thomas H. Jordan,et al.  Full 3D Tomography for the Crustal Structure of the Los Angeles Region , 2007 .

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

[15]  En-Jui Lee,et al.  Perturbation kernels for generalized seismological data functionals (GSDF) , 2010 .

[16]  Andreas Fichtner,et al.  Full waveform tomography for radially anisotropic structure: New insights into present and past states of the Australasian upper mantle , 2009 .

[17]  Michael A. Saunders,et al.  Algorithm 583: LSQR: Sparse Linear Equations and Least Squares Problems , 1982, TOMS.

[18]  Gregory C. Beroza,et al.  Testing Community Velocity Models for Southern California Using the Ambient Seismic Field , 2008 .

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

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

[21]  Douglas S. Dreger,et al.  Source parameters of the Sierra Madre Earthquake from regional and local body waves , 1991 .

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

[23]  Location of shallow subduction-zone earthquakes inferred from teleseismic body waveforms , 1994 .

[24]  D. Dreger,et al.  Monitoring of strain release in central and northern California using broadband data , 1993 .

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

[26]  T. Okamoto.  Teleseismic synthetics obtained from 3‐D calculations in 2‐D media , 1994 .

[27]  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.

[28]  M. Sambridge,et al.  Seismic moment tensor inversion using a 3-D structural model: applications for the Australian region , 2011 .

[29]  P. Shearer,et al.  Applying a three-dimensional velocity model, waveform cross correlation, and cluster analysis to locate southern California seismicity from 1981 to 2005 , 2007 .

[30]  Donald V. Helmberger,et al.  Source Estimation from Broadband Regional Seismograms , 1994 .

[31]  Douglas S. Dreger,et al.  85.11 - TDMT_INV: Time Domain Seismic Moment Tensor INVersion , 2003 .

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

[33]  H. Thio,et al.  Moment-tensor inversions for local earthquakes using surface waves recorded at TERRAscope , 1995 .

[34]  Douglas S. Dreger,et al.  Determination of source parameters at regional distances with three‐component sparse network data , 1993 .

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

[36]  John Clinton,et al.  An Evaluation of the SCSN Moment Tensor Solutions: Robustness of the Mw Magnitude Scale, Style of Faulting, and Automation of the Method , 2006 .

[37]  William Hung Kan Lee,et al.  International handbook of earthquake and engineering seismology , 2002 .

[38]  D. Helmberger,et al.  Modeling the long-period body waves from shallow earthquakes at regional ranges , 1980 .

[39]  John H. Woodhouse,et al.  Determination of earthquake source parameters from waveform data for studies of global and regional seismicity , 1981 .

[40]  D. Komatitsch,et al.  Effects of crust and mantle heterogeneity on PP/P and SS/S amplitude ratios , 2002 .

[41]  K. Aki,et al.  Quantitative Seismology, 2nd Ed. , 2002 .

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

[43]  Thomas H. Jordan,et al.  Strain Green's Tensors, Reciprocity, and Their Applications to Seismic Source and Structure Studies , 2006 .

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

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

[46]  P. Shearer,et al.  A three-dimensional crustal seismic velocity model for southern California from a composite event method , 2007 .

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