Next Article in Journal
Radiometric and Polarimetric Quality Validation of Gaofen-3 over a Five-Year Operation Period
Next Article in Special Issue
Remote Sensing Perspectives on Geomorphology and Tectonic Processes
Previous Article in Journal
Study of Land Cover Change in the City with the Fastest Economic Growth in China (Hefei) from 2000 to 2020 Based on Google Earth Engine Platform
Previous Article in Special Issue
DEM Generation from GF-7 Satellite Stereo Imagery Assisted by Space-Borne LiDAR and Its Application to Active Tectonics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quaternary Crustal Shortening of the Houyanshan Structure in the Eastern Chinese Tian Shan: Constrained from Geological and Geomorphological Analyses

1
Guangdong Provincial Key Lab of Geodynamics and Geohazards, School of Earth Sciences and Engineering, Sun Yat-sen University, Guangzhou 510275, China
2
Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519080, China
3
China Institute of Disaster Prevention, Sanhe 065201, China
4
Guangzhou Design Institute Group Co., Ltd., Guangzhou 510620, China
5
School of Civil Engineering and Transportation, South China University of Technology, Guangzhou 510640, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(6), 1603; https://doi.org/10.3390/rs15061603
Submission received: 5 January 2023 / Revised: 17 February 2023 / Accepted: 13 March 2023 / Published: 15 March 2023
(This article belongs to the Special Issue Remote Sensing Perspectives of Geomorphology and Tectonic Processes)

Abstract

:
The Tian Shan is one of the most active intracontinental orogenic belts in the world. It has undergone complex deformation that has resulted in the formation of several fold-and-thrust belts (FTBs) in the piedmonts and intermontane basins. Investigating the deformation histories of these FTBs is important for understanding the tectonic propagation processes of the Tian Shan. Here, we gain insight into these crustal shortening processes by deciphering the Houyanshan structure, a typical fold-thrust belt in the eastern Chinese Tian Shan. We first describe a curved thrust ramp and related fold pairs of the structure using high-resolution remote sensing photography, deformation of fluvial terraces, and field-based geological cross-section. Combined with deformed terrace records and optically stimulated luminescence (OSL) dating results, the kinematic style allows us to yield a geologic shortening rate of 1.6 ± 0.2 mm/a since ~52 ka. Second, to reduce uncertainty in the seismic interpretation and quantify the amount and time of crustal shortening, we interpret three seismic reflection profiles by using the theory of quantitative fault-related fold, area-depth-strain (ADS), and reverse modeling analyses. These profiles provide direct evidence that this structure connects by means of a listric thrust ramp to a shallow detachment level. ADS analysis reveals that the maximum shortening of the Huoyanshan structure is ~4.5 km, which is consistent with the result of quantitative inverse modeling. Each of the structural analysis methods gives similar parameters, and the high consistency of results greatly improves the soundness of a given geologic interpretation. Finally, the shortening rate and total shortening amount suggest that the structure may have formed at 1.8–3.7 Ma, which is nearly synchronous around the Tibetan Plateau. Together, these results indicate that this combined geological and geomorphological analysis provides greater insight into deformation information than can be achieved by any individual technique in studying fold-and-thrust belts worldwide.

1. Introduction

