Abstract

The Qatruyeh iron deposits, located on the eastern border of the NW-SE trending Sanandaj-Sirjan metamorphic zone, southwest of Iran, are hosted by a late Proterozoic to early Paleozoic sequence dominated by metamorphosed carbonate rocks. The magnetite ores occurred as layered to massive bodies, with lesser amounts of disseminated magnetite and hematite-bearing veins. Textural evidences, along with geochemical analyses of the high field strengths (HFSEs), large ion lithophiles (LILEs), and rare earth elements (REEs), indicate that the main mineralization stage occurred as low-grade layered magnetite ores due to high-temperature hydrothermal fluids accompanied by Na-Ca alteration. Most of the main ore-stage minerals precipitated from an aqueous-carbonic fluid (3.5–15 wt.% NaCl equiv.) at temperatures ranging between 300° and 410°C during fluid mixing process, CO2 effervescence, cooling, and increasing of pH. Low-temperature hydrothermal activity subsequently produced hematite ores associated with propylitic alteration. The metacarbonate host rocks are LILE-depleted and HFSE-enriched due to metasomatic alteration.

1. Introduction

The last decade has seen major progress in our understanding of the origin of iron ore deposits worldwide. The majority of interpretations focused on the igneous iron oxide deposits either having formed by magmatic liquid immiscibility [17] or by hydrothermal alteration and replacement [814]. The hydrothermal iron deposits generally are found at many locations around the Pacific basin, Central America, Australia, and Japan [15]. They are commercially far less important as global source of iron than banded iron formations and igneous iron deposits, except for many countries without these types of iron deposits. Genetic models for hydrothermal Fe-oxide deposits associated with hybrid magmatic to nonmagmatic fluids, which locally mix at the trap site [16], suggest that metamorphic processes or a primary felsic to intermediate intrusion could be regarded as source of the majority of hydrothermal fluids and metals [17, 18]. Therefore, metasomatic iron-rich fluid regimes are associated with magmatism or metamorphism, or both. Depending on the physicochemical conditions the protore iron mineralization could consist of magnetite or hematite, or a mix of the two. They occur in different tectonic environments, such as intracontinental terranes associated with anorogenic magmatism, continental arc terranes, and metamorphic belts [1820].

The most favorable area for hydrothermal ore deposits in Iran is the Sanandaj-Sirjan zone (SSZ), which has a NW-SE trend and is characterized by magmatism-metamorphism and an obliquely thrusted wedge, with asymmetrical structures in the HP-LT metamorphic rocks [21]. This metamorphic belt contains abundance of iron ore deposits, including the Gol-Gohar (Kerman), Totak, Dehbid (Fars), Galali, Baba-Ali (Hamedan), and Shams-Abad (Arak). However, distribution of these deposits from NW toward SE Sanandaj-Sirjan is very important but the origin of these deposits has been a subject of heated debate. Eshraghi et al. [21] were the first to present the details of the geology and genesis of the Gol-Gohar and the Qatruyeh deposits and proposed metamorphism as main factor controlling the mineralization. However, based on textural and geochemical evidence, [2224] suggested iron oxide-copper-gold (IOCG) origin for Sanandaj-Sirjan iron deposits. According to [25], the iron deposits from Hamadan to Gol-Gohar situated in the Sanandaj-Sirjan zone about 1200 km apart and show striking mineralogical and textural similarities. They believed that the ore bodies are of magmatic origin and formed at late stages of hydrothermal conditions.

In spite of the existence of the important iron ore deposits in the Qatruyeh area, there are a few reports on the ores with no exploitation activities. This paper presents detailed field observations, mineralogical alteration assemblages, textural relationships, fluid inclusions and geochemical data on the ore bodies, and the associated dolomitic-marble host and plutonic rocks at the Qatruyeh iron ores, southwest of Iran, in order to reveal a genetic model for the deposits in the area.

2. Methods

Major and trace elements analyses were carried out on representative samples using inductively coupled plasma optical emission spectroscopy (ICP-OES), with 0.1 ppm detection limit at Geological Survey of Iran (22 ore samples), and using standard wet chemical analysis by flame atomic absorption spectrophotometry (FAAS), with 0.5 ppm detection limit at the Department of Earth Sciences, at Shiraz University (5 fresh marbles, 11 altered marbles, and 12 diorites). Rare earth elements (REE) were analyzed by inductively coupled plasma mass spectrometry (ICP-MS) with 0.01 ppm detection limit at ALS Chemex Analytical Services, Canada (4 diorites and 22 ore samples). Sample locations are shown in Table 1 and they are described with regard to type of analyses.

3. Geological Setting and Lithostratigraphy

The Qatruyeh area is located about 50 km NE of Neyriz and lies within the Zagros orogenic belt in SW of Iran. The Zagros orogenic belt of Iran consists of three parallel tectonic domains: (1) the Urumieh-Dokhtar magmatic belt (UDMB), (2) the Sanandaj-Sirjan zone (SSZ), and (3) the Zagros folded-thrust belt (ZFTB) [35]. The iron deposits are located in the eastern margin of the NW-SE trending Sanandaj-Sirjan zone (Figure 1).

