Next Article in Journal
Coupling Relationship Analysis of Gold Content Using Gaofen-5 (GF-5) Satellite Hyperspectral Remote Sensing Data: A Potential Method in Chahuazhai Gold Mining Area, Qiubei County, SW China
Previous Article in Journal
UAV Thermal Images for Water Presence Detection in a Mediterranean Headwater Catchment
Previous Article in Special Issue
Analysis of Seismic Deformation from Global Three-Decade GNSS Displacements: Implications for a Three-Dimensional Earth GNSS Velocity Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deformation of the Crust and Upper Mantle beneath the North China Craton and Its Adjacent Areas Constrained by Rayleigh Wave Phase Velocity and Azimuthal Anisotropy

1
Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
2
Key Laboratory of Earthquake Source Physics, China Earthquake Administration, Beijing 100081, China
3
School of Geosciences and Info-Physics, Central South University, Changsha 410083, China
4
School of Geophysics and Information Technology, China University of Geosciences, Beijing 100029, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(1), 110; https://doi.org/10.3390/rs14010110
Submission received: 9 December 2021 / Revised: 24 December 2021 / Accepted: 24 December 2021 / Published: 28 December 2021
(This article belongs to the Special Issue Remote Sensing and Geodynamics)

Abstract

:
The North China Craton (NCC) has experienced strong tectonic deformation and lithospheric thinning since the Cenozoic. To better constrain the geodynamic processes and mechanisms of the lithospheric deformation, we used a linear damped least squares method to invert simultaneously Rayleigh wave phase velocity and azimuthal anisotropy at periods of 10–80 s with teleseismic data recorded by 388 permanent stations in the NCC and its adjacent areas. The results reveal that the anomalies of Rayleigh wave phase velocity and azimuthal anisotropy are in good agreement with the tectonic domains in the study area. Low-phase velocities appear in the rift grabens and sedimentary basins at short periods. A rotation pattern of the fast axis direction of the Rayleigh wave together with a distinct low-velocity anomaly occurs around the Datong volcano. A NW–SE trending azimuthal anisotropy and a low-velocity anomaly at periods of 60–80 s are observed subparallel to the Zhangbo fault zone. The whole lithosphere domain of the Ordos block shows a high-phase velocity and counterclockwise rotated fast axis. The northeastern margin of the Tibetan plateau is dominated by a low-velocity and coherent NW–SE fast axis direction. We infer that the subduction of the Paleo-Pacific plate and eastward material escape of the Tibetan plateau mainly contribute to the deformation of the crust and upper mantle in the NCC.

1. Introduction

The North China Craton (NCC) is the oldest Archean craton in China where the earliest crustal rock can be dated back to 3.8 Ga [1]. In the late Paleoproterozoic (~1.85 Ga), the North China Craton was amalgamated by the collision between the eastern NCC (ENCC) and western NCC (WNCC) along the Trans-North China orogen (TNCO) until its final cratonization [2] (Figure 1a). The NCC, in the eastern Eurasian plate, is composed of several tectonic blocks including horsts and grabens shown in Figure 1b. In addition, the NCC is currently surrounded by orogenic belts formed in different geological episodes. To the north, the central Asian orogenic belt (CAOB) is a product of the closure of the Paleo-Asian ocean in the Paleozoic–Mesozoic [3]. To the southwest, the Qilian orogen (QLO) has been assembled by Paleozoic subduction–accretion [4]. To the south, the Qinling–Dabie orogenic belt (QDO) has been formed by the collision between the Yangtze plate and the NCC in the Triassic [2].
The Paleo-Pacific plate was subducted beneath the Eurasian plate in a NNW direction (Figure 1a), and then turned gradually to a westward direction from ~50–42 Ma describing a counterclockwise rotation [6]. Previous studies have shown that this subduction of the Paleo-Pacific plate below the Eurasian plate transformed the dynamic deformation of the NCC from N–S shortening to almost E–W intraplate extension, resulting in a series of large-scale extensional structures, detachment faults and rift grabens in the NCC [10,11,12,13]. The stagnant slab dehydration of the Paleo-Pacific plate in the mantle transition zone led to an unsteady asthenospheric flow, partial melting of the overlying lithosphere and magmatic activity [12,14,15,16]. Especially in the eastern NCC, the lithospheric thickness has been thinned by at least 100 km, and the original cratonic structures and properties have been significantly modified due to the lithospheric partial melting mentioned above [12,13]. The uplift and outward growth of the Tibetan plateau have been caused by the collision between the Indian and the Eurasian plates since ~50 Ma [17,18]. The Qinling–Dabie orogen between the Sichuan basin and the Ordos block was proposed as acting as an important channel for the eastward material escape of the Tibetan plateau [17,19]. This extrusion of the Tibetan plateau even affected the deformation of the crust and mantle in the Trans-North China orogen [20,21,22]. Therefore, the NCC suffers the westward subduction of the Paleo-Pacific plate in the east and eastward extrusion of the Tibetan plateau in the southwest, leading to significant deformation of the crust and upper mantle structure [10,11,23,24]. In general, cratons are generally considered to be underlain by a thick, cold and refractory mantle lithosphere [13,14,15]. Lithospheric deformation is likely to cause prominent changes in velocity and anisotropy of the crust and mantle. Thus, the velocity structure and anisotropy can effectively constrain the deformation of the crust and mantle structure in turn, which help us discuss the tectonic evolution of the NCC and its adjacent areas.
Recently, many studies revealed the velocity structure and seismic anisotropy beneath the crust and upper mantle of the NCC [23,25,26,27]. However, these results on understanding the lithospheric thinning and deformation of the NCC are still controversial. For example, the low-velocity anomalies of the asthenosphere revealed by Xu et al. [16], related to the lithospheric thinning, are segmented by irregular high-velocity anomalies in the eastern NCC, which are obviously inconsistent with the results obtained by full-waveform tomography [23], indicating that the eastern NCC is dominated by extremely low velocities. The upwelling of hot material beneath the Datong volcano is indicated by low-velocity anomalies in previous results [20,21,25,27], but its origin and mechanism are still under debate. Huang et al. [27] proposed that the asthenospheric upwelling came from the northeastern margin of the Ordos block. Yao et al. [20] confirmed this hypothesis and attributed it to lithospheric stretching. However, the results of S-wave [21,25] suggested that the hot materials originated from the asthenosphere in the vertical direction rather than the northeastern margin of the Ordos block. Chen et al. [25] further speculated that the upwelling of the asthenosphere beneath the Datong volcano would flow into the eastern NCC. Seismic anisotropy is often used to investigate the deformation of the crust and upper mantle. The S-wave splitting studies provide useful information about anisotropy in the crust and upper mantle [24,26]. However, the depth of the anisotropic structure is usually difficult to obtain by S-wave splitting. Moreover, studies on lithospheric deformation are still discussed due to the different interpretations of the NCC. Based on the fast axis directions of the local S-wave and SKS/SKKS phases, Gao et al. [28] suggested that the deformation of the crust and upper mantle was not a simple coupled pattern. The fast axis direction of the SKS/SKKS phases correlates well with the absolute plate motion (APM) direction, implying the APM-driven model of the upper mantle anisotropy in the NCC [29]. This model disagrees with the hypotheses inferred from stress-aligned microcracks, frozen fossil anisotropy and asthenospheric flow [22,25,26]. Compared with S-wave splitting, the surface wave azimuthal anisotropy estimated from different periods has a better vertical resolution and can provide a depth independent anisotropy to constrain the deformation of the crust and upper mantle.
Using teleseismic data recorded by the permanent stations in the NCC and its adjacent areas, the goal of this study is to further improve constraints on the deformation of the crust and upper mantle. We estimate simultaneously Rayleigh wave phase velocity and azimuthal anisotropy at periods of 10–80 s. By comparing the previous results, we focus on investigating the tectonic deformation and geodynamic mechanisms forming the NCC and its adjacent areas.

2. Data and Analysis

2.1. Data Acquisition and Processing

