Multistep modelling of teleseismic receiver functions combined with constraints from seismic tomography: crustal structure beneath southeast China

With a growing number of modern broad-band seismographic stations in Asia, the conditions have improved to allow higher resolution structural studies on regional scales. Here, we perform a receiver-based study of the lithosphere of southeast China using waveform records of excellent quality from 14 Chinese National Digital Seismic Network and four Global Seismic Network stations. Calculating the theoretical receiver functions (RFs) that match the observed RFs from teleseismic waveforms is an established technique for retrieving information about crustal and upper mantle structure beneath a seismic receiver. RFs, however, are predominantly sensitive to the gradients in the lithospheric elastic parameters, and it is impossible to determine a non-unique distribution of seismic parameters such as absolute shear wave speeds as a function of depth unless other geophysical data are combined with RFs. Thus, we combine RFs with independent information from shear and compressional wave speeds above and below the Mohorovicic discontinuity, available from the existing tomographic studies. We introduce a statistical approach for automatically selecting only mutually coherent RFs from a large set of observed waveforms. Furthermore, an interactive forward modelling software is introduced and applied to observed RFs to define a prior, physically acceptable range of elastic parameters in the lithosphere. This is followed by a grid-search for a simple crustal structure. An initial model for a linearized, iterative inversion is constructed from multiple constraints, including results from the grid-search for shear wave speed, the Moho depth versus v(p)/v(s) ratio domain search and tomography. The thickness of the crust constrained by our multistep approach appears to be more variable in comparison with tomographic studies, with the crust thinning significantly towards the east. We observe low values of v(p)/v(s) ratios across the entire region, which indicates the presence of a very silicic crust. We do not observe any correlation between the crustal thickness or age of the crust with v(p)/v(s) ratios, which argues against a notion that there is a simple relationship between mineralogical composition and crustal thickness and age on a global scale.