Comparison of nonlinear least squares and log transformation procedures for calculating volatilization coefficients

A nonlinear least squares procedure and a log transformation procedure for calculating first-order rate coefficients from experimental concentration-versus-time data were compared using laboratory measurements of the volatilization from water of 1,1,1-trichloroethane and 1,2-dichloroethane and the absorption of oxygen by water. Ratios of the nonlinear least squares to log transformation volatilization and absorption coefficients for 77 tests ranged from 0.955 to 1.08 and averaged 1.01. Comparison of the maximum, minimum, and mean root-mean-square errors of prediction for six sets of coefficients showed that the errors for the nonlinear least squares procedure were almost always smaller than the errors for the log transformation procedure.     

Seasonal air-water exchange fluxes of polychlorinated biphenyls in the Hudson River Estuary

Polychlorinated biphenyls (PCBs) were measured in the air and water over the Hudson River Estuary during six intensive field campaigns from December 1999 to April 2001. Over-water gas-phase SigmaPCB concentrations averaged 1100 pg/m3 and varied with temperature. Dissolved-phase SigmaPCB concentrations averaged 1100 pg/L and displayed no seasonal trend. Uncertainty analysis of the results suggests that PCBs with 5 or fewer chlorines exhibited net volatilization. The direction of net air/water exchange could not be determined for PCBs with 6 or more chlorines. Instantaneous net fluxes of SigmaPCBs ranged from +0.2 to +630 ng m(-2) d(-1). Annual fluxes of SigmaPCBs were predicted from modeled gas-phase concentrations, measured dissolved-phase concentrations, daily surface water temperatures and wind speeds. The net volatilization flux was +62 microg m(-2) yr(-1), corresponding to an annual loss of +28 kg/yr of SigmaPCBs from the Hudson River Estuary for the year of 2000.     

Washing of marine coastal sand in a batch reactor: sorption and desorption of BTEX

This study addresses the issues related to decontamination of marine beach sand accidentally contaminated by petroleum products. Sorption and desorption of BTEX (i.e., benzene, toluene, ethylbenzene, and xylene) onto the sand from Uran Beach, located near the city of Mumbai, India, were studied, and isotherms were determined using the bottle point method to estimate sorption coefficients. Alternatively, QSARs (i.e., quantitative structure activity relationships) were developed and used to estimate the sorption coefficients. Experiments for kinetics of volatilization as well as for kinetics of sorption and desorption in the presence of volatilization were conducted in a fabricated laboratory batch reactor. A mathematical model describing the fate of volatile hydrophobic organic pollutants like BTEX (via sorption and desorption in presence of volatilization) in a batch sediment-washing reactor was proposed. The experimental kinetic data were compared with the values predicted using the proposed models for sorption and desorption, and the optimum values of overall mass transfer coefficients for sorption (K(s)a(s)) and desorption (K(d)a(d)) were estimated. This was achieved by minimization of errors while using the sorption coefficients (Kp) obtained from either laboratory isotherm studies or the QSARs developed in the present study. Independent experimental data were also collected and used for calibration of the model for volatilization, and the values of the overall mass transfer coefficient for volatilization (K(g)a(g)) were estimated for BTEX. In these exercises of minimization of errors, comparable cumulative errors were obtained from the use of Kp values derived from experimental isotherms and QSARs.     

Volatility of TCDD and PCB from soil

Volatilisation fluxes of 2,3,7,8-tetrachlorodibenzodioxin and Aroclors 1242, 1254 and 1260 were estimated from the model of Jury et.al., developed for compounds with similar physical properties. The influence of soil water content, soil organic matter and water evaporation was investigated. Increasing soil organic matter decreased the volatilisation flux of all four chemicals. Water evaporation was accompanied by increased fluxes, implying that the chemicals were susceptible to wicking. The trends observed were in qualitative agreement with field and laboratory studies. The predicted fluxes serve as a basis for comparing their volatilisation behaviour against pesticides and other homologues, and for the construction of models describing the environmental fate of these compounds.     

Volatilization of 1,2-dibromo-3-chloropropane (DBCP) from soils.

The volatilization of DBCP from soils, as affected by the soil characteristics and application techniques, was studied in a laboratory experiment. The volatilization rate of DBCP applied in water was higher from sandy and silty loam soils than from clay soil. Water added after DBCP application acted as a soil cover, decreasing the volatilization rate. The results obtained with DBCP application in hexane to air-dry soils, indicate that adsorption could be an important factor in reducing the volatilization losses. Diffusion coefficients were calculated from the volatilization parameters, by using a simplified relationship between volatilization losses and diffusion through soil.     

Effect of water content on transient nonequilibrium NAPL-gas mass transfer during soil vapor extraction

The effect of water content on the volatilization of nonaqueous phase liquid (NAPL) in unsaturated soils was characterized by one-dimensional venting experiments conducted to evaluate the lumped mass transfer coefficient. An empirical correlation based upon the modified Sherwood number, Peclet number, and normalized mean grain size was used to estimate initial lumped mass transfer coefficients over a range of water content. The effects of water content on the soil vapor extraction SVE process have been investigated through experimentation and mathematical modeling. The experimental results indicated that a rate-limited NAPL-gas mass transfer occurred in water-wet soils. A severe mass transfer limitation was observed at 61.0% water saturation where the normalized effluent gas concentrations fell below 1.0 almost immediately, declined exponentially from the initiation of venting, and showed long tailing. This result was attributed to the reduction of interfacial area between the NAPL and mobile gas phases due to the increased water content. A transient mathematical model describing the change of the lumped mass transfer coefficient was used. Simulations showed that the nonequilibrium mass transfer process could be characterized by the exponent beta, a parameter which described the reduction of the specific area available for NAPL volatilization. The nonequilibrium mass transfer limitations were controlled by the soil mean grain size and pore gas velocity, were well described by beta values below 1.0 at low water saturation, and were well predicted with beta values greater than 1.0 at high water saturation.     

Volatilization of 1,2‐dibromo‐3‐chloropropane (DBCP) from soils

Abstract

The volatilization of DBCP from soils, as affected by the soil characteristics and application techniques, was studied in a laboratory experiment. The volatilization rate of DBCP applied in water was higher from sandy and silty loam soils than from clay soil. Water added after DBCP application acted as a soil cover, decreasing the volatilization rate. The results obtained with DBCP application in hexane to air‐dry soils, indicate that adsorption could be an important factor in reducing the volatilization losses.

Diffusion coefficients were calculated from the volatilization parameters, by using a simplified relationship between volatilization losses and diffusion through soil.     

PCB mass transfer coefficients determined by application of a water surface sampler

A water surface sampler (WSS) was employed in combination with greased surface deposition plates (GSDPs) to measure the particulate dry deposition and gas exchange of polychlorinated biphenyls (PCBs) in Chicago, IL. Vapor phase PCB fluxes were calculated by subtracting the particulate fluxes obtained from GSDPs from total (particulate+gas) fluxes obtained from the WSS. Vapor phase PCB fluxes were divided by ambient air concentrations measured with a high volume sampler to calculate overall gas phase PCB mass transfer coefficients (K(G)). The calculated average PCB MTC was 0.54+/-0.47 cm s(-1). This experimentally determined average gas phase overall mass transfer coefficient, K(G), agreed well with the ones reported from studies using similar techniques and agreed well with modeled values obtained using MTC correlations developed for the WSS.     

A pesticide emission model (PEM) Part II: model evaluation

In the first part of the paper, the development of a numerical pesticide emission model (PEM) is described for predicting the volatilization of pesticides applied to agricultural soils and crops through soil incorporation, surface spraying, or in the furrow at the time of planting. In this paper the results of three steps toward the evaluation of PEM are reported. The evaluation involves: (i) verifying the numerical algorithms and computer code through comparison of PEM simulations with an available analytical solution of the advection/diffusion equation for semi-volatile solutes in soil; (ii) comparing hourly heat, moisture and emission fluxes of trifluralin and triallate modeled by PEM with fluxes measured using the relaxed eddy-accumulation technique; and (iii) comparison of the PEM predictions of persistence half-life for 29 pesticides with the ranges of persistence found in the literature. The overall conclusion from this limited evaluation study is that PEM is a useful model for estimating the volatilization rates of pesticides from agricultural soils and crops. The lack of reliable estimates of chemical and photochemical degradation rates of pesticide on foliage, however, introduces large uncertainties in the estimates from any model of the volatilization of pesticide that impacts the canopy.     

Washing of Marine Coastal Sand in a Batch Reactor: Sorption and Desorption of BTEX

ABSTRACT

This study addresses the issues related to decontamination of marine beach sand accidentally contaminated by petroleum products. Sorption and desorption of BTEX (i.e., benzene, toluene, ethylbenzene, and xylene) onto the sand from Uran Beach, located near the city of Mumbai, India, were studied, and isotherms were determined using the bottle point method to estimate sorption coefficients. Alternatively, QSARs (i.e., quantitative structure activity relationships) were developed and used to estimate the sorption coefficients. Experiments for kinetics of volatilization as well as for kinetics of sorption and desorption in the presence of volatilization were conducted in a fabricated laboratory batch reactor. A mathematical model describing the fate of volatile hydrophobic organic pollutants like BTEX (via sorption and desorption in presence of volatilization) in a batch sediment-washing reactor was proposed. The experimental kinetic data were compared with the values predicted using the proposed models for sorption and desorption, and the optimum values of overall mass transfer coefficients for sorption (K_sa_s) and desorption (K_da_d) were estimated.This was achieved by minimization of errors while using the sorption coefficients (K_p) obtained from either laboratory isotherm studies or the QSARs developed in the present study. Independent experimental data were also collected and used for calibration of the model for volatilization,and the values of the overall mass transfer coefficient for volatilization (K_ga_g) were estimated for BTEX. In these exercises of minimization of errors, comparable cumulative errors were obtained from the use of K_p values derived from experimental isotherms and QSARs.     

Volatilisation of o-xylene from sandy soil

