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: Submarine groundwater discharge (SGD) is an important source of nutrients and metals to the coastal ocean, affects coastal ecosystems, and is gaining recognition as a relevant water resource. SGD is usually quantified using geochemical tracers such as radon or radium. However, a few studies have also used dissolved silicon (DSi) as a tracer for SGD, as DSi is usually enriched in groundwater when compared to surface waters. In this study, we discuss the potential of DSi as a tracer in SGD studies based on a literature review and two case studies from contrasting environments. In the first case study, DSi is used to calculate SGD fluxes in a tropical volcanic-carbonate karstic region (southern Java, Indonesia), where SGD is dominated by terrestrial groundwater discharge. The second case study discusses DSi as a tracer for marine SGD (i.e., recirculated seawater) in the tidal flat area of Spiekeroog (southern North Sea), where SGD is dominantly driven by tidal pumping through beach sands. Our results indicate that DSi is a useful tracer for SGD in various lithologies (e.g., karstic, volcanic, complex) to quantify terrestrial and marine SGD fluxes. DSi can also be used to trace groundwater transport processes in the sediment and the coastal aquifer. Care has to be taken that all sources and sinks of DSi are known and can be quantified or neglected. One major limitation is that DSi is used by siliceous phytoplankton and therefore limits its applicability to times of the year when primary production of siliceous phytoplankton is low. In general, DSi is a powerful tracer for SGD in many environments. We recommend that DSi should be used to complement other conventionally used tracers, such as radon or radium, to help account for their own shortcomings.