Optimization of Reclamation Measures in a Mining Area by Analysis of Variations in Soil Nutrient Grades under Different Types of Land Usage—A Case Study of Pingshuo Coal Mine, China

: The development of reclaimed mine soils is normally spatially heterogeneous, making the fine management and utilization of reclaimed mined lands difficult. Soil nutrient grading can provide a scientific basis for the precise regulation of soil nutrients, but few related studies are available in reclaimed mined areas. This study aimed to quantify the spatiotemporal variations in soil nutrient grades under different land-use types in a reclaimed mined area on the Loess Plateau, China. The study area was graded by four soil nutrients (soil available potassium (SAK), soil available phosphorus (SAP), soil total nitrogen (STN), and soil organic matter (SOM)), and the variation features of soil nutrient grades in the initial stage of reclamation under four land-use types (i.e., cultivated land, grassland, forestland, and barren land) were systematically characterized by geostatis-tical analysis, pedodiversity analysis, and correspondence analysis. The results show that during the initial five years after reclamation, the soil nutrient grades of most reclaimed areas increased from Grade V and VI to Grade I–IV, while the improvements were significantly heterogeneous. Notably, the four land-use types had distinct variation characteristics. The barren land had the lowest SAP level, whereas it had the highest proportion, and medium–high grades of SAK, STN, and SOM (88.3, 100.0, and 100.0%, respectively). In terms of quantitative structure, it had the lowest richness index ( � � , 2.5) and Shannon’s entropy index ( � � , 0.7) and the highest evenness index ( � � , 0.8). These results suggest that the barren land had relatively high and balanced nutrients, with the highest homogeneity among the four land-use types. The grassland had considerable improvement in all nutrients (especially SAP; 95.6% of the area had high SAP grades); however, its improvement was the most heterogeneous ( � � = 4.5, � � = 0.7). As the second-most heterogeneous land-use type ( � � = 4.0, � � = 0.8), the forestland had relatively low STN, SAP, and SAK levels due to high nutrient uptake and storage by tree species, but it had the highest proportion of area that reached high SOM grades (36.4%) and medium to high SOM grades (100.0%) due to its high community productivity. The cultivated land, which received fertilization for an additional three years, was the most imbalanced in terms of nutrients. It had the highest proportion of area that reached high SAP grades (98.0%); in contrast, its area proportions of low-grade SAK and SOM (69.0 and 32.9%, respectively) were the highest among the four land-use types. Based on the above comprehensive characterization of soil nutrient grade variation, guidance was given for fine management of reclaimed mined land and the optimization of reclamation measures.


