Groundwater is a crucial source of fresh water for domestic, agricultural and industrial uses, and for maintaining aquaticand groundwater-dependent ecosystems. It is also of importance for identifying and securing sustainable use of energy resources and subsurface storage sites. Groundwater interacts with different parts of system earth, including the atmosphere, surface water, soil, geological environment and with many geological processes. This chapter provides an overview of groundwater systems in the Netherlands and the Dutch continental shelf at depths from the surface down to about 5 km. The overview shows how the groundwater systems are shaped by the interplay between natural mechanisms operating on various geological timescales and the impact of anthropogenic activities. Important mechanisms involved in the development of groundwater flow systems include differences in the elevation of the groundwater table (its topographic relief), contrasts in groundwater density and hydromechanical interaction of groundwater with geologic media. The topography-driven groundwater flow systems in the coastal dunes, Pleistocene ice-pushed ridges, and the southeastern part of the country contain important fresh groundwater resources of meteoric origin. These resources occur largely in unconsolidated sedimentary sequences of Holocene and Pleistocene to Neogene age. Natural and anthropogenic factors explain the Holocene history of salinization and seepage in the coastal zone. The large transboundary topography-driven groundwater flow system in the southeast of the Netherlands has developed since Miocene times. It induced freshening of groundwater to relatively great depths and cooling of subsurface temperatures. Case studies show the effects of shallow and deep fault zones on flow and chemical conditions of groundwater. Groundwater in older, pre-Paleogene to Carboniferous units outside the realm of topography-driven flow mostly consists of highly saline brines. Groundwater in these units also shows high overpressures in the northern offshore and northern and northeastern part of the Netherlands, while close to hydrostatic pressures prevail in the southern onshore and adjacent offshore area.This spatial difference reflects the differences in burial history and hydrogeological framework.

Deltas are complex and are among the most vulnerable landforms under climate change. Studying them collectively highlights common stressors that drive their most significant challenges. A holistic conceptual framing of a delta and its feeding river basin is fundamental to effective adaptation planning.

We discuss the various performance aspects of parallelizing our transient global-scale groundwater model at 30′′ resolution (30arcsec; °1/41km at the Equator) on large distributed memory parallel clusters. This model, referred to as GLOBGM, is the successor of our 5′ (5arcmin; °1/410km at the Equator) PCR-GLOBWB 2 (PCRaster Global Water Balance model) groundwater model, based on MODFLOW having two model layers. The current version of GLOBGM (v1.0) used in this study also has two model layers, is uncalibrated, and uses available 30′′ PCR-GLOBWB data. Increasing the model resolution from 5′ to 30′′ creates challenges, including increased runtime, memory usage, and data storage that exceed the capacity of a single computer. We show that our parallelization tackles these problems with relatively low parallel hardware requirements to meet the needs of users or modelers who do not have exclusive access to hundreds or thousands of nodes within a supercomputer. For our simulation, we use unstructured grids and a prototype version of MODFLOW 6 that we have parallelized using the message-passing interface. We construct independent unstructured grids with a total of 278 million active cells to cancel all redundant sea and land cells, while satisfying all necessary boundary conditions, and distribute them over three continental-scale groundwater models (168 million - Afro-Eurasia; 77 million - the Americas; 16 million - Australia) and one remaining model for the smaller islands (17 million). Each of the four groundwater models is partitioned into multiple non-overlapping submodels that are tightly coupled within the MODFLOW linear solver, where each submodel is uniquely assigned to one processor core, and associated submodel data are written in parallel during the pre-processing, using data tiles. For balancing the parallel workload in advance, we apply the widely used METIS graph partitioner in two ways: it is straightforwardly applied to all (lateral) model grid cells, and it is applied in an area-based manner to HydroBASINS catchments that are assigned to submodels for pre-sorting to a future coupling with surface water. We consider an experiment for simulating the years 1958-2015 with daily time steps and monthly input, including a 20-year spin-up, on the Dutch national supercomputer Snellius. Given that the serial simulation would require °1/44.5 months of runtime, we set a hypothetical target of a maximum of 16h of simulation runtime. We show that 12 nodes (32 cores per node; 384 cores in total) are sufficient to achieve this target, resulting in a speedup of 138 for the largest Afro-Eurasia model when using 7 nodes (224 cores) in parallel. A limited evaluation of the model output using the United States Geological Survey (USGS) National Water Information System (NWIS) head observations for the contiguous United States was conducted. This showed that increasing the resolution from 5′ to 30′′ results in a significant improvement with GLOBGM for the steady-state simulation when compared to the 5′ PCR-GLOBWB groundwater model. However, results for the transient simulation are quite similar, and there is much room for improvement. Monthly and multi-year total terrestrial water storage anomalies derived from the GLOBGM and PCR-GLOBWB models, however, compared favorably with observations from the GRACE satellite. For the next versions of GLOBGM, further improvements require a more detailed (hydro)geological schematization and better information on the locations, depths, and pumping rates of abstraction wells.

