Determination and validation of soil thresholds for cadmium based on food quality standard and health risk assessment

Soil environmental carrying capacity and its spatial high-precision accounting framework

Human activity intensity should be controlled within the carrying capacity of soil units, which is crucial for environmental sustainability. However, the existing assessment methods for soil environmental carrying capacity (SECC) rarely consider the relationship between human activity intensity and pollutant emissions, making it difficult to provide effective early warning of human activity intensity. Moreover, there is a lack of spatial high-precision accounting methods for SECC. This study first established a spatial soil environmental capacity (SEC) model based on the pollutant thresholds corresponding to the specific protection target. Next, a spatial net-input flux model was proposed based on soil pollutants' input/output fluxes. Then, the quantitative relationship between human activity intensity and pollutant emissions was established and further incorporated into the SECC model. Finally, the spatial high-precision accounting framework of SECC was proposed. The methodology was used to assess the SECC for the copper production capacity in a typical copper smelting area in China. The results showed that (i) the average SECs for Cu, Cd, Pb, Zn, As and Cr are 427.89, 16.84, 306.41, 376.8, 71.63, and 392.7 kg hm-2, respectively; (ii) heavy metal (HM) concentrations and land-use types jointly influence the spatial distribution pattern of SEC; (iii) atmospheric deposition is the dominant HM input pathway and the high net-input fluxes are mainly located in the southeast of the study area; (iv) with the current human activity intensity for 50 years, the average SECs for Cu, Cd, Pb, Zn, As and Cr are 202.31, 1.71, 20.9, 66.15, 36.73, and 3 kg hm-2, respectively; and (v) to maintain the protection target at the acceptable risk level within 50 years, the SECC for the increased copper production capacity is 1.53 × 106 t. This study provided an effective tool for early warning of human activity intensity.

Development of an effective QSAR-based hazard threshold prediction model for the ecological risk assessment of aromatic hydrocarbon compounds

The insufficient hazard thresholds of specific individual aromatic hydrocarbon compounds (AHCs) with diverse structures limit their ecological risk assessment. Thus, herein, quantitative structure-activity relationship (QSAR) models for estimating the hazard threshold of AHCs were developed based on the hazardous concentration for 5% of species (HC5) determined using the optimal species sensitivity distribution models and on the molecular descriptors calculated via the PADEL software and ORCA software. Results revealed that the optimal QSAR model, which involved eight descriptors, namely, Zagreb, GATS2m, VR3_Dzs, AATSC2s, GATS2c, ATSC2i, ω, and Vm, displayed excellent performance, as reflected by an optimal goodness of fit (R2adj = 0.918), robustness (Q2LOO = 0.869), and external prediction ability (Q2F1 = 0.760, Q2F2 = 0.782, and Q2F3 = 0.774). The hazard thresholds estimated using the optimal QSAR model were approximately close to the published water quality criteria developed by different countries and regions. The quantitative structure-toxicity relationship demonstrated that the molecular descriptors associated with electrophilicity and topological and electrotopological properties were important factors that affected the risks of AHCs. A new and reliable approach to estimate the hazard threshold of ecological risk assessment for various aromatic hydrocarbon pollutants was provided in this study, which can be widely popularised to similar contaminants with diverse structures.

Ecological and health risk assessments of heavy metals and their accumulation in a peanut-soil system

Heavy metals pollution is a notable threat to environment and human health. This study evaluated the potential ecological and health risks of heavy metals (Cu, Cr, Cd, Pb, Zn, Ni, and As) and their accumulation in a peanut-soil system based on 34 soil and peanut kernel paired samples across China. Soil As and Cd posed the greatest pollution risk with 47.1% and 17.6% of soil samples exceeding the risk screen levels, respectively, with 26.5% and 20.6% of the soil sites at relatively strong potential ecological risk level, respectively, and with the geo-accumulation levels at several soil sites in the uncontaminated to moderately contaminated categories. About 35.29% and 2.94% of soil sites were moderately and severely polluted based on Nemerow comprehensive pollution index, respectively, and a total of 32.4% of samples were at moderate ecological hazard level based on comprehensive potential ecological risk index values. The Cd, Cr, Ni, and Cu contents exceeded the standard in 11.76, 8.82, 11.76 and 5.88% of the peanut kernel samples, respectively. Soil metals posed more health risks to children than adults in the order As > Ni > Cr > Cu > Pb > Zn > Cd for non-carcinogenic health risks and Ni > Cr ≫ Cd > As > Pb for carcinogenic health risks. The soil As non-cancer risk index for children was greater than the permitted limits at 14 sites, and soil Ni and Cr posed the greatest carcinogenic risk to adults and children at many soil sites. The metals in peanut did not pose a non-carcinogenic risk according to standard. Peanut kernels had strong enrichment ability for Cd with an average bio-concentration factor (BCF) of 1.62. Soil metals contents and significant soil properties accounted for 35-74% of the variation in the BCF values of metals based on empirical prediction models.

Factors driving antimony accumulation in soil-pakchoi and wheat agroecosystems: Insights and predictive models

The accumulation of antimony (Sb) in plants and its potential effects on human health are of increasing concern. Nevertheless, only a few countries or regions have established soil Sb thresholds for agricultural purposes, and soil properties have not been taken into account. This study investigated the accumulation of Sb in the edible parts of pakchoi and wheat grain by adding exogenous Sb to 21 soils with varying properties. The results revealed a positive correlation between bioavailable Sb (Sbava, extracted by 0.1 M K2HPO4) in soil and Sb in the edible parts of pakchoi (R2 = 0.77, p < 0.05) and wheat grain (R2 = 0.54, p < 0.05). Both machine learning and traditional multiple regression analysis indicated Sbava was the most critical feature and the main soil properties that contributed to Sb uptake by pakchoi and wheat were CaCO3 and clay, respectively. The advisory food limits for Sb in pakchoi and wheat were estimated based on health risk assessment, and used to derive soil thresholds for safe pakchoi and wheat production based on Sbtot and Sbava, respectively. These findings hold potential for predicting Sb uptake by crops with different soil properties and informing safe production management strategies.

