Next Article in Journal
Assessment and Application of EPIC in Simulating Upland Rice Productivity, Soil Water, and Nitrogen Dynamics under Different Nitrogen Applications and Planting Windows
Previous Article in Journal
Identification of Soybean Germplasm and Associated Molecular Markers with Resistance to Fusarium graminearum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bioavailability of Cd in Agricultural Soils Evaluated by DGT Measurements and the DIFS Model in Relation to Uptake by Rice and Tea Plants

1
School of Geographical Science, Nantong University, Nantong 226019, China
2
Technology Innovation Center for Ecological Monitoring & Restoration Project on Land (Arable), Ministry of Natural Resources, Geological Survey of Jiangsu Province, Nanjing 210018, China
3
Geological Survey of Anhui Province (Anhui Institute of Geological Sciences), Hefei 230001, China
4
State Environmental Protection Key Laboratory of Soil Environmental Management and Pollution Control, Nanjing Institute of Environmental Sciences, Ministry of Ecology and Environment, Nanjing 210042, China
5
Nanjing Center, China Geological Survey, Nanjing 210016, China
6
State Key Laboratory of Lake Science and Environment, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing 210008, China
*
Author to whom correspondence should be addressed.
Agronomy 2023, 13(9), 2378; https://doi.org/10.3390/agronomy13092378
Submission received: 21 August 2023 / Revised: 5 September 2023 / Accepted: 12 September 2023 / Published: 13 September 2023
(This article belongs to the Section Soil and Plant Nutrition)

Abstract

:
The elevated accumulation of cadmium (Cd) in rice (Oryza sativa L.) and tea (Camellia sinensis L.) grown in agricultural soils may lead to a variety of adverse health effects. This study collected and analyzed crop samples along with paired rhizosphere soil samples from 61 sites in Cd-contaminated regions in Anhui Province, China. The findings revealed that both the diffusive gradients in thin-films (DGT) and soil solution were capable of effectively predicting Cd contents in crops. Conventional chemical extraction methods were inappropriate to evaluate the bioavailability of Cd. However, the effective concentrations (CE) corrected by the DGT-induced fluxes in soils (DIFS) model exhibited the strongest correlation with crop Cd contents. Except for CE, various measurement methods yielded better results for predicting Cd bioavailability in tea compared to rice. Pearson’s correlation analysis and the random forest (RF) model identified the key influencing factors controlling Cd uptake by rice and tea, including pH, soil texture, and contents of zinc (Zn) and selenium (Se) in soils, which antagonize Cd. To reduce the potential health risk from rice and tea, the application of soil liming and/or Se-oxidizing bacteria was expected to be an effective management strategy.

1. Introduction

Cadmium (Cd), which has been identified as the primary contaminant in Chinese agricultural land, poses a substantial threat to both food safety and public health due to its widespread presence in China’s agricultural ecosystems [1]. Elevated Cd accumulation in Chinese soils is the result of a complex interplay between natural geological factors and human activities, including industrialization and agricultural practices. While geological bedrock and mineral deposits can naturally contribute to Cd content in certain areas, anthropogenic activities, especially industrial pollution and agricultural practices, have played a significant role in increasing Cd levels in many regions [2]. In comparison to other toxic elements, Cd exhibits a heightened capacity for absorption by edible plants from soils characterized by high transfer rates. Given its status as a noxious and carcinogenic metal for humans, exposure to Cd can induce kidney damage and various other detrimental impacts, and long-term exposure to Cd has been associated with an increased risk of certain cancers, particularly lung cancer [3]. Some studies have suggested that adequate selenium (Se) intake may mitigate the adverse health effects of Cd exposure. Se can also enhance the activity of certain enzymes that are involved in detoxifying and eliminating heavy metals, including Cd, from the body. This is particularly relevant in regions where there is a risk of Cd contamination in the food supply or environment [4].
Rice (Oryza sativa L.), which serves as a pivotal cereal crop in eastern China, holds a propensity to accumulate more Cd than other cereals, establishing itself as a principal contributor to human Cd exposure through dietary intake. Correspondingly, tea (Camellia sinensis L.), a widely consumed beverage worldwide, particularly in Asia, is susceptible to Cd contamination. In some Asian countries like India, China, and Japan, Cd levels in tea infusions may surpass permissible daily limits, thereby giving rise to cumulative chronic and acute health complications through oral ingestion [5,6]. Moreover, the health risks stemming from Cd exposure via crops hinge not only upon its concentrations within edible portions but also on consumption rates [7]. Consequently, an accurate assessment of Cd bioavailability in rice and tea cultivated in Anhui, a typical province in China where the local populace predominantly relies on rice for caloric intake and tea as a favored hot beverage, becomes an exigent necessity.
Soil characteristics such as pH, organic matter content, and soil texture can exert influence over the bioavailability of Cd [8,9,10]. For example, in acidic conditions, Cd ions are less likely to bind to soil particles, making them more mobile and potentially more harmful to plants. The dissolution of Cd-containing minerals can also increase under acidic conditions; for example, Cd sorption on kaolinite reduced significantly at low soil pHs (pH < 6) since Cd is sorbed as outer-sphere complexes. Organic matter can form complexes with Cd ions, reducing their mobility and uptake by plants. This is due to the formation of stable complexes between Cd and organic ligands. Further, Cd tends to have higher retention and lower availability in clayey soils due to the higher cation exchange capacity (CEC) and surface area for adsorption. Cd ions can bind tightly to clay particles, reducing their movement in the soil [11,12]. However, the specific mechanisms by which these influencing factors affect the bioavailability of Cd remain insufficiently elucidated. The classification and regression tree (CART) methodology stands as a widely employed machine learning approach known for its adeptness in unveiling relationships. In tandem, the random forest (RF) technique, an extension of CART, represents a relatively novel method aimed at enhancing prediction accuracy [13,14,15].
Furthermore, it is widely acknowledged that the uptake of Cd by plants is contingent not solely upon Cd speciation within soils but also upon the dynamic interplay of Cd resupply from the solid phase to the soil solution [16]. To gauge labile Cd concentrations, encompassing both dissolved fractions within the soil solution and a portion of resupplied fractions from the solid phase, the diffusive gradients in thin films (DGT) technique is employed. Previous investigations have attested to the efficacy of the DGT technique in evaluating Cd bioavailability in soils [17,18]. Facilitating quantitative descriptions of chemical exchange between the solid phase and solution under DGT-induced perturbation, the DGT-induced fluxes in soils (DIFS) model provides a dynamic framework for the DGT-medium system. Additionally, this model permits the calculation of diffusion coefficients to rectify DGT measurements [19,20].
The evaluation of Cd bioavailability in rice and tea, along with its governing factors within the soil properties, is imperative within this study area. This study encompasses the pursuit of several objectives: (i) to appraise the predictive efficacy of Cd accumulation in rice grains and tea leaves, employing diverse measurement approaches encompassing conventional chemical extraction, the DGT technique, soil solution analysis, and the DIFS model; (ii) to discern the factors that wield influence over Cd uptake by rice and tea; (iii) to offer recommendations that can steer secure local cultivation management.

