Identifying and determining hydraulic parameters of physically heterogeneous aquifers is pivotal for flow field analysis, contaminant migration and risk assessment. In this research, we applied a novel uniquely sequenced DNA tagged superparamagnetic silica microparticles (SiDNAmag) to quantify hydraulic parameters and associated uncertainties of a heterogeneous sand tank. In the sand tank with lens shaped heterogeneity, we conducted three sets of multi – point injection experiments in unconsolidated (1) homogeneous (zone 0), (2) heterogeneous with a no-conductivity-zone (zone 1), and (3) heterogeneous with a high-conductive-zone (zone 2). From the breakthrough curves (BTC), we estimated the parameters distributions of hydraulic conductivity (k), effective porosity (n_e), longitudinal dispersivity (α_L), transverse vertical (α_TV), and transverse horizontal dispersivities (α_TH) applying Monte Carlo simulation approach for BTC fitting. The estimated parameters and associated uncertainties for each of the heterogeneous sections were further statistically compared (distribution non-specific Mann Whitney U test) these parameter distributions with parameter distributions estimated from the conservative salt tracer. While the time of arrival and time to peak concentration of SiDNAmag and salt in effluent were comparable, peak concentration of SiDNAmag was 1–3 log reduced as compared to the salt tracer due to first order kinetic attachment. Nonetheless, the parameters and associated uncertainty distributions (5 %–95 %) of K, n_e, α_L, α_TV, and α_TH, determined from SiDNAmag BTCs were statistically equivalent to the salt tracer in all three experiment systems. Through our experimental and modelling approach, our work demonstrated that in a coarse to very coarse grain sand medium, with lens shaped heterogeneity, the uniquely sequenced SiDNAmag were a promising tool to identify heterogeneity and determine hydraulic parameters and associated uncertainty distributions.

In this study, we focused on the 3D dispersion of colloids. To our knowledge, we were the first to do so. Thereto, we injected silica encapsulated DNA tagged superparamagnetic particles (SiDNAmag) in a homogeneous coarse grain sand tank. At four downstream locations, SiDNAmag concentrations were determined as a function of time. Longitudinal and transverse dispersivity values and associated uncertainties of SiDNAmag were determined using Monte Carlo modelling approach. The parameter associated uncertainties of hydraulic conductivity as well as of the effective porosity estimated from SiDNAmag breakthrough curves were statistically similar to those estimated from salt tracer breakthrough curves. Further, the SiDNAmag dispersivity uncertainty ranges were then statistically compared with the salt tracer (NaCl, and fluorescein) dispersivities. Our results indicated that time to rise, time of peak concentration and shape of the breakthrough curves of SiDNAmag were similar to those of the salt tracer breakthrough curves. Despite the size difference between the salt tracer molecules and SiDNAmag, size exclusion did not occur, probably due to the large pore throat diameter to SiDNAmag diameter ratio. The median longitudinal dispersivity (α_L) of salt tracer and SiDNAmag were 4.9 and 5.8 × 10⁻⁴ m, respectively. The median ratio of horizontal and vertical transverse dispersivities to α_L, (α_TH /α_L and α_TV /α_L, respectively), for salt tracer and SiDNAmag ranged between 0.52 and 0.56. Through the statistical tests, we concluded that the longitudinal and traverse dispersivities of SiDNAmag were not statistically significantly different from salt tracer in 3 dimensions and could be used to characterize the dispersive properties of the medium we used. Our work contributes to a better understanding of 3D dispersion of SiDNAmag in saturated porous media.

We investigated the applicability of Silica encapsulated, superparamagnetic DNA particles (SiDNAmag) in determining aquifer hydraulic parameters at different ionic strengths (1 mM, 5 mM, and 20 mM phosphate buffer) of injection suspension. Thereto, in a homogeneous, unconsolidated sand tank we pulse - injected two uniquely sequenced SiDNAmag at two injection points. At 0.5 m and 0.8 m downstream from the injection points, we measured the concentration of SiDNAmags at three vertically distributed and two horizontally distributed sampling locations. We estimated the hydraulic parameter distributions from the SiDNAmag breakthrough curves through a Monte – Carlo approach and compared the parameter distributions with salt tracer breakthrough curves. Our results indicated that at all the ionic strengths, the times of peak concentrations, and the shapes of the breakthrough curves were similar to the salt tracer. As compared to the salt, a 1 – 3 log units reduction in the maximum effluent concentration of SiDNAmag was due to kinetic attachment. The attachment rate reduced from 1 mM to 5 mM phosphate buffer possibly due to competitive adsorption of phosphate onto the favourable attachment sites. SiDNAmag attachment rate further increased in 20 mM buffer suspension, possibly due to the compression of electric double layer and reduction in energy barrier for attachment. The parameter distributions of hydraulic conductivity (k), effective porosity (n_e), longitudinal dispersivity (α_L), vertical transverse dispersivity (α_TV /α_L) and horizontal transverse dispersivity (α_TH /α_L) estimated from the SiDNAmag and the salt tracer breakthrough curves were statistically similar. Our work contributes to the applicability of colloidal SiDNAmags for determining hydraulic parameters at different ionic strength conditions.

Agricultural aquifer storage recovery and transfer (ASTR) stores excess fresh water for later reuse in irrigation. Moreover, water quality improves because chemical pollutants and pathogens will be removed by degradation and attachment to the aquifer material. The source water may contain the bacterial plant pathogen Ralstonia solanacearum which causes plant infections and high yield losses. We used quantitative microbial risk assessment (QMRA) to investigate the removal of R. solanacearum during ASTR to predict infection risks of potato plants after irrigation with the recovered water. Laboratory experiments analyzed the ASTR treatment by investigating the bacterial die-off in the water phase and the removal by attachment to the aquifer sediment. Die-off in the water phase depends on the residence time and ranged between 1.3 and 2.7 log₁₀ after 10 or 60 days water storage, respectively. A subpopulation of the bacteria persisted for a prolonged time at low concentrations which may pose a risk if the water is recovered too early. However, the natural aquifer sand filtration proofed to be highly effective in removing R. solanacearum by attachment which depends on the distance between injection and abstraction well. The high removal by attachment alone (18 log₁₀ after 1 m) would reduce bacterial concentrations to negligible numbers. Upscaling to longer soil passages is discussed in the paper. Infection risks of potato plants were calculated using a dose-response model and ASTR treatment resulted in negligible infection risks of a single plant, but also when simulating the irrigation of a 5 ha potato field. This is the first QMRA that analyzed an agricultural ASTR and the fate of a plant pathogen focusing on plant health. QMRA is a useful (water) management tool to evaluate the treatment steps of water reclamation technologies with the aim to provide safe irrigation water and reduce risks disseminating plant diseases.

Irrigation with surface water carrying plant pathogens poses a risk for agriculture. Managed aquifer recharge enhances fresh water availability while simultaneously it may reduce the risk of plant diseases by removal of pathogens during aquifer passage. We compared the transport of three plant pathogenic bacteria with Escherichia coli WR1 as reference strain in saturated laboratory column experiments filled with quartz sand, or sandy aquifer sediments. E. coli showed the highest removal, followed by Pectobacterium carotovorum, Dickeya solani and Ralstonia solanacearum. Bacterial and non-reactive tracer breakthrough curves were fitted with Hydrus-1D and compared with colloid filtration theory (CFT). Bacterial attachment to fine and medium aquifer sand under anoxic conditions was highest with attachment rates of max. k _att1 = 765 day ^-1 and 355 day ^-1, respectively. Attachment was the least to quartz sand under oxic conditions (k _att1 = 61 day ^-1). In CFT, sticking efficiencies were higher in aquifer than in quartz sand but there was no differentiation between fine and medium aquifer sand. Overall removal ranged between < 6.8 log ₁₀ m ⁻¹ in quartz and up to 40 log ₁₀ m ⁻¹ in fine aquifer sand. Oxygenation of the anoxic aquifer sediments for two weeks with oxic influent water decreased the removal. The results highlight the potential of natural sand filtration to sufficiently remove plant pathogenic bacteria during aquifer storage.

Managed aquifer recharge (MAR) can provide irrigation water and overcome water scarcity in agriculture. Removal of potentially present plant pathogens during MAR is essential to prevent crop diseases. We studied the die-off of three plant pathogenic bacteria in water microcosms with natural or filtered tile drainage water (TDW) at 10 and 25◦C and with natural anoxic aquifer water (AW) at 10◦C from a MAR site. These bacteria were: Ralstonia solanacearum (bacterial wilt), and the soft rot Pectobacteriaceae (SRP) Dickeya solani and Pectobacterium carotovorum sp. carotovorum (soft rot, blackleg). They are present in surface waters and cause destructive crop diseases worldwide which have been linked to contaminated irrigation water. Nevertheless, little is known about the survival of the SRP in aqueous environments and no study has investigated the persistence of R. solanacearum under natural anoxic conditions. We found that all bacteria were undetectable in 0.1 mL samples within 19 days under oxic conditions in natural TDW at 10◦C, using viable cell counting, corresponding to 3-log₁₀ reduction by die-off. The SRP were no longer detected within 6 days at 25◦C, whereas R. solanacearum was detectable for 25 days. Whereas in anoxic natural aquifer water at 10◦C, the bacterial concentrations declined slower and the detection limit was reached within 56 days. Finally, we modelled the inactivation curves with a modified Weibull model that can simulate different curve shapes such as shoulder phenomena in the beginning and long tails reflecting persistent bacterial populations. The nonlinear model was shown to be a reliable tool to predict the die-off of the analysed plant pathogenic bacteria, suggesting its further application to other pathogenic microorganisms in the context of microbial risk assessment.

Soil passage of (pretreated) surface water to remove pathogenic microorganisms is a highly efficient process under oxic conditions, reducing microorganism concentrations about 8 log¹⁰ within tens of meters. However, under anoxic conditions, it has been shown that removal of microorganisms can be limited very much. Setback distances for adequate protection of natural groundwater may, therefore, be too short if anoxic conditions apply. Because removal of microorganisms under suboxic conditions is unknown, this research investigated removal of bacteriophage MS2 and PRD1 by soil passage under suboxic conditions at field scale. At the field location (dune area), one injection well and six monitoring wells were installed at different depths along three suboxic flow lines, where oxygen concentrations ranged from 0.4 to 1.7 mg/l and nitrate concentrations ranged from 13 to 16 mg/L. PRD1 and MS2 were injected directly at the corresponding depths and their removal in each flow line was determined. The highest bacteriophage removal was observed in the top layer, with about 9 log removal of MS2, and 7 log removal of PRD1 after 16 meters of aquifer transport. Less removal was observed at 12 m below surface, probably due to a higher groundwater velocity in this coarser grained layer. MS2 was removed more effectively than PRD1 under all conditions. Due to short travel times, inactivation of the phages was limited and the reported log removal was mainly associated with attachment of phages to the aquifer matrix. This study shows that attachment of MS2 and PRD1 is similar for oxic and suboxic sandy aquifers, and, therefore, setback distances used for sandy aquifers under oxic and suboxic conditions provide a similar level of safety. Sticking efficiency and the attachment rate coefficient, as measures for virus attachment, were evaluated as a function of the physico-chemical conditions.