Predictive modeling of methylmercury in rice (Oryza sativa L.) and species-sensitivity-distribution-based derivation of the threshold of soil mercury in karst mountain areas

The bioavailable mercury (Hg) in the soil is highly active and can affect the formulation of methyl-Hg (MeHg) in soil and its accumulation in rice. Herein, we predicted the concentration of MeHg in rice using bioavailable Hg extracted from soils; additionally, we determined the threshold value of soil Hg in karst mountain areas based on species sensitivity distribution. The bioavailable Hg was extracted using calcium chloride, hydrochloric acid (HCl), diethylenetriaminepentaacetic acid mixture, ammonium acetate, and thioglycolic acid. Results showed that HCl is the best extractant, and the prediction model demonstrated good predictability of the MeHg concentration in rice based on the HCl-extractable Hg, pH, and soil organic matter (SOM) data. Compared with the actual MeHg concentration in rice, approximately 99% of the predicted values (n = 103) were within the 95% prediction range, indicating the good performance of the rice MeHg prediction model based on soil pH, SOM, and bioavailable Hg in karst mountain areas. Based on this MeHg prediction model, the safety threshold of soil Hg was calculated to be 0.0936 mg/kg, which is much lower than the soil pollution risk screening value of agricultural land (0.5 mg/kg), suggesting that a stricter standard should be applied regarding soil Hg in karst mountain areas. This study presents the threshold of soil Hg pollution for rice safety in karst mountain areas, and future studies should target this threshold range.

Comprehensive assessment of health and ecological risk of cadmium in agricultural soils across China: A tiered framework

Soil cadmium (Cd) contamination has been increasingly serious in agricultural land across China, posing unexpected risks to human health concerning crop safety and terrestrial ecosystems. This study collected Cd concentration data from 3388 soil sites in agricultural regions. To assess the Cd risk to crop safety, a comprehensive sampling investigation was performed to develop reliable Soil Plant Transfer (SPT) model. Eco-toxicity tests with representative soils and organism was conducted to construct the Species Sensitivity Distribution (SSD) for ecological risk assessment. Then, a tiered framework was applied based on Accumulation index, deterministic method (Hazard quotient), and probabilistic assessment (Monte Carlo and Joint Probability Curve). The results revealed the widespread Cd enrichment in agricultural soils, mainly concentrated in Central, Southern, and Southwest China. Risk assessments demonstrated the greater risks related to crop safety, while the ecological risks posed by soil Cd were manageable. Notably, agricultural soils in southern regions of China exhibited more severe risks to both crop safety and soil ecosystem, compared to other agricultural regions. Furthermore, tiered methodology proposed here, can be adapted to other trace elements with potential risks to crop safety and terrestrial ecosystem.

Thallium content in vegetables and derivation of threshold for safe food production in soil: A meta-analysis

Soil thallium (Tl) pollution is a serious environmental problem, and vegetables are the primary pathway for human exposure to Tl. Therefore, it is important to investigate the characteristics of soil Tl uptake by vegetables. In this study, the meta-analysis approach was first applied to explore the relationship between Tl content in vegetables and soil environment, as well as key factors influencing soil physical-chemical properties, and to derive soil thresholds for Tl. The results indicated that various types of vegetables have different capabilities for Tl accumulation. Vegetables from contaminated areas showed high Tl accumulation, and the geomean Tl content in different types of vegetables was in the following order: leafy > root-stalk > solanaceous vegetables. Taro and kale had significantly higher capability to accumulate soil Tl among the 35 species studied, with Tl bioconcentration factor values of 0.060 and 0.133, respectively. Pearson correlation analysis and meta-analysis revealed that the Tl content in vegetables was significantly correlated with soil pH and Tl content in soil. The linear predictive model for Tl accumulation in vegetables based on soil Tl content described the data well, and the fitting coefficient R2 increased with soil pH value. According to potential dietary toxicity, the derived soil Tl thresholds for all, leafy and root-stalk vegetables increased with an increase in soil pH, and were in the range of 1.46-6.72, 1.74-5.26 and 0.92-6.06 mg/kg, respectively. The soil Tl thresholds for kale, lettuce and carrot were in the range of 0.24-4.89, 2.94-3.32 and 3.77-14.43 mg/kg, respectively. Ingestion of kale, beet, sweet potato, potato, taro, pepper, turnip, Chinese cabbage, eggplant and carrot poses potential health risks. The study provides scientific guidance for vegetable production in Tl-contaminated areas and can help with the selection of vegetable species suitable for avoiding the absorption of Tl from contaminated soil.

Evaluation of non-cancer risk owing to groundwater fluoride and iron in a semi-arid region near the Indo-Bangladesh international frontier

Groundwater quality in Hili, a semi-arid border region at Indo-Bangladesh border, was investigated in the post-monsoon season of 2021, succeeded by assessment of probabilistic health risk arising from fluoride (F-) and iron (Fe) intake, with the hypothesis that groundwater quality of the region was not satisfactory for human consumption and health, considering earlier reports on high groundwater F- and Fe in few of the neighboring districts. All water samples were found to be potable in terms of Ca2+, Mg2+, Cl-, SO42- and NO3-, , but F- and Fe exceeded prescribed safe limits for drinking water in about 48% and 7% samples. Almost all water samples were found to be good for irrigation in terms of sodium adsorption ratio (SAR), soluble sodium percentage (SSP), Kelly's index (KI), %Na and magnesium ratio (MR). The principal component analysis (PCA) identified three major factors influencing groundwater quality, explaining about 71.8% of total variance and indicated that groundwater quality was primarily influenced by geochemical factors. Carbonate and silicate weathering were mainly responsible for dissolution of minerals in groundwater. Non-carcinogenic risk due to cumulative impact of F-and Fe intake was in the order of THIChildren > THIInfant > THIAdult. As per Monte Carlo simulation run with 5000 trials to ascertain the order of probabilistic health risk, the most dominant governing factors behind non-carcinogenic risk caused by F-and Fe intake were their concentration (Ci) followed by ingestion rate (IR), and exposure duration (ED).

