Meso-Cenozoic Tectonic Evolution of the Kexueshan Basin, Northwestern Ordos, China: Evidence from Palaeo-Tectonic Stress Fields Analyses

The Kexueshan Basin in the western Ordos Block is characteristic of typical regional superimposed folds in the Jurassic strata, documenting the main tectonic shortening since the Mesozoic. Field structural investigation, fault kinematic analysis, together with regional structural analysis, we reconstruct two-phase tectonic stress fields related to crustal shortening. Synthesizing the regional structural analysis and geochronology, two stages of tectonic evolution are proposed in the Kexueshan Basin since the late Mesozoic. In the mid-late Middle Jurassic (J2), controlled by NW-SE compaction, NW-SE structures formed in the Kexueshan Basin on the western margin of the Ordos Block, marked by NE-to NNE-trending folds and thrust faults. Growth strata found in the study area constrain the starting time of the NW-SE structures. This shortening is likely related to a combination of the southward convergence of the Siberian plate and the northwestward subduction of the Pacific plate. In the end Cenozoic (N1 2-N2), the tectonic stress field was transformed to NE-SW compression, resulting in the development of NW- to NNW-trending fold structures, which were superimposed on early NE-to NNE-trending structures to produce regional superimposed folds. The tectonic shortening coincided with the northeastward growth of the Tibetan Plateau.


