Search

Recent observations beneath Kanto, Japan have shown that seismic activity and seismic attenuation within the overlying continental plate change with time due to drainage caused by slow‐slip events (SSEs) along the upper boundary of the Philippine Sea plate. However, associated changes in rock properties have not been investigated. In this study, we estimate frequency‐dependent shear‐wave anisotropy to provide a detailed insight into the structural change associated with drainage. We perform shear‐wave splitting analysis in frequency ranges of 1–4, 2–6, and 4–8 Hz for 306 earthquakes that occur during September 2009–August 2021 and recorded at the Metropolitan Seismic Observation network. Obtained time differences between fast and slow S waves (delay time) range from almost zero to 0.16–0.18 s, exhibiting spatio‐temporal variation and frequency dependence. The fast S‐wave polarization directions are generally consistent with the direction of the maximum horizontal compressional axis in the study region, which suggests that the observed anisotropy is probably caused by the NE–SW‐oriented fractures developed under the regional stress field. The temporal variation in delay times is correlated with SSEs activity with a lag time of 0.0–0.1 year. Furthermore, comparisons between observed frequency‐dependent delay times and numerical calculation of fracture‐induced anisotropy suggest that the average fracture radius is almost constant (0.30–0.35 m) over time but fracture density temporally varies from 0.025 to 0.035. We infer that the fracture density is probably enhanced by opening of the NE–SW‐oriented fractures during the upward migration of fluids that are expelled from the plate interface. Plain Language Summary: Measurement of shear‐wave polarization anisotropy characterizes the intensity and orientation of fractures within the rocks beneath the Earth's surface. Furthermore, its frequency dependence provides important information for determining the fracture size and density. This study reveals that the strength of anisotropy shows temporal variations associated with the inferred periodic drainage from the subducting slab. We discuss a possible mechanism of the structural change focusing on the frequency dependence of the observed anisotropy and propose that the opening of fractures occurs immediately above the subducting slab following the periodic drainage. The obtained results will provide important constraints on fluid‐rock interaction above the subducting plate interface. Key Points: Shear‐wave splitting analysis in Kanto, Japan suggests that seismic anisotropy is orientated sub‐parallel to the regional stressAnisotropy above the plate interface is enhanced with a lag of 0.0–0.1 years from the occurrence of slow‐slip eventsThe observed frequency‐dependent anisotropy suggests that fractures open responding to drainage from plate interface [ABSTRACT FROM AUTHOR]

The major- and trace-element compositions of amphiboles in andesite from Quaternary Yufu Volcano, northeastern Kyushu, Japan were analysed to investigate the generation processes of andesitic magma from Yufu Volcano. The amphiboles in andesite from Yufu volcano can be divided into two groups based on major-element composition: pargasite and magnesio-hornblende. To estimate temperature, pressure, and major- and trace-element compositions of melts in equilibrium with amphiboles, we used the recently proposed methods that can calculate temperature, pressure, major element compositions, and partition coefficients of trace-element between amphibole and melt using only the major-element compositions of amphibole. The estimated temperature, pressure, and major-element composition of melt in equilibrium with the amphibole phenocrysts indicate that each group crystallised under different conditions. These differences suggest that two magma chambers at different depths existed beneath Yufu Volcano and that the andesitic magma of Yufu Volcano was formed by mixing of the two magmas. The trace-element compositions of melts in equilibrium with the pargasite and magnesio-hornblende, estimated by applying the partition coefficients calculated from major-element compositions of amphibole to trace-element compositions of amphiboles, indicate magma derived from slab melt and the partial melting of crustal material, respectively. Because magma is a mixture of minerals and melt, we estimate the chemical compositional ranges of the two end-member magmas on the Y versus SiO2 diagram from the mixing relationship between amphibole and estimated melt, as well as phenocrysts of plagioclase, clinopyroxene, and orthopyroxene. The overlap of the estimated compositional range with the trend of whole-rock composition represents the chemical compositions of the end-members of magma mixing, yielding estimates of the mafic (SiO2 ≈ 45 wt %) and felsic (SiO2 ≈ 68 wt %) end-member magmas. Furthermore, we estimate the concentrations of other elements in the end-member magmas by substituting the estimated SiO2 concentrations of the magmas into linear regression equations between the whole-rock contents of other elements and SiO2. The trace-element compositions of the mafic and felsic end-member magmas, as estimated in this study, have similar features to those of gabbroids and Cretaceous granitic rocks, respectively, that are presumed to lie beneath Yufu Volcano. These similarities could be explained by the possibility that the compositions of the end-member magmas were influenced by basement rocks. [ABSTRACT FROM AUTHOR]

The stability and breakdown of mineral phases in subducting slabs control the cycling of trace elements through subduction zones. Stability of key minerals and the partitioning of trace elements between these minerals and liquid phases of interests have been charted by natural sample analysis and experimental constraints. However, systematic study from arc front to far back arc has rarely shown that the expected geochemical variations of the slab liquid are actually recorded by natural samples. Complexities arise by uncertainties on the nature of the slab component (melts, fluids and supercritical liquids), source heterogeneities and transport processes. Using data from olivine-hosted melt inclusions sampled along and across the NE Japan and southern Kurile arcs, we demonstrate that experimentally and thermodynamically constrained phase stabilities in subducted materials indeed control the trace element signatures as predicted by these models and experiments. The main reactions that can be traced across arc are progressive breakdown of light rare earth element-rich accessory phases (e.g. allanite), enhanced dehydration of the lithospheric mantle (serpentine breakdown) and changes in the nature of the slab component. This work elucidates subduction zone elemental cycling in a well-characterized petrogenetic setting and provides important constraints on the interpretation of trace element ratios in arc magmas in terms of the prograde metamorphic reactions within the subducting slab. [ABSTRACT FROM AUTHOR]

