Next Article in Journal
Will Customers’ Understanding of the Trolley Dilemma Hinder Their Adoption of Robotaxi?
Previous Article in Journal
Manufacturing Execution System Application within Manufacturing Small–Medium Enterprises towards Key Performance Indicators Development and Their Implementation in the Production Line
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Incorporating the Results of Geological Disaster Ecological Risk Assessment into Spatial Policies for Ecological Functional Areas: Practice in the Qilian Mountains of China

Chengdu University of Technology, Chengdu 610059, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(7), 2976; https://doi.org/10.3390/su16072976
Submission received: 1 February 2024 / Revised: 16 March 2024 / Accepted: 29 March 2024 / Published: 3 April 2024

Abstract

:
Geological hazards cause changes in the quality of the ecological environment, affect the function and stability of ecosystems, and negatively impact the maintenance and restoration of ecological functions in ecological functional areas (EFAs). This study integrates machine learning, geographic information technology, and multivariate statistical analysis modeling to develop a technical framework for quantitative analysis of ecological risk assessment (ERA) based on the causal logic between geological hazards and ecosystems. The results of the geological disaster ERA are mapped to EFAs, effectively identifying and quantifying the risk characteristics of different EFAs. The results show that: (1) The hazard–vulnerability–exposure ERA framework effectively identifies the distribution characteristics of high ecological risk around the Qilian Mountains, with high risk in the east and low risk in the west. (2) In high ecological risk areas, high hazard–high vulnerability–low exposure is the main combination pattern, accounting for 83.3%. (3) Overall, hazard and vulnerability have a greater impact on geological disaster ecological risk than exposure, with path coefficients of 0.802 (significant at p = 0.01 level) and 0.438 (significant at p = 0.05 level), respectively, in SEM. The random forest model (R2 = 0.748) shows that social factors such as human density and road density contribute significantly more to extreme high risk than other factors, with a contribution rate of up to 44%. (4) Thirty-five ecological functional units were systematically grouped into four clusters and used to formulate a “layered” spatial policy for EFAs. The results of the research are expected to provide support for maximizing the policy impact of EFAs and formulating management decisions that serve ecological protection.

1. Introduction

Ecological risk assessment (ERA) is a methodology for assessing and predicting the extent of ecosystem loss and degradation caused by human activities or natural hazards, in order to determine the state of ecosystem quality and the level of risk [1,2,3]. ERA provides an integrated assessment of human activities and environmental systems by identifying and linking the possible impacts of risks over long time periods and large regional scales in environmental management [4,5,6]. Ultimately, ecological environmental risk management and environmental monitoring provide decision support. ERA of geohazards is one of the hotspots of current research [7,8,9], focusing on the impacts and extent of damage to the ecological environment caused by geohazards [10,11,12,13] and the impacts and threats to ecosystem stability and function [14,15,16]. The main research contents of geohazard ERA include hazard assessment, exposure assessment, receptor analysis, and risk characterization [17,18,19,20]. ERA of geohazards is mostly performed by selecting multiple risk sources [21,22,23] and using ecological index and vulnerability index as assessment indicators [24,25,26,27], but studies on ERA of multiple risk sources at landscape scale are still relatively rare.
Ecological function areas (EFAs) is a type of geospatial zoning based on ecosystem characteristics, stress processes and effects, the importance of ecological service functions, and ecological sensitivity [28]. It is the basis and prerequisite for the implementation of regional ecological environment zoning management, which aims to clarify the important areas of regional ecological security and key areas for protection, analyze the existence of ecological problems and vulnerable areas, and provide a scientific basis for ecological protection and construction planning. It is the basis and prerequisite for the implementation of regional ecological environment zoning management [29]. There is a close relationship between ERA and EFAs, where the receptor of ERA is usually the ecosystem service function, and EFAs are areas whose main function is to provide ecological products [30,31]. We found that current research results on ERA focus more on the exposure of ecosystem services as a whole, and lack the differentiation of ecosystem services in terms of type and space, which makes it difficult to use the risk assessment results in the management and policy formulation of EFAs.
The Qilian Mountain region is important for China’s ecological security and protection. However, there are prominent problems in the maintenance and restoration of the regional ecological function of the Qilian Mountains due to geological disasters and anthropogenic factors. Therefore, there is a current need to understand how to more effectively use the EFAs for policy impact and the development of services for ecological protection and management to protect the stability of the ecosystem services in the Qilian Mountain ecosystem.
Therefore, the objectives of this study are: (1) to develop a framework for geological disaster ERA; (2) to analyze the spatial distribution, patterns, and mechanisms of geological disaster ecological risk in the Qilian Mountains region of China; (3) to integrate the results of geological disaster ERA into policies for EFAs. The aim is to provide scientific references for the ecological protection and management of the Qilian Mountains region.

2. Materials and Methods

2.1. Study Area

In this study, Qilian Mountain was selected as the study area, the region is one of China’s national nature reserves and is a world-class wetland ecological reserve [32] (Figure 1), which is an important place for plateau birds and fowl to live and reproduce, and at the same time, it is the only habitat for the endangered Przewalski’s gazelle [33]. On the other hand, various geological disasters such as landslides and rockfalls frequently hit the area [34,35], posing a serious threat to the ecological security of the region. The Qilian Mountains serve as an important ecological barrier and water conservation area in the northern Qinghai–Tibet Plateau [36,37]. For the ecological protection and restoration of the ecological barrier area in northern Qinghai, it is urgent to propose spatial policies to address the ecological risks of geological hazards and curb the trend of ecological degradation [38,39].

2.2. Data Sources and Data Pre-Processing

The main data sources in the research area are (1) compilation maps of historical collapses and landslides, as well as relevant field survey data; (2) digital elevation model (DEM) obtained from the Geographic Spatial Data Cloud (http://www.gscloud.cn/ (accessed on 10 January 2022)), mainly used for acquiring terrain and landform basic environmental factors such as elevation, slope, and aspect; (3) rock type data obtained from the Chinese Academy of Sciences Resource and Environment Science Data Center (http://www.resdc.cn (accessed on 10 January 2022)); (4) Landsat-8 remote sensing imagery used to obtain land use data and the normalized difference vegetation index (NDVI); (5) data on rivers, roads, population density, and fault zones obtained from the National Geographic Information Resource Catalog Service System (https://www.webmap.cn (accessed on 10 January 2022)).

2.3. Methodology

Based on the causal logic between external disturbances and ecosystems, a framework for ERA was constructed (Figure 2). Hazard reflects the probability and severity of geological hazards, while vulnerability reflects the ability of landscape patterns to withstand geological hazards. Exposure reflects the sensitivity of ecosystem functions to geological hazards. Based on the evaluation of hazards, vulnerability, and ecosystem function values, the calculation formulas are as follows [40]:
E R = H × V × E
where E R is ecological risk, H is hazard, V  is vulnerability, and E is exposure.

2.3.1. Hazard Calculation

Step 1: First, we selected the common driving factors for different subclasses of geological hazards, including elevation, slope, aspect, terrain wetness index, lithology, soil type, fault, annual precipitation, temperature, NDVI, roads, population density, and 12 other factors (Figure 3). After resampling, the spatial resolution of all factor data was set to 1000 m. Different factor layers were stacked in a specific order to form multi-channel raster data.
Step 2: Random non-geological hazard points were generated outside the 5–10 km buffer distance from geological hazard points to construct data point samples. Geological hazard and non-geological hazard block datasets were then created based on the multi-channel raster data overlaid in Step 1. In addition, the training and test samples were split in a 7:3 ratio.
Step 3: The entire sample set was fed into the SVM model. The radial basis function (polynomial kernel) was used to evaluate the risk of geological hazards. The grid search algorithm was used to optimize the parameters and find the optimal gamma and penalty factor C [41]. The SVM formula is as follows:
y i ( w T x i + b ) 1
max 1 w
Step 4: The receiver operating characteristic (ROC) curve is a comprehensive indicator of sensitivity and specificity as continuous variables [42,43]. In this study, it was chosen to assess the accuracy of the predictions.

2.3.2. Calculation of Ecosystem Vulnerability

Ecosystem vulnerability refers to the sensitivity and resilience of an ecosystem to external disturbance. From the perspective of the landscape background, the degree of disturbance of landscape types, the landscape vulnerability index, and the landscape adaptive capacity index can reflect the vulnerability of ecosystems. The ecosystem vulnerability index is as follows [40]:
L V I = U × V / ( 1 + L A I )
where L V I represents ecosystem vulnerability, U represents the degree of disturbance of landscape types, V represents landscape vulnerability, and L A I represents the landscape adaptability index (Table 1).

2.3.3. Calculation of Ecosystem Exposure

In this study, the value of ecosystem services is considered as the receptor of ecological risk. Based on the equivalent table of the value of ecosystem services per unit area in China estimated by Xie et al. (2015) [44] through a survey of 500 ecologists, it is adjusted in combination with the land use types in the Qilian Mountains. Therefore, the ecosystem service value formula is as follows:
E i = j = 1 9 k i j A j
E i represents the value of ecosystem services for a certain landscape type i . k i j represents the value-weighted factor of j types of ecosystem service functions for landscape type i . A j  represents the area of landscape type i.

3. Results

3.1. Hazard Prediction Results

There were 2252 pixels in the training set and 965 pixels in the test set. The results show that when gamma is 0.5 and C is 5, the AUC value is relatively high with minimal loss of accuracy, making the model optimal (Figure 3). As shown in Figure 4, the AUC value for disaster assessment based on the SVM model reached 0.8898.

3.2. Results of ERA

As shown in Figure 5a, the calculation results of the ERA are presented using a color gradient from red to blue, where red represents a higher ecological risk and blue represents a lower ecological risk. It can be seen that the red areas are concentrated in the eastern region and scattered in the western and peripheral areas. The forests and grasslands of the Qilian Mountains are mainly distributed in the east, where ecosystem services are concentrated, and the ecological service values are generally higher. In contrast, the western region is dominated by large areas of bare land, sandy areas, and deserts, resulting in sparse and fragmented ecological resources and relatively lower ecosystem service values compared to the east. Therefore, we consider that the eastern part of the Qilian Mountains has a higher risk from human activities, coupled with the concentration of ecological functions in forests and grasslands, resulting in a high aggregation of ecological risk. Meanwhile, the western part of the Qilian Mountains has a dispersed distribution of potential ecological risk due to the relative dispersion of ecological sources and the high vulnerability of sandy and desert areas.
Taking the example of the eastern region with the highest risk, we analyze the indicators. The highest-risk areas are mainly in the low-lying river valleys of the eastern region, where vulnerability and exposure are significantly higher than in other areas. This area has a higher population density and more urbanization, with significant impacts from human activities such as deforestation, overgrazing, logging for fuelwood, and mining. Therefore, Figure 5e–g show higher values for landscape fragmentation, landscape fractal dimension, and isolation in this region. In addition, the predominant landscape types in the east are forests and grasslands, which generally have higher ecosystem service values, leading to higher vulnerability in this area.

3.3. Spatial Patterns of Ecological Risk

We used spatial analysis to discuss the spatial autocorrelation patterns of ecological risk and its assessment indicators. On a global scale, the Moran indices of ecological risk and hazard are 0.288 (significant at the p = 0.05 level) and 0.287 (significant at the p = 0.05 level), respectively, both showing an agglomeration distribution pattern. The Moran index of vulnerability and exposure are 0.079 and 0.124, respectively, but do not pass the significance test and both show a random distribution pattern. Figure 6a–d shows the local spatial autocorrelation patterns of various indicators of ecological risk, where ecological risk shows a high-value agglomeration in the eastern part of the Qilian Mountains and a low-value agglomeration in the central part. The spatial distribution of hazard is consistent with that of ecological risk. Vulnerability shows a high-value agglomeration in the eastern part of the Qilian Mountains and is not significant in other areas. Exposure shows a high-value agglomeration in the southern part of the Qilian Mountains and a low-value agglomeration in the western part.
Following the natural breakpoint method, the ecological risk and its three-dimensional assessment indicators are classified into low risk, medium risk, and high risk. With the ecological risk classification as the base map, the classification combination patterns of hazard, vulnerability, and exposure are shown. From Figure 6e, there are mainly six combination patterns: low hazard–low vulnerability–low exposure, low hazard–low vulnerability–high exposure, low hazard–high vulnerability–low exposure, low hazard–high vulnerability–high exposure, high hazard–low vulnerability–low exposure, and high hazard–high vulnerability–low exposure. Table 2 summarizes the spatial combination patterns distributed in different levels of ecological risk. In low ecological risk, low hazard–low vulnerability–low exposure is the main combination pattern, accounting for 50%. In medium ecological risk, high hazard–high vulnerability–low exposure is the main combination pattern, accounting for 46.2%, followed by low hazard–high vulnerability–low exposure, accounting for 30.8%. At high ecological risk, high hazard–high vulnerability–low exposure is the main combination pattern, accounting for 83.3%.

3.4. Sources of Ecological Risk

We performed a linear fit of the risk values of ecosystem functions with geological disaster hazard, vulnerability, and exposure, respectively (Figure 7a), with R2 values of 0.37, 0.19, and 0.17, respectively, passing the significance test. This indicates that ecological risk is positively correlated with geological disaster hazard, vulnerability, and exposure overall, and is more affected by geological disaster hazard. To further reveal the causal mechanism of geological disaster ecological risk in the Qilian Mountains, we selected population density and road density as social factors and elevation and slope as natural environmental factors, combined with the three indicators of ecological risk to construct an SEM model of influencing factors and ecological risk. As shown in Figure 7c, the path coefficients of hazard, vulnerability, and exposure with ecological risk are 0.802 (significant at the p = 0.01 level), 0.438 (significant at the p = 0.05 level), and 0.286 (significant at the p = 0.01 level), respectively. From the path coefficients of social factors and natural environmental factors with hazard, vulnerability, and exposure, the path coefficient of hazard with social factors is higher than that with natural factors, indicating that geological disasters in the Qilian Mountains are greatly affected by human activities. The path coefficient of vulnerability with social factors is 5.275 (significant at the p = 0.1 level), while its path coefficient with natural environmental factors is 1.269 (significant at the p = 0.1 level), indicating that vulnerability is jointly affected by human activities and natural environmental factors. The path coefficient of exposure with social factors is also higher than that with natural factors. Therefore, we believe that geological disaster hazards affected by human activities are the main source of ecological risk, followed by vulnerability distribution affected by both human activities and the natural environment as a secondary source of ecological risk.
We divided the geological disaster ecological risk into five levels: extremely low, low, medium, high, and extremely high. Using the random forest model, we calculated the importance of influencing factors at each level. From Figure 7d–h, it can be seen that in the extremely low ecological risk zone, the contribution of natural environmental factors and vulnerability is relatively high, both exceeding 18%. In low-risk zones, the contribution rate of vulnerability is significantly higher than other factors, exceeding 40%. In medium-risk zones, the contribution of hazard is the highest. In high- and extremely high-risk zones, the contribution of social factors such as human density and road density is significantly higher than other factors. The importance of factors calculated by the random forest model is consistent with the results of the SEM model.

4. Discussion

4.1. Spatial Policies for Addressing Ecological Risks of Geological Disasters

Based on the calculation results of ecological risk levels, spatial combination patterns, and source factors in the Qilian Mountains, a total of 35 functional units were grouped into four clusters using a systematic clustering method to formulate policy zones for managing geological hazard risks. Cluster 1 includes 15 ecological functional units located in the outer ring of the Qilian Mountains, which mainly function in water conservation and biodiversity. Cluster 2, located at the outermost periphery, consists of six ecological functional units, mainly focusing on soil conservation. Cluster 3 consists of 13 ecological functional units located in the eastern and central parts, which are mainly dedicated to biodiversity protection. Cluster 4 consists of one ecological functional unit. The distribution of ecological risk results across the different clusters shows a clear concentric structure (Figure 8). Cluster 4, which is the innermost cluster, has the lowest risk, followed by clusters 3 and 1, while cluster 2, which is the outermost cluster, has the highest ecological risk.
Cluster 2 areas are characterized by high risk levels and high exposure, making them priority areas for hazard prevention and management. Since this region is also a primary area for human activities, we suggest strengthening geological hazard monitoring [35] and enhancing landscape restoration, coupled with conducting environmental impact assessments for land use planning.
Compared to Cluster 2, Cluster 1 exhibits relatively higher average vulnerability but lower hazard and exposure levels. Additionally, this region shows a belt-like distribution, where potential ecological risks from geological hazards directly impact ecological processes such as species migration and population distribution among ecological source areas [45]. In geological hazard prevention and control, enhancing ecosystem resilience can be achieved through measures such as improving vegetation coverage and diversity [46].
In cluster 3 areas, geological hazard ecological risks primarily stem from natural environments. Compared to Clusters 1 and 2, these areas experience relatively low precipitation, sparse populations, economic underdevelopment, fewer objects threatened by disasters, and predominantly desert landscapes. Focus in this region should be on protecting mountain vegetation, controlling soil erosion, enhancing local soil and water conservation capacities, and strengthening geological hazard risk monitoring for important species habitats.
Cluster 4 areas exhibit relatively low potential ecological risks from geological hazards. Due to the low geological hazard risks and good ecological resilience foundation, prevention of ecological risks takes precedence in this region.

4.2. Advantages of the Technical Framework Proposed in This Study

Natural disasters can cause changes in the quality of the ecological environment, resulting in some damage to ecosystems [47]. The management of ecosystems under disaster risk pressure has become an increasingly pressing issue [48,49]. Effective maintenance of existing healthy ecosystems and comprehensive management and restoration of degraded ecosystems have become critical issues that require urgent attention [50,51]. In this study, we developed a quantitative analytical framework to assess regional ecological security based on the causal logic between geological hazard disturbances and ecosystems. We assigned the results of geological hazard ERAs to EFAs, thereby effectively identifying and quantifying the risk characteristics of different types of EFAs, and providing scientifically based optimization solutions. Compared to previous studies (2, 3, 31), the technical framework proposed in this study establishes links between ERA and EFAs (Figure 9), which can provide references for decision-making in ecosystem protection and management of EFAs under geological hazard pressure.

4.3. Limitations and Prospects

This study establishes a technical framework for ecosystem risk management by linking geological disaster ERA with EFAs, which can help to better understand the safety status of EFAs and formulate targeted restoration and repair plans according to the degree of risk, thus providing a scientific basis for geological disaster prevention and ecological restoration. However, this study also has some limitations. Disaster risk is mainly focused on geological disaster risk. In the future, multiple indicators and complex models can be integrated to introduce more risk indicators related to ecological security, such as drought and flood, and to consider more factors related to risk management, such as climate change and socio-economics. More complex models that truly reflect the dynamic changes and mechanisms of ecosystems under risk disturbance can be used.

5. Conclusions

In this study, we conducted a comprehensive assessment of ecological risks associated with geological disasters in the Qilian Mountains region. Using a support vector machine model for hazard prediction, we determined optimal parameter configurations and delineated the spatial distribution patterns of ecological risks using geographic information systems. Using structural equation modeling and spatial analysis, we elucidated the mechanisms through which factors such as geological hazards, vulnerability, and exposure influence ecological risks. We also proposed spatial policies and management strategies tailored to different risk zones and discussed the advantages and limitations of our technical framework. Our research provides a scientific basis for ecosystem protection and management while advocating for future studies to deepen our understanding of ecological security through multi-factor, multi-model approaches to better address environmental challenges.

Author Contributions

X.L.: writing—original draft preparation, writing—review and editing, and project administration. Q.X.: methodology, software, investigation, data curation, conceptualization, validation and resources. H.H.: formal analysis, visualization and supervision. R.Z.: funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant no. 42177466).

Institutional Review Board Statement

Not applicable, because this article does not contain any studies with human or animal subjects.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Ju, H.; Niu, C.; Zhang, S.; Jiang, W.; Zhang, Z.; Zhang, X.; Yang, Z.; Cui, Y. Spatiotemporal patterns and modifiable areal unit problems of the landscape ecological risk in coastal areas: A case study of the Shandong Peninsula, China. J. Clean. Prod. 2021, 310, 127522. [Google Scholar] [CrossRef]
  2. Kapustka, L.A.; Galbraith, H.; Luxon, M.; Yocum, J. Using landscape ecology to focus ecological risk assessment and guide risk management decision-making. Toxicol. Ind. Health 2001, 17, 236–246. [Google Scholar] [CrossRef] [PubMed]
  3. Ran, P.; Hu, S.; Frazier, A.E.; Qu, S.; Yu, D.; Tong, L. Exploring changes in landscape ecological risk in the Yangtze River Economic Belt from a spatiotemporal perspective. Ecol. Indic. 2022, 137, 108744. [Google Scholar] [CrossRef]
  4. Formenti, S.; Peace, A.; Eyles, C.; Lee, R.; Waldron, J.W. Fractures in the Niagara Escarpment in Ontario, Canada: Distribution, connectivity, and geohazard implications. Geol. Mag. 2022, 159, 1936–1951. [Google Scholar] [CrossRef]
  5. Liu, S.; Yin, Y.; Liu, X.; Cheng, F.; Yang, J.; Li, J.; Dong, S.; Zhu, A. Ecosystem Services and landscape change associated with plantation expansion in a tropical rainforest region of Southwest China. Ecol. Model. 2017, 353, 129–138. [Google Scholar] [CrossRef]
  6. Fan, X.; Rong, Y.; Tian, C.; Ou, S.; Li, J.; Shi, H.; Qin, Y.; He, J.; Huang, C. Construction of an Ecological Security Pattern in an Urban–Lake Symbiosis Area: A Case Study of Hefei Metropolitan Area. Remote Sens. 2022, 14, 2498. [Google Scholar] [CrossRef]
  7. Zheng, Q.; Lyu, H.M.; Zhou, A.; Shen, S.L. Risk assessment of geohazards along Cheng-Kun railway using fuzzy AHP incorporated into GIS. Geomatics. Geomat. Nat. Hazards Risk 2021, 12, 1508–1531. [Google Scholar] [CrossRef]
  8. Arnous, M.O.; Green, D.R. GIS and remote sensing as tools for conducting geo-hazards risk assessment along Gulf of Aqaba coastal zone, Egypt. J. Coast. Conserv. 2011, 15, 457–475. [Google Scholar] [CrossRef]
  9. Wen, H.J.; Xie, P.; Xiao, P.; Hu, D.P. Rapid susceptibility mapping of earthquake-triggered slope geohazards in Lushan County by combining remote sensing with the AHP model developed for the Wenchuan earthquake. Bull. Eng. Geol. Environ. 2017, 76, 909–921. [Google Scholar] [CrossRef]
  10. Leuven, R.S.; Poudevigne, I. Riverine landscape dynamics and ecological risk assessment. Freshw. Biol. 2002, 47, 845–865. [Google Scholar] [CrossRef]
  11. Ratter, B.M. Surprise and uncertainty—Framing Regional Geohazards in the theory of complexity. Humanities 2013, 2, 1–19. [Google Scholar] [CrossRef]
  12. Gao, M.; Hu, Y.; Bai, Y. Construction of ecological security pattern in national land space from the perspective of the community of life in mountain, water, forest, field, lake and grass: A case study in Guangxi Hechi, China. Ecol. Indic. 2022, 139, 108867. [Google Scholar] [CrossRef]
  13. Lai, Y.; Huang, G.; Chen, S.; Lin, S.; Lin, W.; Lyu, J. Land Use Dynamics and Optimization from 2000 to 2020 in East Guangdong Province, China. Sustainability 2021, 13, 3473. [Google Scholar] [CrossRef]
  14. Huang, S.; Li, X.; Wang, Y. A New model of geo-environmental impact assessment of mining: A multiple-criteria assessment method integrating Fuzzy-AHP with fuzzy synthetic ranking. Environ. Earth Sci. 2012, 66, 275–284. [Google Scholar] [CrossRef]
  15. He, F.; Liu, C.X.; Liu, H.J. Integration and Fusion of Geologic Hazard Data under Deep Learning and Big Data Analysis Technology. Complexity 2021, 2021, 2871770. [Google Scholar] [CrossRef]
  16. Fang, J.P.; Quan, H.C. Risk Analysis of Debris Flow hazardsin Yanbian Prefecture Based on Information Quantity Method. IOP Conf. Ser. Earth Environ. Sci. 2019, 300, 022051. [Google Scholar] [CrossRef]
  17. Van Den Eeckhaut, M.; Hervás, J. State of the art of national landslide databases in Europe and their potential for assessing landslide susceptibility, hazard and risk. Geomorphology 2012, 139, 545–558. [Google Scholar] [CrossRef]
  18. Zhao, J.; Zhang, Q.; Wang, D.; Wu, W.; Yuan, R. Machine learning-based evaluation of susceptibility to geological hazards in the Hengduan mountains region, China. Int. J. Disaster Risk Sci. 2022, 13, 305–316. [Google Scholar] [CrossRef]
  19. Zou, Q.; Jiang, H.; Cui, P.; Zhou, B.; Jiang, Y.; Qin, M.; Liu, Y.; Li, C. A new approach to assess landslide susceptibility based on slope failure mechanisms. Catena 2021, 204, 105388. [Google Scholar] [CrossRef]
  20. Achour, Y.; Saidani, Z.; Touati, R.; Pham, Q.B.; Pal, S.C.; Mustafa, F.; Balik Sanli, F. Assessing landslide susceptibility using a machine learning-based approach to achieving land degradation neutrality. Environ. Earth Sci. 2021, 80, 575. [Google Scholar] [CrossRef]
  21. Ma, Z.J.; Mei, G. Deep learning for geological hazards analysis: Data, models, applications, and opportunities. Earth-Sci. Rev. 2021, 223, 103858. [Google Scholar] [CrossRef]
  22. Dikshit, A.; Pradhan, B.; Alamri, A.M. Pathways and challenges of the application of artificial intelligence to geohazards modelling. Gondwana Res. 2021, 100, 290–301. [Google Scholar] [CrossRef]
  23. Pavlova, I.; Makarigakis, A.; Depret, T.; Jomelli, V. Global overview of the geological hazard exposure and disaster risk awareness at world heritage sites. J. Cult. Herit. 2017, 28, 151–157. [Google Scholar] [CrossRef]
  24. Lin, J.H.; Chen, W.H.; Qi, X.H.; Hou, H.R. Risk assessment and its influencing factors analysis of geological hazards in typical mountain environment. J. Clean. Prod. 2021, 309, 127077. [Google Scholar] [CrossRef]
  25. Quanyuan, W.; Jiewu, P.; Shanzhong, Q.; Yiping, L.; Congcong, H.; Tingxiang, L.; Limei, H. Impacts of coal mining subsidence on the surface landscape in Longkou city, Shandong Province of China. Environ. Earth Sci. 2009, 59, 783–791. [Google Scholar] [CrossRef]
  26. Duan, Y.; Di, B.; Ustin, S.L.; Xu, C.; Xie, Q.; Wu, S.; Li, J.; Zhang, R. Changes in ecosystem services in a montane landscape impacted by major earthquakes: A case study in Wenchuan earthquake-affected area, China. Ecol. Indic. 2021, 126, 107683. [Google Scholar] [CrossRef]
  27. Wu, J.; Zhu, Q.; Qiao, N.; Wang, Z.; Sha, W.; Luo, K.; Wang, H.; Feng, Z. Ecological risk assessment of coal mine area based on “source-sink” landscape theory–A case study of Pingshuo mining area. J. Clean. Prod. 2021, 295, 126371. [Google Scholar] [CrossRef]
  28. Zhang, M.; Bao, Y.; Xu, J.; Han, A.; Liu, X.; Zhang, J.; Tong, Z. Ecological security evaluation and ecological regulation approach of East-Liao River basin based on ecological function area. Ecol. Indic. 2021, 132, 108255. [Google Scholar] [CrossRef]
  29. Hong, W.; Jiang, R.; Yang, C.; Zhang, F.; Su, M.; Liao, Q. Establishing an ecological vulnerability assessment indicator system for spatial recognition and management of ecologically vulnerable areas in highly urbanized regions: A case study of Shenzhen, China. Ecol. Indic. 2016, 69, 540–547. [Google Scholar] [CrossRef]
  30. Li, S.; Zhao, X.; Pu, J.; Miao, P.; Wang, Q.; Tan, K. Optimize and control territorial spatial functional areas to improve the ecological stability and total environment in karst areas of Southwest China. Land Use Policy 2021, 100, 104940. [Google Scholar] [CrossRef]
  31. Chen, D.; Pan, Y.; Jin, X.; Du, H.; Li, M.; Jiang, P. The delineation of ecological redline area for catchment sustainable management from the perspective of ecosystem services and social needs: A case study of the Xiangjiang watershed, China. Ecol. Indic. 2021, 121, 107130. [Google Scholar] [CrossRef]
  32. Tian, H.; Yang, T.; Liu, Q. Climate change and glacier area shrinkage in the Qilian mountains, China, from 1956 to 2010. Ann. Glaciol. 2014, 55, 187–197. [Google Scholar] [CrossRef]
  33. Han, Z.; Song, W.; Deng, X. Responses of Ecosystem Service to Land Use Change in Qinghai Province. Energies 2016, 9, 303. [Google Scholar] [CrossRef]
  34. Liu, J.; Xu, X.; Shao, Q. Grassland degradation in the “Three-River Headwaters” region, Qinghai Province. J. Geogr. Sci. 2008, 18, 259–273. [Google Scholar] [CrossRef]
  35. Bai, A.J.; Zhu, K.X.; Guan, Q. Characteristics of Meteorological Geology Disasters and Rainfall in Qinghai Plateau. Plateau Mt. Meteorol. Res. 2021, 41, 119–124. [Google Scholar] [CrossRef]
  36. Zongxing, L.; Qi, F.; Zongjie, L.; Xufeng, W.; Juan, G.; Baijuan, Z.; Yuchen, L.; Xiaohong, D.; Jian, X.; Wende, G.; et al. Reversing conflict between humans and the environment-The experience in the Qilian Mountains. Renew. Sustain. Energy Rev. 2021, 148, 111333. [Google Scholar] [CrossRef]
  37. Chuanyan, Z.; Zhongren, N.; Guodong, C. Methods for modelling of temporal and spatial distribution of air temperature at landscape scale in the southern Qilian mountains, China. Ecol. Model. 2005, 189, 209–220. [Google Scholar] [CrossRef]
  38. Sun, F.; Lyu, Y.; Fu, B.; Hu, J. Hydrological services by mountain ecosystems in Qilian Mountain of China: A review. Chin. Geogr. Sci. 2016, 26, 174–187. [Google Scholar] [CrossRef]
  39. Deng, S.F.; Yang, T.B.; Zeng, B.; Zhu, X.F.; Xu, H.J. Vegetation cover variation in the Qilian Mountains and its response to climate change in 2000–2011. J. Mt. Sci. 2013, 10, 1050–1062. [Google Scholar] [CrossRef]
  40. Liu, Q.; Yang, Z.; Shi, H.; Wang, Z. Ecological risk assessment of geohazards in Natural World Heritage Sites: An empirical analysis of Bogda, Tianshan. Open Geosci. 2019, 11, 327–340. [Google Scholar] [CrossRef]
  41. Huang, W.; Ding, M.; Li, Z.; Zhuang, J.; Yang, J.; Li, X.; Meng, L.; Zhang, H.; Dong, Y. An Efficient User-Friendly Integration Tool for Landslide Susceptibility Mapping Based on Support Vector Machines: SVM-LSM Toolbox. Remote Sens. 2022, 14, 3408. [Google Scholar] [CrossRef]
  42. Cao, J.; Zhang, Z.; Du, J.; Zhang, L.; Song, Y.; Sun, G. Multi-geohazards susceptibility mapping based on machine learning—A case study in Jiuzhaigou, China. Nat. Hazards 2020, 102, 851–871. [Google Scholar] [CrossRef]
  43. Huang, Y.; Zhao, L. Review on landslide susceptibility mapping using support vector machines. Catena 2018, 165, 520–529. [Google Scholar] [CrossRef]
  44. Xie, G.D.; Zhang, C.X.; Zhang, L.M.; Chen, W.H.; Li, S.M. Improvement of the evaluation method for ecosystem service value based on per unit area. J. Nat. Resour 2015, 30, 1243. [Google Scholar]
  45. Hu, C.; Wang, Z.; Huang, G.; Ding, Y. Construction, Evaluation, and Optimization of a Regional Ecological Security Pattern Based on MSPA–Circuit Theory Approach. Int. J. Environ. Res. Public Health 2022, 19, 16184. [Google Scholar] [CrossRef] [PubMed]
  46. Wang, S.; Wu, M.; Hu, M.; Xia, B. Integrating ecosystem services and landscape connectivity into the optimization of ecological security pattern: A case study of the Pearl River Delta, China. Environ. Sci. Pollut. Res. 2022, 29, 76051–76065. [Google Scholar] [CrossRef] [PubMed]
  47. Grêt-Regamey, A.; Weibel, B. Global assessment of mountain ecosystem services using earth observation data. Ecosyst. Serv. 2020, 46, 101213. [Google Scholar] [CrossRef]
  48. Brardinoni, F.; Grischott, R.; Kober, F.; Morelli, C.; Christl, M. Evaluating debris-flow and anthropogenic disturbance on 10Be concentration in mountain drainage basins: Implications for functional connectivity and denudation rates across time scales. Earth Surf. Process. Landf. 2020, 45, 3955–3974. [Google Scholar] [CrossRef]
  49. Abdullah, S.A.; Hezri, A.A. From forest landscape to agricultural landscape in the developing tropical country of Malaysia: Pattern, process, and their significance on policy. Environ. Manag. 2008, 42, 907–917. [Google Scholar] [CrossRef]
  50. Singh, K.; Byun, C.; Bux, F. Ecological restoration of degraded ecosystems in India: Science and practices. Ecol. Eng. 2022, 182, 106708. [Google Scholar] [CrossRef]
  51. Mir, M.; Maleki, S.; Rahdari, V. Ecosystem restoration and degradation monitoring using ecological indices. Int. J. Environ. Sci. Technol 2022, 20, 1713–1724. [Google Scholar] [CrossRef]
Figure 1. Overview of the study area: (a) location; (b) distribution of landslides and rockfall; (c) land use in 2020.
Figure 1. Overview of the study area: (a) location; (b) distribution of landslides and rockfall; (c) land use in 2020.
Sustainability 16 02976 g001
Figure 2. The technical process of incorporating the results of geological disaster ERA into policies for EFAs.
Figure 2. The technical process of incorporating the results of geological disaster ERA into policies for EFAs.
Sustainability 16 02976 g002
Figure 3. SVM parameter gamma and c selection.
Figure 3. SVM parameter gamma and c selection.
Sustainability 16 02976 g003
Figure 4. Accuracy verification.
Figure 4. Accuracy verification.
Sustainability 16 02976 g004
Figure 5. Results of ERA (see Table 1 for explanation of variables).
Figure 5. Results of ERA (see Table 1 for explanation of variables).
Sustainability 16 02976 g005
Figure 6. (ad) Spatial hot and cold spot analysis of geological disaster ecological risk, hazard, vulnerability, and exposure. (e) Combination pattern of geological disaster ecological risk.
Figure 6. (ad) Spatial hot and cold spot analysis of geological disaster ecological risk, hazard, vulnerability, and exposure. (e) Combination pattern of geological disaster ecological risk.
Sustainability 16 02976 g006
Figure 7. Identification of ecological risk sources. (a) the correlation between hazard, vulnerability, and exposure; (b) the average value of hazard, vulnerability, and exposure of different ecological risk levels. (c) SEM of Ecological Risk. (dh) the contribution of different environmental factors to ecological risk.
Figure 7. Identification of ecological risk sources. (a) the correlation between hazard, vulnerability, and exposure; (b) the average value of hazard, vulnerability, and exposure of different ecological risk levels. (c) SEM of Ecological Risk. (dh) the contribution of different environmental factors to ecological risk.
Sustainability 16 02976 g007
Figure 8. Spatial strategies for addressing ecological risks of geological disasters. (a) systematic clustering of ecological functional units; (b) spatial characteristics of ecological risk; (c) space policy map.
Figure 8. Spatial strategies for addressing ecological risks of geological disasters. (a) systematic clustering of ecological functional units; (b) spatial characteristics of ecological risk; (c) space policy map.
Sustainability 16 02976 g008
Figure 9. Framework for establishing the connection between ERA and EFAs.
Figure 9. Framework for establishing the connection between ERA and EFAs.
Sustainability 16 02976 g009
Table 1. Indicators used to assess ecosystem vulnerability.
Table 1. Indicators used to assess ecosystem vulnerability.
IndicatorCalculationLandscape Ecology Significance
Interference degree of landscape (U)U = a × FN + b × FI + c × FDThe degree of landscape disturbance indicates the degree of loss in the area after being disturbed. Based on expert opinions, take values of 0.5, 0.3, and 0.2 for a, b, and c, respectively.
Fragmentation (FN) F N = M P S × ( N f 1 ) / N c The degree to which the landscape is fragmented and broken reflects the extent of human interference with the landscape to a certain extent.
Isolation of landscape (FI) F I = 1 2 n A × A i A Used to reveal the degree of separation of individual patches in landscape types. The higher the detachment of patches, the smaller their ability to resist risks, and the lower their landscape safety.
Fractal dimension (FD) F D = 2 ln ( P / 4 ) / ln A The smaller the index value, the less likely it is to be interfered with by human activities, while the larger the index value, the more likely it is to be interfered with by human activities.
Landscape vulnerability (V)Experts knowledge acquiredBarren, Sonw/Ice, Water, Wetland, Grassland, Forest, Shrub, Cropland, and Impervious are divided into 9, 8, 7, 6, 5, 4, 3, 2, and 1, respectively.
Landscape adaptation index (LAI)LAI = PRD + SHDI + SHEILAI is determined by referencing other people’s research and is composed of the patch richness index, Shannon diversity index, and Shannon evenness index.
Patch richness index (PRD) P R D = m / A PRD is the number of patches per unit area; it reflects the dispersion of patches in a landscape type.
Shannon diversity index (SHDI) S H D I = i = 1 m ( p i × ln p i ) Reflects the richness and complexity of the landscape.
Shannon evenness index (SHEI) S H E I = i = 1 m ( p i   ×   ln p i ) ln m Indicates the maximum possible diversity of a given landscape’s richness. Reflects the stability of the ecosystem.
Table 2. Statistics on the ecological risk of geological disasters and the combination patterns.
Table 2. Statistics on the ecological risk of geological disasters and the combination patterns.
Combination PatternLow Risk Medium RiskHigh Risk
Functional Unit/NumberShare/%Functional Unit/NumberShare/%Functional Unit/NumberShare/%
Low hazard–low vulnerability–low exposure8500000
Low hazard–low vulnerability–high exposure0017.700
Low hazard–high vulnerability–low exposure425430.8116.7
Low hazard–high vulnerability–high exposure16.30000
High hazard–low vulnerability–low exposure318.8215.400
High hazard–high vulnerability–low exposure00646.2583.3
Total16100131006100
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Long, X.; Xiang, Q.; Zhang, R.; Huang, H. Incorporating the Results of Geological Disaster Ecological Risk Assessment into Spatial Policies for Ecological Functional Areas: Practice in the Qilian Mountains of China. Sustainability 2024, 16, 2976. https://doi.org/10.3390/su16072976

AMA Style

Long X, Xiang Q, Zhang R, Huang H. Incorporating the Results of Geological Disaster Ecological Risk Assessment into Spatial Policies for Ecological Functional Areas: Practice in the Qilian Mountains of China. Sustainability. 2024; 16(7):2976. https://doi.org/10.3390/su16072976

Chicago/Turabian Style

Long, Xu, Qing Xiang, Rongguang Zhang, and Hong Huang. 2024. "Incorporating the Results of Geological Disaster Ecological Risk Assessment into Spatial Policies for Ecological Functional Areas: Practice in the Qilian Mountains of China" Sustainability 16, no. 7: 2976. https://doi.org/10.3390/su16072976

APA Style

Long, X., Xiang, Q., Zhang, R., & Huang, H. (2024). Incorporating the Results of Geological Disaster Ecological Risk Assessment into Spatial Policies for Ecological Functional Areas: Practice in the Qilian Mountains of China. Sustainability, 16(7), 2976. https://doi.org/10.3390/su16072976

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop