Quantifying effects of soil heterogeneity on groundwater pollution at four sites in USA

Four sites located in the north-eastern region of the United States of America have been chosen to investigate the impacts of soil heterogeneity in the transport of solutes (bromide and chloride) through the vadose zone (the zone in the soil that lies below the root zone and above the permanent saturated groundwater). A recently proposed mathematical model based on the cumulative beta distribution has been deployed to compare and contrast the regions' heterogeneity from multiple sample percolation experiments. Significant differences in patterns of solute leaching were observed even over a small spatial scale, indicating that traditional sampling methods for solute transport, for example the gravity pan or suction lysimeters, or more recent inventions such as the multiple sample percolation systems may not be effective in estimating solute fluxes in soils when a significant degree of soil heterogeneity is present. Consequently, ignoring soil heterogeneity in solute transport studies will likely result in under- or overprediction of leached fluxes and potentially lead to serious pollution of soils and/or groundwater.The cumulative beta distribution technique is found to be a versatile and simple technique of gaining valuable information regarding soil heterogeneity effects on solute transport. It is also an excellent tool for guiding future decisions of experimental designs particularly in regard to the number of samples within one site and the number of sampling locations between sites required to obtain a representative estimate of field solute or drainage flux.     

Quantifying effects of soil heterogeneity on groundwater pollution at four sites in USA

Four sites located in the north-eastern region of the United States of America have been chosen to investigate the impacts of soil heterogeneity in the transport of solutes (bromide and chloride) through the vadose zone (the zone in the soil that lies below the root zone and above the permanent saturated groundwater). A recently proposed mathematical model based on the cumulative beta distribution has been deployed to compare and contrast the regions' heterogeneity from multiple sample percolation experiments. Significant differences in patterns of solute leaching were observed even over a small spatial scale, indicating that traditional sampling methods for solute transport, for example the gravity pan or suction lysimeters, or more recent inventions such as the multiple sample percolation systems may not be effective in estimating solute fluxes in soils when a significant degree of soil heterogeneity is present. Consequently, ignoring soil heterogeneity in solute transport studies will likely result in     

The effect of permafrost on stream biogeochemistry: A case study of two streams in the Alaskan (U.S.A.) taiga

Understanding interactions between permanently frozen soils and stream chemistry is important in predicting the effects of management, natural disturbance and changing permafrost distribution on stream ecosystems and nutrient budgets in subarctic watersheds. Chemical measurements of groundwater, soil water and stream water were made in two watersheds in the taiga of interior Alaska. One watershed (HiP) had extensive permafrost and the other (LoP) had limited permafrost. Soil water collected within the rooting zone (0.3--0.5 m) in both watersheds was high in dissolved organic carbon (DOC), dissolved organic nitrogen (DON) and dissolved inorganic nitrogen (DIN) but low in dissolved minerals (dominantly Ca, Mg and Na) and conductivity. The reverse was true for groundwater from springs and wells. Permafrost in the HiP basin prevented deep percolation of water and generated stormflows rich in DOC. The presence of permafrost in HiP resulted in higher fluxes of DOC, DON and DIN into stream water from upland soils.     

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.     

The influence of tidal marshes on upland groundwater discharge to estuaries

We investigated subsurface hydrology in two fringing tidal marshes and in underlying aquifers in the coastal plain of Virginia. Vertical distributions of hydraulic conductivity, hydraulic head and salinity were measured in each marsh and a nearby subtidal sediment. Discharge of hillslope groundwater into the base of the marshes and subtidal sediment was calculated using Darcy's law. In the marshes, fluxes of pore water across the sediment surface were measured or estimated by water balance methods. The vertical distribution of salt in shoreline sediments was modeled to assess transport and mixing conditions at depth. Hydraulic gradients were upward beneath shoreline sediments; indicating that groundwater was passing through marsh and subtidal deposits before reaching the estuary. Calculated discharge (6 to 10 liters per meter of shoreline per day) was small relative to fluxes of pore water across the marsh surface at those sites; even where discharge was maximal (at the upland border) it was 10 to 50 times less than infiltration into marsh soils. Pore water turnover in our marshes was therefore dominated by exchange with estuarine surface water. In contrast, new interstitial water entering subtidal sediments appeared to be primarily groundwater, discharged from below. The presence of fringing tidal marshes delayed transport and increased mixing of groundwater and solute as it traveled towards the estuaries. Soil-contact times of discharged groundwater were up to 100% longer in marshes than in subtidal shoreline sediments. Measured and modeled salinity profiles indicated that, prior to export to estuaries, the solutes of groundwater, marsh pore water and estuarine surface water were more thoroughly mixed in marsh soils compared to subtidal shoreline sediments. These findings suggest that transport of reactive solutes in groundwater may be strongly influenced by shoreline type. Longer soil-contact times in marshes provide greater opportunity for immobilization of excess nutrients by plants, microbes and by adsorption on sediment. Also, the greater dispersive mixing of groundwater and pore water in marshes should lead to increased availability of labile, dissolved organic carbon at depth which could in turn enhance microbial activity and increase the rate of denitrification in situations where groundwater nitrate is high.     

