Original Research ARTICLE
The Geothermal Resource in the Guanacaste Region (Costa Rica): New Hints From the Geochemistry of Naturally Discharging Fluids
- 1Department of Earth Sciences, University of Florence, Florence, Italy
- 2CNR-IGG Institute of Geosciences and Earth Resources, Florence, Italy
- 3Observatorio Vulcanológico y Sismológico de Costa Rica, OVSICORI-UNA, Heredia, Costa Rica
- 4Sezione di Palermo, Istituto Nazionale di Geofisica e Vulcanologia, Palermo, Italy
The Guanacaste Geothermal Province (GGP) encompasses the three major volcanoes of northern Costa Rica, namely from NW to SE: Rincón de la Vieja, Miravalles, and Tenorio. The dominant occurrence of (i) SO4-rich acidic fluids at Rincón de la Vieja, (ii) Cl-rich mature fluids at Miravalles, and (iii) HCO3--rich and low-temperature fluids at Tenorio was previously interpreted as due to a north-to-south general flow of thermal waters and a magmatic gas upwelling mostly centered at Rincón de la Vieja, whereas Miravalles volcano was regarded as fed by a typical geothermal reservoir consisting of a highly saline Na-Cl aquifer. The uniformity in chemical and isotopic (R/Ra and δ34S) compositions of the neutral Cl-rich waters suggested to state that all the thermal discharges in the GGP are linked at depth to a single, regional geothermal system. In this scenario, the thermal manifestations related to Tenorio volcano were regarded as a distal and diluted fluid outflow. In this study, a new gas geochemical dataset, including both chemical and isotopic (δ13C-CO2 and R/Ra) parameters of fluid discharges from the three volcanoes, is presented and discussed. Particular attention was devoted to the Tenorio thermal manifestations, since they were poorly studied in the past because this area has been considered of low geothermal potential. The aim is to provide insights into the magmatic-hydrothermal fluid circulation and to verify the spatial distribution of the heat fluid source feeding the fluid manifestations. According to this new dataset, CO2, i.e., the most abundant dry gas in the fluid manifestations, is mostly produced by limestone, whereas the mantle CO2 contribution is ≤3.3%. Strongly acidic gas compounds from magma degassing were absent in the discharged fluids, being scrubbed by secondary processes related to prolonged fluid-rock interactions and mixing with shallow aquifers. Our results only partially confirm the previously depicted model, because the geochemical and isotopic features (e.g., relatively high concentrations of temperature-dependent gases and high R/Ra values) shown by fluids seeping out from the southern sector of Tenorio volcano are more representative of medium-to-high enthalpy volcanic systems than those typically occurring in distal areas. This implies that the geothermal potential in the south of the GGP is higher than previously thought.
Actions aimed to the evaluation of the sustainable potential of geothermal fluid reservoirs can provide fundamental information for minimizing the financial risk related to the exploitation activity of such natural resources. Geochemical prospections, coupled with geophysical measurements, are considered useful tools to investigate the chemical-physical features of hot, deep aquifers at relatively low cost (e.g., Giggenbach, 1991), especially for developing countries poor in conventional energy sources, such as carbon, oil, and natural gases. The geothermal potential in the Guanacaste Geothermal Province (GGP), northern Costa Rica, has been investigated since the early 1970's (Fergurson and Alfonso, 1977; Gardner and Corrales, 1977). There are three stratovolcanoes in this region, namely from NW to SE: Rincón de la Vieja, Miravalles, and Tenorio. In the last centuries, phreatic and phreatomagmatic eruptions frequently occurred at Rincón de la Vieja (e.g., Barquero and Segura, 1983; Boudon et al., 1996; Kempter et al., 1996), whereas at Miravalles and Tenorio (Figure 1) no historical eruptions were recorded (Mainieri et al., 1985; Alvarado, 2000). GGP shows numerous thermal fluid discharges, e.g., fumaroles, steam-heated soils, low-to-high temperature SO4- and Cl-rich springs, and CO2-rich bubbling pools (Giggenbach and Corrales, 1992; Gherardi et al., 2002; Tassi et al., 2005), which, as commonly occurs in a volcanic environment, are related to a hydrothermal-magmatic source. The spatial distribution of neutral, geothermal-type Cl−-rich water discharges, locally known as salitrales, and acidic rich springs produced by interaction of meteoric water and hydrothermal gases, was interpreted as an evidence for the presence of a single contiguous geothermal reservoir beneath the GGP (Giggenbach and Corrales, 1992). The Cl− concentration and temperature of this deep aquifer were estimated at ~3,500 mg/L and ~250°C, respectively (Giggenbach and Corrales, 1992; Gherardi et al., 2002). A similar temperature was also found in fluids from other exploited geothermal systems of Central America, such as those located in Salvador (Ahuachapan and Berlin), whose Cl− concentrations at the beginning of the exploitation where up to 9,000 and 11,000 mg/L, respectively (D'Amore and Mejia, 1999). The temperature hypothesized for the GGP reservoir, supported by calculations based on water isotopes and the contents of the main anions (Cl− and ), were consistent with those measured on the fluids discharged from the geothermal wells along the SW slope of the Miravalles volcano (Gherardi et al., 2002). According to this scenario, the area named Las Pailas (on the NW flank of Rincón de la Vieja volcano; #6 in Figure 1) and the inner zone of the Miravalles caldera were regarded as the areas most affected by deep fluid upwelling, whereas the thermal manifestations located within the SE portion of GGP, corresponding to Tenorio volcano, were considered as distal fluid outflow (Figure 1). These indications gave a pulse to the development of the geothermal fluid exploitation at Miravalles (Vallejos-Ruiz et al., 2005), as well as the prospection at Rincón de la Vieja (Molina and Martì, 2016, and references therein) where the geothermal fluid exploitation is still at the preliminary stage.
Figure 1. Map of the Guanacaste region with the location of sampling sites. Green square: Rincón de la Vieja (R); blue square: Miravalles (M); red square: Tenorio (T).
In the present study, the chemical and isotopic (δ13C-CO2 and R/Ra) composition of the GGP fluid discharges, with a special focus on those related to Tenorio volcano, is reported. At the light of this robust geochemical database, integrated with data from literature, the primary sources and secondary chemical-physical processes controlling the fate of the hydrothermal fluids are investigated. A critical revision of the hydro-geochemical conceptual models proposed for GGP by the previous authors is carried out, providing new insight into the geothermal potential of this region.
Regional Setting and Volcanic Activity
Costa Rica has a complex geology, related to the presence of three tectonic plates (Caribbean, Cocos and Nazca), the Panama microplate and an uncertain number of tectonic terranes (Alvarado and Cárdenes, 2016). These slabs are separated by a zone of diffusive deformation named Central Costa Rica Deformed Belt (CCRBD) and are bound to the southwest by the Middle America Trench (MAT; Montero et al., 2017). The magmatic arc, which is located ~150 km northeast of MAT, consists of (i) two active volcanic ranges, namely Cordillera Central (CC) and Cordillera de Guanacaste (CG), and (ii) two extinct ranges, i.e., Cordillera de Tilarán-Aguacate (CAg) and Cordillera de Talamanca (CT) (Alvarado and Cárdenes, 2016).
The Cordillera de Guanacaste (NW Costa Rica) was formed about 0.6 Ma and comprises four stratovolcanic complexes (1,500–2,000 m high): Orosí-Cacao, Rincón de la Vieja-Santa María, Miravalles-Zapote, and Tenorio-Montezuma (Chiesa et al., 1994). A series of E-W-trending lineament were identified between Orosí-Cacao and Rincón de la Vieja volcanoes, whereas secondary lineaments showing a different trend (NNE-SSW, N-S-, and NNW-SSE) were recognized along the eastern flank of Miravalles volcano, as well as at Rincón de la Vieja and Tenorio volcanoes (Chavarría et al., 2010) (Figure 1).
Rincón de la Vieja is the largest volcano of the Cordillera de Guanacaste, with an estimated volume of 130 km3 (Carr, 1984). It is a complex andesitic stratovolcano consisting of an elongated, arcuate NW-SE trending of nine coalescing pyroclastic cones, spreading over an 8-km long axis (Chiesa et al., 1994). The maximum elevation (Santa Maria cone) is 1916 m a.s.l. The last plinian eruption occurred 3.5 ka (Alvarado et al., 1992), whereas the historical activity was characterized by the occurrence of frequent phreatic and phreato-magmatic eruptions from the Active Crater (Tristan, 1921; Barquero and Segura, 1983). Since 1991, 20 phreatic eruptions, producing ash deposits and debris-flows, occurred (OVSICORI-UNA, 1998), while a transition to a phreato-magmatic activity took place in 1995 (OVSICORI-UNA, 1995). Intense phreatic activity was recorded in August-September 2011, February-April 2012, September 2014, February 2015, March 2016, February-October 2017 (OVSICORI-UNA, 2011, 2012, 2014, 2015, 2016, 2017).
Miravalles is the highest volcano of the Cordillera de Guanacaste (2028 m a.s.l.). This andesitic stratovolcano pertains to the Miravalles-Zapote complex, whose activity started with the edification of the Guayabo volcano (1.17–1.0 Ma; Alvarado and Gans, 2012). The latter was destroyed by the following explosive activity (1.5 and 0.6 Ma) giving rise to the formation of the 15 km wide Guayabo caldera (Chiesa et al., 1992). The latter hosts Cabro Muco volcano (0.57–0.28 Ma), in the eastern sector, and Zapote and Miravalles volcanoes (<0.28 Ma) in the northern sector (Chiesa et al., 1992; Alvarado and Gans, 2012). Tephrochronology data suggested that an important eruptive event occurred at about 5050 BC, whilst local witnesses reported a volcanic event on the 14th of September 1946 in the SW flank of Miravalles.
Tenorio-Montezuma, in the south-easternmost part of the Cordillera de Guanacaste, is a basaltic-andesitic volcanic complex, consisting of NNW-ESE-oriented volcanic structures, including pyroclastic cones and two twin-craters (Tenorio and Montezuma) (Alvarado, 1993). An andesitic dome field (Bijagua) grew over a relatively old volcanic center in the northern flank of Tenorio (Chiesa et al., 1994). Little is known about the geology of this intensely forested area. The andesitic and basaltic lavas of Monteverde (2.17–1.92 Ma) are the basement of the Tenorio edifice that developed during two main eruptive phases, dated 0.74–0.54 and 0.37–0.26 Ma, respectively (Alvarado and Gans, 2012). The occurrence of a caldera in the southern flank was hypothesized by ICE (Istituto Costaricence de Electricidad) and ENEL (Italian Board of Electricity) (ICE, 1990), although Alvarado and Gans (2012) suggested that the depression should be regarded as a northernmost branch of the Arenal graben. No historical volcanic activity has likely occurred at Tenorio-Montezuma. Notwithstanding, a legend, not supported by scientific investigation, speaks of an eruption that occurred in 1816 (Alvarado, 1993).
Previous Geothermal Prospection and Exploitation
A pre-feasibility study carried out in 1976 in the Miravalles area suggested the presence of a geothermal reservoir with temperature ~240°C (Gardner and Corrales, 1977). From 1977 to 1986, new exploratory activities, including the drilling of eight wells, provided the first insights into the geochemical features of the liquid-dominated system, i.e., neutral pH and a typical Na+-Cl− composition (Giggenbach and Corrales, 1992; Gherardi et al., 2002). The chemical and isotopic composition of fluid discharges from the GGP sampled between 1982 and 1988, allowed the construction of a hydro-geochemical model of the hydrothermal system (Giggenbach and Corrales, 1992). A geophysical prospection carried out by ICE and ELC (Electroconsult) between 1989 and 1991 (ICE, 1991), revealed areas of great interest for geothermal fluid exploitation along the slopes of Rincón de la Vieja, Miravalles, and Tenorio volcanoes. From 1994 to 2004, 5 geothermal power plants were installed in the Miravalles area, amounting for a maximum sustainable capacity of 163.5 MW (Sànchez-Rivera and Vallejos-Ruiz, 2015). Then, ICE focused its attention on the Rincón de la Vieja area, especially on the westernmost flank of this volcano (Las Pailas and Borinquen hydrothermal areas), where the occurrence of a number of fumaroles, silicified rocks, soils and thermal springs clearly indicated a significant thermal anomaly (Molina and Martì, 2016). At Las Pailas, a first geothermal plant (Unit I; 41.5 MW) was installed in 2011, whereas a second unit (55 MW) is planned to be placed to the east of Unit I, in 2018. At Borinquen, the construction of two units (55 MW each) is planned for 2023-2024 (ICE, 2014; Sànchez-Rivera and Vallejos-Ruiz, 2015), following the indication of a recent study (Molina and Martì, 2016) suggesting the occurrence of an aquifer at a temperature ≥220°C located at <3 km depth.
In the Tenorio area, preliminary geothermal prospection, carried out by ICE in 1987-1991, recognized an area of geothermal interest on the southern flank of the volcano (Mayorga, 2009). According to this indication, in 1999, two exploratory wells, PGT-1 and PGT-2, were drilled down to 2,473 m and 1,345 m depth, respectively (Mayorga, 2009). Unfortunately, the low temperature of the discharged water (<160°C), the low permeability and the proximity of a protected area (Tenorio Volcano National Park) discouraged the development of geothermal exploitation (Mayorga, 2009). New preliminary geochemical prospecting activities, carried out since 2008 in a 198 km2 wide area that includes the northern Tenorio Volcano National Park section and a zone located east of Bijagua (Figure 1), did not provide promising results (Sànchez-Rivera and Vallejos-Ruiz, 2015).
Materials and Methods
Gases were collected during 4 campaigns (April 2000, March 2005, February 2009, and February 2017) from thermal manifestations located in the surroundings of Rincón de la Vieja, Miravalles, and Tenorio volcanoes (Figure 1).
Gases from fumaroles were collected using a sampling line consisting of a 1 m-long titanium tube (Ø = 2.5 cm) and pyrex glass dewared pipes connected to pre-weighted and pre-evacuated 50-mL glass Thorion®-tapped flasks filled with 20 mL of a 4N NaOH and 0.15M Cd(OH)2 suspension (Montegrossi et al., 2001). Bubbling gases were sampled using a funnel up-side-down positioned above the bubbles and connected to the soda flasks using a silicon tube. Water vapor and CO2, dissolved in the alkaline solution, whereas H2S formed insoluble CdS. Low-solubility gas species (N2, O2, CO, H2, He, Ar, CH4, and light hydrocarbons) were trapped in the flask headspace. Dry gases (for the analysis of the 13C/12C ratios in CO2 and He isotopes) were collected using a water-cooled condenser connected to the soda flask sampling line.
Inorganic (N2, Ar, O2, H2, He, and CO) and organic (CH4, C2H6, C3H8, and C6H6) gases were analyzed at the Laboratory of Fluid Geochemistry of the University of Florence (Italy) by gas chromatography (GC; Shimadzu 15A, Shimadzu 14, and Thermo Focus). The liquid and the solid phases of the suspension were separated by centrifuge, to analyze for CO2 in the form of by acidimetric titration (AT; Metrohm Basic Titrino). The CdS precipitate was dissolved by oxidation with H2O2 in order to analyze H2S as by ion chromatography (IC; Methrom 761). The analytical error for GC, AT and IC analyses was <5%.
The analysis of 13C/12C of CO2 (expressed as δ13C-CO2 ‰ vs. V-PDB) was carried out at the laboratory of Stable Isotopes of CNR-IGG (Pisa, Italy) using a Finningan MAT252 mass spectrometer after standard extraction and purification procedures of the gas mixtures (Evans et al., 1988) performed at the Department of Earth Sciences of Florence (Italy). Carrara and San Vincenzo marbles, as well as international NBS18 and NBS19 standards, were used to estimate the external precision. The analytical error and the reproducibility were ±0.05 and ±0.1‰, respectively.
The R/Ra (where R is the 3He/4He measured ratio and Ra is the 3He/4He ratio in the air: 1.39 × 10−6; Mamyrin and Tolstikhin, 1984) and 4He/20Ne ratios were determined at the INGV laboratories in Palermo (Italy) by separately introducing He and Ne into a split-flight-tube mass spectrometer (GVI Helix SFT) after performing standard purification procedures (Rizzo et al., 2015). The analytical error was ±1%. The measured R/Ra values were corrected for air contamination using the 4He/20Ne ratios (Poreda and Craig, 1989), as follows:
where r = (4He/20Ne)air/(4He/20Ne)meas, the (4He/20Ne)air ratio being that in the atmosphere (0.318; Ozima and Posodek, 1983) and the (4He/20Ne)meas ratio that measured in the gas sample.
The chemical compositions of gases collected in 2000, 2009, 2015, and 2017 from the GGP are reported in Table 1, whereas the location of the fluid discharges is shown in Figure 1. Data from literature of some of these thermal fluid discharges from Rincón de la Vieja and Miravalles that were sampled in 1984 and 1987 (BOR1 and PA; Giggenbach and Corrales, 1992), as well as in 1998, 1999, and 2001 (BOR1 and PA; Tassi et al., 2005), are also reported. The dataset also includes the chemical composition of gases collected from different sites within GGP in 1987 and 1990 (MZ, TM, and GB; Giggenbach and Corrales, 1992), 1998 (GB; Gherardi et al., 2002), 1999, 2001, and 2002 (BOR2, BOB, LC, PAB, PAV, VO2, SM, and AC; Tassi et al., 2005), as well as that of gases from exploitation wells of the Miravalles geothermal plant (PGM01, PGM05, PGM10, PGM11, PGM12, PGM17, PGM31, and PGM49; Gherardi et al., 2002).
Table 1. Outlet temperature (°C), chemical and stable isotope (δ13C-CO2, in ‰ vs. V-PDB; Rc/Ra; R/Ra) composition of fumaroles from GGP.
The thermal manifestations consist of bubbling and boiling pools with temperatures ranging from 21 to 99°C. The chemical composition of the dry fraction of the fluids (the concentrations of steam were mostly controlled by the liquid phase at the surface) was dominated by CO2 (up to 996 mmol/mol), followed by N2 (from 2 to 53 and mmol/mol). The concentrations of H2S and CH4 ranged widely (from 0.075 to 17 and from 0.071 to 20, respectively), while those of H2 were relatively high (up to 4.8 mmol/mol) with the exception of those (from 0.0013 to 0.0085 mmol/mol) measured in the gases from the NE flank of Tenorio volcano (BO, RR, RC, PM, PC and BA; Figure 1). Oxygen and Ar varied from 0.019 to 0.79 and from 0.015 to 0.44 mmol/mol, respectively, whereas CO and He were ≤0.054 and ≤0.0068 mmol/mol, respectively. The Miravalles thermal discharges were characterized by relatively high concentrations of ethane (C2H6; up to 67 μmol/mol), propane (C3H8; up to 11 μmol/mol), and benzene (C6H6; up to 75 μmol/mol) with respect to those of the other thermal manifestations (≤8.2, ≤2.0, and ≤2.4 μmol/mol, respectively). Minor amounts (≤0.21 μmol/mol) of propene (C3H6) were also measured.
δ13C-CO2 and Rc/Ra Values
The δ13C-CO2 and Rc/Ra values of gases from the GGP are listed in Table 1. Data from literature (Tassi et al., 2005) of Rincón de la Vieja thermal manifestations are also reported. The δ13C-CO2 values of gases from the N flank of Tenorio (Figure 1) were significantly more negative (from −5.44 to −2.80‰ vs. V-PDB) with respect to those from the other thermal discharges (from −2.71 to −1.16‰ vs. V-PDB). The Rc/Ra values of gases from Rincón de la Vieja volcano ranged from 3.55 to 4.66, whereas that measured in HN1 gas (Miravalles volcano) was slightly lower (3.24). On the contrary, significantly higher Rc/Ra values were measured at Tenorio, especially in the AL1 gases (7.11 and 7.07, respectively).
Gas Sources and Secondary Chemical-Physical Processes
The isotopic ratios of He in hydrothermal gas discharges from different geodynamic settings are helpful to evaluate the contribution of mantle- and crust-related fluids (e.g., Polyak and Tolstikhin, 1985). The R/Ra values of gases released along convergent plate boundaries may vary between 3 and 8, depending on the degree of crustal contamination affecting the mantle source (Poreda and Craig, 1989; Hilton et al., 2002). The Rc/Ra values measured in the GGP gases (Table 1) fall within such range (Table 1). The highest mantle He contribution observed in the study area was found in GB (Miravalles), as already highlighted by previous data (Poreda and Craig, 1989), and, surprisingly, in AL (Tenorio) gases. It is worth noting that gas discharges from the same area, such as Borinquen and Las Pailas (Figure 1), i.e., located close to each other, were characterized by relatively strong differences in Rc/Ra values, suggesting that the hydrothermal fluids uprising from the deep reservoirs were affected by variable inputs of crustal gases at a local scale. The δ13C-CO2 values showed overlapping ranges between mantle and crustal CO2 (Sherwood Lollar et al., 1997). Thus, this geochemical parameter cannot clearly distinguish gases from these two sources. The δ13C-CO2 values of most GGP gases, being higher than those that characterize mantle gases (from −7 to −3‰ vs. V-PDB) (Pineau and Javoy, 1983), may be interpreted as due to interactions of limestone with the hydrothermal/magmatic fluids, as also observed in other volcanic systems such as Vulcano Island (Tedesco and Nagao, 1996), Mt. Vesuvius (Chiodini et al., 2001), and Nysiros (Brombach et al., 2003). The δ13C-CO2 values of the gases from the NE sector of Tenorio volcano (Figure 1), which were slightly more negative, may imply (i) a pure mantle CO2 source or (ii) a contribution of 12C-enriched CO2 (≤ − 20‰ vs. V-PDB) from degradation of organic matter (O'Leary, 1988; Hoefs, 2009). These values may also be related to 13C-12C fractionation processes due to interactions with shallow aquifers and/or carbonate precipitation/dissolution (Venturi et al., 2017, and references therein). Useful insights into the origin of CO2 of the GGP gases are provided by combining δ13C-CO2 values and CO2/3He ratios (Marty and Jambon, 1987; O'Nions and Oxburgh, 1988). As shown in Figure 2, all the gases from the three volcanic systems of GGP displayed a strong CO2-excess (1–2 orders of magnitude) with respect to mantle, mostly due to CO2 production from limestone and/or marine carbonate. As suggested by Sano and Marty (1995), the relative contribution of carbon from mantle degassing (M), limestone (L), and/or organic-rich sediments (S) can be computed, as follows:
Equations (3, 4) were solved considering the following data: (δ13C-CO2)MORB = −6.5‰, (δ13C-CO2)Sed = −30‰, (δ13C-CO2)Lim = 0‰, (CO2/3He)MORB = 1.5 × 109, (CO2/3He)Sed = 1 × 1013 and (CO2/3He)Lim = 1 × 1013. According to these calculations (Table 1), the highest mantle CO2 fraction (3.3%) was measured at Altomassis (AL, Tenorio), whereas those of the other emissions ranged from 0.4 to 1.8%. As already shown in Figure 2, limestone was by far the most abundant crustal carbon source (from 85 to 87%), with variable contribution from organic-rich sediments (from 1.8 to 14%). As already reported by Shaw et al. (2003) on the basis of the CO2-He systematics of various thermal fluid discharges from volcanoes and geothermal systems in the Costa Rica and Nicaragua regions including Rincón de la Vieja and Miravalles volcanoes, the fluids from the study area were characterized by a relatively high L/S average ratio (21.7), i.e., strongly higher than the worldwide arc L/S average value (6.0 ± 0.38; Sano and Williams, 1996), possibly due to the carbonate-rich composition of the Cocos Plate subducting sedimentary material (Patino et al., 2000; Saginor et al., 2013). Similar features were also observed in fluids from the Chilean Southern and Central Volcanic Zones (CVZ and SVZ) (Ray et al., 2009; Benavente et al., 2013, 2016).
Figure 2. CO2/3He vs. δ13C-CO2 binary diagram. End-member compositions of organic matter, limestone and mantle are from Sano and Marty (1995). Mixing curves between the three end-members are reported. Symbols as in Figure 1.
The N2/Ar ratios of most GGP gases were comprised between those of Air Saturated Water (ASW) and air (Figure 3), likely due to atmospheric gas contribution from shallow aquifers in contact with the uprising hydrothermal fluids. In contrast, gases from (i) the Miravalles geothermal gases (GW), (ii) the active crater (AC) of Rincón de la Vieja, and (iii) AL were marked by a significant N2-excess (Figure 3). The GW gases were expected to show a limited interaction with the shallow environment, being preserved by the tubing of the wells, whereas the low air contamination of AC was likely related to a strong magmatic fluid inputs, considering that this gas collected from a high-flux emission at the border of the hyperacidic and extremely active carter lake of Rincón de la Vieja (Tassi et al., 2005). The high N2/Ar ratios of the AL gases were surprising since these gas emissions were apparently similar to the other hydrothermal discharges, i.e., low flux bubbling pools with a relatively low (<100°C) outlet temperature (Table 1). Such a N2-excess, as well as the relatively high Rc/Ra values and H2 and CO concentrations (Table 1), clearly distinguish the AL gases from those emitted on the northern flank of Tenorio volcano (Figure 1). A similar chemical composition was also shown by the TM boiling pools (Table 1; Giggenbach and Corrales, 1992), which are located relatively close to the AL discharges (Figure 1). These compositional features suggest that gases from this area were characterized by a relatively high contribution of a deep gas component, a peculiarity that may be related to the occurrence of tectonic structures able to favor the uprising of deep fluids. Unfortunately, the promising indications provided by gas geochemistry cannot be supported by the poorly developed studies on the local structural setting, since, to date, the intense vegetation and the relatively low temperatures of the fluids exploited during the preliminary geothermal prospection (Mayorga, 2009) discouraged appropriate investigations.
Temperature- and redox-sensitive compounds of waters and gases from volcanic systems are commonly used as geoindicators for fluid reservoirs (e.g., Giggenbach, 1991, and references therein). Geothermometers can provide reliable results when (i) the chemical reactions among the species considered for the geothermometric computations have attained equilibrium within the reservoirs, and (ii) no significant compositional changes have occurred during the uprising of fluids toward the surface. According to these assumptions, the composition of fluids exploited from geothermal wells, rapidly uprising through tubing that prevent the contact between the hydrothermal fluids and the external environment, tends to strictly reflect the chemical-physical conditions at the fluid source. On the contrary, secondary processes, such as cooling, steam condensation, mixing with shallow aquifers, and air contamination may affect fluids feeding the natural emissions, modifying the chemical features they have acquired at depth. For example, water vapor concentrations of deep-originated gases discharged from bubbling and boiling pools, such as those of the present study, are strongly depending on condensation processes occurring close to the surface, to shallow depth. Therefore, gas geothermometry was preferably carried out in the CO2-CO-H2-CH4 system. Among these gas compounds, CO and H2 have a similar rapid kinetics in response to changes of the chemical-physical conditions (Giggenbach, 1987, 1996; Chiodini et al., 1993), e.g., during the fluid uprising toward the surface. Thus, the following pressure-independent chemical reaction was considered:
The dependence on temperature of reaction (5) in the vapor phase is given by:
whereas in the liquid phase is, as follows:
where T is in K, XCO and XCO2 are the molar fractions of CO2 and CO, respectively, BCO2 and BCO are the vapor/liquid distribution coefficients of CO2 and CO, respectively, whereas RH is log(XH2/XH2O) (Giggenbach, 1980, 1987). Assuming that Ar concentrations are fixed by the equilibrium between the atmosphere and air saturated water (ASW) (Giggenbach, 1991), the dependence of H2 on RH in the two phases can also be described by the following equations:
where XH2, , , and BH2 are the molar fractions and the vapor/liquid distribution coefficients of Ar* and H2, respectively. The values, computed by Ar* = XAr − (XO2/22) considering that hydrothermal fluids are O2-free, was used instead of the XAr to minimize the effects on hydrothermal fluids of air contamination at the surface. Assuming that the RH values were controlled by the rock redox buffer system proposed by D'Amore and Panichi (1980) (DP), the combination of Equations (6–9) allowed to construct liquid-vapor equilibrium grid reported in Figure 4. Gases from Rincón de la Vieja and Miravalles plot within and out of the liquid-vapor grid, and their apparent equilibrium temperatures ranges from ~180° to >300°C. Gases from the northern sector of Tenorio seem to indicate lower equilibrium temperatures (130–180°C). The strong scattering of these data suggests that secondary processes likely controlled H2 and CO. Interactions with shallow aquifers, which are typically rich in free O2, may have dramatically changed the redox conditions acting on the uprising fluid, severely affecting the behavior of these reduced gas species. Moreover, CO tends to dissolve in liquid water to form HCOOH (Shock, 1993).
Figure 4. Log(XH2/XAr*) vs. log(XCO/XCO2) binary diagram. Symbols as in Figure 1.
On the contrary, the AL gases, which cluster close to the vapor equilibrium curve at 240–250°C (Figure 4), seem to have maintained the chemical composition attained at depth. This is consistent with the relatively high fraction of mantle He and CO2 (Table 1) and the high N2/Ar ratios (Figure 3) of these gases. The H2/Ar values of the AL gases were similar to those measured by Giggenbach and Corrales (1992) in the TM gas emission (Table 1), located in the same area (Figure 1). However, those authors surprising interpreted these data as evidence of low-temperature conditions.
To complete the geothermometric estimations in the CO2-CO-H2-CH4 system, we hypothesized that the CH4 concentrations were controlled by the Sabatier reaction, as follows:
Considering that log(f H2O) = 4.9 − 1820/T (Giggenbach, 1987), the dependence on temperature and RH of the XCH4/XCO2 log-ratios in the vapor and liquid phases can be expressed, as follows:
where BCH4 is the vapor/liquid distribution coefficient of CH4.
The log(XH2/XAr*) vs. log(XCH4/XCO2) diagram (Figure 5) was constructed on the basis of Equations (8–12) and assuming that the RH values were regulated the DP redox buffer system. It confirms a strong chemical disequilibrium of the gases from the three volcanoes, which is consistent with the large variations of the Rc/Ra values measured in gas emissions from the same volcanic system (i.e., Rincón de la Vieja, Table 1). AL gases are the only exception, their equilibrium temperatures were the same as those indicated by the H2-CO pair. It is worth noting that the AL calculated temperatures were consistent with those hypothesized for the geothermal reservoir associated with the Miravalles volcano (Gardner and Corrales, 1977). It is worth noting that the CH4-CO2 ratios of the AL gases were one order of magnitude lower than those of the TM and MZ gases (Giggenbach and Corrales, 1992) (Figure 1), suggesting a large chemical heterogeneity even at a relatively small spatial scale. The XH2/XAr* and XCH4/XCO2 ratios of the Miravalles geothermal gases show strong variations (~1 log-unit; Figure 5), although supposedly less affected by the influence of the shallow environment. Comparable results were also found by Gherardi et al. (2002), who suggested that the fluids from Miravalles suffered the influence of fluid reinjection.
Figure 5. Log(XH2/XAr*) vs. log(XCH4/XCO2) binary diagram. Symbols as in Figure 1.
The occurrence of a strong heat source, demonstrated by the numerous thermal manifestations of GGP, and the abundant meteoric water recharge characterizing the volcanic chain of Guanacaste explain the high geothermal potential of this region. The chemical and isotopic features of the discharged fluids are dictated by the interaction of the uprising hydrothermal fluids and shallow aquifers. At a first approximation, the highest deep-originated fluid contribution can be recognized at the active crater of the Rincón de la Vieja volcano, where phreatic and phreato-magmatic eruptions frequently occur, and progressively decreases with increasing distance from this site. Carbon dioxide production from limestone seems the main source for this gas, possibly related to carbonate assimilation at the magmatic source and/or due to interaction of carbonate-rich rocks and the hydrothermal fluids. The geochemical features of the fluids exploited from Miravalles correspond to those typically found in the surroundings of active volcanic systems, i.e., Na-Cl brines produced by prolonged fluid-rock interactions and dissolution of magmatic gases into a well-developed hydrothermal aquifer. The concentrations of minor reduced gases, such as CO, H2, and CH4 were likely controlled by chemical reactions occurring within the hydrothermal deep fluid source, although the chemical composition attained at depth was strongly modified by secondary processes resulting in a large scattering of the chemistry of the surficial manifestations. According to the scenario suggested by Giggenbach and Corrales (1992), the Tenorio thermal discharges would represent the distal fluid outflow of the regional geothermal reservoir. The new geochemical evidences reported in the present study significantly modify this schematic model. As shown in Figure 6, the relatively high contribution of CO2 from the mantle (M in %), as well as the high Rc/Ra and N2/Ar ratios shown by the fluid discharges located in the southern sector of Tenorio volcano with respect to those from Miravalles and Rincón de la Vieja, suggest a high-temperature source in this area. This implies that, although the preliminary geochemical prospections carried out in the 1990's produced discouraging results, an important geothermal resource is hosted in this peripheral zone of GGP, where local tectonics favor the uprising of hot fluids minimizing their cooling and mixing with crustal and atmospheric gases. In light of this hypothesis, which needs to be supported by geophysical prospection activities and accurate investigations of the geologic and tectonic settings, promising results are expected to be obtained by extending the exploitation activity to the south of the Miravalles geothermal plant. Accordingly, the geothermal potential of GGP should likely be regarded as higher than previously thought.
Figure 6. Schematic draft of the spatial distribution of Rc/Ra and N2/Ar ratios and mantle CO2 fraction values (M in %) in hydrothermal fluid discharges from the studied volcanic systems of Guanacaste. Symbols as in Figure 1.
FT: field campaigns, gas analysis, manuscript preparation; OV, GB, FC, and JdM: field campaigns, manuscript preparation; GP: gas analysis, manuscript preparation; SV: manuscript preparation.
Conflict of Interest Statement
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.
This work was financially supported by the laboratories of Fluid and Rock Geochemistry and Stable Isotope Geochemistry of the Department of Earth Sciences and the Institute of Geosciences and Earth Resources of the National Research Council of Italy (IGG-CNR) of Florence.
Alvarado, G. E., and Cárdenes, G. (2016). “Geology, tectonics, and geomorphology of costa rica: a natural history approach,” in Costa Rican Ecosystems, ed M. Kappelle (Chicago, IL: The University of Chicago Press), 30–63.
Alvarado, G. E., and Gans, P. B. (2012). Síntesis geocronológica del magmatismo, metamorfismo y metalogenia de Costa Rica, América Central. Rev. Geol. Am. Central 46, 7–122. doi: 10.15517/rgac.v0i46.1836
Alvarado, G. E., Kussmaul, S., Chiesa, S., Gillot, P. Y., Appel, H., Worner, G., et al. (1992). Resumen cronoestratigrafico de las rocas igneas de Costa Rica basado en dataciones radiometricas. J. South Am. Earth Sci. 6, 151–168. doi: 10.1016/0895-9811(92)90005-J
Benavente, O., Tassi, F., Gutiérrez, F., Vaselli, O., Aguilera, F., and Reich, M. (2013). Origin of fumarolic fluids from Tupungatito Volcano (Central Chile): interplay between magmatic, hydrothermal, and shallow meteoric sources. Bull. Volcanol. 75, 1–15. doi: 10.1007/s00445-013-0746-x
Benavente, O., Tassi, F., Reich, M., Aguilera, F., Capecchiacci, F., Gutiérrez, F., et al. (2016). Chemical and isotopic features of cold and thermal fluids discharged in the Southern Volcanic Zone between 32.5°S and 36°S: insights into the physical and chemical processes controlling fluid geochemistry in geothermal systems of Central Chile. Chem. Geol. 420, 97–113. doi: 10.1016/j.chemgeo.2015.11.010
Boudon, G., Rancon, J. P., Kieffer, G., Soto, G. J., Traineau, H., and Rossignol, J. C. (1996). Les eruptions de 1966-1970 et 1991-1992 du volcan Rincón de la Vieja, Costa Rica: exemple d'activité réecurrente d'un systéme hydromagmatique. C. R. Acad. Sci. 322, 101–108.
Brombach, T., Caliro, S., Chiodini, G., Fiebig, J., Hunziker, J., and Raco, B. (2003). Geochemical evidence for mixing of magmatic fluids with seawater, Nisyros hydrothermal system, Greece. Bull. Volcanol. 65, 505–516. doi: 10.1007/s00445-003-0278-x
Carr, M. J. (1984). Symmetrical and segmented variation of physical and geochemical characteristics of the Central America Volcanic front. J. Volcanol. Geotherm. Res. 20, 231–252. doi: 10.1016/0377-0273(84)90041-6
Chavarría, L., Mora, O., Hakanson, E., Galvez, M., Rojas, M., Molina, F., et al. (2010). “Geologic model of the pailas geothermal field, guanacaste, costa rica,” in Proceedings World Geothermal Congress (Bali).
Chiesa, S., Alvarado, G. E., Pecchio, M., Corrella, M., and Zanchi, A. (1994). Contribution to petrological and stratigraphical understanding of the Cordillera de Guanacaste lava flows, Costa Rica. Rev. Geol. Am. Central 17, 19–43.
Chiesa, S., Civelli, G., Gillot, P. Y., Morz, O., and Alvarado, G. E. (1992). Rocas piroclásticas asociadas a la formación de la Caldera de Guayabo, Cordillera de Guanacaste, Costa Rica. Rev. Geol. Am. Central 14, 59–75.
Chiodini, G., Cioni, R., and Marini, L. (1993). Reactions governing the chemistry of crater fumaroles from Vulcano Island, Italy, and implications for volcanic surveillance. Appl. Geochem. 8, 357–371. doi: 10.1016/0883-2927(93)90004-Z
Chiodini, G., Marini, L., and Russo, M. (2001). Geochemical evidence for the existence of high-temperature hydrothermal brines at Vesuvio volcano, Italy. Geochim. Cosmochim. Acta 65, 2129–2147. doi: 10.1016/S0016-7037(01)00583-X
Gherardi, F., Panichi, C., Yock, A., and Gerardo-Abaya, J. (2002). Geochemistry of the surface and deep fluids of the Miravalles volcano geothermal system (Costa Rica). Geothermics 31, 91–128. doi: 10.1016/S0375-6505(01)00030-X
Giggenbach, W. F., and Corrales, R. S. (1992). Isotopic and chemical composition of water and steam discharges from volcanic-magmatic-hydrothermal systems of the Guanacaste Geothermal Province, Costa Rica. Appl. Geochem. 7, 309–332. doi: 10.1016/0883-2927(92)90022-U
Instituto Costarricense de Electricidad (ICE) (1990). Estudios de Reconocimiento y Prefactibilidad Geotérmica en la República de Costa Rica. Fase II, Estudio de Prefactibilidad del área del Tenorio, Informe Geovulcanológico. Proyecto COS/83/T01 [Inf. Interno].
Instituto Costarricense de Electricidad (ICE) (2014). Plan de Expansión de la Generación Eléctrica, Período 2014-2035. Report, Centro Nacional de Planificación Eléctrica, Proceso de Expansión Integrada; Instituto Costarricense de Electricidad; San José, Costa Rica.
Kempter, K. A., Benner, S. G., and Williams, S. N. (1996). Rincón de la Vieja volcano, Guanacaste province, Costa Rica: geology of the south-western flank and hazard implication. J. Volcanol. Geotherm. Res. 71, 109–127. doi: 10.1016/0377-0273(95)00072-0
Mayorga (2009). Mayorga. Available online at: http://www.kinertech.com/yahoo_site_admin/assets/docs/Geotermia__Caso_de_Costa_Rica.27362721.pdf
Patino, L. C., Carr, M. J., and Feigenson, M. D. (2000). Local and regional variations in Central American arc lavas controlled by variations in subducted sediment input. Contrib. Miner. Petrol. 138, 265–283. doi: 10.1007/s004100050562
Ray, M., Hilton, D., Muñoz, J., Fischer, T., and Shaw, A. (2009). The effects of volatile recycling, degassing and crustal contamination on the helium and carbon geochemistry of hydrothermal fluids from the southern volcanic zone of Chile. Chem. Geol. 266, 38–49. doi: 10.1016/j.chemgeo.2008.12.026
Rizzo, A., Barberi, F., Carapezza, M. L., Di Piazza, A., Francalanci, L., Sortino, F., et al. (2015). New mafic magma refilling a quiescent volcano: evidence from He-Ne-Ar isotopes during the 2011-2012 unrest at Santorini, Greece. Geochem. Geophys. Geosyst. 16, 798–814. doi: 10.1002/2014GC005653
Saginor, I., Gazel, E., Condie, C., and Carr, M. J. (2013). Evolution of geochemical variations along the Central American volcanic front. Geochem. Geophys. Geosyst. 14, 4504–4522. doi: 10.1002/ggge.20259
Shaw, A. M., Hilton, D. R., Fischer, T. P., Walker, J. A., and Alvarado, G. E. (2003). Contrasting He-C relationships in Nicaragua and Costa Rica: insights into C cycling through subduction zones. Earth Planet. Sci. Lett. 214, 499–513. doi: 10.1016/S0012-821X(03)00401-1
Sherwood Lollar, B., Ballentine, C. J., and O'Nions, R. K. (1997). The fate of mantle-derived carbon in a continental sedimentary basin: integration of C/He relationships and stable isotope signatures. Geochim. Cosmochim. Acta 61, 2295–2307.
Tassi, F., Vaselli, O., Capaccioni, B., Giolito, C., Duarte, E., Fernandez, et al. (2005). The hydrothermal-volcanic system of Rincón de la Vieja volcano (Costa Rica): a combined (inorganic and organic) geochemical approach to understanding the origin of the fluid discharges and its possible application to volcanic surveillance. J. Volcanol. Geotherm. Res. 148, 315–333. doi: 10.1016/j.jvolgeores.2005.05.001
Tedesco, D., and Nagao, K. (1996). Radiogenic 4He, 21Ne and 40Ar in fumarolic gases at Vulcano island: implication for the subductuted African continental plate beneath Italy. Earth Planet. Sci. Lett. 144, 517–528. doi: 10.1016/S0012-821X(96)00196-3
Venturi, S., Tassi, F., Bicocchi, G., Cabassi, J., Capecchiacci, F., Capasso, G., et al. (2017). Fractionation processes affecting the stable carbon isotope signature of thermal waters from hydrothermal/volcanic systems: the examples of Campi Flegrei and Vulcano Island (southern Italy). J. Volcanol. Geotherm. Res. 345, 46–57. doi: 10.1016/j.jvolgeores.2017.08.001
Keywords: geothermal prospection, gas geothermometry, fluid geochemistry, fluid reservoir, Guanacaste province
Citation: Tassi F, Vaselli O, Bini G, Capecchiacci F, de Moor JM, Pecoraino G and Venturi S (2018) The Geothermal Resource in the Guanacaste Region (Costa Rica): New Hints From the Geochemistry of Naturally Discharging Fluids. Front. Earth Sci. 6:69. doi: 10.3389/feart.2018.00069
Received: 21 March 2018; Accepted: 17 May 2018;
Published: 05 June 2018.
Edited by:Jacob B. Lowenstern, Cascades Volcano Observatory (CVO), Volcano Disaster Assistance Program (USGS), United States
Reviewed by:Loic Peiffer, Centro de Investigación Científica y de Educación Superior de Ensenada, Mexico
Hiroshi Shinohara, Geological Survey of Japan (AIST), Japan
Copyright © 2018 Tassi, Vaselli, Bini, Capecchiacci, de Moor, Pecoraino and Venturi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Franco Tassi, firstname.lastname@example.org