Deep Electrical Structure of the Langshan Mountain-Linhe Basin Area on the Northwest Edge of the Ordos Block, China

A series of fault depression structures have developed around the Ordos Block. The Langshan Mountain-Linhe Basin area (LLA), located on the northwest edge of the Ordos Block, is a typical, normal tension fault system. A geological survey shows that the Langshan Piedmont fault (LPF) in this area has a large slip rate and indicates risk of earthquake preparation. Broadband magnetotelluric (MT) exploration research was recently carried out across the LLA in the NW–SE direction, and the three-dimensional deep electrical structure thus obtained revealed that the LPF in the LLA is an evident electrical boundary zone on the whole crustal scale and is the main boundary fault of the primary structural block of the Alxa and Ordos Blocks. The MT results also show that the Linhe Basin and Ordos Block belong to the same tectonic basement. The Linhe and Dengkou faults belong to the internal faults of the Ordos Block. The upper crust of the Langshan Mountain on the west side of LPF is characterized by high-resistivity, the middle and lower crust have a low-resistivity layer, and the Linhe Basin on the east side has a Cenozoic low-resistivity sedimentary layer of approximately 10 km thick, which reveals that the Linhe Basin is a faulted basin with sedimentary thickness around the Ordos Block. This indicates that the LLA has experienced continuous and strong tension, normal fault depression sedimentary activities since the Cenozoic era. The current Global Positioning System velocity field shows that there is an apparent NW–SE acceleration zone in the LLA. The leveling data indicate that Linhe Basin shows a subsidence trend relative to the Ordos Block, indicating that the area is undergoing continuous NW–SE tension and faulting. It is speculated that there is a risk of earthquake preparation in the LPF.


