Episodic stress tensor and fluid pressure cycling in subducting oceanic crust during Northern Hikurangi slow slip events

Author(s):  
Emily Warren-Smith ◽  
Bill Fry ◽  
Laura Wallace ◽  
Enrique Chon ◽  
Stuart Henrys ◽  
...  

<p><span>The occurrence of slow slip events (SSEs) in subduction zones has been proposed to be linked to the presence of, and fluctuations in near-lithostatic fluid pressures (P</span><sub><span>f</span></sub><span>) within the megathrust shear zone and subducting oceanic crust. In particular, the 'fault-valve' model is commonly used to describe occasional, repeated breaching of a low-permeability interface shear zone barrier, which caps an overpressured hydrothermal fluid reservoir. In this model, a precursory increase in fluid pressure may therefore be anticipated to precede megathrust rupture. Resulting activation of fractures during slip opens permeable pathways for fluid migration and fluid pressure decreases once more, until the system becomes sealed and overpressure can re-accumulate. While the priming conditions for cyclical valving behaviour have been observed at subduction zones globally, and evidence for post-megathrust rupture drainage exists, physical observations of precursory fluid pressure increases, and subsequent decreases, particularly within the subducting slab where hydrothermal fluids are sourced, remain elusive. </span></p><p><span>Here we use earthquake focal mechanisms recorded on an ocean-bottom seismic network to identify changes in the stress tensor within subducting oceanic crust during four SSEs in New Zealand’s Northern Hikurangi subduction zone. We show that the stress, or shape ratio, which describes the relative magnitudes of the principal compressive stress axes, shows repeated decreases prior to, and rapid increases during the occurrence of geodetically documented SSEs. We propose that these changes represent precursory accumulation and subsequent release of fluid pressure within overpressured subducting oceanic crust via a ‘valving’ model for megathrust slip behaviour. Our observations indicate that the timing of slow slip events on subduction megathrusts may be controlled by cyclical accumulation of fluid pressure within subducting oceanic crust.</span></p><p><span>Our model is further supported by observations of seismicity preceding a large SSE in the northern Hikurangi Margin in 2019, captured by ocean-bottom seismometer</span><span>s</span><span> and </span><span>absolute </span><span>pressure </span><span>recorders.</span> <span>O</span><span>bservations of microseismicity </span><span>during this period </span><span>indicate that a stress state conducive to vertical fluid flow was present in the downgoing plate prior to SSE initiation, before subsequently returning to a</span><span> down-dip</span><span> extensional state following the SSE. We propose this precursory seismicity is indicative of fluid migration towards the interface shear zone from the lower plate fluid reservoir, which may have helped triggering slip on the megathrust. </span></p><p><span>We also present preliminary results of a moment tensor study to investigate spatial and temporal patterns in earthquake source properties in SSE regions along the Hikurangi Margin. In particular, earthquakes near Porangahau – a region susceptible to dynamic triggering of tremor and where </span><span>shallow </span><span>SSEs occur every 5 years or so – exhibit distinctly lower double couple components than elsewhere along the margin. We </span><span>attribute this to elevated fluid pressures within the crust here, which is consistent with recent observations of high seismic reflectivity from an autocorrelation study. Such high fluid pressure may control the broad range of seismic and aseismic phenomena observed at Porangahau. </span></p>

2020 ◽  
Author(s):  
Claudio Petrini ◽  
Luca Dal Zilio ◽  
Taras Gerya

