Next Article in Journal
Assessment of the Knowledge, Perceptions, and Reactions towards the African Apefly (Spalgis lemolea lemolea) in Tanzania
Next Article in Special Issue
Effects of Aquatic Macrophytes on Spatial Distribution and Feeding Habits of Exotic Fish Species Lepomis macrochirus and Micropterus salmoides in Shallow Reservoirs in South Korea
Previous Article in Journal
Employee Commitment Matters for CSR Practice, Reputation and Corporate Brand Performance—European Model
Previous Article in Special Issue
The Potential of Gully Erosion on the Yamal Peninsula, West Siberia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Evolution of Landscape Ecological Risk Based on Geomorphological Regionalization during 1980–2017: A Case Study of Shaanxi Province, China

1
College of Urban and Environmental Science, Northwest University, Xi’an 710127, China
2
Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, Xi’an 710127, China
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(3), 941; https://0-doi-org.brum.beds.ac.uk/10.3390/su12030941
Submission received: 20 November 2019 / Revised: 18 January 2020 / Accepted: 21 January 2020 / Published: 28 January 2020
(This article belongs to the Collection Risk Assessment and Management)

Abstract

:
Land surface elements, such as land use, are in constant change and dynamically balanced, driving changes in global ecological processes and forming the regional differentiation of surface landscapes, which causes many ecological risks under multiple sources of stress. The landscape pattern index can quickly identify the disturbance caused by the vulnerability of the ecosystem itself, thus providing an effective method to support the spatial heterogeneity of landscape ecological risk. A landscape ecological risk model based on the degree of interference and fragility was constructed and spatiotemporal differentiation of risk between 1980 and 2017 in Shaanxi Province was analyzed. The spatiotemporal migration of risk was demonstrated from the perspective of geomorphological regionalization and risk gravity. Several conclusions were drawn: The risk of Shaanxi Province first increased and then decreased, at the same time, the spatial differentiation of landscape ecological risk was very significant. The ecological risk presented a significant positive correlation but the degree of autocorrelation decreased. The risk of the Qinba Mountains was low and the risk of the Guanzhong Plain and Han River basin was high. The risk of Loess Plateau and sandstorm transition zone decreased greatly and their risk gravities shifted to the southwest. The gravity of the Guanzhong Plain and Qinling Mountains had a northward trend, while the gravity of the Han River basin and Daba Mountains shifted to the southeast. In the analysis of typical regions, there were different relationships between morphological indicators and risk indexes under different geomorphological features. The appropriate engineering measures and landscape management for different geomorphological regionalization were suggested for effective reduction of ecological risks.

1. Introduction

Land use is a comprehensive reflection of human economic and social activities acting on terrestrial surface resources and the natural environment. It is the result of the combined action of human activities and natural processes [1]. As the most direct manifestation of human utilization of the natural environment [2], land use and other land surface elements (water, soil, atmosphere, biology, etc.) are in constant change and have a dynamic balance; thus, these processes drive the changes in global ecological processes and form the spatial differentiation of surface landscapes. The beginning of the “Anthropocene” is reflected by changes in global landscape patterns and ecological functions [3] and this new period is related to the multi-directional evolution of the ecological environment and has generated many ecological risks. The assessment and prevention of ecological risk has become a popular research topic in geography and ecology in response to the comprehensive management of social ecosystems [4,5,6,7,8].
Ecological risk is the possibility of adverse ecological effects due to the external pressures exerted on ecosystem components [9,10,11]. Ecological risk assessment (ERA) started in the 1980s and initially focused on the impact of chemical pollutants on human health [6,12]. Subsequently, due to the complexity of risk sources and the diversity of risk receptors, ERA has gradually extended to the category of socio-economic-natural complex ecosystems, including multiple risk sources and risk receptors [13]. The current research methods for assessing ecological risk include the perspective of risk source-sink and landscape structure [14,15]. Risk assessment based on source-sink can be used to identify risk sources and receptors in the area threatened by an exposure-response [16,17], and additional methods include the PSR framework [18] and the RRM model [19,20,21]. Risk assessment based on the landscape structure relies on the coupling and correlation perspective of patterns and processes and focuses more attention on the spatiotemporal heterogeneity of risks; this approach often uses the land use intensity parameters [2,22,23] and landscape loss degree [8,24,25] as the main elements. In the absence of ecological monitoring data, the landscape pattern characteristics reflected by land use can clarify the impact of quantitative change and spatial allocation of land use on the ecosystems and processes to reflect the spatial differentiation of the ecological environment [26,27,28]. Therefore, using the landscape of heterogeneous land units composed of different types of ecosystems as the perspective of evaluation, this paper explores the impact of human activities on land use based on the landscape pattern, providing the best generalized approach for further clarifying the evolutionary mechanism of ecological risks [5]. This method has been applied in many areas, such as watersheds [7,29], coastal zones [30] and administrative regions [25].
Land uses have spatiotemporal heterogeneity and are subject to regional topographic and geomorphic characteristics [31], which are associated with spatiotemporal changes in landscape pattern and ecological risks. This spatiotemporal variation is mainly manifested as the directional migration of ecological risks on terrain with different gradients [8]. Correlation studies show that the simple response of ecological risk to terrain is achieved using the methods of ecological risk distribution in the terrain grade area ratio [22,32], ecological risk index calculation on terrain grade [23], and terrain and ecological risk section line [8,33]. However, based on the unique and single geomorphologic features in a study area, it is difficult to discuss the influence of geomorphologic differentiation on ecological risk evolution in a large region. Therefore, this method may be an effective way to identify the variation trend of ecological risk and its related elements from the perspective of geomorphological regionalization.
Using Shaanxi Province, China, which has various geomorphic types, as an example, this paper aims to restore and explain the causal chain of “geomorphologic features-human drivers-pattern change-risk response” under the ecological risk paradigm of landscape pattern. At the same time, the spatiotemporal evolution model of risk is explored based on the perspective of geomorphological regionalization to promote the transition of regional ecosystems from vulnerable to sustainable through ecological restoration and landscape pattern optimization.

2. Study Area and Data Sources

2.1. Study Area

