Next Article in Journal
Spatiotemporal Patterns and Risk Zoning of Wildfire Occurrences in Northeast China from 2001 to 2019
Previous Article in Journal
Responses of Vegetation Phenology to Urbanization in Plateau Mountains in Yunnan, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Landscape Fragmentation and Spatial Autocorrelation of a Typical Watershed in the Wenchuan Earthquake-Affected Area—A Case Study in the Longxi River Basin

1
College of Forestry, Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
College of Natural Resources and Environment, Northwest Agriculture and Forestry University, Xianyang 712100, China
3
Institute of Soil and Water Conservation, Northwest Agriculture and Forestry University, Xianyang 712100, China
4
Key Laboratory for Forest Ecosystem Process and Management of Fujian Province, Fuzhou 350002, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Forests 2023, 14(12), 2349; https://doi.org/10.3390/f14122349
Submission received: 27 September 2023 / Revised: 20 November 2023 / Accepted: 24 November 2023 / Published: 29 November 2023
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)

Abstract

:
The 2008 8.0 Ms Wenchuan earthquake caused serious ecological degradation in the northwest of China’s Sichuan Province. Elucidating landscape fragmentation and spatial pattern of post-earthquake landscape is of great interest in order to improve ecological restoration and predict the spatial distribution of damaged ecosystems in earthquake-affected areas. We used four real-time remote sensing images to analyze the landscape pattern characteristics and spatial autocorrelation of the Longxi River Basin before and after the Wenchuan earthquake during the period from 2005 to 2015. In the study period, the degree of landscape fragmentation increased 1 year after the earthquake in 2009 but decreased in 2011 and 2015. The spatial distribution of forest, farmland, and shrub-grassland had significantly positive spatial correlation and the characteristics of spatial aggregation from 2005 to 2015. Construction land had no significant spatial correlation before the earthquake in 2005, but significantly positive spatial correlation after the earthquake, while traffic land had no significant spatial correlation across the study period. Unlike the other landscape types, geological disasters changed from a significantly negative spatial correlation before the earthquake to a significantly positive spatial correlation after the earthquake. However, the positive autocorrelation of all landscape types decreased with the increase of space distance, but in different distance-decay rates. The High-high spatial aggregation areas of geological disasters and construction land clustered gradually upstream of the basin and distributed in Longchi town, respectively; those of farmland distributed in the southeast of the basin increased in 2009 but then decreased, while those of forests had the opposite pattern as the dominant landscape type after the earthquake. Our results elucidated the spatial structure and distribution features of the Longxi River Basin to give a theoretical foundation for assessing the effects of ecological restoration and reconstruction management in earthquake-affected areas.

1. Introduction

Major earthquakes can produce a large amount of geo-hazards from available land resources. This change was always accompanied by heavy economic losses [1], serious ecological degradation [2], and landscape fragmentation [3]. To improve the restoration process in earthquake-affected areas, many researchers used the integrated techniques of 3S technology combined with geography and landscape ecology theories to study the spatiotemporal dynamic characteristics of land use affected by secondary disasters induced by earthquakes [4,5,6,7,8]. They also reported the automatic detection of earthquake-triggered landscape changes [9], landscape evolution [10], and changes in ecosystem services in the earthquake-stricken area [11]. For example, the Jiuzhaigou earthquake occurred on 8 August 2017, with a magnitude of 7.0, which caused the overall landscape grade of Jiuzhaigou Nature Reserve to decrease by 15% (Luo et al., 2023) [12]. In the 2004 Mw 6.6 Chuetsu earthquake-affected areas in Japan, natural succession improved landscape structure, resulting in good vegetation recovery (85.5%) by 2021 (Xiang et al., 2023) [13]. Their works help us predict geomorphological evolution, landscape patterns, and ecosystem functioning and processes [14,15,16,17]. However, little is known about the consequence of landscape fragmentation caused by the earthquake on a small scale, leaving a great need for an ecological restoration process and landscape pattern optimization after an earthquake.
In order to describe the dynamic of landscape patterns and its control factors and predict possibilities of landscape change trajectory, researchers set up a dynamic model to establish robust diagnostics, prediction, and simulation through mathematical models with regard to logical systems and flow interactions among various spatial elements [18]. At present, cellular automata (CA) modeling, based on multi-criteria evaluation (MCE) and logistic regression (LR), is used widely to investigate landscape dynamics due to its ability to capture dynamic information and predict future scenarios [19,20]. Classic statistical methods including MCE and LR assumed individual observations to be statistically independent, while spatial landscape data have the tendency to be dependent, which was defined as spatial autocorrelation—a phenomenon that random variables affected by geographic proximity would have more similar or less similar values over distance than randomly associated pairs of observations [21,22].
Spatial autocorrelation can be used to quantify the spatial dependency in univariate data, including Moran’s index (Moran, 1950) and Geary’s coefficient (Geary, 1954), and multivariate data, such as a Mantel correlogram (Legendre and Legendre, 1998), across a geo-referenced space. Though some researchers debate the methodological disadvantage of the spatial autocorrelation analysis [21], it was still applied to describe the spatial variability of land use or landscape [21,23,24,25]. For example, to determine the law of landscape pattern evolution, scholars analyzed the geographical correlation of the landscape index [26], the change in land cover characteristics [27], the dynamic evolution characteristics of habitat quality [28], and the correlation features with interference elements [29] using spatial autocorrelation analysis. However, previous studies mainly concerned the spatial autocorrelation of landscape pattern change in a large-scale river basin but did not consider how landscape patterns respond to earthquakes in a small watershed.
The Longxi River Basin located in Dujiangyan City, Sichuan Province is considered a typical basin in Wenchuan earthquake-affected areas [10,30]. The Ms 8.0 Wenchuan earthquake caused frequent secondary disasters, including debris flows, landslides, collapses, and floods, as well as subsequent deterioration in the ecological environment of the Longxi River Basin [31]. The changes caused by secondary disasters in the geomorphology and river system induced a large area of vegetation destruction and reshaped land use patterns [32]. Our previous study showed that unutilized land induced by secondary disasters replaced grassland as the second dominant land type, accounting for 15.99% of the Longxi River Basin’s area in 2009 [8]. However, to our knowledge, the spatial autocorrelation analysis in a typical basin after the earthquake is still poorly understood and undocumented. Moreover, the Longxi River Basin belongs to the Wolong Nature Reserve, which supports approximately 10% of the total wild giant panda population across the globe [3] and plays an important role in the safe operation of the Zipingpu Reservoir and water supply for the Chengdu Plain [31]. Therefore, it is necessary to elucidate landscape fragmentation and the spatial pattern of post-earthquake landscapes.
In this study, we collected four different high-resolution remote-sensing satellite images of interpreted land use types in the Longxi River Basin. The aim of this study was to determine how changes in the spatial structure and distribution features were affected by the earthquake and subsequent recovery. Here, we hypothesized that natural recovery could reduce the degree of landscape fragmentation and change spatial autocorrelation characteristics of landscape types after the earthquake. Our main objectives were to (1) analyze the characteristics of the basin’s landscape fragmentation before and after the earthquake, (2) determine the change of spatial correlation of landscape types under ecological recovery after the earthquake, and (3) demonstrate spatial aggregation characteristics of landscape types. This study could provide a theoretical foundation and scientific recommendations for the restoration and reconstruction of fragile habitats in the Longxi River Basin.

2. Materials and Methods

2.1. Study Area

