Transport and transformation of sulfadiazine in soil columns packed with a silty loam and a loamy sand

Concerning the transport of the veterinary antibiotic sulfadiazine (SDZ) little is known about its possible degradation during transport. Also its sorption behaviour is not yet completely understood. We investigated the transport of SDZ in soil columns with a special emphasis on the detection of transformation products in the outflow of the soil columns and on modelling of the concentration distribution in the soil columns afterwards. We used disturbed soil columns near saturation, packed with a loamy sand and a silty loam. SDZ was applied as a 0.57 mg L(-1) solution at a constant flow rate of 0.25 cm h(-1) for 68 h. Breakthrough curves (BTC) of SDZ and its transformation products 4-(2-iminopyrimidin-1(2H)-yl)aniline and 4-hydroxy-SDZ were measured for both soils. For the silty loam we additionally measured a BTC for an unknown transformation product which we only detected in the outflow samples of this soil. After the leaching experiments the (14)C-concentration was quantified in different layers of the soil columns. The transformation rates were low with mean SDZ mass fractions in the outflow samples of 95% for the loamy sand compared to 97% for the silty loam. The formation of 4-(2-iminopyrimidin-1(2H)-yl)aniline appears to be light dependent and did probably not occur in the soils, but afterwards. In the soil columns most of the (14)C was found near the soil surface. The BTCs in both soils were described well by a model with one reversible (kinetic) and one irreversible sorption site. Sorption kinetics played a more prominent role than sorption capacity. The prediction of the (14)C -concentration profiles was improved by applying two empirical models other than first order to predict irreversible sorption, but also these models were not able to describe the (14)C concentration profiles correctly. Irreversible sorption of sulfadiazine still is not well understood.     

Sorption and desorption of sulfadiazine in soil and soil-manure systems

Sulfadiazine is a widely used veterinary medicine that has high potential to enter the environment, especially the soil compartment by the application of manure on agricultural land and grass land or by the deposition of dung pats on pasture. Once it reaches the soil environments, it may enter into surface and ground water. Therefore, sorption-desorption behavior of sulfadiazine was studied under laboratory conditions in five different soils varying in their physicochemical properties. A batch equilibration technique was used with initial aqueous solution concentration of sulfadiazine at 5, 0.5, 0.05, and 0.005 microg mL(-1). Sorption-desorption data in soils with and without manure were well fitted with Freundlich model in log form (r(2), 0.99). A sorption-desorption hysteresis effect was apparent in all soils. A significant amount of sulfadiazine was found tightly bound to the soil particles and did not desorb after the desorption process. Moreover, presence of manure enhanced hysteresis effect. Hysteresis coefficient (H) value from soils in absence of manure (0.9-1.0) increased to the soils in presence of manure (0.9-1.8). Soils in the absence of manure showed low level of K(D Sorp.) values ranging from 0.1 to 24.3, suggesting low level sorption of sulfadiazine with appreciable risk of run-off and leaching, and in turn, surface and ground water contamination. However, presence of manure increased the sorption tendency of sulfadiazine significantly (K(D Sorp.), 6.9-40.2). K(D) values pertaining to desorption cycle increased from 1.2-90.4 to 10.4-167.3 in absence and presence of manure, respectively.     

Second-order modeling of arsenite transport in soils

Rate limited processes including kinetic adsorption-desorption can greatly impact the fate and behavior of toxic arsenic compounds in heterogeneous soils. In this study, miscible displacement column experiments were carried out to investigate the extent of reactivity during transport of arsenite in soils. Arsenite breakthrough curves (BTCs) of Olivier and Windsor soils exhibited strong retardation with diffusive effluent fronts followed by slow release or tailing during leaching. Such behavior is indicative of the dominance of kinetic retention reactions for arsenite transport in the soil columns. Sharp decrease or increase in arsenite concentration in response to flow interruptions (stop-flow) further verified that non-equilibrium conditions are dominant. After some 40-60 pore volumes of continued leaching, 30-70% of the applied arsenite was retained by the soil in the columns. Furthermore, continued arsenite slow release for months was evident by the high levels of residual arsenite concentrations observed during leaching. In contrast, arsenite transport in a reference sand material exhibited no retention where complete mass recovery in the effluent solution was attained. A second-order model (SOM) which accounts for equilibrium, reversible, and irreversible retention mechanisms was utilized to describe arsenite transport results from the soil columns. Based on inverse and predictive modeling results, the SOM model successfully depicted arsenite BTCs from several soil columns. Based on inverse and predictive modeling results, a second-order model which accounts for kinetic reversible and irreversible reactions is recommended for describing arsenite transport in soils.     

