Shallow Sea Gas Manifestations in the Aegean Sea (Greece) as Natural Analogs to Study Ocean Acidification: First Catalog and Geochemical Characterization

Daskalopoulou, Kyriaki; D’Alessandro, Walter; Longo, Manfredi; Pecoraino, Giovannella; Calabrese, Sergio

doi:10.3389/fmars.2021.775247

REVIEW article

Front. Mar. Sci., 31 January 2022

Sec. Ocean Observation

Volume 8 - 2021 | https://doi.org/10.3389/fmars.2021.775247

This article is part of the Research Topic The Discovery of the Unknown Planet: The Ocean View all 17 articles

Shallow Sea Gas Manifestations in the Aegean Sea (Greece) as Natural Analogs to Study Ocean Acidification: First Catalog and Geochemical Characterization

$\r\nKyriaki Daskalopoulou,*$ Kyriaki Daskalopoulou^1,2*

Walter D’Alessandro³

Manfredi Longo³

Giovannella Pecoraino³

Sergio Calabrese^3,4

¹Institute of Geosciences, University of Potsdam, Potsdam Golm, Germany
²GFZ German Research Centre for Geosciences, Potsdam, Germany
³Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Palermo, Palermo, Italy
⁴Dipartimento di Scienze della Terra e del Mare, Università degli Studi di Palermo, Palermo, Italy

The concepts of CO₂ emission, global warming, climate change, and their environmental impacts are of utmost importance for the understanding and protection of the ecosystems. Among the natural sources of gases into the atmosphere, the contribution of geogenic sources plays a crucial role. However, while subaerial emissions are widely studied, submarine outgassing is not yet well understood. In this study, we review and catalog 122 literature and unpublished data of submarine emissions distributed in ten coastal areas of the Aegean Sea. This catalog includes descriptions of the degassing vents through in situ observations, their chemical and isotopic compositions, and flux estimations. Temperatures and pH data of surface seawaters in four areas affected by submarine degassing are also presented. This overview provides useful information to researchers studying the impact of enhanced seawater CO₂ concentrations related either to increasing CO₂ levels in the atmosphere or leaking carbon capture and storage systems.

Introduction

The concentration of carbon dioxide (CO₂) in the atmosphere is increasing mainly due to fossil fuel combustion and industrial processes. Since the beginning of the industrial revolution at the end of the eighteenth century, its level increased from about 280 ppm and exceeded the average yearly value of 413 ppm during the year 2021 (NOAA, 2021). Being one of the major greenhouse gases, such rapid increase has severe consequences on earth’s climate (IPCC, 2021). About one third of the anthropogenic CO₂ released into the atmosphere in the past two centuries has been taken up by the ocean (Sabine et al., 2004; Bindoff et al., 2019; Gruber et al., 2019). In aquatic systems CO₂ gas dissolves, hydrates and dissociates to form weak carbonic acid (Drever, 1997), and the pH is lowered according to the following reaction:

C O_{2 (g)} + H_{2} O \Leftrightarrow H_{2} C O_{3 (aq)} \Leftrightarrow H^{+} + H C O_{3}^{-} (1)

Current CO₂ emission rates exceed the buffering capacity of the oceans and cause a shift of marine carbonate chemistry and a decrease of pH that has been quantified in 0.1 units compared to the pre-industrial period (Haugan and Drange, 1996; Doney et al., 2009; Gattuso and Hansson, 2011; Jiang et al., 2019). Depending on different emission scenarios, models predicted that further CO₂ increase would cause an additional reduction of pH between 0.3 and 0.5 units by the end of the century (Caldeira and Wicket, 2005; Joos et al., 2011; Jiang et al., 2019). Business-as-usual CO₂ emission scenarios predict that atmospheric CO₂ will reach 750 ppm and pH levels will decrease to 7.8 by the year 2100 (Jiang et al., 2019). In addition to this, both surface temperature and heat content of the ocean have increased. Specifically at the ocean surface, temperature increased by 0.88°C on average from 1850–1900 to 2011–2020. Possible future scenarios anticipate that it will arrive at 0.86°C from 1995–2014 to 2081–2100 (IPCC, 2021). Similarly, ocean heat content increased by 0.28–0.55 YJ between 1971 and 2018 and will probably continue to increase until at least 2300 (IPCC, 2021). This projection applies also for low emission scenarios due to the slow circulation of the deep ocean.

Many studies evidenced that ocean acidification (OA) will exert significant and sometimes unexpected effects on marine ecosystems (Jiang et al., 2019). Because these changes decrease the saturation state of the carbonate ion (CO₃^2–) in seawater, organisms relying on calcification for growth or protection are assumed to be most severely affected (Doney et al., 2012). On the contrary, photosynthetic organisms, such as seagrass and algae, may benefit from the increasing pCO₂ which is an essential resource for their photosynthesis and survival (Fabricius et al., 2011; Koch et al., 2012; Russell et al., 2013). It should be mentioned that even though laboratory experiments documented the benefits of OA on seagrass growth, anthropogenic stressors might counterbalance positive effects of increased CO₂ and have likely blocked potential beneficial responses of OA (Koch et al., 2012; Doo et al., 2020). To face the problem of atmospheric CO₂ increase, apart from the most logical solution remaining the strong reduction of anthropogenic CO₂ emission, one remedy proposed is the geologic carbon sequestration. CO₂ capture and storage (CCS) systems concentrate and transfer liquid CO₂ into storage sites, including sub-seabed deep geological formations such as exhausted oil or gas reservoirs. This approach is considered promising, since technically feasible (IPCC, 2005), but as with all other human technologies, it is not exempt from drawbacks. One of these drawbacks is the possibility that the chosen reservoir is not perfectly sealed and undergoes CO₂ leakage (Monastersky, 2013). If these reservoirs are offshore, CO₂ leakages from CCS can drive strong local seawater acidification (Blackford et al., 2014), exceeding the values predicted by the worst scenario of climate change. Moreover, in the case of a CO₂ leak from a storage site, the gas will also acidify the pore water in the sediments surrounding the storage site (Millero et al., 2009). This may also increase the release of harmful elements from the sediments creating an additional negative impact on the marine environment (De Orte et al., 2014; Foo et al., 2018). Flohr et al. (2021) simulated CO₂ leakage from an offshore CO₂ storage site in the British sector of the central North Sea. The CO₂ release experiment (Flohr et al., 2021; Gros et al., 2021) lasted for 1 month and the authors illustrated that different approaches can detect, attribute and quantify the release.

