We present a new digital model (NCcrust) of the seismic crustal structure of the Neoarchean North China Craton (NCC) and its surrounding Paleozoic-Mesozoic orogenic belts (30 degrees-45 degrees N, 100 degrees -130 degrees E). All available seismic profiles, complemented by receiver function interpretations of crustal thickness, are used to constrain a new comprehensive crustal model NCcrust. The model, presented on a 0.25 degrees x 0.25 degrees grid, includes the Moho depth and the internal structure (thickness and velocity) of the crust specified for four layers (the sedimentary cover, upper, middle, and lower crust) and the Pn velocity in the uppermost mantle. The crust is thin (30-32 km) in the east, while the Moho depth in the western part of the NCC is 38-44 km. The Moho depth of the Sulu-Dabie-Qinling-Qilian orogenic belt ranges from 31 km to 51 km, with a general westward increase in crustal thickness. The sedimentary cover is 2-5 km thick in most of the region, and typical thicknesses of the upper crust, middle crust, and lower crust are 16-24 km, 6-24 km, and 0-6 km, respectively. We document a general trend of westward increase in the thickness of all crustal layers of the crystalline basement and as a consequence, the depth of the Moho. There is no systematic regional pattern in the average crustal Vp velocity and the Pn velocity. We examine correlation between the Moho depth and topography for seven tectonic provinces in the North China Craton and speculate on mechanisms of isostatic compensation.