Premium
Test of statistical means for the extrapolation of soil depth point information using overlays of spatial environmental data and bootstrapping techniques
Author(s) -
Dahlke Helen E.,
Behrens Thorsten,
Seibert Jan,
Andersson Lotta
Publication year - 2009
Publication title -
hydrological processes
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.222
H-Index - 161
eISSN - 1099-1085
pISSN - 0885-6087
DOI - 10.1002/hyp.7413
Subject(s) - kriging , soil map , transect , digital soil mapping , mean squared error , digital elevation model , range (aeronautics) , statistics , spatial analysis , extrapolation , bootstrapping (finance) , multivariate interpolation , interpolation (computer graphics) , sampling (signal processing) , data set , mathematics , soil science , soil water , environmental science , remote sensing , geology , computer science , animation , oceanography , materials science , computer graphics (images) , filter (signal processing) , bilinear interpolation , computer vision , composite material , econometrics
Abstract Hydrological modelling depends highly on the accuracy and uncertainty of model input parameters such as soil properties. Since most of these data are field surveyed, geostatistical techniques such as kriging, classification and regression trees or more sophisticated soil‐landscape models need to be applied to interpolate point information to the area. Most of the existing interpolation techniques require a random or regular distribution of points within the study area but are not adequate to satisfactorily interpolate soil catena or transect data. The soil landscape model presented in this study is predicting soil information from transect or catena point data using a statistical mean (arithmetic, geometric and harmonic mean) to calculate the soil information based on class means of merged spatial explanatory variables. A data set of 226 soil depth measurements covering a range of 0–6·5 m was used to test the model. The point data were sampled along four transects in the Stubbetorp catchment, SE‐Sweden. We overlaid a geomorphology map (8 classes) with digital elevation model‐derived topographic index maps (2–9 classes) to estimate the range of error the model produces with changing sample size and input maps. The accuracy of the soil depth predictions was estimated with the root mean square error (RMSE) based on a testing and training data set. RMSE ranged generally between 0·73 and 0·83 m ± 0·013 m depending on the amount of classes the merged layers had, but were smallest for a map combination with a low number of classes predicted with the harmonic mean (RMSE = 0·46 m). The results show that the prediction accuracy of this method depends on the number of point values in the sample, the value range of the measured attribute and the initial correlations between point values and explanatory variables, but suggests that the model approach is in general scale invariant. Copyright © 2009 John Wiley & Sons, Ltd.