Abstract

Volume calculation is important for quantifying the erosion driven by coseismic landslides in geomorphology. With the advent of digital elevation models (DEMs), quantifying features of landslide bodies have become possible, permitting to calculate the landslide volume in terms of elevation changes. To further test this approach, this work calculates the volume of landslides near the epicenter of the 2008 Mw 7.9 Wenchuan earthquake in the Yingxiu region, Sichuan, China, by comparing pre- and postearthquake DEMs. Results suggest that effective application of this method needs to consider the DEM resolution and eliminates background errors of individual landslides. The volume of coseismic landslides calculated by the proposed method may represent a minimum value compared to that from the existing empirical V-A formulas. Considering that it is difficult to quantify the coseismic landslide volume throughout in a broad region, this method can be applied to the preliminary stage of characterizing coseismic landslides quantitatively for some key localities of the affected area of major earthquakes.

1. Introduction

Frequency, area, and volume are three main parameters to measure seismically induced landslides. Determination of these parameters is important for geologic hazard assessment and examining the role of landslides in landscape evolution [16]. For example, by comparing the estimated coseismic landslide volume with the volume of material added to the orogen through coseismic rock uplift, Parker et al. argued that mass wasting triggered by the 2008 Wenchuan earthquake was greater than orogenic growth, leading to a net material deficit in the Longmenshan [6]. In general, the total number and area of landslides can be estimated based on aerial photography or remote sensing data, coupled with field investigations [711]. In contrast, estimating the landslide volume remains difficult because the depth of landslide bodies is typically unknown [1012].

To solve this problem, researchers use a power-scaling relation between the volume and area to estimate the volume of individual landslides or total volume of multiple landslides. Such a relation is expressed aswhere V is the volume of landslides, A is the area of landslides, and α and γ are coefficients. This relation was originally proposed by Simonett [12] which assumes that the volume of a landslide is positively correlated with its area, determined using regression analysis of morphometric data. This formula has been applied to many case studies [5, 6, 13]. However, variations in landslide distribution and size, which are influenced by environmental factors (e.g., tectonic setting, lithology, grain size, and slope), can create uncertainties in the result of this approach, especially for large-scale landslides [10, 1317]. Sometimes, the landslide volume derived from the area-volume equation is much larger than that estimated based on the relation between the landslide volume and moment magnitude of the earthquake [1, 3, 18].

Recently, high-resolution and multitemporal DEMs have proven valuable to identify patterns of landscape evolution with coseismic topographic changes [1820]. It is possible to estimate the landslide volume in terms of elevation changes from DEMs before and after earthquakes. For example, in their study on the landslides triggered by the 2008 Wenchuan M7.9 earthquake, Xu et al. (2016) calculated volumes of 1,415 landslides based on differences between pre- and postearthquake DEMs with 20 m resolution. They proposed three optimized models to calculate volumes of these coseismic landslides, depending on multiple morphometric attributes (area, length, width, height etc.). Unfortunately, uncertainties remain a problem when calculating differences between DEMs, as they often are collected from various data sources [20]. Both pre- and postseismic DEMs contain elevation errors, so comparison of continuous GPS data is required to yield relatively accurate elevation data. Besides, the elevation error is not a constant, instead varying with locations [18], which makes it somewhat difficult using DEM data to measure the volume of landslides.

Here, we propose a method to more accurately calculate the volume of coseismic landslides by removing measurement error caused by differences in DEM resolution and other irregularities using pre- and postearthquake DEMs. The test of this approach is performed based on the data from the Yingxiu area, Sichuan, China, near the epicenter of the 2008 Mw 7.9 Wenchuan earthquake which triggered massive coseismic slope failures. We then compare calculated coseismic landslide volumes with the values from different area-volume relations. Finally, we discuss the advantages of this method and provide a linear relation between log(V) and log(A).

2. Data and Methods

2.1. Study Area

The 2008 Mw 7.9 Wenchuan earthquake triggered thousands of coseismic landslides over an expansive region [21]. The study area of this work is located near the epicenter of this earthquake covering 66 km2. This region has the highest density of coseismic landslides caused by the event [7, 21] and experienced a Modified Mercalli Intensity (MMI) of XI [22]. The exposed strata in the study area are mainly sourced from the Precambrian Pengguan massif [23]. Topographically, the area is featured by an approximately 2,350 meter deep canyon, where the Minjiang River and its two tributaries converge. The Wenchuan earthquake surface rupture zone passes through the southeastern border of the study area near the epicenter, and the majority of landslides resides within the hanging wall of the seismogenic fault (Figures 1(a) and 1(b)).

