This repository contains all datasets, input files, and results required to reproduce the dynamic rupture scenarios of the 2025 Mw 7.7 Myanmar earthquake, including outputs from a comprehensive parameter grid search. For questions, please contact: - **Zhongqiu He** (hezq2021@mail.sustech.edu.cn) - **Zhenguo Zhang** (zhangzg@sustech.edu.cn) ## 📁 Repository Structure ### 📂 Folders #### `Mod1_preferred` - Preferred Model (Mod1) - `params.json` - Configuration parameters - `init_stress.nc` - Initial stress conditions - `mod_*_mrf.mat` - Modeled moment rate functions - `mod_*_ge_npw_syn.mat` - Modeled velocity waveforms at NPW station - `mod_*_surf_slip.txt` - Modeled fault-parallel on-fault slip #### `Mod2_no_barrier` - Model without Localized Barrier (Mod2) - Contains same file types as `Mod1_preferred` #### `Mod3_uniform_dc` - Uniform D_c Model (Mod3) - Contains same file types as `Mod1_preferred` #### `insar_pot` - Sentinel-1 azimuth pixel offset tracking (POT) data *Source: Wang, X., & Xu, X. (2025). Coseismic Surface Deformation of the 2025 Myanmar Mw 7.7 Earthquake Revealed by Sentinel-1 Pixel Offset Tracking and DInSAR Analyses. Zenodo. https://doi.org/10.5281/zenodo.15465249* #### `sm` - NPW velocity waveforms ### 📄 Files #### Core Model Files - `media1D.dat` - 1D velocity structure - `fault_coord.nc` - Fault geometry - `interp_sd_0510.mat` - Interpolated stress drop derived from USGS kinematic source model - `us7000pn9s_moment_rate_v4.txt` - Source time function from USGS - `fig_s01_mrf_comp.png` - Comparison of moment rate functions from different models #### Observation Data - `S2_parallel_displacement_20250401.txt` - Fault-parallel on-fault displacements from Sentinel-2 optical imagery - `Myanmar_Sagaing_fault.txt` - Fault rupture trace - `myanmar_aftershocks_tmd.txt` - 56-day aftershocks from Thai Meteorological Department (TMD) catalog #### Grid Search Results - `grid_inv_scaled_dc_all.txt` - Grid search results for slip-scaled D_c scenarios - **Columns**: `no.`, `C`, `alpha`, `S`, `Mw`, `M0`, `var_sm`, `var_surf`, `var_pot` - **Total misfit**: χ² = `var_sm` × 0.5 + `var_surf` × 0.25 + `var_pot` × 0.25 - `grid_inv_uniform_dc_all.txt` - Grid search results for uniform D_c scenarios - **Columns**: `no.`, `D_c`, `alpha`, `S`, `Mw`, `M0`, `var_sm`, `var_surf`, `var_pot` - **Total misfit**: χ² = `var_sm` × 0.5 + `var_surf` × 0.25 + `var_pot` × 0.25 ## 🔗 Related Resources ### Code - **GPU-accelerated CGFDM dynamic simulation code (`dynfd3dcuda`)**: [https://github.com/wenqiang-geophys/dynfd3dcuda](https://github.com/wenqiang-geophys/dynfd3dcuda) ### Data Sources - **USGS Finite Fault Model**: [https://earthquake.usgs.gov/earthquakes/eventpage/us7000pn9s/finite-fault](https://earthquake.usgs.gov/earthquakes/eventpage/us7000pn9s/finite-fault) - **USGS Emergency Public Data**: [https://earthquake.usgs.gov/static/Emergency%20Public%20Data/Burma/](https://earthquake.usgs.gov/static/Emergency%20Public%20Data/Burma/) - **TMD Earthquake Catalog**: [https://earthquake.tmd.go.th/](https://earthquake.tmd.go.th/) ## 📝 Note A detailed description of each file and folder is provided in this README.md file.