INTRODUCTION
As an important tectonic boundary in the western part of the North China Craton (NCC), the western marginal structural belt of the Ordos Block (OB) has attracted the attention of many researchers for a long time (Figure 1). In previous studies, this belt has been named the "Longxi System" and "the spine of the Qilian-Lvliang-Helan mountain-shaped structure" (Li, 1954), the Helan-Liupan fold belt (Ma et al., 1961), and the Ordos platform fold belt (Huang, 1983). Some later studies indicated that this region was a plate suture zone formed in the early Paleozoic, a Meso-Neoproterozoic rift trough (Zhang, 1983;Zhang, 1989), or a "thrust structural belt on the western edge of Ordos" (Chen, 1986;Tang et al., 1988;Liu et al., 2000). Based on the basin-mountain coupling mechanism, combined with geodynamic analysis, Liu et al. (2000) proposed that the southwestern edge of the Ordos Basin was a peripheral foreland basin formed by Late Triassic-Cretaceous arc-continent collision and continental collision orogeny. U-Pb and Hf isotope studies of Middle-Late Devonian detrital zircons in the southeastern part of the Alxa Block and paleomagnetic analysis show that the Indosinian movement caused the Alxa Block to rotate 32°c ounterclockwise relative to the North China Block, followed by convergence and amalgamation to form a Triassic continental collision zone (Yang et al., 2014;Yuan and Yang, 2015). Detailed structural geology and sedimentological analysis showed that the northern part of the western margin of the Ordos Block formed a rift basin in the Late Triassic, while the southern part developed a foreland basin controlled by thrust faults (Liu et al., 2013).  proposed an unrestricted lateral extrusion structure model through the north-south structural difference in the western margin, which was reasonably explained from the perspective of structural geology and dynamic mechanism (Liu and Yang, 1996;Liu and Yang., 1996;Liu and Yang, 1997;Liu et al., 2013).
Sedimentary analysis shows that the western margin of the Ordos Block belonged to a residual inner cratonic basin in the Mesozoic; the foreland basins began to form in the Late Jurassic and further developed in the Early Cretaceous, and the current foreland basin pattern was basically formed by the middle of the Cenozoic (Zhao, 2003;Liu et al., 2005;Liu et al., 2006a;Wang et al., 2006). Darby and Ritts (2002) proposed that the Helan Mountains and Zhuozishan areas on the western margin of the Ordos Block represent an intracontinental deformation zone, in which deformation and structural shortening of at least 30% occurred from the Middle-Early Jurassic to the Late Jurassic. Its tectonic shorting and deformation mechanism may have been related to the subduction of the Paleo-Pacific plate. Zhang et al. (2008) concluded that the northern part of the western margin of the Ordos Block is a large-scale continental margin thrust nappe structural system through regional structural investigation and analysis; this system consists of a Mesozoic sedimentary strata and an early Paleozoic continental margin sedimentary cap layer. The main development stage of the thrust structure was the Jurassic, and the cumulative displacement can reach 60-80 km. A recent analysis of the superimposed deformation of the Helan Mountains shows that the western margin of the Ordos Block experienced two main phases of intracontinental shortening and deformation during the end of Triassic in the NNE-SSW direction and during the Middle Jurassic in the NW-SE direction, forming a regional superimposed fold, which is typical of an intracontinental tectonic deformation zone (Cheng et al., 2019;Huang et al., 2015). These studies show that the western margin of the Ordos Block has experienced complicated intracontinental deformation since the Mesozoic, and the tectonic evolution process remains very controversial. In addition to the controversy over the structures in the western margin of the Ordos Block, the kinematic model and the geodynamic origin of this area still remain unclear due to a lack of detailed structural investigations and the fact that there is no comprehensive understanding of regional structures within the western of the Ordos Block, This paper selects the late Mesozoic basin in the Kexueshan basin on the western margin of the Ordos Block (Figure 1), presents a detailed structural analysis, proposes the main tectonic evolution processes since the Mesozoic, and provides a basis for determining the late Mesozoic-Cenozoic tectonic framework of the western margin of the Ordos Block.
Liumugao fault to the east. The central part of the KB is flat and surrounded by mountains. The Mesozoic sedimentary characteristics and stratigraphic contact relationships in the basin are consistent with those in the Ordos Block (Wang et al., 2006;Zhao et al., 2006). Ordovician, Devonian, Jurassic, Cretaceous and a few Cenozoic strata are exposed in the study area (Figure 1; Geological Bureau and Mineral Resources of Ningxia Hui Autonomous Region, 1983;Region, 1990;Li et al., 1999). The Mesozoic stratigraphic distribution is characterized by the Jurassic strata exposed in the central part of the KB (Ningxia Hui Autonomous Region Geological Survey Institute, 2017), including the Yan'an Formation (J 1 -2 y), Zhiluo Formation (J 2 z), and Anding Formation (J 3 a); the eastern margin is bounded by the Early Cretaceous Miaoshanhu Formation (K 1 ms), and older strata are exposed south, west and north of the KB.
The Jurassic Yan'an Formation (J 1-2 y) is mainly distributed in the Wanbulegou, Tujuan, Shuangjuan, Dazhanchang, Heigeda, and Weigou areas; the lithology is characterized by purple, purple-brown, gray-yellow moderately calcareous breccia and conglomerate at the bottom of the Yan'an Formation, with purple-brown calcareous feldspar quartz sandstone, quartz sandstone, gravel-bearing coarse sandstone, and a coal-bearing seam in the middle of the Yan'an Formation. The Yan'an Formation (J 1-2 y) is the most important coal-bearing layer in the Ordos Basin. In the KB, it mainly comprises sandstones and mudstones intercalated with coal beds (Figure 2). The Xinjing coal mine produces coal mainly from the Yan'an Formation, and the Yan'an Formation can be further divided into five members representing five sedimentary cycles (Liu et al., 2006b). The results of previous studies show that the Yan'an Formation (J 1-2 y) records the complete evolution of a continental lake from its occurrence and development to its shrinkage and extinction, and the middle three members represent the phases of lake expansion . During the deposition of the last member, along with the evident uplift of the Ordos Basin, the lake shrank considerably, and fluvial facies dominated the western margin of the Ordos Block. There are many feldspar quartz sandstones in the Yan'an Formation, indicative of recycled sedimentation. Then, a short-lived sedimentary hiatus occurred, and the upper part of the Yan'an Formation was exhumed and eroded to varying degrees (Liu et al., 2006a).
The Zhiluo Formation (J 2 z) is widely distributed in the study area and is found from north of the Xinjing coal mine to the Dazhanchang. The Zhiluo Formation (J 2 z) and the underlying Yan'an Formation (J 1-2 y) are in parallel unconformity contact, and the Zhiluo Formation (J 2 z) and the Anding Formation (J 3 a) are in conformable contact.
The Zhiluo Formation is a significant uranium-bearing layer in the Ordos Basin (Liu et al., 2006b). Its lithology is comparatively uniform and mainly comprises sandstones, with a small proportion of mudstones. Conglomeratic sandstone is developed locally in the Zhiluo Formation ( Figure 2). Braided rivers occurred in the early phase, and meandering and anastomosed rivers then dominated in the later phase. At its base, a unit of thick-bedded conglomerate was deposited, which can be correlated across the whole basin . Many large lithic fragments are found in the Zhiluo Formation sandstones, including extrusive rocks and quartzite grains. Compared with the Yan'an Formation, the coarser sediments with more lithic fragments in the Zhiluo Formation may indicate tectonic disturbance in the source areas.
The Anding Formation (J 3 a) has a wider distribution range than the Yan'an Formation (J 1-2 y) and Zhiluo Formation (J 2 z), with obvious lithological and lithofacies differentiation from south to north. The sedimentary rocks south of Tujuan are relatively coarse in grain size, and the proportion of coarse clastic rocks is higher than that in the northern sedimentary rocks. The lithologies are mainly purple and gray-brown conglomerate and sandstone, with interposed gray-purple feldspar quartz sandstone. The upper part often contains purple-red calcareous siltstone, argillaceous siltstone and a small number of glutenite lenses, and the sedimentary rocks are mainly purple-red. During its deposition, the color of the sediments changed from gray and grayish green to aubergine and olive, possibly indicating a change in climate . The main sedimentary characteristics are fluvial facies and bank-flood floodplain sedimentary facies. The sediments of the Anding Formation are much finer than those of the Zhiluo Formation.
The Anding Formation (J 3 a) is in conformable contact with the underlying Zhiluo Formation (J 2 z), and the top is in unconformable contact with the clastic rocks of the Lower Cretaceous Miaoshanhu Formation (K 1 ms). The eastern mountain is dominated by Early Cretaceous Miaoshanhu Formation (K 1 ms) strata comprising a set of sedimentary rocks dominated by clastic rocks, with an overall upward thinning sequence structure belonging to alluvial fanlacustrine deposits. The lower part of the Miaoshanhu Formation (K 1 ms) has brown-red, massive medium-coarse conglomerates and coarse boulders, intercalated conglomerate sandstone lenses and a small amount of silty mudstone. The upper part is a light brown-red and light graygreen thick layer of coarse boulders, glutenite and interbedded gray-green mudstone, siltstone, and sandy limestone. The southern and western mountains are composed of the Miboshan Formation (O 2-3 mb) limestone and Abuqiehai Formation (Є 2 -O 1 a) moderately thin microcrystalline limestone with argillaceous bands. Triassic granodiorite intrusions (Tγδ) in the western part of the study area are produced in the form of stock (Bai and Lu, 2009). The pluton is in intrusive contact with the surrounding limestone of the Abuqiehai Formation (Є 2 -O 1 a) and the feldspar quartz sandstone of the Jurassic Yan'an Formation (J 1-2 y). The northern part of the study area is dominated by the Ordovician Majiagou Formation (O 1-2 m) with gray-yellow fragmented crystalline limestone and the Devonian Laojunshan Formation (D 3 l) with gray-purple and graygreen calcareous siltstone. The lower part of this formation is dominated by purple-red medium-thick calcareous conglomerate, and the upper part is gray-purple calcareous siltstone with a small amount of glutenite. The lithological combination is characterized by piedmont accumulationfluvial lacustrine deposits.

