Next Article in Journal
Urban Link Travel Time Prediction Based on a Gradient Boosting Method Considering Spatiotemporal Correlations
Next Article in Special Issue
Detection of Catchment-Scale Gully-Affected Areas Using Unmanned Aerial Vehicle (UAV) on the Chinese Loess Plateau
Previous Article in Journal
A Long Baseline Three Carrier Ambiguity Resolution with a New Ionospheric Constraint
Previous Article in Special Issue
A GIS Study of the Influences of Warm Ocean Eddies on the Intensity Variations of Tropical Cyclones in the South China Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gully Erosion Mapping and Monitoring at Multiple Scales Based on Multi-Source Remote Sensing Data of the Sancha River Catchment, Northeast China

1
College of Earth Science, Jilin University, Changchun 130061, China
2
Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun 130102, China
3
College of Arts & Sciences, Beijing Union University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2016, 5(11), 200; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi5110200
Submission received: 19 July 2016 / Revised: 25 October 2016 / Accepted: 31 October 2016 / Published: 4 November 2016

Abstract

:
This research is focused on gully erosion mapping and monitoring at multiple spatial scales using multi-source remote sensing data of the Sancha River catchment in Northeast China, where gullies extend over a vast area. A high resolution satellite image (Pleiades 1A, 0.7 m) was used to obtain the spatial distribution of the gullies of the overall basin. Image visual interpretation with field verification was employed to map the geometric gully features and evaluate gully erosion as well as the topographic differentiation characteristics. Unmanned Aerial Vehicle (UAV) remote sensing data and the 3D photo-reconstruction method were employed for detailed gully mapping at a site scale. The results showed that: (1) the sub-meter image showed a strong ability in the recognition of various gully types and obtained satisfactory results, and the topographic factors of elevation, slope and slope aspects exerted significant influence on the gully spatial distribution at the catchment scale; and (2) at a more detailed site scale, UAV imagery combined with 3D photo-reconstruction provided a Digital Surface Model (DSM) and ortho-image at the centimeter level as well as a detailed 3D model. The resulting products revealed the area of agricultural utilization and its shaping by human agricultural activities and water erosion in detail, and also provided the gully volume. The present study indicates that using multi-source remote sensing data, including satellite and UAV imagery simultaneously, results in an effective assessment of gully erosion over multiple spatial scales. The combined approach should be continued to regularly monitor gully erosion to understand the erosion process and its relationship with the environment from a comprehensive perspective.

1. Introduction

Gully erosion is the process in which soil is removed according to the concentration of surface water and subsurface water in narrow flow paths, resulting in the formation of incised channels that may grow into gullies deeper than 30 cm over short time periods [1]. This causes significant soil loss and soil degradation, and it is generally considered as an indicator of desertification and land degradation [2]. Several studies proposed that gully erosion is the major soil erosion type in the black soil region of Northeast China, which suffered from severe soil erosion [3,4,5]. Despite being a natural process, gully erosion is generally initiated and accelerated by inappropriate agricultural and human activities [6]. Gully erosion mapping is a crucial step to monitor land degradation and study its current and future local impacts.
To date, numerous gully erosion studies have been performed at different scales and with different objectives [7]. For instance, gully location maps at the national scale [8] and gully risk map at regional scale have been obtained [9], the effect of topography on gully distributions [5] and the spatial variations of gully distributions over a large scale have been analyzed [10], gully networks have been defined for small catchments at a medium scale [11], and gully headcut retreat rates at a small scale have been estimated [12]. Various techniques have been applied for measuring, mapping and monitoring gully erosion accordingly. Compared with conventional techniques (e.g., ruler, tape, microtopographic profilers [13], poles, total stations [14], pins [15] and differential GPS [16,17]), which are time-consuming for achieving high accuracy in field surveys at a small scale, the development of remote sensing techniques has provided an efficient method to obtain spatially continuous gully information over large scales for different time periods. Classical aerial photography has been effectively employed for large-scale areas and long-term investigations, but with low spatial and temporal resolution [18,19]. Airborne LiDAR [20,21,22] and terrestrial laser scanner (TLS) [23,24,25] derive high resolution and precision data products for precise modeling of land surface. They offer the ability to measure gullies contactless with a high spatial resolution but require significant expertise. Analysis of satellite images integrated with field surveys has proven to be an effective and practical approach for mapping gully erosion over large areas [5]. However, the spatial resolution of satellite images in previous studies was mostly higher than 5 m, which results in omitting a considerable quantity of ephemeral gullies [26] and does not allow for a detailed analysis [27]. In addition, small changes of an individual gully within short time periods are difficult to identify using satellite images [27].
To increase ground resolution for the monitoring of short-term process, lower altitude unmanned aerial platforms, e.g., blimps or kites, are starting to be employed [28,29]. Unmanned Aerial Vehicles (UAVs) have integrated autopilot technology, which allows for semi- or fully autonomous navigation, flight control and image acquisition capabilities [30]. Researchers are beginning to explore the use of UAVs for the acquisition of remote sensing imagery and application in Earth surface analysis with a higher level of detail [31,32]. Several authors have already employed high resolution DSMs of cm-accuracy obtained from UAV-based imagery to monitor and quantify gully erosion [27,33,34]. Moreover, automatic 3D photo-reconstruction techniques have shown great advantages over traditional photogrammetry techniques by making image processing remarkably easier when images are derived from uncalibrated and nonmetric cameras [35]. Compared with laser scanners and rigorous photogrammetric methods, Structure-from-motion (SfM) and multi-view-stereo (MVS) algorithms, which are the generally proposed reconstruction methods, just use consumer-grade cameras and require less expertise [36,37]. Previous work has compared 3D photo-reconstruction with existing 2-D and 3-D field methods (e.g., LiDAR, laser profilometer, total station) for assessing gully erosion with regard to cost, accuracy and effectiveness, revealing the superiority of 3D photo-reconstruction [14,38]. By adopting 3D photo-reconstruction techniques, numerous researchers used images derived from UAV to generate high resolution topographic data and ortho-image for monitoring gully erosion [34,39].
It is still a major challenge in gully erosion research to develop and apply appropriate measuring techniques for monitoring or better understanding the initiation and change of various gully types at various spatial and temporal scales [40]. Most of the previous gully erosion studies were conducted at a single spatial scale. However, rare studies focused on the mapping and monitoring of gullies over multiple spatial scales which would contribute to a better understanding of gully erosion dynamics and their interact with change of gully environments from a comprehensive perspective [41]. This study aims to map and monitor gully erosion over multiple spatial scales using multi-source remote sensing data. High resolution satellite image visual interpretation with field verification was used to map gully geometric features and evaluate gully erosion at a catchment scale. UAV imagery combined with the 3D photo-reconstruction method was employed for gully mapping at a site scale in further detail.

