Slow Earthquakes Research Articles

Seismicity distribution in the Tonankai and Nankai seismogenic zones and its spatiotemporal relationship with interplate coupling and slow earthquakes

We conducted seismic tomography to estimate the seismic velocity structure and to evaluate the spatiotemporal distribution of interplate earthquakes of the Kii Peninsula, central Honshu, Japan, where the Tonankai and Nankai megathrusts are located. Microearthquakes were quantitatively detected by using the data from a cable-type seafloor seismic observation network, completed in 2015. Our velocity model was consistent with the previous 2-D active-source surveys, which reported the areal extent of key structures: a high-velocity zone beneath Cape Shionomisaki, a subducted seamount off Cape Muroto, and the subducted Paleo-Zenith Ridge. The absence of any other subducted seamount with the same or larger spatial scale, than the identified key structures, was confirmed. Our velocity model also revealed that there was not a simple relationship between areas of large coseismic slip or strong interplate coupling and areas of high velocity in the overriding plate. Relocated hypocenters widely ranged from the upper plate to within the slab, while the most active region was attributed to the oceanic crust in the aftershock region of 2004 off-Kii earthquake. Compared with the results from the land-based observation network, the accuracy of the focal depth estimation was substantially improved. Furthermore, we identified the seismic activity in the vicinity of the plate boundary and determined 14 locations for interplate seismicity areas. They were primarily distributed in the range of seismogenic zone temperature (150–350 °C) along the plate boundary and were located outside of the strong interplate coupling zone. Several active areas of interplate earthquakes exhibited clustered activity during the periods of slow-slip events, observed and accompanied with shallow very-low-frequency earthquakes. Thus, regular interplate microearthquakes became active at the plate boundary in the conjunction with slow slip. In summary, as regular earthquakes provide a more accurate source location than slow earthquakes and can detect events of smaller magnitude, monitoring such interplate earthquakes may reveal spatiotemporal variations in the stick–slip conditions on the plate boundary.

Depth profile of frictional properties in the inner Nankai accretionary prism using cuttings from IODP Site C0002

We conduct frictional experiments using cuttings collected at Nankai Trough IODP Site C0002 over 980.5–3262.5 mbsf (meters below seafloor) depth interval to better understand the frictional properties through the accretionary prism. The experiments are conducted at the in situ effective normal stresses (9–37 MPa) under brine-saturated conditions, and the slip velocity is abruptly changed in a stepwise manner to either of 0.3, 3, or 33 µm/s after the steady-state friction is reached. The friction coefficient (μ) of the cuttings samples ranges from 0.45 to 0.60, with a slight increase in μ with increasing depth, related to decreasing smectite content. The velocity dependence of friction (a − b) is positive at all depths and ranges from 0.001 to 0.006, which indicates a velocity-strengthening behavior; these values are consistent with relatively homogeneous deformation microstructures. The critical slip distance (Dc) ranges from 0.5 to 123 μm, with relatively large values obtained for the smectite-rich samples. The changes in both the friction coefficient and rate- and state-friction parameters are likely associated with mineralogical change and consolidation with increasing depth. Although all of the cuttings samples collected from Site C0002 exhibit a velocity-strengthening behavior, a slight decreasing trend in a − b with increasing depth indicates either a nearly neutral velocity dependence or a possible transition to velocity-weakening behavior at greater depths, which may be attributed to the occurrence of slow earthquakes in the Nankai accretionary prism.

The High-Frequency Signature of Slow and Fast Laboratory Earthquakes.