Stress Inversion of Fault-Slip Vectors
The reconstruction of the paleotectonic stress field is one of the important topics in the study of structural deformation, providing a basis for regional structural evolution and structural reconstruction (Zhang et al., 2003a;2003b). Interlayer sliding and fold-related fault kinematics, as an important method for the restoration of the paleotectonic stress field in a fold deformation zone, are gradually applied in structural analysis to study the stress mechanism of crustal shortening (Shi et al., 2012;Shi et al., 2013a;Shi et al., 2013b;Shi et al., 2015b). Mesoscale shear faults, notably slickensides, are widespread in geological terrains that are affected by brittle deformation. Shear faults permit the acquisition of so-called fault-slip vector in the field. These data consist of orientations of faults and slip sense or direction from kinematic indicators such as slickensides. By now, fault-slip vectors are almost routinely acquired in structural, mostly neo-tectonic, field studies (Shi et al., 2015b). In this paper, we using the Fault software from Géophysique et Géodynamique Laboratore, Université de Paris-Sud, France (Carey, 1979;Angelier, 1984), yielding three principal stress axes (σ1, σ2, σ3). The reliability and regional significance of a stress field have to be considered, including data collection, separation, calculation so on (see Mercier et al. (1991), and Sperner and Zweigel (2010) for the details). The usage of such data sets for regional tectonic analyses is motivated by the availability of software packages enabling rapid processing of strain increments or states of paleo-stress from fault-slip data, known as fault-slip inversions. To define the formation mechanism of the superposed folds, we reconstructed their paleo-stress fields in this work based on the syn-folding fault kinematic methodology. Better definition of a paleo-stress field can be achieved by calculating fault slip vectors comprising dip directions and angles of fault planes, pitch directions and angles of striations, and strike-slip directions.
In this paper, we dominantly introduce two approaches to determine a paleo-stress direction regarding syn-tectonic shortening by measuring folding-related striations based on the previous studies Shi et al., 2015a). One of represent the growth strata development process and the fold axis change; (D) is the final form of the growth strata, in which the growth axis under the pivot migration mode is always in a constant state, and the growth triangle can be clearly identified, In the wing rotation mode, the growth axis is active and migrates to the front edge, and the growth triangle is not obvious. the approaches is to measure the striations on bedding (S0) that are largely perpendicular to the fold hinge as the result of the mutual slipping along bedding planes within sediments during folding deformation in the study area. Another method is to measure the striations presenting on fault planes and coevally parallel to bedding, which generally indicate the orientation of initial strike-slip movements on fault planes, coincident with the shortening of the original units (see Shi et al., 2015b for the details).

