Register      Login
Soil Research Soil Research Society
Soil, land care and environmental research
RESEARCH ARTICLE

Simulating the vertical transition of soil textural layers in north-western China with a Markov chain model

Danfeng Li A C and Ming'an Shao B D
+ Author Affiliations
- Author Affiliations

A State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Chinese Academy of Sciences and Ministry of Water Resources, Yangling 712100, Shaanxi, P.R. China.

B Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, P.R. China.

C University of Chinese Academy of Sciences, Beijing 100049, P.R. China.

D Corresponding author. Email: shaoma@igsnrr.ac.cn

Soil Research 51(3) 182-192 https://doi.org/10.1071/SR12332
Submitted: 14 November 2012  Accepted: 27 May 2013   Published: 18 June 2013

Abstract

The heterogeneity of textures in soil profiles is important for quantifying the movement of water and solutes through soil. Soil-profile textures to a depth of 300 cm were investigated at 100 sites in a 100-km2 area in the central region of the Heihe River system, where oases coexist with widespread deserts and wetland. The probability distribution of textural-layer thickness was quantified. The vertical transition of the soil textural layers was characterised by a Markov chain–log-normal distribution (MC-LN) model based on the probability of one textural type transitioning to another. Nine types of textural layers were observed: sand, loamy sand, sandy loam, silt loam, loam, clay loam, silty clay loam, silty clay, and clay. Sand was the most frequent in the profiles, whereas silt loam and clay were rare. The layers of sand and silty clay were relatively thick, and the layers of loam and clay were relatively thin. The coefficients of variation ranged from 36–87%, indicating moderate variation in the layer thickness of each textural type. The soil profile was characterised as a log-normal distribution. A χ2 test verified the Markov characteristic and the stability of the vertical change of soil textural layers. Realisations of the soil textural profiles were generated by the MC-LN model. A Monte Carlo simulation indicated that the simulated mean layer thickness of each textural type agreed well with the corresponding field observations. Element values of the transition probability matrix of the textural layers simulated by the MC-LN model deviated <12.6% from the measured values, excluding the data from the layers of clay and silt loam. The main combinations of upper to lower textural layers in the study area were loamy sand and sand (or sandy loam), sandy loam and sand (or loamy sand and loam), loam and clay loam, clay loam (or silty clay) and silty clay loam, and silty clay loam and silty clay. The MC-LN model was able to accurately quantify the vertical changes of textures in the soil profiles. This study will aid in quantification of water and solute transport in soils with vertical heterogeneity of soil textural layers.

Additional keywords: geostatistics, Markov chain, profile textural layers, stochastic simulation, transition probability matrix.


References

Aaviksoo K (1995) Simulating vegetation dynamics and land use in a mire landscape using a Markov model. Landscape and Urban Planning 31, 129–142.
Simulating vegetation dynamics and land use in a mire landscape using a Markov model.Crossref | GoogleScholarGoogle Scholar |

Anderson TW, Goodman LA (1957) Statistical inference about Markov chains. Annals of Mathematical Statistics 28, 89–110.
Statistical inference about Markov chains.Crossref | GoogleScholarGoogle Scholar |

Atkinson PM (2001) Geographical information science: GeoComputation and nonstationarity. Progress in Physical Geography 25, 111–122.

Aurbacher J, Dabbert S (2011) Generating crop sequences in land-use models using maximum entropy and Markov chains. Agricultural Systems 104, 470–479.
Generating crop sequences in land-use models using maximum entropy and Markov chains.Crossref | GoogleScholarGoogle Scholar |

Balzter H (2000) Markov chain models for vegetation dynamics. Ecological Modelling 126, 139–154.
Markov chain models for vegetation dynamics.Crossref | GoogleScholarGoogle Scholar |

Benabdellah B, Albrecht KF, Pomaz VL, Denisenko EA, Logofet DO (2003) Markov chain models for forest successions in Erzgebirge, Germany. Ecological Modelling 159, 145–160.
Markov chain models for forest successions in Erzgebirge, Germany.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD38XpsFSksrc%3D&md5=0066444cd625d33b05f47d2f177c26faCAS |

Bogaert P (2002) Spatial prediction of categorical variables: the Bayesian maximum entropy approach. Stochastic Environmental Research and Risk Assessment 16, 425–448.
Spatial prediction of categorical variables: the Bayesian maximum entropy approach.Crossref | GoogleScholarGoogle Scholar |

