Advanced Search



ISSN  2096-3955

CN  10-1502/P

Citation: XinYan Zhang, ZhiMing Bai, Tao Xu, Rui Gao, QiuSheng Li, Jue Hou, José Badal, 2018: Joint tomographic inversion of first-arrival and reflection traveltimes for recovering 2-D seismic velocity structure with an irregular free surface, Earth and Planetary Physics, 2, 220-230. doi: 10.26464/epp2018021

2018, 2(3): 220-230. doi: 10.26464/epp2018021


Joint tomographic inversion of first-arrival and reflection traveltimes for recovering 2-D seismic velocity structure with an irregular free surface


Key Laboratory of Deep-Earth Dynamics of Ministry of Land and Resources, Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China


State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China


CAS Center for Excellence in Tibetan Plateau Earth Sciences, Beijing 100101, China


School of Earth Science and Engineering, Sun Yat-sen University, Guangzhou 510275, China


Institute of Geophysics, China Earthquake Administration, Beijing 100081, China


Physics of the Earth, Sciences B, University of Zaragoza, Spain

Corresponding author: XinYan Zhang, zhangxinyana@163.comZhiMing Bai,

Received Date: 2018-03-02
Web Publishing Date: 2018-05-01

Irregular surface flattening, which is based on a boundary conforming grid and the transformation between curvilinear and Cartesian coordinate systems, is a mathematical method that can elegantly handle irregular surfaces, but has been limited to obtaining first arrivals only. By combining a multistage scheme with the fast-sweeping method (FSM, the method to obtain first-arrival traveltime in curvilinear coordinates), the reflected waves from a crustal interface can be traced in a topographic model, in which the reflected wave-front is obtained by reinitializing traveltimes in the interface for upwind branches. A local triangulation is applied to make a connection between velocity and interface nodes. Then a joint inversion of first-arrival and reflection traveltimes for imaging seismic velocity structures in complex terrains is presented. Numerical examples all perform well with different seismic velocity models. The increasing topographic complexity and even use of a high curvature reflector in these models demonstrate the reliability, accuracy and robustness of the new working scheme; checkerboard testing illustrates the method’s high resolution. Noise tolerance testing indicates the method’s ability to yield practical traveltime tomography. Further development of the multistage scheme will allow other later arrivals to be traced and used in the traveltime inversion.

Key words: irregular surface flattening, boundary conforming grid, multistage scheme, traveltime tomography

Afnimar, and Koketsu, K. (2000). Finite difference traveltime calculation for head waves travelling along an irregular interface. Geophys. J. Int., 143(3), 729–734.

Appelö, D., and Petersson, N. A. (2009). A stable finite difference method for the elastic wave equation on complex geometries with free surfaces. Commun. Comput. Phys., 5(1), 84–107.

Artemieva, I. M. (2003). Lithospheric structure, composition, and thermal regime of the East European Craton: implications for the subsidence of the Russian platform. Earth Planet. Sci. Lett., 213(3-4), 431–446.

Bai, C. Y., Huang, G. J., and Zhao, R. (2009). 2-D/3-D irregular shortest-path ray tracing for multiple arrivals and its applications. Geophys. J. Int., 183(3), 1596–1612.

Bai, C. Y., Tang, X. P., and Zhao, R. (2009). 2-D/3-D multiply transmitted, converted and reflected arrivals in complex layered media with the modified shortest path method. Geophys. J. Int., 179(1), 201–214.

Bai, Z. M., Zhang, Z. J., and Wang, Y. H. (2007). Crustal Structure across the Dabie-Sulu orogenic belt revealed by seismic velocity profiles. J. Geophys. Eng., 4(4), 436–442.

Benamou, J. D. (1996). Big ray tracing: Multivalued travel time field computation using viscosity solutions of the Eikonal equation. J. Comput. Phys., 128(2), 463–474.

Cao, S. H., and Greenhalgh, S. (1994). Finite-difference solution of the eikonal equation using an efficient, first-arrival, wavefront tracking scheme. Geophysics, 59(4), 632–643.

Carbonell, R., Simancas, F., Juhlin, C., Pous, J., Pérez-Estaún, A., Gonzalez-Lodeiro, F., Muñoz, G., Heise, W., and Ayarza, P. (2004). Geophysical evidence of a mantle derived intrusion in SW Iberia. Geophys. Res. Lett., 31(11), L11601.

