Abstract: Mathematical groundwater modelling with homogeneous permeability zones has been used for decades to manage water resources in the Almonte-Marismas aquifer (southwest Spain). This is a highly heterogeneous detrital aquifer which supports valuable ecological systems in the Doñana National Park. The present study demonstrates that it is possible to better characterize this heterogeneity by numerical discretization of the geophysical and lithological data available. We identified six hydrofacies whose spatial characteristics were quantified with indicator variogram modelling. Sequential Indicator Simulation then made it possible to construct a 3D geological model. Finally, this detailed model was included in MODFLOW through the Model Muse interface. This final process is still a challenge due to the difficulty of downscaling to a handy numerical modelling scale. New piezometric surfaces and water budgets were obtained. The classical model with zones and the model with 3D simulation were compared to confirm that, for management purposes, the effort of improving the geological heterogeneities is worthwhile. This paper also highlights the relevance of including subsurface heterogeneities within a real groundwater management model in the present global change scenario.

Abstract: The current global intrusion of seawater into coastal aquifers causes salinization of groundwater and thus significant degradation of its quality. This study quantified the effect of seawater intrusion and freshening events in coastal aquifers on trace elements (Li, B, Mn and Ba) across the fresh-saline water interface (FSI) and their possible use as indicators for these events. This was done by combining field data and column experiments simulating these events. The experiments enabled quantification of the processes affecting the trace element composition and examination of whether salinization and freshening events are geochemically reversible, which has been seldom investigated. The dominant process affecting trace element composition during salinization and freshening is ion exchange. The results of the experiments show that the concentrations of major cations and Li+ were reversible during salinization and freshening, whereas B, Mn2+ and Ba2+ were not. During salinization, Li+ and B were depleted due to sorption by 10 and 100μmol·L−1, respectively, to about half of their expected conservative concentrations. The relative depletion of Li+ increased with distance from the shore, representing the propagation of salinization. Ba2+ and Mn2+ were desorbed from the sediment during salinization and enriched by tenfold in the aqueous phase compared to their concentration in seawater ( 0.1 μeq·L−1). During freshening both were depleted by almost tenfold compared to their concentration in fresh groundwater ( 0.7 μeq·L−1). The depletion of Mn2+ is a sensitive marker for freshening because Mn2+ has a strong affinity to the solid phase. Moreover, this study shows that both Mn2+ and Ba2+ can be used as sensitive hydrogeochemical tools to distinguish between salinization and freshening events in the FSI zone in coastal aquifers.