Open access

Terrain Analysis for Locating Erosion Channels: Assessing LiDAR Data and Flow Direction Algorithm

Written By

Adam Pike, Tom Mueller, Eduardo Rienzi, Surendran Neelakantan, Blazan Mijatovic, Tasos Karathanasis and Marcos Rodrigues

Submitted: 10 February 2012 Published: 21 November 2012

DOI: 10.5772/51526

From the Edited Volume

Research on Soil Erosion

Edited by Danilo Godone, Silvia Stanchi

Chapter metrics overview

2,978 Chapter Downloads

View Full Metrics

1. Introduction

Terrain analysis can be used to locate concentrated flow erosion (e.g., ephemeral gully erosion) across landscapes. For example, studies have found that ephemeral gullies were likely to occur when field specific thresholds were exceeded for the following terrain attributes: the product of upslope area, slope, and plan curvatures [1]; topographic wetness index, upslope area, and slope [2]; and the topographic wetness index and the product of the upslope area and slope [3]. Another study used a cartographic classification and threshold procedure for erosion channel identification [4].

An alternative approach utilized logistic regression and artificial neural network procedures to predict where erosion channels would appear in agricultural fields based on digital terrain attributes [5]. With leave-one-field-out validation, it was determined that the more simple logistic regression was more appropriate because it performed as well as the non-linear neural network procedure. In a follow up study, erosion channels predicted from terrain attributes derived from 10-m US Geological Survey (USGS) digital elevation models (DEMs) were compared to those derived from DEMs created with survey-grade real-time kinematic (RTK) Global Positioning System (GPS) data [6]. The USGS models identified most eroded features but the RTK analyses delineated them more clearly. The authors concluded that the USGS predictions were adequate for many agricultural applications because creating DEMs with RTK was relatively costly while USGS data was freely available on the Internet for most of the United States. A graphical representation illustrates how the logistic regression analysis [5,6] can be fit (step 1) and then applied (step 2) in Figure 1.

In another study, the model from reference [5] predicted where eroded waterways would occur 223 km away in a different physiograpic region of Kentucky where there were marked textural and parent material differences between the soils [7]. The locations of the Outer Bluegrass [6,7] and Western Coal Fields studies are shown in Figure 2.

1.1. RTK GPS, LiDAR, and USGS DEM Accuracy

The performance of terrain analysis applications depends on the source and quality of the elevation information. Very accurate elevation measurements can be obtained with survey-grade RTK GPS (e.g., horizontal rmse < 2.2 cm).

Figure 1.

Description of the basic two-step logistic regression analysis procedure used in references [5,6]. The first step involves fitting a model with field observations of soil erosion as a function of terrain attributes. The second step involves applying the model. In this example, the model was applied to the same field as it was fit. However, references [5,6] used a leave one-field out validation and fit the model to four fields and applied it in a fifth. See equation 3 in methods section for a mathematical definition of the Length-Slope terrain attribute.

Figure 2.

Locations of the Pike et al. (2009 and 2010) [5,6] and Luck et al. (2010) [7] studies. The Kentucky Physiographic regions are labelled and colored.

This usually involves setting up one RTK receiver as a base station, to broadcast via radio corrections to a mobile receiver, which calculates elevation on-the-fly. Elevation data can also be obtained with light detecting and ranging (LiDAR) systems mounted at the bottom of airplanes with a pulsing laser scanning rapidly from side to side. Detectors determine the time required for laser pulses to bounce back in order to calculate the distance between the aircraft and the ground surface. With simultaneous in-flight RTK GPS measurements, estimates of ground elevation are obtained at a high spatial intensity (e.g., several points per square meter). In one study [8], RMSEs for LIDAR were determined to be 33.3 cm in short grass. The USGS National Elevation Dataset (NED) includes Level-2 DEMs for most of the United States, which in many cases were reinterpolated from USGS topographic contours. These contour maps were originally created by the USGS using stereo orthophotogrammetry. Level-2 DEMs have an accuracy of one-half the contour width (e.g., the contour width was for the studies in [5,6] was 304 cm).