The FTB is a significant deformation feature within compressional regimes [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. With the growth of collisional orogens, deformation progrades and disrupts adjacent foreland basins, which often form several FTBs. These FTBs record pre- and synorogenic deformation history between an orogenic core and its undeformed foreland and meanwhile act as an important coseismic deformation process [15,16,17,18], which is manifested by large magnitude earthquakes that occur within FTBs. These earthquakes have led to tens of thousands of injuries and deaths and huge property losses [19,20,21,22,23,24,25]. As a consequence, knowledge about the development of these structures, including style, shortening estimates, calculation of deformation rate, and the specific timing of initiation of individual structures, not only provides detailed temporal and spatial contexts to improve our understanding of the deformation propagation process in a compressional orogenic belt over geologic time scales but is also of great economic interest because it represents crucial parameters in seismic hazard analysis [26,27,28].
The styles, amounts, and rates as well as the deformation processes of the Tian Shan provide important meaning for studying and understanding the geodynamics of other intracontinental mountains [29]. GPS measurements indicate that the shortening rate across the Tian Shan is ~20 mm/a, accommodating up to 40% of the total present-day plate convergence between India and Asia [30,31,32,33,34], testifying that tectonics of the Tian Shan remain very active. Crustal shortening of the Tian Shan is characterized by a series of approximately E-W trending FTBs within the northern and southern range fronts and intermontane basins [35,36,37,38,39]. These structural belts provide an abundance of deformation records for studying the temporal and spatial evolution of tectonic deformation in the regional range. In recent decades, many earth scientists have investigated the late Cenozoic deformation process of the Tian Shan [35,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55], and efforts have concentrated on constraining the rate, amount, and onset timing of late Cenozoic deformation along the Tian Shan foreland basins that bound the northern and southern piedmonts. In contrast, in the interior, how the shortening is accommodated has not been adequately quantified, especially in the eastern Chinese Tian Shan. Despite the growing number of studies, very few studies have concentrated on this range [37,56,57,58]. Efforts to constrain shortening amounts and rates of individual structures in this range are still needed.
Here, we document the structural deformation across the central part of the Turpan Basin. At the latitude of this range, analyses were focused on the Huoyanshan anticline because it is the most remarkable structure in the Turpan Basin and defines the ∼90 km wide central belt, which likely underwent most of the deformation accommodation of this range [35,59,60]. The Huoyanshan anticline is built by a flight of uplifted terraces in its cut-through valleys. Therefore, surface deformation of the anticline can be documented from the deformation of differential tilting river terraces as well as the well-exposed stratigraphy. Benefiting from the development of oil fields, the subsurface geometry and deformation can be documented from seismic reflection profiles. However, to date, only a few studies were reported on this structure. By using the elevation of the vertical offset of the fault scarp and the corresponding age, the slip rate has been estimated to be 1.8–2.4 mm/a (since 12 ka) in the western part of this structure [37]. The subsurface geometry of the fold, the surface deformation pattern, and the ages of the studied geomorphic markers are barely known. Currently, several new well-imaged seismic lines that reveal the subsurface structure of the fold are available, the surface syntectonic deformation has been constrained by unmanned aerial vehicle (UVA) photography, and the age of the warped terrace across this fold has been better constrained. These new data enable the deformation process of the Huoyanshan anticline to be better estimated. Therefore, it is timely to evaluate how the Huoyanshan fold evolved in the shortening-accommodating stage in the eastern Chinese Tian Shan over the late Quaternary.
Specifically, in this study, a combined approach is used to evaluate the deformation history of the Huoyanshan structure. We document the deformation of fluvial terraces along the Shengjingou Valley in the digital elevation model (DEM) obtained from UVA. Then, we link the surface deformation pattern with bedrock attitudes data to reconstruct the subsurface geometry of the fold. A kinematic model is used to describe the trajectory of thrust propagation and quantitively constrain the deformation history. We describe how these surface records can be used to interpret the subsurface structure. Moreover, based on a new terrace age constraint, a ten thousand years timescale crustal shortening rate of this fold is recalculated in its western segment. Using subsurface high-resolution seismic reflection imagery, the subsurface geometry of the Huoyanshan structure is explored in detail. In addition, ADS analyses are carried out on three seismic lines across the eastern part of the Huoyanshan anticline in order to further constrain the choices of both a kinematic model and the depth of detachment, validate the interpreted detachment depth, and reveal the total shortening. Then, we match the results of seismic interpretation and ADS analysis using the structural reverse modeling method. A quantitative inverse model that is implemented by incorporating trishear in front of a listric thrust and inclined shear in the backlimb accurately restores the structure and gives the optimal solution of the best-fit kinematic model [61,62]. We show that the model-predicted shortenings are consistent with those from ADS analyses. Collectively, these constraints have implications for our understanding of regional tectonism in the Turfan Basin. In conclusion, we discuss the performance of the restoration, the uncertainty of the solution, and the synchronous deformation stage around the Tibetan Plateau.

2. Geologic Setting

2.1. The Tian Shan

The Tian Shan, a mountain belt more than 2500 km long, dominates the topography of central Asia (Figure 1a). From west to east, it extends through the Republics of Uzbekistan, Kazakhstan, Tajikistan, Kyrgyzstan, and China, generally rising to approximately 4000 m in altitude and with peaks exceeding 7000 m (Figure 1b). The Chinese Tien Shan Mountains are divided into two segments along 88°E longitude: the eastern Chinese Tian Shan and the western Chinese Tian Shan. The western Chinese Tian Shan is further divided into the southern Tian Shan, central Tian Shan, and northern Tian Shan [63]. The Tian Shan has undergone complex deformation. The overall range was formed in the late Paleozoic [63,64,65,66,67,68], and was reactivated and intensely uplifted during the Oligocene and Miocene in response to the collision between India and Asia [42,69,70,71,72]. As a consequence of regional S-N trending compression active since the late Cenozoic, deformation propagated basinward and several FTBs composed of subparallel structural belts formed in the range fronts. Different places have different stages of tectonic development. At the northern Tian Shan piedmont, at least three E-W trending structural belts have been formed; in Kuqa, there are four structural belts; Kashi has three structural belts; and in Turfan, there is mainly one structural belt. In general, the present-day Tian Shan topography, composed of a series of ranges with elevations >4000 m separated by large E-W-striking intermontane basins that are bounded by several thrust systems, reflects both its Paleozoic basement architecture and its Cenozoic reactivation history [73,74,75]. According to the results of GPS measurements, present-day geodetic shortening rates decrease from 20 mm/a at ~75° longitude [30,31], to less than 10 mm/a at 84° longitude [76], and then to less than 5 mm/a at 87° longitude [77].
The exact timing of the reactivation of the Tian Shan during the late Cenozoic remains controversial. On the basis of the unconformity at the base of the Oligocene conglomerates in the eastern Tarim and Turfan basins, Allen et al. [63] and Windley et al. [65] suggested that the late Cenozoic uplift of the Tian Shan was induced by the Indian-Asia collision. According to Avouac et al. [35], the uplift of the Tian Shan began in the early to middle Miocene (16 + 22/−9 Ma). Based on sediment flux acceleration detected from changes in the sedimentary logs of drill holes and cross sections from the Tarim and Jungar basins, Metivier and Gaudemer [78] found that sediment flux acceleration began at ~17 Ma and enhanced again at ~5 Ma. Abdrakhmatov et al. [30] estimated the shortening and shortening rate at 76°E and concluded that the deformation initiated at ~10 Ma. However, these pioneering studies contained a degree of uncertainty due to the absence of geologic cross sections across the range and precise chronological constraints. By utilizing low-temperature thermochronology, the uplift history of the Tian Shan has been better constrained. According to fission-track analyses, the unroofing commenced at 25–24 Ma at both the northern and southern margins of the Tian Shan [42,68]. Sobel and Dumitru [71] and Sobel et al. [72] studied sections of the northwestern margin of the Tarim Basin and indicated that exhumation started at ~25–18 Ma, and pulsed events took place along the southern margin of the Tian Shan. Similarly, magnetostratigraphy also provides chronologic analyses of the Cenozoic deformation of the Tian Shan. Bullen et al. [44,79] developed magnetostratigraphy and detrital fission-track studies on the Kyrgyz Tian Shan and discovered pulsed deformation that occurred at ~11 Ma and ~3 Ma. Magnetostratigraphy of late Cenozoic sediments in the Kuqa Depressions, southern Tian Shan piedmont, revealed exhumation ages of 24–20 Ma, and deformation episodes around 16–15 Ma and 7 Ma were also found [41,80]. Based on the magnetostratigraphy of a section in the northern Tian Shan piedmont, Ji et al. [81] concluded that increases in the sedimentation rate at ~26–22 Ma marked the onset of uplift [81]. However, based on increased sedimentation rates, changes in magnetic parameters and thermochronology along the northern and southern Tian Shan range fronts, Charreau et al. [82,83] implied that rock uplift did not begin until ~11 Ma, which is a much younger age. Sun et al. [50,84] commented on a magnetostratigraphy study in the northern Tian Shan piedmont and concluded that the mountain building started at 7 Ma and lasted to ~2.6 Ma. The controversy in these conclusions has demonstrated the ongoing challenges associated with assessing the exact initial uplift. Currently, the deformation study of the Tian Shan is focused on several subparallel folds within the fold-and-thrust belt of the northern and southern margins, which documented the most recent records of the growth of the Tian Shan from the Miocene to the Holocene. Multiple research methods, including morphological and chronologic investigations of deformed geomorphological surfaces, structural analysis, and low-temperature thermochronology as well as magnetostratigraphy of terrestrial sediments of the fold in foreland basins, have been employed to constrain the deformation process [32,39,47,49,85,86,87,88,89,90].

2.2. The Turpan Basin and the Huoyanshan Anticline

The Turpan Basin (Figure 1c), an intermontane basin located in the eastern Chinese Tian Shan, is considered one of the most hydrocarbon-prone areas in the world [91,92]. It extends approximately 500 km from east to west and 60–100 km from south to north and covers an area of approximately 53,000 km2. A minimum altitude of 154 m below sea level is located southwest of the Turfan Basin, known as Aiding Lake, which is the second lowest exposed land surface on Earth. Boundaries of the basin are defined by the Bogda Mountains to the north and the Juelotage Mountains to the south. To the north of the Turpan Basin, seismicity, active deformation slip estimates, and evidence of surface faulting and folding all demonstrate active deformation. This range has a history of infrequent but significant contemporary earthquakes (Figure 1c). The rate of shortening at the northern piedmont of the Bogda Mountains was calculated to be 1 mm/a [58]. At the southern margin of the Bogda Mountains, overthrusting of the Turpan Basin manifests as E-W alignments of active thrust-related folds, where a W-N-W tending FTB can be observed (Figure 1c). This foreland FTB includes (from west to east) the Hongshankou (HSK), Yanshan (YS), Kendeke (KDK), Huoyanshan (HYS), and Qiketai anticlines (QKT). These anticlines generally formed above the detachment in the early Jurassic sediments [37] and are thought to be a long-distance manifestation of the collision of India and Asia [63]. In contrast, Quaternary deformation in the southern margin of the basin is weak, and crustal shortening is not prominent in this range [57,60,93,94].
The Huoyanshan anticline, 90 km long west-east, and 10 km wide north-south, is the most active structure in the basin (Figure 1c and Figure 2). The outcropped Jurassic, Cretaceous, Palaeogene, and Neogene deposits as well as the overlying late Quaternary alluvial terraces all have been warped by folding. The thrust fault lies under the Huoyanshan anticline outcrops at the south margin of the anticline, and the dip angle of the thrust is ~30–50° [59,95]. Based on deformed geomorphological features along the fold, Yang et al. [95] suggested that the Huoyanshan structure accommodates 2.0–3.2 mm/a of crustal shortening. Based on the surface deformation and bedding attitudes, Yang et al. [95] also established a kinematic model and indicated that the Huoyanshan structure is a fault-propagation fold over a listric thrust rooting into a flat detachment layer at depth. However, due to the lack of subsurface structure deformation, it is difficult to further retrieve the complete history of fold growth.

2.3. Stratigraphy

In the Turpan Basin, the sediment fill contains more than 7000 m of continental sediments, extending from the late Permian to the Quaternary, and the Huoyanshan anticline mainly exposes Mesozoic and Cenozoic deposits. From the Permian to the Jurassic, the Jueluotage Mountains were the primary provenance area for the basin, as compared to the locally uplifted Bogda Mountains, which provide little source rocks, forming sedimentary wedges with layers thinning southward. Then, with the gradual uplifting and folding of the Bogda Mountains since the latest Jurassic, the Bogda area built another provenance area for the basin. Until the late Cenozoic caused by the India–Asia collision, the intense uplift of the Bogda Mountains made it the dominant sediment source area [96,97], forming sedimentary wedges with layers thinning northward. Currently, the morphology pattern of stratigraphic units resembles a kind of seesaw (Figure 3).
Figure 4 shows deposited sequences of the foreland FTB. Unconformably above the Permian, the Triassic units are formed of continental sediments with thicknesses ranging from 0.6 to 0.8 km. The Lower to Middle Jurassic units overlie the Triassic units with thicknesses ranging from 1.0 to 1.7 km. The stratigraphic successions of these sections are the Badaowan Formation (J1b), Sangonghe Formation (J1s), Xishanyao Formation (J2x), Sjianfang Formation (J2s), and Qiketai Formation (J2q). A strong, high-amplitude reflection presence at the bottom of the J2x serves as a significant marker for defining the detachment level of the Huoyanshan fold. The Upper Jurassic units, dominated by red and purple-red coarse clastic rocks, crop out on the southern Huoyanshan anticline with a thickness of ~0.6 km, and are divided into the Qigu Formation (J3q) and Kalaza Formation (J3k). Cretaceous deposits overlie unconformably with the Jurassic units and consist of the Sanshilidadun Formation (K1s), Shengjinkou Formation (K1sh), Lianmuqin Formation (K1l), Kumutage Formation (K2k), and Subaishi Formation (K2s). The maximum thickness of the Cretaceous is ~1.0 km. Paleogene and Neogene units outcrop in the northern part of the Huoyanshan anticline and are composed of the Shanshan Group (Esh), Taoshuyuan Formation (E3–N1t), and Putaoyuan Formation (N2p), with thicknesses of 1.6–2.2 km. The primary components of the Paleogene deposits are coarse clastic sediments, while those of the Neogene are fine clasts. The Quaternary deposits are just above the Neogene, including the Xiyu Formation (Q1x) and Middle Pleistocene to Holocene sediments (Q2–4), which are 0.2 to 1.0 km thick.

3. Materials and Methods

3.1. Geomorphology and Quaternary Chronology

Deformed fluvial terraces preserved along the Shengjingou Valley have recorded localized folding deformation across the Huoyanshan anticline [37,56,95]. In this study, the surficial features of these fluvial terraces have been extensively investigated through mapping and GPS surveying. High-resolution DEMs (~1 m) across the Shengjingou Valley were acquired by using UAV photography to correlate surface geomorphology across the entire transect of the Huoyanshan anticline (Figure 5). A total of 5112 photos were captured, and 75 ground control points were placed by using a real-time kinematic GPS (RTK-GPS) with centimeter vertical and horizontal precision (Figure 5) that covered all surveying areas to further improve the precision of the DEM. Details of this procedure can be found in Yang et al. [95]. Based on DEM interpretation and field surveying, we classified and mapped eight generations of terraces across the Shengjingou Valley (T1–T8, from young to old, Figure 6). Then, to obtain the spatial distribution of terraces and riverbeds, profiles along the tops of terrace treads and riverbeds coupled with the coordinates and elevations of survey points (Figure 7a) were extracted from the DEM. To erase the riverbed gradient from the terrace profiles, elevations above the riverbed of the terraces were recalculated by subtracting the riverbed gradient using a polynomial that closely approximates the current riverbed profile (Figure 7b). To reflect the vertical uplift caused by folding, the undeformed terrace elevation was removed from the maximum terrace elevation (Figure 8). Table 1 displays our evaluation of the uncertainties associated with these measures.
OSL dating of deformed surface sediments is proven suitable for constraining the chronology, yielding reliable absolute ages for terrace abandonment through the Holocene, under conditions of sufficient exposure of material being dated before deposition. We collected our sample SJK-T4 from a layer of interbedded silts and sands lying ~3 m beneath outwash gravel on a prominent terrace (T4) in the middle part of the Shengjingou Valley (Figure 5a). Two cores were collected by driving 50 cm aluminum tubes into the fluvial deposit, and light shielding materials were capped on the ends. Samples were processed at the Research Laboratory of Luminescence Dating at the Institute of Geology, China Earthquake Administration. Sample preparation and analysis followed the standard procedures of Aitken [100,101]. The equivalent dose (De) was obtained from the Simple-aliquot Regenerative Dose (SMAR) protocol of fine particle quartz samples (4~11 μm). The environmental dose rate was determined through activity concentration measurements of the U, Th, and K. The dose rate was computed from the measured activity concentration and attenuation values [100]. Both the cosmic ray contribution and the sample water content were used to assess the environmental dose [102]. The simple and weighted means and standard errors are at a 68% confidence level (Table 1).

3.2. Structural Interpretation from Outcrops and Seismic Profiles

The subsurface geometry of the Huoyanshan structure is constrained by the bedrock geology and seismic reflection profiles. Outcrops of the folded strata are exposed by the entrenchment along the Shengjingou Valley. Thus, 21 bedding attitudes were measured across the Shengjingou Valley (Figure 2b and Figure 6). The results reveal a very uniform strike and a continuous and progressive variation in dip (Figure 8), indicating fold growth by limb rotation over listric faults. We acquired 2D pre-stack time migrated seismic reflection data along 3 lines across the eastern part of the Huoyanshan anticline (Figure 2a), which constrains the subsurface geometries to ~7 km below sea level in detail (Figure 9, Figure 10 and Figure 11). These seismic profiles were converted into depth by utilizing a multilayer velocity model [103] in the Move 2013 software (© 2011 Midland Valley Exploration Ltd., 144 West George Street, Glasgow, G2 2HG, UK). Field observations, regional geologic maps, and well data were used to constrain boundaries of horizons in seismic profiles. Several prominent reflections in the seismic sections, including the bottoms of lithological units J2x, J2q, J3k, K1s, Esh, N1t, N2p, and Q1x, can be readily traced and correlated, allowing a robust structural interpretation and analysis with little uncertainty. The kinematic model was inferred from the characteristics of surface and subsurface geometries. The crustal shortening was calculated using the trigonometric relationship in the kinematic model (Figure 8b), and when paired with the terrace abandonment age, the shortening rate was computed (Table 1).

3.3. ADS Analysis

Area-depth-strain (ADS) analysis is a very useful technique for analyzing structural interpretations in contractional systems [104,105,106,107,108,109,110,111,112,113]. ADS analysis provides a simple method of quantitatively constraining the structural evolution history that would not be apparent from the interpretation alone, which in turn validates interpretations [114,115]. Based on the assumption of a constant change in the area during deformation, this method allows the determination of the magnitude of horizontal shortening for growth and pre-growth strata and detachment depth independent of kinematic assumptions. These parameters are determined by measuring the excess area (s) (against the regional level) and the depths (h) (against the reference level) of multiple marker horizons and plotting the corresponding area-depth graph. In this study, ADS analyses followed the method of Wang et al. [111], a generalized ADS method that can be used to handle contractional systems with strata thinning in one direction. The first step was loading three interpreted seismic sections into the program Move and then defining the boundaries where the marker horizons remain undeformed as the pin line and the moving line (regional limits) in each seismic section. Excess areas and depths against the reference level (zero-area line) were measured for each horizon within the program. In the area-depth diagram, x-axis represents the depth against the reference level, y axis represents the excess area, and each horizon has an area-depth value. Balanced interpretations within pre-growth strata yield lines with good linear trends in the area-depth diagram, and the slope of the line equals the overall boundary shortening. The x-intercept equals the distance of the predicted detachment against the reference level. Positive values indicate upward change, while negative values indicate downward change. Note that if the trend of the undeformed marker between two regional limits is not horizontal, the detachment depths at the two regional limits are distinct. Together, the ADS analysis method allows us to constrain the detachment depth and the overall shortening amount more precisely and meanwhile facilitate assessing the structural integrity in the seismic interpretation and geological modeling (Figure 9, Figure 10 and Figure 11) [108,109].

3.4. Inverse Modeling

We adopted a listric fault-related fold kinematic model, including folding over a listric fault connected to a detachment layer at depth. Attention is only paid to the backlimb deformation in this model because this area is often incompletely imaged in seismic lines and it is difficult to model folding in front of propagating thrusts [116]. However, this problem can be solved by using trishear modeling [7,117,118], as it has the ability to explain many of the details of forelimb deformation and mimic those complexities [119]. Thus, in our study, a symmetric trishear was implemented in front of the listric fault to further optimize the kinematic model. The trajectory of listric fault propagation was determined by using the velocity description as described by Cardozo and Brandenburg [62]. The listric fault was modeled as a circular arc with a given center of curvature (CC), the radius of curvature (RC), and the maximum central angle (θmax). In the backlimb, parallel shear takes place and results in constant fault slip behind the trishear zone as well as backlimb geometry concentric to the listric fault. In the forelimb, linear and symmetric trishear operates in front of the fault tip [62,117], and two important parameters, the fault propagation to fault slip ratio (P/S) and trishear angle, mainly control the geometry of the forelimb (Figure 9, Figure 10 and Figure 11). Several model parameters are relatively easy to determine if the final geometry of the structure is well interpreted, such as the CC, RC, θmax, and fault silp (S). However, the quantification of the trishear angle and P/S ratio has proven elusive [119], because trishear is an incremental rather than a graphical model and there is no mathematical method to quantify these two parameters from the final fold geometry [62]. The difficulty of defining these two parameters has been one of the greatest barriers to the widespread acceptance of trishear. To obtain the trishear angle and P/S ratio, evaluate the fitness of a given model parameter combination, and search for the optimal model parameter combination, the primary vehicle is to carry out inverse modeling. The inversion was carried out by using a script offered by Cardozo (https://thecdem.com/ (accessed on 1 March 2022)) that utilizes the simulated annealing MATLAB function (simulannealbnd in the MATLAB Global Optimization Toolbox) [120].
This inversion method is much faster in comparison to the standard grid search inversion proposed by Allmendinger [118]. Additionally, rather than providing a single solution, the algorithm gives a range of models that can fit the structure (Figure 12, Figure 13 and Figure 14) [62].

4. Results

4.1. Late Pleistocene Shortening Rate

4.1.1. Surface Geometry and OSL Age Constraint

Along the Shengjingou Valley, eight fluvial terraces were identified, T1 (youngest) to T8 (oldest) (Figure 5 and Figure 6). Figure 7 shows the profiles of the terrace surfaces and riverbeds. The upper older suite (T5–T8) is not very continuous and is preserved primarily in both the terminal parts of the valley. T5–T8 are significantly elevated (~80–200 m) above the riverbed (Figure 7b). Among the lower set of terraces (T1–T4, <100 m above the modern riverbed), T4 is the most continuous terrace and forms an extensive surface along the valley; therefore, we focus on this prominent terrace (T4) in this study. The vertical uplift (subtracted undeformed terrace surface) of this terrace at the peak of the fold is ~56 m (Table 1). In general, terraces in the Shengjingou Valley are syntectonic sediments of folding. These geomorphic markers provide the most recent record of the folding deformation and fault slip activity of the Huoyangshan structure. Consistent and progressive rotation of the surfaces resulted in no parallelism between the younger and older terraces, which is consistent with backlimb deformation resulting from listric thrust-related folding. Further evidence from the outcrop of bedrock along the Shengjingou Valley is consistent with this observation. The attitudes of bedding along the Shengjingou Valley reveal a very uniform strike and a continuous and progressive shift in the dip. The gentler backlimb dips from ~40° to 9° with northeast dip directions and no dip domain (Figure 8a). This characteristic also implies fold growth by limb rotation over listric faults.
Sample SJK-T4 was collected from fine silts within a 4.5 m section of fluvial gravel cap on the base rock of terrace T4 at the middle of the Shengjingou Valley (Figure 5a), which yields an age of 51.7 ± 5.1 ka. The position of this sample lies at the bottom of the outwash gravel section, indicating that this age is a maximum estimate of the time of terrace abandonment (Table 1).

4.1.2. Subsurface Geometry

Three seismic reflection profiles, Line 1, Line 2, and Line 3, located in the eastern segment of the Huoyanshan anticline (Figure 2a), were used to investigate the subsurface geological structure. These seismic lines have similar geological structures and sediment properties. In all three seismic reflection profiles, above ~4.5 km depth, on the northern part of the profiles, we observed fairly continuous reflectors that are subhorizontal. In the central part, these reflectors gently dip northward and are subhorizontal again in the southern part (Figure 9, Figure 10 and Figure 11). The detachment depth was estimated from the transition between undeformed and folded reflectors (a better constraint of the detachment depth will be presented in Section 4.2). The fault geometry was determined by the discontinuities of the reflectors, delineating a listric thrust boundary of the fold. These geometries outline that the Huoyanshan anticline is an asymmetric fold controlled by a listric thrust rooted in a subhorizontal detachment level. In great detail, as expected from the listric thrust-related fold model characterized by limb rotation, smoothly continuous progressively rotated backlimb concentric to the listric thrust. The area near the surface is incompletely imaged; therefore, it is difficult to estimate the exact position of the fault tip. The thrust was interpreted to be modeled with a circular geometry, with a center (CC) and a radius (RC) of curvature (Figure 9, Figure 10 and Figure 11). The interpreted maximum central angles are 55° in Line 1, 54° in Line 2, and 55° in Line 3. Beyond the maximum central angles, the interpreted thrusts are planar and outcrop the ground with dip angle values the same as the maximum central angle values. In Section 4.3, the position of the fault tip in each profile will be expanded on through restorations of deformed reflectors. The locations of axial surfaces were identified by the transitions from the listric ramp to the upper planar ramp and detachment and were used to define the southern and northern fold limits, respectively. Collectively, the seismic interpretation results show direct and robust evidence that the Huoyanshan anticline is a listric fault-related fold, and that the thrust is mainly composed of a listric fault that connects the upper planar ramp with a detachment located ~4.8 km under the surface (near the base of J2x). This interpretation is consistent with a previous study utilizing the combination of surface terrace deformation and bedding attitudes to infer the kinematics of the Huoyanshan structure [95].

4.1.3. Kinematic Model and Shortening Rate

The features of both the surface and subsurface geometries of the Huoyanshan structure are inconsistent with the traditional fault-related folding models [7,116,121,122,123,124,125]. We adopted a listric thrust folding theory for the Huoyanshan anticline in this research [98,99,126]. In our listric thrust-related folding model, the thrust grows listrically by limb rotation and connects the upper planar ramp with a detachment located at depth (Figure 8b), and the dip angle of the planar segment equals the maximum central angle values. Shortening (S) was estimated from the dip angle of the thrust at the surface (θ) and the maximal vertical uplift (dz) of the terrace [95,99]. The trigonometric relationship between these parameters is as follows:
S = dz   sin θ
This theoretical geometric relationship allows estimation of the horizontal shortenings associated with the different vertical uplifts and the surface fault dip angle (42°). The model-predicted horizontal shortening is 81.7 ± 2.3 m for T4. Combining the horizontal shortening value with the abandoned age of T4, the horizontal shortening rate was estimated to be 1.6 ± 0.2 mm/a since ~52 ka (Table 1, uncertainties were calculated by using Monte Carlo simulation) [127]. The terrace abandonment age is a maximum (described in Section 4.1.1), which results in the shortening rate becoming a maximum.

4.2. ADS Results

The upper horizons of the fold suffer from erosion, which in turn affects the integrity of these horizons. Therefore, in each seismic profile, four horizons, numbered 1–4, from the bottom up, which can be traced across the fold, were used to perform ADS analysis in this study (Figure 9, Figure 10 and Figure 11). The regional level of each horizon was defined as the northward extrapolation of the undeformed horizon imaged on the southern part of the seismic lines, considering its dip as the regional level dip. Where the horizon intersects the move line and the pin line mark, the regional level depths for each horizon are set. The reference level was defined as the interpreted, predicted detachment level. Following this strategy, we measured the excess areas (S) for the four horizons and the distances (h) of the regional levels against the reference levels at the pin line as well as the moving line in each seismic profile (Figure 9, Figure 10 and Figure 11). Figure 9c and Figure 11c show the area-depth diagrams. Area-depth relationships on both move lines and pin lines all show prominent linear trends with R2 values of >0.98 in all three seismic profiles, demonstrating that interpretations for all seismic reflection profiles are regarded as reasonably balanced and that the four horizons used to perform ADS analysis are pre-growth strata. Intercepts of the linear fits on the x-axis give the detachment adjustive values. For Line 1, the area-depth diagram gives the vertical distance between the intersections on the regional limits to be ~0.4 km. The horizontal distance between the regional limits is ~12 km. Thus, the angle between the reference level and the best-fit detachment level is calculated at 2°, and the maximum central angle is adjusted to 52° (Figure 9). For Line 2, the horizontal distance between the regional limits is ~14 km, and the vertical distance between the intersections on the regional limits is ~1.2 km. The angle between the reference level and the best-fit detachment level is calculated at 5°. The maximum central angle is adjusted to 50° (Figure 10). For Line 3, the vertical distance between the intersections on the regional limits is ~1.3 km, and the horizontal distance between the regional limits is ~15 km. The angle between the reference level and the best-fit detachment level is calculated at 5°. The maximum central angle is adjusted to 49° (Figure 11). According to Wang et al. [111], the crustal shortening determined from the area-depth diagram at the moving line has a very small error (~3%). Thus, in this study, we followed this strategy. The horizontal crustal shortenings for Line 1, Line 2, and Line 3 are 4.50, 4.44, and 3.10 km, respectively. The profile along the strike of the Huoyanshan fold shows that the displacement along the fold crest displays an asymmetric bow shape, and maximal displacement exists in locations around Line 1, with displacement decreasing eastward and westward (Figure 2c). In addition, from observation of the geological map of the Huoyanshan anticline (Figure 2b), it can be seen that the oldest strata are exposed at the locations of these three seismic lines, indicating that these locations have accumulated maximal shortening. Together, these two aspects demonstrate that the shortening of Line 1 can represent the overall shortening of the Huoayanshan anticline.

4.3. Best-Fit Model Parameters

Building on the information given by geometry interpretation and ADS analysis, the interpreted pre-growth horizons 1–4 and listric thrust geometries in Figure 9, Figure 10 and Figure 11 were read back and digitized as a deformed section to invert for the combinations of best-fit model parameters. Six parameters, the CC, RC, θmax, P/S ratio, trishear angle, and fault slip, were searched in a constrained parameter space in each seismic line. The CC, RC, and θmax were well interpreted in the seismic line and the fault slip was constrained by the ADS analysis. Therefore, the main objective of the inversion is to search over reasonable ranges of P/S and trishear angles. The search space of these parameters has lower limits [CCx−10, CCy−10, RC−10, 1.0 (P/S), 40 (trishear angle), 1 (fault slip)] and upper limits [CCx + 10, CCy + 10, RC + 10, 3.0, 80, 3], noting that CCx, CCy, RC, and fault slip are digitized values. The inversion was run horizon-by-horizon, for one horizon, the search was run for 10,000 iterations, and reannealing was applied every 100 iterations. The solutions of the best CC, RC, and θmax are very robust, and the parameters found are the same as, or only slightly different from, those of the interpreted values. The distributions of the fault slip, trishear angle, and P/S for the three seismic lines are shown in Figure 12, Figure 13 and Figure 14. For Line 1, the P/S values of the four horizons are between 1.9 and 2.4. The trishear angles for horizon 1, horizon 2, and horizon 4 are ~40°, and the trishear angle is ~50° in horizon 3. The ranges of fault slip values are 4.3–4.7 km on horizons 1–4 (Figure 12). For Line 2, the P/S is less spread, with values of 2.6–2.8 in all four horizons. The trishear angle, with values ranging from 40° to 45° in horizon 1, horizon 2, and horizon 3, and in horizon 2, the value is ~50°. Fault slip values can be divided into two sets. The lower horizons (horizons 1 and 2) have relatively smaller fault slips, with values between 3.2–3.7 km, and the upper horizons (horizons 3 and 4) have relatively larger fault slips, with values of approximately 4.2–4.3 km (Figure 13). For Line 3, the P/S values are ~3.0 on horizons 1, 3, and 4, and the P/S for horizon 2 is ~2.6. The trishear angle is less spread in horizon 1, horizon 2, and horizon 4, with values ranging from 40° to 44°, and in horizon 3, it is slightly higher, with values of ~48°. Fault slip has a very small spread, with values of approximately 3.0–3.1 km in all four horizons (Figure 14). As noted above, small spreads exist in the fault slip, trishear angle, and P/S from horizon to horizon, but in general, the values are centralized in each horizon. This is because there are uncertainties in the folded area where steeply dipping horizons are difficult to image and thus affect the solutions of the trishear angle and P/S. In addition, in comparison to the lower horizons, the upper horizons have larger extents in the backlimb which may have better constraints on the kinematics in the backlimb. Taking the latter aspect into account, we take the parameters inverted from horizon 4 as the best-fit model parameters in each seismic line. The inversion of horizon 4 in Line 1 results in a fault slip of 4.41 km, trishear angle of 40°, and P/S of 2.42. For Line 2, the best-fit model parameters are fault slip as 4.31 km, trishear angle as 40°, and P/S as 2.69. For Line 3, the best-fit model parameters are fault slip as 3.12 km, trishear angle as 40°, and P/S as 2.99. The displacement identified by the inversion matches the shortening estimated from ADS analysis in each seismic line. These consistent results confirm the reliability of the adopted inverse procedure.

5. Discussion

5.1. Estimation of Uncertainties in Structural Interpretation

Although high-quality seismic profiles provide an unprecedented level of detail in thrust settings, a good structural interpretation involves not only a well-imaged seismic profile but also a reliable depth conversion and an applicable theoretical model [128,129,130]. In our case, based on the features present in both surface and subsurface tectonic structures, we adopted a listric thrust-fold theory in our structural interpretation. The model we chose better accounts for features present in nature and avoids oversimplified kinematic interpretation. Then, we combined ADS analysis with reverse modeling of the interpreted structure to obtain a maximum amount of structural information. These structural analysis techniques quantified displacement and detachment depth, validated the original seismic interpretation, and furthermore gave the best-fit kinematic model parameters. However, there are still some uncertainties in each independent technique. For instance, despite the ADS analysis showing a good linear correlation in area-depth values for pre-growth horizons, which reveals a balanced interpretation, the detachment depth and overall shortening derived from ADS analysis may vary with the choice of seismic interpretation because seismic reflection profiles may not image well at the folded area where steeply dipping horizons exist. Similarly, with the inversion procedure, the question judging procedure works well or does not arise from the factor of the uniqueness of the optimal parameter solution. Evaluating the consistency of results obtained from ADS analysis and reverse modeling can help to answer this question and meanwhile estimate structural uncertainties present in interpretations. The inverse model-predicted displacement for each seismic line matches the ADS estimated displacement within a 3% deviation, which is likely within the uncertainty associated with seismic interpretation. In addition, when implementing the inverse modeling procedure, trishear is a continuous deformation. However, the computer models it as a series of discrete displacement steps, and this discrete displacement is very sensitive to inversion [119]. In this paper, we have used a sufficiently small slip increment (less than 0.8% of the total slip) to ensure that the inversion works well. In sum, this combined structure analysis approach provides greater insight into structural information than can be achieved by any independent technique and realizes mutual verification. This process has implications for guiding the development of fully quantified and robust constrained interpretations of complex structures.

5.2. Onset Time of the Huoyanshan Anticline and Its Reliability

For the Huoyanshan anticline, a total shortening of ~4.5 km was constrained in Section 4.2 and Section 4.3. When combined with the shortening rate of 1.6 ± 0.2 mm/a estimated by geomorphology and OSL dating in Section 4.1, it was inferred that the initiation of folding started at 1.8–3.7 Ma. This onset time of folding was determined from two aspects: the total shortening amount and the shortening rate. The total shortening value is quite robust because the combined structural analysis approach has quantified similar results, and the extent to which the results agree raises confidence in the reliability of this value. Whereas the use of the shortening rate is attached to a range of uncertainties, one must question the extent this rate can be linearly extrapolated in time. Although the assumption that geodetic slip rates are representative of Quaternary slip rates is challenged, there is still a good agreement between GPS and Quaternary slip rates [131]. The consistency between the long- and short-term horizontal shortening rates is broadly shown by the Kyrgyz Tian Shan [38], Pamir Frontal Thrust [53], and Longmen Shan piedmont [132,133]. Therefore, we begin with the comparison of GPS and geologic rates (here referred to as the N-S direction shortening rate) for faults in our study longitude. Results of present-day GPS measurements show that the N-S shortening rate across the eastern Chinese Tian Shan is 3–5 mm/a [76,77]. In the northern Bogda range, an ~1 mm/a shortening rate since the late Quaternary was estimated on the northern Bogda fault [58]. To the south of our study range, Wang et al. [94] constrained the shortening rate of the Kumysh fault to ~0.3 mm/a during the late Quaternary. Further south, in the Yanqi Basin, Huang et al. [93] estimated that the shortening rate of the north-edge fault is ~0.4–0.5 mm/a, and Li et al. [134] suggested that the shortening rate in the Hejing thrust-fold belt has been ~0.3 mm/a since the late Pleistocene. These available shortening rates of individual structures, taken together, indicate that the long-term rate of shortening across the eastern Chinese Tian Shan (except for the Houyanshan anticline) is ~2 mm/a. If we use the total GPS rates minus this long-term rate, the left rate is ~1–3 mm/a, which is a rough estimate. However, it is still consistent with the estimated shortening rate of the Huoyanshan anticline. Relying on this comparison, we conclude that the assumption of constant deformation rates since folding is acceptable. The onset time of the Huoyanshan anticline will be better constrained if we have the chronological constraint of the growth strata and/or long-term shortening rate. However, while these data are absent, our approach to the estimation of the onset time could still provide insight into the deformation process of the eastern Chinese Tian Shan range.

5.3. The 2 Ma Synchronous Growth Stage around the Tibetan Plateau

In addition, there seems to be a synchronous growth stage at ~2 Ma at the margins around the Tibetan Plateau. For instance, the Dushanzi anticline, located in the third structural belt of the northern Tian Shan piedmont, and an angular unconformity exists between the early Pleistocene and the middle Pleistocene strata [50]. Previous studies all suggested that the folding started at 2.6–1.5 Ma [39,50,52]. Likewise, according to analyses of the growth strata in the Longmen Shan piedmont, middle to late Pleistocene activity at 2.5–1.8 Ma was revealed [133,135,136,137]. At the deformed margin of the Himalayan foreland basin of northern India and Pakistan, the application of magnetic stratigraphy to constrain the chronology of sediment accumulation and unconformities shows that the thrust sheet began to form from 2.1 to 1.6 Ma [138,139]. At the Himalayan front of northwest India, the surface geology, interpretation of seismic reflection profiles, and magnetostratigraphic data suggest that a minimum shortening of 23 km has occurred since 1.9–1.5 Ma. At the Yarlung Tsangpo River, the eastern end of the Himalayan region, the chronology study of drill core sediments shows that the steepening of the Tsangpo Gorge started at approximately 2.5 to 2 Ma [140]. At the northern Tibetan Plateau, detailed magnetostratigraphy in the western Kunlun Mountains reflects a stage of uplift that has occurred since 1.6 Ma [141]. In northeastern Tibet, the uplifts of the Heli Shan and Jintannan Shan were reported at ~2 Ma [142,143,144]. As noted above, these studies indicate the ~2 Ma synchronous growth stage around the Tibetan Plateau and further confirm the reliability of our estimation of the deformation onset time of the Huoyanshan anticline.

6. Conclusions

The key findings of this paper are as follows:
(1)
Based on the features present in both surface and subsurface deformation, we first describe a curved thrust ramp and related fold pairs of the structure.
(2)
By using the kinematic model, deformed record, and OSL dating results of terrace T4, we estimate a geologic shortening rate of 1.6 ± 0.2 mm/a since ~52 ka.
(3)
We interpret three seismic reflection profiles using the theory of quantitative fault-related fold, ADS, and reverse modeling analyses. These profiles provide direct evidence that this structure connects by a listric thrust ramp to a shallow detachment level. ADS analysis reveals that the maximum shortening of the structure is ~4.5 km, which is consistent with a quantitative inverse model.
(4)
The above shortening rate and total shortening amount suggest that the structure may have formed at 1.8–3.7 Ma, a nearly synchronous growth stage around the Tibetan Plateau.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/rs15061603/s1.

Author Contributions

Conceptualization, Z.L. and X.Y.; methodology, Z.L., X.Y., W.W., and P.Z.; software, Z.L. and X.Y.; validation, Z.L., X.Y., W.W., and P.Z.; formal analysis, Z.L. and X.Y.; investigation, Z.L., X.Y., W.W., C.W., G.C., L.D., and K.L.; writing—original draft preparation, X.Y. and Z.L.; writing—review and editing, X.Y.; Z.L. and X.W.; supervision, P.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

This research is partially funded by the National Natural Science Foundation of China (42172233 and 42020301) and Guangdong Province Introduced Innovative R&D Team of Geological Processes and Natural Disasters around the South China Sea (2016ZT06N331).

Data Availability Statement

The seismic profiles presented in this study can be found in Supplementary Materials.

Acknowledgments

We thank Néstor Cardozo for providing trishear inverse modeling scripts, which are available on his web page: http://homepage.mac.com/nfcd/work/matlabScripts.html.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Armstrong, F.C.; Oriel, S.S. Tectonic development of Idaho-Wyoming thrust belt. AAPG Bull. 1965, 49, 1847–1866. [Google Scholar]
  2. Bally, A.W.; Gordy, P.; Stewart, G.A. Structure, seismic data, and orogenic evolution of southern Canadian Rocky Mountains. Bull. Can. Pet. Geol. 1966, 14, 337–381. [Google Scholar]
  3. Dahlstrom, C.D. Structural geology in the eastern margin of the Canadian Rocky Mountains. Bull. Can. Pet. Geol. 1970, 18, 332–406. [Google Scholar]
  4. Boyer, S.E. Styles of folding within thrust sheets: Examples from the Appalachian and Rocky Mountains of the USA and Canada. J. Struct. Geol. 1986, 8, 325–339. [Google Scholar] [CrossRef]
  5. Vann, I.; Graham, R.; Hayward, A. The structure of mountain fronts. J. Struct. Geol. 1986, 8, 215–227. [Google Scholar] [CrossRef]
  6. Erslev, E.A.; Mayborn, K.R. Multiple geometries and modes of fault-propagation folding in the Canadian thrust belt. J. Struct. Geol. 1997, 19, 321–335. [Google Scholar] [CrossRef]
  7. Hardy, S.; Ford, M. Numerical modeling of trishear fault propagation folding. Tectonics 1997, 16, 841–854. [Google Scholar] [CrossRef]
  8. Homza, T.X.; Wallace, W.K. Detachment folds with fixed hinges and variable detachment depth, northeastern Brooks Range, Alaska. J. Struct. Geol. 1997, 19, 337–354. [Google Scholar] [CrossRef]
  9. Delcaillau, B.; Deffontaines, B.; Floissac, L.; Angelier, J.; Deramond, J.; Souquet, P.; Chu, H.T.; Lee, J. Morphotectonic evidence from lateral propagation of an active frontal fold; Pakuashan anticline, foothills of Taiwan. Geomorphology 1998, 24, 263–290. [Google Scholar] [CrossRef]
  10. Burbank, D.; McLean, J.K.; Bullen, M.; Abdrakhmatov, K.; Miller, M.M. Partitioning of intermontane basins by thrust-related folding, Tien Shan, Kyrgyzstan. Basin Res. 1999, 11, 75. [Google Scholar] [CrossRef] [Green Version]
  11. Tozer, R.; Butler, R.; Corrado, S. Comparing thin-and thick-skinned thrust tectonic models of the Central Apennines, Italy. EGU Stephan Mueller Spec. Publ. Ser. 2002, 1, 181–194. [Google Scholar] [CrossRef]
  12. Heermance, R.V.; Chen, J.; Burbank, D.W.; Miao, J. Temporal constraints and pulsed Late Cenozoic deformation during the structural disruption of the active Kashi foreland, northwest China. Tectonics 2008, 27, TC6012. [Google Scholar] [CrossRef] [Green Version]
  13. McClay, K. Introduction to Thrust Fault-Related Folding. In Thrust Fault-Related Folding; American Association of Petroleum Geologists: Tulsa, OK, USA, 2011. [Google Scholar]
  14. Calamita, F.; Pace, P.; Satolli, S. Coexistence of fault-propagation and fault-bend folding in curve-shaped foreland fold-and-thrust belts: Examples from the Northern Apennines (Italy). Terra Nova 2012, 24, 396–406. [Google Scholar] [CrossRef]
  15. King, G.; Vita-Finzi, C. Active folding in the Algerian earthquake of 10 October 1980. Nature 1981, 292, 22–26. [Google Scholar] [CrossRef]
  16. Philip, H.; Meghraoui, M. Structural analysis and interpretation of the surface deformations of the El Asnam earthquake of October 10, 1980. Tectonics 1983, 2, 17–49. [Google Scholar] [CrossRef]
  17. Stein, R.S.; King, G.C. Seismic potential revealed by surface folding: 1983 Coalinga, California, earthquake. Science 1984, 224, 869–872. [Google Scholar] [CrossRef]
  18. Nicol, A.; Campbell, J.K. The impact of episodic fault-related folding on late Holocene degradation terraces along Waipara River, New Zealand. N. Z. J. Geol. Geophys. 2001, 44, 145–156. [Google Scholar] [CrossRef]
  19. Kao, H.; Chen, W.-P. The Chi-Chi earthquake sequence: Active, out-of-sequence thrust faulting in Taiwan. Science 2000, 288, 2346–2349. [Google Scholar] [CrossRef] [Green Version]
  20. Chen, Y.-G.; Chen, W.-S.; Lee, J.-C.; Lee, Y.-H.; Lee, C.-T.; Chang, H.-C.; Lo, C.-H. Surface rupture of 1999 Chi-Chi earthquake yields insights on active tectonics of central Taiwan. Bull. Seismol. Soc. Am. 2001, 91, 977–985. [Google Scholar] [CrossRef]
  21. Johnson, K.; Segall, P. Viscoelastic earthquake cycle models with deep stress-driven creep along the San Andreas fault system. J. Geophys. Res. Solid Earth 2004, 109, B10403. [Google Scholar] [CrossRef] [Green Version]
  22. Guzofski, C.A.; Shaw, J.H.; Lin, G.; Shearer, P.M. Seismically active wedge structure beneath the Coalinga anticline, San Joaquin basin, California. J. Geophys. Res. Solid Earth 2007, 112, B03S05. [Google Scholar] [CrossRef] [Green Version]
  23. Lin, A. Fossil Earthquakes: The Formation and Preservation of Pseudotachylytes; Springer: Berlin/Heidelberg, Germany, 2007; Volume 111. [Google Scholar]
  24. Allmendinger, R.W.; Shaw, J.H. Estimation of fault propagation distance from fold shape: Implications for earthquake hazard assessment. Geology 2000, 28, 1099–1102. [Google Scholar] [CrossRef]
  25. DeVecchio, D.E.; Keller, E.A.; Fuchs, M.; Owen, L.A. Late Pleistocene structural evolution of the Camarillo fold belt: Implications for lateral fault growth and seismic hazard in Southern California. Lithosphere 2012, 4, 91–109. [Google Scholar] [CrossRef] [Green Version]
  26. Saint-Carlier, D.; Charreau, J.; Lavé, J.; Blard, P.-H.; Dominguez, S.; Avouac, J.-P.; Wang, S.; Team, A. Major temporal variations in shortening rate absorbed along a large active fold of the southeastern Tianshan piedmont (China). Earth Planet. Sci. Lett. 2016, 434, 333–348. [Google Scholar] [CrossRef]
  27. Yang, H.; Yang, X.; Zhang, H.; Huang, X.; Huang, W.; Zhang, N. Active fold deformation and crustal shortening rates of the Qilian Shan Foreland Thrust Belt, NE Tibet, since the Late Pleistocene. Tectonophysics 2018, 742, 84–100. [Google Scholar] [CrossRef]
  28. Charreau, J.; Sartégou, A.; Saint-Carlier, D.; Lavé, J.; Blard, P.H.; Dominguez, S.; Wang, S.L.; Rao, G.; Team, A.; Aumaître, G. Late Miocene to Quaternary slip history across the Qiulitag anticline in the southern Tianshan piedmont. Terra Nova 2020, 32, 89–96. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, P.Z.; Deng, Q.D.; Yang, X.P.; Peng, S.Z.; Xu, X.W. Late Cenozoic tectonic deformation and mechanism along the Tianshan Mountain, northwestern China. Earthq. Res. China 1996, 12, 127–140, (In Chinese with English abstract). [Google Scholar]
  30. Abdrakhmatov, K.Y.; Aldazhanov, S.; Hager, B.; Hamburger, M.; Herring, T.; Kalabaev, K.; Makarov, V.; Molnar, P.; Panasyuk, S.; Prilepin, M. Relatively recent construction of the Tien Shan inferred from GPS measurements of present-day crustal deformation rates. Nature 1996, 384, 450–453. [Google Scholar] [CrossRef]
  31. Reigber, C.; Michel, G.; Galas, R.; Angermann, D.; Klotz, J.; Chen, J.; Papschev, A.; Arslanov, R.; Tzurkov, V.; Ishanov, M. New space geodetic constraints on the distribution of deformation in Central Asia. Earth Planet. Sci. Lett. 2001, 191, 157–165. [Google Scholar] [CrossRef]
  32. Yang, S.; Li, J.; Wang, Q. The deformation pattern and fault rate in the Tianshan Mountains inferred from GPS observations. Sci. China Ser. D Earth Sci. 2008, 51, 1064–1080. [Google Scholar] [CrossRef]
  33. Zubovich, A.V.; Wang, X.Q.; Scherba, Y.G.; Schelochkov, G.G.; Reilinger, R.; Reigber, C.; Mosienko, O.I.; Molnar, P.; Michajljow, W.; Makarov, V.I. GPS velocity field for the Tien Shan and surrounding regions. Tectonics 2010, 29, TC6014. [Google Scholar] [CrossRef]
  34. Zhang, P.-Z.; Shen, Z.; Wang, M.; Gan, W.; Burgmann, R.; Molnar, P.; Wang, Q.; Niu, Z.; Sun, J.; Wu, J. Continuous deformation of the Tibetan Plateau from global positioning system data. Geology 2004, 32, 809–812. [Google Scholar]
  35. Avouac, J.-P.; Tapponnier, P.; Bai, M.; You, H.; Wang, G. Active thrusting and folding along the northern Tien Shan and late Cenozoic rotation of the Tarim relative to Dzungaria and Kazakhstan. J. Geophys. Res. Solid Earth 1993, 98, 6755–6804. [Google Scholar] [CrossRef] [Green Version]
  36. Burchfiel, B.; Brown, E.; Qidong, D.; Xianyue, F.; Jun, L.; Molnar, P.; Jianbang, S.; Zhangming, W.; Huichuan, Y. Crustal shortening on the margins of the Tien Shan, Xinjiang, China. Int. Geol. Rev. 1999, 41, 665–700. [Google Scholar] [CrossRef]
  37. Deng, Q.D.; Feng, X.Y.; Zhang, P.Z.; Xu, X.W.; Yang, X.P.; Peng, S.Z.; Li, J. Active Tectonics of the Tian Shan Mountains; Seismology Press: Beijing, China, 2000; p. 399. (In Chinese) [Google Scholar]
  38. Thompson, S.C.; Weldon, R.J.; Rubin, C.M.; Abdrakhmatov, K.; Molnar, P.; Berger, G.W. Late Quaternary slip rates across the central Tien Shan, Kyrgyzstan, central Asia. J. Geophys. Res. Solid Earth 2002, 107, ETG 7-1–ETG 7-32. [Google Scholar] [CrossRef]
  39. Lu, H.; Li, B.; Wu, D.; Zhao, J.; Zheng, X.; Xiong, J.; Li, Y. Spatiotemporal patterns of the Late Quaternary deformation across the northern Chinese Tian Shan foreland. Earth-Sci. Rev. 2019, 194, 19–37. [Google Scholar] [CrossRef]
  40. Brown, E.T.; Bourlès, D.L.; Raisbeck, G.M.; Yiou, F.; Clark Burchfiel, B.; Molnar, P.; Qidong, D.; Jun, L. Estimation of slip rates in the southern Tien Shan using cosmic ray exposure dates of abandoned alluvial fans. Geol. Soc. Am. Bull. 1998, 110, 377–386. [Google Scholar] [CrossRef]
  41. Yin, A.; Nie, S.; Craig, P.; Harrison, T.; Ryerson, F.; Xianglin, Q.; Geng, Y. Late Cenozoic tectonic evolution of the southern Chinese Tian Shan. Tectonics 1998, 17, 1–27. [Google Scholar] [CrossRef]
  42. Hendrix, M.S.; Dumitru, T.A.; Graham, S.A. Late Oligocene-early Miocene unroofing in the Chinese Tian Shan: An early effect of the India-Asia collision. Geology 1994, 22, 487–490. [Google Scholar] [CrossRef]
  43. Molnar, P.; Brown, E.T.; Burchfiel, B.C.; Deng, Q.; Feng, X.; Li, J.; Raisbeck, G.M.; Shi, J.; Zhangming, W.; Yiou, F. Quaternary climate change and the formation of river terraces across growing anticlines on the north flank of the Tien Shan, China. J. Geol. 1994, 102, 583–602. [Google Scholar] [CrossRef]
  44. Bullen, M.; Burbank, D.; Garver, J.; Abdrakhmatov, K.Y. Late Cenozoic tectonic evolution of the northwestern Tien Shan: New age estimates for the initiation of mountain building. Geol. Soc. Am. Bull. 2001, 113, 1544–1559. [Google Scholar] [CrossRef]
  45. Lu, H.; Burbank, D.W.; Li, Y. Alluvial sequence in the north piedmont of the Chinese Tian Shan over the past 550 kyr and its relationship to climate change. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2010, 285, 343–353. [Google Scholar] [CrossRef]
  46. Thompson, J.A.; Burbank, D.W.; Li, T.; Chen, J.; Bookhagen, B. Late Miocene northward propagation of the northeast Pamir thrust system, northwest China. Tectonics 2015, 34, 510–534. [Google Scholar] [CrossRef] [Green Version]
  47. Thompson Jobe, J.A.; Li, T.; Chen, J.; Burbank, D.W.; Bufe, A. Quaternary tectonic evolution of the Pamir-Tian Shan convergence zone, Northwest China. Tectonics 2017, 36, 2748–2776. [Google Scholar] [CrossRef]
  48. Chen, J.; Burbank, D.; Scharer, K.; Sobel, E.; Yin, J.; Rubin, C.; Zhao, R. Magnetochronology of the Upper Cenozoic strata in the Southwestern Chinese Tian Shan: Rates of Pleistocene folding and thrusting. Earth Planet. Sci. Lett. 2002, 195, 113–130. [Google Scholar] [CrossRef]
  49. Chen, J.; Heermance, R.; Burbank, D.W.; Scharer, K.M.; Miao, J.; Wang, C. Quantification of growth and lateral propagation of the Kashi anticline, southwest Chinese Tian Shan. J. Geophys. Res. Solid Earth 2007, 112, B03S16. [Google Scholar] [CrossRef] [Green Version]
  50. Sun, J.; Zhu, R.; Bowler, J. Timing of the Tianshan Mountains uplift constrained by magnetostratigraphic analysis of molasse deposits. Earth Planet. Sci. Lett. 2004, 219, 239–253. [Google Scholar] [CrossRef]
  51. Charreau, J.; Avouac, J.-P.; Chen, Y.; Dominguez, S.; Gilder, S. Miocene to present kinematics of fault-bend folding across the Huerguosi anticline, northern Tianshan (China), derived from structural, seismic, and magnetostratigraphic data. Geology 2008, 36, 871–874. [Google Scholar] [CrossRef]
  52. Charreau, J.; Saint-Carlier, D.; Lavé, J.; Dominguez, S.; Blard, P.-H.; Avouac, J.-P.; Brown, N.D.; Malatesta, L.C.; Wang, S.; Rhodes, E.J. Late Pleistocene acceleration of deformation across the northern Tianshan piedmont (China) evidenced from the morpho-tectonic evolution of the Dushanzi anticline. Tectonophysics 2018, 730, 132–140. [Google Scholar] [CrossRef] [Green Version]
  53. Li, T.; Chen, J.; Thompson, J.A.; Burbank, D.W.; Xiao, W. Equivalency of geologic and geodetic rates in contractional orogens: New insights from the Pamir Frontal Thrust. Geophys. Res. Lett. 2012, 39, L15305. [Google Scholar] [CrossRef] [Green Version]
  54. Lu, H.; Zhang, T.; Zhao, J.; Si, S.; Wang, H.; Chen, S.; Zheng, X.; Li, Y. Late Quaternary alluvial sequence and uplift-driven incision of the Urumqi River in the north front of the Tian Shan, northwestern China. Geomorphology 2014, 219, 141–151. [Google Scholar] [CrossRef]
  55. Fu, X.; Li, S.-H.; Li, B.; Fu, B. A fluvial terrace record of late Quaternary folding rate of the Anjihai anticline in the northern piedmont of Tian Shan, China. Geomorphology 2017, 278, 91–104. [Google Scholar] [CrossRef]
  56. Peng, S. Active Tectonics and Earthquake Hazards of the Turfan Basin, Northwestern China. Sci. China Ser. D Earth Sci. 1995, 46, 13–24, (In Chinese with English abstract). [Google Scholar]
  57. Ren, G.; Li, C.; Wu, C.; Wang, S.; Zhang, H.; Ren, Z.; Li, X. The late Quaternary activity and formation mechanism of Baoertu Fault zone, eastern Tian Shan segment. Seismol. Geol. 2019, 41, 856–871, (In Chinese with English abstract). [Google Scholar]
  58. Wu, C.; Wu, G.; Shen, J.; Dai, X.; Chen, J.; Song, H. Late Quaternary tectonic activity and crustal shortening rate of the Bogda mountain area, eastern Tian Shan, China. J. Asian Earth Sci. 2016, 119, 20–29. [Google Scholar] [CrossRef] [Green Version]
  59. Deng, Q.; Zhang, P.; Xu, X.; Yang, X.; Peng, S.; Feng, X. Paleoseismology of the northern piedmont of Tianshan Mountains, northwestern China. J. Geophys. Res. Solid Earth 1996, 101, 5895–5920. [Google Scholar] [CrossRef]
  60. Fu, B.; Lin, A.; Kano, K.-I.; Maruyama, T.; Guo, J. Quaternary folding of the eastern Tian Shan, northwest China. Tectonophysics 2003, 369, 79–101. [Google Scholar] [CrossRef]
  61. Cardozo, N.; Jackson, C.A.-L.; Whipp, P.S. Determining the uniqueness of best-fit trishear models. J. Struct. Geol. 2011, 33, 1063–1078. [Google Scholar]
  62. Cardozo, N.; Brandenburg, J. Kinematic modeling of folding above listric propagating thrusts. J. Struct. Geol. 2014, 60, 1–12. [Google Scholar] [CrossRef]
  63. Allen, M.; Windley, B.; Zhang, C. Palaeozoic collisional tectonics and magmatism of the Chinese Tien Shan, central Asia. Tectonophysics 1993, 220, 89–115. [Google Scholar] [CrossRef]
  64. Burtman, V. Structural geology of variscan Tien Shan, USSR. Am. J. Sci. 1975, 275, 157–186. [Google Scholar]
  65. Windley, B.; Allen, M.; Zhang, C.; Zhao, Z.; Wang, G. Paleozoic accretion and Cenozoic redeformation of the Chinese Tien Shan range, central Asia. Geology 1990, 18, 128–131. [Google Scholar] [CrossRef]
  66. Hendrix, M.S.; Graham, S.A.; Carroll, A.R.; Sobel, E.R.; McKNIGHT, C.L.; Schulein, B.J.; Wang, Z. Sedimentary record and climatic implications of recurrent deformation in the Tian Shan: Evidence from Mesozoic strata of the north Tarim, south Junggar, and Turpan basins, northwest China. Geol. Soc. Am. Bull. 1992, 104, 53–79. [Google Scholar] [CrossRef]
  67. Carroll, A.; Graham, S.; Hendrix, M.; Ying, D.; Zhou, D. Late Paleozoic tectonic amalgamation of northwestern China: Sedimentary record of the northern Tarim, northwestern Turpan, and southern Junggar basins. Geol. Soc. Am. Bull. 1995, 107, 571–594. [Google Scholar] [CrossRef]
  68. Dumitru, T.A.; Zhou, D.; Chang, E.Z.; Graham, S.A.; Hendrix, M.S.; Sobel, E.R.; Carroll, A.R. Uplift, exhumation, and deformation in the Chinese Tian Shan. Mem.-Geol. Soc. Am. 2001, 194, 71–100. [Google Scholar]
  69. Tapponnier, P.; Molnar, P. Active faulting and tectonics in China. J. Geophys. Res. 1977, 82, 2905–2930. [Google Scholar] [CrossRef]
  70. Gao, J.; Long, L.; Klemd, R.; Qian, Q.; Liu, D.; Xiong, X.; Su, W.; Liu, W.; Wang, Y.; Yang, F. Tectonic evolution of the South Tianshan orogen and adjacent regions, NW China: Geochemical and age constraints of granitoid rocks. Int. J. Earth Sci. 2009, 98, 1221–1238. [Google Scholar] [CrossRef]
  71. Sobel, E.R.; Dumitru, T.A. Thrusting and exhumation around the margins of the western Tarim basin during the India-Asia collision. J. Geophys. Res. Solid Earth 1997, 102, 5043–5063. [Google Scholar] [CrossRef]
  72. Sobel, E.R.; Chen, J.; Heermance, R.V. Late Oligocene–Early Miocene initiation of shortening in the Southwestern Chinese Tian Shan: Implications for Neogene shortening rate variations. Earth Planet. Sci. Lett. 2006, 247, 70–81. [Google Scholar] [CrossRef]
  73. Allen, M.B.; Vincent, S.J. Fault reactivation in the Junggar region, northwest China: The role of basement structures during Mesozoic-Cenozoic compression. J. Geol. Soc. 1997, 154, 151–155. [Google Scholar] [CrossRef]
  74. Charvet, J.; Shu, L.S.; Laurent-Charvet, S. Paleozoic structural and geodynamic evolution of eastern Tianshan (NW China): Welding of the Tarim and Junggar plates. Epis. J. Int. Geosci. 2007, 30, 162–186. [Google Scholar]
  75. Buslov, M.; De Grave, J.; Bataleva, E.; Batalev, V.Y. Cenozoic tectonic and geodynamic evolution of the Kyrgyz Tien Shan Mountains: A review of geological, thermochronological and geophysical data. J. Asian Earth Sci. 2007, 29, 205–214. [Google Scholar] [CrossRef]
  76. 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, e2019JB018774. [Google Scholar] [CrossRef] [Green Version]
  77. Liang, S.; Gan, W.; Shen, C.; Xiao, G.; Liu, J.; Chen, W.; Ding, X.; Zhou, D. Three-dimensional velocity field of present-day crustal motion of the Tibetan Plateau derived from GPS measurements. J. Geophys. Res. Solid Earth 2013, 118, 5722–5732. [Google Scholar] [CrossRef]
  78. Métivier, F.; Gaudemer, Y. Mass transfer between eastern Tien Shan and adjacent basins (central Asia): Constraints on regional tectonics and topography. Geophys. J. Int. 1997, 128, 1–17. [Google Scholar] [CrossRef] [Green Version]
  79. Bullen, M.; Burbank, D.; Garver, J. Building the northern Tien Shan: Integrated thermal, structural, and topographic constraints. J. Geol. 2003, 111, 149–165. [Google Scholar] [CrossRef] [Green Version]
  80. Huang, B.; Piper, J.D.; Peng, S.; Liu, T.; Li, Z.; Wang, Q.; Zhu, R. Magnetostratigraphic study of the Kuche Depression, Tarim Basin, and Cenozoic uplift of the Tian Shan range, western China. Earth Planet. Sci. Lett. 2006, 251, 346–364. [Google Scholar] [CrossRef] [Green Version]
  81. Ji, J.; Luo, P.; White, P.; Jiang, H.; Gao, L.; Ding, Z. Episodic uplift of the Tianshan Mountains since the late Oligocene constrained by magnetostratigraphy of the Jingou River section, in the southern margin of the Junggar Basin, China. J. Geophys. Res. Solid Earth 2008, 113, B05102. [Google Scholar] [CrossRef]
  82. Charreau, J.; Chen, Y.; Gilder, S.; Dominguez, S.; Avouac, J.-P.; Sen, S.; Sun, D.; Li, Y.; Wang, W.-M. Magnetostratigraphy and rock magnetism of the Neogene Kuitun He section (northwest China): Implications for Late Cenozoic uplift of the Tianshan mountains. Earth Planet. Sci. Lett. 2005, 230, 177–192. [Google Scholar] [CrossRef] [Green Version]
  83. Charreau, J.; Gilder, S.; Chen, Y.; Dominguez, S.; Avouac, J.-P.; Sen, S.; Jolivet, M.; Li, Y.; Wang, W. Magnetostratigraphy of the Yaha section, Tarim Basin (China): 11 Ma acceleration in erosion and uplift of the Tian Shan mountains. Geology 2006, 34, 181–184. [Google Scholar] [CrossRef] [Green Version]
  84. Sun, J.; Zhang, Z. Syntectonic growth strata and implications for late Cenozoic tectonic uplift in the northern Tian Shan, China. Tectonophysics 2009, 463, 60–68. [Google Scholar] [CrossRef]
  85. Daëron, M.; Avouac, J.P.; Charreau, J. Modeling the shortening history of a fault tip fold using structural and geomorphic records of deformation. J. Geophys. Res. Solid Earth 2007, 112, 3–13. [Google Scholar] [CrossRef] [Green Version]
  86. Li, T.; Chen, J.; Thompson, J.A.; Burbank, D.W.; Yang, X. Quantification of three-dimensional folding using fluvial terraces: A case study from the Mushi anticline, northern margin of the Chinese Pamir. J. Geophys. Res. Solid Earth 2013, 118, 4628–4647. [Google Scholar] [CrossRef] [Green Version]
  87. Li, T.; Chen, J.; Thompson, J.A.; Burbank, D.W.; Yang, H. Hinge-migrated fold-scarp model based on an analysis of bed geometry: A study from the Mingyaole anticline, southern foreland of Chinese Tian Shan. J. Geophys. Res. Solid Earth 2015, 120, 6592–6613. [Google Scholar] [CrossRef] [Green Version]
  88. Lu, H.; Burbank, D.W.; Li, Y.; Liu, Y. Late Cenozoic structural and stratigraphic evolution of the northern Chinese Tian Shan foreland. Basin Res. 2010, 22, 249–269. [Google Scholar] [CrossRef]
  89. Lu, H.; Wang, Z.; Zhang, T.; Zhao, J.; Zheng, X.; Li, Y. Latest Miocene to Quaternary deformation in the southern Chaiwopu Basin, northern Chinese Tian Shan foreland. J. Geophys. Res. Solid Earth 2015, 120, 8656–8671. [Google Scholar] [CrossRef] [Green Version]
  90. Lü, L.; Sun, J.; Jia, Y.; Wu, L. Late Cenozoic thrust propagation within the Keping fold-and-thrust belt along the southern foreland of Chinese Tian Shan: Evidence from apatite (UTh)/He results. Tectonophysics 2021, 814, 228966. [Google Scholar] [CrossRef]
  91. Allen, M.; Windley, B.; Chi, Z.; Zhong-Yan, Z.; Guang-Rei, W. Basin evolution within and adjacent to the Tien Shan Range, NW China. J. Geol. Soc. 1991, 148, 369–378. [Google Scholar] [CrossRef]
  92. Nishidai, T.; Berry, J. Geological interpretation and hydrocarbon potential of the Turpan Basin (NW China) from satellite imagery. In Proceedings of the Thematic Conference on Geologic Remote Sensing, 8th, Denver, CO, USA, 29 April–2 May 1991. [Google Scholar]
  93. Huang, W.; Yang, X.; Li, S.; Yang, H. The late Quaternary activity characteristics of the strike-slip faults in the Tianshan orogenic belt: A case study of Kaiduhe Fault. Seismol. Geol. 2018, 40, 1040–1058. [Google Scholar]
  94. Wang, S.; Jiao, R.; Ren, Z.; Wu, C.; Ren, G.; Zhang, H.; Lei, Q. Active thrusting in an intermontane basin: The Kumysh Fault, eastern Tian Shan. Tectonics 2020, 39, e2019TC006029. [Google Scholar] [CrossRef]
  95. Yang, X.; Wu, C.; Li, Z.; Wang, W.; Chen, G.; Duan, L. Late Quaternary Kinematics and Deformation Rate of the Huoyanshan Structure Derived from Deformed River Terraces in the South Piedmont of the Eastern Chinese Tian Shan. Front. Earth Sci. 2021, 9, 649011. [Google Scholar] [CrossRef]
  96. Shao, L.; Stattegger, K.; Li, W.; Haupt, B.J. Depositional style and subsidence history of the Turpan Basin (NW China). Sediment. Geol. 1999, 128, 155–169. [Google Scholar] [CrossRef]
  97. Carroll, A.; Graham, S.; Hendrix, M.; Chu, J.; McKnight, C.; Feng, Y.; Liang, Y.; Xiao, X.; Zhao, M.; Tang, Y. Characteristics of sedimentation in Late Paleozoic Junggar basin and its basement. In On Tectonic Evolution of the Southern Margin of the Paleoasian Composite Megasuture Zone; Beijing Science and Technic Publishing House: Beijing, China, 1991; pp. 136–150. [Google Scholar]
  98. Amos, C.B.; Burbank, D.W.; Nobes, D.C.; Read, S.A. Geomorphic constraints on listric thrust faulting: Implications for active deformation in the Mackenzie Basin, South Island, New Zealand. J. Geophys. Res. Solid Earth 2007, 112, B03S11. [Google Scholar] [CrossRef] [Green Version]
  99. Trexler, C.C.; Cowgill, E.; Spencer, J.Q.; Godoladze, T. Rate of active shortening across the southern thrust front of the Greater Caucasus in western Georgia from kinematic modeling of folded river terraces above a listric thrust. Earth Planet. Sci. Lett. 2020, 544, 116362. [Google Scholar] [CrossRef]
  100. Aitken, M.J. Physics and Archaeology; CERN-VIDEO-C-204-B; CERN: Meyrin, Switzerland, 1995. [Google Scholar]
  101. Aitken, M.J. Introduction to Optical Dating: The Dating of Quaternary Sediments by the Use of Photon-Stimulated Luminescence; Clarendon Press: Oxford, UK, 1998. [Google Scholar]
  102. Prescott, J.R.; Hutton, J.T. Cosmic ray contributions to dose rates for luminescence and ESR dating: Large depths and long-term time variations. Radiat. Meas. 1994, 23, 497–500. [Google Scholar] [CrossRef]
  103. Chen, K.G.C.; Augier, R.; Chen, Y.; Mei, Y.H.; Lin, W.; Wang, Q.C. A multidiscipline method of geological survey, seismic line and gravity measurement applied to fold-and-thrust belt: A case study along the Hutubi River in the northern piedmont of Tianshan. Chin. J. Geophys. 2014, 57, 75–87. [Google Scholar]
  104. Gonzalez-Mieres, R.; Suppe, J. Relief and shortening in detachment folds. J. Struct. Geol. 2006, 28, 1785–1807. [Google Scholar] [CrossRef]
  105. Gonzalez-Mieres, R.; Suppe, J. Shortening Histories in Active Detachment Folds Based on Area-of-Relief Methods; American Association of Petroleum Geologists: Tulsa, OK, USA, 2011. [Google Scholar]
  106. Groshong, R.H., Jr.; Withjack, M.O.; Schlische, R.W.; Hidayah, T.N. Bed length does not remain constant during deformation: Recognition and why it matters. J. Struct. Geol. 2012, 41, 86–97. [Google Scholar] [CrossRef]
  107. Moretti, I.; Callot, J.P. Area, length and thickness conservation: Dogma or reality? J. Struct. Geol. 2012, 41, 64–75. [Google Scholar] [CrossRef]
  108. Schlische, R.W.; Groshong Jr, R.H.; Withjack, M.O.; Hidayah, T.N. Quantifying the geometry, displacements, and subresolution deformation in thrust-ramp anticlines with growth and erosion: From models to seismic-reflection profile. J. Struct. Geol. 2014, 69, 304–319. [Google Scholar] [CrossRef]
  109. Groshong, R.H., Jr. Quality control and risk assessment of seismic profiles using area-depth-strain analysis. Interpretation 2015, 3, SAA1–SAA15. [Google Scholar] [CrossRef]
  110. Eichelberger, N.W.; Hughes, A.N.; Nunns, A.G. Combining multiple quantitative structural analysis techniques to create robust structural interpretations. Interpretation 2015, 3, SAA89–SAA104. [Google Scholar] [CrossRef]
  111. Wang, W.; Yin, H.; Jia, D.; Wu, Z.; Wu, C.; Zhou, P. Calculating detachment depth and dip angle in sedimentary wedges using the area–depth graph. J. Struct. Geol. 2018, 107, 1–11. [Google Scholar] [CrossRef]
  112. Qiu, J.; Rao, G.; Wang, X.; Yang, D.; Xiao, L. Effects of fault slip distribution on the geometry and kinematics of the southern Junggar fold-and-thrust belt, northern Tian Shan. Tectonophysics 2019, 772, 228209. [Google Scholar] [CrossRef]
  113. Gao, L.; Rao, G.; Tang, P.; Qiu, J.; Peng, Z.; Pei, Y.; Yu, Y.; Zhao, B.; Wang, R. Structural development at the leading edge of the salt-bearing Kuqa fold-and-thrust belt, southern Tian Shan, NW China. J. Struct. Geol. 2020, 140, 104184. [Google Scholar] [CrossRef]
  114. Epard, J.-L.; Groshong, R.H., Jr. Excess area and depth to detachment. AAPG Bull. 1993, 77, 1291–1302. [Google Scholar]
  115. Groshong, R.H., Jr.; Epard, J.-L. The role of strain in area-constant detachment folding. J. Struct. Geol. 1994, 16, 613–618. [Google Scholar] [CrossRef]
  116. Erslev, E.A. Trishear fault-propagation folding. Geology 1991, 19, 617–620. [Google Scholar] [CrossRef]
  117. Zehnder, A.T.; Allmendinger, R.W. Velocity field for the trishear model. J. Struct. Geol. 2000, 22, 1009–1014. [Google Scholar] [CrossRef]
  118. Allmendinger, R.W. Inverse and forward numerical modeling of trishear fault-propagation folds. Tectonics 1998, 17, 640–656. [Google Scholar] [CrossRef]
  119. Allmendinger, R.W.; Zapata, T.; Manceda, R.; Dzelalija, F. Trishear Kinematic Modeling of Structures, with Examples from the Neuqun Basin, Argentina. AAPG Mem. 2004, 82, 356–371. [Google Scholar]
  120. MathWorks, T. Global Optimization Toolbox User’s Guide; The MathWorks Inc.: Natick, MA, USA, 2013. [Google Scholar]
  121. Suppe, J. Geometry and kinematics of fault-bend folding. Am. J. Sci. 1983, 283, 684–721. [Google Scholar]
  122. Suppe, J.; Medwedeff, D.A. Geometry and kinematics of fault-propagation folding. Eclogae Geol. Helv. 1990, 83, 409–454. [Google Scholar]
  123. Hardy, S.; Poblet, J. Geometric and numerical model of progressive limb rotation in detachment folds. Geology 1994, 22, 371–374. [Google Scholar]
  124. Suppe, J.; Connors, C.D.; Zhang, Y. Shear Fault-Bend Folding; American Association of Petroleum Geologists: Tulsa, OK, USA, 2004. [Google Scholar]
  125. Suppe, J.; Sàbat, F.; Munoz, J.A.; Poblet, J.; Roca, E.; Vergés, J. Bed-by-bed fold growth by kink-band migration: Sant Llorenç de Morunys, eastern Pyrenees. J. Struct. Geol. 1997, 19, 443–461. [Google Scholar]
  126. Seeber, L.; Sorlien, C.C. Listric thrusts in the western Transverse Ranges, California. Geol. Soc. Am. Bull. 2000, 112, 1067–1079. [Google Scholar]
  127. Mooney, C.Z. Monte Carlo Simulation; Sage: Thousand Oaks, CA, USA, 1997. [Google Scholar]
  128. Etris, E.L.; Crabtree, N.J.; Dewar, J.; Pickford, S. True depth conversion: More than a pretty picture. CSEG Rec. 2001, 26, 11–22. [Google Scholar]
  129. Torvela, T.; Bond, C.E. Do experts use idealised structural models? Insights from a deepwater fold–thrust belt. J. Struct. Geol. 2011, 33, 51–58. [Google Scholar] [CrossRef]
  130. Totake, Y.; Butler, R.W.; Bond, C.E. Structural validation as an input into seismic depth conversion to decrease assigned structural uncertainty. J. Struct. Geol. 2017, 95, 32–47. [Google Scholar] [CrossRef] [Green Version]
  131. Mohadjer, S.; Ehlers, T.A.; Bendick, R.; Mutz, S.G. Review of GPS and Quaternary fault slip rates in the Himalaya-Tibet orogen. Earth-Sci. Rev. 2017, 174, 39–52. [Google Scholar] [CrossRef]
  132. Willett, S.D.; Brandon, M.T. On steady states in mountain belts. Geology 2002, 30, 175–178. [Google Scholar] [CrossRef]
  133. Li, Z.; Liu-Zeng, J.; Jia, D.; Sun, C.; Wang, W.; Yuan, Z.; Liu, B. Quaternary activity of the range front thrust system in the Longmen Shan piedmont, China, revealed by seismic imaging and growth strata. Tectonics 2016, 35, 2807–2827. [Google Scholar] [CrossRef] [Green Version]
  134. Li, A.; Yang, X.; Huang, W.; Yiliyaer, Y. Quaternary deformation of the Hejing thrust-fold belt on northern margin of the Yanqi Basin, southern Tianshan. Seismol. Geol. 2012, 3, 240–253. [Google Scholar]
  135. Li, Y.; Zhou, R.-J.; Densmore, A.; Ellis, M.; Li, B. Sedimentary responses to Late Cenozoic thrusting and strike-slipping of Longmen Shan along eastern margin of Tibetan Plateau. Acta Sedimentol. Sin. 2006, 24, 153. [Google Scholar]
  136. Zhao, Z.; Qiao, Y.; Wang, Y.; Fu, J.; Wang, S.; Li, C.; Yao, H.; Jiang, F. Magnetostratigraphic and paleoclimatic studies on the red earth formation from the Chengdu Plain in Sichuan province, China. Sci. China Ser. D Earth Sci. 2007, 50, 927–935. [Google Scholar] [CrossRef]
  137. Kong, P.; Zheng, Y.; Fu, B. Cosmogenic nuclide burial ages and provenance of Late Cenozoic deposits in the Sichuan Basin: Implications for Early Quaternary glaciations in east Tibet. Quat. Geochronol. 2011, 6, 304–312. [Google Scholar] [CrossRef]
  138. Johnson, G.D.; Raynolds, R.G.; Burbank, D.W. Late Cenozoic Tectonics and Sedimentation in the North-Western Himalayan Foredeep: I. Thrust Ramping and Associated Deformation in the Potwar Region. In Foreland Basins; Blackwell Publishing: Oxford, UK, 1986; Volume 8, pp. 273–291. [Google Scholar]
  139. Baker, D.M.; Lillie, R.J.; Yeats, R.S.; Johnson, G.D.; Yousuf, M.; Zamin, A.S.H. Development of the Himalayan frontal thrust zone: Salt Range, Pakistan. Geology 1988, 16, 3–7. [Google Scholar] [CrossRef]
  140. Wang, P.; Scherler, D.; Liu-Zeng, J.; Mey, J.; Avouac, J.-P.; Zhang, Y.; Shi, D. Tectonic control of Yarlung Tsangpo Gorge revealed by a buried canyon in Southern Tibet. Science 2014, 346, 978–981. [Google Scholar] [CrossRef] [Green Version]
  141. Zheng, H.; Powell, C.M.; An, Z.; Zhou, J.; Dong, G. Pliocene uplift of the northern Tibetan Plateau. Geology 2000, 28, 715–718. [Google Scholar] [CrossRef]
  142. Tapponnier, P.; Meyer, B.; Avouac, J.P.; Peltzer, G.; Gaudemer, Y.; Guo, S.; Xiang, H.; Yin, K.; Chen, Z.; Cai, S. Active thrusting and folding in the Qilian Shan, and decoupling between upper crust and mantle in northeastern Tibet. Earth Planet. Sci. Lett. 1990, 97, 382–403. [Google Scholar] [CrossRef]
  143. Métivier, F.; Gaudemer, Y.; Tapponnier, P.; Meyer, B. Northeastward growth of the Tibet plateau deduced from balanced reconstruction of two depositional areas: The Qaidam and Hexi Corridor basins, China. Tectonics 1998, 17, 823–842. [Google Scholar] [CrossRef] [Green Version]
  144. Zheng, W.J.; Zhang, P.Z.; Ge, W.P.; Molnar, P.; Zhang, H.P.; Yuan, D.Y.; Liu, J.H. Late Quaternary slip rate of the South Heli Shan Fault (northern Hexi Corridor, NW China) and its implications for northeastward growth of the Tibetan Plateau. Tectonics 2013, 32, 271–293. [Google Scholar] [CrossRef]
Figure 1. (a) Topography of the Tian Shan, showing its location and geological setting. The blue box shows the location of (b). (b) Simplified tectonic map showing major active structures in the Tian Shan. The blue box shows the location of (c). (c) Simplified tectonic map of the eastern Chinese Tian Shan. Faults are modified by Deng et al. [37] and Wu et al. [58]. Red hollow circles indicate earthquake’s magnitude (Mw or mb) since 1900 (earthquake data are downloaded from https://earthquake.usgs.gov/earthquakes/search/ (accessed on 22 September 2022)).
Figure 1. (a) Topography of the Tian Shan, showing its location and geological setting. The blue box shows the location of (b). (b) Simplified tectonic map showing major active structures in the Tian Shan. The blue box shows the location of (c). (c) Simplified tectonic map of the eastern Chinese Tian Shan. Faults are modified by Deng et al. [37] and Wu et al. [58]. Red hollow circles indicate earthquake’s magnitude (Mw or mb) since 1900 (earthquake data are downloaded from https://earthquake.usgs.gov/earthquakes/search/ (accessed on 22 September 2022)).
Remotesensing 15 01603 g001
Figure 2. (a) Topography of the Huoyanshan structure (location in Figure 1c), lines 1–3 are seismic lines. (b) Simplified geologic map (modified from Peng) [56], location in (a). (c) Strike profile of the Huoyanshan anticline (location of B–B′ in (b)).
Figure 2. (a) Topography of the Huoyanshan structure (location in Figure 1c), lines 1–3 are seismic lines. (b) Simplified geologic map (modified from Peng) [56], location in (a). (c) Strike profile of the Huoyanshan anticline (location of B–B′ in (b)).
Remotesensing 15 01603 g002
Figure 3. Geological cross-section showing the structural style and the morphology pattern of stratigraphic units of the FTB in the Turpan Basin (modified from Nishidai and Berry [92], Fu et al. [60], location of A–A′ in Figure 1c), T: Triassic unis, J1: Lower Jurassic units, J2: Middle Jurassic units, J3: Upper Jurassic units, K: Cretaceous units, E: Paleogene units, N1: Miocene units, N2: Pliocene units, Q: Quaternary units, see text for details.
Figure 3. Geological cross-section showing the structural style and the morphology pattern of stratigraphic units of the FTB in the Turpan Basin (modified from Nishidai and Berry [92], Fu et al. [60], location of A–A′ in Figure 1c), T: Triassic unis, J1: Lower Jurassic units, J2: Middle Jurassic units, J3: Upper Jurassic units, K: Cretaceous units, E: Paleogene units, N1: Miocene units, N2: Pliocene units, Q: Quaternary units, see text for details.
Remotesensing 15 01603 g003
Figure 4. Schematic regional stratigraphic chart illustrating deposition sequences of the FTB in the Turpan Basin (modified from Shao et al. [96] and references therein).
Figure 4. Schematic regional stratigraphic chart illustrating deposition sequences of the FTB in the Turpan Basin (modified from Shao et al. [96] and references therein).
Remotesensing 15 01603 g004
Figure 5. (a) Orthography imaging of the Shengjingou Valley (location in Figure 2b). Blue flags represent the surveyed ground control points. Red stars represent the location of OSL sample. (b) High-resolution DEM surveyed by UVA. (c) DEM shaded topography showing the locations of terraces identified in the Shengjingou Valley. Colors and numbers represent different terrace generations, and black hollow dots show survey points on terraces and river channels.
Figure 5. (a) Orthography imaging of the Shengjingou Valley (location in Figure 2b). Blue flags represent the surveyed ground control points. Red stars represent the location of OSL sample. (b) High-resolution DEM surveyed by UVA. (c) DEM shaded topography showing the locations of terraces identified in the Shengjingou Valley. Colors and numbers represent different terrace generations, and black hollow dots show survey points on terraces and river channels.
Remotesensing 15 01603 g005
Figure 6. (a) Red arrows show the location of the fault scarp. (b,c) Views N and E of the river terraces across the Shengjingou Valley (T1–T6, from young to old). (df) Views S, N, and S of the bedrock outcrop in the Shengjingou Valley (locations in Figure 5a).
Figure 6. (a) Red arrows show the location of the fault scarp. (b,c) Views N and E of the river terraces across the Shengjingou Valley (T1–T6, from young to old). (df) Views S, N, and S of the bedrock outcrop in the Shengjingou Valley (locations in Figure 5a).
Remotesensing 15 01603 g006
Figure 7. (a) Longitudinal profiles of terraces T1–T8 and river channel across the Shengjingou Valley. (b) Topographic profiles of terraces T1–T8 that removed the modern river gradient. The distance is projected along the line C–C′ (location see Figure 5c), where C is the direction upstream. In both profiles, the vertical exaggeration is 10.
Figure 7. (a) Longitudinal profiles of terraces T1–T8 and river channel across the Shengjingou Valley. (b) Topographic profiles of terraces T1–T8 that removed the modern river gradient. The distance is projected along the line C–C′ (location see Figure 5c), where C is the direction upstream. In both profiles, the vertical exaggeration is 10.
Remotesensing 15 01603 g007
Figure 8. (a) Terrace deformation patterns and bedding attitudes were used to deduce the subsurface geometry and correlate the kinematic model, and the profile for T4 was used to quantify the uplift it recorded. (b) Kinematic model modified from Amos et al. [98], Cardozo and Brandenburg [62], and Trexler et al. [99].
Figure 8. (a) Terrace deformation patterns and bedding attitudes were used to deduce the subsurface geometry and correlate the kinematic model, and the profile for T4 was used to quantify the uplift it recorded. (b) Kinematic model modified from Amos et al. [98], Cardozo and Brandenburg [62], and Trexler et al. [99].
Remotesensing 15 01603 g008
Figure 9. Seismic interpretation and ADS analysis of Line 1 (location in Figure 2a). (a) Interpretations of the seismic reflection profile. Vertical exaggeration is 1:1. (b) Areas and heights measured for horizons 1–4. (c) Area-depth diagram illustrating the area-depth relationships on both move lines and pin lines with R2 values of >0.99.
Figure 9. Seismic interpretation and ADS analysis of Line 1 (location in Figure 2a). (a) Interpretations of the seismic reflection profile. Vertical exaggeration is 1:1. (b) Areas and heights measured for horizons 1–4. (c) Area-depth diagram illustrating the area-depth relationships on both move lines and pin lines with R2 values of >0.99.
Remotesensing 15 01603 g009
Figure 10. Seismic interpretation and ADS analysis of Line 2 (location in Figure 2a). (a) Interpretations of the seismic reflection profile. Vertical exaggeration is 1:1. (b) Areas and heights measured for horizons 1–4. (c) Area-depth diagram illustrating the area-depth relationships on both move lines and pin lines, with R2 values of >0.98.
Figure 10. Seismic interpretation and ADS analysis of Line 2 (location in Figure 2a). (a) Interpretations of the seismic reflection profile. Vertical exaggeration is 1:1. (b) Areas and heights measured for horizons 1–4. (c) Area-depth diagram illustrating the area-depth relationships on both move lines and pin lines, with R2 values of >0.98.
Remotesensing 15 01603 g010
Figure 11. Seismic interpretation and ADS analysis of Line 3 (location in Figure 2a). (a) Interpretations of seismic reflection profile. Vertical exaggeration is 1:1. (b) Areas and heights measured for horizons 1–4. (c) Area-depth diagram illustrating the area-depth relationships on both move lines and pin lines, with R2 values of >0.98.
Figure 11. Seismic interpretation and ADS analysis of Line 3 (location in Figure 2a). (a) Interpretations of seismic reflection profile. Vertical exaggeration is 1:1. (b) Areas and heights measured for horizons 1–4. (c) Area-depth diagram illustrating the area-depth relationships on both move lines and pin lines, with R2 values of >0.98.
Remotesensing 15 01603 g011
Figure 12. The inversion results of horizons 1–4 in Figure 9b show the distributions of the (a) Fault slip, (b) Trishear angle, and (c) P/S ratio. Red dashed lines indicate the values of best-fit parameters acquired from inverse modeling. See the text for a detailed analysis.
Figure 12. The inversion results of horizons 1–4 in Figure 9b show the distributions of the (a) Fault slip, (b) Trishear angle, and (c) P/S ratio. Red dashed lines indicate the values of best-fit parameters acquired from inverse modeling. See the text for a detailed analysis.
Remotesensing 15 01603 g012
Figure 13. The inversion results of horizons 1–4 in Figure 10b show the distributions of the (a) Fault slip, (b) Trishear angle, and (c) P/S ratio. Red dashed lines indicate the values of best-fit parameters acquired from inverse modeling. See the text for a detailed analysis.
Figure 13. The inversion results of horizons 1–4 in Figure 10b show the distributions of the (a) Fault slip, (b) Trishear angle, and (c) P/S ratio. Red dashed lines indicate the values of best-fit parameters acquired from inverse modeling. See the text for a detailed analysis.
Remotesensing 15 01603 g013
Figure 14. The inversion results of horizons 1–4 in Figure 11b show the distributions of the (a) Fault slip, (b) Trishear angle, and (c) P/S ratio. Red dashed lines indicate the values of best-fit parameters acquired from inverse modeling. See the text for a detailed analysis.
Figure 14. The inversion results of horizons 1–4 in Figure 11b show the distributions of the (a) Fault slip, (b) Trishear angle, and (c) P/S ratio. Red dashed lines indicate the values of best-fit parameters acquired from inverse modeling. See the text for a detailed analysis.
Remotesensing 15 01603 g014
Table 1. Model associated deformation estimates, age constrains, and calculated shortening rate.
Table 1. Model associated deformation estimates, age constrains, and calculated shortening rate.
TerraceDeformed Height (1σ)Undeformed Height (1σ)Surface Uplift (1σ)Horizontal Shortening (1σ)Age of Terrace (1σ)Horizontal Shortening Rate (1σ)
T493.5 ± 1.2 m40.9 ± 0.9 m56 ± 1.5 m81.7 ± 2.3 m51.7 ± 5.1 ka1.6 ± 0.2 mm/a
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, X.; Li, Z.; Wang, W.; Zhang, P.; Wu, C.; Chen, G.; Duan, L.; Wu, X.; Liu, K. Quaternary Crustal Shortening of the Houyanshan Structure in the Eastern Chinese Tian Shan: Constrained from Geological and Geomorphological Analyses. Remote Sens. 2023, 15, 1603. https://doi.org/10.3390/rs15061603

AMA Style

Yang X, Li Z, Wang W, Zhang P, Wu C, Chen G, Duan L, Wu X, Liu K. Quaternary Crustal Shortening of the Houyanshan Structure in the Eastern Chinese Tian Shan: Constrained from Geological and Geomorphological Analyses. Remote Sensing. 2023; 15(6):1603. https://doi.org/10.3390/rs15061603

Chicago/Turabian Style

Yang, Xue, Zhigang Li, Weitao Wang, Peizhen Zhang, Chuanyong Wu, Gan Chen, Lei Duan, Xiancan Wu, and Kang Liu. 2023. "Quaternary Crustal Shortening of the Houyanshan Structure in the Eastern Chinese Tian Shan: Constrained from Geological and Geomorphological Analyses" Remote Sensing 15, no. 6: 1603. https://doi.org/10.3390/rs15061603

APA Style

Yang, X., Li, Z., Wang, W., Zhang, P., Wu, C., Chen, G., Duan, L., Wu, X., & Liu, K. (2023). Quaternary Crustal Shortening of the Houyanshan Structure in the Eastern Chinese Tian Shan: Constrained from Geological and Geomorphological Analyses. Remote Sensing, 15(6), 1603. https://doi.org/10.3390/rs15061603

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