Stochastic analysis of the relationships between saturated hydraulic conductivity variance and solute dispersion in heterogeneous soils

The effect of three‐dimensional heterogeneity of saturated hydraulic conductivity on the vertical transport of solutes in soils is examined by means of controlled numerical experiments. Saturated hydraulic conductivity, an important transport parameter that controls the dispersion of pollutants in heterogeneous soils, is assumed to be composed of a homogeneous mean value and a perturbation caused by the vertical variability of soil properties, producing a stochastic process in the mean flow direction. The spatial heterogeneity of porous soils is characterized by the variance and correlation scale of the saturated hydraulic conductivity in the transport domain. Numerical experiments are carried out to evaluate the extent of contaminant dispersion in Hawaiian Oxic soils when uncertainty exists as a result of the spatial heterogeneity of saturated hydraulic conductivity. Statistical analysis of the saturated hydraulic conductivity measurements on undisturbed soil cores from two locations in Hawaiian Oxic soils indicated two different soils with the same mean and different variances. The partial differential equations describing three‐dimensional transient flow and solute transport in soils with a random conductivity field were solved to evaluate the effect of these two variance levels on the transport of a contaminant plume originating from the surface. The significance of the variance on the spatial and temporal distribution of tracer concentrations is demonstrated using solute breakthrough curves at various depths in the soil profile. The longitudinal macrodispersivity resulting from tracer spreading in the heterogeneous soils with a finite local dispersivity is also analyzed. The analysis shows a similar solute dispersion behavior for the two variances. However, there is an overall reduction in the dispersion of solutes resulting from a uniform velocity field with the same mean. Macrodispersivity values in heterogeneous soils are proportional to the variance at smaller travel distances but converge to the same value at larger travel distances.     

黄土高原生物结皮覆盖对风沙土和黄绵土溶质运移的影响

干旱和半干旱地区生物结皮的普遍发育显著改变了表层土壤的结构与养分富集特征,但其对土壤养分迁移和淋失的影响目前尚不明确。本研究针对黄土高原风沙土和黄绵土上发育的藓类生物结皮,以Ca²⁺和Cl^-为示踪离子开展溶质穿透试验,对有无生物结皮层及其覆盖下不同深度土壤的溶质运移特征进行了研究。结果表明: 在0～5 cm土层,生物结皮覆盖延缓了风沙土和黄绵土的溶质穿透过程,其Cl^-的穿透时间比无结皮延长了3.83(风沙土)和2.09倍(黄绵土),而Ca²⁺则分别延长了2.50和2.73倍。生物结皮覆盖条件下,表层0～5 cm土壤溶质完全穿透所对应的孔隙体积数比下层5～10 cm土壤更高,且其穿透历时更长;其中,Cl^-的穿透时间分别增加了67.3%(风沙土)和51.8%(黄绵土),Ca²⁺的穿透时间分别增加了8.0%和33.7%。生物结皮覆盖降低了土壤孔隙水流速(37.5%～70.2%);除风沙土的5～10 cm土层外,生物结皮使溶质弥散系数提高了1.73～6.29倍,使溶质弥散度提高了2.77～20.95倍。置换液完全穿透土壤后,风沙土和黄绵土的生物结皮层Ca²⁺含量显著高于无结皮,其分别比无结皮高4.14和2.58倍。研究证实,生物结皮覆盖能够提高表层土壤对养分的吸附与固持能力,从而减少土壤表聚养分的深层渗漏和流失,对干旱和半干旱地区退化土壤的肥力提升与植被恢复具有重要意义。     

Leaching of Cryptosporidium parvum oocysts, Escherichia coli, and a Salmonella enterica serovar Typhimurium bacteriophage through intact soil cores following surface application and injection of slurry

Increasing amounts of livestock manure are being applied to agricultural soil, but it is unknown to what extent this may be associated with contamination of aquatic recipients and groundwater if microorganisms are transported through the soil under natural weather conditions. The objective of this study was therefore to evaluate how injection and surface application of pig slurry on intact sandy clay loam soil cores influenced the leaching of Salmonella enterica serovar Typhimurium bacteriophage 28B, Escherichia coli, and Cryptosporidium parvum oocysts. All three microbial tracers were detected in the leachate on day 1, and the highest relative concentration was detected on the fourth day (0.1 pore volume). Although the concentration of the phage 28B declined over time, the phage was still found in leachate at day 148. C. parvum oocysts and chloride had an additional rise in the relative concentration at a 0.5 pore volume, corresponding to the exchange of the total pore volume. The leaching of E. coli was delayed compared with that of the added microbial tracers, indicating a stronger attachment to slurry particles, but E. coli could be detected up to 3 months. Significantly enhanced leaching of phage 28B and oocysts by the injection method was seen, whereas leaching of the indigenous E. coli was not affected by the application method. Preferential flow was the primary transport vehicle, and the diameter of the fractures in the intact soil cores facilitated transport of all sizes of microbial tracers under natural weather conditions.     

