Next Article in Journal
Editorial of the Special Issue “Geomorphology, Geoheritage, Geoparks and Geotourism in Volcanic Areas”
Previous Article in Journal
Defining the Boundary Conditions for Seismic Response Analysis—A Practical Review of Some Widely-Used Codes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contributions for the Understanding of the São Pedro do Sul (North of Portugal) Geohydraulic and Thermomineral System: Hydrochemistry and Stable Isotopes Studies

1
Departamento de Engenharia Civil e Arquitetura, Faculdade de Engenharia, Universidade da Beira Interior, R. Marquês de Ávila e Bolama, 6201-001 Covilhã, Portugal
2
Geobiociências, Geoengenharias e Geotecnologias (Geobiotec), Universidade de Aveiro, Campus Universitário de Santiago, 3810-193 Aveiro, Portugal
3
Departamento de Geologia, Escola de Ciências da Vida e do Ambiente, Universidade de Trás-os-Montes e Alto Douro, Quinta de Prados, 5001-801 Vila Real, Portugal
4
Centro de Geociências (CGeo), Universidade de Coimbra, Rua Sílvio Lima, Universidade de Coimbra—Polo II, 3030-790 Coimbra, Portugal
5
Centro de Ciências e Tecnologias Nucleares (C2TN), Instituto Superior Técnico, Universidade de Lisboa, Estrada Nacional 10, ao km 139.7, 2695-066 Bobadela LRS, Portugal
*
Author to whom correspondence should be addressed.
Geosciences 2022, 12(2), 84; https://doi.org/10.3390/geosciences12020084
Submission received: 21 December 2021 / Revised: 8 February 2022 / Accepted: 10 February 2022 / Published: 13 February 2022

Abstract

:
São Pedro do Sul thermomineral aquifer system is located in the North of Portugal. Hydrogeochemistry and isotopic studies were conducted in order to improve knowledge of this groundwater system, known since ancient Roman times for their therapeutic properties. One thermomineral spring (NT) and three boreholes (AC1, SDV1 and SDV2) have a mean temperature around 68 °C. Currently, these waters are mainly used in thermal spas and for geothermal energy supply. Major cations and anions were used to identify and characterize different water types and sub-systems. Stable isotopes composition, δ18O and δ2H, have been used to determine the origin and have been used as a first approach to estimate the mean preferential recharge altitude of thermomineral water. The results suggest a single aquifer system with a relative composition of cations and anions and similar mineralization values. δ18O and δ2H values indicate a meteoric origin with no significant evaporation before infiltration. Besides, the isotopic composition points to recharge areas located at altitudes of about 1000 m a.s.l. This suggests a recharge area in the Freita/Arada mountains located NW of the thermomineral poles.

1. Introduction

Since the Roman occupation, spas have been a tradition in Portugal, owing mainly to the great diversity in the chemistry of its natural mineral waters, which reflect the geological and structural diversity of the territory. The research region is located near the village of São Pedro do Sul, district of Viseu, central Portugal (Figure 1), where this thermomineral aquifer system represents one of the most important sites of thermomineral waters of Portugal. The economic importance of spas for regional development can be inferred by the number of visitors every year to these facilities. According to [1], in 1998, the number of spa users amounted to 20,012. Due to the above reasons (health and economics), protection areas are of vital importance to guarantee the quality of the mineral water and to prevent contamination problems.
Two mineral water producing poles can be recognized within the area of the study at a distance of around 1.2 km from one another. Termas Pole consists of a bathhouse for medicinal treatments and a geothermal central for the heating of non-mineral waters for domestic use (heating and water supply). This area is supplied by the borehole AC1, with a 500 m depth, and by an ancient spring (Traditional Spring—NT, used for medical treatments since the Roman empire) which together produce, through artesianism, about 18 L/s at 67.6 °C [1]. According to these authors, a small mineral spring and two wells, the SDV1 with a depth of 216 m and SDV2 with a depth of 151 m, characterize the Vau pole. The water flow of 1.5 L/s at a temperature of 67 °C, mainly produced from well SDV1, was used to heat greenhouses for banana and pineapple production until a few years ago.
The main goals of this investigation were to investigate the characterization of São Pedro do Sul mineral waters though an isotopic and hydrogeochemical approach. A hydrogeochemistry and isotopic study was conducted in order to improve knowledge about the São Pedro do Sul groundwater system: Major cations and anions were used to define different water types and to understand the water–rock interactions of thermal waters; the reservoir temperature was estimated using geothermometers from measured SiO2, K+, and Mg2+ chemical species. In addition, stable isotopes 18O and 2H were used as indicators of the origin of the thermal water and as a first approach to understanding the recharge areas of the thermomineral water. Numerous studies related to chemistry were performed by many scientists [2,3,4,5,6,7,8,9,10,11,12]. Besides, aspects of geology, hydrogeology, hydrochemistry of geothermal water in the study region were performed during the last decades through conference proceedings, scientific articles, theses and internal reports [3,5,13,14,15,16,17,18,19,20].
Among the frequently used techniques in hydrogeology to characterize aquifer dynamics, environmental isotopes such as 18O and 2H have been the most useful in terms of providing new insights to understand the changing aspects within aquifer systems, and they are sometimes used as a background tool in the development of conceptual models (e.g., preferential recharge altitudes). In the last decades, isotopic signatures of meteoric waters have become a key tracer, intensively used in hydrogeological studies as a basis for the elaboration of conceptual hydrodynamic models [4,6,21,22,23,24,25,26,27,28,29,30,31].
In most cases, groundwater recharge results from direct infiltration of precipitation (regional precipitation). In this case, groundwater samples will reflect the mean isotopic composition of the regional precipitation. If changes occur between the groundwater and precipitation, an explanation for this deviation should be sought, and different hypotheses can be formulated, for example, recharge derived from a different source, either from regional precipitation, for example from surface water (rivers, lakes or dams) instead of precipitation; in this case, the groundwater will reflect the mean isotopic composition of the contributing source (river, lake or dam waters), or if evaporation occurs before infiltration, an enrichment in the isotopic composition will be noted.
For this purpose, water samples were collected from different water sources, with respect to different types of rocks and altitudes and were analyzed for important major parameters and isotopic data following standard test procedures. For this purpose, a total of 16 water samples (Figure 2) were collected from a cool spring in June and July 2016 in the surrounding area of São Pedro do Sul site, and a set of major cations and anions results obtained over nine years in NT, AC1, SDV1 and SDV2 was employed in the Piper trilinear diagram plot (Section 4.1). Conversely, 19 samples were used for the isotopic analysis with NT, AC1, SDV1 and SDV2 included (Figure 2). In order to achieve this objective, groundwater samples were collected and analyzed.

2. Geological and Hydrogeological Settings

