Predicting diffuse-source transfers of surfactants to surface waters using SWAT

Chemicals present in domestic wastewater can adsorb to solid phase materials during sewage treatment. If biosolids (or sewage sludge) are applied to land, these chemicals can be transferred to soil. Under some circumstances they can also be transferred to surface waters during storm events either in solution or attached to sediment. In this paper we describe the utility of the SWAT 2000 model to estimate diffuse-source surface water exposure to "down-the-drain" chemicals. The model was applied hypothetically to predict the behaviour of linear alkylbenzene sulphonate (LAS), an anionic surfactant commonly used in household detergents, in a small catchment in Bedfordshire, UK, where it has previously been successfully validated for stream flow, sediments and pesticides. LAS transfers were estimated for two scenarios: (1) realistic and (2) reasonable worst case, based on assumptions on sludge application rates and the concentration of LAS in sludge. In addition, the sensitivity of the model output to the proportion of the catchment to which sludge is applied was established. Soil wetness and the total quantity of biosolids applied were the biggest determinants of chemical transport from the catchment. The potential of SWAT as a higher-tier tool in environmental risk assessments is discussed.     

A screening method for preliminary assessment of risk to groundwater from land-applied chemicals

A simple mathematical model for initial screening is presented that can aid in evaluating the relative risk to groundwater from applying nonpolar synthetic organic chemicals to soil. The basic premise is that the magnitude of the quotient of the chemical concentration of the water entering the aquifer and the maximum allowable concentration (as established by EPA or Health Departments) represents the health risk of a chemical. The chemical concentration of the soil water is estimated based on conservative, simplifying assumptions and requires only readily available data such as: basic soil properties (organic matter and saturated hydraulic conductivity), organic chemical properties (octanol-water partition coefficient and degradation rate) and environmental factors (recharge rate and depth to groundwater).The methodology was applied to assess the relative risk of organic chemicals in municipal sewage sludge and pesticides applied to agricultural land. The results are realistic.     

One-at-a-time sensitivity analysis of pollutant loadings to subsurface properties for the assessment of soil and groundwater pollution potential

Chemical leak was numerically simulated for four chemical substances: benzene (light non-aqueous phase liquid (NAPL)), tetrachloroethylene (dense NAPL), phenol (soluble in water), and pentachlorophenol (white crystalline solid) in a hypothetical subsurface leak situation using a multiphase compositional transport model. One metric ton of chemical substances was assumed to leak at a point 3.51 m above the water table in a homogeneous unconfined aquifer which had the depth to water table of 7.135 m, the hydraulic gradient of 0.00097, the recharge rate of 0.7 mm/day, and the permeability of 2.92?×?10^?10 m². For comparison, surface spill scenarios, which had a long pathway from source to the water table, were simulated. Using the model results, point-source pollutant loadings to soil and groundwater were calculated by multiplying mass, impact area, and duration above and below the water table respectively. Their sensitivity to subsurface properties (depth to water table, recharge rate, porosity, organic carbon content, decay rate, hydraulic gradient, capillary pressure, relative permeability, permeability) was analyzed, with changing each parameter within acceptable ranges. The study result showed that the pollutant loading to groundwater was more sensitive to the subsurface properties than the pollutant loading to soil. Decay rate, groundwater depth, hydraulic gradient and porosity were influential to pollutant loadings. The impact of influential parameters on pollutant loadings was nonlinear. The dominant subsurface properties of pollution loadings (e.g., decay rate, groundwater depth, hydraulic gradient, and porosity for groundwater) also affect the vulnerability, and the subsurface pollutant loadings defined in this study are dependent on chemical properties as well, which indicates that the influential hydrogeological and physicochemical parameters to pollutant loadings can be used for pollution potential assessment. The contribution of this work is the suggestion that the sensitivity of pollutant loadings can be used for pollution potential assessment. Soil and groundwater pollution potential of chemicals are discussed altogether for leak scenarios. A physics-based model is used to understand the impact of subsurface properties on the fate and transport of chemicals above and below the water table, and consequently their impact on the pollutant loading to soil and groundwater.

     

