First U-Pb Detrital Zircon Ages from Kamlial Formation (Kashmir, Pakistan): Tectonic Implications for Himalayan Exhumation

This study reports the first-ever detrital zircon provenance investigation of sandstones of the Kamlial Formation, exposed in Kashmir Basin along the Kohala–Bagh road section (Muzaffarabad, Pakistan). Analysis of probability density plots of detrital U-Pb zircon ages displayed a major age population clustered around ≈400–1200 Ma and a minor age population clustered around ≈1600–1900 Ma. In addition, scattered ages existed between ≈2000 and 3000 Ma. This age pattern resembled strongly the Himalayan sources, including the Tethys Himalaya, Greater/Higher Himalaya, and Lesser Himalaya. The younger ages (<150 Ma) present in the studied samples indicated the Asian provenance. The Lesser Himalayan component (≈166–1900 Ma) was more pronounced in the 2015KM03 and 2015KM04 samples, representing the middle to the upper portion of the formation. The recycled orogen provenance of the Kamlial Formation as deduced from the sandstone petrography supports the mixed detrital zircon provenance. Considering the provenance, we propose a tectonic model that suggests that large-scale exhumation occurred in the Himalaya as a result of Panjal thrust activation during 25–14 Ma (age of Kamlial Formation), which uplifted the hinterland zone that acted as a source area that fed the foreland basin, where the Kamlial Formation deposited.


Introduction
The Himalayan fold-thrust belt and foreland basin are manifestations of the ≈50-55 Ma continental collision between India and Eurasia ( Figure 1A) [1]. The foreland basin sediments contain important information about paleo-tectonic processes that occurred during orogen development. An in-depth examination of sedimentary rocks of the basin enables us to find out the source of the sediments and rebuild their tectonic history [2]. Sedimentary rocks in the basin hold information about fluctuations in sea level, origins of sedimentary rocks, and the tectonic history of the basin [3]. The sedimentary archives of the foreland basins are extensively analyzed using modern techniques such as U-Pb dating [4][5][6][7][8] to estimate the maximum depositional ages and constrain the sediment provenance. This is a powerful tool in identifying the precise information about the sediment provenance, which is very important in explaining the paleogeographic reconstructions of the terranes in geological history [9][10][11]. Furthermore, this technique has also been widely applied in constraining the timing of orogenic events such as Ophiolite obduction and final collision in  Simplified geological map of the study area (adopted after Razzaq [27]) showing the location of the studied stratigraphic section.

