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 surface water is water obtained from surface supplies [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 porous-media (Keyword) returned 29 results for the whole karstbase:
Showing 1 to 15 of 29
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

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

The kinetics of the reaction CO2?>H? as one of the rate limiting steps for the dissolution of calcite in the system H2O-CO2-CaCO3, 1996, Dreybrodt W, Lauckner J, Liu Zh, Svensson U, Buhmann D,
Dissolution of CaCO3 in the system H2O-CO2-CaCO3 is controlled by three rate-determining processes: The kinetics of dissolution at the mineral surface, mass transport by diffusion, and the slow kinetics of the reaction H2O CO2 = H HCO3-. A theoretical model of Buhmann and Dreybrodt (1985a,b) predicts that the dissolution rates depend critically on the ratio V/A of the volume V of the solution and the surface area A of the reacting mineral. Experimental data verifying these predictions for stagnant solutions have been already obtained in the range 0.01 cm < V/A < 0.1 cm. We have performed measurements of dissolution rates in a porous medium of sized CaCO3 particles for V/A in the range of 2 . 10(-4) cm and 0.01 cm in a system closed with respect to CO2 using solutions pre-equilibrated with an initial partial pressure of CO2 of 1 . 10(-2) and 5 . 10(-2) atm. The results are in satisfactory agreement with the theoretical predictions and show that especially for V/A < 10(-3) cm dissolution is controlled entirely by conversion of CO2 into H and HCO3-, whereas in the range from 10(-3) cm up to 10(-1) cm both CO2-conversion and molecular diffusion are the rate controlling processes. This is corroborated by performing dissolution experiments using 0.6 mu molar solutions of carbonic anhydrase, an enzyme enhancing the CO2-conversion rates by several orders of magnitude. In these experiments CO2 conversion is no longer rate limiting and consequently the dissolution rates of CaCO3 increase significantly. We have also performed batch experiments at various initial pressures of CO2 by stirring sized calcite particles in a solution with V/A = 0.6 cm and V/A = 0.038 cm. These data also clearly show the influence of CO2-conversion on the dissolution rates. In all experiments inhibition of dissolution occurs close to equilibrium. Therefore, the theoretical predictions are valid for concentrations c less than or equal to 0.9 c(eq). Summarising we find good agreement between experimental and theoretically predicted dissolution rates. Therefore, the theoretical model can be used with confidence to find reliable dissolution rates from the chemical composition of a solution for a wide field of geological applications

Precipitation kinetics of calcite in the system CaCO3-H2O-CO2: The conversion to CO2 by the slow process H?->CO2? as a rate limiting step, 1997, Dreybrodt W, Eisenlohr L, Madry B, Ringer S,
Precipitation rates of CaCO3 from supersaturated solutions in the H2O - CO2 - CaCO3 system are controlled by three rate-determining processes: the kinetics of precipitation at the mineral surface, mass transport of the reaction species involved to and from the mineral surface, and the slow kinetics of the overall reaction HCO3- H --> CO2 H2O. A theoretical model by Buhmann and Dreybrodt (1985a,b) taking these processes into account predicts that, due to the slow kinetics of this reaction, precipitation rates to the surface of CaCO3 minerals depend critically on the ratio V/A of the volume V of the solution to the surface area A of the mineral in contact with it, for both laminar and turbulent flow. We have performed measurements of precipitation rates in a porous medium of sized particles of marble, limestone, and synthetic calcite, with V/A ratios ranging from 3.10(-4) to 1.2-10(-2) cm at 10 degrees C. Calcite was precipitated from supersaturated solutions with [Ca2] approximate to 4 mmol/L and an initial P-CO2 of 5.10(-3) or 1.10(-3) atm, respectively, using experimental conditions which prevented exchange of CO2 with the atmosphere, i.e., closed system. The results are in qualitative agreement with the theoretical predictions. Agreement with the observed data, however, is obtained by modifying the rate law of Plummer et al. (1978) to take into account surface-controlled inhibition effects. Experiments with supersaturated solutions containing carbonic anhydrase, an enzyme which enhances the conversion of HCO3- into CO2, yield rates increased by a factor of up to 15. This provides for the first time unambiguous experimental evidence that this reaction is rate limiting. We have also measured precipitation rates in batch experiments, stirring sized mineral particles in a solution with V/A ranging from 0.03 to 0.75 cm. These experiments also give clear evidence on the importance of the conversion of HCO3- into CO2 as rate limiting step. Taken together our experiments provide evidence that the theoretical model of Buhmann and Dreybrodt (1985a,b) can be used to predict reliable rates from the composition of CaHCO3- solutions with low ionic strength in many geologically relevant situations. Copyright (C) 1997 Elsevier Science Ltd

The problem of modeling limestone springs: The case of Bagnara (north Apennines, Italy), 1997, Angelini P, Dragoni W,
The Bagnara spring (Central Italy), fed by a fractured, carbonate, and, in some areas, karstic aquifer, was examined. The available information is derived from geological mapping and daily flows over a period of 20 consecutive years. There are no data on the hydrogeological parameters nor on the aquifer hydraulic head, which is known only at the elevation of the spring. The objective of the work was to construct an appropriate mathematical model for the spring despite the scarcity of available information. The MODFLOW code was used to simulate the system following the equivalent porous media approach. The hydraulic conductivity and the specific yield equivalents were estimated by calibrating the model on the master depletion curve and taking into consideration the topographic elevation of the system's surface. The size of the protection area around the spring was investigated on the basis of the isochrons constructed from the results of the model

Particle size distributions in waters from a karstic aquifer: from particles to colloids, 1997, Atteia O, Kozel R,
Waters from the surface hydrologic network and the spring of a karstic aquifer in Switzerland were sampled to analyse their colloidal content. The measurements were done weekly with a single particle counter and were verified by other techniques. The particle size distribution (PSD) was modelled in two portions, below and above 5 mu m, using two types of equation: a power law (Pareto distribution) and an exponential law. The model results matched well with the entire PSD data set by varying the parameter values. The parameters obtained from fitting the measured PSD curves were then interpreted in relation to environmental factors. It appears that the two parts of the curves vary independently. The first part of the PSD curve, relating to the smallest particles, is dependent on the pH value of the spring or the temperature of the surface brook. In contrast, the second part of the curve depends mostly on the spring discharge volume. During high flow events, the major effect of the discharge on particle size occurs during the rising limb of the hydrograph, interpreted as clays deposited in the aquifer and resuspended due to high water velocity. The contrasted behaviour of the two parts of the PSD curves suggested that the break point in the curves represents the limit between colloidal and particulate behaviour. Knowing these dependencies, and the characteristics of the particulate matter, allowed the estimation of the role of the colloids in contaminant transport. Large fluxes of suspended matter, specific to karstic aquifers, demonstrate the critical role of colloids in contaminant transport, which is markedly different from what typically occurs in porous media.

Evolution of size distributions of natural particles during aggregation: modelling versus field results, 1998, Atteia O,
In this paper a discretized model simulating aggregation of size distributions jointly with sedimentation and transport is presented. A review of the current theory provides some helpful hints about the relative importance of each aggregation process, i.e. Brownian motion, shear flow and differential sedimentation, which are tested by using collision efficiency factors. The novel aspect of the model arises from the use of a varying mean particle diameter in each size class. This allows both non-steady-state and steady-state calculations and free choice of size classes. A comparison with a classical approach shows the exactitude of the results and the improvment obtained for several cases. The simulations gave a family of curves characterized by three parts corresponding to peri-, and orthokinetic aggregation and to sedimentation. The role of collision effciency is crucial in the relative extent of each part of the size distribution. The comparison with a series of data from a karst spring showed that the model was able to fit most of the particle size distributions using significant values of each parameter. This allowed information about particle aggregation and transport within a non-accessible aquifer to be inferred.

Influence of aperture variability on dissolutional growth of fissures in karst formations, 1998, Hanna R. B. , Rajaram H. ,
The influence of aperture variability on dissolutional growth of fissures is investigated on the basis of two-dimensional numerical simulations. The logarithm of the fissure aperture before dissolution begins is modeled as a Gaussian stationary isotropic random field. The initial phase of dissolutional growth is studied up to the time when turbulent flow first occurs at a point within the fissure (the breakthrough time). The breakthrough time in variable aperture fissures is smaller than that in uniform fissures and decreases as the coefficient of variation of the aperture field (sigma/mu) increases. In comparing uniform and variable aperture fissures in limestone, the breakthrough time with sigma/mu = 0.1 is about a factor of 2 smaller than that in a uniform fissure. The breakthrough time is reduced by about an order of magnitude with sigma/mu = 2.0. The mechanism leading to reduced breakthrough times is the focusing of flow into preferential flow channels which are enlarged at a faster rate than the surrounding regions of slower flow. Dissolution channels are narrower and more tortuous as sigma/mu. increases. Investigations of the influence of reaction rate reveal that the influence of aperture variability is more pronounced in rapidly dissolving rock. In uniform fissures in rapidly dissolving minerals, breakthrough times are very long since water becomes saturated with respect to the mineral within a short distance of the entrance to the flow path. However, in variable aperture fissures, breakthrough occurs rapidly because of selective growth along preferential flow channels, which progressively capture larger fractions of the total flow. These results partly explain why conduits develop rapidly in gypsum, although previous one-dimensional studies suggest that conduit growth will not occur

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

Transient-state history matching of a karst aquifer ground water flow model, 2000, Larocque M. , Banton O. , Razack M. ,
Ground water flow modeling in a karst aquifer presents many difficulties. In particular, the hydrodynamic properties and the now behavior can vary over time. History matching of transient-state conditions is required to test the accuracy of the model under varying hydrodynamic conditions. The objective of this study was to illustrate how transient-state conditions can be used to history match a ground water flow model of a large aquifer, the La Rochefoucauld karst (Charente, France). The model used a porous medium equivalent and was based on a steady-state calibration of hydraulic conductivities. The history match consisted of studying the simulated heads and spring flow rates to test the capacity of the model to reproduce different aspects of the aquifer behavior, The simulated heads and flow rates were analyzed as new data using correlation and spectral analyses to compare the temporal structures of the measured and simulated time series. The analyses provided information on the storage capacity of the aquifer, the input-output delays, the degree of correlation between input and output, and the length of the impulse response of the aquifer, These data were used to study the impact of the hypotheses underlying the model (hydraulic conductivities, storage coefficient, representation of rivers, use of a porous medium equivalent). The results show that the model adequately simulates the overall behavior of the studied aquifer, The model can be used under variable hydrodynamic conditions to simulate ground water flow on a regional scale. This case study illustrates how a complete history match of a simplified representation of reality can lead to an adequate mathematical tool

A strategy for modeling ground water rebound in abandoned deep mine systems, 2001, Adams R, Younger Pl,
Discharges of polluted water from abandoned mines are a major cause of degradation of water resources worldwide, Pollution arises after abandoned workings flood up to surface level, by the process termed ground water rebound, As flow in large, open mine voids is often turbulent, standard techniques for modeling ground water flow (which assume laminar flow) are inappropriate for predicting ground water rebound. More physically realistic models are therefore desirable, yet these are often expensive to apply to all but the smallest of systems. An overall strategy for ground water rebound modeling is proposed, with models of decreasing complexity applied as the temporal and spatial scales of the systems under analysis increase. For relatively modest systems (area < 200 km(2)), a physically based modeling approach has been developed, in which 3-D pipe networks (representing major mine roadways, etc.) are routed through a variably saturated, 3-D porous medium (representing the country rock). For systems extending more than 100 to 3000 km(2), a semidistributed model (GRAM) has been developed, which conceptualizes extensively interconnected volumes of workings as ponds, which are connected to other ponds only at discrete overflow points, such as major inter-mine roadways, through which flow can be efficiently modeled using the Prandtl-Nikuradse pipe-flow formulation. At the very largest scales, simple water-balance calculations are Probably as useful as any other approach, and a variety of proprietary codes may be used for the purpose

Dispersion, retardation and scale effect in tracer breakthrough curves in karst conduits, 2001, Hauns M. , Jeannin P. Y. , Atteia O. ,
Characteristics of tracer breakthrough curves in karst conduits are examined and compared to results generated using well known equations applied to porous media. The equations of the turbulent dispersion lead to a transport equation similar to the classical advection-dispersion equation for porous media with a slightly different meaning for the dispersion and advection terms. For investigations at the meter length scale, we used a three-dimensional (3-D) computational fluid dynamics (CFD) code to simulate tracer transport in several conduit geometries. The simulations show that turbulent dispersion can be considered as Fickian at a meter length scale of observation and that turbulent dispersivity depends linearly on the average flow velocity in the range of observed velocities. The simulations show that pools induce retardation (tailing of the breakthrough curve) due to flow reversal in eddies. Retardation has a complex relationship with the pool dimensions. Irregularity of the conduit cross-section along the investigated section clearly produces retardation. This is obvious at the meter length scale but may still be visible 10(3) m downstream from the injection point. A transfer function ('black box') approach is used for upscaling from a meter to a 10(3) m length scale. Before applying it to natural examples, the transfer function approach is tested by using the 3-D CFD code and appears to perform well. Several tests, based on numerical, laboratory and held experiments, of conduit segments which includes various dispersive features indicate that retardation tends to be transformed to symmetrical dispersion with distance. At large scale it appears that the dominant dispersion factor is the irregularity of the conduit geometry, which produces an increase in dispersivity with distance ('scale effect'), similar to that observed in porous media. In conclusion this suggests that retardation and high dispersion provide evidence of an irregular conduit, including either numerous dispersive features or large-scale ones (pools for example). Conversely no retardation and moderate dispersion (close to 0.012 m) must result from turbulent Row through a smooth conduit. (C) 2001 Elsevier Science B.V. All rights reserved

Coastal karst springs in the Mediterranean basin : study of the mechanisms of saline pollution at the Almyros spring (Crete), observations and modelling, 2002, Arfib B, De Marsily G, Ganoulis J,
Variations in salinity and flow rate in the aerial, naturally salty spring of Almyros of Heraklion on Crete were monitored during two hydrological cycles. We describe the functioning of the coastal karstic system of the Almyros and show the influence of the duality of the flow in the karst (conduits and fractured matrix) on the quality of the water resource in the coastal area. A mechanism of saltwater intrusion into this highly heterogeneous system is proposed and validated with a hydraulic mathematical model, which describes the observations remarkably well. Introduction. - Fresh groundwater is a precious resource in many coastal regions, for drinking water supply, either to complement surface water resources, or when such resources are polluted or unavailable in the dry season. But coastal groundwater is fragile, and its exploitation must be made with care to prevent saltwater intrusion as a result of withdrawal, for any aquifer type, porous, fractured or karstic. In karstic zones, the problem is very complex because of the heterogeneous nature of the karst, which makes it difficult to use the concept of representative elementary volume developed for porous or densely fractured systems. The karstic conduits focus the major part of the flow in preferential paths, where the water velocity is high. In coastal systems, these conduits have also an effect on the distribution of the saline intrusion. As was shown e.g. by Moore et al. [1992] and Howard and Mullings [1996], both freshwater and salt-water flow along the fractures and conduits to reach the mixing zone, or the zone where these fluids are superposed in a dynamic equilibrium because of their differences in density ; but the dynamics of such a saltwater intrusion are generally unknown and not represented in models. Such coastal karstic systems are intensely studied at this moment in the Mediterranean region [Gilli, 1999], both as above sea-level or underwater springs, for potential use in areas where this resource would be of great value for economic development. This article discusses the freshwater-saltwater exchange mechanisms in the karstic aquifer of the Almyros of Heraklion aquifer (Crete) and explains the salinity variations observed in the spring. First, the general hydrogeology of the study site is described, then the functioning of the spring : a main conduit drains the freshwater over several kilometres and passes at depth through a zone where seawater is naturally present. The matrix-conduit exchanges are the result of pressure differences between the two media. These processes are represented in a mathematical model that confirms their relevance. General hydrogeology of the studied site. - The karstic coastal system of the Almyros of Heraklion (Crete) covers 300 km2 in the Ida massif whose borders are a main detachment fault, and the Sea of Crete in the north, the Psiloritis massif (highest summit at 2,456 m) in the south and west, and the collapsed basin of Heraklion filled in by mainly neo-geneous marl sediments in the east. The watershed basin consists of the two lower units of characteristic overthrust formations of Crete (fig. 1) : the Cretaceous Plattenkalk and the Cretaceous Tripolitza limestones. The two limestone formations are locally separated by interbedded flysch or phyllade units that form an impervious layer [Bonneau et al., 1977 ; Fassoulas, 1999] and may lead to different flow behaviour within the two karstic formations. Neo-tectonic activity has dissected these formations with large faults and fractures. The present-day climate in Crete is of Mediterranean mountain type, with heavy rain storms and snow on the summits in winter. Rainfall is unevenly distributed over the year, with 80 % of the annual total between October and March and a year-to-year average of 1,370 mm. The flow rate of the spring is high during the whole hydrologic cycle, with a minimum in summer on the order of 3 m3.s-1 and peak flow in winter reaching up to 40 m3.s -1. The water is brackish during low flow, up to a chloride content of 6 g.l-1, i.e. 23 % of seawater, but it is fresh during floods, when the flow rate exceeds 15 m3.s-1. During the 1999-2000 and 2000-2001 hydrologic cycles, the water was fresh during 14 and 31 days, respectively. The water temperature is high and varies very little during the year (see table I). In the areas of Keri and Tilissos (fig. 1), immediately south of the spring, the city of Heraklion extracts water from the karstic system through a series of 15 wells with depth reaching 50 to 100 m below sea level. Initially, when the wells were drilled, the water was fresh, but nowadays the salinity rises progressively, but unequally from well to well (fig. 2). The relatively constant temperatures and salinities of the wells, during the hydrological cycle, contrast with the large salinity variations at the spring (fig. 2 and table I). They show that the karstic system is complex and comprises different compartments, where each aquifer unit reacts to its individual pressures (pumping, rainfall) according to its own hydrodynamic characteristics [Arfib et al., 2000]. The Almyros spring seems disconnected from the surrounding aquifer and behaves differently from that which feeds the wells (upper Tripolitza limestone). It is recharged by fresh water from the mountains, which descends to depths where it probably acquires its salinity. The spring would thus be the largest resource of the area, if it was possible to prevent its pollution by seawater. A general functioning sketch is proposed (fig. 3), which includes the different geological units of interest. Identification of the functioning of the Almyros spring through monitoring of physical and chemical parameters. - The functioning of the aquifer system of the Almyros spring was analysed by monitoring, over two hydrological cycles, the level of the spring, the discharge, the electric conductivity and the temperature recorded at a 30 min time interval. In the centre of the watershed basin, a meteorological station at an altitude of 800 m measures and records at a 30 min time interval the air temperature, rainfall, relative humidity, wind velocity and direction ; moreover, an automatic rain gauge is installed in the northern part of the basin at an altitude of 500 m. The winter floods follow the rhythm of the rainfall with strong flow-rate variations. In contrast, the summer and autumn are long periods of drought (fig. 7). The flow rate increases a few hours after each rainfall event ; the water salinity decreases in inverse proportion to the flow rate a few hours to a few days later. Observations showed that the water volume discharged at the Almyros spring between the beginning of the flow rate increase and the beginning of the salinity decrease is quite constant, around 770,000 m3 (fig. 4) for any value of the flow rate, of the salinity and also of the initial or final rainfall rates. To determine this constant volume was of the upmost importance when analyzing the functioning of the Almyros spring. The lag illustrates the differences between the pressure wave that moves almost instantaneously through the karst conduit and causes an immediate flow rate increase after rainfall and the movement of the water molecules (transfer of matter) that arrives with a time lag proportionate to the length of the travel distance. The variation of the salinity with the flow rate acts as a tracer and gives a direct indication of the distance between the outlet and the seawater entrance point into the conduit. In the case of the Almyros, the constant volume of expelled water indicates that sea-water intrusion occurs in a portion of the conduit situated several kilometres away from the spring (table II), probably inland, with no subsequent sideways exchange in the part of the gallery leading up to the spring. As the lag between the flow rate and the salinity recorded at the spring is constant, one can correct the salinity value by taking, at each time step, with a given flow rate, the salinity value measured after the expulsion of 770,000 m3 at the spring, which transforms the output of the system so as to put the pressure waves and the matter transfer in phase [Arfib, 2001]. After this correction, the saline flux at the spring, equal to the flow rate multiplied by the corrected salinity, indicates the amount of sea-water in the total flow. This flux varies in inverse proportion to the total flow rate in the high-flow period and the beginning of the low-flow period, thereby demonstrating that the salinity decrease in the spring is not simply a dilution effect (fig. 5). The relationship that exists between flow rate and corrected salinity provides the additional information needed to build the conceptual model of the functioning of the part of the Almyros of Heraklion aquifer that communicates with the spring. Freshwater from the Psiloritis mountains feeds the Almyros spring. It circulates through a main karst conduit that descends deep into the aquifer and crosses a zone naturally invaded by seawater several kilometers from the spring. The seawater enters the conduit and the resulting brackish water is then transported to the spring without any further change in salinity. The conduit-matrix and matrix-conduit exchanges are governed by the head differences in the two media. Mathematical modelling of seawater intrusion into a karst conduit Method. - The functioning pattern exposed above shows that such a system cannot be treated as an equivalent porous medium and highlights the influence of heterogeneous structures such as karst conduits on the quantity and quality of water resources. Our model is called SWIKAC (Salt Water Intrusion in Karst Conduits), written in Matlab(R). It is a 1 D mixing-cell type model with an explicit finite-difference calculation. This numerical method has already been used to simulate flow and transport in porous [e.g. Bajracharya and Barry, 1994 ; Van Ommen, 1985] and karst media [e.g. Bauer et al., 1999 ; Liedl and Sauter, 1998 ; Tezcan, 1998]. It reduces the aquifer to a single circular conduit surrounded by a matrix equivalent to a homogeneous porous medium where pressure and salinity conditions are in relation with sea-water. The conduit is fed by freshwater at its upstream end and seawater penetrates through its walls over the length L (fig. 6) at a rate given by an equation based on the Dupuit-Forchheimer solution and the method of images. The model calculates, in each mesh of the conduit and at each time step, the head in conditions of turbulent flow with the Darcy-Weisbach equation. The head loss coefficient {lambda} is calculated by Louis' formula for turbulent flow of non-parallel liquid streams [Jeannin, 2001 ; Jeannin and Marechal, 1995]. The fitting of the model is intended to simulate the chloride concentration at the spring for a given matrix permeability (K), depth (P) and conduit diameter (D) while varying its length (L) and its relative roughness (kr). The spring flow rates are the measured ones ; at present, the model is not meant to predict the flow rate of the spring but only to explain its salinity variations. Results and discussion. - The simulations of chloride concentrations were made in the period from September 1999 to May 2001. The depth of the horizontal conduit where matrix-conduit exchanges occur was tested down to 800 m below sea level. The diameter of the conduit varied between 10 and 20 m, which is larger than that observed by divers close to the spring but plausible for the seawater intrusion zone. The average hydraulic conductivity of the equivalent continuous matrix was estimated at 10-4 m/s. A higher value (10-3 m/s) was tested and found to be possible since the fractured limestone in the intrusion zone may locally be more permeable but a smaller value (10-5 m/s) produces an unrealistic length (L) of the saline intrusion zone (over 15 km). For each combination of hydraulic conductivity, diameter and depth there is one set of L (length) and kr (relative roughness) calibration parameters. All combinations for a depth of 400 m or more produce practically equivalent results, close to the measured values. When the depth of the conduit is less than 400 m, the simulated salinity is always too high. Figure 7 shows results for a depth of 500 m, a diameter of 15 m and a hydraulic conductivity of 10-4 m/s. The length of the saltwater intrusion zone is then 1,320 m, 4,350 m away from the spring and the relative roughness coefficient is 1.1. All the simulations (table II) need a very high relative roughness coefficient which may be interpreted as an equivalent coefficient that takes into account the heavy head losses by friction and the variations of the conduit dimensions which, locally, cause great head losses. The model simulates very well the general shape of the salinity curve and the succession of high water levels in the Almyros spring but two periods are poorly described due to the simplicity of the model. They are (1) the period following strong freshwater floods, where the model does not account for the expulsion of freshwater outside the conduit and the return of this freshwater which dilutes the tail of the flood and (2) the end of the low-water period when the measured flux of chlorides falls unexpectedly (fig. 5), which might be explained by density stratification phenomena of freshwater-saltwater in the conduit (as observed in the karst gallery of Port-Miou near Cassis, France [Potie and Ricour, 1974]), an aspect that the model does not take into account. Conclusions. - The good results produced by the model confirm the proposed functioning pattern of the spring. The regulation of the saline intrusion occurs over a limited area at depth, through the action of the pressure differences between the fractured limestone continuous matrix with its natural saline intrusion and a karst conduit carrying water that is first fresh then brackish up to the Almyros spring. The depth of the horizontal conduit is more than 400 m. An attempt at raising the water level at the spring, with a concrete dam, made in 1987, which was also modelled, indicates that the real depth is around 500 m but the poor quality of these data requires new tests to be made before any firm conclusions on the exact depth of the conduit can be drawn. The Almyros spring is a particularly favorable for observing the exchanges in the conduit network for which it is the direct outlet but it is not representative of the surrounding area. To sustainably manage the water in this region, it is essential to change the present working of the wells in order to limit the irreversible saline intrusion into the terrain of the upper aquifers. It seems possible to exploit the spring directly if the level of its outlet is raised. This would reduce the salinity in the spring to almost zero in all seasons by increasing the head in the conduit. In its present state of calibration, the model calculates a height on the order of 15 m for obtaining freshwater at the spring throughout the year, but real tests with the existing dam are needed to quantify any flow-rate losses or functional changes when there is continual overpressure in the system. The cause of the development of this karstic conduit at such a great depth could be the lowering of the sea level during the Messinian [Clauzon et al., 1996], or recent tectonic movements

Eogenetic karst from the perspective of an equivalent porous medium, 2002, Vacher H. L. , Mylroie J. E. ,
The porosity of young limestones experiencing meteoric diagenesis in the vicinity of their deposition (eogenetic karst) is mainly a double porosity consisting of touching-vug channels and preferred passageways lacing through a matrix of interparticle porosity. In contrast, the porosity of limestones experiencing subaerial erosion following burial diagenesis and uplift (telogenetic karst) is mainly a double porosity consisting of conduits within a network of fractures. The stark contrast between these two kinds of karst is illustrated by their position on a graph showing the hydraulic characteristics of an equivalent porous medium consisting of straight, cylindrical tubes (n-D space, where n is porosity, D is the diameter of the tubes, and log n is plotted against log D). Studies of the hydrology of small carbonate islands show that large-scale, horizontal hydraulic conductivity (K) increases by orders of magnitude during the evolution of eogenetic karst. Earlier petrologic studies have shown there is little if any change in the total porosity of the limestone during eogenetic diagenesis. The limestone of eogenetic karst, therefore, tracks horizontally in n-D space. In contrast, the path from initial sedimentary material to telogenetic karst comprises a descent on the graph with reduction of n during burial diagenesis, then a sideways shift with increasing D due to opening of fractures during uplift and exposure, and finally an increase in D and n during development of the conduits along the fractures. Eogenetic caves are mainly limited to boundaries between geologic units and hydrologic zones: stream caves at the contact between carbonates and underlying impermeable rocks (and collapse-origin caves derived therefrom); vertical caves along platform-margin fractures; epikarst; phreatic pockets (banana holes) along the water table; and flank margin caves that form as mixing chambers at the coastal freshwater-saltwater 'interface'. In contrast, the caverns of telogenetic karst are part of a system of interconnected conduits that drain an entire region. The eogenetic caves of small carbonate islands are, for the most part, not significantly involved in the drainage of the island

Efficient hydrologic tracer-test design for tracer-mass estimation and sample-collection frequency, 2. Experimental results, 2002, Field Ms,
Effective tracer-test design requires that the likely results be predicted in advance of test initiation to ensure tracer-test success. EHTD-predicted breakthrough curves (BTCs) for various hydrological conditions were compared with measured BTCs obtained from actual tracer tests. The hydrological conditions for the tracer tests ranged from flowing streams to porous-media systems. Tracer tests evaluated included flowing streams tracer tests conducted in small and large surface-water streams, a karst solution conduit, and a glacial-meltwater stream and porous-media systems conducted as natural-gradient, forced-gradient, injection-withdrawal, and recirculation tracer tests. Comparisons between the actual tracer tests and the predicted results showed that tracer breakthrough, hydraulic characteristics, and sample-collection frequency may be forecasted sufficiently well in most instances as to facilitate good tracer-test design. Comparisons were generally improved by including tracer decay and/or retardation in the simulations. Inclusion of tracer decay in the simulations also tended to require an increase in set average tracer concentration to facilitate matching peak concentrations in the measured BTCs, however. Both nonreactive tracer and reactive tracer predictions produced recommended sample-collection frequencies that would adequately define the actual BTCs, but estimated tracer-mass estimates were less precise

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