首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 25 毫秒
1.
Hydrological and hydrochemical processes in the critical zone of karst environments are controlled by the fracture‐conduit network. Modelling hydrological and hydrochemical dynamics in such heterogeneous hydrogeological settings remains a research challenge. In this study, water and solute transport in the dual flow system of the karst critical zone were investigated in a 73.5‐km2 catchment in southwest China. We developed a dual reservoir conceptual run‐off model combined with an autoregressive and moving average model with algorithms to assess dissolution rates in the “fast flow” and “slow flow” systems. This model was applied to 3 catchments with typical karst critical zone architectures, to show how flow exchange between fracture and conduit networks changes in relation to catchment storage dynamics. The flux of bidirectional water and solute exchange between the fissure and conduit system increases from the headwaters to the outfall due to the large area of the developed conduits and low hydraulic gradient in the lower catchment. Rainfall amounts have a significant influence on partitioning the relative proportions of flow and solutes derived from different sources reaching the underground outlet. The effect of rainfall on catchment function is modulated by the structure of the karst critical zone (e.g., epikarst and sinkholes). Thin epikarst and well‐developed sinkholes in the headwaters divert more surface water (younger water) into the underground channel network, leading to a higher fraction of rainfall recharge into the fast flow system and total outflow. Also, the contribution of carbonate weathering to mass export is also higher in the headwaters due to the infiltration of younger water with low solute concentrations through sinkholes.  相似文献   

2.
Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ18O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.  相似文献   

3.
Karst aquifers are highly productive groundwater systems often associated with conduit flow. These systems can be highly vulnerable to contamination, resulting in a high potential for contaminant exposure to humans and ecosystems. This work develops statistical models to spatially characterize flow and transport patterns in karstified limestone and determines the effect of aquifer flow rates on these patterns. A laboratory‐scale Geo‐HydroBed model is used to simulate flow and transport processes in a karstic limestone unit. The model consists of stainless steel tanks containing a karstified limestone block collected from a karst aquifer formation in northern Puerto Rico. Experimental work involves making a series of flow and tracer injections, while monitoring hydraulic and tracer response spatially and temporally. Statistical mixed models (SMMs) are applied to hydraulic data to determine likely pathways of preferential flow in the limestone units. The models indicate a highly heterogeneous system with dominant, flow‐dependent preferential flow regions. Results indicate that regions of preferential flow tend to expand at higher groundwater flow rates, suggesting a greater volume of the system being flushed by flowing water at higher rates. Spatial and temporal distribution of tracer concentrations indicates the presence of conduit‐like and diffuse flow transport in the system, supporting the notion of both combined transport mechanisms in the limestone unit. The temporal response of tracer concentrations at different locations in the model coincide with, and confirms the preferential flow distribution generated with the SMMs used in the study.  相似文献   

4.
High‐resolution measurements of rainfall, water level, pH, conductivity, temperature and carbonate chemistry parameters of groundwater at two adjacent locations within the peak cluster karst of the Guilin Karst Experimental Site in Guangxi Province, China, were made with different types of multiparameter sonde. The data were stored using data loggers recording with 2 min or 15 min resolution. Waters from a large, perennial spring represent the exit for the aquifer's conduit flow, and a nearby well measures water in the conduit‐adjacent, fractured media. During flood pulses, the pH of the conduit flow water rises as the conductivity falls. In contrast, and at the same time, the pH of groundwater in the fractures drops, as conductivity rises. As Ca2+ and HCO3? were the dominant (>90%) ions, we developed linear relationships (both r2 > 0·91) between conductivity and those ions, respectively, and in turn calculated variations in the calcite saturation index (SIC) and CO2 partial pressure (P) of water during flood pulses. Results indicate that the P of fracture water during flood periods is higher than that at lower flows, and its SIC is lower. Simultaneously, P of conduit water during the flood period is lower than that at lower flows, and its SIC also is lower. From these results we conclude that at least two key processes are controlling hydrochemical variations during flood periods: (i) dilution by precipitation and (ii) water–rock–gas interactions. To explain hydrochemical variations in the fracture water, the water–rock–gas interactions may be more important. For example, during flood periods, soil gas with high CO2 concentrations dissolves in water and enters the fracture system, the water, which in turn has become more highly undersaturated, dissolves more limestone, and the conductivity increases. Dilution of rainfall is more important in controlling hydrochemical variations of conduit water, because rainfall with higher pH (in this area apparently owing to interaction with limestone dust in the lower atmosphere) and low conductivity travels through the conduit system rapidly. These results illustrate that to understand the hydrochemical variations in karst systems, considering only water–rock interactions is not sufficient, and the variable effects of CO2 on the system should be evaluated. Consideration of water–rock–gas interactions is thus a must in understanding variations in karst hydrochemistry. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

