The Geology, Geochemistry, and Origin of the Porphyry Cu-Au-(Mo) System at Vathi, Serbo-Macedonian Massif, Greece

: The Vathi porphyry Cu-Au ± Mo mineralization is located in the Serbo-Macedonian metallogenic province of the Western Tethyan Metallogenic Belt. It is mainly hosted by a latite and is genetically associated with a quartz monzonite intrusion, which intruded the basement rocks of the Vertiskos Unit and the latite, 18 to 17 Ma ago. A phreatic breccia crosscuts the latite. The quartz monzonite was affected by potassic alteration, whereas the latite was subjected to local propylitic alteration. Both styles of alteration were subsequently overprinted by intense sericitic alteration. M-type and A-type veins are spatially associated with potassic alteration, whereas D-type veins are related to the sericitic alteration. Three ore assemblages are associated with the porphyry stage: (1) pyrite + chalcopyrite + bornite + molybdenite + magnetite associated with potassic alteration; (2) pyrite + chalcopyrite related to propylitic alteration; and (3) pyrite + chalcopyrite + native gold ± tetradymite associated with sericitic alteration. A fourth assemblage consisting of sphalerite + galena + arsenopyrite + pyrrhotite + pyrite ± stibnite ± tennantite is related to an epithermal overprint. Fluid inclusion data indicate that the A-type veins and related porphyry-style mineralization formed at 390–540 ◦ C and pressures of up to 646 bars (<2.6 km depth) from boiling hydrothermal ﬂuids. A later condensation of vapor-rich inclusions resulted in a moderately saline ﬂuid (8.4–11.2 wt % NaCl equiv) at temperatures between 311 and 392 ◦ C, which were related to sericitic alteration, D-type veins, and associated metallic mineralization. Subsequent dilution of the moderately saline ﬂuid at lower temperatures (205–259 ◦ C) produced a less saline (1.4–2.9 wt % NaCl equiv.) ﬂuid, which is likely associated with the late epithermal overprint. core Η4, 43 166 Optical microscopy on 44 thin, polished, and polished-thin from the surface along with the samples obtained from drill cores, Leitz Laborlux 11 Pol S dual reflected-transmitted light polarizing microscope. Scanning electron microscope (SEM) of energy-dispersive (EDS). Twenty-one samples of

The Vathi area is located in Kilkis ore district in northern Greece, which is a poorly explored region containing known porphyry-style deposits in the SMMP (Figures 1 and 2). Preliminary estimates of the resources at Vathi are 15 Mt ore grading 0.30 wt % Cu and 0.8 g/t Au [3,6]. This deposit has been exploited for gold and copper probably during the reign of Alexander I of Macedonia (498-454 BC) [18]. Remnants of surface excavations and galleries are scattered in the area, reflecting the ancient mining activity.  Figure 3 (modified after [6,[13][14][15][16]19,20]).
This study focuses on the porphyry system at Vathi and discusses the genetic and spatial relationships among hydrothermal alteration, geochemistry, metallic mineralization, metalliferous distribution, characteristics of the ore-forming fluids, and the structural control on the localization of the deposit. The integrated results of this study are used to discuss the metallogenic processes that formed the mineralization and to establish a genetic model for porphyry-style mineralization in the SMMP.  Figure 3 (modified after [6,[13][14][15][16]19,20]). places as sheeted veins, are associated with the E-and NE-trending group of structures and dip 60° to 76° toward SSW. They occur on the northeastern side of Ragian 1 Hill (Figures 3 and 4). At this site, an ancient gallery, approximately 1 m wide, was opened along these veins possibly for gold and/or copper exploitation ( Figure 4).

Figure 4.
Veins, locally sheeted, containing supergene mineralization with minor quartz and sericite that are associated with E-and NE-trending structures and which crosscut latite at the eastern side of Ragian 1 Hill.

Methods
The studied mineralization (N41.1458°, E22.9639°) is located northwest of the Vathi village in the Kilkis area and spans two hills, Ragian 1 Hill and Ragian 2 Hill (Figure 3a). Surface geological mapping and sampling resulted in the collection of 43 mineralized surface and hydrothermally altered samples from the igneous and surrounding metamorphic rocks. Geological mapping collected information on lithologies, lithological boundaries, and structural relationships. In addition, samples from legacy drill cores were obtained along with previously acquired logging data.
These legacy drill cores were collected in the late 1960s and early 1970s during an exploration program that was targeting uranium-bearing mineralization in the area. The Institute of Geology and Mineral Exploration (IGME) of Greece conducted an exploration drilling project during 1968. Two vertical drill holes (D1, D5; Figure 3a,b) and two inclined drill holes (D2, D3; Figure 3a,b), 512.90 m in length were drilled [46], with the results of geological mapping and geophysical exploration being reported by Melidonis [42], Thanasoulas [43], and Andronopoulos [44]. A further exploration drilling program was carried out by Bauxites Parnasse Mining Co. S.A. to target porphyry copper-style mineralization. This program included 11 vertical drill holes (H1, H4-13). Eight polished sections were made from drill core Η4, between 43 and 166 m, with seven additional polished sections being made from drill core Η5, between 53 and 73 m (H4, H5; Figure  3a).
Optical microscopy was conducted on 44 thin, polished, and polished-thin sections from the surface samples along with the samples obtained from drill cores, using a Leitz Laborlux 11 Pol S dual reflected-transmitted light polarizing microscope. Scanning electron microscope (SEM) studies of oxidized rock samples were conducted at the Faculty of Sciences, Aristotle University of Thessaloniki, using a JEOL 840 SEM equipped with an ISIS 300 OXFORD energy-dispersive spectrometer (EDS). Twenty-one samples of mineralization and the host rocks were analyzed for major and trace elements, including rare

Methods
The studied mineralization (N41.1458 • , E22.9639 • ) is located northwest of the Vathi village in the Kilkis area and spans two hills, Ragian 1 Hill and Ragian 2 Hill (Figure 3a). Surface geological mapping and sampling resulted in the collection of 43 mineralized surface and hydrothermally altered samples from the igneous and surrounding metamorphic rocks. Geological mapping collected information on lithologies, lithological boundaries, and structural relationships. In addition, samples from legacy drill cores were obtained along with previously acquired logging data.
These legacy drill cores were collected in the late 1960s and early 1970s during an exploration program that was targeting uranium-bearing mineralization in the area. The Institute of Geology and Mineral Exploration (IGME) of Greece conducted an exploration drilling project during 1968. Two vertical drill holes (D1, D5; Figure 3a,b) and two inclined drill holes (D2, D3; Figure 3a,b), 512.90 m in length were drilled [46], with the results of geological mapping and geophysical exploration being reported by Melidonis [42], Thanasoulas [43], and Andronopoulos [44]. A further exploration drilling program was carried out by Bauxites Parnasse Mining Co. S.A. to target porphyry copper-style mineralization. This program included 11 vertical drill holes (H1, H4-13). Eight polished sections were made from drill core H4, between 43 and 166 m, with seven additional polished sections being made from drill core H5, between 53 and 73 m (H4, H5; Figure 3a).
Optical microscopy was conducted on 44 thin, polished, and polished-thin sections from the surface samples along with the samples obtained from drill cores, using a Leitz Laborlux 11 Pol S dual reflected-transmitted light polarizing microscope. Scanning electron microscope (SEM) studies of oxidized rock samples were conducted at the Faculty of Sciences, Aristotle University of Thessaloniki, using a JEOL 840 SEM equipped with an ISIS 300 OXFORD energy-dispersive spectrometer (EDS). Twenty-one samples of mineralization and the host rocks were analyzed for major and trace elements, including rare earth elements, by X-ray fluorescence (XRF), and inductively coupled plasma-mass spectrometry (ICP-MS) at ACME Laboratories, Vancouver, Canada.
Fluid inclusions were studied in five doubly-polished wafers of quartz in veins containing metallic minerals. A LINKAM THM-600/TMS 90 heating-freezing stage coupled to a Leitz SM-LUX-POL microscope was used for the microthermometric measurements at the Department of Mineralogy, Petrology, Economic Geology, Aristotle University of Thessaloniki. Ice (H 2 O) and organic reference substances with known melting points were used to calibrate the stage. The upper limit of heating was at 600 • C, with a precision of the measurements ±1 • C during heating and ±0.2 • C during freezing. The software packages FLUIDS [47] and SoWat [48] were used for the estimation of the fluid compositions and properties. Densities were calculated from equations proposed by Driesner [49].

