Boso seismicity swarm propagation driven by slow slip stress change

The interactions between aseismic slip and seismicity is mostly studied during postseismic afterslip associated with the generation of aftershocks following large earthquakes. However because of the large stress perturbations produced by the coseismic rupture, it remains difficult to distinguish the contribution of the coseismic stress perturbation and the effects of stresses induced by the afterslip on the triggering of aftershocks. Studying seismicity triggered by slow slip events enables us to understand the direct effect of aseismic slip on the generation of the seismicity. While most subduction slow slip events are deep-seated, at the down-dip edge of seismogenic zones accompanied by tectonic tremors, some are also observed at shallower depths associated with seismic swarms. Among them are the well-studied Boso slow slip events located on the Sagami trough, between 10 and 20 km depth. Recorded every ~4 years since 1996, they are always accompanied by swarms of Mw 1 to 5 earthquakes on their northern and western flanks. Being located right beneath the Boso Peninsula coastline, the kinematics of these slow slip events is particularly well imaged by dense GNSS and tiltmeter networks. In this study, we model the time dependent aseismic slip of the 2018 Boso slow slip event, with the largest moment released of all Boso slow slip events, by inverting time step by time step the slip on the fault with joint GNSS and tiltmeter data. We do not impose arbitrary temporal smoothing in the inversion and find that the well constrained fault slip is first growing and then migrating southwestward with a migration speed of ~ 2 km/day. In order to model the interaction with the seismicity, we compute the Coulomb stress change due to the transient slip on receiver faults located in 9 cubes centered in the seismicity swarm and parallel to the subduction interface. From June 2nd to June 18th, the seismicity is migrating up-dip at a rate of 1 km/day. This migration period coincides with the peak slip rate and with Coulomb stress produced by the slow slip migrating updip together with the seismicity, indicating a causal relationship. Adopting a rate and state friction formalism to explain the nucleation of the seismicity, we finally investigate the ensemble of parameters, in particular the constitutive parameter that relates changes in stress to logarithmic changes of slip velocity, the effective normal stress and the tectonic stressing rates, that can explain the seismicity rate.