The seismological observation system in China has achieved rapid development [30] and the permanent seismic stations cover densely the NCC and its adjacent areas (Figure 1b). In this study, we collected the three-component waveform data of Rayleigh waves recorded at 388 permanent seismic stations from the China National Seismic Network from May 2013 to December 2019. These stations consist of 318 broadband seismometers, with bandwidths of 0.0167–80 Hz, and 70 very broadband seismometers, with bandwidths of 0.0083–80 Hz, including Guralp CMG-3ESPC and 3T, CTS-1, Geodevice BBVS-60 and Trillium-120 (Figure 1b).
In this study, the classical two-station method is used to measure the Rayleigh wave phase velocity. This method is not affected by errors such as the epicenter, the earthquake origin time and the initial phase, but only requires that the event and two stations are roughly on the same great circle path [25]. In order to measure the high-quality and reliable dispersion curves, we require that: (1) the magnitude of the event is not less than 5.5 and the focal depth is not greater than 50 km, which may ensure the signal-to-noise ratio of the Rayleigh wave; (2) the epicenter distance is between 20° and 120° to reduce the influence of the near-source effect and high-order surface wave interference; (3) the interstation distance should be larger than 3 wavelengths to decrease the errors in the dispersion measurement; (4) calculating the cross-correlation function with the seismograms of the same event recorded by the selected station pair after down sampling and removing the means, trends and instrument response. According to the above processing, a total of 1347 earthquake events were selected in this study (Figure 2). We found that the distribution of these earthquakes has good epicentral distance and azimuth coverage, but a considerable number of earthquakes are concentrated in areas such as the western Paleo-Pacific subduction zone, resulting in multiple great circle paths passing through a certain station pair. Therefore, we adopt the time-frequency phase weighted stacking method [31,32] to stack the waveforms of the cross-correlation function with the same great circle path for the station pair, which can improve the signal-to-noise ratio of the cross-correlation functions and measurement accuracy of Rayleigh wave phase velocity dispersion curve. Then, the automated frequency-time analysis method [33] is used to estimate the Rayleigh wave phase velocity dispersion curves in the period range of 10–80 s for all station pairs. The final dispersion curves are selected by the criteria of signal-to-noise ratios greater than 20 and phase velocities less than twice the standard deviation at each period.

2.2. Surface Wave Tomography Inversion

We use the tomographic method of Barmin et al. [34] to invert Rayleigh wave phase velocity and azimuthal anisotropy on a 0.5° × 0.5° spatial grid beneath the NCC and its adjacent areas. This method is based on the damped least squares inversion to minimize a penalty function composed of a linear combination of data misfit, model smoothness and the perturbation to a reference model. For the weakly anisotropic medium, the relationship between Rayleigh wave phase velocity c(T) with a period of T and the azimuthal angle 𝜃 can be expressed as [35]:
c(T) = 𝑐0(𝑇) + 𝑐1(𝑇)𝑐𝑜𝑠(2𝜃) + 𝑐2(𝑇)𝑠𝑖𝑛(2𝜃) + 𝑐3(𝑇)𝑐𝑜𝑠(4𝜃) + 𝑐4(T)s𝑖𝑛(4𝜃)
Here, 𝑐0 is the isotropic phase velocity of the reference model, and 𝑐1, 𝑐2 and 𝑐3, 𝑐4 are 2𝜃 and 4𝜃 anisotropic components, respectively. Generally, 4𝜃 can often be ignored in weakly anisotropic media, and then the strength of anisotropy, A, and the azimuth angle of the fast wave, 𝜃, can be obtained by the following two formulas:
A = c 1 2 + c 2 2
θ = 1 2 a r c t a n c 1 c 2  
During the inversion, the two parameters specifying the weight of smoothing and the width of the smoothing area control the damping factor, which is determined by a number of tests using different combinations of parameters by considering data misfit, model resolution and model smoothness [36]. Meanwhile, we refer to the processing steps of Guo et al. [37] to perform the surface wave tomography inversion. We first invert Rayleigh wave phase velocity and azimuth anisotropy at each period. We then remove the dispersion curves with large misfits, update the data set and then perform the inversion again to obtain the final model of Rayleigh wave phase velocity and azimuth anisotropy in this study.

2.3. Robustness Analysis

After the quality control of dispersion curves in Section 2.1, we plot Rayleigh wave phase velocity dispersion curves as a function of period in Figure 3a. The fewest dispersion curve numbers at a 10 s period are still as high as 8760. At a period of 55 s, the ray number is higher than 45,000. A considerable number of dispersion curves provide quite dense and crossed ray paths and ensure the high resolution of the final inversion. Here, we also take two station pairs as an example to check the phase velocity dispersion curves (Figure 3b). The paths of the station pairs WUL–LZH and JJ–BJT span different geological blocks in the study area (Figure 1b). The station pair WUL–LZH is from the south of the Qilian orogen, and across the Ordos block and the Trans-North China orogen. The path between stations JJ and BJT lies mostly within the Bohai Bay basin. The Rayleigh wave phase velocities at short (<15 s) and long (>40 s) periods along the path of JJ–BJT are obviously lower than that along the path of WUL–LZH because of the sediments and thin lithosphere in the Bohai Bay basin of the eastern NCC. In contrast, in the period range of 15–40 s, the phase velocities are much higher in the Bohai Bay basin, indicating the colder and stronger lower crust and uppermost mantle between stations JJ and BJT.
Histograms of data misfits to measurements from the final tomography at different periods are presented in Figure 4. The standard deviation displayed in the upper right corner of each panel increases gradually from 1.38 s to 2.43 s. This rising standard deviation is mainly due to the increase of random error and decrease of the signal-to-noise ratio in Rayleigh wave phase velocity measurements. However, the Gaussian distributions of data misfits are nearly concentrated between ±4 s, which implies that the final models obtained by this inversion fit the observed data well.
In principle, the resolution and reliability of tomographic results depend on the dense and crossed paths of Rayleigh wave phase velocity measurements. The study area is dominated by a good path coverage (Figure 5), although the rate of data rejection at each period is more than 30% based on the strict selection criteria in Section 2.1. The resolution maps are simultaneously determined during tomography (Figure 6). As described by Barmin et al. [34], the resolution is estimated from the resolution matrix and defined by a complicated function of path density. In Figure 6, we find that the better resolution gradually shrinks from the edge to the central part of the study area due to the path coverage. Overall, most of the study area is covered by a robust resolution of ~1°–2° (~100–200 km).

3. Results

3.1. Rayleigh Wave Phase Velocity

Figure 7 shows the anisotropic Rayleigh wave phase velocity maps at different periods. Anomalies of the phase velocity at periods of 15, 20, 30, 40, 60 and 80 s are demonstrated by the perturbations relative to the average values at each individual period of 3.29, 3.44, 3.70, 3.84, 3.94 and 4.01 km/s, respectively. The phase-velocity perturbations indicate the lateral heterogeneities beneath the NCC and its adjacent areas.
At the short period of 15 s (Figure 7a), Rayleigh wave phase velocities show low-velocity anomalies in the sedimentary basins, such as the Yinchuan, Hetao, Weihe and Fenhe grabens, and the Bohai Bay basin. In contrast, the Lüliang and Taihang mountains, and the Sulu orogen with old crystalline rocks present high velocities except for the young orogenic belts of the Qilian orogen and Songpan–Ganzi terrane dominated by low-velocity anomalies. The 20–40 s phase velocity maps (Figure 7b–d) roughly exhibit similar distributions and are correlated negatively with the Moho depth (Figure 9a) in the study area. Moreover, prominent low-velocity perturbations can be found beneath the Datong volcano (Figure 7 and Figure S1). This low-velocity anomaly extends slightly northwest (Figure 7c,d), which is also observed by the corresponding low S-wave velocity [20,27]. The low-velocity anomalies beneath the northeastern margin of the Tibetan plateau are obvious but not continuous. At long periods of 60–80 s (Figure 7e,f), the Ordos block and the northeastern part of the Alax block show significant high-velocity anomalies. The NW–SE trending low-phase velocities are estimated in the southwestern Zhangbo fault zone (Figure 7e,f and Figure S2). Similar low-velocity anomalies were revealed by Chen et al. [25] and inferred to be related to the upwelling asthenosphere below the Datong volcano.

