Karst and Cave RSS news feed Like us on Facebook! follow us on Twitter!
Community news

Speleology in Kazakhstan

Shakalov on 04 Jul, 2018
Hello everyone!   I pleased to invite you to the official site of Central Asian Karstic-Speleological commission ("Kaspeko")   There, we regularly publish reports about our expeditions, articles and reports on speleotopics, lecture course for instructors, photos etc. ...

Speleology in Kazakhstan

Shakalov on 04 Jul, 2018
Hello everyone!   I pleased to invite you to the official site of Central Asian Karstic-Speleological commission ("Kaspeko")   There, we regularly publish reports about our expeditions, articles and reports on speleotopics, lecture course for instructors, photos etc. ...

Speleology in Kazakhstan

Shakalov on 11 Jul, 2012
Hello everyone!   I pleased to invite you to the official site of Central Asian Karstic-Speleological commission ("Kaspeko")   There, we regularly publish reports about our expeditions, articles and reports on speleotopics, lecture course for instructors, photos etc. ...

New publications on hypogene speleogenesis

Klimchouk on 26 Mar, 2012
Dear Colleagues, This is to draw your attention to several recent publications added to KarstBase, relevant to hypogenic karst/speleogenesis: Corrosion of limestone tablets in sulfidic ground-water: measurements and speleogenetic implications Galdenzi,

The deepest terrestrial animal

Klimchouk on 23 Feb, 2012
A recent publication of Spanish researchers describes the biology of Krubera Cave, including the deepest terrestrial animal ever found: Jordana, Rafael; Baquero, Enrique; Reboleira, Sofía and Sendra, Alberto. ...

Caves - landscapes without light

akop on 05 Feb, 2012
Exhibition dedicated to caves is taking place in the Vienna Natural History Museum   The exhibition at the Natural History Museum presents the surprising variety of caves and cave formations such as stalactites and various crystals. ...

Did you know?

That water catchment is the intake of water from an aquifer or a surface reservoir [16].?

Checkout all 2699 terms in the KarstBase Glossary of Karst and Cave Terms

What is Karstbase?



Browse Speleogenesis Issues:

KarstBase a bibliography database in karst and cave science.

Featured articles from Cave & Karst Science Journals
Chemistry and Karst, White, William B.
Engineering challenges in Karst, Stevanović, Zoran; Milanović, Petar
See all featured articles
Featured articles from other Geoscience Journals
Geochemical and mineralogical fingerprints to distinguish the exploited ferruginous mineralisations of Grotta della Monaca (Calabria, Italy), Dimuccio, L.A.; Rodrigues, N.; Larocca, F.; Pratas, J.; Amado, A.M.; Batista de Carvalho, L.A.
Karst environment, Culver D.C.
Mushroom Speleothems: Stromatolites That Formed in the Absence of Phototrophs, Bontognali, Tomaso R.R.; D’Angeli Ilenia M.; Tisato, Nicola; Vasconcelos, Crisogono; Bernasconi, Stefano M.; Gonzales, Esteban R. G.; De Waele, Jo
Calculating flux to predict future cave radon concentrations, Rowberry, Matt; Marti, Xavi; Frontera, Carlos; Van De Wiel, Marco; Briestensky, Milos
See all featured articles from other geoscience journals

Search in KarstBase

Your search for numerical model (Keyword) returned 96 results for the whole karstbase:
Showing 1 to 15 of 96
The role of tributary mixing in chemical variations at a karst spring, Milandre, Switzerland, , Perrin J. , Jeannin P. Y. , Cornaton F. ,
SummarySolute concentration variations during flood events were investigated in a karst aquifer of the Swiss Jura. Observations were made at the spring, and at the three main subterraneous tributaries feeding the spring. A simple transient flow and transport numerical model was able to reproduce chemographs and hydrographs observed at the spring, as a result of a mixing of the concentration and discharge of the respective tributaries. Sensitivity analysis carried out with the model showed that it is possible to produce chemical variations at the spring even if all tributaries have constant (but different for each of them) solute concentrations. This process is called tributary mixing. The good match between observed and modelled curves indicate that, in the phreatic zone, tributary mixing is probably an important process that shapes spring chemographs. Chemical reactions and other mixing components (e.g. from low permeability volumes) have a limited influence.Dissolution-related (calcium, bicarbonate, specific conductance) and pollution-related parameters (nitrate, chloride, potassium) displayed slightly different behaviours: during moderate flood events, the former showed limited variations compared to the latter. During large flood events, both presented chemographs with significant changes. No significant event water participates in moderate flood events and tributary mixing will be the major process shaping chemographs. Variations are greater for parameters with higher spatial variability (e.g. pollution-related). Whereas for large flood events, the contribution of event water becomes significant and influences the chemographs of all the parameters. As a result, spring water vulnerability to an accidental pollution is low during moderate flood events and under base flow conditions. It strongly increases during large flood events, because event water contributes to the spring discharge

