A quantitative method for the characterisation of karst aquifers based on spring hydrograph analysis, 2005, Kovacs A. , Perrochet P. , Kiraly L. , Jeannin P. Y. ,

This paper presents a method for characterizing flow systems in karst aquifers by acquiring quantitative information about the geometric and hydraulic aquifer parameters from spring hydrograph analysis. Numerical sensitivity analyses identified two fundamentally different flow domains, depending on the overall configuration of aquifer parameters. These two domains have been quantitatively characterized by deducing analytical solutions for the global hydraulic response of simple two-dimensional model geometries. During the baseflow recession of mature karst systems, the hydraulic parameters of karst conduits do not influence the drainage of the low-permeability matrix. In this case the drainage process is influenced by the size and hydraulic parameters of the low-permeability blocks alone. This flow condition has been defined as matrix-restrained flow regime (MRFR). During the baseflow recession of early karst systems and fissured systems, as well as the flood recession of mature systems, the recession process depends on the hydraulic parameters and the size of the low-permeability blocks, conduit conductivity and the total extent of the aquifer. This flow condition has been defined as conduit-influenced flow regime (CIFR). Analytical formulae demonstrated the limitations of equivalent models. While equivalent discrete-continuum models of early karst systems may reflect their real hydraulic response, there is only one adequate parameter configuration for mature systems that yields appropriate recession coefficient. Consequently, equivalent discrete-continuum models are inadequate for simulating global response of mature karst systems. The recession coefficient of equivalent porous medium models corresponds to the transition between matrix-restrained and conduit-influenced flow. Consequently, equivalent porous medium models yield corrupted hydrographs both in mature and early systems, and this approach is basically inadequate for modelling global response of karst aquifers. (c) 2004 Elsevier B.V. All rights reserved

A dimensionless number describing the effects of recharge and geometry on discharge from simple karstic aquifers, 2009, Covington M. D. , Wicks C. M. , Saar M. O.

The responses of karstic aquifers to storms are often used to obtain information about aquifer geometry. In general, spring hydrographs are a function of both system geometry and recharge. However, the majority of prior work on storm pulses through karst has not studied the effect of recharge on spring hydrographs. To examine the relative importance of geometry and recharge, we break karstic aquifers into elements according to the manner of their response to transient flow and demonstrate that each element has a characteristic response timescale. These fundamental elements are full pipes, open channels, reservoir/constrictions, and the porous matrix. Taking the ratio of the element timescale with the recharge timescale produces a dimensionless number, γ, that is used to characterize aquifer response to a storm event. Using sets of simulations run with randomly selected element parameters, we demonstrate that each element type has a critical value of γ below which the shape of the spring hydrograph is dominated by the shape of the recharge hydrograph and above which the spring hydrograph is significantly modified by the system geometry. This allows separation of particular element/storm pairs into recharge-dominated and geometry-dominated regimes. While most real karstic aquifers are complex combinations of these elements, we draw examples from several karst systems that can be represented by single elements. These examples demonstrate that for real karstic aquifers full pipe and open channel elements are generally in the recharge-dominated regime, whereas reservoir/constriction elements can fall in either the recharge- or geometry-dominated regimes.

Characterisation and modelling of conduit restricted karst aquifers – Example of the Auja spring, Jordan Valley, 2014, Schmidta Sebastian, Geyera Tobias, Guttmanb Joseph, Mareic Amer, Riesd Fabian, Sauter Martin

The conduit system of mature karstified carbonate aquifers is typically characterised by a high hydraulic conductivity and does not impose a major flow constriction on catchment discharge. As a result, discharge at karst springs is usually flashy and displays pronounced peaks following recharge events. In contrast, some karst springs reported in literature display a discharge maximum, attributed to reaching the finite discharge capacity of the conduit system (flow threshold). This phenomenon also often leads to a non-standard recession behaviour, a so called “convex recession”, i.e. an increase in the recession coefficient during flow recession, which in turn might be used as an indicator for conduit restricted aquifers. The main objective of the study is the characterisation and modelling of those hydrogeologically challenging aquifers. The applied approach consists of a combination of hydrometric monitoring, a spring hydrograph recession and event analysis, as well as the setup and calibration of a non-linear reservoir model. It is demonstrated for the Auja spring, the largest freshwater spring in the Lower Jordan Valley. The semi-arid environment with its short but intensive precipitation events and an extended dry season leads to sharp input signals and undisturbed recession periods. The spring displays complex recession behaviour, exhibiting exponential (coefficient α) and linear (coefficient β) recession periods. Numerous different recession coefficients α were observed: ∼0.2 to 0.8 d−1 (presumably main conduit system), 0.004 d−1 (fractured matrix), 0.0009 d−1 (plateau caused by flow threshold being exceeded), plus many intermediate values. The reasons for this observed behaviour are the outflow threshold at 0.47 m3 s−1 and a variable conduit–matrix cross-flow in the aquifer. Despite system complexity, and hence the necessity of incorporating features such as a flow threshold, conduit–matrix cross-flow, and a spatially variable soil/epikarst field capacity, the developed reservoir model is regarded as relatively simplistic. As a number of required parameters were calculated from the hydrogeological analysis of the system, it requires only six calibration parameters and performs well for the highly variable flow conditions observed. Calculated groundwater recharge in this semi-arid environment displays high interannual variability. For example, during the 45-year simulation period, only five wet winter seasons account for 33% of the total cumulative groundwater recharge.