3.2. Azimuthal Anisotropy

We denote the azimuthal anisotropy of the Rayleigh wave phase velocity with a black error bar, where the length and orientation stand for the anisotropic magnitude and direction of the fast axis, respectively. The statistical direction of the fast axis in each tectonic block, e.g., the Alax block, the Qilian orogen and Songpan–Ganzi terrane, the south Ordos block, the north Ordos block, the Lüliang and Taihang mountain, the Bohai Bay basin, and the Datong volcano region, are represented by the red rose diagram at each period (Figure 7).
Two prominent rotations of azimuthal anisotropy are observed at periods of 15–40 s in the study area (Figure 7a–d). One is centered on the Ordos block, where its surrounding fast axis direction rotated counterclockwise. Another is present under the Datong volcano region (Figure 7c,d and Figure S3). Bounded by the latitude of 38°, the azimuthal anisotropy of Rayleigh wave phase velocity in the north and south side of the Ordos block is obviously different. The direction of the fast axis of the Ordos block is NE–SW in the north, while it turns to NW–SE in the south. At periods of 60–80 s (Figure 7e,f), the fast axis in the Ordos block seems to be parallel to the direction of the SKS/SKKS splitting [24,38] and the APM [39] (Figure 9). In the northeastern margin of the Tibetan plateau, the fast direction is NW–SE and parallel to the strike of the Qilian orogen (Figure 7), which is consistent with the azimuthal fast direction of the Ps and SKS/SKKS phases splitting [24,40]. The direction of the fast axis in the Trans-North China orogen is NE–SW at periods of 10–30 s and parallel to the strike of the Lüliang and Taihang mountains. At periods of 40–80 s, the fast axis of the southern Trans-North China orogen aligns along a NWW–SEE direction that is different from the direction of the northern Trans-North China orogen where the fast axis trends almost N–S at 40 s period and NNW–SSE at 60–80 s periods. The direction of the fast axis is NW–SE in most of the Bohai Bay basin region. Especially, the magnitude of anisotropy in the southwestern Zhangbo fault zone is very strong at periods of 60–80 s (Figure 7e,f).

4. Discussion

The Rayleigh wave phase velocity is most sensitive to the S-wave velocity structure of the Earth’s medium in the depth range of about 1/3 to 1/2 the wavelength, although it is also affected by the temperature and mineral composition of the medium. The sensitivity kernels of Rayleigh wave phase velocity are shown in Figure 8. According to the sensitivity kernel curves, the short-period Rayleigh waves have a narrow sensitive depth range and a high resolution while long-period Rayleigh waves shows the opposite pattern. The depth of the sensitive kernel overlaps with the increase of the period. The Rayleigh wave phase velocity of 80 s is mainly sensitive to the S-wave velocity structure at 100–150 km. Therefore, maps of Rayleigh wave phase velocity and azimuthal anisotropy (Figure 7) in this study can reliably characterize the structure and deformation within the depth range of 150 km.
The cause of azimuth anisotropy of the Rayleigh wave phase velocity is relatively complex. It is generally believed that the shape preferred orientation and the lattice preferred orientation mainly contribute to the seismic anisotropy in the crust and lithospheric mantle [41,42]. The shape preferred orientation is considered to be caused by geological structures such as oriented microcracks, faults, rifts and orogenic belts [41]. The lattice preferred orientation is produced by the alignment of anisotropic minerals due to structural deformation [42]. Under the regional tectonic stress field, a series of faults and orogenic belts will be formed on the surface of the Earth, resulting in the alignment of microcracks in the superficial crust. Rayleigh wave phase velocity azimuthal anisotropy of the upper crust is controlled by the surface tectonic evolution so that the direction of the fast axis is parallel to the strike of faults or orogenic belts [41,43]. In the mid-lower crust and lithospheric mantle, the Rayleigh wave phase velocity azimuthal anisotropy is the “fossil” anisotropy frozen by the preferentially aligned anisotropic minerals during the recent tectonic movement [44,45]. Quartz, mica, amphibole and other anisotropic minerals are widely distributed in the mid-lower crust. Olivine and pyroxene are the main mineral constituents of the upper mantle. These minerals are anisotropic. Especially olivine has the strongest ability of alignment, which is considered to be the main cause of upper mantle anisotropy [24,42]. The Rayleigh wave phase velocity azimuthal anisotropy of the asthenospheric mantle may be related to the flow of mantle material and is usually parallel to the flow direction [46,47]. Therefore, measurements of Rayleigh wave phase velocity azimuthal anisotropy are helpful to study the regional tectonic evolution and interaction among adjacent tectonic blocks and possible perturbations of the mantle flow.

4.1. Lithospheric Thinning and Volcanism of the NCC

Previous studies have shown that the present lithospheric thickness of the NCC gradually decreases from west to east and the lithospheric mantle gradually varies from an old cratonic type to modified and young types [48]. The composition and nature of the Archean mantle lithosphere beneath the eastern NCC has transformed since the Mesozoic, leading to a significant craton rejuvenation [12,48]. This evolution of the mantle lithosphere is thought to cause the geological deformation responses in the crustal and upper mantle, such as the extensional basins and faults in the upper crust, the lithospheric deformation and the intracontinental volcanism [12,13,48].
Most of the eastern NCC is occupied by the Bohai bay basin. The Taihang mountain and the Sulu orogen serve as the boundaries of the Bohai bay basin along its western and eastern edges, respectively (Figure 1). The Bohai bay basin is an extensional basin with several secondary depocenters. The extension of the Bohai bay basin reached a climax between 42 and 32.8 Ma, with NNW–SSE-directed extensional rates of 0.2–0.4 cm/yr [6]. In this study, we observe a low-velocity anomaly at a period of 15 s below the Bohai bay basin that reflects the low-velocity structure of the sedimentary basin. The 15 s fast axis under the Bohai bay basin trends NNW–SSE and is roughly consistent with the direction of the present GPS velocity [9] (Figure 1a, magenta arrows) and the past extensional velocity of the basin [6].
The Trans-North China orogen and western NCC are interpreted to represent the old cratonic stable lithosphere since the Precambrian based on the thick crust [49,50] and lithosphere thickness [51] and low heat flow [52]. However, the recent high-resolution full-waveform images [23], the detailed lithospheric thickness by S-wave receiver functions [15] and the geochemical study on Cenozoic basalts [53] suggest that the lithosphere was thinned in the eastern NCC, and this thinning may also affect the Trans-North China orogen. The rejuvenated craton is generally characterized by a shallow lithosphere–asthenosphere boundary (LAB) and a low-velocity below the lithosphere–asthenosphere boundary depth. The Rayleigh wave phase velocities of 60–80 s correspond to the S-wave velocity structures at depths of 70–150 km (Figure 8). Referring to the lithosphere–asthenosphere boundary depth (Figure 9b), Figure 7e,f demonstrates the structure of the mid-lower lithosphere and the uppermost asthenosphere beneath the NCC. The distinct high-velocity anomalies beneath the Ordos block and the northeastern part of the Alax block indicate that the thick lithospheric root and stable craton are still preserved in the western NCC. In contrast, the low-velocity anomalies on both sides of the North–South Gravity Line beneath the northern Trans-North China orogen support the assumption that the lithosphere would be modified locally beneath the Trans-North China orogen [15,23,53]. Our results provide a new understanding that this lithospheric thinning may occur only in the northern Trans-North China orogen.
The widespread Early Cretaceous magmatism associated with the Paleo-Pacific plate subduction in the NCC, especially in the Bohai bay basin and its surrounding regions, implies that the thinning of the eastern NCC lithosphere mostly took place during this period [11,12,13]. Beneath the Datong volcano, the prominent low-velocity zone indicates the upwelling of hot materials (Figure S1). The asthenospheric upwelling probably originates from the mantle transition zone beneath the Datong volcano based on the prominent low P- and S-wave velocity anomalies from the surface down to approximately 660 km [23,54]. Another important phenomenon is that the upwelling asthenosphere would flow horizontally after encountering the upper lithospheric layer, flowing southeastward into the eastern NCC (Figure 7e,f and Figure S2). The hot materials continue to rise by thermally eroding the overlying lithosphere and extending in the lower crust (Figure 7c,d), driving surface extension and volcanism (Figure 10). This upwelling and its horizontal flow probably causes the observed azimuthal anisotropy [25,55]. The low-velocity anomalies beneath the Datong volcano at periods of 30–40 s (Figure 7c,d) shows the intrusion northwestward in the uppermost mantle with weaker azimuthal anisotropy due to the small range expansion in the strong lithosphere. However, the rotating fast axis around the Datong volcano (Figure 7c,d and Figure S3) agrees well with the results obtained from Ps and SKS/SKKS splitting [22,24], affecting the whole lithosphere, which may represent stress-induced anisotropy caused by the mantle upwelling [22,55]. The low-velocity anomalies of 60–80 s along the southwest side of the Zhangbo fault zone (Figure 7e,f and Figure S2), a noticeable seismic zone with strong azimuthal anisotropy (Figure 7e,f), suggest that the large-scale flow of the upwelling materials enhanced the observed azimuthal anisotropy. Previous studies speculated that the Zhangbo fault zone plays a role as the channel for the upwelling asthenosphere [25,56]. Our results provide new insights into the range of the low-velocity anomaly and the corresponding azimuthal anisotropy, supporting that the upwelling materials of the Datong volcano flow into the low-velocity asthenosphere of the eastern NCC along the southwestern Zhangbo fault zone (Figure 7e,f and Figure S2). Combined with the previous hypotheses [12,13,14,15], that the thinning of the cratonic lithosphere in the eastern NCC resulted from slab rollback related to the subduction of the Paleo-Pacific plate, we further infer that the asthenospheric flow from the mantle upwelling of the Datong volcano also contributed to this thinning (Figure 10).