Introduction
Coal resources are known as "black gold" and "industrial food" and are one of the main energy sources supporting China's social and economic development.In the past 30 years, China has been the country with the highest coal production globally.According to the National Bureau of Statistics, China's raw coal production reached 3.75 billion tons in 2019, an annual increase of 4.2% [1].While promoting the rapid development of China's economy, coal mining has inevitably destroyed the ecological environment of the mining area and damaged the mine's land resources.This is especially pronounced in opencast coal mining areas, where the original vegetation is devastated, the natural soil structure is disordered and the terrain is completely reformed [2][3][4].Land reclamation, as an effective measure for remediating the postmining environment, guiding ecosystem recovery into a desirable trajectory and recovering the utility value of postmining lands, has become an important means of coordinating coal resource extraction and land resource protection [5,6].
Soil provides nutrients, water, and physical support for vegetation growth, and soil reconstruction is a prerequisite for ecosystem restoration in mining areas.The quality of reconstructed soil largely determines the rate and degree of vegetation restoration [7,8].How to construct reclaimed soil with optimal soil environmental quality in the process of "geomorphologic remodeling, soil reconstruction and vegetation restoration" has always been a key scientific issue for scholars and managers in the field of land reclamation [9].
In the mineral extraction processes of opencast mining, the original topsoil is stripped and stored and then respread over mineral waste in land reclamation to support vegetation recovery.The soil properties undergo considerable changes during these intensive artificial disturbances [10][11][12].Current studies on reclaimed soil have mostly concentrated on the surface or subsurface layer of the soil [13], and their research fields mainly cover the physical and chemical properties [2,14], nutrient content [15], heavy metal pollution and evaluation [16], organic carbon storage rate [10,17], microorganisms [18], and so on.The research methods for evaluating reclaimed soil mainly include principal component analysis [19], correlation coefficients [20], fuzzy comprehensive evaluation [21], regressionkriging approaches [22] and the modified Nemerow index [23].This series of studies has evaluated reclaimed soil quality variation following mined land reclamation from different angles and analyzed the evolution characteristics of reclaimed soils under different conditions, providing guidance for constructing ecologically stable reclaimed mine soil.
The land-use type can adequately reflect the effects on land resources imposed by humans [24].In reclaimed mining areas, different land-use types lead to varying soil nutrient evolution characteristics due to differences in vegetation selection and artificial management [25,26].Land-use conversion can change the soil nutrient evolution track [27].Studying the spatial distribution of soil nutrients under different land-use types can provide a reference for the selection and adjustment of land types in mined land reclamation.Using descriptive statistics, two-way multivariate analysis of variance (two-way MANOVA), spatial statistics, and geostatistical analysis, existing studies have confirmed that soil nutrients such as total nitrogen (TN), total phosphorus (TP), available phosphorus (AP), and available potassium (AK) have significantly different characteristics under different land-use types [26,28].Soil nutrient grading quantifies the soil nutrient level according to measurement criteria, which can provide a scientific basis for the precise regulation of soil nutrients in a region [29].It is widely applied in agricultural land management and protection in China [30].In the Second National Soil Census, China developed a national standard of nutrient grading for soils of natural ecosystems or agricultural lands, i.e., the Second National Soil Census Nutrient Grading Standard.Currently, soil nutrient grading research either focuses on a specific soil nutrient element [29,31] or comprehensively evaluates various soil nutrients.Frequently used research methods of the latter usually include the comprehensive quality index [32], fuzzy evaluation [33], and the matter-element model [34].
However, little research has been conducted on nutrient grading of reclaimed mine soils.As a category of technosols [35,36], reclaimed mine soils have great spatial heterogeneity due to heterogeneous parent material, microtopography, and reclamation measures, and different land-use types could further increase this trend [26,28].The high spatial heterogeneity of reclaimed mine soils makes effective management and rational utilization of reclaimed lands difficult.Nutrient grading of reclaimed mined lands and its spatial-temporal variation can provide a scientific basis for fine management measures and the adjustment of reclamation orientation.When combined with geostatistical analysis, the variation of the characteristics of soil nutrients at different levels can be quantitatively characterized and visually displayed [37].Pedodiversity analysis has always been a hot topic in soil science research, and many soil scientists have widely discussed the concepts and methods of pedodiversity [38][39][40][41][42]. Characterization of soil nutrient grade diversities can intuitively display not only the richness of soil nutrient grades in each region but also the dispersion degree of soil nutrient distribution at each grade in a certain area.Correspondence analysis (CA) is a visualizable data dimensionality reduction method that can reveal the differences between different categories of the same variable and the corresponding relationships with different variables [43].Through application of the theory of correspondence analysis to the quantitative characterization of soil nutrient grades, two-or three-dimensional correlation diagrams can be used to clearly express the complex relationship between soil nutrient grades under different conditions.Overall, a comprehensive understanding of the characteristics and applicability of geostatistics, pedodiversity, and correspondence analysis can provide a new idea for the quantitative characterization of the nutrient grades in reclaimed coal mining areas.
This study graded reclaimed mined land in an opencast coal mine according to soil nutrient levels and further comprehensively characterized the variation in soil nutrient grades under different land-use types, aiming to provide a scientific basis for adjusting reclamation measures and optimizing land-use patterns in mining areas.The objectives of this study were: (1) to grade the reclaimed soil nutrients (soil organic matter, total nitrogen, available potassium, and available phosphorus) under different land-use types; (2) to quantify the spatiotemporal variation characteristics of reclaimed soil nutrient grades under different land-use types for enlightening the optimization of land reclamation measures.

Study Area
The Pingshuo opencast coal mine (112°10′-113°30′ E and 39°23′-39°37′ N) is located in the northern part of Shanxi Province and the eastern part of the Loess Plateau within the territory of Pinglu District and Shuocheng District.It is the largest opencast coal mine in China, consisting of three major mines: Antaibao, Anjialing, and Donglutian.The climate type is a typical temperate arid and semiarid continental monsoon climate.The annual rainfall in this region is between 428.2 mm and 449 mm and is mainly distributed in July, August, and September.The annual evaporation is between 1786.6 mm and 2598.0 mm, far exceeding the rainfall amount.The average temperature in the region is between 4.8 °C and 7.8 °C, and the largest temperature difference can be as high as 61.8 °C during the year.The regional soil type is Kastanozems according to World Reference Base for Soil Resources [35], with a 35 m thick Quarternary Loess beneath the soil horizon.The abundant Quarternary Loess is an ideal soil substitute material for land reclamation.According to previous tests by our research group, the loess has low content of total nitrogen (<0.5 g/kg), 0.62 g/kg total phosphorus, and about 20 g/kg total potassium [44].
In Antaibao opencast coal mine, the coal seam is burried over 300 m deep.On the strip profile of the mining pit, overburden materials (i.e., soil and rock stratum laying above the coal seam) were stripped to expose the coal seam, and then transported and dumped on the dump profile of the pit.In this way, the mining pit moves on, leaving vast areas of destructed land, i.e., the inner dump.In land reclamation, a loess layer of 1-1.5 m thick was covered on the overburden materials to support plant growth.As the loess was from the Quarternary Loess layer, buried beneath the original soil horizon, the reclaimed soils are typical anthropic technosols according to World Reference Base for Soil Resources [35].To avoid severe compaction, the rolling process of vehicles and machinery was avoided as much as possible during the restoration process.Our research area is a large platform on the inner dump, which is surrounded by slopes (Figure 1b).The formation of this area began in 2011 and was completed in 2012.The final elevation of the dump was 1474-1480 m.The total study area was 0.44 km 2 , divided into 36 plots each with an average size of 250 × 150 m.Revegetation was implemented in 2013, and the land-use types are cultivated land, forestland, grassland, and barren land (Figure 1).Before vegetation restoration, this area was fertilized in a unified way except for the barren land.A combination of organic and inorganic fertilizers was used, and a one-time base application was applied, which included 750 kg/hm 2 ammonium bicarbonate in autumn, 30.0 t/hm 2 organic fertilizer, 300.0 kg/hm 2 urea, and 300.0 kg/hm 2 diammonium phosphate before sowing.The cultivated vegetation is one crop per year, including buckwheat and maize.Continuous topdressing was carried out within three years of planting, with 10.0 t/hm 2 organic fertilizer, 50.0 kg/hm 2 urea, and 50.0 kg/hm 2 diammonium phosphate.Medicago sativa, a leguminous plant, was planted in the grassland.The grass seeds were sown at 50 kg/hm 2 , and the amount of replanting was 5-10 kg/hm 2 each year for three consecutive years.The mixed planting of trees and shrubs was selected, the area of Pinus sinensis, Populus tomentosa Carr, and Robinia pseudoacacia was 3 × 3 m, and that of Caragana microphylla and Hippophae rhamnoides was 1.5 × 1.5 m.The seedlings were thinned, rearranged, and fixed within three years.Vegetation reconstruction was not carried out on barren land, and the vegetation in this region mainly consists of wild grass seeds carried by the loess parent material covering the surface of the barren land, including Calamagrostis and wheatgrass; additionally, no other management and protection measures have been implemented in the barren land area.

