Potential role of fluids in modulating stable and unstable slip

1. Potential role of fluids in modulating stable and unstable slip Nadaya Cubas GE 277/169b • Liu and Rice (2005), Aseismic slip transients emerge spontaneously in three-dimensional rate and state modeling of subduction earthquake sequences, JGR • Liu and Rice (2007), Spontaneous and triggered aseismic deformation transients in a subduction fault model, JGR

2. Slow slip events Common features: • slip velocity 10 to 100* Vplate(10-9 to 10-7 m/s) • transients migrate along trench at speed from km/day to km/year • limited to a depth range near and below end of the seismogenic zone. • quasiperiodicity of aseismic transients observed in some shallow subduction zones from few months up to several (10) years • Associated with non-volcanic tremors • Evidence for high Pore Fluid Pressure

3. Slow slip events Liu and Rice 2005: Slow slip events can arise spontaneously • rate and state formulation, with T (depth) dependence (a-b) parameter • for certain (low) effective normal stress with the simplest representation of metamorphic fluid release (typically 50 MPa or 100 MPa) in the seismogenic zone and downdip region.

4. Model 3D thrust fault with rate- and state- dependent friction law Dietrich-Ruina

5. Parameters: (1) Temperature (depth) dependence of material parameter (a-b), for granite gouge under hydrothermal conditions Transition weakening/strengthening (2) P lower than lithostatic – constant with time, sn eff from 50 to 150MPa (100MPa) (3) Lf = 120mm

6. Results: in the seismogenic zone with small perturbations of (a-b) parameters Non uniform slip mode, Complex slip, with strongly locked zones

7. Results: near and down dip the transition zone Aseismic transients locked Vplate downdip Transient, Vtransient = 3km/year 10*Vplate shallow aseismic slip during interseismic

8. Common characteristics compared to GPS observations: Transients have typical aseismic slip rates averaging around 10 to 102 times Vpl. Transients migrate along strike at a speed of several to tens of km/yr in simulations. Transients are limited to depths around the downdip end of seismogenic zone defined in the model. They generally cannot rupture into the updip portion because the fault there is very firmly locked. Nor can slip be enhanced in the velocity strengthening region that is further downdip. CCL: Can arise spontaneously near neutral stability Transition of frictional properties from velocity weakening to velocity strengthening, combined with great earthquake sequences that leave zones of heterogenous slip and stress alteration along strike, seem to be the cause of the aseismic slip event occurrence. They may reflect episodic relaxation, during the interseismic period, of stress concentrated around the transition zone.

9. Common characteristics compared to GPS observations: CCL: Can arise spontaneously near neutral stability Transition of frictional properties from velocity weakening to velocity strengthening, combined with great earthquake sequences that leave zones of heterogenous slip and stress alteration along strike, seem to be the cause of the aseismic slip event occurrence. They may reflect episodic relaxation, during the interseismic period, of stress concentrated around the transition zone. However: Long transient: tens of years rather than the several weeks span observed in some natural events - Large characteristic slip distance L, with smaller L : speed increase - The speed with which transient slipping zones propagate along the strike increases with the reduction of s; e.g., that speed increases by a factor of 3 when the assumed s is reduced by a factor of 2, from 200 MPa to 100.

10. Slow slip events • Liu and Rice 2007: • rate and state formulation, with T (depth) dependence (a-b) parameter • 2D to test smaller L • Lithostatic pore pressure • Reproduce oscillatory behavior, ranging from • 3 months beneath the Bungo Channel, southwest Japan • 11-19 months along the Cascadia • 2–3 years near Gisborne, along the Hikurangi subduction zone in New Zealand • 4 years in Guerrero, Mexico • 6 years beneath the Bungo Channel • 10 years or more, Alaska, Sunda Megathrust between the Mentawai

11. Evidences for HPP: 1.petrological analysis of phase equilibrium: dehydration conditions would be met for shallow dipping subduction zones around 350C and above, which means near and downdip of the frictional stability transition Fluids released from such dehydration reactions can significantly increase the pore pressure

12. Here, add one slide about HPP shallow NE japan – thermal pressurization presented before.

13. Model 2D thrust fault with rate- and state- dependent friction law • Realistic model: • effective normal stress is high (finite) in the seismogeniczone, • much lower from stability transition and further downdip. b) Simplified model: • RW patch completely locked (s effinfinitely high) • W extending up dip of the stability transition • Whole down-dip region at a much lower, uniform s eff, due to dehydration.

14. Parameters: (1) Temperature (depth) dependence of material parameter (a-b), for granite gouge under hydrothermal conditions (2) P near lithostatic (3) Lf = 20mm

15. Results: with locked seismogenic zone Oscillation depends on W/h* W: (a-b) < 0 but very high pore pressure: seff very low h* critical cell size [Rice, 1993; Lapustaet al., 2000] Stable slip patch size for steady sliding with rate-weakening friction Shear modulus Poisson’s ratio

16. Results: with locked seismogenic zone Oscillation depends on W/h* With L = 30mm, seff = 2-3 Mpa: 14 months period However, oscillatory slip rate amplitude extremely low case b and c: lower than Vplate

17. Results: with unlocked seismogenic zone Higher amplitude Rate weakening zone as being at a finitely higher seff so that slip is allowed to penetrate naturally into it by a small amount show that the amplitude of the oscillatory slips increases substantially. Average period: 2 years Slip rate: 2 to 4 * Vplate Accumulated slip: 1 to 2 cm

18. Can arise spontaneously near neutral stability Or transient sequences can be triggered • from the stress field alteration due to nearby earthquakes, • or the pore pressure variations due to episodes of metamorphic fluid release along the subducting slab, • or along-strike variations in seismic slips.

19. Transient sequences can be triggered by slab seismicity aseismic transients might act as a spatial-temporal connection between distant seismicity clusters in subduction zones NF TF Transient aseismic slip NF TF + Repetition in 2006

20. Transients response will dependent on: Perturbation Introduction Time (Ttp /Tr0 ) : Accelerate or delay next thrust EQ Relative Location Along the Thrust Fault: Slightly varies interseismic time Magnitude of the Stress Perturbation: Pore pressure rise or drop can trigger transient >> And have direct consequences on te timing of future thrust EQ

21. Ccl Our studies can produce transients in abundance, of a wide range of recurrence intervals, including transients of relatively short periods (1 year) in cases for which there is a zone of near-lithostatic fluid pressurization around and downdip of the frictional stability transition, a situation independently supported by tremor triggering and petrology. However: However, with the simple model of a completely locked seismogenic zone updip of width W, the amount of deformation at the earth’s surface associated with the spontaneous short-term transients is at least an order of magnitude smaller than observed in Cascadia. Further, even with a relatively more realistic model of an unlocked seismogenic zone as in section 3.2, most of the aseismic moment in our modeling is contributed over a zone extending updip and moderately downdip from the stability transition, whereas a much broader downdip zone is thought to be activated during the transient events in the northern Cascadia and Guerrero subduction zones [Dragert et al., 2001; Hirn and Laigle, 2004; Larson et al., 2004 Have we made the downdip rate-strengthening zone too stable? Segall et al. (2010), Dilatant strengthening as a mechanism for slow slip event, JGR