Performance Evaluation of Multiple Groundwater Flow and Nitrate Mass Transport Numerical Models

Environmental Modeling & Assessment - Benchmarking of different numerical models simulating groundwater flow and contaminant mass transport is the aim of the present study, in order to...     

Terrain Discontinuity Effects in the Regional Flow of a Complex Karstified Aquifer

Karst aquifers are characterized by spatial heterogeneity due to the presence of highly permeable channels and conduits in low-permeable fractured rocks (matrix block). Recent studies have reported a close relationship between surface and subsurface water in karstic regions due to the water flow through a complicated network of paths formed by fracture intersections. Subsurface flow in karstified aquifers ranges between conduit flow, in large passages with relatively high flow velocities, and diffuse flow, in the matrix block where Darcy’s law is still valid. In this paper, we present the simulation of a complex karstified aquifer system in Crete, Greece, where the presence of main faults drastically affects the regional flow. A discrete fracture approach in conjunction with an equivalent porous medium approach was adopted to simulate the mixed flow in the area of interest. The simulation results have shown that the length and the orientation of the dominant faults, primarily during the rainy season, affect the flow field.     

Comparison of stochastic and deterministic methods for mapping groundwater level spatial variability in sparsely monitored basins

In sparsely monitored basins, accurate mapping of the spatial variability of groundwater level requires the interpolation of scattered data. This paper presents a comparison of deterministic interpolation methods, i.e. inverse distance weight (IDW) and minimum curvature (MC), with stochastic methods, i.e. ordinary kriging (OK), universal kriging (UK) and kriging with Delaunay triangulation (DK). The study area is the Mires Basin of Mesara Valley in Crete (Greece). This sparsely sampled basin has limited groundwater resources which are vital for the island’s economy; spatial variations of the groundwater level are important for developing management and monitoring strategies. We evaluate the performance of the interpolation methods with respect to different statistical measures. The Spartan variogram family is applied for the first time to hydrological data and is shown to be optimal with respect to stochastic interpolation of this dataset. The three stochastic methods (OK, DK and UK) perform overall better than the deterministic counterparts (IDW and MC). DK, which is herein for the first time applied to hydrological data, yields the most accurate cross-validation estimate for the lowest value in the dataset. OK and UK lead to smooth isolevel contours, whilst DK and IDW generate more edges. The stochastic methods deliver estimates of prediction uncertainty which becomes highest near the southeastern border of the basin.     

Utilizing Successive Linearization Optimization to Control the Saltwater Intrusion Phenomenon in Unconfined Coastal Aquifers in Crete,Greece

In the present work, a simulation-optimization method is employed in order to manage saltwater intrusion in two unconfined coastal aquifers in Crete, Greece. The optimization formulation seeks to maximize groundwater withdrawal rates while maintaining the saltwater intrusion front at the current location or inhibiting it at locations closer to the coast. A combination of a groundwater flow model (MODFLOW) with the Ghyben-Herzberg saltwater front approximation and a sequential implementation of the Simplex algorithm (GWM) are employed. The results show that under the current pumping strategies, the saltwater intrusion front will continue to move inland, posing a serious threat to the groundwater quality of these regions. Optimal groundwater withdrawal scenarios that take into consideration the water needs of the local communities and environmental concerns are presented and discussed. In both case studies, significant reductions in pumping are required in order for the saltwater intrusion front to retract closer to the shoreline.