Sorption kinetics and its effects on retention and leaching.

Sorption of pesticides to substrates used in biopurification systems is important as it controls the system's efficiency. Ideally, pesticide sorption should occur fast so that leaching of the pesticide in the biopurification system is minimized. Although modeling of pesticide transport commonly assumes equilibrium, this may not always be true in practice. Sorption kinetics have to be taken into account. This study investigated the batch sorption kinetics of linuron, isoproturon, metalaxyl, isoxaben and lenacil on substrates commonly used in a biopurification system, i.e. cow manure, straw, willow chopping, sandy loam soil, coconut chips, garden waste compost and peat mix. The first-order sorption kinetics model was fitted to the observed pesticide concentrations versus time resulting in an estimated kinetic rate constant alpha. Sorption appeared to be fast for the pesticides linuron and isoxaben, pesticides which were classified as immobile, while less mobile pesticides displayed an overall slower sorption. However, the substrate does not seem to be the main parameter influencing the sorption kinetics. Coconut chips, which is a substrate with a high organic matter content showed slow sorption for most of the pesticides. The effect of different estimated alpha values on the breakthrough of pesticides through a biopurification system was evaluated using the HYDRUS 1D model. Significant differences in leaching behavior were observed as a result of the obtained differences in sorption kinetics.     

Column- and catchment-scale transport of cadmium: effect of dissolved organic matter

In the Ellen catchment on the Pinjarra Plain, NE of Perth in Western Australia, cadmium from fertilisers is starting to leach from soils. About 70% of surface soils in the Ellen catchment are sandy and often on top of a shallow ephemeral water table. Adsorption of Cd in the sandy soils of the Ellen catchment was studied by batch adsorption and by leaching small columns of soil. Adsorption of Cd increases linearly with increasing soil organic matter content and exponentially with increasing pH. Cadmium is significantly mobilised in the sandy soils by dissolved organic matter.The capacity of most of the sandy soils in the Ellen catchment to adsorb phosphate from fertiliser has been saturated. Resulting concentrations in Ellen Brook average 500 μg L⁻¹ P. Cadmium is adsorbed more strongly in the sandy soils than phosphate and is just starting to leach into Ellen Brook. From a comparison of Cd/P ratios in water, soils and fertiliser, cadmium concentrations in Ellen Brook are estimated to be at 10–30% of their maximum for complete breakthrough from soils. Present concentrations of Cd in Ellen Brook average 0.1 μg L⁻¹ and are estimated to approach the maximum for complete breakthrough in 100 yr. Maximum Cd concentrations in Ellen Brook could range from 0.6 to 2 μg L⁻¹, depending on rates of input with fertiliser and future increases in agricultural land use in the catchment.Breakthrough curves, resulting from leaching Cd through small columns of sandy soil, indicate that adsorption significantly increases the effective hydrodynamic dispersion of Cd. Longitudinal dispersivities, measured at pore-water velocities of 0.7–14 m day⁻¹, were 5 cm for Cd and 0.1–0.2 cm for chloride. The much greater dispersion of Cd in the sandy soils than of chloride is shown not to be caused by non-equilibrium adsorption.     

Soil attenuation of As(III,V) and Se(IV,VI) seepage potential at ash disposal facilities

Leachate from ash landfills is frequently enriched with As and Se but their off-site movement is not well understood. The attenuation potential of As and Se by soils surrounding selected landfills during leachate seepage was investigated in laboratory column studies using simulated ash leachate. As(III, V) and Se(IV, VI) concentrations as well as pH, flow rate, and a tracer were monitored in influent and effluent for up to 800 pore volumes followed by sequential desorption, extraction, and digestion of column segments. Column breakthrough curves (BTCs) were compared to predictions based on previously measured sorption isotherms. Early As(V) breakthrough and retarded As(III) breakthrough relative to predicted BTCs are indicative of oxidative transformation during seepage. For Se(VI), which exhibits linear sorption and the lowest sorption propensity, measured BTCs were predicted fairly well by equilibrium sorption isotherms, except for the early arrival of Se(IV) in one site soil, which in part, may be due to higher column pH values compared to batch isotherms. Most of the As and Se retained by soils during leaching was found to be strongly sorbed (60–90%) or irreversibly bound (10–40%) with <5% readily desorbable. Redox potential favoring transformation to the more sorptive valence states of As(V) and Se(IV) will invoke additional attenuation beyond equilibrium sorption-based predictions. With the exception of Se(IV) on one site soil, results indicate that attenuation by down-gradient soils of As and Se in ash landfill seepage will often be no less than what is predicted by equilibrium sorption capacity with further attenuation expected due to favorable redox transformation processes, thus mitigating contaminant plumes and associated risks.     

