Variations in Soil Physico-Chemical Properties along Slope Position Gradient in Secondary Vegetation of the Hilly Region, Guilin, Southwest China

Understanding the variation of soil physico-chemical properties along slope position gradients is essential for vegetation restoration and reconstruction, but how slope positions impact the soil physico-chemical properties in the secondary vegetation of hilly regions is poorly understood. To address these uncertainties, we examined the changes in soil physico-chemical properties and their relationships along the slope position gradient in secondary vegetation of the hilly region in Guilin, southwest China. The results showed that except for the soil water content, soil total phosphorus and soil total potassium which reached the highest value at the footslope, other soil physico-chemical properties reached the highest content in the middle slope, and most of the soil physico-chemical properties showed the lowest content in the upper slope or footslope. Moreover, Pearson’s correlation analysis revealed that there were no significant correlations between most of the soil physico-chemical properties and that the correlations between soil physico-chemical properties were not consistent across different slope positions. Additionally, the principal component analysis showed that the first 4 principal components together explained 84.32% of the total variation and might be interpreted as the change of soil total nitrogen, soil organic matter, soil available nitrogen, soil available potassium, soil water content and soil total potassium. Overall, our results strongly demonstrated that slope positions showed significant effects on most of the soil physico-chemical properties and would provide an important reference for the formulation of restoration strategies in different slope positions to facilitate vegetation restoration and reconstruction and the sustainable development of the ecological environment in the hilly region.


Introduction
Soil is an important natural resource, as a key component of the mountain ecosystem, and it underpins our very existence through providing a wide range of food, feed, fiber and timber production, as well as through earth system functions that support the delivery of other "ecosystem services" [1,2]. The suitability of the soil for those productions is based on the quality of the soil's physical, chemical and biological properties [3,4]. However, the spatial variability of soil properties in an ecosystem is a key step in the development of revegetation strategies for this hilly region. However, despite the international fame of this region, there are few studies about the soil physico-chemical properties of this region, and most of them are about forest vegetation. Accordingly, how slope positions affect the soil physico-chemical properties in secondary vegetation of the hilly region remains unknown, which greatly limits our understanding of the cooperative restoration process of vegetation and soil. Hence, it is especially important to investigate how slope positions affect the soil physico-chemical properties in the secondary vegetation of the hilly region.
Here, we explored the variation of soil physico-chemical properties and their relationships along the slope position gradient in the secondary vegetation of the hilly region, Guilin, southwest China. Specifically, we hypothesized that soil physico-chemical properties and their relationships will change with the slope position in this secondary vegetation ecosystem. To test our hypotheses, we investigated soil physico-chemical properties along the slope position gradient in secondary vegetation of the hilly region, Guilin, southwest China. Consequently, the main objectives of this study were to: (1) investigate how soil physico-chemical properties change in different slope positions and (2) explore the relationships between soil physico-chemical properties at different slope positions. Our results are expected to provide fundamental information about the soil physico-chemical properties of this area. Such information would be helpful to support the reasonable utilization and management of soil resources, and further to facilitate ecological restoration and reconstruction strategies and the sustainable development of the ecological environment in the hilly region.

