Dilatant strengthening as a mechanism for slow slip events

The mechanics of slow slip events (SSE) in subduction zones remain unresolved. We suggest that SSE nucleate in areas of unstable friction under drained conditions, but as slip accelerates dilatancy reduces pore pressure p quenching instability. Competition between dilatant strengthening and thermal pressurization may control whether slip is slow or fast. We model SSE with 2‐D elasticity, rate‐state friction, and a dilatancy law where porosity evolves toward steady state ss over distance d c and ss = 0 + ln( v / v 0 ); v is slip speed. We consider two diffusion models. Membrane diffusion (MD) is approximated by −( p − p ∞ )/ t f where p and p ∞ are shear zone and remote pore pressure and t f is a characteristic diffusion time. Homogeneous diffusion (HD) accurately models fault‐normal flow with diffusivity c hyd . For MD, linearized analysis defines a boundary ℰ = 1 − a / b between slow and fast slip, where ℰ ≡ f 0 ε/ βb ( σ − p ∞ ), f 0 , a , and b are friction parameters and β is compressibility. When ℰ < 1 − a / b slip accelerates to instability for sufficiently large faults, whereas for ℰ > 1 − a / b slip speeds remain quasi‐static. For HD, E p ≡ ε h /( β ( σ − p ∞ ) ) defines dilatancy efficiency, where h is shear zone thickness and v ∞ is plate velocity. SSE are favored by large ε h and low effective stress. The ratio E p to thermal pressurization efficiency scales with 1/( σ − p ∞ ), so high p ∞ favors SSE, consistent with seismic observations. For E p ∼ 10 −3 transient slip rates, repeat times, average slip, and stress drops are comparable to field observations. Model updip propagation speeds are comparable to those observed along‐strike. Many simulations exhibit slow phases driven by steady downdip slip and faster phases that relax the accumulated stress. Model SSE accommodate only a fraction of plate motion; the remaining deficit must be accommodated during coseismic or postseismic slip.