The studied area is located in the 1:100,000-scale geological map of Songour and the west of Iran )Figure 1).
The some parts of the stratigraphic units in the studied area belong to Sanandaj- Sirjan structural zone and other parts are related to Zagros Thrust Zone, which marked by the large fractures and multiple thrusts, and resulted to structural complexity of the area (Figure 2). This area, by the following of geological events in the area, was under the influence of different deformation phases whether orogenic or epeirogenic, the emplacement of intrusive bodies, metemorphism, and other geological events.
Figure 1. Map of the study population  .
Figure 2. Structural and sedimentary zones of Iran  .
Given that the Sanandaj-Sirjan zone having NW-SE trending is parallel to the Zagros Thrust Faults, and by considering that this zone has severely affected by these faults activities in during different geological periods after the Cretaceous, as a result, in the area investigated the severe effects of tectonic activities could be well observed due to the proximity to the Main Zagros Fault. Among these features, it can be referred to thrust nappes, the presence of thrust faults with Quaternary age, regional metamorphism and intrusion of gabbroic to syenitic bodied and a small amount of granite.
The rocks cropping out in the studied area belong to the Mesozoic and Cenozoic era and the oldest of them is related to the Upper Cretaceous. These intrusions are Tertiary in age and have been injected in the upper Eocene-Oligocene due to the Pyrenean orogenic phase.
This study aims to determine the formation mechanism of gabbros, syenites and granites, discover the genetic relationship among the intrusive bodies of the region, and determine the change and frequency of major and trace elements in the rocks of the region.
2. Geological and Tectonic Setting
The area studied with coordinates of 47˚30' to 47˚40' east longitude and 34˚30' to 34˚40' north latitude is located in the southern part of 1:100,000 geological sheet of Songour and mountainous area with rough topography, and covers an area of over 280 square kilometers. The area of the intrusions outcrops in this region is about 25 square kilometers (Figures 3-5). Many faults in this area are of thrust or high angle reverse type, but there are less normal and strike-slip faults. One of the characteristics of the faults in the area is their seismicity so that until now several major earthquakes have been occurred in the region and caused the destruction of houses. Faults in the area are the valley-building. This means that the valleys and waterways have been created in trend of some faults due to slip, displacement and crushing caused by faulting and subsequent erosion.
Figure 3. Gabbroic rocks of the area.
Figure 4. An outcrop of syenite rock.
Figure 5. Selected microscopic images of rock groups studied (in crosspolarizedlight, XPL) (a) Hornblend (b) Augite (c) Quartz recrystallization (d) Polysynthetic Plagioclase (e) Quartz and (f) Perthite Orthose (XPL. 40×).
In the late Eocene and probably during the Pyrenean orogeny  , governing compressive forces on the region along with other geological processes have caused to intrusion and emplacement of intrusive bodies in the studied area. Under the influence of the same orogenic phase, the contact metamorphism as well as metasomatism resulting from emplacement of intrusion have been occurred and caused the formation of skarn units in Deh Asyab and the position to form in the marble limestones in the studied area as well. This phase also resulted in folding and possibly faulting of the older units.
According to 19 different type of rocks using ICP-MS and alkaline fusion technique the analysis (Table 1 and Table 2) and the geochemical diagrams of rocks nomenclature  and  , most samples in the studied area would be fallen in the range of the gabbro,
Table 1. Chemical analysis of intrusive pluton in south of Songour.
In the diagram AFM  , the samples would be fallen in both of tholeiitic and calc- alkaline range. In the diagrams Na2O/SiO2  , and K2O/SiO2  , most samples are
Table 2. Chemical analysis of intrusive pluton in south of Songour.
plotted in alkaline and sub-alkaline basalt fields. These diagrams represent an almost disorder linear trend for samples in the area studied. These disorders could be resulted from the crustal contamination of some intrusive magmas and increase in the amount of K2O within them, which of course does not reject the same origin of all magmas of
Figure 6. Determination of intrusive rocks type diagram  .
Figure 7. Determination of intrusive rocks type diagram  .
the area (Figures 8-10).
Based on the aluminum saturation diagram  , 18 samples are in metaluminous and peralkaline fields. One example is in field of peraluminous, indicting the igneous origin of all the samples (Figure 11).
Figure 8. Determination of magmatic series diagram using the major element oxides data  .
Figure 9. Na2O/SiO2 diagram to determine magmatic series  .
Figure 10. K2O/SiO2 diagram to determine magmatic series  .
Figure 11. Aluminum saturation diagram to determine magmatic series  .
Based on the major element oxides vs. SiO2  , the oxides Al2O3, MgO, CaO, TiO2, P2O5 and FeO decrease with increasing silica content and being more acidic of rocks, and K2O Na2O oxides show increasing (Figure 12).
Figure 12. Variation plots of the major element oxides vs. silica  .
Based on the magnesium and calcium oxidesvs. silica diagram, the two almost regular trends it could be inferenced so that the basic to intermediate samples and the intermediate and acidic samples would be fallen follow each other, which this relationship between the two series of basicand intermediate magmas, on the one hand, and on the other hand between intermediate and acidic magmas can be deduced, and this would be confirmed by field geological studies so that the basic to intermediate rocks (gabbros and diorites) present with each other and the intermediate to acidic rocks (monzonite, syenite, quartz syenite and granite) could also be found near each other, and a distinct gradual variation among them can be understandable.
On the diagram Na2O vs. SiO2, the two trends can be observed. The increasing Na2O from basic rocks to high-silica acidic rocks (about to 63%) and decreasing trend from rocks with more than 66% silica toward more acidic rocks, which could be attributed to Na2O depletion in magma and its use in the formation of sodium-bearing feldspars.
According to the rare earth elements vs. SiO2 diagrams, the elements rubidium (Rb), yttrium (Y), zirconium (Zr), selenium (Se), nickel (Ni), thorium (Th) and niobium (Nb) increase with increasing SiO2, and the elements strontium (Sr), barium (Ba), chromium (Cr) and cobalt (Co) show decreasing trend (Figure 13).
Figure 13. Variation plots of the trace elements vs. silica  .
The value of Rb increases with increasing SiO2 in the magma, due to consumption of Rb alkali feldspar’s structure. A dispersion could be observed on the diagram Ba vs. SiO2, suggesting the process of partial melting, assimilation and magmatic contamination in the region. Yttrium (Y) is one of the rare earth elements, which concentrates in acidic terms particularly pegmatites and display a systematic increasing trend from basic to acidic rocks. This systematic increasing could be indicating the same source of magmas in the studied area. Dispersion could be seen on the variation plot of Co vs. SiO2, of course, higher Co content would be consumed within initial Mg-bearing minerals especially olivine. The scattered points on the diagrams are another evidence for the existence of the assimilation and magma mixing or only one of these two processes. Cr would be removed from the magma as chromite mineral in the early stages of crystallization. Also, Cr would be accumulated in pyroxenes, especially the ultra-basic pyroxenes. Based on the diagram Cr vs. SiO2, Cr has more concentrated in four samples (normal) and two syenite samples (abnormal). This abnormal concentration could be resulted from the assimilation and magma mixing processes, of course, at the small- scale, because has created anomaly only in two samples.
According to petrographic and mineralogical studies of thin sections of rocks collected from the studied area, the intrusive igneous rocks in the region from basic to acidic terms include gabbro, gabbro-diorite, diorite, quartz-bearing syenodiorite, syenite, quartz-bearing syenite and granite, which gabbroic and syenitic rocks contain the largest volume in the area studied. The most studied samples are saturated with silica and in terms of the aluminum saturation are peralkaline and metaluminous. Generally, almost all the lithological and geochemical diagrams indicate a single source for these samples. It appears that the basaltic magma originated from the mantle, has intruded into the crust by weakness point of the crust namely the Zagros Main Thrust located in the south part of intrusive bodies. This magma has slightly fractionated in situ during solidification and this resulted in the formation of rocks from gabbro to syenite and granite, which of course, from the time of stopping until the magma has been cooled and solidified, the main basaltic intrusive (confirmed by large volume of gabbros) was able to contaminate with the crust and led to increase in some elements, such as potassium. In the Eocene and Oligocene magmatism of the Zagros Thrust and Sanandaj- Sirjan zones including the studied area, the cause of high compression by tensile and opening forces is acceptable.