Contamination Status and Ecological Security Thresholds of Fluoride in Farmland around a Phosphorus Chemical Plant in a Karst Area of Southwestern China

The phosphorus chemical plant (PCP) production process leads to the substantial discharge of fluoride into the surrounding environment. However, there is limited research data regarding the pollution levels and ecological safety thresholds of farmland fluoride around PCPs in karst areas. This study evaluates the local contamination status and health risks by determining the fluoride content in farmland and vegetables in vicinity of a PCP in a karst area of Southwestern China. Farmland ecological security threshold was derived based on the species sensitivity distribution (SSD) model. Results showed that the fluoride contents in the soil ranged from 529.36 to 1496.02 mg kg^-1, with the median value of 823.93 mg kg^-1, which was higher than the national background value in China (478 mg kg^-1) by 172.37%. Fluoride contents in vegetables ranged from 0.01 to 25.21 mg kg^-1, with the median value of 1.29 mg kg^-1, which was higher than the limits of contaminants in food (1 mg kg^-1) by 129.00%, and 53.85% of vegetable samples were contaminated. Leafy vegetables showed a higher tendency of fluoride enrichment compared to non-leafy vegetables. Despite there being no significant health risk to the residents via the intake of vegetables overall, there may be potential health risks from the intake of sweet potato leaves. Based on the SSD model, the ecological safety thresholds of fluoride in local farmland was classified into the suitable for planting category (≤174.13 mg kg^-1), safe plant utilization category (174.13-4005.42 mg kg^-1), and strict control of planting category (≥4005.42 mg kg^-1). These findings will provide valuable insights to facilitate the safe cultivation of local agricultural products and optimize the utilization of soil resources.

Remediation of Cd contaminated paddy fields by intercropping of the high- and low- Cd-accumulating rice cultivars

Intercropping cadmium (Cd) hyperaccumulators with crops have been widely applied in the remediation of contaminated farmland soils. However, most studies were done on drylands since the majority of the hyperaccumulators are susceptible to the aquatic environment, making the remediation of Cd-contaminated paddy fields particularly difficult. Our study attempts to address the issue by intercropping the high-Cd-accumulating (henceforth, "high-Cd") rice cultivars with the low-Cd-accumulating (henceforth, "low-Cd") ones, and to study the Cd removal, uptake and translocation during the remediation process. The results indicated that intercropping mode with 20-cm row spacing (intercropping-20 treatment) performed better than the that with 30-cm row spacing (intercropping-30 treatment), while intercropping had stronger impact on late rice compared to early rice. In general, the physiological condition of rice was stable under the intercropping-20 treatment, suggesting the growth of rice was not impeded. For late rice, as the intercropping-20 treatment can significantly reduce soil pH and increase the diethylenetriaminepentaacetic acid extracted Cd (DTPA-extracted Cd) from the rhizosphere soil, Cd accumulated more in the tissues of the high-Cd rice cultivars (H2), and its dry biomass increased. As a result, a drastic improvement in the total Cd removal rate by 38.55 % was noticed. Therefore, the reduction of total Cd concentration in 0-20 cm profile caused by removal, thus it could provide safer soil environment for the growth of low Cd-rice cultivars (L2), leading to a significant drop in the root Cd concentration and safer production of L2. Interestingly, intercropping had no effect on the yield per plant of low-Cd rice cultivars. For early rice, intercropping-20 treatment exerted trivial effects to all aspects. The intercropping-30 treatment has poor representativeness of all indicators because of the large intercropping distance. Our results demonstrate that intercropping of the high-Cd and the low-Cd rice cultivars is a potential mode for Cd remediation in paddy fields.

Accurate derivation and modelling of criteria of soil extractable and total cadmium for safe wheat production

It is significant to establish an accurate model to predict cadmium (Cd) criteria for safe wheat production. More importantly, for better evaluation of the risk of Cd pollution in high natural background areas, the soil extractable Cd criteria are needed. In the present study, the soil total Cd criteria were derived using the method of cultivars sensitivity distribution integrated with soil aging and bioavailability as affected by soil properties. Firstly, the dataset that meet the requirements was established. Dataset from thirty-five wheat cultivars planted in different soils published in literature of five bibliographic databases were screened using designated search strings. Then, the empirical soil-plant transfer model was used to normalize the bioaccumulation data. Afterwards, the soil Cd concentration for protecting 95 % (HC₅) of the species was calculated from species sensitivity distribution curves, and the derived soil criteria were obtained from HC₅ prediction models that based on pH. The process of derivation for soil EDTA-extractable Cd criteria was the same way as the soil total Cd criteria. Soil total Cd criteria ranged from 0.25 to 0.60 mg/kg and soil EDTA-extractable Cd criteria ranged from 0.12 to 0.30 mg/kg. Both the criteria of soil total Cd and soil EDTA-extractable Cd were further validated to be reliable using data from field experiments. The results suggested that the criteria of soil total Cd and soil EDTA-extractable Cd in the study can ensure the safety of Cd in wheat grains and thereby enable local agricultural practitioners to develop appropriate management for croplands.

Derivation and validation of soil cadmium thresholds for the safe farmland production of vegetables in high geological background area

Excessive dietary intake of cadmium (Cd) poses toxicity risks to human health, and it is therefore essential to establish accurate and regionally appropriate soil Cd thresholds that ensure the safety of agricultural products grown in different areas. This study investigated the differences in the Cd accumulation in 32 vegetable varieties and found that the Cd content ranged from 0.01 to 0.24 mg·kg^-1, and decreased in the order of stem and bulb vegetables > leafy vegetables > solanaceous crops > bean cultivars. A correlation analysis and structural equation model showed that pH, soil organic matter, and the cation exchange capacity had significant effects on Cd accumulation in the vegetables and explained 72.1 % of the variance. In addition, species sensitivity distribution (SSD) curves showed that stem and bulb vegetables were more sensitive to Cd than other types of vegetables. Using the Burr Type III function for curve fitting, we derived Cd thresholds of 6.66, 4.15, and 1.57 mg·kg^-1 for vegetable soils. These thresholds will ensure that 20 %, 50 %, and 95 % of these vegetable varieties were risk-free, respectively. The predicted threshold of soil Cd was more than twice that of China's current National Soil Quality Standard (GB 15618-2018) for Cd values. Therefore, soil scenarios and cultivars should be considered comprehensively when determining farmland soil thresholds. The present results provide a new model for setting soil Cd criteria in high geological background areas.

