Next Article in Journal
The Association Between Fitness Test Scores and Musculoskeletal Injury in Police Officers
Previous Article in Journal
A Game-Theory Method to Design Job Rotation Schedules to Prevent Musculoskeletal Disorders Based on Workers’ Preferences and Competencies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enrichments of Cadmium and Arsenic and Their Effects on the Karst Forest Area

1
College of Resource and Environmental Engineering, Guizhou University, Guiyang 550025, China
2
State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences, Guiyang 550081, China
3
School of Tourism, Historical Culture, Liupanshui Normal University, Liupanshui 553004, China
4
Puding Karst Ecosystem Observation and Research Station, Chinese Academy of Sciences, Puding 562100, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(23), 4665; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16234665
Submission received: 24 October 2019 / Revised: 19 November 2019 / Accepted: 21 November 2019 / Published: 22 November 2019

Abstract

:
An understanding of the enrichment mechanisms of cadmium (Cd) and arsenic (As) in the process of rock weathering and soil formation is essential to develop agriculture according to local conditions. However, the enrichments of soil Cd and As under natural background conditions in karst areas are still uncertain. The enrichment factor, geo-accumulation index, redundancy analysis, and other methods were used to analyze the enrichment degree and the influencing factors of Cd and As on 5 rock–soil profiles and 15 topsoil samples, which were collected from a karst forest area in Libo County, Guizhou Province. The results showed that the enrichment process was divided into three stages. In the first stage, Cd and As were enriched in carbonate rocks, and their mean concentrations were 1.65 and 3.9 times those of the corresponding abundance of the crust. In the second stage, the enrichment of the parent rock into the soil, the enrichment factors of Cd and As in the parent material horizon relative to the bedrock horizon were 9.2 and 2.82, respectively. The third stage refers to the enrichments of Cd and As in the topsoil, where Cd enrichment was more obvious than that of As. Soil organic matter (SOM) and phosphorus (P) are important factors that influenced the enrichments of Cd and As in the topsoil. The functional groups of SOM were complexed with Cd and As; P easily formed precipitates with Cd, and the tree litter was fed back to the topsoil, which may be the reason for the surface enrichment of Cd and As. This study will help the scientific community understand the enrichment mechanisms of soil Cd and As in karst areas.

1. Introduction

Cadmium (Cd) frequently co-exists with arsenic (As) in soils [1]. The enrichments of Cd and As have always been an important scientific issue for environmental ecosystems. The karst area in southwest China is considered a geochemically anomalous area of Cd and As, with Guizhou province in the center. The mean Cd and As concentrations in the soil for this region were 0.37 mg/kg [2] and 19.82 mg/kg [3], respectively. These values were higher than the mean concentrations of Cd (0.27 mg/kg) and As (11.2 mg/kg) in China and considerably higher than the worldwide mean concentrations of 0.36 mg/kg and 6.0 mg/kg, respectively [4,5]. Cadimum and As are non-essential for the human body and are extremely dispersed potentially toxic elements (PTEs) in soils [6,7,8]. High concentrations of Cd and As in soils can seriously affect the growth and development of plants. Cadimum and As can also enter the human body through the food chain, which may affect human health [9,10,11,12,13]. The enrichment of Cd and As in soils, and their resulting ecological risks and potential threats to human health, have attracted widespread attention from local and foreign scholars.
Several studies have suggested that the enrichment of Cd and As in karst areas is due to the high geochemical background values of Cd and As in the soil [14,15]. A series of studies proposed that the enrichment of Cd and As is closely related to carbonate rocks [16,17], which have high contents of Cd and As. Other scholars have argued that a double enrichment mechanism is the source of Cd and As in karst areas. A notable change of rock/soil volume during carbonate rock weathering can easily lead to the relative enrichment of PTEs, such as Cd and As, from rock to soil [18]. In the process of rock weathering, the soil inherits the high background geochemical characteristics of the parent rock, which may form the high background values of Cd and As concentrations in the soil.
Some progress has been made in understanding the sources and enrichment characteristics of Cd and As in the soil. However, the specific enrichment process and formation mechanism of PTEs in karst remain unclear. Second, the many previous studies have only focused on topsoil; thus, research on PTE enrichment of the entire rock–soil profile is lacking. Third, the majority of studies focused on agricultural areas cultivated by humans, polluted industrial areas, and mining areas, and less on forest areas where human activities have had a weak impact. Based on these issues, this study takes a karst forest in Libo County, Guizhou Province, as its study area to explore the enrichment process of Cd and As in the soil throughout the rock–soil profile and its influencing factors on the karst forest ecosystem. Our findings provide a scientific basis for the spatial distribution and enrichment of Cd and As in the process of rock weathering and soil formation in karst areas.

2. Materials and Methods

2.1. Study Area and Samples Collection

The study area is located near the Maolan National Natural Reserve in Libo County, southeastern Guizhou, with a geographical location of 107°54′–108°02′ E, 25°12′–25°18′ N. This area is predominated by Carboniferous and Permian carbonate rocks. This forest area is also affected by anthropogenic activities to a lesser degree. The soil was mainly formed by weathering through limestone, dolomite, and argillaceous carbonate rocks. The study area is located in a mid-subtropical monsoon climate zone, with an average annual temperature of 15.3 °C and an annual precipitation of 1752.5 mm, mainly from April to October. The natural vegetation of the area is evergreen broad-leaved forest. In October 2018, we conducted field sampling and chose the representative sections in the karst mountainous areas, as well as in secondary forests, to excavate the soil profile. According to the principle of soil layering, the soil profile can be divided into the humus layer (Horizon A), the leached layer (Horizon B), the parent material layer (Horizon C), and the bedrock layer (Horizon D). A total of five soil profiles were collected, and three surface soil samples close to each profile were randomly taken at 0–15 cm (Figure 1). Surface soil samples were taken by the plum-shaped point method. Concretely, approximately 1.0 kg of topsoil samples were collected at 5 points, including the 4 vertices and the center of the (10 m × 10 m) square. After mixing, the samples were divided by quarters, and, finally, ~1.0 kg of soil samples were retained as the point. A total of 15 mixed surface soil samples were collected.

