Numerical study of groundwater flow cycling controlled by seawater/freshwater interaction in a coastal karst aquifer through conduit network using CFPv2

Assessment of confined aquifer response to recharge variations and water inflow distributions using analytical approach

Quantification of the amount of the exchanged water between the surface water and a confined aquifer is a basic step in water balance and environmental hydraulics. The hydraulic connection between a surface water and a confined aquifer may occur through different recharge variations. The current research presents new analytical solutions for confined aquifer response to recharge variations and different inflow distributions. Different cases were studied, where a constant piezometric head is applied at the right boundary of the 2D confined aquifer plane and various distributions of water inflow through the recharging windows are considered on a part and/or parts of the left boundary. Finally, a uniform water inflow distribution on parts of the left boundary and a uniform distribution of water outflow at the right boundary was considered. Both steady and unsteady state problems can be solved using proposed equations. The performance of developed analytical solutions was examined compared to the numerical finite difference modeling. The results show reasonable precision of the developed analytical solutions. The developed solutions can be used as a benchmark to verify numerical approaches with similar boundary conditions.

Managing aquifer recharge with multi-source water to realize sustainable management of groundwater resources in Jinan, China

Managed aquifer recharge (MAR) is an important approach to address water security, water quality decline, ground subsidence, and aquifer degradation. In this study, the large-scale recharge experiments were conducted in a natural river with multiple water sources. The MAR with multi-source water was investigated by developing an improved matter-element model under a limited recharged quantity and period in Jinan, China. Results showed that the background levels (BL) of groundwater quality before recharge was relatively good. However, the use of different water sources would cause a significant increase in the content of some groundwater quality indexes, which might further induce deterioration of regional groundwater quality. And the water quality in porous and karst aquifer displayed deteriorating trends during different water source recharge. Additionally, the adverse effects of recharge water sources on regional groundwater quality in turn was South-to-North Water Diversion Project (SN) > Yellow River (YR) > Wohushan Reservoir (WR). Meanwhile, the high-risk indexes in groundwater quality were presented during different water source recharge. Accordingly, relevant suggestions and measures were then put forward to optimize the MAR with multi-source water and explore the high-efficiency and low-risk recharge mode.

Taking into Account both Explicit Conduits and the Unsaturated Zone in Karst Reservoir Hybrid Models: Impact on the Outlet Hydrograph

The main outlets of karst systems are springs, the hydrographs of which are largely affected by flow processes in the unsaturated zone. These processes differ between the epikarst and transmission zone on the one hand and the matrix and conduit on the other hand. However, numerical models rarely consider the unsaturated zone, let alone distinguishing its subsystems. Likewise, few models represent conduits through a second medium, and even fewer do this explicitly with discrete features. This paper focuses on the interest of hybrid models that take into account both unsaturated subsystems and discrete conduits to simulate the reservoir-scale response, especially the outlet hydrograph. In a synthetic karst aquifer model, we performed simulations for several parameter sets and showed the ability of hybrid models to simulate the overall response of complex karst aquifers. Varying parameters affect the pathway distribution and transit times, which results in a large variety of hydrograph shapes. We propose a classification of hydrographs and selected characteristics, which proves useful for analysing the results. The relationships between model parameters and hydrograph characteristics are not all linear; some of them have local extrema or threshold limits. The numerous simulations help to assess the sensitivity of hydrograph characteristics to the different parameters and, conversely, to identify the key parameters which can be manipulated to enhance the modelling of field cases.

Modeling the Matrix-Conduit Exchanges in Both the Epikarst and the Transmission Zone of Karst Systems

Usual conceptual models of karst hydrodynamics highlight the important role of unsaturated subsystems in recharge repartition. However, few of them have been compared with scarce suitable physically-based numerical models. Hybrid models that couple single continuum medium with discrete features promise an improved consideration of karst specificities. Here we evaluate their capability to properly reproduce interactions between a vertical conduit and the surrounding unsaturated matrix. We simulate the response of such a configuration to a single recharge event for various sets of parameters. We show the ability of hybrid models to reproduce the most significant behaviors described in the literature, i.e., transient storage and distribution of recharge, flow concentration towards conduits in the epikarst, and matrix-conduit exchanges varying in time and space. In addition to the explicit conduits, simulating variably saturated flows with the Richards equation and distinguishing the epikarst and the transmission zone are key elements to reproduce most processes. The contrasts between subsystems necessary to observe desired behaviors have been quantified. They are reinforced by the varying matrix saturation that causes realistic competition between matrix and explicit conduits. The study also highlights the need to deepen knowledge of the scaled medium properties we need to know to apply such models to actual cases.

Coupling SKS and SWMM to Solve the Inverse Problem Based on Artificial Tracer Tests in Karstic Aquifers

