Addressing geographic variability in the comparative toxicity potential of copper and nickel in soils

Bioavailability (BA)-based risk assessment of soil heavy metals in provinces of China through the predictive BA-models

The real biological effect is not generated by the total content of heavy metals (HMs), but rather by bioavailable content. A new bioavailability-based ecological risk assessment (BA-based ERA) framework was developed for deriving bioavailability-based soil quality criteria (BA-based SQC) and accurately assessing the ecological risk of soil HMs at a multi-regional scale in this study. Through the random forest (RF) models and BA-based ERA framework, the 217 BA-based SQC for HMs in 31 Chinese provinces were derived and the BA-based ERA was comprehensively assessed. This study found that bioavailable HMs extraction methods (BHEMs) and total HMs content play the predominant role in affecting HMs (As, Cd, Cr, Cu, Ni, Pb, and Zn) bioavailability by explaining 27.55-56.11% and 9.20-62.09% of the variation, respectively. The RF model had accurate and stable prediction ability for the bioavailability of soil HMs with the mean R2 and RMSE of 0.83 and 0.43 for the test set, respectively. The results of BA-based ERA showed that bioavailability could avoid the overestimation of ecological risks to some extent after reducing the uncertainty of soil differences. This study confirmed the feasibility of using bioavailability for ERA and will utilised to revise the soil environmental standards based on bioavailability for HMs.

Environmental risk thresholds and prediction models of Cd in Chinese agricultural soils