2.2. Determination of the Soil and Bedrock Samples

All soil samples were naturally air-dried for one month under room temperature conditions in the laboratory to remove grit and animal and plant residues. The samples were then ground and passed through a nylon sieve with a pore size of 0.25 mm. The bedrock was broken and treated in the same manner. Soil pH was determined using an acidity meter (HI98160, HANNA, USA) (water to soil ratio of 2.5:1), and soil organic matter (SOM) was determined by the potassium dichromate volumetric method [19]. Humic acid (HA) was determined by the NaOH and Na4P2O7 solution extraction and K2CrO4 solution oxidation methods [20]. Total phosphorus (P) was determined by the sulfuric acid–perchloric acid digestion–molybdenum blue colorimetric methods [21]. Soil total nitrogen (N) was determined by the Kjeldahl determination method [22]. Soil particle size distribution was analyzed using an international soil quality classification standard [23]. Thereafter, the 2 mm air-dried soil samples were divided into three types—clay, silt, and sand—using the hydrometer method.
The soil Cd was analyzed using inductively coupled plasma mass spectrometry (NexION300X, PerkinElmer, Waltham, MA, USA), and As was measured using an atomic fluorescence spectrophotometer (AFS-920, Beijing Titan Instruments Co. Ltd., China), with detection limits of 0.01 and 0.05 mg/kg for Cd and As, respectively. Approximately 100 mg of the soil samples were digested in a high-pressure Teflon vessel using a HF + HNO3 (ultra-pure, 0.5 mL +5 mL) mixture for 48 h (150 °C). Exactly 2 mL of the digests were diluted to 50 mL in deuterium depleted water for testing. The quality controls consisted of method blanks, certified reference materials (CRMs), and blind duplicates. The CRMs were in agreement with recommended values (with the recovery rate ranging from 92% to 110%), and the analytical precision was generally better than 5% for all elements.

2.3. Geo-Accumulation Index

In 1986, the German scientist Müller introduced the geo-accumulation index (Igeo) [24]. This method considers human factors, environmental geochemical background values, and the influence of natural diagenesis background values, which can be used to reflect the degree of PTE enrichment in soils or other materials. The calculation formula was presented as follows:
I g e o = log 2 ( C i k B i )
where Ci is the measured concentrations of Cd or As, and Bi is the geochemical background values of Cd and As. k is a constant, which considers that diagenesis may initiate a change in the background values (1.5 was frequently used). The background values of Cd and As in the surface soil of Guizhou were 0.31 and 14.89 mg/kg, respectively, and these background values are representative for Guizhou province [25]. Table 1 presents the enrichment degree corresponding to the Igeo.

2.4. Statistical Analysis

To evaluate the effect of the physicochemical properties of soil on the concentrations of Cd and As, we used the R3.5.3 software (MathSoft Company, Cambridge, Massachusetts, USA) to analyze the correlation among Cd, As, and the influencing factors. Pearson’s correlation is a statistical method for studying the correlation between random variables to reflect the closeness of the relationship between variables. The Kolmogorov–Smirnova and Shapiro–Wilk methods were used to test the normality of Cd and As contents in surface soil based on the IBM SPSS 20 software (IBM Corporation, Armonk, NY, USA). To further explore the relationship between the PTEs and the physicochemical properties of soil within a spatial distribution, a redundant analysis was performed using the Canoco 4.5 software (Microcomputer Power, Ithaca, NY, USA). In addition, to construct a stepwise regression model, we performed a stepwise regression analysis using IBM SPSS 20 software.

3. Results and Discussion

3.1. Distribution Characteristics of Cd and As in Rock–Soil Profiles

The Cd and As contents of the different layers of soil and bedrock in the rock–soil profiles of the study area were 5.44 and 27.74 mg/kg (horizon A) > 4.21, 24.62 mg/kg (horizon B) > 3.09, 24.26 mg/kg (horizon C) > 0.33, and 8.54 mg/kg (horizon D) (Figure 2).
In general, the Cd and As concentrations were maximal and minimal in Horizons A (surface soil) and D (lower bedrock), respectively. The mean concentrations of As and Cd in the limestone were 3.9 and 1.65 times the abundance of crust at 2.2 and 0.2 mg/kg, respectively [26]. This finding indicates that the geological background of the limestone in the study area has a high geochemical background of Cd and As. A previous study showed that the Cd concentration of braided carbonate rocks of the Middle Jurassic in Switzerland points to extreme enrichment, as the highest concentration of Cd was 4.9 mg/kg [27]. The Cd concentration of the Permian limestone in central Japan was relatively high, ranging from 0.17 to 2.41 mg/kg, with an average concentration of 0.5 mg/kg [28]. The results obtained in this study are similar to those of previous studies.
The Cd and As contents of the soil samples from layers A, B, and C far exceed the soil limit values (0.20 and 15 mg/kg) of soil environmental quality standards (GB15618-1995) in China. Furthermore, these values are higher than the background values for Cd and As in Guizhou soil [25], thereby indicating that the soil in the study area has high Cd and As geochemical backgrounds. The study showed that the Cd concentrations of soils in cultivated land (1.33 mg/kg) and forest land (1.57 mg/kg) in the karst distribution area of Luodian County, southern Guizhou, were significantly greater than those of non-karst rock areas (0.27 mg/kg) and woodland soil (0.22 mg/kg). The high concentration of Cd in the soil was clearly controlled by the distribution of carbonate rocks [29]. In a small area of southwest Guizhou, the concentration of As in agricultural land in the karst area (47.9 mg/kg) was also significantly greater than those in the sub-karst area (36.8 mg/kg) and in non-karst area (12.7 mg/kg) [30]. In addition, the concentrations of As in the overlying soil in Guizhou and carbonate rock were high, whereas the concentration of As in clastic rocks was low. Furthermore, the background value of carbonate rocks for As was high, and As easily produced relative enrichment in the process of rock weathering and soil formation. These are the main reasons for high As content in carbonate rock distribution areas [3,31]. Cadmium and As in soil developed by carbonate parent materials showed significant enrichment characteristics in southern Guangxi [32]. In summary, the Cd and As contents of the carbonated weathered soil profiles show a high background distribution, where the Cd and As enrichments of the topsoil are notable in forest areas weakly affected by human activities.

