The petrological and geochronological study of the Cenozoic Shaivar Dagh composite intrusion in the Alborz Mountain belt (NW Iran) reveals important clues to decipher complex relations between magmatic and tectonic processes in the central sectors of the Tethyan (Alpine–Himalayan) orogenic belt. This pluton is formed by intrusion at different times of two main magmatic cycles. The older (Cycle 1) is formed by calc-alkaline silicic rocks, which range in composition from diorites to granodiorites and biotite granites, with abundant mafic microgranular enclaves. The younger cycle (Cycle 2) is formed by K-rich monzodiorite and monzonite of marked shoshonitic affinity. The latter form the larger volumes of the exposed plutonic rocks in the studied complex. Zircon geochronology (laser ablation ICP-MS analyses) gives a concordia age of 30.8 ± 2.1 Ma for the calc-alkaline rocks (Cycle 1) and a range from 23.3 ± 0.5 to 25.1 ± 0.9 Ma for the shoshonitic association (Cycle 2). Major and trace element relations strongly support distinct origins for each magmatic cycle. Rocks of Cycle 1 have all the characteristic features of active continental margins. Shoshonitic rocks (Cycle 2) define two continuous fractionation trends: one departing from a K-rich basaltic composition and the other from an intermediate, K-rich composition. A metasomatized-mantle origin for the two shoshonitic series of Cycle 2 is proposed on the basis of comparisons with experimental data. The origin of the calc-alkaline series is more controversial but it can be attributed to processes in the suprasubduction mantle wedge related to the incorporation of subducted mélanges in the form of silicic cold plumes. A time sequence can be established for the processes responsible of the generation of the two magmatic cycles: first a calc-alkaline cycle typical of active continental margins, and second a K-rich cycle formed by monzonites and monzodiorites. This sequence precludes the younger potassic magmas as precursors of the older calc-alkaline series. By contrast, the older calc-alkaline magmas may represent the metasomatic agents that modified the mantle wedge during the last stages of subduction and cooked a fertile mantle region for late potassic magmatism after continental collision.

Quantified balanced and restored crustal cross-sections across the NW Zagros Mountains are presented in this work integrating geological and geophysical local and global datasets. The balanced crustal cross-section reproduces the surficial folding and thrusting of the thick cover succession, including the near top of the Sarvak Formation (~90 Ma) that forms the top of the restored crustal cross-section. The base of the Arabian crust in the balanced cross-section is constrained by recently published seismic receiver function results showing a deepening of the Moho from 42 ± 2 km in the undeformed foreland basin to 56 ± 2 km beneath the High Zagros. The internal parts of the deformed crustal cross-section are constrained by new seismic tomographic sections imaging a ~50° NE-dipping sharp contact between the Arabian and Iranian crusts. These surfaces bound an area of 10800 km2 that should be kept constant during the Zagros orogeny. The Arabian crustal cross-section is restored using six different tectonosedimentary domains according to their sedimentary facies and palaeobathymetries, and assuming Airy isostasy and area conservation. While the two southwestern domains were directly determined from well-constrained surface data, the reconstruction of the distal domains to the NE was made using the recent margin model of Wrobel-Daveau et al. (2010) and fitting the total area calculated in the balanced cross-section. The Arabian continental–oceanic boundary, at the time corresponding to the near top of the Sarvak Formation, is located 169 km to the NE of the trace of the Main Recent Fault. Shortening is estimated at ~180 km for the cover rocks and ~149 km for the Arabian basement, including all compressional events from Late Cretaceous to Recent time, with an average shortening rate of ~2 mm yr−1 for the last 90 Ma.