Of the three elevation datasources discussed, RTK GPS is one of the most accurate methods for creating ground surveys. Unfortunately, these systems are relatively expensive and obtaining ground measurements and data processing may be labor intensive. USGS DEMs are an attractive option because they are free; however, they do not identify erosion pathways as clearly as the RTK data [6]. Elevation data obtained with LiDAR is relatively inexpensive on a per area basis and are currently being obtained by various government agencies not necessarily associated with agriculture, but urban planning, transportation, and geophysics applications. LiDAR has a great potential to aid in the identification of concentrated flow pathways [9,10].

1.2. Flow Direction Algorithm

Terrain attributes (e.g., length–slope, topographic wetness index) require estimates of the upslope contributing area for each cell in the DEM. This necessitates the calculation of single or multiple flow direction for each cell in the DEM. The most common single direction flow model is the deterministic eight-neighbor (D8) procedure. Multiple direction flow models include fractional deterministic eight-neighbor (FD8), digital elevation model networks (DEMON), and deterministic infinity (D∞). There are a number of software programs that can calculate terrain attributes, including ArcGIS, Grid-Based Terrain Analysis Programs for the Environmental Sciences (TAPES-G), and Terrain Analysis Using Digital Elevation Models (TauDEM). ArcGIS, TauDEM, and TAPES-G predict D8 single direction flow. The TAPES-G program can also estimate multidirectional flow with FD8 and DEMON while TauDEM program can estimate flow with D∞.

O’Callaghan and Mark [11] developed the D8 model that routs flow from each cell to one of its single eight neighbors in the cardinal or diagonal direction with the steepest grade. However, the D8 method does not accurately model divergent flow in ridge areas and it produces unrealistic parallel flow lines [12]. The FD8 multidirection method [13] allows flow to be routed to more than one of its eight neighbors in an amount proportional to the slope gradient between the center cell and the adjacent eight neighbors. Once the upslope contributing area for each cell exceeds a threshold (i.e., maximum cross grading area), the FD8 procedure switches to D8 flow, allowing the modeling of both divergent and convergent flow [12]. Lea [14] developed an aspect driven kinematic routing algorithm that was no longer restricted to the eight nearest neighbors (i.e., cardinal and diagonal directions). This algorithm inspired both the D∞ and DEMON methods. The DEMON stream tube method developed by Costa-Cabral and Burges [15] is computationally intensive and involves routing flow downstream along tubes that expand and contract in a way where the tubes do not necessarily coincide with the cell boundaries [12]. The D∞ algorithm first calculates flow direction from the infinite set of possible flow directions around each cell, not just in the 8 cardinal directions.

The terrain attributes in references [5-7] were calculated with TAPES-G utilizing the FD8 method. TAPES-G still exists but it is no longer being supported and updated by the developers. The Windows installation only works on legacy versions of ArcGIS. The preferred approach for terrain analysis today has become TauDEM with the D∞ flow direction algorithm.

1.3. Concentrated Flow Erosion, Grassed Waterways and The Environment

Agriculture contributes 80% of the excessive phosphorus (P) flowing to the Gulf of Mexico from the Mississippi River Basin. Among the states contributing to this environmental impact, Kentucky is the sixth largest contributor of nitrogen and phosphorus to the Mississippi River Basin [16]. Sediment is the leading cause of impairment of Kentucky’s rivers and streams, impacting over 4,329 linear km with agriculture being the primary contributing source [17]. The problem of sedimentation, often underestimated, is enormous. It not only adversely impacts aquatic life [18] and impaires ecological and economic functioning of drainage systems, wetlands, streams, and rivers but also increases the risk of flooding and the need for water treatment. Properly constructed grassed waterways (NRCS Code 412), with appropriately sized vegetative filters on either side reduce considerable runoff, nutrient, and sediment delivery to surface waters. Specifically, they reduce ephemeral gully erosion, which is a substantial but often overlooked problem that accounts for about 40% of all erosion in agricultural fields [19]. In one study, the installation of grassed waterways led to reductions of 39 and 82% in runoff and sediment, respectively [20]. In another study, grassed waterways reduced dissolved reactive phosphorus losses by factors ranging between 4 to 7 and particulate phosphorus by factors ranging between 4 and 10 [21]. Because grassed waterways are so effective, the USDA Conservation Reserve Program (CRP) and Environmental Quality Incentives Program (EQIP) provide funding for this and other conservation practices.

