MWH Global

Enviroscan Ukrainian Institute of Speleology and Karstology

Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 5
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 uniformity coefficient is a numerical expression of the variety in particle sizes in mixed natural soils, defined as the ratio of the sieve size on which 40% (by weight) of the material is retained to the sieve size on which 90% of the material is retained [6].?

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

Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 7
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 flux (Keyword) returned 160 results for the whole karstbase:
Showing 31 to 45 of 160
Tectonic network as the initial factor of karstification of the chalk limestones in the Perche hills (Orne, Normandy, France), 1999,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Rodet J. ,
Tectonics has been known to be the initial factor of karstification for a long time, but this has never been demonstrated because of the progressive nature of karstification, which destroys evidence of the preliminary phases. Because of its dual hydrological quality, chalk limestone has conserved a great variety of specific forms illustrating various steps in evolution, from the initial primokarst to the classic well-developed drain. In the Perche hills, very recent research has shown specific endokarstic galleries, filled up by sands and clays. Physical, chemical and X-ray analyses have shown that such infillings result not from a classical fluviatile deposit process, but on the contrary, from an alteration process, underlining the progress of a weathering front into the limestone mass down to the watertable. The infillings are not sedimentary deposits, but represent a kind of 'shadow rock', a chemical in situ transformation of limestone, without any transport of the solid fraction, The genesis results from the progression of weathering fronts located on the tectonic pattern into the input karst. When the weathering front crosses the watertable, the resulting water mixing produces a chemical reaction capable of opening the original joint. This results in infilled galleries, similar to classical karst drains, which have never known fluviatile drainage. This illustrates the first step in normal karstic evolution, just before a water flux drains the gallery, resulting in a new karstification step: the drained passage. The specificity of the region of the Perche hills is this karat evolution stopped in the first step, illustrating the 'primorkarst' modelling by the researcher, which has never been described before. The conditions for this genesis are a well-developed tectonic pattern, an absence of thick superficial layers, a high soil permeability, and an absence of superficial drainage (lack of sinkholes). (C) Elsevier, Paris

Bridging the gap between real and mathematically simulated karst aquifers, 1999,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Groves C. , Meiman J. , Howard A. D.
Although several numerical codes have been developed to study the patterns of karst aquifer evolution and behavior, in the current generation of models simplifying assumptions must be made because of incomplete quantitative understanding of key processesA one-year, high-temporal-resolution study of carbonate chemistry with Mammoth Cave's Logsdon River, designed to investigate details of these processes, reveals that limestone dissolution rates vary appreciably over storm and seasonal time scales due to variations in the flux of CO2-rich waters that wash through, and flood, conduits during storm eventsThis undersaturated storm water dissolves rock within a flood zone 25-30 m thickThrough the year, waters were undersaturated only 31% of the timeTime scales of actual karst development may thus be impacted by time-varying processes different from the constant-input chemistry assumed in current published numerical codesA dual approach, coupling quantitative modeling and refinement of the models by careful measurement of processes within real karst aquifers, provides a framework for developing a comprehensive understanding of karst system behavior

Conduit hydrogeology of a tropical coastal carbonate aquifer. MSc thesis, 1999,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Beddows, P. A

The aim of this study is to investigate the hydrogeology of the submerged conduit systems of a coastal carbonate aquifer (Caribbean coast, Yucatan Peninsula, Mexico), and thereby better understand their significance as large permeability heterogeneities. A complex spatial trend in conduit flow rates (determined by quantitative fluorescent dye tracing initiated 6 km inland) was found, including significant velocity variation between consecutive conduit segments. Elevated coastal velocities under low tide conditions are shown by salinity profiling, to be induced by the volumetric increase of discharging water from mixing with marine water. Semi-diurnal micro-tidal loading is sufficient to induce flooding from the sea into the conduits at one coastal discharge point, and significantly reduce flow rates at another. Furthermore, a network of four observation sites extending 5 km inland indicates efficient propagation of the ~0.30 m tidal signal through the Nohoch Nah Chich conduit system, a distance several time greater than previously appreciated in this environment. The field results clearly indicate that the hydrogeological flux is dominated by cavernous porosity, and that the aquifer is dynamically responsive to the high-frequency low-magnitude tidal loading to a significant distance inland. Conventional coastal groundwater models such as the Ghyben-Herzberg lens model, assume isotropic homogeneous equivalent-porous-medium conditions. Because the corollaries of the conventional models are inconsistent with the field evidence, they are inapplicable in this environment. It is hoped that these results will aid future modelling efforts, and improve our capacity to manage the valuable groundwater resources which represents the unique source of potable water to the local population.