2. Materials

2.1. Study Area

The study area (126°8′49″E–126°17′19″E, 44°9′18″N–44°22′20″N), including the middle and upper reaches of the Sancha River basin (Figure 1c), is located in Northeast China (Figure 1a,b), covering a total area of 163.83 km2. This area is situated at the southeastern edge of the black soil region (Figure 1b) and is covered by typical black soil, meadow soil, albic soil and dark brown soil. The region has a continental climate with a mean annual temperature of 5.3 °C and annual precipitation of 550–600 mm. Long-term human activities, especially excessive reclamation, resulted in a sloped farmland, with little vegetation coverage in the area. Due to the substantial terrain undulation, this region has suffered from severe soil erosion by water. As one of the most critical and destructive water erosion types, gully erosion is widely distributed in this area. Hence, the study area is one of the typical black soil erosion regions of Northeast China.

2.2. Data Sources

The sources of the remote sensing data included satellite imagery and UAV images. Area-wide gully mapping was based on Pleiades-1A imagery (panchromatic sharpened images at 0.7 m resolution fused with 2.8 m multi-spectral bands). The image was acquired on 12 October 2013 (Figure 1c). The geometric, radiometric and atmospheric corrections were performed by the data provider, the Jilin Province Geomatics Center.
The UAV images were acquired in the airborne campaign, which was conducted in October 2015 by using a battery-powered quadrocopter type Phantom 3 Professional (DJI, Shenzhen, China). The quadrocopter has four rotors, a weight of 1280 g and a maximum flight time of 23 min. The UAV was equipped with an integrated and stabilized camera. The sensor Sony EXMOR 1/2.3″ was utilized for image acquisition. The detailed image acquisition workflow is described in Section 3.2. These UAV images were used for detailed mapping at a small site.
Remote sensing data need to be supported and validated with ground truth data, which is achieved through field surveys. We attempted to achieve the following goals through field investigation: (1) determine the gully interpretation signs (i.e., possible gully appearance, such as gully size and shape) in the Pleiades images; (2) verify the mapping accuracy and further correct errors of omission and commission; and (3) conduct an UAV airborne campaign and install Ground Control Points (Figure 1d). According to the above goals, this study conducted field surveys in May and October 2015. Fieldwork was supported with the Jilin Continuously Operating Reference Station (JLCORS) installed at a Trimble device; the horizontal and vertical positioning accuracy of the device is better than 3 cm and 5 cm, respectively. The gully interpretation signs in the high-resolution satellite images were determined (Table 1). To verify the accuracy of the area-wide mapping, field validation work was conducted, and gullies located adjacent to the road were identified.
Other important geographical condition data (road, water, residential sites, etc.) were also provided by the Jilin Province Geomatics Center. The Digital Elevation Model (DEM) of the 5-m pixel size was produced using 1:10,000 scale topographical maps.

3. Methods

3.1. Pleiades-1A Imagery for Area-Wide Mapping

Visual interpretation and vectorization were used to obtain gully distribution data of the study area. Table 1 shows four types of possible gullies in the high-resolution satellite images and ground photos of the study area. Gully-1 is typical of ephemeral gullies, which generally have long narrow structures and occur on farmland. Ephemeral gullies are small channels eroded by concentrated flow that can be easily filled by normal tillage, but reappear in the same location by additional runoff events [2]. As Gully-2 shows, the further developed linear gullies are much more deeply incised. The channels are too deep and wide to level with ordinary farm tillage equipment and are referred to as permanent gullies [2]. Gully-3 represents a larger and wider permanent gully, in which the gully bottom is sparsely covered by plants. The gully headcut and steeper gully walls further decrease in the adjacent flat lands. Compared with the single gully of Gully-3, a ramified gully system with several branches can be observed in Gully-4. Between the individual branches, smaller flat parts of the original surface remain, and have even been reclaimed. However, gullies in forests were not the object of this study since these gullies are difficult to recognize from optical satellite images.
The gullies are mainly discriminated from the surfaces they incise according to their spatial characteristics, i.e., patterns (drainage), shape and size. However, the gullies are difficult to differentiate from the surrounding environment based on the spectral properties because the spectral reflectance between gullies varies markedly over large areas [8,42]. Meanwhile, the spectral properties of gullies rely on vegetation cover or soil, but the spectral values of the bare soil or vegetation within the gully may be the same as those outside. It is gratifying that the shadow cast by the gully walls contribute to the distinction of gullies from their surroundings. Since the incident angle of sunlight varies relatively from the landform surface, shadow cast is visible in high resolution optical satellite images for most gullies. The shadow of the Sun-diverted gully walls presents a spectral contrast between the surroundings and the sunlit opposite gully walls [43]. For example, Gully-2 is more discernible than Gully-1 due to the shadowing. Accordingly, the shadow helps to present the linear characteristics and shape of the gullies by highlighting their edges. In addition, the topographical map, hydrographical network map, traffic map and related information contribute to recognizing erosion gullies.
The gully vectorization was performed using ArcGIS 10.0 software. The gully distribution data included gully length and area. Gullies with width less than 3 pixels were only outlined for length, not area. According to the field validation result, the gully data were modified. The omitted and committed gullies were included and excluded, respectively. Then, the final gully distribution data were obtained. In the overall assessment and validation of gully erosion, gullies were not classified. To assess gully erosion comprehensively, the indicators of gully density and intensity were adopted. Gully density within the study area was defined as the total gully length per area of the whole study area. The gully density distribution was estimated by using the Focal Statistics tool in ArcGIS. The specific steps are as follows: (1) a fishnet polygon containing cells of 5 m × 5 m was created, then the total length per area of each cell was calculated and a value was assigned to the corresponding cell; (2) the polygon was converted into a grid with a 5-m pixel size; and (3) Focal Statistics was applied to calculate the average value of each input cell within a rectangular neighborhood of 100 m × 100 m; thus, the distribution map of gully density was obtained. Similarly, gully intensity was calculated by using the measurement of total gully area per area. Furthermore, the topographic differentiation of gully erosion was analyzed using the DEM and the zonal statistics tool in ArcGIS.

