Articles | Volume 12, issue 11
Solid Earth, 12, 2479–2501, 2021
Solid Earth, 12, 2479–2501, 2021

Research article 03 Nov 2021

Research article | 03 Nov 2021

Detrital zircon provenance record of the Zagros mountain building from the Neotethys obduction to the Arabia–Eurasia collision, NW Zagros fold–thrust belt, Kurdistan region of Iraq

Detrital zircon provenance record of the Zagros mountain building from the Neotethys obduction to the Arabia–Eurasia collision, NW Zagros fold–thrust belt, Kurdistan region of Iraq
Renas I. Koshnaw1,a, Fritz Schlunegger1, and Daniel F. Stockli2 Renas I. Koshnaw et al.
  • 1Institute of Geological Sciences, University of Bern, Baltzerstrasse 1+3, 3012 Bern, Switzerland
  • 2Department of Geological Sciences, Jackson School of Geosciences, The University of Texas at Austin, Austin, TX 78712, USA
  • apresent address: Department of Structural Geology and Geodynamics, Geoscience Center, University of Göttingen, Goldschmidtstraße 3, 37077 Göttingen, Germany

Correspondence: Renas I. Koshnaw (


Recognition of a new angular unconformity and the synthesis of new detrital zircon U–Pb and (U–Th) / He provenance records, including zircon (U–Th) / He double dating, from the NW Zagros region elucidate the basin dynamics of the foreland wedge-top and intermontane units, as well as the tectonic processes in the source terranes in response to the different geodynamic phases. In this contribution, we present field observations and detrital zircon provenance data from hinterland basins to reconstruct the basin dynamics and the underlying tectonic controls in the NW Zagros area in the Kurdistan region of Iraq. Results reveal that the deposition of the suture zone units of the Red Beds Series (RBS; Suwais Group, Govanda Formation, Merga Group) occurred in an intermontane basin on top of folded Upper Cretaceous units and that the RBS deposits rest with an angular unconformity on the underlying older strata. The RBS provenance data point to the Paleogene Walash–Naopurdan–Kamyaran (WNK) complex as a source area and imply a substantial decrease in magmatism by  36 Ma, as reflected by the youngest age peaks. New detrital zircon provenance data from the hinterland wedge-top units of the proto-Zagros foreland basin (the Tanjero, Kolosh, and Gercus formations) exhibit exclusive derivation from the Upper Cretaceous Neotethys ophiolitic terranes, which differs from the provenance of the older Lower Cretaceous and Paleozoic units that are dominated by the Paleozoic and Neoproterozoic age spectra. These shifts in provenance between different tectonostratigraphic units argue for a sediment source reversal from the west to the east in response to ophiolite obduction, arrival of the WNK complex, and commencement of the Arabia–Eurasia continental collision during the latest Eocene (< 36 Ma). According to the provenance data, the incipient collision was followed by the deposition of the RBS in the hinterland of the proto-Zagros fold–thrust belt as well as the connection of drainages with the collision-related Neogene foreland basin.

1 Introduction

In sedimentary basins adjacent to collisional mountain belts, synorogenic clastic deposits potentially provide information on the geodynamic setting prior to and during the collisional processes. Understanding the causes and consequences of mountain-building processes is critical for the reconstruction of paleogeography, paleoclimate, and for the understanding of the long-term deformation and erosion history (Molnar and England, 1990; Avouac and Burov, 1996; Najman et al., 2010, Norton and Schlunegger, 2011; Avouac et al., 2015). Hinterland basins, such as wedge-top and intermontane basins are valuable archives for the assessment of the exhumation and unroofing history of the adjacent uplifted terranes because of their proximity to the source areas. Nevertheless, well-preserved ancient stratigraphic successions are scarce due to the deformation of sedimentary strata as orogenesis proceeds (Horton, 2012; Orme et al., 2015). A possible approach to overcome this drawback is the utilization of the geochronologic and thermochronologic records, which are preserved by detrital zircons (e.g., Cawood et al., 2012; Webb et al., 2013; Gehrels, 2014; Colleps et al., 2020).

The Zagros fold–thrust belt is a prominent collisional orogen on Earth that stretches for  2000 km across the Middle East (Fig. 1). This orogenic belt formed during the Late Cretaceous and the Cenozoic as a consequence of the convergence and the subsequent collision of the Arabian and Eurasian continental plates (Dewey et al., 1973; Hempton, 1985; Dercourt et al., 1986; Mouthereau et al., 2012). This prolonged history of deformation resulted in an amalgamation of different tectonic terranes that were situated between the Arabian and Eurasian plates, such as the Bisotoun block, the middle-Cretaceous intraoceanic oceanic subduction and back-arc spreading zone, and the early-Tertiary magmatic domain. Because the accretion processes overprinted the preceding tectonic configurations (Wrobel-Daveau et al., 2010; Agard et al., 2011; Vergés et al., 2011; Barber et al., 2018, 2019), any paleotectonic models particularly for the NW Zagros suture zone have been associated with uncertainties. Previous paleotectonic reconstructions in the NW Zagros suture zone involved several subduction zones, a mantle plume, and different loci for specific tectonic terranes and sedimentary successions of suture zones (e.g., Numan et al., 1997; Karim et al., 2008; Al-Qayim et al., 2012; Azizi et al., 2013; Whitechurch et al., 2013; Ali et al., 2014, 2019; Moghadam et al., 2020). Similarly, possibly due to the complexity of the convergence history, various timings have been suggested for the collision (e.g., Fig. 2 from Zhang et al., 2017) with a conceivable minimum age during the late Oligocene (Koshnaw et al., 2019).

Figure 1(a) Regional tectonic map of the Middle East displaying the study area in the Kurdistan region of Iraq, and the suture of the Arabian and Eurasian plates along the Main Zagros fault (Koshnaw et al., 2017, and references therein). (b) Geologic map of the NW Zagros area in the Kurdistan region of Iraq showing new (yellow) and published (orange) sample locations (Koshnaw et al., 2019, 2020a, and references therein).

This research aims at constraining the basin dynamics recorded by the suture zone deposits and the wedge-top units, and the convergence history of the Arabian and Eurasian continental plates based on data established with the detrital zircon U–Pb and (U–Th) / He double-dating methods. The targeted NW Zagros hinterland deposits are the Red Beds Series (potential equivalent of the Razak Formation in Iran; Etemad-Saeed et al., 2020), which is situated below the Main Zagros fault, and the proto-Zagros Tanjero (Amiran), Kolosh (Amiran), and Gercus (Kashghan) formations, in the Kurdistan region of Iraq (Fig. 1b). We will combine the new provenance information with published data from the wedge-top deposits of the Neogene Zagros Basin, including the Fatha (Gachsaran), Injana (Aghajari), Mukdadiya (Lahbari Member?/Lower Bakhtiari, where “?” indicates uncertainty), and Bai-Hasan (Bakhtiari) formations, as well as the Lower Cretaceous Garagu Formation and the Paleozoic Khabour and Pirispiki formations (Fig. 2). The ultimate goal is to track the provenance shifts in the erosional hinterland, which will build the basis to invoke the occurrence of specific geodynamic events throughout the Zagros orogeny, including subduction, collision, and post-collisional uplift and exhumation.

Figure 2Generalized Phanerozoic stratigraphic column and key tectonic events for the NW Zagros fold–thrust belt and foreland basin, including the suture zone units in the Kurdistan region of Iraq (Sissakian et al., 1997; Sharland et al., 2001; English et al., 2015). The names in parentheses represent the potential equivalent formation names in Iran. Red circles represent stratigraphic units that were addressed in this research.


2 Geologic background

2.1 The pre-Zagros orogeny strata

From the Cambrian and throughout most of the Paleozoic (Fig. 2), the NE Arabian stratigraphic succession developed on the Pan-African basement and was part of the northern Gondwana realm along the Paleotethys Ocean. The Paleozoic stratigraphic column contains mostly clastic sedimentary rocks with some marginal marine carbonate successions (Konert et al., 2001). Prior to the onset of the Hercynian orogeny during the Late Devonian, the Arabian Plate was considered as a stable platform. During the Hercynian orogeny, however, the platform experienced intracratonic deformation, which lead to several zones of subsidence and uplift (Sharland et al., 2001). During the Permo-Triassic, the eastern margin of Gondwana experienced rifting and the Neotethys Ocean started to open. In the Kurdistan region of Iraq, the exposed Paleozoic rocks include two main clastic formations, which are referred to as the Khabour (Ordovician) and Pirispiki (Late Devonian) formations (Al-Hadidy, 2007). The Khabour Formation consists of thin-bedded, fine-grained sandstones, quartzites, and shales that were deposited in a shallow to deep marine environment (van Bellen et al., 1959; Al-Bassam, 2010). The Pirispiki Formation includes sandstones, siltstones, and red mudstone beds as well as conglomerates, which have been interpreted as records of a partially marine to a fluvial environment but mostly nonmarine (van Bellen et al., 1959; Al-Hadide et al., 2002; Al-Hadidy, 2007; Al-Juboury et al., 2020).

After the opening of the Neotethys Ocean during the Permo-Triassic, the border of the Arabian Plate developed into a passive margin that accumulated thick sequences of limestones, shales, and evaporite beds during most of the Mesozoic (Fig. 2) (Sharland et al., 2001; Stampfli and Borel, 2002; Ziegler, 2001; English et al., 2015). During the Middle Jurassic, the Gotnia Basin was formed, which hosted deep-water carbonates and shales (Aqrawi et al., 2010). During the Late Jurassic, the west and northwestern part of Iraq experienced uplift, resulting in the formation of the Khlesia and Mosul basement highs, and the Gotnia Basin configuration changed considerably (Ameen, 1992; Aqrawi et al., 2010; English et al., 2015). This variation in basin geometry resulted in the deposition of relatively thin (<100 m) carbonate rocks in association with clastic sandstone beds in the western side of the basin during the Early Cretaceous. Toward the eastern side of the basin, a basinal carbonate succession such as the Balambo Formation ( 2000 m) was deposited (van Bellen et al., 1959; Jassim and Buday, 2006; English et al., 2015). Later, during most of the Late Cretaceous, the basin was a site for further deposition and filled with shelfal and lagoonal carbonates. During the Campanian–Maastrichtian time, the northern Arabian Plate was affected by a rifting phase that lead to the formation of grabens filled with synrift marly deposits. In addition, an environment with shallow-water carbonates established on structural highs (Aqrawi et al., 2010; English et al., 2015).

2.2 The proto-Zagros foreland basin strata

By the Maastrichtian, the proto-Zagros flexural foreland basin started to form in response to the arrival of the Neotethys intraoceanic subduction zone at the Arabian Plate margin, leading to the ophiolite obduction (Stoneley, 1975; Koop and Stoneley, 1982; Alavi, 1994; Homke et al., 2010; Saura et al., 2011; Barber et al., 2019). This process resulted in a basin geometry where the sedimentary infill thinned toward the west with a surge of clastic input sourced mostly from the NW and the NE (Karim and Surdashy, 2005; Aqrawi et al., 2010; Çelik and Salih, 2021). The main clastic formations of the proto-Zagros foreland basin in the Kurdistan region of Iraq are the Tanjero, Kolosh, and Gercus formations with the ages of Maastrichtian, Paleocene, and Eocene, respectively (Fig. 2) (van Bellen et al., 1959; Abdel-Kireem, 1986). The reported nature of the contact between these formations is both conformable and unconformable (van Bellen et al., 1959; Jassim and Sisakian, 1978; Tamar-Agha et al., 1978; Al-Shaibani and Al-Qayim, 1990; Jassim and Buday, 2006; Abdallah and Al-Dulaimi, 2019; Kharajiany et al., 2019). The depositional environment of the proto-Zagros foreland basin strata was interpreted as distal to proximal deep-marine deposits that shallow upward to nonmarine deposits (Al-Qayim, 1994; Al-Qayim et al., 2008; Jassim and Buday, 2006; Kadem, 2006; Al-Mashaikie et al., 2014). By the late Eocene, this clastic succession of the proto-Zagros foreland basin was covered by the shallow lagoonal marine carbonates of the Pila-Spi (Shahbazan?) Formation (van Bellen et al., 1959; Aqrawi et al., 2010). On top of the Pila-Spi Formation, an unconformity has been recorded based on absence of the Oligocene–early-Miocene rocks that are linked to the collision between Arabia and Eurasia (Fig. 2) (Dunnington, 1958; Ameen, 2009; Lawa et al., 2013).

2.3 The Neogene Zagros foreland basin strata

During the Neogene, the continent–continent collision between the Arabian and Eurasian plates became prominent and, consequently, new clastic material was supplied into the new flexural basin (Koshnaw et al., 2017, 2020a). The synorogenic Neogene Zagros foreland basin formations are the Fatha (Gachsaran), Injana (Aghajari), Mukdadiya (Lahbari Member?/Lower Bakhtiari?), and Bai-Hasan (Bakhtiari) formations (Fig. 2). From bottom to top, the stratigraphic succession changes from the mixed clastic–carbonate–evaporite deposits of the Fatha Formation to the meandering and braided fluvial, and alluvial deposits of the Injana, Mukdadiya, and Bai-Hasan formations (Shawkat and Tucker, 1978; Tamar-Agha and Al-Aslami, 2015; Tamar-Agha and Salman, 2015; Koshnaw et al., 2020a).

2.4 The Red Beds Series strata

The Red Beds Series (RBS) is located along the suture zone in the hinterland of the NW Zagros fold–thrust belt on the Arabian Plate and overthrust by the Main Zagros fault (Fig. 1b). The RBS deposits consist of several NW–SE oriented discrete outcrops where alternating mudstones, sandstones, conglomerate beds, and carbonate beds are encountered. They have been grouped into three major units: the Suwais Group, the Govanda Formation, and the Merga Group. The Suwais and Merga groups have been further subdivided into four and two subgroups, respectively, based on their lithofacies. Generally, the thickness of the clastic Suwais Group varies geographically between 1000 and 1500 m (Jassim et al., 2006). In some locations, this unit is even thicker. Overall, the Suwais Group represents a succession of dark-red mudstones, siltstones, and sandstones with frequently occurring conglomeratic beds. Generally, the suite coarsens and thickens upward so that the conglomerate beds with cobbles and boulders become the dominant lithofacies toward the upper part. The Suwais Group is regarded as marine to mostly nonmarine fluvial deposits (Jassim et al., 2006; Alsultan and Gayara, 2016). A Paleocene–Eocene depositional age has been inferred for this unit based on fossil contents; however, no fossil content list has been reported (Jassim et al., 2006). Recently, using detrital zircon ages as criteria, Koshnaw et al. (2019) inferred a maximum depositional age as old as the late Oligocene ( 26 Ma).

The Suwais Group conformably overlies the Maastrichtian Tanjero Formation, but it was considered to follow the Lower Cretaceous Qulqula Formation with an unconformity (Karim et al., 2011; Hassan et al., 2014). However, using the late-Oligocene maximum depositional age of the Suwais Group as criteria, Koshnaw et al. (2019) interpreted the contact with the Tanjero Formation as unconformable. The upper contact of the Suwais Group with the Govanda Formation is reported as unconformable (Jassim et al., 2006). The Govanda Formation comprises  100–120 m thick reefal limestones with a few sandstone and conglomerate interbeds. Based on fossiliferous records, a Burdigalian to Langhian depositional age was proposed (van Bellen et al., 1959; Abdula et al., 2018). The  200–500 m thick Merga unit, which overlies the Govanda Formation (Fig. 2), consists of red mudstones and sandstones, some conglomerate interbeds with cobbles and boulders particularly in the upper  200 m of the Merga Group. This group is considered to have been deposited in a fluvial to alluvial environment (Jassim et al., 2006; Alsultan and Gayara, 2016). The contact between the Govanda Formation and Merga Group is conformable (van Bellen et al., 1959; Jassim et al., 2006).

Along-strike of the Zagros suture zone toward Iran, the potential equivalent of the RBS, the Razak Formation, appears not to be preserved in the Lurestan segment, but in the Dezful and Fars segments it has been documented (James and Wynd, 1965; Alavi, 2004; Khadivi et al., 2010, 2012; Vergés et al., 2019; Etemad-Saeed et al., 2020). The Razak Formation shows a lithostratigraphy and a timing of deposition that are similar to the RBS, particularly in the hinterland. The stratigraphy includes carbonates and calcareous argillites toward the foreland, and shales, siltstones, sandstones and conglomerates toward the suture zone (Alavi, 2004). Available magnetostratigraphic chronologies, between the High Zagros fault and the Mountain Front Flexure but not adjacent to the Main Zagros fault, suggest 19.7–16.6 Ma as the time of deposition of the Razak Formation (Khadivi et al., 2010). The older deposits of the Razak Formation closer to the Main Zagros fault might have been eroded due to deformation and exhumation (Alavi, 2004; Khadivi et al., 2012). These characteristics regarding the lithology and the timing of deposition suggest that the Razak Formation could correspond to the Govanda Formation and the Merga Group of the RBS. However, unlike the RBS in the Kurdistan region of Iraq, the Razak Formation seems to be more widespread geographically from the hinterland toward the foreland, and it appears in direct contact with the Asamri (Jeribe), Gachsaran (Fatha), and Aghajari (Injana) formations (Alavi, 1994; Khadivi, 2010; Etemad-Saeed et al., 2020).

3 Structural and stratigraphic relationships with the Red Beds Series

In the NW Zagros suture zone, the Red Beds Series (RBS) deposits are overthrusted by the Qandeel Series, the Walash Group, and the Naopurdan Group (Fig. 3). Along-strike, these units are interchangeably thrust onto the RBS, and they mark the Main Zagros fault (Fig. 1b). In the vicinity of the Main Zagros fault along the Tanun anticline (Figs. 4, 5), angular unconformities are preserved between the older Upper Cretaceous formations and (i) the middle-Miocene Govanda Formation on the NW plunge (Fig. 5b), and (ii) the upper-Oligocene–lower-Miocene (?) Suwais Group on the SE plunge (Fig. 5c). Along-strike, farther to the SE, an angular unconformity between the Red Beds Series (Suwais Group) and the Lower Cretaceous Qulqula radiolarites is also preserved (Karim et al., 2011). Structurally, the highest unit is the Qandeel Series, which is an ophiolite-bearing Upper Cretaceous metamorphosed klippe that likely disjointed by an out-of-sequence deformation from the adjacent terranes of the Sanandaj–Sirjan zone (Jassim et al., 2006; Ali et al., 2014, 2019). The Walash Group (Figs. 3e and 6a, b), situated below the Qandeel Series thrust sheet (Figs. 3, 4), is a volcaniclastic unit that is made up of dark-red mudstones, shales, slates, relatively limited beds of sandstones and conglomerates, radiolarites and limestone beds, and metasediments near the top. These beds are associated with pyroclasts, basic pillow lavas (diabase rocks, pyroxene-bearing spilitic basalts, and spilites) and basaltic andesite flows (pyroxene, andesite, and amphibole andesites) (Jassim et al., 2006; Ali et al., 2013). Pyroclastic feldspar grains of the Walash Group have been dated by the 40Ar/36Ar method and yielded 43.01 ± 0.15 Ma (Ali et al., 2013). Additionally, the 40Ar/36Ar age of the andesite lava flow is 43.1 ± 0.3 Ma, and the diabase rock is 40.1 ± 0.3–32.3 ± 0.4 Ma (Aswad et al., 2014). The isochronal age of 87Rb/86Sr yielded  44 Ma (Koyi, 2009). Depending on the planktonic and benthic foraminiferal assemblages, a marine environment with a depositional age of middle Eocene (Lutetian) has been proposed (Al-Banna and Al-Mutwali, 2008). The Walash Group itself was thrusted onto the Naopurdan Group, which is mostly made up of metasedimentary rock but includes basalt flows and pillow lavas as well. The Naopurdan Group (Figs. 3e and 6c–e) is composed of sandstones, shales, and conglomerates with basic volcanic, nummulitic limestone, and slate constituents. The upper part consists of an upward coarsening succession from sandstones to conglomerates with basic volcanic clasts (Jassim et al., 2006; Ali et al., 2013). Based on the planktonic foraminiferal assemblage, a marine environment of Eocene age has been suggested (Jassim et al., 2006), but 40Ar/36Ar ages measured on basaltic feldspars yielded younger ages between 24.31 ± 0.60 Ma and 33.42 ± 0.44 Ma (Ali et al., 2013). The Walash and Naopurdan groups were correlated with similar rock units in the adjacent part of the Zagros belt in Iran and were named Walash–Naopurdan–Kamyaran (WNK) (Ali et al., 2014; Moghadam et al., 2020). In the study area, the Naopurdan Group was thrust onto the autochthonous Red Beds Series deposits (Figs. 3, 4, 6) on the Arabian Plate. Along the thrust, the Naopurdan Group shows clear evidence of tectonic deformation, such as slate pencil cleavage (Figs. 3e, 6e). At a different locality below the Naopurdan thrust sheet, the Merga Group shows a well-developed shear zone with C and S planes (Figs. 3e, 6f).

Figure 3Landscape and outcrop photos exhibiting (a–d) the main transported thrust sheets and their interpreted direction of motion (c, top SW sense of shear) in the NW Zagros suture zone. (e) Simplified geologic map showing field photograph locations (yellow star).

Figure 4Geologic map illustrating key structural features, the growth strata location of the Tanjero Formation (Le Garzic et al., 2019), and the angular unconformity locations in the vicinity of the NW and SE plunges of the Tanun anticline. The base maps were sourced from Sissakian et al. (1997), MNR (2016), and Le Garzic et al. (2019).

Figure 5Landscape photographs showing field relationships. (a) The transported Paleogene Walash Group thrusted onto the Miocene units of the Red Beds Series that deposited on the Maastrichtian Tanjero Formation. (b) Angular unconformity between the middle-Miocene Govanda Formation of the Red Beds Series and the folded Cretaceous units. (c) Angular unconformity between the upper-Oligocene–lower-Miocene (?) Suwais Group of the Red Beds Series and the folded Cretaceous units. (d) Thrust fault between the Lower Jurassic strata (Sehkaniyan Formation?) and the Upper Cretaceous strata (Tanjero Formation).


Figure 6Landscape and outcrop photos displaying field examples of (a, b) the Walash Group, (c, d, e) the Naopurdan Group, and the Red Beds Series units of (f, h) the Merga Group and (g) the Suwais Group. The pencil cleavage tectonic fabric (e) developed in the lower part of the Naopurdan Group that thrusted onto the Merga Group. The C–S fabric developed in the upper part of the Merga Group that was overlain by the thrusted Walash Group.


4 Detrital zircon geochronology and thermochronology

4.1 Sampling and methods

In this paper, 1097 new detrital zircon U–Pb ages are presented from 11 samples (Tables S1 and S2 in the Supplement), which include 8 samples from the Red Beds Series and 3 samples from the proto-Zagros formations. These new data are integrated with previously published U–Pb ages in the study area (Koshnaw et al., 2019, 2020a). Additionally, for the purpose of the zircon (U–Th) / He double dating, 74 detrital zircons were selected for conducting new (U–Th) / He analyses from these geochronologically dated grains. These minerals were extracted from five Red Beds Series samples (Tables S1 and S3 in the Supplement).

All sample preparation and analyses on the new samples were conducted at The University of Texas at Austin Geo- and Thermochronometry and Mineral Separation laboratories. Mineral separation was conducted following mineral separation procedure (e.g., Gehrels, 2000), and  120 zircon grains of different size (30–100 µm) and shape were selected randomly from each sample for further analysis. The zircon U–Pb ages were produced using laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) following the procedures outlined in Marsh and Stockli (2015) and Hart et al. (2016). In order to preserve the zircon grains for future (U–Th) / He analysis, the detrital zircon geochronologic ages were obtained by depth-profile U–Pb analysis on unpolished grains (Stockli, 2017). The reported zircon U–Pb ages in this study represent discordance ages < 20 % with 2σ analytical uncertainty. The age concordance was determined based on the age of 206Pb/238U for zircon grains < 950 Ma, and the 207Pb/206Pb age for zircon grains > 950 Ma. After completion of the U–Pb analytical process, 13–21 representative zircon grains (except sample 12KRD-132, n=7) from major age peaks were manually picked from the tape mount and packed in a platinum (Pt) pocket for zircon (U–Th) / He thermochronometric analysis. In this study, the presented zircon helium ages are α-ejection corrected (Farley et al., 1996) with a standard procedure uncertainty of  8 % and 2σ (Reiners et al., 2002). The zircon (U–Th) / He ages were acquired following the methodology described in Wolfe and Stockli (2010). For detailed laboratory and data reduction processes, please refer to Thomson et al. (2017), Xu et al. (2017), and Pujols et al. (2020).

