Determination of Soil Cadmium Threshold for Potato

In this study, pot experiments were conducted to investigate the characteristics of Cd transfer to potato tubers from two types of soil. The results showed that the Log-normal and Burr III functions can be used to determine the sensitivity of different potato varieties to Cd as well as the soil Cd threshold. With regard to the prediction accuracy, the root mean squared error (RMSE) values for the total Cd biaoaccumulation factor (BAF_total) calculated with both functions were smaller than those for exogenous Cd biaoaccumulation (BAF_add) in acidic and alkaline soils, indicating that BAF_total is more appropriate for the calculation of the soil Cd threshold. The average Cd threshold values in acidic soil calculated with the Log-normal and Burr III functions were 0.411 and 0.461 mg kg⁻¹, and the average values in alkaline soil were 0.716 and 0.888, respectively. The Log-normal function can also be applied to fit the sensitivity distributions of different species for the development of appropriate soil Cd threshold values for conservation purposes.

Introduction

Cadmium (Cd) is a toxic heavy metal that is harmful to plants, microorganisms and humans (Zhong, et al., 2015). The evaluation of soil Cd levels and Cd transfer from the soil to plants is therefore a critical issue. Potato is the fourth most important crop in China, behind wheat, rice and corn. The ecology of Karst areas is fragile and human activities are more likely to cause soil Cd pollution in these areas. Therefore, it is necessary to pay more attention to Cd transfer from the soil to crops in Karst areas (Zhang, et al., 2020).

The Cd content in different organs generally increases in the order of leaves > stem/root > tuber (Chen, et al., 2014). The Cd content in the tuber is less than that in other organs and most of the Cd is retained by other organs (Ye, et al., 2020). Soil properties such as pH or organic matter (OM) markedly influence the distribution of metals between the phytoavailable and non-phytoavailable sections of a plant (Pinto, et al., 2015); (Hough, et al., 2004). Other metal ions can also affect plant absorption of Cd from the soil; for instance, correlation analysis of different potato varieties showed a strong positive relation between the bioaccumulation of Cd and the presence of zinc (Zn), copper (Cu) or manganese (Mn) (Ashrafzadeh, et al., 2017). The plant species also affect soil Cd uptake. For example, the tomato genotype affects the nutritional value, fruit quality, yield, etc, and Cd-tolerant cultivars were found to accumulate more Cd than sensitive cultivars (Carvalho, et al., 2018). The results of these studies confirmed that Cd uptake by plants is influenced by various factors including the soil properties and crop varieties. Biaoaccumulation factor (BAF) is used to study the transfer of a pollutant from the soil to plants and to calculate the threshold value for the pollutant in the soil to ensure food safety. A prediction model for Pb transfer to wheat from 17 soil types with different Pb levels in China was fitted using the BAF and soil properties to determine the soil Pb threshold (Liu, et al., 2016). Some scholars developed and compared regression models based on BAF for the uptake of Cd by potatoes through field investigations (Novotna, et al., 2015).

The cumulative probability distribution method have been used to evaluate the risk of pollutants in humans and are helpful to establish the soil threshold for heavy metals as well as for soil screening. The species sensitivity distribution (SSD) is a typical cumulative probability distribution that reflects the differences in the sensitivity of different species to the same pollutant (Liu, et al., 2015). Environmental toxicity data for pollutants such as the concentration necessary to achieve 50% of the maximal effect (EC₅₀), no observed effect concentration (NOEC) and BAF for a limited number of species obtained through biological toxicity tests are sampled based on the SSD. Establishing a threshold that covers the fifth percentile of the SSD protects the majority of the species in an ecosystem (van der Hoeven, 2001). This strategy is applied more frequently to water environments; for example, the SSD for nine compounds indicates that crustaceans are the most sensitive species to anticancer drugs in aquatic ecosystems (Li D et al., 2021). Furthermore, it has been increasingly applied to soil environments in recent years. A soil-plant transfer model was applied to Cd data for 183 soil-plant pairs to fit the SSD and 5 grades of Cd sensitivity were obtained based on the Cd-BAF (Li, et al., 2019). Many function models are available for SSD fitting, including the Burr III, Log-normal, Log-logistic and Gamma functions. The Gaussian distribution function was used to determine the cumulative probability distribution for vegetables and the soil Cd threshold in the Hunan Province of southern China as well as to implement a reasonable soil Cd risk assessment strategy for the locality (Yang et al., 2016a). The potential health risk to humans from Cd exposure was evaluated using the target hazard quotient (THQ) and the probability distribution of THQ for 324 vegetable samples in Shandong, Jiangsu, Yunnan Province in China was identified with the Log-normal function (Hu, et al., 2018). The Cd toxicity threshold for paddy soil was determined based on the variation of BAF in various rice cultivars using SSD fitted with the Burr III function (Song, et al., 2015). The appropriate model depends on the specific situation in terms of the polluted species and the regulation of the environment.