Bouma J, van Alphen BJ, Stoorvogel JJ (2002) Fine tuning water quality regulations in agriculture to soil differences. Environmental Science & Policy 5, 113–120.
Fine tuning water quality regulations in agriculture to soil differences.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD38XjsFOiu7c%3D&md5=f22b75ab4a5faea811ef4d92cc522c74CAS |

Burgess TM, Webster R (1984a) Optimal sampling strategies for mapping soil types. I. Distribution of boundary spacings. Journal of Soil Science 35, 641–654.
Optimal sampling strategies for mapping soil types. I. Distribution of boundary spacings.Crossref | GoogleScholarGoogle Scholar |

Burgess TM, Webster R (1984b) Optimal sampling strategies for mapping soil types. II. Risk functions and sampling intervals. Journal of Soil Science 35, 655–665.
Optimal sampling strategies for mapping soil types. II. Risk functions and sampling intervals.Crossref | GoogleScholarGoogle Scholar |

Cheng GD, Xiao DN, Wang GX (1999) On the characteristics and building of landscape ecology in arid area. Advances in Earth Science 14, 11–15. [in Chinese with English abstract]

Chuang CW, Lin CY, Chien CH, Chou WC (2011) Application of Markov-chain model for vegetation restoration assessment at landslide areas caused by a catastrophic earthquake in Central Taiwan. Ecological Modelling 222, 835–845.
Application of Markov-chain model for vegetation restoration assessment at landslide areas caused by a catastrophic earthquake in Central Taiwan.Crossref | GoogleScholarGoogle Scholar |

D’Or D, Bogaert P (2004) Spatial prediction of categorical variables with the Bayesian maximum entropy approach: the Ooypolder case study. European Journal of Soil Science 55, 763–775.
Spatial prediction of categorical variables with the Bayesian maximum entropy approach: the Ooypolder case study.Crossref | GoogleScholarGoogle Scholar |

Deutsch CV (2006) A sequential indicator simulation program for categorical variables with point and block data: BlockSIS. Computers & Geosciences 32, 1669–1681.
A sequential indicator simulation program for categorical variables with point and block data: BlockSIS.Crossref | GoogleScholarGoogle Scholar |

FAO (1998) World reference base for soil resources. In ‘Key to the reference soil groups’. World Soil Resources Report 84. (FAO: Rome)

Goovaerts P (1996) Stochastic simulation of categorical variables using a classification algorithm and simulated annealing. Mathematical Geology 28, 909–921.
Stochastic simulation of categorical variables using a classification algorithm and simulated annealing.Crossref | GoogleScholarGoogle Scholar |

Han DL, Meng XY (1999) Recent progress of research on oasis in China. Chinese Geographical Science 9, 199–205.
Recent progress of research on oasis in China.Crossref | GoogleScholarGoogle Scholar |

Hu W, Tallon LK, Si BC (2012) Elevation of time stability indices for soil water storage upscaling. Journal of Hydrology 475, 229–241.
Elevation of time stability indices for soil water storage upscaling.Crossref | GoogleScholarGoogle Scholar |

Journel AG (1983) Nonparametric estimation of spatial distributions. Journal of the International Association for Mathematical Geology 15, 445–468.
Nonparametric estimation of spatial distributions.Crossref | GoogleScholarGoogle Scholar |

Kite GW, Kouwen N (1992) Watershed modeling using land classifications. Water Resources Research 28, 3193–3200.
Watershed modeling using land classifications.Crossref | GoogleScholarGoogle Scholar |

Koltermann CE, Gorelick SM (1996) Heterogeneity in sedimentary deposits: a review of structure-imitating, process-imitating, and descriptive approaches. Water Resources Research 32, 2617–2658.
Heterogeneity in sedimentary deposits: a review of structure-imitating, process-imitating, and descriptive approaches.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DyaK28XlvF2lsr4%3D&md5=9c8046b76099004478bb12d9662a7410CAS |

Kuusk A (1995) A Markov chain model of canopy reflectance. Agricultural and Forest Meteorology 76, 221–236.
A Markov chain model of canopy reflectance.Crossref | GoogleScholarGoogle Scholar |