2. Materials and Methods

2.1. Study Area

The research area is located in Shitai County, Anhui Province, China. The area is about 24 km2, with elevation between 50 to 1000 m, and the terrain is generally mountains in the northwest and valley plains in the middle. Limestone, shale, and carbonaceous siliceous rock are the dominant lithology across the study area. The annual precipitation is approximately 1100 mm, and the daily temperature ranges between 5 and 25 °C [10]. The main types of soil-forming parent materials are carbonate rocks, clastic rocks, and alluvial deposits. The major soil types have been identified as Cambisols, Acrisols, and Anthrosols [21]. The soils in Shitai County have been identified as seleniferous and Cd-contaminated [10,22]. Rice and tea are the most extensively grown crops in the region.

2.2. Sample Collection and Pretreatment

In this investigation, a total of 61 agricultural sites were surveyed. We collected 43 samples of rice grains (Indica cultivar) and 18 samples of tea leaves, with principal growth stage of ripening, and their corresponding rhizosphere soil samples were extracted from depths of 0 to 20 cm (Figure 1). The soil samples were collected from 3 to 5 sub-sample points per site and subsequently amalgamated to form a composite representative of each location. Employing a meticulous approach, all collected samples were hermetically sealed within airtight plastic pouches to maintain their integrity during transportation to the laboratory within a time of 48 h. After collection, the soil samples were air-dried at ambient temperature and subsequently sieved through a 10-mesh sieve to eliminate any extraneous matter or pebbles, after which they were finely ground into a powder of uniform consistency. The rice grains, divested of their husks, and the tea leaves were subjected to an oven-drying process at 50 °C for a duration of three days to ensure a consistent weight, following which they too were ground into a finely textured powder.

2.3. Chemical Analysis

Soil pH was measured using the ion selectivity electrode (ISE) method. Soil particle size composition was determined using Malvern Mastersizer 2000, with sand, silt, and clay fractions defined as particles >20 μm, 2~20 μm, and <2 μm, respectively [23]. The oxidation with dichromate procedure is used for the determination of soil organic carbon (SOC). The measurement of soil cation exchange capacities (CEC) was carried out through the ammonium acetate exchange method, as detailed in reference [24]. The determination of total Cd concentrations in soil, rice grains, and tea leaves involved the use of 1.00 g of powdered samples for each matrix. Soil samples underwent digestion with HF, HNO3, and HClO4, while the crop samples were subjected to digestion with a concentrated mixture of HNO3 and HClO4. Inductively coupled plasma mass spectrometry (ICP-MS) was employed to quantify Cd, Zn, and Se concentrations in both soil and crop digestion solutions. Quality assurance and quality control (QA/QC) measures were applied to ensure accuracy and precision, incorporating blank samples, duplicates, and standard reference materials.
The sequential extraction procedure, adapted from the methodologies of Tessier et al. [25] and Quezada-Hinojosa et al. [26], facilitated the fractionation of Cd into seven distinct forms: F1 water-soluble, F2 exchangeable (extracted for 1 h using 1 mol·L−1 MgCl2), F3 carbonate-bound, F4 humic acid-bound, F5 amorphous Fe/Mn oxide-bound, F6 crystalline Fe oxide-bound, and F7 residue [27]. The filtrates obtained at each extraction step were subjected to ICP-MS analysis. The reliability of the sequential extractions was verified by comparing the cumulative concentrations of fractions F1 to F7 with the total Cd concentrations.

2.4. DGT Measurements and the DIFS Model

In this investigation, we employed Chelex-100 resin DGT devices procured from Weishen DGT® Research Ltd. (Nanjing, China). Our procedures adhered to established standards for the deployment and measurement of DGT devices in soil, as outlined in reference [28]. A total of sixty-one rhizosphere soil samples were subjected to measurement via the DGT technique. Preceding device deployment, 40 g of air-dried soil was introduced into a beaker, set at 60% of its maximum water holding capacity (MWHC), for a span of 48 h, followed by a subsequent phase at 80% MWHC for 24 h. After 24 h of deployment at a temperature of 25 °C, the DGT devices were carefully retrieved from the soil, subjected to a thorough rinse with deionized water, and then disassembled. Subsequently, the binding gels were subjected to elution with 1 mL of 1 M HNO3, allowing for a 24-h interval prior to the commencement of analysis. The eluates were diluted by a factor of ten utilizing HNO3 before Cd analysis was conducted using ICP-MS. The computed concentration of Cd (CDGT) through the DGT approach was determined using Equation (1), with a comprehensive elucidation of the calculation provided in reference [20]:
C D G T = M × Δ g D × A × t
In this equation, the symbol M represents the accumulated mass of Cd (in ng) during the deployment interval, while Δg denotes the thickness of the diffusive layer (in cm). The diffusion coefficient of Cd within the resin is represented by D (in cm2·s−1), A signifies the area of the exposure window (in cm2), and t corresponds to the deployment time (in seconds). Following the retrieval of the DGT devices, soil paste was procured, transferred into 50 mL polyethylene tubes, and then subjected to centrifugation at 5000 rpm for a duration of 20 min, aimed at extracting the soil solution. After centrifugation, the supernatant was subjected to filtration through a 0.45 μm filter, with the resultant filtrate utilized for the quantification of the soil solution’s Cd content (Csolu). To assess the resupply originating from the solid phase, the ratio (R) between CDGT and Csolu was computed [29]. To comprehend the kinetics of Cd resupply from the soil’s solid phase to the porewater within the rhizosphere zone or to ascertain the effective concentration (CE) of Cd, the DGT-induced fluxes in soils (DIFS) model for the DGT-soil system was employed. The term CE pertains to the theoretical porewater concentration necessary for the accumulation of the measured Cd quantity on the DGT resin solely via diffusional supply [18]. Equation (2) is employed to convert CDGT to CE:
C E = C D G T R d i f f
Here, Rdiff denotes the CDGT to Csolu ratio when supply to the DGT device solely occurs via diffusion. The DIFS model permits the quantification of the interdependency between Rdiff and the rate and extent of Cd resupply from the solid phase to the solution, encompassing both the soil-to-interface transition and diffusion through the diffusion layer to the binding layer’s sink [20]. The software 2D DIFS (version 1.2.3-3, 2005, Lancaster, UK) [30], available online, was utilized for the calculation of Rdiff values pertinent to the Cd in the DGT-soil system.

