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.

@en

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.

@en

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.

@en

Ralstonia solanacearum is the causative agent of bacterial wilt of potato and other vegetable crops. Contaminated irrigation water contributes to the dissemination of this pathogen but the exact concentration or biological threshold to cause an infection is unknown. In two greenhouse experiments, potted potato plants (Solanum tuberosum) were exposed to a single irrigation with 50 mL water (non-invasive soil-soak inoculation) containing no or 102 – 108 CFU/mL R. solanacearum. The disease response of two cultivars, Kondor and HB, were compared. Disease development was monitored over a three-month period after which stems, roots and tubers of asymptomatic plants were analyzed for latent infections. First wilting symptoms were observed 15 days post inoculation in a plant inoculated with 5x109 CFU and a mean disease index was used to monitor disease development over time. An inoculum of 5x105 CFU per pot (1.3x102 CFU/g soil) was the minimum dose required to cause wilting symptoms, while one latent infection was detected at the lowest dose of 5x102 CFU per pot (0.13 CFU/g). In a second set of experiments, stem-inoculated potato plants grown in vitro were used to investigate the dose-response relationship under optimal conditions for pathogen growth and disease development. Plants were inoculated with doses between 0.5 and 5x105 CFU/plant which resulted in visible symptoms at all doses. The results led to a dose-response model describing the relationship between R. solanacearum exposure and probability of infection or illness of potato plants. Cultivar Kondor was more susceptible to brown-rot infections than HB in greenhouse experiments while there was no significant difference between the dose-response models of both cultivars in in vitro experiments. The ED50 for infection of cv Kondor was 1.1x107 CFU. Results can be used in management strategies aimed to reduce or eliminate the risk of bacterial wilt infection when using treated water in irrigation.@en

We investigated the transport and retention kinetics of silica encapsulated – silica core double stranded DNA particles (SiDNASi) through 15cm saturated quartz sand columns as a function of a wide range of colloid injection concentrations (C0 = 8.7×102 - 6.6×108 particles ml-1). The breakthrough curves (BTCs) exhibited an overall 2-log increase of maximum relative effluent concentration with increasing C0. Inverse curve fitting, using HYDRUS1D, demonstrated that a 1-site first order kinetic attachment (katt) and detachment (kdet) model sufficed to explain the C0-dependent SiDNASi retention behaviour. With increasing C0, katt log-linearly decreased, which could be expressed as an overall decrease in the single-collector removal efficiency (ƞ). The decrease in ƞ was likely due to increased electrostatic repulsion between aqueous phase- solid phase colloids, formation of shadow zones downstream of deposited colloids and removal of weakly attached colloids from the solid phase (quartz sand) attributing to increased aqueous phase-solid phase intercolloidal collisions as a function of increasing SiDNASi concentration. Our results implied, firstly, that the aqueous phase colloid concentration should be carefully considered in determining colloidal retention behaviour in saturated porous media. Secondly, colloidal transport and retention dynamics in column studies should not be compared without considering colloid influent concentration. Thirdly, our results implied that the applicability of SiDNASi as a conservative subsurface tracer was restricted, since transport distance and retention was colloid concentration dependent. However, the uniqueness of the DNA sequences in SiDNASi imparts the advantage of concurrent use of multiple SiDNASi for flow tracking or porous media characterization.@en

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.

@en

While the burden of disease from well-studied drinking water contaminants is declining, risks from emerging chemical and microbial contaminants arise because of social, technological, demographic and climatological developments. At present, emerging chemical and microbial drinking water contaminants are not assessed in a systematic way, but reactively and incidence based. Furthermore, they are assessed separately despite similar pollution sources. As a result, risks might be addressed ineffectively. Integrated risk assessment approaches are thus needed that elucidate the uncertainties in the risk evaluation of emerging drinking water contaminants, while considering risk assessors’ values. This study therefore aimed to (1) construct an assessment hierarchy for the integrated evaluation of the potential risks from emerging chemical and microbial contaminants in drinking water and (2) develop a decision support tool, based on the agreed assessment hierarchy, to quantify (uncertain) risk scores. A multi-actor approach was used to construct the assessment hierarchy, involving chemical and microbial risk assessors, drinking water experts and members of responsible authorities. The concept of value- focused thinking was applied to guide the problem-structuring and model-building process. The development of the decision support tool was done using Decisi-o-rama, an open-source Python library. With the developed decision support tool (uncertain) risk scores can be calculated for emerging chemical and microbial drinking water contaminants, which can be used for the evidence-based prioritisation of actions on emerging chemical and microbial drinking water risks. The decision support tool improves existing prioritisation approaches as it combines uncertain indicator levels with a multi-stakeholder approach and integrated the risk assessment of chemical and microbial contaminants. By applying the concept of value-focused thinking, this study addressed difficulties in evidence-based decision-making related to emerging drinking water contaminants. Suggestions to improve the model were made to guide future research in assisting policy makers to effectively protect public health from emerging drinking water risks. @en

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.

