www.geologia-croatica.hr ABSTRACT Characterization of various fluid parameters in magmatic-hydrothermal ore mineralizations is potentially essential for interpretation of the conditions of formation and therefore for mineral exploration. Fluid inclusions can provide a useful and promising tool in the re- search of the ore forming processes in these systems. This review focuses on the nature, composition and origin of magmatic-hydrothermal ore forming fluids involved in the for- mation of representative Tertiary ore deposits at the Rhodope metallogenic province in NE Greece. These deposits are spatially related to Tertiary magmatism in NE Greece. Case studies are presented here and include an intrusion-hosted sheeted vein system (Kavala), a Au-rich carbonate replacement and quartz-vein mineralization (Asimotrypes), mineralized veins in Eptadendro-Rachi and Thasos island (Kapsalina and Panagia), porphyry Cu-Mo- Re-Au deposits in Pagoni Rachi and Maronia and epithermal Au-Ag mineralizations in Perama and Loutros. Hydrothermal fluids rich in CO2 together with elevated Au and Te content are common and occur at the Kavala intrusion hosted sheeted vein system, at the Asimotrypes Au-rich carbonate replacement mineralization and at the Panagia (Thasos) vein system. We classify all these ore mineralizations as intrusion-related gold systems (IRGS). Transport and precipitation of metals including Au and Te is favoured when CO2 is present. Precipitation of the ore mineralization takes place due to the immiscibility of the carbonic and the aqueous fluids which have a magmatic origin with the contribution of meteoric water. Cooling of magmatic hydrothermal fluids and dilution with meteoric water is a common cause for ore mineral formation in the vein mineralizations of Eptadendro/ Rachi and Kapsalina Thasos. At the Pagoni Rachi and Maronia porphyry deposits, boiling and the high proportion of the vapour phase are the most essential fluid processes which affected ore formation. The epithermal veins overprinting the Pagoni Rachi and the Ma- ronia porphyry systems and the HS-IS epithermal system in Perama Hill and the IS epither- mal mineralization in Loutros are characterized by low to moderate temperatures and low to moderate salinities. Cooling and dilution of the ore fluids are the main process for gold precipitation. We conclude that the different fluid parameters and microthermometric data indicate a variety of fluid origin conditions and sources which can affect the strategy for exploration and prospecting for gold, rare and critical metals. Keywords: Fluid inclusions, metallogeny, Rhodope, intrusion related systems, porphyry systems, epithermal mineralizations Article history: Received December 02, 2015 Revised and accepted January 01, 2016 Avaliable online February 29, 2016 1. INTRODUCTION Large-scale magmatic-hydrothermal fluid circulation in the Earth’s crust is derived from intermediate to felsic hydrous magmas mainly at convergent plate margins, providing endo- genic heat and mass transfer for the formation of ore deposits including base, precious and rare metals (LEHMANN et al., 2000; BLUNDELL et al., 2005; SIMMONS & BROWN, 2006; BORTNIKOV, 2006; HEINRICH, 2007). These sys- tems involve large amounts of volatiles including magmatic, meteoric or metamorphic water in aqueous or vapour form, CO2, NaCl, KCl, CaCl2, MgCl2, H2S, ±CH4, ±N2, as well as metals which originate from the magmas or the leaching pro- cesses from the rocks penetrated by the fluids through fluid- wall-rock interactions. The detection of the various fluid characteristics is poten- tially important for interpretation of the formation of numer- ous ore deposit types and consequently for mineral exploration (BAKER & LANG, 2001; WILKINSON, 2001; HEINRICH, 2005; CHANG et al., 2011; KESLER et al., 2013). For exam- ple boiling has been a valuable guide in gold and copper explo- ration in porphyry systems. It is useful therefore to know how the volatile phases originate, how they dissolve and transport metals, which is the role of temperature and pressure and un- der what conditions the metallic minerals are deposited. A useful fluid parameter which can be achieved by the fluid inclusions study, especially in the systems associated with Au and other sidirophile elements enrichment, is the CO2 content, because there is a strong correlation between CO2 content and gold mineralization, such as in the deep- seated intrusion-related systems (LOWENSTERN, 2001; BAKER 2002). In general, fluid inclusions with elevated ho- mogenization temperatures and high salinities indicate mag- matic productive fluids, while fluid inclusions with low ho- 157-167 3 Figs. 1 Tab. doi: 10.4154/gc.2016.12 Fluid evolution in Tertiary magmatic-hydrothermal ore systems at the Rhodope metallogenic province, NE Greece. A review. Vasilios Melfos1,* and Panagiotis Voudouris2 1 Aristotle University of Thessaloniki, Department of Mineralogy, Petrology and Economic Geology, Thessaloniki, 54124, Greece (*corresponding author, email: melfosv@geo.auth.gr) 2 University of Athens, Department of Mineralogy-Petrology, Athens 15784, Greece 158 Vol. 69/1 mogenization temperature and low salinity show incorpora- tion of meteoric waters (WILKINSON, 2001; KESLER et al., 2013). It is evident that homogenization temperatures, fluid salinity, density and pressure, can be useful in the deter- mination of the fluid sources, the physical and chemical pro- cesses of transport and deposition of metals and also to ex- plain the chemical environments of ore deposition and the possible zonation patterns of some ore systems. This review paper summarizes the nature, composition and origin of magmatic-hydrothermal ore forming fluids which were involved in the formation of representative Tertiary ore deposits at the Rhodope metallogenic province in NE Greece (Fig. 1) which are developed within or around magmatic intru- sive or extrusive rocks. We explore the evolution of high tem- perature and high salinity magmatic brines or CO2-rich fluids, as well high temperature and low salinity magmatic vapours, to low temperature and low salinity meteoric fluids which were involved in the formation of these ore mineralizations. This paper is partly motivated by microthermometric mea- surements of fluid inclusions which have been published in previous works of several researchers (VAVELIDIS et al., 1995; ELIOPOULOS, 2000; MELFOS et al., 2002; VAVELI- DIS & MELFOS, 2004; MELFOS et al., 2008, 2010; FORNA- DEL et al., 2011; ELIOPOULOS & KILIAS, 2011; VOU- DOURIS et al., 2011, 2013b) and include a reduced intrusion- related system (Kavala), a Au-rich carbonate replacement mineralization (Asimotrypes), mineralized veins in Epta- dendro-Rachi and Thasos island (Kapsalina and Panagia), por- phyry Cu±Au±Mo deposits in Pagoni Rachi and Maronia and epithermal mineralizations in Perama and Loutros (Fig. 1). The purpose of this work is to show that detailed exploration for the discovery of new ore deposits requires a thorough fluid inclusion study which will interpret fluid phase relations and evaluate fluid characteristics in view of the ore forming conditions. 2. GEOLOGICAL SETTING OF TERTIARY ORE DEPOSITS Greece is traditionally a mining country with numerous ore deposits which include base (Cr, Ni, Mn, Al, Fe, Pb, Zn, Cu), precious (Au, Ag, Platinum Group Elements-PGE), rare (Bi, Te, Mo, Re, Ti) & critical (Mg, Sb, W, Co, Ga, Ge, In, REE) metals (MELFOS & VOUDOURIS, 2012). Mining dates back to prehistoric times focusing initially on gold, copper and iron mineralizations. Later, in the Classical and Hellenistic times, ore exploitation became more intense especially for Au, Ag, Cu and Fe. Famous ancient mining centres operated in Milos and Sifnos Islands (Aegean Sea), in Lavrion (Attica), in Chalkidiki, Palea Kavala, Pangeon Mt and Thasos Island (N. Greece), where a great number of underground galleries and metallurgical furnaces are still preserved (PERNICKA et al., 1981; GEORGAKOPOULOU et al., 2011). Several metals, such as Ag, Au, Fe, Cu, were also produced during the Byzan- tine and Ottoman periods (NERANTZIS, 2009). Since the be- ginning of the 20th century a great number of ore deposits have been targeted for exploration and exploitation of Cr, Ni, Al, Mn, Fe, Pb, Zn, Cu, Au and Ag. Greece is the most important producing country in the EU for Ni-Fe (from laterites) and Al Figure 1. Simplified geological map of the Rhodope metallogenic province, NE Greece, showing the studied mineralizations. Kavala (Kv) intrusion-related system, Asimotrypes (As) Au-rich carbonate replacement mineralization, Eptadendro-Rachi (ER) mineralized veins, Kapsalina (Kp) and Panagia (Pn) Cu-Bi-Ag±Te vein systems in Thasos island, Pagoni Rachi (PR) and Maronia (Mr) Cu-Mo-Re-Au porphyry deposits, Perama Hill (PH) and Loutros (Lt) epithermal mineralizations. 159 Melfos & Voudouris: Fluid evolution in Tertiary magmatic-hydrothermal ore systems at the Rhodope metallogenic province... (from bauxites). The production of Mg from magnesite is also significant (TSIRAMBIDES & FILIPPIDIS, 2012). The geological setting of the Hellenides, part of the Alpine- Himalayan orogeny, comprises numerous geotectonic terrains, which resulted as a consequence of the complicated geody- namic events in the region of the Palaeotethys and Neotethys Oceans, including ocean rifting, spreading, subduction, colli- sion, post-collision and extensive subduction and post-subduc- tion magmatism (PE-PIPER & PIPER, 2002). After closure of the Neotethys Ocean, syn-orogenic exhumation of HP-LT rocks occurred during the Late Cretaceous-Eocene, before an acceleration of slab retreat changed the subduction regime and caused collapse of the Hellenic mountain belt and thinning of the Aegean Sea from Oligocene to the present (JOLIVET & BRUN, 2010). During this post-orogenic episode large scale detachments formed which exhumed metamorphic core com- plexes in a back-arc setting. Tertiary to Quaternary calc-alka- line to alkaline magmatism in the Aegean region occurred mostly in a post-collisional setting behind the active Hellenic subduction zone (PE-PIPER & PIPER, 2001). Lithospheric delamination, roll-back and/or break-off of the subducting slab were the principal mechanisms resulting in magma generation, which was triggered by asthenospheric mantle upwelling (JO- LIVET & BRUN, 2010). In this geotectonic regime, large ore deposits of economic importance were formed during the Tertiary, in relation to magmatism and the evolution of extensive magmatic-hydro- thermal systems. The ore types include reduced intrusion-re- lated, skarn, porphyry, high-, intermediate- and low-sulfidation epithermal, carbonate-hosted replacement Pb-Zn and polyme- tallic vein deposits (MELFOS et al., 2002; VOUDOURIS & ALFIERIS, 2005; VOUDOURIS, 2006; BONSALL et al., 2011; MELFOS & VOUDOURIS, 2012; VOUDOURIS et al., 2013a). They are concentrated into four regional metallogenic regions: 1) the Rhodope metallogenic province (RMP), 2) the Serbomacedonian metallogenic province (SMMP), 3) the At- ticocycladic metallogenic province (ACMP) and 4) the East- ern Aegean metallogenic province (EAMP). These regions are the most promising targets for the near future exploration and exploitation of precious, rare and critical metals in Greece. Tertiary basins and magmatic settings in the Rhodope mas- sif are associated with significant skarn, polymetallic vein, porphyry and epithermal ore deposits, which are occasionally rich in Au and critical metals (MELFOS et al., 2002; VOU- DOURIS, 2006; VOUDOURIS et al., 2009, 2011, 2013b), which constitute a part of the Rhodope metallogenic province. The magmatic-hydrothermal deposits of the Rhodope meta- llogenic province were formed during the final stage of the Tertiary orogenic collapse, and are related to the emplacement of hot asthenosphere into shallow crustal levels above a de- tached lithospheric plate which generated widespread magma- tism (DE BOORDER et al., 1998; BLUNDELL et al., 2005). The Oligocene-Miocene magmas were intruded to shallow depths along deep seated detachment faults, forming mafic to felsic mantle-derived plutonic, subvolcanic and volcanic rocks of calc-alkaline, high-K calc-alkaline and shoshonitic to ultra- potassic (lamprophyric) affinities (CHRISTOFIDES, 1996; MARCHEV et al., 2005). 3. FLUID INCLUSION CASE STUDIES 3.1. KAVALA INTRUSION-RELATED GOLD SYSTEM The Kavala intrusion-related gold system (IRGS) is part of the Palea Kavala ore system (Fig. 1) which includes ~150 ore mine ralizations occurring either as primary hypogene quartz- calcite-sulfide veins or as supergene oxidized ore bodies, asso- ciated with the ~21–22 Ma granodioritic Kavala pluton of I- type (FORNADEL et al., 2011). This pluton intruded marbles, schists, amphibolites and gneisses of the Palaeozoic Rhodope metamorphic core complex, along the regional E–W trending Kavala-Xanthi-Komotini fault (Fig. 1). The mineralization is considered to be part of an intrusion-related ore system in overlapping zones centred on the pluton. The marble hosts dis- tal oxidized Fe-Mn (Pb±Zn±Ag) and Fe-Mn-Au ore bodies, whereas the gneiss and pluton contain proximal Fe-Cu-Au, Fe- As-Au and Bi-Te-Au vein mineralization. These ore minera- lizations were exploited in ancient times for gold and silver with extensive underground and surface mining operations (VAVELIDIS et al., 1996). Figure 2. Fluid inclusion types from Tertiary ore mineralizations in the Rhodope metallogenic province, NE Greece. a. Type L two phase aque- ous liquid-vapour fluid inclusion coexisting with type C three phase inclusion which contain H2O-CO2, in Kavala sheeted vein system in the pluton; b. Type L two phase, liquid and vapour, fluid inclusion in a quartz vein at Eptadendro/Rachi; c. Type V vapour-rich two-phase aqueous inclusion coexisting with type L-S inclusion which contain three phases: a liquid, a vapour and a halite crystal in Pagoni Rachi por- phyry system, indicating that they originate from a deep near-critical single phase fluid by phase separation; d. Type V two-phase vapour- rich inclusion with an opaque daughter crystal possibly chalcopyrite in Pagoni Rachi porphyry system; e. Type V vapour-rich two-phase aque- ous inclusion coexisting with type L-S inclusion which contain three phases: a liquid, a vapour and a halite crystal, and type L two phase aqueous liquid-vapour fluid inclusion in Maronia porphyry system; f. Type L aqueous liquid rich fluid inclusions which contain mainly NaCl in quartz at Perama Hill epithermal deposit. 160 Vol. 69/1 The SE-NW trending Bi-Te±Pb±Sb±Au quartz sheeted veins crosscut the pluton and the neighbouring metamorphic rocks in the Chalkero area, and contain pyrite, tetradymite, bis- muthinite, cosalite and Sb-lillianite (FORNADEL et al., 2011). The veins consist of quartz, K-feldspar, albite and muscovite, altered to sericite and kaolinite. Fluid inclusions in these veins (Tab. 1) are rounded in shape, consistently distributed through- out the samples with a length of up to 60 μm (with the majo rity in the 25 to 35 μm range). They are classified into three types which coexist in each sample (MELFOS et al., 2008; FORNA- DEL et al., 2011) and therefore belong to the same primary fluid inclusion assemblage based on the criteria of GOLD- STEIN & REYNOLDS (1994). Type L contains two phase aqueous liquid-vapour fluid in- clusions (Fig. 2a) with NaCl and CaCl2 and 20-30 volume % vapour. Homogenization temperatures range between 216° and 420°C, with a maximum at 380°C, and the salinities vary from 15.9 to 22.6 wt% NaCl equiv. Type C inclusions (Fig. 2a) are three phase: aqueous liquid, CO2-liquid, CO2-vapour ap- proximately 50-70 volume %, with minor CH4. They homoge- nize at 255°-414°C, with a peak at 320°C, demonstrating sa- linities from 5.5 to 11.2 wt% NaCl equiv. The third type V consists of two phase aqueous vapour-liquid fluid inclusions with a bubble of 90-95%, which homogenized to the vapour phase at Th between 210° and 323°C (max= 280°C). The three fluid inclusion types are considered to represent coexisting immiscible carbonic and aqueous ore forming fluids (Fig. 3a). Isochores for fluid inclusions within the Kavala veins revealed a depth of formation, between 3 and 5 km (MELFOS et al., 2008; FORNADEL et al., 2011). Carbon dioxide played a significant role in the ore mineral precipitation in the Kavala veins. 3.2. ASIMOTRYPES CARBONATEREPLACEMENT AND VEIN TYPE Au DEPOSIT AT PANGEON MOUNTAIN The Pangeon mountain includes numerous carbonate replace- ment, vein and shear-hosted Au-Ag ore mineralizations, which were mined intensively during ancient times for gold and sil- ver (VAXEVANOPOULOS et al., in press). They are hosted in the marbles and gneisses of the Rhodope massif as well as in the Pangeon granitoids. One of the most important ore mineralizations occurs in Asimotrypes (Fig. 1) which is hosted in the Nikisiani pluton, part of the Pangeon granitoids, and in the marbles of the Rho- dope massif (ELIOPOULOS, 2000; ELIOPOULOS & KIL- Figure 3. Homogenization temperatures versus salinity plots of fluid inclusions from various Tertiary ore mineralizations in the the Rhodope metallogenic province, NE Greece. a. Kavala intrusion related sheeted quartz veins (n=241); b. Eptadendro/Rachi mineralized quartz veins (n=249); c. Pagoni Rachi porphyry sys- tem overprinted by epithermal veins (n=362); d. Maronia porphyry system overprinted by epithermal veins (n=90); e. Perama Hill HS-IS epithermal system (n=141); f. Loutros IS epithermal system (n=36). n=number of measurements. 161 Melfos & Voudouris: Fluid evolution in Tertiary magmatic-hydrothermal ore systems at the Rhodope metallogenic province... IAS, 2011; VAXEVANOPOULOS et al., in press). The Pan- geon granitoids consist of tonalite, granodiorite and granite. They demonstrate characteristics of high-K calc-alkaline met- aluminous to slightly peraluminous post-collisional rocks and have a hornblende 40Ar/39Ar age of 21 to 22 Ma (ELEFTHE- RIADIS & KORONEOS, 2003). The mineralization is associated with structurally controlled quartz veins, pods and lenses, and in gold-bearing replacement bodies, and is partly oxidized. It consists mainly of pyrite and arsenopyrite, with minor galena, sphalerite, chalcopyrite, pyr- rhotite, tetrahedrite-tennantite, marcasite and gold (ELIO- POULOS, 2000; VAXEVANOPOULOS et al., in press). Fluid inclusions are irregular to euhedral in shape with smooth boundaries reaching in size up to 20 μm. Three types of primary fluid inclusions (Tab. 1) were distin- guished by ELIOPOULOS (2000) and coexist in each sample being in equilibrium with the sulfide minerals. Type L contains two phase aqueous liquid-vapour fluid inclusions with NaCl. The bubble occupies 10-20 volume % and the homogenization temperatures range from 135° to 256°C with a peak at 190°C and salinities from 5.62 to 9.05 wt% NaCl equiv. Type C inclu- sions are three phase and contain liquid H2O-rich, liquid CO2 and vapour CO2 (70 to 90% by volume). They homogenize at temperatures between 275° and 340°C, with a maximum at 320°C, and show salinities from 2.61 to 7.31 wt% NaCl equiv. The third type V consists of two phase aqueous vapour-liquid fluid inclusions which show decrepitation or leaking and no microthermometric measurements were achieved. The gold mineralization was formed by coexisting immis- cible carbonic and aqueous fluids at pressures of 1.8 to 2.0 kbars under lithostatic conditions corresponding to depths of 6 to 9 km with the ore fluids at ~270°C (ELIOPOULOS & KIL- IAS, 2011). 3.3. EPTADENDRO-RACHI VEIN TYPE Ag-Au ORE MINERALIZATION The Cu-rich sulfide mineralization at the Eptadendro-Rachi ore district in Thrace, NE Greece, occurs at the contacts be- tween granitoid intrusions (trodhjemites and pegmatites) of Upper Cretaceous to Early Tertiary age, and the meta-ultraba- sic to basic rocks of the eastern Rhodope massif (Fig. 1). Ac- cording to CHATZIKIRKOU & MICHAILIDIS (2004) the ore bodies form quartz veins with pyrite, chalcopyrite, sphalerite, galena, hessite, bismuthinite, emplectite, tetradymite, aikinite, wittichenite, siegenite, millerite, bornite, pyrrhotite, covellite, magnetite, haematite and goethite. Chlorite, calcite and seric- ite are the alteration minerals associated with the vein minera- lization. Fluid inclusions have negative crystal or rounded to elon- gated isometric shapes with a length up to 42 μm, and are dis- tributed along crystal faces of quartz and therefore are sug- gested to have a primary origin (MELFOS et al., 2010). At room temperature they contain two phases, a liquid and a va- pour phase which occupies 20-30 % by volume. Two fluid in- clusion assemblages (FIAs) were identified MELFOS et al. (2010). The first FIA is hosted in quartz and is characterized by two phase, liquid and vapour, inclusions (Type L, Fig. 2b) at room temperature which homogenize to liquid at temperatures from 300° to 380°C, with a peak at 330°C, and exhibit low salinities, 1.6 to 7.2 wt% NaCl equiv, corresponding to the flu- ids of the main ore stage. The second FIA also include two phase fluid inclusions (Type L) and are distinguished by a drop in Th (210° to 260°C, with a peak at 240°C) corresponding to the late ore stage as- sociated with calcite formation, whereas the salinities (3.2 to 6.3 wt% NaCl equiv) are similar to the first FIA. NaCl is the dominant salt component of the fluids. The low to moderate salinities (1.6 to 7.2 wt% NaCl equiv) of the two FIAs are possibly attributed to mixing between as- cending magmatic fluids and downward migrating dilute mete- oric water. According to MELFOS et al. (2010) such a process is very likely in a hydrothermal environment with a fault-con- trolled vein structure and causes ore deposition by diluted flu- ids at decreasing temperatures, from 380° to 210°C (Fig. 3b). 3.4. THASOS (KAPSALINA AND PANAGIA) VEIN TYPE ORE MINERALIZATIONS Thasos has been mined for Au and Ag both during antiquity, and in the 20th century (VAVELIDIS et al., 1988). Silver rich carbonate Pb-Zn replacement deposits occur at the western part of the island, whereas Cu-Au and gold deposits exist at the northeastern and eastern part of the island (VAVELIDIS & AMSTUTZ, 1983). Two Cu-Bi-Ag±Te mineralizations on Thasos island, in Kapsalina and Panagia areas (Fig. 1) have been described by VAVELIDIS et al. (1995) and VAVELIDIS & MELFOS (2004). The Kapsalina copper mineralization is hosted in quartz veins which crosscut the schists of the Rhodope massif. The ore mineralization consists mainly of chalcopyrite and pyrite, with minor bismuthinite, tetrahedrite-tennantite and gersdorf- fite (VAVELIDIS & MELFOS, 2004). Fluid inclusions were studied in quartz from the mineralized veins and were identi- fied as primary and pseudosecondary. Their size ranges from 4 to 45 μm, mainly up to 10 μm. Based on their phase assem- blages at room temperature one type of inclusion (Type L) was distinguished containing a liquid aqueous solution and a va- pour bubble (approximately 10 volume %). The fluid is domi- nated by NaCl with a minor contribution of KCl demonstrating salinities from 1.0 to 3.3 wt% NaCl equiv. Homogenization temperatures of the fluid inclusions vary mainly from 120 to 280°C, with a peak at 200°C (VAVELIDIS & MELFOS, 2004). The mineralization was formed by cooling of the hy- drothermal fluids which were derived mainly by meteoric wa- ters. The mineralization in Panagia is hosted in quartz veins which crosscut the schists, marbles, dolomitic marbles, am- phibolites and gneisses of the Rhodope massif. The minera- lization consists of chalcopyrite, pyrrhotite, pyrite, tetrahe- drite-tennantite, gersdorffite, tetradymite and hessite. Microthermometric measurements were carried out on fluid inclusions of the mineralized quartz (VAVELIDIS et al., 1995). Based on their phase proportions at room temperature, three types of fluid inclusions were distinguished which belong to the same primary fluid inclusion assemblage based on the cri- teria of GOLDSTEIN & REYNOLDS (1994). Fluid inclusions of type L (up to 18 μm in length) consist of liquid rich two- 162 Vol. 69/1 phase inclusions which contained mainly NaCl and the vapour ranges from 10 to 20 % in volume. The homogenization tem- peratures vary between 229° and 305°C, with a maximum at 290°C and the salinities are low, from 3.12 to 3.28 wt% NaCl equiv. The type C fluid inclusions, with a length up to 20 μm, con- tain three phases dominated by aqueous liquid, CO2-liquid and CO2-vapour approximately 60-70 % of the total volume. A to- tal fluid salinity of 3.7 to 6.5 wt% NaCl equiv was estimated. Final homogenization temperatures of this type of inclusions vary from 313 and 360°C, with a maximum at 330°C. A third type of fluid inclusion is made up by large irregularly shaped cavities, which contain pure CO2. The vein type mineraliza- tions in Panagia were formed by two coexisting immiscible carbonic and aqueous fluids and CO2 played a significant role in the ore precipitation. 3.5. PAGONI RACHI PORPHYRY Cu-Mo-Re-Au DEPOSIT The Pagoni Rachi Cu-Mo-Re-Au mineralization (Fig. 1) be- longs to the Kirki mining district and is a telescoped porphy- ry-epithermal system which is spatially related to calc-alka- line granodiorite-tonalite porphyry stock of Oligocene age (VOUDOURIS et al., 2013b). Alteration types and vein rela- tionships suggest that mineralization can be divided into four paragenetic stages, which are, from early to late: (1) sodic/ potassic-calcic alteration with quartz- and magnetite-bearing veins (A- and M-type) and distal propylitic alteration, (2) so- dic/potassic alteration with quartz-pyrite-chalcopyrite-mo- lybdenite veins (B-type), (3) sericitic alteration with “transi- tional” porphyry to epithermal pyrite-chalcopyrite-molybde- nite veins (D-type), and (4) argillic alteration with quartz- calcite base metal and precious metal rich veins (E-type) with epithermal affinity. The mineralization contains molybdenite with up to 4.7 wt% Re and the rare mineral rheniite (VOU- DOURIS et al., 2009, 2013b). Fluid inclusions were studied in the various vein types (A, B, and E) and are in equilibrium with the sulfide minerals and therefore represent the ore forming fluids (VOUDOURIS et al., 2013b). Their maximum size is 30 μm (mainly between 7 and 15 μm). They are isolated or arranged along quartz and calcite growth zones and are considered to have a primary origin based on the criteria of BODNAR (2003). Three types of primary fluid inclusions were identified (VOUDOURIS et al., 2013b). Type L consists of two-phase aqueous inclusions with a vapour bubble occupying 10 to 20 volume %, which homogenize to a liquid phase. Type L-S inclusions contain three phases: a liquid, a vapour (10 to 20 volume %) and a halite crystal (Fig. 2c). Type V inclusions are two-phase vapour-rich (~80 volume %) aqueous inclu- sions (Fig. 2c,d). They homogenize to the vapour phase. They commonly contain an accidental-trapped opaque crys- tal (Fig. 2d) possibly chalcopyrite (VOUDOURIS et al., 2013b). Three distinct primary fluid inclusion assemblages (FIAs) were identified (VOUDOURIS et al., 2009, 2013b). The first two assemblages contain coexisting L, L-S and V inclusion types in A- and B-veins. The third assemblage in- cludes late type L inclusions in the epithermal E-veins. A- and B-veins were deposited at temperatures ranging between 360 and 510 °C and at pressures up to 690 bars (<2 km depth) for A-veins and up to 360 bars (<1.5 km depth) for B-veins, from boiling hydrothermal fluids. This process pro- duced coexisting low to moderately saline (1.7-10.7 wt% NaCl equiv) and highly saline (36->74 wt% NaCl equiv) flu- ids (Fig. 3c). The fluid inclusions in the A- and B-type veins indicate that vapour was the predominant ore-forming fluid. Subsequent dilution of the moderately saline fluid resulted in a lower temperature (210–340 °C), less saline (1.4–2.9 wt% NaCl equiv) mineralizing fluid in late epithermal E-veins (Fig. 3c, VOUDOURIS et al., 2013b). 3.6. MARONIA PORPHYRY Cu-Mo-Re-Au MINERA LIZATION The Maronia Cu-Mo-Re-Au mineralization (Fig. 1) is hosted by a porphyritic microgranite (MELFOS et al., 2002). Three hydrothermal alteration zones have been recognized: an ar- gillic, a sericitic and a propylitic zone. Additionally, three highly silicified zones crop out at the microgranite. Chalco- pyrite-pyrite-molybdenite mineralization is located mainly in the silicified zones and is associated with areas of sericitic and propylitic alteration. The molybdenite revealed unusual- ly high and variable rhenium concentrations (Re from 0.12 to 4.21 wt%). Fluid inclusions which were studied in the ore-related quartz (MELFOS et al., 2002) have rounded to elongated or negative-crystal shapes and are isolated or arranged along quartz demonstrating a primary origin. Their size reaches up to 60 μm. Four types of fluid inclusions were identified (MELFOS et al., 2002). Type L contains 2-phase aqueous liquid-rich inclusions (vapour 20-35 volume %) and type L-S consists of 3-phases, an aqueous liquid phase, a vapour (20- 40 volume %) and a halite (NaCl) crystal (Fig. 2e). Type V inclusions are two-phase vapour-rich (70-90 volume %), whereas type V-S are 3-phase aqueous vapour-rich (60-80 volume %) inclusions with a daughter mineral (halite). Fluid inclusions, mainly types L and L-S contain an accidentally- trapped opaque phase (pyrite and/or chalcopyrite?), which suggests that they represent the mineralizing fluid at the Ma- ronia deposit. Microthermometric results from the porphyry system showed that propylitic alteration and associated magnetite- pyrite-quartz veins are characterized by temperatures be- tween 320 to 410°C, with a mean at 400°C. Trapping tem- peratures in the quartz-pyrite-chalcopyrite veins from seri- citic alteration zone range from 360 to 420°C, with a peak at 380°C, and record the main temperature range of copper de- position under boiling conditions. Pyrite-molybdenite veins in sericitic zone were formed at slightly lower temperatures, from 250 to 360°C, with a maximum at 340°C. The ore-bear- ing fluids were dominated by NaCl, with salinities ranging from 7.0 to 16.4 wt% NaCl equiv in type L inclusions, from 28 to 55 wt% NaCl equiv in L-S fluid inclusions, from 5.5 to 14.0 wt% NaCl equiv in type V inclusions and from 38 to 51 wt% NaCl equiv in V-S inclusions (Fig. 3d). Trapping pres- sures of the ore-forming fluids were estimated to vary from 150 to 510 bar and boiling along with cooling of the fluids are 163 Melfos & Voudouris: Fluid evolution in Tertiary magmatic-hydrothermal ore systems at the Rhodope metallogenic province... considered to be the main processes of ore formation (MEL- FOS et al., 2002). 3.7. PERAMA HILL HS-IS EPITHERMAL Au-Ag DEPOSIT An epithermal mineralization rich in Au, Ag, Te and Se oc- curs in Perama Hill at the Petrota graben, NE Greece (Fig. 1). The ore mineralization is hosted in andesitic rocks and the overlying conglomerates and sandstones and exhibits charac- teristics of a high sulfidation system which was followed by deposition of an intermediate sulfidation mineralization (VOUDOURIS et al., 2011). It forms a sulfide rich vein-type in the andesite and a Au-rich stratabound oxidized type in the overlying sedimentary rocks. The ore mineral paragenesis in- cludes a complex Pb-Zn-Bi-Se-Te-Sn-Au-Ag hypogene ore assemblage with native gold, native tellurium and electrum. Fluid inclusions have negative-crystal or elongated iso- metric shapes with a maximum diameter of 30 μm, and were observed in quartz (Fig. 2f) and barite from the high sulfida- tion veins in andesite at depths of 88 m and 97 m (VOU- DOURIS et al., 2011). They are aqueous liquid rich (Type L) with ~20 volume % vapour and contain mainly NaCl. They are isolated or arranged along the quartz or barite growth zones and consequently are defined as primary. Homogenization temperatures in quartz vary from 211 to 331°C, with a maximum at 240°C, and salinities from 0.7 to 4.9 wt% NaCl equiv. In barite the inclusions display Th val- ues of 222 to 256°C, with a maximum at 240°C, overlapping those of quartz, and very low salinities, between 0.8 and 2.3 wt% NaCl equiv. The mineralization was formed by cooling and dilution of the ore fluids which were possibly derived by condensation of magmatic vapour with limited mixing with meteoric water between a higher temperature (330°C) and a higher salinity (4.9 wt% NaCl equiv) end-member and a low- er temperature (210°C) dilute (0.7 wt% NaCl equiv) compo- nent (Fig. 3e, VOUDOURIS et al., 2011). 3.8. LOUTROS HILL IS EPITHERMAL Ag DEPOSIT The Loutros area belongs at the Feres-Soufli-Dadia volcano- sedimentary basin. An intermediate sulfidation epithermal mineralization (Fig. 1) is associated with a zeolite-altered rhyolitic lava dome which demonstrates an age of 19.5 Ma (whole rock K/Ar geochronology, CHRISTOFIDES et al., 2004). The ore mineralization consists mainly of pyrite and marcasite and occurs either in the form of hydrothermal brec- cias or in NW-trending massive or barite bearing veins. Rare- ly barite-galena veins with sphalerite and minor chalcopyrite crosscut the rhyolite. The mineralized barite veins are en- riched in Mo (up to ≤15.6 g/t), As (up to 1058 g/t) and Ag (up to 31 g/t). Fluid inclusions with irregular elongated shapes and a size up to 25 μm, were observed in barite and are represented by type L which contains an aqueous liquid rich phase (10-20 volume % vapour). Microthermometric results revealed a ho- mogenization temperature range mainly between 220° and 329°C, with a maximum at 310°C, and salinities from 8.65 to 16.14 wt% NaCl equiv (unpublished data). The formation of the mineralization is attributed to cooling of an ore magmatic fluid, dominated by NaCl, and a possible limited mixing with meteoric water at 330 to 220°C (Fig. 3f). 4. DISCUSSION Based on phase composition at room temperature, five types of fluid inclusions were identified in quartz, calcite and barite from magmatic-hydrothermal ore mineralizations in the Rho- dope metallogenic province in NE Greece: (1) Type L - two- phase fluid inclusions consisting of an aqueous (low-salinity) solution with a vapour bubble with NaCl and sometimes minor CaCl2; (2) Type C - three phase inclusions which contain H2O- CO2, sometimes with occasional CH4; (3) Type V - two phase aqueous vapour-liquid fluid inclusions; (4) Type L-S - multi- phase phase inclusions which contain an aqueous solution, a vapour bubble and a halite crystal; (5) Type V-S – three phase aqueous vapour-rich inclusions with a halite crystal. At the porphyry mineralizations L and L-S fluid inclusions contain a pyrite or chalcopyrite crystal which shows that they represent the mineralizing fluid (ROEDDER, 1977; BORTNIKOV, 2006). However not all types of fluid inclusions occur in all ore mineralizations, suggesting that they formed from fluids of dif- ferent composition and under various conditions (Table 1). Fluid inclusions in the Kavala reduced intrusion related sheeted vein system hosting a Bi-Te-Pb-Sb-Au mineralization, display a range of fluid compositions (Types L, C). Inclusions with low to moderate salinities (5.5 to 11.2 wt% NaCl equiv) CO2-rich fluids are common and co-exist with moderate to high salinity fluids (salinities 15.9 to 22.6 wt% NaCl equiv). Methane is detectable in the CO2-rich fluid inclusions. Kavala veins are characterized as a deep seated system (~3 to 5 km) which was formed by two co-existing high temperature (Th=300° and 400°C), moderate to high salinity, immiscible carbonic and aqueous fluids (FORNADEL et al., 2011). Fluid inclusions from Asimotrypes carbonate replacement mineralization at Pangaion mountain demonstrate the co-exis- tence of two immiscible fluids, a low to moderate (2.61 to 7.31 wt% NaCl equiv) carbonic fluid and a moderate salinity (5.62 to 9.05 wt% NaCl equiv) aqueous fluid. ELIOPOULOS & KILIAS (2011) have shown that the main cause of ore precipi- tation should be attributed to unmixing of two fluids, a car- bonic and an aqueous, at depths from 6 to 9 km at ~270°C. In the quartz veins at Panagia of Thasos the fluids are dominated by two immiscible fluids, an aqueous-carbonic and an aqueous fluid, which have low to moderate salinities (3.7 to 6.5 wt% eq. NaCl) and trapping temperatures from 340° to 360°C at 500- 700 kbar (VAVELIDIS et al., 1995). The presence of Te and the elevated temperatures of the Panagia vein system possibly show that it is located proximally to a magmatic rock, which however is not exposed on Thasos island but can only be as- sumed. Kavala, Asimotrypes and Panagia ore systems show many similarities regarding their fluid characteristics. According to the intrusion related model presented by HART (2005) the Ka- vala sheeted veins in the pluton and the adjacent metamorphic rocks represent the most proximal system, whereas Asimo- trypes carbonate replacement mineralization occurs at some distance from the Nikisiani pluton, being part of an intrusion- related gold system (IRGS). Panagia vein mineralization at 164 Vol. 69/1 Ta bl e 1. T er tia ry o re m in er al iz at io ns o f t he R ho do pe m et al lo ge ni c pr ov in ce , N E G re ec e, a nd th ei r fl ui d in cl us io n ch ar ac te ris tic s (ty pe s, te m pe ra tu re s, sa lin iti es , c om po si tio n) . L : l iq ui d ric h aq ue ou s flu id in cl us io ns ; C : t hr ee p ha se fl ui d in cl us io ns w ith a n aq ue ou s l iq ui d, a C O 2-l iq ui d an d a C O 2-v ap ou r p ha se ; V : v ap ou r r ic h aq ue ou s fl ui d in cl us io ns ; L -S : l iq ui d ric h aq ue ou s fl ui d in cl us io ns w ith a da ug ht er so lid p ha se ; L -S : v ap ou r r ic h aq ue ou s fl ui d in cl us io ns w ith a d au gh te r s ol id p ha se . O re d is tr ic t M et al s O re st yl e H os t m in er al FI ty pe s T h° C (m ax ) Sa lin ity w t% N aC l e qu iv Fl ui d co m po si - tio n O re fo rm at io n co nd iti on s R ef er en ce s K av al a Pb , S b, B i, Te , A u In tru si on -h os te d, sh ee te d ve in s Q ua rtz L 21 6- 4 20 (3 80 ) 15 .9 - 22 .6 N aC l a nd C aC l 2 C oe xi st in g im m is ci bl e ca rb on ic a nd a qu e- ou s fl ui ds a t a d ep th o f 3 -5 k m . C O 2 p la ye d a si gn ifi ca nt ro le in th e or e pr ec ip ita tio n M EL FO S et a l. (2 00 8) ; F O R N A D EL et a l. (2 01 1) C 25 5- 41 4 (3 20 ) 5. 5 - 1 1. 2 H 2O -C O 2, m in or C H 4 V 21 0- 32 3 (2 80 ) - N aC l a nd C aC l 2 A si m ot ry pe s (P an ge on m ou n- ta in ) Fe , A s, Pb , Z n, C u, S b, A u C ar bo na te re pl ac e- m en t s ty le , p ar t o f a n IR G S Q ua rtz L 13 5- 25 6( 19 0) 5. 62 - 9. 05 N aC l C oe xi st in g im m is ci bl e ca rb on ic a nd a qu e- ou s fl ui ds a t a d ep th o f 6 -9 k m a t ~ 27 0° C . C O 2 p la ye d a si gn ifi ca nt ro le in th e or e pr ec ip ita tio n EL IO PO U LO S (2 00 0) ; E LI O PO U - LO S & K IL IA S (2 01 1) C 27 5- 34 0 (3 20 ) 2. 61 - 7. 31 H 2O -C O 2 V - - - D ec re pi ta tio n E pt ad en dr o/ R ac hi C u, F e, Z n, P b, Te , B i, C o, N i, A g Ve in st yl e, e ar ly st ag e Q ua rtz L 30 0- 38 0 (3 30 ) 1. 6 - 7 .2 N aC l C oo lin g an d di lu tio n of fl ui ds w ith a m ag - m at ic o rig in m ix ed w ith m et eo ric w at er M EL FO S et a l. (2 01 0) Ve in st yl e, la te st ag e C al ci te L 21 0- 26 0 (2 40 ) 3. 2 - 6 .3 T ha so s ( K ap sa - lin a) C u, F e, A s, Sb , B i, N i, A g, A u Ve in st yl e Q ua rtz L 12 0- 28 0 (2 00 ) 1. 0 - 3 .3 N aC l ( m i- no r K C l) C oo lin g of h yd ro th er m al fl ui ds d om in at ed by m et eo ric w at er VA V EL ID IS & M EL FO S (2 00 4) T ha so s ( Pa na gi a) C u, F e, A s, Sb , B i, N i, Te , A g, A u Ve in st yl e, p ar t o f an in tru si on re la te d sy st em Q ua rtz L 22 9- 30 5 (2 90 ) 3. 12 - 3. 28 N aC l C oe xi st in g im m is ci bl e ca rb on ic a nd a qu e- ou s fl ui ds . C O 2 p la ye d a si gn ifi ca nt ro le in th e or e pr ec ip ita tio n VA V EL ID IS e t a l. (1 99 5) C 31 3- 36 0 (3 30 ) 3. 7 - 6 .5 H 2O -C O 2 Pa go ni R ac hi C u, F e, M o, A s, Sb , P b, Z n, B i, R e, T e, A g, A u Po rp hy ry C u- M o -R e- A u st yl e (A , B v ei ns ) Q ua rtz L A : 2 47 -5 14 (4 50 ) A : 2 .4 – 1 0. 7 N aC l A - a nd B -ty pe p or ph yr y st yl e ve in s w er e de po si te d at a d ep th < 2 km fr om a b oi lin g hy dr ot he rm al fl ui d. S ub se qu en t d ilu tio n of th e m od er at el y sa lin e flu id re su lte d in th e fo rm at io n of a lo w er te m pe ra tu re a nd le ss sa lin e flu id in th e la te e pi th er m al E -v ei ns V O U D O U R IS e t a l. (2 01 3) B : 2 10 -5 13 (3 90 ) B : 1 .7 – 9 .2 L- S A : 2 04 -> 60 0 (3 70 ) A : 3 8. 9 - > 74 N aC l B : 3 61 -> 60 0 (3 90 ) B : 3 6. 4 - > 74 V A : 2 45 -5 22 (3 90 ) A : 4 .9 – 1 0. 4 N aC l B : 2 76 -4 94 (4 00 ) B : 4 .6 – 9 .2 Ep ith er m al o ve rp rin t- in g (E v ei ns ) Q ua rtz L 20 6- 37 7 (2 50 ) 1. 1 – 2. 9 N aC l C al ci te L 13 6- 33 9 (2 54 ) 1. 4 - 2 .9 N aC l M ar on ia C u, F e, M o, A s, Sb , P b, B i, R e, A u Po rp hy ry C u- M o -R e- A u st yl e Q ua rtz L 20 0- 58 0 (3 65 ) 7 – 16 N aC l B oi lin g an d co ol in g of th e flu id s M EL FO S et a l. (2 00 2) L- S 21 0- 51 5 (3 31 ) 28 – 5 5 V 27 5- 59 0 (3 93 ) 5 – 14 V- S 30 7- 45 5 (3 81 ) 38 - 51 Pe ra m a H ill Pb , Z n, B i, Se , H S- IS e pi th er m al st yl e Q ua rtz L 21 1- 33 1 (2 40 ) 0. 7 - 4 .9 N aC l C oo lin g an d di lu tio n of fl ui ds w ith a m ag - m at ic o rig in m ix ed w ith m et eo ric w at er V O U D O U R IS e t a l. (2 01 1) Te , S n, A g, A u B ar ite L 22 2- 25 6 (2 40 ) 0. 8 - 2 .3 N aC l L ou tr os H ill Fe , P b, A s, A g IS e pi th er m al st yl e B ar ite L 22 0- 32 9 (3 10 ) 8. 65 - 16 .1 4 N aC l C oo lin g of a m ag m at ic o re fl ui d po ss ib ly m ix ed w ith m et eo ric w at er Th is st ud y 165 Melfos & Voudouris: Fluid evolution in Tertiary magmatic-hydrothermal ore systems at the Rhodope metallogenic province... Thasos island also demonstrates features which are similar to an IRGS, although there is no plutonic rock occurrence on the island. In these three systems aqueous-carbonic and aqueous fluids coexist in the same samples, demonstrating that carbon diox- ide played a significant role in ore formation. It has been em- phasized that CO2 is the medium which can induce transport and precipitation of some metals, including gold, and aids the process of ore deposition by phase separation in magmatic- hydrothermal systems (LANG & BAKER, 2001; BAKER, 2002). This process is common in many intrusion related gold systems and has been cited as direct evidence for a magmatic fluid source (BAKER & LANG, 2001), although decarboniza- tion reactions in wall-rock and dissolution of CO2 from crustal hydrothermal fluids may be also encountered (LOWEN- STERN, 2001). Vein type ore mineralizations in the Rhodope metallogenic province are fault-controlled vein structures. They show simi- lar formation conditions (temperatures-salinities), but relative- ly different fluid chemical compositions. In the Eptadendro- Rachi vein system the mineralization was caused by mixing between ascending magmatic fluids and downward migrating dilute meteoric water during cooling at homogenization tem- peratures from 380° to 210°C. These fluids with low to mode- rate salinities (1.6-7.2 wt% NaCl equiv) probably show an as- sociation with a granitoid body at depth, of Tertiary age, which is not exposed at the surface. The mineralized veins in Kapsa- lina of Thasos island were formed under similar conditions by low to moderate temperatures and low density fluids (Th=120- 280°C, salinities=1.0-3.3 wt% NaCl equiv) also by cooling of the hydrothermal fluids which were derived mainly by mete- oric waters (VAVELIDIS & MELFOS, 2004). The porphyry Cu-Mo-Re-Au ore systems in Pagoni Rachi and Maronia demonstrate the contribution of the vapour phase and the highly saline fluids under boiling conditions for ore mineral precipitation. Boiling processes resulted in the genera- tion of coexisting low to moderately saline and highly saline fluids at temperatures from 300 to 500°C and the formation of the different hydrothermal alteration and ore zones (MELFOS et al., 2002; VOUDOURIS et al., 2013b). Vapour and brine dominance in both porphyry systems indicate that they origi- nate from a deep near-critical single phase fluid by decompres- sion and subsequent phase separation. In this case a low-densi- ty magmatic vapour phase is the principal agent for metal transportation, according to the model proposed by HEIN- RICH (2005), WILLIAMS-JONES & HEINRICH (2005), HURTIG & WILLIAMS-JONES (2015). The released magmatic derived vapour rich fluids can trans- port high concentrations of metals mainly Cu and Au. Rapid expansion of vapour at a shallow depth favours Cu, Fe, Mo, Re and Au saturation and precipitation of the related ore minerals. WILLIAMS-JONES & HEINRICH (2005) emphasize the role of vapour volatile and precipitation of the bulk of Cu, Au and Mo at temperatures between 320 and 425°C at a shallow depth of about 2 km. At this temperature in the presence of iron, SO2 in the fluid phase is reduced to H2S and H2SO4, which favours the precipitation of sulfides, especially pyrite, chalcopyrite, bornite and molybdenite, and increases gold solubility. The same process decreases drastically the fluid pH and causes the generation of an increasingly acidic fluid, which is associated with the related hydrothermal alteration minerals, mainly ser- icite and after mixing with groundwater, kaolinite (HEIN- RICH, 2005). The shallow epithermal veins in the Pagoni Rachi porphyry system and the epithermal overprinting identified in Maronia (MELFOS et al., 2002; VOUDOURIS et al., 2013b) were formed under lower temperature conditions due to cooling and dilution of the ore fluids, which were evolved from initial high temperature and low salinity fluids towards cooler and very low salinity hydrothermal fluids. The Perama Hill HS-IS and the Loutros IS epithermal mine- ralizations show similar temperature conditions (220-330°C) but slightly different salinities. In Perama Hill the salinities range between 0.7 to 4.9 wt% NaCl equiv indicating conden- sation of a magmatic vapour in the hydrothermal fluid with limited mixing with a fluid that had a major meteoric water component. In Loutros the salinities are higher (8.65 to 16.14 wt% NaCl equiv) which show a significant magmatic contri- bution in the hydrothermal fluid with possibly a limited mixing with meteoric water. The epithermal mineralizations in Perama Hill and Loutros were formed by cooling and dilution of the ore fluids. These fluids may represent late stage magmatic flu- ids that probably mixed with meteoric or seawater fluids at shallow levels. 5. CONCLUSIONS Fluid inclusion studies of Tertiary magmatic-hydrothermal ore mineralizations in the Rhodope metallogenic province, NE Greece, focus mainly on the chemical composition of the ore forming fluids and the formation conditions of the mineraliza- tions, the role of the magmatic and meteoric fluids and their possible interaction for the metallogenesis, the immiscibility and boiling processes and the role of CO2 in ore precipitation. The different fluid characteristics and microthermometric data indicate a variety of fluid origin conditions and sources which can affect the strategy for exploration and prospecting of ore deposits. Carbon dioxide rich hydrothermal fluids and Au and Te mineralization is common and occurs at the Kavala intrusion hosted sheeted vein system, at the Asimotrypes Au-rich car- bonate replacement mineralization and at the Panagia (Thasos) vein system. Transport and precipitation of metals like Au and Te is favoured when CO2 is present. Precipitation of the ore mineralization takes place due to the immiscibility of the car- bonic and the aqueous fluids resulting in their unmixing. These fluids have a magmatic origin with the contribution of mete- oric water and are derived from deep seated systems. Based on geological and mineralogical criteria and the fluid inclusion data, we classify all these three ore mineralizations in Kavala, Asimotrypes and Panagia of Thasos, as parts of IRGS. Cooling of magmatic hydrothermal fluids and dilution with meteoric water is a common cause of ore mineral formation in vein mineralizations (Eptadendro/Rachi, Kapsalina Thasos). Boiling and the vapour phase predominance are the most im- portant processes for metallogenesis in the porphyry deposits (Pagoni Rachi and Maronia). These features are common in 166 Vol. 69/1 porphyry systems and the mining industry should focus on them during early stage exploration work. The epithermal veins overprinting the Pagoni Rachi and the Maronia porphyry systems and the HS-IS epithermal system in Perama Hill and the IS epithermal mineralization in Loutros are characterized by low to moderate temperatures and low to moderate salinities. Cooling and dilution of the ore fluids are the main processes for gold precipitation. Following these criteria of the fluid characteristics, explora- tion for gold, rare and critical metals could have encouraging results to promote the mining industry in the future by using fluid inclusion techniques. Acknowledgment The authors would like to thank the reviewer Ladislav PALINKAŠ and an Anonymous Reviewer for their helpful and constructive comments on the manuscript. Ms Alisa MARTEK is especially thanked for her editorial handling of the paper. REFERENCES BAKER, T. & LANG, J.R. (2001): Fluid inclusion characteristics of intrusion- related gold mineralization, Tombstone–Tungsten magmatic belt, Yukon Territory, Canada.– Mineralium Deposita, 36, 563–582. BAKER, T. (2002): Emplacement depth and carbon dioxide-rich fluid inclu- sions in intrusion-related gold deposits.– Economic Geology, 97, 1111–1117. doi: 10.2113/97.5.1111 BLUNDELL, D., ARNDT, N., COBBLOD P.R. & HEINRICH, C. (2005): Processes of tectonism, magmatism and mineralization: Lessons from Eu- rope.– Ore Geology Reviews, 27, 333–349. doi: 10.1016/j.oregeor- ev.2005.07.003 BODNAR, R.J. (2003): Introduction to fluid inclusions.– In: SAMSON, I.M., ANDERSON, A.J. & MARSHALL, D.D. (eds.): Fluid inclusions: analy- sis and interpretation. Mineralogical Association of Canada Short Course, 32, 1–8. BONSALL, T.A., SPRY, P.G., VOUDOURIS, P.C., TOMBROS, S., SEY- MOUR, K.S. & MELFOS, V. (2011): The geochemistry of carbonate-re- placement Pb-Zn-Ag mineralization in the Lavrion district, Attica, Greece: Fluid inclusion, stable isotope, and rare earth element studies.– Economic Geology, 106, 619–651. doi: 10.2113/eco ngeo.106.4.619 BORTNIKOV, N.S. (2006): Geochemistry and origin of the ore-forming fluids in hydrothermal-magmatic systems in tectonically active zones.– Geolo- gy of Ore Deposits, 48, 1–22. doi: 10.1134/S1075701506010016 CHANG, Z., HEDENQUIST, J.W., WHITE, N.C., COOKE, D.R., ROACH, M., DEYELL, C.L. & CUISON, A.L. (2011): Exploration tools for linked porphyry and epithermal deposits: Example from the Mankayan intru- sion-centered Cu-Au district, Luzon, Philippines.– Economic Geology, 106, 1365–1398. doi: 10.2113/eco ngeo.106.8.1365 CHATZIKIRKOU, A. & MICHAILIDIS K. (2004): The Eptadendro and Ra- chi sulfide mineralization, Eastern Rhodope, Greece.– Bull. Geol. Soc. Greece, 36, 397–405 (in Greek with English abstract). CHRISTOFIDES, G. (1996): Tertiary magmatism in the Greek Rhodope Mas- sif, northern Greece: Granitic plutons.– In: KNEZEVIC, V. & KRSTIC, B. (eds.): Terranes of Serbia: The Formation of the Geologic Framework of Serbia and the Adjacent Regions. University of Belgrade, Belgrade, 155–160. CHRISTOFIDES, G., PECSKAY, Z., ELEFTHERIADIS, G., SOLDATOS, T. KORONEOS, A. (2004): The Tertiary Evros volcanic rocks (Thrace, northeastern Greece): petrology and K/Ar geochronology.– Geologica Carpathica, 55, 397–410. DE BOORDER, H., SPAKMAN, W., WHITE, S.H. & WORTEL, M.J.R. (1998): Late Cenozoic mineralization, orogenic collapse and slab detach- ment in the European Alpine Belt.– Earth and Planetary Science Letters, 164, 569–575. doi: 10.1016/S0012-821X(98)00247-7 ELEFTHERIADIS, G. & KORONEOS, A. (2003): Geochemistry and petro- genesis of post-collision Pangeon granitoids in central Macedonia, North- ern Greece.– Chemie der Erde-Geochemistry, 63, 364–389. doi: 10.1078/0009-2819-00024 ELIOPOULOS, D.G. & KILIAS, S.P. (2011): Marble-hosted submicroscopic gold mineralization at Asimotrypes area, mount Pangeon, southern Rho- dope core complex, Greece.– Economic Geology, 106, 751–780. doi: 10.2113/econgeo.106.5.751 ELIOPOULOS, D.G. (2000): Geochemistry and origin of the Asimotrypes carbonatehosted mesothermal gold deposit, Pangeon mountain, N. Greece.– Unpublished PhD thesis, University of Southampton, UK, 244 p. FORNADEL, A.P., SPRY, P.G., MELFOS, V., VAVELIDIS, M. & VOU- DOURIS, P. (2011): Is the Palea Kavala Bi-Te-Pb-Sb ± Au district, north- eastern Greece, an intrusion-related system?– Ore Geology Reviews, 39, 119–133. doi: 10.1016/j.oregeorev.2011.01.002 GEORGAKOPOULOU, M., BASSIAKOS, Y. & PHILANIOTOU, O. (2011): Seriphos surfaces: a study of copper slag heaps and copper sources in the context of Early Bronze Age Aegean metal production.– Archaeometry, 53, 123–145. doi: 10.1111/j.1475-4754.2010.00529.x GOLDSTEIN, R.H. & REYNOLDS, T.J. (1994): Systematics of fluid inclu- sions in diageneticminerals.– SEPM Short Course 31, Tulsa, 199 p. HART, C.J.R. (2005): Classifying, distinguishing and exploring for intrusion- related gold systems. The Gangue, Geological Association of Canada, Mineral Deposits Division Newsletter, 87, 1 and 4–9. HEINRICH, C.A. (2005): The physical and chemical evolution of low-salinity magmatic fluids at the porphyry to epithermal transition: a thermodyna- mic study.– Mineralium Deposita, 39, 864–889. doi: 10.1007/s00126- 004-0461-9 HEINRICH, C.A. (2007): Fluid-fluid interactions in magmatic-hydrothermal ore formation.– Reviews in Mineralogy and Geochemistry, 65, 363–387. doi: 10.2138/rmg.2007.65.11 HURTIG, N.C. & WILLIAMS-JONES, A.E. (2015): Porphyry-epithermal Au-Ag-Mo ore formation by vapour-like fluids: New insights from geo- chemical modeling.– Geology, G36685-1. JOLIVET, L. & BRUN, J.P. (2010): Cenozoic geodynamic evolution of the Aegean.– International Journal of Earth Sciences, 99, 109–138. doi: 10.1007/s00531-008-0366-4 KESLER, S.E., BODNAR, R.J. & MERNAGH, T.P. (2013): Role of fluid and melt inclusion studies in geologic research.– Geofluids, 13, 398–404. doi: 10.1111/gfl.12055 LANG, J.R. & BAKER, T. (2001): Intrusion-related gold systems: the present level of understanding.– Mineralium Deposita, 36, 477–489. doi: 10.1007/ s001260100184 LEHMANN, B., DIETRICH, A. & WALLIANOS, A. (2000): From rocks to ore.– International Journal of Earth Sciences, 89, 284–294. doi: 10.1007/ s005310000085 LOWENSTERN, J.B. (2001): Carbon dioxide in magmas and implications for hydrothermal systems.– Mineralium Deposita, 36, 490–502. doi: 10.1007/s001260100185 MARCHEV, P., KAISER-ROHRMEIER, M., HEINRICH, C., OVTCHA- ROVA, M., VON QUADT, A. & RAICHEVA, R. (2005): Hydrothermal ore deposits related to post-orogenic extensional magmatism and core complex formation: the Rhodope Massif of Bulgaria and Greece.– Ore Geology Reviews, 27, 53–89. doi: 10.1016/j.oregeorev.2005.07.027 MELFOS, V. & VOUDOURIS, P.C. (2012): Geological, mineralogical and geochemical aspects for critical and rare metals in Greece.– Minerals, 2, 300–317. doi: 10.3390/min2040300 MELFOS, V., CHATZIKIRKOU, A., MICHAILIDIS, K. & VOUDOURIS, P. (2010): Fluids related to remobilization of Mesozoic sulfide mineraliza- tion in the Eptadendro-Rachi region in eastern Rhodope, Thrace, Greece.– In CHRISTOFIDES et al. (eds.): Proceedings of the XIX Congress of the Carpathian-Balkan Geological Association, Thessaloniki. Scientific An- nals of the School of Geology A.U.Th., 100, 343–350. MELFOS, V., VAVELIDIS, M., CHRISTOFIDES, G. & SEIDEL, E. (2002): Origin and evolution of the Tertiary Maronia porphyry copper-molybde- num deposit, Thrace, Greece.– Mineralium Deposita, 37, 648–668. doi: 10.1007/s00126-002-0277-4 MELFOS, V., VOUDOURIS, P., VAVELIDIS, M. & SPRY, P. (2008): Micro- thermometric results and formation conditions of a new intrusion-related Bi-Te-Pb-Sb±Au deposit in the Kavala pluton, Greece.– In: Proceedings of the XIII All-Russian conference on thermobarogeochemistry in con- junction with IV APIFIS symposium, vol 2, 165–168. NERANTZIS, N. (2009): Pillars of Power: Silver and Steel of the Ottoman Empire.– Mediterranean Archaeology and Archaeometry, 9, 71–85. PE-PIPER, G. & PIPER D.J. (2001): Late Cenozoic, post-collisional Aegean igneous rocks: Nd, Pb and Sr isotopic constraints on petrogenetic and tectonic models.– Geological Magazine, 138, 653–668. PE-PIPER, G. & PIPER, D.J. (2002): The Igneous Rocks of Greece.– Ge- broder Borntraeger, Berlin, Germany, 573 p. doi: 10.1017/ S0016756801005957 PERNICKA, E., GENTNER, W., WAGNER, G.A., VAVELIDIS, M. & GALE, N.H. (1981): Ancient lead and silver production on Thasos (Greece).– Revue d’Archéométrie, 1, 227–237. doi: 10.3406/ar- sci.1981.1151 ROEDDER, Ε. (1977): Fluid inclusions as tools in mineral exploration.– Eco- nomic Geology, 72, 503–525. 167 Melfos & Voudouris: Fluid evolution in Tertiary magmatic-hydrothermal ore systems at the Rhodope metallogenic province... SIMMONS, S.F. & BROWN, K.L. (2006): Gold in magmatic hydrothermal solutions and the rapid formation of a giant ore deposit.– Science, 314, 288–290. doi: 10.1126/science.1132866 TSIRAMBIDES, A. & FILIPPIDIS, A. (2012): Exploration key to growing Greek industry.– Industrial Minerals, 533, 44–47. VAVELIDIS, M. & AMSTUTZ, G.C. (1983): New genetic investigations on the Pb-Zn deposits of Thasos (Greece).– In: SCHNEIDER, H.J. (ed.): Mineral Deposits of the Alpine Epoch in Europe, Berlin: Springer-Verlag, 359–365. doi: 10.1007/978-3-642-68988-8_36 VAVELIDIS, M. & MELFOS, V. (2004): Bi-Ag-bearing tetrahedrite-tennan- tite in the Kapsalina copper mineralisation, Thasos island, Northern Greece.– Neues Jahrbuch für Mineralogie – Abhandlungen (Journal of Mineralogy and Geochemistry), 180, 149–169. doi: 10.1127/0077- 7757/2004/0180-0149 VAVELIDIS, M., GIALOGLOU, G., MELFOS, V. & WAGNER, G.A. (1996): Goldgrube in Palaea Kavala-Griechenland: Entdeckung von Skaptehyle?– Erzmetall, 49, 547–554. VAVELIDIS, M., GIALOGLOU, G., WAGNER, G.A. & PERNICKA, E. (1988): Die Bundmetall- und Eisen-Mangan-Lagerstätten von Thasos.– Anschnitt, 6, 40–58. VAVELIDIS, M., SCHMIDT-MUMM, A. & MELFOS, V. (1995): Microther- mometric investigations in the Bi-Te-Ag bearing Cu-mineralization of Panagia area, Thasos Island, Greece.– Boletin Soc. Espan. Mineralogia, 18–1, 261–262. VAXEVANOPOULOS, M., VAVELIDIS, M., MALAMIDOU, D., MELFOS, V. & PAVLIDIS, S. (in press): Excavations in ancient mines and metal- lurgical sites on the Pangaion mountain (Asimotrypes, Valtouda).– The Archaeological Work in Macedonia and Thrace (in press, in Greek). VOUDOURIS, P. & ALFIERIS, D. (2005): New porphyry Cu±Mo occur- rences in the north-eastern Aegean, Greece: Ore mineralogy and epither- mal relationships.– In: MAO, J. & BIERLEIN, F.P. (eds.): Mineral De- posit Research: Meeting the Global Challenge, Springer Berlin Heidel- berg, 473–476. VOUDOURIS, P. (2006): Comparative mineralogical study of Tertiary Te-rich epithermal and porphyry systems in northeastern Greece.– Mineralogy and Petrology, 87, 241–275. VOUDOURIS, P., MELFOS, V., SPRY, P.G., BINDI, L., MORITZ, R., OR- TELLI, M. & KARTAL, T. (2013a): Extremely Re-rich molybdenite from porphyry Cu-Mo-Au prospects in northeastern Greece: Mode of occur- rence causes of enrichment, and implications for gold exploration.– Mine- rals, 3, 165–191. doi: 10.3390/min3020165 VOUDOURIS, P., MELFOS, V., SPRY, P.G., KARTAL, T., SCHLEICHER, H., MORITZ, R. & ORTELLI, M. (2013b): The Pagoni Rachi/Kirki Cu- Mo±Re±Au deposit, Northern Greece: Mineralogical and fluid inclusion constrains on the evolution of a telescoped porphyry-epithermal system.– Canadian Mineralogist, 51, 411–442. VOUDOURIS, P., MELFOS, V., SPRY, P.G., MORITZ, R., PAPAVASSIL- IOU, C. & FALALAKIS, G. (2011): Mineralogy and geochemical envi- ronment of formation of the Perama Hill high-sulfidation epithermal Au- Ag-Te-Se deposit, Petrota Graben, NE Greece.– Mineralogy and Petrolo- gy, 103, 79–100. doi: 10.1007/s00710-011-0160-z VOUDOURIS, P.C., MELFOS, V., SPRY, P.G., BINDI, L., KARTAL, T., ARIKAS, K., MORITZ, R. & ORTELLI, M. (2009): Rhenium-rich mo- lybdenite and rheniite (ReS2) in the Pagoni Rachi-Kirki Mo-Cu-Te-Ag- Au deposit, Northern Greece: Implications for the rhenium geochemistry of porphyry style Cu-Mo and Mo mineralization.– Canadian Mineralo- gist, 47, 1013–1036. WILKINSON, J.J. (2001): Fluid inclusions in hydrothermal ore deposits.– Lithos, 55, 229–272. doi: 10.1016/S0024-4937(00)00047-5 WILLIAMS-JONES, A.E. & HEINRICH, C.A. (2005): 100th Anniversary special paper: vapour transport of metals and the formation of magmatic- hydrothermal ore deposits.– Economic Geology, 100, 1287–1312. doi: 10.2113/gsecongeo.100.7.1287