Leaching of phage from Class B biosolids and potential transport through soil

The objective of this study was to investigate leaching and transport of viruses, specifically those of an indigenous coliphage host specific to Escherichia coli ATTC 15597 (i.e., MS-2), from a biosolid-soil matrix. Serial extractions of 2% and 7% (solids) class B biosolid matrices were performed to determine the number of phage present in the biosolids and to evaluate their general leaching potential. Significant concentrations of coliphage were removed from the biosolids for each sequential extraction, indicating that many phage remained associated with the solid phase. The fact that phage was associated with or attached to solid particles appeared to influence the potential for release and subsequent transport of phage under saturated-flow conditions, which was examined in a series of column experiments. The results indicated that less than 8% of the indigenous coliphage initially present in the biosolids leached out of the biosolid-soil matrix. A fraction of this was subsequently transported through the sandy porous medium with minimal retention. The minimal retention observed for the indigenous phage, once released from the biosolids, was consistent with the results of control experiments conducted to examine MS-2 transport through the porous medium.     

Controlled environment soil‐core microcosm unit for investigating fate,migration, and transformation of chemicals in soils

The controlled environment soil‐core microcosm unit (CESMU) methods embody a collection of techniques that began with soil sampling in the field and continued throughout the laboratory investigation of chemical fate, migration, and transformation in site‐specific soils; it was a cost‐effective investigative methodology that could be used to screen chemical materials before initiating high‐cost environmental field studies. Intact soil cores were collected in the field using a hydraulically controlled probe, delivering intact soil‐cores with minimal disturbance directly into high‐density polyethylene pipe (10.3‐cm ID). The inert polyethylene pipe was an effective hydrophobic barrier that remained an integral part of the soil‐core column, obviating subsequent transfers of soil. In the laboratory, each soil column was fitted with a porous ceramic plate and a polyethylene endcap containing fittings for teflon tubing, so that a tension could be applied at the bottom of each soil column (30–35 kPa) to mimic field conditions, thus preventing the undue buildup of water within columns that otherwise would change the chemical, physical, and biological properties of the soil. The intact soil‐cores were housed in the CESMU chamber, a controlled temperature unit with sufficient capacity for maintaining constant temperature within entire soil‐cores. Synthetic rain was added twice a week by peristaltic pump at rates simulating rainfall. Leachates were collected under tension via teflon tubing into flasks in darkness and kept at soil column temperature inside CESMU until harvested for analyses. Soil columns were harvested at intervals for sectioning by depth, extraction, and soil analyses. CESMU methods are applicable to investigations of water movement, soil chemistry, solute transport/transformation, and effects on plants.     

Estimating hydraulic conductivity of a sandy soil under different plant covers using minidisk infiltrometer and a dye tracer experiment

The objective of this study was to estimate the hydraulic conductivity of sandy soil under different plant cover at the locality Mláky II at Sekule (southwest Slovakia). Two sites were demarcated at the locality, with mainly moss species at glade site, and pine forest at forest site. The estimation of unsaturated hydraulic conductivity was conducted by (a) minidisk infiltrometer and (b) the analysis of a dye tracer total resident concentration. The latter approach assumed the applicability of the stochastic—convective flow theory in the sandy soil. In the dye tracer experiment, two plots (1 × 1 m each) were established in both sites, and 100 mm of dye tracer (Brilliant Blue FCF) solution was applied on the soil surface. Similar results were obtained in both plots, with more than 70 % area of horizons stained in the depth of 30–50 cm. In some cases, the predicted and measured hydraulic conductivity were found within an order of magnitude, thus revealing similar impact of different plant cover on hydraulic properties of sandy soil studied. In contrast to sandy soils used for agriculture, the influence of the plant/surface humus and topsoil interface extended in the form of a highly heterogeneous matrix flow to the depth of 50–60 cm, where it was dampened by horizontal layering.     

Measuring and modelling the transport and root uptake of chemicals in the unsaturated zone

