MASARYK UNIVERSITY FACULTY OF SCIENCE DEPARTMENT OF GEOLOGICAL SCIENCES
Hydrogeochemistry of Dripwaters in Selected Caves of Moravian Karst Ph D Dissertation
Pavel. . Pracný
SUPERVISOR: DOC. ING. JIŘÍ FAIMON, DR. BRNO 2017 BIBLIOGRAPHIC ENTRY
Author Mgr. Pavel Pracný Faculty of Science, Masaryk University Department of Geological Sciences
Title of Dissertation Hydrogeochemistry of Dripwaters in Selected Caves of Moravian Karst
Degree Programme Geology
Field of Study Geological Sciences
Supervisor doc. Ing. Jiří Faimon, Dr. Faculty of Science, Masaryk University, Department of Geological Sciences Faculty of Science, Palacký University, Department of Geology Academic Year 2016/2017
Number of Pages 53+65
Keywords Cave dripwater; Moravian Karst; Anomalous drip; Car- bon dioxide; Mg/Ca ratio; Kinetic modeling; Limestone dissolution; Degassing; Mg-calcite BIBLIOGRAFICKÝ ZÁZNAM
Autor Mgr. Pavel Pracný Přírodovědecká fakulta, Masarykova univerzita Ústav geologických věd
Název práce Hydrogeochemie skapových vod ve vybraných jeskyních Moravského krasu
Studijní program Geologie
Studijní obor Geologické vědy
Školitel doc. Ing. Jiří Faimon, Dr. Přírodovědecká fakulta, Masarykova univerzita, Ústav geologických věd Přírodovědecká fakulta, Univerzita Palackého, Katedra geologie Akademický rok 2016/2017
Počet stran 53+65
Klíčová slova Skapové vody; Moravský kras; Anomální skap; Oxid uhličitý; Poměr Mg/Ca; Kinetické modelování; Rozpouštění vápenců; Odplyňování; Mg-kalcit ABSTRACT
Karst dripwaters are an important factor of speleothem formation. These cave precipitates provide various proxy data (e.g. stable isotopes, minor and trace elements or grow laminae) about paleoenvironment. To better under- stand the interrelationship between proxies and environment, an investigation of recent karst processes is important. A dripwater hydrogeochemistry and cave PCO2 were studied in the dry part of Punkva Caves (Moravian Karst). The sampling was conducted twice per month from February 2012 to March 2013. Additional dripwater samples for stable isotopes analyses were collected in April and November 2014. An anomalous drip was identified showing hydrogeochemical properties significantly different from other regular drips in the cave system as well as other caves in Moravian Karst. The anomalous drip showed a low SIcalcite ~ 0.14±0.11 (standard deviation), low specific conductivity 297±22.2 μS cm−1 and enhanced values of δ13C (−7.85 to −8.35‰ VPDB), Mg/Ca × 1000 ratio (45.7±3.3) and Sr/Ca × 1000 ratio (0.65±0.06). In contrast, the regular drips showed satu- ration SIcalcite in range from 0.83 to 1.07, high specific conductivity (604±32 μS cm−1) and lower Mg/Ca × 1000 (17.0±1.4) and Sr/Ca × 1000 (0.31±0.02) ratios as well as lower δ13C values (−10.34 to −10.94‰). The data analysis supports conclusion that the anomalous drip properties were a consequence of a prior calcite precipitation or/and water mixing. This idea is supported by the position of the drip on a crevice edge.
The partial pressure of CO2 measured in cave air, PCO2(air), was in range from 10−3.31 to 10−2.49 (0.06–0.32 vol%). These values were compared to CO2 par- tial pressures calculated from dripwater hydrogeochemistry as a partial pres- sure of CO2 corresponding to aqueous carbonates, PCO2(W) (10−2.91 to 10−2.35 or
0.12–0.45 vol%), and hypothetical CO2 partial pressure participating on the in- itial dripwater formation, PCO2(H) (10−1.77 to 10−1.49 or 1.7–3.2 vol%). Both the
PCO2(air) and PCO2(W) showed clear seasonal variations with maxima in summer and minima in winter. It seems that the cave air CO2 had been controlled by cave ventilation modes: the higher PCO2(air) were a result of a downward airflow mode during the period of active ventilation with increased influx of CO2 from epikarst and vadose zone. In contrast, the PCO2(H) was very stable without any significant seasonality indicating independence on seasonally changing surface conditions. It could mean that the source of CO2 is deployed deeper in karst profile under the soil. The anomalous drip represented an exception with lower and varying PCO2(W) and PCO2(H) close to PCO2(air) indicating prior CO2 degassing and calcite precipitation. A geochemical model of CO2 degassing shows that the regular dripwaters data are plotted along a degassing line with slope ~ −1 pointing to a unique value of PCO2(H) regardless of season. In addition, it shows that the anomalous drip data are much more scattered and estimated values of PCO2(H) are most probably incorrect due to previous calcite precipitation changing dripwater hydrochemical properties. The possibility of dripwater conversion into solution aggressive with re- spect to calcite due to anthropogenic CO2 influx into cave were studied in Výpustek Cave. The model showed that it is possible to reach sufficient cave
CO2 concentrations during longer events with enhanced attendance (500 peo- ple). Ordinary guided tours (50 people, ca. 0.5 h) seem to be of inconsequential effect. A dynamic model of the Mg/Ca ratio theoretical evolution during lime- stone dissolution in epikarst (T = 10 °C, logPCO2 = –1.5) was designed. The lime- stone was modeled as a dolomite and Mg-calcite mix with various content ra- tios. Two distinct stages of the dissolution were observed: (a) an initial stage with stoichiometric release of Ca and Mg (congruent dissolution) and (b) an advanced stage beginning when the solution reached calcite saturation, char- acterized by a continuous release of Mg and concurrent Ca decrease due to cal- cite precipitation (incongruent dissolution). The overall Mg/Ca ratio evolution, represented by shape of dissolution reaction paths, is determined by the Mg- calcite composition and the ratio of Mg-calcite and dolomite (D/C). The dynam- ics of Mg-calcite dissolution dominates for all ratios under D/C = 1, when the reaction paths divert from pure Mg-calcite paths. A minor factor influencing the reaction path of Mg/Ca evolution was identified in the ratio of limestone surface to water volume ({L}/V). However, the {L}/V ratio controls overall inter- action dynamics. In epikarst, the dissolution dynamics is given by conditions of a system open to gaseous CO2, leading to enhanced epikarst dissolution. The ratio is probably higher deeper in vadose zone, but the dissolution is limited because the system is closed to CO2. The actual Mg/Ca ratio in dripwater de- pends on water residence time (i.e., water-rock interaction time) which controls how far the dissolution proceeds along the reaction path. Dripwaters from Punkva Caves and other sites over the world were compared with the model results. Most of the waters showed Mg/Ca ratios similar to reaction paths for Mg-calcite and low-dolomite limestone, whereas dripwaters from dolostone were similar to evolution during dissolution of pure dolomite. Both the difference in Mg/Ca ratios of the anomalous drip compared to regular drips and the Mg/Ca evolution model show the importance of water flow paths in dripwater formation. However, the water flow paths could change both temporally and spatially with the evolution of karst system, inde- pendently on climate conditions. In addition, the Mg-calcite composition and dissolution dynamics seem to play substantial role in dripwater Mg/Ca ratio evolution. Therefore, it is important to take these factors into consideration in paleoenvironmental studies of karst proxies. ABSTRAKT
Krasové vody představují důležitý činitel při vzniku speleotém. Tyto jes- kynní sedimenty jsou významným zdrojem proxy dat (např. stabilní isotopy, stopové prvky, přírůstkové linie) o environmentálních podmínkách v době svého vzniku. Pro lepší pochopení vazeb mezi proxy daty a klimatem je důležité studium recentních krasových procesů. Hydrogeochemické vlastnosti skapo- vých vod a PCO2 v Punkevních jeskyních (Moravský kras) byly studovány dva- krát měsíčně od února 2012 do března 2013. Dodatečné vzorky pro izotopické analýzy byly odebrány v dubnu a listopadu 2014. V Punkevních jeskyních byl identifikován anomální skap s výrazně odliš- nými vlastnostmi ve srovnání s běžnými skapy v této i jiných jeskyních Morav- ského krasu. Anomální skap vykazuje nízké přesycení ke kalcitu SIcalcite ~ 0,14±0,11 (směrodatná odchylka) a nízkou specifickou vodivost 297±22,2 μS cm−1 a zvýšené hodnoty δ13C (−7,85 až −8,35 ‰ VPDB), poměru Mg/Ca × 1000 (45,7±3,3) a poměru Sr/Ca × 1000 (0,65±0,06). Běžné skapy vykazují hodnoty indexu nasycení SIcalcite rozmezí od 0,83 do 1,07, vysokou specifickou vodivost (604±32 μS cm−1) a nižší hodnoty δ13C (−10,34 to −10,94 ‰) a poměrů Mg/Ca × 1000 (17,0±1,4) a Sr/Ca × 1000 (0,31±0,02). Data naznačují, že vlastnosti ano- málního skapu jsou důsledkem předběžného srážení kalcitu a/nebo mixování vod v nadloží, což podporuje také pozice skapu na hraně komína.
Praciální tlak CO2 v jeskynním vzduchu, PCO2(air), se pohyboval v rozmezí od 10−3,31 do 10−2,49 (0,06–0,32 obj. %). Tyto hodnoty byly porovnány s parciál- ními tlaky vypočítanými na základě hydrogeochemických vlastností skapových vod. Jsou to parciální tlak odpovídající obsahu karbonátů v roztoku, PCO2(W) (10−2,91 až 10−2,35, tj. 0,12–0,45 obj. %), a hypotetický parciální tlak podílející se na formování skapové vody, PCO2(H) (10−1,77 až 10−1,49, tj. 1,7–3,2 obj. %). Jak
PCO2(air) tak i PCO2(W) vykazují jasnou sezónnost s maximy v létě a minimy v zimě. Toto chování je patrně dáno jeskynní ventilací: vyšší PCO2(air) je důsled- kem režimu sestupného proudění v období aktivní jeskynní ventilace se zvýše- ným přínosem CO2 z epikrasu a vadózní zóny. Naopak hodnoty PCO2(H) byly velmi stálé bez výraznější sezónnosti, což naznačuje významnou nezávislost na povrchových podmínkách a původ CO2 formujícího skapové vody ne v půdě, ale hlouběji v krasovém profilu. Anomální skap představuje významnou výjimku s výrazně nižším a variabilnějším PCO2(W) a hodnotami PCO2(H) blízko PCO2(air) což indikuje předběžné odplyňování CO2 a srážení kalcitu. Vývoj skapových vod při odplyňování CO2 a následném srážení kalcitu je ilustrován na modelu. Ten ukazuje, že data odpovídající běžným skapům leží podél linie odplyňování se směrnicí ~ −1 směřující k unikátní hodnotě PCO2(H) nezávisle na ročním období. Dále je vidět, že data z anomálního skapu jsou mnohem rozptýlenější a odhad hodnoty PCO2(H) z linie odplyňování bude nesprávně ukazovat hodnotu danou předchozím srážením kalcitu, které změnilo hydrogeochemické vlastnosti vody. Pří větších kulturních akcích v jeskyni Výpůstek byl posuzován také vliv antropogenního CO2 na krasové vody a možnost jejich konverze na vody agre- sivní ke kalcitu. Modelování ukazuje, že při delších akcích se zvýšenou návštěv- ností (500 lidí) je možné dosáhnout dostatečných koncentrací CO2. Běžné pro- hlídky (50 lidí, cca 0,5 h) se zdají být bez významnějšího účinku. Byl sestaven model teoretického vývoje poměru Mg/Ca v průběhu roz- pouštění vápence v epikrasu (T = 10 °C, logPCO2 = –1,5). Vápenec byl simulován jako mix různých poměrů dolomitu a hořečnatého kalcitu. Byly pozorovány dvě zřetelné fáze rozpouštění: (a) počáteční období stechiometrického uvolňování Ca a Mg (kongruentní rozpouštění) a (b) pokročilá fáze začínající po dosažení nasycení kalcitem, odpovídající nekongruentnímu rozpouštění s nárůstem Mg a souběžným poklesem Ca díky srážení kalcitu. Celkový vývoj poměru Mg/Ca reprezentovaný tvarem reakční cesty rozpouštění je dán složením Mg-kalcitu a poměrem Mg-kalcitu a dolomitu (D/C). Při poměrech D/C menších než 1 domi- nuje dynamika rozpouštění Mg-kalcitu a reakční cesty pro vápenec se překrý- vají s cestami pro čistý Mg-kalcit. Další faktor v omezené míře určující tvar reakční cesty je poměr povrchu vápence ku objemu vody ({L}/V). Poměr {L}/V nicméně určuje celkovou dynamiku interakcí. V epikrasu probíhá rozpouštění za podmínek systému otevřeného vůči CO2, což vede k intenzivnějšímu roz- pouštění. Navzdory tomu, že v puklinách hlouběji ve vadózní zóně budou hod- noty {L}/V vyšší, rozpouštění je limitováno uzavřením systému vůči CO2. Oka- mžitý poměr Mg/Ca ve vodě je dán tím, jak daleko podél reakční cesty postou- pilo rozpouštění, což je určeno dobou zadržení vody (tj. doba interakce voda- hornina). Modelové výsledky byly srovnány se skapovými vodami z Punkevních jeskyní i dalších lokalit ve světě. Většina skapových vod vykazovala poměry Mg/Ca odpovídající reakčním cestám Mg-kalcitu s malou příměsí dolomitu, za- tímco vody z dolomitických hornin se podobaly vývoji rozpouštění čistého dolo- mitu. Jak rozdíl v poměrech Mg/Ca mezi anomálním a běžnými skapy, tak mo- del vývoje poměrů Mg/Ca ukazují na důležitost cest proudění vody v krasu pro vznik skapových vod. Nicméně cesty proudění se mohou měnit v čase i prostoru tak, jak se vyvíjí krasový systém, a to nezávisle na klimatických podmínkách. Navíc se zdá, že složení Mg-kalcitu a dynamika rozpouštění hrají podstatnou roli ve vývoji Mg/Ca poměrů ve vápencích. A proto je při paleoenvironmentál- ních studiích krasových proxy dat důležité brát tyto vlivy v potaz. © Pavel Pracný, Masaryk University, 2017 © John Wiley & Sons, Ltd., 2015 © Springer-Verlag Berlin Heidelberg, 2015 © Springer Science+Business Media Dordrecht, 2017 © 2016 Elsevier GmbH., 2017 DECLARATION
I declare that this PhD dissertation is an original report of my reseach and was composed by myself. I confirm that the work submitted is my own, except where work which has formed part of jointly-authored publications has been included. I confirm that appropriate credit has been given within this dissertation where reference has been made to the work of others. I agree that my dissertation may be available in the library of Masaryk Uni- versity.
Brno, 12. 5. 2017 Pavel Pracný TABLE OF CONTENTS
1 Foreword ...... 13 2 Review of dripwater formation and processes ...... 17 2.1 Karst Hydrogeology ...... 17
2.1.1 Dripwater hydrology ...... 19
2.1.2 Residence times ...... 21
2.2 Hydrogeochemistry of dripwaters ...... 21
2.2.1 Carbonate system ...... 21
2.2.2 Dissolution kinetics of carbonate minerals ...... 24
2.2.3 Karst water properties and evolution ...... 25
2.3 Paleoenvironmental reconstructions ...... 28
3 Results and discussion ...... 29 3.1 Hydrogeochemistry of dripwaters in Moravian Karst ...... 29
3.2 Cave and epikarstic PCO2 ...... 32
3.3 Modeling limestone dissolution and Mg/Ca evolution in epikarst ...... 35
3.4 Anthropogenic CO2 influence on dripwaters and speleothem corrosion 38
4 Conclusions ...... 41 5 References ...... 44 Appendix 1 ...... 54 Appendix 2 ...... 70 Appendix 3 ...... 83 Appendix 4 ...... 105
12 1 FOREWORD
Limestone karst areas form large portion of continental crust surface. Be- sides being an important reservoir in biogeochemical cycle of carbon, they pro- vide various benefits for humankind. Many islands are formed of limestone rocks, karst aquifers constitute major fresh water source and karst regions pro- vide plentiful recreational opportunities – from marvelous beaches of southeast Asia, complex cave system of Mammoth Cave in North America to picturesque landscapes of Mediterranean, to name a few. With increasing awareness of climate changes over the last decades, the importance of various paleoclimatic data archives rapidly increased. Karst sys- tems offer important high-resolution terrestrial archives via both the surface (tufas) and the underground (speleothems) secondary karst sediments. Espe- cially speleothems, conveniently sheltered in mostly inaccessible caves, pre- serve detailed record of environmental conditions at the time of their for- mation. In order to recover this information, we have to understand mecha- nisms of speleothem growth as well as the processes, which formed karst wa- ters precipitating speleothems. Although seemingly elementary in its composition, the carbonate karst environment exhibits intriguing levels of complexity. Phenomena considered in holistic approach to karst include biological activity (esp. CO2 production in soil), processes on water-atmosphere boundary (e.g., CO2 degassing in cave), water-rock interaction during limestone dissolution, petrological and miner- alogical characteristics of bedrock, hydrogeological properties of vadose zone among many other aspects. Such intricacy allows for many climatic dependent variables to be reflected in speleothems, while, on the other hand, a careful separation of paleoenvironmental signals is required. Arguably, the most uti- lized paleoenvironmental proxies are stable isotopes, albeit other parameters, esp. trace elements or growth laminae, provide valuable support for reconstruc- tions. Despite utilization of trace elements in reconstructions, the exact mech- anisms of their transfer into the original karst water are not very well under- stood.
13 Additionally, speleothem destruction repeatedly emerges as a lively dis- cussed topic among environmentalists of all kinds (Baker & Genty 1998; Auler & Smart 2004; Dreybrodt et al. 2005; Martín-Peréz et al. 2012). Although pre- vious studies found solid evidence of structural damage causing the destruction in Moravian Karst (Faimon et al. 2004), the possibility of corrosive effects of dripwater cannot still be ruled out. The systematic long-term research in the Punkva Caves might produce decisive insights. This thesis describes results of a research focused on hydrogeochemistry of dripwaters in Moravian Karst. It was particularly aimed to (1) improve our understanding of processes determining the dripwater parameters, possibly applicable as paleoenvironmental proxies (especially minor (Mg2+) and trace (Sr2+) cations) and to (2) find and characterize possibly corrosive dripwaters. The basis for all considerations is a dataset collected in the Punkva Caves in Moravian Karst from February 2012 to March 2013. Although previous works described dripwater properties in some Moravian Karst caves (Faimon et al. 2004; Faimon et al. 2012), our dataset represents an unprecedentedly detailed and thorough study of one site. From hydrogeochemical properties, the relatively low-saturated anoma- lous drip in the Punkva Caves was identified. Based on differences between the anomalous and the regular drips, we could infer possible implications for pale- oenvironmental reconstructions regarding Mg/Ca ratios. Moreover, the de- tailed dataset allowed a comparison of seasonal variations of cave and drip- water CO2 in addition to reconstruction of a hypothetical epikarst PCO2 showing low seasonality. Finally, we devised a kinetic model of limestone dissolution in epikarst, which provided valuable insights into the Mg/Ca evolution during in- congruent dissolution of carbonate minerals. This PhD dissertation presents (a) an introductory review of karst hydro- geology and hydrogeochemistry with regards to cave dripwater and (b) an over- view of the journals’ papers where I have participated as the main author. Other publications related to this topic with my minor contribution or confer- ence presentations are not included in the thesis, but are referenced in the re- view part.
14 The initial impulse which started the whole research endeavor occurred in early 2011, when Ing. Luděk Kabelka, an analytical chemist whose curiosity was fortuitously triggered by a newspaper article about speleothem destruc- tion, contacted Dr. Jiří Faimon, an associate professor at Department of Geo- logical Sciences, Masaryk University, to discuss some insights into the prob- lem. The result was a serendipitous scientific disagreement, which was later drafted into a dissertation topic and offered to me. Dr. Jiří Faimon became my PhD advisor and I am deeply thankful for this opportunity – he taught me many principles of scientific thought, work and writing, spent numerous hours discussing (and on his part mostly explaining) various geochemical topics and helped me to tame my somewhat chaotic and spontaneous nature in order to become more organized and precise – to name a few things I learned from him. In addition, without assistance from Ing. Luděk Kabelka from GEOtest s.r.o., who participated in large portion of the field campaign, provided essen- tial analyses in hydrogeochemical laboratories of GEOtest s.r.o. and negotiated generous donation, the planned intensity of sampling would be impossible to achieve and process. I would like to thank my colleagues and students from Department of Ge- ological Sciences, who either helped me directly with my research or with whom I had the opportunity to cooperate on other research projects, advancing my skills. Namely to Dr. Marek Lang, Ms. Miluška Hradská, Mr. Tomáš Praj, Ms. Klára Blažková, Mrs. Radka Bodláková, Mr. Erik Rzepiel, Dr. Dalibor Všianský and Dr. Josef Zeman, even though there are many more involved. I am also very thankful to Mr. Pavel Kadlec who helped me with the laboratory work and to prof. Ondra Sracek from Palacký University, who participated in isotopic research of dripwaters. The bulk of field work was conducted in the Punkva Caves and would not be possible without kind cooperation of the Cave Administration of the Czech Republic, namely Mr. Hynek Pavelka and Mr. Jiří Hebelka, who also provided
15 climatic data from administration’s meteo-station, and additional cooperation from the Administration of the Protected Landscape Area Moravský kras. At last but not least, I am wholeheartedly thankful to my whole family for lasting support through my study years and their confidence in happy ending. And finally, I would like to state with profound gratitude that I deeply admire my wife Simona for her immense kindness, patience and encouragement in face of my PhD adventures. Thank you!
16 2 REVIEW OF DRIPWATER FORMATION AND PROCESSES
2.1 KARST HYDROGEOLOGY The most important karst rock is limestone, a rock composed of carbonate minerals, mostly calcite (CaCO3) and dolomite (CaMg(CO3)2) and to a lesser extent of aragonite (CaCO3). A significant component might be Mg-calcite which contains up to 20% of magnesium, although it is less stable than calcite and tends to re-crystallize during diagenesis (Mackenzie et al. 1983; Ford & Williams 2007). Geomorphologically, limestone karst in temperate climate re- gions is composed of wide plateaus intersected by deep valleys with various surface features formed by water (sinkholes, polje, limestone pavement etc.). One of the primary sources of water flowing through karst profile is infil- trated precipitation. The infiltrating water moves (a) by a slow matrix flow through interstitial spaces between soil grains and/or (b) via preferential flow along the macropores, e.g., animal dens/holes, root residue, or mud cracks (Kogovšek & Šebela 2004). Although soil thickness significantly increases wa- ter storage capacity, water percolation into caves is commonly observed even in regions without any soil coverage (Klimchouk 2004). Another substantial source of karst water is infiltrating water flowing from non-karstic terrains into carbonate rocks (Ford & Williams 2007). Due to the high solubility of limestone, any surface water in karst areas rapidly disappears under surface shaping characteristic underground cave phenomena. Therefore, the conditions in karstic aquifers are largely deter- mined by the geological and lithological properties of given karst region. Per some authors (e.g., Ternan 1972; Trček 2003; Smart & Worthington 2004) the underground environment can be divided into two hydrogeological structures: A. A network of interconnected karst conduits with high overall permeabil- ity, fracture porosity and quick preferential flow along these conduits. The water residence time is short. Therefore, the structure serves as a drainage of less permeable rock. B. Rock massive with low permeability, i.e., limestone bedrock with pri- mary porosity and secondary tectonic permeability. In spite of the low
17 permeability, the flow velocity is significantly lower and the residence time is usually very long. Another division of karst underground (Palmer 2006; Ford & Williams 2007) is a zonation based on the hydrogeological properties and describes it in terms of three major zones: 1. The uppermost unsaturated zone of strongly weathered rock underneath soil – epikarst. 2. Lower unsaturated zone (vadose zone) with large underground caverns. The lowest part of vadose zone can be periodically saturated (e.g., during floods) forming an epiphreatic (sub)zone. 3. Saturated zone (phreatic zone), the zone of permanent water saturation under the water table where water flows primarily horizontally via con- duit systems. As the vertical permeability in epikarst zone rapidly diminishes on the boundary with vadose zone, excessive water is stored in pores, fractures and joints forming perched aquifers. From the aquifers water flows along fissures downward into the vadose zone and is eventually released into the cave as drip- water (Perrin et al. 2003a; Ford & Williams 2007; Williams 2008; Jones 2013). Therefore, epikarst is of major hydrologic importance. Its overall water capac- ity is given by i. epikarst zone thickness, ii. porosity, iii. water inflow and outflow balance. Whereas the thickness and porosity determine the sum of space available for the water, the ability to retain water is given by the water inflow/outflow balance. The outflow is in addition to hydrostatic pressure determined by the vertical hydraulic conductivity of the vadose zone, which is dependent on frac- ture porosity and can be very variable. Therefore, the epikarst water capacity might significantly vary even on small scale.
18 2.1.1 Dripwater hydrology Rudimental dripwater hydrology studies are based on discharge response to precipitation (Smart & Friederich 1987; Baker et al. 1997; Genty & Deflan- dre 1998; McDonald & Drysdale 2007), but the research focus is gradually shifted towards analyses of relation between precipitation and geochemical properties of dripwater (Perrin et al. 2003a; Musgrove & Banner 2004; Cruz et al. 2005; Schwarz et al. 2009; Riechelmann et al. 2011; Faimon et al. 2012; Kamas et al. 2015). A few modes of water flow were distinguished in karst un- saturated zone and are consistent with the division described in previous chap- ter: water flows (A) via a system of interconnected fissures, joints and conduits (conduit flow) and/or (B) seeps through the bedrock matrix (matrix flow). These flow types may differ in residence times and hydrogeochemical properties of resulting dripwaters. A classification of karst underground water flow was based on the discharge and discharge variability (Friederich and Smart 1982), which was latter improved by Baker et al. (1997) and is commonly utilized (e.g., Spötl et al. 2005; Baldini et al. 2006a; McDonald & Drysdale 2007; Hartland et al. 2012; Fairchild & Baker 2012; Pracný et al. 2016). Dripwaters show various responses to atmospheric precipitation leading to linear, non-linear, or even no direct correlation (e.g., McDonald et al. 2007; Pronk et al. 2009; Riechelmann et al. 2011; Faimon et al. 2016). The drips with stable discharge during dry periods without atmospheric precipitation might indicate a water source in slowly percolating infiltration or very large perched aquifer. On the contrary, the drips with rapid response to precipitation may have faster connection with the surface. Based on a long-term monitoring Tooth and Fairchild (2003) described four types of precipitation responses: 1. Rapid response without time-lag, after which the discharge slowly de- creases. 2. Rapid response with associated time-lag, after which the discharge slowly decreases. 3. Intermittent response, when a particular threshold of water input must be exceeded, before the drip discharge responses, occasionally even by decreasing.
19 4. No response to particular precipitation; the changes in discharge inten- sity are not correlated with precipitation. The non-linear responses to precipitation are explained by an aquifer with an overflow – it maintains a constant water head and in case of intensive influx the excessive water is drained via otherwise dry outlet (and could feed a sea- sonal drip). Tracer studies (Kogovšek & Šebela 2004; Goldscheider et al. 2008; Pronk et al. 2009; Kogovsek & Petric 2014) showed that the tracer concentration in dripwater diminished exponentially suggesting a simple linear flow. Neverthe- less, in case of some drips, the concentration increased again after another rain- fall event indicating existence of periodically drained primary collectors. Be- sides, significant lateral dispersion during the wet periods was observed in some dripwater tracer experiments (Bottrell & Atkinson 1992). Another pro- cess considerably affecting tracer concentrations is a mixing along the water flow path. It seems that hydraulically the epikarst is composed of numerous flow paths that can be in limited contact and allow partial water mixing (Perrin et al. 2003a; Perrin et al. 2007). The observed lateral dispersion proves hori- zontal spread of water during wet periods, indicating overflow from the pri- mary aquifer into adjacent free spaces enabling increased mixing and feeding of seasonal drips. Conversely, the epikarst aquifer becomes fractionated into separate reservoirs feeding individual drips (or groups of drips) during dry pe- riods. An insightful study by Genty and Deflandre (1998) showed that the vol- ume of drop in dripwater was ~0.15 mL in 94% cases. Nevertheless, the volume varies under extreme discharge. At higher discharges (over ~50 drops/min), the drop volume decreases, possibly due to limitations of fluid physical properties (density, surface tension, water pressure connected to the flow) or development of side-drips. On the contrary, the volume increases at lower discharges (under 1 drop/min). In addition, Fairchild et al. (2006b) report that drop volume for drips with various discharges was in range from 0.146 to 0.154 mL.
20 2.1.2 Residence times An important mechanism participating in dripwater hydrogeochemistry could be the piston flow effect. Piston flow is a water transport mechanism oc- curring when the infiltrating water pushes the water retained in flow paths since previous precipitation event towards the outlets. In this case, water dis- charge in cave increases almost instantly after major precipitation events. In addition, the water pushed by the piston flow is often flushed out from parts of the reservoir with longer residence times as indicated by the increased trace element content and isotope ratios (e.g., Aquilina et al. 2006, Tooth & Fairchild 2003, Emblanch et al. 2003). The overall residence times of water in vadose zone are quite variable and ranging from days (Bottrell & Atkinson 1992; Genty & Deflandre 1998; Perrin et al. 2003b; Kogovšek & Šebela 2004; Kamas et al. 2015; Faimon et al. 2016) to months and even years (Spötl et al. 2005; Aquilina et al. 2006; Kluge et al. 2010; Kogovšek & Petric 2014) as a result of very complicated water flow paths or extremely slow matrix flow. During the wet periods, a continuous flow through all types of hydraulically connected fissures is enabled. Therefore, a relatively fast transport of solution from surface into cave occurs. In contrast, during the dry periods, only a portion of flow paths is active and the major part of infiltrated precipitation is stored in less permeable parts of vadose zone. This water might be flushed out after intense precipitation events (Kogovšek & Šebela 2004; Kogovšek & Petric 2014).
2.2 HYDROGEOCHEMISTRY OF DRIPWATERS The hydrogeochemical properties of dripwaters are given by processes oc- curring along the water flow path from karst surface into a cave and are ulti- mately embedded into the speleothem structure and composition.
2.2.1 Carbonate system Biogeochemically the most important compound of carbon is carbon diox- ide CO2. The solubility of gaseous CO2 in water decreases with rising tempera- ture. The equilibrium between CO2 in atmosphere above the solution and CO2 in the solution is expressed by Henry’s constant (Sander 2015):
21 c Hcp = CO2(푎푞) (1) PCO2
Where cCO2(aq) is the concentration of CO2 in the aqueous phase and PCO2 is the partial pressure of CO2 in the atmosphere. Oftentimes, it is defined as a dimensionless parameter given by the ratio between the aqueous phase con- centration (cCO2(aq)) of a specie and its gas-phase concentration (cCO2(g)) c Hcc = CO2(푎푞) (2) cCO2(푔) The values for standard conditions are presented in Table 1. Dissolved carbon dioxide reacts with water to form carbonic acid (Stumm & Morgan 2012):
CO2(푎푞) + H2O = H2CO3 (3) with the equilibrium constant
푎퐻2퐶푂3(푎푞) K0 = (4) 푎퐶푂2(푎푞)
Although unhydrated CO2(aq) is much more abundant than H2CO3, a con- vention was adopted to include the hydration in the overall dissolution in order to facilitate calculations as the hydration step is effectively instantaneous. The two carbonate species are summarily expressed as H2CO3* and the overall re- action of CO2 dissolution becomes ∗ CO2(푔) + H2O = H2CO3 (5) Carbonic acid further dissociates into two species
* - + H2CO3 = HCO3 + H (6) and − 2− + HCO3 = CO3 + H (7)
with equilibrium constants K1 and K2 which are presented in Table 1 (and
K1 is in fact composite constant including the CO2(aq) hydration). The distribu- tion of carbonate species is determined by pH value. Under acidic conditions of pH < 4 the system contains effectively only carbonic acid. As the pH increases, the acid dissociates into bicarbonate reaching maximum at pH = 8.3. With fur- ther pH increase, the bicarbonate ion dissociates into the carbonate ion, which effectively dominates the carbonate system when pH is above 12. The calcium carbonate (calcite or aragonite) dissolution can be described by simple dissolution equation:
22 2− 2+ CaCO3(s) = CO3 + Ca (8) with solubility products for calcite (Kc) and aragonite (Ka)
K = 푎 2+ 푎 2− (9) 푐/푎 Ca CO3 In karst environment studies, the pinnacle of interest lies in interactions of carbonate minerals with water and carbon dioxide. The dissolution under open system conditions can be expressed by equation − 2+ CaCO3(푠) + CO2(푔) + H2O = 2HCO3 + Ca (10) and total equilibrium constant
2 푎 − 푎 2+ HCO3 Ca KT = (11) PCO2 The value of calcite solubility product is temperature dependent and is also affected by variations in Ca activity caused by presence of aqueous com- plexes such as CaCO30 and CaHCO3+ (Jacobson & Langmuir 1974). Dissolution of dolomite in pure water is analogically described as 2− 2+ 2+ CaMg(CO3)2(푠) = 2CO3 + Ca + Mg (12) with solubility product 2 K푑 = 푎 2+푎 2+푎 2− (13) Ca Mg CO3
The equation for open system in contact with CO2 and dolomite is − 2+ 2+ CaMg(CO3)2(푠) + 2CO2(푔) + 2H2O = 4HCO3 + Ca + Mg (14) Dissolution of Mg-calcite with content of Ca = x and Mg = y where x + y = 1 is given by equation: 2− 2+ 2+ Ca푥Mg푦CO3(s) = CO3 + 푥Ca + 푦Mg (15) with solubility product 푥 푦 K = 푎 2+ 푎 2+ 푎 2− (16) 푀푔−푐푎푙푐푖푡푒 Ca Mg CO3 Due to various sources, genesis and composition of Mg-calcite, the pub- lished solubility products vary in a wide range from ca. 10−8.5 to 10−7.4 (Plummer & Mackenzie 1974; Mackenzie et al. 1983; Morse & Mackenzie 1990).
Solubility product Ksp of a mineral is defined by the equilibrium activities of minerals’ constituents. If the system is not in equilibrium the immediate species’ activities constitute ion activity product IAP. By comparing Ksp to IAP the reaction direction can be assumed – if IAP/Ksp > 1 the reaction will proceed towards the reactants and vice versa. This relation is useful to specify solution
23 saturation with respect to a mineral. The most common expression is the satu- ration index SI = log IAP/Ksp. Supersaturated solutions show positive SI values, whereas not saturated solutions show negative values.
Table 1 Equilibrium constants for the carbonate system − log T [°C] 0 5 10 15 20 25
K0 1.11 1.19 1.27 1.34 1.41 1.47
K1 6.58 6.52 6.46 6.42 6.38 6.35
K2 10.63 10.55 10.49 10.43 10.38 10.33
Kc 8.38 8.39 8.41 8.34 8.45 8.48
Ka 8.22 8.24 8.26 8.28 8.31 8.34
Kd - - - - - 17.2±2 All values from Plummer & Busenberg (1982) except dolomite (Sherman & Barak 2000)
2.2.2 Dissolution kinetics of carbonate minerals Dissolution and precipitation rates of carbonate minerals are determined by various factors and were studied under wide range of conditions with signif- icant variations in resulting rate values (see Morse & Arvidson 2002 for review; Arvidson et al. 2003; Kaufmann & Dreybrodt 2007; Morse et al. 2007; Cubillas et al. 2005; Pokrovsky et al. 2005; Pokrovsky et al. 2009). The basic mechanisms of dissolution on mineral surface were described by Plummer et al. (1978) and expanded by Chou et al. (1989) as follows for calcite:
kc1 HCaCO Ca 2 HCO (17) 3 3 kc1
kc 2 HCaCO CO * Ca 2 2HCO (18) 323 3 kc 2
kc3 CaCO 2 COCa 2 (19) 3 3 kc3 and for dolomite:
kd1 H2)CaMg(CO Ca 2 Mg 2 2H CO (20) 23 3 kd1
k d 2 H2)CaMg(CO CO * Ca 2 Mg 2 4H CO (21) 3223 3 kd 2
24 k d 3 )CaMg(CO Ca 2 Mg 2 2CO 2 (22) 23 3 kd 3 The total forward and backward rates for calcite are than expressed as:
+ ∗ (23) 푅푓(푐푎푙푐푖푡푒) = 푘푐1푎H + 푘푐2푎H2CO3 + 푘푐3
푅 = 푘 푎 2+푎 − + 푘 푎 2+푎 − + 푘 푎 2+푎 2− (24) 푏(푐푎푙푐푖푡푒) −푐1 Ca HCO3 −푐2 Ca HCO3 −푐3 Ca CO3 The forward rate equations for calcite and dolomite differ only in reaction order. Whereas for calcite the reaction order n = 1, experiments show that the reaction order of dolomite dissolution with respect to aH+ is a fractional number. Busenberg and Plummer (1982) found that n = 0.5 (for 25 °C) and the value increases with increasing temperature. Definition of the dolomite backward rate is more complicated. Busenberg & Plummer (1982) determined that the
CO32− is not responsible for the observed backward reaction and it is independ- ent on the activities of Ca2+ and Mg2+ below pH 6 and far from equilibrium, resulting in
− (25) 푅푏(푑표푙표푚푖푡푒) = 푅푓(푑표푙표푚푖푡푒) − 푅푡표푡푎푙 = 푘4푎HCO3
Where k4 is the rate constant of backward reaction. Generally speaking, the dissolution mechanisms of calcite are transferra- ble to other mono-cation carbonate minerals (e.g., magnesite and aragonite) whereas dissolution of dolomite is an example of composite carbonate mineral. The expression for dissolution rate of Mg-calcite is expanded to include magne- sium activity and stoichiometry of Ca2+ and Mg2+.
2.2.3 Karst water properties and evolution Hydrochemical properties of karst water are defined by the content of sub- stances in solution, which are incorporated especially via dissolution of miner- als in contact with infiltrating water and gaseous CO2. The initial composition of infiltrating water is result of complex interactions in atmosphere. The chem- ical constituents dissolved in rainwater are usually connected to source oceanic water and particle/gases in atmosphere. Local sources of pollution in urban ar- eas may participate. They are utilized as an indicator of the extent of anthro- pogenic pollution. Among the most important pollutants emitted by humans are both nitrogen and sulfur oxides, which are converted into the acids that are
25 a major cause of water acidity. In turn, acidic rainwater outwashes heavy met- als present in atmosphere mainly from industrial sources and road transporta- tion emissions (e.g. Eriksson 1952; Carrol 1962; Gatz 1991; Paternoster et al. 2014; Vet et al. 2014). The seawater contribution to precipitation composition is ordinarily determined by comparing ion ratios to the same ratios in marine water (e.g. D'Alessandro et al. 2013). The dominant dissolved ions in carbonate karst waters are calcium and carbonates released from calcite dissolution. Other cations released by dissolu- tion of natural calcite/limestone are Mg2+, Fe2+, Mn2+ and Sr2+. One of the de- termining constituents of karst water is dissolved carbon dioxide and its’ spe- cies (Atkinson 1977). The additional anions (e.g. Cl−; SO42− or NO3−) are either originating from minor minerals present in limestone and soil, from anthropo- genic pollution, or are initially present in precipitated water (e.g., Perrin et al. 2003b). Formation of dripwater begins immediately after infiltration of atmos- pheric precipitation. Infiltrating water dissolves CO2 present in soil/epikarst atmosphere that participates on carbonate mineral dissolution (see equation
3). The source of CO2 is linked to biological activity – e.g., plant roots respira- tion or microbial decomposition of organic matter (Kuzyakov 2006). CO2 con- centrations in soils are climatically and seasonally dependent (esp. in temper- ate climate); they are enhanced in warm and wet regions and seasons (Sanchez- Canete et al. 2011; Plestenjak et al. 2012). Albeit soil is widely accepted as a main source of CO2, enhanced concentrations (2–6 vol%) were measured in va- dose zone of karst in Mediterranean area (Benavente et al. 2010). In karst soils of Moravian Karst, the concentrations usually reach up to 1 vol% with strong seasonal fluctuations leading to highest values in summer and lowest in winter (Faimon & Ličbinská 2010; Blecha & Faimon 2014). The acidic water in epikarstic perched aquifer dissolves limestone under conditions (esp. PCO2) which can be retrospectively estimated from dripwater composition (e.g. Faimon et al. 2012; Peyraube et al. 2013; Milanolo & Gabrovšek 2015; Pracný et al. 2016b). Eventually, the water is drained into
26 joints or fissures and flows downwards into the cave system. The water resi- dence times of perennial drips are usually long enough for the water to be fully equilibrated with respect to calcite (Spötl et al. 2005; Kogovšek & Petric 2014). Then, the residence time usually does not affect the concentrations of major constituents of dripwater. Nevertheless, the limestone bedrock might contain numerous minerals with slower dissolution kinetics (e.g. dolomite or clay min- erals) that release trace elements, which are important in paleoenvironmental analysis of speleothems. One of processes affecting hydrogeochemical properties of dripwater is prior calcite precipitation (PCP, Fairchild et al. 2000). It occurs whenever the aqueous CO2 is able to degas into ventilated spaces within the vadose zone be- tween epikarst aquifer and drip site. During calcite precipitation, aqueous Ca concentration decreases and, thus, the trace element to calcium ratios in water increase, possibly archived in forming speleothems. PCP is believed to be pro- moted by drier climatic conditions (Fairchild et al. 2000, McMillan et al. 2005, Fairchild et al. 2006b, Tremaine & Froehlich 2013). Another process influencing composition of dripwater is mixing along the water flow paths. Waters formed in contact with different PCO2 and distinct bedrock composition, or under various surface conditions (forestation, anthro- pogenic pollution etc.) are mixed in karst profile above the cave (Perrin et al. 2003a; Perrin et al. 2007; Moore et al. 2009; Schwarz et al. 2009; Gabrovšek & Dreybrodt 2010). After reaching the cave, a drip is formed and the dripwater degasses the excess of dissolved CO2. The driving force of the degassing is the difference be- tween CO2 present in water and CO2 in cave air. Because aqueous concentra- tion of CO2 corresponds to a specific value of PCO2 in contact with the solution (see equation 3), the partial pressure is commonly used to express the content of CO2 in water and referenced to as a CO2 partial pressure in water (PCO2(W)). Therefore, the driving force of dripwater degassing can be conveniently de- scribed as the difference between PCO2(W) and PCO2(atmosphere). As the CO2 is re- leased from the water, the driving force diminishes until the partial pressures equalize. Dripwater degassing carbon dioxide is an important source of CO2 in
27 cave air (Bourges et al. 2001; Baldini et al. 2008). As the dripwater degasses, the carbonate equilibria shift and the water becomes supersaturated with re- spect to calcite, which is followed by calcite precipitation and speleothem for- mation.
2.3 PALEOENVIRONMENTAL RECONSTRUCTIONS Cave dripwaters and their precipitates (speleothems) are studied world- wide as a source of information about environment in geological past (see McDermott 2004 or Fairchild et al. 2006a for review). As many of processes participating on water composition are climatically controlled, speleothem com- position can be used to reconstruct paleoenvironmental conditions during its formation (Li et al. 2005; Griffiths et al. 2010; Borsato et al. 2016; Paar et al. 2016). Most common technique is stable isotope analysis (McDermott 2004; Drysdale et al. 2005; Verheyden et al. 2008; Lachniet 2009), nevertheless also trace elements can provide useful data (e.g., Verheyden et al. 2000; Huang et al. 2001; McMillan et al. 2005; Fairchild et al. 2006b; Cruz et al. 2007; Wong et al. 2011; Frisia et al. 2012; Jochum et al. 2012; Sinclair et al. 2012; Meyer et al. 2014; Tan et al. 2014; Orland et al. 2014; Casteel & Banner 2015; Bernal et al. 2016). The enhanced trace element ratio is interpreted as an effect of PCP or longer residence time in arid climate, and therefore is frequently tested as a paleoenvironmental proxy (Verheyden et al. 2000; Fairchild et al. 2000; Fairchild & McMillan 2007; Fairchild & Treble 2009; Tremaine & Froehlich 2013). In fact, Mg/Ca ratios may sometimes paradoxically show positive corre- lation with rainfall (e.g., Baldini et al. 2012). In the last few years, more thor- ough studies indicate that dripwater properties are climate sensitive to lesser extent than previously anticipated (e.g., Baker et al. 2016). Furthermore, the variations in dripwater hydrogeochemistry are not necessarily linked with changes in drip discharge (Musgrave & Banner 2004; Faimon et al. 2016) or the response can be non-linear (Karmann et al. 2007). Additionally, it seems that changes in calcite fabrics and crystal habits might be also a record of sat- uration and water discharge variability (e.g. Genty & Quinif 1996; Frisia et al. 2000; Niggemann et al. 2003 or Riechelmann et al. 2014).
28 3 RESULTS AND DISCUSSION
3.1 HYDROGEOCHEMISTRY OF DRIPWATERS IN MORAVIAN KARST Hydrogeochemical properties of dripwaters directly determine speleo- them growth and composition. In order to study differences between dripwaters in a cave system, a long term monitoring of dripwater properties and cave en- vironment was realized in Punkva Caves in Moravian Karst. Sampling sites were located in a corridor behind Přední Dome (drips PC), in Zadní Dome (drip ZD) and in Tunnel Corridor (drips TC). Collected data were furthermore com- pared with older dripwater research on other sites in Moravian Karst to assess possible similarities/differences between particular caves in one karst region. The study was realized from February 2012 to March 2013 and 126 samples were collected during 26 sampling events (twice per month) and compared with 45 analyses from the archive dataset. Meteorological data were provided by the Cave Administration of the Czech Republic from a station situated above the cave and operated by the Czech Meteorological Institute. The sampled drip- waters were partially analyzed directly on site (volumetric determination of alkalinity and Ca) and in a laboratory (ICP-OES analysis of Mg and Sr). Drip- waters were also twice (April and November) sampled for stable isotope analy- sis (LAS; δ18O, δ2H and δ13C). For detailed description of the site and used methods, see Appendix 1 and Appendix 2. Studied drips showed different flow regimes regardless of drip location – drip CP2 had very stable discharge (variation coefficient 17.6%), whereas drips CP3 and TC1 showed higher variation (v.c. ~50%) and drips CP1, TC2 and ZD showed very wide variations. The drip TC2 eventually even stopped dripping. With exception of the stable drip CP2, the decreasing discharges indicate, that the epikarstic perched aquifers feeding the drips were gradually emptied from May to December. Interestingly, no drips showed response to summer storm events, which might be caused by very wide measurement step. Another con- tributing factor might be evapotranspiration as shown by a model (see Appen- dix 1: Fig. 6). The discharges started to rise again in mid-December and rose
29 until the end of measurements in March. This development indicates a sub- stantial replenishment of epikarst aquifers from snowmelt water. Neverthe- less, drip CP3 did not increase discharge and it seems that the perched aquifer is replenished under specific conditions. In addition, a delay of 3–4 months be- tween modeled infiltration and drip discharges suggests that the water is sea- sonally drained through some preferential pathways. In such case, it would not significantly contribute to epikarst aquifers feeding the studied drips. The studied dripwaters were divided into two groups according to their hydrogeochemical properties. On one side, there is the anomalous drip TC1, on the other side the rest of the drips referred to as the regular drips. These drips showed higher values of EC, significant supersaturation with respect to calcite and low trace elements ratios. Isotopically, the drips were very similar to each other with no signs of additional processes. They demonstrated slight enrich- ment in δ18O and δ2H in November samples and δ13C values corresponding to calcite dissolution under closed system conditions. In contrast, the anomalous drip showed peculiar properties – systematically lower EC, saturation close to equilibrium with respect to calcite and enhanced Mg/Ca and Sr/Ca ratios com- pared to regular drips. These ratios were caused by significantly lower Ca con- centration. The isotopic composition showed strong enrichment in δ13C indicat- ing water degassing, but the δ18O and δ2H values were similar to values in regular drips. The δ18O and δ2H are close to Global Meteoric Water Line (GMWL) and indicate very limited effect of evaporation and fast precipitation infiltration. What is more, the summer/fall difference in δ18O values indicates incomplete mixing in the epikarst aquifer as the isotopically lighter summer rainwater remains at least partially differentiable. A few possible explanations of anomalous drip properties were proposed and discussed. Firstly (a) an enhanced water dynamic in karst profile – a situ- ation when infiltrating dripwater reaches the cave before it could attain equi- librium with soil/epikarstic CO2 and calcite. It is therefore under-saturated with respect to calcite, shows low mineralization and Mg/Ca ratio and very var- iable discharge with strong reaction to rainfall. However, these properties do not fit the anomalous drip. Another possible explanation is (b) water mixing
30 along the water flow path somewhere in vadose zone. Although water mixing can explain low saturation or under-saturation with respect to calcite and lower mineralization, it does not seem possible – considering Moravian Karst lime- stones composition – that it could provide enhanced trace element ratios in Punkva Caves. In conclusion, an explanation of anomalous properties via (c) prior calcite precipitation (PCP) seems as the most plausible. Effects of PCP on dripwater include lower mineralization and saturation with respect to calcite, enhanced trace element ratios and increase in δ13C, while not necessarily in- fluencing the dripwater hydrology. Even the speleological situation on site sup- ports PCP hypothesis, as the water flows through ca. 20 long crevice with many speleothems before being sampled (Glozar 1984). Considering the possible speleothem corrosion, the anomalous dripwater does not seem to pose threat to cave environment. Its properties and supposed origin via PCP suggest water in equilibrium or slight supersaturation with re- spect to calcite and cave CO2. An unrealistically high PCO2 in cave air would be necessary for the water to become aggressive to calcite as demonstrated in Výpustek – another Moravian Karst cave – where influence of ventilation and increased PCO2 during large cultural events was studied (see chapter 3.4 and Appendix 4). The anomalous dripwater geochemical properties would be imprinted into a speleothem growing from the water, raising numerous questions regarding speleothem utilization as climatic proxies. The trace element incorporation in calcite is interpreted as indicator of temperature changes (Roberts et al. 1999; Huang & Fairchild 2001), dry periods (Fairchild and McMillan 2007) or chang- ing elemental sources (Ayalon 1999). In contrast to general understanding (e.g. Fairchild et al. 2000), the correlation of Mg/Ca to the drip discharge in Punkva Caves does not show negative correlation. Although a several years long da- taset would be necessary to assess the effect of long term arid periods, no intra- seasonal changes were observed. Another problematic utilization of a speleo- them precipitated from anomalous dripwater is evaluation of calcite precipita- tion dynamic from crystal habits (e.g., Frisia et al. 2000; Niggemann et al. 2003
31 or Riechelmann et al. 2014). The final problematic interpretation is the utiliza- tion of δ13C isotopic data. The δ13C is supposed to be determined by soil CO2 composition (Dreybrodt & Scholz 2011), but the anomalous drip shows enrich- ment caused by the PCP. An anomalous dripwater speleothem would therefore incorrectly indicate lower biological activity in soil. The properties of the anomalous drip are of permanent character and sup- posedly controlled spatially, implying that the spatial conditions might be an important factor that should be considered with the temporal conditions (e.g. dry/wet seasons). Moreover, these spatial effects might change over long peri- ods as karst/cave conditions develop (e.g. by weathering and karst evolution). It could be concluded that it would be appropriate to use more speleothems from one cave system in any paleoclimatic reconstruction based on speleothem proxies to eliminate possible distortion by anomalous speleothems.
3.2 CAVE AND EPIKARSTIC PCO2 The amount of limestone dissolved in dripwater is reflected by Ca2+ and
CO32− content that determines solution’s undersaturation or supersaturation with respect to calcite. The total quantity of dissolved calcite is limited by the conditions under which the dissolution occurred, especially PCO2 in soil/epikarst aquifer. When the percolating water enters a cave, it immediately starts to de- gas CO2, the pH rises as well as CO32− content as the carbonate balance shifts and the water becomes supersaturated to calcite. This subsequently leads to calcite precipitation and speleothem growth. The intensity of degassing is de- pendent on the difference between PCO2 equivalent present in the water and
PCO2 in the cave’s air: the bigger the difference, the higher supersaturation can be reached. Based on the hydrogeochemical properties of dripwater, the esti- mate of hypothetical CO2 partial pressure under which the water was formed can be made (Faimon et al. 2012; Peyraube et al. 2012; Milanolo & Gabrovšek 2015; Pracný et al. 2016b). A model based on reconstruction from calcium and carbonate species concentrations and pH was composed to study hypothetical
PCO2 participating on formation of dripwaters in Punkva Caves in Moravian
Karst (Czech Republic) and compare them with PCO2 in the cave’s air and in
32 sampled dripwaters to better understand some of the processes determining speleothem growth. Punkva Cave system is developed in central part of Moravian Karst in Devonian limestones of the Macocha formation. Studied dripwaters are situ- ated in upper dry level of the cave under ca. 100 m thick layer of limestone.
Concentrations of CO2 in Moravian Karst reach up to 1% in soils (Faimon & Ličbinská 2010; Blecha & Faimon 2014) and 1–11% in caves (Otava 1995; Fai- mon et al. 2012). Dripwater and speleoclimatic data were collected during 26 measuring campaigns from February 2012 to March 2013. Meteorological data were obtained from measuring station situated on the surface above the cave. All data were statistically analyzed for presence of significant correlations and presence of cycles. For detailed description of the site and used methods, see Appendix 1 and Appendix 2.
The partial pressure of CO2 in the cave’s air (PCO2(air)) showed significant seasonality with calculated period of about 304 days. The minima were meas- ured in winter (10−3.31, i.e., 0.06 vol%) and the maxima in summer (10−2.49, i.e.,
0.32 vol%). Both the partial pressure of the CO2 corresponding to aqueous car- bonates, PCO2(W), and the partial pressure of CO2 participating on the initial water formation, PCO2(H), were calculated from dripwater hydrogeochemistry.
Whereas the PCO2(W) showed clear seasonal variations with minima in winter (10−2.91, i.e., 0.12 vol%) and maxima in summer (10−2.35, i.e., 0.45 vol%) the
PCO2(H) did not show any significant seasonality. The partial pressure was cal- culated in a narrow range from 10−1.77 to 10−1.49 (1.7–3.2 vol%) indicating only slight dependence on surface conditions.
The cave air CO2 seasonality might be driven by the difference between exterior and interior temperature. Whereas the daily temperature maxima be- low mean annual temperature (MAT) indicate upward airflow ventilation mode (UAF mode, totally 141 days), the daily temperature minima above MAT indi- cate daily downward airflow mode (DAF mode, totally 105 days) in the cave and both modes represent the periods of active cave ventilation (Faimon et al.
2012). In UAF mode, the PCO2(air) is systematically lower, as the air flows out from the cave through upper openings, which contrasts with DAF mode, when
33 the PCO2(air) rises probably by influx from cracks and joints leading through epikarst and soil (Lang et al. 2017). Nevertheless, the correlations of external temperature with stationary PCO2(air) data are not significant, which might be caused by time shifts in monitoring steps. In addition, neither the internal tem- peratures correlate with PCO2(air).
The periodicity in PCO2(W) is less pronounced than periodicity of PCO2(air) and with different periods. The raw data are correlated statistically signifi- cantly, which indicates interconnection via degassing. No significant correla- tion was found in the stationary data, which might be caused by additional effect, e.g., drip rate variations or variations in sampling. If the drip rate is slow, the residence time of an individual hanging drop on the speleothem in- creases and it can degas much more CO2 than under higher discharge. The intensity of degassing is also affected by the difference between the initial
PCO2(W) of water entering cave and PCO2(air) (Faimon et al. 2016).
In contrast to PCO2(air) and PCO2(W), the hypothetical PCO2(H) shows very sta- ble values for all drips except for the previously identified anomalous drip. The stable PCO2(H) values in all other drips (the regular drips) indicate independence on surface conditions and suggest that the CO2 source might be situated in deeper parts of epikarst or vadose zone. Nevertheless, the precise source and presumed biogeochemical conditions are unknown. Moreover, the estimated values of PCO2(H) are higher than concentrations measured in soils (Faimon & Ličbinská 2010; Sanchez-Cañete et al. 2011; Plestenjak et al. 2012; Blecha & Faimon 2014).
On the contrary, the PCO2(H) reconstructed for the anomalous drip shows seasonal variations similar to PCO2(air) and PCO2(W). Statistical correlation as well as the calculated periodicity demonstrate a strong connection. The close rela- tion of the partial pressures indicates existence of a mechanism disrupting the original PCO2(H). A model of geochemical evolution (based on Peyraube et al. 2012 and Milanolo & Gabrovšek 2015) illustrates the difference between the regular drips and the anomalous drip (see Appendix 2). The model shows that regular dripwaters have undergone uniform degassing evolution shifted only by different initial conditions, i.e., by different PCO2(H). In addition, the water
34 from anomalous drip seems to have firstly degassed to equilibrium with PCO2(air) and then precipitated calcite. This sequence is especially plausible considering kinetics of these processes. Therefore, the most probable cause of the anoma- lous properties seems to be prior calcite precipitation in some caverns above sampling site.
3.3 MODELING LIMESTONE DISSOLUTION AND MG/CA EVOLUTION IN EPIKARST The main source of minor (e.g. Mg) and trace (e.g. Sr) elements in drip- waters is limestone dissolution. The mechanism and dynamics of Mg release from carbonates is generally neglected in paleoclimatic studies despite using Mg/Ca ratio as an important proxy. In order to examine this assessment, a model of theoretical evolution of Mg/Ca ratio during limestone dissolution un- der conditions expected in epikarst was composed and compared with drip- water data from Punkva Caves in Moravian Karst and other caves worldwide. The study was focused especially on the effect of Mg-calcite and dolomite dis- solution dynamics. Limestone dissolution was studied using a dynamic model of conjoined dissolution of two separate minerals, Mg-calcite and dolomite, while calcite and magnesite were initially not present in solution, but were allowed to precipitate when the solution became supersaturated with respect to the minerals during simulation. Series of model calculations were run for epikarstic conditions of T
= 10 °C and log PCO2 = −1.5, fixed rock/solution ({L}/V) and water/atmosphere ({S}/V) boundary area and a range of dolomite/Mg-calcite ratios (D/C). For de- tailed description, how the ratios were estimated, see Appendix 3. The model is an open system where CO2 from soil/epikarst air enters the solution. The
CO2 exchange is derived from the two-layer model (Liss & Slater 1974; Stumm & Morgan 1996). The detailed descriptions of Punkva Caves limestones and dripwaters compared to the model solutions are presented in Appendix 3. The model showed that the reaction path of dissolution of composite car- bonates (dolomite or Mg-calcite) initially follows a straight line determined by mineral stoichiometry indicating congruent dissolution. When the solution reaches calcite saturation, the character of dissolution changes to incongruent
35 – this is illustrated in the model by the reaction path becoming non-linear and increasing the slope up to negative values. It is due to decline in release of Ca compared to Mg and subsequent decrease of Ca concentration, while Mg con- centration is continuously increasing. The incongruence is caused by precipita- tion of calcite, to which the solution is supersaturated, while still dissolving the Mg-bearing mineral. The composition of Mg-calcite determines shape of the re- action path. Model of concurrent dissolution of Mg-calcite and dolomite mix (representing limestone) shows very similar development with reaction path dependent on percent of Mg in calcite and dolomite/calcite ratio. Nevertheless, the effect of dolomite component is almost indistinguishable from pure Mg-cal- cite dissolution if the D/C ratio is less than 1 as the Mg-calcite kinetics are much faster. Therefore, in most limestones is the Mg-calcite composition a key factor determining Mg/Ca evolution. Compared to real dripwater data, most of the dripwaters follow the reac- tion paths defined by the dissolution model. The most probable cause of the data lying elsewhere is that the site conditions are different from conditions in Moravian Karst, upon which the model is based. For example, if the water was formed under higher PCO2, the resulting solution would show enhanced Ca con- centrations. Another defining parameter is temperature because of the rate constant dependence. The wide range of Ca concentrations and relatively narrow range of Mg concentrations indicates dissolution of very low-Mg calcites and only minor do- lomite component. This seems to be the case not only of Punkva Caves’ drip- waters but also of other cave systems, where the Ca concentrations in drip- waters show much wider range than Mg concentrations (Immenhauser et al. 2010; Riechelmann et al. 2011). In contrast, other sites evolve along the reac- tion paths for dissolution of limestones with a significant dolomite component or dolostone (Fairchild et al. 2000; Wong et al. 2011). These paths are gradually straighter with increasing dolomite component and the reaction path of dolo- mite is a straight line. Furthermore, dynamics of the modeled reaction paths are determined by the {L}/V ratio, the limestone-solution boundary. In addition, the ratio modifies
36 the reaction path shape in transition from congruent to incongruent dissolu- tion. Generally speaking, very high {L}/V ratio can be expected in epikarst, where the water is stored in pores and spaces between grains/clasts and frac- tures, whereas in fissures in deeper parts of vadose zone the {L}/V is much lower. Moreover, the perched aquifer forms an open system, where CO2 con- sumed by carbonates dissolution can be replenished from soil/sediments, whereas fissure flow is more of a closed system with very limited supply of free
CO2 in solution. These conditions seem to contribute to dissolution dominantly occurring in epikarst – as illustrated by the limestone corrosion diminishing downward the vadose zone. Therefore, although the value of {L}/V is generally unknown, it seems to be of high importance. Thus, the actual composition of dripwater is given not only by the reaction path shape, but also by its position along the reaction path, i.e., how far the dissolution proceeded (Appendix 3: Fig. 2 and 3). This underlines the im- portance of residence time in karst water formation. Under expected {L}/V ra- tio, the model of epikarst dissolution of Moravian Karst limestone indicates the residence times for regular dripwaters in the range from 100 to 150 days. These values are very well plausible compared to previously estimated residence times (Kamas et al. 2015; Faimon et al. 2016). Intriguingly, interpretation of the high Mg/Ca ratio in the anomalous drip based on the model would indicate extreme residence times in the range of hundreds of days. Nevertheless, the anomalous properties are most probably a result of PCP as previously ex- plained (Appendix 1). It shows that application of the model on real dripwater data could be limited because additional processes may participate on water hydrogeochemistry evolution in addition to dissolution. The arguably most studied process possibly influencing Mg/Ca ratio is the prior calcite precipita- tion (PCP) (Fairchild et al. 2000). Precipitating calcite consumes dissolved Ca and therefore increases the Mg/Ca ratio. It is believed that PCP is promoted by arid and warm climate, whereas the changes in Mg/Ca might be also caused by rapid rainfall events. During such events, infiltrating water might ‘flush out’ the water stored in semi-isolated reservoirs. Such water could show enhanced Mg/Ca ratio due to prolonged residence time. In addition, PCP leads to decrease
37 in Ca concentration, whereas Mg concentration remains the same, while the dissolution dynamics affect both Ca and Mg concentrations. This might be an important distinguishing factor in dripwater, although it is indistinguishable in precipitated speleothems. Other mechanisms participating on dripwater Mg/Ca ratios might be e.g. preferential Mg leaching from fresh surfaces (McGil- len and Fairchild 2005; Morse et al. 2007; Sinclair 2011), cation capture on mineral surfaces (e.g., on clay minerals present in epikarst) or dissolution of additional Mg-rich minerals (e.g., evaporites present in limestone). In paleoenvironmental studies evaluating the Mg/Ca ratio, the effect of PCP on dripwater formation (and therefore speleothem composition) is usually favored above any dissolution dynamic influence. This study identifies further significant factors besides incongruent dissolution due to different dissolution dynamics of carbonate minerals. These are (1) the limestone composition (Mg- calcite composition and dolomite content), (2) the residence time, (3) the ratio of limestone surface area to water volume (i.e., the karst system structure) and
(4) the epikarstic PCO2. With exception of PCO2, all these factors are intrinsically water-flow-path dependent. Therefore, they may change both spatially and temporarily with the natural evolution of the karst system, independently on climate changes. Finally, it is important to note, that very similar effects can be expected from mixing of waters formed at contact with different limestones.
3.4 ANTHROPOGENIC CO2 INFLUENCE ON DRIPWATERS AND SPELEOTHEM CORROSION Cave air ventilation was monitored in the Bear Chamber of the Výpustek Cave in Moravian Karst during events with large visitor numbers. The natural cave air CO2 levels in the Výpustek Cave (0.05–0.10 vol%, i.e., PCO2 = 10−3.32 to
10−2.98) are comparable with the Punkva Caves (0.06–0.32 vol%, i.e., PCO2 = 10−3.31 to 10−2.49), albeit being lower. The difference might be caused by season- ally more representative dataset from Punkva Caves. The study in the
Výpustek Cave showed that the anthropogenic influx of CO2 could lead to in- creased steady state concentrations up to PCO2(anthrop) = 10−2.22 (0.61 vol%), de- pending on the number of visitors and their duration of stay in the cave (for modeling details see Appendix 4).
38 Additionally, dripwaters were sampled for a study of hydrogeochemical properties. Due to limited occurrence of suitable permanent drips in the cave, only two, D1 and D2, were studied in detail. Both drips showed little seasonal variations. The drip D1 had higher specific conductivity (736–775 µS cm−1) com- pared to drip D2 (360–395 µS cm−1). Both dripwaters showed supersaturation with respect to calcite with higher values for D1 (SI = 1.03–1.19) than for D2
(SI = 0.58–0.64). Partial pressures of CO2 corresponding to aqueous carbonate species, PCO2(W), and hypothetical PCO2 under which the water was formed in soil/epikarst, PCO2(H), were calculated (for calculation background see chapter
3.2 and Appendix 2). Whereas the values for D1 (PCO2(H) = 10−(1.56–1.85)) were very similar to properties of other dripwaters in Moravian Karst (e.g., regular drip- waters in Punkva Caves show PCO2(H) = 10−(1.49–1.77)) the values for D2 (PCO2(H) = 10−(2.27–2.51)) were somewhat lower (Pracný et al. 2016). The ventilation model was compared to hydrogeochemical model and showed that the peak PCO2(anthrop) could principally exceed dripwater PCO2(H) turning the water into solution aggressive with respect to calcite. Further mod- eling and analysis was focused on drip D2 because it requires much lower CO2 concentrations than dripwater D1 to be converted into aggressive solution and is therefore of greater environmental concern. Depending on the ventilation mode the duration of visitors’ stay in cave required for PCO2(anthrop) to outreach
PCO2(H) is between 4.24–7.25 h for the ordinary visitor group (50 people) and 0.42–0.73 h for the enhanced attendance (500 people). For comparison, an in- dividual guided tour usually stays in the chamber for 0.25 h and this time would require thousands of visitors to balance the partial pressures. Such at- tendance is unreal given the cave capacity. Although an event with a large at- tendance and longer duration of stay (e.g., concert, performance or wedding) might cause conversion of some dripwater to a solution aggressive with respect to calcite, it seems that the cave dripwater cannot be converted by ordinary tours under current regime. Nevertheless, if condensation waters on cave walls are considered, the probability of corrosion increases. The condensing water dissolves gaseous CO2
39 and becomes aggressive to calcite. Therefore, any increase in cave air PCO2 re- sults in increasing potential of speleothem corrosion, although not via drip- water conversion.
40 4 CONCLUSIONS
The aim of the presented research of hydrogeochemical properties of drip- waters in Moravian Karst was to contribute to general understanding of condi- tions and mechanisms of dripwater formation. Insight into these processes pro- vides (1) consolidation of basic karst genesis concepts, (2) better paleoenviron- mental analyses and (3) knowledge applicable in more efficient karst/cave pro- tection. The long-term dripwater monitoring in Punkva Caves showed no imme- diate relation of meteorological conditions on karst surface (temperature, pre- cipitation) to dripwater properties. Dripwater discharge showed seasonal vari- ation with minima in summer/fall and maxima in winter/spring. The water amount seemed to be strongly affected by evapotranspiration and significantly recharged by snowmelt. Furthermore, not all drips show the same hydrological behavior, as some dripped only seasonally while other showed only very slight discharge variations. The hydrogeochemical properties of dripwaters in Moravian Karst were a foundation for the drips’ division into two groups: (a) the regular drips (high
SIcalcite and EC, low Mg/Ca and Sr/Ca ratios) and (b) an anomalous drip (low
SIcalcite and EC, high Mg/Ca and Sr/Ca ratios). The δ2H and δ18O isotopic com- positions of drips in both groups were the same, but the δ13C enrichment indi- cates anomalous dripwater degassing. In other properties (esp. discharge) the anomalous drip does not differ from regular drips. Prior calcite precipitation was identified as the most probable cause of the anomalous properties, followed by water mixing in the vadose zone. In addition, although the anomalous drip shows low saturation with respect to calcite, it seems incapable of becoming substantially undersaturated and corroding speleothems.
The study of CO2 in the air and water in Punkva Caves showed (I) season- ality of the cave air PCO2 with maxima in summer and minima in winter with the period slightly shorter than 1 year. The proposed cause is the cave ventila- tion driven by the difference of cave temperature and external temperature. Seasonality, albeit with different periods and less conclusive correlations, was
41 also showed by (II) the PCO2 corresponding to aqueous carbonates. The control- ling factor seems to be water degassing dynamics determined by the CO2 con- centration gradient between the water and cave air. A geochemical model was used to reconstruct (III) the hypothetical PCO2 in epikarst – the partial pressure of CO2 under which the water was formed. The results showed very stable val- ues without significant seasonality, indicating independence on surface condi- tions. Thus, a CO2 source is probably situated deeper in epikarst/vadose zone rather than in karst soils that are controlled seasonally. Possibility of speleothem corrosion via dripwater conversion under in- creased CO2 concentration due to anthropogenic influx was studied in the Výpustek Cave. Modeling showed that although it is possible to reach concen- trations causing undersaturation of dripwater with respect to calcite during events with enhanced attendance, ordinary guided tours seem to be incapable of substantially affecting dripwaters. Nevertheless, the effect on cave conden- sation waters might be of importance. The comparison of dripwater data with results of geochemical modeling of limestone dissolution under epikarstic conditions showed that Mg/Ca ratios are indeed controlled by incongruent dissolution of Mg-calcite, lesser of dolomite. Initially, the dissolution proceeds congruently following a straight path. When the solution reaches saturation with respect to calcite, the dissolution changes to incongruent, as part of material released from Mg-calcite or dolomite is pre- cipitated as calcite, while Mg accumulates and the Mg/Ca ratio slowly in- creases. The slope of reaction path nonlinearly increases and becomes negative, resulting in bent overall reaction path. Most of compared real dripwater data follow these reaction paths with a few exceptions caused probably by distinct conditions in individual karst systems (e.g., temperature and PCO2). The shape of reaction path (and resulting Mg/Ca ratio) primarily depends on dissolved rock composition. Besides, the Mg/Ca ratio is given by the dynamics of system evolution, i.e., by the distance along the reaction path that the system covers during dissolution. This distance is controlled by the extent of limestone/water boundary and by the overall water residence time.
42 It is apparent that the Mg/Ca ratios are dependent on additional factors beside the incongruent dissolution. These are (1) the limestone composition (ra- tio of calcite/dolomite component and Mg-calcite composition), (2) the time of water/rock interaction (water residence time), and (3) the ratio between water volume and rock surface during dissolution (effectively the karst structure). Independently on climatic conditions, all these factors might change with the variations in water flow paths. The results lead to the conclusion that numerous processes and factors usually considered negligible in paleoclimatic reconstructions might be of sub- stantial importance. These processes may disrupt desirable features (e.g., dis- solution dynamics affecting Mg/Ca ratio in the same fashion as PCP, or crystal growth and crystal habitus) and should be at least considered in paleorecon- structions. For example, both the model of Mg/Ca evolution and the difference in hydrogeochemical properties of the anomalous drip compared to the regular drips show much higher importance of water flow path in dripwater formation than anticipated. The flow path not only determines residence times, the lime- stone/water boundary, possibility of prior calcite precipitation or water mixing in vadose zone. More importantly, it can change both spatially and temporally independently on climatic variables as the karst system naturally evolves. A more complex analysis of various factors impacting dripwater forming speleo- them may result in more relevant interpretation of paleoenvironmental data.
43 5 REFERENCES
Aquilina L, Ladouche B, Dörfliger N (2006) Water storage and transfer in the epikarst of karstic systems during high flow periods. J Hydrol 327:472– 485. doi: 10.1016/j.jhydrol.2005.11.054 Arvidson RS, Ertan IE, Amonette JE, Luttge A (2003) Variation in calcite dis- solution rates: A fundamental problem? Geochim Cosmochim Acta 67:1623–1634. doi: 10.1016/S0016-7037(02)01177-8 Atkinson TC (1977) Carbon dioxide in the atmosphere of the unsaturated zone: An important control of groundwater hardness in limestones. J Hydrol 35:111–123. doi: 10.1016/0022-1694(77)90080-4 Auler AS, Smart PL (2004) Speleogenesis and Evolution of Karst Aquifers Rates of condensation corrosion in speleothems of semi-arid northeastern Brazil. Speleogenes Evol Karst Aquifers 2:1–2. Ayalon A, Bar-Matthews M, Kaufman A (1999) Petrography, strontium, bar- ium and uranium concentrations, and strontium and uranium isotope ra- tios in speleothems as palaeoclimatic proxies: Soreq Cave, Israel. The Hol- ocene 9:715–722. doi: 10.1191/095968399673664163 Baker A, Barnes W, Smart P (1997) Variations in the discharge and organic matter content of stalagmite drip waters in Lower Cave, Bristol. Hydrol Process 11:1541–1555. doi: 10.1002/(SICI)1099- 1085(199709)11:11<1541::AID-HYP484>3.0.CO;2-Z Baker A, Genty D (1998) Environmental pressures on conserving cave speleo- thems: Effects of changing surface land use and increased cave tourism. J Environ Manage 53:165–175. doi: 10.1006/jema.1998.0208 Baker A, Flemons I, Andersen MS, et al (2016) What determines the calcium concentration of speleothem-forming drip waters? Glob Planet Change 143:152–161. doi: 10.1016/j.gloplacha.2016.06.001 Baldini JUL, McDermott F, Fairchild IJ (2006) Spatial variability in cave drip water hydrochemistry: Implications for stalagmite paleoclimate records. Chem Geol 235:390–404. doi: 10.1016/j.chemgeo.2006.08.005 Baldini JUL, McDermott F, Baldini LM, et al (2012) Identifying short-term and seasonal trends in cave drip water trace element concentrations based on a daily-scale automatically collected drip water dataset. Chem Geol 330– 331:1–16. doi: 10.1016/j.chemgeo.2012.08.009 Baldini JUL, McDermott F, Hoffmann DL, et al (2008) Very high-frequency and seasonal cave atmosphere PCO2 variability: Implications for stalag- mite growth and oxygen isotope-based paleoclimate records. Earth Planet Sci Lett 272:118–129. doi: 10.1016/j.epsl.2008.04.031 Benavente J, Vadillo I, Carrasco F, et al (2010) Air Carbon Dioxide Contents in the Vadose Zone of a Mediterranean Karst. Vadose Zo J 9:126. doi: 10.2136/vzj2009.0027
44 Bernal JP, Cruz FW, Stríkis NM, et al (2016) High-resolution Holocene South American monsoon history recorded by a speleothem from Botuverá Cave, Brazil. Earth Planet Sci Lett 450:186–196. doi: 10.1016/j.epsl.2016.06.008 Blecha M, Faimon J (2014) Spatial and temporal variations in carbon dioxide (CO2) concentrations in selected soils of the Moravian Karst (Czech Re- public). Carbonates and Evaporites 29:395–408. doi: 10.1007/s13146-014- 0220-7 Borsato A, Johnston VE, Frisia S, et al (2016) Temperature and altitudinal influence on karst dripwater chemistry: Implications for regional-scale pal- aeoclimate reconstructions from speleothems. Geochim Cosmochim Acta 177:275–297. doi: 10.1016/j.gca.2015.11.043 Bottrell, SH. Atkinson, TC (1992) Tracer study of flow and storage in the un- saturated zone of a karstic limestone aquifer. Tracer Hydrology. Rotter- dam: Balkema, 207-211. Bourges F, Mangin A, d’Hulst D (2001) Carbon dioxide in karst cavity atmos- phere dynamics: the example of the Aven d’Orgnac (Ardeche). Comptes Rendus L Acad Des Sci Ser Ii Fasc a-Sciences La Terre Des Planetes 333:685–692. doi: 10.1016/s1251-8050(01)01682-2 Busenberg E, Plummer LN (1982) The kinetics of dissolution of dolomite in CO2-H2O systems at 1.5 to 65 degrees C and O to 1 atm PCO 2. Am. J. Sci. 282:45–78. Carroll D (1962) Rainwater as a Chemical Agent of Geologic Processes-A Re- view. US Geol Surv Water-Supply Pap 1535:1–18. Casteel RC, Banner JL (2015) Temperature-driven seasonal calcite growth and dripwater trace element variations in a well-ventilated Texas cave: Impli- cations for speleothem paleoclimate studies. Chem Geol 392:43–58. doi: 10.1016/j.chemgeo.2014.11.002 Cruz FW, Burns SJ, Jercinovic M, et al (2007) Evidence of rainfall variations in Southern Brazil from trace element ratios (Mg/Ca and Sr/Ca) in a Late Pleistocene stalagmite. Geochim Cosmochim Acta 71:2250–2263. doi: 10.1016/j.gca.2007.02.005 Cruz FW, Karmann I, Viana O, et al (2005) Stable isotope study of cave perco- lation waters in subtropical Brazil: Implications for paleoclimate infer- ences from speleothems. Chem Geol 220:245–262. doi: 10.1016/j.chemgeo.2005.04.001 Cubillas P, Köhler S, Prieto M, et al (2005) Experimental determination of the dissolution rates of calcite, aragonite, and bivalves. Chem Geol 216:59–77. doi: 10.1016/j.chemgeo.2004.11.009 D’Alessandro W, Katsanou K, Lambrakis N, et al (2013) Chemical and isotopic characterisation of bulk deposition in the Louros basin (Epirus, Greece). Atmos Res 132–133:399–410. doi: 10.1016/j.atmosres.2013.07.007 Dreybrodt W, Gabrovšek F, Perne M (2005) Condensation Corrosion: A Theo- retical Approach. Acta Carsologica 34:317–348. doi: 10.3986/ac.v34i2.262
45 Dreybrodt W, Scholz D (2011) Climatic dependence of stable carbon and oxygen isotope signals recorded in speleothems: From soil water to speleothem cal- cite. Geochim Cosmochim Acta 75:734–752. doi: 10.1016/j.gca.2010.11.002 Drysdale RN, Zanchetta G, Hellstrom JC, et al (2005) Stalagmite evidence for the onset of the Last Interglacial in southern Europe at 129 ± 1 ka. Ge- ophys Res Lett 32:1–4. doi: 10.1029/2005GL024658 Emblanch C, Zuppi GM, Mudry J, et al (2003) Carbon 13 of TDIC to quantify the role of the unsaturated zone: The example of the Vaucluse karst sys- tems (Southeastern France). J Hydrol 279:262–274. doi: 10.1016/S0022- 1694(03)00180-X Eriksson, E (1952) Composition of Atmospheric Precipitation. Tellus 4:280– 303. doi: 10.1111/j.2153-3490.1952.tb01014.x Eriksson E (1952) Composition of Atmospheric Precipitation. Tellus 4:215–232. Faimon J, Bodláková R, Pracný P, Hebelka J (2016) Transfer of climatic varia- bles by dripwater: a case study from Kateřinská Cave (Moravian Karst). Environ Earth Sci 75:1151. doi: 10.1007/s12665-016-5982-x Faimon J, Ličbinska M (2010) Carbon dioxide in the soils and adjacent caves of the Moravian Karst. Acta Carsologica 39:463–475. doi: 10.3986/ac.v39i3.76 Faimon J, Ličbinská M, Zajíček P, Sracek O (2012) Partial Pressures of CO2 in Epikarstic Zone Deduced from Hydrogeochemistry of Permanent Drips, the Moravian Karst, Czech Republic. Acta Carsologica 41:47–57. doi: 10.3986/ac.v41i1.47 Faimon J, Zimák J, Zajíček P, et al (2004) Současné destruktivní procesy v jeskyních Moravského krasu. Geographica 38:9–14. Fairchild IJ, Borsato A, Tooth AF, et al (2000) Controls on trace element (Sr- Mg) compositions of carbonate cave waters: Implications for speleothem climatic records. Chem Geol 166:255–269. doi: 10.1016/S0009- 2541(99)00216-8 Fairchild IJ, Smith CL, Baker A, et al (2006) Modification and preservation of environmental signals in speleothems. Earth-Science Rev 75:105–153. doi: 10.1016/j.earscirev.2005.08.003 Fairchild IJ, Treble PC (2009) Trace elements in speleothems as recorders of environmental change. Quat Sci Rev 28:449–468. doi: 10.1016/j.quasci- rev.2008.11.007 Fairchild IJ, Tuckwell GW, Baker A, Tooth AF (2006) Modelling of dripwater hydrology and hydrogeochemistry in a weakly karstified aquifer (Bath, UK): Implications for climate change studies. J Hydrol 321:213–231. doi: 10.1016/j.jhydrol.2005.08.002 Fairchild IJ, McMillan EA (2007) Speleothems as indicators of wet and dry pe- riods. Int J Speleol 36:69–74. doi: 10.5038/1827-806X.36.2.2 Ford, D, Williams, PD (2007) Karst hydrogeology and geomorphology. John Wiley & Sons.
46 Friederich H, Smart PL (1982) The classification of autogenic percolation wa- ters in karst aquifers : a study in G.B. Cave, Mendip Hills, England. Proc Univ Bristol Spelaeol Soc 16:143–159. Frisia S, Borsato A, Drysdale RN, et al (2012) A re-evaluation of the palaeocli- matic significance of phosphorus variability in speleothems revealed by high-resolution synchrotron micro XRF mapping. Clim Past 8:2039–2051. doi: 10.5194/cp-8-2039-2012 Frisia S, Borsato A, Fairchild IJ, McDermott F (2000) Calcite Fabrics, Growth Mechanisms, and Environments of Formation in Speleothems from the Italian Alps and Southwestern Ireland. J Sediment Res 70:1183–1196. doi: 10.1306/022900701183 Gabrovšek F, Dreybrodt W (2010) Karstification in unconfined limestone aqui- fers by mixing of phreatic water with surface water from a local input: A model. J Hydrol 386:130–141. doi: 10.1016/j.jhydrol.2010.03.015 Gatz DF (1991) Urban precipitation chemistry: A review and synthesis. Atmos Environ Part B Urban Atmos 25:1–15. doi: 10.1016/0957-1272(91)90036-E Genty D, Quinif Y (1996) Annually laminated sequences in the internal struc- ture of some Belgian stalagmites - Importance for paleoclimatology. J Sed- iment Res 66:275–288. Genty D, Deflandre G (1998) Drip flow variations under a stalactite of the Pere Noel cave (Belgium). Evidence of seasonal variations and air pressure con- straints. J Hydrol 211:208–232. doi: 10.1016/S0022-1694(98)00235-2 Glozar P, (1984) Průzkum komínových prostor v suché části Punkevních jeskyní (In Czech: Exploration of chimney spaces in dry parts of the Punkva Caves). Československý kras 35: 87-93. Goldscheider N, Meiman J, Pronk M, Smart C (2008) Tracer tests in karst hy- drogeology and speleology. Int J Speleol 37:27–40. doi: 10.5038/1827- 806X.37.1.3 Griffiths ML, Drysdale RN, Gagan MK, et al (2010) Evidence for Holocene changes in Australian-Indonesian monsoon rainfall from stalagmite trace element and stable isotope ratios. Earth Planet Sci Lett 292:27–38. doi: 10.1016/j.epsl.2010.01.002 Hartland A, Fairchild IJ, Lead JR, et al (2012) From soil to cave: Transport of trace metals by natural organic matter in karst dripwaters. Chem Geol 304–305:68–82. doi: 10.1016/j.chemgeo.2012.01.032 Huang Y, Fairchild IJ, Borsato A, et al (2001) Seasonal variations in Sr, Mg and P in modern speleothems (Grotta di Ernesto, Italy). Chem Geol 175:429–448. doi: 10.1016/S0009-2541(00)00337-5 Chou L, Garrels RM, Wollast R (1989) Comparative study of the kinetics and mechanisms of dissolution of carbonate minerals. Chem Geol 78:269–282. doi: 10.1016/0009-2541(89)90063-6 Immenhauser A, Buhl D, Richter D, et al (2010) Magnesium-isotope fractiona- tion during low-Mg calcite precipitation in a limestone cave – Field study
47 and experiments. Geochim Cosmochim Acta 74:4346–4364. doi: 10.1016/j.gca.2010.05.006 Jacobson RL, Langmuir D (1974) Dissociation constants of calcite and CaHC03+ from 0 to 50°C. Geochim Cosmochim Acta 38:301–318. doi: 10.1016/0016-7037(74)90112-4 Jochum KP, Scholz D, Stoll B, et al (2012) Accurate trace element analysis of speleothems and biogenic calcium carbonates by LA-ICP-MS. Chem Geol 318–319:31–44. doi: 10.1016/j.chemgeo.2012.05.009 Jones WK (2013) Physical structure of the epikarst. Acta Carsologica 42:311– 314. Kamas J, Bruthans J, Vysoka H, Kovařík M (2015) Range of horizontal transport and residence time of nitrate in a mature karst vadose zone. Int J Speleol 44:49–59. doi: 10.5038/1827-806X.44.1.5 Karmann I, Cruz FW, Viana O, Burns SJ (2007) Climate influence on geochem- istry parameters of waters from Santana-Pérolas cave system, Brazil. Chem Geol 244:232–247. doi: 10.1016/j.chemgeo.2007.06.029 Kaufmann G, Dreybrodt W (2007) Calcite dissolution kinetics in the system CaCO3-H2O-CO2 at high undersaturation. Geochim Cosmochim Acta 71:1398–1410. doi: 10.1016/j.gca.2006.10.024 Klimchouk AB (2004) Towards defining, delimiting and classifying epikarst: Its origin, processes and variants of geomorphic evolution. Proc Epikarst Conf Spec Publ 9 23–35. Kluge T, Riechelmann DFC, Wieser M, et al (2010) Dating cave drip water by tritium. J Hydrol 394:396–406. doi: 10.1016/j.jhydrol.2010.09.015 Kogovsek J, Petric M (2014) Solute transport processes in a karst vadose zone characterized by long-term tracer tests (the cave system of Postojnska Jama, Slovenia). J Hydrol 519:1205–1213. doi: 10.1016/j.jhy- drol.2014.08.047 Kogovšek J, Šebela S (2004) Water tracing through the vadose zone above Post- ojnska Jama, Slovenia. Environ Geol 45:992–1001. doi: 10.1007/s00254- 003-0958-z Kuzyakov Y (2006) Sources of CO2 efflux from soil and review of partitioning methods. Soil Biol Biochem 38:425–448. doi: 10.1016/j.soilbio.2005.08.020 Lachniet MS (2009) Climatic and environmental controls on speleothem oxy- gen-isotope values. Quat Sci Rev 28:412–432. doi: 10.1016/j.quasci- rev.2008.10.021 Lang M, Faimon J, Pracný P, Kejíková S (2017) A show cave management: Anthropogenic CO2 in atmosphere of Výpustek Cave (Moravian Karst, Czech Republic). J Nat Conserv 35:40–52. doi: 10.1016/j.jnc.2016.11.007 Li HC, Ku TL, You CF, et al (2005) 87Sr/86Sr and Sr/Ca in speleothems for paleoclimate reconstruction in Central China between 70 and 280 kyr ago. Geochim Cosmochim Acta 69:3933–3947. doi: 10.1016/j.gca.2005.01.009
48 Liss PS, Slater PG (1974) Flux of Gases across the Air-Sea Interface. Nature 247:181–184. doi: 10.1038/247181a0 Liu Z, Yuan D, Dreybrodt W (2005) Comparative study of dissolution rate-de- termining mechanisms of limestone and dolomite. Environ Geol 49:274– 279. doi: 10.1007/s00254-005-0086-z Mackenzie FT, Bischoff WD, Bishop FC, et al (1983) Magnesian calcites; low- temperature occurrence, solubility and solid-solution behavior. Rev Min- eral Geochemistry 11:97–144. Martín-Pérez A, Martín-García R, Alonso-Zarza A (2012) Diagenesis of a dra- pery speleothem from Castañar Cave: from dissolution to dolomitization. Int J Speleol 41:251–266. doi: 10.5038/1827-806X.41.2.11 McDermott F (2004) Palaeo-climate reconstruction from stable isotope varia- tions in speleothems: A review. Quat Sci Rev 23:901–918. doi: 10.1016/j.quascirev.2003.06.021 McDonald J, Drysdale R (2007) Hydrology of cave drip waters at varying bed- rock depths from a karst system in southeastern Australia. Hydrol Process 21:1737–1748. doi: 10.1002/hyp McDonald J, Drysdale R, Hill D, et al (2007) The hydrochemical response of cave drip waters to sub-annual and inter-annual climate variability, Wombeyan Caves, SE Australia. Chem Geol 244:605–623. doi: 10.1016/j.chemgeo.2007.07.007 McMillan EA, Fairchild IJ, Frisia S, et al (2005) Annual trace element cycles in calcite-aragonite speleothems: Evidence of drought in the western Med- iterranean 1200-1100 yr BP. J Quat Sci 20:423–433. doi: 10.1002/jqs.943 Meyer KW, Feng W, Breecker DO, et al (2014) Interpretation of speleothem calcite δ13C variations: Evidence from monitoring soil CO2, drip water, and modern speleothem calcite in central Texas. Geochim Cosmochim Acta 142:281–298. doi: 10.1016/j.gca.2014.07.027 Milanolo S, Gabrovšek F (2015) Estimation of carbon dioxide flux degassing from percolating waters in a karst cave: Case study from Bijambare cave, Bosnia and Herzegovina. Chemie der Erde - Geochemistry 75:465–474. doi: 10.1016/j.chemer.2015.10.004 Moore PJ, Martin JB, Screaton EJ (2009) Geochemical and statistical evidence of recharge, mixing, and controls on spring discharge in an eogenetic karst aquifer. J Hydrol 376:443–455. doi: 10.1016/j.jhydrol.2009.07.052 Morse JW, Arvidson RS (2002) The dissolution kinetics of major sedimentary carbonate minerals. Earth-Science Rev 58:51–84. doi: 10.1016/S0012- 8252(01)00083-6 Morse JW, Arvidson RS, Lüttge A (2007) Calcium carbonate formation and dis- solution. Chem Rev 107:342–381. doi: 10.1021/cr050358j Morse JW, Mackenzie FT (1990) Geochemistry of sedimentary carbonates. Elsevier
49 Musgrove M, Banner JL (2004) Controls on the spatial and temporal variability of vadose dripwater geochemistry: Edwards aquifer, central Texas. Geo- chim Cosmochim Acta 68:1007–1020. doi: 10.1016/j.gca.2003.08.014 Niggemann S, Mangini A, Richter DK, Wurth G (2003) A paleoclimate record of the last 17,600 years in stalagmites from the B7 cave, Sauerland, Ger- many. Quat Sci Rev 22:555–567. doi: 10.1016/S0277-3791(02)00143-9 Orland IJ, Burstyn Y, Bar-Matthews M, et al (2014) Seasonal climate signals (1990-2008) in a modern Soreq Cave stalagmite as revealed by high-reso- lution geochemical analysis. Chem Geol 363:322–333. doi: 10.1016/j.chemgeo.2013.11.011 Otava J (1995) Isotopic analyses and origin of CO2 in some moravian caves. Acta Carsologica XXIV:339–446. Palmer, AN (2006) Digital modeling of karst aquifers—successes, failures, and promises. Geological Society of America Special Papers, 404: 243-250 Paar D, Frančišković-Bilinski S, Buzjak N, et al (2016) Geochemical and min- eralogical characterization of speleothems from the karst of Croatia as po- tential sources of data for environmental researches. J Geochemical Ex- plor. doi: 10.1016/j.gexplo.2016.05.004 Paternoster M, Sinisi R, Mancusi C, et al (2014) Natural versus anthropogenic influences on the chemical composition of bulk precipitation in the south- ern Apennines, Italy: A case study of the town of Potenza. J Geochemical Explor 145:242–249. doi: 10.1016/j.gexplo.2014.07.003 Perrin J, Jeannin PY, Zwahlen F (2003) Implications of the spatial variability of infiltration-water chemistry for the investigation of a karst aquifer: A field study at Milandre test site, Swiss Jura. Hydrogeol J 11:673–686. doi: 10.1007/s10040-003-0281-5 Perrin J, Jeannin P-Y, Cornaton F (2007) The role of tributary mixing in chem- ical variations at a karst spring, Milandre, Switzerland. J Hydrol 332:158– 173. doi: 10.1016/j.jhydrol.2006.06.027 Perrin J, Jeannin PY, Zwahlen F (2003) Epikarst storage in a karst aquifer: A conceptual model based on isotopic data, Milandre test site, Switzerland. J Hydrol 279:106–124. doi: 10.1016/S0022-1694(03)00171-9 Peyraube N, Lastennet R, Denis a. (2012) Geochemical evolution of groundwa- ter in the unsaturated zone of a karstic massif, using the P CO2-SIc rela- tionship. J Hydrol 430–431:13–24. doi: 10.1016/j.jhydrol.2012.01.033 Peyraube N, Lastennet R, Denis a., Malaurent P (2013) Estimation of epikarst air PCO2 using measurements of water δ13CTDIC, cave air PCO2 and δ13CCO2. Geochim Cosmochim Acta 118:1–17. doi: 10.1016/j.gca.2013.03.046 Plestenjak G, Eler K, Vodnik D, et al (2012) Sources of soil CO2 in calcareous grassland with woody plant encroachment. J Soils Sediments 12:1327– 1338. doi: 10.1007/s11368-012-0564-3
50 Plummer LN, Busenberg E (1982) The solubilities of calcite, aragonite and va- terite in CO2-H2O solutions between 0 and 90°C, and an evaluation of the aqueous model for the system CaCO3-CO2-H2O. Geochim Cosmochim Acta 46:1011–1040. doi: 10.1016/0016-7037(82)90056-4 Plummer LN, Mackenzie FT (1974) Predicting mineral solubility from rate data; application to the dissolution of magnesian calcites. Am. J. Sci. 274:61–83. Plummer LN, Wigley TML, Parkhurst DL (1978) The kinetics of calcite disso- lution in CO 2 -water systems at 5 degrees to 60 degrees C and 0.0 to 1.0 atm CO 2. Am J Sci 278:179–216. doi: 10.2475/ajs.278.2.179 Pokrovsky OS, Golubev S V., Schott J (2005) Dissolution kinetics of calcite, dolomite and magnesite at 25 °C and 0 to 50 atm pCO2. Chem Geol 217:239–255. doi: 10.1016/j.chemgeo.2004.12.012 Pokrovsky OS, Golubev S V., Schott J, Castillo A (2009) Calcite, dolomite and magnesite dissolution kinetics in aqueous solutions at acid to circumneu- tral pH, 25 to 150 °C and 1 to 55 atm pCO2: New constraints on CO2 se- questration in sedimentary basins. Chem Geol 260:317–329. doi: 10.1016/j.chemgeo.2009.01.013 Pracný P, Faimon J, Sracek O, et al (2016) Anomalous drip in the Punkva caves (Moravian Karst): relevant implications for paleoclimatic proxies. Hydrol Process 30:1506–1520. doi: 10.1002/hyp.10731 Pronk M, Goldscheider N, Zopfi J, Zwahlen F (2008) Percolation and particle transport in the unsaturated zone of a karst aquifer. Ground Water 47:361–9. Riechelmann DFC, Schröder-Ritzrau A, Scholz D, et al (2011) Monitoring Bun- ker Cave (NW Germany): A prerequisite to interpret geochemical proxy data of speleothems from this site. J Hydrol 409:682–695. doi: 10.1016/j.jhydrol.2011.08.068 Riechelmann S, Schröder-Ritzrau A, Wassenburg JA, et al (2014) Physicochem- ical characteristics of drip waters: Influence on mineralogy and crystal morphology of recent cave carbonate precipitates. Geochim Cosmochim Acta 145:13–29. doi: 10.1016/j.gca.2014.09.019 Roberts MS, Hawkesworth CJ, Smart PL, et al (1999) Trace element variations in coeval Holocene speleothems from GB Cave, southwest England. The Holocene 9:707–713. doi: 10.1191/095968399672615014 Sander R (2015) Compilation of Henry’s law constants (version 4.0) for water as solvent. Atmos Chem Phys 15:4399–4981. doi: 10.5194/acp-15-4399- 2015 Sanchez-Cañete EP, Serrano-Ortiz P, Kowalski a. S, et al (2011) Subterranean CO2 ventilation and its role in the net ecosystem carbon balance of a karstic shrubland. Geophys Res Lett 38:8–11. doi: 10.1029/2011GL047077
51 Schwarz K, Barth JAC, Postigo-Rebollo C, Grathwohl P (2009) Mixing and transport of water in a karst catchment: a case study from precipitation via seepage to the spring. Hydrol Earth Syst Sci 13(3):285–292. doi: PNR61 Sherman LA, Barak P (2000) Solubility and Dissolution Kinetics of Dolomite in Ca–Mg–HCO/CO Solutions at 25°C and 0.1 MPa Carbon Dioxide. Soil Sci Soc Am J 64:1959. doi: 10.2136/sssaj2000.6461959x Sinclair DJ (2011) Two mathematical models of Mg and Sr partitioning into solution during incongruent calcite dissolution. Implications for dripwater and speleothem studies. Chem Geol 283:119–133. doi: 10.1016/j.chemgeo.2010.05.022 Sinclair DJ, Banner JL, Taylor FW, et al (2012) Magnesium and strontium sys- tematics in tropical speleothems from the Western Pacific. Chem Geol 294– 295:1–17. doi: 10.1016/j.chemgeo.2011.10.008 Smart PL, Friederich H (1987) Water movement and storage in the unsatu- rated zone of a maturely karstified carbonate aquifer, Mendip Hills, Eng- land. In: Proceedings of the Environmental Problems in Karst Terranes and their Solutions Conference. p pp.57-87. Smart, C, Worthington, SRH (2004) Groundwater in karst. Encyclopedia of cave and karst science. Fitzroy Dearborn, New York, 394-397. Spötl C, Fairchild IJ, Tooth AF (2005) Cave air control on dripwater geochem- istry, Obir Caves (Austria): Implications for speleothem deposition in dy- namically ventilated caves. Geochim Cosmochim Acta 69:2451–2468. doi: 10.1016/j.gca.2004.12.009 Stumm W, Morgan JJ (1996) Aquatic Chemistry: Chemical Equilibria and Rates in Natural Waters. John Wiley & Sons, New York Tan L, Shen CC, Cai Y, et al (2014) Trace-element variations in an annually layered stalagmite as recorders of climatic changes and anthropogenic pol- lution in central China. Quat Res (United States) 81:181–188. doi: 10.1016/j.yqres.2013.12.001 Ternan JL (1972) Comments on the use of a calcium hardness variability index in the study of carbonate aquifers: With reference to the central pennines, England. J Hydrol 16:317–321. doi: 10.1016/0022-1694(72)90136-9 Tooth AF, Fairchild IJ (2003) Soil and karst aquifer hydrological controls on the geochemical evolution of speleothem-forming drip waters, Crag Cave, southwest Ireland. J Hydrol 273:51–68. doi: 10.1016/S0022- 1694(02)00349-9 Trček, B (2003). Epikarst zone and the karst aquifer behaviour: a case study of the Hubelj catchment, Slovenia (p. 100). Geološki zavod Slovenije. Tremaine DM, Froelich PN (2013) Speleothem trace element signatures: A hy- drologic geochemical study of modern cave dripwaters and farmed calcite. Geochim Cosmochim Acta 121:522–545. doi: 10.1016/j.gca.2013.07.026 Verheyden S, Keppens E, Fairchild IJ, et al (2000) Mg, Sr and Sr isotope geo- chemistry of a Belgian Holocene speleothem: Implications for paleoclimate
52 reconstructions. Chem Geol 169:131–144. doi: 10.1016/S0009- 2541(00)00299-0 Verheyden S, Nader FH, Cheng HJ, et al (2008) Paleoclimate reconstruction in the Levant region from the geochemistry of a Holocene stalagmite from the Jeita cave, Lebanon. Quat Res 70:368–381. doi: 10.1016/j.yqres.2008.05.004 Vet R, Artz RS, Carou S, et al (2014) A global assessment of precipitation chem- istry and deposition of sulfur, nitrogen, sea salt, base cations, organic ac- ids, acidity and pH, and phosphorus. Atmos Environ 93:3–100. doi: 10.1016/j.atmosenv.2013.10.060 Williams PW (2008) The role of the epikarst in karst and cave hydrogeology: A review. Int. J. Speleol. 37:1–10. Wong CI, Banner JL, Musgrove M (2011) Seasonal dripwater Mg/Ca and Sr/Ca variations driven by cave ventilation: Implications for and modeling of spe- leothem paleoclimate records. Geochim Cosmochim Acta 75:3514–3529. doi: 10.1016/j.gca.2011.03.025
53 APPENDIX 1
This appendix presents following research paper:
Pracný, P., Faimon, J., Sracek, O., Kabelka, L., & Hebelka, J. (2016). Anomalous drip in the Punkva caves (Moravian Karst): relevant implications for paleoclimatic proxies. Hydrological Processes, 30(10), 1506–1520. http://doi.org/10.1002/hyp.10731
© 2015 John Wiley & Sons, Ltd. The original publication is available at Wiley via http://doi.wiley.com/10.1002/hyp.10731
54 HYDROLOGICAL PROCESSES Hydrol. Process. 30, 1506–1520 (2016) Published online 30 November 2015 in Wiley Online Library (wileyonlinelibrary.com) DOI: 10.1002/hyp.10731
Anomalous drip in the Punkva caves (Moravian Karst): relevant implications for paleoclimatic proxies
Pavel Pracný,1* Jiří Faimon,1,2 Ondra Sracek,2 Ludvík Kabelka3 and Jiří Hebelka4 1 Department of Geological Sciences, Faculty of Science, Masaryk University, Kotlářská 267/2, 611 37, Brno, Czech Republic 2 Department of Geology, Faculty of Science, Palacky University, 17. listopadu 12, 771 46, Olomouc, Czech Republic 3 Hydrochemical Laboratories, GEOtest, a.s., Šmahova 1244/112, 627 00, Brno, Czech Republic 4 Cave Administration of the Czech Republic, Svitavská 11, 678 01, Blansko, Czech Republic
Abstract: The anomalous drip in the Punkva caves (Moravian Karst) shows specific hydrogeochemical properties such as low 3 1 13 SIcalcite ~ 0.14 ± 0.11 (standard deviation), low mineralization (4.53 ± 0.42) × 10 mol l , and enhanced values of δ C( 7.85 to 8.35‰ VPDB), Mg/Ca × 1000 ratio (45.7 ± 3.3), and Sr/Ca × 1000 ratio (0.65 ± 0.06). By these properties, the anomalous drip significantly differs from other regular drips in the same cave and other caves in the region. The study suggests that the anomalous drip properties are a consequence of prior calcite precipitation or/and water mixing along the water flow path. As the former processes are spatially controlled, the knowledge of dripwater flow path seems to be necessary for correct paleoclimatic/paleoenvironmental reconstructions. Copyright © 2015 John Wiley & Sons, Ltd.
13 18 KEY WORDS cave dripwater; stable isotope C and O; Mg/Ca and Sr/Ca ratios; prior calcite precipitation; karst water mixing Received 24 June 2015; Accepted 21 October 2015
INTRODUCTION comprehensive characterization of several drips from the same cave/region is a critical step in defining possible Autochthonous cave deposits (speleothems) are frequent- differences in stalagmite records. ly used as terrestrial archives of paleoclimatic data (e.g. The goal of the study is to introduce an anomalous Fairchild et al., 2006 or Fairchild and Treble, 2009, for a dripwater from the Punkva caves (Moravian Karst), the review). This utilization is based on an assumption that composition of which is strongly influenced by processes external surface conditions are projected/encoded into in the epikarst/vadose zone. Its hydrology and chemistry speleothems (Verheyden et al., 2003; Li et al., 2005; Cruz is compared with the ‘regular dripwaters’ from different et al., 2007; Verheyden et al., 2000; Griffiths et al., sites in the same cave and in other caves in the same 2010). Unfortunately, the link between the surface region. The anomalous properties of water are expected to conditions and speleothems is influenced by many factors change dramatically morphological/geochemical proper- on the water reaction–transport path that may disturb the ties of a potentially formed speleothem: growth fabrics climatic signal (Baldini et al., 2006; Karmann et al., (lamina thickness), trace elements contents (Mg, Sr), and 2007; McDonald et al., 2007; Miorandi et al., 2010; stable isotopes values (δ13C). A prospective Sherwin and Baldini, 2011). Therefore, it is extremely paleoenvironmental study based on such properties might important to study recent processes in detail to see how lead to incorrect interpretations. they influence the climatic proxies. Some studies have shown that coeval stalagmites from the same cave may exhibit different trace element patterns (Roberts et al., SITE OF STUDY 1999; Finch et al., 2003). This either suggests that some stalagmite properties are not representative of climatic The Punkva caves are the best known and most popular conditions or that different drip sites preserve distinct show caves in the Czech Republic. They are situated components of the climate signal. Therefore, a more approximately 20 km northwards from Brno in the Moravian Karst, the largest karst region in the country. The caves have developed in very pure Devonian *Correspondence to: Pavel Pracný, Department of Geological Sciences, limestone of the Macocha formation, specifically in the Faculty of Science, Masaryk University, Kotlářská 267/2, 611 37 Brno, ž Czech Republic. La ánky and Vilémovice Limestones (Faimon et al., E-mail: [email protected] 2012; Blecha and Faimon, 2014; Pracný et al., 2015).
Copyright © 2015 John Wiley & Sons, Ltd. 55 ANOMALOUS DRIP: IMPLICATIONS FOR PALEOCLIMATIC PROXIES 1507
The study was conducted during the period from Meteorological data are from the meteorological station February 2012 to March 2013 in the Punkva caves situated close to the Macocha Abyss Upper Bridge and (Moravian Karst). Dripwaters were sampled in a corridor run by the Czech Meteorological Institute and the Cave behind Přední Chamber (drips CP1, CP2, and CP3), Administration of the Czech Republic (Figure 1). Tunnel Corridor (drips TC1/the anomalous drip, TC2), A set of archive data was provided for a comparison. and Zadní Chamber (drip ZD) in the so-called dry part of The data were collected from June 2003 to May 2004. the caves (Figure 1). Waters in the CP1, CP2, CP3, and The data come from a remote part of the Punkva caves TC2 sites have dripped from small straw stalactites on the and three different caves in the Moravian Karst (Figure 1). corridor ceiling. The water of TC1 has dripped from a Two drips, one from a stalactite, MD1, and one from a curtain about 30 cm wide developed on the edge of a straw stalactite, MD2, are situated in the Masaryk Dome crevasse. The water of ZD dripped on the top of Chamber, Punkva caves. Another drip comes from a stalagmite ‘Vase’ from the height of about 20 m. The straw stalactite in Big Foch’s Dome Chamber in the sampling was conducted twice per month: 126 samples in Balcarka cave (BC). Two drips are in the Amatérská cave; total were collected during 26 sampling courses. Because the first one drips from a straw stalactite in the crossroad of their extremely low drip rate during part of the year, the of the adit and Javor Corridor (AC1), and the second one drips TC2 and ZD were not sampled in all cases. comes from a curtain in the Rozlehlá Corridor (AC2).
Figure 1. Sketch map of the Punkva caves system and its localization within the Moravian Karst and the Czech Republic. For explanation of the drip site acronyms, see the text
Copyright © 2015 John Wiley & Sons, Ltd. Hydrol. Process. 30, 1506–1520 (2016) 56 1508 P. PRACNÝ ET AL.
There were nine sampling courses with chemical analysis regimes: the drip CP2 was very stable (variation for all the archive drips and six courses with just drip coefficient 17.6%; Figure 2e), whereas drips CP3 and discharge measurement. TC1 were more variable (variation coefficients 59% and 43.2%, respectively, Figures 2f and g). More variable rates were found for the drip CP1 (from 6 to MATERIALS AND METHODS 174 drops min 1, variation coefficient 87.1%; Figure 2d) 1 Drip discharge was measured by counting drips during a and ZD (rate from 3.5 to 167 drops min , variation fi given time period. The drip size was found by measuring coef cient 149.8%; Figure 2c). The drip TC2 was fi the weight of a water drop caught into a plastic container extremely variable (variation coef cient 230%). Its rate < (<0.9 g) using digital scales. An average value based on was extremely slow ( 1 drop per 5 m) or zero during a three measurements from individual drip sites was used substantial part of the monitoring period. for drip rate calculations. The archive data (not presented in Figure 2) show the – 1 Immediately in the cave, the basic hydrogeochemical most variable drip AC2 (the rate 50 400 drops min , fi parameters were determined: pH, specific electrical con- variation coef cient 77.6%). Other drips are much slower – 1 ductivity (EC), alkalinity (by acidimetric titration with (the drip BC with rate of 0.5 4 drops min , variation fi evaluation of titration curve via the Gran’s function, coef cient 60.9%; the drips MD1 and AC1 with the rates – 1 fi Stumm and Morgan, 1996), and Ca concentration of 8 18 drops min , variation coef cient below 25%; the – 1 (complexometric microtitration using calcein as indica- drip MD2 with the rate of 3 5 drops min ; variation fi tor). Other analyses (Mg, Sr) were conducted in a coef cient below 25%). laboratory using ICP (iCAP 6000 by Thermo Scientific). Saturation indices were calculated using the PHREEQC Hydrogeochemistry code (Parkhurst and Appelo, 2013). Totally, 126 new Specific ECs of the common drips TC2, CP1, CP2, CP3, water samples were analysed and compared with 45 and ZD (representing water mineralization) varied in the analyses from the archive dataset. range 604 ± 32 μScm 1 (standard deviation) (Figure 3a). Dripwaters for analyses of stable isotopes were Mean values of other parameters were as follows: sampled twice, in spring (April) and fall (November) of pH ~ 8.06 ± 0.13, alkalinity ~ (5.68 ± 1.34) × 10 3 eq l 1, 2014, and stored in airtight sealed plastic bottles. Isotope Ca concentrations ~ (3.38 ± 0.19) × 10 3 mol l 1, Mg con- δ18 δ2 values O and H in water were determined at the centrations ~ (5.71 ± 0.31) × 10 5 mol l 1, and Sr concen- Czech Geological Survey in Prague, Czech Republic, trations ~ (1.03 ± 0.08) × 10 6 mol l 1. The Mg/Ca × 1000 using a Los Gatos Research laser absorption spectrom- and Sr/Ca × 1000 ratios for the drips were 17.0 ± 1.4 eter. The results were normalized to the standard Vienna (Figure 3c) and 0.31 ± 0.02 (Figure 3e), respectively. Standard Mean Ocean Water (VSMOW) and reported in In contrast, the drip TC1 showed lower EC ~ 297 δ the -notation. The reproducibility of measurements was ± 22.2 μScm 1 (Figure 3a), alkalinity ~ (2.14 ‰ δ2 ‰ δ18 δ13 0.5 for H and 0.08 for O. Precipitates for C ±0.28)×10 3 eq l 1, and Ca concentrations ~ (1.47 ± 0.13) analyses were prepared by adding NaOH and BaCl2 and mol l 1. The ratios Mg/Ca × 1000 ~ 45.7 ± 3.3 and then the precipitate was filtered out. In the next step, the 13 Sr/Ca × 1000 ~ 0.62 ± 0.06 were enhanced. The mean Mg BaCO3 precipitate was dissolved by H3PO4, and δ C 5 ‰ and Sr concentrations, (6.71 ± 0.37) × 10 and (9.59 ± was measured with a precision better than 0.05 . Results 0.74) × 10 7 mol l 1, respectively, were roughly compara- were expressed with respect to the Vienna Pee Dee ble with those in other drips. Belemnite (VPDB) standard. Archive data showed the parameters similar to the common drips. The EC values were slightly reduced: the RESULTS dripwaters of the Masaryk Dome Chamber from the Punkva Caves showed EC 437 ± 30 μScm 1. The drips in Hydrology the Amatérská cave showed EC 495 ± 35 μScm 1.Asan During the monitoring period, daily precipitation exception, the occasional drip in the Balcarka cave reached up to 30 mm. Overall precipitation for the entire showed very low EC 311 ± 14 μScm 1, which was near monitoring period was 641.7 mm (Figure 2a). External the values of the drip TC1 (Figure 3b). The Mg/Ca × 1000 daily temperatures ranged from 9.3 °C to 25.9 °C. The ratios varied between 10 and 25 with the only exception mean temperature for the entire period was 7.89 °C for the drip AC1 showing much higher ratio (63 ± 3.9) (Figure 2b). Drip discharges in the cave reached up to (Figure 3d). In this case, however, the drip AC1 showed 180 drops min 1 with the local minima in winter enhanced Mg concentrations, (2.53 ± 0.50) × 10 4 mol l 1. (November to February) and the local maxima in spring For the archive data, the Sr/Ca × 1000 ratios were not (March to May). Studied drips showed different flow available.
Copyright © 2015 John Wiley & Sons, Ltd. Hydrol. Process. 30, 1506–1520 (2016) 57 ANOMALOUS DRIP: IMPLICATIONS FOR PALEOCLIMATIC PROXIES 1509
Figure 2. Meteorological and hydrological data. Daily rainfall (a), external maximum and minimum temperatures (b), drip discharges of ZD (c), CP1, 2, 3 (d, e, f), anomalous TC1 (g), and TC2 (h). The orange columns and the light blue columns indicate the maximum daily temperature over the freezing point and the minimum temperature below the freezing point, respectively, during winter/spring period
Copyright © 2015 John Wiley & Sons, Ltd. Hydrol. Process. 30, 1506–1520 (2016) 58 1510 P. PRACNÝ ET AL.
Figure 3. Selected hydrogeochemical data of drip waters: electrical conductivity, main dataset (a) and archive data (b); (Mg/Ca) × 1000 ratios, main dataset (c) and archive data (d); (Sr/Ca) × 1000 ratios, main dataset (e)
The results for all dripwater samples on δ2H and δ18O from 71.16 to 73.35‰ VSMOW for δ2H and from isotopes are in Figure 4. The spring values ranged from 10.17 to 10.41‰ VSMOW for δ18O. The results for 70.70 to 74.20‰ VSMOW for δ2H and from 10.30 δ13C are in Figure 5. Majority of dripwater samples to 10.60‰ VSMOW for δ18O. The fall values ranged showed the δ13C values in the range from 10.34 to
Copyright © 2015 John Wiley & Sons, Ltd. Hydrol. Process. 30, 1506–1520 (2016) 59 ANOMALOUS DRIP: IMPLICATIONS FOR PALEOCLIMATIC PROXIES 1511
Figure 4. Dripwater isotopes δ2D and δ18O
Figure 5. Dripwater alkalinity versus δ13C