Subsurface flow and contaminant transport in heterogeneous media: concepts, modelling, observations and challenges in applications such as risk assessment and remediation.



This session combines presentations on recent developments in understanding, measuring, and modeling subsurface flow and transport. We aim to include processes in both the saturated and unsaturated zones, as well as across boundaries at different scales. At the same time, we address unsolved problems related groundwater contamination management as risk assessment and remediation.

The correct quantification of transport processes, which occur at different spatial and temporal scales, is challenging. It strongly influences predicted spreading, dilution and mixing rates. However, dispersion, mixing and chemical reactions are local phenomena that strongly depend on the interplay between large-scale system heterogeneity and smaller-scale processes. Much effort has been placed in the fundamental understanding of these processes since they are of practical relevance to identify the fate of contaminants in surface and subsurface water that can affect human health and the environment. Particularly newly emerging contaminants of such as PFASs, pharmaceuticals, pesticides, or nanoparticles are increasingly being detected at low levels in surface- and groundwater. Some of these anthropogenic chemicals are potentially harmful can produce long-term adverse health effects even at very low levels of exposure.

The aim of this session is to discuss the effect of flow heterogeneity on transport at different scales, from pore scale up to catchment scale - including theory, modeling, laboratory and field experiments as well as applications. Our contributions deal with the questions: Is macrodispersivity a meaningful parameter? Under which conditions does spatially variable flow enhance mixing and chemical reactions? What is the role played by diffusive processes in modeling transport in porous media? How to upscale dispersion and reactive transport from pore to field-scale? What is the relation between ADE models and dynamic structures of catchment hydrology like travel time distributions? What are appropriate methods to characterize the relevant aquifer properties? What are the recent improvements in transport measurement technologies? What is the best way to physically and chemically characterize sites contaminated by anthropogenic chemicals? What is their mobility and persistence in both the unsaturated and unsaturated zones? How can we improve remediation through laboratory and field research?

The session is co-sponsored by the Groundwater Commission of IAHS.

Co-sponsored by IAHS
Convener: Alraune Zech | Co-conveners: Fritjof Fagerlund, Felipe de Barros, Antonio Zarlenga, Marco Dentz, Aldo Fiori, Tissa Illangasekare
vPICO presentations
| Fri, 30 Apr, 13:30–15:00 (CEST)

vPICO presentations: Fri, 30 Apr

Chairpersons: Alraune Zech, Fritjof Fagerlund, Felipe de Barros
Heterogeneity from Pore to Lab Scale
Joris Heyman, Daniel R. Lester, and Tanguy Le Borgne

Recent works have shown that laminar flows through porous media generate Lagrangian chaos at pore scale, with strong implications for a range of transport, reactive, and biological processes in the subsurface. The characterization and understanding of mixing dynamics in these opaque environments remains an outstanding challenge. We present a novel experimental technique based upon high-resolution imaging of the scalar signature produced by push-pull flows through various porous materials (beads, gravels, sandstones) at high Péclet number. We show that this method provides a direct image (see below) of the invariant unstable manifold of the chaotic flow, while allowing a precise quantification of the incompleteness of mixing at pore scale. In the limit of large Péclet numbers, we demonstrate that the decay rate of the scalar variance is directly related to the Lyapunov exponent of the chaotic flow. Thus, this new push-pull method has the potential to provide a complete characterization of chaotic mixing dynamics in a large class of opaque porous materials.


How to cite: Heyman, J., Lester, D. R., and Le Borgne, T.: Push-pull flows reveal the scalar signature of chaos in porous media, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-665, https://doi.org/10.5194/egusphere-egu21-665, 2021.

Hugo Sanquer, Joris Heyman, Tanguy Le Borgne, and Khalil Hanna

Solute transport in porous media plays a key role in a range of chemical and biological processes, including contaminant degradation, precipitation, dissolution and microbiological dynamics. Increasing evidences have shown that the conventional complete mixing assumption at the pore scale can lead to a strong overestimation of reaction rates. Recent 3D imaging experiments of mixing in porous media suggest that these pore scale chemical gradients may be sustained by chaotic mixing dynamics. However, the consequences of such chaotic mixing on reactive processes are unknown.

In this work, we use reactive transport experiments coupled to 3D imaging to investigate the impact of micro-scale chaotic flows on mixing-limited reactions in the fluid phase.  We use optical index matching and laser-induced fluorescence to characterize the pore scale distribution of reactive product concentration for a range of Peclet and Damkhöler numbers. We use these measurements to develop a reactive lamellar theory that quantifies the impact of pore scale chemical gradients induced by chaotic mixing on effective reaction rates. These results provide new perspectives for upscaling reactive transport processes in porous media.

How to cite: Sanquer, H., Heyman, J., Le Borgne, T., and Hanna, K.: Impact of chaotic mixing on reactive transport: experiment in porous media at high Pe and Da, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15447, https://doi.org/10.5194/egusphere-egu21-15447, 2021.

Oshri Borgman, Turuban Régis, Baudouin Géraud, Le Borgne Tanguy, and Méheust Yves