To determine the mechanisms prescribing the movement and uptake of chemicals in the soil of the rootzone, controlled experiments were carried out in four lysimeters growing tomatoes. Each lysimeter had a depth-wise array of 9 Time Domain Reflectometry (TDR) probes to monitor the soil's water content. Chloride was used as an inert tracer, and was applied with the nutrient solution used for irrigation. Sulphate was used as a reactive tracer, and was applied as a pulse resident in the upper 100 mm of the soil. The measured water contents and the concentrations of the chemicals in the soil profile at the end of the experiment were compared to a deterministic model based on Richards' equation and the convection–dispersion equation linked with various macroscopic sink terms for root water and chemical uptake. The uptake function based on matric pressure head seems to describe the uptake of water and chemicals of our tomato plants best. At high soil solution concentration chloride and sulphate exclusion occurred. Our simple model could be used to describe the major features of coupled water and chemical uptake. However, our approach of inverse modelling to infer the parameters for solute transport and root uptake could not be used to distinguish between soil-based mechanisms and plant uptake mechanisms. The choice of the root water uptake model had only a small effect on the final water content profiles, but led to differences in the final solute profiles of sulphur and chloride. This indicates that tracers might provide improved determination of the uptake mechanisms.     

Impacts of Environmental Colloids on the Transport of 17 β-estradiol in Intact Soil Cores

Estrogens such as 17-β estradiol (E2) are endocrine-disrupting compounds and can affect the reproductive systems of aquatic organisms. Therefore, it is important to understand the mechanisms of their transport in the environment. E2 and its daughter product estrone (E1) are both strongly adsorbed by soil organic matter and have relatively short half-lives. Reduced contact time with soil makes transport of E2 and E1 in soil more likely. In this study, intact soil cores from three soils representing a range of particle size distribution, structure, and organic matter content were used to compare the transport of E2 with and without the presence of colloidal material fractionated from soil or swine manure. In chemical transport experiments conducted with undisturbed soil columns, E2 and E1 were measured both in solution and attached to suspended solids in column effluent. During the transport experiments, colloids carrying E2 passed through all soils, with the exception of the sandy soil. The presence of colloids decreased the first detection time of E2 in the aqueous phase, was correlated with greater peak E2 concentrations in the effluent of both loamy and clayey soils, but not through the sandy soil, and increased mass fractions of the E2 that was transported.     

Modeling Solute Transport by DLA in Soils of Northeastern Egypt

Arid soils in Egypt display large variability in solute transport properties, causing problems in soil management. To characterize this variability, dye infiltration experiments were conducted on four plots representing three main soil types in northeastern Egypt. The plots represented both cultivated and uncultivated land use. The observed dye patterns displayed a large variability and especially the clay soils indicated a high degree of preferential flow. The loamy sand and sandy soils displayed a more uniform dye distribution indicating more homogeneous soil properties. The observed dye patterns were modeled using a diffusion limited aggregation (DLA) model. The DLA is a random walk model where model parameters can be optimized using genetic algorithms (GA). The DLA model reproduced the observed dye patterns for all soils in an excellent way. The best fit was obtained with a specific combination of directional random walk probabilities P_u, P_d, P_r, and P_l for each plot (correlation 0.97–0.99). To account for soil layers with different hydraulic properties a two layer DLA model was developed. For all plots the P_u (upward random walk probability) was higher for the upper more homogeneous soil layer. The overall results showed that spatial variability resulting from solute transport for the investigated soils can be modeled using a DLA approach.     

Transduction of Escherichia coli in soil

Bacteriophage P1-mediated generalized transduction of Escherichia coli K-12 was assessed in nonsterile soil. Auxotrophic recipient cells (thr- leu- thi- rpsL) were incubated in a sandy and a silty clay loam soil, and the transducing phage lysates from prototrophic strains carrying transposon 10(Tn10) in either purE or aroL regions were added. At intervals, the bacterial populations derived from the soils were plated on selective-differential media to enumerate prototrophic (thr+, leu+, or Tcr) transductants. Of 100 bacterial isolates obtained on the selective-differential media, 58 (14 thr+; 11, leu+; 33 Tcr) were confirmed E. coli transductants. The frequency of transduction in soil was ca. 10(-6). These data demonstrate the potential use of bacteriophage P1 to genetically manipulate E. coli in situ.     

Bacterial transport in volcanic tuff cores under saturated flow conditions

