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:
-
vs=vsv*(1+gamma)/2.0
-
xi=2*(gamma-1)/(gamma+1)*100%
##############
2021/08/08
- check parallel computation can be used
- add roughness computation
- output Rayleigh/Love raypath at the final iteration
- add some useful scripts in utils/
- check noiselevel can be used in the synthetic tests
##############