Shaanxi Province (105°29′ E–111°15′ E, 31°42′ N–39°35′ N) is located in central China (Figure 1). The administrative divisions include 10 cities, which account for a total catchment area of 20.57 × 104 km2. The average annual temperature of the entire province is 13.7 °C, and temperature decreases from south to north. The average annual precipitation is 680 mm, with more precipitation in the south and less precipitation in the north. Southern Shaanxi Province has a humid climate, Guanzhong Plain has a semi-humid climate, and northern Shaanxi Province has a semi-arid climate [34,35]. The Qinba Mountains in southern Shaanxi are adjacent to the Qinling Mountains in the north and the Daba Mountains in the south, and the Han River basin in located in the central part. The Qinba Mountains are an important water source in the South-north Water Diversion Project and an ecological function-restricted development zone. There are serious landslides, floods and other geological disasters in this mountainous area. The central part of the province is the Guanzhong Plain, which is high in the west and low in the east, high in the north and south, and low in the middle. The Guanzhong Plain is a major urban agglomeration area with a strong intensity of human activity. The geomorphologic types of northern Shaanxi include the sandstorm transition zone and the Loess Plateau (including tableland, beam, and gully) geomorphology. Precipitation is scarce and unevenly distributed in this region. It is a typical ecologically fragile region in the world, with a high drought frequency and serious soil erosion. The spatial heterogeneity of human disturbances caused by the various geomorphic differentiation provides an ideal case study for discussing the spatiotemporal differentiation of eco-logical risk from the perspective of geomorphological regionalization.

2.2. Data Sources