Soil Sampling
In this experiment, soil samples were collected twice in 2013 and 2017.The first sampling was conducted in June 2013, and 78 soil samples were collected.On the basis of the initial sampling, 81 soil samples were collected in July 2017 according to the reclamation level of each land-use type.To ensure comparability between the two samplings, the layout of quadrats was consistent with the sampling process.According to the size of the plots and the principle of random sampling, two quadrats were established in each plot (three quadrats were established in large plots, and one quadrat was established in small plots).The distribution of quadrats is shown in Figure 1c.
The sampling process used a five-point sampling method; that is, five sampling points were randomly selected near the location of each quadrat, the soil collected from the five points was mixed, and approximately 1 kg of the sample was retained from the mixture.Since the distribution and change in soil nutrients in the topsoil (0-20 cm) were the most obvious and representative [45], the soil was sampled at 0-20 cm.Before the soil samples were tested and analyzed, the soil samples were tagged, air-dried, screened through a 2 mm sieve, and stored for later use.

Methods
To quantify the variation characteristics of soil nutrient grades under different landuse types in the opencast coal mine area on the Loess Plateau, the Pingshuo opencast coal mining area was selected as the study area.Spatial interpolation of soil nutrients was carried out by geostatistical methods and then graded into six grades for each soil nutrient.The soil nutrient grading referred to the Second National Soil Census Nutrient Grading Standard, and the specific grading methods are described in detail in Section 3.3.Next, pedodiversity theory was introduced to analyze the diversity characteristics of spatial distribution and the quantitative structure of the soil nutrient grades under different landuse types.Moreover, correspondence analysis was used to explore the relationship between the land-use type and the reclaimed soil nutrient grade.According to the quantitative analysis of the variation characteristics of the soil nutrient grades, the mechanisms of variation under different reclamation types were explored, which provided implications for improving the reclamation soil quality and optimizing the land-use types in mining areas.The research framework is shown in Figure 2.

Analytical Sources
Soil nutrients were determined after the soil samples were collected and treated, and soil organic matter (SOM) was determined by the potassium dichromate volumetric method [46].The semimicro Kelvin method was used to determine soil total nitrogen (STN) [46]; available soil potassium (SAK) was determined by the method of NH4OAC extraction and flame photometry at the concentration of 1 mol L −1 [47].The available soil phosphorus (SAP) was determined by NaHCO3 extraction with a concentration of 0.5 mol L −1 and molybdenum-antimony colorimetry [46].

Geostatistical Analysis
Geostatistical analysis was used to analyze the spatial variability and distribution characteristics of surface soil nutrients in the study area in 2013 and 2017.The optimal semivariance model and parameters of soil nutrients were fitted by a semivariogram, and interpolation analysis was used to evaluate the spatial value of soil nutrients in the unsampled part of the study area.The semivariogram model is as follows [48]: where γ(h) represents the semivariogram of h, h represents the spatial interval between two quadrats, (ℎ) represents the total logarithm of the quadrats separated by distance from ℎ, and ( ) and ( + ℎ) are the sample values for ( ) at locations and + ℎ [ = 1,2..., (ℎ)], respectively.
In this study, the semivariogram was calculated by GS + 7.0 (Gamma Design Software, Plainwell, MI, USA).The optimal semivariogram results of soil nutrients have been given in previous research by our group [28].According to the best fitting semivariance model and parameters, and based on the ArcGIS 10.2 interpolation analysis module (ESRI, Redlands, CA, USA), the ordinary kriging method was used to estimate the soil nutrient distribution of unsampled sites in the study area.