The geological and structural regional map (Figure 2) was built by a combination of different data, which can be summarized in the following resources: geological map of Portugal, sheets 13-D—Oliveira de Azeméis, 14-C—Castro Daire and 17-A—Viseu, scale 1:50,000, as well as the written documents—Notícia Explicativa— which provides detailed and additional information about each geological map [32,33,34,35,36,37]; geological map of Portugal at 1:500,000 scale [38]; aerial and satellite images available on Google Earth; digital elevation model (DEM), produced from 1:25,000 Portuguese military topographic maps (sheets number 165, 166, 176 and 177) [39]; studies already performed in the area such as internal reports, geological and structural documents [14,17,40] and field geo-structural survey carried out for this work.
The São Pedro do Sul region is characterized by hercynic granite, emplaced between 290 and 321 Myr, in relation to the third phase of Variscan deformation (F3) occupying extensive areas [19,32,33,35,36,37,38]. These hercynian granites can be separated into two outcrops: (1) Syn-tectonics granites in relation to F3 Variscan deformation [40] peraluminous [41] with two micas. These batholiths occur throughout the Serra da Freita/Arada-São Pedro do Sul corridor, occupying a large area around São Pedro do Sul Thermal Poles (Figure 2). These granites have an anisotropic structure, with a NW–SE direction related to the deformation of the hercynian tectonic axes. (2) Post-tectonics granites in relation to F3 Variscan deformation that are porphyroid granites composed mainly of biotite, quartz, microcline, plagioclase (andesine and albite) and few muscovite and occur mainly in the NE of São Pedro do Sul.
Metasedimentary rocks are present in the region, although with less importance than granitic rocks. These metasedimentary rocks consist of grey metaquartzites, with calco-silicated rocks layers, schist, phylites, mica schist and metagraywackes of the Cambrian age (and Upper Precambrian). These rocks are intruded by the hercynian granites and present a schist structure with a NW–SE direction.
Amphibolitized basic rock veins occur in metasedimentary rock levels according to the NW–SE orientation. Furthermore, quartz veins are found in association with both granitic and sedimentary rocks related to the NNE–SSW to NE–SW and NW–SE fault directions.
In addition, quaternary sedimentary deposits can be found in a few areas: fluvial terraces, mainly with gravel, are located between 100 and 200 m a.s.l.; alluvial, mainly with sand and gravel, are found in the Holocene floodplains.
The fault pattern of the regional complex is genetically correlated to the Hercynian orogeny and has been reactivated by neotectonic activity [20]. The São Pedro do Sul region appears to be highly faulted as indicated in Figure 2. This regional geological map shows three dominant fault families with orientations: (i) N–S to NNE–SSW; (ii) WNW–ESE to NW–SE; and (iii) NE–SW. The current regional tectonic stress field with maximum horizontal compressive stress trends from NW–SE to WNW–ESE [20] can promote depth fluid circulation from WNW–ESE to NW–SE faults.
Representing the region by a local scale (Figure 3), three main granite facies, which belong to the syn-tectonics granites, were identified in the São Pedro do Sul Thermal Poles area; namely [34,40], Vouzela granite, São Pedro do Sul granite and Fataunços granite. Vouzela granite is a monzonitic granite medium to coarse grained, formed mainly by quartz, albite-oligoclase, potassium feldspar, biotite and muscovite; São Pedro do Sul is a fine to medium grained granite composed mainly of microcline, plagioclase, quartz, muscovite and biotite; Fataunços is a granite close to granodiorite, fine- to medium-grained, with two micas. All these granites are slightly porphyroid textured.
From the structural point of view, the local geologic faults system is mainly represented by NNE–SSW, NE–SW and WNW–ESE structural directions, where, respectively, the São Pedro do Sul-Ribamá Fault, Termas Fault and Fataunços Fault are examples (Figure 3a). These hercynian fault systems are active today and represent tectonic structures with high regional importance. The regional or local fracture networks support the main fractured aquifer systems. Studies developed in the region [7,20] show that there are two types of aquifer systems associated with fractured granites and metasedimentary rocks. The shallow cold aquifer systems occur throughout the region with lightly mineralized water, while a deeper aquifer system is characterized by higher temperature and strongly mineralized water, to which the thermal poles are closely linked. In addition, a close relationship between this aquifer and São Pedro do Sul-Ribamá Fault was identified. Differences in hydrochemical characteristics are included in Section 4 (Results and Discussion).
Pereira and Ferreira [40] studied the importance of the structural control of the thermomineral springs at São Pedro do Sul, mentioning the positive relation between the active faults and the occurrence of thermal waters and the importance that these structures have in its mineralization and temperature. Furthermore, these authors proposed a schematic cross section important in the conceptual circulation model of São Pedro do Sul thermal waters (Figure 3b). Thermal poles and their hot springs are linked with local fault orientation and connectivity, namely the Termas Fault oriented NE–SW, the Fataunços Fault oriented WNE–ESE, and the São Pedro do Sul-Ribamá Fault oriented NNE–SSW.

3. Methodology

Groundwater samples were collected from boreholes and springs in order to obtain a representative composition of the thermomineral and cold shallow aquifer systems of the São Pedro do Sul region (Figure 2). The selected cold springs intend to represent the shallow groundwater of the region, characterized by the low mineralization, low temperature, and pH values around 5 to 6 (see Table 1). These selected springs are located at different altitude sites within the surroundings of the research region and represent the initial recharge water composition, “before” water-rock interaction processes, as found in the deep aquifer.
The sampling procedures, such as collection, handling, storage and transport, were carefully planned and applied to ensure the safety and quality of sampling. All the processes took place in line with the relevant specifications or standard practices. Thus, the bottles were cleaned with distilled water before sampling and then pre-rinsed with water from the sampling site two or three times before collecting. The water samples were collected in wide-neck high-density polyethylene (HDPE) bottles (filled without leaving any air bubbles), stored in the absence of light and refrigerated at a constant low temperature of 4 °C until analysis in the laboratory. Water samples for cations were acidified with concentrated nitric acid to keep a pH < 2 and were reserved in 500 mL bottles. For anion analysis, 500 mL bottles were used, while samples for oxigen-18 and deuterium were stored in 50 mL bottles. All water samples were marked, sealed, and stored, and then sent to the laboratory for analysis as soon as possible after collection. In addition, field-measured parameters of pH, electrical conductivity (EC) and temperature were performed in situ using a portable multiparameter WTW Multi 3420 calibrated with a standard solution.
The chemical analyses were performed at LNEG (Laboratório Nacional de Energia e Geologia), Portugal, for major dissolved species determination using the following methods: atomic absorption spectrometry for Ca2+ and Mg2+; emission spectrometry for Na+ and K+; colorimetric methods for SiO2, F; ion chromatography for SO42− and Cl; potentiometry for alkalinity, here referred as HCO3 (Table 1). The stable isotope composition of all groundwater samples was determined at C2TN/IST (Centro de Ciências e Tecnologias Nucleares, Instituto Superior Técnico, Universidade de Lisboa) Portugal. The δ2H and δ18O content is reported in per mil (‰) vs. V-SMOW and was obtained by Laser Spectroscopic (LGR 24d).
As mentioned before, a hydrogeochemistry and isotopic study was conducted in order to improve knowledge about the São Pedro do Sul thermomineral groundwater system. From chemical data, the cations and anions content was used to identify the different water facies and possible evolution trends within the groundwater flow, namely using: (i) the Piper diagram to identify different water types; (ii) Schoeller–Berkaloff diagram to highlight the hydrogeochemistry profile of water samples and to evaluate the main similarities and differences between them.
The groundwater recharge in both aquifer systems (shallow and deep) derives from direct infiltration of precipitation; groundwater samples in both units will reflect the mean isotopic composition of that regional precipitation. Using the “relation”, besides for stable isotopic composition, δ18O and δ2H content were used to determine the origin and a first approach to estimate the mean preferential recharge altitude areas of the thermomineral water.
A total of 12 groundwater samples were collected inside of the investigated region, among these, 7 cold springs (AS1, AS2, AS3, AS4, AS5, AS6, AS7) and 5 wells (AF2, AF3, AF4, AF5, AF1A) with a mean depth of around 70 m were sampled between June and July 2016 in the surrounding areas of the São Pedro do Sul site.