The SSZ consists of Paleozoic to Triassic high-grade metamorphic rocks overlaid by Mesozoic metasedimentary strata. The Mesozoic sediments were deposited along a passive continental margin [3639] and subsequently large-scale Mesozoic plutons ranging from gabbro to granite, emplaced in about  Ma by the K-Ar methods [40]. The ore bodies show a strong deformation in parallel to the regional west-northwest structural trend. Lithostratigraphic relationships in the Qatruyeh area are summarized in Figure 2. Some researchers [40, 41] recognized five rock groups in Neyriz area, from base to top and from low- to high-grade metamorphic rocks, which are unconformably overlaid by Upper Triassic turbiditc rocks and basalt. These groups include the following.

Group A. Homogeneous orthogneiss with some amphibolite lenses and small bodies of anatectic granite. The orthogneiss has been dated by U-Pb method on zircons, at 520 Ma [41].

Group B. Alternation of amphibolite, migmatite, minor micaschist, marble, and metamorphosed mafic and ultramafic rocks. The age of the metamorphic rocks is 310 Ma based on 40K-40Ar dating method on extracted white mica, biotite, plagioclase, and amphibole.

Group C. Micaschist, calc-schist, marble, minor amphibolite, and mafic to ultramafic rocks. Along the Sanandaj-Sirjan zone, comparable rocks yielded deformed Devonian corals [21, 42, 43]. The Chah Pazhan metagabbro contains medium-grained clinopyroxene (e.g., augite, hedenbergite, and pigeonite), plagioclase feldspar, and orthopyroxenes (e.g., enstatite). The clinopyroxenes have been amphibolized and uralitized (actinolitization-tremolitization).

Group D. Schist, slate, calc-schist, metabasite, quartzite, and intercalations of marble. The metacarbonate intercalations are more frequent towards the upper part of the group. Palynomorphs (Galyptosporites, Laevigatosporites, Polyedrxium sp., Ancyspora sp., Geminospora sp., and Relusotrilletes sp.) indicate a late Devonian-early Carboniferous age for this group [21].

Group E. Limestone, dolomite, basalt, and synsedimentary monogenetic breccias. Above the early Cimmerian unconformity, the Upper Triassic and Jurassic sequence were deposited in turbiditic basins, geometry of which is still poorly known [40]. These rocks show a weak but obvious cleavage, indicative of a deformation that affects the upper Triassic and Jurassic conglomerates and sandstones, but not the overlying Cretaceous series in adjacent areas.

The Deh Vazir Jurassic deformed conglomerate is located in the north part of the area and is sandwiched between two thrust sheets which are part of the Main Zagros Thrust zone. Most of the pebbles in this conglomerate consist of micaschist and phyllite [44], but at some localities, the presence of massive magnetite pebbles indicates the ore formation is prior to main metamorphism (Laramide phase) (Figure 3(a)). The conglomerate was occurred extensively in the area which covers metamorphosed limestones, turbidite shales, and sandstones.

4. Deformation and Metamorphism

The SSZ underwent polyphase deformation events during the Mesozoic era. The latest event recorded in the rocks reflects collision of the Afro-Arabian continent and Iranian microcontinents and subsequent southward propagation of the fold-thrust belt [39, 45, 46]. Geological and textural information suggest that there are at least three main deformation phases D1, D2, and D3 [47].

The regional lithologies at the Qatruyeh record a clockwise metamorphic path followed by thrust deformation events. The peak pressure and temperature in the upper Triassic and Jurassic sequences were estimated at  kbar and °C, respectively. Using zircon sensitive high-resolution ion microprobe (SHRIMP) U-Pb and monazite chemical Th-U-total Pb isochron methods (CHIME) presumably regional metamorphic event occurred between  Ma and  Ma, respectively. These rocks experienced clockwise P-T path and thrusted into the imbricate slices during exhumation [48]. The amphibolite facies metamorphism in early Cimmerian orogeny could be correlated with the regional D1 deformation. Exhumation of the metamorphic complexes and the resultant cooling and decompression of the metamorphic rocks followed by development of postearly Cimmerian sedimentary basins occurred in a transpressive regime due to subduction of the Neotethys. As a result, calc-alkaline mafic to felsic intrusive rocks were emplaced in the northeastern part of the region. A high-pressure-low-temperature metamorphic belt was formed during the first stage of the deformation. The result of this deformation is a composite foliation between massive magnetite and dolomitic marble that occurred along a wide ductile shear zone [40].

The D2 deformation is associated with overthrusting of the old units (e.g., the Chah Pazhan ultramafic-mafic rocks) onto the Qatruyeh metacarbonate unit and the development of a mylonitic foliation in the various rocks units. Veinlets of iron oxides display pinch and swell features along the foliation. These relationships indicate that percolating fluids moved through permeable zones of foliation. The resulting boudinage show that ductile deformation went on after formation of iron oxides (hematite) and dolomitic-marble host rock (Figure 3(b)).

