Among the advances made in analytical and numerical analysis methods to quantify groundwater/surface-water interaction, one methodology that stands out is the use of heat as an environmental tracer. A large data set of river and riverbed temperature profiles from the Aa River in Belgium has been used to examine the spatial-temporal variations of groundwater/surface-water interaction. Exchange fluxes were calculated with the numerical heat-transport code STRIVE. The code was applied in transient mode to overcome previous limitations of steady-state analysis, and allowed for the calculation of model quality. In autumn and winter the mean exchange fluxes reached −90 mm d⁻¹, while in spring and early summer fluxes were −42 mm d⁻¹. Predominantly gaining conditions occurred along the river reach; however, in a few areas the direction of flow changed in time. The river banks showed elevated fluxes up to a factor of 3 compared to the center of the river. Higher fluxes were detected in the upstream section of the reach. Due to the influence of exchange fluxes along the river banks, larger temporal variations were found in the downstream section. The exchange fluxes at the river banks seemed more driven by variable local exchange flows, while the center of the river was dominated by deep and steady regional groundwater flows. These spatial and temporal differences in groundwater/surface-water exchange show the importance of long-term investigations on the driving forces of hyporheic processes across different scales.

This study investigates the occurrence of greenhouse gases (GHGs) and the role of groundwater as an indirect pathway of GHG emissions into surface waters in a gaining stretch of the Triffoy River agricultural catchment (Belgium). To this end, nitrous oxide (N₂O), methane (CH₄) and carbon dioxide (CO₂) concentrations, the stable isotopes of nitrate, and major ions were monitored in river and groundwater over 8 months. Results indicated that groundwater was strongly oversaturated in N₂O and CO₂ with respect to atmospheric equilibrium (50.1 vs. 0.55 μg L⁻¹ for N₂O and 14,569 vs. 400 ppm for CO₂), but only marginally for CH₄ (0.45 vs. 0.056 μg L⁻¹), suggesting that groundwater can be a source of these GHGs to the atmosphere. Nitrification seemed to be the main process for the accumulation of N₂O in groundwater. Oxic conditions prevailing in the aquifer were not prone for the accumulation of CH₄. In fact, the emissions of CH₄ from the river were one to two orders of magnitude higher than the inputs from groundwater, meaning that CH₄ emissions from the river were due to CH₄ in-situ production in riverbed or riparian zone sediments. For CO₂ and N₂O, average emissions from groundwater were 1.5 × 10⁵ kg CO₂ ha⁻¹ year⁻¹ and 207 kg N₂O ha⁻¹ year⁻¹, respectively. Groundwater is probably an important source of N₂O and CO₂ in gaining streams but when the measures are scaled at catchment scale, these fluxes are probably relatively modest. Nevertheless, their quantification would better constrain nitrogen and carbon budgets in natural systems.

Faults may impact regional groundwater flow and transport, so it is important to include them during aquifer-system conceptualization and while constructing groundwater flow models. For the Neogene aquifer in Flanders (Belgium), three-dimensional groundwater-flow and solute-transport models have been developed in the framework of safety and feasibility studies for the underlying Boom Clay Formation as a potential host for geological disposal of radioactive waste. The model outcomes are subject to uncertainties as they are typically constrained only by hydraulic heads, and their current conceptualization does not differentiate the fault zones from the undisturbed hydrostratigraphic formations. A groundwater flow model has been developed using MODFLOW-2005 to investigate how groundwater flow in the sedimentary Neogene aquifer may be disturbed by the Rauw Fault—a 55-km-long normal fault—across the Nete catchment. The observed hydraulic gradient across the fault zone appears significant, with head differences of 1.5–2.0 m over a horizontal distance of 60 m. A simulated hydraulic-head difference of 2.4 m was achieved largely corresponding to the observed behavior. The Neogene aquifer, within the Nete catchment, seems to be composed of several local flow systems and potentially with a deeper more semiregional/intermediate flow system. Testing different fault configurations shows the impact on the local/semiregional flow system, with pronounced effects in the fault’s vicinity, and extending or narrowing the flow systems further away. These results demonstrate the importance of considering faults, or any other hydrogeological subsurface barrier/conduit, and suggests they should be accounted for in the general practice of subsurface activity impact assessment.

Temperature distribution and heat transport are studied in a coastal aquifer at De Panne in the western Belgian coastal plain. Field observations include temperature profiles of groundwater in the dunes and temperature measurements at the water table in a profile on the shore. Freshwater–saltwater distribution is known from previous studies. These are used to constrain a density-dependent model simulating the freshwater–saltwater distribution and heat transport using the SEAWAT code. The yearly fluctuation of the groundwater temperature in the phreatic aquifer under the dunes, shore and sea, and the influence of a tidal inlet in the dunes are simulated. The observations show that seawater temperature variations determine the temperature variations on the shore whereas atmospheric temperature changes determine this in the dunes. Yearly temperature fluctuations imposed at the water table propagate mainly vertically in the aquifer with only limited lateral influence. Heat transport is mainly convection dominated. Thickness of the surficial zone is determined by the amplitude of the groundwater temperature at the water table and the groundwater flow. Establishment of a tidal inlet in the dunes results in asymmetric temperature profiles under and in the vicinity of it.

Groundwater vulnerability maps can be combined with pollution hazards to assess risks of groundwater pollution. However, groundwater vulnerability maps are generally difficult to interpret because they differ according to the factors considered and the way they are combined. Here, starting from process-based concepts and criteria, a robust definition for groundwater vulnerability to pollution is discussed. A methodology is developed based on processes governing the fate of pollutants at the land surface (i.e. runoff and infiltration) and below ground (i.e. pollutant transport in the subsurface). Groundwater vulnerability is evaluated based on combination of the land surface hazard and the subsurface attenuation capacity. Land surface hazard is defined to consider direct and lateral infiltration capacity of pollutants, regardless of any subsurface attenuation capacity, which refers to any process that leads to pollutant mass reduction from the infiltration location to the water table. The concept of subsurface attenuation capacity is adapted to the case of groundwater intrinsic vulnerability assessment, considering three process-based vulnerability coefficients, which are the pollutant minimum travel time from the hazard location to the water table, the pollution duration at the water table, and the maximum concentration of pollutant discharging into the groundwater. The concepts are illustrated by applying the developed method (named APSÛ) for intrinsic groundwater vulnerability assessment in the Néblon catchment, a karstified limestone/sandstone aquifer system in Belgium. The APSÛ method results are discussed and the perspectives for generalizing the method to groundwater-specific vulnerability and risk mapping are presented.