The Longxi River Basin is located on the upper reach of the Minjiang River along the eastern margin of the Qinghai–Tibet Plateau in the northwest of Dujiangyan City, Sichuan Province, China (103°30′4.98″–103°36′2.08″ E, 31°1′28.06″–31°10′49.53″ N). Its elevation ranges from 822 m to 3291 m a.s.l. with a maximum east–west distance of 6.88 km and north–south distance of 17.32 km, and an area of 79.24 km2. The topography of the basin is dominated by moderate and high mountains with an intensive surface incision from north to south, which leads to a “V”-shaped canyon [30]. According to WRB classification (Hempel et al.) [33], the basin has three soil types: ali-udic argosols are found in the zone at an elevation of 822 to 1600 m a.s.l., haplic acrisols are found in the zone at an elevation of 1600 to 2200 m a.s.l., and haplic luvisols are found in the zone at an elevation of 2200 to 3291 m a.s.l. The climate is tropical humid monsoon with concentrated precipitation (approximately 1100 mm per year), which occurs between May and October (Ding et al.) [34]. The length of the basin’s main channel is 15.7 km, with an average slope gradient of >20°, and the gully gradient is between 287‰ and 626‰ [30]. Due to frequent floods and debris flows after the Wenchuan earthquake, the basin is one of the three major demonstration sites for typical debris flow monitoring and early warning in Sichuan Province (Figure 1 and Figure 2).

2.2. Data

The remote sensing images used in this study were SPOT-5 in 2009, Quick Bird in 2005 and 2011, and Worldview-2 in 2015, supplied by the Institute of Mountain Hazards and Environment, Chinese Academy of Sciences (CAS). The spectral resolution and information of the images are shown in Table 1. All images were covered less than 2% by clouds. We simulated blue band of the SPOT-5 image in trained fitting model in Erdas 9.2 to receive its natural color image [35]. Then, we used ENVI 5.3 version software to correct spatial incompatibility across four images on the basis of World Geodetic System (WGS) 1984 coordinate system. The second-order polynomial model was used on the basis of primary scale Chinese national topographic maps (1:50,000) to ensure the accuracy of all images’ geometric correction within only one pixel with 24 ground control points. The method of human–computer interactive interpretation was used to interpret land use types with the reference of field survey data and unmanned aerial vehicle (UAV) images in 2018 in Longxi River Basin. Moreover, Shuttle Radar Topographic Mission (SRTM) with 30 m × 30 m resolution was used to produce the Digital Elevation Model (DEM) and extract topographic features.

2.3. Methodology

2.3.1. Landscape Types Classification System

According to the classification system created by Jia et al. (2003) [36], we divided the land use types of Longxi River Basin into six categories as cultivated land, forest, construction land, water, shrub-grassland, and unutilized land (mainly including mountains damaged by natural disasters, bare rock mass, bare land, etc.). The process of land use category identification, field observation, pixel-based supervised image classification, and classification accuracy assessment could be seen in Tian et al. (2022) [8]. The overall classification accuracy and the Kappa coefficient were 95.2% and 0.896, respectively.
Based on the classification of land use types, we classified them into nine landscape types, including forest land, farmland, shrub-grassland, construction land, traffic land, bare land, gully channel, geological hazard, and water according to the standard created by You et al. (2006) [37]. This helped us to analyze the change in each landscape in detail.

2.3.2. Landscape Pattern Index

The landscape pattern index is a simple quantitative indicator and one of the important indicators that highly summarizes the composition of landscape structure and spatial pattern characteristics. Landscape indices can describe the landscape pattern, establish the connection between landscape structure and processes or phenomena, and better explain and understand landscape functions. In this study, we used Fragstats 4.2 to assess the spatial pattern of landscape types [38]. To describe the characteristics of landscape fragmentation, the following landscape metrics, including edge density (ED), Shannon diversity index (SHDI), and landscape division index (DIVISION), were chosen according to the method described in You et al. (2006) [37]. ED represents the edge length of a patch per hectare, which determines the density of the patch. SHDI is a sensitive measurement indicator for the uneven distribution of various patch types in a landscape. If SHDI = 0, it indicates that the entire landscape is composed of one patch. The larger value indicates that there are more patch types with the uniform distribution of each patch type. DIVISION refers to the degree of fragmentation in the distribution of individuals with different patch numbers in a certain landscape type. Their formulas are as follows:
ED = k = 1 m e ik A ( 10000 )
where eik is the total edge length of patch type i, and A is the total landscape area. A smaller ED value indicates that the patch type is only adjacent to a few other types, while the higher value indicates the higher degree of landscape fragmentation, more significant edge effect, and stronger openness.
SHDI = i = 1 n P i ln P i
where Pi is the proportion of patch type i in the landscape. The larger index indicates the richer the types and distribution of various patches in the landscape.
D I V I S I O N = 1 j = 1 n a j A
where aj is the area of a certain type of patch j in the landscape type. The larger division index indicates the more severe the fragmentation of landscape types.

2.4. Landscape Pattern Spatial Autocorrelation

Moran’s I, Moran scatter plot, and Local Indicators of Spatial Association (LISA) were selected in this study to analyze and test spatial autocorrelation at the global and local levels using ArcGIS 10.2 and Geoda 1.6.7 software [39]. In detail, global and local Moran’s I can assess the degree of spatial association, which was visualized by Moran scatter plot and LISA to reveal the correlations and difference among the attributes of study units of the Longxi River Basin [40]. In this study, we used the grid analysis method to divide the Longxi River Basin into 374 grid cells of uniform size of 500 m × 500 m based on the principle of optimal sample square described by Fu et al. (2017) [41]. Each landscape type was superimposed with the grid cells for analysis except water because it occupies a concentrated area, resulting in insufficient patch numbers for autocorrelation analysis.

2.4.1. Global Spatial Autocorrelation Statistics

The global Moran’s I could reveal the spatial structure, distribution pattern, and their causes of regional variables [36]. Its formula is expressed as follows [42]:
Moran s   I = n i = 1 n i = j n W ij ( x i x ¯ ) ( x j x ¯ ) i = 1 n j = 1 n W ij i = 1 n ( x i x ¯ ) 2 ( i j )
where n is the total number of spatial units in the Longxi River Basin. xi and xj represent the attribute values of spatial elements of sites i and j, respectively. x is the average attribute values of spatial elements. Wij is a binary symmetric spatial weight matrix that indicates the distance or closeness between sites i and j. Distance-based spatial weight data were created using the queen’s standard with the distances of 0.5, 1, 2, 3, 4, 5, 6, 8, and 10 km, respectively. The value of Moran’s I generally ranges from −1 to 1. On this basis, the global spatial autocorrelation features and Moran’s I distance effects of each landscape type within the Longxi River Basin were examined between 2005 and 2015. Low negative values represent a negative spatial autocorrelation showing the clustering of dissimilar value, while high positive values represent a positive spatial autocorrelation showing the clustering of similar value. The significance of Moran’s I is usually assessed using the Z-score [43], and a Z-score > 1.96 indicates the significance at the level of p < 0.05.

2.4.2. Local Spatial Autocorrelation

Global autocorrelation could not determine the specific position of clustering locations or dispersed places due to the deficiency in aberrant local spatial pattern [44]. Therefore, we used LISA and the Moran scatter plot outlined by Anselin (1995) to assess the local spatial autocorrelation of attribute values and to identify the spatial distribution of anomalous regions. The formula calculating the local Moran’s I is as follows [45]:
I i = ( x i x ¯ ) i = 1 n ( x i x ¯ ) 2 n j = 1 W ij x i x ¯
where xi and xj are the normalized cellular observations from samples i and j; the other parameters are the same as the global Moran’s I.
Moran scatter diagram shown as Moran scatter plot can be used to analyze the spatial agglomeration characteristics between the attribute value in an observation and the average of its neighbors [46,47]. Moran scatter plot yields five types of spatial clusters, including not significant, High-high (HH), Low-low (LL), Low-high (LH), and High-low (HL) [48]. The type of “High-high” means that the attribute values of a spatial unit and its neighbors are both high, while that of “Low-Low” has the opposite pattern. Moreover, the type “High-Low” means that the attribute value of a spatial unit is high but its neighbors have low values, while that of “Low-High” has the opposite pattern [44,45,46,47,48,49,50,51].

3. Results

3.1. The Changes in Landscape Types

