The role of continental lithospheric thermal structure in the evolution of orogenic systems: application to the Himalayan–Tibetan collision zone

Continental collision is a crucial process in plate tectonics. However, in terms of the evolution and the controlling parameters of its lateral heterogeneity, our understanding of the tectonic complexities at such a convergent plate boundary remains largely unclear. In this study, we conducted a series of two-dimensional numerical experiments to investigate how continental lithospheric thermal structure influences the development of lateral heterogeneity along the continental collision zone. The following two end-members were achieved. First, continuous subduction mode, which prevails when the model has a cold procontinental Moho temperature (≤450 ∘C). In this case, a narrow collision orogen develops, and the subducting angle steepens with the increasing retrocontinental Moho temperature. Second, continental subduction with a slab break-off, which generates a relative wide collision orogen and dominates when the model has a relatively hot procontinental Moho temperature (≥500 ∘C), especially when the Moho temperature ≥ 550 ∘C. Radioactive heat production is the second-order controlling parameter in varying the continental collision mode, while it prefers to enhance strain localization in the upper part of the continental lithosphere and promote the growth of shear zones there. By comparing the model results with geological observations, we suggest that the discrepant evolutionary paths from the continuous subduction underlying the Hindu Kush to the continental subduction after slab break-off beneath eastern Tibet may originate from the inherited lateral inhomogeneity of the Indian lithospheric thermal structure. Besides, the high content of crustal radioactive elements may be one of the most important factors that controls the formation of large thrust fault zones in the Himalayas.