4. Results and Discussion

4.1. Hydrogeochemical Approach

At São Pedro do Sul, the geological and tectonic environment seems to control some of the major hydrogeological features and processes of the thermomineral waters, such as infiltration, aquifer recharge, groundwater flow paths, and hydrogeochemical evolution. In Table 1, the chemical data presented related to the thermomineral water composition represent the mean composition. The set of major cations and anions results was obtained over nine years of survey in NT, AC1, SDV1 and SDV2. This option (mean values) was chosen considering the high stability in the composition of this group of sites along the years. Besides, from Table 1, it is also possible to note the minor dissimilarity of content in the different parameters within the four sampling sites, presenting similar physic-chemical composition.
Considering the data of Table 1, the chemical characteristics of the thermal water of São Pedro do Sul can be classified as sulfureous waters, with high chemical composition stability. The major elements in these waters are hydrogen carbonate, sodium, silica, chloride and fluoride. Hydrolysis of Na-plagioclase should be considered as the main process responsible for the occurrence of the HCO3–Na waters because Na-plagioclase is the main mineral of the dominant granitic rock of the regional hydrogeological system.
The water samples collected in shallow cold aquifer systems are characterized by the low total mineralization content. These waters should be considered as a signature of local recharge associated with a shallow and short flow path regimen, explaining by this assumption, the development of the Na–Cl waters.
The thermomineral system is characterized by a mean output temperature of around 67 to 68 °C (spring and boreholes) and high pH values around 8.8. Besides the total mineralization values in the range of 365–370 mg/L and the waters showing a HCO3–Na facies; the high silica and fluoride contents, with values around 70 mg/L and up to 18 mg/L, respectively, also represent distinct parameters if they are compared with the shallow cold dilute waters. An important contribution of reduced species of sulfur ≈ 10 mg/L was also determined in this group of waters.
On the contrary, the shallow cold waters present low mineralization, low pH (around 6), and a HCO3-Cl-Na facies with lower silica, fluor and sulfur content. Nevertheless, focusing only on the cold groundwaters systems, a difference in the chemical and physical parameters is noted when the springs are compared with the boreholes. With the exception of AF4 (see Table 1), all springs present lower mineralization, ECmean = 34.5 ± 16.6 µS/cm, while for the boreholes, ECmean = 87.3 ± 21.6 µS/cm. This difference probably indicates a not-so-deep circulation of the springs, with shorter residence time allowing for minor water–rock interaction processes.
The chemical signatures found in the deep water system suggest that the thermomineral fluid’s composition is mainly dominated by the hydrolysis of Na plagioclase within the granitic reservoir. Groundwater at this depth shows an increase in salinity and in pH, with an increase in sodium over calcium compared with the shallow groundwaters. Moreover, considering the dominant lithology of the research region (Hercynian granitoids), oxidation of sulfide minerals (such as pyrite) at this depth should explain the sulphate content increase in the thermomineral groundwater system.
Granitic rocks dominate the research area, where the feldspars hydrolysis represents an important water–rock interaction process, contributing to a high percentage of dissolved silica. The silica content represents about 20% of the total mineralization of the thermomineral water samples (see Table 1), and higher than 30% in the shallow cold aquifer systems.
The direct application of chemistry predicted that the water–rock interaction processes involved in the geochemical evolution of the groundwaters and these parameters were applied in the estimation of the reservoir temperature and the depth of the thermomineral waters [42,43]. To estimate this, the chemical geothermometers SiO2 and K2/Mg proposed by Truesdell (1975) and Giggenbach (1988), respectively, (Equations (1) and (2)) [44,45] were selected based on the acceptable and consistent estimation of reservoir temperature in low enthalpy systems:
T (°C) = [1315/5.205 − log SiO2] − 273.15
and
T (°C) = (4410/14 − log (K2/Mg)) − 273.15
In both equations, the concentrations of K, Mg and SiO2 is in mg/L. The results obtained at São Pedro do Sul thermomineral waters point to a mean reservoir temperature between 119 °C (SiO2) and 116 °C (K2/Mg), indicating a maximum depth from 3.10 to 3.19 km. The results obtained are very consistent (Table 2). To calculate the maximum depth “reached” by the thermomineral waters”, the following equation was considered: depth = (Tr − Ta)/gg, where Tr stands for the reservoir temperature (°C), Ta for the mean regional annual temperature (13 °C), and gg the regional geothermal gradient of 30 °C/km [46].
Using the available data of the deep and shallow aquifer systems, the data were plotted in a Piper diagram (Figure 4) and in a Schoeller–Berkaloff diagram (Figure 5).
The Piper diagram [47] shows three different domains: cluster (a), standing for the thermomineral waters that are sodium bicarbonate; the second and third groups of water are composed of the water samples from the shallow cold aquifer systems. The non-mineral waters are represented in this diagram by two clusters: the first is cluster (b), composed of sodium chlorinated waters, mainly associated with water–rock interaction processes with granitic rocks, while cluster (c) is composed of sodium, calcium or magnesium bicarbonate waters. The increase in Ca2+ and/or in Mg2+ is associated with interaction processes with granitic and metasedimentar rocks or cover deposits where there is agricultural land use. The borehole AF2 reveals a chloride value significantly distinct from the other samples, probably resulting from local anthropogenic activities.
From Figure 5 (Schoeller–Berkaloff diagram), the first impression is that all thermomineral and cold waters samples seem to present the same evolution pattern, although with different contents. However, the deviation in Mg content observed in NT, AC1, SDV1 and SDV2 thermomineral waters can most probably be ascribed to lithology heterogeneities, such as absence or presence along the groundwater flow of amphibolite veins, in which its existence can increase Mg content.
The hydrogeochemical signatures identified in the two water systems (thermomineral and shallow cold aquifer systems) are ascribed to the regional/local lithological differences, together with the anthropogenic action (pollution) in the region. The different chemical type of waters is well observed when the Na + K or Ca + Mg content is plotted as a function of the hydrogen carbonate (Figure 6a,b). Besides, the possible existence of mixing between the thermal fluids with the local shallow cold dilute groundwater systems or even with surface waters (river Vouga that cross the area) was not identified.
The same conclusion of the absence of mixing with the thermomineral waters is also supported as shown by Figure 7; the deeper waters are plotted apart from all the other samples. No transitional water sample compositions were found in the region, corroborating the hypothesis of no mixing between the different aquifers. Conversely, the dispersion within the shallow cold diluted groundwater samples is enhanced when the Ca + Mg content is plotted as a function of HCO3 or as a function of SiO2 content, which is probably related to metasedimentary formations. The same dispersion is noted in the Cl content. (Figure 7a), however this increase is associated to anthropogenic activities in the surrounding area (agricultural activities). Conversely, it is also noted that the majority of the springs exhibits lower mineralization when compared with the boreholes. This suggests shallow circulation pathways and smaller residence time, which is important in water-rock interactions processes.
Furthermore, Figure 6a and Figure 7a,b suggest an evolution of the hydrochemical process from the non-mineral waters to thermomineral waters, dominated by hydrolysis of Na plagioclase present in the regional granitic rocks.
The thermomineral waters together with the non-mineral water samples exhibit identical hydrochemical profiles, but they present an enrichment of the component concentrations, with the exception of magnesium concentration, which decreases in absolute terms (Figure 5). There is a considerable increase in sodium and hydrogen carbonate concentrations in the mineral water composition, which may be associated with hydrolysis of Na plagioclase of granitic rocks of the São Pedro do Sul region [7]. The concentrations of sodium, hydrogen carbonate and SiO2 results from a major residence time of these mineral waters. In relation to non-mineral waters, in general, the hydrochemical profile is similar, except for the samples whose composition is influenced by anthropic activities.