According to the actual situation of the Longxi River Basin, we divided the basin into nine landscape types, including forest, farmland, shrub-grassland, construction land, traffic land, bare land, gully channel, geological hazards, and water. As seen in Table 2, forest dominated all landscape types but decreased after the earthquake in 2009 and then recovered from 2009 to 2015, while shrub-grassland had a similar pattern but decreased during the period from 2011 to 2015. The area of farmland continually decreased, while that of traffic land had an opposite pattern during the study period. Bare land, gully channel, and geological hazards had similar patterns, in that their areas increased after the earthquake in 2009 but decreased from 2009 to 2015.

3.2. Landscape Fragmentation in the Longxi River Basin

Three landscape pattern indices, including edge density (ED), Shannon diversity (SHDI), and landscape isolation index (DIVISION), were selected to quantify the landscape fragmentation in the Longxi River Basin.

3.2.1. The Change in ED

ED is used to reveal the degree which a landscape or type is segmented by boundaries, showing a direct reflection of the degree of landscape fragmentation.
The change in ED can be intuitively seen from Figure 3. High-value ED areas were distributed along the Longxi River channel in 2005, while they spread over the basin after the Wenchuan earthquake in 2009. After natural recovery, the high-value ED area gradually reduced in 2011 and 2015, but it was still higher than in 2005.

3.2.2. The Change in SHDI

SHDI is a measurement index based on information theory and is widely used in ecology, reflecting landscape heterogeneity. It is particularly sensitive to the uneven distribution of various types in the landscape, emphasizing the contribution of rare types to information. Therefore, it is also used as a sensitive indicator when comparing and analyzing the diversity and heterogeneity changes of different landscapes or different periods of the same landscape.
As seen in Figure 4, the changes in SHDI in the Longxi River Basin also showed a similar trend from 2005 to 2015. The SHDI value along both sides of the Longxi River main channel was comparatively greater than the other parts in 2005, while it clustered in the tributary of the Longxi River after the Wenchuan earthquake in 2009. After natural recovery, the high-value SHDI area also gradually reduced in 2011 and 2015, but it was still higher than in 2005.

3.2.3. The Change in DIVISION

The DIVISION reflects the separation degree of patches in the landscape, indicating the complex degree of the landscape.
As seen in Figure 5, the high-value DIVISION areas increased dramatically after the Wenchuan earthquake in 2009 compared with those in 2005. After natural recovery, the high-value DIVISION areas still clustered in the upper reach of the Longxi River during the period from 2011 to 2015. In general, the high-value DIVISION areas changed significantly on the slopes of both sides of the Longxi River gully and its tributaries across the study period.
Overall, the degree of landscape fragmentation had a tendency of “increasing-reducing-reducing” in the whole basin across the study period, but it remained comparatively stable in the southwest of the basin.

3.3. Global Spatial Autocorrelation of Landscape Types

As shown in Table 3, Moran’s I values of forest, farmland, and shrub-grassland in the Longxi River Basin were all positive during the period from 2005 to 2015, indicating significantly positive spatial correlation characteristics (p < 0.01). Construction land had no significant spatial autocorrelation before the earthquake but had significantly positive spatial autocorrelation after the earthquake (2009–2015). Across the study period, traffic land had no significant spatial autocorrelation (p > 0.05). Bare land had significantly positive spatial autocorrelation except in 2009. Gully channel only had significantly positive spatial autocorrelation in 2005 before the earthquake. Geological hazards had significantly negative spatial autocorrelation in 2005 before the earthquake but had significantly positive spatial autocorrelation after the earthquake (2009–2015).
The global Moran’s I of landscape types except the distance of 0.5 km and geological hazards in 2005 in the Longxi River Basin gradually decreased with increasing distance from 2005 to 2015 (Figure 6). Within the lag distance of 10 km, the spatial autocorrelation of different landscape types had different characteristics across the lag distances before and after the earthquake. Before the earthquake in 2005, unlike the other landscape types, geological hazard had a significant negative spatial autocorrelation with the lowest Moran’s I (p < 0.01) at a distance of 0.5 km but had no significant spatial autocorrelation with the highest Moran’s I (p > 0.05) at a distance of 2 km. In addition, Moran’s I of construction land and traffic land had a lower value at a distance of 0.5 km than at a distance of 1 km. In 2009 (1 year after the earthquake), Moran’s I of traffic land, bare land, and gully channel had a lower value at a distance of 0.5 km than at a distance of 1 km. In 2011 (3 years after the earthquake), Moran’s I of traffic land and gully channel had a lower value at a distance of 0.5 km than at a distance of 1 km. In 2015 (7 years after the earthquake), only gully channel had a lower Moran’s I at a distance of 0.5 km than a distance of 1 km. Across the study period, Moran’s I of all landscape types had no significant spatial autocorrelation when the lag distance was >3 km.

3.4. Local Autocorrelation Analysis of Landscape Types

3.4.1. Moran Scatter Plot of Local Spatial Autocorrelation

Moran scatter plots were drawn using GeoDa, as seen in Figure 7 (in 2005), Figure 8 (in 2009), Figure 9 (in 2011), and Figure 10 (in 2015). In 2005, only geological hazards had a negative slope, and most of their points were located in quadrants II (LH) and IV (HL) in the scatter plot, indicating a negative spatial correlation. However, geological hazards had positive slopes and most of their points were located in quadrants I (HH) and III (LL), indicating positive spatial correlation after the earthquake (from 2009 to 2015). Across the study period, traffic land had no significant spatial correlation. Construction land had a low positive spatial correlation before the earthquake but a high positive spatial correlation after the earthquake, while bare land had the opposite pattern. The other four landscape types, forest, farmland, shrub-grassland, and gully channel, had comparatively stable Moran scatter plots with aggregation in HH and LL quadrants during the period from 2005 to 2015.

3.4.2. LISA Distribution of Local Spatial Autocorrelation

As seen in Table 4 and Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18, the spatial aggregation or anomaly distribution of landscape types in the Longxi River Basin changed dramatically during the period from 2005 to 2015. As the substrate landscape in the Longxi River Basin, forest had the highest and fluctuating HH aggregation area across the study period (15.64 km2 in 2005 accounting for 21.56%, 7.04 km2 in 2009 accounting for 11.33%, 10.76 km2 in 2011 accounting for 15.67%, and 12.71 km2 in 2011 accounting for 17.69%, respectively), but comparatively stable LL aggregation, LH anomaly, and HL anomaly area (Table 3 and Figure 11). Farmland had stable LL aggregation, LH anomaly, and HL anomaly area across the study period, except that the area of HH aggregation decreased to 0.04 km2, accounting for only 8.04% of the total area in 2015 (Table 3 and Figure 12). However, the HH aggregation area of construction land, gully channel, and geological hazards increased in 2009 after the earthquake, while the LL aggregation, LH anomaly, and HL anomaly area of those changed comparatively slightly (Table 3; Figure 14, Figure 17 and Figure 18). In addition, shrub-grassland and bare land had continuously decreasing HH aggregation area but also comparatively stable LL aggregation, LH anomaly, and HL anomaly area across the study period (Table 3; Figure 13 and Figure 16). Compared to the other landscape types, the spatial aggregation or anomaly distribution of traffic land slightly changed across the study period (Table 3 and Figure 15).

4. Discussion

4.1. Fragmentation of Landscape Patterns

Major earthquakes and their subsequent geological hazards caused huge damage to natural landscapes and ecosystems, resulting in heavy loss of life and property, ecological degradation, as well as landscape fragmentation [12,52,53]. It has been proven that landscape fragmentation had a high importance degree in assessing the extent of earthquake-induced landscape destruction due to its negative effects on biodiversity and the ecological process of species [12]. Our result showed that human activities and gully erosion across the Longxi River Basin were the main factors controlling landscape fragmentation in the basin before the earthquake in 2005. They caused high-value areas of edge density (ED), Shannon diversity (SHDI), and landscape isolation index (DIVISION), distributed along the banks of the Longxi River and concentrated in the southeast of the downstream due to human settlements. However, high-value areas of ED, SHDI, and DIVISION increased and spread over the basin due to landslides, collapses, and debris flows after the earthquake in 2009, indicating that a high degree of landscape fragmentation occurred in the Longxi River Basin. This study showed that high-value areas of ED, SHDI, and DIVISION decreased after 3 and 7 years of natural recovery (in 2011 and 2015), suggesting that the landscape of the basin had a great recovery potential due to the dominant natural extant forests managed strongly by the Longxi Hongkou National Nature Reserve administration [8]. This result may prove the hypothesis that natural extant forests are easier to restore after disturbances in comparison to plantations because of the even distribution of species, complicated community structure, and stable ecological equilibrium [54].

