Slow Slip Research Articles

Juxtaposed slab dehydration, decarbonation and seismotectonic variation beneath the Philippine subduction zone based on 3-D modeling.

Largescale volcanic eruptions and earthquakes are occurring frequently in the Philippines, and research has shown that slab metamorphism and diversity alter the impacts of subducted oceanic plates by changing water‒carbon productivity and interplate stability. Within the framework of the thermal evolution history of subducting slabs, the relationships between subduction zone seismicity characterized by both regular megathrust earthquakes and slow slip events of various magnitudes and long-term slab dehydration-decarbonation evolution in the Philippines remain poorly understood. Here, we constructed a comprehensive thermal model incorporating 3-D slab geometric data for the incoming plate and a 3-D subduction velocity field based on the MORVEL plate motion dataset for the Philippine subduction zone with high spatial and temporal resolutions. Our findings reveal that subduction seismicity and arc volcanism are prominent in belt-shaped regions with high thermal gradients (> 5°C/km) and large-scale slab dehydration (> 0.05wt%/km). Dehydration of serpentinite in ultramafic rocks in the subducting slab and decarbonation of carbonate minerals preferentially contribute to the generation and transport of fluids and carbonate melts, thus facilitating seismicity and carbon-rich magmatism. Our results suggest that slab geometry diversity-induced juxtaposed slab dehydration-decarbonation processes play a vital role in the generation of megathrust earthquakes below the forearc.

Rheological Structure and Stress Triggered Megathrust Slip Constrained From the 2016 Mw 7.8 Kaikōura Crustal Earthquake

AbstractUnderstanding postseismic processes following the 2016 Mw 7.8 Kaikōura earthquake remains challenging due to the time‐dependent afterslip over the complex forearc system including crustal faults and megathrust, and the viscoelastic relaxation of the upper mantle. How the 2016 Mw 7.8 Kaikōura crustal earthquake interacts with the megathrust has yet to be better understood. Here we have derived the first 5‐year postseismic displacements from Global Positioning System (GPS) time series of 75 stations to study postseismic processes through a three‐dimensional viscoelastic finite element model. The optimal steady state viscosities of the crustal shear zone, megathrust shear zone, Australian upper mantle and Pacific upper mantle in the lowest‐misfit model among test models are 1018 Pa s, 4 × 1017 Pa s, 2 × 1019 Pa s and 1020 Pa s, respectively. The stress‐driven afterslip within the first 5 years after the earthquake is up to 80 cm over crustal faults, and up to 70 cm over the megathrust. A Kapiti slow slip sequence is probably promoted with a shorter interval by the 2016 earthquake, and is up to ∼11 cm within the first year after the earthquake. Afterslip over crustal faults and the megathrust are both required to reproduce the first‐order pattern of horizontal GPS observations. Coseismic rupture over the megathrust enhances shallow megathrust afterslip, which better fit the eastward postseismic displacement of sites near the rupture area. The southern end of Hikurangi megathrust may be activated during the 2016 earthquake and undergo continuous aseismic slip after the event.

Estimating Vertical Movement and Slip Distribution During the 2018 Boso, Japan, Slow Slip Event From Ocean Bottom Pressure Gauge Data and an Oceanic Model

AbstractMany slow slip events (SSEs) occur beneath the ocean, and continuous ocean‐bottom pressure gauge (OBP) observations provide useful data. OBPs record both oceanic variations and crustal movements, so we developed a multi‐channel singular spectrum analysis method to remove oceanic variations and applied our method to OBPs and oceanic model data. Then components of the oceanic model with good correlations to the OBP data were subtracted from the observed data. This method compensates for the incompleteness of the oceanic model and removes oceanic variations better than use of the original model. We applied the method to OBP data for the 2018 Boso, Japan, SSE to estimate its slip distribution. Comparing slip distributions obtained with and without the OBP data, we found that the distribution obtained using OBP data extended further offshore, and the offshore estimation error was smaller. Our study shows that offshore observations using OBPs are important for characterizing SSEs.

Effective Bulk Rheology of a Two‐Phase Subduction Shear Zone: Insights From Micromechanics‐Based Modeling and Implications for Subduction Interface Slow Slip Events

