Reconstruction of the proto-type basin and tectono-paleographical evolution of Tarim in the Cenozoic

The Tarim Basin is the largest sedimentary basin in China, and it has experienced a complex tectonic evolutionary history. Reconstruction of the proto-type basin and tectono-paleogeography is helpful to understand the different stages of evolution of the sedimentary basin and basin-mountain relationship. It is significant to combine the basin with the regional tectonic background to discuss the process of basin-mountain coupling and the tectonic evolution of the peripheral orogenic belts. With a reliable residual thickness map and lithofacies map of the Tarim Basin in the Cenozoic, based on the amount of shortening we quantified from previous works and 81 balanced cross-sections, we restored the original range and compiled the proto-type basin map of Tarim Basin. From a compilation of previous studies on the lithofacies of peripheral blocks, the tectono-paleogeography of the Tarim Basin in Cenozoic has been reconstructed. The Indian Plate collided with the Eurasian continent at ∼45–40 Ma. The remote effect of the collision led to the resurrection and reactivation of the Kunlun and Tianshan Mountains. The Southwest Tarim and Kuqa rejuvenated foreland basins separately developed along the north front of the Kunlun Mountains and the south front of the Tianshan Mountains. The tectonic evolution process of the Tarim Basin in the Cenozoic was divided into two stages: 1) in the Paleogene, the Neo-Tethys Ocean retreated stepwise westward from the Southwest of the Tarim Basin, and the sedimentary lithofacies of the Southwest Tarim Depression were bay lagoon facies and lake facies; 2) the Neo-Tethys Ocean retreat finally occurred in the Tarim Basin during the Late Oligocene to Early Miocene, and it became an almost closed terrestrial basin, with the deposition of fluvial facies and lacustrine facies. The Cenozoic tectono-paleogeography of the Tarim Basin is closely related to the closure of the Neo-Tethys Ocean and the reactivation of the Kunlun and Tianshan Mountains.


Introduction
The Tarim Basin is the largest intracratonic basin in China and covers an area of 560,000 km 2 . Source rocks for oil and gas were deposited widely in the Tarim Basin during different geological periods, including the Paleogene, and excellent oil and gas can usually be found in its marine rocks. The Tarim Basin is a natural laboratory with unique geological and geomorphological units, and can be divided into seventeen tectonic units in the Cenozoic (Dong et al., 2016;He et al., 2016) ( Figure 1). Geographically bounded by the Western Kunlun Mountains, the Altyn Tagh mountain belt, and the Tianshan Mountains, the Tarim Basin has experienced a complicated tectonic history (Allen et al., 1993;Zhang et al., 2011;He et al., 2016;Laborde et al., 2019) (Figure 1).
In the early Cenozoic, the India-Asia collision caused largescale intra-continental deformations (Liu et al., 2009;Fang et al., 2020). Under the far-field effect of the collision, a prominent tectonic salient formed on the Pamir plateau and the Pamir Salient experienced some northward translation (Lin et al., 2019;Li Y. P. et al., 2020). A major reactivation of the Western Kunlun Range occurred during the Cenozoic in response to the India-Asia collision (Tapponnier and Molnar, 1977;Laborde et al., 2019). A major reactivation of the Tian Shan also occurred during the Cenozoic in response to the India-Asia collision and this reactivation seems to have been initiated in the Oligocene with subsequent Miocene accelerations, which is supported by evidence including sediment provenance and growth strata (Izquierdo-Llavall et al., 2018;Li et al., 2019). The Altyn Tagh Range was also reactivated during the Cenozoic due to the India-Asia collision and this reactivation seems to have been initiated in the Eocene with subsequent Oligocene to Miocene accelerations (Cheng et al., 2015;Zhao et al., 2016). The maximum principal stress direction in different areas of the Tarim Basin is diverse. In the north of the Tarim Basin, the maximum principal stress direction is mainly N-S and NNW-SSE (Huang et al., 2013). In the southwest, the maximum principal stress direction is mainly NE-SW; and in the southeast, the maximum principal stress direction is mainly NW-SE (Wang et al., 2007).
A schematic structural map and seismic cross-sections was constructed to document the Cenozoic deformation through the Tarim area (Figures 1, 2). In the Cenozoic, the Tarim Basin was in a stage of strong compression, and the continuous collision around the Tarim Block accommodated a large amount of Cenozoic compressive deformation (Li et al., 2007;Laborde et al., 2019). Ancient orogens around the Tarim Basin were rejuvenated, and a thick (up to~10 km) and continuous Cenozoic sedimentary series, mainly supplied by the erosion of the surrounding mountain ranges (Wei et al., 2000;Jia, 2009;He et al., 2013), formed the rejuvenated foreland basin in Kuqa and Southwest Tarim (Laborde et al., 2019) and the complex fault system in Bachu Uplift (Ning et al., 2021;Zhang et al., 2021) (Figure 2). In the case of the Western Kunlun and Southwestern Tian Shan compressive systems, these foothills mainly formed the Basin-Mountain system at the edge of the Tarim Basin ( Figure 2). Wide (from~50 tõ 150 km) and long (100 s of km) uplifts, such as the Bachu uplift, are visible to the west of the basin (Figure 2).
The far-field effect caused the reactivation of the Paleozoic and Mesozoic tectonic belts, changing the sea-land configuration in eastern Asia (Tapponnier et al., 1981;Wang et al., 2014;Lin et al., 2019;Li Y. P. et al., 2020). The Neo-Tethys Ocean retreated from Western Tarim in response to a combination of eustatic sea-level lowering, plate tectonics, and the autocyclic Four seismic cross-sections of the Tarim Basin (see Figure 1 for locations).