4.2. Stable Isotopes 2H and 18O

At the research area, 19 groundwater samples were collected for stable isotopic determinations (see Table 1), from AS1 to AS15 cold shallow diluted springs and from NT, AC1, SDV1 and SDV2, thermomineral waters. The isotopic values were plotted in the orthogonal diagram (Figure 8); besides, in this diagram, the Global Meteoric Water Line (GMWL) defined by Craig (1961) [48] and the Local Meteoric Water Line are plotted.
The isotopic composition of all groundwater samples (thermomineral and cold diluted groundwater samples) were used to establish the Local Meteoric Water Line (Local-MWL). The obtained equation (δ2H = 6.5δ18O + 5.6; n = 19; r2 = 0.95) was used since no isotopic precipitation record of this region is available and to observe the difference between the Local-MWL and the Global MWL established by Craig (1961). The local δ18O-δ2H relationship for precipitation often varies from the G-MWL in the slope and the intercept, since they are a function of meteorological and geographic characteristics of the region.
From the distribution of the groundwater samples, a common origin of the water vapor masses responsible for the precipitation and aquifers recharge seems to be present, based on strong correlations between 18O and 2H content and similar deuterium excess content. Yet, it is well noted that there is an isotopic gap between the cold water mean composition (δ18Omean = −5.25 ± 0.26‰ and δ2Hmean = −28.6 ± 1.6‰) and the thermomineral waters (δ18Omean = −6.5 ± 0.06‰ and δ2Hmean = −36.8 ± 0.3‰). This difference is in the order of 1.2‰ in 18O and around 8.2‰ in 2H. No evidence of evaporation previous to infiltration is noted in the isotopic content, based on the graphic distribution of all water samples along the Local Meteoric Water Line and in the deuterium excess values (d = 12.1‰ in AS10 to 15.6‰ in AS8); in addition, a common water vapor mass is assumed to be present. In the estimation of the mean isotopic composition of the cold groundwater samples, the spring AS1 was excluded. The isotopic difference observed between these two groups should be an altitude effect.
Several authors [4,6,23,24,27,49,50,51] use this isotopic shift as a natural tracer for the natural recharge of the water systems and in the establishment of the preferential recharge altitudes of the aquifer systems. Since one of the aims of this study is the identification of any differences/similarities in isotopic input to the water systems, thermomineral and cold dilute aquifers would lead to characterization of its isotopic signatures (mean isotopic composition). In São Pedro do Sul, this “gap” can be used to estimate the preferential recharge altitude of the thermomineral water system.
In general, the distribution of δ18O and δ2H mimics the topography of the continents; the waters of mountain chains is marked by more negative δ values. This feature is known as the “altitude effect” and was first mention by Dansgaard (1964) [52,53]. The lowering of temperature with increasing elevation in mountain regions usually leads to enhanced condensation and therefore to a progressive depletion in heavy isotopes in precipitation with altitude (Figure 9) [54].
The altitude effect has been used in numerous hydrological conceptual models all over the world, including on continental Portugal and Madeira Island. The groundwater isotopic composition has been applied as a basis of preferential recharge areas and definition of protection of aquifer systems, and to investigate the origin and to identify the interconnection of water bodies [6,22,25,26,55,56]. The vertical isotope gradient in mid-latitudes generally varies between 0.15‰ and 0.50‰ per each 100 m of altitude gained for oxygen-18, and about 1‰ to 4‰ per each 100 m for deuterium [57,58,59].
At São Pedro do Sul, the assessment of the elevation of the recharge area of the thermomineral system was established based on the δ18O values of the shallow cold aquifer systems that are characterized by a short residence time and short circulation pathway. The several springs sampled in the region are located at different altitudes. The results supported the relationship between δ18O of these cold shallow springs and the elevation of the sampling site, bearing in mind the conservative behavior of the oxygen isotopic composition of the groundwater from recharge to discharge. The δ18O isotopic gradient obtained for the research area was 0.19‰/100 m of altitude.
Taking into account the isotopic gradient of the region and the isotopic signatures (δ18O) of the thermomineral waters located at altitudes higher than 1000 m a.s.l., it is suggested that the Freita/Arada mountains, with a maximum altitude of 1120 m a.s.l. and located NW of the thermomineral poles, represent the preferential recharge area of this thermomineral system (Figure 2). This compatibility is corroborated by the close relationship between the isotopic composition (δ18O ) of São Pedro do Sul thermomineral waters and AS1 sample at an altitude of 1000 m a.s.l. (Figure 2 and Table 1).
In addition to the preferential recharge altitude of the thermomineral water, the isotopic composition of these samples also gives information concerning the water flux and potential water–rock interaction processes present along the groundwater flow, such as, for example, dissolution/precipitation processes. Formulating the hypothesis of the isotopic gap observed in the two groups of samples is due to dissolution/precipitation processes. The relation between oxygen-18 and deuterium should not be observed, and the water samples should not be plotted along the Local Meteoric Water Line. This isotopic pattern corroborates that the altitude effect is the main factor responsible for the observed isotopic distribution.

5. Final Remarks

The study performed at the São Pedro do Sul thermomineral pole is the continuity of previous works that have been performed in the region by different authors [1,7,40] in order to understand the conceptual model of these particular waters and their protection. Their knowledge and protection is important since the spa has a strong influence in the regional economy, inferred by the numbers of visitors and by direct and indirect revenue [1].
From the geochemical point of view, according to the Piper diagram, the groundwaters can be classified into two main groups: HCO3-Na and Cl-Na types. Cl–Na water types are generally associated with less mineralized waters (shallow cold diluted systems of local circulation), while HCO3-Na water types are associated with more mineralized waters (São Pedro do Sul thermomineral waters). In the evolutionary context of the typology of Cl-Na to HCO3-Na, the sulphureous thermomineral waters of the region represent the maximum of this evolution where HCO3 and Na ions are the most representative of ionic composition. This evolution results from hydrolysis of silicate mineral processes that occur over time. The four sources of the Termas and Vau poles are similar waters, not only by the mineralization values of around 365 mg/L, but also in relation to their relative composition of cations and anions. These results suggest a single aquifer system, a hypothesis that is corroborated by the reservoir temperature using the K2/Mg and SiO2 geothermalmeters (119 °C—maximum temperature), and homogenous maximum depths of the aquifer system (3.18 km—maximum depth average).
The δ18O and δ2H values indicate that the thermal sources mainly originated from meteoric water, with no significant evaporation before infiltration, and they are mainly recharged from precipitation at altitudes higher than 1000 m a.s.l. In the geohydraulic flow, this suggests a recharge area in the Freita/Arada mountains located NW of the thermomineral poles. The circulation of fluid is believed to happen through the WNW–ESE to NW–SE faults along several kilometers, which flow down into the crust where the water is warmed by the increasing natural temperature and then seeking pathways to reach the surface at the springs. The emergence is controlled mainly by the São Pedro do Sul-Ribamá, Termas and Fataunços faults. The flow-up process in the hydraulic circuit appears to be promoted by the difference in density between cold water inflow and hot water upflow and due to the difference in hydraulic charge between the recharge and discharge zones. Further studies are required to confirm the conclusion we have reached in this work.

