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

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

Research output: Contribution to journalArticleResearchpeer-review

Abstract

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 is considered as 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 torques. 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
Volume40
DOIs
Publication statusPublished - 1 Dec 2016
Externally publishedYes

Keywords

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

Cite this

Huangfu, Pengpeng ; Wang, Yuejun ; Cawood, Peter A. ; Li, Zhong Hai ; Fan, Weiming ; Gerya, Taras V. / Thermo-mechanical controls of flat subduction : Insights from numerical modeling. In: Gondwana Research. 2016 ; Vol. 40. pp. 170-183.
@article{81bffba87f684ea9944acbd667fd2f36,
title = "Thermo-mechanical controls of flat subduction: Insights from numerical modeling",
abstract = "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 is considered as 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 torques. 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.",
keywords = "External kinematic conditions, Flat subduction, Numerical modeling, Slab buoyancy, Viscous interplate coupling",
author = "Pengpeng Huangfu and Yuejun Wang and Cawood, {Peter A.} and Li, {Zhong Hai} and Weiming Fan and Gerya, {Taras V.}",
year = "2016",
month = "12",
day = "1",
doi = "10.1016/j.gr.2016.08.012",
language = "English",
volume = "40",
pages = "170--183",
journal = "Gondwana Research",
issn = "1342-937X",
publisher = "Elsevier",

}

Thermo-mechanical controls of flat subduction : Insights from numerical modeling. / Huangfu, Pengpeng; Wang, Yuejun; Cawood, Peter A.; Li, Zhong Hai; Fan, Weiming; Gerya, Taras V.

In: Gondwana Research, Vol. 40, 01.12.2016, p. 170-183.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

T1 - Thermo-mechanical controls of flat subduction

T2 - Insights from numerical modeling

AU - Huangfu, Pengpeng

AU - Wang, Yuejun

AU - Cawood, Peter A.

AU - Li, Zhong Hai

AU - Fan, Weiming

AU - Gerya, Taras V.

PY - 2016/12/1

Y1 - 2016/12/1

N2 - 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 is considered as 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 torques. 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.

AB - 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 is considered as 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 torques. 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.

KW - External kinematic conditions

KW - Flat subduction

KW - Numerical modeling

KW - Slab buoyancy

KW - Viscous interplate coupling

UR - http://www.scopus.com/inward/record.url?scp=84995948198&partnerID=8YFLogxK

U2 - 10.1016/j.gr.2016.08.012

DO - 10.1016/j.gr.2016.08.012

M3 - Article

VL - 40

SP - 170

EP - 183

JO - Gondwana Research

JF - Gondwana Research

SN - 1342-937X

ER -