Next Article in Journal
Tree Rings Elucidate Differential Drought Responses in Stands of Three Mexican Pines
Next Article in Special Issue
Organic Mulching Alters the Soil Microclimate, Increases Survival and Growth of Tree Seedlings in Restoration Planting
Previous Article in Journal
Estimating the Vertical Distribution of Biomass in Subtropical Tree Species Using an Integrated Random Forest and Least Squares Machine Learning Mode
Previous Article in Special Issue
Influence of Date Palm-Based Biochar and Compost on Water Retention Properties of Soils with Different Sand Contents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil Quality Variation under Different Land Use Types and Its Driving Factors in Beijing

1
Jianshui Research Station, School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China
2
College of Life Science, Yan’an University, Yan’an 716000, China
3
State Key Laboratory of Efficient Production of Forestry Resources, Beijing Forestry University, Beijing 100083, China
4
Engineering Research Center of Forestry Ecological Engineering, Ministry of Education, Beijing Forestry University, Beijing 100083, China
5
China National Botanical Garden, Beijing 100093, China
6
Academy of Forestry Inventory and Planning, National Forestry and Grassland Administration, Beijing 100714, China
*
Author to whom correspondence should be addressed.
Forests 2024, 15(6), 993; https://doi.org/10.3390/f15060993
Submission received: 9 May 2024 / Revised: 3 June 2024 / Accepted: 4 June 2024 / Published: 6 June 2024

Abstract

:
With the advancement of urbanization, land resources are becoming increasingly strained, particularly for urban greening purposes. In this context, a large number of newly cultivated lands dominated by construction waste and backfill soil are emerging in cities. Assessing the soil quality of these newly cultivated lands and achieving their rational utilization accurately and quantitatively has become an urgent issue. In this study, soil samples of five land use types, namely newly cultivated land (NCL, control), adjacent cropland (CL), arbor–shrub mixed forest (ASF), arbor forest (AF), and shrubland (SL) were selected around Beijing, China. ASF, AF, and SL are also newly cultivated lands composed of construction waste and backfill before greening. Based on principal component analysis (PCA), a total data set (TDS) and a minimum data set (MDS) were used to construct the soil quality index (SQI) model. Soil quality indicators covering the physical and chemical characteristics of the soil and their relationships with land use types were studied with the Partial Least Squares Path Model (PLS-PM). The results were summarized as follows: (1) The soil quality index under different land use types in the Beijing plain area were in the order of arbor–shrub mixed forest (ASF) > arbor forest (AF) > shrubland (SL) > cropland (CL) > newly cultivated land (NCL). (2) Soil organic carbon (SOC), soil water content (SWC), maximum water-holding capacity (MWHC), capillary water-holding capacity (CWHC), Pb, and Cd were identified as the MDS. The MDS of the soil quality assessment model showed a linear relationship with the TDS (y = 0.946x + 0.050, R2 = 0.51). (3) Land use types have an indirect impact on soil quality by changing the content of Pb. The chemical indicators’ coefficient (0.602) contributed more to the SQI than did the physical indicators’ (0.259) and heavy metal elements’ (−0.234). In general, afforestation and agricultural production could improve the newly cultivated lands’ soil quality, but afforestation is much better than agricultural production. These results will help to evaluate the SQI in the Beijing plain area objectively and accurately, and they have significant implications for soil restoration and management.

1. Introduction

Soil is an important carrier and the major source of nutrients for plant growth. Soil is also a critical factor to support ecosystem biodiversity, structure, and function [1]. Soil quality is defined as the capacity of soil to sustain plant and animal productivity and maintain human and environment health in natural or human ecosystems [2,3]. Reasonable land use types have been considered as an effective method to improve soil structure, contributing to ensure crop yields and improve environmental quality [4]. Nevertheless, unreasonable land use types will degrade the soil and decrease the soil quality, lead to slow vegetation growth, and finally even vegetation death [5]. It is difficult to evaluate soil quality accurately. Numerous soil quality evaluation methods have been developed to comprehensively assess soil quality, such as the matter–element method and the grey correlation analysis method. Yu et al. evaluated the land eco-security of nine cities with the theory of entropy weight and the matter–element model [6]. Tang et al. evaluated the effects of land use type on soil characteristics by using grey correlation analysis [7]. The matter–element method has a relatively complicated calculation process, and the grey correlation analysis method is prone to lose information with discrete values of indicators. However, the soil quality index (SQI) is an effective approach which integrates multiple soil properties into a comprehensive index and is widely used to assess the impact of natural and anthropic factors on soil quality because of its flexibility and simplicity [8,9]. There are numerous factors that affect soil quality; thus, when evaluating soil quality, it is necessary to identify indicators that can truly and effectively reflect the soil conditions. [10]. Principal component analysis (PCA), reducing the complexity of the total data sets by reducing redundancy and improving the reliability and usefulness of the results, has become the most used method for selecting soil quality evaluation indicators.
With the expansion of cities and the increase in population, urban land resources are becoming increasingly strained, leading to more severe ecological security issues. Since 2012, in the plain area of Beijing, a lot of newly cultivated land containing construction waste and backfill has appeared in the urban and suburban areas, and afforestation projects have been implemented on this area, forming different land use types. These newly cultivated lands’ soil has problems, such as unreasonable structure, nutrient deficiency, and heavy metal pollution [11]. Different land use types have a profound impact on soil and vegetation. Therefore, how to evaluate soil quality and rationally utilize it is an urgent problem to be solved at present. In most studies of soil quality assessment, the characteristics of soil quality have been emphasized, with attention being given to soil physical indicators, nutrient indicators, and microbial indicators [12]. It is reported that construction waste may cause heavy metal pollution in soil. Heavy metal pollution of urban soil has been a hot topic in recent research [13]. Wang et al. [14] stated that construction waste pollution in Beijing was at a moderate level, and the main pollution factors were Pb and Cd. Soil heavy metals, which are persistent, difficult to degrade, and easy to accumulate, will lead to the deterioration of soil environmental quality. In this paper, we choose soil samples of five different land use types in Fangshan District, Beijing for study, including newly cultivated land (NCL, control), adjacent cropland (CL), arbor–shrub mixed forest (ASF), arbor forest (AF), and shrubland (SL). ASF, AF, and SL are also newly cultivated lands composed of construction waste and backfill before greening. We measured 15 indicators, including physical, nutrient, and heavy metal elements. Principal component analysis was used to screen soil indexes, and the nonlinear scoring method was used to calculate soil quality. The main objectives of this study were to (1) compare the differences in soil properties among different land use types, (2) screen the soil quality evaluation indicators and establish a soil quality evaluation model based on the minimum data set to evaluate the soil quality of different land use types, and (3) analyze the relationship between soil quality, soil indicators, and land use types. The results are expected to provide a scientific basis and data support for quantitative land assessment and rational planning and utilization in the Beijing Plain and similar areas.

2. Materials and Methods

2.1. Study Area

