Journal of Geosciences and Geomatics
ISSN (Print): 2373-6690 ISSN (Online): 2373-6704 Website: Editor-in-chief: Maria TSAKIRI
Open Access
Journal Browser
Journal of Geosciences and Geomatics. 2015, 3(4), 109-115
DOI: 10.12691/jgg-3-4-3
Open AccessArticle

Delineation of Crustal Structure at VLC Seismographic Station Using Joint Inversion of Receiver Functions and Dispersion Data

Ohaegbuchu H. E.1, , Igboekwe M. U.1 and Chukwu G. U.1

1Department of Physics, College of Physical and Applied Sciences, Michael Okpara University of Agriculture, Umudike, Abia State, Nigeria

Pub. Date: August 25, 2015

Cite this paper:
Ohaegbuchu H. E., Igboekwe M. U. and Chukwu G. U.. Delineation of Crustal Structure at VLC Seismographic Station Using Joint Inversion of Receiver Functions and Dispersion Data. Journal of Geosciences and Geomatics. 2015; 3(4):109-115. doi: 10.12691/jgg-3-4-3


Joint inversion of body wave receiver functions and dispersion data was used to model the shear wave velocity distribution of the crust and upper mantle below VLC (44.16° N, 10.39° E), a broadband seismographic station in Italy. Receiver functions are primarily sensitive to the shear wave velocity contrasts and vertical travel times, whereas the surface wave dispersion measurements are sensitive to absolute vertical shear wave velocity averages and changes as function of depth. Each data set has inherent lapses but by jointly inverting both we are able to draw on the capabilities of one to compensate the imperfections of the other and this provides better S-wave velocity constraints than we would obtain by inverting either data set individually. The receiver functions were computed from the teleseismic earthquakes recorded by VLC station between 2005 and 2012, while the dispersion curves at regional scale were determined by the Frequency Time Analysis and have been used to obtain tomography maps, using the two-dimensional tomography algorithm developed by Ditmar and Yanovskaya in 1987. The inversion results include a crust with a sharp gradient near the surface (shear velocity changing from 2.15 to 3.4 kms-1 in 5 km) underlain by a 13-km-thick layer with a shear velocity of and 3.4 kms-1 another 15-km- thick layer with a shear velocity of 3.72 kms-1, and an upper mantle with an average shear velocity of 4.4 kms-1. The crust–mantle transition has a significant gradient, with velocity values varying from 3.72 to 4.4 kms-1 at about 32 km depth. This result is also in agreement with shear wave velocity cross-section of the area obtained from ITA-LSO data sets.

delineation joint inversion receiver functions dispersion data dispersion curves crustal structure

Creative CommonsThis work is licensed under a Creative Commons Attribution 4.0 International License. To view a copy of this license, visit


Figure of 8


[1]  A. L. Levshin, T. B. Yanovskaya, A. V. Lander, B. G. Bukchin, M. P. Barmin, and L. I. Ratnicova (1989) Surface Seismic Waves in Laterally Inhomogenous Earth. In: Keilis-Borok VI (ed).
[2]  C. J. Ammon, G. E. Randall, and G. Zandt (1990) On the nonuniqueness of receiver function inversions. J. Geophys. Res., 95:15303-15318.
[3]  C. A. Langston (1979) Structure under Mount Rainier, Washington, inferred from the teleseismic body wages. J. Geophys.Res., 84:4749-4762.
[4]  Enrico Brandmayr, Reneta Blagoeva Raykova, Marco Zuri, Fabio Romanelli, Carlo Doglioni, and Giuliano Francesco Panza (2010) The Lithosphere in Italy: Structure and seismicity. Journal of the Virtual Explorer., 36:20-44.
[5]  G. Boyandzhiev, E. Brandmayr, T. Pinat, and G. F. Panza (2008) Optimization for nonlinear inversion. Earth and Environmental Science, 19:17-43.
[6]  G. F. Panza (1981) The resolving power of seismic surface waves with respect to crust and upper mantle structural models. In: Cassinis R (ed) The solution of the inverse problem in geophysical interpretation. Plenum Press New York.
[7]  G. F. Panza, A. Peccerillo, A. Aoudia, and B. M. Farina (2007) Geophysical and petrological modelling of the structure and composition of the crust and upper mantle in complex geodynamic settings: The Tyrrhenian Sea and surroundings. Earth Science Reviews, 80:1-46.
[8]  J. Juli`a, C. J. Ammon, R. B. Hermann, and A. M. Correig. (2000) Joint inversion of receiver function and surface wave dispersion observation. Geophys. J. Int., 143:99-112.
[9]  O. Gonz ́alez, L. Alvarez, M. Guidarelli, and G. F. Panza (2007) Crust and upper mantle structure in the Caribbean region by group velocity tomography and regionalization. Pure Appl. Geophys., 164:1985-2007.
[10]  O’Leary Gonz ́alez, Bladimir Moreno, Fabio Romanelli, and Giuliano F. Panza (2012) Lithospheric structure below seismic stations in Cuba from the joint inversion of Rayleigh surface waves dispersion and receiver functions. Geophys. J. Int., 189:1047-1059.
[11]  R. B. Hermann and C. J. Ammon (2002) Computer programs in seismology: surface waves, receiver functions and crustal structure.
[12]  T. B. Yanovskaya and G. M. Molchan (2001) Development of methods for surface wave tomography based on Backus–Gilbert approach. Computational Seismology. In: Keilis-Boro, VI (ed).