The land use data of Shaanxi Province in 1980, 2000 and 2017 were derived from the 1:100,000 vector data set of land use provided by the Cold and Arid Regions Sciences Data Center of China (http://westdc.westgis.ac.cn). These data were based on Landsat MSS, TM and ETM images, which were generated by combining visual interpretation of terrain and cognition of image spectrum. The accuracy of data interpretation was 95%, which was the highest precision land use data product in China. According to the classification standard of land use of China (GB/T21010-2017), land use types were divided into six categories: farmland, forest, grassland, water body, construction land, and unused land. Land use data were resampled to 30-m grids for subsequent analysis. The digital elevation model (DEM) was constructed from data sourced from the Geospatial Data Cloud (http://www.gscloud.cn/), and the elevation of the study area was extracted. The boundary of geomorphological regionalization was vectorized from the geographic atlas of Shaanxi Province, and the subarea types included the sandstorm transition zone, Loess Plateau, Guanzhong Plain, Qinling Mountains, Han River basin and Daba Mountains. Then, the boundary vector files were generated.

3. Research Methods

3.1. The Spatial Display Approach

Regarding ecological risk spatialization, the equidistant sampling grid method was the most common type [15]. In landscape ecology, when the area of a landscape sample (that is, the sampling unit) is 2–5 times the area of a landscape type patches, it can comprehensively reflect the landscape pattern information around the sampling point [36,37]. Based on this criterion, a grid size of 12 km was selected as a suitable scale for the spatialization of ecological risks. The study regions were divided into 1471 cell grids and the risk values were assigned to the corresponding risk point sets. As a regional variable, the ecological risk index can be measured by spatial evolution trends using the geostatistical semi-variance function [28,30]. The formula was given as follows:
γ ( h ) = 1 2 N ( h ) i = 1 N ( h ) [ Z ( x i + h ) Z ( x i ) ] 2
where γ(h) is the semi-variance function; h is the step length; Z(xi + h) and Z(xi) are the observed value of the ecological risk index at spatial location xi and xi + h, respectively; and N(h) is the number of samples within the interval distance h.
In this paper, GS+7.0 software was used to achieve the best fit of the semi-variance function [30]. Combined with the model fitting effect, the spatial analysis of ecological risk in 1980 was based on the spherical model for fitting and the spatial analysis of ecological risk in 2000 and 2017 was based on the exponential model. Finally, a spatial distribution map of the ecological risk index was created by ordinary Kriging interpolation using the Geostatistics Analyst Tools of ArcGIS 10.2.

3.2. Landscape Index Method

The landscape index model with the landscape disturbance indexes and the landscape vulnerability indexes coupled with the area weighting function as the main components method is widely used to quantify the characteristics of landscape patterns, providing an effective method to clarify the spatial heterogeneity of landscape ecological risk; this method is the main approach of landscape ERA [38,39]. In this paper, the ecological risk index model was established to describe the relative value of the ecological risk index in the evaluation unit. The formula is as follows:
E R I k = i = 1 n S k i S k E i × F i
where ERIk denotes the landscape ecological risk index of an unit k; n is the number of landscape types; Ei and Fi are the disturbance indexes and fragility indexes of landscape i, respectively; Ski is the area of landscape component i in unit k; and Sk is the total area of unit k. By referring to previous studies, the detailed calculation methods of Ei are given in Table 1 [40].
Landscape fragility represented the vulnerability of the ecosystem structure represented by each land use type. The level of vulnerability can reflect the sensitivity and resistance of the landscape to interference by external risk. The greater the vulnerability of a landscape is, the greater the ecological risk is, and vice versa. By referring to relevant studies [8,41] and using the expert scoring method, the fragility of land use was divided into 6 grades (Table 2).

3.3. Spatial Autocorrelation Method

Spatial autocorrelation is the basic form of spatial dependence, and this value reflects the degree of spatial correlation of the risk variables. Global spatial autocorrelation is used to measure the global spatial distribution characteristics of risk value. The Moran’s I index can measure the risk value from the degree of global spatial correlation and spatial difference [2,42]. The formula is as follows:
I ( d ) = n i = 1 n j = 1 n w i j ( x i x ¯ ) ( x j x ¯ ) i = 1 n ( x i x ¯ ) i = 1 n j = 1 n w i j
where xi and xj are the values of variable x at adjacent space points, respectively; x ¯ is the average of the variables; wij is the adjacent weight; and N is the total number of risk points. The value of Moran’s I index is [−1, 1]. When Moran’s I < 0, there is a negative correlation. When Moran’s I > 0, there is a positive correlation. When Moran’s I = 0, there is no correlation. Finally, the autocorrelation is stronger as the absolute value approaches 1.
Local anomalies can be introduced into the local spatial correlation index LISA to identify the spatial aggregation pattern of cold spots and hot spots. The formula is as follows:
I i = Z i j = 1 n w i j Z j , ( i j )
where Ii is the unit value of the LISA index, Zi and Zj are the standardized risk values of units i and j, respectively, and wij is the spatial weight matrix.
In this paper, GeoDA software was used to measure the spatial distribution pattern of the ecological risk index in three periods in Shaanxi Province. A Moran’s I scatter plot was generated, and the Moran’s I values and Z values of the standardized statistics were extracted (p = 0.01).

3.4. Standard Deviational Ellipse and Risk Gravity

The standard deviation ellipse (SDE) is widely used in the field of spatial statistics, and it is a common method used to quantitatively analyze the overall characteristics of the spatial distribution of point elements. As spatial variables, the global characteristics of the spatial distribution of risk points are presented from multiple angles using the SDE and its gravity tool [27,43]. The expression of gravity is as follows:
X = i = 1 n w i x i i = 1 n w i ; Y = i = 1 n w i y i i = 1 n w i
where (X, Y) is the weighted average gravity, and (xi, yi) represents the coordinates of risk points. wi represents the weight. In this paper, the direction distribution module of ArcGIS10.2 was used to calculate the risk gravity of geomorphological regionalization, with the risk value as the weight, and to visualize the space.

4. Results

4.1. Spatial and Temporal Differentiation of Landscape Ecological Risk

The ecological risk value was classified using the ArcGIS10.2 natural breakpoint method [44]: grade I (ERI < 0.1190), grade II (0.1190 ≤ ERI < 0.1271), grade III (0.1271 ≤ ERI < 0.1353), grade IV (0.1353 ≤ ERI < 0.1445), and grade V (ERI > 0.1445); this method enabled us to generate the spatial distribution of ecological risk in Shaanxi Province from 1980 to 2017 (Figure 2). The landscape ecological risk index of Shaanxi Province was 0.1282, 0.1287 and 0.1273, respectively, and the risk value first increased and then decreased, with an overall decrease of 0.78%. The spatiotemporal differentiation of ecological risk was very significant in Shaanxi Province. An analysis of Figure 2 shows that the grade I-II risk area in 1980 accounted for approximately half of the study area, and this area was mainly distributed in the Qinba Mountains and the southern part of the Loess Plateau. The grade IV-V risk area accounted for 19.25% of the study area, and this area was mainly distributed in the northern Loess Plateau, Guanzhong Plain and Han River basin. Compared with 1980, the grade IV risk area in the northern part of the Loess Plateau increased significantly. At the same time, the grade V risk in the Guanzhong Plain extended to the north and south. The urban built-up areas of the Han River basin formed the core of the grade V risk area. Compared with 2000, the ecological risk in 2017 decreased significantly. The grade IV risk area in the north of the Loess Plateau was cut off by the grade III risk area. The grade V risk area in the south of Xi’an has been cut off, which formed a double core area of high risk in Xianyang city and Weinan city. The grade IV-V risk area of the Han River basin has been greatly reduced, and the grade V risk area has been withdrawn from the Hanzhong urban construction area.

4.2. Spatial Autocorrelation Analysis

From 1980 to 2017, the Moran’s I index of ecological risk in Shaanxi Province was 0.7279, 0.7277 and 0.7053, respectively, and the Z-values of the random distribution test were 37.9032, 38.2940 and 35.5271, respectively. This finding indicates that the ecological risk presented a significant positive correlation. The Moran’s I index showed a downward trend over time, indicating that the spatial autocorrelation degree of ecological risk had been reduced. The landscape pattern and ecological process were affected by the evolution of land use types, and these changes affected the landscape ecological risk, leading to the gradual transition of its spatial distribution from aggregation to uniformity.
Based on the LISA index, GeoDA software was used to analyze the clustering characteristics of landscape ecological risk in Shaanxi Province (Figure 3), and the result was significant at the p ≤ 0.05 level. There is a large number of “high–high” and “low–low” agglomeration areas in Figure 3a. Risk hotspots were mainly concentrated in the north of the Loess Plateau, the sandstorm transition zone and the central Guanzhong Plain. The risk hotspots in 2000 showed an increasing trend, especially in the Loess Plateau. However, the area of risk hotspots in this region decreased to some extent in 2017, and the significance of spatial agglomeration weakened. The risk cold spots were mainly located in the Qinling Mountains and the eastern region of Yan’an, and there was little change in the risk cold spots during the study period, which was closely related to the fact that most of the cold spot areas were located in semi-natural and natural areas. In addition, there were some spatial “singularity spots” in the research region. The “low–high” spots were located in the periphery of risk hotspots, and these sites were greatly disturbed by human activities and had farmland that was easily converted into other land use types. The change in land use types affected the habitat vulnerability and affected the regional landscape pattern, which made it more possible for this part of the region to transform into high risk. The “high–low” spots were located around the risk cold spots, which may be caused by the unused land that was surrounded by forest and grassland. In general, the singularity spots were located in the transitional region of land use change and landscape pattern evolution and they were also the areas where ecological risk was most likely to spread.

4.3. Temporal Variation of Ecological Risk Based on Geomorphologic Regionalization

The mean value of ecological risk in geomorphologic regionalization was extracted using the Partition Statistical Tool in ArcGIS10.2, and the rate of the temporal variation in ecological risks was calculated (Figure 4). Overall, the changes in ecological risk values were obvious in the different geomorphological regionalizations.
(1)
The Qinling Mountains (ranked sixth) and the Daba Mountains (ranked fifth) were always low-value areas of ecological risk. The degree of risk evolution of the Daba Mountains was more obvious than that of the Qinling Mountains. Mountainous areas with higher topographical relief greatly restrict farming behaviors while increasing the expansion cost of urban construction, which results in a low intensity of human disturbance. This scenario led to a high degree of dominance of the regional forest and grassland, forming a relatively closed and safe biological refuge, thus maintaining low ecological risk during the study period.
(2)
The internal spatial difference of ecological risk in the Loess Plateau (ranked fourth) was very obvious. The north of the Loess Plateau is the hilly and gully region with higher topographical relief. The disordered landscape structures and broken ecological corridors were formed in this region, resulting in a high level of ecological risk. The forest coverage of the southern part of the Loess Plateau is second only to that of the Qinba Mountains; thus, the ecological risk remained relatively low. Overall, the ecological risk of the Loess Plateau increased from 1980 to 2000, which may reflect the changes in land use types and spatial patterns in the context of agricultural farming and township construction. The rate of decline of ecological risk in the Loess Plateau ranked third in terms of the risk change from 2000 to 2017, which was related to the national policy of returning farmland to forests implemented in 1999. The ecological policy of China increased the area of forests and grassland, reduced the degree of separation of forests and grassland landscapes, greatly reducing the spread trend of ecological risk.
(3)
The sandstorm transition zone was located in the northwest of Shaanxi Province (ranked second), where natural vegetation is rare and wind erosion and desertification of soil are very serious. This condition led to a high concentration of unused land, high habitat fragility and high ecological risk. The ecological risk in this region tended to increase from 1980 to 2000. However, a large amount of unused land was converted into forest and grassland after the policy of returning farmland to forest was implemented and this policy enhanced the regional ability to resist sand erosion and effectively reduced the regional habitat sensitivity. The risk reduction rate in this region was the first in terms of the risk change in the corresponding period, with the risk decreasing by 3.22%.
(4)
The Han River basin (ranked third) and the Guanzhong Plain (ranked first) were always high-risk areas, but the rates of risk change are fourth and third, respectively. The Guanzhong Plain and Han River basin were the main urban construction areas and important irrigated agricultural areas. Urban construction land and farmland were the most typical man-made landscapes and their land use patterns and boundary fractal characteristics were more marked by human interference, affecting the evolution direction of the landscape pattern. In the early period of disturbance (1980–2000), urban construction land patches were abundant, forming fragmented and disordered landscape areas, which increased the ecological risk. In the later period of interference (2000–2017), the space of construction land expanded rapidly. The construction land with a large area gradually swallowed up the surrounding smaller construction land patches, or many small patches gradually merged to form larger construction land patches, which greatly enhanced the stability of the regional landscape and clearly reduced the ecological risk.

4.4. Spatial Evolution of Ecological Risk Based on Geomorphologic Regionalization

The SDE of geomorphological regionalization was generated and the risk gravity was extracted based on the direction distribution tool of ArcGIS10.2. The first-order standard deviation was selected as the model parameter and the ellipse size covered 68% of the risk value (Figure 5). As shown in Figure 5, the shifting direction of ecological risk gravity can be roughly divided into three categories: the southwest, northeast and southeast directions.
(1)
The areas where the risk gravity shifted to the southwest included the sandstorm transition zone and the Loess Plateau. Agricultural production, the policy of returning farmland to forests, and the development of mineral resources should be considered as key factors in the ecological risk changes. This region is important energy accumulation places for coal, oil resources in China. The ecological risk was higher under the influence of a fragile natural basement and excessive energy exploitation. With the transformation of energy structure, the petroleum exploitation in the west of Yulin city strengthened, while the traditional coal cities in the north were at a disadvantage during the industrial transformation, which directly affected the evolution direction of local landscape structure. In addition, as a result of the policy implementation of the three-north shelterbelt construction, a large amount of forest and grassland formed stable habitats in the northern area, which led to the southward shift of risk gravity.
(2)
The areas where the risk gravity shifted to the northeast included the Guanzhong Plain and the Qinling Mountains. Population expansion brought about by urban expansion and agricultural production should be a core factor in risk changes in the Guanzhong Plain. The risk changes in Qinling Mountains was related to the establishment of protected areas. The Wei River ran roughly from southwest to northeast, which was basically consistent with the direction of the risk gravity shift. The upper reaches of the Wei River have complicated topography, limited farming activities and urban expansion; thus, human activities in the Guanzhong Plain tended to expand to the lower reaches of the Wei River due to improved transportation during the study period. As the national central park, human disturbance acted on the northern and southern slopes of the Qinling Mountains to varying degrees. The northern Qinling Mountains faced a higher degree of human disturbance due to their proximity to the Guanzhong Plain [45], which caused the risk gravity of the Qinling Mountains to move northward.
(3)
The areas where the risk gravity shifted to the southeast included the Han River basin and the Daba Mountains. Urban expansion and agricultural production were also considered in risk changes of Han River basin and the risk changes in Daba Mountains were related to the policy of returning farmland to forests. During the research period, the construction land of Hanzhong tended to expand to the downstream regions of the Han River, which destroyed the integrity of the farmland patch. The urban progress of Ankang city was slower than that of Hanzhong city, but the risk increased due to urban expansion. This pattern shifted the risk gravity to the southeast. The increase in the large amounts of forest and grassland areas in the western part of the Daba Mountains further stabilized the ecological process, which significantly reduced the risk in the western part of the Daba Mountains and shifted the risk gravity to the southeast.

5. Discussion

5.1. Geomorphologic Features, Human Drivers, Pattern Change, Risk Response

The restriction of geomorphologic features on human activities and the natural selectivity of the spatial distribution of land use made human disturbance vary in different geomorphological regionalization scales and the change in related elements, such as landscape structure and habitat fragility, led to the change in ecological risk value in different geomorphological regionalization. In combination with Figure 2a,c, the spatial distribution of human activities under geomorphologic features and the resulting pattern changes and the risk response were discussed by using the south of Shenmu City and the south of Xi’an city, which have obvious risk changes, as an example (Figure 6).
(1)
The restricted type: Shenmu City is located in the typical loess hilly gully region with large topographic relief. Human activity was forced to separate because of the mixed terrain, which simultaneously formed a singular and repeated landscape patch. The farmland was located in terraces or gullies and most of the abandoned farmland was located in the high-slope area. A large amount of farmland and unused land was converted into grassland during the period of farmland regression, which greatly enhanced the dominance of grassland patches and the surrounding stable habitats, benefiting the restoration of regional habitats and effectively reducing ecological risk. In the ecological risk study of the Loess Plateau, Wang et al. believed that the ecological risk of the resettlement area increased first and then decreased from 1995 to 2015 [25], which is similar to the results of this study. However, Ren et al. concluded that the ecological risk on the Loess Plateau has generally increased and declined locally based on comprehensive risk models since 1990 [46]. In short, it is undeniable that the policy of returning farmland to forests has effectively reduced the severity of ecological risks in this region.
(2)
The beneficial type: As a typical plain city, the topography of Xi’an was less restrictive to human activities and formed diverse and complex landscape patches. The spatial pattern was dominated by farmland at the beginning of the period and changed to become co-dominated by farmland and construction land at the end of the period. Xi’an city mainly expanded to the southwest during the study period. Chang’an District in the south of Xi’an was dominated by “enclave” patches at the beginning of the period and was gradually “submerged” by spreading expansion. The occupation of farmland by construction land gradually evolved into a more stable large-scale landscape base, thereby increasing the stability of the regional landscape. In the study of the risk of plain cities, Shi et al. believes that the risk of the urban-rural junction is gradually decreasing due to the consolidation of building land patches [47], which is consistent with the characteristics of the ecological risk impact of urban expansion in Xi’an. At the same time, Wang et al. suggested that plain urban agglomerations have always had higher ecological risks due to severe human interference [27]. The risk change of plain cities is closely related to the degree of urban development. Human disturbance can influence the evolution direction of landscape patterns in the plain area, thus leading to an increase or decrease in ecological risk.
In order to better reflect the relationship between pattern changes and ecological risks, Patch Density (PD), Edge Density (ED), Fractal Dimension (FD), CONTAG, DIVISION, and SHDI were used based on the Fragstats4.2 software in typical regions [27,48,49]. The calculation result is shown in Table 3. In Table 3, the ecological risks in typical regions both decreased. At the same time, the PD and DIVISION indexes also decreased, while the CONTAG index increased. This is consistent with the results of related ecological risk studies [49,50]. In contrast, ED, FD, and SHDI have different trends in different regions. These morphological indicators reveal the characteristics of urban morphology [48] and landscape resilience [51], but the accuracy of these indicators for ecological risk characterization needs to be further demonstrated.

5.2. Risk Reduction from the Perspective of Geomorphological Regionalization

Based on the perspective of geomorphological regionalization, this can provide a basis for regional ecological risk reduction and sustainable development in the context of habitat restoration and pattern optimization.
(1)
The Loess Plateau and sandstorm transition zone are ecologically fragile regions in China. In the process of farming and urban construction, the landscape layout and three-dimensional structure should be reasonably designed according to the original landform, precipitation condition and distribution of forests and rivers [13]. The naturalness of the basement edge of the ecosystem should be maintained while the fragmentation of the landscape should be reduced. The policy of returning farmland to forests should be further implemented, as this policy promotes the restoration of fragile habitats and reduces the sensitivity of regional habitats, thus slowing the spread of ecological risk.
(2)
The Guanzhong Plain and Han River basin were the main areas of urban concentration in Shaanxi Province. The terrain struggled to restrict the expansion of Xi’an and other plain cities; thus, the subjective initiative of people can be greatly developed. Efforts should be made to build the ecological space structure of urban agglomerations, and the ecological security pattern was constructed through the optimization and reorganization of the regional ecological matrix, patches and corridors based on the perspective of landscape ecology to effectively reduce ecological risk [52].
(3)
The Qinba Mountains are part of a national park and represent an ecological function restriction development zone in China. A new eco-environmental protection regulation on the Qinba Mountains was formulated for vegetation protection, biodiversity maintenance and construction constraints in 2019. Vegetation protection relies on the policy of returning farmland to forests. The expansion of forests and the improvement of connectivity are positive changes in forest landscapes [53], which helps to better exert water conservation and climate regulation. Wildlife habitat reserves have been set up to avoid severe disturbance of biological habitats by human activities. Construction constraints are mainly coping with spatial conflicts with local habitats arising from transportation corridors, rural settlements, and mineral resource development. Landscape management is implemented through corridor design, rural reconstruction, and mineral development governance, thereby reducing ecological risks.

5.3. Limitations of Study

The spatiotemporal variation of ecological risk in Shaanxi Province was explored based on the landscape pattern index as a proxy indicator of ecological risk. This method is effective in the absence of long-term regional monitoring data, but it also has some limitations. Using only landscape pattern indicators to characterize ecological risk was incomplete; therefore, indicators of natural disasters and human disturbance risk sources should be given more attention. At the same time, using landscape value [54] and ecosystem services [13] as indicators of the degree of ecological risk loss can more clearly describe the comprehensive loss of landscape and ecosystem under ecological risk stress. The design and implementation of an integrated framework of risk sources and losses is the focus of future ecological risk research. Scale selection is one of the core issues of landscape ecology. Although 12 km was a suitable scale for workload and ecological risk expression, it may not be the optimal scale. Therefore, analysis on the scale effect of landscape ecology on granularity and amplitude will be one of the tasks for the author in the future research. Moran ’s I and LISA are applied to determine the optimal scale to improve the scientificity of the ecological risk assessment results. The change of ecological risk is affected by many factors and the relationship between these factors and risk should be observed and quantified to explore the main stress factors affecting the change of risk in order to better assist risk management.

6. Conclusions

Landscape ERA models based on landscape patterns indexes can effectively reflect the evolution of land use patterns and changes in ecological environments. Depending on the perspective of geomorphological regionalization, the spatial transition of ecological risk can be further clarified. From 1980 to 2017, the landscape ecological risk index of Shaanxi Province was 0.1282, 0.1287 and 0.1273, respectively, with an overall decrease of 0.78%. The Moran’s I index was 0.7279, 0.7277 and 0.7053, respectively. The ecological risk presented a significant positive correlation, and the degree of autocorrelation decreased to some extent. The mountains maintained a low ecological risk and a small change rate during the study period. Due to the spatial distribution of building land and arable land, the plain area faced a more serious risk level. Benefiting from the national policy of returning farmland to forests, the ecological risk of the Loess Plateau and sandstorm transition zone greatly reduced and the risk gravity of them gradually retreated southward. The risk gravity of the Guanzhong Plain and Qinling Mountains tended to move northward, while the risk gravity of the Han River basin and the Daba Mountains shifted to the southeast. In the analysis of typical regions, different reasons for the reduction of ecological risks were explored. At the same time, there were different relationships between morphological indicators and risk indexes under different geomorphological features. Therefore, the ability of some morphological indicators to characterize risks should be further verified. Finally, appropriate engineering measures and landscape management for different geomorphological areas were suggested for effective reduction of ecological risks.

Author Contributions

D.L., H.C., and H.Z. designed the research and wrote the paper. T.G., Q.S. analyzed the data. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (Grant No.41671086 & 41871185) and the APC was funded by Grant No.41671086.

Acknowledgments

Supports from the National Natural Science Foundation of China are gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, Y.; Xu, X.; Tang, Q. Human activity intensity of land surface: Concept, methods and application in China. J. Geogr. Sci. 2016, 26, 1349–1361. [Google Scholar]
  2. Zhou, R.J.; Zhang, Y.Z.; He, H.C. Ecological risk assessment based on land use changes in the coastal area in Yancheng city. Geogr. Res. 2016, 35, 1017–1028. [Google Scholar]
  3. Zheng, D.; Wu, S.H.; Yin, Y.H.; Yang, Q.Y.; Zhao, D.S.; Zhang, X.Q. Frontiers in terrestrial system research in China under global change. Acta. Geogr. Sin. 2016, 71, 1475–1483. [Google Scholar]
  4. Chang, H.F.; Peng, J.; Wang, Y.L.; Wu, W.B.; Yang, P.; Liu, Y.X.; Song, Z.Q.; Hsueh, I.C. Characteristics of landscape system in the context of global change. Progr. Geogr. 2017, 36, 1167–1175. [Google Scholar]
  5. Peng, J.; Du, Y.Y.; Liu, Y.X.; Wu, J.S.; Wang, Y.L. From natural regionalization, land change to landscape service: The development of integrated physical geography in China. Geogr. Res. 2017, 36, 1819–1833. [Google Scholar]
  6. Kicińska, A.; Mamak, M. Health risks associated with municipal waste combustion on the example of Laskowa commune (Southern Poland). Hum. Ecol. Risk Assess. 2017, 23, 2087–2096. [Google Scholar]
  7. Tian, P.; Li, J.L.; Gong, H.B.; Pu, R.L.; Cao, L.D.; Shao, S.Y.; Shi, Z.Q.; Feng, X.L.; Wang, L.J.; Liu, R.Q. Research on land use changes and ecological risk assessment in Yongjiang River Basin in Zhejiang Province, China. Sustainability 2019, 11, 2817. [Google Scholar]
  8. Liu, D.; Liang, X.Y.; Chen, H.; Zhang, H.; Mao, N.Z. A quantitative assessment of comprehensive ecological risk for a loess erosion gully: A case study of Dujiashi Gully, northern Shaanxi Province, China. Sustainability 2018, 10, 3239. [Google Scholar]
  9. Faber, J.; Van Wensem, J. Elaborations on the use of the ecosystem services concept for application in ecological risk assessment for soil. Sci. Total Environ. 2012, 415, 3–8. [Google Scholar]
  10. 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]
  11. Li, P.; Wu, J. Sustainable living with risks: Meeting the challenges. Hum. Ecol. Risk Assess. 2019, 25, 1–10. [Google Scholar]
  12. Mirzabeygi, M.; Abbasnia, A.; Yunesian, M.; Nodehi, R.N.; Mahvi, A.H. Heavy metal contamination and health risk assessment in drinking water of Sistan-and-Baluchistan, Southeastern Iran. Hum. Ecol. Risk Assess. 2017, 23, 1893–1905. [Google Scholar]
  13. Du, Y.Y.; Peng, J.; Zhao, S.Q.; Hu, Z.C.; Wang, Y.L. Ecological risk assessment of landslide disasters in mountainous areas of Southwest China: A case study in Dali Bai Autonomous Prefecture. Acta Geogr. Sin. 2016, 71, 1544–1561. [Google Scholar]
  14. Peng, J.; Dang, W.X.; Liu, Y.X.; Zong, M.L.; Hu, X.X. Review on landscape ecological risk assessment. Acta Geogr. Sin. 2015, 70, 664–677. [Google Scholar]
  15. Cao, Q.W.; Zhang, X.W.; Ma, H.K.; Wu, J.S. Review of landscape ecological risk and an assessment framework based on ecological services: ESRISK. Acta Geogr. Sin. 2018, 73, 843–855. [Google Scholar]
  16. Yang, P.; Mao, X.L.; Li, T.H.; Gao, X.W. Ecological Risk Assessment of the Shenzhen River-Bay Watershed. Hum. Ecol. Risk Assess. 2011, 17, 580–597. [Google Scholar]
  17. Malekmohammadi, B.; Blouchi, L.R. Ecological risk assessment of wetland ecosystems using multi criteria decision making and geographic information system. Ecol. Indic. 2014, 41, 133–144. [Google Scholar]
  18. Zhou, Q.; Zhang, X.; Wang, Z. Land use ecological risk evaluation in Three Gorges Reservoir Area based on normal cloud model. Trans. Chin. Soc. Agric. Eng. 2014, 30, 289–297. [Google Scholar]
  19. Bartolo, R.; Van Dam, R.; Bayliss, P. Regional ecological risk assessment for Australia’s tropical rivers: Application of relative risk model. Hum. Ecol. Risk Assess. 2012, 18, 16–46. [Google Scholar]
  20. Kanwar, P.; Bowden, W.B.; Greenhalgh, S. A regional ecological risk assessment of the Kaipara Harbour, New Zealand, using a relative risk model. Hum. Ecol. Risk Assess. 2015, 21, 1123–1146. [Google Scholar]
  21. Heenkenda, M.K.; Bartolo, R. Regional ecological risk assessment using a relative risk model: A case study of the Darwin Harbour, Darwin, Australia. Hum. Ecol. Risk Assess. 2016, 22, 401–423. [Google Scholar]
  22. Xu, Y.; Zhong, Y.; Feng, X.H.; Xu, L.T.; Zheng, L. Ecological risk pattern of Poyang Lake basin based on land use. Acta Ecol. Sin. 2016, 36, 7850–7857. [Google Scholar]
  23. Hu, J.L.; Zhou, Z.X.; Teng, M.J.; Luo, N. Ecological risk assessment of typical karst basin based on land use change: A case study of Lijiang River basin, Southern China. Chin. J. Appl. Ecol. 2017, 28, 2003–2012. [Google Scholar]
  24. Liu, Y.C.; Liu, Y.X.; Li, J.L.; Lu, W.Y.; Wei, X.L.; Sun, C. Evolution of landscape ecological risk at the optimal scale: A case study of the open Coastal Wetlands in Jiangsu, China. Int. J. Environ. Res. Public Health 2018, 15, 1691. [Google Scholar]
  25. 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 Ningxia. Acta Ecol. Sin. 2018, 38, 2672–2682. [Google Scholar]
  26. Paukert, C.P.; Pitts, K.L.; Whittier, J.B.; Olden, J.D. Development and assessment of a landscape-scale ecological threat index for the Lower Colorado River Basin. Ecol. Indic. 2011, 11, 304–310. [Google Scholar]
  27. Wang, F.C.; Wang, D.C.; Zhang, L.H.; Liu, J.Y.; Hu, B.X.; Sun, Z.C.; Chen, J.H. Spatiotemporal analysis of the dynamic changes in land use ecological risks in the urban agglomeration of Beijing-Tianjin-Hebei region. Acta Ecol. Sin. 2018, 38, 4307–4316. [Google Scholar]
  28. Zhao, Y.; Luo, Z.J.; Li, Y.T.; Guo, J.Y.; Lai, X.H.; Song, J. Study of the spatial-temporal variation of landscape ecological risk in the upper reaches of the Ganjiang River Basin based on the “production-living-ecological space”. Acta Ecol. Sin. 2019, 39, 4676–4686. [Google Scholar]
  29. Xie, H.L.; Wang, P.; Huang, H.S. Ecological risk assessment of land use change in the Poyang Lake eco-economic zone, China. Int. J. Environ. Res. Public Health 2013, 10, 328–346. [Google Scholar]
  30. Li, J.L.; Pu, R.L.; Gong, H.B.; Ye, M.Y.; Feng, B.X. Evolution characteristics of landscape ecological risk patterns in coastal zones in Zhejiang Province, China. Sustainability 2017, 9, 584. [Google Scholar]
  31. Ha, K.; Ding, Q.L.; Men, M.X.; Xu, H. Spatial distribution of land use and its relationship with terrain factors in hilly area. Geogr. Res. 2015, 34, 909–921. [Google Scholar]
  32. Zhao, Y.J.; Li, Y.B.; Shao, J.A. Ecological risk assessment of small watershed of the Three Gorges reservoir area based on land use change: A case study of Caotang River. J. Nat. Resour. 2013, 28, 944–956. [Google Scholar]
  33. Peng, J.; Xie, P.; Liu, Y.X.; Hu, X.X. Integrated ecological risk assessment and spatial development trade-offs in low-slope hilly land: A case study in Dali Bai Autonomous Prefecture, China. Acta Geogr. Sin. 2015, 70, 1747–1761. [Google Scholar]
  34. Chen, H.; Lo’pez-Carr, D.; Tan, Y.; Liang, X.Y. China’s Grain for Green policy and farm dynamics: Simulating household land-use responses. Reg. Environ. Chang. 2016, 16, 1147–1159. [Google Scholar]
  35. Chen, H.; Marter-Kenyon, J.; Lo’pez-Carr, D.; Liang, X.Y. Land cover and landscape changes in Shaanxi Province during China’s Grain for Green Program (2000–2010). Environ. Monit. Assess. 2015, 187, 644. [Google Scholar]
  36. Su, H.M.; He, A.X. Analysis of land use based on RS and geostatistics in Fuzhou City. J. Nat. Resour. 2010, 25, 91–99. [Google Scholar]
  37. Qing, 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. [Google Scholar]
  38. Focks, A.; Ter Horst, M.; Van Den Berg, E.; Baveco, H. Integrating chemical fate and population-level effect models for pesticides at landscape scale: New options for risk assessment. Ecol. Model. 2014, 280, 102–116. [Google Scholar]
  39. Zhang, T.; Liu, Y.X.; Peng, J.; Wang, Y.L. Correlation of the landscape ecological risk on multi-scales in Shenzhen City. Chin. J. Ecol. 2016, 35, 2478–2486. [Google Scholar]
  40. Chang, X.Y.; Li, X.J.; Li, X.C.; Guo, P.; Gao, F. Spatial-temporal heterogeneity of ecological risk of land use in mining areas. Acta Ecol. Sin. 2019, 39, 3075–3088. [Google Scholar]
  41. Hua, L.Z.; Liao, J.F.; Chen, H.X.; Chen, D.K. Assessment of ecological risks induced by land use and land cover changes in Xiamen City, China. Int. J. Sustain. Dev. World Ecol. 2018, 25, 439–447. [Google Scholar]
  42. Cao, Y.H.; Chen, C.; Zhang, D.P.; Liu, M.Y.; Dong, S.S. Evolution of ecological risk pattern of land use change in Wanjiang City Belt. Acta Ecol. Sin. 2019, 39, 4773–4781. [Google Scholar]
  43. 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. [Google Scholar]
  44. Wang, F.; Ye, C.S.; Hua, J.Q.; Li, X. Coupling relationship between urban spatial expansion and landscape ecological risk in Nanchang City. Acta Ecol. Sin. 2019, 39, 1248–1262. [Google Scholar]
  45. Guo, S.Z.; Bai, H.Y.; Meng, Q.; Huang, X.Y.; Qi, G.Z. Landscape pattern change and its response to anthropogenic disturbance in the Qinling Mountains during 1980 to 2015. Chin. J. Appl. Ecol. 2018, 29, 4080–4088. [Google Scholar]
  46. Ren, Z.Y.; Sun, Y.J.; Wu, L.X. Comprehensive evaluation of ecological risk of land during 1990–2012 in Yulin City. Chin. J. Eco-Agric 2017, 25, 656–664. [Google Scholar]
  47. Shi, H.P.; Yu, K.Q.; Feng, Y.J. Ecological risk assessment of rural-urban ecotone based on landscape pattern: A case study in Daiyue District of Tai’an City, Shandong Province of East China. Chin. J. Appl. Ecol. 2013, 24, 705–712. [Google Scholar]
  48. Song, Y.Z.; Long, Y.; Wu, P.; Wang, X.Y. Are all cities with similar urban form or not? Redefining cities with ubiquitous points of interest and evaluating them with indicators at city and block levels in China. Int. J. Geogr. Inf. Sci. 2018, 32, 2447–2476. [Google Scholar]
  49. Pei, H.; Wei, Y.; Wang, X.Y.; Qin, Z.H.; Hou, C.L. Method of cultivated land landscape ecological security evaluation and its application. Trans. Chin. Soc. Agric. Eng. 2014, 30, 212–219. [Google Scholar]
  50. Yu, X.; Wu, K.N.; Yun, W.J.; Yun, W.J.; Wei, H.B.; Liu, L.; Song, Y.H.; Gao, X. Analysis on temporal and spatial variation of landscape ecological security in modern agricultural area. Trans. Chin. Soc. Agric. Eng. 2016, 32, 253–259. [Google Scholar]
  51. Zhang, H.; Liang, X.Y.; Liu, D.; Shi, Q.Q.; Chen, H. The resilience evolution and scenario simulation of social-ecological landscape in the fragile area. Acta Geogr. Sin. 2019, 74, 1450–1466. [Google Scholar]
  52. Yang, T.R.; Kuang, W.H.; Liu, W.D.; Liu, A.L.; Pan, T. Optimizing the layout of eco-spatial structure in Guanzhong urban agglomeration based on the ecological security pattern. Geogr. Res. 2017, 36, 441–452. [Google Scholar]
  53. Solecka, I.; Raszka, B.; Krajewski, P. Landscape analysis for sustainable land use policy: A case study in the municipality of Popielów, Poland. Land Use Policy 2018, 75, 116–126. [Google Scholar]
  54. Sylla, M.; Solecka, I. Highly valued agricultural landscapes and their ecosystem services in the urban-rural fringe–an integrated approach. J. Environ. Plan. Manag. 2019, 62, 1–29. [Google Scholar]