Soil Nutrient Grading
The Second National Soil Census Nutrient Grading Standard is the normative criterion for soil nutrient grading of natural ecosystems or agricultural lands in China, but it cannot adequately characterize the gradients of the infertile reclaimed mine soil in our study area.As the nutrient levels of the underdeveloped reclaimed mine soil were much lower than those of natural or agricultural soils, nearly all of the study areas corresponded to the lowest or lowest two grades according to the Second National Soil Census Nutrient Grading Standard, and therefore this standard cannot adequately distinguish the heterogeneity of soil nutrients in different parts of the study area (Table 1).Therefore, referring to this standard, we established soil nutrient grades for the study area by dividing the nutrient range of the soil samples into six equal intervals; Grade I represented the highest nutrient level, while Grade VI represented the lowest (Table 1).The study area was then graded according to these grading ranges.

Pedodiversity Analysis
Based on the soil nutrient grading of the study area, the richness index (S), Shannon's entropy index (H), and the evenness index (E) of soil nutrient grades in different land-use types were calculated.In this study, the modified Shannon entropy equation [49], which was improved from mathematical ecology [50], was used to calculate the diversity index: In Equations ( 2) and ( 3), H represents the diversity of the spatial distribution of soil nutrient grades, S is the number of spatial grids, pi is the ratio of the specific soil area to the total soil area in the spatial grids, and E is the evenness of each soil nutrient grade within the study area.In Equations ( 4) and ( 5), ′ represents the diversity of the quantitative structure of soil nutrient grades, ′ is the number of soil grades, ′ is the ratio of soil nutrient grades to the total soil area, and ′ is the evenness of various soil nutrient grades within the study area.

Correspondence Analysis
Correspondence analysis (CA) is a visualizable data dimensionality reduction method that can reveal the differences between different categories of the same variable and the corresponding relationships with different variables [43].The key step of CA is to sum the original data matrix X, which contains M samples and N variables, according to rows and columns, respectively, to obtain rows , columns , and sum T, and to calculate the probability matrix P and data transformation matrix Z of the original data.The calculations are as follows: Then, the covariance matrices A and B of the samples and variables were calculated: Factor analysis of A and B was performed, and the two most important common factors ( , , , ) were extracted.Thus, the data points ( , ) and ( , ) were plotted in the same Cartesian coordinate system, consisting of dimension 1 and dimension 2. In this way, the complex relationship between samples and variables could be analyzed intuitively.
In this study, the land-use types and soil nutrient grades were taken as samples, and the corresponding area of the samples was taken as a variable.SPSS 21.0 statistical software (SPSS, Chicago, IL, USA) was used to analyze the corresponding relationships between the four types of land use and the 24 soil nutrient grades (four types of soil nutrients, and each type of soil nutrient was divided into six grades).

Spatial and Temporal Patterns
By comparing the spatial distribution of soil nutrients at all grades in 2013 and 2017 (Figure 3), it can be seen that, in most parts, the grades of all four soil nutrient elements increased.However, the variation in the distribution characteristics was different between the nutrient elements.The relative difference in the spatial distribution of the SAP grade was almost unchanged from 2013 to 2017.The high-grade SAP was mainly distributed in cultivated land and grassland, while the low-grade SAP was mainly distributed in forestland and barren land.In 2013, there was little difference in the distribution of SAK among different land-use types, but it had a significant agglomeration effect, with higher SAK in the southeastern part of the study area than in the other areas.After five years of reclamation, the distribution of the SAK grade was relatively dispersed, and the areas with high-grade soil nutrient quality were mainly located in grassland and forestland.STN and SOM had similar spatial distribution characteristics.In 2013, low-grade STN and SOM were concentrated in the northwestern part of the study area, which was mainly composed of cultivated land.By 2017, the contents of STN and SOM showed a patchy distribution, and the high-grade STN and SOM were mainly distributed in the northern and southwestern parts of grassland, while the low-grade STN and SOM were still mainly distributed in cultivated land.

