This repository contains data and code used in the seismic and petrological analysis described in our study.

READMEThis repository contains data and code used in the seismic and petrological analysis described in our study.1. Receiver Function (RF) WaveformsRaw receiver functions (RFs) were obtained using Gaussian low-pass filters with different width parameters (a). The filtered RFs are organized in the following directories:• ./RFwaveforms/Deconv-1.0 – RFs with Gaussian filter parameter a = 1.0• ./RFwaveforms/Deconv-2.0 – RFs with Gaussian filter parameter a = 2.0• ./RFwaveforms/Deconv-3.0 – RFs with Gaussian filter parameter a = 3.0• ./RFwaveforms/Sta.info – Station location information2. Trans-Dimensional Bayesian Joint Inversion CodeThe joint inversion was performed using a trans-dimensional Markov Chain Monte Carlo (McMC) Bayesian framework.Code is provided in: ./McMC-jointinvPlease check the Makefile for installing.3. Obtained Velocity ModelsFinal Vs and Vp/Vs models obtained from the joint inversion are stored in: ./VelocityModels/${STA}/${STA}-${backazimuth}Each folder corresponds to a specific station and back-azimuth combination.4. Petrological ModelingPetrological modeling was conducted using a modified version of the ElasticC code, originally available at: https://github.com/michpaulatto/ElasticC.The code has been modified to support three-phase modeling (solid + melt + fluid).本仓库包含本研究中地震学与岩石学分析所用的数据与代码。 1. 接收函数（Receiver Function, RF）波形 原始接收函数（RF）通过配置不同宽度参数（a）的高斯低通滤波器获取。滤波后的RF按如下目录结构组织： • ./RFwaveforms/Deconv-1.0：高斯滤波参数a=1.0的接收函数 • ./RFwaveforms/Deconv-2.0：高斯滤波参数a=2.0的接收函数 • ./RFwaveforms/Deconv-3.0：高斯滤波参数a=3.0的接收函数 • ./RFwaveforms/Sta.info：台站位置信息 2. 跨维贝叶斯联合反演代码 本次联合反演采用跨维马尔可夫链蒙特卡洛（Trans-dimensional Markov Chain Monte Carlo, McMC）贝叶斯框架实现。相关代码存放于./McMC-jointinv目录，请查阅该目录下的Makefile完成安装。 3. 反演得到的速度模型 联合反演得到的最终横波速度（Vs）与纵波与横波速度比（Vp/Vs）模型存储于./VelocityModels/${STA}/${STA}-${backazimuth}路径下，每个文件夹对应一组特定的台站与反方位角组合。 4. 岩石学建模 岩石学建模采用修改版的ElasticC代码完成，该代码的原始版本可从https://github.com/michpaulatto/ElasticC获取。修改后的代码支持三相建模（固态+熔体+流体）。