Decarbonation of subducting slabs

Insight from petrological–thermomechanical modeling

C. M. Gonzalez, W. Gorczyk, T. V. Gerya

Research output: Contribution to journalArticleResearchpeer-review

11 Citations (Scopus)

Abstract

Subduction of heterogeneous lithologies (sediments and altered basalts) carries a mixture of volatile components (H2O ± CO2) into the mantle, which are later mobilized during episodes of devolatilization and flux melting. Several petrologic and thermodynamic studies investigated CO2 decarbonation to better understand carbon cycling at convergent margins. A paradox arose when investigations showed little to no decarbonation along present day subduction geotherms at subarc depths despite field based observations. Sediment diapirism is invoked as one of several methods for carbon transfer from the subducting slab. We employ high-resolution 2D petrological–thermomechanical modeling to elucidate the role subduction dynamics has with respect to slab decarbonation and the sediment diapirism hypothesis. Our thermodynamic database is modified to account for H2O–CO2 binary fluids via the following lithologies: GLOSS average sediments (H2O: 7.29 wt.% & CO2: 3.01 wt.%), carbonated altered basalts (H2O: 2.63 wt.% & CO2: 2.90 wt.%), and carbonated peridotites (H2O: 1.98 wt.% & CO2: 1.50 wt.%). We include a CO2 solubility P–x[H2O wt.%] parameterization for sediment melts. We parameterize our model by varying two components: slab age (20, 40, 60, 80 Ma) and convergence velocity (1, 2, 3, 4, 5, 6 cm year− 1). 59 numerical models were run and show excellent agreement with the original code base. Three geodynamic regimes showed significant decarbonation. 1) Sedimentary diapirism acts as an efficient physical mechanism for CO2 removal from the slab as it advects into the hotter mantle wedge. 2) If subduction rates are slow, frictional coupling between the subducting and overriding plate occurs. Mafic crust is mechanically incorporated into a section of the lower crust and undergoes decarbonation. 3) During extension and slab rollback, interaction between hot asthenosphere and sediments at shallow depths result in a small window (~ 12.5 Ma) of high integrated CO2 fluxes (205 kg m− 3 Ma− 1).

Original languageEnglish
Pages (from-to)314-332
Number of pages19
JournalGondwana Research
Volume36
DOIs
Publication statusPublished - 1 Aug 2016
Externally publishedYes

Keywords

  • Decarbonation
  • HO–CO fluids
  • Melting
  • Numerical modeling
  • Subduction

Cite this

@article{c6638af67ddf4687beae30f9e68b8a10,
title = "Decarbonation of subducting slabs: Insight from petrological–thermomechanical modeling",
abstract = "Subduction of heterogeneous lithologies (sediments and altered basalts) carries a mixture of volatile components (H2O ± CO2) into the mantle, which are later mobilized during episodes of devolatilization and flux melting. Several petrologic and thermodynamic studies investigated CO2 decarbonation to better understand carbon cycling at convergent margins. A paradox arose when investigations showed little to no decarbonation along present day subduction geotherms at subarc depths despite field based observations. Sediment diapirism is invoked as one of several methods for carbon transfer from the subducting slab. We employ high-resolution 2D petrological–thermomechanical modeling to elucidate the role subduction dynamics has with respect to slab decarbonation and the sediment diapirism hypothesis. Our thermodynamic database is modified to account for H2O–CO2 binary fluids via the following lithologies: GLOSS average sediments (H2O: 7.29 wt.{\%} & CO2: 3.01 wt.{\%}), carbonated altered basalts (H2O: 2.63 wt.{\%} & CO2: 2.90 wt.{\%}), and carbonated peridotites (H2O: 1.98 wt.{\%} & CO2: 1.50 wt.{\%}). We include a CO2 solubility P–x[H2O wt.{\%}] parameterization for sediment melts. We parameterize our model by varying two components: slab age (20, 40, 60, 80 Ma) and convergence velocity (1, 2, 3, 4, 5, 6 cm year− 1). 59 numerical models were run and show excellent agreement with the original code base. Three geodynamic regimes showed significant decarbonation. 1) Sedimentary diapirism acts as an efficient physical mechanism for CO2 removal from the slab as it advects into the hotter mantle wedge. 2) If subduction rates are slow, frictional coupling between the subducting and overriding plate occurs. Mafic crust is mechanically incorporated into a section of the lower crust and undergoes decarbonation. 3) During extension and slab rollback, interaction between hot asthenosphere and sediments at shallow depths result in a small window (~ 12.5 Ma) of high integrated CO2 fluxes (205 kg m− 3 Ma− 1).",
keywords = "Decarbonation, HO–CO fluids, Melting, Numerical modeling, Subduction",
author = "Gonzalez, {C. M.} and W. Gorczyk and Gerya, {T. V.}",
year = "2016",
month = "8",
day = "1",
doi = "10.1016/j.gr.2015.07.011",
language = "English",
volume = "36",
pages = "314--332",
journal = "Gondwana Research",
issn = "1342-937X",
publisher = "Elsevier",

}