2.2. Data

Two DEMs from different sources are used in this study. The pre-earthquake DEM is downloaded from 30 m Shuttle Radar Topography Mission (SRTM) elevation data (http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp). The postearthquake DEM is derived from GF2 remote images after radiometric and orthophoto correction and panchromatic and multispectral fusion processing and has a spatial resolution of 10 meters.

Based on aerial photography (with a resolution of 2 m) in June 2008, we manually interpret coseismic landslides in the study area following the criteria suggested by Dai et al., namely, (a) landslide deposits that can be clearly identified at the foot of the slope covering roads or other infrastructures or extending into rivers or streams; (b) landslides that show newly denuded vegetation on the slope; (c) landslides that show bright white or dark brown contrast to surrounding slopes; and (d) landslide debris movement paths that can be clearly observed [7]. Because it is difficult to make a distinction between the source area and deposit for most landslides on aerial photos, each mapped polygon includes both source and deposit areas. Also, some small landslides are considered as single landslides when they appear as partly merged. As a result, more than 304 landslides are mapped with areas ranging 1800 to over 160,000 m2. Although the study area is a very small part of the Wenchuan earthquake affected area, the relationship between the landslide area and frequency density, which is defined as the ratio between the landslide number within the corresponding landslide area class and the total number of landslides in the whole inventory, is consistent with that of Malamud et al. [3] and others [11, 24, 25] (Figure 1(c)).

2.3. Methods

Landsliding can change elevation of the affected area (Figure 2). In general, when subtracting the pre-earthquake DEM from the postearthquake DEM at a pixel of the grid model, the negative values and the positive values represent the elevation changes in the source region and the deposition region, respectively. Then, based on the pre- and postearthquake DEM data, the difference of elevation () for each pixel can be calculated aswhere Dpre and Dpost are elevations before and after the earthquake, respectively (Figure 2). Assuming each pixel represents a unit area (A) in the DEM, the changes in volume driven by landsliding can be calculated using the following equation:where N is the number of pixels within a single landslide polygon, A is the area of a pixel, and is the elevation difference of ith pixel. This equation means that the volume of a landslide can be calculated by DEM changes of all pixels within the landslide polygon.

In general, when interpreting landslides from remote imagery, both source and deposit areas are included in one polygon. Therefore, the mass of material evacuated from the source area should be equal to the mass gained in the deposition area (i.e., net gain is equal to net loss). Thus, when equation (2) is used, we see the result of V = 0. As A is a constant value of the DEM pixel area, the value of V depends on . Within a landslide,  = 0 leads to V = 0. Although some studies show the expansion rate of a landslide as about 20% in cohesive soil-mantled regions due to long run-out distances [26, 27], we temporally neglect such mass expansion given that the study region is bedrock-mantled by the Pengguan massif rather than sedimentary rocks.

Before the calculation of elevation changes, the pre-earthquake DEM, which has a resolution of 30 m, is resized to the same resolution of 10 m similar to the postearthquake DEM. However, applying the DEM differential method faces a problem of background errors in data, mainly due to different DEM resolutions [10, 1820, 28]. For example, when summing up using the number of pixels within the extent of a landslide, the value of is always greater or less than zero. As shown in the map of coseismic elevation change, which is obtained by applying the cut-fill tool of ArcGIS, the topographic change is not consistent with the distribution of coseismic landslides (Figure 3). The errors here are obvious because most of the presented values are unreasonable in reality. For example, even in areas not impacted by landslides, the results still illustrate either net gain or net loss, implying the changes of elevation. This phenomenon differs from our assumptions, which should place this value near zero following the net loss from the source area which is equal to the net gain in the deposition area.

Based on the idea aforementioned, we use the following procedures to correct the errors in DEMs firstly and then calculate the landslide volume.

2.3.1. Calculating the Initial Difference between Pre- and Postearthquake DEMs

For a given landslide, the initial difference of elevation for each pixel, , is defined as . Here, represents the deposition area with net gain, while represents the source area with net loss (Figure 2). Due to systematic errors, is always far from zero, displaying both greater than and less than zero values. In this case, it is assumed that this nonzero value comes from all the pixels in the landslide.

