Holocene marine transgressions are often put forward to explain observed groundwater salinities that extend far inland in deltas. This hypothesis was also proposed in the literature to explain the large land-inward extent of saline groundwater in the Nile Delta. The groundwater models previously built for the area used very large dispersivities to reconstruct this saline and brackish groundwater zone. However, this approach cannot explain the observed freshening of this zone. Here, we investigated the physical plausibility of the Holocene-transgression hypothesis to explain observed salinities by conducting a palaeohydrogeological reconstruction of groundwater salinity for the last 32 ka with a complex 3-D variable-density groundwater flow model, using a state-of-the-art version of the SEAWAT computer code that allows for parallel computation. Several scenarios with different lithologies and hypersaline groundwater provenances were simulated, of which five were selected that showed the best match with the observations. Amongst these selections, total freshwater volumes varied strongly, ranging from 1526 to 2659 km 3 , mainly due to uncertainties in the lithology offshore and at larger depths. This range is smaller (1511-1989 km 3 ) when we only consider the volumes of onshore fresh groundwater within 300 m depth. In all five selected scenarios the total volume of hypersaline groundwater exceeded that of seawater. We also show that during the last 32 ka, total freshwater volumes significantly declined, with a factor ranging from 2 to 5, due to the rising sea level. Furthermore, the time period required to reach a steady state under current boundary conditions exceeded 5.5 ka for all scenarios. Finally, under highly permeable conditions the marine transgression simulated with the palaeohydrogeological reconstruction led to a steeper fresh-salt interface compared to its steady-state equivalent, while low-permeable clay layers allowed for the preservation of fresh groundwater volumes. This shows that long-term transient simulations are needed when estimating present-day fresh-salt groundwater distributions in large deltas. The insights of this study are also applicable to other major deltaic areas, since many also experienced a Holocene marine transgression.

No abstract

Floating fen ecosystems are home to several protected habitats and species. Their development and conservation require special attention regarding water management. Although they are known to be heterogeneous and partially buoyant, their root mats are simulated in hydrological models as homogeneous, static systems. The objective of this study is to quantify root mat heterogeneity and buoyancy and to assess their effects on groundwater flow and transport, and to determine if these factors need to be taken into account in modelling. We conducted field measurements of root mat heterogeneity and buoyancy in the 'Nieuwkoopse Plassen', the Netherlands. We found that hydraulic conductivity varied over four orders of magnitude and negatively correlated with degree of decomposition, resulting in a zonation of high conductivity near the surface and low conductivity in the deeper layers. Also, we found that the root mat moved vertically with the surface water. It became more buoyant with higher temperatures, but less buoyant with increasing groundwater levels relative to the surface. We implemented the findings in a semi-steady state hydrological model of a floating fen to compare the effects with other parameters. The profound heterogeneity had a limited effect on the water budget, but a clear effect on the flow lines and thus should be taken into account when modelling transport processes in floating fens. Although buoyancy affected the relative groundwater level near the root mat edge, it did not affect the water budget or the flow lines and may therefore be neglected in water budget modelling.