STRUCTURAL ANALYSIS IN THE KEXUESHAN BASIN Structural Deformation and Fault Kinematic Analysis
Relevant structural analysis was carried out in the KB on the western margin of the Ordos Block. The results of regional geological mapping showed that two sets of fold structures have developed in this area, namely, fold structures extending in the NW-SE and NE-SW directions, forming superimposed folds ( Figure 2). Based on the structural maps and field observations, the main superimposed fold structure styles in this area are "T-shaped", "crescent-shaped", "dome-shaped", etc. (Figure 2; Ramsay, 1987). The edge of the KB is mainly controlled by thrust faults. The faults mostly cut the late Mesozoic-Cenozoic strata and locally cut the Cambrian-Ordovician strata, indicating that the fold-thrust deformation has mainly developed since the late Mesozoic. Regional structural analysis shows that most NE-SW-trending thrust faults are cut by NW-SE-trending faults, which roughly indicates that the NE-SW-trending faults in this area were active earlier than the NW-SE-trending faults. The KB is preliminarily confirmed to have experienced two main strong tectonic shortening and deformation events since the late Mesozoic. In this research work, some typical structures in the field were investigated, and structural analysis was carried out from three aspects, namely, NE-SW-oriented structures, NW-SE-oriented structures and superimposed structures, to obtain the twophase structural shortening and deformation events in the KB.

NW-SE Structural Shortening Deformation
A large Jurassic stratum is exposed in the study area, which was subjected to strong tectonic compressive stress, forming a series of NE-SW-trending anticline and syncline structures with various fold styles. According to the size of the limb angle, open folds ( Figure 3E, Figure 4A), medium-normal folds, and closed folds are found in the field; according to the differences in fold axis morphology, there are arc folds, box-shaped folds, and flexural structures (Song et al., 2013). The strata of the southwestern Anding Formation (J 3 a) at Guodikeng in the northern part of the study area are strongly deformed, some strata are nearly vertical or even overturned, and some folds retain a series of simultaneous thrust faults and associated secondary faults ( Figure 4B).
At the northern margin of the KB, the Laojunshan Formation (D 3 l) is unconformably thrust over the gray-green sandstone of the Yan'an Formation (J 1-2 y); the limestone of the Majiagou Formation (O 1-2 m) is thrust over the Yan'an Formation (J 1-2 y) ( Figure 3A), and the structure has been shortened and strongly deformed in this area. The Laojunshan Formation (D 3 l) has developed a fan-shaped compound anticline. The left limb of the fold is an "S"-shaped fold, and the right limb is a "Z"-shaped fold. The fold axis extends in the NE-SW direction. The results of statistical analysis of the fold hinge roughly indicate the effect of NW-SE compression ( Figure 3A). Near Tujuan in the middle of the KB, a flexural structure has developed in the Anding Formation (J 3 a) (observation point K06 in Figure 3B). The gentle strata suddenly became steep, the fold surface is bent The field measurement of limb attitudes yield 9 pivot data, and the polar density maps of the lower hemisphere equal-area stereographic projection are arranged in the NE-SW direction ( Figures 4A,B), indicating the existence of a NW-SE compressive paleostress field. Due to the differences in lithology in various regions and the diverse compressive deformation coefficients of similar rocks, a thrust fault has developed on the right flank of the fold ( Figure 4B) with a fault throw of approximately 1 m. Observation point K73 in the Anding Formation (J 3 a) has developed a recumbent fold with two limbs spaced less than 10 m apart. The inversion results for the fold hinge are scattered in the NE-SW direction also indicates the occurrence of shortening stress in the NW-SE direction. In summary, the analysis of folds with different shapes shows that this area is controlled by compressive stress in the NW-SE direction and is mainly characterized by the development of NE-SW-trending fold structures.
At observation point K05 east of Dazhanchang, the moderately thick sandstone of the Anding Formation (J 3 a) has developed an outcrop-scale thrust fault ( Figure 4C), and fault slip has produced obvious calcite lineations on the fault plane ( Figure 4C). The statistical result for the slip vector indicates the effect of NW-SE compression ( Figure 4C). In the sandstone of the Anding Formation (J 3 a) near Tujuan (observation point K04), "X-shaped" conjugate shear joints have developed ( Figure 4D). The purple-red pebbly sandstone is cut into a diamond-shaped checkerboard, and the two sets of joints are arranged nearly equidistantly; part of the gravel is cut off, the joint surfaces are straight and clean, and the data projection results for the joint surfaces show that joints are controlled by the compressive stress in the NW-SE direction.
In recent years, some researchers have combined growth strata with tectonic chronology to further reveal the time limits of regional tectonic events and constrain the dynamics and kinematics of the basin (Wang et al., 2017a;Liu et al., 2018) and have made good research progress. There are two types of growth strata: fold growth strata and thrust fault growth strata. Growth strata are often found in orogenic foreland basins. In the process of structural deformation, newly deposited strata from nearby provenances form growth folds and growth strata related to thrusting. At present, scholars recognize two main types of deformation mechanism models of growth strata ( Figure 5): a. Hinge migration and b. Limb rotation. The characteristics of growth strata controlled by the two different deformation mechanisms are quite distinct. The characteristics of growth strata formed by hinge migration are as follows: the limb inclination angle of the growth strata remains unchanged, and the fold axis of the growth strata continues to change with structural deformation. Obvious growth triangles can be observed on the vertical section or on the seismic section. In growth strata controlled by the rotation of the limbs, the strata on the side of the rotation are subject to continuous changes in the mechanical direction, resulting in the continuous active structural deformation of the axial plane with tectonic activity ( Figure 5).
The KB is located in the northern part of the thrust tectonic belt on the western margin of the Ordos Block. On the eastern side of Dazhanchang in the KB, synsedimentary structural growth strata developed in the Jurassic Zhiluo Formation (J 2 z) and Anding Formation (J 3 a) at Songjigou. Two profiles were measured in the field, section G-H (Figure 7) and section I-J ( Figure 6). The locations are shown in Figure 2. The two measured sections record the dynamic process of shortening and deformation of the Middle-Late Jurassic structures and the related thrust faults and thrust fold structural styles. Through the analysis of growth strata, the time limits of structural deformation can be well defined (Ford et al., 1997;Zhang P. Z et al., 2006;Liu et al., 2018). Growth strata are a special

