Multi-station joint inversion of receiver function and surface-wave phase velocity data for exploration of deep sedimentary layers

Abstract

In this study, we propose a joint inversion method, using genetic algorithms, to estimate an S-wave velocity structure for deep sedimentary layers from receiver functions and surface-wave phase velocity observed at several sites. The method takes layer continuity over a target area into consideration by assuming that each layer has uniform physical properties, especially an S-wave velocity, at all the sites in a target area in order to invert datasets acquired at different sites simultaneously. Numerical experiments with synthetic data indicate that the proposed method is effective in reducing uncertainty in deep structure parameters when modelling only surface-wave dispersion data over a limited period range. We then apply the method to receiver functions derived from earthquake records at one site and two datasets of Rayleigh-wave phase velocity obtained from microtremor array surveys performed in central Tokyo, Japan. The estimated subsurface structure is in good agreement with the results of previous seismic refraction surveys and deep borehole data. We also conclude that the proposed method can provide a more accurate and reliable model than individual inversions of either receiver function data only or surface-wave dispersion data only.