INTRODUCTION
There are a series of extensional basin-mountain faulted structural areas around the Ordos Block, such as the Langshan Mountain-Linhe Basin area (LLA) on the northwest edge of the Ordos Block ( Figure 1B), Helanshan Mountain-Yinchuan Basin faulted structural area on the west edge, Qinling Mountain-Weihe Basin faulted structural area on the south edge, and Lvliangshan Mountain-Linfen Basin on the east edge. These basin-fault rift tectonic areas are often strong earthquake preparedness zones for the Chinese mainland (State Seismological Bureau, 1988). In 1739, the Pingluo M8 earthquake occurred in the faulted structural area of Helanshan Mountain-Yinchuan Basin, causing more than 50,000 casualties . The Huashan earthquake occurred in the faulted structural area of the Qinling Mountain-Weihe Basin in 1,556 and was the most serious natural disaster in Chinese historical records (Hou et al., 1998;Rao et al., 2014;Du et al., 2017;Feng et al., 2020). The Hongdong M8 earthquake occurred in the structural area of the Lvliangshan Mountain-Linfen Basin in 1,303 (Deng et al., 2003;Xu et al., 2018). These areas have been the key research areas in the field of Geosciences in recent years. The LLA is located in the northern part of the north-south seismic belt of the Chinese mainland, and it is also subjected to the northeast (NE) pushing of the Tibetan Plateau and northwest-southeast (NW-SE) tensioning caused by deep subduction of the Pacific plate (Northrup et al., 1995;Zhang et al., 2006). The Langshan Piedmont fault (LPF) in this structural area is a large fault that controls the NW boundary of the Hetao Basin (the Linhe Basin belongs to the western region of the Hetao Basin. The Hetao Basin can be divided into three secondary depressions from west to east: Linhe Basin, Baiyanhua Basin, and Huhe Basin), and may also be the boundary fault between the Alxa and Ordos Blocks (State Seismological Bureau, 1988;Rao et al., 2016;Zhang et al., 2020). Geological surveys show that activity of the LPF, especially in the middle section, is very intense (Dong et al., 2016;Rao et al., 2016;Dong et al., 2018a;2018b). Dong et al. (2018a) and Rao et al. (2016) revealed that the fault is an important seismic zone in the Hetao Basin through trenching, dating, and other methods. The average recurrence period of strong earthquakes is about 2,500 A, and the average magnitude of each earthquake reaches 7-8 (Dong et al., 2018a(Dong et al., , 2018bLiang et al., 2021). At present, the recurrence period of strong earthquakes at the fault has been exceeded, so there is a great possibility of destructive, strong earthquakes occurring again (Dong et al., 2018a). Li et al. (2015) and Rao et al. (2016) believe that the latest earthquake rupture event of LPF may be related to the great earthquake in 7 BC. Liang et al. (2021) believe that the magnitude of the earthquake in 7 BC reached 8.1, and there is still a risk of magnitude 7.4-8.0 earthquakes in this area. Due to the lack of large earthquakes recorded by modern instruments in the fault depression structural area of LLA, the existing research results in this area focus on geomorphology and paleoseismology (Jia et al., 2016;Rao et al., 2016;He et al., 2017He et al., , 2018Dong et al., 2018aDong et al., , 2018bDong et al., , 2018cHe et al., 2020;Zhang et al., 2020;Li J. et al., 2021;Liang et al., 2021;Xu et al., 2022). The results obtained lack the support of deep geophysical data. The Linhe Basin has had a relatively developed economy and society since ancient times; there are important towns such as Linhe (LH), Hangjinhouqi (HJHQ), and Wulatehouqi (WLTHQ) with large population densities. Therefore, the deep distribution of faults and seismic risk in the structural area of LLA is of great concern. Wu et al. (2022) recently obtained the large-scale, threedimensional (3D) S-wave velocity structure of the crust and upper mantle of the Ordos Block and its surrounding areas and found that there is a thickening of the Cenozoic sedimentary layer in the Linhe Basin. Chen et al. (2020) revealed that the Moho depth of Alxa Block is 38-45 km, that of Ordos Block is 38-42 km, and that of LLA in the transition zone is 52 km. Seismology also shows that there are obvious differences in the crustal structure between the Alxa and Ordos Blocks . It is speculated that the Alxa and Ordos Blocks belong to different tectonic units. Chang et al. (2011) found that the anisotropic fast wave direction in LLA was in a NW-SE direction when summarizing the anisotropic characteristics of the Ordos Block and surrounding upper mantle, which indicates that there may be NW-SE stretching in this area. Numerical simulation shows that the lateral horizontal movement of the stable Ordos Block has good integrity (Chen et al., 2011). Scholars have carried out Coulomb stress statistics , anisotropy analysis (Hu et al., 2011), and crustal strain analysis  in the LLA on the northwest edge of the Ordos Block and have concluded that the LLA is in a state of tensile strain. The results of multi-scale tomography reveal that the northwest edge of the Ordos Block may have experienced mantle thermal material upwelling, transforming the lithospheric structure (Gao et al., 2018). These research results focus on crustal deformation and large-scale geophysical imaging; there is a lack of research on the fine structure of deep faults, especially fine, deep, geophysical exploration near the LPF.
The electrical structural information of the crust obtained using the magnetotelluric (MT) method can provide a basis for the division of block units (Chen., 1990). It is an important deep geophysical method used to detect the deep extension characteristics of the fault zone and has been widely used (Zhao et al., 2012;Cai et al., 2017;Yu et al., 2017;Wang et al., 2018;Hu et al., 2020;Sun et al., 2020;Li L. et al., 2021;Ye et al., 2021;Zhan et al., 2021;Zhao et al., 2021;Chang et al., 2022;Zhao et al., 2022). Predecessors have completed much MT exploration work in the Ordos Block and its surrounding areas. The Ordos Block belongs to a stable block unit and is in contact with the surrounding blocks (Zhao et al., , 2010Dong et al., 2014;Li et al., 2017;Xu et al., 2017;Zhan et al., 2017;Han et al., 2022). The LLA in the contact area between the Ordos and Alxa Blocks lacks MT exploration results. Therefore, we completed a broadband MT profile with a length of more than 270 km passing through the LLA. The MT profile started from the Alxa Block in the west and extended to the Ordos Block in the east. The data were inversed and used to calculate the deep electrical structural information of the area. From the perspective of electrical structure, this paper discusses the deep structural characteristics of LLA and its adjacent areas, and further evaluates the contact relationship between the Ordos and Alxa Blocks, and then discusses the dynamic background of LLA and the earthquake breeding risk of LPF in combination with the regional 3D deformation field data. The results obtained can advocate the seismic risk assessment of LLA, and can also provide deep geophysical data for the formation mechanism of faulted basins around the Ordos Block.

GEOLOGIC SETTING
LLA is located in the contact area between the Alxa and Ordos Blocks, which is a part of Hetao Graben (Rao et al., 2016;Zhang et al., 2020), The LPF is located at the southeast foot of the NEtrending Langshan Mountain. It intersects with the EW-trending Sertengshan Piedmont fault in the north and terminates in Ulanbuhe Desert in the southwest. The total length of the LPF is approximately 160 km and the overall strike is NE 55°. The overall tendency of the fault is SW (Dong et al., 2018a;2018b). According to the block division theory of Deng et al. (2003), LPF is the boundary fault between the Ordos and Alxa Blocks. The Langshan Mountain on the west side of the fault belongs to the Alxa Block, and the Linhe Basin on the east side of the fault belongs to the Ordos Block (Rao et al., 2016;Zhang et al., 2020). Previous researchers have summarized (State Seismological Bureau, 1988;Dong et al., 2016;Dong et al., 2018aDong et al., , 2018bHe et al., 2020;Liang et al., 2021) the evolutionary history of multiple stages in the LLA: 1) Oligocene stage: This was the initial developmental stage of LLA. Due to the activity of LPF and the relative uplifting of the mountain, the Linhe Basin began to form. The sedimentary thickness of the subsidence center was 2,400 m and the subsidence rate was 0.10 mm/a (State Seismological Bureau, 1988). The maximum thickness of the stratum was mainly distributed along the LPF. 2) Late tertiary stage: This was the most active stage in the development of the LLA. The Linhe Basin continued to decline strongly under the control of the Piedmont faults, with a sedimentary thickness of 3,800 m and a subsidence rate of 0.20 mm/a. 3) In the Pliocene stage, the faulting was more intense, the sedimentary range was further expanded, the thickness reached 6,000 m, and the settlement rate was 0.64 mm/a. 4) In the Quaternary period, the Piedmont Fault activity was still relatively strong, generally faulting the Quaternary system, and the fault activity has maintained its characteristics since the late Tertiary. The mountain rose slowly and the basin continued to settle. The thickness of the Quaternary system was 2,400 m and the settlement rate was 1.00 mm/a. Drilling and petroleum geological profile data (State Seismological Bureau, 1988;Liang et al., 2021) show that Linhe Basin is a sedimentary basin with the largest sedimentary layer thickness around the Ordos Block and is much higher than that of other basins around the Ordos Block over the same period, possibly exceeding 10 km near HJHQ. This shows that the Linhe Basin has a high subsidence rate, which should correspond to very strong fault depression sedimentary activities.