3.2. Enrichments of Cd and As Across Soil Layers

As mentioned, Cd and As exhibit relative enrichment characteristics in carbonate rocks relative to the crust, which provides a high background value for weathered soil.
The enrichment factor of Cd and As contents in the A, B, and C layers of the five soil profiles with respect to bedrock horizon (D) was calculated to explore the relative enrichment degrees of Cd and As during soil formation through the weathering of carbonate rocks (Figure 3). Figure 3 shows that the enrichment factor of Cd in Horizon C relative to Horizon D was between 4.58 and 12.9 with a mean value of 9.2, whereas that of As was between 1.9 and 3.64 with a mean value of 2.82. This finding indicates that Cd and As were relatively enriched and showed “second-stage enrichment” when Horizon C was formed by the weathering of carbonate rocks.
Cadmium was more enriched than As. Compared with Horizon C, Horizon B had an enriched Cd. In profiles one and five, the Cd concentration in Horizon B was slightly greater than that in Horizon A. On the other hand profiles two, three, and four showed that Horizon A was significantly greater than horizon B.
In soil profiles one and three, the concentration of As in Horizon B was slightly lower than that in Horizon C, whereas the opposite was true for profiles two, four, and five—that is, Horizon B was higher than Horizon C. The enrichment degree of As in Horizon A was higher than that in Horizons B and C. In summary, Cd and As in surface soil were relatively enriched with respect to their parent material layer, and the average relative enrichment factors were 1.78 and 1.14, respectively. Results indicate that Cd and As were further enriched during the development of humus (Horizon A) from the parent material layer (Horizon C).
The results of the Igeo method show that the soil in each layer of the five profiles has greater than moderate enrichment (Figure 3), with moderate to heavy enrichment accounting for 33.3% and heavy enrichment accounting for 60%. The ratio of As in the soil for each layer with non-enrichment to moderate enrichment was 66.7%, and the remainder pertained to non-enrichment. In summary, humus (Horizon A) is more concentrated than the parent material layer (Horizon C), and Cd and As are concentrated in the surface soil, which indicates “third-stage enrichment”.

3.3. Effects of the Physicochemical Properties of Topsoil on Cd and As Contents

Cadmium and As show notable enrichment in Horizon A. Several studies have shown that PTEs exhibit surface aggregation in the soil profile [33,34]. Similarly, PTE enrichment and pollution mainly affect the topsoil. Therefore, the effects of the physicochemical properties of the surface on soil Cd and As enrichment were analyzed and are provided in Table 2, as well as Figure 4 and Figure 5.
The total concentrations of Cd and As in the surface soil ranged from 2.67 to 5.92 mg/kg (n = 15) and from 14.90 to 36.20 mg/kg, respectively. The contents of both Cd and As obeyed normal distribution, which is based on the Kolmogorov–Smirnova and Shapiro–Wilk of methods. The average concentration for Cd was 3.83 ± 0.93 mg/kg, while the mean value for As was 23.8 ± 5.6 mg/kg. The average concentrations of Cd and As in the topsoil were 12.4 and 1.6 times the background values of soil Cd and As in Guizhou, respectively. The average pH value was 7.54 ± 0.35, which is slightly alkaline. The SOM concentration reached an average value of 67.35 ± 17.39 g/kg. The mean values for total P and N concentrations were 0.61 ± 0.20 and 4.02 ± 1.20 mg/kg, respectively. Total HA concentration ranged from 9.61 to 28.65 mg/kg. The soil particle size distribution showed that the surface soil was composed of clay loam.
Figure 4 depicts the correlation analysis among the eight indicators of the physicochemical properties of the topsoil and concentrations of Cd and As. The figure shows that Cd in the topsoil was significantly positively correlated with N, P, HA, and SOM, with correlation coefficients of 0.86, 0.66, 0.85, and 0.86, respectively. Total Cd was positively correlated with sand, with a correlation coefficient of 0.49. The mineral composition of sand was the residual particles of limestone weathering. In addition, the mineral composition of the grains was mainly calcite with a chemical composition of mainly CaCO3. Cd could easily replace Ca in the form of an isomorphism and enter the crystal lattice of calcite minerals, which would be beneficial for the accumulation of Cd.
Second, the topsoil was alkaline; thus, the migration was weak, and Cd easily accumulated in the topsoil. However, the total Cd in the topsoil was negatively correlated with pH and the proportion of clay in the soil particulate distribution; the correlation coefficients were −0.42 and −0.38. It can be seen that the soil alkalinity and soil aggregate particulate distribution were not unique factors affecting PTE concentrations. At the same time, no significant correlation was observed between the total As and soil physicochemical parameters, thereby indicating that Cd was more likely to accumulate in the topsoil compared with As.
The results of redundancy analysis (RDA) showed a notable separation among the 15 sample points on the surface (Figure 5), which were distributed in the four quadrants and indicated that, although the soil type remained the same, the topsoil exhibited a certain degree of spatial heterogeneity.
Figure 5 shows that the concentrations of Cd and As in the topsoil were in the same quadrants as N, P, HA, and SOM. The direction of the arrow remains the same, and the angle between each factor is acute, which shows that surface soil concentrations of Cd and As were mainly influenced by N, P, HA, and SOM. In addition, the angle between the sand was acute, which indicates that it was also affected by the composition of sand in the soil. PTEs are present in the soil particles. Their amount in the soil has a considerable influence on the distribution and migration of heavy metals, which indicates that the concentration of heavy metals in soil was closely related to soil texture. In summary, the enrichment of Cd in surface soil was mainly influenced by SOM, N, HA, and P. As is mainly affected by SOM and P.