Cassell, B. R. (1982). A method for calculating synthetic seismograms in laterally varying media. Geophys. J. R. Astron. Soc., 69(2), 339–354.

Fornberg, B. (1988). The pseudospectral method: Accurate representation of interfaces in elastic wave calculations. Geophysics, 53(5), 625–637.

Gao, R., Lu, Z. W., Li, Q. S., Guan, Y., Zhang, J. S., He, R. Z., and Hang, L. Y. (2005). Geophysical survey and geodynamic study of crust and upper mantle in the Qinghai-Tibet Plateau. Episodes, 28(4), 263–273.

Haines, A. J. (1988). Multi-source, multi-receiver synthetic seismograms for laterally heterogeneous media using F-K domain propagators. Geophys. J. Int., 95(2), 237–260.

Hole, J. (1992). Nonlinear high‐resolution three‐dimensional seismic travel time tomography. J. Geophys. Res., 97(B5), 6553–6562.

Huang G. J., Bai C. Y., Zhu D. I., and Greenhalgh S. (2012). 2D/3D seismic simultaneous inversion for the velocity and interface geometry using multiple classes of arrivals. Bull. Seismol. Soc. Am., 102(2), 790–801.

Hvid, S. (1994). Three Dimensional Algebraic Grid Generation. Copenhagen, Denmark: Technical University of Denmark.222

Julian, B. R., and Gubbins, D. (1977). Three-dimensional seismic ray tracing. J. Geophys. Res., 43, 95–113.

Kaila, K. L., and Krishna, V. G. (1992). Deep seismic sounding studies in India and major discoveries. Curr. Sci., 62(1-2), 117–154.

Kao, C. Y., Osher, S., and Qian, J. L. (2008). Legendre-transform-based fast sweeping methods for static Hamilton–Jacobi equations on triangulated meshes. J. Comput. Phys., 227(24), 10209–10225.

Kim, S., and Cook, R. (1999). 3-D traveltime computation using second-order ENO scheme. Geophysics, 64(6), 1867–1876.

Kimmel, R., and Sethian, J. A. (1998). Computing geodesic paths on manifolds. Proc. Natl. Acad. Sci. USA, 95(15), 8431–8435.

Knapp, C. C., Knapp, J. H., and Connor, J. A. (2004). Crustal-scale structure of the South Caspian Basin revealed by deep seismic reflection profiling. Mar. Pet. Geol., 21(8), 1073–1081.

Koketsu, K., and Sekine, S. (1998). Pseudo-bending method for three-dimensional seismic ray tracing in a spherical earth with discontinuities. Geophys. J. Int., 132(2), 339–346.

Lan, H. Q., and Zhang, Z. J. (2011a). Comparative study of the free-surface boundary condition in two-dimensional finite-difference elastic wave field simulation. J. Geophys. Eng., 8(2), 275–286.

Lan, H. Q., and Zhang, Z. J. (2011b). Three-dimensional wave-field simulation in heterogeneous transversely isotropic medium with irregular free surface. Bull. Seismol. Soc. Am., 101(3), 1354–1370.

Lan, H. Q., Zhang, Z., Xu, T., and Bai, Z. M. (2012). Influences of anisotropic stretching of boundary conforming grid on traveltime computation by topography-dependent eikonal equation. Chin. J. Geophys., 55(5), 564–579.

Lan, H. Q., and Zhang, Z. J. (2013a). Topography-dependent eikonal equation and its solver for calculating first-arrival traveltimes with an irregular surface. Geophys. J. Int., 193(2), 1010–1026.

Lan, H. Q, and Zhang, Z. J. (2013b). A high-order fast-sweeping scheme for calculating first-arrival travel times with an irregular surface. Bull. Seismol. Soc. Am., 103(3), 2070–2082.

Lelièvre, P. G., Farquharson, C. G., and Hurich, C. A. (2011). Computing first-arrival seismic traveltimes on unstructured 3-D tetrahedral grids using the fast marching method. Geophys. J. Int., 184(2), 885–896.

Li, S. L., and Mooney, W. D. (1998). Crustal structure of China from deep seismic sounding profiles. Tectonophysics, 288(1-4), 105–113.

Ma, T., and Zhang, Z. J. (2014a). Calculating ray paths for first-arrival travel times using a topography-dependent eikonal equation solver. Bull. Seismol. Soc. Am., 104(3), 1501–1517.