Knowledge of the state of stress in subducting slabs is essential for understanding their mechanical behavior and the physical processes that generate earthquakes. Here, we develop a framework which uses a high-resolution focal mechanism catalog to determine the change in the position of the neutral plane before and after the M9 Tohoku-oki earthquake to determine that the deviatoric stress within the slab at intermediate depths must be very low (~1 MPa). We show that by combining the static stress calculated from coseismic slip distributions with the stress orientations before and after the mainshock, we can determine the full deviatoric stress tensor within the subducting slab at intermediate depths. These results preclude earthquake source mechanisms that require large background driving stresses, favoring a mechanically weak subducting slab, thus providing quantitative constraints on the physical processes that generate intermediate-depth earthquakes. [ABSTRACT FROM AUTHOR]

The subduction of the Philippine Sea (PHS) plate along the Nankai Trough in in southwest Japan is a relatively recent process compared with subduction along the Japan Trench in northeast Japan. However, the tectonic evolution of the PHS plate along the Nankai Trough is still controversial and not fully understood. There are several competing hypotheses based on different estimates for the time variations of convergence rate and plate age. Our study employs numerical modelling of subduction in order to evaluate the slab evolution for the last 15 Myr and aims to evaluate each tectonic scenario against the present-day slab geometry along a profile passing through the Shikoku and Chugoku regions. The modelling strategy involves a parameter study where subduction initiation and various subduction parameters are analyzed in terms of subduction geometry evolution. Two-dimensional visco-elasto-plastic numerical simulations of spontaneous bending subduction predict that convergence rate and plate age variations play an important role in the evolution of subduction geometry. Modeling results after 15 Myr of evolution reveal that the tectonic model based on a high convergence rate between ~ 15 Ma and ~ 3 Ma produces a slab geometry that agrees well with the observed present-day slab shape specific for the Shikoku and Chugoku regions. [ABSTRACT FROM AUTHOR]

Boron is an effective tracer of fluid processes in subduction zones. High B and δ11B in arc magmas require efficient B transfer from the slab to magma source regions. The Higashi-akaishi metaperidotite body in the Sanbagawa high-pressure belt, Japan, is composed of locally serpentinized mantle wedge peridotites exhumed in a subduction channel. Cores of coarse-grained primary mantle olivine have 1-4 µg/g B, enriched compared to typical mantle olivine, and δ11B of --10%c to --1%o, consistent with incorporation of fluids from dehydrating slab at --90-120 km depth. Rims of primary mantle olivine as well as olivine neoblasts have even higher B (5-20 µg/g) and higher δ11B (--8%c to +2%o) due to incorporating slab fluids at depths of --70-100 km. Antigorite, formed below 650 °C, shows comparable δ11B and B contents as olivine rims. The data show that olivine is capable of scavenging significant amounts of B from fluids by diffusion and recrystallization at sub-arc pressures and temperatures. Considering the large amount of olivine in the mantle wedge, transport of slab-derived material to magma sources requires processes with minimal interaction with mantle peridotite, such as intensely channelized fluid flow or ascent of mélange diapirs, and limited porous fluid flow. [ABSTRACT FROM AUTHOR]

The geometry of the Philippine Sea slab (PHS) subducting beneath the Japanese islands has been imaged to 400 km depth beneath the Kyushu and Chugoku regions, whereas the PHS slab geometry beneath the Hokuriku region has only been determined to ~ 140 km depth, thereby indicating a large east–west asymmetry in the slab subduction. However, geologic evidence suggests that there was symmetrical east–west seafloor spreading along the axis of the Kinan seamount chain when the Shikoku basin was an active spreading center in the PHS plate. This inconsistency suggests that the PHS slab should be present beneath the Hokuriku region. Here we perform P-wave travel-time tomography across central Japan and conduct a two-dimensional plate subduction numerical simulation that reproduces the dual subduction of the PHS and Pacific (PAC) plates to elucidate the PHS slab geometry beneath central Japan. The tomography results reveal a high-velocity anomaly at ~ 150–250 km depth that extends from Wakasa bay to Noto peninsula and a slab window beneath the Hokuriku region. The numerical simulation results suggest that the PHS slab may have torn when it collided with the PAC slab, with the once leading edge of the PHS slab now present along the upper surface of the PAC slab beneath Noto peninsula. These results indicate that the PHS slab exists at ~ 250 km depth beneath the Hokuriku region, although it has been torn owing to its collision with the PAC slab, with this tear propagating westward to form a triangular slab window beneath the Hokuriku region. [ABSTRACT FROM AUTHOR]