Soil environmental risk threshold of cadmium (Cd) is an important index in formulating soil protection policy. Environmental risk threshold refers to the maximal allowable critical concentration of hazardous substances in the environment. Although there is less study on how to determine soil Cd environmental risk threshold, it is a crucial indicator in formulating soil conservation policies and a key factor in assessing soil environmental quality. The main research content of the study is deducing the environmental risk threshold, aiming to provide scientific basis for the study of environmental quality standards of agricultural land and provide technical support for the protection of Cd pollution of agricultural land. The hazard concentration of 5 % species (HC5, which protects 95 % of species) was determined here using different toxicological data of Cd from 23 test endpoints, interspecific extrapolation using the species sensitivity distribution (SSD) method, and a prediction model was created on the basis of several soil parameters. According to the findings, Cd effective concentration (EC10) (Cd concentration which blocks 10 % of an endpoint's bioactivity) varied from 0.109 to 221 mg·kg-1, and the hormetic response induced by Cd reached 118 % displaying in the dose-response curve of Lolium perenne L.. Toxicology data was rectified by the aging factor considering biogeochemical processes of the newly added pollutants prior to SSD curves fitting. After that, the prediction model was created with the equation of LogHC5 = 0.147 pH + 0.067 OC -1.616. The field test properly validated the prediction model, demonstrating its ability to forecast Cd toxicity levels for various soil conditions. This study offers a scientifically sound methodology for determining the environmental risk limitation for Cd and identifies critical paths for the preservation of environmental species.

Ecotoxicity characterization of chemicals: Global recommendations and implementation in USEtox

Chemicals emitted to the environment affect ecosystem health from local to global scale, and reducing chemical impacts has become an important element of European and global sustainability efforts. The present work advances ecotoxicity characterization of chemicals in life cycle impact assessment by proposing recommendations resulting from international expert workshops and work conducted under the umbrella of the UNEP-SETAC Life Cycle Initiative in the GLAM project (Global guidance on environmental life cycle impact assessment indicators). We include specific recommendations for broadening the assessment scope through proposing to introduce additional environmental compartments beyond freshwater and related ecotoxicity indicators, as well as for adapting the ecotoxicity effect modelling approach to better reflect environmentally relevant exposure levels and including to a larger extent chronic test data. As result, we (1) propose a consistent mathematical framework for calculating freshwater ecotoxicity characterization factors and their underlying fate, exposure and effect parameters; (2) implement the framework into the USEtox scientific consensus model; (3) calculate characterization factors for chemicals reported in an inventory of a life cycle assessment case study on rice production and consumption; and (4) investigate the influence of effect data selection criteria on resulting indicator scores. Our results highlight the need for careful interpretation of life cycle assessment impact scores in light of robustness of underlying species sensitivity distributions. Next steps are to apply the recommended characterization framework in additional case studies, and to adapt it to soil, sediment and the marine environment. Our framework is applicable for evaluating chemicals in life cycle assessment, chemical and environmental footprinting, chemical substitution, risk screening, chemical prioritization, and comparison with environmental sustainability targets.

Influence of metal speciation on soil ecotoxicity impacts in life cycle assessment

It is unknown whether metallic elements remain important contributors to terrestrial ecotoxicity impact scores in life cycle assessment (LCA) when solid- and liquid-phase speciation are considered in environmental fate, exposure and effects. Here, a new speciation-based method for calculating comparative toxicity potentials (CTP) of 23 metallic elements in soils was compared with two other widely used methods which do not consider speciation (i.e., IMPACT, 2002+ and ReCiPe 2008) for nearly 13,000 life cycles of unit processes taken from different sectors. Differences in impact scores between method were driven either by differences in characterization models (ReCiPe 2008) or both by differences in characterization models and substance coverage (IMPACT, 2002+). Strong correlations (r > 0.98) and seemingly constant shifts in impact scores were found for those processes where one or few substances (usually metals) contributed most to total impact and there were large differences in CTPs between methods for these substances. Weaker correlations but often better agreement in impact scores were found for those processes where organic substances were dominant contributors to total impact. Our results suggest that metals are expected to remain important contributors to soil ecotoxicity impacts in LCA when speciation is considered.

Regionalizing eco-toxicity characterization factors for copper soil emissions considering edaphic information for Northern Spain and Portuguese vineyards

The management of vineyards depends on the use of plant protection agents. Regardless of the numerous environmental impacts that these pesticides generate during their production, their dosage as pest control agents in vineyards causes an important toxic effect that must be monitored. Copper-based inorganic pesticides are the most widely used agents to control fungal diseases in humid wine-growing regions. It is, however, significant that the environmental analysis of their use through the Life Cycle Assessment (LCA) methodology does not provide detailed information on the potential toxicity of this type of pesticides. Hence, most studies report average values for copper characterization factors (CFs), excluding local soil characteristics. The objective of the study was the spatial characterization of the ecotoxicity factors of copper soil emissions as a function of the chemical characteristics of vineyard soils located in Portugal and Galicia (NW Spain). A multiple linear regression model was applied to calculate the comparative toxic potential. Subsequently, CFs for copper were calculated based on spatial differentiation considering the variable properties of the soil within each wine appellation. The CFs obtained for the area evaluated ranged from 141 to 5937 PAF·m³·day/kg_{Cu emitted}, for fibric histosols (HSf) and dystic cambisols (CMd), respectively. Moreover, the average values obtained for Galician and Portuguese soils were 1145 and 2274 PAF·m³·day/kg_{Cu emitted}, respectively. The results obtained illustrate the high variability of CF values as a function of the chemical characteristics of each type of soil. For example, Cu soil mobility was linked to organic carbon content and pH. Finally, to validate the representativeness of the calculated CFs, these were applied to the results of 12 literature life cycle inventories of grape production in the area evaluated, revealing that impact scores associated with Cu emissions can considerably vary when spatially-differentiated CFs are implemented.

Overview and recommendations for regionalized life cycle impact assessment

PURPOSE

Regionalized life cycle impact assessment (LCIA) has rapidly developed in the past decade, though its widespread application, robustness, and validity still faces multiple challenges. Under the umbrella of UNEP/SETAC Life Cycle Initiative, a dedicated cross-cutting working group on regionalized LCIA aims to provides an overview of the status of regionalization in LCIA methods. We give guidance and recommendations to harmonize and support regionalization in LCIA for developers of LCIA methods, LCI databases, and LCA software.

METHOD

A survey of current practice among regionalized LCIA method developers was conducted. The survey included questions on chosen method spatial resolution and scale, the spatial resolution of input parameters, choice of native spatial resolution and limitations, operationalization and alignment with life cycle inventory data, methods for spatial aggregation, the assessment of uncertainty from input parameters and model structure, and variability due to spatial aggregation. Recommendations are formulated based on the survey results and extensive discussion by the authors.

RESULTS AND DISCUSSION

Survey results indicate that majority of regionalized LCIA models have global coverage. Native spatial resolutions are generally chosen based on the availability of global input data. Annual modelled or measured elementary flow quantities are mostly used for aggregating characterization factors (CFs) to larger spatial scales, although some use proxies, such as population counts. Aggregated CFs are mostly available at the country level. Although uncertainty due to input parameter, model structure, and spatial aggregation are available for some LCIA methods, they are rarely implemented for LCA studies. So far, there is no agreement if a finer native spatial resolution is the best way to reduce overall uncertainty. When spatially differentiated models CFs are not easily available, archetype models are sometimes developed.

CONCLUSIONS

Regionalized LCIA methods should be provided as a transparent and consistent set of data and metadata using standardized data formats. Regionalized CFs should include both uncertainty and variability. In addition to the native-scale CFs, aggregated CFs should always be provided, and should be calculated as the weighted averages of constituent CFs using annual flow quantities as weights whenever available. This paper is an important step forward for increasing transparency, consistency and robustness in the development and application of regionalized LCIA methods.

Improving the Life Cycle Impact Assessment of Metal Ecotoxicity: Importance of Chromium Speciation, Water Chemistry, and Metal Release

Investigations of metal ecotoxicity in life cycle assessment (LCA) and life cycle impact assessment (LCIA) are becoming important tools for evaluating the environmental impact of a product or process. There is, however, improvement needed for LCIA of metal ecotoxicity in order to make this assessment more relevant and robust. In this work, three issues within the LCIA of metal ecotoxicity are investigated, mainly focusing on topics related to stainless steel manufacturing. The first issue is the importance of considering regional water chemistry when constructing the characterization factor (CF). A model freshwater of relevance for stainless steel manufacturing in a region of Sweden was created with chemistry different from available options. The second issue is related to the lack of consideration on changes in speciation of Cr(VI) in freshwater for a given emission, as Cr(VI) to some extent will be reduced to Cr(III). Two new options are suggested based on relationships between the Cr(VI)–total Cr ratio as a way to improve the relevancy of LCIA for Cr(VI) in freshwater. The last issue is how to treat metal release from slags in LCIA. Metal release from slags was shown to vary significantly between different ways of modelling slag emissions (differences in total metal content, slag leaching tests, estimated emissions to groundwater).

Toward harmonizing ecotoxicity characterization in life cycle impact assessment

Ecosystem quality is an important area of protection in life cycle impact assessment (LCIA). Chemical pollution has adverse impacts on ecosystems on a global scale. To improve methods for assessing ecosystem impacts, the Life Cycle Initiative hosted by the United Nations Environment Programme established a task force to evaluate the state-of-the-science in modeling chemical exposure of organisms and the resulting ecotoxicological effects for use in LCIA. The outcome of the task force work will be global guidance and harmonization by recommending changes to the existing practice of exposure and effect modeling in ecotoxicity characterization. These changes will reflect the current science and ensure the stability of recommended practice. Recommendations must work within the needs of LCIA in terms of 1) operating on information from any inventory reporting chemical emissions with limited spatiotemporal information, 2) applying best estimates rather than conservative assumptions to ensure unbiased comparison with results for other impact categories, and 3) yielding results that are additive across substances and life cycle stages and that will allow a quantitative expression of damage to the exposed ecosystem. We describe the current framework and discuss research questions identified in a roadmap. Primary research questions relate to the approach toward ecotoxicological effect assessment, the need to clarify the method's scope and interpretation of its results, the need to consider additional environmental compartments and impact pathways, and the relevance of effect metrics other than the currently applied geometric mean of toxicity effect data across species. Because they often dominate ecotoxicity results in LCIA, we give metals a special focus, including consideration of their possible essentiality and changes in environmental bioavailability. We conclude with a summary of key questions along with preliminary recommendations to address them as well as open questions that require additional research efforts. Environ Toxicol Chem 2018;37:2955-2971. © 2018 SETAC.