In the present study, we aimed to: 1) explore the transfer of Cd from two types of soil with different pH values to potato tubers of different cultivars; 2) compare the fit of SSD curves using the Log-normal and Burr III functions for different potato cultivars in the two soil types; and 3) obtain a soil threshold value for the safe production of potato using different calculation methods for BAF.

Materials and Methods

Pot Experiments

Soil was collected from two Karst geological regions in China, specifically from Weining and Guiyang cities in Guizhou Province. Soil samples were collected from the surface of a field (0–20 cm) in these regions, air dried and passed through a 2 mm sieve to remove large pieces of gravel prior to the pot experiments. The base Cd concentrations were 0.61 and 0.37 mg kg⁻¹ for the Weining (yellow) and Guiyang (calcareous) soil, respectively. The two soils are the main soil types in Guizhou province. Soil properties were determined using standard methods. Soil pH was determined using a glass electrode with a soil:water ratio of 1:2.5 (g ml⁻¹). The pH values for the Weining and Guiyang soil were 4.8 and 7.8. The measured OM values were 24.5 and 48.5 g kg⁻¹, respectively, using the potassium dichromate oxidation method. Ten commercial potato varieties commonly grown in the local area, namely Eeshu 5 (E5), Lishu 15 (L15), Qianyu 8 (QY8), Qingshu 9 (QS9), Weiyu 5 (W5), Lishu 13 (L13), Hongbaoshi (HBS), Heimeiren (HMR), ChuangyuA5e (CA5e) and Weiyu 7 (W7), were used in the experiments.

The potato varieties were planted in the two soil types in the way of pot experiment (greenhouse). Six kilograms of soil was placed in a plastic pot to ensure the normal growth of the potatoes. Exogenous Cd in solution (3CdSO₄ 8H₂O) was added to the soil and mixed uniformly. Two treatments were applied, with three replicates for each pot. The control treatment (CK) lacked additional Cd, and 0.9 mg kg⁻¹ Cd was mixed with the two soil types for the Cd treatment. The amount of Cd added was based on the risk value for Chinese Cd soil contamination in agricultural land (0.3 mg kg⁻¹ for soil pH≦7.5, 0.6 mg kg⁻¹ for soil pH > 7.5). Additionally, 0.15 g kg⁻¹N (CO(NH₂)₂), 0.05 g kg⁻¹ P(Ca (H₂PO₄)₂) and 0.1 g kg⁻¹ K (K₂SO₄) were added to each pot as the base fertilizer. The treated soil was left to equilibrate for 40 days to allow the behavior of Cd ions in the soil to stabilize. Potatoes with similar sizes and weights were selected for planting in the pot experiments. The potatoes were watered with deionized water to ensure that the Cd absorbed by the plants came from the soil. The soil moisture content was maintained at about 70% of the water holding capacity.

Detection of Cd in Soil and Plants

After passing the soil through a 100 mesh sieve, 0.1 g was weighed (accurate to 0.0001) and transferred to a polytetrafluoroethylene (PTFE) tank. Afterward, 3 ml nitric acid, 1 ml hydrochloric acid and 1 ml hydrofluoric acid were added and the capped tank was placed in a matched steel pipe before digesting for 20 h at 180°C on an electric heating plate. Perchloric acid (1 ml) was added as the acid driver. After digestion of the soil sample, the digestion solution was transferred to a 50 ml volumetric flask and a constant volume was maintained with ultrapure water. Finally, 10 ml of the supernatant from the digested sample solution was collected and the Cd concentration was measured using inductively coupled plasma mass spectrometry (ICP-MS; Thermo Fisher Scientific, Waltham, MA, United Sates, x2). The basic test conditions was injection speed: 0.8 ml min⁻¹; atomization gas flow rate: 0.92 ml min⁻¹; auxiliary gas flow rate: 0.7 ml min⁻¹; cooling gas: 13 L min⁻¹; the detection limit of the instrument is less than 1 ppb.

