Refine
Has Fulltext
- no (2)
Year of publication
- 2020 (2)
Document Type
- Article (2)
Language
- English (2)
Is part of the Bibliography
- yes (2)
Keywords
- earthquake modeling (1)
- giant earthquakes (1)
- multi-slab (1)
- numerical model (1)
- plume (1)
- singleslab (1)
- subduction (1)
- subduction zone (1)
Institute
- Institut für Geowissenschaften (2) (remove)
Initiation of subduction following the impingement of a hot buoyant mantle plume is one of the few scenarios that allow breaking the lithosphere and recycling a stagnant lid without requiring any preexisting weak zones. Here, we investigate factors controlling the number and shape of retreating subducting slabs formed by plume-lithosphere interaction. Using 3-D thermomechanical models we show that the deformation regime, which defines formation of single-slab or multi-slab subduction, depends on several parameters such as age of oceanic lithosphere, thickness of the crust and large-scale lithospheric extension rate. Our model results indicate that on present-day Earth multi-slab plume-induced subduction is initiated only if the oceanic lithosphere is relatively young (<30-40 Myr, but >10 Myr), and the crust has a typical thickness of 8 km. In turn, development of single-slab subduction is facilitated by older lithosphere and pre-imposed extensional stresses. In early Earth, plume-lithosphere interaction could have led to formation of either episodic short-lived circular subduction when the oceanic lithosphere was young or to multi-slab subduction when the lithosphere was old.
Giant earthquakes with magnitudes above 8.5 occur only in subduction zones. Despite the developments made in observing large subduction zone earthquakes with geophysical instruments, the factors controlling the maximum size of these earthquakes are still poorly understood. Previous studies have suggested the importance of slab shape, roughness of the plate interface contact, state of the strain in the upper plate, thickness of sediments filling the trenches, and subduction rate. Here, we present 2-D cross-scale numerical models of seismic cycles for subduction zones with various geometries, subduction channel friction configurations, and subduction rates. We found that low-angle subduction and thick sediments in the subduction channel are the necessary conditions for generating giant earthquakes, while the subduction rate has a negligible effect. We suggest that these key parameters determine the maximum magnitude of a subduction earthquake by controlling the seismogenic zone width and smoothness of the subduction interface. This interpretation supports previous studies that are based upon observations and scaling laws. Our modeling results also suggest that low static friction in the sediment-filled subduction channel results in neutral or moderate compressive deformation in the overriding plate for low-angle subduction zones hosting giant earthquakes. These modeling results agree well with observations for the largest earthquakes. Based on our models we predict maximum magnitudes of subduction earthquakes worldwide, demonstrating the fit to magnitudes of all giant earthquakes of the 20th and 21st centuries and good agreement with the predictions based on statistical analyses of observations.