Springflow hydrographs: eogenetic vs. telogenetic karst

Matrix permeability in the range of 10(-11) to 10(-14) m(2) characterizes eogenetic karst, where limestones have not been deeply buried. In contrast, limestones of postburial, telogenetic karst have matrix permeabilities on the order of 10(-15) to 10(-20) m(2). Is this difference in matrix permeability paralleled by a difference in the behavior of springs draining eogenetic and telogenetic karst? Log Q/Q(min) flow duration curves from 11 eogenetic-karst springs in Florida and 12 telogenetic-karst springs in Missouri, Kentucky, and Switzerland, plot in different fields because of the disparate slopes of the curves. The substantially lower flow variability in eogenetic-karst springs, which results in the steeper slopes of their flow duration curves, also makes for a strong contrast in patterns (e.g., "flashiness") between the eogenetic-karst and telogenetic-karst spring hydrographs. With respect to both spring hydrographs and the flow duration curves derived from them, the eogenetic-karst springs of Florida are more like basalt springs of Idaho than the telogenetic-karst springs of the study. From time-series analyses on discharge records for 31 springs and published time-series results for 28 additional sites spanning 11 countries, we conclude that (1) the ratio of maximum to mean (Q(max)/Q(mean)) discharge is less in springs of eogenetic karst than springs of telogenetic karst; (2) aquifer inertia (system memory) is larger in eogenetic karst; (3) eogenetic-karst aquifers take longer to respond to input signals; and (4) high-frequency events affect discharge less in eogenetic karst. All four of these results are consistent with the hypothesis that accessible storage is larger in eogenetic-karst aquifers than in telogenetic-karst aquifers.     

Vadose CO2 gas drives dissolution at water tables in eogenetic karst aquifers more than mixing dissolution

Most models of cave formation in limestone that remains near its depositional environment and has not been deeply buried (i.e. eogenetic limestone) invoke dissolution from mixing of waters that have different ionic strengths or have equilibrated with calcite at different _pCO₂ values. In eogenetic karst aquifers lacking saline water, mixing of vadose and phreatic waters is thought to form caves. We show here calcite dissolution in a cave in eogenetic limestone occurred due to increases in vadose CO₂ gas concentrations and subsequent dissolution of CO₂ into groundwater, not by mixing dissolution. We collected high‐resolution time series measurements (1 year) of specific conductivity (SpC), temperature, meteorological data, and synoptic water chemical composition from a water table cave in central Florida (Briar Cave). We found SpC, _pCO₂ and calcite undersaturation increased through late summer, when Briar Cave experienced little ventilation by outside air, and decreased through winter, when increased ventilation lowered cave CO_2(g) concentrations. We hypothesize dissolution occurred when water flowed from aquifer regions with low _pCO₂ into the cave, which had elevated _pCO₂. Elevated _pCO₂ would be promoted by fractures connecting the soil to the water table. Simple geochemical models demonstrate that changes in _pCO₂ of less than 1% along flow paths are an order of magnitude more efficient at dissolving limestone than mixing of vadose and phreatic water. We conclude that spatially or temporally variable vadose CO_2(g) concentrations are responsible for cave formation because mixing is too slow to generate observed cave sizes in the time available for formation. While this study emphasized dissolution, gas exchange between the atmosphere and karst aquifer vadose zones that is facilitated by conduits likely exerts important controls on other geochemical processes in limestone critical zones by transporting oxygen deep into vadose zones, creating redox boundaries that would not exist in the absence of caves. Copyright © 2014 John Wiley & Sons, Ltd.     

Infiltration Measured by the Drip of Stalactites