The diffusive release of o-xylene from two soils with different contents of organic carbon (1.1 % and 0.11 % TOC) and with two different water contents (app. 5 % w/w and 15 % w/w was studied in the laboratory. The soils were spiked with o-xylene in the laboratory. The fluxes were measured over a period of 24 hours. The measured fluxes were compared to predictions by two models. Model I, which is an analytical model, assumed instant local equilibrium between soil air, water and solids. The distribution coefficients were measured for the two soils, and Henry's constant and the diffusion coefficient in air were taken from the literature. This model overestimated the flux for o-xylene for all the tested combinations. The ratios between estimated and observed fluxes at 1 h were between 1.7 and 7.3. Model II assumed that the mass transfer of o-xylene between the solids and the water phase was kinetically controlled and was solved numerically. However, the predictions by the more advanced model were not significantly better than the prediction by the simple analytical model. The results indicate that prediction of o-xylene volatilisation from unsaturated soil is associated with substantial uncertainty.     

Volatilization of contaminants from suspended sediment in a water column during dredging

Remedial dredging of contaminated bed sediments in rivers and lakes results in the suspension of sediment solids in the water column, which can potentially be a source for evaporation of hydrophobic organic compounds (HOCs) associated with the sediment solids. Laboratory experiments were conducted in an oscillating grid chamber to simulate the suspension of contaminated sediments and flux to air from the surface of the water column. A contaminated field sediment from Indiana Harbor Canal (IHC) and a laboratory-inoculated University Lake (UL) sediment, Baton Rouge, LA, were used in the experiments, where water and solids concentration and particle size distribution were measured in addition to contaminant fluxes to air. A transient model that takes into account contaminant desorption from sediment to water and evaporation from the water column was used to simulate water and sediment concentrations and air fluxes from the solids suspension. In experiments with both sediments, the total suspended solids (TSS) concentration and the average particle diameter of the suspended solids decreased with time. As expected, the evaporative losses were higher for compounds with higher vapor pressure and lower hydrophobicity. For the laboratory-inoculated sediment (UL), the water concentrations and air fluxes were high initially and decreased steadily implying that contaminant release to the water column from the suspended solids was rapid, followed by evaporative decay. For the field sediments (IHC), the fluxes and water concentrations increased initially and subsequently decreased steadily. This implied that the initial desorption to water was slow and that perhaps the presence of oil and grease and aging influenced the contaminant release. Comparison of the model and experimental data suggested that a realistic determination of the TSS concentration that can be input into the model was the most critical parameter for predicting air emission rates.     

Coupled production and transport of selenium vapor in unsaturated soil: evaluation by experiments and numerical simulation

Volatilization of selenium (Se) from soil to the atmosphere involves several sequential chemical reactions that form volatile Se species, followed by transport of the gaseous Se through the soil. This paper describes a numerical model that simulates the chemical and physical processes governing the production and transport of Se vapor in unsaturated soil. The model couples the four Se species involved in the production of Se vapor through chemical reactions, and allows each to migrate through the soil by advection, liquid or vapor diffusion depending on its affinity for the dissolved or vapor phase. The coupled transformations and transport of the four Se species, i.e., selenate, selenite, elemental and organic Se, and Se vapor, were calculated based on the Crank-Nicolson finite difference method. The model was used to analyze fluxes of Se vapor measured from a soil amended with inorganic Se in the form of selenate and covered with unamended clean soil of various thicknesses. Evolution of Se vapor from the soil was very fast, with measurable amounts of Se detected within 24 h. The peak of Se volatilization, detected at the 6th day, reached 3.31 Se microgram/day for the uncovered soil, but was reduced to near the detection limit (0.05 microgram/day) in the presence of a 8- or 16-cm clean soil cover. With two reaction rate coefficients fitted to the data, the model described Se volatilization very well. The estimated rate coefficient of Se methylation was unexpectedly high, with a value of 0.167/day. The net volatilization of Se, however, was severely inhibited by the fast demethylation, i.e., the reverse reaction which converted volatile Se species back into nonvolatile forms. As a result, Se vapor only penetrated a few centimeters in the soil. The demethylation rate coefficient, assessed by independent transport experiments using dimethyl selenide, was estimated as 186.8/day, corresponding to a half-life of only 5.3 min for Se vapor. Results of this study indicated that rapid demethylation of Se vapor during its diffusive transport through a soil is probably an important limiting factor in the volatilization of Se under natural conditions.     

Volatilization of lindance from water in soil-free and flooded soil systems.

Volatilization of 14C-lindane from water in planchets and under flooded soil ecosystem was investigated. Lindane disappeared faster than parathion from planchets. More rapid loss of both insecticides occurred from water than from chloroform. Loss of lindane and parathion was related to measured losses of water by evaporation. During 5-day incubation under flooded soil conditions, disappearance of lindane was faster from open vials than from sealed vials, whereas in nonflooded soil, no volatile loss of the insecticide was evident despite water evaporation. Over 5 day incubation under flooded conditions, greater volatile loss of lindane occurred in sandy soil than in alluvial soil apparanetly due to greater adsorption to the soil colloids decreasing the insecticide concentration in the standing water on the laterite soil. Under identical conditions of water evaporation, lindane loss was directly proportional to its initial concentration in the water. These results suggest that considerable loss of soil applied pesticides can occur by volatilization from the standing water in flooded rice fields, particularly under tropical conditions.     

Simulating the temporal changes of OCP pollution in Hangzhou, China

A dynamic fugacity model was applied to simulate the changes of contents and transfer fluxes of hexachlorocyclohexane (HCHs) and dichloro-diphenyl-trichloroethane (DDTs) from 1950s in the environment of Hangzhou, China. The receptors are composed of air, surface water, soils, sediment and biota compartments. The model provides a method to combine loadings of HCHs and DDTs from various sources with a series of physical-chemical processes to estimate concentrations and transport fluxes of HCHs and DDTs. Model results suggested that the calculated concentrations were in line with the observed ones. The highest contents of HCH and DDT in the environment of study area were 523 t and 471 t before 1983, among which about 80.7% HCHs and 93.2% DDTs remained in the soil compartment. From 1984 to now, contents of HCHs and DDTs had decreased to about 0.07% and 0.40% of their highest amount (before 1983), and only about 0.001% and 0.014% will expect to be left in 2020 in the study area according to the model prediction. Before 1983, the main transfer fluxes of HCHs were deposition from air to soil, runoff from soil to water and diffusion from soil to air, but for DDTs the main transfer fluxes were deposition from air to soil and water, and transfer from water to sediment. From 1984 to now, runoff from soil to water and transfer from water to sediment became the dominant processes. Although a large amount of HCHs and DDTs had been applied to the study area, their residue levels in the soils were much lower than those in North China (had lesser HCHs and DDTs application than in South China) at present time, and close to other locations of South China (had similar HCHs and DDTs application level). It can be attributed to the high precipitation and temperature that enhances the processes of wet deposition, evaporation and degradation of OCPs. Sensitivities of the input parameters to the calculated concentrations were evaluated using coefficient-of-variation normalized sensitivity coefficients. The model was also subjected to uncertainty analyses using a Monte Carlo simulation.     

Modelling pesticide volatilization after soil application using the mechanistic model Volt'Air

Volatilization of pesticides participates in atmospheric contamination and affects environmental ecosystems including human welfare. Modelling at relevant time and spatial scales is needed to better understand the complex processes involved in pesticide volatilization. Volt'Air-Pesticides has been developed following a two-step procedure to study pesticide volatilization at the field scale and at a quarter time step. Firstly, Volt'Air-NH₃ was adapted by extending the initial transfer of solutes to pesticides and by adding specific calculations for physico-chemical equilibriums as well as for the degradation of pesticides in soil. Secondly, the model was evaluated in terms of 3 pesticides applied on bare soil (atrazine, alachlor, and trifluralin) which display a wide range of volatilization rates. A sensitivity analysis confirmed the relevance of tuning to K_h. Then, using Volt'Air-Pesticides, environmental conditions and emission fluxes of the pesticides were compared to fluxes measured under 2 environmental conditions. The model fairly well described water temporal dynamics, soil surface temperature, and energy budget. Overall, Volt'Air-Pesticides estimates of the order of magnitude of the volatilization flux of all three compounds were in good agreement with the field measurements. The model also satisfactorily simulated the decrease in the volatilization rate of the three pesticides during night-time as well as the decrease in the soil surface residue of trifluralin before and after incorporation. However, the timing of the maximum flux rate during the day was not correctly described, thought to be linked to an increased adsorption under dry soil conditions. Thanks to Volt'Air's capacity to deal with pedo-climatic conditions, several existing parameterizations describing adsorption as a function of soil water content could be tested. However, this point requires further investigation. Practically speaking, Volt'Air-Pesticides can be a useful tool to make decision about agricultural practices such as incorporation or for the estimation of overall pesticide volatilization rates, and it holds promise for time specific dynamics.     

Prediction of pesticide volatilization with PELMO 3.31

In the current EU risk assessment for pesticide registration, the European Community requires prediction of the concentration of each pesticide in air. A number of mathematical models are used to assess the fate of pesticides in groundwater, surface water and soil. PELMO 3.20 calculates the volatilization fluxes from bare soil and was improved in the new version PELMO 3.31 to include the effect of temperature and sorption in dry soil. The objective of this study was to evaluate the new version of PELMO 3.31 in predicting the pesticide volatilization under field conditions. Procymidone, malathion, and ethoprophos were the test compounds in two different seasons (autumn and winter). Comparing simulation results obtained with PELMO 3.31, after calibration, with the previous version PELMO 3.20 shows that the estimated volatilization results seems improved for malathion, similar or slightly overestimating in the warmer season for ethoprophos, and similar or slightly underestimating in the colder season for procymidone. The new release of PELMO allows a more accurate estimation of pesticides volatilization from soil as function of meteorological factors, especially for medium or low volatility pesticides. Some difficulties remain, such as the determination of the active air layer and the sorption increment with the soil drying.     

Pig slurry concentration by vacuum evaporation: influence of previous mesophilic anaerobic digestion process

Water can be removed from pig slurry by evaporation, through the application of wasted heat from a power plant or from other processes. Apart from obtaining a concentrate with an obviously higher nutrient concentration than the original slurry, another objective of water removal is to obtain water as condensate, which could be reused. The objective of this work was to study the vacuum evaporation of pig slurry liquid fraction and to evaluate condensate composition as a function of both pH (4, 5, and 6) and pig slurry type (fresh slurry and anaerobically digested slurry). Batch experiments showed that condensate characteristics, total ammonia nitrogen (NH3-N), volatile fatty acids (VFA), and chemical oxygen demand were strongly dependent on initial slurry pH. In addition to producing part of the required thermal energy, previous anaerobic digestion presented several other clear advantages. The consumption of VFA and other volatile organic compounds during anaerobic digestion reduced the volatilization of organic matter in the evaporation treatment and, consequently, provided a higher quality condensate.     

Volatilization of lindane from water in soil‐free and flooded soil systems

Abstract

Volatilization of ¹⁴C‐lindane from water in planchets and under flooded soil ecosystem was investigated. Lindane disappeared faster than parathion from planchets. More rapid loss of both insecticides occurred from water than from chloroform. Loss of lindane and parathion was related to measured losses of water by evaporation. During 5‐day incubation under flooded soil conditions, disappearance of lindane was faster from open vials than from sealed vials, whereas in nonflooded soil, no volatile loss of the insecticide was evident despite water evaporation. Over 5 day incubation under flooded conditions, greater volatile loss of lindane occurred in sandy soil than in alluvial soil apparently due to greater adsorption to the soil colloids decreasing the insecticide concentration in the standing water of the laterite soil. Under identical conditions of water evaporation, lindane loss was directly proportional to its initial concentration in the water. These results suggest that considerable loss of soil applied pesticides can occur by volatilization from the standing water in flooded rice fields, particularly under tropical conditions.     

Semi-empirical estimation of sorption of hydrophobic pollutants on natural sediments and soils

Methods were developed for estimating the equilibrium sorption behavior of hydrophobic pollutants. At low pollutant concentration (aqueous phase concentration less than half the solubility), sorption isotherms were linear, reversible, and characterized by a partition coefficient, K_p. Partition coefficients normalized to organic carbon, K_OC ($\text{K}{}_{\mathrm{OC}}\text{=}\text{K}{}_{\mathrm{p}}\text{fraction}$ organic carbon), were highly invariant over a set of sediments and soils collected from throughout the nation. Equations for estimating K_OC from water solubility (including crystal energy) and octanol/water partition coefficients were developed. The predictive equations were tested on literature sorption data and found to estimate measured K_OC's generally within a factor of two.