Linear Free Energy Relationship for Predicting the Rate Constants of Munition Compound Reduction by the Fe(II)-Hematite and Fe(II)-Goethite Redox Couples.

Abiotic reduction by iron minerals is arguably the most important fate process for munition compounds (MCs) in subsurface environments. No model currently exists that can predict the abiotic reduction rates of structurally diverse MCs by iron (oxyhydr)oxides. We performed batch experiments to measure the rate constants for the reduction of three classes of MCs (poly-nitroaromatics, nitramines, and azoles) by hematite or goethite in the presence of aqueous Fe2+. The surface area-normalized reduction rate constant (kSA) depended on the aqueous-phase one-electron reduction potential (EH1) of the MC and the thermodynamic state (i.e., pe and pH) of the iron oxide-Feaq2+ system. A linear free energy relationship (LFER), similar to that reported previously for nitrobenzene, successfully captures all MC reduction rate constants that span 6 orders of magnitude: log(kSA)=(1.12±0.04)[0.53EH159mV-(pH+pe)]+(5.52±0.23). The finding that the rate constants of all the different classes of MCs can be described by a single LFER suggests that these structurally diverse nitro compounds are reduced by iron oxide-Feaq2+ couples through a common mechanism up to the rate-limiting step. Multiple mechanistic implications of the results are discussed. This study expands the applicability of the LFER model for predicting the reduction rates of legacy and emerging MCs and potentially other nitro compounds.

Thermodynamic Two-Site Surface Reaction Model for Predicting Munition Constituent Reduction Kinetics with Iron (Oxyhydr)oxides.

Iron (oxyhydr)oxides comprise a significant portion of the redox-active fraction of soils and are key reductants for remediation of sites contaminated with munition constituents (MCs). Previous studies of MC reduction kinetics with iron oxides have focused on the concentration of sorbed Fe(II) as a key parameter. To build a reaction kinetic model, it is necessary to predict the concentration of sorbed Fe(II) as a function of system conditions and the redox state. A thermodynamic framework is formulated that includes a generalized double-layer model that utilizes surface acidity and surface complexation reactions to predict sorbed Fe(II) concentrations that are used for fitting MC reduction kinetics. Monodentate- and bidentate Fe(II)-binding sites are used with individual oxide sorption characteristics determined through data fitting. Results with four oxides (goethite, hematite, lepidocrocite, and ferrihydrite) and four nitro compounds (NB, CN-NB, Cl-NB, and NTO) from six separate studies have shown good agreement when comparing observed and predicted surface area-normalized rate constants. While both site types are required to reproduce the experimental redox titration, only the monodentate site concentration controls the MC reaction kinetics. This model represents a significant step toward predicting the timescales of MC degradation in the subsurface.

Electron Transfer Energy and Hydrogen Atom Transfer Energy-Based Linear Free Energy Relationships for Predicting the Rate Constants of Munition Constituent Reduction by Hydroquinones.

No single linear free energy relationship (LFER) exists that can predict reduction rate constants of all munition constituents (MCs). To address this knowledge gap, we measured the reduction rates of MCs and their surrogates including nitroaromatics [NACs; 2,4,6-trinitrotoluene (TNT), 2,4-dinitroanisole (DNAN), 2-amino-4,6-dinitrotoluene (2-A-DNT), 4-amino-2,6-dinitrotoluene (4-A-DNT), and 2,4-dinitrotoluene (DNT)], nitramines [hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) and nitroguanidine (NQ)], and azoles [3-nitro-1,2,4-triazol-5-one (NTO) and 3,4-dinitropyrazole (DNP)] by three dithionite-reduced quinones (lawsone, AQDS, and AQS). All MCs/NACs were reduced by the hydroquinones except NQ. Hydroquinone and MC speciations were varied by controlling pH, permitting the application of a speciation model to determine second-order rate constants (k) from observed pseudo-first-order rate constants. The intrinsic reactivity of MCs (oxidants) decreased upon deprotonation, while the opposite was true for hydroquinones (reductants). The rate constants spanned ∼6 orders of magnitude in the order NTO ≈ TNT > DNP > DNT ≈ DNAN ≈ 2-A-DNT > DNP- > 4-A-DNT > NTO- > RDX. LFERs developed using density functional theory-calculated electron transfer and hydrogen atom transfer energies and reported one-electron reduction potentials successfully predicted k, suggesting that these structurally diverse MCs/NACs are all reduced by hydroquinones through the same mechanism and rate-limiting step. These results increase the applicability of LFER models for predicting the fate and half-lives of MCs and related nitro compounds in reducing environments.

