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

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.

Vitamin C Alleviates the Risk of Gestational Diabetes Mellitus Associated With Exposure to Metals

Background: Exposure to heavy metals has been suggested to increase the risk of gestational diabetes mellitus (GDM) through the oxidative stress pathway. The study is aimed at examining whether vitamin C could modify the association between exposure to heavy metals and risk of GDM. Methods: We conducted a case-control study in Taiyuan, China, with 776 GDM cases and 776 controls. Data on vitamin C intake from diet and supplements were collected through questionnaires. Concentrations of metals in participants' blood were measured using inductively coupled plasma-mass spectrometry (ICP-MS). Unconditional logistic regression models were applied to estimate effect modification of vitamin C on the association between heavy metals and GDM. Results: Women with higher blood levels of mercury (Hg) (odds ratio (OR) = 2.36, 95% confidence interval (CI): 1.43, 3.92 and 2.04, 95% CI: 1.20, 3.46 for the second and third vs. the first tertile) and arsenic (As) (OR = 2.46, 95% CI: 1.37, 4.43 and 2.16, 95% CI: 1.12, 4.17 for the second and third vs. the first tertile) exposure were associated with increased risk of GDM among women without vitamin C supplement use and having dietary vitamin C intake < 85 mg/day. We found no significant association with metals among women who took vitamin C supplements and/or dietary vitamin C ≥ 85 mg/day. Significant interactions were observed between vitamin C and exposures to metals (i.e., Hg and As) on the risk of GDM (P interaction = 0.048 and 0.045, respectively). Conclusions: Our study, for the first time, suggests that vitamin C supplement use or higher dietary vitamin C intake during preconception and early pregnancy could alleviate the risk of GDM associated with exposure to As and Hg. The results warrant further investigation.

Deriving ecological risk thresholds for soil molybdenum in China based on interspecies correlation estimation and quantitative ion character-activity relationship models

Soil molybdenum (Mo) levels can reach ecologically hazardous levels. China has not yet established the relevant thresholds, posing challenges for environmental management. Therefore, we present our data relevant to Mo toxicity for several important species. By normalizing soil properties, we obtained a correlation model of Mo toxicity to Hordeum vulgare, as well as 31 models for the toxicity of other elements including Cu and Ni to invertebrates and microbial processes. Using interspecies correlation estimation (ICE) extrapolation, the sensitivity coefficient (0.12-0.71) for five plants were found. For invertebrates and microbial processes lacking Mo data, we used regression analysis to establish Mo toxicity models based on the soil quantitative ion character-activity relationships (s-QICAR; R2 =0.70-0.95) and known toxicities of other metal elements to invertebrate and microbial processes. Furthermore, combining species sensitivity distribution calculations, the HC5 values for protecting 95% of soil species from Mo in three typical soil scenarios in China were calculated. After correction, the predicted no-effect concentrations were 6.8, 4.8, and 3.4 mg/kg, respectively. This study innovatively combined ICE and s - QICAR to derive soil Mo thresholds. Our results can provide a basis for decision-making in the assessment and management of soil Mo pollution.

Efficient reduction of Cr(VI) by guava (Psidium guajava) leaf extract and its mitigation effect on Cr toxicity in rice seedlings

