A two-step multi-frequency receiver function inversion method for shallow crustal S-wave velocity structure and its application across the basin-mountain range belts in Northeast China

A shallow crustal velocity structure(above 10 km depth) is essential for understanding the crustal structures and deformation and assessing the exploration prospect of natural resources, and also provides priori information for imaging deeper crustal and mantle structure. Passive-source seismic methods are cost-effective and advantageous for regional-scale imaging of shallow crustal structures compared to active-source methods. Among these passive methods, techniques utilizing receiver function waveforms and/or body-wave amplitude ratios have recently gained prominence due to their relatively high spatial resolution. However, in basin regions, reverberations caused by near-surface unconsolidated sedimentary layers often introduce strong non-uniqueness and uncertainty, limiting the applicability of such methods. To address these challenges, we propose a two-step inversion method that uses multi-frequency P-RF waveforms and P-RF horizontal-to-vertical amplitude ratios. Synthetic tests indicate that our two-step inversion method can mitigate the non-uniqueness of the inversion and enhance the stability of the results. Applying this method to teleseismic data from a linear seismic array across the sedimentary basins in Northeast China, we obtain a high-resolution image of the shallow crustal S-wave velocity structure along the array. Our results reveal significant differences between the basins and mountains. The identification of low-velocity anomalies(<2.8 km s^(-1)) at depths less than 1.0 km beneath the Erlian Basin and less than 2.5 km beneath the Songliao Basin suggests the existence of sedimentary layers. Moreover, the high-velocity anomalies(~3.4–3.8 km s^(-1)) occurring at depths greater than 7 km in the Songliao Basin may reflect mafic intrusions emplaced during the Early Cretaceous. Velocity anomaly distribution in our imaging result is consistent with the location of the major faults, uplifts, and sedimentary depressions, as well as active-source seismic results. This application further validates the effectiveness of our method in constraining the depth-dependent characteristics of the S-wave velocity in basins with unconsolidated sedimentary cover.