MWH Global

Enviroscan Ukrainian Institute of Speleology and Karstology

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 meteoric water is water recently involved in atmospheric circulation [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 simulation (Keyword) returned 143 results for the whole karstbase:
Showing 31 to 45 of 143
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

Computersimulationen zum Schichtenbild der Stalagmiten., 2000, Franke, H(erbert). W(erner).
[mit Computerbilder]

Charakterisierung von Karstgrundwasserleitern durch Simulation der Aquifergenese und des W_-rmetransports, 2000, Liedl Rudolf, Sauter Martin,

Simulated effect of vadose infiltration on water levels in the Northern Guam Lens Aquifer, 2000, Contractor Dn, Jenson Jw,
Regional-scale hydrology of the fresh water lens in the Northern Guam Lens Aquifer has been simulated in the past using a finite element, sharp interface computer model, SWIG2D. Systematic differences exist between observed and computed water levels. Computed seasonal peak water levels are higher, and the computed seasonal lows are lower than the respective observed levels. It is hypothesized that vadose storage must store a substantial amount of water during the wet season and release it gradually into the lens during the dry season. Flow through the vadose zone was simulated with a one-dimensional finite element, unsaturated flow program UNSATID, in which the van Genuchten model is used to characterize unsaturated diffuse flow through the matrix of the vadose zone. An additional parameter (SINK) was added to the van Genuchten set to account for rapid infiltration down open pathways (fractures) associated with the closed depressions of the karst terrain. A global-optimization technique (Shuffled Complex Evolution or SCE-UA Method) was used to obtain the parameters that minimized the difference between simulated and observed water levels. Simulations incorporating the van Genuchten model were accomplished by combining the two programs, UNSATID and SWIG2D, into a single program. The sum-of-squared-errors (SSE) between computed and observed water levels in four observation wells was minimized using SCE-UA, reducing the arithmetically averaged SSE of the four wells by 30% compared with the SSE obtained when the vadose zone was not modeled. These results suggest that vadose storage is significant. On the other hand, the fact that the best fit obtained with an optimum parameter set was able to reduce the SSE by no more than 30% suggests that additional phenomena have yet to be accounted for to mon fully explain differences between simulated and observed well water levels. (C) 2000 Elsevier Science B.V. All rights reserved

Rainfall-runoff relations for karstic springs. Part I: convolution and spectral analyses, 2000, Labat D. , Ababou R. , Mangin A. ,
Karstic basins contain large reserves of subsurface water. In this paper, three karstic systems located in the Pyrenees Mountains (Ariege, France) are studied. Long records of rainfall and discharge rates for these karstic springs are available, sampled at different rates: daily, hourly and half-hourly. This study aims at illustrating and assessing the capabilities and limitations of linear black-box methods for analysing rainfall-runoff type relationships and reconstructing runoffs from rainfall rate data using such systems. In this study, precipitation and discharge rates are considered as two autocorrelated and cross-correlated stochastic processes. A Linear and stationary rainfall-runoff model is adopted, which is used for identification and simulation purposes. Different versions are analysed, including a model based on a convolution integral between the precipitation rate P(tau) and a transfer function h(t - tau) which can be thought of as the unit impulse response of the system. It is shown that this linear stochastic model (i.e. the statistical version), although accurate in some respects, does not represent the hydraulic behaviour of the system very well during low flow episodes and floods. It is also shown that the use of Fourier analysis, alone, does not lead to a satisfactory reconstitution of observed runoff sequences. For these reasons, the use of non-linear random process input-output models based on Volterra integral series is proposed and discussed. (C) 2000 Elsevier Science B.V. All rights reserved

Climate-change impacts in a regional karst aquifer, Texas, USA, 2000, Loaiciga H. A. , Maidment D. R. , Valdes J. B. ,
Climate-change scenarios were created from scaling factors derived from several general circulation models to assess the likely impacts of aquifer pumping on the water resources of the Edwards Balcones Fault Zone (BFZ) aquifer, Texas, one of the largest aquifer systems in the United States. Historical climatic time series in periods of extreme water shortage (1947-1959), near-average recharge (1978-1989), and above-average recharge (1975-1990) were scaled to 2 x CO2 conditions to create aquifer recharge scenarios in a wanner climate. Several pumping scenarios were combined with 2 x CO2 climate scenarios to assess the sensitivity of water resources impacts to human-induced stresses on the Edwards BFZ aquifer. The 2 x CO2 climate change scenarios were linked to surface hydrology and used to drive aquifer dynamics with alternative numerical simulation models calibrated to the Edwards BFZ aquifer, Aquifer simulations indicate that, given the predicted growth and water demand in the Edwards BFZ aquifer region, the aquifer's ground water resources appear threatened under 2 x CO2 climate scenarios. Our simulations indicate that 2 x CO2 climatic conditions could exacerbate negative impacts and water shortages in the Edwards BFZ aquifer even if pumping does not increase above its present average level. The historical evidence and the results of this article indicate that without proper consideration to variations in aquifer recharge and sound pumping strategies, the water resources of the Edwards BFZ aquifer could be severely impacted under a warmer climate. (C) 2000 Elsevier Science B.V. All rights reserved

Simulation of daily and monthly stream discharge from small watersheds using the SWAT model, 2000, Spruill C. A. , Workman S. R. , Taraba J. L. ,
The Soil and Water Assessment Tool (SWAT) was evaluated and parameter sensitivities were determined while modeling daily streamflows in a small central Kentucky watershed over a two-year period. Streamflow data from 1996 were used to calibrate the model and streamflow data from 1995 were used for evaluation. The model adequately predicted the trends in daily streamflow during this period although Nash-Sutcliffe R-2 values were -0.04 and 0.19 for 1995 and 1996, respectively The model poorly predicted the timing of some peak flow values and recession rates during the last half of 1995. Excluding daily peak flow values from August to December improved the daily R-2 to 0.15, which was similar to the 1996 daily R2 value. The Nash-Sutcliffe R-2 for monthly total flows were 0.58 for 1995 and 0.89 for 1996 which were similar to values found in the literature. Since very little information was available on the sensitivity of the SWAT model to various inputs, a sensitivity analysis/calibration procedure was designed to evaluate parameters that were thought to influence stream discharge predictions. These parameters included, drainage area, slope length, channel length, saturated hydraulic conductivity, and available water capacity. Minimization of the average absolute deviation between observed and simulated streamflows identified optimum values/ranges for each parameter. Saturated hydraulic conductivity alpha baseflow factor; drainage area, channel length, and channel width were the most sensitive parameters in modeling the karst influenced watershed. The sensitivity analysis process confirmed die trace studies in the karst watershed that a much larger area contributes to streamflow than can be described by the topographic boundaries. Overall, the results indicate that the SWAT model can be an effective tool for describing monthly, runoff from small watersheds in central Kentucky that have developed on karat hydrology however calibration data are necessary to account for solution channels draining into or out of the topographic watershed

Effects of nearshore recharge on groundwater interactions with a lake in mantled karst terrain, 2000, Lee T. M. ,
The recharge and discharge of groundwater were investigated for a lake basin in the mantled karst terrain of central Florida to determine the relative importance of transient groundwater inflow to the lake water budget. Variably saturated groundwater flow modeling simulated water table responses observed beneath two hillsides radiating outward from the groundwater flow-through lake. Modeling results indicated that transient water table mounding and groundwater flow reversals in the nearshore region following large daily rainfall events generated most of the net groundwater inflow to the lake. Simulated daily groundwater inflow was greatest following water table mounding near the lake, not following subsequent peaks in the water level of upper basin wells. Transient mounding generated net groundwater inflow to the lake, that is, groundwater inflow in excess of the outflow occurring through the deeper lake bottom. The timing of the modeled net groundwater inflow agreed with an independent lake water budget; however, the quantity was considerably less than the budget-derived value

Groundwater Flow and Contaminant Transport in Carbonate Aquifers, 2000, Sasowsky I. D. , Wicks C. M.
Carbonate aquifers are an important source of water throughout the world. They are complicated systems and not always easy to interpret. Caves and channels form in the rock, leading to complex flow pathways and unpredictable contaminant behaviour. This volume covers the range of techniques used to analyse groundwater flow and contaminant transport in carbonate aquifers. The book opens with a review of thoughts and methods, and continues by discussing the use of tracers, hydrograph and hydrochemograph evaluation, estimation of aquifer properties from outcrop studies, numerical simulation, analogue simulation, and 3-D visualization of conduits. Other papers address the critical evaluation of matrix, fracture and conduit components of flow and storage. An understanding of these approaches is important to engineers or hydrogeologists working in carbonate aquifers.

Modelling karst aquifer genesis using a coupled continuum-pipe flow model, 2000, Sauter M. , Liedl R.
This contribution presents a modeling approach for the simulation of karst aquifer genesis during geologic time periods. The approach used here is based upon a hydraulic coupling of a pipe network to a continuum in order to represent the well-known duality of karst aquifer flow systems. In addition, the growth of the karst conduits due to carbonate dissolution can be simulated so that a feedback mechanism with respect to the hydraulic properties of the pipe network is established. The model is applied to some typical scenarios of karst aquifer genesis, i. e. the evolution of a dendritic tube system and of a maze cave network. It is found that the model results can successfully reproduce the basic features of karst development, e. g. the development of maze caves in a carbonate aquifer which is covered by a sandstone caprock.

PCB Pollution of the Karstic Environment (Krupa River, Slovenia), 2000, Polič, S. , Leskovš, Ek H. , Horvat M.

The PCB (Polychlorinated biphenyls) pollution problems in Krupa River are related to sinking surficial streams that mix with the regional groundwater supply, thus endangering the quality of the groundwater reservoirs. The last state of contamination with PCBs and their fate and exposure in the polluted environment were studied in the period 1995 to 1998 within the Remedial Environmental Programme that began in 1984. The concentrations of PCBs in the environmental compartments (air, water, sediment, and soil) of the Krupa were determined. A simple two-compartment (air/water) fugacity mass balance model was applied to these figures, and yielded a reasonable simulation of long-term trends in concentration. The net flux of PCBs from water to air was determined (0.3 ng/m2s). Research into levels of pollution in the Krupa area showed the situation and trends regarding pollution after ten years of remedial measures, and produced a quantitative assessment of the emission of PCBs from the underground catchment area of the source of the river. Comparisons between PCB concentrations in the polluted water measured from 1986 to 1988 and the concentrations measured in 1995 to 1997 show a decreasing trend (i.e. from a mean of 380 ng/l to a mean of 100 ng/l in the Krupa's water). The emission flux of PCBs from water to air for the entire River Krupa (2.6 km long) was estimated at 200 to500 g/year, i.e. approximately one-tenth of the figure at the time the remedial measure was initiated.

Characteristics of the recharge-discharge relation of the karst aquifer in the background of the Vipava springs (Slovenia), 2000, Petrič, Metka

In the presented research I tried to find answers to the following questions: what the characteristics of the relation between recharge and discharge of the karst aquifer in the background of the Vipava springs are and which conclusions about the functioning of the karst system can be made based on stated relations between its input and output signal. Several different models of the system recharge-discharge were set. For each the conversion of the input signal to the system response was defined by transfer functions, which express conditions and processes in the karst system in mathematical form. Based on the comparison of the suitability of these different models an important influence of vegetation and processes in atmosphere and soil on the quantity of water that actually enters the karst aquifer and is expressed as recharge function was proved. An increase in the accuracy of the simulation was obtained also by the temporal distribution of recharge in fast and slow component, which is in agreement with the double mechanism of the functioning of the system: fast flow through the karst drainage net and longer retardation in the system with storage of water in less permeable zones.

Computersimulationen zum Schichtenbild der Stalagmiten, 2000, Franke, H. W.

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

Results 31 to 45 of 143
You probably didn't submit anything to search for