Author Contributions

Conceptualization, S.A., L.G., A.O. and P.C.; methodology, S.A., L.G., A.O. and P.C.; software, S.A., A.O. and P.C.; validation, S.A., L.G., A.O. and P.C.; formal analysis, S.A., L.G., A.O. and P.C.; investigation, S.A., L.G., A.O. and P.C; resources, S.A., L.G., A.O. and P.C; data curation, S.A. and L.G.; writing—original draft preparation, S.A., L.G., A.O. and P.C; writing—review and editing, S.A., L.G., A.O. and P.C.; visualization, S.A., L.G., A.O. and P.C; supervision, S.A., L.G., A.O. and P.C.; project administration, S.A. and L.G.; funding acquisition, S.A., L.G., A.O. and P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research did not received specific funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge the Termalistur—Termas de São Pedro do Sul, E.M., S.A. and to the Geobiociências, Geoengenharias e Geotecnologias (Geobiotec) for the financial support in this study. The CGeo/UTAD author acknowledges the FCT (Portuguese Foundation for Science, I. P.) support through the UIDB/00073/2020 and UIDP/00073/2020 projects of the I&D Unit Geosciences Center (CGEO), and the C2TN/IST author acknowledges the FCT support through the UIDB/04349/2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ferreira Gomes, L.M.; Albuquerque, F.; Fresco, H. Protection areas of the São Pedro do Sul Spa, Portugal. Eng. Geol. 2001, 60, 341–349. [Google Scholar] [CrossRef]
  2. Domenico, P.A.; Schwartz, F.W. Physical and Chemical Hydrogeology; John Wiley & Sons: New York, NY, USA, 1990; p. 824. [Google Scholar]
  3. Calado, C. A ocorrência de água sulfúrea alcalina no maciço hespérico: Quadro hidrogeológico e quimiogénese. Ph.D. Thesis, Lisbon University, Lisbon, Portugal, 2001. [Google Scholar]
  4. Sousa Oliveira, A. Hidrogeologia dos Sistemas Gasocarbónicos da Província Hidromineral Transmontana: Ribeirinha (Mirandela, Sandim (Vinhais), Segirei e Salgadela (Chaves). Ph.D. Thesis, Trás os Montes e Alto Douro University, Vila Real, Portugal, 2001. [Google Scholar]
  5. Morais, M. Sistemas hidrominerais nos terrenos graníticos da zona centro-ibérica em Portugal central: Perspetivas químicas, isotópicas e genéticas sobre as águas sulfúreas. Ph.D. Thesis, Coimbra University, Coimbra, Portugal, 2012. [Google Scholar]
  6. Marques, J.M.; Graça, H.; Eggenkamp, H.G.M.; Neves, O.; Carreira, P.M.; Matias, M.J.; Mayer, B.; Nunes, D.; Trancoso, V.N. Isotopic and hydrochemical data as indicators of recharge áreas, flow paths and water-rock interaction in the Caldas da Rainha —Quinta das Janelas thermomineral carbonate rock aquifer (Central Portugal). J. Hydrol. 2013, 476, 302–313. [Google Scholar] [CrossRef]
  7. Almeida, S.; Ferreira Gomes, L.M.; Carvalho, P.; Oliveira, A. Contribution to the knowledge about hydrogeochemistry of São Pedro do Sul thermomineral system. Comunicações Geológicas 2015, 102, 125–128. [Google Scholar]
  8. Taillefer, A.; Guillou-Frottier, L.; Soliva, R.; Magri, F.; Lopez, S.; Courrioux, G. Topographic and faults control of hydrothermal circulation along dormant faults in an orogen. Geochem. Geophys. Geosystems 2018, 19, 4972–4995. [Google Scholar] [CrossRef] [Green Version]
  9. Martínez-Florentino, A.K.; Esteller, M.; Domínguez-Mariani, E.; Expósito, J.L.; Paredes, J. Hydrogeochemistry Isotopes and Geothermometry of Ixtapan de la Sal–Tonatico Hot Springs, Mexico. Environ. Earth Sci. 2019, 78, 600. [Google Scholar] [CrossRef]
  10. Barberio, M.D.; Gori, F.; Barbieri, M.; Boschetti, T.; Caracausi, A.; Cardello, G.L.; Petitta, M. Understanding the Origin and Mixing of Deep Fluids in Shallow Aquifers and Possible Implications for Crustal Deformation Studies: San Vittorino Plain, Central Apennines. Appl. Sci. 2021, 11, 1353. [Google Scholar] [CrossRef]
  11. Fronzi, D.; Mirabella, F.; Cardellini, C.; Caliro, S.; Palpacelli, S.; Cambi, C.; Valigi, D.; Tazioli, A. The Role of Faults in Groundwater Circulation before and after Seismic Events: Insights from Tracers, Water Isotopes and Geochemistry. Water 2021, 13, 1499. [Google Scholar] [CrossRef]
  12. Li, C.; Zhou, X.; Yan, Y.; Ouyang, S.; Liu, F. Hydrogeochemical Characteristics of Hot Springs and Their Short-Term Seismic Precursor Anomalies along the Xiaojiang Fault Zone, Southeast Tibet Plateau. Water 2021, 13, 2638. [Google Scholar] [CrossRef]
  13. Soares, C.A.R. Estudo Hidrogeológico das Termas de S. Pedro do Sul; Internal Report; Direção Geral de Geologia e Minas: Lisbon, Portugal, 1981.
  14. Haven, H.; Konings, R.; Schoonen, M.; Jansen, J.; Vriend, S.; Van Der Weijden, C.; Buitenkamp, J. Geochemical studies in the drainage basin of the Rio Vouga (Portugal): II. A model for the origin of hydrothermal water in the Vouzela region. Chem. Geol. 1985, 51, 225–238. [Google Scholar] [CrossRef]
  15. Cavaco, A. Projecto geotérmico de S. Pedro do Sul: Fase 1, Estudo Prévio de viabilidade. I item, Final Report; Internal Report; Câmara Municipal de São Pedro do Sul: São Pedro do Sul, Portugal, 1991. [Google Scholar]
  16. Cavaco, A. Definição do Perímetro de Proteção ao Aquífero de S. Pedro do Sul, Final Report; Internal Report Vol. I; Câmara Municipal de São Pedro do Sul: São Pedro do Sul, Portugal, 1995. [Google Scholar]
  17. Nolasco, S.M. Relatório de Trabalhos Geológicos e Geofísicos Efectuados nos Polos do Vau e das Termas; Projecto Geotérmico de S. Pedro do Sul, D.G.G.M.-Direção Geral de Geologia e Minas: Lisboa, Portugal, 1991; Volume II.
  18. Serviço de Fomento Mineiro. Projecto Geotérmico de S. Pedro do Sul: Relatório dos Trabalhos Geológicos e Geofísicos Efectuados nos Polos do Vau e das Termas. In Protocolo entre a Direcção-Geral de Geologia e Minas e a Câmara Municipal de S. Pedro do Sul; Serviço de Fomento Mineiro: Lisboa, Portugal, 1991. [Google Scholar]
  19. Ferreira, N.; Castro, P.; Godinho, M.M.; Neves, L.; Pereira, A.; Ferreira Pinto, A.; Simões, L.; Silva, F.; Aguado, B.; Azevedo, M.; et al. Folha 17-A Viseu. Carta geológica de Portugal à escala 1:50000; LNEG-Laboratório Nacional de Energia e Geologia, I. P.: Lisboa, Portugal, 2010. [Google Scholar]
  20. Ferreira Gomes, L.M.; Trota AP, N.; Sousa Oliveira, A.; Almeida, S.M.S. Reservoir considerations and direct uses of São Pedro do Sul hydromineral and geothermal field, northern Portugal. In Proceedings of the IOP Conference Series: Earth and Environmental Science, Prague, Czech Republic, 11–15 September 2017; p. 9. [Google Scholar] [CrossRef]
  21. Cox, S.F.; Knackstedt, M.A.; Braun, J. Principles of structural control on permeability and fluid flow in hydrothermal systems. Rev. Econ. Geol. 2001, 14, 1–24. [Google Scholar]
  22. Carreira, P.M.; Marques, J.M.; Monteiro Santos, F.A.; Andrade, M.; Matias, H.; Luzio, R.; Nunes, D. Role of geophysics, geochemistry and environmental isotopes in the assessment of Caldas de Monção low-temperature geothermal system. Geotherm. Resour. Counc. Trans. 2005, 29, 209–213. [Google Scholar]
  23. Carreira, P.M.; Marques, J.M.; Espinha Marques, J.; Chaminé, H.I.; Fonseca, P.E.; Monteiro Santos, F.; Moura, R.M.; Carvalho, J.M. Defining the dynamics of groundwater in Serra da Estrela Mountain area, central Portugal: An isotopic and hydrogeochemical approach. Hydrogeol. J. 2011, 19, 117–131. [Google Scholar] [CrossRef]
  24. Matiatos, I.; Alexopoulos, A.; Zouridakis, N. Use of stable isotopes in the determination of the mean altitude of recharge and the investigation of function mechanism of spring waters in Argolis Peninsula (Greece). Bull. Geol. Soc. Greece 2010, 43, 1792–1801. [Google Scholar] [CrossRef] [Green Version]
  25. Marques, J.M.; Matos, C.; Carreira, P.M.; Espinha Marques, J.; Teixeira, J.; Chaminé, H.I. Assessment of mixing between shallow and thermal waters using geochemical and environmental isotope tracers (N Portugal): A review and reinterpretation. Environ. Earth Sci. 2014, 72, 2557–2567. [Google Scholar] [CrossRef]
  26. Prada, S.; Virgílio Cruz, J.; Figueira, C. Using stable isotopes to characterize groundwater recharge sources in the volcanic island of Madeira, Portugal. J. Hydrology 2016, 536, 409–425. [Google Scholar] [CrossRef]
  27. Bahir, M.; Carreira, P.M.; Ouhamdouch, S.; Chamchati, H. Recharge conceptual model and mineralization of groundwater in a semi-arid region; Essaouira basin (Morocco). Procedia Earth Planet. Sci. 2017, 17, 69–72. [Google Scholar] [CrossRef]
  28. Scott, B.E.; Newell, D.L.; Jessup, M.J.; Grambling, T.A.; Shaw, C.A. Structural Controls on Crustal Fluid Circulation and Hot Spring Geochemistry Above a Flat-Slab Subduction Zone, Peru. Geochem. Geophys. Geosystems 2020, 21, e2020GC008919. [Google Scholar] [CrossRef]
  29. Apollaro, C.; Tripodi, V.; Vespasiano, G.; De Rosa, R.; Dotsika, E.; Fuoco, I.; Critelli, S.; Muto, F. Chemical Isotopic and geotectonic relations of the warm and cold waters of the Galatro and Antonimina thermal areas, southern Calabria, Italy. Mar. Pet. Geol. 2019, 109, 469–483. [Google Scholar] [CrossRef]
  30. Aydin, H.; Karaku¸s, H.; Mutlu, H. Hydrogeochemistry of Geothermal Waters in Eastern Turkey: Geochemical and Isotopic Constraints on Water-Rock Interaction. J. Volcanol. Geotherm. 2020, 390, 106708. [Google Scholar] [CrossRef]
  31. Caine, J.S.; Evans, J.P.; Forster, C.B. Fault zone architecture and permeability structure. Geology 1996, 24, 1025–1028. [Google Scholar] [CrossRef]
  32. Pereira, E.; Moreira, A.; Gonçalves, S.; Rodrigues, J.; Silva, A. Folha 13-D Oliveira de Azeméis. Carta geológica de Portugal à escala 1:50000; INETI-Instituto Nacional de Engenharia, Tecnologia e Inovação, I. P.: Lisboa, Portugal, 2006. [Google Scholar]
  33. Schermerhorn, L.; Sluijk, D.; Ramos, J.; Ávila Martins, J.; Fernandes, A.; Farinha Ramos, J. Folha 14-C Castro Daire, Carta Geológica de Portugal à Escala 1:50000; Serviços Geológicos de Portugal: Lisboa, Portugal, 1977. [Google Scholar]
  34. LNEG Laboratório Nacional de Energia e Geologia. Carta Geológica de Portugal, Folha 17-A (Viseu), Escala 1:50000; LNEG-Laboratório Nacional de Energia e Geologia, I. P.: Lisboa, Portugal, 2009. [Google Scholar]
  35. Pereira, E.; Rodrigues, J.; Gonçalves, L.; Moreira, A.; Silva, A. Noticia Explicativa da Folha 13-D Oliveira de Azeméis; INETI-Instituto Nacional de Engenharia, Tecnologia e Inovação, I. P.: Lisboa, Portugal, 2007. [Google Scholar]
  36. Schermerhorn, L. Notícia explicativa da folha 14 C Castro Daire; Serviços Geológicos de Portugal: Lisboa, Portugal, 1980. [Google Scholar]
  37. Ferreira, N.; Godinho, M.M.; Neves, L.; Pereira, A.S.; Castro, P.; Santos, T.B. Notícia Explicativa da Folha 17-A Viseu; LNEG-Laboratório Nacional de Energia e Geologia, I. P.: Lisboa, Portugal, 2010; ISBN 978-989-675-010-7. [Google Scholar]
  38. Serviços Geológicos de Portugal. Carta Geológica de Portugal à Escala 1:500000; Serviços Geológicos de Portugal: Lisboa, Portugal, 1992. [Google Scholar]
  39. IGeoE Instituto Geográfico do Exército. Cartas Militares de Portugal à Escala 1:25000, nº 165, 166, 176, 177; Série M888, Instituto Geográfico do Exército: Lisbon, Portugal, 1995. [Google Scholar]
  40. Pereira, E.; Ferreira, N. Geologia Regional e Controlo Estrutural das Nascentes de S. Pedro do Sul. Comun. Serv. Geol. Portugal 1985, 71, 17–25. [Google Scholar]
  41. Clarke, D.B. The mineralogy of peraluminous granites: A review. Can. Mineral. 1981, 19, 3–17. [Google Scholar]
  42. Marques, J.M.; Matias, M.J.; Basto, M.J.; Carreira, P.M.; Aires-Barros, L.A.; Goff, F.E. Hydrothermal alteration of Hercynian granites, its significance to the evolution of geothermal systems in granitic rocks. Geothermics 2010, 39, 152–160. [Google Scholar] [CrossRef]
  43. Barnicoat, A.; Sheldon, H.; Ord, A. Faulting and fluid flow in porous rocks and sediments: Implications for mineralisation and other processes. Miner. Depos. 2009, 44, 705–718. [Google Scholar] [CrossRef]
  44. Truesdell, A.H. Summary of Section III. Geochemical techniques in exploration. In Proceedings of the Second United Nations Symposium on the Development and Use of Geothermal Resources, San Francisco, CA, USA, 20 May 1975; pp. 53–79. [Google Scholar]
  45. Giggenbach, W.F. Geothermal solute equilibria-Derivation of Na-K-Ca-Mg geoindicators. Geochim. Et Cosmochim. Acta 1988, 52, 2759–2765. [Google Scholar] [CrossRef]
  46. IGM Instituto Geológico Mineiro. Recursos Geotérmicos em Portugal Continental: Baixa Entalpia. Obtido de LNEG. 1998. Available online: http://www.lneg.pt/CienciaParaTodos/edicoes_online/diversos/rec_geotermicos/texto (accessed on 15 April 2016).
  47. Piper, A.M. A Graphic Procedure in the Geochemical Interpretation of Water-Analyses. Trans. Am. Geophys. Union 1944, 914–928. [Google Scholar] [CrossRef]
  48. Craig, H. Isotopic Variations in Meteoric Waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef]
  49. Koeniger, P.; Toll, M.; Himmelsbach, T. Stable isotopes of precipitation and spring waters reveal an altitude effect in the Anti-Lebanon Mountains, Syria. Hydrol. Process. 2016, 30, 2851–2860. [Google Scholar] [CrossRef]
  50. Steinfeld, B.; Scott, J.; Vilander, G.; Marx, L.; Quirk, M.; Lindberg, J.; Koerner, K. Revision of the hydrogeological conceptual models of two Portuguese thermomineral water systems: Similarities and differences. Sustain. Water Resour. Manag. 2019, 5, 117–133. [Google Scholar] [CrossRef]
  51. Yang, J.; Large, R.R.; Bull, S.W. Factors controlling free thermal convection in faults in sedimentary basins: Implications for the formation of zinc– lead mineral deposits. Geofluids 2004, 4, 237–247. [Google Scholar] [CrossRef]
  52. Rozanski, K.; Araguás-Araguás, L.; Gonfiantini, R. Relation between long-term trends of oxygen-18 isotope composition of precipitation and climate. Science 1992, 258, 981–985. [Google Scholar] [CrossRef]
  53. Rozanski, K.; Araguás-Araguás, L.; Gonfiantini, R. Isotopic patterns in modern global precipitation. In Geophysical Monograph 78, Climate Change in Continental Isotopic Records; American Geophysical Union: Washington, DC, USA, 1993; pp. 1–36. [Google Scholar]
  54. Hoefs, J. Stable Isotope Geochemistry; Springer: Berlin/Heidelberg, Germany, 1987. [Google Scholar]
  55. Carreira, P.M.; Marques, J.M.; Carvalho, M.R.; Nunes, D.; Antunes da Silva, M. Carbon isotopes and geochemical processes in CO2-rich cold mineral water, N-Portugal. Environ. Earth Sci. 2014, 71, 2941–2953. [Google Scholar] [CrossRef]
  56. Carreira, P.M.; do Rosário Carvalho, M.; Nunes, J.C.; Grassa, F.; Capasso, G.; Marques, J.M. Preliminary Geochemical and Isotopic Results in Thermal and Cold Waters of Graciosa Volcanic Island (Azores). Procedia Earth Planet. Sci. 2017, 17, 630–633. [Google Scholar] [CrossRef]
  57. Araguás-Araguás, L.; Froehlich, K.; Rozanski, K. Deuterium and oxygen-18 isotope composition of precipitation and atmospheric moisture. Hydrol. Process 2000, 14, 1341–1355. [Google Scholar] [CrossRef]
  58. Gonfiantini, R.; Roche, M.A.; Olivry, J.C.; Fontes, J.C.; Zuppi, G.M. The altitude effect on the isotopic composition of tropical rains. Chem. Geol. 2001, 181, 147–167. [Google Scholar] [CrossRef]
  59. Carreira, P.M.; Nunes, D.; Valerio, P.; Araujo, M.F. A 15-year record of seasonal variation in the isotopic composition of precipitation water over continental Portugal. J. Radioanal. Nucl. Chem. 2009, 281, 153–156. [Google Scholar] [CrossRef]