Assessing the importance of conduit geometry and physical parameters in karst systems using the storm water management model (SWMM), , Peterson Eric W. , Wicks Carol M. ,
SummaryQuestions about the importance of conduit geometry and about the values of hydraulic parameters in controlling ground-water flow and solute transport through karstic aquifers have remained largely speculative. One goal of this project was to assess the role that the conduit geometry and the hydraulic parameters have on controlling transport dynamics within karstic aquifers. The storm water management model (SWMM) was applied to the Devil's Icebox-Connor's Cave System in central Missouri, USA. Simulations with incremental changes to conduit geometry or hydraulic parameters were performed with the output compared to a calibrated baseline model. Ten percent changes in the length or width of a conduit produced statistically significant different fluid flow responses. The model exhibited minimal sensitivity to slope and infiltration rates; however, slight changes in Manning's roughness coefficient can highly alter the simulated output.Traditionally, the difference in flow dynamics between karstified aquifers and porous media aquifers has led to the idea that modeling of karst aquifers is more difficult and less precise than modeling of porous media aquifers. When evaluated against models for porous media aquifers, SWMM produced results that were as accurate (10% error compared to basecase). In addition, SWMM has the advantage of providing data about local flow. While SWMM may be an appropriate modeling technique for some karstic aquifers, SWMM should not be viewed as a universal solution to modeling karst systems

Anomalous behaviour of specific electrical conductivity at a karst spring induced by variable catchment boundaries: the case of the Podstenjšek spring, Slovenia, , Ravbar, N. , Engelhardt, I. , Goldscheider, N.

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.

Hydrogeology of the Umm Er Radhuma aquifer, Saudi Arabia, with reference to fossil gradients, 1982, Bakiewicz W, Milne Dm, Noori M,
Much of North Africa and the Arabian peninsula, lying in the Saharan climate zone, are underlain by huge tabular sandstone and carbonate aquifers, ranging in age from Cambrian to Tertiary. These are often saturated with water of reasonable quality and form very valuable resources in an area often desperately short of water. The Palaeocene Umm Er Radhuma carbonate aquifer is one such formation which has been the subject of intensive recent investigation. The formation contains groundwater of a reasonable quality, has adequate transmission and storage characteristics and hence considerable potential for future development. The origin of the water in such aquifers is the subject of continuing controversy. It is not disputed that the water is moving under the influence of regional groundwater gradients but origins of these gradients are the subject of considerable argument. On the one hand, there are those who hold that the presently observed gradients are fossil remnants of conditions created by a much wetter climatic regime prevalent some thousands of years ago. Against this are those who maintain that the gradients, at least in part, reflect a present day system with groundwater discharge in approximate dynamic equilibrium with recharge. This paper examines the hydrogeology of a typical Middle Eastern formation of the disputed kind, the Umm Er Radhuma aquifer in Saudi Arabia, and, with the aid of analytical and numerical models, attempts to resolve the problem of the origin of the observed groundwater gradients and to discover the extent to which the past must influence present day plans for future development

The Hydrology of a Glacierised Alpine Karst Castlegaurd Mountain, Alberta, PhD Thesis, 1983, Smart, Charles Christopher

Alpine karst throughout the world has been affected by past glaciation, and yet little is known of the interactions between glacier ice and karst. This dissertation attempts to gain some understanding of the problem through the study of the Castleguard Area, Alberta, where a karst aquifer is presently overlain by temperate glacier ice.
Quantitative fluorometric tracing and hydrometric measurements generated a broad data base on aquifer behaviour. Tracer breakthrough curves were interpreted using a new systematic approach which considers an explicit set of processes likely to affect the particular tracer under the given experimental conditions. Non-linearity in aquifer behaviour and rapid groundwater velocities demonstrated the aquifer to be an extreme conduit type Conduit springs are elements in a vertical hierarchy in which the topmost springs are "overflows" and exhibit greater flow variability than their associated "underflows". A numerical model was developed to simulate a conduit aquifer. It demonstrated that pulse train and recession analysis widely accepted methods of karst aquifer investigation, could be rather misleading when applied to conduit aquifers.
Interactions between ice and groundwater were observed at two scales: regulation water appeared to feed a diffuse percolation system and supraglacial melt passed into subglacial conduits which entered open vadose shafts. Karst is unlikely to be entirely subglacial in origin because of the limited aggressiveness of subglacial waters.
The Castlegaurd karst appeared to have originated preglacially in response to the breaching of impermeable caprock. Glaciation re-ordered the landscape and produced abundant clastic debris which subsequently blocked or obstructed karst conduits. Much of the resulting karst is paragenetic and comparatively immature due to glacial disruption and slow growth rates. Geomorphic and hydrologic interactions between ice and karst depend intimately upon the relationship between the geographic zones of the glacier and the aquifer.

Thermal aspects of the East Midlands aquifer system, 1987, Wilson N. P. , Luheshi M. N. ,
A case study of a heat flow anomaly in the E Midlands of England is reported. The anomaly has been suggested to be an effect of water movement at depth within the E Midlands basin, with recharge to the Lower Carboniferous limestones in their outcrop, eastward movement and ascent of water up a steep faulted anticline at Eakring where the heat flow measurements were made. Numerical modelling of heat and fluid flow has been undertaken for a section running from the Peak District through Eakring to the coast. The results indicate that, although an anomaly is expected for reasonable values of hydrological parameters, its magnitude is less than that observed. The geological structure at Eakring is such that three-dimensional flow is likely to be important, and this could easily account for the discrepancy between the modelling results and the observations. The regional water flow regime has other effects on heat flow, notably the depression of heat flow above the Sherwood Sandstone aquifer

Principles of early development of karst conduits under natural and man-made conditions revealed by mathematical analysis of numerical models., 1996, Dreybrodt W.

Principles of early development of karst conduits under natural and man-made conditions revealed by mathematical analysis of numerical models, 1996, Dreybrodt W,
Numerical models of the enlargement of primary fissures in limestone by calcite aggressive water show a complex behavior. If the lengths of the fractures are large and hydraulic heads are low, as is the case in nature, dissolution rates at the exit of the channel determine its development by causing a slow increase of water flow, which after a long gestation time by positive feedback accelerates dramatically within a short time span. Mathematical analysis of simplified approximations yields an analytical expression for the breakthrough time, when this happens, in excellent agreement with the results of a numerical model. This expression quantifies the geometrical, hydraulic, and chemical parameters determining such karat processes. If the lengths of the enlarging channels are small, but hydraulic heads are high, as is the case for artificial hydraulic structures such as darns, it is the widening at the entrance of the flow path which determines the enlargement of the conduit. Within the lifetime of the dam this can cause serious water losses, This can also be explained by mathematical analysis of simplified approximations which yield an analytical threshold condition from which the safety of a dam can be judged. Thus in both cases the dynamic processes of karstification are revealed to gain a deeper understanding of the early development of karst systems. As a further important result, one finds that minimum conditions, below which karstification cannot develop, do not exist

The applicability of numerical models to adequately characterize ground-water flow in karstic and other triple-porosity aquifers, 1996, Quinlan J. F. , Davies G. J. , Jones S. W. , Huntoon P. W.

Structure et comportement hydraulique des aquifers karstiques, DSc thesis, 1996, Jeannin, P. Y.

