scholarly journals Fluid pressurisation and earthquake propagation in the Hikurangi subduction zone

2021 ◽  
Vol 12 (1) ◽  
Author(s):  
S. Aretusini ◽  
F. Meneghini ◽  
E. Spagnuolo ◽  
C. W. Harbord ◽  
G. Di Toro
Keyword(s):  
Shear Strength ◽  
Subduction Zone ◽  
Subduction Zones ◽  
Seismic Velocity ◽  
Fluid Pressure ◽  
Velocity Shear ◽  
Low Permeability ◽  
Earthquake Rupture ◽  
Seismic Slip ◽  
Saturated Clay

AbstractIn subduction zones, seismic slip at shallow crustal depths can lead to the generation of tsunamis. Large slip displacements during tsunamogenic earthquakes are attributed to the low coseismic shear strength of the fluid-saturated and non-lithified clay-rich fault rocks. However, because of experimental challenges in confining these materials, the physical processes responsible for the coseismic reduction in fault shear strength are poorly understood. Using a novel experimental setup, we measured pore fluid pressure during simulated seismic slip in clay-rich materials sampled from the deep oceanic drilling of the Pāpaku thrust (Hikurangi subduction zone, New Zealand). Here, we show that at seismic velocity, shear-induced dilatancy is followed by pressurisation of fluids. The thermal and mechanical pressurisation of fluids, enhanced by the low permeability of the fault, reduces the energy required to propagate earthquake rupture. We suggest that fluid-saturated clay-rich sediments, occurring at shallow depth in subduction zones, can promote earthquake rupture propagation and slip because of their low permeability and tendency to pressurise when sheared at seismic slip velocities.

scholarly journals Fluid pressurisation and earthquake propagation in the Hikurangi subduction zone

2021 ◽  
Author(s):  
Stefano Aretusini ◽  
Francesca Meneghini ◽  
Elena Spagnuolo ◽  
Christopher Harbord ◽  
Giulio Di Toro
Keyword(s):  
Shear Strength ◽  
Pore Pressure ◽  
Subduction Zone ◽  
Subduction Zones ◽  
Fluid Pressure ◽  
Low Permeability ◽  
Earthquake Rupture ◽  
Fault Rocks ◽  
Seismic Slip ◽  
Saturated Clay

<p>In subduction zones, seismic slip at shallow crustal depths can lead to the generation of tsunamis. Large slip displacements during tsunamogenic earthquakes are attributed to the low coseismic shear strength of the fluid-saturated and non-lithified clay-rich fault rocks. However, because of experimental challenges in confining these materials, the physical processes responsible of the coseismic reduction in fault shear strength are poorly understood. Using a novel experimental setup, we measured pore fluid pressure during simulated seismic slip in clay-rich materials sampled from the deep oceanic drilling of the Pāpaku thrust (Hikurangi subduction zone, New Zealand). Here we show that seismic slip is characterized by an initial decrease followed by an increase of pore pressure. The initial pore pressure decrease is indicative of dilatant behavior. The following pore pressure increase, enhanced by the low permeability of the fault, reduces the energy required to propagate earthquake rupture. We suggest that thermal and mechanical pressurisation of fluids facilitates seismic slip in the Hikurangi subduction zone, which was tsunamigenic about 70 years ago. Fluid-saturated clay-rich sediments, occurring at shallow depth in subduction zones, can promote earthquake rupture propagation and slip because of their low permeability and tendency to pressurise when sheared at seismic slip velocities.</p>

scholarly journals Seismic evidence for megathrust fault-valve behavior during episodic tremor and slip

2020 ◽  
Vol 6 (4) ◽  
pp. eaay5174 ◽  
Author(s):  
Jeremy M. Gosselin ◽  
Pascal Audet ◽  
Clément Estève ◽  
Morgan McLellan ◽  
Stephen G. Mosher ◽  
...  
Keyword(s):  
Subduction Zone ◽  
Subduction Zones ◽  
Seismic Velocity ◽  
Fluid Pressure ◽  
Pore Fluid ◽  
Cascadia Subduction Zone ◽  
Low Velocity ◽  
Seismic Scattering ◽  
Fault Valve ◽  
Episodic Tremor