3.2. UAV-Based Remote Sensing for Detailed Mapping

A small study site (Figure 1d) in the Sancha River catchment was selected for detailed mapping using UAV, as it includes typical gullies (for example, Gully-4) and is located near the road. The image acquisition workflow is as follows: Firstly, Ground Control Points (GCPs) were distributed over the chosen study site prior to the airborne campaign (Figure 1d). The GCPs were installed and precisely measured with JLCORS, as well as used for georeferencing during further image processing. The geographic coordinate system and projection were set to China Geodetic Coordinate System 2000 (CGCS 2000) and Gauss Kruger, respectively. Secondly, the UAV was operated in the field by a trained operator and controlled autonomously by the autopilot software DJI GO (DJI, Shenzhen, China) and Pix4Dcapture (Pix4d, Lausanne, Switzerland) to follow predefined flight plans for imagery acquisition during takeoff, flight and landing. The flight plans were designed for longitudinal overlap of 80% and lateral overlap of 60% to avoid gaps between single photographs and ensure successful photo-reconstruction. The lens of the camera was set in the 90-degree position, and the flying height was set at 120 m above the ground for the duration of the imaging process.
As one of the main aims of this research, the ortho-image mosaic and Digital Surface Model (DSM) were generated through the 3D photo-reconstruction method by new photogrammetry software in an automated way. 3D reconstruction using structure-from-motion (SfM) and multi-view stereo (MVS) techniques was performed for surface geomorphometry reconstruction to generate 3D data from a series of non-calibrated images. In general, the reconstruction workflow includes the processing steps: (1) identification of homologous image points and consequent image matching; (2) reconstruction of the image acquisition configuration, including camera orientation, position and internal camera parameters, as well as corresponding 3D coordinates, with an iterative bundle block adjustment (BBA); and (3) dense matching of the sparse point cloud by forward intersections in space [34,36,44]. In this study, the commercial photogrammetry software, Pix4DMapper (Pix4D, Lausanne, Switzerland), based on SfM and MVS techniques, was used to automatically generate an ortho-mosaic and DSM. In addition, although each image was recorded with the synchronized GPS position of low accuracy, high-precision GCPs were necessary to be included into the processing by manually adding additional tie points and using an automated tie point generation program within Pix4dmapper software. As an important gully indicator, the volume was calculated based on previous studies [27,33]. The 3D polyline of the gully edge was used to create a 3D polygon. We assumed that the previous surface was completely flat between the current gully edges, so the polygon was then transformed to a raster of the original former surface prior to erosion. The gully DSM was further subtracted. The resulting raster contained the height difference between the state before incision and the current state for each raster cell. These values were then multiplied with the raster resolution and summed up, resulting in the gully volume. The calculation of the gully volume was conducted in ArcGIS software, and the method was described in detail by d’Oleire-Oltmanns et al. [27] and Peter et al. [33].

4. Results

4.1. Spatial Distribution of Erosion Gullies in the Catchment

According to visual interpretation of very high resolution satellite images, the status of gully erosion in the whole catchment area was obtained. The statistical results of gully erosion are shown in Table 2. In 2013, there were 1190 erosion gullies in the study area, among which 572 (nearly half of the total) were erosion gullies, with width less than 3 pixels. This result indicated that the gully erosion is still relatively active in the study area. We assumed that these narrow gullies would further develop with intensified erosion without any prevention. The gully density of the whole area was 1.37 km/km2 (gully length per unit area) and the gully intensity was 5434.84 m2/km2 (gully area per unit area). The total gully area and length were 89.039 hm2 and 224.978 km, respectively. According to the standard for classification and gradation of soil erosion (SL190-2007) issued by the Ministry of Water Resources of China, the gully density of the whole area was classified as slight erosion.
Figure 2 displays the spatial distribution of gully erosion across the Sancha River catchment. Figure 2a illustrates that gullies are widespread across the catchment, especially in the northeastern part. Figure 2b,c shows that the gully density and gully intensity ranged from 0 to 71.99 km/km2 and from 0 to 792,896 m2/km2, respectively. The results all indicate that gully erosion in the northeastern part was most serious. However, there are differences between the spatial distribution of gully density and intensity. The high gully intensity generally indicates that the lateral erosion of the gully erosion is significant, which results in notable widening of ditch sides and an increase in the area.
In the field validation work (Table 3), 20 gullies were observed and later recorded by JLCORS. The indoor study showed that 18 gullies were captured correctly, while a roadside gully was omitted and another field path was mistaken for an ephemeral gully. Nevertheless, the overall accuracy of gully interpretation reached 90%, and the very high resolution (sub-meter) satellite image showed a strong ability in the detection of various gully types.

4.2. High-Resolution Ortho-Image Mosaic and DSM of a Small Site

A total number of 256 JPEG photos were acquired during the field airborne campaign; 255 out of 256 images were calibrated (99%) and all images were enabled. The ortho-image mosaic (Figure 3) and DSM (Figure 4) were generated for the small study site, covering a total of 13.85 ha. The average ground sampling distance (GSD) of pixels was 4.42 cm. Forty installed GCPs were used for ground control. The root mean square errors (RMSE) of the GCPs were between 0.026 and 0.016 m in the horizontal directions, and 0.024 m in the vertical direction. A total of 2,622,204 2D keypoint observations and 899,037 3D points were used for the bundle block adjustment, while the mean reprojection error was evaluated at approximately 0.26 pixel. Table 4 shows the summary of these outputs for the results obtained in Pix4dmapper software using automatic mode.
The resulting ortho-image mosaics and DSM reveal the area of agricultural utilization and its shaping by human agricultural activities and water erosion in detail. The maize straws with accumulation on the flat farmland are visible along with landfill within the headcut. The cultivated land in the bottom of the gully is even distinguishable. In addition, the depression of existing gully erosion is highly visible. In addition to the three main gullies, a roadside gully and small rills are visible. In Figure 4, the light blue part illustrates the lowest area at the Wadi bottom. However, the maximum incision depth of the gullies reaches 10.32 m in the upper leftmost main gully (i.e., Gully-4 mentioned above). The gully volumes were calculated as 68,145 m3, 12,240 m3 and 5405 m3 for the three main gullies (from north to south in Figure 4, respectively). Meanwhile, a precise and comprehensive 3D model of the gullies was calculated (Figure 5). The enlarged inset shows the 3D structure of the gully sidewall and headcut in detail. The results shown in Figure 3, Figure 4 and Figure 5 illustrate the high degree of detail achievable with the method of UAV-based remote sensing and 3D photo-reconstruction.