Notwithstanding the increasing number of studies on the ecological consequences of OA and CCS leakage, many issues remain unexplored and, until now, the vast majority of them have been performed in laboratories mainly as short-term and univariate experiments (Cornwall and Hurd, 2016). To have a more realistic picture, experiments should be made on marine organisms in their natural ecosystems. In this sense, areas with natural CO₂ vents represent useful experimental locations to investigate the impact of OA on entire ecosystems (Hall-Spencer et al., 2008). Natural underwater vents of volcanic origin release gases composed mainly of CO₂ and may therefore represent a natural analog to study the impact of seawater acidification. The CO₂ vent areas are also perfect natural laboratories to study the impact of CO₂ leakage from CCS systems.

Few of these “natural laboratories” have already been used to study the effects of elevated CO₂ on ecosystems (Vizzini et al., 2010; Lauritano et al., 2015; Linares et al., 2015) sometimes evidencing the adaptation of complex ecosystems such as coral reefs (Golbuu et al., 2016; Teixidó et al., 2020).

These vent sites allow to study different habitats, including shallow coral reefs in Papua New Guinea, Japan, and Northern Mariana Islands (Enochs et al., 2016; Golbuu et al., 2016); seagrass meadows, macroalgae stands, and coralligenous in the Mediterranean Sea (Columbretes Islands, Spain—Linares et al., 2015; Ischia, Italy—Hall-Spencer et al., 2008; Vulcano, Italy—Boatta et al., 2013; Panarea, Italy—Rogelja et al., 2016; Methana, Greece—Baggini et al., 2014); as well as in the subtropical North East Atlantic reefs (La Palma, Canary Islands—Hernández et al., 2016). However, a larger representation of environments is needed to predict the biological and ecological consequences of OA.

Our study will give a first catalog of the gas vents within the Aegean Sea comprising a description of the areas. It will provide important information to researchers who study the impact of enhanced seawater CO₂ concentrations related to increasing CO₂ levels in the atmosphere or even to leaking CCS systems like (i) extension and morphology of the exhaling area; (ii) preliminary gas flux estimations and geochemical characterization of the gases; (iii) presence of possible confounding factors as for example emission of thermal waters and/or hydrogen sulfide, iron oxi-hydroxide flocculation. The geochemical characterization is based almost exclusively on literature data that are gathered together with some new results and are made available to the reader in Supplementary Table 1, while a description of the degassing sites and a rough estimation of the gas fluxes are available in Table 1. We also present unpublished data on pH and temperature measurements of surface seawaters in four areas affected by the submarine degassing (Supplementary Table 2).

TABLE 1

Table 1. List and general characteristics of the underwater degassing areas.

Study Area

The Aegean Sea (Figure 1) is located in the eastern Mediterranean and is a rift formed in a “backarc” setting. It is situated in the upper plate of the Hellenic subduction zone and west of Anatolia, where active tectonics is observed. In fact, the northern Aegean Sea is a part of the Eurasian plate and the boundary with the Aegean microplate is called the North Anatolian Trough (NAT). The latter is the continuation of the North Anatolian Fault Zone (NAFZ) and is a ∼300 km long system of tectonically active marine basins, up to 1,000 m deep (Le Pichon et al., 1987; Taymaz et al., 1991; Kreemer et al., 2004).

FIGURE 1

Figure 1. Map of the Aegean sea where the south Aegean active volcanic arc (SAAVA) is drawn with a dashed red line. Study areas are delimited in black squares, while Kolumbo is marked with a light blue triangle. Details on the study areas are found in section “Geological and Geochemical Description of the Submarine Degassing Areas and in situ Observations”.

The thinning of various tectonic units mainly emplaced during the Upper Cretaceous– Paleocene convergence–collision processes has resulted in the creation of the basin (Boccaletti et al., 1974; Robertson et al., 1991). It should be noted that the Hellenic subduction system was active since at least the Late Cretaceous, while the “backarc” rift was developed during Eocene-Early Miocene (Agostini et al., 2010; Jolivet et al., 2013). Despite the long-lasting formation of the Aegean basin (∼40 Ma), the extension rate is relatively low, so that the oceanic crust was not generated (Agostini et al., 2010).

Nowadays, the extension is seemingly localized around the Corinth-Patras rift (southern Greece), however; it was widespread during the Miocene (Sébrier, 1977; Mercier et al., 1979). Oligo-Miocene extensional metamorphic complexes outcrop in the Cyclades archipelago and the northern Aegean Sea (Lister et al., 1984; Gautier et al., 1993). The extension has proceeded from north to south, while the subduction front was retreating southward (Lauritano et al., 2015).

This geodynamically active regime is also characterized by intense seismic activity (Taymaz et al., 2007), by the presence of the south Aegean active volcanic arc (SAAVA) (Fytikas et al., 1984) and anomalous geothermal gradients (Fytikas and Kolios, 1979). Similar to other regions of intense geodynamic activity, extensive geogenic degassing takes place (Daskalopoulou et al., 2018a, 2019a) with gas manifestations being widespread both on land and underwater.

Geological and Geochemical Description of the Submarine Degassing Areas and in situ Observations

A total of 10 areas characterized by submarine degassing were documented and sampled along the Aegean Sea. Table 1 summarizes the general characteristics of the underwater sampling sites. A brief description of the degassing sites and in situ observations documented during the field campaigns are presented in this paragraph. Where possible, underwater filming allowed us to document and describe the degassing areas, and estimate the gas fluxes (as described in Supplementary Material).

Samothraki Island

The island of Samothraki is located at the NE part of the Aegean Sea of Greece (Figure 1) and belongs to the Circum Rhodope Zone (Kauffmann et al., 1976). It comprises five lithological units, which include: (i) low-grade metamorphic rocks (basement unit), (ii) an ophiolitic complex, (iii) a granite intrusion with biotite and a contact metamorphic event, (iv) Cenozoic volcanic rocks, and (v) Quaternary clastic sedimentary rocks (Kotopouli et al., 1989; St. Seymour et al., 1996). The rough relief with steep slopes characterizing the SSE part of the island is the result of the tectonic uplift movements, whereas natural weathering and erosion are responsible for the geomorphology (Pavlidis et al., 2005).

Sparse emission points characterized by ambient temperatures are found within the fisherman port of Therma (Figure 2A). The manifestations are rich in CH₄ (72.7% on average-Supplementary Table 1), while CO₂ is also present (23.8% on average-Supplementary Table 1). The flux of the bubbles is low (Table 1), and the gas manifestations, which are spread on an area of a few hundred m², are permanent.

FIGURE 2

