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.

Abstract. Management of coastal fresh groundwater reserves requires a thorough understanding of the present-day groundwater salinity distribution and its possible future development. However, coastal groundwater often still reflects a complex history of marine transgressions and regressions, and is only rarely in equilibrium with current boundary conditions. In addition, the distribution of groundwater salinity is virtually impossible to characterize satisfactorily, complicating efforts to model and predict coastal groundwater flow. A way forward may be to account for the historical development of groundwater salinity when modeling present-day coastal groundwater flow. In this paper, we construct a palaeo-hydrogeological model to simulate the evolution of groundwater salinity in the coastal area of the Netherlands throughout the Holocene. While intended as a perceptual tool, confidence in our model results is warranted by a good correspondence with a hydrochemical characterization of groundwater origin. Model results 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. Not once reaching steady-state throughout the Holocene, our results demonstrate the long-term dynamics of salinity in coastal aquifers. This stresses the importance of accounting for the historical evolution of coastal groundwater salinity when modeling present-day coastal groundwater flow, or when predicting impacts of e.g. sea level rise on coastal aquifers. Of more local importance, our findings suggest a more significant role of pre-Holocene groundwater in the present-day groundwater salinity distribution in the Netherlands than previously recognized. The implications of our results extend beyond understanding the present-day distribution of salinity, as the proven complex history of coastal groundwater also holds important clues for understanding and predicting the distribution of other societally relevant groundwater constituents.

No abstract