4.2 Detrital zircon U–Pb geochronology results

The detrital zircon U–Pb age distributions from three new proto-Zagros foreland basin samples manifest a unimodal U–Pb age spectra of Early Cretaceous to Late Cretaceous time ( 120–85 Ma) throughout the Tanjero (12KRD-142:  120–90), Kolosh (12KRD-143:  105–85), and Gercus (12KRD-146:  110–85) formations, with an age peak of  100 Ma (Fig. 7b). Among these, the Kolosh sample 12KRD-143 contains two relatively young grains of middle-Eocene ( 40 Ma) age. Interestingly, the detrital zircon U–Pb age populations that are determined for the Upper Cretaceous–Eocene proto-Zagros foreland basin fill are very different compared with the older Paleozoic and Lower Cretaceous pre-Zagros orogeny deposits and the younger Neogene Zagros foreland fill (Fig. 7a, c) (Koshnaw et al., 2020a). The detrital zircon U–Pb age results from the eight new Red Beds Series samples display a bimodal signature with late-Eocene ( 35–45 Ma) and Paleocene ( 55–65 Ma) ages with peaks of  40 Ma and  60 Ma, respectively (Fig. 8). Less commonly, an older age peak of  100 Ma is also found, except for a sample from the Suwais Group (SH17S2) that shows a dominant  100 Ma peak.

Figure 7Detrital zircon U–Pb age distribution plots for the (a) Neogene Zagros foreland basin, (b) proto-Zagros foreland basin, and (c) pre-Zagros orogeny deposits depicted as probability density plot (PDP; bandwidth =20, bin width =20) and age histograms (Vermeesch, 2012). Colored shading represents the key source terrane age components. Samples with an asterisk are new samples. Published samples are from Koshnaw et al. (2017, 2020a). Sample locations are shown in Fig. 1b.


Figure 8Detrital zircon U–Pb age distribution plots for the Red Beds Series units of the Suwais and Merga groups depicted as probability density plot (PDP; bandwidth =20, bin width =20) and age histograms (Vermeesch, 2012). Dark-red vertical bars delineate the youngest age peak ( 3 grains). Colored shading represents key source terrane age components (see Fig. 7 for legend). Samples with an asterisk are new samples. Published samples are from Koshnaw et al. (2019). Sample sets are separated by a white space according to their localities. Sample locations are shown in Fig. 1b.


4.3 Detrital zircon (U–Th) / He thermochronology and (U–Th) / He double-dating results

