2024 JCR Q1
X
Advanced Search
Zhi Wang (2010). Interplate coupling and seismotectonics under the fore-arc regions of Japan. Earthq Sci 23(6): 555-565. DOI: 10.1007/s11589-010-0755-1
Citation: Zhi Wang (2010). Interplate coupling and seismotectonics under the fore-arc regions of Japan. Earthq Sci 23(6): 555-565. DOI: 10.1007/s11589-010-0755-1

Interplate coupling and seismotectonics under the fore-arc regions of Japan

More Information
  • Corresponding author:

    Zhi Wang, e-mail: mike-wang@sohu.com

  • Received Date: 27 Jun 2010
  • Accepted Date: 09 Oct 2010
  • Available Online: 30 May 2022
  • Published Date: 09 Dec 2010
  • Three-dimensional P- and S-wave velocities (vP, vS), Poisson's ratio (σ), crack-density (ε) and bulk-sound velocity (vϕ) structures along the slab upper boundary beneath the fore-arc regions were determined using a large number of high-quality P-wave and S-wave arrival times from both onshore and offshore earthquakes in Japan. The velocity and Poisson's ratio images provide a compelling evidence for a highly hydrated and serpentinized fore-arc mantle and fluid-bearing anomalous low velocity and high Poisson's ratio associated with slab dehydration under the fore-arc areas. Most great thrust earthquakes (M > 7.5) occurred at or close to the high-velocity areas along the slab interface under the fore-arc areas, suggesting strong interplate coupling (asperities) with slab subduction. On the other hand, prominent low-velocity anomalies were revealed along the slab upper boundary, which may reflect weak coupled or decoupled patches (aseismicity) of the plates due to serpentinization of the fore-arc mantle wedge. The crack-density and bulk-sound velocity images, calculated from the corresponding velocity models, indicate that the interplate coupling in northeastern Japan is different from that under central and southwestern Japan owing to differences between the tectonic backgrounds of the subduction system, such as the geological age, thermal regime and dipping angle of the oceanic plates. A comparison between fluid-related anomalies of Japan, Cascadia, Chile, and Costa Rica subduction zones suggests that seismic mantle may be common in fore-arc settings and these reflect similar 3-D seismic structures relatively to fluid liberating processes. We consider that the fluid-bearing anomalies along the interface of the subducting slab, attributing to processes such as slab dehydration and serpentinization of the fore-arc mantle, are mainly contributed to the interplate coupling and the repeated generation of the great thrust earthquakes under the fore-arc regions in Japan.

  • The Islands of Japan is a distinct region characterized by the presence of three main tectonic plates: the Pacific plate, the Eurasian plate and the Philippine Sea (PHS) plate (Figure 1). The Pacific plate subducts WNW under the PHS plate and the Eurasian plates with a ~30° angle of dip and a rate of ~10 cm/a from the Japan-Kuril Trench in Northeastern Japan (; ). The PHS plate is subducting to the WNW under the Eurasian plate at a relative rate of ~5.5 cm/a from the Nankai Trough in Southwestern Japan (; ; , ). In central Japan, the PHS plate is sandwiched by the overlying Eurasian plate and the underlying Pacific slab (; ). The fore-arc regions of Hokkaido, Northeastern Japan, Kanto district and Kyushu in Japan are the sites of important processes associated with the subduction of the Pacific and PHS plates, including serpentinization of the fore-arc mantle, repeated genesis of megathrust earthquakes, slab dehydration and the interplate coupling. Many previous studies were conducted to investigate the structural heterogeneities and their implications for the seismotectonics and arc volcanisms under the land areas of Japan (e.g., ; ; ; ). Few studies of the structural heterogeneities under the fore-arc areas in the Pacific Ocean and the Philippine Sea were done because the offshore earthquakes are poorly located by the routine procedure of the seismic network on the land areas owing to few seismic stations located there. However, more great earthquakes occurred under the offshore areas than onshore. The frequent occurrence of great thrust earthquakes in the fore-arc regions has caused widespread damage to the coastal areas through strong shakings and tsunamis. The spatial distribution of the megathrust earthquakes indicates significant variations under the fore-arc areas along the Islands of Japan (; ; ). We consider that the seismicity and interplate coupling or decoupling along the slab's upper boundary might be correlated with the structural heterogeneities. Therefore, information on the detailed velocity structures together with other physical parameters along the upper boundary of the subducting slab in the fore-arc regions is crucial for a better understanding of the initiation of the interplate coupling, the rupture nucleation of destructive interplate earthquakes and the slab dehydration. However, such information is scanty because the fore-arc areas are generally covered under oceans, and consequently there are few seismic stations available there to record the offshore earthquakes.

    Figure 1. Map showing the tectonic framework around the Islands of Japan (Modified from Kamata and Kodama, 1999). Curved lines show the trenches that represent the major plate boundaries in Japan. Three boxes show the study regions, the Northeastern Japan (Figure 3), the Kanto district in Central Japan (Figure 4) and Kyushu in Southwestern Japan (Figure 5).
    Figure  1.  Map showing the tectonic framework around the Islands of Japan (Modified from ). Curved lines show the trenches that represent the major plate boundaries in Japan. Three boxes show the study regions, the Northeastern Japan (Figure 3), the Kanto district in Central Japan (Figure 4) and Kyushu in Southwestern Japan (Figure 5).
    Figure 3. P-wave, S-wave, Poisson's ratio and crack-density as well as bulk sound velocity images along the upper boundary of the subducting Pacific slab in Northeastern Japan. The location of the plan views along the line shown in the insert map. Blue color denotes high velocity and red color indicates low velocity. Opened circles represent epicenters of large historic interplate earthquakes (M > 6.0) occurred between 830.A.D and 2001. Red triangles denote active volcanoes along the Kuril-Northeastern Japan arc. The contour lines show distribution of rupture areas of the great earthquakes (epicenters as red circles, M > 7.5) estimated by Nagai et al.  (2001). The color scales of the velocity, Poisson's ratio, crack-density and bulk sound velocity are shown horizontally.
    Figure  3.  P-wave, S-wave, Poisson's ratio and crack-density as well as bulk sound velocity images along the upper boundary of the subducting Pacific slab in Northeastern Japan. The location of the plan views along the line shown in the insert map. Blue color denotes high velocity and red color indicates low velocity. Opened circles represent epicenters of large historic interplate earthquakes (M > 6.0) occurred between 830.A.D and 2001. Red triangles denote active volcanoes along the Kuril-Northeastern Japan arc. The contour lines show distribution of rupture areas of the great earthquakes (epicenters as red circles, M > 7.5) estimated by Nagai et al. (2001). The color scales of the velocity, Poisson's ratio, crack-density and bulk sound velocity are shown horizontally.
    Figure 4. P-wave, S-wave, Poisson's ratio, crack-density and bulk sound velocity images along the upper boundary of the subducting Pacific slab under Central Japan. Black triangles indicate active volcanoes and red circles denote the large thrust earthquakes (M > 6.0) occurred during the period from 1923 to 1998 (Utsu, 1999). The color scales of the velocity, Poisson's ratio, crack-density and bulk sound velocity are shown at the bottom of each map.
    Figure  4.  P-wave, S-wave, Poisson's ratio, crack-density and bulk sound velocity images along the upper boundary of the subducting Pacific slab under Central Japan. Black triangles indicate active volcanoes and red circles denote the large thrust earthquakes (M > 6.0) occurred during the period from 1923 to 1998 (). The color scales of the velocity, Poisson's ratio, crack-density and bulk sound velocity are shown at the bottom of each map.
    Figure 5. P-wave, S-wave, Poisson's ratio, crack-density and bulk sound velocity images along the upper boundary of the subducting Philippine Sea slab in Kyushu. Red triangles denote active volcanoes and red circles denote the large thrust earthquakes (M > 6.0) occurred during the period from 1923 to 1998 (Utsu, 1999; Yamaoka et al., 2002). The color scales of the velocity, Poisson's ratio, crack-density and bulk sound velocity are shown at the base of each map.
    Figure  5.  P-wave, S-wave, Poisson's ratio, crack-density and bulk sound velocity images along the upper boundary of the subducting Philippine Sea slab in Kyushu. Red triangles denote active volcanoes and red circles denote the large thrust earthquakes (M > 6.0) occurred during the period from 1923 to 1998 (; ). The color scales of the velocity, Poisson's ratio, crack-density and bulk sound velocity are shown at the base of each map.

    detected the sP depth phase in the offshore area of Northeastern Japan, and then they relocated 690 small offshore earthquakes using theoretical sP-P travel times calculated with a two-dimensional ray-tracing method. In their study only the focal depths of the offshore earthquakes were determined, while the epicenters were fixed. and , , ) applied the same procedure as the previous study to have relocated another 92 offshore earthquakes and then they estimated the velocity structures under the offshore region in Tohoku. However, the seafloor topography of the Pacific Ocean, which has strong influence on the calculation of the ray paths and the travel times, was not considered in their studies. This is because the depth of the seafloor varies from tens to thousands of meters under the offshore regions from the trench to the eastern coast line of Japan.

    Based on the above reasons, we presented a new method to find the bounce point of the sP depth phase using a three-dimensional perturbation method and then relocated 385 offshore earthquakes including both focal depths and epicenters under the offshore regions of Northeast Japan and Hokkaido (). Meanwhile, the influence of the seafloor topography was considered for the first time in their study. Furthermore, taking advantage of a new hypocenter location techniques (; ), we presented an effective hypocenter location method (Master-Event location: MEL) for the offshore earthquakes, outside the areas covered by the seismic networks (, , ). This method makes use of the sP depth phase combined with the extended double-difference location method. To demonstrate the effectiveness of the MEL method, they applied it to both a real data set in Central Japan () and a synthetic data set in Southwestern Japan (). The testing results showed that the hypocenters of the offshore earthquakes that were relocated with the MEL method are consistent with those determined by using the arrivals recorded by the offshore stations.

    In this study the Poisson's ratio (σ), crack-density (ε) and bulk-sound velocity (vϕ) images were calculated from the corresponding velocity models. The calculated parameters are less reliable than the inverted velocity structures. The velocity images and the calculated physical parameters revealed not only the common consensus on interpretational framework of seismic tomography but also the difference of the structural heterogeneities under the fore-arc regions from Northeastern to Southwestern Japan. Noted that the calculated crack-density images are reliable feature at shallow depth under the fore-arc region but would be unreliable at deep depth under the mantle wedge. Here Ⅰ only discuss the crack-density structures beneath the fore-arc region. Those observations suggest that most of the seismological processes, such as the slab dehydration, interplate coupling and the genesis of megathrust earthquakes, are mainly attributed to the subduction conditions, for instance, the geological ages, the dipping angle and temperature of the subducting slab as well as the convergence of the Pacific and the PHS plates.

    The high-resolution 3-D velocity models shown in this study (Figures 35) were calculated using a combined data set including two groups of earthquakes (, , ). One group consists of 27 388 shallow and intermediate-depth earthquakes (depth < 200 km) that occurred beneath the land areas of Japan. These events were located within the seismic network and so they have reliable hypocentral locations with errors smaller than 0.5 km in epicentral location and 2.0 km in depth. The other group includes 3 728 offshore earthquakes at depths of 10-80 km within the fore-arc regions of Japan from Hokkaido to Kyushu. These offshore events have been relocated using the sP depth-phases that were identified from the three-component seismograms recorded by the Hi-net using the MEL method (, ). Figure 2 shows the ray path of the sP depth phase and four examples of the seismograms generated from the offshore earthquakes in Japan. Hypocenters of the offshore events are well constrained by the use of sP depth phase data and the MEL method. Their mislocation errors are generally smaller than 3.0 km in focal depth; comparable to that of the earthquakes occurred within the seismic network (). A total of 645 018 P-wave, 615 943 S-wave and 7 643 sP depth phase arrivals from the selected 27 388 land earthquakes and 3 728 offshore earthquakes have been collected. The picking accuracy is 0.05−0.1 s for P-wave data, 0.05−0.2 s for S-wave data, and 0.1−0.8 s for the sP depth phase data.

    Figure 2. (a) Schematic map showing the ray path of the sP depth phase from an offshore earthquake to the land-based station. (b) Locations of the offshore events and the stations that recorded the arrival times from the offshore events. (c) Four examples of the three-component seismograms that were generated from the offshore events recorded by the land-based stations. From the left to right in each seismogram, three vertical lines indicate the arrivals of P-wave, sP-depth phase and S-wave arrival times, respectively. The epicentral distance of the corresponding stations (alphabetical codes) are shown of the left of each seismogram. The hypocentral parameters are listed above each seismogram.
    Figure  2.  (a) Schematic map showing the ray path of the sP depth phase from an offshore earthquake to the land-based station. (b) Locations of the offshore events and the stations that recorded the arrival times from the offshore events. (c) Four examples of the three-component seismograms that were generated from the offshore events recorded by the land-based stations. From the left to right in each seismogram, three vertical lines indicate the arrivals of P-wave, sP-depth phase and S-wave arrival times, respectively. The epicentral distance of the corresponding stations (alphabetical codes) are shown of the left of each seismogram. The hypocentral parameters are listed above each seismogram.

    The regional tomography method () was applied to the arrival times from the well relocated earthquakes to determine the 3-D velocity images (, , ). We then applied the crack theory (; ; , ) and method of the bulk-sound velocity (; ; ; ) to calculating the crack-density and bulk-sound velocity structures from the tomography velocity models along the slab upper boundary, respectively. Estimated crack parameter together with the Poisson's ratio and bulk-sound velocity images can significantly improve our understanding of structural heterogeneities under the fore-arc regions of Japan.

    Figures 3, 4 and 5 show the plan views of P- and S-wave velocity structures together with the Poisson's ratio, crack-density and bulk-sound velocity images for a two-km-thick layer right above the upper boundaries of the Pacific plate (Figures 3 and 4) and the PHS slab (Figure 5). The spatial distributions and perturbations of P- and S-wave velocity anomalies are similar to one another, although the amplitudes of the P-wave velocity anomalies are less than the S-wave under Northeastern Japan (Figure 3) and the Kanto district (Figure 4). However, this feature of the velocity anomalies is not observed under the Kyushu subduction zone (Figure 5). Beneath the offshore area in Northeastern Japan, two Pand S-wave low-velocity areas with high Poisson's ratio and high crack-density anomalies are imaged (Figure 3). Similarly, high vP/vS zones above the subducting Pacific slab were revealed by double-difference tomography in northern Honshu, which was interpreted as the serpentine dehydration (). The features of the structural heterogeneities under the offshore region are consistent with the previous studies (; ; ). Two slow velocity zones (39°N−41.2°N and 42.0°N−43.0°N) (with low-velocity perturbations of ~2%−6% for P wave and 4%-6% for S wave) are imaged clearly, where no great interplate earthquakes (M > 7.5) had occurred. The low-velocity anomalies are consistent with a proposed seismic moment-release gap (). Strong heterogeneity is revealed at the shallow depths above the subducting slab under the Kanto district (Figure 4). Similar features are observed from the inverted velocity structures in the Kyushu subduction zone. Those low-velocity anomalies could be caused by fluids that are released upwards to the surface by porosity collapse or dehydration of the subducting sediments and the oceanic crust.

    Many large earthquakes occurred between the subducting slab and the overriding continental plate along the interplate thrust zones in Hokkaido, Tohoku, Kanto district and Kyushu of Japan, for example, the 1936 Miyagi earthquake (M7.5), the 1994 Shanriku earthquake (M7.5) and the 2003 Tokachi-oki earthquake (M8.0) in Northeast Japan, the 1941 (M7.2) and 1968 (M7.5) and 1984 (M7.1) Hyuga-nada earthquakes in Kyushu. Figures 3, 4 and 5 demonstrate that more than 92% of the great interplate earthquakes (M≥7.5) occurred around the low-velocity zones, while only a few large interplate earthquakes (M6.0−6.5) occurred in the low-velocity areas, which inferred the existence of weakly coupled patches along the upper boundary of the subducting PHS slab. Similar features were revealed under the fore-arc area of Northeastern Japan (; ), which may reflect a common feature on the spatial distribution of thrust earthquakes and the structural heterogeneities in subduction zones. The mechanism of large thrust earthquakes (Figure 6) calculated by using the Harvard CMT and Hi-Net catalogs shows similar generated mechanism of the large interplate earthquakes in the Northeast Japan and Kyushu but the Kanto district where the Philippine Sea plate is sandwiched by the underlying Pacific and the overlying Eurasian plates (Figure 1). For the subducting slab, the down-dip limit of the great thrust earthquakes often corresponded to the intersection of the thrust with the fore-arc mantle (), which might be explained by aseismic hydrous minerals, such as serpentinite, which are present in the fore-arc mantle wedge and exhibit stable sliding (, , ). Our results show that these interplate thrust earthquakes were located in areas with seismic velocities that are higher than those of their surrounding areas (Figures 3, 4 and 5), which might reflect the futures of strong-coupled asperities.

    Figure 6. Focal mechanism of the thrust earthquakes occurred in the fore-arc regions of Northeast Japan (a), the Kanto district (b) and Kyushu (c). Open stars denote the large thrust earthquakes (Mb > 6.0) during the period from 830 to 2007. Black dots present the earthquakes (M > 5.0) from 1970 to 2007. Gray thin lines indicate the active faults. The gray bold line shows the Japan Trench (a). Gray solid lines show the upper boundaries of the Pacific plate (a and b). The magnitude of the large thrust earthquakes is shown at the bottom.
    Figure  6.  Focal mechanism of the thrust earthquakes occurred in the fore-arc regions of Northeast Japan (a), the Kanto district (b) and Kyushu (c). Open stars denote the large thrust earthquakes (Mb > 6.0) during the period from 830 to 2007. Black dots present the earthquakes (M > 5.0) from 1970 to 2007. Gray thin lines indicate the active faults. The gray bold line shows the Japan Trench (a). Gray solid lines show the upper boundaries of the Pacific plate (a and b). The magnitude of the large thrust earthquakes is shown at the bottom.

    We consider that one possible explanation for the spatial distribution of the large interplate earthquakes is that the lateral heterogeneities, including strong coupling sections or asperities and weak coupled or decoupled patches located along the upper slab boundary, might control the degree and spatial extent of the interplate seismic coupling and the rupture nucleation of the interplate earthquakes. The occurrence of slow and ultra-slow great interplate earthquakes that ruptured the oceanic lithosphere from the seafloor to a depth of several tens of kilometers (; ) may be related to the fluid contents along the slab upper boundary. It is considered that the fluids are released from free pore water contained in the upper oceanic crust and sediments at shallow depth, because of the increasing temperature and compressional pressure, and/or are produced by the extensive dehydration of the subducting slab. We suggest that the low-velocity zones correspond to the weakly coupled or decoupled patches on the slab boundary (; ), which are caused by the presence of fluid-related anomalies.

    The subduction process is also found to be significant in illustrating the role of water and gas (e.g. carbon dioxide, sulphur dioxide) or fluids by influencing the thermo-petrological model of both cool and warm subduction metamorphism (, ). The velocity images along with the Poisson's ratio, crack-density and bulk-sound velocity structures (Figures 35) show that considerable lateral heterogeneities exist on the upper boundary of the Pacific and the Philippine Sea slabs in the fore-arc regions, which might reflect serpentinization of the fore-arc mantle associated with dehydration of the subducting slab. The fore-arc mantle is probably aseismic because of the stable-sliding serpentinite hydrated by water that is expelled upward from the under-thrusting oceanic crust and sediments. Laboratory studies have indicated that serpentinite generally exhibits stable-sliding behavior (e.g., ). Serpentinite is believed to exhibit stable-sliding behavior at the plate subducting rate, thus impeding rupture into the fore-arc mantle (). The mechanism for generating stable-sliding behavior in a subducting slab is the serpentinization of the mantle wedge. Water liberated by the dehydration reactions within the subducting plate might intrude into the overlying fore-arc mantle wedge, resulting in the formation of serpentine minerals and possibly other hydrous minerals, such as talc and brucite (). Because serpentinization will dramatically reduce the seismic velocity, the mechanical strength and density in the fore-arc mantle while increase the Poisson's ratio, magnetization and electrical conductivity (), the serpentinized fore-arc mantle is probably highly heterogeneous and some fluids escape to the surface. Such fluids will migrate into the overlain mantle and decrease the seismic velocity in the fore-arc mantle by the serpentinization process. If the fore-arc mantle is serpentinized, the pattern of convection will be altered in the mantle wedge because of the weak rheology, low seismic velocity and rock density. Such a proposal is evidenced from our estimated seismic velocity, Poisson's ratio and crack-density structures (Figures 35). Strong low-velocity and high Poison's ratio anomalies under the fore-arc regions along the upper slab boundaries are expected to correspond to the serpentinization of the corearc mantle associated with the slab dehydration (Figures 35). Therefore, it is believed that large volumes of aqueous fluids is released upwards by dehydration reactions in the subducting oceanic crust and sediments; the subduction of the oceanic lithosphere will then cool the overlying fore-arc mantle such that low-temperature hydrous serpentine minerals are stable in the fore-arc mantle.

    Several researchers estimated the extent of hydration and serpentinization of mantle wedge for different subduction zones using laboratory experiments, such as Chile margin (serpentinization, 0−12%; water, 0−3.6% by weight), central Japan (serpentinization, 30%−40%, water, 4%−5% by weight), SW Japan (serpentinization, 0−15 %, water, 0−2% by weight), Cascadia (serpentinization, 25%−30%, water, 3%−4% by weight), and Costa Rica (serpentinization, 15%−25 %, water, 3%−4% by weight) (e.g., ; ; ; ; ). Serpentinized mantle wedge near the subducting Cocos plate in Costa Rica implies that rupture during large magnitude earthquakes may not initiate but may propagate along the oceanic slab/continental mantle interface making the downdip limit of the seismic zone beneath northern Costa Rica strain-rate dependent (). The crack-density and bulk-sound velocity images, calculated from the corresponding velocity models, indicate that the interplate coupling in Northeastern Japan is different from that under central and Southwestern Japan owing to differences between the tectonic backgrounds of the subduction system, such as the geological age, thermal regime and dipping angle of the oceanic plates. The geological age of the subducting Pacific plate under northeastern Japan is more than 140 Ma with plate thickness of 80 km, while the geological age of the Philippine Sea plate in southwest Japan is about 50 Ma with plate thickness of 30 km. The dipping angle is about 30° of the Pacific plate in northeast Japan and about 60° of the Philippine Sea plate in southwest Japan. The thermal-petrologic model predicts that the oceanic crust of the Philippine Sea plate under southwest Japan is 300℃−500℃ which is warmer than that beneath Northeast Japan (). These features result in different coupling models of the Pacific plate in northeast Japan and the Philippine Sea plate in southwest Japan.

    For the subducting slab, the down-dip limit of the great thrust earthquakes generally corresponds to the intersection of the thrust zone with the fore-arc mantle (; ), which was explained by the distribution of aseismic hydrous minerals such as serpentinite (; ; ). Our results show that large earthquakes that occurred along the PHS-Eurasian and the upper Pacific boundaries are located in areas with seismic velocities higher than that of surrounding areas (Figures 35), reflecting strongly coupled asperities. We conclude that the spatial distribution of large interplate earthquakes is possibly controlled by lateral heterogeneities, including strongly coupled sections (asperities) and weakly coupled or decoupled patches along the upper slab boundary due to the extensive dehydration process.

    In this study some new features of the structural heterogeneities have been revealed from the 3-D seismic velocity, Poisson's ratio, crack-density and bulk-sound velocity images along the slab upper boundaries under the fore-arc regions of Japan from Hokkaido to Kyushu. The velocity structures together with other physical parameters exhibit strong lateral heterogeneities under the fore-arc regions of Japan from northeast to southwest. On the upper boundary of the subducting Pacific and the PHS slabs, 92% large interplate earthquakes (mb > 7.5) are located at the high-velocity, low Poisson's ratio, high crack-density and high bulk-sound velocity zones, which may represent strong-coupled portions or asperities on the slab upper boundary. Few large interplate earthquakes occur in areas of low velocity, high Poisson's ratio, low crack-density and low bulk-sound velocity areas, representing decoupled areas associated with fore-arc mantle serpentinization caused by fluids from the slab dehydration. The lateral heterogeneities along the slab upper boundary may control the degree and spatial extent of the interplate seismic coupling and the rupture nucleation of the large thrust earthquakes. The crack-density and bulk-sound velocity images indicate that the interplate coupling varies from Tohoku-Hokkaido subduction zone to the regions of Kanto district region and southwest Japan owing to different conditions of the tectonic setting and the subducting slab. We consider that the features of the seismotectonics and interplate coupling are mainly attributed to the slab dehydration controlled by the plate geological ages, the dipping angle and temperature of the subducting slab, and the convergence of the plates. Our study suggests that fluid plays an important role on the interplate coupling and seismotectonics in the subduction zone.

    The author appreciate Prof. John Chen and two anonymous referees who provided thoughtful review comments which improved the manuscript. The waveform and arrival time data used in this work are downloaded from the Hi-net data center of Japan via Internet. Most of the figures were made using GMT (). This work was partially supported by the research grants of the National Natural Science Foundation of China (Nos.40872148 and 40974024) and the New Century Excellent Talents in University, Ministry of Education (NCET2009).

  • Bostock M G., Hyndman R D, Rondenay S and Peacock S M (2002). An inverted continental Moho and serpentinization of the forearc mantle. Nature 417: 536-538. doi: 10.1038/417536a
    Deshon H R, Schwartz S Y, Biklek S L, Dorman L M, Gonzalez V, Protti J M, Dixon T H and Flueh E R (2003). Seismogenic zone structure of the southern middle American Trench, Costa Rica. J Geophys Res 108, doi 10.1029/2002JB002294. doi: 10.1029/2002JB002294
    Hacker B R, Abers G A and Peacock S M (2003a). Subduction factory 1. theoretical mineralogy, densities, seismic wave speeds, and H2O contents. J Geophys Res 108: 2 029, doi: 10.1029/2001JB001127.
    Hacker B R, Peacock S M, Abers G A and Holloway S D (2003b). Subduction factory 2. Are intermediate-depth earthquakes in subducting slabs linked to metamorphic dehydration reactions? J Geophys Res 108: 2 030, doi: 10.1029/2001JB001129, 2003b.
    Hasegawa A, Horiuchi S and Umino N (1994). Seismic structure of the northeastern Japan convergent margin: A synthesis. J Geophys Res 99: 22 195-22 311. doi: 10.1029/93JB02797
    Hyndman R D, Yamano M and Oleskevich D A (1997). The seismogenic zone of subduction thrust faults. The Island Arc 6: 244-260. doi: 10.1111/iar.1997.6.issue-3
    Hyndman R D and Peacock S M (2003). Serpentinization of the forearc mantle. Earth Planet Sci Lett 212: 417-432. doi: 10.1016/S0012-821X(03)00263-2
    Ishida M (1992). Geometry and relative motion of the Philippine Sea plate and Pacific plate beneath the Kanto-Tokai district, Japan. J Geophys Res 97: 489-513. doi: 10.1029/91JB02567
    Iio Y, Sagiya T, Kobayashi Y and Shiozaki I (2002). Water-weakened lower crust and its role in the concentrated deformation in the Japanese Islands. Earth Planet Sci Lett 203: 245-253. doi: 10.1016/S0012-821X(02)00879-8
    Kamiya S and Kobayashi Y (2000). Seismological evidence for the existence of serpentinized wedge mantle. Geophys Res Lett 27: 819-822. doi: 10.1029/1999GL011080
    Kawakatsu H and Seno T (1983). Triple seismic zone and the regional variation of seismicity along the northern Honshu arc. J Geophys Res 88: 4 215-4 230. doi: 10.1029/JB088iB05p04215
    Kanamori H (1971). Seismological evidence for a lithosphere normal faulting the Sanriku earthquakes of 1933. Phys Earth Planet Inter 4: 289-300. doi: 10.1016/0031-9201(71)90013-6
    Kawasaki I, Asai Y and Tamura Y (2001). Space-time distribution of interplate moment release including slow earthquakes and seismogeodetic coupling in the Sanriku-oki region along the Japan trench. Tectonophysics 330: 267-283. doi: 10.1016/S0040-1951(00)00245-6
    Kamata H and Kodama K (1999). Volcanic history and tectonics of the southwest Japan arcs. Isl Rac 8: 393-403.
    Kennett B L N, Wdiyantoro S and van der Hilst R D (1998). Joint seismic tomography for bulk sound and shear wave speed in the Earth's mantle. J Geophys Res 103: 12 469-12 493. doi: 10.1029/98JB00150
    Kennett B L N and Cummins P R (2005). The relationship of the seismic source and subduction zone structure for the 2004 December 26 Sumatra-Andaman earthquake. Earth Planet Sci Lett 239: 1-8. doi: 10.1016/j.epsl.2005.08.015
    Long M and Van der Hilst R D (2006). Shear wave splitting from local events beneath the Ryukyu Arc: Trench-parallel anisotropy in the mantle wedge. Phys Earth Planet Inter 155: 300-312. doi: 10.1016/j.pepi.2006.01.003
    Matsuzawa T, Igrashi T and Hasegawa A (2002). Characteristic small-earthquakes sequence off Sanriku, Northeastern Honshu, Japan. Geophys Res Lett 29(11): 1543, doi: 10.1029/2001GL014632.
    Mishra O P, Zhao D, Umino N and Hasegawa A (2003). Tomography of the northeast Japan forarc and its implications for interplate seismic coupling. Geophys Res Lett 30(16): 1 850, doi: 10.1029/2003GL017736.
    Nakajima J, Matsuzawa T and Hasegawa A (2002). Moho depth variation in the central part of northeastern Japan estimated from reflected and converted waves. Phys Earth Planet Inter 130: 31-47. doi: 10.1016/S0031-9201(01)00307-7
    Nakajima J and Hasegawa A (2007). Subduction of the Philippine Sea slab beneath southwestern Japan: Slab geometry and its relationship to arc magmatism. J Geophys Res 112: B08306, doi: 10.1029/2006JB004770.
    O'Connell R J and Budiansky B (1974). Seismic velocities in dry and saturated cracked solids. J Geophys Res 79: 5 412-5 426. doi: 10.1029/JB079i035p05412
    Oleskevich D A, Hyndman R D and Wang K (1999). The updip and downdip limits to subduction earthquakes: Thermal and structure models of Cascadia, SW Japan, Alaska and Chile. J Geophys Res 104: 14 965-14 991. doi: 10.1029/1999JB900060
    Peacock M S and Hyndman R D (1999). Hydrous minerals in the mantle wedge and the maximum depth of subduction thrust earthquakes. Geophys Res Lett 26: 2 517-2 520. doi: 10.1029/1999GL900558
    Peacock S M (1990). Fluid processes in subduction zones. Science 248: 329-337. doi: 10.1126/science.248.4953.329
    Peacock M S and Wang K (1999). Seismic consequences of warm versus cool subduction metamorphism: Example from southwest and northeast Japan. Science 286: 937-939. doi: 10.1126/science.286.5441.937
    Reinen L A (2000). Seismic and aseismic slip indicators in serpentinite gouge. Geology 28: 135-138. doi: 10.1130/0091-7613(2000)28<135:SAASII>2.0.CO;2
    Schmid C, van der Lee S and Giardini D (2006). Correlated shear and bulk moduli to 1400 km beneath the Mediterranean region. Phys Earth Planet Inter 159: 213-224. doi: 10.1016/j.pepi.2006.07.003
    Seno T, Zhao D, Kobayashi Y and Makamura M (2001). Dehydration in serpentinized slab mantle: Seismic evidence from southwest Japan. Earth Planets Space 53: 861-871. doi: 10.1186/BF03351683
    Seno T (2005). Variation of downdip limit of the seismogenic zone near the Japanese islands: Implications for the serpentinization mechanism of the forearc mantle wedge. Earth Planet Sci Lett 231: 249-262. doi: 10.1016/j.epsl.2004.12.027
    Springer M and Forster A (1998). Heat flow density across the central Andean subduction zone. Tectonophyics 291: 123-139. doi: 10.1016/S0040-1951(98)00035-3
    Umino N, Hasegawa A and Matsuzawa T (1995). sP depth phase at small epicentral distances and estimated subducting plate boundary Geophys J Int 120: 356-366. doi: 10.1111/gji.1995.120.issue-2
    Utsu T (1999). Seismic study: A comprehensive review. University of Tokyo Press, Tokyo, 769-829 (in Japanese).
    Waldhauser F and Ellsworth W L (2000). A double-difference earthquake location algorithm: Method and application to the Northern Hayward Fault, California. Bull Seismol Soc Am 90: 1 353-1 368. doi: 10.1785/0120000006
    Wang Z and Zhao D (2005). Seismic imaging of the entire arc of Tohoku and Hokkaido in Japan by using P-wave, S-wave and sP depth-phase data. Phys Earth Planet Inter 152: 144-162. doi: 10.1016/j.pepi.2005.06.010
    Wang Z and Zhao D (2006a). νP and νS tomography of Kyushu, Japan: New insight into arc magmatism and forearc seismotectonics. Phys Earth Planet Inter 157: 269-285. doi: 10.1016/j.pepi.2006.04.008
    Wang Z and Zhao D (2006b). Suboceanic earthquake location and seismic structure in the Kanto district, central Japan. Earth Planet Sci Lett 241: 789-803. doi: 10.1016/j.epsl.2005.11.015
    Wang Z and Zhao D (2006c). Seismic evidence for the influence of fluids on the 2005 west off Fukuoka prefecture earthquake in southwest Japan. Phys Earth Planet Inter 155: 313-324. doi: 10.1016/j.pepi.2006.01.006
    Wang Z, Zhao D, Mishra O P and Yamada A (2006). Lithospheric heterogeneity and its implications for the lowfrequency tremor in Southwest Japan. Earth Planet Sci Letts 251: 66-78. doi: 10.1016/j.epsl.2006.08.025
    Wang Z, Huang R, Huang J and He Z (2008). P velocity and gradient images beneath the Okinawa Trough. Tectonophysics 455: 1-13. doi: 10.1016/j.tecto.2008.03.004
    Wang Z, Zhao D, Huang R, Tang X and Mishra O P (2009). Structural heterogeneity in Northeast Japan and its implications for the genesis of the 2004 and 2007 Niigata earthquakes in Northeast Japan. Bull Seismol Soc Am 99(6): 3 355-3 373. doi: 10.1785/0120080281
    Wessel P and Smith W (1998). New improved version of the Generic Mapping Tools released. Eos 79: 579. doi: 10.1029/98EO00426
    Yamaoka K, Churei M and Ando M (2002). Inland earthquakes precede the interplate earthquakes in Hyuga-nada, Kyushu, Japan. J Geogr 111: 185-191 (in Japanese with English abstrct). doi: 10.5026/jgeography.111.2_185
    Zhang H, Thurber C H, Shelly A, Ide S, Beroza C and Hasegawa A (2004). High-resolution subducting-slab beneath northern Honshu, Japan, revealed by double-difference tomography. Geology 32(4): 361-364. doi: 10.1130/G20261.2
    Zhao D, Hasegawa A and Horiuchi S (1992). Tomographic imaging of P and S wave velocity structure beneath northeastern Japan. J Geophys Res 97: 19 909-19 928. doi: 10.1029/92JB00603
    Zhao D, Mishra O P and Sanda R (2002). Influence of fluids and magma on earthquakes: Seismological evidence. Phys Earth Planet Int 132: 249-267. doi: 10.1016/S0031-9201(02)00082-1
    Zhao D, Wang Z, Umino N and Hasegawa A (2007). Tomographic imaging outside a seismic network: Application to the northeast Japan arc. Bull Seismol Soc Am 97: 1 121-1 132. doi: 10.1785/0120050256
    Zhao D, Wang Z, Umino N and Hasegawa A (2009). Mapping the mantle wedge and interpolate thrust zone of the northeast Japan arc. Tectonophysics 467: 89-106. doi: 10.1016/j.tecto.2008.12.017
  • Related Articles

Catalog

    Corresponding author: Zhi Wang, mike-wang@sohu.com

    1. On this Site
    2. On Google Scholar
    3. On PubMed

    Figures(6)

    Article views (367) PDF downloads (2) Cited by()

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return