A significant increase in surface water salinization in low-lying deltas is expected globally due to saline groundwater exfiltration driven by rising sea levels and decreasing freshwater availability. Sustaining fresh water-dependent agriculture in such areas will entail an increased demand for fresh water flushing. Unfortunately, the flushing of surface water is not operationally optimised and results in excessive use of scarce freshwater. To meet the increased demand for flushing, while minimizing the need for diverted freshwater, new operational designs are required. This paper presents a novel network model based approach that uses De Saint Venant (SV) and Advection Dispersion (AD) equations to optimize multiple objectives on water level and salinity control using a Nonlinear Model Predictive Control (NMPC). The resulting NMPC problem is solved with a receding horizon implementation, where the nonlinear program (NLP) at each iteration is solved using state-of-the-art large scale interior point solver (IPOPT). We evaluate the performance of the proposed approach and compare it to the traditional fixed flushing for a representative Dutch polder. Firstly, the approach is shown to be capable of controlling the water level and salinity level in the polder. Secondly, the results highlight that the network of canals, which were originally made for drainage, could not be made sufficiently fresh with current intake capacity. A simple design approach was used to identify appropriate new capacities for two of the gates that allow optimal flushing to guarantee the required water level and salinity constraints.

The Nile Delta Aquifer (NDA) is threatened by salt water intrusion (SWI). This article demonstrates an approach for identifying critical salinity concentration zones using a three-dimensional (3D) variable-density groundwater flow model in the NDA. An innovative procedure is presented for the delineation of salinity concentration in 2010 by testing different simulation periods. The results confirm the presence of saline groundwater caused by SWI in the north of the NDA. In addition, certain regions in the east and southwest of the NDA show increased salinity concentration levels, possibly due to excessive groundwater extraction and dissolution of marine fractured limestone and shale that form the bedrock underlying the aquifer. The research shows that the NDA is still not in a state of dynamic equilibrium. The modeling instrument can be used for simulating future scenarios of SWI to provide a sustainable adaptation plan for groundwater resource.

As Egypt's population increases, the demand for fresh groundwater extraction will intensify. Consequently, the groundwater quality will deteriorate, including an increase in salinization. On the other hand, salinization caused by saltwater intrusion in the coastal Nile Delta Aquifer (NDA) is also threatening the groundwater resources. The aim of this article is to assess the situation in 2010 (since this is when most data is sufficiently available) regarding the available fresh groundwater resources and to evaluate future salinization in the NDA using a 3D variable-density groundwater flow model coupled with salt transport that was developed with SEAWAT. This is achieved by examining six future scenarios that combine two driving forces: increased extraction and sea level rise (SLR). Given the prognosis of the intergovernmental panel on climate change (IPCC), the scenarios are used to assess the impact of groundwater extraction versus SLR on the seawater intrusion in the Delta and evaluate their contributions to increased groundwater salinization. The results show that groundwater extraction has a greater impact on salinization of the NDA than SLR, while the two factors combined cause the largest reduction of available fresh groundwater resources. The significant findings of this research are the determination of the groundwater volumes of fresh water, brackish, light brackish and saline water in the NDA as a whole and in each governorate and the identification of the governorates that are most vulnerable to salinization. It is highly recommended that the results of this analysis are considered in future mitigation and/or adaptation plans.