Ma, T., and Zhang, Z. J. (2014b). A model expansion criterion for treating surface topography in ray path calculations using the eikonal equation. J. Geophys. Eng., 11(2), 025007.

Ma, T., and Zhang, Z. J. (2015). Topography-dependent eikonal traveltime tomography for upper crustal structure beneath an irregular surface. Pure Appl. Geophys., 172(6), 1511–1529.

Podvin, P., and Lecomte, I. (1991). Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools. Geophys. J. Int., 105(1), 271–284.

Qian, J. L., and Symes, W. W. (2002). An adaptive finite-difference method for traveltimes and amplitudes. Geophysics, 67(1), 167–176.

Qian, J. L., Zhang, Y. T., and Zhao, H. K. (2007a). A fast sweeping method for static convex Hamilton-Jacobi equations. J. Sci. Comput., 31(1-2), 237–271.

Qian, J. L., Zhang, Y. T., and Zhao, H. K. (2007b). Fast sweeping methods for Eikonal equations on triangular meshes. SIAM J. Numer. Anal., 45(1), 83–107.

Qin, F. H., Luo, Y., Olsen, K. B., Cai, W. Y., and Schuster, G. T. (1992). Finite-difference solution of the eikonal equation along expanding wavefronts. Geophysics, 57(3), 478–487.

Rawlinson, N., and Sambridge, M. (2004a). Multiple reflection and transmission phases in complex layered media using a multistage fast marching method. Geophysics, 69(5), 1338–1350.

Rawlinson, N., and Sambridge, M. (2004b). Wave front evolution in strongly heterogeneous layered media using the fast marching method. Geophys. J. Int., 156(3), 631–647.

Rawlinson, N., and Goleby, B. R. (2012). Seismic imaging of continents and their margins: New research at the confluence of active and passive seismology. Tectonophysics, 572-573, 1–6.

Reshef, M. (1991). Depth migration from irregular surfaces with depth extrapolation methods. Geophysics, 56(1), 119–122.

Riahi, M. A., and Juhlin, C. (1994). 3-D interpretation of reflected arrival times by finite-difference techniques. Geophysics, 59(5), 844–849.

Sambridge, M. S. (1990). Non-linear arrival time inversion: Constraining velocity anomalies by seeking smooth models in 3-D. Geophys. J. Int., 102(3), 653–677.

Scarascia, S., and Cassinis, R. (1997). Crustal structures in the central-eastern Alpine sector: a revision of the available DSS data. Tectonophysics, 271(1-2), 157–188.

Sethian, J. A. (1999). Fast marching methods. SIAM Rev., 41(2), 199–235.

Sethian, J. A., and Vladimirsky, A. (2000). Fast methods for the Eikonal and related Hamilton–Jacobi equations on unstructured meshes. Proc. Natl. Acad. Sci. USA, 97(11), 5699–5703.

Symes, W. W., and Qian, J. L. (2003). A slowness matching eulerian method for multivalued solutions of eikonal equations. J. Sci. Comput., 19(1-3), 501–526.

Teng, J. W., Wei, S. Y., Sun, K. Z., and Xue, C. S. (1987). The characteristics of the seismic activity in the Qinghai-Xizang (Tibet) Plateau of China. Tectonophysics, 134(1-3), 129–144.

Teng, J. W., Zeng, R. S., Yan, Y. F., and Zhang, H. (2003). Depth distribution of Moho and tectonic framework in eastern Asian continent and its adjacent ocean areas. Sci. China Ser. D Earth Sci., 46(5), 428–446.

Thompson, J. F., Warsi, Z. U. A., and Mastin, C. W. (1985). Numerical Grid Generation: Foundations and Applications. Amsterdam: North-Holland.222

Tian, X. B., Teng, J. W., Zhang, H. S., Zhang, Z. J., Zhang, Y. Q., Yang, H., and Zhang K. K. (2011). Structure of crust and upper mantle beneath the Ordos Block and the Yinshan Mountains revealed by receiver function analysis. Phys. Earth Planet. Inter., 184(3-4), 186–193.

Um, J., and Thurber, C. (1987). A fast algorithm for two-point seismic ray tracing. Bull. Seismol. Soc. Am., 77(3), 972–986.

van Trier, J., and Symes, W. W. (1991). Upwind finite-difference calculation of traveltimes. Geophysics, 56(6), 812–821.

