New Detrital Apatite Fission Track Thermochronological Constraints on the Meso-Cenozoic Tectono-Thermal Evolution of the Micangshan-Dabashan Tectonic Belt, Central China

The Micangshan-Dabashan tectonic belt, located in the southern Qinling-Dabie Orogen near the northeastern Tibetan Plateau, is a crucial area for understanding the processes and mechanisms of orogenesis. Previous studies have been focused on the cooling process via thermochronology and the mechanism and process of basement uplift have been investigated. However, the coupling process of basement exhumation and sedimentary cap cooling is unclear. The tectono-thermal history constrained by the detrital apatite fission track (AFT) results could provide valuable information for understanding crustal evolution and the coupling process. In this study, we provided new detrital AFT thermochronology results from the Micangshan-Dabashan tectonic belt and obtained nine high-quality tectono-thermal models revealing the Meso-Cenozoic cooling histories. The AFT ages and lengths suggest that the cooling events in the Micangshan area were gradual from north (N) to south (S) and different uplift occurred on both sides of Micangshan massif. The cooling in Dabashan tectonic zone was gradual from northeast (NS) to southwest (SW). The thermal histories show that a relatively rapid cooling since ca. 160 Ma occurred in the Micangshan-Dabashan tectonic belt, which was a response to the event of Qinling orogenic belt entered the intracontinental orogenic deformation. This cooling event may relate to the northeastward dextral compression of the Yangtze Block. The sedimentary cap of Cambriano-Ordovician strata responded positively to this rapid cooling event and entered the PAZ since ca. 63 Ma. The deep buried samples may be limited affected by climate and water erosion and the accelerated cooling was not obvious in the Late Cenozoic. Collectively, the cooling processes of basement and sedimentary cap in Micangshan-Dabashan tectonic belt were inconsistent. The uplift of the sedimentary area is not completely consistent with that of the basement under thrust and nappe action. The rigid basement was not always continuous and rapidly uplifted or mainly showed as lateral migration in a certain stage because of the different intensities and modes of thrust and nappe action, and the plastic sedimentary strata rapidly uplifted due to intense folding deformation.

The Micangshan-Dabashan tectonic belt, located in the southern Qinling-Dabie Orogen near the northeastern Tibetan Plateau, is a crucial area for understanding the processes and mechanisms of orogenesis. Previous studies have been focused on the cooling process via thermochronology and the mechanism and process of basement uplift have been investigated. However, the coupling process of basement exhumation and sedimentary cap cooling is unclear. The tectono-thermal history constrained by the detrital apatite fission track (AFT) results could provide valuable information for understanding crustal evolution and the coupling process. In this study, we provided new detrital AFT thermochronology results from the Micangshan-Dabashan tectonic belt and obtained nine high-quality tectono-thermal models revealing the Meso-Cenozoic cooling histories. The AFT ages and lengths suggest that the cooling events in the Micangshan area were gradual from north (N) to south (S) and different uplift occurred on both sides of Micangshan massif. The cooling in Dabashan tectonic zone was gradual from northeast (NS) to southwest (SW). The thermal histories show that a relatively rapid cooling since ca. 160 Ma occurred in the Micangshan-Dabashan tectonic belt, which was a response to the event of Qinling orogenic belt entered the intracontinental orogenic deformation. This cooling event may relate to the northeastward dextral compression of the Yangtze Block. The sedimentary cap of Cambriano-Ordovician strata responded positively to this rapid cooling event and entered the PAZ since ca. 63 Ma. The deep buried samples may be limited affected by climate and water erosion and the accelerated cooling was not obvious in the Late Cenozoic. Collectively, the cooling processes of basement and sedimentary cap in Micangshan-Dabashan tectonic belt were inconsistent. The uplift of the sedimentary area is not completely consistent with that of the basement under thrust and nappe action. The rigid basement was not always continuous and rapidly uplifted or mainly showed as lateral migration in a certain stage because of the different intensities and