This thesis aims to provide a better knowledge of karst flow systems, from a functional point of view (behaviour with time), as well as from a structural one (behaviour in space). The first part of the thesis deals with the hydrodynamic behaviour of karst systems, and the second part with the geometry of karstic networks, which is a strong conditioning factor for the hydrodynamic behaviour.
Many models have been developed in the past for describing the hydrodynamic behaviour of karst hydrogeological systems. They usually aim to provide a tool to extrapolate, in time and/or space, some characteristics of the flow fields, which can only be measured at a few points. Such models often provide a new understanding of the systems, beyond what can be observed directly in the field. Only special field measurements can verify such hypotheses based on numerical models. This is an significant part of this work. For this purpose, two experimental sites have been equipped and measured: Bure site or Milandrine, Ajoie, Switzerland, and Holloch site, Muotathal, Schwyz, Switzerland. These sites gave us this opportunity of simultaneously observe hydrodynamic parameters within the conduit network and, in drillholes, the "low permeability volumes" (LPV) surrounding the conduits.
These observations clearly show the existence of a flow circulation across the low permeability volumes. This flow may represent about 50% of the infiltrated water in the Bure test-field. The epikarst appears to play an important role into the allotment of the infiltrated waters: Part of the infiltrated water is stored at the bottom of the epikarst and slowly flows through the low permeability volumes (LPV) contributing to base flow. When infiltration is significant enough the other part of the water exceeds the storage capacity and flows quickly into the conduit network (quick flow).
For the phreatic zone, observations and models show that the following scheme is adequate to describe the flow behaviour: a network of high permeability conduits, of tow volume, leading to the spring, is surrounded by a large volume of low permeability fissured rock (LPV), which is hydraulically connected to the conduits. Due to the strong difference in hydraulic conductivity between conduits and LPV, hydraulic heads and their variations in time and space are strongly heterogeneous. This makes the use of piezometric maps in karst very questionable.
Flow in LPV can be considered as similar to flow in fractured rocks (laminar flow within joints and joints intersections). At a catchment scale, they can be effectively considered as an equivalent porous media with a hydraulic conductivity of about 10-6 to 10-7 m/s.
Flow in conduits is turbulent and loss of head has to be calculated with appropriate formulas, if wanting any quantitative results. Our observations permitted us to determine the turbulent hydraulic conductivity of some simple karst conduits (k', turbulent flow), which ranges from 0.2 to 11 m/s. Examples also show that the structure of the conduit network plays a significant role on the spatial distribution of hydraulic heads. Particularity hydraulic transmissivity of the aquifer varies with respect to hydrological conditions, because of the presence of overflow conduits located within the epiphreatic zone. This makes the relation between head and discharge not quadratic as would be expected from a (too) simple model (with only one single conduit). The model applied to the downstream part of Holloch is a good illustration of this phenomena.
The flow velocity strongly varies along the length of karst conduits, as shown by tracer experiments. Also, changes in the conduit cross-section produce changes in the (tow velocity profile. Such heterogeneous flow-field plays a significant role in the shape of the breakthrough curves of tracer experiments. It is empirically demonstrated that conduit enlargements induce retardation of the breakthrough curve. If there are several enlargements one after the other, an increase of the apparent dispersivity will result, although no diffusion with the rock matrix or immobile water is present. This produces a scale effect (increase of the apparent dispersivity with observation scale). Such observations can easily be simulated by deterministic and/or black box models.
The structure of karst conduit networks, especially within the phreatic zone, plays an important role not only on the spatial distribution of the hydraulic heads in the conduits themselves, but in the LPV as well. Study of the network geometry is therefore useful for assessing the shape of the flow systems. We further suggest that any hydrogeological study aiming to assess the major characteristics of a flow system should start with a preliminary estimation of the conduit network geometry. Theories and examples presented show that the geometry of karst conduits mainly depends on boundary conditions and the permeability field at the initial stage of the karst genesis. The most significant boundary conditions are: the geometry of the impervious boundaries, infiltration and exfiltration conditions (spring). The initial permeability field is mainly determined by discontinuities (fractures and bedding planes). Today's knowledge allows us to approximate the geometry of a karst network by studying these parameters (impervious boundaries, infiltration, exfiltration, discontinuity field). Analogs and recently developed numerical models help to qualitatively evaluate the sensitivity of the geometry to these parameters. Within the near future, new numerical tools will be developed and will help more closely to address this difficult problem. This development will only be possible if speleological networks can be sufficiently explored and used to calibrate models. Images provided by speleologists to date are and will for a long time be the only data which can adequately portray the conduit networks in karst systems. This is helpful to hydrogeologists. The reason that we present the example of the Lake Thun karst system is that it illustrates the geometry of such conduits networks. Unfortunately, these networks are three-dimensional and their visualisation on paper (2 dimensions) is very restrictive, when compared to more effective 3-D views we can create with computers. As an alternative to deterministic models of speleogenesis, fractal and/or random walk models could be employed.

Mixed transport reaction control of gypsum dissolution kinetics in aqueous solutions and initiation of gypsum karst, 1997, Raines M. A. , Dewers T. A. ,
Experiments with gypsum in aqueous solutions at 25 degrees C, low ionic strengths, and a range of saturation states indicate a mixed surface reaction and diffusional transport control of gypsum dissolution kinetics. Dissolution rates were determined in a mixed flow/rotating disc reactor operating under steady-state conditions, in which polished gypsum discs were rotated at constant speed and reactant solutions were continuously fed into the reactor. Rates increase with velocity of spin under laminar conditions (low rates of spin), but increase asymptotically to a constant rate as turbulent conditions develop with increasing spin velocity, experiencing a small jump in magnitude across the laminar-turbulent transition. A Linear dependence of rates on the square root of spin velocity in the laminar regime is consistent with rates being limited by transport through a hydrodynamic boundary layer. The increase in rate with onset of turbulence accompanies a near discontinuous drop in hydrodynamic boundary layer thickness across the transition. A relative independence of rates on spinning velocity in the turbulent regime plus a nonlinear dependence of rates on saturation state are factors consistent with surface reaction control. Together these behaviors implicate a 'mixed' transport and reaction control of gypsum dissolution kinetics. A rate law which combines both kinetic mechanisms and can reproduce experimental results under laminar flow conditions is proposed as follows: R = k(t) {1 - Omega(b)() zeta [1 - (1 2(1 - Omega(b)())(1/2)]} where k(t) is the rate coefficient for transport control, and Omega(b)() is the mean ionic saturation state of the bulk fluid. The dimensionless parameter zeta(=Dm(eq)()/2 delta k() where m(eq)() = mean ionic molal equilibrium concentration, D is the diffusion coefficient through the hydrodynamic boundary layer, delta equals the boundary layer thickness and k() is the rate constant for surface reaction control) indicates which process, transport or surface reaction, dominates, and is sensitive to the hydrodynamic conditions in the reactor. For the range of conditions used in our experiments, zeta varies from about 1.4 to 4.5. Rates of gypsum dissolution were also determined in situ in a cavern system in the Permian Blaine Formation, southwestern Oklahoma. Although the flow conditions in the caverns were not determinable, there is good agreement between lab- and field-determined rates in that field rate magnitudes lie within a range of rates determined experimentally under zero to low spin velocities A numerical model coupling fluid flow and gypsum reaction in an idealized circular conduit is used to estimate the distance which undersaturated solutions will travel into small incipient conduits before saturation is achieved. Simulations of conduit wall dissolution showed-member behaviors of conduit formation and surface denudation that depend on flow boundary conditions (constant discharge or constant hydraulic gradient and initial conduit radius. Surface-control of dissolution rates. which becomes more influential with higher fluid flow velocity, has the effect that rate decrease more slowly as saturation is approached than otherwise would occur if rates were controlled by transport alone. This has the effect that reactive solutions can penetrate much farther into gypsum-bearing karst conduits than heretofore thought possible, influencing timing and mechanism of karst development as well as stability of engineered structures above karst terrain

Numerical simulation as a tool for checking the interpretation of karst spring hydrographs, 1997, Eisenlohr L, Kiraly L, Bouzelboudjen M, Rossier Y,
A schematic representation of karst aquifers may be that of a high hydraulic conductivity channel network with kilometre-wide intervals, surrounded by a low hydraulic conductivity fractured limestone volume and connected to a local discharge area, the karst spring, The behaviour of the karst spring (hydrographs, chemical or isotopic composition, etc.) represents the global response of the karst aquifer to input events. The available data an karst aquifer hydraulic parameters are limited, Global response is therefore more easily obtained and is commonly used to make inferences on the recharge and groundwater How processes, as well as on the hydraulic parameter fields. Direct verification of these interpretations is, obviously, very difficult. We have used an indirect method of verification, consisting of introducing well-defined theoretical karst structures into a finite element model and then analysing the simulated global response according to presently accepted interpretation schemes. As we know what we put into the numerical model, the validity of any interpretation may be checked. The first results indicate that some of the generally accepted interpretations are not necessarily true. In particular: (i) separation of simulated recession hydrographs into several components shows that different exponential components do not necessarily correspond to aquifer volumes with different hydraulic conductivities: (ii) non-exponential parts of recession hydrographs do not always give information about the infiltration process: and (iii) the recession coefficient of the baseflow (i.e. the last, nearly exponential part of the recession hydrograph) depends on the global configuration of the whole karst aquifer, not just on the hydraulic properties of the low hydraulic conductivity volumes. (C) 1997 Elsevier Science B.V

Simulation of the evolution of maze caves., 1997, Clemens T. , Hiickinghaus D. , Sauter M. , Liedl R. , Teutsch G.
The development of cave systems in carbonate rocks depends on a variety of boundary and initial conditions. Among the cave systems, two main types of geometries can be distinguished: the dendritic and maze pattern. A numerical model has been developed capable of modeling the genesis of karst systems in complex geological environments. It is applied to simulate the development of the above mentioned two different types of cave geometries. The results confirm that a prerequisite for the development of maze caves is evenly distributed recharge (White 1969). However more important for the development of maze caves is that flow through the system is restricted by an overlying less conductive horizon, e.g. a sandstone caprock. Thus the feed back mechanism of higher flow rates leading to higher dissolution rates and therefore the preferential development of a small number of tubes does not dominate the evolution of the karst aquifer. This hydraulic restriction furthers the development of other conduits also to achieve a significant diameter.

Parameter identification in double-continuum models applied to karst aquifers., 1997, Mohrlok U. , Kienie J. , Teutsch G.
One modelling approach which proved successful in describing the groundwater flow within karst terraines is based on the double-continuum concept. This concept was first introduced by Teutsch (1988) and subsequently used by Teutsch & Sauter (1991), Sauter (1992), Lang (1995), Mohrlok (1996) and others to describe the ambivalent characteristics of karst aquifers. However, the approach has the drawback that the double-continuum model parameters can be determined only through model calibration (inverse approach), i.e. so far the model parameters cannot be related directly to physical field measurements. Therefore, in order to develop a better understanding of the physical significance of hydraulic parameters within double-continuum systems, a detailed numerical modeling study was conducted. For this purpose a number of synthetic but realistic karst aquifer network geometries were generated and analysed. The response of the karst network to recharge events was simulated using a detailed discrete fracture flow model with the resulting head and spring flow variographs being subsequently assumed as field measurements. This 'measured' data was then used for the calibration of a double-continuum model and the resulting parameters were compared to the original karst network geometry data. This comparison was used to develop mathematical/physical relationships between the discrete karst network geometry representing reality and the double-continuum parameter representation of it.

Physical response of a karst drainage basin to flood pulses: Example of the Devil's Icebox cave system (Missouri, USA), 1998, Halihan T. , Wicks C. M. , Engeln J. F. ,
In karst aquifers, water moves from the recharge area (sinkhole plains and swallets) to the discharge area (springs), traveling kilometers through the groundwater system in a period of hems to days. Transit rimes through karst aquifers are a function of the conduit geometry and connectedness, intensity and duration of the recharge event, and antecedent soil moisture. Often many of these factors are unknown or difficult to quantify. Therefore, predicting the response of a karst basin to recharge is difficult. Numerous researchers have attempted to understand the response of a karst basin, but a good understanding of whether the response is dependent on local features or regional effects is currently lacking. From April 1994 to May 1995, flood pulse hydrographs from a karst aquifer with well-developed and well-documented conduits (Devil's Icebox cave system) were obtained from a gaging station near the spring of the karst basin. Data were also collected from within the conduit system in an attempt to determine whether flow was locally controlled by constrictions in the conduits. Based on an application of Bernoulli's equation, analyses of the changes in kinetic head and potential head over time indicated local control during storm events. The observed sediment patterns and water level variations also support localized flow control during storm events. A numerical model of the constrictions was rested and reproduced the responses observed at the spring during initial periods of storm events. The model illustrated that the constricted flow was very sensitive to recharge. It also illustrated the transition from local control due to constriction to regional controls due to the aquifer matrix. (C) 1998 Elsevier Science B.V

Results 1 to 15 of 96
You probably didn't submit anything to search for