Tectonic faults fail through a spectrum of slip modes, ranging from slow aseismic creep to rapid slip during earthquakes. Understanding the seismic radiation emitted during these slip modes is key for advancing earthquake science and earthquake hazard assessment. In this work, we use laboratory friction experiments instrumented with ultrasonic sensors to document the seismic radiation properties of slow and fast laboratory earthquakes. Stick‐slip experiments were conducted at a constant loading rate of 8 μm/s and the normal stress was systematically increased from 7 to 15 MPa. We produced a full spectrum of slip modes by modulating the loading stiffness in tandem with the fault zone normal stress. Acoustic emission data were recorded continuously at 5 MHz. We demonstrate that the full continuum of slip modes radiate measurable high‐frequency energy between 100 and 500 kHz, including the slowest events that have peak fault slip rates <100 μm/s. The peak amplitude of the high‐frequency time‐domain signals scales systematically with fault slip velocity. Stable sliding experiments further support the connection between fault slip rate and high‐frequency radiation. Experiments demonstrate that the origin of the high‐frequency energy is fundamentally linked to changes in fault slip rate, shear strain, and breaking of contact junctions within the fault gouge. Our results suggest that having measurements close to the fault zone may be key for documenting seismic radiation properties and fully understanding the connection between different slip modes.

Source Characteristics and Along‐Strike Variations of Shallow Very Low Frequency Earthquake Swarms on the Nankai Trough Shallow Plate Boundary

AbstractWe detected shallow very low frequency earthquakes (VLFEs) off the Cape Muroto and Kii Channel in the Nankai subduction zone and estimated their moment rate functions. Combining the new and previously estimated catalogs, we obtained the comprehensive catalog of shallow VLFE moment rate functions along the Nankai Trough. We defined the shallow VLFE swarms and investigated the scaling relationships of their cumulative moments, activity area, and durations in each region. Detected swarms were considered candidates for shallow slow slip events. A similar scaling relationship was observed between the cumulative moments and activity areas, irrespective of regions. It indicates similar stress drops in each region. However, the relationship between the cumulative moments and durations varied. This difference was explained by the along‐strike variations in the faulting conditions of shallow slow earthquakes, such as material or hydrological properties.

Frictional stability and hydromechanical coupling of serpentinite-bearing fault gouge

SUMMARY Observations of slow earthquakes and tremor have raised fundamental questions about the physics of quasi-dynamic rupture and the underlying fault zone processes. The presence of serpentinite at P-T conditions characteristic of deep tremor and slow earthquakes suggests that it plays an important role in controlling complex fault slip behaviour. Here, we report on experiments designed to investigate the frictional behaviour of serpentinite sampled from outcrop exposures (SO1 and SO2) of altered ultramafic rocks present at depth, and recovered from the SAFOD borehole (G27). XRD analyses reveal the presence of chrisotyle, lizardite, kaolinite, talc in SO1; lizardite, clinochlore and magnetite in SO2; and lizardite, quartz and calcite in G27. We sheared fault gouge in a double-direct shear configuration using a true triaxial deformation apparatus. The effective normal stress was varied from 2 to 40 MPa. We conducted velocity stepping tests and slide-hold-slide (SHS) tests in each experiment to characterize frictional stability and healing. At the end of each experiment, post-shear permeability was measured and the samples were recovered for microstructural analysis. The steady-state friction coefficient was μ = 0.17 for SO1, μ = 0.33 for SO2 and μ = 0.53 for G27. Overall, the gouges exhibit velocity strengthening behaviour, and become nearly velocity neutral at 40 MPa effective normal stress. SHS tests show positive healing rates for SO2 and G27, whereas SO1 exhibits zero or negative healing rates. Permeability decreases with increasing σn’, with SO1 (k = 10–20 m2) showing the lowest values. Microstructural observations reveal a well-developed R-Y-P fabric in SO1, which is not observed in SO2 and G27. We posit that the development of shear fabric controlled by mineralogy governs frictional and hydrological properties. In this context, when serpentinite is associated with other weak phyllosilicate minerals, frictional stability and hydrological properties can vary greatly, with a potential control on the mode of fault failure.

Structural Anomaly at the Boundary Between Strong and Weak Plate Coupling in the Central‐Western Nankai Trough

AbstractTo investigate structural characteristics associated with variations in slip behavior on the plate boundary in the Nankai subduction zone, we conducted seismic reflection surveys in the central‐western Nankai Trough. Data were processed by pre‐stack depth migration. The resulting depth map of the top of the subducting Philippine Sea plate shows detailed topographic features of the plate boundary, including subducting seamounts and large‐scale undulations of the plate surface. We report the presence of a structurally anomalous region, potentially with low velocity, in the overriding plate at the boundary between the zone of large coseismic slip during the 1946 Nankai earthquake and the area producing slow earthquakes. This anomaly appears to be related to depth‐dependent variations in plate boundary slip style and plate coupling in the central‐western Nankai Trough.

