Eogenetic karst hydrology: insights from the 2004 hurricanes, peninsular Florida

Eogenetic karst lies geographically and temporally close to the depositional environment of limestone in warm marine water at low latitude, in areas marked by midafternoon thunderstorms during a summer rainy season. Spring hydrographs from such an environment in north-central Florida are characterized by smooth, months-long, seasonal maxima. The passage of Hurricanes Frances and Jeanne in September 2004 over three field locations shows how the eogenetic karst of the Upper Floridan Aquifer responds to unequivocal recharge events. Hydrographs at wells in the High Springs area, Rainbow Springs, and at Morris, Briar, and Bat Caves all responded promptly with a similar drawn-out rise to a maximum that extended long into the winter dry season. The timing indicates that the typical hydrograph of eogenetic karst is not the short-term fluctuations of springs in epigenic, telogenetic karst, or the smoothed response to all the summer thunderstorms, but rather the protracted response of the system to rainfall that exceeds a threshold. The similarity of cave and noncave hydrographs indicates distributed autogenic recharge and a free communication between secondary porosity and permeable matrix-both of which differ from the hydrology of epigenic, telogenetic karst. At Briar Cave, drip rates lagged behind the water table rise, suggesting that recharge was delivered by fractures, which control the cave's morphology. At High Springs, hydrographs at the Santa Fe River and a submerged conduit apparently connected to it show sharp maxima after the storms, unlike the other cave hydrographs. Our interpretation is that the caves, in general, are discontinuous.     

Natural responses of Neoproterozoic dynamic karst springs to rainfall events,São Miguel Watershed,Minas Gerais,Brazil

Karst aquifers consist of complex networks of conduits in which groundwater flows and recharge/discharge processes are generally more dynamic than in other types of aquifers. Due to their intrinsic heterogeneity and anisotropy, monitoring, quantifying, and analysing natural responses of karst springs is an efficient tool. Unlike Cenozoic and Mesozoic rocks, in Neoproterozoic karst systems, groundwater circulates and stores generally in dissolution features known as tertiary porosity, as the rock's primary porosity is recrystallized, considered negligible. This article studies the hydrodynamics of a karst portion of the São Miguel River basin, southwest of the state of Minas Gerais, Brazil. The region is predominantly composed of Neoproterozoic carbonate rocks, dating from about 570 to 540 million years ago. During a hydrological year (2019–2020), three karst springs (S1, S2, and S3) were daily monitored through their natural responses (variations of electrical conductivity, EC, temperature, T, and discharge, Q) to rainfall episodes. The data were interpreted based on the analysis of spring hydrographs, time series, recession curves (seasonal and intra-annual), and statistics of EC, T, and Q variations. The results show the three springs generally exhibit quick flow, typically karstic, in the case of hydrosystems with a well-structured and functional underground drainage network. The time series indicate the hydrosystem drained by S1 presents slower circulation and a lower degree of linearity, resulting from the higher sinuosity of the system, while the hydrosystems of S2 and S3 have similar behaviours, of quick water circulations immediately after a rainy episode. The degrees of karstification classify S1 and S2 as complex and extensive karst systems consisting of several subsystems, and S3 as a system in which the conduit network is more developed at the upper epiphreatic zone than near the outlet.     

Formation of phreatic caves in an eogenetic karst aquifer by CO2 enrichment at lower water tables and subsequent flooding by sea level rise