Assessment of soil heavy metal pollution in provinces of China based on different soil types: From normalization to soil quality criteria and ecological risk assessment

Ecological risks can vary dramatically depending on abiotic factors, such as soil properties and the background values of elements. This study developed a framework for an integrated risk assessment system to derive soil quality criteria (SQC) for heavy metals (HMs) applicable to different soil types and to assess ecological risks at a multi-regional scale. Through the construction of normalization and species sensitivity distribution models, 248 SQC values for Cd, Pb, Zn, As, Cu, Cr, Sb, and Ni in 31 Chinese provinces were derived. These SQC considered the soil types and background values of the elements and effectively reduced the uncertainty caused by spatial heterogeneity. Using the derived SQC values, the qualitative and quantitative ecological risks of HMs in the terrestrial environment of China were comprehensively assessed using a three-level ecological risk assessment (ERA) approach. Compared to traditional ERA methods, the new methodology reached a more quantitative conclusion. The mean overall probabilities of ecological risk in China were 2.42 % (Cd), 2.82 % (Pb), 12.17 % (Zn), 14.89 % (As), 10.42 % (Cu), 32.20 %(Cr), and 8.88 % (Ni). The new framework could be useful for the ERA of various soil types.

Establishment and Optimization of Soil Cd Risk Threshold in Typical Karst Area with Potato Production, China

The threshold is key to risk assessment of soil cadmium (Cd) pollution. However, there is limited research on the soil Cd risk threshold of potatoes. Soil and potato samples (n = 256) were used to establish and optimize the Cd risk threshold by using relative cumulative frequency, bioconcentration factor, and regression model. The results showed that suggested risk screening values (SRSVs) for soil Cd were divided into 2.465 (pH ≤ 5.5), 2.564 (5.5 < pH ≤ 6.5), 2.778 (6.5 < pH ≤ 7.5), and 4.348 mg kg^-1 (pH > 7.5). SRSVs were applied to classify soil Cd risk assessment by collecting soil samples (n = 100). Low-risk areas only comprised 0.98% of the total area using risk screening values (RSVs) (GB15618-2018), and risk areas comprised as much as 99.02%. Low-risk area and risk area comprised 97.75% and 2.25% of the total area based on SRSVs. SRSVs are appropriate for potato production in typical karst areas.

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.

Selection of rice and maize varieties with low cadmium accumulation and derivation of soil environmental thresholds in karst

Cadmium (Cd) is considered the primary dietary toxic element. Previous studies have demonstrated significant differences in heavy metal accumulation among crop species. However, this information in karst areas with low heavy metal activity is missing. In this study, the uptake and accumulation characteristics of cadmium in soil-crop samples of group 504 in the core karst region of East Asia were analyzed. Cadmium low-accumulating maize and rice were screened using cluster and Pareto analytic methods. In addition, a new method, the species-sensitive distribution model (SSD), was proposed, which could be used to estimate the environmental threshold for cadmium in regional cropland. The results showed that both maize and rice soils in the research area were contaminated with varying degrees of cadmium. The total concentrations of cadmium ω(T-Cd) in maize and rice fields are 0.18-1.32 and 0.20-4.42 mg kg^-1, respectively. The ω(T-Cd) of heavy metals in maize kernels and rice grains is 0.002-0.429 and 0.003-0.393 mg kg^-1, respectively. The bioaccumulation factor (BCF) of cadmium in maize ranged from 0.0079 to 0.9701, with a coefficient of variation of 1.71; the BCF of cadmium in rice ranged from 0.0074 to 0.1345, with a coefficient of variation of 0.99. According to cluster and Pareto analyses, the maize crop varieties with low cadmium accumulation suitable for local cultivation were screened as JHY809, JDY808, AD778, SN3H and SY13, and the rice varieties were DMY6188, GY725, NY6368, SY451 and DX4103. In addition, the environmental cadmium threshold ranges of 0.30-10.05 mg kg^-1 and 0.89-24.39 mg kg^-1 for maize and rice soils, respectively, were deduced in this study. This threshold will ensure that 5-95% of maize and rice will not be contaminated with cadmium in the soil.

Ecological risk assessment for perfluorohexanesulfonic acid (PFHxS) in soil using species sensitivity distribution (SSD) approach

Perfluorohexanesulfonic acid (PFHxS) is one of the persistent organic pollutants that has been recommended to be listed in Annex A of the Stockholm Convention. It has gained increasing attention in recent years due to its toxic effects. The guideline values of PFHxS are commonly associated with PFOS in various countries and regulatory agencies. In this study, multispecies bioassays were conducted to determine the ecological toxic effects of PFHxS, including plants, soil invertebrates, and soil microorganisms, which indicated the EC10/NOEC values ranged from 2.9 to 250 mg/kg. Where possible, logistic models were used to calculate the EC30 values for various endpoints. The species sensitivity distributions were employed to estimate the ecological investigation levels for PFHxS contamination in soils using toxicity results from literature and this study. The calculation using EC10/NOEC values from both literature and this study indicated a most conservative HC5 as 1.0 mg/kg (hazardous concentration for 5 % of the species being impacted). However, utilisation of EC30 values derived from this study resulted in a much higher HC5 for PFHxS in contaminated soils (13.0 mg/kg) which is at the higher end of the existing guideline values for PFOS for protecting ecological systems. The results obtained in this study can be useful in risk assessment processes to minimize any uncertainty using combined values with PFOS.

Ecological risk assessment of toxic metal(loid)s for land application of sewage sludge in China