3.4. Regression Analysis

Studies have shown that correlation may not necessarily indicate a significant causal association. The stepwise regression method gradually selects variables according to their contribution values and obtains the important influencing variables by eliminating factors with non-significant contribution values. The regression equations for Cd and As constructed by the stepwise regression model were derived as follows:
Y1 = 6.298 + 1.679X2 + 0.032X4 − 0.037X6 − 0.566X (F = 14.78, p < 0.01; R = 0.925; Ry1 = 0.617, Ry2 = 0.761, Ry3 = −0.424, Ry4 = −0.41; t1 = 2.477, t2 = 3.711, t3 = 1.481, t4 = 1.431; P1 = 0.031, P2 = 0.003, P3 = 0.167, P4 = 0.180)
Y2 = 0.416 + 1.395X2 + 0.038X4 (F = 25.1; R = 0.899, Ry1 = 0.511, Ry2 = 0.812; t1 = 2.058, t2 = 4.831; P1 = 0.060, P2 = 0.003).
The F test results show that the regression equation of Y1 reached a very significant level. The significance test of the partial regression coefficient indicated that P and SOM were the most important factors that influenced Cd concentration in soil. Surface Cd concentration increases with an increase in the P and SOM concentrations in soil. The regression equation of Y2 did not reach a significant level. Although the concentration of As showed a non-significant linear relationship with SOM and P, As was apparently mainly influenced by SOM and P. Therefore, SOM and P were the key factors for Cd and As enrichment in forest soil in karsts.

3.5. Analysis of Cd and As Enrichment Mechanisms

3.5.1. Factors of Geological Background

The concentrations of Cd and As in the study area exceeded the standards for soil environmental quality in China and the corresponding PTE background values in Guizhou Province. The research area is a forest without mining activities and is less affected by human activities. Therefore, the enrichments of both Cd and As were not a result of human factors. We speculate that enrichment could be mainly controlled by the natural environment.
The concentrations of Cd and As in the bedrock were 1.65 and 3.9 times the crustal abundance in the study area. Both values were greater than one, thereby indicating that Cd and As had high background values in carbonate rocks. This finding further indicates that the enrichments of Cd and As in the karst area may be restricted by carbonate lithology [16]. Rambeau et al. [35] argued that the concentration of Cd in Jurassic limestone in the Lower Burgundy region of France reached as high as 2.6 mg/kg, which indicates that Cd enrichment in the karst region is closely related to the high background values of carbonate rocks. He et al. [36] pointed out that the high concentration of Cd in the area was mainly distributed in exposed carbonate rocks. Ni et al. [37] noted the enrichment trend of As in carbonate rocks of the Lower Paleozoic in Beijing by the symbiotic soil covering the rocks. In addition, the results of the 1:200,000 regional geochemical sweeping in Yunnan Province showed that the regions with geochemically high contents of As and Cd had good spatial coupling with widely distributed carbonates [14].
The concentrations of Cd and As in Horizon C of the five rock–soil profiles in the study area were relatively enriched, thereby indicating second-stage enrichment. Previous studies have shown that overlying soils and underlying carbonate rocks in Switzerland clearly have a shared provenance with the underlying carbonate rocks. The high Cd concentration in the soil was closely related to the high concentrations of Cd in the bedrock [38]. In addition, the ratio of Cd to limestone bedrock concentration in the Lower Burgundy region of France was between 4.6 and 5.7, which indicates that soil developed by carbonate weathering exhibits obvious Cd enrichment characteristics [35]. The weathering of carbonate rocks is an important source of Cd and As enrichment for soil in the karst area of Yunnan. The accumulation of As and Cd from rock to soil was approximately 4 to 5 times, and the enrichment mechanism was mainly due to the adsorption of clay minerals and organic matter [39]. A study also showed that soil formed by the weathering of carbonate rocks leads to the relative enrichment characteristics of Fe, Al, Mn, and As compared with the bedrock. Conversely, the concentrations of Ca and Mg are relatively deficient in the bedrock [40].
On the one hand, the high concentrations of Cd and As in soil are related to the high background values of carbonate rocks. On the other hand, this phenomenon may be due to clay mineral enrichment, such as the oxides and hydroxides of Fe and Al, which may cause the enrichment of Cd and As concentrations in the soil parent material compared to the bedrock.

3.5.2. Effects of Plant Growth Processes