AbstractSubduction interfaces exhibit various slip styles, including slow slip events (SSEs). We use a micromechanics‐based approach to calculate the effective rheology of a shear zone containing ellipsoidal amphibolite clasts deforming by dislocation creep within an interconnected linear‐viscous phyllosilicate‐dominated matrix. Frictional failure occurs if local stress exceeds Mohr‐Coulomb yield strength. At moderate fluid overpressure, mixed‐frictional‐viscous behavior emerges at 350–560C, consistent with a broad zone of mixed fault slip behavior without requiring extreme fluid overpressures. Increasing stress in this transition zone promotes local frictional failure and raises bulk strain rate. If, however, the bulk strain rate increases by more than one order of magnitude, system‐wide frictional sliding becomes preferable. This strain rate increase is insufficient to explain the slip rates observed in geophysically detectable SSEs. Therefore, viscous matrix flow as modeled here cannot explain SSEs without either invoking dynamic weakening within a frictional‐viscous flow or a mechanism switch to dominantly frictional sliding.

Seasonal slow slip in landslides as a window into the frictional rheology of creeping shear zones.

Whether Earth materials exhibit frictional creep or catastrophic failure is a crucial but unresolved problem in predicting landslide and earthquake hazards. Here, we show that field-scale observations of sliding velocity and pore water pressure at two creeping landslides are explained by velocity-strengthening friction, in close agreement with laboratory measurements on similar materials. This suggests that the rate-strengthening friction commonly measured in clay-rich materials may govern episodic slow slip in landslides, in addition to tectonic faults. Further, our results show more generally that transient slow slip can arise in velocity-strengthening materials from modulation of effective normal stress through pore pressure fluctuations. This challenges the idea that episodic slow slip requires a narrow range of transitional frictional properties near the stability threshold, or pore pressure feedbacks operating on initially unstable frictional slip.

Detection of Hidden Low-Frequency Earthquakes in Southern Vancouver Island with Deep Learning

Low-frequency earthquakes (LFEs) are small-magnitude earthquakes that are depleted in high-frequency content relative to traditional earthquakes of the same magnitude. These events occur in conjunction with slow slip events (SSEs) and can be used to infer the space and time evolution of SSEs. However, because LFEs have weak signals, and the methods used to identify them are computationally expensive, LFEs are not routinely cataloged in most places. Here, we develop a deep-learning model that learns from an existing LFE catalog to detect LFEs in 14 years of continuous waveform data in southern Vancouver Island. The result shows significant increases in detection rates at individual stations. We associate the detections and locate them using a grid search approach in a 3D regional velocity model, resulting in over 1 million LFEs during the performing period. Our resulting catalog is consistent with a widely used tremor catalog during periods of large-magnitude SSEs. However, there are time periods where it registers far more LFEs than the tremor catalog. We highlight a 16-day period in May 2010, when our model detects nearly 3,000 LFEs, whereas the tremor catalog contains only one tremor detection in the same region. This suggests the possibility of hidden small-magnitude SSEs that are undetected by current approaches. Our approach improves the temporal and spatial resolution of the LFE activities and provides new opportunities to understand deep subduction zone processes in this region.

Parallel dynamics of slow slips and fluid-induced seismic swarms

Earthquake swarms may be driven by fluids, through hydraulic injections or natural fluid circulation, but also by slow and aseismic slip transients. Understanding the driving factors for these prolific sequences and how they can potentially develop into larger ruptures remains a challenge. A notable and almost ubiquitous feature of swarms is their hypocenters migration, which occurrence is closely related to the processes driving the observed seismicity, in a similar way as seismicity accompanies slow-slip events at subduction zones. Here, we analyze global data on migrating sequences, and identify scaling laws for migration velocity, moment and duration measured on natural and injection-induced swarms, foreshock sequences, and slow slip events. We highlight two different behaviors among these sequences: one linked to slow slips, with elevated migration velocities and moments, and the other related to fluid-induced processes, featuring lower velocities and moments. These results provide metrics for distinguishing between the drivers of earthquake swarms, fluid or slow-slip related, and prompt a reevaluation of scaling laws of fault slip transients, especially for swarms.