Whereas D1 and D2 developed under amphibolite and greenschist facies metamorphic conditions, D3 represents an event leading to the development of kink bands, localized shear zones, and slickensides along the thrust sheet of the Qatruyeh area [47]. Abundant kink bands of several centimeter thicknesses overprinted all other deformation structures such as both magnetite and hematite ores.

5. Ore Bodies’ Geology

There are twelve ore bodies in the Qatruyeh area. Kuh-Sorkh and Kuh-Ghermez were chosen for detailed field investigation, petrography, and chemical analyses. Both ore bodies occur within the metasedimentary rocks of the early Paleozoic which consists predominantly of dolomitic limestone, greenschist, and quartzite (group B). The iron oxide ores can be divided into massive and layered magnetite ores with minor hematite. The ore bodies were not exploited due to their low tonnages in the area.

The massive magnetite ore bodies are often in contact with micaschist and dolomitic marble. They occurred often as very small sized masses with maximum 5 m2 area which replaced metacarbonates. The ore bodies have variable orientation and shape and commonly have been fractured and cut by quartz-bearing hematite veins. The layered magnetite ores, however, have few centimeters to several meters width (>10 m) and are found throughout the area interbedded with the stratified metacarbonate rocks (Figure 3(c)). These host rocks were strongly affected by sodic-calcic alteration. These sodic-calcic alteration zones are spatially associated with layered magnetite bodies that appear to have replaced marbles, dolomitic marbles, and chlorite-schist rocks. Pervasive sodic-calcic alteration grades into zones of layered, actinolite crystals intergrown with magnetite. Most intense sodic-calcic alteration is dominated by pervasive rock replacement by siderite + paragonite with associated layered magnetite ores. The ores are most often horizontal to subhorizontal and elongated WNW-ESE trending, concordant with the regional structures.

Minor disseminated subhedral to euhedral magnetites are common in the area, dominantly appears in chlorite-schist (Figure 3(d)). Minor hematite occurrences are found in the upper parts of the brecciated and mylonitized metacarbonates (Figure 3(e)) that display propylitic alteration near the fault contact.

The metacarbonates and chlorite-schists adjacent to intrusive rocks (intrusions) are crossed by later veins of hematite + quartz ± pyrite ± chalcopyrite (Figure 3(f)). Such veins are generally poorly mineralized and appear to represent a late, high-level zone of hydrothermal alteration.

6. Petrography

The metamorphic rocks are extensively formed in the Qatruyeh area which is composed of metacarbonate rocks (marble and dolomitic marble), chlorite-schist, micaschist, and metagabbroic schist with small amounts of amphibolite. A heterogeneous intrusion with about 1 km in diameter is located in the south of the area showing a chilled margin. It has been emplaced within the calcareous sandstone and shale whose Dictyoconous sp. microfossils indicate late Jurassic to early Cretaceous age (Figure 4(a)). The plutonic rocks range from dioritic to monzodioritic in composition.

The dolomitic-marble rocks appear white (fresh) to yellowish-brown (altered) in color due to the absence or presence of siderite, respectively (Figure 4(b)). The rocks are considered to be extensively replaced by magnetite and hematite. The rhythmic banding in dark altered and light colored fresh rocks is easily distinguished in the field.

Microscopic observations indicate several replacement textures due to alteration. These include widening of a fracture filling where a fracture crosses a chemically reactive rock (Figure 5(a)), irregular or vermicular intergrowths of hematite and calcite-dolomite (Figure 5(b)), remnant islands of nonreplaced host rock, rims of one mineral penetrating another along its crystallographic directions (Figure 5(c)), nonmatching walls or borders of a fracture, and cusp and caries textures (Figure 5(d)) [49].

7. Alteration and Mineral Assemblages

Metasomatism accompanied by intrusions, and/or metamorphism, may be responsible for iron oxide mineralization in the Qatruyeh area. Microscopic observations and X-ray diffraction patterns on ore samples and metacarbonates have led to the recognition of two main alteration stages in the layered magnetite- and hematite-bearing ores (Table 2).

The first stage of alteration includes Na-Ca metasomatism (actinolite + titanite + epidote + calcite + albite ± scapolite ± quartz), which is characteristic of diorite-monzodiorite plutonic rocks (Figure 6(a)). Epidote and actinolite replaced plagioclase; magnetite inclusions were found in some actinolites (Figure 6(b)). Magnetite, actinolite, tremolite, tourmaline, titanite, and quartz are the major minerals; pyrite and chalcopyrite occur as minor phases in the accompanying layered ore bodies (Figure 6(c)). Layered magnetite ores replaced dolomitic marble and the calc-silicate minerals that are typical of skarns are virtually absent from the ore zones [50]. Paragonite and siderite are exclusively identified within the halos of 1–3 m in diameter close to ore bodies that were only recorded at XRD data. Paragonite was formed during the early metasomatic stage. The same data records were observed in Lena Gold District, Siberia [51]. Siderite is produced by replacement of calcite or dolomite by an iron-rich solution [52] and paragonite is formed by solutions during regional metamorphism [51].

