Next Article in Journal
Burned Area Detection Using Multi-Sensor SAR, Optical, and Thermal Data in Mediterranean Pine Forest
Next Article in Special Issue
Effects of Soil Amelioration and Vegetation Introduction on the Restoration of Abandoned Coal Mine Spoils in South Korea
Previous Article in Journal
Research on the Temporal and Spatial Distributions of Standing Wood Carbon Storage Based on Remote Sensing Images and Local Models
Previous Article in Special Issue
Soil Organic Carbon Pools and Associated Soil Chemical Properties under Two Pine Species (Pinus sylvestris L. and Pinus nigra Arn.) Introduced on Reclaimed Sandy Soils
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Distribution and Regulating Factors of Soil Nutrient Stocks in Afforested Dump of Pingshuo Opencast Coalmine, China

1
School of Land Science and Technology, China University of Geosciences (Beijing), Beijing 100083, China
2
Key Lab of Land Consolidation, Ministry of Natural Resources of the People’s Republic of China, Beijing 100035, China
*
Author to whom correspondence should be addressed.
Submission received: 17 January 2022 / Revised: 13 February 2022 / Accepted: 14 February 2022 / Published: 18 February 2022
(This article belongs to the Special Issue Forest Restoration and Secondary Succession)

Abstract

:
Determining the regulating factors of soil nutrient variations can guide the implementation of land reclamation measures in opencast coalmine regions. In this study, 132 soil samples were collected at 22 sample sites in the South Dump of Pingshuo opencast coalmine, and soil physicochemical properties were separately measured to obtain the related soil information. Geostatistical analyses were employed to analyze the spatial distribution patterns of soil organic carbon stocks (SOCD), total nitrogen stocks (TND), available phosphorus stocks (APD), and available potassium stocks (AKD) at 0–60 cm. The results showed that the spatial distributions of these soil nutrient stocks were characterized by moderate (TND) to strong (SOCD, APD, and AKD) spatial dependence. Meanwhile, the values of SOCD (16.4–60.1 Mg ha−1) and TND (1.9–15.5 Mg ha−1) were much higher than those of APD (0.022–0.095 Mg ha−1) and AKD (0.31–1.40 Mg ha−1). The statistical analyses indicated that the influence of afforestation on SOCD, TND, APD, and AKD was not significant, and the dynamic variations of soil nutrient contents were mainly regulated by soil pH in the South Dump. The findings of this study can provide some scientific guidance for soil nutrient management in the opencast coalmine regions of similar ecosystems.

1. Introduction

The cost-effective opencast mining is a widely employed exploitation method for mineral resources all over the world, the high-intensity and large-scale mining activities of which have caused severe damages to the ecosystems of original mining regions [1,2,3]. Moreover, most of China’s opencast coalmine is distributed in the ecologically fragile regions of Loess Plateau, where the local soils are mainly characterized by poor structure and low soil fertility [4,5,6]. Land reclamation has been proven to be the most effective method for the restoration of the ecological environment in opencast coalmine regions, which mainly comprises revegetation, soil reconstruction, and landform remodeling [7]. Soil fertility is an essential indicator for land reclamation quality, which can be directly reflected by soil nutrient contents [1,8]. Moreover, soil nutrient conditions also determine the spatial vegetation configurations in degraded opencast coalmine regions [1,9]. Therefore, it is of great significance to explore the spatial distribution patterns and regulating factors of soil nutrient stocks in the reclaimed land of opencast coalmine regions.
Soil organic carbon (SOC), total nitrogen (TN), available phosphorus (AP), and potassium (AK) are necessary macronutrients for the growth and productivity of vegetation, which can significantly regulate terrestrial functions by influencing soil physicochemical properties and the activities of soil microorganisms [10,11,12]. Generally, soil nutrient stocks are determined by the balance between organic matter input from regional biota and the losses from vegetation utilization and microbial respiration [13,14]. Previous studies have reported that the biogeochemical cycling of SOC, TN AP, and AK was greatly affected by natural processes (e.g., decomposition, leaching, and adsorption/desorption) and anthropogenic activities (e.g., fertilization, weeding, and afforestation) [15,16,17]. For instance, He et al. [18] explored the succession of soil C, N, P, and K stocks in Masson pine ecosystems and they found that the plantations of Masson pine significantly altered soil nutrient stocks and that soil physicochemical properties were the dominant influencing factors for these variations. Li et al. [19] found that the contents of SOC, TN, AN, and AK in woodland were much higher than those in grassland, shrubland, and cropland, but land-use types did not significantly affect other soil nutrient contents (e.g., TP and TK); they also drew the conclusion that climate factors, soil texture, and pH induced the variations of soil nutrient contents. The conclusions of regional-scale studies might contrast with each other due to the difference in climate, geological background, and soil conditions [20], which highlights the significance of regional studies into different ecosystems. Different from natural soils, reclaimed soils are much more heterogeneous and compact due to the irregular mixing of different materials (e.g., plant residues, fly ash, and soil particles) and the compaction of it by heavy machinery [21,22]. However, most studies were carried out in natural ecosystems, whose findings cannot directly guide the land reclamation practice in the ecosystems of opencast coalmine. Moreover, most studies only estimated the stocks of SOC (SOCD), TN (TND), AP (APD), and AK (AKD) at 0–20 cm, while recent studies revealed that soil nutrient stocks in deep soil horizon are also obviously influenced by anthropogenic and natural disturbances [23,24,25]. Therefore, more studies are required to explore the variations and influencing factors of the SOCD, TND APD, and AKD in deep soil horizons.
The soil nutrient management in reclaimed lands requires accurate spatial distribution patterns of soil nutrient stocks and deep understandings of the regulating factors for soil nutrient content variations [10]. Afforestation has been frequently employed for the reclamation of opencast coalmine regions because trees own a strong ability to increase soil organic matter (SOM) content and improve the activities of soil microorganisms [26]. The effects of afforestation modes and soil physicochemical properties on the recovery of reclaimed soil nutrients remain controversial and should be separately examined to provide some insight for the land reclamation in opencast coalmine regions. Therefore, the main objectives of this study were to: (1) estimate the nutrient stocks based on soil nutrient contents, soil bulk density, and soil layer thickness in different soil horizons; (2) present and analyze the spatial distribution characteristics of SOCD, TND, APD, and AKD at 0–60 cm based on the ordinary Kriging interpolation method and geostatistical analyses; (3) explore how afforestation and soil properties regulate the variations of soil nutrients in reclaimed soil profiles.

2. Materials and Methods

2.1. Study Area

