Japan Subduction Zone Research Articles

Three-dimensional ray tracing in P-wave azimuthal anisotropic media

SUMMARY A new 3-D ray-tracing method is developed for P-wave azimuthal anisotropic (AAN) media. We assume anisotropic media with hexagonal symmetry and take advantage of the property that the AAN symmetry axis, the phase velocity vector and the group velocity vector are located in the same plane. The 3-D ray-tracing method that combines the pseudo-bending technique and Snell's law is improved for the AAN media. We compute isotropic (ISO) and AAN rays in synthetic models and an actual 3-D P-wave AAN model of the East Japan subduction zone. The accuracy of our ray-tracing code is evaluated by comparing the ray-path and travel-time differences between the ISO and AAN rays. Our results show that the AAN rays in each model bend in the right direction and satisfy Fermat's principle, so the theory and approximations adopted in the calculations are reasonable. For long rays (>350 km), the ray-path difference between the ISO and AAN rays is > 20 km, and the travel-time difference is > 0.1 s, suggesting that it is necessary and important to take azimuthal anisotropy into account in the 3-D ray tracing.

The effect of antigorite dehydration on velocity structure and water migration in subduction zones

The water migration in subduction zones, primarily driven by the phase transition in hydrous minerals, can give rise to hydrated regions with reduced velocity. A fundamental element in comprehending and deciphering these low-velocity zones revolves around acquiring insights into the stability and elasticity of relevant hydrous minerals. As one of the main water carriers in shallow areas, antigorite can dehydrate to form talc, forsterite, and fluid (talc–bearing peridotites) in deep areas of subduction zones, and then the talc thus serves as one of the minerals that can bring water to the deep Earth. Here, the elasticity of talc up to 24 GPa and forsterite up to 12 GPa are calculated by using the first principles method. The result supposes that the talc structure transforming from talc I to talc II is at a pressure between 6 GPa and 8 GPa, impacting the trend of elastic wave velocity in response to pressure. Furthermore, the elastic wave velocity of forsterite can be significantly affected by iron concentration. Meanwhile, a variation velocity model with antigorite consumption and talc content is set up for talc-bearing serpentinized peridotite based on the elastic properties of talc and forsterite in this study, and antigorite in Wang et al. (2022). The results of our model demonstrate a decrease in the low-velocity anomaly in subduction zones, particularly in deep regions or areas with higher initial serpentinization degrees. The results also suggest that the mode of antigorite dehydration can diminish the estimation of water content transported to depths of subduction zones, such as the Mariana Trench and Northern Japan subduction zones. The mode of antigorite dehydration thus provides a useful tool for constraining the composition, seismic velocity structure, and water migration in subduction zones.

Geodetic Matched Filter Slow Slip Event Detection Along the Northern Japan Subduction Zones

AbstractWe apply a template matching method on GNSS data for stations located in Honshu, Japan, to detect slow slip events associated with the subducting Philippine Sea and Pacific plates during the period from 1997 to 2020. A measure of the minimum detectable moment magnitude is proposed, from which we infer that the method could potentially detect SSEs as small as Mw 5.2 on the westernmost part of the Philippine Sea plate and Mw 6 on the Pacific plate below Honshu eastern coastline. We find 12 slow slip events on the Philippine Sea plate, among which eight are located on the known Boso slow slip event asperity and the four others are located offshore north‐east relative to the Boso SSEs, at the transition with the Pacific plate. We find 9 SSEs on the Pacific plate, mainly on the northern section, offshore Iwate prefecture. A clear gap with no SSEs coincides with the main asperity that broke during the 2011 Tohoku earthquake. Most event locations correlate with low locking areas. We do not find any clear temporal pattern apart from the regular occurrence of the largest Boso SSEs.

Frontal Thrust Ramp‐Up and Slow Earthquakes Due To Underthrusting of Basement High in the Nankai Trough