NW-SE Structural Shortening Deformation and its Tectonic Stress Field
A series of NW-SE-striking folds and thrust structures are also developed in the Jurassic strata in the study area. This paper analyzes some typical outcrop-scale structures. The moderately thin sandstone of the Anding Formation (J 3 a) from east of Dazhanchang in the northern part of the KB to Guodikeng is strongly deformed, and large-scale vertical folds can be seen in many places, which may indicate superimposed deformation. Observation point K58 ( Figure 8A) has developed superimposed folds. The inclination angle of the fold pivot is between 14°and 21°, and the axial inclination angle is in the range of 47°-70°. According to the classification criteria of the fold position, the fold is described as a plunging and inclined fold (Song et al., 2013). The stereographic projection results for the fold hinge indicate that it was controlled by the stress field of NW-SE compression. The results of the interlayer sliding vector analysis reveal the existence of the NE-SW compression stress field, indicating that the vertical fold was controlled by a twophase compressive stress ( Figure 8A). An outcrop-scale vertically dipping fold ( Figure 8B) can be seen in the sandstone of the Anding Formation (J 3 a) at observation point K91. The lower hemisphere equal-area stereographic projection of the fold hinge data indicates that it is controlled by a NE-SW-directed compressive stress field. At the top of the northern mountain in Dazhanchang (observation point K95), the Anding Formation (J 3 a) exposes a broad and gentle syncline ( Figure 8C). The results for the interlayer slip vector and the stereographic projection of the fold hinge all indicate this NE-SW direction of tectonic compressive stress. Well-developed "X"shaped conjugate shear joints are present in the Anding Formation (J 3 a) near Tujuan in the KB ( Figure 8D). The joint surfaces are straight and smooth, with large inclination angles and long extent. The joint spacing is small and equidistant. The projection results for the joint surfaces show that they are affected by the NE-SW-directed tectonic stress field. At observation point K01 of the Xinjing coal mine in the north, the Ordovician Miboshan Formation (O 2-3 mb) has been thrust onto the Yan'an Formation (J 1-2 y) sandstone; the fault plane has developed fine rough grooves, and the sliding vector projection results indicate that the thrust extrusion in this period was controlled by the NE-SW extrusion stress ( Figure 8E). These structural analyses confirm that strong NE-SW crustal shortening occurred in the KB.