Petrography
The latite shows a porphyritic, fine to medium grained (locally coarse grained,~1 cm in size) texture and moderate to strong hydrothermal alteration (Figure 5a). It contains up to 40 vol % phenocrysts of quartz, K-feldspar, and minor plagioclase, and clinopyroxene, in a fine-grained ground mass dominated by quartz, K-feldspar, and sericite (Figure 5b,c). The main accessory minerals are magnetite, apatite, monazite, and zircon. The hydrothermal alteration assemblage includes biotite, K-feldspar, quartz, and minor titanite and magnetite related to potassic zone at depth; epidote, clinozoisite, chlorite, and albite related to propylitic zone; and sericite, quartz, pyrite, kaolinite, dolomite and rutile related to sericitic alteration. The surface of the latite bears a significant amount of totally oxidized metallic mineralization.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 39 earth elements, by X-ray fluorescence (XRF), and inductively coupled plasma-mass spectrometry (ICP-MS) at ACME Laboratories, Vancouver, Canada. Fluid inclusions were studied in five doubly-polished wafers of quartz in veins containing metallic minerals. A LINKAM THM-600/ΤΜS 90 heating-freezing stage coupled to a Leitz SM-LUX-POL microscope was used for the microthermometric measurements at the Department of Mineralogy, Petrology, Economic Geology, Aristotle University of Thessaloniki. Ice (H2O) and organic reference substances with known melting points were used to calibrate the stage. The upper limit of heating was at 600 °C, with a precision of the measurements ±1 °C during heating and ±0.2 °C during freezing. The software packages FLUIDS [47] and SoWat [48] were used for the estimation of the fluid compositions and properties. Densities were calculated from equations proposed by Driesner [49].

Petrography
The latite shows a porphyritic, fine to medium grained (locally coarse grained, ~1 cm in size) texture and moderate to strong hydrothermal alteration (Figure 5a). It contains up to 40 vol % phenocrysts of quartz, K-feldspar, and minor plagioclase, and clinopyroxene, in a fine-grained ground mass dominated by quartz, K-feldspar, and sericite ( Figure  5b,c). The main accessory minerals are magnetite, apatite, monazite, and zircon. The hydrothermal alteration assemblage includes biotite, K-feldspar, quartz, and minor titanite and magnetite related to potassic zone at depth; epidote, clinozoisite, chlorite, and albite related to propylitic zone; and sericite, quartz, pyrite, kaolinite, dolomite and rutile related to sericitic alteration. The surface of the latite bears a significant amount of totally oxidized metallic mineralization.  The quartz monzonite is a micro-to equi-granular, locally porphyritic intrusion grey to light brown in color (Figure 5d). It consists of plagioclase, quartz, biotite, clinopyroxene, and hornblende, with accessory zircon and apatite crystals less than 2 mm in size (Figure 5e,f). The hydrothermal alteration assemblage in the quartz monzonite consists of biotite, K-feldspar, magnetite, titanite (potassic alteration), chlorite and sericite (sericitic alteration). Metallic minerals include pyrite and chalcopyrite. Mineralized and unmineralized quartz veinlets (up to 2 cm wide) crosscut the rock (Figure 5g).
The phreatic breccia consists of medium to strongly hydrothermally altered fragments of gneiss, schist, and latite of various sizes (up to 30 cm long) and shapes (Figure 6a,b). The breccia is extensively oxidized and well cemented by an argillaceous matrix. Towards the margins of the breccia body the fragments are more abundant, coarse, and angular in shape. The matrix is composed of illite, kaolinite, sericite, quartz, and oxidized metallic mineralization. Textural and mineralogical characteristics such as the heterolithic wellcemented fragments, their sizes and shapes, the fine argillaceous matrix, the oxidized mineralization, and the absence of cavities, classify this breccia as phreatic, according to the criteria of Sillitoe [50].
The cataclasite is composed of angular quartz fragments (up to 2 cm) and of a dark colored, fine-grained, silicified matrix ( Figure 6c). Quartz is the main mineral while the hydrothermal alteration assemblage includes microcrystalline quartz and pyrophyllite.
The quartz monzonite is a micro-to equi-granular, locally porphyritic intrusion grey to light brown in color (Figure 5d). It consists of plagioclase, quartz, biotite, clinopyroxene, and hornblende, with accessory zircon and apatite crystals less than 2 mm in size (Figure 5e,f). The hydrothermal alteration assemblage in the quartz monzonite consists of biotite, K-feldspar, magnetite, titanite (potassic alteration), chlorite and sericite (sericitic alteration). Metallic minerals include pyrite and chalcopyrite. Mineralized and unmineralized quartz veinlets (up to 2 cm wide) crosscut the rock (Figure 5g).
The phreatic breccia consists of medium to strongly hydrothermally altered fragments of gneiss, schist, and latite of various sizes (up to 30 cm long) and shapes ( Figure  6a,b). The breccia is extensively oxidized and well cemented by an argillaceous matrix. Towards the margins of the breccia body the fragments are more abundant, coarse, and angular in shape. The matrix is composed of illite, kaolinite, sericite, quartz, and oxidized metallic mineralization. Textural and mineralogical characteristics such as the heterolithic well-cemented fragments, their sizes and shapes, the fine argillaceous matrix, the oxidized mineralization, and the absence of cavities, classify this breccia as phreatic, according to the criteria of Sillitoe [50].
The cataclasite is composed of angular quartz fragments (up to 2 cm) and of a dark colored, fine-grained, silicified matrix ( Figure 6c). Quartz is the main mineral while the hydrothermal alteration assemblage includes microcrystalline quartz and pyrophyllite.

Lithogeochemistry
The compositions of 15 samples of latite, three samples of quartz monzonite, and single samples of supergene mineralization, phreatic breccia, and cataclasite are shown in