<p>Slow slip events (SSEs) are part of a spectrum of aseismic processes that relieve tectonic stress on faults. Their occurrence in subduction zones have been suggested to trigger megathrust earthquakes due to perturbations in fluid pressure. However, examples to date have been poorly recorded and physical observations of temporal fluid pressure fluctuations through slow slip cycles remain elusive. Here, we use a newly developed two-phase flow numerical model — which couples solid rock deformation and pervasive fluid flow — to show how crustal stresses and fluid pressures within subducting megathrust evolve before and during slow slip and regular events. This unified 2D numerical framework couples inertial mechanical deformation and fluid flow by using finite difference methods, marker-in-cell technique, and poro-visco-elasto-plastic rheologies. Furthermore, an adaptive time stepping allows the correct resolution of both long- and short-time scales, ranging from years to milliseconds during the dynamic propagation of earthquake rupture.</p><p>Here we show how permeability and its spatial distribution control the degree of locking along the megathrust interface and the interplay between seismic and aseismic slip. While a constant permeability leads to more regular seismic cycles, a depth dependent permeability contributes substantially to the development of two distinct megathrust zones: a shallow, locked seismogenic zone and a deep, narrow aseismic segment characterized by SSEs. Furthermore, we show that without requiring any specific friction law, our model shows that permeability, episodic stress transfer and fluid pressure cycling control the predominant slip mode along the subduction megathrust. Specifically, we find that the up-dip propagation of episodic SSEs systematically decreases the fault strength due to a continuous accumulation and release of fluid pressure within overpressured subducting interface, thus affecting the timing of large megathrust earthquakes. These results contribute to improve our understanding of the physical driving forces underlying the interplay between seismic and aseismic slip, and demonstrate that slow slip events may prove useful for short-term earthquake forecasts.</p>


Geosphere ◽  
2019 ◽  
Vol 15 (4) ◽  
pp. 1038-1065 ◽  
Author(s):  
Alissa J. Kotowski ◽  
Whitney M. Behr

Abstract We use structural and microstructural observations from exhumed subduction-related rocks exposed on Syros Island (Cyclades, Greece) to provide constraints on the length scales and types of heterogeneities that occupy the deep subduction interface, with possible implications for episodic tremor and slow slip. We selected three Syros localities that represent different oceanic protoliths and deformation conditions within a subduction interface shear zone, including: (1) prograde subduction of oceanic crust to eclogite facies; (2) exhumation of oceanic crust from eclogite through blueschist-greenschist facies; and (3) exhumation of mixed mafic crust and sediments from eclogite through blueschist-greenschist facies. All three localities preserve rheological heterogeneities that reflect metamorphism of primary lithological, geochemical, and/or textural variations in the subducted protoliths and that take the form of brittle pods and lenses within a viscous matrix. Microstructural observations indicate that the matrix lithologies (blueschists and quartz-rich metasediments) deformed by distributed power-law viscous flow accommodated by dislocation creep in multiple mineral phases. We estimate bulk shear zone viscosities ranging from ∼1018 to 1020 Pa-s, depending on the relative proportion of sediments to (partially eclogitized) oceanic crust. Eclogite and coarse-grained blueschist heterogeneities within the matrix preserve multiple generations of dilational shear fractures and veins formed under high-pressure conditions. The veins commonly show coeval or overprinting viscous shear, suggesting repeated cycles of frictional and viscous strain. These geologic observations are consistent with a mechanical model of episodic tremor and slow slip (ETS), in which the deep subduction interface is a rheologically heterogeneous distributed shear zone comprising transiently brittle (potentially tremor-genic) sub-patches within a larger, viscously creeping interface patch. Based on our observations of outcrop and map areas of heterogeneous patches and the sizes, distributions, and amounts of brittle offset recorded by heterogeneities, we estimate that simultaneous brittle failure of heterogeneities could produce tremor bursts with equivalent seismic moments of 4.5 × 109–4.7 × 1014 N m, consistent with seismic moments estimated from geophysical data at active subduction zones.


2020 ◽  
Author(s):  
Cailey Condit ◽  
◽  
Victor Guevara ◽  
Jonathan R. Delph ◽  
Melodie French

2015 ◽  
Vol 16 (11) ◽  
pp. 4006-4020 ◽  
Author(s):  
E. K. Mitchell ◽  
Y. Fialko ◽  
K. M. Brown

2020 ◽  
Author(s):  
Rebecca Bell

