Thermal springs and balneology in the Peri-Adriatic area: geochemical status and prospects Riccardo Petrini, Pisa University (with contributions by R. Cataldi and B. Della Vedova) 1
The Peri-Adriatic Region (Picha F.J. 2002) Epicenters of seismic events of magnitude > 2.0, from 1975 to 2005, U.S. Geological Survey Earthquake Database (Del Gaudio et al., 2007) 2
Thermal Spas in Peri-Adriatic Region COUNTRY Number of active SPAs Water volume used in 2010 (m 3 /year, excluding pool uses) Albania 5 140.000 Bosnia Herzegovina 20 1.290.000 Croatia 18 1.345.000 Greece > 23 (up to 50?) > 500.000 Italy 149 48.000.000 Montenegro 1 60.000 Slovenia 25 1.780.000 T O T A L > 260 > 53.115.000 Including pool uses, a flow of about 112x10 6 m 3 /year (3.55 m 3 /s) of thermal waters can be estimated to feed Spas of the Adriatic-Jonian Region 3
Thermal Spas in Peri-Adriatic Region: Q,T, V, MW, TJ Country N. of SPAs Total average flow rate (kg/s) Annual Total average water used (10 6 *m 3 ) T ( C) Average T ( C) Resource potential capacity (MWt) Total geothermal energy used (TJ/yr) Estimated number of tourists (1) Source (2) Slovenia 26 115* 3.6 23-63 33-38 25 311 751000 (2010)? WGC 2010 Croatia 18 96* 3.0 25-85 42-47 33 216 400000 (2010)? WGC 2010 Bosnia and Herzegovina 20 79* 2.5 21-76 39-44 13 100 192000 (2009)? WGC 2010 Montenegro 1?????? WGC 2010 Albania 6 21 0.7 27-60 45-50 12 8? WGC 2010 Greece >60 >1500 47 18-100? 39?** 238?** >100000 (2010)? WGC 2010 Serbia 59 700-800? 23 29-96 45-50 40 647? WGC 2010 TOTAL > 191 > 2.561 79,8 162 1520 >1.443.000 (1) Data taken from different sources, they likely contain large uncertainties. (2) Proceedings World Geothermal Congress 2010 (*) Water volumes have been estimeted starting from the Annual Energy Use of each country. For each spring we estimated the average mean annual energy (based on the average flow) for the effective temperature intervals (T in T out ), using the Annual Energy Use (TJ/yr) = Ave. flow rate (kg/s) x [inlet temp. ( C) - outlet temp. ( C)] x 0.1319 (**) Figures largely underestimated 4
Energy Potential of thermal springs and waters Conceptual flow diagram for water and energy utilization in spa systems, excluding pools Conceptual flow diagram for water and energy utilization of thermal pools associated to spa systems Mass and energy balances are crucial for sustainable use of excess heat, up- and down-stream Spa uses, and are tools to develop the balneotherapy sector. 5
Thermal springs and waters in peri-adriatic context Geothermal reservoirs are often related to areas of active or recent magmatism, where heat is transferred from a magmatic source to the circulating fluids However, low-enthalpy waters in different tectonic settings, including aquifers in deep carbonate-rocks and metamorphic basements, are receiving much attention, in particular for bathing and swimming The Peri-Adriatic area has many favorable geothermal characteristics; however, long-term abstraction of resources should be carefully planned to avoid overexploitation 6
Water geochemistry: source aquifers for a sustainable exploitation policy The major ions chemistry gives the basic information on water-type and processes cations Ca 2+, Mg 2+, Na +, K + anions HCO 3-, Cl -, SO 4 2- Piper diagram 7
Geothermal manifestations in the Peri-Adriatic Region belong to different hydrofacies, reflecting different origin and nature of aquifers They include: Thermal waters in Mesozoic carbonate-rock aquifers Thermal waters in aquifers within the metamorphic basement Thermal waters in porous media in sedimentary basins Thermal waters of marine origin in coastal environments Each requires to enhance knowledge to reduce future quantitative and qualitative threats 8
Hydrofacies: mostly of the Ca-HCO 3 type Bosnia Some of the Ca-HCO 3 waters are displaced towards Na-Cl marine components seawater 9
Hydrofacies: waters range from the Ca- HCO 3 to the Na-Cl, Ca-SO 4 type Croatia limited ion-exchanges The data suggest mixing processes between the different water-types seawater 10
Hydrofacies: the Na-Cl type dominate; highly variable composition, including the Na-HCO 3 type. Greece The data indicate that a number of different processes are active in the different environments, and that gases (in particular CO 2 ) have an active role in driving equilibria Ca-HCO 3 type SO 4 reduction CaCO 3 deposition CaCO 3 solution seawater 11
Hydrofacies: waters range from Ca-HCO 3, Na-Cl and the Na-HCO 3 type Macedonia 12
Most waters range from the Ca-HCO 3 to the Na-HCO 3 type Serbia Ca-HCO 3 type 13
Hydrofacies: mostly waters range from the Ca- HCO 3 and Na-HCO 3 type Na-Cl waters are limited Slovenia Conservative mixing between Ca-HCO 3 and Na-Cl waters are absent or negligible, indicating non-interacting reservoirs. Possible processes of freshwater intrusion are highlighted Ca-HCO 3 type seawater fresh water intrusion? 14
dd Waters origin: oxygen and hydrogen stable isotopes 10 Seawater source -12-10 -8-6 -4-2 0 2-10 -30-50 Serbia Greece Bosnia Slovenia d 18 O -70-90 Global Meteoric Water Line Central Mediterranean Water Line Most water-types have a meteoric origin. Na-Cl waters have a marine origin. Modern or ancient seawater? 15
Strontium isotopes: the modern or ancient origin of the marine Na-Cl type waters 0.7095 The 87 Sr/ 86 Sr ratio in seawater changes through geological time, allowing the age of the saline reservoir to be inferred 87 Sr / 86 Sr 0.7090 0.7085 0.7080 0.7075 0.7070 0.7065 0 100 200 300 400 500 Numerical Age, Ma Mc Arthur et al, 2001, 2007 16
HCO3 (mg/l) The Na-HCO 3 waters: extreme compositions through water-rock interaction 10000 1000 100 10 Bosnia Croatia Greece Macedonia 1 0.1 0.1 1 10 100 1000 10000 100000 Na (mg/l) Serbia Slovenia The waters from the Peri-Adriatic region show a Na vs. HCO 3 scattered distribution, in some cases with high HCO 3. Most waters from Greece are displaced towards low HCO 3 17
Some of the Na HCO 3 thermal waters likely reflect silicate weathering: silicates have low solubilities (NaAlSi 3 O 8 : 6x10-7 mol/l) and low dissolution rates the Na-HCO 3 signature reflects long residence times of waters in deep environments These thermal reservoirs are of particular interest 18
TDS (mg/l) Cl (mg/l) Temperature 3500 3000 (Miosic et al., 2013) 2500 2000 1500 1000 500 0 0 20 40 60 80 Temperature ( C) 400 350 300 250 200 150 100 50 0 0 20 40 60 80 Temperature ( C) Bosnia the trends of increasing temperature with increasing salinity and chloride content suggest, at least in some cases, the role of a high-thermal component of marine origin 19
Ca (mg/l) SiO2 (mg/l) Temperature 70 60 50 40 30 (Milenic et al., 2012) Serbia 20 10 0 0 20 40 60 80 Temperature ( C) 140 120 100 80 60 40 20 0 0 10 20 30 40 Temperatura ( C) Despite some scatter, the trend of increasing temperature with increasing silica likely indicate the role of high-thermal components related to silicate dissolution. The pattern of decreasing calcium with increasing temperature could reflect the saturation state of carbonates, yielding cement precipitation 20
SiO2 (mg/l) Temperature 100 90 80 70 60 50 40 30 20 10 0 (Lambrakis & Stamatis, 2008) 0 20 40 60 80 100 Temperature ( C) Greece Also in this case, some of the waters likely reflect the role of high-thermal, deep components related to silicate dissolution. Trends of decreasing calcite and dolomite saturation states with decreasing temperature are also observed, indicating the role of CO 2 -mediate carbonate equilibria 21
A variety of hydrofacies characterize the active geothermal manifestations used for thermal balneology in the Peri-Adriatic region Waters have both a meteoric origin or maintain a marine signature, likely attributable to modern or to remnants of ancient seawater. Isotopic studies would help to clarify this point Some of the waters are diagenetically modified by long-lasting interactions with the aquifer lithologies, suggesting high residence times and/or slow flowpaths. In particular, the interactions with low-solubility silicate minerals is likely responsible for the particular Na-HCO 3 signature of the thermal, deeper reservoirs. These waters would deserve further studies. For these waters, long-term abstraction should be carefully planned In some cases, mixing of waters with different geochemical signature occurs 22
Example of application of geochemical tools The Monfalcone (Grado Lignano) thermal waters in the coastal area of the Friuli Plain (NE-Italy): a case history 23
Sites of study Monfalcone: natural springs (34-40 C) Grado: Grado-1 well (1100 m depth; 43 C) Lignano: SIL well ( 2000 m depth; 53 C) 24
Geological framework of the Trieste Gulf Lignano Monfalcone B Grado A Mesozoic-Cenozoic carbonate Units Paleogene Terrigenous Units Quaternary alluvial sediments After Carulli (2011) J. Geodyn. 51, 156-165 25
Multichannel seismic profile across the Trieste Gulf A Lignano Grado-1 Monfalcone B Plio-Quaternary marine sediments Eocene Flysch terrigenous sequence Friuli Carbonate Platform Busetti et al 2010 26
Hydrochemistry The Monfalcone, Grado-1 and Lignano deep waters have all the chemical signature like modern seawater NOTE Some of the waters from geothermal wells in the FGV Low Plain are of the Na-HCO 3 type Lignano Monfalcone Grado-1 Wells in the FVG Low Plain freshwater Modern seawater 27
Na (mg/l) The marine signature 12000 10000 Modern seawater 8000 6000 4000 Monfalcone Grado-1 Lignano 2000 0 0 5000 10000 15000 20000 Cl (mg/l) Waters from the Lignano well are the most representative of the marine end-member Dilution processes are active at the Monfalcone springs and Grado-1 well to a lower extent A marine component also contributes to waters from geothermal wells in the FVG Plain (crosses) 28
dd Dilution processes and origin of non-marine component 10 0-13 -11-9 -7-5 -3-1 1-10 -20-30 -40 Monfalcone Grado-1 Lignano Karst-type waters -50-60 FVG Low Plain geothermal wells d 18 O -70-80 The freshwater component might be represented by karst-type groundwaters 29
SO 4 (mg/l) Evolutionary patterns of the marine component 3000 2500 seawater 2000 Sulfate excess 1500 1000 Sulfate depletion Monfalcone Grado-1 Lignano 500 0 Sulfate excess and mixing corrosion 0 5000 10000 15000 20000 H H H 2 2 2 S 2O SO 4 2 S 2O Cl (mg/l) 2CaCO 2 H 2 SO 3 4 2CaCO 2Ca 3 2 2Ca SO 2 2 4 SO 2HCO 2 4 3 2HCO 3 30
Ca (mg/l) 1600 1400 1200 1000 Calcium excess 800 600 Monfalcone Grado-1 Lignano 400 seawater 200 0 0 5000 10000 15000 20000 Cl (mg/l) Diagenetic reactions of seawater with the hosting carbonates 31
Crustal gas-water interactions Samples deviate from the composition of air saturated water towards the admixing of crustal gases. This has implication on the volume of the thermal reservoir, since the original atmospheric signature has been replaced by crustal components: Sulfur springs in Eastern Alps Alpi orientali Monfalcone Grado-1 Lignano LOW WATER VOLUME OR DISCONTINUOUS WATER BODIES (?) 32
Modern or ancient seawater? 0,7084 Sr-isotope ratio of modern seawater 0.70918 0,7083 0,7082 87 Sr/ 86 Sr 0,7081 0,7080 0,7079 0,7078 0,7077 lignano monfalcone Grado-1 The thermal waters do not have the isotopic signature of present-day seawater! A possible common marine source for Monfalcone, Grado-1 and Lignano 33
87Sr / 86 Sr Deconvolving the diagenetic chemical changes and age of the thermal component Mass-balance equation 87 Sr Sr 86 87 86 excess 87 / Srreservoir Sr/ Srmeasured 1 Sr reservoir Sr/ 86 Sr excess Sr Sr excess reservoir 0,7095 0,7090 Monfalcone, Grado-1, Lignano Sr isotopic changes in seawater through time 0,7085 0,7080 0,7075 0,7070 0 20 40 60 80 100 Numerical Age, Ma the marine reservoir would be represented by ancient seawaters dated at Miocene 34
EC (us/cm) Temperature 45000 40000 35000 30000 25000 20000 15000 Monfalcone Grado-1 Lignano 10000 5000 Karst waters? 0 0 10 20 30 40 50 60 T ( C) The thermal reservoir is best represented by the saline ancient seawater Cooling by dilution is observed; however the karst-type waters of superficial aquifers are unsuitable as cold term 35
Deconvolution of the diagenetic processes of ancient seawater in the deep carbonate reservoir to estimate temperature of the geothermal component 60 70 C 80 Lignano Grado-1 65 C is the T at about 2.0-2.5 km in carbonate reservoir Modern seawater 36
Concluding remarks Friuli thermal waters represent the ancient, diagenetically modified seawater of a low temperature geothermal system Deep saline reservoir might represent remnants of seawater entrapped in Mesozoic carbonates during the late Oligocene Miocene sea transgression In some parts the deep thermal reservoir mixes and homogeneizes with colder kast(?)-type waters flowing in, through a deep and long-lasting circuit Thermal reservoir might be composed by discrete aquifers, or might not be widely extended Hydrothermal diagenesis including dolomite cementation might change the water flow dynamics at depth 37
Generalized geothermal model? 38
Thank you for your attention 39
Some of the waters show a Na vs. Cl linear correlation, suggesting the role of seawater or dissolution of marine salts Cl (mg/l) 100000 10000 1000 100 10 1 0.1 0.1 1 10 100 1000 10000 100000 Na (mg/l) On this ground, the contribution of the marine component can be removed 40
HCO3 (mg/l) 9000 8000 7000 6000 5000 4000 3000 2000 1000 0 0 500 1000 1500 2000 2500 3000 Na (mg/l) After correction a more clear correlation in observed, following the equimolar distribution of Na and HCO 3 ions (solid line) This is consistent with the (incongruent) silicate dissolution to generate the Na- HCO 3 waters: 2NaAlSi SiO as an example 3O8 3H2O 2CO2 2Na 2HCO3 Al2Si2O5 ( OH) 4 4 2 41