2.2.1. Ecological Sensitivity Assessment
- (1)
Indicator Selection
Drawing insights from the relevant literature, specific factors related to terrain, soil, and climate were identified [
13,
14,
22,
23,
24,
25,
26,
27,
28,
29]. The intermediate layer of abandoned mining land is divided into topographical conditions, surface water systems, and vegetation landscapes. On one hand, the topography of abandoned mining land is unique, requiring the consideration of factors such as elevation, slope, aspect, land use, and other evaluation factors. On the other hand, abandoned mining land experiences severe ecological damage, with vegetation being the most affected. Therefore, it is necessary to investigate factors such as vegetation coverage and vegetation types. In the research area, the prominent feature of abandoned mining land is water accumulation forming a lake, which makes the surface water system a decisive factor. Initial screening resulted in 56 assessment factors across terrain conditions (
Table 1), surface water systems (
Table 2), and plant landscapes (
Table 3) [
30]. After eliminating duplicates, 9 factors remained. Applying the Delphi method for questionnaire selection, quantified scoring was conducted on the factors influencing the initial selection set. The extracted influencing factors were then subjected to secondary selection. The process involved integration, refinement, and other enhancement steps for comprehensive improvement. Considering the current site conditions, the final factors closely related to ecological sensitivity suitability were identified as elevation, slope, aspect, soil texture, land use, runoff buffer zone, water body buffer zone, vegetation coverage, and vegetation type (
Figure 3).
- (2)
Determination of Indicator Weights
- (1)
Determination of Weights by Individual Experts
Applying the Delphi method, a combination of online and offline questionnaires was administered to experts with postgraduate qualifications or above in fields such as landscape ecology, landscape architecture, design, urban and rural planning, and architecture. The questionnaire theme was “Analysis of the Importance of Ecological Sensitivity in Abandoned Mining Land.” Following a pairwise comparison approach, questions were framed such as “Compared to ecological sensitivity in abandoned mining land, which is more important: topographical conditions or surface water systems?” The importance levels were rated on a scale from 1 to 9, where a higher number indicated a higher degree of importance (
Table 2).
- (2)
Analytic Hierarchy Process (AHP)
The obtained expert index scores were imported into the Yaahp V11.3 software to construct a hierarchical structure model. The m indicators to be analyzed were then transformed into a judgment matrix
A based on the ratings provided by each expert (
Table 1).
The elements in set
A satisfy the following:
Constructing a judgment matrix involves pairwise comparisons between various elements and determining the weights of each criterion layer on the target layer. In simple terms, this means evaluating the indicators of the criterion layer via pairwise comparisons, typically using Saaty’s 1–9 scale method.
- B.
Hierarchical single sorting
Hierarchical single sorting refers to comparing each element in the current layer pairwise with all other elements in the same layer, conducting hierarchical ranking, and arranging the order of importance. The specific calculations can be based on the judgment matrix A, ensuring that it satisfies the conditions of the characteristic roots and eigenvectors of AW. Here, the maximum characteristic root of A is denoted as λmax, and the normalized eigenvector corresponding to λmax is represented as W. Each component of W represents the weight and corresponds to the single sorting of the respective element. Utilizing the judgment matrix, the weights (coefficients) of each factor on the target layer can be calculated. The calculation steps for the weight vector (W) and the maximum characteristic root (λmax) using the root method or summation method are shown:
Multiply the elements along each row, and then take the
nth root.
Normalizing
(making the sum of its elements equal to 1) results in the ranking weight vector. Denote this normalized vector as
W (where the elements of
W represent the relative importance ranking weights of factors in the same hierarchy concerning a factor in the previous hierarchy). Therefore,
= (
W1,
W2, …,
Wn).
- C.
Calculating the maximum characteristic root and CI value. Based on the CI and RI values, calculate the CR value to determine whether the consistency is acceptable.
Calculate the maximum eigenvalue of the matrix:
In the formula,
is the order of the judgment matrix, and
represents the weight of each factor (
i = 1, 2, …,
n).
In the formula, CI is the consistency index. CI = 0 indicates complete consistency in the judgment matrix. The larger the CI, the more severe the inconsistency in the judgment matrix; CR is the consistency ratio; RI is the random index; the data can be obtained via the random index (R.I.) value table (
Table A1). Obtaining the consistency test result CR < 0.1 indicates that the judgment matrix has satisfactory consistency. If CR is greater than or equal to 0.1, it is advisable to consider making adjustments to the judgment matrix
A.
In the formula, represents the normalized subjective assessment weight values.
- (3)
Classification of Evaluation Indicators
The impact of different indicator factors on the ecological environment varies, referring to the fact that different levels and types have different degrees of impact on the ecological environment. Therefore, refining the indicator factors based on the general phenological characteristics and objective principles of each factor, and combining the characteristics of the study area, a hierarchical assignment is carried out. Ecological sensitivity is divided into 5 levels, and the sensitivity levels of each indicator factor are divided into 9, 7, 5, 3, and 1, corresponding to highly sensitive, moderately sensitive, sensitive, less sensitive, and insensitive ecological sensitivities, respectively. Ecological suitability is categorized as most unsuitable, less unsuitable, suitable, less suitable, and most suitable. Via a literature review, field surveys, expert consultations, etc., the definition and criteria for the division of each indicator factor are clarified (
Table 3).
Soil texture refers to the particle composition in the soil, closely related to soil permeability and nutrient content. The growth and development of vegetation are influenced by soil texture. The higher the nutrient content and permeability of the soil, the more luxuriant the vegetation, making it less suitable for development, and the ecological sensitivity is higher [
31]. This article primarily refers to the “Green Mining Construction Standard for Non-metallic Mines” (DZ/T 0312-2018) [
32] and the “Technical Guidelines for Environmental Impact Assessment—Soil Environment (Trial)” (HJ 964-2018) [
33]. It explores the efficacy of different soil textures in serving ecological restoration and evaluates and classifies them accordingly. The soil is categorized as Level I Soil Texture Sensitivity, sandy loam as Level II Soil Texture Sensitivity, clay as Level III Soil Texture Sensitivity, heavy clay as Level IV Soil Texture Sensitivity, and sandy soil as Level V Soil Texture Sensitivity. The slope represents the inclination angle of the ground. The steeper the slope, the more difficult it is for vegetation to recover after damage, making it less suitable for development, and ecological sensitivity is higher. According to the “Technical Regulations for Forest Resource Planning and Design Survey” (GB-T 26424-2010) [
34], different slope intervals are classified as follows: flat slope (slope range 0°–5°), gentle slope (slope range 6°–15°), moderate slope (slope range 16°–25°), steep slope (slope range 26°–35°), and steep incline (slope range 36°–45°). Based on the current situation of the study area, the steepest slope is 35.2°, and the gentlest slope is 8.8°. Based on the current conditions in the study area, with the steepest slope being 35.2° and the gentlest slope being 8.8°, an equal interval reclassification was conducted with adjusted intervals, essentially meeting the requirements of the “Technical Regulations for Forest Resource Planning and Design Survey”. Slopes above 35.2° are classified as Slope Sensitivity Level I, 26.4°–35.2° as Slope Sensitivity Level II, 17.6°–26.4° as Slope Sensitivity Level III, below 8.8° as Slope Sensitivity Level V, and 8.8°–17.6° as Slope Sensitivity Level IV. Land-use type is classified based on the current state of the land, determining the difficulty of land development and affecting the ability of ecological restoration [
35]. The higher the land-use grade, the greater the response of the land ecosystem to human activities and natural environmental changes, making it less suitable for development, and ecological sensitivity is higher [
35,
36]. According to the “Land Quality Classification for Cultivated Land” (GB/T 33469-2016) [
37], forests are classified as Land Sensitivity Zone I, water bodies as Land Sensitivity Zone II, grasslands as Land Sensitivity Zone III, cultivated land as Land Sensitivity Zone IV, and construction land and bare land as Land Sensitivity Zone V. Elevation refers to the height above sea level. Higher elevations are associated with lower biodiversity, indicating higher ecological sensitivity [
38]. Based on an elevation difference of 109.5 m in the study area, using equal interval reclassification, with the highest altitude being 93 m and the lowest being −16.5 m, a reclassification with adjusted intervals was conducted. Elevations above 71.2 m up to 93 m are classified as Altitude Sensitivity Level I, 49.4–71.2 m as Altitude Sensitivity Level II, 27.6–49.4 m as Altitude Sensitivity Level III, −16.5 m to 5.8 m as Altitude Sensitivity Level IV, and 5.8–27.6 m as Altitude Sensitivity Level V. The aspect refers to the direction a slope faces, and it influences the reception of solar radiation on the ground, indirectly affecting the photosynthesis of plants [
39]. The study area is located north of the Tropic of Cancer, and throughout the year, the noon sunlight predominantly shines on the south-facing slopes, with the north-facing slopes receiving less solar radiation. Therefore, the Aspect I sensitive area is situated in the northern slope region, the Aspect II sensitive area is in the northwest and northeast slope regions, the Aspect III sensitive area is in the east and west slope regions, the Aspect IV sensitive area is in the southeast, and southwest slope regions and the Aspect V sensitive area is on the south-facing slopes and flat terrain.
Due to the unique terrain of the mining pit, precipitation gathers to form site runoff, which can serve as a crucial water resource reserve for future ecological restoration [
40,
41]. Utilizing a buffer zone analysis based on Euclidean distance measurement, the water body itself is classified as extremely highly sensitive. The closer the buffer zone is to the water body, the less suitable it is for development, indicating higher ecological sensitivity. Therefore, employing the equal interval method, the water body buffer zone is divided into Sensitivity Level I for distances less than 10 m, Sensitivity Level II for distances between 10 m and 30 m, Sensitivity Level III for distances between 30 m and 50 m, Sensitivity Level IV for distances between 50 m and 70 m, and Sensitivity Level V for distances above 70 m. Surface runoff, as an important means of water supply, impacts the growth and development environment of plants. When surface runoff increases, it can erode the soil environment essential for plant survival, which is unfavorable for plant growth. Therefore, after consulting relevant research and seeking advice from an expert in soil and water conservation (at the associate professor level), the sensitivity classification of surface runoff was conducted based on an understanding of its impact range. This classification was carried out in a manner that divides it into suitable intervals [
40,
41,
42,
43,
44]. Among them, the 5–10 m range is classified as a Level I runoff buffer zone sensitivity area, 0–5 m and 10–15 m are Level II runoff buffer zone sensitivity areas, 15–20 m is Level III runoff buffer zone sensitivity area, 20–30 m is Level IV runoff buffer zone sensitivity area, and above 30 m is Level V runoff buffer zone sensitivity area.
Vegetation type refers to the plant communities covering a specific region on the Earth’s surface, categorized based on different plant communities. Vegetation type is a common qualitative factor that reflects the ecological quality of site vegetation communities. Generally, mixed forests with a multi-layered structure exhibit higher ecological service capabilities and stability compared to ordinary forests, shrublands, grasslands, etc., making them relatively more sensitive. Therefore, mixed forests are classified as Vegetation Type Sensitivity Level I, broadleaf pure forests and coniferous pure forests are classified as Vegetation Type Sensitivity Level II, shrub forests are classified as Vegetation Type Sensitivity Level III, crops and grasslands are classified as Vegetation Type Sensitivity Level IV, and others are classified as Vegetation Type Sensitivity Level V. Vegetation coverage indicates the percentage of the ground area occupied by vegetation (including leaves, stems, branches) about the total area of the surveyed region. Vegetation coverage can be categorized into four types: high, moderately high, moderate, and low. The higher the vegetation coverage and the richer the vegetation types in a certain area, the better the geographical condition of the land unit in that region. Consequently, it is less suitable for development, and the ecological sensitivity is higher [
45,
46]. Based on a 30% difference in vegetation coverage in the study area, an equal interval reclassification was applied. The highest vegetation coverage is 30%, and the lowest is no vegetation coverage. With equal interval reclassification, areas with vegetation coverage above 30% are classified as Vegetation Coverage Sensitivity Level I, 20–30% as Vegetation Coverage Sensitivity Level II, 10–20% as Vegetation Coverage Sensitivity Level III, 0–10% as Vegetation Coverage Sensitivity Level IV, and areas with no vegetation coverage as Vegetation Coverage Sensitivity Level V.
2.2.3. Data Sources and Processing
- (1)
Data sources
To make the obtained conclusions more convincing, the selected data sources were all the latest ones available at the time of the study. Elevation, slope, aspect, and runoff factor data for the ecological sensitivity assessment system of abandoned mining land were derived from the DEM (Digital Elevation Model) of the Geographic Spatial Data Cloud (
http://gscloud.cn, accessed on 23 March 2023); water body data were obtained using Baidu API crawling in 23 March 2023 (
https://lbsyun.baidu.com/, accessed on 23 March 2023); soil texture data were sourced from the 30 m precision soil type dataset of the Chinese Academy of Sciences in 2022 (
http://www.csdn.store/, accessed on 23 March 2023); land-use data were derived from the 30 m precision land-use dataset of the Chinese Academy of Sciences in 2022, manually calibrated (
https://www.resdc.cn/, accessed on 23 March 2023); vegetation type data were sourced from the 1:1,000,000 vegetation dataset (
http://www.ncdc.ac.cn/, accessed on 23 March 2023); and vegetation cover data were calculated from the 30 m precision NDVI (Normalized Difference Vegetation Index) in 2022 (
http://www.ncdc.ac.cn/, accessed on 23 March 2023).
The landscape pattern indices, including NP, CA, PLAND, LPI, TE, ED, SHAPE_MD, AREA_MN, CONTAG, and PR, were calculated based on the reclassified land-use types from the Chinese Academy of Sciences in 2022 (
https://www.resdc.cn/, assessed in 23 March 2023). These indices were imported into Fragstats 4.2 for computation.
- (2)
Data Processing
Using GIS reclassification tools to obtain suitability ratings for individual indicators, a comprehensive ecological sensitivity analysis for the scenic area was conducted using a multi-criteria weighted overlay method [
54]. The evaluation model is as follows:
In the formula, represents the comprehensive score for a specific indicator; is the weight of the -th indicator; is the score for the -th indicator; and is the number of indicators.