The influence of tidal marshes on upland groundwater discharge to estuaries

We investigated subsurface hydrology in two fringing tidal marshes and in underlying aquifers in the coastal plain of Virginia. Vertical distributions of hydraulic conductivity, hydraulic head and salinity were measured in each marsh and a nearby subtidal sediment. Discharge of hillslope groundwater into the base of the marshes and subtidal sediment was calculated using Darcy's law. In the marshes, fluxes of pore water across the sediment surface were measured or estimated by water balance methods. The vertical distribution of salt in shoreline sediments was modeled to assess transport and mixing conditions at depth. Hydraulic gradients were upward beneath shoreline sediments; indicating that groundwater was passing through marsh and subtidal deposits before reaching the estuary. Calculated discharge (6 to 10 liters per meter of shoreline per day) was small relative to fluxes of pore water across the marsh surface at those sites; even where discharge was maximal (at the upland border) it was 10 to 50 times less than infiltration into marsh soils. Pore water turnover in our marshes was therefore dominated by exchange with estuarine surface water. In contrast, new interstitial water entering subtidal sediments appeared to be primarily groundwater, discharged from below. The presence of fringing tidal marshes delayed transport and increased mixing of groundwater and solute as it traveled towards the estuaries. Soil-contact times of discharged groundwater were up to 100% longer in marshes than in subtidal shoreline sediments. Measured and modeled salinity profiles indicated that, prior to export to estuaries, the solutes of groundwater, marsh pore water and estuarine surface water were more thoroughly mixed in marsh soils compared to subtidal shoreline sediments. These findings suggest that transport of reactive solutes in groundwater may be strongly influenced by shoreline type. Longer soil-contact times in marshes provide greater opportunity for immobilization of excess nutrients by plants, microbes and by adsorption on sediment. Also, the greater dispersive mixing of groundwater and pore water in marshes should lead to increased availability of labile, dissolved organic carbon at depth which could in turn enhance microbial activity and increase the rate of denitrification in situations where groundwater nitrate is high.     

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

干旱和半干旱地区生物结皮的普遍发育显著改变了表层土壤的结构与养分富集特征,但其对土壤养分迁移和淋失的影响目前尚不明确。本研究针对黄土高原风沙土和黄绵土上发育的藓类生物结皮,以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倍。研究证实,生物结皮覆盖能够提高表层土壤对养分的吸附与固持能力,从而减少土壤表聚养分的深层渗漏和流失,对干旱和半干旱地区退化土壤的肥力提升与植被恢复具有重要意义。     

Transport of bacterial inoculants through intact cores of two different soils as affected by water percolation and the presence of wheat plants

Abstract Water flow-innduced transport of Burkholderia cepacia strain P2 and Pseudomonas fluorescens strain R2f cells through intact cores of loamy sand and silt loam field soils was measured for two percolation regimes, 0.9 and 4.4 mm h⁻¹, applied daily during 1 hour. For each strain, transport was generally similar between the two water regimes. Translocation of B. cepacia , with 4.4 mm h⁻¹, did occur initially in both soils. In the loamy sand soil, no change in the bacterial distribution occurred during the experiment (51 days). In the silt loam, B. cepacia cell numbers in the lower soil layers were significantly reduced, to levels at or below the limit of detection. Transport of P. fluorescens in both soils also occurred initially and was comparable to that of B. cepacia . Later in the experiment, P. fluorescens was not detectable in the lower soil layers of the loamy sand cores, due to a large decrease in surviving cell numbers. In the silt loam, the inoculant cell distribution did not change with time. Pre-incubation of the inoculated cores before starting percolation reduced B. cepacia inoculant transport in the loamy sand soil measured after 5 days, but not that determined after 54 days. Delayed percolation in the silt loam soil affected bacterial transport only after 54 days. The presence of growing wheat plants overall enhanced bacterial translocation as compared to that in unplanted soil cores, but only with percolating water. Percolation water from silt loam cores appeared the day after the onset of percolation and often contained inoculant bacteria. With loamy sand, percolation water appeared only 5 days after the start of percolation, and no inoculant bacteria were found. The results presented aid in predicting the fate of genetically manipulated bacteria in a field experiment.     

