Petrogenesis and Metallogenesis of Malek Siah Kuh Adakite-Like Rocks and Associated Hydrothermal Mineralization ( Sistan Area , Iran )

The granitoids-diroretic complex and associated volcanic suits of Malek Siah Kuh (eastern Iran) were examined by RS, XRD, XRF, AAS and ICP-MS methods. The Multispectral data were used for differentiating of lithological units, main alterations, structural trends, and hydrothermally altered minerals. Based on the results of remote sensing, field studies and petrographic data, the complex classified into three main groups: 1) mafic unit (gabbro-basalt), 2) intermediate rocks (diorite-andesite suite) and 3) felsic unit (granitoide rocks). The complex is characterized by high Al2O3 and Sr contents as well as high Sr/Y and La/Yb ratios, while low in Y and Yb contents. The geochemical characteristics of these rocks indicate that they are similar to adakite like rocks. The complex is high calc-alkaline to shoshonitic nature. The results indicated derivation from similar silicate melt source followed by differentiation and assimilation of crustal rocks. Based on set of petrological and geochemical relations, 3 distinct stages of mineralization process can be recognized. 1) It is differentiation during magma evolution as well as metal enrichment and subsolidus leaching by hot fluids. 2) It occurred by ascending of the gas-rich magma and release of S and As. 3) It was accompanied by hydrothermal activities leading to the generation of sulfidation and associated polymetallic (Cu, Au, Fe) mineralization.


Introduction
Granitoid rocks are characterized by a large compositional diversity arising from different source compositions, variable melting conditions, chemical complexity, magma mixing, fractional crystallization and crustal contamination [1] [2] [3].
The study area shows geological and structural complexity due to its location between major structural zones of the Lut block, Afghan block, central Iranian blocks and several magmatic episodes.As a result an important mineralization is originated from this complexity.Their complex lithological units and strong re-mobilization cause ambiguity on the primary setting of the ore deposition.
The aim of this study was to investigate the mechanism of magma differentiation processes, parental magma compositions, changes and processes involved in the formation of the complex and related ore deposits.

Geological Setting
The study area is located between 60˚50' -60˚57' East longitudes and 29˚53' -29˚55' North latitudes (40 km north of Zahedan) in Sistan and Baluchestan Province (Iran).The Sistan region is located in southeastern Iran close to the Iranian borders with Pakistan and Afghanistan.Late Cretaceous to Tertiary rocks of the Sistan suture zone separate two blocks of continental character (Figure 1): the Lut block to the west and the Afghan block to the east [19].The Sistan suture zone (SSZ) represent a north-south trending segment of an ocean (the Sistan Ocean) in Early Cretaceous [20].Since Early Cretaceous the Sistan suture zone (SSZ) has undergone a complex history marked by rifting, subduction, ophiolite emplacement, continent-trench collision, uplift, at least three phases of Cenozoic deformation and several magmatic episodes have been recorded throughout its dynamic history [21].
The study area is mainly comprised of the plutonic and volcanic rocks (from Cretaceous to Eocene volcanic group).The granitoid plutons are elongated in a NW-SE direction and have been called the Zahedan-Saravan granitoids [21].
According to the nature of the components and degree of mixing, the Sistan Suture Zone rocks have been grouped into 3 basic lithotectonic units: 1) ophiolites in various degrees of disarray (with Senonian to Maastrichtian ages); 2) phyllite and phyllonite of late Cretaceous to Eocene age; and 3) terrigenous marine sedimentary rocks [20].The phyllite and phyllonite are weakly metamorphosed turbidite successions constituting almost half of the Sistan Suture Zone [20].
K-Ar dating of biotite samples from different areas in the intrusive belt indicates solidification ages in the range of 31 -33 Ma for the main phase intrusions [21].From the middle Cretaceous to recent (Figure 2), the area of the Sistan suture zone has undergone a number of important tectonic events and continued convergence between the Lut and Afghan blocks resulted in widespread folding and conjugate strike-slip faulting [20].Several mineral deposit occurrences with variable Au and Cu anomalous are hosted in the vicinity of granitoid, mafic and ultramafic rocks and iron oxide veins.