Formation of extensive phreatic caves in eogenetic karst aquifers is widely believed to require mixing of fresh and saltwater. Extensive phreatic caves also occur, however, in eogenetic karst aquifers where fresh and saltwater do not mix, for example in the upper Floridan aquifer. These caves are thought to have formed in their modern settings by dissolution from sinking streams or by convergence of groundwater flow paths on springs. Alternatively, these caves have been hypothesized to have formed at lower water tables during sea level low‐stands. These hypotheses have not previously been tested against one another. Analyzing morphological data and water chemistry from caves in the Suwannee River Basin in north‐central Florida and water chemistry from wells in the central Florida carbonate platform indicates that phreatic caves within the Suwannee River Basin most likely formed at lower water tables during lower sea levels. Consideration of the hydrological and geochemical constraints posed by the upper Floridan aquifer leads to the conclusion that cave formation was most likely driven by dissolution of vadose CO₂ gas into the groundwater. Sea level rise and a wetter climate during the mid‐Holocene lifted the water table above the elevation of the caves and placed the caves tens of meters below the modern water table. When rising water tables reached the land surface, surface streams formed. Incision of surface streams breached the pre‐existing caves to form modern springs, which provide access to the phreatic caves. Phreatic caves in the Suwannee River Basin are thus relict and have no causal relationship with modern surficial drainage systems. Neither mixing dissolution nor sinking streams are necessary to form laterally extensive phreatic caves in eogenetic karst aquifers. Dissolution at water tables, potentially driven by vadose CO₂ gas, offers an underappreciated mechanism to form cavernous porosity in eogenetic carbonate rocks. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of karst aquifer spring flows with a gray system decomposition model

There are approximately 470,000 km(2) of karst aquifers that feed many large springs in North China. Turbulent flow often exists in these karst aquifers, which means that the classical ground water model based on Darcy's law cannot be applied here. Ground water data are rare for these aquifers. As a consequence, it is difficult to quantitatively investigate ground water flow in these karst systems. The purpose of this study is to develop a parsimonious model that predicts karst spring discharge using gray system theory. In this theory, a white color denotes a system that is completely characterized and a black color represents a system that is totally unknown. A gray system thus describes a complex system whose characteristics are only partially known or known with uncertainty. Using this theory, we investigated the karst spring discharge time series over different time scales. First, we identified three specific components of spring discharge: the long-term trend, periodic variation, and random fluctuation. We then used the gray system model to simulate the long-term trend and obtain periodic variation and random fluctuation components. Subsequently, we developed a predictive model for karst spring discharge. Application of the model to Liulin Springs, a representative example of karst springs in northern China, shows that the model performs well. The predicted results suggest that the Liulin Springs discharge will likely decrease over time, with small fluctuations.     

Quantifying the role of karstic groundwater in a snowmelt-dominated hydrologic system

River discharge in mountainous regions of the world is often dominated by snowmelt, but base flows are sustained primarily by groundwater storage and discharge. Although numerous recent studies have focused on base-flow discharge in mountain systems, almost no work has explicitly investigated the role of karst groundwater in these systems across a full range of flow conditions. We directly measured groundwater discharge from 48 karst springs in the Kaweah River and its five forks in the Sierra Nevada mountains, California, United States. Relationships between spring and river discharge showed that karst aquifers and springs provide significant storage and delayed discharge to the river. Regression models showed that, of all potential seasonal groundwater storage compartments in the river basin, the area of karst (0.1–4.4%) present provides the best explanation of base-flow recession in each fork of the Kaweah River (directly measured contributions from karst springs ranged from 3.5 to 16% during high-flow to 20 to 65% during base-flow conditions). These results show that, even in settings where karst represents a small portion of basin area, it may play an over-sized role in seasonal storage and water resources in mountain systems. Karst aquifers are the single most important non-snow storage component in the Kaweah River basin, and likely provide similar water storage capacities and higher base flows in other mountain river systems with karst when compared with systems without karst.     

Distribution functions of spring discharges according to their lithologies and the influence of lower limit to flow: an example from Spain

Specific functions of distribution of the number of springs and of their contributions according to flow for each of nine lithological groups are established. These functions confirm the suitability of a general type of distribution (Sanz 1996), which is a borderline case of the lognormal distribution. The statistical methods proposed for the lithological types presented can be applied to any region if data on flow and lithology of the aquifers drained by springs are available. These methods have been applied to Spain, a representative region with varied geology, climate and topography; 71.2% of spring flow is supplied by limestones, 19.17% by alluvial sediments and marls, 6.7% by conglomerates and sandstones, and 3% by slates, plutonic rocks, quartzites, and other groups. Springs with discharge rates exceeding 2000 L/s exist only in limestones. The majority of springs with low flow occur in marls. If we consider springs with flow greater than 0.5 L/s, alluvial sediments, sandstones, marls, and limestones have the greatest density of springs per surface area, although the average flow varies greatly from one lithology to the next. The cumulative estimated discharge and the estimated number of springs depend upon the definition of the lower limit of flow. For example, total discharge is reduced by 30% if we eliminate all springs with flow lower than 0.5 L/min.     

