The Neotethyan Sanandaj-Sirjan zone of Iran as an archetype for passive margin-arc transitions

The Sanandaj-Sirjan zone of Iran is a northwest trending orogenic belt immediately north of the Zagros suture, which represents the former position of the Neotethys Ocean. The zone contains the most extensive, best preserved record of key events in the formation and evolution of the Neotethys, from its birth in Late Paleozoic time through its demise during the mid-Tertiary collision of Arabia with Eurasia. The record includes rifting of continental fragments off of the northern margin of Gondwanaland, formation of facing passive continental margins, initiation of subduction along the northern margin, and progressive development of a continental magmatic arc. The latter two of these events are critical phases of the Wilson Cycle that, elsewhere in the world, are poorly preserved in the geologic record because of superimposed events. Our new synthesis reaffirms the similarity between this zone and various terranes to the north in Central Iran. Late Paleozoic rifting, preserved as A-type granites and accelerated subsidence, was followed by a phase of pronounced subsidence and shallow marine sedimentation in Permian through Triassic time, marking the formation and evolution of passive margins on both sides of the suture. Subduction and arc magmatism began in latest Triassic/Early Jurassic time, culminating at ~170 Ma. The extinction of arc magmatism in this zone, and its shift northeastward to form the subparallel Urumieh-Dokhtar arc, occurred diachronously along strike, in Late Cretaceous or Paleogene time. Post-Cretaceous uplift transformed the zone from a primarily marine borderland into a marine archipelago that persisted until mid-Tertiary time

Comment on "Neotethyan subduction ignited the Iran arc and back-arc differently" by Shafaii Moghadam et al. (2020)

Shafaii Moghadam et al. (2020) contribute important new data on Late Cretaceous-Tertiary subduction- related magmatism in Iran, but their plate convergence model, wherein Neotethyan subduction begins in mid-Cretaceous time (c. 100 Ma), overlooks well established facts relating to the tectonic history of Neotethys, in regard to global plate reconstructions, paleolatitude data, the regional stratigraphy, geochronology and geochemistry, and metamorphic history. Based on their model, Neotethys subduction beneath Eurasia began at ~100 Ma, meaning that the Neotethys was spreading and bounded by opposing passive margins during Jurassic and Early Cretaceous time, for ~100 Ma prior to their proposed onset of Neotethyan convergence. Consequently, their subduction model contradicts (1) the Indian Ocean spreading history derived from magnetic anomalies; (2) continental paleolatitude data from paleomagnetism; (3) sedimentary and igneous evolution of the Mesozoic continental margins in Arabia and southern Asia, (4) the age and geochemistry of Jurassic igneous rocks in southernmost Eurasia; and (5) the preservation of Early to Middle Jurassic eclogite metamorphism and exhumation on the northern side of the Arabia-Eurasia suture. Reconciliation of each of these omissions and contradictions of their model would be welcome, and perhaps an advisory that readers may wish to evaluate their concept of Cretaceous subduction initiation with due circumspection

Comment on "Neotethyan subduction ignited the Iran arc and back-arc differently" by Shafaii Moghadam et al. (2020)

Shafaii Moghadam et al. (2020) contribute important new data on Late Cretaceous-Tertiary subduction- related magmatism in Iran, but their plate convergence model, wherein Neotethyan subduction begins in mid-Cretaceous time (c. 100 Ma), overlooks well established facts relating to the tectonic history of Neotethys, in regard to global plate reconstructions, paleolatitude data, the regional stratigraphy, geochronology and geochemistry, and metamorphic history. Based on their model, Neotethys subduction beneath Eurasia began at ~100 Ma, meaning that the Neotethys was spreading and bounded by opposing passive margins during Jurassic and Early Cretaceous time, for ~100 Ma prior to their proposed onset of Neotethyan convergence. Consequently, their subduction model contradicts (1) the Indian Ocean spreading history derived from magnetic anomalies; (2) continental paleolatitude data from paleomagnetism; (3) sedimentary and igneous evolution of the Mesozoic continental margins in Arabia and southern Asia, (4) the age and geochemistry of Jurassic igneous rocks in southernmost Eurasia; and (5) the preservation of Early to Middle Jurassic eclogite metamorphism and exhumation on the northern side of the Arabia-Eurasia suture. Reconciliation of each of these omissions and contradictions of their model would be welcome, and perhaps an advisory that readers may wish to evaluate their concept of Cretaceous subduction initiation with due circumspection

