Numerical study on multi-scaling earthquake rupture

A new numerical scheme using a renormalization and a 3D boundary integral equation method is proposed to simulate a multi-scaling dynamic rupture of earthquakes: How a small earthquake grows up to a large one in spatially heterogeneous field of critical slip-weakening distance D-c ( fracture energy G(c))? We examine the case where D-c grows according to a hypocentral distance L (D-c proportional to L-beta). When beta = 1, we succeed to show numerically that a rupture propagates at a constant rupture speed in uniform initial stress field. This result still keeps the scaling relation of Gc and Dc inferred for earthquake size, however no scale-dependent initial process is required. The break of the proportional relation changes rupture speed as well as slip velocity to keep the energy balance. The rupture is accelerated up to a speed even faster than the shear wave velocity (beta < 1) or naturally arrested (β > 1).