To further constrain the potential source regions for the Red Beds Series deposits based on the zircon (U–Th) / He double dating, five samples were chosen for the zircon (U–Th) / He thermochronometric analysis (ZHe). Two samples were collected from the NW of the study area and three samples from the SE of the study area, including the Suwais and Merga groups in each area (Fig. 9). Detrital zircon grains were selected from each main U–Pb age peak including  40 Ma (late Eocene),  60 Ma (middle Paleocene),  75 Ma (Late Cretaceous),  100 Ma (late Early Cretaceous),  170 Ma (Jurassic),  300 Ma (Carboniferous–Permian), and  600 Ma (Neoproterozoic). Generally, the bulk ZHe age results display notable ranges: in the  66–55 Ma range, the most dominant age is  60 Ma; in the  50–35 Ma range, the most dominant ages are  40 Ma (for the Suwais Group) and  45 and  37 Ma (for the Merga Group) (Fig. 9a). Two grains from the Merga Group that were extracted from samples MT17M3 and SH17M5 chronicle the youngest exhumation at  21 Ma. From the Suwais Group, three ZHe ages that overlap within error mark the youngest peak of  29 Ma. The detrital zircon (U–Th) / He double-dating results reveal cooling ages of  40–30 Ma for the late-Eocene grains and  60–40 Ma for the middle-Paleocene grains (Fig. 9b, c). The Late and late Early Cretaceous zircon grains have cooling ages of  70–35 Ma. The Jurassic zircon grains show cooling ages of  60–40 Ma, which are comparable to the middle-Paleocene zircon grains. The Neoproterozoic and Carboniferous–Permian zircon grains manifest a wide range of cooling ages of  150–30 and  120–50 Ma, respectively.

Figure 9(a) Detrital zircon (U–Th) / He age distribution plots for the Red Beds Series units of the Suwais and Merga groups, depicted as kernel density estimation (KDE; bandwidth =2, bin width =5) and age histograms (Vermeesch, 2012). (b) Detrital zircon (U–Th) / He double-dating plots of individual zircon analysis. (c) Expanded view of the 0–120 Ma period to show detrital zircon (U–Th) / He ages for the late-Mesozoic and Cenozoic times.


5 Discussion

5.1 Sedimentary provenance evolution

5.1.1 Provenance of the NW Zagros wedge-top deposits

The new detrital zircon U–Pb age results form this study highlight distinctive signatures for each of the various tectonic phases including (i) the passive margin of the pre-Zagros orogeny, (ii) the proto-Zagros foreland basin stage related to the obduction of the Neotethys oceanic crust, and (iii) the Neogene Zagros basins related to the Arabia–Eurasia continental collision (Fig. 10a, b). In particular, the provenance data from the Maastrichtian Tanjero (Amiran), the Paleocene–Eocene Kolosh (Amiran), and the Eocene Gercus (Kashghan) formations, all of which were deposited in the proto-Zagros foreland basin, show a strong influence of Neotethys-related terranes and a striking difference from the older and younger formations, implying a switch of the main sediment sources through time. The key source terranes in the NW Zagros region and the surrounding areas are (i) the Paleocene–Eocene Walash–Naopurdan–Kamyaran complex (WNK), (ii) the Cretaceous Neotethys oceanic crust and the island arc ophiolitic terranes, (iii) the Jurassic and Triassic igneous rocks within the Sanandaj–Sirjan zone (SSZ), (iv) the Carboniferous–Permian northern Gondwana and Variscan rocks, and (v) the Neoproterozoic Pan-African/Arabian–Nubian shield and older rocks (Fig. 7).

Figure 10Comparison plots of the pre-Zagros orogeny, proto-Zagros foreland basin, Neogene Zagros foreland basin, and Red Beds Series deposits in the present-day Zagros hinterland: (a) detrital zircon U–Pb cumulative age distribution (Saylor and Sundell, 2016), and (b) two-dimensional multidimensional scaling plot depicting relative statistical similarity/dissimilarity based on the Kolmogorov–Smirnov test D value (Vermeesch, 2013) of the detrital zircon U–Pb ages for the same sample sets as in the cumulative age distribution plot. Dimension 1 and 2 are unitless.


Samples from the proto-Zagros foreland basin formations that were obtained from the NW of the study area (Fig. 1b) have a sole age peak of  100 Ma. This age peak points to the Neotethys oceanic crust and island arc ophiolitic terranes as principal sources during the Maastrichtian and the Paleogene times. The pre-Zagros passive margin strata (Fig. 7) from the northern part of the study area have detrital zircon U–Pb age distributions dominated by the early Neoproterozoic, Pan-African, Cadomian (peri-Gondwana), and other Gondwana-related sources (Stern and Johnson, 2010; Alirezaei and Hassanzadeh, 2012; Stern et al., 2014; Avigad et al., 2016, 2017; Golan et al., 2018; Koshnaw et al., 2017; Al-Juboury et al., 2020; Omer et al., 2021). These age components are also found in the detrital zircon populations of the Lower Cretaceous formations in southern Iraq (Wells et al., 2017). The late-Paleozoic age components from the pre-Zagros formations are likely Gondwana-related, unlike the comparable age components from the younger formations that likely involve Variscan-derived detritus (e.g., Barber et al., 2019). In general, the sediment of the Paleozoic and Lower Cretaceous rocks largely originated from (i) the Arabian–Nubian shield; (ii) the Gondwana super fan system, which itself was sourced from the East African Orogen; or from (iii) the SE Saharan Metacraton (Koshnaw et al., 2017; Meinhold et al., 2021).

Overall, the contrast in the detrital zircon provenance signature between the pre-Zagros and proto-Zagros deposits reflects a fundamental reversal in the sediment source from the west to the east as the basin started to flex toward the east. This flexural downbending most likely reflects the arrival of the Neotethys oceanic crust and ophiolitic terranes, which were associated with numerous Late Cretaceous magmatism events, including granitoid intrusions with zircons of a middle-Cretaceous age (Delaloye and Desmons, 1980; Ali et al., 2012; Nouri et al., 2016; Al Humadi et al., 2019; Ismail et al., 2020). In the SE of the study area, detrital zircons from the Tanjero Formation with Paleozoic and Proterozoic ages as well as a few grains with Triassic and latest-Paleocene ages have been reported (Jones et al., 2020). Farther SE, in Iran, the proto-Zagros strata of the Amiran and Kashghan formations have  240 Ma Triassic and  100 Ma middle-Cretaceous age signatures in their detrital records (Zhang et al., 2017; Barber et al., 2019). These additional old detrital zircon age components potentially imply a local variation in the paleogeography along-strike. Even though zircon is not a common mineral in mafic rocks, it has been recognized (Grimes et al., 2007, and references therein). The  240 Ma Triassic age signature in the detrital zircon records of the Amiran and Kashghan formations has been attributed to the mid-oceanic ridge based on the trace element data measured in zircon grains (Barber et al., 2019). Additionally, sandstone petrographic analysis on the Tanjero Formation revealed the occurrence of limestone and chert grains that were particularly sourced from the older Qulqula radiolarite and Balambo formations in the NW (Aziz and Sadiq, 2020; Jones et al., 2020; Çelik and Salih, 2021). Such variation in the sediment source for the proto-Zagros foreland basin is in line with (i) the destruction of the Gotnia Basin architecture, (ii) an input of carbonate material into the new flexural basin, and (iii) the occurrence of some recycled Paleozoic and older zircon grains in the Tanjero Formation (Aqrawi et al., 2010; English et al., 2015; Aziz and Sadiq, 2020; Jones et al., 2020). The younger Neogene Zagros foreland basin strata also have distinct U–Pb age spectra that are dissimilar to the proto-Zagros and the pre-Zagros U–Pb age distributions. In the same locality where the new proto-Zagros samples were collected (Fig. 1b), published data show that most of the Neogene formations, except the Injana Formation, are characterized by an Eocene ( 40 Ma) age peak with subordinate contributions of older ages to the age spectra (Fig. 7). However, in different localities, the Neogene deposits show a more mixed signature (Koshnaw et al., 2020a). These sediments with zircon grains recording a dominant Eocene ( 40 Ma) age peak were interpreted to have been supplied by transverse rivers from the NE during a late stage of the Zagros orogeny (Koshnaw et al., 2020a). In addition to the Eocene age peak, the Injana Formation records a wider range of source terranes where the detrital zircons are characterized by younger early-Oligocene ages and older age components. This suggests that the sedimentary material was likely recycled from older uplifted strata and delivered by axial rivers from the NW (Koshnaw et al., 2020a).

Ultimately, such variation in the detrital zircon U–Pb provenance between the different hinterland wedge-top units of the NW Zagros Basin denote that (i) the pre-Zagros passive margin deposits are diagnosed by a wide age range of old zircon grains, (ii) the proto-Zagros sediments that are related to the Neotethyan obduction record a unique zircon age component of Albian–Cenomanian, and (iii) the Neogene Zagros foreland basin deposits are dominated by zircon grains with Eocene and early-Oligocene ages.

5.1.2 Provenance of the NW Zagros Red Beds Series deposits

The Red Beds Series (RBS) provenance data put an age limit on the diminution of the magmatic activities and on the exhumation of the source terranes. The detrital zircon U–Pb age populations from the RBS depict two characteristic age components, which are Paleocene ages with a peak of  60 Ma and late-Eocene ages with a peak of  40 Ma (Fig. 8). Additionally, older subordinate age peaks similar to the pre- and proto-Zagros age components are also present. Zircon grains from the Paleocene and Eocene age components record identical (U–Th) / He and U–Pb ages, pointing to a volcanic origin (Fig. 9b, c). Furthermore, their cooling over  60–40 Ma for the middle-Paleocene grains and  40–30 Ma for the late-Eocene grains suggest the occurrence of a systematic exhumation and unroofing of magmatic plutons through time.

The older subordinate age components possibly reflect the combined effect of recycling of the previously deposited strata as shortening continued and local variations in rock outcrops along-strike, including potential exposure of the ophiolitic rocks. The potential source terrane candidates that have ages comparable to the Paleocene and Eocene age components are the Walash–Naopurdan volcaniclastic complex and Urumieh-Dokhtar magmatic zone (UDMZ), which both are located to the NE of the study area. The Paleogene Walash–Naopurdan volcaniclastic complex extends along-strike into Iran and was previously termed the Walash–Naopurdan–Kamyaran (WNK) (Ali et al., 2014; Moghadam et al., 2020). The WNK complex potentially extends into Turkey as well and is known as the Maden–Hakkari complex (Braud and Ricou, 1975; Yılmaz, 1993; Robertson et al., 2007; Oberhänsli et al., 2010). In Iran, the Walash–Naopurdan complex is equivalent to the volcaniclastic deposits near Kamyaran and Kermanshah, which have been named differently as Gaveh-Rud domain (Braud and Ricou, 1975; Homke et al., 2010; Saura et al., 2015), early-Tertiary magmatic domain (Agard et al., 2005, 2011), and Kamyaran Paleocene–Eocene arc (Whitechurch et al., 2013). These Paleogene rocks represent arc-related (defined differently as arc, back-arc, and forearc) volcanic and volcaniclastic units (Homke et al., 2010; Ali et al., 2013; Whitechurch et al., 2013; Saura et al., 2015), which likely developed in the proximity of the SW margin of the Eurasian Plate in association with the subduction of the Neotethys Plate beneath Eurasia (Robertson et al., 2007; Barrier and Vrielynck, 2008; Homke et al., 2010; Agard et al., 2011; Whitechurch et al., 2013; Saura et al., 2015). The dated magmatic bedrock within the WNK complex and part of the adjacent Sanandaj–Sirjan zone (SSZ) show both age components that are found in the RBS, i.e., Paleocene (Braud, 1987 in Homke et al., 2010; Azizi et al., 2011; Whitechurch et al., 2013) and Eocene (Mazhari et al., 2009, 2020; Azizi et al., 2011, 2019; Ali et al., 2013; Aswad et al., 2014; Moghadam et al., 2020). The Eocene intrusions are also found within the Upper Cretaceous ophiolitic rocks (Leterrier, 1985; Aswad et al., 2016; Ali et al., 2017; Ismail et al., 2020). Younger Oligocene magmatic rocks within the WNK complex have also been reported (Whitechurch et al., 2013; Ali et al., 2013). The UDMZ contains large Eocene and Oligocene igneous bodies marking a prolonged magmatic flare-up (Verdel et al., 2011; Chiu et al., 2013) with a peak magmatism age of  40 Ma (van der Boon et al., 2021). Accordingly, both the WNK complex and the UDMZ could have served as source areas for the Red Beds Series. However, the WNK complex was thrusted onto the RBS and is adjacent to the SSZ, whereas the SSZ is situated between the RBS and the WNK complex to the west and the UDMZ to the east, which makes any sedimentary pathways from the UDMZ more complicated. Additionally, the regional uplift of the SSZ occurred earlier than that of the UDMZ; this makes it very unlikely that large volumes of UDMZ material were supplied to the depocenter of the Red Beds Series since uplift of the SSZ disconnected the sedimentary pathways (François et al., 2014; Barber et al., 2018). Moreover, the majority of the double-dated detrital zircons from the RBS that have Ediacaran, Carboniferous–Permian, and Jurassic ages show an exhumation age that corresponds to the time span between the Late Cretaceous and the Eocene (Fig. 9b), which is comparable to the exhumation age of the Sanandaj–Sirjan zone (e.g., Homke et al., 2010; Khadivi et al., 2012; Mouthereau et al., 2012; Barber et al., 2018). Such an exhumation pattern of the old zircon grains that are found in the WNK-derived RBS supports the interpretation that the development of the WNK complex occurred in the vicinity of the Sanandaj–Sirjan zone.