Terrestrial Ecotoxic Impacts Stemming from Emissions of Cd, Cu, Ni, Pb and Zn from Manure: A Spatially Differentiated Assessment in Europe

Metallic elements present in livestock manure as co-contaminants have the potential to cause terrestrial ecotoxic impacts when the manure is used as fertilizer on agricultural soils. The magnitude of this impact at country scale in Europe has, to date, not been quantified. Here, we address this knowledge gap by combining recently developed national emission inventories of Cd, Cu, Ni, Pb and Zn releases from manure with metal- and soil-specific comparative toxicity potentials (CTP) calculated for cropland grid cells at 1 × 1 km resolution for 33 European countries. The CTPs account for speciation in environmental fate, exposure and effects, including reduction in the solid-phase reactivity of a metal when it is associated with organic carbon present in the manure. Given the scarcity of inventory data at sub-national level, it was assumed that each unit area of cropland within a given country has the same probability to receive manure. The resulting CTPs span a range of several orders of magnitude reflecting the influence of soil type and properties on the speciation patterns and resulting CTP values. However, when combined with the use of manure in each European country, the resulting national impact scores were mainly explained by the total mass input of metal released to soil rather than by geographic variability in the CTP values. Simple linear regression is then sufficient to predict terrestrial ecotoxic impacts from input mass. Although some changes in ranking of metals and countries were observed, both mass- and impact-based comparisons between metals agreed that Zn and Cu are dominant contributors to total impacts, and that top contributing countries were those emitting the largest amounts of metals. Our findings show that spatially differentiated impact assessment is important for ranking of countries when differences in national emission inventories between countries are smaller than a factor of two (Ni), a factor of three (Cd, Cu, Zn) or a factor of four (Pb). In other cases, ranking of countries can be based on national emission inventories.

