Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

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 δ¹⁸O) 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.     

Characterizing the variability of transit time distributions and young water fractions in karst catchments using flux tracking

Hydrological and biogeochemical processes in karst environments are strongly controlled by heterogeneous fracture-conduit networks. Quantifying the spatio-temporal variability of water transit time and young water fractions in such heterogeneous hydrogeological systems is fundamental to linking discharge and water quality dynamics in the karst critical zone. We used a tracer-aided conceptual hydrological model to track the fate of each hour of rain input individually. Using this approach, the variability of transit time distributions and young water fraction were estimated in the main landscape units in a karst catchment of Chenqi in Guizhou Province, Southwest China. The model predicted that the mean young water (i.e., <~2 months old) fraction of ground conduit flow is 0.31. Marked seasonal variabilities in water storage and hydrological connectivity between the conduit network and fractured matrix, as well as between hillslopes and topographic depression, drive the dynamics of young water fraction and travel time distributions in each landscape unit. Especially, the strong hydrological connectivity between the land surface and underground conduits caused by the direct infiltration through large fractures and sinkholes, leads the drastic increasement in young water fraction of runoff after heavy rain. Even though the contribution of young water to runoff is greater, the strong mixing and drainage of small fractures accelerate the old water release during high flows during the wet season. It is notable that the young water may sometimes be the most contaminated component contributing to the underground conduit network in karst catchments, because of the direct transfer of contaminants from the ground surface with rain water via large fractures and sinkholes.     

Modelling hydrological processes influenced by soil,rock and vegetation in a small karst basin of southwest China

Hydrological processes in karst basins are controlled by permeable multimedia, consisting of soil pores, epikarst fractures, and underground conduits. Distributed modelling of hydrological dynamics in such heterogeneous hydrogeological conditions is a challenging task. Basing on the multilayer structure of the distributed hydrology‐soil‐vegetation model (DHSVM), a distributed hydrological model for a karst basin was developed by integrating mathematical routings of porous Darcy flow, fissure flow and underground channel flow. Specifically, infiltration and saturated flow movement within epikarst fractures are expressed by the ‘cubic law’ equation which is associated with fractural width, direction, and spacing. A small karst basin located in Guizhou province of southwest China was selected for this hydrological simulation. The model parameters were determined on the basis of field measurement and calibrated against the observed soil moisture contents, vegetation interception, surface runoff, and underground flow discharges from the basin outlet. The results show that due to high permeability of the epikarst zone, a significant amount of surface runoff is only generated after heavy rainfall events during the wet season. Rock exposure and the epikarst zone significantly increase flood discharge and decrease evapotranspiration (ET) loss; the peak flood discharge is directly proportional to the size of the aperture. Distribution of soil moisture content (SMC) primarily depends on topographic variations just after a heavy rainfall, while SMC and actual ET are dominated by land cover after a period of consecutive non‐rainfall days. The new model was able to capture the sharp increase and decrease of the underground streamflow hydrograph, and as such can be used to investigate hydrological effects in such rock features and land covers. Copyright © 2011 John Wiley & Sons, Ltd.     

Role of epikarst in near‐surface hydrological processes in a soil mantled subtropical dolomite karst slope: implications of field rainfall simulation experiments

Epikarst exerts a strong control on run‐off generation in karst regions, but it is still unclear in karst regions. Our study aimed to demonstrate the effect of epikarst on near‐surface hydrological processes in a subtropical cockpit karst region of southwest China, using plot‐scale rainfall simulation experiments with different rainfall intensities (low and high) and antecedent moisture conditions (dry and wet). A trench excavated to the epikarst lower boundary allowed identification of flow pathways in the entire soil–epikarst architecture system, thus facilitating the water balance calculations using a conceptual model with the assumption of a two‐stage hydrological evolution. More than 70% of the total rainfall water moved vertically through the shallow soil layer and then was redistributed by the epikarst as subsurface flow occurring on the soil–epikarst interface, depression filling on epikarst surface, water held by epikarst and deep percolation. Epikarst water regulation capacity, defined as the sum of depression filling on epikarst surface, water held by epikarst, epikarst seepage flow and deep percolation, was 58 mm (wet antecedent condition) and 223 mm (dry antecedent condition). Total run‐off from the soil–epikarst system was dominated by saturated subsurface flow showing a threshold process controlled by epikarst storage capacity (storing as much as 181 mm of rainfall water under dry antecedent condition). Our study proved that despite the epikarst being relatively poorly developed and covered by a soil mantle, it still exerted a strong influence on near‐surface hydrological processes and thus should be adequately considered in future modelling of water recharge and depletion dynamics in this integrated soil–epikarst system. Copyright © 2015 John Wiley & Sons, Ltd.     