<p>The discovery of slow slip events (SSEs) at subduction margins in the last two decades has changed our understanding of how stress is released at subduction zones. Fault slip is now viewed as a continuum of different slip modes between regular earthquakes and aseismic creep, and an appreciation of seismic hazard can only be realised by understanding the full spectrum of slip. SSEs may have the potential to trigger destructive earthquakes and tsunami on faults nearby, but whether this is possible and why SSEs occur at all are two of the most important questions in earthquake seismology today. Laboratory and numerical models suggest that slow slip can be spontaneously generated under conditions of very low effective stresses, facilitated by high pore fluid pressure, but it has also been suggested that variations in frictional behaviour, potentially caused by very heterogeneous fault zone lithology, may be required to promote slow slip.</p><p>Testing these hypotheses is difficult as it requires resolving rock properties at a high resolution many km below the seabed sometimes in km’s of water, where drilling is technically challenging and expensive. Traditional geophysical methods like travel-time tomography cannot provide fine-scale enough velocity models to probe the rock properties in fault zones specifically. In the last decade, however, computational power has improved to the point where 3D full-waveform inversion (FWI) methods make it possible to use the full wavefield rather than just travel times to produce seismic velocity models with a resolution an order of magnitude better than conventional models. Although the hydrocarbon industry have demonstrated many successful examples of 3D FWI the method requires extremely high density arrays of instruments, very different to the 2D transect data collection style which is still commonly employed at subduction zones.</p><p> The north Hikurangi subduction zone, New Zealand is special, as it hosts the world’s most well characterised shallow SSEs (<2 km to 15 km below the seabed).  This makes it an ideal location to collect 3D data optimally for FWI to resolve rock properties in the slow slip zone. In 2017-2018 an unprecedentedly large 3D experiment including 3D multi-channel seismic reflection, 99 ocean bottom seismometers and 194 onshore seismometers was conducted along the north Hikurangi margin in an 100 km x 15 km area, with an average 2 km instrument spacing. In addition, IODP Expeditions 372 and 375 collected logging-while drilling and core data, and deployed two bore-hole observatories to target slow slip in the same area. In this presentation I will introduce you to this world class 3D dataset and preliminary results, which will enable high resolution 3D models of physical properties to be made to bring slow slip processes into focus.  </p>


2020 ◽  
Author(s):  
Mathilde Radiguet ◽  
Ekaterina Kazachkina ◽  
Louise Maubant ◽  
Nathalie Cotte ◽  
Vladimir Kostoglodov ◽  
...  

<p>Slow slip events (SSEs) represent a significant mechanism of strain release along several subduction zones, and understanding their occurrence and relations with major earthquake asperities is essential for a comprehensive understanding of the seismic cycle. Here, we focus on the Mexican subduction zone, characterized by the occurrence of recurrent large slow slip events (SSEs), both in the Guerrero region, where the SSEs are among the largest observed worldwide, and in the Oaxaca region, where smaller, more frequent SSEs occur. Up to now, most slow slip studies in the Mexican subduction zone focused either on the detailed analysis of a single event, were limited to a small area (Guerrero or Oaxaca), or were limited to data before 2012 [e.g.1-4]. In this study, our aim is to build an updated and consistent catalog of major slow slip events in the Guerrero-Oaxaca region.</p><p>We use an approach similar to Michel et al. 2018 [5]. We analyze the GPS time series from 2000 to 2019 using Independent Component Analysis (ICA), in order to separate temporally varying sources of different origins (seasonal signals, SSEs and afterslip of major earthquakes). We are able to isolate a component corresponding to seasonal loading, which matches the temporal evolution of displacement modeled from the GRACE data. The sources (independent components) identified as tectonic sources of deep origin are inverted for slip on the subduction interface. We thus obtain a model of the spatio-temporal evolution of aseismic slip on the subduction interface over 19 years, from which we can isolate around 30 individual slow slip events of M<sub>w </sub>> 6.2.</p><p> The obtained catalog is coherent with previous studies (in terms of number of events detected, magnitude and duration) which validates the methodology. The observed moment-duration scaling is close to M<sub>0</sub>~T<sup>3 </sup>as recently suggested by Michel [6] for Cascadia SSEs, and our study extends the range of magnitude considered in their analysis. Finally, we also investigate the spatio-temporal relations between the SSEs occurring in the adjacent regions of Guerrero and Oaxaca, and their interaction with local and distant earthquakes.</p><p> </p><p>References:</p><ol><li>Kostoglodov, V. et al. A large silent earthquake in the Guerrero seismic gap, Mexico. Geophys. Res. Lett <strong>30</strong>, 1807 (2003).</li> <li>Graham, S. et al. Slow Slip History for the Mexico Subduction Zone: 2005 Through 2011. Pure and Applied Geophysics 1–21 (2015). doi:10.1007/s00024-015-1211-x</li> <li>Larson, K. M., Kostoglodov, V. & Shin’ichi Miyazaki, J. A. S. The 2006 aseismic slow slip event in Guerrero, Mexico: New results from GPS. Geophys. Res. Lett. <strong>34</strong>, L13309 (2007).</li> <li>Radiguet, M. et al. Slow slip events and strain accumulation in the Guerrero gap, Mexico. J. Geophys. Res. <strong>117</strong>, B04305 (2012).</li> <li>Michel, S., Gualandi, A. & Avouac, J.-P. Interseismic Coupling and Slow Slip Events on the Cascadia Megathrust. Pure Appl. Geophys. (2018). doi:10.1007/s00024-018-1991-x</li> <li>Michel, S., Gualandi, A. & Avouac, J. Similar scaling laws for earthquakes and Cascadia slow-slip events. Nature <strong>574, </strong>522–526 (2019) doi:10.1038/s41586-019-1673-6</li> </ol><p> </p>