AbstractRecently, integrated geophysical‐geological surveys in the Nankai subduction zone in Japan have revealed that slow earthquakes repeatedly occur beneath the outer wedge of the forearc. During December 2020 to February 2021, clustered slow earthquakes propagated around the frontal thrust of the accretionary wedge. The frontal thrust ramps up from the basal décollement and slips over trench‐filling sediment along the landward edge of the Nankai trough floor. Here, the Paleo‐Zenisu ridge has been subducted beneath the inner‐outer slope border. In addition, ocean floor topography and geologic structure revealed by seismic reflection surveys completed before 2022 document that the basement of the Philippine Sea Plate beneath the frontal thrust has a seamount and a horst‐like basement high. The northern edge of the basement high is located at the ramp‐up position of the frontal thrust. The 2020–2021 clustered slow earthquakes started at the Paleo‐Zenisu ridge and propagated to the topographic highs beneath the deformation front. Considering that the relative plate convergence between the upper Amurian Plate of the Nankai forearc and the subducting Philippine Sea Plate is ∼6.0 cm/year, the basement high at the deformation front has uplifted the frontal crest of the wedge at an average rate of 2.7–5.7 mm/year for several tens to hundred thousand years. These rates are among some of the highest rock uplift rates measured in the world. The slow earthquakes in the off‐Kumano Nankai Trough in 2020–2021 are a snapshot of a “living” Nankai frontal thrust during the megathrust interseismic period.

Nonergodic Ground-Motion Models for Subduction Zone and Crustal Earthquakes in Japan

ABSTRACT We investigate the nonergodic behavior of ground motions from subduction zone earthquakes and crustal earthquakes in Japan using the National Research Institute for Earth Science and Disaster Resilience strong-motion flatfile together with various reference ergodic ground-motion models (GMMs). For the nonergodic path effects, our nonergodic GMM has two path effects terms: a cell-specific linear-distance scaling, denoted as δP2PQ, that mimics the effects of a 3D Q structure, and a source- and site-specific term related to the effects of the 3D velocity structure, denoted δP2PV. The resulting model provides spatially varying nonergodic source, site, and path terms. The GMM smoothly interpolates and extrapolates the nonergodic terms in space so that the GMM can be applied to any combination of source and site locations in the region. In regions where data are sparse, the median nonergodic terms will approach zero but with large epistemic uncertainty. Over the period range of 0–10 s, the ranges of the standard deviations of the nonergodic source, path, and site terms are 0.2–0.65, 0.25–0.40, and 0.15–0.40 natural log units, respectively. The fully nonergodic models lead to a 40%–55% decrease in the aleatory standard deviation compared to the reference ergodic GMMs. This large reduction in the aleatory standard deviation combined with the change in the median given by the nonergodic terms can exert a significant impact on the computed seismic hazard for the Japan region.

Understanding subduction infancy to mature subduction in Southwest Japan via the self-consistent formation of a weak slab interface

The weak slab interface controls long-term subduction dynamics. A weak hydrous layer at the slab interface promotes mechanical decoupling between the forearc mantle and the subducting slab and converts a hot forearc mantle to a cold mantle. Often referred to as a cold nose, the cold forearc mantle, plays a key role in the transition from subduction infancy to mature subduction. This study was the first to numerically demonstrate the self-consistent formation of a weak hydrous layer with permeability anisotropy based on the Southwest Japan subduction zone case, where transition-related geological features were present. Our models showed that mechanical decoupling by spontaneous downdip growth of the weak hydrous layer created a cold nose by converting a hot forearc mantle to a cold mantle. The emergence of a cold nose explained the migration of the forearc-to-arc volcanic front, expressed as the formation of mid-Miocene forearc high-magnesium andesite and Quaternary arc adakite. Furthermore, the weak hydrous layer providing a pathway for free-water transport toward the mantle wedge tip elucidates slab/mantle-derived geochemical components in deep groundwater as well as large S-wave delay times and non-volcanic seismic tremors in the forearc.