This study was carried out in the South Dump of Pingshuo opencast coalmine that is located in the north of Shanxi Province and east of the Loess Plateau (39°23′–39°37′ N, 112°11′–113°30′ E) (Figure 1). The Loess Plateau has suffered from extensive mining, cultivation, and deforestation in recent years, where the ecological environment is increasingly threatened by desertification, declining soil fertility, land degradation, and high-intensity soil erosion [27]. Moreover, more than two-thirds of China’s coal resource is distributed across the Loess Plateau and opencast mining has inevitably caused damages to local ecosystems [28], for which the Chinese Government has initiated many ecological restoration or land reclamation programs to improve the natural environment in these regions. It has been more than 30 years since the South Dump was reclaimed, and the local ecosystem has been greatly improved. More studies should be carried out to extract the effective land reclamation measures used (e.g., vegetation configuration modes, soil reconstruction methods, and landform remodeling schemes) and fully understand the mechanisms by which these measures influence the evolution path of mining ecosystems.
The South Dump is dominated by a typical semi-arid and arid continental monsoon climate, where the average annual precipitation is about 450 mm (196–757 mm) and the average annual temperature is 4.8–7.8 °C, with the highest temperature of 34.5 °C and lowest temperature of −27.3 °C [1]. Moreover, most rainfall (>65%) is unevenly distributed from June to September, and the average annual evaporation (2160 mm) is much higher than rainfall [21,28]. The main reconstructed soil is classified as chestnut brown soil with low-content organic matter and poor soil structure, making it difficult for the land reclamation in the study area [29]. Afforestation has been initiated since the 1990s to increase the surface vegetation coverage rate in the South Dump and a multilevel ecological landscape has been constructed by planting various vegetation species, such as Ulmus pumila, Robinia pseudoacacia, Pinus tabuliformis, and Elymus Dahuricus Turcz [30].

2.2. Sampling and Analysis

In May and August 2018, 22 representative soil profiles were randomly dug across the South Dump (Figure 2) and soil samples were collected in different soil horizons (0–10, 10–20, 20–30, 30–40, 40–50, and 50–60 cm). Overall, these soil profiles were dug in the forest land (sites S1, S2, S4, S5, S6, S8, S14, S18, S20, S21, S22, S32, S33, and S34) and the non-forest land (sites S3, S13, S16, S17, S23, S28, S31, and S35), respectively. Before soil physicochemical properties were measured, all soil samples were sufficiently air-dried in an oven and additional materials (e.g., plant residues and unwanted gravels) were removed by a 2 mm sieve. Soil texture was quantitively determined using a laser particle size analyzer with the precision of 1% (i.e., Longbench Mastersizer 2000; Malvern Instrument, Malvern, England) while soil pH was measured by a pH meter with the precision of 0.05 [31,32,33]. The soil organic matter (SOM) and TN contents were determined based on the wet oxidation method (KMnO4) and the Kjeldahl method [28,34,35], respectively. Soil AP was extracted by the sodium bicarbonate solution (NaHCO3) and measured through Olsen’s bicarbonate method, while soil AK was extracted by the ammonium acetate solution (NH4OAc) and measured through the atomic absorption spectroscopic analysis [36]. Soil bulk density and herbaceous biomass have been reported in our previous studies [21,30]. Most experiments were performed in Beijing Academy of Agriculture and Forestry Sciences.

2.3. Calculation of Reclaimed Soil Nutrient Stock

The stocks of soil nutrients (i.e., SOC, TN, AP, and AK) were calculated based on the following equation [27,37,38]:
N T S D h = i = 1 n N T S i × β i × t i × ( 1 δ i × 10 2 ) 10
where NTSDh represents the total amount of soil nutrients at 0-h cm per unit area (Mg ha−1); n refers to the number of soil layers and i represents the ith soil layer; NTSi, βi, ti, and δi refer to soil nutrient contents (g kg−1), soil bulk density (g cm−3), thickness (cm) of the ith soil layer, and the percentage (%) of coarse particles (>2 mm), respectively. Since the soil particles of loess in China are always smaller than 2 mm, the value of δ is normally considered as 0 [27,39]. In this study, the reclaimed soil nutrient stocks at 0–60 cm were calculated based on these measured soil properties.

2.4. Geostatistical Analysis

Geostatistical methods can be employed to present the spatially continuous distribution of soil nutrients based on some discrete points through optimal interpolation, which was established on the spatial autocorrelation principle [27,40,41,42]. In this study, the preferentially considered ordinary Kriging interpolation (OK) was adopted to describe the spatial distribution characteristics of different soil nutrient stocks at 0–60 cm. The optimal input parameters of OK should be determined by the application of a semivariogram function, the formula of which is shown as follows [1,21]:
γ ( h ) = 1 2 N ( h ) i = 1 N ( h ) [ H ( x i )   H ( x i + 1 ) ] 2  
where γ(h) is the semivariogram function value, N(h) denotes the total number of data pairs at different sample sites, and H(xi) represents the value of variable H at the sample site xi. Moreover, the C0/(C + C0) ratio (i.e., Nugget/Sill) is a widely used indicator to quantify the spatial autocorrelation or dependency degree of studied variables [38,43]. According to the C0/(C + C0) ratio, the spatial autocorrelation degree could be classified into strong spatial autocorrelation (<25%), moderate spatial autocorrelation (25–75%), or weak spatial autocorrelation (>75%) [15,44].
The cross-validation method has been widely employed to evaluate the reliability of OK interpolation results [38,45,46]. One datum would be considered unknown to construct a variation model based on the data left, and the omitted value would be estimated based on the constructed model. The procedure would be repeated many times to obtain the estimated values at different sample sites and the errors between these actual and estimated values would be calculated to assess the accuracy of the selected interpolation models (e.g., spherical model, linear model, Gaussian model, and exponential model). Here, the absolute mean error (AME), mean error (ME), and root mean square error (RMSE) were selected as the evaluation indices, whose calculation formulae are shown as follows:
M E = 1 n i = 1 n ( P i A i )
A M E = 1 n i = 1 n | P i A i |
R M S E = 1 n i = 1 n ( P i A i ) 2
where n is the number of soil samples, Pi and Ai represent the predicted and actual values, respectively. The geostatistical analyses were achieved in GS+ 9.0 and the evaluation indices were further calculated in Microsoft Excel 2016.

3. Results and Discussion

3.1. Distribution Characteristics of Measured Soil Properties, and Nutrient Stocks

3.1.1. Soil Texture, pH, and Herbaceous Biomass at Different Sample Sites

The soil texture composition, soil pH and herbaceous biomass distributions at different sample sites are shown in Figure 3. The sand contents varied dramatically (12.11–81.42%) when compared with clay or silt contents and most points were distributed at the bottom of the USDA soil texture triangle (Figure 3). Based on the American soil classification scheme, most soil samples can be classified as sandy loam (24.24%), silt loam (26.52%), and loam (40.91%) in the study area. On the other hand, the herbaceous biomass was extremely unevenly distributed at different sample sites, whose coefficient of variation (CV) reached 75.78%. After the long-term afforestation, large amounts of suitable woody species were planted across the South Dump, which greatly improved local ecosystems. The dominant species in the forest land comprised Robinia pseudoacacia, Populus alba, Ulmus pumila, and Pinus tabulaeformis, while those in the non-forest land mainly included Caragana Korshinskii and Artemisia Capillaris Thunb. The mean herbaceous biomass in forest land was 150.56 g m−2, while it was 88.47 g m−2 in non-forest land. No significant difference was found between them based on the independent-sample T test, which might result from the high varibaility of herbaceous biomass at different sample sites. Overall, soil pH ranged from 5.7 to 8.4 and most soil samples (96.2%) were alkaline in the study area. Moreover, the vertical variations of soil pH were quite limited in most soil profiles except for sites S1, S2, and S22. Previous studies reported that mining and land reclamation operations would trigger the increase of reclaimed soil pH [47,48], which could account for the alkaline soil environment in this study.

