Search

Marginal parts of a plate and subducting slab can play important roles in geodynamics. This is because in areas where a plate interacts with other plates or with the mantle thermal, geochemical, and mechanical interactions are expected. The Philippine Sea (PHS) slab that subducts beneath the Japan arcs has such an edge. To examine the relationship between arc magmatism and the slab edge in the transition zone from Northeast Japan to Central Japan, we investigated isotopic systematics of the regional volcanic rocks, incorporating data from literature and new data for five isotopic ratios of Sr, Nd, and Pb. The new data included major element compositions of 22 samples from the back-arc area, and 5 isotopic ratios for 6 samples selected from Pleistocene to early Quaternary epochs. Consequently, several findings were determined based on the spatial variation of the isotopic ratios and the estimated amount of slab-derived fluid: (1) the amount of fluid derived from the two subducting slabs (i.e., the Pacific slab and the PHS slab) decreases northward from a significantly high value (~ 5 wt% fluid added to the source mantle), away from the seismically determined edge of the PHS slab; (2) the proportion of the PHS component in the total slab-derived fluid also decays northward; and (3) the PHS component spreads to the north beyond the seismically determined edge of the PHS slab. These observations strongly suggest that the existence of an aseismic PHS slab beneath southernmost parts of Northeast Japan delivers the PHS component to the arc magmatism. As was indicated by previous geodynamical studies, subduction of the PHS and PAC slabs may generate suction force towards the corner of mantle wedge, which might account for the large amount of fluid near the seismically determined slab edge as described in (1) above. [ABSTRACT FROM AUTHOR]

The Pacific and Okhotsk plates are converging at a rate of about 9 cm/a, thereby forming the Japan Trench subduction zone. Two large earthquakes occurred in the offshore region of Fukushima on 13 February 2021 (M w 7.1) and 16 March 2022 (M w 7.3). The occurrences of the two earthquakes are close in space and time, thus constituting the Fukushima-Oki earthquake doublet. The joint inversions of the onshore and offshore strong-motion data and the teleseismic P waves reveal that both earthquakes belong to intraslab events, and the unilateral ruptures of the 2021 and 2022 events propagate southwestward and northeastward, respectively. Both events are dominated by thrust motions, and major fault slips are mainly localized away from the rupture initiating area. The geometries and locations of the two seismogenic faults suggest that the Fukushima-Oki doublet arises from the reactivations of two preexisting bending-related faults inside the subducting Pacific plate, which is driven by the downdip compressional stress regime caused by plate unbending. Furthermore, the regional tectonic settings indicate a complicated seismogenic environment for the Fukushima-Oki doublet, in which the dehydration embrittlement and seamount subduction are likely to play significant roles. Additionally, the distinctive rupture patterns of both events are also associated with the variation of fluid content and the small-scale structure of a slab slice enclosed by the two seismogenic faults within the subducting plate, which may control the heterogenous distributions of fault slips and aftershocks. The analysis of Coulomb failure stress changes suggests that the 2021 event and the 2022 foreshock collectively increase the stress loading, which may facilitate the occurrence of the 2022 mainshock, and the seismicity of the seismic belt consisting of the three intraslab sequences following the great 2011 Tohoku-Oki earthquake is also likely to remain active in the future. • Source rupture processes of two Fukushima-Oki earthquakes in northeastern Japan. • Joint inversions of onshore and offshore strong-motion and teleseismic waveforms. • Two unilateral ruptures separately propagating southwestward and northeastward. • Both intraslab events caused by reactivations of preexisting outer-trench faults. • Seismogenesis attributed to heterogenous fluid content and regional tectonics. [ABSTRACT FROM AUTHOR]

We investigate the 3-D P and S wave velocity structure of the East Japan forearc using a large number of high-quality arrival-time data of local earthquakes recorded by a permanent seafloor seismic network (S-net) and the land-based Kiban seismic network. Geometries of the subducting slab boundary and the Conrad and Moho discontinuities under the forearc are refined by referring to results of previous studies including active-source seismic soundings and seismicity located by the S-net. The obtained P and S wave velocity images are quite similar to each other, suggesting that they reflect robust features of the 3-D seismic structure. The subducting Pacific slab exhibits high velocities, whereas significant low-velocity (low-V) anomalies are revealed in the overlying Okhotsk plate beneath the forearc, which may reflect accretionary materials containing abundant fluids. The subducting oceanic crust atop the Pacific slab is revealed clearly as a thin low-V layer extending from the trench axis down to ~100 km depth under the volcanic front. Our results of earthquake relocation using the 3-D velocity model indicate that the double seismic zone in the subducting Pacific slab occurs from the trench axis down to ~180 km depth. The occurrence and focal mechanisms of the double seismic zone are caused by slab deformation during its subduction, such as slab bending or unbending, as well as other processes including hydration, dehydration, compositional variations, and phase changes in the subducting slab. [Display omitted] • Seafloor S-net and Kiban network data are used to study the 3-D forearc structure. • The Conrad, Moho and Pacific slab geometries in the East Japan forearc are refined. • P and S wave velocity images are quite similar, reflecting robust structural features. • The subducting oceanic crust is revealed down to ~100 km depth under the volcanic front. • The intraslab double seismic zone occurs from the trench axis down to ~180 km depth. [ABSTRACT FROM AUTHOR]