Introduction to the MT Method
The core of the MT method is to determine the underground conductivity structure by using the Earth's external magnetic field and electric field, which vary with time. Specifically, conductivity or resistivity can be determined by the dependence of amplitude and the phase relationship between the electrical field and the magnetic field vector on the Earth's surface. These relationships are expressed in the frequency domain by the correlation transfer function between the surface components of the electromagnetic field, which normalizes the induced response of the magnetic field fluctuation intensity. The transfer function is represented by the impedance tensor, Z, which is defined as: where E is the electric field vector and H is the magnetic field vector (Chen et al., 1990).
In recent years, with the application of supercomputers and the increasingly mature and practical MT 3D inversion algorithms (Siripunvaraporn et al., 2005;Egbert et al., 2012;Kelbert et al., 2014), 3D inversion technology has been widely used in the practical research of various MT detections (Dong et al., 2014;Sun et al., 2019Sun et al., , 2020Zhao et al., 2022).

MT Profile and Data Acquisition
The field data acquisition of LLA was carried out from September to October 2021 and took approximately 35 days. Four sets of MTU-5A equipment (Phoenix Company, Canada) were used for data acquisition at the same time (frequency band range 320-0.0005 Hz). The completed MT profile started from Alxa Block in the border area between China and Mongolia in the northwest, passed through LLA in the southeast, exited the Langshan Mountain near Qingshan (QS), entered Linhe Basin, passed through the Yellow River, entered the hinterland of the Ordos Block, and finally terminated near HJQ ( Figure 1B). The contact area between the Langshan Mountain and the Linhe Basin is the area through which LPF passes, and it is also the key area of this study (as shown by the blue empty frame in Figure 1B). A total of 34 MT sites were obtained along the profile. The profile was~270 km long and densified to 2-3 km near the LPF, which was mainly used to obtain the fine electrical structural information of the LPF area. Due to the relatively simple structure in the peripheral Alxa and Ordos Blocks, the MT site spacing was increased to 10-20 km, which was mainly used to obtain the large-scale, regional electrical structural information ( Figure 1B).