Regional Geology
The Himalaya appeared on the earth's crust in response to the continent-continent collision of Indian and Asian plates in the Cenozoic [1]. The interaction of these plates occurred in response to the closure of the intervening ocean basin along the Indus suture zone/Main Mantle Thrust (MMT). As a result, a foreland basin was developed on the northern Indian margin in the footwall of MMT [28]. This basin received the detritus from the Asian plate and the Indian plate in response to the evolution of the mountain system. The Kohistan-Ladakh arc (KLA) is a magmatic arc representing the southern margin of the Asian plate, which is accreted to the Karakoram block along MKT. The Early Cretaceous to Paleocene igneous rocks belonging to KLA represent a typical island arc sequence [29]. The northern and western part of the KLA consist of Late Cretaceous-Paleocene volcanic and sedimentary rocks [30]. The Early Cretaceous (110-85 Ma) age Kohistan batholith is cropped out in the central part of the arc [31]. The Kohistan batholith is intruded by younger ≈34-29 Ma leucogranitic sills and dikes [32]. The southern part of the arc is comprised of ultramafic-mafic plutonic complexes named as Jijal complex, Spat complex, and Chilas complex of Middle Cretaceous (120-80 Ma) age [33][34][35], which intruded into a pre-early-arc metavolcanic basement complex called Kamila amphibolites of Early Cretaceous age [29,36]. The southern part of the KLA is obducted onto the northern Indian margin along the MMT and marks the collision zone between the Indian Plate and Kohistan-Karakoram terrane [37].
Immediate to the south of the collision zone (MMT), Indian Plate rocks can be described by dividing them into three tectonostratigraphic zones, which are from north to south as (1) internal metamorphosed zone, (2) external low-grade metamorphic zone, and (3) foreland basin sediments [38]. The internal metamorphosed zone comprises an intensely deformed imbricate stack of middle Proterozoic to Archean metamorphic rocks, which are dominantly high-grade gneisses and pelitic schists intruded by Cambrian-Ordovician granitic intrusions [39]. The internal zone is dominantly characterized by Barrovian-type high-grade metamorphism, but locally, high-pressure metamorphism is represented by eclogite such as those reported from the upper Kaghan valley [40]. The Panjal Thrust separates the internal metamorphosed zone from the external low-grade metamorphic zone. This external zone consists of weakly metamorphosed late Proterozoic rocks and Mesozoic to Eocene sedimentary deposits.
Further to the south, the MBT marks the boundary between the external zone and the foreland basin sequence. It consists of Paleocene-Eocene marine sediments, overlain by Oligocene-Miocene and younger siliciclastic strata marking regional unconformity [41,42]. In terms of the Himalayan classification [43], part of the internal metamorphosed zone, north of MCT (Laut Thrust/Batal Fault), marks the Greater/Higher Himalaya ( Figure 1A,B). The Lesser Himalayan zone is marked between the MBT and MCT (Laut Thrust/Batal Fault). The Hazara-Kashmir syntaxial region in Pakistan represents the foreland basin's northern extreme, currently occupying the sub-Himalayan region [8]. It contains marine sediments (Paleocene-Eocene) overlain by molasse sediments (Oligocene-younger) [44]. The Kamlial Formation lies toward the top of the stratigraphy of the molasse sequence in the Sub Himalaya [7] and is positioned stratigraphically below the Siwaliks, around 35 km to the southeast of the Muzaffarabad ( Figure 1B).
The Kamlial Formation is widespread in the sub-Himalaya in Pakistan with variable thickness, ranging from 100 to ≈580 m [25]. The Kohala section represents its northernmost exposure within the sub-Himalaya. The thickness of the Kamlial Formation in this section is ≈102 m (Figure 2A) and exposed in a synclinal structure ( Figures 1C and 2B). The dominant component of the formation is greenish-gray sandstone, which is ≈74 % and portrays spheroidal weathering. The sandstone is medium to coarse-grained and well compacted. It also consists of wood ( Figure 2C) and leaf imprint fossils. The second major portion of the formation consisted of maroon siltstone and mudstone ( Figure 2D), which is ≈20 % of the total thickness. The rest of the formation consisted of intraformational conglomerates and paleosol layers ( Figure 2E,F). The bioturbation is also observed in the sandstone ( Figure 2G). In addition, volcanic clast ( Figure 2H) is also observed in the sandstone, possibly associated with the Panjal volcanics. It is also observed that the lower portion of the Kamlial Formation is dominated in siltstone and mudstone, while the upper part is dominated in sandstone ( Figure 2B). The Kamlial Formation consisted of a cyclic deposition of mudstone, siltstone, and sandstone in a fining upward sequence. This cyclic sequence and sedimentary structures represent the deposition by braided streams. ( Figure 2G). In addition, volcanic clast ( Figure 2H) is also observed in the sandstone, possibly associated with the Panjal volcanics. It is also observed that the lower portion of the Kamlial Formation is dominated in siltstone and mudstone, while the upper part is dominated in sandstone ( Figure 2B). The Kamlial Formation consisted of a cyclic deposition of mudstone, siltstone, and sandstone in a fining upward sequence. This cyclic sequence and sedimentary structures represent the deposition by braided streams. The Kamlial Formation has transitional contact with the underlying Murree Formation ( Figure 2I), while the upper contact is not exposed in the study area. However, in the vicinity of the studied area, it is conformable with the overlying Chinji Formation of the Siwalik group.

