The central Iranian volcano-plotonic belt (CIVB), composed on tertiary volcanic and plutonic rock that shows Eocene early Oligocene, middle miocen, Pliocene and quaternary magmatism  . The urge Eocene volcanism    composed of basic to acidic calcalkaline magmatic activity in subduction zone   . The Oligocene and Miocene plutonic rocks in CIVB introduced to Eocene volcanic rocks  . The CIVB is an active volcanic zone in upper cretaceous and Eocene  . The CIVB composed of mafic-felsic volcanic and plutonic igneous rocks which are introduced by subvolcanic Oligocene Miocene dykes and stocks  -  . The country rocks of porphyry copper deposits in CIVB are subvolcano plutonic rocks (Granodiorite, diorite, monzonite, quartz onzonite, …) occurred in subduction zones and similar to Andian copper deposits     .
The CIVB has situated between Sanandaj-Sirjan and central Iranian zone. The salafchegan studied area is located at North and North-west of CIVB and shows adakitic and alkaline sequence of colligenal and post colligenal magmas in Pliocene and Quaternary  . The magmatic activity in Salafchegan similar to CIVB, is due to subduction of Arabian plate and Neogene magmatism in cretaceous and Oligocene periods  -  (Figure 1). The host rock of most purified copper which has the position of magmatic arc is shallow intrusive in CIVB. These intrusive have been formed in a magmatic arc and convergent plate margins, which can also be seen in copper mine in Andes     . The history of magmatism area in Salafchegan in central Iran represents the subduction of Tethys ocean lithosphere in Eocene below the continent of central Iran which is associated with the rift phenomenon of the Red Sea, Gulf of Aden    .
After the Oligocene, the rate of convergence increases between Arabic plate and continent of central Iran. The ocean began to be closed and collisions were occurred       . Consequently, in Salafchegan region, a number of dextral strike-slip faults with direction of North West-South East were formed as a result of oblique convergence   .
Traditionally, this area has attracted the attention of geologists by various magmatic and mineral phenomena. In geological map of Emami & Hajian, 1991, Qom 1,250,000  ; igneous masses of region have been studied. Also in geological map of Tafresh (1999), Salafchegan-Khorhe (2000), Alaee Mahabadi and Hajian   have studied the mentioned masses.
Figure 1. Positions of study area in CIVB magmatic belt (Mirnejad et al., 2013).
Emami  was studied the petrology of Khastak Caldera in Zavarian region in his doctoral thesis titled Geology of Qom-Aran region. Monsef (2001)  studied the geochemistry and petrogenesis of Neogene volcanic and semi-volcanic rocks in West of Salafchegan up to North of Delijan. In this study, we attempt to evaluate geochemical features, formation and origin of exterior igneous rocks of the region according to the field data, petrographic studies, and the results of chemical analysis of exposed volcanic rocks.
2. Research Method
During the Field Operations, 100 different rock samples were selected with chip sampling method. Samples with minimum weight of 1kg placed inside of plastic bags. In order to identify lithological composition and crosscutting relationships, thin microscopic section were prepared and studied. In order geochemical studies, 29 volcanic rocks samples as much as possible unaltered and minimum weathering and secondary streak which covered characteristics of the area rocks were selected and transferred to Geological Survey & Mineral Explorations of Iran (GSI) laboratory for XRF and ICP-MS analysis. The limit of oxygen index (LOI) was also measured.
3. Results and Discussion
3.1. Region General Geology
The study area is located in Salafchegan district in 40 km South West of Qom and near the village of Zvaryan (Figure 1). The study main outcrop was in 34˚15' to 34˚40' North latitude and 50˚15' to 50˚28' East latitude and it is located in an area of about 200 square kilometers.
The Salafchegan region is located in CIVB magmatic zone with direction of North West-South East and composed of Tertiary volcanic intrusive rocks with a length of 1700 km and an average width of about 150 km. there weren’t any outcrops for Precambrian basement in the region and the oldest rock units in Salafchegan area, was sedimentary rocks in Jurassic era (Figure 2). The magmatic activity in this region started in early Eocene and continued until late Eocene  . This volcanic complex in the hole of region consists of basaltic to rhyolitic lavas and associated pyroclastic with sedimentary rocks among these layers (Figure 2). The Paleogene volcano-sedimentary sequence has been covered discontinuously with unconformity layers of red and calcareous units  . These sequences invaded by granitoid masses with calc-alkaline nature. Shallow intrusive activity continued until the late Miocene and tectonic settings were occurred in Pliocene after the collision. Within the volcanic masses of the area is seen the andesitic dikes to be seen which has cut the volcanic masses that indicates extensive stretching at the end of the volcanic stage.
In addition to nomination of initial modal and characterization of Mineralogical and textural thin sections, analysis of events affecting the formation of minerals
Figure 2. Geology map of study area the resulting Merj (Taheri, 2016) of 1:100,000 Salafchegan-Khorhe geology map (Alaee Mahabadi, 2000) and Tafresh (Hajian, 1999) and field observations.
and identification of texture, mineral relationship (relationship of minerals), reactive margins, Xenoliths and other notable features of were evaluated. Igneous rocks of studied area were divided in two main intrusive and extrusive categories. The volcanic rocks include: olivine basalt, basalt, andesite basalt, andesite, trachyandesite, quartz Latite, Latite, trachyte and dacite and member of intrusive igneous rocks include: diorite, monzonite and quartz monzonite. Based on main purpose of this article, the main volcanic rocks of region are described below.
3.2.1. Basaltic Lavas
These rock groups consist of Olivine basalt, basalt, andesite basalt. In terms of texture, composed of porphyric texture with microcrystalline pulp (Figure 3) and phenocrysts include of plagioclase, olivine and clinopyroxene. In some samples chlorite and opacities amphibole were observed. Alkali feldspar and quartz can be seen as a rarely sub-grade mineralization. Pennine mineral as result of serve chloritisation were seen in some samples (Figure 4). OPEC minerals were present primarily and secondarily in rocks and rarely tiny cavities which filled with calcite were observed. In some cases, olivine has been serpantnised and formed skeletal structure. Olivine idnegization phenomena startss from margins and fractures and continues to the center of crystals (Figure 5).
3.2.2. Andesite, Latite, Trachite Lavas
The most extensive outcrops of volcano rocks in this region are andesitic lavas. These rocks consists of porphyric texture with microlitic pulp stream or microcrystalline, which in some cases formed as a micro porphyric amygdaloid (Figure 6). Disequilibrium textures such as sieve texture (Figure 7) and embayment in these crystals are rarely seen. Presence of sieve texture in plagioclase phenocrysts of andesitic lavas can be due to an imbalance in the magmatic system. Also,
Figure 3. Porphyric texture with microcrystalline pulp (field length: 4.5 mm, 40X) (XPL).
Figure 4. Presence of pennine minerals as result of serve chloritisation in basaltic andesite (field length: 4.5 mm, 40X) (XPL).
Figure 5. Skeletal structure as result of epastrization of grains of olivine in basalt rocks with serpantinization phenomena (field length: 4.5 mm, 40X) (XPL).
Figure 6. Micro-porphyric amigdual texture with cavities filled by calcite (field length: 4.5 mm, 40X) (XPL).
Figure 7. Sieve texture of plagioclase in andesitic rocks (field length: 4.5 mm, 40X) (XPL).
sieve texture in some cases indicates magmatic mixing and sudden fall of pressure and shows unbalanced conditions in magmatic reservoir in general.
Andesitic rocks is consists of fine microliters-oriented plagioclase with some glass between them. Plagioclase which has zoning is the main mineral in these rocks (Figure 8). Twill twinning (Figure 9), crusader twinning and pericline twinning (Figure 10) can be seen in these rocks. Coarse crystals are seen in different sizes in the fine pulps consists of plagioclase, clinopyroxene, burned hornblende (Figure 11), opaque and glass minerals and in some cases alkali feldspar (sanidine) and rarely quartz. Hornblende type of amphibole are mainly altered and replaced with chlorite, calcite and opaque minerals. Presence of rim reaction in amphibole crystals depends on rate of ascent of magma. Secondary minerals include sericite, epidote, silica, iron oxide and hydroxide mineral.
Phenocrysts are visible in Latite of plagioclase, potassium feldspar, approximately in the same values in the crystallized in fine to glasses. Plagioclase is the most abundant Connie Stone which is largely altered. Amphibole macro crystals were seen as automorphe and sub-automorphe with burned margins are often seen as chloritisate, opacitiesate, epidotisate and carbonatizate. The dominant
Figure 8. Zoning of plagioclase phenocrysts in the andesite (field length: 4.5 mm, 40X) (XPL).
Figure 9. Twill twining of plagioclase in andesite (field length: 4.5 mm, 40X) (XPL).
Figure 10. Pericline twinning of plagioclase in andesite (field length: 1.8 mm, 10X) (XPL).
Figure 11. Amphibole phenocrysts in the andesite with burned margins (field length: 4.5 mm, 40X) (XPL).
process in region hornblendes is apasytisation of hornblendes. In some samples, hornblendes have maintained their hexagonal or diamond shape but completely oxidized. This phenomenon is result of formaion OPEC minerals which seems to be magnetite and hematite. Alkali feldspar is formed as macro crystals up to 40% of phenocrysts. Minor minerals consist of epidote, chlorite, epic minerals and quartz which includes about 2 percent of pulp. Micro enclaves with granitic composition can be seen in a number of andesitic sections (Figure 12).
3.2.3. Dacitic Lavas
Plagioclase, alkali feldspar, quartz amphibole have created the mineralogical composition of the rocks. Plagioclase as the most abundant phenocrysts has polysynthetic twinning and clear zoning combination. Presence of zoning combination and consequently regional blackouts represent complex process in the liquid-crystal interface during rapid cooling in lower temperature  . In the other hand, the zoning of plagioclase indicates an imbalance during the crystallization and solidification which can be the result of changes in crystal local conditions like thermal turbulence caused by the entry fresh magma in crystallization of magma chamber, rapid decrease in temperature and changes in the water vapor pressure and dissolved gases during crystallization of plagioclase in the magma  .
Amphibole, is often seen as shaped crystals (Figure 13) to hypidiomorphic in fine pulps. In some amphibole crystals, burned margins are observed. In most cases these burned margins are the result of oxidation reactions which are caused by heat from sudden pressure drop in magma during ascent.
Alkali feldspar is another minerals of these rocks which seen mostly as fine crystals in stone pulp and it is difficult to be diagnosed. Epic minerals were seen primarily and secondarily in stone pulp. Secondary epic minerals mostly formed as result of hornblende alteration and they are seen in the margin and center of minerals.
Porphyric texture with micro crystal pulp (Figure 14), is the most abundant texture of region dacite rocks. Plagioclase and alkali feldspar, shaped feldspars with eroded margins, are the most important coarse crystals. In the other hand, plagioclase coarse crystals were seen with swinging zoning, which it is the sign of crystallization in unstable physicochemical conditions. Erosion of a solid phase in a fluid phase is a sign of fluid saturation which means that the initial
Figure 12. Granite micro enclaves in the andesitic rocks (field length: 4.5 mm, 40X) (XPL).
Figure 13. Amphibole phenocrysts in dacite (field length: 1.8 mm, 100X) (XPL).
phase, which has crystallized in fluid due to saturation, is undergoing liquidation due to changing circumstances. Sudden change in temperature, pressure, oxygen fugacity, vapor pressure, digestion stones stuck inside or mixing with other liquid cause changing in the composition of the fluid  .
3.3.1. Classification Extrusive Igneous Rocks and Magma
In this study we used diagram presented by Cox et al. diagram  for nomination of volcanic rocks. According to the results the rocks of region are within the basalt, andesite basalt, andesite and dacite (Figure 15) with sub-alkaline nature and acid intermediate composition. According to the results of geochemical analysis of samples (Table 1), these samples based on the Winchester & Floyd
Figure 14. Porphyric texture in dacite (field length: 4.5 mm, 40X) (XPL).
Figure 15. Location of studied volcanic rocks in volcanic rocks nomination charts based on TAS method (Cox et al., 1979).
Table 1. The result of samples analyses from the studied area by XRF and ICP-OES methods (major elements in % and trace elements in ppm).
charts  , are within andesite, dacite and rhyodacite (Figure 16). For determining the magmatic series of volcanic rocks, Irvina & Baragar charts  were used and samples were within sub-alcaline, calc-alkaline and rarely tholeiitic because of weathering and lose of potassium oxide (Figure 17 and Figure 18). According to the Luba charts  (Figure 19) region rocks are sub-alcaline and
Figure 16. Location of region samples in volcanic rocks nomination charts by using SiO2, Zr/TiO2.
Figure 17. FeO* − (Na2O + K2O) − MgO for separation of tholeiitic and calc-alkaline range (Irvina & Baragar, 1971).
Figure 18. Samples location in the volcanic rocks segmentation chart to alkaline and sub-alkaline according to the total alkali versus silica (Irvina & Baragar, 1971).
Figure 19. Location of studied samples in volcanic rocks nomination chart according to the TAS method (Le Bas, 1986).
within andesite, basaltic andesite basaltic and dacite. According to the K2O-SiO2 (Figure 20)  most of the region volcanic rocks are in intermediate groups and fewer in acid range. All of the samples are in low-potassium tholeiitic range.
3.3.2. Magmatic Tectonic Pattern
For recognizing tectonic environment of volcanic rocks within studied area different charts were used. According to the Muller & Groves (1977)  charts all samples are in range of volcanic rocks of volcanic arcs after collision (Figure 21 and Figure 22). Moein Vaziri  believes that Neotethys Ocean due to movement of afro-arabian plate towards the northeast has subducted (drived) to the
Figure 20. Changes in SiO2, K2O in volcanic rocks of studied region (Muller et al., 1998).
Figure 21. Y/Zr diagram for determination of tectonic setting of volcanic rocks of studied area (Muller and Groves, 1977).
Figure 22. Zr/TiO2-Ce/P2O5 diagtam (Muller and Groves, 1977) for separation of active continental margin magmatic arc (CAP) and magmatic arc after hitting (PAP).
continent of Central Iran. The matching of Subduction zones and dealing with old hydrothermal blades which had created Mesozoic rift within the continent, could cause disturbance in classical magmatism in mentioned tectonic setting. A break of Neotethys oceanic crust in the upper Cretaceous allows that depressed and stopped part of the shell get uniformly warm (hot) and causes the creation of severe volcanoes in the Eocene period. These phenomena continued with less intensity in other parts of the Tertiary  . Some researchers   believed that most intrusive and semi-volcanic masses of Iran’s tertiary are characterized by Oligo-Miocene ages which formed as result of pressure relief deep magmatic reservoirs of Eocene volcano in orogenic period of Pyrenees and astrains era  .
In Figure 23, distribution pattern of rare earth elements (REE) in region samples has been ruled against chondrites  . According to the diagrams, all samples are richer than chondrites in REE. Among these, Light rare earth elements (LREE) are richer than High rare earth elements (HREE). As seen in figure, behavior of Eu is almost like other REE and do not represent any particular anomaly. Rollinson,  mentioned that, this phenomena happens when fugacity of oxygen is high and distribution coefficient of Eu is low. In spider diagram normative to MORB  , region rocks show somewhat depletion in HFSE elements (Figure 24), While enriched from K, Th and Rb elements which has features of magmatic arc rocks or volcanic arc rocks  . These elements are special for crustal rocks and may indicate the digestion of crustal rocks by the initial molten material.
In subduction area, fluid released from the upper part of subducted lithosphere which is rich of LILE and poor of Nb are added to the mantle wedge  .
Figure 23. The normative pattern of rare earth elements compared to chondrites in studied samples (Boyton, 1984).
Figure 24. Spider diagram normalized to MORB in studied area (Pearce, 1983).
The high amount of Ce, K, Rb and Ba elements in these rocks can indicate: 1) Digestion continental crust rocks in the direction of the upper mantle and molten rock is swept; 2) The involvement of magmatic origin of continental crust (partial melting) or deeper magmatic origin (magma mixing); 3) The effect of alkaline metasomatism on the area rocks during or after crystallization and positioning the rocks and 4) Extensive hydrothermal alteration in area rocks (kaolinisation with quartz-sericite alteration) which is a result of magmatic fluids rich shell elements.
In the primitive mantle normalized charts (Figure 25) Enrichment from LREE and HREE were observed in all samples. Enrichment of LREE elements is higher than HREE and negative gradient of the graph indicates the higher amount of LREE than HREE that conforms with the overall pattern of spider diagrams subduction zones. Defant & Drummond  mentioned partial enrichment of P in basaltic and andesite samples caused due to presence high amount of apatite in these samples. Provatke and Klemme  believed that presence of apatite in rocks in addition to creating a positive anomaly in P can cause Sr enrichment in rocks. Sr plays as a highly consistent element in apatite. Positive anomaly of Rb May be happened due to contamination magma with crustal (due to the high concentration of these elements in the continental crust). Because of the continental crust enriched with Nb and Ta, so suspension of these rocks magmas (although short-term), in the lower continental crust with penetration of crustal fluids into the magma or by digestion of crustal magma rocks cause this type of anomaly in studied samples   .
In conjunction with enriched potassium relative to primitive mantle, Watson  believed that contamination of magma with continental crust, Even if the other elements remain unchanged, significant amounts of potassium finds a way into the basaltic magma. Also the positive anomaly of Sr in basaltic rocks is justified by the presence in these rocks plagioclase  . Similar and parallel trends drawn in all spider diagrams can introduce a single source for these rocks.
Meanwhile, negative anomalies of Ti element in lava can be formed one hand to participate in the partial melting of lithospheric mantle linked metasomatized and on the other hand crustal contamination resulting from digestion are attributed zinocryst  .
Using SiO2 alternation plot against Sr/Y  (Figure 26) volcanic rocks of studied area are not located in adakitic rocks groups. Also the Y against Sr/Y diagram  (Figure 27) indicate location samples within the volcanic arc.
Figure 25. Spider diagram normalized to primitive mantle (Sun & McDonough, 1989).
Figure 26. SiO2/Sr/Y diagram (Yofodzinski et al., 2001).
Figure 27. Y alternation against Sr/Y (Defant & Drommond, 1990).
Lithological composition of volcanic rocks in studied area consists lava basalt, andesite and dacite. The most abundant volcanic rock outcrop is andesite. The dominant texture of these rocks is porphyric with microcrystalline pulp. The magmatic nature of sub-alcaline volcanic rocks is generally alkaline and less tholeiitic. The region rocks enriched with LILE and HFSE elements were less depleted in comparison to MORB which depends on volcanic arcs. Volcanic rocks of studied area are located in tectonic setting of volcanic arc after the collision (PAP). Also presence of negative anomaly in P, Nb, Ta and Ti elements could indicate the formation of rocks in subduction zone. According to the Zr/Y diagram, all samples are located in extent of volcanic arc. Volcanic rocks of region are not in adactic rock groups. The tectonic evidence of region shows that Neogene magmatism of the region is connected with a collision area with Tethys Ocean subduction under Magmatic Belt of CIVB and collision between the Arabic plate central Iran during the Neogene era.
 Emami, M.H. (1981) Géologie de la région de Qom-Aran (Iran): Contribution a l’étude dynamique et géochimique du volcanisme Te rtiaire de l’I-ran Central. Ph.D. Thèse, Sciences naturelles Univ. Sc. Et Medicale de Grenoble., 489 p.
 Jung, D., Kursten, M. and Tarakian, M. (1976) Post-Mesozoic Volcanism in Iran and Its Relation to the Sub-Duction of the Afro-Arabian under the Eurasian Plate. In: Pilger, A. and Rosler, A., Eds., Afar between Continental and Oceanic Rifting, Volume II, Stuttgart, Schweizerbatsche Verlagsbuch-Handlung, 175-181.
 Berberian, F., Muir, I.D., Pankhurst, R.J. and Berberian, M. (1982) Late Cretaceous and Early Miocene Andean-Type Plutonic Activity in Northern Makran and Central Iran. Journal of the Geological Society, 139, 605-614.
 Boomeri, M., Nakashima, K. and Lentz, D.R. (2010) The Sarcheshmeh Porphyry Copper Deposit, Kerman, Iran: A Mineralogical Analysis of the Igneous Rocks and Alteration Zones Including Halogen Element Systematics Related to Cu Mineralization Processes. Ore Geology Reviews, 38, 367-381.
 Hosseini, M., Zarasvandi, A.R. and Liaghat, S. (2010) Tectonic Setting and Petrogenesis of Darreh-Zereshk Granitoids (SW Yazd, Iran) and Comparison with Some World Skarn Granitoids. Abstract for the 1st International Applied Geological Congress, Department of Geology, Islamic Azad University-Mashad Branch, Iran, 26-28 April 2010.
 Barzegar, H. (2007) Geology, Petrology and Geochemical Characteristics of Alteration Zones within the Seridune Prospect, Kerman, Iran. Ph.D. Thesis, Georessourcen und Material technik der Rheinisch Westfalischen Technischen Hochschule Aachen, 202 p.
 Kirkham, R.V. and Dunne, K.P.E. (2000) World Distribution of Porphyry, Porphyry-Associated Skarn, and Bulk-Tonnage Epithermal Deposits and Deposits. Geological Survey of Canada Open File Report 3792a, 26 p.
 McInnes, B.I.A., Evans, N.J., Belousova, E., Griffin, W.T. and Andrew, R.L. (2003) Timing of Mineralization and Exhumation Processes at the Sarcheshmeh and Meiduk Porphyry Cu Deposits, Kerman Belt, Iran: Mineral Exploration and Sustainable Development. Millpress, Rotterdam, 1197-1200.
 Sillitoe, R.H. (1997) Characteristics and Controls of the Largest Porphyry Copper-Gold and Epithermal Gold Deposits in the Circum-Pacific Region. Australian Journal of Earth Sciences, 44, 373-388.
 Ghadami, G., Shahre-Babaki, A.M. and Mortazavi, M. (2008) Post-Collisional Plio-Pleistocene Adakitic Volcanism in Central Iranian Volcanic Belt: Geochemical and Geodynamic Implications. Journal of Sciences Islamic Republic of Iran, 19, 223-235.
 Ricou, L.E. (1994) Tethys Reconstructed: Plates Continental Fragments and Their Boundaries Since 260 Ma from Central America to South-Eastern Asia. Geodinamica Acta, 7, 169-218.
 Mohajjel, M., Fergusson, C.L. and Sahandi, M.R. (2003) Cretaceous-Tertiary Convergence and Continental Collision, Sanandaj-Sirjan Zone, Western Iran. Journal of Asian Earth Sciences, 21, 397-412.
 Agard, P., Omrani, J., Jolivet, L. and Mouthereau, F. (2005) Convergence History across Zagros (Iran): Constraints from Collisional and Earlier Deformation. International Journal of Earth Sciences, 94, 401-419.
 McClusky, S., Balassanian, S., Barka, A., Demir, C., Ergintav, S., Georgiev, I., Gurkan, O., Hamburger, M., Hurst, K., Kahle, H., Kastens, K., Kekelidze, G., King, R., Kotzev, V., Lenk, O., Mahmoud, S., Mishin, A., Nadariya, M., Ouzounis, A., Paradissis, D., Peter, Y., Prilepin, M., Reilinger, R., Sanli, I., Seeger, H., Tealeb, A., Toksoz, M.N. and Veis, G. (2003) Global Positioning System Constraints on Plate Kinematics and Dynamics in the Eastern Mediterranean and Caucasus. Journal of Geophysical Research: Solid Earth, 105, 2695-5719.
 Fowler, A.D. and Ach, D.E. (1996) A Model and Simulation of Branching Mineral Growth from Cooling Contacts and Glasses. Mineralogical Magazine, 60, 595-601.
 Le Bas, M.J., Le Maitre, R.W., Streckeisen, A. and Zanettin, B. (1986) A Chemical Classification of Volcanic Rocks Based on the Total Alkali-Silica Diagram. Journal of Petrology, 27, 745-750.
 Muller, A., et al. (1998) New Aspects of Electron Transfer Revealed by the Crystal Structure of a Truncated Bovine Adrenodoxin, Adx(4-108). Structure, 6, 269-280.
 Pearce, J.A. (1983) Role of the Sub-Continental Lithosphere in Magma Genesis at Active Continental Margins. In: Hawkesworth, C.J. and Norry, M.J., Eds., Continental Basalts and Mantle Xenoliths, Shiva Cheshire, UK, 230-249.
 Wu, B.-M., Jahn, F.-Y., Wilde, S.A., Lo, C.-H., Yui, T.-F., Lin, Q., Ge, W.-C. and Sun, D.-Y. (2003) Highly Fractionated I-Type Granites in NE China (II): Isotopic Geochemistry and Implications for Crustal Growth in the Phanerozoic. Lithos, 67, 191-204.
 Wang, Y., Bryan, C., Xu, H. and Gao, H. (2003) Nanogeochemistry: Geochemical Reactions And Mass Transfers in Nanopores. Geology, 31, 387-390.
 Moghazi, A.M. (2003) Geochemistry and Petrogenesis of a High-K Calc-Alkaline Dokhan Volcanic Suite, South Safaga Area, Egypt: The Role of Late Neoproterozoic Crustal Extension. Precambrian Research, 125, 161-178.
 Yogodzinski, G.M., Lees, J.M., Churikova, T.G., Dorendorf, F., Woerner, G. and Volynets, O.N. (2001) Geochemical Evidence for the Melting of Subducting Oceanic Lithosphere at Plate Edges. Nature, 409, 500-504.