4.2. Landscape Spatial Global Autocorrelation

Our previous study showed that all vegetation types had significantly positive spatial autocorrelation in the most severely earthquake-affected areas of west Sichuan Province, China [55]. However, the spatial autocorrelation analysis was seldom used to characterize the spatial structure of landscape types on a small scale, such as a watershed or a basin.
In this study, whether the earthquake occurred or not, forest, farmland, and shrub-grassland had a significantly positive spatial correlation, showing spatial aggregation characteristics. Forest, as the first dominant landscape type, accounted for 91.53% before the earthquake in 2005, and still 78.42% after the earthquake in 2009. Its large area and wide-spread range caused positive spatial correlation across the study period. Shrub-grassland, as the second dominant landscape type, had a similar pattern. Farmland is concentrated in the area with flat terrain, an abundant source of water, and a dense population, which is mainly located in the southeast of the downstream of the Longxi River, so it also had a positive spatial correlation. In addition, though secondary disasters caused the reduction of farmland area, the demand for crop harvest and economic development with convenience resulted in the cluster of farmland.
Construction land had a significantly positive spatial correlation after the earthquake because the local government carried out the building reconstruction for the relocation of victims, resulting in an intensive distribution of settlements in the flat downstream alluvial fan, which is suitable for reconstruction in the basin. Though the construction of residential areas and the laying of cement land interrupted the natural environment and caused ecological damage, the centralized reconstruction after the earthquake greatly reduced the damage to the surrounding environment compared to the scattered distribution before the earthquake, resulting in the improvement of the local environment. This result is consistent with the findings of Shi et al. (2020) [56], who suggested the positive effect of reconstruction on construction land. It also proved that flat land in a mountainous area had the convenience of supporting a dense population (Zou et al., 2021) [57].
However, there is no significant spatial autocorrelation in traffic land across the study period because it is mainly affected by nonrandom human activities. Bare land had significantly positive spatial autocorrelation across the study period except 1 year after the earthquake in 2009. This may be due to the transition from bare land to geological disasters 1 year after the earthquake that caused the independent tendency over distance, resulting in less similar values for randomly associated pairs of observations [21]. This result is consistent with the findings of Li et al. (2012) [58], who observed that earthquake-induced disasters caused the random distribution of bare land.
Gully channels had a significantly positive spatial correlation before the earthquake in 2005 but not after the earthquake. This was due to frequent debris flows burying the original channels and forming headfall after the earthquake, resulting in the dispersed distribution of channels. In addition, the changes in river channels also caused road realignment and even flush back into roads and buildings. Therefore, an early warning system should be set up for disaster prevention and control. However, geological disasters changed from significantly negative spatial correlation before the earthquake to significantly positive spatial correlation after the earthquake, suggesting that they had high spatial dependency controlled by topographical factors (including elevation, slope gradient, distance to river) and seismic factors (including seismic intensity and distance to fault). This result is consistent with other studies demonstrating the positive relationship between topographical and seismic factors and geological disasters after an earthquake [59].
The Moran’s I of all landscape types decreased with distance and had no significant spatial autocorrelation when the lag distance was > 3 km. This means that the spatial autocorrelation of all landscape types could not be detected by aggregating the data when their aggregation level exceeded the level of increasing spatial scale (Overmars et al., 2003) [21]. This result is consistent with other studies demonstrating significant impacts of spatial scale on spatial structure [55].

4.3. Landscape Spatial Local Autocorrelation

Cluster and significance maps of LISA can show the interactions between sites and their neighbors in close proximity to identify the local spatial structure and instability by comparing their values [44]. In this study, forest had a high HH aggregation area in 2005 but the proportion of forest aggregation area decreased after the earthquake in 2009, then increased in 2011 and 2015, indicating that the spatial aggregation degree of forest decreased due to the disturbance of geological disasters but gradually recovered after natural restoration and succession.
Topographical factors are the main controlling factors of farmland distribution, while socioeconomic factors, including human tillage management, fertilization, and construction control the direction, quantity, and speed of farmland transition [17]. In this study, the distribution of farmland showed a trend of HH and LL aggregation but with a low area before the earthquake, suggesting that farmland resources in the Longxi River Basin had comparatively fewer spatial aggregation areas, distributed discretely. The HH aggregation area of farmland changed slightly one year after the earthquake but then decreased in 2011 and 2015, indicating that the spatial dependency of farmland diminished. This may be due to the following two reasons. First, frequent geological disasters such as landslides and debris flows induced by heavy rainfall in rainy seasons buried farmland distributed along the river [30,58]. Second, the local government converted farmland to forest and shrub-grassland or made a requisition of farmland due to its convenience for transportation, construction, and development for reconstruction for the relocation of victims, resulting in the reduction of the HH aggregation area of farmland [8].
The HH aggregation area of construction land increased from 2005 to 2009 but changed slightly during the period from 2009 to 2015. This may be due to the immediate reconstruction in disaster-affected areas after the earthquake and the stable state of reconstruction in subsequent years. This result is consistent with a previous study, which highlighted the positive effects of reconstruction on the stability of construction land after the earthquake [31]. Unlike its spatial global autocorrelation, the HH aggregation area of bare land continually decreased from 44.18% to 5.21% during the study period, indicating that natural disasters induced dispersed distribution of bare land, which were mainly distributed in the banks of the middle and lower stream of the Longxi River. In addition, the HH aggregation area of geological disasters increased intensively after the earthquake in 2009 but decreased with a higher percentage from 2011 to 2015, demonstrating a concentrated distribution of geological disasters in the north of the Longxi River Basin after natural recovery (Figure 18). This may be due to steep landforms and difficult human accessibility to the north of the basin that caused no treatment and frequent occurrence of geological disasters around the Longchi Lake distributed there [10,30].

4.4. The Boundedness and Prospect

The overall landscape pattern of the Longxi River Basin has undergone strong changes, resulting in fragmented landscapes. After 7 years, the damaged ecology has gradually recovered. The formation and evolution of landscape patterns are influenced by multiple factors. For example, from 2005 to 2009, due to the impact of the Wenchuan earthquake, the forest landscape area significantly decreased, and geological disasters and gully landscape area significantly increased, resulting in serious damage to the ecological environment, the decrease in wildlife habitats and resources, the reduction of water conservation capacity, and a high incidence of debris flow. In the future, closing mountains to facilitate afforestation; protecting natural forest, animal, and plant resources; and vigorously developing ecotourism will be important ways to reduce the impact of disasters and protect the environment of the basin. From 2009 to 2015, although the population in the watershed increased and geological disasters such as debris flows and landslides occurred frequently under heavy rainfall conditions, the ecological environment in the watershed gradually recovered. This is due to the policy of returning farmland to grassland and forest, as well as ecological and natural restoration. In addition, migration and resettlement, new town construction, and industrial structure adjustment also have a certain impact on human life and production activities in the basin. Affected by urbanization and human activities, the area of farmland and construction land in the Longxi River Basin has decreased, while the area of transportation land has increased.
In this study, we used four different remote-sensing satellite images to analyze the landscape pattern characteristics and spatial autocorrelation of the typical earthquake-affected basin. However, we only conducted the analysis using the 500 m × 500 m grid as a research cell. It is well known that the occurrence of spatial autocorrelation in geological elements can disappear or emerge going from one scale to the other [21]. Therefore, elucidating multi-scale effects on spatial autocorrelation of landscape would require further research. In addition, we were not able to consider potential driving factors of landscape pattern (e.g., elevation; slope gradient and aspect; distances to road, river, and fault; and rainfall). A spatial autoregressive model combining potential driving factors should be carried out to explain the correlation between landscape patterns and their driving factors. Human intervention, including reconstruction, disaster treatment, reforestation, and so on, can also change landscape patterns [60]. Understanding how human interventions combined with other driving factors control the process of landscape types could improve our knowledge of the ecological restoration of earthquake-affected areas.