Potato tubers were collected at maturity and washed three times successively with tap water and purified water to remove soil particles from the surface. The potato tubers were placed in a 105°C electric oven for 15 min to remove green sections and then dried to a constant weight at 75°C. The plant samples were ground and stored for determination of the Cd content. We weighed 0.2 g (accurate to 0.0001) of edible plant samples into a PTFE tank, then added 5 ml nitric acid and two to three drops of hydrofluoric acid for 12 h for pre-digestion. The initial digestion solution was heated to 160 °C and digested for 8 h to remove OM. After cooling, a constant volume of 50 ml was maintained. The determination method for the Cd content in the plants was the same as that for the soil. For quality control, all the acid reagents used were of high grade purity. The glassware and PTFE tube used for digestion were soaked in 25% nitric acid solution for 12 h before use. The soil reference material GSS-5 and plant reference material GBW10021 were adopted for quality control.

Calculation Method for BAF

BAF represents the ratio of the heavy metal concentration in plants to the heavy metal concentration in the soil environment. The BAF for the CK treatment (BAF_ck) was calculated using the equation:

${\text{BAF}}_{\text{ck}}=\raisebox{1ex}{${\text{C}}_{\text{CK}-\text{plant}}$}\!\left/ \!\raisebox{-1ex}{${\text{C}}_{\text{CK}-\text{soil}}$}\right.(1)$

Where, C_CK-plant and C_CK-soil are the heavy metal concentrations in the plants and soil for the CK treatment, respectively.

The BAF for total Cd (BAF_total) was calculated with the equation:

${\text{BAF}}_{\text{total}}=\raisebox{1ex}{${\text{C}}_{\text{plant}}$}\!\left/ \!\raisebox{-1ex}{${\text{C}}_{\text{soil}}$}\right.(2)$

Where, C_plant is the concentration of heavy metals in the plants and C_soil is the concentration of heavy metals in the soil for the Cd treatments.

Some scholars (Ding, et al., 2013) proposed a BAF calculation method for exogenous pollutants that only considers the effects of exogenous pollutants on plants. The BAF for exogenous Cd (BAF_add) was calculated with the equation:

${\text{BAF}}_{\text{add}}=\raisebox{1ex}{$\left({\text{C}}_{\text{plant}}{-\text{C}}_{\text{CK}-\text{plant}}\right)$}\!\left/ \!\raisebox{-1ex}{$\left({\text{C}}_{\text{soil}}{-\text{C}}_{\text{CK}-\text{soil}}\right)$}\right.(3)$

Where, C_plant and C_soil are the concentrations of heavy metals in the plants and soil for the Cd treatments, respectively.

Results and Discussion

BAF for Cd in Two Soil Types

The variability of BAF is important for the assessment of Cd exposure for plants in contaminated soil and represents the current contamination situation (Augustsson, et al., 2015). As seen in Figure 1, the total Cd-BAF for tubers grown in Weining soil varied between 0.09 (HMR) - 0.33 (QS9) for BAF_ck and 0.07 (HMR) - 0.25 (CA5e) for BAF_total. The average BAF_ck and BAF_total values for the 10 potato varieties evaluated in this study were 0.23 and 0.17, respectively. Twenty-one potato samples were collected from a field in the Czech Republic and the average BAF was 0.34, which was greater than the values obtained in the current study because of differences in the soil properties (Novotna, et al., 2015). The BAF_ck for each potato variety was slightly higher than the corresponding BAF values for the Cd treatments. In contrast, the BAF_ck for carrot in 21 soil types from China was less than the BAF_total for the Cd treatment (soil spiked with 0.6 mg kg⁻¹ Cd) (Ding, et al., 2013). The difference in the results may be attributed to the use of different crops and the crop growth conditions. For both the CK and Cd treatments, the HMR potato variety accumulated the least Cd while the CA5e variety accumulated the most. In a previous study, three potato cultivars that easily uptake Cd from the soil, known as Cd accumulators, were identified from 10 cultivars based on the BAF and plant Cd concentration (Ashrafzadeh, et al., 2017). These studies demonstrate that identifying crop genotypes with low or high Cd accumulation based on the enrichment factor is an effective method to avoid Cd pollution in food.

FIGURE 1

FIGURE 1. BAF_ck, BAF_total and BAF_add for different potato varieties in Weining soil.