The hydrodynamic processes and mechanisms involved in rain infiltration and recharge in local areas of karst terrain can be identified and quantified by using measurements of the seepage of cave stalactites. Detailed measurements of the seepage of stalactites in seven caves located in an area close to the land surface, or the subcutaneous area of the karst, show a diversity of complex factors involved in infiltration: type of precipitation (rain or snow), air temperature, soil type and thickness, etc., which give rise to larger or smaller variations of flow in the espeleothem hydrographs. In some cases, no explanation can be found for the response of stalactites to rainfall, while in others there is a relationship between outer atmospheric parameters and the recharge represented by the stalactite drip. Romperopas Cave (Spain) has both a rapid and a basic flow, with hydrograph recessions similar to those observed in other caves. Water seepage in this cave varies greatly both in space and in time. The infiltration in Altamira Cave (Spain) was calculated and a multiple regression was found between infiltration, rain and outside air temperature. In other cases, the balance of the water on the soil is responsible for the seepage. Thus, a precipitation runoff numerical model that simulated the stalactite hydrographs could be applied to the Baradla and Beké Caves (Hungary). The complex properties of the ground, which are required for other flow numerical models for the unsaturated zone, were not taken into consideration.     

Formation of phreatic caves in an eogenetic karst aquifer by CO2 enrichment at lower water tables and subsequent flooding by sea level rise

Formation of extensive phreatic caves in eogenetic karst aquifers is widely believed to require mixing of fresh and saltwater. Extensive phreatic caves also occur, however, in eogenetic karst aquifers where fresh and saltwater do not mix, for example in the upper Floridan aquifer. These caves are thought to have formed in their modern settings by dissolution from sinking streams or by convergence of groundwater flow paths on springs. Alternatively, these caves have been hypothesized to have formed at lower water tables during sea level low‐stands. These hypotheses have not previously been tested against one another. Analyzing morphological data and water chemistry from caves in the Suwannee River Basin in north‐central Florida and water chemistry from wells in the central Florida carbonate platform indicates that phreatic caves within the Suwannee River Basin most likely formed at lower water tables during lower sea levels. Consideration of the hydrological and geochemical constraints posed by the upper Floridan aquifer leads to the conclusion that cave formation was most likely driven by dissolution of vadose CO₂ gas into the groundwater. Sea level rise and a wetter climate during the mid‐Holocene lifted the water table above the elevation of the caves and placed the caves tens of meters below the modern water table. When rising water tables reached the land surface, surface streams formed. Incision of surface streams breached the pre‐existing caves to form modern springs, which provide access to the phreatic caves. Phreatic caves in the Suwannee River Basin are thus relict and have no causal relationship with modern surficial drainage systems. Neither mixing dissolution nor sinking streams are necessary to form laterally extensive phreatic caves in eogenetic karst aquifers. Dissolution at water tables, potentially driven by vadose CO₂ gas, offers an underappreciated mechanism to form cavernous porosity in eogenetic carbonate rocks. Copyright © 2012 John Wiley & Sons, Ltd.     

Integrated cave drip monitoring for epikarst recharge estimation in a dry Mediterranean area,Sif Cave,Israel

Understanding recharge mechanisms and controls in karst regions is extremely important for managing water resources because of the dynamic nature of the system. The objective of this study was to evaluate water percolation through epikarst by monitoring water flow into a cave and conducting artificial irrigation and tracer experiments, at Sif Cave in Wadi Sussi, Israel from 2005 through 2007. The research is based on continuous high‐resolution direct measurements of both rainfall and water percolation in the cave chamber collected by three large PVC sheets which integrate drips from three different areas (17, 46, and 52 m²). Barrels equipped with pressure transducers record drip rate and volume for each of the three areas. The combined measured rainfall and cave data enables estimation of recharge into the epikarst and to better understand the relationship of rainfall‐recharge. Three distinct types of flow regimes were identified: (1) ‘Quick flow’ through preferential flow paths (large fractures and conduits); (2) ‘Intermediate flow’ through a secondary crack system; and (3) ‘Slow flow’ through the matrix. A threshold of ～100 mm of rain at the beginning of the rainy season is required to increase soil water content allowing later rainfall events to percolate deeper through the soil and to initiate dripping in the cave. During winter, as the soil water content rises, the lag time between a rain event and cave drip response decreases. Annual recharge (140–160 mm in different areas in the cave) measured represents 30–35% of annual rainfall (460 mm). Copyright © 2011 John Wiley & Sons, Ltd.     

