Data products for research paper titled: "Assessing occurrence patterns of shallow Hikurangi slow slip events using renewal processes"
收藏资源简介:
Content Figures_SSEcatalog.fig contains two files: SSEnum_**.fig: Figure of east component GPS time series showing start and end time of each SSE, labelled by event number. SSEnum_**_map.fig: Figure showing location of GPS stations that detected a given SSE, labelled by event number. files_SSEcatalog.zip contains three files: SSEcatalog.mat: Structure containing the start and end times of all SSEs identified and the stations that detected each event. StationID.mat: Structure containing the start and end time of all SSEs identified at the StationID. SSEcatalog_simplified_refmodel.Rda: R-data frame containing the simplified SSE catalog information shown in Figure 2a. Note that start and end times of SSEs are given as dates (i.e., yyyy-mm-dd) without including time. waveletDetail_R1.zip: Contains the stacked wavelet details for R1, as in Figures S2 to S4. waveletDetail_R2.zip: Contains the stacked wavelet details for R2, as in Figures S5 to S7. Scripts.zip contains 13 scripts: denslines.R: Function used to plot posterior distributions. fit_renewal_process_case1.R: R-script to model SSE occurrence as a marked renewal process for Case 1. fit_renewal_process_case2.R: R-script to model SSE occurrence as a marked renewal process for Case 2. plot_posterior_distributions.R: R-script used to plot posterior distributions obtained from model results (Figure 3). ovlp.chain.R: Function that calculates the overlap between posterior distributions. hpd.interval.R: Function that calculates the highest posterior density (HPD) interval for each model parameter. forecast_after_Kaikoura.R: R-script to forecast SSE occurrence time after Kaikoura (Figure S12). plot_wavelet_details.m: Script to plot the wavelet details from Figures S2 to S7. residual_analysis.R: R-script to carry out the residual analysis shown in Figures S10 and S11. waic.R: Script to calculate the Watanabe-Akaike Information Criterion. get_tocc_Kaik.R: Function used to get occurrence time of SSEs after Kaikoura. get_tocc.R: Function used to get ocurrence time of SSEs.
### 数据集内容说明 1. **Figures_SSEcatalog.fig** 包含两个文件: - `SSEnum_**.fig`:展示各慢滑移事件(Slow Slip Event, SSE)东向分量GPS时间序列的图表,标注了每个SSE的起始与结束时间,并以事件编号进行标记。 - `SSEnum_**_map.fig`:展示了探测到指定SSE的GPS测站的空间分布,同样以事件编号进行标记。 2. **files_SSEcatalog.zip** 包含三个文件: - `SSEcatalog.mat`:存储了所有已识别SSE的起始与结束时间,以及每个事件对应的探测测站信息的结构体数据。 - `StationID.mat`:存储了各测站(StationID)所识别到的全部SSE的起始与结束时间的结构体数据。 - `SSEcatalog_simplified_refmodel.Rda`:R语言数据框文件,包含了图2a中展示的简化版SSE目录信息。 > 备注:SSE的起始与结束时间以日期格式(即`yyyy-mm-dd`)给出,未包含具体时刻。 3. **waveletDetail_R1.zip**:包含区域1(Region 1, R1)的叠加小波细节数据,对应图S2至图S4。 4. **waveletDetail_R2.zip**:包含区域2(Region 2, R2)的叠加小波细节数据,对应图S5至图S7。 5. **Scripts.zip** 包含13个脚本文件: - `denslines.R`:用于绘制后验分布的函数脚本。 - `fit_renewal_process_case1.R`:用于将SSE发生过程建模为标记更新过程(marked renewal process)的R脚本,适用于案例1。 - `fit_renewal_process_case2.R`:用于将SSE发生过程建模为标记更新过程的R脚本,适用于案例2。 - `plot_posterior_distributions.R`:用于绘制模型结果得到的后验分布的R脚本,对应图3。 - `ovlp.chain.R`:用于计算后验分布之间重叠程度的函数脚本。 - `hpd.interval.R`:用于计算每个模型参数的最高后验密度(Highest Posterior Density, HPD)区间的函数脚本。 - `forecast_after_Kaikoura.R`:用于预测凯库拉(Kaikoura)地震后SSE发生时间的R脚本,对应图S12。 - `plot_wavelet_details.m`:用于绘制图S2至S7中的小波细节的MATLAB脚本。 - `residual_analysis.R`:用于开展残差分析的R脚本,对应图S10与S11。 - `waic.R`:用于计算渡边赤池信息准则(Watanabe-Akaike Information Criterion, WAIC)的脚本。 - `get_tocc_Kaik.R`:用于获取凯库拉地震后SSE发生时间的函数脚本。 - `get_tocc.R`:用于获取SSE发生时间的函数脚本。