INTRODUCTION
The Micangshan-Dabashan tectonic belt, located at the northeastern margin of the Tibetan Plateau, is a key "basinmountain coupling" transitional orogenic belt that separates the Qinling orogenic belt to the north and the Sichuan Basin to the south ( Figure 1). Its tectonic evolution and thermal histories have attracted much attention due to its relevance in understanding the onset of the Asian monsoon, drainage adjustment of the Yangtze River system, growth of the Tibetan Plateau, evolution of the Qinling-Dabie Orogen, erosional lowering of the Sichuan Basin as well as for evaluating the new Frontier exploration potential of lower Cambrian and Neoproterozoic hydrocarbon reservoirs (including shale gas) (Clark et al., 2004;Enkelmann et al., 2006;Shen et al., 2007;Richardson et al., 2008;Tian et al., 2010;Tian et al., 2012;Yang et al., 2013;Yang et al., 2017). The Micangshan-Dabashan tectonic belt formed after collision with the South China Block (SCB) and North China Block (NCB) in the Indochina period and intracontinental orogeny in the Yanshanian-Himalayan period (Zhang et al., 1989;Meng and Zhang, 2000;Meng et al., 2005;Tan et al., 2007;Xu et al., 2009;Tian et al., 2010;Tian et al., 2012). The internal structure is complex, and the deformation characteristics are obviously different among different units. The tectonothermal evolutionary histories need to be addressed in a timely manner, especially given that the evolution of the Micangshan-Dabashan tectonic belt since the Meso-Cenozoic has remained controversial.
As one of the most important thermochronological methods, apatite fission track (AFT) parameters provide quantitative information on cooling and record the thermal history of rocks Green, 1988;Gallagher, 2003), and these parameters are typically interpreted as a result of surface uplift and erosion (Tian et al., 2010;Tian et al., 2012;Tian et al., 2013;Lei et al., 2012;Ren et al., 2015;Qi et al., 2016;Powell et al., 2017;. Previous work that has focused on the thermal history derived from a variety of thermochronometry dating methods, such as U-Pb, 40 Ar/ 39 Ar, (U-Th)/He and AFT, in the Longmenshan, Micangshan and Dabashan areas has been published, and the dating methods and results have mainly been applied to investigate the cooling events and the coupling process of the Tibetan Plateau uplift (Arne et al., 1997;Enkelmann et al., 2006;Shen et al., 2007;Richardson et al., 2008;Li ZW et al., 2010;Li ZW et al., 2012;Li JH et al., 2010;Li et al., 2011;Chang et al., 2010;Tian et al., 2010;Tian et al., 2012;Tian et al., 2013;Tian et al., 2016;Tian YT et al., 2018;Yang et al., 2013;Yang et al., 2017;Shen et al., 2019). The thermochronometry dating and thermal histories results above and others from the, Tongbai-Dabieshan (Webb et al., 1999;Reiners et al., 2003;Hu et al., 2006), Taibaishan  and Western Qinling (Zheng et al., 2004;Enkelmann et al., 2006) show that the regional uplift occurred from the Late Jurassic to the Early Cretaceous and an accelerated uplift since the Late Cenozoic which related to the growth of the northeastern Tibetan Plateau. The Longmenshan tectonic defines the eastern margin of the eastern part of the Tibetan Plateauis and the western boundary of the Micangshan tectonic zone. The Longmenshan experienced multiple phases of intracontinental deformation during Mesozoic-Cenozoic (Tian et al., 2013;Tian et al., 2016;Tian YT et al., 2018;Shen et al., 2019), such as the Early Cretaceous (ca.119-131 Ma) and the Late Cenozoic (ca.30-25 Ma, 20 Ma, and 12-5 Ma, Arne et al., 1997;Kirby et al., 2002;Godard et al., 2009;Wang et al., 2012), and these tectonic activities may relate to the uplift and denudation in Micangshan tectonic zone.
Generally, the previous low-temperature thermochronological data was focused on the uplift of basement and show that the timing of the uplift in different areas of the Micangshan-Dabashan tectonic belt during the Mesozoic to Cenozoic was different. Furthermore, the apatite fission track (AFT) ages in the same area were inconsistently affected by factors such as experimental conditions and the nature of the samples (Figure 1). In contrast to the abundant thermochronology studies of the basement uplift, few studies have been conducted in the thermal evolution of sedimentary cap and its relationship with the uplift of basement. During the orogenic process, the uplift and denudation of the mountain were closely related to deposition in the nearby basin, which was a coupling process. The basin sediments adjacent to the orogenic belt recorded abundant uplift and denudation information during the late orogenic process. Therefore, the detrital AFT collected from the sedimentary basin could reveal its thermal history and the exhumation cooling history of the orogenic belt in the provenance area (Shen et al., 2005;Homke et al., 2010;Lin et al., 2015;Zhang et al., 2016;Du et al., 2018).
Detrital AFT analysis in the hinterland of the Micangshan and Dabashan areas is rare, and there is a lack of information regarding the tectono-thermal history of the sedimentary area. There are some differences in the results of previous thermochronometry studies related to the timing of cooling events in Micangshan and Dabashan. Therefore, new detrital AFT data and analyses need to be provided as key evidence for the timing of cooling events. In this study, we present a lowtemperature thermochronology study of the Micangshan and Northern Dabashan, employing apatite fission track (AFT) data from nine detrital samples to assess the regional cooling history and providing information about the coupling relationship of uplift in basement and sedimentary areas.

GEOLOGICAL SETTING
The Micangshan-Dabashan tectonic belt is located at the junction of the Qinling orogenic belt and Sichuan Basin (Figures 2A,B), which was a passive continental margin in the early Paleozoic (Tian et al., 2010;Yang et al., 2013). The South China Block subducted under the North China Block beginning in the late Caledonian, underwent full collision until the Late Triassic, and experienced continuous compression and reconstruction in the Yanshan-Himalayan. Structural phenomena such as continuous fold deformation, differential uplift and thrust-slip faults are widely developed (Shi and Shi, 2014, Figures 2C,D).
The Dabashan tectonic belt is an important part of the multilayered thrust nappe tectonic system on the southern margin of the Qinling Orogen. The Late Triassic collision of the South Qinling and Yangtze Block along the Mianlue suture and the large-scale intracontinental subduction of the Yangtze Block beneath the Qinling Orogen during the Yanshanian period (Cheng et al., 2004;Huang and Wu, 1992) resulted in extensive fold and thrust deformation. Usually, the Dabashan tectonic belt is divided into two different tectonic units and sedimentary systems, the Southern Dabashan (SDBS) and the Northern Dabashan (NDBS), and they are bounded by the Chengkou fault (He et al., 1997;Li et al., 2011Yue, 1998, Figure 2E). The Dabashan thrust belt initiated and propagated southwards into the northern Sichuan Basin, which resulted in the formation of the Dabashan foreland basin . A number of NW-SE-trending faults are developed, and the Chengkou fault is the main active surface stacking the thrusts in the SW direction Dong et al., 2008). The thrust belt mainly comprises Proterozoic basement and Sinian-Silurian limestone-dominated cover, which are all intruded by Silurian mafic Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 754137 4 dikes (Zou et al., 2011;Li RX. et al., 2012). Devonian to Middle Triassic strata are exposed only north of the Ankang fault, and lower Paleozoic strata occur along faults and unconformably contact rare Lower to Middle Jurassic strata. Micangshan is located on the northwestern margin of the Upper Yangtze, adjacent to the Longmenshan tectonic belt in the west and to the Southern Dabashan foreland fold-thrust belt in the east. The main structural lines in the area spread out in the NE direction, and most of them behave similarly to thrust faults and slip faults ( Figure 2). The intersection with the South Qinling block shows a strong dextral arc, with a nearly N-S-trending compound anticline. The Archean-early Proterozoic granitic basement and a small amount of basic intrusive rocks are exposed in the central uplift area of Micangshan, in which there are relatively intact residual synclines, mainly the Permian-Triassic system and the Sinian to Silurian system on the two wings. The overlying strata experienced deformation during the later stage. The lithology of outcrops in the regions varies after uplift and denudation Li et al., 2008;Tian YT et al., 2018;Tian et al., 2019). From north to south in the Micangshan tectonic zone, a distinct succession of Neogene-Quaternary sedimentary (Hanzhong Basin), crystalline basement (Hannan and Micangshan massifs), Sinian-Triassic sedimentary (Huijunba syncline) and Sinian-Cretaceous sedimentary cover (Norther Sichuan Basin) can be identified ( Figure 2D).