Electrical Orientation Analysis
The phase tensor decomposition technique can be used to obtain the electrical orientation information of the MT profile (Caldwell et al., 2004). The rose-petal diagram of the electrical orientation of the full frequency band and sub-band of all 34 MT sites along the profile was calculated using the phase tensor decomposition method in the MTP software developed by Chen et al. (2004). The rose petal diagram of electrical orientation indicates that the electrical orientation of the full frequency band of all sites along the profile is NE 60°or NW 30° (  Figure 2A). Combined with the tectonic characteristics of the Langshan mountain, LPF, LHF, and DKF in this area, which mainly show a NE 60°orientation, it is reasonable to judge that the electrical orientation of this area is NE 60°on the whole. The rose-petal diagram of the electrical orientation of the subband shows similar characteristics to the full-band, which indicates that the electrical orientation of the study area is relatively uniform from shallow to deep and the linear characteristics are more evident.

Dimensional Analysis
The phase tensor invariant, φ 2 , of each MT site was obtained using the phase tensor decomposition technique (Caldwell et al., 2004). The distribution characteristics of φ 2 with frequency can indicate the variation pattern of the apparent resistivity value below that MT site, which is actually the variation relationship of the apparent resistivity value with depth due to the correspondence between frequency and depth. When φ 2 is less than 45°, the apparent resistivity value below this MT site increases with depth, and when φ 2 is greater than 45°, the apparent resistivity value below this MT site decreases with depth (Heise et al., 2008). Figure 2B shows the distribution of the phase invariance, φ 2 , with frequency for 34 sites along the profile. The electrical variation trends of the plots along the profile are not consistent, the electrical variation trends of the Ordos and Alxa Blocks are relatively simple, and the variation of φ 2 is mainly from high to low and then to high values, which shows that the stratigraphic characteristics of these two blocks are more clear. The value of φ 2 fluctuates more on both sides of the LPF and the high frequency part of the Langshan Mountain is dominated by a low phase, while the middle and low frequency parts are dominated by high frequency, and the φ 2 is larger than 70°. This shows that there may be an evident low-resistivity structure below the Langshan Mountain. In the Linhe Basin east of the LPF, the middle-high frequency part of φ 2 is dominated by high values (black dashed part in the figure), which indicates that there may be an existing, large-scale, low resistivity structure in the shallow part of the area, which corresponds to the existence of a thick, low-resistivity sedimentary layer in the Linhe Basin . The φ 2 in the low-frequency part of the Linhe Basin is dominated by low values, indicating that the resistivity values in the Linhe Basin show a significant increasing trend with increasing depth.
The phase tensor decomposition technique can also obtain the two-dimensional (2D) skewness angle (|β|) with a frequency distribution for each MT site. Considering the interference in the tectonic zone of LLA, it was determined that the subsurface medium could be approximated as a onedimensional (1D) or 2D case when |β| is less than 6, and the subsurface medium can be considered as a 3D case when |β| is greater than 6. The larger the value of |β|, the stronger the three-dimensionality of MT data (Caldwell et al., 2004;Bibby et al., 2005;Cai et al., 2017). The high-frequency part of |β| in the Linhe Basin is smaller (black dashed part in the figure), close to 0, and extends all the way to the vicinity of the Ordos Block, indicating that the shallow part of the region is an approximately 1D or 2D laminar structure, which is presumed to be due to the thicker isotropic Cenozoic sedimentary layers, which is more similar to the distribution characteristics of φ 2 ( Figure 2C). The overall |β| indicates that the Langshan Mountain is more complex than the Linhe Basin, while the smaller values of |β| in the Ordos and Alxa Blocks indicate that the tectonics of these two regions are relatively simple.