Petrography
The traditional petrography is widely applied to determine the composition of the sandstones and their provenance [45,46]. The sandstone petrography was used to supplement the detrital zircon U-Pb age data in this study. A total of five sandstone samples from the entire section were selected for sandstone petrography. Thin sections of the samples were prepared in the Rock Cutting and Thin section lab at the Department of Earth Sciences, COMSATS University Islamabad, Abbottabad Campus, Pakistan.
Each thin section was examined under a polarized petrographic microscope, and >300 framework grains (Quartz, feldspars and lithics) were counted from the different views of the thin section following the Gazzi-Dickinson point-counting method [46].

Zircon Imaging
The detrital zircons extracted from the selected studied samples represent the multiple sources from the different blocks, consisting of sedimentary, metamorphic, and igneous rocks. It is important to image the detrital zircons before the in situ U-Pb analysis [47,48]. The cathode luminescence (CL) images were recorded for all three samples, which clearly show the internal structure of the detrital zircons (Supplementary Figure S1). The analysis spots are placed at the outer part in zircons, which reflects a complex rim to core internal structure (Supplementary Figure S1).

U-Pb Detrital Zircon Dating
U-Pb dating of the detrital zircon is a modern technique extensively used in clastic sediments and provides significant information about the sediment source region by comparing ages with the ages of the adjacent source terranes. In this study, the samples of the Kamlial Formation were selected for detrital zircon U-Pb dating to constrain the provenance. A total of three samples were collected representing the bottom, mid, and upper part of the formation from the southwestern limb of Chamankot syncline ( Figures 1C and 2A).
The detrital zircon grains were isolated using heavy liquids and the magnetic separation method from the crushed sandstone samples. In a second stage, these isolated zircon grains were mounted on a glue strip, and then, the epoxy resin was poured on the strip to prepare the sample for polishing. Finally, the prepared samples were polished to make the surface of the grains smooth.
In the analysis stage, to avoid lead contamination, the mounted detrital grains' surface was cleaned with dilute nitric acid and pure alcohol before U-Pb analysis. The U-Pb age dating was conducted at the Key Laboratory of Continental Collision and Plateau Uplift, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing, China. From each sample, 100 individual detrital grains were analyzed using an Agilent 7500a Laser Ablation Inductively Plasma Mass Spectrometer (LAICPMS). Each analyzed spot has a diameter of 35 µm. The ages of the unknown detrital grains were calibrated using the Plesovice standard (average age of 337 + 0.37 Ma) [49]. The final ages of the detrital grains were selected using two age systems. The 207 Pb/ 206 Pb age system was adopted for the detrital zircon grains with ages > 1400 Ma, whereas the 206 Pb/ 238 U age system was adopted for the zircon grain with ages < 1400 Ma [50]. The grains with >10% discordance in each age system were excluded in the final data interpretation. The raw data were processed using Glitter 4.0 software. The age data were displayed in probability density plots using Isoplot [51] and provided in Table S1.

Petrography
Five thin sections were prepared from the collected samples covering the entire measured section. The samples 2015KM01 and 2015KM02 represented the lower part and were collected from the medium to coarse-grained greenish color thickly bedded sandstone ( Figure 2A). The samples 2015KM03 and 2015KM05 represented the middle part and were collected from the greenish-gray sandstone. Similarly, the sample 2015KM04 represented the upper portion of the Kamlial Formation and was collected from the very coarse-grained greenish sandstone (Figure 2A).
The grains of samples are sub-rounded to angular, with calcite and siltstone cementing materials. The framework grains consisted of quartz, feldspar, and lithics ( Figure 3A). The monocrystalline quartz is abundant ( Figure 3B,G). However, few polycrystalline quartz grains were observed in all thin sections ( Figure 3C). The feldspar includes both the Kfeldspar and Na-feldspar. The albite is prominent in the thin section ( Figure 3B,G). The lithics consisted of all three rocks ( Figure 3B-E). Igneous and metamorphic lithics are very common (Figure 3B,E,G). The sedimentary lithics were mostly the shales and siltstone ( Figure 3C,D). However, few sedimentary carbonate lithics were also observed in the thin section ( Figure 3F,G). In addition, mica was also observed in the thin sections, including biotite and muscovite ( Figure 3D,G). the detrital zircon grains with ages >1400 Ma, whereas the 206 Pb/ 238 U age system was adopted for the zircon grain with ages <1400 Ma [50]. The grains with >10% discordance in each age system were excluded in the final data interpretation. The raw data were processed using Glitter 4.0 software. The age data were displayed in probability density plots using Isoplot [51] and provided in Table S1.

