Arsenian Pyrite and Cinnabar from Active Submarine Nearshore Vents, Paleochori Bay, Milos Island, Greece

Active, shallow-water (2–10 m below sea level) and low temperature (up to 115 °C) hydrothermal venting at Paleochori Bay, nearshore Milos Island, Greece, discharges CO2 and H2S rich vapors (e.g., low-Cl fluid) and high-salinity liquids, which leads to a diverse assemblage of sulfide and alteration phases in an area of approximately 1 km2. Volcaniclastic detritus recovered from the seafloor is cemented by hydrothermal pyrite and marcasite, while semi-massive to massive pyrite-marcasite constitute mounds and chimney-like edifices. Paragenetic relationships indicate deposition of two distinct mineralogical assemblages related to the venting of high-Cl and low-Cl fluids, respectively: (1) colloform As- and Hg-bearing pyrite (Py I), associated with marcasite, calcite, and apatite, as well as (2) porous and/or massive As-rich pyrite (Py II), associated with barite, alunite/jarosite, and late-stage hydrous ferric oxides. Mercury, in the form of cinnabar, occurs within the As-rich pyrite (Py I) layers, usually forming distinct cinnabar-enriched micro-layers. Arsenic in colloform pyrite I shows a negative correlation with S indicating that As1− dominates in the pyrite structure suggesting formation from a relatively reducing As-rich fluid at conditions similar to low-sulfidation epithermal systems. On the contrary, As3+ dominates in the structure of porous to massive pyrite II suggesting deposition from a sulfate-dominated fluid with lower pH and higher fO2. Bulk sulfide data of pyrite-bearing hydrothermal precipitates also show elevated As (up to 2587 ppm) together with various epithermal-type elements, such as Sb (up to 274 ppm), Tl (up to 513 ppm), and Hg (up to 34 ppm) suggesting an epithermal nature for the hydrothermal activity at Paleochori Bay. Textural relationships indicate a contemporaneous deposition of As and Hg, which is suggested to be the result of venting from both high-salinity, liquid-dominated, as well as CO2- and H2S-rich vapor-dominated fluids that formed during fluid boiling. The CO2- and H2S-rich vapor that physically separated during fluid boiling from the high-salinity liquid led to calcite formation upon condensation in seawater together with the precipitation of As- and Hg-bearing pyrite I. This also led to the formation of sulfuric acid, thereby causing leaching and dissolution of primary iron-rich minerals in the volcaniclastic sediments, finally resulting in pyrite II precipitation in association with alunite/jarosite. The Paleochori vents contain the first documented occurrence of cinnabar on the seafloor in the Mediterranean area and provide an important link between offshore hydrothermal activity and the onshore mercury and arsenic mineralizing system on Milos Island. The results of this study therefore demonstrate that metal and metalloid precipitation in shallow-water continental arc environments is controlled by epithermal processes known from their subaerial analogues.