Symposium Abstract: Carbon flux and aquifer evolution in the South-Central Kentucky karst [USA], 2000,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Groves C. , Meiman J.

Hydrochemical Interpretation of Cave Patterns in the Guadalupe Mountains, New Mexico, 2000,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Palmer, A. N. , Palmer, M. V.
Most caves in the Guadalupe Mountains have ramifying patterns consisting of large rooms with narrow rifts extending downward, and with successive outlet passages arranged in crude levels. They were formed by sulfuric acid from the oxidation of hydrogen sulfide, a process that is now dormant. Episodic escape of H2S-rich water from the adjacent Delaware Basin, and perhaps also from strata beneath the Guadalupes, followed different routes at different times. For this reason, major rooms and passages correlate poorly between caves, and within large individual caves. The largest cave volumes formed where H2S emerged at the contemporary water table, where oxidation was most rapid. Steeply ascending passages formed where oxygenated meteoric water converged with deep-seated H2S-rich water at depths as much as 200 m below the water table. Spongework and network mazes were formed by highly aggressive water in mixing zones, and they commonly rim, underlie, or connect rooms. Transport of H2S in aqueous solution was the main mode of H2S influx. Neither upwelling of gas bubbles nor molecular diffusion appears to have played a major role in cave development, although some H2S could have been carried by less-soluble methane bubbles. Most cave origin was phreatic, although subaerial dissolution and gypsum-replacement of carbonate rock in acidic water films and drips account for considerable cave enlargement above the water table. Estimates of enlargement rates are complicated by gypsum replacement of carbonate rock because the gypsum continues to be dissolved by fresh vadose water long after the major carbonate dissolution has ceased. Volume-for-volume replacement of calcite by gypsum can take place at the moderate pH values typical of phreatic water in carbonates, preserving the original bedrock textures. At pHs less than about 6.4, this replacement usually takes place on a molar basis, with an approximately two-fold volume increase, forming blistered crusts.

Is it always dark in caves?, 2000,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Badino Giovanni
Underground natural sources of visible light are considered. The main light producer is Cerenkov radiation emitted in air, water and rock by cosmic ray muons, that depends, in a complex way, on shape of mountain and of caves. In general the illumination increases linearly with the cavity dimensions. Other light sources are from secondary processes generated by radioactive decays in rock from minerals luminescence. The natural light fluxes in caves are in general easy to detect but are not used from underground life.

Carbonate platform systems: components and interactions -- an introduction, 2000,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Insalaco Enzo, Skelton Peter, Palmer Tim J. ,
Carbonate platforms are open systems with natural boundaries in space and time. Across their spatial boundaries there are fluxes of energy (e.g. light, chemical energy in compounds, and kinetic energy in currents and mass flows) and matter (e.g. nutrients, dissolved gases such as CO2, and sediment -- especially, of course, carbonates). Internally, these fluxes are regulated by myriads of interactions and feedbacks (Masse 1995), and the residue is consigned to the geological record. The most distinctive aspect of carbonate platforms is the predominant role of organisms in producing, processing and/or trapping carbonate sediment, even in Precambrian examples. Because of evolutionary changes in this strong biotic input, it is harder to generalize about carbonate platforms than about most other sedimentary systems. Evolution has altered both the constructive and destructive effects of platform-dwelling organisms on carbonate fabrics, with profound consequences for facies development. Moreover, changing patterns in the provision of accommodation space (e.g. between greenhouse and icehouse climatic regimes) have also left their stamp on facies geometries, in turn feeding back to the evolution of the platform biotas. Hence simplistic analogies between modern and ancient platforms may give rise to misleading interpretations of what the latter were like and how they formed. Although a number of carbonate platform and reef specialists have warned of the dangers of such misplaced uniformitarianism (e.g. Braithwaite 1973; Gili et al. 1995; Wood 1999), it remains depressingly commonplace in the literature on ancient carbonate platforms. The endless quest in the literature for an allpurpose definition of reefs' ... This 250-word extract was created in the absence of an abstract

PCB Pollution of the Karstic Environment (Krupa River, Slovenia), 2000,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
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.

Groundwater pollution by contaminant transport from soil to fractured rock, 2000,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Witthü, Ser Kai, Č, Enč, Ur Curk Barbara