Crust-mantle interaction inferred from the petrology and Sr-Nd-Pb isotope geochemistry of Eocene arc lavas from the Kahrizak Mountains, north-Central Iran

The Eocene volcanic rocks from the Kahrizak Mountains in north-central Iran are part of the Urumieh-Dokhtar magmatic arc, which runs parallel to the Main Zagros Thrust segment of the Neo-Tethys suture. These volcanic rocks, similar to those from eastern Pontides and northern Anatolia, Turkey, were mainly produced during the Eocene magmatic flare-up associated with the Arabia-Eurasia convergence. The rock suite includes basalt, trachyandesite/andesite and trachydacite/rhyolite lavas and pyroclastic deposits that evolved compositionally from calc-alkalic to shoshonitic. Their normalized trace element concentration patterns are moderately enriched in light rare earth element and depleted in high field-strength elements (HFSE; e.g., Nb, Ta, Ti). They have narrow ranges of initial Pb isotopic ratios and ^(143)Nd/^(144)Nd_i, but highly variable ^(87)Sr/^(86)Sr_i. The new analyses indicate that the parental magmas of the volcanic rocks were derived from a mantle source that had been enriched by fluids released from a subducted oceanic slab. The fluids introduced significant amounts of large ion lithophile elements, but negligible HFSE to the source. The parental magmas underwent fractional crystallization and assimilation of upper crustal materials to produce the range of volcanic rocks. Integration of new analyses with regional data suggests that the Eocene volcanic rocks from north-central Iran, together with ~coeval volcanic rocks in eastern Pontides and northern Anatolia, were most probably derived from a lithospheric mantle source that had been previously metasomatized by fluids derived from a subducted slab before and during the Arabia-Eurasia collision

Crystal size and shape distribution systematics of plagioclase and the determination of crystal residence times in the micromonzogabbros of Qisir Dagh, SE of Sabalan volcano (NW Iran)

The Qisir Dagh igneous complex occurs as a combination of volcanic and intrusive rocks to the south-east of the Sabalan volcano, north-western Iran. Micromonzogabbroic rocks in the region consist of plagioclase, alkaline feldspar and clinopyroxene as the major mineral phases and orthopyroxene, olivine, apatite and opaque minerals as the accessory minerals. Microgranular and microporphyritic textures are well developed in these rocks. Considering the importance of plagioclase in reconstructing magma cooling processes, the size and shape distribution and chemical composition of this mineral were investigated. Based on microscopic studies, it is shown that the 2-dimensional size average of plagioclase in the micromonzogabbros is 538 micrometers and its 3-dimensional shape varies between tabular to prolate. Crystal size distribution diagrams point to the presence of at least two populations of plagioclase, indicating the occurrence of magma mixing and/or fractional crystallization during magma cooling. The chemical composition of plagioclase shows a wide variation in abundances of Anorthite-Albite-Orthoclase (An = 0.31—64.58, Ab = 29.26—72.13, Or = 0.9—66.97), suggesting a complex process during the crystal growth. This is also supported by the formation of antiperthite lamellae, which formed as the result of alkali feldspar exsolution in plagioclase. The calculated residence time of magma in Qisir Dagh, based on 3D crystal size distribution data, and using growth rate G = 10—10 mm/s, varies between 457 and 685 years, which indicates a shallow depth (near surface) magma crystallization and subvolcanic nature of the studied samples

Lithospheric structure beneath the Zagros collision zone resolved by non-linear teleseismic tomography