The second stage of alteration is shown by hematitization in brecciated metacarbonate (Figure 6(d)). Hydrothermal activity has produced a specific alteration around hematite ores, with a more regionally developed propylitic zone. This alteration style is very common in hematite ores and indicates the activity of a volatile-rich hydrothermal solution [6]. The alteration that is associated with oxidizing fluids often results in the formation of minerals with a high Fe3+/Fe2+ ratio and, in particular, hematite associated with sericite + quartz + epidote + chlorite ± magnetite [53]. Hematite + quartz ± pyrite ± chalcopyrite veins accompanied the late stage of each main ore mineralization in the Qatruyeh area.

8. Results

8.1. Geochemistry

In order to study the major, trace, and REE compositions of the Chah Anjir plutonic rocks, eleven samples were collected from surficial exposures. The emphasis during sampling was on fresh rocks, so as to avoid metasomatic changes during hydrothermal alteration.

Results of the whole rocks chemical analyses of major, trace, and REE in the Chah Anjir migmatite are summarized in Table 3. According to - discrimination diagram [26] and the classification diagram for plutonic rocks [28], the Chah Anjir intrusive rocks can be classified as diorite, monzodiorite (Figure 7(a)), and diorite to syenodiorite, respectively (Figure 7(b)), with subalkaline to alkaline affinity (Figure 8(a)). As indicated in Figure 8(b), the plutonic rocks plot within the calc-alkaline to tholeiitic fields. In order to determine the tectonic setting for the plutonic rocks and its relation to ore mineralization, Y − Nb and (Y + Nb)-Rb diagrams [30] were used. On the Y − Nb diagram, the samples plot in the volcanic arc to syn-collisional granite fields (Figure 8(c)). On the Rb-(Y + Nb) diagram the samples fall in the field of volcanic arc granitoids (Figure 8(d)). The tectonic setting is confirmed to geological characteristics of the Qatruyeh area.

In the A-B diagram of [31], the samples show a negative slope in IV sectors of metaaluminous area that is such a trend characteristic of cafemic (CAFEM) rocks due to presence of hornblende, ortho- and clinopyroxene, biotite, epidote, and titanite in origin (Figure 9(a)). According to the A/CNK <1.1 ratio [54] and Zn versus SiO2 diagrams (Figure 9(b)) [32] and SiO2-Fe2O3/FeO diagram [33] the Qatruyeh plutonic rocks are considered I-type and magnetite series granitoids (Figure 9(c)). These results associated with those of the intrusions fall in the hydrothermal iron deposits field.

REEs from the samples were normalized to chondrites using the values of Wakita et al. [55]. Enrichment coefficients of La (), La (), Gd (), and [Eu]/[Eu*] as Eu anomaly were examined. The results indicate that the Chah Anjir intrusive rocks were enriched in light rare earth elements (LREEs) compared to heavy rare earth elements (HREEs) (Figure 9(d)). Moderate to negative Eu anomalies, high total REE contents ( with average of 168), and high LREE/HREE ratios () are also observed. These evidences suggest continental or continental-margin geological settings for the intrusive rocks in the Qatruyeh area. The negative Eu anomalies can be explained by the abundance of plagioclase feldspar in the melt. As proposed by Henderson [56], small amounts of garnet, amphibole, or pyroxene can reduce the HREE relative to LREE in the melt.

As mentioned, iron ores could be divided into high-grade magnetite ores (consisting of massive magnetite with minor silicate minerals) and low-grade magnetite ores (layered magnetite with silicate minerals) and hematite ores. The major, trace, and rare earth elements of the Qatruyeh iron oxide ores are presented in Table 4. The high-grade magnetite ores typically have high contents of TiO2 (0.09 to 1.08 wt.%), whereas the low-grade ores indicate the high amounts Na2O + K2O (avg. 1.08 wt.%) and CaO (avg. 3.8 wt.%).

Enrichment of sulfur in the layered magnetite is reflected by the presence of sulfides such as pyrite and chalcopyrite. There are strong positive correlations among Fe2O3, TiO2 and HFSEs (Hf, Sc, Th, U, Nb, and Zr), LREEs, and transition metals (Cr, Co, and Ni) in layered magnetite. In contrast, the strong negative correlations between the two oxides with HREEs and LILEs (Ba, K, and Rb) indicate that layered magnetite acts as a scavenger for HFSEs, LREEs, and intermediate elements (Figure 10). It was inferred that the relationship was due to the similarity of the behavior of those elements in spinel group minerals.

The positive correlation between LREEs and Al2O3, CaO, Na2O, and MgO indicate that the LREEs are collected in the silicate phases in the ore bodies.

High values of chalcophile elements (Hg, Bi, Cu, and Ag) were observed in layered magnetite as a below-ore halo and deep-seated iron mineralization.