Modelling of atrazine transport in the presence of surfactants

Laboratory experiments were conducted to examine the effect of detergents on transport of atrazine [2-chloro-4-(ethylamino)-6-(isopropylamino)-s-triazine] through loam and sandy loam soils under saturation conditions. The Convection Dispersion Equation (CDE) was used to model and quantify the effects of detergents on atrazine model parameters: the retardation factor (R), pore velocity (v) and dispersion coefficient (D). The transport parameters were estimated using moment technique and partition coefficient obtained from batch experiments and compared with best-fitted parameters, R and D, keeping pore velocity constant. Results indicated the CDE model was not successful in predicting atrazine transport in the presence of surfactants at high concentrations. In the case of anionic surfactant with Elora loam, the average predicted R and D from moment technique of 3.4 and 11.1 cm2/h, respectively were significantly different than fitted parameters (R = 39 and D = 227 cm2/h). The poor performance of CDE in the presence of surfactants results from physiochemical changes in herbicide solubility and retention to the soil matrix rather than changes in soil hydraulic properties since the predicted pore water velocities from moment technique were similar to those measured during leaching experiments. Nevertheless, BTC analysis with CDE showed that land application of anionic surfactant (sulphonic) significantly increased R and D and decrease v for both soils. Addition of sulphonic increased R of atrazine by 12 and 26 folds for loam and sandy loam soils, respectively. On the other hand non-ionic surfactants seemed to decrease R, especially in sandy loam soil, thus facilitating atrazine leaching through soil. Non-equilibrium conditions seemed to govern atrazine transport in the presence of surfactants; double peaks in breakthrough curves were observed, indicating a need for mathematical models to account for such phenomena. Atrazine dispersion and tailing seemed to be higher through Elora loam compared to Caledon sandy loam due to higher aggregation of the Elora soil.     

Tylosin sorption to silty clay loam soils, swine manure, and sand

The objectives of this study were to assess sorption and desorption of tylosin, a macrolide antimicrobial chemical used in swine, cattle, and poultry production, in three silty clay loam soils of South Dakota and compare soil sorption to sand and manure sorption. The silty clay loam soils, from a toposequence in eastern South Dakota, standardized sand samples, and swine manure were used in 24-h batch sorption studies with tylosin concentrations ranging from 25 to 232 micro mole/L. Desorption from soil was conducted over a four-day period. Partition coefficients, based on the Freundlich isotherm (K(f)) or K(d) values, were calculated. K(f) values for the silty clay loams were similar, not influenced by landscape position, and averaged 1,350 with isotherm slopes ranging from 0.85 to 0.93. K(f) values for sand were dependent on solution/sand ratios and pH, ranging from 1.4 to 25.1. K(d) values of manure were dependent on the solution type and ranged from 840 L/kg with urine to about 175 L/kg when sorbed from water. Desorption of tylosin from each soil over the four-day period was < 0.2% of the amount added. The soils' high K(f) values and low desorption amounts suggest that once tylosin is in these soils, leaching to lower depths may not occur. However, this does not preclude runoff with soil eroded particles. If tylosin reaches a sand aquifer, through bypass flow or other mechanism(s), movement in the aquifer most likely would occur.     

Evaluating equilibrium and non-equilibrium transport of bromide and isoproturon in disturbed and undisturbed soil columns