The new data presented in this paper in combination with the previous data by Koshnaw et al. (2019) of the RBS illustrate (Fig. 8) that the youngest zircon age peak component of the Suwais Group samples is between  37 and  41 Ma. In the same sense, the youngest zircon age peak component of the Merga Group samples is between  36 and  45 Ma; however, a younger  33 Ma peak (three grains) also exists. Along-strike, both groups show older ages toward the SE, similar to the age trend of the Walash and Naopurdan rocks,  24–33 Ma in the NW and  32–43 Ma in the SE (Ali et al., 2019). Published detrital zircon U–Pb data (Jones et al., 2020) of three samples collected from the Tanjero Formation (T, n=28) and the Red Beds Series (LR, n=5; UR, n=39) but farther to the SE of our study area show disparity with the provenance results from this study. This mainly concerns the samples from the Tanjero Formation and lower Red Beds (LR), as they lack the middle-Cretaceous and the Paleogene age components, respectively. This variation is likely due to the limited number of the analyzed grains, as dating less than 117 detrital zircon grains for provenance studies increases the chance of not detecting age components that consist of more than 5 % of the population at a 95 % confidence level (Vermeesch, 2004).

So far, in the vicinity of the study area, late-Oligocene magmatism has only been documented in two localities in the NW Zagros area, including the Kurdistan region of Iraq and Iran (Ali et al., 2013; Whitechurch et al., 2013); however, in Iran (and thus farther NE of the study area), magmatism is more abundant (Agard et al., 2011; Chiu et al., 2013). In general, magmatic activity appears to be considerably lower during the late Oligocene than during the Eocene (van der Boon et al., 2021). In the Kurdistan region of Iraq, earlier partial cessation of magmatism in the SE part is potentially due to the presence of the Bisotoun block (Agard et al, 2005; Wrobel-Daveau et al., 2010; Sasvari et al., 2015), which implies that the thrusting of the WNK terranes occurred earlier farther SE of our study area. As the geographic extension of the Bisotoun block is limited, this local difference in the timing of the WNK terrane thrusting cannot be in contradiction to the regional oblique Arabia–Eurasia collision that potentially commenced first farther to the NW, away from the study area (e.g., Rolland et al., 2012; Darin et al., 2018).

Overall, the U–Pb age populations of the Red Beds Series deposits show consistent Paleocene and Eocene age signatures from the bottom to the top, which imply that the WNK complex was a major source area and that the sediment sources have not changed significantly during deposition of the RBS.

5.2 Basin dynamic of the NW Zagros Red Beds Series deposits

The angular unconformity between the Suwais Group and the Upper Cretaceous strata, and the detrital zircon U–Pb and (U–Th) / He double-dating results advocate for deposition of the RBS in a distinctive hinterland basin.

On the basis of the documented growth strata in the Maastrichtian Tanjero Formation of the Tanun anticline (Le Garzic et al., 2019) and the lack of growth strata in the Suwais Group (Figs. 4, 5c), the Suwais Group deposition must have taken place after the deposition of the Tanjero Formation and the initial growth of the Tanun anticline. This period of no stratigraphic record is also evident from the  26 Ma detrital zircon maximum depositional age of the Suwais Group (Koshnaw et al., 2019) and the Maastrichtian age of the Tanjero Formation (Abdel-Kireem, 1986). However, no direct angularity is manifested between the Tanjero Formation and the Suwais Group (Le Garzic et al., 2019), possibly because of erosion of the tilted beds of the Tanjero Formation where the Suwais Group is deposited and the preservation of the semi-horizontal beds. The angular unconformity between the Suwais Group and the Upper Cretaceous beds (Fig. 5c) along with other angular unconformities between the Suwais Group and the Qulqula Formation (Karim et al., 2011), and the Govanda Formation and the Upper Cretaceous beds (Fig. 5b) suggest that the RBS unit was deposited in an intermontane setting after the development of the proto-Zagros fold–thrust belt (Lawa et al., 2013; Koshnaw et al., 2017, 2020b; Le Garzic et al., 2019).

In this study, the comparison of the detrital zircon U–Pb provenance data established for the RBS, pre-Zagros, proto-Zagros, and Neogene Zagros sediments reveals a robust dissimilarity of the RBS deposits with the pre-Zagros and proto-Zagros units but a possible linkage with the Neogene Zagros units (Fig. 10). On the basis of the sandstone petrography, the Suwais Group and the Tanjero Formation are considered as litharenites, although a key difference is the presence of the magmatic and metamorphic grains in the Suwais Group and their scarcity in the Tanjero Formation (Hassan et al., 2014; Çelik and Salih, 2018, 2021).

To further constrain the tectonic setting of the RBS, three principal hypotheses can be tested: (a) the RBS was deposited in a foreland basin, (b) the RBS was deposited in a foreland basin but later developed into an intermontane basin, (c) and the RBS was deposited in an intermontane basin (Fig. 11). To render a unique solution, the detrital zircon (U–Th) / He double-dating method has been utilized. In scenario (a) (Fig. 11a), the RBS deposits are expected to have a similar signature to the proto-Zagros foreland basin where the detrital zircon U–Pb age spectra is dominated by a  100 age peak (this study) and characterized by a  70 Ma zircon (U–Th) / He age peak (Barber et al., 2019). For scenario (b) (Fig. 11b), where the RBS deposition would take an early-stage foreland and a late-stage intermontane basin position, we would expect zircon grains with a mixed U–Pb age distribution, which should reflect a mixture of (i) recycled grains with  100 and  70 Ma zircon (U–Th) / He age peaks, and (ii) new  60 and  40 Ma grains with their corresponding zircon (U–Th )/ He ages. In scenario (c) (Fig. 11c), where the tectonic setting of the RBS is unrelated to the foreland basin, a mixed U–Pb age spectra similar to scenario (b) is predicted, although with no inherited  70 Ma zircon (U–Th) / He age because no recycling of foreland basin strata is expected. Scenario (a) can be ruled out, as the RBS detrital zircon U–Pb age spectra is substantially different (Figs. 7, 8, 10). Scenarios (b) and (c) offer valid alternatives, but there are two key differences in the U–Pb and (U–Th) / He ages. In scenario (b), because the RBS basin involves an early-stage foreland basin, no detrital zircons with Paleocene and Eocene ages should be expected in its lower strata. However, samples from the Suwais Group are dominated by the Paleocene and Eocene age components (Fig. 8). Additionally, the detrital zircon (U–Th) / He double-dating results (Fig. 9) show that most of the  100 Ma U–Pb age grains have a  65 Ma and  40 Ma exhumation age that is comparable to the Paleocene and Eocene age grains. These ages are still younger than the  70 Ma zircon (U–Th) / He age of the proto-Zagros strata. Moreover, due to the presence of the  40 Ma zircon (U–Th) / He age and insufficient burial in the proto-Zagros basin to reset the zircon (U–Th) / He ages, they are unlikely to record an origin from the deformed proto-Zagros strata. Indeed, the pre-deformational thickness of the collective proto-Zagros formations (Tanjero, Kolosh, and Gercus formations) is about 3–3.6 km (van Bellen et al., 1959). Assuming 20  for the surface temperature and 25 C km−1 for the geothermal gradient, about 6–7 km thickness would be required to reset the zircon (U–Th) / He ages. Accordingly, the  100 Ma U–Pb ages encountered in detrital zircon grains rather reflect the supply of material from the newly eroded ophiolitic terranes where igneous intrusions have a similar age. Therefore, on the basis of the detrital zircon U–Pb signature and the (U–Th) / He double-dating results, we argue that scenario (c), the intermontane basin hypothesis, is more plausible and compatible with the provenance characteristics of the RBS deposits.

Figure 11Predicted detrital zircon (U–Th) / He double-dating ages for the competitive hypotheses of the Red Beds Series depositional setting: (a) a foreland basin; (b) a foreland basin, but later partitioned and evolved into an intermontane basin; (c) an intermontane basin.


All things considered, the angular unconformities and the provenance information suggest neither a significant connection of the RBS intermontane basin with the proto-Zagros foreland basin nor a potential connection during the Miocene with the Neogene Zagros foreland basin, possibly as a result of a broad Zagros collisional zone where uplift and exhumation occurred. This interpretation about the position of the RBS in the NW Zagros area in the Kurdistan region of Iraq is most likely also valid for the Razak Formation in the Dezful and Fars areas of Iran, although considering a possible variation in the structural architecture of the basin. In the study area, the RBS is structurally bound by the allochthonous WNK complex toward the NE and by the anticlines of the proto-Zagros fold–thrust belt toward the SW (e.g., Figs. 4 and 5), and the RBS deposits overlay the Paleogene and older rocks in direct contact. However in the Dezful and Fars segments of the Zagros in Iran, where the development of the Paleogene proto-Zagros fold–thrust belt is limited (Hessami et al., 2001), the Razak Formation deposits occur in direct contact with the Neogene Zagros foreland basin deposits with no documented unconformities (James and Wynd, 1965; Khadivi et al., 2010; Vergés et al., 2019; Etemad-Saeed et al., 2020). Such a basin setting suggests that the Razak Formation could be interpreted as an early foreland basin unit, at least during the Neogene. In either case, we emphasize that both the RBS and the Razak Formation are represent the deposits of the Arabia–Eurasia early-stage collision (Khadivi et al., 2012; Koshnaw et al., 2019; Etemad-Saeed et al., 2020).

5.3 Geodynamic evolution

Combination of new field evidence with the provenance data and with the existing knowledge about the magmatic and tectonic history of the Arabian and Eurasian margins allows us to unravel the Cenozoic convergence history between the two plates in the present-day NW Zagros region (Figs. 12, 13).

After the initiation of the Neotethys intraoceanic subduction zone during the Aptian and its development until the middle Late Cretaceous, including the growth of ophiolitic island arcs (Aswad and Elias, 1988; Aziz et al., 2011; Ali et al., 2012; Barber et al., 2019), the ophiolitic terranes obducted on the NE margin of the Arabian Plate during the Campanian (Figs. 12, 13a). Sediments were shed from the uplifted terranes on the newly formed proto-Zagros foreland basin, as evidenced by the detrital zircon provenance data from this study and the detrital zircon (U–Th) / He ages from the Amiran and Kashghan formations from Lurestan Province in Iran (Barber et al., 2019). In the Upper Cretaceous Bekhma Formation (late-Campanian–early-Maastrichtian age; van Bellen et al., 1959), the U–Pb dating of hydrothermal cementation yielded  74 Ma (Salih et al., 2020), consistent with the detrital geochronological and thermochronological constraints (Fig. 12). These lines of evidence suggest that the proto-Zagros fold–thrust belt growth initiated as a result of the ophiolite obduction. During the Paleocene (Figs. 12, 13a), the WNK complex possibly originated as a back-arc to arc system, associated with magmatic activity (Braud, 1987 in Homke et al., 2010; Azizi et al., 2011, 2019; Ali et al., 2013; Whitechurch et al., 2013). The signature of the Paleocene magmatism is also reflected in the detrital zircon U–Pb and (U–Th) / He records of the Red Beds Series that originated from the WNK complex. The locality of this arc system is contested. In particular, it is not clear whether it developed in the vicinity of the Eurasian Plate (Braud and Ricou, 1975; Agard et al., 2005, 2011; Robertson et al., 2007; Barrier and Vrielynck, 2008; Oberhänsli et al., 2010; Al-Qayim et al., 2012; Whitechurch et al., 2013) or close to the Arabian Plate (Aswad et al., 2014; Ali et al., 2017, 2019). However, due to presence of Eurasian derived detritus (Homke et al., 2010; Whitechurch et al., 2013) and the occurrence of low-grade metamorphic grains (Ali et al., 2017), which are abundant in the Sanandaj–Sirjan zone (SSZ) (Hassanzadeh and Wernicke, 2016), a Eurasian-related setting that resulted from the oceanic slab rollback is more likely (Robertson et al., 2007; Agard et al., 2011; Whitechurch et al., 2013; Abdollahi et al., 2020).

Figure 12Time–space chart synthesizing hinterland basin stratigraphy and structural zone thermochronological and geochronological constraints on key tectonic and geodynamic, metamorphic, magmatic, and exhumational events in the broader NW Zagros sector. Numbered arrows represent the estimated relative timing of the thrusting of their respective terranes. Foreland basin stratigraphy is based largely on van Bellen et al. (1959), Aqrawi et al. (2010), and Al-Husseini (2000). The abbreviations used in the figure are as follows: MDA, maximum depositional age; AHe, apatite (U–Th) / He; ZHe, zircon (U–Th) / He; AFT, apatite fission track; ZFT, zircon fission track; MFF, mountain front flexure; HZF, High Zagros fault; MZF, Main Zagros fault; WNK, Walash–Naopurdan–Kamyaran; UDMZ, Urumieh-Dokhtar magmatic zone.


During the Eocene (Fig. 13b), the magmatism in the WNK complex and the SSZ edge continued; however, the SSZ magmatism was mostly associated with intrusions rather than volcanism (Moritz et al., 2006; Mazhari et al., 2009, 2020; Azizi et al., 2011, 2019; Ali et al., 2013; Whitechurch et al., 2013; Aswad et al., 2014; Abdollahi et al., 2020; Moghadam et al., 2020). This period of magmatism in the WNK complex was coeval with the nearby voluminous UDMZ magmatic flare-up (Verdel et al., 2011; Chiu et al., 2013; van der Boon et al., 2021), implying a genetic relationship. The downgoing oceanic slab was possibly subducted deeper and farther beneath the Eurasian Plate with a relatively shallower angle that led to the generation of magmas in the WNK and UDMZ at a comparable period (Whitechurch et al., 2013; Hassanzadeh and Wernicke, 2016). Additionally, in the proximity of the Arabian Plate, the Cretaceous ophiolitic terranes also contain Eocene intrusions (Fig. 13b; Leterrier, 1985; Aswad et al., 2016; Ali et al., 2017; Ismail et al., 2020), with an origin that is likely different from that of the WNK complex (Aswad et al., 2016). These Eocene intrusions in the Cretaceous ophiolitic terranes are likely the consequence of the break-off of the subducted oceanic slab. Tomographic data for the Middle East show the occurrence of cold material at deep levels that has been considered as a remnant from an older early-Paleogene slab break-off event (Agard et al., 2011; van der Meer et al., 2018). This inferred cold material is different from the remnant of a younger (Miocene) and shallower slab break-off (Mesopotamia and Zagros slabs; van der Meer et al., 2018).