Fault slip behavior during episodic tremor and slow slip (ETS) events, which occur at the deep extension of subduction zone megathrust faults, is believed to be related to cyclic fluid processes that necessitate fluctuations in pore-fluid pressures. In most subduction zones, a layer of anomalously low seismic wave velocities [low-velocity layer (LVL)] is observed in the vicinity of ETS and suggests high pore-fluid pressures that weaken the megathrust. Using repeated seismic scattering observations in the Cascadia subduction zone, we observe a change in the seismic velocity associated with the LVL after ETS events, which we interpret as a response to fluctuations in pore-fluid pressure. These results provide direct evidence of megathrust fault-valve processes during ETS.

Subduction zone heterogeneity observed across the magnitude spectrum for megathrust earthquakes

2021 ◽  
Author(s):  
Susan Bilek ◽  
Emily Morton
Keyword(s):  
Spatial Heterogeneity ◽  
Subduction Zone ◽  
Stress Drop ◽  
Subduction Zones ◽  
Cascadia Subduction Zone ◽  
Ocean Bottom ◽  
Small Earthquake ◽  
Megathrust Earthquakes ◽  
Seismic Slip ◽  
Subduction Zone Earthquakes

<p>Observations from recent great subduction zone earthquakes highlight the influence of spatial geologic heterogeneity on overall rupture characteristics, such as areas of high co-seismic slip, and resulting tsunami generation.  Defining the relevant spatial heterogeneity is thus important to understanding potential hazards associated with the megathrust. The more frequent, smaller magnitude earthquakes that commonly occur in subduction zones are often used to help delineate the spatial heterogeneity.  Here we provide an overview of several subduction zones, including Costa Rica, Mexico, and Cascadia, highlighting connections between the small earthquake source characteristics and rupture behavior of larger earthquakes.  Estimates of small earthquake locations and stress drop are presented in each location, utilizing data from coastal and/or ocean bottom seismic stations.  These seismicity characteristics are then compared with other geologic and geophysical parameters, such as upper and lower plate characteristics, geodetic locking, and asperity locations from past large earthquakes.  For example, in the Cascadia subduction zone, we find clusters of small earthquakes located in regions of previous seamount subduction, with variations in earthquake stress drop reflecting potentially disrupted upper plate material deformed as a seamount passed.  Other variations in earthquake location and stress drop can be correlated with observed geodetic locking variations. </p>

Bringing slow slip processes into focus

2020 ◽  
Author(s):  
Rebecca Bell
Keyword(s):  
High Resolution ◽  
Subduction Zones ◽  
Seismic Velocity ◽  
Numerical Models ◽  
Fluid Pressure ◽  
Geophysical Methods ◽  
Rock Properties ◽  
Slow Slip ◽  
Velocity Models ◽  
The North

<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>

scholarly journals Monitoring transient changes within overpressured regions of subduction zones using ambient seismic noise

2016 ◽  
Vol 2 (1) ◽  
pp. e1501289 ◽  
Author(s):  
Esteban J. Chaves ◽  
Susan Y. Schwartz
Keyword(s):  
Seismic Noise ◽  
Subduction Zones ◽  
Seismic Velocity ◽  
Fluid Pressure ◽  
Pore Fluid ◽  
Ambient Seismic Noise ◽  
Pore Fluid Pressure ◽  
Near Surface ◽  
Velocity Reduction ◽  
Pressurized Fluids

In subduction zones, elevated pore fluid pressure, generally linked to metamorphic dehydration reactions, has a profound influence on the mechanical behavior of the plate interface and forearc crust through its control on effective stress. We use seismic noise–based monitoring to characterize seismic velocity variations following the 2012 Nicoya Peninsula, Costa Rica earthquake [Mw(moment magnitude) 7.6] that we attribute to the presence of pressurized pore fluids. Our study reveals a strong velocity reduction (~0.6%) in a region where previous work identified high forearc pore fluid pressure. The depth of this velocity reduction is constrained to be below 5 km and therefore not the result of near-surface damage due to strong ground motions; rather, we posit that it is caused by fracturing of the fluid-pressurized weakened crust due to dynamic stresses. Although pressurized fluids have been implicated in causing coseismic velocity reductions beneath the Japanese volcanic arc, this is the first report of a similar phenomenon in a subduction zone setting. It demonstrates the potential to identify pressurized fluids in subduction zones using temporal variations of seismic velocity inferred from ambient seismic noise correlations.

Geometric controls on megathrust earthquakes

