Petrology , Geochemistry and Tectonomagmatic Setting of Neshveh Intrusion ( NW Saveh )

Neshveh intrusion is located in the NW of Saveh City and is a part of Orumieh-Dokhtar magmatic arc. This intrusion consists of quartz monzodiorite, granodiorite and granite that have intruded into the Eocene volcano-sedimentary rocks. This intrusion is high-K calk-alkaline and metaluminous and is classified as I-type granitoids. Field investigations along with petrographic and geochemical studies indicate that all phases of Neshveh intrusion are derived from a common magma source as a result of mineral differentiation. Different phases of this intrusion have low Mg#, Ni, Cr, Co and V which are indicative for higher evolution of magma during the magma ascent and before complete crystallization. All phases of Neshveh granitoid are characterized by LREE-rich patterns with high LREE/HREE ratio and negative Eu anomalies. Similarity of the mentioned patterns suggests a comagmatic source for these rocks and demonstrates the role of magmatic differentiation in their evolution. There are negative anomalies in the Nb and Ti along with positive anomalies of Rb, Ba, K and Pb on the spider diagrams. These anomalies are indicative for a subduction setting for magma source of these rocks. Geochemical studies indicate that the Neshveh granitoid is formed in a volcanic arc and active continental margin. In this base, it is assumed that this intrusion is formed as a result of Neo-Tethys oceanic lithosphere subduction beneath the Central Iran zone which is replaced in the Orumieh-Dokhtar magmatic arc.


Introduction
Neshveh intrusion is located in the north west of Saveh and is a small part of the Orumieh-Dokhtar magmatic arc in the Alpine-Himalayan orogenic belt.Volcanic rocks are so common in the Orumieh-Dokhtar magmatic arc.These rocks have varieties in composition and tectonic setting which are varying from acidic to basic and continental to shallow marine environments.There are some intrusions also within the volcanic rocks.Acidic-intermediate volcanic and intrusive rocks are widespread in contrast to basic rocks that are probably indicative for the effect of continental crust in their formation.
Several studies had been carried out in the Saveh region, such as studies of [1] and [2] which were focused on the petrography and Petrochemistry of igneous rocks in the south west of Saveh, volcanic and plutonic rocks in the south of studied area.According to [3], volcanic rocks around the Neshveh village consist of three parts: the oldest volcanic rocks include rhyolite, dacite and some andesite; pyroclastic rocks and lava flows along with intercalation of marl, sandstone and limestone with Lutetian age; and finally andesites and latites with porphyritic textures in the upper part with Priabonian age.The Eocene volcanic flows are widespread in the studied area, which are the result of hybrid crustal-mantle source.The hybrid source of these rocks is attributed to the ascending of basic magma with mantlic source and melting of the continental crust.
Neshveh intrusion located between 35˚10' to 35˚14' latitudes and 50˚10' to 50˚14' longitude at the 40 km north west of Saveh.This intrusion was intruded into the volcanic-sedimentary rocks of Eocene (Figure 1).This intrusion didn't investigate for detail petrology and geochemistry, yet.Also, no reliable analyses of trace elements have been published from this intrusion.This research is a detailed study on the petrology and geochemistry of major and trace elements and tectonomagmatic setting of the Neshveh intrusion.

Geology
Neshveh region which is located in the north west of Saveh is a small part of the Orumieh-Dokhtar magmatic arc within the Central Iran zone.All rocks crops out in this region have Cenozoic age and older rocks don't exposed.Based on Saveh geological map 1:100,000 [4], the main rocks in this area consist of sedimentary, volcanic and intrusive rocks have Paleocene-Eocene age and later (Figure 2).These rocks composed of conglomerate equal to Fajan formation, volcanic-sedimentary rocks of Middle-Upper Eocene consist of intermediateacidic tuffs and lavas, agglomerate and intermediate-basic lavas along with the tuffaceous and sandy limestone and finally Upper Eocene ignimbrite, Lower Red Formation (conglomerate, sandstone, siltstone and red-brown marl), Qom Formation (limestone, marly limestone and marl), Upper Red Formation (conglomerate with red sandstone intercalation), Mio-Pliocene marl, conglomerate equal to Hezar Dareh Formation and quaternary alluvial deposits.
Several granitoid intrusions were intruded into the Eocene volcanic-sedimentary rocks.Dating of some of these intrusions by [1] by K-Ar method on the amphibole, biotite, K-feldspar and whole rock indicate 28 -42 Ma  (Table 1).The Silijard granite, Kalkhab tonalite-diorite, Neshveh granitoid (this study) and Niousht micro-granite are the most important intrusions in the northwest of Saveh.
In addition to mentioned intrusives, some sub-volcanic domes also exist in different parts of the northwest Saveh that have acidic-intermediate composition.These domes intruded into the Upper Eocene-Oligocene intrusives.In this base, these sub-volcanic intrusions attributed to the Post Oligocene.
There are some acidic and basic dykes in this area.Basic dykes with andesite-andesitic basalt composition are widespread in comparison to acidic dike with dacite-andesite composition.These dikes cut the Eocene volcanic rocks and Oligocene intrusives.In this base, their intrusions have been attributed to the post Oligocene.