4.2. Deformation Patterns of the Northeastern Margin of the Tibetan Plateau and the Western NCC

The continental collision between the Indian and Eurasian plates has resulted in the uplift and outward growth of the Tibetan plateau for ~50 Ma [17,18]. Counteracted by the western NCC, the northeastern margin of the Tibetan plateau has a steep topography and strong lithospheric deformation. Thus, it is essential to propose a suitable model for understanding the mechanism of the growth and deformation in the northeastern margin of the Tibetan plateau. The crustal channel flow [5,57] and vertical coherent deformation of the entire lithosphere [58,59] have been debated recently to predict the very different seismic anisotropy and velocity structures in the northeastern margin of the Tibetan plateau [24,40,60,61].
According to the continuous low-velocity layer in the mid-lower crust associated with the observed seismic anisotropy, the crustal channel-flow model may explain the lateral crustal deformation along the southeastern margin of the Tibetan plateau [62,63,64]. In the Songpan–Ganzi terrane, the crustal thickening and deformation were also considered to be driven by the crustal channel flow based on the existence of regional low-velocity zones [60,61]. However, the low-velocity anomalies of Rayleigh wave phase velocity at periods of 30–40 s (Figure 7c,d), sampling the mid-lower crust (Figure 8 and Figure 9), are discontinuous, coinciding with the distribution of S-wave velocity revealed from ambient noise tomography [65]. The discontinuous intracrustal low-velocity anomalies are most likely caused by partial melting and may facilitate the shortening between the Tibetan plateau and the western NCC [65,66,67]. In particular, the coherent crustal and upper mantle anisotropy (Figure 9) estimated from the Ps and SKS/SKKS phases splitting contradicts the crust-mantle decoupling due to crustal channel flow [24,40]. The statistical rose diagrams of the northeastern margin of the Tibetan plateau at periods of 15–80 s (Figure 7) show a NW–SE trending azimuthal anisotropy, which agrees well with the S-wave splitting in the crust and upper mantle (Figure 9). Our results in this study conform to the prediction of a vertical coherent deformation of the entire lithosphere. Moreover, the prominent low-velocity anomalies at 20–40 s periods are not interrupted along the collision belt between the northeastern margin of the Tibetan plateau and the western NCC, and invade partially the interior of the Alax block. The encroachment of the Tibetan plateau depends on the strength of the surrounding blocks. This may explain the observed low-phase velocity anomalies in this study. Therefore, we propose a scenario where the vertical coherent lithospheric deformation in the northeastern margin of the Tibetan plateau is produced by the northeastward expansion of this domain upon the western NCC.

4.3. Rotation of the Ordos Block

The Ordos block located at the western NCC is a stable block but its peripheral areas with high-strain rate and seismicity have been active since the Cenozoic [19,68,69]. The Ordos block is surrounded by four extensional rift systems, e.g., the Weihe graben to the southeast, the Fenhe graben to the east, the Hetao graben to the north and the Yinchuan graben to the northwest, and the Liupan fold and thrust belt to the southwest (Figure 1). Therefore, the Ordos block and its adjacent regions play a significant role in the present tectonic deformation. A kinematic counterclockwise rotation of the Ordos block was proposed to accommodate tectonic driving forces from its surrounding regions [19,24,68].
The recent field investigations [70], GPS observations [68,69] and high-resolution remote sensing imagery [71] reported that a series of right-lateral strike-slip faults, with a rate of about 0.8–2.6 mm/yr, occur along the extensional grabens surrounding the Ordos block. In addition, the earthquake focal mechanism solutions in these grabens support the pattern of superficial crustal deformation mentioned above [8]. The rotated direction of Paleo-Pacific plate subduction (Figure 1a) might be coupled to the crustal deformation pattern that provide the extensional–transtensional stress environment for promoting the counterclockwise rotation of the Ordos block [6,10]. This counterclockwise rotation pattern was also observed by the fast axis directions of the Ps and SKS/SKKS phases splitting (Figure 9) in the mid-lower crust and upper mantle [22,24,40]. It is concluded that the fast axis directions from S-wave splitting are produced by the extrusion of the lithospheric material of the Tibetan plateau [22,24,40]. The azimuthal anisotropy of the Rayleigh wave phase velocity around the Ordos block in the crustal and mantle lithosphere, corresponding to the periods of 15–40 s (Figure 7a–d), is quite consistent with the proposed model of rotation. The low Rayleigh wave phase velocity (Figure 7a,b) as well as the S-wave velocity [27] anomalies along the periphery of the Ordos block are the manifestation of the weak crustal structure that is more likely to produce the transtensional shear deformation inside the crust of the extensional grabens. Our results seem to support the previous hypothesis that the eastward extrusion of the Tibetan plateau and the subduction of the Paleo-Pacific plate caused the shortening rate as high as 3.0 mm/yr in the Liupan thrust-fold belt and extensional rate of 1.3–1.7 mm/yr across the Weihe and Fenhe grabens [69], respectively, and jointly contributed to the rotation of the Ordos block [20,21,22].

5. Conclusions

