Mid-Cretaceous extensional magmatism in the Alborz Mountains, north Iran; geochemistry and geochronology of Gasht-Masuleh gabbros

In the Gasht-Masuleh area in the Alborz Mountains, gabbroic magma intruded Palaeozoic metasediments and Mesozoic sediments and crystallised as isotropic and cumulate gabbros. LREE enrichment points to relatively low degrees of mantle melting and depletion of Ti, Nb and Ta relative to primitive mantle points to an arc related component in the magma. Clinopyroxene compositions indicate MORB to arc signatures. U–Pb zircon crystallisation ages of 99.5 ± 0.6 Ma and 99.4 ± 0.6 Ma and phlogopite 40Ar/39Ar ages of 97.1 ± 0.4 Ma, 97.5 ± 0.4 Ma, 97.1 ± 0.1 Ma, within 2σ error, indicate that gabbro intrusion occurred in the (Albian-)Cenomanian (mid-Cretaceous). As active subduction did not take place in the Cretaceous in North Iran, the small volume mafic magmatism in the Gasht-Masuleh area must be due to local, extension-related mantle melting. Melting was most likely caused by far field effects triggered by roll-back of the Neo-Tethys subducting slab. As subduction took place at a distance of ~ 400 km (present distance) from the Alborz Mountains, the observed arc geochemical signatures must be inherited from a previous subduction event and concomitant mantle metasomatism, possibly in combination with contamination of the magma by crustal material. Gasht-Masuleh gabbros calc-alkaline, extensional setting, inherited arc signature (Albian-) Cenomanian, 99 Ma zircon ages 97–98 Ma 40Ar/39Ar phlogopite ages Mantle upwelling and formation of gabbros triggered by Neo-Tethys slab roll-back causing extension Gasht-Masuleh gabbros calc-alkaline, extensional setting, inherited arc signature (Albian-) Cenomanian, 99 Ma zircon ages 97–98 Ma 40Ar/39Ar phlogopite ages Mantle upwelling and formation of gabbros triggered by Neo-Tethys slab roll-back causing extension


Introduction
The two most important orogenies in Iran are related to the subduction and closure of the Palaeo-and Neo-Tethys Oceans in the late Triassic and late Cenozoic, respectively (e.g.Stampfli et al., 1991).Palaeo-Tethys closure was associated with accretion of Gondwana-derived Cimmerian terranes or micro-continents, such as the Central Iran micro-continent, to the southern margin of the Eurasian continent resulting in the approximately E-W trending Alborz orogen (Zanchi et al., 2009).Opening of the Neo-Tethys Ocean in the mid-Permian (Stampfli & Kozur, 2006) was followed by Late Triassic to Early Jurassic subduction initiation, the formation of a Jurassicage magmatic arc (Hassanzadeh & Wernicke, 2016 and references therein), and late Oligocene-Miocene collision of Arabia forming the NW-SE trending Zagros orogen to the SW (Agard et al., 2011).
In the Alborz Mountains in the Late Triassic, extension associated with continental alkaline magmatic activity (e.g.Doroozi et al., 2016Doroozi et al., , 2018) ) was followed by deposition of shallow-marine siliciclastic rocks of the Late Triassic-Early Jurassic Shemshak Group.A continuous sequence of Early Jurassic to Cretaceous marine sediments overlies both basement rocks and the Shemshak Group with an unconformity that represents a long-lasting depositional hiatus (Fürsich et al., 2005).
Small volumes of mid-Cretaceous volcanic and plutonic rocks of intermediate to mafic composition and volcano-clastic sediments have been described in the central Alborz Mountains (Agard et al., 2011;Monsef et al., 2022), such as the Gasht-Masuleh gabbros (this study).In late Cretaceous times, thick sequences of limestones and carbonate-rich sediments were deposited in the north-central Alborz.Late Palaeocene-Eocene extension-related volcanism was widespread in the Alborz Orogen (e.g.Karaj Formation, Asiabanha & Foden, 2012) as well as in the adjacent areas of Turkey, Armenia, Georgia and the Talesh Mountains in NW Iran, and was widespread, showing a large basic to acid, and sub-alkaline to alkaline compositional range.
Compared to the widespread and compositionally varied Palaeocene-Eocene volcanism (Verdel et al., 2011), magmatic rocks of Cretaceous age in the Alborz Mountains are rare and relatively small in volume.They are mostly exposed in the central Alborz Mountains, south-east of the Gasht-Masuleh area (Fig. 1).Rezaei et al., (2019Rezaei et al., ( , 2020) ) and Monsef et al. (2022) suggested that the Cretaceous, largely mafic magmas in the Alborz Mountains (and other regions such as the Iran-Armenia-Turkey Plateau; Monsef et al., 2022) were sourced from mantle material that was modified during subduction of the Palaeo-Tethys Ocean, and subsequently underwent fractional crystallisation.Mantle doming and far-field rift magmatism probably resulted from the Neo-Tethys subduction that affected Central Iran as far north as the South Caspian Basin (Monsef et al., 2022).
The occurrence of late Mesozoic mafic mantle-derived magmatic activity remote from active subduction zones and in an environment of quiet carbonate sedimentation raises questions about its age, source, cause and tectonic setting.This paper presents whole rock and mineral compositions, U-Pb zircon and 40 Ar/ 39 Ar phlogopite ages for gabbros in the Gasht-Masuleh area in order to determine the geochemical trends and tectonic setting of their magma.It focuses on the whole-rock composition of the isotropic gabbros (being the closest approximation Fig. 1 Geological map of the Gasht-Masuleh area (simplified after Clark et al., 1975) showing sample locations.GMC = Gasht-Masuleh Metamorphic Complex Mid-Cretaceous extensional magmatism in the Alborz Mountains to magma) and mineral compositions of both isotropic and cumulate samples.In addition, their geochemistry is compared with published data from other Cretaceous mafic magmatic rocks in the Alborz Mountains in order to establish if they were derived from similar mantle sources and resulted from the same tectonic event(s).

