Unique composition and evolution histories of low velocity mantle domains: Data and analysis

Dataset and Jupyter notebook accompanying 'Unique composition and evolution histories of low velocity mantle domains'. 

Dataset includes:

• Present day properties of simulated mantle for simulations RCY, B=0.22, B=0.44, visc2, visc3, CMB2600, CMB2800, COMP, PRM, MER.
• Present day predicted seismic properties for simulations RCY, B=0.22, B=0.44, visc2, visc3, CMB2600, CMB2800, COMP, MER.
• Present day predicted seismic properties for simulation PRM assuming 'primordial' material to be i) basaltic oceanic crust ii) chondrite enriched basalt (CEB).
• Present day delta Vs for simulations RCY, B=0.22, B=0.44, visc2, visc3, CMB2600, CMB2800, COMP, PRM, MER, filtered using the resolution of seismic tomography model S40RTS.
• Properties of simulated mantle at 100 Myr intervals from 900 Ma - 100 Ma inclusive for simulation RCY. 
• P-T tables with predicted abundance of post-perovskite for different mantle lithologies (harzburgite, lherzolite and basalt - as defined in the paper).

Jupyter notebook `s-llvps.ipynb` contains code for identifying simulated large low-velocity provinces (S-LLVPs), extracting assoicated model properties and plotting results. Python module files terra_utils.py and ppv.py are also included and required by the code in the notebook. 

There are a number of pre-requisite packages that will need to be installed in order to run the Jupyter notebook, including terratools, a software package written specifically for reading and postprocessing outputs from TERRA simulations. Installation instructions can be found on the GitHub repository. 

Due to the TERRA code pre-dating open source licensing, we do not currently have permission to publicly share all aspects of the code. In code_pieces.F90 we include code snippets which were implemented for this study. 

Simulations were conducted using ARCHER2, the UK's national super-computing service. 

RCY.mp4 is a movie produced for simualtion RCY, visualising the evolution of temperature (right panels) and bulk composition (left panels). Hot iso-surface (red) drawn at +500 K and cold iso-surface (blue) drawn at -400 K, composition iso-surface drawn at C=0.6. Red and blue lines indicate overlying ridges / subduction zones taken from the plate motion reconstructions of Müller et al (2022).