Full 3D Tomography for the Crustal Structure of the Los Angeles Region

We apply full 3D tomography (f3dt) to a regional dataset in the Los Angeles area. In f3dt, the starting model as well as the model perturbation is three- dimensional and the sensitivity (Frechet) kernels are computed using numerical schemes that account for the full physics of 3D wave propagation. f3dt can account for the nonlinearity of structural inverse problem through iteration, thus providing the most efficient means for assimilating seismic observations into dynamic ground- motion models. We have successfully applied a scattering-integral (si) formulation of f3dt to improve a 3D elastic structure model, the Southern California Earthquake Center (scec) Community Velocity Model version 3.0 (CVM3.0) with lowest shear velocity clamped at 1 km/sec, in the Los Angeles region. Our data are time- and frequency-localized measurements of phase-delay anomalies relative to synthetics computed from the 3D elastic starting model. The Frechet kernels for our measurements were computed by convolving the earthquake wave fields generated by point earthquake sources with the receiver Green tensors (rgts), which are the spatial- temporal fields produced by three orthogonal unit impulsive point forces acting at the receiver locations. We inverted 7364 phase-delay measurements of P and S body waves by using the lsqr method. The revised 3D model, LAF3D, provides substantially better fit to the observed waveform data than the 3D starting model. To our knowledge, this study is the first successful application of f3dt using real data in structural seismology.

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

[2]  Paul G. Richards,et al.  Quantitative Seismology: Theory and Methods , 1980 .

[3]  Mostafa A. Bassiouni,et al.  Data Compression in Scientific and Statistical Databases , 1985, IEEE Transactions on Software Engineering.

[4]  G. Backus.  Long-Wave Elastic Anisotropy Produced by Horizontal Layering , 1962 .

[5]  Remo Guidieri.  Res , 1995, RES: Anthropology and Aesthetics.

[6]  F. Dahlen,et al.  Fréchet kernels for body-wave amplitudes , 2001 .

[7]  Stephen S. Gao,et al.  Localized amplification of seismic waves and correlation with damage due to the Northridge earthquake: Evidence for focusing in Santa Monica , 1996, Bulletin of the Seismological Society of America.

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

[9]  Gao,et al.  Northridge earthquake damage caused by geologic focusing of seismic waves , 2000, Science.

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

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

[12]  Jeroen Tromp,et al.  Finite-frequency kernels for wave propagation in porous media based upon adjoint methods , 2009 .

[13]  Kazuki Koketsu,et al.  Visualization of 3D Wave Propagation from the 2000 Tottori-ken Seibu, Japan, Earthquake: Observation and Numerical Simulation , 2003 .

[14]  Robert W. Graves,et al.  The seismic response of the Los Angeles basin, California , 1998, Bulletin of the Seismological Society of America.

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

[16]  Guust Nolet,et al.  Three-dimensional waveform sensitivity kernels , 1998 .

[17]  Robert L. Woodward,et al.  Global upper mantle structure from long-period differential travel times , 1991 .

[18]  J. Tromp,et al.  Finite-Frequency Kernels Based on Adjoint Methods , 2006 .

[19]  Robert W. Graves,et al.  Ground-motion amplification in the Santa Monica area: Effects of shallow basin-edge structure , 1998, Bulletin of the Seismological Society of America.

[20]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[21]  P. Maechling,et al.  Strong shaking in Los Angeles expected from southern San Andreas earthquake , 2006 .

[22]  Jeroen Tromp,et al.  Measurements and global models of surface wave propagation , 1997 .

[23]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[24]  Thomas H. Jordan,et al.  Sensitivity of frequency-dependent traveltimes to laterally heterogeneous, anisotropic Earth structure , 1998 .

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

[26]  Guust Nolet,et al.  Fréchet kernels for finite‐frequency traveltimes—II. Examples , 2000 .

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

[28]  Edward H. Field,et al.  A Modified Ground-Motion Attenuation Relationship for Southern California that Accounts for Detailed Site Classification and a Basin-Depth Effect , 2000 .

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

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

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

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

[33]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[34]  Gerard T. Schuster,et al.  Wave-equation traveltime inversion , 1991 .

[35]  T. Jordan,et al.  FAST TRACK PAPER: Full three-dimensional tomography: a comparison between the scattering-integral and adjoint-wavefield methods , 2007 .

[36]  T. Furumura,et al.  Visualization of 3 D Wave Propagation from the 2000 Tottori-ken Seibu , Japan , Earthquake : Observation and Numerical Simulation , 2003 .

[37]  Roel Snieder,et al.  Wavefield smoothing and the effect of rough velocity perturbations on arrival times and amplitudes , 1996 .

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

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

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

[41]  Paul M. Davis,et al.  Localized Amplification of Seismic Waves and Correlation with Damae due to the Northridge Earthquake: Evidence for Focusing in Santa Monica , 1996 .

[42]  Guust Nolet,et al.  Fréchet kernels for finite-frequency traveltimes—I. Theory , 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]  Joab R Winkler,et al.  Numerical recipes in C: The art of scientific computing, second edition , 1993 .

[47]  Richard A. Silverman,et al.  Wave Propagation in a Random Medium , 1960 .

[48]  O. Savasta,et al.  Wave-equation Traveltime Inversion Using 2D Cross-correlation Function , 1997 .

[49]  G. Backus,et al.  The Resolving Power of Gross Earth Data , 1968 .

[50]  Malcolm Sambridge,et al.  Non-linear arrival time inversion: constraining velocity anomalies by seeking smooth models in 3-D , 1990 .

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