Geological background
The Gasht-Masuleh area is located in the western Alborz Mountains to the west of Rasht city (Gilan, also Guilan, Province).Davies et al. (1972) and Clark et al. (1975) mapped the Masuleh area at the 1:100,000 scale and the Bandar-e-Anzali (formerly: Bandar-e-Pahlavi) area at the 1:250,000 scale.Later, the Geological Survey of Iran published the 1:100,000 scale map of the Bandar-e-Anzali area (Nazari et al., 2004).
The Alborz Mountains have a basement of predominantly metasediments of the Gasht Metamorphic Complex and small volumes of late Ediacaran granite (Hassanzadeh et al., 2008).The low-grade part of the Gasht Metamorphic Complex comprises phyllites, greenschist-facies metapelites and metasandstones with subordinate calc-silicates.The higher-grade part consists of amphibolite-facies metapelites (Fig. 1).Protolith sedimentation ages of the metasediments are probably Neoproterozoic (Clark et al., 1975).Amphibolite boulders in river-beds reveal the presence of nearby metabasites which, however, have not been found in exposure due to dense vegetation.The granite from Lahijan east of Rasht has a 551 ± 9 Ma U-Pb zircon age (Guest et al., 2006; hereafter, all age uncertainties are cited at 2-sigma level).In the Talesh Mountains in the West Alborz, lower Palaeozoic volcanic rocks occur interbedded with Silurian limestones and are overlain by Early Devonian limestones (Wendt et al., 2005;Zanchi et al., 2009).
In the study area, the Gasht Metamorphic Complex is overlain by unmetamorphosed Carboniferous black shales and Permian limestones and marls (Clark et al., 1975).The Late Triassic to Middle Jurassic (Bajocian) Shemshak Group overlies these rocks with an angular unconformity and mainly comprises shallow-marine siliciclastic rocks with coal bearing layers, changing to carbonates in Bajocian times (Fürsich et al., 2009).The unconformity is thought to be related to the closure of the Palaeo-Tethys Ocean and terrane accretion (Guest et al., 2006).The younger, mid-to late Jurassic Shal Formation consists of glauconitic and calcareous sandstones that were deposited in a shallow marine environment under oxidizing conditions.Shallow marine carbonates and volcanoclastic rocks of Cretaceous age cover the Jurassic formations (Clark et al., 1975).A granite that intrudes metamorphic basement rocks has an early Jurassic 180.5 ± 1.3 Ma U-Pb zircon crystallisation age (sample AB422B of Chu et al., 2021).Further published U-Pb zircon ages are 99.1 ± 0.7 Ma for a small diorite intrusion in the western Gasht-Masuleh area (Chu et al., 2021), and 97.4 ± 1.8 Ma for the diorite phase and 96.9 ± 2.4 Ma for Nusha granite of a composite intrusion ESE of the Gasht-Masuleh area (Guest et al., 2006).Monsef et al., (2022) reported 99.6 ± 0.5 and 98.5 ± 0.7 Ma U-Pb zircon ages for Javaherdeh-Javaherdasht isotropic gabbros ~ 110 km to the E of Gasht-Masuleh and a 108.1 ± 1.4 Ma zircon U-Pb age for pillow lavas of the Shanderman-Asalem to the NW of the study area.
In the Gasht-Masuleh area, small, ca.1.5 km 2 sized bodies of cumulates and isotropic gabbros are exposed near the villages of Zudel, Gilvanderud and Chapul (Fig. 1).These may the exposed parts of a single, larger gabbro intrusion, but this remains conjectural due to the poor degree of exposure.Their magma intruded the basement metasediments and, locally, Shemshak Group and Shal Formation limestones and marls.

Gasht-Masuleh gabbro intrusions
The gabbros are coarse-to medium-grained clinopyroxene-olivine to plagioclase-rich cumulates and isotropic gabbros.Compared to the isotropic gabbros, the cumulate gabbros volumetrically predominate in the study area.Due to the dense vegetation cover, extensive weathering and thick soil in the area, contacts between gabbros and host rocks are not exposed.Large, loose and locally-derived blocks of the Chapul intrusion in the main valley dissecting the western part of the body are clinopyroxene-and olivine-rich cumulates (Fig. 2A).The Gilvanderud intrusion cuts metasediments of the Gasht Metamorphic Complex and consists of coarse-to medium-grained isotropic gabbros and coarse-grained clinopyroxene-olivine cumulates.The intrusion is crosscut by centimetre to decimetre thick dolerite and diorite dykes (Fig. 2B).The Zudel gabbro intruded the Gasht Metamorphic Complex and sediments of the Shemshak Group and Shal Formation.It consists of coarse-grained isotropic gabbros, and of plagioclase-rich cumulates showing modal layering of plagioclase-rich light-coloured layers and dark-coloured layers composed of pyroxene, phlogopite and opaque minerals (Fig. 2C).The gabbros are crosscut by undated dolerite, diorite dykes and finegrained gabbro dykes of up to 50 cm thick.
Magmatic textures (e.g.layering) and minerals are well preserved and the rocks are not metamorphosed or folded.Brittle deformation and faulting locally resulted in the formation of slickensides on fault planes and tilting of the magmatic layering, the latter feature is especially clear in the Zudel intrusion (Fig. 2D).These features indicate that the intrusions were affected by younger active tectonic events probably during the exhumation of the rocks during the Cenozoic.
The gabbros were considered to be Palaeozoic in age because they intruded the Gasht Metamorphic Complex (Nazari et al., 2004).However, map information indicates that the gabbros intruded Jurassic sediments showing that they must have Jurassic or younger crystallisation ages (Fig. 1).Field relations between the gabbros and the early Jurassic granites are unknown because of lack of exposure.