Dehydration of subducting oceanic lithosphere (slab) induces Li‐isotope fractionation between the fluid and the slab, suggested by the δ7Li variation (∼10‰) in exhumed subduction complexes. Given that arc magmas represent melt of the supraslab mantle, a large δ7Li variation is anticipated for arc volcanic rocks. However, the δ7Li values in these rocks are mostly homogeneous within the range of mid‐ocean ridge basalts (+1.6 to +5.6‰). The lack of a subduction‐related δ7Li signature has been explained by (a) homogenization by mixing of different magma sources, (b) loss of Li from the slab via dehydration, or (c) homogenization by diffusive exchange of slab‐derived Li and the mantle. The Chugoku district in SW Japan is an ideal place to study the process responsible for Li‐isotope variation in arc magmas, since the Chugoku volcanic rocks show large δ7Li variation (−1.9 to +7.4‰). High δ7Li values (+6.3 to +7.4‰) are found in some high‐Sr andesites and dacites (adakites) whereas low δ7Li values (−1.0 to −0.1‰) are found in high‐Mg andesites. The parental magmas of these rocks have been sourced from subducted oceanic crust and sediments, respectively, with various extents of the interaction with wedge mantle. The limited extents of Li isotope modification are indicated by the similarity of the δ7Li values of these rocks and their supposed sources. The models for a slab dehydration and a diffusive exchange between slab‐derived melt and mantle demonstrate that the δ7Li signatures of the sources can be preserved in the adakites if they ascent rapidly in mantle. Plain Language Summary: Many of Earth's volcanoes occur near the trenches, where one plate subducts beneath another. These volcanoes are formed as a result of the melting of mantle that contains materials of subducted plates and the transportation of melts through the mantle. A quantitative understanding of these important processes has been hampered by the lack of geochemical tools to examine consecutive chemical reactions beneath the volcanoes. Here we conducted a lithium isotope investigation of the volcanic rocks from the Chugoku district in Southwest Japan. This district hosts the volcanoes of andesites and dacites which have been derived from the melting of materials in the subducted plate. We found that these volcanic rocks preserve Li‐isotope compositions of the original materials when they melted. This in turn suggests that these magmas ascent rapidly (<300,000 years) in the mantle. Key Points: Volcanic rocks in the Chugoku district in Southwest Japan show large δ7Li variation (−1.9‰ to +7.4‰)Li isotopes fractionate to the limited extent during arc magma productionSlab‐derived δ7Li features can be preserved if slab‐derived fluids/melts have distinct δ7Li values from mantle and ascend rapidly [ABSTRACT FROM AUTHOR]

During moderate to deep (35–260 km) earthquakes within the Philippine-sea slab along the Ryukyu arc, distinctive later phases after S are observed across the Japanese archipelago for epicentral distances from 1500 to 2200 km, producing anomalous amplification of ground motion in central and northern Japan. Broad-band observations show that these later phases have a faster apparent wave speed (7–9.5 km s−1) than S , and are dominant at low frequencies (0.05–1 Hz), indicating S -wave reflections returned from the upper mantle with strong attenuation for high frequencies. Numerical simulation of seismic wave propagation in 3-D models including subduction of both the Pacific (PAC) and Philippine-sea (PHS) plates reveal the origin of these reflections. The triplicated S wave front from the 410 km discontinuity undergoes wide-angle reflection at the top of the Pacific slab bringing strong amplitudes in a narrow band of epicentral distance near 1500 km. Also, wide-angle S reflections at the 660 km discontinuity are reinforced by refraction when travelling through the high-wave speed Pacific slab; these arrive at the surface beyond 2000 km epicentral distance. The characteristics of the deep mantle reflections and the distribution of large ground motion across Japan due to the reflections are strongly dependent on the source depth and distance to the Ryukyu earthquakes. Tomographic imagery for western Japan indicates loss of the high-wave speed signal of the PAC slab below 200 km depth; the character of the large S reflections and pattern of enhanced ground motions favours a model in which the PAC slab is thinned rather than fully broken. The peculiar pattern of ground motion from Ryukyu earthquakes can be useful for constraining deep slab structures that are difficult to identify based on tomography. [ABSTRACT FROM AUTHOR]

Repeating earthquakes, or repeaters, which repeatedly rupture the same location and release the strain energy, are interpreted as repeated ruptures of an isolated asperity patch surrounded by a stable sliding regime. While repeaters are frequently observed along plate‐boundary faults, the occurrence of repeater sequences are reported in various tectonic settings. Here we systematically investigate whether repeaters occur in the continental crust and subducting slabs beneath the Japanese Islands and show a prevalence of repeaters in both locations. The crustal and intraslab repeaters show the power‐law decay of seismicity rates that is identical to that along plate‐boundary faults and occur on well‐defined fault planes coincidentally with non‐repeating earthquakes. These observations suggest that repeaters share a common generation process independent of the tectonics regime and are not distinct from non‐repeaters in their source locations and occurrence times. We thus infer that repeaters and non‐repeaters are both generated by ruptures of locked asperities in response to stress loading by aseismic slip in the surrounding stable regime. Repeater sequences can be observed when aseismic slip is large enough to rupture an isolated asperity more than twice, while no repeaters are generated when the amount of aseismic slip is not sufficient to cause multiple ruptures at an isolated asperity. Since an earthquake cannot be triggered unless the amount of the aseismic stress loading exceeds the strength of asperities at the fault, hidden aseismic slip probably occurs repeatedly and frequently during an interseismic period. Plain Language Summary: The majority of earthquakes occur at three places in the Earth; along plate boundaries where two plates are adjacent, in downgoing plates (slabs), and in continental plates. However, whether the same generation mechanism works for the earthquakes in the different tectonic regimes is still debatable. Earthquakes with very similar waveforms have been recognized since ∼1900, and we recognize at present that many of such similar earthquakes are repeating earthquakes, or repeaters. The repeaters are caused by repeated ruptures of an isolated, frictionally locked asperity patch by stress loading due to aseismic slip in the surrounding stable regime. This study reveals that these repeaters occur prevalently in both the continental crust and slabs beneath the Japanese Islands. We find that there are no substantial differences in the hypocenter locations and occurrence time between repeaters and non‐repeaters. These observations suggest that all earthquakes share a common generation mechanism. We hypothesize that aseismic slip plays an essential role in stress loading to frictionally locked areas and an earthquake occur when the amount of aseismic slip exceeds the strength of faults. Our idea, although it needs further verification, will help better understanding of how frequent aseismic slip occurs and how earthquakes are triggered. Key Points: Repeating earthquakes are more common in the continental crust and subducting slabs than previously consideredRepeaters are not distinct from non‐repeaters in their source locations and occurrence timesRepeaters and non‐repeaters are triggered by ruptures of locked asperities by stress loading of aseismic slip [ABSTRACT FROM AUTHOR]

