2025-01-14 12:02:34 +08:00
|
|
|
Direct surface wave radial anisotropy tomography package (DSurfRTomo)
|
|
|
|
|
|
|
|
|
|
please refer to:
|
|
|
|
|
|
|
|
|
|
Hu, S., H. Yao, and H. Huang (2020), Direct surface wave radial
|
|
|
|
|
anisotropy tomography in the crust of the eastern Himalayan
|
|
|
|
|
syntaxis. Journal of Geophysical Research: Solid Earth, 125,
|
|
|
|
|
e2019JB018257. https://doi.org/10.1029/2019JB018257
|
|
|
|
|
|
|
|
|
|
for details of the algorithm
|
|
|
|
|
|
|
|
|
|
The code is based on previous studies, especially on the implementation
|
|
|
|
|
of DSurfTomo:
|
|
|
|
|
|
|
|
|
|
Fang, H., H. Yao, H. Zhang, Y-C Huang, and R. D. van der Hilst (2015)
|
|
|
|
|
Direct inversion of surface wave dispersion for three-dimensional
|
|
|
|
|
shallow crustal structure based on ray tracing methodology and
|
|
|
|
|
application. Geophysical Journal International, 201, 1251-1263.
|
|
|
|
|
|
|
|
|
|
Please also refer to:
|
|
|
|
|
|
|
|
|
|
Rawlinson, N. and M. Sambridge (2004) Wave front evolution in
|
|
|
|
|
strongly heterogeneous layered media using the fast marching method,
|
|
|
|
|
Geophysical Journal International, 156(3), 631-647
|
|
|
|
|
|
|
|
|
|
for implementation of the fast marching method, and
|
|
|
|
|
|
|
|
|
|
Herrmann, R. B. (2013) Computer programs in seismology: An evolving
|
|
|
|
|
tool for instruction and research. Seismological Research Letter,
|
|
|
|
|
84(6), 1081-1088
|
|
|
|
|
|
|
|
|
|
for implementation of the 1-D surface wave dispersion kernel.
|
|
|
|
|
|
|
|
|
|
The dispersion data (ALLR.dat for Rayleigh wave
|
|
|
|
|
and ALLT.dat for Love wave), resulting model (DSurfRTomo.inMeasurement.dat)
|
|
|
|
|
in the crust of the eastern Himalayan syntaxis is provided in example/
|
|
|
|
|
|
|
|
|
|
#############
|
|
|
|
|
|
|
|
|
|
2019/05/18
|
|
|
|
|
The code may still need some modification
|
|
|
|
|
|
|
|
|
|
##############
|
|
|
|
|
|
|
|
|
|
output (default DRadiSurfTomo.inMeasure.dat) is in the format
|
|
|
|
|
: lon lat dep vsv gamma
|
|
|
|
|
|
|
|
|
|
To compute average shear wave velocity (vs) and radial anisotropy (xi),
|
|
|
|
|
use the following equaitons:
|
|
|
|
|
|
|
|
|
|
1. vs=vsv*(1+gamma)/2.0
|
|
|
|
|
|
|
|
|
|
2. xi=2*(gamma-1)/(gamma+1)*100%
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
##############
|
|
|
|
|
|
|
|
|
|
2021/08/08
|
|
|
|
|
1. check parallel computation can be used
|
|
|
|
|
2. add roughness computation
|
|
|
|
|
3. output Rayleigh/Love raypath at the final iteration
|
|
|
|
|
4. add some useful scripts in utils/
|
|
|
|
|
5. check noiselevel can be used in the synthetic tests
|
|
|
|
|
|
|
|
|
|
##############
|
2025-01-14 12:01:47 +08:00
|
|
|
|