Using periodic hydrologic and geochemical sampling with limited continuous monitoring to characterize remote karst aquifers in the Kaweah River Basin,California, USA

Hydrogeologic field work in remote settings is often challenging: assessing spring behaviour and aquifer characteristics can be expensive in both time commitment and resources needed to assess these systems. In this study, we document the hydrology and geochemistry of 47 perennial karst springs in the Kaweah River, a mountain river basin in the Sierra Nevada, California. After preliminary hydrogeochemical characterization and grouping, selected springs were continuously monitored to further assess aquifer characteristics in each group. Later, in areas without previous dye‐tracing work, traces were conducted to establish connections between large sinking streams and springs. The springs have a wide range of inter‐spring and intra‐spring variability in discharge and geochemistry. We assessed this variability by performing statistical comparisons with spring chemistry and principal components analysis of all measured variables. Results show that springs can be divided into two distinct groups: high elevation springs of the Mineral King Valley and lower elevation springs throughout the rest of the basin. Continuous discharge, temperature and specific conductivity data from four springs (two from each group) were then used to characterize the hydrograph recession behaviour of springs in each group. Both groups showed statistically similar baseflow recession slopes, suggesting that both groups contain baseflow storage compartments with similar hydrogeologic properties. The biggest difference between each group is the variability in amount of water remaining in the aquifer during baseflow conditions. High elevation springs have lower baseflow discharges, relative to peak flow, than lower elevation springs, despite the fact that more precipitation falls at higher elevation. This is likely caused by differences in the amount of soil and epikarst storage, which are related to recent geomorphic events: high elevation aquifers were glaciated as recent as 41 thousand years ago (kya), while there is no evidence that low elevation aquifers were glaciated. As a result, lower elevations have developed thicker soils, weathered bedrock and epikarst. Copyright © 2016 John Wiley & Sons, Ltd.     

Interpretation of groundwater level monitoring results in karst aquifers: examples from the Dinaric karst

The paper presents an attempt to determine the characteristics of karst aquifers using information on groundwater level (GWL) in natural holes and boreholes with different data quantity and time resolution of GWL measurements. In this paper the particulars of karst aquifers were analysed for four examples from the Dinaric karst. In all four study areas, aquifers are formed in bare, deep and well‐developed Dinaric karst consisting of Cretaceous limestones. The first example represents a wide area of Imotsko polje in the karst. The aquifer was analysed on the basis of infrequent water level monitoring in natural karst water features (jamas, lakes, wells) and discharges of springs and rivers. The karst aquifer in this example is complex, non‐homogenous and variable in space and time, which is frequent in the Dinaric karst. Regardless of the aforementioned it was possible to determine its elementary characteristics. The second example represents 10 wells used for the water supply for the city of Pula. The GWL and salinity were measured once a week in the period between 1981 and 1996. Even though these measurements were relatively infrequent in space and time, they served as bases for assessment of average and maximum aquifer conditions as well as boundaries of saltwater intrusion. In the third example only a portion of aquifer of the karst spring Blaz, which is in the contact with the Adriatic Seas, has been analyzed. It is a spring with an intrusion of salt water. For purposes of study of saltwater intrusion, 26 piezometers were drilled in its vicinity in which GWL, salinity and temperature were measured once a day during 168 days, a period comprising one complete cycle of seawater intrusion and retreat. These measurements proved the existence of dispersed discharge from the aquifer into the sea and its non‐homogeneity in space. In the fourth example GWL was measured continuously in 10 deep (up to 300 m) piezometers in the hinterland of the Ombla Spring catchment. The measurement period lasted 2 years (January 1988 to December 1989). The analyses are made with hourly data. The results made it possible to determine numerous characteristics of the karst aquifer and a significant non‐homogeneity of groundwater distribution in karst aquifers, depending more on the underground karst phenomena than the surface karst forms. Copyright © 2000 John Wiley & Sons, Ltd.     

