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 with a "dynamic" slab (Stokes equations solved in the slab portion of the domain as in van Keken & Wilson, PEPS, 2023 and `D80new` below) and performed 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.

`vankeken_wilson_peps_2023_TF_lowres_minimal.zip`: contains a set of low resolution vtu files from TerraFERMA with the final potential temperature output from steady state and time dependent simulations.  Intended for fast downloading (compared to `vankeken_wilson_peps_2023.zip`) and comparison.

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