Development and verification of a screening model for surface spreading of petroleum

Overflows and leakage from aboveground storage tanks and pipelines carrying crude oil and petroleum products occur frequently. The spilled hydrocarbons pose environmental threats by contaminating the surrounding soil and the underlying ground water. Predicting the fate and transport of these chemicals is required for environmental risk assessment and for remedial measure design. The present paper discusses the formulation and application of the Oil Surface Flow Screening Model (OILSFSM) for predicting the surface flow of oil by taking into account infiltration and evaporation. Surface flow is simulated using a semi-analytical model based on the lubrication theory approximation of viscous flow. Infiltration is simulated using a version of the Green and Ampt infiltration model, which is modified to account for oil properties. Evaporation of volatile compounds is simulated using a compositional model that accounts for the changes in the fraction of each compound in the spilled oil. The coupling between surface flow, infiltration and evaporation is achieved by incorporating the infiltration and evaporation fluxes into the global continuity equation of the spilled oil. The model was verified against numerical models for infiltration and analytical models for surface flow. The verification study demonstrates the applicability of the model.     

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.     

Non-steady state partitioning of dry cleaning surfactants between tetrachloroethylene (PCE) and water in porous media

Trapped organic solvents, in both the vadose zone and below the water table, are frequent sources of environmental contamination. A common source of organic solvent contamination is spills, leaks, and improper solvent disposal associated with dry cleaning processes. Dry cleaning solvents, such as tetrachloroethylene (PCE), are typically enhanced with the addition of surfactants to improve cleaning performance. The objective of this work was to examine the partitioning behavior of surfactants from PCE in contact with water. The relative rates of surfactants partitioning and PCE dissolution are important for modeling the behavior of waste PCE in the subsurface, in that they influence the interfacial tension of the PCE, and how (or if) interfacial tension changes over time in the subsurface. The work described here uses a flow-through system to examine simultaneous partitioning and PCE dissolution in a porous medium. Results indicate that both nonylphenol ethoxylate nonionic surfactants and a sulfosuccinate anionic surfactant partition out of residual PCE much more rapidly than the PCE dissolves, suggesting that in many cases interfacial tension changes caused by partitioning may influence infiltration and distribution of PCE in the subsurface. Non-steady-state partitioning is found to be well-described by a linear driving force model incorporating measured surfactant partition coefficients.     

Evaluation and comparison of multimedia mass balance models of chemical fate: application of EUSES and ChemCAN to 68 chemicals in Japan

The European Union System for Evaluation of Substances (EUSES) and the ChemCAN chemical fate model are applied to describe the fate of 68 chemicals on two spatial scales in Japan. Emission information on the chemicals has been obtained from Japan's Pollutant Release and Transfer Registry and available monitoring data gathered from government reports. Environmental concentrations calculated by the two models for the four primary environmental media of air, water, soil and sediment agree within a factor of 3 for over 70% of the data, and within a factor of 10 for over 87% of the data. Reasons for certain large discrepancies are discussed. Concentrations calculated by the models are generally consistent with the lower range of concentrations that are observed in the environment. Agreement between modeled and observed concentrations is considerably improved by including an estimate of the advective input of chemicals in air from outside Japan. The agreement between the EUSES and ChemCAN models suggests that results of individual chemical assessments are not likely to be significantly affected by the choice of chemical fate model. Primary sources of discrepancy between modeled and observed concentrations are believed to be uncertainties in emission rates, degradation half-lives, and the lack of data on advective inflow of contaminants in air.     

Environmental impact assessment (EIA): Equilibrium fugacity calculations and worst case approach