Seismic anisotropy and geodynamics of the East Japan subduction zone

Seismic anisotropy in the East Japan arc has been extensively investigated by conducting shear-wave splitting measurements, receiver-function analyses, and tomographic inversions of body-wave travel times and surface-wave dispersion data, which have provided a wealth of information on dynamic processes associated with active subduction of the Pacific plate. Measuring shear-wave splitting is popular and effective to detect seismic anisotropy, but it has poor depth resolution. This problem has been overcome by conducting 3-D anisotropic tomography, which has high-resolution in both lateral and vertical directions. Both P and S wave anisotropies are revealed in the crust, which are caused by alignment or preferred orientation of crustal minerals and stress-induced microcracks related to active faults. Trench-normal fast-velocity directions (FVDs) of azimuthal anisotropy are revealed in the back-arc mantle wedge, reflecting subduction-driven convection there. Trench-parallel FVDs appear in the forearc mantle wedge under the land area, which may reflect deformation that results in B-type olivine fabric. The forearc mantle wedge offshore may lack anisotropy, suggesting that it is stagnant and decoupled from the subducting slab and does not participate in the viscous flow, in sharp contrast with the rest of the mantle wedge. The most significant findings of the body-wave anisotropic tomography are its constraints on the slab anisotropy. The subducting Pacific slab exhibits mainly trench-parallel FVDs, which reflect shape-preferred orientation of crystals and cracks related to normal faults produced in the outer-rise area before the plate subduction, overprinting the fossil anisotropy that the Pacific plate gained when it was produced at the mid-ocean ridge. Trench-parallel intraslab fast velocity planes of anisotropy intersect the slab upper surface at high angles (∼45–90°), reflecting aligned hydrated faults in the slab. Ruptures of the hydrated faults in the upper part of the slab may cause large intraslab earthquakes (M ≥7.0) that take place frequently beneath the forearc area. Trench-normal FVDs also appear in the subslab mantle, which may reflect asthenospheric shear deformation associated with the overlying slab subduction.

Long-term slow slip events with and without tremor activation in the Bungo Channel and Hyuganada, southwest Japan

Slow slip events (SSEs) lasting for approximately 1 year occur every 6–8 years around the Bungo Channel in the southwest Japan subduction zone. The slip time evolution of the latest Bungo Channel SSE that occurred in 2018–2019 has been studied; however, the detailed spatial and temporal relationship between the slip process and other nearby phenomena, such as tectonic tremors and SSEs, is not well understood. Moreover, the migration of such long-term SSEs from Hyuganada to Shikoku through the Bungo Channel has been suggested, but a slip process connecting the SSEs has not been observed. In this study, we utilized 21 continuous global navigation satellite system (GNSS) stations around the Bungo Channel and Hyuganada that have been installed by us since 2014 in addition to GNSS Earth Observation System (GEONET) stations to improve the spatial resolution of such interplate slip. Based on these data, we estimate the spatial and temporal slip evolutions of the major SSE in 2018–2019, which was accompanied by tremor activity in the deep episodic tremor and slip (ETS) zone, and a smaller SSE in 2015–2016 without tremor activity. We show that the slip area of the major SSE overlaps the ETS zone, whereas that of the smaller SSE does not. This strongly suggests that synchronized tremor activity with an SSE requires a slip close to or overlapping the ETS zone. We also show two distinct slip propagation paths from the Oita area during the 2018–2019 sequence: one is a southward propagation to the Miyazaki area, leading to an SSE around the Miyazaki Plain, and the other is an eastward propagation to an area close to Cape Ashizuri, where “invading slip” is proposed to propagate from the ETS zone to a shallower megathrust source area. These slip propagations may be two of fundamental slip modes that connect slow-slip patch-like areas around the Bungo Channel and Hyuganada.Graphic

Anisotropic tomography of the East Japan subduction zone: influence of inversion algorithms