5. Conclusions

We used the landscape index method to analyze the landscape pattern characteristics and the succession characteristics of landscape patterns in the Longxi River Basin. The following results are obtained:
(1)
In the study period, due to the earthquake, the degree of landscape fragmentation increased in 2009 compared with 2005 but decreased after 3 and 7 years of natural recovery.
(2)
Forest, farmland, and shrub-grassland in the basin had significantly positive spatial autocorrelation across the study period, while construction land only had significantly positive spatial autocorrelation after the earthquake. Bare land had significantly positive spatial autocorrelation except in 2009, while gully channel only had significantly positive spatial autocorrelation in 2005. Geological hazards shifted from significantly negative in 2005 to positive spatial autocorrelation from 2009 to 2015.
(3)
The Moran’s I of all landscape types decreased with the increase of space distance but in different distance-decay rates. In general, all landscape types had no significant spatial autocorrelation when the lag distance was >3 km.
(4)
Forest had the highest and fluctuating HH spatial aggregation area across the study period. Geological disasters had the clustered HH spatial aggregation areas upstream of the basin, while construction land was gradually distributed in Longchi town. The HH spatial aggregation areas of farmland distribution in the southeast of the basin increased in 2009 but then decreased. In addition, shrub-grassland and bare land had a continuously decreasing HH aggregation area.
Despite some limitations, our research contributes to improving landscape patterns, human intervention, and ecological restoration in earthquake-affected areas. Especially for post-disaster ecological restoration and reconstruction work, data support can be provided. For example, based on the results in this article, the degree of disaster in ecological areas can be graded and evaluated, and different ecological protection measures can be implemented for different regions, such as closing mountains to facilitate afforestation and returning farmland to forests. For post-disaster reconstruction, various reconstruction and relocation strategies can be implemented according to local conditions. Future research should extend this work by combining more influencing factors, and using more high-resolution remote-sensing satellite images to predict the spatial distribution of landscape types during ecological recovery.

Author Contributions

X.T., conception and design; fieldwork; data analysis and interpretation; writing the paper; submission preparation. L.Y., idea discussion; data analysis and interpretation; writing the paper. X.W., idea discussion; analysis and interpretation; writing the paper. J.W., analysis and interpretation; critical revision of the paper; idea discussion. Y.G. (Yiting Guo), paper revision; idea discussion. Y.G. (Yuhao Guo), fieldwork; data analysis. H.C., fieldwork; data analysis. J.L., fieldwork; data analysis. Y.L., idea discussion; fieldwork; writing the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Fujian Provincial Natural Science Foundation (2021J01060) and the National Natural Science Foundation of China (42071132).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to confidentiality requirements of the corresponding author’s university.

Acknowledgments