Kyriakidis PC, Dungan JL (2001) A geostatistical approach for mapping thematic classification accuracy and evaluating the impact of inaccurate spatial data on ecological model predictions. Environmental and Ecological Statistics 8, 311–330.
A geostatistical approach for mapping thematic classification accuracy and evaluating the impact of inaccurate spatial data on ecological model predictions.Crossref | GoogleScholarGoogle Scholar |

Li WD (2007) Markov chain random fields for estimation of categorical variables. Mathematical Geology 39, 321–335.
Markov chain random fields for estimation of categorical variables.Crossref | GoogleScholarGoogle Scholar |

Li WD, Zhang CR (2006) A generalized Markov chain approach for conditional simulation of categorical variables from grid samples. Transactions in GIS 10, 651–669.
A generalized Markov chain approach for conditional simulation of categorical variables from grid samples.Crossref | GoogleScholarGoogle Scholar |

Li WD, Li BG, Shi YC, Tang DY (1997) Application of the Markov chain theory to describe spatial distribution of textural layers. Soil Science 162, 672–683.
Application of the Markov chain theory to describe spatial distribution of textural layers.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DyaK2sXmsV2rsLc%3D&md5=437dadf35c16e2aad2887b2f659a1038CAS |

Li WD, Li BG, Shi YC (1999) Markov-chain simulation of soil textural profiles. Geoderma 92, 37–53.
Markov-chain simulation of soil textural profiles.Crossref | GoogleScholarGoogle Scholar |

Li WD, Zhang CR, Burt JE, Zhu AX, Feyan J (2004) Two-dimensional Markov chain simulation of soil type spatial distribution. Soil Science Society of America Journal 68, 1479–1490.
Two-dimensional Markov chain simulation of soil type spatial distribution.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD2cXns1Smsb4%3D&md5=0d88af9f9f907d4b126e6fc1a6a1d73aCAS |

Li WD, Zhang CR, Burt JE, Zhu AX (2005) A Markov chain-based probability vector approach for modeling spatial uncertainties of soil classes. Soil Science Society of America Journal 69, 1931–1942.

Liu YP, Tong J, Li XN (2005) Analysing the silt particles with the Malvern Mastersizer 2000. Water Conservancy Science and Technology and Economy 11, 329–331. [in Chinese with English abstract]

Liu B, Zhao WZ, Chang XX, Li SB, Zhang ZH, Du MW (2010) Water requirements and stability of oasis ecosystem in arid region, China. Environmental Earth Science 59, 1235–1244.
Water requirements and stability of oasis ecosystem in arid region, China.Crossref | GoogleScholarGoogle Scholar |

Lobe I, Amelung W, Du Preez CC (2001) Losses of carbon and nitrogen with prolonged arable cropping from sandy soils of the South Africa Highveld. European Journal of Soil Science 52, 93–101.
Losses of carbon and nitrogen with prolonged arable cropping from sandy soils of the South Africa Highveld.Crossref | GoogleScholarGoogle Scholar |

Logofet DO, Lesnaya EV (2000) The mathematics of Markov models: what Markov chains can really predict in forest successions. Ecological Modelling 126, 285–298.
The mathematics of Markov models: what Markov chains can really predict in forest successions.Crossref | GoogleScholarGoogle Scholar |

Miller J, Franklin J (2002) Modeling the distribution of four vegetation alliances using generalized linear models and classification trees with spatial dependence. Ecological Modelling 157, 227–247.
Modeling the distribution of four vegetation alliances using generalized linear models and classification trees with spatial dependence.Crossref | GoogleScholarGoogle Scholar |

Nielsen DR, Bouma J (1985) Soil spatial variability. In ‘Proceedings of Workshop of the ISSS and the SSSA’. 30 Nov.–1 Dec. 1984, Las Vegas, NV. (Pudoc: Wageningen, the Netherlands)

Pan XL (2001) A preliminary study on the stability of oasis ecosystem in arid area. Quaternary Sciences 21, 345–351. [in Chinese with English abstract]

Shi RJ (1992) ‘The basic theory of Markov chain and its application.’ (Xi’an Electronic and Technical University Publishing: Xi’an, China) [in Chinese].