Quantitative Structure
The quantitative structure of soil nutrient classification under different land-use types is shown in Figure 4.In 2013, different land-use types had similar soil nutrient grades.The majority of areas of all land-use types were at the lowest two grades for all soil nutrients except for SAP.Due to nonuniform fertilization, cultivated land had relatively higher SAP grades than the other land-use types, with 53.7 and 12.7% of the areas at Grade III and Grade II, respectively.
In 2017, the soil nutrient grades of each land-use type increased; however, different land-use types showed varying quantitative structures (Figure 4).The disparity in the SAP grades between the four land-use types was obvious, following the order of cultivated land > grassland > forestland > barren land.For the other three soil nutrients, only a small proportion of the area reached high grades (Grade I and II) in each land-use type, while most areas were at medium to high grades (Grade III and IV).Hence, for SAK, STN, and SOM, the four different reclamation orientations were ranked by the proportion of area that reached Grade IV, Grade III, and Grade III, respectively.The ranking of SAK was barren land (88.3%) > grassland (66.0%) > forestland (38.5%) > cultivated land (31.0%); the ranking of STN was barren land (90.4%) > forestland (73.1%) > grassland (71.0%) > cultivated land (37.3%); the order of SOM was forestland (93.6%) > barren land (89.2%) > grassland (45.2%) > cultivated land (32.7%).In addition, it is worth mentioning that barren land had the lowest and the most homogeneous improvement of soil nutrient grades among the four land-use types, with its whole area being inferior to Grade II in terms of SAK, SAP, and SOM, and only 30.1% of the area reached STN II.

Diversity Characteristics of Soil Nutrient Grades
The study area was divided into 36 grids by land parcels, and the spatial distribution diversity index of each soil nutrient grade in the 36 grids was calculated (Figure 5).The S value represents the number of grids occupied by this type of soil nutrient grade in the study area.In 2013, the high values of S were mainly distributed in Grades V and VI, and a few of them were scattered in other soil nutrient grades.By 2017, high values of S were located in Grades II, III, and IV, showing that, at five years after reclamation, the soil nutrient grades had significantly increased in most areas.The E value represents the spatial dispersion degree of soil nutrients of various grades.The larger the E value is, the more discrete is the spatial distribution of such soil nutrients, and vice versa [42].In 2013, the soil nutrient grades with the highest E values were SAK-VI, SOM-V, STN-V, and SAP-IV.By 2017, the soil nutrient grades with the highest E values were SAK-V, SAK-IV, STN-III, SOM-III, and SAP-II.As the number of grids remained unchanged, the variation of the H value had a similar trend with that of the E value, i.e., the larger the H value is, the larger the E value is.The diversity of the quantitative structure of soil nutrient grades is shown in Table 2.
In the whole study area, the ′ values of SAK, STN, and SOM increased by three units from 2013 to 2017, while that of SAP decreased by one unit.The ′ values of the four kinds of soil nutrients all increased to a certain extent.The ′ value of all soil nutrient grades increased, except for that of STN, which remained the same.The diversity indices of each land-use type also changed.In 2013, the average ′ values of the soil nutrient grades of the four land-use types were in the order of barren land < cultivated land < forestland < grassland.The ranking of the average ′ value was barren land = forestland < cultivated land < grassland.The ranking of the average ′ value was forestland < barren land < cultivated land < grassland.In 2017, the average ′ value of the soil nutrient grades of the four land-use types was in the order of barren land < cultivated land < forestland < grassland.The ranking of the average ′ value was barren land < cultivated land = forestland = grassland.The average ′ value was in the order of grassland < forestland = barren land = cultivated land.Meanwhile, the average ′ of each land-use type increased, indicating that each of them became heterogeneous in terms of the soil nutrient grades.

Land-Use Type
The Whole Study Area Cultivated Land Forestland Grassland Barren Land

Correspondence Analysis between Soil Nutrient Grades and Land-Use Types
The correspondence analysis between land-use types and soil nutrient grades (Figure 6) showed that the cumulative variance contribution rates in 2013 and 2017 were 97.7 and 95.2%, respectively, indicating that most information could be explained through these two dimensions.In 2013, the angle between vectors of forestland and barren land (from the origin to FL and BL, respectively) was the smallest, while the angles between the vector of cultivated land and vectors of the other three land-use types were all obtuse angles, indicating that forestland and barren land were the most similar land-use types in terms of the soil nutrient grades, and cultivated land differed widely from the other three landuse types.In 2017, the four land-use types were distributed within four different quadrants, showing that after five years of reclamation, the characteristics of the four land-use types were all different.High grades of SOM-I, STN-I, STN-II, SAK-I, SAK-II, and SAP-II were all nearest to grassland, signifying that grassland had the highest proportion of highgrade soil nutrients, and it had a balanced improvement in different nutrients.This result may relate to the cultivation of alfalfa in grassland and indicate the potential lack of fertilization in alfalfa culture in post-mining areas [36].Cultivated land was surrounded by SAP-I and SOM-V, meaning that cultivated land had the highest proportion of high-grade SAP and low-grade SOM, which mainly affected by multiple effects of fertilization supplementation and crop uptake [51,52].