The study was conducted at the southeast of Fangshan District (115°25′~116°15′ E and 39°30′~39°55′ N), Beijing, which is dominated by plains with an elevation of approximately 20–60 m (Figure 1). The region belongs to a continental monsoon climate with a mean annual temperature of 10–12 °C and a mean annual frost-free period of 202 days. The climate is characterized by four distinct seasons (spring, summer, autumn, and winter). The mean annual rainfall is approximately 587 mm, where about 85% falls from June to August. Fangshan District is rich in plant resources, with natural vegetation dominated by temperate deciduous broad-leaved forests, followed by temperate coniferous forests. The dominant greening tree species include Populus L., Salix babylonica L., Ulmus pumila L., Styphnolobium japonicum (L.) Schott, and Larix gmelinii (Rupr.) Kuzen. Platycladus orientalis (L.) Franco, and the main shrubs are Vitex negundo var. heterophylla (Franch.) Rehd. and Ziziphus jujuba var. spinosa (Bunge) Hu ex H.F.Chow. The soil type is classified as Anthrosols [15], and the average soil thickness is 60 cm.

2.2. Soil Sampling and Laboratory Analyses

Through comprehensive field research and investigation in the afforestation area of the Beijing Plain, we selected five different land use types in the plain area of Fangshan District, including newly cultivated land mainly composed of construction waste and backfill, as well as arbor–shrub mixed forest, arbor forest, and shrubland planted on the newly cultivated land, as well as nearby cropland (Table 1). An area of 20 m × 20 m was selected on each field plot for investigation and sampling. At each sampling spot, we excavated three soil profiles with a depth of 60 cm, which were divided into four layers vertically from the surface, namely 0–10 cm, 10–20 cm, 20–40 cm, and 40–60 cm. Soil samples were collected using the cutting-ring method. Three ring knife samples were taken from each layer for replication to determine the physical properties of the soil. Meanwhile, 500 g mixed soil samples were taken from the same soil layer, packed into the cloth bags, and taken back to the laboratory for wind drying to test the chemical characteristics and heavy metals.
The detailed methods for measuring and analyzing soil indicators are shown in Table 2 [16,17,18].

2.3. Soil Quality Evaluation Methods

This study used the SQI method to evaluate the soil quality of different land use types in the Beijing plain afforestation area. The total data set and the minimum data set are used to identify appropriate indicators. One-way analysis of variance (ANOVA) was performed on the 15 indicators to assess the influence of different land uses on soil indicators. Only the indicators showing significant treatment differences (p < 0.05) were chosen as members of the total dataset [4]. The indicators of the minimum data set are screened out through principal component analysis of the standardized data matrix of the total dataset [4,19]. The calculation process of soil quality includes the following four steps: (1) identification of the minimum data set, (2) calculation of membership degrees, (3) determination of indicator weights, and (4) calculation of the SQI [19,20,21].

2.3.1. Identification of the Minimum Data Set (MDS)

To reduce redundancy and information overlap, this study established a minimum dataset to screen out the most representative indicators. Firstly, the loadings of each soil indicator on the principal components (PCs) with eigenvalues ≥ 1 are calculated, and the explained total variance is greater than 85%. For each PC, the indicators with loadings within 10% of the maximum weighted loading are retained to represent that PC. If multiple indicators are retained for a PC, the correlation between the indicators needs to be considered to determine whether they should be included in the minimum dataset. If the remaining indicators are positively correlated, the indicator with the highest correlation coefficient with the rest of the indicators is selected, otherwise, all indicators in that group are included in the minimum dataset [22].
The calculation formula for the norm value is as follows:
N i k = i k ( U i k 2 λ k )
where Nik is the comprehensive loading of the i-th variable on the first k principal components with eigenvalues ≥ 1; Uik is the loading of the i-th variable on the k-th principal component; and λk is the eigenvalue of the k-th principal component.

2.3.2. Calculation of Membership Degrees

The data are converted to values between 0 and 1 using a nonlinear scoring function method, where 1 represents a high level of the indicator and 0 represents a low level of the indicator.
According to the contribution of indicators to soil quality, they are divided into “more is better” and “less is better” types. The formula is as follows:
F ( x ) = 1 1 + ( x / x m a x ) b
where F(x) is the nonlinear score of the indicator, x is the observed value of the indicator, and xmax is the average observed value of the indicator. b is the slope of the equation. For “more is better” type curves, b is −2.5, and for “less is better” type curves, b is 2.5.

2.3.3. Determination of Indicator Weights

In order to objectively evaluate soil quality, the weights of each soil indicator in this study are determined using PCA. First, the indicators are analyzed using PCA to obtain the common factor variance of each indicator. Then, the weight of each indicator is calculated by dividing its common factor variance by the sum of common factor variances of all indicators.

2.3.4. Calculation of the SQI

The calculation formula for the SQI is as follows:
S Q I = i = 1 n R i × F ( x i )
where Ri is the weight of the soil indicator, F(xi) is the membership degree of the soil indicator, and n is the number of soil indicators.

2.4. PLS-PM Method

The PLS-PM includes the measurement model and the structural models. The measurement model is used to define the latent variables. The structural model is used to reflect the relationship between the latent variables [23,24]. The formulas are as follows:
x = Λxξ + δ
y = Λyη + ε
where x represents exogenous indicators and y represents endogenous indicators. Λx and Λy are the matrix of factor loadings of ξ to x and η to y. δ and ε refer to the measurement error. The formula is as follows:
η = + Γξ + ζ
where η is the endogenous latent variable, ξ is the exogenous latent variable, B and Γ refer to the regression path coefficient of the effect between different η and the effect of ξ on η. ζ represents the regression residuals.

2.5. Statistical Analysis

The data were processed and calculated using Excel 2010. The statistical analysis of the data was performed using SPSS 26.0, including one-way analysis of variance and principal component analysis. The driving factors of soil quality differences were analyzed using SmartPLS 4 for PLS-PM analysis. The figures were plotted using Origin 2018.

3. Results

3.1. Changes in Measured Soil Indicators

Significant differences were found for all soil properties, except for noncapillary porosity (NCP). A total of 14 soil indicators were chosen as members of the total data set. Then, we used PCA to reduce redundant indexes of the SQI calculation, including four chemical indicators, five physical indicators and five trace element indicators.
As shown in Table 3, the NCL had a higher pH (8.04) and lower SOC (0.45 g∙kg−1), TN (0.05 g∙kg−1), and TP (0.18 g∙kg−1) compared with those of the other four land use types, and showed a significant difference in TN. ASF was characterized by higher SOC (1.79 g∙kg−1) and TN (0.23 g∙kg−1) compared with those of the other land use types, and showed a significant difference in SOC, except for CL. For TP, CL was the highest (1.52) and had a significant difference with other land use types.
In terms of physical properties, NCL has the largest MWHC and CWHC among the five land use types, which are 38.98% and 31.29%, respectively, while its soil BD is the smallest (1.05 g∙cm−3). The SWC of ASF (15.69%) is significantly higher than that of CL (11.67%), AF (8.46%), and SL (8.16%), but there is no significant difference between ASF and NCL (14.93%). AF has the largest CP (36.55%) without a significant difference from SL (36.01%) and ASF (33.63%). There is no significant difference in NCP among the five land use types.
ASF has significantly higher Cu (25.77 mg∙kg−1) and Fe (26.70 mg∙kg−1) contents than NCL, AF, and SL, but there is no significant difference between ASF and CL. CL has significantly higher Cd (0.05 mg∙kg−1) and Pb (19.53 mg∙kg−1) contents than the other four land use types, which are approximately five times and 1.2–1.8 times higher, respectively. CL’s Zn (106.28 mg∙kg−1) content is about 1.7 times higher than AF and SL, but there is no significant difference between CL and NCL (94.32 mg∙kg−1) or ASF (104.73 mg∙kg−1).