Water percolation process studies in a Mediterranean karst area

To investigate processes of water percolation, the drip response of stalactites in a karstic cave below a 143 m² sprinkling plot was measured. The experiment was conducted in Mount Carmel, Israel, at the end of the dry season and intended to simulate a series of two high‐intensity storms on dry and wet soils. In addition to hydrometric measurements (soil moisture, surface runoff, stalactite dripping rates), two types of tracers (electrical conductivity and bromide) were used to study recharge processes, water origin and mixing inside a 28‐m vadose zone. Results suggested that slow, continuous percolation through the rock matrix is of minor importance and that percolating water follows a complicated pattern including vertical and horizontal flow directions. While bromide tracing allowed identification of quick direct flow paths at all drips with maximum flow velocities of 4·3 m/h, mixing analysis suggested that major water fractions were mobilized by piston flow, pushing out water stored in the unsaturated zone above the cave. Under dry preconditions, 80 mm of artificial rainfall applied in less than 7 h was not enough to initiate significant downward water percolation. Most water was required to fill uppermost soil and rock storages. Under wet preconditions during the second day sprinkling, higher water contents in soils and karst cavities facilitated piston flow effects and a more intense response of the cave drips. Results indicate that in Mediterranean karst regions, filling of the unsaturated zone, including soil and rock storages, is an important precondition for the onset of significant water percolation and recharge. This results in a higher seasonal threshold for water percolation than for the generation of surface runoff. Copyright © 2010 John Wiley & Sons, Ltd.     

Modeling the Usefulness of Spatial Correlation Analysis on Karst Systems

Cross-correlation analyses on field data collected in karst aquifer systems can be used to develop a conceptual understanding of the aquifer. This includes the use of many data sets from the same aquifer to develop an understanding of how properties vary spatially. We focus on a method for characterizing the distribution of recharge, which is becoming increasingly important in regions where urban development encroaches on these important sources of water. Spatially varying precipitation data and cross-correlation analysis provide a means of spatially characterizing recharge locations on a karst aquifer. Our work expands on the numerical experiments conducted by Padilla and Pulido-Bosch (1995) using the numerical ground water model MODFLOW to introduce spatially varying parameters. The numerical experiments include conduit-controlled, matrix-controlled, and mixed karst systems with more than one precipitation time series input. The results show that spatially varying parameters can be inferred based on the cross-correlation of precipitation data and spring discharge. Simulations were completed using aquifer parameters derived from studies of the Barton Springs segment of the Edwards Aquifer. The simulations indicate that spatial variability within an aquifer can be inferred using cross-correlation analysis. A field study using these methods is summarized for Barton Springs near Austin, Texas.     

Evidence of karst from behaviour of the Asmari limestone aquifer at the Khersan 3 Dam site,southern Iran

Abstract

The Asmari limestone formation is the major aquifer system at the Khersan 3 Dam site, Zagros, Iran. Characterization of the aquifer system and study of karst development are essential for forecasting leakage potential and to plan remediation works. The aquifer functioning and karst structure were evaluated by geology, well hydrodynamics and natural tracing studies, showing one unconfined and three artesian sub-aquifers, the last being characterized by rapid flow, with contributing old recharge water, and a recharge at higher elevation than the unconfined sub-aquifer. The anticline structure outcrops the artesian sub-aquifers downstream of the site. The confining layers disappear and the aquifer system discharges as a mix of all groundwater. Artesian groundwater is estimated to contribute about 80% (dry season) and 50% (wet season) of flow in the mixing zone. Very few karst features are observed in cores and galleries at the site, despite some karst landforms in the reservoir area. However the aquifer exhibits fast flow of karst type in the artesian sub-aquifers.     

Using nitrate to quantify quick flow in a karst aquifer