Petrography
Five thin sections were prepared from the collected samples covering the entire measured section. The samples 2015KM01 and 2015KM02 represented the lower part and were collected from the medium to coarse-grained greenish color thickly bedded sandstone (Figure 2A). The samples 2015KM03 and 2015KM05 represented the middle part and were collected from the greenish-gray sandstone. Similarly, the sample 2015KM04 represented the upper portion of the Kamlial Formation and was collected from the very coarsegrained greenish sandstone (Figure 2A).
The grains of samples are sub-rounded to angular, with calcite and siltstone cementing materials. The framework grains consisted of quartz, feldspar, and lithics ( Figure 3A). The monocrystalline quartz is abundant ( Figure 3B,G). However, few polycrystalline quartz grains were observed in all thin sections ( Figure 3C). The feldspar includes both the K-feldspar and Na-feldspar. The albite is prominent in the thin section ( Figure 3B,G). The lithics consisted of all three rocks ( Figure 3B-E). Igneous and metamorphic lithics are very common (Figure 3B,E,G). The sedimentary lithics were mostly the shales and siltstone ( Figure 3C,D). However, few sedimentary carbonate lithics were also observed in the thin section ( Figure 3F,G). In addition, mica was also observed in the thin sections, including biotite and muscovite ( Figure 3D,G).  The percentage composition of the framework grains calculated through the pointcounting method suggests that the quartz concentration ranges between ≈52 and 75%, the feldspar concentration ranges between 4 and 6.92%, and the lithics concentration ranges between 20 and 41.42% ( Figure 3A). The data plotted on the QtFLt diagram of the Dickinson [46] placed the samples of the Kamlial Formation within the recycled orogen provenance ( Figure 3H).

Detrital Zircon U-Pb Dating Results
Detrital zircon geochronology has been the most accurate and reliable method for determining the provenance of sediments for the last few decades because detrital zircon is chemically stable and mechanically robust in various depositional environments and weathering conditions. According to previous provenance studies based on U-Pb detrital zircon dating, the recommended number of selected zircon grains is 80-100 grains per sample to precisely determine the source components of the sediments [52]. In this research, we adopted U-Pb age dating of the detrital zircons, extracted from sandstones of the Kamlial Formation exposed along Kohala-Bagh road section, sub-Himalaya, Azad Kashmir, Pakistan, for constraining the provenance. The results of the analyzed samples are described below in detail.
Sample 2015KM01 was collected above the lower contact of the Kamlial Formation with Murree Formation from a coarse sandstone horizon (Figure 2A). One hundred detrital zircons were analyzed for U-Pb age dating. One hundred analyses yielded 98 concordant ages ( Figure 4, Table S1). The highest concentration of detrital zircons occurred between ≈400 and ≈1000 Ma, with peak ages of ≈432 Ma, ≈484 Ma, and ≈875Ma ( Figure 4). There were a few subordinate peaks of Paleoproterozoic ages (≈1600-2500 Ma) and Mesoproterozoic ages (≈1000-1600 Ma) age. In addition, a younger age cluster existed between ≈31 and ≈134 Ma, with peak ages at ≈72 Ma and ≈93 Ma (Figure 4).
The sample 2015KM03 represents the middle part of the formation. All the 100 detrital zircons yielded concordant ages. The major proportion of the detrital zircon ages existed between ≈453 and ≈1150 Ma, which is ≈70% of the total proportion ( Figure 4). The age peaks were recorded in this broad range at ≈495 Ma, ≈525 Ma, ≈615 Ma, and ≈978 Ma (Figure 4). The second cluster existed between ≈1500 and ≈2000 Ma, which is ≈10% of the total proportion. The rest of the detrital ages were scattered between ≈2100 and ≈3300 Ma (Figure 4). Only two younger detrital ages were reported from the 2015KM03 sample, which were 95 Ma and 118 Ma.

