Seismic ray path variations in a 3D global velocity model
|
|
- Lilian Beverley Porter
- 6 years ago
- Views:
Transcription
1 Physics of the Earth and Planetary Interiors 141 (2004) Seismic ray path variations in a 3D global velocity model Dapeng Zhao, Jianshe Lei Geodynamics Research Center, Ehime University, Matsuyama , Japan Received 5 December 2002; accepted 22 November 2003 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 Elsevier B.V. All rights reserved. Keywords: Seismic tomography; Ray tracing; Seismic discontinuity; Subducting slabs; Mantle plumes 1. Introduction During the last two decades, seismic tomography has become the most powerful tool to explore the heterogeneous structure of the Earth. Many three-dimensional (3D) Earth models have been produced that have greatly advanced our knowledge about the structure and dynamics of the Earth s interior. Continuing efforts, however, are needed to improve the theoretical aspects of seismic tomography as well as its applications. Seismic ray tracing, Corresponding author. Tel.: ; fax: addresses: zhao@sci.ehime-u.ac.jp (D. Zhao), leijs@sci.ehime-u.ac.jp (J. Lei). the computation of ray paths and travel times in a seismic velocity model, is one of the key elements in traveltime tomographic imagings. The accuracy of ray tracing directly affects the final results of a tomographic inversion. Most of local and regional tomographic studies have already used 3D ray tracing schemes to compute travel times and ray paths (see a recent review by Zhao, 2001b). Because usually millions of travel time data are used, most of the global tomographic studies are still using one-dimensional (1D) ray tracing methods, that is, ray paths are only calculated for a 1D velocity model where velocity changes with depth alone, and the ray paths are fixed during the inversion process; ray path variations due to the lateral velocity heterogeneity are not considered /$ see front matter 2003 Elsevier B.V. All rights reserved. doi: /j.pepi
2 154 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) With the remarkable advances in computer technology, 3D ray tracing methods have become feasible in large tomographic computations. Recently, a few researchers have started to apply 3D ray tracing schemes in their global tomographic studies (e.g., Bijwaard and Spakman, 2000; Gorbatov et al., 2001; Zhao, 2001a). Bijwaard and Spakman (1999) examined ray path variations of first P waves in a 3D mantle velocity model. They found that for a direct P wave with an epicentral distance of 87, the ray path is displaced by nearly 100 km and the travel time change is larger than 2 s in their 3D velocity model. The changes in trajectory for other well used later phases, such as pp, PP, and PcP, in a 3D global velocity model are still unclear. Studying this subject is important for not only theoretical aspects of seismology but also practical imagings of the Earth s heterogeneous structure. The Moho depth ranges from about 10 km under oceans to km under the continents (Mooney et al., 1998). The 410 and 660 km discontinuities exhibit depth variations of up to 36 km on a global scale (Flanagan and Shearer, 1998). Such large depth changes of the discontinuities would certainly affect the travel times and ray paths of seismic waves. However, the changes in travel time and ray path due to the discontinuity topography have not been reported in literature, including Bijwaard and Spakman (1999). In this work, we have attempted to use a 3D ray tracing method to investigate ray path variations of P, pp, PP, and PcP phases due to the lateral velocity variations (Zhao, 2001a) as well as the topography of the Moho, 410 and 660 km discontinuities (Mooney et al., 1998; Flanagan and Shearer, 1998). 2. Three-dimensional ray tracing In this work, we have used a modified version (Zhao, 2001a) of the 3D ray tracing technique originally developed by Zhao et al. (1992). Its principle is shown schematically in Fig. 1. This ray tracing scheme can deal with a velocity model that contains several velocity discontinuities of complex geometry and with 3D velocity variations everywhere in the model. Consider a case as shown in Fig. 1a where a complexly shaped seismic velocity discontinuity, MM,exists and velocities on both sides of MM are continuous and change in three dimensions. A and B are two Fig. 1. (a) Schematic illustration of the scheme using Snell s law to find the location of intersection C between a seismic ray (ACB) and a discontinuity MM. (b) Schematic illustration of the 3D ray tracing method used in this study. The thick straight line denotes the initial ray path. The dashed, dotted and thin solid lines denote ray tracing processes using Snell s law, the pseudo-bending scheme, and Snell s law again, respectively. See text for details. points on different sides of MM. Velocities at A and B are V a and V b, respectively. C is the intersection between a seismic ray AB and the discontinuity MM. The velocities at C on the two sides of MM are V c1 and V c2, respectively. We take the arithmetic average V 1 of V a and V c1 and arithmetic average V 2 of V b and V c2 : V 1 = V a + V c1 2 V 2 = V b + V c2 2 (1) (2)
3 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) to represent the mean velocity around AC and BC, respectively. This approximation may be poor when A and B are far away from MM. When A and B are close to MM, however, the approximation will become sufficiently good. In order to find the location of point C, we first estimate the locations of two points, A and B, where AA and BB are perpendicular to MM. The point C should exist within the section of MM between A and B. Then we use the bisection method to gradually narrow the range in which C exists until we find the location of C which satisfies Snell s law: sin θ 1 = sin θ 2 (3) V 1 V 2 with sufficient accuracy. Fig. 1b shows the ray tracing algorithm schematically for the case with three discontinuities. The triangle denotes a station and the star a hypocenter. The straight line A1 A5 connecting the station and hypocenter is assumed to be the initial ray path. The points on the ray are divided into two types: one is the intersections between the ray and the discontinuities (A2, A3, A4), which we call discontinuous points (DPs); the other is the points in the continuous medium bounded by two adjacent discontinuities (B1, B2, B3, B4), which we call continuous points (CPs). The principle of the ray tracing algorithm is that Snell s law is used to perturb the DPs (Fig. 1a), while the pseudo-bending technique (Um and Thurber, 1987) is adopted to perturb the CPs. In Fig. 1b, we first use Snell s law to find new DPs A 2 from A1 and A3, A 3 from A2 and A4, and A 4 from A3 and A5. Then we use the pseudo-bending to find CPs B1 from A1 and A 2, B2 from A 2 and A 3, B3 from A 3 and A 4, and B4 from A 4 and A5. Then again using Snell s law we find A 2 from B1 and B2, A 3 from B2 and B3, and A 4 from B3 and B4, and so on. After a number of iterations the ray path converges to its true location. It is found that this algorithm is very efficient and works very well for 3D seismic ray tracing in local to regional scales (epicentral distance less than about 1000 km) (Zhao et al., 1992, 1994). In Zhao (2001a) and the present study, we found that the original routine of this algorithm developed by Zhao et al. (1992) gives less accurate results for mantle rays in a global scale (P, PP, PcP and pp, etc.); the epicentral distance is as large as 100 (1 = km). After careful examinations we found that the problem was caused by the enhancement factor used in the pseudo-bending technique (Um and Thurber, 1987). The enhancement factor is adopted to quicken the convergence of the ray tracing process. We found that the optimal value of the enhancement factor is dependent on the length of seismic rays. The optimal enhancement factor for the local and regional scale is inadequate for seismic rays in the global scale. After numerous tests we determined the optimal values of the enhancement factor for different lengths of the mantle seismic rays with different focal depths. Our modified version of the 3D ray tracing algorithm gives very accurate results with the absolute computational error in travel time smaller than 0.05 s, which is considered to be accurate enough because the uncertainty of arrival time pickings is generally 0.1 s or greater for the mantle rays in the global scale. Our new 3D ray tracing code can deal with not only the first P and S rays but also the later phases such as pp, PP, PcP and Pdiff, etc. 3. 3D global velocity model We examined seismic rays along many cross sections using the 3D ray tracing technique and the global P wave tomography model of Zhao (2001a). Here we show examples of the 3D ray tracing for three regions: Western Pacific to East Asia, South Pacific, and Tibet (Fig. 2). The subducting Pacific slab is clearly imaged in the upper mantle beneath the Japan Island and the Sea of Japan (Fig. 3a). Intermediate-depth and deep earthquakes form a clear Wadati Benioff zone within the high-velocity (high-v) subducting slab. The slab becomes stagnant in the mantle transition zone under Northeast China. Moderately fast anomalies in the lower mantle represent pieces of the old Pacific slab which are collapsing down to core mantle boundary (CMB) after stagnating in the transition zone for a long time (Maruyama, 1994; Zhao, 2004). Prominent low-velocity (low-v) anomalies are visible under the active Changbai volcano and the Datong Quaternary volcano, which may be caused by the deep dehydration process of the subducting Pacific slab and the convective circulation process in the mantle wedge (Tatsumi et al., 1990; Zhao et al., 1994, 1997; Zhao,
4 156 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 2. Map showing locations of three vertical cross sections (AB, CD and EF) in which seismic ray paths are investigated in this work. Stars and open triangles denote hypothesized earthquake hypocenters and seismic stations, respectively. Solid triangles denote the surface hotspot volcanoes except for those in East China (WV, Wudalianchi volcano; DV, Datong volcano; CV, Changbai volcano). TA, Tahiti hotspot. 2004). Hence the active intraplate volcanism in Northeast China is not hotspot like Hawaii but a kind of back-arc volcanism. Slow anomalies are also visible under the Pacific slab in the transition zone and the lower mantle, which represent upwellings either related to a mantle plume or caused by a local mantle convection associated with the subduction of the Pacific slab (Zhao, 2004). Several hotspot volcanoes are located on the South Pacific Superswell (McNutt, 1998) (Fig. 2). A huge slow anomaly with a lateral extent of over 1000 km is visible in the entire mantle from the surface down to the CMB (Fig. 3b), which represents the Pacific superplume as has been detected by many previous studies with various geophysical and geochemical approaches (Su et al., 1994; Maruyama, 1994; McNutt, 1998; Condie, 2001; Zhao, 2001a; Tanaka, 2002). Detailed resolution analyses confirmed that the major structural features in Fig. 3 as mentioned above are reliable (Zhao, 2004). The geodynamic implications of the tomographic results (Fig. 3) are discussed in detail by Zhao (2004). 4. Influence of lateral velocity variations Fig. 3 shows the 1D and 3D ray paths of P, pp, PP and PcP phases along the AB and CD profiles. A focal depth of 300 km is adopted, though any other depth values are feasible. Here 3D ray path means the ray trajectory that is determined by applying the 3D ray tracing method (Fig. 1) to a velocity model with 3D P wave velocity variations determined by Zhao (2001a), while all the velocity discontinuities (the Moho at 35, 410 and 660 km) are spherical interfaces without lateral depth variations. The 1D ray path denotes that determined by using the ray tracing scheme to a 1D velocity model that is the average of the 3D model for each depth. The average 1D velocity model is nearly the same as the iasp91 Earth model (Kennett and Engdahl, 1991), the difference between them is less than 1%. To see the ray path variations more clearly, we plot the differences between the 1D and 3D rays in the horizontal and vertical directions as well as the total displacement in Figs The straight lines show the 1D rays, while the curved lines show the 3D rays. P wave velocity images along the 1D ray paths in the horizontal and vertical directions are also shown (Figs. 4 7). The changes in travel time and ray path for the rays along profile AB are listed in Table 1. The results for those in profile CD are listed in Table 2. In addition to the travel time (T11) spent by a 1D ray path in the 1D velocity model and that (T33) taken by a 3D ray path in the 3D velocity model, we also computed the travel times (T13) of 1D rays with velocities of the 3D
5 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 3. 1D and 3D ray paths of P, pp, PP and PcP phases and P wave velocity images (Zhao, 2001a) along the profile AB (a, c) and CD (b) in Fig. 2. Blue and red colors denote fast and slow velocities, respectively. The dotted and solid lines denote 1D and 3D rays, respectively. The two dashed lines show the 410 and 660 km discontinuities. The Moho is taken as a spherical surface at 35 km depth. Solid triangles denote active volcanoes, small white dots in (a) and (c) show earthquakes that occurred within a 150 km width of profile AB. Stars and open triangles denote hypothesized earthquake hypocenters and seismic stations, respectively. JA, Japan arc; JT, Japan trench. See text for details. model and that (T31) of 3D rays with velocities of the 1D model. In most of the previous global tomographic inversions (e.g., Inoue et al., 1990), ray paths are fixed to those in a 1D velocity model, while travel times are updated by using the 3D velocity distribution obtained in each of the iterations during the inversion process. The differences between the four travel times (T11, T33, T13 and T31) for the four types of phases (P, pp, PP, and PcP) are shown in Tables 1 and 2. For the rays along profile AB (Fig. 3a), the travel time differences are relatively small (<0.2 s), and the ray path deviations are generally less than 33 km. Table 1 Travel time and ray path changes for rays along profile AB Phase T11 T13 (s) T11 T31 (s) T13 T33 (s) T11 T33 (s) DMH (km) DMV (km) DMT (km) DS3D-1D (km) P PcP pp PP T ij denotes travel time of an id ray path in a jd velocity model. DMH, DMV and DMT denote the maximum horizontal, vertical and total ray path deviations, respectively. DS3D-1D denotes the difference in ray length between 3D and 1D ray paths. See text for details.
6 158 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 4. Difference in ray path (in km) between the 1D and 3D rays vs. epicentral distance (in degree) for direct P phase (a c) and pp phase (d f) along the profile AB as shown in Fig. 3a. The straight and curved lines denote the 1D and 3D rays, respectively. (a) Horizontal, (b) vertical and (c) total ray deviations for direct P wave. (d) Horizontal, (e) vertical and (f) total ray deviations for pp wave. P wave velocity images along the 1D ray paths are shown in gray scale. The velocity perturbation scale is shown at the bottom. Table 2 Travel time and ray path changes for rays along profile CD Phase T11 T13 (s) T11 T31 (s) T13 T33 (s) T11 T33 (s) DMH (km) DMV (km) DMT (km) DS3D-1D (km) P PcP pp PP T ij denotes travel time of an id ray path in a jd velocity model. DMH, DMV and DMT denote the maximum horizontal, vertical and total ray path deviations, respectively. DS3D-1D denotes the difference in ray length between 3D and 1D ray paths. See text for details.
7 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 5. The same as Fig. 4 but for PcP (a c) and PP rays (d f). These relatively small differences in both travel time and ray path are related to the alternate low- and high-velocity anomalies along the rays (see Figs. 3a, 4 and 5). For the rays along profile CD, however, the travel time differences become much larger, amounting to 3.1 s, because of the huge slow anomalies associated with the Pacific superplume. The differential times (T13 T33) amount to 0.7 s, exceeding the picking accuracy of arrival time data recorded by short-period instruments ( s), as is the case for the International Seismological Center (ISC) data sets which were used by many global tomographic studies including Zhao (2001a). These results suggest that 1D ray tracing could cause large calculation errors for rays which pass through strong velocity anomalies. We can also see in Tables 1 and 2 that T33 is smaller than T13 for all the cases, indicating that 3D ray tracing is always better than 1D ray tracing. Ray trajectories are sensitive to lateral and vertical velocity gradients, and they are apt to move toward
8 160 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 6. The same as Fig. 4 but for direct P (a c) and pp rays (d f) along the profile CD in Fig. 3b. faster areas in a 3D velocity field. This is clearly visible in Fig. 3. Because of the slow anomaly under the Changbai volcano, the PP bounce point at the surface is moved by 16 km toward the west (Fig. 3a). The easternmost segment of PP is elevated by 27 km to escape the slow anomaly under the subducting Pacific slab (Figs. 3a and 5e). Because of the high-v stagnant slab in the mantle transition zone, the western portion of the long segment of pp is elevated by 12 km (Figs. 3a and 4e). The PcP bounce point at the CMB is moved by 18 km westward from a low-v zone to a high-v patch in the lowermost mantle (Fig. 3a). The most significant displacement in ray path occurs for the northern segment of PcP under South Pacific (Figs. 3b and 7c). The horizontal, vertical and total deviations in ray path are 26, 48, and 77 km, respectively (Fig. 7a c). The ray is displaced from the large slow Pacific superplume toward fast areas in the north (Figs. 3b and 7c). The northern segment of pp is moved by 32 km toward a high-v zone in the north
9 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 7. The same as Fig. 6 but for PcP (a c) and PP rays (d f). (Figs. 3b and 6e). In addition, the PP bounce point at the surface is displaced by 20 km toward the north to escape the slow anomaly under Tahiti. Most of the PP and pp rays propagate a long distance in and around the mantle transition zone and their paths show complex changes of significant amounts because of the existence of high-v stagnant slabs and low-v mantle plumes in the region. Thus, calculating travel times and ray paths of PP and pp phases precisely using a 3D ray tracing method becomes particularly important and necessary. Fig. 3c shows five 1D and 3D rays for direct P waves passing through the upper mantle and transition zone under the Western Pacific region; their epicentral distances are 14,17,20,25 and 29, respectively. The hypothesized hypocenter is located near the Japan trench with a focal depth of 75 km. The high-v subducting slab and the low-v anomalies in the mantle
10 162 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 8. Ray paths of P, pp and PP phases along the profiles AB (a) and CD (b) in Fig. 2. The dotted lines denote rays in the 1D iasp91 model. The solid lines show the rays in a heterogeneous model with lateral depth variations of the Moho discontinuity. The 410 and 660 km discontinuities are assumed to be spherical interfaces without lateral depth variations. The dashed lines denote the 410 and 660 km discontinuities and the global-average Moho discontinuity at 35 km depth. The Moho depth variations (Mooney et al., 1998) along the profiles AB and CD are shown in thin solid lines. The iasp91 Earth model is used for the 1D velocities in each layer. Lateral variations in seismic velocity are not considered. wedge and below the slab cause significant changes in the travel times and ray paths. For the rays P1 to P5 (Fig. 3c), the travel time difference (T11 T33) is 0.64, 0.11, 1.58, 2.24 and 1.11 s, respectively, and the maximum displacement in ray path is 6, 19, 10, 9 and 8 km, respectively. The rays bend toward the high-v slab. 5. Influence of the discontinuity topography The changes in seismic ray paths caused by the depth variations of the Moho, 410 and 660 km discontinuities along the three profiles in Fig. 2 are shown in Figs The models by Mooney et al. (1998) and Flanagan and Shearer (1998) are used for the Fig. 9. The same as Fig. 8 but for three direct P phases along the profile EF in Fig. 2. The focal depth is 80 km.
11 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 10. The same as Fig. 8 but for rays (solid lines) in a velocity model with a flat Moho at 35 km depth and with lateral depth variations of the 410 and 660 km discontinuities from the model of Flanagan and Shearer (1998). Fig. 11. The same as Fig. 8 but for rays (solid lines) in a velocity model with lateral depth variations of the Moho, 410 and 660 km discontinuities from the models of Mooney et al. (1998) and Flanagan and Shearer (1998).
12 164 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) topography of the three discontinuities. The iasp91 Earth model (Kennett and Engdahl, 1991) is used for the 1D P wave velocities in the crust and mantle; lateral velocity variations are not considered. In Figs. 8 11, ray path changes in the models with depth variations of the discontinuities are compared with those in the 1D iasp91 Earth model where the Moho (35 km depth), 410 and 660 km discontinuities are taken as spherical interfaces without lateral depth variations. In Figs. 8 and 9, ray path changes caused only by the Moho topography are examined. The Moho depth changes from 12 to 38 km along the profile AB, and from 11 to 14 km along the profile CD. The Moho topography caused the surface bounce point of the PP ray displaced by 34 km in Fig. 8a and 40 km in Fig. 8b. The maximum change in the pp ray path is 17 km along profile AB and 11 km along profile CD. The P rays are little affected by the Moho topography, though their travel times are changed significantly. In Fig. 8a, the travel time changes for P, pp and PP rays are 0.91, 0.91, and 1.10 s, respectively. In Fig. 8b, the corresponding travel time changes are 0.91, 3.02, and 3.87 s, respectively. Along the profile EF passing through the Himalaya mountains (Fig. 2), the Moho depth changes from 39 to 60 km (Fig. 9). For rays P1, P2 and P3 in Fig. 9, the travel time changes are 1.43, 2.13 and 1.16 s, and the maximum deviations in ray paths are 10, 17 and 9 km, respectively. Ray path changes caused by the topography of the 410 and 660 discontinuities are shown in Fig. 10. The Moho is flat with a depth of 35 km. The ray path of direct P wave is only slightly affected by the 410 and 660 km topography; the path deviations are less than 6 km in both AB and CD profiles. The maximum path deviations of pp and PP rays are 9 and 34 km in Fig. 10a, and 13 and 11 km in Fig. 10b, respectively. The travel time changes of P, pp and PP rays are 0.27, 0.42 and 1.00 s in Fig. 10a, and 0.12, 0.10 and 0.44 s in Fig. 10b, respectively. Fig. 12. The same as Fig. 3 but for rays (solid lines) in a velocity model with lateral depth variations of the Moho, 410 and 660 km discontinuities from the models of Mooney et al. (1998) and Flanagan and Shearer (1998).
13 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Fig. 11 shows the ray path changes caused by the lateral depth variations of all the three velocity discontinuities (the Moho, 410 and 660). The maximum path changes of P, pp and PP rays are 4, 19 and 28 km in profile AB, and 6, 6 and 30 km in profile CD, respectively. The corresponding changes in travel times are 0.64, 0.52 and 0.29 s in Fig. 11a, and 0.80, 2.95 and 3.42 s in Fig. 11b, respectively. 6. Composite effect of velocity changes and discontinuity topography Fig. 12 shows P, pp and PP rays in the model of Zhao (2001a) where both P wave velocity changes in three dimensions and the Moho, 410 and 660 km discontinuities have lateral depth variations as revealed by Mooney et al. (1998) and Flanagan and Shearer (1998). These 3D ray paths (solid lines in Fig. 12) are compared with those in the 1D iasp91 model (dotted lines in Fig. 12). The maximum path differences between the 1D and 3D rays for P, pp and PP phases are 12, 16 and 26 km in Fig. 12a, and 46, 38 and 32 km in Fig. 12b, respectively. The corresponding travel time changes are 0.59, 0.73 and 0.14 s in Fig. 12a, and 0.95, 1.62 and 0.45 s in Fig. 12b, respectively. 7. Discussion and conclusions Tracing seismic rays in a heterogeneous Earth model precisely and quickly is essential in tomographic studies because usually a great number of arrival times are used and, for each arrival, ray tracing is employed several times to conduct hypocenter relocation and structure determination. Although several efficient algorithms for 3D ray tracing have been developed (see Zhao, 2001b) and there have been great advances in computer technology, continuing efforts are still needed to search better ray tracing schemes in order to use more data to image the Earth structure in greater details and in higher accuracy. Our present results suggest that although the maximal velocity anomalies of the 3D global model are only 1 2%, rays passing through those very heterogeneous regions (such as subducting slabs and mantle plumes) can have deviations from the 1D rays amounting to 77 km because of the very long ray trajectories in the global case. Since damping and smoothing are usually applied in the tomographic inversions, the amplitude of velocity anomalies in the tomographic models is generally underestimated. Actual amplitudes of mantle heterogeneities would be greater. For example, high-resolution (5 30 km) tomographic imagings of local and regional scales produce velocity anomalies of up to 6% for P wave and 10% for S-wave (Zhao et al., 1992, 1997). Unlike the global case, in local or regional scale studies, arrival time data have high quality and ray path coverage is good, thus heavy damping and smoothing as in the global tomography can be avoided, and so not only the pattern but also the amplitude of velocity heterogeneity can be recovered well. With the improvement of spatial resolution of global tomography, greater amplitudes of velocity anomalies are expected. Thus ray paths in the mantle would have larger deviations. The significant lateral depth variations of the Moho, 410 and 660 km discontinuities (e.g., Mooney et al., 1998; Flanagan and Shearer, 1998) make the changes in ray path and travel time more complex. All these necessitate the use of 3D ray tracing in the future generation of global tomographic inversions. Recently, the computation of frequency-dependent travel times and Frechet kernels in a heterogeneous Earth model has been paid much attention (e.g., Dahlen et al., 2000; Hung et al., 2000; Zhao et al., 2000). Although it will be a challenge to apply the new approach to conduct tomographic inversions with a great number of travel time data, it may enable us to better understand the propagational behaviors of seismic waves in the heterogeneous interior of the Earth. Acknowledgements This work was partially supported by a grant to D. Zhao (Kiban-B No ) and a grant to E. Takahashi (Tokyo Institute of Technology) and D. Zhao (Special grant No ) from Ministry of Education and Science, Japan. K. Creager and two anonymous referees provided thoughtful comments which improved the manuscript. References Bijwaard, H., Spakman, W., Fast kinematic ray tracing of first- and later-arriving global seismic phases. Geophys. J. Int. 139,
14 166 D. Zhao, J. Lei / Physics of the Earth and Planetary Interiors 141 (2004) Bijwaard, H., Spakman, W., Non-linear global P-wave tomography by iterated linearized inversion. Geophys. J. Int. 141, Condie, K., Mantle Plumes and their Record in Earth History. Cambridge University Press, Cambridge, UK, p Dahlen, F., Hung, S., Nolet, G., Frechet kernels for finite-frequency traveltimes. I. Theory. Geophys. J. Int. 141, Flanagan, M., Shearer, P., Global mapping of topography on transition zone velocity discontinuities by stacking SS precursors. J. Geophys. Res. 103, Gorbatov, A., Fukao, Y., Widiyantoro, S., Application of a three-dimensional ray-tracing technique to global P, PP and Pdiff traveltime tomography. Geophys. J. Int. 146, Hung, S., Dahlen, F., Nolet, G., Frechet kernels for finite-frequency traveltimes. I. Examples. Geophys. J. Int. 141, Inoue, H., Fukao, Y., Tanabe, K., Ogata, Y., Whole mantle P wave travel time tomography. Phys. Earth Planet. Inter. 59, Kennett, B., Engdahl, E., Traveltimes for global earthquake location and phase identification. Geophys. J. Int. 105, Maruyama, S., Plume tectonics. J. Geol. Soc. Jpn. 100, McNutt, M., Superswells. Rev. Geophys. 36, Mooney, W., Laske, G., Master, T., CRUST 5.1: a global crustal model at 5 5. J. Geophys. Res. 103, Su, W., Woodward, R., Dziewonski, A., Degree 12 model of shear velocity heterogeneity in the mantle. J. Geophys. Res. 99, Tanaka, S., Very low shear wave velocity at the base of the mantle under the South Pacific Superswell. Earth Planet. Sci. Lett. 203, Tatsumi, Y., Maruyama, S., Nohda, S., Mechanism of backarc opening in the Japan Sea: role of asthenospheric injection. Tectonophysics 181, Um, J., Thurber, C., A fast algorithm for two-point seismic ray tracing. Bull. Seism. Soc. Am. 77, Zhao, D., 2001a. Seismic structure and origin of hotspots and mantle plumes. Earth Planet. Sci. Lett. 192, Zhao, D., 2001b. New advances of seismic tomography and its applications to subduction zones and earthquake fault zones: a review. The Island Arc 10, Zhao, D., Global tomographic images of mantle plumes and subducting slabs: insight into deep Earth dynamics. Phys. Earth Planet. Inter., in press. Zhao, D., Hasegawa, A., Horiuchi, S., Tomographic imaging of P and S wave velocity structure beneath northeastern Japan. J. Geophys. Res. 97, Zhao, D., Hasegawa, A., Kanamori, H., Deep structure of Japan subduction zone as derived from local, regional and teleseismic events. J. Geophys. Res. 99, Zhao, D., Xu, Y., Wiens, D., Dorman, L., Hildebrand, J., Webb, S., Depth extent of the Lau back-arc spreading center and its relation to subduction processes. Science 278, Zhao, L., Jordan, T., Chapman, C., Three-dimensional Frechet differential kernels for seismic delay times. Geophys. J. Int. 141,
Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima
1 2 3 Auxiliary materials for Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima nuclear power plant area 4 5 6 7 8 9 Ping Tong 1,2, Dapeng Zhao 1 and Dinghui Yang 2 [1] {Department of Geophysics,
More informationSmall scale hot upwelling near the North Yellow Sea of eastern China
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L20305, doi:10.1029/2008gl035269, 2008 Small scale hot upwelling near the North Yellow Sea of eastern China Yinshuang Ai, 1 Tianyu Zheng, 1 Weiwei Xu, 1 and Qiang
More informationTomographic imaging of P wave velocity structure beneath the region around Beijing
403 Doi: 10.1007/s11589-009-0403-9 Tomographic imaging of P wave velocity structure beneath the region around Beijing Zhifeng Ding Xiaofeng Zhou Yan Wu Guiyin Li and Hong Zhang Institute of Geophysics,
More information27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
P AND S WAVE VELOCITY STRUCTURE OF THE CRUST AND UPPER MANTLE UNDER CHINA AND SURROUNDING AREAS FROM BODY AND SURFACE WAVE TOMOGRAPHY M. Nafi Toksöz, Robert D. Van der Hilst, Youshun Sun, and Chang Li
More informationEstimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
More informationC3.4.1 Vertical (radial) variations in mantle structure
C3.4 Mantle structure Mantle behaves as a solid on short time scales (seismic waves travel through it and this requires elastic behaviour). Over geological time scales the mantle behaves as a very viscous
More informationChangbaishan volcanism in northeast China linked to subduction-induced mantle upwelling
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2166 Changbaishan volcanism in northeast China linked to subduction-induced mantle upwelling Youcai Tang, Masayuki Obayashi, Fenglin Niu, Stephen P. Grand, Yongshun
More informationFenglin Niu and Hitoshi Kawakatsu. Earthquake Research Institute, University of Tokyo, Yayoi, Bunkyo-ku, Tokyo , Japan
Earth Planets Space, 50, 965 975, 1998 Determination of the absolute depths of the mantle transition zone discontinuities beneath China: Effect of stagnant slabs on transition zone discontinuities Fenglin
More information3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES
3D IMAGING OF THE EARTH S MANTLE: FROM SLABS TO PLUMES Barbara Romanowicz Department of Earth and Planetary Science, U. C. Berkeley Dr. Barbara Romanowicz, UC Berkeley (KITP Colloquium 9/11/02) 1 Cartoon
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationSeismogram Interpretation. Seismogram Interpretation
Travel times in the Earth Ray paths, phases and their name Wavefields in the Earth: SH waves, P-SV waves Seismic Tomography Receiver Functions Seismogram Example Long-period transverse displacement for
More informationSome aspects of seismic tomography
Some aspects of seismic tomography Peter Shearer IGPP/SIO/U.C. San Diego September 7, 2009 Earthquake Research Institute Part 1: Global Tomography P velocity perturbations 200 km 1000 km 2700 km MIT 2006
More informationThe Earth s Structure from Travel Times
from Travel Times Spherically symmetric structure: PREM - Crustal Structure - Upper Mantle structure Phase transitions Anisotropy - Lower Mantle Structure D D - Structure of of the Outer and Inner Core
More informationMantle dynamics of Western Pacific and East Asia: Insight from seismic tomography and mineral physics
Gondwana Research 11 (2007) 120 131 www.elsevier.com/locate/gr Mantle dynamics of Western Pacific and East Asia: Insight from seismic tomography and mineral physics Dapeng Zhao a, Shigenori Maruyama b,
More informationTracing rays through the Earth
Tracing rays through the Earth Ray parameter p: source receiv er i 1 V 1 sin i 1 = sin i 2 = = sin i n = const. = p V 1 V 2 V n p is constant for a given ray i 2 i 3 i 4 V 2 V 3 V 4 i critical If V increases
More informationDeep structure and origin of the Baikal rift zone
Earth and Planetary Science Letters 243 (2006) 681 691 www.elsevier.com/locate/epsl Deep structure and origin of the Baikal rift zone Dapeng Zhao a,, Jianshe Lei a, Tomofumi Inoue a, Akira Yamada a, Stephen
More informationCrust and upper mantle P- and S-wave delay times at Eurasian seismic stations
Physics of the Earth and Planetary Interiors 123 (2001) 205 219 Crust and upper mantle P- and S-wave delay times at Eurasian seismic stations E.R. Engdahl, M.H. Ritzwoller Center for Imaging the Earth
More informationANEWJOINTP AND S VELOCITY MODEL OF THE MANTLE PARAMETERIZED IN CUBIC B-SPLINES
ANEWJOINTP AND S VELOCITY MODEL OF THE MANTLE PARAMETERIZED IN CUBIC B-SPLINES Michael Antolik, Göran Ekström, Adam M. Dziewonski, Yu J. Gu, Jian-feng Pan, and Lapo Boschi Department of Earth and Planetary
More informationContinent-sized anomalous zones with low seismic velocity at the base of Earth s mantle
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2733 Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle Edward J. Garnero 1, Allen K. McNamara 1, and Sang-Heon D. Shim 1
More informationCOMPUTATION OF REGIONAL TRAVEL TIMES AND STATION CORRECTIONS FROM THREE-DIMENSIONAL VELOCITY MODELS
COMPUTATION OF REGIONAL TRAVEL TIMES AND STATION CORRECTIONS FROM THREE-DIMENSIONAL VELOCITY MODELS A. Villaseñor 1,2, M.P. Barmin 1, M.H. Ritzwoller 1, and A.L. Levshin 1 1 Department of Physics, University
More information27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
CRUSTAL AND MANTLE STRUCTURE BENEATH EASTERN EURASIA FROM FINITE FREQUENCY SEISMIC TOMOGRAPHY (FFST) Ting Yang 2, Yang Shen 2, and Xiaoping Yang 1 Science Applications International Corporation 1 and University
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 188, 293 300 doi: 10.1111/j.1365-246X.2011.05256.x An analysis of SS precursors using spectral-element method seismograms L. Bai, Y. Zhang and
More informationMorphology of the distorted subducted Pacific slab beneath the Hokkaido corner, Japan
Physics of the Earth and Planetary Interiors 156 (2006) 1 11 Morphology of the distorted subducted Pacific slab beneath the Hokkaido corner, Japan M.S. Miller a,, B.L.N. Kennett a,1, A. Gorbatov b,2 a
More informationBasic Ray Tracing. Rick Aster and Sue Bilek. October 3, 2003
Basic Ray Tracing Rick Aster and Sue Bilek October 3, 3 A key observation that we can make about a seismic signal is its arrival time. From systematic observations of arrival times, we can deduce useful
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Seismic evidence for a global low velocity layer within the Earth s upper mantle SUPPLEMENTARY MATERIAL Benoît Tauzin 1, Eric Debayle 2 & Gérard Wittlinger 3 1 Department of Earth
More informationContents of this file
Geophysical Research Letters Supporting Information for Intraplate volcanism controlled by back-arc and continental structures in NE Asia inferred from trans-dimensional ambient noise tomography Seongryong
More informationAPPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY
APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY Timur TEZEL Supervisor: Takuo SHIBUTANI MEE07169 ABSTRACT In this study I tried to determine the shear wave velocity structure in the crust
More informationImaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges
2017 SCEC Proposal Report #17133 Imaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges Principal Investigator Zhongwen Zhan Seismological Laboratory, California
More informationSeismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone
IJMS 2017 vol. 4 (2): 49-54 International Journal of Multidisciplinary Studies (IJMS) Volume 4, Issue 2, 2017 DOI: http://doi.org/10.4038/ijms.v4i2.22 Seismic Activity near the Sunda and Andaman Trenches
More informationGEOPHYSICAL RESEARCH LETTERS, VOL. 35, L14308, doi: /2008gl034461, 2008
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35,, doi:10.1029/2008gl034461, 2008 Tomographic evidence for hydrated oceanic crust of the Pacific slab beneath northeastern Japan: Implications
More informationSeismic Velocity Structure in the Crust and Upper Mantle beneath Northern Japan
J. Phys. Earth, 42, 269-301, 1994 Seismic Velocity Structure in the Crust and Upper Mantle beneath Northern Japan Hiroki Miyamachi,l,* Minoru Kasahara,2 Sadaomi Suzuki,2,** Kazuo Tanaka,3 and Akira Hasegawa
More informationMagnitude 8.3 SEA OF OKHOTSK
A powerful earthquake in Russia's Far East was felt as far away as Moscow, about 7,000 kilometers (4,400 miles) west of the epicenter, but no casualties or damage were reported. The epicenter was in the
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
TRANSITION ZONE WAVE PROPAGATION: CHARACTERIZING TRAVEL-TIME AND AMPLITUDE INFORMATION Peter M. Shearer and Jesse F. Lawrence University of California San Diego, Institute of Geophysics and Planetary Physics
More informationUSU 1360 TECTONICS / PROCESSES
USU 1360 TECTONICS / PROCESSES Observe the world map and each enlargement Pacific Northwest Tibet South America Japan 03.00.a1 South Atlantic Arabian Peninsula Observe features near the Pacific Northwest
More informationANOTHER MEXICAN EARTHQUAKE! Magnitude 7.1, Tuesday Sept. 19, 2017
ANOTHER MEXICAN EARTHQUAKE! Magnitude 7.1, Tuesday Sept. 19, 2017 Why is there no oceanic crust older than 200 million years? SUBDUCTION If new oceanic crust is being continuously created along the earth
More informationWhole Mantle Convection
Whole Mantle Convection Overview 1. Evidence for whole mantle convection 2. Model of whole mantle convection reconciling geophysical and geochemical data Transition Zone Water Filter Model 3. Evidence
More informationSupplementary Online Material for. Seismic evidence for a chemically distinct thermochemical reservoir in Earth s deep mantle beneath Hawaii
Supplementary Online Material for Seismic evidence for a chemically distinct thermochemical reservoir in Earth s deep mantle beneath Hawaii Authors: Chunpeng Zhao 1, Edward J. Garnero 1,*, Allen K. McNamara
More informationTrans-Pacific temperature field in the mantle transition region from seismic and electromagnetic tomography
Trans-Pacific temperature field in the mantle transition region from seismic and electromagnetic tomography Yoshio Fukao 1,3, Takao Koyama 2, Masayuki Obayashi 1 and Hisashi Utada 3 1 Research Program
More informationBroadband converted phases from midmantle discontinuities
Earth Planets Space, 50, 987 997, 1998 Broadband converted phases from midmantle discontinuities Lev Vinnik 1,2, Fenglin Niu 1, and Hitoshi Kawakatsu 1 1 Earthquake Research Institute, University of Tokyo,
More informationGEOPHYSICAL RESEARCH LETTERS, VOL. 31, L19604, doi: /2004gl020366, 2004
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L19604, doi:10.1029/2004gl020366, 2004 Characteristic seismic activity in the subducting plate boundary zone off Kamaishi, northeastern Japan, revealed by precise
More informationThe influence of short wavelength variations in viscosity on subduction dynamics
1 Introduction Deformation within the earth, driven by mantle convection due primarily to cooling and subduction of oceanic lithosphere, is expressed at every length scale in various geophysical observations.
More informationSUPPLEMENTARY INFORMATION
The major uncertainties in our model predictions arise from the input parameters, which include mantle density models (i.e. seismic tomography and choices about scaling velocities to temperature), crustal
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/326/5949/112/dc1 Supporting Online Material for Global Surface Wave Tomography Using Seismic Hum Kiwamu Nishida,* Jean-Paul Montagner, Hitoshi Kawakatsu *To whom correspondence
More informationSupplementary Figure 1. Distribution of seismic event locations determined using the final 3-D velocity model. We separate the crust-related
Supplementary Figure 1. Distribution of seismic event locations determined using the final 3-D velocity model. We separate the crust-related seismicity at depths of less than 40 km (panel A) from the deeper
More informationScattering objects in the lower mantle beneath northeastern China observed with a short-period seismic array
Physics of the Earth and Planetary Interiors 138 (2003) 55 69 Scattering objects in the lower mantle beneath northeastern China observed with a short-period seismic array Tadashi Kito a,, Takuo Shibutani
More informationGlobal Tectonics. Kearey, Philip. Table of Contents ISBN-13: Historical perspective. 2. The interior of the Earth.
Global Tectonics Kearey, Philip ISBN-13: 9781405107778 Table of Contents Preface. Acknowledgments. 1. Historical perspective. 1.1 Continental drift. 1.2 Sea floor spreading and the birth of plate tectonics.
More informationTime dependence of PKP(BC) PKP(DF) times: could this be an artifact of systematic earthquake mislocations?
Physics of the Earth and Planetary Interiors 122 (2000) 221 228 Time dependence of PKP(BC) PKP(DF) times: could this be an artifact of systematic earthquake mislocations? Xiaodong Song Department of Geology,
More informationLab 1: Plate Tectonics April 2, 2009
Name: Lab 1: Plate Tectonics April 2, 2009 Objective: Students will be introduced to the theory of plate tectonics and different styles of plate margins and interactions. Introduction The planet can be
More informationRelocation of a global earthquake data set with a 3-D velocity model
Chapter 8 Relocation of a global earthquake data set with a 3-D velocity model We have relocated a global earthquake data set of 450,000 events contained in the International Seismological Centre (ISC),
More information9th Workshop on Three-Dimensional Modelling of Seismic Waves Generation, Propagation and their Inversion
1965-36 9th Workshop on Three-Dimensional Modelling of Seismic Waves Generation, Propagation and their Inversion 22 September - 4 October, 2008 Tomography and Active Tectonics in Kanto, Japan Francis T.
More informationFigures S1 S4 show the measurements taken from the synthetic vespagrams where a)
Figures S1 S4 show the measurements taken from the synthetic vespagrams where a) is the differential travel time versus the Dʺ discontinuity height, b) is the differential travel time versus δv S, c) is
More informationMantle plume tomography
Chemical Geology 241 (2007) 248 263 www.elsevier.com/locate/chemgeo Mantle plume tomography Guust Nolet a, Richard Allen b, Dapeng Zhao c, a Department of Geosciences, Princeton University, Princeton,
More informationMantle and crust anisotropy in the eastern China region inferred from waveform
Earth Planets Space, 53, 159 168, 2001 Mantle and crust anisotropy in the eastern China region inferred from waveform splitting of SKS and PpSms Takashi Iidaka 1 and Fenglin Niu 2 1 ERI, University of
More informationData Repository Item
Data Repository Item 2009003 An abrupt transition from magma-starved to magma-rich rifting in the eastern Black Sea Donna J. Shillington, Caroline L. Scott, Timothy A. Minshull, Rosemary A. Edwards, Peter
More information4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution
4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution PI: Lijun Liu Department of Geology, University of Illinois at Urbana-Champaign Corresponding author: Lijun Liu, ljliu@illinois.edu
More informationTravel time tomography of the uppermost mantle beneath Europe
Chapter 3 Travel time tomography of the uppermost mantle beneath Europe We have obtained a detailed P and S model of the uppermost mantle beneath Europe using regional travel time data based on the ISC
More informationMapping the Upper Mantle Discontinuities beneath China with Teleseismic Receiver Functions
Mapping the Upper Mantle Discontinuities beneath China with Teleseismic Receiver Functions Xuzhang Shen 1, 2, 3 Huilan Zhou 1 and Hitoshi Kawakatsu 2 1 Laboratory of Computational Geodynamics, Graduate
More informationEarthq Sci (2011)24:
Earthq Sci (2011)24: 27 33 27 doi:10.1007/s11589-011-0766-6 Receiver function study of the crustal structure of Northeast China: Seismic evidence for a mantle upwelling beneath the eastern flank of the
More informationHigh-precision location of North Korea s 2009 nuclear test
Copyright, Seismological Research Letters, Seismological Society of America 1 High-precision location of North Korea s 2009 nuclear test Lianxing Wen & Hui Long Department of Geosciences State University
More informationMECHANISM OF THE 2011 TOHOKU-OKI EARTHQUAKE: INSIGHT FROM SEISMIC TOMOGRAPHY
Proceedings of the International Symposium on Engineering Lessons Learned from the 2011 Great East Japan Earthquake, March 1-4, 2012, Tokyo, Japan MECHANISM OF THE 2011 TOHOKU-OKI EARTHQUAKE: INSIGHT FROM
More informationSearch for deep slabs in the Northwest Pacific mantle (subduction/lithosphere/residual spheres/lower mantle)
Proc. Nati. Acad. Sci. USA Vol. 86, pp. 8602-8606, November 1989 Geophysics Search for deep slabs in the Northwest Pacific mantle (subduction/lithosphere/residual spheres/lower mantle) HUA-WE ZHOU AND
More information4 Deforming the Earth s Crust
CHAPTER 7 4 Deforming the Earth s Crust SECTION Plate Tectonics BEFORE YOU READ After you read this section, you should be able to answer these questions: What happens when rock is placed under stress?
More informationData Repository Item For: Kinematics and geometry of active detachment faulting beneath the TAG hydrothermal field on the Mid-Atlantic Ridge
GSA Data Repository Item: 2007183 Data Repository Item For: Kinematics and geometry of active detachment faulting beneath the TAG hydrothermal field on the Mid-Atlantic Ridge Brian J. demartin 1*, Robert
More informationSCIENCE CHINA Earth Sciences
SCIENCE CHINA Earth Sciences RESEARCH PAPER September 2011 Vol.54 No.9: 1386 1393 doi: 10.1007/s11430-011-4177-2 Crustal P-wave velocity structure of the Longmenshan region and its tectonic implications
More informationPlate Tectonics. entirely rock both and rock
Plate Tectonics I. Tectonics A. Tectonic Forces are forces generated from within Earth causing rock to become. B. 1. The study of the origin and arrangement of Earth surface including mountain belts, continents,
More informationHaruhisa N. (Fig. + ) *+ Graduate School of Environmental Studies, Nagoya University, Furo-cho, Chikusa-ku, Nagoya.0. 20*+ Japan.
/- (,**2) 0,+/,,+ Source Mechanism and Seismic Velocity Structure of Source Region of Deep Low-frequency Earthquakes beneath Volcanoes: Case Studies of Mt Iwate and Mt Fuji Haruhisa N AKAMICHI + +3 (Fig
More informationMapping the upper mantle discontinuities beneath China with teleseismic receiver functions
Earth Planets Space, 60, 713 719, 2008 Mapping the upper mantle discontinuities beneath China with teleseismic receiver functions Xuzhang Shen 1,2,3, Huilan Zhou 1, and Hitoshi Kawakatsu 2 1 Laboratory
More information1. I can describe evidence for continental drift theory (e.g., fossil evidence, mountain belts, paleoglaciation)
Science 10 Review Earth Science Vocabulary asthenosphere continental drift theory converging plates diverging plates earthquakes epicentre fault hot spot inner core lithosphere mantle mantle convection
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
MODELING TRAVEL-TIME CORRELATIONS BASED ON SENSITIVITY KERNELS AND CORRELATED VELOCITY ANOMALIES William L. Rodi 1 and Stephen C. Myers 2 Massachusetts Institute of Technology 1 and Lawrence Livermore
More informationDetailed structure of the upper mantle discontinuities around the Japan subduction zone imaged by receiver function analyses
Earth Planets Space, 57, 5 14, 25 Detailed structure of the upper mantle discontinuities around the Japan subduction zone imaged by receiver function analyses Takashi Tonegawa 1, Kazuro Hirahara 1, and
More informationWhole Earth Structure and Plate Tectonics
Whole Earth Structure and Plate Tectonics Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 4/5/2017 14:45 We Discuss Whole Earth Structure and Plate
More informationFull file at
Chapter 2 PLATE TECTONICS AND PHYSICAL HAZARDS MULTIPLE-CHOICE QUESTIONS 1. What direction is the Pacific Plate currently moving, based on the chain of Hawaiian Islands with only the easternmost island
More informationWhat Forces Drive Plate Tectonics?
What Forces Drive Plate Tectonics? The tectonic plates are moving, but with varying rates and directions. What hypotheses have been proposed to explain the plate motion? Convection Cells in the Mantle
More informationEarth as a planet: Interior and Surface layers
Earth as a planet: Interior and Surface layers Bibliographic material: Langmuir & Broecker (2012) How to build a habitable planet Internal structure of the Earth: Observational techniques Seismology Analysis
More informationStructural heterogeneity in the megathrust zone and mechanism of the 2011 Tohoku oki earthquake (Mw 9.0)
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl048408, 2011 Structural heterogeneity in the megathrust zone and mechanism of the 2011 Tohoku oki earthquake (Mw 9.0) Dapeng Zhao, 1 Zhouchuan
More informationAn Introduction of Aleutian Subduction Zone. Chuanmao Yang, Hong Yang, Meng Zhang, Wenzhong Wang 2016/04/29
An Introduction of Aleutian Subduction Zone Chuanmao Yang, Hong Yang, Meng Zhang, Wenzhong Wang 2016/04/29 Outline General Introduction Formation history Structure from seismic study geochemical features
More informationImaging the Gutenberg Seismic Discontinuity beneath the Oceanic Crust of the North American Plate
Imaging the Gutenberg Seismic Discontinuity beneath the Oceanic Crust of the North American Plate Robbie Burgess 11-25-15 Dr. Nicholas Schmerr GEOL 394 1 1. Abstract: The lithosphere-asthenosphere boundary
More informationVertical coherence in mantle heterogeneity from global seismic data
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl049281, 2011 Vertical coherence in mantle heterogeneity from global seismic data L. Boschi 1,2 and T. W. Becker 3 Received 11 August 2011; revised
More information5: ABSOLUTE PLATE MOTIONS & HOTSPOTS
5-1 5: ABSOLUTE PLATE MOTIONS & HOTSPOTS 1 Relative motions between plates are most important In some applications important to consider absolute plate motions, those with respect to the deep mantle ABSOLUTE
More informationPlate Tectonics. Structure of the Earth
Plate Tectonics Structure of the Earth The Earth can be considered as being made up of a series of concentric spheres, each made up of materials that differ in terms of composition and mechanical properties.
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
SASIA3D: A 3-D CRUST AND UPPER-MANTLE VELOCITY MODEL OF SOUTH ASIA DERIVED FROM JOINT INVERSION OF P-WAVE TRAVEL TIMES AND SURFACE-WAVE DISPERSION DATA William L. Rodi 1 and Delaine T. Reiter 2 Massachusetts
More informationShear-wave anisotropy beneath the Ryukyu arc
LETTER Earth Planets Space, 61, 1197 1202, 2009 Shear-wave anisotropy beneath the Ryukyu arc Mamoru Nakamura and Ai Yamamoto Department of Physics and Earth Sciences, Faculty of Science, University of
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/1131692/dc1 Supporting Online Material for Localized Temporal Change of the Earth s Inner Core Boundary This PDF file includes: Materials and Methods Figs. S1 to S3
More informationThorne, Garnero, Jahnke, Igel, McNamara Supplementary Material - 1 -
Supplementary Material S1. Bin Location Map Location of 2.5 2.5 bins for which data was obtained to infer the ULVZ distribution. The preferred ULVZ model is drawn (yellow areas). Gray shaded regions indicate
More informationBRIEF ARTICLE THE AUTHOR
BRIEF ARTICLE THE AUTHOR 1 2 THE AUTHOR S Pd K Pd S mantle OC IC CMB Figure 1 Figure 1. Illustration of the SPdKS / SKPdS ray-paths with sub-segments labeled. SPdKS is an SKS that intersects the source-side
More informationDOUBLE-PLANED STRUCTURE OF INTERMEDIATE- DEPTH SEISMIC ZONE AND THERMAL STRESS IN THE DESCENDING PLATE. (Received December 20, 1983)
J. Phys. Earth, 31, 329-347, 1983 DOUBLE-PLANED STRUCTURE OF INTERMEDIATE- DEPTH SEISMIC ZONE AND THERMAL STRESS IN THE DESCENDING PLATE Hiroyuki HAMAGUCHI, Kazuhiko GOTO, and Ziro SUZUKI Geophysical Institute,
More informationThree-dimensional structure of the African superplume from waveform modelling
Geophys. J. Int. (5) 161, 283 294 doi: 1.1111/j.1365-X.5.258.x Three-dimensional structure of the African superplume from waveform modelling Sidao Ni, 1,2 Don V. Helmberger 2 and Jeroen Tromp 2 1 School
More informationGeological structure analysis in Central Java using travel time tomography technique of S waves
Journal of Physics: Conference Series PAPER OPEN ACCESS Geological structure analysis in Central Java using travel time tomography technique of S waves To cite this article: I. R. Palupi et al 2016 J.
More informationAverage shear-wave velocity structure of the Kamchatka peninsula from the dispersion of surface waves
Earth Planets Space, 52, 573 577, 2000 Average shear-wave velocity structure of the Kamchatka peninsula from the dispersion of surface waves N. M. Shapiro 1, A. V. Gorbatov 2, E. Gordeev 3, and J. Dominguez
More informationGeneral Oceanography Geology 105 Expedition 10 - Whole Lotta Shakin' Goin' On
General Oceanography Geology 105 Expedition 10 - Whole Lotta Shakin' Goin' On Name Not attempting to answer questions on expeditions will result in point deductions on course workbook (two or more blank
More informationA magnitude 7.8 earthquake has occurred km (63.3 mi) ESE of Suva, Fiji at a depth of km (378 miles).
A magnitude 7.8 earthquake has occurred 101.8 km (63.3 mi) ESE of Suva, Fiji at a depth of 608.6 km (378 miles). There is no risk of a tsunami from an earthquake at this depth. Images courtesy of Google
More informationMagnitude 7.0 PERU. This region of the Andes is a sparsely populated area, there were no immediate reports of injuries or damage.
A magnitude 7.0 earthquake occurred in southeastern Peru on Friday about 27 kilometers northeast of the town of Azángaro, Peru, near the border with Bolivia. The earthquake occurred at a depth of 257.4
More informationScience 10 PROVINCIAL EXAM STUDY BOOKLET. Unit 4. Earth Science
Science 10 PROVNCAL EXAM STUDY BOOKLET Unit 4 Earth Science Student nstructions 1. Ensure that you have blank paper and a Data Booklet. 2. Record all answers on a separate piece of paper. 3. Answer keys
More informationSuperplume Project: Towards a new view of whole Earth dynamics
RESEARCH NEWS Earth Planets Space, 51(1), i v, 1999 Superplume Project: Towards a new view of whole Earth dynamics Mizuho Ishida 1, Shigenori Maruyama 2, Daisuke Suetsugu 3, Shigeru Matsuzaka 4, and Takao
More informationSelected Seismic Observations of Upper-Mantle Discontinuities
Selected Seismic Observations of Upper-Mantle Discontinuities Peter Shearer IGPP/SIO/U.C. San Diego August 31, 2009 Earthquake Research Institute Interface Depth vs. Publication Date Most depths are sampled
More informationPlate Tectonics. By Destiny, Jarrek, Kaidence, and Autumn
Plate Tectonics By Destiny, Jarrek, Kaidence, and Autumn .The Denali Fault and San Andreas Fault - The San Andreas Fault is a continental transform fault that extends roughly 1300 km (810 miles) through
More informationThe Structure of the Earth and Plate Tectonics
The Structure of the Earth and Plate Tectonics Agree or Disagree? 1. The Earth if made up of 4 different layers. 2. The crust (where we live) can be made of either less dense continental crust or the more
More informationIcelandic Lithosphere & Upper Mantle Structure. Natalia Solomatova
Icelandic Lithosphere & Upper Mantle Structure Natalia Solomatova Location of Iceland maps.google.com Lithosphere Thickness Anomalously thick lithosphere beneath Greenland and Baltic shield may be due
More informationOnline material: Updated 3-D tomography model and plotting script with examples
Model Update January 2010: Upper Mantle Heterogeneity beneath North merica from Traveltime Tomography with Global and usarray Transportable rray Data Scott urdick, Robert D. van der Hilst, Frank L. Vernon,
More informationConvergent plate boundaries. Objective to be able to explain the formation and key features of these zones.
Convergent plate boundaries Objective to be able to explain the formation and key features of these zones. Destructive plate margins When plates collide due to convection currents/slab pull in the asthenosphere
More information