In karst aquifers, contaminated recharge can degrade spring water quality, but quantifying the rapid recharge (quick flow) component of spring flow is challenging because of its temporal variability. Here, we investigate the use of nitrate in a two-endmember mixing model to quantify quick flow in Barton Springs, Austin, Texas. Historical nitrate data from recharging creeks and Barton Springs were evaluated to determine a representative nitrate concentration for the aquifer water endmember (1.5 mg/L) and the quick flow endmember (0.17 mg/L for nonstormflow conditions and 0.25 mg/L for stormflow conditions). Under nonstormflow conditions for 1990 to 2005, model results indicated that quick flow contributed from 0% to 55% of spring flow. The nitrate-based two-endmember model was applied to the response of Barton Springs to a storm and results compared to those produced using the same model with δ¹⁸O and specific conductance (SC) as tracers. Additionally, the mixing model was modified to allow endmember quick flow values to vary over time. Of the three tracers, nitrate appears to be the most advantageous because it is conservative and because the difference between the concentrations in the two endmembers is large relative to their variance. The δ¹⁸O-based model was very sensitive to variability within the quick flow endmember, and SC was not conservative over the timescale of the storm response. We conclude that a nitrate-based two-endmember mixing model might provide a useful approach for quantifying the temporally variable quick flow component of spring flow in some karst systems.     

Isotopic variation in groundwater across the conterminous United States – Insight into hydrologic processes

Groundwater supplies a significant proportion of water use in the United States and is critical to the maintenance of healthy ecosystems and environmental processes, thus characterizing aquifer hydrology is important to managing and preserving these resources. While groundwater isotopes provide insight into hydrologic and ecologic processes, their application is limited to where measurements exist. To help overcome this limitation, we used the random forest algorithm to develop a predictive model for shallow groundwater isotopes in the conterminous United States. Our model uses environmental variables (e.g. temperature, elevation, precipitation isotopes) as predictors. We used our model to develop the first shallow groundwater isoscape of δ²H and δ¹⁸O for the conterminous United States. We describe the patterns in groundwater isotopes using both observations and our modelled isoscape. We find that throughout much of the Eastern United States, groundwater isotopes are close to annual amount weighted precipitation, while groundwater isotopes are significantly depleted relative precipitation across much of the High Plains and Western United States. Furthermore, we compare the observations compiled for this study to isotopes of precipitation, which allows us to determine the relative recharge efficiency (i.e. ratio of groundwater recharge to precipitation) between seasons and the proportion of annual recharge that occurs in a given season. Our findings suggest that winter recharge is generally more efficient than summer recharge; however, the dominant recharge season is more varied as it is the product of both seasonal recharge efficiency and the seasonal timing of precipitation. Parts of the central United States have summer dominant recharge, which is likely the result of heavy summer precipitation/nocturnal summer precipitation. Interestingly, parts of coastal California appear to have summer dominant recharge, which we suggest could be due to recharge from fog-drip. Our results summarize spatial patterns in groundwater isotopes across the conterminous United States, provide insight into the hydrologic processes affecting shallow groundwater, and are valuable information for future ecologic and hydrologic studies.     

Classification of Thermal Patterns at Karst Springs and Cave Streams

Thermal patterns of karst springs and cave streams provide potentially useful information concerning aquifer geometry and recharge. Temperature monitoring at 25 springs and cave streams in southeastern Minnesota has shown four distinct thermal patterns. These patterns can be divided into two types: those produced by flow paths with ineffective heat exchange, such as conduits, and those produced by flow paths with effective heat exchange, such as small fractures and pore space. Thermally ineffective patterns result when water flows through the aquifer before it can equilibrate to the rock temperature. Thermally ineffective patterns can be either event‐scale, as produced by rainfall or snowmelt events, or seasonal scale, as produced by input from a perennial surface stream. Thermally effective patterns result when water equilibrates to rock temperature, and the patterns displayed depend on whether the aquifer temperature is changing over time. Shallow aquifers with seasonally varying temperatures display a phase‐shifted seasonal signal, whereas deeper aquifers with constant temperatures display a stable temperature pattern. An individual aquifer may display more than one of these patterns. Since karst aquifers typically contain both thermally effective and ineffective routes, we argue that the thermal response is strongly influenced by recharge mode.     

