Skip to content

Research at St Andrews

Thermo-mechanical controls of flat subduction: insights from numerical modeling

Research output: Contribution to journalArticle


Pengpeng Huangfu, Yuejun Wang, Peter A. Cawood, Zhong-Hai Li, Weiming Fan, Taras V. Gerya

School/Research organisations


Numerical experiments are used to investigate the thermo- mechanical controls for inducing flat subduction and why flat subduction is rare relative to normal/steep subduction. Our modeling results demonstrate that flat subduction is an end-member of a steady state subduction geometry and is characterized by a curved slab with a nearly-horizontal slab section. Intermediate cases between normal/steep and flat subduction appear to be transient in origin and evolve toward one of the stable end-members. Physical parameters inducing flat subduction can be classified into four categories: buoyancy of the subducting oceanic lithosphere (e.g., slab age, oceanic crustal thickness), viscous coupling between the overriding and downgoing plates (e.g., initial subduction angle), external kinematic conditions, and rheological properties of the subduction zone. On the basis of parameter sensitivity tests and the main characteristics of present-day flat subduction zones, positive buoyancy from either the young slab or the thickened oceanic crust are considered the primary controlling parameter. Our results show that the possibility of flat subduction is directly proportional to oceanic crustal thickness and inversely proportional to the slab age. Furthermore, oceanic crust must be thicker than 8 km to induce flat subduction, when the slab is older than 30 Ma with an initial subduction angle of ≥ 20°, and without absolute trenchward motion of the overriding plate. The lower the initial subduction angle or the thicker the overriding continental lithosphere, the more likelihood for flat subduction. The initial subduction angle is more influential for the development of flat subduction than the overriding lithospheric thickness, and a thick overriding lithosphere induces flat subduction only under the condition of an initial subduction angle of ≤ 25°, with a slab age of ≥ 30 Ma and without absolute trenchward motion of the overriding plate. However, when the initial subduction angle is increased to > 25°, no flat subduction is predicted. All the parameters are evaluated within the constraints of a mechanical framework in which the slab geometry is regarded as a result of a balance between the gravitational and hydrodynamic torque. Any factor that can sufficiently reduce gravitational torque or increase hydrodynamic torque will exert a strong effect on flat subduction development. Our results are consistent with the observations of modern flat subduction zones on Earth.


Original languageEnglish
Pages (from-to)170-183
Number of pages14
JournalGondwana Research
Early online date29 Sep 2016
Publication statusPublished - Dec 2016

    Research areas

  • Flat subduction, Numerical modeling, Slab buoyancy, Viscous interplate coupling, External kinematic conditions

Discover related content
Find related publications, people, projects and more using interactive charts.

View graph of relations

Related by author

  1. Neoarchean and Paleoproterozoic K-rich granites in the Phan Si Pan Complex, north Vietnam: constraints on the early crustal evolution of the Yangtze Block

    Zhao, T., Cawood, P. A., Wang, K., Zi, J-W., Feng, Q., Nguyen, Q. M. & Tran, D. M., 15 Sep 2019, In : Precambrian Research. 332, 105395.

    Research output: Contribution to journalArticle

  2. Jiangnan Orogen, South China: a ~970–820 Ma Rodinia margin accretionary belt

    Yao, J., Cawood, P. A., Shu, L. & Zhao, G., 3 Jun 2019, In : Earth Science Reviews. In press

    Research output: Contribution to journalArticle

Related by journal

  1. Evolution of the Mozambique Belt in Malawi constrained by granitoid U-Pb, Sm-Nd and Lu-Hf isotopic data

    Manda, B. W. C., Cawood, P. A., Spencer, C. J., Prave, T., Robinson, R. & Roberts, N. M. W., 29 Nov 2018, In : Gondwana Research. In press

    Research output: Contribution to journalArticle

  2. The crustal architecture of Myanmar imaged through zircon U-Pb, Lu-Hf and O isotopes: tectonic and metallogenic implications

    Gardiner, N. J., Searle, M. P., Morley, C. K., Robb, L. J., Whitehouse, M. J., Roberts, N. M. W., Kirkland, C. L. & Spencer, C. J., Oct 2018, In : Gondwana Research. 62, p. 27-60 34 p.

    Research output: Contribution to journalReview article

  3. Mapping temporal and spatial patterns of zircon U-Pb disturbance: a Yilgarn Craton case study

    Kirkland, C. L., Abello, F., Danišík, M., Gardiner, N. J. & Spencer, C., Dec 2017, In : Gondwana Research. 52, p. 39-47 9 p.

    Research output: Contribution to journalArticle

  4. From convergent plate margin to arc-continent collision: formation of the Kenting Mélange, southern Taiwan

    Zhang, X., Cawood, P. A., Huang, C-Y., Wang, Y., Yan, Y., Santosh, M., Chen, W. & Yu, M., Oct 2016, In : Gondwana Research. 38, p. 171-182 12 p.

    Research output: Contribution to journalArticle

ID: 246389865