Search in KarstBase
![]() |
![]() |
Most concepts of conduit development have focused on telogenetic karst aquifers, where low matrix permeability focuses flow and dissolution along joints, fractures, and bedding planes. However, conduits also exist in eogenetic karst aquifers, despite high matrix permeability which accounts for a significant component of flow. This study investigates dissolution within a 6-km long conduit system in the eogenetic Upper Floridan aquifer of north-central Florida that begins with a continuous source of allogenic recharge at the Santa Fe River Sink and discharges from a first-magnitude spring at the Santa Fe River Rise. Three sources of water to the conduit include the allogenic recharge, diffuse recharge through epikarst, and mineralized water upwelling from depth. Results of sampling and inverse modeling using PHREEQC suggest that dissolution within the conduit is episodic, occurring only during 30% of 16 sampling times between March 2003 and April 2007. During low flow conditions, carbonate saturated water flows from the matrix to the conduit, restricting contact between undersaturated allogenic water with the conduit wall. When gradients reverse during high flow conditions, undersaturated allogenic recharge enters the matrix. During these limited periods, estimates of dissolution within the conduit suggest wall retreat averages about 4 × 10−6 m/day, in agreement with upper estimates of maximum wall retreat for telogenetic karst. Because dissolution is episodic, time-averaged dissolution rates in the sink-rise system results in a wall retreat rate of about 7 × 10−7 m/day, which is at the lower end of wall retreat for telogenetic karst. Because of the high permeability matrix, conduits in eogenetic karst thus enlarge not just at the walls of fractures or pre-existing conduits such as those in telogenetic karst, but also may produce a friable halo surrounding the conduits that may be removed by additional mechanical processes. These observations stress the importance of matrix permeability in eogenetic karst and suggest new concepts may be necessary to describe how conduits develop within these porous rocks.
Anomalous behaviour of specific electrical conductivity (SEC) was observed at a karst spring in Slovenia during 26 high-flow events in an 18-month monitoring period. A conceptual model explaining this anomalous SEC variability is presented and reproduced by numerical modelling, and the practical relevance for source protection zoning is discussed. After storm rainfall, discharge increases rapidly, which is typical for karst springs. SEC displays a first maximum during the rising limb of the spring hydrograph, followed by a minimum indicating the arrival of freshly infiltrated water, often confirmed by increased levels of total organic carbon (TOC). The anomalous behaviour starts after this SEC minimum, when SEC rises again and remains elevated during the entire high-flow period, typically 20–40 µS/cm above the baseflow value. This is explained by variable catchment boundaries: When the water level in the aquifer rises, the catchment expands, incorporating zones of groundwater with higher SEC, caused by higher unsaturated zone thickness and subtle lithologic changes. This conceptual model has been checked by numerical investigations. A generalized finite-difference model including high-conductivity cells representing the conduit network (“discrete-continuum approach”) was set up to simulate the observed behaviour of the karst system. The model reproduces the shifting groundwater divide and the nearly simultaneous increase of discharge and SEC during high-flow periods. The observed behaviour is relevant for groundwater source protection zoning, which requires reliable delineation of catchment areas. Anomalous behaviour of SEC can point to variable catchment boundaries that can be checked by tracer tests during different hydrologic conditions.
Karst aquifers represent dual flow systems consisting of a highly conductive conduit system embedded in a less permeable rock matrix. Hybrid models iteratively coupling both flow systems generally consume much time, especially because of the nonlinearity of turbulent conduit flow. To reduce calculation times compared to those of existing approaches, a new iterative equation solver for the conduit system is developed based on an approximated Newton–Raphson expression and a Gauß–Seidel or successive over-relaxation scheme with a single iteration step at the innermost level. It is implemented and tested in the research code CAVE but should be easily adaptable to similar models such as the Conduit Flow Process for MODFLOW-2005. It substantially reduces the computational effort as demonstrated by steady-state benchmark scenarios as well as by transient karst genesis simulations. Water balance errors are found to be acceptable in most of the test cases. However, the performance and accuracy may deteriorate under unfavorable conditions such as sudden, strong changes of the flow field at some stages of the karst genesis simulations.
Do you want to read the rest of this publication?
After brief descriptions of the geomorphology of the Cooleman Plain karst and in particular of the Blue Waterholes, the methods adopted to analyse the functioning of these major risings are detailed. The discharge regime of Cave Creek below them is oceanic pluvial in type perturbed by drought and snow. There is much annual variation both in seasonal incidence and total amount, with catchment efficiency correspondingly variable. Suspended sediment concentration is even more erratic and monthly determinations are inadequate for calculating corrasional denudation rates. Mean concentrations of suspended solids are about 1/18th of solute load. Total dissolved salts have a strong inverse relationship with discharge, and mean values are high compared with those for other catchments in eastern Australia but none of these determinations are from limestone catchments. Sodium, potassium, and chlorine contents are low compared with the same catchments but silica is relatively high. The ratio of alkaline earths to alkalis indicate that Cave Creek carries carbonate waters and there is an inverse regression of the ratio on discharge. There is inverse correlation of total hardness on discharge likewise due to concentration of surface waters by evaporation in dry periods, together with reduced underground solution rate at times of large, rapid flow. The spring waters remain aggressive. Close regressions of hardness on specific conductivity now permit the latter to be determined in the place of the former. Much evidence converges to indicate that all the springs at the Blue Waterholes are fed from the same conduit. The intermittent flow which comes down the North Branch on the surface to the Blue Waterholes differs significantly in many characters from the spring waters. Rates of Ca + M carbonate equivalent removal vary directly with discharge since hardness varies much less than does water volume. These gross rates have to be adjusted for (a) atmospheric salts entering the karst directly, (b) peripheral solute inputs from the non-karst two-thirds of the catchment and (c) subjacent karst solution before they can be taken as a measure of exposed karst denudation. The methods for achieving this are set out. The total corrections amount to about one third of the total hardness, though the correction for subjacent karst on its own lies within the experimental error of the investigation. The residual rate of limestone removal from the exposed karst also shows a winter/spring high rate and a summer/autumn low rate but the seasonal incidence and annual total varied very much from year to year. In comparison with results from karsts in broadly similar climate, the seasonal rhythm conforms and so does the high proportion (78%) of the solution taking place at or close to the surface. This reduces the importance of the impounded condition of this small karst but supports the use of karst denudation rate as a measure of surface lowering. Cave passage solution may however be more important in impounded karst than its absolute contribution might suggest, by promoting rapid development of underground circulation. The mean value of limestone removal is low for the climatic type and this is probably due to high evapotranspirational loss as well as to the process of eliminating atmospheric, peripheral non-karst and subjacent karst contributions. The difficulties of applying modern solution removal rate to the historical geomorphology of this karst are made evident; at the same time even crude extrapolations are shown to isolate problems valuably.
This study was undertaken to gain a better understanding of karst hydrology. To do this, the present day hydrology and the paleohydrology were determined in three karst basins. The basins chosen were the Swago, Locust and Spring Creek basins in Pocahontas and Greenbrier Counties, West Virginia. A number of conventional field techniques were used successfully in this study, including the following: current meter and dye dilution gauging; dye and lycopodium stream tracing; geological and cave mapping; the setting up of stage recorders; geochemistry; and limestone erosion measurements. The climate of the region was investigated to obtain realistic precipitation, temperature and potential evaporation data over the study basins.
It was found that the mean precipitation over two of the basins was 30% higher than recorded data in the valleys. The karst development of the basins was found to take place in four major stages. These were: A) initial surficial flow, B) strike controlled drainage, C) major piracies from one sub-basin to another, and D) shortening of the flow routes. The major controls on the karst development were found to be: A) the Taggard shale, B) the strike direction, which controlled early basin development, and C) the hydraulic gradient from the sink to rising, which controlled later basin development.
To better assess the quantitative hydrology, and to assist in determining the type of unexplorable flow paths, a watershed model was developed. This modelled the streamflow from known climatic inputs using a number of measured or optimized parameters. The simulation model handled snowmelt, interception, infiltration, interflow, baseflow, overland flow, channel routing, and evaporation from the interception, soil water, ground water, snowpack and channel water. The modelled basin could be split up into 20 segments, each with different hydrological characteristics, but a maximum of 3 segments was used in this study.
A total of 29 parameters was used in the model although only 10 (other than those directly measurable) were found to be sensitive in the three basins. The simulated streamflow did not match the real flows very well due to errors in the data input and due to simplifications in the model. It was found, however, that as the proportion of the limestone in a segment increased the overland flow decreased, the interflow increased, the baseflow and interflow recessions were faster, the soil storages were smaller and the infiltration rate was higher, than in segments with a larger proportion of exposed clastics. The flow characteristics of the inaccessible conduits were inferred from the channel routing parameters and it was postulated that the majority of the underground flow in the karst basins was taking place under vadose conditions.
This is the first detailed examination of the karst geomorphology of the Bruce Peninsula. It attempts to review all aspects including pavement phenomena and formation (microkarst features), surface and subsurface karst hydrology (meso to macro scale) and water chemistry. The latter is based on over 250 samples collected in 1973 and 1974.
The dolomite pavement is the best example of its kind that has been described in the literature. It covers much of the northern and eastern parts of the peninsula and can be differentiated into three types based on karren assemblages. Two of these are a product of lithology and the third reflects local environmental controls. The Amabel Formation produces characteristic karren such as rundkarren, hohlkarren, meanderkarren, clint and grike, kamentizas and rillenkarren on glacially abraded biohermal structures. The Guelph Formation develops into a very irregular, often cavernous surface with clint and grike and pitkarren as the only common recognizable karren. The third assemblage is characterized by pitkarren and is found only in the Lake Huron littoral zone. Biological factors are believed to have played a major role in the formation of the pavement. Vegetation supplies humic acids which help boost the solution process and helps to maintain a wet surface. This tends to prolong solution and permit the development of karren with rounded lips and bottoms.
Three types of drainage other than normal surface runoff are found on the Bruce. These are partial underground capture of surface streams, complete underground capture (fluvio-karst), and wholly vertical drainage without stream action (holokarst). Holokarst covers most of the northern and eastern edge of the peninsula along the top of the escarpment. Inland it is replaced by fluvial drainage, some of which has been, or is in the process of being captured. Four perennial streams and one lake disappear into sinkholes. These range from very simple channel capture and resurgence, as shown by a creek east of Wiarton, to more mature and complex cave development of the St. Edmunds cave near Tobermory. Partial underground capture represents the first stage of karst drainage. This was found to occur in one major river well inland of the fluvio-karst and probably occurs in other streams as well. This chapter also examines the possible future karst development of the Bruce and other karst feature such as isolated sinks and sea caves.
The water chemistry presented in Chapter 5 represents the most complete data set from southern Ontario. It is examined on a seasonal basis as well as grouped into classes representing water types (streams, Lake Huron and Georgian Bay, inland lakes, swamps, diffuse springs and conduit springs). The spring analyses are also fitted into climatic models of limestone solution based on data from other regions of North America. It was found that solution rates in southern Ontario are very substantial. Total hardness ranges from 150 to 250 ppm (expressed as CaCO3) in most lakes and streams and up to 326 ppm in springs. These rates compare with more southerly latitudes. The theoretical equilibrium partial pressure of CO2 was found to be the most significant chemical variable for comparing solution on different kinds of carbonates and between glaciated and non-glaciated regions. Expect for diffuse flow springs and Lake Huron, the Bruce data do not separate easily into water types using either graphical or statistical (i.e. Linear Discriminant Analysis) analyses. This is partly because of the seasonality of the data and because of the intimate contact all waters have with bedrock.
![]() |
![]() |