Integrated rupture mechanics for slow slip events and earthquakes

Author(s):  
Huihui Weng

Abstract Slow slip events usually occur downdip of seismogenic zones in subduction megathrusts and crustal faults, with rupture speeds much slower than earthquakes. The empirical moment-duration scaling relation can help constrain the physical mechanism of slow slip events, yet it is still debated whether this scaling is linear or cubic and a fundamental model unifying slow slip events and earthquakes is still lacking. Here I present numerical simulations that show that slow slip events are regular earthquakes with negligible dynamic-wave effects. A continuum of rupture speeds, from arbitrarily-slow speeds up to the S-wave speed, is primarily controlled by the stress drop and a transition slip rate above which the fault friction transitions from rate-weakening behaviour to rate-strengthening behaviour. This continuum includes tsunami earthquakes, whose rupture speeds are about one-third of the S-wave speed. These numerical simulation results are predicted by the three-dimensional theory of dynamic fracture mechanics of elongated ruptures. This fundamental model unifies slow slip events and earthquakes, reconciles the observed moment-duration scaling relations, and opens new avenues for understanding earthquakes through investigations of the kinematics and dynamics of frequently occurring slow slip events.

2021 ◽  
Author(s):  
Huihui Weng

Abstract Slow slip events usually occur downdip of seismogenic zones in subduction megathrusts and crustal faults, with rupture speeds much slower than earthquakes. The empirical moment-duration scaling relation can help constrain the physical mechanism of slow slip events, yet it is still debated whether this scaling is linear or cubic and a fundamental model unifying slow slip events and earthquakes is still lacking. Here I present numerical simulations that show that slow slip events are regular earthquakes with negligible dynamic-wave effects. A continuum of rupture speeds, from arbitrarily-slow speeds up to the S-wave speed, is primarily controlled by the stress drop and a transition slip rate above which the fault friction transitions from rate-weakening behaviour to rate-strengthening behaviour. This continuum includes tsunami earthquakes, whose rupture speeds are about one-third of the S-wave speed. These numerical simulation results are predicted by the three-dimensional theory of dynamic fracture mechanics of elongated ruptures. This fundamental model unifies slow slip events and earthquakes, reconciles the observed moment-duration scaling relations, and opens new avenues for understanding earthquakes through investigations of the kinematics and dynamics of frequently occurring slow slip events.


2019 ◽  
Vol 219 (3) ◽  
pp. 2074-2096 ◽  
Author(s):  
Kazuro Hirahara ◽  
Kento Nishikiori

Summary A variety of slow slip events at subduction zones have been observed. They can be stress meters for monitoring the stress state of megathrust faults during their earthquake cycles. In this study, we focus on long-term slow slip events (LSSEs) recurring at downdip portions of megathrust faults among such slow earthquakes. Data analyses and simulation studies of LSSEs have so far been executed independently. In atmosphere and ocean sciences, data assimilations that optimally combine data analyses and simulation studies have been developed. We develop a method for estimating frictional properties and monitoring slip evolution on an LSSE fault, with a sequential data assimilation method, the ensemble Kalman filter (EnKF). We executed numerical twin experiments for the Bungo Channel LSSE fault in southwest Japan to validate the method. First, based on a rate- and state-dependent friction law, we set a rate-weakening circular LSSE patch on the rate-strengthening flat plate interface, whose critical nucleation size is larger than that of the patch, and reproduced the observed Bungo Channel LSSEs with recurrence times of approximately 7 yr and slip durations of 1 yr. Then, we synthesized the observed data of surface displacement rates at uniformly distributed stations with noises from the simulated slip model. Using our EnKF method, we successfully estimated the frictional parameters and the slip rate evolution after a few cycles. Secondly, we considered the effect of the megathrust fault existing in the updip portion of the LSSE fault, as revealed by kinematic inversion studies of Global Navigation Satellite System (GNSS) data and added this locked region with a slip deficit rate in the model. We estimated the slip rate on the locked region only kinematically, but the quasi-dynamic equation of motion in each LSSE fault cell includes the stress term arising from the locked region. Based on this model, we synthesized the observed surface displacement rate data for the actual distribution of GNSS stations and executed EnKF estimations including the slip rate on the locked region. The slip rate on the locked region could be quickly retrieved. Even for the actual distribution of GNSS stations, we could successfully estimate frictional parameters and slip evolution on the LSSE fault. Thus, our twin numerical experiments showed the validity of our EnKF method, although we need further studies for actual GNSS data analyses.


2020 ◽  
Vol 13 (10) ◽  
pp. 705-710
Author(s):  
Kyungjae Im ◽  
Demian Saffer ◽  
Chris Marone ◽  
Jean-Philippe Avouac