In this study, Rayleigh wave phase velocity and azimuthal anisotropy at periods of 10–80 s was estimated from the teleseismic data recorded by 388 permanent stations in the NCC and its adjacent areas. Our results reveal the detailed crustal and upper mantle structures of the NCC and its adjacent areas, which provides a new reference for understanding the lithospheric deformation. The prominent low-velocity anomalies with strong NW–SE azimuthal anisotropy from the Datong volcano to the eastern NCC indicate that the lithospheric thinning of the eastern NCC may be also caused by the asthenospheric upwelling of the Datong volcano, except for the subduction of the Paleo-Pacific plate. Our observed coherent crustal and upper-mantle azimuthal anisotropy beneath the northeastern margin of the Tibetan plateau conforms to the prediction of a vertical coherent deformation of the entire lithosphere. Combined with the previous geological and geophysical studies, the eastward extrusion of the Tibetan plateau and the subduction of the Paleo-Pacific plate may contribute to the counterclockwise rotation of the Ordos block.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/rs14010110/s1, Figure S1: Schematic of the asthenospheric upwelling beneath the Datong volcano. (a) Map of the study area, indicating the location of the profile shown in (b). The abbreviations of geological blocks are defined in the caption of Figure 1. (b) Vertical profile of Rayleigh wave phase velocity. Figure S2: Refined Rayleigh wave phase velocity and azimuthal anisotropy at 60 and 80 s periods indicating the horizontal flow of the asthenosphere. The abbreviations of geological units are defined in the caption of Figure 1. The black error bar denotes the magnitude and fast wave direction of azimuthal anisotropy. The average Rayleigh wave phase velocity (Ave.) at each period is shown in the corresponding panel. Figure S3: Refined Rayleigh wave phase velocity and azimuthal anisotropy at 30 and 40 s periods indicating the rotating fast axis around the Datong volcano. The abbreviations of geological blocks are defined in the caption of Figure 1. The black error bar denotes the magnitude and fast wave direction of azimuthal anisotropy. The average Rayleigh wave phase velocity (Ave.) at each period is shown in the corresponding panel. Figure S4: Refined Rayleigh wave phase velocity and azimuthal anisotropy at 30 period indicating the discontinuous low velocities in the northeastern margin of the Tibetan plateau. The abbreviations of geological blocks are defined in the caption of Figure 1. The black error bar denotes the magnitude and fast wave direction of azimuthal anisotropy. The average Rayleigh wave phase velocity (Ave.) at 30 s period is shown in the figure.

Author Contributions

X.X. and D.Z. conceived and designed the paper; X.X. inversed the results and wrote the manuscript; X.H. and X.C. participated in the data processing; D.Z. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China, grant number 41974100.

Data Availability Statement

The model of Rayleigh wave phase velocity and azimuthal anisotropy in current research are available from Xiaoming Xu ([email protected]) upon request.

Acknowledgments

