Hydrologic functioning of the deep critical zone and contributions to streamflow in a high‐elevation catchment: Testing of multiple conceptual models

High‐elevation mountain catchments are often subject to large climatic and topographic gradients. Therefore, high‐density hydrogeochemical observations are needed to understand water sources to streamflow and the temporal and spatial behaviour of flow paths. These sources and flow paths vary seasonally, which dictates short‐term storage and the flux of water in the critical zone (CZ) and affect long‐term CZ evolution. This study utilizes multiyear observations of chemical compositions and water residence times from the Santa Catalina Mountains Critical Zone Observatory, Tucson, Arizona to develop and evaluate competing conceptual models of seasonal streamflow generation. These models were tested using endmember mixing analysis, baseflow recession analysis, and tritium model “ages” of various catchment water sources. A conceptual model involving four endmembers (precipitation, soil water, shallow, and deep groundwater) provided the best match to observations. On average, precipitation contributes 39–69% (55 ± 16%), soil water contributes 25–56% (41 ± 16%), shallow groundwater contributes 1–5% (3 ± 2%), and deep groundwater contributes ~0–3% (1 ± 1%) towards annual streamflow. The mixing space comprised two principal planes formed by (a) precipitation‐soil water‐deep groundwater (dry and summer monsoon season samples) and (b) precipitation‐soil water‐shallow groundwater (winter season samples). Groundwater contribution was most important during the wet winter season. During periods of high dynamic groundwater storage and increased hydrologic connectivity (i.e., spring snowmelt), stream water was more geochemically heterogeneous, that is, geochemical heterogeneity of stream water is storage‐dependent. Endmember mixing analysis and ³H model age results indicate that only 1.4 ± 0.3% of the long‐term annual precipitation becomes deep CZ groundwater flux that influences long‐term deep CZ development through both intercatchment and intracatchment deep groundwater flows.     

Catchment‐scale sensitivity and uncertainty in water quality modelling

In this paper, we assess the performance of the catchment model SIMulated CATchment model (SIMCAT), to predict nitrate and soluble reactive phosphorus concentrations against four monitoring regimes with different spatial and temporal sampling frequencies. The Generalised Likelihood Uncertainty Estimation (GLUE) uncertainty framework is used, along with a general sensitivity analysis to understand relative parameter sensitivity. Improvements to model calibration are explored by introducing more detailed process representation using the Integrated Catchments model (INCA) water quality model, driven by the European hydrological predictions for the environment model. The results show how targeted sampling of headwater watercourses upstream of point discharges is essential for calibrating diffuse loads and can exert a strong influence on the whole‐catchment model performance. Further downstream, if the point discharges and loads are accurately represented, then the improvement in the catchment‐scale model performance is relatively small as more calibration points are added or frequency is increased. The higher‐order, dynamic model integrated catchments model of phosphorus dynamics, which incorporates sediment and biotic interaction, resulted in improved whole‐catchment performance over SIMCAT, although there are still large epistemic uncertainties from land‐phase export coefficients and runoff. However, the very large sampling errors in routine monitoring make it difficult to invest confidence in the modelling, especially because we know phosphorous transport to be very episodic and driven by high flow conditions for which there are few samples. The environmental modelling community seems to have been stuck in this position for some time, and whilst it is useful to use an uncertainty framework to highlight these issues, it has not widely been adopted, perhaps because there is no clear mechanism to allow uncertainties to influence investment decisions. This raises the question as to whether it might better place a cost on uncertainty and use this to drive more data collection or improved models, before making investment decisions concerning, for example, mitigation strategies. Copyright © 2016 John Wiley & Sons, Ltd.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

The effect of density-driven flow on the transport of high-concentration solutes in the hyporheic zone