3D Inversion
In this study, the resistivity result model of the study area was calculated using the 3D inversion algorithm to invert the original data collected in the field. The latest version of the ModEM program (Egbert et al., 2012;Kelbert et al., 2014) was used and the full impedance tensor (four components: Zxx, Zxy, Zyx, and Zyy) was used as the input data for the 3D inversion calculation. A total  Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 916044 5 of 35 frequency points were used in the 3D inversion calculated from high to low frequency (Supplementary Figure S5). Before inversion calculation, data preparation and unified interpolation were carried out using MTP software  and Topeak software (Liu et al., 2017). The initial model grid was constructed using 3D inversion and was divided into horizontal and vertical grids (Figure 3). The horizontal grid was divided into the core part of the MT sites involved in the calculation and the expanded boundary grid. The horizontal grid of the core part was uniformly dissected using 3.0 × 3.0 km spacing, and a total of 60 (NS) × 60 (EW) grids were obtained. The extended boundary grid was constructed outside the core part by expanding 10 grids in four directions with 1.5 time multipliers. Finally, a total of 80 × 80 grids were constructed in the horizontal direction. The vertical grid design was similar to that of Cai (2017), Sun (2020) and , and a total of 75 layers were obtained in the vertical direction.
In the 3D inversion calculation, the selection of data threshold error was consistent with the previous work of Sun (2020) and . Since the diagonal components (Zxx, Zyy) were too small and the error was too large, we enlarged the error of the diagonal components in the actual 3D inversion and decreased the inversion weighting. The data threshold errors of diagonal components (Zxx, Zyy) and off-diagonal components (Zxy, Zyx) were set to 10% and 5%, respectively. A 100 Ω m uniform half space was used as the initial model and an automatically updated regularization factor was adopted. The initial value of the regularization factor was 5,000. When the inversion was no longer convergent, the regularization factor was updated to one 10th of the original value and the inversion was continued. The inversion was calculated many times and all data were carefully selected. After 65 iterations to the end of the inversion, the root mean square error (RMS) was reduced from 32 to 3.1. The RMS distribution map of each component for all 34 MT sites is shown in Figure 4 and the detailed inversion process is shown in Supplementary Figures S5-S7; Table S1. Figure 5 shows the 3D inversion result within the depth range of 0-60 km along the profile. The topographical changes along the profile and the distribution of blocks and faults are drawn above the electrical structure image. The study area was divided into the Alxa Block, Langshan Mountain, Linhe Basin, and Ordos Block (Dong et al., 2016;Dong et al., 2018aDong et al., , 2018bHe et al., 2020;Liang et al., 2021). Among them, the terrain of the Langshan Mountain fluctuates greatly, with an altitude of 1,500-2,300 m. The terrain of the Linhe Basin is relatively flat, and the average altitude is 1,000 m (Supplementary Figure S2). The maximum altitude difference between the two blocks is up to thousands of meters, which is the location of LPF at the boundary of the two blocks. To show the fine structure of LPF more intuitively and the sedimentary layer thickness of the Linhe Basin, the electrical structure and geomorphic distribution map within 30 km of the core area of LLA are drawn in Figure 6 (the range is shown in the blue dotted box in Figure 1 and Figure 5). The LPF is an obvious electrical boundary zone between the Langshan Mountain and Linhe Basin, whether it is 10 km above or below the area ( Figure 5). It is not only the boundary fault between the Langshan Mountain area and Linhe Basin, but also the largest main boundary fault in the study area. This corresponds to the view of Wen et al. (2014) that the structural deformation in the LLA is mainly distributed on the LPF. Dong et al. (2018a) revealed that the maximum sliding rate and maximum displacement of the fault are located in the middle of the fault, which is also where the MT profile passes through at this time. The upper crust in the Langshan Mountain area west of the LPF is a horizontally layered, high-resistivity body, with a depth of nearly 10 km as a whole, showing relatively complete high-resistivity characteristics, corresponding to the Paleozoic bedrock strata exposed in the Langshan Mountain area above (Rao et al., 2016;Dong et al., 2018a;2018b;Zhang et al., 2020;Liang et al., 2021). The middle and lower crust of the Langshan Mountain area is dominated by a low-resistivity/highconductivity layer (HCL), and there is a clear high and low electrical boundary zone between the middle and lower crust, and the high resistivity body structure of the upper crust at the bottom to 10 km, which indicates that the deep medium characteristics of the Langshan Mountain have a predominantly layered distribution. The HCL in Figures 5, 6 is a large-scale, low resistivity structure in the middle and lower  crust below the Langshan Mountain. Based on other geophysical data, such as the results of multi-scale tomography, we speculate that this may be due to the upwelling of mantle thermal materials and the transformation of the lithospheric structure of Langshan Mountain (Gao et al., 2018). There is a low-resistivity layer in the shallow part of the Linhe Basin to the east of LPF. The distribution of the low-resistivity layer is relatively flat on the profile, and the bottom depth is more than 5 km, becoming shallower and shallower to the east and extending near DKF.