The authors thank C.H. Yu, and the staff of the Institute of Mountain Hazards and Environment, CAS for help in collecting and analyzing the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cui, P.; Lin, Y.M.; Chen, C. Destruction of vegetation due to geo-hazards and its environmental impacts in the Wenchuan earthquake areas. Ecol. Eng. 2012, 44, 61–69. [Google Scholar] [CrossRef]
  2. Allen, R.B.; Bellingham, P.J.; Wiser, S.K. Immediate damage by an earthquake to a temperate montane forest. J. Ecol. 1999, 80, 708–714. [Google Scholar] [CrossRef]
  3. Zhang, J.D.; Hull, V.; Xu, W.H.; Liu, J.G.; Ouyang, Z.Y.; Huang, J.Y.; Wang, X.Z.; Li, R.G. Impact of the 2008 Wenchuan earthquake on biodiversity and giant panda habitat in Wolong Nature Reserve, China. Ecol. Res. 2011, 26, 523–531. [Google Scholar] [CrossRef]
  4. Zhang, B.; Jiao, Q.; Wu, Y.; Zhang, W. Estimating soil erosion changes in the Wenchuan earthquake disaster area using geo-spatial information technology. J. Appl. Remote Sens. 2009, 3, 107–121. [Google Scholar] [CrossRef]
  5. He, J.; Tang, C.; Liu, G.; Li, W.L. Effect of landslides on the structural characteristics of land-cover based on complex networks. Int. J. Mod. Phys. B 2017, 31, 1750156. [Google Scholar] [CrossRef]
  6. Sheykhmousa, M.; Kerle, N.; Kuffer, M.; Ghaffarian, S. Post-disaster recovery assessment with machine learning-derived land cover and land use information. Remote Sens. 2019, 11, 1174. [Google Scholar] [CrossRef]
  7. Qiu, H.J.; Zhu, Y.R.; Zhou, W.Q.; Sun, H.S.; He, J.Y.; Liu, Z.J. Influence of DEM resolution on landslide simulation performance based on the Scoops3D model. Geomat. Nat. Hazards Risk 2022, 13, 1663–1681. [Google Scholar] [CrossRef]
  8. Tian, X.; Ma, X.; Huang, M.; Guo, Y.; Yang, H.; Yang, L.; Chen, H.; Gao, R.; Li, J.; Lin, Y. Spatiotemporal Dynamic Characteristics of Land Use in the Typical Watershed of Wenchuan Earthquake-Affected Areas—A Case Study in the Longxi River Basin. Sustainability 2022, 14, 15937. [Google Scholar] [CrossRef]
  9. Jelenek, J.; Kopačková-Strnadová, V. Synergic use of Sentinel-1 and Sentinel-2 data for automatic detection of earthquake-triggered landscape changes: A case study of the 2016 Kaikoura earthquake (Mw 7.8), New Zealand. Remote Sens. Environ. 2021, 265, 112634. [Google Scholar] [CrossRef]
  10. Li, C.; Wang, M.; Liu, K.; Coulthard, T.J. Landscape evolution of the Wenchuan earthquake-stricken area in response to future climate change. J. Hydrol. 2020, 590, 125244. [Google Scholar] [CrossRef]
  11. Duan, Y.; Di, B.; Ustin, S.L.; Chong, X.; Qiang, X.; Wu, S.; Li, J.; Zhang, R. Changes in ecosystem services in a montane landscape impacted by major earthquakes: A case study in Wenchuan earthquake-affected area, China. Ecol. Indic. 2021, 126, 107683. [Google Scholar] [CrossRef]
  12. Luo, X.Y.; Li, X.; Liu, W.; Wang, M.X.; Liu, C.H. Comprehensive evaluation of natural landscapes before and after earthquake based on GST-FAHP-GIS analytical method: A case study of Jiuzhaigou Nature Reserve, China. Ecol. Indic. 2023, 154, 110850. [Google Scholar] [CrossRef]
  13. Xiang, Z.L.; Dou, J.; Yunus, A.P.; Zhang, L.L.; Wang, X.K.; Luo, W.Q. Vegetation-landslide nexus and topographic changes post the 2004 Mw 6.6 Chuetsu earthquake. Catena 2023, 223, 106946. [Google Scholar] [CrossRef]
  14. Singh, R.P.; Bhoi, S.; Sahoo, A.K. Changes observed on land and ocean after Gujarat earthquake 26 January 2001 using IRS data. Int. J. Remote Sens. 2002, 23, 3123–3128. [Google Scholar] [CrossRef]
  15. Dey, S.; Singh, R.P. Surface latent heat flux as an earthquake precursor. Nat. Hazard Earth Syst. 2003, 3, 749–755. [Google Scholar] [CrossRef]
  16. Singh, V.P.; Singh, R.P. Changes in stress pattern around epicentral region of Bhuj earthquake of 26 January 2001. Geophys. Res. Lett. 2005, 32, 1–4. [Google Scholar] [CrossRef]
  17. Nath, B.; Niu, Z.; Singh, R.P. Land Use and Land Cover changes, and environment and risk evaluation of Dujiangyan city (SW China) using remote sensing and GIS techniques. Sustainability 2018, 10, 4631. [Google Scholar] [CrossRef]
  18. Pinheiro, E.; Camini, N.A.; Soares, M.; Sumida, S.S. Cellular automata model-landscape dynamics simulation tool in the process of change in land use and cover in the city of gacha do norte-Mt. In Proceedings of the 2020 IEEE Latin American GRSS & ISPRS Remote Sensing Conference (LAGIRS) 2020, Santiago, Chile, 22–26 March 2020; pp. 225–229. [Google Scholar] [CrossRef]
  19. Feng, Y.; Yang, Q.; Hong, Z.; Cuia, L. Modelling coastal land use change by incorporating spatial autocorrelation into cellular automata model. Geocarto. Int. 2018, 33, 470–488. [Google Scholar] [CrossRef]
  20. Liu, Y.; Feng, Y. Simulating the impact of economic and environmental strategies on future urban growth scenarios in Ningbo, China. Sustainability 2016, 8, 1045. [Google Scholar] [CrossRef]
  21. Overmars, K.; Koning, G.; Veldkamp, A. Spatial autocorrelation in multi-scale land use models. Ecol. Model. 2003, 164, 257–270. [Google Scholar] [CrossRef]
  22. Dormann, C.F. Effects of incorporating spatial autocorrelation into the analysis of species distribution data. Glob. Ecol. Biogeogr. 2007, 16, 129–138. [Google Scholar] [CrossRef]
  23. Caruso, T.; Migliorini, M.; Bucci, C.; Bargagli, R. Spatial patterns and autocorrelation in the response of microarthropods to soil pollutants: The example of oribatid mites in an abandoned mining and smelting area. Environ. Pollut. 2009, 157, 2939–2948. [Google Scholar] [CrossRef]
  24. Piotrowska, A.; Długosz, J. Spatio-temporal variability of microbial biomass content and activities related to some physicochemical properties of Luvisols. Geoderma 2012, 173, 199–208. [Google Scholar] [CrossRef]
  25. Sun, Y.Q.; Li, Y.H.; Li, Q.J.; Chen, Y.Z.; Liu, S.X.; Li, Y.L.; Li, J.H. Analysis of ecological risk evolution in the “Three Lakes” basin of Yuxi City based on landscape pattern. J. Hydroecol. 2022, 1–15. [Google Scholar] [CrossRef]
  26. Li, Y.J.; Shen, B.X.; Li, J.K.; Li, Y.; Zhou, X.; Liu, Y.W. Evolution and correlation analysis of habitat quality and landscape pattern in Dan-Hanjiang river basin. Environ. Sci. Technol. 2022, 45, 206–216. [Google Scholar] [CrossRef]
  27. Chang, Y.Y.; Gao, Y.; Xie, Z.; Zhang, T.Z.; Yu, X.Z. Spatiotemporal evolution and spatial correlation of habitat quality and landscape pattern over Beijing-Tianjin-Hebei region. China Environ. Sci. 2021, 41, 848–859. [Google Scholar]
  28. He, B.; Chang, J.X.; Guo, A.J.; Wang, Y.M.; Wang, Y.; Li, Z.H. Assessment of river basin habitat quality and its relationship with disturbance factors: A case study of the Tarim River Basin in Northwest China. J. Arid Land 2022, 14, 167–185. [Google Scholar] [CrossRef]
  29. Dong, J.H.; Zhang, Z.B.; Liu, B.T.; Zhang, X.H.; Zhang, W.B.; Chen, L. Spatiotemporal variations and driving factors of habitat quality in the loess hilly area of the Yellow River Basin: A case study of Lanzhou City, China. J. Arid Land 2022, 14, 637–652. [Google Scholar] [CrossRef]
  30. Zou, Q.; Cui, P.; He, J.; Lei, Y.; Li, S. Regional risk assessment of debris flows in China-An HRU-based approach. Geomorphology 2019, 340, 84–102. [Google Scholar] [CrossRef]
  31. Kang, D.; Yin, C.; Zhu, D.; Zou, S. Altitude and landslide scale regulated the assembly of grassland communities on landslides during the recovery process after the magnitude 8.0 Wenchuan earthquake, China. Ecol. Eng. 2021, 172, 106413. [Google Scholar] [CrossRef]
  32. Liu, G.Q.; Jiao, X. Ecological problems caused by Wenchuan earthquake and countermeasures. Soil Water Conserv. China 2008, 11, 11–13. [Google Scholar] [CrossRef]
  33. Hempel, J.; Michéli, E.; Owens, P.; McBratney, A. Universal soil classification system report from the International Union of Soil Sciences Working Group. Soil Horiz. 2013, 54, 1–6. [Google Scholar] [CrossRef]
  34. Ding, M.T.; Cheng, Z.L.; Wang, Q. Coupling mechanism of rural settlements and mountain disasters in the upper reaches of Min river. J. Mt. Sci. 2014, 11, 66–72. [Google Scholar] [CrossRef]
  35. Yang, Y.D.; Yan, X.S.; Zhang, H.B. Characteristics of SPOT5 remote sensing images of debris-flow hazards in Dongchuan County, Yunnan Province. J. Catastrophol. 2010, 25, 59–62+67. [Google Scholar] [CrossRef]
  36. Jia, L.; Du, J.K.; Zhao, P.; Zhang, Y.S. Land use dynamic monitoring in Hainan by TM data. Remote Sens. Inform. 2003, 22–25+56. [Google Scholar] [CrossRef]
  37. You, C.; Zhou, Y.B.; Yu, S.F. An Introduction of Quantitative Methods in Landscape Pattern Fragmentation. Chin. Agric. Sci. Bull. 2006, 22, 146–151. [Google Scholar] [CrossRef]
  38. Lei, J.R.; Chen, Z.Z.; Wu, T.T.; Li, Y.L.; Yang, Q.; Chen, X.H. Spatial autocorrelation pattern analysis of land use and the value of ecosystem services in northeast Hainan island. Acta Ecol. Sin. 2019, 39, 2366–2377. [Google Scholar] [CrossRef]
  39. Anselin, L.; Syabri, I.; Kho, Y. GeoDa: An introduction to spatial data analysis. Geogr. Anal. 2006, 38, 5–22. [Google Scholar] [CrossRef]
  40. Chen, J.Y.; Luo, H.X.; Hu, L.L.; Li, Q.H.; Jiang, W.J. Correlation between social economic development and ecosystem service value–A Case of Chongqing metropolitan area. Res. Soil Water Conserv. 2017, 24, 330–336. [Google Scholar] [CrossRef]
  41. Fu, J.X.; Zheng, F.L.; Li, Y.Y. Analysis of Land Use Spatial Autocorrelation Patterns and Influence Factors of Xiaolihe Watershed. China Acad. J. Electr. Pub. House 2017, 48, 128–138. [Google Scholar] [CrossRef]
  42. Legendre, P.L.; Legendre, L.F.J. Numerical Ecology. Ecology 1983, 63, 853. [Google Scholar] [CrossRef]
  43. Meng, B.; Wang, J.F.; Zhang, W.Z.; Liu, X.H. Evaluation of Regional Disparity in China Based on Spatial Analysis. Sci. Geogr. Sin. 2005, 25, 11–18. [Google Scholar]
  44. Wang, H.; Zhang, J.F.; Zhu, F.B.; Zhang, W.W. Analysis of spatial pattern of aerosol optical depth and affecting factors using spatial autocorrelation and spatial autoregressive model. Environ. Earth Sci. 2016, 75, 822–838. [Google Scholar]
  45. Anselin, L. Local indicators of spatial association-LISA. Geogr. Anal. 1995, 27, 93–115. [Google Scholar] [CrossRef]
  46. Ma, Y.J.; Xu, J.G. A Study on the Spatial-Temporal Development of Regional Economic Disparities in Henan Province Based on ESDA. Soft Sci. 2006, 20, 51–54. [Google Scholar] [CrossRef]
  47. Zhang, J.X.; Gong, J.; Ma, X.C.; Liu, D.Q. GeoDA-based spatial correlation analysis of landscape fragmentation in Bailongjiang Watershed of Gansu. Chin. J. Ecol. 2018, 37, 1476–1483. [Google Scholar]
  48. Dong, Y.H.; Peng, F.L.; Li, H.; Men, Y.Q. Spatial autocorrelation and spatial heterogeneity of underground parking space development in Chinese megacities based on multisource open data. Appl. Geogr. 2023, 153, 102897. [Google Scholar] [CrossRef]
  49. Gong, J.; Xie, Y.C. Spatiotemporal Changes of Watershed Landscape Pattern and Ecosystem Services: A Case Study of Bailongjiang Watershed in Gansu, 1st ed.; Science Press: Beijing, China, 2018; pp. 143–146. [Google Scholar]
  50. Wang, L.; Li, T.; Cao, X.S. Spatial Difference of Urban and Rural Harmonious Development with ESDA-GIS Framework in Guangdong Province. Econ. Geogr. 2012, 32, 44–50. [Google Scholar]
  51. Chen, Y.S. Adaptive Study on Local Spatial Autocorrelation Index. Master’s Thesis, East China Normal University, Shanghai, China, 1 June 2009. [Google Scholar]
  52. Aamir, M.; Ali, T.; Irfan, M.; Shaf, A.; Azam, M.Z.; Glowacz, A.; Brumercik, F.; Glowacz, W.; Alqhtani, S.; Rahman, S. Natural disasters intensity analysis and classification based on multispectral images using multi-layered deep convolutional neural network. Sensors 2021, 21, 2648. [Google Scholar] [CrossRef]
  53. Bian, Z.; Zhang, D.; Xu, J.; Tang, H.; Bai, Z.; Li, Y. Study on the evolution law of surface landscape pattern in earthquake-stricken areas by remote sensing: A case study of Jiuzhaigou County. Sichuan Province. Sustainability 2022, 14, 13032. [Google Scholar] [CrossRef]
  54. Ren, Y.; Lv, Y.; Fu, B.J.; Zhang, K. Biodiversity and ecosystem functional enhancement by forest restoration: A meta-analysis in China. Land Degrad. Dev. 2017, 28, 2062–2073. [Google Scholar] [CrossRef]
  55. Li, J.; He, J.; Liu, Y.; Wang, D.; Rafay, L.; Chen, C.; Hong, T.; Fan, H.; Lin, Y. Spatial Autocorrelation Analysis of Multi-Scale Damaged Vegetation in the Wenchuan Earthquake-Affected Area, Southwest China. Forests 2019, 10, 195. [Google Scholar] [CrossRef]
  56. Shi, Y.F.; Di, B.F.; Zuo, Q.; Wu, S.L.; Duan, Y.N. Interactive relationship between industrial structure and land use pattern in the Wenchuan earthquake stricken area, China. Mt. Res. 2020, 38, 916–925. [Google Scholar] [CrossRef]
  57. Zou, Y.D.; He, L.; Zhang, X.P.; Ma, B.Y.; Wang, H.J.; Wang, M.Q.; Xue, F.; He, J. Characteristics of Land Use Structure Change in Beiluo River Basin During 1970–2019 Based on Google Erath Engine. Bull. Soil Water Conserv. 2021, 41, 209–219. [Google Scholar] [CrossRef]
  58. Li, Z.; He, Z.W.; Xu, H.X. The research of land use dynamic change before and after the earthquake in Dujiangyan City. J. South Chin. Norm. Univ. Nat. Sci. Ed. 2012, 35, 412–417. [Google Scholar] [CrossRef]
  59. Wang, X.; Mao, H. Spatio-temporal evolution of post-seismic landslides and debris flows: 2017 Ms 7.0 Jiuzhaigou earthquake. Environ. Sci. Pollut. Res. 2022, 29, 15681–15702. [Google Scholar] [CrossRef]
  60. Piquer-Rodríguez, M.; Butsic, V.; Gärtner, P.; Macchi, L.; Baumann, M.; Pizarro, G.G.; Volante, J.N.; Gasparri, I.N.; Kuemmerle, T. Drivers of agricultural land-use change in the Argentine Pampas and Chaco regions. Appl. Geogr. 2018, 91, 111–122. [Google Scholar] [CrossRef]