5. Discussion

5.1. Visual Pleiades Imagery Interpretation for Gully Erosion Mapping

The Pleiades imagery displayed high feasibility for the recognition of various gully types, especially ephemeral gullies in this study (e.g., Gully-1 in Table 1). However, with the enhancement of satellite image resolution, the phenomena of “noise” and “different objects with the same spectrum” as well as “the same objects with different spectra” increase [45]. In the present study, a small farm path was easily mistaken for an ephemeral gully because their spectral characteristics and geomorphological structures were similar in the Pleiades imagery. Meanwhile, it is difficult to distinguish narrow roadside gullies from the road, as they were too similar to be differentiated from the imagery. Nevertheless, the obtained interpretation result was satisfactory due to the high accuracy (Table 3). However, achieving the satisfactory result was time consuming and required a certain expertise. According to previous studies, with respect to high spatial resolution satellite images (such as Pleiades-1A), although the visual interpretation of gully erosion with field verification achieved much higher accuracies [8,46], (semi-)automatic mapping of gullies at a regional scale or larger scales showed great potential for obtaining a general statement of the distribution of gullies to a wider extent [43,47,48].
The higher spatial resolution of satellite images made it possible to identify gullies with complex surroundings as well as with a wide range of factors that control the development and distribution of gully erosion [49,50]. The topographic differentiation characteristics of gully erosion at a large scale were obvious (see Figure 6). Figure 6a illustrates a layered distribution of gully erosion in the vertical direction. This result indicates that the slope land between geomorphic types of lower elevations (e.g., valley plain) and higher elevations (e.g., hill) is more conducive to the formation and development of gully erosion due to a longer slope and larger catchment area [5]. Figure 6b shows that gully density and intensity all increase first and then decrease as the slope increases. Although the steeper slope will increase the shear force of the soil particles [51], the slope is not the main factor restricting the form of erosion gully, which was more affected by slope length, catchment area, slope shape and other terrain factors [52] when the slope was greater than six degrees. The relation between slope aspect and gully erosion was analyzed (Figure 6c). The results indicate that the gully erosion on windward slopes (S, SSW, SW, WSW, and W) was greater than on leeward slopes because the prevailing southwest winds result in stronger rainfall erosivity on windward slopes than on leeward slopes. Furthermore, the difference between gully density and intensity on various terrain factors should be given further consideration.

5.2. Assessment of UAV Remote Sensing and 3D Reconstruction

Numerous studies already illustrated that the low-altitude UAV photogrammetry and 3D photo-reconstruction method showed the definite advantage in the achievable detail with very high resolution [27,34,39,53]. The extremely detailed measurements were generally obtained on a sub-decimeter level, which is comparable to TLS [14]. TLS is often used as high quality reference spatial dataset to assess different field methods for measuring gully erosion due to its high accuracy [54,55]. In this study, accurate TLS or LiDAR are unavailable. In other words, the spatially continuous precise ground truth data are unavailable as a reference for the chosen study site. Thus, the traditional accuracy assessment is impossible, as already illustrated by Harwin and Lucieer [56]. This study used the high-precision GCPs as independent checkpoints to assess the accuracy automatically (see Table 4) and the same assessment method is also applied by d’Oleire-Oltmanns et al. [27]. Therefore, the resulting accuracy of the derived products is considered to be reliable. Note that the obtained value for the gully volume should not be regarded as a precise absolute erosion amount. The quantification values obtained at multiple periods should be a good estimation of the magnitude of erosion-induced change [34]. Moreover, the manual digitizing of the gully edge influences the accuracy of the volumes and defining the location of the gully edge still remains a challenge [57].

5.3. Comparison and Combination of Satellite Imagery and UAV Remote Sensing

Even though the very high resolution Pleiades-1A satellite imagery has significantly improved the identification of gully erosion, it does not allow for a more detailed analysis of an individual gully compared with UAV imagery. For instance, the use of corn straw mulch buffer within the gully headcut is a simple and common gully prevention measure in sloping farmland in the black soil region of Northeast China [58]. Pleiades-1A with 0.7 m does not provide more detailed information about the corn straws within the gully headcut. However, the corn straw mulch is visible in the ortho-image derived from UAV (0.042 m) (Figure 7). According to the enlarged view of the UAV image (the lower right of Figure 7), we can see the vegetation in the bottom clearly, and even distinguish the active gully headcut, which was eroding the farmland.
In addition to the gully length and width, the gully volume is also an important factor for the determination of sediment yield from gully erosion. The very high resolution DSM derived from UAV-based remote sensing and 3D photo-reconstruction was used for quantification of gully volumes in detail (Figure 4). In contrast, the resolution of the available DEM, i.e., 5 m × 5 m resolution, is too coarse and does not allow this degree of detail. Meanwhile, it is difficult to extract gully volumes from the single Pleiades imagery, and stereoscopic Pleiades images are not available in this study. In addition, the 3D model provided more details about the 3D structure of the gully sidewall and headcut (Figure 5), but 3D information was lacking in the Pleiades image. The high resolution and precision of the data products derived from UAV enable detailed 2D and 3D analysis, which is not possible at the satellite scale level.
The previous gully mapping and monitoring studies were generally conducted at a single spatial scale to achieve a certain goal. Meanwhile, the data types achieved from different scales were too complex to be integrated properly. However, UAV-based remote sensing provides a valuable intermediate scale between the satellite scale and field scale [27]. The UAV-based remote sensing data can provide detailed 2D and 3D analysis of an individual gully at a small site scale, which is not possible at the satellite scale level. The satellite image can provide the surrounding environment of the individual gully. The obtained distribution data of various types of gullies at a large scale can be used for analyzing how the topographic factors control the formation and development of gully erosion. The two sources of remote sensing data can compare and complement each other at multiple scales.
Although gully erosion mapping and monitoring were conducted at a single time in this study, small changes of an individual gully within short time periods can already be monitored by applying short term UAS-based remote sensing [34]. In addition, the regional variation of the spatial distribution of erosion gullies over decades can also be analyzed by using satellite images [10]. The combined approach will be continued to become a regular monitoring for gully erosion. Numerous gully erosion studies around the world have been performed to monitor or better understand the initiation and development of various gully types at various temporal and spatial scales [4,7,20,27,34]. Compared with a single scale, the approach based on multi-source remote sensing data over multiple scale levels to combine the very detailed field scale and the larger scale level can help to elucidate the erosion process and its relationship with the environment at different spatial and temporal scales from a comprehensive perspective.