2.5. Statistical Analysis

Descriptive statistics and Pearson’s correlation analysis were conducted employing IBM SPSS Statistics version 27 (IBM Corp, Armonk, NY, USA). To ascertain and quantify the significance of potential factors influencing the bioavailability of Cd in agricultural soils, the random forest (RF) technique was employed. Implementation of the RF model was carried out in R version 4.2.3 (R Development Core Team, Vienna, Austria), utilizing the random forest package [31]. For graphical representation, box plots and linear fit analyses were generated using Origin 2021 (Origin Lab Corporation, Northampton, MA, USA).

3. Results and Discussion

3.1. Concentrations of Cadmium in Soils, Rice Grains and Tea Leaves

The box plots illustrating the total Cd content in the plow layer soil, rice grains, and tea leaves within the study area are presented in Figure 2. In these box plots, the internal horizontal line represents the median value of Cd, while the cross symbol denotes the mean value. The top and bottom edges of the box correspond to the 25th and 75th percentiles of Cd, and the lines extending from the top and bottom of the box represent the 5th and 95th percentiles of Cd, respectively. In the Cd-contaminated region of Anhui, the variation in Cd total content in rhizosphere soils ranged from 0.117 to 1.029 mg·kg−1, with an average of 0.564 mg·kg−1 in rice grains. The average Cd content was 0.191 mg·kg−1 in rice seeds and 0.038 mg·kg−1 in tea leaves within the study area (Figure 2b).
We employed the modified Tessier sequential extraction method to investigate the various binding forms of Cd in 61 rhizosphere soil samples [25,26], yielding a total of seven forms (F1~F7). Among these forms, the exchangeable Cd fraction exhibited the highest content, constituting 42.4% of the total, indicating notable Cd mobility (Figure 2a). The amorphous Fe/Mn oxide-bound fraction also exhibited considerable content at 21.9% of the total, signifying the influence of Fe/Mn oxides on Cd distribution. The residual form appeared as the most stable, accounting for 16.0% of the total Cd content. The carbonate-bound and humic acid-bound fractions demonstrated similar contents, representing 7.19% and 7.37% of the total Cd content, respectively. However, the carbonate-bound fraction displayed greater variability among different soil samples. The crystalline Fe oxide-bound and water-soluble fractions exhibited the lowest contents at 3.90% and 1.21% of the total Cd content, respectively. In summary, Cd exhibited a substantial proportion of mobile forms. Previous research suggests that Cd in soil is more prone to migration compared to other heavy metals, thereby necessitating heightened attention to its ecological risk [2].
The bioconcentration factor (BCF) is defined as the ratio of the concentration of heavy metals, such as Cd, in crops to the concentration of these elements in the topsoil where the crops are grown. It is commonly used to indicate the capacity of crops to absorb and accumulate heavy metals. The formula for calculating the BCF is provided in Equation (3) [32]:
BCF   = C crop C soil
where BCF represents the value of the bioconcentration factor, Ccrop is the concentration of the heavy metal in the crops, and Csoil is the corresponding concentration of the heavy metal in the rhizosphere soil. In this study, BCFrice and BCFtea represent the BCF for Cd in rice and tea, respectively. The transfer of harmful elements like Cd from the soil to crops is a vital link in the food chain for human exposure to heavy metal pollutants. Investigating the bioaccumulation characteristics of Cd is essential for evaluating the safety risks associated with agricultural products and managing safe crop cultivation. Table 1 presents the Cd content in soil and crops, BCF values, and exceedance rates for the study area as well as other Cd-contaminated regions (Jiangsu and Guizhou) and non-contaminated regions (Henan and Yunnan).
Recent studies have shown that in Cd-contaminated regions in Jiangsu, the average BCFrice is approximately 0.155 [8]. In comparison, in the agricultural fields of Henan, the average BCFrice is around 0.145, which is at a similar level [33]. In the study area, the average BCFrice is higher than those in both the contaminated and non-contaminated regions. This indicates a greater Cd activity in the Cd-contaminated areas of Anhui. As discussed earlier, the sum of the water-soluble, exchangeable, and carbonate-bound forms (F1+F2+F3) exceeds 50% of the total Cd content in the soil, highlighting its potential biological availability [34,35]. However, according to the national standards [36,37], the exceedance rate of rice (32.5%) is at a comparable level to that of Jiangsu (30.6%), whereas it is more than 10 times higher than the non-contaminated area (2.5%). In contrast to rice, the ability of tea leaves to accumulate Cd shows notable differences between contaminated and non-contaminated areas. In Guizhou, the average BCFtea (0.303) is ten times higher than that in Yunnan (0.038), while in this study, BCFtea falls between these two values [38,39]. The underlying cause of this divergence in rice and tea lies potentially in the significantly lower pH of the rhizosphere soil for tea plants compared to rice, resulting in a greater release of labile forms for Cd into the soil solution, thereby facilitating their uptake by crops.
Table 1. The median contents of total soil Cd, crop Cd, BCF, and exceedance rate according to national thresholds.
Table 1. The median contents of total soil Cd, crop Cd, BCF, and exceedance rate according to national thresholds.
RegionAnhuiAnhuiJiangsuGuizhouHenanYunnan
Cropricetearicetearicetea
n43187022408
pH6.335.396.104.305.904.48
Total soil Cd0.7370.4170.9700.3300.1240.079
Crop Cd0.1910.0380.1500.1000.0180.003
BCF0.2590.0910.1550.3030.1450.038
Thresholds0.2 a1.0 b0.200 a1 b0.200 a1 b
Exceedance rate (%)32.50.030.60.02.50.0
Referencethis studythis study[8][38][33][39]
a National Standards for Food Safety and Limits of Contaminants for Food [36]. b Residue Limits for Chromium, Cadmium, Mercury, Arsenic, and Fluoride in Tea [37].

