Calculate PRFs of the station CB.NJ2¶
Here, we provide an example for calculating PRFs of CB.NJ2, which is a permanent station at Nanjing China mainland.
Download this example¶
Download link: ex-prf.tar.gz
unzip the file to any directory.
wget https://osf.io/dxcfz/download -O ex-prf.tar.gz tar -xzf ex-prf.tar.gz
the package include the seismic records of the station CB.NJ2 and configure files for RF calculation and H-k stacking.
The seismic records of CB.NJ2 were cut out from origin time of seismic events to 1000 s after the origin time, which contained the direct P wave and its conversion phases.
Estimate the Moho depth and the crustal Vp/Vs using the H-k stacking method¶
The PRFs had been saved to
./RFresult/CB.NJ2, and the
CB.NJ2finallist.dat had been created via the virtual check. The Moho depth can be estimated using the the H-k stacking method.
The Seispy has provided a command
hk to do this process.
usage: hk [-h] [-v] cfg_file HK stacking for single station positional arguments: cfg_file Path to HK configure file optional arguments: -h, --help show this help message and exit -v Display results to standard output
Prepare a configure file for
Here we provide an
hk.cfg in this example. See Templates of configurations in detail.
Run this example¶
The figure would be generated to
hkpath, which shows the stacking energy of difference phases and the final result of the estimation.
The result of the Moho depth and Vp/Vs were saved into the path of
hklst (defaults to