Solute mixing mediated by flow in porous media plays a significant role in controlling reaction rates in subsurface environments. In many practical cases, incomplete mixing—inhomogeneous solute concentrations—occurs at the pore-scale, limiting local and thus upscaled reaction rates, and renders their prediction based on effective dispersion coefficients derived from dispersion models (or by assuming Taylor-Aris dispersion) inaccurate. We perform solute transport experiments in transparent, quasi-two-dimensional, soil analog models to investigate the relationships between pore-scale solute dispersion and mixing under different flow conditions. We use Fluorescein as a conservative tracer and record its fluorescence intensity in monochrome images at fixed time intervals. We convert the fluorescence intensity to solute concentration fields based on a calibration curve obtained with various homogeneous solute concentrations and subsequently compute concentration gradients. Our images provide evidence for incomplete mixing at the pore-scale and show strong gradients transverse to the overall flow direction. We fit the mean longitudinal concentration profile to an analytical solution of the advection-dispersion equation and compute the effective longitudinal dispersion coefficient. Based on the lamellar mixing theory, we also infer an effective diffusion coefficient relevant to the mean concentration gradient’s dynamics. By comparing these two diffusion/dispersion coefficients in saturated flow conditions, we show that while their values are similar at low Péclet, their scaling behaviors as a function of Péclet are different. Hence, as pointed out by several previous studies, modeling reactive transport processes requires accounting for a mixing behavior driven by a diffusive process that cannot entirely be described by the solute dispersion coefficient. We extend this work by varying the saturation degree in the experiments and our samples' structural heterogeneity to investigate how flow desaturation and porous medium structure impact solute mixing.

How to cite: Borgman, O., Régis, T., Géraud, B., Tanguy, L. B., and Yves, M.: Impact of flow conditions on pore-scale solute mixing: experiments in heterogeneous 2D porous media, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8636, https://doi.org/10.5194/egusphere-egu21-8636, 2021.

Hiwa Abdullah, Huhao Gao, Martin Sauter, and Alexandru Tatomir

The newly developed kinetic interface sensitive (KIS) tracers have been the focus of research in the past decade, as a new method to determine the mobile interfacial area between immiscible fluids in porous media. An accurate and reliable interfacial area determination is crucial to several industrial applications and the geoscientific research.

In this work we investigate the relationship between the concentration breakthrough curves of the KIS tracer, consequently the specific interfacial area and the evolution of the mobile non-wetting-phase front.

Up to now, such laboratory experiments have been conducted only in columns, quasi-one-dimensional systems. In this study we consider two-dimensional domains filled with porous material where immiscible displacement of water by oil takes place. The presence of heterogenous inclusions leads to perturbations in the fluid interface and causes fingers. By means of numerical modelling we investigate these effects and the results will help as a basis in the design of a new two-dimensional flume setup.

An analysis is performed for different viscosity ratios, capillary numbers corresponding to different capillary pressure-saturation relationships, injection rates and geometrical heterogeneity. We found that the presence of higher or lower permeability inclusions have a significant but clearly distinct impact on the destruction and/or production of the fluid-fluid interfacial area. Lower permeability inclusions increase the overall area of the front, compared to a decrease in the overall area for higher permeability inclusions. By increasing the interfacial area an increase of the reactive tracer concentration is observed. The mobile interfacial area is evaluated at the front of the saturation profile by using a cut-off value from the saturation profile, and then the area of the mobile concentration of the reactive tracer is calculated.

How to cite: Abdullah, H., Gao, H., Sauter, M., and Tatomir, A.: Numerical simulation of kinetic interface sensitive tracers in two-phase flow in porous media applied in two-dimensional experimental setups, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-5242, https://doi.org/10.5194/egusphere-egu21-5242, 2021.

Ester Marafini, Michele La Rocca, Aldo Fiori, Ilenia Battiato, and Pietro Prestininzi

Limitations stemming from the employment of 2D models to investigate the properties of 3D flows in porous media are generally overlooked. In this study, the extent to which 2D modelling can be employed for the representation of genuinely 3D flows in porous media is quantified. To this scope, Representative Elementary Volume (REV) sizes of 2D and 3D media sharing the same porosity are compared. The spatial stationarity of several Quantities of Interest (QoIs) namely, porosity, permeability, mean and variance of velocity, is numerically evaluated. In order to extend conclusions to transport phenomena, the analysis of the velocity variance, which is closely associated to the hydrodynamic dispersion process, is included. Porous media adopted in this study are composed by spheres and disks in 3D and 2D domains respectively, where both 2D and 3D geometries are characterized by random locations. Specifically, for 3D random packings creation, a sphere packing generator program is used. Pore scale flow is simulated by means of the Lattice Boltzmann Model (LBM): the LBM is employed as a numerical flow solver to reproduce the Darcy's experiment through the aforementioned domains. The LBM represents a powerful tool to model flow in porous media and it is able to accurately predict flow paths, permeability and hydraulic conductivity. Hydraulic QoIs are analysed at steady state conditions. To this purpose, the flow velocity field is used to inspect stationarity. The quantitative approaches adopted in the REV assessment procedure allow one to determine the residual variability of the quantity associated to the REV and consequently the level of accuracy that the modeller wants to achieve with respect to the QoIs. Such criteria show that REV estimations through 2D models are much larger than their 3D counterparts. In conclusion, pore scale LBM simulations highlight that the 2D approach leads to inconsistent results, due to the profound difference between 2D and 3D porous flows.

How to cite: Marafini, E., La Rocca, M., Fiori, A., Battiato, I., and Prestininzi, P.: Suitability of 2D modelling to evaluate flow properties in 3D porous media, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8710, https://doi.org/10.5194/egusphere-egu21-8710, 2021.

Heterogeneity from Lab to Field Scale
Eugenio Pescimoro, Matteo Icardi, Federico Municchi, and Marco Bianchi

Due to the increasing challenges to preserve water quality and supply at global scale, groundwater flow modeling has become a tool of pivotal relevance for remediation, implementation of policies, and design of applications for recharge management. The strain towards faster and more reliable hydrogeological simulations triggered the development of upscaled and multi-scale approaches employing different diffusion and dispersion models that are still the object of much debate in the community. Our ongoing study focuses on the up scaling of solute transport through heterogeneous geological domains by means of an extensive three-dimensional simulation study, based on a new open-source C++ library, built on top of the well-know finite-volume library OpenFOAM®. We integrate the whole workflow, from geostatistical random field generators to flow and transport solvers with integrated post-processing capabilities. The robustness, scalability and flexibility of the library makes it suitable framework for the development, testing, and application of upscaling techniques.  

