The solid-solution partitioning of heavy metals (Cu, Zn, Cd, Pb) in upland soils of England and Wales

Quantitative analysis of the compliance of EU Sewage Sludge Directive by using the heavy metal concentrations from LUCAS topsoil database

In the European Union (EU), a common understanding of the potential harmful effect of sewage sludge (SS) on the environment is regulated by the Sewage Sludge Directive 86/278/EEC (SSD). Limit values (LVs) for concentrations of heavy metals in soil are listed in Impact Assessment of this directive, and they were transposed by EU member states using different criteria. Member states adopted either single limit values or based on soil factors such as pH and texture to define the maximum limit values for concentrations of heavy metals in soils. Our work presents the first quantitative analysis of the SSD at the European level by using the Land Use and Coverage Area Frame Survey (LUCAS) 2009 topsoil database. The reference values at the European level were arranged taking into account the upper value (EU_UL) and the lower value (EU_LL) for each heavy metal (arsenic, cadmium, copper, chromium, mercury, nickel, lead, and zinc) as well as taking into account the pH of the soil (cadmium, copper, mercury, nickel, lead, and zinc) as introduced in the SSD Annex IA. Single and integrated contamination rate indices were developed to identify those agricultural soils that exceeded the reference values for each heavy metal. In total, 10%, 36%, and 19% of the LUCAS 2009 topsoil samples exceeded the limit values. Additionally, 12% and 16% of agricultural soils exceeded the concentration of at least one single heavy metal when European LVs were fixed following the soil pH in Strategy II compared to those national ones in Strategy I. Generally, all member states apply similar or stricter limit values than those laid down in the SSD. Our work indicates that choosing LVs quantitatively affects further actions such as monitoring and remediation of contaminated soils. The actual soil parameters, such as heavy metal concentrations and soil pH values from the LUCAS 2009 topsoil database, could be used by SSD-involved policy stakeholders not only to lay down the LVs for concentrations of heavy metal in soils but also for monitoring the SSD compliance grade by using the LUCAS surveys over time (past and upcoming LUCAS datasets).

Aging rice straw reduces the bioavailability of mercury and methylmercury in paddy soil

Straw amendment is a prevalent agricultural practice worldwide, which can reduce air pollution and improve soil fertility. However, the impact of aging straw amendment on the bioavailability of mercury (Hg) and methylmercury (MeHg) in paddy soil remains unclear. To investigate this, incubation experiments were conducted using the diffusive gradient in thin-film technique. Results showed that amendments of dry-wet aging (DRS), photochemical aging (LRS), and freeze-thaw aging rice straw (FRS) reduced the bioavailable MeHg in paddy soil by 2.2-27.6%, 13.5-69.8%, and 23.5-86.1%, respectively, compared to fresh rice straw (RS) amendment. This result could be due to changes in soil properties such as soil pH and overlying water Fe and Mn as well as microbial abundance (including Clostridiaceae, Firmicutes, and Actinobacteriota). Simultaneously, The LRS and FRS amendments reduced bioavailable Hg in paddy soil by 20.0-40.8% and 17.1-48.6%, respectively, while DRS increased the bioavailable Hg by 15.8-120.0%. This could be attributed to changes in soil oxidation-reduction potential and overlying water SO42- content. Additionally, the results of sand culture experiments showed that the concentrations of Hg uptake by rice seedlings were 97.1-118.2%, 28.1-35.6%, and 198.0-217.1% higher in dissolved organic matter (DOM) derived from DRS, LRS, and FRS than RS, indicating that aging straw leached DOM may promote the Hg bioavailable when straw amendment. This result could be due to lower molecular weight and higher CO functional group content. These results provide new insight into how aging straw amendment affects the bioavailability of Hg and MeHg in paddy soil under different climates.

Mapping risks associated with soil copper contamination using availability and bio-availability proxies at the European scale

Soil contamination by trace elements like copper (Cu) can affect soil functioning. Environmental policies with guidelines and soil survey measurements still refer to the total content of Cu in soils. However, Cu content in soil solution or free Cu content have been shown to be better proxies of risks of Cu mobility or (bio-)availability for soil organisms. Several empirical equations have been defined at the local scale to predict the amount of Cu in soil solution based on both total soil Cu content and main soil parameters involved in the soil/solution partitioning. Nevertheless, despite the relevance for risk assessment, these equations are not applied at a large spatial scale due to difficulties to perform changes from local to regional. To progress in this challenge, we collected several empirical equations from literature and selected those allowing estimation of the amount of Cu in solution, used as a proxy of available Cu, from the knowledge of both total soil Cu content and soil parameters. We did the same for the estimation of free Cu in solution, used as a proxy of bio-available Cu. These equations were used to provide European maps of (bio-)available Cu based on the one of total soil Cu over Europe. Results allowed comparing the maps of available and bio-available Cu at the European scale. This was done with respective median values of each form of Cu to identify specific areas of risks linked to these two proxies. Higher discrepancies were highlighted between the map of bio-available Cu and the map of soil total Cu compared to the Cu available map. Such results can be used to assess environmental-related issues for land use planning.

Unified Modeling Approach for Quantifying the Proton and Metal Binding Ability of Soil Dissolved Organic Matter

Soil dissolved organic matter (DOM) is composed of a mass of complex organic compounds in soil solutions and significantly affects a range of (bio)geochemical processes in soil environment. However, how the chemical complexity (i.e., heterogeneity and chemodiversity) of soil DOM molecules affects their proton and metal binding ability remains unclear, which limits our ability for predicting the environmental behavior of DOM and metals. In this study, we developed a unified modeling approach for quantifying the proton and metal binding ability of soil DOM based on Cu titration experiments, Fourier transform ion cyclotron resonance mass spectrometry data, and molecular modeling method. Although soil DOM samples from different regions have enormously heterogeneous and diverse properties, we found that the molecules of soil DOM can be divided into three representative groups according to their Cu binding capacity. Based on the molecular models for individual molecular groups and the relative contributions of each group in each soil DOM, we were able to further develop molecular models for all soil DOM to predict their molecular properties and proton and metal binding ability. Our results will help to develop mechanistic models for predicting the reactivity of soil DOM from various sources.

Environmental quality standards for agricultural land in China: What should be improved on derivation methodology?

Soil pollution has caused increasingly widespread attention in China. The environmental risk threshold of pollutants is a yardstick to measure soil environmental quality. For decades, plenty of research on soil environmental quality standards (SEQSs) has been carried out, providing scientific basis for the investigation and supervision of soil environmental quality. This paper summaries the development of SEQSs in China, the corresponding influencing factors and methodology of SEQSs derivation. In the current version of SEQSs (GB15618-2018), the thresholds of soil pollutants are derived by the methods of environmental risk assessment, which are more methodologically scientific than geochemical method and ecological effect method used in the previous version (GB15618-1995). Abundant toxicology data on related species is required for risk assessment of soil pollution using extrapolation; however, basic toxicological data is insufficient and few valid data is available at present. Besides, the inadequate consideration on influencing factors for the derivation of soil pollutant threshold would affect the scientificity and rationality of SEQSs, such as biotic factors (species type, test endpoint etc.) and abiotic factors (aging effect, leaching effect, synergistic or antagonistic effects of elements etc.). These problems should be paid close attention in future research on soil environmental quality standards. The contents summarized in this review may provide reference for decision-making on supervision of soil environmental quality and point out important directions for future studies.

Heavy metals speciation in surface sediments of the Cross River Estuary, Gulf of Guinea, South East Nigeria

The speciation of heavy metals (Ni, Cr, Cu, Zn, Pb, and Cd) was studied in surface sediments of the Cross River Estuary (CRE), Gulf of Guinea, South East Nigeria. Pb (~56 %), Cd (~71 %), Zn (~67 %), and Cr (~76 %) were mainly available in non-residual phases, suggesting potential bioavailability. High contents of Ni and Cu in residual phase indicated immobilization of these metals in aluminosilicate minerals. Cd was the most polluted heavy metal with the highest bioavailability risk. Bayesian Network model results revealed that sedimentary organic carbon (OC) from terrestrial C₃ plants controlled the contents and variability of Pb and Zn, while the input of terrestrial soil OC strongly influenced Cu and Ni. However, Cd and Cr were dominantly influenced by sediment pH, while Ni was mainly influenced by sediment salinity. Strong interdependency between Cd and total nitrogen (TN) suggested that nitrogen might increase Cd bioavailability upon release from sediments.

Effect of land use pattern on the bioavailability of heavy metals: A case study with a multi-surface model