A mathematical model UNITRISK which can be used for screening purposes and risk ranking was set up to calculate relative risk values of single chemicals and mixtures of chemicals towards single or mixtures of organisms via contamination of air, water, soil and food dose. The concentration values are calculated by applying the fugacity concept or a dose must be defined. The dose-response is assumed to be linear versus concentration and degradative processes are not considered which is the worst case. The assumption that equilibria between the environmental phases exist is simplifying the model and is also representing the worst case. A mean risk factor is calculated which is 1 if the endpoint values (LC₅₀, LD₅₀, ADI, etc.) are exceeded for the investigated organisms and man.     

Impact of varying soil structure on transport processes in different diagnostic horizons of three soil types

When soil structure varies in different soil types and the horizons of these soil types, it has a significant impact on water flow and contaminant transport in soils. This paper focuses on the effect of soil structure variations on the transport of pesticides in the soil above the water table. Transport of a pesticide (chlorotoluron) initially applied on soil columns taken from various horizons of three different soil types (Haplic Luvisol, Greyic Phaeozem and Haplic Cambisol) was studied using two scenarios of ponding infiltration. The highest infiltration rate and pesticide mobility were observed for the Bt₁ horizon of Haplic Luvisol that exhibited a well-developed prismatic structure. The lowest infiltration rate was measured for the Bw horizon of Haplic Cambisol, which had a poorly developed soil structure and a low fraction of large capillary pores and gravitational pores. Water infiltration rates were reduced during the experiments by a soil structure breakdown, swelling of clay and/or air entrapped in soil samples. The largest soil structure breakdown and infiltration decrease was observed for the Ap horizon of Haplic Luvisol due to the low aggregate stability of the initially well-aggregated soil. Single-porosity and dual-permeability (with matrix and macropore domains) flow models in HYDRUS-1D were used to estimate soil hydraulic parameters via numerical inversion using data from the first infiltration experiment. A fraction of the macropore domain in the dual-permeability model was estimated using the micro-morphological images. Final soil hydraulic parameters determined using the single-porosity and dual-permeability models were subsequently used to optimize solute transport parameters. To improve numerical inversion results, the two-site sorption model was also applied. Although structural changes observed during the experiment affected water flow and solute transport, the dual-permeability model together with the two-site sorption model proved to be able to approximate experimental data.     

Developing a multimedia model of chemical dynamics in an urban area.

A multimedia model has been developed to account for the movement of semi-volatile organic compounds (SOCs) in an urban environment. The model, based on a Level III fugacity model of D. Mackay (Multimedia Environmental Models: The Fugacity Approach, Lewis Publishers, Boca Raton, FL, 1991), consists of six compartments: air, surface water, sediment, soil, vegetation, and an organic film that coats impervious surfaces. The latter is a newly identified compartment into which gas-phase SOCs partition and particle-phase SOCs are believed to be efficiently captured (M.L. Diamond, S.E. Gingrich, K. Fertuck, B.E. McCarry, G.A. Stern, B. Billeck, B. Grift, D. Brooker, T.D. Yager, Environ. Sci. Technol., 34 (2000a), 2900-2908). The model, parameterized for downtown Toronto, Ontario, Canada, and run with an illustrative emission rate for selected polycyclic aromatic hydrocarbons and homologues of polychlorinated dibenzodioxins, indicates that the film achieves the highest concentrations among media but that soils are the greatest sink for all but the least hydrophobic chemicals. The film "reflects" the more volatile chemicals into air, facilitates removal to surface waters by wash-off, and provides a surface on which photolytic degradation can occur. As such, the film is a transient sink that increases chemical mobility in urban areas by increasing air concentrations and the cycling of these compounds between air and urban surfaces and increasing water concentrations. Vegetation also accumulates SOCs, a portion of which is transferred to soil that reduces chemical mobility.     

Field calibration of surface: a model of agricultural chemicals in surface waters