3.1.2. Statistical Characteristics of SOCD, TND, APD, and AKD in Different Soil Horizons

The statistical results of the studied soil nutrient stocks in different soil horizons are shown in Table 1. The SOCD showed extremely strong variability (CV > 250%) at 20–60 cm and seemed to increase with soil depth, which should be mainly attributed to the existence of coal gangue in soil profiles at several sample sites. During the early stages of land reclamation, some coal gangue was inappropriately backfilled in reclaimed soils due to the limitation of theories and techniques [49]. Based on our field observations, the abnormally high SOCD values were always found in the reclaimed soils where coal gangue was distributed, which changed the original distribution characteristics of SOCD and cannot reflect the actual conditions of reclaimed soils [50]. Therefore, the abnormal SOCD data at sample sites S2, S22, and S28 were omitted in the following analyses to avoid the interference of coal gangue. Meanwhile, the average values of TND, APD, and AKD in different soil horizons did not present clear vertical distribution patterns. Moreover, TND was characterized by moderate (0–20 cm) to strong (20–60 cm) variability while APD and AKD were characterized by moderate variability in all soil horizons, indicating that available soil nutrients were less variable than SOC and TN in the reclaimed soils.

3.2. Geostatistical Analyses of the Calculated Nutrient Stocks

3.2.1. Semi-Variogram Analyses

In this study, the stocks of SOC, TN, AP, and AK at 0–60 cm were calculated to examine their spatial distribution patterns across the South Dump. Overall, SOCD and TND conformed to a lognormal distribution while APD and AKD were normally distributed, which met the requirement of semi-variogram analyses. The log-transformed SOCD and TND have been transformed back to their actual values in the following analyses. Semivariograms have always been considered as the optimal analytical method for the spatial variability of soil properties [44], whose parameters are listed in Table 2. The optimal model for SOCD, APD, and AKD is the Gaussian model while that for TND is the Spherical model. The results of R2 (≥0.58) and RSS (≤ 2.65 × 10−2) indicated that these models were reliable to analyze the spatial variability of the soil nutrient stocks in this study. Meanwhile, the C0/(C + C0) ratios of SOCD, APD, and AKD were lower than 18% while that of TND was higher than 25%, which demonstrated that these soil nutrient stocks were characterized by moderate (TND) to strong (SOCD, APD, and AKD) spatial dependence. The results indicated that the spatial distributions of these soil nutrient stocks were mainly influenced by intrinsic factors (e.g., heterogeneity of reclaimed soil, soil pH, surface vegetation configuration, and soil texture) [21,51,52], although extrinsic factors could slightly affect the spatial distribution of TND. The cross-validation was carried out to examine the accuracy of the optimal models, which demonstrated that the interpolated results were reliable and the errors were acceptable (Table 3). Moreover, it should be noted that all RMSE values were greater than AME values for these models, indicating that the predicted nutrient stocks were slightly underestimated [21].

3.2.2. Spatial Distribution Patterns of SOCD, TND, APD, and AKD across South Dump

The spatial distribution maps of SOCD, TND, APD, and AKD are displayed in Figure 4. The stocks of SOC at 0–60 cm ranged from 16.4 to 60.1 Mg ha−1 and the high-SOCD regions were mainly distributed in the northeastern areas of the South Dump, which was roughly characterized by a striped distribution pattern (Figure 4). The TND values were lower than 3.7 Mg ha−1 in most areas while the high-TND areas were mainly concentrated in the northwestern regions of the South Dump. Available soil nutrients can be directly utilized by vegetation because they are always dissolved in soil solutions [53], for which it was of great significance to explore their spatial distribution characteristics. In this study, the APD values ranged from 0.022 to 0.095 Mg ha−1, while AKD values ranged from 0.31 to 1.40 Mg ha−1, much lower than SOCD and TND. Moreover, most areas were characterized by low APD (<0.051 Mg ha−1) and AKD values (<0.82 Mg ha−1), while several high-APD or high-AKD “points” were sporadically scattered across the South Dump (Figure 4).

3.3. Influencing Factors of SOCD, TND, APD, and AKD in Reclaimed Land

3.3.1. Afforestation

Numerous studies have explored the effects of afforestation on soil nutrient stocks in natural ecosystems [14,18,20,54,55]. It has remained controversial about the specific influence of afforestation on the evolution of soil nutrient stocks and different studies might draw different conclusions due to the comprehensive influence of climate, topography, and geological background [20]. Zhao et al. [56] explored the spatial distribution characteristics of SOC across the Wangmaogou Watershed on the Loess Plateau and they found that the SOC content at 0–20 cm was significantly influenced by land-use types. However, few studies have been carried out to examine how afforestation regulates the sequestration of both soil nutrients and available nutrients in the reclaimed land of opencast coalmines. As mentioned above, the forest land was mainly dominated by Robinia pseudoacacia, Pinus tabuliformis, and Ulmus pumila, while non-forest land was mainly covered by Caragana korshins and Artemisia Capillaris Thunb. In this study, soil samples were categorized into two groups based on the land use types (forest land vs. non-forest land) and the difference of the nutrient stocks between the two groups was compared by the Mann-Whitney U test (Figure 5). The results demonstrated that there was no significant difference between these soil nutrient stocks in forest land and non-forest land, indicating that the influence of afforestation on soil nutrient stocks was positive but not significant. The phenomenon could be attributed to that soil nutrient stocks were affected by numerous factors and the influence of single factors was difficult to be directly examined. Moreover, the influence of afforestation was also determined by reclamation time because the improvement of soil physicochemical properties and the accumulation of SOM was gradually achieved.
Understanding the mechanisms by which afforestation regulated nutrient stocks in reclaimed soils was of great significance for the optimization of land reclamation measures in opencast coalmine. Previous studies reported that afforestation can increase aboveground biomass and provide abundant raw materials for microorganisms; on the other hand, afforestation can improve soil properties and benefit the accumulation of soil nutrients [57,58,59]. For instance, our previous study found that surface runoff or soil erosion intensity in forest land was much lower than that in non-forest land [21], indicating that the forest land owned a stronger ability to retain soil nutrients when compared with non-forest land during soil erosion process. Moreover, afforestation can reduce soil bulk density and improve the infiltration of dissolved soil nutrients (e.g., AP and AK) along soil profiles [21], which can also benefit the retention of soil nutrients in reclaimed soils. More studies are required to explore the detailed mechanisms by which afforestation affects the balance of soil nutrient stocks. Overall, afforestation can be identified as an effective pathway to increase the stocks of SOC, TN, AP, and AK in the South Dump.

3.3.2. Soil Texture and pH