Improved understanding of dynamic water and mass budgets of high-alpine karst systems obtained from studying a well-defined catchment area

Large areas of Europe, especially in the Alps, are covered by carbonate rocks and in many alpine regions, karst springs are important sources for drinking water supply. Because of their high variability and heterogeneity, the understanding of the hydrogeological functioning of karst aquifers is of particular importance for their protection and utilisation. Climate change and heavy rainfall events are major challenges in managing alpine karst aquifers which possess an enormous potential for future drinking water supply. In this study, we present research from a high-alpine karst system in the UNESCO Biosphere Reserve Großes Walsertal in Austria, which has a clearly defined catchment and is drained by only one spring system. Results show that (a) the investigated system is a highly dynamic karst aquifer with distinct reactions to rainfall events in discharge and electrical conductivity; (b) the estimated transient atmospheric CO₂ sink is about 270 t/a; (c) the calculated carbonate rock denudation rate is between 23 and 47 mm/1000a and (d) the rainfall-discharge behaviour and the internal flow dynamics can be successfully simulated using the modelling package KarstMod. The modelling results indicate the relevance of matrix storage in determining the discharge behaviour of the spring, particularly during low-flow periods. This research and the consequent results can contribute and initiate a better understanding and management of alpine karst aquifers considering climate change with more heavy rainfall events and also longer dry periods.     

A high‐altitude temporary spring in a compartmentalized carbonate aquifer: the role of low‐permeability faults and karst conduits

The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd.     

Interpretation of spring recession curves

Recession curves contain information on storage properties and different types of media such as porous, fractured, cracked lithologies and karst. Recession curve analysis provides a function that quantitatively describes the temporal discharge decay and expresses the drained volume between specific time limits (Hall 1968). This analysis also allows estimating the hydrological significance of the discharge function parameters and the hydrological properties of the aquifer. In this study, we analyze data from perennial springs in the Judean Mountains and from others in the Galilee Mountains, northern Israel. All the springs drain perched carbonate aquifers. Eight of the studied springs discharge from a karst dolomite sequence, whereas one flows out from a fractured, slumped block of chalk. We show that all the recession curves can be well fitted by a function that consists of two exponential terms with exponential coefficients alpha1 and alpha2. These coefficients are approximately constant for each spring, reflecting the hydraulic conductivity of different media through which the ground water flows to the spring. The highest coefficient represents the fast flow, probably through cracks, or quickflow, whereas the lower one reflects the slow flow through the porous medium, or baseflow. The comparison of recession curves from different springs and different years leads to the conclusion that the main factors that affect the recession curve exponential coefficients are the aquifer lithology and the geometry of the water conduits therein. In normal years of rainy winter and dry summer, alpha1 is constant in time. However, when the dry period is longer than usual because of a dry winter, alpha1 slightly decreases with time.     

Temporal cycles of karst denudation in northwest Georgia,U.S.A.

Time patterns of karst denudation in northwest Georgia (U.S.A.) were investigated at three spring sites for 12 months and at five stream sites for 10 years. Rainfall was evenly distributed and showed no significant seasonality. At the springs, as well as the streams, water hardness was largely controlled by discharge. At the springs, soil pCO₂ and water pH were strongly correlated (r + -0·69 to -0·83). Solute transport in spring waters was highly seasonal, with two conduit flow springs removing more limestone in the winter, and the diffuse flow spring removing more during the growing season. At the stream sites, most denudation occurred during the winter and spring seasons, and least during the summer. Fourier analysis showed that variations in denudation occur on deterministic (long-wave) as well as stochastic (shortwave) time scales. As contributing variables, discharge varied in short-wave and long-wave cycles, whereas soil pCO₂ showed only a long-wave cycle. The 12 month deterministic cycles were the most important, with changes in discharge taking precedence over soil pCO₂. Time series regression explains up to 69 per cent of changes in denudation through rain and soil pCO₂. Time cycles in available water are the key controlling factor of denudation, and amounts of available soil CO₂ may not be as important in the temporal patterns of karst downwearing as has been believed previously.     

Variation of karst spring discharge in the recent five decades as an indicator of global climate change: A case study at Shanxi,northern China