Waveform data for this study are provided by Data Management Centre of China National Seismic Network at Institute of Geophysics (SEISDMC, doi:10.11998/SeisDmc/SN). All figures were plotted by Generic Mapping Tools (GMT).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, D.Y.; Nutman, A.P.; Compston, W.; Wu, J.S.; Shen, Q.H. Remnants of ≥3800 Ma crust in the Chinese part of the Sino-Korean craton. Geology 1992, 20, 339–342. [Google Scholar] [CrossRef]
  2. Zhao, G.; Wilde, S.A.; Cawood, P.A.; Sun, M. Archean blocks and their boundaries in the North China Craton: Lithological, geochemical, structural and P–T path constraints and tectonic evolution. Precambrian Res. 2001, 107, 45–73. [Google Scholar] [CrossRef]
  3. Xiao, W.; Windley, B.F.; Hao, J.; Zhai, M. Accretion leading to collision and the Permian Solonker suture, Inner Mongolia, China: Termination of the central Asian orogenic belt. Tectonics 2003, 22, 1069. [Google Scholar] [CrossRef] [Green Version]
  4. Xiao, W.; Windley, B.F.; Yong, Y.; Yan, Z.; Yuan, C.; Liu, C.-Z.; Li, J. Early Paleozoic to Devonian multiple-accretionary model for the Qilian Shan, NW China. J. Asian Earth Sci. 2009, 35, 323–333. [Google Scholar] [CrossRef]
  5. Clark, M.; Royden, L. Topographic ooze: Building the eastern margin of Tibet by lower crustal flow. Geology 2000, 28, 703–706. [Google Scholar] [CrossRef]
  6. Zhu, Y.; Liu, S.; Zhang, B.; Gurnis, M.; Ma, P. Reconstruction of the Cenozoic deformation of the Bohai Bay Basin, North China. Basin Res. 2020, 33, 364–381. [Google Scholar] [CrossRef]
  7. Liu, J.; Xie, F.; Lv, Y. Seismic hazard assessments for the Ordos Block and its periphery in China. Soil Dyn. Earthq. Eng. 2016, 84, 70–82. [Google Scholar] [CrossRef]
  8. Guo, X.; Jiang, C.; Wang, X.; Tian, X. Characteristics of small to moderate focal mechanism solutions stress field of the circum-Ordos block. J. Geod. Geodyn. 2017, 37, 675–685. [Google Scholar] [CrossRef]
  9. Wang, M.; Shen, Z.-K. Present-Day Crustal Deformation of Continental China Derived From GPS and Its Tectonic Implications. J. Geophys. Res. Solid Earth 2020, 125. [Google Scholar] [CrossRef] [Green Version]
  10. Ren, J.; Tamaki, K.; Li, S.; Junxia, Z. LateMesozoic and Cenozoic rifting and its dynamic setting in Eastern China and adjacent areas. Tectonophysics 2002, 344, 175–205. [Google Scholar] [CrossRef]
  11. Xu, Y.; Chung, S.; Ma, J.; Shi, L. Contrasting Cenozoic Lithospheric Evolution and Architecture in the Western and Eastern Sino-Korean Craton: Constraints from Geochemistry of Basalts and Mantle Xenoliths. J. Geol. 2004, 112, 593–605. [Google Scholar] [CrossRef] [Green Version]
  12. Zhu, R.X.; Xu, Y.G.; Zhu, G.; Zhang, H.F.; Xia, Q.K.; Zheng, T.Y. Destruction of the North China Craton. Sci. China Earth Sci. 2012, 55, 1565–1587. [Google Scholar] [CrossRef]
  13. Wu, F.; Xu, Y.; Zhu, R.; Zhang, G. Thinning and destruction of the cratonic lithosphere: A global perspective. Sci. China Earth Sci. 2014, 57, 2878–2890. [Google Scholar] [CrossRef]
  14. Zheng, Y.; Wu, F. Growth and reworking of cratonic lithosphere. Chin. Sci. Bull. 2009, 54, 3347–3353. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, L. Concordant structural variations from the surface to the base of the upper mantle in the North China Craton and its tectonic implications. Lithos 2010, 120, 96–115. [Google Scholar] [CrossRef]
  16. Xu, X.; Zhao, L.; Wang, K.; Yang, J. Indication from finite-frequency tomography beneath the North China Craton: The heterogeneity of craton destruction. Sci. China Earth Sci. 2018, 61, 1238–1260. [Google Scholar] [CrossRef]
  17. Molnar, P.; Tapponnier, P. Cenozoic Tectonics of Asia: Effects of a Continental Collision: Features of recent continental tectonics in Asia can be interpreted as results of the India-Eurasia collision. Science 1975, 189, 419–426. [Google Scholar] [CrossRef]
  18. Yin, A.; Harrison, T.M. Geologic evolution of the Himalayan–Tibetan orogeny. Annu. Rev. Earth Planet. Sci. 2000, 28, 211–280. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, Y.Q.; Mercier, J.L.; Vergély, P. Extension in the graben systems around the Ordos (China), and its contribution to the extrusion tectonics of south China with respect to Gobi-Mongolia. Tectonophysics 1998, 285, 41–75. [Google Scholar] [CrossRef]
  20. Yao, Z.; Eric, S.; Wang, C.; Ding, Z.; Chen, Y. Asthenospheric upwelling beneath northeastern margin of Ordos Block: Constraints from Rayleigh surface-wave tomography. Tectonophysics 2020, 790, 228548. [Google Scholar] [CrossRef]
  21. Cai, Y.; Wu, J.; Rietbrock, A.; Wang, W.; Fang, L.; Yi, S.; Liu, J. S wave Velocity Structure of the Crust and Upper Mantle Beneath Shanxi Rift, Central North China Craton and its Tectonic Implications. Tectonics 2021, 40, e2020TC006239. [Google Scholar] [CrossRef]
  22. Xu, X.; Ding, Z.; Li, L.; Niu, F. Crustal Anisotropy Beneath the Trans-North China Orogen and its Adjacent Areas From Receiver Functions. Front. Earth Sci. 2021, 9, 753612. [Google Scholar] [CrossRef]
  23. Tao, K.; Grand, S.P.; Niu, F. Seismic Structure of the Upper Mantle Beneath Eastern Asia From Full Waveform Seismic Tomography. Geochem. Geophys. Geosystems 2018, 19, 2732–2763. [Google Scholar] [CrossRef]
  24. Chang, L.J.; Ding, Z.F.; Wang, C.Y. Upper mantle anisotropy and implications beneath the central and western North China and the NE margin of Tibetan Plateau. Chin. J. Geophys. 2021, 64, 114–130. [Google Scholar] [CrossRef]
  25. Chen, H.; Zhu, L.; Ye, Q.; Wang, Q.; Yang, Y.; Zhang, P. Azimuthal anisotropy of the crust and uppermost mantle in northeast North China Craton from inversion of Rayleigh wave phase velocity. Geophys. J. Int. 2015, 202, 624–639. [Google Scholar] [CrossRef]
  26. Zheng, T.; Ding, Z.; Ning, J.; Liu, K.H.; Gao, S.S.; Chang, L.; Kong, F.; Fan, X. Crustal Azimuthal Anisotropy Beneath the Central North China Craton Revealed by Receiver Functions. Geochem. Geophys. Geosystems 2019, 20, 2235–2251. [Google Scholar] [CrossRef]
  27. Huang, X.; Ding, Z.; Ning, J.; Niu, F.; Li, G.; Wang, X.; Xu, X. Sedimentary and crustal velocity structure of Trans-North China Orogen from joint inversion of Rayleigh wave phase velocity and ellipticity and some implication for Syn-rift volcanism. Tectonophysics 2021, 819, 229104. [Google Scholar] [CrossRef]
  28. Gao, Y.; Wu, J.; Fukao, Y.; Shi, Y.; Zhu, A. Shear wave splitting in the crust in North China: Stress, faults and tectonic implications. Geophys. J. Int. 2011, 187, 642–654. [Google Scholar] [CrossRef] [Green Version]
  29. Chang, L.; Wang, C.-Y.; Ding, Z. Upper mantle anisotropy beneath North China from shear wave splitting measurements. Tectonophysics 2012, 522-523, 235–242. [Google Scholar] [CrossRef]
  30. Zheng, X.-F.; Yao, Z.-X.; Liang, J.-H.; Zheng, J. The Role Played and Opportunities Provided by IGP DMC of China National Seismic Network in Wenchuan Earthquake Disaster Relief and Researches. Bull. Seismol. Soc. Am. 2010, 100, 2866–2872. [Google Scholar] [CrossRef]
  31. Schimmel, M.; Stutzmann, E.; Gallart, J. Using instantaneous phase coherence for signal extraction from ambient noise data at a local to a global scale. Geophys. J. Int. 2010, 184, 494–506. [Google Scholar] [CrossRef] [Green Version]
  32. Li, G.; Niu, F.; Yang, Y.; Xie, J. An investigation of time–frequency domain phase-weighted stacking and its application to phase-velocity extraction from ambient noise’s empirical Green’s functions. Geophys. J. Int. 2017, 212, 1143–1156. [Google Scholar] [CrossRef]
  33. Levshin, A.L.; Ritzwoller, M.H. Automated Detection, Extraction, and Measurement of Regional Surface Waves. Pure Appl Geophys. 2001, 158, 1531–1545. [Google Scholar] [CrossRef]
  34. Barmin, M.P.; Ritzwoller, M.H.; Levshin, A.L. A Fast and Reliable Method for Surface Wave Tomography. Pure Appl. Geophys. 2001, 158, 1351–1375. [Google Scholar] [CrossRef]
  35. Smith, M.L.; Dahlen, F.A. The azimuthal dependence of Love and Rayleigh wave propagation in a slightly anisotropic medium. J. Geophys. Res. 1973, 78, 3321–3333. [Google Scholar] [CrossRef]
  36. Yang, Y.; Zheng, Y.; Chen, J.; Zhou, S.; Ceylan, S.; Sandvol, E.; Tilmann, F.; Priestley, K.; Hearn, T.M.; Ni, J.F. Rayleigh wave phase velocity maps of Tibet and the surrounding regions from ambient seismic noise tomography. Geochem. Geophys. Geosystems 2010, 11, Q08010. [Google Scholar] [CrossRef]
  37. Guo, Z.; Chen, Y.J.; Ning, J.; Yang, Y.; Afonso, J.C.; Tang, Y. Seismic evidence of on-going sublithosphere upper mantle convection for intra-plate volcanism in Northeast China. Earth Planet. Sci. Lett. 2016, 433, 31–43. [Google Scholar] [CrossRef]
  38. Zhu, T.; Ma, X. Uppermantle shear-wave splitting measurements in Mainland China: A review. Earth-Sci. Rev. 2020, 212, 103437. [Google Scholar] [CrossRef]
  39. Kreemer, C.; Blewitt, G.; Klein, E.C. A geodetic plate motion and Global Strain Rate Model. Geochem. Geophys. Geosystems 2014, 15, 3849–3889. [Google Scholar] [CrossRef]
  40. Xu, X.; Niu, F.; Ding, Z.; Chen, Q. Complicated crustal deformation beneath the NE margin of the Tibetan plateau and its adjacent areas revealed by multi-station receiver-function gathering. Earth Planet. Sci. Lett. 2018, 497, 204–216. [Google Scholar] [CrossRef]
  41. Crampin, S.; Peacock, S. A review of the current understanding of seismic shear-wave splitting in the Earth’s crust and common fallacies in interpretation. Wave Motion 2008, 45, 675–722. [Google Scholar] [CrossRef]
  42. Silver, P.G. Seismic anisotropy beneath the continents: Probing the depth of geology. Annu. Rev. Earth Planet. Sci. 1996, 24, 385–432. [Google Scholar] [CrossRef]
  43. Crampin, S.; Lovell, J. A decade of shear-wave splitting in the Earth’s crust: What does it mean? What use can we make of it? and what should we do next? Geophys. J. Int. 1991, 107, 387–407. [Google Scholar] [CrossRef] [Green Version]
  44. Savage, M.K. Lower crustal anisotropy or dipping boundaries? Effects on receiver functions and a case study in New Zealand. J. Geophys. Res. 1998, 103, 15069–15087. [Google Scholar] [CrossRef]
  45. Ko, B.; Jung, H. Crystal preferred orientation of an amphibole experimentally deformed by simple shear. Nat. Commun. 2015, 6, 6586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Long, M.D.; Becker, T.W. Mantle dynamics and seismic anisotropy. Earth Planet. Sci. Lett. 2010, 297, 341–354. [Google Scholar] [CrossRef]
  47. Takeo, A.; Kawakatsu, H.; Isse, T.; Nishida, K.; Sugioka, H.; Ito, A.; Shiobara, H.; Suetsugu, D. Seismic azimuthal anisotropy in the oceanic lithosphere and asthenosphere from broadband surface wave analysis of OBS array records at 60 Ma seafloor. J. Geophys. Res. Solid Earth 2016, 121, 1927–1947. [Google Scholar] [CrossRef]
  48. Tang, Y.; Ying, J.; Zhao, Y.; Xu, X. Nature and secular evolution of the lithospheric mantle beneath the North China Craton. Sci. China Earth Sci. 2021, 64, 1492–1503. [Google Scholar] [CrossRef]
  49. Li, S.; Mooney, W.D.; Fan, J. Crustal structure of mainland China from deep seismic sounding data. Tectonophysics 2006, 420, 239–252. [Google Scholar] [CrossRef]
  50. Li, Y.; Gao, M.; Wu, Q. Crustal thickness map of the Chinese mainland from teleseismic receiver functions. Tectonophysics 2014, 611, 51–60. [Google Scholar] [CrossRef]
  51. Huang, Z.-X. Velocity Anisotropy in the Crust and Upper Mantle of North China. Chin. J. Geophys. 2011, 54, 681–691. [Google Scholar] [CrossRef]
  52. Hu, S.; He, L.; Wang, J. Heat flow in the continental area of China: A new data set. Earth Planet. Sci. Lett. 2000, 179, 407–419. [Google Scholar] [CrossRef]
  53. Qian, S.-P.; Ren, Z.-Y.; Richard, W.; Zhang, L.; Zhang, Y.-H.; Hong, L.-B.; Ding, X.-L.; Wu, Y.-D. Petrogenesis of Early Cretaceous basaltic lavas from the North China Craton: Implications for cratonic destruction. J. Geophys. Res. Solid Earth 2017, 122, 1900–1918. [Google Scholar] [CrossRef]
  54. Tian, Y.; Zhao, D.; Sun, R.; Teng, J. Seismic imaging of the crust and upper mantle beneath the North China Craton. Phys. Earth Planet. Inter. 2009, 172, 169–182. [Google Scholar] [CrossRef]
  55. Johnson, J.H.; Savage, M.K.; Townend, J. Distinguishing between stress-induced and structural anisotropy at Mount Ruapehu volcano, New Zealand. J. Geophys. Res. Earth Surf. 2011, 116, B12303. [Google Scholar] [CrossRef]
  56. Huang, J.; Zhao, D. Seismic imaging of the crust and upper mantle under Beijing and surrounding regions. Phys. Earth Planet. Inter. 2009, 173, 330–348. [Google Scholar] [CrossRef]
  57. Royden, L.H.; Burchfiel, B.C.; King, R.W.; Wang, E.; Chen, Z.; Shen, F.; Liu, Y. Surface Deformation and Lower Crustal Flow in Eastern Tibet. Science 1997, 276, 788–790. [Google Scholar] [CrossRef]
  58. England, P.; Houseman, G. Finite strain calculations of continental deformation: 2. Comparison with the India-Asia Collision Zone. J. Geophys. Res. Space Phys. 1986, 91, 3664–3676. [Google Scholar] [CrossRef]
  59. Dewey, J.F.; Shackleton, R.M.; Chengfa, C.; Yiyin, S. The tectonic evolution of the Tibetan Plateau. Philos. Trans. R. Soc. London. Ser. A Math. Phys. Sci. 1988, 327, 379–413. [Google Scholar] [CrossRef]
  60. Sun, A.; Zhao, D. Anisotropic Tomography Beneath Northeast Tibet: Evidence for Regional Crustal Flow. Tectonics 2020, 39, e2020TC006161. [Google Scholar] [CrossRef]
  61. Zhao, P.; Chen, J.; Li, Y.; Liu, Q.; Chen, Y.; Guo, B.; Yin, X. Growth of the Northeastern Tibetan Plateau Driven by Crustal Channel Flow: Evidence From High-Resolution Ambient Noise Imaging. Geophys. Res. Lett. 2021, 48, e2021GL093387. [Google Scholar] [CrossRef]
  62. Yang, Y.; Ritzwoller, M.H.; Zheng, Y.; Shen, W.; Levshin, A.L.; Xie, Z. A synoptic view of the distribution and connectivity of the mid-crustal low velocity zone beneath Tibet. J. Geophys. Res. Space Phys. 2012, 117, B04303. [Google Scholar] [CrossRef]
  63. Bao, X.; Sun, X.; Xu, M.; Eaton, D.W.; Song, X.; Wang, L.; Ding, Z.; Mi, N.; Li, H.; Yu, D. Two crustal low-velocity channels beneath SE Tibet revealed by joint inversion of Rayleigh wave dispersion and receiver functions. Earth Planet. Sci. Lett. 2015, 415, 16–24. [Google Scholar] [CrossRef]
  64. Qiao, L.; Yao, H.; Lai, Y.-C.; Huang, B.-S.; Zhang, P. Crustal Structure of Southwest China and Northern Vietnam From Ambient Noise Tomography: Implication for the Large-Scale Material Transport Model in SE Tibet. Tectonics 2018, 37, 1492–1506. [Google Scholar] [CrossRef]
  65. Li, H.; Shen, Y.; Huang, Z.; Li, X.; Gong, M.; Shi, D.; Sandvol, E.; Li, A. The distribution of the mid-to-lower crustal low-velocity zone beneath the northeastern Tibetan Plateau revealed from ambient noise tomography. J. Geophys. Res. Solid Earth 2014, 119, 1954–1970. [Google Scholar] [CrossRef]
  66. Zheng, D.; Li, H.; Shen, Y.; Tan, J.; Ouyang, L.; Li, X. Crustal and upper mantle structure beneath the northeastern Tibetan Plateau from joint analysis of receiver functions and Rayleigh wave dispersions. Geophys. J. Int. 2016, 204, 583–590. [Google Scholar] [CrossRef] [Green Version]
  67. Li, Y.; Pan, J.; Wu, Q.; Ding, Z. Lithospheric structure beneath the northeastern Tibetan Plateau and the western Sino-Korea Craton revealed by Rayleigh wave tomography. Geophys. J. Int. 2017, 210, 570–584. [Google Scholar] [CrossRef]
  68. Qu, W.; Lu, Z.; Zhang, M.; Zhang, Q.; Wang, Q.; Zhu, W.; Qu, F. Crustal strain fields in the surrounding areas of the Ordos Block, central China, estimated by the least-squares collocation technique. J. Geodyn. 2017, 106, 1–11. [Google Scholar] [CrossRef]
  69. Zhao, B.; Zhang, C.; Wang, D.; Huang, Y.; Tan, K.; Du, R.; Liu, J. Contemporary kinematics of the Ordos block, North China and its adjacent rift systems constrained by dense GPS observations. J. Asian Earth Sci. 2017, 135, 257–267. [Google Scholar] [CrossRef] [Green Version]
  70. Lin, A.; Rao, G.; Yan, B. Flexural fold structures and active faults in the northern–western Weihe Graben, central China. J. Asian Earth Sci. 2015, 114, 226–241. [Google Scholar] [CrossRef] [Green Version]
  71. Rao, G.; Lin, A.; Yan, B.; Jia, D.; Wu, X. Tectonic activity and structural features of active intracontinental normal faults in the Weihe Graben, central China. Tectonophysics 2014, 636, 270–285. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Tectonic setting and seismic station coverage in this study. (a) Schematic diagram showing the collision between the Indian and Eurasian plates, and the subduction of Paleo-Pacific plate (modified after Clark et al. [5] and Zhu et al. [6]). The thick blue rectangle indicates the study area. The bold black lines mark the boundaries of the major tectonic units. The cyan, blue and red circles are the large historic earthquakes (M ≥ 6, [7,8]). The black arrows represent the movements of the different blocks or plates. The magenta arrows show the GPS velocity [9]. WNCC: western North China Craton, TNCO: Trans-North China orogen, ENCC: eastern North China Craton, ZBZ: Zhangbo fault zone, SB: Sichuan basin. (b) Topographic map with the seismic stations used in this study. The thick gray line denotes the North–South Gravity Line (NSGL). The black triangles show the location of the stations. The red triangle is the Datong volcano (DTV). The blue lines indicate the paths of station pair WUL–LZH and BJT–JJ. SGT: Songpan–Ganzi terrane, QLO: Qilian orogen, AB: Alax block, OB: Ordos block, QDO: Qinling–Dabie orogen, SLO: Sulu orogen, LPM: Liupan mountain, YM: Yin mountain, LM: Lüliang mountain, TM: Taihang mountain, YG: Yinchuan graben, HG: Hetao graben, WG: Weihe graben, FG: Fenhe graben, BBB: Bohai bay basin, BS: Bo sea.