In modern agricultural practice, the land use pattern has been changing due to economic reasons and related policies, which significantly affects the basic physical and chemical properties of soils, thereby influencing the speciation and distribution of heavy metals (HMs) in soils. In this study, we selected three typical types of land use patterns (vegetable field, paddy field and forest field) in Shaoguan City, Guangdong Province, to analyze the content and distribution of HMs, screen the sensitive physicochemical properties, and predict the phytoavailability of HMs under different land use patterns with the multi-surface model (MSM). The forest field had relatively lower levels of labile and free HM ions than both paddy and vegetable fields, which may be attributed to the lower HM content in forest field. The modeling results revealed that organic matter (OM) is the primary carrier of HMs, accounting for 0.19%-97.92% of labile HMs. The sensitivity of soil physicochemical properties to free HM ions followed the order of pH > SOM > goethite > clay. Besides, the conversion of paddy field into vegetable or forest field increased the environmental risk of HMs. Our results may help better decision making in agricultural restructuring to reduce the risk of HM-contaminated soils, as well as give a demonstration for the application of the MSM in predicting the phytoavailability of HMs as a powerful technique.

Estimating remobilization of potentially toxic elements in soil and road dust of an industrialized urban environment

The mobility of potentially toxic elements (PTEs) is of paramount concern in urban settings, particularly those affected by industrial activities. Here, contaminated soils and road dusts of the medium-size, industrialized city of Volos, Central Greece, were subjected to single-step extractions (0.43 M HNO₃ and 0.5 M HCl) and the modified BCR sequential extraction procedure. This approach will allow for a better understanding of the geochemical phase partitioning of PTEs and associated risks in urban environmental matrices. Based on single extraction procedures, Pb and Zn exhibited the highest remobilization potential. Of the non-residual phases, the reducible was the most important for Pb, and the oxidizable for Cu and Zn in both media. On the other hand, mobility of Ni, Cr, and Fe was low, as inferred by their dominance into the residual fraction. Interestingly, we found a significant increase of the residual fraction in the road dust samples compared to soils. Carbonate content and organic matter controlled the extractabilities of PTEs in the soil samples. By contrast, for the road dust, magnetic susceptibility exerted the main control on the geochemical partitioning of PTEs. We suggest that anthropogenic particles emitted by heavy industries reside in the residual fraction of the SEP, raising concerns about the assessment of this fraction in terms of origin of PTEs and potential environmental risks. Conclusively, the application of sequential extraction procedures should be complemented with source identification of PTEs with the aim to better estimate the remobilization of PHEs in soil and road dust influenced by industrial emissions.

Competitive kinetics of Ni(II)/Co(II) and Cr(VI)/P(V) adsorption and desorption on goethite: A unified thermodynamically based model

Predicting the dynamic behavior of coexisting ions on mineral interface is essential to understanding their lability in soil matrix, but a mechanical kinetic model for predicting competitive adsorption is lacking. In this study, the thermodynamic and kinetic adsorption behaviors of Ni(II), Co(II), Cr(VI), and P(V) on goethite under various condition were investigated by batch and stirred-flow experiments, respectively. The equilibrium model CD-MUSIC was developed to describe their equilibrium behavior, followed by the development of a multi-rate kinetic model constrained by the equilibrium model to describe their kinetic behavior. Ni(II) and Co(II) exhibit similar adsorption affinities, while the adsorption of P(V) was stronger and faster than that of Cr(VI). The two surface species of Cr(VI) and P(V) differed in dynamic features, a finding confirmed by in-situ ATR-FTIR spectroscopy. The kinetic model was successfully used to predict the binary competitive adsorption of Ni(II)-Co(II) and Cr(VI)-P(V), and especially the overshooting of Cr(VI) induced by P(V). Our results showed that an integrated thermodynamic-kinetic model obtained from a single-ion experiment can be extended to describe complex multi-ion interactions, indicating the robustness and scalability of the model's parameters. This approach can be used to construct more comprehensive equilibrium and dynamic models of the actual soil environment.

Liming legacy effects associated with the world's largest soil liming and regreening program in Sudbury, Ontario, Canada

Limestone is a common amendment used to counteract soil acidity and metal pollution. Understanding the legacy effects of a one-time soil limestone application and subsequent afforestation is needed to evaluate the long-term success of remediation efforts. In this study, soil and tree chemistry were measured across 15 limed sites that were treated and planted 14 to 37 years ago in Sudbury, Ontario, along with two untreated sites. Soil pH and exchangeable base cation (calcium (Ca) and magnesium (Mg)) concentrations were generally elevated especially in surface organic [FH] horizons up to 37-years post limestone treatment. High site-to-site variation however, obscured clear patterns over time and base cation mass budgets were generally unable to account for the mass of added Ca and Mg. Metal partitioning (K_d) in soil was most influenced by soil pH rather than organic matter (OM) showing that metal availability increases as liming effects fade. This study shows that the legacy effects of soil liming can persist for several decades and are most apparent in the forest floor (FH), but legacy effects are quite modest, and it is likely that a considerable amount of limestone has been lost through erosion.

Simple models efficiently predict free cadmium Cd2+ in the solutions of low-contaminated agricultural soils

Speciation of Cd in soil solutions strongly determines the fate of this toxic metal in the environment. Generally, in soil solutions, Cd predominantly binds to the dissolved organic matter (DOM). The determination of the quantity and reactivity of DOM that actually complexes Cd in soil solutions is challenging for operational purposes. Therefore, this study tested whether Cd²⁺ concentration in soil solutions could effectively be predicted by considering complexation with a single mean organic ligand assumed to be a fraction of DOM of unspecified nature or assumed to be purely fulvic acids (FA) with reactivity as described in WHAM VII. The reactivity of the unspecified ligand and the concentration of FA were modelled and fitted to experimental data from 76 agricultural soils with low Cd contents. The optimal reactivity and FA concentration that minimized the relative error (RE) of predictions of the concentration of Cd²⁺ in soil solutions were either considered constant across soils or modelled from soil properties by multiple linear regressions (MLR) or random forests (RF), giving 6 models, the predictive value of which was assessed by 10-folds cross-validation. When the reactivity of the mean ligand and the optimal FA concentration were considered constant across soils, the models were biased and 66.9% of predictions had relative errors below a factor of 2. By contrast, if the reactivity of the mean ligand or the optimal FA concentration were allowed to vary with soil characteristics, these performances increased to 95.5%, soil pH being the main predictor and RF being slightly more efficient than MLR. With more than 95% of the relative errors of prediction below a factor of 2, the models developed in this work could be valuable for assessing Cd speciation in the solution of soils having a low Cd content.

Proposal of new distribution coefficients (Kd) of potentially toxic elements in soils for improving environmental risk assessment in the State of São Paulo, southeastern Brazil

Soil solid-solution distribution coefficients (K_d) are used in predictive environmental models to assess public health risks. This study was undertaken to determine K_d for potentially toxic elements (PTE) Cd, Co, Cr, Cu, Ni, Pb, and Zn in topsoil samples (0-20 cm) from 30 soils in the State of São Paulo, southeastern Brazil. Batch sorption experiments were carried out, and PTE concentrations in the equilibrium solution were determined by High Resolution Inductively Coupled Plasma Mass Spectrometry (HR-ICPMS). Sorption data was fitted to the Freundlich model. The K_d values were either obtained directly from the slope coefficients of C-type isotherms or derived from the slope of the straight line tangent to the non-linear L-type and H-type isotherms. Stepwise multiple regression models were used to estimate the K_d values through the combined effect of a number of soil attributes [pH_H2O, effective cation exchange capacity (ECEC) and contents of clay, organic carbon, and Fe (oxy)hydroxides]. The smallest variation in K_d values was recorded for Cu (105-4598 L kg^-1), Pb (121-7020 L kg^-1), Ni (6-998 L kg^-1), as variation across four orders of magnitude was observed for Cd (7-14,339 L kg^-1), Co (2-34,473 L kg^-1), and Cr (1-21,267 L kg^-1). The K_d values for Zn were between 5 and 123,849 L kg^-1. According to median values of K_d, PTE were sorbed in the following preferential order: Pb > Cu > Cd > Ni > Zn > Cr > Co. The K_d values were best predicted using metal-specific and highly significant (p < 0.001) linear regressions that included pH_H2O, ECEC, and clay contents. The K_d values reported in this study are a novel result that can help minimize erroneous estimates and improve both environmental and public health risk assessments under humid tropical edaphoclimatic conditions.

Interactions of arsenic, copper, and zinc in soil-plant system: Partition, uptake and phytotoxicity

