Supplementary material for "Dynamic restrengthening and fault heterogeneity explain megathrust earthquake complexity"

This repository contains the supplementary files of the "Coseismic restrengthening and fault heterogeneity explain megathrust earthquake complexity" by Jeremy Wing Ching Wong (wcwong@ucsd.edu), Alice-Agnes Gabriel, and Wenyuan Fan. [Code availability] - The open-source SeisSol is available at [https://github.com/SeisSol/SeisSol]- We use the commit code '8cddb43' on the master branch to perform all dynamic rupture simulations presented in this study. - Instructions for downloading, installing, and running the code are available in the SeisSol documentation. - Downloading and compiling instructions are here. - Instructions for setting up and running simulations are here. - Quickstart containerized installations and introductory materials are provided in the Docker container and Jupyter Notebooks at the SeisSol Training Repository. - Example problems and model configuration files are provided at the SeisSol Example repository, many of which reproduce the SCEC 3D Dynamic Rupture benchmark problems. [File structure]Preferred dynamic rupture model (in zip file):├── Tohoku_FL103_plasticity.par : SeisSol parameter file for fast-velocity weakening rate-state friction├── Tohoku_fault.yaml : yaml file defining frictional properties├── Tohoku_initial_stress.yaml : yaml file defining initial stress conditions├── Tohoku_lithostress_gamma.yaml : yaml file defining lithostatic stress and pore-pressure level├── Tohoku_material_1D_C0_1MPa.yaml : yaml file for 1D velocity structure and plasticity ├── Tohoku_nucleation_stress_R0.8.yaml : yaml file for nucleation stress├── FL33_stress_change_kinematic_model.yaml : yaml file for stress heterogeneity based on imposed median slip distribution (Wong et al., 2024)├── M_stress_change_1km_1Hz_1000.nc : NC grid file for 3D stress change tensor based on imposed median slip distribution (Wong et al., 2024)├── mesh : mesh files │ ├── JIVSM_v7_15d_1000_1D_1Hz.puml.h5 │ ├── JIVSM_v7_15d_1000_1D_1Hz.xdmf │ └── JIVSM_v7_15d_1000_1D_1Hz.xml├── fault_receivers : on-fault and off-fault receivers │ ├── Fault_10km_node.txt │ ├── Tohoku_receivers.csv │ └── Tohoku_receivers.dat└── run_Seissol_general_o4.sh : SuperMUC job submission file Additional model comparison in the Discussion section (in zip file): - Fw_heterogeneity_model - heterogeneous fully-weakened friction with heterogeneity pattern derived from the median slip distribution and with homogeneous prestress - Sl0 - Dynamic rupture model with multiscale asperities in state-evolution distance - Kubota_FFM_DR_model - Model using stress heterogeneity pattern derived from Kubota et al., 2022 finite-fault slip model - Melgar_FFM_DR_model - Model using stress heterogeneity pattern derived from Melgar & Bock, 2016 finite-fault slip model - Yamazaki_FFM_DR_model - Model using stress heterogeneity pattern derived from Yamazaki et al., 2018 finite-fault slip model [Computational Resources] - SeisSol simulations were conducted on the supercomputing cluster SuperMUC-NG, hosted at the Leibniz Supercomputing Center (https://www.lrz.de/), Frontera, hosted at the Texas Advanced Computing Center (https://tacc.utexas.edu/), and Lenoardo, hosted at the CINECA datacenter in Bologna, Italy.