Journal of Geographical Sciences ›› 2015, Vol. 25 ›› Issue (2): 196210.doi: 10.1007/s1144201511622
• Orginal Article • Previous Articles Next Articles
BERGONSE Rafaello(), REIS Eusébio
Received:
20130605
Accepted:
20140420
Online:
20150215
Published:
20150624
About author:
Author: Rafaello Bergonse, Email:
BERGONSE Rafaello, REIS Eusébio. Reconstructing preerosion topography using spatial interpolation techniques: A validationbased approach[J].Journal of Geographical Sciences, 2015, 25(2): 196210.
Table 1
Contexts and methodologies of some published preerosion surface reconstructions"
Author  Purpose of reconstruction  Interpolation method 

Wells and Gutiérrez (1982)  Estimation of eroded volumes and combination of results with current mean erosion rates in order to date badland initiation  Undefined 
Daba et al. (2003)  Estimation of eroded volume in a large gully system for two different dates; comparison of results in order to quantify temporal evolution  Undefined^{1} 
Alexander et al. (2008)  Understanding the geomorphic evolution of a badland site from a set of remnant surfaces, and estimating rates of denudation  Linear interpolation (Triangulated Irregular Networks) 
Perroy et al. (2010)  Estimating volumetric soil loss from a set of gully channels  Linear interpolation (grid based) 
Buccolini et al. (2012)  Estimating the volume eroded by a set of gully systems (calanchi), and relating preerosion topography to gully system properties  Linear interpolation (manual^{2}) 
Table 2
General properties of the most common commercially available exact spatial interpolation methods"
Method  General features  Smoothing  Proximity  Geostatistical assumptions 