The influence of aging on the comparative terrestrial ecotoxicity potential of copper and zinc in soils

Metal exposure to terrestrial organism is influenced by the reactivity of the solid-phase metal pool. Aging is one of the important factors that control the reactivity of the solid-phase metal pool in soil. In this study, the selected 13 soils were collected from different locations of China, representing different soil types. The reactivity variation of spiked Cu and Zn with aging was assessed in these 13 soils, and their comparative toxicity potentials (CTPs) were also calculated. The median reactive fractions (f_reactive) of Cu and Zn with 95% confidence intervals were 1.6 × 10^-2 (3.5 × 10^-6 to 2.2 × 10^-1) and 0.10 (9.1 × 10^-4 to 0.44) kg_reactive/kg_total, and the median CTPs for Cu and Zn were 2.09 (8.1 × 10^-4 to 2.2 × 10⁴) and 0.85 (8.5 × 10^-4 to 7.2 × 10²) m³/kg day, respectively. The statistical analysis indicated that aging variability in the CTP of Cu and Zn was mainly associated with the variability in soil organic carbon and pH. These results stress the importance of dealing with aging in the calculation of CTPs for terrestrial ecotoxicity of metals.

Regionalized Terrestrial Ecotoxicity Assessment of Copper-Based Fungicides Applied in Viticulture

