A viscoplastic shear-zone model for deep (15-50 km) slow-slip events at plate convergent margins

A key issue in understanding the physics of deep (15-50 km) slow-slip events (D-SSE) at plate convergent margins is how their initially unstable motion becomes stabilized. Here we address this issue by quantifying a rate-strengthening mechanism using a viscoplastic shear-zone model inspired by recent advances in field observations and laboratory experiments. The well-established segmentation of slip modes in the downdip direction of a subduction shear zone allows discretization of an interseismic forearc system into the (1) frontal segment bounded by an interseismically locked megathrust, (2) middle segment bounded by episodically locked and unlocked viscoplastic shear zone, and (3) interior segment that slips freely. The three segments are assumed to be linked laterally by two springs that tighten with time, and the increasing elastic stress due to spring tightening eventually leads to plastic failure and initial viscous shear. This simplification leads to seven key model parameters that dictate a wide range of mechanical behaviors of an idealized convergent margin. Specifically, the viscoplastic rheology requires the initially unstable sliding to be terminated nearly instantaneously at a characteristic velocity, which is followed by stable sliding (i.e., slow-slip). The characteristic velocity, which is on the order of <10−7 m/s for the convergent margins examined in this study, depends on the (1) effective coefficient of friction, (2) thickness, (3) depth, and (4) viscosity of the viscoplastic shear zone. As viscosity decreases exponentially with temperature, our model predicts faster slow-slip rates, shorter slow-slip durations, more frequent slow-slip occurrences, and larger slow-slip magnitudes at warmer convergent margins.