Sewage sludge, including those after biological or thermochemical treatments, has the potential to be used as fertilizers for recycle of resources. However, its potential ecological risk is also of great concern to policy making. This study employed comprehensive ecological risk assessment (ERA) methods to evaluate the risk caused by the toxic metal(loid)s in sewage sludge throughout China. The conventional geo-accumulation index and potential ecological risk index revealed that cadmium (Cd) and mercury (Hg) were of significant concern in treating sewage sludge before land application, but chromium (Cr) and zinc (Zn) were preferred by potential affected proportion (PAF) and overall risk probability (ORP) of species sensitivity distribution (SSD). Because SSD considered both the community and the ecotoxicity of toxic metal(loid)s, it was more advantageous and promising in assessing ecological risks caused by land application of sewage sludge. Based on the predicted no-effect concentration (PNEC) of toxic metal(loid) calculated by hazardous concentration that cause death of 50% of species (HC50) by SSD, the maximum allowable disposal amount (MADA) of sewage sludge in the whole China indicated that chromium (Cr) should be totally eliminated because of its high risks in the present background soil. After excluding Cr, the MADA of sewage sludge in China was 3.24 × 10⁶ t and 6.47 × 10⁷ t under land application scenarios with high and low ecological risks, respectively. Additionally, the MADA could be increased by mixing sewage sludge with deeper soil in wider areas. This study emphasized that local laws and regulations on land application of sewage sludge and the subsequent ERA system need to be addressed in the future.

Ecological risk thresholds for Zn in Chinese soils

The environmental risk threshold of a pollutant is a yardstick to measure soil environmental quality. The derivation of ecological risk thresholds of the heavy metal zinc (Zn) in soil environments based on up-to-date ecological risk assessments plays an important role in soil protection policy. According to regional soil classification, different representative soils with various degrees of acidity and alkalinity were selected, and a data set comprising ecotoxicities of Zn to 21 different test endpoints (plants, soil fauna, microorganisms, etc.) found in representative farmland soils of China was compiled based on new and published data to determine toxicological limits of Zn effects on endpoints. These limits were derived from fitted dose-response model parameters and indicated by EC₁₀ values (the effective concentrations of Zn that inhibit 10% of endpoint bioactivity and also represents the toxicity threshold of Zn in this study) ranging from 36 mg·kg^-1 to 682 mg·kg^-1. The hormesis effect appeared in the dose-response curve of Zn, for example, the relative Chinese cabbage growth reached more than 120% at most. Zn concentrations added in toxicity tests were also corrected for aging and leaching effects in order to more accurately reflect field conditions. The hazardous concentrations for 5% of the species affected (HC₅) were derived by the species sensitivity distribution (SSD) approach for four major types of Chinese soils: acidic (38 mg·kg^-1), neutral (106 mg·kg^-1), alkaline (217 mg·kg^-1), and alkaline calcareous soils (155 mg·kg^-1). Prediction models of ecological risk thresholds for Zn based on soil properties were generated, such as logHC₅ = 0.564 + 0.218pH + 0.097OC (R² = 0.790,p < 0.001). The predicted models based on lab test data were verified in the field, and the measured field data fell within two-fold of the prediction intervals. This work provides a scientific framework for developing soil-specific guidance on Zn toxicity thresholds.

Derivation of Soil Criteria of Cadmium for Safe Rice Production Applying Soil–Plant Transfer Model and Species Sensitivity Distribution

Widespread soil contamination is hazardous to agricultural products, posing harmful effects on human health through the food chain. In China, Cadmium (Cd) is the primary contaminant in soils and easily accumulates in rice, the main food for the Chinese population. Therefore, it is essential to derive soil criteria to safeguard rice products by assessing Cd intake risk through the soil–grain–human pathway. Based on a 2-year field investigation, a total of 328 soil–rice grain paired samples were collected in China, covering a wide variation in soil Cd concentrations and physicochemical properties. Two probabilistic methods used to derive soil criteria are soil–plant transfer models (SPT), with predictive intervals, and species sensitivity distribution (SSD), composed of soil type-specific bioconcentration factor (BCF, Cd concentration ratio in rice grain to soil). The soil criteria were back-calculated from the Chinese food quality standard. The results suggested that field data with a proper Cd concentration gradient could increase the model accuracy in the soil–plant transfer system. The derived soil criteria based on soil pH were 0.06–0.11, 0.33–0.59, and 1.51–2.82 mg kg−1 for protecting 95%, 50% and 5% of the rice safety, respectively. The soil criteria with soil pH further validated the soil as being safe for rice grains.

Bibliometric Analysis of the Influencing Factors, Derivation, and Application of Heavy Metal Thresholds in Soil

The study of threshold levels of heavy metals in soil is essential for the assessment and management of soil environmental quality. This study reviewed the influencing factors, the derivation, and application aspects of heavy metals' threshold values comprehensively by a combination of bibliometric analysis and scientific knowledge mapping. A total of 1106 related studies were comprehensively extracted from the Web of Science database during the period from 2001 to 2020. The results showed that the publication output has been growing strongly. An analysis on the subject, journal, country, and institution was carried out to demonstrate the development and evolution of this research branch during the two decades. According to high-frequency keywords analysis, external factors (e.g., soil physicochemical properties) and internal factors (e.g., crop genotype) can affect heavy metal threshold values in the soil-crop system. The current methods mainly include the Point model (e.g., evaluation factor method), the Probability model (e.g., species sensitivity distribution method), and the Empirical model (e.g., ecological environment effect method). A threshold study can be applicable to the risk assessment for soil heavy metal contamination in order to determinate the soil pollution degree and its spatial and temporal distribution characteristics. Moreover, challenges and prospects of the study of heavy metal threshold values are proposed, indicating that research should focus on the relationships between human health risks and the established threshold values of heavy metals in the soil, long-term field trials and bioavailability of heavy metals for the derivation of the thresholds, and the establishment of more scientific and rational soil environmental benchmarks.