PHYTOREMEDIATION OF BTEX HYDROCARBONS: POTENTIAL IMPACTS OF DIURNAL GROUNDWATER FLUCTUATION ON MICROBIAL DEGRADATION

Volatile hydrocarbons have multiple potential fates in phytoremediation. This research investigated the relationship between biodegradation and plant uptake of BTEX compounds in laboratory and field settings. At a phytoremediation site, preliminary studies revealed minimal uptake into trees and enhanced degradation potential in the rhizosphere and in the bulk soil. Increased oxygen transport to the vadose zone caused by diurnal rise and fall of the water table was hypothesized to enhance degradation in the bulk soil. A detailed greenhouse study was then conducted to investigate potential bioremediation impacts using field-site soil and DN34 hybrid poplar trees.

In rhizosphere soils, the contaminated-planted reactor had significantly higher BTEX degrader populations versus the uncontaminated-planted reactor, as was anticipated. The bulk soil in the planted-contaminated reactor had increased degrader populations than the unplanted-contaminated soil or planted-uncontaminated soil, and planting increased degradation throughout the soil profile, not just in the limited volume of rhizosphere soils. Oxygen diffusive and advective transport into reactors was modeled and calculated. Oxygen input in planted reactors was at least 3 to 5 times higher than in unplanted reactors, and increasing oxygen input lead to increased degrader populations in a linear manner. These results combined with the knowledge that high-transpiration trees draw the contaminated groundwater to the capillary fringe and the rhizosphere indicate that phytoremediation can aid microbial degradation via multiple mechanisms: increasing degrader populations, increasing oxygen input via groundwater diurnal fluctuations, and transporting contaminants to the biologically-enriched soil profile.     

A Study of Vadose Zone Transport Model VLEACH

Application of vadose zone transport models has been hampered by lack of model validation. Difficulties to validate vadose zone models using field data not only come from model assumptions that are uncertain to the subsurface transport processes but also from the uncertainties associated with soil contaminants’ release time and quantity, soil sampling, sample transport, and analytical procedures. This article first conducts a test of a popularly used vadose zone transport VLEACH by comparing model results with a set of laboratory soil column infiltration and volatilization study data. The comparison shows a close agreement between the VLEACH model results and the laboratory data. Second, the sorption coefficient K_d calculated in VLEACH is compared with field data. The comparison indicates that VLEACH may overestimate the mass leached from soil to groundwater. The article also discusses the selection of the model simulation timestep, the vertical dimension increment, the Courant criterion, and the lower boundary condition using the sensitivity analysis method based on a case study of soil remediation for trichloroethylene. The procedures presented in this paper are important to practical model application and modification. This level of work should be routinely conducted for any new or modified version of vadose zone models.     

Optimization of a simple field method to determine mercury volatilization from soils—Examples of 13 sites in floodplain ecosystems at the Elbe River (Germany)

Mercury fluxes between soil and atmosphere have often been determined by using dynamic flux chambers and micrometeorological methods to assess ecological risks. However, both systems are complex, stationary, and expensive impeding measurements of Hg emissions at various field sites.The mobile, easy to handle, and cost-effective field method to determine total gaseous mercury (TGM), according to [Böhme, F., Rinklebe, J., Stärk, H.-J., Wennrich, R., Mothes, S., Neue, H.-U., 2005. A simple field method to determine mercury volatilisation from soils. Environ. Sci. Pollut. Res. (ESPR), 12: 133–135] creates a drop in air pressure that enhance the Hg emission. We optimized the sampling set-up using an air circulation system resulting in a continuous air flow over the soil surface. Thus, a drop in air pressure can be avoided and the detected TGM emissions are closer to reality. Additional benefits are an in-ground cylinder which inhibits lateral flow of gaseous mercury and the reduced size of the glass socket facilitating handling.To test the suitability of the optimized method, TGM emissions have been quantified on a set of Hg-contaminated riverine soils. Compared with non-polluted soils, mean Hg fluxes were strongly increased (between 138 and 711 ng m^?2 h^?1) and showed high spatial heterogeneity. Due to impacts of multiple environmental conditions that affect TGM emissions, no significant correlations have been found between Hg stocks in bulk soils and Hg fluxes.     

Effect of Soil Fabric on Transport of a LNAPL through Unsaturated Fine-Grained Soils: A Centrifugal Model Study