5.
Comparison of flowpaths to a well and spring in a karst aquifer   总被引：3，自引：0，他引：3  
Toran L  Herman EK  White WB 《Ground water》2007,45(3):281-287
The permeability of some karst aquifers consists of networks of poorly integrated conduits and dissolution-widened fractures. The flow includes conduit flow, especially during storm recharge, but lacks the focused recharge into single master conduits that occurs in more highly developed karst systems. The proportions of conduit and dispersed flow are difficult to quantify in such systems. This study examines the flowpaths in a small karst watershed, based on comparing the physical and chemical response to storm flow at both a spring and a well. By conducting continuous monitoring at both locations, a better understanding of the flowpaths in a poorly integrated network was obtained. A more permeable flowpath to the spring leads to faster storm response and lower ion concentrations. The flowpath to and from the well is more complicated. The higher ion content and slower storm response suggest slower, more dispersed flowpaths. However, the well has greater variation in ion chemistry. Periodic recharge may dilute well concentrations due to faster (conduit or fracture) flowpaths. Although karst systems such as this are difficult to characterize, applying a variety of geochemical and physical monitoring techniques at multiple locations illustrates that the flowpaths can vary in both space and time.  相似文献   

6.
Karst Aquifer GIS‐based model (KAGIS model) is developed and applied to Mela aquifer, a small karst aquifer located in a Mediterranean region (SE Spain). This model considers different variables, such as precipitation, land use, surface slope and lithology, and their geographical heterogeneity to calculate both, the run‐off coefficients and the fraction of precipitation which contributes to fill the soil water reservoir existing above the aquifer. Evapotranspiration uptakes deplete water, exclusively, from this soil water reservoir and aquifer recharge occurs when water in the soil reservoir exceeds the soil field capacity. The proposed model also obtains variations of the effective porosity in a vertical profile, an intrinsic consequence of the karstification processes. A new proposal from the Nash–Sutcliffe efficiency index, adapted to arid environments, is presented and employed to evaluate the model's ability to predict the water table oscillations. The uncertainty in the model parameters is determined by the Generalized Likelihood Uncertainty Estimation method. Afterwards, when KAGIS is calibrated, wavelet analysis is applied to the resulting data in order to evaluate the variability in the aquifer behaviour. Wavelet analysis reveals that the rapid hydrogeological response, typical of a wide variety of karst systems, is the prevailing feature of Mela aquifer. This study proves that KAGIS is a useful tool to quantify recharge and discharge rates of karst aquifers and can be effectively applied to develop a proper management of water resources in Mediterranean areas.  相似文献   

7.
The modeling of groundwater flow in karst aquifers is a challenge due to the extreme heterogeneity of its hydraulic parameters and the duality in their discharge behavior, that is, rapid response of highly conductive karst conduits and delayed drainage of the low‐permeability fractured matrix after recharge events. There are a number of different modeling approaches for the simulation of the karst groundwater dynamics, applicable to different aquifer as well as modeling problem types, ranging from continuum models to double continuum models to discrete and hybrid models. This study presents the application of an equivalent porous model approach (EPM, single continuum model) to construct a steady‐state numerical flow model for an important karst aquifer, that is, the Western Mountain Aquifer Basin (WMAB), shared by Israel and the West‐Bank, using MODFLOW2000. The WMAB was used as a catchment since it is a well‐constrained catchment with well‐defined recharge and discharge components and therefore allows a control on the modeling approach, a very rare opportunity for karst aquifer modeling. The model demonstrates the applicability of equivalent porous medium models for the simulation of karst systems, despite their large contrast in hydraulic conductivities. As long as the simulated saturated volume is large enough to average out the local influence of karst conduits and as long as transport velocities are not an issue, EPM models excellently simulate the observed head distribution. The model serves as a starting basis that will be used as a reference for developing a long‐term dynamic model for the WMAB, starting from the pre‐development period (i.e., 1940s) up to date.  相似文献   