2019 ◽  
Vol 124 (12) ◽  
pp. 13250-13268 ◽  
Author(s):  
Bunichiro Shibazaki ◽  
Laura M. Wallace ◽  
Yoshihiro Kaneko ◽  
Ian Hamling ◽  
Yoshihiro Ito ◽  
...  

1974 ◽  
Vol 64 (6) ◽  
pp. 1789-1808 ◽  
Author(s):  
Robert Burridge ◽  
Cesar Levy

abstract It has recently been shown (Burridge, 1973) that in two dimensions plane-strain shear cracks lacking cohesion may run at speeds ranging from the Rayleigh-wave to the S-wave speed or possibly at the P-wave speed. On the other hand, it has long been known that in antiplane strain, cracks lacking cohesion must run at least at the S-wave speed. Since locally at the edge of a three-dimensional crack there is a combination of plane and antiplane strain, we have sought and found solutions for circular shear cracks expanding at the S-wave speed and at the P-wave speed. These have finite shear tractions ahead of the crack and so may correspond to frictional sliding in the absence of cohesion. The analysis combines the method of Kostrov (1964b) with the results of Burridge (1973). We carry out a complete evaluation for the displacements and tractions on the fault plane, and far-field radiation for the S-wave-speed crack. The corresponding evaluations for the P-wave speed are not elementary and are not attempted here. As far as the authors are aware, these are the first analytic solutions of three-dimensional crack problems which satisfy a physically plausible fracture criterion for failure under shear.


2018 ◽  
Vol 4 (10) ◽  
pp. eaat8472 ◽  
Author(s):  
N. Voss ◽  
T. H. Dixon ◽  
Z. Liu ◽  
R. Malservisi ◽  
M. Protti ◽  
...  

Slow slip events have been suggested to trigger subduction earthquakes. However, examples to date have been poorly recorded, occurring offshore, where data are sparse. Better understanding of slow slip events and their influence on subsequent earthquakes is critical for hazard forecasts. We analyze a well-recorded event beginning 6 months before the 2012 Mw (moment magnitude) 7.6 earthquake in Costa Rica. The event migrates to the eventual megathrust rupture. Peak slip rate reached a maximum of 5 mm/day, 43 days before the earthquake, remaining high until the earthquake. However, changes in Mohr-Coulomb failure stress at the hypocenter were small (0.1 bar). Our data contradict models of earthquake nucleation that involve power law acceleration of slip and foreshocks. Slow slip events may prove useful for short-term earthquake forecasts.


2021 ◽  
Author(s):  
Mattia Pizzati ◽  
Fabrizio Balsamo ◽  
Fabrizio Storti

<p>Valuable information concerning the seismic cycle are mainly provided by the study of exposed fossil subduction-accretionary complexes and by coring and probing through present-day active major plate boundary interfaces. Subduction zone investigation and monitoring allowed to comprehend the mechanics of thrust-related faulting and to discern seismic events with different slip rate (coseismic events, slow slip events and tremor). While subduction zones received particular attention especially following the M<sub>w</sub> 9 Tohoku-Oki earthquake in Japan, relatively small-scale extensional faults affecting the uppermost portion of seismogenic zone of the Earth’s crust are still less studied.</p><p>Here, we present a field and laboratory study of meso-scale structures recorded within the fault core of an extensional fault zone (Rocca di Neto fault, offset < 100 m) affecting Pleistocene siliciclastic sediments in the Crotone Basin, Calabria, Southern Italy. Due to shallow burial conditions experienced by deformed sediments (< 400-500 m), the fault zone structure is characterised by deformation features typical of high-porosity granular rocks, with extensive occurrence of deformation bands, subsidiary faults and gouges. The 1 m-thick fault core displays a complex network of mutually cross-cutting black gouges and deformation bands developed in foliated sand. Some black gouges have straight pattern parallel to the master fault surface, while others are displaced and dragged along the deformation bands (mm-offset). Black gouges, previously interpreted as coseismic events due to moderate to high-temperature mineral assemblage, are characterised by cm-offset and extreme grain comminution via severe cataclasis (mean grain size of 20-30 μm and fractal dimension from 3.0 to 3.3); clast preferred orientation is almost parallel to the gouge outer boundaries, thus resulting in a well-developed foliation. Deformation bands are organised in two conjugate sets and display moderate to intense cataclasis depending on the accommodated displacement (mean grain size of 80-170 μm and fractal dimension from 2.4 to 2.8), with preferred orientation of clasts describing an angle of 30-45° from the band surface. Within deformation bands the foliation is less defined compared to black gouges. At the intersections between gouges and deformation bands, the resulting tectonic fabric is given by the superposition of different deformation events overprinting the original one.</p><p>The difference in grain size distribution, fractal dimension, clast shape preferred orientation (i.e., foliation) and mineral composition between black gouges and deformation bands supports the hypothesis of different slip rates causing their development. In particular, black gouges are interpreted to develop during coseismic slip (~0.1-1 m/s), while deformation bands formed during interseismic intervals (slip rate from nm/s to μm/s). The cross-cutting relationship between gouges and deformation bands, combined with the overprinting of different tectonic fabrics along the intersections, suggests they formed as a result of repeating coseismic (fast slip) and aseismic (slow slip) events occurring at shallow burial-near surface conditions. This feature could be a key point to evaluate the deformation style (fast vs slow slip) and to estimate the potential seismic hazard of superficial faults affecting high-porosity sediments.</p>