Long‐Term Weakening Processes and Short‐Term Seismic Slip Behavior of an Intraplate Mature Fault Zone: A Case Study of the Yangsan Fault, SE Korea

AbstractAn intraplate tectonic fault evolves with repeated seismic cycles that may exhibit multi‐mode slip behaviors (continuous aseismic sliding, slow slip with silent earthquakes, and earthquakes resulting from stick‐slip frictional instability). Here, we show long‐term (geological time scale) fault zone weakening processes and short‐term (seismic cycle scale) fault behavior of the Yangsan Fault, SE Korea, based on structural and mineralogical observations of an across‐strike fault zone sequence retrieved by inclined borehole drilling. Since the Late Cretaceous, physical and chemical weakening processes have formed a wide and complex fault zone characterized by interconnected clay‐rich fault core strands and fractured lenses. The deformation is more concentrated in the eastern block (Cretaceous sedimentary rock with felsic dikes) of the fault zone than in the western block (Paleogene granite) of the fault zone. The asymmetrical deformation pattern is controlled by the distinct properties of the different host rocks and differential fluid‐assisted reaction softening involving hydrothermal processes. Of the five main strands and numerous subsidiary strands in the wide fault zone, strand S1 between the sedimentary rock and granite has accommodated the largest displacements. In particular, the smectite‐rich, &lt;2‐cm‐wide principal slip zone within S1 has acted as the major pathway for past seismic slip. Our findings consider that slow deformation in the sedimentary rocks formed the asymmetrical wide fault zone over long geological time, whereas that the preferential propagation of repeated seismic slip events propagated preferentially along the narrow principal slip zone of S1.

Mechanical behaviors of intact antigorite as functions of temperature: Faulting, slow stick-slip and stable sliding

Triaxial compression experiments on intact antigorite were performed under undrained conditions with a confining pressure of 100–400 MPa, temperature of 25–700 °C and strain rate of ∼1.5 × 10−6−10−5 s−1. We report a sequential transition in mechanical behavior from faulting to slow stick-slip behavior and then to stable sliding with increasing temperature. Syndeformational dehydration of antigorite at 100 MPa and 550–650 °C produces slow stick-slip, and the slip velocity and shear stress drop generally decrease with increasing temperature. At 100 MPa and 650 °C, the slip velocity and shear stress drop of slow stick-slip gradually decrease, implying a trend toward stable sliding, which might be related to a continuous increase in fluid pressure induced by the accumulation of fluid over time in the process of antigorite dehydration. Our results indicate that dehydration makes the sample weaker but more stable and does not cause embrittlement; therefore, an intermediate-depth earthquake in a subduction zone may be triggered by brittle deformation in strong peridotite surrounding dehydrating antigorite. Similar to the antigorite fault gouge previously studied, intact antigorite can also induce slow stick-slip under partial antigorite dehydration, indicating that slow earthquakes could be triggered in the mantle wedge of a hot subduction zone.

Cover Image

The cover image is based on the Research Article Thermal regime and slab dehydration beneath the Izu-Bonin arc: Implications for fast and slow subduction earthquakes by Yingfeng Ji et al., https://doi.org/10.1111/ter.12568.

Stick-slip behavior of a clayey crustal fault

Recent studies have revealed that slip on clay-rich fault planes, which develop in the shallow crust and can cause huge earthquakes, is controlled by their flow properties rather than their frictional properties. Here, we simulate a seismogenic crustal fault and show that such a fault plane may behave as a thixotropic yield stress fluid. We then conduct numerical experiments using a simple spring-slider model that incorporates this thixotropic property, and we successfully reproduce spontaneous stick-slip behaviors that correspond to a variety of seismogenic processes, ranging from regular earthquakes with high slip rates to slow earthquakes with lower slip rates. This finding suggests that the seismic activity on a shallow clay-rich crustal fault may be governed by the time-evolving microstructure of the clay-water system that exists along the fault plane.