Figure 2. (A) Map of the north section of Samothraki Island that shows the location and name of the emission point. The circled area corresponds to the degassing points of Pigi A and B, and Giotrisi at the hydrothermal system of Therma (Dotsika, 2012; Daskalopoulou et al., 2018a). (B) Map of the Kassandra Peninsula that shows the locations and names of the emission points. (C) Map of the Ilion area (Euboea) that shows the location of the sampling points.

This manifestation seems to have no relationship with the nearby (800 m south) on-land hydrothermal system of Therma (Figure 2A) that reaches emission temperatures up to 74°C (Dotsika, 2012) and whose bubbling gases have CO₂-rich composition (Daskalopoulou et al., 2018a). On the contrary, its CH₄-rich composition points toward a hydrocarbon reservoir like those that are widespread and exploited in the north Aegean Sea (Rigakis et al., 2001).

Chalkidiki Peninsula

Chalkidiki peninsula is located at NNW part of the Aegean Sea (Figure 1) and is a part of the Vardar-Axios Zone, and the Serbomacedonian and Rhodope Massif (from west to the east) (Kauffmann et al., 1976). The area of interest is situated at Kassandra peninsula in Vardas-Axios zone, with the latter being considered as a narrow fragment of the Serbomacedonian Massif (Kockel et al., 1977). Despite the various metamorphic facies, the zone mainly comprises granitic intrusions of Upper Jurassic age and carbonates of a similar age enclosing bauxite horizons (Mountrakis, 1985). The tectonic regime of the area is mainly influenced by Thermaikos gulf, which is the relic of an older larger elongated tectonic depression trending from NNW to SSE.

Underwater degassing takes place in two areas at the Kassandra Peninsula (Figure 2B). The first emission site is found in Agia Paraskevi in front of “Halkidiki Thermal Spa” hotel. The two neighboring main degassing points can be visually recognized from the hotel due to a lighter color with respect to the main sea body; they appear like large stains in the sea. One of the main degassing areas is next to the coast close to the thermal springs on land. The springs are at sea level within some small caves (Lazaridis et al., 2011). Bubbling gases sometimes occurred also inside the caves but are more widespread in the sea. In this area, the gases come up from a very shallow depth (<1 m) between the boulders that form the shore. Another degassing area is at some tens of meters from the coast. Here intense degassing occurs at about 5 m depth. Gas emission vents form some recognizable alignments and are probably accompanied by thermal water emission because most of the orifices are surrounded by white deposits. It is worth noting that H₂S is present both within the cave and in the two underwater degassing areas (Supplementary Table 1). The degassing is constant and the flux elevated (Table 1), especially at the area further from the coast. Here H₂S, although below the analytical detection limit (<10 ppm) in the sample collected close to the sea surface, could still be smelled in the atmosphere above the bubbling site. Considering that H₂S is highly soluble in water only where bubbling is very intense it may reach the surface after crossing 5 m of seawater.

The second degassing spot is found on the eastern coast of the Chalkidiki peninsula in Xyna. It is near the shoreline at the eastern end of a 3 km long sandy beach at the border with a private luxury resort. The gas flux is very low (Table 1). In correspondence to the bubbling site, on the beach (5 m from the shore) there is a small hypothermal spring (23°C) captured with a shallow well.

Euboea Island

The island of Euboea is found at the western part of the Aegean Sea (Figure 1) and is the second largest island of Greece. It consists of formations from the Sub-Pelagonian structural zone, while its southern part belongs to the Atticocycldic massif. Volcanism of Pliocene and Quaternary age took place in the area (Fytikas et al., 1976; Pe-Piper and Piper, 1989, 2002) contributing to the formation of geothermal fields. The major fault structures of the North Euboean Gulf, where the underwater vents are found (Figure 2C), comprise several segments of normal faults, trending about NW-SE and dipping NE with a total length of about 20–30 km (Pavlides et al., 2004).

Widespread underwater manifestations are found a few meters by the coast in the area of Ilion (Figure 2C). The widespread bubbling is constant and the flux, according to our estimation, is classified intense (Table 1). Hydrogen sulfide is present in minor concentrations (Supplementary Table 1), while the rusty color of the sediments suggests the existence of iron oxides deposition. It is worth noting that low pH values have been documented along the coast, with the lowest values being found in front of a high temperature and intensely degassing spring on land 10 m from the sea (Supplementary Table 2).

Methana Peninsula

Methana peninsula is located in Saronikos Gulf and represents the northwestern, still active part of the SAAVA (Figure 1). It belongs to the Atticocycladic zone and consists of Quaternary calc-alkaline volcanic rocks (andesites to dacites; Fytikas et al., 1986). Two different Pliocene-Quaternary tectonic domains within the Aegean plate are affecting the tectonic regime; a rapid N–S extension to the north and an E–W extension to the south (Jolivet et al., 2013).

Around the peninsula, two areas of submarine gas emissions have been recognized so far (Figure 3A). The first one is located in the northern part of the peninsula, where the “Pausanias baths” are found. These baths are associated with hydrothermal degassing emissions composed of almost pure CO₂ (D’Alessandro et al., 2008). They are found at about 3 m from the coast in a depth of 2 m. The seabed is formed of boulders and the flux is considered low (Table 1). Baggini et al. (2014) state that along the entire northern coast, the measured pH values of seawater are lower and highly variable with respect to their reference site. Local fishermen reported intermittent gas bubbling close to the coast about 3 km west of the Pausanias baths. This degassing site has not been confirmed and no gas sample has been collected. It may be related either to the historical eruption of 230 BCE, whose lava flow entered the sea in that area, or to the nearby active submarine volcanic edifice called Pausanias (Foutrakis and Anastasakis, 2018). In any case, the sea close to the Pausanias baths has been the site of many studies to investigate the effect of higher pCO₂ values on the local ecosystem (Baggini et al., 2014, 2015; Bray et al., 2014; Triantaphyllou et al., 2018; Patoucheas et al., 2021).

FIGURE 3

Figure 3. (A) Map of the Methana Peninsula that shows the locations and names of the emission points. (B) Map of the western section of Milos Island that shows the locations and names of the emission points. (C) Map of the Kammeni Islands that shows the locations and names of the emission points.

The second submarine vent is situated in the eastern part of the peninsula in a small bay called by the locals “Thiafi bay.” This area is nearly 300 m long and is demonstrating on the beach widespread alteration from recent fumarolic activity. The alteration is particularly evident at its northern and southern ends, where it is expressed as native S and sulfates (alunite, gypsum, and alunogen) (Rahders et al., 1997). CO₂ fluxes on land are sometimes elevated and account for the whole area for about 500 t/a (D’Alessandro et al., 2008). The underwater gas vents show very low fluxes and are visible at shallow depths (<3 m) (Table 1) close to the coast (from the shoreline to distances of a few tens of meters). Few of these vents have a distinguishable orifice, and only some of them are associated with obvious deposits (possibly amorphous silica), whereas no thermal anomaly has been found at the gas vents (D’Alessandro et al., 2008).

