Kriging soil texture under different types of nonstationarity

Abstract A geostatistical analysis assumes some form of stationarity of the variable under study, but different types of stationarity exist and often spatial data exhibit some form of nonstationarity. However, most studies assume one type of nonstationarity and consequently apply one type of interpolation method within the study area. A study area of 8×18-km area was selected because it was expected to contain complex nonstationary conditions in soil texture. Therefore, four geostatistical interpolation methods were evaluated in their ability to account for different types of nonstationarity in the topsoil silt content: two univariate interpolation methods, ordinary kriging (OK) and universal kriging (UK), and two bivariate methods, simple kriging with varying local means (SKlm) and ordinary cokriging (OCK). A digital elevation model (DEM) was used as the exhaustive secondary information for the bivariate methods. Two kinds of nonstationary conditions were identified inside the study area: (1) a large-scale trend in both the silt content and elevation, with a strong correlation between them, and (2) a very strong local fluctuation around a mean value, representing a local nonstationarity. Consequently, different techniques were applied in different parts inside the study area: the global trend was best accounted for by OCK and UK could best account for the local nonstationarity. After combining the results of the two prediction methods, it was found that the overall estimation of the silt content was more precise than when any single method was used over the entire study area.