Very Low Frequency Earthquakes in Between the Seismogenic and Tremor Zones in Cascadia?

AbstractMegathrust earthquakes and their associated tsunamis cause some of the worst natural disasters. In addition to earthquakes, a wide range of slip behaviors are present at subduction zones, including slow earthquakes that span multiple orders of spatial and temporal scales. Understanding these events may shed light on the stress or strength conditions of the megathrust fault. Out of all types of slow earthquakes, very low frequency earthquakes (VLFEs) are most enigmatic because they are difficult to detect reliably, and the physical nature of VLFEs are poorly understood. Here we show three VLFEs in Cascadia that were dynamically triggered by a 2009 Mw 6.9 Canal de Ballenas earthquake in the Gulf of California. The VLFEs likely locate in between the seismogenic zone and the Cascadia episodic tremor and slip (ETS) zone, including one event with a moment magnitude of 5.7. This is the largest VLFE reported to date, causing clear geodetic signals. Our results show that the Cascadia megathrust fault might slip rapidly at some spots in this gap zone, and such a permissible slip behavior has direct seismic hazard implications for coastal communities and perhaps further inland. Further, the observed seismic sources may represent a new class of slip events, whose characteristics do not fit current understandings of slow or regular earthquakes.

Fluid Migration Before and During Slow Earthquakes in the Shallow Nankai Subduction Zone

AbstractFluid migration in subduction zones is a key controlling factor of slow and megathrust earthquakes at plate boundaries. During the migration, seismic velocity and heterogeneous structures in its pathways may be temporarily perturbed, preferably triggering slow earthquakes. Here, we show that transient changes of seismic heterogeneity occurred 0–9 months before shallow slow earthquakes in the Nankai subduction zone, Japan, using very long‐term (6–10 years) records of ambient seafloor noise. The heterogeneity changes preceding to shallow slow earthquakes were observed near the margin of the source region, while concurrent changes primarily occurred in the source region. We propose that the heterogeneity changes are attributed to dynamic fluid migration, and the difference in timings reflects the pore pressure level in the corresponding source region. When fluids are supplied to a source region under relatively low pressure, fluids are leaked out from its downdip or updip side, and slow earthquakes occur not immediately but with a time delay of at most 9 months. In the high pore pressure case, slow earthquakes occur immediately with fluid migration from the source region. This study suggests that the heterogeneous seismic structure is possibly changed by fluid migration before slow earthquakes in the Nankai subduction zone.

Spatiotemporal Variations of Shallow Very Low Frequency Earthquake Activity Southeast Off the Kii Peninsula, Along the Nankai Trough, Japan

AbstractCross‐correlation analysis was applied to long‐term onshore broadband records from April 2004 to March 2021 to detect and relocate shallow very low frequency earthquakes (VLFEs) southeast off the Kii Peninsula, along the Nankai Trough, Japan. We then determined the moment rate functions of detected shallow VLFEs using the Monte Carlo‐based simulated annealing method. According to this new comprehensive catalog, shallow VLFEs are widespread beneath the accretionary prism toe, but shallow VLFEs with large cumulative moments are localized around the western edge of the paleo‐Zenisu Ridge, which is subducted beneath southeast off the Kii Peninsula. Our results from the long‐term shallow VLFE catalog are well consistent with previous studies in this region, suggesting that heterogeneous structures and stress conditions due to the subducted paleo‐Zenisu Ridge promote the occurrence of shallow slow earthquakes. The relocated shallow VLFE epicenters illustrated three major episodes characterized by a similar activity area and five minor episodes characterized by different areas. The three major episodes exhibited slow frontal migration with different initiation locations, directions, and speeds, as well as several rapid reverse migrations. Episodes of minor activity were distributed in different locations within part of the area of major activity. Different patterns of shallow VLFE migration could reflect temporal changes in the pore‐fluid distribution or stress conditions of the plate boundary.