Research Methods
Our researches in this study consist of two parts: field and laboratory studies.Field studies include identifying the different phases of intrusion, relationship between them and host rock and finally sampling of different phases for laboratory studies.Laboratory investigations include preparing of 90 thin section and petrographic studies, analyzing of 9 samples by ICP-AES and ICP-MS for major, trace and rare earth elements at the ALS Chemex (Table 2).Finally, these analyses were processed by using Excel, Minpet, GCDKit and Igpet programs.

Petrography
Neshveh granitic intrusion with an area about 22 km 2 , consist of intrusive rocks in the mountains of Shahpasand and Zarou, north of Jafarabad Village, north of Aftabroo Qeshlaq village and south of Neshveh village.This intrusion crops out in the east-west trend and intruded into the Upper Eocene volcanic-sedimentary rocks (Figure 3).There are chloritic and argillic alterations present at the contact of intrusion with host rocks.Major faults include Kooshk Nosrat, Neshveh and QermezAghaj with NW-SE trend cut the Neshveh granitic intrusion and caused to crushing and fracturing the intrusion (Figure 3).Moreover, this intrusion cut by a series of andesitic-  basaltic dikes.Some xenoliths of the Eocene volcanic rocks present at the marginal parts of the intrusion, which their dimension and frequency decrease to the interior part.In addition, there are some mafic micro-granular enclaves in different parts of the Neshveh intrusion which are characterized by rounded shapes.These enclaves mineralogically composed of the monzodiorite and diorite.
Based on Petrographic studies of 90 samples from different parts of the Neshveh intrusion, it was subdivided into two compositional zones: quartz monzogabbro in the West and quartz monzodiorite, granodiorite and granite in the East (Figure 4).The contacts are gradational but the rock types are distinct and easily recognized both in outcrop and in aerial photos and satellite images.Compositionally, dikes include andesite-basaltic andesite.Samples Neshveh intrusion has medium to coarse grained granular texture which is locally accompanied by porphyroidic texture with plagioclase phenocrysts, micro-granular and graphic texture.
In the all type of rocks at the studied area, plagioclase present as euhedral and subhedral crystals and are of andesine and oligoclase type.Some of them show zoning (Figure 5).These minerals underwent low grade alte-  ration to sericite and sometimes epidote and calcite (Figure 5).Alkali feldspars (orthoclase) indicate argillic alteration in some samples.At the granitic samples, perthite texture was formed as a result of sub-solidus processes.Quartz presents as fine to medium-grained xenomorphic crystals between the other minerals.Graphic texture was formed as a result of quartz and orthoclase in the granites and occasionally in granodiorites (Figure 5).In some samples, quartz shows wavy extinction indicating tectonic stresses on the intrusion.Mafic minerals in the Neshveh intrusion are mainly hornblende and biotite is less.Hornblende present as euhedral and subhedral crystals at the all phases of Neshveh intrusion (Figure 5).Some crystals of hornblende, altered to actinolite which sometimes accompanied by biotite and iron hydroxides.Biotite present at the all phases of this intrusion which its content increase at the granites.Clinopyroxene is observed with low frequency in the quartz-monzodiorite which mainly has been replaced by actinolite.The rests of the clinopyroxene are found as islands form within the actinolite (Figure 5).Sphene, opaque, zircon and apatite are the minor minerals in the Neshveh intrusion and secondary minerals include sericite, epidote, actinolite, chlorite and clay minerals.As mentioned before, there are some medium-basic dykes within the Neshveh intrusion that compositionally are andesite and basaltic andesite.They have porphyry, microlitic porphyry, inter-granular and intersertal textures phenocrysts of plagioclase, clinopyroxene and hornblende in a fine-grained to microlitic matrix (Figure 5).Some plagioclase crystals show zoning and classified as andesine-labradorite.Some plagioclase crystals altered to sericite and some clinopyroxene crystals show alteration to chlorite.Fine-grained matrix underwent alteration to chlorite.

