Mapping soil texture classes using field texturing, particle size distribution and local knowledge by both conventional and geostatistical methods

We investigated the utility of three interpolation techniques that ignored descriptive 'soft' information and one that used it for mapping topsoil tex ture classes: re-coding of soil map units within Geographical Information S ystems (GIS), Thiessen polygons, and classification of probability Vectors estimated by ordinary indicator kriging and simple indicator kriging with l ocal prior means. The results were compared with texture maps based on a cl assification of kriged maps of particle size distribution. The methods were applied to two distinct regions, which represent large areas in rain-fed r ice ecosystems and irrigated rice ecosystems. The 'hard' databases for both areas contained soil information needed for mapping at regional scales (1: 100 000-1:150 000). These data were complemented with 'soft' information de rived from farmers and soil experts (Northeast Thailand) and soil maps (Nue va Ecija, Philippines). Interpolated maps agreed with the texture map based on interpolation of par ticle size distribution, and field estimates of soil texture proved to be v aluable surrogates for laboratory measurements of soil texture classes. The interpolation of categorical data such as soil texture classes allows for upgrading and increasing the resolution of maps in sparsely sampled regions by using simple held measurements. Validation using independent test sets showed that indicator kriging with l ocal prior means performed best in the rain-fed lands, whereas soft informa tion did not improve the predictions in Nueva Ecija. Local knowledge in a f ormalized form was valuable in Northeast Thailand and the interpolated soil texture map for this area had an accuracy and resolution to support agrono mic decisions at the Village scale. The poor quality of the soil map and th e fact that the gradually changing variability in young alluvial soils can be modelled with fewer data explained the lower accuracy of simple indicato r kriging with local prior means in Nueva Ecija. Thiessen polygons performe d well in the undulating rain-fed lands but were not as reliable as indicat or kriging in the gradually changing irrigated lands.