Abstract

The ecological environment is fragile in Xinjiang, so it is necessary to carry out land reclamation for mines to restore its ecology. The premise of mines land reclamation is to determine the direction of land reclamation, which requires the suitability evaluation for land reclamation. In this paper, the evaluation index system and suitability evaluation model for land reclamation of nonmetallic mines in Xinjiang Uygur Autonomous Region were established. This model was established by using factor analysis, cluster analysis, and discriminant analysis and tested by back-substitution. First, using 149 units of 21 nonmetallic mines as research samples, the samples were divided into 4 categories by a combination of factor and cluster analysis. Then, the samples were trained using a discriminant analysis method to establish the corresponding land reclamation suitability evaluation model. This model was verified by back-substitution with an accuracy of 98.7%, and only 2 of 149 samples were misclassified. Finally, the evaluation model was applied to the Dabancheng Toga Solo limestone mine in Urumqi. Evaluation analysis of 15 land reclamation units of this mine showed satisfactory results. The evaluation model developed here could serve as a powerful complement to the evaluation of land reclamation suitability in Xinjiang.

1. Introduction

Exploitation of mineral resources has always been the pillar industry of Xinjiang's economic development. Xinjiang is rich in nonmetallic mining resources, such as limestone, dolomite, marble, granite, jade, vermiculite ore, sand, and gravel [1, 2]. Up to now, just the nonmetallic mining permits issued by the Land and Resources Department of Xinjiang Uygur Autonomous Region have amounted to more than 450, with a combined area of over 650 km2. Nonmetallic mining methods are mainly open-pit mining (Figure 1), and large-scale mining is bound to severely damage the local ecological environment. Because the area of Xinjiang has an arid or semi-arid climate, the ecological environment is extremely fragile [3]. Land reclamation must be carried out after the pits are closed. An important prerequisite task and a specific challenge of land reclamation of nonmetallic mines in Xinjiang is to determine the direction of land types for land reclamation, that is, to evaluate the suitability results for land reclamation. Different land reclamation suitability evaluation results directly determine the future land reclamation measures [4]. It is necessary to establish an effective evaluation method for land reclamation evaluation of nonmetallic mines in Xinjiang.

In the past, many studies have produced fruitful results on the suitability evaluation of mining land reclamation using a variety of evaluation methods [5]. Previously established methods of evaluation mainly include the extremum value evaluation method [68], the GIS method [918], the discriminant analysis method used in the suitability analysis of coal mine reclamation [19], and the extension evaluation method [20]. Also the corresponding index system was established for evaluation [21, 22], fuzzy comprehensive evaluation method [2325], clustering analysis methods [26], analytic hierarchy analysis [27, 28], evaluation using neural network models [29], and evaluation using spatial decision support systems [30]. However, the above evaluation methods were mostly aimed at the suitability evaluation of coal mines, and for specific evaluation areas. Because the land reclamation suitability evaluation method of coal mines is mainly aimed at the land subsidence area formed by underground mining of coal mines, the way of land destruction is very different from that of nonmetallic mines, and the construction facilities of coal mines are also different from that of nonmetallic mines, so the above land reclamation suitability evaluation methods are obviously not suitable for the land reclamation work of nonmetallic mines. The evaluation indexes selected in the evaluation process are also different from place to place. To date, very few studies of land reclamation have taken into account the characteristics of damaged land in nonmetallic mines in Xinjiang and considered the special ecological environmental nature of the area.

At present, the suitability evaluation of mines land reclamation in Xinjiang is mainly carried out according to the relevant regulations and standards (Table 1) [31, 32]. Table 1 lists the results of the evaluation, which are categorized into four categories of very suitable (1), somewhat suitable (2), not very suitable (3), and unsuitable (4). The evaluation indexes include mainly four aspects: terrain slope, soil quality, irrigation and drainage conditions, and degree of pollution. Some indexes are difficult to quantify. The severity of land damage was inadequately considered, and the same evaluation indexes level may lead to two different results. For example, as shown in Table 1, for a reclamation unit with a terrain slope greater than 35°, the current land reclamation suitability evaluation grade would be “not very suitable” or “unsuitable” for reclamation as grassland. Such an evaluation result cannot determine in the first place whether the land reclamation unit is suitable or unsuitable for reclamation as grassland, and hence impossible to determine the land reclamation method and the reclamation workload. If it is unsuitable for reclamation as grassland, the land reclamation direction of this unit is bare land. A final land reclamation plan is only possible based on expert evaluation of the original land type, land use plan, and actual mine conditions.