Frontiers in Earth Science
frontiersin.org 03 sedimentary infilling of the basin (Bosboom et al., 2011;Bosboom et al., 2014). At least five different seawater transgression-regression cycles occurred and the sea retreated from different parts of the Tarim Basin (Sun and Jiang, 2013;Bosboom et al., 2014;Sun et al., 2016a). The last regression marking the Tarim Basin generated a huge transformation from marine facies to terrestrial facies.
As a result of the collision of the Indian Plate and the Eurasian continent, the Tarim Block has different types of sedimentary facies from different stages of the Cenozoic. There are two issues to be addressed in this study: the tectono-paleogeographic evolution of the Tarim Basin and the tectonic interrelationships between the basin and the orogenic belts during the different stages. There are some studies that have mapped the lithofacies paleogeography of the Tarim Basin; however, several studies considered the whole basin and its periphery, and few compiled the prototype basin by considering the India-Asia collision and the huge shortening of the basin. In this study, we restored the proto-type basin of the Tarim Basin and reconstructed the Tectonopaleogeographic evolution around the Tarim Basin in the Cenozoic, analyzing the changes of the paleo-geographical environment and Cenozoic deformation of the Tarim Basin and its periphery. When and how the retreat of the Neo-Tethys Ocean occurred were specially considered to reconstruct the proto-type basin and tectono-paleogeography of the Tarim Basin and its periphery.

Methodology and database
Guided by three steps of map compilation, we reconstructed the proto-type basin and tectono-paleogeography of the Tarim Basin and its periphery in the Cenozoic.
First, we restored the original thickness map (Figure 3). Based on the analysis of wells and outcrops in the whole basin, the lithofacies-paleographic maps in the range of the present basin were reliable. If there were border facies in the basin, they could be used to help delineate the scope of the basin. If there were none, to recover the denudation area, it was necessary to reconstruct the original thickness map based on the Cenozoic residual thickness map with thickness trend analysis ( Figure 3). For an intracratonic basin, the thickness line needed to be concentric and closed. While for the foreland basin, the thickness line needed to be asymmetrical and the distance of each thickness line outside the basin at one-half of the inside. Therefore, the scope of the foreland basin was relatively narrow, and the terrane was near the Tarim Block.
Second, we restored the range of the proto-type basin, which meant restoring it to what it was before shortening or extension. The balanced cross-section method (Lou et al., 2016;Laborde et al., 2019) is one of the most practical ways to obtain the proto-type basin boundary. Based on previous balanced sections (Laborde et al., 2019), 71 balanced sections provided by Tarim Oilfield, and 10 balanced sections restored by us (Figure 4), we quantified the shortening and reconstructed the boundary of the pre-shortening proto-type basin.
Third, we restored the tectonic configuration of the Tarim Block and its adjacent plates during the development of the proto-type basin. After clarifying the sedimentation and deformation of the Tarim Basin, it was essential to clarify the peripheral margin of the Tarim Basin and the tectonic relationship between the Tarim Basin and its surrounding plates. To compile the proto-type basin maps and the tectono-paleogeographic maps, it was essential to figure out the answers to some questions: 1) the lithofacies distribution of the Tarim Basin in the Paleogene and Neogene (Guan and Guan, 2002;Shao et al., 2006;Yue et al., 2017;Ma CM. et al., 2020); 2) the paleogeographic distribution in its adjacent plates Ma T. et al., 2020;Song et al., 2022); 3) the sea retreat process of Neo-Tethys and its effect on the Tarim Basin (Bosboom et al., 2011;Bosboom, 2013;Sun and Jiang, 2013;Bosboom et al., 2014;Sun et al., 2016b;Li Q. et al., 2020;Sun et al., 2021).

