We present a digital crustal model in North China Craton (NCC). The construction of crustal model is based on digitization of original seismic sounding profiles, and new results of three-dimensional structure images of receiver functions. The crustal model includes seismic velocity and thickness of crustal layers. The depths to Moho indicate a thinning crust ~30 km in the east areas and a general westward deepening to more than 40 km in the west. The P wave velocity varies from 2.0 to 5.6 km/s in the sedimentary cover, from 5.8 to 6.4 km/s in the upper crust, and from 6.5 to 7.0 km/s in the lower crust. By analyzing regional trends in crustal structure and links to tectonic evolution illustrated by typical profiles, we conclude that: (1) The delimited area by the shallowing Moho in the eastern NCC represents the spatial range of the craton destruction. The present structure of the eastern NCC crust retains the tectonic information about craton destruction by extension and magmatism; (2) The tectonic activities of the craton destruction have modified the crustal structure of the convergence boundaries at the northern and southern margin of the NCC; (3) The Ordos terrene may represent a relatively stable tectonic feature in the NCC, but with the tectonic remnant of the continental collision during the assembly of the NCC in the north-east area and the response to the lateral expansion of the Tibetan Plateau during the Cenozoic in the south-west.
A method for reconstructing crustal velocity structure using the optimization of stacking receiver function amplitude in the depth domain, named common conversion amplitude (CCA) inversion, is presented. The conversion amplitude in the depth domain, which represents the impedance change in the medium, is obtained by assigning the receiver function amplitude to the corresponding conversion position where the P-to-S conversion occurred. Utilizing the conversion amplitude variation with depth as an optimization objective, imposing reliable prior constraints on the structural model frame and velocity range, and adopting a stepwise search inversion technique, this method efficiently weakens the tendency of easily falling into the local extremum in conventional receiver function inversion. Synthetic tests show that the CCA inversion can reconstruct complex crustal velocity structures well and is especially suitable for revealing crustal evolution by estimating diverse velocity distributions. Its performance in reconstructing crustal structure is superior to that of the conventional receiver function imaging method.