Abstract:Receiver function techniques are widely used in imaging crustal and mantle structure beneath a seismic station. The weak P-to-S conversions at deep seismic structures are usually masked by strong shallow reverberations when unconsolidated sediments are present below the station, making it nearly impossible to utilize receiver function techniques. We develop a method to estimate sediment and crustal structures beneath a seismic station based on wavefield downward continuation and decomposition method. The method parametrizes velocity structure beneath the station with a stack of constant velocity layers overlying a homogeneous half-space, and approximates the teleseismic P wave and its coda by the structural response to an incoming plane P wave. Our method is based on the principle that the upgoing S wavefield is absent in the half-space, and searches for the optimum velocity and thickness of the layers that give the minimum S-wave energy flux from the half-space to the layers. An iterative grid-search algorithm from the top to the bottom layers is employed to implement the search. In this study, we only use models comprising either only one crust layer or two layers (sediment + crust) with a half-space mantle, although models with more layers are also implementable. The method is especially useful in resolving seismic structure beneath a station sitting on unconsolidated sediments. It not only can be used to determine the sediment thickness and velocity structure, but also provides an effective way to generate subsurface receiver functions, which are formed by deconvolving the upgoing P wavefield from the upgoing S waves at the top of hardrock crust, and thus are free from shallow reverberations. The technique is applied to various synthetic data generated with different types of velocity model and noise levels, and appears to have good capability in recovering the input models. We further applied this method to teleseismic data recorded at a station inside the Songliao Basin in northeast China. The estimated sediment thickness and velocity agrees well with the results of previous active-source studies. The subsurface receiver functions also show a superior power in exposing the Moho Ps conversions, resulting in a well-defined peak in the H-kappa domain, which are absent in the regular receiver function data.

P-Sv Wave-Field Connection Technique for Regional Wave Propagation Simulation

THE EFFECT OF 2-DIMENSIONAL TERRAIN WITH POINT-SOURCE ON RESISTIVITY SURVEYS

Oscillatory Responses of a Spatially Random Seabed in the Vicinity of Structures: Validation and Application of the Integrated CFD-SFEM

Impact of Multiphase Flow Simulation of Breaking Waves on a Vertical Seawall with a Recurved Parapet

Investigation and Use of Surface-wave Characteristics for Near-surface Applications

Modeling Of Rayleigh-Wave Propagation In Three-Dimensional Media With Topography Using Staggered-Grid Finite Difference Scheme

Numerical Investigation of Rayleigh-wave Propagation on Topography Surface

Study of Transmission Characteristics of Rayleigh Waves on an Elastic Wedge

Seismic wave-field modeling of 3D irregularities in layered half-space under P-, S-, and Rayleigh waves with arbitrary incident azimuth

Effects of topographic amplification induced by a U-shaped canyon on seismic waves

Pseudospectral modeling and dispersion analysis of Rayleigh waves in viscoelastic media

Multidomain Finite-Difference and Chebyshev Pseudospectral Hybrid Method for Elastic Wave Simulation: Two-Dimensional Case

Viscoelastic wave finite-difference modeling in the presence of topography with adaptive free-surface boundary condition

Simulate the elastic wavefields in media with an irregular surface topography based on staggered grid finite difference

A Flexible Wavefield Simulation Method for Layered Viscoelastic Media with Dipping Interfaces

Estimating Sedimentary and Crustal Structure Using Wavefield Continuation: Theory, Techniques and Applications

The seismic response analysis of three-dimensional Canyon complex topography under incident SV seismic waves

Enhancing Crustal Velocity Structure in Sedimentary Basin by joint inversion of Teleseismic P-Wave Reverberations and Surface Wave Dispersion

Time-domain Method for Free Field in Layered Half Space under P-SV Waves of Oblique Incidence

A New Method to Constrain Shallow Crustal S-wave Velocities Based on Direct P-wave Amplitudes in Receiver Functions and Its Application in Northeastern Tibet

Viscous-spring artificial boundary and its application to study on seismic effect of basin