An antibiotic‐resistant bacterium was tested for transport through volcanic tuff and sandstone cores. Tuff cores were representative of the geology of Rainier Mesa located on the Nevada Test Site (NTS). Rapid bacterial transport occurred in some of the tuff cores and all sandstone cores under the hydraulic heads used (5–500 cm). Hydraulic conductivity of the tuff cores ranged widely, 9.6 × 10^‐5 to 7.2 x 10^‐3 cm h^‐1. A much narrower range was observed for sandstone cores, 1.6 × 10^‐2 to 5.9 X 10^‐2 cm h^‐1, which served as experimental controls. The percentage of the initial bacterial inoculum recovered within 3 pore volumes from tuff and sandstone cores ranged from 9.4 to 54.7% and 0.20 to 2.9%, respectively. Bacterial recovery appeared to be controlled by the structure of the flow paths in rock cores and not by overall hydraulic conductivity. Saturated clay‐infiltrated and unfractured zeolitized tuff cores were impermeable to water flow, and therefore bacterial transport was not detected. Three routes of bacterial transport were discerned in permeable rock cores by comparison of the breakthrough patterns of bacteria and tracer solution (chloride ions) in cores of differing lithologies. In sandstone cores, where water flowed evenly through the matrix, bacteria were transported in a dispersed manner throughout the sandstone, whereas bacteria were transported primarily along preferred flow paths (fractures or macropores) in permeable tuff cores.     

森林生态系统林木根系对优先流的影响

土壤水和溶质运移是土壤学和环境科学研究的难点和热点,优先流是一种常见的土壤溶质运移形式,绕过土壤基质而优先运移至地下水源,造成土壤养分的流失和水质的恶化。林木根系是土壤层的重要部分,其结构形态影响着优先流过程,为量化林木根系结构对土壤优先流的影响,以首都圈森林生态系统鹫峰定位监测站为研究区域,利用野外染色示踪与室内分析相结合的方法,定量分析根长密度和根系生物量在优先流区和基质流区的变化。结果表明:1)随着土层深度的增加,根长密度表现为减小的趋势,对径级d1 mm,1d3 mm和3d5 mm根系而言,根长密度在优先流区大于基质流区发生概率分别为66.7%,88.9%和83.3%;2)根系d1 mm对优先流贡献度最大,均值为94.8%,1d3 mm和3d5 mm根系对优先流贡献度较小,均值分别为4.3%和0.9%;3)研究点根系生物量进行统计,66.7%优先流区根系生物量大于基质流区根系生物量。开展根系对优先流的影响研究,有助于探明土壤水分运移规律,分析地表地下水质恶化根源,为生态环境安全提供理论指导和技术支持。     

Validation of theoretical framework explaining active solute uptake in dynamically loaded porous media

Solute transport in biological tissues is a fundamental process necessary for cell metabolism. In connective soft tissues, such as articular cartilage, cells are embedded within a dense extracellular matrix that hinders the transport of solutes. However, according to a recent theoretical study (Mauck et al., 2003, J. Biomech. Eng. 125, 602–614), the convective motion of a dynamically loaded porous solid matrix can also impart momentum to solutes, pumping them into the tissue and giving rise to concentrations which exceed those achived under passive diffusion alone. In this study, the theoretical predictions of this model are verified against experimental measurements. The mechanical and transport properties of an agarose–dextran model system were characterized from independent measurements and substituted into the theory to predict solute uptake or desorption under dynamic mechanical loading for various agarose concentrations and dextran molecular weights, as well as different boundary and initial conditions. In every tested case, agreement was observed between experiments and theoretical predictions as assessed by coefficients of determination ranging from R²=0.61 to 0.95. These results provide strong support for the hypothesis that dynamic loading of a deformable porous tissue can produce active transport of solutes via a pumping mechanisms mediated by momentum exchange between the solute and solid matrix.     

Transport of iodide in structured soil under spring barley during irrigation experiment analyzed using dual-continuum model

Transport of radioactive iodide ¹³¹I^? in a black clay loam soil under spring barley in an early ontogenesis phase was monitored during controlled field irrigation experiment. It was found that iodide bound in the soil matrix could be mobilized by the surface leaching enhanced by mechanical impact of water drops and transported below the root zone of crops via soil cracks. The iodide transport through structured soil profile was simulated by the one-dimensional dual-continuum model, which assumes the existence of two inter-connected flow domains: the soil matrix domain and the preferential flow domain. The model predicted relatively deep percolation of iodide within a short time, in a good agreement with the observed vertical iodide distribution in soil. The dual-continuum approach proved to be an adequate tool for evaluation of field irrigation experiments conducted in structured soils.     

A new approach to model the spatiotemporal development of biofilm phase in porous media

