Evolutionary Model of the Sepid Sarve Manto Type Copper Mineralization, Doruneh Fault Volcanic Plutonic Belt (NW Bardeskan, Central Iran): An Integrated Geological, Geochem-ical, Fluid Inclusion and Stable O–S Isotope Study

1. Introduction

In the uppermost levels of the Eocene volcanic sequences in Iran, numerous copper deposits and mineral occurrences have been identified (e.g., Govin Cu-Au deposit [1], Kuh Pang Copper deposit [2], Masjed Daghi Porphyry Cu-Au deposit [3], Gheysari Copper deposit [4]). These occurrences have been variously described in previous studies as hydrothermal, vein-type, stockwork, volcanic-related, and others.

Mesozoic manto-type Cu deposits are classified into two categories: sed-imentary-hosted deposits (e.g., Cerro Negro) and volcanic- hosted deposits (e.g., Mantos Blancos and Buena Esperanza) [6]. Three main theories have been proposed for the genesis of these ore deposits: (a) syngenetically formed volcanogenic deposits [5,7,8]; (b) epigenetic deposits formed by magmatic emanations from associated plutonic intrusions [9,10]; and (c) epigenetic deposits formed by fluid–rock interaction involving Cu-bearing host rocks [11,12,13,14,15].

Most Manto-type copper deposits are hosted within stratigraphic sequences predominantly composed of andesitic to basaltic rocks. These sequences correspond to the Mesozoic era (in Chile and Iran), the Proterozoic to Triassic (in North America), the Late Silurian to Early Devonian (in China), and the Cenozoic (in Iran) (Figure 1). The presence of Manto-type copper mineralization in Iran was first identified and reported by Samani at the Veshnaveh deposit (Qom Mountains) in 2000. Since then, the possibility of similar occurrences elsewhere in the country has received increased attention [16]. A comparative analysis of copper metallogeny in Chile and its correlation with the geological characteristics of known Iranian deposits suggests that a specific type of copper mineralization, known as Mantos, may exist in regions such as Torud–Abbasabad, Golestan, Veshnaveh (Qom Mountains), Anarak, and other localities [16]. Accordingly, the proper identification and understanding of these types of mineralization can offer a new perspective in the study of copper geology and resource potential in Iran.