Mass transport triggered by heavy rainfall: the role of endorheic basins and epikarst in a regional karst aquifer

This study aims at recognizing the mechanisms of mass transport between the karst surface and the saturated zone in a morphostructural relief of the Mesozoic karst carbonate platform of Murgia (Puglia, Southern Italy). The large dimension of the karst aquifer, the regional scale of the flow system, the boundary condition constituted by the sea, and the lack of freshwater springs constrain to the use of wells as monitoring points and limit the study area to the recharge area comprising 986 endorheic basins. The concentrations of non‐reactive tracers (nitrates) in the waters of autogenic recharge (from endorheic basins) have been modeled through the evaluation of effective infiltration, land use, and nitrogen surplus, with reference to a time window, which includes a low precipitation period followed by significant rainfall events. The comparison between the modeled nitrate concentrations and the nitrate concentrations measured in ground waters, coupled with the analysis of groundwater chemograms and records of hydraulic heads (all referred to the same time window), allows inferring the mechanism of mass transport between the karst surface and the groundwater table. The mass transport conceptual model requires the presence of the epikarst. The infiltration of significant rainfall in the endorheic basins after a low precipitation period displaces waters stored in the epikarst toward the saturated zone. Ground waters in the post‐event period show higher concentrations of nitrates, lower concentrations of total organic carbon, and higher Mg/Ca ratios than both those of the pre‐event period and the autumn‐winter recharge period. The post‐event recharge from epikarst storage determines a transient hazard of groundwater pollution with a time lag from the occurrence of the heavy rainfall.     

Anomalous behaviour of specific electrical conductivity at a karst spring induced by variable catchment boundaries: the case of the Podstenjšek spring,Slovenia

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.     

Hydrochemical variations during flood pulses in the south‐west China peak cluster karst: impacts of CaCO3–H2O–CO2 interactions

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 Ca²⁺ and HCO₃^? were the dominant (>90%) ions, we developed linear relationships (both r² > 0·91) between conductivity and those ions, respectively, and in turn calculated variations in the calcite saturation index (SI_C) and CO₂ 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 SI_C is lower. Simultaneously, P of conduit water during the flood period is lower than that at lower flows, and its SI_C 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 CO₂ 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 CO₂ 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.     

Contributing sources to baseflow in pre‐alpine headwaters using spatial snapshot sampling

Mountainous headwaters consist of different landscape units including forests, meadows and wetlands. In these headwaters it is unclear which landscape units contribute what percentage to baseflow. In this study, we analysed spatiotemporal differences in baseflow isotope and hydrochemistry to identify catchment‐scale runoff contribution. Three baseflow snapshot sampling campaigns were performed in the Swiss pre‐alpine headwater catchment of the Zwäckentobel (4.25 km²) and six of its adjacent subcatchments. The spatial and temporal variability of δ²H, Ca, DOC, AT, pH, SO₄, Mg and H₄SiO₄ of streamflow, groundwater and spring water samples was analysed and related to catchment area and wetland percentage using bivariate and multivariate methods. Our study found that in the six subcatchments, with variable arrangements of landscape units, the inter‐ and intra catchment variability of isotopic and hydrochemical compositions was small and generally not significant. Stream samples were distinctly different from shallow groundwater. An upper spring zone located near the water divide above 1,400 m and a larger wetland were identified by their distinct spatial isotopic and hydrochemical composition. The upstream wetland percentage was not correlated to the hydrochemical streamflow composition, suggesting that wetlands were less connected and act as passive features with a negligible contribution to baseflow runoff. The isotopic and hydrochemical composition of baseflow changed slightly from the upper spring zone towards the subcatchment outlets and corresponded to the signature of deep groundwater. Our results confirm the need and benefits of spatially distributed snapshot sampling to derive process understanding of heterogeneous headwaters during baseflow. Copyright © 2015 John Wiley & Sons, Ltd.     

Using Stormwater Hysteresis to Characterize Karst Spring Discharge

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.     

Defining hydrochemical evolution of streamflow through flowpath dynamics in Kawakami headwater catchment,Central Japan