Hexavalent chromium (Cr(VI)) is a toxic element that has negative impacts on crop growth and yield. Using plant extracts to convert toxic Cr(VI) into less toxic Cr(III) may be a more favorable option compared to chemical reducing agents. In this study, the potential effects and mechanisms of using an aqueous extract of Psidium guajava L. leaves (AEP) in reducing Cr(VI) toxicity in rice were comprehensively studied. Firstly, the reducing power of AEP for Cr(VI) was confirmed by the cyclic voltammetry combined with X-ray photoelectron spectroscopy (XPS) assays. The highest Cr(VI) reduction efficiency reached approximately 78% under 1.5 mg gallic acid equivalent (GAE)/mL of AEP and 10 mg/L Cr(VI) condition. Additionally, Cr(VI) stress had a significant inhibitory effect on rice growth. However, the exogenous application of AEP alleviated the growth inhibition and oxidative damage of rice under Cr(VI) stress by increasing the activity and level of enzymatic and non-enzymatic antioxidants. Furthermore, the addition of AEP restored the ultrastructure of root cells, promoted Cr adsorption onto root cell walls, and limited the translocation Cr to shoots. In shoots, AEP application also triggered the expression of specific genes involved in Cr defense and detoxification response, including photosynthesis pathways, antioxidant systems, flavonoids biosynthesis, and plant hormone signal transduction. These results suggest that AEP is an efficient reduction agent for Cr(VI), and exogenous application of AEP may be a promising strategy to mitigate the harm of Cr(VI) on rice, ultimately contributing to improved crop yield in Cr-contaminated environments.

Resistance mechanisms of cereal plants and rhizosphere soil microbial communities to chromium stress

