10Vadose zone oxygen dynamics control all subsurface redox reactions and play a decisive role in 11 maintaining groundwater quality. Although drying and wetting events are common in artificial 12 recharge, their effects on subsurface oxygen distribution are poorly documented. We monitored 13 oxygen concentration in the unsaturated zone in a mid-scale (1 m high) laboratory soil lysimeter, 14which was subjected to short wetting and drying cycles that simulated a highly permeable 38• Surface scraping results in an immediate but temporary increase in the infiltration rate 39• Quantifying small changes in space and time is vital for accurate pore-scale O 2 mapping 40 3

Surface water-groundwater interactions were studied in a coastal lagoon performing 180 seepage meter measurements and using heat as a tracer in 30 locations along a lagoon inlet. The direct seepage meter measurements were compared with the results from analytical solutions for the 1D heat transport equation in three different scenarios: (1) Homogeneous bulk thermal conductivity (Ke); (2) horizontal heterogeneity in Ke; and (3) horizontal and vertical heterogeneity in Ke. The proportion of fresh groundwater and saline recirculated lagoon water collected from the seepage experiment was used to infer the location of the saline wedge and its effect on both the seepage meter results and the thermal regime in the lagoon bed, conditioning the use of the thermal methods. The different scenarios provided the basis for a better understanding of the underlying processes in a coastal groundwater-discharging area, a key factor to apply the best-suited method to characterize such processes. The thermal methods were more reliable in areas with high fresh groundwater discharge than in areas with high recirculation of saline lagoon water. The seepage meter experiments highlighted the importance of geochemical water sampling to estimate the origin of the exchanged water through the lagoon bed.

Submarine groundwater discharge (SGD)—including terrestrial freshwater, density-driven flow at the saltwater–freshwater interface, and benthic exchange—can deliver nutrients to coastal areas, generating a negative effect in the quality of marine water bodies. It is recognized that water stable isotopes (18O and 2H) can be helpful tracers to identify different flow paths and origins of water. Here, we show that they can be also applied when assessing sources of nutrients to coastal areas. A field site near a lagoon (Ringkøbing Fjord, Denmark) has been monitored at a metric scale to test if stable isotopes of water can be used to achieve a better understanding of the hydrochemical processes taking place in coastal aquifers, where there is a transition from freshwater to saltwater. Results show that 18O and 2H differentiate the coastal aquifer into three zones: Freshwater, shallow, and deep saline zones, which corresponded well with zones having distinct concentrations of inorganic phosphorous. The explanation is associated with three mechanisms: (1) Differences in sediment composition, (2) chemical reactions triggered by mixing of different type of fluxes, and (3) biochemical and diffusive processes in the lagoon bed. The different behaviors of nutrients in Ringkøbing Fjord need to be considered in water quality management. PO4 underneath the lagoon exceeds the groundwater concentration inland, thus demonstrating an intra-lagoon origin, while NO3, higher inland due to anthropogenic activity, is denitrified in the study area before reaching the lagoon.

Abstract. Essentially all hydrogeological processes are strongly influenced by the subsurface spatial heterogeneity and the temporal variation of environmental conditions, hydraulic properties, and solute concentrations. This spatial and temporal variability generally leads to effective behaviors and emerging phenomena that cannot be predicted from conventional approaches based on homogeneous assumptions and models. However, it is not always clear when, why, how, and at what scale the 4D (3D + time) nature of the subsurface needs to be considered in hydrogeological monitoring, modeling, and applications. In this paper, we discuss the interest and potential for the monitoring and characterization of spatial and temporal variability, including 4D imaging, in a series of hydrogeological processes: (1) groundwater fluxes, (2) solute transport and reaction, (3) vadose zone dynamics, and (4) surface–subsurface water interactions. We first identify the main challenges related to the coupling of spatial and temporal fluctuations for these processes. We then highlight recent innovations that have led to significant breakthroughs in high-resolution space–time imaging and modeling the characterization, monitoring, and modeling of these spatial and temporal fluctuations. We finally propose a classification of processes and applications at different scales according to their need and potential for high-resolution space–time imaging. We thus advocate a more systematic characterization of the dynamic and 3D nature of the subsurface for a series of critical processes and emerging applications. This calls for the validation of 4D imaging techniques at highly instrumented observatories and the harmonization of open databases to share hydrogeological data sets in their 4D components.