Crustal structure study based on principal component analysis of receiver functions

作     者：Jianyong ZHANG Ling CHEN Xu WANG 

作者机构：State Key Laboratory of Lithospheric Evolution Institute of Geology and Geophysics Chinese Academy of Sciences University of Chinese Academy of Sciences CAS Center for Excellence in Tibetan Plateau Earth Sciences 

出 版 物：《Science China Earth Sciences》 (中国科学（地球科学英文版）)

年 卷 期：2019年第62卷第7期

页      面：1110-1124页

核心收录：

学科分类：07[理学] 0708[理学-地球物理学] 0704[理学-天文学] 

基　　金：supported by the National Natural Science Foundation of China (Grant No. 41688103) the independent project of the State Key Laboratory of the Lithospheric Evolution, IGGCAS (Grant No. SKL-Z201704-11712180) 

主　　题：Principal component analysis Receiver function Crustal structure Dipping discontinuity Anisotropy Sichuan Basin 

摘      要：The receiver function(RF) technique is an effective method for studying crustal structure. For a single station, the average 1-D crustal structure is usually derived by stacking the radial RFs from all back-azimuths, whereas structural variations(such as dipping discontinuities or anisotropy) can be constrained through analysis of waveform dependence on the backazimuth of both the radial and tangential RFs. However, it is often difficult to directly extract information about structural variations from the waveform of RF, due to the common presence of noise in real data. In this study, we proposed a new method to derive structural variation information for individual stations by applying principal component analysis(PCA) to RFs sorted by back-azimuth. In this method(termed as RF-PCA), a set of principal components(PCs), which are uncorrelated with each other and reflect different characteristics of the RF data, were extracted and utilized separately to reconstruct new RFs. Synthetic tests show that the first PC of the radial RFs contains the average structural information of the crust beneath the corresponding station, and the second PC of the radial RFs and the first PC of the tangential RFs both reflect the variations of the crustal structure. Our synthetic modeling results indicate that the new RF-PCA method is valid for a variety of synthetic models with intra-crustal dipping discontinuities and/or anisotropy. We applied this method to the real data from a broadband temporary seismic station(s233) in the central part of the Sichuan Basin. The results suggest that the RF data can be best explained by the presence of two nearly parallel dipping discontinuities within the crust. Combining with previous logging data, seismic exploration and deep sounding observations, we interpret the shallow dipping discontinuity as the top boundary of the Precambrian crystalline basement of the Sichuan Basin and the deep one corresponding to the Conrad interface between the upper and lower crust, consistent with the geological feature of the study area. In this work, both synthetic tests and real data application results demonstrate the effectiveness of the RF-PCA method for studying crustal structures.