Strain waves, earthquakes, slow earthquakes, and afterslip in the framework of the Frenkel-Kontorova model

The one-dimensional Frenkel-Kontorova (FK) model, well known from the theory of dislocations in crystal materials, is applied to the simulation of the process of nonelastic stress propagation along transform faults. Dynamic parameters of plate boundary earthquakes as well as slow earthquakes and afterslip are quantitatively described, including propagation velocity along the strike, plate boundary velocity during and after the strike, stress drop, displacement, extent of the rupture zone, and spatiotemporal distribution of stress and strain. The three fundamental speeds of plate movement, earthquake migration, and seismic waves are shown to be connected in framework of the continuum FK model. The magnitude of the strain wave velocity is a strong (almost exponential) function of accumulated stress or strain. It changes from a few km/s during earthquakes to a few dozen km per day, month, or year during afterslip and interearthquake periods. Results of the earthquake parameter calculation based on real data are in reasonable agreement with measured values. The distributions of aftershocks in this model are consistent with the Omori law for temporal distribution and a 1/r for the spatial distributions.