In both degassing sites, CO₂ is the prevailing gas component (up to ∼98%), while H₂S was documented only in Thiafi (Supplementary Table 1).

Milos Island

Milos Island is found in the center of SAAVA (Figure 1) in the convergence zone between the African and the Aegean plates. It belongs in the Atticocycladic zone and comprises Upper Pliocene submarine and Upper Pleistocene to Holocene submarine-to-subaerial calc-alkaline volcanic domes, lavas, and pyroclastic deposits (andesites, dacites, and rhyolites; Fytikas et al., 1986; Stewart and McPhie, 2006). The youngest (Upper Pleistocene) volcanic activity is located in the volcanic centers of Fyriplaka in the south and Trachilas in the north, which are also present-day exhaling areas (Fytikas et al., 1986). However, the seabed around the island is hydrothermally very active. This is particularly noticeable in eight areas, located along the eastern part of the island (Dando et al., 1995; Cronan and Varnavas, 1999; Daskalopoulou et al., 2018b; Ivarsson et al., 2019). Gas emissions have been studied in the following sites (clockwise starting from the north): Mandrakia, Voudia, Paleochori, Spathi Bay, Agia Kyriaki, DEH (Kanavas), Skinopi (Figure 3B).

The little center of Mandrakia stands just on the prominent central part of a gulf oriented to the north. In the eastern part of the gulf, just close to some boulders, in a sandy seabed with a large Posidonia grassland, some sparse emissions are present tens of meters away from the shoreline, at a depth of about 3 m (Table 1). The vents stand exactly on the direction of the impluvium present on land and emit mainly CO₂ (about 98%).

Voudia bay has two submarine hydrothermal vents. One of them is situated few meters away from the shoreline aligned with altered rocks on the beach, while the other one is found in the southern part of the bay. According to Megalovasilis (2020), the temperatures of the vents range from 28 to 78°C. Daskalopoulou et al. (2019b) collected samples from the emanations spot and reported a CO₂-rich gas (96.6%), with minor contents of H₂S (3,100 μmol/mol; Supplementary Table 1). Evidence of sulfur yellow-concretions were noted on the walls of some degassing vents (Supplementary Figure 2). While sampling, the same authors have observed a relatively low flow (Table 1), which was afterwards confirmed by the documentation of Megalovasilis (2020).

Paleochori Bay is an 800 m long bay with apparent fumarolic activity at its eastern and western parts. Numerous intensively degassing seeps of elevated temperatures (up to 122°C; Dando et al., 1995; Khimasia et al., 2020) occur in the area. The emanating gases are rich in CO₂ (up to 93%; Supplementary Table 1) with minor enrichments in N₂ (up to 14%; Supplementary Table 1). H₂S is also present in concentrations up to 3.5% (Supplementary Table 1). The area is characterized by Fe- and S-alteration products (Cronan and Varnavas, 1999; Baltatzis et al., 2001; Voudouris et al., 2021). According to various authors (Dando et al., 1995; Yücel et al., 2013; Khimasia et al., 2020), the sea-bottom of the bay is dominated by areas of white and brown bacterial mat (Supplementary Figures 2C,D). Areas of gray/yellow sand with encrustations and of light brown sand with numerous burrows are also present. Godelitsas et al. (2015) suggested that the reddish or yellow-colored sediment patches in the center of the white mats might have been caused by elemental sulfur (Supplementary Figure 2C) and arsenic sulfides precipitation. It is worth noting that the gas ascends from rock fissures at a depth of about 2–3 m, resulting in bacteria-dominated outlets known as “White Smoker” (Supplementary Figure 2E). This degassing area is distinguished by the highest gas fluxes estimated in this study (about 1.8 L/min) (Table 1). It is important to note that the hydrothermal fluids release to the seawater huge quantities of potentially toxic elements. For example, several studies revealed concentrations of the order of thousands of μg/L of As, Ba, Fe, and Mn, and up to 1 μg/L of Hg (Price et al., 2013; Roberts et al., 2021) leading to the formation of hydrothermal precipitates rich in these and other (Sb, Tl) elements (Voudouris et al., 2021).

Spathi bay is located in the south-eastern sector of the island and extends for about 400 m. The coastline is characterized by pebbles and is bordered to the west by an imposing promontory and to the east by striking stacks. We do not have detailed information on the underwater emissions as the gas sample was kindly taken by colleagues.

Agia Kyriaki is located in the southern sector of Milos, in a bay that stretches for about 500 m. The coast is predominantly sandy and the seabed is characterized by alternating sandy areas and reefs. The only underwater gas manifestation in this area is found at less than 10 m from the shoreline, where a few isolated bubble-trains with a low flux outcome from a substrate of rock blocks, at a depth of about 2–3 m (Table 1). H₂S was present in concentrations lower than 10 μmol/mol (Supplementary Table 1).

DEH is located in the sea along Kanavas coast close to the power plant of the Hellenic Public Power Corporation. The degassing area is about 400 m² with a CO₂ output estimated at 1.06 t/d (Daskalopoulou et al., 2018b). The numerous and widespread gas manifestations have low temperatures with H₂S being less than 16 μmol/mol (Supplementary Table 1). The sea bottom is sandy and very shallow (about 1–2 m) and clear alignments of high-flux degassing vents from fissures were noted (Supplementary Figure 2A and Table 1). The environment is disturbed by many human activities, i.e., a shipyard 500 m north, abandoned salt flats about 200 m south, and the power plant less than 100 m away which discharges its cooling water about 50 m from the bubbling area.

Skinopi is a small bay with a 100 m long pebble beach on which stand some characteristic houses and small jetties for fishermen’s boats. It is located 1,500 m west of Adamas, the main port of the island. Close to the coast (tens on meters) there are many bubbling areas. The degassing area is shallow (<2 m) and not very active (low fluxes) (Table 1), while the sandy and boulders sea bottom is mostly covered by Posidonia grasslands. The gases are mainly composed of CO₂ (>82%), and H₂S is undetectable (Supplementary Table 1).

Santorini Island