2.3.2. Calculating the Background Error

In our assumption, for each landslide which includes both source and deposit areas, the net gain should be equal to the net loss when the DEMs are applied. Then, the error (i.e., ) is evenly distributed among the pixels within a landslide and forms a background error: . The following equation shows the procedure for calculating :

2.3.3. Calculating True Elevation Change

This means that the true elevation change for each pixel, namely, corrected difference () between the pre- and postearthquake DEM, is obtained by calculating the background error and subtracting this value from the initial difference in elevation, (Equation (5)). Figure 4 shows the process and the comparison between landslide values containing the background error and after removal of this error.

2.3.4. Calculating the Landslide Volume

After correcting background errors, the mass removed from the source area of a landslide is equal to the mass accumulated in the deposition area. From the corrected data, a landslide volume can be calculated by summing all pixels with positive values within the landslide area and multiplying this value by the area of one pixel (100 m2) (represented by purple dots in Figure 4), as shown in the following equation:

After calculation of volumes of all individual landslides, it is easy to obtain the volume of all landslides in the study area.

3. Results

Using the approach above, we calculate the volume of coseismic landslides in the study area. Figure 5(a) shows coseismic landslides with (before correction), and Figure 5(b) shows coseismic landslides with (after correction). Compared to Figure 5(a), results of Figure 5(b) are reasonable because they show that each landslide satisfies equation , namely, in a landslide, the mass moved from the source area is equal to the mass gained in the deposit area. It should be noted that, in Figure 5(b), as the same background error is subtracted by each pixel, it is not the reflection of real landslides (e.g., there are some net gain pixels appearing in the net loss areas).

Finally, the volume of coseismic landslides in the study area is calculated to be 1.765 × 108 m3.

In addition, based on the calculated volume and known area of individual landslides, this work uses equation (7) to perform a regression analysis:

The relation of log(A) versus log(V) has an R2 value of 0.8739 (Figure 6). Thus, the relation of V-A iswhich is largely consistent with those in the previous work (see Table 1).

Combining the result in Figure 6 with the landslide area-frequency relation shown in Figure 1(c), it demonstrates that most coseismic landslides in the study area are within the size range of 104–105 m2. There are 210 landslides falling in this area range, which account for 66% of the total number of the landslides in the study area. However, these landslides have a total volume of 0.677 × 108 m3, which accounts for only 38.4% of the total landslide volume. This implies that large-scale landslides are the primary contributor to the volume of coseismic landslides and have major influence on the calculation of this parameter. Considering these big slope failures are typically deep-seated, of which the depths are generally greater than the errors in elevations from DEMs, the method using elevation variations to calculate the landslide volume is more suitable for large slope failures, as they can produce surface changes larger than errors of DEMs.

4. Discussion

As a comparison, we use the empirical V-A formulas with α and γ from different studies to estimate the total landslide volume in the study area. Results are displayed in Table 1.

It can be seen that different α and γ values result in a wide range of calculated landslide volumes. Even for the formula derived from this work, the difference exists between volumes obtained from the V-A formula (1.388  108 m3 in Table 1) and the DEMs method (1.765 × 108 m3).

In fact, in the conventional relation , the functions of α and γ are analogous to landslide depths. With parameters A and V of landslides, using the various fitting methods to solve for depth is a basic principle. Therefore, the landslide area and quantity are important for obtaining a suitable regression formula. However, due to the subjective mapping by researchers, which defines the existence and/or extent of landslides [7, 8, 29], there exist obvious discrepancies in the quantity of landslides or the landslide boundaries, such as different landslide inventories compiled for the Wenchuan earthquake [7, 8, 29, 30]. In addition, coseismic landslide occurrence is subject to many factors including local topography, geology, and even local rock features. All these factors will greatly affect the coefficients α and γ when applying various fitting methods. For these reasons, it is not surprising that there are different α and γ in different studies (Table 1). Some researchers suggested that the existing V-A relations may overestimate the volumes of shallow landslides [18]. As for very large landslides, like the Daguangbao landslide triggered by the 2008 Wenchuan M7.9 earthquake, these relations can result in highly variable estimation values of landslide volumes [10]. Therefore, results of landslide morphometrics derived from various V-A relations are generally not comparable.