Earthquake‐cycle deformation at subduction zones is not only induced by cyclical faulting processes but also impacted by the Earth structure and rheology. Oversimplifying the first‐order Earth structure in models leads to biased estimations of fault kinematics and Earth rheology and hence misinterpretation of subduction‐zone dynamics. Here we first use synthetic models to investigate the impacts of an elastic slab on viscoelastic postseismic and interseismic deformation and on determining afterslip, locking distribution, and mantle viscosity. We confirm that the presence of a slab causes anchored viscous flow in the mantle wedge and less extensive viscoelastic deformation in the overriding plate than the depth‐dependent layered models without a slab during the interseismic and postseismic phases. By assuming the predicted viscoelastic deformation from the model with a slab as the synthetic true observations, layered models obtain a higher viscosity and more slip on the downdip of the seismogenic zone than the true model during both phases. These biased estimates are found varying with the displacement component, spatial and temporal coverages of the used observations. We then conduct two real‐world case studies at the Japan subduction zone and find that the layered model obtains the biased fault kinematics and mantle viscosity that are similar to the synthetic modeling results, indicating that the slab signature may have been observed in the geodetic data. Our findings highlight the importance of including a slab in modeling the earthquake‐cycle deformation and thoroughly exploring the sensitivities of surface deformation, especially the verticals, on the deep Earth structures. Plain Language Summary: A descending oceanic plate (i.e., slab) is much colder and stronger than its surrounding mantle rocks at depth. However, such differences between the slab and its surroundings are often not considered in computational subduction models. In this research, we investigate how different the surface motion is predicted by models with and without a slab and what are the consequences in quantifying the mantle properties at depth as well as the slip on the subduction fault. These two questions are studied by first running purely synthetic models and then applying the models to the real motion data from Japan. We find that if using models without a slab, the mantle is predicted to be stronger and the fault to slip more at depth than the true models. The extents of the deviations depend partially on when and what locations that the data are collected for the study. To better understand the fault kinematics and mantle properties at subduction zones, we suggest to revisit the previous no‐slab models and to correctly include a slab in future models. Key Points: Effects of an elastic slab on the viscoelastic postseismic and interseismic deformation are explored using synthetic models and real casesIgnoring the slab in models leads to a broader horizontal deformation pattern and a contrasting, even reversed, vertical patternLayered models without a slab prefer a higher mantle viscosity and more slip near downdip seismogenic edge in post‐ and interseismic phases [ABSTRACT FROM AUTHOR]

It is commonly assumed that intermediate-depth seismicity is in some way linked to dehydration reactions inside subducting oceanic plates. Although there is growing evidence that the hydration state of an oceanic plate is controlled by its structure and degree of faulting, we do not have a quantitative understanding of this relationship. Double seismic zones offer the possibility of investigating changes in oceanic-plate hydration not only along strike but also with depth beneath the slab surface. To quantify the impact of oceanic-plate structure and faulting on slab hydration and intermediate-depth seismicity, with a focus on the genesis of double seismic zones, we correlate high-resolution earthquake catalogs and seafloor maps of ship-based bathymetry for the northern Chilean and Japan Trench subduction zones. The correlations show only a weak influence of oceanic-plate structure and faulting on seismicity on the upper plane of the double seismic zone, which may imply that hydration is limited by slow reaction kinetics at low temperatures 5-7 km below the seafloor and by the finite amount of exposed wall rock in the outer-rise region. These factors seem to limit hydration even if abundant water is available. Seismicity in the lower plane is, in contrast, substantially enhanced where deformation of the oceanic plate is high and distributed across intersecting faults. This likely leads to an increase in the volume of damaged wall rock around the faults, thereby promoting the circulation of water to mantle depths where serpentinization is faster due to elevated temperatures. Increased lower-plane seismicity around subducting oceanic features such as seamounts or fracture zones may also be caused by enhanced faulting around these features. Our results provide a possible explanation for the globally observed presence of rather homogeneous upper-plane seismicity in double seismic zones as well as for the commonly patchy and inhomogeneous distribution of lower-plane seismicity. [ABSTRACT FROM AUTHOR]