The BAF_add of Cd for the 10 potato varieties evaluated in this study ranged from 0.04 (HMR) to 0.19 (CA5e), and the average BAF_add was 0.13, which was slightly less than the average BAF_total and BAF_ck. The BAF_total and BAF_ck indicated that the CA5e potato variety could easily uptake soil Cd, while the HMR and HBS varieties absorbed less Cd from the soil. A prediction model was established based on soil Cd content, plant Cd content and soil properties (pH and organic carbon) and the results showed that the BAF_total and BAF_add can be used to develop a prediction model for Cd transfer from the soil to vegetables as well as for identifying soil thresholds for food safety (Liang, et al., 2013). Empirical soil-carrot transfer models were developed for 21 soil types in China and the prediction model for added Cd was found to be better than the model for total Cd (R²: 0.77 > 0.73) (Ding, et al., 2013). These results indicated that BAF_total and BAF_add are both valuable parameters for studying soil Cd uptake in crops and calculating the soil threshold.

The Cd-BAF of tuber grown in Guiyang soil ranged from 0.03 (W5) to 0.08 (L15) for BAF_ck and from 0.05 (L15) to 0.13 (QY8) for BAF_total. The average BAF_ck and BAF_total of the tubers were 0.06 and 0.09 for different potato varieties, respectively. Some researchers studied the accumulation of Cd in potato with different soil Cd levels (soil Cd < 25 mg kg⁻¹) and found that the BAF_total of the tubers ranged from 0.20 to 0.81. The BAF_total values in these studies were larger than our results, probably due to higher Cd pollution (<25 mg kg⁻¹) and low OM content (14.6 g kg⁻¹) in the soil compared to our soil conditions (soil Cd: 0.9 mg kg⁻¹, OM: 48 g kg⁻¹) (Chen, et al., 2014). Most of the BAF_total values were higher than the BAF_ck values for specific varieties in the alkaline Guiyang soil and the opposite results were obtained in the acidic Weining soil. The biggest difference between the BAF_ck and BAF_total was found in the W5 variety, which had a BAF_total value that was 2.74 times greater than the BAF_ck. The BAF_add ranged from 0.02 (L15) to 0.17 (QY8), with an average BAF_add of 0.11, which was close to the BAF_total for the Cd treatments.

The transfer of Cd from the soil to the potato tubers was limited with the HMR and HBS varieties in both the Weining soil and Guiyang soil for the CK and Cd treatments. The QY8 and CA5e varieties were sensitive to Cd uptake from the Guiyang soil, while QY8, QS9 and CA5e were sensitive to Cd uptake from the Weining soil. QY8 and CA5e were sensitive to Cd uptake with the CK and Cd treatments in the Weining and Guiyang soils, as indicated by the high BAF_total and BAF_add. The BAF_add and BAF_total values demonstrated that Cd-sensitive and insensitive potato varieties were similar in the two soil types with different pH values and the genotype affected food safety in the polluted soils.

The BAF_ck and BAF_total values in the alkaline Guiyang soil differed from those in the acidic Weining soil. Figures 1, 2 show that the ability of the potato varieties to absorb Cd in acidic soil was stronger than that in alkaline soil for both the CK and Cd treatments based on the average BAF_total and BAF_ck. pH is an important factor that affects the adsorption, desorption, precipitation, dissolution and complexation of heavy metals in the soil and can directly or indirectly affect the retention of heavy metals in the soil. In a previous study in which the soil pH was increased by 0.50 units through lime addition, the uptake of Cd by rice grain decreased by 35.3% and the Cd transfer ratio increased as the soil pH decreased in acidic and alkaline soils (Zhu, et al., 2016). The results of the current study demonstrated that the soil type and genotype significantly affect Cd uptake by potato and so does the interaction between these two factors. Luyin No. One was identified from four potato cultivars as a suitable cultivar for planting in slightly Cd-contaminated soils (<0.6 mg kg⁻¹). The average reported BAF_add for different cultivars is 2.3, which is less than the reported BAF_total value (2.76) in soil with pH 6.9 and is similar to the results for our acidic Weining soil (BAF_add 0.13 < BAF_total 0.17) (Ding, et al., 2014). In general, there was significant difference between the Weining soil and Guiyang soil in both the BAF_ck and BAF_total, however, BAF_add was different from BAF_ck and BAF_total (Table 1 in the Supplementary Data S1). The BAF_ck, BAF_total and BAF_add of some varieties have significant differences, while other varieties have no significant differences among each other in the two soils, moreover, the significant difference of the BAF_total was similar to that of BAF_add (Table 2 in the Supplementary Data S1).

FIGURE 2

FIGURE 2. BAF_ck, BAF_total and BAF_add for different potato varieties in Guiyang soil.

TABLE 1

TABLE 1. RMSE of soil Cd threshold calculated with two functions.

TABLE 2

TABLE 2. Soil Cd threshold calculated with two functions for two soil types.

Comparison of Log-Normal and Burr III Functions