Geochemistry
For studying of geochemical characteristics of Neshveh intrusion, 8 samples from plutonic rock and one sample from basaltic andesite dikes were analyzed at the ALS Chemex laboratories for major, trace and rare earth elements by ICP-MS and ICP-AES analysis method.The results of this analysis are presented in Table 1.
Based on K 2 O + Na 2 O versus SiO 2 diagram [5] and [6] and the AFM triangular diagram [5], all rocks units of the Neshveh intrusion are calc-alkaline character (Figure 6).At the SiO 2 against K 2 O diagram [7], they are classified as high potassium calc-alkaline rocks (Figure 7).Based on calcic-alkali index [8], all samples Neshveh intrusion have Calc-alkali nature (Figure 7).On the ASI versus A/NK diagram [9], all samples are classified as metaluminous rocks (Figure 7).K 2 O/Na 2 O ratio is more than one in the studied rocks.
On the variation diagrams, with increasing SiO 2 , contents of Al 2 O 3 , P 2 O 5 , CaO, TiO 2 , MgO, MnO and Fe 2 O 3 are decreased and K 2 O and Na 2 O demonstrate increasing trend.These trends may reflect the crystal fractionation process in the evolution of the Neshveh intrusion.The Nd/Sr vs. La/Nb diagram [20], is useful indicator for fractionation of plagioclase in the magma.Based on Figure 8, crystallization and fractionation of plagioclase led to increasing of the Nd/Sr ratio in the magma while the La/Nb ratio will be constant.
[21] exhibited a method to recognition of the three main processes of partial melting, fractional crystallization and magma mixing in the magmas.In this method, an incompatible element is used against a consistent element.According to this method, the consistent element put in the horizontal axis and the incompatible element is in the vertical axis.If the trend of samples is a smooth and straight and almost horizontal line, fractional crystallization is the probable process in the magma evolution.If the trend is a smooth and straight and almost vertical line, partial melting is the effective process in the evolution of magma.In the cases that the trends of samples show a concave curved, this trend indicating that the magma mixing is the most probable process in the evolution of magma.In the V vs. Rb and Ba diagrams (Figure 9), despite the very small scattering, samples of the Neshveh intrusion has relatively horizontal trend indicating that the main mechanism for evolution of the Neshveh intrusion was the fractional crystallization.
Based on geochemical studies, Neshveh granitoid can be classified as sodic granite (Na 2 O > K 2 O) and has low values of Mg#, Ni, Cr, Co and V. Lower contents of these elements indicative for evolution of magma during ascending and before complete crystallization [22].
The primitive mantle normalized spider diagrams for different units of Neshveh intrusion demonstrate negative anomalies of HFSEs including Nb, Ti and P, along with positive anomalies in LILEs including Rb, Ba, Th, K and Pb (Figure 10).The negative anomalies of Ti and Nb are controlled by titanium bearing minerals such as titanite, ilmenite, rutile, garnet and some amphiboles.With increasing pressure, the solubility of titanium bearing minerals reduces in aqueous fluids [24].Therefore, theses minerals that are rich in HFSEs will be stable during partial melting process at depths greater than 30 km that is led to negative anomalies in the melt [24].The negative anomaly in P indicative for low content of apatite in the studied rocks.As mentioned before, in the I-type granitoids, phosphorus act as a consistent element and fractionation this element in the early stages of magma-   crystallization, led to negative anomalies in the crystalized rocks.The sample from dikes has negative Th anomalies (Figure 9).Based on mineralogical and geochemical composition of dykes; it seems that their parental magma source is the mantle.Since the Th content in the mantle is higher than the crust, this negative anomaly indicating that this magma rise up rapidly and underwent very low crustal contamination.
Fluids and melts derived from the subducted oceanic slab, led to the metasomatism of upper mantle wedge and thereupon, negative anomalies of Nb and Ta [25]- [27].Positive anomalies in Pb, K, Rb and Ba are attributed to metasomatism of mantle wedge by fluids derived from the subducted slab or/and contamination with continental crust [28] [29].
The Chondrite-normalized REE patterns in the different phases of Neshveh intrusion show a LREE rich pattern with a high LREE/HREE ration (Figure 10).All phases of Neshveh intrusion indicate negative anomalies Eu while sample from the dikes show a positive anomaly in Eu (Figure 10).Negative anomalies in Eu are related to the feldspar fractionation during magma crystallization or remaining of feldspar at the source.If these negative anomalies are associated with negative anomalies in Sr, fractional crystallization of plagioclase is re-sponsible for these anomalies.While, the negative anomalies in Eu are associated with negative anomalies in Ba, K-feldspar fractionation has the main role in this respect [25].In the studied rocks, negative anomalies of Eu occur along with Sr, which refers to the plagioclase fractionation.Positive anomaly in Eu from the dike sample is indicative for the absence of plagioclase fractionation during the magma crystallization and high concentration of this mineral in these rocks.

