An introductory review of thermal structure of subduction zones: II & III

This repository contains the output from a global suite of subduction zone models used in: "An introductory review of the thermal structure of subduction zones: II—numerical approach and validation", Wilson & van Keken, PEPS, 2023 and "An introductory review of the thermal structure of subduction zones: III—Comparison between models and observations", van Keken & Wilson, PEPS, 2023.

Directories & Files

The subdirectories of `vankeken_wilson_peps_2023.zip` are:

• `D80`: a reproduction of the "D80" Syracuse et al., PEPI, 2010 suite of subduction zones using Sepran
• `D80new`: a new global suite based on Syracuse et al., PEPI, 2010 but using a lower mantle potential temperature and an error function boundary condition and performed using both Sepran and TerraFERMA.  Input files provided for TerraFERMA.
• `benchmark`: a new subduction benchmark performed using both Sepran and TerraFERMA.  Input files provided for TerraFERMA.
• `Figure_III_4_models`: the results used in Figure 4, Part III.

Files include:

• `Table2_Syracuse2010_updated.pdf` and `Table2_Syracuse2010_updated.xlsx`: updated parameter values for the  Syracuse et al., PEPI, 2010 global suite.
• `Table3_part_Syracuse2010_updated.pdf` and `Table3_part_Syracuse2010_updated.xlsx`: updated parameter values for the  Syracuse et al., PEPI, 2010 global suite.
• `Figure_comp_D80_WW09_all.pdf`: a comparison between the Sepran models presented in `D80` and `D80new` and those in Wada & Wang, G-cubed, 2009.
• `Figure_comp_D80new_TF-Sepran_all.pdf`: comparison between the TerraFERMA and Sepran models presented in `D80new`.
• `Figure_III_6_all_D80.pdf`: complete set of model results for Figure 6, Part III.
• `Subduction_parameters.pdf` and `Subduction_parameters.json`: a description of the parameters used in the global suite.
• `PEPS_part_II_preprint.pdf`: preprint of part II.
• `PEPS_part_III_preprint.pdf`: preprint of part III.

Running a simulation

Input files are provided for the TerraFERMA (the Transparent Finite ELement Rapid Model Assembler) simulations in this repository (`TF` subdirectories).  Running them requires a working installation of TerraFERMA.  If one is not available then consider using the docker image provided for this paper at:

https://github.com/users/cianwilson/packages/container/package/vankeken_wilson_peps_2023

This docker image contains a complete installation of TerraFERMA and its dependencies, PETSc, FEniCS and SPuD, within an Ubuntu 20.04LTS OS.  For a full description of TerraFERMA please refer to the webpage:

http://terraferma.github.io