Antimony, beryllium, cobalt, and vanadium in urban park soils in Beijing: Machine learning-based source identification and health risk-based soil environmental criteria

The pollution situation of antimony (Sb), beryllium (Be), cobalt (Co), and vanadium (V) is poorly understood, although they are widely used in daily life and production processes. Moreover, threshold levels ("soil environmental criteria", SEC) for these pollutants are lacking in China, which impedes effective soil quality management. This study explored pollution characteristics for park soils in urban area of Beijing, China at first. Then multivariate statistical analysis and machine learning model were used to identify the main sources of pollutants. Additionally, probabilistic health risk and SEC were studied to assess the risks of pollutants and manage soil pollutants. The results revealed that the overall pollution levels of Be, Co, Sb, and V were low, but Be and Sb were enriched to varying degrees. Source apportionment showed that Sb (85.5%) was mainly derived from fuel combustion and industrial legacy, Co (66.7%) and V (82.5%) from natural processes, and Be from the natural background (39.3%) and anthropogenic sources (53.8%). Risk assessment indicated that the pollutants' carcinogenic and noncarcinogenic risks were negligible. Exposure frequency and soil ingestion rate were the most important parameters affecting health risks. The SEC of Be, Co, Sb, and V were 31, 39.7, 41.3, and 348 mg/kg, respectively, all of which are higher than the corresponding soil quality standards in China, indicating that current soil quality standards may be too conservative for urban park land. This study provides a reference for the management of soil pollutants in Beijing's urban parks, and the formulation of soil environmental quality standards.

Theobroma cacao L. agricultural soils with natural low and high cadmium (Cd) in Santander (Colombia), contain a persistent shared bacterial composition shaped by multiple soil variables and bacterial isolates highly resistant to Cd concentrations

•

This study report statistically significant differences in beta-diversity correlated with soil properties, including Cd concentrations.

•

Culture-dependent techniques allowed the isolation of bacterial strains tolerating high Cd concentrations up to 120 mgL⁻¹ for potencial Cd biosorption or intracellular sequestration.

•

The combination of different approaches provides a baseline about the bacterial composition and Cd tolerant strains found in these soils influence for lower Cd accumulation in cocoa crops.

Heavy metals can be found in soil as natural components or as product of contaminations events; plants growing in soils are prone to bioaccumulate heavy metals on their biomass. Theobroma cacao L. can bioaccumulate cadmium (Cd) in the seed and could be in derived food products, it considered a human health risk; therefore, removal of Cd is desirable but not vet technically and economically feasible; only to avoid Cd in cocoa is by selecting lands plots exhibiting lower Cd concentrations in soils, imposing a serious limitation to farmers and regulators. The study of bacterial communities and isolation bacteria with tolerance and mechanisms to counteract the translocation of Cd to the parts of cocoa plant exhibits high relevance in Colombia economy and especially to companies producing chocolate and derivatives. Here, we explore bacterial communities associated with soils having relatively high natural Cd concentrations in a large agricultural cocoa plot located in the Santander region. We characterized the bacterial communities’ compositions by amplicon 16S rRNA sequencing from metagenomics soil DNA and by culturing-based enumeration and isolation approaches.

Culture-dependent techniques allowed the isolation of bacteria tolerant to Cd concentration, complement the information for Colombia, and expand the number of strains characterized with adaptive capacity against Cd with tolerance in a concentration of 120 mg/L, which represents the first capacity for Exiguobacterium sp., Ralstonia sp., Serratia sp., Dermacoccus sp., Klebsiella sp., Lactococcus sp. and Staphylococcus sp. In addition to confirming that there is a greater diversity of Cd-tolerant bacteria present in soils of farms cultivated with cocoa in Colombia. As for the results of new generation sequencing, they revealed that, the alpha-diversity in bacterial composition, according to the ANOVA, there are statistically significant differences of the bacterial communities present in the samples. Regarding Pearson correlation analysis, it was found the Shannon Simpson indices, have a positive correlation against OM, C, pH, Mn, C.E.C.I., Ca, P and negatively correlated with S; respect to bacterial community structure, a principal component analysis, which revealed that independent of the concentration of Cd present in soil samples, separates them according to pH value. Phyla to high abundance relative in all samples were Proteobacteria, Acidobacteriota, Actinobacteriota, Verrucomicrobiota, Myxococcota, Chloroflexi, Plactomycetota, Bacteroidota, Gemmatimonadota, Nitrospirota, Firmicutes and NB1_J; the bacteria genera with higher relative abundance (>0.5%) Nitrospira, candidatus Udaeobacter, Haliangium, Cupriavidus, MND1, Bacillus, Kitasatospora, Niveibacterium, Acidothermus, Burkholderia, Acidibacter, Terrimonas, Gaiella, candidatus Solibacter, Kitasatospora, Sphingomonas, Streptomyces, this genus with a relationship with the Cd tolerance process. After it, redundancy analysis was performed between the variation of the bacterial communities identified by dependent and independent techniques and edaphic soil variables, where their positive correlation was found against K, OM, C, Ca, pH (p<0.01) and P, C.E.C.I (p<0.05). For soil samples, the bacterial genera that make up the core community were identified, which are present in all samples as Nitrospira sp., Cupriavidus sp., Burkholderia sp., Haliangium sp., candidatus Udaeobacter, MND1, Kitasatospora, Acidothermus, Acidibacter, Streptomyces, Gaiella, candidatus Solibacter and Terramonas; the genera identified has a different and fundamental role in ecosystem functioning. The combination of different approaches offers new clues regarding the assessment of bacterial communities in soils cultivated with cocoa in soils with elevated Cd content in Colombia, and the ecological role and interplay of soil components and bacterial communities that contribute to modulate the effect of bioaccumulation in products.

Derivation and validation of thresholds of cadmium, chromium, lead, mercury and arsenic for safe rice production in paddy soil

