Lackovic.indd 213 Lacković et al.: Evolution of the Veternica cave (Medvednica Mountain, Croatia) drainage system... � AB STRA CT Detailed fi eld examination, U-Th age dating, and stable isotope analysis of Veternica fl owstone and shelfstone de- posits have been employed in order to unravel details about the geological history of Veternica Cave (Medvednica Mountain, Croatia). The study was carried out in the hydrologically inactive part of the main cave channel, which is developed mainly along the unconformity between Triassic dolostone and Miocene limestone. For 180 m from the cave entrance, (located at 320 metres above sea level (asl)), the morphology of the main channel refl ects exclusively phreat- ic conditions in the cave until the end of its hydrological activity. From 180 to 390 m, the phreatic channel has a second- ary vadose entrenchment in the bottom part as marked by massive fl owstones at elevations from 306 to 313.5 m asl. From 390 m farther inside the cave (in the upstream direction), the main channel has a tall, narrow cross-section and is of mainly vadose origin. In this part of the cave shelfstone precipitates are observed at 9 different levels ranging from 318.8 to 320.2 m asl, indicating the water palaeolevels in the cave. U-Th dating revealed the age of the highest shelfstone (320.2 m asl) of ~380 kyrs BP, and the age of the lowermost analyzed shelfstone (318.9 m asl) of ~245 kyrs BP. Dating of fl owstone deposits, located below the shelfstone level, revealed their formation from ~235 to 205 kyrs BP. A relatively rapid water table lowering and transition from phreatic to vadose conditions occurred within the cave, from ~245 to 235 ka BP, between formation of the youngest shelfstone (representing phreatic conditions below their level) and the oldest fl owstone (marking the beginning of vadose conditions). The results provide unique new information about Veternica Cave genesis and geomorphologic evolution of the Medvednica Mountain area. Keywords: speleogenesis, shelfstone, fl owstone, U-Th dating, water table lowering, Veternica Cave, Croatia Evolution of the Veternica Cave (Medvednica Mountain, Croatia) drainage system: insights from the distribution and dating of cave deposits � Damir Lacković1, Bosiljka Glumac2, Yemane Asmerom3 & Andrej Stroj4 1Croatian Natural History Museum, Zagreb, Croatia; (damir.lackovic@hpm.hr) 2Department of Geosciences, Smith College, Northampton, Massachusetts 01063, USA; (bglumac@smith.edu) 3Deptartment of Earth & Planetary Sciences, University of New Mexico, Albuquerque, New Mexico, USA; (asmerom@unm.edu) 4Croatian Geological Survey, Zagreb, Croatia; (andrej.stroj@hgi-cgs.hr) doi: 104154/gc.2011.18 Geologia Croatica 64/3 213–221 6 Figs. 2 Tabs. Zagreb 2011 Geologia CroaticaGeologia Croatica 1. INTRODUCTION Veternica Cave is located in the southwestern part of the Medvednica Mountain in the vicinity of Zagreb, the capital of Croatia (Fig. 1). The cave is located in a small isolated karst area, separated from the main Dinaric karst region. The cave has been declared a Geomorphologic Natural Monu- ment and is protected within the Medvednica Nature Park. Veternica was discovered in the 19th century (GORJANOVIĆ- KRAMBERGER, 1899), and is named after a local term for wind (veter or vjetar; i.e. the wind cave), because of the con- stant air circulation through the cave channels that can be felt near the cave entrance. The fi rst detailed speleological and geological investigations of the Veternica Cave were Geologia Croatica 64/3Geologia Croatica 214 performed by POLJAK (1934), and later on by MALEZ (1965), ČEPELAK (1977), PANUŠKA & MARJANAC (1977), SUTLOVIĆ (1991), MARJANAC et al. (2005), and MARJANAC (2007), among others. MALEZ (1965) (a pa- laeontologist), conducted systematic research on the cave entrance sediments, during which the natural entrance was widened from 0.45 x 0.29 m (HIRC, 1903) to about 4 x 3 m. This research revealed a rich Late Pleistocene fauna and Ne- anderthal artifacts (MALEZ, 1965). More recent research by MIRACLE & BRAJKOVIĆ (2010) provided a revision of Ple istocene faunal characteristics from the Veternica en- trance sediments and their palaeoclimatic implications. MALEZ (1965) gave a detailed description of clastic deposits at the cave entrance. Analysis of associated fossils suggested the Late Pleistocene age of the deposits, ranging from the end of the Riss glacial to the Holocene. Deposits from the deeper parts of the cave were not described in de- tail, but MALEZ (1965) documented the channel morpho- logy and made a map and cross-section of the channels up to about 2000 m from the cave entrance. MALEZ (1965) also produced a general interpretation of the cave formation through multiple phases of its geological history. He con- cluded that during the oldest phase there was an active spring at the present-day cave entrance, and he described a gradual retreat of the water towards the cave interior to its current sink. Nowadays, there is a great public interest in Veternica. It is the 6th longest cave in Croatia with a total cave channel length of 7128 m (CROATIAN SPELEOLOGICAL SER- VER, 2011), and is easily accessible from Zagreb, the lar gest city in Croatia. The fi rst 380 m of the cave are open for tour- ists and many speleological schools are held in the cave. The current study aimed at gaining detailed information about the history of the cave that could then be incorporated into tourist guides and shared with visitors, such as the nature and timing of palaeoenvironmental changes associated with a water level lowering in the karst drainage system of the cave. This information was primarily gathered through an in-depth fi eld investigation and a detailed study of the cave fl owstone and shelfstone deposits. Field investigation focused on the morphology of the Veternica Cave conduits and associated alluvial sediments in order to determine hydrological conditions in the cave during its hydrological activity (e.g., BÖGLI, 1980; FORD & WILLIAMS, 1989; PALMER, 2007). Phreatic conduits usually have elliptical to round cross-sections, and as they are formed by pressure fl ow, their inclination direction is not constrained by fl ow direction. Vadose channels, on the other hand, are typically vertical shafts or meandering canyons formed by gravity fl ows, and consequently their inclination direction is uniform and governed by fl ow direction. Multi- ple phase channels commonly have phreatic upper parts with vadose trenching in their lower parts (i.e., a typical “keyhole” profi le). This study also focused on the Veternica Cave fl owstones and shelfstones because they can provide important informa- tion for the detailed reconstruction of cave formation and its palaeoenvironments. The cave fl owstones precipitate from thin fi lms of water fl owing over subaerially exposed surfaces, and shelfstones, also known as water-table speleothems, are calcite deposits precipitated at the water surface as ledges attached to the cave walls (HILL & FORTI, 1997). These precipitates can be very good indicators of former water le- vels in caves, as shown by FORD et al. (1993), SZABO et al. (1994), and AULER & SMART (2001). Veternica shelf- stones from the “Majmunski prolaz” channel, at about 450 m from the entrance, and the “Kameni slap” fl owstones in the “Touristic” part of the cave about 250 m from the en- trance (sample locations 1–2 and 3–5; Fig. 2B), were fi rst studied by LACKOVIĆ and HORVATINČIĆ in 1999 (un- publ.), who collected 4 shelfstone samples and 1 sample of the “Kameni slap” fl owstone for radiocarbon dating at the Ruđer Bošković Institute in Zagreb. Since the resulting ages were greater than the limit of the 14C technique at that time (>35 kyrs) no further analyses were conducted. The current study refl ects renewed interest in the geological history of the Veternica Cave made possible through the use of the modern U-Th radiometric dating method. 2. GEOLOGICAL AND HYDROGEOLOGICAL SETTING Unlike the rest of the Medvednica Mountain region where carbonate rocks are associated with siliciclastic, metamor- phic and volcanic rocks, the western part of Medvednica that hosts the Veternica Cave is made up mainly of carbonate rocks. The oldest rocks in the area are Triassic dolostones, dolomitic limestones, and limestones with interbedded chert, shale and pyroclastic deposits, which are overlain by Mio- cene dolomitic breccias and conglomerates as well as bio- clastic limestones (ŠIKIĆ, 1995). The Veternica Cave pas- sages are of varying morphology and are located mainly Figure 1: Location of Veternica Cave in the southwestern part of Medved- nica Mountain. Geologia Croatica 215 Lacković et al.: Evolution of the Veternica cave (Medvednica Mountain, Croatia) drainage system... along the unconformable contact between Triassic dolos- tones and Miocene limestones (MALEZ, 1965; ČEPELAK 1977; Fig. 2A). Some of the passages also follow fractures in both of these deposits. The cave is part of a karst channel system that drains water from sinkholes of the “Ponikve” karst polje and the surrounding karst terrain, through the subsurface to the “Du- bravica” spring (MALEZ, 1963; BOŽIĆEVIĆ, 1976; ČEPELAK, 1977; Fig. 2A), as confi rmed by dye testing (BOŽIČEVIĆ, 1976). There are 9 sinking streams that drain into Veternica. Two of them are located in the Ponikve polje, and the rest are located west of Veternica in the extension direction of the cave lateral channels (Fig. 2B; ČEPELAK, 1977; 1979). There are currently 14 active streams in the cave (ČEPELAK, 1977). The cave entrance is located 320 m asl, and the cave channels range in elevation from 265 m asl (at the bottom of “Ponor 16 m” channel; Fig. 2B) to 435 m asl, yielding a to- tal elevation difference of 170 m. The cave channels up to approximately 900 m from the cave entrance are hydrologi- cally inactive, but their phreatic and/or vadose morphology is very well preserved owing to an absence of collapse pro- cesses and rare speleothem precipitates. From 900 m farther inside the cave, there is an active stream that ends in the “Main siphon” located below the main cave channel (Fig. 2B). Be- sides this main channel, there are several generally narrower lateral side channels (e.g., “Velebitaški” channel; Fig. 2B). Above the main cave channel with active water fl ow (ap- proximately at 950 to 1350 m from the cave entrance), dry chambers and passages occur at higher elevation. Steep side channels with vadose morphology are present below the level of the main channel in several places at 180 to 900 m from the cave entrance (e.g., “Ponor 16 m”, “Vele bitaški si- phon” below “Velebitaški” side channel, “Kukušni” channel, channel towards the “Main siphon”; Fig. 2B). The main sink- holes of Ponikve karst polje are at 470 m asl, while the pre- sent-day “Main” cave siphon is at 295 m asl, and the cave water surfaces at the Dubravica spring at 220 m asl (Fig. 2A). 3. METHODS Fieldwork for this study was carried out in 2008 and 2009 in the hydrologically inactive part of the main Veternica Cave channel. The cave shelfstone and fl owstone deposits Figure 2: Cross-sections of the main Veternica Cave passages. A) Geological profi le illustrating the water drainage pattern from sinkholes along the Ponikve karst polje and through the cave to the Dubravica spring. B) Enlarged area to show detail of the main cave passage and parts of the side channels. Shelfstone (1 & 2) and fl owstone (3, 4 & 5) sampling locations are indicated. Reconstructed palaeolake are also shown. Note vertical exaggeration. Geologia Croatica 64/3Geologia Croatica 216 were described and documented in detail and samples for U-Th age dating and stable isotope (oxygen and carbon) analyses were collected. To precisely measure the elevation of various shelfstone and fl owstone deposits, a detailed top- ographic survey from the cave entrance to the “Main siphon” that marks the present-day water level and the contact be- tween the vadose and phreatic zone in the cave was carried out (Fig. 2B). The topographic survey was conducted using a precise compass and clinometer (Suunto types) and cali- brated rods. Every measurement of polygonal segments was executed both forward and backward to obtain the maximum survey precision. The 234U–230Th dating was done at the Radiogenic Iso- tope Laboratory, the University of New Mexico, USA. Sam- ples were spiked with a mixed 229Th–233U–236U spike. U and Th were separated using conventional anion exchange chro- matography. U and Th isotopes were measured using a Thermo Neptune multi-collector inductively coupled plasma mass spectrometer (MC-ICPMS) which was optimized for U-series analytical work as described by ASMEROM et al. (2006). 234U was measured on a secondary electron multiplier with high abundance fi lter, while the other isotopes of ura- nium were measured on Faraday cups with amplifi ers that had mixed 1010, 1011 and 1012 ohm resistors for 233U and 236U, 235U and 238U, respectively. Mass fractionation was monitored using the 236U/233U ratio, while SEM/Faraday gain was set us- ing sample standard bracketing. A similar procedure was used for Th isotope measurements. We used an initial 230Th/232Th atomic ratio of 4.4x10–6 assuming a source of Th with a bulk earth 232Th/238U ratio of 3.8. The age errors in Table 1 refl ect analytical errors and uncertainties in the value of the initial ratio (± 50%). The laboratory U and Th procedural blanks range from 10–40 pg and 5–20 pg, respectively and were not analytically signifi cant. The CRM145 U isotope standard was measured with the samples obtaining the conventionally accepted δ234U value of –36.5±0.5‰ (CHENG et al., 2000). Samples for stable isotope analysis represent a small amount of carbonate powder collected from fl owstone and shelfstone specimens using a microdrill mounted on a bin- ocular microscope. The fl owstone samples were mainly dril- led along stratigraphic horizons that represent a series of thin laminae made up of fi nely to medium crystalline calcite. The shelfstone samples were also drilled along distinct strati- graphic horizons composed of coarse to very coarse crystal- line calcite with only faint laminations visible in some spec- imens. All samples were heated at 400°C for one hour to remove volatile organic components. At the University of Massachusetts at Amherst, USA, samples were reacted with 100% anhydrous phosphoric acid (H3PO4) for ten minutes and analyzed using an on-line automated carbonate prepara- tion system (Kiell III) linked to a Finnigan-MAT DeltaXL+ ratio mass spectrometer. Standard isobaric and phosphoric acid fractionation corrections were applied to all data. Inter- nal analytical precision, monitored through daily analysis of carbonate standards, is better than or equal to 0.1‰ for both carbon and oxygen isotope values. Stable isotope results are expressed as δ13C and δ18O values in ‰ relative to the Vi- enna PeeDee Belemnite standard (VPDB). 4. RESULTS 4.1. Field investigation The fi rst part of the main Veternica channel, from the cave entrance to 180 m (Fig. 2A) is characterized by wide, circu- lar to elliptical conduits with basal fi lls of clastic alluvial sediments (Fig. 3B). Alluvial sediments at the bottom of the cave channels have been partly dug out by cavers during early cave explorations and also later on to ease tourist ac- cess (Fig. 3A). The overall channel shape and the gently un- dulating surface morphology of the basal alluvial sediments characterized by the absence of uniform slope direction and the lack of evidence for any vadose downcutting, all indicate the presence of phreatic or epiphreatic fl ow conditions in this part of the cave until the end of its hydrological activity. At 180 to 390 m from the cave entrance, the initially phreatic channel is characterized by secondary vadose en- trenchment with a uniform downstream slope direction. This vadose entrenchment ends at the beginning of the “Ponor 16 m” side channel at 180 m from the cave entrance, and the main cave channel contains fl owstone deposits at 250 m from the entrance (Fig. 2B; sampling sites 3, 4 and 5). Figure 3: Characteristic morphology of the cave channels. A) Photograph of a phreatic channel (located 75 m from the cave entrance), covered with clastic sediment in the basal part (the cut through the sediment is man made); B) Sketch of a phreatic channel (located about 130 m from the cave entrance; modifi ed after MALEZ, 1965); C) Photograph of a typical vadose morphology channel (located at 550 m from the cave entrance); D) Sketch of a channel with phreatic roof and tall, narrow vadose trenching morpho- logy (located about 500 m from the cave entrance; modifi ed after MALEZ, 1965). Geologia Croatica 217 Lacković et al.: Evolution of the Veternica cave (Medvednica Mountain, Croatia) drainage system... The third part of the main channel extends farther inside the cave (beyond 390 m), and has a tall, narrow shape with phreatic roof and vadose trenching morphology (Fig. 3D). In this part of the cave, called the “Majmunski prolaz” channel, at 410 to 475 m from the cave entrance there are 9 different levels of shelfstone deposits (Fig 2B; sampling sites 1 and 2). The shelfstones are located along the channel sidewalls and extend towards the middle of the channel from barely notice- able 0.1 cm thick precipitates (Fig. 4A) to a 105 cm long shelf- stone bridging the entire width of the channel. The thickness of individual shelfstones varies from 1 to 30 mm, with an av- erage of about 5 mm. The highest and the lowermost shelfs- tone levels are 20 cm above and 120 cm below the bottom of the present-day cave entrance, respectively. The lateral extent and level of the Veternica shelfstones mark the position of a free surface palaeolake. According to the observed channel elevation as well as a major (vadose to phreatic) change in channel morphology, the palaeolake end ed in the downstream direction at a distance of 390 m from the cave entrance. From that point, all the way to the present en- trance, the cave channels were phreatic during the existence of the palaeolake. Therefore, the formation of the vadose en- trenchment at 180 to 390 m as well as the fl owstone deposition at 250 m from the cave entrance both postdate the palaeolake period. In the upstream direction, the palaeolake extended to 475 m from the cave entrance (Fig. 2B). From that point up- stream and farther inside the cave, the channel bottom is at higher elevation than the highest observed shelfstone level. The fl owstones at the focus of this study are located ap- proximately 250 m from the cave entrance, at an elevation of 306 to 313.5 m asl (Fig. 5). They are below the entrance to a fracture-related side channel (“Velebitaški”), which was the source of the palaeowater fl ow responsible for precipita- tion of the fl owstone known as “Kameni slap” at the eleva- tion of 306 to 312 m asl (Fig. 2B: sampling sites 3 and 4; Fig. 5B). The known thickness of the “Kameni slap” fl ow- stone is 55 cm and its aerial extent is about 20 m2. The sur- face of the fl owstone is covered with karren and small dis- solution scallops (Fig. 5B), and there is a 40 cm deep erosional cut in the middle part of the fl owstone (Fig. 5C). Clastic deposits cover the lower part of the fl owstone and its base. Excavations of these deposits during cave maintenance in 1990 revealed bones of a cave bear (Ursus spelaeus; MA- LINAR & LACKOVIĆ, unpublished data). During the win- ter of 2008/2009 a detailed fi eld study of the “Kameni slap” fl owstone was carried out, and specimens for U-Th age dat- ing were collected by drilling a 15 cm deep core in the deep- est part of the naturally eroded cut, through the fl owstone at sampling site 3 (Figs. 2B, 5D and E). The core, however, did not reach the oldest “Kameni slap” fl owstone precipitate. Sampling site 4 is on the surface of the fl owstone (Figs. 2B and 5C). Besides the “Kameni slap” fl owstone, a sample from an eroded fl owstone above “Kameni slap” (sampling site 5; Figs. 2B and 5A), from the entrance to the “Velebitaški” channel at 313.5 m asl was also collected. However, it was not possible to correlate the erosional remnants of this fl ow- stone to the nearby “Kameni slap” fl owstone deposits. 4.2. U-Th age dating Shelfstone samples from the highest level, or 20 cm above the bottom of the cave entrance (sampling site 1 at 320.2 m asl; Figs. 2B and 4B), and from 110 cm below the bottom of the Figure 4: Field photographs of the Veternica Cave shelfstone deposits. A) Speleologist pointing at multiple shelfstone deposits, about 1 cm thick, along the sidewalls of the “Majmunski prolaz” channel. B) Close-up of the 1 cm thick shelfstone shown in (A) and the location of sample VET08-4 at sampling site 1 (Fig. 2B). C) A ~1 cm thick shelfstone remnant and location of sample VET09-4 at sampling site 2 (Fig. 2B). D) Close-up photograph showing the dis- tribution of shelfstones at 4 diff erent levels. Each shelfstone is about 1 cm thick. Geologia Croatica 64/3Geologia Croatica 218 cave entrance (sampling site 2 at 318.9 m asl; Figs. 2B and 4C) were dated, because the lowermost shelfstone (10 cm be- low the analyzed one) is relatively small and poorly preserved. The highest shelfstone gave an age of 379.8 ± 15.5 kyrs BP, and the age of the lowermost analyzed shelfstone is 244.6 ± 3.6 kyrs BP (Table 1). The fl owstone at the elevation of 313.5 m asl (sampling site 5; Figs. 2B and 5A) gave an age of 234.6 ± 9.1 kyrs BP (Table 1). The oldest “Kameni slap” fl owstone sample from the bottom of the drill core, at the elevation of 306 m asl, formed at 212.4 ± 8.1 kyrs BP (sampling site 3; Figs. 2B, 5D and 5E), whereas the youngest sample (sampling site 4; Fig. 2B and 5C) formed at 205.7 ± 10.9 kyrs BP (Table 1). 4.3. Stable isotope analysis The results of preliminary stable isotope analysis of 7 fl ow- stone and 3 shelfstone samples are included in Table 2 and illustrated in Figure 6. In general, the samples analyzed have Figure 5: Veternica Cave fl owstone deposits. A) Location of sample VET08-5/1 at sampling site 5 from an eroded fl owstone at the entrance to the “Velebitaški” channel (Fig. 2B). B) “Kameni slap” fl owstone in the “Touristic” part of the cave (sampling sites 3 and 4). Karren and small dissolution scallops characterize the surface of the fl owstone. C) Location of sample VET08-5/2 at sampling site 4 (Fig. 2B) on fl owstone surface on top of the erosional cut through the “Kameni slap”. D) Arrow marks the location of sampling site 3 (Fig. 2B), which is a coring site on the bottom of the erosional cut through the “Kameni slap” fl owstone. E) Photograph of the core from sampling site 3, with an arrow indicating location of sample VET08-5/3/1. Table 1: Uranium-series data for selected Veternica Cave fl owstone and shelfstone samples. Material: Site (sample) 238U (ng/g) 232Th (pg/g) 230Th/232Th activity ratio 230Th/238U activity ratio measured δ234U (‰) Initial δ234U (‰) uncorrected age (yrs BP) corrected age (yrs BP) Flowstone: 4 (VET08-5/2) 828 ± 2.4 345737 ± 1236 7.30 ± 0.04 0.996 ± 0.005 125±1 223±7 216073 ± 3411 205692 ± 10880 3 (VET08-5/3/1) 699 ± 1.9 211121 ± 498 9.98 ± 0.04 0.986 ± 0.004 110±1 200±5 219968 ± 3009 212400 ± 8105 5 (VET08 5/1) 787 ± 1.9 318073 ± 813 8.60 ± 0.03 1.138 ± 0.004 216±1 419±11 243227 ± 3273 234644 ± 9144 Shelfstone: 2 (VET09-4) 924 ± 2.2 35260 ± 102 87.25 ± 0.39 1.089 ± 0.004 171±1 341±4 245477 ± 3549 244646 ± 3623 1 (VET08-4) 651 ± 1.5 284977 ± 662 8.15 ± 0.03 1.168 ± 0.005 149±1 435±20 389005 ± 13183 379776 ± 15528 Yrs BP = years before present, where present = AD 2009. Geologia Croatica 219 Lacković et al.: Evolution of the Veternica cave (Medvednica Mountain, Croatia) drainage system... relatively invariable oxygen, and highly variable carbon iso- tope values, characteristic of meteoric precipitation. Five of the fl owstone samples (three from sampling site 3 – deeper section of “Kameni slap” fl owstone; one from sampling site 5 – fl owstone at the entrance to “Velebitaški” channel; and one from sampling site 4 – surface layers of “Kameni slap” fl owstone), cluster towards the most negative carbon and oxygen values measured, and partially overlap with the com- positional fi eld defi ned by shelfstone samples, which on ave- rage have slightly less negative oxygen isotope values (Fig. 6). The remaining two fl owstone samples (both from sampling site 4 – surface layers of “Kameni slap” fl owstone) have sub- stantially less negative carbon isotope values (Fig. 6). 5. DISCUSSION Field observations and age determination of the “Majmunski prolaz” shelfstone (sampling sites 1 and 2; Fig. 2B) yielded information about the timing and lateral extent of vadose and phreatic conditions in the Veternica Cave during the time period of shelfstone formation. Since shelfstones generally form in quiet subsurface lakes and pools, their presence in Veternica is related to the occurrence of standing water mark ing a signifi cant decrease in water fl ow intensity. As the age of both analyzed shelfstone samples (~380 and ~245 kyrs BP) coincides with the Riss glacial maxima (OIS 8 and 10; SIEGERT, 2001), it is likely that the decrease in water fl ow through the cave was related to cold and dry glacial cli- matic periods. Since the study area was located in a perigla- cial environment without any substantial evidence for major glaciations of the Medvednica Mountain and nearby areas (ŠIKIĆ, 1995), the main infl uence of Pleistocene climatic oscillations on Veternica Cave hydrology was probably in terms of changing precipitation rates. It is also possible that during shelfstone formation, the “Majmunski prolaz” chan- nel was a hydrologically less active part of the cave system, and that the main water circulation was taking place in other lateral side channels or in channels at different elevations. Given that the elevation difference between the 9 observed shelfstone levels is only 1.4 m, but their age range is at least 135 kyrs (age difference between the second lowest and the highest shelfstone), the period of shelfstone formation rep- resents a relative stagnation in water level within this karst terrain. Particular shelfstone levels probably refl ect water levels during periods of diminished water circulation alter- nating with more hydrologically active periods. The eleva- tion difference between the shelfstones could be a conse- quence of the changing morphology of a palaeospring, or of very slow gradual karstifi cation. Precipitation from stagnant palaeolake water of rather constant composition is also supported by the relatively uni- form stable isotope values of the examined shelfstone sam- ples (Fig. 6). These values overlap in part with a group of fl owstone samples suggesting similar environmental condi- tions for precipitation of some of these cave speleothems. The slightly less negative oxygen isotope values of one of the shelfstone samples may be related to evaporation or tem- perature variations of the palaeolake waters (Fig. 6). On the other hand, the substantially less negative carbon isotope values of the two examined fl owstone samples (Fig. 6) sug- gest a greater variability in the composition of the water per- colating through the cave under vadose conditions. During the formation of shelfstone in the vadose envi- ronment of the “Majmunski prolaz” channel, the lower part of the cave (Fig. 2B) was characterized by phreatic conditions. Timing of the onset of the fl owstone formation in the upper part of the channel at 250 m from the entrance (Figs. 2B and 5A, sampling site 5), prior to ~235 kyrs BP, indicates that the transition from phreatic into vadose conditions happened during a relatively short time period of less than ~10 kyrs. This relatively fast lowering of the water level in Veternica could have been caused by rapid advance of downward karstifi cation as a consequence of the base level lowering due to enhanced tectonism and uplift of the western Medved- nica Mountain relative to the Sava river valley (regional base level). The present-day distribution of Quaternary deposits Table 2: Carbon and oxygen stable isotope data for selected Veternica Cave fl owstone and shelfstone samples. Material: Sampling site (sample) δ13C (‰ VPDB) δ18O (‰ VPDB) Flowstone: 5 (VET 08 5/1) 4 (VET 08 5/2 a) 4 (VET 08 5/2 b) 4 (VET 08 5/2 c) 3 (VET 08 5/3/1 a) 3 (VET 08 5/3/1 b) 3 (VET 08 5/3/1 c) –11.78 –11.43 –2.38 –7.07 –11.83 –10.85 –11.35 –10.31 –9.94 –8.55 –8.89 –9.12 –9.23 –9.48 Shelfstone: 1 (VET 08-4) 2 (VET 09-4 b) 2 (VET 09-4 c) –11.11 –11.26 –11.23 –9.35 –8.45 –9.06 Figure 6: Results of preliminary stable isotope analysis of the Veternica Cave fl owstone and shelfstone deposits. Geologia Croatica 64/3Geologia Croatica 220 in the area suggests a total Quaternary uplift of about 350 m along the southern fl anks of Medvednica Mt., including the late Pleistocene and Holocene (the last 130 kyrs) uplift of about 60 m (KUK et al., 2000). This supports active tecton- ism during the late Quaternary, when Veternica experienced rapid water table fall. During the period of enhanced verti- cal karstifi cation, those formerly phreatic channels (“Tour- istic” part, Fig. 2B) fi rst became epiphreatic or only occa- sionally flooded, and then completely dried out as the drainage capacity of the channels in dolostone deposits at lower elevations increased. The ages of “Kameni slap” fl owstones from the lower- most part of the “Touristic” part of Veternica indicate the end of signifi cant hydrologic activity at the entrance area of the cave. As it was not possible to sample the oldest “Kameni slap” fl owstone deposits, we could not determine the exact onset of fl owstone formation, but we established that it start ed before ~212 kyrs BP (sampling sites 3 and 4; Fig. 2B). The tectonic uplift fi rst reactivated a vadose fl ow through the “Majmunski prolaz” channel. This fl ow eroded most of the sediment from the lower parts of the channel, as supported by the presence of vadose downcutting of the channel bot- tom in the “Touristic” part of the cave, from the “Majmunski prolaz” channel downstream to “Ponor 16 m” (Fig. 2B). This downcutting was likely coeval with the formation and sub- sequent partial erosion of the “Kameni slap” fl owstone (sam- pling sites 3 and 4; Figs. 2B and 5C). At that time, the “Ponor 16 m” was probably the main drainage channel for vadose fl ow, and part of the cave from “Ponor 16 m” downstream to the cave entrance became completely dry as it is today. The preserved phreatic shape of the channel and the bottom sediments surface morphology between the cave entrance and “Ponor 16 m” (Fig. 2B), indicate phreatic conditions in this part of the channel until the cessation of hydrologic ac- tivity (e.g., PALMER, 2007). The reconstructed timing of the end of hydrologic activity at the cave entrance is in agree- ment with previous studies indicating that the oldest clastic sediments at the cave entrance were deposited after hydro- logic activity ceased at the end of the Riss glaciation (prior to 130 to 150 kyrs BP; MALEZ, 1963). After formation of the youngest “Kameni slap” fl ow- stone (sampling site 4; ~205 kyrs BP), water fl ow from the “Velebitaški” channel temporarily increased. This is refl ected in both erosional and corrosional forms present throughout the entire fl owstone surface and in a 40 cm deep erosional cut through the fl owstone (Fig. 5). Today, water fl ows only in the deeper part of the “Velebitaški” channel, where it sinks into the “Velebitaški siphon” (Fig. 2B). The results of detailed GPS measurements suggest that the Medvednica region is still being actively uplifted (PRIBIČEVIĆ et al., 2007). Present-day cave morphology indicates that lowering of the groundwater circulation level is an ongoing process, gradually advancing upward from the downstream end of the karst system. The morphology of younger generation Veternica Cave channels (Fig. 2B) re- fl ects this advance: steep inclination of the hydrologically active channel towards the “Main siphon”, and the presence of abandoned generations of channels and siphons, which were, during past karstifi cation phases, draining water from the main channel to the lower levels (e.g., “Ponor 16 m”, “Kukušni channel”, Fig. 2B). To improve reconstruction of the geological history of the Veternica Cave, future studies should include the dating of the oldest fl owstones as well as shelfstones from all dif- ferent levels, in addition to detailed sedimentological analy- ses and dating of clastic alluvial deposits from the cave chan- nels. Such multifaceted studies can provide unique and interesting information about the cave history that can then be shared with the public. Numerous visitors are attracted annually to the Veternica Cave, not only because of its easy access and proximity to the large metropolitan area of Za- greb, but also because this is one of the longest caves in Croatia, rich in cave formations, while its deposits have yielded cave bear remains and Neanderthal artifacts. This public interest in the cave is a perfect opportunity to dissem- inate the results of up-to-date scientifi c research on Veter- nica, but it also necessitates the need for preservation and protection of this delicate cave system. 6. CONCLUSIONS 1) The study of cave deposits from the hydrologically inactive part of the main channel in Veternica Cave docu- mented the morphological evolution of the cave and sur- rounding area. 2) The cave entrance is at 320 metres asl. For 180 m from the entrance, there were phreatic conditions in the main cave channel until the end of its hydrological activity, while the main channel at 180 to 390 m from the cave entrance shows phreatic morphology with later vadose entrenchment. In this part of the channel there are fl owstone precipitates at 250 m from the cave entrance at elevations ranging from 306 to 313.5 m asl. 3) Part of the channel that is farther inside the cave and higher upstream, has a tall, narrow morphology of vadose origin. Nine different shelfstone remnants were observed in this part of the channel at elevations ranging from 318.8 to 320.2 m asl. These shelfstones mark the position of a free surface palaeolake that extended from 390 to 475 m from the present cave entrance. 4) U-Th dating revealed the age of the highest shelfstone (from 320.2 m asl) of ~380 kyrs BP, and the age of the low- est analyzed shelfstone (from 318.9 m asl) of ~245 kyrs BP. Dating of the fl owstones in the downstream part of the cave indicates their formation from ~235 to 205 kyrs BP. 5) The time interval between formation of the youngest shelfstone (representing phreatic conditions below their level), and the oldest fl owstone (marking the beginning of vadose conditions) represents a period of relatively fast wa- ter table fall between ~245 and 235 kyrs BP and the conse- quent transition from phreatic to vadose conditions in the “Touristic” part of the cave. This fast water table lowering resulted in the preserved phreatic characteristics of the chan- nel between the cave entrance and “Ponor 16 m” (or for about 180 m from the cave entrance) according to the major change in channel and sediment morphology. Geologia Croatica 221 Lacković et al.: Evolution of the Veternica cave (Medvednica Mountain, Croatia) drainage system... AKNOWLEDGEMENT We would like to thank Dr. Nenad BUZJAK and our caving friends from the Velebit and Samobor caving clubs for helping with fi eld work, the Medvednica Nature Park and the Kahn Liberal Arts Institute at Smith College for fi nancial support, and the Croatian Natural History Museum for general support. Dr. Philippe AUDRA (University of Nice, France) provided useful comments on an earlier version of this manu- script. We also thank two anonymous reviewers for their detailed and thoughtful reviews. REFERENCES ASMEROM, Y., POLYAK, V., SCHWIETERS, J. & BOUMAN, C. (2006): Routine high-precision U-Th isotope analysis for paleo- climate chronology.– Geochimica et Cosmochimica Acta, 70, p. A24. AULER, A.S. & SMART, P.L. (2001): Late Quaternary paleoclimate in semiarid northeastern Brazil from U-series dating of travertine and water-table speleothems.– Quaternary Research, 55, 159–167. BÖGLI, A. (1980): Karst Hydrology and Physical Speleology.– Spring- er-Verlag, Berlin, 284 p. BOŽIČEVIĆ, S. (1976): Hidrologija pećine Veternice kraj Zagreba na planini Medvednici [Hydrogeology of the Veternica Cave near Za- greb on the Medvednica Mt. – in Croatian].– 8. jugosl. geol. kon- gres, Ljubljana, 4, 5–12. CHENG, H., EDWARDS, R.L., HOFF, J., GALLUP, C.D., RICHARDS, D.A. & ASMEROM, Y. (2000): The half-lives of uranium-234 and thorium-230.– Chemical Geology, 169, 17–33. CROATIAN SPELEOLOGICAL SERVER (2011): The caves in Croatia longer than 1000 meters.– [Online] http://www.speleologija.hr/kar- ta.html [Updated 11th July 2011]. ČEPELAK, M. (1977): Novija i buduća speleološka istraživanja u Vet- ernici [New and future speleological explorations in Veternica – in Croatian].– Speleolog, 24–25, 1–8. ČEPELAK, M. (1979): Objašnjenje uz nacrt špilje Veternice [Explana- tion for the topographic survey of Veternica Cave – in Croatian].– Speleolog, Zagreb, 26–27, 33–34. FORD, D.C., LUNDBERG, J., PALMER, A.N., PALMER, M.V., DREYBRODT, W. & SCHWARCZ, H.P. (1993): Uranium–series dating of the draining of an aquifer: The example of Wind Cave, Black Hills, South Dakota.– Geological Society of America Bul- letin, 105, 241–250. FORD, D.C. & WILLIAMS, P.W. (1989): Karst Geomorphology and Hydrology.– Chapman & Hall, London, 601 p. GORJANOVIĆ–KRAMBERGER, D. (1899): Krš Zagrebačke gore [Karst of the Medvednica Mt. – in Croatian].– Hrvatski planinar, 11, 161–164. HILL, C. & FORTI, P. (1997): Cave Minerals of the World, 2nd ed.– Na- tional Speleological Society, Huntsville, 463 p. HIRC, D. (1903): U zapadnom prigorju Zagrebačke gore [Western slopes of the Medvednica Mt. – in Croatian].– Hrvatski planinar, 1–2, 7–10. KUK, V., PRELOGOVIĆ, E., SOVIĆ, I., KUK, K. & ŠARIRI, K. (2000): Seizmološke i seizmotektonske značajke šireg zagrebačkog područja [Seismic and seismotectonic characteristics of the greater Zagreb area – in Croatian].– Građevinar, 52/11, 647–653. MALEZ, M. (1963): Stratigrafska i paleontološka proučavanja diluvijal- nog nalazišta u pećini Veternici (Medvednica) [Stratigraphic and paleontological research of a diluvial site in Veternica Cave (Med- vednica Mt.) – in Croatian].– Unpubl. PhD thesis, University of Zagreb, Croatia, 372 p. MALEZ, M. (1965): Pećina Veternica u Medvednici. I opći speleološki pregled, II stratigrafi ja kvartarnih taložina [Veternica Cave in the Medvednica Mt. I General speleological description, II Stratigraphy of Quaternary deposits – in Croatian].– Acta Geol., 5, 175–229. MARJANAC, T. (2007): Špilja Veternica – geološka riznica [Veternica Cave – a geological treasure – in Croatian].– Speleolog, Zagreb, 55, 40–44. MARJANAC, T., MARJANAC, LJ. & BERMANEC, V. (2005): Veter- nica on Medvednica Mt. near Zagreb – A cave developed in dolo- mites, clastics and pyroclastics. 13th International Karstological School “Classical karst”: Karst in various rocks, Postojna, Sloven- ija. Abstracts for posters, 32–33. MIRACLE., P.T. & BRAJKOVIĆ, D. (2010): The palaeoecological sig- nifi cance of the Pleistocene mammalian fauna from Veternica Cave, Croatia. Revision of the lagomorpha, canidae, mustelidae and feli- dae.– Geologia Croatica, 63/2, 207–224. PALMER, A. (2007): Cave geology.– Cave books, Cave Research Foun- dation, USA, 454 p. PANUŠKA, LJ. & MARJANAC, T. (1977): Rezultati nekih novijih geoloških istraživanja u spilji Veternici [Results of new geological research in Veternica Cave – in Croatian].– Speleolog, 24–25, 13–14. POLJAK, J. (1934): Pećina Veternica u Zagrebačkoj gori [Veternica Cave on the Medvednica Mt. – in Croatian].– Priroda, 5, 133–139. PRIBIČEVIĆ, B., MEDAK, D. & ĐAPO, A. (2007): Progušćenje točaka geodinamičke mreže grada Zagreba u podsljemenskoj zoni [Addi- tion to geodynamical monitoring points for the city of Zagreb in the Medvednica foothill area – in Croatian].– Geod. list, 4, 247–258. SIEGERT, M.J. (2001): Ice Sheets and Late Quaternary Environmental Change.– John Wiley & Sons Ltd, London, 231 p. SUTLOVIĆ, A. (1991): Novi metri u staroj Veternici [New measure- ments in the Veternica Cave – in Croatian].– Velebiten, 8, 34–36. SZABO, B.J., KOLESAR, P.T., RIGGS, A.C., WINOGRAD, I.J. & LUDWIG, K.R. (1994): Paleoclimatic inferences from a 120,000- yr calcite record of water-table fl uctuation in Browns Room of Dev- ils Hole, Nevada.– Quaternary Research, 41, 59–69. ŠIKIĆ, K. (1995): Prikaz geološke građe Medvednice [Geological over- view of the Medvednica Mt. – in Croatian].– In: ŠIKIĆ, K (ed.): Geološki vodič Medvednice. Hrvatski Geološki Institut, Zagreb, 7–30. Manuscript received July 11, 2011 Revised manuscript accepted September 01, 2011 Available online October 26, 2011