Small‐scale intraslab heterogeneity is well documented seismically in multiple subduction zones, but its nature remains elusive. Previous efforts have been mostly focusing on the scattering strength at intermediate depth (<350 km), without constraining its evolution as a function of depth. Here, we illustrate that the inter‐source interferometry method, which turns deep earthquakes into virtual receivers, can resolve small‐scale intraslab heterogeneity in the mantle transition zone. The interferometric waveform observations in the Japan subduction zone require weak scattering (<1.0%) within the slab below 410 km. Combining with previous studies that suggest high heterogeneity level (∼2.5%) at intermediate depth, we conclude that the small‐scale intraslab heterogeneity weakens as slabs subduct. We suggest that the heterogeneities are caused by intraslab hydrous minerals, and the decrease in their scattering strength with depth reveals processes associated with dehydration of subducting slabs. Plain Language Summary: Deep earthquakes often generate surprisingly lengthy surface shaking, due to wave trapping and guiding along laminar small‐scale scatterers within subducting slabs. The nature of these scatterers is not well understood, although they have been observed in multiple subduction zones above 350 km. Here, by turning some deep earthquakes into virtual sensors closer to our target, we find that the small‐scale scatterers within the slab core fade substantially as slab subducts, from ∼2.5% above 350 km to <1% below 410 km. The fading signal suggests that the scatterers are caused by heterogeneous hydration of the slab in the outer rise that decreases as the slab core dehydrates. Key Points: We apply inter‐source interferometry to probe the small‐scale intraslab heterogeneity in the mantle transition zone beneath the Japan SeaThe deep intraslab scattering level needs to be weaker than at intermediate depth to fit the high‐frequency interferometric waveformsThe decrease of intraslab scattering with depth reveals the processes associated with slab dehydration [ABSTRACT FROM AUTHOR]

Ar-Ar ages, and petrographical and geochemical characteristics of pyroclastics and an overlying lava from Teshima Island, southwest Japan are presented. Although previous geological and age data suggested Teshima pyroclastics were products of magmatism > 3 my prior to lava flows of Setouchi volcanic rocks generated in association with southward migration of the southwest Japan arc sliver during opening of the Sea of Japan backarc basin at 15 Ma, the present results led to the conclusion that a sequence of Setouchi volcanism, induced by slab melting and subsequent melt-mantle reactions, produced both pyroclastics and lava at 14.6-14.8 Ma. This age is oldest among those reported so far and may represent the timing of onset of characteristic Setouchi magmatism immediately posterior to and hence as a result of the mega-tectonic event including rotation of the southwest Japan arc sliver. [ABSTRACT FROM AUTHOR]

Volcanism following the initiation of subduction is vital to our understanding of this specific magma-generation environment. This setting is represented by the first development of the Izu–Bonin–Mariana arc system as subduction commenced along the Western Pacific margin in the Eocene. A new collection of volcanic rocks recovered from the islands and exposed crustal sections of the Bonin Ridge spans the first 10 Myr of arc evolution. An elemental and radiogenic isotope dataset from this material is presented in conjuction with new 40Ar/39Ar ages and a stratigraphic framework developed by a detailed mapping campaign through the volcanic sections of the Bonin Islands. The dating results reveal that both the locus and type of magmatism systematically changed with time in response to the progressive sinking of the slab until the establishment of steady-state subduction at around 7–8 Ma. Following initial mid-ocean ridge basalt (MORB)-like spreading-related basalt magmatism, volcanic centres migrated away from the trench and changed from high-Si boninite to low-Si boninite or high-Mg andesite, then finally tholeiitic or calcalkaline arc magma. Subducting pelagic sediment combined with Pacific-type igneous ocean crust dominates the slab input to the shallow source of high-Si boninites at 49 Ma, but high-precision Pb isotope data show that this sediment varies in composition along the subducting plate. At around 45 Ma, volcanism switched to low-Si boninite and the pelagic sediment signature was almost entirely replaced by volcanic or volcaniclastic material originating from a HIMU ocean island source. These low-Si boninites are isotopically consistent with a slab component comprising variable proportions of HIMU volcaniclastic rocks and Pacific MORB. In turn, this signature was replaced by a Pacific MORB-dominated flux in the post 45 Ma tholeiite and calcalkaline volcanic rocks. Notably, each change in slab-derived flux coincided with a change in the magma type. Fluctuations in the slab-derived geochemical signature were superimposed on a change in the mantle wedge source from highly depleted harzburgite to a depleted MORB-type mantle-type source. In turn, this may correspond to the increasing depth of the leading edge of the slab through this 5 Myr period. [ABSTRACT FROM AUTHOR]

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]