Figure 1. Location and DEM of Shaanxi Province, China.
Figure 1. Location and DEM of Shaanxi Province, China.
Sustainability 12 00941 g001
Figure 2. Spatial distribution of landscape ecological risk in Shaanxi Province during 1980–2017.
Figure 2. Spatial distribution of landscape ecological risk in Shaanxi Province during 1980–2017.
Sustainability 12 00941 g002
Figure 3. Spatial distribution of LISA in Shaanxi Province during 1980–2017.
Figure 3. Spatial distribution of LISA in Shaanxi Province during 1980–2017.
Sustainability 12 00941 g003
Figure 4. Temporal changes in landscape ecological risk based on geomorphological regionalization in Shaanxi Province during 1980–2017.
Figure 4. Temporal changes in landscape ecological risk based on geomorphological regionalization in Shaanxi Province during 1980–2017.
Sustainability 12 00941 g004
Figure 5. Temporal changes in gravity of landscape ecological risk based on geomorphological regionalization in Shaanxi Province during 1980–2017.
Figure 5. Temporal changes in gravity of landscape ecological risk based on geomorphological regionalization in Shaanxi Province during 1980–2017.
Sustainability 12 00941 g005
Figure 6. The response of landscape ecological risk of typical regions to land use change.
Figure 6. The response of landscape ecological risk of typical regions to land use change.
Sustainability 12 00941 g006
Table 1. Calculation of landscape disturbance indexes.
Table 1. Calculation of landscape disturbance indexes.
Landscape Disturbance IndexesCalculation MethodMeaning
Landscape fragmentation index (Ci) C i = n i / A i Ci refers to the process in which landscape types change from single, homogeneous, and continuous to complex, heterogeneous and discontinuous due to the interference of nature or human beings. The increase in Ci indicates a decrease in system stability and an aggravation of ecological risk.
Landscape isolation index (Ni) N i = A 2 A i n i A Ni refers to the degree of separation of different patch numbers in landscape types. The increase in the Ni index indicates that the landscape is more scattered, the landscape distribution is more complex, and the fragmentation degree is higher.
Landscape dominant index (DOi) D O i = Q i + M i 4 + L i 2 DOi is an indicator used to measure the importance of patches in a landscape, which directly reflects the impact of this landscape type on landscape pattern.
Landscape disturbance index (Ei) E i = a C i + b N i + c D i Ei reflects the disturbance degree of an ecosystem represented by different landscape types.
Note: ni is the number of patches in the ith landscape; Ai is the area of the ith landscape; A is the landscape total area; Qi is the sample number of the ith pattern/total number of samples; Mi is the number of ith patches/total number of patches; Li is the area of the ith patches/area of samples; a, b and c are the weighted coefficients, and a + b + c = 1. As suggested in the literature [28,30], a = 0.5, b = 0.3 and c = 0.2 in this study.
Table 2. The grades of fragility and corresponding values of fragility for the different land use types.
Table 2. The grades of fragility and corresponding values of fragility for the different land use types.
Land Use TypesGrades of FragilityValues of Fragility
Construction land10.0476
Forest20.0952
Grassland30.1429
Farmland40.1905
Waterbody50.2381
Unused land60.2857
Table 3. The relationship between morphological indicators and ecological risks in typical regions.
Table 3. The relationship between morphological indicators and ecological risks in typical regions.
Typical RegionsYearMorphological IndicatorsEcological Risk
PDEDFDCONTAGDIVISIONSHDI
Shenmu City19800.679558.12401.127059.06540.93741.13140.1364
20170.657253.64291.125761.07000.87021.07740.1302
Xi’an City19800.794423.54211.052959.39620.87021.06460.1129
20170.745924.95871.059360.31340.82701.11870.1044

Share and Cite

MDPI and ACS Style

Liu, D.; Chen, H.; Zhang, H.; Geng, T.; Shi, Q. Spatiotemporal Evolution of Landscape Ecological Risk Based on Geomorphological Regionalization during 1980–2017: A Case Study of Shaanxi Province, China. Sustainability 2020, 12, 941. https://0-doi-org.brum.beds.ac.uk/10.3390/su12030941

AMA Style

Liu D, Chen H, Zhang H, Geng T, Shi Q. Spatiotemporal Evolution of Landscape Ecological Risk Based on Geomorphological Regionalization during 1980–2017: A Case Study of Shaanxi Province, China. Sustainability. 2020; 12(3):941. https://0-doi-org.brum.beds.ac.uk/10.3390/su12030941

Chicago/Turabian Style

Liu, Di, Hai Chen, Hang Zhang, Tianwei Geng, and Qinqin Shi. 2020. "Spatiotemporal Evolution of Landscape Ecological Risk Based on Geomorphological Regionalization during 1980–2017: A Case Study of Shaanxi Province, China" Sustainability 12, no. 3: 941. https://0-doi-org.brum.beds.ac.uk/10.3390/su12030941

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