The Zagros orogen in the Iranian Plateau is a natural laboratory in which to investigate the tectonic evolution of the transition from oceanic subduction to continental collision. However, where slab detachment occurred and how slab detachment affects shallower continental subduction remain poorly understood. The formation mechanism of the post-collision magmatism is also unclear. Here, we construct a high-resolution Pn-wave attenuation model for the uppermost mantle beneath the Iranian Plateau and surrounding areas using a newly compiled dataset. Weak Pn attenuation delineates the Arabian Plate front near the Moho discontinuity, extending further toward the subduction direction in the northwestern and southeastern parts, possibly due to the up-bending and underplating of the Arabian lithosphere after the loss of slab pull. The correlations among the surface Miocene-Quaternary volcanism and strong Pn attenuation in the uppermost mantle suggest that asthenospheric materials escaped from slab windows and further feed the post-collision volcanoes.

Tethyan evolution is characterized by cyclical continent-transfer from Gondwana to the continents in the Northern Hemisphere, similar to a “one-way” train. Subduction has been viewed as the primary driver of transference. Therefore, it is crucial to understand the tectonic evolution of all past subduction zones that occurred along Eurasia’s southern margin. We studied the earliest known eclogite located at the Neo-Tethyan suture in the Iranian segment. A prograde-E-MORB-like eclogite reached a peak metamorphic condition of 2.2 GPa and 560°C, at 190 ± 11 Ma (1 rutile U-Pb ages), which constrains the youngest age for subduction initiation of the Neo-Tethyan slab. Combined with regional magmatic and structural data, the oldest age for Neo-Tethys subduction initiation is 210–192 Ma, which is younger than the Paleo-Tethyan closure time of 228–209 Ma. These data, used with previous numerical modeling, supports collision-induced subduction initiation. The collision-induced force, together with the Paleo-Tethyan subduction driven-mantle flow, is likely to have exploited weak inherited structures from earlier Neo-Tethyan rifting, resulting in a northward directed subduction zone along the southern margin of Central Iran Block.

Removal and thinning of cratonic lithosphere is believed to have occurred under different tectonic settings, for example, near subduction zones and above mantle plumes. Subduction-induced cratonic modification has been widely discussed; however, the mechanisms and dynamic processes of plume-induced lithospheric removal remain elusive and require further systematic investigation. In this study, we conduct a series of 2-D thermo-mechanical models to explore the dynamics of the removal and thinning of cratonic lithosphere due to the interaction between a mantle plume and a weak mid-lithosphere discontinuity (MLD) layer. Our modeling results suggest that the interaction between a mantle plume and weak MLD layer can lead to a large-scale removal of the cratonic lithosphere as long as the connection between the hot upwelling and weak MLD layer is satisfied. The presence of a vertical lithospheric weak zone and its closeness to the plume center play critical roles in creating a connection between the weak MLD and hot plume/asthenosphere. Furthermore, delamination of cratonic lithosphere is favored by a larger plume radius/volume, a higher plume temperature anomaly, and a lower viscosity of the MLD layer. A systematic comparison between subduction-induced and plume-induced lithospheric thinning patterns is further conducted. We summarize their significant differences on the origin and migration of melt generation, the water content in melts, and topographic evolution. The combination of numerical models and geological/geophysical observations indicates that mantle plume-MLD interaction may have played a crucial role in lithospheric removal beneath South Indian, South American and North Siberian Cratons.

Felsic continental crust is unique to the Earth in the solar system, but it still remains controversial regarding its formation, accretion and reworking. The plate tectonics theory has been significantly challenged in explaining the origin of continents as Archean continents rarely preserve hallmarks of plate tectonics. In contrast, growing evidence emerges to support mantle plume-derived oceanic plateau models as the models can reasonably explain the origin of bimodal volcanic assemblages and nearly coeval emplacement of tonalite-trondjhemite-granodiorite (TTG) rocks, presence of ~1600ºC komatiites and dominant dome structures, and lack of ultra-high-pressure rocks, paired metamorphic belts and ophiolites in Archean continents. Although plate tectonics seems to fail in explaining the origin of continents, it has been successfully applied to interpret the accretion or outgrowth of continents along subduction zones where new mafic crust is generated at the base of continental crust through partial melting of the mantle wedge with addition of H2O-dominant fluids from the subducted oceanic slabs, and partial melting of the juvenile mafic crust results in the formation of new felsic continental crust, leading to the outside accretion of continents. Subduction processes also cause the softening, thinning and recycling of continental lithosphere due to the vigorous infiltration of volatile-rich fluids and melts especially along weak layers or weak belts, leading to the widespread reworking and even destruction of continental lithosphere. Reworking of continents also occurs in continental interiors due to plume-lithosphere interactions, which, however, leads to much less degrees of lithospheric modification than subduction-induced craton destruction.

The subducting slab morphology in the mantle transition zone (MTZ) is strongly affected by the mantle viscosity and density variations at the 660-km discontinuity (D660). Besides the negative Clapeyron slope of phase transition and the viscosity increase, a possible thin weak layer at around D660 is proposed to play a key role in the slab stagnation, which is however not well constrained. In this study, a series of numerical models are systematically conducted, which reveal that a weak layer beneath D660 does not change the slab mode selection (penetration versus stagnation). However, it will contribute to longer slab flattening at the bottom of the MTZ, when slab sinking is strongly resisted by either the viscosity increase or a large Clapeyron slope at D660. The role of a weak layer on slab flattening is dependent on the lubrication effect that promotes sub-horizontal slab movement at the bottom of the MTZ.