Figure 13Schematic illustrations portraying the late-Mesozoic–Cenozoic history of the Arabia–Eurasia convergence and the NW Zagros suture zone evolution. The possible configuration of the subducted Neotethys oceanic slabs, and the Arabian and Eurasian Plate margins are constrained by the results of this study as well as previous data and interpretations (Agard et al., 2005; Barrier and Vrielynck, 2008; Omrani et al., 2008; Agard et al., 2011; Moghadam and Stern, 2011; Agard et al., 2011; Mouthereau et al., 2012; Whitechurch et al, 2013; Ali et al., 2014; Jagoutz et al., 2015; Hassanzadeh and Wernicke, 2016; and the publications cited in Fig. 12).


For this period, the Eocene, several models suggested multiple subduction zones with 200–300 km distance between the zones across the Neotethys to explain the documented magmatic activities (Aswad et al., 2016; Ali et al., 2017, 2019). However, due to a space problem and difficulties with subduction inception, a multiple subduction system has been considered unlikely (Whitechurch et al., 2013). Present-day natural examples of one single complete suite including a trench, a forearc, an arc, and a back-arc necessitate a horizontal width of 200–500 km, depending on the subduction angle (e.g., Cloos et al., 1993; Stern, 2002, 2010). In this study, the GPlates 2.2.0 plate kinematic reconstruction software (Müller et al., 2018) has been utilized, and the published plate circuit model of Arabia–Eurasia (McQuarrie and Hinsbergen, 2013) was used as a base for the reconstruction, with the assumption of  35 Ma as an onset of collision. Based on this setting, the calculated distance between the Arabian and Eurasian margins was estimated to be 283 ± 66 km at  40 Ma. For this calculation, the length of the Harsin (35 km) and Bisotoun (49 km) domains were excluded (Vergés et al., 2011), as the study area is farther NW of Lurestan Province and no evidence has been found regarding their occurrence. Therefore, we hypothesize that the 283 ± 66 km distance between the Arabian and Eurasian plates is more compatible with a convergence model that involves a subduction of the Neotethys slab below Eurasia as well as an older subducted slab below the obducted ophiolite on Arabia. This conclusion is in line with the plate kinematic reconstruction models in the broader Middle East (Dewey et al., 1973; Barrier and Vrielynck, 2008; Jagoutz et al., 2015; van Hinsbergen et al., 2020).

The latest Eocene to the late Oligocene ( 36–26 Ma) marks a period of significant reduction of magmatism in the UDMZ, SSZ, WNK, and the Late Cretaceous ophiolitic terranes (Figs. 12, 13c). During the same time, a widespread unconformity was recorded in the NW Zagros area in the Kurdistan region of Iraq (Dunnington, 1958; Ameen, 2009; Lawa et al., 2013). Furthermore, the development of the syndepositional compressional joints within the upper-Eocene Pila-Spi Formation has also been documented (Numan et al., 1998) in the Zagros fold–thrust belt. Lastly, fractures that were filled by hydrothermal cementation in the Zagros fold–thrust belt hint at a period of fluid flow during the early Oligocene ( 30 Ma), as constrained by the U–Pb method (Salih et al., 2020). Among the presented RBS samples, there are four young detrital zircon age peaks of  36–37 Ma (Fig. 8). Such a pattern in the detrital zircon ages signifies a remarkable attenuation of magmatism in the source terranes by  36 Ma. However, the occurrence of a single younger peak of  33 Ma, and limited younger magmatic activity of  26 Ma (Whitechurch et al., 2013) and  24 Ma (Ali et al., 2013) imply that magmatic activities continued locally. Samples from the RBS, particularly from the bottom section of the Suwais Group in the study area (e.g., Fig. 5c), include the Eurasian detrital zircon age signature and contain zircon grains as young as  26 Ma (Koshnaw et al., 2019). We argue that these data set the constraints on the timing of the Arabia–Eurasia continental collision between  36 and  26 Ma, which includes the entire geodynamic development from the commencement to the termination in the NW Zagros region, a conclusion which is in agreement with previous studies (Hempton, 1987; Agard et al., 2005, 2011; Horton et al., 2008; Allen and Armstrong, 2008; Ballato et al., 2011; Mouthereau et al., 2012; Mohammad et al., 2014; Barber et al., 2018).

After the continental collision, the Zagros suture as well as the SSZ and the UDMZ underwent notable uplift and exhumation starting in the earliest Miocene (Figs. 12, 13d). The youngest detrital zircon grains from the RBS deposits (Merga Group) record a  21 Ma zircon (U–Th) / He age, comparable to (i) the  23 Ma zircon (U–Th) / He age of an Eocene intrusion within the Cretaceous ophiolite (Ismail et al., 2020); (ii) the period of accelerated exhumation in the SSZ and UDMZ as documented by apatite (U–Th) / He (AHe), apatite fission track (AFT), and zircon fission track (ZFT) ages (François et al., 2014; Behyari et al., 2017; Barber et al., 2018); and (iii) a  22 Ma AFT detrital age from the Aghajari Formation of the Neogene Zagros foreland basin in Lurestan Province in Iran (Homke et al., 2010). Additionally, the earliest Miocene is also the time of alkaline magmatism initiation in the UDMZ (Homke et al., 2010, and references therein; Mouthereau et al., 2012, and references therein; van der Boon, 2021). Furthermore, the present-day field relationships in some localities reveal that the Cretaceous ophiolitic terranes, including the associated Eocene intrusions that have a different petrochemistry than the Eocene WNK, were thrusted onto the WNK (e.g., Figs. 3, 4) likely by an out-of-sequence mechanism (Ali et al., 2014; Aswad et al., 2016). Such an event likely occurred at a similar time to the documented  21–23 Ma exhumation age across the NW Zagros collisional zone as a consequence of the terminal collision. Moreover, in the SE Fars segment of the Zagros belt, the AFT detrital age from the potential stratigraphic equivalent of the RBS, the Razak Formation, yielded  25 Ma, and it has been interpreted as recording the exhumation of the Main Zagros hanging wall in the suture zone (Khadivi et al., 2012). Such a similarity in the suture zone exhumation during the latest Oligocene and the earliest Miocene across the Zagros  2000 km long orogen may suggest a uniform suturing between Arabia and Eurasia, and possibly reduction in the oblique collision effect and onset of the northward motion of Arabia (McQuarrie et al., 2003; Navabpour et al., 2014). Later, during the middle and late Miocene ( 14–4 Ma), further enhanced exhumation took place from the suture zone in the foreland direction, with sediment inputs from the NW and the NE of the study area that were possibly triggered by a new slab break-off event (Koshnaw et al., 2020a, b).

6 Conclusions

An angular unconformity between the Red Beds Series (RBS) units, the Suwais Group and the Govanda Formation, and the deformed Upper Cretaceous formations in the NW Zagros hinterland denotes deposition of the RBS in an intermontane basin after the development of the proto-Zagros fold–thrust belt. The detrital zircon U–Pb age data of the RBS units (Suwais and Merga groups) and the proto-Zagros foreland basin strata (Tanjero, Kolosh, and Gercus formations) are strikingly different, yet the RBS data are partially commensurate with the Neogene Zagros foreland basin fill, and all data sets are unrelated to the pre-Zagros Paleozoic and Lower Cretaceous record. The detrital zircon (U–Th) / He double-dating result pattern of the RBS further highlights that these sediments were deposited in an intermontane basin. The dating results additionally demonstrate a signal related to the supply of volcanic material and to the unroofing of potential plutons, dissimilar to the proto-Zagros detrital record. The evaluation of the provenance data shows a derivation of the RBS deposits from the Paleogene Walash–Naopurdan–Kamyaran complex that was developed along the Eurasian margin, whereas the proto-Zagros deposits were derived from the Upper Cretaceous ophiolitic terranes. Significant cutback of magmatism by  36 Ma, as reflected in the RBS youngest U–Pb age peaks, points to the latest Eocene (< 36 Ma) as the onset of the Arabia–Eurasia continental collision in the NW Zagros region. Altogether, the provenance shifts recorded in the pre-Zagros, proto-Zagros, and Neogene Zagros units, and in the WNK-derived RBS intermontane basin deposits reflect different phases of the Cenozoic Zagros orogeny.

Data availability

The new data used in this paper can be found online at (Koshnaw, 2021).


The supplement related to this article is available online at:

Author contributions

RK was responsible for conceptualizing the study, carrying out the investigation and data collection, writing original draft of the paper, reviewing and editing the paper, creating the visualizations, and acquiring funding. FS contributed to the investigation, resource acquisition, funding acquisition, and reviewing and editing the paper. DS was responsible for carrying out the investigation, data collection, and resource acquisition.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


This research was partially funded by the State Secretariat for Education, Research and Innovation of Switzerland via the Swiss Government Excellence Scholarship awarded to Renas Koshnaw. We would like to thank Idrees Ahmed, Varoujan Sissakian, Mazin Tamar-Agha, Waleed Shingaly, Nihad Karo, and Mjahid Zebari for helpful discussions and logistical support. We extend gratitude to Rudra Chatterjee, Desmond Patterson, and Lisa Stockli at The University of Texas at Austin UTChron Laboratories for their assistance. We thank Frédéric Mouthereau and an anonymous reviewer for their insightful reviews.

Financial support

This research has been partially supported by the State Secretariat for Education, Research and Innovation of Switzerland via the Swiss Government Excellence Scholarship.

Review statement

This paper was edited by Mark Allen and reviewed by Frédéric Mouthereau and one anonymous referee.


Abdallah, F. T. and Al-Dulaimi, S. I.: Biostratigraphy of the Upper Cretaceous for selected sections in northern Iraq, Iraqi Journal of Science, 60, 545–553, 2019. 

Abdel-Kireem, M. R.: Planktonic foraminifera and stratigraphy of the Tanjero Formation (Maastrichtian), northeastern Iraq, Micropaleontology, 32, 215–231, 1986. 

Abdollahi, F., Nabatian, G., Li, J. W., Honarmand, M., and Ebrahimi, M.: Saheb Granitoid Batholith, North of Kurdistan: An Evidence of Cretaceous-Paleocene Magmatism in the Sanandaj-Sirjan Zone, Journal Of Economic Geology, 12, 359–376, 2020. 

Abdula, R. A., Chicho, J., Surdashy, A., Nourmohammadi, M. S., Hamad, E., Muhammad, M. M., Smail, A. A., and Ashoor, A.: Sedimentology of the Govanda Formation at Gali Baza locality, Kurdistan region, Iraq, Iraqi Bulletin of Geology and Mining, 14, 1–12, 2018. 

Agard, P., Omrani, J., Jolivet, L., and Mouthereau, F.: Convergence history across Zagros (Iran): constraints from collisional and earlier deformation, Int. J. Earth Sci., 94, 401–419, 2005. 

Agard, P., Omrani, J., Jolivet, L., Whitechurch, H., Vrielynck, B., Spakman, W., Monié, P., Meyer, B., and Wortel, R.: Zagros orogeny: a subduction-dominated process, Geological Magazine, 148, 692–725, 2011. 

Al Humadi, H., Väisänen, M., Ismail, S. A., Kara, J., O'Brien, H., Lahaye, Y., and Lehtonen, M.: U–Pb geochronology and Hf isotope data from the Late Cretaceous Mawat ophiolite, NE Iraq, Heliyon, 5, e02721,, 2019. 

Al-Banna, N. Y. and Al-Mutwali, M. M.: Microfacies and age determination of the sedimentary sequences within Walash volcano-sedimentary Group, Mawat Complex, northeast Iraq, Tikrit Journal of Pure Science, 13, 46–52, 2008. 

Al-Bassam, K. S.: Petrology, geochemistry and depositional environment of the Khabour Formation in Ora and Khabour localities, Northern Iraq, Iraqi Bulletin of Geology and Mining, 6, 71–94, 2010. 

Al-Hadide, A. H., Al-Juboury, A. I., and Al-Eisa, M. I.: Stratigraphic sequence of late Paleozoic era in Iraq, Iraqi National Journal of Earth Sciences, 2, 28–42, 2002. 

Al-Hadidy, A. H.: Paleozoic stratigraphic lexicon and hydrocarbon habitat of Iraq, GeoArabia, 12, 63–130, 2007. 

Al-Husseini, M.: Preliminary chronostratigraphic chart of the Gulf Region, Gulf Petrolink, Manama, Bahrain, poster, 2000. 

Al-Juboury, A. I., Morton, A., Shingaly, W. S., Howard, J., Thusu, B., Vincent, S., and Fanning, M.: Stratigraphy and age revision of the Pirispiki Formation, Kurdistan Region, Northern Iraq, Arab. J. Geosci., 13, 1–15, 2020. 

Al-Mashaikie, S. Z. A. K., Al-Azzawi, A. T. A., and Kadum, A. K.: Depositional Environment of the Gercus Formation in Jabal Haibat Sultan, NE Iraq, New Sedimentological Approach, Iraqi Journal of Science, 55, 471–483, 2014. 