Use of tsunami observations for resolving coseismic slip of recent and historic large submarine earthquakes

Tsunamis generated by seafloor displacements accompanying large submarine earthquakes provide sensitivity to absolute slip position and distribution for offshore faulting analogous to that of geodetic observations for landward faulting. Tsunami recordings at deep‐water and near‐shore ocean bottom pressure sensors and tide gauges, along with runup and inundation measurements, can now be reliably modeled using detailed bathymetric structures and robust numerical codes. As a result, tsunami observations now play an important role in quantifying coseismic slip distributions for large submarine earthquakes in subduction zones and other tectonic environments. Applications of joint modeling or inversion of seismic, geodetic and tsunami observations for recent major earthquakes are described, highlighting the specific contributions of the tsunami observations to source model resolution. Tsunami observations provide unique information on the up‐dip extent of earthquake coseismic slip on subduction zone megathrust faults and occurrence of near‐trench slip, which are usually not well constrained by seismic and land‐based geodetic signals. Tsunami signals also help to detect offshore slow slip that is not evident in seismic or land‐based geodetic data and to balance geophysical constraints on ruptures that extend from on‐shore to off‐shore. Tsunami runup measurements and stratigraphic deposits further provide unique constraints on large earthquake ruptures that occurred prior to modern geophysical instrumentation.

Locating Boundaries Between Locked and Creeping Regions at Nankai and Cascadia Subduction Zones

AbstractInterseismic coupling maps and, especially, estimates of the location of the fully coupled (locked) zone relative to the trench, coastline, and slow slip events are crucial for determining megathrust earthquake hazard at subduction zones. We present an interseismic coupling inversion that estimates the locations of the upper and lower boundaries of the locked zone, the lower boundary of the deep transition zone, and downdip gradient of creep rate in the transition from locked to freely creeping in the downdip transition zone. We show that the locked zone at Cascadia is west of the coastline and 10 km updip of the slow slip zone along much of the margin, widest (25–125 km, extending to ∼19 km depth) in northern Cascadia, narrowest (0–70 km) in central Cascadia, with moment accumulation rate equivalent to a Mw 8.71 and Mw 8.85 earthquake for 300‐ and 500‐year earthquake cycles. We find a steep gradient in creep immediately below the locked zone, indicative of propagating creep, along the entire margin. At Nankai, we find three distinct zones of locking (offshore Shikoku, offshore southeast Kii peninsula, and offshore Shima peninsula) with a total moment accumulation rate equivalent to a Mw 8.70 earthquake for a 150‐year earthquake cycle. The bottom of the locked zone is nearly under the coastline for all three locked regions at Nankai and is positioned 0–5 km updip of the slow slip zone. In contrast with Cascadia, creep rate gradients below the locked zone at Nankai are generally gradual, consistent with stationary locking.

Frictional Strength and Frictional Instability of Glaucophane Gouges at Blueschist Temperatures Support Diverse Modes of Fault Slip From Slow Slip Events to Moderate‐Sized Earthquakes

AbstractFluid overpressure from the water released by subducted sediments and oceanic crust is an important mechanism for generating earthquakes via brittle failure and frictional instability. If unstable, such fault materials may also host diverse fault reactivation mechanisms from slow slip events to moderate‐sized earthquakes in cold subduction zones. We examine this hypothesis for glaucophane gouge ‐ a key index minerals for blueschist facies ‐ at lower confining stresses where behavior is poorly understood. We measure friction and stability at temperatures of 100°C–500°C and effective normal stresses of 50–200 MPa, to explore the controls of temperature, stresses and excess pore fluid pressures on fault friction. The frictional coefficient of glaucophane at representative temperatures and stresses is ∼0.70, insensitive to temperature but with a slight increase at lower effective stresses. Elevating temperature promotes a transition from velocity‐strengthening to weak velocity‐weakening behavior, indicating the destabilizing effect of high temperature downdip in subduction zones. Reducing effective normal stress or concomitantly elevating pore fluid pressure further strengthens the velocity‐weakening response and would be manifest as moderate‐sized earthquakes. Our results support the potential for enhanced unstable sliding of glaucophane gouges at lower effective stresses and blueschist facies temperatures ‐ with weak to moderate velocity‐weakening responses upon fluid pressurization vital for understanding the abundance of slow slip to moderate‐sized earthquakes apparent in cold subduction zones.

