[1] Climate change in combination with increased anthropogenic activities will affect coastal groundwater systems throughout the world. In this paper, we focus on a coastal groundwater system that is already threatened by a relatively high seawater level: the low-lying Dutch Delta. Nearly one third of the Netherlands lies below mean sea level, and the land surface is still subsiding up to 1 m per century. This densely populated delta region, where fresh groundwater resources are used intensively for domestic, agricultural, and industrial purposes, can serve as a laboratory case for other low-lying delta areas throughout the world. Our findings on hydrogeological effects can be scaled up since the problems the Dutch face now will very likely be the problems encountered in other delta areas in the future. We calculated the possible impacts of future sea level rise, land subsidence, changes in recharge, autonomous salinization, and the effects of two mitigation countermeasures with a three-dimensional numerical model for variable density groundwater flow and coupled solute transport. We considered the effects on hydraulic heads, seepage fluxes, salt loads to surface waters, and changes in fresh groundwater resources as a function of time and for seven scenarios. Our numerical modeling results show that the impact of sea level rise is limited to areas within 10 km of the coastline and main rivers because the increased head in the groundwater system at the coast can easily be produced though the highly permeable Holocene confining layer. Along the southwest coast of the Netherlands, salt loads will double in some parts of the deep and large polders by the year 2100 A.D. due to sea level rise. More inland, ongoing land subsidence will cause hydraulic heads and phreatic water levels to drop, which may result in damage to dikes, infrastructure, and urban areas. In the deep polders more inland, autonomous upconing of deeper and more saline groundwater will be responsible for increasing salt loads. The future increase of salt loads will cause salinization of surface waters and shallow groundwater and put the total volumes of fresh groundwater volumes for drinking water supply, agricultural purposes, industry, and ecosystems under pressure.Citation: Oude Essink, G. H. P., E. S. van Baaren, and P. G. B. de Louw (2010), Effects of climate change on coastal groundwater systems: A modeling study in the Netherlands, Water Resour. Res., 46, W00F04,

The complex conductivity of soils remains poorly known despite the growing importance of this method in hydrogeophysics. In order to fill this gap of knowledge, we investigate the complex conductivity of 71 soils samples (including four peat samples) and one clean sand in the frequency range 0.1 Hz to 45 kHz. The soil samples are saturated with six different NaCl brines with conductivities (0.031, 0.53, 1.15, 5.7, 14.7, and 22 S m−1, NaCl, 25°C) in order to determine their intrinsic formation factor and surface conductivity. This data set is used to test the predictions of the dynamic Stern polarization model of porous media in terms of relationship between the quadrature conductivity and the surface conductivity. We also investigate the relationship between the normalized chargeability (the difference of in‐phase conductivity between two frequencies) and the quadrature conductivity at the geometric mean frequency. This data set confirms the relationships between the surface conductivity, the quadrature conductivity, and the normalized chargeability. The normalized chargeability depends linearly on the cation exchange capacity and specific surface area while the chargeability shows no dependence on these parameters. These new data and the dynamic Stern layer polarization model are observed to be mutually consistent. Traditionally, in hydrogeophysics, surface conductivity is neglected in the analysis of resistivity data. The relationships we have developed can be used in field conditions to avoid neglecting surface conductivity in the interpretation of DC resistivity tomograms. We also investigate the effects of temperature and saturation and, here again, the dynamic Stern layer predictions and the experimental observations are mutually consistent.

Abstract. In deltaic areas with saline seepage, freshwater availability is often limited to shallow rainwater lenses lying on top of saline groundwater. Here we describe the characteristics and spatial variability of such lenses in areas with saline seepage and the mechanisms that control their occurrence and size. Our findings are based on different types of field measurements and detailed numerical groundwater models applied in the south-western delta of the Netherlands. By combining the applied techniques we could extrapolate measurements at point scale (groundwater sampling, temperature and electrical soil conductivity (TEC)-probe measurements, electrical cone penetration tests (ECPT)) to field scale (continuous vertical electrical soundings (CVES), electromagnetic survey with EM31), and even to regional scale using helicopter-borne electromagnetic measurements (HEM). The measurements show a gradual mixing zone between infiltrating fresh rainwater and upward flowing saline groundwater. The mixing zone is best characterized by the depth of the centre of the mixing zone D mix , where the salinity is half that of seepage water, and the bottom of the mixing zone B mix , with a salinity equal to that of the seepage water (Cl-conc. 10 to 16 g l −1 ). D mix is found at very shallow depth in the confining top layer, on average at 1.7 m below ground level (b.g.l.), while B mix lies about 2.5 m b.g.l. The model results show that the constantly alternating upward and downward flow at low velocities in the confining layer is the main mechanism of mixing between rainwater and saline seepage and determines the position and extent of the mixing zone (D mix and B mix ). Recharge, seepage flux, and drainage depth are the controlling factors.

No abstract

Abstract. Coastal groundwater reserves often reflect a complex evolution of marine transgressions and regressions, and are only rarely in equilibrium with current boundary conditions. Understanding and managing the present-day distribution and future development of these reserves and their hydrochemical characteristics therefore requires insight into their complex evolution history. In this paper, we construct a paleo-hydrogeological model, together with groundwater age and origin calculations, to simulate, study and evaluate the evolution of groundwater salinity in the coastal area of the Netherlands throughout the last 8.5 kyr of the Holocene. While intended as a conceptual tool, confidence in our model results is warranted by a good correspondence with a hydrochemical characterization of groundwater origin. Throughout the modeled period, coastal groundwater distribution never reached equilibrium with contemporaneous boundary conditions. This result highlights the importance of historically changing boundary conditions in shaping the present-day distribution of groundwater and its chemical composition. As such, it acts as a warning against the common use of a steady-state situation given present-day boundary conditions to initialize groundwater transport modeling in complex coastal aquifers or, more general, against explaining existing groundwater composition patterns from the currently existing flow situation. The importance of historical boundary conditions not only holds true for the effects of the large-scale marine transgression around 5 kyr BC that thoroughly reworked groundwater composition, but also for the more local effects of a temporary gaining river system still recognizable today. Model results further attest to the impact of groundwater density differences on coastal groundwater flow on millennial timescales and highlight their importance in shaping today's groundwater salinity distribution. We found free convection to drive large-scale fingered infiltration of seawater to depths of 200 m within decades after a marine transgression, displacing the originally present groundwater upwards. Subsequent infiltration of fresh meteoric water was, in contrast, hampered by the existing density gradient. We observed discontinuous aquitards to exert a significant control on infiltration patterns and the resulting evolution of groundwater salinity. Finally, adding to a long-term scientific debate on the origins of groundwater salinity in Dutch coastal aquifers, our modeling results suggest a more significant role of pre-Holocene groundwater in the present-day groundwater salinity distribution in the Netherlands than previously recognized. Though conceptual, comprehensively modeling the Holocene evolution of groundwater salinity, age and origin offered a unique view on the complex processes shaping groundwater in coastal aquifers over millennial timescales.