PHYSICAL AND CHEMICAL RESEARCH IN VELEBITA PIT (CROATIA) FIZIKALNO KEMIČNE RAZISKAVE V BREZNU VELEBITA (HRVAŠKA) Dalibor PAAR1-5, Magdalena UJEVIČ2, Darko BAKŠIČ3'5, Damir LACKOVIČ4'5, Ana ČOP5 & Vanja RADOLIČ6 Abstract UDC 551.44(497.5):546.296 Dalibor Paar, Magdalena Ujevič, Darko Bakšič, Damir Lackovič, Ana Čop & Vanja Radolič: Physical and Chemical Research in Velebita pit (Croatia) We have performed measurements of the cave microclimate, water quality parameters and radon concentration to a depth of -1000 m in Velebita Cave system (Northern Velebit, Croatia). The results were analyzed as a function of the cave depth and geomorphological characteristics. Two different air temperature gradients were obtained, which can be attributed to a cave morphology and air circulation in the upper part of the pit. The water quality parameters show that the studied waters are poorly mineralized and are of weakly alkaline type. Water chemistry is probably predominantly controlled by the petrog-raphy of the bedrock (limestone) and the cave morphology. Water in the cave is not affected by pollution. The average value of radon concentration is rather low, much lower than in some other Croatian caves. Keywords: Cave climate, pH, TDS, conductivity, ORP, radon-concentration, Mt.Velebit, Croatia. Izvleček UDK 551.44(497.5):546.296 Dalibor Paar, Magdalena Ujevič, Darko Bakšič, Damir Lackovič, Ana Čop & Vanja Radolič: Fizikalno kemične raziskave v breznu Velebira (Hrvaška) V jamskem sistemu Velebita (Severni Velebit, Hrvaška), smo do globine 1000 m merili klimatske in hidrološke parametre ter koncentracijo radona. Rezultate smo interpretirali z ozirom na globino in morfološke značilnosti jame. V jami smo izmerili dva različna temperaturna gradienta, ki jih povezujemo z jamsko morfologijo in kroženjem zraka v zgornjem delu brezna. Fizikalno-kemični parametri vode kažejo na nizko mineralizacijo sicer šibko alkalne vode. Kemija vode je v glavnem odvisna od petrografskih značilnosti kamnine (apnenec) in morfologijo jame oz. vodnih poti. Voda v jami ni onesnažena. Povprečna vrednost koncentracije radona je nizka, vsekakor nižja od drugih jam na Hrvaškem.. Ključne besede: Jamska klima, pH, TDS, specifična električna prevodnost, koncentracija radona, Velebit, Hrvaška. introduction Northern Velebit is a mountainous region in Croatia between the Adriatic Sea and Ličko-Gacko Polje. The characteristics of this region place Mt. Velebit to alpine karst. Alpine karst is an important, but also very vulnerable natural resource with large amounts of groundwater stored in its depths. Alpine karst aquifers are complex and very vulnerable systems. To asses the vulnerability of alpine karst water resources and plan their protection, it is necessary to analyze water dynamics, the quality and the diverse characteristics of these systems (Petrič 1 Faculty of Science and Mathematics, Bijenička 32, Zagreb, Croatia, dpaar@phy.hr 2 Croatian National Institute of Public Health, Rockefellerova 7, Zagreb, Croatia 3 Faculty of Forestry, Svetošimunska 25, Zagreb, Croatia 4 Croatian Natural History Museum, Demetrova 1, Zagreb, Croatia 5 Speleological Section Velebit, Radiceva 23, Zagreb, Croatia 6 Department of Physics, University of Osijek, P.O. Box 125, Osijek, Croatia Received/Prejeto: 9.9.2008 2004, Ravbar 2007). As in Slovenia, alpine karst waters in Croatia can be a very important water resource in the future. The amount and quality of groundwater stored in its depths should be explored. The quality of groundwater is largely affected by processes occurring in the vadose zone of karstic aquifers. Deep and accessible pits provide excellent natural laboratories for these investigations. geomorphological characteristics of velebita pit The Northern Velebit karst area is conducive to the formation of deep pits, three of which are deeper than 1 km - Lukina Jama Cave (-1392 m), Slovačka Jama (-1320 m) and the recently discovered Cave system Velebita (-1026 m). So far, 259 pits have been discovered in the Northern Velebit area, which extends over some 25 square kilometers (Bakšic&Paar 2006). This area is composed of lithostratigraphic units ranging from the Middle Triassic to the Paleogene age. The largest part of Northern Velebit is composed of Jurassic sediments which contain carbonate rocks only. Limestone prevails in the composition of deposits, but dolomites are also present. Jurassic deposits are approximately 2,850 m thick. However, speleologi-cal explorations have revealed more complex geology of Northern Velebit. For example, Lukina Jama contains deposits of carbonate breccias from -450 to -700 m and from -750 to -950 m. A similar situation occurs in Slovačka Jama as well. This poses a number of questions related to their stratigraphy and tectonic movements (Lackovic et al. 1999). A significant part of the subject area is covered with the Jelar formation of the Upper Paleogene age. Its origin is closely related to strong tectonic movements affecting the area at that time (Bahun 1974). In the hinterland (Lika) calcareous breccias are partially permeable but in the higher positions on Northern Velebit they are highly permeable. Calcareous breccias are up to 300 m thick (Kuhta & Bakšic 2001). The morphology of the Cave system Velebita is complex to the depth of -228 m. The entrance to Velebita pit descends between lodged boulders to the entrance chamber. The cave continues with the shaft after a passage between the boulders and the compact rock on the floor. The entrance to Dva Javora (Two Maples) pit is a vertical crack of 5x0.5 m situated 66 m NW and 7.3 m below the entrance to Velebita pit. At a depth of -74 m, in the southern part of the main passage, there is a large passage to a parallel meander which goes back south-wards to the neighboring pit Breskva. In Dva Javora there is also another vertical part down to -228 m. The entrance parts of Velebita are in the S-N and S/SW - N/NE directions, while the entrance parts of Dva Javora are in the W-E direction. At the contact of the passages (meanders) in Dva Javora and Velebita stretching through high and wide cracks in the W-E and the S-N directions, there is an entrance to the shaft of Divka Gromovnica (the total shaft length is 513 m). The bottom part of the shaft at -574 m ends with a sloped chamber covered with big boulders. In the southern part of the chamber there is a 10 m wide, shallow drip-water pool, and in the south-eastern part there is another small shallow drip-water pool with a small stream flowing to the bottom of the chamber. In the southern rock of the Divka Gromovnica shaft at -534 m there is a meander in the south direction going down to Bivouac 2 at -786 m. In the vicinity of this point there is an incoming shaft with a stream of water which is strong during rainy weather. The stream continues to the main passage, which goes down to -941 m. After the "Ljutiša" waterfall at -941 m, it joins the SW passage with another incoming water stream from the north direction. Down to the stream is a meander finishing with the „Podzemljar" Chamber of 70 x 60 m, whose bottom is at a depth -1026 m. materials and methods In situ measurements of water pH, specific electric conductivity (COND), total dissolved solids (TDS), oxidation reduction potential (ORP) and temperature were performed in 14 measuring points in Velebita pit using the Ultrameter II 6P Myron device. Measurements were performed during an expedition in Velebita pit (5 - 12 Aug 2007). Two water samples were taken for laboratory analy-sis at the Croatian National Institute of Public Health. Bicarbonate was determined by titration with HCl and Cave system VeEebita North Voíeblt, Croate Hali kova igranja dcplh; -10H ai horiznnlal IcnyLh SG5 ni laiÉl lerliflh: ZSfiO m ExpleníJ 2003-Z0P?. Üy: SO PDS Velebit, SO HPD DuJiovac, SO HP D Mosor, SO HPD Željezničar, Ttie Spo Ico logicai Committse of the Croaüan Mountarneermg Association Tapo: Matija Čepelak. Daiibor Paar, Ronald îeleinjak. Darko Bafcšič. Lovro Čepe Lik. Tihana Dasmič, Vesna HrdNika, Luka Mudronja. Filip HNpovIt, Donatjcj lauti t Ediled by: Darko Bak£ič and Lovro Čepela k Enlrjncc Dva javor* \ ^ dvoran» Enlrance Velebila *= 4557,021 y = S499,2B9 z = 15S7 m ne Id mark coseros Etili sut c Dva jsvera r » 5494,244 i - tSsO m P Jt3 Oltt loo FU Enlrance Dva javora Bivak ; Entrane* V»l»bH« . Pifuri&v odvajali Vflrtlkals ■ - Olvj« Gromovfiic« ■"J Bivak i1 0 m IDO m 200 m 100 m 400 m 500 m £00 m -1026 m Fig. 1: Velebita Cave system. Insert: Location of the cave. methyl orange as an indicator. An ICS 500 ion chroma-tography system (Dionex) was used to determine major anions (sulphate, chloride and nitrate) and cations (sodium, potassium, calcium and magnesium) concentrations. Microclimatic parameters were measured with Onset Hobo Temp/RH and Oakton RH/Temp data loggers. Airflow was measured with Silva Alba wind watch. Integrated measurements of radon and its short-lived progenies in the air were performed by means of the passive track etching method with the LR-115 SSNT detector, type II (Kodak-Pathé, France) at 10 measuring sites. The cylindrical detector cup, with the diameter and length of 11 cm and 7 cm respectively, was either covered with a paper filter of 0.078 kg/m2 surface density (diffusion detector), or was open. Radon concentration in the air was determined as a product of the sensitiv-ity coefficient and track density of the diffusion detector. The measurement method with two detectors (diffusion and open one) enables determination of the equilibrium factor for radon and its progeny in the air (Planinic et al. 1997). results of measurements and discussion MICROCLIMATE CONDITIONS In order to obtain a complete picture of microclimatic parameters, monitoring should be undertaken over an extended period of time (Cigna 2002). However, deep pits such as Velebita do not allow this due to a very com-plex access to the measuring points. Consequently, the cave microclimatic parameters were monitored during Fig. 2: Air and water temperature variation with the cave depth (5.-12.8.2007.). Dashed lines are two different linear regressions of air temperature between -30 and -100 m, and -100 and -800 m. Insert: Relative air humidity RH% as a function of the cave depth. speleological exploration. The temperature and relative humidity of the cave atmosphere were measured at 14 locations throughout the cave (Fig. 2). Relative air humidity in the entrance chamber of Velebita was 95.6 ± 6.3 %, and at other measuring points it was close to 100%. The highest air temperature with the average value 10.5 ± 1.1°C was recorded in the entrance chamber. The temperature falls to 3.4°C at a depth of -95 m. In the shaft of Divka Gromovnica the temperature rises to 4.3°C at the bottom of the shaft at -570 m and continues to rise at a rate given in Table 1. Water temperature also rises at a lower rate. At the entrance to the first shaft in Velebita pit, the air flow changes with a maximum speed of 2.5 ± 0.2 m/s. At -35 m, air flow with change of direction was detected with period between 10 and 20 seconds. The maximum speed was 5.6 ± 0.2 m/s (direction: cold air going out) and 3.5 ± 0.1 m/s (direction: hot air going in). WATER QUALITY Water quality parameters were measured in situ at 14 locations in the cave (Table 2). The pH of meteoric water (collected at the surface) is 6.98 ± 0.01, and in the cave it slightly decreases with the depth, and it is between 8.0 and 8.35 at maximum. The oxidation reduction potential (ORP) value decreases with the depth at a rate -5.5 ± 0.7 mV/100 m. The lowest value of 186 ± 1 mV is at -930 m. The ORP value of meteoric water 206 ± 1 mV is comparable with the value in the cave at -760 m. The specific electrical conductivity (COND) gives general information on water chemistry. The average value of COND is 182 ± 25 ^S/cm and does not show any sig- Table 1: Air and water temperature gradients Depth Temperature gradient Air -30 to -100 m dT/dh = -3.5 ± 0.2 oC / 100 m Air -100 to -800 m dT/dh = +0.25 ± 0.03 oC / 100 m Water -100 to -900 m dT/dh = +0.19 ± 0.02 oC / 100 m. Table 2: In-situ measurements of water quality parameters At depth 0 the meteoric water sample is collected on the surface (rainwater). Location: pH COND ORP TDS T Depth (m) (yS/cm) mV PPm (°C) 0 6.98 13.12 436 8.33 15.9 -35 7.95 183 492 120.9 8.1 -65 8.22 178.5 459 118.7 4.0 -95 8.35 206.2 465 137.8 3.6 -390 8.11 178.5 466 118.2 5.9 -570 8.19 156.3 467 103.5 4.3 -570 8.22 153.1 460 101.4 4.2 -570 8.28 177.8 463 118.3 4.2 -675 8.20 192.4 469 128.0 4.7 -760 8.14 181.5 434 120.6 4.8 -785 8.10 260.6 430 173.4 4.8 -895 8.09 171.9 429 113.9 5.0 -915 8.07 166.8 423 110.4 5.1 -930 8.19 169.5 416 112.1 5.2 Table 3: Chemical analysis of two water samples: a) water sample taken at -584 m, b) at -920 m in velebita pit. Location: Depth (m) hco3 Cl- N03 O2- Na+ K+ Mg2+ Ca2+ Y (mg/L) Y (mg/L) Y (mg/L) Y (mg/L) Y (mg/L) Y (mg/L) Y (mg/L) Y (mg/L) a) -584 92.4 1.56 3.58 2.74 0.5 0.05 0.12 33 b) -920 105.2 1.71 4.16 2.92 0.77 0.08 0.19 36.7 saturated with calcite (SIc= -0.04) while water at -920 m is lightly saturated with calcite (SIc=0.03). The observed SIc values are very close to equilibrium (SI=0). Waters at -584 m and -920 m were unsaturated regarding to dolomite (SId = -2.50 and -2.21 respectively). RADON CONCENTRATION Radon is a naturally occur-ring radioactive gas that is produced in the radioactive decay of uranium and other elements. Radon concentration measurements were performed in the upper part of Velebita pit and. The results are given in Table 4. The measurement time was de-fined by cave exploration, 9 days for the locations down to 100 m, and 6 days for deeper locations. We can interpret these results as measures of lower radon concentrations in the cave. nificant différences between waters from différent locations. However, electric conductivity has a lower value 153 ± 0.1 ^S/cm at a depth -570 m and the highest value 260.6 ± 0.1 ^S/cm at -785 m. The observed changes re-flect the impact of cave litology, geomorphology and the mixing of waters from different flows on water chemis-try. Total dissolved solids (TDS) also have the lowest value of 101.4 ± 0.1 ppm at -570 m, and the highest value of 173.4 ± 0.1 ppm at -785 m. The COND and TDS values of meteoric water are about 10 to 20 times lower than the values in the cave. Data related to the chemical content for only two samples are available (Table 3); however, some general statements can be pointed out. The concentration in Ca2+ largely exceeds that of Mg2+, Na+ and K+. The predominant anion is HCO3- so water in the cave is of Ca-HCO3 type. This water type reflects the predominant rock type in the cave - limestone. Even if dolomite can be present too due to short water-rock interaction time, Mg concentrations are low. Water at -584 m is lightly un- Table 4: Measurements of radon concentration c and equilibrium constant F at various depths in velebita pit Location: c ± ac (Bq m-3) F ± ar Depth (m) -5 96 ± 36 - 38 112 ± 39 0.62 ± 0.38 -68 421 ± 53 -86 255 ± 45 0.20 ± 0.13 -94 24 ± 35 -157 319 ± 67 0.56 ± 0.22 -207 269 ± 64 0.43 ± 0.22 -227 408 ± 69 0.36 ± 0.15 -294 192 ± 61 0.91 ± 0.40 -370 456 ± 71 conclusions The results were analyzed as a function of cave depth and geomorphological characteristics. Two air temperature gradients were obtained: dT/dh = -3.9 ± 0.2°C/100 m and +0.25 ± 0.03°C/100 m. Since at about -100 m there is an entrance to a 513 m long shaft, this change can be attrib-uted to a cave morphology and detected air circulation in the upper part of the pit. In the upper part of the pit, air circulation changes the temperature gradient by condensation or evaporation that is releasing or absorbing the heat. In the lower part, the temperature increases with depth because of Earth's geothermal gradient. It is also shown, that small water masses flowing into a mountain are able to peturb the rock temperature (Badino 2005). The water quality parameters measured in situ at 14 locations in the cave show that the studied waters are poorly mineralized and weakly alkaline. Oxidation reduction potential (ORP) shows clear reduction in the value with depth, while the calculated rate is -5.5 ± 0.7 mV/100 m. Water chemistry is probably predomi-nantly controlled by the petrography of the bedrock (limestone) and the cave morphology. Chemical analyses of water performed on two samples show that waters are Ca-HCO3 type. Low concentrations of NO3-, Cl-, SO42-, K and Na suggest that the water in the cave is not affected by pollution. However, more frequent and continuous research should be undertaken in order to determine fac-tors which control the chemical composition of water. The average value of radon concentration from the entrance down to a depth of -100 m is 181.6 ± 41.6 Bq/ m3, and between -100 m and -400 m is 328.8 ± 66.4 Bq/m3. Radon concentrations are much lower than in some other Croatian caves, where concentrations up to 14.400 Bq/m3 were measured (Paar et al. 2005). acknowledgements The authors would like to thank to North Velebit National Park and Croatian waters for the support of presented work. Their thanks also go to the members of the Spe- leological Section Velebit, Zagreb, for their assistance in field work. references Badino, G., 2005: Underground drainage systems and geothermal flux. Acta carsologica, 34/2, 277-316, Ljubljana. Bahun, S., 1974: The tectogenesis of Mt. Velebit and the formation of Jelar deposits. Geološki vjesnik 27, Zagreb : 35-51. Bakšič, D. & Paar, D., 2006: Croatia and the Deep Caves of Northern Velebit. Alpine Karst, vol 2., ed. J. & T. Oliphant, 105-124, Cave books, Dayton, USA. Cigna, A.A., 2002: Modern trend in Cave monitoring. Acta carsologica 31/1, 35-54, Ljubljana. Kuhta, M. & Bakšič, D., 2001: Karstification Dynamics and Development of the Deep Caves on the North Velebit Mt. - Croatia. 13th Internatinal Congress of Speleology, 1-4, Brazil. Lackovič, D., Šmida, B., Horvatinčic, N. & Tibljaš, D., 1999: Some geological observations in Slovačka jama cave (-1268 m) in Velebit mountain, Croatia. Acta carsologica 28/2, 113-120, Ljubljana. Paar, D., Radolič, V. & Buzjak, N., 2005: Measurements of concentration of radioactive gas radon in the caves of Žumberak. Znanstveno-stručni skup istraživača krša Žumberačke gore, 20, Bregana. Petrič, M., 2004: Alpine karst waters in Slovenia. Acta carsologica 33/1, 11-24, Ljubljana. Planinič, J., Radolič, V., Faj, Z. & Šuveljak, B., 1997: Radon equilibrium factor and aerosols. Nuclear Instruments and Methods A, 396, 414-417. Ravbar, N., 2007: The protection of karst waters - a com-prehensive Slovene Approach to vulnerability and contamination risk mapping. Carsologica 6 - series ed. F. Gabrovšek. Postojna-Ljubljana.