The REE distribution patterns in high-grade massive magnetite and low-grade layered magnetite ores are very different (Figure 9(d)). These patterns can show differential origins of both ore bodies. An average ratio of 6.3 in the layered ores indicates a high degree of fractionation. A weak negative Eu anomaly and LREE enrichment are observed in the layered magnetite. In contrast, high-grade ores are depleted in LREE, enriched in HREE, and display a strong negative Eu anomaly. The depleted LREE pattern is interpreted as high-temperature leaching of LREE or the occurrence of a lower temperature metamorphism [57]. It is suggested that Ce depletion in high-grade magnetite may be the result of a late stage metasomatic alteration in the Qatruyeh area. This indicates that the two ore occurrences, layered and massive, are different in origin.

The bulk compositions of the altered and fresh metacarbonates are given in Table 5. Altered rocks are depleted in CaO and MgO relative to fresh metacarbonate samples whereas the Fe2O3, Al2O3, K2O, and Na2O contents increased in the rocks during metasomatic alteration. An increase in the K2O and Al2O3 contents is the result of muscovite abundance in the altered rocks. TiO2, Ba, Rb, Cr, and V can also be present within muscovite. The difference in the position of data points of altered and fresh rocks is explained in Figure 11.

High MnO and Fe2O3 values in altered dolomitic marbles indicate that Fe and Mn replace Mg in dolomite lattice. In order to evaluate the compositional change which accompanies hydrothermal alteration, the concentrations of elements of interest are compared to those of immobile elements such as Zirconium. Typically, Zr is used as immobile and incompatible element during hydrothermal alteration due to the very high ionic radius [58]. In the case of chemically incompatible elements, the regression lines REEs and HFSEs which are considered to be immobile with each other in the Qatruyeh altered dolomitic marbles (Figures 12(a)12(d)). The REEs are less well correlated with the above HFSEs and the regression lines shift from the origin, which suggests that the REEs were somewhat mobile during alteration. The interelement correlation among the REEs is significant (generally with coefficients higher than 0.7), which is consistent with their similar geochemical behavior. However, the interelement correlation among LREEs and HREEs treated separately is considerably stronger, which indicates fractionation between the two subgroups. LILEs (Ba, Sr, Rb, and Cs) are more mobile relative to Zr during hydrothermal alteration as shown by their poor correlations in Figures 12(e) and 12(f).

8.2. Fluid Inclusions Study

Microthermometric measurements were made using a Linkam THMSG600 combined heating and freezing stage with a temperature range of −196 to +600°C and a reproducibility of ±0.1°C, at the fluid research laboratory at the Geological Survey of Iran by optical microscopy and microthermometry, attached to an Olympus petrographic microscope with LinkSys software (version 1.83). Ten doubly polished quartz wafers were prepared using the procedure of Shepherd et al. [59]. The thickness of the wafers varied between 150 and 300 μm, depending upon the transparency of the quartz crystals. Sample selection was biased to quartz containing an abundance of Na-Ca minerals in low-grade magnetite type. The fluid inclusions have a variety of shapes (mainly rounded, elongate, or square and, less commonly, irregular) and typically are of low to medium size (5 to 50 μm). Quartz in the samples contains inclusions >10 μm which allows for more confident thermometric analysis. Based on petrographic and microthermometric studies and present phases at temperature room, inclusions are divided into the following types (Figures 13(a)–13(d)).

(1) Two-phase liquid + vapor (type A), (2) three-phase vapor + liquid () + CO2 liquid (L2) (type B), (3) monophase liquid (type C), and (4) monophase vapor (type D). Types A and B constitute the majority of fluid inclusions in the Qatruyeh quartz veins (proximal to the Chah Anjir intrusive unit) and account for 80% of the inclusion population and were chosen for microthermometric studies for two important reasons: (1) the inclusions are mostly intimately associated with layered magnetite mineralization, (2) these types contain inclusions >5 μm which allows for more confident thermometric analysis. In addition, primary fluid inclusions, with no evidence of necking down, were selected for microthermometric analyses.

Type A inclusions are dominated by more than 80 vol.% of H2O at room temperature. These inclusions occur in a variety of forms: in growth zones, small clusters, and as dispersed arrays. Measurements on trapped phases were not attempted. The first ice temperature of melting () often occurs around −24.5 to −19.5°C, corresponding to stable eutectic melting in the NaCl-MgCl2-CaCl2-H2O system (Figure 14(a)). Salinities determined by last ice temperature of melting () were 3.5 to 15 wt.% NaCl equivalent (size of inclusions between 10 and 30 μm) for inclusions trapped in whole quartz samples (Figure 14(b)). The homogenization temperature () values change between 280 and 385°C (Figure 14(c)) and Tm measurements range from −11.3 to −3.5°C (Figure 14(d)). Homogenization temperature existing for primary types (A) of fluid inclusion may indicate episodes of hydrothermal activities in mesothermal systems for Fe deposits in this region.