Continuous Tremor Activity With Stable Polarization Direction Following the 2014 Large Slow Slip Event in the Hikurangi Subduction Margin Offshore New Zealand

AbstractMany types of slow earthquakes have been discovered at subduction zones around the world. However, the physical process of these slow earthquakes is not well understood. To monitor offshore slow earthquakes, a marine seismic and geodetic experiment was conducted at the Hikurangi subduction margin from May 2014 to June 2015. During this experiment, a large slow slip event (Mw 6.8) occurred directly beneath the ocean bottom seismometer (OBS) network. In this study, S‐wave splitting and polarization analysis methods, which have been previously used on onshore data to investigate tremor and anisotropy, are applied to continuous OBS waveform data to identify tremors that are too small to detect by the envelope cross correlation method. Continuous tremor activity with stable polarization directions is detected at the end of the 2014 slow slip event and continued for about 2 weeks. The tremors are generated around a southwest bend in the slow slip contours and at the landward edge of a subducted seamount. Our findings corroborate a previous interpretation, based on burst‐type repeating earthquakes and intermittent tremor, that localized slow slip and tremor around the seamount was triggered by fluid migration following the large plate boundary slow slip event and indicate tremor occurred continuously rather than as isolated and sporadic individual events.

Model for tectonic tremors: Enduring events, moment rate spectrum, and moment-duration scaling.

Numerous attempts have been made to obtain earthquake statistics from a theoretical-physics perspective, but these studies mostly involve regular earthquakes. In recent years, a new category of earthquakes, referred to as slow earthquakes, has been discovered. Slow earthquakes emit only weak or no seismic signals and have different statistics than regular earthquakes. Here we propose a physical model for the tremor, which is a type of slow earthquake, introducing two competing timescales in a cellular automaton model. The proposed model reproduces some observation results for tremors, such as enduring events, moment-duration scaling, size distribution, and the power spectrum of the moment rate function.

Cryptoexplosive breccia: A genetic mark of the paleoearthquakes

全球到处发育的隐爆角砾岩表明隐爆作用的普遍性，并标志着古地震的一种形成机制。隐爆角砾岩具有特殊的角砾形状和构造。地震地质资料表明大地震产生的地表破裂带没有肉眼可观测到的断层位移。材料聚集应变能估算值表明岩石聚集的应变能不足以产生大地震。携带高能量的地核-地幔流体不仅为自己逸散提供了能量，而且为产生地震提供了充足的能量。深部流体以运移-聚集-爆炸-运移-聚集-爆炸的循环方式向地表逃逸，在地球内部不同深度产生了不同震级的地震。因此，超高压流体隐爆应该是地震形成的主要机制。地震隐爆成因机制不仅能够更好地解释各种天然地震（中深源地震、震群、慢地震和非双力偶性地震等）以及全球地震、火山和地热带的空间吻合，而且能够解释隐爆角砾岩及其伴生矿的形成。;The cryptoexplosive breccias found all over the world indicated that cryptoexplosion was pervasive in the geological history. Cryptoexplosive breccias are characterized by the special structure, which indicates a genetic mark of paleoearthquakes. The seismic-geological data indicated that the great earthquakes produced the ruptures with evident surface displacement, but on visible fault displacement. Estimation of strain energy accumulated in the materials demonstrates that rocks in the crust and mantle cannot accumulate huge strain energy to produce great earthquake. High energy density fluids derived from the core and mantle can provide energy for their migration, as well as enough energy for producing great earthquakes. The deep fluids migrate upward in the way of migration-accumulation-explosion-migration-accumulation-explosion, producing earthquakes with the different magnitudes at different depths. Therefore, the crypto-explosion can be considered as the predominant mechanism of earthquake generation. The proposed genetic mechanism not only can be used to much better explain the generation of different kinds of earthquakes, such as the moderate and deep focus earthquakes, earthquake swarm, slow earthquakes and no-dipole-couple earthquakes, as well as special overlap of the seismic, volcanic and geothermal zones in the world, but also to explain the formation of ore deposits companied with cryptoexplosive breccias.

Tremor Waveform Extraction and Automatic Location With Neural Network Interpretation