Cadmium (Cd), chromium (Cr), lead (Pb), mercury (Hg) and arsenic (As) are potent toxicants to human health via dietary intake. It is imperative to establish accurate soil thresholds based on soil-plant transfer models and food safety standards for safe agricultural production. This study takes rice genotypes and soil properties into account to derive soil thresholds for five heavy metal(loid)s using the bioconcentration factors (BCF) and species sensitivity distribution (SSD) based on the food safety standard. The BCF generated from two paddy soils was calculated to investigate the sensitivity of heavy metal accumulation in nine rice cultivars in a greenhouse pot experiment. Then, empirical soil-plant transfer models were developed from a middle-sensitivity rice cultivar (Denong 2000, one selected from nine rice) grown in nineteen paddy soils with various soil properties under a proper exogenously metal(loid)s concentration gradient. After normalization, hazardous concentrations from the fifth percentile (HC5) were calculated from the SSD curves, and the derived soil thresholds were obtained from HC5 prediction models that based on the combination of pH and organic carbon (OC) or cation exchange capacity (CEC). The soil Cd threshold derived based on pH and organic carbon (pH < 7.5, OC ≥ 20 g kg^-1) was 1.3-fold of those only considering pH, whereas the Pb threshold (pH > 6, CEC ≥ 20 cmol_c kg^-1) was 3.1 times lower than the current threshold. The derived thresholds for five elements were validated to be reliable through literature data and field experiments. The results suggested that deriving soil heavy metal(loid)s threshold using SSD method and local food safety standards is feasible and also applicable to other crops as well as other regions with potential health risks of toxic elements contamination in agricultural production.

Integrating bioavailability and aging in the criteria derivation of cadmium for the safe production of rice in paddy soils

It is significant to derive an appropriate cadmium (Cd) threshold for the rice to ensure that the Cd concentration of rice grains meets the food safety standards. In the present study, soil thresholds for Cd were derived using the method of species sensitivity distribution based on aging time, cultivars, and soil properties. Dataset from thirty-nine rice cultivars planted in different soils published in literature of five bibliographic databases were screened using designated search strings to explore their differences in Cd accumulation capacity in paddy soil. The empirical soil-plant transfer model was used to normalize the bioaccumulation data. Later, the soil Cd concentration for protecting 95% (HC5) of the cultivars was calculated using species sensitivity distribution curves fitted by the Burr III function. The soil Cd criteria derived from the added approach for risk assessment were proposed as continuous criteria based on the combination of soil pH and organic carbon in soil. Soil Cd criteria ranged from 0.34 to 1.18 mg/kg. The prediction model for HC5 was applied to field experimental data to validate its validity and applicability. The predicted Cd thresholds were less than the field experiment Cd thresholds. In conclusion, this study provided valuable and scientific bases for setting soil Cd criteria for paddy soils.

A phytoremediation coupled with agro-production mode suppresses Fusarium wilt disease and alleviates cadmium phytotoxicity of cucumber (Cucumis sativus L.) in continuous cropping greenhouse soil

Cadmium (Cd) contamination and continuous cropping obstacle often coexist in greenhouse soil and seriously restrict cucumber production. In this study, hyperaccumulator Sedum alfredii Hance was intercropped with spring cucumber (Cucumis sativus L.), then rotated with low accumulator water spinach and autumn cucumber under rational water regime, composited amendment was applied to soil before transplanting autumn cucumber. The results showed that, compared with conventional crop rotation system (Chinese cabbage and cucumber rotation), superposition management practice suppressed Fusarium wilt disease by 28.4 and 57.4% and increased yield by 35.2 and 383% for spring and autumn cucumbers, respectively. Meanwhile, photosynthetic characteristics, antioxidant system and fruit quality were significantly improved. Furthermore, this mode modified soil microbial community structure, enhanced soil enzyme activities, and simultaneously reduced soil total and phytoavailable Cd by 30.3 and 47.7%, respectively. These results demonstrated a feasible technical system to achieve phytoremediation coupled with argo-production in Cd contaminated greenhouse soil with continuous cropping obstacle and provided useful information for further revelation of interaction mechanisms between multicropping and comprehensive biofortification measurements.

A new perspective on the toxicity of arsenic-contaminated soil: Tandem mass tag proteomics and metabolomics in earthworms

The toxicity of low-level arsenic (As)-contaminated soil is not well understood. An integrated proteomic and metabolomic approach combined with morphological examination was used to investigate the potential biological toxicity of As-contaminated soil based on an exposure experiment with the earthworm Eisenia fetida. The results showed that the earthworm hindgut accumulated high As concentrations resulting in injury to the intestinal epithelia, chloragogenous tissues and coelom tissues. Furthermore, As-contaminated soil induced a significant increase in betaine levels and a decrease in dimethylglycine and myo-inositol levels in the earthworms, suggesting that the osmoregulatory metabolism of the earthworms may have been disturbed. The significantly altered levels of asparagine and dimethylglycine were proposed as potential biomarkers of As-contaminated soil. The upregulation of soluble calcium-binding proteins and profilin, the downregulation of sodium/potassium-transporting ATPase, and the proteins changes identified by gene ontology enrichment analysis confirmed that the earthworms suffered from osmotic stress. In addition, the significant changes in glycine-tRNA ligase activity and coelomic tissue injury revealed that As accumulation may disturb the earthworm immune system. This work provided new insight into the proteomic and metabolic toxicity of low-level As-contaminated soil ecosystems in earthworms, extended our knowledge of dual omics and highlighted the mechanisms underlying toxicity.

Metals source apportionment in farmland soil and the prediction of metal transfer in the soil-rice-human chain

