Dynamic earthquake sequence simulation with an SBIEM accounting for interseismic poroelastic rebound

Abstract

AbstractAfterslip inside a coseismic slip patch is rarely observed, though some previous studies suggest that it is driven by poroelastic rebound (PER). These studies assume constant frictional strength, whereas time-dependent strengthening (healing) of a fault is expected from laboratory experiments, which provide a basis for a rate- and state-dependent friction law (RSF). In this study, quasistatic poroelasticity (PE) was implemented in a dynamic earthquake sequence simulation using a spectral boundary integral equation method, and the effect of PER on the behavior of a fault governed by RSF was examined. Spatio-temporal convolution for PE would significantly affect the resolution of the numerical simulation affordable. This problem has been resolved by numerical approximation of the time dependency of Green’s function of PE in the wavenumber domain, definition of memory variables, and reformulation of the temporal convolution into ordinary differential equations of them. In the novel method, the additional numerical costs due to PE are negligible. A planar fault with a rate-weakening patch embedded in the rate-strengthening region was simulated. Because it is the healing of the fault that competes against PER, both the aging law and slip law were examined, which have different characteristics in the evolution of the fault strength. The simulation results indicate that PER causes postseismic loading to the patch, but the healing efficiently suppresses afterslip not only for the aging law, but also for the slip law. When cases with different friction laws are compared, the healing is more significant for the aging law, which has log-$$t$$ t strengthening at a limit of $$V\to 0$$ V → 0 . However, the effect of PER on the slip rate is minor for the slip law. The slip law yields additional healing if the fault is accelerated by loading owing to PER. The simulation results are consistent with the absence of afterslip within the coseismic slip patches in the observations. Graphical Abstract