Agricultural soils contaminated with heavy metals poison crops and disturb the normal functioning of rhizosphere microbial communities. Different crops and rhizosphere microbial communities exhibit different heavy metal resistance mechanisms. Here, indoor pot studies were used to assess the mechanisms of grain and soil rhizosphere microbial communities on chromium (Cr) stress. Millet grain variety 'Jingu 21' (Setaria italica) and soil samples were collected prior to control (CK), 6 hours after (Cr_6h), and 6 days following (Cr_6d) Cr stress. Transcriptomic analysis, high-throughput sequencing and quantitative polymerase chain reaction (qPCR) were used for sample determination and data analysis. Cr stress inhibited the expression of genes related to cell division, and photosynthesis in grain plants while stimulating the expression of genes related to DNA replication and repair, in addition to plant defense systems resist Cr stress. In response to chromium stress, rhizosphere soil bacterial and fungal community compositions and diversity changed significantly (p < 0.05). Both bacterial and fungal co-occurrence networks primarily comprised positively correlated edges that would serve to increase community stability. However, bacterial community networks were larger than fungal community networks and were more tightly connected and less modular than fungal networks. The abundances of C/N functional genes exhibited increasing trends with increased Cr exposure. Overall, these results suggest that Cr stress primarily prevented cereal seedlings from completing photosynthesis, cell division, and proliferation while simultaneously triggering plant defense mechanisms to resist the toxic effects of Cr. Soil bacterial and fungal populations exhibited diverse response traits, community-assembly mechanisms, and increased expression of functional genes related to carbon and nitrogen cycling, all of which are likely related to microbial survival during Cr stress. This study provides new insights into resistance mechanisms, microbial community structures, and mechanisms of C/N functional genes responses in cereal plants to heavy metal contaminated agricultural soils. Portions of this text were previously published as part of a preprint (https://www.researchsquare.com/article/rs-2891904/v1).

The environmental risk threshold (HC5) for Cd remediation in Chinese agricultural soils

Given the increasing concern over Cd contamination of agricultural soils in China, reducing the availability of the toxic metal has become an important remedial strategy. However, the lack of a unified evaluation framework complicates the assessment of remediation efficiency of different practices. Here, we evaluated the general extraction method (GEM) of available Cd in nine typical soil types by comparing extraction agents, including CaCl2, EDTA, Mehlich-Ⅲ, HCl and DTPA. The safe grain concentration of different agricultural products from National Food Safety Standards Limits of Contaminants in Food (GB 2762-2022) was then applied to understand soil limited available Cd concentration based on dose-response curves. We also derived environmental risk threshold (HC5) values for Cd remediation in agricultural soils by constructing species sensitivity distribution (SSD) curves. The results showed that Mehlich-Ⅲ best predicted Cd accumulation in crops (with 76.5% of explanation of grain Cd) and was selected as the GEM of soil available Cd for subsequent analyses. The regression coefficient (R2) of dose-response curves fitting between Cd absorption in crop tissues and soil available Cd extracted by GEM based on 30 different crop species varied from 51.0% to 79.5%, and the derived limit concentration of soil available Cd based on standard GB 2762-2022 was 0.18-0.76 mg‧kg-1. An HC5 of 0.19 mg‧kg-1 was then calculated, meaning that a concentration of available Cd in agricultural soil below 0.19 mg‧kg-1 ensures that 95% of agricultural products meet the quality and safety requirements of standard GB 2762-2022. The prediction model was well verified in the field test, indicating that can correctly estimate the soil available Cd based on the content of Cd in plant. This study provides a robust scientific framework for deriving the risk threshold for Cd remediation in agricultural soils and could be quite useful for establishing soil remediation standards.

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.

Spectra metrology for interaction of heavy metals with extracellular polymeric substances (EPS) of Pseudomonas aeruginosa OMCS-1 reveals static quenching and complexation dynamics of EPS with heavy metals

The adsorption behavior and interaction mechanisms of extracellular polymeric substances (EPS) of Pseudomonas aeruginosa OMCS-1 towards chromium (Cr), lead (Pb), and cadmium (Cd) were investigated. EPS-covered (EPS-C) cells exhibited significantly higher (p < 0.0001; two-way ANOVA) removal of Cr (85.58 ± 0.39%), Pb (81.98 ± 1.02%), and Cd (73.88 ± 1%) than EPS-removed (EPS-R) cells. Interactions between EPS-heavy metals were spontaneous (ΔG<0). EPS-Cr(VI) and EPS-Pb(II) binding were exothermic (ΔH<0), while EPS-Cd(II) binding was endothermic (ΔH>0) process. EPS bonded to Pb(II) via inner-sphere complexation by displacement of surrounding water molecules, while EPS-Cr(VI) and EPS-Cd(II) binding occurred through outer-sphere complexation via electrostatic interactions. Increased zeta potential of Cr (29.75%), Pb (41.46%), and Cd (46.83%) treated EPS and unchanged crystallinity (CIXRD=0.13), inferred EPS-metal binding via both electrostatic interactions and complexation mechanism. EPS-metal interaction was predominantly promoted through hydroxyl, amide, carboxyl, and phosphate groups. Metal adsorption deviated EPS protein secondary structures. Strong static quenching mechanism between tryptophan protein-like substances in EPS and heavy metals was evidenced. EPS sequestered heavy metals via complexation with C-O, C-OH, CO/O-C-O, and NH/NH2 groups and ion exchange with -COOH group. This study unveils the fate of Cr, Pb, and Cd on EPS surface and provides insight into the interactions among EPS and metal ions for metal sequestration.

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.

Spatial and variety distributions, risk assessment, and prediction model for heavy metals in rice grains in China

In this study, 6229 brown rice grains from three major rice-producing regions were collected to investigate the spatial and variety distributions of heavy metals in rice grains in China. The potential sources of heavy metals in rice grains were identified using the Pearson correlation matrix and principal component analysis, and the health risks of dietary exposure to heavy metals via rice consumption were assessed using the hazard index (HI) and total carcinogenic risk (TCR) method, respectively. Moreover, 48 paired soil and rice samples from 11 cities were collected to construct a predicting model for Cd accumulation in rice grains using the multiple linear stepwise regression analysis. The results indicated that Cd and Ni were the main heavy metal pollutants in rice grains in China, with approximately 10% of samples exceeding their corresponding maximum allowable limits. The Yangtze River basin had heavier pollution of heavy metals than the Southeast Coastal Region and Northeast Plain, and the indica rice varieties had higher heavy metal accumulation abilities compared with the japonica rice. The Cu, Pb, and Cd mainly originated from anthropogenic sources, while As, Hg, Cr, and Ni originated from both natural and anthropogenic sources. The mean HI and TCR values of dietary exposure to heavy metals via rice consumption ranged from 2.92 to 4.31 and 9.74 × 10-3 to 1.44 × 10-2, respectively, much higher than the acceptable range, and As and Ni were the main contributor to the HI and TCR for Chinese adults and children, respectively. The available Si (ASi), total Cd (TCd), available Mo (AMo), and available S (AS) were the main soil factors determining grain Cd accumulation. A multiple linear stepwise regression model was constructed based on ASi, TCd, AMo, and AS in soils with good accuracy and precision, which could be applied to predict Cd accumulation in rice grains and guide safe rice production in contaminated paddy fields.

Screen for low-arsenic-risk rice varieties based on environment-genotype interactions by using GGE analysis

Arsenic (As) accumulation in rice is a global health concern that has received increased attention in recent years. In this study, 12 rice genotypes were cultivated at four As-contaminated paddy sites in Taiwan. According to the different crop seasons and As levels in the soil, the sites were further divided into 18 environmental conditions. For As in soils, results showed that 67% of the studied environments were likely to represent As contamination. For As in rice, the mean total As concentration in brown rice grains ranged from 0.17 to 0.45 mg kg-1. The analysis of variance for the environment effect indicated that grain As concentration was mainly affected by the environmental conditions, suggesting that there was a remarkable degree of variation across the trial environments. According to the combination of the GGE biplot and cumulative distribution function of order statistics (CDFOS) analysis, five genotypes-TCS17, TCS10, TT30, KH139, and TC192-were regarded as stable, low-risk genotypes because the probability of grain As concentration exceeding the maximum permissible concentration (MPC) was lower for these genotypes across all environmental conditions. Particularly, TCS17 was recommended to be the safest rice genotype. Thus, grain As levels in the selected genotypes were applied to assess the health risk to Taiwanese residents associated with As exposure through rice consumption. Results showed that the upper 75th percentile values of the hazard quotient were all less than unity. This suggested that the health risk associated with consuming the selected rice genotypes was acceptable for most of the residents. The methodology developed here would be applicable to screen for stable, low-As-risk rice genotypes across multiple field environments in other regions or countries.

Health risk assessment of trace metal(loid)s in agricultural soil using an integrated model combining soil-related and plants-accumulation exposures: A case study on Hainan Island, South China

Traditional health risk assessment of trace metal(loid)s (TMs) in agricultural soil exclusively considers direct soil-related exposure and may underestimate the health risks they pose. In this study, the health risks of TMs were evaluated using an integrated model that combined soil-related and plant-accumulating exposures. A detailed investigation of common TMs (Cr, Pb, Cd, As, and Hg) coupled with probability risk analysis based on a Monte Carlo simulation was conducted on Hainan Island. Our results showed that, except for As, the non-carcinogenic risk (HI) and carcinogenic risk (CR) of the TMs were all within the acceptable ranges (HI < 1.0, and CR < 1E-06) for direct soil-related exposure to bio-accessible fractions and indirect exposure via plant accumulation (CR substantially lower than the warning threshold 1E-04). We identified crop food ingestion as the essential pathway for TM exposure and As as the critical toxic element in terms of risk control. Moreover, we determined that RfDo and SFo are the most suitable parameters for assessing As health risk severity. Our study demonstrated that the proposed integrated model combining soil-related and plant-accumulating exposures can avoid major health risk assessment deviations. The results obtained and the integrated model proposed in this study can facilitate future multi-pathway exposure research and could be the basis for determining agricultural soil quality criteria in tropical areas.

Toxicity of chromium to wheat (Triticum aestivum L.) in two soils: influence of soil properties and chromium form

The toxicity of Cr to plants depends on Cr form and soil properties. Currently, the phytotoxicity differences of Cr(VI) and Cr(III) in different soils are not clear. In this study, the toxicity of Cr(VI) and Cr(III) to root growth and root morphology of wheat (Triticum aestivum L.) were compared in Shandong fluvo-aquic soil (SD soil) and Jiangxi red soil (JX soil) that is differing in soil properties. The toxicity thresholds of Cr(VI) and Cr(III) on wheat root elongation were determined by fitting the dose-effect curves. Results showed that the 10% and 50% root length inhibitory concentrations (EC10 and EC50) of Cr(III) were 53.1 and 125 times of Cr(VI) in SD soil and 8.11 and 1.36 times of Cr(VI) in JX soil, indicating that Cr(VI) was more toxic to wheat roots than Cr(III) in both soils and the toxicity discrepancy of the two forms of Cr was more prominent in SD soil. Cr(VI) exhibited higher toxicity in SD soil (alkaline) than in JX soil (acidic), whereas Cr(III) showed the opposite pattern. In addition, the ethylene diamine tetraacetic acid extractable Cr (EDTA-Cr) concentrations in soils were correlated well with the relative wheat root elongation (R2=0.854, P<0.01), indicating that soil EDTA-Cr concentration can be used as a predictor of Cr phytotoxicity. Both Cr(VI) and Cr(III) showed significant biphasic dose effects on wheat root morphology (root length, root surface area, root volume, and root tip number) in JX soil. These findings are helpful for the risk evaluation of Cr contamination in agricultural soils.

In-situ conversion of geopolymer into novel floral magnetic sodalite microspheres for efficient removal of Cd(II) from water

In the present work, a novel, floral-like, magnetic sodalite microsphere (SODM) was synthesized in situ by using fly ash (FA) and metakaolin (MK) as raw materials and was used to remove Cd(II) from water. Its magnetism can solve the problems of adsorbent recovery and possible secondary pollution. During the static adsorption, SODM shows a maximum adsorption capacity of 245.17 mg/g. The adsorption of Cd(II) on the SODM surface is spontaneous, exothermic, and physicochemical adsorption, which was evaluated by thermodynamics, kinetics, and isotherm studies. During dynamic adsorption, SODM shows a maximum adsorption capacity of 342.74 mg/g in the simulated solution prepared by the deionized water, compared to 215.88 mg/g in the simulated solution prepared using Xiangsi Lake water from Guangxi Minzu University. At 0.5 g SODM dosage in the dynamic adsorption, the adsorption capacity could rise to 632.81 mg/g. These results demonstrated the excellent Cd (II) adsorption performance of the SODM. The adsorption of cadmium on the SODM surface includes the synergistic effects of electrostatic attraction, ion exchange, and surface coordination reaction. Besides, the SODM shows good regeneration performance in both the deionized water and Xiangsi Lake water. The present study explores SODM as an adsorbent for the Cd (II) removal from wastewater and unbolts the industrial applicability of the SODM in the field of wastewater purification.

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.

Rapid identification of high and low cadmium (Cd) accumulating rice cultivars using machine learning models with molecular markers and soil Cd levels as input data

Excessive accumulation of cadmium (Cd) in rice grains threatens food safety and human health. Growing low Cd accumulating rice cultivars is an effective approach to produce low-Cd rice. However, field screening of low-Cd rice cultivars is laborious, time-consuming, and subjected to the influence of environment × genotype interactions. In the present study, we investigated whether machine learning-based methods incorporating genotype and soil Cd concentration can identify high and low-Cd accumulating rice cultivars. One hundred and sixty-seven locally adapted high-yielding rice cultivars were grown in three fields with different soil Cd levels and genotyped using four molecular markers related to grain Cd accumulation. We identified sixteen cultivars as stable low-Cd accumulators with grain Cd concentrations below the 0.2 mg kg^-1 food safety limit in all three paddy fields. In addition, we developed eight machine learning-based models to predict low- and high-Cd accumulating rice cultivars with genotypes and soil Cd levels as input data. The optimized model classifies low- or high-Cd cultivars (i.e., the grain Cd concentration below or above 0.2 mg kg^-1) with an overall accuracy of 76%. These results indicate that machine learning-based classification models constructed with molecular markers and soil Cd levels can quickly and accurately identify the high- and low-Cd accumulating rice cultivars.

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.

Accumulation of Cr in different vegetables and derivation of soil Cr threshold using the species sensitivity distribution method

Due to its high mobility and bioavailability, hexavalent chromium [Cr(VI)] in agricultural soil can be taken up by crops and pose threat to human being. In this study, two soils (Jiangxi red soil and Shandong fluvo-aquic soil) spiked with Cr(VI) and 8 common vegetable varieties were used to conduct the pot experiment. The bioconcentration factor (BCF) values based on the tetraacetic acid extractable Cr (EDTA-Cr) in soils were used to construct the species sensitivity distribution (SSD) curve. Afterwards, the soil Cr threshold was derived based on the critical BCF value and the permissible limit of Cr for vegetables. The results showed that when spiked with 5.6 mg kg^-1 of Cr(Ⅵ), the soil EDTA-Cr concentrations were significantly increased compared with the control except Jiangxi red soil planted with carrot and radish, while the Cr concentrations in the edible parts of vegetables in both soils were below the permissible limit (0.5 mg kg^-1 FW). However, there are dramatic differences in the accumulation of Cr by different varieties of vegetables. Apparent discrepancy was observed between the two soils for the bioconcentration of Cr by carrot. Among the leafy vegetables, lettuce and oilseed rape are the most and the least sensitive to Cr pollution, respectively. The safety threshold values of EDTA-Cr were 0.70 mg kg^-1 for Shandong fluvo-aquic soil and 0.85 mg kg^-1 for Jiangxi red soil, respectively. This study provides information on the safety production of vegetable products in Cr(Ⅵ) polluted soils and is helpful to the revision of soil quality standards of Cr.

Human health risks associated with metals in paddy plant (Oryza sativa) based on target hazard quotient and target cancer risk

Paddy plants (Oryza sativa) contaminated with metals could be detrimental to human health if the concentrations of metals exceed the permissible limit. Thus, this study aims to assess the risk of the concentrations of As, Se, Cu, Cr, Co, and Ni and their distributions in various parts (roots, stems, leaves, and grains) of paddy plants collected from Sekinchan, Malaysia. Both soil and plant samples were digested according to the United States Environmental Protection Agency (USEPA) Method 3050B and the metal concentrations were determined by the Inductively Coupled Plasma-Mass Spectrometry (ICP-MS). The highest mean translocation factor (TF) was from soil to roots (TF roots/soil ranged from 0.12 to 6.15) and the lowest was from leaves to grain (TF grain/leaves ranged from 0.06 to 0.87). Meanwhile, the bioaccumulation factor (BAF) for all metals was less than 1.0 indicating that paddy plants only absorb metals from the soil but do not accumulate in the grains. The average daily intake for As (1.15 ± 0.25 µg/kg/day) has exceeded the limit proposed by ATSDR and IRIS USEPA (0.30 µg/kg/day). Target cancer risk (TR) of 1.10 × 10^-3 for As through rice consumption indicates that the potential cancer risk exists in one out of 1000 exposed individuals. The results from this study could serve as a reference for researchers and policymakers to monitor and formulate strategies in managing As and other metals in paddy plants, especially in Southeast Asian countries.

A New Method for Ecological Risk Assessment of Combined Contaminated Soil

Ecological risk assessment of combined polluted soil has been conducted mostly on the basis of the risk screening value (RSV) of a single pollutant. However, due to its defects, this method is not accurate enough. Not only were the effects of soil properties neglected, but the interactions among different pollutants were also overlooked. In this study, the ecological risks of 22 soils collected from four smelting sites were assessed by toxicity tests using soil invertebrates (Eisenia fetida, Folsomia candida, Caenorhabditis elegans) as subjects. Besides a risk assessment based on RSVs, a new method was developed and applied. A toxicity effect index (EI) was introduced to normalize the toxicity effects of different toxicity endpoints, rendering assessments comparable based on different toxicity endpoints. Additionally, an assessment method of ecological risk probability (RP), based on the cumulative probability distribution of EI, was established. Significant correlation was found between EI-based RP and the RSV-based Nemerow ecological risk index (NRI) (p < 0.05). In addition, the new method can visually present the probability distribution of different toxicity endpoints, which is conducive to aiding risk managers in establishing more reasonable risk management plans to protect key species. The new method is expected to be combined with a complex dose-effect relationship prediction model constructed by machine learning algorithm, providing a new method and idea for the ecological risk assessment of combined contaminated soil.

Heavy metal(loid)s in agriculture soils, rice, and wheat across China: Status assessment and spatiotemporal analysis

Heavy metal(loid)s (HMs) accumulation in agricultural soils, rice, and wheat is of particular concern in China, while the status and spatio-temporal distribution of HMs in the soil-crops system have been rarely reported at the national scale. This study aimed to summarize the overall pollution status, spatiotemporal patterns, and drivers of HMs in agricultural soil, rice, and wheat nationwide. The metal-polluted data from 1030 agricultural soils, rice, and wheat in China were collected from the literature published from 2000 to 2022. The results showed that Cd was the most prevailing contaminant in soils based on its spatiotemporal distribution and accumulation. The pollution cases and severe pollution percentage of Cd (103 %) and Hg (128 %) show an increasing trend pattern. Mining activities are the main anthropogenic sources of agricultural soil HMs in China. Cd and Pb had the highest exceedance rate in rice (33.5 and 32.2 %) and wheat (25.8 and 30.3 %). The rice from Hunan, Fujian, and Guangxi showed the highest average concentration of Cd and Pb, respectively, while wheat samples from Hubei had the greatest exceedance rate of Pb. Besides, HMs in crops was not usually corresponding to soil HMs but increased gradually from north to south areas. Several mitigation strategies and accurate health risk assessments model of HMs based on bioavailability were also proposed and recommended. Collectively, this review provides valuable information to improve the management of farmland nationwide, optimize the accurate risk assessment, and reduce HMs pollution.

Coupling phytotoxicity and human health risk assessment to refine the soil quality standard for As in farmlands

In the present study, a field experiment was conducted to investigate arsenic (As) concentrations in soils and in grains of 15 rice varieties in a contaminated site in Taiwan. The studied site was divided into two experimental units, namely plot A and plot B. The results showed that mean total As concentrations were 70.94 and 61.80 mg kg^-1 in plot A and plot B, respectively, and thus greater than or approximate to the soil quality standard for total As in Taiwan (60 mg kg^-1). The As levels in rhizosphere soil in plot A (19.71-32.33 mg kg^-1) were much higher than in plot B (6.41-8.60 mg kg^-1); however, As accumulation in brown rice did not significantly differ between the plots. These results implied that a significant variation in the bioconcentration factor (BCF) value of As existed among different rice genotypes, and a negative correlation was observed between BCF value and rhizosphere As level in the soil. In phytotoxicity, the median values of the ecological risk indicator were 104.85 and 103.89 in plot A and plot B, respectively, indicating considerable risk. In human health risk assessment, the median and 97.5%-tile values for cancer risk for both male and female residents were markedly higher than the acceptable risk (1 × 10^-4). Furthermore, non-cancer and cancer risks were higher for males than females, mainly due to the greater rice ingestion rate of males. Sensitivity analysis showed that total As concentration in soil was the main factor affecting health risks, suggesting that priority should be given to the reduction of soil As levels. To better manage the phytotoxicity of As on rice, as well as the health risk to residents resulting from exposure to As-contaminated soils, the soil quality standard for As in farmland soils should be set between 5 and 10 mg kg^-1. The methodology developed in this study could also be applied to provide the basis for refining and revising the soil quality standard for heavy metals in farmland in other regions and countries.

Ecological risk threshold for Pb in Chinese soils

Derivation of ecological risk threshold (the threshold concentration value that protect a certain proportion of species within the acceptable hazard level) of lead (Pb) is a yardstick and plays a key role in formulating soil protection policies, while the research about deducing soil Pb ecological risk threshold is still limited. In this study, toxicological data of Pb based on 30 different test endpoints was collected from our experiment and literature, and applied into interspecific extrapolation by species sensitivity distribution (SSD) method to derive the hazard concentration for 5% of species (HC₅, that can protect 95% of species), the prediction models according to different soil properties were established. The results showed that EC₁₀ (the effective concentrations of Pb that inhibit 10% of endpoint bioactivity) ranged from 205.6 to 1596.3 mg kg¹, and hormesis induced by Pb were up to 118%. Toxicity data were corrected by leaching and aging process before SSD curves fitting. HC₅ was then derived and prediction model was developed, as LogHC₅ = 0.134 pH + 0.315 LogOC + 0.324 LogCEC + 1.077. The prediction model was well verified in the field test, indicating that can correctly estimate Pb ecotoxicity thresholds in different soils. This study provides a scientific frame for deriving the ecological risk threshold of Pb and is of great significance for ecological species protection.

Prediction model for Cd accumulation of rice (Oryza sativa L.) based on extractable Cd in soils and prediction for high Cd-risk regions of southern Zhejiang Province, China

Soil environmental quality in China for agricultural land always considers the effect of total cadmium (Cd) in soil, ignoring the bioavailability of soil Cd. The 139 paired rice (Oryza sativa L.) and soil samples were collected from the Cd-contaminated paddy fields of southern Zhejiang Province, China. The objectives of this study were to establish accurate prediction models for Cd accumulation in brown rice based on bioavailable Cd and physiochemical properties of soils and to evaluate the safety of rice production in Cd-contaminated paddy. The bioavailable Cd in soils was extracted and evaluated by using CaCl₂, HNO₃, diethylenetriamine pentaacetic acid (DTPA), diffusive gradients in thin-films technique (DGT), and sequential extraction method proposed by the European Community Bureau of Reference; 100 pairs of data were used as training sets, and the remaining 39 sets were used as validation sets. Stepwise multiple linear regression analysis and random forest analysis showed that total Cd in soil could roughly indicate the content of Cd in rice, while extractable Cd could better explain the accumulation of Cd in rice grain and DTPA and DGT extractive technologies are the most evaluative. The validation sets also showed that the prediction model has a good fit. Based on the prediction model for Cd in brown rice based on soil pH and DGT extractive Cd, the Monte Carlo simulation showed that 74.32% and 89.35% of the estimated brown rice hazard quotient (HQ) of the daily Cd intake of adults and children in safe utilization paddy sites could exceed the safe level of 1, respectively. Additionally, the threshold values for extractable Cd by DGT or DTPA for rice safe production were 3.4 μg/kg or 0.13 mg/kg when the pH in soils was below 5.5. The results further proved the threshold concentration of extractable Cd for predicting high-risk soils of Cd contamination in brown rice.

Assessment soil cadmium and copper toxicity on barley growth and the influencing soil properties in subtropical agricultural soils

Evaluation joint cadmium (Cd) and copper (Cu) phytotoxicity in wide range of subtropical agricultural soils is highly vital for phytoremediation of soils contaminated with Cd and Cu. In this study, barley root elongation assays were performed in 30 representative soils in response to single and combined Cd and Cu inhibition. The single Cd caused nearly 50% inhibition of barley root elongation, and Cu induced more than 50% inhibition in most soils. Mixed Cd + Cu caused significant inhibition on barley growth with average relative root elongation values of 20.0% and 30.4% in soil with a pH < 7 and pH > 7, respectively. An antagonistic interaction was evaluated in combined Cd + Cu toxicity, which was strong in soils containing low soluble Cu and Cd contents. Soil pH was the controlling factor in predicting single and mixed Cd and Cu phytotoxicity, which could explain 44% and 46% variation of single Cd and Cu toxicity, respectively. Soil organic carbon and effective cation exchange capacity were another important factor positively influencing metal toxicity, which further improved empirical prediction models accuracy, with determined coefficient (r²) values of 0.44-0.84. These results provide a theoretical basis for soils Cd and Cu pollution control.

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.

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 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.