Al-Qayim, B. A., Al-Mutwali, M. M., and Nissan, B. Y.: Flysch–Molasse sediments of the Paleogene foreland basin of north Arabia, Shiranish area, north Iraq, Iraqi Bulletin of Geology and Mining, 4, 1–20, 2008. 

Al-Qayim, B. and Saadallah, A.: Sedimentary faces and evolution of an active margin buildup, Aqra-Bekhme Formation, N. E. Iraq, Iraqi Geological Journal, 27, 1–28, 1994. 

Al-Qayim, B., Omer, A., and Koyi, H.: Tectonostratigraphic overview of the Zagros suture zone, Kurdistan region, Northeast Iraq, GeoArabia, 17, 109–156, 2012. 

Al-Shaibani, S. and Al-Qayim, B.: Micropaleontology of Tertiary-Cretaceous transect, Dokan area, northeastern Iraq, Iraqi Journal of Science, 31, 353-396, 1990. 

Alavi, M.: Tectonics of the Zagros orogenic belt of Iran: new data and interpretations, Tectonophysics, 229, 211–238, 1994. 

Alavi, M.: Regional stratigraphy of the Zagros fold-thrust belt of Iran and its proforeland evolution, Am. J. Sci., 304, 1–20, 2004. 

Ali, S. A., Buckman, S., Aswad, K. J., Jones, B. G., Ismail, S. A., and Nutman, A. P.: Recognition of Late Cretaceous Hasanbag ophiolite-arc rocks in the Kurdistan Region of the Iraqi Zagros suture zone: A missing link in the paleogeography of the closing Neotethys Ocean, Lithosphere-US, 4, 395–410, 2012. 

Ali, S. A., Buckman, S., Aswad, K. J., Jones, B. G., Ismail, S. A., and Nutman, A. P.: The tectonic evolution of a Neo-T ethyan (Eocene-Oligocene) island-arc (Walash and Naopurdan groups) in the Kurdistan region of the Northeast Iraqi Zagros Suture Zone, Island Arc, 22, 104–125, 2013. 

Ali, S. A., Mohajjel, M., Aswad, K., Ismail, S., Buckman, S., and Jones, B.: Tectono-stratigraphy and structure of the northwestern Zagros collision zone across the Iraq-Iran border, Journal of Environment and Earth Science, 4, 92–110, 2014. 

Ali, S. A., Nutman, A. P., Aswad, K. J., and Jones, B. G.: Overview of the tectonic evolution of the Iraqi Zagros thrust zone: Sixty million years of Neotethyan ocean subduction, J. Geodyn., 129, 162–177, 2019. 

Ali, S. A., Sleabi, R. S., Talabani, M. J., and Jones, B. G.: Provenance of the Walash-Naopurdan back-arc-arc clastic sequences in the Iraqi Zagros Suture Zone, J. Afr. Earth Sci., 125, 73–87, 2017. 

Alirezaei, S. and Hassanzadeh, J.: Geochemistry and zircon geochronology of the Permian A-type Hasanrobat granite, Sanandaj–Sirjan belt: A new record of the Gondwana break-up in Iran, Lithos, 151, 122–134, 2012. 

Allen, M. B. and Armstrong, H. A.: Arabia-Eurasia collision and the forcing of mid-Cenozoic global cooling, Palaeogeogr. Palaeocl., 265, 52–58, 2008. 

Alsultan, H. A. A. and Gayara, D. A.: Basin Development of the Red Bed Series, Ne Iraq, Journal of University of Babylon, 24, 2, 435–555, 2016. 

Ameen, B. M.: Lthological indicators for the Oligocene unconformity, NE Iraq, Iraqi Bulletin of Geology and Mining, 5, 25–34, 2009. 

Ameen, M. S.: Effect of basement tectonics on hydrocarbon generation, migration, and accumulation in northern Iraq, AAPG Bull., 76, 356–370, 1992. 

Aqrawi, A. A., Goff, J. C., Horbury, A. D., and Sadooni, F. N.: The Petroleum Geology of Iraq, Scientific Press Ltd, Beaconsfield, UK, 424 pp., ISBN 978-0-901360-36-8, 2010. 

Aswad, K. J. and Elias, E. M.: Petrogenesis, geochemistry and metamorphism of spilitized subvolcanic rocks of the Mawat Ophiolite Complex, NE Iraq, Ofioliti, 13, 95–109, 1988. 

Aswad, K. J., Al-Samman, A. H., Aziz, N. R., and Koyi, A. M.: The geochronology and petrogenesis of Walash volcanic rocks, Mawat nappes: constraints on the evolution of the northwestern Zagros suture zone, Kurdistan Region, Iraq, Arab. J. Geosci., 7, 1403–1432, 2014. 

Aswad, K. J., Ali, S. A., Sheraefy, R. M. A., Nutman, A. P., Buckman, S., Jones, B. G., and Jourdan, F.: 40Ar/39Ar hornblende and biotite geochronology of the Bulfat igneous complex, Zagros suture zone, NE Iraq: new insights on complexities of Paleogene arc magmatism during closure of the Neotethys Ocean, Lithos, 266, 406–413, 2016. 

Avigad, D., Abbo, A., and Gerdes, A.: Origin of the Eastern Mediterranean: Neotethys rifting along a cryptic Cadomian suture with Afro-Arabia, Bulletin, 128, 1286–1296, 2016. 

Avigad, D., Morag, N., Abbo, A., and Gerdes, A.: Detrital rutile U-Pb perspective on the origin of the great Cambro-Ordovician sandstone of North Gondwana and its linkage to orogeny, Gondwana Res., 51, 17–29, 2017. 

Avouac, J. P. and Burov, E. B.: Erosion as a driving mechanism of intracontinental mountain growth, J. Geophys. Res.-Sol. Ea., 101, 17747–17769, 1996. 

Avouac, J. P., Meng, L., Wei, S., Wang, T., and Ampuero, J. P.: Lower edge of locked Main Himalayan Thrust unzipped by the 2015 Gorkha earthquake, Nat. Geosci., 8, 708–711, 2015. 

Aziz, N. R. H. and Sadiq, D. M.: U-Pb Zircon Dating of Upper Cretaceous Siliciclastic Rocks from the Tanjero Flysch, NE Iraq: New Constraints on their Provenance, and Tectonic Evolution, Kuwait J. Sci., 47, 4, 2020. 

Aziz, N. R., Elias, E. M., and Aswad, K. J.: Rb−Sr and Sm−Nd isotopes study of serpentinites and their impact on the tectonic setting of Zagros suture zone, NE Iraq, Iraqi Bulletin of Geology and Mining, 7, 67–75, 2011. 

Azizi, H., Tanaka, T., Asahara, Y., Chung, S. L., and Zarrinkoub, M. H.: Discrimination of the age and tectonic setting for magmatic rocks along the Zagros thrust zone, northwest Iran, using the zircon U-Pb age and Sr−Nd isotopes, J. Geodyn., 52, 304–320, 2011. 

Azizi, H., Hadi, A., Asahara, Y., and Mohammad, Y.: Geochemistry and geodynamics of the Mawat mafic complex in the Zagros Suture zone, northeast Iraq, Open Geosci., 5, 523–537, 2013. 

Azizi, H., Hadad, S., Stern, R. J., and Asahara, Y.: Age, geochemistry, and emplacement of the  40-Ma Baneh granite–appinite complex in a transpressional tectonic regime, Zagros suture zone, northwest Iran, Int. Geol. Rev., 61, 195–223, 2019. 

Ballato, P., Uba, C. E., Landgraf, A., Strecker, M. R., Sudo, M., Stockli, D. F., Friedrich, A., and Tabatabaei, S. H.: Arabia-Eurasia continental collision: Insights from late Tertiary foreland-basin evolution in the Alborz Mountains, northern Iran, Bulletin, 123, 106–131, 2011. 

Barber, D. E., Stockli, D. F., Horton, B. K., and Koshnaw, R. I.: Cenozoic exhumation and foreland basin evolution of the Zagros orogen during the Arabia-Eurasia collision, western Iran, Tectonics, 37, 4396–4420, 2018. 

Barber, D. E., Stockli, D. F., and Galster, F.: The Proto-Zagros foreland basin in Lorestan, western Iran: Insights from multimineral detrital geothermochronometric and trace elemental provenance analysis, Geochem. Geophy. Geosy., 20, 2657–2680, 2019. 

Barrier, E. and Vrielynck, B.: MEBE atlas of paleotectonic maps of the Middle East, Commission for the geological map of the world, Paris, France, 2008. 

Behyari, M., Mohajjel, M., Sobel, E. R., Rezaeian, M., Moayyed, M., and Schmidt, A.: Analysis of exhumation history in Misho Mountains, NW Iran: Insights from structural and apatite fission track data, Neues Jahrb. Geol. P.-A., 283, 291–308, 2017. 

Braud, J.: La Suture Du Zagros Au Niveau De Kermanshah (Kurdistan Iranien): Reconstitution Paleogeographique, Evolution Geodynamique, Magmatique Et Structurale, PhD Thesis, Universtie Paris-Sud, pp. 488, unpublished, 1987. 

Braud, J. E. A. N. and Ricou, L. E.: Eléments de continuité entre le Zagros et la Turquie du Sud-Est, B. Soc. Geol. Fr., 7, 1015–1023, 1975. 

Cawood, P. A., Hawkesworth, C. J., and Dhuime, B.: Detrital zircon record and tectonic setting, Geology, 40, 875–878, 2012. 

Çelik, H. and Salih, T. M. H.: Provenance investigation from sedimentary petrography of the Upper Cretaceous deep marine low density turbidites of the Tanjero Formation around Arbat, northeastern Iraq, Turk. J. Earth Sci., 27, 432–459, 2018. 

Çelik, H. and Salih, T. M. H.: Petrographic characteristics of deep marine turbidite sandstones of the Upper Cretaceous Tanjero Formation, Northwestern Sulaimaniyah, Iraq: implications for provenance and tectonic setting, Bulletin of the Mineral Research and Exploration, 164, 11–38, 2021. 

Chiu, H. Y., Chung, S. L., Zarrinkoub, M. H., Mohammadi, S. S., Khatib, M. M., and Iizuka, Y.: Zircon U-Pb age constraints from Iran on the magmatic evolution related to Neotethyan subduction and Zagros orogeny, Lithos, 162, 70–87, 2013. 

Cloos, M.: Lithospheric buoyancy and collisional orogenesis: Subduction of oceanic plateaus, continental margins, island arcs, spreading ridges, and seamounts, Geol. Soc. Am. Bull., 105, 715–737, 1993. 

Colleps, C. L., McKenzie, N. R., Horton, B. K., Webb, A. A. G., Ng, Y. W., and Singh, B. P.: Sediment provenance of pre-and post-collisional Cretaceous-Paleogene strata from the frontal Himalaya of northwest India, Earth Planet. Sc. Lett., 534, 116079, 2020. 

Darin, M. H., Umhoefer, P. J., and Thomson, S. N.: Rapid late Eocene exhumation of the Sivas Basin (Central Anatolia) driven by initial Arabia-Eurasia collision, Tectonics, 37, 3805–3833, 2018. 

Delaloye, M. and Desmons, J.: Ophiolites and mélange terranes in Iran: a geochronological study and its paleotectonic implications, Tectonophysics, 68, 83–111, 1980. 

Dercourt, J., Zonenshain, L. P., Ricou, L. E., Kazmin, V. G., Le Pichon, X., Knipper, A. L., Grandjacquet, C., Sbortshikov, I. M., Geyssant, J., Lepvrier, C., and Pechersky, D. H.: Geological evolution of the Tethys belt from the Atlantic to the Pamirs since the Lias, Tectonophysics, 123, 241–315, 1986. 

Dewey, J. F., Pitman III, W. C., Ryan, W. B., and Bonnin, J.: Plate tectonics and the evolution of the Alpine system, Geol. Soc. Am. Bull., 84, 3137–3180, 1973. 

Dunnington, H. V.: Generation, migration, accumulation and dissipation of oil in northern Iraq, in: Habitat of oil, a symposium, edited by: Weeks, L. G., Am. Assoc. Petr. Geol., Tulsa, 1149–1251, 1958. 

English, J. M., Lunn, G. A., Ferreira, L., and Yacu, G.: Geologic evolution of the Iraqi Zagros, and its influence on the distribution of hydrocarbons in the Kurdistan region, AAPG Bull., 99, 231–272, 2015. 

Etemad-Saeed, N., Najafi, M., and Vergés, J.: Provenance evolution of Oligocene-Pliocene foreland deposits in the Dezful embayment to constrain Central Zagros exhumation history, J. Geol. Soc., 177, 799–817, 2020. 

Farley, K. A., Wolf, R. A., and Silver, L. T.: The effects of long alpha-stopping distances on (U-Th) / He ages, Geochim. Cosmochim. Ac., 60, 4223–4229, 1996. 

Francois, T., Agard, P., Bernet, M., Meyer, B., Chung, S. L., Zarrinkoub, M. H., Burov, E., and Monie, P.: Cenozoic exhumation of the internal Zagros: first constraints from low-temperature thermochronology and implications for the build-up of the Iranian plateau, Lithos, 206, 100–112, 2014. 

Gehrels, G.: Detrital zircon U-Pb geochronology applied to tectonics, Annu. Rev. Earth Pl. Sc., 42, 127–149, 2014. 

Gehrels, G. E.: Introduction to detrital zircon studies of Paleozoic and Triassic strata in western Nevada and northern California, Special Paper of the Geological Society of America, 347, 1–17, 2000. 

Golan, T., Katzir, Y., and Coble, M. A.: Early Carboniferous anorogenic magmatism in the Levant: implications for rifting in northern Gondwana, Int. Geol. Rev., 60, 101–108, 2018. 