2020 ◽  
Author(s):  
Giuseppe Costantino ◽  
Mauro Dalla Mura ◽  
David Marsan ◽  
Sophie Giffard-Roisin ◽  
Mathilde Radiguet ◽  
...  

<p>The deployment of increasingly dense geophysical networks in many geologically active regions on the Earth has given the possibility to reveal deformation signals that were not detectable beforehand. An example of these newly discovered signals are those associated with low-frequency earthquakes, which can be linked with the slow slip (aseismic slip) of faults. Aseismic fault slip is a crucial phenomenon as it might play a key role in the precursory phase before large earthquakes (in particular in subduction zones), during which the seismicity rate grows as well as does the ground deformation. Geodetic measurements, e.g. the Global Positioning System (GPS), are capable to track surface deformation transients likely induced by an episode of slow slip. However, very little is known about the mechanisms underlying this precursory phase, in particular regarding to how slow slip and seismicity relate.</p><p>The analysis done in this work focuses on recordings acquired by the Japan Meteorological Agency in the Boso area, Japan. In the Boso peninsula, interactions between seismicity and slow slip events can be observed over different time spans: regular slow slip events occur every 4 to 5 years, lasting about 10 days, and are associated with a burst of seismicity (Hirose et al. 2012, 2014, Gardonio et al. 2018), whereas an accelerated seismicity rate has been observed over decades that is likely associated with an increasing shear stress rate (i.e., tectonic loading) on the subduction interface (Ozawa et al. 2014, Reverso et al. 2016, Marsan et al. 2017).</p><p>This work aims to explore the potential of  Deep Learning  for better characterizing the interplay between seismicity and ground surface deformation. The analysis is based on a data-driven approach for building a model for assessing if a link seismicity – surface deformation exists and to characterize the nature of this link. This has potentially strong implications, as (small) earthquakes are the prime observable, so that better understanding the seismicity rate response to potentially small slow slip (so far undetected by GPS) could help monitoring those small slow slip events. The statistical problem is expressed as a regression between some features extracted from the seismic data and the GPS displacements registered at one or more stations.</p><p>The proposed method, based on a Long-Short Term Memory (LSTM) neural network, has been designed in a way that it is possible to estimate which features are more relevant in the estimation process. From a geophysical point of view, this can provide interesting insights for validating the results, assessing the robustness of the algorithms and giving insights on the underlying process. This kind of approach represents a novelty in this field, since it opens original perspectives for the joint analysis of seismic / aseismic phenomena with respect to traditional methods based on more classical geophysical data exploration.</p>