Centrifugal model tests were performed to study the impact of the fabric of a fine-grained soil on transport of a light non-aqueous phase liquid (LNAPL). An image processing technique was developed to extract contaminant transport and fate data from the centrifugal model. Two unconsolidated sites with different moisture contents and a saturated site consolidated due to self-weight were simulated using the centrifuge. The LNAPL migrated in the vertical direction as a narrow plume and formed a free product pool above the saturated zone in unsaturated and unconsolidated soils. However, the LNAPL migrated in the horizontal direction before moving in the vertical direction as a broad plume in the consolidated site. The test results showed that the final width of the plume in the unsaturated zone of the consolidated site was nearly two times as large as that for the unconsolidated sites. In addition, the rate of leak from the underground storage tanks (USTs) on consolidated soils was substantially higher when compared with those on the unconsolidated state. The comparison of LNAPL saturation profiles at the centerline of the centrifugal models during leakage showed that, depending on the soil fabric at a given time and depth, the LNAPL phase would be different; i.e., mobile or immobile (residual) in the same soil type. The test results provided additional insight into contribution of soil fabric on transport and fate of contaminants. The soil fabric controls the geological and hydro-geological properties of fine-grained soils and hence the contamination plume.     

Modeling jet fuel (JP-8) fate and transport in soils with plants

Vegetation is often used to clean up soils and groundwater contaminated with organic contaminants. Plant-induced upward water movement may draw organic contaminants spilled near the watertable to the more aerated near-surface soil. The objective of this study was to develop and verify a 1-D model of fate and transport of JP-8, a kerosene-based jet fuel, in soil. The modeling approach considered the advective and dispersive transport of jet fuels dissolved in groundwater, which may undergo simple first-order decay or linear adsorption. The governing partial differential advection dispersion equation was solved in one dimension. Data from an experiment of fate and transport of JP-8 with plant-induced upward water movement were used to verify the model. Simulated results with different scenarios described the experimental results well for different depths above the contaminated zone in both vegetated and unvegetated columns. Advection was the dominant mechanism near the contaminated zone and advection with retardation and decay was used to fit the data away from the contaminated zone. Results indicated that the soil water movement impacted the transport and concentration of JP-8 in the soil columns. This model can be used to simulate the fate of JP-8 associated with phytoremediation and evapotranspiration.     

A Comparison of Hydrocarbon Vapor Attenuation in the Field with Predictions from Vapor Diffusion Models

This study used field data from three sites in Southern California to evaluate vapor phase transport from: (1) free product (die-sel and gasoline spill) on groundwater; (2) dissolved benzene (gasoline spill) in groundwater; and (3) hydrocarbon-impacted soil (gasoline spill) in the vadose zone. A sampling program to evaluate the vapor pathway included the following: vertical profile data, minimal purging prior to sample collection, field analysis of data, confirmation of field data using a fixed laboratory analysis, and soil physical property data. Comparison of hydrocarbon vapor concentrations measured in this field study with those calculated using vapor diffusion models suggest that an additional attenuation factor of between 500 and 35,000 is needed to account for observed concentrations. Comparison of hydrocarbon profiles with oxygen, carbon dioxide, and methane values is consistent with the interpretation that biodegradation is primarily responsible for the observed attenuation. Therefore, vapor pathway models that do not account for bioattenuation will result in a large overesti-mation of the risk at spill sites and will not be consistent with field data.     

Leaching of Cd,Zn, Pb,and Cu in Packed and Undisturbed Columns of Soils Affected by the Spill from a Pyrite Mine in the South of Spain

A study has been made of the leaching of Cd, Zn, Pb, and Cu in three representative soils within the zone affected by the spill from a pyrite mine in Aznalcollar (Sevilla, Spain) employing packed soil columns. According to the breakthrough and cumulative leaching curves, the relative mobilities of the different toxic elements in the columns are as follows: Cd> Zn> Cu> Pb. The effect of leaching on the distribution of metals as a function of depth using intact soil cores was also studied. The results showed that the soils themselves have a good capacity for immobilizing the soluble fraction of the elements from the spilled mud. This capacity varied as follows: clayey soil with a high carbonate content > clayey soil with a moderate carbonate content > sandy-clay loam soil with a low carbonate content. However, sandy soils with a low carbonate content could pose a risk to groundwater if initial contamination was high. These results could be considered during the evaluation of remedial technologies for the immobilization of soil metals.     

Bioventing soils contaminated with petroleum hydrocarbons