In view of the above-mentioned research status and shortcomings of land reclamation suitability evaluation, we made a comprehensive study of the characteristics of damaged land and the ecological environmental conditions of nonmetallic mines in Xinjiang and determined the corresponding land reclamation suitability evaluation indexes system which is different from the previous research results according to relevant regulations and standards. The 21 nonmetallic mines are divided into 149 units according to the characteristics of land damage (e.g., the depth of pit excavation). The 149 land reclamation units were divided into four categories using combined cluster analysis method and classification factors without prior knowledge of specific classification results. Based on the accurate classification, a discriminant model was established for suitability evaluation of nonmetallic mines reclamation in Xinjiang. The land reclamation suitability evaluation model established in this paper could determine the land reclamation direction of land reclamation units for nonmetallic mines in Xinjiang at the first time, which is conducive to guiding land reclamation work. The evaluation model can be widely applied to the land reclamation of nonmetallic mines in Xinjiang.

2. Suitability Evaluation Indexes System of Mines Land Reclamation

The types of land damaged by nonmetallic mines in Xinjiang are mainly grassland and bare land. This land reclamation suitability evaluation is also aimed at grassland reclamation. A system of indexes in which the indexes are both mutually independent and complementary was constructed for evaluating the suitability of land reclamation. The indexes chosen in this study consider mainly the following factors: first, measurability; that is, the indexes contributing to the evaluation must be measurable and numerically quantifiable. Second, dominant; that is, the increase and decrease of the indexes affect the rise and fall of the quality of the land reclamation units. Third, nonoverlapping; that is, the boundaries between the evaluating indexes are clear and do not overlap. Given the above requirements, considerations were given to the unique geographical characteristics of land damaged by nonmetallic mines in Xinjiang. Combined with relevant standards and characteristics of nonmetallic mines in Xinjiang, the selection of evaluation indexes was focused mainly on four major aspects: the severity of land damage, the soil conditions, the original vegetation conditions of the land reclamation units, and the irrigation and drainage conditions. These are described as follows:(1)The severity of land damage includes the following aspects: (a) the dominant slope of the surface after excavation or the deposition of debris, which is denoted as X1. This refers to the post-excavation dominant slope or the post-deposit dominant slope of the land. The greater the slope is, the more severe the damage to the land is. (b) The depth of the excavation (m) or the flatness of the deposited debris (m3/m2), which is denoted as X2. This refers to the dominant depth of the ground relative to the surrounding ground after the excavation or the amount of leveling per unit area of the land after deposit. The greater the value of this index is, the more severe the land damage is, and the less suitable the land is for reclamation. (c) Yield reduction caused by the damaged land (%), which is denoted as X3. This refers to the decrease in productivity due to land damage, mainly referring to the land loss due to excavation in the mining area, and the reduction of productivity due to waste deposit and building construction. Similarly, the greater the value of X3, the more severe the damage to the land, and the less suitable the land for reclamation.(2)The soil conditions are mainly considered in four aspects, namely, clay content (%) X4, organic matter content (%) X5, soil thickness before damage (cm) X6, and soil conservation rate available for land reclamation (%) X7.(3)The original vegetation condition of the land reclamation units, that is, the vegetation coverage rate (%) of the land reclamation units before mining, X8. The more vegetation there is, the more the original ecological environment is conducive to vegetation growth, and the more favorable to land reclamation.(4)Irrigation and drainage conditions: according to the requirements of the standards, the irrigation and drainage conditions of the land reclamation units must be considered. Due to the dry climate in Xinjiang, water resources are extremely scarce, and the surface basically does not retain water. For this reason, drainage conditions are not considered and only irrigation conditions are considered in the selection of indexes. That is, the surface river runoff which is available for irrigation in the land reclamation units (m3/s) is denoted as X9.

3. Suitability Evaluation Units of Mines Land Reclamation

