  We discover a remarkable correlation between the inter-tremor time interval
and the slenderness ratio of the overriding plate in subduction zones all over
the world. In order to understand this phenomenon better, we perform numerical
simulations of 3D deformation. The numerical buckling studies show that
critical load and slenderness ratio indeed have an inverse nonlinear relation
between them -- identical to the classical Euler's critical load relation, and
closely resemble the relationship observed between the inter-tremor time
interval and the slenderness ratio of the overriding plate. From the above
analysis, we conclude that the observed relation is the result of buckling of
the overriding continental plate. In addition to the above numerical analysis,
we analyze the surficial 3D spatio-temporal displacements of the overriding
plates in Cascadia and Alaska using 3-component GPS data. We find that these
deformations are consistent with the buckling of the overriding continental
crust. Based on these novel observations and guided by numerous existing
scientific observations and findings, we propose an Episodic Buckling and
Collapse model of subduction zones, wherein periodic geodetic changes and
tectonic-tremor activity, result from the episodic buckling of the overriding
continental crust and its rapid collapse on the subducting oceanic slab.
According to this model, geodetic measurements, previously inferred as slow
slip, are the surficial expressions of slowly-evolving buckling and rapid
collapse of the overriding plate, while tremor swarms result from the striking
of the collapsing overriding plate on the subducting slab (as opposed to
slipping or shearing).
