Dynamic Emergence of Plate Motions and Great Megathrust Earthquakes Across Length and Time Scales

The slow motion of tectonic plates over thousands of kilometers is intermittently interrupted by great earthquakes with sudden slips localized near convergent plate boundaries. We developed a subduction model that self-consistently integrates buoyancy forces, diffusion and dislocation creep, and inter-plate friction. From the nonlinear dynamics emerge long-term plate motions that achieve velocities of approximate to 5 ${\approx} 5$ cm/year, effective viscosities of approximate to 1019 ${\approx} 1{0}<^>{19}$ Pa & sdot; $\cdot $s below plates, and sudden slips up to approximate to 10 ${\approx} 10$ m repeating every several hundred years. Along-strike resistance arising from long-wavelength variation of coseismic slip is naturally incorporated with a rupture length scale, L similar to $\tilde{L}$. Computations with L similar to similar to 103 $\tilde{L}\sim 1{0}<^>{3}$ km generate events with Mw approximate to 9 ${M}_{w}\approx 9$. When L similar to $\tilde{L}$ decreases, there is a commensurate decrease in the effective moment of rupture events. Predicted long-term plate velocities, mantle viscosities, cycles of stress loading and release, and rupture event size and magnitude all show good agreement with observations.