SSD is used to create a cumulative probability distribution curve with toxicity data that follow certain probability distributions such as the Log-normal distribution. A cumulative probability of “p” indicates that (100-p)% of organisms in the ecosystem are relatively safe in an environment containing the pollutant at the given concentration (Maltby, et al., 2009). It is worth noting that the “p” value is determined by the local environment safety management department rather than science, and 5% is often used as the threshold value for ecological safety. The soil Cd threshold for each potato variety was obtained based on the BAF (C_{soil threshold} = 0.1/BAF, 0.1 mg kg⁻¹ is the limit for the Chinese food safety standard (GB 2762–2017); BAF can refer to BAF_ck, BAF_total or BAF_add). The C_{soil threshold} of 10 potato varieties was applied to calculated the final soil Cd threshold in order to protect 95% potato and the value of soil Cd threshold was equaled to the 5% cumulative probability. The Burr III or Log-normal function was used to generate an SSD curve and a 5% cumulative probability to protect 95% of species was set as the final soil threshold to protect most potato varieties from Cd damage. A similar method was applied to determine the Cd soil threshold for wheat in different soils using the SSD fitted with the Log-normal function (Liu, et al., 2015). The soil Cd threshold increased with the increase of the cumulative probability. When the cumulative probability was less than 30%, the soil threshold calculated with the Log-normal function was less than that calculated with the Burr III function. In contrast, when the cumulative probability was greater than 30%, the opposite results were obtained for the soil threshold calculated with the two functions. The BAF_ck indicated that the sensitivity levels of the QS9 and HMR varieties to soil Cd were significantly different. QY8 and QS9 easily absorbed Cd from the soil and other potato varieties were relatively insensitive to Cd compared to these two varieties (Figure 3A). The BAF_total values for the potato varieties were evenly distributed on the two SSD curves. The soil Cd threshold fitted with the Log-normal function was less than that obtained with the Burr III function at a low cumulative probability (<30%) and was similar to the BAF_ck. The HMR variety was the most insensitive to Cd, evidenced by both the BAF_ck and BAF_total. The soil Cd threshold was obtained based on the BAF_total and the SSD curves fitted with the Burr III and Log-normal functions were “S-shaped” (Figure 3B). When we focused on the effect of Cd addition on the bioavailability of potato tubers (BAF_add), we found that the CA5e potato variety was relatively sensitive to soil Cd while the HMR variety had a strong resistance to soil Cd uptake. The curves calculated with the Burr III and Log-normal functions were almost the same when the cumulative probability was less than 5%; therefore, the soil Cd thresholds calculated with the two functions were relatively similar when the protection of more than 95% of potato varieties was taken into consideration (Figure 3C). Based on a combination of the BAF_total and BAF_add values, the CA5e and HMR potato varieties were the most sensitive and insensitive to soil Cd in Weining soil in Karst areas, respectively. The identification of Cd-resistant potato varieties is helpful for local food safety. In a study using SSD, Xiangzao 17 was identified as the most sensitive rice cultivar to soil Cd among 20 rice cultivars based on the cumulative probability distribution curve calculated with 1/BAF. The results demonstrated that the SSD method was helpful to screen sensitive and insensitive varieties for pollutants (Song, et al., 2015).

FIGURE 3

FIGURE 3. SSD curves of different potato varieties in Weining soil. (A). BAF_ck, (B). BAF_total, (C). BAF_add.

The soil threshold for Cd established with the Burr III function was higher than that calculated with the Log-normal function at a low cumulative probability, and the opposite was true for a high cumulative probability (Figures 4A,B). This result was similar to the results for the SSD curves of BAF_ck and BAF_total for the Cd treatments in Weining soil. However, the SSD curve fitting for BAF_add with the Burr III function was quite different from that with the Log-normal function, and unlike the SSD curve in acidic Weining soil (Figure 4C). QY8 was sensitive and L15 was insensitive to Cd uptake, as indicated by their respective curves for BAF_total and BAF_add. The Burr III function is considered the best function for soil Cd threshold calculations and has been used to determine the soil Cd threshold for vegetables including potato in acidic and neutral soils based on the BAF_add. Cd-sensitive species and hazardous concentrations within the fifth percentile, reflecting the protection of 95% of species, can also be determined using the SSD (Ding, et al., 2018). The sensitivity levels of the potato varieties were similar in the acidic Weining soil and the alkaline Guiyang soil, as indicated by the BAF_total and BAF_add, but differed from the sensitivity levels with CK treatment, as indicated by BAF_ck. The toxicity sensitivity of different potato varieties to the soil Cd can be directly observed from the SSD curves based on the BAF_ck, BAF_total and BAF_add, which reflect the transfer capacity of heavy metals.