As one of the largest international scientific pro- grams in geoscience and environmental science, global change studies were initiated in the early 1980s[1,2]. Noticeable achievements have been made in the stud- ies using indicators such as loess, marine sediment, permafrost, vermicular red earth, and even magmatic activity[2―6]. In recent years, the importance of ground- water as a new type of global change indicators has caused wide attention[7]. Stochastic, isotopic and hy- drochemical st…     

Sustained post-seismic effects on groundwater flow in fractured carbonate aquifers in Central Italy

A sustained increase in spring discharges was monitored after the 2016 Central Italy seismic sequence in the fractured carbonate aquifer of Valnerina–Sibillini Mts. The groundwater surplus recorded between August 2016 and November 2017 was determined to be between 400 and 500 × 10⁶ m³. In fractured aquifers, the post-seismic rise in spring discharges is generally attributed to an increase in bulk permeability caused by the fracture cleaning effect, which is induced by pore pressure propagation. In the studied aquifers, the large amount of additional discharge cannot only be attributed to the enhanced permeability, which was evaluated to be less than 20% after each main seismic event. A detailed analysis of the spring discharge hydrographs and of the water level at five gauging stations was carried out to determine the possible causes of this sudden increase in groundwater outflow. Taking into account the geological and structural framework, a conceptual model of a basin-in-series has been adopted to describe the complex hydrogeological setting, where the thrusts and extensional faults have clearly influenced the groundwater flow directions before and after the seismic sequence. The prevalent portion of the total post-seismic discharge surplus not explained by the increase in permeability has been attributed to changes in the hydraulic gradient that caused seismogenic fault rupture and the disruption in the upgradient sector of the aquifer. The additional flow calculated through the breach of the pre-existing hydrostructural barrier corresponds to approximately 470 × 10⁶ m³. This value is consistent with the total discharge increase measured in the whole study area, validating the proposed conceptual model. Consequently, a shift in the piezometric divide of the hydrogeological system has been induced, causing a potentially permanent change that lowers the discharge amount of the eastern springs.     

Transport and Variability of Fecal Bacteria in Carbonate Conglomerate Aquifers

Clastic sedimentary rocks are generally considered non‐karstifiable and thus less vulnerable to pathogen contamination than karst aquifers. However, dissolution phenomena have been observed in clastic carbonate conglomerates of the Subalpine Molasse zone of the northern Alps and other regions of Europe, indicating karstification and high vulnerability, which is currently not considered for source protection zoning. Therefore, a research program was established at the Hochgrat site (Austria/Germany), as a demonstration that karst‐like characteristics, flow behavior, and high vulnerability to microbial contamination are possible in this type of aquifer. The study included geomorphologic mapping, comparative multi‐tracer tests with fluorescent dyes and bacteria‐sized fluorescent microspheres, and analyses of fecal indicator bacteria (FIB) in spring waters during different seasons. Results demonstrate that (1) flow velocities in carbonate conglomerates are similar as in typical karst aquifers, often exceeding 100 m/h; (2) microbial contaminants are rapidly transported toward springs; and (3) the magnitude and seasonal pattern of FIB variability depends on the land use in the spring catchment and its altitude. Different groundwater protection strategies that currently applied are consequently required in regions formed by karstified carbonatic clastic rocks, taking into account their high degree of heterogeneity and vulnerability.     

Method for calibrating a theoretical model in karst springs: an example for a hydropower station in South China

A theoretical, dimensionless rainfall–runoff model was used to simulate the discharge of Wulongdong spring in western Hubei province, South China. The single parameter (time constant τ) in the model is easy to obtain by fitting the recession rate of the observed hydrographs. The model was scaled by simply matching the total annual flow volume of the model to the observed value. Annual distribution of actual evapotranspiration was embedded in the model input to calculate the accumulated deficit of soil moisture before each rain event. Hourly precipitation input data performed better than daily data, defining τ of 0.85 days and returning a Nash–Sutcliffe efficiency of 0.89 and the root mean square error of 0.07. This model offers an effective way to simulate the discharge of karst springs that respond sensitively to rainfall events. The model parameters of a successful simulation can be used to estimate the recharge area and indicate the intrinsic response time of the basin. Copyright © 2016 John Wiley & Sons, Ltd.     