Santorini volcanic complex is found in the center of the SAAVA (Figure 1). It comprises the islands of Thera, Thirasia, Palea Kammeni, Nea Kammeni, and Aspronisi, and belongs to the Atticocycladic zone. The complex consists of volcanic rocks (mainly pumice and glass) and metamorphic formations (mainly marbles and phyllites) (Druitt et al., 1989, 1999; Oikonomidis and Pavlides, 2017). The evolution of the volcanic centers is associated with two NE-SW faults; the Kammeni Line (Heiken and McCoy, 1984; Druitt et al., 1989, 1999; Parks et al., 2013) and the Columbo Fault Zone (Druitt et al., 1989, 1999; Mountrakis et al., 1998). In addition to this, Tzanis et al. (2020) demonstrated that both volcanism and the shape of the volcanic center are controlled by the tectonics.

Submarine gas vents are located in Palea and Nea Kammeni islets (Figure 3C). The emission point in Palea Kammeni is in the bay of Agios Nikolaos and is located close to the coast in the eastern part of the island. Two emissions have been documented in the Nea Kammeni island. One is called Agios Giorgios and is found on the western side of the island, while Irinia is on the eastern side where most of the island visitors are disembarked. Both Agios Nikolaos and Agios Giorgios are CO₂ dominated (Chiodini et al., 1998; Tassi et al., 2013; Daskalopoulou et al., 2018a) and present H₂S content up to 26 and 415 μmol/mol, respectively (Supplementary Table 1). The outlet temperatures of the two points are 36 (Agios Nikolaos) and 40°C (Agios Giorgios) (Böstrom and Widenfalk, 1984; Dotsika et al., 2009). On the other hand, Irinia shows a mixed CO₂-N₂ composition (Daskalopoulou et al., 2018b; Tarchini et al., 2019). All three sites are in protected coves where the seawater is heavily stained by iron oxides due to the input of the hydrothermal fluids. In the deepest parts of the coves, iron concentrations in seawater exceed 1 mg/L (up to 13.7 mg/L; Smith and Cronan, 1983). The iron, solubilized by the low pH CO₂-rich hydrothermal fluids, becomes oxidized on mixing with seawater and precipitates accumulating abundantly at the bottom of the exhaling areas (Smith and Cronan, 1983).

Kolumbo

Kolumbo is a submarine volcano found 7 km northeast off Santorini island (Figure 1). It is a high-temperature hydrothermal field (Sigurdsson et al., 2006) characterized by numerous vents of CO₂-rich gases (<97%) and fluids of ∼220°C (Carey et al., 2011). Despite its proximity to Santorini, volcanological and petrological evidence suggest the existence of two separate plumbing systems beneath the two volcanic edifices (Francalanci et al., 2005; Dimitriadis et al., 2009; Kilias et al., 2013). In addition to this, Rizzo et al. (2016) has documented a more than 85% mantle contribution for He (the highest across SAAVA). Gases collected in Kolumbo are CO₂ dominated (>98% on average; Carey et al., 2011; Rizzo et al., 2016, 2019; Supplementary Table 1).

Kos Island

Kos Island is located in the eastern part of the SAAVA (Figure 1). It comprises alluvial deposits with greenschists and flysch in the north, lacustrine and terrestrial deposits of the Pliocene age in the central part, while tuffs and ignimbrites of the Quaternary age are found in the south (La Ruffa et al., 1999). Faults of WNW-ESE and NE-SW orientation seem to control the tectonic evolution of the island and to be related to extensional processes and volcanic activity during the Pleistocene and Pliocene (Lagios et al., 1998). In the area, four submarine degassing centers have been recognized (Daskalopoulou et al., 2019b). Paradise beach and Kefalos are located in the western part of the island, while Therma and Agia Irini in the eastern (Figure 4A).

FIGURE 4

Figure 4. (A) Map of the southern sector of Kos Island that shows the locations and names of the emission points. (B) Map of the NE sector of Nisyros Island that shows the locations and names of the emission points. (C) Map of Gyali Island that shows the locations and names of the emission points.

The submarine emissions of Paradise beach are rich in CO₂. H₂S is always below detection limit (Supplementary Table 1). The vents are found at approximately 20 m from the coast at 1–1.5 m depth. They are widespread and having elevated gas flows (Table 1).

The marine area of Kefalos at the SW of the island is interested by diffuse degassing, with a lot of bubble streams mainly concentered just to the east of the harbor area. Hundreds of little vents that emit trains of little bubbles are present in a sandy seabed at a depth of few meters.

Submarine gases of Therma present similar chemical characteristics to the gases of Paradise. However, they are found by the coast and are characterized by elevated temperatures (up to 45°C).

Two degassing vents have been recognized in the sea in front of the Agia Irini church. Despite the vicinity of the emission points, their prevailing gas components differ significantly. One (Agia Irini 1) is rich in CO₂ (>95%), while the other (Agia Irini 2) is rich in N₂ (34–99%). At both sites, H₂S is generally below detection limit (Supplementary Table 1).

Nisyros Island

Nisyros Island is found at the eastern end of the SAAVA (Figure 1) and is a quiescent active stratovolcano with intense fumarolic activity that is generated by the presence of a high enthalpy geothermal system (Marini et al., 1993). It belongs to the Atticocycladic unit and consists of Quaternary volcanic rocks and alternations of lava flows, pyroclastic deposits and lava domes. The island has an area of 47 km² and forms a truncated cone with a base diameter of 8 km and a 4 km wide central caldera (Hunziker and Marini, 2005), known as the Lakki Caldera. Numerous, mostly sub-vertical faults crosscut the island and the caldera. The vertical offsets for the majority of these faults decrease from the caldera rim toward the coast, where they practically disappear (“scissors-type” faults); something that evidences their association with volcano-tectonic effects (Stiros, 2000).

Two points of submarine vents have been recognized in Nisyros (Figure 4B). Lies and Katsouni are two gas manifestations rich in CO₂ (Supplementary Table 1). They are found few meters from the coast in < 2 m depth, where the sea bottom is mostly covered by boulders and pebbles. Both sites are characterized by low temperatures and medium to low gas fluxes.

Gyali Island

Gyali Island is located between the Islands of Kos and Nisyros (Figure 1). The small island is uninhabited except by workers for the extraction of pumice and occasional tourists visiting the picturesque bays. It consists of a thick rhyolitic pumice succession to the south (Gyali pumice breccia and overlying units) and rhyolitic lava to the north. These two formations are separated by an isthmus, which is found in the center of the island. According to Allen and McPhie (2000), the pumice breccia of Gyali is the result of a submarine phreatomagmatic eruption.

Three submarine vents have been recognized in the island (Figure 4C), in the area where the fault zones are located (Daskalopoulou et al., 2021a). The gas manifestations are rich in CO₂, with H₂S being found in minor content (Supplementary Table 1). Intense degassing activity is observed in both Gyali West and Gyali South, where the emissions are widespread at some tens of meters from the shore. In both areas the sea bottom is sandy, but while the former is at about 10 m depth the latter is shallower (1–3 m depth). The third site is less active in terms of degassing and is found very close to the shore. The bubbles rise at shallow depth (about 1 m) between large boulders.