The study area is a forest with vigorous vegetation. The Cd and As in the soil absorbed by plants were returned to the surface soil in the form of litter, which provides a source for the accumulation of Cd and As on the soil’s surface [41].
Excess phosphates in the study area may induce the precipitation of Cd3(PO4)2 and Cd2+ [42,43]. The competitive adsorption between P and As in soil can reduce the adsorption of As in soil with an increase in P concentrations [44]. Therefore, Cd in the topsoil is more concentrated than As. The forest area is flourishing, and SOM is largely accumulated under biological actions [45,46]. A large number of studies have shown that SOM has specific complexation and accumulation capacities for heavy metals in soil because of its abundant functional groups [47,48,49]. SOM is an important adsorbent of heavy metals in soil because it can retain heavy metals in soil by forming metal–SOM complexes to reduce the mobility of heavy metals [50]. Lin et al. [51] showed that removing SOMs resulted in an average decrease of 20% in the adsorption of Cd2+. Luo et al. [52] pointed out that concentrations of Cd in soil generally decreased with a deepening of stony desertification in a karst, which indicates that the concentration of heavy metals decreased with a decrease in SOM.
Therefore, topsoil Cd and As were mainly accumulated due to the large amount of SOM complexed with heavy metals formed by plant growth and the precipitation of Cd by P and feedback litter.

4. Conclusions

We used the enrichment factor, geo-accumulation index, redundancy analysis, and other methods to analyze the enrichment mechanism of Cd and As in the process of rock weathering and soil formation in the karst area. We drew the following conclusions:
(1) The mean concentrations of Cd and As in the lime soil profile were characterized as follows: horizon A > horizon B > horizon C > horizon D. The entire enrichment process was mainly divided into three stages. In the first stage, carbonate was enriched in the Cd and As concentrations relative to the crust. The second stage featured the enrichment of Cd and As in the process of rock weathering and soil formation. The third stage entailed the topsoil enrichment of Cd and As.
(2) The enrichment process and its influencing factors differed across stages. The first stage was mainly influenced by lithology. In the second stage, the adsorption of clay minerals, such as the oxides and hydroxides of Fe and Al, led to the enrichment of Cd and As concentrations in the soil parent material in comparison to the bedrock. Cadimum and As surface aggregation, especially the enrichment of Cd, was more obvious in the third stage. Complexation and precipitation could be the main enrichment mechanisms for Cd and As in topsoil.
(3) This study revealed the enrichment process of soil Cd and As in karst areas. The results could provide a valuable reference for the source apportionment of soil Cd and As in the karstic geochemical anomalous area of Southwest China.

Author Contributions

F.L. designed the experiments; F.C., Q.L., and J.W. did the field sampling; X.B. and J.Z. analyzed the data; J.W. wrote paper; all authors read and approved the final version.

Funding