碳酸盐岩溶洞成像要素分析与研究

碳酸盐岩油气藏储集空间以裂缝和溶洞为主,在地震剖面上主要表现为＂串珠状＂。为了研究溶洞成像规律,本文通过大量的数值模拟,分析溶洞的成像分辨率,探讨溶洞成像与上覆地层、信噪比（SNR）、覆盖次数、偏移速度和偏移方法之间的定量关系。研究认为：①无论是溶洞的纵向分辨率还是横向分辨率,均满足1/4波长调谐规律,可以将薄层机理转用至溶洞分辨率解释;②溶洞的存在会对上覆地层产生下拉效应,使溶洞在地震剖面上呈塌陷等＂非串珠状＂特征;③高覆盖和小道距（面元）可提高串珠能量,压制偏移附带的背景噪声;④叠前偏移能改善资料的SNR,提高成像精度,特别是低SNR资料;⑤偏移速度偏小,串珠边缘下拉,反之会上翘,速度偏大带来的成像质量变差的程度要大于速度偏小;⑥不同偏移方法对缝洞成像效果存在差异,对于低SNR资料,溶洞成像拟选取频率域或频波域方法。     

Heterogeneous distributions of CO2 may be more important for dissolution and karstification in coastal eogenetic limestone than mixing dissolution

Mixing dissolution, a process whereby mixtures of two waters with different chemical compositions drive undersaturation with respect to carbonate minerals, is commonly considered to form cavernous macroporosity (e.g. flank margin caves and banana holes) in eogenetic karst aquifers. On small islands, macroporosity commonly originates when focused dissolution forms globular chambers lacking entrances to the surface, suggesting that dissolution processes are decoupled from surface hydrology. Mixing dissolution has been thought to be the primary dissolution process because meteoric water would equilibrate rapidly with calcium carbonate as it infiltrates through matrix porosity and because pCO₂ was assumed to be homogeneously distributed within the phreatic zone. Here, we report data from two abandoned well fields in an eogenetic karst aquifer on San Salvador Island, Bahamas, that demonstrate pCO₂ in the phreatic zone is distributed heterogeneously. The pCO₂ varied from less than log ?2.0 to more than log ?1.0 atm over distances of less than 30 m, generating dissolution in the subsurface where water flows from regions of low to high pCO₂ and cementation where water flows from regions of high to low pCO_2. Using simple geochemical models, we show dissolution caused by heterogeneously distributed pCO₂ can dissolve 2.5 to 10 times more calcite than the maximum amount possible by mixing of freshwater and seawater. Dissolution resulting from spatial variability in pCO₂ forms isolated, globular chambers lacking initial entrances to the surface, a morphology that is characteristic of flank margin caves and banana holes, both of which have entrances that form by erosion or collapse after cave formation. Our results indicate that heterogeneous pCO₂, rather than mixing dissolution, may be the dominant mechanism for observed spatial distribution of dissolution, cementation and macroporosity generation in eogenetic karst aquifers and for landscape development in these settings. Copyright © 2015 John Wiley & Sons, Ltd.     

Conceptualization of Karstic Aquifer with Multiple Outlets Using a Dual Porosity Model

In this study, two conceptual models, the classic reservoir (CR) model and exchange reservoirs model embedded by dual porosity approach (DPR) are developed for simulation of karst aquifer functioning drained by multiple outlets. The performances of two developed models are demonstrated at a less developed karstic aquifer with three spring outlets located in Zagros Mountain in the south‐west of Iran using 22‐years of daily data. During the surface recharge, a production function based on water mass balance is implemented for computing the time series of surface recharge to the karst formations. The efficiency of both models has been assessed for simulation of daily spring discharge during the recession and also surface recharge periods. Results indicate that both CR and DPR models are capable of simulating the ordinates of spring hydrographs which drainage less developed karstic aquifer. However, the goodness of fit criteria indicates outperformance of DPR model for simulation of total hydrograph ordinates. In addition, the DPR model is capable of quantifying hydraulic properties of two hydrologically connected overlapping continua conduits network and fissure matrix which lays important foundations for the mining operation and water resource management whereas homogeneous model representations of the karstic subsurface (e.g., the CR) do not work accurately in the karstic environment.     