FIGURE 4

FIGURE 4. SSD curves of different potato varieties in Guiyang soil. (A). BAF_ck, (B). BAF_total, (C). BAF_add.

The RMSE represents the deviation between the prediction and the measured value, and a smaller RMSE indicates better prediction accuracy. The RMSE between the predicted and measured soil Cd threshold values for the 10 potato varieties established with BAF_add was larger than that established with BAF_ck and BAF_total in the two soil types using the two functions (Table 1). This suggested that BAF_total and BAF_ck were more suitable for studying the transfer of Cd from the soil to plants. The predictive accuracy of the BAF_ck and BAF_total values was better in Weining and Guiyang soils, respectively. The RMSE values for the Burr III function with the BAF_ck and BAF_total values were less than those for the Log-normal function in acidic Weining soil, indicating that the Burr III function was better suited for determining the soil Cd threshold for this soil type. Log-normal was more appropriate for BAF_add and BAF_ck in alkaline Guiyang soil. The high prediction accuracy of the soil threshold for BAF_ck, BAF_total and BAF_add confirmed that the Burr III and Log-normal functions can both be applied to obtain the soil Cd threshold. However, the optimal function needs to be selected depending on the soil type and BAF calculation method. BAF_add was reported to be better for the determination of the soil Cd threshold for the safe production of rice in paddy soils (pH: 4.94–6.72) and the Burr III function is commonly used to fit SSD curves (Li L et al., 2021). The results of these studies are different from ours, in which BAF_total was better than BAF_add, probably due to different experimental conditions and species.

As shown in Table 2, the soil Cd threshold for the Weining soil calculated with the Log-normal function was a little less than that calculated with the Burr III function using the BAF_ck, BAF_total and BAF_add values. The soil Cd threshold obtained with BAF_add was higher than those obtained with BAF_ck and BAF_total. The average soil Cd threshold was slightly higher than the Chinese national soil standard for Cd (0.3 mg kg⁻¹). Similarly, the soil Cd threshold in Guiyang soil calculated with the Burr III function was higher than that calculated with the Log-normal function and was slightly higher than the Chinese soil screening value for agriculture (0.6 mg kg⁻¹). The soil Cd threshold for BAF_add was lower than those for BAF_ck and BAF_total. These findings differed from the results obtained for Weining soil. Some researchers have applied similar methods to determine the soil threshold for other pollutants. The soil Pb threshold was obtained for 12 root vegetable cultivars based on the SSD and the threshold values based on BAF_add were lower than the standard for Pb in the soil (Ding, et al., 2016). Risk evaluation of Cd in soil-rice systems (n = 124) was conducted using SSD curves from rice Cd-BAF_total in Hunan Province, China (Yang et al., 2016b). The soil Cd threshold obtained in this study is conducive to the protection of local potatoes from Cd pollution, the evaluation of Cd pollution and the development of soil quality guidelines.

Conclusion

The BAF values for different potato varieties in acidic soil were higher than those in alkaline soil for the CK and Cd treatments. The HMR and CA5e potato varieties were quite insensitive and sensitive in acidic soil, whereas QY8 and L15 were sensitive and insensitive in alkaline soil.

The soil Cd threshold values obtained with the Log-normal function using BAF_ck, BAF_total and BAF_add were less than those obtained with the Burr III function in the Weining and Guiyang soils. The soil Cd threshold values were in the order BAF_add > BAF_total > BAF_ck in Weining soil for the Log-normal and Burr III functions. The soil Cd threshold values were in the order BAF_ck > BAF_total > BAF_add for the Guiyang soil. The RMSE for BAF_total was less than that for BAF_add based on the two functions in the Weining and Guiyang soils, indicating that BAF_total was more appropriate for calculating the soil threshold of Cd compared to BAF_add.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

WL: Writing—original draft preparation, Investigation and formal analysis. XH: Investigation, Writing and Editing. JZ: Experimental formal analysis and software analysis. KL: Methodology, Writing—review and editing, Formal analysis.

Funding

This research was supported by the National Natural Science Foundation of China (21866008) and the Science and Technology Planning Project of Guizhou Province, China (Qiankehejichu no (2019)1103, Qiankehehoubuzu no. (2020)3001).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenvs.2022.808362/full#supplementary-material

References

