The 2008 Wenchuan earthquake destroyed large areas of vegetation. Presently, these areas of damaged vegetation are at various
stages of recovery. In this study, we present a probabilistic approach for slope stability analysis that quantitatively relates
data on earthquake-damaged vegetation with slope stability in a given river basin. The Mianyuan River basin was selected for model
development, and earthquake-damaged vegetation and post-earthquake recovery conditions were identified via the normalized
difference vegetation index (NDVI), from multi-temporal (2001–2014) remote sensing images. DSAL (digital elevation model, slope,
aspect, and lithology) spatial zonation was applied to characterize the survival environments of vegetation, which were used to
discern the relationships between successful vegetation regrowth and environmental conditions. Finally, the slope stability
susceptibility model was trained through multivariate analysis of earthquake-damaged vegetation and its controlling factors
(i.e. topographic environments and material properties). Application to the Subao River basin validated the proposed model, showing
that most of the damaged vegetation areas have high susceptibility levels (88.1 %
On 12 May 2008, the high-magnitude (Ms 8.0) Wenchuan earthquake struck China, and the large release of energy triggered many geo-hazards at different scales; catastrophic landslides and collapses were particularly widespread in middle- and high-elevation mountainous areas (Cui et al., 2009; Huang and Li, 2009a). These geo-hazards resulted in massive movements of surface material that led to large areas of damaged vegetation, and in some places the vegetation was completely destroyed (Cui et al., 2012). Presently, the areas of damaged vegetation are at various stages of recovery and some have fully recovered within only a few years of the earthquake (Liu et al., 2010). Nevertheless, the potential for future slope instability still exists in this region (Khan et al., 2013). Few studies have focused explicitly on the recovery process of earthquake-damaged vegetation in this region. For example, Zhang et al. (2014) documented the natural recovery of forests after the earthquake and revealed that factors including soil cover and slope were correlated with successful vegetation recovery (Zhang et al., 2014). Soil moisture is one of the most important environmental factors for vegetation recovery in the landslide sites (Lin et al., 2005).
The occurrence and frequency of landslides in an area depend fundamentally on the interaction between triggering mechanisms and natural conditions (Lee et al., 2008; Guzzetti et al., 2012; Peruccacci et al., 2012; Borgomeo et al., 2014). Events such as the Wenchuan earthquake not only trigger serious co-seismic landslides, but can also lead to increased long-term post-seismic slope instability (Koi et al., 2008; Tang et al., 2011). We consider that the locations of vegetation damage are related to slope stability, and that recovery conditions can represent its suffered landslides. Establishing statistical relationships between damaged vegetation and the environmental conditions under which vegetation will recover, hereafter called the “survival environment,” is an important task for predicting and mapping areas susceptible to future slope failure.
This paper presents a general framework for determining the relationships between landslide-damaged vegetation and the topographical and vegetative factors associated with survival environments. The Mianyuan River basin was selected for model development, and we identified earthquake-damaged vegetation and post-earthquake recovery conditions using the normalized difference vegetation index (NDVI) from multi-temporal (2001–2014) remote sensing images. We then applied DSAL (digital elevation model, slope, aspect, and lithology) spatial zonation to characterize the survival environments of vegetation, which were used to discern the relationships between environmental conditions and successful vegetation regrowth. Finally, the slope susceptibility model was trained through multivariate analysis of earthquake-damaged vegetation and its controlling factors (i.e. topographic environments and material properties). The proposed model was then applied to the Subao River basin to evaluate whether or not the data could be used to assess future susceptibility to slope failure in regions affected by the Wenchuan earthquake. In addition, the recovery capacity model of the damaged vegetation was trained through multivariate analysis of recovery vegetation and its survival factors (i.e. topographic environments and material properties), and was used to assess the landslide processes of the damaged vegetation areas under certain conditions.
The study area is located in the upstream reaches of the Mianyuan River watershed (31
The Mianyuan River is situated in the transitional mountainous belt between the Sichuan basin and the western Sichuan Plateau,
which is characterized by rugged mountains with deeply incised valleys. Active geotectonic movements induced by the complicated
fault system occur frequently in this region. The Longmenshan thrust belt, which ruptured during the 2008 Wenchuan earthquake,
runs through the central part of the Mianyuan River basin (Fig. 1c). The rocks of Mianyuan River basin primarily consist of Sinian
sandstone and siltstone (Z); Cambrian sandstone, siltstone, and slate (E); Silurian phyllite, schist, and slate (S); Devonian
dolomite limestone and sandstone (D); Carboniferous limestone (C); Permian limestone and shale (P); Triassic sandstone, dolomite,
limestone, siltstone, and shale (T); Quaternary deposits (Q); and magmatic rocks, granite, and diorite (r). These lithologies are
based on
The basin is influenced by a subtropical moist climate and monsoonal rains that start in early June and continue until September.
Mean annual precipitation is approximately 1500–1700
During the Wenchuan earthquake, the earthquake-induced geo-hazards resulted in massive movement of surface material that
diminished and destroyed large areas of vegetation (Fig. 2c and d). These damaged vegetation areas alter the spectral signatures
and NDVI values recorded on remote sensing images. These data sets are derived from both the red and near-infrared spectral bands
and are sensitive to changes in biophysical conditions of vegetation, and can therefore be used to detect these damaged
areas. Following the earthquake, obvious changes in NDVI values indicated areas (Fig. 2e and f) of damaged vegetation (Liu et al.,
2012), which were distributed along the stream network of the basin (Fig. 2g). Statistical analysis shows that the total size of
the damaged area was about 95.4
More recently, 17 June 2014 Landsat-8 images (Fig. 2h) of the basin indicate that most areas of damaged vegetation have
recovered, whereas other areas have not. To analyse this recovery process in detail, the 17 June 2014 NDVI image was processed to
extract the unrecovered vegetation areas. The earthquake-damaged and recovered vegetation areas are shown in Fig. 2i. The
results indicate that these areas of damaged vegetation exhibit various stages of recovery depending on their environments
(particularly local topographical conditions). Statistical analysis shows that approximately 57.5
Vegetation growth primarily depends on sunshine, water, temperature and nutrients, and these inputs are directly related to environmental conditions. For example, topography influences vegetation growth in a variety of ways: elevation influences temperature, aspect determines sunshine, and slope affects hydrological conditions. In addition, the geological setting (especially lithology) controls soil properties. River basins represent natural hydrological units for which it is possible to determine balances between the major constituent fluxes of rainfall, evaporation, and river discharge, along with groundwater storage (Bathurst et al., 2010). Additionally, river basins are a natural result of changes in geomorphic processes. Hence, river basins can be characterized by their natural features, particularly in terms of topographical environments and vegetation growth.
On the basis of these natural features, the DEM and slope gradient data were mapped into four classes (1:
According to the above considerations, DSAL spatial zonation was defined as one spatial zonation for one given river basin, based on the natural features of the DEM, slope data, aspect data, and lithology conditions. In theory, one river basin can be mapped into 256 DSAL classes (Zhang et al., 2015).
To determine the key topographic characteristics of Mianyuan River basin, digital elevation model (DEM) data (resolution:
25
To determine the relationships between damaged vegetation and survival environments, the vegetation damage probability
(
All of the above results indicate that the probabilities of vegetation damage and its recovery capacity are both strongly correlated with topographical characteristics. Additionally, these results demonstrate that DSAL spatial zonation is reasonable for detecting the relationships between earthquake-damaged vegetation and survival environments.
Vegetation cover, especially of a woody type with strong and large root systems, helps to improve the stability of slopes (Gray, 1996; Wu and Sidle, 1995; Cammeraat et al., 2005; Greenwood et al., 2004). Importantly, vegetation helps to stabilize slope materials by two main mechanisms: (1) above-ground biomass changes the soil hydrology extensively through evaporation and transpiration processes (Marston, 2010; Haneberg, 1991; Harden, 2006), and leaves and litter intercept raindrops, thereby dissipating erosive energy (Parsons et al., 1996; Marston and Dolan, 1999; Keim and Skaugset, 2003), (2) the development of root systems affects the mechanical and hydrological properties of soils. Thus, vegetation improves the resistance of slopes to both surficial erosion and mass wasting. Conversely, the removal of slope vegetation tends to accelerate or increase slope failures (Gray, 1996). Although the impacts of vegetation on slope stability in mountainous regions are reasonably well understood and documented, it is still difficult to predict exactly how vegetation will impact mass movement processes including debris slides, mudflows, and rockfalls.
Endo and Tsuruta (1968) determined the reinforcement effect of a root system based on shear strength, and showed that increased strength is directly proportional to root density (Endo and Tsuruta, 1968). Roots are the primary pathway for water and nutrient uptake by plants, and the surface material along the slope creates a root ecosystem conducive to vegetation growth. The nutrient and water inputs, which control vegetation growth, are positively correlated with surface material properties. Meanwhile, surface materials make up the bulk of landslides, and landslide processes are largely determined by the slope gradient, slope materials, elevation, and hydrological conditions. All these causative factors influence vegetation growth and the vegetation survival environment. Thus, there are some coupling relationships between landslide processes and vegetation growth, and spatial analysis of landslide-damaged vegetation is important for understanding which terrain and surface materials are susceptible to landslide processes. Although there is not yet a complete understanding of the interaction of vegetation and landslide processes, these coupling relationships can be confirmed by the use of statistical analyses or other methods.
The Wenchuan earthquake triggered many geo-hazards, especially landslides and collapses in the middle- and high-mountainous areas of China. The reduction of shear strength and overall soil cohesion, combined with ground acceleration from seismic waves, was responsible for the failure of many shallow slopes in the affected region (Yang and Chen, 2010; Zhang et al., 2012). In essence, the earthquake-induced geo-hazards are strongly controlled by the three most important factors, i.e. topographical conditions, slope material properties, and ground shaking. For it is one quite complex process that an earthquake produces a range of ground shaking levels at sites throughout the region, the ground shaking condition is more difficult to handle in modelling practice. We consider that earthquake-induced geo-hazards are much influenced by topographical conditions and slope materials properties at the basin scale.
Previous studies in the area affected by the Wenchuan earthquake have demonstrated that earthquake-induced geo-hazards are
strongly correlated with elevation, slope gradient and aspect, geological units, distances to the epicenter and active faults, and
seismic intensity (Huang and Li, 2009b; Qi et al., 2010; Cui et al., 2009; Dai et al., 2011; Xu et al., 2013a, b; Chen et al.,
2012a, c; Zhang et al., 2011c). Our spatial distribution analysis of earthquake-damaged vegetation indicates that the probability
of vegetation damage
The model results show that high-susceptibility areas were located along the steep slopes of main streams. Predicted error was defined as the difference between the mean of the values predicted by the model and the actual values in each DSAL region. The results (Fig. 4b) show good agreement between predicted susceptibility and actual conditions, except for mountaintop areas and some deep-seated landslide areas.
The above results indicate that the quantitative slope stability susceptibility should be regarded as providing a general assessment of landslide susceptibility rather than the precise probability of future landslides, owing to the complexity of the landslide processes. Actually, the predicted slope stability susceptibility at a given basin is relative, and it is no mean to compare two predicted values at different basins. In practice, it is more useful to divide the predicted slope stability susceptibility at one given basin into different levels, identifying those areas with high susceptibility to future landslides.
In addition, the recovery capacity of damaged vegetation is strongly influenced by its survival environments (especially soil moisture) (Lin et al., 2005) and the landslide processes that occurred. In the Mianyuan River basin, these unrecovered areas had suffered more serious landslide processes than the recovered areas in similar survival environments, such as the deep-seated landslides that occurred in Wenjia gully. Hence, monitoring the recovery conditions of landslide-damaged vegetation can provide useful information about the landslide processes that occurred.
The Subao River basin (31
Furthermore, Subao River basin is under the influence of a humid subtropical monsoon climate, which is warm and wet. The mean
annual temperature is 15.6
The areas of damaged vegetation in Subao River basin were used to verify the predicted slope stability susceptibility map. To assess the accuracy of the proposed method, the predicted slope stability susceptibilities were divided into five levels and mapped (Fig. 6); the results were then compared to known landslide-prone areas, by means of high-precision imagery and field survey data (Fig. 6). Areas predicted to have high susceptibility overlapped with many of the landslide-prone areas. Statistical analysis indicated that 88.1 % of the total damaged vegetation areas have high susceptibility levels (exceeding level 3), and that 61.5 % of the total damaged vegetation areas in the Subao River basin exceeded susceptibility level 4.
In addition, the recovery conditions of damaged vegetation areas and the predicted recovery probability can be used to analyze the landslide processes that occurred during earthquake events. Statistical analysis shows that 55.4 % of the total unrecovered area has high recovery capacity (exceeding 0.60) and 11.4 % of the total recovered area has low recovery capacity (less than 0.45). These unrecovered areas with high recovery capacity had suffered more serious landslide processes than the recovered areas with low recovery capacity.
We know that one key factor for earthquake-induced landslide is the grounding shaking experienced at that location during the
earthquake. For example, Minayuan River is mainly located in the seismic intensity
The results of the application indicate that vegetative features and data on damaged vegetation can be used to assess future landslide susceptibilities in areas recovering from an earthquake. In addition, the recovery conditions can represent the landslide processes that occurred in some conditions, and these detailed relationships require further analysis and research. Although our modeling approach would benefit from further development, the data collected to date demonstrate that this approach can be valuable for disaster mitigation efforts.
Our study presents a probabilistic approach for slope stability analysis that quantitatively relates earthquake-damaged vegetation data with slope stability susceptibilities in a given river basin. The results indicate that spatial analysis of earthquake-damaged vegetation can provide useful information about the types of terrain and surface materials that are susceptible to landslide processes, and that historical earthquake-damaged vegetation data can be used for predicting future slope instability.
Presently, the areas of damaged vegetation are in various stages of recovery and a few areas of the Wenchuan earthquake-hit area have fully recovered. In essence, these recovered areas returned to geomorphic equilibrium within a few years after the earthquake, possibly because most of the landslides that damaged vegetated areas were of shallow type, and because the high recovery capacity of the vegetation helped stabilize the slopes. Nevertheless, there remains potential for future slope instability in the region. Hence, monitoring the recovery processes of earthquake-damaged vegetated areas is helpful for understanding landslide processes and predicting future landslide hazards.
This research was financially supported by the National Science and Technology Support Program (No: 2012BAC06B02), the key Projects of the Chinese Academy of Sciences (No: KZZD-EW-05-01-04), and the Beijing Natural Science Foundation (No: 4144088).
Parameters of the function (Eq. 1) in the Mianyuan River basins.
Parameters of the function (Eq. 3) in the Mianyuan River basins.
The study area:
2001–2014 vegetation changes in Mianyuan River basin:
Predicted slope susceptibility in Mianyuan River basin
Subao River basin:
Predicted slope stability