3-D Simulation of Earthquake Generation Cycles and Evolution of Fault Constitutitve Properties

The earthquake generation cycle consists of tectonic loading, quasi-static rupture nucleation, dynamic rupture propagation and stop, and subsequent stress redistribution and fault restrengthening. From a macroscopic point of view, the entire process of earthquake generation cycles should be consistently described by a coupled nonlinear system of a slip-response function, a fault constitutive law and a driving force. On the basis of such a general idea, we constructed a realistic 3-D simulation model for earthquake generation cycles at a transcurrent plate boundary by combining the viscoelastic slip-response function derived for a two-layered elastic-viscoelastic structure model, the slipand time-dependent fault constitutive law that has an inherent mechanism of fault restrengthening, and the steady relative plate motion as a driving force into a single closed system. With this model we numerically simulated the earthquake generation cycles repeated in a seismogenic region on a plate interface, and examined space-time changes in shear stress, slip deficits and fault constitutive properties during one complete cycle in detail. The occurrence of unstable dynamic slip brings about decrease both in fault strength and shear stress to a constant residual level. After the arrest of dynamic slip, the breakdown strength drop Drp of fault is restored rapidly and the process of stress accumulation resumes in the seismogenic region. On the other hand, the restoration of the critical weakening displacement Dc proceeds gradually with time through the interseismic period. The restoration of Dc can be regarded as the macroscopic manifestation of the microscopic recovery process of fractal fault surface structure. Through numerical simulation with a multi-segmented fault model, we examined the effects of viscoelastic faultto- fault interaction. The effect of transient viscoelastic stress transfer through the asthenosphere is significant as well as the direct effect of elastic stress transfer, and it possibly explains the time lag of the sequential occurrence of large events along a plate boundary.