8.
Groundwater is the principal water resource in semi‐arid and arid environments. Therefore, quantitative estimates of its replenishment rate are important for managing groundwater systems. In dry regions, karst outcrops often show enhanced recharge rates compared with other surface and sub‐surface conditions. Areas with exposed karst features like sinkholes or open shafts allow point recharge, even from single rainfall events. Using the example of the As Sulb plateau in Saudi Arabia, this study introduces a cost‐effective and robust method for recharge monitoring and modelling in karst outcrops. The measurement of discharge of a representative small catchment (4.0 · 104 m2) into a sinkhole, and hence the direct recharge into the aquifer, was carried out with a time‐lapse camera. During the monitoring period of two rainy seasons (autumn 2012 to spring 2014), four recharge events were recorded. Afterwards, recharge data as well as proxy data about the drying of the sediment cover are used to set up a conceptual water balance model. The model was run for 17 years (1971 to 1986 and 2012 to 2014). Simulation results show highly variable seasonal recharge–precipitation ratios between 0 and 0.27. In addition to the amount of seasonal precipitation, this ratio is influenced by the interannual distribution of rainfall events. Overall, an average annual groundwater recharge for the doline (sinkhole) catchment on As Sulb plateau of 5.1 mm has estimated for the simulation period. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

9.
Karst spring responses examined by process-based modeling   总被引：8，自引：0，他引：8  
Birk S  Liedl R  Sauter M 《Ground water》2006,44(6):832-836
Ground water in karst terrains is highly vulnerable to contamination due to the rapid transport of contaminants through the highly conductive conduit system. For contamination risk assessment purposes, information about hydraulic and geometric characteristics of the conduits and their hydraulic interaction with the fissured porous rock is an important prerequisite. The relationship between aquifer characteristics and short-term responses to recharge events of both spring discharge and physicochemical parameters of the discharged water was examined using a process-based flow and transport model. In the respective software, a pipe-network model, representing fast conduit flow, is coupled to MODFLOW, which simulates flow in the fissured porous rock. This hybrid flow model was extended to include modules simulating heat and reactive solute transport in conduits. The application of this modeling tool demonstrates that variations of physicochemical parameters, such as solute concentration and water temperature, depend to a large extent on the intensity and duration of recharge events and provide information about the structure and geometry of the conduit system as well as about the interaction between conduits and fissured porous rock. Moreover, the responses of solute concentration and temperature of spring discharge appear to reflect different processes, thus complementing each other in the aquifer characterization.  相似文献   

10.
Tracer breakthrough curves provide valuable information about the traced media, especially in inherently heterogeneous karst aquifers. In order to study the effect of variations in hydraulic gradient and conduit systems on breakthrough curves, a bench scale karst model was constructed. The bench scale karst model contains both matrix and a conduit. Eight tracing tests were conducted under a wide range of hydraulic gradients from 1 to greater than 5 for branchwork and network-conduit systems. Sampling points at varying distances from the injection point were utilized. Results demonstrate that mean tracer velocities, tracer mass recovery and linear rising slope of the breakthrough curves were directly controlled by hydraulic gradient. As hydraulic gradient increased, both one half the time for peak concentration and one fifth the time for peak concentration decreased. The results demonstrate the variations in one half the time for peak concentration and one fifth the time for peak concentration of the descending limb for different sampling points under differing hydraulic gradients are mainly controlled by the interactions of advection with dispersion. The results are discussed from three perspectives: different conduit systems, different hydraulic-gradient conditions, and different sampling points. The research confirmed the undeniable role of hydrogeological setting (i.e., hydraulic gradient and conduit system) on the shape of the breakthrough curve. The extracted parameters (mobile-fluid velocity, tracer-mass recovery, linear rising limb, one half the time for peak concentration, and one fifth the time for peak concentration) allow for differentiating hydrogeological settings and enhance interpretations the tracing tests in karst aquifers.  相似文献   