6. Conclusions

The present study was one of the first to attempt to map and monitor gully erosion over multiple scales by combining satellite and UAV imagery. At the catchment scale, the sub-meter images show strong ability in the recognition of various gully types and obtained satisfactory results. The gully distribution of the whole Sancha River catchment was derived, and the topographic differentiation characteristics of gully erosion over a large scale suggested that the gully erosion was stronger on windward sunny slopes than on leeward slopes and mainly occurred at middle elevations. When the slope was greater than six degrees, the slope was not the main factor restricting gully erosion. At a more detailed site scale, UAV imagery combined with 3D photo-reconstruction provided DSM and ortho-image data at a centimeter level as well as a 3D model. The resulting products revealed the area of agricultural utilization and its shaping by human agricultural activities and water erosion in detail and also provided the gully volume. However, detailed 2D and 3D analysis were lacking at the satellite scale level.
Using the Sancha River catchment as the study area, we illustrated that using multi-source remote sensing data, including satellite and UAV imagery simultaneously, results in an effective assessment of gully erosion over multiple scales. The combined mapping approach will be continued to regularly monitor the area and aid in further understanding the erosion process and its relationship with the environment from a comprehensive perspective. The work offers guidance for policy makers to develop sustainable comprehensive soil and water conservation management plans in the study area.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (No. 41271416; 41301467, and 41601458) and the Project 2016213 supported by the Graduate Innovation Fund of Jilin University. The authors would like to thank the Jilin Province Geomatics Center and its director Zheng Wang for their support with the fieldwork.

Author Contributions