Arsenic, copper, and zinc are common elements found in contaminated soils but little is known about their combined effects on plants when presented simultaneously. Here, we systematically investigated the phytotoxicity and uptake of binary and ternary mixtures of As, Cu, and Zn in a soil-plant system, using wheat (Triticum aestivum) as model species. The reference models of concentration addition (CA) and response addition (RA) coupled with different expressions of exposure (total concentrations in soil ([M]_tot, mg/kg), free ion activities in soil solution ({M}, μM), and internal concentrations in plant roots ([M]_int, μg/g)), were selected to assess the interaction mechanisms of binary mixtures of AsCu, AsZn, and CuZn. Metal(loid) interactions in soil were estimated in terms of solution-solid partitioning, root uptake, and root elongation effects. The partitioning of one metal(loid) between the soil solution and solid phase was most often inhibited by the presence of the other metal(loid). In terms of uptake, inhibitory effects and no effects were observed in the mixtures of As, Cu, and Zn, depending on the mixture combinations and the dose metrics used. In terms of toxicity, simple (antagonistic or synergistic) and more complex (dose ratio-dependent or dose level-dependent) interaction patterns of binary mixtures occurred, depending on the dose metrics selected and the reference models used. For ternary mixtures (As-Cu-Zn), nearly additive effects were observed irrespective of dose descriptors and reference models. The observed interactions in this study may help to understand and predict the joint toxicity of metal(loid)s mixtures in soil-plant system. Mixture interactions and bioavailability should be incorporated into the regulatory framework for accurate risk assessment of multimetal-contaminated sites.

Evaluation of single-extraction methods to estimate the oral bioaccessibility of metal(loid)s in soils

Incidental ingestion of polluted soil particles exposes the population to toxic metal(loid)s. To refine the methods of exposure and risk assessment, it is relevant to use bioaccessible concentrations of metal(loid)s determined via in vitro digestion methods. However, some validated methods are complex and costly, involving high technical skills and numerous reagents. The objective of the present study was to evaluate the suitability of four simple chemical extractions to mimic the bioaccessible fraction of As, Cd, and Pb in the gastric (G) and gastrointestinal (GI) phases obtained using the validated UBM (unified bioaccessibility method) test. Acetic acid (0.11 M), citric acid (0.11 M), EDTA (0.16 M), and hydrochloric acid (HCl, 0.65%) were separately tested in 201 soil samples with a wide range of physicochemical parameters and metal(loid)s concentrations. Significant linear relationships were observed with HCl, EDTA, and to a lesser extent with citric acid. For the cheaper HCl method, correlations with the UBM ranged from 0.91 to 0.99 for the G phase and from 0.72 to 0.97 for the GI phase. This test can be used at least as a first-tier screening to assess the oral bioaccessibility of As, Cd, and Pb.

Single metal and metal mixture toxicity of five metals to Oppia nitens in five different Canadian soils

Metal mixture toxicity across soil types is a daunting challenge to risk assessment. Here, we evaluated metal mixture toxicity in Oppia nitens, using ten fixed metal mixture ratios in five Canadian soils that closely matched some of the EU PNEC reference soils. Soils were dosed with five metals (Cu, Zn, Pb, Co, Ni) as single metals (ten concentrations) and as mixtures (eight concentrations). Synchronized adult mites were exposed to metals, with survival and reproduction assessed after 28 days. We found out that (i) the differences among soils in mite sensitivity and single metals were not consistent when mites were exposed to metal mixtures, (ii) assuming concentration addition, the mixture interaction factor (MIF) showed that single metal low effect levels excessively underestimated low level metal mixture effects (iii) Zn emerged as a protective metal in most mixtures, and (iv) Soil properties such as CEC, independent of effects on metal speciation, explained more of the variation than measured metals. This study suggests that metal risk assessment should be done on a case by case basis. Further work is needed to ensure that by protecting soil-dwelling organisms from single metals, the risk from metal mixtures is appropriately protected for.

Review: mine tailings in an African tropical environment-mechanisms for the bioavailability of heavy metals in soils

Heavy metals are of environmental significance due to their effect on human health and the ecosystem. One of the major exposure pathways of Heavy metals for humans is through food crops. It is postulated in the literature that when crops are grown in soils which have excessive concentrations of heavy metals, they may absorb elevated levels of these elements thereby endangering consumers. However, due to land scarcity, especially in urban areas of Africa, potentially contaminated land around industrial dumps such as tailings is cultivated with food crops. The lack of regulation for land-usage on or near to mine tailings has not helped this situation. Moreover, most countries in tropical Africa have not defined guideline values for heavy metals in soils for various land uses, and even where such limits exist, they are based on total soil concentrations. However, the risk of uptake of heavy metals by crops or any soil organisms is determined by the bioavailable portion and not the total soil concentration. Therefore, defining bioavailable levels of heavy metals becomes very important in HM risk assessment, but methods used must be specific for particular soil types depending on the dominant sorption phases. Geochemical speciation modelling has proved to be a valuable tool in risk assessment of heavy metal-contaminated soils. Among the notable ones is WHAM (Windermere Humic Aqueous Model). But just like most other geochemical models, it was developed and adapted on temperate soils, and because major controlling variables in soils such as SOM, temperature, redox potential and mineralogy differ between temperate and tropical soils, its predictions on tropical soils may be poor. Validation and adaptation of such models for tropical soils are thus imperative before such they can be used. The latest versions (VI and VII) of WHAM are among the few that consider binding to all major binding phases. WHAM VI and VII are assemblages of three sub-models which describe binding to organic matter, (hydr)oxides of Fe, Al and Mn and clays. They predict free ion concentration, total dissolved ion concentration and organic and inorganic metal ion complexes, in soils, which are all important components for bioavailability and leaching to groundwater ways. Both WHAM VI and VII have been applied in a good number of soils studies with reported promising results. However, all these studies have been on temperate soils and have not been tried on any typical tropical soils. Nonetheless, since WHAM VII considers binding to all major binding phases, including those which are dominant in tropical soils, it would be a valuable tool in risk assessment of heavy metals in tropical soils. A discussion of the contamination of soils with heavy metals, their subsequent bioavailability to crops that are grown in these soils and the methods used to determine various bioavailable phases of heavy metals are presented in this review, with an emphasis on prospective modelling techniques for tropical soils.

Prediction of free metal ion activity in contaminated soils using WHAM VII, baker soil test and solubility model

Bioavailability and ecotoxicity of metals in contaminated soils depend largely on their solubility. The present investigation was carried out to predict the free ion activity of Zn²⁺, Cu²⁺, Ni²⁺, Pb²⁺ and Cd²⁺ in contaminated soils as a function of pH, organic carbon content and extractable metal concentration. Twenty-five composite soil samples were collected from various locations which had a history of receiving sewage sludge (Keshopur and IARI, Delhi), municipal solid waste (Kolkata, West Bengal), polluted river water (Madanpur, Delhi) and industrial effluents (Debari, Rajasthan and Sonepat, Haryana). Four composite soil samples were also collected from adjacent fields which had not received contaminated amendments. Free ion activities (-log₁₀ values), viz. pZn²⁺, pCu²⁺, pNi²⁺, pPb²⁺ and pCd²⁺ as measured by the Baker soil test, were 10.1 ± 1.12, 13.4 ± 1.23, 12.9 ± 0.85, 11.6 ± 0.74 and 12.6 ± 2.26, respectively. Free metal ion activities were also determined using the geochemical speciation model WHAM-VII following extraction of soil solution with porous Rhizon samplers from the rhizosphere of growing plants. pH dependent Freundlich model based on soil properties could explain the variation in pZn²⁺, pCu²⁺, pNi²⁺, pPb²⁺ and pCd²⁺ to the extent of 84, 52, 73, 60 and 70%, respectively, in the case of data from Rhizon samplers coupled with speciation modelling. Whereas, C-Q model could explain 84, 57, 82, 72 and 74% variability in pZn²⁺, pCu²⁺, pNi²⁺, pPb²⁺ and pCd²⁺, respectively, based on soil properties and free metal ion activity as determined with integrated use of Rhizon-WHAM-VII. Modelling approach was superior compared to that based on the Baker soil test solution.

Lead was mobilized in acid silty clay loam paddy soil with potassium dihydrogen phosphate (KDP) amendment

The immobilization effectiveness between Pb and phosphorus in soil varies with soil types. To clarify the effect of phosphate on the availability of Pb in agricultural soil, a culture experiment with three types of paddy soil was performed with potassium dihydrogen phosphate (KDP) added. EDTA, DGT and in-situ solution extraction methods were used to represent different available Pb content. Results showed that the concentration of EDTA-Pb in HN soil was slightly elevated after exogenous KDP added. The supplement of 300 mg/kg KDP significantly increased the content of soluble Pb in both acid silty clay loam soil and neutral silty loam soil (increased by 104.65% and 65.12%, respectively). However, there was no significant influence of KDP on the concentration of DGT extracted Pb. XANES results showed that Pb(OH)₂, PbHPO₄, humic acid-Pb and GSH-Pb were the major speciation of Pb in soil colloids. The proportion of Pb(OH)₂ and humic acid-bounded Pb in soil colloids were elevated after exogenous KDP added. Our results indicated that there was a mobilization effect of KDP on Pb by increasing the amount of colloidal Pb in soil solution, especially in acid silty clay loam paddy soil. Such colloid-facilitated transport might promote the uptake of Pb in rice and pose a potential threat to human health.

Effects of Soil Properties on the Toxicity and Bioaccumulation of Lead in Soil Invertebrates

The present study examined the effects of soil physical and chemical properties on the toxicity of lead (Pb) to earthworms (Eisenia fetida) and collembolans (Folsomia candida), and on bioaccumulation of Pb by earthworms, in soils amended with Pb salts. Toxicity tests were conducted in 7 soils varying in soil properties (pH 4.7-7.4, effective cation exchange capacity [eCEC] 4-42 cmol_c /kg, organic carbon 10-50 g C/kg) that were leached and pH corrected after spiking with PbCl₂ . The median effect concentrations (EC50s) based on total soil Pb concentrations ranged from 35 to 5080 mg Pb/kg for earthworms and 389 to >7190 mg/kg for Collembola. Significant positive correlations were observed between log (EC50) for earthworm reproduction and log (eCEC, total C, exchangeable Ca and Mg, or clay content), but no significant correlations were observed between Pb toxicity to Collembola and soil properties. Expressing Pb dose as either the free ion (Pb²⁺ ) activity in porewater or as the measured dissolved porewater concentration of Pb did not explain differences in toxicity among soils. The bioaccumulation factors (BAFs) for Pb in earthworms ranged up to >10-fold across 6 soil treatments, with a median of 0.16, and the BAF was significantly correlated with eCEC (p = 0.038, r = -0.84), but not with any other soil properties. Soil properties related to eCEC (total C, exchangeable Ca and Mg, clay content) had a significant effect on Pb toxicity and bioaccumulation in earthworms, but no relationship was found for Collembola. As a major soil property affecting the bioavailability of Pb, CEC should be incorporated into any soil hazard assessment of Pb as a modifying factor of toxicity and bioaccumulation for earthworms. Environ Toxicol Chem 2019;38:1486-1494. © 2019 SETAC.

The forgotten role of toxicodynamics: How habitat quality alters the mite, Oppia nitens, susceptibility to zinc, independent of toxicokinetics

Soil habitat quality is thought to influence metal toxicity via changes in speciation and thereby toxicokinetics. Here, we assessed the toxicokinetic and toxicodynamic effects of habitat quality on mite, Oppia nitens when exposed to zinc (Zn) contaminated soils. Forty-seven soils were ranked into three habitat qualities; high, medium, and low based on biological reproduction of Folsomia candida, Enchytraeus crypticus, and Elymus lanceolatus. From the 47 soils, eighteen soils (comprising of six soils from each habitat quality) were randomly selected and dosed with field relevant concentrations of Zn. Mite survival and reproduction were assessed after 28 days. Total Zn, bioaccessible Zn, Zn bioavailability, Zn body burden, lactate dehydrogenase activity (LDH) and glucose-6-phosphate dehydrogenase (G6PDH) activities of the mites were determined. Zinc toxicity and potency were much less in the high compared to low quality soils and the mites in the high habitat quality soils tolerated higher zinc body burdens (2040 ± 130 μg/g b.w) than the lower habitat quality (1180 ± 310 μg/g b.w). Lower LDH activity (20 ± 2 μU mg^-1) in the high quality soils compared to lower quality soils (50 ± 8 μU mg^-1) suggested that there was less stress in the high habitat quality mites. Despite changes in speciation across habitat qualities, bioavailability of zinc was similar (∼20%) irrespective of habitat quality. Our results suggest that the influence of soil properties on survival is modulated by toxicodynamics rather than toxicokinetics. Restoring habitat quality may be more important for soil invertebrate protection than metal concentration at contaminated sites.

Increased Methylmercury Accumulation in Rice after Straw Amendment

Consumption of rice has been shown to be an important route of dietary exposure to methylmercury (MeHg, a neurotoxin) for Asians having a low fish but high rice diet. Therefore, factors that increase MeHg production and bioaccumulation in soil-rice systems, could enhance the risk of MeHg exposure. On the basis of a national-scale survey in China (64 sites in 12 provinces) and rice cultivation experiments, we report that straw amendment, a globally prevalent farming practice, could increase MeHg concentrations in paddy soils (11-1043%) and rice grains (95%). By carrying out a series of batch incubation, seedling uptake and sand culture experiments, we demonstrate that these increases could be attributed to (1) enhanced abundances/activities of microbial methylators and the transformation of refractory HgS to organic matter-complexed Hg, facilitating microbial Hg methylation in soils; (2) enhanced MeHg mobility, and increased root lengths (35-41%) and tip numbers (60-105%), increasing MeHg uptake by rice roots; and (3) enhanced MeHg translocation to rice grains from other tissues. Results of this study emphasize fresh organic matter-enhanced MeHg production and bioaccumulation, and highlight the increased risk of MeHg after straw amendment and thus the need for new policies concerning straw management.

A Comparison between the Oil Removal Capacity of Polymer-Coated Magnetic Nanoparticles in Natural and Synthetic Environmental Samples

Oil spills can have dramatic impacts on the environment. The limitations of current oil remediation techniques have inspired researchers to study the application of nanotechnology for oil cleanup. Previously, we reported essentially 100% removal of a reference MC252 oil using polyvinylpyrrolidone (PVP)-coated iron oxide nanoparticles (NPs) from oil-water mixtures under a wide range of environmentally relevant conditions. Our previous results showed that in the synthetic water samples, the concentration of cations and natural organic macromolecules (NOM) can significantly affect the oil removal efficiency of NPs. Here, we studied the application of these NPs for oil removal from natural freshwater samples and compared the results with the synthetic water samples with the same concentrations of major ions and NOM. For both natural and synthetic samples, concentrations of NOM, calcium, and magnesium were positively correlated with oil removal ( p-value <0.05). NPs show an average of 30% higher oil removal efficiency from natural samples compared to the synthetic samples. Using up to 50 ppm of NP, essentially 100% oil removal was observed under most conditions specially hardwater samples (initial 0.15 g L^-1 oil concentration). Results show that these NPs are a facile and reliable technique for removing oil under realistic conditions.

Trace metal availability in soil horizons amended with various urban waste composts during 17 years - Monitoring and modelling

Recycling organic residues in agrosystems presents several benefits but faces the question of contaminants, among them a few trace metals which eventually accumulate in soils following regular applications of organic waste products (OWP) and represent an ecological risk. The increase of total trace metal contents in amended topsoils can be predicted by a mass balance approach, but the evolution of their available fractions is a more intricate issue. We aimed at modelling this evolution by using the dataset of a long-term field experiment of OWP applications (manure and three urban waste composts). Two operationally-defined fractions of 6 trace metals have been quantified in the OWP and amended topsoils between 2002 and 2015: the soluble and potentially available metals, extracted in 0.01 M CaCl₂ and 0.05 M EDTA solutions, respectively. The potentially available metals have progressively increased in amended topsoils, at rates depending on elements and types of OWP. For Zn, these increases corresponded in average to inputs of potentially available Zn from OWP. But the soil stocks of potentially available Cu increased faster than from the inputs of EDTA-extractable Cu, showing linear regression slopes between 1.4 and 2.5, depending on OWP type. The influence of OWP has been provisionally interpreted in the light of their efficiency to increase soil organic matter and their inputs of reactive oxides. Soluble copper has increased with repeated amendments. But soluble cadmium, nickel and zinc have generally decreased, as they are influenced by changing soil variables such as pH and organic matter. Statistic models were used to unravel the relationships between soluble and EDTA-extractable metals and other soil variables. For Cu, the most satisfactory models just relate soluble and potentially available Cu. Developing such models could contribute to predict the long-term effects of a precise scenario of agricultural OWP recycling upon available trace metals in soils.

Arsenic Distribution Assessment in a Residential Area Polluted with Mining Residues

Mining is a major source for metals and metalloids pollution, which could pose a risk for human health. In San Guillermo, Chihuahua, Mexico mining wastes are found adjacent to a residential area. A soil-surface sampling was performed, collecting 88 samples for arsenic determination by atomic absorption. Arsenic concentration data set was interpolated using the ArcGis models: inverse distance weighting (IDW), ordinary kriging (OK), and radial basis function (RBF). For method validation purposes, a set of the data was selected and two tests were performed (P1 and P2). In P1 the models were processed without the validation data; in P2 the validation data were removed one by one, models were processed every time that a data point was removed. An arsenic concentration range of 22.7 to 2190 mg/kg was reported. The 39% of data set was classified as contaminated soil and 61% as industrial land use. In P1 the method of interpolation with the lowest RMSE was RBF (0.80), the highest coefficient of E was RBF (46.25), and the highest Ceff value was with RBF (0.48). In P2 the method with the lowest RMSE was OK (0.76), the highest E value was 50.65 with OK, and the Ceff reported the highest value with OK (0.52). The high arsenic contamination in soil of the site indicates an abundant dispersion of this metalloid. Furthermore, the difference between the models was not very wide. The incorporation of more parameters would be of interest to observe the behavior of interpolation methods.

Geochemical cadmium anomaly and bioaccumulation of cadmium and lead by rapeseed (Brassica napus L.) from noncalcareous soils in the Guizhou Plateau

The cruciferous crop, oil rapeseed (Brassica napus L.), may bioaccumulate excessive cadmium (Cd) and lead (Pb) as well, from Cd-enriched noncalcareous soils in Guizhou province of southwestern China. Field paired soil-rapeseed sampling and greenhouse experiment were performed to characterize the Cd anomaly in rapeseed-planting soils and to predict the bioaccumulation of Cd and Pb in raw seeds using soil variables. The results indicated that total soil-Cd concentrations averaged 0.43 mg kg^-1 (range 0.11-1.41 mg kg^-1) from field investigation; and a soil type dependent Cd anomaly was observed. Besides, cumulative frequency of total soil-Cd was plotted to be helpful in delimitation of regional Cd anomalies. Rapeseeds readily bioaccumulated Cd from soils as validated by greenhouse experiment and field data. Contrary results were observed in relation to rapeseed Pb levels measured from greenhouse experiment (very low) and field (very high) which was likely due to soil particle contamination as indicated by the considerably higher ratio of Pb to Cd level in seeds harvested from fields. Based on multiple stepwise regression analysis, reliable Cd soil-rapeseed relationships, but less reliable for Pb, were derived using either total or (bio)available metal concentrations and were further inversely used to derive local soil Cd criteria (e.g., total soil-Cd based, 2.5 mg kg^-1) based on Hygienic Standard for Feeds (GB13082-2001). Although seed Cd levels (<feed standard) observed in field data indicated a least human dietary risk, however, high (bio)availability of Cd, but not Pb, in Cd-enriched acid soil still poses high environmental risks and may threaten food safety of other crops.

Influence of electrode placement for mobilising and removing metals during electrodialytic remediation of metals from shooting range soil

Electrodialytic remediation was applied to a shooting range soil to investigate the influence of electrode placement on the removal and binding of metals during the treatment. The set-up was based on a 2-compartment cell, in which the cathode was separated from the soil by a cation exchange membrane and the anode was placed directly in the soil, thereby introducing protons and oxygen directly in the soil. Mobilisation of metals from less available fractions (oxidisable and residual) in the soil occurred, due to oxidation/dissolution of insoluble/soluble organic matter and possibly metal oxides in the residual fraction. The transport via electromigration out of the soil and/or re-precipitation in other fractions of the soil (oxidisable, reducible, exchangeable) depended on the metal. More than 30% of the initial content of Mn, Cd, Cu, Pb and Zn and less than 20% of the initial content of Al, Fe, K, Mg, As, Cr and Ni was transported out of the soil. By decreasing the distance between the electrodes from 3.0 to 1.5 cm, the removal of the targeted metal for remediation, Pb, was improved by more than 200%, from 14 to 31%. A similar removal could be achieved in experiments with long distance between electrodes (3.0 cm) by increasing the current intensity from 4 to 10 mA and/or the remediation time from 7 to 35 d. The experiments showed that the design and optimisation of electrodialytic remediation depends on the targeted metal and metal partitioning.

Soil pH effects on the toxicity of zinc oxide nanoparticles to soil microbial community

We conducted an experiment with two agricultural soils with acidic and alkaline pH levels to assess the effects of zinc oxide nanoparticles (nZnO) on the bacterial community. The effect of the nZnO concentrations (0, 0.1, 1, 10, 100, 1000 mg Zn/kg soil) and contact time between nanoparticles and soil (180 days) was considered. We measured the microbial respiration rate, nitrogen transformation, enzymatic activities (dehydrogenase (DH), acidic phosphatase (ACP), and alkaline phosphatase (ALP)), and the community-level physiological profile (CLPP) soil parameters. Respiration potential and nitrogen transformation were negatively affected only at the highest nZnO concentration. The changes in enzymatic activities were very variable with time and between both soils. A stimulating effect of the nanoparticles on microbial activity was clearly shown at 30 days after the nZnO application in both soils, except for the 1000 mg/kg in calcareous soil, after which time in the latter, the functional richness of the bacterial community was reduced to virtually zero. However, values of the enzymatic activities demonstrated that there was self-adaptation of microbial communities over the study period (180 days). The nZnO 1000 mg/kg dose produced an increase in bacterial growth in the acidic soil, without apparent changes in their metabolic profiles over time. A good correlation was found between microbial respiration rates (calcareous and acidic soils) and microbial metabolic activity (acidic soil) based on extracted Zn concentrations. Our findings suggest the necessity of additional studies to examine the effects of nZnO in natural microorganism populations in soil with different pH levels for extended periods of time.

Controls on accumulation and soil solution partitioning of heavy metals across upland sites in United Kingdom (UK)

A significant body of knowledge suggests that soil solution pH and dissolved organic carbon (DOC) strongly influence metal concentrations and speciation in porewater, however, these effects vary between different metals. This study investigated the factors influencing soil and soil solution concentrations of copper (Cu), lead (Pb), nickel (Ni) and zinc (Zn) under field conditions in upland soils from UK having a wide range of pH, DOC and organic matter contents. The study primarily focussed on predicting soil and soil solution metal concentrations from the data on total soil metal concentrations (HNO₃ extracts) and soil and soil solution properties (pH, DOC and organic matter content). We tested the multiple regression models proposed by Tipping et al. (2003) to predict heavy metal concentrations in soil solutions and the results indicated a better fit (higher R² values) in both studies for Pb compared to the Zn and Cu concentrations. Both studies observed consistent negative relationships of metals with pH and loss on ignition (LOI) suggesting an increase in soil solution metal concentrations with increasing acidity. The positive relationship between Pb concentrations in porewater and HNO₃ extracts was similar for both studies, however, similar relationships were not found for the Zn and Cu concentrations because of the negative coefficients for these metals in our study. The results of this study conclude that the predictive equations of Tipping et al. (2003) may not be applicable to the field sites where the range of DOC and metal concentrations is much lower than their study. Our study also suggests that the extent to which metals are partitioned into soil solution is lower in soils with a higher organic matter contents due to binding of these metals to soil organic matter.

Proper management of lead-contaminated agricultural lands against the exceedance of lead in agricultural produce: Derivation of local soil criteria

The Measures for Management of Soil Environment in Agricultural Land (Trial, Nov. 01, 2017, China) recently came into effect and highlighted the proper management of contaminated croplands to lower risks of exceedances of contaminants, especially toxic trace metals in agricultural produce. We aimed to develop local soil criteria for lead (Pb) in Hezhang county of southwestern China by the inverse use of reliable models linking Pb contamination levels between soils and vegetables. Dilute nitric acid (0.43 M) extraction, a new ISO standard (ISO-17586:2016) for extracting the geochemically reactive Pb fraction (Pb_NA), and calcium chloride (0.01 M) extraction (ISO-14255: 1998) for estimating the plant-available Pb (Pb_CC) were performed in fifty historically polluted and newly Pb-spiked soils with differing soil types, properties (pH 4.1-8.0), and total soil Pb levels (Pb_T, 20-6153 mg kg^-1). Greenhouse experiments for Brassica pekinensis L., and in-situ soil porewater measurement for Pb were conducted to investigate the mechanism of Pb uptake, and to establish reliable Pb soil-plant relationships. The results indicated that about 83% of the variation for Pb concentrations in vegetable (Pb_CL, 0.009-1.06 mg kg^-1) was contributable to free Pb²⁺ activity in soil porewater, which was mainly influenced by pH and dissolved organic matter. Pb_CL was satisfactorily predicted using Pb_NA and key soil properties (adj. R² 0.852). Soil Pb criteria for Pb_T and Pb_NA are then derived based on food standard. The full implementation of criteria derived for Pb_NA (i.e., 27-127 mg kg^-1, soil pH 5.5-8.0) can avoid the exceedance of Pb in 95% of cabbage samples in this study, 95% of cabbage cultivars by model extrapolation, and one widely cultivated root vegetable, radish, in the study region. We provide a successful case study that has effectively tackled the challenge for the complexity of the soil management in contaminated croplands.

Assessing bioavailability of complex chemical mixtures in contaminated soils: Progress made and research needs