The 1891 Nobi earthquake was the largest historic intraplate earthquake in Japan. The rupture origin corresponds to a zone of high strain rate inferred from geodetic data. To understand the geologic setting of this event, we deployed temporary seismic stations in the area. Shear-wave splitting analysis was performed using data from temporary and permanent seismic stations, revealing significant lateral variations in polarization directions. Polarization directions of NE-SW, ESE-WNW, and ENE-WSW were observed in the northeastern, central, and southwestern parts of the study area, respectively. The NE-SW- and ENE-WSW-aligned polarizations are consistent with the subduction directions of the Philippine Sea plate and Pacific plate, respectively; thus, shear-wave splitting in the northeastern and southwestern regions of the study area is likely caused by mantle wedge anisotropy, a consequence of mantle flow caused by the subducting oceanic slabs. However, the ESE-WNW orientations observed in the central Chubu Region are inconsistent with the subduction direction of either slab. Regions of low seismic velocity and low resistivity have been reported in the inferred position of the mantle wedge; these heterogeneities are thought to be caused by fluid rising from the dehydrated oceanic slabs. Thus, the ESE-WNW polarization in central Chubu could be a consequence of structural heterogeneities created by fluid to the crust from the mantle. The presence of crustal fluid is closely related to weakening, and the faults responsible for the 1891 Nobi earthquake are located just above the anisotropic region. Because fluids in the crust weaken the surrounding rock, this could explain the occurrence of the 1891 Nobi earthquake. [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]

Abstract: Oceanic plates subducted at trenches penetrate into the deep mantle, and encounter a structural boundary at a depth of 410km where olivine, the dominant element of mantle rocks, transforms into a higher density form wadsleyite. This transformation may be delayed within the coldest core of subducting plates (slabs) due to kinetic effects, and it has been suggested that metastable olivine may persist deeper than 410km. Using high density seismic array data in Japan, we show the direct image of the structure corresponding to this metastable olivine wedge (MOW) beneath southwest Japan. Numerical simulation of a subducting slab, including the kinetic effect of water (H2O) on the olivine–wadsleyite transformation, indicates that the presence of the imaged MOW requires an insignificant amount of water (less than 100wt. ppm) be present in the slab mantle, thus a deep dry cold slab. We infer that the transportation of water into the deep mantle occurs along the top surface of the subducting slab, but no significant amount within the slab itself. We also demonstrate that a numerical simulation including the kinetics of 660-km phase transformation can reconcile the observed deep depression of the 660-km discontinuity with a gentle Clapeyron slope. [Copyright &y& Elsevier]

Abstract: The petrological consequences of deep subhorizontal deflection (“stagnation”) of subducting slabs should affect both apparent seismic velocity structures and slab morphology. We construct kinematic thermal models of stagnant slabs and perform thermodynamic modeling of the consequent perturbation of high-pressure phase transitions in mantle minerals, focusing upon Japan as our study area. We calculate associated thermo-petrological buoyancy forces and bending moments which (along with other factors such as viscosity variations and rollback dynamics) may contribute to slab deformation. We consider effects of variations in depth of stagnation, post-stagnation dip angle, phase transition sharpness, transition triplication due to multiple intersection of geotherms with phase boundaries, and potential persistence of metastable phases due to kinetic hindrance. We also estimate seismic velocity anomalies, as might be imaged by seismic tomography, and corresponding seismic velocity gradients, as might be imaged by receiver-function analysis. We find that buoyant bending moment gradients of petrological origin at the base of the transition zone may contribute to slab stagnation. Such buoyancy forces vary with the depth at which stagnation occurs, so that slabs may seek an equilibrium slab stagnation depth. Metastable phase bending moment gradients further enhance slab stagnation, but they thermally decay after ∼600–700km of horizontal travel, potentially allowing stagnant slabs to descend into the lower mantle. Stagnant slabs superimpose zones of negative velocity gradient onto a depressed 660-km seismic discontinuity, affecting the seismological visibility of such features. Seismologically resolvable details should depend upon both stagnation depth and the nature of the imaging technique (travel-time tomography vs. boundary-interaction phases). While seismic tomography appears to yield images of stagnant slabs, discontinuity topography beneath Japan resolved by migrated receiver functions appears to be consistent with slab penetration of the transition zone. However, model slabs which bottom around ∼780–810km and then bend upwards by a few degrees can match both the tomographic and receiver-function images. [Copyright &y& Elsevier]

Abstract: The metaperidotites found in the Happo-O''ne region of central Japan are mostly lherzolites to harzburgites with subordinate dunites. Their protoliths originated as a series of refractory residues by near-fractional melting, 15–25% melting in the spinel field, which is confirmed by a HREE melting model of the apparently primary clinopyroxene. These metaperidotites display U-shaped primitive mantle (PM)-normalized REE patterns (0.02–0.5 times PM) and are highly enriched in LILE (0.2–20 times PM) relative to HFSE (<0.2 times PM), providing evidence for mantle-wedge metasomatism. In-situ analyses confirmed that clinopyroxenes and tremolite are enriched in fluid-mobile elements (B, Sr, Pb, Li, Cs, Ba and Rb; 0.1–100 times PM) coupled with depletion of HFSE (Ta, Hf, Th, Zr, Ti and Nb; <0.7 times PM)+U; these chemical features of clinopyroxenes and tremolite are similar to those of their whole rocks, and reflect slab-fluid metasomatism. We recognized three stages of clinopyroxenes that are different in morphology and REE patterns, but show the same major-element chemistry. Nearly euhedral clinopyroxene (Cpx1) displays a spoon-shaped REE pattern (0.2–3 times PM), reflecting the first stage (stage 1) of metasomatism. Acicular clinopyroxene (Cpx2) after tremolite+olivine shows a U-shaped REE pattern (0.06–1 times PM) (stage 2), whereas fine-grained clinopyroxene (Cpx3) derived from orthopyroxene shows very low concentrations of REE (stage 3). These three successive stages of retrogressive clinopyroxene formation possibly denote a retrogressive chemical change of the involved fluids during multi-stage metasomatism. The fluid of stage 1 is highly enriched in LREE (0.2–2 times PM), MREE, Pb, Sr, Li and Rb relative to the stage-3 fluid that is very low in MREE, Pb, LREE and Sr, but has high Na, Mn, Ba, B and Cs. The stage-2 fluid, which is low in LREE (0.06–1 times PM), Li and Pb, and very low in Cs, Rb and Hf relative to the fluid of stage 1, is considered to be a transitional fluid between stages 1 and 3. The retrogressive change of the fluid composition is possibly equivalent to a transversal change of the slab-derived fluid within the mantle wedge due to a change of metamorphic condition, depth of the subducting slab and an episode of metasomatism during the exhumation of the Happo-O''ne peridotites. [Copyright &y& Elsevier]

