We analyzed continuous GPS data collected from 2002–2020 to characterize slow slip events (SSEs) in and near the Nicoya Peninsula, Costa Rica. These data are bisected by the 5 September 2012
Slow slip events (SSEs) have been observed in spatial and temporal proximity to megathrust earthquakes in various subduction zones, including the 2014
- NSF-PAR ID:
- 10497763
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- AGU Advances
- Volume:
- 5
- Issue:
- 2
- ISSN:
- 2576-604X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract M w 7.6 earthquake. The displacement time series contain multiple signals, including plate convergence, plate interface locking, coseismic and postseismic deformation, seasonal oscillations, SSEs, and noise. GPS‐measured coseismic and postseismic displacements associated with theM w 7.6 earthquake are modeled and removed by a step function plus multiple timescale relaxation processes with four characteristic times: 11, 94, 470, and 1,865 days. Seasonal oscillations are eliminated using a multichannel singular spectrum analysis (M‐SSA). Ten major SSEs (M w > 6.6) are observed in the remaining time series, with a constant recurrence interval of 21.7 ± 2.6 months. SSEs occur in both shallow (~10 km) and deep (~35 km) portions of the plate interface, but the latter last longer and have larger magnitudes. There is minimum to no slow slip in theM w 7.6 seismic rupture area and a persistent slow slip patch beneath the Nicoya Gulf entrance. Despite strong earthquake‐related stress perturbations, the inter‐SSE locking status on the megathrust is very similar between the late and early stages of the earthquake cycle and includes locked patches that ruptured in the 2012 earthquake or continue to rupture via SSEs. Some locked patches offshore south of the Nicoya Peninsula did not rupture in 2012, do not participate in SSEs, and may be indicative of supercycle behavior, that is, strain accumulation over several seismic cycles. These areas warrant heightened monitoring. -
more » « less
Abstract Fluids influence fault zone strength and the occurrence of earthquakes, slow slip events, and aseismic slip. We introduce an earthquake sequence model with fault zone fluid transport, accounting for elastic, viscous, and plastic porosity evolution, with permeability having a power‐law dependence on porosity. Fluids, sourced at a constant rate below the seismogenic zone, ascend along the fault. While the modeling is done for a vertical strike‐slip fault with 2D antiplane shear deformation, the general behavior and processes are anticipated to apply also to subduction zones. The model produces large earthquakes in the seismogenic zone, whose recurrence interval is controlled in part by compaction‐driven pressurization and weakening. The model also produces a complex sequence of slow slip events (SSEs) beneath the seismogenic zone. The SSEs are initiated by compaction‐driven pressurization and weakening and stalled by dilatant suctions. Modeled SSE sequences include long‐term events lasting from a few months to years and very rapid short‐term events lasting for only a few days; slip is ∼1–10 cm. Despite ∼1–10 MPa pore pressure changes, porosity and permeability changes are small and hence fluid flux is relatively constant except in the immediate vicinity of slip fronts. This contrasts with alternative fault valving models that feature much larger changes in permeability from the evolution of pore connectivity. Our model demonstrates the important role that compaction and dilatancy have on fluid pressure and fault slip, with possible relevance to slow slip events in subduction zones and elsewhere.
-
Abstract Shallow slow-slip events (SSEs) contribute to strain release near the shallow portions of subduction interfaces and may contribute to promoting shallow subduction earthquakes. Recent efforts in offshore monitoring of shallow SSEs have provided evidence of possible interactions between shallow SSEs and megathrust earthquakes. In this study, we use a dynamic earthquake simulator that captures both quasi-static (for SSEs) and dynamic (for megathrust earthquakes) slip to explore their interactions and implications for seismic and tsunami hazards. We model slip behaviors of a shallow-dipping subduction interface on which two locally locked patches (asperities) with different strengths are embedded within a conditionally stable zone. We find that both SSEs and earthquakes can occur, and they interact over multiple earthquake cycles in the model. Dynamic ruptures can nucleate on the asperities and propagate into the surrounding conditionally stable zone at slow speeds, generating tsunami earthquakes. A clear correlation emerges between the size of an earthquake and SSE activities preceding it. Small earthquakes rupture only the low-strength asperity, whereas large earthquakes rupture both. Before a large earthquake, periodic SSEs occur around the high-strength asperity, gradually loading stress into its interior. The critically stressed high-strength asperity can be ruptured together with the low-strength one in the large earthquake, followed by a relatively quiet interseismic period with very few SSEs and then a small earthquake. An SSE may or may not directly lead to nucleation of an earthquake, depending on whether a nearby asperity is ready for spontaneously dynamic failure. In addition, because of different SSE activities, the coupling degree may change dramatically between different interseismic periods, suggesting that its estimate based on a short period of observation may be biased.more » « less
-
more » « less
Abstract We study the effects of pore fluid pressure (
P f ) on the pre‐earthquake, near‐fault stress state, and 3‐D earthquake rupture dynamics through six scenarios utilizing a structural model based on the 2004M w 9.1 Sumatra‐Andaman earthquake. As pre‐earthquakeP f magnitude increases, effective normal stress and fault shear strength decrease. As a result, magnitude, slip, peak slip rate, stress drop, and rupture velocity of the scenario earthquakes decrease. Comparison of results with observations of the 2004 earthquake support that pre‐earthquakeP f averages near 97% of lithostatic pressure, leading to pre‐earthquake average shear and effective normal tractions of 4–5 and 22 MPa. The megathrust in these scenarios is weak, in terms of low mean shear traction at static failure and low dynamic friction coefficient during rupture. Apparent co‐seismic principal stress rotations and absolute post‐seismic stresses in these scenarios are consistent with the variety of observed aftershock focal mechanisms. In all scenarios, the mean apparent stress rotations are larger above than below the megathrust. Scenarios with largerP f magnitudes exhibit lower mean apparent principal stress rotations. We further evaluate pre‐earthquakeP f depth distribution. IfP f follows a sublithostatic gradient, pre‐earthquake effective normal stress increases with depth. IfP f follows the lithostatic gradient exactly, then this normal stress is constant, shifting peak slip and peak slip rate updip. This renders constraints on near‐trench strength and constitutive behavior crucial for mitigating hazard. These scenarios provide opportunity for future calibration with site‐specific measurements to constrain dynamically plausible megathrust strength andP f gradients. -
more » « less
Abstract Subduction zones host some of Earth's most damaging natural hazards, including megathrust earthquakes and earthquake‐induced tsunamis. A major control on the initiation and rupture characteristics of subduction megathrust earthquakes is how the coupled zone along the subduction interface accumulates elastic strain between events. We present results from observations of slow slip events (SSEs) in Cascadia occurring during the interseismic period downdip of the fully coupled zone, which imply that the orientation of strain accumulation within the coupled zone can vary with depth. Interseismic GPS motions suggest that forces derived from relative plate motions across a shallow, offshore locked plate interface dominate over decadal timescales. Deeper on the plate interface, below the locked (seismogenic) patch, slip during SSEs dominantly occurs in the updip direction, reflecting a dip‐parallel force acting on the slab, such as slab pull. This implies that in subduction zones with obliquely convergent plate motions, the seismogenic zone of the megathrust is loaded by forces acting in two discrete directions, leading to a depth‐varying orientation of strain accumulation on the plate interface.
