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 peat is decomposed matter, mainly vegetable [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 16 to 30 of 96
Modeling of storm responses in conduit flow aquifers with reservoirs, 1998, Halihan Todd, Wicks Carol M. ,
In aquifers containing large voids, such as karst aquifers with caves or basaltic aquifers with lava tubes, hydrographs at wells or springs are used to analyze the structure and response of the hydrogeological system. Numerical modeling of hydrograph response is commonly based on either inverse techniques or postulated flow geometries. However, the range of mechanisms for generating hydrograph responses have not been fully investigated.Physical modeling of these complex non-Darcian systems permits better understanding of the storm responses that conduit systems may generate. Using a numerical model of conduit flow systems which incorporates turbulent flow, some of the mechanisms that can alter storm pulses were investigated by treating them as combinations of pipes that connect reservoirs.The results indicate that the response of a conduit-flow aquifer can range from what has been called 'diffuse' or 'steady' to 'conduit' or 'flashy', without employing a diffusive component. A full range of behavior can be the result of changes from phreatic to epiphreatic conditions in a conduit, changes in conduit geometry, or multiple springs draining the same system. The results provide a quantitative tool to assess spring and well hydrographs, and illustrate mechanisms that can generate observed responses, which have previously been qualitatively interpreted

Structure et comportement hydraulique des aquifers karstiques, DSc. Thesis, faculte des Sciences de l'Universite de Neuchatel., 1998, Jeannin Py.
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.

On predicting contaminant transport in carbonate terrains: Behavior and prediction, 1999, Annable W. K , Sudicky E. A.
A three-dimensional numerical model was used to quantify the fate of conservative transport in carbonate terrainsNumerical flow and transport experiments were conducted in proto-conduit scale limestone terrains (conduits less than 10cm) which determined that a priori information on the 'spill' and/or 'tracer injection location' and discharge locations provided little insight in characterizing the complexity of the internal labyrinth of interconnected conduitsScaling, aside from the characterization of the geologic media, was one of the most limiting factors in quantifying recharge tracers or contaminant distributionHowever, if sufficient numbers of discharge locations (springs) are known, the extent of downstream contaminant migration can be characterized

Solutionally enhanced leakage rates of dams in karst regions, 1999, Bauer S. , Birk S. , Liedl R. , Sauter M.
This paper presents numerical model studies regarding the development of leakage rates of dams in karst regions due to solutional enlargement of conduits in the rock beneath the reservoirA hybrid continuum-discrete flow model (CAVE) is used for the modelingThe fractured carbonate rock beneath the reservoir is represented by a network of initial conduitsThe effects of a preferential flowpath and a grout curtain on leakage rates are studied in different scenariosFor the parameters considered in this paper, the simulation results suggest that dissolutional widening of the network leads in all scenarios to leakage rates endangering the reservoir within 100 years

The initiation of hypogene caves in fractured limestone by rising thermal water: investigation of a parallel series of competing fractures, 1999, Dumont K. A. , Rajaram H. , Budd D. A.
Integrated cave systems can either form at or near the surface of the earth (epigenic) or at some depth below the earth's surface (hypogenic)For caves that form in fractured limestone, the two most common types of cave-system morphologies are branchwork and mazeworkBranchwork caves are composed of tributaries that coalesce in the downstream direction, similar to surface streamsMazework caves exhibit two or more sets of parallel passages intersecting in a grid-like patternThe majority of epigenic caves exhibit branchwork morphologies, which represent the dominance of individual flow pathsIn contrast, mazework caves develop when dissolution occurs along numerous flow pathsWhereas most epigenic caves are related to surficial meteoric flow systems, some mazework caves are thought to have formed in hypogene environments where rising thermal water cools in response to the geothermal gradientOur objective is to examine the fundamental cause for the difference in morphology between epigenic and thermal hypogenic cave systems using numerical modelsIn particular, we are examining the competition between different flow paths in fractured limestone undergoing dissolutional enlargementAs noted in previous numerical studies, epigenic systems are characterized by the dominance of a single flow path, which is consistent with the structure of epigenic cavesSo, in order to explain the structure of maze caves, one has to explain why no single flow path attains dominanceThe retrograde solubility of calcite coupled with heat transfer from the fluid to the rock is hypothesized to provide the mechanism by which dissolutional power is distributed among all competing flow pathsNumerical models of fluid flow, heat transfer, and calcite dissolution chemistry are integrated to develop a model of hypogene cave initiation in fractured limestoneFlow is assumed to occur in the presence of a spatially variable rock temperature field that is constant through timePreliminary numerical modeling results for a system of parallel fractures demonstrate the differences in the nature of competition between flow paths in epigenic (constant temperature) and hypogenic systems (flow in the presence of a negative thermal gradient)Differences in results using various kinetic models for calcite dissolution are also presentedThe role of aperture variation and distribution in a parallel set of fractures is also examined

Structural effects on carbonate aquifers, 1999, Sasowsky I. D
Structural geology affects the behavior of karst aquifers by controlling the overall placement and orientation of the limestone and through fracturesThe placement and orientation affect the position of recharge and discharge boundaries to the system, while the fractures serve as pathways for water movementWhen creating a conceptual or numerical model of a karst site, it is useful and cost-efficient to consider all of these effects, as well as the geologic and geomorphic history of the areaBy understanding structural controls on the genesis of the aquifer, predictions can be made regarding current-day behavior in terms of heterogeneity and anisotropy of flowBecause conduits and fissures mainly form along structurally created discontinuities, structural data can be very useful for understanding aquifer behavior, and determining specific high-conductivity flowpaths

Infiltration measured by the drip of stalactites, 2000, Sanz E. , Lopez J. J. ,
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 Beke 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

Influence of fracture roughness on karstification times, 2000, Dreybrodt W. , Gabrovsek F.
Karst aquifers develop along water-transmitting fractures which are rough by nature. Most numerical models on the evolution of the widths of those fractures by limestone dissolution, however, approximate such rough fractures by two parallel smooth planes. Here we investigate the influence of fracture roughness on the breakthrough time of conduits. All results show that even when the flow through a fracture is reduced by its roughness by a factor of ten in comparison to a corresponding smooth fracture, breakthrough time is increased only by a factor of 4. This gives a first answer to the influence of natural fracture roughness on the evaluation of karst and shows that extreme roughness, unlikely in nature, is necessary to exert significant influence on karstification time.

Genesis of a large cave system: the case study of the North of Lake Thun system (Canton Bern, Switzerland), 2000, Jeannin Py. , Bitterli T. , Hauselmann P.
The genesis of the cave system in the region Hohgant-Sieben Hengste-Lake of Thun (more than 250 km of surveyed passage) has been reconstructed based on speleomorphological observations (mainly by observing where the morphology changes from vadose to phreatic). Eight flow systems (phases) and their respective conduit networks have been distinguished so far. The oldest had a phreatic level at an altitude of 1950 m a.s.l. The last corresponds to today's phreatic zone located at 658 m a.s.l. Between each system, the water table dropped several hundred meters. This appears to be a consequence of changes in boundary conditions, mainly the springis position, which moved down as a tectonic uplift and deepening of the nearby valleys occured. Observations demonstrate that phreatic conduits are sometimes developed close to the ancient water table, but often much deeper, down to 200 to 400 m below this level. The change from one phase to the next seems to have been quick. This stepwise evolution is compatible with the results of computer models which give durations of 10'000 to 30i000 years for conduits systems to develop. Analysis of the conduit networks of each flow system shows that their geometry is mainly influenced by the hydraulic gradients and the overall geometry of the aquifer. The orientation of discontinuity surfaces (fractures and bedding planes) and/or their intersections, play a subordinate role. This is also supported by numerical models found in the literature. As, despite a high fracture density, we observe deep rather than shallow phreatic loops, we assume that the heterogeneity of the discontinuity openings plays a more important role in the depth of karstification than the frequency of the discontinuities.

On flow and transport of miscible tracers in the vadose zone, 2000, Loboferreira J. P. , Leitã, O Teresa, Quigley Seá, N P. , Theves Thomas

The paper presents a synthesis of LNEC's contribution to the European Commission (DGXII) Contract n° CI1*-CT94-0014 (DG 12 HSMU), based on the results reported by Lobo-Ferreira et al. (1998). The main contribution of LNEC to the project was the development of flow and transport experiments, for different flow type conditions (i.e. saturated and unsaturated). These included experiments at: (1) different scales (two scales in the laboratory and one in a medium scale - artificial aquifer), (2) different saturation conditions, and (3) different tracers (MgCl2, CaCl2, NO3-, Ni, Cd). The goals of the experiments were threefold: (1) analysis and quantification of the physical parameters that control circulation in this zone (hydraulic conductivity correspondent to the degree of saturation etc.); (2) determination and quantification of the more important processes that control the chemical behaviour of heavy metals and nitrate, and (3) to obtain data for calibrating the numerical models.

HJWFTAC: software for Hantush-Jacob analysis of variable-rate, multiple-extraction well pumping tests, 2002, Fleming Sw, Ruskauff Gj, Adams A,
Analytical well test solutions are a powerful approach to aquifer characterization and the parameterization of comprehensive numerical models. In addition, wellfield drawdown tests, which consist of coordinated pumping and data collection at a suite of monitoring and operating production wells, are of growing significance due to increasing pressures upon groundwater resources and the consequent management and planning requirement for superior hydrogeologic characterization of existing production wellfields. However. few pumping test analysis codes accommodate the multiple extraction wells involved, particularly for more sophisticated analytic aquifer test solutions. We present and demonstrate here a FORTRAN code for analysis of drawdown at a monitor well due to simultaneous variable-rate pumping at multiple independent production wells, which we developed in response to a need to refine an existing numerical, coupled groundwater/surface water resource management model, Spatial and temporal superposition are used to accommodate the typical operational properties of wellfield pumping tests, The software invokes the well-accepted Hantush-Jacob method for semiconfined or 'leaky' aquifers in a forward simulation procedure and effectively assumes homogeneity in applicable aquifer parameters (transmissivity, coefficient of storage, and leakance). Intended for both professionals and students, the code is widely applicable and straightforward to use as written. However, it can be modified with relative ease to use alternative well test solutions and/or formal inverse modeling techniques, or to accommodate spatial hydrogeologic variability. An application to a pumping test conducted in a karst limestone aquifer at the Cross Bar Ranch wellfield in Tampa Bay, Florida, demonstrates the utility of the software. (C) 2002 Elsevier Science Ltd. All rights reserved

Is the water still hot? Sustainability and the thermal springs at Bath, England, 2002, Atkinson Tc, Davison Rm,
The hot springs at Bath are the largest natural thermal source in Britain. Sustainable use of the waters for a spa requires maintenance of their temperature and flow rate. Together with smaller springs at Hotwells, Bristol, they form the outflow from a regional thermal aquifer that occurs where the Carboniferous Limestone is buried at depths > 2.7 km in the Bristol-Bath structural basin. The aquifer is recharged via limestone outcrops forming the south and west portions of the basin rim. Current knowledge of the basin's structure is reviewed, and important uncertainties identified concerning the hydrogeological role of thrust faults which may cut the limestone at depth. A simple numerical model is used to determine the possible influence of thrusts upon groundwater flow within the thermal aquifer. Comparison of the modelled flow patterns with geochemical data and structure contours eliminates the hypothesis that thrusts completely disrupt the continuity of the aquifer. The most successful model is used to simulate the possible impact of dewatering by large quarries at the limestone outcrops north and south of Bath. Substantial reductions in modelled flow at Bath result from proposed dewatering in the eastern Mendips, although the steady-state approach adopted has severe limitations in that it does not take account of the incremental staging of actual dewatering, nor allow for partial restitution of groundwater levels. The geological uncertainties highlighted by the modelling could be addressed by future research into the effect of thrusts on the continuity of the Carboniferous Limestone. More refined modelling to predict the timing of possible impacts of quarry dewatering will require measurements of the storativity of the thermal aquifer

Assessments of the sensitivity to climate change of flow and natural water quality in four major carbonate aquifers of Europe, 2002, Younger P. L. , Teutsch G. , Custodio E. , Elliot T. , Manzano M. , Sauter M. ,
A numerical modelling approach has been developed to predict the vulnerability of aquifers to future climate change. This approach encompasses changes in recharge regime, dynamics of flow and storage patterns within aquifers, and natural hydrochemical changes. An application of the approach has been made to four hypothetical spring catchments representative of major carbonate aquifers in three European climatic zones. Since prolific carbonate aquifers typically combine a high transmissivity with a low specific yield, they can be expected to be more sensitive than clastic aquifers to changes in recharge patterns. Simulations of the study systems to the middle of the 21st century predict different outcomes in the three different climate zones: (1) in the northern maritime zone (UK) recharge (and therefore discharge) is predicted to increase by as much as 21 0n response to anticipated increases in precipitation; (2) in the continental zone (Germany) recharge in winter is predicted to remain approximately the same as at present, but summer recharge will decline dramatically (by as much as 32%), so that a net decrease in aquifer discharge is predicted; and (3) in the Mediterranean zone (Spain) recharge is predicted to decrease by as much as 160f the present-day values. For all three systems, increases in water hardness in response to rising CO2 are predicted, but are expected to be negligible in water resources terms

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico, 2002, Gonzalezherrera R. , Pinto I. , Gamboavargas J. ,
The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7-10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity

Analytical 1D dual-porosity equivalent solutions to 3D discrete single-continuum models. Application to karstic spring hydrograph modelling, 2002, Cornaton F, Perrochet P,
One-dimensional analytical porosity-weighted solutions of the dual-porosity model are derived, providing insights on how to relate exchange and storage coefficients to the volumetric density of the high-permeability medium. It is shown that porosity-weighted storage and exchange coefficients are needed when handling highly heterogeneous systems-such as karstic aquifers-using equivalent dual-porosity models. The sensitivity of these coefficients is illustrated by means of numerical experiments with theoretical karst systems. The presented ID dual-porosity analytical model is used to reproduce the hydraulic responses of reference 3D karst aquifers, modelled by a discrete single-continuum approach. Under various stress conditions, simulation results show the relations between the dual-porosity model coefficients and the structural features of the discrete single-continuum model. The calibration of the equivalent 1D analytical dual-porosity model on reference hydraulic responses confirms the dependence of the exchange coefficient with the karstic network density. The use of the analytical model could also point out some fundamental structural properties of the karstic network that rule the shape of the hydraulic responses, such as density and connectivity. (C) 2002 Elsevier Science B.V. All rights reserved

Results 16 to 30 of 96
You probably didn't submit anything to search for