This research work was supported jointly by the national key research program of China (No. 2016YFC0502300 and 2016YFC0502102), National Natural Science Foundation (No. U161244), Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDA23060100), Western Light Talent Program (Category A) (No. 2018-99), United fund of karst science research center (No. U1612441), Science and Technology Plan of Guizhou Province of China (2017–2966), Chinese academy of science and technology services network program (No. KFJ-STS-ZDTP-036), and international cooperation agency international partnership program (No.132852KYSB20170029, No. 2014-3), and International cooperation research projects of the national natural science fund committee (No. 41571130074 & 41571130042).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yuan, Y.; Chai, L.; Yang, Z.; Yang, W. Simultaneous immobilization of lead, cadmium, and arsenic in combined contaminated soil with iron hydroxyl phosphate. J. Soil. Sediments 2017, 17, 432–439. [Google Scholar] [CrossRef]
  2. Luo, H.; Liu, Q.M.; Wang, S.J. Pollution characteristics and sources of cadmium insoils of the karst area in South China. Chin. J. Ecol. 2018, 37, 1538–1544. [Google Scholar]
  3. Wang, H.; Liu, X.M.; Liu, F.; Wang, S.J.; Tang, Q.L. Constraints of geological lithology in the spatial distribution of arsenic. Ceochimica 2019, 48, 114–125. [Google Scholar]
  4. Zhang, X.; Chen, D.; Zhong, T. Assessment of cadmium (Cd) concentration in arable soil in China. Environ. Sci. Pollut. Res. 2015, 22, 4932–4941. [Google Scholar] [CrossRef] [PubMed]
  5. Kubier, A.; Wilkin, R.T.; Pichler, T. Cadmium in soils and groundwater: A review. Appl. Geochem. 2019, 108, 104–107. [Google Scholar] [CrossRef]
  6. Gurung, B.; Race, M.; Fabbricino, M.; Komínková, D.; Libralato, G.; Siciliano, A.; Guida, M. Assessment of metal pollution in the Lambro Creek (Italy). Ecotoxicol. Environ. Saf. 2018, 148, 754–762. [Google Scholar] [CrossRef] [PubMed]
  7. Race, M.; Ferraro, A.; Fabbricino, M.; La Marca, A.; Panico, A.; Spasiano, D.; Tognacchini, A.; Pirozzi, F. Ethylenediamine-N,N′-Disuccinic Acid (EDDS)—Enhanced Flushing Optimization for Contaminated Agricultural Soil Remediation and Assessment of Prospective Cu and Zn Transport. Int. J. Environ. Res. Public Health 2018, 15, 543. [Google Scholar] [CrossRef] [PubMed]
  8. Liang, Y.; Yi, X.; Dang, Z.; Wang, Q.; Luo, H.; Tang, J. Heavy Metal Contamination and Health Risk Assessment in the Vicinity of a Tailing Pond in Guangdong, China. Int. J. Environ. Res. Public Health 2017, 14, 1557. [Google Scholar] [CrossRef] [PubMed]
  9. Obiri, S.; Yeboah, P.O.; Osae, S.; Adu-Kumi, S. Levels of arsenic, mercury, cadmium, copper, lead, zinc and manganese in serum and whole blood of resident adults from mining and non-mining communities in Ghana. Environ. Sci. Pollut. Res. Int. 2016, 23, 16589–16597. [Google Scholar] [CrossRef] [PubMed]
  10. Pan, J.; Plant, J.A.; Voulvoulis, N.; Oates, C.J.; Ihlenfeld, C. Cadmium levels in Europe: Implications for human health. Environ. Geochem. Health 2010, 32, 1–12. [Google Scholar] [CrossRef]
  11. Rangkadilok, N.; Siripriwon, P.; Nookabkaew, S.; Suriyo, T.; Satayavivad, J. Arsenic, Cadmium, and Manganese Levels in Shellfish from Map Ta Phut, an Industrial Area in Thailand, and the Potential Toxic Effects on Human Cells. Arch. Environ. Contam. Toxicol. 2015, 68, 169–180. [Google Scholar] [CrossRef] [PubMed]
  12. Kumasaka, M.Y.; Yamanoshita, O.; Shimizu, S.; Ohnuma, S.; Furuta, A.; Yajima, I.; Nizam, S.; Khalequzzzaman, M.; Shekhar, H.U.; Nakajima, T.; et al. Enhanced carcinogenicity by coexposure to arsenic and iron and a novel remediation system for the elements in well drinking water. Arch. Toxicol. 2013, 87, 439–447. [Google Scholar] [CrossRef] [PubMed]
  13. Yanez, L.; Garcia-Nieto, E.; Rojas, E.; Carrizales, L.; Mejia, J. DNA damage in blood cells from children exposed to arsenic and lead in a mining area. Environ. Res. 2003, 93, 231–240. [Google Scholar] [CrossRef] [PubMed]
  14. Li, H.; Wang, B.L. Geochemical characteristics of as and Cd elements in soils of YunnanProvince. Geophys. Geochem. Explor. 2008, 21, 245–250. [Google Scholar]
  15. Mandal, B.K.; Suzuki, K.T. Arsenic round the world: A review. Talanta 2002, 58, 226–235. [Google Scholar] [CrossRef]
  16. Rambeau, C.; Föllmi, K.B.; Matera, V.; Adatte, T.; Steinmann, P. Cadmium enrichments in jurassic carbonates: Causes and mechanisms. Geochim. Et Cosmochim. Acta Suppl. 2017, 69, 406–412. [Google Scholar]
  17. Shao, S.; Ye, L. Geochemical characteristics of carbonate rock weathering and forming soil in the center area of guizhou province. Chin. Sci. Bull. 1999, 20, 218–219. [Google Scholar]
  18. Wang, S.J.; Ji, H.B.; Ouyang, Z.Y.; Zheng, H.P.; Ni, T.Y. Preliminary study on weathering and pedogenesis of carbonate rock. Sci. China Ser. D Earth Sci. 1999, 42, 572–581. [Google Scholar] [CrossRef]
  19. Zhang, N.; He, H.-D.; Wu, W. Studies on the characteristics of soil organic matter and pedogenic calcium carbonate for three kinds of soil in the tengri desert. Acta Ecol. Sin. 2009, 29, 4094–4101. [Google Scholar]
  20. Tregubova, P.N.; Koptsik, G.N.; Stepanov, A.A.; Smirnova, I.E.; Kadulin, M.S.; Turbaevskaya, V.V.; Zakharenko, A.L. Application of humic substances in the remediation of heavy-metal-polluted soils of the subarctic zone of the kola peninsula. Mosc. Univ. Soil Sci. Bull. 2017, 72, 207–214. [Google Scholar] [CrossRef]
  21. Fang, L.; Teng, B.H.; Qian, X.G.; Hai, B.L.; Ying, G.S. The status of phosphorus and its environmental effects analysis under various land use patterns in yellow soil of hilly areas. Soil Environ. Sci. 2002, 11, 232–236. [Google Scholar]
  22. Bremner, J.M. Determination of nitrogen in soil by the kjeldahl method. J. Agric. Sci. 1960, 55, 11–33. [Google Scholar] [CrossRef]
  23. Wilson, M.J. The importance of parent material in soil classification: A review in a historical context. Catena 2019, 182, 104131. [Google Scholar] [CrossRef]
  24. Müller, V.G. Schadstoffe in sedimenten–Sedimenteals Schadstoffe. Mitt Österr. Geol. Ges 1986, 79, 107–126. [Google Scholar]
  25. He, S.L. Geochemical background of supergene sediments in Guizhou. Guizhou Geol. 1998, 15, 149–156. [Google Scholar]
  26. Li, T. Geochemical abundance of chemical elements. Geochemistry 1976, 3, 167–174. [Google Scholar]
  27. Quezada-Hinojosa, R.P.; Matera, V.; Adatte, T.; Rambeau, C.; Föllmi, K.B. Cadmium distribution in soils covering jurassic oolitic limestone with high cd contents in the swiss jura. Geoderma 2009, 150, 287–301. [Google Scholar] [CrossRef]
  28. Aizawa, S.; Akaiwa, H. Cadmium contents of Triassic and Permian limestones in central Japan. Chem. Geol. 1992, 98, 103–110. [Google Scholar] [CrossRef]
  29. Tang, Q.L.; Liu, X.M.; Liu, F.; Wang, H.; Wang, S.J. Accumulation of Cd and Risk Assessment in Arable Soil in Karst Region of Northern Luodian in Guizhou. Environ. Sci. 2019, 10, 1–13. [Google Scholar]
  30. Wang, H.; Liu, X.M.; Liu, F.; Tang, Q.L.; Wang, S.J. Spatial Distribution and Pollution Assessment of As at a Small Scale in Agricultural Soils of the Karst Region. Environ. Sci. 2019, 40, 2895–2902. [Google Scholar]
  31. Vital, M.; Martínez, D.E.; Babay, P.; Quiroga, S.; Clément, A.; Daval, D. Control of the mobilization of arsenic and other natural pollutants in groundwater by calcium carbonate concretions in the Pampean Aquifer, southeast of the Buenos Aires province, Argentina. Sci. Total Environ. 2019, 15, 532–543. [Google Scholar] [CrossRef] [PubMed]
  32. Ma, H.H.; Peng, M.; Liu, F.; Guo, F.; Tang, S.Q.; Liu, X.J.; Zhou, Y.L.; Yang, K.; Li, K.; Yang, Z.; et al. Bioavailability, Translocation, and Accumulation Characteristic of Heavy Metals in Soil-Crop System from a Typical Carbonate Rocks Area in Guangxi, China. Environ. Sci. 2019. [Google Scholar] [CrossRef]
  33. Emurotu, J.E.; Onianwa, P.C. Bioaccumulation of heavy metals in soil and selected food crops cultivated in Kogi State, north central Nigeria. Environ. Syst. Res. 2017, 6, 21. [Google Scholar] [CrossRef]
  34. Novak, J.M.; Watts, D.W.; Stone, K.C. Copper and zinc accumulation, profile distribution, and crop removal in coastal plain soils receiving long-term, intensive applications of swine manure. Trans. ASABE 2004, 47, 1513. [Google Scholar] [CrossRef]
  35. Rambeau, C.M.C.; Baize, D.; Saby, N.; Matera, V.; Adatte, T.; Föllmi, K.B. High cadmium concentrations in jurassic limestone as the cause for elevated cadmium levels in deriving soils: A case study in lower burgundy, france. Environ. Earth Sci. 2010, 61, 1573–1585. [Google Scholar] [CrossRef]
  36. He, S.L.; Long, C.L.; Liu, Y.Z. Geochemistry and environment of the cadmium in the soil and Sediments at the surface of Guizhou Province. Guizhou Geol. 2004, 21, 245–250. [Google Scholar]
  37. Ni, S.Q.; Hou, Q.L.; Ju, Y.W.; Wang, S.; Liu, Q.; Wu, Y.; Xiao, L. Migration and enrichment of trace elements in the Lower Paleozoic carbonate strata in Beijing. Prog. Nat. Sci. 2009, 39, 488–496. [Google Scholar]
  38. Wang, J.; Wang, S.J. The sources and crop effects of heavy metal elements of contamination in soil. J. Guizhou Normal. Uni. 2005, 2, 117–124. [Google Scholar]
  39. Wang, Y.; Peng, S.H.; Yang, S.L. The anomaly features of as and Cd in the karst area in Yunan province. Carsol. Sin. 2012, 31, 377–381. [Google Scholar]
  40. Goh, K.H.; Lim, T.T. Arsenic fractionation in a fine soil fraction and influence of various anions on its mobility in the subsurface environment. Appl. Geochem. 2005, 20, 231–239. [Google Scholar] [CrossRef]
  41. Sung, U.K.; Vancen, O.; Yong, G.K.; Sang, M.L. Effect of phosphate addition on Cadmium precipitation and adsorpti on incontaminated arable soil with a low concentration of Cadmium. Bull. Environ. Contam. Toxicol. 2015, 95, 675–679. [Google Scholar]
  42. Gramlich, A.; Tandy, S.; Gauggel, C.; López, M.; Perla, D.; Gonzalez, V. Soil cadmium uptake by cocoa in Honduras. Sci. Total Environ. 2018, 612, 370–378. [Google Scholar] [CrossRef] [PubMed]
  43. Penido, E.S.; Melo, L.C.A.; Guilherme, L.R.G.; Bianchi, M.L. Cadmium binding mechanism sand adsorption capacity by novel phosphorus/magnesium-engineered biochars. Sci. Total Environ. 2019, 25, 1134–1143. [Google Scholar] [CrossRef]
  44. Deng, T.T.; Duan, H.X.; Liu, B.B.; Ma, P.; Sun, X.M.; Fan, J.H. Effects of Four Typical Heavy Metal Pollutions on Soil Phosphorus Adsorption. Chin. J. Soil Sci. 2019, 50, 463–473. [Google Scholar]
  45. Lal, R. Forest soil sand carbon sequestration. For. Ecol. Manag. 2005, 220, 242–258. [Google Scholar] [CrossRef]
  46. Sedjo, R.; Sohngen, B. Carbon sequestration in forests and soils. Ann. Rev. Resour. Econ. 2012, 4, 127–144. [Google Scholar] [CrossRef]
  47. De Mora, A.P.; Ortega-Calvo, J.J.; Cabrera, F.; Madejón, E. Changes in enzyme activities and microbial biomass after ‘’in situ’’ remediation of a heavy metal-contaminated soil. Appl. Soil Ecol. 2005, 28, 125–137. [Google Scholar] [CrossRef]
  48. Kleber, M.; Mikutta, R.; Torn, M.S.; Jahn, R. Poorly crystalline mineral phases protect organic matter in acid subsoil horizons. Eur. J. Soil Sci. 2005, 56, 717–725. [Google Scholar] [CrossRef]
  49. Tiberg, C.; Sjöstedt, C.; Gustafsson, J.P. Metal sorption to spodosol bs horizons: Organic matter complexes predominate. Chemosphere 2018, 196, 556–565. [Google Scholar] [CrossRef]
  50. Matijevic, L.; Romic, D.; Romic, M. Soil organic matter and salinity affect copper bioavailability in root zone and uptake by Vicia fabaL.plants. Environ. Geochem. Health 2014, 36, 883–896. [Google Scholar] [CrossRef]
  51. Liang, J.; Xu, R.K.; Jiang, X.; Bian, Y.R.; Tan, W.F. Adsorption and Desorption of Cu(II), Pb(II) and Cd(II) in Two Variable Charge Soils Different in pH. Soils 2007, 39, 992–995. [Google Scholar]
  52. Luo, X.Q.; Wang, S.J.; Liu, X.M.; Wang, C.Y.; Yu, H. Geochemical characteristics of the soil cadmium In the process of karst rocky desertification. Ecol. Environ. Sci. 2009, 18, 160–166. [Google Scholar]