Ashrafzadeh, S., Gaw, S., Genet, R., Glover, C. N., and Leung, D. W. M. (2017). Natural Variation in Correlations between Cadmium and Micronutrients in Potato Tubers. J. Food Compost. Anal. 59, 55–60. doi:10.1016/j.jfca.2017.02.008

CrossRef Full Text | Google Scholar

Augustsson, A. L. M., Uddh-Söderberg, T. E., Hogmalm, K. J., and Filipsson, M. E. M. (2015). Metal Uptake by Homegrown Vegetables - the Relative Importance in Human Health Risk Assessments at Contaminated Sites. Environ. Res. 138, 181–190. doi:10.1016/j.envres.2015.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, M. E. A., Piotto, F. A., Gaziola, S. A., Jacomino, A. P., Jozefczak, M., and Cuypers, A. etal. (2018). New Insights about Cadmium Impacts on Tomato: Plant Acclimation, Nutritional Changes, Fruit Quality and Yield. Food Energy. Secur. 7. doi:10.1002/fes3.131

CrossRef Full Text | Google Scholar

Chen, Z., Zhao, Y., Gu, L., Wang, S., Li, Y., and Dong, F. (2014). Accumulation and Localization of Cadmium in Potato (Solanum tuberosum) under Different Soil Cd Levels. Bull. Environ. Contam. Toxicol. 92, 745–751. doi:10.1007/s00128-014-1252-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, C., Ma, Y., Li, X., Zhang, T., and Wang, X. (2016). Derivation of Soil Thresholds for lead Applying Species Sensitivity Distribution: A Case Study for Root Vegetables. J. Hazard. Mater. 303, 21–27. doi:10.1016/j.jhazmat.2015.10.027

CrossRef Full Text | Google Scholar

Ding, C., Ma, Y., Li, X., Zhang, T., and Wang, X. (2018). Determination and Validation of Soil Thresholds for Cadmium Based on Food Quality Standard and Health Risk Assessment. Sci. Total Environ. 619-620, 700–706. doi:10.1016/j.scitotenv.2017.11.137

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, C., Zhang, T., Wang, X., Zhou, F., Yang, Y., and Huang, G. (2013). Prediction Model for Cadmium Transfer from Soil to Carrot (Daucus Carota L.) and its Application to Derive Soil Thresholds for Food Safety. J. Agric. Food Chem. 61, 10273–10282. doi:10.1021/jf4029859

CrossRef Full Text | Google Scholar

Ding, C., Zhang, T., Wang, X., Zhou, F., Yang, Y., and Yin, Y. (2014). Effects of Soil Type and Genotype on Cadmium Accumulation by Rootstalk Crops: Implications for Phytomanagement. Int. J. Phytoremediation 16, 1018–1030. doi:10.1080/15226514.2013.810581

CrossRef Full Text | Google Scholar

Hough, R. L., Breward, N., Young, S. D., Crout, N. M. J., Tye, A. M., au, A. M., et al. (2004). Assessing Potential Risk of Heavy Metal Exposure from Consumption of home-produced Vegetables by Urban Populations. Environ. Health Perspect. 112, 215–221. doi:10.1289/ehp.5589

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, W., Huang, B., Borggaard, O. K., Ye, M., Tian, K., au, H., et al. (2018). Soil Threshold Values for Cadmium Based on Paired Soil-Vegetable Content Analyses of Greenhouse Vegetable Production Systems in China: Implications for Safe Food Production. Environ. Pollut. 241, 922–929. doi:10.1016/j.envpol.2018.06.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Chen, H., Liu, H., Schlenk, D., Mu, J., au, S., et al. (2021). Anticancer Drugs in the Aquatic Ecosystem: Environmental Occurrence, Ecotoxicological Effect and Risk Assessment. Environ. Int. 153, 106543. doi:10.1016/j.envint.2021.106543

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Cao, C., Ma, Y., Su, D., and Li, J. (2019). Identification of Cadmium Bioaccumulation in rice (Oryza Sativa L.) by the Soil-Plant Transfer Model and Species Sensitivity Distribution. Sci. Total Environ. 692, 1022–1028. doi:10.1016/j.scitotenv.2019.07.091

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Jiang, B., Wan, Y., Li, J., and Ma, Y. (2021). Integrating Bioavailability and Aging in the Criteria Derivation of Cadmium for the Safe Production of rice in Paddy Soils. Ecotox. Environ. Safe. 219. doi:10.1016/j.ecoenv.2021.112356

CrossRef Full Text | Google Scholar

Liang, Z., Ding, Q., Wei, D., Li, J., Chen, S., and Ma, Y. (2013). Major Controlling Factors and Predictions for Cadmium Transfer from the Soil into Spinach Plants. Ecotoxicology Environ. Saf. 93, 180–185. doi:10.1016/j.ecoenv.2013.04.003