11.
Karst aquifers represent dual flow systems consisting of a highly conductive conduit system embedded in a less permeable rock matrix. Hybrid models iteratively coupling both flow systems generally consume much time, especially because of the nonlinearity of turbulent conduit flow. To reduce calculation times compared to those of existing approaches, a new iterative equation solver for the conduit system is developed based on an approximated Newton–Raphson expression and a Gauß–Seidel or successive over‐relaxation scheme with a single iteration step at the innermost level. It is implemented and tested in the research code CAVE but should be easily adaptable to similar models such as the Conduit Flow Process for MODFLOW‐2005. It substantially reduces the computational effort as demonstrated by steady‐state benchmark scenarios as well as by transient karst genesis simulations. Water balance errors are found to be acceptable in most of the test cases. However, the performance and accuracy may deteriorate under unfavorable conditions such as sudden, strong changes of the flow field at some stages of the karst genesis simulations.  相似文献   

12.
Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.  相似文献   

13.
A specific characteristic of karst systems is the occurrence of time variant recharge areas. In our study we present a new type of hydrological karst model and a new calibration approach both considering this specific characteristic. The new model type considers the spatial variability of karst system properties by distribution functions, and is compared to a simple reservoir model. Both models are applied to a karst system in Southern Spain where objective functions applied on hydrodynamic and hydrochemical information helped to determine model parameters playing a role for hydrodynamic response. Thereafter, the recharge area is determined separately for individual hydrological years and for the entire time series by calibrating the model to match the water balance. We show that hydrochemical information is crucial to find a reasonable set of parameters for both models. Considering different hydrological years, we find that the recharge area is changing significantly (from 28 to 53 km2). The newly developed model is able to reproduce this variation and provide acceptable simulation results for the entire time series of available data. The classic reservoir model shows inferior performance concerning hydrodynamics and fails to reproduce the water balance because it does not consider variations of recharge area. Our calibration approach allows identifying a variable recharge area and our new model is able to reproduce its variability. Hence we obtain a more realistic system representation, which can be of high significance when models are used for prediction, i.e. beyond the conditions they were calibrated, e.g. for land-use or climate change scenarios.  相似文献   

14.
Grasso DA  Jeannin PY 《Ground water》2002,40(6):608-617
Investigation techniques for karst flow systems are based mainly on the study of different signals leaving the system caused by natural or induced external influences. Each signal represents one of the systems outputs (e.g., hydraulic, chemical, physical, or isotopic responses) that reflect the characteristics of the entire system. In this paper, we present a method to infer information about the structure of karst systems. It is based on a simultaneous analysis of chemical and hydraulic responses. Beside the classical piston flow at the beginning of a flood pulse, we define a chemically based recession flow phase. During this phase, field data show that the concentration of total dissolved solids can be considered as an exponential function of the logarithm of flow. This relationship allows two parameters to be defined, one of which is dependent on the structure and degree of development of the karst conduit network, the other is dependent mainly on bioclimatic factors. Data collected from seven karst springs are used to support ideas introduced in the paper.  相似文献   

15.
South China karst aquifer storm-scale hydrochemistry   总被引：7，自引：0，他引：7  
Liu Z  Groves C  Yuan D  Meiman J 《Ground water》2004,42(4):491-499
  相似文献   

16.
Karst is regarded as a geological and physicochemical process of interaction between the natural water and soluble rock. Its development is possible only under non-equilibrium conditions in the “groundwater–congruently soluble rock” system. It is shown that intense karst development, leading to the formation of underground cavities and passages, is possible in the presence of advective transport of dissolved substance, which can redistribute its considerable amounts in a relatively short time. The specific features of karst development are considered in the case of its manifestation in carbonate, sulfate, and chloride rocks.  相似文献   

17.
Rapid changes in spring water quality in karst areas due to rapid recharge of bacterially contaminated water are a major concern for drinking water suppliers and users. The main objective of this study was to use field experiments with fecal indicators to verify the vulnerability of a karst spring to pathogens, as determined by using a numerical modeling approach. The groundwater modeling was based on linear storage models that can be used to simulate karst water flow. The vulnerability of the karst groundwater is estimated using such models to calculate criteria that influence the likelihood of spring water being affected by microbial contamination. Specifically, the temporal variation in the vulnerability, depending on rainfall events and overall recharge conditions, can be assessed and quantified using the dynamic vulnerability index (DVI). DVI corresponds to the ratio of conduit to diffuse flow contributions to spring discharge. To evaluate model performance with respect to predicted vulnerability, samples from a spring were analyzed for Escherichia coli, enterococci, Clostridium perfringens, and heterotrophic plate count bacteria during and after several rainfall events. DVI was shown to be an indication of the risk of fecal contamination of spring water with sufficient accuracy to be used in drinking water management. We conclude that numerical models are a useful tool for evaluating the vulnerability of karst systems to pathogens under varying recharge conditions  相似文献   