Ranghu Wang, Shuwen Zhang and Jiuchun Yang designed the research; Ranghu Wang and Luoman Pu processed and analyzed the satellite data; Ranghu Wang, Luoman Pu, Chaobin Yang, Jing Chen, Cong Guan, Qing Wang and Dan Chen collected the UAV data; Ranghu Wang and Xuejia Sang completed the 3D reconstruction; and Ranghu Wang wrote the paper and Bolin Fu assisted the revision.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Luffman, I.E.; Nandi, A.; Spiegel, T. Gully morphology, hillslope erosion, and precipitation characteristics in the Appalachian Valley and Ridge province, southeastern USA. Catena 2015, 133, 221–232. [Google Scholar] [CrossRef]
  2. Poesen, J.; Nachtergaele, J.; Verstraeten, G.; Valentin, C. Gully erosion and environmental change: Importance and research needs. Catena 2003, 50, 91–133. [Google Scholar] [CrossRef]
  3. Liu, B.; Yan, B.; Shen, B.; Wang, Z.; Wei, X. Current status and comprehensive control strategies of soil erosion for cultivated land in the northeastern black soil area of China. Sci. Soil Water Conserv. 2008, 6, 1–8. (In Chinese) [Google Scholar]
  4. Hu, G.; Wu, Y.; Liu, B.; Yu, Z.; You, Z.; Zhang, Y. Short-term gully retreat rates over rolling hill areas in black soil of Northeast China. Catena 2007, 71, 321–329. [Google Scholar] [CrossRef]
  5. Zhang, S.; Li, F.; Li, T.; Yang, J.; Bu, K.; Chang, L.; Wang, W.; Yan, Y. Remote sensing monitoring of gullies on a regional scale: A case study of Kebai region in Heilongjiang Province, China. Chin. Geogr. Sci. 2015, 25, 602–611. [Google Scholar] [CrossRef]
  6. Kakembo, V.; Rowntree, K.M. The relationship between land use and soil erosion in the communal lands near Peddie town, Eastern Cape, South Africa. Land Degrad. Dev. 2003, 14, 39–49. [Google Scholar] [CrossRef]
  7. Gomez-Gutierrez, A.; Conoscenti, C.; Angileri, S.E.; Rotigliano, E.; Schnabel, S. Using topographical attributes to evaluate gully erosion proneness (susceptibility) in two Mediterranean basins: Advantages and limitations. Nat. Hazards 2015, 79, S291–S314. [Google Scholar] [CrossRef]
  8. Roux, N.M.; Le, J.J. Gully location mapping at a national scale for South Africa. S. Afr. Geogr. J. 2012, 94, 208–218. [Google Scholar]
  9. Eustace, A.H.; Pringle, M.J.; Denham, R.J. A risk map for gully locations in central Queensland, Australia. Eur. J. Soil Sci. 2011, 62, 431–441. [Google Scholar] [CrossRef]
  10. Yan, Y.C.; Zhang, S.W.; Yue, S.P. Application of corona and spot imagery on erosion gully research in typical black soil regions of Northeast China. Resour. Sci. 2006, 27, 650–656. (In Chinese) [Google Scholar]
  11. Perroy, R.L.; Bookhagen, B.; Asner, G.P.; Chadwick, O.A. Comparison of gully erosion estimates using airborne and ground-based LiDAR on Santa Cruz Island, California. Geomorphology 2010, 118, 288–300. [Google Scholar] [CrossRef]
  12. Wu, Y.; Zheng, Q.; Zhang, Y.; Liu, B.; Cheng, H.; Wang, Y. Development of gullies and sediment production in the black soil region of Northeastern China. Geomorphology 2008, 101, 683–691. [Google Scholar] [CrossRef]
  13. Casali, J.; Loizu, J.; Campo, M.A.; De Santisteban, L.M.; Alvarez-Mozos, J. Accuracy of methods for field assessment of rill and ephemeral gully erosion. Catena 2006, 67, 128–138. [Google Scholar] [CrossRef]
  14. Castillo, C.; Perez, R.; James, M.R.; Quinton, J.N.; Taguas, E.V.; Gomez, J.A. Comparing the accuracy of several field methods for measuring gully erosion. Soil Sci. Soc. Am. J. 2012, 76, 1319–1332. [Google Scholar] [CrossRef] [Green Version]
  15. Ionita, I. Gully development in the Moldavian Plateau of Romania. Catena 2006, 68, 133–140. [Google Scholar] [CrossRef]
  16. Wu, Y.; Cheng, H. Monitoring of gully erosion on the Loess Plateau of China using a global positioning system. Catena 2005, 63, 154–166. [Google Scholar] [CrossRef]
  17. Hu, G.; Wu, Y.Q.; Liu, B.Y.; Zhang, Y.G.; You, Z.M.; Yu, Z.T. The characteristics of gully erosion over rolling hilly black soil areas of Northeast China. J. Geogr. Sci. 2009, 19, 309–320. [Google Scholar] [CrossRef]
  18. Vandekerckhove, L.; Poesen, J.; Govers, G. Medium-term gully headcut retreat rates in Southeast Spain determined from aerial photographs and ground measurements. Catena 2003, 50, 329–352. [Google Scholar] [CrossRef]
  19. Bouchnak, H.; Felfoul, M.S.; Boussema, M.R.; Snane, M.H. Slope and rainfall effects on the volume of sediment yield by gully erosion in the Souar lithologic formation (Tunisia). Catena 2009, 78, 170–177. [Google Scholar] [CrossRef]
  20. James, L.A.; Watson, D.G.; Hansen, W.F. Using LiDAR data to map gullies and headwater streams under forest canopy: South Carolina, USA. Catena 2007, 71, 132–144. [Google Scholar] [CrossRef]
  21. Jackson, T.J.; Ritchie, J.C.; White, J.; Leschack, L. Airborne laser profile data for measuring ephemeral gully erosion. Photogramm. Eng. Remote Sens. 1988, 54, 1181–1185. [Google Scholar]
  22. Armston, J.; Disney, M.; Lewis, P.; Scarth, P.; Phinn, S.; Lucas, R.; Bunting, P.; Goodwin, N. Direct retrieval of canopy gap probability using airborne waveform LiDAR. Remote Sens. Environ. 2013, 134, 24–38. [Google Scholar] [CrossRef]
  23. Goodwin, N.R.; Armston, J.; Stiller, I.; Muir, J. Assessing the repeatability of terrestrial laser scanning for monitoring gully topography: A case study from Aratula, Queensland, Australia. Geomorphology 2016, 262, 24–36. [Google Scholar] [CrossRef]
  24. Kociuba, W.; Janicki, G.; Rodzik, J. 3D laser scanning as a new tool of assessment of erosion rates in forested loess gullies (case study: Kolonia Celejów, Lublin Upland). Ann. UMCS 2014, 69, 107–116. [Google Scholar] [CrossRef]
  25. Kociuba, W.; Kubisz, W.; Zagorski, P. Use of terrestrial laser scanning (TLS) for monitoring and modelling of geomorphic processes and phenomena at a small and medium spatial scale in polar environment (Scott River—Spitsbergen). Geomorphology 2014, 212, 84–96. [Google Scholar] [CrossRef]
  26. Li, H.; Zhang, X.; Liu, S.; Yury, K.; Kateryna, I.; Li, X.; Chen, Q. Evolvement of gully erosion in village scale in the typical black soil area. Sci. Soil Water Conserv. 2012, 10, 21–28. (In Chinese) [Google Scholar]
  27. D’Oleire-Oltmanns, S.; Marzolff, I.; Peter, K.D.; Ries, J.B. Unmanned aerial vehicle (UAV) for monitoring soil erosion in Morocco. Remote Sens.-Basel 2012, 4, 3390–3416. [Google Scholar] [CrossRef]
  28. Marzolff, I.; Poesen, J. The potential of 3d gully monitoring with GIS using high-resolution aerial photography and a digital photogrammetry system. Geomorphology 2009, 111, 48–60. [Google Scholar] [CrossRef]
  29. Ries, J.B.; Marzolff, I. Monitoring of gully erosion in the Central Ebro Basin by large-scale aerial photography taken from a remotely controlled blimp. Catena 2003, 50, 309–328. [Google Scholar] [CrossRef]
  30. Hugenholtz, C.H.; Whitehead, K.; Brown, O.W.; Barchyn, T.E.; Moorman, B.J.; LeClair, A.; Riddell, K.; Hamilton, T. Geomorphological mapping with a small unmanned aircraft system (sUAS): Feature detection and accuracy assessment of a photogrammetrically-derived digital terrain model. Geomorphology 2013, 194, 16–24. [Google Scholar] [CrossRef]
  31. Haas, F.; Hilger, L.; Neugirg, F.; Umstädter, K.; Breitung, C.; Fischer, P.; Hilger, P.; Heckmann, T.; Dusik, J.; Kaiser, A. Quantification and analysis of geomorphic processes on a recultivated iron ore mine on the Italian island Elba using long-time ground-based LiDAR and photogrammetric data by an UAV. Nat. Hazards Earth Syst. Sci. 2015, 3, 6271–6319. [Google Scholar] [CrossRef]
  32. D’Oleire-Oltmanns, S.; Eisank, C.; Dragut, L.; Blaschke, T. An object-based workflow to extract landforms at multiple scales from two distinct data types. IEEE Geosci. Remote Sens. 2013, 10, 947–951. [Google Scholar] [CrossRef]
  33. Peter, K.D.; d’Oleire-Oltmanns, S.; Ries, J.B.; Marzolff, I.; Hssaine, A.A. Soil erosion in gully catchments affected by land-levelling measures in the Souss Basin, Morocco, analysed by rainfall simulation and UAV remote sensing data. Catena 2014, 113, 24–40. [Google Scholar] [CrossRef]
  34. Stocker, C.; Eltner, A.; Karrasch, P. Measuring gullies by synergetic application of UAV and close range photogrammetry—A case study from Andalusia, Spain. Catena 2015, 132, 1–11. [Google Scholar] [CrossRef]
  35. Snavely, N.; Seitz, S.M.; Szeliski, R. Photo tourism: Exploring photo collections in 3D. ACM Trans. Graph. 2006, 25, 835–846. [Google Scholar] [CrossRef]
  36. Furukawa, Y.; Ponce, J. Accurate, dense, and robust multiview stereopsis. IEEE Trans. Pattern Anal. Mach. Intell. 2007, 32, 1362–1376. [Google Scholar] [CrossRef] [PubMed]
  37. Castillo, C.; James, M.R.; Redelmacías, M.D.; Pérez, R.; Gómez, J.A. SF3M software: 3-D photo-reconstruction for non-expert users and its application to a gully network. Soil 2015, 1, 583–594. [Google Scholar] [CrossRef]
  38. Gomez-Gutierrez, A.; Schnabel, S.; Berenguer-Sempere, F.; Lavado-Contador, F.; Rubio-Delgado, J. Using 3D photo-reconstruction methods to estimate gully headcut erosion. Catena 2014, 120, 91–101. [Google Scholar] [CrossRef]
  39. Eltner, A.; Baumgart, P.; Maas, H.G.; Faust, D. Multi-temporal UAV data for automatic measurement of rill and interrill erosion on loess soil. Earth Surf. Process. Landf. 2015, 40, 741–755. [Google Scholar] [CrossRef]
  40. Poesen, J. Challenges in gully erosion research. Landf. Anal. 2011, 17, 5–9. [Google Scholar]
  41. Smith, M.W.; Vericat, D. From experimental plots to experimental landscapes: Topography, erosion and deposition in sub-humid badlands from structure-from-motion photogrammetry. Earth Surf. Process. Landf. 2015, 40, 1656–1671. [Google Scholar] [CrossRef]
  42. King, C.; Baghdadi, N.; Lecomte, V.; Cerdan, O. The application of remote-sensing data to monitoring and modelling of soil erosion. Catena 2005, 62, 79–93. [Google Scholar] [CrossRef]
  43. D’Oleire-Oltmanns, S.; Marzolff, I.; Tiede, D.; Blaschke, T. Detection of gully-affected areas by applying object-based image analysis (OBIA) in the region of Taroudannt, Morocco. Remote Sens.-Basel 2014, 6, 8287–8309. [Google Scholar] [CrossRef]
  44. Eltner, A.; Kaiser, A.; Castillo, C.; Rock, G.; Neugirg, F.; Abellan, A. Image-based surface reconstruction in geomorphometry–merits, limits and developments of a promising tool for geoscientists. Image 2015, 3, 1445–1508. [Google Scholar] [CrossRef]
  45. Wang, Q.; Xu, J.; Chen, Y.; Li, J.; Wang, X. Influence of the varied spatial resolution of remote sensing images on urban and rural residential information extraction. Resour. Sci. 2012, 34, 159–165. (In Chinese) [Google Scholar]
  46. Shruthi, R.B.V.; Kerle, N.; Jetten, V.; Abdellah, L.; Machmach, I. Quantifying temporal changes in gully erosion areas with object oriented analysis. Catena 2015, 128, 262–277. [Google Scholar] [CrossRef]
  47. Wang, T.; He, F.H.; Zhang, A.D.; Gu, L.J.; Wen, Y.M.; Jiang, W.G.; Shao, H.B. A quantitative study of gully erosion based on object-oriented analysis techniques: A case study in Beiyanzikou catchment of Qixia, Shandong, China. Sci. World J. 2014, 2014, 417325. [Google Scholar] [CrossRef] [PubMed]
  48. Hughes, A.O.; Prosser, I.P. Gully erosion prediction across a large region: Murray-Darling Basin, Australia. Soil Res. 2012, 50, 267–277. [Google Scholar] [CrossRef]
  49. Mousazadeh, F.; Salleh, K.O. The influence of lithology and soil on the occurrence and expansion of gully erosion, Toroud Basin—Iran. Procedia Soc. Behav. Sci. 2014, 120, 749–756. [Google Scholar] [CrossRef]
  50. Mousazadeh, F.; Salleh, K.O. Factors controlling gully erosion development in Toroud Basin—Iran. Procedia Soc. Behav. Sci. 2014, 120, 506–512. [Google Scholar] [CrossRef]
  51. Ranzi, R.; Le, T.H.; Rulli, M.C. A RUSLE approach to model suspended sediment load in the Lo river (Vietnam): Effects of reservoirs and land use changes. J. Hydrol. 2012, 422–423, 17–29. [Google Scholar] [CrossRef]
  52. Wang, W.; Zhang, S.; Fang, H. Coupling mechanism of slope-gully erosion in typical black soil area of Northeast China. J. Nat. Resour. 2012, 27, 2113–2122. (In Chinese) [Google Scholar]
  53. Mancini, F.; Dubbini, M.; Gattelli, M.; Stecchi, F.; Fabbri, S.; Gabbianelli, G. Using unmanned aerial vehicles (UAV) for high-resolution reconstruction of topography: The structure from motion approach on coastal environments. Remote Sens.-Basel 2013, 5, 6880–6898. [Google Scholar] [CrossRef] [Green Version]
  54. Hofle, B.; Griesbaum, L.; Forbriger, M. GIS-based detection of gullies in terrestrial LiDAR data of the Cerro llamoca peatland (Peru). Remote Sens.-Basel 2013, 5, 5851–5870. [Google Scholar] [CrossRef]
  55. Kociuba, W.; Janicki, G.; Rodzik, J.; Stepniewski, K. Comparison of volumetric and remote sensing methods (TLS) for assessing the development of a permanent forested loess gully. Nat. Hazards 2015, 79, S139–S158. [Google Scholar] [CrossRef]
  56. Harwin, S.; Lucieer, A. Assessing the accuracy of georeferenced point clouds produced via multi-view stereopsis from unmanned aerial vehicle (UAV) imagery. Remote Sens.-Basel 2012, 4, 1573–1599. [Google Scholar] [CrossRef]
  57. Evans, M.; Lindsay, J. High resolution quantification of gully erosion in upland peatlands at the landscape scale. Earth Surf. Process. Landf. 2010, 35, 876–886. [Google Scholar] [CrossRef]
  58. Wen, L.; Zheng, F.; Shen, H.; Gao, Y. Effects of corn straw mulch buffer in the gully head on gully erosion of sloping cropland in the black soil region of Northeast China. J. Sediment. Res. 2014, 6, 73–80. (In Chinese) [Google Scholar]