SAMPLES AND METHODS
A total of nine outcrop samples were collected for apatite fission track analysis, three of which were from the vicinity of the Huijunba syncline in the Micangshan tectonic belt, and six were from the Dabashan tectonic belt ( Figure 2C). All the samples were from lower Cambrian to Lower-Middle Ordovician siltstone and sandstone.
Apatite was separated from each sample using standard heavy liquid and magnetic separation techniques. Apatite grains were mounted in epoxy resin on glass slides and then ground and polished to an optical finish to expose internal grain surfaces. Spontaneous fission tracks in the apatite were etched in 6.6% HNO 3 at 25°C for 30 s to reveal the spontaneous fission tracks. Low-U muscovite in close contact with these grains served as an external detector during irradiation (Yuan et al., 2003;Yuan et al., 2006). Neutron fluence was monitored in CN5 U dosimeter glasses (Bellemans et al., 1995). After irradiation in the reactor, the external muscovite detector was detached and etched in 40% HF for 20 min at room temperature to reveal induced tracks. Track densities of both spontaneous and induced fission track populations were measured with a dry objective at ×10015 magnification. Fission track ages were calculated using the IUGS-recommended zeta calibration approach (Hurford, 1990), which is shown in Eq. 1. Zeta values used in this study were determined from repeated measurements of standard apatites (Hurford and Green, 1983). The weighted mean zeta value for apatite used by the fission track operator was 391 ± 17.8 a/cm 2 according to the calibration of sample standards. The lengths of the horizontally confined fission tracks were measured exclusively in prismatic apatite crystals because of the anisotropy of annealing of fission tracks in apatite .
where T sample is the AFT age of the samples; λ d , which is 1.55125 × 10 −10 a −1 , is the decay constant of 238 U; ρ s is the track density for spontaneous fission tracks; ρ i is the track density for induced fission tracks; ρ d is the track density for the CN5 glass standard; and ξ is the Zeta constant, which is 391 ± 17.8 in this study. Vitrinite reflectance can provide an accurate temperature range for fission track annealing and evidence for thermal history modelling (Sweeney and Burnham, 1990;George et al., 2001). The equivalent vitrinite reflectance calculated by bitumen reflection (Feng and Chen, 1988) in the Lower Camburian Niutitang shale is in the range of 1.86-2.54% , and the range of 2.98-3.03% in the Lower Cambrian Jianzhuba mudstone which is adjacent to the Lower-Middle Ordovician in NDBS (Table 1), reflecting exmaximum paleotemperatures of 233.4°C-273.3°C and 294°C-295°C according to Barker and Pawlewicz (1986). These results demonstrate that fission tracks in apatite were completely annealed in a geological episode and that the apatite fission tracks were all new tracks through the partial annealing zone (PAZ).
The apparent age of the AFT has no direct geological relevance (Gleadow et al., 2002). However, the thermal history simulation combined with the length distribution of the fission track and the annealing kinetic parameters of apatite itself can well reveal the thermal history experienced by the sample. To achieve robust thermal histories, inverse modelling of AFT data was carried out using the programme HeFTy (version 1.8.3, Ketcham, 2014). Model paths between specified time-temperature constraints were chosen to be monotonic rather than monotonically consistent in HeFTy. This choice provides more freedom in model paths between constraint points. Dpar was used as a kinetic parameter, and thermal histories were calculated using a multikinetic annealing model. The Kolmogorov-Smirnov test function was used to compare measured length histograms to model results (Ketcham et al., 2007). The initial mean track length was calculated based on Dpar (Carlson et al., 1999). The number of fitting curves was set to 20,000. The model quality was judged by the goodness of fit (GOF). Generally, a GOF value greater than 0.05 indicates acceptable simulation results, while a GOF value greater than 0.5 indicates high-quality simulation results (Ketcham, 2005).