Figure 1. Tectonic setting and seismic station coverage in this study. (a) Schematic diagram showing the collision between the Indian and Eurasian plates, and the subduction of Paleo-Pacific plate (modified after Clark et al. [5] and Zhu et al. [6]). The thick blue rectangle indicates the study area. The bold black lines mark the boundaries of the major tectonic units. The cyan, blue and red circles are the large historic earthquakes (M ≥ 6, [7,8]). The black arrows represent the movements of the different blocks or plates. The magenta arrows show the GPS velocity [9]. WNCC: western North China Craton, TNCO: Trans-North China orogen, ENCC: eastern North China Craton, ZBZ: Zhangbo fault zone, SB: Sichuan basin. (b) Topographic map with the seismic stations used in this study. The thick gray line denotes the North–South Gravity Line (NSGL). The black triangles show the location of the stations. The red triangle is the Datong volcano (DTV). The blue lines indicate the paths of station pair WUL–LZH and BJT–JJ. SGT: Songpan–Ganzi terrane, QLO: Qilian orogen, AB: Alax block, OB: Ordos block, QDO: Qinling–Dabie orogen, SLO: Sulu orogen, LPM: Liupan mountain, YM: Yin mountain, LM: Lüliang mountain, TM: Taihang mountain, YG: Yinchuan graben, HG: Hetao graben, WG: Weihe graben, FG: Fenhe graben, BBB: Bohai bay basin, BS: Bo sea.
Remotesensing 14 00110 g001
Figure 2. Distribution of 1347 earthquakes used in this study. Although some events are located inside 20° and outside 120°, all three component waveforms have the epicentral distance between 20° and 120°.
Figure 2. Distribution of 1347 earthquakes used in this study. Although some events are located inside 20° and outside 120°, all three component waveforms have the epicentral distance between 20° and 120°.
Remotesensing 14 00110 g002
Figure 3. (a) Number of paths at each period; (b) Rayleigh wave phase velocity dispersion curves of the station pairs WUL–LZH and JJ–BJT. The paths of station pair WUL–LZH and JJ–BJT are shown in Figure 1b.
Figure 3. (a) Number of paths at each period; (b) Rayleigh wave phase velocity dispersion curves of the station pairs WUL–LZH and JJ–BJT. The paths of station pair WUL–LZH and JJ–BJT are shown in Figure 1b.
Remotesensing 14 00110 g003
Figure 4. The misfit histograms for the final data set of Rayleigh wave phase velocity at periods of (a) 15, (b) 20, (c) 30, (d) 40, (e) 60 and (f) 80 s, respectively. The standard deviation (Std) of each misfit distribution is estimated in each panel.
Figure 4. The misfit histograms for the final data set of Rayleigh wave phase velocity at periods of (a) 15, (b) 20, (c) 30, (d) 40, (e) 60 and (f) 80 s, respectively. The standard deviation (Std) of each misfit distribution is estimated in each panel.
Remotesensing 14 00110 g004
Figure 5. Path coverage for Rayleigh wave phase velocity at periods of (a) 15, (b) 20, (c) 30, (d) 40, (e) 60, and (f) 80 s, respectively.
Figure 5. Path coverage for Rayleigh wave phase velocity at periods of (a) 15, (b) 20, (c) 30, (d) 40, (e) 60, and (f) 80 s, respectively.
Remotesensing 14 00110 g005
Figure 6. Resolution maps of Rayleigh wave tomography at periods of (a) 15, (b) 20, (c) 30, (d) 40, (e) 60 and (f) 80 s, respectively. We mask the area where the resolution is higher than 200 km.
Figure 6. Resolution maps of Rayleigh wave tomography at periods of (a) 15, (b) 20, (c) 30, (d) 40, (e) 60 and (f) 80 s, respectively. We mask the area where the resolution is higher than 200 km.
Remotesensing 14 00110 g006
Figure 7. Examples of Rayleigh wave phase velocity and azimuthal anisotropy at 15, 20, 30, 40, 60 and 80 s periods. The abbreviations of geological units are defined in the caption of Figure 1. The black error bar denotes the magnitude and fast wave direction of azimuthal anisotropy. The average Rayleigh wave phase velocity (Ave.) at each period is shown in the corresponding panel. The red rose diagrams indicate the statistical direction of the fast axis in the Alax block (AB), Qilian orogen and Songpan–Ganzi terrane (QS), south Ordos block (SOB), north Ordos block (NOB), Lüliang and Taihang mountain (LT), Bohai bay basin (BBB) and Datong volcano (DTV).
Figure 7. Examples of Rayleigh wave phase velocity and azimuthal anisotropy at 15, 20, 30, 40, 60 and 80 s periods. The abbreviations of geological units are defined in the caption of Figure 1. The black error bar denotes the magnitude and fast wave direction of azimuthal anisotropy. The average Rayleigh wave phase velocity (Ave.) at each period is shown in the corresponding panel. The red rose diagrams indicate the statistical direction of the fast axis in the Alax block (AB), Qilian orogen and Songpan–Ganzi terrane (QS), south Ordos block (SOB), north Ordos block (NOB), Lüliang and Taihang mountain (LT), Bohai bay basin (BBB) and Datong volcano (DTV).
Remotesensing 14 00110 g007
Figure 8. Sensitivity kernels of Rayleigh wave phase velocity at different periods.
Figure 8. Sensitivity kernels of Rayleigh wave phase velocity at different periods.
Remotesensing 14 00110 g008
Figure 9. Azimuthal anisotropy of Rayleigh wave at period 20 s vs. the Ps phase splitting of receiver function [22,26,40] and at period 40 s vs. the SKS/SKKS phases splitting [24,38] (a) Distribution of Moho depth (after Li et al. [50] and Xu et al. [22,40]). The black arrows indicate the absolute plate motion (APM) direction of the GSRMv2.1 model [39]. The abbreviations of geological units are defined in the caption of Figure 1. (b) Distribution of lithosphere–asthenosphere boundary depth (after Chen et al. [15]).
Figure 9. Azimuthal anisotropy of Rayleigh wave at period 20 s vs. the Ps phase splitting of receiver function [22,26,40] and at period 40 s vs. the SKS/SKKS phases splitting [24,38] (a) Distribution of Moho depth (after Li et al. [50] and Xu et al. [22,40]). The black arrows indicate the absolute plate motion (APM) direction of the GSRMv2.1 model [39]. The abbreviations of geological units are defined in the caption of Figure 1. (b) Distribution of lithosphere–asthenosphere boundary depth (after Chen et al. [15]).
Remotesensing 14 00110 g009
Figure 10. Tectonic sketch map for the lithospheric thinning of the NCC [12,16]. The abbreviations of geological units are defined in the caption of Figure 1.
Figure 10. Tectonic sketch map for the lithospheric thinning of the NCC [12,16]. The abbreviations of geological units are defined in the caption of Figure 1.
Remotesensing 14 00110 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, X.; Zhang, D.; Huang, X.; Cao, X. Deformation of the Crust and Upper Mantle beneath the North China Craton and Its Adjacent Areas Constrained by Rayleigh Wave Phase Velocity and Azimuthal Anisotropy. Remote Sens. 2022, 14, 110. https://doi.org/10.3390/rs14010110

AMA Style

Xu X, Zhang D, Huang X, Cao X. Deformation of the Crust and Upper Mantle beneath the North China Craton and Its Adjacent Areas Constrained by Rayleigh Wave Phase Velocity and Azimuthal Anisotropy. Remote Sensing. 2022; 14(1):110. https://doi.org/10.3390/rs14010110

Chicago/Turabian Style

Xu, Xiaoming, Dazhou Zhang, Xiang Huang, and Xiaoman Cao. 2022. "Deformation of the Crust and Upper Mantle beneath the North China Craton and Its Adjacent Areas Constrained by Rayleigh Wave Phase Velocity and Azimuthal Anisotropy" Remote Sensing 14, no. 1: 110. https://doi.org/10.3390/rs14010110

APA Style

Xu, X., Zhang, D., Huang, X., & Cao, X. (2022). Deformation of the Crust and Upper Mantle beneath the North China Craton and Its Adjacent Areas Constrained by Rayleigh Wave Phase Velocity and Azimuthal Anisotropy. Remote Sensing, 14(1), 110. https://doi.org/10.3390/rs14010110

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop