DOI: 10.5281/zenodo.15117293
title: Supplementary material for "Dynamic restrengthening and fault heterogeneity explain megathrust earthquake complexity"

description: 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.