Geochemistry of Submarine Gas Vents

In the current study, a total of 122 data from submarine gas manifestations are presented. This dataset comprises both literature (Chiodini et al., 1998; Shimizu et al., 2005; Kyriakopoulos, 2010; Carey et al., 2013; Tassi et al., 2013; Rizzo et al., 2016, 2019; Daskalopoulou et al., 2018a, b, 2019b, 2021a; Tarchini et al., 2019) and 7 unpublished results from submarine degassing areas distributed in ten areas of the Aegean Sea; seven of which belong to the SAAVA. Gases were collected using the inverted funnel method (Figure 5). Sample IDs, coordinates, references, gas flux estimations and their chemical and isotopic content are found in Supplementary Table 1. pH and temperature data of seawater affected by the submarine degassing are also presented in Supplementary Table 2. Details on the sampling techniques and the laboratory methods used are found in Methods section in Supplementary Data Sheet 1.

FIGURE 5

Figure 5. Gas sampling with the use of the inverted funnel method at Kefalos (Kos Island). Note that the funnel was constructed at the mechanical laboratory of Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Palermo (INGV-Palermo). Additional photos regarding the sampling technique are provided in Supplementary Figure 1. Photo courtesy of SC.

No samples plot close to the atmospheric point (Kipfer et al., 2002) in the CO₂-N₂-O₂ ternary diagram (air in Figure 6A), thus excluding important air contaminations. The vast majority of the samples present N₂/O₂ ratios higher than the ratios of air saturated waters (ASW; Kipfer et al., 2002), suggesting that the atmospheric component deriving from meteoric recharge has been modified by microbial or inorganic redox reactions that took place in the subterranean circuit. In their majority, gases are rich in CO₂, while few samples have N₂ or CH₄ as the prevailing gas species (Figure 6B). Geographically, gases rich in N₂ (Santorini, Nisyros, Kos, Methana) are distributed in the SAAVA and just after SAAVA to the north (Euboea), while CH₄ prevails in the gas vents of north eastern Aegean Sea (Samothraki) (Daskalopoulou et al., 2018a, 2019a). Figures 6A,B show the occurrence of dissolution processes, with the latter being likely responsible for the CO₂ loss and the enrichment of the less soluble gases (CH₄, O₂, and N₂).

FIGURE 6

Figure 6. Ternary plot of (A) CO₂-N₂-O₂ and (B) CH₄-N₂-CO₂. Processes impacting the gases are drawn with an arrow. The abbreviation “ASW” stands for air saturated water. Values of air and air saturated water (ASW) after Kipfer et al. (2002). Literature data from Chiodini et al. (1998), Shimizu et al. (2005), Kyriakopoulos (2010), Carey et al. (2013), Tassi et al. (2013), Rizzo et al. (2016, 2019), Daskalopoulou et al. (2018a, b, 2019b, 2021a), and Tarchini et al. (2019).

An important atmospheric contribution for He is noticed for the gases of Samothraki Island as they plot close to the atmospheric point (Figure 7A, after Sano and Wakita, 1985). As expected, the vents located at the SAAVA present an enhanced MORB-type mantle contribution arriving up to ∼90% (Kolumbo), while the gases of non-volcanic areas show a more crustal origin for He (up to ∼95%). Figure 7B (Sano and Marty, 1995) reveals a mixed mantle-limestone origin for C for the great majority of the gases. The contribution of the organic sediment is relatively negligible. Some gases present CO₂/³He ratios that fall below the Mantle field, indicating CO₂ loss. This is likely due to the dissolution of CO₂ in water or to the precipitation of carbonates (Kanellopoulos, 2012; Winkel et al., 2013; Stefánsson et al., 2016, 2017; Kanellopoulos et al., 2017). A change in the flow as observed from Daskalopoulou et al. (2019b) for the sample of Kos presenting a strong decrease in δ¹³C-CO₂ may have also contributed to the dissolution processes. Methane for most submarine gases plots in the field ascribed to abiogenic hydrocarbons emitted from volcanic-geothermal systems (Figure 7C, after Etiope and Schoell, 2014). Contribution from biogenic sources cannot be excluded due to the wide range (from ∼4 to 244) of CH₄/(C₂H₆+C₃H₈) ratios (Bernard et al., 1978). Some gases exhibit high δ¹³C-CH₄ and δ²H-CH₄ values. This points to either organic or inorganic CH₄ oxidation processes. It is worth mentioning that isotope fractionation for organic and inorganic processes follow different fractionation paths (details in Daskalopoulou et al., 2018a, 2019b). Low δ²H-CH₄ values of Milos can be explained by non-equilibrium fractionation of CH₄-H with either H₂O or H₂ (Botz et al., 1996). Biogenic origin is attributed to CH₄ for the area of Samothraki. In particular, low δ¹³C-CH₄ and δ²H-CH₄ values as well as intermediate CH₄/(C₂H₆+C₃H₈) ratio indicate mixing between thermogenic and microbially-derived gases. The latter has been attributed to CO₂-reduction by Daskalopoulou et al. (2018a) on the basis of carbon isotope fractionation factor between coexisting CO₂ and CH₄ (Whiticar et al., 1986).

FIGURE 7

Figure 7. Binary plot of (A) R/R_A vs. ⁴He/²⁰Ne of the Hellenic gas emissions. The mixing lines between Atmosphere and Mantle and between Atmosphere and Crust are also plotted. Dashed lines represent mixing between atmosphere and end-members with different percentages of mantle contribution (after Sano and Wakita, 1985); (B) CO₂/³He vs. δ¹³C-CO₂. The composition for Sediments, MORB-like Mantle and Limestones end-members are, as follows: δ¹³C-CO₂ = −30‰, −5‰ and 0‰ and CO₂/³He = 1 × 10¹³, 2 × 10⁹ and 1 × 10¹³, respectively (after Sano and Marty, 1995); and (C) modified Schoell binary diagram (Etiope and Schoell, 2014) between δ²H-CH₄ and δ¹³C-CH₄ ratios for the Aegean submarine gas discharges. Slopes of biogenic and abiogenic oxidation of CH₄ are also plotted. Literature data from Chiodini et al. (1998), Shimizu et al. (2005), Kyriakopoulos (2010), Carey et al. (2013), Tassi et al. (2013), Rizzo et al. (2016, 2019), Daskalopoulou et al. (2018a, b, 2019b, 2021a), and Tarchini et al. (2019).

