SSSAJ
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 5 April 2007
Published in Soil Sci Soc Am J 71:656-668 (2007)
DOI: 10.2136/sssaj2006.0173
© 2007 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Full Text Free
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Li, W.
Right arrow Articles by Zhang, C.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Li, W.
Right arrow Articles by Zhang, C.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Li, W.
Right arrow Articles by Zhang, C.
Related Collections
Right arrow Geostatistics
Right arrow Stochastic Processes
Right arrow Predictive soil mapping

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

Weidong Li* and Chuanrong Zhang

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


Figure 1
View larger version (9K):
[in this window]
[in a new window]

 
Fig. 1. Possible neighborhoods within a search radius in a Markov chain random field if only considering nearest known neighbors in four cardinal directions: (a) four known neighbors; (b) three known neighbors; (c) two known neighbors; (d) one known neighbor; and (e) no known neighbors. The white cell represents the unknown one to be estimated and black cells represent known data locations in cardinal directions.

 

Figure 2
View larger version (9K):
[in this window]
[in a new window]

 
Fig. 2. The search circle, search sectors, and the considered neighborhoods of nearest known neighbors. The search circle is divided into four equal search sectors. Left: four nearest known neighbors are found, one from each sector. Right: two nearest known neighbors are found because the other two sectors have no known data locations.

 

Figure 3
View larger version (35K):
[in this window]
[in a new window]

 
Fig. 3. The reference soil map of seven soil classes and randomly distributed sample data sets (dense: 646 points, 2.9% of total pixels; medium: 179 points, 0.8% of total pixels; sparse: 50 points, 0.2% of total pixels). The seven soil classes (soil series) include: 1, loamy alluvium; 2, Sparta; 3, Plainfield; 4, Dakota; 5, Richwood; 6, others (e.g., water); and 7, peat and muck.

 

Figure 4
View larger version (23K):
[in this window]
[in a new window]

 
Fig. 4. Experimental transiograms headed by Soil Class 1, estimated from each of the three sample data sets: (a) from the dense data set; (b) from the medium data set; and (c) from the sparse data set. Note that estimated points for each experimental transiogram are connected as a line so that different experimental transiograms can be readily differentiated.

 

Figure 5
View larger version (21K):
[in this window]
[in a new window]

 
Fig. 5. A subset of experimental transiograms (dots) headed by Class 1, estimated from the dense data set, and interpolated models (solid lines).

 

Figure 6
View larger version (50K):
[in this window]
[in a new window]

 
Fig. 6. Simulated results by Markov chain sequential simulation, conditioned on the dense data set: (a) optimal prediction map; (b) and (c) two simulated realizations; (d) maximum occurrence probability map; (e) occurrence probability map of Class 1; (f) occurrence probability map of Class 4.

 

Figure 7
View larger version (55K):
[in this window]
[in a new window]

 
Fig. 7. Simulated results by Markov chain sequential simulation, conditioned on the medium data set: (a) optimal prediction map; (b) and (c) two simulated realizations; (d) maximum occurrence probability map; (e) occurrence probability map of Class 1; (f) occurrence probability map of Class 4.

 

Figure 8
View larger version (55K):
[in this window]
[in a new window]

 
Fig. 8. Simulated results by Markov chain sequential simulation, conditioned on the sparse data set: (a) optimal prediction map; (b) and (c) two simulated realizations; (d) maximum occurrence probability map; (e) occurrence probability map of Class 1; (f) occurrence probability map of Class 4.

 

Figure 9
View larger version (26K):
[in this window]
[in a new window]

 
Fig. 9. Proportions of classes in sample data sets and simulated realizations (averaged from 100 realizations) generated by Markov chain sequential simulation.

 

Figure 10
View larger version (32K):
[in this window]
[in a new window]

 
Fig. 10. Transiograms estimated from single realizations generated by Markov chain sequential simulation, conditioned on the dense data set.

 

Figure 11
View larger version (36K):
[in this window]
[in a new window]

 
Fig. 11. Transiograms estimated from single realizations generated by Markov chain sequential simulation, conditioned on the medium data set.

 

Figure 12
View larger version (20K):
[in this window]
[in a new window]

 
Fig. 12. The interclass relationships between Class 4 and Class 7, estimated from the dense data set. They have no or little chances to be neighbors.

 

Figure 13
View larger version (59K):
[in this window]
[in a new window]

 
Fig. 13. Simulated results by sequential indicator simulation with ordinary kriging, conditioned on the dense dataset: (a) optimal prediction map; (b) and (c) two simulated realizations; (d) maximum occurrence probability map; (e) occurrence probability map of Class 1; (f) occurrence probability map of Class 4.

 

Figure 14
View larger version (61K):
[in this window]
[in a new window]

 
Fig. 14. Simulated results by sequential indicator simulation with ordinary kriging, conditioned on the medium dataset: (a) optimal prediction map; (b) and (c) two simulated realizations; (d) maximum occurrence probability map; (e) occurrence probability map of Class 1; (f) occurrence probability map of Class 4.

 

Figure 15
View larger version (61K):
[in this window]
[in a new window]

 
Fig. 15. Simulated results by sequential indicator simulation with ordinary kriging, conditioned on the sparse dataset: (a) optimal prediction map; (b) and (c) two simulated realizations; (d) maximum occurrence probability map; (e) occurrence probability map of Class 1; (f) occurrence probability map of Class 4.

 

Figure 16
View larger version (12K):
[in this window]
[in a new window]

 
Fig. 16. Average maximum probabilities estimated from maximum occurrence probability maps generated by Markov chain sequential simulation (MCSS) and sequential indicator simulation with ordinary kriging (SISoik).

 





HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Journal of Natural Resources
and Life Sciences Education
Vadose Zone Journal
Journal of Plant Registrations Journal of
Environmental Quality
The Plant Genome
Copyright © 2007 by the Soil Science Society of America.