Compared to scaling relations derived from a limited number of measurements in traditional methods [5, 6], using DEMs is more powerful to develop methods to quantify changes in elevation [1820]. However, as mentioned previously, the elevation errors are always present in DEM data, and different DEM resolutions or DEM source data can cause errors when they are used for comparative analysis [10, 28]. It is difficult to have perfect matchable DEM data before and after a major earthquake [10]. Then, how to effectively apply the available DEM seems more important.

This study focuses on the landslide volume calculation, exploring how to eliminate the background errors in DEMs. However, individual landslides existing in areas both classified as net gain and net loss (Figure 3) mean that the background error is not a constant and cannot be generalized over an entire study region. So, it is relatively reasonable to consider the background error of each slope for calculation of the landslide volume. Figure 5(b) shows the changes after eliminating background errors with regard to individual landslides. When compared with Figure 3, it can be noted that both areas of net gain and net loss appear within a delineated landslide.

Unlike previous studies (e.g., Xu et al.) based on the deposition area of landslides determined by pre- and postearthquake DEM, this work takes the whole landslide (including both source and deposit area) as a target. Then, the principle that, in a landslide, the mass moved from the source area is equal to the mass gained in the deposit area can be applied effectively. This application permits to overlook unattainable data, such as landslide slip surface, clear boundaries between the source and deposit areas, and negate the issue of discrepancies within landslides inventories.

The limits of the method are mainly from the resolution of DEMs, and we do notice that this is still an open question to be addressed. Using pre- and postearthquake DEMs to calculate the landslide volume requires that the elevation difference should be greater than the depth of landslides, which ensures that the precision is sufficient for measuring the landslide volume by subtracting the pre-earthquake DEM from the postearthquake one [18]. This is why, this method is more suitable for large-scale landslides [20, 28], although both pre- and postearthquake DEMs inevitably still have elevation errors in this case. Unfortunately, the study area lacks postearthquake GPS data, so it is difficult to estimate the elevation error. Besides, this method removes the background error based on the difference of changes of the slope volume between the source and deposit areas, while material left at the source area after the landsliding is not taken into consideration because it did not move to the deposit area. Thus, the volume of the whole landslide may be underestimated. Nevertheless, many landslide cases show that only a small part of slope material was left at the source area after strong seismic shaking, so the overlay area has a little effect on the volume calculation for coseismic landslides.

In the case that deposit expansion rates of up to 20% of landslide masses as observed in some studies [26, 27], the DEM-based method may underestimate the deposit volume because of background error (), which is subtracted from the deposit area leading to a smaller volume. However, because the background error is evenly distributed through the whole landslide polygon, the impact of the background error () on the volume is very limited (<10%) compared to the error inherent to the DEM, which is around several meters [18]. In other words, the volume of coseismic landslides calculated by the DEM-based method may represent a minimum value ((Figure 7).

5. Conclusions

Calculating the volume of earthquake-triggered landslides is an important but difficult issue. Comparing pre- and postearthquake DEMs permits to calculate the elevation change by landsliding and then calculate the volume of landslides in a chosen area. This work tests this approach based on the data from the Yingxiu area, Sichuan, China, near the epicenter of the 2008 Mw 7.9 Wenchuan earthquake, which triggered massive coseismic slope failures. The conclusions of this study are as follows:(1)To use effectively the DEMs-based method to calculate the landslide volume, it is necessary to eliminate background errors of elevations in individual landslides firstly, which are associated with different sources and resolutions of pre- and postearthquake DEMs. With such a preprocessing, the calculated landslide volume may be more reasonable.(2)Using the DEMs-based method, the total landslide volume in the study area is 1.765 × 108 m3. Considering the simplification in this approach, the volume of coseismic landslides calculated may represent a minimum value.(3)By a log-linear regression, a relation between the landslide volume and area for the study area is derived, largely consistent with that in previous work. Combining the frequency-area and volume-area relations in the study area, it is noted that the large-scale landslides are the major contributor to the calculated volume. And, the DEMs-based method is more suitable for the calculation of large-scale earthquake-triggered landslides.

Data Availability

The data used to support the findings of this study are available at http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Key Research and Development Program (no. 2017YFC1501004) and Basic Scientific Research Project of Institute of Geology, China Earthquake Administration (no. IGCEA1604).