Modeling the Reduction Kinetics of Munition Compounds by Humic Acids.

Dissolved organic matter (DOM) comprises a sizeable portion of the redox-active constituents in the environment and is an important reductant for the abiotic transformation of nitroaromatic compounds and munition constituents (NACs/MCs). Building a predictive kinetic model for these reactions would require the energies associated with both the reduction of the NACs/MCs and the oxidation of the DOM. The heterogeneous and unknown structure of DOM, however, has prohibited reliable determination of its oxidation energies. To overcome this limitation, humic acids (HAs) were used as model DOM, and their redox moieties were modeled as a collection of quinones of different redox potentials. The reduction and oxidation energies of the NACs/MCs and hydroquinones, respectively, via hydrogen atom transfer (HAT) reactions were then calculated quantum chemically. HAT energies have been used successfully in a linear free energy relationship (LFER) to predict second-order rate constants for NAC reduction by hydroquinones. Furthermore, a linear relationship between the HAT energies and the reduction potentials of quinones was established, which allows estimation of hydroquinone reactivity (i.e., rate constants) from HA redox titration data. A training set of three HAs and two NACs/MCs was used to generate a mean HA redox profile that successfully predicted reduction kinetics in multiple HA/MC systems.

Reductive Transformation of 3-Nitro-1,2,4-triazol-5-one (NTO) by Leonardite Humic Acid and Anthraquinone-2,6-disulfonate (AQDS).

3-Nitro-1,2,4-triazol-5-one (NTO) is a major and the most water-soluble constituent in the insensitive munition formulations IMX-101 and IMX-104. While NTO is known to undergo redox reactions in soils, its reaction with soil humic acid has not been evaluated. We studied NTO reduction by anthraquinone-2,6-disulfonate (AQDS) and Leonardite humic acid (LHA) reduced with dithionite. Both LHA and AQDS reduced NTO to 3-amino-1,2,4-triazol-5-one (ATO), stoichiometrically at alkaline pH and partially (50-60%) at pH ≤ 6.5. Due to NTO and hydroquinone speciation, the pseudo-first-order rate constants (kObs) varied by 3 orders of magnitude from pH 1.5 to 12.5 but remained constant from pH 4 to 10. This distinct pH dependency of kObs suggests that NTO reactivity decreases upon deprotonation and offsets the increasing AQDS reactivity with pH. The reduction of NTO by LHA deviated continuously from first-order behavior for >600 h. The extent of reduction increased with pH and LHA electron content, likely due to greater reactivity of and/or accessibility to hydroquinone groups. Only a fraction of the electrons stored in LHA was utilized for NTO reduction. Electron balance analysis and LHA redox potential profile suggest that the physical conformation of LHA kinetically limited NTO access to hydroquinone groups. This study demonstrates the importance of carbonaceous materials in controlling the environmental fate of NTO.

A Unified Linear Free Energy Relationship for Abiotic Reduction Rate of Nitroaromatics and Hydroquinones Using Quantum Chemically Estimated Energies.

Determining the fate of nitroaromatic compounds (NACs) in the environment requires the use of predictive models for compounds and conditions for which experimental data are insufficient. Previous studies have developed linear free energy relationships (LFERs) that relate the thermodynamic energy of NAC reduction to its corresponding rate constant. We present a comprehensive LFER that incorporates both the reduction and oxidation half-reactions through quantum chemically calculated energies. Environ Toxicol Chem 2020;39:2389-2395. © 2020 SETAC.