2021 ◽  
Author(s):  
Olivia Pita Sllim

<p><b>The Hikurangi subduction zone beneath the eastern North Island, New Zealand exhibits a variety of fault-slip related phenomena including tsunami earthquakes, non-volcanic tremor, low-frequency earthquakes, episodic slow slip, and repeating earthquakes. The northern Hikurangi margin hosts shallow slow-slip and is weakly coupled to shallow depths. In contrast, the southern Hikurangi margin is strongly coupled, and only deep slow-slip has been observed. The transition in coupling occurs beneath the township of Porangahau, and is an exemplary focus region for studying how this change in locking is accommodated. </b></p><p>To examine slip processes beneath Porangahau, we have constructed and analysed a catalogue of repeating earthquakes that occurred between 2004 and 2018. Repeating earthquakes are thought to re-rupture the same fault patch at different times, and thus have nearly identical waveforms, locations and magnitudes. Because repeating earthquakes represent cyclic loading, they can be used to detect temporal and spatial changes of slip-rate at depth and hence monitor how stress is transferred to seismogenic zones. </p><p>To build a catalogue of repeating earthquakes we first clustered the GeoNet earthquake catalogue by distance and correlation to identify potential repeating events. We then used a stronger cross-correlation threshold of at least 0.95 normalised cross-correlation value at three or more stations to identify repeating earthquakes from the initial clusters. This threshold was determined by our group's previous work on the northern Hikurangi margin. We identified 225 families of repeating earthquakes, with each family having two or more earthquakes in the 14-year study period from 2004 to 2018. </p><p>We carried out manual phase picking and polarity identification for the most recent event in each family and computed absolute locations, local magnitudes calibrated with moment magnitude, and high-quality focal mechanisms. For the rest of the events in each family, we conducted cross-correlation re-picking to obtain precise relative locations and relative magnitudes. With precise locations and well-constrained focal mechanisms, we determined whether the repeating earthquake families originated within the Pacific Plate, Australian Plate or on the subduction interface. Most of the families are located within the Pacific Plate, and the majority of families that originate on the subduction interface are located near the township of Porangahau. At least 220 of the 532 identified repeating earthquakes locate at the transition from strong- to weak-coupling of the subduction interface near the township of Porangahau. </p><p>A variety of slow slip events have been detected near Porangahau in the last two decades. Even though some repeating earthquakes correlate spatially and temporally with slow slip events, temporal and spatial correlations between slow slip events and repeating earthquakes are scarce and sparse. The majority of repeating earthquakes are located up-dip or down-dip of modelled slow slip patches, with very few families having spatial correlation with slow slip events. We obtained a moment-recurrence interval relationship for the catalogue of repeating earthquakes near Porangahau and compared it to the relationship obtained by Nadeau and Johnson (1998) at Parkfield, California. Finally, we computed slip-rates using the families located on the subduction interface and obtained an average slip-rate of 13 mm/yr. The insights gained from this study lay the groundwork for future work constraining processes of strain accumulation at the creeping-to-locked transition zone near Porangahau.</p>


2021 ◽  
Author(s):  
Romain Jolivet ◽  
Bertrand Rouet-Leduc ◽  
Jorge Jara ◽  
Manon Dalaison ◽  
Claudia Hulbert ◽  
...  