SUMMARY An important element of seismic tomography is the inversion process. In this work, we use P-wave arrival times of local earthquakes recorded at onshore and offshore seismic stations in East Japan to investigate the influence of two well-known inversion algorithms (LSQR and L-BFGS-B) on anisotropic tomography. Our synthetic tests show that a large damping parameter in the LSQR algorithm can lead to a stable and fast convergence, but it can result in many small value disturbances. The L-BFGS-B algorithm, which has second-order convergence, could converge fast to the optimal solution without damping regularization, but an inappropriate bound on the unknown parameters makes them hard to be recovered fully and causes strong trade-off between isotropic velocity and azimuthal anisotropy. If appropriate control parameters are adopted, the two inversion algorithms lead to almost the same results, though the L-BFGS-B provides a more efficient convergence and leads to a slightly better fit to the data than LSQR does. The two algorithms are applied to investigate the 3-D P-wave velocity (Vp) structure and azimuthal anisotropy of the East Japan subduction zone. Our results show that high-Vp anomalies and trench-normal fast-velocity directions (FVDs) exist in the forearc crust beneath the Pacific Ocean off South Hokkaido, which may reflect a cold and hydrated forearc crust with aligned microcracks or fractures. Significant low-Vp anomalies and trench-parallel FVDs exist at 40–80 km depths beneath Hokkaido, reflecting a water-rich mantle wedge with aligned B-type olivine. In the subducting Pacific slab, strong anisotropy with trench-parallel FVDs is revealed, reflecting localized horizontal bending of the slab.

Controls on melt focusing beneath old subduction zones: A case study of northeast Japan

Arc volcanoes in global subduction zones are geographically focused regardless of subduction parameters, and the Japan subduction zone is an excellent natural laboratory to examine arc volcanoes because young and old oceanic plates subduct along southwest and northeast Japan, respectively. Compared with the arc volcanoes in southwest Japan, which are formed by restricted sub-arc melting in the mantle wedge, those in northeast Japan should result from melt transport from the sub-backarc to the sub-arc mantle (i.e., melt focusing) because the liberated water from the sub-arc to the sub-backarc slab generates flux melting in the sub-arc to sub-backarc mantle. In this study, we quantitatively evaluated the partial melting and melt transport in the mantle wedge beneath northeast Japan using a series of two-dimensional numerical models with controlling factors, such as melt viscosity, melt density, mechanical decoupling depth at the slab interface, and melt freezing. Model calculations showed that a large melt viscosity and/or density allow corner-flow dominant melt behavior, which focuses the melt from the sub-backarc to the sub-arc mantle. A greater mechanical decoupling depth prevents the melt from being dragged toward the forearc mantle, enhancing the sub-arc melt focusing. Melt freezing prevents the melt from being dragged toward the back-arc mantle, further enhancing the sub-arc melt focusing. Calculated melt distribution was consistent with geographically focused arc volcanoes, seismic tomography, and geochemically estimated melt production in the mantle beneath northeast Japan; that is, the melt focusing in the old subduction zones is attributed to the corner flow and melt freezing in the mantle wedge.

Frictional Behavior Downdip Along the Subduction Megathrust: Insights From Laboratory Experiments on Exhumed Samples at In Situ Conditions

AbstractSubduction megathrusts exhibit a spectrum of slip modes, including catastrophic earthquakes. Although the mechanical and frictional properties of materials sampled from subduction zones have been studied extensively, few datasets have been collected for compositions and at pressure and temperature conditions representative of those in situ. The Nankai subduction zone in southwest Japan is a well‐studied margin, and abundant data provide an opportunity to advance our understanding of fault and earthquake physics. Here, we use samples exhumed in the Shimanto and Sanbagawa Belts on Shikoku Island of southwest Japan that represent analogs for materials along the present‐day megathrust at depths of ∼5–>25 km, and we shear these at their peak in situ pressure‐temperature (P‐T) conditions. Effective normal stresses range from 28 to 192 MPa, and temperatures from 105°C to 470°C. We used pore fluid pressures of 45–240 MPa, corresponding to fluid overpressure ratios λ of 0.65 and 0.90. Slip velocities of 0.1–100 μm/s were used, in order to focus on the nucleation of instability and earthquakes. We found predominantly velocity‐strengthening (inherently stable) behavior under all conditions for λ = 0.65. For λ = 0.90, velocity‐weakening behavior was observed at 350°C, with velocity‐strengthening behavior at lower and higher temperatures. The rate/state frictional stability parameter (a‐b) increases with slip velocity at temperatures up to ∼200°C and remains constant or decreases with slip velocity at higher temperatures. Overall, our results demonstrate the potentially important roles of both temperature and slip velocity in controlling the distribution of stress and frictional rheology along subduction thrusts.