18.
Anomalous behaviour of specific electrical conductivity (SEC) was observed at a karst spring in Slovenia during 26 high‐flow events in an 18‐month monitoring period. A conceptual model explaining this anomalous SEC variability is presented and reproduced by numerical modelling, and the practical relevance for source protection zoning is discussed. After storm rainfall, discharge increases rapidly, which is typical for karst springs. SEC displays a first maximum during the rising limb of the spring hydrograph, followed by a minimum indicating the arrival of freshly infiltrated water, often confirmed by increased levels of total organic carbon (TOC). The anomalous behaviour starts after this SEC minimum, when SEC rises again and remains elevated during the entire high‐flow period, typically 20–40 µS/cm above the baseflow value. This is explained by variable catchment boundaries: When the water level in the aquifer rises, the catchment expands, incorporating zones of groundwater with higher SEC, caused by higher unsaturated zone thickness and subtle lithologic changes. This conceptual model has been checked by numerical investigations. A generalized finite‐difference model including high‐conductivity cells representing the conduit network (“discrete‐continuum approach”) was set up to simulate the observed behaviour of the karst system. The model reproduces the shifting groundwater divide and the nearly simultaneous increase of discharge and SEC during high‐flow periods. The observed behaviour is relevant for groundwater source protection zoning, which requires reliable delineation of catchment areas. Anomalous behaviour of SEC can point to variable catchment boundaries that can be checked by tracer tests during different hydrologic conditions. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

19.
This paper describes how climate influences the hydrology of an ephemeral depressional wetland. Surface water and groundwater elevation data were collected for 7 years in a Coastal Plain watershed in South Carolina USA containing depressional wetlands, known as Carolina bays. Rainfall and temperature data were compared with water‐table well and piezometer data in and around one wetland. Using these data a conceptual model was created that describes the hydrology of the system under wet, dry, and drought conditions. The data suggest this wetland operates as a focal point for groundwater recharge under most climate conditions. During years of below‐normal to normal rainfall the hydraulic gradient indicated the potential for groundwater recharge from the depression, whereas during years of above‐normal rainfall, the hydraulic gradient between the adjacent upland, the wetland margin, and the wetland centre showed the potential for groundwater discharge into the wetland. Using high‐resolution water‐level measurements, this groundwater discharge condition was found to hold true even during individual rainfall events, especially under wet antecedent soil conditions. The dynamic nature of the hydrology in this Carolina bay clearly indicates it is not an isolated system as previously believed, and our groundwater data expand upon previous hydrologic investigations at similar sites which do not account for the role of groundwater in estimating the water budget of such systems. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

20.
Discharge from karst springs contains a mixture of conduit and matrix water, but the variations in groundwater mixing are poorly known. Storm events present an opportunity to try to map flow components because water entering during storms is more dilute and provides a tracer as it mixes with pre‐event water along the flowpath from the recharge area to discharge at a spring. We used hysteresis plots of Mg/Ca ratios in a spring in the Cumberland Valley of Pennsylvania to map conduit (higher Ca) vs. diffuse (higher Mg) sources of recharge. We observed two types of temporal heterogeneity: within a storm event and from storm to storm. The timing of the variation in Mg/Ca suggested sources of mixing waters. An increase in the Mg/Ca ratio at the beginning of some storms while conductivity declined suggested diffuse recharge through the epikarst. The rapid changes in Mg/Ca ratios for low‐intensity events probably occurred as the rainfall waxed and waned and illustrate that a variety of flowpaths are available at this spring because additional flushing of Mg occurred. In contrast, the conductivity hysteresis began with dilute water initially and rotation was similar from storm to storm. Hysteresis plots of the Mg/Ca ratio have the potential of revealing more of the complexity in discharge than conductivity alone. A better understanding of flow components in karst is needed to protect these aquifers as a groundwater resource.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号