Frontiers in Earth Science
frontiersin.org 04 Based on previously published data (Laborde et al., 2019) and our new data, we obtained the proto-type basin maps and the tectonopaleogeographic maps in the Cenozoic. It is noteworthy that the biggest difference between the proto-type basin map and the tectonopaleogeographic map is the latter includes larger geographical range, which provides more information with the tectonic configuration of the Tarim Basin and its adjacent plates.

Analysis
Before compiling proto-type basin maps, we clarified three issues: firstly, the division of the evolution stages in the Cenozoic; secondly, the marine sedimentation problem of the Tarim Basin in the Cenozoic; and thirdly, the quantified shortening amount of the Tarim Basin.   (Ding et al., 1993;Guo, 1994;Jia et al., 2004;Shao et al., 2006;Li et al., 2017).

Strata
Southwest Tarim

Division of the evolution stages in the Cenozoic
The strata of the Cenozoic Tarim Basin were mainly exposed around the basin, and most of the basin was covered by desert. According to the sedimentary characteristics of Cenozoic strata in the basin, the Tarim Basin could be divided into three regions: Southwest Tarim, North Tarim, and Southeast Tarim (Table 1). Although there were different division schemes for Cenozoic stratigraphic subdivision in this area, in this study, the evolution stages in the Cenozoic were based on the stratigraphic subdivision unified as follows (Table 1).
In Southwest Tarim, the Paleogene strata included the Bashibulake Formation, Zhuoyoulegansu Formation, Wulagen Formation, Kalatar Formation, Gaijitage Formation, Qimugen Formation, and Aertashi Formation in the Kashi Group; the Neogene strata included the Pakabulake Formation, Anjuan Formation, Keziluoyi Formation in the Wuqia Group and the Atushi Formation (Table 1). In North Tarim, the Paleogene strata included the Awate Formation, Xiaokuzibai Formation, and Talak Formation in the Kulugeliemu Group; the Neogene strata included the Kuche Formation, Kangcun Formation, Jidike Formation, and Suweiyi Formation (Table 1). In Southeast Tarim, the Paleogene strata included the Kulugeliemu Group; the Neogene strata included the Kuche Formation, Kangcun Formation, Jidike Formation, and Suweiyi Formation (Table 1).
Since the Paleogene, the tectonic compression caused the uplift of the surrounding mountains around the Tarim Basin, while the basin itself was relatively stable Laborde et al., 2019;Li Y. P. et al., 2020). During the Neogene, with the uplift of the Qinghai -Tibet Plateau, there are long-term growth and erosion in the mountain ranges and caused progressive thrusting over the Tarim basin (Jolivet et al., 2010;Wang et al., 2014;Laborde et al., 2019). Since the Quaternary, the Tarim Basin was surrounded by high mountains, progressively soil drying and forming China's largest inland water system and the largest desert.

Marine sedimentation problem of the Tarim Basin in the Cenozoic
According to the long-term eustatic trend of the Tarim Basin, there were at least five episodes of marine incursion in the Tarim Basin from the Cretaceous to the Paleogene ( Figure 5), which are supported by Cretaceous-Paleogene sedimentary records (Bosboom et al., 2011;Bosboom, 2013;Bosboom et al., 2014;Li et al., 2017;Zhang et al., 2018). There are reliable pieces of evidence to restrict the time of the fourth and fifth incursions and the sea retreat to the Cenozoic (Bosboom et al., 2011;Bosboom et al., 2014;Zhang et al., 2018) ( Figure 5).
The palaeomagnetic chronology and paleoenvironmental records of Cenozoic sections indicate that the fourth sea retreat from Aertashi occurred at~41 Ma (Bosboom et al., 2014). Research on marine-terrestrial sediments at Oytag in the Early Cenozoic also suggests that the sea retreat from the Tarim Basin may have been diachronous (Sun and Jiang, 2013). The sea retreat from Oytag occurred in the Early Eocene at~47 Ma and from Keliyang at 45 Ma (Sun and Jiang, 2013;Sun et al., 2016a); the~45 Ma peak in detrital zircon U/Pb was considered as an indication of a new sedimentary source (Wang et al., 2021).