Internal Zoning Pattern and Th/U ratio
The CL images of the detrital zircons show different internal zoning patterns [47], which indicates the derivation from different sources, including igneous, metamorphic, and sedimentary rocks [54,55]. The size of the zircon grains ranges from 50 to >200 μm (Supplementary Figure S1). The analyzed zircons abundantly displayed the oscillatory zoning pattern, indicating the igneous originʹs derivation. In addition, the zircons with xenocrystic cores are also present [47]. A fair number of the detrital zircons displayed the metamorphic textures with a well-developed core and younger rims ( Supplementary Figure S1). However, a sufficient number of zircons also display plane texture without any zonation. Some analyzed zircons possessed sector zoning patterns (Supplementary Figure  S1).
In addition, the zircon Th/U ratio is also an important elemental ratio, which is used to discriminate between the metamorphic and igneous origin zircon [56,57]. The detrital zircon grains derived from the igneous origin possessed a Th/U ratio value >0.3, whereas

Internal Zoning Pattern and Th/U ratio
The CL images of the detrital zircons show different internal zoning patterns [47], which indicates the derivation from different sources, including igneous, metamorphic, and sedimentary rocks [54,55]. The size of the zircon grains ranges from 50 to >200 µm (Supplementary Figure S1). The analyzed zircons abundantly displayed the oscillatory zoning pattern, indicating the igneous origin's derivation. In addition, the zircons with xenocrystic cores are also present [47]. A fair number of the detrital zircons displayed the metamorphic textures with a well-developed core and younger rims (Supplementary Figure S1). However, a sufficient number of zircons also display plane texture without any zonation. Some analyzed zircons possessed sector zoning patterns (Supplementary Figure S1).
In addition, the zircon Th/U ratio is also an important elemental ratio, which is used to discriminate between the metamorphic and igneous origin zircon [56,57]. The detrital zircon grains derived from the igneous origin possessed a Th/U ratio value > 0.3, whereas the detrital zircons with a Th/U ratio <0.3 indicate the metamorphic origin [57]. The U-Pb ages of the detrital zircons and Th/U ratio plot reflect the detrital zircons' dominant population from the igneous origin ( Figure 5). However, a significant population between ≈400 and ≈1200 Ma reflects the metamorphic origin. This metamorphic episode is most possibly associated with the Pan-African orogeny, which is well documented in the northwestern Himalayas [58]. At the same time, the igneous zircons in this similar age domain might represent the lesser Himalayan granites, which possessed similar ages [59,60]. The younger population of detrital ages < 100 Ma reflects the igneous origin dominantly ( Figure 5). In contrast, the younger grains reflect the metamorphic origin, which is more likely associated with Himalayan metamorphism [61,62]. northwestern Himalayas [58]. At the same time, the igneous zircons in this similar ag domain might represent the lesser Himalayan granites, which possessed similar age [59,60]. The younger population of detrital ages <100 Ma reflects the igneous origin dom inantly ( Figure 5). In contrast, the younger grains reflect the metamorphic origin, which is more likely associated with Himalayan metamorphism [61,62].