Agricultural chemicals sporadically occur at detectable levels in the surface waters of intensively farmed watersheds. HSPF, a previously released model of agricultural chemicals in surface water, had been used to predict concentrations which were much higher (10 X) than those actually observed during monitoring studies. A new model, SURFACE, is described here which is much simpler than HSPF and gives better predictions of surface water concentrations. SURFACE uses PRZM, an EPA model, to calculate edge-of-field runoff losses and simple hydraulic routing algorithms to determine concentrations at the bottom of large river basins. In water systems sampled during 1985 and 1986, SURFACE predictions of annualized mean concentrations for alachlor, atrazine, cyanazine and metolachlor were within 0.09 ppb half of the time.     

The effect of export to the deep sea on the long-range transport potential of persistent organic pollutants

BACKGROUND: Export to the deep sea has been found to be a relevant pathway for highly hydrophobic chemicals. The objective of this study is to investigate the influence of this process on the potential for long-range transport (LRT) of such chemicals. METHODS: The spatial range as a measure of potential for LRT is calculated for seven PCB congeners with the multimedia fate and transport model ChemRange. Spatial ranges for cases with and without deep sea export are compared. RESULTS AND DISCUSSION: Export to the deep sea leads to increased transfer from the air to the surface ocean and, thereby, to lower spatial ranges for PCB congeners whose net deposition rate constant is similar to or greater than the atmospheric degradation rate constant. This is fulfilled for the PCB congeners 101, 153, 180, and 194. The spatial ranges of the congeners 8, 28, and 52, in contrast, are not affected by deep sea export. With export to the deep sea included in the model, the spatial ranges of the heavier congener are similar to those of the lighter ones, while the intermediate congeners 101 and 153 have the highest potential for long-range transport. CONCLUSIONS: Transfer to the deep ocean affects the mass balance and the potential for LRT of highly hydrophobic chemicals and should be included in multimedia fate models containing a compartment for ocean water.     

Electrical resistivity tomography as monitoring tool for unsaturated zone transport: an example of preferential transport of deicing chemicals

Non-invasive spatially resolved monitoring techniques may hold the key to observe heterogeneous flow and transport behavior of contaminants in soils. In this study, time-lapse electrical resistivity tomography (ERT) was employed during an infiltration experiment with deicing chemical in a small field lysimeter. Deicing chemicals like potassium formate, which frequently impact soils on airport sites, were infiltrated during snow melt. Chemical composition of seepage water and the electrical response was recorded over the spring period 2010. Time-lapse electrical resistivity tomographs are able to show the infiltration of the melt water loaded with ionic constituents of deicing chemicals and their degradation product hydrogen carbonate. The tomographs indicate early breakthrough behavior in parts of the profile. Groundtruthing with pore fluid conductivity and water content variations shows disagreement between expected and observed bulk conductivity. This was attributed to the different sampling volume of traditional methods and ERT due to a considerable fraction of immobile water in the soil. The results show that ERT can be used as a soil monitoring tool on airport sites if assisted by common soil monitoring techniques.     

Field calibration of surface: A model of agricultural chemicals in surface waters

Abstract

Agricultural chemicals sporadically occur at detectable levels in the surface waters of intensively farmed watersheds. HSPF, a previously released model of agricultural chemicals in surface water, had been used to predict concentrations which were much higher (10 X) than those actually observed during monitoring studies. A new model, SURFACE, is described here which is much simpler than HSPF and gives better predictions of surface water concentrations. SURFACE uses PRZM, an EPA model, to calculate edge‐of‐field runoff losses and simple hydraulic routing algorithms to determine concentrations at the bottom of large river basins. In water systems sampled during 1985 and 1986, SURFACE predictions of annualized mean concentrations for alachlor, atrazine, cyanazine and metolachlor were within 0.09 ppb half of the time.     

Modeling dissolution and volatilization of LNAPL sources migrating on the groundwater table