Soil texture is an essential property determining the mechanical composition of soil particles, which can significantly affect the cycling and storage of soil nutrients [20,32,60]. Generally, finer soil particles owned a stronger ability to absorb soil nutrients and can prevent them from being leached along soil profiles [60,61,62]. However, the effects of soil texture on SOCD, TND, APD, and AKD were very limited in this study. As shown in Figure 3, the main soil types (91.67%) were silt loam, loam, and sandy loam, which were characterized by similar clay contents. The similar clay contents might make it difficult to distinguish the effects of different soil texture on soil nutrient stocks, which supported the idea that the influence of soil texture varied regionally based on local soil conditions [63].
On the other hand, soil pH has been considered as the dominant factor regulating the bioavailability and turnover of soil nutrients [32,64], which highlights the significance of examining the interactions between soil pH and soil nutrients in reclaimed soils. To minimize the potential interference of surface vegetation and backfilled coal gangue, soil samples at 40–60 cm were adopted for regression analyses to explore the pure influence of soil pH on soil nutrient accumulation in reclaimed soils. Since soil nutrient stocks were significantly affected by soil bulk density that was not correlated with soil pH, soil nutrient contents instead of stocks were employed as dependent variables in the following analyses. Overall, the results indicated that SOC and TN contents were negatively associated with soil pH, while AP and AK contents were positively associated with soil pH. Meanwhile, the relationships between SOC, TN, and soil pH can be appropriately described by unary quadratic equations, while the relationships between AP, AK, and soil pH can be roughly described by sigmoid fitting equations (Figure 6). The ratio of total carbon to total nitrogen (TC/TN) was an essential parameter for the evaluation of SOM quality [32,65,66], which controlled the mineralization/immobilization processes of SOM. Since the soil inorganic carbon contents were quite low and can be negligible, SOC/TN was adopted to denote the decomposition rates of SOM in this study. Generally, a lower C/N ratio can denote the quicker depletion of SOC and TN in the process of decomposition by microorganisms [67]. As shown in Figure 6, the C/N ratio can be estimated by soil pH based on a logarithmic equation (Figure 6).
The relationships between SOC, TN, and soil pH are very complex, and it was difficult to identify the dominant influencing factor among them. On the one hand, the breakdown of SOM would lower soil pH by excreting organic acids accompanied by the increase of SOC and TN contents, which can account for the negative relationships between SOC, TN, and soil pH. On the other hand, soil pH was an essential parameter regulating the activity and community structure of microorganisms that were reported to be the main reason for the decomposition of SOM in terrestrial ecosystems [32,67,68,69]. In this study, soil pH was negatively associated with C/N ratio, denoting that relatively higher soil pH (5.7–8.4) would accelerate the utilization of SOC and TN by microorganisms and reduce SOC and TN contents [32,67,70]. Therefore, the negative association between soil pH and C/N ratio can support the conclusion that soil pH influenced SOC and TN contents by regulating microorganisms’ decomposition process. By contrast, available soil nutrients (AP and AK) were positively correlated with soil pH, indicating that a relatively higher pH benefited the accumulation of AP and AK in reclaimed soils. Since these soil samples were collected at 40–60 cm, where few plant roots can reach and utilize soil nutrients, the influence of plant roots on AP and AK would not be considered as the dominant factor for their negative relationships with soil pH. Previous studies also found that the adsorption of AP and AK by clay is strongly regulated by soil pH, and that alkaline soils own a stronger ability to absorb and accumulate available soil nutrients [71,72,73]. The results of this study supported the assumption that soil pH can regulate the adsorption/desorption processes of AP and AK in soil environment.

3.4. Suggestion for Nutrient Management in Reclaimed Soils of Opencast Coalmine

Soil fertility maintenance was essential for the recovery of the natural resilience of local soil ecosystems in the process of land reclamation [1,74]. The findings of this study can provide some inspiration for the soil nutrient management in similar ecosystems. The spatial distribution patterns of SOCD were much more complex than those of TND, APD, and AKD (Figure 4), indicating that the SOC stocks in the reclaimed land were more sensitive to environmental factors (e.g., surface vegetation and soil properties). Regular monitoring of soil nutrient stocks was required, and land reclamation measures should differ regionally based on soil conditions. Afforestation or revegetation was an effective pathway to improve soil nutrient stocks and the optimal vegetation configuration should be determined based on local soil conditions. The Robinia pseudoacacia, Pinus tabuliformis, and Ulmus pumila species were the suitable vegetation species in the South Dump, which can provide some reference for the vegetation configuration in the opencast coalmine of similar ecosystems. Moreover, soil pH was a vital regulating factor for the accumulation of SOC, TN, AP, and AK. Numerous studies have reported that the availability of soil nutrients to plant roots would reach maximum when soil pH ranged from 5.5–7 [72,73]. Meanwhile, this study showed that alkaline soils owned a stronger ability to absorb or accumulate AP and AK in deep soil horizons, but the utilization of SOC and TN would be accelerated in an alkaline soil environment. Measures should be taken to adjust the pH values of reclaimed soils, but specific decisions should be made depending on the stages of land reclamation, and different measures should be taken in surface soil and deep soil.

4. Conclusions

This study investigated the spatial distributions of SOCD, TND, APD, and AKD at 0–60 cm across the South Dump and explored the influence of afforestation and soil physicochemical properties on the variations of these soil nutrients. The stocks of these soil nutrients did not show clear vertical distribution patterns, although the SOCD in deep soil was much higher than that in topsoil due to the existence of backfilled coal gangue. Horizontally, the results of C0/(C + C0) ratios demonstrated that soil nutrient stocks at 0–60 cm were characterized by moderate (TND) to strong (SOCD, APD, and AKD) spatial dependence, indicating that their spatial distributions were mainly regulated by intrinsic factors, although extrinsic factors also slightly influenced the spatial distribution of TND. Moreover, the spatial distribution of SOCD was characterized by a striped distribution pattern while those of TND, APD, and AKD were characterized by “point distribution patterns”. On the other hand, the effect of afforestation on soil nutrient stocks appeared to be positive but not significant in our study, most likely due to the growth of some species in the non-forest land. Moreover, soil pH was the dominant controlling factor for the dynamic balance of SOC, TN, AP, and AK in reclaimed soils and their relationships can be roughly described by regression equations, while the influence of soil texture on these soil nutrients was quite limited in the study area. These results indicated that more measures should be taken to appropriately adjust reclaimed soil pH according to land reclamation stages and soil nutrient conditions.

Author Contributions

Conceptualization, M.Q., Y.C. and W.Z.; methodology, M.Q., Y.C. and W.Z.; software, W.Z. and Y.C.; validation, W.Z., M.Q. and Y.C.; formal analysis, W.Z.; investigation, S.W., S.L., M.Q., Y.C. and W.Z.; resources, M.Q. and Y.C.; data curation, W.Z., S.W., S.L. and M.Q.; writing—original draft preparation, W.Z.; writing—review and editing, W.Z., M.Q. and Y.C.; visualization, W.Z., M.Q., S.W., S.L. and Y.C.; supervision, M.Q. and Y.C.; project administration, M.Q. and Y.C.; funding acquisition, Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