In this study, both laboratory experiments and numerical simulations were conducted to investigate the effect of density-driven flow on the transport of high-concentration pollutants in the hyporheic zone. The results show that the density gradient can change the flow of pore water and the strong density-driven flow can lead to an unstable flow, which increases the effect of preferential flow and thus causes the appearance of solute fingers in the hyporheic zone. Notably, these solute fingers become more obvious with the increase of depth. The appearance of solute fingers depends on the relative strength of the pumping exchange and density gradient, which are represented by the dimensionless number M* and N* respectively. Finger flows appear near the interface when M^* is less than 0.5 N^*. This study may contribute to better understanding the transport and destination of solutes and thus may provide some insights into the assessment on pollution incidents.     

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.     

A novel platform to evaluate the dampening of water and solute transport in an experimental channel under unsteady flow conditions

This paper presents a novel platform to study the dampening of water and solute transport in an experimental channel under unsteady flow conditions, where literature data are scarce. We address the question about what could be the smallest size of experimental platform that is useful for research, project studies, and teaching activities and that allows to do rational experiments characterized by small space occupation, short experimental duration, high measurement precision, high quality and reproducible experimental curves, low water and energy consumption, and the possibility to test a large variety of hydrograph scenarios. Whereas large scale hydraulic laboratories have focused their studies on sediment transport, our platform deals with solute transport. The objectives of our study are (a) building a platform that allows to do rational experiments, (b) enriching the lack of experimental data concerning water and solute transport under unsteady state conditions, and (c) studying the dampening of water and solute transport. We studied solute transport in a channel with lateral gain and lateral loss under different experimental configurations, and we show how the same lateral loss flow event can lead to different lateral loss mass repartitions under different configurations. In order to characterize water and solute dampening between the input and the output of the channel, we calculate dampening ratios based on peak coordinates of time flow curves and time mass curves and that express the decrease of peak amplitude and the increase of peak occurrence time between the input and output curves. Finally, we use a solute transport model coupling the diffusive wave equation for water transfer and the advection–diffusion equation for solute transport in order to simulate the experimental data. The simulations are quite good with a Nash–Sutcliffe efficiency NSE > 0.98 for water transfer and 0.84 < NSE < 0.97 for solute transport. This platform could serve hydrological modellers because it offers a variety of measured parameters (flow, water height, and solute concentration), at a fine time step under unsteady flow conditions.     

Time‐lapse electrical resistivity imaging of solute transport in a karst conduit

The use of electrical resistivity surveys to locate karst conduits has shown mixed success. However, time‐lapse electrical resistivity imaging combined with salt injection improves conduit detection and can yield valuable insight into solute transport behaviour. We present a proof of concept above a known karst conduit in the Kentucky Horse Park (Lexington, Kentucky). A salt tracer solution was injected into a karst window over a 45‐min interval, and repeat resistivity surveys were collected every 20 min along a 125‐m transect near a monitoring well approximately 750 m downgradient from the injection site. In situ fluid conductivity measurements in the well peaked at approximately 25% of the initial value about 3 h after salt injection. Time‐lapse electrical resistivity inversions show two broad zones at the approximate conduit depth where resistivity decreased and then recovered in general agreement with in situ measurements. Combined salt injection and electrical resistivity imaging are a promising tool for locating karst conduits. The method is also useful for gaining insight into conduit geometry and could be expanded to include multiple electrical resistivity transects. Copyright © 2015 John Wiley & Sons, Ltd.     

Improvement of the solute transfer in a conceptual unsaturated zone scheme: a case study of the Seine River basin