Zones of fluid-bearing anomalies occur in volcanic regions, on a plate interface, or in the mantle wedge of a subduction zone. The various causes include slab dehydration, serpentinization or metamorphic reactions in different thermal regimes. To reveal the impact factors that might have contributed to slab melting and arc magmatism behind the Japan Trench, we determined high-resolution three-dimensional (3-D) structures of seismic velocities (Vp, Vs), Vp/Vs ratio and temperatures. The temperatures (T) were calculated from the inverted seismic structures, and two petrological models of the upper mantle: a peridotite assemblage and a pyrolite assemblage, which differ from the results of the thermal simulation. The multi-parameter structures obtained for Vp, Vs, Vp/Vs and T show good consistency in the upper mantle behind the Japan Trench. The temperatures derived from the peridotite assemblage were found to match the mean geotherm and simulated temperature of the upper mantle more closely than those from the pyrolite model. A layer approximately 10 km thick with low-Vp and Vs, high-Vp/Vs, and slightly high-T anomalies was observed on the upper boundary of the subducting Pacific slab. This distinctive anomalous layer is interpreted as partial melting of the oceanic crust due to deep-seated metamorphic reactions whose characteristics are dependent on the source of fluids, mineral composition, and the thermal regime. Such a process enriches the peridotite content of the basalt underneath the island arc in the mantle wedge. Localized zones of significantly low-Vp and low-Vs, high-Vp/Vs and high-T perturbations were imaged in the mantle wedge under the active volcanoes, suggesting partial melting of peridotite-rich mantle material to produce tholeiitic magma. The present study demonstrates that fluids released from slab dehydration, mineral composition and the thermal regime play crucial roles in both arc magmatism, and slab melting in the subduction zone. [Display omitted] • High resolution seismic and thermal structures are imaged in the NE Japan subduction zone. • Slab melting is controlled by the source of fluids and thermal regime behind the Japan Trench. • Fluids and mineral composition play a key role in slab melting and arc magmatism. [ABSTRACT FROM AUTHOR]

The motivation for this work was to model the crustal thickness, and the geometry of the Mohorovičić (Moho) discontinuity, the crust-upper mantle boundary around the Japanese archipelago. The western Pacific Ocean is one of the most eminent subduction zones in the world. The subduction zone complexity is due to alteration within the geometry of the stagnant slab. We compiled a three-dimensional density distribution model for subduction plates around the Japan Islands, and estimated the depth to Moho using the constrained non-linear inversion of gravity data derived from the latest Gravity field and Ocean Circulation Explorer, or GOCE, model known as GOCO06s. The study had three steps, starting with the data corrections followed by the inversion, and ending with the modeling and comparing of the estimated results with previous models and studies. In this study, we used available receiver function results, seismic tomography models, global crustal thickness, and geological and tectonic information as constraints. The model estimated the crustal thickness, surface geometry of the Moho, and geometry of the stagnant slab in the Kuril Trench, Japan Trench, Nankai Trough, and Izu-Bonin Trough. In gravity inversion, we used 34 km and 600 kg m−3 as the Moho depth reference and density contrast, respectively. Results showed that the crustal thickness ranged from 14 to 43 km with a − 0.78 mean misfit for Moho depths derived from seismic P-wave analysis. The calculated gravity data had a 0.16 mGal misfit with the observed data. The normalized analytical signal helped to sharply delineate the different plates and their edges. Three cross-sections were conducted to model the stagnant slab along the subduction zone of the Eastern Eurasian plate, western Pacific plate, Okhotsk plate, and Philippine Sea plate. The eastern area of the entire region (Pacific plate) had a low crust thickness and high density compared to the western area (Eurasian plate), which had a high crust thickness and low density. • 3D gravity inversion of the Northwestern Pacific Plate at the Japan Islands. • Estimating the Moho depths around the Japan Islands. • 3D density modeling of the Stagnant Slabs and Subduction plates. [ABSTRACT FROM AUTHOR]

• Effects of anisotropic thermal conductivity on the thermal structure are investigated. • Two regions with a high degree of anisotropy are observed in the mantle wedge. • Slab temperature changes only slightly due to the anisotropic thermal conductivity. The effects of anisotropic thermal conductivity due to the crystallographic-preferred orientation (CPO) of olivine in the mantle wedge on the steady-state thermal structure of the Tohoku subduction zone, Northeast Japan, are investigated using geodynamic modeling. Four different olivine crystal fabrics are modeled based on theoretical CPO development in the mantle wedge: random (i.e., isotropic thermal conductivity), A-type, C-type, and E-type fabrics. The thermal conductivity in the mantle wedge is highly anisotropic in the regions just above the descending slab and near the base of the overriding plate. The obtained temperature differences between the anisotropic and isotropic cases are only up to 30 °C within the slab. Therefore, the impact of anisotropic thermal conductivity on the slab temperature may be limited when anisotropy is only considered in the mantle wedge. [ABSTRACT FROM AUTHOR]