Palmer.indd 143 �AB STRA CT Determining the nature of water fl ow and contaminant dispersion in karst requires far more information than can be provided by simple dye traces. Tracing can delineate drainage divides, fl ow directions, and fl ow velocities at various stages, but from water management purposes it is also important to determine such variables as groundwater storage, retention times, patterns of convergence and divergence, and response to wet-dry cycles in the soil. These are most signifi cant in the non-conduit portions of the karst aquifer, which supply most wells. Dye tracing can be augmented by hydrograph analysis at various stages, tracing with tagged solid particles or microbes, evaluation of dissolved sol- ids and chemical equilibria, and isotopic analysis. This paper concentrates on some of the uses of chemical equilibria and isotopes. Stable isotopes (e.g. 18O and deuterium) and the various radium isotopes are among the most useful. Ratios among the four radium isotopes (228Ra and 224Ra, with half-lives in years; and 223Ra and 226Ra with half-lives in days) are well suited to karst studies. These techniques are time-consuming and costly, so a full analysis of a karst aquifer is rarely feasible. Instead, it is recommended that selective analyses be made of representative parts of the aquifer, and that they be applied as follows: (1) Develop conceptual models based on fi eld observation, which allow one to anticipate a range of probable scenarios of contaminant transport and remediation. (2) If digital models are used, it is most effective to design simple generalized models in which the boundary conditions are clearly defi ned, and then to gain insight into real aquifers by noting the differences between the model and fi eld observations. (3) Use fi eld techniques to become familiar with the local hydrology and then apply hydraulic and chemical principles to an- ticipating contaminant behaviour, rather than reacting only to emergencies. These approaches encourage the growth of interpretive skills based on the same scientifi c principles that govern the origin of caves and karst. Keywords: hydrologic models, geochemistry, isotopes, contaminant tracking Understanding the hydrology of karst � Arthur N. Palmer Department of Earth Sciences, State University of New York, Oneonta, NY 13820, USA; (palmeran@oneonta.edu) doi: 104154/gc.2010.11 Geologia Croatica 63/2 143–148 2 Figs. Zagreb 2010 Geologia CroaticaGeologia Croatica 1. INTRODUCTION Maintaining a safe water supply is diffi cult in karst, where the physical setting and modes of water movement are complex. This paper addresses two issues: fi rst, to introduce the great variety of interpretive fi eld techniques that are available to the karst hydrologist; and second, to suggest a realistic strategy for applying the resulting information to karst. Effective water management in karst requires several types of fi eld data: (1) Delineation of drainage basins and drainage divides. This task is complicated by the overlap of catchment areas where perched vadose fl ow crosses phreatic divides, and by the shifting of divides with time as fl ow stage varies. (2) Flow directions, velocities, and discharge in solu- tion conduits and in the dispersed fl ow outside the conduits. (3) Patterns of convergence and divergence of groundwater, and their role in chemical dispersion. (4) Effective pore and fi ssure sizes and their interconnectivity. This informati on helps the anticipation of likely paths of pathogen migration. (5) Volume and duration of water retained in storage betwe en infiltration events (e.g., capillary water, perched pools, etc.). Geologia Croatica 63/2Geologia Croatica 144 ence of dispersion of water into the aquifer. During a fl ood pulse, the infl ow at fi rst typically exceeds the outfl ow, indi- cating accumulation of storage in conduits, abandoned up- per-level cave passages, and neighbouring fi ssures and pores. During the waning phase of a fl ood pulse, the outfl ow ex- ceeds the infl ow as storage is released. A long tail with a roughly logarithmic rate of decrease suggests the release of water as laminar fl ow, after the main turbulent-fl ow pulse has passed (ATKINSON, 1977). The volume released is found by integrating the area under the hydrograph (volume = ∑QΔt, where Q = discharge and t = time. Response to storm events will vary with antecedent pre- cipitation – for example, after a prolonged dry period when the infi ltration capacity of soil and fi ll material in epikarst fi ssures is dominated by desiccation fi ssures, in comparison with a prolonged wet period when soils are moist and have expanded to fi ll available bedrock fi ssures. Much of the wa- ter that reaches the phreatic zone during an infi ltration event is composed of capillary water that has been held in storage from prior events and displaced by the new water pulse (shown by PITTY (1966)). 2.3. Measuring travel times: quantitative tracer tests Estimating the position of drainage divides with dye traces is a long-established technique. Detection of fl ow diver- gence, crossovers, and shifting of divides with fl ow stage are also important outcomes. These complexities increase with structural deformation. Dominantly convergent and well-de- fi ned drainage patterns are typical of undisturbed strata, as in the central U.S.A. (e.g. QUINLAN & RAY, 1981). In areas of greater structural complexity, such as the Dinaric chain (e.g., BAUČIĆ, 1968), features including divergent fl ow paths, crossovers, and overlapping groundwater divides are far more likely. Dye tracing is discussed in many other pub- lications (e.g., KÄSS, 1998). Much information can be gained by repeating quantita- tive dye traces at a wide variety of discharges. Dye concen- trations at detection points are measured at short time inter- vals so the exact arrival times, peak dye concentrations, and mean travel times (to centres of mass) can be detected. Mean breakthrough times to each monitoring point are plotted as a function of discharge. A logarithmic scale allows coverage of a broad range of values and emphasizes differences in the shape of the plots (PALMER, 2007). Conduits that are mostly water-fi lled will show linear plots with slopes of –1 (Fig. 1). Paths composed mainly of open channels, such as canyons, produce concave-upward curves: as discharge rises, water depth also increases, and the increase in wetted perimeter, (p, length of contact of water with solid surfaces in a cross section), can increase almost as fast as the cross-sectional area (A). Flow velocity is a function of the ratio A/p, so ve- locity in a fl ooding vadose channel tends to increase more slowly than in a water-fi lled conduit. A combination of open channels and closed conduits will give linear plots at low discharge that change to concave-upward curves at high dis- charge. Complex patterns, such as overfl ow and divergence, 2. ANALYTICAL TECHNIQUES The variables listed above can be investigated in many dif- ferent ways. Major categories of hydrologic fi eld data are outlined here, mainly to demonstrate their utility and how they relate to each other. Most are time-consuming, and some are costly as well, so it is rarely feasible to apply them in suffi cient detail to truly understand the behaviour of a karst aquifer. However, by considering them, we are forced to think about how karst aquifers behave. For further informa- tion see also MILANOVIĆ, 1981; BONACCI, 1987, and FORD & WILLIAMS, 2007. 2.1. Measurement of hydraulic head and gradients The most direct way to measure local hydrologic conditions is to measure heads in water wells and to determine their spatial and temporal variation. Contouring of static water levels can give a good approximation of drainage divides and fl ow patterns, if the data are abundant enough and the geology is relatively simple (e.g., QUINLAN & RAY, 1989). Natural variations in head with time, provide much in- formation on proximity to turbulent-fl ow conduits, as well as fl ow convergence and divergence at different stages. Rap id rise and fall of head in response to rainfall or snow-melt events facilitates the delineation of major fl ow routes, with turbulent-fl ow solution conduits having the most extreme rates of response. Laminar-fl ow zones within the aquifer re- spond more slowly, with broad, long-term rises and declines in head. Volume of groundwater storage and rates of release can be estimated. Variations in head with time can demon- strate the dispersal of water during fl ood pulses. As in other aquifers, pumping tests in karst provide con- siderable information on fl ow patterns, effective hydraulic conductivity, and anisotropy. Spatial variability tends to be much greater than in other aquifer types. Flow directions and anisotropy can also be anticipated from the geologic set- ting. Hydraulic gradient is a critical component in all fl ow equations, whether laminar or turbulent. Discharge and ve- locity are directly proportional to gradient in laminar fl ow but are related more or less to the square root of the gradient in turbulent fl ow. 2.2. Discharge measurements Combining discharge measurements with head data provides insight into the structural nature of the aquifer. Discharge at springs and in cave passages is proportional to catchment area. It is more useful in determining aquifer properties to measure variations in discharge with time, especially in re- sponse to storms and/or snowmelt. Differing responses of springs to what appear to be similar storm events can usu- ally be accounted for by evapotranspiration loss and ante- cedent soil-moisture conditions, both of which show a strong seasonal effect. Where the hydrology is simple, as where ponors feed a single spring, the comparison of discharge in vs. discharge out is a simple method for assessing the pres- Palmer: Understanding the hydrology of karst Geologia Croatica 145 tive enrichment in the soil by calcite precipitation. A variety of de-icers on highways (NaCl, CaCl2, etc.) can be tracked and used as a proxy for contaminant fl ow paths from acci- dental spills. 2.6. Equilibrium chemistry The degree of saturation of dissolved minerals in wells, springs, and caves helps to determine the fl ow history of the water and the nature of the openings through which it has passed. Drip-waters in caves can also provide estimates of PCO2 in the overlying soil. Samples of infi ltration water must be col- lected where it fi rst emerges from the rock or sediment ma- trix so there is no loss of dissolved gases. The pH is mea- sured in situ, and precipitation of dissolved solids is pre- vented by HCl acidifi cation of one of a pair of samples. Most seepage water is near saturation with calcite. Saturation with dolomite generally requires lengthy residence time in an aquifer (typically years). Signifi cant supersaturation with ei- ther mineral generally indicates loss of CO2 from the water, e.g., by air exchange through open cave entrances or fi s- sures. Dripwater chemistry in caves, combined with discharge of the drips, can clarify the nature of vadose fl ow routes. Is the water seeping slowly, with large areas of contact with fi ssure walls? Is it following discrete solution channels? Some of this information can be obtained from rate equa- tions (e.g., PLUMMER et al., 1978, DREYBRODT, 1988). A generic rate equation for dissolution is dC/dt = (A’ k / V) (1 – C/Cs)n (1) where dC/dt = change in concentration with time, k = reac- tion coeffi cient, n = reaction order, Cs = saturation concen- tration (and therefore C/Cs = degree of saturation, where 1.0 = saturation), A’ = surface area in contact with the water, and V = water volume. It is rarely possible to fi nd a unique solu- tion to this equation, but the range of possibilities becomes clear. An integral solution, as well as values for n and k for typical karst situations, are given by PALMER (1991). Note that t/V is dimensionally equivalent to Q, and A’ is proportional to fl ow length, so C/Cs relates directly to discharge and to the network of openings through which the water has passed. A PCO2 signifi cantly greater than that of the soil suggests a hypogenic CO2 source. For example, at Saratoga Springs, New York, water is actively degassing as it emerges from the springs. The fl ow is discontinuous, with spurts and bub- bles. With the assumption that the water at depth was at equi- librium with calcite, back-calculation indicates that its PCO2 was more than 6 atmospheres. Reaction-path software is use- ful in the calculations (seewww.usgs.gov/software). PCO2 signifi cantly less than atmospheric is common but diffi cult to measure. Seepage into carbonate rock through an insoluble cap-rock allows dissolution to take place under nearly closed conditions, because the reaction is isolated from the reservoir of soil CO2. PCO2 can drop to as little as 10–4 to 10–5 atm. This is diffi cult to measure, because local fl ow rates are very small, and CO2 is absorbed rapidly from may complicate the analysis, but they can also be revealed by this method. Approximate conduit size and water volume within it can also be crudely estimated. It is useful to know the approximate percentage of a con- duit that is vadose. Low-density contaminants that are mainly insoluble in water will tend to fl oat and accumulate at sumps. During fl ood pulses, much of a normally vadose conduit will fi ll with water and drive the contaminants upward into many narrow fi ssures. Leakage of volatile gases to the surface has been documented (e.g., in Bowling Green, Kentucky; CRAW- FORD, 1984), where volatile gases have accumulated in buildings to near explosive limits. Even if the entire inter- pretive procedure is not feasible, it is helpful merely to an- ticipate the fi eld conditions that lead to this problem. 2.4. Use of tagged tracers The ability of pathogens to be carried by groundwater can be estimated by tracing with solid particles of discrete size. An ideal tracer consists of natural indigenous microbes that have been tagged with an identifi er that is incorporated into the living cells. BRAHANA, (2009) describes this technique in a karst basin in Arkansas that is isolated for hydrologic research. The bacteria are exposed to a europium solution and reintroduced into sinking streams and other infi ltration points. The tagged bacteria are transported, temporarily de- posited, and re-suspended during fl ood pulses in the karst conduits. Their movement through the system is monitored in caves, wells, and springs. Such studies can determine the fl ow conditions in which they are most prevalent at the sam- pling points. Use of local microbes gives a more realistic re- sponse in terms of transmission, retention, and delay than artifi cial particles. 2.5. Analysis of dissolved solids Measurement of dissolved solids at different fl ow stages can clarify fl ow paths where a variety of porous media are pres- ent. Interpretations are fairly straightforward and are not elaborated here. For example, high values of sulphates and chlorides suggest deep fl ow components (most abundant during low fl ow); and Mg/Ca ratios >1.0 suggest evapora- Fi gu re 1: Variation in tracer breakthrough time vs. discharge in (a) a 1000- m-long conduit 1 m in diameter; (b) a vadose canyon 1 m wide and 1000 m long, with slope = 0.01 and Manning friction factor = 0.05; and (c) a 500-m- long cayon as in b, leading to a 500-m-long tube as in a. Geologia Croatica 63/2Geologia Croatica 146 the cave air. Results include rills, etching, and weathering of cave walls where seepage water suddenly becomes highly aggressive toward carbonates (PALMER, 2007). Discharge per unit area is low, but the total over large areas can be sub- stantial. Infi ltration through insoluble rock can be verifi ed, and contaminant paths and fi ltering can be anticipated. 2.7. Stable-isotope chemistry Flow rates and storage in non-conduit parts of a karst system can be assessed through the use of isotope chemistry. 18O and deuterium (2H) in precipitation vary seasonally, and when plotted on a graph of δ18O vs. δ2H they tend to fall on a straight line. This “local meteoric water line” is typically close to, and nearly parallel to, that of the global mean data, defi ned by δ2H = 10 + 8 [δ18O], as shown in Fig. 2. During the cold season, the points plot at the lower left (more de- pleted in the heavier isotopes); and during the warm season they plot at the upper right. As this water infi ltrates and trav- els through the aquifer, it retains most of its initial isotopic signature. Changes due to CO2 uptake and dissolution of car- bonates have little effect, because their total molar contribu- tion to the water is only a tiny fraction. The discrepancy in isotopic signature between local precipitation and of a ground water sample increases with time for the fi rst half year, and then diminishes. This discrepancy forms a sinusoi- dal curve with time (Fig. 3), and therefore the “age” (resi- dence time) of the water since it fi rst precipitated can be es- timated in most cases. A single measurement is not suffi cient; a time series needs to be developed that shows the seasonal variation over a long time (ideally several years). A ground- water sample that tracks almost exactly with that of the pre- cipitation is likely to have spent less than a single season reaching its present point in the aquifer. Large discrepancies indicate longer residence times in the ground. An erratic cor- relation requires more care in interpretation. Comparison of values in a variety of places in the aquifer is likely to clarify the overall picture of residence times. Groundwater that is supplied by a wide range of water sources, typical of bathyphreatic fl ow, tends to have a clus- tering of oxygen-deuterium values somewhat low on the me- teoric water line, regardless of time of year (Fig. 1). This is because most of the infi ltration takes place in the transition from cold to warm, when snow melt is high and evapotran- spiration is low. Points to the right of the local meteoric water line nor- mally indicate evaporative enrichment of the heavy isotopes at the surface. (It may also indicate long-term inheritance from local bedrock, but this appears to be minor in karst aq- uifers.) For example, McFail’s Cave, New York receives many drips along its length of about 11 km. Some of the land above it is swampy, so there is great opportunity for evapo- ration. Some of the drips show a signifi cant shift to the right, away from the meteoric water line, and these are interpreted as being fed by wetlands. The position of the drips relative to the wetlands supports this idea, but there are local discrep- ancies that appear to indicate lateral movement of vadose water down the dip of the strata. Sulphate content can be used as an indicator of depth of groundwater fl ow. For example, in the same fi eld area there are sulphate-rich springs that show a sulfur isotopic signa- ture (δ34S) that is identical to that of the solid sulphates in a local Silurian dolomite, which in this area lies at a depth of up to 100–200 m below the surface. This same water shows a clustering of oxygen-deuterium values that does not vary signifi cantly with time. The δ34S values plot exactly in the range of values for Silurian marine sulphate rocks. 2.8. Use of radioactive isotopes: the example of radium One of the most promising avenues for investigating karst hydrology is the interpretation of radium (Ra) isotopes. Lit- tle has been done in this fi eld, and certain details have yet to be understood. Bedrock contains trace amounts of uranium and thorium: 238U, 235U, and 232Th. Decay of 238U produces 226Ra (half-life of 1601 yr), and 235U produces 223Ra (half-life of 11.1 days). 232Th produces both 228Ra (half-life of 5.7 yr) and 224Ra (half-life of 3.54 days). All Ra isotopes decay to radon (Rn) and eventually to lead (Pb). Water in contact with bedrock or sediment acquires these isotopes by dissolution. The main governing variables are shown by Equation 1. With four radium isotopes of varied sources and half-lives, they provide a potentially great amount of information about fl ow patterns in karst. For gen- eral reference, see KRAEMER & GENEREUX, 1998. Those with long half-lives (years) are 226Ra and 228Ra. Those with short half-lives (days) are 223Ra and 224Ra. As water travels through narrow fi ssures or fi ne-grained material with large areas of contact and long residence times, it equilibrates with the local bedrock and sediment values. The ratio 235U/238U is about 0.0466, so that is the equilibrium ratio for the decay products 223Ra/226Ra. But 224Ra and 228Ra come from the same parent (232Th), so 224Ra/228Ra ideally approaches 1.0. Fi gu re 2: Oxygen-deuterium graph in the eastern New York karst. A = ear- ly spring snowmelt; B = rapid cave drips, late spring; C = cave streams, late spring; D = rapid cave drips, late summer; E = springs fed by deep ground- water, year-round; F = springs fed by ponors, late spring; G = cave drips beneath 60 m of low-permeability cover, fall; H = surface lake, summer; J = cave drips beneath 20 m of high-permeability cover, fall; L (to right of meteoric water line) = cave drip fed by overlying swamp, fall; M = springs fed by ponors, fall; N = rainfall, summer. Data from TERRELL et al. (2005), SIEMION (2006) and the author. Palmer: Understanding the hydrology of karst Geologia Croatica 147 High activities of the long-lived Ra isotopes indicate lengthy and intimate contact with bedrock or sediment. High activities of the short-lived Ra (224Ra and 223Ra) indicate re- cent release of water from pores and narrow fi ssures into larger bodies of water, where they decay faster than they can be replenished from the solids. Zones of seepage into streams (at the surface or underground) can be identifi ed in this way. The time since the seepage has entered the main body of wa- ter can be estimated by the amount of disequilibrium. For example, with time, the 224Ra/228Ra ratio decreases. This de- crease can be caused by simple mixing between the incom- ing seepage and the stream or lake water; but if the seepage input retains its identity as a plume without substantial mix- ing (as in lakes fed by springs), the ratio decreases along an exponential decay curve instead of linearly. Marine rocks (e.g., carbonates) have a low Th/U activ- ity ratio (<0.1), so they produce low 228Ra/226Ra ratios in groundwater. Siliciclastic rocks have a higher Th/U ratio (about 1.0), so 228Ra/226Ra also approaches 1.0, but slowly, and most water will have completed its tour of the aquifer long before it reaches equilibrium. Individual strata can have their own 228Ra/226Ra identity. For example, shaly carbonates tend to have a higher 228Ra/226Ra than more pure carbonates. Laboratory analysis of rock samples can clarify the expected ratios. In the Madison carbonate aquifer of South Dakota, recharge through overlying siliciclastic beds can be distin- guished from recharge directly into the aquifer by their con- trast in 228Ra and 226Ra, even where there is no difference in oxygen-deuterium content (KRAEMER & GENEREUX, 1998). Water can also acquire Ra isotopes by alpha recoil: re- lease of an alpha particle from a radioactive nucleus causes the source particle to recoil up to several hundred nanome- tres. If the source is located near the surface of a solid, it can be ejected into the adjacent fl uid. Alpha recoil has a signifi - cant effect only where there is a very large surface area, as in suspended fi ne-grained sediment. High 228Ra/226Ra ratios can result: up to at least 6 in slow-moving groundwater. This suggests both considerable alpha-recoil and also long resi- dence times (tens or hundreds of years; Tom KRAEMER, Reston, Virginia, personal communication). Residence time of groundwater is often estimated with tritium, radiocarbon methods, etc., but the utility of those produced by human activity (e.g., tritium from atmospheric testing of nuclear devices) is diminishing with time. 3. CONCLUSIONS The interpretive techniques described above are all interre- lated by several major concepts: the mass balance, fl uid me- chanics, dissolution kinetics, and chemical equilibrium. They share the same physical laws, concepts, and variables as those that govern the origin of karst. All require a funda- mental understanding of the local geological setting. By ap- plying these and similar techniques, one can understand the origin and development of karst in addition to how it be- haves hydrologically. These techniques require considerable time and expense, so they are feasible only for long-term academic or govern- ment-sponsored research projects. It is important to develop suitable strategies for contaminant remediation by anticipat- ing problems, rather than by reacting after a spill or leakage has already taken place. Today the most common approach to groundwater man- agement is digital modelling, in which the goal is to predict outcomes on the basis of a few simple fi eld measurements. This is rarely successful. Instead, it is more appropriate to design simple digital models that are not expected to be cor- rect, but the output of which can easily be compared to fi eld observations. Confronted with the differences between the ideal model and reality, we are forced to visualize the fi eld conditions that can account for those differences. Another approach to digital modelling is to develop interactive mod- els of specifi c aspects of karst that allow one to explore the effects of varied fi eld conditions (e.g., DREYBRODT et al., 2005). These methods encourage professional growth, in contrast to the typical accumulation of digital outputs that involve only repetitive data types and thought processes. A recommended strategy is fi rst to use fi eld techniques such as those described in this paper to gain insight into the behaviour of karst aquifers, and only then to progress to aq- uifer modelling, whether it be conceptual, analytical, statis- tical, or digital. Field information should ideally be held in a central repository, (e.g., in a government-sponsored or ac- ademic karst centre), rather than scattered widely in the lit- erature or in personal fi les. Hydrologic and chemical fi eld measurements give personal insight into the internal work- ings of karst – an important step in the training and advance- ment of karst scientists. ACKNOWLEDGMENT Many thanks to Margaret PALMER for fi eld and laboratory assistance, Tom KRAEMER for interpretive insight into radium geochemistry and graduate students Elias MOSKAL, Levia TERRELL and Jason SIE- MION for collaboration in the fi eld. REFERENCES ATKINSON, T.C. (1977): Diffuse fl ow and conduit fl ow in limestone terrain in the Mendip Hills, Somerset (Great Britain).– Journal of Hydrology, 35, 93–110. BAUČIĆ, I. (1968): Subterranean connections in the hydrological drain- age basin of the Cetina River.– Proceedings of 4th International Congress of Speleology, Ljubljana, 271–277. BONACCI, O. (1987): Karst Hydrology.– Springer-Verlag, Berlin, 184 p. BRAHANA, V. (2009) The Savoy Experimental Watershed, Arkansas: a long-term research site for karst hydrogeology.– In: PALMER, A.N. & PALMER, M.V. (eds.): Caves and Karst of the USA. Na- tional Speleological Society, Huntsville, Alabama, USA, 179 p. CRAWFORD, N.C. (1984): Toxic and explosive fumes rising from car- bonate aquifers: a hazard for residents of sinkhole plains.– In: BECK, B.F. (ed.): Sinkholes: Their Geology, Engineering, and En- vironmental Impact. Balkema, Rotterdam, 297–304. DREYBRODT, W. (1988): Processes in Karst Systems: Physics, Chem- istry, and Geology.– Springer-Verlag, Berlin, 288 p. Geologia Croatica 63/2Geologia Croatica 148 DREYBRODT, W., GABROVŠEK, F. & ROMANOV, D. (2005): Pro- cesses of Speleogenesis: A Modeling Approach: Carstologia, Ljublja- na, ZRC Publishing, 376 p. + CD FORD, D.C. & WILLIAMS, P.W. (2007): Karst hydrology and geomor- phology.– Wiley, Chichester, UK, 601 p. KÄSS, W. (1998): Tracing technique in geohydrology.– Balkema, Rotter dam, 581 p. KRAEMER, T.F. & GENEREUX, D.P. (1998): Applications of uranium- and thorium-series radionuclides in catchment hydrology studies.– In: KENDALL, C. & McDONNELL, J.J. (eds.): Isotope Tracers in Catchment Hydrology. Elsevier, Amsterdam, 679–722. MILANOVIĆ, P.T. (1981): Karst Hydrogeology.– Water Resources Pub- lications, Littleton, Colorado, 434 p. PALMER, A.N. (1991): Origin and morphology of limestone caves.– Geological Society of America Bulletin, 105/1, 1–21. PALMER, A.N. (2007): Cave Geology.– Cave Books, Dayton, Ohio, 454 p. PITTY, A.F. (1966): An Approach to the Study of Karst Water.– Univer- sity of Hull, U.K., Occasional Papers in Geography, 5, 70 p. PLUMMER, L.N., WIGLEY, T.M.L. & PARKHURST, D.L. (1978): The kinetics of calcite dissolution in CO2-water systems at 5–60°C and 0.0–1.0 atm CO2.– American Journal of Science, 278, 179–216. QUINLAN, J.F. & RAY, J.A. (1989): Groundwater basins in the Mam- moth Cave region, Kentucky.– In: WHITE, W.B. & WHITE, E.L. (eds.): Karst Hydrdology: Concepts from the Mammoth Cave Re- gion. Van Nostrand Reinhold, New York, one-page map insert. SIEMION, J. (2006): Use of isotopes, tracer tests, and dissolved solids to characterize a karst aquifer.– M.A. thesis, State University of New York, Oneonta, NY, 49 p. TERRELL, L., PALMER, A. & KRAEMER, T. (2008): Isotopic and chemical analysis of sulfate-rich springs along the Helderberg Escarp ment of east-central New York State.– Geological Society of America, Northeastern Section, annual meeting, abstracts with programs, 37/1, 69. Manuscript received January 04, 2010 Revised manuscript accepted March 19, 2010 Available online May 31, 2010