The Origin of Dissolved Sulphate in the Thermal Waters of Budapest Inferred from Stable S and O Isotopes
by
,
István Fórizs
1,*
,
Viktória Rita Szabó
2,
József Deák
3,
Stanisław Hałas
4,†,
Andrzej Pelc
4,
Andrzej Trembaczowski
4 and
Árpád Lorberer
5,† 1
Institute for Geological and Geochemical Research, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences, H-1112 Budapest, Hungary
2
Department of Physical and Applied Geology, Eötvös Loránd University, H-1053 Budapest, Hungary
3
GWIS (Ground Water Isotope Studies) Ltd., H-8200 Veszprém, Hungary
4
Institute of Physics, M. Curie-Skłodowska University, 20-031 Lublin, Poland
5
Retired Hydrogeologist, Budapest, Hungary
*
Author to whom correspondence should be addressed.
†
Passed away.
Geosciences 2019, 9(10), 433; https://doi.org/10.3390/geosciences9100433
Submission received: 16 August 2019
/
Revised: 21 September 2019
/
Accepted: 23 September 2019
/
Published: 8 October 2019
(This article belongs to the Special Issue Isotope Geochemistry of Meteoric Waters)
Abstract
:The thermal waters produced by wells and springs from the Buda Thermal Karst in Budapest and its surroundings are rich in dissolved sulphate. Radiocarbon ages indicate that waters of T >45 ℃ were infiltrated during the Ice Age (more than 11 thousand years ago), on the higher elevations of the Buda-Pilis Hills, whereas waters of lower temperatures were infiltrated during the Holocene. For the origin of dissolved sulphate, two hypotheses can be set up: (1) the dissolved sulphate originates from the oxidation of the sulphide (pyrite) of Oligocene Clay Formation; (2) it is the dissolution product of the sulphate minerals (gypsum and anhydrite) of older carbonaceous rocks (limestone and dolomite). The isotopically stable sulphur composition of the dissolved sulphate in the thermal water (δ34S = 9.7‰ to 17.7‰) indicates its marine origin, so likely it dissolved from the older Permian evaporites.
1. Introduction
Budapest (Figure 1), known as city of spas, is very rich in thermal waters of various chemical compositions [1]. Thermal baths are an important treasure of the Hungarian capital city. The thermal waters are used mainly for balneological purposes; however, they are utilized as bottled water and for heating purposes as well. Most of these thermal waters are rich in dissolved sulphate [1,2]. In spite of the fact that the thermal waters of Budapest have been studied for more than one hundred years, nobody has examined in detail the origin of the dissolved sulphate, although some ideas were published about its origin [1,2,3]. One of the most reliable tools for the determination of the origin of sulphate is the application of environmentally stable isotopes; several papers demonstrated their usefulness (e.g., [4,5,6,7,8]). In our study, we would like to define the origin of the dissolved sulphate by means of the stable isotope measurements of dissolved sulphate (S and O) and water (O).
2. Hydrogeological Framework
The study area is the main discharge area of the Transdanubian Range. Both shallow (cold) and deep (warm) karstic flow systems end in the study area (Figure 2), and mixing between these two flow systems occurs at several places. At the Southern and Eastern foots of the Buda and Pilis Mountains (Figure 1), near the base level at the River Danube, there are many thermal (>30 °C) and lukewarm (from 15 °C to 30 °C) springs, which came from karstified Triassic and Eocene carbonaceous rocks [2]. The warmest component of water flow is in the basement of the Pest side (Figure 2), where karstified rocks are covered by thick sediment. The main mass of Buda and Pilis Mountains consist of well karstified Triassic platform carbonates. The reservoir of thermal waters of Budapest is mainly made of the Triassic karstified and fissured carbonaceous bedrock (limestone and dolomite). This Triassic bedrock is covered somewhere by Eocene carbonaceous rocks (Nummulite bearing limestone) [3]. The Upper Triassic bedrock is in hydraulic connection with the covering Upper Eocene formations in northern and southeastern Budapest [2]. The recharge area of Buda Thermalkarst is not defined exactly, but the discharge area can be delineated more exactly, because of its connection to tectonic lines. This concentrated discharge area is the Buda Spring-line at the right bank of the River Danube.
The springs emerge mainly along a north–south oriented path at the Buda side; some of them are in the bed of River Danube. The springs can be divided into two groups based on their chemistry and temperature: (1) the southern group of higher temperature and higher TDS (total dissolved solids), and (2) the northern group of lower temperature and lower TDS [2].
The Upper Triassic formations of the main thermal aquifer is deposited over older South-West Lower Triassic and Upper Permian rocks which contain evaporites (gypsum and anhydrite). These formations are known only from the Southwestern flank of the Budapest area, extending farther towards the Southwest. These limestones, marls and sandstones have low transmissivity and they are barely karstified, rather they are fissured crack water horizons. The Budapest flow system can probably reach only the very end of the known evaporite sediments. However, in the SW border of Budapest, at the Alcsútdoboz-Biatorbágy-Törökbálint-Budaörs-Budatétény line, they can probably communicate through a big horizontal surface with the more karstified Middle and Upper Triassic dolomites and limestones. The higher sulphate content of thermal water in South Budapest can be explained with that counterclockwise water flow system which contacts the Permian—Triassic evaporitic rocks (Figure 3). A more detailed description of the hydrogeological setting was published in 2017 by Mádl-Szőnyi et al. [10].
3. Samples and Analytical Methods
Twelve thermal water samples (>30 °C), four lukewarm water samples (between 15 °C and 30 °C) and eight cold-water samples (<15 °C) were taken from springs and wells in Budapest and its surroundings (Figure 1) The dissolved sulphate was precipitated from 1–3 L of water on the spot as BaSO4, by applying BaCl2. Beside the BaSO4, some BaCO3 was precipitated as well. The latter was removed by hydrochloric acid in the Institute for Geological and Geochemical Research, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences; meanwhile, the pH was strictly kept above 3, avoiding any oxygen isotopic exchange between water and sulphate. The pure BaSO4 samples were rinsed 5–8 times by deionized water, and then dried in oven at 60 °C. These samples of dry BaSO4 powder were mailed to the Mass Spectrometry Laboratory of Institute of Physics of Maria Curie-Skłodowska University (Lublin, Poland) for the determination of their δ34S and δ18O values. The oxygen isotopic composition of the water samples (0.6 dL) was determined in the Institute for Geological and Geochemical Research, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences (Budapest, Hungary). The isotopic compositions are given in the traditional delta notation as δ = (Rs/Rst − 1) × 1000‰, where Rs and Rst are the 18O/16O or 34S/32S isotopic ratios in the sample and standards, respectively. International standards: VSMOW (Vienna Standard Mean Ocean Water) for oxygen and CDT (Canon Diablo Troilite) for sulphur.
Both the sampling, and chemical analyses of the samples, were performed in the Laboratory of the Spas and Thermal Waters of Budapest Ltd, according to the related ISO standards.
4. Results
The temperature of the studied water samples varies from 8.9 °C to 76.8 °C. The lukewarm and cold-water samples were taken for the sake of comparison. The stable oxygen isotopic compositions of the waters (δ18O, Table 1) were in the range of −12.50‰ to −9.38‰ against VSMOW (mean value of all samples is −10.82‰). The mean value of the thermal samples is −11.68‰, which is a more negative value than that of lukewarm and cold-waters. The stable sulphur isotopic composition of the dissolved sulphate (δ34S, Table 1) was from −5.66‰ to 17.72‰ versus VCDT (mean value of all samples is 7.08‰). The δ34S values of sulphate varied from 9.68‰ to 17.72‰. Their mean value is 13.22‰, which is more positive than that of colder waters. The stable oxygen isotopic compositions of the dissolved sulphate (δ18O, Table 1) were from 0.68‰ to 6.39‰ (mean value of all samples is 4.18‰). The mean value of the thermal samples is 5.13‰, which is more positive than that of colder waters.
5. Discussion
The thermal wells and springs in Budapest and its surroundings are rich in dissolved sulphate. For the origin of this sulphate two hypotheses exist: (1) The dissolved sulphate originates from the oxidation of the sulphide (pyrite) of the Oligocene Kiscelli and Tard Clay Formations. This formation contains significant amounts of organic matter and sulphide; the latter is mostly in pyrite form. This pyrite can be oxidized where thermal water makes a contact with the clay. (2) The dissolved sulphate is a dissolution product of the sulphate minerals (gypsum, anhydrite) found in the Southwestern extension of the main Triassic karst limestone and dolomite aquifer [2].
5.1. Checking the Two Hypotheses Using Sulphate’s Sulphur Isotopes
5.1.1. The Oxidation of Sulphide
Vető et al. [11] published the δ34S values of total sulphur in the Tard Clay Formation (Alcsútdoboz-3 borehole). The mean value of the total sulphur of the clay was 2.6‰. Because sulphur in the clay is nearly 100% pyrite, the value δ34S = 2.6‰ can be regarded as that of pyrite. When pyrite is oxidized in abiogenic processes, the δ34S value of the sulphate formed is close to the value of original pyrite. When pyrite is oxidized in biogenic processes, the δ34S value of the sulphate formed is lower than that of pyrite. Since the mean δ34S value of the sulphate dissolved in the thermal waters studied is much more positive (average 13.6‰) than that of the pyrite in the Tard Clay Formation, (2.6‰), the first hypothesis that the dissolved sulphate originates from the oxidation of the sulphide (pyrite) in the Tard Clay Formation, is unlikely.
5.1.2. The Dissolution of Sulphate
Hámor [12], in his study about sedimentary sulphates, measured sulphur’s stable isotopic compositions mainly in Upper Permian and Lower Triassic evaporites. In these two groups, the isotopic compositions of the sulphur are significantly different. The δ34S values of the Permian samples were from 9.51‰ to 20.93‰ (mean value is 12.78‰; [13]). These Permian samples had been taken from several boreholes in Hungary (Bódvaszilas-1, Bódvarákó-4, Dinnyés-3, Hídvégardó-3, Jósvafő-2, Nagybátony-324, Nagyvisnyó-4, Nagyvisnyó-18, Perkupa-mine, Recsk-136 and Tornakápolna-3). The δ34S values measured for the samples mentioned roughly correspond with world values of Permian marine sulphates [14,15,16]. The δ34S values of the Triassic samples were from 15.98‰ to 33.01‰, with wide dispersion (mean δ34S value is 24.28‰), and are in good agreement with the global values [16]. These Triassic samples had been taken from several boreholes in Hungary (Alcsútdoboz-2, Alsószalmavár-1, Bakonyszűcs-1 and 3, Barcs-West-6, 8, Báta-3, Cún-1, Diósviszló-3, Forráskút-4, Iszkaszentgyörgy-4, Kiskunhalas-North-3, Kiskunhalas-South-5, Máriakéménd-3, Nagykozár-2, Ortaháza-West-1, Öttömös-2, Petend-1, Somberek-1, Tét-2, Turony-1, Vajta-3 and XIV.).
The results obtained prove our second hypothesis, that the dissolved sulphate is a dissolution product of evaporites, because in our study, the δ34S values of the thermal waters were between 9‰ and 18‰ (average is 13.22‰), which practically overlaps with the range of the above Permian evaporites; the latter is within the range of “Devonian to Lower Triassic” in Figure 4. There is a small overlap between the range of the measured δ34S values of thermal waters (9–18‰) and the above Triassic suphates (16–33‰), so some of the dissolved sulphate could originate from Triassic sources, but the amount must be inferior.
5.2. An Evaluation of the Chemical Data
The water chemistry may provide information about the sources of chemical species dissolved in the water (e.g., [18]). Here we provide a simple evaluation of chemical data just for checking whether the conclusion made on the basis of stable isotope data above is supported by the chemical data or not. In the previous chapter, we concluded that the dissolved sulphate in the thermal waters originates from the dissolution of sulphate minerals. Figure 5a indicates that there is a significant amount of surplus calcium compared to the congruent dissolution of sulphate (all the points are over the 1:1 dissolution line), necessarily so, because carbonate minerals have dissolved as well increasing the amount of Ca. But, if we compare the Ca with SO4 + HCO3 (Figure 5b) it is obvious that there is a shortage in Ca, because not only calcite dissolved, but dolomite as well. Regarding the model of dissolution of gypsum (and/or anhydrite), calcite and dolomite as sources for the Ca, Mg, SO4 and HCO3 species in the water, Figure 5c gives a support for it, showing a fairly good correlation between Ca + Mg and SO4 + HCO3. We may suppose that dissolution of halite is the source of chloride. Although the Na–Cl correlation is fairly good (Figure 5d) in many cases there is some surplus Na over the 1:1 line. This can be explained as the result of the Ca ↔ (Na,K) cation exchange between water and clay minerals. Clay releases Na or K and receives Ca ([19] chapter 4). It means that the Ca + Mg should be corrected with this surplus Na + K, which is over the dissolution of halite (Na + K–Cl). Figure 5e shows the corrected diagram, where the correlation is a little bit better than in the Figure 5c. Because not all the points are on the 1:1 line, we can infer that additional processes played roles, although minor ones. If we group the waters according to temperature into three intervals (Table 1), then there is no clear separation between the groups (Figure 5f), which can be the result of the mixing of cold and warm waters.
As we can see, the conclusion drawn by the stable isotope data of dissolved sulphate is supported by the chemical data.
5.3. Water Oxygen Isotopic Composition
The oxygen isotopic composition of the infiltrated water in Hungary is very close to that of the local precipitation [20]. Since the multi-annual stable isotopic characteristics of precipitation depends on the climate, the δ18O value of the infiltrated water reflects the climate in the time of infiltration. The δ18O values of the studied water samples (Table 1) show that the thermal waters (δ18O from −12.50‰ to −10.95‰) were infiltrated under a cooler climate during the Ice Age [21], while the cold-waters (δ18O from −10.67‰ to −9.38‰) were infiltrated during the Holocene [21]. The lukewarm waters (δ18O is from −10.58‰ to −10.24‰) are mixture of thermal and cold-waters. This is in harmony with δ34S values.
5.4. Oxygen Isotope, Bisulphate–Water Geothermometer
In our thermal water samples, δ18O values of the dissolved sulphate (bisulphate) were lower by 5–10‰ than Permian and Triassic marine sulphate values (12–18‰, Figure 3) [15,18,22]. There is an isotope exchange between the bisulphate and water which is rather fast over 200 ℃ [23] and significantly slower in lower temperatures; however, both pH and chemistry play important roles in determining the exchange rate. McKenzie and Truesdell [24] gave a formula for calculating the half time of oxygen exchange between bisulphate and water at different pHs and temperatures. We have calculated the time needed for 99.9% exchange of oxygen in the bisulphate-water system at pH 7 for different temperatures; see Table 2. Because the formula was given for pHs 3.8, 7 and 9, we made the calculations for pH 7, because it is the closest value to the pH we measured in the thermal waters studied (Table 3). Although the 99.9% exchange does not mean isotopic equilibrium, it can be regarded as a good approximation of the equilibrium.
Comparing the data in Table 2 with the water temperatures measured and the calculated radiocarbon ages in Table 3, we can see that the oxygen isotopic equilibrium in the bisulphate-water system has supposedly been attained in the majority of the waters except in four wells: Lukács IV and V, and Margitsziget III and IV wells, which are younger than nine-thousand years (ka), while the time needed for equilibrium is more than 9 ka at the outflow temperature.
Several authors dealt with the temperature dependence of oxygen isotope fractionation factor in the bisulphate-water system [23,25,26,27,28]. We tried the equations of Lloyd [25], Mizutani and Rafter [26] and Kusakabe and Robinson [23], and found that Equation (1), published by Kusakabe and Robinson [23], gave temperatures closest to the measured ones. They were the only authors who applied a correction for ion hydration effects, which is significant, and this can be a reason that their formula gave the closest temperatures.
Although Kusakabe and Robinson determined the O isotope fractionation factor of the BaSO4-H2O system, it can be regarded as the fractionation factor of the HSO4-H2O system, because the isotope fractionation between the solid BaSO4 and the dissolved HSO4 is negligible (e.g., [18] (p. 141)) both for S and O.
Comparing the measured outflow temperature with the temperature determined by the bisulphate-water geothermometer at waters of radiocarbon age >9 ka, we could see that thermometer T in three wells was equal to the measured T within ±2 °C (Széchenyi I and II wells; Zugló and Paskál malom wells); and in five wells it was significantly higher (Csepel II. well; VITUKI, Kvassay-well IX/38; Pesterzsébet, spa well; Margitsziget, Magda-well; Göd, strand thermal well) (Table 3). When the measured and the calculated geothermometer temperatures were ca. equal, then we may conclude that isotopic equilibrium was attained and neither mixing nor cooling took place; i.e. the water came from one flow path and its temperature was never significantly higher than the present one. Those wells are spatially close to each other and get their water from the ascending structures in faults close to each other. In the case that the calculated geothermometer T is significantly higher than the outflow temperature, theoretically, there are two explanations for this phenomenon. (1) Water comes from deep aquifers, where the temperature is high, and cools down during the ascension. The geothermometer gives the temperature of the deep aquifer, because the time of ascension is too short for re-equilibration. (2) The other explanation can be the mixing of small amount of cool water of low bisulphate content to the ascending hot water. In this case, the vast majority of the bisulphate originates from the hot component, but a small portion may come from the cool component, slightly modifying the δ18OHSO4. It is true for the δ18OH2O as well; it can be slightly modified by the cool component. The δ18O values of cool waters were higher than of the hot waters (Table 1), while the δ18O values of bisulphate in cool waters were lower than in the hot waters (Table 1). Both an increase in δ18OH2O and a decrease in δ18OHSO4 increase the geothermometer temperature calculated. At the same time, the effect of mixing of a small amount of young water with old water is very remarkable on the calculated radiocarbon age, resulting in a significantly younger age. (It is because a small increase in the 14C content results a great decrease in the radiocarbon age; the correlation between the age and 14C content is exponential.) Keeping in mind the above statements and having a look at the data in Table 3, we can see the following: in the case when the geothermometer T and the outflow T are close to each other, the T is in the range of 70–76 °C. The same geothermometer T characterizes some other wells: Csepel II well and Margitsziget Magda-well. Both are rather old (23.5 and 15.5 ka), and the outflow Ts are significantly lower than the geothermometer Ts, although for the latter, the difference is smaller. At these two wells, we may interpret the data such that the ascending water cooled down without mixing; therefore, the geothermometer gave the highest temperatures these waters had along their flow path.
There are three wells where the geothermometer T is just a little bit higher (79.7–79.9 °C) than the previously mentioned ones. These are: the VITUKI, Kvassay-well IX/38; the Pesterzsébet, spa well; the strand thermal well (Göd). All these waters are rather old waters, 14.6, 15.8 and 17.1 ka cal. 14C. It may be that a small amount of cool water mixes into the hot water, but according to the data it was insignificant, if it took place. In all the other cases, the geothermometer T is over 80 °C, and at the same time the radiocarbon age is under 7 ka: Lukács IV–V wells; Margitsziget III well. Although the geothermometer temperatures can be realistic, because at Tura (in the most eastern part of the Buda Thermal Karst) ca. 100 °C was measured in a lately drilled well, the parallel increase in geothermometer T and significant decrease in the radiocarbon age clearly indicates the mixing of young and cool water to the old and hot water.
5.5. The Possible Effect of S2− Content on the δ18OHSO4
In six thermal water samples, we detected a little amount of S2− (from 0.2 mg/L to 0.98 mg/L, Table 4), which may indicate bacterial reduction of sulphate, which involves isotope fractionation. However, the bacterial reduction of sulphate has not been proven, and other process, such as inorganic sulphate reduction, could produce that sulphide content [18] (p. 146); but because the latter is very rare, and the most common origin of S2- is the bacterial reduction of sulphate, we calculated the potential effect of bacterial sulphate reduction on the dissolved sulphate. We do not go into the details of the hypothetical process.
The bacterial reduction (for example, by Desulfovibrio desulfuricans, which is a natural anaerobic bacterium species) of sulphate removes, preferentially, the light isotopes of oxygen and sulphur from the sulphate in the reservoir and causes isotopic enrichments for both 18O and 34S in the residual sulphate [29]. The typical values for the ratio of enrichment factors ε34S/ε18O ranged from 2.5 to 4.5 [29]. Supposing that the bacterial reduction took place at enrichment factors of ε34S = 70‰ and ε18O = 20‰, we calculated what the δ18O and δ34S values of bisulphate were before the reduction. As we can see from Table 4, the possible effect of bacterial reduction on the δ values is insignificant; practically, this effect is in the range of the analytical uncertainty.
6. Conclusions
Thermal waters produced by wells and springs from the Buda Thermal Karst Aquifer of Budapest and its surroundings are rich in dissolved sulphate. The origin of this bisulphate has been determined by means of isotope geochemical methods.
Based on the stable sulphur isotopic composition of the dissolved sulphate we can infer that bisulphate in the thermal wells (>30 °C) originates from dissolution of Permian evaporates containing gypsum and anhydrite. These waters are old and were infiltrated during the Ice Age (>11 ka BP) with few exceptions, to which a small amount of cool water mixed, mostly affecting their radiocarbon age. The latter is below 9 ka (14C calculated age).
Luke-warm waters (T = 15–30 °C) are probably a mixture of cool and hot karstic water components based on the oxygen isotopic composition of the waters. The bisulphate originates supposedly from the oxidation of pyrite from the Tard and Kiscelli Clay Formations based on the stable sulphur isotopic composition.
An evaluation of water chemistry of thermal karstic wells suggest that bisulphate is a dissolution product of sulphate minerals; the bicarbonate is a dissolution product of calcite and dolomite; the chloride comes from NaCl; some portion of Na content comes from cation exchange between clay minerals and water [Ca ↔ (Na,K)].
In the waters of >60 °C the bisulphate-water oxygen isotope thermometer gives the same temperature as the outflow temperature, indicating that these waters are undisturbed waters and their temperature was never higher than the present one. In the thermal waters with lower temperatures, the bisulphate-water oxygen isotope thermometer outputs a slightly higher temperature those of the previous ones, indicating that a small amount of cool water could mix with this water. In the thermal waters younger than 9 ka, the bisulphate-water oxygen isotope thermometer gives temperatures of over 80 °C, which can be real, but their age and δ18OH2O indicates the mixing of young cool water to old hot water.
Author Contributions
Conceptualization, I.F., J.D. and S.H.; field work, J.D.; isotope analysis, A.P., A.T. and I.F.; data curation, V.R.S., J.D. and I.F.; writing—original draft preparation, V.R.S. and I.F.; writing—review and editing, A.P., Á.L. and A.T.; visualization, V.R.S.; supervision of sulphate isotope measurements, S.H.; funding acquisition, I.F.
Funding
This research was funded by the Hungarian Scientific Research Fund, grant number K-60921, and the National Research, Development and Innovation Office, grant number SNN118205.
Acknowledgments
The authors would like to thank, for the water oxygen isotope measurements of a set of the samples, Krisztina Kármán.
Conflicts of Interest
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
References
- Szalontai, G. A budapesti hévizek kémiai tulajdonságai (Chemistry of the Budapest Thermal Waters). In Budapest hévizei (Thermal Waters of Budapest); Alföldi, L., Bélteky, L., Böcker, T., Horváth, J., Korim, K., Liebe, P., Rémi, R., Eds.; VITUKI: Budapest, Hungary, 1968; pp. 59–71. (In Hungarian) [Google Scholar]
- Lorberer, Á. Budapest hévizei mérnökgeológiai szemmel (Thermal waters of Budapest from the aspect of engineering geology). In Alagút- és mélyépítő szakmai napok, Proceedings of the a Millenium Után, Európával, jövőnk környezetéért, Eger, Hungary, 27–28 May 2002; MTESZ Közlekedéstudományi Egyesület—Deutsche Ingeinurkammer: Budapest, Hungary, 2002; pp. 71–78. (In Hungarian) [Google Scholar]
- Alföldi, L. Budapesti hévizek (Thermal waters of Budapest)—VITUKI Közlemények, 20.; VITUKI: Budapest, Hungary, 1979; pp. 1–102. (In Hungarian) [Google Scholar]
- Krouse, H.R.; Mayer, B. Sulphur and oxygen isotopes in sulphate. In Environmental Tracers in Subsurface Hydrology; Cook, P.G., Herczeg, A.L., Eds.; Springer: Boston, MA, USA, 2019; pp. 195–231. [Google Scholar]
- Apollaro, C.; Dotsika, E.; Marini, L.; Barca, D.; Bloise, A.; De Rosa, R.; Doveri, M.; Lelli, M.; Muto, F. Chemical and isotopic characterization of the thermo mineral water of Terme Sibarite springs (Northern Calabria, Italy). Geochem. J. 2012, 46, 117–129. [Google Scholar] [CrossRef]
- Michalik, A.; Migaszewski, Z.M. Stable sulfur and oxygen isotope ratios of the Świętokrzyski National Park spring waters generated by natural and anthropogenic factors (south-central Poland). Appl. Geochem. 2012, 27, 1123–1132. [Google Scholar] [CrossRef]
- Migaszewski, Z.; Galuszka, A.; Michalik, A.; Dolegowska, S.; Migaszewski, A.; Halas, S.; Trembaczowski, A. The use of stable sulfur, oxygen and hydrogen isotope ratios as geochemical tracers of sulfates in the podwisniowka acid drainage area (South-Central Poland). Aquat. Geochem. 2013, 19, 261–280. [Google Scholar] [CrossRef]
- Porowski, A.; Porowska, D.; Halas, S. Identification of sulfate sources and biogeochemical processes in an aquifer affected by Peatland: Insights from monitoring the isotopic composition of groundwater sulfate in Kampinos National Park, Poland. Water 2019, 11, 1388. [Google Scholar] [CrossRef]
- Vendel, M.; Kisházi, P. Összefüggések melegforrások és karszt-vizek között a Dunántúli-középhegységben megfigyelt viszonyok alapján I and II (Correlations between thermal springs and Karstic waters in the Transdanubian Ranges I and II). MTA Műszaki Tudományok Osztályának Közleményei 1964, 32–33, 205–234, 393–417. (In Hungarian) [Google Scholar]
- Mádl-Szőnyi, J.; Erőss, A.; Tóth, Á. Fluid flow systems and hypogene karst of the transdanubian range, Hungary—With special emphasis on buda thermal karst. In Hypogene Karst Regions and Caves of the World; Klimchouk, A., Palmer, A.N., De Waele, J., Auler, A.S., Audra, P., Eds.; Springer: Cham, Switzerland, 2017; pp. 267–278. [Google Scholar]
- Vető, I.; Nagymarosy, A.; Brukner-Wein, A.; Hetényi, M.; Sajgó, C. Salinity changes control isotopic composition and preservation of the organic matter: The Oligocene Tard clay, Hungary revisited. In Proceedings the of 19th International Organic Geochemistry, Istanbul, Turkey, 6–10 September 1999; Abstracts Part I, TÜBITAK; Marmara Research Center: Istanbul, Turkey, 1999; pp. 411–412. [Google Scholar]
- Hámor, T. Az anoxikus üledékképződés és a korai diagenezis vizsgálata stabil izotóp mérések alkalmazásával (Study of anoxic sedimentation and early diagenesis by means of stable isotopes). In Országos Tudományos Kutatási Alap 259. sz. téma zárójelentése (Final Report of Project no. 259 Financed by the National Scientific Research Fund); Hungarian Geological Institute (Now Mining and Geological Survey of Hungary): Budapest, Hungary, 1991. (In Hungarian) [Google Scholar]
- Hámor, T. Stabil izotóp mérések alkalmazása szedimentológiai és környezetvédelmi problémák megoldásában (Application of stable isotope measurement for sedimentological and environmental problems). In Az Országos Tudományos Kutatási Alap F007373. sz. kutatásának zárójelentése ((Final Report of Project no. F007373 Financed by the National Scientific Research Fund); Mining and Geological Survey of Hungary: Budapest, Hungary, 1997. (In Hungarian) [Google Scholar]
- Newton, R.J.; Pevitt, E.L.; Wignall, P.B.; Bottrell, S.H. Large shifts in the isotopic composition of seawater sulphate across the Permo-Triassic boundary in northern Italy. Earth Planet. Sci. Lett. 2004, 218, 331–345. [Google Scholar] [CrossRef]
- Claypool, G.E.; Holser, W.T.; Kaplan, I.R.; Sakai, H.; Zak, I. The age curves of sulfur and oxygen isotopes in marine sulfate and their mutual interpretation. Chem. Geol. 1980, 28, 199–260. [Google Scholar] [CrossRef]
- Cortecci, G.; Reyes, E.; Berti, G.; Casati, P. Sulfur and oxygen isotopes in Italian marine sulfates of Permian and Triassic ages. Chem. Geol. 1981, 34, 65–79. [Google Scholar] [CrossRef]
- Clark, I.D.; Fritz, P. Environmental Isotopes in Hydrogeology; Lewis Publishers: Boca Raton, FL, USA, 1997. [Google Scholar]
- Dogramaci, S.; McLean, L.; Skrzypek, G. Hydrochemical and stable isotope indicators of pyrite oxidation in carbonate-rich environment; the Hamersley Basin, Western Australia. J. Hydrol. 2017, 545, 288–298. [Google Scholar] [CrossRef]
- Drever, J.I. The Geochemistry of Natural Waters. Surface and Groundwater Environments, 3rd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 1997. [Google Scholar]
- Fórizs, I.; Czuppon, G.; Kern, Z.; Kohán, B.; Deák, J. Relationship between the isotopic characteristics of local precipitation and groundwater in Hungary. In Book of Abstracts—44th Annual Congress of International Association of Hydrogeologists. “Groundwater Heritage and Sustainability”, Proceedings of the International Association of Hydrogeologists (IAH), Dubrovnik, Croatia, 25–29 September 2017; Posavec, K., Marković, T., Eds.; Croatian National Chapter of the International Association of Hydrogeologists: Dubrovnik, Croatia, 2017; p. 258. [Google Scholar]
- Deák, J.; Coplen, T.B. Identification of Holocene and Pleistocene groundwaters in Hungary using oxygen and hydrogen isotopic ratios. In Isotopes in Water Resources Management (Symposium Proceedings, Vienna, 1995); International Atomic Energy Agency: Vienna, Austria, 1996; Volume 1, p. 438. [Google Scholar]
- Longinelli, A.; Flora, O. Isotopic composition of gypsum samples of Permian and Triassic age from the north-eastern Italian Alps: Paleoenvironmental implications. Chem. Geol. 2007, 245, 275–284. [Google Scholar] [CrossRef]
- Kusakabe, M.; Robinson, B.W. Oxygen and sulfur isotope equilibria in the BaSO4-HSO4−-H2O system from 110 to 350 °C and applications. Geochim. et Cosmochim. Acta 1977, 41, 1033–1040. [Google Scholar] [CrossRef]
- McKenzie, W.F.; Truesdell, A.H. Geothermal reservoir temperatures estimated from the oxygen isotope compositions of dissolved sulfate and water from hot springs and shallow drillholes. Geothermics 1977, 5, 51–61. [Google Scholar] [CrossRef]
- Lloyd, R.M. Oxygen isotope behavior in the sulfate-water system. J. Geophys. Res. 1968, 73, 6099–6110. [Google Scholar] [CrossRef]
- Mizutani, Y.; Rafter, T.A. Oxygen isotopic composition of sulphates, 3. Oxygen isotopic fractionation in the bisulphate ion-water system. N. Zeal. J. Sci. 1969, 12, 54–59. [Google Scholar]
- Mizutani, Y. Isotopic composition and underground temperature of the Otake geothermal water, Kyushu, Japan. Geochem. J. 1972, 6, 67–73. [Google Scholar] [CrossRef]
- Chiba, H.; Sakai, H. Oxygen isotope exchange rate between dissolved sulphate and water at hydrothermal temperatures. Geochim. et Cosmochim. Acta 1985, 49, 993–1000. [Google Scholar] [CrossRef]
- Fritz, P.; Basharmal, G.M.; Drimmie, R.J.; Ibsen, J.; Qureshi, R.M. Oxygen isotope exchange between sulphate and water during bacterial reduction of sulphate. Chem. Geol. (Isotope Geosci. Sect.) 1989, 79, 99–105. [Google Scholar] [CrossRef]
Figure 1.
The location of the studied area in Hungary, and the sampling locations (blue dots) in Budapest and its surroundings.
Figure 1.
The location of the studied area in Hungary, and the sampling locations (blue dots) in Budapest and its surroundings.
Figure 2.
Conceptual flow model of the Buda Thermal Karst system from Vendel–Kisházi (1964) [9]. Abbreviations: T—Triassic dolomite; E—Eocene limestone and Eocene–Oligocene marl; O—Oligocene clay.
Figure 2.
Conceptual flow model of the Buda Thermal Karst system from Vendel–Kisházi (1964) [9]. Abbreviations: T—Triassic dolomite; E—Eocene limestone and Eocene–Oligocene marl; O—Oligocene clay.
Figure 3.
The distribution of the dissolved sulphate (mg/L) in the thermal waters of Budapest. Isolines (light blue) are interpolated from the point data by manual interpretation based on linear estimation (Á. Lorberer co-author). Key: blue dots—measured karst wells; blue circles—other karst wells having SO4-measurements; green lines—border between the unconfined mountainous part and the confined thermal part of the aquifer; grey lines—district and city borders; orange lines—main roads.
Figure 3.
The distribution of the dissolved sulphate (mg/L) in the thermal waters of Budapest. Isolines (light blue) are interpolated from the point data by manual interpretation based on linear estimation (Á. Lorberer co-author). Key: blue dots—measured karst wells; blue circles—other karst wells having SO4-measurements; green lines—border between the unconfined mountainous part and the confined thermal part of the aquifer; grey lines—district and city borders; orange lines—main roads.
Figure 4.
The isotopic composition of the dissolved sulphate in the thermal waters of Budapest and its surroundings on the δ34S–δ18O diagram in comparison with terrestrial, marine and atmospheric sulphates [17].
Figure 4.
The isotopic composition of the dissolved sulphate in the thermal waters of Budapest and its surroundings on the δ34S–δ18O diagram in comparison with terrestrial, marine and atmospheric sulphates [17].
Figure 5.
Graphs showing the relationship between the chemical components of the studied waters. Every component is in milliequivalent/litre (meq/L). The 1:1 line indicates the congruent dissolution(s) involved. (a) Ca versus SO4; (b) Ca versus SO4 + HCO3; (c) Ca + Mg versus SO4 + HCO3; (d) Na versus Cl; (e) and (f) Ca + Mg + (Na + K–Cl) versus SO4 + HCO3.
Figure 5.
Graphs showing the relationship between the chemical components of the studied waters. Every component is in milliequivalent/litre (meq/L). The 1:1 line indicates the congruent dissolution(s) involved. (a) Ca versus SO4; (b) Ca versus SO4 + HCO3; (c) Ca + Mg versus SO4 + HCO3; (d) Na versus Cl; (e) and (f) Ca + Mg + (Na + K–Cl) versus SO4 + HCO3.
Table 1.
Water temperature [°C], δ34S and δ18O values [‰] of dissolved sulphate and water, grouped into temperature ranges.
Table 1.
Water temperature [°C], δ34S and δ18O values [‰] of dissolved sulphate and water, grouped into temperature ranges.
| Well/Spring | δ34S [‰] Sulphate | δ18O [‰] Water | δ18O [‰] Sulphate | Twater [°C] |
|---|---|---|---|---|
| Thermal Waters (>30 °C) | ||||
| Csepel II well (Budapest) | 15.61 | −12.01 | 5.45 | 45.2 |
| Széchenyi II well (Budapest) | 17.18 | −12.46 | 4.91 | 76.8 |
| Széchenyi I well (Budapest) | 17.72 | −12.50 | 5.05 | 73.0 |
| Lukács IV well (Budapest) | 11.63 | −11.59 | 4.22 | 52.5 |
| Lukács V well (Budapest) | 10.28 | −11.36 | 5.18 | 50.9 |
| VITUKI, Kvassay-well IX/38 (Budapest) | 11.65 | −11.53 | 5.27 | 46.2 |
| Pesterzsébet, spa well (Budapest) | 11.32 | −11.60 | 5.21 | 41.8 |
| Zugló, Paskál malom well (Budapest) | 13.35 | −11.75 | 6.39 | 67.7 |
| Margitsziget, Magda-well (Budapest) | 16.12 | −11.88 | 5.36 | 69.0 |
| Margitsziget, IV well (Budapest) | 9.68 | −10.95 | 38.1 | |
| Margitsziget, III well (Budapest) | 10.19 | −10.99 | 4.47 | 37.8 |
| Strand thermal well (Göd) | 14.01 | −11.87 | 4.94 | 50.4 |
| Lukewarm Waters (15–30 °C) | ||||
| Lukács, Boltív spring (Malom Lake) (Budapest) | −4.91 | −10.48 | 2.64 | 22.0 |
| Lukács, Római spring (Budapest) | −5.66 | −10.30 | 2.87 | 20.8 |
| Lukács, Török spring (Budapest) | −3.60 | −10.58 | 23.9 | |
| Óbuda, Bründl spring (Budapest) | −2.47 | −10.24 | 18.5 | |
| Cold Waters (<15 °C) | ||||
| Vízművek-1 (3. works) well (Pilisborosjenő) | −0.22 | −9.38 | 4.12 | 14.4 |
| Vízművek-2 (4. works) well (Pilisborosjenő) | −1.15 | −9.96 | 13.3 | |
| Vízművek, Ipari well (Pilisszentiván) | 7.51 | −9.99 | 3.87 | 11.0 |
| Szentkút spring (Csobánka) | 1.32 | −9.81 | 9.4 | |
| Lajos spring (Szentendre) | 8.08 | −10.67 | 2.10 | 9.4 |
| Diósvölgyi spring (Süttő) | 8.32 | −9.70 | 12.4 | |
| Vízművek 2. well (Tarján) | 1.71 | −10.10 | 12.1 | |
| Török-well (Zsámbék) | 7.04 | −10.24 | 2.59 | 12.5 |
| Szentkút spring (Csobánka) | 2.37 | −9.74 | 0.68 | 8.9 |
Table 2.
Time needed for 99.9% exchange of O in the bisulphate-water system calculated at pH 7 using the formula of McKentzie and Truesdell [24]. a = year.
Table 2.
Time needed for 99.9% exchange of O in the bisulphate-water system calculated at pH 7 using the formula of McKentzie and Truesdell [24]. a = year.
| Temperature [°C] | 30 °C | 40 °C | 50 °C | 60 °C | 70 °C | 80 °C | 90 °C |
|---|---|---|---|---|---|---|---|
| Time needed for 99.9% exchange (10 × half time) [a] | 18,445 | 9962 | 5589 | 3246 | 1946 | 1201 | 761 |
Table 3.
The temperatures measured (Twater), thermometer temperatures (Twater calculated), oxygen isotopic composition of dissolved sulphate, pH and 13C-corrected radiocarbon ages of thermal waters.
Table 3.
The temperatures measured (Twater), thermometer temperatures (Twater calculated), oxygen isotopic composition of dissolved sulphate, pH and 13C-corrected radiocarbon ages of thermal waters.
| Thermal Wells | Twater [°C] | Twater [°C] Calculated | δ18O [‰] Sulphate | pH | 14C age (13C-corr.) [103 × a] |
|---|---|---|---|---|---|
| Csepel II well (Budapest) | 45.2 | 75.1 | 5.45 | 6.6 | 23.5 |
| Széchenyi II well (Budapest) | 76.8 | 75.7 | 4.91 | 6.4 | 18.9 |
| Széchenyi I well (Budapest) | 73.0 | 74.4 | 5.05 | 6.4 | 19.0 |
| Lukács IV well (Budapest) | 52.5 | 87.3 | 4.22 | 6.9 | 3.0 |
| Lukács V well (Budapest) | 50.9 | 81.8 | 5.18 | 6.9 | 4.7 |
| VITUKI, Kvassay-well IX/38 (Budapest) | 46.2 | 79.9 | 5.27 | 6.8 | 14.6 |
| Pesterzsébet, spa well (Budapest) | 41.8 | 79.8 | 5.21 | 6.7 | 15.8 |
| Zugló, Paskál malom well (Budapest) | 67.7 | 70.5 | 6.39 | 6.5 | 15.2 |
| Margitsziget, Magda-well (Budapest) | 69.0 | 76.6 | 5.36 | 6.7 | 15.5 |
| Margitsziget, IV well (Budapest) | 38.1 | 6.9 | 4.9 | ||
| Margitsziget, III well (Budapest) | 37.8 | 90.1 | 4.47 | 6.8 | 7.0 |
| Strand thermal well (Göd) | 50.4 | 79.7 | 4.94 | 6.4 | 17.1 |
Table 4.
The calculated δ18O and δ34S values of bisulphate before bacterial reduction in six sulphate samples. The enrichment factors we used in this calculation are ε34S = 70‰ and ε18O = 20‰. The ε34S/ε18O ratio is 3.5, which is the mean value of the range 2.5–4.5 defined by Fritz et al. (1989) [29].
Table 4.
The calculated δ18O and δ34S values of bisulphate before bacterial reduction in six sulphate samples. The enrichment factors we used in this calculation are ε34S = 70‰ and ε18O = 20‰. The ε34S/ε18O ratio is 3.5, which is the mean value of the range 2.5–4.5 defined by Fritz et al. (1989) [29].
| Thermal Wells | δ34S [‰] Sulphate (Measured) | δ18O [‰] Sulphate (Measured) | S2− [mg/L] | δ34S [‰] Sulphate (Calculated) | δ18O [‰] Sulphate (Calculated) |
|---|---|---|---|---|---|
| Csepel well II (Budapest) | 15.61 | 5.45 | 0.4 | 15.58 | 5.43 |
| Széchenyi well II (Budapest) | 17.18 | 4.91 | 0.2 | 17.16 | 4.90 |
| Széchenyi well I (Budapest) | 17.72 | 5.05 | 0.3 | 17.69 | 5.03 |
| Lukács well IV (Budapest) | 11.63 | 4.22 | 0.98 | 11.50 | 4.14 |
| Lukács well V (Budapest) | 10.28 | 5.18 | 0.64 | 10.16 | 5.11 |
| Strand thermal well (Göd) | 14.01 | 4.94 | 0.34 | 13.95 | 4.90 |
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Fórizs, I.; Szabó, V.R.; Deák, J.; Hałas, S.; Pelc, A.; Trembaczowski, A.; Lorberer, Á. The Origin of Dissolved Sulphate in the Thermal Waters of Budapest Inferred from Stable S and O Isotopes. Geosciences 2019, 9, 433. https://doi.org/10.3390/geosciences9100433
AMA Style
Fórizs I, Szabó VR, Deák J, Hałas S, Pelc A, Trembaczowski A, Lorberer Á. The Origin of Dissolved Sulphate in the Thermal Waters of Budapest Inferred from Stable S and O Isotopes. Geosciences. 2019; 9(10):433. https://doi.org/10.3390/geosciences9100433
Chicago/Turabian StyleFórizs, István, Viktória Rita Szabó, József Deák, Stanisław Hałas, Andrzej Pelc, Andrzej Trembaczowski, and Árpád Lorberer. 2019. "The Origin of Dissolved Sulphate in the Thermal Waters of Budapest Inferred from Stable S and O Isotopes" Geosciences 9, no. 10: 433. https://doi.org/10.3390/geosciences9100433
APA StyleFórizs, I., Szabó, V. R., Deák, J., Hałas, S., Pelc, A., Trembaczowski, A., & Lorberer, Á. (2019). The Origin of Dissolved Sulphate in the Thermal Waters of Budapest Inferred from Stable S and O Isotopes. Geosciences, 9(10), 433. https://doi.org/10.3390/geosciences9100433
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