Being the subsurface inaccessible by nature, the appeal to geostatistical techniques is a well-established approach to construct a realistic domain for flow and transport simulations. However,  additional challenges are posed by the numerical simulation of highly heterogeneous materials. Indeed, the problem is twofold: on one side it is not always possible to characterize the heterogeneity in a deterministic way, while on the other side numerical methods which are effective for elliptic and parabolic equations solved over homogeneous domains might suffer in heterogeneous media. Both challenges were effectively tackled using the open-soruce library OpenFOAM whose implementation and capabilities will be illustrated. Preliminary results on flow and transport simulations performed on truncated pluri-Gaussian permeability fields will be shown and the influence of geostatistical metrics (e.g. correlation lengths, variance, geological entropy) on the flow and transport results (e.g. average velocity and breakthrough curves) analysed.

Extensions to variable-density, mobile-immobile, and multi-rate mass transfer formulations are also presented in the context of the EU project “SECURe”.

How to cite: Pescimoro, E., Icardi, M., Municchi, F., and Bianchi, M.: An open-source integrated solution for flow and transport in random heterogeneous porous media, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2290, https://doi.org/10.5194/egusphere-egu21-2290, 2021.

Qinzhuo Liao, Gang Lei, Dongxiao Zhang, and Shirish Patil

We present a new method to estimate the displacement covariance and macrodispersivity for solute transport in bounded formations. Here we use circulant embedding, which is based on the fast Fourier transform and is much more efficient than eigen-decomposition for the factorization of random spatial fields. We compute the displacement covariances using the analysis of variance approach and introduce an interpolation process to significantly reduce the number of forward simulations. Once the effect of each eigenvector on the displacement covariance is obtained, it is unnecessary to rerun the simulator for different spatial covariance functions or anisotropy ratios, which saves a large amount of computational effort. The proposed method is validated in various tests in two-dimensional and three-dimensional examples and accurately matches the results from the Monte Carlo simulation. It is found that the longitudinal dispersivity is not sensitive to the boundaries, while the transverse and vertical dispersivities are greatly affected. The method is applied to the Borden site and provides a better explanation of the observed data after considering the effect of vertical boundaries. These results show that our method could serve as a promising tool for studying and predicting the characteristics of solute transport in heterogeneous media.

How to cite: Liao, Q., Lei, G., Zhang, D., and Patil, S.: Scale-dependent dispersion for solute transport in bounded formations, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-1746, https://doi.org/10.5194/egusphere-egu21-1746, 2021.

Maria Prieto Espinoza, Sylvain Weill, and Gwenaël Imfeld

Groundwater quality is of increasing concern due to the ubiquitous release of micropollutants, often originating from surface water. Micropollutants comprise a wide range of substances such as pesticides, pharmaceuticals and personal care products (PCPs), and pose risks to groundwater contamination due to their high persistence and toxicity. Although biodegradation is a major process for the removal of organic contaminants in aquifers, the interplay of hydrogeochemical conditions, microbial diversity and micropollutant dissipation at low concentrations remains yet poorly understood. We developed here an integrative approach to understand and predict the factors affecting micropollutant dissipation within the surface-/ground-water transition zone. Compound-specific Isotope Analysis (CSIA) was used to evaluate micropollutant transformation based on changes in the ratio of stable isotopes (i.e., 13C/12C and 15N/14N). The responses of aquifer microbes – the key players during contaminant transformation – to micropollutant exposures was examined through biomolecular approaches, proving advantageous in combination with CSIA.


We examined the dissipation of a micropollutant mixture in two lab-scale aquifer systems fed with river water collected from an agricultural area, thus representing the highly reactive surface-/ground-water transition zone. The micropollutant mixture included legacy and currently used pesticides such as atrazine, terbutryn, S-metolachlor and metalaxyl. Caffeine and metformin were also examined as anthropogenic compounds with physico-chemical properties close to currently used pesticides. Changes in bacterial diversity was examined in both aquifer systems during variations of micropollutant exposures under static hydrological conditions. It is hypothesized that such variations may induce bacterial changes and thus alter micropollutant transformation pathways. To this end, three periods of micropollutant injections during 140 d were induced as follow: (i) a first pulse (about 25 µM) to identify dissipation processes and bacterial adaptation to micropollutants, (ii) a constant injection (2 pore volumes) at 10 fold lower concentrations (chronic exposure phase), and (iii) a second pulse injection (25 µM) to examine whether transformation of micropollutants was enhanced. Concentration breakthrough curves (BTCs) of atrazine, terbutryn and metaxyl showed sorption as the main dissipation process for the three periods, whereas both sorption and degradation were observed for caffeine and S-metolachlor. Carbon and nitrogen CSIA further supported the in situ transformation of caffeine and S-metolachlor (Δδ13C of ≥ 4‰ and ≥ 2‰, respectively), while no significant enrichment of 13C and 15N were observed for atrazine, terbutryn and metalaxyl (Δδ13C < 2‰). In parallel, surface-water microcosm experiments showed half-life times of atrazine, terbutryn and metalaxyl of >200 days. Microbial diversity is currently examined in pore water and sand samples. A numerical model is under development to improve the interpretation of micropollutant dissipation in the highly reactive surface-/ground-water transition zone based on concentrations, CSIA and bacteria diversity data obtained in this study. Altogether, our results demonstrated degradation capacity within the laboratory systems, mainly for caffeine and S-metolachlor, and highlight the persistence and risk to long-term groundwater contamination of both legacy and currently used pesticides.

How to cite: Prieto Espinoza, M., Weill, S., and Imfeld, G.: Reactive transport of micropollutants in laboratory aquifers: combining Compound-Specific Isotope Analysis (CSIA) and biomolecular approaches, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-7864, https://doi.org/10.5194/egusphere-egu21-7864, 2021.

Mónica Basilio Hazas, Francesca Ziliotto, Giorgia Marcolini, Massimo Rolle, and Gabriele Chiogna

Hydropeaking, an artificial flow regime consisting on strong and frequent river stage fluctuations, is known to have important effects on groundwater-surface water interaction. It influences the transient dynamics of water flow and also of solute and energy fluxes between aquifers and rivers. In this work, we focus on the effects of hydropeaking at multiple spatial and temporal scales. We start the investigation at the laboratory scale using quasi-two-dimensional flow-through experiments in which we can characterize  flow and transport mechanisms, as well as the topology of the flow field, at high spatial and temporal resolution. We measure and model the spatial moments, the dilution index and the Okubo-Weiss parameter of a transient plume, and find a correlation between changes in flow topology and mixing enhancement. We then investigate a two-dimensional field scale cross section representative of the Adige aquifer in North-East Italy, where two rivers differently affected by hydropeaking influence groundwater flow, and we investigate the system considering hourly and mean daily fluctuations in the river stage. We characterize the transient groundwater dynamics for this and for other aquifers affected by hydropeaking using the Townley number, analyzing the potentiality of such systems for chaotic advection. Finally, at regional scale we use a three-dimensional transient model to show how the Adige aquifer is differently affected by hydropeaking depending on dry and wet years. Moreover, we apply the continuous wavelet transform to identify the main temporal scales of variability detected in the groundwater fluctuations and how they change with time. Our work therefore highlights the relevance of the effect of hydropeaking on groundwater flow and transport processes, and its impact on flow topology and mixing enhancement at multiple spatial and temporal scales.

How to cite: Basilio Hazas, M., Ziliotto, F., Marcolini, G., Rolle, M., and Chiogna, G.: Effects of hydropeaking on groundwater mixing: from laboratory experiments to field scale observations, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-1533, https://doi.org/10.5194/egusphere-egu21-1533, 2021.

Seonggan Jang and Minjune Yang

A flow chamber experiment was conducted to investigate solute diffusion in a multi-layered aquifer-aquitard system. The two-dimensional flow chamber consisted of a finite thickness aquitard (kaolinite, 2 cm) bounded by two parallel aquifers at the top (unconfined aquifer, 2 cm) and bottom (confined aquifer, 4 cm). New Coccine (red dye) of 500 mg/L in the top aquifer and distilled water in the bottom aquifer were injected with 0.02 mL/min for 13 days. One-dimensional analytical solutions were developed for advection and diffusion through a finite aquitard and compared with the measured data to evaluate experimental validation. The simulated aquitard concentration profiles (E > 0.97) and breakthrough curve (E = 0.95) showed good agreement with the measured data. During the experiment, the penetration distance in the aquitard increased over time and the vertical concentration distribution showed a linear profile through the aquitard after 7 days of loading in the top aquifer, indicating steady-state diffusion. The New Coccine diffused across the aquitard to the bottom aquifer after 1 day of loading. The bottom aquifer concentration increased at early times and was maintained after 7 days of loading (11 ~ 12 mg/L). This study provides experimental validation of the developed analytical solutions and quantitatively evaluates contaminant occurrences of the confined aquifer through the aquitard.

How to cite: Jang, S. and Yang, M.: Diffusive contaminant transport in a multi-layered aquifer-aquitard system, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3782, https://doi.org/10.5194/egusphere-egu21-3782, 2021.

Kan Bun Cheng, Gedeon Dagan, and Avinoam Rabinovich

Characterization of spatially variable aquifer properties is a necessary first step towards modeling flow and transport. An emerging technique in hydraulic tomography, known as diffusivity tests, consist of injecting (or pumping) a volume of water through short segments of a well for a short time and measuring the travel time of the peak of the head signal at different points in the surrounding aquifer volume. In our stochastic model, the specific storage is assumed to be constant, while the hydraulic conductivity of the heterogeneous aquifer is modeled as a random lognormal field. The axi-symmetric anisotropic structure is characterized by a few parameters (logconductivity mean and variance and horizontal and vertical integral scales). The mean and variance of the peak travel time are then determined as a function of distance from an instantaneous source by solving the flow equation using a first-order approximation in the logconductivity variance. The mean travel time is recast in terms of the equivalent conductivity, which decreases from the harmonic mean near the source to the effective conductivity in uniform flow for a sufficiently large distance. Similarly, the variance drops from its maximum near the source to a small value.

A different type of tomographic test is the constant-rate pumping one. We propose to apply the first order stochastic approach to the data from the Boise Hydrogeophysical Research site (BHRS) to characterize the aquifer properties by estimating heterogeneity statistical parameters. Equivalent properties are first calculated by matching a homogeneous aquifer solution to the pointwise data to obtain a spatially varying hydraulic conductivity (Keq) and storativity (Ss,eq). Then the statistical properties of K and Ss are to be computed by a best fit between the theoretically derived statistical moments of the equivalent random properties (Keq, Ss,eq) and those from field measurements. Our preliminary results indicate that the proposed stochastic methodology is robust and reliable as well as computationally more efficient than the conventional hydraulic tomography techniques.

How to cite: Cheng, K. B., Dagan, G., and Rabinovich, A.: Stochastic Modeling of Diffusivity and Constant-Rate Pumping Tests in Heterogeneous Aquifers in a Tomographic Setup and Its Application to Field Measurements, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2680, https://doi.org/10.5194/egusphere-egu21-2680, 2021.

Field Scale Heterogeneity and Applications
Alraune Zech, Sabine Attinger, Alberto Bellin, Vladimir Cvetkovic, Gedeon Dagan, Marco Dentz, Peter Dietrich, Aldo Fiori, and Georg Teutsch

Six conceptually different transport models are applied to the MADE-1 field tracer experiement as a first major attempt for model comparison. The objective was to show that complex mass distributions in heterogeneous aquifers can be predicted without calibration of transport parameters - solely making use of structural and flow data.

The models differ in their conceptualization of the heterogeneous aquifer structure, computational complexity, and use of conductivity data obtained from various observation methods (Direct Push Injection Logging - DPIL, Grain Size Analysis, Pumping Tests and Flowmeter). They agree in the underlying physical transport processes, none of them considering mass transfer. Predictive capability is assessed by comparing results to observed longitudinal mass distributions of the MADE-1 experiment. We deal with data uncertainty indicated by decreasing rates of recovered mass by focusing the comparison on measures, such as peak location, position and shape of bulk mass and leading tail, and we do not normalize observation data.

Comparison of models reveals that the predictions of the solute plume agree reasonably well with observations if the models are underlined by a few parameters of close values: mean velocity, a parameter reflecting log-conductivity variability and a horizontal length scale related to conductivity spatial correlation. The robustness of the results implies that conservative transport models with appropriate conductivity upscaling strategies of various observation data provide reasonable predictions of plumes longitudinal mass distribution as long as key features are taken into account.

How to cite: Zech, A., Attinger, S., Bellin, A., Cvetkovic, V., Dagan, G., Dentz, M., Dietrich, P., Fiori, A., and Teutsch, G.: A Comparison of Six Transport Models of the MADE-1 Experiment Implemented with Different Types of Hydraulic Data, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9469, https://doi.org/10.5194/egusphere-egu21-9469, 2021.

Antonio Zarlenga, Maria Ines Di Dato, Claudia D'Angelo, and Alessandro Casasso

Open-loop shallow geothermal systems, which exploit shallow aquifers as a heat source or sink, have a great potential to reduce greenhouse gas emissions related to the heating and cooling of buildings. In order to limit the depletion of groundwater resources water is generally reinjected into the same aquifer after the heat exchange, as a consequence a thermal plume develops within the aquifer. Furthermore a share of the reinjected water may come back to the abstraction wells, inducing a progressive thermal alteration of the abstracted water temperature that may even result in the plant failure. This phenomenon, known as thermal recycling, strongly depends on the hydraulic conductivity of the aquifer. The design models commonly adopted in the practice assume a homogeneous domain with constant hydraulic conductivity, this assumption, however, is not realistic: neglecting the natural heterogeneity of hydraulic properties of the porous medium may result in large prediction errors.

In this study, we aim to quantify the impact of the different heat transport dynamics in aquifers on the thermal plume development. A stochastic model, which explicitly considers the spatial variability of the hydrological properties, such as the hydraulic conductivity, is developed for low enthalpy geothermal systems. The thermal breakthrough curve at the extraction well is obtained by applying a Lagrangian model and assuming a steady state velocity field. Relevant quantities of thermal recycling, such as the thermal breakthrough time, are adopted for the evaluation of the effects of the hydrogeological and geometrical parameters of the systems.

The results of our study emphasize how the correct representation of the aquifer heterogeneity is fundamental in the design of shallow geothermal systems and in the correct heat plume assessment.

How to cite: Zarlenga, A., Di Dato, M. I., D'Angelo, C., and Casasso, A.: Effects of the aquifer heterogeneity in heat production from low enthalpy geothermal systems, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-4417, https://doi.org/10.5194/egusphere-egu21-4417, 2021.

Tinesh Pathania, T Iype Eldho, and Andrea Bottacin-Busolin

The use of optimization approaches for designing in-situ groundwater bioremediation systems has been demonstrated in a number of previous studies under the assumption of homogenous soil. However, in real applications the soil is typically heterogeneous and knowledge of the spatial conductivity distribution is, to some degree, uncertain. Here, a systematic attempt is made to quantify the effect of soil heterogeneity on the optimal design of in-situ bioremediation systems. To determine the optimal placement of injection and extraction wells within the computational domain, the meshless element-free Galerkin method (EFGM) was coupled with particle swarm optimization (PSO), resulting in a simulation-optimization (S/O) model which is referred to as BIOEFGM-PSO. A hypothetical case study is considered where the design of an in-situ bioremediation system is optimized considering different degrees of heterogeneity of the porous medium. Heterogeneous conductivity fields are generated using a pseudo-random field generator with same mean and varying variance and correlation lengths. The BIOEFGM-PSO model was then applied to the different soil scenarios, and the resulting bioremediation costs were compared. Results show that the optimal placement of injection and extraction wells depends on the soil properties and, on average, heterogeneous soils have higher in-situ bioremediation costs compared with a homogeneous soil with the same mean conductivity. This highlights the importance of considering soil heterogeneity in designing cost-effective in-situ bioremediation systems, and further demonstrates the general applicability of the BIOEFGM-PSO model.

How to cite: Pathania, T., Eldho, T. I., and Bottacin-Busolin, A.: Effect of soil heterogeneity on the optimal design of in-situ groundwater bioremediation systems, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-1922, https://doi.org/10.5194/egusphere-egu21-1922, 2021.

Dilshan Bandara, Thomas Heinze, Mahinda Premathilake, and Stefan Wohnlich

Hydro-geochemistry of groundwater plays an important role in understanding the characteristics of a geothermal system. Mixing zones of geothermal deep waters and shallow groundwater can be identified through chemical distribution maps and help identify geothermal flow paths. The flow paths can be used to calculate the chemical values of the geothermal water leading to a characterization of the heat source. In combination with knowledge about regional structures, the geochemical distribution can further reveal unknown geothermal zones.

In the present study, the geochemical distribution of the groundwater is studied from samples collected from shallow and deep wells, with special reference to the regional structures present in the terrain. The study area was selected as a 20 • 20 km area centered around the Padiyathalawa hot spring field in Sri Lanka. From the results, two main geochemical anomalous zones are identified, especially with the increased values of electrical conductivity (EC), total dissolved solids (TDS), and Sulphate distribution maps. Those two zones include the hot spring itself as well as an area in ~10 km distance in the NE direction from the hot spring. Both zones are characterized by crosscutting structures of dolerite dykes and shear zones. Due to the shear zones, there are deep-seated fractures facilitating water flow from deeper layers towards the surface. This uprising water mixes with the shallow groundwater, affecting the general geochemical values of the shallow groundwater system.

Common minerals in Dolerite in Sri Lanka are Pyroxenes, Feldspar, Ilmenite, Magnetite, and Pyrite with minor amounts of other minerals. The increased EC values in both before mentioned zones relate with higher amounts of iron due to dissolution and mixing processes in regions with fractured Dolerite. Similarly, the increased concentration of Sulphates in the groundwater can be related to Pyrite from the fractured Dolerite, as microbial oxidization of Pyrite leads to origin of Sulphates. The increase of TDS can be interpreted as shallow water mixing with deep geothermal water, which contains a higher amount of minerals from the fractured dolerites.

The similar geochemical anomalies in those two zones can be associated with cross-cutting Dolerite dikes and existing faults in the shear zones at greater depth, subsequently mixing uprising deep geothermal water with shallow groundwater. A similar geochemistry and tectonic setting suggest similar flow paths from the underground and therefore also similar geothermal conditions at both spots. However, due to the rural and remote region, only one of the two before mentioned areas is known as a hot spring field. Thermal signatures dissipate much more quickly in the shallow groundwater than the mineral composition and might not be significant for measurement. Geochemical signatures of groundwater can therefore be a substantial help to locate geothermal springs, identify source mechanisms and characterize fluid flow paths.

How to cite: Bandara, D., Heinze, T., Premathilake, M., and Wohnlich, S.: Regional structures influencing the groundwater geochemistry around geothermal springs: A case study from Padiyathalawa, Sri Lanka, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2378, https://doi.org/10.5194/egusphere-egu21-2378, 2021.

PFAS and Groundwater Contamination
Eva Weidemann, René Lämmer, Thorsten Stahl, Bernd Göckener, Mark Bücking, Jörn Breuer, Janine Kowalczyk, Hildegard Just, and Matthias Gassmann

Per- and polyfluoroalkyl substances (PFASs) are fluorinated anthropogenic compounds of which numerous have a high thermal stability, high surface tension, water and oil repellence. PFASs are used e.g. as surfactants, in agrochemical products and in fire extinguishing foams. Due to negative effects to the environment and human health some compounds are already declared as Persistent Organic Pollutants (POPs) or Substances of Very High Concern (SVHC). The behaviour of PFASs in environmental matrices can differ dependent on the molecular structure and environmental conditions. Sorption, for example, depends on soil characteristics, compound properties such as the carbon chain length and the functional group. When PFASs are found in soils, the knowledge of factors affecting the leaching behaviour in soil can be helpful to assess plant availability and risk of groundwater contamination.

In some contamination cases, e.g. Rastatt/Baden-Baden, Germany, PFASs precursors such as polyfluoroalkyl phosphate diesters (diPAPs) are present in soil and are transformed into PFAAs (perfluoroalkyl acids), a sub-group of PFASs, which are stable in the environment. To study the behaviour of diPAPs in soil, 8:2 diPAP and 6:2 diPAP were applied to 50 cm soil-filled columns with a concentration of 1 mg per kg soil. For two years, the columns are watered on 3-5 days per week with an average weekly natural rainfall of Baden-Württemberg. The leaching water is collected and analysed every two weeks. After the first year, there was no detection of precursors in the leaching water. The main products were, in accordance with literature, PFPeA (Perfluoropentanoic acid) and PFHxA (Perfluorohexanoic acid) for the precursor 6:2 diPAP and PFHpA (Perfluoroheptanoic acid) and PFOA (Perfluorooctanoic acid) for the precursor 8:2 diPAP. The breakthrough peak of PFPeA (864 µg/l) and PFHxA (487 µg/l) was found in week nine after the start of the experiment, the major peak of PFHpA (124 µg/l) in week 15 and the peak of PFOA (303 µg/l) in week 25. Even though the major breakthrough of PFAAs was over after the first year, all of them are still found in leaching water in relatively constant concentrations: ~40 µg PFPeA/l, ~20 µg PFHxA/l (6:2 diPAP) and ~5 µg PFHpA/l, ~50 µg PFOA/l (8:2 diPAP). Thus, as a preliminary result after the first year, there is a higher mass of PFPeA and PFHxA in the leachate compared to PFHpA and PFOA, which can be led back to higher transformation rates of 6:2 diPAP than 8:2 diPAP.

After two years, the soil columns will be cut horizontally in sections, and the soil will be analysed for precursors and PFAAs in order to calculate a mass balance. This will provide more insight into the transformation and mobilisation behaviour of investigated diPAPs and PFAAs.

How to cite: Weidemann, E., Lämmer, R., Stahl, T., Göckener, B., Bücking, M., Breuer, J., Kowalczyk, J., Just, H., and Gassmann, M.: Soil Column Experiments to Study Leaching and Transformation Behaviour of 8:2 diPAP and 6:2 diPAP, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2700, https://doi.org/10.5194/egusphere-egu21-2700, 2021.

Georgios Niarchos, Dan Berggren Kleja, Lutz Ahrens, and Fritjof Fagerlund