The upper-mantle structure across the Zagros collision zone, in southwest Iran, is investigated using a non-linear weighted damped least-squares teleseismic tomography approach. The resolution of the structures/transitions in the upper mantle is enhanced significantly by correcting the teleseismic relative arrival time residuals for an a priori crustal velocity model and then performing the inversion with fixed crustal blocks. To investigate whether or not the lithospheric blocks and major transitions in the resulting model are required by the data or are artefacts of the inversion, the data were inverted using two different inverse methods (singular value decomposition and a quadratic programming method). New high-quality seismic velocity models show apparent correlation between surface geological features and seismic velocity structures at lithospheric depth across the Zagros collision zone. The image shows a sharp lithospheric boundary at the Main Zagros Thrust between 100 km and 250 km depth with P-wave velocity about 3 per cent faster within the Arabian Shield to the south. A step-like increase in lithospheric thickness across the Zagros collision zone is assumed to separate two different mantle structures namely the Arabian (to the south) and the Eurasian (to the north) domains. The most striking feature resolved is a north-dipping slab-like positive velocity anomaly

The Neotethyan Sanandaj-Sirjan zone of Iran as an archetype for passive margin-arc transitions

The Sanandaj-Sirjan zone of Iran is a northwest trending orogenic belt immediately north of the Zagros suture, which represents the former position of the Neotethys Ocean. The zone contains the most extensive, best preserved record of key events in the formation and evolution of the Neotethys, from its birth in Late Paleozoic time through its demise during the mid-Tertiary collision of Arabia with Eurasia. The record includes rifting of continental fragments off of the northern margin of Gondwanaland, formation of facing passive continental margins, initiation of subduction along the northern margin, and progressive development of a continental magmatic arc. The latter two of these events are critical phases of the Wilson Cycle that, elsewhere in the world, are poorly preserved in the geologic record because of superimposed events. Our new synthesis reaffirms the similarity between this zone and various terranes to the north in Central Iran. Late Paleozoic rifting, preserved as A-type granites and accelerated subsidence, was followed by a phase of pronounced subsidence and shallow marine sedimentation in Permian through Triassic time, marking the formation and evolution of passive margins on both sides of the suture. Subduction and arc magmatism began in latest Triassic/Early Jurassic time, culminating at ~170 Ma. The extinction of arc magmatism in this zone, and its shift northeastward to form the subparallel Urumieh-Dokhtar arc, occurred diachronously along strike, in Late Cretaceous or Paleogene time. Post-Cretaceous uplift transformed the zone from a primarily marine borderland into a marine archipelago that persisted until mid-Tertiary time

Flexural bending of the Zagros foreland basin

We constrain and model the geometry of the Zagros foreland to assess the equivalent elastic thickness of the northern edge of the Arabian plate and the loads that have originated due to the Arabia–Eurasia collision. The Oligo-Miocene Asmari formation, and its equivalents in Iraq and Syria, is used to estimate the post-collisional subsidence as they separate passive margin sediments from the younger foreland deposits. The depth to these formations is obtained by synthesizing a large database of well logs, seismic profiles and structural sections from the Mesopotamian basin and the Persian Gulf. The foreland depth varies along strike of the Zagros wedge between 1 and 6 km. The foreland is deepest beneath the Dezful embayment, in southwest Iran, and becomes shallower towards both ends. We investigate how the geometry of the foreland relates to the range topography loading based on simple flexural models. Deflection of the Arabian plate is modelled using point load distribution and convolution technique. The results show that the foreland depth is well predicted with a flexural model which assumes loading by the basin sedimentary fill, and thickened crust of the Zagros. The model also predicts a Moho depth consistent with Free-Air anomalies over the foreland and Zagros wedge. The equivalent elastic thickness of the flexed Arabian lithosphere is estimated to be ca. 50 km. We conclude that other sources of loading of the lithosphere, either related to the density variations (e.g. due to a possible lithospheric root) or dynamic origin (e.g. due to sublithospheric mantle flow or lithospheric buckling) have a negligible influence on the foreland geometry, Moho depth and topography of the Zagros. We calculate the shortening across the Zagros assuming conservation of crustal mass during deformation, trapping of all the sediments eroded from the range in the foreland, and an initial crustal thickness of 38 km. This calculation implies a minimum of 126 ± 18 km of crustal shortening due to ophiolite obduction and post-collisional shortening

