Data and Results for Rupture Dynamics of the 2025 Mw 7.7 Myanmar Earthquake

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** ([email protected]) - **Zhenguo Zhang** ([email protected]) ## 📁 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.