Summary Bioventing combines the capabilities of soil venting and enhanced bioremediation to cost-effectively remove light and middle distillate hydrocarbons from vadose zone soils and the groundwater table. Soil venting removes the more volatile fuel components from unsaturated soil and promotes aerobic biodegradation by driving large volumes of air into the subsurface. In theory, air is several thousand times more effective than water in penetrating and aerating fuel-saturated and low permeability soil horizons. Aerobic microbial degradation can mitigate both residual and vapor phase hydrocarbon concentrations. Soil venting is being evaluated at a number of U.S. military sites contaminated with middle distillate fuels to determine its potential to stimulate in situ aerobic biodegradation and to develop techniques to promote in situ vapor phase degradation. In situ respirometric evaluations and field pilot studies at sites with varying soil conditions indicate that bioventing is a cost-effective method to treat soils contaminated with jet fuels and diesel.     

A field lysimeter study of heavy metal movement down the profile of soils with multiple metal pollution during chelate-enhanced phytoremediation

The agricultural soils near a copper smelter in southeast China were found to be highly contaminated with Cu, Pb, Zn, and Cd. Metal migration from the soil to groundwater presents an environmental risk that depends on the physicochemical properties of the contaminated soils. Soil solution samples were obtained using lysimeters from a loam soil with multiple metal pollutions over a period of about 1 yr. A field lysimeter study was also conducted to examine the potential use of (S, S')-ethylenediamine-N, N'-disuccinic acid trisodium salt (EDDSNa3) in chelate-enhanced phytoremedation and to evaluate the leaching of heavy metals. The average heavy metal concentrations in the soil solution (without the addition of EDDS) were high (e.g., 0.15 mg Pb L(-1) at a 50-cm depth) compared to the upper limit for protection of groundwater in China, but varied during the sampling period. Cu concentrations were not correlated with pH or dissolved organic carbon (DOC), but Zn and Cd concentrations were related to soil solution pH. EDDS enhanced metal solubility in the soil, but plant metal uptake by Elsholtzia splendens Nakai did not increase accordingly. There may be an increasing risk of groundwater pollution by Cu and the EDDS enhanced phytoremediation technique needs to be carefully applied to minimize this side effect.     

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.     

Local and Regional Effects on Community Structure of Dung Beetles in a Mainland-Island Scenario

Understanding the ecological mechanisms driving beta diversity is a major goal of community ecology. Metacommunity theory brings new ways of thinking about the structure of local communities, including processes occurring at different spatial scales. In addition to new theories, new methods have been developed which allow the partitioning of individual and shared contributions of environmental and spatial effects, as well as identification of species and sites that have importance in the generation of beta diversity along ecological gradients. We analyzed the spatial distribution of dung beetle communities in areas of Atlantic Forest in a mainland-island scenario in southern Brazil, with the objective of identifying the mechanisms driving composition, abundance and biomass at three spatial scales (mainland-island, areas and sites). We sampled 20 sites across four large areas, two on the mainland and two on the island. The distribution of our sampling sites was hierarchical and areas are isolated. We used standardized protocols to assess environmental heterogeneity and sample dung beetles. We used spatial eigenfunctions analysis to generate the spatial patterns of sampling points. Environmental heterogeneity showed strong variation among sites and a mild increase with increasing spatial scale. The analysis of diversity partitioning showed an increase in beta diversity with increasing spatial scale. Variation partitioning based on environmental and spatial variables suggests that environmental heterogeneity is the most important driver of beta diversity at the local scale. The spatial effects were significant only at larger spatial scales. Our study presents a case where environmental heterogeneity seems to be the main factor structuring communities at smaller scales, while spatial effects are more important at larger scales. The increase in beta diversity that occurs at larger scales seems to be the result of limitation in species dispersal ability due to habitat fragmentation and the presence of geographical barriers.     

VIRTUS, a model of virus transport in unsaturated soils.

As a result of the recently proposed mandatory groundwater disinfection requirements to inactivate viruses in potable water supplies, there has been increasing interest in virus fate and transport in the subsurface. Several models have been developed to predict the fate of viruses in groundwater, but few include transport in the unsaturated zone and all require a constant virus inactivation rate. These are serious limitations in the models, as it has been well documented that considerable virus removal occurs in the unsaturated zone and that the inactivation rate of viruses is dependent on environmental conditions. The purpose of this research was to develop a predictive model of virus fate and transport in unsaturated soils that allows the virus inactivation rate to vary on the basis of changes in soil temperature. The model was developed on the basis of the law of mass conservation of a contaminant in porous media and couples the flows of water, viruses, and heat through the soil. Model predictions were compared with measured data of virus transport in laboratory column studies and, with the exception of one point, were within the 95% confidence limits of the measured concentrations. The model should be a useful tool for anyone wishing to estimate the number of viruses entering groundwater after traveling through the soil from a contamination source. In addition, model simulations were performed to identify parameters that have a large effect on the results. This information can be used to help design experiments so that important variables are measured accurately.