Figure 1. Location of Longxi River Basin.
Figure 1. Location of Longxi River Basin.
Forests 14 02349 g001
Figure 2. Distribution map of main channels in Longxi River Basin.
Figure 2. Distribution map of main channels in Longxi River Basin.
Forests 14 02349 g002
Figure 3. Spatial distribution of edge density (ED) of Longxi River Watershed from 2005 to 2015.
Figure 3. Spatial distribution of edge density (ED) of Longxi River Watershed from 2005 to 2015.
Forests 14 02349 g003
Figure 4. Spatial distribution of Shannon diversity (SHDI) of Longxi River Watershed from 2005 to 2015.
Figure 4. Spatial distribution of Shannon diversity (SHDI) of Longxi River Watershed from 2005 to 2015.
Forests 14 02349 g004
Figure 5. Spatial distribution of landscape division (DIVISION) of Longxi River Watershed from 2005 to 2015.
Figure 5. Spatial distribution of landscape division (DIVISION) of Longxi River Watershed from 2005 to 2015.
Forests 14 02349 g005
Figure 6. The correlograms of the global Moran’s I of all landscape types with a lag distance of 10 km.
Figure 6. The correlograms of the global Moran’s I of all landscape types with a lag distance of 10 km.
Forests 14 02349 g006
Figure 7. Moran scatter plots of various landscapes in the Longxi River Watershed in 2005.
Figure 7. Moran scatter plots of various landscapes in the Longxi River Watershed in 2005.
Forests 14 02349 g007
Figure 8. Moran scatter plots of various landscapes in the Longxi River Watershed in 2009.
Figure 8. Moran scatter plots of various landscapes in the Longxi River Watershed in 2009.
Forests 14 02349 g008
Figure 9. Moran scatter plots of various landscapes in the Longxi River Watershed in 2011.
Figure 9. Moran scatter plots of various landscapes in the Longxi River Watershed in 2011.
Forests 14 02349 g009
Figure 10. Moran scatter plots of various landscapes in the Longxi River Watershed in 2015.
Figure 10. Moran scatter plots of various landscapes in the Longxi River Watershed in 2015.
Forests 14 02349 g010
Figure 11. LISA maps of forest land for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Figure 11. LISA maps of forest land for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Forests 14 02349 g011
Figure 12. LISA maps of farmland for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Figure 12. LISA maps of farmland for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Forests 14 02349 g012
Figure 13. LISA maps of shrub-grassland for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Figure 13. LISA maps of shrub-grassland for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Forests 14 02349 g013
Figure 14. LISA maps of construction land for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Figure 14. LISA maps of construction land for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Forests 14 02349 g014
Figure 15. LISA maps of traffic land for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Figure 15. LISA maps of traffic land for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Forests 14 02349 g015
Figure 16. LISA maps of bare land for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Figure 16. LISA maps of bare land for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Forests 14 02349 g016
Figure 17. LISA maps of the gully channel for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Figure 17. LISA maps of the gully channel for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Forests 14 02349 g017
Figure 18. LISA maps of geological hazard for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Figure 18. LISA maps of geological hazard for Longxi River Watershed from 2005 to 2015. (HH: High-high cluster; LL: Low-low cluster; HL: High-low anomaly; LH: Low-high anomaly; NS: Not significant).
Forests 14 02349 g018
Table 1. Satellite imagery data information in the study area.
Table 1. Satellite imagery data information in the study area.
Satellite ImageryDateIntervalsResolution/mCoverageData Quality
QuickBird26 June 2005Before the earthquake0.61Entire study areaFew Clouds
SPOT-510 February 20091 year after the earthquake2.50Entire study area Cloudless
QuickBird26 April 20013 years after the earthquake0.61Entire study areaFew Clouds
Worldview-215 April 20157 years after the earthquake0.46Entire study areaCloudless
Table 2. Changes in landscape area from 2005 to 2015.
Table 2. Changes in landscape area from 2005 to 2015.
Year ForestFarmlandShrub-GrasslandConstructionTraffic LandBare LandGully ChannelGeological Hazards
Land
2005area/hm27253.63136.43336.9841.0122.6262.1751.306.18
percentage/%91.701.724.260.520.290.790.650.08
2009area/hm26213.5585.11269.9436.0038.0575.18175.941016.76
percentage/%78.551.083.410.460.480.952.2212.85
2011area/hm26864.1066.15320.5629.9541.0644.04175.78377.42
percentage/%86.680.844.050.380.520.562.224.77
2015area/hm27186.1349.67220.7030.5748.6437.51165.97179.25
percentage/%90.750.632.790.390.610.472.102.26
Table 3. Global spatial autocorrelation characteristics of various landscapes in Longxi River Watershed from 2005 to 2015 at 0.5 km spatial aggregation distance. z-score is the critical value of the test statistic; p-value indicates the significance level; |z| > 1.96, p < 0.05, i.e., confidence level > 95%; |z| > 2.56, p < 0.01, i.e., confidence level > 99%.
Table 3. Global spatial autocorrelation characteristics of various landscapes in Longxi River Watershed from 2005 to 2015 at 0.5 km spatial aggregation distance. z-score is the critical value of the test statistic; p-value indicates the significance level; |z| > 1.96, p < 0.05, i.e., confidence level > 95%; |z| > 2.56, p < 0.01, i.e., confidence level > 99%.
YearParametersForestFarmlandShrub-GrasslandConstruction
Land
Traffic LandBare LandGully ChannelGeological Hazards
2005Moran’s I0.2340.3000.2600.085−0.0210.3750.535−0.748
z5.1482.5923.5100.879−0.0333.1085.934−7.397
p-value<0.001<0.01<0.0010.3800.974<0.01<0.001<0.001
2009Moran’s I0.4780.4940.4740.412−0.1030.0870.1870.482
z9.8493.2947.4544.355−0.8621.0142.3258.680
p-value<0.001<0.001<0.001<0.0010.3890.3100.020<0.001
2011Moran’s I0.2290.5160.2710.3270.0880.2670.1320.458
z4.8653.9824.3993.0900.9982.9791.6327.408
p-value<0.001p < 0.001<0.001<0.010.318<0.010.103<0.001
2015Moran’s I0.2430.3650.3640.3760.0360.3520.1820.531
z5.1672.7315.9133.4900.4652.7772.2827.889
p-value<0.001<0.01<0.001<0.0010.642<0.010.023<0.001
Table 4. Statistics of spatial clustering and abnormal area of various landscapes in Longxi River Watershed.
Table 4. Statistics of spatial clustering and abnormal area of various landscapes in Longxi River Watershed.
ParametersArea/km2 (Proportion/%)
ForestFarmlandShrub-GrasslandConstruction
Land
Traffic LandBare LandGully ChannelGeological Hazards
2005
HH15.64(21.56)0.29(21.55)0.74(22.10)0.05(11.05)0.04(15.72)0.27(44.18)0.13(24.85)-
LL0.89(1.22)0.05(3.50)0.09(2.72)0.001(0.31)0.001(0.38)0.01(2.04)0.01(2.71)-
HL2.51(3.46)-0.08(2.35)0.005(1.11)0.02(9.71)---
LH-0.05(3.55)0.08(2.49)0.01(2.40)-0.01(0.87)0.003(0.60)-
NS53.50(73.76)0.97(71.41)2.37(70.34)0.35(85.13)0.17(74.20)0.33(52.91)0.37(71.85)0.06(100)
2009
HH7.04(11.33)0.32(37.69)0.72(26.69)0.15(40.92)-0.08(10.18)0.37(20.82)2.29(22.57)
LL0.49(0.79)0.05(6.33)0.05(1.77)0.01(2.42)-0.02(2.02)0.03(1.93)0.29(2.82)
HL2.07(3.33)0.03(3.42)0.01(0.43)-0.02(6.40)0.03(3.45)--
LH0.35(0.57)0.02(2.31)0.03(1.21)0.02(4.56)0.01(2.21)0.01(1.19)0.02(1.40)0.16(1.61)
NS52.19(83.99)0.43(50.26)1.89(69.91)0.19(52.10)0.35(91.40)0.62(83.16)1.33(75.84)7.42(73.00)
2011
HH10.76(15.67)0.22(33.74)0.54(16.75)0.13(42.93)0.01(2.75)0.06(13.71)0.24(13.86)1.17(31.02)
LL1.00(1.45)0.03(4.59)0.06(1.91)0.002(0.80)0.005(1.18)0.001(0.28)0.03(1.45)0.09(2.32)
HL2.14(3.11)-0.04(1.12)-0.02(5.99)--0.03(0.79)
LH-0.01(1.37)0.10(3.19)0.004(1.43)0.01(3.19)0.004(0.90)0.01(0.66)0.03(0.81)
NS54.76(79.77)0.40(60.30)2.47(77.03)0.16(54.84)0.36(86.88)0.37(85.11)1.48(84.03)2.45(65.05)
2015
HH12.71(17.69)0.04(8.04)0.52(23.61)0.15(47.74)-0.02(5.21)0.34(20.49)0.65(36.11)
LL0.81(1.12)0.03(6.61)0.04(1.82)0.004(1.33)0.01(1.77)0.002(0.48)0.03(1.70)0.04(2.17)
HL2.28(3.17)0.02(3.16)--0.02(4.78)0.02(4.88)-0.04(2.47)
LH-0.02(3.19)0.05(2.20)0.004(1.30)0.01(2.88)0.003(0.85)0.03(1.64)0.03(1.58)
NS56.07(78.02)0.39(79.00)1.60(72.37)0.15(49.63)0.44(90.57)0.33(88.58)1.26(76.17)1.03(57.67)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tian, X.; Yang, L.; Wu, X.; Wu, J.; Guo, Y.; Guo, Y.; Chen, H.; Li, J.; Lin, Y. Landscape Fragmentation and Spatial Autocorrelation of a Typical Watershed in the Wenchuan Earthquake-Affected Area—A Case Study in the Longxi River Basin. Forests 2023, 14, 2349. https://doi.org/10.3390/f14122349

AMA Style

Tian X, Yang L, Wu X, Wu J, Guo Y, Guo Y, Chen H, Li J, Lin Y. Landscape Fragmentation and Spatial Autocorrelation of a Typical Watershed in the Wenchuan Earthquake-Affected Area—A Case Study in the Longxi River Basin. Forests. 2023; 14(12):2349. https://doi.org/10.3390/f14122349

Chicago/Turabian Style

Tian, Xue, Liusheng Yang, Xuan Wu, Jianzhao Wu, Yiting Guo, Yuhao Guo, Hui Chen, Jian Li, and Yongming Lin. 2023. "Landscape Fragmentation and Spatial Autocorrelation of a Typical Watershed in the Wenchuan Earthquake-Affected Area—A Case Study in the Longxi River Basin" Forests 14, no. 12: 2349. https://doi.org/10.3390/f14122349

APA Style

Tian, X., Yang, L., Wu, X., Wu, J., Guo, Y., Guo, Y., Chen, H., Li, J., & Lin, Y. (2023). Landscape Fragmentation and Spatial Autocorrelation of a Typical Watershed in the Wenchuan Earthquake-Affected Area—A Case Study in the Longxi River Basin. Forests, 14(12), 2349. https://doi.org/10.3390/f14122349

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