U-Pb Ages of Source Terranes
The interpretation of the detrital zircon provenance for the middle to late Miocen Kamlial Formation, exposed along the Kohala-Bagh road section, requires the representa tive U-Pb age dataset for the adjacent litho-tectonic terranes, which were possible contrib uting sediment sources during Himalayan mountain evolution. These litho-tectonic ter ranes include Himalayan tectonostratigraphic zones (Tethyan Himalaya, Higher/Greate Himalaya (GH), and Lesser Himalaya (LH)) and Asian sources (Kohistan-Ladakh ar (KLA), Karakoram Block (KB), and Lahasa Block (LB)) ( Figure 4) [53]. These terranes hav distinct age spectra, which discriminate them from one another. The Asian sources wer grouped for simplicity, as these terranes were already accreted to the Indian plate follow ing the final collision [7].
The Asian terranes possessed a distinct age pattern, with a majority of the zircon age clustered between ≈20 and ≈120 Ma (Figure 4). These ages were representative of th younger batholiths associated with the magmatism during the collision process. The pres ence of these younger ages in the evolving basins indicates these terranesʹ contribution The Himalayan tectonostratigraphic zones are important in understanding the evolution of the mountain belt. These tectonostratigraphic zones possessed a unique age pattern representing the specific blockʹs sources. The zircons from the LH typically yielded ag clusters around 1700-1900 Ma with subordinate minor age clusters around 2400-2600 M The interpretation of the detrital zircon provenance for the middle to late Miocene Kamlial Formation, exposed along the Kohala-Bagh road section, requires the representative U-Pb age dataset for the adjacent litho-tectonic terranes, which were possible contributing sediment sources during Himalayan mountain evolution. These litho-tectonic terranes include Himalayan tectonostratigraphic zones (Tethyan Himalaya, Higher/Greater Himalaya (GH), and Lesser Himalaya (LH)) and Asian sources (Kohistan-Ladakh arc (KLA), Karakoram Block (KB), and Lahasa Block (LB)) ( Figure 4) [53]. These terranes have distinct age spectra, which discriminate them from one another. The Asian sources were grouped for simplicity, as these terranes were already accreted to the Indian plate following the final collision [7].
The Asian terranes possessed a distinct age pattern, with a majority of the zircon ages clustered between ≈20 and ≈120 Ma (Figure 4). These ages were representative of the younger batholiths associated with the magmatism during the collision process. The presence of these younger ages in the evolving basins indicates these terranes' contribution. The Himalayan tectonostratigraphic zones are important in understanding the evolution of the mountain belt. These tectonostratigraphic zones possessed a unique age pattern representing the specific block's sources. The zircons from the LH typically yielded age clus-ters around 1700-1900 Ma with subordinate minor age clusters around 2400-2600 Ma [53]. The detrital zircons from the GH sequence displayed age clusters around 540-750 Ma, 800-1200 Ma, 1600-1900 Ma, and 2400-2600 Ma (Figure 4). The dominant ages of the GH sequence clustered around 900-1100 Ma [28,53,63]. The detrital zircon record of the intruding granitic body depicted age clusters between 470 and 550 Ma with a prominent peak at 485 Ma [64]. The TH detrital age spectra are similar to the GH sequence. However, the 480-570 Ma age cluster is more distinct than other clusters such as 700-1200 Ma and 2430-2560 Ma (Figure 4) [65,66]. The Tethyan strata also contain younger detrital zircon grains ranging between 110 and 140 Ma, which are representative of Indian Plate volcanic rocks [15,16,28,53].

Provenance of Kamlial Formation
The detrital zircon U-Pb plots of the studied samples of Kamlial Formation portrayed the age pattern mainly clustered around ≈400-1200 Ma with major peaks at ≈550 Ma, ≈650 Ma, ≈800 Ma, and ≈978 Ma (Figure 3). In addition to this prominent component, a significant age pattern existed between ≈1600 and ≈1900 Ma. This age component is more distinct in the samples representing middle and upper parts ( Figure 4). In comparison, minor scattered ages existed between ≈2000 and ≈3000 Ma. The younger detrital zircon ages reported in the samples were low in percentage (≈8%). However, these younger ages suggested the Asian provenance contribution. The Th/U ratio of the younger grains are mostly higher than 0.3, which supports the igneous origin. However, few younger grains with Th/U ratio lower than 0.3 reflect Himalayan metamorphism's imprints.
The major age cluster (≈400-1200 Ma) resembled the Tethyan and Greater Himalayan sequences ( Figure 4). It suggests that the major contributing tectonostratigraphic zones during the deposition of the Kamlial Formation were Tethyan and Greater Himalayas. However, the increased input of the detritus portraying the ages between ≈1600 and ≈1900 Ma reflected the contribution of the Lesser Himalayan tectonostratigraphic zone. This pattern is comparatively more pronounced in the samples representing the middle and upper portions of the formation. Considering the contribution from the Himalayan sources and Asian sources, it is suggested that the Kamlial Formation received the mixed detritus indicative of the sediment recycling and southward propagation of the fold-thrust belt. The sandstone petrography also supports the sediment recycling reflected by "Recycled orogen" provenance ( Figure 3H). This sediment recycling also supports the uplifting of the blocks and erosion and depositing in a new site.