To date, more than 30 major Manto-type copper deposits and numerous copper occurrences ranging in age from Cretaceous to Eocene have been identified in Iran(e.g., Veshnaveh deposit[17], Keshtmahaki[[18], Koshkoiye[19], Vorezg-Qayen copper deposit[20], Abbas Abad copper deposit[21], Mari Cu(Ag) Manto copper deposit[22], Kahak Manto-Type copper deposir[23], Narbaghi deposit[24], Yamaghan Manto-Type Cu(Ag) deposit[25] and Madan Bozorg Cu deposits,Abbas Abad[26]). Some of these deposits were historically exploited in ancient times—for example, Veshnaveh deposit [17] and Narbaghi deposit [24] (Figure 1).

Manto-type deposits in Iran are typically associated with andesitic lava flows. The Eocene host rocks (Eocene volcano-sedimentary sequences) are the most common host lithologies for manto-type copper–silver deposits in Iran. These deposits are mainly concentrated within the Urumieh–Dokhtar Magmatic Arc (UDMA) [24], the Alborz Magmatic Arc (AMA), the Sabzevar Zone, and the Lut Block. Approximately 85% of Iran's manto-type copper deposits occur within the UDMA, the Sabzevar Zone (north of the Doruneh Fault), and the AMA. This zone is composed of extensive magmatic rocks and includes the largest Iranian ophiolite m´ elange that is considered as remnants of a Cretaceous back-arc oceanic basin that formed by hy perextension in the Iranian plateau due to the initiation of NeoTethys subduction [27,28]. the ore deposits of the Sabzevar Zone can be separated into two groups, based on their different metal association, spatial distribution, and geodynamic setting. The Cretaceous mineralized orebodies consist of Cr deposits associated with serpentinized peridotites, VMS deposits, and Mn ore deposits in volcano-sedimentary sequences [29]. The Cenozoic mineralization in the Sabzevar Zone includes Cu- Au porphyry, red-bed, and Manto-type Cu deposits. The stratigraphic sequence and host rock lithotypes of the Manto-type deposits of the Sabzevar Zone are very similar to the Cu-(Ag) deposits of the AMA, which are also hosted by the Eocene volcanic rocks. Notable examples include the Abbasabad deposit in the Alborz range, Keshkuiyeh in the UDMA, and Sepid-Sarve and (e.g., Abri, Rahbari, and Cheshmeh Marzeyeh cu deposit[30], Bornaward plutonic complex[31], Baharieh copper deposit[32], Kal Abri copper deposit[33], Nasim copper deposit[34], Sharifabad copper deposit[35], Madan Bozorg[36], Grik and Gukhab[37], Golcheshmeh[38], Cheshmeh-Hadi[39]) in the southern Sabzevar Zone (north of the Doruneh Fault). The majority of the deposits are being actively explored. The Cheshmeh Marzeyeh cu deposit [30], Kal Abri copper deposit [33], Nasim copper deposit [34], Madan Bozorg[36] and Cheshmeh-Hadi[39] deposits are presently mined (Figure 1).

Vein-type mineralization within andesitic lava flows highlights the critical role of stratigraphic control in the formation of Iran’s manto-type copper deposits. In the Sepid-Sarve copper deposit, located ~67 km west of Bardeskan and 17 km from Doruneh village in southwestern Khorasan Razavi Province, mineralization is hosted by Eocene volcanic rocks and restricted to specific stratigraphic horizons. The deposit lies within the Doruneh sheet, along the northern margin of Central Iran and the southern Sabzevar Zone, an area shaped by intense Cenozoic magmatism and activity along the Doruneh Fault. This tectonic setting has exposed lithological units ranging from Precambrian to Quaternary. Despite the high potential of volcano-sedimentary sequences for manto-type copper mineralization, limited understanding of their genetic models and mineralization controls has led to the abandonment of many deposits before achieving economic results.

This study characterizes the geological setting and mineralization style of the Sepid-Sarve deposit, and proposes a large-scale mineralization model integrated with the tectono-magmatic evolution of the region, based on newly acquired field and analytical data. The results provide a refined exploration framework for identifying new manto-type copper targets in similar volcano-sedimentary settings, improve the understanding of ore-forming processes in the Sabzevar Zone, and offer practical guidelines for more effective exploration and resource development strategies in Iran’s metallogenic provinces

2. Materials and Methods

2.1. Geochemical analysis

As part of this study, and following an extensive literature review, a total of 134 samples were collected during fieldwork from intrusive and volcanic igneous units, as well as from mineralized veins and altered zones of the Sepid-Sarve copper deposit. Sampling was conducted from the surface down to a depth of 56.2 meters to facilitate petrographic and lithological investigations. Corresponding thin sections were prepared for detailed analysis. Among these, 20 intrusive rock samples exhibiting minimal alteration were selected for whole-rock geochemical (lithochemical) analysis. Major elements in these samples were measured using X-ray fluorescence (XRF), following the processes of crushing, pulverizing, drying, homogenizing, powder preparation, pressing with boric acid, and finally fusion with lithium metaborate or tetraborate. For the analysis of rare earth elements (REEs) and trace elements, 0.2 grams of each sample was digested with diluted nitric acid and subsequently analyzed using inductively coupled plasma mass spectrometry (ICP-MS), based on ionization and separation of ions according to their mass-to-charge ratios. The analytical procedures were carried out at the Zar Azma Laboratory in Tehran.

2.2. Fluid inclusions measurements

Five double-polished thin sections with a thickness of 80–100 µm were prepared from quartz and calcite (carbonate) veins and veinlets. Petrographic observations of fluid inclusion evidence were conducted, followed by microthermometric analyses on 40 fluid inclusions. These studies were carried out using a Linkam THMSG600 heating–freezing stage with a temperature range of −196°C to +600°C and a precision of ±1°C, at the Department of Geology, Tarbiat Modares University, Tehran. Temperature calibration was performed using Synflinc synthetic standards. During the freezing and heating procedures, the rates were controlled within the range of 1–10 ◦C/min, before being decreased to 0.5–1 ◦C/min in the vicinity of phase transition points. Salinities (wt.% NaCl equiv.) and densities (g/cm³) of the aqueous (NaCl-H₂O) inclusions were determined using the computer program HOKIEFLINCS_H₂O-NaCl Excel developed by [42] based on the final ice melting temperatures.

2.3. Stable isotope measurements

2.3.1. Oxygen Isotope Analysis

Two quartz and calcite samples were selected from the same vein locations where fluid inclusion samples were collected, for stable oxygen isotope (δ¹⁸O) analysis. The samples were first crushed to mesh 60 at the University of Tehran and then manually purified under a binocular microscope to achieve a purity of over 99%. After separation, the samples were powdered, sieved to remove contaminants, and rinsed with ethanol to ensure cleanliness. Isotopic analyses were performed at the Stable Isotope Research Center, Arak University, Iran. For δ¹⁸O analysis in solid samples, high-temperature pyrolysis was conducted at 1450°C using an elemental analyzer, where all oxygen atoms in the molecular structure were converted to carbon monoxide (CO) gas in the presence of carbon black and glassy carbon. The resulting CO gas passed through a series of purification systems, including a water trap to remove moisture and a purge-and-trap column to eliminate CO impurities. After complete removal of interfering gases, the desorption column was heated to 150°C, and the purified CO gas was transferred through a secondary water trap into the isotope ratio mass spectrometer (IRMS). For calibration and calculation, the isotopic composition of the reference CO gas was determined by comparison with international or laboratory secondary standards. The mass ratio of 30/28 was measured to determine the δ¹⁸O values. In this process, δ¹⁸O is measured by converting the sample to CO gas via pyrolysis at 1450°C, followed by purification and subsequent analysis by IRMS, using the 30/28 mass ratio and calibrated against international standards. δ¹⁸O values exhibit a precision of 0.1‰ and are reported in reference to V-SMOW (Vienna Standard Mean Ocean Water, δ¹⁸O = 0, with uncertainties of 0.02‰; International Atomic Energy Agency, (2017).

2.3.2. Sulfur Isotope Analysis

Three representative chalcocite samples from mineralized veinlets were selected for sulfur isotope (δ³⁴S) analysis. The samples were first powdered and sieved to obtain grain sizes between 0.1–0.5 mm. Clean mineral grains were handpicked under a binocular microscope to ensure purity. For isotopic analysis, the samples were combusted at 1150°C in an elemental analyzer, producing sulfur dioxide (SO₂) gas. The generated SO₂ gas was directed through a purification system consisting of a water trap to remove moisture and a purge-and-trap column to eliminate residual impurities. The adsorption column was subsequently heated to 220°C to release the accumulated SO₂ gas, which was then transferred into the isotope ratio mass spectrometer (IRMS). Within the IRMS, the sulfur isotopic composition was determined by measuring the mass ratio 66/64, corresponding to the ³⁴S/³²S ratio of the sample. To ensure the accuracy and reproducibility of the measurements, repeated calibrations were performed using the IAEA-S-4 international standard (δ³⁴S = +16.9 ± 0.2‰ VCDT) and a laboratory secondary standard (δ³⁴S = -6.3 ± 0.2‰ VCDT).

3. Geology Setting

3.1. Regional Geology and Tectonic Setting

According to the tectonostratigraphic framework of Iran [41,43], the study area lies within the Sabzevar Zone of the Central Iranian terrane. The Bardeskan region can be subdivided into two distinct domains: north of the Doruneh Fault, the terrain is mountainous and morphologically diverse; south of the fault, except for the southeastern hills associated with the Ozbak-Kuh outcrops, lies a subsiding basin with an average elevation of ~851 m above sea level. This basin is extensively covered by Quaternary deposits and exhibits geomorphic features such as alluvial terraces, alluvial fans, mudflats, and salt flats.

The most prominent geotectonic feature is the east–west-trending Doruneh Fault, which separates the northern mountainous domain from the southern plains. The Sepid-Sarve area lies within the 1:250,000 Bardeskan quadrangle [44] and the central portion of the 1:100,000 Doruneh geological map [45]. Two major geological characteristics define this region: (1) intense Mesozoic folding, metamorphism, and plutonism; and (2) extensive Cenozoic volcanic activity. Central Iran, however, is not a homogeneous terrane but a mosaic of distinct structural blocks, each with unique geological features resulting from extensive faulting.

Eocene volcano-sedimentary rocks are widely exposed throughout the region, forming a belt with substantial potential for manto-type copper mineralization. In the Sepid-Sarve deposit, mineralization is hosted by Eocene volcanic rocks and restricted to specific stratigraphic horizons. The Doruneh quadrangle, along the northern margin of Central Iran and the southern Sabzevar Zone [40,41], is particularly significant due to intense Cenozoic magmatism and the activity of the Doruneh Fault. These tectono-magmatic processes have exposed a wide range of lithologies, from Precambrian to Quaternary.

The occurrence of vein-type mineralization within andesitic lava flows underscores the critical role of stratigraphic control in the localization and formation of metals in Iran’s manto-type copper deposits. Stratigraphic positioning is therefore considered a key factor in ore formation. Furthermore, the lateral and vertical continuity of mineralization is strongly influenced by host-rock lithology, extensional fault systems, and the degree of fracturing.

3.2. Ore Deposit Geological Characteristics

The Sepid-Sarve copper deposit is hosted within an Eocene volcano–sedimentary sequence in the southern Sabzevar Zone. This sequence comprises intermediate to mafic volcanic units (porphyritic to megaporphyritic andesite, pyroxene andesite, and basalt), pyroclastic rocks (tuff and agglomerate), and sedimentary strata including fossiliferous limestone (bearing Nummulites, Plesiopoda, and brachiopods), calcareous sandstone, detrital limestone, conglomerate, siltstone, shale, gypsiferous marl, and alluvial terrace deposits.

A distinctive feature of the region is the widespread development of Cretaceous volcanic–tuffaceous rocks and Eocene volcanic and pyroclastic units, reflecting intense volcanic activity. In the central part of the study area, Middle Eocene sediments host a thick volcanic–pyroclastic complex dominated by andesite, pyroxene andesite, basaltic andesite, and basalt, typically dark green to black in color and often displaying bedding and banding. In some sections, green to red agglomerates are interbedded within this volcanic package. Copper mineralization, mainly malachite, occurs along fractures within the andesites (Figure 2A). Overlying this complex is a pyroclastic sequence of agglomerate, tuff, lapilli tuff, and andesitic tuff with colors ranging from red to green.

Upper Eocene sediments consist of alternating sandstone, siltstone, silty marl, and thin-bedded limestone, generally light green to cream in color. The Upper Eocene volcanic–pyroclastic succession in the Sepid-Sarve deposit area can be divided into four units. The first unit is composed of andesite, porphyritic andesite, and pyroxene andesite (green to dark gray), with amygdaloidal textures formed by alteration of phenocrysts (Figure 2A). Vesicles are typically filled with secondary calcite. Above this lies a volcanic complex of basalt and basaltic andesite (black), overlain by Oligo–Miocene sedimentary deposits, representing the final phase of Eocene volcanism. These volcanic rocks exhibit gentle morphology and are locally vesicular, with vesicles filled by calcite, silica, and occasionally zeolite.

Several fault systems are present in the study area, including E–W and N–S trending faults, as well as dextral strike-slip faults trending NW–SE and NE–SW, the latter being dominant. In places, these faults displace mineralized zones. In the easternmost part of the deposit, overturned limestone units indicate intense tectonic activity. Here, multiple NW–SE strike-slip faults have disrupted the limestone layers and displaced copper mineralization hosted at the contact between limestone and volcanic units (Figure 2B).

3.2.1. Conglomerate

The EOL^cv unit, representing the oldest lithological formation exposed in the Sepid-Sarve copper deposit area, is widely exposed in the southern, central, and northern parts of the region and constitutes one of the most prominent geological units with significant thickness (Figure 3H,I). One of the most important features of this unit is the occurrence of ore zones at the contact between this unit and a Nummulitic Eocene limestone (Figure 3A,B,D). Outcrops of this unit exhibit a hilly to rugged topography, particularly in the southern part, where it forms the highest elevations of the area and structurally defines the outer limb of an anticline.

Field exposures display varied coloration, ranging from brown and red to dark gray, and are composed of a volcano-sedimentary conglomerate sequence interbedded with volcaniclastic rocks and volcanic lava flows. The conglomerates, forming the main component of this unit, consist of coarse to medium-grained, poorly sorted clasts with moderate to good roundness. The clasts are predominantly composed of older volcanic rocks, including andesitic and basaltic-andesitic fragments (Figure 3E,H).

Intermediate to mafic lava flows, including andesite, basaltic andesite, and pyroxene andesite, are present within this unit, particularly thickened in the eastern parts (Figure 3I).

3.2.2. Andesitic and Basaltic Andesite

The Eocene andesitic and basaltic andesite volcanic units show mineralogical, textural, and lithological diversity. These lava flows are conformable and parallel to adjacent units. Fresh andesites are dark gray, weathering to grayish-brown, with some areas stained reddish by iron oxides (Figure 3F,G). Hand samples reveal porphyritic, megaporphyritic, vesicular, amygdaloidal, and flow textures, featuring phenocrysts of plagioclase, pyroxene, and opaque minerals. Vesicles are often filled with secondary silica, calcite, and chlorite minerals, giving amygdaloidal and vesicular textures.

3.2.3. Pyroxene Andesite

The pyroxene andesite unit appears as low-relief hills with a light gray color and displays porphyritic and amygdaloidal textures. Hand specimens contain rounded pyroxene crystals and coarse plagioclase crystals (27–36%). Secondary minerals like chlorite and calcite are present as veinlets and vesicle-fillings, with carbonate veins measuring 1.5–2.5 mm clearly visible.

3.2.4. Hornblende-Pyroxene Andesite

The hornblende-pyroxene andesite is also light gray, with vesicles filled by celadonite, chlorite, calcite, and minor zeolites, and shows porphyritic and amygdaloidal textures. Plagioclase phenocrysts (43–45%) along with hornblende and pyroxene are easily observed. The hornblende basaltic andesite is dark gray, porphyritic, and amygdaloidal, with visible plagioclase (22–26%) and hornblende (4–6%) phenocrysts. The hornblende-pyroxene basaltic andesite, found in mountainous terrain in the eastern area, also displays porphyritic texture with plagioclase (27–29%), pyroxene (5–6%), and hornblende (4–6%). The pyroxene basaltic andesite, appearing as low-relief hills, contains plagioclase (24–26%) and pyroxene (5–6%). The olivine-pyroxene basalt is dark green to black, fine-grained, and includes plagioclase (16–18%), pyroxene (5–7%), and olivine (12–14%) phenocrysts.

3.2.5. Basalt, Olivine Basalt, and Microgabbro

Basalt, olivine basalt, and microgabbro units occur as low-relief hills with smooth and eroded surfaces across the region (Figure 5A,C). In hand specimens, these rocks appear dark green, dark gray, to black in color. They are primarily composed of phenocrysts of olivine and pyroxene, and commonly display porphyritic, microporphyritic, amygdaloidal, and vesicular textures (Figure 5D,K). Most olivine crystals in these basalts exhibit iron oxide staining and iddingsitization along grain boundaries and fracture zones.

In certain areas of the Sepid-Sarve deposit, this unit contains weak mineralization. The gabbro-basalt rocks have medium-grained textures and consist of plagioclase, olivine, pyroxene, and minor alkali feldspar. Olivine crystals in the groundmass are frequently altered to iron oxides. The plagioclase composition is sodic, and in association with alkali feldspar, suggests disequilibrium with olivine (Figure 5E,I), indicating evidence of magma mixing. Chlorite appears as minor veinlets in parts of the rock. Partial alteration of olivine to iddingsite and secondary iron oxides is evident.

Microscopically, the gabbro-basalt rock contains 45–50% phenocrysts. Of the total phenocryst volume, 60–65% are pyroxene and 35–40% are olivine. Pyroxene crystals range from 1 to 4 mm, and olivine crystals from 0.5 to 3 mm in size. Plagioclase crystals, occasionally observed, measure about 0.5–1 mm and are finer than the mafic phenocrysts (Figure 5). Both orthopyroxene and clinopyroxene are present (Figure 5). Some vesicles are filled with secondary quartz and carbonate minerals. Occasionally, pyroxene appears as large megacrysts (Figure 6E,E-1,E-4).

3.2.6. Agglomerate

The volcaniclastic units are represented by andesitic agglomerates that occupy a significant portion of the area and are stratigraphically interbedded with lava flows. These rocks show weak to moderate hydrothermal alteration, including chloritic, hematitic, and, in some locations, carbonate and argillic alterations. Agglomerates are brown to reddish in color and cover extensive areas. In hand specimens, some show pyroclastic textures and are composed of andesite to pyroxene andesite.

The majority of the copper mineralization in the Sepid-Sarve deposit occurs at the contact of these units with the overlying limestone (Figure 4F,H). Copper mineralization appears as both oxide (malachite) and sulfide (chalcocite) phases. The crystal tuff has an amygdaloidal texture and is composed of plagioclase and pyroxene, with secondary minerals including iron oxides, carbonates, and chalcedonic quartz. The groundmass is iron-rich. Intense alteration is evident, with iron oxides forming a substantial portion of the rock and the vesicles being filled with secondary silica and carbonate minerals. This rock likely represents an oxidized crystal tuff. Pyroxene and scattered plagioclase are observed. A large portion of the rock consists of porous iron oxides, with vesicles filled by secondary minerals.

3.2.7. Limestone

Thin-bedded, cream-colored limestones (OLM^lm), interbedded with marl and siltstone and rich in Nummulite and bivalve fossils, occur in multiple synclines across the Sepid-Sarve area (Figure 4B,C,I). These carbonate units include silty, sandy, and detrital limestones. The main limestone body, trending east-west in the central region, is conformably interbedded with volcanic and volcaniclastic rocks.

This limestone is a fossiliferous micritic packstone with abundant fossils and a fine-crystalline matrix. Copper mineralization often occurs at its contact with volcanic units. Carbonate breccias here consist of angular fossil-rich clasts in a calcitic matrix (Figure 6A), with some fragments containing quartz and plagioclase grains. Large benthic foraminifera like Alveolina are common. Iron oxides sometimes replace the margins of carbonate clasts (Figure 6A,A4).

This limestone is a key exploration target, as major copper mineralization in Sepid-Sarve is associated with its contact with overlying agglomerates.

Petrographically, the limestone includes lumachelle, fossiliferous micritic (biomicrite), sparitic (biosparite), and microsparitic (biomicrosparite) types, showing coarse textures and abundant fossil debris. Minor feldspar and silica detritus are present, with a matrix of micrite, sparite, and iron oxides. Secondary silica fills fractures, and abundant carbonate fossils (mostly Nummulites) are present alongside minor feldspar and quartz grains (Figure 6D,D4).

3.2.8. Marl and Siltstone

The OLM^sm unit consists of thinly bedded marl and siltstone layers, ranging from light gray to cream, overlying the fossiliferous limestone. It forms hilly terrain, with greater thickness and lateral extent in the western sector compared to the east (Figure 4A–C, E, I). This unit is conformably overlain by a red marl layer. Notably, limestone beds analogous to OLMlm occur within the northern part of the western sector, hosting malachite mineralization, and have been structurally emplaced into the marl via thrust faulting.

The OLM^m unit comprises thinly bedded red to gray marl and siltstone overlying the gray-cream marl. It forms hilly topography, is relatively thin but more extensive in the eastern region, and is unconformably overlain by a brown conglomerate (Figure 4A).

The M^c unit is characterized by gray to brown polymictic conglomerates with clasts from pebbles to boulders. Sedimentologically, it is poorly sorted with moderate clast roundness. Interbedded sandstone layers (Mts unit), more prevalent in the west, accompany this unit. The Mc forms hilly relief and exhibits significant thickness, particularly in the east. In some areas, it is unconformably overlain by a gypsum-bearing marl (Figure 4A, D).

The M^ts unit comprises gray to green sandstones interlayered with tuffs, occurring interbedded within the Mc conglomerate. It is more extensive and thicker in the western sector, exhibiting low-relief hills and relatively minor thickness (Figure 3A, B, D).

The M^gy unit includes marls with transparent gypsum layers, gray to white in color, present in both eastern and western sectors (Figure 3C, D). It locally overlies older units discontinuously, with wider distribution in the east and beyond the study area. The unit hosts multiple high-quality gypsum layers currently under extraction, displaying rugged topography and limited thickness (Figure 4B, E).

The M^gy unit includes marls containing transparent gypsum layers, gray to white in color, and is observed in both the eastern and western sectors (Figure 3C,D). It locally and discontinuously overlies older units. The unit is more widespread in the eastern part and is also extensively developed outside the study area. It contains multiple high-quality, bright gypsum layers that are currently being mined. This unit displays relatively rugged topography and has limited thickness (Figure 4B,E).

Figure 3. (A) General view of the lithological units in the Sepid-Sarve copper deposit area. (B) View of the eastern part of the deposit showing Nummulitic limestone, conglomerate, andesite, agglomerate, siltstone, marl, and gypsum units. (C) View of the marl and gypsum units in the central part of the deposit. (D) View of the eastern part of the deposit showing Nummulitic limestone, conglomerate, andesite, agglomerate, siltstone, marl, and gypsum units. (E) Outcrop view of Nummulitic limestone associated with conglomerate, andesite, and agglomerate. (F) View of the marl and gypsum units alongside conglomerate, andesite, and agglomerate. (G) Mineralized zone located at the contact between Nummulitic limestone and andesite. (H) View of the eastern part of the deposit showing Nummulitic limestone, conglomerate, andesite, agglomerate, siltstone, marl, and gypsum units. (I) Outcrop of EOL^cv(Eocene Limestone with Volcanic Intercalations) and OLM^lm (Oligo–Miocene Marl and Limestone)units with andesite in the eastern part of the Sepid-Sarve copper deposit.

Figure 3. (A) General view of the lithological units in the Sepid-Sarve copper deposit area. (B) View of the eastern part of the deposit showing Nummulitic limestone, conglomerate, andesite, agglomerate, siltstone, marl, and gypsum units. (C) View of the marl and gypsum units in the central part of the deposit. (D) View of the eastern part of the deposit showing Nummulitic limestone, conglomerate, andesite, agglomerate, siltstone, marl, and gypsum units. (E) Outcrop view of Nummulitic limestone associated with conglomerate, andesite, and agglomerate. (F) View of the marl and gypsum units alongside conglomerate, andesite, and agglomerate. (G) Mineralized zone located at the contact between Nummulitic limestone and andesite. (H) View of the eastern part of the deposit showing Nummulitic limestone, conglomerate, andesite, agglomerate, siltstone, marl, and gypsum units. (I) Outcrop of EOL^cv(Eocene Limestone with Volcanic Intercalations) and OLM^lm (Oligo–Miocene Marl and Limestone)units with andesite in the eastern part of the Sepid-Sarve copper deposit.

Figure 4. (A) General view of the limestone units (OLM^lm), marl and siltstone (OLM^m), and conglomerate (M^c). (B) Outcrop view showing limestone (OLM^lm), marl and siltstone (OLM^m–OLM^sm), conglomerate (M^c), and gypsiferous marl (M^gy). (C) Outcrop view of limestone (OLM^lm), marl and siltstone (OLM^m), conglomerate (M^c), and gypsiferous marl (M^gy). (D) View of conglomerate (Mc) and gypsiferous marl (M^gy) units. (E) Outcrop of gypsiferous marl (M^gy) and marl–siltstone (OLM^sm) units. (F) Nummulitic limestones emplaced into marl units along thrust faults at the contact with the mineralized zone. (G,H) View of copper mineralization zones at the contact between Nummulitic limestone and conglomerate, andesite, and agglomerate units. (I) View of limestone (OLM^lm), marl and siltstone (OLM^m–OLM^sm), conglomerate (M^c), and gypsiferous marl (M^gy) units.

Figure 4. (A) General view of the limestone units (OLM^lm), marl and siltstone (OLM^m), and conglomerate (M^c). (B) Outcrop view showing limestone (OLM^lm), marl and siltstone (OLM^m–OLM^sm), conglomerate (M^c), and gypsiferous marl (M^gy). (C) Outcrop view of limestone (OLM^lm), marl and siltstone (OLM^m), conglomerate (M^c), and gypsiferous marl (M^gy). (D) View of conglomerate (Mc) and gypsiferous marl (M^gy) units. (E) Outcrop of gypsiferous marl (M^gy) and marl–siltstone (OLM^sm) units. (F) Nummulitic limestones emplaced into marl units along thrust faults at the contact with the mineralized zone. (G,H) View of copper mineralization zones at the contact between Nummulitic limestone and conglomerate, andesite, and agglomerate units. (I) View of limestone (OLM^lm), marl and siltstone (OLM^m–OLM^sm), conglomerate (M^c), and gypsiferous marl (M^gy) units.

Figure 5. (A–C) Outcrop views of the gabbro–basalt unit in the Sepid-Sarve copper deposit area. (D) Olivine phenocryst associated with carbonate minerals. (E) Groundmass composed of alkali feldspar and plagioclase. (F) Pyroxene and olivine phenocrysts. (G) Pyroxene phenocryst within a groundmass that commonly contains quartz. (J) Clinopyroxene and orthopyroxene phenocrysts in the rock. (K) Olivine with fractures partially replaced by chlorite. (H) Fractured olivine phenocryst within the rock matrix. (I) Olivine phenocryst embedded in a groundmass consisting of feldspar, pyroxene, and olivine.

Figure 5. (A–C) Outcrop views of the gabbro–basalt unit in the Sepid-Sarve copper deposit area. (D) Olivine phenocryst associated with carbonate minerals. (E) Groundmass composed of alkali feldspar and plagioclase. (F) Pyroxene and olivine phenocrysts. (G) Pyroxene phenocryst within a groundmass that commonly contains quartz. (J) Clinopyroxene and orthopyroxene phenocrysts in the rock. (K) Olivine with fractures partially replaced by chlorite. (H) Fractured olivine phenocryst within the rock matrix. (I) Olivine phenocryst embedded in a groundmass consisting of feldspar, pyroxene, and olivine.

Figure 6. Overview of the Sepid-Sarve copper deposit area: (A) Packstone unit. (B) View of the mineralized zone within the conglomerate, andesite, and agglomerate units. (C) Eastern part of the Sepid-Sarve copper deposit. (A-1) Coarse pyroxene crystal set in an iron oxide-rich matrix. (A-2) Dispersed carbonate grains in an iron oxide matrix. (A-3) Vesicles filled with carbonate and chalcedonic minerals. (A-4) Cluster of several pyroxene crystals in the host rock. (E) View of the gabbro–basalt unit in the Sepid-Sarve area. (E-1) Olivine and pyroxene phenocrysts. (E-2) Groundmass with medium-sized crystals composed of plagioclase and alkali feldspar. (E-3) Secondary quartz and carbonate minerals filling vesicles. (E-4) Section of a very large (megacryst) pyroxene crystal. (D) View of the nummulitic limestone overlying the mineralized zone. (D-1) Coarse fragment of an alveolinid fossil. (D-2) Rock fragments with packstone composition. (D-3) Coarse calcite clast and carbonate fragments with iron oxides and carbonate minerals. (D-4) Alkali feldspar and plagioclase in the carbonate matrix and lithic fragments.

Figure 6. Overview of the Sepid-Sarve copper deposit area: (A) Packstone unit. (B) View of the mineralized zone within the conglomerate, andesite, and agglomerate units. (C) Eastern part of the Sepid-Sarve copper deposit. (A-1) Coarse pyroxene crystal set in an iron oxide-rich matrix. (A-2) Dispersed carbonate grains in an iron oxide matrix. (A-3) Vesicles filled with carbonate and chalcedonic minerals. (A-4) Cluster of several pyroxene crystals in the host rock. (E) View of the gabbro–basalt unit in the Sepid-Sarve area. (E-1) Olivine and pyroxene phenocrysts. (E-2) Groundmass with medium-sized crystals composed of plagioclase and alkali feldspar. (E-3) Secondary quartz and carbonate minerals filling vesicles. (E-4) Section of a very large (megacryst) pyroxene crystal. (D) View of the nummulitic limestone overlying the mineralized zone. (D-1) Coarse fragment of an alveolinid fossil. (D-2) Rock fragments with packstone composition. (D-3) Coarse calcite clast and carbonate fragments with iron oxides and carbonate minerals. (D-4) Alkali feldspar and plagioclase in the carbonate matrix and lithic fragments.

4. Discussion

4.1. Ore Mineralogy and Paragenetic Sequence

Intense compressional tectonics formed anticlines and synclines, with mineralization localized along their limbs. High-permeability conglomerates facilitated fluid flow and metal deposition. Mineralized zone thickness varies from <1 m to >10 m, occurring as veinlets, disseminations, cavity fillings, and impregnation textures (Figure 7D–F). Veinlet-style mineralization mainly occurs in andesitic and pyroclastic rocks, with limestone less involved (Figure 7G). Conglomerates, rich in volcanic and agglomerate fragments, provided porous pathways for hydrothermal fluids or metamorphic brines, enabling open-space and disseminated mineral precipitation. Copper mineralization is concentrated at contacts between pyroclastic/volcanic units and overlying carbonate sequences (marly, fossiliferous, and clastic limestones). Main minerals include sulfides (chalcocite, covellite, bornite, pyrite, rare chalcopyrite), copper carbonates (azurite, malachite), oxides (hematite, limonite, goethite, cuprite), and gangue (calcite, zeolite, quartz, clays). Secondary copper minerals and iron oxides formed from sulfide oxidation under surface conditions. Brecciated fragments (>1 cm) result primarily from penetration of chalcocite-bearing veinlets, which control mineralization patterns. These veinlets reach 1–2 mm thickness and are partly replaced by chalcocite throughout the host rock.

The veinlets are often discontinuous, with chalcocite mineralization frequently terminating within non-metallic veinlets (Figure 8D–F). Chalcocite occurs as both primary and secondary forms, the latter showing a bluish tint (Figure 9C-1-a). Thin, short chalcocite veinlets are observed locally, and fine disseminated chalcocite particles may be present but are difficult to distinguish from hematite due to similar coloration (Figure 9B-1-a, B-1-b, B-2-a, B-2-b).

Mineralization mainly occurs below the carbonate horizon and within upper volcanic and volcaniclastic units—conglomerate, andesite, and agglomerate (Figure 9A). In some areas, oxide copper minerals appear at the limestone’s upper boundary in contact with shale and marl (Figure 9B,C).

Hematite is found as disseminated granules and short veinlets. Malachite fragments and veinlets also occur, identifiable by their green color. Magnetite crystals, partially or fully altered to hematite, vary in abundance; hematitization intensity is generally low but locally more pronounced (Figure 8H–I, 8L). Magnetite appears as fine dispersed particles and as coarser subhedral to anhedral grains undergoing partial to complete hematite replacement. Brecciated samples show rock fragments with variably oxidized magnetite, where crystals range from anhedral to euhedral, often with marginal or extensive hematite replacement (Figure 9A,A-1-a). Some magnetite crystals are fully transformed to hematite, losing their original phase entirely (Figure 9A-2-a).

Paragenetic analysis indicates that chalcocite formed after magnetite. Magnetite crystals in the matrix are mostly anhedral and commonly partially replaced by hematite, with disseminated magnetite also found within malachite veinlets (Figure 9A-2b). Well-formed euhedral magnetite crystals occur in brecciated fragments (Figure 9A-2-a, A-2-b), while anhedral grains replace portions of the host rock. Magnetite crystals occasionally appear disseminated within non-metallic, chalcocite-bearing veinlets and show variable hematitization. These relationships suggest magnetite crystallized early, followed by the intrusion of chalcocite-bearing veinlets (Figure 9A).

Bornite is primarily observed as thick veinlets and massive zones. Secondary chalcocite and covellite replace bornite margins, indicating alteration of bornite to these copper sulfides (Figure 9C-2-a, C-2-b). Covellite–secondary chalcocite veinlets are sometimes present within large bornite grains along fractures. Fine-grained chalcopyrite occurs at bornite margins, often intergrown (Figure 9C-2-b).

The ore assemblage comprises chalcocite, minor chalcopyrite and bornite, covellite, pyrite, malachite, and chrysocolla (Figure 9H–K). Primary sulfides (chalcocite, bornite, pyrite) were subsequently altered to secondary minerals (malachite, covellite, iron hydroxides) during supergene enrichment (Figure 7, Figure 8 and Figure 9). Alteration is dominated by chloritic, carbonatic, and zeolitic assemblages, with limited argillic zones. Chalcocite, characterized by high copper and lack of iron, is the key primary copper mineral in mantle-type deposits. In contrast, chalcopyrite—common in porphyry and VMS deposits—contains less copper and more iron. The absence of iron-bearing gangue and quartz indicates mineralizing fluids were copper-rich, iron-poor, and silica-deficient, reflecting highly reducing conditions, further supported by the lack of epidote (Figures 8, 9).

Although volcanic host rocks contain up to 5% primary magnetite, no significant iron-rich mineralization (e.g., chalcopyrite, pyrite) or alteration is present, suggesting these rocks are unlikely metal sources. Similarly, intrusive sources would typically yield silica-, iron-, and aluminum-rich fluids; however, the absence of quartz associated with chalcocite argues against this.

Recent studies highlight the crucial role of framboidal pyrite and bitumen in the genesis of mantle-type copper deposits [46,47,48]. In Iranian stratabound deposits, framboidal pyrite commonly acts as a pre-ore reducing agent. For instance, in the Keshkuiyeh deposit, peak copper grades correlate with veins and zones enriched in pyrobitumen and framboidal pyrite. These pyrite textures, occurring as fine-grained aggregates near copper mineralization (Figure 8H, J, K), suggest a pivotal role in ore formation. The pre-ore assemblage in Iranian mantle-type deposits typically comprises framboidal pyrite alongside pervasive propylitic alteration of andesitic host rocks, interpreted as early diagenetic or pre-hydrothermal features [22]. Framboidal pyrite formed during early diagenesis may be partially replaced by Cu-(Fe) sulfides, such as chalcocite and bornite, during burial diagenesis. The presence of zeolite minerals during subsequent burial indicates ongoing diagenetic evolution [22]. These data imply that metal-bearing hydrothermal fluids interacted with pre-existing pyrite, which served as a reducing agent to induce copper sulfide precipitation [48,49]. Similar pyrite replacement processes have been reported for stratabound copper deposits in central Chile [50,51].

At Sepid-Sarve, chalcocite— the hallmark mineral of Manto-type deposits dominates, with lesser bornite, pyrite, digenite, magnetite, and hematite observed as primary sulfides. Secondary minerals, including malachite, azurite, secondary chalcocite, covellite, and hematite, occur throughout the deposit (Figure 10). These exhibits disseminated, cavity-filling (amygdaloidal), martitized, replacement, and skeletal textures, predominantly as veinlets of varying lengths and thicknesses, as well as disseminated grains.

4.2. Alteration

Magmatic and hydrothermal fluids related to volcanic activity caused generally weak alteration in the host rocks, mainly resulting in chloritic and iron oxide alterations (Figure 11A). Alteration in the Sepid-Sarve deposit, consistent with the broader northwest Bardeskan region, occurs in two stages: pre-mineralization and syn-mineralization.

Pre-mineralization alteration comprises carbonate (Figure 11B), silicification (Figure 11C), and propylitic assemblages characterized by epidote and chlorite (Figure 11D, E). Syn-mineralization alteration mainly includes carbonate (Figure 11E, F), zeolite (Figure 11F, G), calcite (Figure 11G, H), minor chlorite, and subordinate argillic alteration. The pre-mineralization stage corresponds to regional propylitic alteration, while syn-mineralization alteration is localized and marked by carbonate, calcite, zeolite, and minor argillic alteration (Figure 11A).

Chlorite forms via alteration of mafic minerals or hydrothermal introduction of Fe and Mg. Zeolitic alteration occurs locally and is linked to hydrothermal alteration, low-grade metamorphism, and meteoric water influence in intermediate to basic rocks. Iron oxide staining, chloritization, and iddingsitization are weak to moderate in andesitic units, with weak carbonate and clay alteration. Iron oxide alteration likely results from pyrite oxidation to goethite–hematite, facilitating sulfur release and sulfide mineral formation such as chalcocite and digenite. Secondary copper oxides like malachite derive from primary sulfides (Figure 11B).

Pervasive propylitic alteration of variable intensity affects the Eocene volcanic–sedimentary complex, including andesitic, basaltic, and microgabbro units, with more intense alteration in olivine basalt and microgabbro (Figure 11D, F). Alteration minerals mainly include chlorite, calcite, and epidote, with minor jarosite, adangsite, iron oxides, and magnetite appearing as veinlets, cavity fillings, replacements, and disseminations. Plagioclase and mafic minerals are locally replaced by up to 20% chlorite and calcite and up to 5% epidote (Figure 11J).

Carbonate alteration, with 30–35% intensity, occurs in andesitic units and Nummulitic limestone near mineralized zones, both at surface and drill cores. Calcite dominates, replacing plagioclase and mafic minerals and filling veinlets and cavities in andesite and limestone (Figure 11B, E, F). Plagioclase alteration reaches 20%, mafic minerals up to 90%. Vesicles in andesite are partially or fully filled with calcite, which commonly grows with colloform texture up to 3 cm in size and density of ~25 per m². Limestone mineralization involves calcium leaching and calcite veinlet precipitation.

Silicic alteration is mainly observed in andesitic tuff, olivine basalt, and andesite units, often coexisting with carbonate and propylitic alterations. Cryptocrystalline quartz is the principal mineral, occurring as cavity fillings and veinlets, sometimes bordered by calcite. Veinlet thickness ranges from 0.1 to 1 mm, with densities up to 10 per m² (Figure 11G).

Weak argillic alteration involves 10–15% feldspar replacement by clay minerals in some volcanic samples. Zeolitic alteration is prominent in megaporphyritic andesite near mineralization, mainly filling vesicles (Figure 11G, H, I).

4.3. Fluid Inclusions

In the Sepid-Sarve copper deposit, fine-crystalline calcite is the most abundant and significant mineral associated with the main mineralizing fluid. The studied fluid inclusions occur as Mono-phase (liquid, L; vapor, V) and two-phase liquid-rich plus vapor-rich (L+V) types (Figure 12B,F), mostly irregular, rectangular, spherical, or elliptical in shape, and dark in color. Mono-phase vapor inclusions range from 5 to 15 microns, predominantly irregular in shape. Two-phase liquid-rich inclusions constitute 5–40% of the inclusion volume, mostly irregular to elliptical, with sizes between 2 and 10 microns. Two-phase vapor-rich inclusions measure about 14 microns, with spherical to flattened shapes. Overall, fluid inclusion sizes in the deposit range from a few microns upward. Mono-phase vapor inclusions occur sporadically and individually within the host (Figure 12D,E), while liquid-rich two-phase inclusions are mostly isolated (Figure 12A,C,G,H).

Homogenization temperatures (Th) of vapor + liquid (V+L) inclusions range from 90 to 386 °C, averaging 193.75 °C, while for liquid + vapor (L+V) inclusions, Th ranges from 265 to 341 °C, averaging 302.6 °C (Figure 13A). The first melting temperature is directly related to the salt composition of the hydrothermal fluid [52] (Figure 13B). Salinity varies from 3.08 to 13.38 wt.% NaCl, with an average of 7.56 wt.% (Figure 13C).

Based on fluid inclusion studies, the formation temperature of the Sepid-Sarve copper deposit ranged between 90 and 386 °C, originating from a solution containing 3.08 to 13.38 wt.% NaCl. The final ice melting temperatures of fluid inclusions vary from –6.9 to –8.1 °C. Salinity, a key factor influencing fluid density, was evaluated using homogenization temperature versus salinity plots for mineralizing inclusions from Sepid-Sarve [53]. Fluid densities vary between 0.7 and 1 g/cm³ (Figure 14A,B). A decrease in homogenization temperature from 386 to 90 °C corresponds to a density increase from below 0.7 to 1 g/cm³, resulting in reduced fluid velocity and enhanced mineral precipitation.

The increase in salinity in fluid inclusions over the 90 to 386 °C range is accompanied by rising density. This process is linked to boiling phenomena whereby, due to cooling or pressure drop, fluids evolve toward a denser, more saline phase compared to the initial fluid. However, salinity increase with decreasing temperature cannot be solely explained by simple boiling. [54] argued that boiling-induced salinity increases are relatively small compared to the larger salinity variations observed here (Figure 14B). Meteoric waters with low temperature and salinity infiltrate faults and fractures, mixing with magmatic fluids and producing convective fluid mixing between deep intrusive bodies and shallower levels.

The homogenization temperature–salinity diagram reflects a progressive cooling and dilution of the mineralizing fluid, with both temperature and salinity declining during late-stage deposition (Figure 14C). Hydrothermal systems dominated by hydrostatic pressure commonly experience boiling, and vapor pressure can be used to estimate trapping depths. Such systems generally relate to near-surface environments within the upper ~3 km of the crust [55,56].

The coexistence of liquid-rich, vapor-rich, two-phase inclusions with varying fill volumes, as well as some homogenized vapor inclusions, confirms boiling processes in the study area. Using the [57] diagram, the depth of mineralization in Sepid-Sarve was estimated based on homogenization temperatures ranging from 90 to 386 °C. The most frequent homogenization temperature range (100–300 °C) corresponds to trapping depths between ~100 and 1550 meters (Figure 16B).

Figure 15. (A) Homogenization temperature versus salinity diagram for determining dominant metal complexes [58]. (B) Homogenization temperature versus salinity plot of various deposit types [56], with the Sepid-Sarve copper deposit samples falling within the epithermal deposit field. (C) Salinity versus homogenization temperature diagram of fluid inclusions showing different water-type fields [62].

Figure 15. (A) Homogenization temperature versus salinity diagram for determining dominant metal complexes [58]. (B) Homogenization temperature versus salinity plot of various deposit types [56], with the Sepid-Sarve copper deposit samples falling within the epithermal deposit field. (C) Salinity versus homogenization temperature diagram of fluid inclusions showing different water-type fields [62].

Figure 16. (A) Pressure vs. salinity diagram showing fluid inclusion pressures of 26/5–28/5 MPa in the Sepid-Sarve copper deposit [60]. (B) Homogenization temperature versus depth diagram for estimating trapping depth of fluids in the Sepid-Sarve copper deposit [57].

Figure 16. (A) Pressure vs. salinity diagram showing fluid inclusion pressures of 26/5–28/5 MPa in the Sepid-Sarve copper deposit [60]. (B) Homogenization temperature versus depth diagram for estimating trapping depth of fluids in the Sepid-Sarve copper deposit [57].

To identify the dominant metal-bearing complex, fluid inclusion data were plotted on the [58] diagram, showing that Sepid-Sarve fluids fall within the sulfide complex field (Figure 15A). This suggests sulfate (SO₄²⁻) or bisulfide (HS⁻) as the dominant anions, consistent with the presence of minerals such as chalcocite and bornite. Using homogenization temperatures and boiling point data, and applying [57], the trapping pressure of fluid inclusions was calculated via P = Hρg, yielding pressures between 240 and 260 bar (Figure 16A). The relationship between homogenization temperature and pressure is illustrated in Figure 16.

According to [56] classification of ore-forming systems based on salinity and homogenization temperature, the Sepid-Sarve samples fall within the epithermal deposit range (Figure 15B,C). Origin determination using the [59] diagram indicates that the fluids derive from magmatic waters, metamorphic waters, and meteoric waters, which likely reacted with host rocks during ascent, migrating along faults and fractures to deposit metals at the contact between conglomerate, andesite, agglomerate, and nummulitic limestone.

Table 1. Microthermometric data of fluid inclusions in the Sepid-Sarve copper deposit.

Sample N.

Mineralization Type

Mineral

phase

Tice

TICE+

salinity

Th V-L

Th L-V

Tutec

THH

THT

S4

Mlc+Az+Cal

Calcite

L+V

(-2)-(-5/5)

2/0-5/5

3/38-9/09

135-210

-36/3

135-210

S7

Mlc+Az+Cal

Calcite

L+V

(-2/7)-(9/2)

2/7-9/2

4/48-12/63

176-330

-43

176-330

S7

Mlc+Az+Cal

Calcite

V+L

(-7/8)-(-8/20)

7/8-8/2

11/19-11/6

265-341

265-341

S6

Cct+Cal

Calcite

L+V

(-1/8)-(-9/4)

1/8-9/4

3/08-13/17

90-386

-49/6

90-386

S9

Cct+Cal

Calcite

L+V

(-2)-(-5)

2/0-5/0

3/39-7/86

149-275

-37/2

149-275

S9

Cct+Cal

Calcite

V+L

(-5/7)-(-5/8)

5/7-5/8

8/76-8/98

280-300

280-300

S15

Cct+Bn+Ccp+Mag+Qtz

Quartz

L+V

(-2)-(-9/6)

2/0-9/6

3/39-13/38

100-234

-39

100-234

S15

Cct+Bn+Ccp+Mag+Qtz

Quartz

L+V+S

-8

8

11/55

280

280

Table 1. Microthermometric data of fluid inclusions in the Sepid-Sarve copper deposit.

Sample N.

Mineralization Type

Mineral

phase

Tice

TICE+

salinity

Th V-L

Th L-V

Tutec

THH

THT

S4

Mlc+Az+Cal

Calcite

L+V

(-2)-(-5/5)

2/0-5/5

3/38-9/09

135-210

-36/3

135-210

S7

Mlc+Az+Cal

Calcite

L+V

(-2/7)-(9/2)

2/7-9/2

4/48-12/63

176-330

-43

176-330

S7

Mlc+Az+Cal

Calcite

V+L

(-7/8)-(-8/20)

7/8-8/2

11/19-11/6

265-341

265-341

S6

Cct+Cal

Calcite

L+V

(-1/8)-(-9/4)

1/8-9/4

3/08-13/17

90-386

-49/6

90-386

S9

Cct+Cal

Calcite

L+V

(-2)-(-5)

2/0-5/0

3/39-7/86

149-275

-37/2

149-275

S9

Cct+Cal

Calcite

V+L

(-5/7)-(-5/8)

5/7-5/8

8/76-8/98

280-300

280-300

S15

Cct+Bn+Ccp+Mag+Qtz

Quartz

L+V

(-2)-(-9/6)

2/0-9/6

3/39-13/38

100-234

-39

100-234

S15

Cct+Bn+Ccp+Mag+Qtz

Quartz

L+V+S

-8

8

11/55

280

280

4.4. Stable Isotopes

4.4.1. Sulfur Isotopes

To investigate the origin and nature of ore-forming and alteration-related fluids in the Sepid-Sarve copper deposit, sulfur and oxygen isotope analyses were conducted. One sulfide mineral sample and two silica mineral samples were selected. Sample preparation was performed at a laboratory in Tehran, after which the specimens were sent to the Stable Isotope Research Laboratory at Arak University for isotopic analysis (Table 2). Sulfur isotope fractionation occurs via two principal mechanisms that result in natural variations of δ³⁴S values in different environments and materials [63]: (a) Abiotic isotopic fractionation, which takes place through various chemical exchange reactions between sulfate and sulfide species, or among different types of sulfides;(b) Biological (kinetic) isotopic fractionation, which occurs through microbial metabolic processes, leading to the preferential incorporation of lighter sulfur isotopes (³²S) into biologically-derived products. In hydrothermal systems, sulfur is primarily present in reduced forms such as H₂S, HS⁻, and S²⁻, which play key roles in ore deposition processes [64]. The most stable sulfur species in reducing and alkaline environments are H₂S and HS⁻. Under these conditions, the oxygen fugacity (fO₂) and pH of the hydrothermal fluids tend to be low, and the δ³⁴S values of deposited sulfide minerals generally reflect those of the total sulfur in the fluid [65].

The The sulfide minerals in the Sepid-Sarve copper deposit exhibit a relative depletion in the heavy sulfur isotope (³⁴S). One plausible explanation for this isotopic depletion is that hydrothermal sulfur was derived from the leaching of sulfide minerals in reduced sedimentary environments [63]. The δ³⁴S values of H₂S in the ore-forming hydrothermal fluid are estimated to range from –22.7‰ to –2.7‰. For comparison, δ³⁴S values in peridotitic kimberlites of mantle origin are approximately 0 ± 1‰, and in basaltic rocks typically range between 0‰ and 5‰ [63]. However, recycling of crustal materials (e.g., sulfides and marine sulfates) in subduction zones, as well as mantle contamination, can cause significant heterogeneity in mantle sulfur isotopic values, ranging from –10‰ to +10‰ [66]. Therefore, the isotopic composition of the hydrothermal fluid may suggest a magmatic origin.

Sulfur in the Sepid-Sarve deposit may have been sourced either from the leaching of primary volcanic-hosted sulfide minerals or directly from magmatic-hydrothermal fluids. One contributing factor to the ³⁴S depletion could be the preferential leaching of ³²S (due to its higher reactivity compared to ³⁴S; [63]) from volcanic and pyroclastic host rocks by meteoric-dominated hydrothermal fluids. Another factor may be isotopic fractionation between sulfate and sulfide species in the hydrothermal fluid. Following isotopic equilibrium fractionation and subsequent ³⁴S depletion in sulfides, two potential scenarios could have occurred:

(a) During sulfide mineralization, the hydrothermal fluid expelled and dispersed the sulfate phase away from the vein system. The scarcity of sulfate minerals alongside sulfides may support this removal of sulfate-bearing fluids during or after mineralization. This scenario is consistent with acidic and oxidizing conditions that facilitate sulfide precipitation from sulfur complexes and promote sulfate phase separation.

(b) Prior to the current stage of sulfide deposition, the sulfate phase may have precipitated at greater depths, allowing the upward migration of a sulfide-dominant hydrothermal fluid into the vein system, where it precipitated ³⁴S-depleted sulfide minerals.

Another possible mechanism responsible for the observed ³⁴S depletion is the derivation of sulfur from sedimentary sequences containing biogenic sulfides, which were leached by meteoric-dominated hydrothermal fluids. However, considering the fluid inclusion temperatures and the absence of sulfur-bearing sedimentary host rocks, a biogenic sulfur source and bacterial reduction processes are unlikely. Instead, the mixing of meteoric and magmatic waters likely produced a system in which both sulfate and reduced sulfur were generated.

The δ³⁴S values of sulfide minerals in the Sepid-Sarve deposit, when compared to various geological environments and sulfur sources, closely resemble those from ilmenite-series granites and hydrothermal sulfides (Figure 17A,C). Additionally, Figure 18 compares the δ³⁴S values of Sepid-Sarve with those from other deposits in Central Iran, the Urumieh–Dokhtar belt, and other global analogues. These comparisons suggest that sulfide mineralization in the Sepid-Sarve deposit was formed from fluids that derived sulfur predominantly via the reduction of seawater sulfate.

4.4.2. Oxygen Isotopes

The oxygen isotope compositions (δ¹⁸O) of hydrothermal fluids associated with altered and mineralized carbonate samples from the Sepid-Sarve copper deposit are presented in Table 3. The δ¹⁸O values of mineralized samples fall within or close to the compositional fields of meteoric water, basalts, and sedimentary rocks, suggesting that the hydrothermal fluids involved in ore deposition were primarily of magmatic origin but underwent mixing with meteoric waters. The δ¹⁸O values of altered but barren carbonate samples plot entirely within the meteoric water range, slightly shifting toward magmatic compositions, indicating that carbonate alteration may have resulted from interaction between magmatic and meteoric fluids.

This isotopic transition is consistent with the fluid inclusion microthermometry data, which indicate decreasing fluid temperatures due to fluid mixing and dilution. The evolution from a magmatic-dominated fluid to cooler, oxidizing, and more diluted meteoric waters is characteristic of hydrothermal systems where ore precipitation is triggered by fluid mixing, cooling, and changes in redox conditions (Figure 19).

4.5. Geochemistry of Major and Trace Elements

Twenty samples from volcanic outcrops and drill cores were selected for petrographic and geochemical analyses of major and trace elements, and 75 samples were analyzed specifically for trace and minor elements (Table 4). For rock classification and to minimize alteration effects, immobile elements were used. Based on the Nb/Y vs. Zr/Ti diagram [96], the volcanic rocks are mainly classified as andesite, basaltic andesite, basalt, and alkali basalt (Figure 20A). Using the Nb/Y vs. Zr/TiO₂ plot [96], the samples are categorized as sub-alkaline basalt, basaltic andesite, and andesite (Figure 20B). The Zr/TiO₂ vs. SiO₂ plot indicates compositions ranging from andesite to dacite and rhyodacite (Figure 20C). Based on Co vs. Th [97], the samples fall within the calc-alkaline basalt–andesite–andesite series (Figure 20D). To distinguish between continental arc and island arc settings, Ta/Yb vs. Th/Yb (Figure 20E) and Zr vs. Zr/Y (Figure 20F) diagrams [98] were used. The samples plot within the continental arc field, indicating a source in an enriched mantle influenced by subduction-zone fluids. High Th contents likely result from crustal contamination and fluids derived from the subducted slab. Magmas from active continental margins typically show stronger enrichment in incompatible LREEs compared to island arc magmas, likely due to a combination of enriched mantle sources and crustal assimilation. The Zr vs. Zr/Y plot further confirms their placement within the continental volcanic arc field.

The SiO₂ content in the volcanic rocks ranges from 57.85 to 64.1 wt%, with an average of 61.95 wt%. Al₂O₃ varies between 16.65 and 18.70 wt%, CaO from 2.18 to 6.82 wt%, Na₂O from 3.76 to 5.97 wt%, K₂O from 3.17 to 4.97 wt%, MgO from 0.01 to 1.66 wt%, TiO₂ from 0.55 to 0.74 wt%, and total Fe₂O₃ (Fe₂O₃T) from 3.18 to 7.07 wt%. To understand the processes influencing magma evolution, such as fractional crystallization, binary plots of oxide–oxide, oxide–element, and element–element (Harker diagrams; [99]) were utilized. Generally, linear or negatively correlated trends suggest the influence of processes such as fractional crystallization, magma mixing, or crustal assimilation. However, significant scatter may reflect crystal accumulation, alteration, or crustal contamination [100]. Among the major oxides, Na₂O and Al₂O₃ show increasing trends (Figure 21C,E). K₂O versus SiO₂ lacks a clear trend and shows scattered values (Figure 21F). In calc-alkaline suites, the concentrations of MgO, CaO, Fe₂O₃, and TiO₂ decrease with increasing SiO₂, which aligns with the volcanic rocks of the Sepid-Sarve deposit, reflecting the crystallization of ferromagnesian minerals as well as plagioclase, amphibole, and clinopyroxene fractionation [101] (Figure 21A,B). The decreasing CaO trend corresponds to plagioclase crystallization from anorthite to albite (Figure 21D). The increasing trends in Na₂O and possibly K₂O may be attributed to the presence of alkali feldspar in the final stages of crystallization or potassium addition during hydrothermal alteration. These trends overall highlight the role of continental crust in the magmatic evolution.

In the SiO₂ vs. Sr diagram (Figure 21G), Sr shows a decreasing trend consistent with plagioclase fractionation at low pressure. Since Sr is compatible in plagioclase but not in clinopyroxene, this trend also corresponds to the decreasing CaO pattern due to plagioclase evolution. Europium (Eu) first increases with increasing SiO₂ and then decreases, reflecting plagioclase involvement in magma evolution, further supported by a positive correlation between Sr and Eu (Figure 21H). The Dy/Yb ratio shows a slight decrease with increasing SiO₂, indicating hornblende or clinopyroxene fractionation (Figure 21I). Sr content also decreases with MgO, suggesting plagioclase removal (Figure 21J). The Rb/Sr ratio first decreases and then increases with increasing Sr, displaying an almost linear trend consistent with plagioclase, hornblende, and clinopyroxene fractionation [102] (Figure 21L).

The Th/Ti vs. Th (Figure 21K), La vs. La/Yb (Figure 21M), and Th vs. Ti/Yb (Figure 21N) diagrams [103] indicate crustal contamination during magmatic differentiation in the Sepid-Sarve region.

Spider diagrams normalized to primitive mantle [104] show negative slopes with enrichment of LILEs over HFSEs, and depletion in HREEs — features characteristic of calc-alkaline magmas. The observed positive Pb anomaly is indicative of crustal contamination during magma evolution (Figure 22A). In the chondrite-normalized REE diagram, volcanic rocks from Sepid-Sarve exhibit moderate LREE enrichment, HREE depletion, and weak negative Eu anomalies, further reflecting plagioclase fractionation and a continental arc signature (Figure 22B).

Table 4. Chemical Analysis Results of the Studied Samples from the Sepid-Sarve Copper Deposit.

Sample

SiO2

Al2O3

Na2O

MgO

K2O

TiO2

CaO

Fe2O3

Ag

Pb

Cu

Mo

Bi

Fe

Co

Ni

Ti

Sb

CH.M-12A

64/05

18

3/9

0/46

3/66

0/58

2/95

5/17

0/6

15

1068

2/2

0/05

35794

8/2

29

4456

1/1

CH.M-13A

63/9

17/69

3/76

0/26

3/17

0/63

2/68

7/07

31/9

23

>5%

3/8

0/05

35658

4/2

18

3139

2/2

CH.M-15A

57/85

16/75

3/84

0/01

4/12

0/74

2/31

5/38

>10

38

>1000

6/82

0/46

41000

4/9

13/8

3800

1/33

CH.M-16A

61/25

17/16

5/38

1/46

3/67

0/63

5/71

4/2

0/42

22

541

3/01

0/16

37900

8/4

19/6

4000

0/58

CH.M-17A

60/75

16/85

4/99

0/05

4/66

0/63

6/81

4/92

0/43

20

575

3/75

0/12

34700

6/2

16/2

3900

0/84

CH.M-17B

59/64

18/25

5/15

1/66

3/75

0/66

6/01

4/69

0/44

22

458

3/25

0/14

32000

4/2

17/3

3400

0/84

CH.M-18A

61/35

18/7

5/97

0/05

4/97

0/71

5/23

3/18

0/5

14

1068

2/1

0/05

3421

2/3

29

3215

1/5

CH.M-18B

63/85

18/63

4/75

0/35

4/08

0/55

2/68

6/07

0/14

17

397

2/79

0/11

32200

9/9

18/7

0/44

1/84

Se.Sa-01

63/3

18/51

4/65

0/43

4/01

0/61

2/64

5/64

0/1

2

54

2

0/1

38519

20

31

3776

<0.5

Se.Sa-02

61/9

17/36

5/05

0/27

3/56

0/58

4/85

4/15

0/1

4

40

0/1

0/1

47492

23/5

37

4245

<0.5

Se.Sa-03

60/85

16/65

4/86

0/05

4/61

0/6

6/45

4/73

0/1

4

63

0/1

0/1

48663

24/6

38

3605

1

Se.Sa-04

60/55

16/85

4/91

0/06

4/15

0/62

6/15

4/68

0/1

9

6186

0/1

0/1

37322

10

46

2650

<0.5

Se.Sa-05

59/85

18/15

4/94

1/45

3/74

0/65

6/82

4/72

0/7

32

284

0/1

0/1

45855

25/4

52

3746

<0.5

Se.Sa-06

62/75

16/95

4/7

0/85

3/64

0/66

2/95

4/33

0/1

16

153

0/1

0/1

45710

26

54

3712

<0.5

Se.Sa-07

63/9

17/35

3/84

0/46

4/01

0/6

2/63

6/1

0/1

3

40

0/1

0/1

27859

5/1

5

2572

<0.5

Se.Sa-08

64/1

18/24

3/95

0/48

3/59

0/58

2/24

5/08

0/1

4

138

0/1

0/1

54893

40/6

299

3802

<0.5

Se.Sa-09

64

18/52

4/15

0/49

3/68

0/57

2/18

5/04

0/1

7

29

1

0/1

54821

39/5

282

3446

<0.5

Se.Sa-10

62/8

17/15

4/86

0/64

3/62

0/63

2/9

4/3

0/1

6

19

0/1

0/1

53398

41/8

329

3458

<0.5

Se.Sa-11

61/65

17/25

5/03

1/15

3/61

0/61

4/86

4/72

0/1

5

25

0/1

0/1

48511

44/4

467

2298

<0.5

Se.Sa-12

60/75

17/45

5/14

1/04

3/74

0/59

4/74

4/66

0/1

3

16

0/1

0/1

52377

49/6

506

2671

<0.5

Table 4. Chemical Analysis Results of the Studied Samples from the Sepid-Sarve Copper Deposit.

Sample

SiO2

Al2O3

Na2O

MgO

K2O

TiO2

CaO

Fe2O3

Ag

Pb

Cu

Mo

Bi

Fe

Co

Ni

Ti

Sb

CH.M-12A

64/05

18

3/9

0/46

3/66

0/58

2/95

5/17

0/6

15

1068

2/2

0/05

35794

8/2

29

4456

1/1

CH.M-13A

63/9

17/69

3/76

0/26

3/17

0/63

2/68

7/07

31/9

23

>5%

3/8

0/05

35658

4/2

18

3139

2/2

CH.M-15A

57/85

16/75

3/84

0/01

4/12

0/74

2/31

5/38

>10

38

>1000

6/82

0/46

41000

4/9

13/8

3800

1/33

CH.M-16A

61/25

17/16

5/38

1/46

3/67

0/63

5/71

4/2

0/42

22

541

3/01

0/16

37900

8/4

19/6

4000

0/58

CH.M-17A

60/75

16/85

4/99

0/05

4/66

0/63

6/81

4/92

0/43

20

575

3/75

0/12

34700

6/2

16/2

3900

0/84

CH.M-17B

59/64

18/25

5/15

1/66

3/75

0/66

6/01

4/69

0/44

22

458

3/25

0/14

32000

4/2

17/3

3400

0/84

CH.M-18A

61/35

18/7

5/97

0/05

4/97

0/71

5/23

3/18

0/5

14

1068

2/1

0/05

3421

2/3

29

3215

1/5

CH.M-18B

63/85

18/63

4/75

0/35

4/08

0/55

2/68

6/07

0/14

17

397

2/79

0/11

32200

9/9

18/7

0/44

1/84

Se.Sa-01

63/3

18/51

4/65

0/43

4/01

0/61

2/64

5/64

0/1

2

54

2

0/1

38519

20

31

3776

<0.5

Se.Sa-02

61/9

17/36

5/05

0/27

3/56

0/58

4/85

4/15

0/1

4

40

0/1

0/1

47492

23/5

37

4245

<0.5

Se.Sa-03

60/85

16/65

4/86

0/05

4/61

0/6

6/45

4/73

0/1

4

63

0/1

0/1

48663

24/6

38

3605

1

Se.Sa-04

60/55

16/85

4/91

0/06

4/15

0/62

6/15

4/68

0/1

9

6186

0/1

0/1

37322

10

46

2650

<0.5

Se.Sa-05

59/85

18/15

4/94

1/45

3/74

0/65

6/82

4/72

0/7

32

284

0/1

0/1

45855

25/4

52

3746

<0.5

Se.Sa-06

62/75

16/95

4/7

0/85

3/64

0/66

2/95

4/33

0/1

16

153

0/1

0/1

45710

26

54

3712

<0.5

Se.Sa-07

63/9

17/35

3/84

0/46

4/01

0/6

2/63

6/1

0/1

3

40

0/1

0/1

27859

5/1

5

2572

<0.5

Se.Sa-08

64/1

18/24

3/95

0/48

3/59

0/58

2/24

5/08

0/1

4

138

0/1

0/1

54893

40/6

299

3802

<0.5

Se.Sa-09

64

18/52

4/15

0/49

3/68

0/57

2/18

5/04

0/1

7

29

1

0/1

54821

39/5

282

3446

<0.5

Se.Sa-10

62/8

17/15

4/86

0/64

3/62

0/63

2/9

4/3

0/1

6

19

0/1

0/1

53398

41/8

329

3458

<0.5

Se.Sa-11

61/65

17/25

5/03

1/15

3/61

0/61

4/86

4/72

0/1

5

25

0/1

0/1

48511

44/4

467

2298

<0.5

Se.Sa-12

60/75

17/45

5/14

1/04

3/74

0/59

4/74

4/66

0/1

3

16

0/1

0/1

52377

49/6

506

2671

<0.5

Table 4. (continued)

Tl	Li	Be	Na	K	Ca	Mg	Al	Ga	Cs	Rb	Ba	Sr	Th	U	Y	Zr	Hf	Nb	Ta
0/55	11	4/9	31348	34626	24966	2871	7465	15/2	4/8	94	613	316	4/35	1/2	14/4	119	3/08	5/5	0/4
0/41	10	3/9	27389	28248	18024	1856	6497	15/2	3/9	79	420	315	4/29	1/5	14/5	18	3/12	4/6	0/31
0/6	15/8	1/1	32600	35200	13500	2400	2/66	15/2	5/39	74/3	254	132	2/7	0/3	8/7	109	3/78	4/7	0/46
0/51	22/6	1/4	36600	36200	47700	3900	4/17	18/6	5/67	78/2	368	223	3/9	0/3	18	128	3/65	5/4	0/43
0/58	18/5	1/2	35700	38100	34500	1400	2/56	17/3	6/69	101	400	184	3/2	0/2	15	110	3/15	4/9	0/37
0/56	17/3	1/5	35723	37512	35000	1290	2/36	16/2	6/78	122	359	190	3/2	0/2	15	110	3/14	5/3	0/25
0/39	10	3/9	30221	36620	23230	2781	7365	15	4	802	512	312	3/9	1/3	14	100	3/04	5/2	0/26
0/56	20/6	1/3	29000	42100	32220	4100	3450	17/8	5/51	86/1	398	173	4/5	0/2	20/6	142	4/15	6	0/48
<0.1	14	0/2	22946	7106	47451	>2%	75884	15/1	0/7	13	126	197	0/91	0/3	15/5	48	1/75	8/1	0/83
<0.1	12	0/2	23126	6382	53087	>2%	81673	15/4	1	10	105	272	0/84	0/4	15/5	49	1/81	4/6	0/57
0/16	3	0/2	19532	2373	55661	>2%	82422	15/8	1/7	6	61	141	0/28	0/2	17/2	42	1/64	4/9	0/6
<0.1	<1	0/2	1761	244	>10%	7241	89852	15/7	<0/5	<1	14	43/7	0/68	0/9	9/3	28	1/29	6/4	0/58
0/15	6	0/2	25505	3146	48663	>2%	86485	16	0/9	3	60	196	<0.1	0/2	11/8	33	1/32	5/4	0/53
<0.1	7	0/2	24429	3515	54628	>2%	88529	16/1	0/9	4	61	182	<0.1	0/2	12/9	33	1/36	4/5	0/5
0/1	18	0/2	26564	16425	17194	9009	70541	15/9	1/8	42	92	169	0/33	0/4	35/7	62	2/16	4/1	0/39
<0.1	17	0/2	15347	8660	66205	>2%	68084	15/6	2/5	17	110	330	0/71	0/3	13/1	37	1/54	5/6	0/62
<0.1	32	0/2	12006	7428	48237	>2%	62404	16/3	1/1	9	105	264	0/58	0/34	11/6	34	1/34	5/2	0/66
<0.1	33	0/2	12023	7004	64240	>2%	59380	16/7	1/4	13	89	274	0/47	0/2	11/3	33	1/33	4/9	0/49
<0.1	39	0/2	2871	3508	42747	>2%	34053	16/6	0/7	<1	55	121	0/2	0/3	6/9	21	1/02	3/7	0/34
0/1	40	0/2	4524	1614	40212	>2%	42514	16/3	4/7	2	20	89/4	0/39	0/3	8/7	25	1/09	3/8	0/36

Table 4. (continued)

Tl	Li	Be	Na	K	Ca	Mg	Al	Ga	Cs	Rb	Ba	Sr	Th	U	Y	Zr	Hf	Nb	Ta
0/55	11	4/9	31348	34626	24966	2871	7465	15/2	4/8	94	613	316	4/35	1/2	14/4	119	3/08	5/5	0/4
0/41	10	3/9	27389	28248	18024	1856	6497	15/2	3/9	79	420	315	4/29	1/5	14/5	18	3/12	4/6	0/31
0/6	15/8	1/1	32600	35200	13500	2400	2/66	15/2	5/39	74/3	254	132	2/7	0/3	8/7	109	3/78	4/7	0/46
0/51	22/6	1/4	36600	36200	47700	3900	4/17	18/6	5/67	78/2	368	223	3/9	0/3	18	128	3/65	5/4	0/43
0/58	18/5	1/2	35700	38100	34500	1400	2/56	17/3	6/69	101	400	184	3/2	0/2	15	110	3/15	4/9	0/37
0/56	17/3	1/5	35723	37512	35000	1290	2/36	16/2	6/78	122	359	190	3/2	0/2	15	110	3/14	5/3	0/25
0/39	10	3/9	30221	36620	23230	2781	7365	15	4	802	512	312	3/9	1/3	14	100	3/04	5/2	0/26
0/56	20/6	1/3	29000	42100	32220	4100	3450	17/8	5/51	86/1	398	173	4/5	0/2	20/6	142	4/15	6	0/48
<0.1	14	0/2	22946	7106	47451	>2%	75884	15/1	0/7	13	126	197	0/91	0/3	15/5	48	1/75	8/1	0/83
<0.1	12	0/2	23126	6382	53087	>2%	81673	15/4	1	10	105	272	0/84	0/4	15/5	49	1/81	4/6	0/57
0/16	3	0/2	19532	2373	55661	>2%	82422	15/8	1/7	6	61	141	0/28	0/2	17/2	42	1/64	4/9	0/6
<0.1	<1	0/2	1761	244	>10%	7241	89852	15/7	<0/5	<1	14	43/7	0/68	0/9	9/3	28	1/29	6/4	0/58
0/15	6	0/2	25505	3146	48663	>2%	86485	16	0/9	3	60	196	<0.1	0/2	11/8	33	1/32	5/4	0/53
<0.1	7	0/2	24429	3515	54628	>2%	88529	16/1	0/9	4	61	182	<0.1	0/2	12/9	33	1/36	4/5	0/5
0/1	18	0/2	26564	16425	17194	9009	70541	15/9	1/8	42	92	169	0/33	0/4	35/7	62	2/16	4/1	0/39
<0.1	17	0/2	15347	8660	66205	>2%	68084	15/6	2/5	17	110	330	0/71	0/3	13/1	37	1/54	5/6	0/62
<0.1	32	0/2	12006	7428	48237	>2%	62404	16/3	1/1	9	105	264	0/58	0/34	11/6	34	1/34	5/2	0/66
<0.1	33	0/2	12023	7004	64240	>2%	59380	16/7	1/4	13	89	274	0/47	0/2	11/3	33	1/33	4/9	0/49
<0.1	39	0/2	2871	3508	42747	>2%	34053	16/6	0/7	<1	55	121	0/2	0/3	6/9	21	1/02	3/7	0/34
0/1	40	0/2	4524	1614	40212	>2%	42514	16/3	4/7	2	20	89/4	0/39	0/3	8/7	25	1/09	3/8	0/36

Table 4. (continued).

P

La

Ce

Pr

Nd

Sm

Eu

Er

Gd

Tb

Dy

Ho

Tm

Yb

Lu

1090

13

30

3/93

13/7

2/98

1

1/66

3/26

0/47

3/07

0/39

0/27

1/8

0/29

1237

11

26

3/7

12/6

2/89

0/95

1/61

3/26

0/46

3/05

0/39

0/27

1/7

0/26

1000

4/6

9/91

1/35

5/8

1/5

0/6

1/3

0

0/27

1/85

0/39

0/21

1/4

0/23

1000

8/9

19/5

2/6

11

3/02

1/04

2/34

0

0/5

3/37

0/68

0/33

2/3

0/36

800

5/6

13/5

1/82

8/02

2/35

0/92

1/77

0

0/42

2/78

0/57

0/3

1/8

0/3

900

6/1

12/5

1/73

9/26

2/36

0/77

1/8

0

0/32

2/45

0/58

0/21

2

0/23

1032

12

32

2/9

12/8

2/7

1

1/56

0

0/41

2/8

0/34

0/24

1/6

0/29

1000

7/4

18

2/31

10/8

3/24

1/02

2/5

3/39

0/52

3/3

0/74

0/32

2/5

0/41

407

6

7

1

6/9

0/22

0/67

1/34

2/16

0/39

2/79

0/39

0/3

1/9

0/31

523

7

7

1/33

7/4

0/21

0/66

1/29

2/23

0/39

3/07

0/39

0/29

2/1

0/3

278

4

4

0/66

4/7

0/02

0/43

1/57

2/1

0/42

3/33

0/65

0/31

2/1

0/32

397

8

7

1/1

5

0/02

0/31

0/81

1/97

0/31

2/1

0/58

0/2

1/4

0/18

195

3

1

0/28

2/5

0/02

0/35

1/19

1/79

0/32

2/41

0/57

0/22

1/7

0/21

205

3

2

0/28

2/4

0/02

0/33

1/13

1/78

0/32

2/46

0/58

0/21

1/8

0/17

590

6

8

1/46

7/8

0/48

0/8

2/86

2/93

0/67

5/06

0/64

0/54

3/3

0/53

656

7

8

1/43

7/2

0/26

0/7

1

2/36

0/39

2/88

0/43

0/2

1/8

0/23

619

7

8

1/44

7/4

0/12

0/66

1/13

2/13

0/36

2/51

0/42

0/23

1/6

0/17

594

7

7

1/43

7/9

0/02

0/6

1

2/29

0/36

2/5

0/71

0/17

1/6

0/21

373

4

3

0/66

3/7

0/02

0/36

0/67

1/76

0/24

1/97

0/67

0/1

1

0/15

432

5

5

0/86

4/3

0/02

0/23

0/76

1/92

0/32

2/18

0/65

0/12

1/3

0/24

Table 4. (continued).

P

La

Ce

Pr

Nd

Sm

Eu

Er

Gd

Tb

Dy

Ho

Tm

Yb

Lu

1090

13

30

3/93

13/7

2/98

1

1/66

3/26

0/47

3/07

0/39

0/27

1/8

0/29

1237

11

26

3/7

12/6

2/89

0/95

1/61

3/26

0/46

3/05

0/39

0/27

1/7

0/26

1000

4/6

9/91

1/35

5/8

1/5

0/6

1/3

0

0/27

1/85

0/39

0/21

1/4

0/23

1000

8/9

19/5

2/6

11

3/02

1/04

2/34

0

0/5

3/37

0/68

0/33

2/3

0/36

800

5/6

13/5

1/82

8/02

2/35

0/92

1/77

0

0/42

2/78

0/57

0/3

1/8

0/3

900

6/1

12/5

1/73

9/26

2/36

0/77

1/8

0

0/32

2/45

0/58

0/21

2

0/23

1032

12

32

2/9

12/8

2/7

1

1/56

0

0/41

2/8

0/34

0/24

1/6

0/29

1000

7/4

18

2/31

10/8

3/24

1/02

2/5

3/39

0/52

3/3

0/74

0/32

2/5

0/41

407

6

7

1

6/9

0/22

0/67

1/34

2/16

0/39

2/79

0/39

0/3

1/9

0/31

523

7

7

1/33

7/4

0/21

0/66

1/29

2/23

0/39

3/07

0/39

0/29

2/1

0/3

278

4

4

0/66

4/7

0/02

0/43

1/57

2/1

0/42

3/33

0/65

0/31

2/1

0/32

397

8

7

1/1

5

0/02

0/31

0/81

1/97

0/31

2/1

0/58

0/2

1/4

0/18

195

3

1

0/28

2/5

0/02

0/35

1/19

1/79

0/32

2/41

0/57

0/22

1/7

0/21

205

3

2

0/28

2/4

0/02

0/33

1/13

1/78

0/32

2/46

0/58

0/21

1/8

0/17

590

6

8

1/46

7/8

0/48

0/8

2/86

2/93

0/67

5/06

0/64

0/54

3/3

0/53

656

7

8

1/43

7/2

0/26

0/7

1

2/36

0/39

2/88

0/43

0/2

1/8

0/23

619

7

8

1/44

7/4

0/12

0/66

1/13

2/13

0/36

2/51

0/42

0/23

1/6

0/17

594

7

7

1/43

7/9

0/02

0/6

1

2/29

0/36

2/5

0/71

0/17

1/6

0/21

373

4

3

0/66

3/7

0/02

0/36

0/67

1/76

0/24

1/97

0/67

0/1

1

0/15

432

5

5

0/86

4/3

0/02

0/23

0/76

1/92

0/32

2/18

0/65

0/12

1/3

0/24

4.6. Tectono-Magmatic Model

Most Most of the known Manto-type copper deposits in Iran formed during the Eocene, a period characterized by intense tectonic and magmatic activity, especially within the Urumieh-Dokhtar Magmatic Arc (UDMA), Alborz Magmatic Arc (AMA), and the Sabzevar region. The Eocene is considered a key metallogenic epoch due to the closure of the Neo-Tethys Ocean. Many of Iran’s mineral deposits, including Manto-type copper deposits, are closely related to the evolution of the Neo-Tethys Ocean. These deposits typically form in arc-related or subduction zone settings and are generally associated with volcano-sedimentary environments along plate boundaries. They are often hosted by basaltic to andesitic rocks.

The general formation process of Manto-type copper deposits involves the generation of calc-alkaline magma in the subduction zone, heat-induced alteration caused by intrusive bodies, and the precipitation of copper along fractures and rock porosity through hydrothermal fluids derived from the mantle and lower crust. The Khaf–Doruneh Belt (in the southern Sabzevar Zone, northwest of Bardeskan) hosts Manto-type copper mineralization due to its tectonically active setting related to the subduction of Neo-Tethyan oceanic crust beneath the Central Iranian continental crust and the concurrent Eocene magmatic activity (Figure 23A). These deposits are typically stratabound and hosted in volcanic and sedimentary sequences.

The Sabzevar Ocean, situated between the Lut Block to the south and the Turan Block to the north, remained open since at least the mid-Cretaceous. Intra-oceanic subduction in the Sabzevar Basin likely began before the Albian, as evidenced by the presence of high-pressure metamorphic rocks [28]. Similar to UDMA magmatism, the Eocene–Oligocene magmatic rocks in the Kashmar region of the Sabzevar area formed during Eocene extensional tectonics [109]. The Eocene volcanic activity in Davarzan–Abbasabad (Sabzevar) began with continental to shallow marine Paleogene–Eocene sediments, followed by widespread subaqueous to subaerial basaltic and andesitic volcanism, volcaniclastic rocks, and sedimentary sequences during the middle to late Eocene [110]. These E–W-trending magmatic rocks were generated in response to the subduction of Neo-Tethyan oceanic lithosphere beneath the eastern Alborz (Figure 23B).

The mid-Eocene volcanic rocks of the Davarzan–Abbasabad area formed in a back-arc extensional setting and are associated with post-Eocene subduction-related magmatism (Figure 23B). This region experienced collisional and post-collisional tectonic events during the Oligocene to Miocene [110,111,112,113,114]. In western Sabzevar (Davarzan–Abbasabad area), these Eocene volcanic rocks host Manto-type Cu–Ag mineralization, including the Bozorg Mine, Damanjala, Asiadeyu, Labkalu, North Choghandarsar, South Choghandarsar, Elhak, and Dochileh deposits.

This model is based on the subduction of the Neo-Tethyan oceanic crust during the Eocene, which led to the development of calc-alkaline magmatism (andesite, basalt, trachyandesite) and hydrothermal activity within the Khaf–Doruneh Belt. Copper mineralization occurs at the contact between volcanic rocks (e.g., andesite) and fossiliferous limestones (mostly conglomerates), where copper-rich hydrothermal fluids precipitated in fractures and cavities. The subduction of the oceanic plate beneath the continental crust during the Neo-Tethys era triggered partial melting in the mantle wedge, producing calc-alkaline magma. The intrusion of this magma and the release of metal-bearing hydrothermal fluids caused mineralization in fractures and stress zones. This magmatic activity was accompanied by regional propylitic and argillic alteration (Figure 23C).

The primary ore minerals include chalcocite, bornite, and pyrite, while the supergene secondary minerals include malachite, azurite, and hematite. Regional zeolitic–carbonatic and local chloritic and hematitic alteration are present. Observed mineral textures include veinlets, cavity fillings, replacement, and stratabound forms. Calc-alkaline volcanic rocks, shear zones, deep faults, dykes, and intense alteration are characteristic features of this model (Figure 23D). The formation model of the Sepid-Sarve copper deposit closely aligns with this model, and the Kal Abri deposit, located east of the Sepid-Sarve copper deposit, represents another similar example.

5. Conclusion

The Sepid-Sarve copper deposit is a strata-bound (Manto-type) sulfide system hosted in volcanic-sedimentary sequences formed in a subduction-related magmatic arc setting. The volcanic succession, predominantly comprising basalt, basaltic andesite, and andesite, developed in a non-marine environment, followed by deposition of sedimentary conglomerates and limestone. Mineralization is mainly associated with chlorite, zeolite, and calcite alteration zones. Copper mineralization occurs predominantly as fine-grained chalcocite within permeable conglomerates and, to a lesser extent, in limestones. The ore-forming fluids, characterized by low silica content and reducing conditions, ascended along faults, depositing copper in fractures and pore spaces.

Geochemical and isotopic data suggest that the mineralizing fluids originated primarily from a magmatic source but underwent mixing with meteoric waters, promoting sulfide precipitation. Fluid inclusions indicate moderate salinities (3.08–13.38 wt.% NaCl eq.) and highlight boiling and fluid mixing as key mineralization mechanisms. Sulfur isotope values (δ³⁴S between –2.7‰ and –24.4‰) and oxygen isotope signatures support a model of magmatic-meteoric fluid interaction.

Geochemical analyses confirm a calc-alkaline to shoshonitic affinity of volcanic rocks derived from an enriched mantle source influenced by subduction-related fluids and crustal contamination. Mineralogy is dominated by primary chalcocite, bornite, and pyrite, with supergene malachite, azurite, and hematite. Alteration includes regional zeolite–carbonate and localized chlorite and hematite zones.

The Sepid-Sarve deposit exemplifies copper mineralization in an active continental margin, where subduction-related magmatism, structural controls (notably the Doruneh Fault), and hydrothermal processes converge to form significant Manto-type deposits in volcano-sedimentary environments