Dense non-aqueous phase liquids, i.e., DNAPLs and the evolving contaminant plumes in aquifers provide significant potential to pose hazards affecting both environment and human health. Therefore, a proper assessment of contaminant spreading within the subsurface is critical. This includes a sufficient characterization of governing parameters describing both the subsurface and the contaminant itself. Thereby, knowledge on the contaminant source zone and especially the source zone geometry, i.e., SZG is critically required, yet very uncertain. This study identifies current limitations and open research questions in the formation and shape determination of source zone geometry, as well as its relevance for contaminant plumes. Our literature review reveals that existing characterization methods are subject to data interpretation uncertainties, while the application of these methods on field scale is limited by technical demands and accompanied efforts. In a next step, methods to implement increased source zone information into calculation methods are discussed. By means of an exemplary application of selected assessment tools, i.e., plume response models, results clearly proof the relevance of SZGs for site assessment. However, existing plume response models consider over-simplified geometries that may compromise their suitability. Our findings identify the demand for improved characterization of complex SZGs and the need to better evaluate the dependency of DNAPL migration on system properties and external influences. With emphasized knowledge on the most relevant SZG features, the delineation of "effective" SZGs allowing for straightforward implementation into plume response models and an adaption of the latter to incorporate more information on SZGs should be possible.

We conducted multiple laboratory trials in a robust and repeatable experimental layout to study dense non-aqueous phase liquid (DNAPL) source zone formation. We extended an image processing and analysis framework to derive DNAPL saturation distributions from reflective optical imaging data, with volume balance deviations < 5.07%. We used a multiphase flow model to simulate source zone formation in a Monte Carlo approach, where the parameter space was defined by the variation of retention curve parameters. Integral and geometric measures were used to characterize the source zones and implemented into a multi-criteria objective function. The latter showed good agreement between observation data and simulation results for effective DNAPL saturation values > 0.04, especially for early stages of DNAPL migration. The common hypothesis that parameters defining the DNAPL-water retention curves are constant over time was not confirmed. Once DNAPL pooling started, the optimal fit in the parameter space was significantly different compared to the earlier DNAPL migration stages. We suspect more complex processes (e.g., capillary hysteresis, adsorption) to become relevant during pool formation. Our results reveal deficits in the grayscale-DNAPL saturation relationship definition and laboratory estimation of DNAPL-water retention curve parameters to overcome current limitations to describe DNAPL source zone formation.

Subsurface DNAPL (dense non-aqueous phase liquid) contamination from (un-) intentional spilling typically leads to severe environmental hazards. A large number of studies have demonstrated the relevance of DNAPL source zone geometry for the determination of contaminant plume propagation in groundwater. Optical imaging represents a promising non-invasive method for identifying DNAPL saturation without disturbing multiphase flow dynamics. However, workflow and image analysis methodologies have not been sufficiently developed or described for general application to related experimental efforts. For example, the choice of dye(s) used for phase colorization affects image processing and can bias final estimations of DNAPL saturations. In this study, we perform a series of DNAPL migration and entrapment studies in transparent tanks that are filled with three different types of porous media. Different dyes are used and raw images are acquired. Subsequently, these are used to evaluate a suite of image processing and analysis approaches, which are organized into a workflow. Our approach allows for us to identify key image processing and analysis steps that introduce the most error. Applicable dye configurations led to uncertainties of up to 41% depending on the selection of processing steps. Based on these findings, it was possible to delineate a flexible framework for image processing and analysis that has the potential for transfer and application in other tank experiment setups.

Leaky urban drainage networks (UDNs) exfiltrating wastewater can contaminate aquifers. Detailed knowledge on spatiotemporal distributions of water-dissolved, sewer-borne contaminants in groundwater is essential to protect urban aquifers and to optimize monitoring systems. We evaluated the effect of UDN layouts on the spreading of sewer-borne contaminants in groundwater using a parsimonious approach. Due to the UDN’s long-term leakage behavior and the existence of non-degradable sewer-borne contaminants (equivalent to a conservative and constant contaminant source), we employed a concept of horizontal line sources to mimic the UDN layout. This does not require the consideration of bio-degradation processes or temporal delay and effectively bypasses the vadose zone, thus reducing computational requirements associated with a full simulation of leakages. We used a set of synthetic leakage scenarios which were generated using fractals and are based on a real-world UDN layout. We investigated the effects of typical leakage rates, varying groundwater flow directions, and UDN’s layouts on the shape of the contaminant plume, disregarding the resulted concentration. Leakage rates showed minimal effects on the total covered plume area, whereas 89% of the variance of the plume’s geometry is explained by both the UDN’s layout (e.g., length and level of complexity) and groundwater flow direction. We demonstrated the potential of applying this approach to identify possible locations of groundwater observation wells using a real UDN layout. This straightforward and parsimonious method can serve as an initial step to strategically identify optimal monitoring systems locations within urban aquifers, and to improve sewer asset management at city scale.

In this study, the influence of varying temperatures on the transport behavior of conservative solutes in heterogeneous porous media has been investigated. Column flow experiments employing potassium bromide as tracer were conducted at four temperature levels (3°C, 10°C, 20°C, and 30°C) and the measured electrical conductivity (EC) signal was interpreted through inverse modeling. Additional experiments were performed with deuterium‐enriched water at 10°C and 30°C. For those experiments, deuterium isotope ratios were analyzed alongside anion and cation concentrations. Obtained EC‐based breakthrough curves showed measurable differences in both the observed peak values and tailing intensities that could be clearly attributed to variations in the experimental temperatures. The EC‐based results were further corroborated by the measured isotope ratios and corresponding anion/cation concentrations, although measured differences were less pronounced. The model‐based interpretation of the results employed the standard advection‐dispersion equation as well as three alternative variants that were all based on the single‐rate dual‐domain mass transfer (DDMT) approach, but embracing varying coupling intensities between the experiments. For all variants transport parameters were determined for EC, bromide, and deuterium, respectively. The estimated ranges of the transport parameters all point toward a direct correlation between the effective DDMT rates and the experimental temperature. The observed correlation directly follows the Arrhenius relationship, but is weaker than the one describing the temperature dependence of the molecular diffusion coefficients, therefore pointing to a contribution of non‐diffusive components.