In this study, displacement experiments of isoproturon were conducted in disturbed and undisturbed columns of a silty clay loam soil under similar rainfall intensities. Solute transport occurred under saturated conditions in the undisturbed soil and under unsaturated conditions in the sieved soil because of a greater bulk density of the compacted undisturbed soil compared to the sieved soil. The objective of this work was to determine transport characteristics of isoproturon relative to bromide tracer. Triplicate column experiments were performed with sieved (structure partially destroyed to simulate conventional tillage) and undisturbed (structure preserved) soils. Bromide experimental breakthrough curves were analyzed using convective-dispersive and dual-permeability (DP) models (HYDRUS-1D). Isoproturon breakthrough curves (BTCs) were analyzed using the DP model that considered either chemical equilibrium or non-equilibrium transport. The DP model described the bromide elution curves of the sieved soil columns well, whereas it overestimated the tailing of the bromide BTCs of the undisturbed soil columns. A higher degree of physical non-equilibrium was found in the undisturbed soil, where 56% of total water was contained in the slow-flow matrix, compared to 26% in the sieved soil. Isoproturon BTCs were best described in both sieved and undisturbed soil columns using the DP model combined with the chemical non-equilibrium. Higher degradation rates were obtained in the transport experiments than in batch studies, for both soils. This was likely caused by hysteresis in sorption of isoproturon. However, it cannot be ruled out that higher degradation rates were due, at least in part, to the adopted first-order model. Results showed that for similar rainfall intensity, physical and chemical non-equilibrium were greater in the saturated undisturbed soil than in the unsaturated sieved soil. Results also suggested faster transport of isoproturon in the undisturbed soil due to higher preferential flow and lower fraction of equilibrium sorption sites.     

Behaviour of simazine in soil amended with the final residue of the olive-oil extraction process

Addition of organic wastes to agricultural soils is becoming a common practice as a disposal strategy and to improve the physical and chemical soil properties. However, in order to optimise the use of organic wastes as soil amendments, their effect on the behaviour of other compounds that are also used in agriculture, such as pesticides, needs to be assessed. In this work, we have investigated the effects of the addition of the final solid residue of the new technology of olive-oil extraction (extracted alperujo or solid olive-mill waste, SOMW2) on the sorption, degradation and leaching of the herbicide simazine in a sandy loam soil. The results are compared with those of a previous study, where the intermediary by-product of the olive-oil processing technology (unextracted alperujo or SOMW1) was applied to the same soil. The soil was amended in the laboratory with SOMW2 at two different rates (5% and 10% w/w). Simazine sorption isotherms showed a great increase in herbicide sorption after SOMW2 addition to soil. SOMW2 addition also increased sorption irreversibility. Incubation studies revealed reduced biodegradation of simazine in the soil amended with SOMW2 compared to the unamended soil. Breakthrough curves of simazine in handpacked soil columns showed that SOMW2 addition retarded the vertical movement of the herbicide through the soil and greatly reduced the amount of herbicide available for leaching. Interestingly, the results were quantitatively different from those obtained for the intermediary by-product SOMW1, illustrating the importance of the specific characteristics of the organic amendment in determining its effect on pesticide behaviour.     

Influence of surfactant sorption on the removal of phenanthrene from contaminated soils

Laboratory column flushing experiments were conducted to remove phenanthrene from contaminated soils by Triton X-100 (TX100) with an aim to investigating the effect of surfactant sorption on the performance of surfactant-enhanced remediation process. The effluent concentration of phenanthrene from soil columns showed strong dependence on the sorption breakthrough curves of TX100. The removal of phenanthrene from contaminated soils was enhanced only when the sorption breakthrough of TX100 occurred and the influent concentration of TX100 was greater than the critical enhanced flushing concentration (CEFC). The sorption of surfactant onto soils and the subsequent partitioning of contaminants into soil-sorbed surfactant had a significant effect on the solute equilibrium distribution coefficient (KD) and thus the flushing efficiency for phenanthrene. A model was developed to predict KD and CEFC values for simulating the performance of surfactant-enhanced flushing for contaminated soils. These results are of practical interest in developing effective and safe surfactant-enhanced remediation technologies.     

Influence of ageing of residues on the availability of herbicides for leaching

Losses by leaching of chlorotoluron, isoproturon and triasulfuron from small intact columns of a structured clay loam and an unstructured sandy loam soil were measured in five separate field experiments. In general, losses of all three herbicides were greater from the clay loam than from the sandy loam soil and the order between herbicides was always triasulfuron > isoproturon > chlorotoluron. Differences between experiments were also consistent for every soil/herbicide combination. There was no relationship between total loss and either total rainfall or cumulative leachate volume. When weighting factors were applied to the rainfall data to make early rainfall more important than later rainfall, there were significant positive relationships between cumulative weighted rainfall and total losses. Also, there were significant negative correlations between total losses and the delay to accumulation of 25 mm rainfall (equivalent to one pore volume of available water) in the different experiments. In laboratory incubations, there was a more rapid decline in aqueous (0.01 M calcium chloride) extractable residues than in total solvent extractable residues indicating increasing sorption with residence time. However, the rate of change in water extractable residues could not completely explain the decrease in leachability with ageing of residues in the field. Short-term sorption studies with aggregates of the two soils indicated slower sorption by those of the clay loam than by those of the sandy loam suggesting that diffusion into and out of aggregates may affect availability for leaching in the more structured soil. Small scale leaching studies with aggregates of the soils also demonstrated reductions in availability for leaching as residence time in soil was increased, which could not be explained by degradation. These results therefore indicate that time-dependent sorption processes are important in controlling pesticide movement in soils, although the data do not give a mechanistic explanation of the changes in leaching with ageing of residues.     