Analysis of Superimposed Deformation
Typical small regional superimposed folds have developed near Dazhanchang in the middle of the KB, showing a relatively complete dome-basin structure (Figure 2, Figure 9A). The latest strata affected by this fold deformation are the Zhiluo Formation (J 2 z) and Anding Formation sandstone (J 3 a).
At observation point K13 ( Figure 9B), small, superimposed syncline basins can be seen spreading in the north-south direction. The graphical interpretation of remote sensing data and the analysis of stratigraphic attitudes show that the basin is caused by shortening in the NW-SE and NE-SW directions, and the topography of superimposed folds resembles a "dumbbell shape" with a clear structure and is called an effective superimposed fold (Ramsay, 1987). NW-SE compression controls the syncline structure oriented in the north-south direction in the basin. Due to shortening in the NE-SW direction, the anticline crosses and overlaps, the dorsal hinge is tilted, and the syncline hinge has risen, forming a saddle structure (Song et al., 2013). At observation point K15, the Zhiluo Formation (J 2 z) thick sandstone is nearly horizontal; two dense and nearly orthogonal vertical cleavages have developed, and partial replacement has occurred. The bedding attitudes are is 033°<70°and 120°<70° (Figures 9B,C), and the bedding attitude is 335°<13°, which indicates that the basin was controlled by crustal shortening in both the NW-SE and NE-SW directions. At observation points K14 and K16 on the western edge of the basin ( Figures 9D,E), a nearly upright cleavage has developed in the Anding Formation (000°<15°), and the attitudes of cleavage are 250°<85°and 248°<70°, indicating the existence of the NE-SW-directed extrusion stress field. Another stage of cleavage developed at observation point K20 ( Figure 9I) in the central part of the superimposed deformed basin, and the attitudes of cleavage are 301°<42°and 330°<59°, indicating that the cleavage in this period was formed by NW-SE structural extrusion. Two sets of diamond-shaped intersecting cleavages have developed at observation points K17, K18, and K19 in the Anding Formation (J 3 a) in the northern part of the basin ( Figures  9F-H); the stratum attitude is 110°<25°, and cleavage planes and stratum planes intersect obliquely, NW-SE compression- related interlayer cleavage attitudes are 301°<42°, 330°<59°, and 332°<60°, and the attitudes of interlayer cleavage related to NE-SW structural shortening are 237°<54°, 240°<60°, and 242°<58°. The two sets of cleavages at observation point K19 ( Figure 9H) have an obvious crosscutting relationship. The cleavage controlled by NW-SE compression is cut by the cleavage controlled by NE-SW crustal shortening, which shows that the structural extrusion deformation in the NE-SW direction occurred later than that in the NW-SE direction. In addition, there are slip vectors at observation points K21 and K23 at the northern and southern ends of the superimposed deformation basin. The inversion results for the slip vectors on the northern margin indicate NW-SE structural compression ( Figure 9J), and the results for the slip vectors on the southern margin indicate NE-SW structural compression ( Figure 9K). A small-scale superimposed deformation syncline has developed at K22 in the Anding Formation (J 3 a) at the northern observation point of the small basin ( Figure 9L). In the early stage, the basic outline of synclinal folds was formed due to extrusion in the NW-SE direction. Later, due to extrusion in the NE-SW direction, the syncline hinge was bent into a wide and gentle "Z" shape, which produced a superimposed fold at the outcrop scale. This shows that the NW-SE tectonic stress field occurred earlier than the NE-SW tectonic stress field.
In summary, the analysis of the deformation of superimposed folds and the crosscutting relationship between the cleavages associated with tectonic compression show that since the Middle-Late Jurassic, the KB has experienced two main periods of structural shortening and deformation oriented NW-SE and NE-SW, and the NW-SE crustal shortening and deformation occurred earlier than the NE-SW deformation; the later NE-SW crustal shortening event basically shaped the tectonic pattern of the KB.

Superimposed Deformation Tectonic Stress Field
Regarding the tectonic stress field in the superimposed fold area, it is necessary to consider the deformation relationship between the earlier and the later fold deformation and fold-related faults. According to the field measurement data, a structure screen is applied to extract the structural deformation patterns at different stages, including the interlayer slip vector, to measure and process tectonic elements at different stages and to restore the tectonic stress field in the superimposed fold area Shi et al., 2012;Shi et al., 2015b). In accordance with the abovementioned ideas and methods, this paper restores the paleotectonic stress field in the superimposed fold deformation zone of the Kexueshan Basin. This work mainly involved fold and fault kinematic analysis of the Mesozoic and Cenozoic strata in the KB, combined with the two phases of structural deformation analysis, and restored the two phases of the tectonic stress field (Figures 10A,B; Table 1).
Observation point K38 at the northern edge of the study area ( Figures 10A,B; Figure 11) is located in the superimposed deformation zone of the Anding Formation in southwestern Guodikeng. The stratum attitude is 145°<25°, and the stratum has not been overturned. There are many secondary calcite stretching lineations developed in the rock layer. The calcite layer is approximately 1-2 cm thick and retains the sliding vector of the shortening deformation that produced the two-stage structure. The lower hemisphere equal-area stereographic projection of the sliding vector retained by the bottom calcite indicates NW-SE tectonic compression, and the overlying late calcite retains the sliding vector formed by NE-SW compression.
Although the early slip vector was destroyed by late tectonic action, it can still be clearly distinguished. According to the crosscutting relationship of the scratch slip vector, the stress field of the NW-SE compression structure in the early stage and the compression structure stress field in the NE-SW direction in the late stage are shown.
Early near-NW-SE compression is found in the Ordovician, Devonian, and Jurassic strata in the study area (Table 1; Figure 10A). The study area has large and wide folds, and due to the uneven local forces on the right limbs of the folds, right limb thrust fault folds developed at Dazhanchang in the Anding Formation ( Figure 4B). Under the control of the early NW-SEtrending paleotectonic stress field, various structural styles, such as flexural folds and box folds, formed in the study area ( Figure 3). Interlayer cleavage and "X"-shaped conjugate shear joints developed in some areas. The abovementioned various structural styles are the result of the NW-SE compressive structural stress field ( Figure 10C). Both the Jurassic Yan'an Formation (J 1-2 y) and Zhiluo Formation (J 2 z) reveal compressive tectonic stress during this period (Table 1; Figure 10A). In summary, the above studies indicate that the starting time of the NW-SE compression paleotectonic stress field was at least after the Middle Jurassic or Late Jurassic. At the same time, the analysis of the statistical results for the slip vector between the folds in the field shows that the NE-SW compression tectonic stress field is well reflected in the strata from the Middle Jurassic Zhiluo Formation (J 2 z) to the end of the Jurassic (Table 1; Figure 10B), forming NW-SE-striking folds, thrust faults, "X"type conjugate joints, and interlayer cleavage and reforming the early folds to produce cross-or diagonal-span superimposed folds. Most importantly, contemporaneous folds developed in the Hongliugou (N 1 h) and Anding Formation (J 3 a) strata (Figure 2), and the pivot direction is oriented NW-SE. Syncline folds are developed in the Hongliugou Formation (N 1 h) northwest of Guodikeng, with the fold axes in the NW-SE direction (Geological Bureau andMineral Resources of Ningxia Hui Autonomous Region, 1983, 1990). This is consistent with fold deformation controlled by the NE-SWtrending paleotectonic stress field ( Figure 10D) in the study area, indicating that NE-SW-trending tectonic compression occurred after the deposition of the Hongliugou (N 1 h) strata in the later period. The Liumugao fault has been an active fault controlled mainly by the NE-SW compressional tectonic stress field since the late Cenozoic (Chen et al., 2013;Gong et al., 2016). After analyzing the two-stage sliding vector lineations in the Jurassic Zhiluo Formation (J2z) sedimentary strata ( Figure 10) and the structural results for the Dazhanchang superimposed deformation basin (Figure 9), the above two lines of evidence support the interpretation that the NW-SE compressional tectonic deformation event occurred earlier than the NE-SW compressional tectonic deformation event. Through statistical analysis of all fold hinges in the study area, the stereographic projection results also indicate the existence of the two-phase paleotectonic stress field ( Figure 10E). Therefore, the above analytical results imply that there have been at least two important tectonic compression and deformation events since the Mesozoic in the KB: the early tectonic stress field was dominated by NW-SE compression ( Figure 10C), and the late tectonic stress field was dominated by NE-SW compression ( Figure 10D).