Abstract: Hypocenter relocations of earthquakes in the Pacific slab have shown an anomalous deepening of a seismic belt in the upper-plane of the double seismic zone at depths of 80–120km beneath the Hokkaido corner, while it is located at depths of 70–90km in the surrounding Tohoku and eastern Hokkaido areas. Seismic tomographic inversions performed beneath the Hokkaido corner have shown that a low-velocity zone having seismic velocities of crust materials exists in the mantle wedge just above the Pacific slab and makes direct contact with the upper surface of the Pacific slab. These observations suggest that: 1) the low-velocity zone just above the Pacific slab is part of the subducted Kuril forearc sliver. 2) The contact with the subducted, and so relatively cold, sliver materials prevents the mantle wedge from heating the Pacific slab and causes a lower temperature condition in the Pacific slab crust beneath the Hokkaido corner than in the surrounding areas. 3) As a result, a delay of eclogite-forming phase transformations occurs and enhances the local deepening of the seismic belt in the slab crust there. [Copyright &y& Elsevier]

Abstract: We examine the upper mantle and transition zone beneath the western Pacific using multiple ScS reverberations. A low-velocity zone (LVZ) is found above the 410-km discontinuity oceanward of the subducting Honshu slab at an average depth of 356km, with a thickness that ranges from 50 to 75km assuming the LVZ continues to the 410-km discontinuity, which is locally elevated. The low-velocity region is evident in previous tomographic studies, and our results suggest that the anomaly is best explained by a layer of partial melt. The layer may be entrained from above by subduction or produced in situ by the combined effects of water and temperature. A self-consistent model that explains local P-wave velocities (Obayashi, M., Sugioka, H., Yoshimitsu, J., Fukao, Y., 2006. High temperature anomalies oceanward of subducting slabs at the 410-km discontinuity. Earth and Planetary Science Letters 243, 149–158) and our observations calls for a maximum temperature anomaly of ∼150°C and a resulting maximum olivine water content after melting of 0.100wt. %. [Copyright &y& Elsevier]

Abstract: Both the global and regional P wave tomographic studies have revealed significant deep structural heterogeneities in subduction zone regions. In particular, low-velocity anomalies have been observed beneath the descending high-velocity slabs in a number of subduction zones. The limited resolution at large depths and possible trade-off between the high and low velocities, however, make it difficult to substantiate this feature and evaluate the vertical extent of the low-velocity structure. From broadband waveform modeling of triplicated phases near the 660-km discontinuity for three deep events, we constrained both the P and SH wave velocity structures around the base of the upper mantle in northeast Asia. For the two events beneath the Southern Kurile, the rays traveled through the lowermost transition zone and uppermost lower mantle under the descending Pacific slab. Our preferred models consistently suggest normal-to-lower P and significantly low SH velocities above and below the 660-km discontinuity extending to about 760-km depth compared with the global IASP91 model, corroborating previous observations for a slow structure underneath the slab. In contrast, both high P and SH velocity anomalies are shown in our preferred models for the Japan subduction zone region, likely reflecting the structural feature of a slab stagnant above the 660-km discontinuity. The velocity jumps across the 660-km discontinuity were found to be on average 4.5% and 7% for P and S waves under the south Kurile, and 3% and 6% under the Japan subduction zone. The respective velocity contrasts in the two regions are consistent with mineralogical models for colder slab interior and hotter under-slab areas. Based on mineral physics data, the depth-averaged ∼1.5% P and ∼2.5% SH velocity differences in the depth range of 560–760km between the two regions could be primarily explained by a 350–450K temperature variation, although the presence of about 0.5–1wt.% water might also contribute to the subtle velocity variations near the base of the transition zone in the Southern Kurile. From our modeling results, we speculate that the slow structure in the Southern Kurile may be correlated to the low-velocity zone observed previously around the 410-km discontinuity under Northern Honshu. If this is the case, both may be associated with a thermal anomaly rooted in the lower mantle beneath the subduction zone in northeast Asia. [Copyright &y& Elsevier]