Figure 1.
Map showing the location of hydrothermal deposits discovered on the seafloor including locations of some very shallow and shallow-water hydrothermal vents discussed in this study (adapted from Jamieson et al. [42]).
Compared to their deeper arc-related settings (200 to 1500 meters below sea level, i.e. mbsl), the shallow water depth (<200 mbsl) of these systems implies that they typically occur nearshore to a volcanic center characterized by a thick volcaniclastic cover, which has a significant imprint on the fluid composition and the conditions of formation of the hydrothermal precipitates [30,32,41,[43][44][45][46][47][48][49][50]. These shallow aqueous hydrothermal systems are characterized by fluid boiling caused by the increasing density contrast of vapor and liquid with decreasing temperature and pressure conditions [51]. This leads to the discharge of both gas bubbles (vapor-dominated fluids) and liquids at "gasohydrothermal" shallow marine vent sites [28,40,52,53].
Milos Island, Greece is known for its mineral wealth, exploited from ancient to recent times. Industrial minerals, as well as base and precious metal epithermal deposits, are related to long-standing hydrothermal activity since the Pliocene [54][55][56][57][58][59][60][61]. The active nearshore expressions (e.g., Paleochori Bay) of the Milos epithermal/geothermal system have attracted research, mainly with respect to metal and metalloid enrichment and fractionation in the hydrothermal fluids (e.g., [28,[44][45][46]). However, only little attention has been Map showing the location of hydrothermal deposits discovered on the seafloor including locations of some very shallow and shallow-water hydrothermal vents discussed in this study (adapted from Jamieson et al. [42]).
Compared to their deeper arc-related settings (200 to 1500 m below sea level, i.e. mbsl), the shallow water depth (<200 mbsl) of these systems implies that they typically occur nearshore to a volcanic center characterized by a thick volcaniclastic cover, which has a significant imprint on the fluid composition and the conditions of formation of the hydrothermal precipitates [30,32,41,[43][44][45][46][47][48][49][50]. These shallow aqueous hydrothermal systems are characterized by fluid boiling caused by the increasing density contrast of vapor and liquid with decreasing temperature and pressure conditions [51]. This leads to the discharge of both gas bubbles (vapor-dominated fluids) and liquids at "gasohydrothermal" shallow marine vent sites [28,40,52,53].
Milos Island, Greece is known for its mineral wealth, exploited from ancient to recent times. Industrial minerals, as well as base and precious metal epithermal deposits, are related to long-standing hydrothermal activity since the Pliocene [54][55][56][57][58][59][60][61]. The active nearshore expressions (e.g., Paleochori Bay) of the Milos epithermal/geothermal system have attracted research, mainly with respect to metal and metalloid enrichment and fractionation in the hydrothermal fluids (e.g., [28,[44][45][46]). However, only little attention has been given to the subaqueous sulfide precipitation, despite its profound importance for providing a better understanding of the mineralization processes in shallow submarine and onshore environments [40,41,43,[47][48][49][50]53,[62][63][64][65][66]. The link between the nearshore and onshore mineralization at Milos is still poorly investigated [57,60], and thus improving our understanding of the seafloor hydrothermal processes at Paleochori Bay will also help to better constrain the formation conditions of the onshore epithermal mineralization of Milos Island.
This study expands on previous work of Kati et al. [67,68], and presents new data from Paleochori Bay reporting on the occurrence of Hg and Tl enrichments associated with As-rich pyrite mineralization in actively forming hydrothermal vents in a very shallow marine continental volcanic arc setting. Active cinnabar formation is a rare feature of submarine hydrothermal systems, and this study reports for the first time on the occurrence of actively formed cinnabar on the sea floor in the Aegean and the Mediterranean area. This study aims to present an integrated view of the atypical mineralization of Paleochori Bay focusing on the mineral constituents, bulk and mineral chemistry, mineralization sequence, and physicochemical formation conditions. This provides new insights into the hydrothermal metallogeny of Milos Island and demonstrates that metal and metalloid precipitation in shallow-water continental arc environments is controlled by epithermal processes known from their subaerial analogues.

Materials and Methods
Sulfide grab samples were recovered from ten sites in an area of about 450 m × 250 m on the seafloor of Paleochori Bay (Table 1 and Figures 2 and 3). Sampling took place over four field seasons (1996)(1997)(1998)(1999) by scuba diving. The collected material was taken from sulfide edifices (mounds and chimney-like structures) and metalliferous sediment in the immediate vicinity of areas of active discharge, where hydrothermal fluids were sampled and described in detail by Valsami-Jones et al. [41]. The samples are considered to be representative, adequately describing the situation at Paleochori Bay.
A total of 31 mineralized samples were examined by optical microscopy, X-ray diffraction (XRD), scanning electron microscopy (SEM), electron probe micro-analysis (EPMA), and cathodoluminescence. X-ray diffraction measurements (n = 26) were carried out in the Department of Geology and Geoenvironment at the National and Kapodistrian University of Athens, using a Siemens (Bruker) D5005 X-ray diffractometer at the following conditions: CuKa radiation (40 kV, 40 mA), 3-65 • 2θ range, 0.02 • step size, and 1 s step time. For the mineralogical identification, the DIFFRAC v. 2.41 and EVA v. 10.0 software package was used.
Cathodoluminescence images were obtained at the Natural History Museum, London, using a CITL CL8200 Mk4 cold cathode system (350 µA and 17 kV gun current) mounted on a Leica DM LM microscope.
Bulk samples were digested in 1:1:1 aqua regia and then analyzed for selected trace elements by ICP-MS at ACME Analytical Laboratories. Detection limits were: 0.2 ppb for Au; 2 ppb for Ag; 5 ppb for Hg; 0.01 ppm for Mo, Cu, Pb, and Cd; 0.02 ppm for Bi, Sb, Tl, and Te; 0.1 ppm for Ni, Ga, Ge, Sn, Co, Zn, As, and Se; 1 ppm for Mn; and 0.01 wt % for Fe.

Regional and Local Geology
Milos Island is a dormant volcanic center located in the central part of the Lower Pliocene to recent South Aegean active volcanic arc, which also includes the volcanic centers of Methana and Poros in its western part, and Santorini, Kos, and Nisyros in its eastern part (Figure 2a). The Aegean volcanic arc is generally considered a continental arc formed by the subduction of the African plate beneath the continental Aegean microplate [70][71][72]. Tertiary to Quaternary magmatism in the Aegean region occurred mostly in a post-collisional setting behind the active Hellenic subduction zone [71,73].
Volcanic activity in western Milos spans a period from~3.5 Ma to present and originated from several emergent eruptive centers characterized by both explosive and effusive activities [74][75][76][77] (Figure 2b). Volcanic activity led to the deposition of shallow marine pumiceous tuffs and a calc-alkaline sequence of basalts, andesites, dacites, and rhyolites covered by alluvial sediments that both overlie an Upper Miocene to Lower Pliocene sedimentary sequence and a Mesozoic metamorphic basement of the Cycladic blueschist unit (Figure 2b) [74,78,79].
The volcanic units can be subdivided into four main phases and are ( Figure 2b): (1) Lower to Upper Pliocene felsic submarine units including pumice-rich pyroclastics, intrusive rhyolitic subvolcanic bodies, and lavas in the western part of the island; (2) Upper Pliocene to Middle Pleistocene submarine pyroclastics and andesitic-dacitic flow domes and lavas, and submarine to subaerial dacitic-rhyodacitic flow banded domes associated with lithic-and pumice-rich pyroclastics in the western part of the island; (3) Lower to Middle Pleistocene submarine felsic domes-cryprodomes and pyroclastics, and subaerial rhyolitic lavas mainly in the eastern and northeastern parts of the island; and (4) Late Pleistocene subaerial explosive eruptions followed by rhyolitic lavas and accompanied by widespread phreatic activity, concentrated along the eastern coast of Milos Gulf on a well-defined NW-trending lineament.
Miocene to Pliocene extensional tectonics resulted in four main fault trends (NW-SE, N-S, NE-SW, and E-W) creating a series of horsts and graben structures controlling the volcanic and hydrothermal activity of the island until today (Figure 2b). Paleochori Bay is located nearshore the southeast coast of Milos ( Figure 2b). The socalled "green lahar" exposed onshore in the Paleochori area consists of volcanic ash with fragments of metamorphic, sedimentary, and volcanic rocks (Figure 2c). The seafloor at Paleochori Bay is covered by sandy volcaniclastic material with larger rock fragments towards the eastern and western coast. Paleochori Bay is covered by sandy volcaniclastic material with larger rock fragments towards the eastern and western coast.

Hydrothermal Venting
Venting activity around the eastern part of Milos island, including Paleochori Bay, occurs over an area of approximately 35 km 2 from onshore to about 100 m water depth, and is characterized by gas-and liquid-dominated fluids that are emitted through the unconsolidated volcaniclastic cover [40,41,[45][46][47]50,53,63,81] (Figure 3). The hydrothermal fluids are highly sulfidic (up to 3 mM H 2 S) and reach temperatures of up to 122 • C [41,82]. The active submarine hydrothermal system of Paleochori Bay discharges vapor-and liquiddominated aqueous fluids with average temperatures (91 ± 23 • C). Vapors released from both subaerial fumaroles and submarine vents at the Paleochori Bay consist mainly of CO 2 (55-92 vol. %), H 2 (<3 vol. %), CH 4 (<10 vol. %), and H 2 S (<8 vol. %) [40]. Valsami-Jones et al. [41] and Wu et al. [28,44,45] identified two types of aqueous fluids at Paleochori Bay: (1) low-Cl fluids containing low concentrations of alkalis (K, Li, Na) and Ca, and high concentrations of silica and sulfate, as well as (2) high-Cl fluids containing high concentrations of alkalis and Ca, and lower concentrations of silica and sulfate compared to seawater. High contents of Zn, Cd, Mn, and Pb in the high-Cl fluids suggest that these metals are likely transported as chloride-complexes [41]. Both types of fluids occur next to each other in vents that are only a few meters apart discharging fluids with As concentrations of up 39 µM in the high-Cl fluids and up to 78 µM in the low-Cl fluids [46]. These two fluid types at Paleochori Bay are interpreted to be vapor-and liquid-rich fractions that formed by boiling of deep reservoir (~1-2 km) and shallow seawater-like fluids [28,41,45]. The low-Cl fluids are likely the result of vapor condensation into seawater in a steam-heated environment, whereas the high-Cl fluids are diluted in Cl as a result of seawater mixing [46].
trations of alkalis and Ca, and lower concentrations of silica and sulfate compared to seawater. High contents of Zn, Cd, Mn, and Pb in the high-Cl fluids suggest that these metals are likely transported as chloride-complexes [41]. Both types of fluids occur next to each other in vents that are only a few meters apart discharging fluids with As concentrations of up 39 μM in the high-Cl fluids and up to 78 μM in the low-Cl fluids [46]. These two fluid types at Paleochori Bay are interpreted to be vapor-and liquid-rich fractions that formed by boiling of deep reservoir (~1-2 km) and shallow seawater-like fluids [28,41,45]. The low-Cl fluids are likely the result of vapor condensation into seawater in a steamheated environment, whereas the high-Cl fluids are diluted in Cl as a result of seawater mixing [46]. Hydrothermal activity with associated sulfide mineralization is concentrated in the center of the bay (Figures 2c and 3). Intense streams of gas bubbles (Figure 3a-c) and high salinity liquid discharge ( Figure 3d) emerge through the sandy seafloor or between sulfide mounds and conical edifices. According to Valsami-Jones et al. [41], venting associated with the formation of sulfide mounds from the hydrothermal fallout is also related to intermittent gas bubbling through the mound structure. Occasionally, the intensity of venting associated with gas bubbles is such that it raises grey to black sand plumes (fluid plumes with entrained sand) ( Figure 3b).
Yellow, orange, white (Figure 3c,f), and brown patches on the sandy seafloor of several meters in diameter are associated with diffuse venting and the deposition of S-and As-rich substances and/or microbial mats [83][84][85][86]. The different colors correlate well to different seafloor temperatures and precipitate compositions, with the hottest areas (>90°C ) covered mostly by yellow native sulfur and orange-colored sulfur/arsenic compounds and biogenic material [85]. White microbial mats occur throughout the bay and are associated with amorphous silica and native sulfur forming at lower temperatures (~45-85 °C) Hydrothermal activity with associated sulfide mineralization is concentrated in the center of the bay (Figures 2c and 3). Intense streams of gas bubbles (Figure 3a-c) and high salinity liquid discharge ( Figure 3d) emerge through the sandy seafloor or between sulfide mounds and conical edifices. According to Valsami-Jones et al. [41], venting associated with the formation of sulfide mounds from the hydrothermal fallout is also related to intermittent gas bubbling through the mound structure. Occasionally, the intensity of venting associated with gas bubbles is such that it raises grey to black sand plumes (fluid plumes with entrained sand) ( Figure 3b).
Yellow, orange, white ( Figure 3c,f), and brown patches on the sandy seafloor of several meters in diameter are associated with diffuse venting and the deposition of S-and As-rich substances and/or microbial mats [83][84][85][86]. The different colors correlate well to different seafloor temperatures and precipitate compositions, with the hottest areas (>90 • C) covered mostly by yellow native sulfur and orange-colored sulfur/arsenic compounds and biogenic material [85]. White microbial mats occur throughout the bay and are associated with amorphous silica and native sulfur forming at lower temperatures (~45-85 • C) [53]. Brown colored manganese and/or hydrous ferric oxide layers are related to the lowest hydrothermal temperature regime (~30-35 • C) [46,85]. Microbial mats host chemolithotrophic sulfur oxidizing and sulfate-reducing microbes that facilitate the precipitation of various S-bearing mineral phases, including elemental S and Fe-sulfides (e.g., [50,86,87]. The shallow marine hydrothermal fluids discharging at Paleochori Bay are related to the active Milos geothermal field [67,88], possibly representing a modern analogue to the onshore Profitis Ilias-Chondro Vouno epithermal Au-Ag mineralization [57]. Stable (δD and δ 18 O) and radiogenic isotope data ( 87 Sr/ 86 Sr), as well as Br/Cl and I/Cl ratios, suggest a seawater and meteoric water signature modified by water-rock interaction and fluid boiling in the geothermal fluids [45,46,57,89].

Hydrothermal Precipitates
The investigated samples represent different types of mineralization, including sulfide edifices with mound-like morphology and metalliferous sediments with sulfideencrustations in the immediate vicinity of areas with active fluid discharge (Figures 3 and 4, Table 1). Precipitate mineralogy and hydrothermal fluid composition (all fluid data from Valsami-Jones et al. [41]) vary between the sampling sites (n = 10) and with the mineralization-type (Table 1). Fluid temperature and pH range from 49 • C to 110 • C and 3.5 to 7.6, respectively, at a water depth of 2.0 m to 4.5 m (Table 1). Table 1. Characteristics of hydrothermal precipitates (this study) and fluids at Paleochori Bay. All fluid-and water-depth data are from Valsami-Jones et al. [41] and Wu et al. [44,45]. The precipitates are associated with different styles of fluid discharge including venting associated with high-salinity liquids and/or gas bubbles and leading to the formation of small mound-like structures ( Figure 3). The sulfide mounds form solid irregular edifices up to 0.5 m high and 1 m wide, locally including chimney-like structures (e.g., conical edifices) (Figure 3a-e). The mounds are breccia structures consisting of detrital siliceous grains (e.g., sand grains of volcaniclastic sands), cemented by hydrothermal minerals, mostly a pyrite matrix (Figure 4a,b). Locally, a good development of colloform bands of pyrite and/or marcasite occurs ( Figure 4c). Small cracks act as pathways for vigorous bubbling and fluid discharge, some of which are sealed by sulfide deposition.

Site
The metalliferous sediments occur on the seafloor in areas of active fluid discharge through volcaniclastic material (Figure 3a,b,f), where they locally also cover the mounds. The metalliferous sediments form individual hydrothermal blackish precipitates on the sandy seafloor, and/or discontinuous sulfide-bearing encrustations on the surfaces of loose sand fragments (Figure 4d,e). Similarly, larger detrital rock fragments (up to several cm in size) host sulfides (mostly pyrite and marcasite) as coatings and encrustations (~1-3 mm) ( Figure 4f). Based on in situ observations, all or parts of these structures are annual features, which grow during the summer, and collapse in the winter as a result of rough weather [41].
The precipitates are associated with different styles of fluid discharge including venting associated with high-salinity liquids and/or gas bubbles and leading to the formation of small mound-like structures ( Figure 3). The sulfide mounds form solid irregular edifices up to 0.5 m high and 1 m wide, locally including chimney-like structures (e.g., conical edifices) (Figure 3a-e). The mounds are breccia structures consisting of detrital siliceous grains (e.g., sand grains of volcaniclastic sands), cemented by hydrothermal minerals, mostly a pyrite matrix (Figure 4a,b). Locally, a good development of colloform bands of pyrite and/or marcasite occurs ( Figure 4c). Small cracks act as pathways for vigorous bubbling and fluid discharge, some of which are sealed by sulfide deposition.
The metalliferous sediments occur on the seafloor in areas of active fluid discharge through volcaniclastic material (Figure 3a,b,f), where they locally also cover the mounds. The metalliferous sediments form individual hydrothermal blackish precipitates on the sandy seafloor, and/or discontinuous sulfide-bearing encrustations on the surfaces of loose sand fragments (Figure 4d,e). Similarly, larger detrital rock fragments (up to several cm in size) host sulfides (mostly pyrite and marcasite) as coatings and encrustations (~1-3 mm) (Figure 4f). Based on in situ observations, all or parts of these structures are annual features, which grow during the summer, and collapse in the winter as a result of rough weather [41].

Bulk Ore Geochemistry
Bulk sulfide analyses of the different mineralization types at Paleochori Bay (Table  2) generally indicate high concentrations in elements like As (up to 2587 ppm), Mn (up to 2360 ppm), Co (up to 813 ppm), Tl (up to 513 ppm), Sb (up to 274 ppm), and Hg (up to 34 ppm). Chimney-like structures tend to be enriched in most trace elements including Fe,

Bulk Ore Geochemistry
Bulk sulfide analyses of the different mineralization types at Paleochori Bay (Table 2) generally indicate high concentrations in elements like As (up to 2587 ppm), Mn (up to 2360 ppm), Co (up to 813 ppm), Tl (up to 513 ppm), Sb (up to 274 ppm), and Hg (up to 34 ppm). Chimney-like structures tend to be enriched in most trace elements including Fe, Mn, As, Hg, Cu, Pb, and Zn compared to the metalliferous sediments (Table 2 and Figure 5). This is likely due to higher sulfide abundances in the chimneys or a dilution by the more abundant detrital material in the sediments. Table 2. Bulk ore analyses (ICP-MS) of sulfide samples from Paleochori Bay. Iron in wt %, Au and Ag in ppb, and all other elements in ppm. Sample numbers, as in Figure 4 and Table 1. Mn, As, Hg, Cu, Pb, and Zn compared to the metalliferous sediments (Table 2 and Figure  5). This is likely due to higher sulfide abundances in the chimneys or a dilution by the more abundant detrital material in the sediments.

Mineralogy
The mound-like structures and the metalliferous sediments are heterogeneous on a small scale (~1 cm). However, the sulfides of both mineralization styles display similar textural and temporal relationships. The main framework material consists of volcaniclastic sandstones including fragments of hydrothermal altered and mineralized volcanic and metamorphic country rock, which are mostly cemented by sulfides like pyrite and marcasite (Figure 6a-c) or by other hydrothermal precipitates like barite and calcite (Figure 6a,d).

Mineralogy
The mound-like structures and the metalliferous sediments are heterogeneous on a small scale (~1 cm). However, the sulfides of both mineralization styles display similar textural and temporal relationships. The main framework material consists of volcaniclastic sandstones including fragments of hydrothermal altered and mineralized volcanic and metamorphic country rock, which are mostly cemented by sulfides like pyrite and marcasite (Figure 6a  The detrital material mostly consists of quartz, feldspars (commonly intensely altered), muscovite, biotite, chlorite, rutile, and titanite, as well as subordinate clasts of mica schists, volcaniclastics, and microcrystalline silica reflecting the country rock of Paleochori Bay. The iron sulfide cements are texturally diverse occurring mostly as massive or laminated colloform layers of pyrite/marcasite intergrowth locally with high porosity, cementing detrital material ( Figure 6). In addition, minor euhedral pyrite associated with calcite in the pore space locally surrounds colloform pyrite/marcasite and framboidal pyrite. Two The detrital material mostly consists of quartz, feldspars (commonly intensely altered), muscovite, biotite, chlorite, rutile, and titanite, as well as subordinate clasts of mica schists, volcaniclastics, and microcrystalline silica reflecting the country rock of Paleochori Bay. The iron sulfide cements are texturally diverse occurring mostly as massive or laminated colloform layers of pyrite/marcasite intergrowth locally with high porosity, cementing detrital material ( Figure 6). In addition, minor euhedral pyrite associated with calcite in the pore space locally surrounds colloform pyrite/marcasite and framboidal pyrite. Two types of pyrite mineralization generally occur at Paleochori Bay (Figures 6 and 7), including colloform pyrite (Py I), which is associated with calcite and apatite (Figure 7), as well as massive pyrite (Py II) that is associated with barite, alunite, and amorphous silica (Figure 6c,g,h).
In the absence of colloform pyrite, calcite appears to be the first mineral cementing the sediments (Figure 6d), which is locally followed by gypsum mineralization. In contrast, in the absence of calcite, pyrite cement is followed by yellow-colored spherulitic aggregates of zeolites associated with clays and Fe-hydroxides. Silica commonly occurs as amorphous to cryptocrystalline quartz together with cinnabar forming microlayers (up to 3 µm) alternating with colloform pyrite/marcasite (Py I) bands (Figure 6e,f). In addition, native sulfur was identified in the chimney-and mound-like structures, which also occurs together with orpiment in microbial mats on the seafloor [46,48]. Microbe cells have been found attached to the sulfide precipitates (Figure 6i). The cells appear to be dividing and excreting a framework of honeycomb-like structures with conditioning polymers around them, both indicating that the microbes are thriving in the vent water and pyritic depositional environment. These microbes may contain high levels of arsenic [85].
types of pyrite mineralization generally occur at Paleochori Bay (Figures 6 and 7), including colloform pyrite (Py I), which is associated with calcite and apatite (Figure 7), as well as massive pyrite (Py II) that is associated with barite, alunite, and amorphous silica (Figure 6c,g,h).
In the absence of colloform pyrite, calcite appears to be the first mineral cementing the sediments (Figure 6d), which is locally followed by gypsum mineralization. In contrast, in the absence of calcite, pyrite cement is followed by yellow-colored spherulitic aggregates of zeolites associated with clays and Fe-hydroxides. Silica commonly occurs as amorphous to cryptocrystalline quartz together with cinnabar forming microlayers (up to 3 μm) alternating with colloform pyrite/marcasite (Py I) bands (Figure 6e,f). In addition, native sulfur was identified in the chimney-and mound-like structures, which also occurs together with orpiment in microbial mats on the seafloor [46,48]. Microbe cells have been found attached to the sulfide precipitates (Figure 6i). The cells appear to be dividing and excreting a framework of honeycomb-like structures with conditioning polymers around them, both indicating that the microbes are thriving in the vent water and pyritic depositional environment. These microbes may contain high levels of arsenic [85].

Hydrothermal Mineral Paragenesis
At Paleochori Bay, the sulfide mineralization is marked by two distinct mineral assemblages ( Figure 8). The formation of As-rich colloform pyrite/marcasite (Py I) is related to the occurrence of cinnabar micro-layers in the colloform textures also suggesting a contemporaneous deposition. Manganese-bearing calcite and apatite also form together with early pyrite/marcasite I, filling interstitial spaces. Massive pyrite II represents the second pyrite generation that contemporaneously formed with barite, gypsum, alunite-jarosite, hydrothermal silica, and Fe-hydroxides. Pyrite II could be additionally divided in subassemblages, as indicated by textural and morphological features shown in Figure 6c,g, and also supported by Figure 9. However, they are unified here for simplicity reasons, and will be studied in detail in a future paper related to their forming conditions. Although not observed, it is most likely that cinnabar also associates with massive pyrite II together with alunite.

Hydrothermal Mineral Paragenesis
At Paleochori Bay, the sulfide mineralization is marked by two distinct mineral assemblages ( Figure 8). The formation of As-rich colloform pyrite/marcasite (Py I) is related to the occurrence of cinnabar micro-layers in the colloform textures also suggesting a contemporaneous deposition. Manganese-bearing calcite and apatite also form together with early pyrite/marcasite I, filling interstitial spaces. Massive pyrite II represents the second pyrite generation that contemporaneously formed with barite, gypsum, alunitejarosite, hydrothermal silica, and Fe-hydroxides. Pyrite II could be additionally divided in sub-assemblages, as indicated by textural and morphological features shown in Figure 6c,g, and also supported by Figure 9. However, they are unified here for simplicity reasons, and will be studied in detail in a future paper related to their forming conditions. Although not observed, it is most likely that cinnabar also associates with massive pyrite II together with alunite.

Mineral-Chemistry
Pyrite is the dominant sulfide and accounts for >70 volume % of the total sulfide mineral content in all types of hydrothermal precipitates, followed by marcasite as the second most abundant sulfide. Electron microprobe analyses revealed a large variation in the chemical composition of pyrite/marcasite ranging from near stoichiometric compositions to elevated As (up to 3.2 wt %) and Mn (up to 0.4 wt %) in colloform pyrite/marcasite (Py I) and up to 1.7 wt % As and 0.6 wt % Mn in massive pyrite (Py II) (Table 3 and Figure 7). Cobalt and Ni occur in minor amounts reaching 0.1 wt % (Table 3). Arsenic-rich colloform pyrite/marcasite (Py I) exhibits a negative correlation between arsenic and sulfur in sample ML38A, suggesting that arsenic substitutes for sulfur in the pyrite lattice (Figure 9a) [90][91][92]. In contrast, a negative correlation exists between arsenic and iron in alunite-associated pyrite II of sample ML22B (Figure 9b). Although pyrite from sample ML11A is also alunite-associated, the negative Fe-As relationship is less well developed and As contents are generally very low (<0.2 wt %; Table 3). Cinnabar (HgS) typically occurs as micro-layers in colloform pyrite/marcasite I (Figure 6e,f), and hence elevated Fe contents in cinnabar (up to 7.0 wt %) are likely caused by an interaction of the electron beam with the pyrite/marcasite host.

Mineral-Chemistry
Pyrite is the dominant sulfide and accounts for > 70 volume % of the total sulfide mineral content in all types of hydrothermal precipitates, followed by marcasite as the second most abundant sulfide. Electron microprobe analyses revealed a large variation in the chemical composition of pyrite/marcasite ranging from near stoichiometric compositions to elevated As (up to 3.2 wt %) and Mn (up to 0.4 wt %) in colloform pyrite/marcasite (Py I) and up to 1.7 wt % As and 0.6 wt % Mn in massive pyrite (Py II) ( Table 3 and Figure  7). Cobalt and Ni occur in minor amounts reaching 0.1 wt % (Table 3). Arsenic-rich colloform pyrite/marcasite (Py I) exhibits a negative correlation between arsenic and sulfur in sample ML38A, suggesting that arsenic substitutes for sulfur in the pyrite lattice ( Figure  9a) [90][91][92]. In contrast, a negative correlation exists between arsenic and iron in aluniteassociated pyrite II of sample ML22B (Figure 9b). Although pyrite from sample ML11A is also alunite-associated, the negative Fe-As relationship is less well developed and As contents are generally very low (<0.2 wt %; Table 3). Cinnabar (HgS) typically occurs as microlayers in colloform pyrite/marcasite I (Figure 6e,f), and hence elevated Fe contents in cinnabar (up to 7.0 wt %) are likely caused by an interaction of the electron beam with the pyrite/marcasite host.

Interpretation on Mineral Paragenesis
The two distinct pyrite assemblages at Paleochori Bay reflect distinct precipitation conditions with respect to temperature, pH and fO2. The formation of As-rich colloform Figure 11. Chemical variation of alunite-jarosite solid solution minerals from Paleochori vents as a function of Na/(Na + K) relative to Fe 3+ .

Interpretation on Mineral Paragenesis
The two distinct pyrite assemblages at Paleochori Bay reflect distinct precipitation conditions with respect to temperature, pH and fO 2 . The formation of As-rich colloform pyrite/marcasite (Py I) (Figures 6-8) is typically related to disequilibrium conditions reflecting rapid crystallization due to abundant fluid-seawater mixing [9,10,95]. This early stage process cements the detrital material (Figure 6a-d) and insulates the fluid pathways from ambient seawater for subsequent precipitation in a more stable fluid environment, as reflected by later minor euhedral pyrite. These early processes are comparable to many other submarine hydrothermal systems and seem to be unaffected by water depth [9,10,39,95,96].
The alternation of pyrite/marcasite in the early colloform textures indicates fluctuations in the physicochemical composition of the hydrothermal fluids, as indicated by marcasite precipitation at lower pH and temperature relative to pyrite [97]. At Luise Harbor, Lihir Island, pH variations due to mixing between vapor-dominated fluids and seawater in variable proportions are likely to be responsible for the alternating deposition of pyrite and marcasite, where marcasite requires lower pH conditions than pyrite [30]. Alternatively, the pyrite/marcasite alternations at Paleochori Bay may be caused by fluid temperature fluctuations as a result of variable proportions of fluid-seawater mixing leading to marcasite deposition at lower temperatures [9,97]. Relatively large fluctuations of vent temperatures in Paleochori Bay occur at tidal frequencies and seismic events [98]. Abrupt fluctuations in fluid chemistry may be responsible for the alternation of the As-rich (up to 3.2 wt %) and As-poor layers in the colloform pyrite/marcasite I (Figure 7) indicating changes of the As content of the precipitating fluids.
The second pyrite generation (Py II) formed contemporaneously with alunite/natroalunite, jarosite/natrojarosite, hydrothermal silica, and Fe-hydroxides (Figures 6, 8 and 10). The latter may be related to a supergene oxidation of the Fe-sulfides by seawater [9,95,99]. The paragenetic sequence (Figure 8) and the observed As variations in pyrite relative to Fe and S ( Figure 9) probably reflect temporal or local changes in the chemistry of the hydrothermal fluids [9,100].

Formation of Arsenian Pyrite
Arsenian pyrite is a common constituent in epithermal systems and formation of arsenian pyrite may be the key to the formation of many Au deposits, in which pyrite can be an important and often undervalued host for Au [92,[100][101][102][103][104][105][106][107]. Two main types of arsenian pyrite are usually distinguished: (1) As 1− pyrite, in which As 1− substitutes for S reflected by a negative correlation of As and S typically found in Carlin-type and low-sulfidation epithermal deposits (e.g., [90]), and (2) As 2+/3+ pyrite, in which oxidized As species substitute for Fe leading to a negative As-Fe relationship commonly related to high-sulfidation epithermal conditions, due to their more oxidized character [92,100,106,108,109]. Alternatively, As can occur as amorphous As-Fe-S nanoparticles in As-rich pyrite, suggesting the presence of As 0 species [105]. Accordingly, the As-Fe-S ternary diagram provides information on the As incorporation and speciation in pyrite either as As 1− , As 2+/3+ , or As 0 [91,92,109]. The As-Fe-S ternary diagram (Figure 12), suggests that As 1− substitutes for S in colloform pyrite/marcasite (Py I) associated with calcite at Paleochori Bay, in accordance to the As vs. S scatterplot (Figure 9). The As-Fe-S ternary diagram also suggests the possible occurrence of massive As 3+ pyrite (Py II), as well as a substitution of Fe 2+ by Me 2+ leading to the incorporation of other divalent trace elements in pyrite II that is associated with alunite at Paleochori (Figure 12). The different As incorporation in pyrite II could reflect changes in forming conditions. However the exact substitution mechanisms and the specific geochemical conditions of formation at Paleochori Bay, cannot be resolved based on available data, and will be the subject of a future study.
The As 1− incorporation in colloform pyrite/marcasite I suggests more reducing lowsulfidation hydrothermal conditions (e.g., [110]) and is compatible with its calcite and apatite association in the chimney-and mound-like structures at Paleochori. Low sulfidation conditions are usually also characterized by lower fluid temperatures and neutral to alkaline pH conditions [111] indicating that marcasite precipitation at Paleochori Bay is thought to be due to temperature fluctuation, rather than being caused by decreasing to more acidic pH [97]. In contrast, massive As 3+ -pyrite II formation in the mounds, implies more oxidizing and lower pH conditions, as also reflected by its association with alunite (e.g., [111]). The lower pH conditions suggested for the pyrite II-related fluids could also explain the As depletion in these pyrites relative to the colloform As 1− pyrite ( Figure 12). Arsenic typically forms OHcomplexes in hydrothermal fluids [112,113] and the fact that acidic fluids have less OHthan their higher pH counterparts leads to a lower As solubility, which may be preserved by pyrite II with characteristically lower As contents (Figure 12).
Reducing and oxidizing fluid conditions at Paleochori Bay correlate with the occurrence of high-Cl liquid-and low-Cl vapor-dominated fluids, respectively. This indicates that high-Cl liquid-dominated fluids were richer in As 1− , Hg, and CO 2 leading to formation of As 1− pyrite, cinnabar, and calcite upon mixing with seawater at the seafloor. On the contrary, the low-Cl vapor-dominated fluids were more oxidizing, lower pH, and sulfaterich, leading to formation of As 3+ pyrite, barite, and alunite. This chemical distinction reflects the different reactions of vent brine-rich solutions with seawater and of vent vapor-rich solutions with seawater, resulting in a mineralization assemblage similar to low-and high-sulfidation epithermal systems, respectively.
Minerals 2021, 11, x FOR PEER REVIEW 16 of 25 and the specific geochemical conditions of formation at Paleochori Bay, cannot be resolved based on available data, and will be the subject of a future study. Figure 12. Ternary diagram of As-S-Fe in pyrite from Paleochori Bay (adopted from Deditius et al. [91] and Román et al. [109]). The different trends show substitution of (1) As for S (colloform As 1− pyrite I; ML38A), (2) As 3+ for Fe (massive As 3+ pyrite II, ML22B); and (3) substitution of divalent metals (Me 2+ ) for Fe (massive pyrite II, ML11A).
The As 1− incorporation in colloform pyrite/marcasite I suggests more reducing lowsulfidation hydrothermal conditions (e.g., [110]) and is compatible with its calcite and apatite association in the chimney-and mound-like structures at Paleochori. Low sulfidation conditions are usually also characterized by lower fluid temperatures and neutral to alkaline pH conditions [111] indicating that marcasite precipitation at Paleochori Bay is thought to be due to temperature fluctuation, rather than being caused by decreasing to more acidic pH [97]. In contrast, massive As 3+ -pyrite II formation in the mounds, implies more oxidizing and lower pH conditions, as also reflected by its association with alunite (e.g., [111]). The lower pH conditions suggested for the pyrite II-related fluids could also explain the As depletion in these pyrites relative to the colloform As 1− pyrite ( Figure 12). Arsenic typically forms OHcomplexes in hydrothermal fluids [112,113] and the fact that acidic fluids have less OHthan their higher pH counterparts leads to a lower As solubility, which may be preserved by pyrite II with characteristically lower As contents ( Figure  12).
Reducing and oxidizing fluid conditions at Paleochori Bay correlate with the occurrence of high-Cl liquid-and low-Cl vapor-dominated fluids, respectively. This indicates that high-Cl liquid-dominated fluids were richer in As 1− , Hg, and CO2 leading to formation of As 1− pyrite, cinnabar, and calcite upon mixing with seawater at the seafloor. On the contrary, the low-Cl vapor-dominated fluids were more oxidizing, lower pH, and sulfaterich, leading to formation of As 3+ pyrite, barite, and alunite. This chemical distinction reflects the different reactions of vent brine-rich solutions with seawater and of vent vaporrich solutions with seawater, resulting in a mineralization assemblage similar to low-and high-sulfidation epithermal systems, respectively.
We propose that the sulfate-dominated fluid and its associated pyrite-alunite mineralization is caused by the reaction of H2S originally dissolved in the deep fluid but separated upon boiling in the vapor phase with oxygenated seawater, producing H2SΟ4-rich acid-sulfate fluids in a steam-heated environment, as known from subaerial epithermal systems like those on Milos island (e.g., [57,111,114]). These acid fluids could leach and Figure 12. Ternary diagram of As-S-Fe in pyrite from Paleochori Bay (adopted from Deditius et al. [91] and Román et al. [109]). The different trends show substitution of (1) As for S (colloform As 1− pyrite I; ML38A), (2) As 3+ for Fe (massive As 3+ pyrite II, ML22B); and (3) substitution of divalent metals (Me 2+ ) for Fe (massive pyrite II, ML11A).
We propose that the sulfate-dominated fluid and its associated pyrite-alunite mineralization is caused by the reaction of H 2 S originally dissolved in the deep fluid but separated upon boiling in the vapor phase with oxygenated seawater, producing H 2 SO 4rich acid-sulfate fluids in a steam-heated environment, as known from subaerial epithermal systems like those on Milos island (e.g., [57,111,114]). These acid fluids could leach and dissolve primary iron-rich grains from the volcaniclastic components of the sediments, and result in precipitation of pyrite and alunite. This mechanism is already proposed by Pichler et al. [30] for Lihir Island shallow submarine vents. Both alunite and jarosite often form in steam-heated environments, and jarosite may be stable relative to alunite under exceptionally low pH (~1-2) and high Fe 3+ conditions [115,116]. These low pH conditions are not documented by the fluid data from Paleochori Bay (Table 1), which is likely due to seawater mixing during the precipitation process and/or during fluid sampling. Natrojarosite surrounding alunite at Paleochori Bay therefore suggests a rapid oxidation process of Fe 2+ to Fe 3+ upon mixing of the vapor-rich fluids with seawater.
Alternatively, the pyrite II-alunite association at Paleochori Bay could be the result of an acidic and oxidized magmatic fluid that was contributed by the seawater-dominated hydrothermal system of Paleochori Bay, which was suggested for several arc-related submarine hydrothermal systems like Brothers volcano, Kermadec arc [21,92], Panarea, Aeolian arc [117], and high-sulfidation epithermal deposits, such as Pueblo Viejo, Dominican Republic [106], and El Indio, Chile [118].
Although a minor contribution of magmatic fluids at Paleochori Bay cannot be excluded (as evidenced from the δ 11 B values of the vent fluids; [28]), a magmatic-hydrothermal scenario is most unlikely, because no SO 2 was detected in the gas phase [40,84], the δ 34 S value of H 2 S discharged to the seafloor demonstrates that any volcanic inputs were buffered by subsurface anhydrite veins [47], and finally, mineralogical evidence, such as enargite deposition is also lacking [102].

Effects of Biological Activity on Sulfide Deposition
Previous studies also suggested that biological activity may have enhanced sulfide deposition at the vent sites of Paleochori Bay [43,46,48,50,53,83,85]. We propose that a combination of microbial reduction of seawater sulphate (Figure 6i) together with sulfide deposition from H 2 S-rich fluids are responsible for the texturally diverse mineralization at Paleochori Bay ( Figure 6). This indicates a biotic and abiotic S cycle at Paleochori Bay, as also suggested by Houghton et al. [50] based on S and O isotope data of hydrothermal sulfides and pore water sulfates providing evidence for microbial sulfate reduction and sulfide oxidation. Furthermore, the S isotope composition of pyrite in seafloor sediments records a microbially influenced signature relative to the hydrothermal endmember [50].
The formation of framboidal pyrite ( Figure 6) together with orpiment-type As sulfides at Paleochori Bay has also been related to biological activity [48,119,120]. Framboidal pyrite formation requires a very high FeS 2 supersaturation in the hydrothermal fluids resulting in nucleation rates that significantly exceed crystallization rates [121]. Such conditions can be achieved in highly oxidizing milieus [122], which occur in the discharge area where mixing between hydrothermal fluids and oxygen-rich seawater is abundant. Beyond this, sulfate reducing microbes could also account for the formation of As 1− in the fluids from which As 1− pyrite precipitates. An oxidized vapor phase, on the other hand, that condensates in seawater in a steam heated environment (cf. Section 6.2) would rather carry As 3+ , from which As 3+ pyrite may precipitate ( Figure 12). However, whether biomineralization processes at Paleochori Bay, such as biogenic reduction of seawater-derived sulfate, prevail over abiotic hydrothermal fluid related precipitation processes is a matter of further investigation.

Comparison with Other Shallow and Very Shallow Marine Vent Sites in the Mediterranean and Elsewhere
Active shallow (e.g., <200 mbsl) seafloor (and sub-seafloor) sulfide deposition in the Mediterranean Sea occur in the Aeolian arc, Italy, and more especially at Baia di Levante, Vulcano Island, at a water depth of up to 18 m with temperatures <100 • C [123], as well as on near-shore Panarea Island (8-80 mbsl) with maximum temperatures for the submarine vents up to 135 • C [83,117,[124][125][126][127] (Figure 1). At Baia di Levante, shallow water fumarolic activity results in the formation of pyrite and marcasite cementing the sand grains in the surface sediments [117]. At Panarea area, gas and thermal water venting associated with sulfide precipitation occurs at several places among others: (a) at the so-called Black Point NW of Lisca Nera (water depth of 23.5 m, T < 137 • C), (b) in the NE of the Secca dei Panarelli (water depth of~55-85 m, T < 180 • C), and (c) N of Lisca Bianca (water depth of 25 m, T < 150 • C). Chimneys at Black Point consist of galena, pyrite, marcasite, sphalerite, and barite [126]. According to Monecke et al. [117], this site currently represents the shallowest site of polymetallic sulfide mineralization discovered in the worlds' oceans. Polymetallic sulfides NE of the Secca dei Panarelli include pyrite and marcasite, sphalerite and barite, and intergrowths of gypsum-anhydrite and of pyrite with alunite [117,124,125]. Finally, at the Lisca Bianca venting site, marcasite co-precipitates with alunite and opal [127]. Compared to the above-mentioned very shallow vent systems in the Mediterranean, the pyrite and marcasite depositional system at Paleochori Bay mineralogically strongly resembles those at Baia di Levante, Vulcano and Lisca Bianca, Panarea.
Punta Mita, Punta Banda, Bahía Concepción (Mexico) and Luise Harbor (Lihir, Papua New Guinea), in addition to Baia di Levante, Vulcano and Lisca Bianca, Panarea (see above), may also be considered as analogs to the Paleochori Bay system (Figure 1). At Punta Mita, sulfide deposition takes place from a mixture of liquid and gas that discharges on the seafloor at similar temperature (~85 • C) and water depth (<10 mbsl) leading to a comparable association of pyrite, cinnabar, and calcite. Similar to Paleochori Bay, the Lihir (e.g., Luise Harbor) hydrothermal field is characterized by abundant Fe sulfides deposited due to the interaction of phase-separated fluid with seawater and Fe-rich sediment at 10 to 50 m water depth and temperatures between 60 • C and 96 • C [30]. The northern Baja California hydrothermal system also vents similar temperature fluids (~102 • C) depositing pyrite and gypsum [36,37]. Generally, in these shallow-water hydrothermal systems, the seafloor sulfide deposition is limited because extensive sub-seafloor boiling or mixing with seawater results in sub-seafloor metal precipitation [19,35].

Enrichment in an Epithermal Suite of Elements and Cinnabar Deposition
The Paleochori Bay precipitates display enrichment in an epithermal suite of elements such as As, Hg, Sb, and Tl, and relatively low contents of Cu, Pb, and Zn, varying between those reported for other very shallow systems depositing iron sulfides (pyrite and marcasite) on the seafloor, such as Punta Mita, Punta Banda, Bahía Concepción, Luise Harbor, Baia di Levante, and Panarea (Lisca Bianca) ( Table 5 and text above). Arsenic, Hg, Sb, and Tl contents at Paleochori Bay are also comparable to those from arc volcanoes including Panarea (Black Point and Secca dei Panarelli). However, precipitates from arc volcanoes are much more enriched in Cu, Pb, and Zn compared to Paleochori due to their emplacement at greater depths and/or at higher temperatures [19,129].
Active cinnabar formation is a rare feature of submarine hydrothermal systems, and was only observed at a few other localities like Punta Mita, central Mexico [32], at Punta Banda and Bahía Concepción of northern Baja California, Mexico [34][35][36], and Bay of Plenty (Calypso), New Zealand [38,130]. Vent precipitates with anomalously high concentrations of Hg have been also reported in the Lihir hydrothermal system in Papua New Guinea [30], but without any report on cinnabar deposition (Table 5). However, although the Punta Mita, the Lihir field, and the northern Baja California systems discharge fluids at about 10-50 m water depth similar to Paleochori Bay, the vents at Bay of Plenty (~200 mbsl) are deeper and higher in temperature (up to 200 • C) [38,130] compared to Paleochori Bay.

Comparison with Other On-Shore Milos Mineralization and Sources of As and Hg in the Fluids
The Paleochori precipitates contain the first documented occurrence of actively formed cinnabar on the sea floor in the Aegean and the Mediterranean area and provide an important link between offshore hydrothermal activity and the mercury-(and As-) depositing mineralizing system on western Milos Island. Mercury and arsenic are enriched (up to 24 and 4350 ppm, respectively) in the northwestern part of Milos Island, where shallow submarine epithermal mineralization at Triades-Galana, Profitis Ilias, AgathiaKondaros, and Katsimouti-Vani is hosted by volcanic rocks [55,58,60] (Figure 2b). Paleochori Bay also shares similarities with As-and Hg-enriched subaqueous cherty silica present offshore at Kondaros-Agathia area where barite and natroalunite are intergrown with iron sulfides [58] ( Figure 2b).
The origin of As and Hg in the Milos mineralization can only be speculated: Wind et al. [131] proposed that enrichments in As, Tl, and Hg, which were observed in deposits that occur in the hanging wall of major detachment faults (e.g., Milos), may be related to leaching from associated supra-detachment sediments of the Miocene conglomeratic sequence. However, according to Katsouri et al. [128], the volcanic arc crust is probably the source of As for the Paleochori Bay hydrothermal system, and Price et al. [46] suggested that the extremely high As contents of the Paleochori Bay fluids are due to leaching of the As-rich metamorphic basement. By expanding to the above-mentioned sources, we suggest here that underlying sulfide mineralization at Paleochori Bay, similar to those present at western Milos island, may represent a potential source for As and Hg in the Paleochori fluids and therefore its precipitates.

Conclusions
(1) We herein report the discovery of active cinnabar-depositing hydrothermal vents in a submarine setting at Paleochori Bay, within the offshore southeastern extension of the Milos Island Geothermal Field, South Aegean Active Volcanic Arc. Active, shallow-water (2-10 m) and low temperature (up to 115 • C) hydrothermal venting at Paleochori Bay discharges CO 2 and H 2 S gas and liquid (pH~3.5-7.6) and has led to an assemblage of sulfide and alteration mineral phases in an area of approximately 1 km 2 . Hydrothermal edifices recovered from the seafloor are composed of volcaniclastic detrital material cemented by pyrite and marcasite in variable proportions. (2) Paragenetic relations indicate deposition of two distinct mineral assemblages: (1) colloform As 1− and Hg-bearing pyrite/marcasite (Py I) together with Mn-bearing calcite and apatite, and (2) massive As 3+ -rich pyrite (Py II) associated with alunite/natroalunitejarosite/natojarosite solid solution and barite. Amorphous silica occur throughout the paragenetic sequence. Mercury, in the form of cinnabar, occurs in up to 3 µm grains within arsenian pyrite (Py I) layers, usually forming distinct cinnabar-enriched micro-layers. (3) A negative correlation between As and S in pyrite I suggests that As 1− substitutes for sulfur, a feature consistent with reducing, higher pH, and low-sulfidation conditions prevailing during colloform pyrite/marcasite (Py I) and cinnabar deposition associated with calcite and apatite. Overgrowths of arsenian pyrite layers (up to 3.2 wt % As) onto As-free pyrite within the colloform pyrite/marcasite bands indicate fluctuation in the As content of the hydrothermal fluid. The deposition of As 3+ pyrite in association with alunite is consistent with low pH, low Cl, vapor-dominated, and higher fO 2 fluids as reflected by the positive correlation between As and Fe, where As substitutes for Fe in the octahedral site of pyrite. (4) The pyrite-bearing hydrothermal precipitates at Paleochori Bay are enriched in epithermal-type elements like As, Sb, Tl and Hg. In addition to boiling, mixing with seawater at the seafloor and biological activity are responsible mechanisms for this epithermal-like mineralization. (5) The Paleochori vents contain the first documented occurrence of cinnabar on the seafloor of the Aegean Sea and the Mediterranean area and provide an important link between offshore hydrothermal activity and the onshore mercuryarsenic mineralizing system on Milos Island. This study demonstrates that metal and metalloid precipitation in shallow-water continental arc environments is controlled by epithermal processes known from their subaerial analogues. Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available in minerals-1046095.pdf and the Supplementary Tables S1 and S2.