Transport modes and pathways of the strongly sorbing pesticides glyphosate and pendimethalin through structured drained soils

Leaching of the strongly sorbing pesticides glyphosate and pendimethalin was evaluated in an 8-month field study focussing on preferential flow and particle-facilitated transport, both of which may enhance the leaching of such pesticides in structured soils. Glyphosate mainly sorbs to mineral sorption sites, while pendimethalin mainly sorbs to organic sorption sites. The two pesticides were applied in equal dosage to a structured, tile-drained soil, and the concentration of the pesticides was then measured in drainage water sampled flow-proportionally.The leaching pattern of glyphosate resembled that of pendimethalin, suggesting that the leaching potential of pesticides sorbed to either the inorganic or organic soil fractions is high in structured soils. Both glyphosate and pendimethalin leached from the root zone, with the average concentration in the drainage water being 3.5 and 2.7 μg L⁻¹, respectively. Particle-facilitated transport (particles >0.24 μm) accounted for only a small proportion of the observed leaching (13-16% for glyphosate and 16-31% for pendimethalin). Drain-connected macropores located above or in the vicinity of the drains facilitated very rapid transport of pesticide to the drains. That the concentration of glyphosate and pendimethalin in the drainage water remained high (>0.1 μg L⁻¹) for up to 7 d after a precipitation event indicates that macropores between the drains connected to underlying fractures were able to transport strongly sorbing pesticides in the dissolved phase. Lateral transport of dissolved pesticide via such discontinuities implies that strongly sorbing pesticides such as glyphosate and pendimethalin could potentially be present in high concentrations (>0.1 μg L⁻¹) in both water originating from the drainage system and the shallow groundwater located at the depth of the drainage system.     

Adsorption and transport of arsenate in carbonate-rich soils: coupled effects of nonlinear and rate-limited sorption

The transport and fate of arsenate in carbonate-rich soil under alkaline conditions was investigated with multiple approaches combining batch, sequential extraction and column experiments as well as transport modeling studies. Batch experiments indicated that sorption isotherm was nonlinear over a wide range of concentration (0.1-200 mg L(-1)) examined. As(V) adsorption to the calcareous soil was initially fast but then continued at a slower rate, indicating the potential effect of rate-limited sorption on transport. Column experiments illustrated that transport of As(V) was significantly retarded compared to a non-reactive tracer. The degree of retardation decreased with increasing As(V) concentration. As(V) breakthrough curves exhibited nonideal transport behavior due to the coupled effects of nonlinear and rate-limited sorption on arsenate transport, which is consistent with the results of modeling studies. The contribution of nonlinear sorption to the arsenate retardation was negligible at low concentration but increased with increasing As(V) concentration. Sequential extraction results showed that nonspecifically sorbed (easily exchangeable, outer sphere complexes) fraction of arsenate is dominant with respect to the inner-sphere surface bound complexes of arsenate in the carbonate soil fraction, indicating high bioavailability and transport for arsenate in the carbonate-rich soils of which Fe and Al oxyhydroxide fractions are limited.     

Measurement and analysis of non-Fickian dispersion in heterogeneous porous media

Contaminant breakthrough behavior in a variety of heterogeneous porous media was measured in laboratory experiments, and evaluated in terms of both the classical advection-dispersion equation (ADE) and the continuous time random walk (CTRW) framework. Heterogeneity can give rise to non-Fickian transport patterns, which are distinguished by "anomalous" early arrival and late time tails in breakthrough curves. Experiments were conducted in two mid-scale laboratory flow cells packed with clean, sieved sand of specified grain sizes. Three sets of experiments were performed, using a "homogeneous" packing, a randomly heterogeneous packing using sand of two grain sizes, and an exponentially correlated structure using sand of three grain sizes. Concentrations of sodium chloride tracer were monitored at the inflow reservoir and measured at the outflow reservoir. Breakthrough curves were then analyzed by comparison to fitted solutions from the ADE and CTRW formulations. In all three systems, including the "homogeneous" one, subtle yet measurable differences between Fickian and non-Fickian transport were observed. Quantitative analysis demonstrated that the CTRW theory characterized the full shape of the breakthrough curves far more effectively than the ADE.     