Reduction of 3-Nitro-1,2,4-Triazol-5-One (NTO) by the Hematite-Aqueous Fe(II) Redox Couple.

3-Nitro-1,2,4-triazol-5-one (NTO) is an insensitive munition compound (MC) that has replaced legacy MC. NTO can be highly mobile in soil and groundwater due to its high solubility and anionic nature, yet little is known about the processes that control its environmental fate. We studied NTO reduction by the hematite-Fe2+ redox couple to assess the importance of this process for the attenuation and remediation of NTO. Fe2+(aq) was either added (type I) or formed through hematite reduction by dithionite (type II). In the presence of both hematite and Fe2+(aq), NTO was quantitatively reduced to 3-amino-1,2,4-triazol-5-one following first-order kinetics. The surface area-normalized rate constant (kSA) showed a strong pH dependency between 5.5 and 7.0 and followed a linear free energy relationship (LFER) proposed in a previous study for nitrobenzene reduction by iron oxide-Fe2+ couples, i.e., log kSA = -(pe + pH) + constant. Sulfite, a major dithionite oxidation product, lowered kSA in type II system by ∼10-fold via at least two mechanisms: by complexing Fe2+ and thereby raising pe, and by making hematite more negatively charged and hence impeding NTO adsorption. This study demonstrates the importance of iron oxide-Fe2+ in controlling NTO transformation, presents an LFER for predicting NTO reduction rate, and illustrates how solutes can shift the LFER by interacting with either iron species.

Hydrogen Atom Transfer Reaction Free Energy as a Predictor of Abiotic Nitroaromatic Reduction Rate Constants: A Comprehensive Analysis.

A linear free energy model is presented that predicts the second-order rate constant for the abiotic reduction of nitroaromatic compounds (NACs). Previously presented models use the one-electron reduction potential EH1(ArNO2) of the NAC reaction ArNO2+e-→ArNO2•- . If EH1(ArNO2) is not available, it has been proposed that EH1(ArNO2) be computed directly or estimated from the gas-phase electron affinity (EA). The model proposed uses the Gibbs free energy of the hydrogen atom transfer (HAT) reaction ArNO2+H•→ArNOOH• as the parameter in the linear free energy model. Both models employ quantum chemical computations for the required thermodynamic energies. The available and proposed models are compared using experimentally determined second-order rate constants from 5 investigations from the literature in which a variety of NACs were exposed to a variety of reductants. A comprehensive analysis utilizing all the NACs and reductants demonstrate that the HAT energy model and the experimental one-electron reduction potential model have similar root mean square errors and residual error probability distributions. In contrast, the model using the computed EA has a more variable residual error distribution with a significant number of outliers. The results suggest that a linear free energy model utilizing computed HAT reaction free energy produces a more reliable prediction of the NAC abiotic reduction second-order rate constant than previously available methods. The advantages of the proposed HAT energy model and its mechanistic implications are discussed as well. Environ Toxicol Chem 2020;39:1678-1684. © 2020 SETAC.

Predicting Cr(vi) adsorption on soils: the role of the competition of soil organic matter.

Cr(vi) has posed a serious risk for the environment and human beings because of its pollution and toxicity. It is essential to understand the equilibrium behavior of Cr(vi) in soils. In this study, the adsorption of Cr(vi) on fourteen soils was studied with batch experiments and quantitative modeling. The batch experiments included the adsorption edge and adsorption isotherm experiments, investigating the adsorption of Cr(vi) with varying soil properties, solution pH, and initial Cr(vi) concentrations. The experimental data were then modeled using the surface complexation models in Visual MINTEQ of CD-MUSIC by considering the adsorption of Cr(vi) and ions onto Fe (hydr)oxides and Al (hydr)oxides, and the Stockholm Humic Model and the fixed charge site model by accounting for the adsorption of the cations to soil organic matter and clay, respectively. Particularly, the modeling method of this study introduced an important parameter RO- to account for the amount of soil organic matter irreversibly adsorbed on soil minerals. Overall, the model predicted reasonably well for the equilibrium partition of Cr(vi) under various conditions with a root-mean-square-error of 0.35 for the adsorption edge data and 0.19 for the adsorption isotherm data. According to the model calculations, ferrihydrite dominated the binding of Cr(vi) at pH of 3.0-7.0. The content of ferrihydrite and reactive soil organic matter was found to be the main factor influencing RO-. The modeling results help to understand and predict Cr(vi) adsorption on different soils and are beneficial to environmental risk assessment and pollution remediation.

Experimental Validation of Hydrogen Atom Transfer Gibbs Free Energy as a Predictor of Nitroaromatic Reduction Rate Constants.

Nitroaromatic compounds (NACs) are a class of prevalent contaminants. Abiotic reduction is an important fate process that initiates NAC degradation in the environment. Many linear free energy relationship (LFER) models have been developed to predict NAC reduction rates. Almost all LFERs to date utilize experimental aqueous-phase one-electron reduction potential ( EH1) of NAC as a predictor, and thus, their utility is limited by the availability of EH1 data. A promising new approach that utilizes computed hydrogen atom transfer (HAT) Gibbs free energy instead of EH1 as a predictor was recently proposed. In this study, we evaluated the feasibility of HAT energy for predicting NAC reduction rate constants. Using dithionite-reduced quinones, we measured the second-order rate constants for the reduction of seven NACs by three hydroquinones of different protonation states. We computed the gas-phase energies for HAT and electron affinity (EA) of NACs and established HAT- and EA-based LFERs for six hydroquinone species. The results suggest that HAT energy is a reliable predictor of NAC reduction rate constants and is superior to EA. This is the first independent, experimental validation of HAT-based LFER, a new approach that enables rate prediction for a broad range of structurally diverse NACs based solely on molecular structures.

Sorption and desorption kinetics of nitroglycerin and 2,4-dinitrotoluene in nitrocellulose and implications for residue-bound energetic materials.

Energetic materials (EMs) bound to propellant residues can contribute to environmental risk and public health concerns. This work investigated how nitrocellulose, a common binding material in propellants, may control the release dynamics of nitroglycerin (NG) and 2,4-dinitrotoluene (2,4-DNT) from propellant residues. Batch adsorption/desorption experiments on nitrocellulose and re-interpretation on results from past leaching studies involving propellant-bound EMs were conducted. Mechanistic modeling of adsorption/desorption kinetics based on intra-particle diffusion (IPD) predicted aqueous intrinsic diffusivities (Diw) to within a factor of 2 of expected values. Furthermore, the IPD model was able to predict effective diffusivities (Deff) during the early leaching of NG from propellant residues to within a factor of 2 over a 3-log unit range. Prediction of leaching Deff's associated with fired residues was less successful probably due to the neglect of compositional and morphological heterogeneity within the residues. Close correlations were found between the early and late Deff's of residue-bound NG and between the fast- and slow-domain rate constants for both EMs, suggesting that the late leaching kinetics of bound-EMs may be empirically assessed from the early kinetics. This work illustrates that, in addition to dissolution, retarded diffusion through nitrocellulose matrix may also limit the overall release and transformation of residue-bound EMs in the field. Implications and limitations of the current study, and the steps forward are also presented.

Validation of Cu toxicity to barley root elongation in soil with a Terrestrial Biotic Ligand Model developed from sand culture.