Frontiers in Earth Science
frontiersin.org Pieces of evidence, including the palaeomagnetic chronology and paleoenvironmental records, also restricted the fifth sea retreat occurring from Kuqa, Kezi, and Bashibulake to~37 Ma and 33-25 Ma (Bosboom et al., 2014) ( Figure 5). In the Pamir-Tianshan convergence zone, the successions changed from marine deposits in the lower, continental clay, fine sand in the middle, and molasse in the upper part, recording the climate changes and evolution of the sea retreat . Small-scale sea incursions still occurred when the Neo-Tethys Ocean retreat occupied the Tarim Basin, which is supported by the sulfur isotope of gypsums or anhydrites in Kuqa (δ 34 S from 11.6‰ to 15.9‰ in the Middle Paleocene, from 12.7‰ to 15.9‰ in the Late Eocene) (Zhang et al., 2013). Particle size analysis of the Suweiyi Formation sediments of the Southwest Tarim Depression indicated the sand was similar to near-shore sedimentary environments including lagoons and shallow lake facies (Tang et al., 2014). The lakes in the Jidike Formation were generally broad and shallow, and the lake shoreline was in response to both climate change and tectonic activities (Ma CM. et al., 2020).
Marine and lagoon deposits with fossils were recognized from southwest to the north, which was more than a thousand of kilometers by the outcrop and wells of the Miocene formations of the northern Tarim Basin (Guo et al., 2002). There were thin gypsum interbeds of marine strata in the Late Oliogene Suweiyi Formation in the Kuqa and Late Oliogene Keziluoyi Formation reserved in Southwest Tarim (Guo et al., 2002;Ma T. et al., 2020) (Figure 5). These fossils, including Foraminifer, Dinoflagellates, Ostracods, and Calcareous nannofossils, would seem to suggest that the range of seawater influence in the Early Miocene included two main areas, Kuqa and Southwest Tarim ( Figure 5B). The high saline concentration of saltwater deposits with Foraminifer and Ostracods identified in Kuqa and Southwest Tarim showed that two lagoon bays in the western Tarim Basin were linked to the residual sea in the Late Oliogene or even until the early Miocene (Guo et al., 2002). However, from the Late Miocene to the Pliocene, the marine or lagoon environments vanished completely, and lake and freshwater deposits with Stonewort were common in the Tarim Basin (Guo et al., 2002) ( Figure 5B).
The final time of the retreat of the Neo-Tethys Ocean from the Tarim Basin remains controversial. Some believe it was in the Early Oligocene (Bosboom et al., 2011;Sun and Jiang, 2013;Sun et al., 2016b), while others believe it occurred in the Late Miocene (Guo et al., 2002). Re-evaluating the fossil evidence and sedimentary facies ( Figure 5), it was concluded that the final sea retreat occurred at 37-25 Ma, and the lagoon environments in the southwest and northwest of the Tarim Basin linked to the Neo-Tethys Ocean during the Late Oligocene. Taken together, these studies supported the notion that, as suggested by the presence of marine fossils and gypsum deposits, stepwise and diachronous sea retreat from the southwestern and northwestern margins of the Tarim Basin occurred in the late Oligocene-early Miocene.

Shortening amount of the Tarim Basin
The tectonic configuration of the Tarim Basin and its adjacent plates played an important role in the reconstruction of the proto-type basin because the compressive component of Cenozoic deformation controlled the basin (Wei et al., 2000;Li et al., 2007;Li Y. P. et al., 2020), with up to~66.7% of compressive deformation accommodated in the Cenozoic orogeny (Laborde et al., 2019). Previous studies have explored the quantified compressive component (Li et al., 2007;Laborde et al., 2019). A regional-scale balanced crosssection constructed on Southwest Tarim-Kuqa including foreland fold-thrust showed a 64% shortening rate for Southwest Tarim and 30.6% for Kuqa (Li et al., 2007). Detailed analysis and examination found that most of the Cenozoic compressive deformation (from~94% to 100%) was concentrated in the mountains along the basin margins, which meant it was negligible (up to~6%) within the Tarim Basin (Laborde et al., 2019). Therefore, it was of great importance to restore the original thickness isopaches and determine the pre-shortening tectonic configuration of the Tarim Basin.
We calculated the shortening amount of every cross-section and quantified them into four representative profiles: AA′, BB′, CC′, and DD′ (Table 2) ( Figure 1). According to the quantified data, we enlarged the range of the pre-shortening original thickness map and considered the zero-thickness line as the boundary of the preshortening proto-type basin.

Reconstruction of the proto-type Tarim Basin
Due to the uplift of the mountains around the Tarim Basin in the context of the India-Asian collision, the Tarim Basin gradually became a closed continental sedimentary basin Lin et al., 2019;Li Y. P. et al., 2020). The paleogeography of the Tarim Basin in the Cenozoic was also related to the evolution of the Neo-Tethys Ocean (Shao et al., 2006;Li et al., 2017;Zhang et al., 2018). The Paleogene strata in the Tarim Basin were bay lagoons and shore-shallow lake facies, and the Neogene strata were flood plain-salt lakes and fluvial-alluvial fan facies (Zhuang et al., 2002;Shao et al., 2006;Li et al., 2017;Yue et al., 2017;Wu et al., 2020). This study also mapped the proto-type basin in the Late Eocene to better understand the tectono-paleogeography evolution in the Cenozoic.