Tectonomagmatic Setting
Besides the geochemical characteristics of rare elements such as spider diagrams, tectonomagmatic setting discrimination diagrams also demonstrate a subduction setting relation for the Neshveh granitoid.
Based on the [30] diagram that is offered with considering, two cationic parameters include R2 = 6Ca + 2Mg + Al and R1 = 4Si − 11(Na + K) -2(Fe + Ti), samples of Neshveh intrusion put into the pre-plate collision setting (Figure 11).On the [31] triangular diagrams that are drowned based on immobile elements include Ta and Hf and mobile element of Rb, samples of studied intrusion put into the volcanic arc (VA) setting (Figure 11) Based on the Rb versus Y + Nb and Rb versus Ta + Yb diagrams [32] samples from the Neshveh intrusion put in the volcanic arc granite (VAG) setting (Figure 12).On the Th/Hf versus Ta/Hf, Th/Ta versus Yb and Ta/Yb versus Th/Ta diagrams [33], all samples from the Neshveh intrusion put into the active continental margins are (Figure 13).The sample from the dikes is situated in the within plate setting (Figure 13).
Based on the AFM diagram [34], the studied rocks located in the continental arc granites (CAG) setting (Figure 14), indicating that these granitoids originated from the subduction of an oceanic crust under a continental crust at the active continental margins.The parental magma was medium to high K calc-alkaline.
The Th/Yb versus Ta/Yb diagram [35] was used to recognition the source magma for Neshveh intrusion.According to this diagram, the studied samples are located within the calc-alkaline Andean granitoids (Figure 15).Moreover, the variation trend of samples on this diagram indicates the evolution of magma by fractional crystallization.As mentioned in the previous parts, based on [19] classification, Neshveh intrusion can be classified as ACG granitoids.According to him, ACG granitoids are hybrid type (mantle + crust) and are related with subduction settings.In this base, it is supposed a subduction setting to Neshveh intrusion.

Conclusion
Field studies, petrographic and geochemical investigations indicate that the Neshveh granitoid consists of the quartz-monzodiorite, granodiorite and granite.Field relations and geochemical studies demonstrate that all of these rocks are related to a single parent magma which underwent fractional crystallization.The Neshveh granitic intrusion has high-K calc-alkaline nature.This intrusion is metaluminous and classified as I-type granitoids.Negative anomalies of Nb and Ti along with positive anomalies of Rb, Ba, K, and Pb in the spider diagrams are indicative for the formation of this intrusion in relation with subduction setting.In terms of geodynamic setting, this intrusion is classified as volcanic arc granites (VAG) and active continental margin granites (CAG) which have been formed as a result of Neotethys oceanic crust subduction beneath the Central Iran continental crust.

Figure 1 .
Figure 1.Geographic location of the Neshveh intrusion in the northwest Saveh.

Figure 3 .
Figure 3. (a) View from the Neshveh intrusion and its contact with Eocene volcanic-sedimentary rocks (view to the south); (b) Close view from the fracturing in the Neshveh intrusion as a result of tectonic.

Figure 4 .
Figure 4. Simplified geological map of NW Saveh intrusions (Based on satellite data and the geological map of Saveh in 1:100,000 scale; Galamgash and Fonudi, 1998).

Figure 9 .
Figure 9. (A) Increasing trend in the Nd/Sr versus La/Nb which is indicative for plagioclase fractionation; (B) and (C) Decreasing trend of Rb and Ba versus V which is indicative for fractional crystallization on the evolution of Neshveh intrusion.

Figure 11 .
Figure 11.(A)Location of the Neshveh intrusion samples on the two cationic parameters of R1 and R2 [29].(B) and (C) Location of the Neshveh intrusion samples on the Hf-Ta-Rb triangular diagrams [30].

Figure 14 .
Figure 14.Location of the studied samples on the AFM diagram [33].

Figure 15 .
Figure 15.Location of the Neshveh intrusion sample on the Th/ Yb versus Ta/Yb diagram [34].

Table 1 .
[1] result of dating analysis on the some intrusions around the Saveh[1].

Table 2 .
The result of samples analysis from the studied area by ICP-AES and ICP-MS methods (major elements in % and trace elements in ppm).