The evaluation unit is the basic physical unit for the suitability evaluation of land reclamation. The appropriateness of the evaluation unit has a direct effect on the quality of the land reclamation suitability evaluation, the size of the reclamation project, and the result of the reclamation. The nonmetallic mines in Xinjiang are mined in an open pit and the facilities of mines include open pits, waste rock yards, industrial squares, living quarters, blasting equipment storage areas, and landfills. Open-pit mining methods include side-hill cut mining and deep-concaved open-pit mining (Figure 2). For side-hill cut mining, the waste rocks cannot be backfilled, so that a high steep slope will be formed. The open pits formed in deep-concaved open-pit mining can be backfilled with waste rocks during or after the mining process, so the terrain is gentler. Hence, when dividing the open pits into land reclamation units, the open pits may be divided into just one land reclamation unit if the pits used deep-concaved open-pit mining. If the pits used side-hill cut mining, they would be divided into two reclamation units: an open-pit slope area and a flat bottom area, or they can be classified as one land reclamation unit of open mining slope area. If the mining method is a combination of side-hill cut mining and deep-concaved open-pit mining, then the open pit is divided into two land reclamation units: the open-pit slope area and the bottom pit mining area. The research conducted here is a study on 21 nonmetallic mines in Xinjiang, and the 21 mines are divided into 149 land reclamation units, as shown in Supplementary Table 1.

4. Classification Results of Land Reclamation Suitability Evaluation

4.1. Factor Analysis
4.1.1. Introduction of Factor Analysis

Factor analysis [3335] is a linear function of the common factor and the special factor that each original variable plays a dominant role in the data of the sample, so as to better explain the correlation of the original variable and reduce its dimension. The basic idea of factor analysis is to group variables according to the size of correlation, so that the correlation between variables in the same group is higher, and the correlation between variables in different groups is lower. Each group of variables represents a basic structure, which is represented by an unobservable comprehensive variable called a common factor. Another function of factor analysis is to classify variables (or samples) [36].

4.1.2. Results of Factor Analysis

The factor analysis data are listed in Supplementary Table 1. The factor analysis extracted the factor by the principal factor standard, and the criterion was that the first k principal factors were extracted when the cumulative contribution rate reached 85%. The load matrix of the factor was rotated according to the Kaiser-standardized orthogonal rotation method. After analysis, it was concluded that the measured KMO value was 0.781, the statistical value of the Bartlett spherical test was 1528.242, and the phase accompanying probability Sig. = 0.000 < 0.01, indicating that the chosen variables were suitable for factor analysis. The calculation process of factor analysis was carried out by SPSS software. The results of specific factor analysis are shown in Tables 2 and 3.

It can be seen from Table 2 that the selection of three factors could explain 85.91% of the variance of the original variables (accumulated contribution rate), indicating that all nine indexes variables could be well explained by the three principal factors. A matrix of rotational components (Table 3) could be obtained by rotating the extracted principal component matrix. Table 3 shows that the variables represented by the rotated principal factors were obviously polarized. The principal factor 1 could be interpreted as soil and vegetation factors; it could also be interpreted as the original ecological environmental conditions of the land reclamation units. Principal factor 2 represented the severity of mining damage of the land, and principal factor 3 represented the irrigation conditions. The three principal factors were consistent with the classification results of evaluation indexes, which showed that the selection of evaluation indexes was reasonable. The individual scores of the three principal factors (Supplementary Table 2) were calculated by SPSS software [37]. The composite score was equal to the individual scores of the three principal factors multiplied by the corresponding variance percentage. The degree of suitability for land reclamation (Supplementary Table 2) could be assessed from the magnitude of the composite score.

4.2. Cluster Analysis
4.2.1. Introduction of Cluster Analysis

Cluster analysis [38] is a statistical method for studying classification problem with the principle of multivariate statistical analysis. The idea of this method is to consider all samples as one class and then define the distance between classes. The two classes with the closest selectivity (e.g., the smallest distance) are merged into a new class. Then the distance between the new class and other classes is calculated, and then the two classes with the closest distance are merged, so that one class is reduced until all samples are merged into one class [36].

4.2.2. Results of Cluster Analysis

In cluster analysis, merging of categories is done using the average distance between categories and the distance between categories is measured by the square of the Euclidean distance. The calculation process of cluster analysis is carried out by SPSS software [37]. Figure 3 is the cluster analysis spectrum (tree diagram) and displays intuitively the gradual process of merging samples. A total of 148 steps were performed in this cluster. The land reclamation units No. 147 and No. 148 were clustered first which in turn was clustered with No. 149, and so the process went on. Finally, after 148 steps of clustering, 149 samples were grouped into one large category.

Figure 3 shows that it was most desirable to classify the 149 land reclamation units into five categories. In order to be consistent with the classification results from regulations and standards (Table 1), the land reclamation suitability evaluation results were also classified into four categories. Since the cluster analysis merely classifies the samples, the classification of the land reclamation suitability evaluation results must be done in combination with the composite scores from factor analysis (Supplementary Table 2). For the seven land reclamation units (unit designations 1∼7) of the limestone mine in the Aktas mining area in Chabuchar County, Xinjiang, category “1” has the highest composite score, and category “4” has the lowest composite score. According to the composite scores of the land reclamation units, categories IV and V of the cluster analysis were finally combined into one category. The results of land reclamation suitability evaluation are shown in Figure 3 and Supplementary Table 2. The results in Supplementary Table 2 show that the evaluation results using combined cluster and factor analysis were consistent with and more accurate than the standard classification results. Such evaluation results were more conducive to guiding the later stage of the land reclamation work.

4.3. Analysis of Classification Results

Through statistical analysis of the land reclamation suitability results of the above 149 land reclamation units, 46 land reclamation units had a reclamation suitability of “4,” which were unsuitable. The land reclamation direction was bare land, and the limiting factors of reclamation included extremely severe damage to the land, no irrigation conditions, poor soil quality, and undeveloped vegetation on the original surface. There were 35 land reclamation units whose reclamation suitability was “3,” which was not very suitable. The restrictions on land reclamation for these units included severe damage to the land, poor soil quality on the surface, and little vegetation on the surface. For land reclamation units with a suitability level of “very suitable” or “somewhat suitable,” the severity of land damage was generally light, the soil quality was also good, the surface vegetation was well developed, and the irrigation conditions were sufficient.

4.4. Analysis of Selected Evaluation Indexes

In order to analyze the correlation between the nine selected evaluation indexes and the land reclamation suitability level, the average of the nine classification indexes of 149 land reclamation units divided into 4 categories was calculated (Table 4 and Figure 4). Table 4 and Figure 4 show that, among the nine classification indexes, three indexes values (X1X3) gradually increased as the suitability of land reclamation became worse, and the values of the remaining six indexes gradually decreased as the suitability of land reclamation became worse. The changing behavior of the index values was consistent with factor analysis results and agreed with the actual situation of the land reclamation work. The selection of the nine classification indexes was reasonable.

5. Land Reclamation Suitability Evaluation Model

5.1. Establishment of Evaluation Model

The land reclamation suitability evaluation model described here adopted the Fisher discriminant analysis method. Fisher discriminant analysis method [39] is to establish a scientific classification rule, i.e., discriminant rule, to classify any new observation sample according to the known characteristics. The basic idea of Fisher discriminant is projection, which projects k-group m metadata to a certain direction, so that the post-projection group and the group can be opened as far as possible. Then, the appropriate discriminant rules are selected to classify and discriminate the new samples.

According to the principle and calculation steps of Fisher discriminant analysis method, nine evaluation indexes were used as independent variables of the discriminant analysis, and the suitability categories determined by combined cluster analysis and factor analysis were used as a grouping variable. The data of Supplementary Table 1 were analyzed using SPSS analysis software [37], and the land reclamation suitability evaluation model of nonmetallic mines reclamation in Xinjiang was established. The coefficients of the discriminant function were obtained, as shown in Table 5.

Using the coefficients of the Fisher discriminant function given in Table 5, the following three linear discriminant function expressions could be established:

The above three linear discriminant functions were three independent calculation formulas. The values of the corresponding sample functions could be calculated using functional formulas. To determine which group a new sample should belong to, the samples function values of the mines evaluation units were compared with the squared distance of the central value (Table 6) of the four classification groups. The sample belonged to the group which had the least squared distance.

5.2. Testing of Evaluation Model

According to the established Fisher discriminant model, the 149 land reclamation units were tested using a backward substitution estimation method. Only two misclassified samples appeared in the analysis, the correct rate was 98.7%, and the two misclassified land reclamation units were No. 12 and No. 139. According to the classification, No. 12 is the waste rock yard of the Dolomite Mine of Tianhu No. 1 in Hami City, Xinjiang, and No. 139 is the waste rock dumping site of the vermiculite mine of 32nd regimen. The classification result was “3” or “4” from relevant regulations and standards; according to the combined analysis of factors and clusters, the result was “4.” In terms of the severity of land damage, the land reclamation suitability should be “3.” However, from the viewpoint of irrigation conditions, the original vegetation development, and the soil quality, it is difficult to reclaim it into grassland. The land reclamation suitability rating is “4.” Figure 5 shows that level “3” and “4” samples are poorly clustered with the central points, with a poorly defined boundary. As a result, some of the samples are subject to misclassification.

5.3. Application of Evaluation Model

Using the land reclamation suitability evaluation model established above, the land reclamation units of the Toga Solo limestone mine in Dabancheng, Urumqi, were evaluated for land reclamation suitability. The limestone mine had an area of 1.5576 km2, with a total of 5 ore bodies. Among them, the No. 2 and No. 5 ore bodies adopted the side-hill cut mining, and the No.1, No. 3, and No. 4 ore bodies adopted the side-hill cut mining and deep-concaved open-pit mining method. The evaluation divided the entire mining area into 15 land reclamation units (Figure 6). The specifically measured indexes of each land reclamation unit are shown in Table 7.

The indexes values of the 15 land reclamation units in Table 7 were substituted, respectively, into three discriminant formula, and each land reclamation unit can obtain three discriminant function values. The squared distance from the discriminant function value to the central value of each category was calculated for each evaluation unit and compared. The results are shown in Table 8. It could be seen from Table 8 that the classification results of the land reclamation suitability of the 15 land reclamation units were consistent with the classification results from regulations and standards. Furthermore, the classification results were unique, the conclusions were more accurate, and the classification could more favorably guide the task of reclaiming the mine land.

6. Conclusions

(1)In this work, the evaluation indexes system of land reclamation suitability evaluation for nonmetallic mines in Xinjiang was established. The evaluation indexes fit the actual situation of nonmetallic mines in Xinjiang [21, 22]. They have included the factors affecting the suitability of land reclamation and the direction of land reclamation to the greatest extent possible. The evaluation results are highly consistent with the actual situation and the analysis of the classification indexes data showed that all nine indexes selected have affected the land reclamation suitability evaluation results to varying degrees.(2)According to the mining mode of nonmetallic mines and the layout of mine facilities, the mines are reasonably divided into several land reclamation units. For open pits mined with the side-hill cut mining method, the pits are divided into the slope region and the bottom flat area. Open pits mined with deep-concaved open-pit mining method are divided into just one open-pit mining area. Open pits mined with the combined side-hill cut and deep-concaved open-pit mining method are divided into slope region and deep-concaved open-pit mining part.(3)The 149 land reclamation units are classified into 4 categories using a combination of factor analysis and cluster analysis. The classification results are compared with that obtained using the regulations and standards; the comparison showed that the classification was unique and consistent with the regulations and standards.(4)In the evaluation model established here, all nine indexes can be obtained quantitatively through actual measurement, experiment, and relevant reports, without resorting to weight factors and human-assigned values [68, 27, 28]. It avoided the effects of human factors on the results of the evaluation model. In actual application of the evaluation model, it is only necessary to collect the relevant classification indexes data.(5)In this work, the discriminant analysis method was applied to the suitability evaluation of land reclamation for mines, and a land reclamation suitability evaluation model was established for land reclamation of nonmetallic mines in Xinjiang. The model was tested and the results of real case calculation proved that the calculation process of the discriminant analysis method was simple to implement and the model structure was stable. The model was shown to have superior discrimination capability for suitability evaluation results of mines reclamation and the evaluation model may serve as a useful reference in the effort of reclaiming for nonmetallic mines in Xinjiang.(6)The biggest difference between the classification results of this model and standard classification results is to distinguish the classification of “not very suitable” and “unsuitable,” which is often the most critical difference affecting the mines land reclamation direction in Xinjiang.(7)The model is established on the basis of analyzing 149 land reclamation units in 21 mines. It is representative, but it still has certain limitations due to the limited samples. In future work, the evaluation model may be further refined by including even more representative nonmetallic mines.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the Natural Science Foundation of Xinjiang (grant no. 2018D01C061).

Supplementary Materials

Samples of suitability evaluation of nonmetallic mine land reclamation in Xinjiang. (Supplementary Materials)