@en

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.

@en

Though coliform bacteria are used worldwide to indicate faecal pollution of groundwater, the parameters determining the transport of Escherichia coli in aquifers are relatively unknown. To investigate the effect of dissolved organic carbon (DOC) on the attachment of E. coli to saturated goethite-coated sand, we carried out column experiments with E. coli with and without humic acid (HA) in monovalent and divalent salt solutions. To characterize sorption of DOC and attachment of E. coli, we measured the pH of the influent and effluent, the cation concentrations and the zeta potential of particles. Depending on the chemistry of the E. coli suspension, the normalized breakthrough concentrations were over 80 times higher in columns treated with HA compared with columns not treated with HA. However, this difference was not constant: there were time-dependent variations in attachment of E. coli to the collector surface, and in the chemical composition of the bacterial suspension. Reduction in removal occurred because HA altered the surface charge of the collector and also sterically hindered E. coli. In addition, reduction of removal in a CaCl₂ bacterial suspension was probably caused by site-blocking mechanisms between HA and Ca²⁺ ions. Our results indicate that in the presence of DOC, the concept of geochemical heterogeneity in explaining attachment of biocolloids has limited relevance.

@en

Though coliform bacteria are used worldwide to indicate fecal pollution of groundwater, the parameters determining the transport of Escherichia coli in aquifers are relatively unknown. We evaluated the occurrence of both straining and attachment of E. coli ATCC25922 in columns of ultra-pure, angular, saturated quartz sand. The column experiments were conducted over a wide range of porous medium sizes, column heights, input concentrations, and pore water flow velocities. Straining and attachment were examined by modelling the breakthrough curves (with HYDRUS 1D). In addition, model output was compared with measured strained and attached bacteria via column extrusion experiments (in which sand was extruded from the column and placed in excess water) and flow reversal experiments (in which the pore water flow direction was reversed, thereby dislodging strained bacteria). Our model consisted of an attachment rate coefficient and a straining rate coefficient; both of these decreased with transport distance. The straining rate coefficient also decreased in a Langmuirian way, in response to the filling of available pore space, which in turn depended on influent bacteria concentration, quartz grain diameter, and transport distance. The maximum strained fraction was 25-30% of total bacteria mass applied to the column; the maximum attached fraction was 30-35%. The fit between modelled and measured (strained and attached) bacteria masses was acceptable, as was the sensitivity of the model output to fitted parameter values. Our results lead to a new description for the time-dependent mass balance of strained bacteria, which entails using three fitting parameters. The results also imply that column experiments in combination with retention profiles (or various column lengths) are not enough to explain the retention processes in a column. Column extrusion and flow reversal experiments provide vital additional information on the occurrence and magnitude of straining. Our straining model could be of assistance in evaluating the importance of straining and in incorporating the straining process in bacteria transport modelling.

@en

Because of heterogeneity among members of a bacteria population, deposition rates of bacteria may decrease upon the distance bacteria are transported in an aquifer. Such deposition rate decreases may result in retained bacteria concentrations, which decrease hyper-exponentially as a function of transport distance, and may therefore significantly affect the transport of colloids in aquifers. We investigated the occurrence of hyper-exponential deposition of Escherichia coli, an important indicator for fecal contamination, and the causes for such behavior. In a series of column experiments with glass beads of various sizes, we found that attachment of E. coli decreased hyper-exponentially, or, on logarithmic scale in a bimodal way, as a function of the transported distance from the column inlet. From data fitting of the retained bacteria concentration profiles, the sticking efficiency of 40% of the E. coli population was high (α=1), while the sticking efficiency of 60% was low (α=0.01). From the E. coli total population, an E. coli subpopulation consisting of slow attachers could be isolated by means of column passage. In subsequent column experiments this subpopulation attached less than the E. coli total population, consisting of both slow and fast attachers. We concluded that the main driver for the observed dual mode deposition was heterogeneity among members of the bacteria population. Intra-population may result in some microbes traveling surprisingly high distances in the subsurface. Extending the colloid filtration theory with intra-population variability may provide a valuable framework for assessing the transport of bacteria in aquifers.

@en