No abstract

Cenozoic exhumation patterns in the internal and external Zagros reveal a long‐term deformation record associated with geodynamic restructuring of Arabia‐Eurasia collisional zone from continental subduction to plate suturing, which can be evaluated from thermochronometric, provenance, and subsidence analyses. Thermal modeling of zircon and apatite (U‐Th)/He ages and apatite fission track data from the Sanandaj‐Sirjan Zone (SSZ) indicates exhumation and inferred uplift along the leading edge of Eurasia starting in the Late Eocene (~35 Ma), coeval with initial foreland flexural subsidence of Arabia. Together with deceleration in Arabia‐Eurasia convergence and diminished subduction‐related magmatism, these events signal the final Neotethys closure and onset of long‐term (15–20 Myr) Arabian continental subduction beneath Eurasia, facilitated by the attenuated architecture of the precollisional Arabian margin. From 35 to 20 Ma, crustal shortening was relatively subdued and restricted to areas along the Arabia‐Eurasia plate boundary and diffuse inversion structures within continental interiors. Acceleration in SSZ cooling/exhumation rates from 19 to 16 Ma was synchronous with rapid basin subsidence and clastic progradation in the Zagros foreland. These events were contemporaneous with 20‐ to 16‐Ma surge in calc‐alkaline magmatism in central Iran and may have been linked to reorganization/deflection of Arabian plate vectors during the main phase of Red Sea rifting at 19–18 Ma. Transition from continental subduction to Arabia‐Eurasia suturing by ~12 Ma forced a transfer of strain from the subduction zone to intraplate deformational structures. This was marked by rapid outward expansion of the Zagros orogen, involving a shift in exhumation from the SSZ to Zagros fold‐thrust belt and Iranian plate interior.

In the northwestern sector of the Zagros foreland basin, axial fluvial systems initially delivered fine‐grained sediments from northwestern source regions into a contiguous basin, and later transverse fluvial systems delivered coarse‐grained sediments from northeastern sources into a structurally partitioned basin by fold‐thrust deformation. Here we integrate sedimentologic, stratigraphic, palaeomagnetic and geochronologic data from the northwestern Zagros foreland basin to define the Neogene history of deposition and sediment routing in response to progressive advance of the Zagros fold‐thrust belt. This study constrains the depositional environments, timing of deposition and provenance of nonmarine clastic deposits of the Injana (Upper Fars), Mukdadiya (Lower Bakhtiari) and Bai‐Hasan (Upper Bakhtiari) Formations in the Kurdistan region of Iraq. Sediments of the Injana Formation (~12.4–7.75 Ma) were transported axially (orogen‐parallel) from northwest to southeast by meandering and low‐sinuosity channel belt system. In contrast, during deposition of the Mukdadiya Formation (~7.75–5 Ma), sediments were delivered transversely (orogen‐perpendicular) from northeast to southwest by braided and low‐sinuosity channel belt system in distributive fluvial megafans. By ~5 Ma, the northwestern Zagros foreland basin became partitioned by growth of the Mountain Front Flexure and considerable gravel was introduced in localized alluvial fans derived from growing topographic highs. Foredeep accumulation rates during deposition of the Injana, Mukdadiya and Bai‐Hasan Formations averaged 350, 400 and 600 m/Myr respectively, suggesting accelerated accommodation generation in a rapidly subsiding basin governed by flexural subsidence. Detrital zircon U‐Pb age spectra show that in addition to sources of Mesozoic‐Cenozoic cover strata, the Injana Formation was derived chiefly from Palaeozoic‐Precambrian (including Carboniferous and latest Neoproterozoic) strata in an axial position to the northwest, likely from the Bitlis‐Puturge Massif and broader Eastern Anatolia. In contrast, the Mukdadiya and Bai‐Hasan Formations yield distinctive Palaeogene U‐Pb age peaks, particularly in the southeastern sector of the study region, consistent with transverse delivery from the arc‐related terranes of the Walash and Naopurdan volcano‐sedimentary groups (Gaveh‐Rud domain?) and Urumieh‐Dokhtar magmatic arc to the northeast. These temporal and spatial variations in stratigraphic framework, depositional environments, sediment routing and compositional provenance reveal a major drainage reorganization during Neogene shortening in the Zagros fold‐thrust belt. Whereas axial fluvial systems initially dominated the foreland basin during early orogenesis in the Kurdistan region of Iraq, transverse fluvial systems were subsequently established and delivered major sediment volumes to the foreland as a consequence of the abrupt deformation advance and associated topographic growth in the Zagros.