Numerical estimation of a tsunami source at the flexural area of Kuril and Japan Trenches in the fifteenth to seventeenth century based on paleotsunami deposit distributions in northern Japan

Paleotsunami deposit investigations and numerical tsunami computations have been performed to elucidate the source and size of large tsunamis along the Kuril to Japan Trenches, particularly for unusual tsunamis that occurred in the seventeenth century, the 1611 CE Keicho tsunami (M 8.1) along the Japan Trench and seventeenth-century tsunami (> Mw 8.8) along the Kuril Trench, which caused serious damages on the coastal residents and environments. Moreover, several paleotsunami deposits dating from the thirteenth to eighteenth centuries have been reported along the area between the Kuril and Japan subduction zones, but their sources have not been clarified. In this study, we estimated the tsunami sources from numerical simulations using the distribution of fifteenth- to seventeenth-century tsunami deposits at Sekinehama along the coast of the Shimokita Peninsula. Based on numerical simulations with previously proposed fault models, the tsunami deposits showing similar ages at Sekinehama and another site on the coast of Shimokita Peninsula, which are within 50 km apart, could not be explained except with the huge earthquake models (> Mw 9.1), whose rupture zones extend to not only the Kuril or Japan Trenches but also their flexural area. Thus, we modified or newly proposed twelve fault models located in the flexural area between the two trenches to explain tsunami deposits possibly around the seventeenth century at the above-mentioned two sites on the coast of Shimokita Peninsula. Simulations using these models elucidated that the rupture in the shallow or deep plate boundaries with > 14–32 m slip (> Mw 8.55–8.76) is necessary. If the tsunami deposits around the seventeenth century along the Iburi–Hidaka coast in Hokkaido and those at the two sites mentioned above might be left by an identical event, an interplate earthquake with > 18–40 m slip (> Mw 8.62–9.2) in the flexural area is needed. Moreover, this interplate earthquake might have occurred in the deep plate boundary than in the shallower plate boundary based on slip deficit and slow earthquake distribution data. Our results offer significant insights into a large earthquake (> M 8) along the Kuril and Japan Trenches in the fifteenth to seventeenth century.

Elastic Slab in Viscoelastic Mantle: Effects on Determining Megathrust Slip and Mantle Viscosity During Postseismic and Interseismic Phases

AbstractEarthquake‐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.

Attenuation Characteristics of High-Frequency Ground Motions from Local Sources Caused by Great Subduction Zone Earthquakes in Northeast Japan

