An introductory review of thermal structure of subduction zones: II. Numerical approach, validation, and comparison

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

Directories & Files

The subdirectories 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.

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.

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