Figure 1. Location of the study area: (a) map of Northeast China; (b) map of the study area, situated at the southeastern edge of the black soil region; (c) overview of the study area with a detailed view of the survey site—the detailed RGB-321 Pleiades-1A satellite image includes the Sancha River catchment (12 October 2013); and (d) the investigated site was used for detailed mapping and the Ground Control Points (GCPs) were distributed over the site.
Figure 1. Location of the study area: (a) map of Northeast China; (b) map of the study area, situated at the southeastern edge of the black soil region; (c) overview of the study area with a detailed view of the survey site—the detailed RGB-321 Pleiades-1A satellite image includes the Sancha River catchment (12 October 2013); and (d) the investigated site was used for detailed mapping and the Ground Control Points (GCPs) were distributed over the site.
Ijgi 05 00200 g001
Figure 2. The distribution map of: (a) erosion gullies; (b) gully density; and (c) gully intensity in the study area.
Figure 2. The distribution map of: (a) erosion gullies; (b) gully density; and (c) gully intensity in the study area.
Ijgi 05 00200 g002aIjgi 05 00200 g002b
Figure 3. Orthorectified image mosaic of the small study site (31 October 2015). Darker areas within the gully system are withered vegetation. The main gully in the northern section is Gully-4, as mentioned above.
Figure 3. Orthorectified image mosaic of the small study site (31 October 2015). Darker areas within the gully system are withered vegetation. The main gully in the northern section is Gully-4, as mentioned above.
Ijgi 05 00200 g003
Figure 4. Digital Surface Model (DSM) of the small study site. It was applied by hill-shading.
Figure 4. Digital Surface Model (DSM) of the small study site. It was applied by hill-shading.
Ijgi 05 00200 g004
Figure 5. (ac). 3D model of the small study site (c). The blue (a); and red (b) enlarged insets show the 3D structure of the gully sidewall and headcut in detail, respectively.
Figure 5. (ac). 3D model of the small study site (c). The blue (a); and red (b) enlarged insets show the 3D structure of the gully sidewall and headcut in detail, respectively.
Ijgi 05 00200 g005
Figure 6. Topographic differentiation of gully erosion: (a) elevation; (b) slope; and (c) slope aspect.
Figure 6. Topographic differentiation of gully erosion: (a) elevation; (b) slope; and (c) slope aspect.
Ijgi 05 00200 g006aIjgi 05 00200 g006b
Figure 7. Comparison of the satellite image (12 October 2013; left) and UAV image (31 October 2015; right).
Figure 7. Comparison of the satellite image (12 October 2013; left) and UAV image (31 October 2015; right).
Ijgi 05 00200 g007
Table 1. Interpretation signs: four types of gully appearance in the high-resolution satellite images and ground photographs of the study area.
Table 1. Interpretation signs: four types of gully appearance in the high-resolution satellite images and ground photographs of the study area.
Types of Gully AppearanceGround PhotosPleiades Images
Gully-1: Ephemeral gully Ijgi 05 00200 i001 Ijgi 05 00200 i002
Gully-2: Small permanent gully Ijgi 05 00200 i003 Ijgi 05 00200 i004
Gully-3: Large permanent gully Ijgi 05 00200 i005 Ijgi 05 00200 i006
Gully-4: Ramified permanent gully Ijgi 05 00200 i007 Ijgi 05 00200 i008
Table 2. Status of gully erosion in the Sancha River catchment.
Table 2. Status of gully erosion in the Sancha River catchment.
Numbers of GulliesThe Total Length (km)The Total Area (hm2)Gully Density of Whole Area (km/km2)Gully Intensity of Whole Area (m2/km2)
1190224.97889.0391.375434.84
Table 3. Field validation results.
Table 3. Field validation results.
Validated GulliesCorrectFalseAccuracy of Visual Interpretation
2018290%
Table 4. Summary of the outputs obtained in Pix4dmapper.
Table 4. Summary of the outputs obtained in Pix4dmapper.
Imagery and ResolutionGeoreferencingBundle Block Adjustment
Images calibratedArea covered (ha)GSD (cm)Number of GCPsRMSE (X/Y/Z) in (m)Number of 2D keypoint observationsNumber of 3D pointsMean reprojection error (pixels)
25513.854.42400.026/0.016/0.0242,622,204899,0370.26