In the Late Eocene
The landform pattern of the Tarim Basin in the Paleocene was mainly characterized by uplift in the center with depressions in the north and south (Shao et al., 2006). A short-term invasion of seawater occurred in the Kuqa Depression and Southwest Tarim Depression from the Paleocene to the Late Eocene (Guo et al., 2002;Li et al., 2017;Zhang et al., 2018;Lin et al., 2019) (Figure 6), where there were marine sediments from tens to hundreds of meters thick deposited in fluvial-deltaic and shallow-water lacustrine settings (Shao et al., 2006).
The degree of seawater influence was decreasing in the Late Eocene due to the small-scale, stepwise Neo-Tethys Ocean retreat (Shao et al., 2006). When seawater came from the west, sabkha or the Frontiers in Earth Science frontiersin.org Note: Shortening in brackets are speculated.

FIGURE 6
The Tarim proto-type basin and its periphery in the Late Eocene.
Frontiers in Earth Science frontiersin.org 08 open platform type of depositional environment was the most important sedimentary facies (Yue et al., 2017). Turning to the period of sea retreat, the basin fill was mainly composed of several typical facies, representing: bay lagoon facies and tidal facies dominated by carbonates and gypsum in the southwestern part; shore-shallow lacustrine facies dominated by sandstone, mudstone, carbonates, and gypsum in the northern part; and fluvial fan facies and fan delta facies dominated by conglomerate, sandstone, and siltstone along the basin margin (Shao et al., 2006;Li et al., 2013) ( Figure 6).
The deformation and northward indentation of the Pamir Salient could have played an important role in controlling the lithofacies and sedimentary facies in the Tarim Basin. The provenance of the Cretaceous-Oligocene strata of the Southwest Tarim Depression was mainly from the Central Pamir (Sun et al., 2016a;Wang et al., 2021) (Figure 6), and the deposition center was near the West Kunlun Mountains (Shao et al., 2006). Regarding the Kuqa Depression, the provenance was mainly from the South Tianshan Mountains (He et al., 2004;Shao et al., 2006) (Figure 6), with the deposition center moving southward gradually (He et al., 2004).

In the Late Paleogene
The India-Asian collision had a remote effect on structural deformation and tectonic activity (Li et al., 2007;Jia, 2009;Laborde et al., 2019). For example, the West Kunlun reactivation was initiated in the Eocene with subsequent Oligocene to Miocene accelerations (Figures 2A, B), the Southwestern Tian Shan thrust wedge was formed above a decollement level (Figures 2A-C), and the Altyn Tagh strike-slip wedge was formed above a deep decollement level (Figures 2C, D).
In order to reconstruct the proto-type basin map of the Tarim Basin and its adjacent plates, we considered the structures activated during the Cenozoic and quantified the compressive component of this deformation. To the west, about 12.5 km of crustal shortening were accommodated in the Kuqa Depression, while to the east, about 20 km of crustal shortening were accommodated in the Southwest Tarim Depression (Table 2). Then we enlarged the range of the pre-shortening proto-type basin and combined the overall tectonic background. Separating the South Tianshan Mountains and the Pamir Salient was the intervening Alai Valley. The The Tarim proto-type basin and its periphery in the Late Oligocene. Thickness map contours show the depth (m) of the upper Eocene to Oligocene in the Tarim Basin.

Frontiers in Earth Science
frontiersin.org history of the Pamir Salient has been argued but data from new studies in the geology of Northern Pamir suggest that subsequent Cenozoic northward translation of the Pamir was interpreted to be only~50 km and no more than~80-100 km (Li Y. P.et al., 2020), establishing a greater degree of accuracy on the compressive component we quantified.
During the Late Eocene to Oligocene, the western part of the Tarim Basin contained an open marine, lagoon environment, with a continental environment in the eastern part   (Figure 7). Eight sedimentary facies were divided in the study area, including lagoon facies, lacustrine facies, plain braided channel facies, braided delta facies, braided delta plain facies, alluvial plain facies, and alluvial fan facies, showing that the characteristics of spatial distribution in this area were from marine to continental (Figure 7). The Southwest Tarim Depression and Kuqa Depression could be identified as major deposition centers, where the maximum sediment deposition were both over eight-hundred-meters thick (He et al., 2004;Shao et al., 2006;Zhang et al., 2022). In Kuqa, the Suweiyi Formation was mainly composed of mudstones, siltstones, sandstones, and gypsum interbeds with nannofossils ( Figure 5). In Southwest Tarim, mudstones with interbedded siltstones and laminated gypsum with nannofossils occurred in the Bashibulake Formation ( Figure 5), which is evidence of the depositional environment of a lagoon.
However, whether the lagoon facies developed its range in the Late Oligocene is controversial. Based on the latest lithofacies-paleographic map of the Tarim Basin in the Suweiyi Formation provided by Tarim Oilfield, and combined with the analysis of the Neo-Tethys Ocean retreat ( Figure 5), the outside seawater was linked to the western Tarim Basin via the Alai Valley, and the range of seawater influence in the Late Oligocene included two main areas, north to Kuqa and south to Southwest Tarim (Figure 7). In Kuqa, gypsum interbeds of marine strata and marine fossils developed in the Suweiyi Formation and Jidike Formation, which provides some support that lagoon facies developed in the Late Oligocene or even early Miocene. It could be concluded that the range of lagoon facies in the Late Oligocene was larger than the Late Eocene (Figures 6, 7). In the eastern part of the Tarim Basin, a continental sedimentary system existed from the basin margins to the center, from alluvial fan to braided river delta plain and braided river delta front (Shao et al., 2006;Zhang et al., 2022) (Figure 7).

