 Geochemical characteristics of thermal waters of Hrvatsko zagorje  Tamara Marković, Staša Borović and Ozren Larva Croatian Geological Survey, Department of hydrogeology and engineering geology, Sachsova, 2, Zagreb, Croatia; (tmarkovic@hgi-cgs.hr; sborovic@hgi-cgs.hr; olarva@hgi-cgs.hr) doi: 10.4154/gc.2015.05 ABSTRACT In the region of Hrvatsko zagorje there are a lot of spa localities. The study area is in the northern Pannonian part of Croatia which is characterized by a high geothermal gradient (0.049 °C/m) and surface heat flow (76 mW/m2). Al- though there are a lot of thermal occurrences in the study area, only ten were taken into consideration because there was a lack of geochemical data at other locations. The thermal springs considered are: Harina Zlaka, Krapinske top- lice, Tuheljske toplice, Stubičke toplice, Sutinske toplice, Šemničke toplice, Topličica (Mađarevo), Topličica (Gota- lovec), Podevčevo, Vraždinske toplice. A compilation of geochemical data from different sources together with our own measurements has been used in this study. The aim of this paper is to review the thermal waters’ geochemical characteristics and demonstrate how these features can be used to discern their origin and the aquifer they equilibrat- ed with. The geochemical characteristics of thermal waters of the study area suggest that the water is in equilibrium with dolomite which means that dolomite is the thermal aquifer in the study area. Keywords: thermal water, groundwater temperature, chemical composition, water type, Hrvatsko zagorje 1. INTRODUCTION Hot water occurring at the earth’s surface has been a subject of awe since the dawn of humankind. Ancient civilizations revered thermal springs because they were believed to have supernatural and healing powers (ROUTH et al., 1996; LAMOREAUX & TANNER, 2001). Thermal waters are in- creasingly being used in the world for industrial processing, agriculture production, heating, electricity production and the extraction of rare elements (ATKINSON & DAVIDSON, 2002; HELLMAN & RAMSEY, 2004; PETRACCIA et al., 2005; LUND, 2006; HAEHNLEIN et al., 2010). Moreover, with the increasing popularity of spas and the growing im- portance attached to the ‘natural’ health industry, thermal waters are becoming a valuable asset for centres of hydro- therapy (ROUTH et al., 1996; LUND, 1996; MAROVIĆ et al., 1996; HARVEY, 2007; BITUH et al., 2009). Thermal waters are thus a natural resource that could make consider- able contribution to the local and regional economy. How- ever, care should be exercised when making decisions re- garding the appropriate use of thermal waters. Croatia is relatively richly endowed with thermal waters and their uti- lization has been a tradition for millennia. Hrvatsko zagorje is a region situated in the north-west- ern part of central Croatia (Fig. 1). It stretches between the Croatian-Slovenian borderline and the Medvednica and Kalnik Mountains. Most of the region is characterized by undulating relief, especially the foothills of Ivanščica Mt., populated since antiquity. A “corrugated relief” with multi- ple interconnected mountains and valleys has always been an obstacle for the construction of meaningful transport net- works in the area. The region sees a daily workforce migra- tion to the national capital Zagreb, but it is also a traditional recreational area for many citizens of Zagreb, who visit its numerous spa resorts (ROGLIĆ, 1974). The number of spa localities in Hrvatsko zagorje is high, and they are densely distributed across the region. The sim- ple indicator of this situation is the sheer number of topo- nyms „toplica“ (singular) or „toplice“ (plural), meaning ther- Geologia Croatica 68/1 67–77 7 Figs. 2 Tabs. Zagreb 2015 Geologia CroaticaGeologia Croatica Geologia Croatica 68/1Geologia Croatica 68 mal water spring(s), and „topličica“ meaning little thermal water spring (with little referring to the lower temperature of the water, not its quantity). Some of the localities were utilized even in prehistoric time, and in ancient Roman time they were already established as curative destinations, espe- cially Aquae Iassae (Varaždinske Toplice) with it’s nume- rous upright archaeological remains and Aquae Vitae (Kra- pinske Toplice). Even the Roman name of the Pannonian tribe populating this area (lat. Iassi) is derived from the sup- posedly curative properties of thermal waters because the root ias- or iat-came from the ancient Greek language mean- ing to heal or cure (SCHEjBAL, 2003). The temperatures of the thermal waters in Hrvatsko za- gorje vary from 16 to 65 °C. There are many available scales to classify thermal resources and thermal springs. According to all of the geothermal resource classifications (e.g. AXELS- SON & GUNLAUGSSON, 1978; MUFFLER & CATALDI, 1978; UNGEMACH & ECONOMIDES, 1987; BENDERIT- TER & CORMY, 1990; HOCHSTEIN, 1990; NICHOLSON, 1993; SANYAL, 2005) these resources at Hrvatsko zagorje are considered low enthalpy resources. Due to their relatively low temperatures, the classifications for thermal springs are more appropriate. These classifications vary significantly, because they use different benchmark values for identifying cold, warm and hot springs. It is common to classify the springs as warm if their temperature is lower than that of the human body, and hot if it is higher (e.g. DEMING, 2002; KREŠIĆ, 2010). Some authors have also considered the lower boundary for classifying waters as thermal, taking e.g. 20 °C which is a limit for balneological purposes (HARA- MUSTEK et al., 1952; IVEKOVIĆ & PEROŠ, 1981) or the temperature which disables ice formation at the rim of the spring, around 15 °C (ĐUROVIĆ, 1960). From the hydro- geological point of view however, all groundwaters with temperatures higher than the average rainwater temperature (approximately the mean annual temperature of the locality) are considered thermal, although they cannot be used in bal- neology. Since the climate of the investigated region is Cfwbx’’ according to the Köppen classification (GAjIĆ-ČAPKA & ZANINOVIĆ, 2008) – a temperate (C) humid (f) climate with no extremely dry periods throughout the year, and the annual mean temperature of the region varies between 7–8 °C at the highest mountains and 10–11 °C in the valleys (GAjIĆ- ČAPKA & ZANINOVIĆ, 2008) – in a hydrogeological sense all groundwaters with temperatures higher than this average are considered thermal. In this paper the waters have been categorized according to the modified balneological classi- fication (Tab. 1). It was created on the basis of the traditional Croatian balneological classification (HARAMUSTEK et al., 1952), which actually stemmed from much earlier con- siderations of German and Swiss medical balneology experts (jACOBj, 1907; HINTZ & GRÜNHUT, 1916; HART- MANN, 1925). The modification was applied to the lower part of the scale, thus including the waters where the tem- perature is higher than the annual mean, but lower than that required in balneology. Figure 1: The location of the Hrvatsko zagorje region in Croatia. Table 1: Categorization of thermal localities in Croatia on the basis of water temperature (adopted from KOVAČIĆ & PERICA, 1998). Category SUBTHERMAL HYPOTHERMAL HOMEOTHERMAL HYPERTHERMAL Temperature 13–20 °C >20–34 °C >34–38 °C >38 °C Marković et al.: Geochemical characteristics of thermal waters of Hrvatsko zagorje Geologia Croatica 69 The waters considered here range from subthermal to hyperthermal. The modes of utilization vary according to the temperatures, e.g. waters of the lowest temperatures are used for the public water supply, bottling and fish farming, while waters of the highest temperatures (40–65 °C) are utilized for water and space heating. According to its geothermal characteristics, Croatia can be divided into two parts where the mean values of the im- portant geothermal parameters vary significantly. The north- ern Pannonian part is characterized by a high geothermal gradient (0,049 °C/m) and surface heat flow (76 mW/m2), while the southern Dinaric part has a low average geother- mal gradient (0,018 °C/m) and surface heat flow (29 mW/ m2) (BOŠNjAK et al., 1998). In accordance with these char- acteristics, many locations where thermal springs exist are located in the northern part of Croatia where the groundwa- ter receives enough heat from its surroundings. Attempts to explain the appearance and origin of the thermal springs in Croatia date back as far as the 18th cen- tury. Many geologists at those times were studying the warm springs of Hrvatsko zagorje owing to their unusual density in the region, and they postulated a number of different hy- potheses to explain their occurrence. For example, in the be- ginning of 20th century Dragutin Gorjanović-Kramberger was convinced that the thermal waters of north-western Croatia were of volcanic origin with magma chambers heat- ing the waters and deep faults, so called “thermal lines”, bringing them to the surface. The main shortcoming of this hypothesis is that there are no recent magmatic bodies to supply heat, but because of Gorjanović-Kramberger’s au- thority it still has a fair number of supporters (ŠIMUNIĆ, 2008). The idea that thermal waters could be of meteoric origin and then heated during circulation through the sub- surface first appeared in the works of Gjuro Pilar at the end of 19th century and is still supported by all the evidence. Un- til now, it has been proven that all significant thermal springs in the Pannonian part of Croatia are located in the tops of anticlines which were fractured by faults of SW – NE or N – S strike. The intersections of faults and folds usually represent a strongly fractured environment of high perme- ability which enables the upwelling of heated water at depth to the surface (ŠIMUNIĆ, 2008). The geochemical studies of HORVATINČIĆ et al. (1991; 1996) proved that the ther- mal waters of the Hrvatsko zagorje area are of meteoric ori- gin. In other papers, the geochemical characteristics of some springs in the research area have been studied and all data point to a meteoric origin (MIHOLIĆ, 1940a, 1952, 1959; MAROVIĆ et al., 1996; MARKOVIĆ & KOVAČIĆ, 2006; BITUH et al., 2009; POLANČEC, 2011). The aim of this paper is to review the thermal waters’ geochemical characteristics and demonstrate how geochem- ical features of thermal waters can be used to discern their origin and the aquifer they equilibrated with. Although there are lots of thermal occurrences in the study area, only ten were taken into consideration because there was a lack of geochemical data at other locations. The thermal springs con- sidered are: Harina Zlaka, Krapinske toplice, Tuheljske top- lice, Stubičke toplice, Sutinske toplice, Šemničke toplice, Topličica (Mađarevo), Topličica (Gotalovec), Podevčevo and Vraždinske toplice (Fig. 2). 2. HYDROGEOLOGICAL SETTINGS In Harina Zlaka (Fig. 2) there was a natural thermal spring with a water temperature of 32.8 °C (MIHOLIĆ, 1940a). The spring dried up after 1966, following exploitation of thermal water from the well on the opposite bank of the Sutla River for the Atomske toplice spa in Podčetrtek, Slovenia. The Figure 2: Locations of the thermal springs described in the paper. Geologia Croatica 68/1Geologia Croatica 70 well is drilled in the carbonate sediments (dolomites). In 1998, near the dried up spring, the well HZL-1 was drilled, also in dolomites. It has been tested with a yield of 45 L/s, and the water temperature at the wellhead was 19–20 °C (MRAZ & KRSNIK, 2005). Thermal springs of Krapinske toplice (Fig. 2) occur in a narrow valley of the Topličica stream and there are three main springs and a few springs of lower yield. The main oc- currence of thermal springs is in the area of Pučke and ja- kobove kupelji. The total measured yield of all springs in Krapinske toplice, according to BAĆ & HERAK (1962), ranges from 69 to 81 L/s. The majority of springs occur along the boundary between the dolomite and limestone. In 1985 a deep well was drilled about 250 m north from the thermal springs. The well passed through limestone and dolomite before it ended up in calcarenites and shale at a depth of 861 m (ŠIMUNIĆ, 1986). The well was tested with a yield of 30 L/s and the water temperature was 45 °C. Thermal wa- ters in Krapinske Toplice are used for recreation, balneo- therapy, space heating and sanitary water. Tuheljske toplice springs (Fig. 2) are located nearby at the Krapinske toplice thermal water source, 4.5 km to the southeast. The surrounding area of Tuheljske toplice consists of limestone, dolomite, sandstone, marl, clay, gravel and sand. Thermal springs are ascending springs and they occur in several places, from which the most important are the Dadino vrelo and Vrelo u bari springs. The total yield of the thermal springs was 85 L/s (BAĆ & HERAK, 1962) and the temperature varies from 32.5 to 33.1 °C. These thermal wa- ters have so far only been used for recreational purposes. Stubičke toplice springs (Fig. 2) are located on the north- ern side of Medvednica Mt. They had consisted of two ma- jor and several smaller springs, but they dried up when deep wells were drilled and put into operation. Spring water tem- perature varied from 30 to 49.8 °C and the water in the wells has a temperature of 65 °C. The yield of the springs was 18 L/s (BAĆ & HERAK, 1962). The surrounding area of the Stubičke toplice springs comprises clastic carbonate sedi- ments, dolomite, alluvial, proluvial-deluvial deposits and sil- iciclastic rocks. Today thermal water from the wells is used for recreation, balneotherapy, water and space heating, san- itary water and greenhouse heating, and it is the locality with most diversified utilization of geothermal waters in Croatia (BOROVIĆ & MARKOVIĆ, 2015). Thermal springs of Sutinske toplice (Fig. 2) are located in the canyon of the Sutinska stream. Water temperature ranges from 30 to 37.4 °C and the water emerges from frac- tured dolomites. In 1988 a well was drilled in the southern exit of Sutinska canyon and a dolomite aquifer was detected in the interval between 75 – 385 m. A pumping test was car- ried out resulting in a yield of 112 L/s and the water tem- perature at the wellhead varied from 38 to 39 °C (BRITVIĆ, 1988). The water had been used for a swimming pool during the summer, but it was closed in the summer of 2013. The Šemničke toplice springs (Fig. 2) are located near to the Krapinske toplice thermal water source, 9 km to the northwest. The yield of the thermal spring is 0.5-1.0 L/s and the water temperature ranges from 31 to 33.1 °C. The well was drilled in the vicinity of the springs and carbonates were found. The well was pumped with the yield of 20 L/s, and the water temperature varied from 33 to 36oC (CAPAR, 1982). The thermal water of Šemničke toplice is not currently being utilized. The Topličica springs in Mađarevo (Fig. 2) occur at the foot of the steep slopes of the Ivanščica massif near Novi Marof. There are four major springs and several smaller ones. The water temperature is inconsistent: during the sum- mer it varies from 18 to 22.5 °C, and during winter temper- atures are lower. Significant fluctuations in temperature and yield indicate mixing of cold and thermal water. The total yield of all the springs has not been measured. Thermal wa- ter is used for heating a fish farm nearby and for recreation in the summer. The Topličica spring in Gotalovec (Fig. 2) is located on the southern slopes of the central part of the Ivanščica mas- sif. The yield of the spring ranges from 15 to 18 L/s and the water temperature is from 25 to 26 °C. The spring occurs at the contact between the fractured dolomite aquifer and erup- tives and clastic sediments. The water is being used for bot- tling as natural spring water. In Podevčevo (Fig. 2) there are two thermal springs, with variable temperatures, from 16.3 to 20.0 °C depending on the hydrological cycle. The aquifer is fractured dolomite. The total yield of both springs does not exceed 1 L/s. The springs are not used. Varaždinske toplice springs (Fig. 2) are located in the north-eastern corner of the study area, nearby the town of Varaždin. They are among the best known, are the largest spas in Croatia and have been used for the longest time. The water temperature varies from 56.5 to 57.5 °C and the yield ranges from 45 to 50 L/s. The thermal waters occur in the brecciated dolomites which are covered by clastic sediments (ŠIMUNIĆ, 2008). The largest spring in Varažidnske toplice is Klokot, accompanied by three smaller hypothermal springs with temperatures of 24–25 °C. There are also a lot of wells. Thermal waters from the spring and wells are used for rec- reation, balneotherapy, water and space heating (BOROVIĆ & MARKOVIĆ, 2015). 3. MATERIALS AND METHODS A compilation of geochemical data from different sources has been used in this study, including chemical analyses by MIHOLIĆ (1952) and jURŠIĆ-MITROVIĆ (2001), and chemical analyses performed during research for the Basic Hydrogeological Map sheets Rogatec and Varaždin in Octo- ber 2009. From a total of 32 analyses, 10 analyses were per- formed by MIHOLIĆ (1952), 10 by jURŠIĆ-MITROVIĆ (2001) and 12 analyses were performed in the scope of the research for the Basic Hydrogeological Map. The samples taken during research for hydrogeological map were anal- ysed at the Hydrochemical laboratory of the Department of Hydrogeology and Engineering Geology of Croatian Geo- logical Survey. Prior to sampling, the electrical conductivity (EC), pH and temperature (T) were determined using por- table WTW probes. Total alkalinity was measured by titra- Marković et al.: Geochemical characteristics of thermal waters of Hrvatsko zagorje Geologia Croatica 71 Table 2: Mean, minimum and maximum values of physicochemical parameters in observed thermal waters. Harina Zlaka Krapinske toplice Tuheljske topliice Stubičke toplice Sutinske toplice Šemničke toplice Topličica Mađarevo Topličica Gotalovec Podevčevo Varaždinske toplice EC(mS/cm) mean min max 635 597 655 519 438 566 623 620 627 627 610 653 515 510 520 532 522 545 568 539 596 422 403 440 675 672 678 1206 1176 1235 T(°C) mean min max 24.2 19.8 32.8 41.8 40.2 45.0 32.8 32.5 33.1 53.9 47.0 65.0 34.9 34.0 35.7 32.7 32.5 33.1 23.0 22.5 23.4 25.6 25.5 25.7 17.7 16.0 19.4 54.2 51.2 57.2 pH mean min max 7.4 7.4 7.5 7.1 7.0 7.3 7.3 7.3 7.4 7.0 6.7 7.3 7.2 6.9 7.5 7.4 7.4 7.5 7.3 7.3 7.3 7.3 7.2 7.3 7.8 7.6 8.0 6.7 6.6 6.8 Na+(mg/l) mean min max 7.2 2.8 14.0 11.4 9.1 17.0 11.3 10.9 11.8 25.3 20.8 28.5 5.0 3.8 6.2 9.7 6.3 11.8 3.2 2.4 4.0 4.3 3.4 5.2 26.1 22.6 29.6 86.2 76.8 95.6 K+(mg/l) mean min max 2.3 2.0 2.9 3.3 2.8 3.5 3.0 2.9 3.1 6.3 4.6 7.3 1.7 1.6 1.8 3.0 2.9 3.1 1.6 0.9 2.2 1.9 1.0 2.8 9.4 9.3 9.6 32.1 23.2 40.9 Ca2+(mg/l) mean min max 68.7 61.9 75.3 54.9 52.4 58.2 63.8 63.2 64.6 66.9 61.0 71.7 65.2 53.7 76.6 61.5 57.6 63.7 62.3 53.7 70.9 60.9 55.8 65.9 58.8 56.6 60.9 126.2 125.1 127.3 Mg2+(mg/l) mean min max 36.0 31.8 44.1 30.7 27.8 35.1 36.9 34.9 38.0 25.4 23.1 30.2 29.8 26.6 33.0 39.0 37.9 41.0 30.4 26.8 33.9 24.7 21.6 27.8 32.8 29.6 36.0 28.3 27.3 29.2 HCO3 – (mg/l) mean min max 362.0 333.8 407.5 296.3 280.6 330.9 367.9 364.8 371.0 281.8 262.3 300.4 288.1 288.0 288.2 355.6 330.9 371.0 328.2 298.0 358.3 292.1 288.2 296.0 341.3 338.0 344.7 454.4 445.3 463.5 SO4 2– (mg/l) mean min max 25.7 15.6 39.7 37.3 35.3 39.7 37.6 35.2 39.4 88.3 80.2 97.9 31.3 27.9 34.8 40.6 38.2 44.2 20.6 15.8 25.4 14.7 13.1 16.3 71.0 61.9 80.0 168.6 155.8 181.5 Cl–(mg/l) mean min max 7.7 3.4 10.7 3.1 2.5 3.7 3.5 2.6 5.0 12.1 9.6 13.9 3.9 1.5 6.2 3.6 2.6 5.3 3.4 3.1 3.7 6.1 5.9 6.4 14.3 12.5 16.2 82.8 79.0 86.6 Number of samples 3 5 3 4 3 3 3 3 2 3 tion with 1.6 N H2SO4 using phenolphthalein and bromcre- sol green – methyl red as indicators. Anions (SO4 2- and Cl-) were measured by ion chromatograph (LabAlliance Ltd) and dissolved cations (Na+, K+, Mg2+ and Ca2+) by atomic ab- sorption spectrophotometer (Perkin Elmer Analyst 700). The NETPATH-WIN program (EL-KADI et al., 2010) was used to calculate saturation indices of the observed wa- ters with respect to the calcite and dolomite mineral phases and CO2 partial pressure (pCO2). 4. GEOCHEMICAL CHARASTERISTICS In Table 2 the minimum, maximum and mean values of ob- served geochemical parameters in thermal waters of ten in- vestigated water sources are given. The mean water temperatures ranged from 19.2 to 54.2 °C. The highest value is observed in Varaždinske top- lice and the lowest in the water of Podevčevo spring (Tab. 2). According to the temperature categorization by KOVAČIĆ & PERICA (1998), the spring water in Podevčevo is sub- thermal; water from Sutinske toplice is homeothermal; wa- ters from Harina Zlaka, Tuheljske toplice, Šemničke toplice, Topličica (Mađarevo) and Topličica (Gotalovec) are hypo- thermal and waters from Stubičke toplice, Krapinske toplice and Varaždinske toplice are hyperthermal. The mean pH values ranged from 6.7 to 7.8 (Tab. 2). The lowest value is observed in the waters from Varaždinske top- lice (acid) and the highest is observed in the water of Podevčevo spring (alkaline). In Stubičke and Krapinske toplice lower pH values in waters are also observed. The pH is controlled by dissolution of gases (such as CO2, H2S, NH3, CH4) which are entering the aquifer system and geochemical processes such as the synthesis and mineralisation of biomass, redox proc- esses etc. in the aquifer. The higher the amount of dissolved gases in the aquifer system, the lower the pH value will be. From the water quality monitoring of the Special Hospital for Medical Rehabilitation Varaždinske toplice (SPECIAL HOS- PITAL, 2014) a high amount of H2S (10.4 mg/l) in thermal waters has been observed. Furthermore, in the waters from Stubičke toplice a high content of CO2 is observed - 101.49 mg/l (HOSPITAL STUBIČKE TOPLICE, 2014) and H2S is also detected in the water from Krapisnke toplice (HOSPITAL KRAPINSKE TOPLICE, 2014). Also, the contribution of H+ ions during the geochemical processes in the aquifer should not be dismissed. For example during the redox process H+ ions occur, e.g. FeS2 + 3.75 O2 + 1.5 H2O= Fe(OH)3 + 2SO4 2- + 4H+ (pyrite oxidation). Geologia Croatica 68/1Geologia Croatica 72 The highest EC value is observed at the Varaždinske toplice springs, with a mean value of 1205 mS/cm, while the lowest value is observed at the Topličica (Gotalovec) spring, with a mean 422 mS/cm (Fig. 3). It is visible from Fig. 3 that some of the waters have similar EC values, e.g. waters from Harina Zlaka, Tuheljske and Stubičke toplice. In the next group there are waters from Sutinske, Krapinske and Šemničke toplice, with Topličica (Mađarevo) close to them. Podevčevo spring waters lie between the last group and the highest EC value. The EC value is influenced by the amount of dissolved solids in the water. If the water contains a high amount of dissolved solids, the EC values will be higher, and vice versa. The amount of dissolved solids in the water is controlled by numerous factors. The type of rock mass which constitutes the catchment area of the spring is one of the im- portant factors. If a rock mass of high solubility (e.g. lime- stone, volcanic ash) constitutes the catchment area, then the EC will be higher. Waters which contain high concentrations Figure 3: Distribution of mean EC values. Figure 4: Piper diagram of major constituents in thermal waters. Marković et al.: Geochemical characteristics of thermal waters of Hrvatsko zagorje Geologia Croatica 73 of CO2 and H2S gases, resulting in high concentrations of dissolved CO2 and H2S, facilitate the dissolution of carbon- ate and silicate rocks (YOSHIMURA et al., 2004) and con- sequently, such waters have higher EC values. Moreover, higher water temperature and a more fractured aquifer in combination with low pH values of the water also facilitate dissolution of carbonate rocks and sediments (BERTRAM et al., 1991; GAUTELIER et al., 1999; MORSE & ARVID- SON, 2002) resulting in high EC values. The thermal water source of Varaždinske toplice has the highest concentrations of sodium, potassium, sulphate and chloride (Tab. 2). It is followed by the thermal springs of Stubičke toplice, the spring in Podevčevo and the thermal water source of Krapinske toplice (Tab. 2). The lowest val- ues are observed in the springs of Topličica in Gotalovec and Mađarevo. Higher concentrations of sodium, potassium, sul- phate and chloride ions are the consequence of the dissolu- tion of magmatic rocks and siliciclastic (including pyroclas- tic) sediments in the areas surrounding the thermal water sources. While some 55 ionic species have been reported in ash leachates (WITHAM et al., 2005) the most abundant species usually found are the cations Na+, K+, Ca2+ and Mg2+ Figure: 5 a) Molar ratios of magnesium and calcium ions in thermal waters; b) Bivariate mixing diagram of Na+-normal- ized Ca2+ versus Na+-normal- ized Mg2+. a) b) Geologia Croatica 68/1Geologia Croatica 74 and the anions Cl-, F- and SO4 2- (jONES & GISALSON, 2008; WITHAM et al., 2005). According to the major ionic composition, water from Varaždinske toplice belongs to a NaCaMg-HCO3SO4 mixed type; Podevčevo and Stubičke toplice waters belong to a CaMgNa-HCO3SO4 mixed type; while the rest of thermal waters (Harina Zlaka, Tuheljske toplice, Krapinske toplice, Topličica (Mađarevo), Topličica (Gotalovec), Sutinske top- lice and Šemničke toplice) belong to the CaMg-HCO3 type (Fig. 4). Hydrochemical facies, as recognised from water chem- istry data, are a consequence of the chemistry of the recharg- ing water, and water–aquifer matrix interactions (e.g. cation exchange), as well as groundwater residence time within the aquifer. Two general processes contribute solutes to ground- water: carbonate dissolution and silicate weathering (GAR- RELS & MACKENZIE, 1971). The chemistry of the evolv- ing water depends not only on the bulk chemistry of the matrix, but also on the rate of weathering. MEYBECK (1987) commented that weathering rates of limestone and dolomite are up to 80 times and 12 times, respectively, faster than silicate weathering rates. From Fig. 5a it can be ob- served that the groundwater chemistry is primarily control- led by the dissolution of carbonate rocks (dolomite and lime- Figure: 6 Saturation index and logpCO2 calculations a) in waters of Harina Zlaka, Tuheljske toplice, Krapinske toplice, Šemničke toplice and Topličica Mađarevci b) in water of Topličica (Gota- lovec), Podevčevo spring, Varaždinske toplice, Sutinske toplice and Stubičke toplice a) b) Marković et al.: Geochemical characteristics of thermal waters of Hrvatsko zagorje Geologia Croatica 75 stone). In addition, the saturation indices of the observed thermal waters with respect to calcite and dolomite mineral phases indicate the dissolution of carbonate rocks (Figs. 6 a, b). Generally, it can be stated that the waters of all observed thermal sources were saturated with respect to calcite and dolomite (Figs. 6 a, b). It was observed that the saturation index with respect to dolomite was higher than with respect to calcite. The chemistry of waters from Varaždinske toplice, Podevčevo and Stubičke toplice is also influenced by the dis- solution of carbonates. However, on the bivariate mixing diagrams (Fig. 5b) of Na+-normalized Ca2+ versus Na+-nor- malized Mg2+, the samples from Varaždinske toplice, Podevčevo and Stubičke toplice tend to fall within the sili- cate weathering domain. This suggests that incongruent leaching of the siliceous volcanic rocks is present in the aq- uifer. The diagram (Fig. 5b) also shows that carbonate dis- solution is dominant in waters of the thermal sources of the Harina Zlaka, Krapinske toplice, Tuheljske toplice, Sutinske toplice, Šemničke toplice, Topličica (Mađarevo) and Topličica (Gotalovec). Also, partial pressure of CO2 (pCO2) of observed waters were calculated. It was observed that the pCO2 values in all sampled waters, were in disequilibrium with respect to the atmosphere (10-3.5 or logpCO2= -3.5). The highest partial pressure was observed in the thermal waters of Varaždinske toplice and Krapinske toplice (Figs. 6 a, b) as a consequence of gas influence. The influence of cation exchange was evaluated by an equivalent bivariate plot of corrected bivalent cations versus corrected Na+ (Fig. 7). Concentrations of bivalent cations (Ca2+ and Mg2+) that may have been involved in exchange reactions were corrected by subtracting equivalent concen- trations of associated anions (HCO3 - and SO4 2- ) that would be derived from other processes (e.g. carbonate or silicate weathering). Similarly, Na+ that may be derived from the aquifer matrix can be accounted for by assuming that Na+ contributions of meteoric origin would be balanced by equiv- alent concentrations of Cl- (MCLEAN & jANKOWSKI, 2000). For active cation exchange taking place in the aqui- fer, the slope of this bivariate plot should be -1 (i.e. y = -x) (jANKOWSKI et al., 1998). Due to small number of sam- ples, all were included in a single diagram. The slope of -1.35 for the observed samples indicates cation exchange (Fig. 7). The mixed cation exchange-mineral weathering might be a reflection of the longer residence time of thermal water in the subsurface. Although a bivariate plot (Fig. 7) shows influence of exchange processes and the bivariate mixing diagram (Fig. 5b) indicates silicate weathering, the majority of the analyzed samples tend to fall closer to the carbonate composition than to the silicate composition which means that dolomite is the dominant geothermal aquifer in the study area. 5. CONCLUSION In the region of Hrvatsko zagorje there are a lot of spa locali- ties. In this part of Croatia geothermal characteristics are fa- vourable. The water temperatures of thermal springs vary from 16 to 65 °C, i.e. from subthermal to hyperthermal. The highest water temperatures are observed in Varaždinske and Stubičke toplice water sources and the lowest in the Podevčevo spring. Lower pH values are observed in the waters from Varaždinske, Stubičke and Krapinske toplice and the highest is observed in the water of the Podevčevo spring. Low pH values are the consequence of H2S and CO2 gas dissolution. According to the major ionic composition, the waters belong to types ranging from mixed NaCaMg-HCO3SO4 to the CaMg-HCO3 type. The geochemical characteristics of the observed thermal waters indicate that the waters are in equilibrium with dolomite which means that dolomite is the Figure 7: Bivariate plot of (Ca2+ + Mg2+) - (HCO3 - + SO4 2-) against (Na+ - Cl-). Geologia Croatica 68/1Geologia Croatica 76 thermo-mineral and mineral waters – in Croatian and French].− In- stitute for geological and geophysical research, Belgrade, 69 p. EL-KADI A.I., PLUMMER L.N., & AGGARWA L.P. (2010): NET- PATH-WIN: An Interactive User Version of the Mass-Balance Model, NETPATH.− Ground Water, 49/4, 593−599. GAjIĆ-ČAPKA, M. & ZANINOVIĆ, K. (2008): KLIMA HRVATSKE [Climate of Croatia – in Croatian]. – In: ZANINOVIĆ, K. (ed.): Klimatski atlas Hrvatske 1961. −1990. i 1971.−2000. [Climate at- las of Croatia 1961−1990 and 1971−2000 – in Croatian].− Mete- orological and Hydrological Service of Croatia, Zagreb, 200 p. GARRELS, R.M. & MACKENZIE, F.T. (1971): Evolution of sedimen- tary rocks.− Norton, New York, 303 p. GAUTELIER, M., OELKERS, E.H. & SCHOTT, j. (1999): An experi- mental study of dolomite dissolution rates as a function of pH from 0.5 to 5 and temperature from 25 to 80 oC.− Chem. Geol., 157, 13–26. HAEHNLEIN, S., BAYER, P. & BLUM, P. (2010): International legal status of the use of shallow geothermal energy.− Renewable and Sustainable Energy Reviews, 14, 2611–2625. HARAMUSTEK, B., MIHOLIĆ, S. & TRAUNER, L. (1952): Godišnjak Baleološko-klimatološkog instituta NR Hrvatske [Annual report of the Balneology and Climatology Institute of PR Croatia – in Croatian].− Balneology and climatology institute, Zagreb, 133 p. HARTMANN, A. (1925): Die Mineral- und Heil-Quellen des Kantons Aargau: ihre geologischen und chemischen Verhältnisse [Mineral and curative springs of Aargau canton: their geological and chem- ical characteristics – in German].− Mitteilungen der aargauischen Naturforschenden Gesellschaft, 17, 255−320. HARVEY, K. (2007): Healing Touch.− Sawubona, 75−76. HELLMAN, M.j. & RAMSEY, M.S. (2004): Analysis of hot mineral springs and associated deposits in Yellowstone National Park using ASTER and AVIRIS remote sensing.− j. Volcanol. Geotherm. Res., 134/1−2, 195−219. HINTZ, E. & GRÜNHUT, L. (1916): Die Mineralwässer, Moore und Mineralschlame [Mineral waters, moors and mineral muds – in German].– In: DIETRICH, E. & KAMINER, D. (eds.): Handbuch der Balneologie, medizinischen Klimatologie und Balneogeogra- phie.− Georg Thieme, Leipzig, 134−178. HOCHSTEIN, M.P. (1990): Classification and assessment of geothermal resources.– In: DICKSON, M.H. & FANELLI, M. (eds.): Small geothermal resources: A guide to development and utilization, UNI- TAR, New York, 31−57. HORVATINČIĆ, N., SRDOČ, D., PEZDIČ, j., CHAFETZ, H., SLIEPČEVIĆ, A. & KRAjCAR BRONIĆ, I. (1991): Determina- tion of the origin of geothermal waters in NW Yugoslavia by iso- topic methods (Abstract).− XIV International Radiocarbon Confer- ence, Tucson, AZ, USA, 20.05.1991.−24.05.1991., Radiocarbon, 209 p. HORVATINČIĆ, N., SRDOČ, D., KRAjCAR BRONIĆ, I., PEZDIČ, j., KAPELj, S. & SLIEPČEVIĆ, A. (1996): A study of geothermal waters in Northwest Croatia and East Slovenia. Isotopes in water resources management.− Vienna, IAEA, Vol. 2, 470−474. HOSPITAL KRAPINSKE TOPLICE (2014): http://www.sbkt.hr/termal- na.php HOSPITAL STUBIČKE TOPLICE (2014): http://www.bolnicastubick- etoplice.com IVEKOVIĆ, H. & PEROŠ, R. (1981): Mineralne i termalne vode Hrvat- ske [Mineral and thermal waters in Croatia – in Croatian].− Works of Yugoslav Academy of Sciences and Arts, Dpt. of math., phys. and techn. sci. and Dpt. of med. sci., 190 p. jACOBj, C. (1907): Einfache warme Quellen − Akratothermen [Regu- lar warm springs - Akratotherms – in German].– In: Keiserliche Gesundheitsamt (eds): Deutsche Bäderbuch. Verlag von j. j. Weber, Leipzig, 24−40. geothermal aquifer in the study area. However, during the movement of thermal water toward the surface, it comes into contact with limestone and volcanic clastic sediments and processes of cation exchange and silicate weathering take placet. The influence of those processes is recognized in the waters of the Gotalovec, Varaždinske toplice, Podevčevo and Stubičke toplice thermal water sources. ACKNOWLEDGEMENT The authors wish to express their gratitude to Vlasta jURŠIĆ- MITROVIĆ and Vinko MRAZ. This paper represents part of the hydrogeological research conducted within the project “Basic Hydrogeological Map of the Republic of Croatia” (project no. 181-1811096-3165) funded by the Ministry of Science, Education and Sports, Republic of Croatia. REFERENCES ATKINSON, T.C. & DAVIDSON, N. (2002): Is the water still hot? Sus- tainability and the thermal springs at Bath, England.− Geol. Soc., London, Special Publications 193, 15−40. AXELSSON, G. & GUNNLAUGSSON, E. (2000): Background: geo- thermal utilization, management and monitoring.− In: Long-term monitoring of high- and low-enthalpy fields under exploitation, World Geothermal Congress Short Course, Kyushu-Tohoku, May 28−june 10, 2000, 3−10. BAĆ, j. & HERAK, M. (1962): Prijedlog za određivanje užih i širih zaštitnih zona termomineralnih izvora u Hrvatskoj [Proposal to de- termine the immediate and wider protection zone of thermo-miner- al resources in Croatia – in Croatian].– Unpubl. Report, Croatian Geological Survey, Zagreb, 40 p. BENDERITTER, Y. & CORMY, G. (1990): Possible approach to geo- thermal research and relative costs.– In: DICKSON, M.H. & FANELLI, M. (eds.): Small geothermal resources: A guide to de- velopment and utilization, UNITAR, New York, 59−69. BERTRAM, M.A., MACKENZIE, F.T., BISHOP, F.C. & BISCHOFF, W.D. (1991): Influence of temperature on the stability of magnesian calcite.− Am. Mineral., 76, 1889–1896. BITUH, T., MAROVIĆ, G., PETRINEC, B., SENČAR, j. & FRANULOVIĆ, I. (2009): Natural Radioactivity of 226Ra and 228Ra in Thermal and Mineral Waters in Croatia.− Radiation Protection Dosimetry, 133/2, 119−123. BOROVIĆ, S. & MARKOVIĆ, I. (2015): Utilization and tourism val- orisation of geothermal waters in Croatia.− Renewable and Sustain- able Energy Reviews, 44, 52−63. BOŠNjAK, R., ČUBRIĆ, S., GOLUB, M., GRABOVSKI, K. & jELIĆ, K. (1998): Geoen – program korištenja geotermalne energije [Geoen – Programme of geothermal energy utilization – in Croatian].− Energy institute “Hrvoje Požar”, Zagreb, 124 p. BRITVIĆ, V. (1988): Tehnički izvještaj o bušenju i ispitivanju bušotine SUT-3 u Sutinskim Toplicama [Technical report on the drilling and testing of well SUT-3 Sutinskim Toplicama − in Croatian].− Unpubl. Report, Geotehnika, Zagreb, 32 p. CAPAR, A. (1982): Vodoistražni radovi Šemničke Toplice [Water re- search of Šemničke Toplice – in Croatian].− Unpubl. Report, Ge- ofizika, Zagreb, 55 p. DEMING, D. (2002): Introduction to hydrogeology.– McGraw Hill, New York, 468 p. ĐUROVIĆ, R. j. (1960): Definicija i klasifikacija termalnih, termo-min- eralnih i mineralnih voda [Definition and classification of thermal, Marković et al.: Geochemical characteristics of thermal waters of Hrvatsko zagorje Geologia Croatica 77 jANKOWSKI, j., ACWORTH, R.I. & SHEKARFOROUSH, S. (1998): Reverse ion exchange in deeply weathered porphyritic dacite frac- tured aquifer system, Yass, NewSouth Wales, Australia.− In: ARE- HART, G.B. & HULSTON, j.R. (eds.): 9th Internat. Symp. water– rock interaction. Rotterdam: AA Balkema, 243–246. jONES, M. T. & GISALSON, S.R. (2008): Rapid releases of metal salts and nutrients following the deposition of volcanic ash into aqueous environments.− Geochimica et Cosmochimica Acta, 72, 3661–3680. jURIŠIĆ MITROVIĆ, V. (2001): Izvještaj o rezultatima kemijskih anal- iza uzoraka voda za zadatak "Monografija termalnih i mineralnih voda RH"[Report on the results of chemical analysis of water sam- ples for the task "Monography of thermal and mineral waters RC" − in Croatian].− Croatian Geological Survey, Zagreb. KOVAČIĆ, M. & PERICA, R. (1998): Stupanj korištenja geotermalnih voda u Hrvatskoj [The level of thermal water utilization in Croatia – in Croatian].– Hrvatske vode, 25, 355−361. KREŠIĆ, N. (2010): Types and classifications of springs.– In: KREŠIĆ, N. & STEVANOVIĆ, Z. (eds.): Groundwater hydrology of springs. Butterworth-Heinemann, Burlington, 31−85. LAMOREAUX, P.E. & TANNER, j.T. (2001): Springs and Bottled Wa- ters of the World. Ancient History, Source, Occurrence, Quality and Use.− Springer Verlag, Berlin, 250 p. LUND, j.W. (1996): Balneological use of thermal and mineral waters in the USA, Geotherm., 25/1, 103−147. LUND, j.W. (2006): Geothermal energy focus: Tapping the earth's natu- ral heat, Refocus, 7/6, 48−51. MARKOVIĆ, S. & KOVAČIĆ, M. (2006): Izvorište geotermalne vode u Stubičkim Toplicama [The source of geothermal water of Stubičke Toplice – in Croatian].− Hrvatske vode, 55, 173−181. MAROVIĆ, G., SENČAR, j., FRANIĆ, Z. & LOKOBAUER, N. (1996): Radium-226 in Thermal and Mineral Springs of Croatia and Associated Health Risk.− j. Environ. Radioactivity, 33/3, 309−317. MCLEAN, W. & jANKOWSKI, j. (2000): Groundwater quality and sustainability in an alluvial aquifer, Australia.− In: SILILO, A. (ed): XXX IAH congress on groundwater: past achievements and future challenges. A.A. Balkema, Rotterdam. MEYBECK, M. (1987): Global chemical weathering from surficial rocks estimated from river dissolved loads.− Am. j. Sci., 287, 401–428. MIHOLIĆ, S. (1940a): Kemijska analiza termalnih voda u Hrvatskom zagorju [Chemical analysis of thermal waters of Hrvatsko zagorje – in Croatian].− Rad jugosl. akad. znan. umjet., 267, 195−215. MIHOLIĆ, S. (1952): Kemijski sastav i svojstva mineralnih voda [Chem- ical compostion and charachteristics of mineral waters – in Croatian].− God. balneo. klimat. insti. NR Hrvatske, 1, 7−18. MIHOLIĆ, S. (1959): Istraživanje termalnih vrela hrvatskog zagorja [In- vestigation of thermal waters of Hrvatsko zagorje – in Croatian].− Ljetopis jugosl. akad. znan. umjet., 63, 326−328. MORSE, j.W. & ARVIDSON, R.S. (2002): The dissolution kinetics of major sedimentary carbonate minerals.− Earth-Science Reviews, 58, 51–84. MRAZ, V. & KRSNIK, M. (2005): Utjecaj temperature podzemne vode na geoelektričnu otpornost u karbonatnim vodonosnicima na primjeru Harine Zlake – Hrvatsko zagorje [Effect of temperature on groundwater geophysical resistivity in carbonate aquifers in the example of Harina Zlake – Hrvatsko zagorje – in Croatian].− 3. Hrvatski geološki kongres, Opatija. MUFFLER, P. & CATALDI, R. (1978): Methods for regional assessment of geothermal resources.− Geothermics, 7, 53−89. NICHOLSON, K. (1993): Geothermal fluids.− Springer Verlag, Berlin, 264 p. PETRACCIA, L., LIBERATI, G. & MASCIULLO, S.G. (2005): Water, mineral waters and health.− Clin. Nutrit., 25/3, 377−385. POLANČEC, M. (2011): Primjena prirodnih hidrogeokemijskih i izo- topnih trasera u istraživanju geotermalnih vodonosnih sustava [Ap- plication of hydrogeochemical and isotopic tracers in the study of geothermal aquifer systems – in Croatian].− Graduate Work, Geo- technical Faculty in Varaždin, University of Zagreb, 56 p. ROGLIĆ, j. (1974): Obilježja prirodne osnove [Physical geography fea- tures].− In: CRKVENČIĆ, I. (ed.): Geografija Hrvatske – Knjiga I., Školska knjiga, 44–92. ROUTH, H.B., BHOWMIK, K.R., PARISH, L.C. & WITKOWSKI, j.A. (1996): Balneology, mineral water, and spas in historical perspec- tive.− Clinics in Dermatology, 14/6, 551−554. SANYAL, S.K. (2005): Classification of geothermal systems – a possi- ble scheme.− Proceedings of the 13th Workshop on Geothermal Reservoir Engineering, Stanford, California, 85−88. SCHEjBAL, B. (2003): Nova razmatranja o Aquae Balisse i narodu jaza [New considerations about Aquae Balissae and the Iassi people − in Croatian].− Opvuscula archeologica, 27, 393−416. SPECIAL HOSPITAL (2014): http://www.minerva.hr ŠIMUNIĆ, AN. (1986): Detaljna geološka istraživanja bliže okolice Krapinskih Toplica u svrhu lociranja istražne bušotine na termalnu vodu i izvještaj o rezultatima bušenja [Detailed geological investi- gations in the vicinity of Krapinske Toplice to locate exploration well of thermal water and report on the results of drilling − in Croatian].− Unpubl. report, Croatian Geological Survey, Zagreb, 46 p. ŠIMUNIĆ, A. (2008): Geotermalne i mineralne vode RH [Geothermal and mineral waters of RC – in Croatian].− Croatian Geological Sur- vey, Zagreb, 343 p. UNGEMACH, P. & ECONOMIDES, M. j. (1987): Applied geother- mics.− john Wiley & Sons, Chichester, 238 p. WITHAM, C.S., OPPENHEIMER, C. & HORWELL, C.j. (2005): Vol- canic ash-leachates: a review and recommendations for sampling methods.– j. Volcanol. Geoth. Res., 141/3–4, 299–326. YOSHIMURA, K., LIU, Z., CAO, j., YUAN, D., INOKURA, Y. & NOTO, M. (2004): Deep source CO2 in natural waters and its role in extensive tufa deposition in the Huanglong Ravines, Sichuan, China.− Chemical Geology, 205, 141–153. Manuscript received November 14, 2014 Revised manuscript accepted January 24, 2015 Available online February 27, 2015