Petrography
The Gasht-Masuleh gabbros are mainly coarse-grained grey gabbros with isotropic and cumulate textures.The petrography of the cumulate samples is described in Additional file 1; this study focuses on the isotropic gabbros.Samples with ophitic texture are medium-to coarse-grained and contain clinopyroxene up to 2 mm (30-40 vol%), euhedral to subhedral laths of plagioclase up to 1 mm (30-50 vol%), olivine up to 0.5 mm (5-10 vol%), phlogopite (5-10 vol%) up to 1 mm in diameter, and magnetite-ilmenite exsolutions (5-10 vol%) and accessory apatite (up to 2 vol%) (Fig. 3A).Orthopyroxene is rare and occurs as small grains on olivine rims and between olivine and clinopyroxene, and as inclusions in clinopyroxene.The approximate crystallisation order is olivine → olivine ± orthopyroxene → clinopyroxene → plagioclase ± clinopyroxene → opaque mineral(s) → phlogopite (based on petrographic observation listed in Additional file 8).Phlogopite is presents as late an-to subhedral grains, as rims overgrowing opaque minerals, and as anhedral patches locally replacing clinopyroxene crystal margins.Clinopyroxene rims partially reacted to brown and green amphibole; phlogopite and amphibole appear in micro-fractures in clinopyroxene.Olivine is partially altered to serpentine and iddingsite developing along cracks and at the rims.Phlogopite in some samples is altered to secondary unidentified opaque minerals and chlorite along cleavages and at the rims.
Opaque minerals in gabbro samples of especially the Zudel body are inhomogeneous and show extensive evidence of magnetite-ilmenite exsolution (Fig. 3B).The amount of opaque minerals is low, ranging between absent to 5 vol%, in some samples reaching10 vol%.In some samples from the Zudel and Gilvanderud intrusions they are relatively large and range from 0.2 to 1 mm.The largest grains are found in plagioclaserich gabbros, all from the Zudel body.The exsolution textures as observed in back-scatter electron images can be divided into four major textural types following Buddington and Lindsley (1964) and Haggerty (1991) (see the figure caption for the textural types).

Methods and materials
The analytical methods are described in detail in Additional file 1, which also specifies the instrumental parameters, experimental conditions, reference materials, data reduction methods and plateau criteria.18 samples were processed for whole rock major (XRF) and trace element (ICP-OES, ICP-MS) analyses.Wavelength-dispersive electron probe microanalyses (EPMA) were carried out on olivine, ortho-and clinopyroxene, plagioclase and phlogopite in 9 samples.Ce, Pr, Nd, Sr and Sc concentrations of clinopyroxenes, Ba, Rb, F and Cl concentrations in phlogopites, and Rb, Ba, Sr and Eu concentrations of plagioclases in selected samples were determined by field-emission microprobe.Zircons from one isotropic gabbro sample and one gabbroic cumulate sample, and phlogopite from three cumulate samples of isotropic gabbro and gabbroic cumulate were dated by LA-ICP-MS on zircon and 40 Ar/ 39 Ar step-heating of phlogopite.

Whole rock geochemistry
Since the whole rock chemistry of cumulates is controlled by modal proportions of cumulate minerals and not by melt, we only focus on the chemistry of isotropic gabbros, considering these to be the closest representative of magma compositions.These we compared with the chemistry of other mid-Cretaceous basic magmatic rocks in the west Alborz Mountains, especially the basalts from Abrud village in the Gasht-Masuleh area (Rezaei et al., 2020) and the Talesh-Lahijan basaltic lavas (Monsef et al., 2022).
The primitive mantle normalised multi-element patterns (Fig. 5B) are sub-parallel and resemble that of OIB, except that some samples have higher Rb, Ba and K concentrations.Nb, Ta, P and Ti are generally depleted causing negative anomalies to varying degrees.

Mineral chemistry
Selected analyses of olivine, orthopyroxene, clinopyroxene, plagioclase and phlogopite are listed in Additional file 3, Additional file 4, Additional file 5, Additional file 6 and Additional file 7, respectively.

Olivine
Individual crystals lack compositional variations from rim to core and are mainly homogenous.Olivine in Fig. 4 A-I Variation diagrams of MgO vs. SiO 2 , Al 2 O 3 , CaO, FeO t , TiO 2 , Na 2 O, Ni, Sr and Zr for the Gasht-Masuleh gabbros.The chemistry of the gabbros is compared with the basalts from Abrud (Rezaei et al., 2020) and Talesh-Lahijan (Monsef et al., 2022) Mid-Cretaceous extensional magmatism in the Alborz Mountains cumulate samples is relatively magnesium-rich with Mg# between 69.1 and 75.3, compared to olivine in isotropic gabbros (65.7-69.5).NiO contents are very low and range from below detection limit to 0.16 wt% in both cumulates and isotropic gabbros, and do not correlate with forsterite (Fo) content.Cr 2 O 3 abundances are below 0.03 wt% in all olivine-bearing samples.MnO ranges between 0.28 and 0.61 wt% and shows a slight inverse correlation with Fo contents.Olivine in cumulate samples has Fo 68.7-75 , whereas Fo 65-69.1 in isotropic gabbros.Furthermore, Fo component in all samples is lower than those in mantle peridotites (Fo > 83%, Takahashi, 1986).Olivine in Chapul cumulate samples is richer in Fo 71.9-75 compared to those in the Gilvanderud and Zudel cumulates (71.1-73.4% resp.68.7-72.2%).