Analytical Methods
Remote sensing techniques have been used for identification of lithological rock units and altered zones and were confirmed on the field before sampling.Multispectral and hyperspectral remote sensing has been used extensively around the globe for lithological and mineral mapping [22] [23] [24] [25].With the advent of multispectral sensors, lithological discrimination and mineral potential mapping were possible from space/airborne platforms [26] [27].Hyperspectral Open Journal of Geology sensors provide a unique combination of both spatially and spectrally contiguous images that allow precise identification of minerals [28].
A total of 35 samples (from volcanic and plutonic rocks) were collected from outcrops and 15 samples were selected for determination of major and selected trace and rare element analysis.Major element contents were determined by XRF whereas trace and rare elements were determined by inductively coupled plasma mass spectrometer (ICP-MS) at Geological Survey of Iran.
Detailed mineralogical studies and modal analysis were performed by X-ray Diffraction (XRD) but textural analyses were carried out by Petrographic microscope and scanning electron microscope (SEM) equipped with an EDS system, in the electron microscopy laboratory of Ferdowsi university of Mashhad.X-ray diffraction (XRD) analysis of the crystallographic phases was performed on over 15 samples.Quantitative XRD analysis was accomplished in order to accurately determine the mineralogical composition of fine-grain volcanic rocks and altered products.All samples were crushed to a fine powder and placed on single silicon crystal sample holder.XRD patterns were collected from the powdered samples using a CuKα source.All samples were analyzed from 5˚ to 75˚ 2θ.Operating conditions were 40 kV X-ray beam and 35 mA X-ray beam current.
Mineral phases present were identified by computer-aided searches using Modern SaxAcq acquisition software.

Remote Sensing
The remote sensing data plus field and petrological studies have provided possible separation of lithological units with high accuracy in the study area.The Multispectral data was used for differentiating of lithological units (Figure 3), main alterations (Figure 4), structural trends, hydrothermally altered minerals (Figure 5), surface mineral identification and mapping (Figure 6).The main litholigcal units and structural trends were identified as a linear body trending NW-SE direction parallel to the major shear zone of the area.The results of remote sensing (RS) and petrographic data provided a means for classification of lithological units into 3 units.The granitoid stock consists of monzogranite, granodiorite unites in the south but quartz-diorites and quartz-monzogranites in the north of the area.Analysis of remote sensing data and field studies suggested that mafic to intermediate and acidic intrusive-extrusive sequence is present in the study area (Figure 6).Moreover, hyperspectral data was used to distinct the location of hydrothermally alterated minerals in Malek Siah Kuh area.Hematite and arsenopyrite with minor quantity of limonite and goethite are found in the study area.