1.4. Research Objective

The objective of this study was to compare how predictions of concentrated flow erosion performed with LiDAR, survey grade RTK GPS, and 10-m USGS DEMs, and using the D8, FD8, DEMON, and D∞ flow direction models.


2. Materials and Methods

This study was conducted in five fields in Shelby County located in the Outer Bluegrass physiographic region of Kentucky. Field A (38°17´ N, 85°9´ W), B (38°18´ N, 85°11´ W), C (38°20´ N, 85°12´ W), D (38°20´ N, 85°11´ W), and Field E (38°20´ N, 85°14´ W), were 23, 36, 11, 57, and 33 ha in size, respectively. The fields had been in a no-till, corn (Zea mays L.), wheat (Triticum aestivum L.), and double-crop soybean [Glycine max (L.) Merr.] or corn-wheat rotation for more than 20 yr. Soils in this region developed primarily from limestone residuum overlain with pedisediment from limestone weathered materials and loess [22].

2.1. Field observations

One of the co-owners of the Worth and Dee Ellis Farms was trained by the NRCS to identify eroded zones from concentrated water flow. The farmer delineated and mapped the eroded ephemeral gullies in the five study fields using a DGPS system. Subsequently, the farmer installed grassed waterways in these areas but did not reshape them as is typically the case when installing these conservation structures. This is an important distinction because reshaping these locations would have changed the terrain attributes, making the analyses presented in this chapter difficult or impossible to interpret.

Because these field observations were conducted by the farmer, we asked the NRCS district conservationist to validate the erosion delineations in the field. The conservationist visited the fields in 2008 and examined 42% of the eroded channels delineated by the farmer in the five study fields. He determined that all the channels would have been eligible for cost support for grassed waterways through the continuous USDA Conservation Reserve Program CRP program if they did not already have waterways installed.

2.2. Acquisition of Elevation Data, DEM Creation, and Terrain Analysis