FIGURE 8
The Tarim proto-type basin and its periphery in the Late Neogene.

Frontiers in Earth Science
frontiersin.org

In the Late Neogene
During the Neogene, the lasting tectonic compression caused the mountains surrounding the Tarim Basin to uplift and the basin closed, and the compressive component of the Neogene deformation could also be quantified. To the west, about 11 km of crustal shortening were accommodated in the Kuqa Depression, while to the east, about 10 km of crustal shortening were accommodated in the Southwest Tarim Depression (Table 2). Then we restored the range of the pre-shortening proto-type basin and combined the overall tectonic background. Compared with the Late Paleogene, the mountains, such as the Tianshan Mountains and West Kunlun Mountains, were reactivated by the ongoing structure deformation and tectonic activity (Figures 7, 8). A northward indentation of the Pamir Salient in the Cenozoic, causing a lower amount of shortening documented along the Alai Valley, and the basins, such as the Tarim Basin and the Qaidam Basin, were considered to have undergone continuous aridification.
Compared with the proto-type basin map of the Tarim Basin and its adjacent plates in the Late Eocene and Late Oligocene, the depositional environment of the Tarim Basin significantly changed in the Oligocene-Miocene: marine facies disappeared, and continental sedimentation patterns formed (Figures 6-8). The sea retreated from the whole Tarim Basin in the Late Neocene, and the Tarim Basin gradually graded to a closed inland environment (Figure 8). The Tarim Basin contained terrestrial depositional environments, such as alluvial fan groups, deltas, and flood plain facies with salt lakes (Figure 8). There were large thick molasse units near the deposition center in the Southwest Tarim Depression and Kuqa Depression, while toward the Bachu Uplift the thickness decreased gradually (Wu et al., 2020).
In the Neogene, it was still active on the south edge of the Tianshan mountain, the deformation was influencing the northern edge of the Tarim Basin (Jolivet et al., 2010). The deposition center of the Kuqa Depression gradually moved north to south and the thickness of the sediment increased (He et al., 2004). Several tectonic stages occurred in the Pamir-Tianshan convergence zone in the Cenozoic, and alluvial deposits dominated the sequences from the early Miocene   (Figure 8). The reactivation event of the Bachu Uplift could be linked to the India-Asian collision, where there were thrust fault systems may related to a deep decollement within the basement below the Tarim Block (Laborde et al., 2019) (Figures 2A, B).
Based on the lithofacies-paleographic map of the proto-type basin and our calculations of the quantified shortening (Table 2), and combined with the regional tectonic configuration, the Tarim proto-type basin and its periphery in the Late Eocene (Figure 6), Late Oligocene (Figure 7), and Late Neogene (Figure 8) was obtained. When comparing the proto-type basin and its periphery in the Late Eocene ( Figure 6) and in Late Oligocene (Figure 7), the biggest difference that could be seen was that the range of seawater influence became larger and the range of lagoon facies increased from the Southwest Tarim Depression to the Kuqa Depression (Figure 7). This variation indicated that there was still an intermittent connection between the Tarim Basin and the Neo-Tethys Ocean that occurred separately during the third, fourth, and fifth transgression-regression events in the Paleogene (Figure 5), as well as the final sea retreat which occurred at 25-30 Ma in the Late Oligocene ( Figure 5) or even the early Miocene.

Reconstruction of the tectonopaleogeography around the Tarim Basin
To obtain the tectono-paleogeographic map around the Tarim Basin in the Cenozoic, more studies were considered, such as the configuration of the plates around the Tarim Block in Asia, the lithofacies and paleogeography of basins in the adjacent plates, especially in the central Asian plates Zhang et al., 2019;Ma T. et al., 2020;Song et al., 2022), and the proto-type basin map of the Tarim Basin and its periphery in the Cenozoic we compiled (Figures 7, 8).

Continued subduction of the Indian ocean in the Late Paleogene
In the long-term paleo-geographical evolution in the Cenozoic, as the global sea levels fell, the Neo-Tethys Ocean retreat gradually occurred from the Junggar Basin and the Tarim Basin, and the Fergana Basin and the Alai Valley in southern Kyrgyzstan, and the Afghan-Tajik Basin (Bosboom et al., 2011;Carrapa et al., 2015;Li et al., 2017;Li Q. et al., 2020). In the Paleogene, the Neo-Tethys Ocean entered the western Tarim Basin via the present-day Alai Valley (Bosboom et al., 2011;Zhang et al., 2013;Bosboom et al., 2014) (Figure 9). The sea retreated westward in at least five transgressions in the west of the Tarim Basin during the Cretaceous-Cenozoic (Bosboom et al., 2011;Bosboom, 2013;Zhang et al., 2013;Bosboom et al., 2014), where tidal flat and lagoon facies associations were extended and carbonate reservoirs and mudstone covers were formed (Shao et al., 2006;Yue et al., 2017). The sedimentary lithofacies were reflected by the Neo-Tethys Ocean seawater, and developed lagoon facies in the Alai Valley and shoreshallow lacustrine facies in the east of the Tarim Basin (Guan and Guan, 2002;Guo et al., 2002;Ma CM. et al., 2020) (Figure 5).
For many decades, the time of the India-Asia continental collision has been debated. But it is accepted that the final collision was completed at~45-40 Ma, followed by continued post-collisional convergence to the present (Xia et al., 2009;Bouilhol et al., 2013;Todrani et al., 2022). The Tibetan Plateau is still experiencing deformation (Liu et al., 2009), distributing a lot of Cenozoic basin groups (Song et al., 2022) (Figure 9). The northeastern margin of the Tibetan Plateau, which includes the North Qiangtang terrane and the South Qiangtang terrane as well as the Kunlun Mountains and the Qaidam Basin, continues to deform in response to the ongoing India-Asia collision (Liu et al., 2009;Fang et al., 2020). As a huge intermontane basin, the Qaidam Basin was dominated by lacustrine facies during the Oligocene (Wu et al., 2020) (Figure 9).
The tectonic evolution of the Pamir Salient in the Cenozoic has played an important role in controlling the shift from marine to continental sedimentation in the Tarim and Tajik Basins (Bosboom Frontiers Sun and Jiang, 2013;Carrapa et al., 2015;Sun et al., 2016b). There were also basin groups distributed in Central Asia that experienced Paleozoic-Cenozoic sedimentary evolution including, for example, the Fergana Basin, Afghanistan-Tajik Basin, and Turgay Basin in the west, and the Junggar Basin and Santanghu Basin in the east (Figure 9) Ma T. et al., 2020). Several intermontane basins, such as the Yili basin and the Turpan basin, were preserved within the interior of the TianShan Mountains ( Figure 9) (Jolivet et al., 2010). In the west of the tectonopaleogeography map around the Tarim Basin in the Late Paleogene, a remnant Neo-Tethys Ocean branched via Alai Valley into a remnant marine embayment, causing the lagoon facies to develop in the Southwest Tarim Depression and Kuqa Depression in the northern Tarim Basin (Figure 9).

Basin-mountain system developed in the Late Neogene
The continental collision between the Arabian and Eurasian plates following the subduction of the Neo-Tethys Ocean beneath Eurasia resulted in a complex deformation across the Zagros Mountains ( Figure 10). There was permanent closure of the northwestern segment of the Tethyan Seaway by 12.8 Ma in the Early-Middle Miocene (Homke et al., 2004;Pirouz et al., 2015;Sun et al., 2021). Both the complex mountain building deformation across the Zagros Mountains and the further northward indentation of the Pamir Salient resulted in the closure of the Alai Valley and a complete sea retreat from the Tarim Basin ( Figure 10). The sea retreat was the one of the most significant events during the Late Oligocene-Early Miocene, causing the biggest differences in the southwest of the Tarim Basin (Figures 9, 10).
As the India-Asia continental collision is still ongoing, the Pamir Salient is still experiencing deformation. Intraplate orogeny occurred in the Kunlun Mountains and Tianshan Mountains due to the strong remote effect in the Miocene. In the mountains surrounding the Tarim Basin, mainly at the front of the Western Kunlun Mountains and Southwestern Tianshan Mountains, a significant fraction of the deformation was accommodated by new-formed structures during the Cenozoic (Li et al., 2007;Laborde et al., 2019). Accordingly, the margins of the Tarim Basin, mainly in the Southwest Tarim Depression and the Kuqa Depression, corresponded to the basin-mountain system and became rejuvenated foreland basins due to the intense regional Tectono-paleogeography map around the Tarim Basin in the late paleogene.