Understanding the distribution, behaviour and interactions of complex chemical mixtures is key for providing the evidence necessary to make informed decisions and implement robust remediation strategies. Much of the current risk assessment frameworks applied to manage land contamination are based on total contaminant concentrations and the exposure assessments embedded within them do not explicitly address the partitioning and bioavailability of chemical mixtures. These oversights may contribute to an overestimation of both the eco-toxicological effects of the fractions and the mobility of contaminants. In turn, this may limit the efficacy of risk frameworks to inform targeted and proportionate remediation strategies. In this review we analyse the science surrounding bioavailability, its regulatory inclusion and the challenges of incorporating bioavailability in decision making process. While a number of physical and chemical techniques have proven to be valuable tools for estimating bioavailability of organic and inorganic contaminants in soils, doubts have been cast on its implementation into risk management soil frameworks mainly due to a general disagreement on the interchangeable use of bioavailability and bioaccessibility, and the associated methods which are still not standardised. This review focuses on the role of biotic and abiotic factors affecting bioavailability along with soil physicochemical properties and contaminant composition. We also included advantages and disadvantages of different extraction techniques and their implications for bioavailability quantitative estimation. In order to move forward the integration of bioavailability into site-specific risk assessments we should (1) account for soil and contaminant physicochemical characteristics and their effect on bioavailability; (2) evaluate receptor's potential exposure and uptake based on mild-extraction; (3) adopt a combined approach where chemical-techniques are used along with biological methods; (4) consider a simplified and cost-effective methodology to apply at regulatory and industry setting; (5) use single-contaminant exposure assessments to inform and predict complex chemical mixture behaviour and bioavailability.

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⁺, Ca²⁺, and Mg²⁺ to Cu²⁺ toxicity were studied individually and independently, and linear relationships between EC50 free Cu²⁺ and H⁺, Ca²⁺, and Mg²⁺ activities were found, meaning that the cations H⁺, Ca²⁺, and Mg²⁺ will alleviate the toxicity of Cu²⁺ in solutions. Toxicity accompanying increasing concentration of solution ions other than Cu²⁺ 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.

Evaluating water quality and ecotoxicology assessment techniques using data from a lead and zinc effected upland limestone catchment

Point and diffuse sources associated with historical metal ore mining are major causes of metal pollution. The understanding of metal behaviour and fate has been improved by the integration of water chemistry, metal availability and toxicity. Efforts have been devoted to the development of efficient methods of assessing and managing the risk posed by metals to aquatic life and meeting national water quality standards. This study focuses on the evaluation of current water quality and ecotoxicology techniques for the metal assessment of an upland limestone catchment located within a historical metal (lead ore) mining area in northern England. Within this catchment, metal toxicity occurs at circumneutral pH (6.2-7.5). Environmental Quality Standards (EQSs) based on a simple single concentration approach like hardness based EQS (EQS-H) are more overprotective, and from sixteen sites monitored in this study more than twelve sites (>75%) failed the EQSs for Zn and Pb. By increasing the complexity of assessment tools (e.g. bioavailability-based (EQS-B) and WHAM-F_TOX), less conservative limits were provided, decreasing the number of sites with predicted ecological risk to seven (44%). Thus, this research supports the use of bioavailability-based approaches and their applicability for future metal risk assessments.

Predicting trace metal solubility and fractionation in Urban soils from isotopic exchangeability

Metal-salt amended soils (MA, n = 23), and historically-contaminated urban soils from two English cities (Urban, n = 50), were investigated to assess the effects of soil properties and contaminant source on metal lability and solubility. A stable isotope dilution method, with and without a resin purification step, was used to measure the lability of Cd, Cu, Ni, Pb and Zn. For all five metals in MA soils, lability (%E-values) could be reasonably well predicted from soil pH value with a simple logistic equation. However, there was evidence of continuing time-dependent fixation of Cd and Zn in the MA soils, following more than a decade of storage under air-dried conditions, mainly in high pH soils. All five metals in MA soils remained much more labile than in Urban soils, strongly indicating an effect of contaminant source on metal lability in the latter. Metal solubility was predicted for both sets of soil by the geochemical speciation model WHAM-VII, using E-value as an input variable. For soils with low metal solution concentrations, over-estimation of Cd, Ni and Zn solubility was associated with binding to the Fe oxide fraction while accurate prediction of Cu solubility was dependent on humic acid content. Lead solubility was most poorly described, especially in the Urban soils. Generally, slightly poorer estimation of metal solubility was observed in Urban soils, possibly due to a greater incidence of high pH values. The use of isotopically exchangeable metal to predict solubility is appropriate both for historically contaminated soils and where amendment with soluble forms of metal is used, as in toxicological trials. However, the major limitation to predicting solubility may lie with the accuracy of model input variables such as humic acid and Fe oxide contents where there is often a reliance on relatively crude analytical estimations of these variables. Trace metal reactivity in urban soils depends on both soil properties and the original source material; the WHAM geochemical model predicts solubility using isotopically exchangeable metal as an input.

Inhibitory effects of rice residues amendment on Cd phytoavailability: A matter of Cd-organic matter interactions?

Incorporating crop residues into soils, a most common way of organic input into farmland soils, is being encouraged in many parts of the world, while its potential impacts on Cd phytoavailability are not well understood. Here, a Cd-contaminated soil was amended with rice residues (RR, i.e., straw + root mixture) or not (Control) and incubated for 81 days under laboratory-controlled conditions. During the incubation, key soil parameters (e.g., dissolved organic carbon and pH), Cd geochemical fractionation (by BCR sequential extraction), and CaCl₂ extracted Cd in soils (by 0.01 M CaCl₂ extraction) were quantified to explain the effects of RR amendment on Cd phytoavailability (assessed by 7 day-cultivation of rice seedling in soils). Besides, hydroponic experiments were designed to explore the effects of D-RR-OM (dissolved-RR-organic matter) on the uptake of Cd by rice seedlings (quantified by uptake constant rate, k_u, using stable isotope tracing technique). Our results demonstrated that RR amendment reduced Cd phytoavailability by 17-92% compared with Control during incubation, which might be explained by the interactions between Cd and RR-OM (RR-organic matter) in soil or porewater: (1) Cd immobilization due to its association with solid-RR-OM in soils, (2) Cd mobilization by D-RR-OM, and (3) Cd complexation with D-RR-OM in porewater, and thus reduced k_u of Cd. Our results suggested that dynamics of RR-OM (e.g., dissolution, decomposition and transformation) in soils, and thus interactions between Cd and solid/dissolved-RR-OM may control Cd phytoavailability under RR amendment. Information gained in this study would further our understanding about Cd phytoavailability in farming soils.

Experimentally derived acute and chronic copper Biotic Ligand Models for rainbow trout

We evaluated the effects of varying water chemistry ([Ca²⁺]=0.2-3mM, [Mg²⁺]=0.05-3mM, dissolved organic matter (DOM, natural, from maple leaves)=0.3-10mg of CL^-1, pH=5.0-8.5) on the acute (96-h, unfed fish) and chronic (30-d, fed fish) toxicity of waterborne Cu to juvenile rainbow trout (Oncorhynchus mykiss) exposed in flow-through conditions. Acute and chronic Biotic Ligand Models (BLMs) were developed from the obtained toxicity data-sets, using the Visual MINTEQ software. Our results indicate that Cu is predominantly an acute toxicant to rainbow trout, as there were no observable growth effects and the 96-h and 30-d LC50 values were similar, with mortality mostly occurring within the first few days of exposure. Calcium and DOM were greatly protective against both acute and chronic Cu toxicity, but Mg seemed to only protect against chronic toxicity. Additional protection by pH 5.0 in acute exposure and by pH 8.5 in chronic exposure occurred. In the range of conditions tested, the observed 96-h LC50 and 30-d LC20 values varied by a factor of 39 and 27 respectively. The newly developed acute and chronic BLMs explained these variations reasonably well (i.e. within a 2-fold error), except at pH≥8 where the high observed acute toxicity could not be explained, even by considering an equal contribution of CuOH⁺ and Cu²⁺ to the overall Cu toxicity. The 96-h LC50 values of 59% of 90 toxicity tests from 19 independent studies in the literature were reasonably well predicted by the new acute BLM. The LC20 predictions from the new chronic BLM were reasonable for 7 out of 14 toxicity tests from 6 independent chronic studies (with variable exposure durations). The observed deviations from BLM predictions may be due to uncertainties in the water chemistry in these literature studies and/or to differences in fish sensitivity. A residual pH effect was also observed for both the acute and the chronic data-sets, as the ratio of predicted vs. observed LC values generally increased with the pH. Additional mechanistic studies are required to understand the influence of pH, Na, and Mg on Cu toxicity to trout. The present study presents the first experimentally developed chronic Cu BLM for the rainbow trout. To the best of our knowledge, it also presents the first acute Cu BLM that is based on a published data-set for trout. These newly developed BLMs should contribute to improving the risk assessment of Cu to fish in freshwater.

Nanospecific Phytotoxicity of CuO Nanoparticles in Soils Disappeared When Bioavailability Factors Were Considered

Bioavailability-modifying factors such as soil type and aging have only rarely been considered in assessing toxicity of metal-containing nanoparticles in soil. Here, we examined the toxicity to barley (Hordeum vulgare) of CuO nanoparticles (CuO-NPs) relative to CuO bulk particles (CuO-BPs) and Cu acetate (Cu(OAc)₂) in six different soils with or without aging. The set up allows identifying whether or not NPs-derived colloidal Cu in soil porewater contributes to toxicity. Ultrafiltration (50 kDa) was performed together with geochemical modeling to determine {Cu²⁺} (free Cu²⁺ activity in soil porewater). Based on total soil Cu concentration, toxicity measured with seedling root elongation ranked Cu(OAc)₂ > CuO-NPs > CuO-BPs in freshly spiked soils. The differences in toxicity among the three toxicants became smaller in soils aged for 90 days. When expressing toxicity as {Cu²⁺}, there was no indication that nanoparticulate or colloidal Cu enhanced toxicity. A calibrated bioavailability-based model based on {Cu²⁺} and pH successfully explained (R² = 0.78, n = 215) toxicity of all Cu forms in different soils with and without aging. Our results suggest that toxicity predictions and risk assessment of CuO-NPs can be carried out properly using the bioavailability-based approaches that are used already for their non-nano counterparts in soil.