Apatite Fission Track Ages and Lengths
The fission track test results of all samples are shown in Table 2. The AFT ages of the sandstone and siltstone range from 34 Ma to 50 Ma, which are much younger than the ages of the corresponding sedimentary formations. These results indicate that the samples were all affected by thermal events and annealed after formation. Except for the F-B-4 sample, the other eight samples all pass the age χ 2 -test (P (χ 2 )>5%), Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 754137 5 Tian et al.
Tectono-Thermal Evolution of Micangshan-Dabashan indicating that the age difference of a single particle is within the statistical error range and is the same age component. The pooled age is generally used as the AFT age. The sample age χ 2 -test of F-B-4 is less than 5%, representing the mixed age of different provenance minerals, and the central age can be used to interpret the cooling event (Galbraith and Laslett, 1993;Sobel et al., 2006). The samples collected from the Micangshan area have a short north-south axial distance, and the AFT ages mainly range from 40 Ma to 49 Ma. The ages of the Z-7, Z-8, P-1, and L-1 samples in the northeastern part of the Dabashan belt range from 44 Ma to 50 Ma, which is considerably older than that of sample Z-20 in the southwest with an AFT age of 34 Ma. This difference may indicate that the Dabashan thrust nappe belt was uplifted from NE to SW. The average measured AFT lengths of the samples are in the range of 11.7-13.2apatite fission trackμm, all of which are less than the original track lengths (generally exceeding 16 μm, Carlson et al., 1999), indicating that the samples underwent strong annealing. The distribution pattern of track lengths is generally unimodal (Figure 3A), and the peak values are mainly distributed in the range of 13-15 μm, indicating that the samples did not experience a multistage thermal history but were in the process of monotonic cooling , which may reflect rapid stripping and cooling during the later orogeny.
N, induced number of apatite grains; ρ s, ρ i and ρ d , track densities for spontaneous fission tracks, induced fission tracks and CN5 glass standard tracks, respectively; Ns, Ni and Nd, number of spontaneous fission tracks counted, induced fission tracks counted and CN5 glass standards, respectively; P (χ 2 ), chisquare probability; L, mean track length; n, number of tracks measured.
Collectively, Dpar is a good potential indicator of fission track annealing kinetic properties in natural apatite grains (Carlson et al., 1999). The Dpar content is a routinely used parameter, exhibits a strong and informative correlation with apatite fissiontrack ages and lengths ( Figure 3B). It should be noted that the lithological heterogeneity and chemical composition of samples might be different, for example the Dpar with diameters of 1.8-5.0 μm (Richardson et al., 2008), resulting in uncertainty and a negative impact on the basin-scale exhumation and tectono-thermal processes. In this study, the Dpar values (mean etch pit size parallel to the c axis, Donelick, 1993) of the samples range from 1.41 to 1.79 μm ( Figure 3C). No significant Dpar variation was observed among the samples reported in this work, in which the single-grain ages and length came from the same kinetic population. The Dpar results suggest that the differences in AFT age and length are due to the differences in thermal history and may not influenced by nature of apatite itself.
There is a positive correlation between the fission track age and the elevation of the samples, which reflects that the samples in high-altitude localities always passed through the annealing zone of apatite and produced tracks (Shen et al., 2007). However, the study area at the junction of the Qinling orogenic belt and Sichuan Basin, coupled with continuing subduction collision since the Indo-Chinese epoch and intracontinental orogeny, caused the large-scale thrust nappe Shen et al., 2007;Chang et al., 2010). The elevation and distribution of apatite often change due to thrust nappes, uplift, denudation and faulting. Therefore, there is often no correlation between the fission track age and elevation in complex tectonic belts Shen et al., 2007;Qi et al., 2016).
Statistical analysis of fission track ages and elevations in the study area and its surrounding structural belt indicates that it underwent complex tectonic processes (Figure 4). The fission track ages of some sections have a good corresponding relation with altitude, but the slopes are quite different, which shows that the tectonic uplift effect is obviously different in the area. The different uplift effects may be related to the peripheral structural belt of acute Himalayan activity, especially with the far-field effect of the rapid expansion of regional topography and strong reformation from the northeastern margin of the Qinghai-Tibet Plateau (Tian et al., 2010;Tian et al., 2012;Yang et al., 2013, Yang et al., 2017. In this study, there was no linear correlation between AFT age and elevation. In addition to the complex geological process, this lack of correlation may also be related to the concentrated distribution of the topographic elevation of the samples in the smaller ranges of 391-658.2 m and 817-1442.2 m, which are insufficient to reflect the linear relationship between AFT ages and elevation. In terms of the relationship between sample lithology and age, the ages of the samples from the outcropping area of the bedrock are generally higher than those of the detrital rocks in the sedimentary area (Figure 4). The AFT ages of the detrital rocks are mainly distributed in the range of 28.3-85 Ma, and those of the bedrock are mainly distributed in the range of 71-123.5 Ma. As old rigid blocks with weak detachment, the bedrock outburst area in the study area and its surrounding area responds more quickly to tectonic activities such as thrust nappes or intracontinental orogenies that cause formation uplift, and the apatite distributed in the area enters the PAZ earlier and resets its age. Micangshan-Dabashan was a passive continental margin in the early Paleozoic, and the samples from the formation were typical marine sediments. The Micangshan massif formed in the Late Triassic (Tian et al., 2010), and the provenance of the samples was not from this massif. Even if there are other provenance inputs, thermal relocation occurs with increasing burial depth and temperature after the detrital particles are transported and deposited. Therefore, the clastic apatite fission tracks reflect the late tectonic and thermal evolution of the sedimentary area rather than the provenance area. The simulation of the thermal evolutionary history of wells HY-1 and SND-1 in the study area indicates that the maximum paleogeotemperature exceeded 220°C at 160 Ma. The Cambrian-Ordovician system at 123.5-71 Ma was buried at a depth of more than 4,000 m and had high paleo-geotemperatures ranging from 120°C to 190°C (Tian et al., 2020), and apatite did not enter the annealing zone. Then, the apatite from the Cambrian-Ordovician system entered the annealing zone with subsequent rapid uplift, so the uplift and cooling events recorded by the fission track ages were relatively late. Overall, the Micangshan-Dabashan tectonic belt experienced uplift and denudation from at least 123.5 Ma and affected the surrounding sedimentary strata, entering the PAZ from 85 Ma to 28.3 Ma.
The AFT age and length are positively correlated to a certain extent according to the data in the research area ( Figure 5). The AFT ages gradually become younger from the Micangshan-Hannan massif to the northern Sichuan Basin, indicating a gradual thrust expansion from north to south in the region as piggyback expansion (Tian et al., 2010;Yang et al., 2013;Wang et al., 2003). In this tectonic framework, the AFT ages of the Cambrian strata samples in the Huijunba syncline should be larger than those in the northern Sichuan foreland area. However, the AFT ages of the Cretaceous-Cambrian stratum samples in the northern Sichuan foreland area range from 61.6 to 85 Ma with the track lengths between 10.6 and 12.2 μm and that in Huijunba syncline at the range of 40-49 Ma with track lengths between 11.7 and 13.2 μm ( Figure 5A). This distribution of AFT ages and FIGURE 4 | Correlation between AFT ages and elevation.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 754137 8 lengths seems inconsistent with regional tectonic evolution history. Micangshan tectonic zone is a large complex anticline formed by South Qinling collision orogeny since the Triassic and is characterized by the asymmetric structure with a gentle north wing orienting Huijunba syncline and steep south wing orienting north Sichuan foreland (Zhang and Dong, 2009;Zhang, 2019). The samples in the Sichuan foreland area was deeper than samples in the north Huijunba syncline. Therefore, these abnormal AFT ages may be caused by differential basement uplift, as the sedimentary cover in the northern Sichuan foreland area uplifted and entered the PAZ earlier than that in the Huijunba syncline during the overall N-S trending uplift process.
Generally, the AFT age from Dabashan area gradually becomes younger from NE to SW, indicating that the timing of uplift in the NE region was earlier than that in the SW (Shen et al., 2007;Li JH. et al., 2010;Yang et al., 2017, Figure 5B). In this study, the AFT age of the NDBS ranges from 35 Ma to 50 Ma, the youngest age of NDBS sample Z-20 is located near the Chengkou fault. Collectively, the AFT ages from the NDBS are older than most ages from the SDBS ( Figure 5B). This distribution of AFT ages indicates that the uplift process of the sedimentary cap is characterized by earlier to NE and later to SW under the influence of regional uplift. This interpretation is consistent with the thrust nappe of Dabashan from NE to SW (Shen et al., 2007) and the NE-SW trending paleotectonic stress during the Late Jurassic to the Early Cretaceous Dong et al., 2014). The analysis of AFT provides a new thermochronological constraint for the progressive extension deformation of the thrust nappe structure in the Dabashan area.