Modeling depth-variant and domain-specific sorption and biodegradation in dual-permeability media

A dual-permeability model (S_1D_DUAL) was developed to simulate the transport of land-applied pesticides in macroporous media. In this model, one flow domain was represented by the bulk matrix and the other by the preferential flow domain (PFD) where water and chemicals move at faster rates. The model assumed the validity of Darcian flow and the advective-dispersive solute transport in each of the two domains with inter-domain transfer of water and solutes due to pressure and concentration gradients. It was conceptualized that sorption and biodegradation rates vary with soil depth as well as in each of the two flow domains. In addition to equilibrium sorption, kinetic sorption was simulated in the PFD. Simulations were conducted to evaluate the combined effects of preferential flow, depth- and domain-variant sorption, and degradation on leaching of two pesticides: one with strong sorption potential (trifluralin) and the other with weak sorption potential (atrazine). Simulation results for a test case showed that water flux in the PFD was three times more than in the matrix for selected storm events. When equilibrium sorption was considered, the simulated profile of trifluralin in each domain was similar; however, the atrazine profile was deeper in the PFD than in the bulk matrix under episodic storm events. With an assumption of negligible sorption in the PFD, both the atrazine and the trifluralin profiles moved twice deeper into the PFD. The simulated concentrations of the chemicals were several orders higher in the PFD than in the matrix, even at deeper depths. The volume fraction of the macropores and the sorption and biodegradation properties of the chemicals could also affect the amount of pesticides leaving the root zone. For an intense storm event, slow sorption reaction rates in the PFD produced higher breakthrough concentrations of atrazine at the bottom of the simulated soil profile, thus posing the risk for breakthrough of chemicals from the root zone.     

Multi-process herbicide transport in structured soil columns: experiments and model analysis

Model predictions of pesticide transport in structured soils are complicated by multiple processes acting concurrently. In this study, the hydraulic, physical, and chemical nonequilibrium (HNE, PNE, and CNE, respectively) processes governing herbicide transport under variably saturated flow conditions were studied. Bromide (Br-), isoproturon (IPU, 3-(4-isoprpylphenyl)-1,1-dimethylurea) and terbuthylazine (TER, N2-tert-butyl-6-chloro-N4-ethyl-1,3,5-triazine-2,4-diamine) were applied to two soil columns. An aggregated Ap soil column and a macroporous, aggregated Ah soil column were irrigated at a rate of 1 cm h(-1) for 3 h. Two more irrigations at the same rate and duration followed in weekly intervals. Nonlinear (Freundlich) equilibrium and two-site kinetic sorption parameters were determined for IPU and TER using batch experiments. The observed water flow and Br- transport were inversely simulated using mobile-immobile (MIM), dual-permeability (DPM), and combined triple-porosity (DP-MIM) numerical models implemented in HYDRUS-1D, with improving correspondence between empirical data and model results. Using the estimated HNE and PNE parameters together with batch-test derived equilibrium sorption parameters, the preferential breakthrough of the weakly adsorbed IPU in the Ah soil could be reasonably well predicted with the DPM approach, whereas leaching of the strongly adsorbed TER was predicted less well. The transport of IPU and TER through the aggregated Ap soil could be described consistently only when HNE, PNE, and CNE were simultaneously accounted for using the DPM. Inverse parameter estimation suggested that two-site kinetic sorption in inter-aggregate flow paths was reduced as compared to within aggregates, and that large values for the first-order degradation rate were an artifact caused by irreversible sorption. Overall, our results should be helpful to enhance the understanding and modeling of multi-process pesticide transport through structured soils during variably saturated water flow.     

Kinetic modeling of virus transport at the field scale