For predicting the evolution of solute concentrations in groundwater and testing the impact of remediation policies, a coupling between the agronomical model STICS and the hydrogeological model MODCOU was implemented. When applied to the Seine River basin, this model accurately represents the temporal evolution of average nitrate concentrations in the aquifer, but with large local errors. We propose an improvement to the simple unsaturated zone (UZ) scheme NonsatSW used in STICS–MODCOU. The modifications are based on a comparison with the mechanistic model Metis considered as a reference as it solves Richards' equation. A more realistic saturation profile and a varying percolation rate are integrated in NonsatSW. This new model, named NonsatVG, is assessed by comparing it with NonsatSW and Metis. In an ideal case, NonsatVG generates a solute transfer and a dispersion closer to that of Metis than of NonsatSW. In real cases, without additional calibration, NonsatVG and Metis simulate better the average transfer velocities of the observed nitrate profiles. Furthermore, modifications in NonsatVG give a direct relationship between the depth of the water table and the saturation profile. We obtain, therefore, as in Metis, an evolution of the solute transfer velocity depending on the piezometric level. These dynamics are not simulated in NonsatSW. Despite a modified water transfer through the UZ, NonsatVG is also as valid as NonsatSW in the modelling of water transfer to the saturated zone. Finally, an application to the Seine basin shows that solute transfer velocities are lower with NonsatVG than with NonsatSW, but are in better agreement with literature. Copyright © 2010 John Wiley & Sons, Ltd.     

Experimental and computational validation and verification of the Stokes–Darcy and continuum pipe flow models for karst aquifers with dual porosity structure

In our previous study, we developed the Stokes–Darcy (SD) model was developed for flow in a karst aquifer with a conduit bedded in matrix, and the Beavers–Joseph (BJ) condition was used to describe the matrix–conduit interface. We also studied the mathematical well‐posedness of a coupled continuum pipe flow (CCPF) model as well as convergence rates of its finite element approximation. In this study, to compare the SD model with the CCPF model, we used numerical analyses to validate finite element discretisation methods for the two models. Using computational experiments, simulation codes implementing the finite element discretisations are then verified. Further model validation studies are based on the results of laboratory experiments. Comparing the results of computer simulations and experiments, we concluded that the SD model with the Beavers–Joseph interface condition is a valid model for conduit–matrix systems. On the other hand, the CCPF model with the value of the exchange parameter chosen within the range suggested in the literature perhaps does not result in good agreement with experimental observations. We then examined the sensitivity of the CCPF model with respect to the exchange parameter, concluding that, as has previously been noted, the model is highly sensitive for small values of the exchange parameter. However, for larger values, the model becomes less sensitive and, more important, also produces results that are in better agreement with experimental observations. This suggests that the CCPF model may also produce accurate simulation results, if one chooses larger values of the exchange parameter than those suggested in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Depth‐averaged,spatially distributed flow dynamic and solute transport modelling of a large‐scaled,subtropical constructed wetland

Constructed wetlands are being utilized worldwide to effectively reduce excess nutrients in agricultural runoff and wastewater. Despite their frequency, a multi‐dimensional, physically based, spatially distributed modelling approach has rarely been applied for flow and solute transport in treatment wetlands. This article presents a two‐dimensional hydrodynamic and solute transport modelling of a large‐scaled, subtropical, free water surface constructed wetland of about 8 km² in the Everglades of Florida, USA. In this study, MIKE 21 was adopted as the basic model framework. Field monitoring of the time series hydrological and chloride data, as well as spatially distributed data such as bathymetry and vegetation distribution, provided the necessary model input and testing data. Simulated water level profiles were in good agreement with the spatio‐temporal variations of measured ones. On average, the root‐mean‐square error of model calibration on annual water level fluctuations was 0·09 m. Manning's roughness coefficients for the dense emergent and submerged aquatic vegetation areas, which were estimated as a function of vegetation type, ranged from 0·67 to 1·0 and 0·12 to 0·15 s/m^1/3, respectively. The solute transport model calibration for four monitoring sites agreed well with the measured annual variations in chloride concentration with an average percent model error of about 15%. The longitudinal dispersivity was estimated to be about 2 m and was more than an order of magnitude higher than the transverse one. This study is expected to play the role of a stepping stone for future modelling efforts on the development and application of more advanced flow and transport models applicable to a variety of constructed wetland systems, as well as to the Everglades stormwater treatment areas in operation or in preparation. Copyright © 2010 John Wiley & Sons, Ltd.     

Mobility of major ions and nutrients in the unsaturated zone during paddy cultivation: a field study and solute transport modelling approach