Characteristics and impacts of trace elements in atmospheric deposition at a high-elevation site, southern China

To investigate the regional background trace element (TE) level in atmospheric deposition (dry and wet), TEs (Fe, Al, V, Cr, Mn, Ni, Cu, Zn, As, Se, Mo, Cd, Ba, and Pb) in 52 rainwater samples and 73 total suspended particles (TSP) samples collected in Mt. Lushan, Southern China, were analyzed using inductively coupled plasma-mass spectrometry (ICP-MS). The results showed that TEs in wet and dry deposition of the target area were significantly elevated compared within and outside China and the volume weight mean pH of rainwater was 4.43. The relative contributions of wet and dry depositions of TEs vary significantly among elements. The wet deposition fluxes of V, As, Cr, Se, Zn, and Cd exceeded considerably their dry deposition fluxes while dry deposition dominated the removal of pollution elements such as Mo, Cu, Ni, Mn, and Al. The summed dry deposition flux was four times higher than the summed wet deposition flux. Prediction results based on a simple accumulation model found that the content of seven toxic elements (Cr, Ni, Cu, Zn, As, Cd, and Pb) in soils could increase rapidly due to the impact of annual atmospheric deposition, and the increasing amounts of them reached 0.063, 0.012, 0.026, 0.459, 0.076, 0.004, and 0.145 mg kg^-1, respectively. In addition, the annual increasing rates ranged from 0.05% (Cr and Ni) to 2.08% (Cd). It was also predicted that atmospheric deposition induced the accumulation of Cr and Cd in surface soils. Cd was the critical element with the greatest potential ecological risk among all the elements in atmospheric deposition.

Development of electrostatic-based bioavailability models for interpreting and predicting differential phytotoxicity and uptake of metal mixtures across different soils

Metals are ubiquitous and normally co-occur as mixtures in soil, but there remains much to do regarding the development of appropriate models which incorporate mixture interactions and bioavailability to estimate their phytotoxicity and phytoaccumulation. Here, we developed a probability-based electrostatic toxicity model (ETM) and a Langmuir-type electrostatic uptake model (EUM) to predict and normalize toxicity and uptake of zinc-copper mixtures in Hordeum vulgare L. in different soils. For model development, the electrical potential (ψ₀) and metal ion activities ({M²⁺}₀) at the cell-membrane surface was computed based on plant physiological properties and soil solution chemistry. Single metal toxicity correlated more closely to their corresponding {M²⁺}₀ than to ion activities in soil solution or total soil metal concentrations. The ETM explained up to 89% of the variance in mixture toxicity across different soils. Incorporation of ψ₀ into the EUM improved the model's ability for predicting metal uptake. Besides, cell-surface H⁺ appeared to significantly inhibit copper uptake via competition or other mechanisms, beyond its effect upon ψ₀. Our results for the first time demonstrate that electrostatic theory can be used to predict and reconcile mixture toxicity and uptake data in different soils, indicating the potential of electrostatic-based models in risk assessment of multimetal-contaminated soils.

Speciation and reactivity of lead and zinc in heavily and poorly contaminated soils: Stable isotope dilution, chemical extraction and model views

Correct characterization of metal speciation and reactivity is a prerequisite for the risk assessment and remedial activity management of contaminated soil. To better understand the intrinsic reactivity of Pb and Zn, nine heavily and poorly contaminated soils were investigated using the combined approaches of chemical extractions, multi-element stable isotopic dilution (ID) method, and multi-surface modelling. The ID results show that 0.1-38% of total Pb and 3-45% of total Zn in the studied soils are isotopically exchangeable after a 3-day equilibration. The intercomparison between experimental and modelling results evidences that single extraction with 0.43 M HNO₃ solubilizes part of non-isotopically exchangeable fraction of Pb and Zn in the studied soils, and cannot be used as a surrogate for ID to assess labile Pb and Zn pools in soil. Both selective sequential extraction (SSE) and modelling reveal that Mn oxides are the predominant sorption surface for Pb in the studied soils; while Zn is predicted to be mainly associated with soil organic matter in the soil with low pH and Fe/Mn oxides in the soils with high pH. Multi-surface modelling can provide a reasonable prediction of Pb and Zn adsorption onto different soil constituents for the most of the studied soils. The modelling could be a promising tool to decipher the underlying mechanism that controls metal reactivity in soil, but the submodel for Mn oxides should be incorporated and the model parameters, especially for the 2-pK diffuse layer model for Mn oxides, should be updated in the further studies.

Parameterizing the binding properties of dissolved organic matter with default values skews the prediction of copper solution speciation and ecotoxicity in soil

Parameterizing speciation models by setting the percentage of dissolved organic matter (DOM) that is reactive (% r-DOM) toward metal cations at a single 65% default value is very common in predictive ecotoxicology. The authors tested this practice by comparing the free copper activity (pCu²⁺  = -log₁₀ [Cu²⁺ ]) measured in 55 soil sample solutions with pCu²⁺ predicted with the Windermere humic aqueous model (WHAM) parameterized by default. Predictions of Cu toxicity to soil organisms based on measured or predicted pCu²⁺ were also compared. Default WHAM parameterization substantially skewed the prediction of measured pCu²⁺ by up to 2.7 pCu²⁺ units (root mean square residual = 0.75-1.3) and subsequently the prediction of Cu toxicity for microbial functions, invertebrates, and plants by up to 36%, 45%, and 59% (root mean square residuals ≤9 %, 11%, and 17%), respectively. Reparametrizing WHAM by optimizing the 2 DOM binding properties (i.e., % r-DOM and the Cu complexation constant) within a physically realistic value range much improved the prediction of measured pCu²⁺ (root mean square residual = 0.14-0.25). Accordingly, this WHAM parameterization successfully predicted Cu toxicity for microbial functions, invertebrates, and plants (root mean square residual ≤3.4%, 4.4%, and 5.8%, respectively). Thus, it is essential to account for the real heterogeneity in DOM binding properties for relatively accurate prediction of Cu speciation in soil solution and Cu toxic effects on soil organisms. Environ Toxicol Chem 2017;36:898-905. © 2016 SETAC.

Redistribution of elements between wastes and organic-bearing material in the dispersion train of gold-bearing sulfide tailings: Part I. Geochemistry and mineralogy