Constants for a Terrestrial Biotic Ligand Model (TBLM) to predict the Cu toxicity to barley root elongation (RE) were developed from controlled sand culture experiments. These constants were used to predict RE in soil culture. The competition of H+, Ca2+, and Mg2+ to Cu2+ toxicity were studied individually and independently, and linear relationships between EC50 free Cu2+ and H+, Ca2+, and Mg2+ activities were found, meaning that the cations H+, Ca2+, and Mg2+ will alleviate the toxicity of Cu2+ in solutions. Toxicity accompanying increasing concentration of solution ions other than Cu2+ was observed and modeled as an osmotic effect which improved soil culture toxicity prediction. The Root Mean Square Error (RMSE) of %RE and EC50 (50% effective concentration) for soil toxicity prediction using TBLM parameters developed from sand culture are 13.0 and 0.23 respectively, which are as good as that of 14.0 and 0.24 using parameters that developed from soil culture itself. A model including the activity at the root plasma membrane surface was tested and found not to provide improvement over the use of bulk solution activity to predict metal toxicity. TBLM parameters obtained from water solution culture were unable to accurately predict the EC50s in soils whereas the parameters obtained from sand culture were able to predict the toxicity in soils. Including the toxicity of CuOH+ was found to improve the toxicity prediction slightly.

Bioconcentration factors and plant-water partition coefficients of munitions compounds in barley.

Plants growing in the soils at military ranges and surrounding locations are exposed, and potentially able to uptake, munitions compounds (MCs). The extent to which a compound is transferred from the environment into organisms such as plants, referred to as bioconcentration, is conventionally measured through uptake experiments with field/synthetic soils. Multiple components/phases that vary among different soil types and affect the bioavailability of the MC, however, hinder the ability to separate the effects of soil characteristics from the MC chemical properties on the resulting plant bioconcentration. To circumvent the problem, this work presents a protocol to measure steady state bioconcentration factors (BCFs) for MCs in barley (Hordeum vulgare L.) using inert laboratory sand rather than field/synthetic soils. Three MCs: 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (2,4-DNT), and 2,4-dinitroanisole (2,4-DNAN), and two munition-like compounds (MLCs): 4-nitroanisole (4-NAN) and 2-methoxy-5-nitropyridine (2-M-5-NPYNE) were evaluated. Approximately constant plant biomass and exposure concentrations were achieved within a one-month period that produced steady state log BCF values: 0.62 ± 0.02, 0.70 ± 0.03, 1.30 ± 0.06, 0.52 ± 0.03, and 0.40 ± 0.05 L kgplant dwt-1 for TNT, 2,4-DNT, 2,4-DNAN, 4-NAN, and 2-M-5-NPYNE, respectively. Furthermore, results suggest that the upper-bounds of the BCFs can be estimated within an order of magnitude by measuring the partitioning of the compounds between barley biomass and water. This highlights the importance of partition equilibrium as a mechanism for the uptake of MCs and MLCs by barley from interstitial water. The results from this work provide chemically meaningful data for prediction models able to estimate the bioconcentration of these contaminants in plants.

Leaching of propellant compounds from munition residues may be controlled by sorption to nitrocellulose.

Sustainable management of military ranges requires effective assessment of surface mobility and leaching potential of propellant compounds (PCs). Previous studies have focused mostly on PCs' dissolution from fired residues and their sorption to soil components. This work investigated the potential role of nitrocellulose, a major component in propellants, in the binding of PCs to propellant residues. Sorption isotherms of military grade nitrocellulose for dissolved nitroglycerine (NG) or 2,4-dinitrotoluene (2,4-DNT) was measured in batch experiments and were determined to be SNG=102.39(±0.05)CNG0.916(±0.032) and S2,4-DNT=103.08(±0.01)C2,4-DNT0.668(±0.010) (S and C in mg/kgnitrocellulose and mg/Lwat, respectively). Solid-to-water partitioning for NG and 2,4-DNT was 100 times greater in propellant residues than in typical military ranges soils. Since nitrocellulose can sorb NG and 2,4-DNT up to 23 and 5% of its mass, respectively, it can slow down, through retarded diffusion, the leaching of PCs from fired residues over the typical composition ranges of common propellants. The slow leaching of PCs from propellant grains in column studies can be better interpreted by considering their sorptive interaction with nitrocellulose in addition to dissolution kinetics. With nitrocellulose as the carrying matrix, residue-bound PCs may migrate farther and persist longer in subsurface environment.

