Scripts to reproduce the results of "Did the short-term precursory accelerating phase precede the 2023 Mw 7.8 Kahramanmaraş earthquake?"

This repository contains the scripts to reproduce the results of Guo and Wang. Directories had to be zipped to be downloaded. Please start by unzipping them. If you wish to run the scripts and reproduce the results, you will first download the PRIDE PPP-AR software (https://github.com/PrideLab/PRIDE-PPPAR) to process RINEX observation files. 1. Data The data files obtained using single-day data processing and continuous processing strategies are converted into the same format as the data files used in the Bletery program. 1.1 TSfile(gapdata)_singleday: Data concatenation files obtained using the single-day processing strategy. 1.2 TSfile(gapdata)_continueday: Data concatenation files obtained using the continuous processing strategy. 1.3 TSfile(fulldata)_singleday: Missing values are interpolated by averaging the preceding and following data points for TSfile(gapdata)_singleday. 1.4 TSfile(fulldata)_continueday: Missing values are interpolated by averaging the preceding and following data points for TSfile(gapdata)_continueday. 1.5 station_single_200km: On the basis of interpolation, the influence of common mode error (CME) is removed for TSfile(fulldata)_singleday. 1.6 station_continue_200km: On the basis of interpolation, the influence of CME is removed for TSfile(fulldata)_continueday. 2. different location of hypocenter Investigating the impact of different source locations on the stacking results. 2.1 The folders 1-55 contain results from different source locations, including time series and computed dot product stack. 2.2 cal_stack2h.m is a MATLAB script used to compute the normalized stacking of the average values of the time series during the two-hour period prior to the earthquake. 2.3 LB.txt is the file that contains the hypocenter locations used for testing purposes. 3. different-distance The different-distance folder contains data used to assess the effect of CME on the results at various distances from the hypocenter. 3.1 serlect.bat is a batch file used to compute and extract stations within specified epicentral distance ranges. 3.2 The station_single_x00km folder contains the results of calculations after removing the CME at various distances. 4. Uniqueness_signals The Uniqueness_signals folder contains data utilized to assess whether the observed signals are attributable to noise. You need to change the absolute paths of all matlab files to your own!本仓库包含用于复现Guo与Wang研究成果的脚本。为便于下载，所有目录均已打包压缩，请先解压。若需运行脚本并复现研究结果，您需首先下载PRIDE PPP-AR软件（https://github.com/PrideLab/PRIDE-PPPAR）以处理RINEX观测文件。 1. 数据 采用单日数据处理与连续数据处理策略得到的数据文件，已转换为与Bletery程序所用数据文件一致的格式。 1.1 TSfile(gapdata)_singleday：采用单日处理策略得到的数据拼接文件。 1.2 TSfile(gapdata)_continueday：采用连续处理策略得到的数据拼接文件。 1.3 TSfile(fulldata)_singleday：针对TSfile(gapdata)_singleday，通过对前后数据点取平均的方式插值补全缺失值后得到的文件。 1.4 TSfile(fulldata)_continueday：针对TSfile(gapdata)_continueday，通过对前后数据点取平均的方式插值补全缺失值后得到的文件。 1.5 station_single_200km：在插值补全的基础上，针对TSfile(fulldata)_singleday去除共模误差（Common Mode Error, CME）的影响后得到的结果。 1.6 station_continue_200km：在插值补全的基础上，针对TSfile(fulldata)_continueday去除共模误差（Common Mode Error, CME）的影响后得到的结果。 2. 不同震源位置 探究不同震源位置对堆叠结果的影响。 2.1 编号1至55的文件夹包含不同震源位置对应的结果，包括时间序列与计算得到的点积堆叠结果。 2.2 cal_stack2h.m为MATLAB脚本，用于计算地震发生前两小时内时间序列平均值的归一化堆叠结果。 2.3 LB.txt为包含测试所用震源位置信息的文件。 3. 不同震中距 不同震中距文件夹包含用于评估共模误差对震源不同距离处结果影响的数据。 3.1 serlect.bat为批处理文件，用于计算并提取指定震中距范围内的测站数据。 3.2 station_single_x00km文件夹包含不同震中距下去除共模误差后的计算结果。 4. 信号唯一性 Uniqueness_signals文件夹包含用于评估观测信号是否由噪声引起的数据。 请将所有MATLAB文件的绝对路径修改为您本地的实际路径！