AbstractGround-motion attenuation characteristics are examined for the peak ground accelerations (PGAs) and peak ground velocities (PGVs) from strong-motion generation areas (SMGAs), which emit strong high-frequency waves during great subduction zone earthquakes. A conventional ground-motion prediction equation (GMPE) for earthquakes is designed based on the magnitude and distance from the source fault to predict the peak ground motion amplitude. For great subduction zone earthquakes, significant wavetrains of high-frequency ground motions are often separately observed in seismograms, and the corresponding rupture areas are estimated as SMGAs along the plate interface. In this case, although the advantages of using the shortest distances measured from the closest SMGAs rather than the shortest fault distance have been confirmed in previous studies, it is more physically reasonable to examine the ground-motion attenuation characteristics of individual SMGAs based on their magnitudes and locations. Therefore, we examine the attenuation characteristics of the PGAs and PGVs of individual SMGAs of the 2003 Mw 8.2 Tokachi–Oki earthquake and 2011 Mw 9.1 Tohoku earthquake in the Japan subduction zone considering data at stations within an SMGA distance of 300 km and the SMGA moment magnitude. The model, exhibiting the same functional form as a conventional GMPE, is fitted to the PGA and PGV data pertaining to each SMGA that are normalized to a bedrock site with VS30=760 m/s at any distance. The uncertainties in the obtained PGAs and PGVs are 3.6% and 13.5% lower, respectively, in terms of the residual in logarithmic units than those in the results of previous approaches considering only the SMGA distance. This result could help develop GMPEs for SMGAs to more appropriately predict the strong motions generated during great subduction zone earthquakes.

Role of fluids and seamount subduction in interplate coupling and the mechanism of the 2021 Mw 7.1 Fukushima-Oki earthquake, Japan

The frequent occurrence of large thrust earthquakes and subsequent tsunamis in a subduction zone are attributed to multiple factors, including structural heterogeneity, fluid saturation, and topographic variations of the subducting oceanic plate. To investigate these impacts on the generation of the 2021 Mw 7.1 Fukushima-Oki earthquake and interplate seismic coupling of the northeastern (NE) Japan subduction zone, high resolution three-dimensional (3-D) seismic structures of velocities (Vp, Vs) and Poisson's ratio (σ) were determined through a multi-parameter joint inversion method using 328,625 travel times of P- and S-wave receiver-source pairs from the earthquakes recorded by onshore seismic stations and ocean bottom seismometers. We found that more than 92.5% of M 6+ thrust earthquakes occurred in and around high-V zones, while less than 7.5% of the quakes were located in low-V patches. The high-V zones with high density M ≥ 6.0 thrust earthquakes (including the 2011 Mw 9.0 Tohoku mainshock) were thought to indicate strong interplate seismic coupling (asperities), while the low-V patches lacking large thrust-earthquakes reflected weak interplate coupling or decoupling. The present results revealed that seamount (including the Joban seamount chain and the bulgy Kitamachi volcanic batholith) subduction had a significant influence on seismography and the mechanical strength variation of interplate coupling in the forearc seismogenic zone. It was demonstrated that structural heterogeneities associated with fluids and seamount subduction played a fundamental role in interplate seismic coupling and large thrust-earthquake generation along the slab upper boundary in the NE Japan subduction zone.

3-D synthetic modelling and observations of anisotropy effects on SS precursors: implications for mantle deformation in the transition zone

SUMMARY The Earth's mantle transition zone (MTZ) plays a key role in the thermal and compositional interactions between the upper and lower mantle. Seismic anisotropy provides useful information about mantle deformation and dynamics across the MTZ. However, seismic anisotropy in the MTZ is difficult to constrain from surface wave or shear wave splitting measurements. Here, we investigate the sensitivity to anisotropy of a body wave method, SS precursors, through 3-D synthetic modelling and apply it to real data. Our study shows that the SS precursors can distinguish the anisotropy originating from three depths: shallow upper mantle (80–220 km), deep upper mantle above 410 km, and MTZ (410–660 km). Synthetic resolution tests indicate that SS precursors can resolve $\ge $3 per cent azimuthal anisotropy where data have an average signal-to-noise ratio (SNR = 7) and sufficient azimuthal coverage. To investigate regional sensitivity, we apply the stacking and inversion methods to two densely sampled areas: the Japan subduction zone and a central Pacific region around the Hawaiian hotspot. We find evidence for significant VS anisotropy (15.3 ± 9.2 per cent) with a trench-perpendicular fast direction (93° ± 5°) in the MTZ near the Japan subduction zone. We attribute the azimuthal anisotropy to the grain-scale shape-preferred orientation of basaltic materials induced by the shear deformation within the subducting slab beneath NE China. In the central Pacific study region, there is a non-detection of MTZ anisotropy, although modelling suggests the data coverage should allow us to resolve at least 3 per cent anisotropy. Therefore, the Hawaiian mantle plume has not produced detectable azimuthal anisotropy in the MTZ.

Small‐Scale Intraslab Heterogeneity Weakens Into the Mantle Transition Zone

AbstractSmall‐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.

Receiver Function Imaging of the Amphibious NE Japan Subduction Zone—Effects of Low‐Velocity Sediment Layer

AbstractThis study presents reflectivity images of the northeast (NE) Japan subduction zone continuous across the ocean and land. As nearly half of its forearc region is under the ocean, data from ocean bottom seismometers (OBSs) must be utilized to fully image the region by passive seismic analysis. The use of OBS data has been a challenge due to inherent characters of the ocean bottom observations: high noise level and effects of seafloor sediment. Now, decent imaging is possible in NE Japan overcoming the high level noise due to the accumulated data set of the OBSs. The low‐velocity of seafloor sediment significantly delays and amplifies S waves passing through them, and thus complicates teleseismic waveforms. We identify and correct these effects to produce coherent receiver function images throughout the amphibious subduction zone. Our images provide a potential for discussing new structural features and will help better understanding of the dynamics of the NE Japan subduction zone.

Attenuation contrast in mantle wedge across the volcanic front of northeastern Japan that controls propagations of high-frequency S-wave later phases

Distinct later phases of waves with rich high-frequency (> 8 Hz) components were observed for intraslab earthquakes that occurred at intermediate depths, particularly at depths exceeding 100 km, in the northeastern (NE) Japan subduction zone. These high-frequency later phases (HFLPs) showed anomalously large peak-amplitude delays, up to ~ 50 s after direct S-wave arrivals at stations in the backarc region. Using a source-scanning algorithm, we investigated the locations of passing points affecting the propagation of HFLPs. The passing points were estimated to be in the forearc region in the entire NE Japan, indicating that HFLPs are scattered waves that pass through the forearc region. The propagating HFLPs seem to bypass the backarc mantle wedge, as a consequence of the distinct attenuation contrast in the mantle wedge across the volcanic front in NE Japan. These HFLP observations suggest that the high-attenuation zone in the backarc mantle wedge controls propagations of the high-frequency waves of intraslab earthquakes, in addition to the scatterers possibly locate in the forearc region.

New Insights Into the Slip Budget at Nankai: An Iterative Approach to Estimate Coseismic Slip and Afterslip

AbstractThe slip budget at a subduction zone megathrust is met by a spectrum of slip behaviors including coseismic slip from great earthquakes, afterslip, episodic tremor and slip (ETS), long‐term slow slip events (SSEs), and quasisteady interseismic creep. Long‐term geodetic monitoring at the Nankai subduction zone in southwest Japan has recorded deformation observations spanning nearly a complete earthquake cycle, making this region the best study area for obtaining a full understanding of how the slip budget is partitioned along the interface and how the rheological conditions evolve through space and time. In this work, we conduct an iterative inversion of vertical surface displacement data to jointly estimate coseismic slip during the 1944/46 Tonankai and Nankai earthquakes, spatiotemporal afterslip history following those earthquakes, and postseismic mantle flow. A suite of postseismic models with variable elastic slab and plate geometry and mantle viscosity all show at least 30 years of sustained afterslip within the ETS zone and the gap between the top of the ETS zone and the fully coupled zone. This result demonstrates that stable, long‐duration afterslip is possible in the ETS zone. Summing the estimated cumulative slip from this study with estimates from previous studies of cumulative slip during the interseismic period, including episodic SSEs, we find that eastern Shikoku island has nearly met all of its slip budget for a 150‐year interseismic cycle, whereas western Shikoku island has a slip deficit of nearly half the budget (4 out of 8.25 meters) above 20 km depth on the subduction interface.