Crust-mantle interaction inferred from the petrology and Sr-Nd-Pb isotope geochemistry of Eocene arc lavas from the Kahrizak Mountains, north-Central Iran

The Eocene volcanic rocks from the Kahrizak Mountains in north-central Iran are part of the Urumieh-Dokhtar magmatic arc, which runs parallel to the Main Zagros Thrust segment of the Neo-Tethys suture. These volcanic rocks, similar to those from eastern Pontides and northern Anatolia, Turkey, were mainly produced during the Eocene magmatic flare-up associated with the Arabia-Eurasia convergence. The rock suite includes basalt, trachyandesite/andesite and trachydacite/rhyolite lavas and pyroclastic deposits that evolved compositionally from calc-alkalic to shoshonitic. Their normalized trace element concentration patterns are moderately enriched in light rare earth element and depleted in high field-strength elements (HFSE; e.g., Nb, Ta, Ti). They have narrow ranges of initial Pb isotopic ratios and ^(143)Nd/^(144)Nd_i, but highly variable ^(87)Sr/^(86)Sr_i. The new analyses indicate that the parental magmas of the volcanic rocks were derived from a mantle source that had been enriched by fluids released from a subducted oceanic slab. The fluids introduced significant amounts of large ion lithophile elements, but negligible HFSE to the source. The parental magmas underwent fractional crystallization and assimilation of upper crustal materials to produce the range of volcanic rocks. Integration of new analyses with regional data suggests that the Eocene volcanic rocks from north-central Iran, together with ~coeval volcanic rocks in eastern Pontides and northern Anatolia, were most probably derived from a lithospheric mantle source that had been previously metasomatized by fluids derived from a subducted slab before and during the Arabia-Eurasia collision

Middle to late Cenozoic basin evolution in the western Alborz Mountains: Implications for the onset of collisional deformation in northern Iran

Oligocene-Miocene strata preserved in synclinal outcrop belts of the western Alborz Mountains record the onset of Arabia-Eurasia collision-related deformation in northern Iran. Two stratigraphic intervals, informally named the Gand Ab and Narijan units, represent a former basin system that existed in the Alborz. The Gand Ab unit is composed of marine lagoonal mudstones, fluvial and alluvial-fan clastic rocks, fossiliferous Rupelian to Burdigalian marine carbonates, and basalt flows yielding ^(40)Ar/^(39)Ar ages of 32.7 ± 0.3 and 32.9 ± 0.2 Ma. The Gand Ab unit is correlated with the Oligocene–lower Miocene Qom Formation of central Iran and is considered a product of thermal subsidence following Eocene extension. The Narijan unit unconformably overlies the Gand Ab unit and is composed of fluvial-lacustrine and alluvial fan sediments exhibiting contractional growth strata. We correlate the Narijan unit with the middle to upper Miocene Upper Red Formation of central Iran on the basis of lithofacies similarities, stratigraphic position, and an 8.74 ± 0.15 Ma microdiorite dike (^(40)Ar/^(39)Ar) that intruded the basal strata. Deformation timing is constrained by crosscutting relationships and independent thermochronological data. The Parachan thrust system along the eastern edge of the ancestral Taleghan-Alamut basin is cut by dikes dated at 8.74 ± 0.15 Ma to 6.68 ± 0.07 Ma (^(40)Ar/^(39)Ar). Subhorizontal gravels that unconformably overlie tightly folded and faulted Narijan strata are capped by 2.86 ± 0.83 Ma (^(40)Ar/^(39)Ar) andesitic lava flows. These relationships suggest that Alborz deformation had migrated southward into the Taleghan-Alamut basin by late Miocene time and shifted to its present location along the active range front by late Pliocene time. Data presented here demonstrate that shortening in the western Alborz Mountains had started by late middle Miocene time. This estimate is consistent with recent thermochronological results that place the onset of rapid exhumation in the western Alborz at ∼12 Ma. Moreover, nearly synchronous Miocene contraction in the Alborz, Zagros Mountains, Turkish-Iranian plateau, and Anatolia suggests that the Arabia-Eurasia collision affected a large region simultaneously, without a systematic outward progression of mountain building away from the collision zone