Fault‐Valve Instability: A Mechanism for Slow Slip Events

AbstractGeophysical and geological studies provide evidence for cyclic changes in fault‐zone pore fluid pressure that synchronize with or at least modulate slip events. A hypothesized explanation is fault valving arising from temporal changes in fault zone permeability. In our study, we investigate how the coupled dynamics of rate and state friction, along‐fault fluid flow, and permeability evolution can produce slow slip events. Permeability decreases with time, and increases with slip. Linear stability analysis shows that steady slip with constant fluid flow along the fault zone is unstable to perturbations, even for velocity‐strengthening friction with no state evolution, if the background flow is sufficiently high. We refer to this instability as the “fault valve instability.” The propagation speed of the fluid pressure and slip pulse, which scales with permeability enhancement, can be much higher than expected from linear pressure diffusion. Two‐dimensional simulations with spatially uniform properties show that the fault valve instability develops into slow slip events, in the form of aseismic slip pulses that propagate in the direction of fluid flow. We also perform earthquake sequence simulations on a megathrust fault, taking into account depth‐dependent frictional and hydrological properties. The simulations produce quasi‐periodic slow slip events from the fault valve instability below the seismogenic zone, in both velocity‐weakening and velocity‐strengthening regions, for a wide range of effective normal stresses. A separation of slow slip events from the seismogenic zone, which is observed in some subduction zones, is reproduced when assuming a fluid sink around the mantle wedge corner.

Unstable Sliding of Plagioclase Gouge and Deformation Mechanisms Under Hydrothermal Conditions With Effective Normal Stresses of 100–300 MPa

AbstractPlagioclase feldspar is a major mineral in mafic crustal rocks. To better understand the deformation mechanism of plagioclase feldspar during frictional faulting, we conducted shearing experiments on simulated plagioclase gouge in a wide range of effective normal stress of 100–300 MPa, pore‐water pressure of 30–100 MPa, and temperatures ranging from 100°C to 600°C. The coefficient of friction is found to range from 0.65 to 0.74 across the entire temperature range, showing no significant thermal weakening process. Except for a case at 200°C with an effective normal stress of 300 MPa, the frictional sliding is velocity weakening over the whole temperature range, showing a steady‐state rate dependence (a−b) ranging from −0.5 × 10−3 to −8.6 × 10−3. This property facilitates nucleation of unstable slips in frictional faulting. Above 200°C, the direct rate effect parameter (a) and the evolution effect parameter (b) of friction increase with temperature up to a threshold of 400°C or 500°C, depending on the effective normal stress. This thermal enhancement suggests thermally activated creep at contact junctions governed by intergranular pressure solution, as evidenced by microstructural signatures indicating the prevalence of very fine precipitates formed at the surfaces of gouge particles as a result of pressure solution. In frictional sliding of plagioclase, a low effective normal stress of 100 MPa corresponds to a higher degree of velocity weakening and tends to facilitate seismic slip rather than slow slips, whereas the high effective normal stress of 300 MPa corresponds to a minor velocity weakening which may cause slow‐slip events in faults of limited size.

Experimental Study on Acoustic Emission Features and Energy Dissipation Properties during Rock Shear-Slip Process.

The features of rock shear-slip fracturing are closely related to the stability of rock mass engineering. Granite, white sandstone, red sandstone, and yellow sandstone specimens were selected in this study. The loading phase of "shear failure > slow slip > fast slip" was set up to explore the correlation between fracture type, acoustic emission (AE) features, and energy dissipation during the rock fracturing process. The results show that there is a strong correlation between fracture type, energy dissipation, and AE features. The energy dissipation ratio of tension-shear (T-S) composite, shear, and tensile types is 10:100:1. The fracture types in the shear failure phase are mainly tensile and TS composite types. The differential mechanism of energy dissipation of different rocks during the shear-slip process is revealed from the physical property perspectives of mineral composition, particle size, and diagenetic mode. These results provide a necessary research basis for energy dissipation research in rock failure and offer an important scientific foundation for analyzing the fracture propagation problem in the shear-slip process. They also provide a research basis for further understanding the acoustic emission characteristics and crack type evolution during rock shear and slip processes, which helps to better understand the shear failure mechanism of natural joints and provides a reference for the identification of precursors of shear disasters in geotechnical engineering.