Tectonic Implications for Himalayan Exhumation
The Himalayas resulted from the India-Asia collision that occurred ≈56 Ma [1,4,7]. Following the collision, a foreland basin was developed in the footwall of the MMT, which received detritus from the accreted Asian terranes, mainly the Kohistan-Ladakh arc ( Figure 6A). As a result of the continued subduction of the Indian plate in the north underneath, the Asian plate gave rise to Himalayan belt evolution by southward propagation of the fold-thrust belt.
In the first developed foreland basin, the shallow marine carbonate sequence was formed in Margalla Hill Limestone and Chorgali formations [7]. The shallow marine conditions prevailed during the foreland basin's initial development, such as the case of the Australia and Papua New Guinea collision [67]. With the emergence of Batal  In the first developed foreland basin, the shallow marine carbonate sequence w formed in Margalla Hill Limestone and Chorgali formations [7]. The shallow marine co ditions prevailed during the foreland basin's initial development, such as the case of t Australia and Papua New Guinea collision [67]. With the emergence of Batal Fault/MC Relying on the U-Pb age data and sandstone petrography of the Kamlial Formation suggests that the provenance is from the Himalayan sources, including Tethys Himalaya, Greater/Higher Himalaya, and Lesser Himalaya. These Himalayan tectonic terranes consisted of particular detrital zircon age signatures. These terranes uplifted with the development of new faults in response to fold-thrust belt propagation. Particularly, the U-Pb ages compared with the Th/U ratio indicates the derivation from the igneous and metamorphic sources of the Proterozoic and Cambro-Ordovician rocks. These Proterozoic and Cambro-Ordovician age metamorphic rocks are widely exposed in the Higher Himalayan crystalline and Lesser Himalayan sequences. The U-Pb ages together with the Th/U ratio suggest the source of the Kamlial Formation from the rocks exposed in Tethyan, Greater/Higher, and Lesser Himalaya. This mixed source suggests that the foreland basin during the middle to late Miocene time existed in the footwall of the Panjal fault that received the source from the Lesser Himalayan block ( Figure 6C). Considering the age of the Kamlial Formation, it can be deduced that the Panjal fault activated during 25-14 Ma. The felsic lithics and increased Cambrian-Ordovician ages in the detrital record suggest the exhumation of the Mansehra granite, which is mainly exposed in the hanging wall of the Panjal fault.

Conclusions
This study presented first-time detrital zircon U-Pb geochronology supported with sandstone petrography of the middle to late Miocene Kamlial Formation to constrain the provenance and its implications for Himalayan Exhumation in the western Himalayas, Pakistan. Our results suggest that the sediments of the Kamlial Formation were mainly sourced from the Himalayan tectonostratigraphic zones, including Tethyan, Greater and lesser Himalaya, with a contribution from the Asian sources as reflected by younger ages (<150 Ma). It is worth mentioning that the Lesser Himalayan contribution (≈1600-1900 Ma) is more pronounced in the middle and upper part of the formation, as reflected in samples 2015KM03 and 2015KM04. Considering this provenance, we proposed a tectonic model that explains the exhumation of the Himalayan tectonostratigraphic zones. This model suggests that the Tethyan and Greater Himalayan tectonostratigraphic zones were already uplifted during the Eocene to Early Miocene time and feeding the foreland basin in the footwall of MCT. However, during the deposition of the Kamlial Formation, the Lesser Himalayan tectonostratigraphic zone uplifted in response to the activation of the Panjal Fault at least around ≈24-14 Ma. With this exhumation, the contribution from the Lesser Himalayan block started and was reflected in the sediments of the Kamlial Formation.