Figure 1. (a) Geographic location of Portugal in Europe (Adapted from: https://commons.wikipedia.org (accessed on 3 April 2016)). (b) Geographic sketch location of the thermomineral poles of São Pedro do Sul.
Figure 1. (a) Geographic location of Portugal in Europe (Adapted from: https://commons.wikipedia.org (accessed on 3 April 2016)). (b) Geographic sketch location of the thermomineral poles of São Pedro do Sul.
Geosciences 12 00084 g001
Figure 2. Simplified regional geological and structural map of the study area, with main lithology units, faults and sampling locations.
Figure 2. Simplified regional geological and structural map of the study area, with main lithology units, faults and sampling locations.
Geosciences 12 00084 g002
Figure 3. (a) Extract from the geologic map of Portugal, 17-A sheet (Viseu), 1:50,000 scale, with the localization of the São Pedro do Sul thermomineral poles; (b) A–B and C–D cross sections.
Figure 3. (a) Extract from the geologic map of Portugal, 17-A sheet (Viseu), 1:50,000 scale, with the localization of the São Pedro do Sul thermomineral poles; (b) A–B and C–D cross sections.
Geosciences 12 00084 g003
Figure 4. Piper diagram for the 16 water samples of the São Pedro do Sul region. (a) Standing for the thermomineral waters; (b) Composed of sodium chlorinated waters; (c) Composed of sodium, calcium or magnesium bicarbonate waters.
Figure 4. Piper diagram for the 16 water samples of the São Pedro do Sul region. (a) Standing for the thermomineral waters; (b) Composed of sodium chlorinated waters; (c) Composed of sodium, calcium or magnesium bicarbonate waters.
Geosciences 12 00084 g004
Figure 5. Schoeller–Berkaloff diagram for the 16 water samples of São Pedro do Sul region. Parameters units: conductivity—μs cm–1; total mineralization and SiO2—mg/L; other parameters—meq/L.
Figure 5. Schoeller–Berkaloff diagram for the 16 water samples of São Pedro do Sul region. Parameters units: conductivity—μs cm–1; total mineralization and SiO2—mg/L; other parameters—meq/L.
Geosciences 12 00084 g005
Figure 6. (a) Na + K vs. HCO3; (b) Ca + Mg vs. HCO3 of São Pedro do Sul groundwater. (), thermomineral waters; (), shallow cold aquifer system boreholes; (), shallow cold aquifer system springs.
Figure 6. (a) Na + K vs. HCO3; (b) Ca + Mg vs. HCO3 of São Pedro do Sul groundwater. (), thermomineral waters; (), shallow cold aquifer system boreholes; (), shallow cold aquifer system springs.
Geosciences 12 00084 g006
Figure 7. (a) Cl vs. SiO2; (b) HCO3 vs. SiO2 and (c) Ca + Mg vs. SiO2 of São Pedro do Sul groundwaters. (), thermomineral waters; (), shallow cold aquifer system boreholes; (), shallow cold aquifer system springs.
Figure 7. (a) Cl vs. SiO2; (b) HCO3 vs. SiO2 and (c) Ca + Mg vs. SiO2 of São Pedro do Sul groundwaters. (), thermomineral waters; (), shallow cold aquifer system boreholes; (), shallow cold aquifer system springs.
Geosciences 12 00084 g007
Figure 8. δ18O vs. δ2H for the São Pedro do Sul thermomineral and cold dilute shallow aquifers. (), thermomineral waters; () shallow cold aquifer system springs.
Figure 8. δ18O vs. δ2H for the São Pedro do Sul thermomineral and cold dilute shallow aquifers. (), thermomineral waters; () shallow cold aquifer system springs.
Geosciences 12 00084 g008
Figure 9. Schematic variation of fractionation in δ18O during the hydrologic cycle; the altitude and continental effect in precipitation (water vapor masses) from coastal regions to inland.
Figure 9. Schematic variation of fractionation in δ18O during the hydrologic cycle; the altitude and continental effect in precipitation (water vapor masses) from coastal regions to inland.
Geosciences 12 00084 g009
Table 1. Summary of physical, chemical and isotopic parameters of shallow springs and well water samples. Alt., for altitude; EC, electrical conductivity; Temp., temperature values.
Table 1. Summary of physical, chemical and isotopic parameters of shallow springs and well water samples. Alt., for altitude; EC, electrical conductivity; Temp., temperature values.
SiteAlt.ECpHTemp.Total min.SiO2Na+Ca2+K+Mg2+HCO3ClSO42−Fδ18Oδ2H
mμs/cm °Cmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/L
Thermal Springs
NT *156400.288.8165.0364.5976.6090.903.013.250.03122.0527.4910.2417.98−6.44−36.80
AC1 *154399.298.8465.4365.4968.5591.152.953.240.02122.0927.719.4918.13−6.53−37.10
SDV1 *142401.118.8156.0371.7470.6291.283.063.230.02124.7628.5110.1118.09−6.45−36.60
SDV2 *142396.068.7862.7367.7370.5591.593.083.170.02123.9828.1910.2718.08−6.38−36.50
Cold Springs
AS1100014.505.8411.412.324.562.070.370.100.191.712.83<0.72<0.25−6.22−35.40
AS242529.505.6014.428.6410.784.700.260.750.485.256.35<0.72<0.25−5.13−28.50
AS335054.806.1715.851.9317.778.720.870.940.9211.5311.09<0.72<0.25−5.26−28.80
AS485025.205.3612.120.977.363.510.400.080.384.394.78<0.72<0.25−5.60−30.00
AS550559.906.0113.961.0320.899.261.571.550.7813.3610.361.48<0.25−5.43−28.40
AS644025.205.1715.621.218.673.630.330.090.203.175.11<0.72<0.25−5.16−28.40
AS772532.605.2413.331.4412.324.790.740.300.255.556.07<0.72<0.25−5.81−32.30
AS885012.7−5.51−28.50
AS970012.9−5.01−25.50
AS1027014.8−4.87−26.90
AS1127514.7−5.09−28.20
AS1232013.9−5.03−27.80
AS1344513.3−5.37−29.50
AS1451015.0−5.19−29.30
AS1525016.7−5.00−27.70
Boreholes
AF2 115.804.7018.081.5117.0912.551.853.652.20<0.5015.84<0.72<0.25
AF3 81.005.9822.784.7327.249.871.661.363.0027.579.213.92<0.25
AF4 33.105.2114.326.659.154.920.530.220.304.586.94<0.72<0.25
AF5 63.905.7512.852.069.167.103.750.901.2015.867.416.24<0.25
AF1A 88.305.4315.564.7714.3211.822.060.621.706.3414.058.40<0.25
* Mean values obtained over nine years (2005–2013) by physico-chemical and chemical analyses programs established by legal authorities, except for δ18O and δ2H values.
Table 2. Reservoir temperature using the K2/Mg and SiO2 geothermalmeters and maximum depth reached by São Pedro do Sul thermomineral waters.
Table 2. Reservoir temperature using the K2/Mg and SiO2 geothermalmeters and maximum depth reached by São Pedro do Sul thermomineral waters.
Ref.Geother. K2/Mg
(°C)
Maximum Depth (km)Geother. SiO2(°C)Maximum Depth (km)Average
Maximum Depth (km)
NT111.93.00122.83.333.17
AC1117.83.18117.23.163.17
SDV1117.73.17118.73.203.19
SDV2117.23.16118.63.203.18
Average116.13.13119.33.223.18
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Almeida, S.; Gomes, L.; Oliveira, A.; Carreira, P. Contributions for the Understanding of the São Pedro do Sul (North of Portugal) Geohydraulic and Thermomineral System: Hydrochemistry and Stable Isotopes Studies. Geosciences 2022, 12, 84. https://doi.org/10.3390/geosciences12020084

AMA Style

Almeida S, Gomes L, Oliveira A, Carreira P. Contributions for the Understanding of the São Pedro do Sul (North of Portugal) Geohydraulic and Thermomineral System: Hydrochemistry and Stable Isotopes Studies. Geosciences. 2022; 12(2):84. https://doi.org/10.3390/geosciences12020084

Chicago/Turabian Style

Almeida, Solange, Luís Gomes, Alcino Oliveira, and Paula Carreira. 2022. "Contributions for the Understanding of the São Pedro do Sul (North of Portugal) Geohydraulic and Thermomineral System: Hydrochemistry and Stable Isotopes Studies" Geosciences 12, no. 2: 84. https://doi.org/10.3390/geosciences12020084

APA Style

Almeida, S., Gomes, L., Oliveira, A., & Carreira, P. (2022). Contributions for the Understanding of the São Pedro do Sul (North of Portugal) Geohydraulic and Thermomineral System: Hydrochemistry and Stable Isotopes Studies. Geosciences, 12(2), 84. https://doi.org/10.3390/geosciences12020084

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop