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.     

杀虫剂十氯酮的多介质环境行为模拟

应用EQC模型模拟十氯酮在多介质环境中的归宿和迁移通量.结果表明:土壤是十氯酮最大的贮存库,在稳态平衡条件下,残留率达到95.0%;在稳态非平衡条件下,十氯酮单独排放到水体,有37.5%残留在于排放的水体中,其在大气的浓度水平和质量分布均很低,在沉积物中的质量则来自于水体向沉积物的沉降迁移;十氯酮主要通过水体的水平迁移和土壤的厌氧降解输出;十氯酮的主要界面迁移过程是大气向土壤的迁移,其次是水体向沉积物的沉降和大气向水体的迁移.     

PCB in soils and estimated soil-air exchange fluxes of selected PCB congeners in the south of Sweden

PCB concentrations were studied in different soils to determine the spatial variation over a region of approximately 11 000 km(2). PCB congener pattern was used to illustrate the spatial differences, as shown by principal component analysis (PCA). The relationship to different soil parameters was studied. PCB concentrations in soil showed a large variation between sampling-areas with median concentrations ranging between 2.3 and 332 ng g(-1) (dw). Highest concentrations were found at two sites with sandy soils, one with extremely high organic carbon content. Both sites were located on the west coast of southern Sweden. Soils with similar soil textures (i.e. sandy silt moraine) did not show any significant differences in PCB concentrations. PCB congener composition was shown to differ between sites, with congener patterns almost site-specific. PCB in air and precipitation was measured and the transfer of chemicals between the soil and air compartments was estimated. Soil-air fugacity quotient calculations showed that the PCBs in the soil consistently had a higher fugacity than the PCBs in the air, with a median quotient value of 2.7. The gaseous fluxes between soil and air were estimated using standard modelling equations and a net soil-air flux estimated by subtracting bulk deposition from gaseous soil-air fluxes. It was shown that inclusion of vertical sorbed phase transport of PCBs in the soil had a large effect on the direction of the net soil-air exchange fluxes.     

Simulation of long-term environmental dynamics of polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans using the dynamic multimedia environmental fate model and its implication to the time trend analysis of dioxins

Long term environmental fate of polychlorinated dibenzo-p-dioxins (PCDD) and polychlorinated dibenzofurans (PCDF) were studied using the fugacity-based dynamic multimedia environmental fate model. New formulation about soil wind erosion into air was included into the model. Effect of process assumption, sensitivity analysis for parameters, and tentative validation against the measured sediment core analysis was performed. Mass fluxes between compartments were estimated by the dynamic modeling. From the model estimates, major mass fluxes coming from emission source were shown to go to the soil and water through wet/dry deposition, then go to degradation mainly in the soil and sediment. Major mass fluxes of TCDD and OCDD come from the impurities in CNP (Chlornitrofen) and PCP (Pentachlorophenol) directly into the soil. Consideration about multimedia environmental dynamics using the modeled mass fluxes was shown in the discussion.     

Application of a level IV fugacity model to simulate the long-term fate of hexachlorocyclohexane isomers in the lower reach of Yellow River basin, China

A level IV multimedia fugacity model was established to simulate the fate and transfer of hexachlorocyclohexane (HCH) isomers in the lower reach of the Yellow River basin, China, during 1952-2010. The predicted concentrations of HCHs are in good agreement with the observed ones, as indicated by the residual errors being generally lower than 0.5 logarithmic units. The effects of extensive agricultural application and subsequent prohibition of HCHs are reflected by the temporal variation of HCHs predicted by the model. It is predicted that only 1.8 tons of HCHs will be left in 2010, less than 0.06% of the highest contents (in 1983) in the study area, and about 99% of HCHs remain in soil. The proportions of HCH isomers in the environment also changed with time due to their different physicochemical properties. Although beta-HCH is not the main component of the technical HCHs, it has become the most abundant isomer in the environment because of its persistence. The dominant transfer processes between the adjacent compartments were deposition from air to soil, air diffusion through the air-water interface and runoff from soil to water. Sensitivity analysis showed that degradation rate in soil, parameters related to major sources, and thickness of soils had the strongest influence on the model result. Results of Monte Carlo simulation indicated the overall uncertainty of model predictions, and the coefficients of variation of the estimated concentrations of HCHs in all the compartments ranged from 0.5 to 5.8.     

Dynamics of PAH deposition, cycling and storage in a mixed-deciduous (Quercus-Fraxinus) woodland ecosystem