Lithogeochemistry
The compositions of 15 samples of latite, three samples of quartz monzonite, and single samples of supergene mineralization, phreatic breccia, and cataclasite are shown in Tables 1 and 2. Due to the localized but extensive hydrothermal alteration of the magmatic rocks, the geochemical affinity and the classification were based on carefully selected fresh or least propylitized samples. Samples Vath 41 (quartz monzonite) and Vath 42 (latite; slightly propylitized), as well as unaltered quartz monzonite sample reported by Filippidis et al. [45] were evaluated (Figure 7a,b).
Samples Vath 41 and Vath 42 have similar silica contents of 63.05 wt % and 61.57 wt % SiO 2 , respectively (Table 1). They have a high total alkali content, ranging from 6.55 to 8.70 wt % K 2 O + Na 2 O, and plot in the alkali-calcic to slightly alkali fields in Figure 7a. The two magmatic rock types are classified as quartz monzonite and as latite (Figure 7b).
The altered latite samples are characterized by wide variations in major-element composition (

Hydrothermal Alteration
Three main types of hydrothermal alteration affected the magmatic and adjacent metamorphic rocks: potassic, propylitic, and sericitic ( Figure 9). No particular patterns are ascribed to the distribution of the hydrothermal alteration, which is mostly spatially associated with the quartz monzonite and regional structures.
Potassic alteration, spatially and temporally related to the quartz monzonite, is generally restricted in space due to the subsequent sericitic alteration overprint. Drill core samples show that it mainly occurs at depth but it is also weakly exposed at the surface [45] ( Figure 9). It was previously described by Frei [16] and Filippidis et al. [45] who recognized biotite-K-feldspar-pyrite veinlets in the drill core from the quartz monzonite stock, which is located 2 km north of Ragian 1 Hill (Figure 3; inset map). These veinlets are approximately 1 cm wide and consist of biotite, K-feldspar, magnetite, titanite, and quartz ( Figure 10a). Fine-grained hydrothermal biotite was also found in A-type veinlets and disseminated in the groundmass of quartz monzonite.
The propylitic alteration is generally weak and affects both the quartz monzonite and the latite (Figure 9). The characteristic hydrothermal minerals, also reported by Frei [16], are epidote, clinozoisite, chlorite, minor sericite, pyrite, and traces of albite ( Figure 10b).
The sericitic alteration affects both the magmatic and metamorphic rocks, and the phreatic breccia ( Figure 9). In the latite, it overprints the propylitic alteration and locally results in the complete replacement of feldspar, especially at the eastern part of the porphyry system. In addition, sericite, pyrite, quartz, kaolinite, dolomite, and rutile are also present in the altered latite ( Figure 10c). Sericite forms acicular to prismatic aggregates, which are accompanied by disseminations and masses of hematite. The K-feldspar replacement by sericite is generally controlled by the distribution of micro-fractures. In these cases, only the relic crystal habit of K-feldspar remains, forming vugs of different sizes. An epithermal overprint occurs locally within the sericitic alteration and is manifested as veinlets and as euhedral crystals (<30 µm in size) of hydrothermal quartz with fine-grained pyrophyllite in the groundmass of the latite. This epithermal event is mainly traced along an E-W trending structure on the northern side of Ragian 1 Hill into a latite cataclasite, which includes angular quartz fragments and pyrophyllite in a silicified groundmass (Figure 10d).
Appl. Sci. 2021, 11, x FOR PEER REVIEW 16 of 39 Figure 9. Spatial distribution of the hydrothermal alterations of the Vathi porphyry system. The locations of the drill holes are also shown (advice Figure 3). Figure 9. Spatial distribution of the hydrothermal alterations of the Vathi porphyry system. The locations of the drill holes are also shown (advice Figure 3).

Ore Mineralization
Ore mineralization at the Vathi porphyry system is hosted by quartz monzonite, latite, phreatic breccia, and surrounding metamorphic rocks. At the surface, the metallic mineralization is almost totally oxidized. The oxidation zone has a maximum depth of 90 to 100 m and is associated with a narrow secondary enrichment zone, 2 m deep [46] (Figure 3c). The hypogene mineralization at the surface is restricted to the quartz monzonite, and in the drill core samples, it is well exposed below the secondary enrichment zone to an unknown depth.

Ore Mineralization
Ore mineralization at the Vathi porphyry system is hosted by quartz monzonite, latite, phreatic breccia, and surrounding metamorphic rocks. At the surface, the metallic mineralization is almost totally oxidized. The oxidation zone has a maximum depth of 90 to 100 m and is associated with a narrow secondary enrichment zone, 2 m deep [46] ( Figure 3c). The hypogene mineralization at the surface is restricted to the quartz monzonite, and in the drill core samples, it is well exposed below the secondary enrichment zone to an unknown depth.
The mineralization associated with the potassic alteration occurs in quartz monzonite and is present mainly along drill hole D1 (66-75 m interval) [46]. It forms A-and M-type veins and disseminations, and it consists mostly of pyrite and chalcopyrite with minor molybdenite, bornite, and magnetite. A-type quartz veins (up to 2 cm wide) are rare but were observed in the drill core samples (Figure 12a). The M-type veins, found on the surface and at the depth, contain magnetite with quartz and titanite (Figure 12b). At a depth of 140 m, pyrite and chalcopyrite are present in approximately equal amounts, whereas at 166 m, only chalcopyrite and molybdenite occur. Restricted amounts of pyrite
The mineralization associated with the potassic alteration occurs in quartz monzonite and is present mainly along drill hole D1 (66-75 m interval) [46]. It forms A-and M-type veins and disseminations, and it consists mostly of pyrite and chalcopyrite with minor molybdenite, bornite, and magnetite. A-type quartz veins (up to 2 cm wide) are rare but were observed in the drill core samples (Figure 12a). The M-type veins, found on the surface and at the depth, contain magnetite with quartz and titanite (Figure 12b). At a depth of 140 m, pyrite and chalcopyrite are present in approximately equal amounts, whereas at 166 m, only chalcopyrite and molybdenite occur. Restricted amounts of pyrite and chalcopyrite were observed in the propylitic alteration that extends in both magmatic rock types at the surface and at a depth of 40 to 140 m.
Pyrite and chalcopyrite, as well as minor amounts of rutile and ilmenite, occur in sericitic alteration mainly in latite but also in the phreatic breccia and surrounding metamorphic rocks [46]. The mineralization forms disseminations, aggregates, and D-type veins and veinlets and associated halos along cracks in the latite rock mass (Figure 12c). In some cases, D-type sheeted veins trending E-W form a dense system crosscutting the latite ( Figure 6). They are usually less than 1 cm wide, locally reaching up to 5 cm in width.A subsequent epithermal overprint is documented locally within the potassic and sericitic alteration. It mainly occurs as E-type veins at depth and consists of sphalerite, galena, arsenopyrite, pyrrhotite, and traces of tennantite, stibnite, rutile, and ilmenite. matic rock types at the surface and at a depth of 40 to 140 m.
Pyrite and chalcopyrite, as well as minor amounts of rutile and ilmenite, occur in sericitic alteration mainly in latite but also in the phreatic breccia and surrounding metamorphic rocks [46]. The mineralization forms disseminations, aggregates, and D-type veins and veinlets and associated halos along cracks in the latite rock mass (Figure 12c). In some cases, D-type sheeted veins trending E-W form a dense system crosscutting the latite ( Figure 6). They are usually less than 1 cm wide, locally reaching up to 5 cm in width.A subsequent epithermal overprint is documented locally within the potassic and sericitic alteration. It mainly occurs as E-type veins at depth and consists of sphalerite, galena, arsenopyrite, pyrrhotite, and traces of tennantite, stibnite, rutile, and ilmenite. Pyrite is the most common sulfide in the Vathi porphyry system, and it is intergrown with chalcopyrite, bornite, pyrrhotite, sphalerite, and galena. Locally, it hosts rounded inclusions of chalcopyrite, as well as euhedral inclusions of magnetite and mo- Pyrite is the most common sulfide in the Vathi porphyry system, and it is intergrown with chalcopyrite, bornite, pyrrhotite, sphalerite, and galena. Locally, it hosts rounded inclusions of chalcopyrite, as well as euhedral inclusions of magnetite and molybdenite. It occurs as disseminations or intergrowths of anhedral to euhedral grains, with sizes varying from 25 to 200 µm and from 400 µm to 2 cm, respectively (Figure 13a). Pyrite also forms veinlets (up to 3 mm in size) and minor dispersions, in which galena locally appears in the interstices (Figure 13b).
Chalcopyrite is the second most abundant sulfide and is found mainly in the potassic alteration where it is intergrown with pyrite, magnetite, molybdenite, and bornite (Figure 13c-f). It occurs as disseminated grains (<200 µm in size), irregular aggregates (<300 µm in size), and rounded inclusions (<20 µm in diameter) in pyrite. Chalcopyrite also forms thin veinlets (<0.5 cm in width) along with pyrite inclusions (10-20 µm in size). Below a depth of 166 m, mineralized stockworks with chalcopyrite and molybdenite are present (Figure 13e).
with sizes varying from 25 to 200 μm and from 400 μm to 2 cm, respectively (Figure 13a). Pyrite also forms veinlets (up to 3 mm in size) and minor dispersions, in which galena locally appears in the interstices (Figure 13b).
Chalcopyrite is the second most abundant sulfide and is found mainly in the potassic alteration where it is intergrown with pyrite, magnetite, molybdenite, and bornite (Figure 13c-f). It occurs as disseminated grains (<200 μm in size), irregular aggregates (<300 μm in size), and rounded inclusions (<20 μm in diameter) in pyrite. Chalcopyrite also forms thin veinlets (<0.5 cm in width) along with pyrite inclusions (10-20 μm in size). Below a depth of 166 m, mineralized stockworks with chalcopyrite and molybdenite are present (Figure 13e). Magnetite is associated with the potassic alteration forming M-type veinlets (up to 1 cm wide) crosscutting the quartz monzonite (Figure 12b). It also occurs as disseminated crystals or intergrowths with chalcopyrite. In latite, massive magnetite is intergrown with pyrite, chalcopyrite, and minor titanite (Figure 13d). It is also found in chalcopyrite as subhedral to euhedral inclusions (Figure 13e).
Molybdenite occurs with pyrite and chalcopyrite in potassic alteration as small (<30 μm in length), disseminated, lamellar grains in latite (Figure 13e). In some cases, molybdenite also occurs as disseminations intergrown with prismatic rutile. Bornite in potassic alteration locally forms rounded inclusions, ranging between 15 and 100 μm in size in pyrite and chalcopyrite (Figure 13f). Magnetite is associated with the potassic alteration forming M-type veinlets (up to 1 cm wide) crosscutting the quartz monzonite (Figure 12b). It also occurs as disseminated crystals or intergrowths with chalcopyrite. In latite, massive magnetite is intergrown with pyrite, chalcopyrite, and minor titanite (Figure 13d). It is also found in chalcopyrite as subhedral to euhedral inclusions (Figure 13e).
Pyrrhotite is present at a depth below 166 m in E-type veins and is intergrown with pyrite and chalcopyrite (1 cm wide) (Figure 14a). Galena occurs in sericitic alteration and is intergrown with pyrite, chalcopyrite (<50 µm in width), and sphalerite. It appears in the interstices in pyrite veinlets and disseminations in quartz monzonite (Figure 13b), or it forms E-type veins (1 to 3 mm wide) with sphalerite and quartz (Figure 14b,c). Pyrrhotite is present at a depth below 166 m in Ε-type veins and is intergrown with pyrite and chalcopyrite (1 cm wide) (Figure 14a). Galena occurs in sericitic alteration and is intergrown with pyrite, chalcopyrite (<50 μm in width), and sphalerite. It appears in the interstices in pyrite veinlets and disseminations in quartz monzonite (Figure 13b), or it forms E-type veins (1 to 3 mm wide) with sphalerite and quartz (Figure 14b,c). Sphalerite, similar to galena, occurs as large anhedral grains intergrown with pyrite and chalcopyrite in sericitic alteration. In E-type veins, sphalerite contains numerous chalcopyrite inclusions (so-called chalcopyrite disease) and is intergrown with chalcopyrite, galena, arsenopyrite, minor pyrite, and quartz (Figure 14b,c). Arsenopyrite was observed only in the E-type veins along with sphalerite and galena. It mainly forms prismatic and rhombohedral crystals, <500 µm in length (Figure 14d), and it occurs in micro-cavities in latite along with quartz and pyrite. Sphalerite, similar to galena, occurs as large anhedral grains intergrown with pyrite and chalcopyrite in sericitic alteration. In E-type veins, sphalerite contains numerous Stibnite is present in traces (<50 µm in size) in sericitic alteration in latite (Figure 14e). Additionally, stibnite occurs in the E-type veins as tiny radiating and elongated crystals in cavities or as intergrowths with sphalerite and galena (Figure 14c). Tennantite is spatially related to sphalerite and galena in sericitic alteration, where it appears as tiny inclusions in sphalerite with galena and chalcopyrite as well as overgrowths on pyrite, chalcopyrite, and bornite (Figure 13c). Native gold occurs as inclusions 10 to 50 µm in size in pyrite and chalcopyrite and as infills between grains in pyrite veinlets related to sericitic alteration (Figure 14f). It is also found in goethite. In D-type pyrite veinlets, native gold (up to 20 µm in length) appears as infills along with chalcopyrite and tetradymite (Figure 14g). Veinlets of gold in goethite were reported previously by Markoulis (1970).
Tetradymite is a minor mineral up to 7 µm in size along with chalcopyrite and native Au in D-type pyrite veins in quartz monzonite where potassic alteration is overprinted by sericitic alteration (Figure 14h). Rutile occurs as prismatic crystals up to 50 µm in size in pyrite, chalcopyrite, and molybdenite in potassic and sericitic alterations. Ilmenite forms lamellar exsolutions in magnetite in the potassic alteration zone. It is also present in pyrite and chalcopyrite as euhedral crystals 20 to 40 µm in length.

Supergene Mineralization
The oxidation zone in latite contains masses, encrustations, and impregnations of secondary Cu-and Fe-minerals, which are characterized by traces of native Au and Cu. The supergene enrichment zone hosts mainly chalcocite and covellite as well as minor cuprite and native copper ( Figure 11).
Hematite and goethite are associated to pyrite and chalcopyrite oxidation (Figure 13a). Hematite is found only at the surface where it forms sheets and botryoidal to granular aggregates. In rock fractures or micro-cavities, hematite appears as radial and bushy aggregates of platy sheets, commonly with sericite and rutile. Goethite is the dominant mineral in the fault zones where it forms massive aggregates. In the oxidation zone, goethite replaces pyrite and chalcopyrite, forming pseudomorphs after pyrite and aggregates in micro-cavities, respectively.
Malachite and minor azurite are mostly present in the oxidation zone of the latite and in the phreatic breccia. Malachite appears as massive, botryoidal aggregates and acicular to tabular-prismatic crystals filling interstices and fractures in the rock, whereas azurite forms encrustations.
Cuprite increases with depth in the oxidation zone [46]. It is mainly found near the supergene enrichment zone at a depth of 90 to 92 m in veinlets with minor chalcocite and covellite. Native copper forms small disseminated grains (up to 50 µm in size) in cuprite. Chalcocite and covellite occur in the secondary enrichment zone where they are most abundant at depths between 90 and 100 m.
Although relatively rare in Greece, turquoise occurs on the eastern side of Ragian 1 Hill [57,58] (Figure 15a,b). It occurs in veinlets and encrustations exhibiting characteristic pale green-blue and forms hexagonal pseudomorphs after apatite [58]. In addition to turquoise, various minor occurrences of woodhouseite and phosphate carbonates, rich in large-ion lithophile elements (LILE) and HFSE (mainly Sr, Nb, Ce, Th, La, U), as well as senaite, are intergrown with sericite or they form in the rock interstices (Figure 15b).
The autunite group minerals, torbernite and meta-torbernite, occur as disseminations and encrustations in the oxidation zone of latite. On the eastern side of Ragian 1 Hill, pale green tabular torbernite and dark-green flat plates of meta-torbernite (crystals up to 3 mm in size) are locally present as disseminations in latite [44,59,60]. A significant concentration of torbernite occurs in drill hole D2 between 72 and 75 m [46] (Figure 3). Torbernite appears as disseminations, up to 3 mm in sizes, and as foliated to earthy aggregations of tabular crystals up to 3 cm in size [46].  (Figure 3). Torbernite appears as disseminations, up to 3 mm in sizes, and as foliated to earthy aggregations of tabular crystals up to 3 cm in size [46].

Geochemistry of Ore Mineralization
Based on the chemical analyses of 21 mineralized surface samples (Table 1), the Vathi porphyry system is locally enriched in Cu (avg: 2313 ppm), Mo (avg: 152 ppm), Au (avg: 0.46 ppm), U (avg: 31 ppm), Pb (avg: 120 ppm), Zn (avg: 55 ppm), and Ag (avg: 1 ppm). This enrichment in base and precious metals is spatially related to hydrothermal alteration and to the distribution of metallic minerals. On the surface, the highest grades of Au, Ag, Cu, Mo, Pb, Zn, and U are mainly associated with the intense sericitic alteration and oxidized D-type veins on the eastern side of the latite at Ragian 1 Hill ( Figure  16). Weak anomalies of gold are present along and around stocks of quartz monzonite.
The propylitic alteration in latite generally contains a paucity of metal (only up to 116.8 ppm Cu, 5.16 ppm Mo, 36.3 ppb Au, 11.15 ppm U, 306.7 ppm Pb, 54.0 ppm Zn, and 0.7 ppm Ag; Table 1). Samples affected by sericitic alteration exhibit higher metal concentrations. The highest metal concentrations in samples VTH 6 and VTH 11 reach up to 9297 ppm Cu, 319.1 ppm Mo, 6996 ppb Au, 262.6 ppm U, 977.3 ppm Pb, 212 ppm Zn, and 4.2 ppm Ag (Table 1). This enrichment is associated with the sericitic alteration and secondary oxidation.
The phreatic breccia is also characterized by high U (40.5 ppm) and REE contents especially La (613.

Geochemistry of Ore Mineralization
Based on the chemical analyses of 21 mineralized surface samples (Table 1) (Table 1). This enrichment is associated with the sericitic alteration and secondary oxidation.
The phreatic breccia is also characterized by high U (40.5 ppm) and REE contents especially La (613.

Fluid Inclusion Study
A fluid inclusion study was conducted on quartz from A-and E-type veins and on hydrothermal quartz grains related to sericitic alteration associated with the porphyry mineralization and the epithermal overprint, respectively. Samples of A-type veins crosscutting quartz monzonite were obtained from the surface (Vath 31c) and from drill core H4

Fluid Inclusion Study
A fluid inclusion study was conducted on quartz from A-and E-type veins and on hydrothermal quartz grains related to sericitic alteration associated with the porphyry mineralization and the epithermal overprint, respectively. Samples of A-type veins crosscutting quartz monzonite were obtained from the surface (Vath 31c) and from drill core H4 (Vath 40). Samples of E-type veins were collected from cataclasite at the surface (Vath 33) and from drill cores H4 (sample Vath 41) and H5 (sample Vath 42). Hydrothermal quartz in the sericitic alteration was derived from the surface (sample Vath 12a). D-type veins are devoid of fluid inclusions.
Quartz in the veins and latite is mainly equigranular, anhedral, and it possesses a polygonal shape. The contacts among quartz grains form triple junctions with angles at ≈120° (Figure 17a). This implies that the quartz was extensively recrystallized so that the fluid inclusions largely migrated to grain boundaries [61]. As a result, fluid inclusions at Vathi are scarce, and only five of eight samples contained workable inclusions.
Where present, fluid inclusions are concentrated in clusters ( Figure 17b) and are considered to be primary, on the basis of the criteria proposed by Roedder [62] and Goldstein and Reynolds [63]. These criteria include constant liquid to vapor ratios and similar shapes and sizes of the inclusions. Primary fluid inclusions found along growth zones of quartz crystals are rare, and their size is too small for microthermometry. Inclu- Quartz in the veins and latite is mainly equigranular, anhedral, and it possesses a polygonal shape. The contacts among quartz grains form triple junctions with angles at ≈120 • (Figure 17a). This implies that the quartz was extensively recrystallized so that the fluid inclusions largely migrated to grain boundaries [61]. As a result, fluid inclusions at Vathi are scarce, and only five of eight samples contained workable inclusions.
Where present, fluid inclusions are concentrated in clusters ( Figure 17b) and are considered to be primary, on the basis of the criteria proposed by Roedder [62] and Goldstein and Reynolds [63]. These criteria include constant liquid to vapor ratios and similar shapes and sizes of the inclusions. Primary fluid inclusions found along growth zones of quartz crystals are rare, and their size is too small for microthermometry. Inclusions that occur in planes along healed fractures (Figure 17c) are considered to be secondary and possibly reflect fluid circulation associated with a late stage of fluid entrapment. The primary inclusions are elliptical to rounded and sometimes irregular in shape. Their maximum diameter is 35 µm, but the vast majority range in size from 5 to 15 µm.
On the basis of phase ratios observed at room temperature and phase transitions during homogenization, four types of fluid inclusions are present in the early-stage A-type veins and the hydrothermal quartz of the host rock: two-phase liquid-rich aqueous (Type 1), multiphase halite-bearing (Type 2), two-phase vapor-rich aqueous (Type 3), and one-phase vapor only (Type 4) inclusions.
Their maximum diameter is 35 μm, but the vast majority range in size from 5 to 15 μm.
On the basis of phase ratios observed at room temperature and phase transitions during homogenization, four types of fluid inclusions are present in the early-stage A-type veins and the hydrothermal quartz of the host rock: two-phase liquid-rich aqueous (Type 1), multiphase halite-bearing (Type 2), two-phase vapor-rich aqueous (Type 3), and one-phase vapor only (Type 4) inclusions.  Type 1 two-phase aqueous inclusions contain a vapor bubble occupying 20-30 vol.% (Figure 17d), which homogenize into the liquid phase during heating. These inclusions occur in A-type veins and are rare and secondary in origin, although some primary inclusions occur in the hydrothermal quartz grains in latite.
Type 2 inclusions are multiphase and consist of liquid, vapor, and one or more daughter crystals consisting mostly of halite and rare sylvite, or other unidentified transparent daughter minerals or an opaque (Figure 17e-h). The vapor phase occupies 20-40 vol.% of the inclusions, which homogenize into the liquid phase. Halite and sylvite, where halite is generally larger than sylvite, were identified by their colorless, isotropic cubic crystal morphology (Figure 17e,f). Halite dissolved at a higher temperature than sylvite during heating. The opaque phases (Figure 17f) are either cubic or triangular in shape and are possibly pyrite or chalcopyrite, respectively, but they do not dissolve during heating up to 600 • C.
Type 3 two-phase inclusions are dominated by a vapor bubble comprising 70-90 vol.% of the inclusions (Figure 17b,g,h) and homogenize into the vapor phase. They are the most common type of inclusion. Type 3 fluid inclusions commonly locally coexist with type 2 inclusions (Figure 17h). Type 4 inclusions are usually dark in color and consist mainly of vapor. They occur along healed fractures and are considered to be secondary in origin (Figure 17c,i).
None of the inclusions analyzed at the Vathi porphyry system contain visible liquid CO 2 , and clathrate formation was not observed during freezing experiments. This excludes the presence of significant quantities of CO 2 (<2 mol.%).
Quartz in the epithermal veins contains few fluid inclusions, but where present, they are <20 µm in length. They are two-phase aqueous liquid-vapor inclusions (Type 1) containing a vapor bubble occupying 10-20 vol.% (Figure 17j). These inclusions homogenized into the liquid phase upon heating. Post-entrapment modifications, such as necking down or leakage, have affected this type of fluid inclusions (Figure 17k,l). They commonly occur along micro-cracks and demonstrate a variable liquid-to-vapor ratio. This is likely due to "stretching" processes during mechanical intracrystalline straining [64]. In order to eliminate the effects of the post-entrapment modifications, only inclusions with constant liquid-to-vapor ratios were selected for the microthermometric measurements.
Microthermometric data from 78 fluid inclusions in hydrothermal quartz and A-type veins from the Vathi porphyry system are summarized in Table 3. Eutectic ice-melting temperatures (Te) obtained from a few type 1 inclusions are concentrated around −22 • C, which suggests that NaCl is the dominant salt species in the inclusions [65]. Eutectic temperatures of type 3 inclusions were not possible to determine due to their small size.
Final ice melting temperatures of type 1 inclusions vary from −7.6 to −5.4 • C, indicating salinities of 8.4 to 11.2 wt % NaCl equiv [66]. Due to the small size (mainly 5-11 µm) of the meniscus liquid phase in type 3 inclusions, which made observation of the phase changes difficult, only a few freezing data could be obtained. On the basis of the freezing point depression (−14.1 to −6.4 • C), the salinities ranged from 9.7 to 18.5 wt % NaCl equiv.
Halite in type 2 fluid inclusions of A-type quartz veins dissolved at temperatures between 220 and 401 • C, corresponding to salinities of 32.7-47.2 wt % NaCl equiv. Salinity estimates based on the dissolution temperature of halite and sylvite in three type 2 inclusions yielded values between 28.4 and 36.7 wt % NaCl and 18.8 to 20.3 wt % KCl, which corresponds to a salinity of between 38.4 and 46.1 wt % NaCl equiv, respectively.
Homogenization temperatures (Th) for type 1 liquid-rich inclusions in quartz range from 311 to 392 • C ( Figure 18). Type 2 liquid-rich halite-bearing fluid inclusions homogenized by disappearance of the vapor bubble after the halite crystal had disappeared at temperatures from 394 to >600 • C (the upper operating limit of the fluid inclusion stage). Homogenization temperatures of type 3 vapor-rich inclusions range from 335 to >600 • C. The majority of values of Th for type 2 and 3 fluid inclusions overlap between 390 and 540 • C and have a distinct maximum at ≈510 • C ( Figure 18).
Microthermometric results for quartz in E-type epithermal veins were obtained from type 1 fluid inclusions in quartz (Table 3). Eutectic ice-melting temperatures range between −22 and −21 • C, suggesting that NaCl is the only dissolved salt in the fluid [65]. Final ice melting temperatures (−1.7 to −0.8 • C) correspond to salinities of 1.4 to 2.9 wt % NaCl equiv. Homogenization temperatures range from 205 to 259 • C, with a distinct maximum at ≈220 • C (Figure 18).  Microthermometric results for quartz in E-type epithermal veins were obtained from type 1 fluid inclusions in quartz (Table 3). Eutectic ice-melting temperatures range between −22 and −21 °C, suggesting that NaCl is the only dissolved salt in the fluid [65]. Final ice melting temperatures (−1.7 to −0.8 °C) correspond to salinities of 1.4 to 2.9 wt % NaCl equiv. Homogenization temperatures range from 205 to 259 °C, with a distinct maximum at ≈220 °C ( Figure 18).

Geotectonic Setting, Structural Control, and Magma Emplacement
The magmatic complex at Vathi was emplaced during post-subduction tectonic extension. During the Miocene, several deep to shallow subvolcanic stocks, plugs, and dikes, which vary from intermediate to felsic compositions, intruded the metamorphic basement rocks in the Vathi area [16,45]. Their emplacement was controlled by regional structures including E-W, NW-SE, and NE-SW trending faults, which also favored the formation of the porphyry-style mineralization by localizing the magmatic-hydrothermal fluid circulation (Figure 3; inset map).
The mineralization at Vathi is mainly hosted in latite, but it is genetically related to quartz monzonite, which intruded the older latite and the crystalline rocks of the Vertiskos Unit. The subvolcanic intrusion formed between 18 ± 0.5 and 17 ± 1 Ma, based on U-Pb zircon ages of Frei [16]. The microgranular and porphyritic texture and the elongate shape of the ore-related quartz monzonite are ascribed to fast ascending magma and subsequent crystallization at shallow depths set in an extensional regime. Melfos et al. [34] describe a similar mechanism for the origin of microgranite porphyry associated with the Maronia Cu-Mo deposit, NE Greece.
The presence of the phreatic breccia intruding the latite at Ragian 1 Hill, its ESE-WNW orientation, and the minor presence of quartz monzonite fragments suggest that quartz monzonite intruded the latite. This breccia was possibly formed when magmatic fluids were mixed with surface meteoric water during the late stages of the quartz monzonite intrusion, at the upper and marginal parts of the porphyry system [50,56].

Hydrothermal Alteration
Potassic and sericitic alteration is spatially related to the main stage of sulfide mineralization, whereas propylitic alteration hosts insignificant amounts of pyrite and chalcopyrite. The potassic alteration formed during the emplacement of the quartz monzonite and is characterized by the presence of quartz, hydrothermal biotite and K-feldspar, magnetite, chalcopyrite, bornite, molybdenite, and pyrite.
Propylitic alteration (epidote and chlorite) was likely coincident with the late stage of formation of the potassic alteration. Pyrite and chalcopyrite are found in trace amounts in the propylitic alteration, suggesting a slight increase in the sulfur fugacity of the hydrothermal fluids and/or a rise of the fluid temperature [67].
Sericitic alteration is expressed by the almost complete replacement of magmatic K-feldspar, clinopyroxene, hornblende, and biotite and by the extensive presence of sericite and quartz in the ground mass. Kaolinite, dolomite, and rutile constitute minor phases. Acicular to prismatic rutile crystals were formed as a result of Ti removal from precursor biotite during alteration [68]. The sericitic alteration overprints the potassic and propylitic alteration. The areas of sericitic alteration with the intense feldspar destructive textures possess relatively high sulfur content (S = ≈1.30 wt %). The magmatic fluid activity that resulted in the formation of this alteration zone during the late phases of the porphyry system was mainly controlled by the major E-W trending fault at the eastern part of Ragian 1 Hill. In this location, minor quartz veins and abundant mineralized D-type veins crosscut the host rocks.

Physicochemical Parameters of the Hydrothermal Fluids
Exsolution of hydrous fluids and silicate melts during the crystallization of the quartz monzonite at a shallow depth favored the hydraulic fracturing of the quartz monzonite and latite, which enhanced the permeability of the host rocks and the deposition of magnetite, chalcopyrite, bornite, molybdenite, and pyrite in A-type veins, fractures, and as disseminations in the host rock. This stage of mineralization is associated with the potassic alteration.
Fluid inclusions in A-type quartz veins at Vathi show that quartz and related metallic mineralization formed as a result of phase separation from an initial high temperature (at least >600 • C), featuring moderate-salinity supercritical magmatic fluid (Figures 19 and 20) under boiling conditions during the ascent and cooling of the fluids. The rapid drop in temperature and pressure favored fluid immiscibility in these veins, resulting in a saline brine (33-47 wt % NaCl equiv; type 2) and a co-existing lower salinity (10-19 wt % NaCl equiv; type 3) low-density vapor-rich fluid at lithostatic pressures that homogenize over the same temperature range, mainly between 390 and 540 • C. This association, plus the presence of type 4 vapor-rich inclusions, suggests that the potassic alteration was associated with fluid boiling (Figure 19) as has been described elsewhere [62,66,69,70].
Fluid boiling can cause extensive hydraulic fracturing of the wall rocks and result in mixing of ascending boiling magmatic fluid with meteoric water at the ductile to brittle transition under hydrostatic pressures [72]. The lower temperature (310-390 °C) and less saline (8-11 wt % NaCl equiv) type 1 liquid-rich fluid inclusions in the hydrothermal quartz of latite at Vathi (Figure 19) likely represents a late fluid that was associated with decompression under hydrostatic pressures ( Figure 20). This decompression was caused by extensive hydrofracturing and generated the fractures associated with the metallic mineralization at Vathi [69], which is associated with the sericitic alteration.  The presence of halite and sylvite in fluid inclusions indicates that the ore-bearing fluid in A-type veins was saturated in NaCl and KCl, as well as the likely presence of chalcopyrite and pyrite in type 2 inclusions, shows that the fluids were hypersaline and that Cu at Vathi was likely precipitated by the condensation of metal-rich brines.
Fluid inclusions trapped from boiling fluids do not require pressure corrections, and the measured homogenization temperatures represent trapping temperatures [62,70]. Trapping pressures of co-existing type 2 and 3 fluid inclusions in A-type veins at Vathi range from 248 to 646 bars, which were estimated using the equation of state of Zhang and Frantz [71] and homogenization temperatures. This corresponds to a depth of ore formation between 1.0 and 2.6 km assuming lithostatic conditions ( Figure 20). The estimated densities of the hydrothermal fluids are >0.77 g cm −3 for the highly saline liquid-rich (type 2) fluids, and <0.78 g cm −3 for the moderate-salinity vapor-rich (type 3) fluids.
Fluid boiling can cause extensive hydraulic fracturing of the wall rocks and result in mixing of ascending boiling magmatic fluid with meteoric water at the ductile to brittle transition under hydrostatic pressures [72]. The lower temperature (310-390 • C) and less saline (8-11 wt % NaCl equiv) type 1 liquid-rich fluid inclusions in the hydrothermal quartz of latite at Vathi (Figure 19) likely represents a late fluid that was associated with decompression under hydrostatic pressures ( Figure 20). This decompression was caused by extensive hydrofracturing and generated the fractures associated with the metallic mineralization at Vathi [69], which is associated with the sericitic alteration.
Fluid inclusions from the E-type epithermal veins are characterized by low salinity (1.4-2.9 wt % NaCl equiv) and relatively low Th ranging between 205 and 259 • C, with a maximum at ≈220 • C (Figures 19 and 20). These fluid inclusions formed during the final stage of the magmatic-hydrothermal evolution of the system and represent a dilute fluid resulting from the mixing of a moderate saline fluid with circulating meteoric water in the epithermal stage. These non-boiling fluids are associated with the assemblage sphalerite + galena + arsenopyrite ± stibnite ± tennantite ore. Therefore, pressures obtained for non-boiling fluid inclusion assemblages represent minimum trapping values, and the pressure correction is estimated to be <10 • C, since these veins probably formed at shallow depths (e.g., [73]). The estimated minimum trapping pressures of the Vathi epithermal veins range from 100 to 244 bars, indicating a depth between 110 and 250 m under hydrostatic conditions. The estimated densities of the hydrothermal fluids in these veins are >0.79 g cm −3 .
(1.4-2.9 wt % NaCl equiv) and relatively low Th ranging between 205 and 259 °C, with a maximum at ≈220 °C (Figures 19 and 20). These fluid inclusions formed during the final stage of the magmatic-hydrothermal evolution of the system and represent a dilute fluid resulting from the mixing of a moderate saline fluid with circulating meteoric water in the epithermal stage. These non-boiling fluids are associated with the assemblage sphalerite + galena + arsenopyrite ± stibnite ± tennantite ore. Therefore, pressures obtained for non-boiling fluid inclusion assemblages represent minimum trapping values, and the pressure correction is estimated to be <10 °C, since these veins probably formed at shallow depths (e.g., [73]). The estimated minimum trapping pressures of the Vathi epithermal veins range from 100 to 244 bars, indicating a depth between 110 and 250 m under hydrostatic conditions. The estimated densities of the hydrothermal fluids in these veins are >0.79 g cm −3 .

Source of Metals
Phase separation promotes Au deposition in porphyry systems, which follows the deposition of Cu at high temperatures [76,77]. In addition, Au precipitation in phreatic breccias is favored by the permeable textures created by the periodic catastrophic release of the build-up of volatiles [56]. This could also be the case for the Au enrichment during the formation of the Vathi breccia (724 ppb Au). At Vathi, gold is also related to the oxidized D-type veins, which occur in the sericitic alteration (6996 ppb Au). A similar relationship between high Au contents and D-type mineralized veins was reported by Voudouris et al. [78], where gold is found as inclusions in pyrite in the Pagoni Rachi Cu-Mo-Re-Au deposit in northern Greece.
Molybdenum at Vathi was precipitated near and below a depth of 166 m after Cu was precipitated, while Fe, Pb, Zn, and As are enriched in the upper, outer, and adjacent areas of the quartz monzonite and was related to the epithermal stage. Frei [16], based on Pb isotope compositions, has shown that the enrichment of Cu, Pb, and Zn at Vathi is associated with crustal rocks assimilation from magma.
The  Table 2). These elements are lithophile and incompatible in the mantle, and the enrichment of these elements could be related to differentiation of the parental magma and the subsequent interaction with crustal rocks, as well as to the onset of significant plagioclase fractionation [79,80]. In this regard, Sr values are considered more enriched, with respect to reported values from intermediate to evolved magmas (<100 ppm) [80]. Moreover, the enrichment in high field strength elements (HFSE; Zr, Hf, Nb, and Ta) is likely the result of limited assimilation of crustal rocks as the magma differentiated [79].
A significant enrichment is also observed in the REE: La (up to 613.1 ppm), Ce (up to 894 ppm), Pr (up to 79.1 ppm), Nd (up to 211.1 ppm), and Gd (up to 12.7 ppm), as well as in Ti (up to 7040 ppm) and P (up to 1340 ppm) ( Table 2). These geochemical anomalies could suggest that the Vathi deposit shares several characteristics in common with alkaline igneous rock-related porphyry deposits [81].
Uranium exhibits moderately to highly enriched concentrations in the magmatic rocks in the Vathi area (14 to 106 ppm). In the phreatic breccia and along the mineralized tectonic structures, the high U content (up to 328 ppm) is manifested by the presence of autunite-group minerals, torbernite and meta-torbernite. Uranium enrichment is rather an uncommon characteristic of porphyry systems associated with calc-alkaline igneous rocks [82]. On the contrary, U tends to be enriched in porphyry deposits related to alkaline intrusions [82,83].
At the Vathi area, the enrichment in U was likely associated with the final stages of crystallization of the mildly alkaline magma combined with movement along deep-seated faults. Rupturing of the intrusive cupola was likely associated with fluid exsolution and the accumulation of a magmatic fluid enriched in U. Fluid exsolution was favored by deep-seated faults, while a network of fractures enhanced rock permeability resulting in the formation of secondary torbernite. Rice et al. [83] described a similar mechanism that resulted in the enrichment of U in the alkaline igneous rock-related porphyry molybdenum system at Central City, Colorado.

Origin of the Vathi Deposit
In the early Cenozoic, melts enriched in volatiles, Cu, and other chalcophile metals were extracted from the subducted oceanic crust in northern Greece [84]. These melts ascended into the overlying delaminated crust, inducing partial melting. This melting resulted in significant crustal rock assimilation and produced magmas of calc-alkaline to high K-calc-alkaline and of shoshonitic affinities, which then intruded the upper crust, resulting in post-collisional back-arc magmatism [26]. An extensional tectonic regime, characterized by detachment and strike-slip faults in transpressional and transtensional tectonic settings, characterizes the Oligocene-Miocene deformation of the upper crust [19]. This deformation structurally controlled the formation of the I-type plutonic rocks and the related intrusive and extrusive magmatic rocks [84].
The oxidizing character of the magma that formed the ore deposit and the quartz monzonite association with porphyry copper-style mineralization have been previously mentioned for Vathi [16]. The Al 2 O 3 /TiO 2 and V/Sc ratios for unaltered and altered samples, respectively, can be applied to assess the magmatic fertility and prospectivity of Cu and Au [85]. Using the criteria of Loucks [85], samples of quartz monzonite and latite at Vathi were plotted as functions of Al2O3/TiO2 versus SiO2 and V/Sc versus SiO 2 (Figure 21a,b). The V/Sc ratio of samples from Vathi ranges from 10 to 50.91 and 11.84 to 17.08 for quartz monzonite and latite, respectively (Table 2), and they fall in the category of prospective Cu-Au deposits worldwide as characterized by V/Sc >10 [85]. In both Figure 21a,b, the Vathi samples plot in the fields of more productive Cu and Cu + Au intrusive rocks.
At the Vathi area, the ore-forming magma was focused at the intersection of major regional structures in the upper crust at a depth of >3 km where it crystallized to form stocks, plugs, and dikes of intermediate to felsic composition ( Figure 3; inset map). The major and minor structural features provided a vertically extensive dilatant network, which favored the evolution of the magmatic-hydrothermal porphyry system at Vathi (Figure 22a).
Disseminations, veins, veinlets, and small aggregates of magnetite + chalcopyrite + pyrite + bornite + molybdenite formed in the potassic alteration at Vathi. Volatile fluids migrated upward from the monzonitic magma through fractures. The initial high-temperature parental magmatic fluid cooled from 540 to 390 • C and decompressed at approximately 646 to 248 bars, which resulted in phase separation at a depth of 1 to 2.6 km (assuming lithostatic conditions). The resulting moderate-density vapor and saline brine fluids interacted with the surrounding host rocks to form potassic alteration, A-type veins, and associated Cu-Mo-Au mineralization (Figure 22a). In addition, the mixing of magmatic fluids and meteoric water during this stage resulted in boiling and phreatic brecciation. During the late stages of potassic alteration, propylitic alteration was developed, overprinting both the latite and the quartz monzonite ( Figure 22b).
Later, the vapor phase formed synchronously with the potassic alteration stage and continued to move upwards, cooling and mingling with meteoric water to form an acidic moderate-salinity (8-11 wt % NaCl equiv) aqueous fluid at lower temperatures (310-390 • C). Sericitic alteration developed from these fluids along brittle fault zones (Figure 22c). The secondary porosity of the latite was further enhanced by extended hydraulic fracturing favoring the formation of a network of cracks in the latite and the quartz monzonite. Within this network, the pyritic mineralization and D-veins associated with gold were developed.
During the late epithermal stage, a low-salinity fluid (1.4-2.9 wt % NaCl equiv) was generated at relatively low Th (205 and 259 • C). This fluid represents a mixture of a moderately saline fluid with circulating meteoric waters at a depth of 110 to 250 m under hydrostatic conditions, which formed the E-type veins and the attendant sphalerite + galena + arsenopyrite ± stibnite ± tennantite mineralization (Figure 22d). These minerals are common in the late stages of many porphyry systems and are the result of Zn, Pb, As, and Sb accumulation in the residual ore-forming solution, commonly as late-stage epithermal E-type veins [55,56].
The assemblage pyrite + pyrrhotite + chalcopyrite at Vathi is a rather uncommon paragenesis in porphyry systems. In reduced porphyry systems, the low solubility of S and Cu, the predominance of S 2-, and the depletion in Cl − favor the formation of massive pyrrhotite in temperatures between 350 and 600 • C [86][87][88]. However, pyrrhotite related to oxidized Cu-Au-porphyry systems has been reported from the biotite-actinolite alteration zones of alkaline igneous rock-related porphyry systems, peripheral to porphyry-Au deposit quartz stockwork veins, and they are caused by the interaction between late chlorite-illite-quartz alteration and carbonaceous wall rocks [56,88]. At Vathi, it is likely that the formation of pyrrhotite is due to a lower sulfur fugacity of the ore fluid. During the late epithermal stage, a low-salinity fluid (1.4-2.9 wt % NaCl equiv) was generated at relatively low Th (205 and 259 °C). This fluid represents a mixture of a moderately saline fluid with circulating meteoric waters at a depth of 110 to 250 m under hydrostatic conditions, which formed the E-type veins and the attendant sphalerite + galena + arsenopyrite ± stibnite ± tennantite mineralization (Figure 22d). These minerals are common in the late stages of many porphyry systems and are the result of Zn, Pb, As,  and Sb accumulation in the residual ore-forming solution, commonly as late-stage epithermal E-type veins [55,56].

Figure 22.
Schematic model of the formation of magmatic rocks, hydrothermal alteration, and ore mineralization at Vathi along an N-S cross-section: (a) Regional structural settings promote the fast ascent of the magmas (the position of faults is indicative) to form the quartz monzonite intrusions. Magmatic fluid exsolves from magma and flows upwards undergoing phase separation to form co-existing vapor and brine. As these fluids cool, they result in the potassic alteration and the related ore mineralization; (b) At the late stage of formation of the potassic alteration, magmatic fluids propylitize the latite and the quartz monzonite; (c) As the system cools, the vapor phase condenses and mixes with meteoric water, resulting in a late sericitic alteration that overprints potassic and propylitic alteration; (d) A new pulse of a low-salinity fluid of meteoric origin flows along the brittle fault zone to form the epithermal alteration and the late-stage metallic mineralization.
The assemblage pyrite + pyrrhotite + chalcopyrite at Vathi is a rather uncommon paragenesis in porphyry systems. In reduced porphyry systems, the low solubility of S and Cu, the predominance of S 2-, and the depletion in Clfavor the formation of massive pyrrhotite in temperatures between 350 and 600 °C [86][87][88]. However, pyrrhotite related to oxidized Cu-Au-porphyry systems has been reported from the biotite-actinolite alter- Figure 22. Schematic model of the formation of magmatic rocks, hydrothermal alteration, and ore mineralization at Vathi along an N-S cross-section: (a) Regional structural settings promote the fast ascent of the magmas (the position of faults is indicative) to form the quartz monzonite intrusions. Magmatic fluid exsolves from magma and flows upwards undergoing phase separation to form co-existing vapor and brine. As these fluids cool, they result in the potassic alteration and the related ore mineralization; (b) At the late stage of formation of the potassic alteration, magmatic fluids propylitize the latite and the quartz monzonite; (c) As the system cools, the vapor phase condenses and mixes with meteoric water, resulting in a late sericitic alteration that overprints potassic and propylitic alteration; (d) A new pulse of a low-salinity fluid of meteoric origin flows along the brittle fault zone to form the epithermal alteration and the late-stage metallic mineralization.

Summary and Conclusions
The Vathi porphyry Cu-Au ± Mo mineralization occurs in the northern part of the Oligocene-Miocene Serbo-Macedonian metallogenic province. The mineralization is mainly hosted in latite and is genetically associated with the intrusion of quartz monzonite (18 to 17 Ma). The quartz monzonite appears as fault-controlled, narrow, elongated stocks, while latite forms a plug-shaped body. These two rocks are silica-saturated and exhibit alkali-calcic to slightly alkaline chemical affinities. A phreatic breccia crosscuts the latite and was likely formed during the intrusion of the quartz monzonite. Three major fault groups trending E-W, NW-SE, and NE-SW characterize the regional structures, which controlled the emplacement of the intrusions in the crystalline basement of the Vertiskos Unit. The quartz monzonite and latite were affected by potassic and propylitic alteration that were subsequently overprinted by sericitic alteration. The potassic alteration is associated with the ore assemblage pyrite + chalcopyrite + bornite + molybdenite + magnetite in M-and A-type veins, and the propylitic alteration is related to pyrite + chalcopyrite, while the later sericitic alteration is associated with the assemblage pyrite + chalcopyrite + native gold ± tetradymite in D-type veins. The assemblage sphalerite + galena + arsenopyrite + pyrrhotite + pyrite ± stibnite ± tennantite is related to a subsequent epithermal overprinting event. The metallic mineralization is enriched in Cu (average 2313 ppm), Mo (average 152 ppm), Au (average 0.46 ppm), and U (average 31 ppm). A local enrichment in some REEs (up to 613 ppm La, 894 ppm Ce, 79 ppm Pr, 211 ppm Nd and 23 ppm Sm) is attributed to the phreatic breccia and to secondary supergene mineralization. Fluid inclusion data indicate that A-type porphyry style veins were deposited at 390->600 • C and at pressures up to 646 bars (<2.6 km depth) from boiling hydrothermal fluids. This process resulted in the generation of coexisting moderately saline (9.7-18.5 wt % NaCl equiv.) and highly saline (32.7-47.2 wt % NaCl equiv.) fluids. Sericitic alteration, D-type veins, and related pyrite + pyrrhotite + chalcopyrite mineralization are associated with a moderately saline fluid (8.4-11.2 wt % NaCl equiv) that formed as a result of condensation of vapor-rich fluids between 311 and 392 • C. At lower temperatures (205-259 • C), a less saline (1.4-2.9 wt % NaCl equiv.) fluid formed, which deposited late epithermal E-type veins and attendant mineralization. The anomalous enrichments in U and REE and the alkaline character of magmatism are possible indicators that Vathi shares several characteristics in common with alkaline igneous rock-related porphyry deposits.