A New Method for Estimating Recharge to Unconfined Aquifers Using Differential River Gauging

In semiarid and arid environments, leakage from rivers is a major source of recharge to underlying unconfined aquifers. Differential river gauging is widely used to estimate the recharge. However, the methods commonly applied are limited in that the temporal resolution is event‐scale or longer. In this paper, a novel method is presented for quantifying both the total recharge volume for an event, and variation in recharge rate during an event from hydrographs recorded at the upstream and downstream ends of a river reach. The proposed method is applied to river hydrographs to illustrate the method steps and investigate recharge processes occurring in a sub‐catchment of the Murray Darling Basin (Australia). Interestingly, although it is the large flood events which are commonly assumed to be the main source of recharge to an aquifer, our analysis revealed that the smaller flow events were more important in providing recharge.     

Extreme precipitation drives groundwater recharge: the Northern High Plains Aquifer,central United States, 1950–2010

Future extreme precipitation (EP, daily rainfall amount over certain thresholds) is projected to increase with global climate change; however, its effect on groundwater recharge has not been fully explored. This study specifically investigates the spatiotemporal dynamics of groundwater recharge and the effects of extreme precipitation (daily rainfall amount over the 95th percentile, which is tagged by ranking the percentiles in each season for a base period) on groundwater recharge from 1950 to 2010 over the Northern High Plains (NHP) Aquifer using the Soil Water Balance Model. The results show that groundwater recharge significantly (p < 0.05) increased in the eastern NHP from 1950 to 2010, where the highest annual average groundwater recharge occurs compared to the central and the western NHP. In the eastern NHP, 45.1% of the annual precipitation fell as EP, which contributed 56.8% of the annual total groundwater recharge. In the western NHP, 30.9% of the annual precipitation fell as extreme precipitation, which contributed 62.5% of the annual total groundwater recharge. In addition, recharge by extreme precipitation mainly occurred in late spring and early summer, before the maximum evapotranspiration rate, which usually occurs in mid‐summer until late fall. A dry site in the western NHP and a wet site in the eastern NHP were analysed to indicate how recharge responds to EP with different precipitation regimes. The maximum daily recharge at the dry site exceeded the wet site when there was EP. When precipitation fell as non‐extreme rainfall, most recharge was less than 5 mm at both the dry and wet sites, and the maximum recharge at the dry site became lower than the wet site. This study shows that extreme precipitation plays a significant role in determining groundwater recharge. © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

Seasonal δ13CDIC sourcing and geochemical flux in telogenetic epikarst of south-central Kentucky

Telogenetic epikarst carbon sourcing and transport processes and their associated hydrogeochemical responses are complex and dynamic. Carbon dioxide (CO₂) transport rates in the epikarst zone are often driven by hydrogeochemical responses, which influence carbonate dissolution and conduit formation. This study examines the influence of land use on carbon sourcing and carbonate dissolution kinetics through a comparative analysis of separate, but similar, epikarst systems in south-central Kentucky. The use of high-resolution hydrogeochemical data from multiple data loggers and isotope analysis from collected water samples reflects the processes within these epikarst aquifers, which are estimated to contribute significantly to bedrock dissolution. Results indicate that, in an agricultural setting, long-term variability and dissolution is governed by seasonal production of CO₂ . In a more urbanized, shallower epikarst system, land cover may affect CO₂ transport between the soil and underlying bedrock. This concentration of CO₂ potentially contributes to ongoing dissolution and conduit development, irrespective of seasonality. The observed responses in telogenetic epikarst systems seem to be more similar to eogenetic settings, which is suggested to be driven by CO₂ transport occurring independent of high matrix porosity. The results of this study indicate site-specific responses with respect to both geochemical and δ¹³C_DIC changes on a seasonal scale, despite regional geologic similarities. The results indicate that further comparative analyses between rural and urban landscapes in other karst settings is needed to delineate the impact of land use and seasonality on dissolution and carbon sourcing during karst formation processes. © 2019 John Wiley & Sons, Ltd.     

