Estimation of CO2 release from thermal springs to the atmosphere

Abstract

Introduction Geodynamically active regions have long been recognized as areas of anomalous Earth degassing [Irwin and Barnes, 1980]. Areas found at plate boundaries are characterized by seismic, volcanic and geothermal activity as well as ore deposition. These processes are enhanced by the circulation of hydrothermal fluids in the crust, which transport volatiles from the deep crust or mantle to the surface [King, 1986]. Kerrick and Caldera, [1998], were the first to indicate the significant contribution of the CO2 degassing by extensional tectonic and hydrothermal activity in metamorphic belts during the Phanerozoic. Moreover, further studies concerning gas emissions from diffuse degassing tectonic structures on various geological regimes suggested in their majority elevated CO2 concentrations [Klusman, 1993]. In fact, it is worth noting that the estimated global hydrothermal CO2 flux from subaerial geothermal environments can be comparable to that of direct volcanic discharges [Kerrick et al., 1995; Seward and Kerrick, 1996]. Study Area The back-arc geothermal fields of Greece include, among others, the Tertiary sedimentary basins of both Sperchios Basin and north Euboea, which are located in central Greece. Their tectonic activity contributes in crust thinning [Papadakis at al., 2016 and references therein] and elevated heat flow values [Fytikas and Kolios, 1979]. These geothermal anomalies due to the tectonic activity and the geological and volcanic regime are expressed as hot springs (Ypatis, Psoroneria, Thermopyles and Kamena Vourla in Sperchios Basin and Edipsos and Ilion in north Euboea). Tectonics of central Greece seems to be of particular interest as major fault structures are found in the area. Sperchios Basin was formed through the activity of WNW-ESE trending faults [Georgalas and Papakis, 1966; Marinos et al., 1973], whilst the Sperchios tectonic graben itself is considered to be the extension of the North Anatolia strike-slip fault. Moreover, in the north Euboean Gulf, the major fault structures are those of the Atalanti Fault Zone (AFZ) that consist of several segments of normal faults, trending about NW-SE [Pavlides et al., 2004]. Materials and Methods Six groups of springs (Ypatis, Psoroneria, Thermopyles, Kamena Vourla, Edipsos and Ilion) were investigated in this study. Bubbling gases were sampled using an inverted funnel positioned above the bubbles and stored in glass flasks equipped with two stopcocks until analysis. Samples for dissolved gas analyses were collected in glass vials and were sealed underwater. In the laboratory, the concentrations of He, H2, H2S, O2, N2, CO2 and CH4, on the samples were analysed by an Agilent 7890B gas chromatograph with Ar as carrier. The total CO2 emitted through bubbling was measured at 6 different pools (Psoroneria, Psoroneria 2,Thermopyles, Leonidas, Kamena Vourla and Ilion), whereas at other springs (Koniavitis-Sperchios Basin, Edipsos-Damaria and Edipsos-Thermopotamos) an estimation of the release was made by visual inspection. The CO2 fluxes were measured using the floating chamber method [Mazot and Bernard, 2015] that was equipped with a portable fluxmeter (WEST Systems, Italy) based on the accumulation chamber method as suggested by Chiodini et al., [1998]. The flux data were processed with both the Graphical Statistical Approach (GSA) and the Stochastic Simulation Approach (SSA), with the latter being based on the algorithm of sequential Gaussian simulation [Deutsch and Journal 1998; Cardellini et al., 2003]. Zonal Statistics on the final CO2 flux maps was obtained using the ArcMap 10.3 (ESRI) Spatial Analyst tool and were used to estimate the total CO2 output to the atmosphere. Results and Conclusions Carbon dioxide is the prevailing gas species for the great majority of the under investigation sites, with only gases collected in the area of Kamena Vourla (Kamena Vourla and Koniavitis) being rich in N2. The total bubbling CO2 emission from the pools to the atm