Life cycle assessment has been recognized as an important decision-making tool to improve the environmental performance of agricultural systems. Still, there are certain modelling issues related to the assessment of their impacts. The first is linked to the assessment of the metal terrestrial ecotoxicity impact, for which metal speciation in soil is disregarded. In fact, emissions of metals in agricultural systems contribute significantly to the ecotoxic impact, as do copper-based fungicides applied in viticulture to combat downy mildew. Another issue is linked to the ways in which the intrinsic geographical variability of agriculture resulting from the variation of management practices, soil properties, and climate is addressed. The aim of this study is to assess the spatial variability of the terrestrial ecotoxicity impact of copper-based fungicides applied in European vineyards, accounting for both geographical variability in terms of agricultural practice and copper speciation in soil. This first entails the development of regionalized characterization factors (CFs) for the copper used in viticulture and then the application of these CFs to a regionalized life-cycle inventory that considers different management practices, soil properties, and climates in different regions, namely Languedoc-Roussillon (France), Minho (Portugal), Tuscany (Italy), and Galicia (Spain). There are two modelling alternatives to determine metal speciation in terrestrial ecotoxicity: (a) empirical regression models; and (b) WHAM 6.0, the geochemical speciation model applied according to the soil properties of the Harmonized World Soil Database (HWSD). Both approaches were used to compute and compare regionalized CFs with each other and with current IMPACT 2002+ CF. The CFs were then aggregated at different spatial resolutions—global, Europe, country, and wine-growing region—to assess the uncertainty related to spatial variability at the different scales and applied in the regionalized case study. The global CF computed for copper terrestrial ecotoxicity is around 3.5 orders of magnitude lower than the one from IMPACT 2002+, demonstrating the impact of including metal speciation. For both methods, an increase in the spatial resolution of the CFs translated into a decrease in the spatial variability of the CFs. With the exception of the aggregated CF for Portugal (Minho) at the country level, all the aggregated CFs derived from empirical regression models are greater than the ones derived from the method based on WHAM 6.0 within a range of 0.2 to 1.2 orders of magnitude. Furthermore, CFs calculated with empirical regression models exhibited a greater spatial variability with respect to the CFs derived from WHAM 6.0. The ranking of the impact scores of the analyzed scenarios was mainly determined by the amount of copper applied in each wine-growing region. However, finer spatial resolutions led to an impact score with lower uncertainty.

Modeling ecotoxicity impacts in vineyard production: Addressing spatial differentiation for copper fungicides

Application of plant protection products (PPP) is a fundamental practice for viticulture. Life Cycle Assessment (LCA) has proved to be a useful tool to assess the environmental performance of agricultural production, where including toxicity-related impacts for PPP use is still associated with methodological limitations, especially for inorganic (i.e. metal-based) pesticides. Downy mildew is one of the most severe diseases for vineyard production. For disease control, copper-based fungicides are the most effective and used PPP in both conventional and organic viticulture. This study aims to improve the toxicity-related characterization of copper-based fungicides (Cu) for LCA studies. Potential freshwater ecotoxicity impacts of 12 active ingredients used to control downy mildew in European vineyards were quantified and compared. Soil ecotoxicity impacts were calculated for specific soil chemistries and textures. To introduce spatial differentiation for Cu in freshwater and soil ecotoxicity characterization, we used 7 European water archetypes and a set of 15,034 non-calcareous vineyard soils for 4 agricultural scenarios. Cu ranked as the most impacting substance for potential freshwater ecotoxicity among the 12 studied active ingredients. With the inclusion of spatial differentiation, Cu toxicity potentials vary 3 orders of magnitude, making variation according to water archetypes potentially relevant. In the case of non-calcareous soils ecotoxicity characterization, the variability of Cu impacts in different receiving environments is about 2 orders of magnitude. Our results show that Cu potential toxicity depends mainly on its capacity to interact with the emission site, and the dynamics of this interaction (speciation). These results represent a better approximation to understand Cu potential toxicity impact profiles, assisting decision makers to better understand copper behavior concerning the receiving environment and therefore how restrictions on the use of copper-based fungicides should be considered in relation to the emission site.

Limitations of experiments performed in artificially made OECD standard soils for predicting cadmium, lead and zinc toxicity towards organisms living in natural soils