Global subduction slow slip events and associated earthquakes.

Three decades of geodetic monitoring have established slow slip events (SSEs) as a common mode of fault slip, sometimes linked with earthquake swarms and in a few cases escalating to major seismic events. However, the connection between SSEs and earthquake hazard has been difficult to quantify and contextualize beyond regional studies. We aggregate a geodetic record of SSEs from subduction zones in the circum-Pacific region. In aggregate, earthquake rates increase up to threefold concurrent with and proximal to SSEs. The relative amplitude of this increase is correlated with the SSE size and, to a lesser extent, their depth and region. The subdued and coincident earthquake response to SSE stress transfer suggests a more limited role of static stress transfer and a very short relaxation timescale for the triggered seismicity. The observed range of behavior does not support a major connection between SSEs and earthquake hazard.

Non-Stationary ETAS Model: How It Works for External Forcing

Abstract The stationary epidemic-type aftershock sequence (ETAS) model is based on the important empirical laws in aftershock statistics with a self-similar feature, and is therefore useful for the statistical analysis of many common earthquake occurrence series. However, because earthquake catalogs become richer, the non-stationarity arising from geophysical heterogeneity becomes more pronounced. This article discusses the utility of a model that assumes time dependence of the first two ETAS parameters, which are sensitive to the short-term prediction. The inversion analyses of the non-stationary ETAS model in this article show heuristic results for swarm earthquakes and complex mainshock–aftershock-type seismicity. The case studies demonstrate how these parameters change quantitatively in swarm seismicity associated with slow slips, increases in pore-fluid pressure such as magma and hydrothermal fluids, stress changes associated with an earthquake motion, and interseismic-induced effects due to geological properties.

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.

Tectonic geomorphology and paleoseismology of the Angelochori fault segment of the Anthemountas extensional detachment fault, Central Macedonia, Greece: Paleoseismic evidence from the 1677 CE earthquake

This study presents the results of a geomorphological and paleoseismological investigation on the Angelochori normal fault. This fault controls the southern boundary of the Anthemountas basin and constitutes the western segment of the Anthemountas extensional detachment fault (AEDF), a ∼48 km-long, E-W-trending, north-dipping low-angle normal fault in Central Macedonia, Greece. To assess the relative tectonic activity of the Angelochori fault and to determine the number, timing, and displacement of surface-faulting earthquakes, we conducted detailed field mapping and paleoseismological investigations. Our approach involved utilizing geomorphic indices, analyzing stratigraphic and structural relations, considering luminescence ages, and employing OxCal modeling techniques. The acquired data reveal, for the first time, the Middle-Late Holocene activity of the Anghelochori fault, providing constraints on the age of three strong (M > 6) surface-rupturing earthquakes over the last ∼7.7 ka. The youngest of these events aligns with the historical earthquake of the Anthemountas basin in 1677 CE. The observed minimum displacement per event ranges from 0.48 m to 0.65 m, consistent with a fault length estimated at about 24–28 km, suggesting an earthquake magnitude on the order of 6.6. The cumulative slip rate over the recorded period averages 0.21 ± 0.02 mm/a, with an average open recurrence interval of 2.6 ± 0.2 ka. Nonetheless, the displacement pattern displays alternating intervals of rapid slip at 0.60 ± 0.16 mm/a and slower slip at 0.08 ± 0.01 mm/a. Classified as Class 1, the Angelochori fault signifies an active tectonic structure posing a significant seismic hazard for the Thessaloniki metropolitan area.

Along-strike forearc and subducted upper slab structure beneath north Chile: Slow slip implications