Clinopyroxene
Mg# of clinopyroxene in cumulate samples varies between 67.4 and 92.3, and in isotropic gabbros between 69.3 and 91.5 (Additional file 5: Table S5).Clinopyroxene in cumulate samples comprises Wo 41.2-51.6En 41.1-51 Fs 0.9-12.7 end member component, whereas clinopyroxene in isotropic gabbros contains Wo 42.5-50.5En 37.5-49.9Fs 2.4-16.5 .Clinopyroxene in plagioclase-bearing cumulate samples has higher amounts of Fs compared to those in olivinebearing cumulates.Clinopyroxene in all three intrusions lacks significant compositional major and trace element zoning.They are diopside to augite according to the classification of Morimoto (1988).

Plagioclase
Plagioclase in most samples is altered to various degrees, and only unaltered parts were analysed.Plagioclase in both cumulate and isotropic samples has similar compositions and is zoned.Anorthite (An) end-member contents range widely from very low 7.1 (probably secondary plagioclase) to 89.2 mol% for cumulate and 11.8 to 92.4 mol% for isotropic samples, with cores having high and rims low An contents.Selected trace element concentrations for plagioclases in three cumulate and two isotropic samples are at or just above detection limits: < 0.02-0.1 wt% Rb 2 O, < 0.04-0.1 wt% BaO, 0.1-0.3wt% SrO and < 0.03-0.3wt% Eu 2 O 3 .Sr and Ba concentrations for low-An rims and high-An cores in both cumulate and isotropic samples overlap and do not correlate with the respective anorthite end-member compositions, but Rb concentrations in rims are at or below the detection limit, and Eu concentrations at the rims are lower on average than in cores.