DISCUSSION: TECTONIC EVOLUTION PROCESS
In summary, the KB on the western margin of the Ordos Block has experienced two main periods of structural shortening events oriented NW-SE and NE-SW. Combining regional structural analysis and deformation-affected stratigraphic characteristics, the two periods of deformation events can be determined to have occurred since the late Mesozoic. This paper proposes a tectonic evolution model of the western margin of the Ordos Block since the late Mesozoic ( Figure 12). The specific analysis is as follows: The Jurassic strata in the Ordos Basin can be divided into five formations, from lower to upper, the Lower Jurassic Fuxian  Figure 2).  (Liu et al., 2006a). Although there existed no sedimentary hiatus between the Zhiluo and Anding formations, U-Pb ages of detrital zircons, lithology, and sedimentary facies show great difference between these two formations (Guo et al., 2017), the sediment colour of the Anding Formation is mainly aubergine and olive and that of the underlying Yan'an and Zhiluo formations is mainly greyish green, the climate of the Anding Formation period was interpreted as cold and dry and that of the Yan'an and Zhiluo formations was warm and humid . All of those indicate that there was a tectonic deformation event in the study area during the Middle Jurassic ( Figure 2C).

Frontiers in Earth
The above structural analysis shows that the early NW-SEtrending tectonic compression affected all the Mesozoic Jurassic strata, forming a regional NE-SW-trending fold in the KB and the basin edge thrust structure involved in the base (Figure 2; Figure 10A). Existing research data indicate that the western margin of the Ordos Block was compressed and uplifted during the Middle-Late Jurassic. The analysis of regional tectonic geology shows that the Ordos Basin and its periphery were compressed in multiple directions under stress field control during the Middle-Late Jurassic . The closure of the Paleo-Tethys Ocean in the southern Ordos Basin resulted in strong intracontinental deformation and a strong northward compressive force. The Mongolia-Okhotsk tectonic belt to the north was in a closure stage, and the Siberian plate was thrust southward, forming a southward compressive force. The eastern circum-Pacific belt, as an active continental margin, formed a strong NW-trending compressive force. Sedimentary analysis shows that the paleocurrent direction on the western margin of the Ordos Basin changed during the Middle Jurassic (J 2 ), indicating a strong intracontinental deformation event (Darby and Ritts, 2002). These studies suggest that the shortening time of the NW-SE-trending structure in the KB should have been in the Middle-Late Jurassic. This period of tectonic shortening and deformation events resulted in a regional angular unconformity contact between the Late Jurassic and Early Cretaceous strata in this area, and thick bottom conglomerates developed (Zhao, 2003;Zhang Y. Q et al., 2006;Zhang Y. Q et al., 2007). The Helan Mountains at the northern margin of the KB experienced superimposed deformation during the Middle-Late Jurassic; that is, the shortening associated with the NW-SEtrending structure was superimposed on the early NWtrending structure (Huang et al., 2015). These studies indicate that the KB also experienced a NW-SE tectonic compression event during this period ( Figure 10A). New progress has been made in the chronological data of the growth strata of the Middle Jurassic Changhangou Formation (J 2 c) in the Daqingshan area on the northern margin of the Ordos Block. The zircon U-Pb ages of the middle-lower volcanic ash interlayers in the growth strata have been obtained, and the structure of this period formed by compression. The onset of compression is limited to the late Middle Jurassic (~170 Ma) (Wang et al., 2017b). The dynamics of this tectonic shortening event are generally attributed to the control of the Middle-Late Jurassic (~165 Ma) East Asian multidirectional convergence system Dong et al., 2008), which may have been related to the joint effect of southward convergence of the Siberian plate and northwestern subduction of the Pacific plate.
Another strong NE-SW-trending tectonic shortening event is identified in the KB, which led to the widespread development of fold structures and thrust faults in the Jurassic strata, and the main structures are NW-SE-trending folds. This period of deformation strongly reshaped the structure related to the early NW-SE compressive stress. In the area, superimposed structural deformation described as "folding", "T-shaped", "crescent" and "dome-shaped" formed, and recumbent folds and overturned folds developed in local areas ( Figure 10B). This tectonic shortening event also affected the Paleogene-Neogene strata on the eastern margin of the study area. Previous studies have shown that a set of stable red clastic sediments in the Eocene-Miocene lake environment developed around the Haiyuan fault zone on the western margin of the Ordos Basin. This indicates that the crust was relatively stable during this period, and the northeastward expansion of the Qinghai-Tibet Plateau did not affect the Haiyuan fault zone on the western margin of the Ordos Basin (Shi et al., 2013c).
A large amount of research shows that since the late Miocene, strong tectonic uplift has occurred on the western margin of the Ordos Block Wang et al., 2011). Lowtemperature thermochronological data also show that a strong NE-SW-trending tectonic shortening event occurred on the western margin of the Ordos Block in the late Miocene (10-8 Ma) Bai and Lu, 2009). From the analysis of the Cenozoic structures on the western margin of the Ordos Block, NE-SW tectonic compression can be seen to have controlled the deposition of the Ganhegou Formation, which unconformably overlies the Qingshuiying Formation (Zhang et al., 2005;Shi et al., 2012); as a result, the sedimentary facies of this group change from underlying riverlake facies to fluvial facies (Wang et al., 2011), and regional basinmountain arc structures formed (Shi et al., 2013a;Shi et al., 2013b;Chen et al., 2015;Shi et al., 2015a). The latest strata involved in the deformation are the late Miocene-Pliocene Ganhegou Formation (Figure 13), and the structural deformation age may be 16.7-5.4 Ma (Wang et al., 2011). Since the end of the Pliocene, against the background of northeastward extrusion of Qinghai-Tibet Plateau material, the Ordos Block has rotated counterclockwise (Zhang et al., 1998). The tectonic stress field in the study area has been mainly transformed to NE compression . A recent Cenozoic paleomagnetic study on the western margin of the Ordos Block determined the age of the Ganhegou Formation to be 9.5-2.7 Ma (Liu et al., 2019), which indicates that the NE-SW-trending structure shortened from the late Miocene to the end of the Pliocene.
In summary, the KB on the western margin of the Ordos Block has experienced two main crustal shortening events. The early NW-SE shortening deformation occurred in the Middle-Late Jurassic, and the late NE-SW shortening deformation occurred Frontiers in Earth Science | www.frontiersin.org March 2022 | Volume 10 | Article 845475 from the late Miocene to the end of the Pliocene. The abovementioned two tectonic deformation events jointly shaped the superimposed fold structure and geomorphic characteristics of the KB. The late dynamic background should be related to the uplift and northeastward material expansion of the Qinghai-Tibet Plateau since the late Cenozoic (Shi et al., 2015a;Shi et al., 2015a).