&lt;p&gt;While some faults remain locked for tens to hundreds of years, some active faults slip slowly, either continuously or episodically. The discovery of slow, generally silent, slip at the turn of the century led to a profound modification of our understanding of the mechanics of faulting, shedding light on the dynamics of fault slip. Such dynamics areis controlled by the past history of stress along the fault plane (i.e. historical ruptures), fluids circulating in the crust and the rheology of the crust and fault plane. Understanding the influence of these different factors requires dense observations, as suggested by the large range of spatial and temporal scales involved in the control of the slip velocity along a fault. Specifically, the smallest scales of slow slip have beenwere inferred by the observation of tremors or low frequency events, interpreted as the chatter of a fault plane while it slips slowly. We are missing direct observations of such kilometer-scale slow slip events and continental creeping faults are an obvious target for such observationsfor such observations.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;Aseismic slip along the North Anatolian Fault was recognized in the 1960&amp;#8217;s by the observation of offset man-made features without earthquakes recorded. Following these early observations, multiple geodetic studies focused on recording aseismic slip and analyzed the average rate of shallow slow slip in the vicinity of the town of Ismetpasa. GPS, InSAR and creepmeter data all converge toward an aseismic slip rate reaching 1 cm/yr in places, with significant along- strike variations. Furthermore, earlyHowever, creepmeter measurements in the 80&amp;#8217;s, confirmed by records from a more recent instrument, suggest aseismic slip is currently episodic, occurring in bursts of slip. Recent InSAR data from the Cosmo-SkyMed constellation captured a month-long slow slip event with a maximum of 2 cm/yr of slip.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;We propose to analyze the geodetic record to search for slow slip events over the 2015-2020 period. We take advantage of a dense network of continuous GNSS stations installed in 2017 and of time series of Sentinel 1 SAR data to identify at least 3 slow slip events along the North Anatolian Fault. Thanks to the dense temporal sampling of the GNSS records, we describe faithfullyobserve the onset of slow slip. We use a deep learning algorithm to extract the surface signature of the slow slip events from the InSAR time series, highlighting a slow rupture front propagating along strike. We compare the occurrences of slow slip events with the local fault geometry, the average distribution of kinematic coupling and the historical seismicity. We discuss the mechanical implications of such detailed description of slow slip along an active fault. In conclusion, while slow slip rate averaged over periods longer than 2-3 years seems constant over the last 40 years, identification of slow slip events suggests this apparently constant rate results from slow slip events over multiple spatial and temporal scales.&lt;/p&gt;


2021 ◽  
Author(s):  
Olivia Pita Sllim

<p><b>The Hikurangi subduction zone beneath the eastern North Island, New Zealand exhibits a variety of fault-slip related phenomena including tsunami earthquakes, non-volcanic tremor, low-frequency earthquakes, episodic slow slip, and repeating earthquakes. The northern Hikurangi margin hosts shallow slow-slip and is weakly coupled to shallow depths. In contrast, the southern Hikurangi margin is strongly coupled, and only deep slow-slip has been observed. The transition in coupling occurs beneath the township of Porangahau, and is an exemplary focus region for studying how this change in locking is accommodated. </b></p><p>To examine slip processes beneath Porangahau, we have constructed and analysed a catalogue of repeating earthquakes that occurred between 2004 and 2018. Repeating earthquakes are thought to re-rupture the same fault patch at different times, and thus have nearly identical waveforms, locations and magnitudes. Because repeating earthquakes represent cyclic loading, they can be used to detect temporal and spatial changes of slip-rate at depth and hence monitor how stress is transferred to seismogenic zones. </p><p>To build a catalogue of repeating earthquakes we first clustered the GeoNet earthquake catalogue by distance and correlation to identify potential repeating events. We then used a stronger cross-correlation threshold of at least 0.95 normalised cross-correlation value at three or more stations to identify repeating earthquakes from the initial clusters. This threshold was determined by our group's previous work on the northern Hikurangi margin. We identified 225 families of repeating earthquakes, with each family having two or more earthquakes in the 14-year study period from 2004 to 2018. </p><p>We carried out manual phase picking and polarity identification for the most recent event in each family and computed absolute locations, local magnitudes calibrated with moment magnitude, and high-quality focal mechanisms. For the rest of the events in each family, we conducted cross-correlation re-picking to obtain precise relative locations and relative magnitudes. With precise locations and well-constrained focal mechanisms, we determined whether the repeating earthquake families originated within the Pacific Plate, Australian Plate or on the subduction interface. Most of the families are located within the Pacific Plate, and the majority of families that originate on the subduction interface are located near the township of Porangahau. At least 220 of the 532 identified repeating earthquakes locate at the transition from strong- to weak-coupling of the subduction interface near the township of Porangahau. </p><p>A variety of slow slip events have been detected near Porangahau in the last two decades. Even though some repeating earthquakes correlate spatially and temporally with slow slip events, temporal and spatial correlations between slow slip events and repeating earthquakes are scarce and sparse. The majority of repeating earthquakes are located up-dip or down-dip of modelled slow slip patches, with very few families having spatial correlation with slow slip events. We obtained a moment-recurrence interval relationship for the catalogue of repeating earthquakes near Porangahau and compared it to the relationship obtained by Nadeau and Johnson (1998) at Parkfield, California. Finally, we computed slip-rates using the families located on the subduction interface and obtained an average slip-rate of 13 mm/yr. The insights gained from this study lay the groundwork for future work constraining processes of strain accumulation at the creeping-to-locked transition zone near Porangahau.</p>


Sign in / Sign up

Export Citation Format

Share Document