Linear interpolation  May be based on a previous Delaunay triangulation, with the value for each cell being defined by the linear surface of the triangle it overlays (e.g. Surfer 10^{1}, ArcGIS 9.1). In other cases, estimations are obtained simply as a function of the nearest known values and the respective distances (e.g. IDRISI Andes^{2}: Eastman, 2006)  None  Local  No 
Inverse Distance Weighted  Interpolated values are a function of the values of the nearest points (quantity is userdefined), with the weight of each in the result being a function of distance.  None  Local to Global  No 
Splines  Generated surface results from fitting a polynomial to a quantity of userdefined known values, subjected to two constraints: (1) surface passes exactly through the known data points; (2) curvature of generated surface is minimized. Has problems representing discrete transitions (e.g. limits of flood plains, slope breaks), sometimes ‘overshooting’ the true surface (  Elevated  Local to Global  No 
Topo to Raster  Similar to Spline, but modified in order to produce a hydrologically correct surface and incorporate slope breaks. Conceived to use points, lines and polygons as input.  Elevated  Local to Global^{3}  No 
Ordinary Kriging  Based on preliminary analysis and statistical modelling of the variation of differences between all known values with spatial distance and/or direction. For each location, the functions thus defined are used to estimate values from surrounding data points of known value.  Medium  Local to Global  Yes 
Natural neighbour  Based on the construction of a network of Voronoy polygons incorporating all known data points. Each point to be estimated is inserted on the network, and the latter is modified in order to incorporate it. Each estimated value is the average of all known surrounding points of known value, weighted by the proportion of the new Voronoy polygon overlaying each of the initial polygons.  None  Local  No 
Table 3
The interpolation methods and parameterizations adopted"
Method (different parameter sets used)  Parameters 

Linear interpolation (1)  Obtained through triangulation and conversion of a TIN model (points as input) 
Topo to Raster (2)  Two parameterizations: points as input and contours as input. Further parameters were set as default. 
Spline (10)  Spline Regularized: w = 0; 0.001; 0.01; 0.1; 0.5 
Spline Tension: w = 0, 1, 4, 7, 10. 
Table 4
Characteristics of the distributions of absolute error obtained for each interpolation method (i.e. square root of the square of the difference between real and interpolated values). Methods are ordered by ascending mean absolute error (MAE). Spline Reg and Spline Ten respectively identify the Regularized and Tension methods; w = weight parameter; P50 and P80 are the 50th and 80th percentiles; SD  standard deviation. All values are in metres."
Method  MAE  Min  Max  P50  P80  SD 

Topo to Raster (contours)  0.752  0.000  5.399  0.440  1.203  0.872 
Spline Reg w=0.01  0.767  0.000  5.001  0.463  1.176  0.889 
Spline Reg w=0.1  0.771  0.000  5.395  0.470  1.218  0.913 
Spline Reg w=0.001  0.810  0.000  5.338  0.493  1.239  0.904 
Spline Reg w=0.5  0.813  0.000  5.827  0.473  1.242  1.000 
Spline Reg w=0  0.834  0.000  5.456  0.526  1.288  0.912 
Spline Ten w=1  0.887  0.000  5.375  0.540  1.496  0.939 
Spline Ten w=4  0.954  0.000  5.354  0.587  1.654  0.976 
Spline Ten w=7  0.998  0.000  5.349  0.619  1.758  1.004 
Spline Ten w=10  1.034  0.000  5.350  0.639  1.828  1.028 
Linear  1.214  0.000  6.908  0.815  1.962  1.239 
Topo to Raster (points)  1.589  0.000  12.773  1.094  2.383  1.795 
Spline Ten w=0  3.463  0.001  72.642  1.084  3.494  7.500 
Table 5
Characteristics of the distributions of absolute error obtained during parameter optimization. Methods are ordered by ascending mean absolute error (MAE).; w = weight parameter in the regularized (Reg) Spline method; R = roughness penalty in the Topo to Raster method; P50 and P80 are the 50th and 80th percentiles; SD  standard deviation. All values are in metres."
Method  Mean  Min  Max  P50  P80  SD 

Spline Reg w=0.033  0.758  0.001  5.206  0.458  1.160  0.890 
Spline Reg w=0.055  0.762  0.001  5.291  0.468  1.191  0.896 
Spline Reg w=0.078  0.766  0.000  5.349  0.466  1.190  0.905 
Spline Reg w=0.008  0.771  0.000  5.042  0.463  1.171  0.890 
Spline Reg w=0.006  0.776  0.000  5.095  0.466  1.175  0.892 
Spline Reg w=0.003  0.791  0.000  5.210  0.480  1.191  0.897 
Topo to Raster R=0.4  0.825  0.000  5.499  0.496  1.318  0.922 
Topo to Raster R=0.3  0.871  0.000  4.920  0.547  1.442  0.923 
Topo to Raster R=0.2  0.938  0.000  5.585  0.618  1.519  0.980 
Topo to Raster R=0.1  0.987  0.000  5.634  0.664  1.599  1.007 
Topo to Raster R=0.5  1.020  0.000  5.716  0.701  1.635  1.029 
1  Alexander R W, CalvoCases A, ArnauRosalén Eet al., 2008. Erosion and stabilization sequences in relation to base level changes in the El Cautivo badlands, SE Spain.Geomorphology, 100: 8390. 
2  Anderson M G, Burt T P, 1978. The role of topography in controlling through flow generation.Earth Surface Processes, 3: 331344. 
3  Bergonse R, Reis E, 2011. Theoretical constraints to gully erosion research: Time for a reevaluation of concepts and assumptions?Earth Surface Processes and Landforms, 36: 15541557. 
4  Betts H D, Trustrum N A, De Rose R C, 2003. Geomorphic changes in a complex gully system measured from sequential digital elevation models, and implications for management.Earth Surface Processes and Landforms, 28: 10431058. 
5  Bou Kheir R, Wilson J, Deng Y, 2007. Use of terrain variables for mapping gully erosion susceptibility in Lebanon.Earth Surface Processes and Landforms, 32: 17701782. 
6  Buccolini M, Coco L, 2010. The role of the hillside in determining the morphometric characteristics of “calanchi”: The example of Adriatic central Italy.Geomorphology, 123: 200210. 
7  Buccolini M, Coco L, Cappadonia Cet al., 2012. Relationships between a new slope morphometric index and calanchi erosion in northern Sicily, Italy.Geomorphology, 149/150: 4148. 
8  Daba S, Rieger W, Strauss P, 2003. Assessment of gully erosion in eastern Ethiopia using photogrammetric techniques.Catena, 50: 273291. 
9  De Rose R C, Gomez B, Marden Met al., 1998. Gully erosion in Mangatu Forest, New Zealand, estimated from digital elevation models.Earth Surface Processes and Landforms, 23: 10451053. 
10  Dunne T, Moore T R, Taylor C H, 1975. Recognition and prediction of runoffproducing zones in humid regions.Hydrological Sciences Bulletin, XX(3): 305327. 
11  Eastman J R, 2006. IDRISI Andres Guide to GIS and Image Processing. Clark University. 
12  ESRI (Environmental Systems Research Institute), 2005 
13  Hengl T, Evans I S, 2009. Mathematical and digital models of the land surface. In: Hengl T, Reuter H I. Geomorphometry: Concepts, Software, Applications. Elsevier, 3163. 
14  Heywood I, Cornelius S, Carver S, 2006. An Introduction to Geographic Information Systems. 3rd ed. Pearson Education. 
15  Hutchinson M F, 1988. Calculation of hydrologically sound digital elevation models. Paper presented at the Third International Symposium on Spatial Data Handling. Sydney, Australia. 
16  Hutchinson M F, 1989. A new procedure for gridding elevation and stream line data with automatic removal of spurious pits.Journal of Hydrology, 106: 211232. 
17  Kakembo V, Xanga W W, Rowntree K, 2009. Topographic thresholds in gully development on the hillslopes of communal areas in Ngqushwa Local Municipality, Eastern Cape, South Africa.Geomorphology, 110: 188194. 
18  MartínezCasasnovas J A, Ramos M C, GarciaHernández D, 2009. Effects of landuse changes in vegetation cover and sidewall erosion in a gully head of the Penedès region (northeast Spain).Earth Surface Processes and Landforms, 34: 19271937. 
19  Neteler M, Mitasova H, 2008. Open Source GIS: A GRASS GIS Approach. Springer. 
20  Olaya V, 2009. Basic landsurface parameters. In: Hengl T, Reuter H (eds.). Geomorphometry: Concepts, Software, Applications. Elsevier, 141169. 
21  Parkner T, Page M, Marutani Tet al., 2006. Development and controlling factors of gullies and gully complexes. East Coast, New Zealand.Earth Surface Processes and Landforms, 31: 187199. 
22  Perroy R L, Bookhagen B, Asner G Pet al., 2010. Comparison of gully erosion estimates using airborne and ground based LiDAR on Santa Cruz Island, California.Geomorphology, 118: 288300. 
23  Poesen J, Nachtergaele J, Verstraeten Get al., 2003. Gully erosion and environmental change: Importance and research needs.Catena, 50: 91133. 
24  Reuter H I, Nelson A, 2009. Geomorphometry in ESRI packages. In: Hengl T, Reuter H I (eds.). Geomorphometry: Concepts, Software, Applications. Elsevier, 270291. 
25  Summerfield M, 1991. Global Geomorphology. Pearson Education. 
26  Weibel R, Heller M, 1991. Digital terrain modelling. In: Maguire D J, Goodchild M F, Rhind D W (eds.). Geographical Information Systems: Principles and Applications: 269297. 
27  Wells N, Andriamihaja B, Rakotovololona H F, 1991. Patterns of development of lavaka, Madagascar’s unusual gullies.Earth Surface Processes and Landforms, 16: 189206. 
28  Wells S G, Gutierrez A A, 1982. Quaternary evolution of badlands in the southeastern Colorado plateau, USA. In: Bryan R, Yair A (eds.). Badland Geomorphology and Piping. GeoBooks. 
29  Wood A L, Simon A, Downs P Wet al., 2001. Banktoe processes in incised channels: the role of apparent cohesion in the entrainment of failed bank materials.Hydrological Processes, 15: 3961. 
30  Yao C, Lei T, Elliot W Jet al., 2008. Critical conditions for rill initiation.Transactions of the ASABE, 51(1): 107114. 
31  Yue Tianxiang, Chen Chuanfa, Li Bailian, 2010a. An adaptive method of high accuracy surface modeling and its application to simulating elevation surfaces.Transactions in GIS, 14(5): 615630. 
32  Yue Tianxiang, Chen Chuanfa, Li Bailian, 2012. A highaccuracy method for filling voids and its verification.International Journal of Remote Sensing, 33(9): 28152830. 
33  Yue Tianxiang, Song Dunjiang, Du Zhengpinget al., 2010b. Highaccuracy surface modelling and its application to DEM generation.International Journal of Remote Sensing, 31(8): 22052226. 
34  Zucca C, Canu A, Della Peruta R, 2006. Effects of land use and landscape on spatial distribution and morphological features of gullies in an agropastoral area in Sardinia (Italy).Catena, 68: 8795. 
[1]  YAO Yonghui, ZHANG Baiping. MODISbased estimation of air temperature of the Tibetan Plateau [J]. , 2013, 23(4): 627640. 
[2]  SHI Wenjiao, LIU Jiyuan, DU Zhengping, YUE Tianxiang. Development of a surface modeling method for mapping soil properties [J]. Journal of Geographical Sciences, 2012, 22(4): 752760. 
[3]  ZHU Guofeng, HE Yuanqing, PU Tao, WANG Xufeng, JIA Wenxiong, LI Zongxing, XIN Huijuan. Spatial distribution and temporal trends in potential evapotranspiration over Hengduan Mountains region from 1960 to 2009 [J]. Journal of Geographical Sciences, 2012, 22(1): 7185 . 
[4]  HU Gang, WU Yongqiu, LIU Baoyuan, ZHANG Yongguang, YOU Zhimin,YU Zhangtao. The characteristics of gully erosion over rolling hilly black soil areas of Northeast China [J]. Journal of Geographical Sciences, 2009, 19(3): 309320. 
[5]  LI Xiaoyan, ZHANG Shuwen, WANG Zongming, ZHANG Huilin. Spatial variability and pattern analysis of soil properties in Dehui city, Jilin province [J]. Journal of Geographical Sciences, 2004, 14(4): 503511. 