Figure 1. Geographical location and sampling point distribution in the study area.
Figure 1. Geographical location and sampling point distribution in the study area.
Ijerph 16 04665 g001
Figure 2. The distribution characteristics of Cd and As content across soil layers. The number of statistical samples per layer is 5. The bars are the mean standard deviation.
Figure 2. The distribution characteristics of Cd and As content across soil layers. The number of statistical samples per layer is 5. The bars are the mean standard deviation.
Ijerph 16 04665 g002
Figure 3. Geo-accumulation indexes of Cd and As content across soil layers and the bedrock ratios in five soil profiles.
Figure 3. Geo-accumulation indexes of Cd and As content across soil layers and the bedrock ratios in five soil profiles.
Ijerph 16 04665 g003
Figure 4. Correlation analysis of soil Cd, As, and soil physicochemical parameters (* manifest p < 0.05, ** manifest p < 0.01, *** manifest p <0.001).
Figure 4. Correlation analysis of soil Cd, As, and soil physicochemical parameters (* manifest p < 0.05, ** manifest p < 0.01, *** manifest p <0.001).
Ijerph 16 04665 g004
Figure 5. Bioplot map of hybrid redundancy analysis (RDA) for potentially toxic elements (PTEs) and environmental factors.
Figure 5. Bioplot map of hybrid redundancy analysis (RDA) for potentially toxic elements (PTEs) and environmental factors.
Ijerph 16 04665 g005
Table 1. Geo-accumulation index.
Table 1. Geo-accumulation index.
ClassIgeo ValueEnrichment levelClassIgeo ValueEnrichment Level
0Igeo ≤ 0Non-enrichment43 < Igeo ≤ 4Heavy enrichment
10 < Igeo ≤ 1Non-enrichment to Moderate enrichment54 < Igeo ≤ 5Heavy to Extreme enrichment
21 < Igeo ≤ 2Moderate enrichment6Igeo > 5Extreme enrichment
32 < Igeo ≤ 3Moderate to heavy enrichment
Table 2. Cadmium (Cd) and Arsenic (As) contents and physicochemical properties of soil.
Table 2. Cadmium (Cd) and Arsenic (As) contents and physicochemical properties of soil.
TopsoilAsCdNPHASOMClaySlitSandpH
121.44.095.080.2823.1782.321335527.54
226.34.274.520.7321.3182.821936457.89
325.03.944.190.4624.5280.062737367.24
430.92.713.500.6916.8654.641941407.94
536.24.033.690.5716.5555.592435417.74
629.93.063.470.5913.5553.002635397.72
726.15.925.741.0028.6592.441442447.18
824.85.236.000.9628.1790.721140497.62
922.54.295.340.6426.6987.571044467.06
1020.54.133.850.6417.3863.642933387.56
1119.83.983.940.6618.4164.491731527.08
1222.33.253.450.4320.2763.751550357.01
1317.52.843.240.5613.8651.622039417.90
1418.22.671.210.299.6235.852540357.59
1514.92.973.150.6013.2451.781839438.03
Max36.2 5.92 6.00 1.00 28.65 92.44 29 50 52 8.03
Min14.9 2.67 1.21 0.28 9.62 35.85 10 31 35 7.01
Mean23.8 3.83 4.02 0.61 19.48 67.35 19 38 42 7.54
SD5.6 0.93 1.20 0.20 5.85 17.39 6 5 6 0.35
Note: The units for As, Cd, N, P, and HA are mg/kg; the unit for SOM is g/kg; the units for Clay, Slit, and Sand are %. N (Nitrogen); P (Phosphorus); HA (Humic acid); SOM (Soil organic matter). Max (Maximum); Min (Minimum); SD (Standard deviation).

Share and Cite

MDPI and ACS Style

Wang, J.; Bai, X.; Liu, F.; Zhang, J.; Chen, F.; Lu, Q. Enrichments of Cadmium and Arsenic and Their Effects on the Karst Forest Area. Int. J. Environ. Res. Public Health 2019, 16, 4665. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16234665

AMA Style

Wang J, Bai X, Liu F, Zhang J, Chen F, Lu Q. Enrichments of Cadmium and Arsenic and Their Effects on the Karst Forest Area. International Journal of Environmental Research and Public Health. 2019; 16(23):4665. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16234665

Chicago/Turabian Style

Wang, Jinfeng, Xiaoyong Bai, Fang Liu, Jian Zhang, Fei Chen, and Qian Lu. 2019. "Enrichments of Cadmium and Arsenic and Their Effects on the Karst Forest Area" International Journal of Environmental Research and Public Health 16, no. 23: 4665. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16234665

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