The CO₂ loss is more evident in the binary plots of Figure 8. The positive correlation between He and N₂ indicates the impact of the CO₂ dissolution on the gases (Figure 8A). He and N₂, as well as CH₄ (Figure 8B), are less soluble respect to CO₂, hence the strong solubility difference between the gases in the marine environment may has resulted in CO₂ loss (Reid et al., 1987). In fact, D’Alessandro et al. (2014) showed that when a gas mixture ascends through non-saturated waters, solubility contrasts might enrich the less soluble gases. This is specifically applicable in reduced gas upflow conditions (Daskalopoulou et al., 2021b).

FIGURE 8

Figure 8. Binary plot of (A) He-N₂ and (B) CO₂/He-CH₄. Processes impacting the gases are drawn with an arrow. Literature data from Chiodini et al. (1998), Shimizu et al. (2005), Kyriakopoulos (2010), Carey et al. (2013), Tassi et al. (2013), Rizzo et al. (2016, 2019), Daskalopoulou et al. (2018a, b, 2019b, 2021a), and Tarchini et al. (2019).

Figure 9 further evidences the impact of solubility-related processes on gas composition. The elevated CO₂/H₂S ratios of the gases in Kos, Gyali and Santorini demonstrate the interaction between magmatic gases and hydrothermal systems (Supplementary Table 1). This process, known as magmatic scrubbing (Symonds et al., 2001), occurs when ascending gases encounter any aquifer interposed between the source magma stored at depth and the surface (Di Napoli et al., 2016). The more water-soluble gas species dissolve due to the gas-water-rock interactions, thus modifying the composition of the primary magmatic gas phase. In their great majority, gases collected along SAAVA are poor in CH₄. This scarcity evidences the lack of underlying organic-rich source rocks in these areas. Few gases of Milos, Kos and Santorini are virtually enriched in CH₄ due to loss of CO₂ by dissolution as evidenced also in Figure 8B. Finally, some gases of Chalkidiki, Milos, and Euboea plot close to the atmospheric point, further demonstrating the CO₂ loss due to dissolution processes, but also indicating some contribution of the atmospheric component.

FIGURE 9

Figure 9. Ternary plot of CO₂-H₂S-CH₄. Processes impacting the gases are drawn with arrows. Values of air after Kipfer et al. (2002). Literature data from Chiodini et al. (1998), Shimizu et al. (2005), Kyriakopoulos (2010), Carey et al. (2013), Tassi et al. (2013), Rizzo et al. (2016, 2019), Daskalopoulou et al. (2018a, b, 2019b, 2021a), and Tarchini et al. (2019).

In order to further demonstrate the impact of water-gas interactions on the gas content, two samples for each site were collected at the degassing centers of Agia Paraskevi (Chalkidiki) and Agia Irini 2 (Kos). One of the two samples was taken in the usual manner from the emission site at sea bottom (for details see “Material and Methods” section in Supplementary Data Sheet 1), while the second at the sea surface after the gas bubbles have risen through the entire water column. Results evidence that sea bottom samples have CO₂ as the major component. In one case they show the presence of some H₂S. The superficial samples have much lower CO₂ concentrations, H₂S always below the detection limit, and become enriched in N₂, O₂, He, and CH₄ (Figure 10 and Supplementary Table 1). These, sometimes very strong changes, can be explained by two processes that drive the gas exchanges between the rising bubbles and the seawater: mixing between two end-members and fractionation due to different solubility. The first process accounts for the virtual enrichment of the less soluble gases of hydrothermal origin (He and CH₄) with respect to CO₂ (Figure 10A). The second process is responsible for the decrease of CO₂ and H₂S (hydrothermal end-member) and the increase in O₂ and N₂ (Air-saturated seawater end-member) (Figure 10B). The extent of the changes suffered by the ascending gases depends on many conditions, which are mainly temperature, area of the interaction surface and interaction time (i.e. distance to be covered from the sea bottom to the surface). In the case of dry gases (no water vapor) the temperature is generally that of seawater because even if the emitted gases are hot they rapidly equilibrate with the seawater temperature due to the water/gas mass ratio and the thermal inertia of water. Both the interaction surface area and interaction time strongly depend on gas flux, bubble dimension, and depth of the water column. Higher gas fluxes, greater bubble dimensions, and lower emission depths all reduce gas exchange between bubbles and water, limiting the changes in gas composition. In the case of the above mentioned sites, the high gas flux and shallow depth of Agia Paraskevi prevents strong compositional changes like those registered in the site of Agia Irini 2. For example, while in the first case 78% of the initial CO₂ content arrives at the sea surface in the second case less than 1% does (Supplementary Table 1).

FIGURE 10

Figure 10. Binary plot of (A) CO₂-CH₄ and (B) CO₂−−N₂. Lines in (A) indicate the course toward less soluble gas species, while the line in (B) indicates mixing from the deeper to the shallower gas sample and then to the atmospheric end member. Values of air after Kipfer et al. (2002).

Changes on water characteristics were documented for the areas of Kanavas and Paleochori (Milos), Therma (Kos), and Ilion (Euboea) (Figures 11A–D). These areas comprise widespread degassing vents that are characterized by elevated CO₂ contents (>90%) and presence of H₂S (Supplementary Table 1). pH transects that were performed along the coast revealed pH values lower than the value of average seawater (Supplementary Table 2). Such lowering of the pH is driven by the dissolution of CO₂ in seawater.

FIGURE 11

Figure 11. pH transects for the areas of (A) Paleochori and (B) Kanavas at Milos, (C) Therma at Kos, and (D) Ilion at Euboea. Points of particular interest are marked with an arrow. The morphology of the seashore at Paleochori (A) was diverse when the measurements took place.

Daskalopoulou et al. (2018b) identified areas (on-land and in the sea) of intense degassing and anomalous CO₂ flux at both Kanavas and Paleochori that correspond to the sites where the lowest pH values are found (Figures 11A,B). These areas are connected to the main fault structures recognized in the island (Stewart and McPhie, 2006) that have likely taken part in the volcanic activity of Milos, acting as pathways for the ascending magma (Kokkalas and Aydin, 2013) and the uprising gases (Dando et al., 1995; Daskalopoulou et al., 2018b). Especially for Paleochori Bay, Aliani et al. (2004) and Khimasia et al. (2021) mapped minor fault structures in the bay with the latter using microbial mats and high temperatures (Khimasia et al., 2020) for fault identification. It is worth mentioning that the temperature profiles by Khimasia et al. (2021) took place at hydrothermal vents within the bacterial mats (Figure 11A). Intense on land and submarine degassing is also identified in the island of Kos at Therma with the CO₂ upflow being completely addressed to the hydrothermal component (Daskalopoulou et al., 2019b). Megalovasilis (2020) suggested that the low pH values (Figure 11C) are the result of either fluid mixing with seawater in the substrate or water-gas interactions. It is important to note that the lowest pH values in Milos and Kos are found in emission points characterized by high temperatures (Figure 12). This, combined with the waters’ enriched SO₄-Cl content (Li Vigni et al., 2021), evidences the impact of hydrothermal activity on the systems.