Survey grade RTK GPS equipment was used to collect elevation data from Field D in 2000 [23], Fields E and B in 2004 [24], and in Fields A and C in 2007 [5]. Dual frequency Trimble AgGPS 214 (base station) and Trimble 5800 (rover) RTK receivers were used to create surveys in Fields A, B, C, and E. The survey for Field D was created with two single-frequency Trimble 4600 receivers and elevation was determined during post processing. The surveys were created at varying intensities with elevation measurements spaced approximately every 3 (Field D) and 4 m (Fields A, B, C, and E) along the direction of travel during the survey. There were 7.5 (Field D) and 12 m (Fields A, B, C, and E) between survey passes. According to the Trimble data specifications, vertical errors for all receivers were expected to be <2.2 cm because the baselines were <1200m [25,26]. The LiDAR survey was created (courtesy of Photo Science) in November of 2009 after harvest with soybean residue and stubble on the ground. The airplane height above ground was 1.219 km and speed was 185 km hr-1. The estimated horizontal and vertical accuracies were 16 and 11 cm, respectively. The average elevation point density was 2.82 points m-2. PhotoScience converted the raw LiDAR data into a triangulated irregular network (TIN). Then they converted the TIN into a 1-m raster file using natural neighbor interpolation. After the 1-m LiDAR DEM was delivered to the University of Kentucky from PhotoScience, the grid was sampled to a 4 by 4 meters. Elevation data on 9.1-m grids for these fields were obtained from Kentucky Division of Geographic Information (KDGI) ( The USGS DEMs had been previously re-interpolated by KDGI from 10-m to 9.1 meters so that they could distribute in the Kentucky State Plane coordinate system. The 9.1-m USGS raster was converted to a point file format. The raw RTK GPS and USGS points data were converted to 4 by 4-m grids with the ArcGIS (ESRI, Redlands, CA) TOPOTORASTER command (drainage enforcement option was not used for the reasons described in reference [5]). To smooth the data, 1-m contour maps were created from the 4 by 4-m RTK, LiDAR, and USGS DEMs with the ArcGIS spatial analyst extension. Then the TOPOTORASTER command was used to convert the contours back into new 4-m DEMs.

Next, TAPESG was used to remove sinks (depressions) and calculate terrain attributes for the RTK and USGS datasets using the FD8 and DEMON flow direction algorithms. TauDEM was used to remove pits and calculate terrain attributes for the RTK, LiDAR, and USGS datasets using the D8 and D∞ flow direction algorithms. Primary terrain attributes included slope (β), plan curvature, upslope contributing area, and specific catchment area. Slope was calculated using the finite distance algorithm in TAPES and the D∞ algorithm (i.e., slope could be in any direction, not just in the cardinal and diagonal directions) with TauDEM. For the FD8 method, the maximum cross-grading area [12], was hardcoded to a value of 50,000 m2 by the authors of the TAPESG for Windows software. Each cell in the DEM had a flow width dependent on the direction of flow entering that cell from the eight neighbors. For the TauDEM D8 and D∞ procedures, flow width was assumed to be 1 grid increment in all directions. The flow width for FD8 and DEMON algorithms is described in detail in reference [12]. The specific catchment area adjusts the upslope contributing area to account for flow width and was calculated as follows:

Specific Catchment Area= Upslope Contributing Area Flow Width E1

Secondary terrain attributes (i.e., those computed from two or more primary attributes) were also determined with TAPESG. This included the topographic wetness index calculated as

Topographic wetness index=ln Specific catchment area tan β E2

and the length-slope factor calculated as

Length-Slope=1.4 Specific Catchment   22.13 0.4 sinβ   0.0896 1.3 . E3

The length-slope terrain attribute (Eq. [3]) was developed to be an estimate of the length-slope factor from the Revised Universal Soil Loss Equation (RUSLE), and index of potential erosion [27]. Unfortunately, this index does not estimate the true length-slope well on longer and steeper slopes [22] but still has adequate utility.

2.3. Statistical analyses

A new variable, “ErWater,” was created that was assigned a value of 0 if the grid point observations were not from areas with evidence of concentrated flow and a value of 1 if there was evidence. The terrain attributes for each 4 x 4 grid from all five fields (n= 99,505) were exported from ArcGIS in a database format so they could be read by SAS. The sample module in SAS Enterprise Miner 4.3 were used to random subsample the dataset with an equal number of observations from each field and areas within the fields with and without concentrated flow erosion. So after subsampling, the final dataset included 1450 observations (145 from each of the eroded areas and non eroded areas in each of the five fields studied). SAS PROC LOGISTIC was used for logistic regression with the 1450 point subset used as the input dataset. The SCORE option was used to obtain full dataset predictions with all 99,505 records. The topographic wetness index, the estimated length-slope for the Universal Soil Loss Equation, and plan curvature were included as predictor variables. The dependent variable was the newly created ErWater variable. Proc FREQUENCY was used to create confusion tables for the score datasets and the results were reported in percentages. The confusion tables presented in this chapter can be interpreted according to the guide presented in Table 1.

Actual Field Status Predicted
NE (non-eroded) E (eroded)
NE (non-eroded) Correctly classified
non-eroded areas
Incorrectly classified un-eroded areas (Type 1 errors)
E (eroded) Incorrectly classified eroded Areas (Type II errors) Correctly classified
eroded areas

Table 1.

Confusion table interpretation guide.

The scored logit data (i.e., 0’s and 1’s) were then exported from SAS into ArcGIS where the point data were converted to raster format for display in GIS.


3. Results and Discussion

3.1. Performance of LiDAR data with D∞ and D8 Flow Direction Models

Logistic regression models and statistical tests are given in Table 2. The topographic wetness index, length-slope, and plan curvature values were all highly significant in the models. The discretized output of the LiDAR D8 (Figure 3) and D∞ (Figure 4) models indicated a high correspondence between the eroded waterway boundaries and the black shaded areas (i.e., areas with probability of concentrated flow erosion > 0.5 and ≤ 1.0). The average of the combined type 1 and type 2 error rates across the five fields was 8% for D8 and 9% for D∞ (Table 2). This data can be interpreted using Table 1 provided in the methods section.

The LiDAR predictions with D8 and D∞ were excellent (Figure 3 and 4). Interestingly, the D8 approach differentiated between the eroded and non eroded areas in the lower central portion of Field A (Figure 3). In this area, an eroded feature that appears to be an island can be observed by examining the polygon boundaries. The water flowing across this island area disappeared underground and reemerged in the waterway below. The terrain modeling predicted high upslope contributing area and therefore high topographic wetness and length-slope index indices below this island area; however, plan curvature values did not have large negative values (indicating concavity) because there was no erosion area as determined during field observations with the NRCS conservationist. The D8 model differentiated this area because most of the prediction weight from this model came from plan curvature values. This was apparent from the Wald Chi Square statistic which was much greater for the D8 (180) than D∞ (71) analyses.

Flow Direction Method Variable Parameter Estimate Wald ChiSquare
D8 Intercept -3.58 68 **
Topographic Wetness Index 0.303 24 **
Length-Slope 0.758 45 **
Plan Curvature -6.81 180 **
D∞ Intercept -5.97 111 **
Topographic Wetness Index 0.573 59 **
Length-Slope 1.06 75 **
Plan Curvature -4.75 71 **

Table 2.

Logistic regression model parameters and tests for the LiDAR data analyses for the TauDEM D8 and D∞ flow direction models.

3.2. Impact of data cleaning

We considered what would have been the result if we had not contoured and rasterized the data. The average combined type-1 and type-2 error rate for D8 and D∞ was 8% for the smoothed data (Table 2) and 12% for the unsmoothed data. The logistic regression analyses did not differentiate as clearly eroded and non-eroded areas when the data were not smoothed. This is very apparent when comparing the unsmoothed analyses of Fields D and E (Figure 5) with smoothed analyses (Figure 3).

Figure 3.

Discretized erosion probability maps derived from LiDAR measurements overlain by the boundaries of the observed concentrated flow pathways. Terrain attributes were calculated with TauDEM using the D 8 flow direction algorithm.

Figure 4.

Discretized erosion probability maps derived from LiDAR measurements overlain by the boundaries of the observed concentrated flow pathways. Terrain attributes were calculated with TauDEM using the D ∞ flow direction algorithm.

Flow Direction Algorithm D8 FD8
Field Actual Field Status NE E Ne E
A NE 81 16 80 17
E 1 3 1 3
B NE 83 11 81 13
E 1 5 2 5
C NE 79 18 79 18
E 0 2 0 2
D NE 78 10 76 13
E 4 7 4 8
F NE 75 17 73 18
E 3 6 3 6

Table 3.

Confusion table for the TauDEM D8 and D∞ flow direction models. The values are given in percentages of the total observations from each field. This table can be interpreted using the guide presented in Table 1.

Figure 5.

Discretized erosion probability maps derived from unsmoothed LiDAR data overlain by field boundaries. Terrain attributes were calculated with TauDEM using the D ∞ flow direction algorithm.

One problem with LiDAR data is that it contains systematic artifacts resulting from differential plant stubble heights in the direction perpendicular to crop rows associated with the direction of movement of farm machinery. With terrain modeling, they can create artificial flow lines along pathways of farm equipment. These lines are apparent in the unsmoothed but not smoothed LiDAR slope data (Figure 6). Unfortunately, the method of smoothing presented in this chapter introduced new artifacts along the contours because the procedure involved

  1. calculating contours from DEMs and

  2. rasterizing the contours.

This smoothing is very computationally resource intensive, and potentially problematic for professional conservation planners. It may be necessary for GIS experts to smooth the LiDAR data during preprocessing and make the analyses available to planners. More computationally efficient smoothing techniques should be considered that minimize new artifacts.

Figure 6.

Comparison of unsmoothed and smoothed LiDAR slope data for Field D.

3.3. Comparison of LiDAR with RTK and USGS

The type-1 and -2 average misclassification errors for TauDEM D8 and D∞ output for LiDAR (8%) were similar in size as errors for the RTK (8%) dataset and lower than those for the USGS data (12%). The predictions were not very different for any of the datasets as shown in Figure 7. In many areas throughout the United States, LiDAR is being purchased for various applications (e.g., agriculture, soil mapping, transportation, land use planning). In those situations, it would be advantageous to use the LiDAR data for identifying eroded waterways. In areas where LiDAR is not available, USGS 10-m grids may be adequate for conservation planning [6].

3.4. Impact of Flow Direction Algorithm

The FD8 and DEMON methods were not used with the LiDAR data because TAPES no longer works with the latest version of ArcGIS (i.e., version 10.0). The regression parameters and Wald Chi Square test for the RTK (Table 4) and USGS (Table 5) analyses could be compared with those for the LiDAR models (Table 2). All of the parameters were statistically significant except for plan curvature in the USGS FD8 analyses (Table 5) which was consistent with [6]. The concentrated flow prediction models can be tested in other areas across the country with the same data source (RTK, LiDAR, and USGS) and flow direction (D8, D∞, FD8, and DEMON). The users should note that use of these may only be valid when similar smoothing techniques and DEM resolutions are used. In our case, all of our analyses were made with 4 by 4-m rasters. Grid scale is a very important factor to consider because at a small scale, there may be no relationship between land forms and curvature values; however, at a large increment, landscapes could have a profoundly large impact on an analysis similar to the one used in this study.

Flow Direction Method Variable Parameter Estimate Wald ChiSquare
D8 Intercept -3.14 57 **
Topographic Wetness Index 0.251 18 **
Length-Slope 0.681 42 **
Plan Curvature -7.979 188 **
D∞ Intercept -6.17 119 **
Topographic Wetness Index 0.603 65 **
Length-Slope 1.01 77 **
Plan Curvature -5.73 82 **
FD8 Intercept -9.88 161 **
Topographic Wetness Index 0.882 88 **
Length-Slope 1.88 137 **
Plan Curvature -5.54 39 **
DEMON Intercept -5.94 79 **
Topographic Wetness Index 8.510 34 **
Length-Slope 1.06 70 **
Plan Curvature -10.1 161 **

Table 4.

Logistic regression parameters for the RTK dataset using the D8, D∞, FD8, and DEMON flow direction models.

Figure 7.

Discretized erosion probability maps for the logistic regression analysis of the LiDAR, RTK, and USGS datasets using the D∞ flow direction algorithm.

Flow Direction Method Variable Parameter Estimate Wald ChiSquare
D8 Intercept -3.08 87 **
Topographic Wetness Index 0.351 56 **
Length-Slope 0.323 17 **
Plan Curvature -4.46 99 **
D∞ Intercept -7.19 157 **
Topographic Wetness Index 0.829 123 **
Length-Slope 0.790 70 **
Plan Curvature -2.11 19 **
FD8 Intercept -7.63 200 **
Topographic Wetness Index 0.827 147 **
Length-Slope 1.02 124 **
Plan Curvature -0.565 2 ns
DEMON Intercept -8.23 123 **
Topographic Wetness Index 0.919 96 **
Length-Slope 0.895 72 **
Plan Curvature -5.69 60 **

Table 5.

Logistic regression parameters for the USGS dataset using the D8, D∞, FD8, and DEMON flow direction models.

The average Type 1 and Type 2 error rates for the RTK dataset were ranked in the follow order (from lowest to highest): FD8 (6%), DEMON (7%), D8 (9%), and D∞(11%). The differences were smaller for the USGS dataset but procedures were ranked in a similar order: FD8 (10%), DEMON (10%), D∞(11%), D8 (12%). The map analysis for the RTK data (Figure 8) demonstrates that differences between methods were small with the FD8 model showing least noise in the southern part of the field.

Figure 8.

Discretized erosion probability maps for the logistic regression analysis of the RTK dataset comparing the D8, D∞, FD8, and DEMON flow direction algorithm.


4. Conclusion and Recommendations

The findings of this study indicate that LiDAR data can be used to clearly identify eroded features in agricultural landscapes with a level of accuracy that is similar to RTK GPS and better than USGS DEMs. It is critical that LiDAR data are smoothed prior to modeling erosion channels. Smoothing removes artifacts resulting from differences in plant residue heights perpendicular to the direction of travel by farm machinery. These differences can actually cause the terrain analysis flow models to incorrectly rout water along the direction of travel. While smoothing produced better results in this chapter, the method of smoothing with ArcGIS TopoToRaster was not efficient and cannot be used over large areas. More work is needed to determine computationally efficient smoothing algorithms for terrain analysis that minimize artifacts.

Conservation planners and GIS analysts should be able to accurately identify erosion features with the D8, D∞, FD8, and DEMON flow direction algorithm. This is important because previous work was based on TAPES G which is no longer being supported. Further TauDEM can utilize very large blocks of memory to cover extensive land areas, and can also operate on high performance computers. It is important to note that the choice flow algorithm will change the model parameters so it is important that they use the correct model. We recommend the TauDEM software program which uses the D8 and D∞ procedures because it works on 64 bit machines, allows the use of multiple core processer, and works with DEMs up to 4 GB in size.

All analyses performed in this study were based on with 4-m DEMs. Efforts are necessary to better understand the impact of the scale of terrain models on the quality of erosion model predictions. It may also be possible to expand the inference space of these models by including erosion parameters in the analyses obtained from soil surveys.


5. Nomenclature

D8, deterministic eight-neighbor; FD8, fractional deterministic eight-neighbor; DEMON, digital elevation model networks; DEM, Digital Elevation Model; TauDEM, Terrain Analysis Using Digital Elevation Models; TAPES, Grid-Based Terrain Analysis Programs for the Environmental Sciences; GIS, geographic information systems; RTK, Real Time Kinematic GIS, Geographic Information System.



We appreciate the generosity of Mike Ellis in providing us access to his grassed waterway datasets and Mike, Bob, and Jim Ellis for allowing us to conduct this research on their farm. We are also grateful for the assistance of Randall Rock, Jack Kuhn, Danny Hughes from the NRCS. We would like to thank Steve Workman from the UK college of Agriculture for support from the SB-271 Water Quality Funding. We also wish to express our appreciation to Photo Science president Mike Ritchie for providing us the complementary LiDAR data used in this chapter.


  1. 1. Thorne C. R. Zezenbergen L. W. Grissinger E. H. Murphey J. B. 1986 Ephemeral Gullies as Sources of Sediment. In: Proceeding of the 4th Federal Interagency Sedimentation Conference 24-27 March 1986, Las Vegas, United States. Washington, DC, United States: Gov. Print. Office
  2. 2. Moore I. D. Burch G. J. Mackenzie D. H. 1988 Topographic Effects on the Distribution of Surface Soil Water and the Location of Ephemeral Gullies Transactions of the ASABE 31 4 1098 1107
  3. 3. Srivastava K. P. Moore I. D. 1989 Application of Terrain Analysis to Land Resource Investigations of Small Catchments in the Caribbean. In: Proceeding of the 20th Int. Conf. of the Erosion Control Association 15-18 Feb. 1989, Vancouver, BC, Canada. Streamboat Springs: International Erosion Control Association
  4. 4. Berry J. K. Delgado J. A. Pierce F. J. Khosla R. 2005 Applying Spatial Analysis for Precision Conservation across the Landscape. Journal of Soil and Water Conservation 60 6 363 370
  5. 5. Pike A. C. Mueller T. G. Schörgendorfer A. Shearer S. A. Karathanasis A. D. 2009 Erosion Index Derived from Terrain Attributes using Logistic Regression and Neural Networks Agronomy Journal 101 5 1068 1079
  6. 6. Pike A. C. Mueller T. G. Schörgendorfer A. Luck J. D. Shearer S. A. Karathanasis A. D. 2010 Locating Eroded Waterways with United States Geologic Survey Elevation Data. Agronomy Journal 102 4 1269 1273
  7. 7. Luck J. D. Mueller T. G. Shearer S. A. Pike A. C. 2010 Grassed Waterway Planning Model Evaluated for Agricultural Fields in the Western Kentucky Coal Field Physiographic Region of Kentucky. Journal of Soil and Water Conservation 65 5 280 288
  8. 8. Hodgson M. E. Jensen J. R. Schmidt L. Schill S. Davis B. 2003 An Evaluation of LiDAR and IFSAR-Derived Digital Elevation Models in Leaf-On Conditions with USGS Level 1 and Level 2 DEMs. Remote Sensing of Environment 84 2 295 308
  9. 9. Evans M. Linsay J. 2010 High Resolution Quantification of Gully Erosion in Upland Peatlands at the Landscape Scale Earth Surface Processes and Landforms 35 8 876 866
  10. 10. Perroy R. L. Bookhagen B. Asner G. P. Chadwick O. A. 2010 Comparison of Gully Erosion Estimates Using Airborne and Ground-Based LiDAR on Santa Cruz Island, California Geomorphology 118 3-4 288 300
  11. 11. O’Callaghan J. F. Mark D. M. 1984 The Extraction of Drainage Networks from Digital Elevation Data. Computer Vision, Graphics, and Image Processing 28 3 323 344
  12. 12. Gallant J. C. Wilson J. P. 2000 Primary Topographic Attributes. In: Wilson JP., Gallant JC. (ed.) Terrain Analysis: Principles and Applications New York: John Wiley & Sons 51 85
  13. 13. Moore I. D. Gessler P. E. Nielsen G. A. Peterson G. A. 1993 Soil Attribute Prediction Using Terrain Analysis Soil Science Society of America Journal 57 2 443 452
  14. 14. Lea N. J. 1992 An Aspect-Driven Kinematic Routing Algorithm. In: Parsons AJ., Abrahams AD. (ed.) Overland Flow: Hydraulics and Erosion Mechanics London: UCL Press 393 407
  15. 15. Costa-Cabral M. C. Burges S. J. 1994 Digital Elevation Model Networks (DEMON): a Model of Flow Over Hillslopes for Computation of Contributing and Dispersal Areas Water Resources Research 30 6 1681 1692
  16. 16. Alexander R. B. Smith R. A. Schwarz G. E. Boyer E. W. Nolan J. V. Brakebill J. W. 2008 Differences in Phosphorus and Nitrogen Delivery to the Gulf of Mexico from the Mississippi River Basin Environmental Science Technology 42 3 822 830
  17. 17. EPA. 2008 Kentucky water quality assessment report. accessed 22 June 2012
  18. 18. Bilotta G. S. Brazier R. E. 2008 Understanding the Influence of Suspended Solids on Water Quality and Aquatic Biota Water Research 42 12 2849 2861
  19. 19. Bennett S. J. Casali J. Robinson K. M. Kadavy K. C. 2000 Characteristics of Actively Eroding Ephemeral Gullies in an Experimental Channel Transactions of the American Society of Agricultural Engineers 43 3 641 649
  20. 20. Fiener P. Auerswald K. 2003 Effectiveness of Grassed Waterways in Reducing Runoff and Sediment Delivery from Agricultural Watersheds. Journal of Environmental Quality 32 3 927 936
  21. 21. Fiener P. Auerswald K. 2009 Effects of Hydrodynamically Rough Grassed Waterways on Dissolved Reactive Phosphorus Loads Coming from Agricultural Watersheds Journal of Environmental Quality 38 2 548 559
  22. 22. Soil Conservation Service. 1980 Soil survey of Shelby County, KY. USDA, NRCS, Fort Worth, TX.
  23. 23. Mueller T. G. Hartsockc N. J. Stombaugh T. S. Shearerb S. A. Cornelius P. L. Barnhisel R. I. 2003 Soil Electrical Conductivity Map Variability in Limestone Soils Overlain by Loess Agronomy Journal 95 3 496 507
  24. 24. Pike A. C. Mueller T. G. Mijatovic B. Koostra B. K. Poulette M. M. Prewitt R. M. Shearer S. A. 2006 Topographic Indices: Impact of Data Source. Soil Science 171 10 800 809
  25. 25. Trimble. 1997 4800 specifications. Trimble Navigation Ltd., Sunnyvale, CA. Available online at accessed 22 June 2012
  26. 26. Trimble. 2008 Trimble 5800 GPS System Datasheet. Trimble Navigation Ltd., Sunnyvale, CA. Available online at accessed 22 June 2012
  27. 27. Moore I. D. Wilson J. P. 1992 Length-Slope Factors for the Revised Universal Soil Loss Equation: Simplified Method of Estimation Journal of Soil and Water Conservation 47 5 423 428
  28. 28. Wilson J. P. Gallant J. C. 2000 Digital Terrain Analysis. In: Wilson JP., Gallant JC. (ed.) Terrain analysis: Principles and Applications New York: John Wiley & Sons 1 28

Written By

Adam Pike, Tom Mueller, Eduardo Rienzi, Surendran Neelakantan, Blazan Mijatovic, Tasos Karathanasis and Marcos Rodrigues

Submitted: 10 February 2012 Published: 21 November 2012