Decarbonation of subducting slabs : Insight from petrological–thermomechanical modeling. / Gonzalez, C. M.; Gorczyk, W.; Gerya, T. V.

In: Gondwana Research, Vol. 36, 01.08.2016, p. 314-332.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

T1 - Decarbonation of subducting slabs

T2 - Insight from petrological–thermomechanical modeling

AU - Gonzalez, C. M.

AU - Gorczyk, W.

AU - Gerya, T. V.

PY - 2016/8/1

Y1 - 2016/8/1

N2 - Subduction of heterogeneous lithologies (sediments and altered basalts) carries a mixture of volatile components (H2O ± CO2) into the mantle, which are later mobilized during episodes of devolatilization and flux melting. Several petrologic and thermodynamic studies investigated CO2 decarbonation to better understand carbon cycling at convergent margins. A paradox arose when investigations showed little to no decarbonation along present day subduction geotherms at subarc depths despite field based observations. Sediment diapirism is invoked as one of several methods for carbon transfer from the subducting slab. We employ high-resolution 2D petrological–thermomechanical modeling to elucidate the role subduction dynamics has with respect to slab decarbonation and the sediment diapirism hypothesis. Our thermodynamic database is modified to account for H2O–CO2 binary fluids via the following lithologies: GLOSS average sediments (H2O: 7.29 wt.% & CO2: 3.01 wt.%), carbonated altered basalts (H2O: 2.63 wt.% & CO2: 2.90 wt.%), and carbonated peridotites (H2O: 1.98 wt.% & CO2: 1.50 wt.%). We include a CO2 solubility P–x[H2O wt.%] parameterization for sediment melts. We parameterize our model by varying two components: slab age (20, 40, 60, 80 Ma) and convergence velocity (1, 2, 3, 4, 5, 6 cm year− 1). 59 numerical models were run and show excellent agreement with the original code base. Three geodynamic regimes showed significant decarbonation. 1) Sedimentary diapirism acts as an efficient physical mechanism for CO2 removal from the slab as it advects into the hotter mantle wedge. 2) If subduction rates are slow, frictional coupling between the subducting and overriding plate occurs. Mafic crust is mechanically incorporated into a section of the lower crust and undergoes decarbonation. 3) During extension and slab rollback, interaction between hot asthenosphere and sediments at shallow depths result in a small window (~ 12.5 Ma) of high integrated CO2 fluxes (205 kg m− 3 Ma− 1).

AB - Subduction of heterogeneous lithologies (sediments and altered basalts) carries a mixture of volatile components (H2O ± CO2) into the mantle, which are later mobilized during episodes of devolatilization and flux melting. Several petrologic and thermodynamic studies investigated CO2 decarbonation to better understand carbon cycling at convergent margins. A paradox arose when investigations showed little to no decarbonation along present day subduction geotherms at subarc depths despite field based observations. Sediment diapirism is invoked as one of several methods for carbon transfer from the subducting slab. We employ high-resolution 2D petrological–thermomechanical modeling to elucidate the role subduction dynamics has with respect to slab decarbonation and the sediment diapirism hypothesis. Our thermodynamic database is modified to account for H2O–CO2 binary fluids via the following lithologies: GLOSS average sediments (H2O: 7.29 wt.% & CO2: 3.01 wt.%), carbonated altered basalts (H2O: 2.63 wt.% & CO2: 2.90 wt.%), and carbonated peridotites (H2O: 1.98 wt.% & CO2: 1.50 wt.%). We include a CO2 solubility P–x[H2O wt.%] parameterization for sediment melts. We parameterize our model by varying two components: slab age (20, 40, 60, 80 Ma) and convergence velocity (1, 2, 3, 4, 5, 6 cm year− 1). 59 numerical models were run and show excellent agreement with the original code base. Three geodynamic regimes showed significant decarbonation. 1) Sedimentary diapirism acts as an efficient physical mechanism for CO2 removal from the slab as it advects into the hotter mantle wedge. 2) If subduction rates are slow, frictional coupling between the subducting and overriding plate occurs. Mafic crust is mechanically incorporated into a section of the lower crust and undergoes decarbonation. 3) During extension and slab rollback, interaction between hot asthenosphere and sediments at shallow depths result in a small window (~ 12.5 Ma) of high integrated CO2 fluxes (205 kg m− 3 Ma− 1).

KW - Decarbonation

KW - HO–CO fluids

KW - Melting

KW - Numerical modeling

KW - Subduction

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

U2 - 10.1016/j.gr.2015.07.011

DO - 10.1016/j.gr.2015.07.011

M3 - Article

VL - 36

SP - 314

EP - 332

JO - Gondwana Research

JF - Gondwana Research

SN - 1342-937X

ER -