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. ...

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 albedo is the ratio of reflected radiation to total radiation on a natural surface [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.
See all featured articles
Featured articles from other Geoscience Journals
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
Microbial mediation of complex subterranean mineral structures, Tirato, Nicola; Torriano, Stefano F.F;, Monteux, Sylvain; Sauro, Francesco; De Waele, Jo; Lavagna, Maria Luisa; D’Angeli, Ilenia Maria; Chailloux, Daniel; Renda, Michel; Eglinton, Timothy I.; Bontognali, Tomaso Renzo Rezio
Evidence of a plate-wide tectonic pressure pulse provided by extensometric monitoring in the Balkan Mountains (Bulgaria), Briestensky, Milos; Rowberry, Matt; Stemberk, Josef; Stefanov, Petar; Vozar, Jozef; Sebela, Stanka; Petro, Lubomir; Bella, Pavel; Gaal, Ludovit; Ormukov, Cholponbek;
See all featured articles from other geoscience journals

Search in KarstBase

Your search for flow simulation (Keyword) returned 9 results for the whole karstbase:
Paleohydrology and Streamflow Simulation of three Karst Basins in Southeastern West Virginia, U.S.A., PhD Thesis, 1975, Coward, Julian Michael Henry

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.

Rospo Mare field is located in the Adriatic Sea, 20 km of the Abruzzes coast, at an average depth of 80 m. The reservoir is a karst which is essentially conductive; yet unlike a conventional porous medium, it cannot be simulated by the usual tools and techniques of reservoir simulation. Therefore, several approaches were used to describe the flow mechanism during the production period in greater detail. The first approach consisted of generating three-dimensional images which were constrained by both petrophysical and geological factors and then, using up-scaling techniques, obtaining the equivalent permeabilities (scalar or tensorial) of grid blocks located in different zones within the karst. This approach shows that within the infiltration zone it is possible, whatever the scale, to find an equivalent homogeneous porous medium; on the other hand, within the epikarst this equivalent medium does not exist below pluridecametric dimensions. Thus it is impossible to study the sweeping mechanism on a small scale, so we must use a deterministic model which describes the network of pipes in the compact matrix, in which a waterflood is simulated by means of a conform finite-element model. This constituted the second approach. The third and final approach consisted of inventing a system of equations to analytically solve the pressure field in a network of vertical pipes which are intersected by a production drain and submitted to a strong bottom water-drive. This model allows us to simulate the water-oil contact rise within the reservoir and study the flows depending on the constraints applied to the production well. It appears that cross flows occur in the pipes even during the production period

Determining karst transmissivities with inverse modeling and an equivalent porous media, 1999, Larocque M. , Banton O. , Ackerer P. , Razack M. ,
Flow simulation is difficult to implement in heterogeneous media such as karst aquifers, primarily because the structure of the rock is extremely complex and usually unknown. The aim of this study was to verify the possibility of using inverse modeling and an equivalent porous media to identify transmissivities in a slightly karstified aquifer, the La Rochefoucauld karst (Charente, France), Different simulation scenarios were tested: using two spatial discretizations with different finite-element cell sizes and using measured or interpolated heads. The inverse modeling was performed with the downscaling parameterization procedure, using a finite-element representation of bidimensional ground water flow. The inverse modeling converged satisfactorily with all scenarios: head residuals were small and spring flow rates and the river/aquifer exchanges were adequately stimulated. The scenario using small cells and measured heads generated a highly heterogeneous transmissivity field, indicating an overparameterization of the problem. The calibrated transmissivities and simulated heads of this scenario proved less reliable overall than those of the other scenarios. The use of interpolated heads generated more uniform transmissivities as a result of the head smoothing. A rotation of the initial parameter mesh showed that the scenarios using interpolated heads generate the most stable and reliable results. The scenarios with interpolated heads could therefore be used when head measurements are limited or are unevenly distributed over the aquifer. Overall, the calibrated transmissivities reproduced the entire range of transmissivities measured in the field using different methods. The results indicate that inverse modeling and an equivalent porous media can be used to determine transmissivities in a moderately karstified aquifer

Rainfall-runoff relations for karstic springs: multifractal analyses, 2002, Labat D. , Mangin A. , Ababou R. ,
Karstic watersheds appear as highly as non-linear and non-stationary systems. The behaviour of karstic springs has been previously studied using non-linear simulation methods (Volterra expansion) and non-stationary analyses methods based on wavelet transforms. The main issue of karstic spring behaviour consists of the presence and the identification of characteristic time-scales. In order to highlight more precisely the scale-properties of the rainfall-runoff relations for karstic springs, the multifractal analysis is introduced. These methods are applied daily and half-hourly rainfall rates and runoffs measured on a three French karstic springs located in the Pyrenees Mountains (Ariege, France): Aliou, Baget and Fontestorbes. They are characterised by a variable development of the drainage systems. We have at our disposal long and uninterrupted series of data over period of several years, which constitute a high quality bank data. Multifractal analyses of both daily and half-hourly rainfall rates and runoffs give evident a scale-dependant behaviour. Effectively, it highlights the presence of different multifractal processes at each sampling rate. Using a universal class of multifractal models based on cascade multiplicative processes, the identified multifractal sub-processes are characterised by the classical parameters alpha and C-1. All these results should lead to several improvements in karstic springflow simulation models. (C) 2002 Elsevier Science B.V. All rights reserved

On the importance of geological heterogeneity for flow simulation, 2006, Eaton Tt,
Geological heterogeneity is recognized as a major control on reservoir production and constraint on many aspects of quantitative hydrogeology. Hydrogeologists and reservoir geologists need to characterize groundwater flow through many different types of geological media for different purposes. In this introductory paper, an updated perspective is provided on the current status of the long effort to understand the effect of geological heterogeneity on flow using numerical simulations. A summary is given of continuum vs. discrete paradigms, and zonal vs. geostatistical approaches, all of which are used to structure model domains. Using these methods and modern simulation tools, flow modelers now have greater opportunities to account for the increasingly detailed understanding of heterogeneous aquifer and reservoir systems.One way of doing this would be to apply a broader interpretation of the idea of hydrofacies, long used by hydrogeologists. Simulating flow through heterogeneous geologic media requires that numerical models capture important aspects of the structure of the flow domain. Hydrofacies are reinterpreted here as scale-dependent hydrogeologic units with a particular representative elementary volume (REV) or structure of a specific size and shape. As such, they can be delineated in indurated sedimentary or even fractured aquifer systems, independently of lithofacies, as well as in the unlithified settings in which they have traditionally been used. This reconsideration of what constitutes hydrofacies, the building blocks for representing geological heterogeneity in flow models, may be of some use in the types of settings described in this special issue

Evidence for hydraulic heterogeneity and anisotropy in the mostly carbonate Prairie du Chien Group, southeastern Minnesota, USA, 2006, Tipping Robert G. , Runkel Anthony C. , Alexander Jr. E. Calvin, Alexander Scott C. , Green Jeffery A. ,
In southeastern Minnesota, Paleozoic bedrock aquifers have typically been represented in groundwater flow simulations as isotropic, porous media. To obtain a more accurate hydrogeologic characterization of the Ordovician Prairie du Chien Group, a new approach was tested, combining detailed geologic observations, particularly of secondary porosity, with hydraulic data. Lithologic observations of the depositional and erosional history of the carbonate-dominated bedrock unit constrained characterization of both primary (matrix) and secondary porosity from outcrops and core. Hydrostratigraphic data include outcrop and core observations along with core plug permeability tests. Hydrogeologic data include discrete interval aquifer tests, borehole geophysics, water chemistry and isotope data, and dye trace studies. Results indicate that the Prairie du Chien Group can be subdivided into the Shakopee aquifer at the top, consisting of interbedded dolostone, sandstone and shale, and the underlying Oneota confining unit consisting of thickly bedded dolostone. The boundary between these two hydrogeologic units does not correspond to lithostratigraphic boundaries, as commonly presumed. Groundwater flow in the Shakopee aquifer is primarily through secondary porosity features, most commonly solution-enlarged bedding planes and sub-horizontal and vertical fractures. Regional scale preferential development of cavernous porosity and permeability along specific stratigraphic intervals that correspond to paleokarst were also identified, along with a general depiction of the distribution of vertical and horizontal fractures. The combination of outcrop and core investigations, along with borehole geophysics, discrete interval aquifer tests, water chemistry and isotope data and dye trace studies show that the Prairie du Chien Group is best represented hydrogeologically as heterogeneous and anisotropic. Furthermore, heterogeneity and anisotropy within the Prairie du Chien Group is mappable at a regional scale (> 15,000 km2)


There are several numerical modelling approaches of speleogenesis in existence today. They take into account physical and chemical laws for flow and dissolution in fractured carbonate aquifers. Nevertheless, the initial void networks considered by these models generally do not correspond to the fracturing reality. The approach proposed here aims to simulate speleogenesis in an aquifer characterized by a fracture network, while matching field reality as closely as possible and respecting geometrical properties. Using statistical and geometrical parameters obtained by field observations and analogue experiments, it is possible to generate 3-D realistic networks in terms of the relative position of joints that control the overall network connectivity. Once the fracture networks are generated, they are adapted and incorporated in a 3-D ground water flow and transport finite element model. The flow simulations in the fracture networks allow determination of the spatial distribution of flow velocities for the initial configuration. This distribution, added to other information such as age and travel time, is used to simulate the evolution of the apertures of the different elements. This paper mainly presents the theoretical basis for the proposed method, from the fracturing model to the incorporation of the generated network in the flow model. Then, it describes the principles leading to forthcoming first benchmark simulations which will be used to develop the analogical rules concerning karst aquifer evolution, and for lead sensibility analysis.

Evaluation of permeability and non-Darcy flow in vuggy macroporous limestone aquifer samples with lattice Boltzmann methods, 2013, Sukop M. C. , Huang H. , Alvarez P. F. , Variano E. A. , Cunningham K. J.

Lattice Boltzmann flow simulations provide a physics-based means of estimating intrinsic permeability from pore structure and accounting for inertial flow that leads to departures from Darcy’s law. Simulations were used to compute intrinsic permeability where standard measurement methods may fail and to provide better understanding of departures from Darcy’s law under field conditions. Simulations also investigated resolution issues. Computed tomography (CT) images were acquired at 0.8 mm interscan spacing for seven samples characterized by centimeter-scale biogenic vuggy macroporosity from the extremely transmissive sole-source carbonate karst Biscayne aquifer in southeastern Florida. Samples were as large as 0.3 m in length; 7–9 cm-scale-length subsamples were used for lattice Boltzmann computations. Macroporosity of the subsamples was as high as 81%. Matrix porosity was ignored in the simulations. Non-Darcy behavior led to a twofold reduction in apparent hydraulic conductivity as an applied hydraulic gradient increased to levels observed at regional scale within the Biscayne aquifer; larger reductions are expected under higher gradients near wells and canals. Thus, inertial flows and departures from Darcy’s law may occur under field conditions. Changes in apparent hydraulic conductivity with changes in head gradient computed with the lattice Boltzmann model closely fit the Darcy-Forchheimer equation allowing estimation of the Forchheimer parameter. CT-scan resolution appeared adequate to capture intrinsic permeability; however, departures from Darcy behavior were less detectable as resolution coarsened.

Nonlinear Flow Process: A New Package to Compute Nonlinear Flow in MODFLOW, 2014, Mayaud C. , Walker P. , Hergarten S. , Birk S.

A new MODFLOW package (Nonlinear Flow Process; NLFP) simulating nonlinear flow following the Forchheimer equation was developed and implemented in MODLFOW-2005. The method is based on an iterative modification of the conductance calculated and used by MODFLOW to obtain an effective Forchheimer conductance. The package is compatible with the different layer types, boundary conditions, and solvers as well as the wetting capability of MODFLOW. The correct implementation is demonstrated using four different benchmark scenarios for which analytical solutions are available. A scenario considering transient flow in a more realistic setting and a larger model domain with a higher number of cells demonstrates that NLFP performs well under more complex conditions, although it converges moderately slower than the standard MODFLOW depending on the nonlinearity of flow. Thus, this new tool opens a field of opportunities to groundwater flow simulation with MODFLOW, especially for core sample simulation or vuggy karstified aquifers as well as for nonlinear flow in the vicinity of pumping wells.

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