Water flow and contaminant transport from soil to underlying fractured rock is mainly controlled by the hydraulic conditions of the soil-bedrock boundary. In respect to the necessary understanding of contaminant transport at the soil-bedrock boundary the identification of flow paths within both the soil cover and the fractured media is decisive on the one side. On the other hand substance-specific behaviour of the often reactive pollutants compared to water flow has to be known in detail. Field scale tracer tests with different tracers (uranine and salts) and a potential pollutant as a reactive tracer (nitrate) were performed at the IRGO field research facility Sinji Vrh (SI). Injection points are located on the surface, in the soil, at the soil-rock interface and in the fractured rock; water is sampled in an underground tunnel with the help of two subhorizontal boreholes equipped with sampling devices and a special construction for collecting water seeping from the ceiling. The goal of these experiments is to identify the flow paths of solutes to the underground tunnel and to estimate their residence time dependent on the injection point. So far only some conclusions regarding the waterflux into the tunnel could be drawn.

Seasonal variations in Sr, Mg and P in modern speleothems (Grotta di Ernesto, Italy), 2001,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Huang Yiming, Fairchild Ian J. , Borsato Andrea, Frisia Silvia, Cassidy Nigel J. , Mcdermott Frank, Hawkesworth Chris J. ,
Sub-annual variations in trace element chemistry and luminescence have recently been demonstrated from speleothems and offer the potential of high-resolution palaeoclimatic proxies. However, no studies have yet examined microscopic trace element variations in relation to modern cave conditions. In this study, the spatial variations in trace element (Sr, Mg and P) concentrations in speleothems (a stalagmite and a soda straw stalactite) from the alpine Ernesto cave (temperature 6.60.1[deg]C) in a forested catchment in NE Italy have been studied using secondary ion mass spectrometry (SIMS) and compared with environmental parameters and waters in the modern cave. An annual lamination exists in the stalagmite and soda straw stalactite in the form of clear calcite with narrow visible layers, which are UV-fluorescent and interpreted to contain soil-derived humic/fulvic acids washed into the cave during autumn rains. Microanalyses were undertaken of seven annual laminae, probably deposited during the 1960s in the stalagmite, and seven laminae in the 1990s for the stalactite.The analysis results show that Sr consistently has a trough and P, a peak centred on the inclusion-rich layer. Mg shows mainly a negative covariation with Sr in laminae formed in the 1990s, but a positive covariation in the stalagmite formed in 1960s. The spatial scale of the main geochemical variations is the same as that of annual laminae of inclusion-poor and inclusion-rich couplets. Mass balance arguments are used to show that the P is inorganic in form and presumably occurs as individual phosphate ions within the calcite.Most drip waters show limited chemical variations, but a summer peak in trace elements in 1995 and a decrease in Mg/Ca in the following winter are notable. More pronounced covariations in Mg/Ca and Sr/Ca are shown by a site with highly variable drip rates where ratios increase at slow drip rates. The strongest seasonal variations are found in pool waters, where ratios increase reflecting significant Ca removal from the water into the calcite during the winter in response to seasonal PCO2 variations in cave air. Thus, the cave waters' compositions tend to reflect climate conditions, such that Mg/Ca and Sr/Ca are tentatively interpreted to be higher when climate conditions are dry.Combining results from the speleothems and cave water along with the behaviour of each trace species, Mg/Ca variations in the speleothems are considered to reflect their variation in the cave waters, whereas, Sr incorporation is also dependent on precipitation rate, in this case, mainly controlled by temporal variations in PCO2 in the cave (and conceivably, also by inhibitors such as phosphate). P adsorption (a fraction of which is subsequently incorporated within calcite) depends on aqueous phosphate concentration and water flux, both of which should increase during the autumn. Therefore, multiple trace element profiles in speleothems reflect multiple aspects of environment seasonality and conditions, and hence, a calibration against weather records is desirable to establish their palaeoclimatological meaning. The strong annual variation of trace elements, and particularly P, can provide chronological markers for high-resolution studies of other climate proxies, such as stable isotopes