3.2. Evaluation of Cd Bioavailability

In this study, conventional chemical extraction methods were inappropriate to evaluate the bioavailability of Cd in the soils. The MgCl2 extractable fraction (F1 + F2 from sequential extraction) showed no significant correlation with Cd content in rice grains and tea leaves (Table 2). In this study, the DGT-measured Cd concentration (CDGT) exhibited a positive correlation with Cd content in rice grains (r = 0.61) and tea leaves (r = 0.66), though the predictive performance was superior for tea compared to rice (Table 2). Surprisingly, a strong positive correlation was found between Cd in soil solution (Csolu) and Cd content in crops, with Pearson’s correlation coefficients surpassing those of DGT. However, the predictive performance for rice remained inferior to that for tea (Figure 3a). This phenomenon is primarily attributed to the formation of iron plaque on the root surfaces of many aquatic plants, including rice, within the aquatic rhizosphere microenvironment. This iron plaque is an adaptation mechanism to aquatic environments. Elevated levels of soil iron oxides lead to an increase in the formation of the iron plaque on rice roots. This root-bound iron plaque is predominantly composed of amorphous or crystalline Fe (hydroxy) oxides and tends to immobilize a portion of Cd as oxide-bound forms on rice roots, reducing the Cd uptake by rice plants. Consequently, the DGT measurement of the bioavailable fraction of Cd performs less effectively for rice compared to upland crops like tea [40,41].
The supply capacity of Cd from the solid phase to the soil solution can be quantified using the R value (R = CDGT/Csolu, 0 < R < 1), where a higher R value (approaching 1) indicates a stronger ability of soil particles to resupply Cd into the soil solution [28]. When the supply of Cd from the solid phase to the solution is limited or approaches negligible supply (mainly diffusion), the R value approaches its minimum possible value (close to 0). In this study, the sequence of different soil R values from high to low is as follows: soils with pH 4.5~5.5 > soils with pH 5.5~6.5 > soils with pH above 6.5. The average R value for soils with pH above 6.5 is 0.39, notably lower than other soils, indicating a continuous supply process of solid-phase Cd into the soil solution in slightly alkaline soils. However, for other acidic soils, especially those with pH 4.5~5.5, the R value is higher (average 0.62), suggesting a relatively abundant supply of solid-phase Cd to the soil solution in more acidic conditions. In these more acidic soils, as plants absorb Cd from the soil solution, the solid phase of the soil can release a greater amount of bioavailable Cd into the liquid phase.
In order to investigate the concentrations of analytes that are effectively utilizable in both the soil solid phase and soil solution, the concept of effective concentration (CE) is commonly utilized when employing DGT for studying soil environments. The calculation method for CE is outlined in Equation (2). The diffusion coefficient (Rdiff) values calculated using the DIFS model are used to correct CDGT, resulting in the determination of analyte concentrations that can be effectively utilized within the soil. When the correlation between CDGT and Cd content in crops is suboptimal, employing CE for predicting Cd accumulation in crops often yields improved results [18]. In this study, the predictive performance of CDGT for Cd accumulation in crops was not entirely satisfactory. Consequently, a correction was applied using the CDGT and Rdiff to derive the effective concentration (CE) for a refined prediction. This approach demonstrated enhanced predictive accuracy (Figure 3b), especially the correlation coefficient (r = 0.84) between CE and Cd accumulation in rice, surpassing the correlation coefficient (r = 0.79) between CE and Cd accumulation in tea leaves. This observation underscores the capacity of the combined DGT technique and DIFS model to better simulate the dynamic migration of Cd in the soil-crop system.

3.3. Factors Controlling Cd Uptake by Rice and Tea

