Gualbert Oude Essink

Improving fresh groundwater supply—problems and solutions

Many coastal regions in the world experience an intensive salt water intrusion in aquifers due to natural and anthropogenic causes. The salinisation of these groundwater systems can lead to a severe deterioration of the quality of existing fresh groundwater resources. In this paper, the characteristics of variable-density groundwater flow and the concept of an interface between fresh and saline groundwater are introduced. Problems associated with these groundwater resources are discussed, such as upconing of saline groundwater caused by excessive overpumping and salt water intrusion caused by global mean sea level rise. Possible human activities to compensate and control the salinisation of coastal aquifers are presented. Most countermeasures appear to be expensive, laborious and should be taken in time. Finally, the effect of lowering piezometric heads on the salinisation of the subsoil is illustrated in a Dutch aquifer system.

Salt Water Intrusion in a Three-dimensional Groundwater System in The Netherlands: A Numerical Study

Salt water intrusion is investigated in a coastal groundwater system in the northern part of the province Noord-Holland, The Netherlands. Density dependent groundwater flow is modeled in three-dimensions with MOCDENS3D. This computer code is a version of MOC3D (Konikow et al., 1996) that has been adapted to simulate transient density-driven groundwater flow. Results from the model suggests that in this Dutch hydrogeologic system a severe and irreversible salinisation is already occurring. Within a few tens to hundreds of years, the salinity of the shallow aquifer is estimated to increase substantially. This salinisation process is a result of human activities such as the reclamation of the low-lying areas during the past centuries. Without changing the present boundary conditions, seepage into the low-lying areas will decrease slightly because of predicted increases in groundwater salinity. However, the rate in salt load through the Holocene aquitard into the low-lying areas will increase significantly due to an increase in salinity in the shallow aquifer. In addition, a relative sea level rise of 0.5 m per century will intensify the salinisation process, causing an enormous increase in salt load in all low-lying areas in this part of The Netherlands.Salt water intrusion is investigated in a coastal groundwater system in the northern part of the province Noord-Holland, The Netherlands. Density dependent groundwater flow is modeled in three-dimensions with MOCDENS3D. This computer code is a version of MOC3D (Konikow et al., 1996) that has been adapted to simulate transient density-driven groundwater flow. Results from the model suggests that in this Dutch hydrogeologic system a severe and irreversible salinisation is already occurring. Within a few tens to hundreds of years, the salinity of the shallow aquifer is estimated to increase substantially. This salinisation process is a result of human activities such as the reclamation of the low-lying areas during the past centuries. Without changing the present boundary conditions, seepage into the low-lying areas will decrease slightly because of predicted increases in groundwater salinity. However, the rate in salt load through the Holocene aquitard into the low-lying areas will increase significantly due to an increase in salinity in the shallow aquifer. In addition, a relative sea level rise of 0.5 m per century will intensify the salinisation process, causing an enormous increase in salt load in all low-lying areas in this part of The Netherlands.

An operational, multi-scale, multi-model system for consensus-based, integrated water management and policy analysis: The Netherlands Hydrological Instrument

Water management in the Netherlands applies to a dense network of surface waters for discharge, storage and distribution, serving highly valuable land-use. National and regional water authorities develop long-term plans for sustainable water use and safety under changing climate conditions. The decisions about investments on adaptive measures are based on analysis supported by the Netherlands Hydrological Instrument NHI based on the best available data and state-of-the-art technology and developed through collaboration between national research institutes. The NHI consists of various physical models at appropriate temporal and spatial scales for all parts of the water system. Intelligent connectors provide transfer between different scales and fast computation, by coupling model codes at a deep level in software. A workflow and version management system guarantees consistency in the data, software, computations and results. The NHI is freely available to hydrologists via an open web interface that enables exchange of all data and tools.

Global sampling to assess the value of diverse observations in conditioning a real‐world groundwater flow and transport model

