DOI: 10.5281/zenodo.18457982
title: Data and Results for Dynamics Rupture Simulation Co-constrained by Interseismic Locking Models and Regional Stress Field: A Case Study of the Luding Earthquake

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