Share and Cite

MDPI and ACS Style

Wang, R.; Zhang, S.; Pu, L.; Yang, J.; Yang, C.; Chen, J.; Guan, C.; Wang, Q.; Chen, D.; Fu, B.; et al. Gully Erosion Mapping and Monitoring at Multiple Scales Based on Multi-Source Remote Sensing Data of the Sancha River Catchment, Northeast China. ISPRS Int. J. Geo-Inf. 2016, 5, 200. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi5110200

AMA Style

Wang R, Zhang S, Pu L, Yang J, Yang C, Chen J, Guan C, Wang Q, Chen D, Fu B, et al. Gully Erosion Mapping and Monitoring at Multiple Scales Based on Multi-Source Remote Sensing Data of the Sancha River Catchment, Northeast China. ISPRS International Journal of Geo-Information. 2016; 5(11):200. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi5110200

Chicago/Turabian Style

Wang, Ranghu, Shuwen Zhang, Luoman Pu, Jiuchun Yang, Chaobin Yang, Jing Chen, Cong Guan, Qing Wang, Dan Chen, Bolin Fu, and et al. 2016. "Gully Erosion Mapping and Monitoring at Multiple Scales Based on Multi-Source Remote Sensing Data of the Sancha River Catchment, Northeast China" ISPRS International Journal of Geo-Information 5, no. 11: 200. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi5110200

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