Imaging evolution of Cascadia slow-slip event using high-rate GPS

The slip history of short-term slow slip event (SSE) is typically inferred from daily Global Positioning System (GPS) data, which, however, cannot image the sub-daily processes, leaving the underlying mechanisms of SSEs elusive. To address the temporal resolution issue, we attempted to employ the kinematic subdaily GPS analysis, which has never been applied to SSE studies because its signal-to-noise ratio has been believed too low. By carefully post-processing sub-daily positions to remove non-tectonic position fluctuation, our 30-minute kinematic data clearly exhibits the transient motion of a few mm during one Cascadia SSE. A spatiotemporal slip image by inverting the 30-minute data exhibits a multi-stage evolution; it consists of an isotropic growth of SSE followed by an along-strike migration and termination within the rheologically controlled down-dip width. This transition at the slip growth mode is similar to the rupture growth of regular earthquakes, implying the presence of common mechanical factors behind the two distinct slip phenomena. The comparison with a slip inversion of the daily GPS demonstrates the current performance and limitation of the subdaily data in the SSE detection and imaging.Better understanding of the non-tectonic noise in the kinematic GPS analysis will further improve the temporal resolution of SSE.

Slow Slip Events Following the 2002 Mw 7.1 Hualien Offshore Earthquake Afterslip

The recurrence intervals of slow slip events may increase gradually after a large earthquake during the afterslip. Stress perturbations during coseismic and postseismic periods may result in such an increase of intervals. However, the increasing recurrence intervals of slow slip events are rarely observed during an afterslip. The evolution process along with the afterslip remains unclear. We report an observation of slow slip events following the 2002 Mw 7.1 Hualien offshore earthquake afterslip in the southernmost Ryukyu subduction zone. Slow slip events in 2005, 2009, and 2015 are adjacent to the Mw 7.1 earthquake hypocenter. An increasing slow-slip interval of 3.1, 4.2, and 6.2 years has been observed after the earthquake. We calculated coseismic and postseismic slips from the Mw 7.1 earthquake and then estimated the Coulomb stress changes in the slow slip region. The Mw 7.1 earthquake has contributed positive Coulomb stresses to both the 2005 slow-slip region and 2009/2015 repeating slow-slip region. The coseismic and postseismic Coulomb stress change on the 2005 slow-slip region is approximately 0.05 MPa and 0.035 MPa, respectively. However, both Coulomb stress changes on the 2009/2015 repeating slow-slip region are not over 0.03 MPa. The ongoing afterslip following the Mw 7.1 earthquake last for at least five years, evolving with a decaying stress rate with time. The long-term stress perturbations may be able to trigger the 2005 slow slip event during the afterslip. The 2009 slow slip event seems to be influenced by the afterslip as well. Postseismic stress evolution and frictional and stressed conditions of the slow-slip region can be a reason to affect the evolution process of slow slip events intervals.

Adjoint slip inversion under a constrained optimization framework: revisiting the 2006 guerrero slow slip event

Summary To shed light on the prevalently slow, aseismic slip interaction between tectonic plates, we developed a new static slip inversion strategy, the ELADIN (ELastostatic ADjoint INversion) method, that uses the adjoint elastostatic equations to compute the gradient of the cost function. ELADIN is a 2-step inversion algorithm to efficiently handle plausible slip constraints. First it finds the slip that best explains the data without any constraint, and then refines the solution by imposing the constraints through a Gradient Projection Method. To obtain a selfsimilar, physically-consistent slip distribution that accounts for sparsity and uncertainty in the data, ELADIN reduces the model space by using a von Karman regularization function that controls the wavenumber content of the solution, and weights the observations according to their covariance using the data precision matrix. Since crustal deformation is the result of different concomitant interactions at the plate interface, ELADIN simultaneously determines the regions of the interface subject to both stressing (i.e., coupling) and relaxing slip regimes. For estimating the resolution, we introduce a mobile checkerboard analysis that allows to determine lower-bound fault resolution zones for an expected slip-patch size and a given stations array. We systematically test ELADIN with synthetic inversions along the whole Mexican subduction zone and use it to invert the 2006 Guerrero Slow Slip Event (SSE), which is one of the most studied SSEs in Mexico. Since only 12 GPS stations recorded the event, careful regularization is thus required to achieve reliable solutions. We compared our preferred slip solution with two previously published models and found that our solution retains their most reliable features. In addition, although all three SSE models predict an upward slip penetration invading the seismogenic zone of the Guerrero seismic gap, our resolution analysis indicates that this penetration might not be a reliable feature of the 2006 SSE.

