Relation between rheological properties and the stress state in subducting slabs

Abstract

Abstract The distribution of different stress states in the subducting slab indicated by centroid moment tensor solutions for intra-slab earthquakes can help constrain the rheological properties of the slabs. A comparison of slabs in the western and eastern Pacific realms shows contrasting patterns in the stress states down to depths of ~ 350 km. The majority of slabs in the western Pacific show a pair of down-dip compression (DC) and down-dip tension (DT) domains in the upper and lower parts of the slab reflecting the effects of the slab unbending during progressive subduction. In contrast, slabs in the eastern Pacific show predominantly in-plane DT stress irrespective of slab geometry. Two-dimensional numerical simulations assuming constant slab thickness and viscosity indicate that the development of slabs with in-plane DT stress at depths of 100–300 km requires the slabs to be thin and have a low viscosity (1023 Pa s). Weak slabs bend easily and tend to fold when they encounter increased resistance to downward movement at the 660-km boundary. The associated DC stresses are not transmitted up the slab so negative buoyancy of the slab and DT stress dominates at intermediate depths for this type of slab. Most experimentally derived rheological parameters predict a high viscosity (> 1024 Pa s) for such slabs. However, two-dimensional numerical simulations using temperature- and pressure-dependent viscosity show that a relatively low activation energy (~ 110 kJ/mol) for diffusion creep is a possible explanation for the observed distribution of stresses in the slabs. Such low activation energies are compatible with recent experimental work on diffusion creep of polyphase mantle materials in which a low effective activation energy for creep results from a slow grain growth due to pinning effect of the secondary phase. The simulations provide a mechanical explanation for the observed dominantly DT stress state at 100–300 km depths for young slabs and paired DT and DC stress states at the same depth range for old slabs. Graphical Abstract

Keywords

• Stress in slab
• Intermediate-depth earthquake
• Slab rheology
• Numerical model
• The eastern and western Pacific subduction margins