Local earthquake (LE) tomography with joint inversion for P‐ and S‐wave velocities using structural constraints

In local earthquake (LE) tomography for P‐ and S‐wave velocities (or Vp/Vs ratios), constraints punishing deviations from a given Vp/Vs ratio are often used. In a synthetic model test we demonstrate that though such constraints inhibit unreasonable Vp/Vs ratio variations, they bias the resulting Vp/Vs ratios. As an alternative approach, structural constraints based on the cross‐product of the model gradients can be used to constrain the joint P‐ and S‐wave inversion. We show that the resulting models are as low in artifacts as if a Vp/Vs damping was used, but the resulting Vp/Vs ratios are less biased, which is important for a quantitative interpretation of physical properties and processes in the Earth.

[1]  G. Pavlis,et al.  The mixed discrete‐continuous inverse problem: Application to the simultaneous determination of earthquake hypocenters and velocity structure , 1980 .

[2]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[3]  H. M. Iyer,et al.  Vp/Vs ratios in the Yellowstone National Park region, Wyoming , 1985 .

[4]  Joan S. Gomberg,et al.  The effect of S-wave arrival times on the accuracy of hypocenter estimation , 1990, Bulletin of the Seismological Society of America.

[5]  D. Eberhart‐Phillips,et al.  Three‐dimensional P and S velocity structure in the Coalinga Region, California , 1990 .

[6]  Viale Romolo Gessi,et al.  Testing the reliability of Vp/ Vs anomalies in traveltime tomography , 1993 .

[7]  A. Michelini.  Testing the reliability of Vp/Vs anomalies in traveltime tomography , 1993 .

[8]  Jens Havskov,et al.  Earthquake prediction research in the South Iceland seismic zone and the SIL project , 1993 .

[9]  C. Sanders,et al.  Seismological evidence for magmatic and hydrothermal structure in Long Valley Caldera from local earthquake attenuation and velocity tomography , 1995 .

[10]  J. Lees,et al.  Poisson's ratio and porosity at Coso geothermal area, California , 2000 .

[11]  Ari Tryggvason,et al.  Three-dimensional imaging of the P- and S-wave velocity structure and earthquake locations beneath Southwest Iceland , 2002 .

[12]  C. Chiarabba,et al.  Tomographic images and 3D earthquake locations of the seismic swarm preceding the 2001 Mt. Etna eruption: Evidence for a dyke intrusion , 2002 .

[13]  M. Meju,et al.  Characterization of heterogeneous near‐surface materials by joint 2D inversion of dc resistivity and seismic data , 2003 .

[14]  M. Meju,et al.  Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints , 2004 .

[15]  A. Zollo,et al.  The deep structure of the Larderello‐Travale geothermal field from 3D microearthquake traveltime tomography , 2004 .