Shallow very-low-frequency earthquakes accompany slow slip events in the Nankai subduction zone

Link between the Nankai underthrust turbidites and shallow slow earthquakes

Trench sediments such as pelagic clay or terrigenous turbidites have long been invoked to explain the seismogenic behavior of the megathrust fault (i.e., décollement). Recent numerous studies suggest that slow earthquakes may be associated with huge megathrust earthquake; however, controls on the slow earthquake occurrence remain poorly understood. We investigate seismic reflection data along the Nankai Trough subduction zone to understand the correlations between the spatial distribution of the broad turbidites and along-strike variations in shallow slow earthquakes and slip-deficit rates. This report presents a unique map of regional distribution of the three discrete Miocene turbidites that underthrust apparently along the décollement beneath the Nankai accretionary prism. A comparison of distributions of the Nankai underthrust turbidites, shallow slow earthquakes, and slip-deficit rates enables us to infer that the underthrust turbidites may cause primarily low pore-fluid overpressures and high effective vertical stresses across the décollement, leading to potentially inhibiting the slow earthquake occurrence. Our findings provide a new insight into potential role of the underthrust turbidites for shallow slow earthquakes at subduction zone.

Décollement geometry controls on shallow very low frequency earthquakes

Recent studies have documented the occurrence of shallow very low frequency earthquakes (VLFE) in subduction zones. The heterogeneity of the materials or stresses that act on the plate interface results in the variable slip rate. Stress on the décollement can be controlled by the décollement geometry and the regional stress, which is also able to control the material properties. We determined the distribution of stress along the shallow portion of the décollement in the Nankai Trough using a three-dimensional (3D) seismic survey and regional stress analysis to construct maps of normalized slip tendency (T_s′) and dilation tendency (T_d). Alignments of VLFEs trend parallel to the trends of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{s}^{^{\prime}}$$\end{document}Ts′ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{d}$$\end{document}Td. On the other hand, very low \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{s}^{^{\prime}}$$\end{document}Ts′ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{d}$$\end{document}Td areas probably act as barriers that limit the number of VLFEs that can migrate towards the trench. Because the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{s}^{^{\prime}}$$\end{document}Ts′ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{d}$$\end{document}Td distributions are derived only from the décollement geometry and the regional stress without incorporating any data on sediment properties, the consistency between the trends suggests that the décollement geometry is the primary control on VLFE activity.

Shallow slow earthquakes to decipher future catastrophic earthquakes in the Guerrero seismic gap

The Guerrero seismic gap is presumed to be a major source of seismic and tsunami hazard along the Mexican subduction zone. Until recently, there were limited observations at the shallow portion of the plate interface offshore Guerrero, so we deployed instruments there to better characterize the extent of the seismogenic zone. Here we report the discovery of episodic shallow tremors and potential slow slip events in Guerrero offshore. Their distribution, together with that of repeating earthquakes, seismicity, residual gravity and bathymetry, suggest that a portion of the shallow plate interface in the gap undergoes stable slip. This mechanical condition may not only explain the long return period of large earthquakes inside the gap, but also reveals why the rupture from past M < 8 earthquakes on adjacent megathrust segments did not propagate into the gap to result in much larger events. However, dynamic rupture effects could drive one of these nearby earthquakes to break through the entire Guerrero seismic gap.

Tracking Deep Sediment Underplating in a Fossil Subduction Margin: Implications for Interface Rheology and Mass and Volatile Recycling

The architecture and mechanical properties of the subduction interface impact large-scale subduction processes, including mass and volatile recycling, upper-plate orogenesis, and seismic behavior. The nature of the deep subduction interface, where a dominantly frictional megathrust likely transitions to a distributed ductile shear zone, is poorly understood, due to a lack of constraints on rock types, strain distribution, and interface thickness in this depth range. We characterized these factors in the Condrey Mountain Schist, a Late Jurassic to Early Cretaceous subduction complex in northern California that consists of an upper and lower unit. The Lower Condrey unit is predominantly pelagic and hemipelagic metasediment with m-to km-scale metamafic and metaserpentinitic ultramafic lenses all deformed at epidote blueschist facies (0.7-1.1 GPa, 450°C). Major and trace element geochemistry suggest tectonic erosion of the overriding plate sourced all ultramafic and some mafic lenses. We identified two major ductile thrust zones responsible for Lower Condrey unit assembly, with earlier strain distributed across the structural thickness between the ductile thrusts. The Lower Condrey unit records distributed deformation across a sediment-dominated, 2+ km thick shear zone, possibly consistent with low velocity zones observed in modern subduction zones, despite subducting along a sediment poor, tectonically erosive margin. Periodic strain localization occurred when rheological heterogeneities (i.e., km-scale ultramafic lenses) entered the interface, facilitating underplating that preserved 10%-60% of the incoming sediment. Modern mass and volatile budgets do not account for erosive margin underplating, so improved quantification is crucial for predicting mass and volatile net flux to Earth's interior.

Shallow slow slip events along the Nankai Trough detected by GNSS-A

Various slow earthquakes (SEQs), including tremors, very low frequency events, and slow slip events (SSEs), occur along megathrust zones. In a shallow plate boundary region, although many SEQs have been observed along pan-Pacific subduction zones, SSEs with a duration on the order of a year or with a large slip have not yet been detected due to difficulty in offshore observation. We try to statistically detect transient seafloor crustal deformations from seafloor geodetic data obtained by the Global Navigation Satellite System-Acoustic (GNSS-A) combination technique, which enables monitoring the seafloor absolute position. Here, we report the first detection of signals probably caused by shallow large SSEs along the Nankai Trough and indicate the timings and approximate locations of probable SSEs. The results show the existence of large SSEs around the shallow side of strong coupling regions and indicate the spatiotemporal relationship with other SEQ activities expected in past studies.

Characteristic activities of slow earthquakes in Japan

Slow earthquakes are a recently discovered phenomenon that mainly occur updip and downdip of the seismogenic zones of great earthquakes along the subducting plate interface. The spatiotemporal activity of various slow earthquakes occurring in the Nankai subduction zone is characterized by along-strike heterogeneity and along-dip systematic changes. Various slow earthquakes are horizontally distributed at their own depths and along-strike segments can be observed with respect to this distribution downdip of the locked zone; however, slow and great earthquakes occur in the same depth range near the Nankai Trough and Japan Trench axes. The frequently observed spatiotemporal interactions between different slow earthquakes can be attributed to their sensitivity and the stress transfer of the surrounding areas. This stress transfer is expected to extend to the adjacent sections in the seismogenic zone. Therefore, precise monitoring of slow earthquakes is important for future evaluations of great earthquakes, which requires the long-term maintenance and continuous improvement of the high-quality observation networks.