HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Free Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Li, W. Articles by Zhang, C. Search for Related Content PubMed Articles by Li, W. Articles by Zhang, C. GeoRef GeoRef Citation Agricola Articles by Li, W. Articles by Zhang, C. Related Collections Geostatistics Stochastic Processes Predictive soil mapping

SOIL PHYSICS

A Random-Path Markov Chain Algorithm for Simulating Categorical Soil Variables from Random Point Samples

Dep. of Geography, Kent State Univ., Kent, OH 44242

^* Corresponding author (weidong6616{at}yahoo.com).

Quantitative prediction and simulation of categorical soil variables from limited samples are crucial for cost-effectively acquiring exhaustive area-class maps. Conventional methods, however, usually cannot meet all of the requirements for class simulation in incorporating interclass relationships and generating polygon features. We developed a random-path sequential simulation algorithm based on the Markov chain random field theory. Our objective was to find a suitable method for predictive area-class soil mapping from irregularly distributed point samples, and thus extend Markov chains into practical nonlinear geostatistics. The algorithm was used to simulate soil type maps conditioned on three different sample data sets, and was compared with the widely used indicator kriging simulation algorithm–sequential indicator simulation with ordinary indicator kriging (SISoik). Results show that the algorithm works well with both dense and sparse random samples in reproducing all classes and input statistics. Compared with SISoik, the algorithm has the following advantages: (i) it more effectively captures complex patterns of soil classes and obeys their interclass relationships; (ii) it generates lower spatial uncertainty and more accurate realizations—for example, the relative increases in average percentage of correctly classified locations of realizations for the sparse, medium, and dense data sets are 5.0, 9.9, and 8.5%, respectively; (iii) it generates polygon features in realizations in accordance with the style of area-class soil maps; and (iv) it can generate classes missed in sampling but confirmed by experts. We concluded that the algorithm provides a practical spatial statistical tool for prediction and simulation of categorical soil spatial variables.

Abbreviations: AMP, average maximum probability • CPD, conditional probability distribution • IK, indicator kriging • MCG, Markov chain geostatistics • MCRF, Markov chain random field • MCSS, Markov chain sequential simulation • OIK, ordinary indicator kriging • PCC, percentage of correctly classified locations • SIK, simple indicator kriging • SIS, sequential indicator simulation • SISoik, sequential indicator simulation with ordinary kriging • TPM, transition probability matrix