Dissolution in a variably confined carbonate platform: effects of allogenic runoff,hydraulic damming of groundwater inputs,and surface–groundwater exchange at the basin scale

In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SI_CAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.     

Inverse Modeling Approach to Allogenic Karst System Characterization

Allogenic karst systems function in a particular way that is influenced by the type of water infiltrating through river water losses, by karstification processes, and by water quality. Management of this system requires a good knowledge of its structure and functioning, for which a new methodology based on an inverse modeling approach appears to be well suited. This approach requires both spring and river inflow discharge measurements and a continuous record of chemical parameters in the river and at the spring. The inverse model calculates unit hydrographs and the impulse responses of fluxes from rainfall hydraulic head at the spring or rainfall flux data, the purpose of which is hydrograph separation. Hydrograph reconstruction is done using rainfall and river inflow data as model input and enables definition at each time step of the ratio of each component. Using chemical data, representing event and pre-event water, as input, it is possible to determine the origin of spring water (either fast flow through the epikarstic zone or slow flow through the saturated zone). This study made it possible to improve a conceptual model of allogenic karst system functioning. The methodology is used to study the Bas-Agly and the Cent Font karst systems, two allogenic karst systems in Southern France.     

Analytical and numerical models to explain steady rates of spring flow

Flow from some springs in former glacial lakebeds of the Upper Midwest is extremely steady throughout the year and does not increase significantly after precipitation events or seasonal recharge. Analytical and simplified numerical models of spring systems were used to determine whether preferential ground water flow through high-permeability features in shallow sandstone aquifers could produce typical values of spring discharge and the unusually steady rates of spring flow. The analytical model is based on a one-dimensional solution for periodic ground water flow. Solutions to this model suggest that it is unlikely that a periodic forcing due to seasonal variations in areal recharge would propagate to springs in a setting where high-permeability features exist. The analytical model shows that the effective length of the aquifer, or the length of flowpaths to a spring, and the total transmissivity of the aquifer have the greatest potential to impact the nature of spring flow in this setting. The numerical models show that high-permeability features can influence the magnitude of spring flow and the results demonstrate that the lengths of ground water flowpaths increase when high-permeability features are explicitly modeled, thus decreasing the likelihood for temporal variations in spring flow.     

Hydrodynamic behaviour of karstic aquifers in Boroujerd,western Iran

Abstract

Considering the geological conditions of the southwest of Boroujerd and northwest of Doroud, Iran, karst development is analysed with respect to the hydrodynamic behaviour of the main draining springs of the units and the karstic aquifers are classified as either those developed in Cretaceous limestone or those developed in Oligomiocene limestone. For this purpose, the yields of the main karstic springs of the region—Absardeh and Zoorabad (Cretaceous karstic limestone aquifer), Kalamsooz and Azizabad (Oligomiocene karstic limestone aquifer)—were measured and analysed. Analysis of the recession curve is used for hydrodynamical analysis and to construct the conceptual model for estimation of karst development in the selected aquifers. Based on the results, the dynamic storage capacity of the saturated zone in Cretaceous limestone is evaluated as low to medium and that in Oligomiocene limestone as medium to high. The dynamic storage capacity of the unsaturated zone in Cretaceous limestone is evaluated as high and that in Oligomiocene limestone as low to medium. Moreover, the contribution of quickflow in karstic aquifers developed in the Cretaceous limestone drained by the Absardeh and Zoorabad springs is 23.5 and 82.2%, respectively, and that for the Kalamsooz and Azizabad springs (Oligomiocene limestone) is 5.7 and 22.5%, respectively. Flow in the Cretaceous limestone aquifer drained by the Zoorabad Spring is of concentrated type and the main flow occurs in the well-developed karstic conduits. The main flow in the Oligomiocene limestone aquifer, drained by the Kalamsooz Spring, occurs in a network of joints and fractures and the contribution of concentrated flow is very low. The transmissivity of the saturated zone in the karstic aquifer drained by the Zoorabad and Absardeh springs is medium to high and that for the Kalamsooz and Azizabad springs is found to be low.