Tectono-Thermal History of the Micangshan-Dabashan Tectonic Belt
Geological background constraints are very important for the inversion of the AFT thermal history (Tian et al., 2010). In the study area, many major geological events occurred during the Jurassic and Cretaceous, and the Qinling orogenic belt truly entered the intracontinental orogenic deformation period at approximately 160 Ma . The Micangshan-Dabashan tectonic belt, as the southern margin of the Qinling orogenic belt, likely responded to this tectonic event. Previous AFT data and thermal simulation results showing that Dabashan area experienced uplift at 180-165 Ma, 170-160 Ma and ca.100 Ma (Shen et al., 2007;Li JH et al., 2010;Xu et al., 2010), and the Micangshan area experienced uplift at 140-100 Ma and 152-150 Ma (Tian et al., 2010;Tian et al., 2012;Xu et al., 2010), which have also provided important constraints for the thermal simulation.
To constrain the maximum paleo-geotemperature experienced by the samples, paleotemperature scales such as vitrinite reflectance are always used along with homogenization temperatures of fluid inclusions in or adjacent to the strata where the samples were located ( Table 1). In the Micangshan area, the inclusion samples collected from the Lower Cambrian Niutitang Formation were emplaced at 212-166 Ma (Tian et al., 2020). The organic fluids in the Dabashan areas were captured in the Late Triassic to Early-Middle Jurassic (Li RX. et al., 2012). Therefore, the vitrinite reflectance and homogenization temperature of organic fluids were used in this study to constrain the maximum paleogeotemperature in the Late Triassic to Early-Middle Jurassic. The F-B-3 and F-B-4 samples are not representative because of the effect of regional thermal anomalies (Tian et al., 2020), and the AFT age χ 2 -test is less than 5%, which reflects the thermal history of different provenances. Thermal history inversion simulations were carried out for samples other than F-B-3 and F-B-4 under the geological constraints in this paper.
The inversely modelled t-T paths for the samples collected from the Micangshan-Dabashan tectonic zone are shown in Figure 6.  Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 754137 9 and Dabashan tectonic zone at 153 Ma , which is a precise response to the event of Qinling orogenic belt entered the intracontinental orogenic deformation period Li et al., 2013). However, this cooling time is much earlier than that in the southwestern Qinling (Late Cretaceous, Enkelmann et al., 2006), east Qinling-Dabieshan (ca.100 Ma, Hu et al., 2006) and Longmenshan tectonic zone (119-131 Ma, Arne et al., 1997). The previous work shows that it is possibly a response to the combined effects of the Mesozoic tectonism in nearby regions, which probably reactivated pre-existing structures through strike-slip faulting (Webb et al., 1999;Hu et al., 2006;Tian et al., 2012). The Yangtze Block rotated clockwise relative to the north China block through Mesozoic time (Sun et al., 2006) may be an important inducement. The Micangshan-Dabashan area is the core area of the deformation of rotated clockwise that uplift and cooling earlier due to the dextral extrusion. And then spread to the further and more marginal areas such as southwestern Qinling, eastern Qinling-Dabie and Longmenshan tectonic zone. Furthermore, the rapid cooling of Micangshan-Dabashan area at ca.160 Ma was probably the northern source for the early Cretaceous deposition in the northern Sichuan Basin (Meng et al., 2005;Tian et al., 2012). On the background of uplift in the region, the basement and sedimentary cap in Micangshan-Dabashan area also underwent different uplifting processes.
The Micangshan area was uplifted as a whole during 160-90 Ma, and the sedimentary strata and basement experienced a relatively consistent regional uplift process. The sedimentary strata uplifted differently from the basement since 90 Ma. Under the action of thrust and nappe tectonics, the rigid basement uplifted slowly and horizontally, and the plastic strata were squeezed and deformed, which produced the characteristics of folding and uplift since 90 Ma and then entered the PAZ at 63 Ma. Previous study concluded that the displacement of upper structure in Micangshan area is larger than the lower during the progressive deformation process of contraction orogenesis from north to south based on the analysis of structural types and geometric evolution characteristics of folds (Ramsay and Lisle, 2000;Zhang, 2019). And our finding of the different uplifting processes of basement and sedimentary cap is comparable with previous results. In the Dabashan area, the uplift of sedimentary strata was different from that of the basement during 165-95 Ma. The basement uplifted relatively slowly, while the sedimentary strata uplifted rapidly due to strong fold deformation. The sedimentary strata and the basement experienced a relatively consistent regional uplift process from 95 Ma to 65 Ma. After 65 Ma, the strata experienced differential uplift again ( Figure 6). In the process of Dabashan nappe extrusion, the rapid uplift of the lower Paleozoic strata in the NDBS since 160 Ma occurred much earlier than that of the Paleozoic strata in the SDBS since ca. 100 Ma (Shen et al., 2007), which reflects the process of nappe uplift from north to south in the Dabashan area. The occurrence of differential uplift between sedimentary strata and basement may reflect the different intensities of thrust and nappe in different stages.
Abundant AFT ages of Late Cenozoic have been used to explain the rapid northeastward growth of Tibetan Plateau, such as ca.20 Ma in northern Longmenshan, 13-15 Ma in central and south Longmenshan and ca.13 Ma in the eastern Tibetan Plateau. Most of the samples which revealed the accelerated cooling were exposed basement samples and shallow buried Cretaceous-Permian samples. However, the simulation results in this study do not show the obvious accelerated cooling in the Late Cenozoic. The accelerated cooling of northeastern Tibet Plateau since Late Cenozoic was the result of multiple factors, which was not only the eastward growth of the Tibetan Plateau, but also intensified by climate change in Southeast Asia (Molnar and Tapponnier, 1975;Tian et al., 2012) and the erodibility of the drainage systems that accelerated the denudation (Richardson et al., 2008;Richardson et al., 2010). Our new samples were from the deep Cambrian-Ordovician strata and they are still buried deeply lager than 2 km (Tian et al., 2020) and was less affected by climate and water erosion in the Late Cenozoic. These may be the reasons why the cooling process was not significantly accelerated in this thermal simulation since the Late Cenozoic.

CONCLUSION
The detrital AFT thermochronological dataset provide new insight into the Meso-Cenozoic tectono-thermal evolution of the Micangshan-Dabashan tectonic belt, central China. New evidence leads us to draw the following conclusions: 1) The AFT ages and lengths suggest that the cooling events in the Micangshan area were gradual from north to south and different uplift occurred on both sides of Micangshan massif. The cooling in Dabashan tectonic zone was from northeast to southwest. 2) A relatively rapid cooling since ca.160 Ma occurred in the Micangshan-Dabashan area which was a precise response to the event of Qinling orogenic belt entered the intracontinental orogenic deformation period. The cooling event may relate to the northeastward dextral compression of the Yangtze Block.

3) The cooling processes of basement and sedimentary cap in
Micangshan-Dabashan tectonic belt were inconsistent under the thrust and nappe action. The rigid basement was not always continuous and rapidly uplifted or mainly showed as lateral migration in a certain stage and the plastic sedimentary strata rapidly uplifted due to intense folding deformation and entered the PAZ since ca. 63 Ma. 4) The deep buried samples in Micangshan-Dabashan tectonic belt may be less affected by climate and water erosion what intensified the cooling so that the accelerated cooling was not obvious in the Late Cenozoic as usual.

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