To identify the major factors influencing the uptake of Cd by rice and tea from the soil, we conducted a correlation analysis between BCFrice, BCFtea, and various soil physicochemical parameters. The Pearson’s correlation coefficients between BCF values and soil properties for rice and tea samples are presented in Table 2. Overall, significant correlations were observed between BCFrice and soil properties such as pH, Se, Zn, clay, and sand content, while BCFtea was influenced by soil pH, Se, sand, and clay content. Given the multitude of factors affecting the migration and accumulation of Cd in soil-crop systems, we introduced the random forest (RF) method as an advanced classification approach to further verify the influencing factors on Cd absorption by rice and tea. The RF model, an extension of the classification and regression trees (CART) model, is described in detail by Breiman [13]. Its advantages include sensitivity to noise and reduced bias in error estimation compared to other models [10]. Notably, soil properties including pH, Zn, sand, Se, and clay content account for 23.6%, 19.6%, 14.7%, 13.5%, and 11.4%, respectively, in influencing BCFrice (Figure 4a). Similarly, sand, pH, clay, and Se account for 22.6%, 20.9%, 19.5%, and 16.5%, respectively, in affecting BCFtea (Figure 4b).
Soil pH has been identified as a crucial factor influencing the uptake of Cd by rice and tea, displaying significant negative correlations with both BCFrice and BCFtea (Table 2). Soil pH could affect the speciation of Cd in soil. As soil pH increases, the reduction in H+ leads to a heightened negative surface charge on the soil, thereby promoting the adsorption of Cd2+ by various minerals in the soil. Additionally, pH affects the reactivity of Cd, compound solubility, and the extent of hydroxylation of Cd2+ on mineral surfaces [42,43]. Soil acidification, on the other hand, reduces the adsorption capacity and affinity of minerals for Cd2+ in the soil, converting Cd2+ adsorption from specific to non-specific ion exchange adsorption. Therefore, the application of soil liming can reduce the bioavailability of Cd in the soils, making it less prone to uptake by crops like rice and tea.
Soil texture is commonly recognized as a significant factor influencing the uptake of heavy metals such as Cd by crops. Typically, finer soil particles tend to accumulate more Cd compared to coarser soil particles [44]. This outcome can be attributed to the enhanced water retention, nutrient-holding capacity, and support for above- or below-ground organisms provided by finer soil particles. Moreover, finer particles exhibit greater capability in adsorbing key carriers of Cd, such as organic matter and soil minerals. In this study area, there is a notable negative correlation between sand content and both BCFrice and BCFtea, while clay content exhibits a significant positive correlation, and the influence of silt content is marginal (Table 2). The impact of sand content on Cd uptake by crops, particularly tea (Figure 4b), is substantial. In soils with high sand content, the ability of rice and tea to absorb Cd is diminished due to their poorer water-holding capacity and weaker adsorption capabilities towards Cd-bearing minerals like iron and manganese oxides [45].
Soil Zn content is an additional factor affecting the uptake of Cd by rice and tea. The soil-crop system constitutes a complex arrangement, wherein alterations in the content of one element can influence the functionality of others, ultimately impacting the migration of heavy metals like Cd within the soil-crop system. There exists a certain antagonistic relationship between Zn and Cd since previous studies have found that applying low doses of ZnO nanoparticles to soil can reduce the toxicity of Cd on crop growth [46]. Given the structural similarity between Zn and Cd ions, higher Zn content in the soil can mitigate Cd toxicity to crops and promote crop growth. As the bioavailable concentration of Zn increases, the bioavailable concentration of Cd in the soil decreases significantly [47].
Furthermore, prior research indicates an antagonistic interaction between Se and Cd. Elevated Se levels in soils can facilitate the deposition of Cd in root cell walls or vacuoles, consequently reducing its transportation to other plant organs, thus diminishing the bioavailability of Cd [48,49]. Soil Se content ranked fourth among the factors influencing the uptake of Cd by both rice and tea (Figure 4). Since Se-oxidizing bacteria can enhance the bioavailability of Se in soils [50]. In this study, the application of Se-oxidizing bacteria in the seleniferous soil was thought to be an effective approach to diminish the uptake of Cd by crops such as rice and tea and its subsequent transport to the aboveground parts.

4. Conclusions

We conducted an analysis of the Cd content characteristics in soil, rice grains, and tea leaves within the study area. Our findings indicate that the cumulative content of water-soluble (F1), exchangeable (F2), and carbonate-bound (F3) Cd fractions collectively accounts for more than 50% of the total Cd content in the soil. This suggests a high potential for Cd bioavailability in the soil. However, this situation is also associated with an elevated rate of Cd accumulation in rice, surpassing 32.5%. When assessing Cd bioavailability using the DGT method, we observed favorable and positive correlations between CDGT, Csolu, and Cd content in crops. The most accurate predictions were obtained by utilizing the corrected CE values based on the Rdiff calculated using the DIFS model.
More alkaline and sandy agricultural soils are considered more suitable for growing rice and tea plants with low Cd content. In addition to soil pH and texture, Se and Zn levels in the soil play crucial roles in influencing the uptake of Cd by rice and tea crops. Despite the known antagonistic effect of Se on Cd uptake by crops, the co-enrichment of Se and Cd in the soil often results in the concurrent accumulation of both elements in crops. This situation may lead to the production of Se-enriched crops that also exceed permissible Cd levels. Therefore, enhancing the bioavailability of Se in the soil and reducing the bioavailability of Cd through soil liming and/or the application of Se-oxidizing bacteria can facilitate the safe cultivation of locally Se-enriched agricultural products while mitigating the risks of Cd contamination.

Author Contributions

Y.W. (Yubo Wen): conceptualization, investigation, visualization, writing—original draft, writing—review and editing; Y.W. (Yuanyuan Wang): conceptualization, investigation, writing—original draft, writing—review and editing; C.T.: conceptualization, investigation, funding acquisition; W.J.: conceptualization, investigation; S.H. and M.Z.: methodology, investigation; X.M.: methodology, investigation, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Natural Resource Technology Project of Anhui Province (no. 2020-K-7), National Natural Science Foundation of China (no. 42207236), Fundamental Science Project of Nantong (no. JC12022075), Special Fund for the Development of Natural Resources in Jiangsu Province “Territory Spatial Ecological Monitoring in Jiangsu Province”, and Fund from the Land (Arable Land) Ecological Monitoring and Restoration Engineering Technology Innovation Center of the Ministry of Natural Resources (No. GTST2021-002).

Data Availability Statement

The data is unavailable due to confidentiality agreement.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have influenced the work reported in this paper.