The research and APC were funded by the National Natural Science Foundation of China (U1810107 and 41701607), and the Fundamental Research Funds for the Central Universities of China (2-9-2018-025 and 2-9-2019-307).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge Gubai Luo, Xinyu Kuang, Yuting Li and Bingjin Geng for their assistance with sample collection and laboratory experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guan, Y.; Zhou, W.; Bai, Z.; Cao, Y.; Huang, Y.; Huang, H. Soil nutrient variations among different land use types after reclamation in the Pingshuo opencast coal mine on the Loess Plateau, China. CATENA 2020, 188, 104427. [Google Scholar] [CrossRef]
  2. Zhao, Z.Q.; Shahrour, I.; Bai, Z.K.; Fan, W.X.; Feng, L.R.; Li, H.F. Soils development in opencast coal mine spoils reclaimed for 1–13 years in the West-Northern Loess Plateau of China. Eur. J. Soil Biol. 2013, 55, 40–46. [Google Scholar] [CrossRef]
  3. Maiti, S.K. Ecorestoration of the Coalmine Degraded Lands; Springer: New Delhi, India, 2013. [Google Scholar]
  4. Cao, Y.G.; Bai, Z.K.; Sun, Q.; Zhou, W. Rural settlement changes in compound land use areas: Characteristics and reasons of changes in a mixed mining-rural-settlement area in Shanxi Province, China. Habitat Int. 2017, 61, 9–21. [Google Scholar] [CrossRef]
  5. Zhang, L.; Wang, J.; Bai, Z.; Lv, C. Effects of vegetation on runoff and soil erosion on reclaimed land in an opencast coal-mine dump in a loess area. CATENA 2015, 128, 44–53. [Google Scholar] [CrossRef]
  6. Guimarães, D.V.; Silva, M.L.N.; Beniaich, A.; Pio, R.; Gonzaga, M.I.S.; Avanzi, J.C.; Bispo, D.F.A.; Curi, N. Dynamics and losses of soil organic matter and nutrients by water erosion in cover crop management systems in olive groves, in tropical regions. Soil Tillage Res. 2021, 209, 104863. [Google Scholar] [CrossRef]
  7. Hu, Z.Q.; Atkinson, K. Principle and method of soil reconstruction for coal mine land reclamation. J. China Coal Soc. 1998, 6, 761–768. [Google Scholar]
  8. Simmons, J.A.; Currie, W.S.; Eshleman, K.N.; Kuers, K.; Monteleone, S.; Negley, T.L.; Pohlad, B.R.; Thomas, C.L. Forest to reclaimed mine land use change leads to altered ecosystem structure and function. Ecol. Appl. 2008, 18, 104–118. [Google Scholar] [CrossRef] [PubMed]
  9. Vasu, D.; Singh, S.K.; Sahu, N.; Tiwary, P.; Chandran, P.; Duraisami, V.P.; Ramamurthy, V.; Lalitha, M.; Kalaiselvi, B. Assessment of spatial variability of soil properties using geospatial techniques for farm level nutrient management. Soil Tillage Res. 2017, 169, 25–34. [Google Scholar] [CrossRef]
  10. Guan, F.; Xia, M.; Tang, X.; Fan, S. Spatial variability of soil nitrogen, phosphorus and potassium contents in Moso bamboo forests in Yong’an City, China. CATENA 2017, 150, 161–172. [Google Scholar] [CrossRef]
  11. Hati, K.M.; Swarup, A.; Mishra, B.; Manna, M.C.; Waniari, R.H.; Mandal, K.G.; Misra, A.K. Impact of long-term application of fertilizer, manure and lime under intensive cropping on physical properties and organic carbon content of an Alfisol. Geoderma 2008, 148, 173–179. [Google Scholar] [CrossRef]
  12. Quilchano, C.; Maranon, T.; Perez-Ramos, I.M.; Noejovich, L.; Valladares, F.; Zavala, M.A. Patterns and ecological consequences of abiotic heterogeneity in managed cork oak forests of Southern Spain. Ecol. Res. 2008, 23, 127–139. [Google Scholar] [CrossRef] [Green Version]
  13. Peng, Y.Y.; Thomas, S.C.; Tian, D.L. Forest management and soil respiration: Implications for carbon sequestration. Environ. Rev. 2008, 16, 93–111. [Google Scholar] [CrossRef]
  14. Wellock, M.L.; Rafique, R.; LaPerle, C.M.; Peichl, M.; Kiely, G. Changes in ecosystem carbon stocks in a grassland ash (Fraxinus excelsior) afforestation chronosequence in Ireland. J. Plant Ecol. 2014, 7, 429–438. [Google Scholar] [CrossRef] [Green Version]
  15. Tang, X.; Xia, M.; Guan, F.; Fan, S. Spatial Distribution of Soil Nitrogen, Phosphorus and Potassium Stocks in Moso Bamboo Forests in Subtropical China. Forests 2016, 7, 267. [Google Scholar] [CrossRef] [Green Version]
  16. Marklein, A.R.; Houlton, B.Z. Nitrogen inputs accelerate phosphorus cycling rates across a wide variety of terrestrial ecosystems. New Phytol. 2012, 193, 696–704. [Google Scholar] [CrossRef] [PubMed]
  17. Liu, M.; Zhang, W.; Wang, X.; Wang, F.; Dong, W.; Hu, C.; Liu, B.; Sun, R. Nitrogen leaching greatly impacts bacterial community and denitrifiers abundance in subsoil under long-term fertilization. Agric. Ecosyst. Environ. 2020, 294, 106885. [Google Scholar] [CrossRef]
  18. He, J.; Dai, Q.; Xu, F.; Yan, Y.; Peng, X. Variability in Soil Macronutrient Stocks across a Chronosequence of Masson Pine Plantations. Forests 2022, 13, 17. [Google Scholar] [CrossRef]
  19. Li, Y.Q.; Ma, J.W.; Xiao, C.; Li, Y.J. Effects of climate factors and soil properties on soil nutrients and elemental stoichiometry across the Huang-Huai-Hai River Basin, China. J. Soils Sediments 2020, 20, 1970–1982. [Google Scholar] [CrossRef]
  20. Li, X.; Li, Y.; Peng, S.; Chen, Y.; Cao, Y. Changes in soil phosphorus and its influencing factors following afforestation in Northern China. Land Degrad. Dev. 2019, 30, 1655–1666. [Google Scholar] [CrossRef]
  21. Huang, Y.; Cao, Y.; Pietrzykowski, M.; Zhou, W.; Bai, Z. Spatial distribution characteristics of reconstructed soil bulk density of opencast coal-mine in the loess area of China. CATENA 2021, 199, 105116. [Google Scholar] [CrossRef]
  22. Pietrzykowski, M. Soil quality index as a tool for Scots pine (Pinus sylvestris) monoculture conversion planning on afforested, reclaimed mine land. J. For. Res. 2014, 25, 63–74. [Google Scholar] [CrossRef]
  23. Xu, H.W.; Qu, Q.; Li, P.; Guo, Z.Q.; Wulan, E.; Xue, S. Stocks and Stoichiometry of Soil Organic Carbon, Total Nitrogen, and Total Phosphorus after Vegetation Restoration in the Loess Hilly Region, China. Forests 2019, 10, 27. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, X.; Ma, J.; Ma, Z.W.; Li, L.H. Soil nutrient contents and stoichiometry as affected by land-use in an agro-pastoral region of northwest China. CATENA 2017, 150, 146–153. [Google Scholar] [CrossRef]
  25. Li, C.Z.; Zhao, L.H.; Sun, P.S.; Zhao, F.Z.; Kang, D.; Yang, G.H.; Han, X.H.; Feng, Y.Z.; Ren, G.X. Deep Soil C, N, and P Stocks and Stoichiometry in Response to Land Use Patterns in the Loess Hilly Region of China. PLoS ONE 2016, 11, e0159075. [Google Scholar] [CrossRef]
  26. Korkanç, S.Y. Effects of afforestation on soil organic carbon and other soil properties. CATENA 2014, 123, 62–69. [Google Scholar] [CrossRef]
  27. Liu, Z.P.; Shao, M.A.; Wang, Y.Q. Effect of environmental factors on regional soil organic carbon stocks across the Loess Plateau region, China. Agric. Ecosyst. Environ. 2011, 142, 184–194. [Google Scholar] [CrossRef]
  28. Zhou, W.; Yang, K.; Bai, Z.K.; Cheng, H.X.; Liu, F. The development of topsoil properties under different reclaimed land uses in the Pingshuo opencast coalmine of Loess Plateau of China. Ecol. Eng. 2017, 100, 237–245. [Google Scholar] [CrossRef]
  29. Cao, Y.; Wang, J.; Bai, Z.; Zhou, W.; Zhao, Z.; Ding, X.; Li, Y. Differentiation and mechanisms on physical properties of reconstructed soils on open-cast mine dump of loess area. Environ. Earth Sci. 2015, 74, 6367–6380. [Google Scholar] [CrossRef]
  30. Wang, S.; Cao, Y.; Pietrzykowski, M.; Zhou, W.; Zhao, Z.; Bai, Z. Spatial distribution of soil bulk density and its relationship with slope and vegetation allocation model in rehabilitation of dumping site in loess open-pit mine area. Environ. Monit. Assess. 2020, 192, 740. [Google Scholar] [CrossRef]
  31. Liu, M.; Han, G.; Li, X. Comparative analysis of soil nutrients under different land-use types in the Mun River basin of Northeast Thailand. J. Soils Sediments 2021, 21, 1136–1150. [Google Scholar] [CrossRef]
  32. Zhou, W.; Han, G.; Liu, M.; Li, X. Effects of soil pH and texture on soil carbon and nitrogen in soil profiles under different land uses in Mun River Basin, Northeast Thailand. PeerJ 2019, 7, e7880. [Google Scholar] [CrossRef] [Green Version]
  33. Liu, M.; Han, G.L.; Li, X.Q. Contributions of soil erosion and decomposition to SOC loss during a short-term paddy land abandonment in Northeast Thailand. Agric. Ecosyst. Environ. 2021, 321, 107629. [Google Scholar] [CrossRef]
  34. Parkinson, J.A.; Allen, S.E. A wet oxidation procedure suitable for the determination of nitrogen and mineral nutrients in biological material. Commun. Soil Sci. Plant Anal. 1975, 6, 1–11. [Google Scholar] [CrossRef]
  35. Nelson, D.W.; Sommers, L.E. Total Carbon, Organic Carbon, and Organic Matter. In Methods of Soil Analysis; Soil Science Society of America: Madison, WI, USA, 1996; pp. 961–1010. [Google Scholar]
  36. Long, O.H.; Seatz, L.F. Correlation of Soil Tests for Available Phosphorus and Potassium with Crop Yield Responses to Fertilization. Soil Sci. Soc. Am. J. 1953, 17, 258–262. [Google Scholar] [CrossRef]
  37. Fernández-Romero, M.L.; Lozano-García, B.; Parras-Alcántara, L. Topography and land use change effects on the soil organic carbon stock of forest soils in Mediterranean natural areas. Agric. Ecosyst. Environ. 2014, 195, 1–9. [Google Scholar] [CrossRef]
  38. Yao, X.; Yu, K.; Deng, Y.; Zeng, Q.; Lai, Z.; Liu, J. Spatial distribution of soil organic carbon stocks in Masson pine (Pinus massoniana) forests in subtropical China. CATENA 2019, 178, 189–198. [Google Scholar] [CrossRef]
  39. Zhao, F.; Kang, D.; Han, X.; Yang, G.; Yang, G.; Feng, Y.; Ren, G. Soil stoichiometry and carbon storage in long-term afforestation soil affected by understory vegetation diversity. Ecol. Eng. 2015, 74, 415–422. [Google Scholar] [CrossRef]
  40. Webster, R. Quantitative spatial analysis of soil in the field. In Advances in Soil Science; Springer: New York, NY, USA, 1985; pp. 1–70. [Google Scholar]
  41. Marchant, B.P.; Lark, R.M. Robust estimation of the variogram by residual maximum likelihood. Geoderma 2007, 140, 62–72. [Google Scholar] [CrossRef]
  42. Carlon, C.; Critto, A.; Marcomini, A.; Nathanail, P. Risk based characterisation of contaminated industrial site using multivariate and geostatistical tools. Environ. Pollut. 2001, 111, 417–427. [Google Scholar] [CrossRef]
  43. Zheng, H.; Wu, J.; Zhang, S. Study on the Spatial Variability of Farmland Soil Nutrient Based on the Kriging Interpolation. In Proceedings of the 2009 International Conference on Artificial Intelligence and Computational Intelligence, Shanghai, China, 7–8 November 2009; pp. 550–555. [Google Scholar]
  44. Cambardella, C.A.; Moorman, T.B.; Novak, J.M.; Parkin, T.B.; Karlen, D.L.; Turco, R.F.; Konopka, A.E. Field-scale variability of soil properties in central Iowa soils. Soil Sci. Soc. Am. J. 1994, 58, 1501–1511. [Google Scholar] [CrossRef]
  45. Addis, H.K.; Klik, A. Predicting the spatial distribution of soil erodibility factor using USLE nomograph in an agricultural watershed, Ethiopia. Int. Soil Water Conserv. Res. 2015, 3, 282–290. [Google Scholar] [CrossRef] [Green Version]
  46. Bogunovic, I.; Mesic, M.; Zgorelec, Z.; Jurisic, A.; Bilandzija, D. Spatial variation of soil nutrients on sandy-loam soil. Soil Tillage Res. 2014, 144, 174–183. [Google Scholar] [CrossRef]
  47. Shrestha, R.K.; Lal, R. Changes in physical and chemical properties of soil after surface mining and reclamation. Geoderma 2011, 161, 168–176. [Google Scholar] [CrossRef]
  48. Juwarkar, A.A.; Mehrotraa, K.L.; Nair, R.; Wanjari, T.; Singh, S.K.; Chakrabarti, T. Carbon sequestration in reclaimed manganese mine land at Gumgaon, India. Environ. Monit. Assess. 2010, 160, 457–464. [Google Scholar] [CrossRef] [PubMed]
  49. Zhou, W.; Cao, Y.; Wang, S.; Huang, Y.; Zhou, W.; Bai, Z. Deciphering the origin and controlling factors of mercury in reclaimed soils: A case study in Pingshuo opencast coalmine of China. Environ. Sci. Pollut. Res. 2022, 1–13. [Google Scholar] [CrossRef] [PubMed]
  50. Das, R.; Maiti, S.K. Importance of carbon fractionation for the estimation of carbon sequestration in reclaimed coalmine soils—A case study from Jharia coalfields, Jharkhand, India. Ecol. Eng. 2016, 90, 135–140. [Google Scholar] [CrossRef]
  51. Fu, W.; Jiang, P.; Zhao, K.; Zhou, G.; Li, Y.; Wu, J.; Du, H. The carbon storage in moso bamboo plantation and its spatial variation in Anji County of southeastern China. J. Soils Sediments 2013, 14, 320–329. [Google Scholar] [CrossRef]
  52. Bandyopadhyay, S.; Novo, L.A.B.; Pietrzykowski, M.; Maiti, S.K. Assessment of Forest Ecosystem Development in Coal Mine Degraded Land by Using Integrated Mine Soil Quality Index (IMSQI): The Evidence from India. Forests 2020, 11, 1310. [Google Scholar] [CrossRef]
  53. Wang, H.; Zhang, G.-h.; Li, N.-n.; Zhang, B.-j.; Yang, H.-y. Soil erodibility as impacted by vegetation restoration strategies on the Loess Plateau of China. Earth Surf. Process. Landf. 2019, 44, 796–807. [Google Scholar] [CrossRef]
  54. Varnagirytė-Kabašinskienė, I.; Žemaitis, P.; Armolaitis, K.; Stakėnas, V.; Urbaitis, G. Soil Organic Carbon Stocks in Afforested Agricultural Land in Lithuanian Hemiboreal Forest Zone. Forests 2021, 12, 1562. [Google Scholar] [CrossRef]
  55. Ahirwal, J.; Maiti, S.K. Development of Technosol properties and recovery of carbon stock after 16 years of revegetation on coal mine degraded lands, India. CATENA 2018, 166, 114–123. [Google Scholar] [CrossRef]
  56. Zhao, B.; Li, Z.; Li, P.; Xu, G.; Gao, H.; Cheng, Y.; Chang, E.; Yuan, S.; Zhang, Y.; Feng, Z. Spatial distribution of soil organic carbon and its influencing factors under the condition of ecological construction in a hilly-gully watershed of the Loess Plateau, China. Geoderma 2017, 296, 10–17. [Google Scholar] [CrossRef]
  57. Graham, J.; Voroney, P.; Coleman, B.; Deen, B.; Gordon, A.; Thimmanagari, M.; Thevathasan, N. Quantifying soil organic carbon stocks in herbaceous biomass crops grown in Ontario, Canada. Agrofor. Syst. 2019, 93, 1627–1635. [Google Scholar] [CrossRef]
  58. Agostini, F.; Gregory, A.S.; Richter, G.M. Carbon Sequestration by Perennial Energy Crops: Is the Jury Still Out? BioEnergy Res. 2015, 8, 1057–1080. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Felten, D.; Emmerling, C. Accumulation of Miscanthus-derived carbon in soils in relation to soil depth and duration of land use under commercial farming conditions. J. Plant Nutr. Soil Sci. 2012, 175, 661–670. [Google Scholar] [CrossRef]
  60. Gonçalves, D.R.P.; Sá, J.C.d.M.; Mishra, U.; Cerri, C.E.P.; Ferreira, L.A.; Furlan, F.J.F. Soil type and texture impacts on soil organic carbon storage in a sub-tropical agro-ecosystem. Geoderma 2017, 286, 88–97. [Google Scholar] [CrossRef]
  61. Gami, S.K.; Lauren, J.G.; Duxbury, J.M. Soil organic carbon and nitrogen stocks in Nepal long-term soil fertility experiments. Soil Tillage Res. 2009, 106, 95–103. [Google Scholar] [CrossRef]
  62. Tian, L.; Zhao, L.; Wu, X.; Fang, H.; Zhao, Y.; Yue, G.; Liu, G.; Chen, H. Vertical patterns and controls of soil nutrients in alpine grassland: Implications for nutrient uptake. Sci. Total Environ. 2017, 607–608, 855–864. [Google Scholar] [CrossRef]
  63. Wang, S.L.; Huang, M.; Shao, X.M.; Mickler, R.A.; Li, K.; Ji, J.J. Vertical distribution of soil organic carbon in China. Environ. Manag. 2004, 33, S200–S209. [Google Scholar] [CrossRef]
  64. Kemmitt, S.J.; Wright, D.; Goulding, K.W.T.; Jones, D.L. pH regulation of carbon and nitrogen dynamics in two agricultural soils. Soil Biol. Biochem. 2006, 38, 898–911. [Google Scholar] [CrossRef]
  65. Dignac, M.F.; Kogel-Knabner, I.; Michel, K.; Matzner, E.; Knicker, H. Chemistry of soil organic matter as related to C : N in Norway spruce forest (Picea abies(L.) Karst.) floors and mineral soils. J. Plant Nutr. Soil Sci. 2002, 165, 281–289. [Google Scholar] [CrossRef] [Green Version]
  66. Baldock, J.A.; Oades, J.M.; Nelson, P.N.; Skene, T.M.; Golchin, A.; Clarke, P. Assessing the extent of decomposition of natural organic materials using solid-state C-13 NMR spectroscopy. Aust. J. Soil Res. 1997, 35, 1061–1083. [Google Scholar] [CrossRef]
  67. Zhou, W.X.; Han, G.L.; Liu, M.; Zeng, J.; Liang, B.; Liu, J.K.; Qu, R. Determining the Distribution and Interaction of Soil Organic Carbon, Nitrogen, pH and Texture in Soil Profiles: A Case Study in the Lancangjiang River Basin, Southwest China. Forests 2020, 11, 532. [Google Scholar] [CrossRef]
  68. Dai, X.L.; Zhou, W.; Liu, G.R.; Liang, G.Q.; He, P.; Liu, Z.B. Soil C/N and pH together as a comprehensive indicator for evaluating the effects of organic substitution management in subtropical paddy fields after application of high-quality amendments. Geoderma 2019, 337, 1116–1125. [Google Scholar] [CrossRef]
  69. Gregorich, E.G.; Carter, M.R.; Angers, D.A.; Monreal, C.M.; Ellert, B.H. Towards a minimum data set to assess soil organic-matter quality in agricultural soils. Can. J. Soil Sci. 1994, 74, 367–385. [Google Scholar] [CrossRef] [Green Version]
  70. Wan, X.H.; Huang, Z.Q.; He, Z.M.; Yu, Z.P.; Wang, M.H.; Davis, M.R.; Yang, Y.S. Soil C:N ratio is the major determinant of soil microbial community structure in subtropical coniferous and broadleaf forest plantations. Plant Soil 2015, 387, 103–116. [Google Scholar] [CrossRef]
  71. Shaheen, S.; Tsadilas, C. Phosphorus Sorption and Availability to Canola Grown in an Alfisol Amended with Various Soil Amendments. Commun. Soil Sci. Plant Anal. 2013, 44, 89–103. [Google Scholar] [CrossRef]
  72. Zhao, J.; Dong, Y.; Xie, X.; Li, X.; Zhang, X.; Shen, X. Effect of annual variation in soil pH on available soil nutrients in pear orchards. Acta Ecol. Sin. 2011, 31, 212–216. [Google Scholar] [CrossRef]
  73. Barrow, N.J. The effects of pH on phosphate uptake from the soil. Plant Soil 2016, 410, 401–410. [Google Scholar] [CrossRef]
  74. Kumar, S.; Maiti, S.K.; Chaudhuri, S. Soil development in 2–21 years old coalmine reclaimed spoil with trees: A case study from Sonepur-Bazari opencast project, Raniganj Coalfield, India. Ecol. Eng. 2015, 84, 311–324. [Google Scholar] [CrossRef]