A vertically averaged two-dimensional model was developed to describe areal spreading and migration of light nonaqueous-phase liquids (LNAPLs) introduced into the subsurface by spills or leaks from underground storage tanks. The NAPL transport model was coupled with two-dimensional contaminant transport models to predict contamination of soil gas and groundwater resulting from a LNAPL migrating on the water table. Numerical solutions were obtained by using the finite-difference method. Simulations and sensitivity analyses were conducted with a LNAPL of pure benzene to study LNAPL migration and groundwater contamination. The model was applied to subsurface contamination by jet fuel. Results indicated that LNAPL migration were affected mostly by volatilization. The generation and movement of the dissolved plume was affected by the geology of the site and the free-product plume. Most of the spilled mass remained as a free LNAPL phase 20 years after the spill. The migration of LNAPL for such a long period resulted in the contamination of both groundwater and a large volume of soil.     

Numerical methods for describing chemical transport in the unsaturated zone of the subsurface

Finite-difference and finite-element methods of approximation have been extended to solve the one-dimensional nonlinear partial differential equations that describe the simultaneous transport of heat, moisture and chemical in the unsaturated zone. Especially for chemical transport, nodal spacing criteria are required to minimize numerical dispersion and oscillatory behavior in the solution vector for chemical concentration. Conservative criteria for nodal spacing for saturated flow can be used to set nodal spacing for unsaturated zone transport. When nodal spacing criteria are satisfied, for the same set of transport and boundary conditions, chemical concentration profiles calculated by the two numerical methods will be almost the same. A situation that is simulated very well with one-dimensional models, is the application of chemicals to land surfaces. To compare and contrast the characteristics of solutions given by the two numerical methods, moisture content, temperature and chemical concentration profiles for a 75-day period after application in the unsaturated zone are calculated for two representative types of organic chemicals. In the first, the chemical is very slowly degraded in the subsurface environment but strongly sorbed to soil surfaces. In the second, the chemical is rapidly degraded but weakly sorbed to soil surfaces. Because of differences in sorption coefficients and mechanisms of degradation, for the same set of hydrodynamic properties of the subsurface, the weakly sorbed chemical is more widely distributed throughout the unsaturated zone, whereas the strongly sorbed chemical stays very close to where it is put initially with little penetration into the subsurface. Satisfying nodal spacing criteria minimizes the impact of the method of approximation on the calculated solutions of the transport equations. For better model predictive performance, however, there are needs for more fundamental information on processes governing transport in the subsurface.     

Effects of temperature and water content on degradation of isoproturon in three soil profiles

The phenylurea herbicide isoproturon, 3-(4-isopropylphenyl)-1,1-dimethylurea (IPU), is widely used to control pre- and post-emergence of grass and broad-leaved weeds in cereal crops. Its degradation in soils is a key process for assessing its leaching risk to groundwater resources. The degradation properties of various samples from surface and subsurface soil (down to 1m depth) of a heterogeneous agricultural field were studied using (14)C-IPU. Laboratory incubations were carried out at 22 and 10 degrees C and at water contents 90% and 50% of the estimated water holding capacity (eWHC) corresponding to water potentials between -56 kPa and -660 MPa. Degradation was found to be more sensitive to water content variations than to temperature variations in the ranges that we used. For surface layers, at 10 and 22 degrees C, the degradation half-life increased by a factor 10 and 15, respectively, when water content decreased from 90% to 50% eWHC. Under optimal degradation conditions (i.e. 22 degrees C and 90% eWHC), 3-(4-isopropylphenyl)-1-methylurea (MDIPU) was the main metabolite in surface samples. At subsurface depths, IPU half-lives were larger than 100 d, IPU was the main compound after 92 d of incubation and the main metabolite was an unidentified polar metabolite. These results suggest a metabolic pathway involving hydroxylations for subsurface materials. IPU degradation was largely affected by water availability in both surface and subsurface horizons. Clay content seemed to play a major role in degradation processes in subsurface soil by determining through sorption IPU availability in soil solution and/or by limiting water availability for microorganisms.