Simulation of saturated and unsaturated flow in karst systems at catchment scale using a double continuum approach

The objective of this work is the simulation of saturated and unsaturated flow in a karstified aquifer using a double continuum approach. The HydroGeoSphere code (Therrien et al., 2006) is employed to simulate spring discharge with the Richards equations and van Genuchten parameters to represent flow in the (1) fractured matrix and (2) conduit continuum coupled by a linear exchange term. Rapid vertical small-scale flow processes in the unsaturated conduit continuum are accounted for by applying recharge boundary conditions at the bottom of the saturated model domain. An extensive sensitivity analysis is performed on single parameters as well as parameter combinations. The transient hydraulic response of the karst spring is strongly controlled by the matrix porosity as well as the van Genuchten parameters of the unsaturated matrix, which determine the head dependent inter-continuum water transfer when the conduits are draining the matrix. Sensitivities of parameter combinations partially reveal a non-linear dependence over the parameter space. This can be observed for parameters not belonging to the same continuum as well as combinations, which involve the exchange parameter, showing that results of the double continuum model may depict a certain degree of ambiguity. The application of van Genuchten parameters for simulation of unsaturated flow in karst systems is critically discussed

Joint inversion of groundwater flow, heat, and solute state variables: A multipurpose approach for characterization and forecast of karst systems

Characterization of karst systems and forecast of their state variables are essential for groundwater management and engineering in karst regions. These objectives can be met by the use of process-based discrete-continuum models (DCMs). However, results of DCMs may suffer from inversion nonuniqueness. It has been demonstrated that the joint inversion of observations regulated by different natural processes can tackle the nonuniqueness issue in groundwater modeling. However, this has not been tested for DCMs thus far. This research proposes a methodology for the joint inversion of hydro-thermo-chemo-graphs, applying to two small-scale sink-to-spring experiments at Freiheit Spring, Minnesota, USA. In order to address conceptual uncertainty, a multimodel approach was implemented, featuring seven mutually exclusive variants. Spring hydro-thermo-chemo-graphs, for all the variants simulated by MODFLOW-CFPv2, were jointly inverted using a weighted least squares algorithm. Subsequently, models were compared in terms of inversion and forecast performances, as well as parameter uncertainties. Results reveal the suitability of the DCM approach for simultaneous inversion and forecast of hydro-physico-chemical behavior of karst systems, even at a scale of meters and seconds. The estimated volume of the tracer conduit passage ranges from approximately 46–51 m³, which is comparable to the estimate from the flood-pulse method. Moreover, it was demonstrated that the thermograph and hydrograph contain more information about aquifer characteristics than the chemograph. However, this finding can be site-specific and should depend on the analysis scale, the considered conceptual models, and the hydrological state, which are potentially affected by minor unaccountable processes and features

Global karst springs hydrograph dataset for research and management of the world's fastest-flowing groundwater [Data paper]

Karst aquifers provide drinking water for 10% of the world's population, support agriculture, groundwater-dependent activities, and ecosystems. These aquifers are characterised by complex groundwater-flow systems, hence, they are extremely vulnerable and protecting them requires an in-depth understanding of the systems. Poor data accessibility has limited advances in karst research and realistic representation of karst processes in large-scale hydrological studies. In this study, we present World Karst Spring hydrograph (WoKaS) database, a community-wide effort to improve data accessibility. WoKaS is the first global karst springs discharge database with over 400 spring observations collected from articles, hydrological databases and researchers. The dataset's coverage compares to the global distribution of carbonate rocks with some bias towards the latitudes of more developed countries. WoKaS database will ensure easy access to a large-sample of good quality datasets suitable for a wide range of applications: comparative studies, trend analysis and model evaluation. This database will largely contribute to research advancement in karst hydrology, supports karst groundwater management, and promotes international and interdisciplinary collaborations

Assessing geotechnical centrifuge modelling in addressing variably saturated flow in soil and fractured rock

The vadose zone typically comprises soil underlain by fractured rock. Often, surface water and groundwater parameters are readily available, but variably saturated flow through soil and rock are oversimplified or estimated as input for hydrological models. In this paper, a series of geotechnical centrifuge experiments are conducted to contribute to the knowledge gaps in: (i) variably saturated flow and dispersion in soil and (ii) variably saturated flow in discrete vertical and horizontal fractures. Findings from the research show that the hydraulic gradient, and not the hydraulic conductivity, is scaled for seepage flow in the geotechnical centrifuge. Furthermore, geotechnical centrifuge modelling has been proven as a viable experimental tool for the modelling of hydrodynamic dispersion as well as the replication of similar flow mechanisms for unsaturated fracture flow, as previously observed in literature. Despite the imminent challenges of modelling variable saturation in the vadose zone, the geotechnical centrifuge offers a powerful experimental tool to physically model and observe variably saturated flow. This can be used to give valuable insight into mechanisms associated with solid–fluid interaction problems under these conditions. Findings from future research can be used to validate current numerical modelling techniques and address the subsequent influence on aquifer recharge and vulnerability, contaminant transport, waste disposal, dam construction, slope stability and seepage into subsurface excavations.The Water Research Commission of South Africahttp://link.springer.com/journal/113562018-05-31hb2017Geolog