The distribution of earthquakes with depth and stress in subducting slabs

The global variation of Benioff zone seismicity with depth and the orientation of stress axes of deep and intermediate earthquakes is explained using numerical models of subducting slabs. Models that match the seismicity and stress require a barrier to flow at the 670 km seismic discontinuity. The barrier may be a viscosity increase of at least an order of magnitude or a chemical discontinuity. Instantaneous flow is subparallel to the slabs for models with a viscosity increase but contorted for models with a chemical barrier. Log N (number of earthquakes) decreases linearly to 250-300 km depth and increases thereafter. Stress magnitude in the models shows the same pattern, in accord with experiments showing N proportional to e(k-sigma), with k a constant and sigma the stress magnitude. The models predict downdip compression in the slabs at depths below 300-400 km, as observed for earthquake stress axes.