Development of Soil Nutrient Grades in the Loess Plateau Coal Mine Dump
According to the spatial distribution of soil nutrient grades in 2013 and 2017 (Figure 3), we found that all nutrient grades in the study area had significantly increased.The structural characteristics of the soil nutrient grades (Figures 4 and 5) showed that soil nutrients in most parts of the study area were at Grades V and VI, with only a small proportion at medium soil nutrient grades in 2013.By 2017, soil nutrients reached grades II, III, and IV in most areas and even reached Grade I in some areas, except the barren land, showing that the reconstructed soil nutrients of the mining dumps have been greatly improved through the active guidance of the soil formation process through four years of land reclamation.That is to say, the development of reconstituted soil nutrients in reclaimed dumps in coal mining areas have mainly come from the induced pedogenic processes.Moreover, the ′ and ′ values in the whole study area increased from 2013 to 2017 (Table 2), indicating that the soil nutrients in the study area were becoming heterogeneous.These results suggested that land reclamation promoted soil nutrient levels and increased the diversity of the soil nutrient grades.
However, reclaimed mine soil is still infertile when compared with natural and agricultural soils.According to the Second National Soil Census Nutrient Grading Standard (Table 1), SOM and STN were at Grade VI (the lowest grade), and SAP and SAK were mainly at Grade IV and Grade III.The low nutrient grades of reclaimed mine soil were mainly due to the short restoration period and the poor properties of the soil substitutes used for reclamation.In the area of the Loess Plateau, abundant loess is the soil parent material and is an ideal soil substitute for land reclamation.According to the Completion Standards on Land Reclamation Quality of China, in this area, the original soil can be abandoned in mining, and the reclaimed land should be covered with a loess layer of 1~1.5 m to reconstruct the soil [53].This standard also requires that the land productivity of the reclaimed mined land should reach the level of unmined land in five years.However, our results suggest that the soil nutrients remained at low levels and could not reach this standard in such a short period, even with intensive manual intervention such as fertilization and revegetation.For example, by 2017, the STN and SOM levels for reclaimed mined land were 0.18-0.40 and 2.3-5.2g kg −1 , respectively, much lower than those of surrounding unmined land (0.6 and 7.0 g kg −1 , respectively [26]).Obviously, a longer reclamation period is needed for reclaimed mine soils to reach the background nutrient levels in this region.

Spatial and Quantitative Variation
In the first year of reclamation, different land-use types had broadly similar soil nutrient grades (Figures 3 and 4).This was mainly due to the consistent reclamation technology (soil reconstruction and fertilization) in the whole study area.In the soil reconstruction process, the covering material and reclamation operation were homogenized [3,54].However, since differences in reclamation orientation, labor management, and vegetation types exist among different land-use types, by 2017, the soil nutrient grades significantly varied with the land-use types.
The cultivated land, on the one hand, had the highest area proportion of high-grade SAP (98.0%).On the other hand, it had the highest proportion of low-grade SAK (69.0%) and low-grade SOM (32.9%), and its STN level was obviously the lowest among all landuse types (Figure 4).Moreover, its SAP level was negatively correlated with the STN level (r = −0.715,Table 3).These results show that the cultivated land had imbalanced soil nutrients, i.e., high overall SAP level, but, in contrast, low overall SAK, STN, and SOM levels.As the whole above-ground part of crop (buckwheat and silage maize) was harvested, the main soil nutrient source and sink were fertilization and crop uptake, respectively.Consequently, the nutrient imbalance was probably due to the fact that sole plant species (monoculture of buckwheat or maize) uptake nutrient elements unevenly when compared to vegetation constituted by multiple species.
Grassland had the second highest proportion of high-grade SAP (95.6%) and highgrade SOM (18.8%).Its STN level was also relatively high and was distributed with the only patch of STN-I in the study area.Moreover, the STN level was positively correlated with the SOM level (r = 0.748, Table 3).These results show that revegetation with Medicago sativa can efficiently improve soil nutrients.Existing research has found that Medicago sativa has strong nitrogen-fixing capacity and decomposes easily, which can rapidly increase STN and SOM [55,56].Our results agree with this conclusion; furthermore, our results also suggest that this reclamation mode can effectively improve SAP.
Forestland had relatively low SAP, SAK, and STN levels, and a portion of the area remained unchanged at SAP-IV, SAK-VI, SAK-V, and STN-V from 2013 to 2017.A patch of forestland even decreased from SAK-V to SAK-VI (Figure 3).However, it had the highest proportion of high-grade SOM (36.4%), with another 57.2% of the area at SOM-III, which was much higher than that of the other land-use types (Figure 4).The slow improvement (even decrease) in the SAP, SAK, and STN levels and the high SOM level were mainly because the growth of tree species consumed a large amount of nutrients that were converted and stored in their branches and trunk; meanwhile, the forestland had the highest productivity among all land-use types, and the decay of litter increased the SOM [57].Note: * represents significant differences at the 0.05 confidence level, and ** represent significant differences at the 0.01 confidence level.
The barren land had the lowest improvement in the SAP level, with no area reaching high-grade SAP (Figure 4).However, it had the highest SAK and STN levels, with 88.3% of the area reaching SAK-IV and 90.4% of the area reaching STN-III, both of which ranked first among the four land-use types.Although it had no area that reached SOM-II, 89.2% of the area reached SOM-III, which was second only to that of the forestland (93.6%).In addition, it had only two or three grades for each soil nutrient, which were the lowest among all land-use types, indicating that its nutrient distribution was the most homogeneous.