The hydrochemical behaviour of catchments is often investigated by inferring stream chemistry through identification of source areas involved in hydrograph separation analysis, yet its dynamic evolution of hydrologic pathways has received little attention. Intensive hydrometric and hydrochemical measurements were performed during two different storms on March 29, 2001 and August 21–22, 2001 to define hydrochemical evolution under the dynamic of flow pathways in a 5·2 ha first‐order drainage of the Kawakami experimental basin (KEB), Central Japan, a forested headwater catchment with various soil depths (1·8 to 5 m) overlying late Neogene of volcanic bedrocks. The hydraulic potential distribution and flow lines data showed that the change in flow direction, which was controlled by rainfall amount and antecedent wetness of the soil profile, agreed well with the hydrochemical change across the slope segment during the storm. Hydrograph separation predicted by end‐member mixing analysis (EMMA) using Ca²⁺ and SiO₂ showed that near surface riparian, hillslope soil water and deep riparian groundwater were important in stream flow generation. The evidence of decrease in solutes concentration at a depth of 1 m in the hillslope and 0·6 m in the near surface riparian during peak storm suggested a flushing of high solutes concentration. Most of the solutes accumulated in the deep riparian groundwater zone, which was due to prominent downward flow and agreed well with the residence time. The distinct flow pathways and chemistry between the near surface riparian and deep riparian groundwater zones and the linkage hillslope aquifer and near surface riparian reservoir, which controls rapid flow and solutes flushing during the storm event, are in conflict with the typical assumption that the whole riparian zone resets flow pathways and chemical signature of hillslope soil water, as has been reported in a previous study. Copyright © 2005 John Wiley & Sons, Ltd.     

Bridging the gap between numerical solutions of travel time distributions and analytical storage selection functions

Recent advancements in analytical solutions to quantify water and solute travel time distributions (TTDs) and the related StorAge Selection (SAS) functions synthesize catchment complexity into a simplified, lumped representation. Although these analytical approaches are efficient in application, they require rarely available long‐term and high‐frequency hydrochemical data for parameter estimation. Alternatively, integrated hydrologic models coupled to Lagrangian particle‐tracking approaches can directly simulate age under different catchment geometries and complexity, but at a greater computational expense. Here, we bridge the two approaches, using a physically based model to explore the uncertainty in the estimation of the SAS function shape. In particular, we study the influence of subsurface heterogeneity, interactions between distinct flow domains (i.e., the vadose zone and saturated groundwater), diversity of flow pathways, and recharge rate on the shape of TTDs and the SAS functions. We use an integrated hydrology model, ParFlow, linked with a particle‐tracking model, SLIM, to compute transient residence times (or ages) at every cell in the domain, facilitating a direct characterization of the SAS function. Steady‐state results reveal that the SAS function shape shows a wide range of variation with respect to the variability in the structure of subsurface heterogeneity. Ensembles of spatially correlated realizations of hydraulic conductivity indicate that the SAS functions in the saturated groundwater have an overall weak tendency toward sampling younger ages, whereas the vadose zone gives a strong preference for older ages. We further show that the influence of recharge rate on the TTD is tightly dependent on the variability of subsurface hydraulic conductivity.     

HYDROCHEMICAL DECIPHERING OF STREAMFLOW GENERATION IN SEMI-ARID EAST AFRICA

There is a critical lack of knowledge regarding the dynamics of streamflow generation in the semi-arid tropics, particularly in Africa. In this project runoff mechanisms in forested and non-forested degraded catchments in northern Tanzania were studied using combined hydrometrical and hydrochemical methods. Following the hydrochemical identification of several flowpaths contributing to runoff, hydrograph separation by an end-member model based on K and Ca was undertaken. Results from the forested catchment indicate that stormflow was dominated by event water (about 75%), via overland flow and throughflow. The proportion of pre-event water (groundwater) displaced into the stream by a suggested riparian groundwater ridge mechanism varied, depending on the rainfall characteristics. In the non-forested, degraded catchment, nearly all stormflow was event water, and groundwater discharge was unaffected by rainfall. It is suggested that macropore flow is pivotal to the transmission of rainfall to runoff via throughflow, particularly in semi-arid tropical areas.     

Nitrate fluctuations at the water table: implications for recharge processes and solute transport in the Chalk aquifer

This study investigates fluctuations in nitrate concentration at the water table to improve understanding of unsaturated zone processes in the Chalk aquifer. Sampling was conducted using a novel multi‐level sampler during periods of water table rise over 5 years at a vertical resolution of 0.05 m. Nitrate concentration increased as the water table seasonally recovered, with similar inter‐annual trends with depth. The rising water table activated horizontal fractures facilitating the delivery of water elevated by up to 10 mg/l of nitrate with respect to the adjacent groundwater below. These fractures are considered to activate via piston displacement of water from the adjoining matrix. Hydrograph analysis identified 16 events which significantly perturbed the water table within 24–48 h of rainfall. Consistent nitrate concentrations indicate recharge through persistent fracture flow from the surface was not generally the primary driver of the rapid water table response during these events. Instead, the response was attributed to the piston displacement of porewater immediately above the water table. However, a single event in November 2012 delivered relatively dilute recharge indicating rapid persistent fracture flow following rainfall was possible to a depth of 14–15 m. Decreases in porewater nitrate concentration around fracture horizons and the dilution of many groundwater samples with respect to porewaters indicate a fresher source of water at depth. This was considered most likely to be a result of near surface water bypassing the matrix because of widespread mineralization on fracture surfaces, which retard water and solute exchange. Therefore, persistent fracture flow maybe considered a frequent process, operating independently of the matrix, and is not necessarily event driven. Copyright © 2015 John Wiley & Sons, Ltd.     

基于修正SWAT模型的岩溶地区非点源污染模拟初探——以横港河流域为例

岩溶流域含水系统的主要特征之一是连通地表的落水洞等垂直管道将近水平的地下暗河联系起来,降水及其形成的地表径流可以通过这些管道迅速地灌入地下河系,从而改变了水及其所携带的非点源污染物质在垂直与水平方向的传输速度与数量,使岩溶流域内地表-地下之间的物质交换与传输过程变得比较复杂;应用广泛的SWAT模型在模拟岩溶地区的水文、水质时会存在一些不足与局限.为此,本文针对岩溶水系统特征,引入落水洞、伏流、暗河的水文过程以及主要营养盐的输移过程,修正SWAT模型原有的水文循环过程及相关算法,改变其只适用于松散均匀介质流域非点源污染模拟的单一特征,并研究建立适合于岩溶流域的非点源污染模型和相应的模拟方法.选取横港河流域岩溶地区作为非点源污染的对象,应用修正后的模型通过控制性的模拟方法和敏感性性分析,定量评估落水洞、伏流、暗河等岩溶特征对氮、磷等主要非点源污染物质输移的影响及其带来的时空效应,并进一步探讨落水洞、伏流、暗河等对地表-地下水文与营养盐的交互作用及转换机理.结果表明,岩溶特征对流域的氮、磷负荷有增加作用,其中总磷的增加明显大于总氮的增加,总磷和总氮的增量分别为0.86%和2.12%;植被岩溶指数的增加会导致流域可溶性磷、有机磷的产出量增加,有机氮、地表产流中硝酸氮和沉积磷的产出量则居其次,落水洞改变了降雨的产流方式,增加了落水洞所在流域的有机磷和有机氮的产出,其增量变化在0~0.7和0~0.3 kg/hm²之间.     

旅游酒店排污影响下的岩溶地下水水化学变化

岩溶地下水系统具有高度的开放性和脆弱性,使得地下水极易遭受污染.为探究旅游酒店排污对岩溶地下水水化学变化的影响,以金佛山世界自然遗产地水房泉流域为例,对水房泉地下水的水化学进行自动化监测,对流域内雨水、土壤水、某酒店自来水、污废水进行定期采集,结合流域硝酸盐氮氧同位素分析.结果表明,监测期间水房泉水化学变化随酒店入住游客量总体表现为3个阶段：前期受降雨影响显著,旅游高峰期间水质急剧恶化,后期水质明显好转.酒店生活污废水的排放成为水房泉水化学演变的重要因素,硝酸盐氮氧同位素表明水房泉的NO₃^-主要来自粪便、污废水以及土壤N的混合.H₂SO₄及污水中HCl、有机酸等可能参与了碳酸盐岩的溶蚀,使水房泉Ca²⁺、Mg²⁺、HCO₃^-浓度增加明显.流域岩溶管道发育,地下水流速快,使污染物质扩散迅速,故在研究期间水房泉主要离子的浓度高峰对污废水排放高峰的响应仅滞后约4 d.     

Application of a Discrete-Continuum Model to Karst Aquifers in North China

A generalized discrete-continuum model is developed to simulate ground water flow in the karst aquifers of North China. The model is a hybrid numerical flow model, which takes into account both quick conduit flow and diffusive fissure flow. The conduit flow is represented by a discrete network model, and the fissure flow is modeled by a continuum approach. The developed model strongly emphasizes the function of the conduits in the flow fields. They control the general drainage pattern, as demonstrated in the simulation of a complex karst aquifer in North China. The model reproduces reasonably well the flow field in response to an unanticipated discharge of ground water from the karst aquifer into an underground mine based on the aquifer parameters that are manually calibrated from a multiple-well pumping test. Sensitivity of the model to the aquifer parameters was evaluated in the context of the case study.     