Bacteria can exist within biofilms that are attached to the solid matrix of a porous medium. Under certain conditions, the biomass can fully occupy the pore space leading to reduced hydraulic conductivity and mass transport. Here, by treating biofilm as a growing, high-viscosity phase, a novel macroscopic approach to model biofilm spatial expansion and its corresponding effects on porous medium hydraulic properties is presented. The separate yet coupled flow of the water and biofilm phases is handled by using relative permeability curves that allow for biofilm movement within the porous medium and bioclogging effects. Fluid flow is governed by Darcy's law and component transport is set by the convection-diffusion equation reaction terms for each component. Here, the system of governing equations is solved by using a commercial multiphase flow reservoir simulator, which is used to validate the model against published laboratory experiments. A comparison of the model and experimental observations reveal that the model provides a reasonable means to predict biomass development in the porous medium. The results reveal that coupled flow of water and movement of biofilm, as described by relative permeability curves, is complex and has a large impact on the development of biomass and consequent bioclogging in the porous medium.     

Transport of natural soil nanoparticles in saturated porous media: effects of pH and ionic strength

To understand the effects of ionic strength and pH on the transport of natural soil nanoparticles (NS) in saturated porous media, aeolian sandy soil nanoparticles (AS), cultivated loessial soil nano particles (CS), manural loessial soil nanoparticles (MS) and red soil nanoparticles (RS) were leached with solutions of varying pH and ionic strength. The recovery rate of soil nanoparticles decreased in the order AS > RS > MS > CS. Transport of soil nanoparticles was enhanced with increasing pH and decreasing ionic strength and was attributable to changes in the Zeta potential of NS. Deposition of NS was also affected by the composition of soil nanoparticles and the surface charge. Column experiments showed that the interaction between soil nanoparticles and saturated quartz sand was mainly due to the physical and chemical properties of soil nanoparticles. The Derjaguin–Landau–Verwey–Overbeek interaction energies between NS and sand were affected by pHs and ionic strengths. Soil nanoparticles transport through saturated porous media could be accurately simulated by the one-dimensional advection-dispersion-reaction equation.     

Further Observations on Asymmetrical Solute Movement across Membranes

The permeability of frog skin under the influence of urea hyperosmolarity has been studied. Flux ratio asymmetry has been demonstrated again for tracer mannitol. The inhibitors DNP, CN^-, and ouabain have been used to eliminate active sodium transport and it was found that urea hyperosmolarity produces asymmetrical mannitol fluxes on frog skins having no short-circuit current. These findings suggest that flux ratio asymmetry is due to solute interaction and is unrelated to sodium transport. Studies with a synthetic membrane show clearly that bulk flow of fluid can produce a "solvent drag" effect and change flux ratios. When bulk flow is blocked and solute gradients allowed their full expression, then solute interaction "solute drag" is easily demonstrable in a synthetic system.     

Influence of precipitation and soil on transport of fecal enterococci in fractured limestone aquifers

Limestone aquifers provide the main drinking water resources of southern Italy. The groundwater is often contaminated by fecal bacteria because of the interaction between rocks having high permeability and microbial pollutants introduced into the environment by grazing and/or manure spreading. The microbial contamination of springwater in picnic areas located in high mountains can cause gastrointestinal illness. This study was carried out in order to analyze the interaction between Enterococcus faecalis and the soil of a limestone aquifer and to verify the influence of this interaction on the time dependence of groundwater contamination. E. faecalis was chosen because, in the study area involved, it represents a better indicator than Escherichia coli. The research was carried out through field (springwater monitoring) and laboratory experiments (column tests with intact soil blocks). The transport of bacterial cells through soil samples was analyzed by simulating an infiltration event that was monitored in the study area. Comparison of laboratory results with data acquired in the field showed that discontinuous precipitation caused an intermittent migration of microorganisms through the soil and produced, together with dispersion in the fractured medium (unsaturated and saturated zones), an articulated breakthrough at the spring. The short distances of bacterial transport in the study area produced a significant daily variability of bacterial contamination at the field scale.     

Simulated cadmium transport in macroporous soil during heavy rainstorm using dual-permeability approach

Numerical modelling is used to analyze the transport of cadmium in response to an extreme rainfall event. The cadmium transport through the soil profile was simulated by the one-dimensional dual-permeability model, which assumes the existence of two mutually communicating domains: the soil matrix domain and the preferential flow domain. The model is based on Richards’ equation for water flow and advection-dispersion equation for solute transport. A modified batch technique allowed us to consider domain specific sorption, i.e. each of the domains has its own distribution coefficient. The dual-permeability model predicts that the cadmium can be transported substantially below the root zone after the storm. On the other hand, classical single permeability approach predicted that almost all applied cadmium stays retained near the soil surface.