Changes in Diversity Characteristics
Different land-use types also lead to varying structural characteristics of soil nutrient grades.According to the diversity indices of soil nutrient grades under different land-use types in 2017 (Table 2), the cultivated land, forestland, and grassland had the same average ′ of 1.0, which was higher than that of the barren land, suggesting that all three artificial restoration methods would lead to higher heterogeneity in soil nutrient grades than natural restoration.Moreover, the ′ value is not only affected by the richness of the soil nutrient grades but is also closely related to the area occupied by the soil nutrients of each grade [58,59].Although the three artificially restored land-use types had the same average ′, grassland had the highest average ′ and lowest ′ among them, indicating that grassland had the largest discrepancy in soil nutrient levels between high-grade and low-grade areas among all land-use types.
In a two-dimensional graph of CA, the further a point is from the center of the origin, the more obvious are the variable's features [60,61].In 2013, the grassland, forestland, and barren land were closer to the origin, while the distance from the cultivated land to the center was longer, indicating that the former three had no significant differences in the structure of the soil nutrient grades but were relatively different from the cultivated land.In 2017, the four land-use types were distributed separately in the four quadrants of the Cartesian coordinate system, and all grades of soil nutrients were scattered around the coordinate plane, indicating that the four land-use types had different variation trends and structure characteristics in terms of the soil nutrient grades.Moreover, high-quality soil nutrients, such as SAK-I, STNI, SAK-II, SOM-I, and SAP-II, were mainly distributed around the grassland, suggesting that the grassland had the highest proportion of highgrade soil nutrients, which demonstrates that Medicago sativa is an ideal species for reclamation in this area.This result is consistent with previous research in this area [62,63].As the land-use type subjected to the most intensive artificial management, cultivated land, it received extra fertilization in the first three years of reclamation but was mostly surrounded by low nutrient grades (such as SOM-VI, STN-V, and SOM-V), except for SAP-I, indicating that the overall nutrient level was lower than that in the grassland despite the high fertilization input.

Implications for Reclaimed Land Management and Reclamation Technology Optimization
As a typical arid and semiarid region, the Loess Plateau in China has eight national coal mine bases.Mining activities destroy the original ecosystem and utilizability of land resources, intensifying the contradiction between people and land.In this area, the mined lands must be reclaimed, and the land productivity must reach the medium level of the surrounding area in five years.Due to the abundance of soil parent material (i.e., loess), the original topsoil is commonly abandoned, and loess is used as a soil substitute in reclamation [64].Additionally, as the loess is abundant in total phosphorus (0.62 g/kg) and total potassium (about 20 g/kg) [44], even without fertiliztion, the soil nutrients increase with years of reclamation [63,65].However, our study has provided strong evidence that nutrient levels of the reconstructed soils cannot reach this standard in such a short period.Despite the limited improvement on the whole, different reclamation orientations, human interventions, and vegetation types have different influences on the accumulation or decomposition of soil nutrients [25,26,65,66].Analysis of the nutrient variation in reclaimed mine soils under different land-use types can provide guidance for the fine management of reclaimed mined land and the optimization of reclamation technology [26,64,67].
There has been a long-standing controversy over whether reclaimed mined land should be artificially revegetated or left for natural restoration, and the conclusion varies with bioclimate zones [68,69].In our study area, the grassland, forestland, and cultivated land were all artificially revegetated, while the barren land was spontaneously restored.From the perspective of soil nutrient development, the barren land, which received the lowest manual intervention, had much higher overall nutrient grades than the cultivated land five years after reclamation (Figure 4), although the cultivated land had been fertilized for an additional three years.Its STN level even exceeded that of the grassland, which was seeded with Medicago sativa, a leguminous species with a high nitrogen-fixing capacity.Moreover, unlike the artificially revegetated land-use types, whose soil nutrient levels showed high heterogeneity, the improvement of all nutrients was homogeneous on the barren land.However, the vegetation cover developed much slower than the artificially revegetated land-use types, leaving the reconstructed soil exposed to a high risk of erosion.Hence, spontaneous restoration could be a good option in reclamation on the Loess Plateau, provided that erosion prevention measures are implemented.
Among the artificially revegetated land-use types, grassland was the most economically important.With modest human intervention (reseeding but no extra fertilization), there was considerable improvement in all soil nutrients (especially in SAP), although this improvement was the most uneven (average ′ = 4.5, Table 2).As another human-induced revegetated land-use type, forestland had a lower area proportion, reaching medium-high SAK and SAP grades, than the grassland.The low nutrient level occurred because the development of a tree requires a large amount of nutrients, and this situation will be mitigated in the middle stage of reclamation [63,65].Moreover, owing to high community productivity, the forestland had the highest SOM level among all land-use types, and its area proportion of SOM-III (93.6%, Figure 4) was much higher than that of the grassland (45.2%).These two types of human-induced revegetation, namely, grassland and forestland reclamation, have their own advantages in promoting soil nutrient accumulation, and both of their vegetation develops quickly, providing effective protection for soils from erosion.
The cultivated land differed from other land-use types in terms of soil nutrient sources and cycling.It was continuously fertilized in the first three years after reclamation, providing a much more artificial input of nutrients.However, substantial nutrients were absorbed and assimilated by crops and then lost through harvesting, with little litter left to decay in the soil [51,52].Hence, the cultivated land had the lowest overall grades of SAK, STN, and SOM (Figure 4).This slow accumulation of STN and SOM on reclaimed cultivated land is consistent with the findings of existing research [5,70,71].In contrast, due to the imbalanced fertilization and nutrient uptake by monoculture crops, the cultivated land had the highest SAP grade.These results suggest that the improvement of soil nutrients on the cultivated land was limited and highly imbalanced, although the fertilization and manual inputs were highest among all land-use types.Consequently, for cultivated land reclamation, it is not cost-effective to plant crops directly and apply a high amount of fertilization.An economic strategy is to plant Medicago sativa with seeds first, in the first 3~5 years (the same management as that used for grassland reclamation), which can increase soil nutrient levels rapidly and with the lowest cost.

Limitations and Prospects
The initial stage of land reclamation is the fastest and most critical recovery stage for the reconstruction of soil properties.Studying the reclaimed soil nutrient variation in this stage is of great significance for the timely optimization of land reclamation measures.Therefore, this study focused only on the initial stage of land reclamation (the first and fifth years) to conduct reclaimed soil sampling design, testing, and statistical analysis.However, existing studies have shown that ecological restoration and nutrient accumulation on reclaimed mined land usually require a much longer time span [72][73][74].Therefore, long-term monitoring of soil nutrients is still needed in further studies.

Conclusions
Taking a representative opencast coal mine on the Loess Plateau as an example, this study applied the theories of geostatistics, pedodiversity, and correspondence analysis to grade the reclaimed mined land according to nutrient levels and to quantitatively characterize the spatiotemporal variation in nutrient grades under different land-use types.The results show that the overall nutrient level of reclaimed mine soil increased significantly five years after reclamation, but it remained inferior to that of the surrounding undisturbed soil and needed a longer period to reach the background level.The variation characteristics of soil nutrient grades under different land-use types were obviously different in terms of their spatiotemporal patterns, quantitative structures, and indices.The barren land, the only spontaneously revegetated land-use type, had the lowest area proportion of high nutrient grades, while the majority of this land-use type was at medium to high grades, and it had the most homogeneous and balanced nutrients.With proper erosion prevention measures (vegetation recovery by spontaneous restoration occurs slowly), this can be a good reclamation orientation on the Loess Plateau.As a human-induced revegetated land-use type, the grassland reached a relatively high level in all soil nutrients five years after reclamation, although it had the highest heterogeneity in the soil nutrient grades ( ′ = 4.5, ′ = 0.7).Considering the moderate manual input it received, this is an economical and practical reclamation orientation.The forestland had lower SAK, SPA, and STN levels than the grassland, but it had the highest SOM level, with as much as 36.4% of the area reaching high SOM grades.Although forestland requires large amounts of N, P, and K elements in the early stage of reclamation, it has the greatest community productivity, ecosystem service function, and ecosystem stability.The cultivated land had the highest area proportion of high-grade SAP (98.0%), low-grade SAK (69.0%), and lowgrade SOM (32.9%), and its STN level was obviously lower than that of the other land-use types.The low and imbalanced soil nutrients, regardless of the highest fertilization and manual input being implemented, suggest that it is not economical to plant crops with continuous fertilization in the early stage of reclamation.Instead, the more economical and effective strategy for improving land productivity is to plant Medicago sativa for 3-5 years.These results comprehensively reflected the variation characteristics of soil nutrient grades under different land-use types in the initial stage of reclamation.Based on these results, suggestions for optimizing reclamation technology and developing fine management of reclaimed mined land were proposed.Our study shows that soil nutrient grading and deep analysis of soil nutrient grades under different land-use types are good tools for adjusting artificial intervention and land-use management of reclaimed mined lands.

Figure 1 .
Figure 1.Geographical location and distribution of quadrats in the study region ((a) Pingshuo coal mine in Shouzhou city in Shanxi Province; (b) the land reclamation design diagram of the inner dump; (c) the distribution map of soil sampling sites; (d1) field pictures in the cultivated land; (d2) field pictures in the grassland; (d3) field pictures in the forestland; and (d4) field pictures in the barren land).

Figure 3 .
Figure 3. Maps of soil nutrient grades in the study area (CL.Cultivated land; GL.Grassland; FL.Forestland; and BL.Barren land).

Figure 5 .
Figure 5. Distribution diagram of the diversity index of different soil nutrient grades in the study area.

Table 1 .
The soil nutrient grades for the study area and comparison with the Second National Soil Census Nutrient Grading Standard.

Table 3 .
Pearson correlation coefficients between soil nutrients under different land-use types in 2017.