Type B inclusions, characterized by much more gas bubble within an aqueous liquid, are commonly found in all selected samples. The inclusions are dominated by more than 45 vol.% H2O-CO2 at room temperature. Carbonic liquid is a visible phase as immiscible solution within H2O liquid. But, vapor bubbles of the inclusions do not contain enough CO2 to form clathrate or a visible phase of liquid CO2 during cooling. This suggests that the CO2 content of the bubble is less than 0.85 mol.% [60]. of type B fluid inclusions (Figure 14(a)) is in the range of −39.5 to −37.6°C. The first ice melting temperatures indicate that the mineralizing fluid contained a complex mixture of metal chlorides, particularly NaCl-MgCl2-CaCl2-FeCl2. In this type, liquid-CO2 homogenizes to vapor in the temperature range 396 to 410°C (Figure 14(c)). A general increase in values is observed for type B inclusions in comparison to type A. Values of of type B inclusions fall within a narrow range from −2.4 to −1.8°C (Figure 14(d)) corresponding to low salinity, between 3.9 and 2.9 wt.% NaCl equivalent. These fluids consist of mixtures of H2O and CO2, the most common fluids released during metamorphic reactions [61].

9. Discussion

These deposits have a number of distinct features [15, 1820, 62]: (1) an epigenetic origin, formed at moderate to high temperatures (200°–350°C), (2) volcanic, sedimentary, and metamorphosed host rocks, (3) widespread alkali-rich alteration halos, (4) the presence of hydrothermal-metasomatic magnetite and/or hematite ores, (5) a common association with calc-alkaline, I-type granitoid, (6) the presence of deep seated hydrothermal replacement and shallower brecciated textures, and (7) hydrothermal fluids with low to high salinity and CO2 and other gaseous components and also chloride complexes for iron mobilization. This indicates that no single genetic model can be developed for hydrothermal Fe-oxide deposits [17, 18, 2225, 63, 64].

There are some similarities between the Qatruyeh iron ores and hydrothermal Fe-oxide deposits such as Andean iron oxide and Iron springs Utah deposits [15, 18]; these are (1) occurrence of ore bodies along a subduction-related active continental margin, (2) abundances of magnetite with lesser hematite, (3) the presence of breccia and replacement textures, (4) iron oxides with low Ti contents, (5) uranium anomalies, and (6) chemical analyses results of magmatic rock samples of the Qatruyeh similar to those of calc-alkaline plutonic rocks in Andean iron belt (Table 1). At the Qatruyeh area, titanium in both types of magnetite ores was depleted like to other hydrothermal iron deposits. It was concentrated in early magmatic ferromagnesian minerals in mafic rocks and hence was depleted in the residual late magmatic hydrothermal fluids. Layered magnetite and hematite ores related to massive magnetite in the Qatruyeh area generally have anomalous, though uneconomic, uranium grades. The relatively high values of uranium could be related to uranium-rich magmatic fluids from causative intrusions, particularly magnetite-series intrusions such as the Chah Anjir plutonic rocks at the Qatruyeh area. Alternatively, uranium could be contributed by a second fluid mixing with hydrothermal fluids at the site of mineralization or it may have been leached from host rocks [65]. The hydrothermal systems within intermediate to mafic calc-alkaline sequences, such as the subduction-related deposits of the Andean iron belt, would be expected to have low uranium contents. Therefore, uranium-rich iron oxide systems with relatively low copper values (e.g., the Qatruyeh area) suggest the presence of iron oxide uranium deposits and may be important exploration targets [65]. On this basis, the intermediate intrusions would be considered as a potential source of mineralizing components for the iron ore deposits.

HFSEs were immobile in the metacarbonate rocks, whereas large-ion lithophile elements were mobile during hydrothermal alteration and associated iron mineralization.

According to major oxide analyses of magnetite, the high Na2O content in magnetite is called alkali iron effect [15] and the system, such as FeO + Na2O + K2O + Al2O3 + MgO and FeO + MnO + TiO2 + CaO + P2O5 + SiO2, kept iron in solution at temperatures significantly lower than the melting point of magnetite at rather shallow depth [66]. The role of other volatile constituents such as CO2 and F is also important in this process [67]. The PCO2 at the Qatruyeh may have been raised by the assimilation of dolomitic-marble host rock by iron-rich hydrothermal fluids, as is evidenced by the presence of carbonate xenoliths in the magnetite ore.

The REE distribution in hydrothermal minerals is often used for discussing ore genesis. It is influenced by chemical-crystallographical constraints, kinetics of crystallization, - conditions, and the composition of the hydrothermal fluids [68]. In addition, it is necessary to take into account the effect of subsequent metamorphism and hydrothermal alteration when using REE data in ore genetic discussions. REEs geochemical investigations on the metacarbonate host rocks were shown that these elements were moderately mobilized due to metasomatic-alteration in the study area. Alderton et al. [69] stated that REE mobilization in fluids is controlled largely by (1) the REE concentrations of the reacting minerals; (2) the relative stability of these minerals against the fluid; (3) the availability of suitable crystallographic sites within the secondary minerals to accommodate released REE during alteration; (4) the REE concentrations of the hydrothermal fluid; and (5) the ability of the fluid phase to remove REE from the system.