Detrital zircon (DZ) U‐Pb geochronology is a widely used provenance tool that leverages bedrock age signatures of hinterland source terranes. However, complex sediment recycling of multicycle zircon and hinterland provinces with nondiagnostic U‐Pb ages represent possible pitfalls for provenance reconstructions. Additional biases pertain to source rock zircon fertilities and insensitivity to low‐ and medium‐grade tectonothermal events that do not result in zircon generation. To bridge these inherent biases and gaps in DZ U‐Pb provenance data sets and derive more comprehensive tectonic reconstruction, this study combined U‐Pb and trace element analyses on DZ, apatite, and rutile as well as zircon (U‐Th)/He analyses from the Proto‐Zagros foreland basin in western Iran to shed light on Late Cretaceous tectonic accretion along Arabia. Integrated multimineral, multimethod data sets record formation of a 110‐ to 85‐Ma island arc within Triassic mid‐oceanic ridge crust of the Neotethys ocean, simultaneous obduction starting in Santonian‐Early Campanian times, and inversion of the Arabian rift margin. Overall, integration of these techniques constrains provenance based on multiple independent criteria including crystallization age, cooling history, and petrogenic‐geodynamic environment. This approach not only more completely described provenance signatures but also helped avoid significant pitfalls. For example, while Triassic DZ U‐Pb ages might be mistaken as input from Eurasia, zircon trace element analysis reveals a MORB signature and attributes these DZ to Neotethyan oceanic rather than Eurasian continental origin, having fundamentally different paleogeographic/tectonic implications for the Arabia‐Eurasia collision. Moreover, detrital rutile and apatite resolve and characterize multiple Paleozoic‐Mesozoic thermotectonic events not recorded by DZ U‐Pb due to their largely amagmatic nature.

Zircon (U‐Th)/He (ZHe) and zircon fission track thermochronometric data for 47 samples spanning the areally extensive Willard thrust sheet within the western part of the Sevier fold‐thrust belt record enhanced cooling and exhumation during major thrust slip spanning approximately 125–90 Ma. ZHe and zircon fission track age‐paleodepth patterns along structural transects and age‐distance relations along stratigraphic‐parallel traverses, combined with thermo‐kinematic modeling, constrain the fault slip history, with estimated slip rates of ~1 km/Myr from 125 to 105 Ma, increasing to ~3 km/Myr from 105 to 92 Ma, and then decreasing as major slip was transferred onto eastern thrusts. Exhumation was concentrated during motion up thrust ramps with estimated erosion rates of ~0.1 to 0.3 km/Myr. Local cooling ages of approximately 160–150 Ma may record a period of regional erosion, or alternatively an early phase of limited (<10 km) thrust slip. Propagation of the Sevier wedge front and major thrust slip during the late Early to mid‐Cretaceous were synchronous with increasing subsidence and deposition of thick synorogenic strata in the foreland, crustal thickening in the hinterland, growing igneous activity in the Sierran magmatic arc, and increasing plate convergence rates. Along‐strike, other parts of the Cordilleran retroarc fold‐thrust belt also experienced major shortening during the late Early to mid‐Cretaceous, following a period of earlier Cretaceous quiescence. Late Jurassic shortening was concentrated nearer the arc, and thus mostly in the hinterland at the latitude of northern Utah, related to width and location of the passive‐margin sedimentary wedge relative to the plate margin.