Modelling field-scale cadmium transport below the root zone of a sewage sludge amended soil in an arid region in Central Iran

Addition of trace metals such as cadmium to soils in metal-rich sewage sludge may result in contamination of soil and groundwater. This study addresses the plot-scale transport of Cd derived from sewage sludge in a layered clay soil in an arid region of central Iran. Sewage sludge was enriched by Cd at rates of 38 and 80 mg kg(-1) and applied to experimental soil plots using a complete random block design with three replicates. Cadmium concentration was measured as a function of depth after 185 and 617 days. HYDRUS-1D and MACRO codes were calibrated for Cd transport in the site treated with 80 mg kg(-1) sewage sludge. Model parameters were estimated by inverse modelling using the SUFI-2 procedure. The site treated with 38 mg kg(-1) cadmium was used to test the calibrated models. Both convection-dispersion equation (CDE) and non-equilibrium CDE in HYDRUS-1D produced reasonable calibration results. However, the estimated Freundlich sorption constants were significantly smaller than those measured in a batch study. A site tracer experiment revealed the existence of substantial macropore flow. For this reason we applied MACRO to account for this process. The calibration and test results with MACRO were as good as those obtained by HYDRUS-1D with the difference that adsorption constants were much closer to the measured ones. This indicates that in HYDRUS-1D, the adsorption parameters were underestimated in order to allow a deeper transport of Cd which had actually occurred due to macropore flow. A 20-year simulation scenario depicting the long-term effect of sludge application indicated small risk of groundwater contamination. However, high concentration of Cd near the soil surface raises a concern about the crop Cd uptake which should be further investigated.     

Ambient concentrations and source apportionment of PCBs and trace elements around an industrial area in Izmir, Turkey

Atmospheric concentrations of polychlorinated biphenyls (PCBs) and trace elements were measured at two sites (Industrial and Urban) located around the Aliaga industrial region, Izmir, Turkey. Average sigma 36PCB concentrations were 3136+/-824 and 1371+/-642 pg m(-3) for summer and winter periods, respectively in the Industrial site and they were 314+/-129 and 847+/-610 pg m(-3) in the Urban site. The elemental content showed that the PM(10) measured at the Industrial site was dominated by terrestrial elements and trace elements emitted by the iron-steel plants (Fe, Zn and Pb). The elemental profile at the Urban site was typical for Aegean Region that was dominated by terrestrial elements (Ca, Al, Mg) and sea salt (Na). Sources of particle-phase PCBs and trace elements were identified using factor analysis (FA) and were apportioned by chemical mass balance (CMB) model. FA suggested that the steel industry, fuel oil combustion, or the nearby vinyl chloride process in the petrochemical plant, and soil were significant PCB sources. CMB results showed that at the Industrial site, the contribution of steel industry and soil to particle-phase PCBs were 71% and 22%, respectively, while at the Urban site, the contributions were 33% and 49%, respectively. Steel industry was also the dominant contributor for trace elements around the site. Fugacity calculations in air and soil showed that the soil acts as a secondary source to the atmosphere for low molecular weight PCBs especially in summer and as a sink for the higher molecular weight ones.     

Modelling trace metal partitioning in forest floors of northern soils near metal smelters

Trace metal (TM) mobility and toxicity varies with changing soil conditions. Geochemical models can account for the influence of soil characteristics on TM behaviour. We tested the effectiveness of the Stockholm humic model (SHM), and the NICA-Donnan model (NDM) to estimate partitioning coefficients (logKd) in 26 forest floor horizons of podzolic soils enriched in trace metals from deposition by metal smelters. We wanted to know if a consistent approach could be applied to model metal partitioning in forest floors without optimizing each individual soil. When optimized, the SHM reproduced the partitioning of Cd, Cu and Zn but not Pb. It was necessary to revise the affinity constants for the NDM to simultaneously simulate the partitioning of the four metals. Revised affinity constants for the NDM model based on a fixed definition of soil organic carbon, i.e., a fixed ratio of fulvic and humic acids per unit carbon, reproduced metal partitioning more effectively in an independent data set of 16 soils than the use of generic affinity constants available for these models. From the perspective of the applicability of these models to risk assessment, this result suggests geochemical models using affinity constants that have been verified and/or modified against multiple soils from a region can provide good estimates of metal partitioning on a regional scale.     

Trace elements in major solid-pesticides used in the Gaza Strip