FIGURE 12

Figure 12. Binary plot of pH-T for the areas where the transects took place.

In the area of Ilion, the low pH values are concentrated in front and east of a thermal (>62°C; D’Alessandro et al., 2014; Li Vigni et al., 2021) spring at the road (Figure 11D). The spring, as well as many of the thermal water emissions in Euboea, is connected to one of the tectonic structures in the border with the Sperchios Basin-Evoikos Gulf graben. The spring water ascents from the deep and hot geothermal system of the area (Li Vigni et al., 2021), and according to D’Alessandro et al. (2014) is affected by the Quaternary volcanic system. The volcanic impact though is also evident at the eastern side of the coast, where the degassing is diffuse. There, rocks present a rusty color, which is indicative of the emission of iron-rich (probably thermal) groundwater in the area. It is worth mentioning that the mobility of Fe and other trace metals (e.g., Cd, As, Pb) is enhanced by reducing conditions (Tarasov et al., 2005). The main driver is the dissolution of reactive gases (e.g., CO₂, H₂S, H₂) that results in intense rock leaching (Aiuppa et al., 2000).

Overview

Submarine degassing may have an impact on marine environments through ocean acidification. Hence, there is a necessity to study and better understand the ocean and its components. The current work reviews all known submarine gas manifestations of the Aegean Sea and summarizes the geochemical processes taking place in the individual areas.

All in all, degassing occurs in both volcanic and non-volcanic areas and is associated with the complex tectonics of the individual systems. Carbon dioxide is the dominant gas species for most vents and is often related to volcanism, geothermal energy, and elevated heat flow (Fytikas and Kolios, 1979). On the other hand, sites where CH₄ is the dominant gas species are likely related to hydrocarbon reservoirs (Rigakis et al., 2001).

The isotope signatures of He for gases found in SAAVA yield an important mantle contribution, while a dominant crustal origin characterizes gases in non-volcanic areas. Carbon dioxide derives from mixed mantle-limestone sources for most samples and in cases exhibits unimportant contributions from organic sediment sources. Methane is attributed to abiogenic hydrocarbons discharged from volcanic-geothermal systems. Inorganic and organic CH₄ oxidation processes resulting in isotope fractionation have also been identified (Daskalopoulou et al., 2018a). Only at Samothraki, where it is the main gas species, CH₄ is of biogenic origin.

The impact of water-gas-rock interactions on the initial gas phase is evident as soluble gas species dissolve in the water. This results in their depletion and the consequent enrichment of less soluble gas species. This phenomenon was also noticeable while comparing the composition of the gases at the emission point on the seafloor and the sea surface after its rising through the entire sea column. In addition to gas content variations, pH transects were performed in 4 sites. These are characterized by volcanic/geothermal activity, have CO₂ as the dominant gas species, and presented lower pH respect to the average marine value.

Even though the impact of gases on marine flora and fauna was not investigated in the current study, it shouldn’t be disregarded. Various researchers (e.g., Boatta et al., 2013; Price and Giovannelli, 2017; Aiuppa et al., 2021; Caramanna et al., 2021) have already highlighted that shallow marine vents can provide us with an accessible and economic way to investigate the effects of CO₂ on the whole marine ecosystems. In addition to this, parameters like the presence of light, wave action, tides, the input of meteoric water, salinity variations, etc., can significantly influence the geochemistry of the vents and the microbial diversity and distribution (Giovannelli and Price, 2018).

It is important to note that this is a preliminary catalog of shallow submarine vents found in the Aegean Sea. Springs found in tectonic structures on-land close to the coast (Li Vigni et al., 2021) characterized by strong degassing (Daskalopoulou et al., 2019a) and soil alterations (D’Alessandro et al., 2020) can be good indications of nearby submarine degassing underscoring that the catalog has still to be completed. Nevertheless, we aim that this study will initiate further research on the OA in Greece and in other countries. Following the identification of new emission sites and the quantification of gas flow, research should move towards a bio-, hydro-, and geochemical monitoring direction. As a next step, the anthropogenic input has to be taken into consideration. Understanding and defining the impact of both geogenic and anthropogenic processes affecting the various spheres in systems such as the Mediterranean, will result in the improvement of not only the carbon cycle knowledge but also of the dynamics and vulnerability of individual systems. This research and our group invite researchers from related disciplines to use multiple approaches and investigate other aspects of this problem. Hence, we make results accessible as electronic supplements and we aim to publish them as stand-alone datasets with their own doi.

Author Contributions

SC and KD contributed to the conception and design of this study. ML, GP, and SC collected the samples by diving. ML contributed to the laboratory analyses. KD and WD’A wrote the first draft while all authors contributed to manuscript revision, read and approved the submitted version. All authors were involved in the sampling campaigns.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

This manuscript comprises literature data collected principally by the authors’ working group. We kindly acknowledge all the friends and colleagues who helped us either in the field or with precious information about the sampling sites. We would like to thank once again: the LAVA Mining and Quarrying SA Company that allowed us to access Gyali Island, offered us lunch and accompanied us around the island with their boat looking for underwater gas emanations (Mr. Diamantis); Sabina Morici who helped us with the pH measurements at Paleochori and DEH (Milos Island) and Roy Price and Thomas Pichler, who collected gas samples at Spathy Bay (Milos Island); Marine warrant officer Ilias Simadakis who gave us the permission to collect samples at Agia Irini (Kos Island), the owner of “Kardamena Watersports Center” Mr. Nikos Nikolakopoulos who gave us one of his boats, and Mr. Yannis Limperis who accompanied us to the points of interest at Agia Irini. For the analyses made at the laboratories of the INGV of Palermo, we are indebted to the heads of the laboratories and technicians: G. Capasso, F. Grassa, M. Martelli, Y. Oliveri, A. Rizzo, F. Salerno, A. Sollami, and M. Tantillo. We are grateful for the insightful comments of the two reviewers and of the editor HR that helped us to improve the manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2021.775247/full#supplementary-material

References

Agostini, S., Doglioni, C., Innocenti, F., Manetti, P., and Tonarini, S. (2010). On the geodynamics of the Aegean rift. Tectonophysics 488, 7–21. doi: 10.1016/j.tecto.2009.07.025