Keyboard Model of Seismic Cycle of Great Earthquakes in Subduction Zones: Simulation Results and Further Generalization

Catastrophic megaearthquakes (M > 8) occurring in the subduction zones are among the most devastating hazards on the planet. In this paper we discuss the seismic cycles of the megathrust earthquakes and propose a blockwise geomechanical model explaining certain features of the stress-deformation cycle revealed in recent decades from seismological and satellite geodesy (GNSS) observations. Starting with an overview of the so-called keyboard model of the seismic cycle by L. Lobkovsky, we outline mathematical formalism describing the motion of seismogenic block system assuming viscous rheology beneath and between the neighboring elastic blocks sitting on top of the subducting slab. By summarizing the GNSS-based evidence from our previous studies concerning the transient motions associated with the 2006-2007 Simushir earthquakes, 2010 Maule earthquake, and 2011 Tohoku earthquake, we demonstrate that those data support the keyboard model and reveal specific effect of the postseismic oceanward motion. However, since the seismogenic blocks in subduction systems are mostly located offshore, the direct analysis of GNSS-measured displacements and velocities is hardly possible in terms of the original keyboard model. Hence, the generalized two-segment keyboard model is introduced, containing both frontal offshore blocks and rear onshore blocks, which allows for direct interpretation of the onshore-collected GNSS data. We present a numerical computation scheme and a series of simulated data, which exhibits the consistency with measured motions and enables estimating the seismic cycle characteristics, important for the long-term earthquake forecasting.