This study describes the purity of pesticides used in Gaza in terms of trace elements. A semi-quantitative EMMA-XRF technique and quantitative ICP/OES was used to determine the concentrations of Al, As, Ba, Br, Ca, Cd, Co, Cr, Cu, Fe, K, Mn, Ni, Pb, Rb, Sc, Se, Sr, Ti and Zn in 50 of the most commonly used solid pesticides collected from the five central shops in the Gaza Strip. The results revealed that the pesticides contain considerable amounts of trace elements and do not comply with the expected-theoretical structure of each species. Moreover, they do not reflect the actual constituents mentioned in the trade labels. Interviews with market owners and field surveys confirmed that the pesticides were not pure. In some cases they have been mixed in local markets with minor inorganic species without a scientific basis. They may also have been smuggled into Gaza with differing impurities. The results indicate that pesticides should be considered as a source of certain trace metals (particularly Cu, Mn and Zn) and other elements (Br, Sr and Ti), which may affect their mass balances in soil and groundwater as well as their plant uptake. Different scenarios and calculation models of the transport of trace elements in soil and groundwater of the Gaza Strip should include pesticides as an additional source.     

Reliable QSAR for estimating Koc for persistent organic pollutants: correlation with molecular connectivity indices

Several recent studies have shown that n-octanol/water partition coefficients may not be a good predictor for estimating soil sorption coefficients of persistent organic pollutants (POPs), defined here as chemicals with log Kow greater than 5. Thus, an alternative QSAR model was developed that seems to provide reliable estimates for the soil sorption coefficients of persistent organic pollutants. This model is based on a set of calculated molecular connectivity indices and evaluated soil sorption data for 18 POPs. The chemical's size and shape, quantified by 1chi, 3chiC and 4chiC(v) indices, have a dominant effect on the soil sorption process of POPs. The developed QSAR model was rationalized in terms of potential hydrophobic interactions between persistent organic pollutants and soil organic matrix. Its high predictive power has been verified by an extensive internal and external validation procedure.     

Application of a new generation of complexing agents in removal of heavy metal ions from different wastes

Complexing agents are extensively applied in many fields of industry. They are used to provide effective controlling trace metal ions in cleaning industries, textile, pulp and paper production, water treatment, agriculture, food industries, etc. Recently, the low biodegradability of these ligands and their accumulation in the environment has become a cause for concern. Therefore, replacement of ethylenediaminetetraacetic acid and diethylenetriaminepentaacetic acid by more environmentally friendly chelating agents is highly desirable. So far, these acids and their salts have been applied as components of household chemistry, cosmetics, modern microelement fertilizers and agrochemicals. This paper reviews the sorption of heavy metal ions such as Cu(II), Zn(II), Cd(II) and Pb(II) in the presence of the above-mentioned complexing agents on commercially available anion exchangers of different matrix. The obtained sorption results were fitted using the Langmuir and Freundlich sorption isotherm models. The kinetic data were also analysed using the Lagergren, Ho and McKay sorption kinetic equations. The studies were carried out considering the effects of such important parameters as phase contact time, initial concentration, pH and temperature.     

Comparison of nonideal sorption formulations in modeling the transport of phthalate esters through packed soil columns

Sorption of dimethyl phthalate (DMP), diethyl phthalate (DEP) and dipropyl phthalate (DPP) to two soil materials that vary in organic matter content was investigated using miscible displacement experiments under saturated flow conditions. Generated breakthrough curves (BTCs) were inversely simulated using linear, equilibrium sorption (LE), nonlinear, equilibrium sorption (NL), linear, first-order nonequilibrium sorption (LFO), linear, radial diffusion (LRD), and nonlinear, first-order nonequilibrium sorption (NFO) models. The Akaike information criterion was utilized to determine the preferred model. The LE model could not adequately describe phthalate ester (PE) BTCs in higher organic matter soil or for more hydrophobic PEs. The LFO and LRD models adequately described the BTCs but a slight improvement in curve-fitting was gained in some cases when the NFO model was used. However, none of the models could properly describe the desorptive tail of DPP for the high organic matter soil. Transport of DPP through this soil was adequately predicted when degradation or sorption hysteresis was considered. Using the optimized parameter values along with values reported by others it was shown that the organic carbon distribution coefficient (K(oc)) of PEs correlates well with the octanol/water partition coefficient (K(ow)). Also, a strong relationship was found between the first-order sorption rate coefficient normalized to injection pulse size and compound residence time. A similar trend of timescale dependence was found for the rate parameter in the radial diffusion model. Results also revealed that the fraction of instantaneous sorption sites is dependent on K(ow) and appears to decrease with the increase in the sorption rate parameter.     