ELECTRICAL STRUCTURE IN THE LANGSHAN MOUNTAIN-LINHE BASIN AREA
To determine the specific location of LPF more accurately, intensive observation was carried out in the contact area between the Langshan Mountain and the Linhe Basin ( Figure 6). It was determined that the accurate location of LPF is between the 2 MT sites, LSL1-16A and LSL1-17B (see Supplementary Figure S4 for the apparent resistivity phase curve of all MT sites and the red triangle in Figure 6). The shallow part of the Linhe Basin is dominated by a low-resistivity layer, combined with 2D deviation angle |β| distribution ( Figure 6). Previous geological survey results (Dong et al., 2016;He et al., 2017; speculate that these low-resistivity layers represent the Cenozoic sediments in the area, with evident 1D or 2D characteristics. The thickness of the low resistivity sedimentary layer near HJHQ in the hinterland of the Linhe Basin exceeds 10 km, which confirms the early geological survey that the Linhe Basin is a sedimentary basin with the largest thickness of sediment around the Ordos Block and that the Cenozoic sedimentary layer may reach 5-6 km (State Seismological Bureau, 1988). The latest research of Liang et al. (2021) found that the sedimentary layer of the Linhe Basin may be close to 10 km in the area crossed by the MT profile, which is similar to our results. Wu et al. (2022) recently obtained the large-scale, three-dimensional (3D) S-wave velocity structure of the crust and upper mantle of the Ordos Block and its surrounding areas, and found that there is a huge thickness of Cenozoic sedimentary layer in the Linhe Basin. The extremely thick Cenozoic sedimentary layer further reflects that the LLA has experienced continuous and strong tension, and normal fault depression sedimentary activities since the Cenozoic stage (State Seismological Bureau, 1988;Dong et al., 2018a;2018b). It is one of the areas with the strongest tension in the Ordos Block. There is an evident high-resistivity body under the low-resistivity sedimentary layer in Linhe Basin. The high-resistivity basement is connected with the Ordos Block and extends to the hinterland of the Ordos Block, indicating that Linhe Basin and Ordos Block may belong to the same structural unit. The Alxa Block to the west of Langshan Mountain has high-resistivity and is a relatively simple structure. Ordos Block was also shown to be a relatively complete block unit, which is similar to the previous MT exploration results (Zhao et al., 2010;Dong et al., 2014;Li et al., 2017;Xu et al., 2017;Han et al., 2022). Through the two north-south MT profiles Li et al., 2017) in Hetao Basin, it was found that the scale of the crust in the north of the Ordos Block is mainly a layered high resistivity body, which corresponds to our MT results. There is a low resistivity layer near the mantle. Our MT results this time did not find a similar situation, which may be due to the short profile. Our results reveal that the electrical structural characteristics of the Ordos Block are similar to those of Han et al. (2022) in the Helanshan Mountain-Yinchuan Basin. In addition, the electrical structural image in Figure 5 shows that LHF and DKF in Linhe Basin are small in scale and do not pass through the high-resistivity body below the Linhe Basin and Ordos Block. These two faults belong to the internal faults of the Ordos Block. It is speculated that it may be an associated fault in the scale of upper crust formed during the process of tension in the structural area of LLA.