Bacteriophage removal by soil passage in two field studies was re-analyzed with the goal to investigate differences between one- and two-dimensional modeling approaches, differences between one- and two-site kinetic sorption models, and the role of heterogeneities in the soil properties. The first study involved removal of bacteriophages MS2 and PRDI by dune recharge, while the second study represented removal of MS2 by deep well injection. In both studies, removal was higher during the first meters of soil passage than thereafter. The software packages HYDRUS-ID and HYDRUS-2D, which simulate water flow and solute transport in one- and two-dimensional variably saturated porous media, respectively, were used. The two codes were modified by incorporating reversible adsorption to two types of kinetic sites. Tracer concentrations were used first to calibrate flow and transport parameters of both models before analyzing transport of bacteriophages. The one-dimensional one-site model did not fully describe the tails of the measured breakthrough curves of MS2 and PRD1 from the dune recharge study. While the one-dimensional one-site model predicted a sudden decrease in virus concentrations immediately after the peaks, measured data displayed much smoother decline and tailing. The one-dimensional two-site model simulated the overall behavior of the breakthrough curves very well. The two-dimensional one-site model predicted a more gradual decrease in virus concentrations after the peaks than the one-dimensional one-site model, but not as good as the one-dimensional two-site model. The dimensionality of the problem hence can partly explain the smooth decrease in concentration after peak breakthrough. The two-dimensional two-site model provided the best results. Values for k(att2) and k(det2) could not be determined at the last two of four monitoring wells, thus suggesting that either a second type of kinetic sites is present in the first few meters of dune passage and not beyond the second monitoring well, or that effects of soil heterogeneity and dimensionality of the problem overshadowed this process. Variations between single collector efficiencies were relatively small, whereas collision efficiencies varied greatly. This implies that the nonlinear removal of MS2 and PRD1 is mainly caused by variations in interactions between grain and virus surfaces rather than by physical heterogeneity of the porous medium. Similarly, a two-site model performed better than the one-site model in describing MS2 concentrations for the deep well injection study. However, the concentration data were too sparse in this study to have much confidence in the fitted parameters.     

Transport of silver nanoparticles (AgNPs) in soil

The effect of soil properties on the transport of silver nanoparticles (AgNPs) was studied in a set of laboratory column experiments, using different combinations of size fractions of a Mediterranean sandy clay soil. The AgNPs with average size of ∼30 nm yielded a stable suspension in water with zeta potential of −39 mV. Early breakthrough of AgNPs in soil was observed in column transport experiments. AgNPs were found to have high mobility in soil with outlet relative concentrations ranging from 30% to 70%, depending on experimental conditions. AgNP mobility through the column decreased when the fraction of smaller soil aggregates was larger. The early breakthrough pattern was not observed for AgNPs in pure quartz columns nor for bromide tracer in soil columns, suggesting that early breakthrough is related to the nature of AgNP transport in natural soils. Micro-CT and image analysis used to investigate structural features of the soil, suggest that soil aggregate size strongly affects AgNP transport in natural soil. The retention of AgNPs in the soil column was reduced when humic acid was added to the leaching solution, while a lower flow rate (Darcy velocity of 0.17 cm/min versus 0.66 cm/min) resulted in higher retention of AgNPs in the soil. When soil residual chloride was exchanged by nitrate prior to column experiments, significantly improved mobility of AgNPs was observed in the soil column. These findings point to the importance of AgNP-soil chemical interactions as a retention mechanism, and demonstrate the need to employ natural soils rather than glass beads or quartz in representative experimental investigations.     

Environmental transport of endogenous dairy manure estrogens

Although estrogens originating from dairy manure applied to agricultural soils as a fertilizer can potentially contaminate surface water and groundwater, the variables that control transport are poorly understood. Our objective was to assess the potential for off-site movement of endogenous dairy cattle estrogens when manure is applied on fields at agronomically relevant fertilization rates. Estrone (E1), 17α-estradiol (α-E2), and 17β-estradiol (β-E2) were used in laboratory sorption, desorption, and transformation incubations with both manure and an agriculturally relevant soil. Sorption on manure containing 44% organic carbon exceeded sorption on soil containing 0.8% organic carbon by 20 to 150 times, following the pattern of β-E2 > α-E2 > E1. Approximately 20% of E1 and 17% of α-E2 were desorbed from manure, whereas only about 4% of β-E2 was desorbed. Thirty to seventy percent of α-E2 and β-E2 were converted to E1 in soil and manure, making it imperative that transformation reactions be considered when predicting transport and potential biological effects in the environment. Overall results indicate that high organic carbon concentrations and relatively low amounts of desorption inhibit the potential for off-site transport of endogenous dairy manure estrogens.