Experimental determination of solvent-water partition coefficients and Abraham parameters for munition constituents.

There is concern about the environmental fate and effects of munition constituents (MCs). Polyparameter linear free energy relationships (pp-LFERs) that employ Abraham solute parameters can aid in evaluating the risk of MCs to the environment. However, poor predictions using pp-LFERs and ABSOLV estimated Abraham solute parameters are found for some key physico-chemical properties. In this work, the Abraham solute parameters are determined using experimental partition coefficients in various solvent-water systems. The compounds investigated include hexahydro-1,3,5-trinitro-1,3,5-triazacyclohexane (RDX), octahydro-1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), hexahydro-1-nitroso-3,5-dinitro-1,3,5-triazine (MNX), hexahydro-1,3,5-trinitroso-1,3,5-triazine (TNX), hexahydro-1,3-dinitroso-5- nitro-1,3,5-triazine (DNX), 2,4,6-trinitrotoluene (TNT), 1,3,5-trinitrobenzene (TNB), and 4-nitroanisole. The solvents in the solvent-water systems are hexane, dichloromethane, trichloromethane, octanol, and toluene. The only available reported solvent-water partition coefficients are for octanol-water for some of the investigated compounds and they are in good agreement with the experimental measurements from this study. Solvent-water partition coefficients fitted using experimentally derived solute parameters from this study have significantly smaller root mean square errors (RMSE = 0.38) than predictions using ABSOLV estimated solute parameters (RMSE = 3.56) for the investigated compounds. Additionally, the predictions for various physico-chemical properties using the experimentally derived solute parameters agree with available literature reported values with prediction errors within 0.79 log units except for water solubility of RDX and HMX with errors of 1.48 and 2.16 log units respectively. However, predictions using ABSOLV estimated solute parameters have larger prediction errors of up to 7.68 log units. This large discrepancy is probably due to the missing R2NNO2 and R2NNO2 functional groups in the ABSOLV fragment database.

Barley root hair growth and morphology in soil, sand, and water solution media and relationship with nickel toxicity.

Barley, Hordeum vulgare (Doyce), was grown in the 3 media of soil, hydroponic sand solution (sand), and hydroponic water solution (water) culture at the same environmental conditions for 4 d. Barley roots were scanned, and root morphology was analyzed. Plants grown in the 3 media had different root morphology and nickel (Ni) toxicity response. Root elongations and total root lengths followed the sequence soil > sand > water. Plants grown in water culture were more sensitive to Ni toxicity and had greater root hair length than those from soil and sand cultures, which increased root surface area. The unit root surface area as root surface area per centimeter of length of root followed the sequence water > sand > soil and was found to be related with root elongation. Including the unit root surface area, the difference in root elongation and 50% effective concentration were diminished, and percentage of root elongations can be improved with a root mean square error approximately 10% for plants grown in different media. Because the unit root surface area of plants in sand culture is closer to that in soil culture, the sand culture method, not water culture, is recommended for toxicity parameter estimation. Environ Toxicol Chem 2016;35:2125-2133. © 2016 SETAC.

A metabolomic study on the responses of daphnia magna exposed to silver nitrate and coated silver nanoparticles.

