Abstract
Molnar and England (1990, https://doi.org/10.1029/JB095iB04p04833) introduced equations using a semianalytical approach that approximate the thermal structure of the forearc regions in subduction zones. A detailed new comparison with high-resolution finite element models shows that the original equations provide robust predictions and can be improved by a few modifications that follow from the theoretical derivation. The updated approximate equations are shown to be quite accurate for a straight-dipping slab that is warmed by heat flowing from its base and by shear heating at its top. The approximation of radiogenic heating in the crust of the overriding plate is less accurate but the overall effect of this heating mode is small. It is shown that the previous and updated approximate equations become increasingly inaccurate with decreasing thermal parameter and increasing variability of slab dip. It is also shown that the approximate equations cannot be extrapolated accurately past the brittle-ductile transition. Conclusions in a recent paper (Kohn et al., 2018, https://doi.org/10.1073/pnas.1809962115) that modest amount of shear heating can explain the thermal conditions of past subduction from the exhumed metamorphic rock record are invalid due to a number of compounding errors in the application of the Molnar and England (1990, https://doi.org/10.1029/JB095iB04p04833) equations past the brittle-ductile transition. The use of the improved approximate equations is highly recommended provided their limitations are taken into account. For subduction zones with variable dip and/or low thermal parameter finite element modeling is recommended.
Original language | English (US) |
---|---|
Pages (from-to) | 3268-3288 |
Number of pages | 21 |
Journal | Geochemistry, Geophysics, Geosystems |
Volume | 20 |
Issue number | 7 |
DOIs | |
State | Published - 2019 |
Bibliographical note
Funding Information:Many of the FEM models presented here can be reproduced with the open-source software of the FEniCS project using the distribution discussed in supporting information S1. We thank Scott King and an anonymous reviewer for their constructive reviews. We thank Rob Harris, Kelin Wang, Brad Hacker, and Xiang Gao for discussions. We thank Cian Wilson, Sarah Penniston-Dorland, Philippe Agard, and Peter Molnar for their comments on an earlier version of this manuscript. The research presented here was partially funded by the National Science Foundation through Grant EAR-1347262 to G.A. and EAR-1850634 to P.v.K. I.W. acknowledges support from the University of Minnesota through startup funds.
Publisher Copyright:
©2019. American Geophysical Union. All Rights Reserved.