CONCLUSION
In summary, through analyzing the structural deformation of the KB on the western margin of the Ordos Block and combining this information with the regional structural analysis, the following preliminary interpretations have been obtained: 1) Growth strata were discovered in the Jurassic Zhiluo Formation (J 2 z) and Anding Formation (J 3 a) in the KB; these growth strata indicate that the initial period of early NW-SE structural shortening occurred in the late Middle Jurassic-Late Jurassic, and the structural shortening and deformation during this period controlled the deposition of the Late Jurassic Anding Formation. 2) Since the late Mesozoic Era, the KB has experienced two main phases of intracontinental shortening and deformation. The early period (J 2 ) was compressed in the NW-SE direction, resulting in strong NE-SW-directed fold shortening; the late period (N 1 2 -N 2 ) was controlled by a NE-SW-trending structural compressive stress field, and a series of NW-SEdirected fold structures developed, strongly transforming the early NE-SW structure to form a regional superimposed structure. These results show that the western margin of the Ordos Block is a Meso-Cenozoic intracontinental deformation zone.
3) In terms of the dynamic background of the two phases of crustal shortening and deformation in the KB, the early intracontinental shortening and deformation may have been related to the combined effect of the southward convergence of the Siberian plate and the northwestward subduction of the Pacific plate. The late intracontinental deformation was mainly due to the strong northeastward extrusion of the Qinghai-Tibet Plateau.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.