Frontiers in Earth Science
frontiersin.org

FIGURE 10
Tectono-paleogeography map around the Tarim Basin in the Late Neogene.

FIGURE 11
Tectonic evolution of the Tarim Basin from the Late Cretaceous to Cenozoic and its response to adjacent blocks and ocean evolution.
The present-day topography of the Tarim Basin and its surrounding mountains and adjacent blocks is considered to be related to the ongoing India-Asia collision (Figures 9-11). In the Cenozoic, the Western Kunlun Mountains and the South Tianshan Mountains were rejuvenated, and the Southwest Tarim Depression and Kuqa Depression formed rejuvenated foreland basins (Figure 11). Associated with these compressive structures, the Bachu Uplift also corresponded to the rejuvenated Western Kunlun Mountains, forming steep bounding thrusts, in particular deep fault systems (Zheng et al., 2014;Laborde et al., 2019) (Figures 2A, B).

The Tarim Block in global plate tectonics and its evolution history in the Cenozoic
The Tarim Block occupies an important position in Central Asia and the tectonic configuration around the Tarim Block has been evolving in the Cenozoic. Based on studies by previously published papers about the Tarim Block in global plate tectonics and the regional geological background in the Cenozoic (Scotese, 2004;Ma CM. et al., 2020;Wang, et al., 2020;Song et al., 2022), we combined our newly compiled proto-type basin maps (Figures 6-8) and tectono-paleogeography maps (Figures 9, 10) to modify the

FIGURE 12
Global plate reconstruction in the late Paleogene (A) and Late Neogene (B) (modified after Scotese, 2004).

Frontiers in Earth Science
frontiersin.org global plate reconstruction in the Late Paleogene ( Figure 12A) and Late Neogene ( Figure 12B). In our study, the Tarim Block was considered important and its response to adjacent blocks and ocean evolution were highlighted.
To the south of the Tarim Block, some terranes and blocks continue to deform, including the Tianshuihai, North Qiangtang, South Qiangtang, and Lhasa blocks, in response to the ongoing India-Asia collision (Figure 12). Since the collisions of the Lhasa Block and the India Plate with the Asian Plate in the Late Cretaceous, the mountains surrounding the Tarim Basin have rapidly uplifted, with a large amount of compressive deformation at the basin edges (Li et al., 2007;Jia, 2009;Laborde et al., 2019). To the north of the Tarim Block, the North and the South Tianshan Mountains were reactivated, and fold-and-thrust belts localized the topographic front of the mountain range (Jolivet et al., 2010;Jolivet et al., 2013;Laborde et al., 2019). The basins that formed on the Central Asian Orogenic Belt were small-sized continental basin groups (Song et al., 2022). To the Northwest of the Tarim Block, there were also some basin groups formed on the Kazakhstan Block, including the Syr Darya, Turgay, Chu-Sarysu, and Balkhash Basin Ma T. et al., 2020).
Besides the long-lasting plate collision, global climate change was another key on the paleogeography and tectonic evolution of the Tarim Block and its adjacent blocks. As the global sea levels fell, the Neo-Tethys Ocean gradually retreated from Asian Plate in the Cenozoic (Sun and Jiang, 2013;Carrapa et al., 2015;Sun et al., 2016a;Sun et al., 2021). The Tarim Basin has a paleogeographical evolution characterized by a long-term stepwise Neo-Tethys sea retreat punctuated by short-term shallow-marine incursions (Bosboom et al., 2011;Bosboom, 2013;Bosboom et al., 2014). The last incursion extended to the southwest and the north of the Tarim Basin with its marine deposits and fossil assemblages, and the sea retreated at~37-25 Ma ( Figure 12).

Conclusion
1. Several localized tectonic activities around the Tarim Basin occurred due to the ongoing India-Asia collision in the Cenozoic, causing tectono-paleogeography changes in the basin.
2. As the Neo-Tethys Ocean retreat finally occurred in the Tarim Basin during the Late Oligocene -Early Miocene, the proto-type basin of the Tarim Basin evolved from lagoon facies in the west and terrestrial facies in the east in the Paleogene to the evolution of terrestrial depositional environments in the Neogene. 3. By quantifying the compressive component of the Cenozoic deformation, it could be seen that the main part of the deformation was accommodated by the reactivated basinmountain system, mainly in the edges of the Tarim Basin and surrounding ranges.

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.

Author contributions
LW did the main work and wrote the paper; SH, CL and YD provided the data; JX, ZZ, XL, and HC contributed to the mapcompiling. All authors read and approved the final manuscript.

Conflict of interest
SH, CL, and YD were employed by Tarim Oilfield Company. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.