Besides, the REE patterns for magnetite, as well as their geochemical character, show a clear difference between the high-grade and low-grade ores. The similarity of REE patterns in the layered magnetite ores and magmatic rocks is consistent with this conclusion that the layered magnetite was formed from magmatic derived hydrothermal fluids. On the contrary, different REE distribution patterns in massive magnetite indicate that these ores formed by dissimilar agents and processes at different times (e.g., by metamorphic fluids derived from dehydration and decarbonation reactions during prograde metamorphism). The metamorphic grade of the Qatruyeh area has been determined to epidote-amphibolite facies, which according to Apted and Liou [70] took place between 575 and 675°C (at 7 kb) and could consequently have produced hot fluids enough for magnetite mineralization. Also, Magnetite is not normally a major REE-bearing mineral; due to the small ionic radius of Fe therefore it has a very little influence on the REE budget in the hydrous systems. In contrast, the silicate partition coefficients for the REE are highest in such systems [71]. LREEs enrichment in low-grade ores is opposite to what one should expect based on the small ionic radius of the iron site. Schock [72] cited analyses of magnetite contained Ca2+ ions in its structure. He further noted the plausibility that REE ions with comparable ionic radii to Ca2+ (7- and 9-fold coordination with ionic radii of 1.06 and 1.18 Å, resp.) could also enter the same lattice sites. Frietsch and Predhal [73] according to de Sitter et al. [74] have showed substitution of Ca2+ ions for ferrous ion in the magnetite structure, which can contain up to 6.67 wt.% Ca2+, corresponding to the composition (Fe2.8Ca0.2)O4. The large ionic radius of Eu2+ (7- and 9-fold coordination with ionic radii of 1.20 and 1.30 Å, resp.) discriminates against, is in contradiction to incorporation into the magnetite. Accordingly, the negative Eu anomaly in low-grade magnetite indicates Eu enrichment in the remaining hydrothermal fluids after iron ore formation. Divalent Eu will thus be stable over a wide range of pH values at the magnetite-hematite buffer and moderate to high temperature hydrothermal systems, due to its large ionic radius. But, trivalent Eu will become the stable species when the system cools or becomes extremely oxidized [75]. The positive Eu anomalies in the Qatruyeh hematite ores reflect cooling of the Eu-enriched hydrothermal fluid subsequent to iron ore deposition.

The difference in Cl contents among layered magnetite, massive magnetite, and hematite can be explained by metal-ligand complexation. Zhu and Sverjensky [76] demonstrated that partitioning of F and Cl between minerals and hydrothermal fluids is a strong function of temperature, pressure, pH, and fluid composition. Chloride solutions are very adequate to transport enough metals to form deposits [77]. According to these authors, the Cl/F ratio increases with increasing temperature and pressure. The Cl ion is a high intermediate Lewis acid ligand and will affect the solubility of divalent metals (e.g., Fe2+) in hydrothermal solutions, forming a range of iron minerals such as sulfides (pyrite and chalcopyrite) and carbonates (siderite) [7880], as observed in the Qatrureh area. Most experimental studies indicate that FeCl+ and FeCl2 are the main complexes involved in the hydrothermal transport of ferrous iron, especially at high temperatures (>300°C). At the Qatruyeh, Cl concentration in the hematite ores relative to layered magnetite decreases with respect to increasing pH and decreasing temperature due to interaction between dolomitic marbles and primitive acidic-high temperature hydrothermal fluids. Therefore, a higher pH of the fluids should indicate higher concentrations of Cl in the layered magnetite relative to hematite. High Cl value records related to Fe in the Qatruyeh iron ores are good evidences for mentioned complexes in Cl-F-Fe triangular diagram (Figure 15(a)). From these data, it is clear that lower pH, higher concentrations of chloride complex as a ligand, or higher temperatures usually increase the capability of hydrothermal fluids to carry sufficient Fe and REE to explain the observed -enrichments in all layered magnetite relative to massive magnetite and hematite ores in Cl-F- trigonal diagram (Figure 15(b)).

Field observations and microscopic studies indicate two different alternative genetic concepts for ore mineralization in the Qatruyeh area. A syn-metamorphic or premetamorphic origin for the massive magnetite is supported by our observations on pebbles of the Jurassic Deh Vazir conglomerate. Magnetite recrystallization and the presence of ductile deformation textures are clear evidence for origin of layered magnetite and hematite ores later to peak of Cimmerian metamorphism followed the continental collision events.