Study of the movement of water and solute within soil profiles is important for a number of reasons. Accumulation of prominent contaminants from agricultural chemicals in the unsaturated zone over the years is a major concern in many parts of the world. As a result, the unsaturated zone has been a subject of great research interest during the past decade. Hence, an intensive field study was conducted in a part of Palar and Cheyyar river basins to understand the variation of major ions and nutrients in the soil zone during paddy cultivation. The chloride and nitrate data were used to model the movement of these chemicals in the unsaturated zone using the HYDRUS‐2D model. The field study shows that fertilizer application and irrigation return flow increases the major ions and nutrients concentration in the unsaturated zone. Further, the nutrient concentrations are regulated by plant uptake, fertilizer application and infiltration rate. Additionally, denitrification and soil mineralization processes also regulate the nitrogen concentration in the unsaturated zone. The solute transport modelling study concluded that the simulated results match reasonably with the observed trends. Simulated concentrations of chloride and nitrate for a 5‐year period indicate that the concentrations of these ions fluctuate in a cyclic manner (from 60 to 68 mg l^?1 and from 3·4 to 3·5 mg l^?1 respectively in groundwater) with no upward and downward trend. The influence of excessive fertilizer application on groundwater was also modelled. The model predicts an increase of about 17 mg l^?1 of chloride and 2·3 mg l^?1 of nitrogen in the groundwater of this area when the application of fertilizers is doubled. The model indicates that the present level of use of agrochemicals is no threat to the groundwater quality. Copyright © 2007 John Wiley & Sons, Ltd.     

A simulation of large‐scale groundwater flow and travel time in a fractured rock environment for waste disposal purposes

Two‐dimensional regional groundwater flow was simulated based on a conceptual model of low‐permeability crystalline rocks of the Whiteshell Research Area (WRA) in south‐eastern Manitoba. The conceptual model consists of fracture zones that strike in different directions and dip at various angles in the background rock mass. The thickness and hydraulic properties of the fracture zones in the conceptual model were varied as were the fluid properties and the boundary conditions of the groundwater flow system. The effects of these variations on the groundwater flow pattern and on the convective travel time along pathways from a hypothetical disposal vault at 500 m depth to discharge locations at the ground surface were evaluated. The vault was located in the regional discharge area of the groundwater system. A homogeneous conceptual model of the WRA, having only freshwater flow, formed a groundwater flow pattern with a regional flow system. Local flow systems developed increasingly with the introduction of fracture zones 20 m and 3 m thick, and depth‐dependent fluid density. This indicates a reduction in groundwater residence time by fracture zones and fluid density. Flow pathways were analysed using both a stream‐function and a particle‐tracking technique. The pathways and their lengths from the location of the vault to the surface varied spatially according to the flow patterns. The minimum travel time along these pathways was less than 150 000 and greater than 4 000 000 years in models with and without fracture zones, respectively, indicating that the presence of fracture zones was the major controlling factor. A precise knowledge and refinement of conceptual model parameters is necessary during site selection for waste disposal purposes. Copyright © 2004 John Wiley & Sons, Ltd.     

Water and air dynamics within a deep vadose zone of a karst massif: Observations from the Lukina jama–Trojama cave system (−1,431 m) in Dinaric karst (Croatia)

Extreme heterogeneity of karst systems makes them very challenging to study. Various processes within the system affect its global response, usually measured at karst springs. Research conducted in caves provides a unique opportunity for in situ analysis of separate processes in karst underground. The aim of the present study was to research the water and air dynamics within a deep karst system. Air and water basic physical parameters across the Lukina jama–Trojama cave system (?1,431 m) were continuously monitored during a 1‐year period. Recorded hydrograph of the siphon lake at the bottom of the cave was used to interpret the characteristics of an unexplored phreatic/epiphreatic conduit network. Water origin in the siphon was determined based on temperature and electrical conductivity. Air temperature and humidity monitoring revealed a strong inflow of air of sub‐zero temperature into the upper portion of the cave during winter. Cave passage morphology was interpreted as the main determinant of air dynamics, which caused ice to accumulate extensively in the upper portions of the cave and caused the temperature on the top of the homothermic zone to be significantly below the mean outside temperature. Air dynamics also lowered the temperature of water flowing through the cave vadose zone and feeding the phreatic zone of the massif. The pronounced temperature difference between the phreatic zone and the top of the homothermic zone probably contributed to the thermal gradient observed in the cave, which is steeper than in ice‐free caves in the area. Our results enabled the development of a conceptual model that describes coupling between air and water dynamics in the cave system and its surroundings.     