Estimating groundwater recharge for an arid karst system using a combined approach of time‐lapse camera monitoring and water balance modelling

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 · 10⁴ m²) 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.     

Conceptual models for carbonate aquifers. 1969

The very diverse types of ground‐water behavior in carbonate terrains can be classified by relating the flow type to a particular hydrogeologic environment each exhibiting a characteristic cave morphology. The ground water may move by diffuse flow, by retarded flow, or by free flow. Diffuse flow occurs in less soluble rocks such as extremely shaley limestones or crystalline dolomites. Integrated conduits are rare. Caves tend to be small, irregular, and often little more than solutionally widened joints. Retarded flows occur in artesian environments and in situations where unfavorable stratigraphy forces ground water to be confined to relatively thin beds. Network cave patterns are characteristic since hydrodynamic forces are damped by the external controls. Solution occurs along many available joints. Free flowing aquifers are those in which solution has developed a subsurface drainage system logically regarded as an underground extension of surface streams. These streams may have fully developed surface tributaries as well as recharge from sinkholes and general infiltration. Characteristic cave patterns are those of integrated conduit systems which are often truncated into linear, angulate, and branchwork caves. Free Flow aquifers may be further subdivided into Open aquifers lying beneath karst plains and Capped aquifers in which significant parts of the drainage net lie beneath an insoluble cap rock. Other geologic factors such as structure, detailed lithology, relief, and locations of major streams, control the details of cave morphology and orientation of the drainage network.     

Assessment of groundwater recharge processes through karst vadose zone by cave percolation monitoring

Recharge processes of karst aquifers are difficult to assess given their strong heterogeneity and the poorly known effect of vadose zone on infiltration. However, recharge assessment is crucial for the evaluation of groundwater resources. Moreover, the vulnerability of karst aquifers depends on vadose zone behaviour because it is the place where most contamination takes place. In this work, an in situ experimental approach was performed to identify and quantify flow and storage processes occurring in karst vadose zone. Cave percolation monitoring and dye tracing were used to investigate unsaturated zone hydrological processes. Two flow components (diffuse and quick) were identified and, respectively, account for 66% and 34% of the recharge. Quickflow was found to be the result of bypass phenomenon in vadose zone related to water saturation. We identify the role of epikarst as a shunting area, most of the storage in the vadose zone occurring via the diffuse flow component in low permeability zones. Relationship between rainfall intensity and transit velocity was demonstrated, with 5 times higher velocities for the quick recharge mode than the diffuse mode. Modelling approach with KarstMod software allowed to simulate the hybrid recharge through vadose zone and shows promising chances to properly assess the recharge processes in karst aquifer based on simple physical models.     

Change in low flows due to catchment management dynamics—Application of a comparative modelling approach

Understanding the natural low flow of a catchment is critical for effective water management policy in semi-arid and arid lands. The Geba catchment in Ethiopia, forming the headwaters of Tekeze-Atbara basin was known for its severe land degradation before the recent large scale Soil and Water conservation (SWC) programs. Such interventions can modify the hydrological processes by changing the partitioning of the incoming rainfall on the land surface. However, the literature lacks studies to quantify the hydrological impacts of these interventions in the semi-arid catchments of the Nile basin. Statistical test and Indicators of Hydrological Alteration (IHA) were used to identify the trends of streamflow in two comparatives adjacent (one treated with intensive SWC intervention and control with fewer interventions) catchments. A distributed hydrological model was developed to understand the differences in hydrological processes of the two catchments. The statistical and IHA tools showed that the low flow in the treated catchment has significantly increased while considerably decreased in the control catchment. Comparative analysis confirmed that the low flow in the catchment with intensive SWC works was greater than that of the control by >30% while the direct runoff was lower by >120%. This implies a large proportion of the rainfall in the treated catchment is infiltrated and recharge aquifers which subsequently contribute to streamflow during the dry season. The proportion of soil storage was more than double compared to the control catchment. Moreover, hydrological response comparison from pre- and post-intervention showed that a drastic reduction in direct runoff (>84%) has improved the low flow by >55%. This strongly suggests that the ongoing intensive SWC works have significantly improved the low flows while it contributed to the reduction of total streamflow in the catchment.