3.2. Selecting MDS Indicators

As shown in Table 4, the first four PCs had eigenvalues > 1.0 and explained >86.65% of the variance of the original data (Table 4). The first PC explained 36.39% of the total variance. The highest loading value is Pb. Cu had a loading value within 10% of the highest loading value. Pb and Cu were significantly (p < 0.01) correlated with each other (Table 5). Therefore, only Pb was selected as the indicator of PC1. The second PC explained 26.15% of the total variance and met the selection requirements in four indicators, which were SOC, N, MWHC, and CWHC. SOC had the highest loading value. SOC and TN were significantly (p < 0.01) correlated with each other, but TN was not chosen. However, the correlations among MWHC, CWHC, and SOC were insignificant, these two indicators were kept in the MDS as the indicator of PC2. The third PC explained 15.62% of the total variance. Only SWC was retained as the indicator of PC3. Similarly, only Cd was reserved as the indicator of PC4. The final indicators of SQI determined by PCA are Pb, SOC, MWHC, CWHC, SWC, and Cd.

3.3. Soil Quality Index under Different Land Uses

We used the nonlinear scoring methods to transform the soil indicators in the TDS and MDS. The parameters of nonlinear equations, and the weights for the selected indicators, are shown in Table 6. Finally, the comparative SQI can be described as follows:
SQI–NLT = 0.072 × C + 0.067 × pH + 0.072 × N + 0.080 × P + 0.071 × SWC + 0.077 × BD + 0.077 × MWHC + 0.052 × CP + 0.073 × CWHC + 0.080 × Cd + 0.076 × Cu + 0.065 × Fe + 0.075 × Pb + 0.063 × Zn
SQI–NLM = 0.170 × C + 0.168 × SWC + 0.159 × MWHC + 0.166 × CWHC + 0.169 × Cd + 0.168 × Pb
As shown in Table 7, the differences in soil quality among different land use types calculated based on the TDS and the MDS are similar. The order of SQI values from high to low for the five land use types is as follows: ASF > AF > SL > CL > NRL. The SQI values of the other four land use types are both higher than that of NCL, indicating that compared to NCL, both afforestation and agricultural production can improve soil quality.
As shown in Figure 2, the linear fitting of the MDA and the TDS gives a fitted equation of y = 0.946x + 0.050 with R2 = 0.51, indicating a significant positive correlation between the two methods.

3.4. Analysis of Driving Factors of Soil Quality under Different Land Uses

To illustrate the relationship between soil quality and soil physicochemical factors and heavy metal elements under different land use types, a Partial Least Squares Path Model (PLS-PM) was constructed (Figure 3). Chemical indexes (SOC, TP), physical indexes (CP, MWHC), and Pb were the key factors contributing to changes in soil quality under different land use types. The meanings of parameters in the model are shown in Table 8 [25,26,27]. The results of the VIF test, the reliability and validity of the model, and the loadings for each of the manifest factors are shown in Table 9. The path validity tests of the model are shown in Figure 3.
In the PM, the path coefficients of SQI with different land use types, chemical indexes, physical indexes, and Pb were −0.18, 0.60, 0.26, and −0.23, respectively.

4. Discussion

4.1. Soil Indicators under Different Land Use Types

Vegetation has a profound influence on soil physicochemical properties [28]. In this study, the soil chemical properties of NCL, which served as a control, were obviously lower than the other four land use types. These results confirmed that vegetation affects soil parameters significantly, which is consistent with the conclusions of a previous study by Wu et al. [29]. The highest SOC and TN contents were found in ASF, which indicated that ASF had a more significant effect on improving soil chemical properties. The cover on land or incorporation of plant litter within the top-soil layer is a major source of soil organic matter. Litter is transformed into humus through microbial decomposition and a slow decomposition process, which can be stably preserved in soil for a long time to maintain carbon storage [30]. In addition, less tillage can reduce the decomposition rate of organic carbon, thereby increasing the soil organic carbon content [31]. In this study, the ASF accumulated more soil organic carbon because of its larger vegetation coverage, higher stand density, and greater amount of litter. Some studies have shown that deciduous broad-leaved tree species have a better effect on soil organic carbon accumulation than conifer species [32]. Research has shown that compared to monoculture forests, mixed forests significantly increase the total organic carbon and active organic carbon content of forest soil [33]. In this study, the ASF mainly consists of deciduous broad-leaved tree species, such as Populus tomentosa and Vitex negundo, which increases the content of soil organic carbon. The TP content of CL is significantly higher than the other four land use types, which is induced by artificial application of P fertilizer. Studies have shown that crops can absorb a small fraction of P. Therefore, long-term fertilization leads to the accumulation of P in CL soil.
Soil physical properties also showed significant differences among different land use types. Soil bulk density and soil porosity can reflect the compactness of soil and can be affected by tillage methods and soil microstructure [34,35]. Most previous studies on soil bulk density under different land use patterns have shown that newly cultivated land has a higher bulk density than forest and grassland [36,37]. However, in this study, the bulk density of NCL is the lowest, where the primary reason for this observation is the lack of human intervention in the land use of NCL compared to other sites. The foreign soil changes the original microstructure of the soil, affecting the aggregate structure of the soil, and thereby reducing the soil bulk density. Higher BD in the cropland can be attributed to frequent cultivation and mechanical compaction, which caused a reduction in soil porosity and deterioration of the soil structure. Soil water-holding capacity is affected by multiple factors, such as soil bulk density, soil porosity, and land use type [38]. In the current study, the water-holding capacity of soil is directly proportional to porosity and inversely proportional to bulk density. Different land use types also affect SWC. There are differences in surface vegetation types and root distribution in soil under different land use types, which in turn affect the process of soil water infiltration and circulation [39]. The highest soil water content in ASF was due to the good soil structure, which was conducive to water content. In addition, the high vegetation coverage has a shading effect, reducing the evaporation of soil moisture. We found that NCL also has a high SWC, which is consistent with the results of Wang et al. [37]. This may be due to the fact that there is little vegetation in NCL, and therefore no vegetation consumption of soil water in NCL.
Previous studies have demonstrated that livestock manure serves as the primary source of trace elements in agricultural soil, particularly Cu and Zn [40,41]. Conversely, mineral fertilizer predominantly contributes to Cd contamination, with phosphate fertilizers generally exhibiting high levels of Cd content [42]. In this study, the higher concentrations of Cu and Zn in CL may be attributed to the application of livestock manure. The TP in CL is significantly higher compared to the others, indicating a greater usage of P fertilizer, which consequently leads to a significantly elevated Cd content. The fertilization also influences soil pH and cation exchange capacity, subsequently affecting the levels of heavy metal elements content [43]. Despite the presence of construction waste and backfill soil in other land use types, variations in plant root absorption and human activities contribute to differences in heavy metal content. For instance, plant extraction is an economical and environmentally friendly method that utilizes special plants to remove heavy metals from contaminated soil, thereby reducing their concentration to a safe level [44]. Additionally, human activities affect both the input and migration redistribution of heavy metals in soil [45].