Abstract: We examined seven ultramafic xenoliths from 1~3 Ma alkali olivine basalt reefs near the Eurasian continent and one sample of the host alkali basalt to identify the mantle wedge material and to constrain the origin and evolution of mantle beneath SW Japan. Six xenoliths are from Kurose and one xenolith is from Takashima, northern part of the Kyushu islands, SW Japan. The Sr and Nd isotopic ratios vary from 0.70416 to 0.70773 and from 0.51228 to 0.51283, respectively. The Kurose and Takashima xenoliths have higher Sr isotopic ratios and lower Nd isotopic ratios than those of the peridotite xenoliths from the other arc settings such as Simcoe and NE Japan. The Kurose xenoliths have less radiogenic Os isotopic ratios (187Os/188Os=0.123–0.129) than the primitive upper mantle (PUM) estimate and limited variation compared to the other arc xenoliths. Their Os isotope compositions are rather similar to the ultramafic xenoliths from NE and east China. In addition, the samples of the Kurose and Takashima xenoliths plot along a mixing line between ultramafic xenoliths from SE and NE China and a slab component in Sr–Nd–Os isotopic space. Our results suggest that fragments of continental lithospheric mantle from the China craton may exist beneath Kurose and Takashima after the Sea of Japan expansion when the Japanese islands were rifted away from the Eurasian continent during Miocene. Later magmatism due to subduction of the Philippine Sea Plate beneath the SW Japan arc around 15 Ma ago may have introduced fluids or melts derived from slab component, interpreted to be oceanic sediments rather than altered oceanic crust, that possibly modified the original composition of the lithospheric mantle sampled by the peridotite xenoliths from Kurose and Takashima. [Copyright &y& Elsevier]

Subducted slabs of oceanic lithosphere below the western Pacific tend to be stagnant in the transition zone with poorly known mechanical properties. Typical examples are the Izu-Bonin and Japan slabs that meet each other to form a cusplike junction beneath southwest Japan. Here, we show that these two slabs are torn apart at their junction when they bend to flatten over the 660-kilometer discontinuity, as is expected from a simple geometric argument. We present three lines of evidence for this ongoing slab tear. [ABSTRACT FROM AUTHOR]

Tungsten is strongly incompatible during magmatic processes and is fluid mobile in subduction zones. Here we show that W isotope fractionation in arc lavas provide a powerful new tool for tracing slab dehydration and melting in subduction zones. Geochemically well characterized, representative arc-lavas from three subduction zones were chosen for this study to evaluate W isotope fractionation under different sub-arc conditions. Arc-lavas from SW Japan are produced by subducting a young, hot slab, and lavas from the volcanic front and rear arc of the Sangihe and Izu arcs are produced during subduction of a cold slab. The heaviest W isotope compositions (δ 184 W ∼ 0.06 ‰) are observed in fluid-rich samples from the volcanic fronts of the Sangihe and Izu arcs. With increasing distance from the volcanic front, more melt-rich samples are characterized by progressively lighter W isotope compositions. Enriched alkali basalts from SW Japan, thought to be the product of mantle melting at a slab tear, and adjacent shoshonites have the lightest W isotope compositions (δ 184 W ∼ 0 ‰). The correlation of W isotope fractionation with various indices of fluid release (e.g., Ce/Pb, Ba/Th) suggests that the heavy W isotope signatures record fluid recycling near the volcanic front due to dehydration of the subducted slab. Upon release of the heavy W, the residual slab preferentially retains isotopically light W, which is released during subsequent melting of drier lithologies in hot subduction zones, such as SW Japan. These data suggest that W isotopes can be used as a tracer of slab dehydration, potentially helping to determine the onset of cold subduction zone magmatism and hence, modern-style plate tectonics. • The W stable isotope composition of bulk earth is equal to chondritic values. • Cold, fluid-rich subduction zones have heavy W stable isotope compositions. • Hot, dry subduction zones have light W stable isotope compositions. • Tungsten stable isotopes have the potential to be used as a tracer for slab dehydration. [ABSTRACT FROM AUTHOR]

Seaward-concave subduction forearcs are controversial settings that exhibit remarkable along-strike variations of surface uplift and rock exhumation, the causes of which are still a matter of debate. While short-term (10-10^2 a) forearc uplift relates to locking of the megathrust, it is unclear which mechanism produces long-term (10^4-10^6 a) permanent surface uplift. In this work we use Coulomb-wedge solutions and background forearc seismicity (not megathrust) in order to investigate permanent deformation of the overriding plate in Peru-Chile and northern Japan. According to the Coulomb-wedge theory, the observed along-strike variations of long-term forearc uplift may reflect changes of (i) slab dip angle, (ii) wedge material strength and/or (iii) plate interface friction. Here, we explore the latter hypothesis through a comparison of the spatial distribution of basal friction (inferred from the taper geometry) and the pattern of background forearc seismicity (obtained from the International Seismological Center global earthquake catalog).In Peru-Chile, concordant patterns of cumulative seismic moment and inferred basal strength suggest that (i) at the curvature apex, low interplate friction is needed to explain the forearc morphology; (ii) permanent vertical displacement of the proximal wedge is likely accomplished through repeated non-megathrust seismic events; (iii) the wedge is largely in equilibrium. Conversely, in northern Japan a complex relationship between long-term deformation and plate-interface processes is deduced. [ABSTRACT FROM AUTHOR]

The article discusses research that provides tomographic evidence for the development of a vertical tear in a slab underlying southwest Japan and shows that the plate is still rupturing. Measurements of stress revealed by active seismic sources were correlated with these images. The concept that slabs rupture as they move through the upper mantle has been used to explain geological observations in the past. In addition to knowing the location of an earthquake, high-resolution tomographic reconstructions and data collection networks are needed to distinguish slab rupture from sharp flexure.