We evaluate the lateral variation of the lithospheric structure beneath the forearc of the North Chile subduction zone along strike, extending from Iquique to Valparaiso, to understand a possible seismic structure that hosts deep slow earthquakes by comparing the deep slow slip areas around 60 km depth, using receiver function and inversion method. Relative differences are investigated in seismic velocities around the plate interface, especially underlying the mantle wedge, the upper part of the subducted oceanic crust, and the oceanic mantle along the subduction strike. We estimate and compare the thickness and velocity of subducted oceanic crust, slab top, subducted oceanic Moho, continental Moho, and wedge mantle velocity beneath the seismic stations. The oceanic crustal thickness range from ∼5 to 10 km respectively, and the oceanic Moho depth is ∼40 km–∼90 km. It was observed that the disparities in velocities were more pronounced near the south of the Iquique zone (18°S-25°S), with high subducting oceanic crust velocities. These discrepancies in velocities could be regarded as the presence of fluids from subducted ridges and seamounts, which might alter seismicity in the area, including slow slip.

Spatiotemporal Evolution of Slow Slip Events at the Offshore Hikurangi Subduction Zone in 2019 Using GNSS, InSAR, and Seafloor Geodetic Data

AbstractDetecting crustal deformation during transient deformation events at offshore subduction zones remains challenging. The spatiotemporal evolution of slow slip events (SSEs) on the offshore Hikurangi subduction zone, New Zealand, during February–July 2019, is revealed through a time‐dependent inversion of onshore and offshore geodetic data that also accounts for spatially varying elastic crustal properties. Our model is constrained by seafloor pressure time series (as a proxy for vertical seafloor deformation), onshore continuous Global Navigation Satellite System (GNSS) data, and Interferometric Synthetic Aperture Radar displacements. Large GNSS displacements onshore and uplift of the seafloor (10–33 mm) require peak slip during the event of 150 to >200 mm at 6–12 km depth offshore Hawkes Bay and Gisborne, comparable to maximum slip observed during previous seafloor pressure deployments at north Hikurangi. The onshore and offshore data reveal a complex evolution of the SSE, over a period of months. Seafloor pressure data indicates the slow slip may have persisted longer near the trench than suggested by onshore GNSS stations in both the Gisborne and Hawkes Bay regions. Seafloor pressure data also reveal up‐dip migration of SSE slip beneath Hawke Bay occurred over a period of a few weeks. The SSE source region appears to coincide with locations of the March 1947 Mw 7.0–7.1 tsunami earthquake offshore Gisborne and estimated great earthquake rupture sources from paleoseismic investigations offshore Hawkes Bay, suggesting that the shallow megathrust at north and central Hikurangi is capable of both seismic and aseismic rupture.

A physical explanation for an unusually long-duration slow slip event in the Nankai Trough

The Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) and the Long Term Borehole Monitoring System (LTBMS), installed above the source region of the 1944 Tonankai earthquake, revealed that crustal deformation is driven by slow slip events (SSEs) in the shallower extension of megathrust earthquakes. However, there are unresolved questions about (A) the duration of the SSE in February 2012, which was longer than expected for SSEs with similar magnitudes, and (B) the relationship of the spatial distribution of fault slip between the SSEs in February and December 2012 under the condition of drilling disturbance. To clarify these questions, we re-analyzed the pore/seafloor pressure data associated with the SSEs. Our refined fault models show that the SSE in February had a significantly slower propagation speed and a longer duration than others, while the SSE in December was comparable to others. We interpret that the difference in duration and propagation speed is related to external and internal stress perturbations, respectively. Using the ocean modeling JCOPE (Japan Coastal Ocean Predictability Experiment), we identified that the decrease and subsequent increase in seafloor pressure due to the passage of the Kuroshio meander coincided with the latter part of the longer duration of the SSE in February and its termination, respectively. This suggests that the Kuroshio meander might affect the duration of SSEs. Our refined fault model also indicates that the amount of shear stress accumulation was small before the occurrence of the SSE in February, which triggered the slow propagation of aseismic slip based on a rate- and state-dependent friction law. These results imply that we need to consider the variety of SSEs from the viewpoint of stress perturbation due to not only interaction between fault segments but also external forces from oceanographic phenomena.