China is facing serious heavy metal pollution in farmland soil, which is a major pressing concern for food safety and human health. This research delivers an integrated methodology for pollution source apportionment and a soil-rice-human model to predict heavy metal transfer in the farmland soil, rice grain and human blood chain. The source identification integrated positive matrix factorization (PMF), cluster analysis (CA) and the life cycle assessment (LCA) survey of agricultural and industrial production and consumption. Based on the case analysis of Shaoxing, this method showed very good performance through the illustration of the source contributions by PMF and LCA at county level and the identification of the pollution sources using CA and LCA at field scale. According to the overall evaluation, the integrated method was superior for the farmland metals pollution source identification comparing to existing source apportionment methods. To predict metal transformation in soil-rice-human chain, a set of models of metals (As, Pb, Cd, Hg, Cr) accumulation ability in rice grain and human blood has been established by literature review and monitoring data. The models showed adequate predictability for the metal content of rice grains at both the field and regional scale, and plausible simulation of the metal concentration in human blood throughout the whole study region. Therefore, this study provides valuable tools for farmland soil heavy metal pollution source identification and for the prediction of heavy metal transformation in soil-rice-human chain; and it can highlight the need to take mitigating action to reduce farmland metal pollution risks in specific regions.

Identification of cadmium bioaccumulation in rice (Oryza sativa L.) by the soil-plant transfer model and species sensitivity distribution

Contamination of agricultural soil with cadmium (Cd) poses a severe threat to food safety and human health, especially for Cd in rice. It is very important to identify Cd bioaccumulation in rice in order to screen Cd-safe cultivars. In the present study, 183 pairs of rice and soil data collected from Cd-contaminated soil were used to investigate the differences of Cd bioaccumulation in grains among rice cultivars. The results showed that the adverse effect on grain Cd accumulation of japonica was less than that of indica under Cd exposure. The percentage of japonica with grain Cd concentration exceeding 0.2 mg/kg reduced 50.3% compared with indica. Partial correlation analyses suggested that lower pH contributed to Cd accumulation in grains, and a significant increase in grain Cd concentration was observed with increasing soil Cd concentration. The bioaccumulation factors (BCF) of Cd in rice grains could be divided into 5 grades by combining an empirical soil-plant transfer model with species sensitivity distribution (SSD). Grades with lower Cd bioaccumulation (grades 1 and 2) were dominated by japonica, and the intrinsic sensitivity index of Cd-enrichment (k value) and straw to grain transfer factors (TF) increased with ascending grades. Average k value and TF of cultivars in grade 5 were 1.4-7.9 and 1.5-5.7 times higher than those of cultivars in grades 1 to 4, which eventually caused the increase of Cd accumulation in grains. The lower level of Cd absorption and translocation contributed to reducing the bioaccumulation of Cd in rice grains had been proved by the classification of rice on Cd accumulation. Considering the influence of soil properties and intrinsic sensitivity of rice, cultivars with grain Cd bioaccumulation controlled at low levels to safe for human consumption could be identified on Cd-contaminated soils.

Uptake effects of toxic heavy metals from growth soils into jujube and persimmon of China

Compared with other agricultural plants, information about uptake effects of toxic heavy metals from growth soils into persimmon and jujube are scarce. In this study, the single and comprehensive uptake effects of five toxic heavy metals (Pb, As, Cd, Cr, and Hg) between the two fruits and their growth soils were investigated. The results showed that the average concentrations of heavy metals in the two fruits were found to be 30 (Pb), 6.6 (As), 2.3 (Cd), 38 (Cr), and 0.33 (Hg) μg/kg, respectively. The average concentrations of heavy metals in their growth soils were 26.31 (Pb), 9.63 (As), 0.12 (Cd), 57.6 (Cr), and 0.049 (Hg) mg/kg, respectively. An uptake effect was observed for the two fruits. The values of Nemerow pollution index (NPI) in the two fruits and their growth soils were 0.10 and 0.55, respectively. The average bioconcentration factor (BCF) values of Pb, As, Cd, Cr, and Hg in the two fruits were 0.0012, 0.00075, 0.021, 0.00077, and 0.012, respectively. Based on the residue levels of toxic heavy metals in the growth soils and soil parameters, the prediction models for NPI and BCF were established, with the adjusted regression coefficients of 0.65 (NPI) and 0.81 (BCF). The contribution rates of different soil parameters to NPI were 21.7% (OC), 16.1% (Pb), 17.1% (Cr), 19.8% (Cd), and 25.4% (As), respectively. The contribution rates of different soil parameters to BCF were 10.2% (OC), 9.4% (Cr), and 80.4% (Cd), respectively.

Cadmium phytoavailability under greenhouse vegetable production system measured by diffusive gradients in thin films (DGT) and its implications for the soil threshold

The diffusive gradients in thin films (DGT) technique is recognized to have advantages over traditional techniques. For example, the passive measurement generally follows the principle of metal uptake by plants, and its result incorporates the influences of soil properties, which may make DGT a good protocol for improving soil quality guidelines (SQGs). However, DGT has rarely been applied to assess Cd phytoavailability in soils under greenhouse vegetable production (GVP) systems. In this study, 29 turnips (Raphanussativus L.), 21 eggplants (Solanum melongena L.) and their corresponding soils were collected from GVP systems in Dongtai and Shouguang, eastern China. Simple linear regression and stepwise regression were performed using the soil Cd content and soil properties to predict the vegetable Cd content. Soil thresholds were derived based on both total and available Cd concentrations. The results showed that total Cd, DGT-measured Cd (DGT-Cd), soil-solution Cd (Soln-Cd) and CaCl₂-extractable Cd (CaCl₂-Cd) were all significantly correlated with vegetable Cd. DGT-Cd had the best correlation with turnip Cd. The total Cd threshold values ranged from 4.87 (pH 6.5) to 5.18 (pH 7.5) mg kg^-1 for turnips and 14.60 (pH 6.5) to 14.90 (pH 7.5) mg kg^-1 for eggplants. These Cd thresholds were higher than the current SQGs. The predicted of turnip Cd by DGT-Cd was not improved significantly by further considering the soil properties. The calculated soil threshold of DGT-Cd was 5.35 μg L^-1 for turnips. However, the predicted soil threshold of DGT-Cd for eggplant was improved by including SOM, with R² values from 0.53 to 0.70. The DGT-Cd threshold was calculated as 1.81 μg L^-1 for eggplant (30.0 g kg^-1 SOM). In conclusion, whether DGT measurements are independent of soil properties and preferable for the evaluation of Cd phytoavailability and the generation of soil thresholds remains to be clarified in future research.