Seismic ray path variations in a 3D global velocity model

Abstract A three-dimensional (3D) ray tracing technique is used to investigate ray path variations of P, PcP, pP and PP phases in a global tomographic model with P wave velocity changing in three dimensions and with lateral depth variations of the Moho, 410 and 660 km discontinuities. The results show that ray paths in the 3D velocity model deviate considerably from those in the average 1D model. For a PcP wave in Western Pacific to East Asia where the high-velocity (1–2%) Pacific slab is subducting beneath the Eurasian continent, the ray path change amounts to 27 km. For a PcP ray in South Pacific where very slow (−2%) velocity anomalies (the Pacific superplume) exist in the whole mantle, the maximum ray path deviation amounts to 77 km. Ray paths of other phases (P, pP, PP) are also displaced by tens of kilometers. Changes in travel time are as large as 3.9 s. These results suggest that although the maximal velocity anomalies of the global tomographic model are only 1–2%, rays passing through regions with strong lateral heterogeneity (in velocity and/or discontinuity topography) can have significant deviations from those in a 1D model because rays have very long trajectories in the global case. If the blocks or grid nodes adopted for inversion are relatively large (3–5°) and only a low-resolution 3D model is estimated, 1D ray tracing may be feasible. But if fine blocks or grid nodes are used to determine a high-resolution model, 3D ray tracing becomes necessary and important for the global tomography.

[1]  Harmen Bijwaard,et al.  Fast kinematic ray tracing of first‐ and later‐arriving global seismic phases , 1999 .

[2]  K. Condie.  Mantle Plumes and Their Record in Earth History: Superplume Events , 2001 .

[3]  Yosihiko Ogata,et al.  Whole mantle P-wave travel time tomography , 1990 .

[4]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

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

[6]  P. Shearer,et al.  Global mapping of topography on transition zone velocity discontinuities by stacking SS precursors , 1998 .

[7]  Gabi Laske,et al.  CRUST 5.1: A global crustal model at 5° × 5° , 1998 .

[8]  Yoshio Fukao,et al.  Application of a three-dimensional ray-tracing technique to global P, PP and Pdiff traveltime tomography , 2001 .

[9]  Hiroo Kanamori,et al.  Deep structure of Japan subduction zone as derived from local, regional, and teleseismic events , 1994 .

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

[11]  John A. Hildebrand,et al.  Depth Extent of the Lau Back-Arc Spreading Center and Its Relation to Subduction Processes , 1997 .

[12]  Harmen Bijwaard,et al.  Non‐linear global P‐wave tomography by iterated linearized inversion , 2000 .

[13]  Dapeng Zhao,et al.  Global tomographic images of mantle plumes and subducting slabs: insight into deep Earth dynamics , 2004 .

[14]  Akira Hasegawa,et al.  Tomographic imaging of P and S wave velocity structure beneath northeastern Japan , 1992 .

[15]  Dapeng Zhao,et al.  New advances of seismic tomography and its applications to subduction zones and earthquake fault zones: A review , 2001 .

[16]  Shigenori Maruyama,et al.  Mechanism of backarc opening in the Japan Sea: role of asthenospheric injection , 1990 .

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

[18]  S. Tanaka.  Very low shear wave velocity at the base of the mantle under the South Pacific Superswell , 2002 .

[19]  Dapeng Zhao,et al.  Seismic structure and origin of hotspots and mantle plumes , 2001 .

[20]  Wei-jia Su,et al.  Degree 12 model of shear velocity heterogeneity in the mantle , 1994 .

[21]  Clifford H. Thurber,et al.  A fast algorithm for two-point seismic ray tracing , 1987 .