CrossRef Full Text | Google Scholar

Liu, K., Lv, J., Dai, Y., Zhang, H., and Cao, Y. (2016). Cross-Species Extrapolation of Models for Predicting lead Transfer from Soil to Wheat Grain. PLoS One 11, e0160552. doi:10.1371/journal.pone.0160552

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, K., Lv, J., He, W., Zhang, H., Cao, Y., and Dai, Y. (2015). Major Factors Influencing Cadmium Uptake from the Soil into Wheat Plants. Ecotoxicology Environ. Saf. 113, 207–213. doi:10.1016/j.ecoenv.2014.12.005

CrossRef Full Text | Google Scholar

Maltby, L., Brock, T. C. M., and van den Brink, P. J. (2009). Fungicide Risk Assessment for Aquatic Ecosystems: Importance of Interspecific Variation, Toxic Mode of Action, and Exposure Regime. Environ. Sci. Technol. 43, 7556–7563. doi:10.1021/es901461c

PubMed Abstract | CrossRef Full Text | Google Scholar

Novotná, M., Mikeš, O., and Komprdová, K. (2015). Development and Comparison of Regression Models for the Uptake of Metals into Various Field Crops. Environ. Pollut. 207, 357–364. doi:10.1016/j.envpol.2015.09.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinto, E., Almeida, A. A., and Ferreira, I. M. P. L. V. O. (2015). Assessment of Metal(loid)s Phytoavailability in Intensive Agricultural Soils by the Application of Single Extractions to Rhizosphere Soil. Ecotoxicology Environ. Saf. 113, 418–424. doi:10.1016/j.ecoenv.2014.12.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, W.-e., Chen, S.-b., Liu, J.-f., Chen, L., Song, N.-n., au, N., et al. (2015). Variation of Cd Concentration in Various rice Cultivars and Derivation of Cadmium Toxicity Thresholds for Paddy Soil by Species-Sensitivity Distribution. J. Integr. Agric. 14, 1845–1854. doi:10.1016/s2095-3119(14)60926-6

CrossRef Full Text | Google Scholar

van der Hoeven, N. (2001). Estimating the 5-percentile of the Species Sensitivity Distributions without Any Assumptions about the Distribution. Ecotoxicology (London, England) 10, 25–34. doi:10.1023/a:1008998405241

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Chen, W. P., Li, Y. L., Wang, M. E., and Peng, C. (2016a). Comprehensive Risk Evaluation of Cadmium in Soil-rice System Based on Uncertainty Analysis. Huan jing ke xue 37, 4800–4805. doi:10.13227/j.hjkx.201606073

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Chen, W., Wang, M., and Peng, C. (2016b). Regional Accumulation Characteristics of Cadmium in Vegetables: Influencing Factors, Transfer Model and Indication of Soil Threshold Content. Environ. Pollut. 219, 1036–1043. doi:10.1016/j.envpol.2016.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, Y., Dong, W., Luo, Y., Fan, T., Xiong, X., au, L., et al. (2020). Cultivar Diversity and Organ Differences of Cadmium Accumulation in Potato (Solanum tuberosum L.) Allow the Potential for Cd-Safe Staple Food Production on Contaminated Soils. Sci. Total Environ. 711, 134534. doi:10.1016/j.scitotenv.2019.134534

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Wu, X., Tu, C., Huang, X., Zhang, J., au, H., et al. (2020). Relationships between Soil Properties and the Accumulation of Heavy Metals in Different Brassica Campestris L. Growth Stages in a Karst Mountainous Area. Ecotoxicol Environ. Saf. 206, 111150. doi:10.1016/j.ecoenv.2020.111150

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, M.-S., Jiang, L., Han, D., Xia, T.-X., Yao, J.-J., au, X.-Y., et al. (2015). Cadmium Exposure via Diet and its Implication on the Derivation of Health-Based Soil Screening Values in China. J. Expo. Sci. Environ. Epidemiol. 25, 433–442. doi:10.1038/jes.2015.5

CrossRef Full Text | Google Scholar

Zhu, H., Chen, C., Xu, C., Zhu, Q., and Huang, D. (2016). Effects of Soil Acidification and Liming on the Phytoavailability of Cadmium in Paddy Soils of central Subtropical China. Environ. Pollut. 219, 99–106. doi:10.1016/j.envpol.2016.10.043

PubMed Abstract | CrossRef Full Text | Google Scholar