A model for estimating the potential biomagnification of chemicals in a generic food web: Preliminary development

GOAL, SCOPE AND BACKGROUND: Bioaccumulation and biomagnification of organic pollutants have been increasingly assessed and modeled during the last years. Due to the complexity of these processes and the large variability of food webs, setting generic assessments for these parameters is really difficult. Equilibrium models, based on a compound's lipophylicity, are the main tool in regulatory proposals, such as for identifying Persistent, Bioaccumulative and Toxic Substances (PBTs), although a refinement has been claimed by the scientific community. Toxicokinetic studies offer an alternative for these estimations, where biomagnification is modeled as a succession of bioaccumulation processes, each one regulated by toxicokinetic parameters. METHODS: A review of kinetic models covering species belonging to different trophic levels and with different ecological behavior has been conducted. The results were employed for setting a conceptual model for estimating the biomagnification potential in a generic food web, which was mathematically implemented through system dynamic models developed under data sheet software. Crystal Ball was then employed for allowing Monte Carlo based probabilistic calculations. Bioaccumulation laboratory assays have been performed to estimate toxicokinetic parameters in mussels (Mytilus edulis) with two PAHs (chrysene and benzo[a]pyrene). The contamination was delivered via food. The exposure period lasted more than one month followed then by a depuration phase. The contaminant content was determined on an individual basis on five replicates. RESULTS AND DISCUSSION:. The reviewed information suggested the development of a tiered conceptual biomagnification model, starting with a simplified food chain which can be refined to more realistic and complex models in successive levels. CONCLUSIONS: The mathematical implementation of the conceptual model offers tools for estimating the potential for bioaccumulation and biomagnification of chemicals under very different conditions. The versatility of the model can be used for both comparative estimations and for validating the model. RECOMMENDATIONS AND PERSPECTIVES: Since bioaccumulation and biomagnification processes are crucial elements for a proper risk assessment of chemicals, their estimation by mathematical models has been widely tested. However, inregulatory assessments, too simplistic models are still being used quite often. The biomagnification model presented in this study should be amore accurate alternative to these models. In comparison to other previously published biomagnification models, the present one covers the time variation of bioaccumulation using just a few toxicokinetic parameters.     

Transport and transformation of sulfadiazine in soil columns packed with a silty loam and a loamy sand

Concerning the transport of the veterinary antibiotic sulfadiazine (SDZ) little is known about its possible degradation during transport. Also its sorption behaviour is not yet completely understood. We investigated the transport of SDZ in soil columns with a special emphasis on the detection of transformation products in the outflow of the soil columns and on modelling of the concentration distribution in the soil columns afterwards. We used disturbed soil columns near saturation, packed with a loamy sand and a silty loam. SDZ was applied as a 0.57 mg L(-1) solution at a constant flow rate of 0.25 cm h(-1) for 68 h. Breakthrough curves (BTC) of SDZ and its transformation products 4-(2-iminopyrimidin-1(2H)-yl)aniline and 4-hydroxy-SDZ were measured for both soils. For the silty loam we additionally measured a BTC for an unknown transformation product which we only detected in the outflow samples of this soil. After the leaching experiments the (14)C-concentration was quantified in different layers of the soil columns. The transformation rates were low with mean SDZ mass fractions in the outflow samples of 95% for the loamy sand compared to 97% for the silty loam. The formation of 4-(2-iminopyrimidin-1(2H)-yl)aniline appears to be light dependent and did probably not occur in the soils, but afterwards. In the soil columns most of the (14)C was found near the soil surface. The BTCs in both soils were described well by a model with one reversible (kinetic) and one irreversible sorption site. Sorption kinetics played a more prominent role than sorption capacity. The prediction of the (14)C -concentration profiles was improved by applying two empirical models other than first order to predict irreversible sorption, but also these models were not able to describe the (14)C concentration profiles correctly. Irreversible sorption of sulfadiazine still is not well understood.     

Enrichment of trace elements in the clay size fraction of mining soils

