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 ([e-mail]), 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.本仓库包含论文《同震强化与断层异质性解释大逆冲地震复杂性》（原英文标题：Coseismic restrengthening and fault heterogeneity explain megathrust earthquake complexity）的补充文件，作者为Jeremy Wing Ching Wong（电子邮箱：[e-mail]）、Alice-Agnes Gabriel以及Wenyuan Fan。 【代码可及性】 - 开源的SeisSol可从[https://github.com/SeisSol/SeisSol]获取——本研究中所有动态破裂模拟均基于master分支的commit代码'8cddb43'完成。 - 代码的下载、安装与运行说明详见SeisSol官方文档。 - 下载与编译说明详见此处。 - 模拟的设置与运行说明详见此处。 - SeisSol训练仓库提供了Docker容器与Jupyter Notebook形式的快速容器化安装指南及入门材料。 - SeisSol示例仓库提供了示例问题与模型配置文件，其中多数可复现SCEC三维动态破裂基准测试问题。 【文件结构】 首选动态破裂模型（压缩包内）： ├── Tohoku_FL103_plasticity.par：用于快速速度弱化速率状态摩擦的SeisSol参数文件 ├── Tohoku_fault.yaml：定义摩擦属性的YAML配置文件 ├── Tohoku_initial_stress.yaml：定义初始应力条件的YAML配置文件 ├── Tohoku_lithostress_gamma.yaml：定义岩石静应力与孔隙压力水平的YAML配置文件 ├── Tohoku_material_1D_C0_1MPa.yaml：用于一维速度结构与塑性行为的YAML配置文件 ├── Tohoku_nucleation_stress_R0.8.yaml：用于定义成核应力的YAML配置文件 ├── FL33_stress_change_kinematic_model.yaml：基于预设中位滑动分布（Wong等，2024）的应力异质性YAML配置文件 ├── M_stress_change_1km_1Hz_1000.nc：基于预设中位滑动分布（Wong等，2024）的三维应力变化张量NC网格文件 ├── mesh：网格文件目录 │ ├── JIVSM_v7_15d_1000_1D_1Hz.puml.h5 │ ├── JIVSM_v7_15d_1000_1D_1Hz.xdmf │ └── JIVSM_v7_15d_1000_1D_1Hz.xml ├── fault_receivers：断层内外监测点目录 │ ├── Fault_10km_node.txt │ ├── Tohoku_receivers.csv │ └── Tohoku_receivers.dat └── run_Seissol_general_o4.sh：SuperMUC集群作业提交脚本 讨论部分的额外模型对比（压缩包内）： - Fw_heterogeneity_model：基于中位滑动分布提取异质性模式、预设均匀预应力的非均匀全弱化摩擦模型 - Sl0：状态演化距离下包含多尺度凹凸体的动态破裂模型 - Kubota_FFM_DR_model：基于Kubota等，2022年有限断层滑动模型提取应力异质性模式的模型 - Melgar_FFM_DR_model：基于Melgar与Bock，2016年有限断层滑动模型提取应力异质性模式的模型 - Yamazaki_FFM_DR_model：基于Yamazaki等，2018年有限断层滑动模型提取应力异质性模式的模型 【计算资源】 SeisSol模拟运行于以下超级计算集群：莱布尼茨超级计算中心（LRZ）托管的SuperMUC-NG（https://www.lrz.de/）、德克萨斯高级计算中心（TACC）托管的Frontera（https://tacc.utexas.edu/），以及意大利博洛尼亚CINECA数据中心托管的Leonardo集群。