Hydrology of cave drip waters at varying bedrock depths from a karst system in southeastern Australia

The study of how cave drip‐water discharge responds to recharge events is fundamental to evaluating the potential of actively forming speleothems as high‐resolution climate archives. Most previous research has focused on caves of the Northern Hemisphere middle latitudes, where recharge is strongly seasonal. Few studies have explored drip‐water behaviour from regions where the expected seasonal rainfall pattern is significantly perturbed on an irregular basis by changing regional atmospheric circulation patterns. Here, we report the results of a 4‐year study of cave drip‐water–climate relationships from two caves in eastern Australia. The discharge of 10 drip sites located beneath bedrock thicknesses of 12, 22 and 45 m was monitored either continuously (using automated infrared sensors) or at discrete approximately monthly intervals and compared with local rainfall and water balance data. The study period traversed two major droughts, including the severe 2002–2003 El Niño. Drips at 12 and 22 m depths responded almost simultaneously to individual recharge events, although the time lag between individual events varied according to the volume of recharge and pre‐event storage. Overall, a steady decline in discharge is evident through the moisture‐deficit period, with increased flows through phases of positive water balance. Speleothems growing at these and similar shallow‐chamber sites have potential for reconstructing palaeo‐rainfall trends at high‐resolution, although the highly variable nature of year‐to‐year recharge would make it difficult to obtain data on a calendrical time‐scale. Drips at 45 m depth did not respond consistently to individual recharge events and displayed hydrological behaviour markedly dissimilar to one another and to the near‐surface drip sites, indicating great complexity in karst architecture and the absence of fissure flow. Although speleothems at this depth may well preserve information on longer‐term rainfall trends, their potential to encode a palaeo‐rainfall variability signal at interannual resolution is poor. Copyright © 2007 John Wiley & Sons, Ltd.     

Carbon flux and landscape evolution in epigenic karst aquifers modeled from geochemical mass balance

This paper considers the contributions of epigenic karst processes as a major element of the carbon cycle and a significant agent of landscape evolution. Geochemical models developed from monitoring data and water samples are used to estimate the variation and magnitude of dissolved inorganic carbon (DIC) flux in karst landscapes at several scales, from local to global. At the local scale, the Cumberland River watershed of southeast Kentucky, these geochemical models are also used to evaluate the potential role of sulfur in the production of DIC and to compute an estimated rate of landscape erosion. Geochemical modeling using ionic species and modeled discharge reveal a variable rate of DIC flux driven by large fluctuations in calcite saturation and discharge. Ratios of reaction products and principal component analyses (PCA) suggest that some bedrock dissolution may be driven by the oxidation of reduced sulfur derived from brines entrained into the karst aquifers. Over the 3730 km² of carbonate exposure in the Cumberland River, 25.8–62.4 Gg/yr of carbon dioxide (CO₂) is conveyed from the atmosphere through the dissolution of carbonate. At the global scale, this translates to 123–296 Tg/yr of CO₂ delivered by karst processes into the aqueous system. The bedrock portion of DIC equates to a flux of 32.6 ± 2.6 m³ – 35.2 ± 2.8 m³ of bedrock during the period of study of which 29% was dolomite. This translates to a landscape erosion rate of 13.1–17.9 mm/ka in the 3.45–4.32 km² of carbonate exposure in the studied watershed. Based upon 16+ km of cave survey data spanning a vertical range of 72 to 75 m above base level, this suggests that cave development in the watershed spans the Plio‐Pleistocene. Using the modeled erosion rates, the ages of cave levels, 4.03–5.71, 3.08–4.56, 1.57–2.43, 1.01–1.67, 0.45–0.91, and < 0.45 Ma, are in good agreement with regional studies of Plio‐Pleistocene landscape evolution in the Appalachian Lowland Plateaus. Copyright © 2015 John Wiley & Sons, Ltd.