4.2. Soil Quality Indexes and Effects of Land Use Types on Soil Quality

The SQI has emerged as a robust and widely adopted approach for assessing soil quality in numerous scientific investigations. TDS and MDS, two commonly used methods for calculating SQI, each have their own advantages and limitations. The TDS can provide a more comprehensive set of soil indicators, while the MDS is simpler and more cost-effective. In this study, six indicators (Pb, SOC, MWHC, CWHC, SWC, and Cd) with highly weighted factors were selected in the MDS for evaluation. The indicators of SOC, MWHC, CWHC, and SWC have been widely adopted as the minimum data set in numerous previous studies [4,22,46,47]. SOC serves as a pivotal indicator for assessing soil quality, exerting positive influences on various physical, chemical, and biological indicators of soil. For instance, the presence of organic carbon enhances soil structure, augments soil fertility, and stimulates the activities of soil fauna and microorganisms [22]. The SWC, MWHC, and CWHC indexes are all closely associated with soil moisture, which is a vital determinant of plant and microorganism survival, directly impacting their ecological conditions. Given the water scarcity in the arid and semi-arid study area, the accurate assessment of soil quality necessitates a comprehensive evaluation of soil moisture indicators. Previous studies on soil quality assessment have rarely included heavy metal indicators; however, given the presence of construction waste in our study site, it is imperative to acknowledge the significant impact of heavy metals associated with such waste [13,14]. Through the principal component analysis and minimum data set screening, Cd and Pb were identified as indicators for evaluating soil quality in this study. The concentrations of Cd and Pb in soil are influenced by both natural geological background and human activities [48]. Wang et al. [14] conducted an analysis on the pollution characteristics of heavy metals in construction waste in Beijing, revealing mild and severe contamination of Pb and Cd, respectively. Consistent with the minimum data set indicator selected in this study, Joimel et al. [49] proposed using Pb as an indicator to assess soil quality. The presence of Cd in agricultural soil in China is highly concerning, with an average annual input of 0.004 mg/kg/a [50]. Furthermore, Cd exhibits pronounced migratory capabilities, facile absorption by plants, and toxicity towards both plants and soil microorganisms. After Cd enters the human body through the food chain, it may lead to deadly diseases such as cancer [11]. Consequently, the higher the Cd concentration in soil, the greater ecological risk it poses.
The soil quality varies among different land use types in this study (Table 6). The SQI results for various land use types indicate that ASF > AF > SL > CL > NCL. Meng et al. [51] have indicated that the growth of vegetation and the application of fertilizers can enhance the physical and chemical properties of soil, thereby improving its quality, which is in line with our findings. In this study, the Vitex negundo in the ASF are characterized as shallow-rooted plants. Previous research has demonstrated that mature trees predominantly utilize deep soil moisture, whereas shrubs primarily absorb surface soil moisture [52]. Consequently, there appears to be minimal water competition among different plant species within the ASF. Additionally, the high vegetation coverage and diverse plant species in ASF contribute to a substantial accumulation of litter, leading to significant nutrient cycling. Consequently, ASF exhibits elevated levels of soil organic matter and nitrogen content, thereby enhancing soil fertility. As a result, the soil quality improvement of ASF is more pronounced compared to AF. However, due to the relatively shallow roots, single plant species, and small biomass of the SL, its soil quality improvement effect is comparatively lower than ASF and AF. The lower soil quality of CL compared to forestland may be attributed to long-term cultivation, which damages the soil structure and mechanically disrupts soil aggregates, thereby accelerating the mineralization rate of organic matter. Additionally, farmers’ harvesting practices also contribute to a decline in soil nutrient content [53]. Fertilization induces alterations in heavy metal content and soil pollution, resulting in lower soil quality on CL compared to forestland. The lowest soil quality observed on NCL suggests that both afforestation and agricultural practices contribute to improving the soil quality, with afforestation exhibiting a superior effect.

4.3. The Dominant Factor Analysis of Soil Quality Change

This study employed PLS-PM to investigate the impact of variables on soil SQI. The findings revealed that variables exerted both direct and indirect effects on soil SQI by influencing other variables. For instance, in this study, the soil chemical indexes (SOC and TP), physical indexes (MWHC and CP), and heavy metal element (Pb) were found to directly influence SQI, with corresponding influence coefficients of 0.602, −0.234, and 0.259, respectively. Moreover, chemical indexes exhibited highly significant influences while physical indexes and heavy metal elements showed significant influences on SQI. However, the direct impact of land use types on SQI was not found to be statistically significant. The land use types, with a path coefficient of −0.752, exert a significant influence on Pb levels, which aligns with the findings reported by Zheng et al. [54] in Beijing. Hence, it is evident that land use types primarily influence SQI through alterations in soil Pb levels, subsequently exerting an indirect impact on SQI. Lead is a highly toxic chemical element that can persist in the environment and has detrimental effects on both animals and humans. It predominantly occurs in soil, where it cannot be biodegraded, posing a significant threat to plant health and human well-being. Lead primarily originates from atmospheric and rainwater deposition, the introduction of lead-containing waste, sludge, and organic fertilizers, as well as the application of pesticides containing lead. The elevated lead content observed in CL may be attributed to the application of chemical fertilizers and pesticides. These findings align with results reported by Li et al. [55], demonstrating that land use types play a significant role in the accumulation and impact of lead in cropland.
In this study, the influence coefficient of chemical indexes on SQI (0.602) surpassed that of physical indexes (0.259) and heavy metal elements (−0.234), aligning with the findings by Zou et al. [56]. This suggests the crucial role of soil chemistry indicators in soil quality restoration. Previous studies have demonstrated that agricultural and forestry practices can effectively enhance soil organic matter content, while the application of fertilizers and improvements in soil structure contribute to increased levels of SOC and TN, thereby augmenting the capacity for carbon and nitrogen sequestration in soils [57]. Furthermore, the augmentation of SOC and TN content can enhance soil structure and improve soil permeability and drainage, thereby mitigating nitrogen loss in the soil [58].

5. Conclusions

Except for NCP, there are differences in soil physical and chemical indicators and heavy metal content among different land use types in the study area. The minimum data set comprised six indicators, including one soil chemical indicator (SOC), three soil physical indicators (MWHC, CWHC, and SWC), and two heavy metal indicators (Cd and Pb). The SQI values of the TDS and the MDS followed the order of ASF > AF > SL > CL > NCL. Afforestation and agricultural practices exhibited improvements in soil quality, with a more pronounced effect observed for afforestation. The land use type has a significant impact on the soil Pb content. Therefore, in the evaluation of urban soil quality, attention should be paid to the heavy metal elements. The future land use types in the study area should be determined according to the needs of production practices, selecting the most rational utilization types. If the primary objective is solely to improve soil quality, then the ASF should be prioritized.

Author Contributions

Conceptualization, J.Z. (Jinxing Zhou); Methodology, F.Q.; Formal analysis, F.Q.; Investigation, C.S.; Resources, F.Q.; Writing—original draft, F.Q.; Writing—review & editing, F.Q.; Supervision, J.Z. (Jinxing Zhou); Project administration, J.Z. (Jiaqi Zhang) and J.Z. (Jinxing Zhou); Funding acquisition, L.J. and J.Z. (Jinxing Zhou). All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the self-developed project of the Academy of Forestry Inventory and Planning, National Forestry and Grassland Administration (LC-6-01), 5·5 Engineering Research & Innovation Team Project of Beijing Forestry University (BLRC2023B09) and Ningxia Key Research and Development Plan Project (2021BEG02005).

Data Availability Statement

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

Acknowledgments

The authors greatly appreciate the assistance of C.S., Qian Dong, Yapei Lu, Gefei Tan, Jianmei Fang, and Liwei Jiang.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Romeo, F.; Settineri, G.; Sidari, M.; Mallamaci, C.; Muscolo, A. Responses of soil quality indicators to innovative and traditional thinning in a beech (Fagus sylvatica L.) forest. For. Ecol. Manag. 2020, 465, 118106. [Google Scholar] [CrossRef]
  2. Karlen, D.L.; Mausbach, M.J.; Doran, J.W.; Cline, R.G.; Harris, R.F.; Schuman, G.E. Soil Quality: A Concept, definition, and framework for evaluation (A guest editorial). Soil Sci. Soc. Am. J. 1997, 61, 4–10. [Google Scholar] [CrossRef]
  3. Guo, L.; Sun, Z.; Ouyang, Z.; Han, D.; Li, F. A comparison of soil quality evaluation methods for Fluvisol along the lower Yellow River. Catena 2017, 152, 135–143. [Google Scholar] [CrossRef]
  4. Raiesi, F. A minimum data set and soil quality index to quantify the effect of land use conversion on soil quality and degradation in native rangelands of upland arid and semiarid regions. Ecol. Indic. 2017, 75, 307–320. [Google Scholar] [CrossRef]
  5. Qiao, X.; Li, Z.; Lin, J.; Wang, H.; Zheng, S.; Yang, S. Assessing current and future soil erosion under changing land use based on InVEST and FLUS models in the Yihe River Basin, North China. Int. Soil Water Conserv. Res. 2023, 12, 298–312. [Google Scholar] [CrossRef]
  6. Yu, J.; Fang, L.; Cang, D.; Zhu, L.; Bian, Z. Evaluation of land eco-security in Wanjiang district base on entropy weight and matter element model. Trans. Chin. Soc. Agric. Eng. 2012, 28, 260–266. [Google Scholar]
  7. Tang, B.; He, B.; Yan, J. Gray correlation analysis of the impact of land use type on soil physical and chemical properties in the hilly area of central Sichuan, China. Chin. J. Appl. Ecol. 2016, 27, 1445–1452. [Google Scholar]
  8. Nortcliff, S. Standardisation of soil quality attributes. Agric. Ecosyst. Environ. 2002, 88, 161–168. [Google Scholar] [CrossRef]
  9. Qiu, X.; Peng, D.; Wang, H.; Wang, Z.; Cheng, S. Minimum data set for evaluation of stand density effects on soil quality in Larix principis-rupprechtii plantations in North China. Ecol. Indic. 2019, 103, 236–247. [Google Scholar] [CrossRef]
  10. Gao, R.; Ai, N.; Liu, G.; Liu, C.; Zhang, Z. Soil C:N:P stoichiometric characteristics and soil quality evaluation under different restoration modes in the loess region of northern Shaanxi Province. Forests 2022, 13, 913. [Google Scholar] [CrossRef]
  11. Yang, J.; Zhang, G. Formation, characteristics, and eco-environmental implications of urban soils—A review. Soil Sci. Plant Nutr. 2015, 61, 30–46. [Google Scholar] [CrossRef]
  12. Maurya, S.; Abraham, J.S.; Somasundaram, S.; Toteja, R.; Gupta, R.; Makhija, S. Indicators for assessment of soil quality: A mini-review. Environ. Monit. Assess. 2020, 192, 604. [Google Scholar] [CrossRef] [PubMed]
  13. Ye, J.; Liu, C.; Zhao, Z.; Li, Y.; Yu, S. Heavy metals in plants and substrate from simulated extensive green roofs. Ecol. Eng. 2013, 55, 29–34. [Google Scholar] [CrossRef]
  14. Wang, H.; Luo, H.; Zhang, W.; Cai, B.; Yang, M.; Qin, T.; Yang, Q.; Jing, Z.; He, B. Research progress on heavy metal leaching characteristics and environmental pollution risk of construction waste. Appl. Chem. Ind. 2023, 52, 2408–2413. [Google Scholar]
  15. lUSS Working Group WRB. World Reference Base for Soil Resource 2014. International Soil Classification System for Naming Soil and Creating Legends for Soil Maps; Food and Agriculture Organization of the United Nations: Rome, Italy, 2014; p. 154. [Google Scholar]
  16. Institute of Soil Science Chinese Academy of Sciences. Soil Physical Properties Determination Method; Science Press: Beijing, China, 1978; pp. 2,4,6. [Google Scholar]
  17. Bao, S. Soil Agrochemical Analysis; China Agriculture Press: Beijing, China, 2000; pp. 2–10. [Google Scholar]
  18. United States Environmental Protection Agency (U.S. EPA). Acid Digestion of Sediments, Sludges, and Soils; Method 3050B.; United States Environmental Protection Agency (U.S. EPA): Washington, DC, USA, 1996. [Google Scholar]
  19. Andrews, S.S.; Mitchell, J.P.; Mancinelli, R.; Karlen, D.L.; Hartz, T.K.; Horwath, W.R.; Munk, D.S. On-farm assessment of soil quality in California’s Central Valley. Agron. J. 2002, 94, 12–23. [Google Scholar]
  20. Ye, C.; Cheng, X.; Zhang, Q. Recovery approach affects soil quality in the water level fluctuation zone of the Three Gorges Reservoir, China: Implications for revegetation. Environ. Sci. Pollut. Res. 2014, 21, 2018–2031. [Google Scholar] [CrossRef] [PubMed]
  21. Yu, P.; Liu, S.; Zhang, L.; Li, Q.; Zhou, D. Selecting the minimum data set and quantitative soil quality indexing of alkaline soils under different land uses in northeastern China. Sci. Total Environ. 2018, 616, 564–571. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, Z.; Ai, N.; Liu, G.; Liu, C.; Qiang, F. Soil quality evaluation of various microtopography types at different restoration modes in the loess area of Northern Shaanxi. Catena 2021, 207, 105633. [Google Scholar] [CrossRef]
  23. Hair, J.F.; Ringle, C.M.; Sarstedt, M. Partial least squares: The better approach to structural equation modeling? Long Range Plan. 2012, 45, 312–319. [Google Scholar] [CrossRef]
  24. Sarstedt, M.; Hair, J.F.; Ringle, C.M.; Thiele, K.O.; Gudergan, S.P. Estimation issues with PLS and CBSEM: Where the bias lies! J. Bus. Res. 2016, 69, 3998–4010. [Google Scholar] [CrossRef]
  25. Sarstedt, M.; Hair, J.F., Jr.; Cheah, J.H.; Becker, J.M.; Ringle, C.M. How to specify, estimate, and validate higher-order constructs in PLS-SEM. Australas. Mark. J. 2019, 27, 197–211. [Google Scholar] [CrossRef]
  26. Turner, G. Organisation in the school dental service. Br. Dent. J. 1971, 131, 33–36. [Google Scholar] [CrossRef] [PubMed]
  27. dos Santos, P.M.; Cirillo, M.Â. Construction of the average variance extracted index for construct validation in structural equation models with adaptive regressions. Commun. Stat.-Simul. Comput. 2023, 52, 1639–1650. [Google Scholar] [CrossRef]
  28. Liu, M.; Chang, Q.; Qi, Y.; Liu, J.; Chen, T. Aggregation and soil organic carbon fractions under different land uses on the tableland of the Loess Plateau of China. Catena 2014, 115, 19–28. [Google Scholar] [CrossRef]
  29. Wu, W.; Chen, G.; Meng, T.; Li, C.; Feng, H.; Si, B.; Siddique, K.H. Effect of different vegetation restoration on soil properties in the semi-arid Loess Plateau of China. Catena 2023, 220, 106630. [Google Scholar] [CrossRef]
  30. Nave, L.E.; Swanston, C.W.; Mishra, U.; Nadelhoffer, K.J. Afforestation effects on soil carbon storage in the United States: A synthesis. Soil Sci. Soc. Am. J. 2013, 77, 1035–1047. [Google Scholar] [CrossRef]
  31. Gao, L.; Wang, B.; Li, S.; Han, Y.; Zhang, X.; Gong, D.; Degré, A. Effects of different long-term tillage systems on the composition of organic matter by 13C CP/TOSS NMR in physical fractions in the Loess Plateau of China. Soil Tillage Res. 2019, 194, 104321. [Google Scholar] [CrossRef]
  32. Du, Z.; Dong, H.; Jing, D.; Ma, B.; Liu, F. Effects of long-term plantations on soil organic carbon pool in Yellow River delta. Bull. Soil Water Conserv. 2016, 36, 056–061. [Google Scholar]
  33. Paul, K.I.; Polglase, P.J.; Nyakuengama, J.G.; Khanna, P.K. Change in soil carbon following afforestation. For. Ecol. Manag. 2002, 168, 241–257. [Google Scholar] [CrossRef]
  34. Guan, J.; Chen, S.; Li, S.; Zhang, Y.; Zhang, X.; Lu, Y.; Yan, Z. Soil tillage practices affecting the soil characteristics and yield of winter wheat and summer maize in North China. Chin. J. Eco-Agric. 2019, 27, 1663–1672. [Google Scholar]
  35. Huang, Y.; Li, Y.; Liu, Y. Effects of soil-layer compounding schemes on the soil fertility of newly-constructed cultivated land. Trans. Chin. Soc. Agric. Eng. 2021, 37, 64–72. [Google Scholar]
  36. Wu, X.; Wang, R.; Gao, C.; Gao, S.; Du, L.; Asif, K.; Sheng, L. Variations of soil properties effect on microbial community structure and functional structure under land uses. Acta Ecol. Sin. 2021, 41, 7989–8002. [Google Scholar]
  37. Wang, M.; Huang, L.; Chen, C. Difference in soil water holding capacity and the influencing factors under different land use types in the alpine region of Tibet, China. J. Appl. Ecol. 2022, 33, 3287–3293. [Google Scholar]
  38. Yang, F.; Zhang, G.; Yang, J.; Li, D.; Zhao, Y.; Liu, F.; Yang, F. Organic matter controls of soilwater retention in an alpine grassland and its significance for hydrological processes. J. Hydrol. 2014, 519, 3086–3093. [Google Scholar] [CrossRef]
  39. Kan, X.; Cheng, J.; Hu, X.; Zhu, F.; Li, M. Effects of Grass and Forests and the Infiltration Amount on Preferential Flow in Karst Regions of China. Water 2019, 11, 1634. [Google Scholar] [CrossRef]
  40. Belon, E.; Boisson, M.; Deportes, I.Z.; Eglin, T.K.; Feix, I.; Bispo, A.O.; Guellier, C.R. An inventory of trace elements inputs to French agricultural soils. Sci. Total Environ. 2012, 439, 87–95. [Google Scholar] [CrossRef] [PubMed]
  41. Santorufo, L.; Memoli, V.; Panico, S.C.; Esposito, F.; Vitale, L.; Di Natale, G.; Maisto, G. Impact of Anthropic activities on soil quality under different land uses. Int. J. Environ. Res. Public Health 2021, 18, 8423. [Google Scholar] [CrossRef] [PubMed]
  42. Nziguheba, G.; Smolders, E. Inputs of trace elements in agricultural soils via phosphate fertilizers in European countries. Sci. Total Environ. 2008, 390, 53–57. [Google Scholar] [CrossRef] [PubMed]
  43. Ning, C.; Gao, P.; Wang, B.; Lin, W.; Jiang, N.; Cai, K. Impacts of chemical fertilizer reduction and organic amendments supplementation on soil nutrient, enzyme activity and heavy metal content. J. Integr. Agric. 2017, 16, 1819–1831. [Google Scholar] [CrossRef]
  44. Si, L.; Peng, X.; Zhou, J. The suitability of growing mulberry (Morus alba L.) on soils consisting of urban sludge composted with garden waste: A new method for urban sludge disposal. Environ. Sci. Pollut. Res. 2019, 26, 1379–1393. [Google Scholar] [CrossRef]
  45. Zhang, Q.; Wang, C. Natural and Human Factors Affect the Distribution of Soil Heavy Metal Pollution: A Review. Water Air Soil Pollut. 2020, 231, 350. [Google Scholar] [CrossRef]
  46. Gong, L.; Ran, Q.; He, G.; Tiyip, T. A soil quality assessment under different land use types in Keriya river basin, Southern Xinjiang, China. Soil Tillage Res. 2015, 146, 223–229. [Google Scholar] [CrossRef]
  47. Tian, Y.; Xu, Z.; Wang, J.; Wang, Z. Evaluation of Soil Quality for Different Types of Land Use Based on Minimum Dataset in the Typical Desert Steppe in Ningxia, China. J. Adv. Transp. 2022, 2022, 1–14. [Google Scholar] [CrossRef]
  48. Jiang, Y.; Yang, Y.; Li, R.; Xi, B.; Li, M.; Hao, Y.; Meng, F.; Gao, S.; Chen, L. Soil Pollution Prevention Strategies and Typical Engineering Cases of Agricultural Products Producing Areas in Beijing–Tianjin–Hebei Region. Strateg. Study Chin. Acad. Eng. 2018, 20, 142–147. [Google Scholar] [CrossRef]
  49. Joimel, S.; Cortet, J.; Jolivet, C.C.; Saby, N.P.A.; Chenot, E.D.; Branchu, P.; Schwartz, C. Physico-chemical characteristics of topsoil for contrasted forest, agricultural, urban and industrial land uses in France. Sci. Total Environ. 2016, 545, 40–47. [Google Scholar] [CrossRef] [PubMed]
  50. Luo, L.; Ma, Y.; Zhang, S.; Wei, D.; Zhu, Y. An inventory of trace element inputs to agricultural soils in China. J. Environ. Manag. 2009, 90, 2524–2530. [Google Scholar] [CrossRef] [PubMed]
  51. Meng, Q.; Yang, J.; Yao, R.; Liu, G. Soil quality in east coastal region of China as related to different land use types. J. Soils Sediments 2013, 13, 664–676. [Google Scholar] [CrossRef]
  52. Jackson, P.C.; Meinzer, F.C.; Bustamante, M.; Goldstein, G.; Franco, A.; Rundel, P.W.; Causin, F. Partitioning of soil water among tree species in a Brazilian Cerrado ecosystem. Tree Physiol. 1999, 19, 717–724. [Google Scholar] [CrossRef] [PubMed]
  53. Zhu, P.; Zhang, G.; Wang, C.; Chen, S.; Wan, Y. Variation in soil infiltration properties under different land use/cover in the black soil region of Northeast China. Int. Soil Water Conserv. Res. 2024, 12, 379–387. [Google Scholar] [CrossRef]
  54. Zheng, Y.; Chen, T.; Chen, H. Lead accumulation in soils under different land use types in Beijing City. Acta Geogr. Sin. 2005, 60, 791–797. [Google Scholar]
  55. Li, L.; Zeng, X.; Bai, L.; Li, S. Characteristics of Lead Accumulation in Soils Under Different Agricultural Utilization Pattern in Shouguang of Shandong Province, China. J. Agro-Environ. Sci. 2010, 29, 1960–1965. [Google Scholar]
  56. Zou, X.; Zhu, X.; Zhu, P.; Singh, A.K.; Zakari, S.; Yang, B.; Liu, W. Soil quality assessment of different Hevea brasiliensis plantations in tropical China. J. Environ. Manag. 2021, 285, 112147. [Google Scholar] [CrossRef] [PubMed]
  57. Cong, W.; Hoffland, E.; Li, L.; Six, J.; Sun, J.; Bao, X.; Van Der Werf, W. Intercropping enhances soil carbon and nitrogen. Glob. Chang. Biol. 2015, 21, 1715–1726. [Google Scholar] [CrossRef] [PubMed]
  58. Sun, R.; Lan, G.; Yang, C.; Wu, Z.; Chen, B.; Fraedrich, K. Soil quality variation and its driving factors within tropical forests on Hainan Island, China. Land Degrad. Dev. 2023, 34, 3418–3432. [Google Scholar] [CrossRef]