2020 ◽  
Vol 222 (2) ◽  
pp. 1270-1282
Author(s):  
Steven M Plescia ◽  
Gavin P Hayes
Keyword(s):  
Subduction Zone ◽  
Subduction Zones ◽  
Earthquake Hazard ◽  
Earthquake Rupture ◽  
Primary Control ◽  
Slab Geometry ◽  
Megathrust Earthquakes ◽  
Maximum Earthquake Magnitude ◽  
Maximum Earthquake

SUMMARY The role of subduction zone geometry in the nucleation and propagation of great-sized earthquake ruptures is an important topic for earthquake hazard, since knowing how big an earthquake can be on a given fault is fundamentally important. Past studies have shown subducting bathymetric features (e.g. ridges, fracture zones, seamount chains) may arrest a propagating rupture. Other studies have correlated the occurrence of great-sized earthquakes with flat megathrusts and homogenous stresses over large distances. It remains unclear, however, how subduction zone geometry and the potential for great-sized earthquakes (M 8+) are quantifiably linked—or indeed whether they can be. Here, we examine the potential role of subduction zone geometry in limiting earthquake rupture by mapping the planarity of seismogenic zones in the Slab2 subduction zone geometry database. We build from the observation that historical great-sized earthquakes have preferentially occurred where the surrounding megathrust is broadly planar, and we use this relationship to search for geometrically similar features elsewhere in subduction zones worldwide. Assuming geometry exerts a primary control on earthquake propagation and termination, we estimate the potential size distribution of large (M 7+) earthquakes and the maximum earthquake magnitude along global subduction faults based on geometrical features alone. Our results suggest that most subduction zones are capable of hosting great-sized earthquakes over much of their area. Many bathymetric features previously identified as barriers are indistinguishable from the surrounding megathrust from the perspective of slab curvature, meaning that they either do not play an important role in arresting earthquake rupture or that their influence on slab geometry at depth is not resolvable at the spatial scale of our subduction zone geometry models.

Mantle dynamics of the Andean Subduction Zone from continent-scale teleseismic S-wave tomography

2020 ◽  
Vol 224 (3) ◽  
pp. 1553-1571
Author(s):  
Emily E Rodríguez ◽  
Daniel Evan Portner ◽  
Susan L Beck ◽  
Marcelo P Rocha ◽  
Marcelo B Bianchi ◽  
...  
Keyword(s):  
South America ◽  
Subduction Zone ◽  
Lower Mantle ◽  
Subduction Zones ◽  
Seismic Velocity ◽  
High Amplitude ◽  
S Wave ◽  
Nazca Plate ◽  
Wave Tomography ◽  
Andean Subduction

SUMMARY The Andean Subduction Zone is one of the longest continuous subduction zones on Earth. The relative simplicity of the two-plate system has makes it an ideal natural laboratory to study the dynamics in subduction zones. We measure teleseismic S and SKS traveltime residuals at >1000 seismic stations that have been deployed across South America over the last 30 yr to produce a finite-frequency teleseismic S-wave tomography model of the mantle beneath the Andean Subduction Zone related to the Nazca Plate, spanning from ∼5°N to 45°S and from depths of ∼130 to 1200 km. Within our model, the subducted Nazca slab is imaged as a fast velocity seismic anomaly. The geometry and amplitude of the Nazca slab anomaly varies along the margin while the slab anomaly continues into the lower mantle along the entirety of the subduction margin. Beneath northern Brazil, the Nazca slab appears to stagnate at ∼1000 km depth and extend eastward subhorizontally for >2000 km. South of 25°S the slab anomaly in the lower mantle extends offshore of eastern Argentina, hence we do not image if a similar stagnation occurs. We image several distinct features surrounding the slab including two vertically oriented slow seismic velocity anomalies: one beneath the Peruvian flat slab and the other beneath the Paraná Basin of Brazil. The presence of the latter anomaly directly adjacent to the stagnant Nazca slab suggests that the plume, known as the Paraná Plume, may be a focused upwelling formed in response to slab stagnation in the lower mantle. Additionally, we image a high amplitude fast seismic velocity anomaly beneath the Chile trench at the latitude of the Sierras Pampeanas which extends from ∼400 to ∼1000 km depth. This anomaly may be the remnants of an older, detached slab, however its relationship with the Nazca–South America subduction zone remains enigmatic.

scholarly journals Cascadia megathrust earthquake rupture model constrained by geodetic fault locking