References

  1. Wang, P.; Chen, H.; Kopittke, P.M.; Zhao, F.-J. Cadmium contamination in agricultural soils of China and the impact on food safety. Environ. Pollut. 2019, 249, 1038–1048. [Google Scholar] [CrossRef]
  2. Zhao, F.-J.; Ma, Y.; Zhu, Y.-G.; Tang, Z.; McGrath, S.P. Soil Contamination in China: Current Status and Mitigation Strategies. Environ. Sci. Technol. 2015, 49, 750–759. [Google Scholar] [CrossRef]
  3. Mao, C.; Song, Y.; Chen, L.; Ji, J.; Li, J.; Yuan, X.; Yang, Z.; Ayoko, G.A.; Frost, R.L.; Theiss, F. Human health risks of heavy metals in paddy rice based on transfer characteristics of heavy metals from soil to rice. Catena 2018, 175, 339–348. [Google Scholar] [CrossRef]
  4. Natasha; Shahid, M.; Niazi, N.K.; Khalid, S.; Murtaza, B.; Bibi, I.; Rashid, M.I. A critical review of selenium biogeochemical behavior in soil-plant system with an inference to human health. Environ. Pollut. 2018, 234, 915–934. [Google Scholar] [CrossRef]
  5. Koupaie, E.H.; Eskicioglu, C. Health risk assessment of heavy metals through the consumption of food crops fertilized by biosolids: A probabilistic-based analysis. J. Hazard. Mater. 2015, 300, 855–865. [Google Scholar] [CrossRef]
  6. Karak, T.; Bhagat, R. Trace elements in tea leaves, made tea and tea infusion: A review. Food Res. Int. 2010, 43, 2234–2252. [Google Scholar] [CrossRef]
  7. Meharg, A.A.; Norton, G.; Deacon, C.; Williams, P.; Adomako, E.E.; Price, A.; Zhu, Y.; Li, G.; Zhao, F.-J.; McGrath, S.; et al. Variation in Rice Cadmium Related to Human Exposure. Environ. Sci. Technol. 2013, 47, 5613–5618. [Google Scholar] [CrossRef]
  8. Li, H.; Luo, N.; Li, Y.W.; Cai, Q.Y.; Mo, C.H.; Wong, M.H. Cadmium in rice: Transport mechanisms, influencing factors, and minimizing measures. Environ. Pollut. 2017, 224, 622–630. [Google Scholar] [CrossRef]
  9. Sarwar, N.; Saifullah Malhi, S.S.; Zia, M.H.; Naeem, A.; Bibi, S.; Farid, G. Role of mineral nutrition in minimizing cadmium accumulation by plants. J. Sci. Food Agric. 2010, 90, 925–937. [Google Scholar] [CrossRef]
  10. Tao, C.; Song, Y.; Chen, Z.; Zhao, W.; Ji, J.; Shen, N.; Frost, R.L. Geological load and health risk of heavy metals uptake by tea from soil: What are the significant influencing factors? Catena 2021, 204, 105419. [Google Scholar] [CrossRef]
  11. Wang, C.; Ji, J.; Yang, Z.; Chen, L.; Browne, P.; Yu, R. Effects of Soil Properties on the Transfer of Cadmium from Soil to Wheat in the Yangtze River Delta Region, China—A Typical Industry–Agriculture Transition Area. Biol. Trace Element Res. 2012, 148, 264–274. [Google Scholar] [CrossRef]
  12. Wu, L.; Zhou, J.; Zhou, T.; Li, Z.; Jiang, J.; Zhu, D.; Hou, J.; Wang, Z.; Luo, Y.; Christie, P. Estimating cadmium availability to the hyperaccumulator Sedum plumbizincicola in a wide range of soil types using a piecewise function. Sci. Total. Environ. 2018, 637–638, 1342–1350. [Google Scholar] [CrossRef]
  13. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  14. Zhang, H.; Yin, A.; Yang, X.; Wu, P.; Fan, M.; Wu, J.; Zhang, M.; Gao, C. Changes in surface soil organic/inorganic carbon concentrations and their driving forces in reclaimed coastal tidal flats. Geoderma 2019, 352, 150–159. [Google Scholar] [CrossRef]
  15. Wen, Y.; Wang, Y.; Ji, W.; Wei, N.; Liao, Q.; Huang, D.; Meng, X.; Song, Y. Influencing Factors of Elevated Levels of Potentially Toxic Elements in Agricultural Soils from Typical Karst Regions of China. Agronomy 2023, 13, 2230. [Google Scholar] [CrossRef]
  16. Luo, J.; Zhang, H.; Zhao, F.-J.; Davison, W. Distinguishing Diffusional and Plant Control of Cd and Ni Uptake by Hyperaccumulator and Nonhyperaccumulator Plants. Environ. Sci. Technol. 2010, 44, 6636–6641. [Google Scholar] [CrossRef]
  17. Wen, Y.; Li, W.; Yang, Z.; Zhuo, X.; Guan, D.-X.; Song, Y.; Guo, C.; Ji, J. Evaluation of various approaches to predict cadmium bioavailability to rice grown in soils with high geochemical background in the karst region, Southwestern China. Environ. Pollut. 2019, 258, 113645. [Google Scholar] [CrossRef]
  18. Williams, P.N.; Zhang, H.; Davison, W.; Zhao, S.; Lu, Y.; Dong, F.; Zhang, L.; Pan, Q. Evaluation of in Situ DGT Measurements for Predicting the Concentration of Cd in Chinese Field-Cultivated Rice: Impact of Soil Cd:Zn Ratios. Environ. Sci. Technol. 2012, 46, 8009–8016. [Google Scholar] [CrossRef]
  19. Almås, R.; Lombnæs, P.; Sogn, T.A.; Mulder, J. Speciation of Cd and Zn in contaminated soils assessed by DGT-DIFS, and WHAM/Model VI in relation to uptake by spinach and ryegrass. Chemosphere 2006, 62, 1647–1655. [Google Scholar] [CrossRef]
  20. Guan, D.-X.; Zheng, J.-L.; Luo, J.; Zhang, H.; Davison, W.; Ma, L.Q. A diffusive gradients in thin-films technique for the assessment of bisphenols desorption from soils. J. Hazard. Mater. 2017, 331, 321–328. [Google Scholar] [CrossRef]
  21. FAO. Food and Agriculture Organization of the United Nations, World reference base for soil resources. In International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; FAO: Rome, Italy, 2015. [Google Scholar]
  22. Long, Z.; Yuan, L.; Hou, Y.; Bañuelos, G.S.; Liu, Y.; Pan, L.; Liu, X.; Yin, X. Spatial variations in soil selenium and residential dietary selenium intake in a selenium-rich county, Shitai, Anhui, China. J. Trace Elements Med. Biol. 2018, 50, 111–116. [Google Scholar] [CrossRef] [PubMed]
  23. Callesen, I.; Keck, H.; Andersen, T.J. Particle size distribution in soils and marine sediments by laser diffraction using Malvern Mastersizer 2000—Method uncertainty including the effect of hydrogen peroxide pretreatment. J. Soils Sediments 2018, 18, 2500–2510. [Google Scholar] [CrossRef]
  24. Tian, Y.; Wang, X.; Luo, J.; Yu, H.; Zhang, H. Evaluation of Holistic Approaches to Predicting the Concentrations of Metals in Field-Cultivated Rice. Environ. Sci. Technol. 2008, 42, 7649–7654. [Google Scholar] [CrossRef]
  25. Tessier, A.; Campbell, P.G.C.; Bisson, M. Sequential extraction procedure for the speciation of particulate trace metals. Anal. Chem. 1979, 51, 844–851. [Google Scholar] [CrossRef]
  26. Quezada-Hinojosa, R.P.; Föllmi, K.B.; Verrecchia, E.; Adatte, T.; Matera, V. Speciation and multivariable analyses of geogenic cadmium in soils at Le Gurnigel, Swiss Jura Mountains. Catena 2015, 125, 10–32. [Google Scholar] [CrossRef]
  27. Wen, Y.; Li, W.; Yang, Z.; Zhang, Q.; Ji, J. Enrichment and source identification of Cd and other heavy metals in soils with high geochemical background in the karst region, Southwestern China. Chemosphere 2019, 245, 125620. [Google Scholar] [CrossRef] [PubMed]
  28. Luo, J.; Cheng, H.; Ren, J.; Davison, W.; Zhang, H. Mechanistic Insights from DGT and Soil Solution Measurements on the Uptake of Ni and Cd by Radish. Environ. Sci. Technol. 2014, 48, 7305–7313. [Google Scholar] [CrossRef]
  29. Tian, K.; Xing, Z.; Liu, G.; Wang, H.; Jia, M.; Hu, W.; Huang, B. Cadmium phytoavailability under greenhouse vegetable production system measured by diffusive gradients in thin films (DGT) and its implications for the soil threshold. Environ. Pollut. 2018, 241, 412–421. [Google Scholar] [CrossRef]
  30. Sochaczewski, Ł.; Tych, W.; Davison, B.; Zhang, H. 2D DGT induced fluxes in sediments and soils (2D DIFS). Environ. Model. Softw. 2007, 22, 14–23. [Google Scholar] [CrossRef]
  31. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013; ISBN 3-900051-07-0. Available online: http://www.R-project.org/ (accessed on 15 March 2023).
  32. Chopra, A.K.; Pathak, C. Accumulation of heavy metals in the vegetables grown in wastewater irrigated areas of Dehradun, India with reference to human health risk. Environ. Monit. Assess. 2015, 187, 1–8. [Google Scholar] [CrossRef]
  33. Liu, W.X.; Shen, L.F.; Liu, J.W.; Wang, Y.W.; Li, S.R. Uptake of toxic heavy metals by rice (Oryza sativa L.) cul-tivated in the agricultural soil near Zhengzhou City, People’s Republic of China. Bull. Environ. Contam. Toxicol. 2007, 79, 209–213. [Google Scholar] [CrossRef] [PubMed]
  34. Ji, W.; Luo, Z.; Huang, J.; Liu, X.; He, H.; Gong, Y.; Chen, M.; Wen, Y.; Ying, R. Geochemical Behaviors of Heavy Metal(loid)s in Soil Ferromanganese Nodules in Typical Karst Areas in Southwest China. Agronomy 2023, 13, 1602. [Google Scholar] [CrossRef]
  35. Agrelli, D.; Caporale, A.G.; Adamo, P. Assessment of the Bioavailability and Speciation of Heavy Metal(loid)s and Hydrocarbons for Risk-Based Soil Remediation. Agronomy 2020, 10, 1440. [Google Scholar] [CrossRef]
  36. GB2762-2017; National Standards for Food Safety and Limits of Contaminants for Food. NHFPC and NMPA, National Health and Family Planning Commission and National Medical Products Administration of the People’s Republic of China: Beijing, China, 2017.
  37. NY 659-2003; Residue Limits for Chromium, Cadmium, Mercury, Arsenic and Fluoride in Tea. MARA, Ministry of Agriculture and Rural Affairs of the People’s Republic of China: Beijing, China, 2003.
  38. Zhang, J.; Yang, R.; Li, Y.C.; Peng, Y.; Wen, X.; Ni, X. Distribution, accumulation, and potential risks of heavy metals in soil and tea leaves from geologically different plantations. Ecotoxicol. Environ. Saf. 2020, 195, 110475. [Google Scholar] [CrossRef]
  39. Cao, H.; Qiao, L.; Zhang, H.; Chen, J. Exposure and risk assessment for aluminum and heavy metals in Puerh tea. Sci. Total Environ. 2010, 408, 2777–2784. [Google Scholar] [CrossRef] [PubMed]
  40. Khan, N.; Seshadri, B.; Bolan, N.; Saint, C.; Kirkham, M.; Chowdhury, S.; Yamaguchi, N.; Lee, D.; Li, G.; Kunhikrishnan, A.; et al. Root Iron Plaque on Wetland Plants as a Dynamic Pool of Nutrients and Contaminants. Adv. Agron. 2016, 138, 1–96. [Google Scholar] [CrossRef]
  41. Chen, Z.; Tang, Y.T.; Yao, A.J.; Cao, J.; Wu, Z.H.; Peng, Z.R.; Wang, S.Z.; Xiao, S.; Baker, A.J.M.; Qiu, L. Mitigation of Cd accumulation in paddy rice (Oryza sativa L.) by Fe fertilization. Environ. Pollut. 2017, 231, 549–559. [Google Scholar] [CrossRef]
  42. Loganathan, P.; Vigneswaran, S.; Kandasamy, J.; Naidu, R. Cadmium Sorption and Desorption in Soils: A Review. Crit. Rev. Environ. Sci. Technol. 2012, 42, 489–533. [Google Scholar] [CrossRef]
  43. Tahervand, S.; Jalali, M. Sorption, desorption, and speciation of Cd, Ni, and Fe by four calcareous soils as affected by pH. Environ. Monit. Assess. 2016, 188, 322. [Google Scholar] [CrossRef]
  44. Liu, G.; Wang, J.; Xue, W.; Zhao, J.; Wang, J.; Liu, X. Effect of the size of variable charge soil particles on cadmium accumulation and adsorption. J. Soils Sediments 2017, 17, 2810–2821. [Google Scholar] [CrossRef]
  45. Kong, L.; Guo, Z.; Peng, C.; Xiao, X.; He, Y. Factors influencing the effectiveness of liming on cadmium reduction in rice: A meta-analysis and decision tree analysis. Sci. Total. Environ. 2021, 779, 146477. [Google Scholar] [CrossRef] [PubMed]
  46. Wang, F.; Adams, C.A.; Shi, Z.; Sun, Y. Combined effects of ZnO NPs and Cd on sweet sorghum as influenced by an arbuscular mycorrhizal fungus. Chemosphere 2018, 209, 421–429. [Google Scholar] [CrossRef] [PubMed]
  47. He, P.P.; Lv, X.Z.; Wang, G.Y. Effects of Se and Zn supplementation on the antagonism against Pb and Cd in vegetables. Environ. Int. 2004, 30, 167–172. [Google Scholar] [CrossRef]
  48. Lin, L.; Zhou, W.; Dai, H.; Cao, F.; Zhang, G.; Wu, F. Selenium reduces cadmium uptake and mitigates cadmium toxicity in rice. J. Hazard. Mater. 2012, 235–236, 343–351. [Google Scholar] [CrossRef] [PubMed]
  49. Li, T.; Chang, Q.; Yuan, X.; Li, J.; Ayoko, G.A.; Frost, R.L.; Chen, H.; Zhang, X.; Song, Y.; Song, W. Cadmium transfer from contaminated soils to the human body through rice consumption in southern Jiangsu Province, China. Environ. Sci. Process. Impacts 2017, 19, 843–850. [Google Scholar] [CrossRef] [PubMed]
  50. Zhu, D.; Niu, Y.; Fan, K.; Zhang, F.; Wang, Y.; Wang, G.; Zheng, S. Selenium-oxidizing Agrobacterium sp. T3F4 steadily colonizes in soil promoting selenium uptake by pak choi (Brassica campestris). Sci. Total. Environ. 2021, 791, 148294. [Google Scholar] [CrossRef]