Shirazi MA, Boersma L (1984) A unifying quantitative analysis of soil texture. Soil Science Society of America Journal 48, 142–147.
A unifying quantitative analysis of soil texture.Crossref | GoogleScholarGoogle Scholar |

Sivakumar B, Harter T, Zhang H (2005) A fractal investigation of solute travel time in a heterogeneous aquifer: transition probability/Markov chain representation. Ecological Modelling 182, 355–370.
A fractal investigation of solute travel time in a heterogeneous aquifer: transition probability/Markov chain representation.Crossref | GoogleScholarGoogle Scholar |

Su YZ, Zhao HL, Zhao WZ, Zhang TH (2004) Fractal features of soil particle size distribution and the implication for indicating desertification. Geoderma 122, 43–49.
Fractal features of soil particle size distribution and the implication for indicating desertification.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD2cXlvVWmtr4%3D&md5=765f5a32befcf2df346c699d1959ed57CAS |

Su YZ, Zhao WZ, Su PX, Zhang ZH, Wang T, Ram R (2007) Ecological effects of desertification control and desertified land reclamation in an oasis-desert ecotone in an arid region: A case study in Hexi Corridor, northwest China. Ecological Engineering 29, 117–124.
Ecological effects of desertification control and desertified land reclamation in an oasis-desert ecotone in an arid region: A case study in Hexi Corridor, northwest China.Crossref | GoogleScholarGoogle Scholar |

Sun DZ, Li XQ (2010) Application of Markov Chain model on environmental fate of phenanthrene in soil and groundwater. Procedia Environmental Science 2, 814–823.
Application of Markov Chain model on environmental fate of phenanthrene in soil and groundwater.Crossref | GoogleScholarGoogle Scholar |

Usher MB (1981) Modelling ecological succession, with particular reference to Markovian models. Vegetation 46–47, 11–18.
Modelling ecological succession, with particular reference to Markovian models.Crossref | GoogleScholarGoogle Scholar |

Weigand H, Totsche KU, Huwe B, Kögel-Knabner I (2001) PAH mobility in contaminated industrial soils: a Markov chain approach to the spatial variability of soil properties and PAH levels. Geoderma 102, 371–389.
PAH mobility in contaminated industrial soils: a Markov chain approach to the spatial variability of soil properties and PAH levels.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD3MXkvF2hsr4%3D&md5=12cb6e3bc8463d39de7d7b75947a6fb0CAS |

Wu KJ, Nunan N, Crawford JW, Young IM, Ritz K (2004) An efficient Markov chain model for the simulation of heterogeneous soil structure. Soil Science Society of America Journal 68, 346–351.
An efficient Markov chain model for the simulation of heterogeneous soil structure.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD2cXitV2ntro%3D&md5=6ec4a54f9087e8e6bc5c0252f0abc5f7CAS |

Ye M, Khaleel R (2008) A Markov chain model for characterizing medium heterogeneity and sediment layering structure. Water Resources Research 44, W09427
A Markov chain model for characterizing medium heterogeneity and sediment layering structure.Crossref | GoogleScholarGoogle Scholar |

Zhang CR, Li WD (2005) Markov chain modeling of multinomial land-cover classes. GIScience & Remote Sensing 42, 1–18.
Markov chain modeling of multinomial land-cover classes.Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD2MXntVCjtLc%3D&md5=6673acec351b0c9a81be427af351fd93CAS |

Zhang H, Wu JW, Zheng QH, Yu YJ (2003) A preliminary study of oasis evolution in the Tarim Basin, Xinjiang, China. Journal of Arid Environments 55, 545–553.
A preliminary study of oasis evolution in the Tarim Basin, Xinjiang, China.Crossref | GoogleScholarGoogle Scholar |

Zhang RQ, Tang CJ, Ma SH, Yuan H, Gao LL, Fan WY (2011) Using Markov chains to analyze changes in wetland trends in arid Yinchuan Plain, China. Mathematical and Computer Modelling 54, 924–930.
Using Markov chains to analyze changes in wetland trends in arid Yinchuan Plain, China.Crossref | GoogleScholarGoogle Scholar |

Zhu AX, Mackay DS (2001) Effects of spatial detail of soil information on watershed modeling. Journal of Hydrology 248, 54–77.
Effects of spatial detail of soil information on watershed modeling.Crossref | GoogleScholarGoogle Scholar |