Figure 1. Location of Pingshuo opencast coalmine (a), current conditions of opencast mining (b), and slope characteristics across the South Dump (c).
Figure 1. Location of Pingshuo opencast coalmine (a), current conditions of opencast mining (b), and slope characteristics across the South Dump (c).
Forests 13 00345 g001
Figure 2. Some field photos of selected quadrat (a), soil profile excavation (b), and sampled soil profiles (c).
Figure 2. Some field photos of selected quadrat (a), soil profile excavation (b), and sampled soil profiles (c).
Forests 13 00345 g002
Figure 3. Characteristics of soil texture (a), pH and herbaceous biomass (b). The soil texture triangle diagram has been reported in our previous study [49]. Soil pH at different depths is displayed by bubble diagrams, where larger bubbles represent the soil samples in deeper horizons. The column in green represents forest land sites, while that in red represents non-forest land sites.
Figure 3. Characteristics of soil texture (a), pH and herbaceous biomass (b). The soil texture triangle diagram has been reported in our previous study [49]. Soil pH at different depths is displayed by bubble diagrams, where larger bubbles represent the soil samples in deeper horizons. The column in green represents forest land sites, while that in red represents non-forest land sites.
Forests 13 00345 g003
Figure 4. Interpolated maps of SOCD (a), TND (b), APD (c), and AKD (d) at 0–60 cm.
Figure 4. Interpolated maps of SOCD (a), TND (b), APD (c), and AKD (d) at 0–60 cm.
Forests 13 00345 g004
Figure 5. The average SOCD (a), TND (b), APD (c), and AKD (d) values in forest land and non-forest land. The same letters indicate that there is no significant difference between the two data groups. The error bars represent the standard errors of different groups.
Figure 5. The average SOCD (a), TND (b), APD (c), and AKD (d) values in forest land and non-forest land. The same letters indicate that there is no significant difference between the two data groups. The error bars represent the standard errors of different groups.
Forests 13 00345 g005
Figure 6. Results of regression analyses between soil pH and SOC (a), TN (b), AP (c), AK (d), and C/N (e). All coefficients and constants in regression equations are significant at 0.01 or 0.001 level.
Figure 6. Results of regression analyses between soil pH and SOC (a), TN (b), AP (c), AK (d), and C/N (e). All coefficients and constants in regression equations are significant at 0.01 or 0.001 level.
Forests 13 00345 g006
Table 1. The calculated SOCD, TND, APD, and AKD in different soil horizons.
Table 1. The calculated SOCD, TND, APD, and AKD in different soil horizons.
Depth (cm)MeanMedianMaximumMinimumCV (%)
SOCD Mg ha−1
0–1015.609.2458.193.4599.14
10–2010.134.2274.882.54169.89
20–3026.563.91304.631.76259.21
30–4031.203.21496.491.90339.75
40–5031.372.82559.121.46378.31
50–6026.723.42397.621.81323.00
TND Mg ha−1
0–100.790.581.970.3557.74
10–200.530.391.620.2568.78
20–300.730.395.110.19147.68
30–400.810.328.650.18219.75
40–500.830.349.920.17247.51
50–600.720.306.670.20192.11
APD Mg ha−1
0–10 0.00560.00500.01160.002639.83
10–200.00530.00420.01490.000463.30
20–300.00620.00470.01490.000463.83
30–400.00740.00610.01650.001863.69
40–500.00910.00720.02550.000560.17
50–600.00970.00860.03030.001062.42
AKD Mg ha−1
0–10 0.170.140.410.0652.12
10–200.100.090.180.0341.63
20–300.090.080.180.0441.98
30–400.100.090.220.0447.09
40–500.110.100.220.0344.37
50–600.110.100.250.0548.10
Table 2. Semi-variogram analyses of SOCD, TND, APD, and AKD at 0–60 cm.
Table 2. Semi-variogram analyses of SOCD, TND, APD, and AKD at 0–60 cm.
Soil nutrient StocksOptimal ModelC0/(C + C0) % aR2 bRSS cA0/m d
SOCD (Mg ha−1)Gaussian17.050.582.65 × 10−2971.68
TND (Mg ha−1)Spherical25.600.821.06 × 10−2242.00
APD (Mg ha−1)Gaussian0.200.668.28 × 10−8171.47
AKD (Mg ha−1)Gaussian0.140.749.33 × 10−4280.59
a Nugget/Sill, spatial heterogeneity ratio; b determination coefficient; c residual sum of squares; d Range, distance of spatial dependency.
Table 3. Cross-validation of OK interpolation for SOCD, TND, APD, and AKD at 0–60 cm.
Table 3. Cross-validation of OK interpolation for SOCD, TND, APD, and AKD at 0–60 cm.
Soil nutrient StocksOptimal ModelMEAMERMSE
SOCD (Mg ha−1)Gaussian−2.56 × 1001.05 × 1011.47 × 101
TND (Mg ha−1)Spherical−1.38 × 1002.69 × 1005.92 × 100
APD (Mg ha−1)Gaussian−1.07 × 10−31.83 × 10−22.25 × 10−2
AKD (Mg ha−1)Gaussian−1.21 × 10−21.81 × 10−12.69 × 10−1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, W.; Qian, M.; Wang, S.; Li, S.; Cao, Y. Spatial Distribution and Regulating Factors of Soil Nutrient Stocks in Afforested Dump of Pingshuo Opencast Coalmine, China. Forests 2022, 13, 345. https://0-doi-org.brum.beds.ac.uk/10.3390/f13020345

AMA Style

Zhou W, Qian M, Wang S, Li S, Cao Y. Spatial Distribution and Regulating Factors of Soil Nutrient Stocks in Afforested Dump of Pingshuo Opencast Coalmine, China. Forests. 2022; 13(2):345. https://0-doi-org.brum.beds.ac.uk/10.3390/f13020345

Chicago/Turabian Style

Zhou, Wenxiang, Mingjie Qian, Shufei Wang, Shengpeng Li, and Yingui Cao. 2022. "Spatial Distribution and Regulating Factors of Soil Nutrient Stocks in Afforested Dump of Pingshuo Opencast Coalmine, China" Forests 13, no. 2: 345. https://0-doi-org.brum.beds.ac.uk/10.3390/f13020345

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