Remediation of sites contaminated with per- and polyfluoroalkyl substances (PFASs) is key to reduce the contamination of drinking water sources and subsequent human exposure. PFAS production and use is increasingly being restricted worldwide with a reduction of point sources; however, legacy plumes are still posing a threat due to the persistence of these chemicals against degradation. One of the most widely studied soil remediation techniques for PFASs is stabilisation, which results in the long-term entrapment of the contaminants with the addition of fixation agents in the subsurface, aiming to prevent their leaching from soil to groundwater. In relation to this, the aim of this study was to identify the leaching behaviour of various PFASs in a treatment scenario using activated carbon. Silt loam soil sampled from central Sweden was used, as well as a mixture of the soil with activated carbon at 0.1% w/w. Spiked artificial groundwater was prepared with a mixture of 21 PFASs, at a total concentration of 1.4 μg mL-1. The sorption of PFASs to the solid phase was investigated using 15 cm-long column experiments under saturated conditions. Uniform packing of the material was validated through non-reactive tracer tests. The desorption behaviour after treatment was also investigated, by switching the inflow from contaminated to clean water after steady state was achieved. Analysis of the compounds was conducted using ultra performance liquid chromatography coupled with tandem mass spectrometry (UPLC-MS/MS). Results have shown significantly increased sorption in soil amended with activated carbons compared to the untreated soil. Additionally, there was a positive correlation between the length of the perfluorocarbon chain and sorption efficiency. The study is a step towards increasing our understanding on the efficiency and longevity of stabilisation with activated carbons as a remediation strategy for PFAS-contaminated soils and groundwater.

How to cite: Niarchos, G., Berggren Kleja, D., Ahrens, L., and Fagerlund, F.: Remediation of per- and polyfluoroalkyl substances (PFASs) contaminated soil and groundwater; evaluating the performance of activated carbon in column tests, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-1117, https://doi.org/10.5194/egusphere-egu21-1117, 2021.

Kevin Pettersson, Mona Mossadeghi Björklund, and Alexandru Tatomir

Per- and poly-fluoroalkyl substances represent a new group of anthropogenic environmental pollutants.
For instance, one of the most commonly used firefighting methods today especially within airports is the use of aqueous film forming foams (AFFFs). AFFFs contain PFAS, which give the foam the ability to easily spread over a large area. Some PFAS species have demonstrated adverse health effects already at low concentrations such as liver damage and reproductive harm. The strong chemical bonding between carbon and fluorine also results in PFAS to being resistant to thermal- and chemical degradation. Thus, they degrade very slowly in nature. In Sweden, the limit for PFAS in drinking water is 90 ng L-1. This means that water needs to be treated for PFAS in the water treatment plants in cases of contamination. where water with a concentration over the limit would be used.

Bålsta is a growing city, with is proximity to Sweden’s largest city and capital Stockholm and Sweden’s fourth largest city Uppsala. Both located about 50 kilometres away, makes Bålsta an attractive city for people to live in and commute to their job in in the larger cities. With this growth new drinking water sources are required in the future Vreta-Bålsta area is an esker with good aquifer properties. However, it is an industrial area, and in the area lies a landfill which has been used as a training ground for firefighting exercises. During these exercises AFFFs have been used and the area is polluted with an unknown amount of PFAS. As the area is planned to be used for artificial infiltration and extraction of groundwater, it needs to be tested for potential PFAS contamination. 

This study aims to investigate by means of numerical modelling different scenarios of PFAS transport from the contaminated. Using data about the topography, soil depth of the area together, hydraulic conductivity, porosity, precipitation and runoff a model of the area was constructed within GMS. The package MODFLOW was used to simulate the groundwater flow. Using the groundwater flow solution, with the PFAS transport is simulated with the package MT3DMS to produce a result of a possible spread of PFAS within the studied area to see which parts that could be contaminated or become contaminated in the future.

The results given from the model was that PFAS from the landfill would transport northeast bound. With the southern part of the area would stay clear from all contamination from the landfill.

How to cite: Pettersson, K., Mossadeghi Björklund, M., and Tatomir, A.: Numerical Modelling of PFAS Transport and Groundwater Flow in the Vreta-Bålsta Esker, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-6542, https://doi.org/10.5194/egusphere-egu21-6542, 2021.

Fritjof Fagerlund, Georgios Niarchos, Lutz Ahrens, Dan Berggren Kleja, Jonny Bergman, Anna Larsson, Gareth Leonard, Jim Forde, Erik Ribeli, Henning Persson, and Lijana Gottby

Due to the exceptional persistence and resistance to degradation of per- and polyfluoroalkyl substances (PFASs), novel technologies for in-situ treatment and remediation of these pollutants are urgently needed. While there is still a need for more evidence from well-documented field applications, a promising technique is the use of activated carbon (AC) sorbents that can immobilize PFASs in groundwater and thereby prevent further spreading of the contaminants.

In Arboga Sweden a small fire-fighting training area connected to aviation industry is contaminated by PFAS from aqueous film forming foams (AFFFs). This site has been characterized for PFAS contamination and hydrogeological parameters affecting the spreading of contaminants with the groundwater in a few smaller site investigations since 2016 and continuous monitoring since 2018. In November 2019 colloidal activated carbon (CAC) was injected in a pilot-scale test to study the capability of CAC to immobilize PFASs in a part of the contamination plume.

The complex geology of the site made the injection of CAC challenging and special measures had to be taken to avoid excessive preferential flow of the CAC particles even at low-pressure injection. The injection pattern was modified and CaCl2 was injected downstream of the CAC injection to reduce CAC mobility and create a defined zone of CAC intercepting the PFAS plume in the groundwater, thus acting like a PFAS-immobilizing permeable barrier.