Reactive waste dumps with sulfide minerals promote acid mine drainage (AMD), which results in water and soil contamination by metals and metalloids. In these systems, contamination is regulated by many factors, such as mineralogical composition of soil and the presence of sorption sites on specific mineral phases. So, the present study dedicates itself to understanding the distribution of trace elements in different size fractions (<2-mm and <2-μm fractions) of mining soils and to evaluate the relationship between chemical and mineralogical composition. Cerdeirinha and Penedono, located in Portugal, were the waste dumps under study. The results revealed that the two waste dumps have high degree of contamination by metals and arsenic and that these elements are concentrated in the clay size fraction. Hence, the higher degree of contamination by toxic elements, especially arsenic in Penedono as well as the role of clay minerals, jarosite, and goethite in retaining trace elements has management implications. Such information must be carefully thought in the rehabilitation projects to be planned for both waste dumps.     

Parameterization and evaluation of sulfate adsorption in a dynamic soil chemistry model

Sulfate adsorption was implemented in the dynamic, multi-layer soil chemistry model SAFE. The process is modeled by an isotherm in which sulfate adsorption is considered to be fully reversible and dependent on sulfate concentration as well as pH in soil solution. The isotherm was parameterized by a site-specific series of simple batch experiments at different pH (3.8-5.0) and sulfate concentration (10-260 micromol 1(-1)) levels. Application of the model to the Lake Gardsj6n roof covered site shows that including sulfate adsorption improves the dynamic behavior of the model and sulfate adsorption and desorption delay acidification and recovery of the soil. The modeled adsorbed pool of sulfate at the site reached a maximum level of 700 mmol/m(2) in the late 1980s, well in line with experimental data.     

Development of a black carbon-inclusive multi-media model: application for PAHs in Stockholm

A multi-media model was developed for predicting the fate of organic chemicals in the Greater Stockholm Area, Sweden, and applied to selected polycyclic aromatic hydrocarbons (PAHs). Although urban models have been previously developed, this model is novel in that it includes sorption to pyrogenically-derived particles, commonly termed "black carbon" (BC), within the model structure. To examine the influence of BC sorption on environmental fate of PAHs, two versions of the model were generated and run: one in which sorption to BC was included and one in which BC sorption was excluded. The inclusion of BC sorption did not cause any significant variations to air levels, but it did cause an average 20-30% increase in sediment concentrations related to increased sediment solids partitioning. The model also predicted reduced advective losses out of the model domain, as well as chemical potential to diffuse from sediments, whilst total chemical inventory increased. In all cases, the lighter PAHs were more affected by BC inclusion than their heavier counterparts. We advocate the addition of sorption to BC in future multi-media fate and exposure models, which as well as influencing fate will also alter (lower) chemical availability and, thus, wildlife exposure to hydrophobic chemicals. A quantification of the latter was derived with the help of the soot-inclusive model version, which estimated a lowering of dissolved water concentrations between five and >200 times for the different PAHs of this study.     

Dynamic studies on the mobility of trace elements in soil and sediment samples influenced by dumping of residues of the flood in the Mulde River region in 2002

In the analysis of soil samples, batch sequential extraction procedures are traditionally used for the fractionation of trace elements to access their mobility and potential risk for the contamination of groundwater. In the present work a continuous-flow technique has been used that enables not only the fast and efficient leaching of trace elements but as well as time-resolved studies on the mobilization of arsenic and selected heavy metals in different forms to be made. Rotating coiled columns (RCC) earlier used mainly in countercurrent chromatography have been successfully applied to the dynamic leaching of heavy metals from soils contaminated by flooding sludge's. The sample was retained in a PTFE rotating column as the stationary phase whereas aqueous solutions were continuously pumped through. The contents of elements were determined by on-line coupling of RCC and inductively coupled plasma atomic emission spectrometry (ICP-AES). This enables real-time data on the leaching process to be obtained. Dynamic and traditional batch procedures were compared. It has been shown that the aqueous elution under centrifugal forced conditions is much more effective for the mobilization of heavy metals. Hence, the dynamic leaching is characterized by a substantially more intensive interaction between solid and water and is besides substantially more time-saving than the conventional batch procedure. The RCC procedure was also employed for preliminary leaching studies with a simulated "acid rain". In comparison with the water leaching, the mobilization of heavy metals and arsenic from soil samples with employment of simulated acid rain as eluent was less effective.     

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.     

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.