Fluid inclusion study has proven to be an important tool for constraining the physicochemical conditions of hydrothermal mineralization [81]. The microthermometric results and the moderate salinity of the studied fluid inclusions suggest that mineralization at the Qatruyeh occurred under mesothermal conditions. At the Qatruyeh area similar to other mesothermal deposits, ore forming system is accompanied by both replacement and brecciated textures with the latter more common in shallower deposits and well developed alteration patterns [62], with sodic-calcite alteration in depth and the propylitic-hematitization haloes in shallower depths [11, 53, 62]. Hematite and quartz (± pyrite ± chalcopyrite) veins accompanied the late stage of each main ore mineralization in the Qatruyeh area.

The fluid inclusion study data also indicate that the ore forming fluids at the Qatruyeh were CO2-bearing, low to moderately saline, NaCl-MgCl2-CaCl2-H2O rich fluids. The source of the fluids (magmatic or metamorphic origin) could not be unequivocally determined from the available data; nevertheless, the distribution of fluid inclusions in ()-salinity diagram [34], with lines of constant fluid density imply cooling and increasing of pH, boiling, and fluid mixing, may play an important role in the iron ore deposition in the area (Figure 16). Cooling and increasing of pH are the main factors affecting precipitation of layered magnetite minerals. However, other processes such as fluid mixing and boiling also seem to be effective on such a precipitation. The precipitation of layered magnetite ores and actinolite was occurred by cooling and an increasing of pH due to reaction of ore bearing fluids with dolomitic-marble rocks. For example, in a simplified reaction:

An increase in pH will move the reaction to the right side, resulting in the precipitation of magnetite and actinolite [78].

Boiling results in the production of vapor and loss of H2O, and other volatile species occur. Salinity variations (from 3.5 to 15 wt.% NaCl equiv.) in the fluid inclusions data can be produced by boiling. This accompaniment to cooling would result in production of a relatively high-density liquid phase. Figure 16 shows that the density of mineralizing fluid at the Qatruyeh ranged from 0.6 to 0.93 g cm−3. In addition, some fluid inclusions that contain relatively high concentrations of CO2 in quartz can be constructed assuming a model involving a fluid rising adiabatically and mixing with metamorphic fluids.

10. Conclusions

Based on ore mineral types and alteration assemblages, the mineralization processes have been divided into the following stages: (I) layered magnetite ore stage, (II) sulfide ore stage, and (III) hematite ore stage. Stage I involved the formation of the layered magnetite ore, high temperature hydrothermal activity, replacement textures, and assemblages of actinolite + quartz + titanite+ siderite + paragonite ± tourmaline. Stage II partially overlapped with stage I and includes hematite ores with epidote + chlorite + sericite + quartz ± magnetite assemblages. Propylitic zones with associated hematite ores and brecciated textures are shown by low-temperature oxidizing- hydrothermal fluids in shallow levels due to cooling and fluid mixing at this stage. Stage III involved the formation of hematite and quartz veins and some minor sulfide minerals as chalcopyrite and pyrite. Therefore, the evolution of hydrothermal fluids varies due to fluid-dolomitic-marble interaction from acidic to alkaline affinity during iron mineralization.

The low-grade magnetite ores and diorites whole rock REE patterns are comparable to each other, indicating that the bulk of REE concentration in low-grade magnetite is inherited from precursor diorite during the replacement process that was obvious in ore textures. Besides, trace elements geochemical investigation of the metamorphosed carbonate host rocks shows that(1)HFSEs, such as Th, Hf, Zr, and Ta, were immobile during hydrothermal alteration and iron mineralization.(2)REEs were moderately mobile, whereas large ion lithophile elements such as Ba, Sr, Rb, and Cs were depleted during the alteration and iron mineralization.

From the available data, it is not obvious whether the massive magnetite formed from an iron-rich melt or through metasomatism but the presence of magnetite pebbles in the Jurassic Deh Vazir conglomerate indicates that a part of ore mineralization occurred syngenetic or premetamorphic.

Fluid inclusion data suggest that low to moderate saline fluid (3.5 to 15 wt.% NaCl equiv.) was possibly responsible for the development of low-grade magnetite ore that indicates episodes of hydrothermal activities in mesothermal systems for Fe occurrences in this region. The microthermometric data reflect variable amounts of magmatic fluids and low salinity metamorphic fluids. The temperature during the main mineralization event was between 300° and 410°C. Increase in pH due to reaction of ore fluids with dolomitic-marble rocks and effervescence can be effective mechanisms of ore precipitation. In addition, mixing with cooler probably metamorphic fluids had the effect of decreasing the salinity and temperature of the equilibrated magmatic fluid, making it more buoyant and capable of moving out from the apophyses into the host rocks, causing metasomatic reactions and precipitating low-grade magnetite ores from chloride-bearing hydrothermal fluids.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publishing of this paper.

Acknowledgments

The authors would like to thank the research council of Shiraz University for financing this research (Grant no. 88-GR-SC-83). They express gratitude to professor Nigel. J. Cook (University of Adelaide) who provided many insightful and constructive comments that greatly improved the paper. They would like also to thank professor Warren. C. Day (US Geological Survey) for criticism of an early version of the paper.