Development of comparative toxicity potentials of cationic metals in soils for applications in hazard ranking and toxic impact assessment is currently jeopardized by the availability of experimental effect data. To compensate for this deficiency, data retrieved from experiments carried out in standardized artificial soils, like OECD soils, could potentially be tapped as a source of effect data. It is, however, unknown whether such data are applicable to natural soils where the variability in pore water concentrations of dissolved base cations is large, and where mass transfer limitations of metal uptake can occur. Here, free ion activity models (FIAM) and empirical regression models (ERM, with pH as a predictor) were derived from total metal EC50 values (concentration with effects in 50% of individuals) using speciation for experiments performed in artificial OECD soils measuring ecotoxicological endpoints for terrestrial earthworms, potworms, and springtails. The models were validated by predicting total metal based EC50 values using backward speciation employing an independent set of natural soils with missing information about ionic composition of pore water, as retrieved from a literature review. ERMs performed better than FIAMs. Pearson's r for log₁₀-transformed total metal based EC50s values (ERM) ranged from 0.25 to 0.74, suggesting a general correlation between predicted and measured values. Yet, root-mean-square-error (RMSE) ranged from 0.16 to 0.87 and was either smaller or comparable with the variability of measured EC50 values, suggesting modest performance. This modest performance was mainly due to the omission of pore water concentrations of base cations during model development and their validation, as verified by comparisons with predictions of published terrestrial biotic ligand models. Thus, the usefulness of data from artificial OECD soils for global-scale assessment of terrestrial ecotoxic impacts of Cd, Pb and Zn in soils is limited due to relatively small variability of pore water concentrations of dissolved base cations in OECD soils, preventing their inclusion in development of predictive models. Our findings stress the importance of considering differences in ionic composition of soil pore water when characterizing terrestrial ecotoxicity of cationic metals in natural soils.

Effect of liming on nickel bioavailability and toxicity to oat and soybean grown in field soils containing aged emissions from a nickel refinery

Remediation of soils elevated in trace metals so that the soils may provide ecosystems services is typically achieved through pH adjustment or addition of sorbents. The present study aimed to generate higher-tier in situ toxicity data for elevated nickel (Ni) in soils with and without lime addition and to explore the effect of liming on soil chemistry and bioavailability of Ni to plants. A multiyear study of agronomic yield of field-grown oat and soybean occurred in 3 adjacent fields that had received air emissions from a Ni refinery for 66 yr. The soil Ni concentration in the plots ranged between 1300 mg/kg and 4900 mg/kg, and each field was amended with either 50 Mg/ha, 10 Mg/ha, or 0 Mg/ha (or tonnes/ha) of crushed dolomitic limestone. As expected, liming raised the pH of the soils and subsequently reduced the plant availability of Ni. Toxicity thresholds (effective concentrations causing 50% reduction in growth) for limed soils supported the hypothesis that liming reduces toxicity. Relationships were found between relative yield and soil cation exchange capacity and between relative yield and soil pH, corroborating findings of the European Union Risk Assessments and the Metals in Asia studies, respectively. Higher tier ecotoxicity data such as these are a valuable contribution to risk assessment for Ni in soils. Environ Toxicol Chem 2017;36:1110-1119. © 2016 SETAC.

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.

Assessing comparative terrestrial ecotoxicity of Cd, Co, Cu, Ni, Pb, and Zn: The influence of aging and emission source

Metal exposure to terrestrial organisms is influenced by the reactivity of the solid-phase metal pool. This reactivity is thought to depend on the type of emission source, on aging mechanisms that are active in the soil, and on ambient conditions. Our work shows, that when controlling for soil pH or soil organic carbon, emission source occasionally has an effect on reactivity of Cd, Co, Cu, Ni, Pb and Zn emitted from various anthropogenic sources followed by aging in the soil from a few years to two centuries. The uncertainties in estimating the age prevent definitive conclusions about the influence of aging time on the reactivity of metals from anthropogenic sources in soils. Thus, for calculating comparative toxicity potentials of man-made metal contaminations in soils, we recommend using time-horizon independent accessibility factors derived from source-specific reactive fractions.