The use of additional types of observational data has often been suggested to alleviate the ill-posedness inherent to parameter estimation of groundwater models and constrain model uncertainty. Disinformation in observational data caused by errors in either the observations or the chosen model structure may, however, confound the value of adding observational data in model conditioning. This paper uses the global generalized likelihood uncertainty estimation methodology to investigate the value of different observational data types (heads, fluxes, salinity, and temperature) in conditioning a groundwater flow and transport model of an extensively monitored field site in the Netherlands. We compared model conditioning using the real observations to a synthetic model experiment, to demonstrate the possible influence of disinformation in observational data in model conditioning. Results showed that the value of different conditioning targets was less evident when conditioning to real measurements than in a measurement error-only synthetic model experiment. While in the synthetic experiment, all conditioning targets clearly improved model outcomes, minor improvements or even worsening of model outcomes was observed for the real measurements. This result was caused by errors in both the model structure and the observations, resulting in disinformation in the observational data. The observed impact of disinformation in the observational data reiterates the necessity of thorough data validation and the need for accounting for both model structural and observational errors in model conditioning. It further suggests caution when translating results of synthetic modeling examples to real-world applications. Still, applying diverse conditioning data types was found to be essential to constrain uncertainty, especially in the transport of solutes in the model.

Fresh Water Lens Persistence and Root Zone Salinization Hazard Under Temperate Climate

In low lying deltaic areas in temperate climates, groundwater can be brackish to saline at shallow depth, even with a yearly rainfall excess. For primary production in horticulture, agriculture, and terrestrial nature areas, the fresh water availability may be restricted to so-called fresh water lenses: relatively thin pockets of fresh groundwater floating on top of saline groundwater. The persistence of such fresh water lenses, as well as the quantity and quality of surface water is expected to be under pressure due to climate change, as summer droughts may intensify in North-West Europe. Better understanding through modelling of these fresh water resources may help anticipate the impact of salinity on primary production. We use a simple model to determine in which circumstances fresh water lenses may disappear during summer droughts, as that could give rise to enhanced root zone salinity. With a more involved combination of expert judgement and numerical simulations, it is possible to give an appraisal of the hazard that fresh water lenses disappear for the Dutch coastal regions. For such situations, we derive an analytical tool for anticipating the resulting salinization of the root zone, which agrees well with numerical simulations. The provided tools give a basis to quantify which lenses are in hazard of disappearing periodically, as well as an impression in which coastal areas this hazard is largest. Accordingly, these results and the followed procedure may assist water management decisions and prioritization strategies leading to a secure/robust fresh water supply on a national to regional scale.

Saltwater Upconing Due to Cyclic Pumping by Horizontal Wells in Freshwater Lenses.

This article deals with the quantification of saltwater upconing below horizontal wells in freshwater lenses using analytical solutions as a computationally fast alternative to numerical simulations. Comparisons between analytical calculations and numerical simulations are presented regarding three aspects: (1) cyclic pumping; (2) dispersion; and (3) finite horizontal wells in a finite domain (a freshwater lens). Various hydrogeological conditions and pumping regimes within a dry half year are considered. The results show that the influence of elastic and phreatic storage (which are not taken into account in the analytical solutions) on the upconing of the interface is minimal. Furthermore, the analytical calculations based on the interface approach compare well with numerical simulations as long as the dimensionless interface upconing is below 1/3, which is in line with previous studies on steady pumping. Superimposing an analytical solution for mixing by dispersion below the well over an analytical solution based on the interface approach is appropriate in case the vertical flow velocity around the interface is nearly constant but should not be used for estimating the salinity of the pumped groundwater. The analytical calculations of interface upconing below a finite horizontal well compare well with the numerical simulations in case the distance between the horizontal well and the initial interface does not vary significantly along the well and in case the natural fluctuation of the freshwater lens is small. In order to maintain a low level of salinity in the well during a dry half year, the dimensionless analytically calculated interface upconing should stay below 0.25.

Fast calculation of groundwater exfiltration salinity in a lowland catchment using a lumped celerity/velocity approach

To support operational water management of freshwater resources in coastal lowlands, a need exists for a rapid, well-identifiable model to simulate salinity dynamics of exfiltrating groundwater. This paper presents the lumped Rapid Saline Groundwater Exfiltration Model (RSGEM). RSGEM simulates groundwater exfiltration salinity dynamics as governed by the interplay between water velocity, gradually adjusting the subsurface salinity distribution, and pressure wave celerity, resulting in a fast flow path response to groundwater level changes. RSGEM was applied to a field site in the coastal region of the Netherlands, parameter estimation and uncertainty analysis were performed using generalized likelihood uncertainty estimation. The model showed good correspondence to measured groundwater levels, exfiltration rates and salinity response. Moreover, RSGEM results were very similar to a detailed, complex groundwater flow and transport model previously applied to this field site.