Study Area
This study was conducted in a hilly region of Yanshan Town, Guilin City, Guangxi Zhuang Autonomous Region, southwest China (25 • Figure 1). The topography of this region was characterized as a hilly region, and the altitude ranges from 100 to 500 m above sea level (a.s.l.). This region has a typical mid-subtropical humid monsoon climate with a climate of hot summers and cold winters, and the mean annual temperature is 17. [8][9][10][11][12][13][14][15][16][17][18][19].1 • C with a frost-free period of 309 d and a mean annual sunshine duration of 1670 h (the China Meteorological Data Service Center; http://data.cma.cn). The mean annual precipitation varies from 1814 to 1941 mm, with a distinct wet season from May to October (more than 80% of precipitation occurs during the wet season) and a dry season from November to April of the following year; the annual average evaporation capacity change ranges from 1490 to 1905 mm. The region's soil is primarily dominated by zonal red soils derived from Quaternary red clay according to the Chinese soil classification system, which was classified as ultisols based on the United States Department of Agriculture (USDA) soil taxonomy [31]. The natural vegetation of this region is characterized by subtropical evergreen and deciduous broadleaved mixed forests. Historically, most of the forests have been severely disturbed with only small tracts of intact or semi-intact forest, and at present, the main vegetation type is secondary vegetation. The dominant woody species include Pinus massoniana var. Massoniana, Melastoma candidum D. Don, Schima superba Gardn. et Champ, Rhododendron molle G. Don, Ficus hirta Vahl.

Vegetation Investigation
The main topographical factor examined in this study was the hilly slope positions, therefore three different slope positions in an altitudinal transect with the same slope aspect were chosen as the sampling sites, namely: footslope, middle slope and upper slope. From September to October 2016, a total of 28 plots (20 m × 20 m) were arranged on random locations from the bottom to the top along the footslope position (10 plots), middle slope position (9 plots) and upper slope position (9 plots) of the Jiangjiaba heaped-up mountain in Yanshan Town, Guilin City, Guangxi Zhuang Autonomous Region, southwest China ( Figure 1). All the plots were on similar northwest aspects, with an altitude range from 195 to 372 m above the sea level (a.s.l.). The plots within our three different slope positions were considered to be true replicates as the distance between them exceeded the spatial dependence of most physico-chemical properties [32]. In each sampling plot, all the individual trees and shrubs were identified and measured for Latin name, diameter at breast height (DBH) or basal diameter, height and layer coverage. The nomenclature of the species followed Flora of China (English edition; http://foc.iplant.cn). Moreover, environmental variables such as geographical position, elevation, topography (such as slope aspect, slope degree and slope position), rock exposure rate and degree of human interference were surveyed and recorded for each sampling plot. The basic description information of the sample sites at each slope position is summarized in Table 1

Vegetation Investigation
The main topographical factor examined in this study was the hilly slope positions, therefore three different slope positions in an altitudinal transect with the same slope aspect were chosen as the sampling sites, namely: footslope, middle slope and upper slope. From September to October 2016, a total of 28 plots (20 m × 20 m) were arranged on random locations from the bottom to the top along the footslope position (10 plots), middle slope position (9 plots) and upper slope position (9 plots) of the Jiangjiaba heaped-up mountain in Yanshan Town, Guilin City, Guangxi Zhuang Autonomous Region, southwest China ( Figure 1). All the plots were on similar northwest aspects, with an altitude range from 195 to 372 m above the sea level (a.s.l.). The plots within our three different slope positions were considered to be true replicates as the distance between them exceeded the spatial dependence of most physico-chemical properties [32]. In each sampling plot, all the individual trees and shrubs were identified and measured for Latin name, diameter at breast height (DBH) or basal diameter, height and layer coverage. The nomenclature of the species followed Flora of China (English edition; http://foc.iplant.cn). Moreover, environmental variables such as geographical position, elevation, topography (such as slope aspect, slope degree and slope position), rock exposure rate and degree of human interference were surveyed and recorded for each sampling plot. The basic description information of the sample sites at each slope position is summarized in Table 1

Soil Sampling and Analysis
Soil samples were also collected from July to September 2016. According to the "S" type, after leaf litter was removed, five topsoils (0-15 cm depth) samples were collected from each 20 m × 20 m plot and then mixed to obtain one sample to represent the soil chemical properties of each plot.
Meanwhile, soil water content was measured by a portable Soil Moisture/Temperature/Salinity Meter (WET-2, Delta-T, England) at the same location where the soil samples were taken. The composite samples were placed in polyethylene bags and immediately transported to the laboratory. After manually removing the fine roots, stones and plant organic residues, each soil sample was air-dried at room temperature. Subsequently, samples were crushed, sieved through a 2 mm mesh sieve, and homogeneously mixed for further soil physico-chemical analysis.
The analyses of soil properties were determined using methods that have been described in soil agricultural and chemical analysis [33]. Briefly, soil pH (pH) was measured with a digital pH meter (FE20K, Mettler-Toledo, Zurich, Switzerland) using a soil: water ratio of 1: 2.5. Soil organic matter (SOM, g kg −1 ) was tested by using the dichromate oxidation method. Soil total nitrogen (TN, g kg −1 ) was determined by using automatic Kjeldahl analysis (KJELTEC TM 8400, FOSS Quality Assurance Co., Ltd., Hillerød, Denmark). Soil total phosphorus (TP, g kg −1 ) was analyzed by acid digestion with an H 2 SO 4 +HClO 4 solution. Soil total potassium (TK, g kg −1 ) was digested by the HF-HClO 4 -HNO 3 acid mixture and determined by the flame photometric method. Soil available nitrogen (AN, mg kg −1 ) was quantified using the alkaline hydrolysis diffusion method. Soil available phosphorus (AP, mg kg −1 ) was measured using the molybdenum blue colorimetric method after samples had been extracted with 0.5 M Na 2 CO 3 whereas the soil available potassium (AK, mg kg −1 ) was shaken for 30 min with 1 M ammonium acetate solution (1:10 w/v) and then analyzed by flame photometry. For each soil sample, measurements were repeated three times, and the average value was taken as the basic data of soil physico-chemical analysis.

Statistical Analysis
All data were tested for normality of distributions and homogeneity of variances prior to analysis. When data were not normally distributed, they were logarithmically transformed before analysis. The main statistical parameters (minimum, maximum, mean, median, variance, standard deviation, skewness coefficient, kurtosis coefficient and coefficient of variance (CV)) of soil physico-chemical properties at three different slope positions were calculated through the "describe" function in the Psych package of the R 3.4.2 software [34]. Differences in soil physico-chemical properties across the three different slope positions were compared by using the one-way analysis of variance (one-way ANOVA) method. Furthermore, the Tukey-Kramer honestly significant difference (HSD) test of the means was used when statistically significant differences (p < 0.05) were observed among slope positions. Subsequently, all relationships between soil variables were evaluated using Pearson's correlation tests, and the size of the correlation coefficient was detected by using the Student's t. In addition, Principal component analysis (PCA) was applied to tease apart differences and examine the contribution degree of the overall variation of each soil's physico-chemical properties. The site distribution map was created using the ArcGIS 10.2 software (ESRI, Inc., Redlands, CA, USA software version 10.2; http://www.esri.com/software/arcgis/arcgis-for-desktop). All the other statistical analyses and graphing were performed using R 3.4.2 software [34] with the additional packages "vegan" [35].

Descriptive Statistics of Soil Physico-Chemical Properties at Different Slope Positions
The descriptive statistics of all soil physico-chemical properties at different slope positions were summarized in Table 2. The coefficients of variation (CV) of soil physico-chemical properties at the footslope, middle slope and upper slope positions were 6.47%-89.51%, 2.15%-41.19%, and 4.71%-155.54%, respectively. The minimum and maximum variable coefficients of soil physico-chemical properties at the three different slope position gradients were soil pH (CV = 6.47%, 2.15% and 4.71%, respectively) and soil available phosphorus (CV = 89.51%, 41.19% and 122.54%, respectively), respectively.

Comparison of Soil Physico-Chemical Properties at Different Slope Positions
In accordance with our predictions, most of the soil physico-chemical properties showed significant differences between the three slope positions. Specifically, the soil water content (SWC) at three slope positions was shown as footslope > middle slope > upper slope, and the difference between each slope position was significant (Figure 2A). Soil pH (pH) was slightly acidic to strongly acidic, and there was no significant difference between the three slope positions ( Figure 2B). Soil organic matter (SOM) showed that the middle slope and upper slope positions were significantly higher than the footslope positions ( Figure 2C). Soil total nitrogen (TN) and available nutrients (AN, AP and AK) at the middle slope positions were significantly higher than at the footslope and upper slope positions ( Figure 2D, G, H and I, respectively). Soil total phosphorus content (TP) was significantly higher at the footslope positions than at the middle slope and upper slope positions ( Figure 2E). Soil total potassium (TK) at the footslope and middle slope positions was significantly higher than the upper slope positions ( Figure 2F). In total, except for the soil water content (SWC), soil total phosphorus (TP) and soil total potassium (TK), which reached the highest value at the footslope positions, all the remaining soil physico-chemical properties showed the highest content in the middle slope positions. In addition, most of soil physico-chemical properties were lowest in the upper slope or footslope positions.

Relationships between Soil Physico-Chemical Properties at Different Slope Positions
Pearson's correlation analysis revealed that there were no significant correlations between most of the soil physico-chemical properties, and that the correlation between soil physico-chemical properties was not consistent across the slope positions ( Figure 3). However, several soil physico-chemical properties in all three slope positions showed significant correlations. For example, the soil organic matter (SOM), soil total nitrogen (TN), soil available nitrogen (AN) and soil available potassium (AK) were significantly positively correlated with each other in all three slope positions ( Figure 3). Moreover, soil pH (pH) had a significantly negative and positive correlation

Relationships between Soil Physico-Chemical Properties at Different Slope Positions
Pearson's correlation analysis revealed that there were no significant correlations between most of the soil physico-chemical properties, and that the correlation between soil physico-chemical properties was not consistent across the slope positions ( Figure 3). However, several soil physico-chemical properties in all three slope positions showed significant correlations. For example, the soil organic matter (SOM), soil total nitrogen (TN), soil available nitrogen (AN) and soil available potassium (AK) were significantly positively correlated with each other in all three slope positions ( Figure 3). Moreover, soil pH (pH) had a significantly negative and positive correlation with soil organic matter (SOM) and soil total phosphorus (TP), respectively. In addition, at the footslope, the soil pH (pH) showed significant negative and positive correlations with soil organic matter (SOM) and soil total phosphorus (TP), respectively, and soil total phosphorus (TP) was significantly negatively correlated with soil organic matter (SOM). At the middle slope, the soil total phosphorus (TP) and soil available nitrogen (AN) were significantly positively correlated with soil total potassium (TK) and soil available potassium (AK), respectively. At the upper slope, soil water content (SWC) was significantly negatively correlated with soil total potassium (TK). Apart from that, no significant correlation was found for other soil physico-chemical properties at the three slope positions.
Sustainability 2020, 12, 1303 10 of 18 phosphorus (TP) and soil available nitrogen (AN) were significantly positively correlated with soil total potassium (TK) and soil available potassium (AK), respectively. At the upper slope, soil water content (SWC) was significantly negatively correlated with soil total potassium (TK). Apart from that, no significant correlation was found for other soil physico-chemical properties at the three slope positions.

Principal Component Analysis of Soil Physico-Chemical Properties
The PCA biplots ordination diagram is called the equilibrium contribution circle, and the radius of this circle represents a variable that would contribute equally to all dimensions of the PCA space. Therefore, for any given pair of axes, the variables that have vectors longer than this radius make a higher contribution than average and can be interpreted with confidence [36]. The weights for nine original variables along the first four principal components are shown in Table 3. According to the results of the PCA ordination, the eigenvalue of the first four ordination axes was greater than 1, which can together explain 84.32% of the total variance in those properties ( Figure 4, Table 3). Specifically, the first principal axis (PC1) accounted for 37.56% of the total variance and was strongly and positively related to soil total nitrogen (TN), soil organic matter (SOM), soil available nitrogen (AN) and soil available potassium (AK). The second principal component (PC2) explained 19.20% of the total variation and was strongly and positively correlated with soil water content (SWC). The third principal component (PC3) accounted for 14.89% of the total variation and showed a significant positive relationship with soil total potassium (TK). In addition, the fourth principal component (PC4) accounted for 12.67% of the total variation.

Principal Component Analysis of Soil Physico-Chemical Properties
The PCA biplots ordination diagram is called the equilibrium contribution circle, and the radius of this circle represents a variable that would contribute equally to all dimensions of the PCA space. Therefore, for any given pair of axes, the variables that have vectors longer than this radius make a higher contribution than average and can be interpreted with confidence [36]. The weights for nine original variables along the first four principal components are shown in Table 3. According to the results of the PCA ordination, the eigenvalue of the first four ordination axes was greater than 1, which can together explain 84.32% of the total variance in those properties ( Figure 4, Table 3). Specifically, the first principal axis (PC1) accounted for 37.56% of the total variance and was strongly and positively related to soil total nitrogen (TN), soil organic matter (SOM), soil available nitrogen (AN) and soil available potassium (AK). The second principal component (PC2) explained 19.20% of the total variation and was strongly and positively correlated with soil water content (SWC). The third principal component (PC3) accounted for 14.89% of the total variation and showed a significant positive relationship with soil total potassium (TK). In addition, the fourth principal component (PC4) accounted for 12.67% of the total variation.   Figure 2. The variables that have vectors longer than the radius of the equilibrium contribution circle are shown by pink line. The intersection angle between vector lines shows significant correlations between its corresponding variables, in which a sharp angle represents positive correlation, an obtuse angle represents negative correlation and a right angle represents no significant correlation relationship. Moreover, the right diagram is a scree plot and broken stick model to help compare the real eigenvalues with null model values generated by the broken stick model (green broken line) and mean eigenvalues (Keiser-Guttman criterion, purple horizontal line).   Figure 2. The variables that have vectors longer than the radius of the equilibrium contribution circle are shown by pink line. The intersection angle between vector lines shows significant correlations between its corresponding variables, in which a sharp angle represents positive correlation, an obtuse angle represents negative correlation and a right angle represents no significant correlation relationship. Moreover, the right diagram is a scree plot and broken stick model to help compare the real eigenvalues with null model values generated by the broken stick model (green broken line) and mean eigenvalues (Keiser-Guttman criterion, purple horizontal line).

Basic Characteristics of Soil Physico-Chemical Properties at Different Slope Positions
According to China's second national soil census nutrient classification standard [37], the average content of soil organic matter, total nitrogen and available nitrogen in this study reached the first level (extremely rich), and soil total potassium and available potassium reached the fifth level (deficit), while the average soil total phosphorus and soil available phosphorus content belong to the sixth level (extremely deficient). In general, the soil pH of our study area was slightly acidic to strongly acidic, and its soil organic matter, total nitrogen and available nitrogen were relatively rich, while soil total phosphorus, available phosphorus, total potassium and available potassium were relatively deficient. Compared with the soil of the subtropical evergreen and deciduous broadleaved mixed forest in the southwestern Hubei Province, China [38], the average soil pH, soil organic matter, soil available phosphorus and soil total potassium in this study area were relatively high, while the other soil physico-chemical properties were relatively low. Similarity, compared with the soil of the evergreen broadleaved forest in Zhejiang Province, China [39], the soil pH and soil organic matter in this region were relatively high, and the soil total nitrogen and total phosphorus were relatively low, while the other soil physico-chemical properties were similar with our study.
Furthermore, descriptive statistics showed that the variable range of physico-chemical properties in the middle slope position was relatively small, while the variation at the upper slope position was relatively large. This smaller variation found in the middle slope position may be attributed to the relatively suitable environmental conditions at this slope position so that the variation of soil physico-chemical properties was relatively low. On the contrary, the environmental conditions of the upper slope position were relatively poor, with a stronger wind, stronger surface evapotranspiration and shallow soil layer, and the soil physico-chemical properties vary considerably. Moreover, soil pH has the smallest variable coefficient among all three slope positions and belongs to weak variation. These results are consistent with some studies [4] that also found that soil pH had no significant difference at different slopes. The possible reason is that soil is a huge buffer body with a self-regulating ability, which can balance the changes of the external acid and alkali environment to some extent [40]. In addition, the variable coefficient of soil available phosphorus was highest among all three slope positions, which is possibly due to the difference in aboveground vegetation composition and related physiological and ecological processes [41].

Responses of Soil Physico-Chemical Properties to Different Slope Positions
Slope position has been regarded as one of the most important abiotic factors that determine the spatial heterogeneity of soil physico-chemical properties by controlling the pedogenic process on a local scale [8]. Different slope positions lead to different influencing factors such as light, heat, water and air, and soil properties may also show differences. Specifically, the soil from the upper slope position was shallow with little capacity for moisture retention and its groundwater level was deeper, and these surfaces are exposed to intense solar radiation and strong wind, which translates into higher evapotranspiration rates; this may result in a higher risk of drought and overheating [42][43][44][45]. By contrast, the footslope position tends to constitute cooler and more humid habitats with lower solar radiation and evapotranspiration rates, and deep soils capable of accumulating surface runoff water, variables translating into habitats with low risk of drought and thermal damage [46][47][48]. For instance, Zhu et al. [49] found that lower and toe slope positions were hot spots receiving lateral surface and subsurface recharge. Our results also found that soil water content at the upper slope position was significantly lower than that at the footslope and middle slope positions. These trends are likely related to the relatively lower altitude of the footslope position, which is related to weak solar radiation and easy-to-maintain relatively high soil water content. Conversely, the upper slope position leads to stronger surface evapotranspiration due to high solar radiation, and the steeper slope contributes to the surface runoff being transported down the slope, and thus the footslope position usually has higher soil water content due to the acceptance of water from the upper slope position. Additionally, our analysis demonstrated that most of the soil physico-chemical properties showed the highest value in the middle slope position, which is consistent with other studies [11,50]. This may be attributed to the fact that the middle slope position is relatively flat among the three slope positions, and the soil from the upper slope position is heavily deposited here, and thus results in the decrease of soil nutrients. Besides, the temperature and moisture of the middle slope position are better than the upper slope position, and the light conditions are better than the footslope position. In general, the habitat conditions of the middle slope position are relatively better. This was evident from the species richness and density of woody species in the middle slope position, which reached the highest value among all three slope positions. Notably, vegetation, which plays important roles not only in reducing soil erosion but also in affecting soil nutrients cycles, is a key medium of the interactions between slope position and soil nutrients [7]. Conversely, soil physico-chemical properties resulting from modifications of slope positions could also influence the redistribution of native vegetation. Therefore, we should pay attention to the interaction between vegetation and soil in future research.

Relationships between Soil Physico-Chemical Properties
Soil physico-chemical properties, as an important indicator of soil quality, possess not only obvious differences but also interactions and close associations with each other, which jointly determine the quality characteristics of the soil [51]. According to the results above, the soil organic matter, soil total nitrogen, soil available nitrogen and soil available potassium were significantly positively correlated with each other in all three slope positions. In agreement with our results, some studies have shown that soil organic matter is a complex mixture, which influences a series of soil properties and nutrient cycling, while soil enzymes are mainly adsorbed on the soil organic matter particles or combined with humus [52]. Moreover, the soil nitrogen and potassium are mainly derived from the decomposition and mineralization of soil organic matter, and the increase of soil organic matter is beneficial to the decomposition of animal and plant residues and humus by soil enzymes, thereby releasing nitrogen and potassium [53]. Besides, soil available nitrogen is the nitrogen element that can be directly absorbed and utilized by plants, which is mainly derived from the mineralization of soil total nitrogen [54]. Therefore, the increase of soil organic matter is conducive to the accumulation of soil total nitrogen, available nitrogen and available potassium.
Furthermore, as one of the most important physico-chemical properties of soil, soil pH directly affects plant growth and microbial activities, and most plants are suitable for growth in an environment where pH is neutral [55]. Conversely, excessive pH will inhibit plant root growth, thereby reducing the accumulation of soil organic matter. Hence, suitable pH is beneficial to the decomposition of soil organic matter by microorganisms and affects the release of phosphorus [56], which may be the reason that the soil pH at the footslope position showed a significant negative and positive association with soil organic matter and soil total phosphorus, respectively. However, soil total phosphorus has a significant negative correlation with soil organic matter at the footslope position, possibly because soil organic matter is mainly wrapped in the soil clay, while the clay has a protective effect on soil organic matter, making it difficult to mineralize and easy to accumulate. Moreover, as a critical element for soil, nitrogen is usually coupled with phosphorus and potassium, so this may be attributed to the positive correlation between nitrogen, phosphorus and potassium, resulting in the soil total phosphorus and soil available nitrogen at the middle slope position showing significant positive correlation with soil total potassium and soil available potassium, respectively. These results also agree with a growing body of research that describes positive correlations between nitrogen, phosphorus and potassium [57,58]. Additionally, the soil total potassium was highly restricted by soil parent material, and the potassium resulting from weathering of parent material was converged in the surface soil through the vegetation [59,60]. Accordingly, our study found that soil total potassium at the upper slope position was negatively correlated with soil water content, which may be due to the fact that higher soil water content might change soil structure by disrupting aggregates and surface crust, thereby leading to increased soil susceptibility to water and wind erosion and stimulating soil total potassium losses.

Relationships between Different Samples and Soil Physico-Chemical Properties
Principal component analysis (PCA), one of several multivariate methods, simplifies the complexity in high-dimensional variables by geometrically projecting them onto lower dimensions called principal components (PCs), in order to find the best summary of the data using a limited number of PCs [36,61]. Our study showed that the eigenvalue of the first four ordination axes is greater than 1, which can together explain 84.32% of the total variance, which was identified as soil total nitrogen, soil available nitrogen, soil organic matter, soil available potassium, soil water content and soil total potassium. Among these ordination axes, the first principal axis accounted for 37.56% of the total variance and represented a multivariate space for soil total nitrogen, soil organic matter, soil available nitrogen and soil available potassium, and the vector lines between them have a sharp angle. This further confirmed the results of the Pearson's correlation analysis, that is, the soil organic matter, soil total nitrogen, soil available nitrogen and soil available potassium were significantly positively correlated in all three slope positions. The second principal component explained 19.20% of the total variation and was strongly and positively correlated with soil water content. The third principal component accounted for 14.89% of the total variation and showed a significant positive relationship with soil total potassium. These findings seem to be consistent with our current research, which found that the Pearson's correlation analysis showed that soil water content was positively correlated with soil total potassium at the upper slope position.
In addition, Principal component analysis not only visually explains the relationship between different soil physico-chemical properties, but also shows the relationship between soil physico-chemical properties and sample plots. PCA ordination results showed that the soil total phosphorus and soil water content have significant effects on the soil physico-chemical properties at the footslope position. In accordance with this result, we also found that soil total phosphorus and soil water content showed the highest content at the footslope. Meanwhile, soil total nitrogen, soil available nitrogen and soil organic matter are important soil physico-chemical properties at the middle slope. This can be clarified by the fact that the highest contents of soil total nitrogen, soil available nitrogen and soil organic matter were found at the footslope. Overall, the results of the principal component analysis are basically similar to the results of the Pearson's correlation analysis and one-way analysis of variance, which further demonstrates the accuracy of the results of this study. As a consequence, the principal component analysis not only better reflects the relationship between the sample plots and the soil physico-chemical properties, but also the interrelationship between the sample plots and soil physico-chemical properties.

Conclusions
In summary, our study clearly demonstrates that the slope position showed significant effects on soil physico-chemical properties, and most of the soil physico-chemical properties showed the highest and lowest content in the middle slope position and upper slope or footslope positions, respectively. These results imply that soil physico-chemical properties should be explicitly considered to determine the impact of management practices on soil quality, and provide a reference for the species selection and vegetation layout planning in vegetation restoration and reconstruction in the hilly region. Specifically, when planting artificial forests in the hilly region, reasonable species collocation, population density and spatial structure allocation should be determined according to the differences of habitat characteristics and growth adaptability strategies of afforestation tree species in different slope positions, so that vegetation can be rationally configured to improve the utilization efficiency of plants for water, heat and nutrient resources. However, soil properties in the forest ecosystem are significantly influenced by multiple environmental factors such as land use. Hence, a further step for future studies would be required in order to clarify the complex interactive relationships among slope positions, soil properties (such as physico-chemical and biological), vegetation, land use and its history and management in this study area, as well as to develop soil physico-chemical and biological quality indicators as a tool to assist in the restoration of degraded ecosystems and sustainable land and vegetation management.