Estimates of standing biomass and fluxes of biomass in a mixed-deciduous woodland were derived, and used with results for concentrations of seven polycyclic aromatic hydrocarbons (PAHs) in different compartments of the woodland system to quantitatively assess some of the key fluxes and burdens of PAHs in this complex system. We quantified PAH burdens in air, in leaves of three deciduous tree species, in leaf litter and in soil, and uptake of PAHs by the tree leaves; PAH fluxes in litterfall, and deposition to the litter layer on the woodland floor during winter were calculated from these data. Air burdens exhibited marked seasonal variations for all compounds, with lowest values in summer when combustion-related emissions were low. Leaves did not accumulate large burdens of PAHs while on the trees and consequently, litterfall-associated fluxes of PAHs were small, representing only a fraction of the burdens in the litter layer to which they were deposited. Higher PAH burdens in air in winter, combined with the organic-matter-rich nature of the litter layer, are thought to be responsible for fluxes of PAHs to the litter layer in winter being 20-170 times the peak litterfall fluxes. The soil compartment was calculated to contain 25 years' worth of deposition of benzo[ghi]perylene, the most recalcitrant PAH in this study. Storage quotients for fluoranthene, pyrene, benzo[k]fluoranthene and benzo[a]pyrene burdens in soil represented 7-10 years' worth of deposition, while fluorene and phenanthrene storage in soil approached unity with inputs (1 and 3 years' worth of deposition, respectively). The relative importance of storage and loss processes was therefore closely related to the physico-chemical properties of the PAH, and is discussed in relation to the cycling of carbon in the woodland.     

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.     

A multimedia environmental model of chemical distribution: fate, transport, and uncertainty analysis

This paper presented a framework for analysis of chemical concentration in the environment and evaluation of variance propagation within the model. This framework was illustrated through a case study of selected organic compounds of benzo[alpha]pyrene (BAP) and hexachlorobenzene (HCB) in the Great Lakes region. A multimedia environmental fate model was applied to perform stochastic simulations of chemical concentrations in various media. Both uncertainty in chemical properties and variability in hydrometeorological parameters were included in the Monte Carlo simulation, resulting in a distribution of concentrations in each medium. Parameters of compartmental dimensions, densities, emissions, and background concentrations were assumed to be constant in this study. The predicted concentrations in air, surface water and sediment were compared to reported data for validation purpose. Based on rank correlations, a sensitivity analysis was conducted to determine the influence of individual input parameters on the output variance for concentration in each environmental medium and for the basin-wide total mass inventory. Results of model validation indicated that the model predictions were in reasonable agreement with spatial distribution patterns, among the five lake basins, of reported data in the literature. For the chemical and environmental parameters given in this study, parameters associated to air-ground partitioning (such as moisture in surface soil, vapor pressure, and deposition velocity) and chemical distribution in soil solid (such as organic carbon partition coefficient and organic carbon content in root-zone soil) were targeted to reduce the uncertainty in basin-wide mass inventory. This results of sensitivity analysis in this study also indicated that the model sensitivity to an input parameter might be affected by the magnitudes of input parameters defined by the parameter settings in the simulation scenario. Therefore, uncertainty and sensitivity analyses for environmental fate models was suggested to be conducted after the model output was validated based on an appropriate input parameter settings.     

The effects of atmospheric nitrogen deposition on the soil chemistry of coniferous forests in The Netherlands

Nitrogen fluxes, particularly those of ammonium, are extremely high in Dutch forests. In soils exposed to high ammonium deposition, acidification, eutrophication or a combination of both processes may occur. In addition to the amounts of ammonium deposited, the rate of soil nitrification determines which process takes place. A nation-wide investigation, in which three coniferous tree species were involved, was carried out to study the relation between deposition fluxes, measured by means of throughfall and bulk samplers, and the chemical composition of the soil. The ammonium deposition accounted directly for the high ammonium content and the high ammonium/cation ratios in the soil. In the top layer of most of the forest soils which were investigated nitrification rates were low. In these stands ammonium/cation ratios in the soil often reflected ammonium/cation ratios in throughfall water. Even in soils with relatively high nitrification rates, ammonium concentrations exceeded those of nitrate in the top layer of the mineral soil, indicating that ammonium deposition was more important than nitrification rate in determining the predominant form of nitrogen.     

Residues of organochlorine pesticides in Alabama soils

A survey was made of 36 Alabama agricultural soils to assess residues of formerly used organochlorine pesticides. Compounds determined comprised alpha- and gamma-hexachlorocyclohexane, heptachlor, heptachlor-exo-epoxide, trans- and cis-chlordane, trans-nonachlor, dieldrin, toxaphene, DDT and DDE. Concentrations varied by several orders of magnitude among farms and appeared to be log-normally distributed. Highest concentrations (ng g(-1) dry soil, arithmetic means) were found for toxaphene (285+/-390) and DDTs (p,p'-DDE, 22.7+/-21.4; p,p'-DDT, 24.6+/-30.5; o,p'-DDT, 4.00+/-5.86; p,p'-DDD, 2.40+/-2.41) which were once heavily used in the southern USA. Pesticide residues were not proportional to soil organic carbon content indicating that residue concentrations were a reflection of pesticide application history and dissipation rates rather than air-soil equilibrium. Mean ratios of DDT/DDE in six regions of the state ranged from 0.39 to 1.5, and compound ratios for chlordanes and toxaphene were different from those in the technical mixtures.     

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.     

Polycyclic aromatic hydrocarbons in lake sediments from the High Tatras

European alpine lake systems are used as indicators of air quality over the continent. Preliminary data showed high polycyclic aromatic hydrocarbons (PAH) loads in the High Tatras (Eastern Europe) in comparison to other mountain regions. Here, insight on the spatial distribution of PAH is provided from analysis of top-core sediments of 27 alpine lakes distributed along the High Tatras.Top-core sediment concentrations were higher than those in deep-cores, and they were higher than those observed in other European high mountain regions. The PAH profiles were uniform and comparable to those observed in aerosols and snow, indicating that atmospheric deposition was the predominant PAH input pathway to the lakes. Good agreement between estimated atmospheric deposition and sedimentation fluxes was observed. However, in several lakes in the western range higher sediment fluxes may correspond to higher PAH depositions levels. The higher concentrations may also reflect inputs from potential emission source areas.     

Analysis of atmospheric behavior of PCDDs/PCDFs by a one-compartment box model

The overall atmospheric behavior of PCDDs/PCDFs in the Kanto region, Japan, was simulated by a one-compartment box model. For each homologue the relative significance and temperature dependences of dry deposition, wet deposition, degradation, and advection in both gas and particulate phases were examined and compared. The results of the model calculation suggested that the rates for dry deposition are comparable to those for wet deposition, and the rates for advection are comparable to those for bulk (dry+wet) depositions in the Kanto region. On the other hand, the rates of degradation for PCDDs/PCDFs in the atmosphere in the Kanto region would be negligible. The emission rates and the bulk deposition fluxes in the entire Kanto region estimated by the model calculation based on observed air concentrations were 0.084-0.90 kg-TEQ/month and 0.045-0.43 kg-TEQ/month, respectively. These estimated emission rates and bulk deposition fluxes were slightly higher than the estimated emission rate based on observed emission concentrations and the estimated bulk deposition fluxes based on observed deposition fluxes collected on water deposition surface, respectively. This study showed the model calculation can be available for understanding of the overall atmospheric behavior, verification of the source inventory, and estimation of deposition flux on the actual environment including various deposition surfaces.     

A GIS-based multimedia watershed model: development and application

A multimedia model was developed using publicly available geographical information system (GIS) data, chemical release information and local monitoring networks to assess the fate of trichloroethene (TCE) within the Passaic River Watershed. Seven environmental media, air, water, sediment, surface soil, terrestrial vegetation, root zone soil and vadose zone soil, were modeled in this study along with their sub-compartments. The Passaic River Watershed is described using the NJDEP geographical information system (GIS) resources, the United States Geological Survey (USGS) and the United States Soil Conservation Services (US SCS) soil data. The introduction of spatial resolution to a multimedia, unsteady state model is performed in this work, and represents an important step in expanding the use of equilibrium models to provide far reaching information on the fate of toxic contaminants within a given environmental unit. The spatial representation of cross-boundary fluxes was successfully demonstrated with the use of sub-watershed as an environmental unit and the direct assessment of TCE for each of the 11 sub-watersheds that make up the Passaic River Basin in northern New Jersey. Important data gaps identified during the development of this model include the lack of comprehensive monitoring data on organic contaminants, and non-uniformity among available physical environmental data from different government agencies.     

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

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

Multimedia fate of petroleum hydrocarbons in the soil: Oil matrix of constructed biopiles

A dynamic multimedia fugacity model was used to evaluate the partitioning and fate of petroleum hydrocarbon fractions and aromatic indicator compounds within the soil: oil matrix of three biopiles. Each biopile was characterised by four compartments: air, water, soil solids and non-aqueous phase liquid (NAPL). Equilibrium partitioning in biopile A and B suggested that most fractions resided in the NAPL, with the exception of the aromatic fraction with an equivalent carbon number from 5 to 7 (EC₅₋₇). In Biopile C, which had the highest soil organic carbon content (13%), the soil solids were the most important compartment for both light aliphatic fractions (EC₅₋₆ and EC₆₋₈) and aromatic fractions, excluding the EC₁₆₋₂₁ and EC₂₁₋₃₅. Our starting hypothesis was that hydrocarbons do not degrade within the NAPL. This was supported by the agreement between predicted and measured hydrocarbon concentrations in Biopile B when the degradation rate constant in NAPL was set to zero. In all scenarios, biodegradation in soil was predicted as the dominant removal process for all fractions, except for the aliphatic EC₅₋₆ which was predominantly lost via volatilization. The absence of an explicit NAPL phase in the model yielded a similar prediction of total petroleum hydrocarbon (TPH) behaviour; however the predicted concentrations in the air and water phases were significantly increased with consequent changes in potential mobility. Further comparisons between predictions and measured data, particularly concentrations in the soil mobile phases, are required to ascertain the true value of including an explicit NAPL in models of this kind.     

Degradation half-life times of PCDDs, PCDFs and PCBs for environmental fate modeling

Literature search of the knowledge on the degradation of persistent organic pollutants (POPs) in environmental compartments air, water, soil and sediment was done in purpose to find properties of POPs of interest for modeling. One degradation process, hydrolysis (chemical degradation), was omitted as negligibly slow for POPs studied. The other two, photolysis and biodegradation processes, were considered separately in purpose to develop estimation procedures. The estimates can be given as pseudo first-order rate constants kP for photolysis and kB for biodegradation. For each compartment, an overall degradation rate is k(tot) = kP + kB and lifetime t(1/2) = ln 2/k(tot). The latter values, lifetimes in each compartment, will be used as input parameters to the Baltic Sea model.     

Estimation of PCDD/F distribution and fluxes in the Venice Lagoon,Italy: combining measurement and modelling approaches

The available experimental information on the occurrence of PCDD/Fs in the Venice Lagoon, Italy, was compiled and used to calculate fugacities for the environmental compartments of sediment, suspended particulate matter (SPM), water and air and then used to estimate fugacity ratios and assess the likely net direction of flux between media. The bottom sediment: SPM fugacity ratios for different PCDD/Fs indicate conditions close to equilibrium, suggestive of the close coupling of SPM with re-suspended sediment. Sediment/water and the sediment/air fugacity ratios suggest that net flux directions vary depending on the congener and the location within the lagoon. Net sediment-water-air movement (i.e. re-mobilisation/volatilisation) is suggested for the lighter congeners from the industrial canals, where the highest PCDD/F concentrations in the lagoon occur. The tendency to volatilise increases with decreasing congener molecular weight. In contrast, net deposition (air-water-sediment) appears to be occurring for the heaviest (hepta- and octa-) substituted PCDD/Fs. OCDF represents a marker of the industrial district of the lagoon, decreasing in concentration and as a fraction of total PCDD/Fs with increasing distance. The fugacity-based quantitative water air sediment interaction (QWASI) mass-balance model was applied to the central part of the lagoon. The key parameters for the determination of the model output, identified by a sensitivity analysis, were: the sediment active depth, the sediment re-suspension and deposition rates, and the total input of PCDD/Fs to the system. The QWASI model also indicates the tendency for the lighter PCDD/Fs to be released from surface sediment to the water column.     

Validation of modeling approach to evaluate congener-specific concentrations of polychlorinated dibenzo-p-dioxins and dibenzofurans in air and soil near a solid waste incinerator

A concise modeling approach using long-term averaged meteorological data was developed to estimate site-specific concentrations of congeners of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) near a solid waste incinerator. This approach consists of calculation of atmospheric dispersion, dry and wet deposition of gaseous and particle-bound congeners, and non-steady-state concentrations in soil. The predictability of this approach was evaluated by comparison of calculated concentrations of congeners in soil with those measured at eight locations near a municipal solid waste incinerator (MSWI). The variation of these concentrations due to variability of meteorological parameters is small. A considerable number of mean values show good agreement with measured concentrations within a factor of three. The reasonable agreement between calculated and measured concentrations indicates that algorithms for the calculation of vapor-phase deposition and non-steady-state concentrations in soil must be included in the modeling approach for an accurate estimation of the concentrations of congeners of PCDD/Fs emitted from MSWIs to the atmosphere. For a detailed estimation of site-specific concentrations, it is important to specify the bulk density of soil in the evaluated area, together with meteorological parameters.