PFAS concentrations were initially reduced by 74% (for a sum of 11 PFASs) directly downstream of the CAC-barrier. However, a few months later PFAS concentrations rebounded to levels equally high or higher than before CAC injection, after which the levels have been going down again. The reasons to the rebound are likely connected to seasonal changes and fluctuations in the groundwater flow directions, causing bypass of the permeable CAC barrier. Lessons learned from applying CAC injections at this field site include the key importance of understanding the groundwater flow patterns and its temporal variability. CAC was able to produce significant reduction in PFAS concentrations (74%), but only when the PFAS plume was properly intercepted. The results illustrate the challenges with application of permeable barrier techniques particularly at geologically complex field sites. At such sites, sorbents for immobilization of PFAS plumes in groundwater should be applied in the most straightforward location where a year-round interception of the plume can be obtained.

How to cite: Fagerlund, F., Niarchos, G., Ahrens, L., Berggren Kleja, D., Bergman, J., Larsson, A., Leonard, G., Forde, J., Ribeli, E., Persson, H., and Gottby, L.: PFAS immobilization using in-situ application of colloidal activated carbon at a geologically complex site, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8592, https://doi.org/10.5194/egusphere-egu21-8592, 2021.

Federica Di Giacinto, Miriam Berti, Luigi Carbone, Riccardo Caprioli, Valentina Colaiuda, Annalina Lombardi, Barbara Tomassetti, Paolo Tuccella, Gianpaolo De Iuliis, Adelina Pietroleonardo, Giuseppina Mascilongo, Ludovica Di Renzo, Nicola D’Alterio, and Nicola Ferri

Emerging contaminants could exert combined toxic effects, including synergetic and antagonistic ones, that cannot be identified by chemical analysis tools. The biological early warning systems (BEWS) perform a real-time and continuous (24 h) monitoring of physiological and/or behavioural parameters of organism alterations, potentially correlated to water pollution. They are based on the response of living sentinels (i.e. molluscs, algae, crustaceans, fish) to a contaminant or mixture of them. Early warnings can be sent by SMS, e-mail, etc. to operators, in order to activate response actions. Belonging to different trophic levels, the crustacean D. magna, the alga C. vulgaris and the mollusc P. casertanum have been used to control Gran Sasso-Sirente aquifer in three different locations. Drinkable water of Teramo province and irrigation water of L’Aquila have been continuously monitored by the commercial tools “bbe® Daphnia Toximeter (DTOX)” and “bbe® Algae Toximeter (ATOX)”, respectively. In Tirino river spring, a novel sensor “SmartShell” has registered for the first time the valve movements of the autochthonous bivalve. After the first testing period, DTOX and ATOX did not register any typology of alarms in the potable and irrigation water. The valve movements of P. casertanum have been examined through spectral analysis in order to evaluate the behavioural rhythms useful for further investigation on their alterations as early warnings. The objective has been to reinforce the aquifer protection by installing instruments internationally recognised as efficient tools and exploring new proposals for guaranteeing human and ecosystem health.

How to cite: Di Giacinto, F., Berti, M., Carbone, L., Caprioli, R., Colaiuda, V., Lombardi, A., Tomassetti, B., Tuccella, P., De Iuliis, G., Pietroleonardo, A., Mascilongo, G., Di Renzo, L., D’Alterio, N., and Ferri, N.: Aquatic organism behaviours at multiple trophic levels for water monitoring: the case of Gran Sasso - Sirente aquifer, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15955, https://doi.org/10.5194/egusphere-egu21-15955, 2021.

Nico Hachgenei, Lorenzo Spadini, Guillaume Nord, Marie-Christine Morel, Anne Lespine, Henri Mora, François Courvoisier, Jean-François Sutra, and Céline Duwig

Ivermectin (IVM) is one of the few pharmaceutics that are still used in a preventive, systematic manner in extensive cattle breeding in our study region in the Ardèche region (France), amongst others. It is an efficient antiparasitic agent with an extreme acute toxicity for most invertebrates, especially aquatic organisms like daphnia (ng/l), and is also highly toxic to different fish species (µg/l). Due to its strong sorption to soil and sediment and quick photodegradation, early environmental risk assessments (ERA) conclude a low risk for aquatic organisms. More recent studies conclude an inacceptable risk for daphnia and dung organisms. One of the critical parameters between these contradictory conclusions is IVM export from cow dung and transfer towards the streams.

The study region is characterized by a Mediterranean climate with a dry summer and intense convective storm events leading to regular flash flood events that coincide with the cattle treatment seasons in spring and autumn. The study region encompasses the Claduègne catchment which is part of the OHMCV observatory and the OZCAR and eLTER research infrastructures.

The key question concerning the risk for aquatic organisms is to what extent and in which conditions IVM is mobilized and transferred from cow dung to soil and river via surface runoff and percolation in this environment prone to rapid flow processes. We approach this question on the scale of 60*30*22 (L*W*D) cm3 intact soil mesocosms, for which we developed an adapted field sampling and laboratory experimentation case. Soil mescosms are collected in the Claduègne catchment. IVM is applied in form of spiked cow dung at realistic environmental concentrations before simulating several rainfall events, representative of this Mediterranean region. Runoff and drainage water are sampled for major anions (including Br- tracer), non-particulate organic carbon and IVM concentrations on a high temporal frequency in order to gain an insight on the intra- and inter-event dynamics of water and IVM transfer. Tested parameters include dung ageing, soils types, initial soil humidity and consecutive rainfall events.

The first results highlight the importance of runoff for the overall export of IVM on the event scale. Concerning the water flux, initial humidity is found to determine the runoff / drainage partitioning as well as the rapidity of percolation through the occurrence of preferential flow. In this context, hydrophobicity seems to play an important role.

How to cite: Hachgenei, N., Spadini, L., Nord, G., Morel, M.-C., Lespine, A., Mora, H., Courvoisier, F., Sutra, J.-F., and Duwig, C.: Ivermectin transfer through percolation and surface runoff from intact soil mesocosms under rainfall simulation, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-12343, https://doi.org/10.5194/egusphere-egu21-12343, 2021.