Water flow and solute transport model of potentially toxic elements through unsaturated zone at regional wellfield Kosnica

The objective of this work was to build a prognostic water flow model and potentially toxic elements (lead, cadmium, zinc) transport model in the unsaturated zone. Research was conducted in the catchment area of Kosnica regional wellfield, where the unsaturated zone is characterised by Fluvisol. Lower sorption capacities were determined in the first horizons for all three potentially toxic elements. Correlation coefficient of the measured and simulated values of tracer concentration is 0.58 for the AC horizon and 0.84 for the 2C/C1 horizon. Based on calibrated water flow and transport parameters, a prognostic water flow model and potentially toxic elements (lead, cadmium, zinc) transport model in the unsaturated zone was built. In case of an accidental spill of potentially toxic elements with concentrations of 1000 mg/l, the risk of contamination of the aquifer is present. Copyright © 2016 John Wiley & Sons, Ltd.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 2. Controls on weathering and permeability profiles,drainage hydraulics,and solute export pathways

The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 1. Hydraulic theory of reactive transport

This is the first of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory (also known as Dupuit or Boussinesq theory). Here, we examine the effect of lateral flow on the downward fluxes of water and solutes through perched groundwater at steady state. We derive analytical expressions describing the decline in the downward flux rate with depth. Using these, we obtain analytical expressions for water age in a number of cases. The results show that when the permeability field is homogeneous, the spatial structure of water age depends qualitatively on a single dimensionless number, Hi. This number captures the relative contributions to the lateral hydraulic potential gradient of the relief of the lower‐most impermeable boundary (which may be below the weathering front within permeable or incipiently weathered bedrock) and the water table. A “scaled lateral symmetry” exists when Hi is low: age varies primarily in the vertical dimension, and variations in the horizontal dimension x almost disappear when the vertical dimension z is expressed as a fraction z/H(x) of the laterally flowing system thickness H(x). Taking advantage of this symmetry, we show how the lateral dimension of the advection–diffusion‐reaction equation can be collapsed, yielding a 1‐D vertical equation in which the advective flux downward declines with depth. The equation holds even when the permeability field is not homogeneous, as long as the variations in permeability have the same scaled lateral symmetry structure. This new 1‐D approximation is used in the accompanying paper to extend chemical weathering models derived for 1‐D columns to hillslope domains.     

Determination of the conduit structure in a karst aquifer based on tracer data—Lurbach system,Austria

A structure model was used to analyse solute‐transport parameter estimates based on tracer breakthrough curves. In the model system, groundwater flow is envisioned to be organised in a complex conduit network providing a variety of short circuits with relative small carrying capacities along different erosion levels. The discharge through the fully filled conduits is limited owing to void geometries and turbulent flow; thus, a hierarchic overflow system evolves where conduits are (re‐)activated or dried up depending on the flow condition. Exemplified on the Lurbach–Tanneben karst aquifer, the applicability of the model approach was tested. Information derived from multi‐tracer experiments performed at different volumetric flow rates enabled to develop a structural model of the karst network, under constraint of the geomorphological and hydrological evolution of the site. Depending on the flow rate, groundwater is divided into up to eight flow paths. The spatial hierarchy of flow paths controls the sequence of flow path activation. Conduits of the topmost level are strongly influenced by reversible alteration processes. Sedimentation or blocking causes an overflow of water to the next higher conduit. Flow path specific dissolutional denudation rates were estimated using the temporal development of the partial discharge. Copyright © 2012 John Wiley & Sons, Ltd.