DYNAMIC ENLIGHTENMENT OF THE LANGSHAN MOUNTAIN-LINHE BASIN AREA
Active normal faults in China mainly develop around the Ordos Block (State Seismological Bureau, 1988). Except that the southwest edge of the Ordos Block is strike slip thrust, the rest are normal fault basins. Hetao Basin, including Linhe basin, is a sinistral shear tension zone (Fan et al., 2003;Chen et al., 2011). Since the Miocene, the Ordos Block has been characterized by NW-SE tensile stress and NE-SW compressive stress (Northrup et al., 1995;Zhang et al., 2006;Cui et al., 2016). Under the current tectonic stress field environment, the normal fault activity of the NE-distributed LPF coordinates the maximum normal strain produced by the regional stress field. In the active fault zone FIGURE 7 | Distribution map of the global positioning system (GPS) velocity field and leveling field in the study area and adjacent areas (A), and GPS velocity decomposition (B) profile (blue arrow indicates the direction of regional stress field (modified from Zhang et al., 2006;Cui et al., 2016).
Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 916044 around Ordos, the LPF may have a large seismic risk (Dong et al., 2018a(Dong et al., , 2018bLiang et al., 2021). The LPF has been located in the transition area between the Pacific extensional tectonic domain and the Tethys compression tectonic domain since the late Mesozoic (Ren et al., 2002;Zhang et al., 2006). It is a preexisting, large-scale fault system. The study of the activity characteristics and trends of the fault is of great significance to earthquake prevention and disaster reduction in the Linhe Basin and even Hetao Basin (Rao et al., 2016;Dong et al., 2018aDong et al., , 2018bZhang et al., 2020;Liang et al., 2021). Therefore, in addition to using the MT method to obtain the crustal scale electrical structural characteristics of LLA where the LPF is located, this study also used the GPS plane motion data of the study area from 2009 to 2015 (Hao et al., 2019) and the precise leveling vertical motion data from 1970 to 2014  recently obtained by the second monitoring center of the China Earthquake Administration. The combination of static deep electrical structure data and dynamic 3D crustal movement field data provides better constraints for the establishment of geodynamic models in the Linhe Basin and the periphery of the Ordos Block, and provides reliable reference data for seismic risk assessment in the LLA area. Figure 7A shows the GPS velocity distribution field of LLA and its adjacent areas relative to the Eurasian block from 2009 to 2015 (blue plane arrow) and the annual vertical deformation field of LLA and its adjacent areas relative to the Ordos Block obtained from precise leveling data from 1970 to 2014 (color vertical arrow) (data quoted from Hao et al., 2016;. The GPS velocity field within the black box is decomposed into the direction along the MT profile ( Figure 7B; the north and east are positive). The GPS velocity field shows an overall trend of SE movement in the study area, and the GPS velocity field values show an increasing trend from the Alxa Block through the LLA to the Ordos Block, increasing from 4 mm/a to 7 mm/a. This shows that the LLA is an evident accelerating zone, indicating that the area is undergoing continuous NW-SE tension. Current research on tectonic stress fields shows that the fault depression system around the Ordos Block, including the structural area of LLA, shows NW-SE horizontal tensile stress (Zhang et al., 2006;Cui et al., 2016), which corresponds to the results obtained in this study. The precise leveling data from 1970 to 2014 indicate that the Linhe Basin is dominated by uplift (colored vertical arrow in Figure 7A), and the uplift rate is approximately 0-2 mm/a, which is less than the annual uplift rate of 2-4 mm/a in the hinterland of the Ordos Block . The Linhe Basin and Ordos Block belong to the same tectonic unit, while the Linhe Basin shows subsidence relative to the Ordos Block, which may be due to the strong tension in LLA. A movement model map of LLA on the northwest edge of the Ordos Block is postulated ( Figure 8). As the Ordos and Alxa Blocks are currently subjected to the NE pushing of the Tibetan Plateau and the deep subduction and extension of the Pacific plate (Northrup et al., 1995;Zhang et al., 2006), the two blocks continue to carry out evident separation activities, resulting in normal tension fault depression sedimentary activities in LLA, forming a huge fault depression basin. The LPF is the boundary fault between the two blocks and there is concomitantly a strong risk of earthquake preparation.

CONCLUSION
Broadband magnetotelluric exploration was carried out across the LLA in the NW-SE direction, and the 3D, deep electrical structure revealed that the LPF in the LLA was an evident electrical boundary zone on the whole crustal scale and is the main boundary fault of the primary structural block of the Alxa and Ordos Blocks. The MT results also show that Linhe Basin and Ordos Block belong to the same tectonic basement. The Linhe fault (LHF), and Dengkou fault (DKF) belong to the internal faults of the Ordos Block. The upper crust of the Langshan Mountain on the west side of LPF is characterized by high-resistivity, the middle and lower crust have low-resistivity, and the Linhe basin on the east side has a Cenozoic, low-resistivity sedimentary layer of~5-10 km thick. This indicates that the Linhe Basin is a faulted basin with sedimentary thickness near the Ordos Block, reflecting that the LLA has experienced continuous and strong tension and normal fault depression sedimentary activities since the Cenozoic stage. The current GPS velocity field shows that there is an evident NW-SE accelerating zone in the LLA. The leveling data indicate that the Linhe Basin shows a subsidence trend relative to the Ordos Block, indicating that the area is undergoing continuous NW-SE tension and faulting. It is speculated that there is a risk of earthquake preparation in the LPF.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
LZ completed the field data acquisition, the data analyses, and wrote the article. QW provided funding for the study. CC assisted in data acquisition. MH provided GPS and leveling data. All authors contributed to manuscript revision and discussion and approved the submitted version.