We examined the short-term toxicity of AgNPs and AgNO3 to Daphnia magna at sublethal levels using (1)H NMR-based metabolomics. Two sizes of polyvinylpyrrolidone-coated AgNPs (10 and 40nm) were synthesized and characterized and their Ag(+) release was studied using centrifugal ultrafiltration and inductively coupled plasma mass spectrometry. Multivariate statistical analysis of the (1)H NMR spectra showed significant changes in the D. magna metabolic profiles following 48h exposure to both AgNP particle sizes and Ag(+) exposure. Most of the metabolic biomarkers for AgNP exposure, including 3-hydroxybutyrate, arginine, lysine and phosphocholine, were identical to those of the Ag(+)-exposed groups, suggesting that the dominant effects of both AgNPs were due to released Ag(+). The observed metabolic changes implied that the released Ag(+) induced disturbance in energy metabolism and oxidative stress, a proposed mechanism of AgNP toxicity. Elevated levels of lactate in all AgNP-treated but not in Ag(+)-treated groups provided evidence for Ag-NP enhanced anaerobic metabolism. These findings show that (1)H NMR-based metabolomics provides a sensitive measure of D. magna response to AgNPs and that further targeted assays are needed to elucidate mechanisms of action of nanoparticle-induced toxicity.

Development and validation of a terrestrial biotic ligand model for Ni toxicity to barley root elongation for non-calcareous soils.

A Terrestrial Biotic Ligand Model (TBLM) for Ni toxicity to barley root elongation (RE) developed from experiments conducted in sand culture was used to predict toxicity in non-calcareous soils. Ca(2+) and Mg(2+) concentrations and pH in sand solution were varied individually and TBLM parameters were computed. EC50 increased as Mg(2+) increased, whereas the effect of Ca(2+) was insignificant. TBLM parameters developed from sand culture were validated by toxicity tests in eight Ni-amended, non-calcareous soils. Additional to Ni(2+) toxicity, toxicity from all solution ions was modelled independently as an osmotic effect and needed to be included for soil culture results. The EC50s and EC10s in soil culture were predicted within twofold of measured results. These are close to the results obtained using parameters estimated from the soil culture data itself.

Soil acidification increases metal extractability and bioavailability in old orchard soils of Northeast Jiaodong Peninsula in China.

The bioavailability of Cu, Zn, Pb and Cd from field-aged orchard soils in a certified fruit plantation area of the Northeast Jiaodong Peninsula in China was assessed using bioassays with earthworms (Eisenia fetida) and chemical assays. Soil acidity increased with increasing fruit cultivation periods with a lowest pH of 4.34. Metals were enriched in topsoils after decades of horticultural cultivation, with highest concentrations of Cu (132 kg(-1)) and Zn (168 mg kg(-1)) in old apple orchards and Pb (73 mg kg(-1)) and Cd (0.57 mg kg(-1)) in vineyard soil. Earthworm tissue concentrations of Cu and Pb significantly correlated with 0.01 M CaCl2-extractable soil concentrations (R(2) = 0.70, p < 0.001 for Cu; R(2) = 0.58, p < 0.01 for Pb). Because of the increased bioavailability, regular monitoring of soil conditions in old orchards and vineyards is recommended, and soil metal guidelines need reevaluation to afford appropriate environmental protection under acidifying conditions.

A general model for kinetics of heavy metal adsorption and desorption on soils.

In this study, we propose a general kinetics model for heavy metal adsorption and desorption reactions in soils when soil organic matter (SOM) is the dominant adsorbent. The kinetics model, integrated with the equilibrium speciation model WHAM VI, specifically considers metal reactions with SOM and dissolved organic matter (DOM) and accounts for the variations of solution chemistry. Metal reactions with SOM are associated with two groups of sites, one from the monodentate sites and another one from the bidentate and tridentate sites. There are three model parameters, desorption rate coefficients of the two groups of SOM sites for each metal and reactive organic carbon (ROC) for each soil. The applicability of the kinetics model was mainly examined with three elements, Cu, Pb, and Zn, which demonstrate different binding ability with organic matter. The kinetic data were collected with a stirred-flow reactor covering a wide range of experimental conditions, including varying SOM, DOM, Ca, and metal concentrations, reaction pHs, and different flow rates. The kinetics model has been successfully applied to describe heavy metal adsorption and desorption on soils under various reaction conditions.