Figure 1. Location of 61 sampling sites in Shitai County, China.
Figure 1. Location of 61 sampling sites in Shitai County, China.
Agronomy 13 02378 g001
Figure 2. Distributions of Cd concentrations in (a) F1–F7 fractions by sequential extractions from soils, and in (b) rice grains and tea leaves.
Figure 2. Distributions of Cd concentrations in (a) F1–F7 fractions by sequential extractions from soils, and in (b) rice grains and tea leaves.
Agronomy 13 02378 g002
Figure 3. Relationships between Cd concentrations in rice grains/tea leaves and (a) soil solution Cd (Csolu), (b) effective concentration of Cd (CE) in soils.
Figure 3. Relationships between Cd concentrations in rice grains/tea leaves and (a) soil solution Cd (Csolu), (b) effective concentration of Cd (CE) in soils.
Agronomy 13 02378 g003
Figure 4. Factor importance derived from the random forest (RF) model for bioconcentration of Cd concentrations in (a) rice grains (BCFrice) and (b) tea leaves (BCFtea).
Figure 4. Factor importance derived from the random forest (RF) model for bioconcentration of Cd concentrations in (a) rice grains (BCFrice) and (b) tea leaves (BCFtea).
Agronomy 13 02378 g004
Table 2. Pearson’s correlation coefficients (r) between Cd accumulations in rice grains (n = 43)/tea leaves (n = 18) and measurements of Cd in soils, between BCFrice/BCFtea and soil properties.
Table 2. Pearson’s correlation coefficients (r) between Cd accumulations in rice grains (n = 43)/tea leaves (n = 18) and measurements of Cd in soils, between BCFrice/BCFtea and soil properties.
MeasurementsRiceTeaSoil PropertiesBCFriceBCFteaSoil PropertiesBCFriceBCFtea
DGT and DIFS (CE)0.84 b0.79 bpH−0.42 b−0.68 bZn−0.49 b0.09 NS
Soil solution (Csolu)0.68 b0.71 bSOC−0.12 NS0.15 NSSand−0.37 a−0.57 a
DGT (CDGT)0.61 b0.66 aCEC−0.09 NS−0.16 NSSilt0.20 NS−0.16 NS
MgCl2 extraction0.28 NS0.48 NSSe−0.41 b−0.51 aClay0.39 a0.49 a
NS, a, b Denote not significant, p < 0.05, p < 0.01, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wen, Y.; Wang, Y.; Tao, C.; Ji, W.; Huang, S.; Zhou, M.; Meng, X. Bioavailability of Cd in Agricultural Soils Evaluated by DGT Measurements and the DIFS Model in Relation to Uptake by Rice and Tea Plants. Agronomy 2023, 13, 2378. https://doi.org/10.3390/agronomy13092378

AMA Style

Wen Y, Wang Y, Tao C, Ji W, Huang S, Zhou M, Meng X. Bioavailability of Cd in Agricultural Soils Evaluated by DGT Measurements and the DIFS Model in Relation to Uptake by Rice and Tea Plants. Agronomy. 2023; 13(9):2378. https://doi.org/10.3390/agronomy13092378

Chicago/Turabian Style

Wen, Yubo, Yuanyuan Wang, Chunjun Tao, Wenbing Ji, Shunsheng Huang, Mo Zhou, and Xianqiang Meng. 2023. "Bioavailability of Cd in Agricultural Soils Evaluated by DGT Measurements and the DIFS Model in Relation to Uptake by Rice and Tea Plants" Agronomy 13, no. 9: 2378. https://doi.org/10.3390/agronomy13092378

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop