1. Introduction
In the context of global change, uncertain factors in the environment have gradually increased, coupled with the increasing degree of interference of human activities on the ecosystem, leading to imbalances in the internal structure of the ecosystem, damage to self-repair capabilities, and deterioration of the quality of the ecological environment [
1,
2,
3,
4,
5,
6]. It is urgent to assess the disturbance of stress factors such as environment, society, and economy on the structure and function of the ecosystem, aim to identify ecological risk (ER) sources, regulate ER processes, and then improve ecosystem service efficiency [
7,
8,
9,
10]. In recent years, ER has become a hot topic of ecological and geographic research [
11,
12,
13]. The ER concept is that an ecosystem, landscape, or region is subject to external threats such as natural- or human-induced, increasing the possibility of reduced productivity, environmental damage, and decreased ecosystem service function. At present, the “source-sink” method and landscape index method are used to evaluate ER [
14,
15]. The “source-sink” method focuses on the identification of “source-sink” properties of landscape risks and is widely used in research of non-point source pollution, soil erosion, and urban heat island effect [
16]. The “source-sink” method applies to the evaluation of ERs with clear threats at the regional level [
17]. Based on land use/cover data, the landscape index method focuses on the ecological environment effects caused by specific landscape assemblages and their spatial heterogeneity. The landscape index emphasizes the interconnectedness of patterns and processes [
18,
19]. Landscape disturbance index and landscape vulnerability index are the main indicators for constructing ER in landscape ecology; while the landscape disturbance indices consist of a landscape fragmentation index, a landscape separation index, and a landscape subdimensionality index [
20,
21]. Compared with the “source-sink” risk assessment, the landscape index method assesses ER at the landscape level. This method can identify the local spatial characteristics of ER, scientifically understand the interrelationship between landscape patterns and ecological processes and diagnose multiple ER sources and their cumulative effects in the landscape mosaic, which provides a theoretical reference for promoting the healthy and sustainable development of ecosystems.
Scale characteristics have been one of the main focuses of ecological and geographic research, and ecosystem structure and function tend to change with scale. Assessing the extent of the disturbance of ecosystems from external threats at different scale units is a major tool for developing ER management measures in the region. Currently, studies on the scale of ERs have been deepened, and ER assessment on different geographical units such as watersheds [
22,
23], mountainous areas [
24], cities [
25,
26] and beaches [
27] have been carried out. In recent years, an increasing number of scholars are beginning to assess ER at a fine scale. Zhang et al. [
28] analyzed the ER of the high background area of heavy metals in the southwest from the township scale, making up for the shortcomings in precision of previous surveys at the regional scale, and meeting the needs of natural resource management. Zhang et al. [
29] evaluated ER in southwest Guangxi-Beibu Bay based on the DPSIR model and OWA-GIS. Previous studies have shown that there are significant differences in the assessment of ER at different scales. If the assessment scale is too large, it is difficult to identify the local characteristics of ER. In contrast, if the assessment scale is too small, the overall characteristics of the region are often ignored. Therefore, the selection of the optimal assessment unit is a prerequisite for the accurate identification of the spatial differentiation characteristics of ERs and a basis for quantifying the impact of natural factors and the regional ecological environment.
The Funiu Mountains are the transition region from the north subtropical to the south warm temperate zone in China, and the area where the north–south divide passes through. It is also an important water-conserving forest area of the Danjiangkou reservoir, with climate regulation, water supply, and ecosystem service function. The ecological quality of the north–south transition region in China plays an important role in conserving biodiversity and enhancing ecosystem service capacity. Since the 21st century, the land-use pattern was profoundly changed by rapid urbanization, resulting in great pressure on the stability of the internal structure in the ecosystem. Given this, we constructed the ecological risk index by landscape disturbance degree and landscape vulnerability degree. The exploratory data analysis and spatial autocorrelation were used to diagnose the difference of ERs at different scales. Furthermore, correlation analysis and geographically weighted regression (GWR) were used to analyze the strength and spatial difference of the impact of various environmental factors on ER in 2015, to provide scientific reference for the ecological environment planning in the study area.
This article aimed to answer three questions: (i) What is the best scale for assessing ER in the Funiu mountain area—watershed, township, or grid? (ii) What changes in the ER and development trends of different grades of ER occurred from 2000 to 2015? (iii) At the scale of a 1 km grid, what are the intensity, and direction of influencing factors of ecological risk? The results of this study can provide relevant references for ecological quality assessment of other geographic transition areas in China and the world.
2. Data and Methodology
2.1. Study Area
Located in the transition zone from the north subtropical to the warm temperate zone in China, the Funiu Mountains are a typical ecological and geographically fragile area with rich biodiversity and important ecosystem services in the region. The vegetation type is the transition type of warm temperate deciduous broadleaf forest and north subtropical evergreen deciduous mixed forest. The temperature range is 13.6–15.1 °C, and the precipitation range is 700–1000 mm. In addition, this area is in the west of Henan Province with a total area of 19,711.75 km
2. The study area ranged from 110°30′ E to 113°30′ E and from 32°45′ N to 34°20′ N, reaching the border of Henan and Shaanxi in the west, north of Fangcheng city in the east, Xionger Mountain and Waifang Mountain in the north, and Nanyang Basin in the south (
Figure 1).
2.2. Data Resources
The land cover/use data (2000–2015) were obtained from the Middle and Lower Yellow River Scientific Data Center, China Earth System Data Sharing Platform (
http://www.geodata.cn/) (accessed on 2 May 2019). According to the “China land use/cover remote sensing monitoring data classification system”, the land-use types in each period were reclassified into six types, including construction land, forest, watershed, grassland, cultivated land, and unused land. The digital elevation model (DEM) data came from the Geospatial Data Cloud Platform (
http://www.gsclod.cn) (accessed on 8 May 2019), with a resolution of 30 m. The slope was obtained by extracting DEM data. Climate data, population density, and watered zoning data were collected from the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (
http://www.resdc.cn/) (accessed on 12 May 2019). The land-use degree and human activity intensity were calculated according to the method of Zhuang et al. [
30] and Yan et al. [
31].
2.3. Methods
In this paper, we constructed the land use risk index using the landscape interference index and landscape vulnerability index in the landscape ecology. The exploratory data analysis was used to evaluate the optimal unit of ER. The GIS software was to analyze the spatio-temporal characteristics. The correlation analysis and geographically weighted regression models (GWR) were applied to estimate the influence degree of different environments on the ER. The overall framework of this paper is shown in
Figure 2.
2.3.1. Construction of Land-Use Risk Index Model
(1) Landscape interference index
The landscape disturbance index (E
i) refers to the degree of external disturbance to different landscape ecosystems. The lower the intensity of landscape disturbance, the more stable the ecosystem, the more beneficial it is to the survival of organisms. In this paper, a landscape fragmentation index, a landscape separation index and a landscape subdimensionality index were used to construct the landscape disturbance index and were measured through the method of Wu et al. [
32] and Xie et al. [
33]:
where C
i was the landscape fragmentation index, which can characterize the degree of ecosystem landscape heterogeneity. S
i was the landscape subdimensionality index used to describe the complexity of the landscape patches. DO
i was the landscape separation index; the larger the value, the more complex the landscape distribution and the greater the degree of external interference.
a,
b, and
c represented the weights of each landscape index, and the sum of the weights was 1. Combined with previous studies, the values of
a,
b, and
c were assigned as 0.502, 0.301, and 0.197, respectively [
34,
35].
(2) Landscape vulnerability index
Landscape vulnerability (
Fi) was the degree of susceptibility in different ecosystems under the stage of the natural succession process of the landscape, which can characterize the integrated effect of human activities and the natural environment on land use. When the ecosystem is in the primary succession stage, the food chain structure is simple and has less ability to resist the risk of external disturbance, while the more complex the ecosystem structure, the greater the ability to resist the risk of external disturbance. According to the research of Zhao et al. [
36] and Xie et al. [
33], unused land was the most vulnerable, followed by watershed, while construction land was the most stable. Therefore, the unused land, grassland, watershed, cultivated land, forest land, and constructed land were assigned values of 6, 5, 4, 3, 2, and 1 in descending sequence, and the
Fi values of each type were normalized to be 0.2857, 0.2381, 0.1905, 0.1429, 0.0952, and 0.0476, respectively [
37].
(3) Land use ecological risk index
The land use ecological risk index (
ERIi) was constructed from the landscape disturbance index and landscape vulnerability index, and grid units were created to calculate the ecological risk value [
33].
ERIi was defined as follows:
where
ERIk was the ecological risk index value of land use on the kth grid.
Ak was the sampling area of the kth grid.
Aki was the area of land-use type
i on the kth grid.
Ei was the disturbance index of land use type I, and
Fi was the corresponding vulnerability index.
2.3.2. Analysis of Ecological Risk Scale Characteristics
(1) Different scales of ecological risk assessment
To evaluate the spatial clustering characteristics of ER at different scales (
Figure 3), the grid-scale (1 km × 1 km, 1.5 km × 1.5 km, 3 km × 3 km, 5 km × 5 km), township administrative unit scale, county administrative unit scale, and watershed-scale were constructed through resampling and magnitude conversion in Geographic Information Systems (GIS). By analyzing the change rule of ER under scale transformation, the best unit for evaluating ER can be determined.
(2) Exploratory data analysis
The exploratory data analysis was based on spatial statistics. The main object of the research was the spatial interaction and changes to the rules between events or phenomena and geospatial characteristic information, to analyze the randomness and structural characteristics, spatial correlation, dependence, and heterogeneity of variables in a specific region [
38]. Spatial autocorrelation is one of the commonly used tools for spatial statistical analysis, which accurately represents the attribute association characteristics of a geographic element on a specific area and adjacent units [
39]. In general, both global autocorrelation and local autocorrelation indicate the degree of spatial autocorrelation. The Global Moran’s I statistic and Getis–Ord Gi* index were used to measure the spatial autocorrelation of ERs in different scales with the software of Geoda and GIS. The Global Moran’s I was defined as follows:
where n was the number of space units, indexed by i and j. Y
i and y
j were the attribute values of unit i and j, respectively. W
ij was the weight matrix. S
2 represented the sample variance. Y
mean was the mean of y. The Global Moran’s I was the spatial autocorrelation value of ER, the value of I ranged from −1 to 1. The Global Moran’s I value greater than zero indicated a spatial positive correlation. The larger the index, the stronger the positive spatial correlation, and an index close to 1 showed a smaller degree of spatial variation. A Moran’s I value less than zero indicated a spatial negative correlation. The lower the index, the stronger the negative spatial correlation, and an index close to −1.0 meant very large spatial differences. A Moran’s I value trending toward zero showed the random spatial distribution of ER. The Getis–Ord Gi* index was calculated as follows [
38]:
where
E(
Gi*) denoted the mathematical expectation of
Gi*.
Var(
Gi*) was the variance of Gi*, and W
ij was the spatial weight. The values of
Z (
Gi*) were classified into four classes: <−1.96, −1.96 to −1.65, −1.65 to 1.65, 1.65 to 1.95, and >1.96, representing the cold spot area, subcold spot area, non-significant change area, subhot spot area, and hot spot area, respectively.
2.4. Correlation Analysis
Correlation analysis was used to measure the correlation between ER and environmental variables and quantify the degree of correlation between variables. If the value was greater than zero, the correlation was positive. Otherwise, it was a negative correlation. The value represents the degree of correlation, and the higher the value, the stronger the correlation [
39]. The definition of correlation analysis was shown as follows:
where
Rxy was the correlation coefficient between two variables.
N was the number of samples.
Xi was the
i value of
x, and
x mean was the average value of
x. Also,
yi and
ymean represented
y and the average value of
y, respectively.
2.5. Geographically Weighted Regression Models
GWR models were improved based on the classical linear regression model (OLS). The spatial scale dependence and heterogeneity of the variables were reduced, and the local regional parameters and their influencing factor were estimated. The non-smoothness of correlation between the variables and the spatial location was found, and the results were more accurate [
40,
41,
42].
where y
i was the dependent variable, and x
ik was the independent variable. (
μi,
vi) was the latitude and longitude of the sampling point i.
βk (
μi,
νi) represented the value of the continuous function
βk (
μ0,
ν0) on the space element of the k sample.
was the regression residual of point i. i = 1, 2, ...,
n, was the number of positions in space.
3. Results
3.1. Scale Characteristics of Ecological Risk in Funiu Mountain Area
3.1.1. Cluster Analysis of Ecological Risk at Different Scales
As shown in
Figure 4, with the evaluation scale gradually pushed down (from 5 km to 1 km), the spatial distribution pattern of ERs became more precise, and the distribution of hot spots and subhot spots become more and more, and the significance level became higher. When the evaluation unit increased from the 5 km grid scale to a watershed scale, the spatial correlation of the ER distribution became weak, and the distribution area of hot spots and subhot spots was insignificant. The results showed that the spatial agglomeration characteristics of ERs were different at various scales, and the larger the evaluation scale, the more obvious the impact of landscape composition and different types of land ecological subsystems on the results of ecological risk assessment. As a result, differences in the composition and structure of the interior landscape were filtered. Therefore, compared with administrative units and watersheds, grid units were more suitable for evaluating the spatial–temporal characteristics and influence mechanisms of ER.
3.1.2. Spatial Autocorrelation Analysis of Ecological Risk at Different Scales
It was clear that the spatial correlation degree of ER varied significantly at different scales (
Table 1). The Global Moran’s I value of ER decreased from 0.53 to 0.04 when the assessment scale moved from the 1 km grid unit to the watershed scale. The correlation of ER decreased with the increase of scale, and the larger the assessment scale, the lower the spatial autocorrelation. The significance test showed that the Z-scores of 1 km, 1.5 km, 3 km, 5 km, and township scale were all greater than 1.96, indicating that both grid-scale and township scale passed the significance test, while the Z-scores of county scale and watershed scale were 1.50 and 1.30, respectively, failing the significance test (
p > 1.96). The results of this study showed that the spatial characteristics and spatial correlation of the 1 km grid were high. Therefore, the spatial and temporal characteristics of the ER and its influencing factors in the study area were analyzed based on the 1 km grid.
3.2. Spatio-Temporal Differentiation Characteristics of Ecological Risk
As can be seen from
Figure 5, the spatial distribution of ER in the study area did not change significantly from 2000 to 2015. The mild ecological risk zone was the most extensive, mainly concentrated in the southwest and central regions. The moderate ecological risk zone and high ecological risk zone were clustered and distributed in the low altitude area, mainly due to the influence of human activities in the low altitude areas, while the distribution of the mild ecological risk zone was small and scattered.
The ER levels in different periods were counted, and the area and proportion were calculated (
Table 2). In this study, the mild ecological risk zone was the most widely distributed and had the highest proportion in 2015 (69.77%) but the lowest proportion in 2010 (65.74%). The moderate ecological risk zone was the largest (2906.53 km
2) in 2015, accounting for 14.75%, and the smallest (2351.22 km
2) in 2005, accounting for 11.93%, reflecting the phase change characteristics of moderate ecological risk decreasing first and then increasing. The distribution area of the slight ecological risk zone was the smallest (2151.08 km
2) in 2000 and the largest (2641.54 km
2) in 2005, which showed the phase characteristics of first increasing and then decreasing. The highest level of ecological risk and ecological stress were found in the high ecological risk zone. In 2010, the high ecological risk zone was the most widely distributed (1896.16 km
2), accounting for 9.62%; in 2015, it was the smallest (755.84 km
2), accounting for 3.83%, and the overall high ecological risk area showed a downward trend. During the study period, the area of different levels of ER changed significantly, where the area of high ecological risk zone decreased by 398.07 km
2, and the proportion dropped by 2.02%, indicating that the high-risk level was controlled and the degree of risk stress on the ecosystem was reduced, which also reflected the trend of transferring from high ecological risk level to other ecological risk levels. The area of mild ecological risk increased by 204.74 km
2, and the area of slight ecological risk and medium ecological risk increased by 144.48 km
2 and 49.33 km
2, respectively.
3.3. Analysis of Ecological Risk Influencing Factors
3.3.1. Correlation Analysis
ArcGIS 10.3 was used to extract the value of ER in the 1 km grid in 2015, and the function of multivalue extraction to a point feature was used to extract the impact of human activities intensity, land-use degree, population density, maximum temperature, slope, and altitude. SPSS software was used for correlation analysis to calculate the pair correlation between the impact factors. The correlation among the influencing factors passed the significance test, except for precipitation and ER and maximum temperature and land-use degree, which did not pass the highly significant test (
p < 0.01) (
Table 3). In addition, the correlation coefficients among the factors were both positive and negative, indicating that the directing of action varied greatly among the influencing factors. Slope and altitude were negatively correlated with other influencing factors, and had the strongest influence on maximum temperature, with coefficients of −0.483 and −0.604, respectively. The correlation coefficients of land-use degree with ER and human activity intensity were 0.768 and 0.696, respectively, suggesting that the change in land-use degree significantly affected the ER and human activity intensity. The ER and human activity intensity were stronger for land-use types with higher development degrees. Both altitude and slope were negatively correlated with ER at the 0.01 level, the correlation coefficients were −0.116 and −0.084, respectively.
3.3.2. Model Selection
To select a suitable model for analyzing the influence of each environmental variable on ER, six environmental variables were selected as independent variables, with ER as the dependent variable in the OLS model to measure the influence of each explanatory variable on the explanatory variable and its significance level (
Table 4). The variance inflation factor (VIF) of each parameter of the regression model was less than 7.5, demonstrating that the equation variables were set reasonably and there was no multicollinearity problem. However, the OLS model only considered the global characteristics of the regression coefficients and did not consider the local effects and spatial spillover of each influencing factor on ER, and there were spatial clustering characteristics of ER. In addition, the whole model had a strong spatial autocorrelation based on the results of Global Moran’s I, so there would be a large error existing in the OLS model which cannot explain the spatial heterogeneity of the risk factors well. Hence, the GWR model was further used to identify the spatial heterogeneity of ER influencing factors, and the parameters of the GWR model were obtained with the OLS model regression results as a reference (
Table 5). Compared with the OLS model results, the GWR model fit was 0.9038, which was 40.25% better than the OLS model; the sigma and residual sum of squares were statistically significantly lower, with a decrease in the value of the deficit pool information criterion. Therefore, we concluded that the GWR model was significantly better than the OLS model.
The diagnostic results of the OLS model showed that all six explanatory variables passed the significance test at 0.01 level. The influence of land-use degree, population density, and maximum temperature on soil moisture was positively correlated, and the explanation level was maximum temperature, land-use degree, and population density from high to low. Though it remained unchanged under other conditions, the ER would increase by 0.0574% when the maximum temperature with the strongest explanatory ability increased by one unit. The values of the correlation coefficients of human activity intensity, slope, and DEM on the ER were negatively correlated, and the absolute value of the correlation coefficients in descending order were human activity intensity, slope, and DEM. The ER decreased by 5.85% with each increase of human activity intensity, and the ER increased slightly with the decrease of slope and altitude.
3.3.3. Spatial Differences of Impacts of Explanatory Variables on Ecological Risk
(1) Descriptive statistics of regression coefficients of GWR model
Descriptive analysis was conducted on the regression coefficients of each specific spatial unit of GWR, and the coefficient values of explanatory variables were obtained statistically, including minimum value, lower quartile value, median value, upper quartile value, maximum value, and average value (
Table 6). It showed that each explanatory variable was positive or negative with a large difference, indicating that each explanatory variable was not simply positively or negatively correlated with risk factors. There was obvious spatial heterogeneity in the degree of explanation, and the impact on the ER may be opposite in different spatial scopes. However, the direction of the median value and the mean value was consistent, suggesting that the GWR model had the same effect on ER in most spatial ranges, and the difference of its value indicated that each explanatory variable had different effects.
(2) Spatial differences of the impact of explanatory variables on ecological risk
The natural breakpoint method in ArcGIS 10.3 was used to visualize the results of model operation and reflect the influence characteristics of explanatory variables on the ER at the local spatial scale.
Figure 6 presents the spatial coefficient distribution of six impact factors generated by the GWR model, indicating that there were clear spatial variations in the degree of influence of the ER impact factors: (1) In terms of human activity intensity, the regression coefficient decreased from central to peripheral, indicating that forest and grassland were the main land-use types in the central part, which were greatly affected by human activity intensity and had a weak ability to resist ER. The northwest and southeast regions were occupied land for urban residents, and the terrain was flat, and the ER decreased when other impact factors remained unchanged, and the human activity intensity increased. (2) In terms of land-use degree, there was a positive correlation between land-use degree and ER, displaying low levels in the middle and high surrounding areas. The high-value areas were mainly concentrated in the northwest, north, and south, mainly in areas with high land development, such as cultivated land and construction land. In addition, the interchange between various land-use types led to alterations in the ecological level. The government departments should implement some measures to protect the ecological land and relieve the external pressure on the ecosystem. (3) In terms of population density, the regression coefficient ranged from −0.61 to −1.05, showing a negative correlation. The low area gathered in an ellipse in the northwest, southwest, northeast, and high-value areas distribution in the northeast and southeast, located in the built-up areas of Lu, Nanzhao, and Zhenping county, the high population density, and flat topography was the main reason. If other factors remain unchanged, the population density increasing would aggravate ER. (4) In terms of elevation, the regression coefficients were negative in most regions, and the absolute values of the regression coefficients increased from the central region to the surrounding areas, indicating that the ER in the central region was less negatively affected by the average elevation than in other regions. (5) In terms of slope, the regression coefficients were concentrated in [−2.08, −1.53], exhibiting negative correlation characteristics, and the high absolute values of the regression coefficients were mainly distributed in the northeast, meaning the ER in the northeast was negatively affected by the average slope to a greater extent. (6) In terms of maximum temperature, the regression coefficients generally showed a trend of increasing from the northwest to the southeast, and a positive correlation was a dominant relationship. The ER in the southeast was more likely to be positively affected by raised maximum temperatures and was significantly affected by monsoonal climate, especially extreme climate changes that increased the threats to ecosystems. In addition, the flat topography of the southeast was less influenced by topographic factors than the northwest, which was a potential cause of increased ER from warmer temperatures in the east.
4. Discussion
4.1. Significance to Management
The Funiu Mountains, situated in the transition area between the north and the south of China, are one of the richest areas in terms of biodiversity conservation and the most sensitive ecological environment in China. We calculated the landscape pattern index by constructing a 1 km grid to evaluate the degree of ER. It found that the ER showed a trend of increasing and then decreasing, with the highest ER value in 2010 and the lowest ER value in 2015. The rapid urbanization in mountainous areas has changed the way and extent of land use, which has a significant impact on ER in the region. In addition, the ecological risk index showed a decreasing trend after 2010, and the ecological environment quality improved. The main reason being that the implementation of the policy of returning farmland to forest and the environmental protection concept of “Green water and green mountains are golden mountains” proposed by Xi Jinping has increased the protection of ecological land and curbed the trend of converting ecological land such as forest land to construction land. In the future, the implementation of the policy of returning farmland to the forest as well as ecological conservation should be implemented in the Funiu Mountains.
From 2000 to 2015, the ER showed an overall spatial characteristic of high and low in the middle, which is mainly because the topography of the Vernon Hills is high in the northwest and low in the southeast. The surrounding terrain is relatively flat, the degree of human activity influence is high. The construction land and cultivated land are widely distributed in low altitude areas, which are subject to a greater degree of human disturbance. The high degree of land fragmentation is an important factor contributing to the formation of the spatial distribution pattern of ER. Hence, government departments should pay attention to the protection of ecological land in the formulation of socio-economic development plans, strengthen the protection of national and regional ecological policies, plan ecological environmental protection well, and balance the relationship between socio-economic development and ecological environmental protection.
In the context of global changes, managers should propose policies that promote sustainable land use, focus on ecological land-use change in sensitive areas, implement red lines for farmland protection and ecological function zones, build buffer zones for the conversion of ecological land and construction land, set biodiversity protection zones, and implement the strictest land control system to promote the sustainable use of land resources in the region.
4.2. Scale Characteristics of Ecological Risk in Funiu Mountains
The selection of suitable risk evaluation units is the precondition for the accurate assessment of ER in the region, and how to determine ER evaluation units on different geographical units is a hot topic of academic interest [
43,
44]. In this study, the response characteristics of ER to scale were analyzed at grid, township, county, and watershed levels. The results showed that the hotspot areas and the Global Moran’s I values of ER increased as the scale of the 5 km grid was pushed down to 1 km, and both the distribution of hotspot areas and the Global Moran’s I values decreased as the scale was pushed from the grid to watershed. Therefore, the 1 km grid scale was determined to be the best unit for assessing ER, which was consistent with the research on the evaluation of ecosystem service values in the Dabie mountains [
45]. In addition, the ecosystem is a very complex giant system, while this study only assessed the ER from the grid unit, which was not sufficient for the formulation of an ecological function zoning scheme and policy implementation. In future research, a region-wide assessment of ecoenvironmental quality and a quantitative analysis of capabilities and trends in ecosystem services, coupled with “water-soil-gas-biology-human” and other elements will be used to measure the external stress factors of ecosystems in a multilevel, comprehensive way.
4.3. Spatial and Temporal Evolution of Ecological Risks and Analysis of Influencing Factors
The landscape index method was used to assess the spatial and temporal characteristics of ER in the Funiu Mountains. The results displayed that the ecological environment quality improved from 2000 to 2015, and the area of high ER decreased the most, while the area of slight ER increased the most, which reflected the mutual conversion and development trend of different ecological risk levels in the region. Overall, due to the implementation of ecological and environmental protection policies such as returning farmland to forest, rural residential renovation, and mining land restoration, the ecological and environmental index of farmers showed a rising trend. From the viewpoint of space, the high ecological risk zone and moderate ecological risk zone were mainly located in the southeast and at the edge of the area, mainly because of the clustered distribution of construction and the strong degree of interference by human activities intensity. Hence, it required further strengthening the regulation and prevention of high-level ecological risk areas, controlling the disorderly expansion of construction land and strictly adhering to the ecological red line.
Topographical and geomorphological factors, construction of the settlement, policymaking, construction of protected areas, and the degree of anthropogenic disturbance can all have important effects on the strength of ER [
46]. Ran et al. [
47] discovered that land-use change is an important cause of ER. Wang et al. [
48] constructed an ecological risk evaluation index system for Dongting lake, incorporating population density and gross regional product into the index system, based on the “pressure-response” model. In this paper, the aim was quantitatively exploring the degree of influence of external stress factors on ER and spatial differences, and indicators were selected from both natural factors and human disturbances. The natural factors were chosen as elevation, slope, and maximum temperature. Human activity intensity, land-use degree, and population density were selected as human disturbance factors. The topography was the key element in determining the quality of the mountain’s ecological environment. As a result, the elevation and slope were chosen as the representative factors. Existing research demonstrated that the maximum temperature has a significantly higher intensity of impact on ER than other climatic elements [
45], and multiple covariances were found in this study when testing for covariance among environmental variables. Therefore, the maximum temperature was considered as a typical climatic element in this paper. Human activity intensity, land-use degree, and population density all indicated the degree of impact of anthropogenic disturbances on ecosystem structure and function, thus identifying the direction and degree of impact of anthropogenic activities on ER under different forms. The correlation analysis and GWR model were used to measure the intensity and spatial variation of each environmental variable on the ER, respectively. We found that the topographic and geomorphic factors had a strong influence on ER, which was consistent with the elevation, slope, and lithology being negatively correlated with the ER of landslide hazards in assessing ER in southwestern mountains [
45]. Ni et al. [
49] discovered that the ecological level was lower in areas with a high intensity of human engineering impacts and high population density, based on the assessment of geological risks in the Changbai mountains. The results of this paper indicate that there was spatial heterogeneity in human activities intensity, land-use degree, and population density with ER, which were manifested in differences in intensity and direction. In some areas, active protection measures may be taken by human beings so that the local ER will be reduced. ER is a joint natural-social-human influence on the ecosystem, and the selection of indicators is of great significance for indepth analysis of the formation mechanism of ER. Vegetation Normalized Difference Vegetation Index, extreme climatic environment (solar radiation, droughts, and floods), soil physical and chemical properties, soil erosion, and socio-economic factors (GDP, pesticides, mechanization level, urbanization level) should be considered to further research ER in the future. The application of high-resolution remote sensing data could improve the accurate analysis of temporal and spatial characteristics of ecological risk and provide the theoretical basis and empirical analysis for improving ecological environment quality.
5. Conclusions
This paper used spatial autocorrelation to identify the spatial aggregation characteristics of ecological risks at the grid, township, county, and watershed scales. The results showed that ecological risk assessment is scale-dependent, and there are differences in the spatial differentiation of ecological risk on different scale units. From the 1 km grid scale to the watershed scale, the correlation of the ecological risk spatial distribution grows weaker and weaker, the hot and cold areas gradually disappear, the Global Moran’s I value also drops from 0.53 to 0.04, and the Z-score is reduced from 109.02 to 1.30. The significance level also shows a downward trend, and it is further proved that the 1 km grid is the most suitable unit for evaluating the ecological risk. Therefore, we applied the landscape index model to calculate the ecological risk in the Funiu mountain area on the 1 km grid from 2000 to 2015. The results indicated that the ecological environment quality showed a good trend, mild ecological risk zone is the main ecological risk level, accounting for more than 60%. During the study period, the area of high ecological risk decreased by 398.07 km2, and the proportion decreased by 2.02%, while the area of slight ecological risk increased the most (204.74 km2), the area of mild ecological risk area increased the most (144.48 km2), and the area of moderate ecological risk area increased the least (49.33 km2). The implementation of the policy measures to return farmland to the forest in mountainous areas is the main reason for the improvement of the ecological environment quality. In 2013, the introduction of the policy that green mountains are golden mountains increased the protection of ecological land and improved the ability of the ecosystem to resist external disturbances. Moreover, we took ER in 2015 as the subject of our investigation and analyzed the direction and degree of influence of different factors on ER. In terms of the correlation among various influencing factors, elevation and slope are negatively correlated with other influencing factors. The correlation coefficient between land-use degree and human activity intensity and ER is the highest, being 0.768 and 0.696, respectively. Human activity intensity, land-use degree, population density, altitude, slope, and the highest temperature had a significant influence on the ER, and there is space differentiation. Land-use degree and human activity intensity were the main obstacles ER to improvement; while population density, altitude, slope, and the highest temperature of ER revealed that there is an influence of different direction and strength.
In this study, we analyzed the spatio-temporal characteristics of ER and driving mechanisms based on multisource data, such as land use/cover data from 2000 to 2015, population density, elevation, slope, and maximum temperature, etc. The timeliness of the data may incur certain errors in the assessment results. Hence, the latest and higher resolution data should be collected to calculate ER, including land use/cover data, population density, and climate data in 2020. The accuracy of the data will be improved and provide suggestions for governmental departments to develop ecological environment planning.
In addition, this paper constructed ER through a limited landscape disturbance index and landscape vulnerability index, which may have some error influence on the results, and more indicators need to be selected to comprehensively evaluate ecological risk. Lastly, there is an urgent need to assess the impact of natural environment-socio-economic-human activities on ecological risk and to identify the driving factors of ecological risk.
Author Contributions
Conceptualization, L.L. and X.Z.; methodology, L.L.; software, L.L.; validation, L.L., L.Y. and J.D.; formal analysis, L.L.; investigation, Z.Z.; resources, X.Z.; data curation, L.L.; writing—original draft preparation, L.Y.; writing—review and editing, X.Z.; visualization, J.D.; supervision, Z.Z. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Not applicable.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Suter, G.W.; Norton, S.B.; Barnthouse, L.W. The evolution of frameworks for ecological risk assessment from the Red, Book ancestor. Hum. Ecol. Risk Assess. 2003, 9, 1349–1360. [Google Scholar] [CrossRef]
- Sahle, M.; Yeshitela, K. Dynamics of land use land cover and their drives study for management of ecosystems in the socio-ecological landscape of Gurage Mountains, Ethiopia. Remote Sens. Appl. Soc. Environ. 2008, 12, 48–56. [Google Scholar] [CrossRef]
- Wang, H.; Liu, X.; Zhao, C.; Chang, Y.; Liu, Y.; Zang, F. Spatial-temporal pattern analysis of landscape ecological risk assessment based on land use/land cover change in Baishuijiang National nature reserve in Gansu Province, China. Ecol. Indic. 2021, 124, 107454. [Google Scholar] [CrossRef]
- Thompson, B.S.; Friess, D.A. Stakeholder preferences for payments for ecosystem services (PES) versus other environmental management approaches for mangrove forests. J. Environ. Manag. 2019, 233, 636–648. [Google Scholar] [CrossRef] [PubMed]
- Mo, W.; Wang, Y.; Zhang, Y.; Zhuang, D. Impacts of road network expansion on landscape ecological risk in a megacity, China: A case study of Beijing. Sci. Total Environ. 2017, 574, 1000–1011. [Google Scholar] [CrossRef] [Green Version]
- Xu, D.; Hou, G. The Spatiotemporal Coupling Characteristics of Regional Urbanization and Its Influencing Factors: Taking the Yangtze River Delta as an Example. Sustainability 2019, 11, 822. [Google Scholar] [CrossRef] [Green Version]
- 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]
- Xie, W.; Huang, Q.; He, C.; Zhao, X. Projecting the impacts of urban expansion on simultaneous losses of ecosystem services: A case study in Beijing, China. Ecol. Indic. 2018, 84, 183–193. [Google Scholar] [CrossRef]
- Hunsaker, C.T.; Graham, R.L.; SuterII, G.W. Assessing ecological risk on a regional scale. Environ. Manag. 1990, 14, 325–332. [Google Scholar] [CrossRef]
- Chen, S.Q.; Chen, B.; Fath, B.D. Ecological risk assessment on the system scale: A review of state-of-the-art models and future perspectives. Ecol. Modeling 2013, 250, 25–33. [Google Scholar] [CrossRef]
- Liang, J.; Li, S.; Li, X.; Li, X.; Liu, Q.; Meng, Q.; Lin, A.; Li, J. Trade-off analyses and optimization of water-related ecosystem services (WRESs) based on land use change in a typical agricultural watershed, southern China. J. Clean. Prod. 2021, 279, 123851. [Google Scholar] [CrossRef]
- Li, J.; Pu, R.; Gong, H.; Luo, X.; Ye, M.; Feng, B. Evolution characteristics of landscape ecological risk patterns in Coastal Zones in Zhejiang Province, China. Sustainability 2017, 9, 584. [Google Scholar] [CrossRef] [Green Version]
- Wang, X.; Zheng, D.; Shen, Y. Land use change and its driving forces on the Tibetan Plateau during 1990–2000. Catena 2008, 72, 56–66. [Google Scholar] [CrossRef]
- Chen, L.D.; Fu, B.J.; Zhao, W.W. Source- sink landscape theory and its ecological significance. Acta Ecol. Sin. 2006, 26, 1444–1449. (In Chinese) [Google Scholar] [CrossRef]
- Fu, B.J. The integrated studies of geography: Coupling of patterns and processes. Acta Geogr. Sin. 2006, 69, 1052–1059. (In Chinese) [Google Scholar]
- Li, W.; Cao, Q.; Lang, K.; Wu, J.S. Linking potential heat source and sink to urban heat island: Heterogeneous effects of landscape pattern on land surface temperature. Sci. Total Environ. 2017, 586, 457–465. [Google Scholar] [CrossRef]
- Fan, J.; Wang, Y.; Zhoum, Z.; You, N.S.; Mengm, J.J. Dynamic ecological risk assessment and management of land use in the middle reaches of the Heihe River based on landscape patterns and spatial statistics. Sustainability 2016, 8, 536. [Google Scholar] [CrossRef] [Green Version]
- Forbes, V.E.; Calow, P. Developing predictive systems models to address complexity and relevance for ecological risk assessment. Integr. Environ. Assess. Manag. 2013, 9, e75–e80. [Google Scholar] [CrossRef] [Green Version]
- Wang, J.; Cui, B.S.; Liu, J.; Yao, H.R.; Cui, H.J. The effect of land use and its change on ecological risk in the Lancang River watershed of Yunan Province at the landscape scale. Acta Sci. Circumstantiae 2008, 28, 269–277. (In Chinese) [Google Scholar]
- Malekmohammadi, B.; Rahimi, B.L. Ecological risk assessment of wetland ecosystems using multi-criteria decision making and geographic information system. Ecol. Indic. 2014, 41, 133–144. [Google Scholar] [CrossRef]
- Liu, S.L.; Cui, B.S.; Dong, S.K.; Yang, M.; Holt. Evaluating the influence of road networks on landscape and regional ecological risk: A case study in Lancang River Valley of Southwest China. Ecol. Eng. 2008, 34, 91–99. [Google Scholar] [CrossRef]
- Stow, C.A.; Borsuk, M.E.; Stanley, D.W. Long-term changes in watershed nutrient inputs and riverine exports in the Neuse River, North Carolina. Water Res. 2001, 35, 1489–1499. [Google Scholar] [CrossRef]
- Xu, X.B.; Yang, G.S.; Tan, Y.; Zhuang, Q.; Li, H.; Wan, R.; Su, W.; Zhang, J. Ecological risk assessment of ecosystem services in the Taihu Lake Basin of China from 1985 to 2020. Sci. Total Environ. 2016, 554, 7–16. [Google Scholar] [CrossRef]
- Gong, J.; Zhao, C.X.; Wang, H.L.; Sun, P.; Xie, Y.C.; Meng, X.M. Ecological risk assessment of Longnan Mountainous Area based on geological disasters: A case study of Wudu. J. Mt. Sci. 2012, 30, 570–577. (In Chinese) [Google Scholar]
- Li, W.; Wang, Y.; Xie, S.; Sun, R.; Cheng, X. Impacts of landscape multifunctionality change on landscape ecological risk in a megacity, China: A case study of Beijing. Ecol. Indic. 2020, 117, 106681. [Google Scholar] [CrossRef]
- Pili, S.; Serra, P.; Salvati, L. Landscape and the city: Agro-forest systems, land fragmentation and the ecological network in Rome, Italy. Urban For. Urban Green. 2019, 41, 230–237. [Google Scholar] [CrossRef]
- Pablo, L.J.; Sarda, R.; Jimenez, J.A. A methodological framework for multi-hazard risk assessment in beaches. Environ. Sci. Policy 2011, 14, 685–696. [Google Scholar] [CrossRef]
- Zhang, F.G.; Peng, M.; Wang, H.Y.; Ma, H.H.; Xu, R.T.; Cheng, X.M.; Hou, Z.L.; Chen, Z.W.; Li, T.; Cheng, H.X. Ecological Risk Assessment of Heavy Metals at Township Scale in the High Background of Heavy Metals, Southwestern, China. Environ. Sci. 2020, 41, 4197–4209. (In Chinese) [Google Scholar]
- Zhang, Y.; Li, Y.Z.; Lv, J.; Wang, J.; Wu, Y. Scenario simulation of ecological risk based on land use/cover change—A case study of the Jinghe county. China Ecol. Indic. 2021, 131, 108176. [Google Scholar] [CrossRef]
- Zhuang, D.F.; Liu, J.Y. A study on the regional variation model of land use degree in China. J. Nat. Resour. 1997, 12, 10–16. (In Chinese) [Google Scholar]
- Yan, E.P.; Lin, H.; Wang, G.X.; Xia, C.Z. Analysis of evolution and driving force of ecosystem service values in the Three Gorges Reservoir region during 1990—2011. Acta Ecol. Sin. 2014, 34, 5962–5973. (In Chinese) [Google Scholar]
- Wu, J.G. Landscape Ecology, Patterns, Processes, Scales and Grades, 2nd ed.; Higher Education Press: Beijing, China, 2007. [Google Scholar]
- Xie, H.L. Spatial characteristic analysis of land use eco-risk based on landscape structure: A case study in the Xingguo County, Jiangxi Province. China Environ. Sci. 2011, 31, 688–695. (In Chinese) [Google Scholar]
- Qin, F.T.; Peng, Y. Temporal and spatial characteristics of ecological risk in Shunyi, Beijing, China based on landscape structure. Chin. J. Appl. Ecol. 2016, 27, 1585–1593. (In Chinese) [Google Scholar]
- Wang, P.; Wang, Y.J.; Liu, X.P.; Chen, X.; Kong, F.X. Ecological risk assessment of an ecological migrant resettlement region based on landscape structure: A case study of Hongsibu in Ningxi. Acta Ecol. Sin. 2018, 38, 2672–2682. (In Chinese) [Google Scholar]
- Zhao, W.Q.; Yang, Z.H.; Su, W.C.; Li, K. Ecological risk assessment and management of watershed based on landscape pattern change- A case study of the Chishui river basin in Guizhou. Resour. Environ. Yangtze Basin 2017, 26, 1218–1227. [Google Scholar]
- Xu, X.G.; Lin, H.P.; Fu, Z.Y. Regional ecological risk assessment of wetland in the Huang River Delta. Acta Sci. Entiarum Natualum Unversit Is Pekinesis 2001, 37, 111–120. [Google Scholar]
- He, S.S.; Li, X.; He, C.L.; Fang, B. Landscape Ecological Risk Assessment in Guangling District of Yangzhou City Based on Land Use Change. J. Nanjing Norm. Univ. Nat. Sci. Ed. 2019, 42, 138–148. (In Chinese) [Google Scholar]
- Lu, W.D. SPSS for Windows Statistical Analysis (Version 2); Electronic Industry Press: Bejing, China, 2002. [Google Scholar]
- Gao, Y.N.; Zhao, M.L.; Wang, H.L.; Xiong, M.; Zhao, T.Q. Analysis of spatial differentiation of grassland ecological quality and its affecting factors based on landscape ecology perspective: A case study of Siziwang Banner. Acta Ecol. Sin. 2019, 39, 5288–5300. (In Chinese) [Google Scholar]
- Fotheringham, A.S.; Brunsdon, C.; Charlton, M. Geographically Weighted Regression: The Analysis of Spatially Varying Relationships; John Wiley & Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
- Brunsdon, C.; Fotheringham, A.S.; Charlton, M.E. Geographically weighted regression: A method for exploring spatial nonstationarity. Geogr. Anal. 1996, 28, 281–298. [Google Scholar] [CrossRef]
- He, L.; Shen, J.; Zhang, Y. Ecological vulnerability assessment for ecological conservation and environmental management. J. Environ. Manag. 2018, 206, 1115–1125. [Google Scholar] [CrossRef]
- Galic, N.; Schmolke, A.; Forbes, V.; Baveco, H.; Van den Brink, P.J. The role of ecological models in linking ecological risk assessment to ecosystem services in agroecosystems. Sci. Total Environ. 2012, 415, 93–100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Huang, M.Y.; Yue, W.Z.; Fang, B.; Feng, S.R. Scale response characteristics and geographic exploration mechanism of spatial differentiation of ecosystem service values in Dabie Mountain area, central China from 1970 to 2015. Acta Geogr. Sin. 2019, 74, 1904–1920. (In Chinese) [Google Scholar]
- Liu, C.Y.; Zhang, K.; Liu, J.P. A long-term site study for the ecological risk migration of landscapes and its driving forces in the Sanjiang Plain from 1976 to 2013. Acta Ecol. Sin. 2018, 38, 3729–3740. (In Chinese) [Google Scholar]
- Rang, S.H.; Tan, M.H.; Lv, C.H. Ecological Risk Analysis of LUCC in Relation to Stakeholders. Prog. Geogr. 2010, 29, 439–444. (In Chinese) [Google Scholar]
- Wang, Y.F.; Ni, Z.K.; Li, X.X.; Wang, S.R. The ecological risk assessment and stage characteristics of Dongting Lake. China Environ. Sci. 2019, 39, 321–329. (In Chinese) [Google Scholar]
- Ni, H.J.; Nan, Y.; Zhu, W.H.; Cui, Y.X.; Liu, G.M.; Liu, C.; Yao, K. Study on comprehensive assessment of ecological security in Changbai Mountain Region based on multi-hazard natural disasters risk. Geogr. Res. 2014, 33, 1348–1360. [Google Scholar]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).