Data and Results for Dynamics Rupture Simulation Co-constrained by Interseismic Locking Models and Regional Stress Field: A Case Study of the Luding Earthquake

This repository contains all datasets, input files, and results required to reproduce the dynamic rupture scenarios and strong ground motion of the 2022 Mw 6.6 Luding earthquake ## 📧 Contact For questions, please contact: - **Yuhao Gu** (guyh1997@mail.ustc.edu.cn) - **Zhengguo Zhang** (zhangzg@sustech.edu.cn) - **Xiaofeng Jia** (xjia@sustech.edu.cn) --- ## 📁 Repository Structure ### 📂 Folders #### `dynamic rupture scenarios` - `params.json` - Configuration parameters - `100.2_22.2.mod` - a case of velocity structure - `fault_coord.nc` - Fault geometry - `source.m` - Generate strong ground motion sources - `getmedia.py` - Interpolation to obtain velocity structure around faults for research ##### `Preferred dynamic rupture scenarios` ###### `Mod1_Co-constrained` - Preferred Model (Mod1) - `init_stress.nc` - Initial stress conditions ###### `Mod2_only_inter` - Model with interseismic locking model (Mod2) - `init_stress.nc` - Initial stress conditions ###### `Mod3_only_field` - Model with regional stress field (Mod3) - `init_stress.nc` - Initial stress conditions ##### `1.15 factor dynamic rupture scenarios` ###### `Mod1_Co-constrained` - Preferred Model (Mod1) - `init_stress.nc` - Initial stress conditions ###### `Mod2_only_inter` - Model with interseismic locking model (Mod2) - `init_stress.nc` - Initial stress conditions ###### `Mod3_only_field` - Model with regional stress field (Mod3) - `init_stress.nc` - Initial stress conditions #### `strong ground motion` - `params.json` - Configuration parameters --- ## 🔗 Related Resources ### Simulation Code - **GPU-accelerated CGFDM dynamic simulation code (`dynfd3dcuda`)**: [https://github.com/wenqiang-geophys/dynfd3dcuda](https://github.com/wenqiang-geophys/dynfd3dcuda) - **strong ground motion code (`CGFDM3D-EQR`)**: [https://github.com/wenqiangwangsustech](https://github.com/wenqiangwangsustech) ### Data Sources - **Xianshuihe fault trace data**: [http://www.cses.ac.cn/sjcp/ggmx/2021/135.shtml](http://www.cses.ac.cn/sjcp/ggmx/2021/135.shtml) - **Intensity distribution data**: [https://www.mem.gov.cn/xw/yjglbgzdt/202209/t20220911_422190.shtml](https://www.mem.gov.cn/xw/yjglbgzdt/202209/t20220911_422190.shtml) - **Stations instrument intensity data**: [https://seismisite.net/index_en.html](https://seismisite.net/index_en.html) - **S-wave velocity data**: [https://zenodo.org/records/11098135#:~:text=This%20high-resolution%20CSRM-1.0%20is%20constructed%20for%20the%20top,of%204435%20seismic%20stations%20during%201990%20and%202020](https://zenodo.org/records/11098135#:~:text=This%20high-resolution%20CSRM-1.0%20is%20constructed%20for%20the%20top,of%204435%20seismic%20stations%20during%201990%20and%202020) --- ## 📝 Note A detailed description of each file and folder is provided in this README.md file. Therefore, the installation and configuration of dynfd3dcuda and CGFDM3D-EQR are not repeated here. In practice, a new directory should be created under ../dynfd3dcuda/jobs/ as the job directory (e.g., ../dynfd3dcuda/jobs/temp). First, download the S-wave velocity data, convert it to the 100.2_22.2.mod format, and name it media02. Then, run getmedia.py to generate the medium file media3d_0.2.nc. Next, copy the init_stress.nc file corresponding to the target model into the same directory as fault_coord.nc (i.e., ../dynfd3dcuda/jobs/temp). Submit and run the dynfd3dcuda program. After the simulation completes, run source.m to obtain the strong ground motion source, move it to the /source/ directory of the CGFDM3D-EQR program, and finally submit and run CGFDM3D-EQR.