Observing seismic signatures of slow slip events with unsupervised learning

<p>Slow slip events are observed in geodetic data, and are occasionally associated with seismic signatures such as slow earthquakes (low-frequency earthquakes, tectonic tremors). In particular, it was shown that swarms of slow earthquake can correlate with slow slip events occurrence, and allowed to reveal the intermittent behavior of several slow slip events. This observation was possible thanks to detailed analysis of slow earthquakes catalogs and continuous geodetic data, but in every case, was limited to particular classes of seismic signatures. In the present study, we propose to infer the classes of seismic signals that best correlate with the observed geodetic data, including the slow slip event. We use a scattering network (a neural network with wavelet filters) in order to find meaningful signal features, and apply a hierarchical clustering algorithm in order to infer classes of seismic signal. We then apply a regression algorithm in order to predict the geodetic data, including slow slip events, from the occurrence of inferred seismic classes. This allow to (1) identify seismic signatures associated with the slow slip events as well as (2) infer the the contribution of each classes to the overall displacement observed in the geodetic data. We illustrate our strategy by revisiting the slow-slip event of 2006 that occurred beneath Guerrero, Mexico.</p>

Slow-slip, earthquake-swarms and fault-interactions at the western-end of the Hellenic Subduction System precede the Mw 6.9 Zakynthos Earthquake, Greece

<p>The month-to-year-long deformation of the Earth’s crust where active subduction zones terminate is poorly explored. Here we report on a multidisciplinary dataset that captures the synergy of slow-slip events, earthquake swarms and fault-interactions during the ~5 years leading up to the 2018 M<sub>w</sub> 6.9 Zakynthos Earthquake at the western termination of the Hellenic Subduction System (HSS). It appears that this long-lasting preparatory phase initiated due to a slow-slip event that lasted ~4 months and released strain equivalent to a ~M<sub>w</sub> 6.3 earthquake. We propose that the slow-slip event, which is the first to be reported in the HSS, tectonically destabilised the upper 20-40 km of the crust, producing alternating phases of seismic and aseismic deformation, including intense microseismicity (M<4) on neighbouring faults, earthquake swarms in the epicentral area of the M<sub>w</sub> 6.9 earthquake ~1.5 years before the main event, another episode of slow-slip immediately preceding the mainshock and, eventually, the large (M<sub>w </sub>6.9) Zakynthos Earthquake. Tectonic instability in the area is evidenced by a prolonged (~4 years) period of overall suppressed b-values (<1) and strong earthquake interactions on discrete strike-slip, thrust and normal faults. We propose that composite faulting patterns accompanied by alternating (seismic/aseismic) deformation styles may characterise multi-fault subduction-termination zones and may operate over a range of timescales (from individual earthquakes to millions of years).</p>

Slip bursts during coalescence of slow slip events in Cascadia

<p>Both laboratory experiments and dynamic simulations suggest that earthquakes can be preceded by a precursory phase of slow slip. Observing processes leading to an acceleration or spreading of slow slip along faults is therefore key to understand the dynamics potentially leading to seismic ruptures. Here, we use continuous GPS measurements of the ground displacement to image the daily slip along the fault beneath Vancouver Island during a slow slip event in 2013. We image the coalescence of three originally distinct slow slip fronts merging together. We show that during coalescence phases lasting for 2 to 5 days, the rate of energy (moment) release significantly increases. This observation supports the view proposed by theoretical and experimental studies that the coalescence of slow slip fronts is a possible mechanism for initiating earthquakes.</p>

Slow slip on the northern Hikurangi subduction interface, New Zealand

In October 2002, a surface displacement episode of 20-30 mm magnitude was observed over a ∼10 day period on two continuous Global Positioning System (GPS) instruments near Gisborne, North Island, New Zealand. We interpret this to result from slow slip on the northern Hikurangi subduction interface. Using ten years of regional campaign GPS (1995-2004) and recent continuous GPS data, we estimate the recurrence interval for similar events to be 2-3 yrs. In November 2004, a similar slow slip event occurred within this recurrence period. The 2002 event can be modeled by ∼18 cm of slow slip near the down-dip end of the seismogenic zone on the subduction interface offshore of Gisborne. The campaign GPS data show that the 2002 slow slip event had little effect on regional strain patterns. Copyright 2005 by the American Geophysical Union.