Artificial tracer tests constitute one of the most powerful tools to investigate solute transport in conduit-dominated karstic aquifers. One can retrieve information about the internal structure of the aquifer directly by a careful analysis of the residence time distribution (RTD). Moreover, recent studies have shown the strong dependence of solute transport in karstic aquifers on boundary conditions. Information from artificial tracer tests leads us to propose a hypothesis about the internal structure of the aquifers and the effect of the boundary conditions (mainly high or low water level). So, a multi-tracer test calibration of a model appeared to be more consistent in identifying the effects of changes to the boundary conditions and to take into consideration their effects on solute transport. In this study, we proposed to run the inverse problem based on artificial tracer tests with a numerical procedure composed of the following three main steps: (1) conduit network geometries were simulated using a pseudo-genetic algorithm; (2) the hypothesis about boundary conditions was imposed in the simulated conduit networks; and (3) flow and solute transport were simulated. Then, using a trial-and-error procedure, the simulated RTDs were compared to the observed RTD on a large range of simulations, allowing identification of the conduit geometries and boundary conditions that better honor the field data. This constitutes a new approach to better constrain inverse problems using a multi-tracer test calibration including transient flow.

Analysis of the Saltwater Wedge in a Coastal Karst Aquifer with a Double Conduit Network, Numerical Simulations and Sensitivity Analysis

We investigate the long-distance salinity in a dual permeability coastal karst aquifer with a double conduit network using a three-dimensional variable-density groundwater flow and multispecies transport SEAWAT model. Sensitivity analyses were used to evaluate the impact of the parameters and boundary conditions on the modeling saltwater wedge in a karstic aquifer situated in the Cuban land territory, including hydraulic conductivity, vertical anisotropy and salinity concentration; both in the conduits network and the fractured medium. These analyses indicated that hydraulic conductivity of the fractured medium and salt concentration were the ones that have a stronger effect on saltwater intrusion in a karstic aquifer. We also show results of the three-dimensional numerical simulations on groundwater salinity for different scenarios with the variabilities of the important parameters and compare results with electric conductivity profiles measured in a well.

Implications of groundwater development and seawater intrusion for sustainability of a Mediterranean coastal aquifer in Tunisia

Tunisia relies extensively on coastal groundwater resources that are pumped at unsustainable rates to support irrigated agriculture, causing groundwater drawdown and water quality problems due to seawater intrusion. It is imperative for the country to regulate future groundwater allocations and implement conservation strategies based on robust hydrogeological assessments to alleviate the adverse impacts of groundwater depletion. We developed a 3D transient density-dependent groundwater model by coupling MODFLOW-2000 and MT3DMS to improve understanding of seawater intrusion into the Korba aquifer in Tunisia. Results indicate that groundwater overexploitation since 1965 induced 5.15 Mm³/year of seawater inflow while reducing submarine discharge into the sea by about 9.74 Mm³/year as compared to the steady state water budget in 1965. Projecting withdrawals from 2014 up to 2050 results in a slow but extensive groundwater table decline forming a cone of depression 15 m below sea level. The seawater wedge under this business-as-usual scenario is expected to reach 1.8 km from the shoreline, causing significant mixing of the TDS-rich seawater in the aquifer system. The cone of depression under a 25% increase in groundwater withdrawal drops to about 20 m below sea level while the saltwater front reaches 2.5 km inland. Countering the seawater intrusion problem requires reducing groundwater pumping by 17 Mm³/year to push back the saltwater front along the coastline by about 25% over a 43-year period. Application of the presented generic groundwater simulation framework guides developing management strategies to mitigate seawater intrusion in the Korba coastal aquifer and similar areas.

Long distance seawater intrusion through a karst conduit network in the Woodville Karst Plain, Florida

Five periods of increased electrical conductivity have been found in the karst conduits supplying one of the largest first magnitude springs in Florida with water. Numerous well-developed conduit networks are distributed in the Woodville Karst Plain (WKP), Florida and connected to the Gulf of Mexico. A composite analysis of precipitation and electrical conductivity data provides strong evidence that the increases in conductivity are directly tied to seawater intrusion moving inland and traveling 11 miles against the prevailing regional hydraulic gradient from from Spring Creek Spring Complex (SCSC), a group of submarine springs at the Gulf Coast. A geochemical analysis of samples from the spring vent rules out anthropogenic contamination and upwelling regional recharge from the deep aquifer as sources of the rising conductivity. The interpretation is supported by the conceptual model established by prior researchers working to characterize the study area. This paper documents the first and longest case of seawater intrusion in the WKP, and also indicates significant possibility of seawater contamination through subsurface conduit networks in a coastal karst aquifer.