Migration and redistribution of elements during prolonged interaction of cyanide wastes with the underlying natural organic-bearing material have been studied in two ~40cm deep cores that sample primary ores and their weathering profile (wastes I and II, respectively) in the dispersion train of gold-bearing sulfide tailings in Siberia. Analytical results of SR-XRF, whole-rock XRF, AAS, CHNS, and SEM measurements of core samples show high K, Sr, Ti, and Fe enrichments and correlation of P₂O₅ and Mn with LOI and C_org. Organic material interlayered or mixed with the wastes accumulates Cu, Zn, Se, Cd, Ag, Au, and Hg. The peat that contacts wastes II bears up to 3wt.% Zn, 1000g/t Se, 100g/t Cd, and 8000g/t Hg. New phases of Zn and Hg sulfides and Hg selenides occur as abundant sheaths over bacterial cells suggesting microbial mediation in sorption of elements. Organic-bearing material in the cores contains 10-30g/t Au in 2-5cm thick intervals, both within and outside the intervals rich in sulfides and selenides. Most of gold is invisible but reaches 345g/t and forms 50nm to 1.5μm Au⁰ particles in a thin 2-3cm interval of organic remnants mixed with wastes I. Vertical and lateral infiltration of AMD waters in peat and oxidative dissolution of wastes within the dispersion train of the Ursk tailings lead to redistribution of elements and their accumulation by combined physical (material's permeability, direction AMD), chemical (complexing, sorption by organic matter and Fe(III) hydroxides) and biochemical (metabolism of sulfate-reducing bacteria) processes. The accumulated elements form secondary sulfates, and Hg and Zn selenides. The results provide insights into accumulation of elements in the early history of coal and black shale deposits and have implications for remediation of polluted areas and for secondary enrichment technologies.

Evaluation of the Single Dilute (0.43 M) Nitric Acid Extraction to Determine Geochemically Reactive Elements in Soil

Recently a dilute nitric acid extraction (0.43 M) was adopted by ISO (ISO-17586:2016) as standard for extraction of geochemically reactive elements in soil and soil like materials. Here we evaluate the performance of this extraction for a wide range of elements by mechanistic geochemical modeling. Model predictions indicate that the extraction recovers the reactive concentration quantitatively (>90%). However, at low ratios of element to reactive surfaces the extraction underestimates reactive Cu, Cr, As, and Mo, that is, elements with a particularly high affinity for organic matter or oxides. The 0.43 M HNO₃ together with more dilute and concentrated acid extractions were evaluated by comparing model-predicted and measured dissolved concentrations in CaCl₂ soil extracts, using the different extractions as alternative model-input. Mean errors of the predictions based on 0.43 M HNO₃ are generally within a factor three, while Mo is underestimated and Co, Ni and Zn in soils with pH > 6 are overestimated, for which possible causes are discussed. Model predictions using 0.43 M HNO₃ are superior to those using 0.1 M HNO₃ or Aqua Regia that under- and overestimate the reactive element contents, respectively. Low concentrations of oxyanions in our data set and structural underestimation of their reactive concentrations warrant further investigation.

Relationship between heavy metals and minerals extracted from soil clay by standard and novel acid extraction procedures

Strong acid digestions are commonly used to determine heavy metal (HM) contents in soils. In order to understand more fully the acid digestion processes, a logical step is to determine the extent of dissolution of mineral phases. The aims of this study were to compare the efficiency of extraction of HM by different acid digestions and to monitor the associated dissolution of the clay fraction. The context of the study was to develop a milder chemical extraction method (microwave-assisted 1 mol L^-1 HNO₃ closed system (NACS)), which recovers more reactive HM and with little dissolution of minerals. The different acid digestion methods dissolved different amounts of minerals from the clay fraction. Both aqua regia (AR) and EPA 3051 dissolved all of the Fe and Al oxides, and the dissolution of kaolin was limited to thinner particles (c dimension), smaller particles in a and b dimensions and grains with lower crystallinity. The lower recovery of HM for AR compared with EPA 3051 was related to the large amount of short-range order phases formed during the AR extraction as these phases have the capacity to re-adsorb HM. The new method (NACS) has the potential to replace other methods of determining bioavailable forms of HM, such as AR and EPA 3051. The contents of Pb, As, Co, Zn, and Cu determined by EPA 3051 and EPA 3052 were quite close.

Temporal variability in trace metal solubility in a paddy soil not reflected in uptake by rice (Oryza sativa L.)

Alternating flooding and drainage conditions have a strong influence on redox chemistry and the solubility of trace metals in paddy soils. However, current knowledge of how the effects of water management on trace metal solubility are linked to trace metal uptake by rice plants over time is still limited. Here, a field-contaminated paddy soil was subjected to two flooding and drainage cycles in a pot experiment with two rice plant cultivars, exhibiting either high or low Cd accumulation characteristics. Flooding led to a strong vertical gradient in the redox potential (Eh). The pH and Mn, Fe, and dissolved organic carbon concentrations increased with decreasing Eh and vice versa. During flooding, trace metal solubility decreased markedly, probably due to sulfide mineral precipitation. Despite its low solubility, the Cd content in rice grains exceeded the food quality standards for both cultivars. Trace metal contents in different rice plant tissues (roots, stem, and leaves) increased at a constant rate during the first flooding and drainage cycle but decreased after reaching a maximum during the second cycle. As such, the high temporal variability in trace metal solubility was not reflected in trace metal uptake by rice plants over time. This might be due to the presence of aerobic conditions and a consequent higher trace metal solubility near the root surface, even during flooding. Trace metal solubility in the rhizosphere should be considered when linking water management to trace metal uptake by rice over time.

A structural equation model of soil metal bioavailability to earthworms: confronting causal theory and observations using a laboratory exposure to field-contaminated soils

Structural equation models (SEM) are increasingly used in ecology as multivariate analysis that can represent theoretical variables and address complex sets of hypotheses. Here we demonstrate the interest of SEM in ecotoxicology, more precisely to test the three-step concept of metal bioavailability to earthworms. The SEM modeled the three-step causal chain between environmental availability, environmental bioavailability and toxicological bioavailability. In the model, each step is an unmeasured (latent) variable reflected by several observed variables. In an exposure experiment designed specifically to test this SEM for Cd, Pb and Zn, Aporrectodea caliginosa was exposed to 31 agricultural field-contaminated soils. Chemical and biological measurements used included CaC12-extractable metal concentrations in soils, free ion concentration in soil solution as predicted by a geochemical model, dissolved metal concentration as predicted by a semi-mechanistic model, internal metal concentrations in total earthworms and in subcellular fractions, and several biomarkers. The observations verified the causal definition of Cd and Pb bioavailability in the SEM, but not for Zn. Several indicators consistently reflected the hypothetical causal definition and could thus be pertinent measurements of Cd and Pb bioavailability to earthworm in field-contaminated soils. SEM highlights that the metals present in the soil solution and easily extractable are not the main source of available metals for earthworms. This study further highlights SEM as a powerful tool that can handle natural ecosystem complexity, thus participating to the paradigm change in ecotoxicology from a bottom-up to a top-down approach.

Kinetics of Heavy Metal Dissociation from Natural Organic Matter: Roles of the Carboxylic and Phenolic Sites

We developed a unifying model for the kinetics of heavy metal dissociation from natural organic matter (NOM) in this study. The kinetics model, integrated with the equilibrium model WHAM 7, specifically considered metal ion reactions with various NOM sites formed by the carboxylic and phenolic sites. The association and dissociation rate coefficients for metal reactions with various NOM sites were constrained by WHAM predicted equilibrium distribution coefficients at specific reaction conditions. We developed the relationship for the dissociation rate coefficients among different binding sites for each metal, which was internally constrained by the metal binding constants. The model had only one fitting parameter, the dissociation rate coefficient for the metal complexes formed with two weak carboxylic sites, and all other parameters were derived from WHAM 7. The kinetic data for metal dissociation from NOM were collected from the literature, and the model was able to reproduce most of relevant data analyzed. The bidentate complexes appeared to be the predominated species controlling metal dissociation under most environmental conditions. The model can help to predict the reactivity and bioavailability of heavy metals under the impact of multiple competing ligands including NOM.

Sorption and pH determine the long-term partitioning of cadmium in natural soils

The bioavailability of metals in soil is a dynamic process. For a proper extrapolation to the field of laboratory studies on fate and effects, it is important to understand the dynamics of metal bioavailability and the way it is influenced by soil properties. The aim of this study was to assess the parallel (concurrent) effect of pH and aging time on the partitioning of cadmium in natural LUFA 2.2 soil. Cadmium nitrate-spiked pH-amended LUFA 2.2 soils were incubated under laboratory conditions for up to 30 weeks. Measured pHpw was lower after 3 weeks and decreased only slightly toward the end of the test. Cadmium concentrations in the pore water increased with time for all soil pH levels, while they decreased with increasing pH. Freundlich kf values ranged between 4.26 and 934 L kg(-1) (n = 0.79 to 1.36) and were highest at the highest pH tested (pH = 6.5). Multiple linear regression analysis, based on a soil ligand modeling approach, resulted in affinity constants of 2.61 for Ca(2+) (log KCa-SL) and 5.05 for H(+) (log KH-SL) for their binding to the active sites on the soil surface. The results showed that pH and aging time are two important factors which together affect cadmium partitioning and mobility in spiked natural soils.

Using isotopic dilution to assess chemical extraction of labile Ni, Cu, Zn, Cd and Pb in soils

Chemical extractants used to measure labile soil metal must ideally select for and solubilise the labile fraction, with minimal solubilisation of non-labile metal. We assessed four extractants (0.43 M HNO3, 0.43 M CH3COOH, 0.05 M Na2H2EDTA and 1 M CaCl2) against these requirements. For soils contaminated by contrasting sources, we compared isotopically exchangeable Ni, Cu, Zn, Cd and Pb (EValue, mg kg(-1)), with the concentrations of metal solubilised by the chemical extractants (MExt, mg kg(-1)). Crucially, we also determined isotopically exchangeable metal in the soil-extractant systems (EExt, mg kg(-1)). Thus 'EExt - EValue' quantifies the concentration of mobilised non-labile metal, while 'EExt - MExt' represents adsorbed labile metal in the presence of the extractant. Extraction with CaCl2 consistently underestimated EValue for Ni, Cu, Zn and Pb, while providing a reasonable estimate of EValue for Cd. In contrast, extraction with HNO3 both consistently mobilised non-labile metal and overestimated the EValue. Extraction with CH3COOH appeared to provide a good estimate of EValue for Cd; however, this was the net outcome of incomplete solubilisation of labile metal, and concurrent mobilisation of non-labile metal by the extractant (MExt<EExt>EValue). The Na2H2EDTA extractant mobilised some non-labile metal in three of the four soils, but consistently solubilised the entire labile fraction for all soil-metal combinations (MExt ≈ EExt). Comparison of EValue, MExt and EExt provides a rigorous means of assessing the underlying action of soil chemical extraction methods and could be used to refine long-standing soil extraction methodologies.