Inverse modeling of the hydrological and the hydrochemical behavior of hydrosystems: Characterization of karst system functioning, 2001,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Pinault J. L. , Plagnes V. , Aquilina L. , Bakalowicz M. ,
Inverse modeling of mass transfer characterizes the dynamic processes affecting the function of karst systems and can be used to identify karst properties. An inverse model is proposed to calculate unit hydrographs as well as impulse response of fluxes from rainfall-runoff or rainfall-flux data, the purpose of which is hydrograph separation. Contrary to what hydrologists have been doing for years, hydrograph separation is carried out by using transfer functions in their entirety, which enables accurate separation of fluxes, as was explained in the companion paper [Pinault et al., this issue]. The unit hydrograph as well as impulse response of fluxes is decomposed into a quick and a slow component, and, consequently, the effective rainfall is decomposed into two parts, one contributing to the quick flow (or flux) and the other contributing to the slow flow generation. This approach is applied to seven French karstic aquifers located on the Larzac plateau in the Grands Causses area (in the south of France). Both hydrodynamical and hydrogeochemical data have been recorded from these springs over several hydrological cycles. For modeling purposes, karst properties can be represented by the impulse responses of flow and flux of dissolved species. The heterogeneity of aquifers is translated to time-modulated flow and transport at the outlet. Monitoring these fluxes enables the evaluation of slow and quick components in the hydrograph. The quick component refers to the 'flush flow' effect and results from fast infiltration in the karst conduit network when connection is established between the infiltration and phreatic zones, inducing an increase in water head. This component reflects flood events where flow behavior is nonlinear and is described by a very short transfer function, which increases and decreases according to water head. The slow component consists of slow and fast infiltration, underground runoff, storage in annex-to-drain systems, and discharge from the saturated zone. These components can be further subdivided by measuring chemical responses at the karst outlet. Using Such natural tracers enables the slow component of the unit hydrograph to be separated into preevent water, i.e., water of the reservoir and event water, i.e., water whose origin can be related to a particular rainfall event. These measurements can be used to determine the rate of water renewal. Since the preevent water hydrograph is produced by stored water when pushed by a rainfall event and the event water hydrograph reflects rainwater transfer, separating the two components can yield insights into the characteristics of karst aquifers, the modes of infiltration, and the mechanisms involved in karstification, as well as the degree of organization of the aquifer

Soil carbon dioxide in a summer-dry subalpine karst, Marble Mountains, California, USA, 2001,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Davis J, Amato P, Kiefer R,
Studies of the seasonality, spatial variation and geomorphic effects of Soil CO2 concentrations in a summer-dry subalpine karst landscape in the Marble Mountains, Klamath National Forest, California, demonstrate the significance of soil moisture as a limiting factor. Modeled actual evapotranspiration (AET) in the four weeks prior to sampling explains 36% of the observed soil-CO2 concentrations, pointing to the importance of root respiration processes in these systems. Late snows are significant in controlling the timing of a snowmelt-initiated pulse of respiration and groundwater. CO2 concentrations were measured at multiple sites in two seasons - 1995 and 1997 - with contrasting patterns of snowmelt. Other than wet-meadow anomalies, where CO2 concentrations reached up to 3.8% in midsummer, alpine meadows on schist were the sites of the highest spring peak concentrations of approximately 1%. Forest sites and sites with thin soils on marble typically peaked at approximately 0.5%, also within a month of snowmelt exposure. Ongoing karstification in the upper bare karst is focused in soil-filled grikes where late-season snowmelt concentrates flow during high-respiration periods, but the lack of active speleothem development suggests that the carbonate solution system is greatly reduced from preglacial periods

Inorganic carbon flux and aquifer evolution in the south central Kentucky karst, 2001,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Groves C. , Meiman J.

Tracer Study on the Tectonic Controll of the Drainage System in the Contact Karst Zone of Lake Voralp (Swiss Alps), 2001,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
Ross Janhenning, Rieg Alfred, Leibundgut Chris

Lake Voralp is a small karst lake, dammed up by a Pleistocene rockslide and without any surface drainage. A dye tracer test under flood conditions shows a widespread drainage system to remote springs and direct afflux to the porous aquifer of the Rhine Valley. These karstwater passages are leading through nonkarstified flysch and a marly layer. Both layers are lithologically rather impermeable and karstwater passages through these layers indicates water leading fault zones. Properties of the drainage patterns strongly depend on the hydrological situation.

Coastal karst springs in the Mediterranean basin : study of the mechanisms of saline pollution at the Almyros spring (Crete), observations and modelling, 2002,
Deprecated: Function get_magic_quotes_gpc() is deprecated in /home/isthin5/public_html/addon-domains/ on line 943
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 160
You probably didn't submit anything to search for