Grimes, C. B., John, B. E., Kelemen, P. B., Mazdab, F. K., Wooden, J. L., Cheadle, M. J., Hanghøj, K., and Schwartz, J. J.: Trace element chemistry of zircons from oceanic crust: A method for distinguishing detrital zircon provenance, Geology, 35, 643–646, 2007. 

Hart, N. R., Stockli, D. F., and Hayman, N. W.: Provenance evolution during progressive rifting and hyperextension using bedrock and detrital zircon U-Pb geochronology, Mauléon Basin, western Pyrenees, Geosphere, 12, 1166–1186, 2016. 

Hassan, M. M., Jones, B. G., Buckman, S., Al-Jubory, A. I., and Al Gahtani, F. M.: Provenance of Paleocene-Eocene red beds from NE Iraq: constraints from framework petrography, Geol. Mag., 151, 1034–1050, 2014. 

Hassanzadeh, J. and Wernicke, B. P.: The Neotethyan Sanandaj-Sirjan zone of Iran as an archetype for passive margin-arc transitions, Tectonics, 35, 586–621, 2016. 

Hempton, M. R.: Structure and deformation history of the Bitlis suture near Lake Hazar, southeastern Turkey, Geol. Soc. Am. Bull., 96, 233–243, 1985. 

Hempton, M. R.: Constraints on Arabian plate motion and extensional history of the Red Sea, Tectonics, 6, 687–705, 1987. 

Hessami, K., Koyi, H. A., Talbot, C. J., Tabasi, H., and Shabanian, E.: Progressive unconformities within an evolving foreland fold–thrust belt, Zagros Mountains, J. Geol. Soc., 158, 969–981, 2001. 

Homke, S., Vergés, J., Van Der Beek, P., Fernàndez, M., Saura, E., Barbero, L., Badics, B., and Labrin, E.: Insights in the exhumation history of the NW Zagros from bedrock and detrital apatite fission-track analysis: evidence for a long-lived orogeny, Basin Res., 22, 659–680, 2010. 

Horton, B. K., Hassanzadeh, J., Stockli, D. F., Axen, G. J., Gillis, R. J., Guest, B., Amini, A., Fakhari, M. D., Zamanzadeh, S. M., and Grove, M.: Detrital zircon provenance of Neoproterozoic to Cenozoic deposits in Iran: Implications for chronostratigraphy and collisional tectonics, Tectonophysics, 451, 97–122, 2008. 

Horton, B. K.: Cenozoic evolution of hinterland basins in the Andes and Tibet, in: Tectonics of Sedimentary Basins: Recent Advances, first edn., edited by: Busby, C. and Azor, A., Blackwell Publishing Ltd., Oxford, UK, 2012. 

Ismail, S. A., Koshnaw, R. I., Barber, D. E., Al Humadi, H., and Stockli, D. F.: Generation and exhumation of granitoid intrusions in the Penjween ophiolite complex, NW Zagros of the Kurdistan region of Iraq: Implications for the geodynamic evolution of the Arabian-Eurasian collision zone, Lithos, 376, 105714,, 2020. 

Jagoutz, O., Royden, L., Holt, A. F., and Becker, T. W.: Anomalously fast convergence of India and Eurasia caused by double subduction, Nat. Geosci., 8, 475–478, 2015. 

James, G. A. and Wynd, J. G.: Stratigraphic nomenclature of Iranian oil consortium agreement area, AAPG Bull., 49, 2182–2245, 1965. 

Jassim, S. Z. and Sisakian, V.: Field guide to the geology of Salah Al-Din-Shaqlawa area, northeastern Iraq, Geological society of Iraq and Union of Iraqi Geologists, Fifth Iraqi geological congress, Baghdad, Iraq, 1978. 

Jassim, S. Z. and Buday, T.: Units of the unstable shelf and the Zagros suture, in: Geology of Iraq, edited by: Jassim, S. Z. and Goff, J. C., Dolin and Moravian Museum, Prague and Brno, Czech Republic, 73–90, 2006. 

Jassim, S. Z., Buday, T., Cicha, I., and Opletal, M.: Tectonostratigraphy of the Zagros suture, in: Geology of Iraq, edited by: Jassim, S. Z. and Goff, J. C., Dolin and Moravian Museum, Prague and Brno, Czech Republic, 276–29, 2006. 

Jones, B. G., Ali, S. A., and Nutman, A. P.: Provenance of Tanjero and Red Bed clastic sedimentary rocks revealed by detrital zircon SHRIMP dating, Kurdistan region, NE Iraq: Constraints on ocean closure and unroofing of Neo-Tethyan allochthons, J. Afr. Earth Sci., 172, 103981, 2020. 

Kadem, L. S.: Mixed carbonate clastic sediment facies of the Paleocene-L, Eocene Kolosh Formation from selected locations in North Iraq, Tikrit Journal of Pure Science, 11, 205–209, (last access: 29 October 2021) 2006. 

Karim, H. K. and Surdashy, A. M.: Paleocurrent Analysis of Upper Cretaceous Zagros Foreland Basin: A Case Study for Tanjero Formation in Sulaimaniya Area NE-Iraq, Iraqi National Journal of Earth Sciences, 5, 30–44, 2005. 

Karim, K. H., Al-Barzinjy, S. T., and Ameen, B. M.: History and geological setting of intermontane basins in the Zagros fold-thrust belt, Kurdistan region, NE Iraq, Iraqi Bulletin of Geology and Mining, 4, 21–33, 2008. 

Karim, K. H., Koyi, H., Baziany, M. M., and Hessami, K.: Significance of angular unconformities between Cretaceous and Tertiary strata in the northwestern segment of the Zagros fold-thrust belt, Kurdistan Region, NE Iraq, Geol. Mag., 148, 925–939, 2011. 

Karo, N. M., Oberhänsli, R., Aqrawi, A. M., Elias, E. M., Aswad, K. J., and Sudo, M.: New 40Ar/39Ar age constraints on cooling and unroofing history of the metamorphic host rocks (and igneous intrusion associates) from the Bulfat Complex (Bulfat area), NE-Iraq, Arab. J. Geosci., 11, 1–11, 2018. 

Khadivi, S., Mouthereau, F., Larrasoaña, J. C., Vergés, J., Lacombe, O., Khademi, E., and Suc, J. P.: Magnetochronology of synorogenic Miocene foreland sediments in the Fars arc of the Zagros Folded Belt (SW Iran), Basin Res., 22, 918–932, 2010. 

Khadivi, S., Mouthereau, F., Barbarand, J., Adatte, T., and Lacombe, O.: Constraints on palaeodrainage evolution induced by uplift and exhumation on the southern flank of the Zagros-Iranian Plateau, J. Geol. Soc., 169, 83–97, 2012. 

Kharajiany, S., Al-Qayim, B., and Wise Jr., S.: Calcareous nannofossil stratigraphy of the Upper Cretaceous-lower Paleocene sequence from the Chinarok section, Sulaimaniah area, Kurdistan region, NE Iraq, Iraqi Bulletin of Geology and Mining, 15, 1–4, 2019. 

Konert, G., Afifi, A. M., Al-Hajri, S. I. A., and Droste, H. J.: Paleozoic stratigraphy and hydrocarbon habitat of the Arabian Plate, GeoArabia, 6, 407–442, 2001. 

Koop, W. J. and Stoneley, R.: Subsidence history of the Middle East Zagros basin, Permian to recent, Philosophical Transactions of the Royal Society of London, Series A, Mathematical and Physical Sciences, 305, 149–168, 1982. 

Koshnaw, R.: Supporting_information_Koshnaw_etal_SolidEarth, V1, Mendeley Data [data set],, 2021. 

Koshnaw, R. I., Horton, B. K., Stockli, D. F., Barber, D. E., Tamar-Agha, M. Y., and Kendall, J. J.: Neogene shortening and exhumation of the Zagros fold-thrust belt and foreland basin in the Kurdistan region of northern Iraq, Tectonophysics, 694, 332–355, 2017. 

Koshnaw, R. I., Stockli, D. F., and Schlunegger, F.: Timing of the Arabia-Eurasia continental collision-Evidence from detrital zircon U-Pb geochronology of the Red Bed Series strata of the northwest Zagros hinterland, Kurdistan region of Iraq, Geology, 47, 47–50, 2019. 

Koshnaw, R. I., Horton, B. K., Stockli, D. F., Barber, D. E., and Tamar-Agha, M. Y.: Sediment routing in the Zagros foreland basin: Drainage reorganization and a shift from axial to transverse sediment dispersal in the Kurdistan region of Iraq, Basin Res., 32, 688–715, 2020a. 

Koshnaw, R. I., Stockli, D. F., Horton, B. K., Teixell, A., Barber, D. E., and Kendall, J. J.: Late Miocene Deformation Kinematics Along the NW Zagros Fold-Thrust Belt, Kurdistan Region of Iraq: Constraints From Apatite (U-Th) / He Thermochronometry and Balanced Cross Sections, Tectonics, 39, e2019TC005865, 2020b. 

Koyi, A. M.: Sr−Nd isotopical significance of Walash volcanic rocks, Mawat area, NE Iraq, Zanco J. Pure Appl. Sci. 21, 39–45, 2009. 

Lawa, F. A., Koyi, H., and Ibrahim, A.: Tectono-stratigraphic evolution of the NW segment OF the Zagros fold-thrust belt, Kurdistan, NE Iraq, J. Petrol. Geol., 36, 75–96, 2013. 

Le Garzic, E., Vergés, J., Sapin, F., Saura, E., Meresse, F., and Ringenbach, J. C.: Evolution of the NW Zagros Fold-and-Thrust Belt in Kurdistan Region of Iraq from balanced and restored crustal-scale sections and forward modelling, J. Struct. Geol., 124, 51–69, 2019. 

Leterrier, J.: Mineralogical, geochemical and isotopic evolution of two Miocene mafic intrusions from the Zagros (Iran), Lithos, 18, 311–329, 1985. 

Marsh, J. H. and Stockli, D. F.: Zircon U-Pb and trace element zoning characteristics in an anatectic granulite domain: Insights from LASS-ICP-MS depth profiling, Lithos, 239, 170–185, 2015. 

Mazhari, S. A., Bea, F., Amini, S., Ghalamghash, J., Molina, J. F., Montero, P., Scarrow, J. H., and Williams, I. S.: The Eocene bimodal Piranshahr massif of the Sanandaj-Sirjan Zone, NW Iran: a marker of the end of the collision in the Zagros orogeny, J. Geol. Soc., 166, 53–69, 2009. 

Mazhari, S. A., Ghalamghash, J., Kumar, S., Shellnutt, J. G., and Bea, F.: Tectonomagmatic development of the Eocene Pasevh pluton (NW Iran): Implications for the Arabia-Eurasia collision, J. Asian Earth Sci., 203, 104551, 2020. 

McQuarrie, N. and van Hinsbergen, D. J.: Retrodeforming the Arabia-Eurasia collision zone: Age of collision versus magnitude of continental subduction, Geology, 41, 315–318, 2013. 

McQuarrie, N., Stock, J. M., Verdel, C., and Wernicke, B. P.: Cenozoic evolution of Neotethys and implications for the causes of plate motions, Geophys. Res. Lett., 30, 2036,, 2003. 

Meinhold, G., Bassis, A., Hinderer, M., Lewin, A., and Berndt, J.: Detrital zircon provenance of north Gondwana Palaeozoic sandstones from Saudi Arabia, Geological Magazine, 158, 442–458, 2021. 

Minstry of Natural Resources (MNR), Kurdistan regional government.: Occurrences of metallic deposits in the kurdistan region-Iraq, Kurdistan region, Iraq, 2016. 

Moghadam, H. S. and Stern, R. J.: Geodynamic evolution of Upper Cretaceous Zagros ophiolites: formation of oceanic lithosphere above a nascent subduction zone, Geological Magazine, 148, 762–801, 2011. 

Moghadam, H. S., Li, Q. L., Stern, R. J., Chiaradia, M., Karsli, O., and Rahimzadeh, B.: The Paleogene ophiolite conundrum of the Iran-Iraq border region, J. Geol. Soc., 177, 955–964, 2020. 

Mohammad, Y. O. and Cornell, D. H.: U-Pb zircon geochronology of the Daraban leucogranite, Mawat ophiolite, Northeastern Iraq: a record of the subduction to collision history for the Arabia-Eurasia plates, Isl. Arc, 26, e12188, 2017. 

Mohammad, Y. O., Cornell, D. H., Qaradaghi, J. H., and Mohammad, F. O.: Geochemistry and Ar−Ar muscovite ages of the Daraban Leucogranite, Mawat Ophiolite, northeastern Iraq: implications for Arabia-Eurasia continental collision, J. Asian Earth Sci., 86, 151–165, 2014. 

Molnar, P. and England, P.: Late Cenozoic uplift of mountain ranges and global climate change: chicken or egg?, Nature, 346, 29–34, 1990. 

Moritz, R., Ghazban, F., and Singer, B. S.: Eocene gold ore formation at Muteh, Sanandaj-Sirjan tectonic zone, Western Iran: a result of late-stage extension and exhumation of metamorphic basement rocks within the Zagros Orogen, Econ. Geol., 101, 1497–1524, 2006. 

Mouthereau, F., Lacombe, O., and Vergés, J.: Building the Zagros collisional orogen: timing, strain distribution and the dynamics of Arabia/Eurasia plate convergence, Tectonophysics, 532, 27–60, 2012. 

Müller, R. D., Cannon, J., Qin, X., Watson, R. J., Gurnis, M., Williams, S., Pfaffelmoser, T., Seton, M., Russell, S. H., and Zahirovic, S.: GPlates: building a virtual Earth through deep time, Geochem. Geophy. Geosy., 19, 2243–2261, 2018. 