2020 ◽  
Author(s):  
Frederique Rolandone ◽  
Jean-Mathieu nocquet ◽  
Patricia Mothes ◽  
Paul Jarrin ◽  
Mathilde Vergnolle

<p>In subduction zones, slip along the plate interface occurs in various modes including earthquakes, steady slip, and transient accelerated aseismic slip during either Slow Slip Events (SSE) or afterslip. We analyze continuous GPS measurements along the central Ecuador subduction segment to illuminate how the different slip modes are organized in space and time in the zone of the 2016 Mw 7.8 Pedernales earthquake. The early post-seismic period (1 month after the earthquake) shows large and rapid afterslip developing at discrete areas of the megathrust and a slow slip event remotely triggered (∼100 km) south of the rupture of the Pedernales earthquake. We find that areas of large and rapid early afterslip correlate with areas of the subduction interface that had hosted SSEs in years prior to the 2016 earthquake. Areas along the Ecuadorian margin hosting regular SSEs and large afterslip had a dominant aseismic slip mode that persisted throughout the earthquake cycle during several years and decades: they regularly experienced SSEs during the interseismic phase, they did not rupture during the 2016 Pedernales earthquake, they had large aseismic slip after it. Four years after the Pedernales earthquake, postseismic deformation is still on-going. Afterslip and SSEs are both involved in the postseimsic deformation. Two large aftershocks (Mw 6.7 & 6.8) occurred after the first month of postseismic deformation in May 18, and later in July 7 2016 two other large aftershocks (Mw 5.9 & 6.3) occurred, all were located north east of the rupture. They may have triggered their own postseismic deformation. Several seismic swarms were identified south and north of the rupture area by a dense network of seismic stations installed during one year after the Pedernales earthquakes, suggesting the occurrence of SSEs. Geodetically, several SSEs were detected during the postseismic deformation either in areas where no SSEs were detected previously, or in areas where regular seismic swarms and repeating earthquakes were identified. The SSEs may have been triggered by the stress increment due to aftershocks or due to afterslip.</p>


2021 ◽  
Author(s):  
Lise Alalouf ◽  
Yajing Liu

<p>Subduction zones are where the largest earthquakes occur. In the past few decades, scientists have also discovered the presence of episodic aseismic slip, including slow slip events (SSEs), along most of the subduction zones. However, it is still unclear how these SSEs can influence megathrust earthquake ruptures. The Costa Rica subduction zone is a particularly interesting area because a SSE was recorded 6 months before the 2012 Mw7.6 earthquake in the Nicoya Peninsula, suggesting a potential stress transfer from the SSE to the earthquake slip zone. SSEs beneath the Nicoya Peninsula were also recorded both updip and downdip the seismogenic zone, making it a unique area to study the complex interaction between SSEs and earthquakes.</p><p>As most of the shallow SSEs were recorded around the Nicoya Peninsula, we chose to start using a 1D planar fault embedded in a homogeneous elastic half-space, with different dipping angles following several geometric profiles of the subduction fault beneath the Nicoya Peninsula section of the Costa Rica margin. This 1D modelling study allows us to better investigate the interaction between shallow and deep SSEs and megathrust earthquakes with high numerical resolution and relatively short computation time. The model provides information on the long-term seismic history by reproducing the different stages of the seismic cycle (interseismic slip, shallow and deep episodic slow slip, and coseismic slip).</p><p>We study the influence of the variation of numerical parameters and frictional properties on the recurrence interval, maximum slip velocity and cumulative slip of SSEs (both shallow and deep) and earthquakes and their interaction with each other. We then compare our results with GPS and seismic observations (i.e. cumulative slip, characteristic duration, moment rate, depth and size of the rupture, equivalent magnitude) to identify an optimal set of model parameters to understand the interaction between various modes of subduction fault deformation.</p>


Sign in / Sign up

Export Citation Format

Share Document