A dataset related to the manuscript entitled "Three stages of plateau evolution manifested in present-day Tibetan Plateau"

This dataset contains the raw and cluset-stacked receiver functions and phase velocities and ZH ratios of Rayleigh wave used in the joint inversion of this study, as well as the inversion codes and the inverted seismic velocity models.

Code

CODE-MCMC_iniversion: codes used for joint inversion

README FILE

# System requirements:

This code is supported for Linux and MacOS.

The code has been tested on CentOS Linux 7 and MacOS Sequoia 15.1.

The code is under the standard Python (>3.5) environment with the following modules installed:

obspy scipy numpy matplotlib dill disba tqdm pykrige patsy rcvFnpy

# Installation guide:

These Python modules (except for the module rcvFnpy) can be installed using Python's pip tool with the command: 

pip3 install module_name

The rcvFnpy can be installed following the step:

Enter the Demo/rcvFnPy directory and run pip3 install ./.

This installation process took ~20 min on our desktop computer.

# Demo and Intructions for use

Enter the Demo/without_Pls directory and run:

python3 CrustUpperMantle_without_Pls.py ./inverted/SC.AXI.baz_60_140

Or enter the Demo/with_Pls directory and run:

python3 CrustUpperMantle_with_Pls.py ./inverted/T1.KJC02.baz_110_120

Afterward, you will get the inversion results for the two demos. This process took ~10 hours on our desktop computer.

If you want to run this code on other data, please format the data according to the demo and run the code.

Data

DATA-RF_raw: raw receiver functions (0 to 12 s after the direct-P wave) used in the joint inversion of this study

Stored in SAC (Seismic Analysis Code) format

DATA-RF_cluster_stacked: cluster-stacked receiver functions (1 to 12 s after the direct-P wave) used in the joint inversion of this study

Stored in SAC (Seismic Analysis Code) format

DATA-PV: phase velocities (8 to 70 s in period) used in the joint inversion of this study

Stored in `.txt` text format

Format: periods (s), phase velocity (km/s), phase velocity uncertainty(km/s)

DATA-ZH: ZH ratios (20 to 80 s in period) used in the joint inversion of this study

Stored in `.txt` text format

Format: periods (s), ZH ratio, ZH ratio uncertainty

RESULT-Models: inversion results of seismic velocity profiles in this study

LVZ-free: inverted models for the "LVZ-free" type profiles

Stored in `.txt` text format

Format: layer thickness (km), Vs (km/s), Vp (km/s)

LVZ-smooth: inverted models for the "LVZ-smooth" type profiles

Stored in `.txt` text format

Format: layer thickness (km), Vs (km/s), Vp (km/s)

LVZ-sharp: inverted models for the "LVZ-sharp" type profiles

Stored in `.txt` text format

Format: layer thickness (km), Vs (km/s), Vp (km/s)

mean_and_uncertainty: average and uncertainty interval of all acceptable models

Stored in `.txt` text format

Format: depth (km), Vs (km/s), uncertainty interval (km/s)

station.list: station list and locations used in the joint inversion of this study

If you face any problem or issue in the usage of this dataset, please feel free to communicate with the dataset creator, Shihua Cheng ([e-mail]).