Petrography
The igneous complex at Malek Siah Kuh area is comprised of a wide range of mafic and intermediate to acidic intrusive-extrusive sequences.The complex can be subdivided into three groups: 1) mafic unit (gabbro-basalt), 2) Intermediate rocks (diorite-andesite suite and much less trachyandesites) and 3) felsic unit (granitoide rocks).The predominant units in plutonic rocks are chains of diorite and granodiorite suits.The diorite rocks in the Malek Siah Kuh batholith are exposed in the north-eastern part of the study area.These rocks are medium to coarse grained, characterized by granular, ophitic to subophitic textures.Although all types of diorites in the study area are mineralogically similar but the dioritic unites comprise varying lithologies of altered gabrodiorite, meladiorite, mesodiorite and leucodiorite.The principal constituent minerals in these rocks are plagioclase (andesine) (44% to 56%), orthoclase (8% to 16%) and β-quartz (5% to 15%), but augite (5% to 15%), hornblende (15% to 20%), epidote, calcite, zoisite and clinozoisite, occur as minor constituents.The XRD analysis of dioritic altered rocks revealed the assemblage of calcite, zoesite, clinozoesite and epidote.In meladioritic rocks, saussuritization products (calcite, clinozoesite and epidote) occur as altered minerals.The most abundant accessory minerals are  The volcanic rocks consist mainly of andesite, trachyandesite and basalt (Figure 7).They are dark-colored to brown with phenocrysts of plagioclase (labradorite to andesine), mafic minerals (hornblende, Pyroxene) and alkali feldspar with porphyritc and gleumeroporphyritc textures.The more Ca-rich members of the plagioclase series (phenocrysts) are completly zoned with inner zones altered to fine cryptocrystalline mixture (sausuritization).The process of sausuritization had been produced calcite, zoesite, prehnite, clinozoesite and epidote.In volcanic rocks (andesites and trachyandesites), the fine groundmass consists of plagioclase (andesine), hornblende, pyroxene (augite) and alkali feldspar.Moreover, titanomagnetite and opaque minerals are present as accessory minerals (Figure 8).
Granitoid unit is the most abundant plutonic rock in the area and consists mainly of granodiorite but includes other rocks such as monzogranite and quartz monzonite.Granodiorite rocks are often found related with typical hydrothermal minerals and they vary from medium to coarse grained.In this group, the dominant mineral assemblage is plagioclase (oligoclase) (40% to

Geochemestry
Geochemical data from the igneous complex have been used to characterize the petrogenesis and tectonic setting of the plutonic and associated volcanic rocks of the study area.Major oxides, trace and rare earth element (REE) abundances of 15 samples from the complex were analyzed by XRF, ICP-MS and AAS (Table 1).The result data indicate that unaltered trachyandesite is silica oversaturated with SiO 2 ranging from 55.46% to 60.45% and alkali contents (K 2 O + N 2 O) ranging from 6% to 10% but the intermediate rocks of the area exhibit a narrow range of SiO 2 content (55% to 60 wt.%).
Moreover, the amount of SiO 2 in granodiorites ranging from 63.31% to 64.1% and total alkali contents from 5.56% to 6.46%.As a result, these granodiorites belong to high-K calcalkaline series.The granodiorite is characterized by high LREE contents and low HREE contents.All of the geological and petrological evidences suggest that quartz-diorite, quartz-monzodiorite and granodiorite rocks in Malek Siah Kuh apparently form as a result of fractional crystallization of source magma.
In the TAS diagram [29], the most rocks plot in the fields of sub-alkaline basalt, andesite and trachyandesite (Figure 9).However, variations of silica versus the FeOt/MgO and standard AFM diagram (Figure 10) as well as silica versus  K 2 O (Figure 11) showed that most of the samples are calc-alkaline in nature.Thus, due to chemical and petrological relationships, these rocks are considered part of the calc-alkaline province which extends from Lut block into Sistan area.Therefore, the geochemical characteristic of igneous complex in Malek Siah Kuh was related strong subduction components.In the complex, TiO 2 increases with decreasing SiO 2 , that consistent with magma evolution during differentiation from a mafic magma [30].Despite the clear link between the compositions of plutonic complex that can be related to metaluminous or peraluminous but some volcanic rocks in the Malek Siah Kuh area are shoshonitic.They are enriched in LREE and depleted in Ti, Nb and Ba.Plot of major elements versus silica (Figure 12) indicates that the rocks of the complex show smooth variations with gradual decreases in CaO, FeO (total) , and TiO 2 with increasing silica.The positive correlation of Na 2 O and incompatible elements (Ni, Zr and La) with SiO 2 (Figure 12 and Figure 13) and the negative  correlations of TiO 2 , MgO, CaO, FeO (total) and P 2 O 5 with SiO 2 (Figure 12 and Figure 13) show the impact of differentiation of magma in the magmatic evolution of igneous complex in the Malek Siah Kuh area.Beside these plots, the result of XRF and ICP-MS in Harker diagrams showed the effect of differentiation process as well as contamination phenomenon.Diorite rocks of the complex are characterized by high SiO 2 and Al 2 O 3 contents but intermediate alkali and MgO contents.Geochemically, diorite-andesite unit is predominantly calc-alkaline and metaluminous in characteristics.These compositional characteristics suggest the parental magma of andesitic composition and this kind of high-aluminum magma is typically produced in subduction related settings particularly at continental margins [31].Although the igneous complex in Malek Siah Kuh, crystallised from basaltic magma, however, granodiorites and trachyandesites are peraluminous and may have formed by assimilation of felsic continental crust.
Although, field relations of intermediate plutonic and acidic rocks do not indicate a direct petrogenetic link, but lithogeochemistry investigations showed some similarities in trace and rare earth element abundances.These evidences may indicate derivation from similar silicate melt source during differentiation from a mafic magma [30] and assimilation of crustal rocks.Immobile, incompatible elements (such as Ti and Zr) exhibit highly significant positive correlation coefficients (r) in cogenetic lavas, even when altered or metamorphosed [32].In another case, the plot in the magnesian granitoid field indicates hydrous, oxidizing condition of the magma.The origin of magnesian granitoids has been interpreted as related to subduction processes [33].
Petrological studies and geochemical data (including REE and trace element analyses) has highlighted the fact that magma mixing was not responsible for the Open Journal of Geology formation of intermediate rocks.In contrast, the trace element analyses and enrichment of LILE of the igneous complex relative to HFSE is suggesting calc-alkaline magmatism and involvement of strong subduction components.
The igneous complex shows a wide range of Zr (70.2).Broadly, the characteristics of these rocks are geochemically similar to adakite.In addition to the petrological aspects, the chemical trend of major and trace elements from intermediate plutonic to acidic rocks showed arc related active continental margins (Figure 15 and Figure 16).The rocks in the study area are complex and in order to take account of this complexity, the diagram of Sr/Y versus Y (Figure 17) was used to distinguish adakite type rocks from normal arc igneous rocks.This is a practical illustration that plots of MgO wt.% vs. SiO 2 , Sr vs. (CaO + Na 2 O) and Sr/Y vs. Y show the chemical differences between low silica adakite and high silica adakite types (Figure 18).Thus, plotting La/Sm ratio vs Sm/Yb ratio diagram (Figure 19) would allow for deducing the possible mantle source in the study area.It is clear that the presence of igneous activities (in the form of diorite-granodioritic complex), hosting the mineralization indicate a role of magma in the generation of ore forming magmatic hydrothermal systems.Adakites in subduction zones are now recognized to be particularly favorable for Cu-Au mineralization [8].Porphyry Cu-Au deposits are generally derived from sulfur-rich, highly oxidized magmatic systems, with oxygen fugacities (fO 2 ) between the nickel-nickel oxide (NNO) and magnetite-hematite oxygen buffers [12].[12] suggested that only slab-derived adakitic magmas or supercritical fluids would have a high oxidation potential to generate epithermal and porphyry Cu-Au deposits.
Cu and Au mineralization in Malek Siah Kuh is mainly associated with diorite-granodiorite complex and by hydrothermal activities.In the study area, wall

Conclusions
The granitoid of Malek Siah Kuh is metaluminous to per aluminous and the complex is high calc-alkaline to shoshonitic nature.

Figure 1 .
Figure 1.(a) Location of the Sistan suture zone, the Lut and Afghan blocks, and surrounding physiographic and geologic features [21]; (b) Major geologic subdivisions of the Sistan suture zone [20].

Figure 3 .
Figure 3. Color view of Malek Siah Kuh area and its surrounding area.

Figure 5 .
Figure 5. (a) Malek Siah Kuh area and its surrounding area; (b) Yellow colours show hydrothermal alteration zones; (c) Distribution of iron minerals (hematite and goethite) on the Malek Siah Kuh area.

Figure 6 .
Figure 6.Simplified geologic map of the study area showing the spatial distribution of igneous complex in the Malek Siah Kuh area.

Figure 12 .
Figure 12.Selected Harker diagrams showing the major element variations of igneous rocks in the Malek Siah Kuh area.

Figure 13 .
Figure 13.Selected Harker diagrams showing the trace element variations of igneous rocks in the Malek Siah Kuh area.

Figure 14 .
Figure 14.(a) Spider trace element; (b) REE elements normalized to MORB in igneous rocks in the Malek Siah Kuh area.

Figure 17 .Figure 18 .
Figure 17.Plot of Sr/Y vs. Y.The diagram typically used to distinguish adakite from normal arc igneous rocks [4].

Figure 19 .
Figure 19.La/Sm vs Sm/Yb diagram to determine the possible mantle sources of the studied rocks.Most of the samples are in the range of amphibolite derived rocks [34].

Figure 20 .
Figure 20.Photographs of typical alterations from Cu-Au mineralization in the Malek Siah Kuh area.
Petrology and geochemical data of mineralized granitoid intrusions of the study area showed a geochemical relationship with adakitic rocks with high Al 2 O 3 (15.5% wt.%), MgO (4 wt.%), Sr (800 ppm) contents as well as high Sr/Y ratios, light rare earth element (LREE) enriched but low in Y (12 ppm), Yb (2 ppm) contents and heavy REE (HREE) depletion.Based on set of petrological and geochemical-mineral deposit relations, three distinct stages of mineralization process can be recognized in the study area.The first stage is related to differentiation of primary magma, metal enrichment and subsolidus leaching by hot fluids.The second phase appeared by ascending of the gas-rich magma (reduction in pressure) and the release of S and As.The final stage was accompanied by hydrothermal activities leading to the generation of sulfidation (formation of Au and Cu sulfide) and associated polymetallic (Cu, Au, Fe) mineralization.

Table 2 .
[35]arison of the geochemical features of the studied rocks with geochemical features of adakite[35].