Najman, Y., Appel, E., Boudagher-Fadel, M., Bown, P., Carter, A., Garzanti, E., Godin, L., Han, J., Liebke, U., Oliver, G., and Parrish, R.: Timing of India-Asia collision: Geological, biostratigraphic, and palaeomagnetic constraints, J. Geophys. Res.-Sol. Ea., 115, , 2010. 

Navabpour, P., Barrier, E., and McQuillan, H.: Oblique oceanic opening and passive margin irregularity, as inherited in the Zagros fold-and-thrust belt, Terra Nova, 26, 208–215, 2014. 

Norton, K. and Schlunegger, F.: Migrating deformation in the Central Andes from enhanced orographic rainfall, Nat. Commun., 2, 1–7, 2011. 

Nouri, F., Azizi, H., Golonka, J., Asahara, Y., Orihashi, Y., Yamamoto, K., Tsuboi, M., and Anma, R.: Age and petrogenesis of Na-rich felsic rocks in western Iran: evidence for closure of the southern branch of the Neo-Tethys in the Late Cretaceous, Tectonophysics, 671, 151–172, 2016. 

Numan, N. M.: A plate tectonic scenario for the Phanerozoic succession in Iraq, Iraqi Geological Journal, 30, 85–110, 1997. 

Numan, N. M. S., Hammoudi, R. A., and Chorowicz, J.: Synsedimentary tectonics in the Eocene Pila Spi limestone formation in Iraq and its geodynamic implications, J. Afr. Earth Sci., 27, 141–148, 1998. 

Oberhänsli, R. O. R. G. A. J. M., Candan, O., Bousquet, R., Rimmele, G., Okay, A., and Goff, J.: Alpine high pressure evolution of the eastern Bitlis complex, SE Turkey, Special Publications, London, Geol. Soc., 340, 461–483, 2010. 

Omrani, J., Agard, P., Whitechurch, H., Benoit, M., Prouteau, G., and Jolivet, L.: Arc-magmatism and subduction history beneath the Zagros Mountains, Iran: a new report of adakites and geodynamic consequences, Lithos, 106, 380–398, 2008. 

Orme, D. A., Carrapa, B., and Kapp, P.: Sedimentology, provenance and geochronology of the Upper Cretaceous-Lower Eocene western Xigaze forearc basin, southern Tibet, Basin Res., 27, 387–411, 2015. 

Omer, M. F., Friis, H., Kokfelt, T. F., and Thomsen, T. B.: Provenance of northern Gondwana Lower Ordovician sandstones (Khabour Formation, northern Iraq) revealed by detrital zircon using LA-ICP-MS dating, Geol. J., 56, 1–18,, 2021. 

Pujols, E. J., Stockli, D. F., Constenius, K. N., and Horton, B. K.: Thermochronological and geochronological constraints on Late Cretaceous unroofing and proximal sedimentation in the Sevier orogenic belt, Utah, Tectonics, 39, e2019TC005794, 2020. 

Reiners, P. W., Farley, K. A., and Hickes, H. J.: He diffusion and (U-Th) / He thermochronometry of zircon: initial results from Fish Canyon Tuff and Gold Butte, Tectonophysics, 349, 297–308, 2002. 

Robertson, A. H. F., Parlak, O., Rizaoğlu, T., Ünlügenç, Ü., İnan, N., Tasli, K., and Ustaömer, T.: Tectonic evolution of the South Tethyan ocean: evidence from the Eastern Taurus Mountains (Elaziğ region, SE Turkey), Special Publications, London, Geol. Soc., 272, 231–270, 2007. 

Rolland, Y., Perincek, D., Kaymakci, N., Sosson, M., Barrier, E., and Avagyan, A.: Evidence for 80–75 Ma subduction jump during Anatolide-Tauride-Armenian block accretion and 48 Ma Arabia-Eurasia collision in Lesser Caucasus-East Anatolia, J. Geodyn., 56, 76–85, 2012. 

Salih, N., Mansurbeg, H., Kolo, K., Gerdes, A., and Préat, A.: In situ U-Pb dating of hydrothermal diagenesis in tectonically controlled fracturing in the Upper Cretaceous Bekhme Formation, Kurdistan Region-Iraq, Int. Geol. Rev., 62, 2261–2279, 2020. 

Sasvari, A., Davies, L., Mann, A., Afzal, J., Vakarcs, G., and Iwaniw, E.: Dachstein-type Avroman Formation: An indicator of the Harsin Basin in Iraq, GeoArabia, 20, 17–36, 2015. 

Saura, E., Vergés, J., Homke, S., Blanc, E., Serra-Kiel, J., Bernaola, G., Casciello, E., Fernández, N., Romaire, I., Casini, G., and Embry, J. C.: Basin architecture and growth folding of the NW Zagros early foreland basin during the Late Cretaceous and early Tertiary, J. Geol. Soc., 168, 235–250, 2011. 

Saura, E., Garcia-Castellanos, D., Casciello, E., Parravano, V., Urruela, A., and Vergés, J.: Modeling the flexural evolution of the Amiran and Mesopotamian foreland basins of NW Zagros (Iran-Iraq), Tectonics, 34, 377–395, 2015. 

Saylor, J. E. and Sundell, K. E.: Quantifying comparison of large detrital geochronology data sets, Geosphere, 12, 203–220, 2016. 

Sharland, P. R., Archer, R., Casey, D. M., Davies, R. B., Hall, S. H., and Heward, A. P.: Arabian plate sequence stratigraphy, GeoArabia Special Publication 2, Gulf PetroLink, Manama, Bahrain, 371 pp. 2001. 

Shawkat, M. G. and Tucker, M. E.: Stromatolites and sabkha cycles from the lower Fars formation (Miocene) of Iraq, Geol. Rev., 67, 1–14, 1978. 

Sissakian, V. K.: Geological Map of Arbeel and Mahabad Quadrangles Sheets NJ-38-14 and NJ-38-15, Geological Survey and Mining, Baghdad, Iraq, Scale 1:250 000, 1997. 

Stampfli, G. M. and Borel, G. D.: A plate tectonic model for the Paleozoic and Mesozoic constrained by dynamic plate boundaries and restored synthetic oceanic isochrones, Earth Planet. Sc. Lett., 196, 17–33, 2002. 

Stern, R. J.: Subduction zones, Rev. Geophy., 40, 3–1, 2002. 

Stern, R. J.: The anatomy and ontogeny of modern intra-oceanic arc systems, Special Publications, London, Geol. Soc., 338, 7–34, 2010. 

Stern, R. J. and Johnson, P.: Continental lithosphere of the Arabian Plate: a geologic, petrologic, and geophysical synthesis, Earth-Sci. Rev., 101, 29–67, 2010. 

Stern, R. J., Ren, M., Ali, K., Förster, H. J., Al Safarjalani, A., Nasir, S., and Romer, R. L.: Early Carboniferous (357 Ma) crust beneath northern Arabia: Tales from Tell Thannoun (southern Syria), Earth Planet. Sci. Lett., 393, 83–93, 2014. 

Stockli, D., Boyd, P., and Galster, F.: Intragrain common Pb correction in apatite by LA-ICP-MS depth profiling and implications for detrital apatite U-Pb dating, Geophysical Research Abstracts, 19th EGU General Assembly, Vienna, Austria, 23–28 April 2017, p. 12225, 2017. 

Stoneley, R.: On the origin of ophiolite complexes in the southern Tethys region, Tectonophysics, 25, 303–322, 1975. 

Tamar-Agha, M. Y. and Al-Aslami, O. J. M.: Facies, Depositional Environment and Cyclicity of the Fatha Formation in East Baghdad Oil Field, Iraq, Iraqi Journal of Science, 56, 2939–2952, 2015. 

Tamar-Agha, M. Y. and Salman, N. A.: Facies and Depositional Environments of Injana Formation in Zawita, Amadia and Zakho Areas, Northern Iraq, Iraqi Bulletin of Geology and Mining, 11, 39–60, 2015. 

Tamar-Agha, M., Numan, N. M. S., and Malala, K.: Field guide for the geology of Dohok and Sinjar areas, Geological society of Iraq and Union of Iraqi Geologists, Fifth Iraqi Geological Congress, Baghdad, Iraq, 1978. 

Thomson, K. D., Stockli, D. F., Clark, J. D., Puigdefàbregas, C., and Fildani, A.: Detrital zircon (U-Th) / He double-dating constraints on provenance and foreland basin evolution of the Ainsa Basin, south-central Pyrenees, Spain, Tectonics, 36, 1352–1375, 2017. 

van Bellen, V. R. C., Dunnington, H. V., Wetzel, R., and Morton, D. M.: Lexique Stratigraphique International 3, Asie, fasc. 10a, Iraq, CNRS (Cent. Natl. Rech. Sci.), Paris, 1959. 

van der Boon, A., Kuiper, K. F., van der Ploeg, R., Cramwinckel, M. J., Honarmand, M., Sluijs, A., and Krijgsman, W.: Exploring a link between the Middle Eocene Climatic Optimum and Neotethys continental arc flare-up, Clim. Past, 17, 229–239,, 2021. 

van der Meer, D. G., Van Hinsbergen, D. J., and Spakman, W.: Atlas of the underworld: Slab remnants in the mantle, their sinking history, and a new outlook on lower mantle viscosity, Tectonophysics, 723, 309–448, 2018. 

van Hinsbergen, D. J., Torsvik, T. H., Schmid, S. M., Maţenco, L. C., Maffione, M., Vissers, R. L., and Spakman, W.: Orogenic architecture of the Mediterranean region and kinematic reconstruction of its tectonic evolution since the Triassic, Gondwana Res., 81, 79–229, 2020. 

Verdel, C., Wernicke, B. P., Hassanzadeh, J., and Guest, B.: A Paleogene extensional arc flare-up in Iran, Tectonics, 30, TC3008,, 2011. 

Vergés, J., Saura, E., Casciello, E., Fernandez, M., Villaseñor, A., Jimenez-Munt, I., and García-Castellanos, D.: Crustal-scale cross-sections across the NW Zagros belt: implications for the Arabian margin reconstruction, Geol. Mag., 148, 739–761, 2011. 

Vergés, J. H. M. E. S. P. A., Emami, H., Garcés, M., Beamud, E., Homke, S., and Skott, P.: Zagros foreland fold belt timing across Lurestan to constrain Arabia-Iran collision, in: Developments in Structural Geology and Tectonics, Elsevier, 3, 29–52, 2019. 

Vermeesch, P.: How many grains are needed for a provenance study?, Earth Planet. Sc. Lett., 224, 441–451, 2004. 

Vermeesch, P.: On the visualisation of detrital age distributions, Chem. Geol., 312, 190–194,, 2012.  

Vermeesch, P.: Multi-sample comparison of detrital age distributions, Chem. Geol., 341, 140–146, 2013. 

Webb, A. A. G., Yin, A., and Dubey, C. S.: U-Pb zircon geochronology of major lithologic units in the eastern Himalaya: Implications for the origin and assembly of Himalayan rocks, Geol. Soc. Am. Bull., 125, 499–522, 2013. 

Wells, M., Morton, A., and Frei, D.: Provenance of Lower Cretaceous clastic reservoirs in the Middle East, J. Geol. Soc., 174, 1048–1061, 2017. 

Whitechurch, H., Omrani, J., Agard, P., Humbert, F., Montigny, R., and Jolivet, L.: Evidence for Paleocene-Eocene evolution of the foot of the Eurasian margin (Kermanshah ophiolite, SW Iran) from back-arc to arc: implications for regional geodynamics and obduction, Lithos, 182, 11–32, 2013. 

Wolfe, M. R. and Stockli, D. F.: Zircon (U-Th) / He thermochronometry in the KTB drill hole, Germany, and its implications for bulk He diffusion kinetics in zircon, Earth Planet. Sci. Lett., 295, 69–82, 2010. 

Wrobel-Daveau, J. C., Ringenbach, J. C., Tavakoli, S., Ruiz, G. M., Masse, P., and de Lamotte, D. F.: Evidence for mantle exhumation along the Arabian margin in the Zagros (Kermanshah area, Iran), Arab. J. Geosci., 3, 499–513, 2010. 

Xu, J., Stockli, D. F. and Snedden, J. W.: Enhanced provenance interpretation using combined U-Pb and (U-Th) / He double dating of detrital zircon grains from lower Miocene strata, proximal Gulf of Mexico Basin, North America, Earth Planet. Sci. Lett., 475, 44–57, 2017. 

Yilmaz, Y.: New evidence and model on the evolution of the southeast Anatolian orogen, Geol. Soc. Am. Bull., 105, 251–271, 1993. 

Zhang, Z., Xiao, W., Majidifard, M. R., Zhu, R., Wan, B., Ao, S., Chen, L., Rezaeian, M., and Esmaeili, R.: Detrital zircon provenance analysis in the Zagros Orogen, SW Iran: implications for the amalgamation history of the Neo-Tethys, Int. J. Earth Sci., 106, 1223–1238, 2017. 

Ziegler, M. A.: Late Permian to Holocene paleofacies evolution of the Arabian Plate and its hydrocarbon occurrences, GeoArabia, 6, 445–504, 2001. 

Short summary
As continental plates collide, mountain belts grow. This study investigated the provenance of rocks from the northwestern segment of the Zagros mountain belt to unravel the convergence history of the Arabian and Eurasian plates. Provenance data synthesis and field relationships suggest that the Zagros Mountains developed as a result of the oceanic crust emplacement on the Arabian continental plate, followed by the Arabia–Eurasia collision and later uplift of the broader region.