Articles | Volume 21, issue 7
Hydrol. Earth Syst. Sci., 21, 3635–3653, 2017
https://doi.org/10.5194/hess-21-3635-2017
Hydrol. Earth Syst. Sci., 21, 3635–3653, 2017
https://doi.org/10.5194/hess-21-3635-2017

Research article 18 Jul 2017

Research article | 18 Jul 2017

Assessing lateral flows and solute transport during floods in a conduit-flow-dominated karst system using the inverse problem for the advection–diffusion equation

Cybèle Cholet1,2, Jean-Baptiste Charlier2, Roger Moussa3, Marc Steinmann1, and Sophie Denimal1 Cybèle Cholet et al.
  • 1Chrono-Environnement, UMR 6249 UBFC/CNRS, University of Burgundy Franche-Comté, Besançon, 25000, France
  • 2BRGM, 1039 rue de Pinville, 34000 Montpellier, France
  • 3INRA, UMR LISAH, 2 Place Pierre Viala, 34060 Montpellier, France

Abstract. The aim of this study is to present a framework that provides new ways to characterize the spatio-temporal variability of lateral exchanges for water flow and solute transport in a karst conduit network during flood events, treating both the diffusive wave equation and the advection–diffusion equation with the same mathematical approach, assuming uniform lateral flow and solute transport. A solution to the inverse problem for the advection–diffusion equations is then applied to data from two successive gauging stations to simulate flows and solute exchange dynamics after recharge. The study site is the karst conduit network of the Fourbanne aquifer in the French Jura Mountains, which includes two reaches characterizing the network from sinkhole to cave stream to the spring. The model is applied, after separation of the base from the flood components, on discharge and total dissolved solids (TDSs) in order to assess lateral flows and solute concentrations and compare them to help identify water origin. The results showed various lateral contributions in space – between the two reaches located in the unsaturated zone (R1), and in the zone that is both unsaturated and saturated (R2) – as well as in time, according to hydrological conditions. Globally, the two reaches show a distinct response to flood routing, with important lateral inflows on R1 and large outflows on R2. By combining these results with solute exchanges and the analysis of flood routing parameters distribution, we showed that lateral inflows on R1 are the addition of diffuse infiltration (observed whatever the hydrological conditions) and localized infiltration in the secondary conduit network (tributaries) in the unsaturated zone, except in extreme dry periods. On R2, despite inflows on the base component, lateral outflows are observed during floods. This pattern was attributed to the concept of reversal flows of conduit–matrix exchanges, inducing a complex water mixing effect in the saturated zone. From our results we build the functional scheme of the karst system. It demonstrates the impact of the saturated zone on matrix–conduit exchanges in this shallow phreatic aquifer and highlights the important role of the unsaturated zone on storage and transfer functions of the system.

Download
Short summary
This paper aims to improve the understanding of transport processes in karst aquifers, which remains a great challenge due to its complex nature. A framework is proposed to identify and quantify the spatio-temporal variability of lateral exchanges along a karst conduit network during flood events, for both flow and solute transport. An inverse problem approach is used on discharge and water mineralization data sets and gives new insights into the hydrogeological behavior of such complex systems.