Active faults release tectonic stress imposed by plate motion through a spectrum of slip modes, from slow, aseismic slip, to dynamic, seismic events. Slow earthquakes are often associated with tectonic tremor, nonimpulsive signals that can easily be buried in seismic noise and go undetected. We present a new methodology aimed at improving the detection and location of tremors hidden within seismic noise. After identifying tremors with a classic convolutional neural network (CNN), we rely on neural network attribution to extract core tremor signatures. We observe that the signals resulting from the neural network attribution analysis correspond to a waveform traveling in the Earth&#x2019;s crust and mantle at wavespeeds consistent with seismological estimates. We then use these waveforms signatures to locate the source of tremors with standard array-based techniques. We apply this method to the Cascadia subduction zone, where we identify tremor patches consistent with existing catalogs. This approach allows us to extract small signals hidden within the noise, and to locate more tremors than in existing catalogs.

Slow earthquakes rumble through the Pacific Northwest

Slow earthquakes rumble through the Pacific Northwest

Segmentation of Shallow Slow Slip Events at the Hikurangi Subduction Zone Explained by Along‐Strike Changes in Fault Geometry and Plate Convergence Rates

Over the last two decades, geodetic and seismic observations have revealed a spectrum of slow earthquakes along the Hikurangi subduction zone in New Zealand. Of those, shallow slow slip events (SSEs) that occur at depths of less than 15 km along the plate interface show a strong along‐strike segmentation in their recurrence intervals, which vary from ∼1 yr from offshore Tolaga Bay in the northeast to ∼5 yr offshore Cape Turnagain ∼300 km to the southwest. To understand the factors that control this segmentation, we conduct numerical simulations of SSEs incorporating laboratory‐derived rate‐and‐state friction laws with both planar and non‐planar fault geometries. We find that a relatively simple model assuming a realistic non‐planar fault geometry reproduces the characteristics of shallow SSEs as constrained by geodetic observations. Our preferred model captures the magnitudes and durations of SSEs, as well as the northward decrease of their recurrence intervals. Our results indicate that the segmentation of SSE recurrence intervals is favored by along‐strike changes in both the plate convergence rate and the downdip width of the SSE source region. Modeled SSEs with longer recurrence intervals concentrate in the southern part of the fault (offshore Cape Turnagain), where the plate convergence rate is lowest and the source region of SSEs is widest due to the shallower slab dip angle. Notably, the observed segmentation of shallow SSEs cannot be reproduced with a simple planar fault model, which indicates that a realistic plate interface is an important factor to account for in modeling SSEs.

Automatic Detection of Slow Slip Events Using the PICCA: Application to Chilean GNSS Data

The detection of transient events related to slow earthquakes in GNSS positional time series is key to understanding seismogenic processes in subduction zones. Here, we present a novel Principal and Independent Components Correlation Analysis (PICCA) method that allows for the temporal and spatial detection of transient signals. The PICCA is based on an optimal combination of the principal (PCA) and independent component analysis (ICA) of positional time series of a GNSS network. We assume that the transient signal is mostly contained in one of the principal or independent components. To detect the transient, we applied a method where correlations between sliding windows of each PCA/ICA component and each time series are calculated, obtaining the stations affected by the slow slip event and the onset time from the resulting correlation peaks. We first tested and calibrated the method using synthetic signals from slow earthquakes of different magnitudes and durations and modelled their effect in the network of GNSS stations in Chile. Then, we analyzed three transient events related to slow earthquakes recorded in Chile, in the areas of Iquique, Copiapó, and Valparaíso. For synthetic data, a 150 days event was detected using the PCA-based method, while a 3 days event was detected using the ICA-based method. For the real data, a long-term transient was detected by PCA, while a 16 days transient was detected by ICA. It is concluded that simultaneous use of both signal separation methods (PICCA) is more effective when searching for transient events. The PCA method is more useful for long-term events, while the ICA method is better suited to recognize events of short duration. PICCA is a promising tool to detect transients of different characteristics in GNSS time series, which will be used in a next stage to generate a catalog of SSEs in Chile.