Figure 1. Location map of the study area.
Figure 1. Location map of the study area.
Forests 15 00993 g001
Figure 2. Linear fitting of soil quality.
Figure 2. Linear fitting of soil quality.
Forests 15 00993 g002
Figure 3. Path model (PM) analysis and validity of the effects of chemical indicators, physical indicators, Pb, and land use types on the soil quality index (SQI). Red arrows indicate negative effects and green arrows represent positive effects. Numbers adjacent to arrows are path coefficients (p-values) indicating the effect size of the relationship. **, and *** represent 0.01 < p < 0.05, and p < 0.01, respectively.
Figure 3. Path model (PM) analysis and validity of the effects of chemical indicators, physical indicators, Pb, and land use types on the soil quality index (SQI). Red arrows indicate negative effects and green arrows represent positive effects. Numbers adjacent to arrows are path coefficients (p-values) indicating the effect size of the relationship. **, and *** represent 0.01 < p < 0.05, and p < 0.01, respectively.
Forests 15 00993 g003
Table 1. Basic information for each sampling site.
Table 1. Basic information for each sampling site.
Land Use TypesSlope/(°)Elevation/(m)Main Plant SpeciesCoverage/(%)
Newly cultivated land
(NCL)
0–540Eleusine indica (L.) Gaertn., Setaria viridis (L.) Beauv., Artemisia caruifolia Buch.-Ham. ex Roxb.10
Cropland (CL)0–5100Helianthus annuus L., Ipomoea batatas (L.) Lam., Artemisia argyi H. Lév. & Vaniot, Setaria viridis (L.) Beauv., Xanthium strumarium L.50
Arbor forest
(AF)
0–550Populus tomentosa Carrière, Larix gmelinii (Rupr.) Kuzen.80
Shrubland
(SL)
0–550Prunus mume ‘Meiren’60
Arbor–shrub mixed forest
(ASF)
0–590Populus tomentosa Carrière, Larix gmelinii (Rupr.) Kuzen., Vitex negundo var. heterophylla (Franch.) Rehd.90
Table 2. Soil indicators measurement.
Table 2. Soil indicators measurement.
Soil IndicatorsMeasurement Methods
Soil water content (SWC)The drying method (oven-drying method (105 °C, 12 h))
Soil bulk density (BD)
Maximum water-holding capacity (MWHC), capillary water-holding capacity (CWHC), capillary porosity (CP), non-capillary porosity (NCP)
The ring knife method
Soil organic carbon content (SOC)The potassium dichromate oxidation method
Soil total nitrogen (TN)Kjeldahl method
Soil total phosphorus (TP)NaOH melting-molybdenum antimony colorimetric method
pHThe PHS-3E meter (INESA, Shanghai, China) (the water–soil ratio was 2.5:1)
The metal elements (Cd, Cu, Pb, Fe, Zn)Inductively coupled plasma–optical emission spectrometry (PRODIGY-XP, Leeman Labs, Hudson, USA) (RF power: 1150 W; cooling air flow: 1.0 L/min; injection cleaning time: 30 s; integration time: 30 s; flushing pump speed: 45 r/min; analyzing pump speed: 45 r/min)
Table 3. The statistical features of soil indicators under different restoration modes.
Table 3. The statistical features of soil indicators under different restoration modes.
Soil
Index
Land Use Types
NCLCLAFSLASF
MeanMinMaxMeanMinMaxMeanMinMaxMeanMinMaxMeanMinMax
SOC
(g∙kg−1)
0.45 ± 0.09b0.100.961.18 ± 0.32ab0.332.980.91 ± 0.23b0.322.180.62 ± 0.13b0.251.371.79 ± 0.42a0.333.38
pH8.04 ± 0.03a7.868.177.26 ± 0.09c6.907.747.90 ± 0.04a7.718.097.61 ± 0.09b7.087.907.34 ± 0.15c6.788.09
TN
(g∙kg−1)
0.05 ± 0.02b0.010.190.20 ± 0.04a0.080.420.15 ± 0.03a0.040.320.15 ± 0.02a0.090.250.23 ± 0.05a0.080.50
TP
(g∙kg−1)
0.18 ± 0.06b0.020.521.52 ± 0.52a0.144.620.37 ± 0.06b0.180.690.72 ± 0.03b0.560.830.79 ± 0.10b0.361.23
SWC
(%)
14.93 ± 0.42a12.8418.5911.67 ± 0.34b10.4213.968.46 ± 0.69c4.6212.288.16 ± 0.67c4.9212.6715.69 ± 0.37a12.8917.00
BD
(g∙cm−3)
1.05 ± 0.02d0.971.161.51 ± 0.05a1.271.781.36 ± 0.02b1.231.461.26 ± 0.02c1.181.381.24 ± 0.03c1.101.41
MWHC
(%)
38.98 ± 0.32a37.1440.9126.17 ± 1.65c18.7136.2531.82 ± 0.87b27.3737.2233.60 ± 1.03b28.7638.8633.19 ± 1.16b26.3940.62
CWHC
(%)
31.29 ± 0.35a29.1233.7221.60 ± 1.44c16.3432.6827.00 ± 0.74b24.0131.5028.63 ± 1.39ab22.0337.7027.37 ± 1.01b20.9934.55
CP
(%)
32.71 ± 0.64bc28.2936.2032.13 ± 1.46c24.7544.6936.55 ± 0.93a32.6442.7036.01 ± 1.49ab27.7744.5233.63 ± 0.92abc28.1738.18
NCP
(%)
8.05 ± 0.36a5.089.536.65 ± 1.12a1.4713.306.43 ± 1.04a1.3411.696.30 ± 0.93a1.3611.597.07 ± 1.04a0.8813.96
Cd
(mg∙kg−1)
0.01 ± 0.00b0.010.010.05 ± 0.02a0.010.150.01 ± 0.00b0.010.020.01 ± 0.00b0.010.020.01 ± 0.00b0.010.02
Cu
(mg∙kg−1)
10.90 ± 1.07b5.3314.8823.61 ± 1.68a15.6530.208.15 ± 0.23b6.888.888.20 ± 0.49b6.6310.6025.78 ± 0.55a22.9027.90
Fe
(mg∙kg−1)
13.32 ± 0.29c12.0414.7724.68 ± 2.00a19.6436.0516.90 ± 0.79b12.4819.2517.22 ± 0.19b16.5618.2826.70 ± 0.47a24.1028.25
Pb
(mg∙kg−1)
10.95 ± 0.12d10.3511.2819.53 ± 1.09a16.1825.6310.87 ± 0.08d10.5811.3012.60 ± 0.56c11.2015.7516.80 ± 0.38b14.6817.85
Zn
(mg∙kg−1)
94.32 ± 3.11a83.10110.58106.28 ± 9.49a66.00152.3359.46 ± 2.32b47.7069.0063.02 ± 3.06b55.1580.38104.73 ± 3.82a89.70134.68
Notes: Different lowercase letters indicate significant differences among different land use types (one-way ANOVA, p < 0.05).
Table 4. Results of principal component analysis of the TDS.
Table 4. Results of principal component analysis of the TDS.
Soil IndicatorsPC1PC2PC3PC4Communalities
SOC0.470.75−0.06−0.290.92
pH−0.70−0.390.140.380.87
TN0.470.71−0.26−0.300.93
TP0.440.48−0.470.570.97
SWC0.270.420.780.110.89
BD0.54−0.60−0.51−0.140.95
MWHC−0.550.720.340.080.94
CP−0.270.50−0.46−0.310.98
CWHC−0.560.740.14−0.130.98
Cd0.350.53−0.470.590.97
Cu0.900.100.290.120.96
Fe0.83−0.270.13−0.080.88
Pb0.93−0.110.09−0.110.92
Zn0.630.040.590.130.77
Eigenvalue5.093.662.191.19
Variance (%)36.3926.1515.628.50
Cumulative variance (%)36.3962.5378.1586.65
Notes: Red highlighting factor loading values are considered highly weighted. Red highlighting and underlined loading values correspond to the soil indicators included in the MDS.
Table 5. Correlation matrix for the TDS.
Table 5. Correlation matrix for the TDS.
SOCpHTNTPSWCBDMWHCCWHCCPCdCuFePbZn
SOC1
pH−0.66 **1
TN0.84 **−0.74 **1
TP0.15−0.31 *0.38**1
SWC0.25−0.070.08−0.121
BD−0.07−0.30 *0.050.26 *−0.52 **1
MWHC0.130.27 *0.00−0.26 *0.46 **−0.91 **1
CWHC0.160.200.08−0.26 *0.33 *−0.75 **0.84 **1
CP0.19−0.050.29 *0.14−0.16−0.010.230.59 **1
Cd0.27 *−0.020.230.48 **−0.180.060.05−0.040.091
Cu0.48 **−0.50 **0.34 **0.130.56 **0.16−0.20−0.26 *−0.28 *0.111
Fe0.29 *−0.50 **0.30 *0.35 **0.200.50 **−0.52 **−0.52 **−0.180.090.78 **1
Pb0.40 **−0.77 **0.47 **0.28 *0.230.33 **−0.37 **−0.35 **−0.210.000.79 **0.81 **1
Zn0.22−0.210.040.020.69 **−0.09−0.04−0.12−0.36 **−0.230.72 **0.48 **0.51 **1
Notes: ** Correlation is significant at p < 0.01 level. * Correlation is significant at p < 0.05 level.
Table 6. The parameters and weights in the TDS and MDS.
Table 6. The parameters and weights in the TDS and MDS.
IndicatorsMeanSlopeWeight 1Weight 2
SOC0.99−2.500.0720.170
pH7.632.500.067
TN0.16−2.500.072
TP0.72−2.500.080
SWC10.89−2.500.0710.168
BD0.122.500.077
MWHC1.28−2.500.0770.159
CP0.33−2.500.052
CWHC0.34−2.500.0730.166
Cd0.072.500.0800.169
Cu0.41−2.500.076
Fe0.27−2.500.065
Pb0.022.500.0750.168
Zn15.32−2.500.063
Table 7. The soil quality evaluation results.
Table 7. The soil quality evaluation results.
Restoration PatternTDSMDS
NCL0.4660.485
CL0.4850.510
AF0.4930.519
SL0.4880.515
ASF0.5030.527
Table 8. The meanings of parameters in PLS-PM.
Table 8. The meanings of parameters in PLS-PM.
ParametersMeanings
Variance inflation factor (VIF)Measure the severity of multicollinearity
Reliability and validityValidate the accuracy of the PLS-PM results
Composite reliability (CR, not <0.6)Evaluate the reliability of internal coherence
Extracted average variance (AVE > 0.5)Verify the convergence validity
Table 9. Modeling (PLS-PM) of the variance inflation factor (VIF) of each indicator and the PLS-PM reliability and validity evaluations.
Table 9. Modeling (PLS-PM) of the variance inflation factor (VIF) of each indicator and the PLS-PM reliability and validity evaluations.
Driving FactorsSoil IndicatorsVariance Inflation Factor (VIF)Outer
Loadings
Cronbach’s AlphaComposite
Reliability (CR)
Average Variance Extracted (AVE)
Chemical indicatorsSOC1.190.930.570.730.69
TP1.190.71
Physical indicatorsCP1.460.810.720.850.77
MWHC1.460.94
Pb11///
Land use types11///
SQI11///
Notes: “/” is meaningless.
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

Qiang, F.; Sheng, C.; Zhang, J.; Jiang, L.; Zhou, J. Soil Quality Variation under Different Land Use Types and Its Driving Factors in Beijing. Forests 2024, 15, 993. https://doi.org/10.3390/f15060993

AMA Style

Qiang F, Sheng C, Zhang J, Jiang L, Zhou J. Soil Quality Variation under Different Land Use Types and Its Driving Factors in Beijing. Forests. 2024; 15(6):993. https://doi.org/10.3390/f15060993

Chicago/Turabian Style

Qiang, Fangfang, Changchang Sheng, Jiaqi Zhang, Liwei Jiang, and Jinxing Zhou. 2024. "Soil Quality Variation under Different Land Use Types and Its Driving Factors in Beijing" Forests 15, no. 6: 993. https://doi.org/10.3390/f15060993

APA Style

Qiang, F., Sheng, C., Zhang, J., Jiang, L., & Zhou, J. (2024). Soil Quality Variation under Different Land Use Types and Its Driving Factors in Beijing. Forests, 15(6), 993. https://doi.org/10.3390/f15060993

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