Vidale, J. (1988). Finite-difference calculation of travel times. Bull. Seismol. Soc. Am., 78(6), 2062–2076.

Wang, C. Y., Zeng, R. S., Mooney, W. D., and Hacker, B. R. (2000). A crustal model of the ultrahigh-pressure Dabie Shan orogenic belt, China, derived from deep seismic refraction profiling. J. Geophys. Res., 105(B5), 10857–10869.

Wang, C. Y., Han, W. B., Wu, J. P., Lou, H., and Chan, W. W. (2007). Crustal structure beneath the eastern margin of the Tibetan Plateau and its tectonic implications. J. Geophys. Res., 112(B7), B07307.

Wu, C. L., Harris, J. M., Nihei, K. T., and Nakagawa, S. (2005). Two-dimensional finite-difference seismic modeling of an open fluid-filled fracture: Comparison of thin-layer and linear-slip models. Geophysics, 70(4), T57-T62.

Xu, T., Xu, G. M., Gao, E. G., Li, Y. C., Jiang, X. Y., and Luo, K. Y. (2006). Block modeling and segmentally iterative ray tracing in complex 3D media. Geophysics, 71(3), T41-T51.

Xu, T., Zhang, Z. J., Gao, E. G., Xu, G. M., and Sun, L. (2010). Segmentally iterative ray tracing in complex 2D and 3D heterogeneous block models. Bull. Seismol. Soc. Am., 100(2), 841–850.

Xu, T., Li, F., Wu, Z. B., Wu, C. L., Gao, E. G., Zhou, B., Zhang, Z. J., and Xu, G. M. (2014). A successive three-point perturbation method for fast ray tracing in complex 2D and 3D geological models. Tectonophysics, 627, 72–81.

Zelt, C. A., and Smith, R. B. (1992). Seismic traveltime inversion for 2-D crustal velocity structure. Geophys. J. Int., 108(1), 16–34.

Zeng, R. S., Ding, Z. F., and Wu, Q. J. (1995). A review on the lithospheric structures in the Tibetan Plateau and constraints for dynamics. Pure Appl. Geophys., 145(3-4), 425–443.

Zhang, W., and Chen, X. F. (2006). Traction image method for irregular free surface boundaries in finite difference seismic wave simulation. Geophys. J. Int., 167(1), 337–353.

Zhang, Z. J., and Klemperer, S. L. (2005). West-east variation in crustal thickness in northern Lhasa block, central Tibet, from deep seismic sounding data. J. Geophys. Res., 110(B9), B09403.

Zhang, Z. J., and Klemperer, S. (2010). Crustal structure of the Tethyan Himalaya, southern Tibet: New constraints from old wide-angle seismic data. Geophys. J. Int., 181(3), 1247–1260.

Zhang, Z. J., Deng, Y. F., Teng, J. W., Wang, C. Y., Gao, R., Chen, Y., and Fan, W. M. (2011). An overview of the crustal structure of the Tibetan plateau after 35 years of deep seismic soundings. J. Earth Sci., 40(4), 977–989.

Zhou, B., Greenhalgh, S. A., and Sinadlnovskl, C. (1992). Iterative algorithm for the damped minimum norm, least-squares and constrained problem in seismic tomography. Explor. Geophys., 23(3), 497–505.


Qing Wang, XiaoDong Song, JianYe Ren, 2017: Ambient noise surface wave tomography of marginal seas in east Asia, Earth and Planetary Physics, 1, 13-25. doi: 10.26464/epp2017003


JinHai Zhang, ZhenXing Yao, 2017: Exact local refinement using Fourier interpolation for nonuniform-grid modeling, Earth and Planetary Physics, 1, 58-62. doi: 10.26464/epp2017008


LiangQuan Ge, JianKun Zhao, QingXian Zhang, YaoYao Luo, Yi Gu, 2018: Mapping of the lunar surface by average atomic number based on positron annihilation radiation from Chang’e-1, Earth and Planetary Physics, 2, 238-246. doi: 10.26464/epp2018023

Article Metrics
  • PDF Downloads()
  • Abstract views()
  • HTML views()
  • Cited by(0)

Figures And Tables

Joint tomographic inversion of first-arrival and reflection traveltimes for recovering 2-D seismic velocity structure with an irregular free surface

XinYan Zhang, ZhiMing Bai, Tao Xu, Rui Gao, QiuSheng Li, Jue Hou, José Badal