2021 ◽  
Vol 379 (2196) ◽  
pp. 20200135 ◽  
Author(s):  
Duo Li ◽  
Yajing Liu
Keyword(s):  
Subduction Zone ◽  
Subduction Zones ◽  
Seismogenic Zone ◽  
Cascadia Subduction Zone ◽  
Earthquake Rupture ◽  
Megathrust Earthquake ◽  
Fracture Dynamics ◽  
Coastal Marine ◽  
Rate And State Friction ◽  
Fault Locking

Paleo-earthquakes along the Cascadia subduction zone inferred from offshore sediments and Japan coastal tsunami deposits approximated to M9+ and ruptured the entire margin. However, due to the lack of modern megathrust earthquake records and general quiescence of subduction fault seismicity, the potential megathrust rupture scenario and influence of downdip limit of the seismogenic zone are still obscure. In this study, we present a numerical simulation of Cascadia subduction zone earthquake sequences in the laboratory-derived rate-and-state friction framework to investigate the potential influence of the geodetic fault locking on the megathrust sequences. We consider the rate-state friction stability parameter constrained by geodetic fault locking models derived from decadal GPS records, tidal gauge and levelling-derived uplift rate data along the Cascadia margin. We incorporate historical coseismic subsidence inferred from coastal marine sediments to validate our coseismic rupture scenarios. Earthquake rupture pattern is strongly controlled by the downdip width of the seismogenic, velocity-weakening zone and by the earthquake nucleation zone size. In our model, along-strike heterogeneous characteristic slip distance is required to generate margin-wide ruptures that result in reasonable agreement between the synthetic and observed coastal subsidence for the AD 1700 Cascadia Mw∼9.0 megathrust rupture. Our results suggest the geodetically inferred fault locking model can provide a useful constraint on earthquake rupture scenarios in subduction zones. This article is part of the theme issue ‘Fracture dynamics of solid materials: from particles to the globe’.

scholarly journals Seismic velocity structure at the southern termination of the 2016 Kumamoto Earthquake rupture, Japan

2020 ◽  
Vol 72 (1) ◽  
Author(s):  
Yasuhira Aoyagi ◽  
Haruo Kimura ◽  
Kazuo Mizoguchi
Keyword(s):  
Fault Zone ◽  
Velocity Structure ◽  
Seismic Velocity ◽  
Earthquake Rupture ◽  
2016 Kumamoto Earthquake ◽  
Seismic Velocity Structure ◽  
Kumamoto Earthquake ◽  
Seismogenic Layer ◽  
The 2016 Kumamoto Earthquake ◽  
Tectonic Line

Abstract The earthquake rupture termination mechanism and size of the ruptured area are crucial parameters for earthquake magnitude estimations and seismic hazard assessments. The 2016 Mw 7.0 Kumamoto Earthquake, central Kyushu, Japan, ruptured a 34-km-long area along previously recognized active faults, eastern part of the Futagawa fault zone and northernmost part of the Hinagu fault zone. Many researchers have suggested that a magma chamber under Aso Volcano terminated the eastward rupture. However, the termination mechanism of the southward rupture has remained unclear. Here, we conduct a local seismic tomographic inversion using a dense temporary seismic network to detail the seismic velocity structure around the southern termination of the rupture. The compressional-wave velocity (Vp) results and compressional- to shear-wave velocity (Vp/Vs) structure indicate several E–W- and ENE–WSW-trending zonal anomalies in the upper to middle crust. These zonal anomalies may reflect regional geological structures that follow the same trends as the Oita–Kumamoto Tectonic Line and Usuki–Yatsushiro Tectonic Line. While the 2016 Kumamoto Earthquake rupture mainly propagated through a low-Vp/Vs area (1.62–1.74) along the Hinagu fault zone, the southern termination of the earthquake at the focal depth of the mainshock is adjacent to a 3-km-diameter high-Vp/Vs body. There is a rapid 5-km step in the depth of the seismogenic layer across the E–W-trending velocity boundary between the low- and high-Vp/Vs areas that corresponds well with the Rokkoku Tectonic Line; this geological boundary is the likely cause of the dislocation of the seismogenic layer because it is intruded by serpentinite veins. A possible factor in the southern rupture termination of the 2016 Kumamoto Earthquake is the existence of a high-Vp/Vs body in the direction of southern rupture propagation. The provided details of this inhomogeneous barrier, which are inferred from the seismic velocity structures, may improve future seismic hazard assessments for a complex fault system composed of multiple segments.

Sign in / Sign up

Export Citation Format

Share Document