Phlogopite
Phlogopites generally have high TiO 2 contents of up to 7.6 wt% and high MgO contents of up to 23.9 wt%.CaO  Rudnick and Gao, (2003) contents reach 0.2 wt% (K 2 O contents vary from 6.8 to 9.9 wt%.Mg# for phlogopite in cumulates are 53.0-85.0,and 60.7-83.3 in isotropic gabbros.Many phlogopites are enriched in Cl (0.1-0.4 wt% in cumulates and 0.2-0.5 wt% in isotropic gabbros) and F (up to 0.3 wt% in cumulates and up to 0.2 wt% in isotropic gabbros).Phlogopites in samples 16FMN45A and 16FMN20A that were dated by 40 Ar/ 39 Ar step-heating show no evidence for Ti, Mg, K and Fe compositional zoning.

LA-ICP MS U-Pb zircon dating results
Two samples, a cumulate (15FMN06) and an isotropic gabbro (GH 16-77) were chosen for LA-ICP-MS U-Pb zircon dating.Samples descriptions are given in Additional file 8.The dated zircons are colourless, mainly euhedral long prismatic crystals of up to 100 µm long and ~ 50 µm wide.Some crystals are fragmented due to sample processing.In CL images most crystals show concentric growth zoning typical for magmatic zircons (Fig. 6A, B).U concentrations vary from 174 to 3002 ppm, Th concentrations from 85 and 8764 ppm, and Th/U ratios range from 0.1 to 11.2 ppm.Analysis number (gh1677_13) displays a Th content of 22280 ppm and U content of 3790 ppm, suggesting it is a potential outlier.These concentrations may be analytical artefacts or due to heterogeneity of the zircon grain.However, the Th, U and Pb concentrations of analysis gh1677-13 do not have a significant influence on the calculated age.For the calculation of Concordia ages, only analyses that are 2% or less discordant were used.The analytical data, the ranked and weighted 238 U-206 Pb age average plots, the Concordia diagrams with all analysed points, and the data reporting sheet are given in Additional file 9.

40 Ar/ 39 Ar phlogopite step-heating dating results
Step-heating dating of magmatic phlogopite of three cumulate gabbro samples yielded plateau ages on 2-sigma level for samples 16FMN20A (Zudel) and 16FMN45A (Gilvanderud), but the age spectrum of sample 15FMN05A (Chapul) only allowed calculation of a forced plateau age. 40Ar/ 39 Ar step-heating data including uncertainties are given in Additional file 10.Phlogopite 16FMN20A (Zudel) yielded a 97.1 ± 0.4 Ma plateau age for five gas fractions 8 to 12 comprising 69% of the total 39 Ar released (Fig. 7A), similar to the 97.2 ± 0.4 Ma total gas age.Phlogopite 16FMN45A (Gilvanderud) yielded a 97.5 ± 0.4 Ma plateau age for six gas fractions 4 to 9 comprising 94% of the 39 Ar released (Fig. 7B), which is indistinguishable from the 97.7 ± 0.4 Ma total gas age.In contrast, apparent ages for 15FMN05A (Chapul) did not yield a plateau age, but noting that gas fractions 4 to 10 have apparent ages that are similar (between 94.7 Ma and 98.9 Ma) and almost within 4 Ma of each other, comprising 95% of the 39 Ar released.However, the longest continuous steps which agree within 2-sigma error by their internal errors (excluding common J value errors) were only obtained for steps 5 to 6 (comprising 34% of the 39 Ar released).These steps are here named a forced plateau, which does not satisfy the typical criteria of the plateau that should comprise more than 50% of the 39 Ar released.The weighted mean age calculated from these forced plateau steps is 97.1 ± 0.1 Ma.This forced plateau age is indistinguishable from the 97.2 ± 0.4 Ma total gas age (Fig. 7C) and both ages also overlap the plateau ages of samples 16FMN20A and 16FMN45A within 1-sigma.

Petrography and crystallisation sequence
The crystallisation sequence of the gabbros is olivine, present as isolated individual crystals or as inclusions in clinopyroxene, followed by clinopyroxene, and clinopyroxene and plagioclase.Fe-Ti oxides occur both as inclusions in phlogopite and clinopyroxene, and as individual grains.Phlogopite appears as the latest phase as both subhedral to anhedral crystals, overgrows Fe-Ti oxides and locally replaces clinopyroxene.
The textural varieties of the Fe-Ti oxides (Fig. 3B) indicate a cooling history under oxidizing conditions of an originally homogeneous magmatic Fe-Ti phase, such as ulvöspinel-rich magnetite.Trellis, sandwich and composite textures are super-solvus exsolution textures that formed at temperatures above the magnetite-ulvöspinel solvus, whereas cloth-like to graphic and vermicular textures formed at sub-solvus conditions (Buddington & Lindsley, 1964).Furthermore, crystallisation of an OHbearing magmatic mineral such as phlogopite points to elevated water contents of the magma during the final stages of crystallisation.We suggest that, following crystallisation and (in cumulate samples) gravitational settling of a homogenous Fe-Ti oxide phase, the trellis, sandwich and composite type textures developed upon cooling of the rocks to temperatures of 800-850 °C (e.g.Buddington & Lindsley, 1964) in the presence of small amounts of H 2 O-rich fluids.The fine-grained character of the minerals involved in the cloth-like ilmenite-magnetite exsolutions indicate formation during increased cooling rates at lower temperatures.

Whole rock compositions
The MgO contents of the isotropic gabbros are between 3.9 and 12.6 wt%.This also explains the dominantly linear trends in the variation diagrams that lack inflections (Fig. 4A-I The multi-element patterns for the Gasht-Masuleh isotropic gabbros overlap that of bulk continental crust of Rudnick and Gao (2003) (Fig. 5A).Compared to primitive mantle, the samples are enriched in LILE (Ba, Rb, K) and depleted in HFSE (Nb, Ta, Zr, Ti and P in Fig. 5B) which points to a subduction-related component, the influence of slab derived fluids in their mantle source, and possibly to contamination of magma by continental crust and/ or low degree of mantle melting.Generally, the patterns gently slope from Rb to Yb, and show positive and negative anomalies.Nb and Ta depletions seen in many of the samples point to magma generation in an arc-related setting, or their contamination by continental crust (White & Patchett, 1984), or may be a feature inherited from the mantle source (see Discussion below).Ti depletion may be a primary feature of the mantle source, but can also be due to removal of Fe-Ti oxides during fractionation.
REE and trace elements are considered to describe the character of the of Gasht-Masuleh gabbro magma.In the Nb/Y vs. Ti/Y diagram of Pearce, (1982) the samples trend from the low Ti and low Nb mid-ocean ridge fields to the volcanic arc field together with the basalts from Abrud and Talesh-Lahijan areas, having sub-alkaline Nb/Y ratios < 1 (Fig. 8A).
Nb/La ratios of the isotropic gabbros are lower than those for MORB and OIB (0.2-0.7) suggesting a lithospheric mantle source, and/or contributions by lithospheric and asthenospheric mantle to the Gasht-Masuleh gabbros (Fig. 8B; Smith et al., 1999).
Taking La/Yb as a proxy for the degree of melting, the high La/Yb indicate low degrees of melting.Non-modal melting curves calculated for garnet-bearing and spinel lherzolite mantle sources of depleted and primitive composition are shown on the primitive mantle normalised (La/Yb) N vs. (Dy/Yb) N plot (Fig. 10).This diagram utilises the fractionation of the middle (Dy) from the heavy (Yb) REE during melting in the garnet peridotite stability field, whereas melting in the spinel-facies produces little change in Dy/Yb ratios (Thirlwall et al., 1994).The Gasht-Masuleh gabbros, Abrud and Talesh-Lahijan basalts plot close to the spinel lherzolite melting curves indicating that mid-Cretaceous gabbros and the basalts in the western Alborz Mountains originate from relatively shallow mantle.This is supported by the relatively flat HREE patterns of the Gasht-Masuleh gabbros that do not indicate the presence of garnet in their source ((Dy/Yb) N = 1.4-1.8;Fig. 10).
Comparing published data for basic rocks from the western, central and eastern Alborz (Fig. 9A-E, references for the published data are given in the figure caption), basic rocks in the west and central Alborz and the Gasht-Masuleh gabbros have very similar compositions, whereas those of basic rocks in the east Alborz (Sabzevar area) differ (Additional file 11 shows the location of compared basic rocks).Those in the West Alborz were generated from a mantle source that was more enriched compared to those in the Sabzevar region.Figure 9 shows that basic rocks in the West-Central Alborz were formed in an OIB-like/arc-like tectonic setting, while those in the East Alborz were mostly formed in a MORB setting.As mentioned above, the evolution of the Gasht-Masuleh gabbros magma involved contamination the crustal material.In almost all comparative diagrams, the compositions of the Gasht-Masuleh gabbros and the West-Central Alborz basic rocks overlap with calc-alkaline rocks in the Basin and Range (USA) area (Hawkesworth et al., 1995;Hooper et al., 1995) and in the Eastern Pontides (Turkey) and the Lesser Caucasus (Armenia; Monsef et al., 2022 and references therein) that formed in an extensional tectonic setting.

Mineral compositions
In the Mg# vs. Al 2 O 3 diagram (DeBari & Coleman, 1989) the orthopyroxene does not plot on the Al 2 O 3 enrichment trend but along a low pressure differentiation trend due to their relatively low Al 2 O 3 contents.This points to crystallisation of the magma under low pressure conditions in a shallow magma chamber at upper crustal levels (Fig. 11A).We have found no petrographic evidence, such as plagioclase crystallisation before olivine, which might have supported shallow levels of crystallisation.Orthopyroxene Al 2 O 3 contents range from 0.1 to 1.2 wt.%.The low Al 2 O 3 contents of the orthopyroxenes (< 0.6 apfu) are an indicator for low pressures of crystallization.Using Emslie et al. (1994), we calculate very low pressures (ca.0.23 kbar) that should probably not be taken literally as medium-grained, doleritic textures should be expected, but as an indication for crystallization at shallow crustal levels.Rezaei et al., (2019) point out that the Al 2 O 3 , Na 2 O and Cr 2 O 3 contents of orthopyroxene and clinopyroxene are low, indicating crystallisation at low pressure at upper crustal levels.They obtained temperatures in the range 800-900 °C, using pyroxene thermometry (Köhler The partial melting curves were calculated using the non-modal aggregated fractional melting model C L = (C 0 /F) (1-(1-PF/D 0 ) 1/P ) of Shaw, (1970).The numbers in italics on the model melting curves indicate degree of melting in percent.The modal mineralogies and mineral/melt partition coefficients are from McKenzie andO'Nions (1991a, 1991b) and the melt proportions are from Thirlwall et al., (1994).The primitive and depleted mantle source compositions are from Hofmann (1988).Primitive Mantle normalisation values for La, Dy and Yb are from Sun and McDonough, (1989) & Brey, 1990;Nimis & Taylor, 2000), that are unusually low for mafic melts and do not reflect crystallisation temperatures.
The TiO 2 and Al 2 O 3 contents of the clinopyroxenes plot in the overlapping fields for arc and mid-ocean ridge gabbros (Beccaluva et al., 1989) (Fig. 11D).In the TiO 2 vs. Al 2 O 3 diagram for clinopyroxenes from isotropic and cumulate gabbros (Parlak et al., 2020), the clinopyroxenes mainly cluster in the low pressure field with some analyses plotting in the moderate pressure field and beyond (Fig. 11E).The relatively low TiO 2 and Al 2 O 3 and high SiO 2 of the analysed clinopyroxenes tend to indicate a calc-alkaline nature for the magma, as clinopyroxene in mid-ocean ridge basalts generally are higher, and those in arc-related basalts lower Ti contents (LeBas, 1962).

Zircon and phlogopite ages
The 99.5 ± 0.6 Ma and 99.4 ± 0.6 Ma U-Pb zircon Concordia ages for the Zudel leucogabbro and Gilvanderud isotropic gabbro overlap within 2σ error and are interpreted as crystallisation ages of the magmas that must have intruded in the mid-Cretaceous around the (Albian-)Cenomanian boundary.The 40 Ar/ 39 Ar phlogopite step-heating ages for the three intrusions fall in the narrow range ca.97-98 Ma and are 1 to 2 Ma younger than the U-Pb zircon ages.The younger phlogopite ages agree with the petrographic observations that it is the last phase to crystallise.We consider that the three intrusions crystallised and cooled coevally.(Cartier, 1971;Wensink & Varekamp, 1980).These erupted in a shallow marine basin mainly filled with limestones (Allen et al., 2003) that overlie the Jurassic Shemshak Group with an angular unconformity and are overlain by Aptian-Cenomanian limestones of the Tizkuh Formation (Ansari, 2013;Ansari et al., 2011;Cartier, 1971).Whole rock geochemistry of the olivine basaltic andesites to olivine andesitic basalts of the Kojur and basanites of the Abbasabad volcanic fields indicate that the mantle-derived primary magma were alkaline and originated from a garnet bearing lherzolite source (Ansari, 2013;Ansari et al., 2011).This, and the fact that they erupted onto continental crust points to an intra-continental setting of the Chalus Formation magmatism, e.g. the central Alborz crust (Ansari, 2013;Ansari et al., 2011).Guest et al., (2006 and references therein), reported a 97.4 ± 1.8 Ma ionprobe 206 Pb/ 238 U zircon age for the diorite phase and a 96.9 ± 2.4 Ma age for the (rapakivi) granite phase of the composite Nusha Pluton 150 km to the ESE of the Gasht-Masuleh area.These zircon ages overlap with ca.97-100 Ma 40 Ar/ 39 Ar phlogopite and U-Pb zircon ages for the Gasht-Masuleh gabbros (this study).This indicates that (Albian-)Cenomanian plutonic activity was not local but occurred over large areas, possibly constituting a long-lasting and widespread mafic volcanoplutonic stage in the Alborz Mountains.120 km to the east of the Gasht-Masuleh area on the Caspian Sea coast, basic volcanic rocks, small gabbro intrusions and ultramafic cumulates of suspected Cretaceous age occur south of Amlash (Annells, 1975).The gabbros and cumulates have a mantle source that was modified by subduction (Salavati et al., 2013).In contrast, the basalts are alkaline and have an OIB-like mantle source (Salavati, 2008).In the Javaherdasht area 20 km SE of Amlash, isotropic and layered gabbros intrude Jurassic clastic sediments and limestone of the Javaherdeh Formation.The isotropic gabbros have high Na 2 O, K 2 O, TiO 2 and FeO t contents compared to the layered gabbros and were probably derived from MORB-like spinel-bearing mantle source(s) and show trace element evidence for contamination with crustal material (Haghnazar et al., 2011).The layered gabbros have high Mg# and high Ni, Cr and Co abundances, both types may have formed in a continental rift setting (Haghnazar et al., 2011).Basalts and gabbros in the Javaherdasht, the Gasht-Masuleh and the Talesh Mountains areas formed in an extensional setting but show an arc geochemical signature that is the result of contamination by crustal components (Haghnazar et al., 2011;Monsef et al., 2022;Rezaei et al., 2019Rezaei et al., , 2020)).Cretaceous magmatic rocks in the Deylaman area (Additional file 11) comprise pillow and sheet lavas, trachytes, alkali rhyolites and olivine gabbro and monzodiorite intrusions with sub-alkaline to alkaline signatures (Akmali et al., 2019).Mid-Cretaceous extensional magmatism in the Alborz Mountains

Comparison with other
Comparison of the gabbros in the Gasht-Masuleh area with published data of other Cretaceous mafic magmatic rocks in the Alborz Mountains shows that more than half have elevated Ba/La (i.e. higher than MORB) at relatively low Th/Yb ratios (Fig. 12; Woodhead et al., 2001).This may be attributed to the effects of fluids derived from a previously subducted oceanic slab (Arculus & Powell, 1986;Woodhead et al., 2001), possibly of the Palaeo-Tethys Ocean because most magmatic activities related to Neo-Tethys subduction are concentrated in the Sanandaj-Sirjan Zone far to the south of the Alborz Mountains (Hassanzadeh & Wernicke, 2016).The volcanic rocks in the Chalus area form an exception, trending to high Th/Yb at low Ba/La ratios reflecting mixing of mantle source material with magma derived from subducted sediments (Ansari, 2013;Ansari et al., 2011).In the Alborz Mountains, tectonic activity was limited in the Cretaceous.Sedimentation was of a quiet nature with low sedimentation rates of predominant shallow marine reef carbonates and marls (Clark et al., 1975;Saidi et al., 1997) and subduction was cetrainly not taking place in this area.Guest et al. (2006) tentatively suggested that the Cretaceous magmatic activity might be related to subduction of the Neo-Tethys Ocean and arc formation in the Zagros Belt at least 380 km to the SW, implying that a very wide area was affected by this subduction system.Thus, the mantle derived magmatism in the Alborz Mountains may be due to extension causing mantle melting, both triggered by far-field effects of the Neo-Tethys subduction zone.A similar setting was proposed for the Permo-Carboniferous mafic magmatism in the foreland of the Variscan Orogen in Europe (e.g.Timmerman, 2004).
Geochemical data for the Gasht-Masuleh gabbros such as Ti, Nb and Ta negative anomalies, low Nb/Th and elevated and variable Ba/La ratios, indicate that their magma was derived from a mantle source marked an arc signature.In general, mantle metasomatism appears to be mainly the result of slab-derived fluids, except for the alkaline basalts of the Chalus Formation that show evidence for contributions from subducted sediments (Fig. 12).The Cenomanian plate tectonic reconstruction of Barrier et al., (2018) (Fig. 13) illustrates that north-dipping Neo-Tethys subduction occurred in the future Zagros Belt and that in the northeastern, overriding plate, extension-related magmatism took place in the Alborz Mountains.Here, magmatism occurred in a shallow marine environment and was associated with faulting that probably facilitated the magma ascent such as that of the Gasht-Masuleh gabbros.As there is no evidence for Jurassic-Cretaceous subduction in the Alborz  Woodhead et al., (2001).Data sources: Abrud basalt (Rezaei et al., 2020), Talash-Lahijan basalt (Monsef et al., 2022), Masuleh gabbro (Khalatbari Jafari et al., 2016), Javaherdasht gabbro (Haghnazar et al., 2011;Monsef et al., 2022), South Amlash basic rocks (Salavati, 2008;Salavati et al., 2013), Marzanabad basalt (Doroozi et al., 2016), Chalus basalt (Ansari, 2013) and Deylaman igneous complex (Akmali et al., 2019) Fig. 13 Palaeotectonic reconstruction of Iran (modified from Barrier et al., 2018) showing the approximate position of the study area (star symbol) in Cenomanian in relation to the zone of active subduction of the Neo-Tethys slab beneath the Central Iranian Platform, causing basin formation and extension-related magmatism in the Alborz Mountains Mid-Cretaceous extensional magmatism in the Alborz Mountains Mountains, these signatures must have been inherited from an older event, probably subduction of the Palaeo-Tethys Ocean in the Carboniferous and its closure in the Triassic (Stampfli & Borel, 2002).The Palaeo-Tethys suture is assumed to lie to the north of the Alborz Mountains (e.g.Omrani et al., 2013;Stampfli & Kozur, 2006).Local normal faulting and the generation of small volumes of magmatic rocks in the Alborz Mountains in the Cretaceous were probably the response to extension and the far field tectonic effects of slab roll-back and back-arc opening due to the north dipping Neo-Tethys subduction zone in central Iran (Agard et al., 2011).The short-lived, local, relatively small volume and often mafic mid-Cretaceous magmatism in the Gasht-Masuleh area may be linked to such short-lived extension stages that caused passive mantle upwelling and low degree mantle decompression melting (Fig. 14).Relatively rapid ascent of the mantle-derived magma probably took place along deepseated fault zones.These may have been newly formed normal faults or reactivated old faults (including thrusts and strike-slip faults) inherited from closure of the Palaeo-Tethys Ocean in the Late Triassic.

Conclusions
Dating of magmatic zircon and phlogopite shows that in the Gasht-Masuleh area, intrusion and crystallisation of mantle derived mafic magma occurred at ca. 100 Ma during the mid-Cretaceous times.The magma shows normal gabbroic fractionation (Ol → Ol + Cpx + magnetite-ilmenite ± Opx → Pl + magnetite-ilmenite).The relatively low Mg# of the olivines indicates that the magma was fractionated before intrusion emplacement.Crystallisation occurred at relatively high crustal levels indicated by the low Al contents of clino-and orthopyroxenes, defining a low pressure differentiation trend.Clinopyroxene compositions point to a low Ti, calc-alkaline magma, which is supported by the transitional to calc-alkaline whole rock major and trace element compositions.
Gasht-Masuleh gabbros are slightly enriched in LREE resulting in sloping REE patterns that generally lack Eu anomalies.Although there is no evidence for active subduction and arc magmatism in Alborz in the Cretaceous, Nb and Ta negative anomalies in multi-element diagrams point to the presence of an arc-related component, as do the variable Ti abundances which, however, may also be due to the fractionation of magnetite-ilmenite exsolution.Overall, the trace element signature points to the presence of an arc-related component inherited from the mantle source, with possible contamination by crustal components.
In summary, trace element compositions indicate that the gabbros in the Gasht-Masuleh area (and many other Cretaceous mafic magmatic rocks in the Alborz Mountains) were generated by partial melting of a metasomatized mantle source.The most likely explanation is that mantle metasomatism occurred during subduction of the Palaeo-Tethys Ocean in north Iran in Palaeozoic-Triassic times, and that the arc/subduction signature of the gabbros is an inherited feature (e.g.Monsef et al., 2022;Woodhead et al., 2001).
U-Pb zircon and 40 Ar/ 30 Ar phlogopite ages show that the Gasht-Masuleh gabbros intruded in (Albian-)Cenomanian times and cooled rapidly, indicating magma generation, intrusion and cooling at high crustal levels all taking place during a short-lived extensional event.This event was not restricted to the Gasht-Masuleh area but also occurred further to the SSE in the area NW of Tehran (ca.97 Ma Nusha diorite-granite composite pluton, Guest et al., 2006).In general, the mid-Cretaceous magmatism may reflect extension of the crust caused by slab roll-back initiation of the northward subducting Neo-Tethys Ocean.

Fig. 2
Fig. 2 Field photographs of the Gasht-Masuleh gabbros, A Cumulate gabbro from the Chapul intrusion.B Typical exposure of gabbros of the Gilvanderud intrusion, the gabbros are jointed and weathered on the surface.C Hand specimen of plagioclase-rich gabbro from the Zudel intrusion showing modal magmatic layering, with accumulation of plagioclase in the light layers and of mainly clinopyroxene, phlogopite and Fe-Ti oxides in dark layers.D Gabbro of the Zudel intrusion with tilted modal magmatic layering (highlighted by white lines) Figure 4A-I shows the variation diagrams for whole rock SiO 2 , Al 2 O 3 , CaO, FeO t , TiO 2 , Na 2 O, Sr, Ni and Zr concentrations against MgO as fractionation index.They show that SiO 2 concentrations stay within the narrow range of 41-48 wt% whereas CaO increases and Al 2 O 3 , TiO 2 and Na 2 O decrease with increasing MgO.FeO t shows no recognizable correlation and, except two samples (16FMN48B and 16FMN70A).Sr, Zr contents decrease and Ni increases with increasing MgO.The Gasht-Masuleh gabbros do not show distinct compositional differences from the basalts from Abrud and Talesh-Lahijan areas.

Fig. 3 A
Fig. 3 A Photomicrographs (crossed nicols) of samples of Gasht-Masuleh isotropic gabbro sample 16FMN22D with ophitic texture from Zudel intrusion.B Back-scattered electron image of Fe-Ti oxides in sample 16FMN22C of the Zudel intrusion showing textures of trellis type: thin ilmenite spindles in magnetite displaying sharp contacts with the magnetite host, sandwich type: thick, lath-shaped ilmenite lamellae in magnetite with sharp mutual boundaries, composite type: aggregates in which ilmenite surrounds or is partly included in magnetite with sharp contact, and cloth-like to graphic and vermicular: fine-grained exsolutions of anhedral ilmenite and magnetite.Mineral name abbreviations in the text, tables and figures are afterWhitney and Evans, (2010)

Fig. 5 A
Fig. 5 A Chondrite normalised rare earth element patterns of the Gasht-Masuleh isotropic gabbros showing moderate LREE enrichment compared to the HREE.B Primitive Mantle normalised multi-element pattern of the gabbros overlapping the bulk continental crust composition.REE and multi-element patterns from the Abrud (Rezaei et al., 2020) and Talesh-Lahijan (Monsef et al., 2022) areas are plotted for chemical comparison.Chondrite and Primitive Mantle normalisation values, and OIB and N-MORB are from Sun and McDonough, (1989); bulk continental crust values are fromRudnick and Gao, (2003)

Fig. 6
Fig. 6 Cathodoluminescence images and Concordia diagrams of magmatic zircons from the Gasht-Masuleh gabbros A Cumulate gabbro 15FMN06 of the Zudel intrusion.B isotropic gabbro GH16-77 of the Gilvanderud intrusion.The spots where laser-ablation analyses were performed are labelled with spot number and 206 Pb/ 238 U ages.Data-point error ellipses are 2σ.(D) denotes the analyses that are more than 2 percent discordant

Fig. 7 .
Fig.7.40 Ar/ 39 Ar incremental step heating age spectra, Ca/K ratios and fractions (%) of the radiogenic.40Ar of magmatic phlogopite of cumulate samples A 16FMN20A (Zudel intrusion) and B 16FMN45A (Gilvanderud intrusion) yielding plateau ages.C A forced plateau age (*) for 2 gas fractions of step 5-6 was calculated for phlogopite from 15FMN05A (Chapul intrusion) ).The variation diagrams show that Al 2 O 3 and Na 2 O decrease and CaO increases with increasing MgO concentrations, overall indicating fractional crystallisation of clinopyroxene before plagioclase.Samples with low MgO have an elevated TiO 2 concentrations caused by the presence of large phlogopite crystals and (especially for the Zudel intrusion) magnetite-ilmenite exsolutions.Ni concentrations increase with increasing MgO due to the presence olivine.Like Al 2 O 3 , Sr concentrations decrease with increasing MgO in gabbros reflecting fractionation of plagioclase.The major oxides and trace element compositions of the isotropic gabbros overlap with the compositions of basalts from Abrud and Talesh-Lahijan areas.