3.1 Multivariate statistical analysis
3.1.1 Physicochemical properties and heavy metal concentrations.
The average contents of heavy metals in plants were compared with the Hygienical standard for feeds (CNSBQTS, 2017). As shown in the Table 1, the concentrations of heavy metals in herbage were far lower than the standard contents, indicating that the concentrations of heavy metals in herbage in this area were at a safe level, which was suitable for grazing.
The pH of soils ranged from 7.9 to 8.2 with an average value of 8.1 (Table 2). EC in soils was within the range of 165 - 282 μs cm-1 and with an average of 215 μs cm-1. The EC in soil reflects the content of soluble salt directly. High EC value could affect microbial community and herbage growing. CaCO3 contents in soils ranged from 43.39 to 183.45 g kg-1, with a mean value of 68.48 g kg-1. Calcium carbonate in soil was one of the main factors of alkalinity (pH > 7). OM in soils ranged within 48.77 - 160.79 g kg-1 with a mean value of 88.68 g kg-1 which was higher than the background.
The mean values of Hg, As, Cu, Mn, Cd, Zn, Pb, Cr were 40.9 ± 10.8 μg kg-1, 13.2 ± 1.3 mg kg-1, 34.0 ± 2.1 mg kg-1, 876 ± 159 mg kg-1, 0.232 ± 0.057 mg kg-1, 81 ± 16 mg kg-1, 34.8 ± 7.3 mg kg-1, 136.6 ± 31.6 mg kg-1 respectively. The mean values of heavy metals were lower than the Environmental Quality Standard for Soil (CNEPA,1995) except Cd and Cr. All of the mean values of metals in soils were higher than their background values of Gansu Province soil (CNEMC, 1990), suggesting that heavy metal concentrations in the studied soils were affected by human activities, but not leading to a significant contamination of soils. The coefficient variations (CV) of heavy metals were comparatively low, suggesting that the concentrations levels of these heavy metals in soils were less variable, and the extrinsic inputs of the metals into the soils were limited.
3.1.2 Correlation analysis
The results of correlation analysis were shown in Table 3. For soil physicochemical properties, pH showed a positive correlation with Pb, indicating that neutral-alkaline soil was good for the accumulation of lead in soil. EC showed a significant positive correlation with mercury (r = 0.528, p<0.05) but negative with arsenic (r = -0.600, p<0.05). The reason might be that the properties of mercury and arsenic were different. Hg was positively charged and As was negatively charged. It was also revealed that soluble salts might reduce arsenic concentration. OM showed a negative correlation with arsenic (r = -0.546, p < 0.05), contrary to the results of some other researchers (Quenea et al., 2009), which might be caused by regional differences. CaCO3 content correlated positively with Cd (r = 0.656, p < 0.01), indicating that calcium carbonate promotes the accumulation of Cd in soil, which also confirmed that cadmium and calcium may have similar properties. Hg, As, Cu, Mn, Cd, Zn, Pb and Cr were not significantly correlated with each other, implying that the sources of these heavy metals were different.
3.1.3 Relationship between elevation and soil physicochemical properties
Elevation was the most important factor affecting the abiotic environment by changing climatic and topography (Holechek et al., 2010, Roukos et al., 2017). Most subalpine grasslands were on steep slopes, both elevation and degree of slope influenced plant diversity and soil properties (Hadjigeorgiou et al., 2005, Roukos et al., 2011).
Linear regression analysis (Figure 2) showed that the contents of CaCO3 and TP decreased significantly with the increase in elevation, while the contents of soil organic matter increased significantly with the increase in elevation. EC and herbage had a positive correlation with elevation, while pH had a negative correlation with elevation, but this trend was not very obvious.
The decrease in CaCO3 content with the increase of elevation could be due to the increase in soil organic matter and water content (Ali et al., 2019). The difference in soil parent material and small area climate might also be the reasons. The change of pH, temperature and precipitation also affected the change of CaCO3 with elevation (Ali et al., 2017), because acidity, temperature and precipitation affected the solubility of CaCO3.
The change of OM with elevation was mainly related to temperature and plant density. Lower temperature delays decomposition of OM (Charan et al., 2012), because low temperature reduces soil enzyme activity. As the main source of organic matter, the vegetation material also had a positive relationship with vegetation density (Williams et al., 2003). The plant yield in the study area increased with the increase in elevation. However, the content of soil organic matter at low elevation was still higher than that at high elevation in this study, which could be because of grazing in the study area. Grazing was more frequent at higher elevation (more than 3000 m asl) in this area, and as a result, the enrichment efficiency of organic matter in soil was greatly reduced. Total Kjeldahl nitrogen (TKN) had almost no linear correlation with elevation, which was consistent with Cao et al (2020). But there was a significant negative correlation between TP and elevation. This could be caused by the dissolution of soil soluble phosphorus into rainwater and its downward migration along the slope (Roberts and Bettany, 1985). In addition, soil water content, phosphorus adsorption capacity, pH value and the interaction among microclimate, topography and vegetation might also cause the change of soil total phosphorus with elevation (He et al., 2016).
3.1.4 Relationship between elevation and soil trace element
Elevation was one of the factors affecting the distribution of heavy metals in subalpine grassland soil. Affected by elevation, the precipitation increased with the increase in elevation. When the elevation reached a certain level, the precipitation decreased and the total deposition might be changed (Salerno et al., 2015; Reiners et al., 1975). The Figure 3 showed that except for Cu and Cd, the content of other trace elements was linear with elevation.
The concentration of Hg, Zn and Cr increased with elevation, while that of As, Pb, Mn decreased. Hg, in the forms of Hg0, gaseous oxidized Hg, and particulate phase mercury (Grigal, 2002), can undergo long-range atmospheric transport, due to its volatilization and persistence in the environment. Then, Hg entered the soil through precipitation or atmospheric deposition, which might be one of the reasons for the positive correlation between Hg and elevation. The change trend of Zn and Cr were similar to that in Magnani et al. (2018), but Mn and Pb were not the same to some studies (Magnani et al., 2018; McGee and Vallejo, 1996; Reiners et al., 1975). As a matter of fact, in addition to air pollution, topographic conditions and microclimate also affect the distribution of heavy metals (Bergamasch et al., 2002; Cong et al., 2015). Several studies in Alps showed that deposition was increased with elevation (Camarero et al., 2009), but it was not found in other studies (Kang et al., 2007). The non-linear changes of Cu and Cd indicated that these two elements might be affected by many factors, such as human activities, atmospheric deposition, topographic conditions and parent materials. As a metalloid, the concentration of arsenic was affected by many factors, such as pH, speciation (arsenate, arsenite), organic matter, etc (Cai et al., 2002; González et al., 2006). In this study, arsenic concentration and organic matter showed a significant negative correlation, indicating that arsenic was more affected by organic matter concentration, which was different from the results of González et al. (2006).
3.1.5 BCFs in herbage
As shown in Table 2, Hg, Cd and Zn were the most translocated metals. BCFs were in the decreasing order of Hg (0.417) > Cd (0.393) > Zn (0.302) > Cu (0.181) > Mn (0.087) > Pb(0.021)> As (0.019) > Cr (0.004). There were many factors to control the accumulation and bioavailability of heavy metals, including: sequestration and speciation, active/passive transfer processes, redox states and the response of plants to elements in relation to seasonal cycles (Badr et al., 2012). Soil structure and texture also affected the absorption of heavy metals by herbage. Unlike other elements, Zn occurs in the soil frequently in easily soluble forms (Gawryluk et al., 2020) and Zn is usually accumulated the most in the aboveground tissues of plants in an ecosystem where this element occurs in the air (Kabata-Pendias and Pendias, 2001).
All BCFs in this study were less than 1, indicating that heavy metals were not accumulated excessively in plants. It was worth noting that the BCF of mercury, cadmium and zinc in herbage in the study area reached 0.3. Although the concentration of heavy metals in herbage was at a safe level, the concentration of these three metals in soil should be monitored.
3.2 Assessment of ecological risk
3.2.1 Geo-Accumulation Index
The Igeo of Hg, Cd and Cr in most samples reached level 1 to 2 (Table 4), indicating that Hg, Cd and Cr in the soil of the study area were more obviously affected by human activities compared with other metals. And 93% samples reached level 1 of Pb, 33% samples reached level 1 of Mn, and the other samples were level 0, implying that Pb, Mn in the soil slightly affected by human activities. As, Zn and Cu in soil were almost not affected by human beings, because almost all samples were evaluated at level 0.
Qilian Mountain was rich in mineral resources (Qin et al., 2016). According to the survey, there were some mining sites around the study area. Although some of the mining sites were closed at the request of the government, the impact on the soil quality directly or through atmospheric subsidence was long-term. The behavior of herdsmen in the process of grazing, other human activities such as traffic in the scenic area may be the main reason that the metal content in the soil in the study area was higher than the background value.
3.2.2 Potential ecological risk index and Nemerow pollution index
The contamination level of single element was generally low, and most CF values were ranged from 1 to 3 (Table 5), which indicated that they were at moderate contamination and the impact of human activities was limited (Williams and Antoine, 2020). Arsenic, zinc, copper and manganese in most sampling sites reached moderate contamination. The sites with high zinc ecological risk were concentrated in areas with high elevation, which were mainly affected by grazing activities, while the sites with high ecological risk of arsenic were concentrated in the areas with low altitude, and the human impact was relatively less frequent. The higher ecological risk of mercury in high elevation areas might be related to the migration characteristics of mercury. It was well known that mercury could migrate with the atmosphere and enter the soil with precipitation and atmospheric deposition. The higher precipitation in high elevation in study area made the soil more likely to enrich mercury. As mentioned above, plant biomass in the study area increased with elevation, and higher plant density might have a positive effect on mercury interception in precipitation.
In terms of RI (Table 6), 5 sampling sites had RI values lower than 150, indicating low ecological risk, and other sampling points were in the range of 150-300, which implying that the ecological risk was moderate. The maximum RI value was 227.31 (S8 sampling site), which was still far below the upper limit value of moderate risk. The results showed that most of the study area was affected by human activities, but it had no serious impact on the ecological environment.
The calculated PN values of trace elements were presented in Table 5. The data indicated that most of the study area were slightly contaminated by trace elements. It was worth noting that S6 and S12 sampling sites were moderately contaminated. The S6 sampling points was located in the buffer zone of the reserve, which was a pasture in autumn and winter. Therefore, grazing was a factor affecting the environmental quality.
From the analysis results of risk assessment, it could be seen that most of the soils in the study area were at or below the moderate risk. Only a few sampling sites reached high risk for Hg or Cr or Cd, which was due to serious human activities in these sites.
3.3 Chemical speciation of heavy metals
The chemical activity, mobility, bioavailability of heavy metals in the environment and their impact on ecosystem or organism couldn’t be well explained only based on the total amount of elements, therefore, the present study examined the different chemical forms of heavy metals in soil samples through the Tessier extraction procedure aside from determining the total concentration. The results were shown in the Figure 4.
Exchangeable state referred to the part of metal that was not specifically adsorbed on the surface of soil colloid, but also easily absorbed by plant roots. Among the eight metals, the largest proportion of exchangeable state to total value was Pb, and they were in the decreasing order of Pb (20.58%) > Mn (9.15%) > Cd (6.59%) > Cu (3.96%) > Cr (3.49%) > Zn (2.50%) > As (0.05%) > Hg (0.03%). The activity of Pb was obviously higher than that of other metals, indicating that Pb had a strong contribution from anthropogenic source in soils and a high probability transferring from soil to plants and underground water (Ma and Rao, 1997; Kaasalainen and Yli-Halla, 2003). The activity of As and Hg was very low, indicating that the impact of the two heavy metals on the environment was limited. The results were different from the BCFs analyzed above, which might be due to the low content of heavy metals in the soil and the limited absorption of plants, so even though the exchangeable metal content was high, the metal content absorbed by plants was still small due to the low enrichment ability of plants. For example, the exchangeable state content of Pb accounted for 20.58% of the total value, but the BCF value was 0.021, while the BCF value of Hg was the largest (0.471), but the exchangeable state content of Hg only accounted for 0.03% of the total value. It made the quality of herbage in the study area not at risk. However, it was still necessary to detect the exogenous introduction of heavy metals, especially Pb and Hg.
As shown in the Figure 4, in terms of chemical form distribution, the metals studied could be divided into three categories. Hg and As were mainly in the form of residue, organic matter and iron manganese oxide (94.75% and 92.63%, respectively), indicating that the bioavailability and mobility of these two heavy metals were limited. The residue, organic matter and iron manganese oxide of Cu, Zn and Cr accounted for 86.57%, 88.8% and 87.01% of the total value, respectively. These three metals mainly existed in stable form, but the proportion of this form was significantly lower than Hg and As. Pb, Cd and Mn were the most environmentally risky metals, and their exchangeable and carbonate bound states accounted for 35.65%, 44.59% and 25.09% of the total values, respectively, and it confirmed the strong contribution of anthropogenic pollution to their accumulation.