Document Type : Research Article

Authors

University of Kurdistan

Abstract

Introduction: Soil organic matter (SOM) is an important soil quality factor that affects physical, chemical and biological properties of soil. Accurate estimation of SOM variability provides critical information especially in precision agriculture. Geostatistics and geographic information system (GIS) are powerful tools for characterizing and mapping the spatial distribution and variability of soil properties. Kriging is a basic geostatistical technique that provides the best linear unbiased estimation for a spatially dependent variable. This method will produce satisfying results if enough sample points are available. Unfortunately, laboratory measurements of the SOM are costly and time-consuming. Artificial neural network-kriging (ANNK) is another geostatistical method that extends kriging of a primary variable to the readily available auxiliary variables based on their relationship with the primary variable. This relationship is captured using an artificial neural network (ANN) model. The residuals of the model were then interpolated using kriging, and added to the prediction obtained from the ANN model. Terrain attributes, derived from digital elevation models (DEMs), are useful for estimating SOM at landscape scale. Topographic indicators including slope, aspect, elevation, and topographic wetness index may be the dominant factors affecting SOM variability in an area with same parent material and climate. Hence, these factors can be used as auxiliary variables for estimating spatial variability of SOM using ANNK. The objective of this study was to estimate SOM spatial variability using ANNK and topographic indices and assess its status in hilly areas of Ghorveh in Kurdistan province (Iran).
Materials and Methods: A total of 150 soil samples from a depth of 0-15 cm were systematically collected in a grid spaced 2 Km × 2 Km. The SOM content of soil samples was measured in the laboratory. Topographic indicators including slope, aspect, elevation, and topographic wetness index were derived from the DEM. ANN was used to predict SOM variability based on topographic index combinations. The feed-forward network consisted of an input layer, one hidden layer with sigmoid neurons, and an output layer with linear neurons. The network was trained with Levenberg-Marquardt backpropagation algorithm. According to the Kolomogrov’s theorem, the number of nodes in the hidden layer was 2n+1, in which n is the number of input neurons. The optimal subset of topographic index combinations correlating best with the SOM was selected as the best ANN model. This model was used to generate an initial SOM surface. The residuals of ANN model were interpolated using ordinary kriging (OK) and combined with the initial SOM surface to produce the final ANNK SOM surface. The SOM status map was derived from overlaying of soil texture and SOM maps in four different levels (very low, low, medium and high).
Results and Discussion: The results of ANN suggested that elevation was the most important variable determining the distribution of SOM across the landscape. Further, aspect was the other variable which had a significant influence on SOM distribution. The selected two inputs ANN model (elevation and aspect) can explain about 33% of total variance of SOM. The cross-validation results indicated that the OK and ANNK techniques can explain about 37 and 89% of total variance of SOM, respectively. The ANNK technique performed better than the OK and ANN techniques since it was able to capture most of the small variations of SOM. The resulting SOM status map indicated a low and very low SOM content in relation with soil texture in most regions surveyed (79%). Low SOM level can be attributed to the erosive processes under Mediterranean climate on hills coupled with intensive and/or inappropriate agricultural practices. Based on the results of this study, proper agronomical and environmental planning such as soil conservation strategy is highly required in this area to restore and increase the SOM content in agricultural soils, combat soil erosion and maintain soil ecological functions and productivity. The SOM replenishment can be achieved in the degraded areas (i.e., low SOM content) by adopting conservative practices such as conservation tillage or no-tillage (e.g., direct seeding), improving land use rotations with forage crops, returning crop residues to soil, growing green manure crops, and supplying the soil with proper exogenous organic matter (compost, manure, sewage sludge, etc.). Furthermore, the results highlighted the potential of ANNK in combination with GIS to provide improved distribution patterns of SOM.

Keywords

1- Beven K.J., and Kirkby M.J. 1979. A physically based, variable contributing area model of basin hydrology. Hydrological Sciences Bulletin, 24:43–69.
2- Bishap C.M. 1995. Neural networks for pattern recogniation. Oxford University Press, Oxford, 482.
3- Burgess T.M., and Webster, R. 1980. Optimal interpolation and isarithmic mapping of soil properties: the semivariogram and punctual kriging. Soil Science, 31:315–331.
4- Burrough P.A. 1986. Principles of geographical information systems for land resources assessment. Oxford university press, New York.
5- Chivenge P.P., Murwira H.K., Giller K.E., Mapfumo P., and Six J. 2007. Long-term impact of reduced tillage and residue management on soil carbon stabilization: implications for conservation agriculture on contrasting soils. Soil and Tillage Research, 94:328–337.
6- Eldeiry A., and Garcia L.A. 2009. Comparison of Regression kriging and cokriging techniques to estimate soil salinity using Landsat images. Hydrology Days, 27:38.
7- Guo P.T., Wu W., and Sheng Q.K. 2013. Prediction of soil organic matter using artificial neural network and topographic indicators in hilly areas. Nutrient Cycling in Agroecosystems, 95:333-344.
8- Ingleby H.R., and Crowe T.G. 2001. Neural network models for predicting organic matter content in Saskatchewan soils. Canadian Biosystems Engineering, 43:71-75.
9- Istok J.D., Smyth J.D., and Flint A.L. 1993. Multivariate geostatistical analysis of groundwater contaminant: a case history. Groundwater, 31:63–74.
10- Li Z.Y. 1998. Supervised classification of multi-spectral remote sensing image using B-P neural network. Journal of Infrared and Millimeter Waves, 17:153-156.
11- Liao K., Xu S., Wu J., and Zhu Q. 2013. Spatial estimation of surface soil texture using remote sensing data. Soil Science and Plant Nutrition, 59(4):488-500.
12- Mahmoodi S., and Hakimian M. 1998. Fundamentals of soil science. Tehran university press, Tehran. (In Persian)
13- Marchetti A., Piccini C., Francaviglia R., and Mabit L. 2012. Spatial distribution of soil organic matter using geostatistics: A key indicator to assess soil degradation status in central Italy. Pedosphere, 22(2):230–242.
14- McBratney A.B., and Webster R. 1986. Choosing functions for semivariograms of soil properties and fitting them to sampling estimates. Journal of Soil Science, 37:617–639.
15- McBratney A.B., Santos M.L.M., and Minasny B. 2003. On digital soil mapping. Geoderma, 117: 3-52.
16- Mirzaie M. 2015. Prediction of soil organic matter based on soil characteristics, topography and remote sensing data using artificial neural networks. M.Sc. thesis, University of Kurdistan, Sanandaj. (In Persian with English abstract)
17- Moore A.D., McLaughlin R.A., Mitasova H., and Line D.E. 2007. Calibrating WEPP model parameters for erosion prediction on construction sites. Transactions of the ASABE, 50(2):507-516.
18- Mulla D.J., and McBratney A.B. 2000. Soil Spatial Variability. p. 343–373. In A.W. Warrick (ed.) Soil Physics Companion. CRC Press, NewYork.
19- Nash J.E., and Sutcliffe J.V. 1970. River flow forecasting through conceptual models: Part I. A discussion of principles. Journal of Hydrology,10(3):282-290.
20- Odeh I.O.A., McBratney A.B., and Chittleborough D.J. 1995. Further results on prediction of soil properties from terrain attributes: Heterotopic cokriging and regression-kriging. Geoderma, 67(3-4):215–226.
21- Quinton J.N. 1997. Reducing predictive uncertainty in model simulations: A comparison of two methods using the European Soil Erosion Model (EUROSEM). Catena, 30(2): 101-117.
22- Shouse P.J., Gerik T.J., Russell W.B., and Cassel D.K. 1990. Spatial distribution of soil particle size and aggregate stability index in a clay soil. Soil Science, 149:351–360.
23- Societ`a Italiana dei Laboratori Pubblici di Agrochimica (SILPA). 1999. From soil analysis to the fertilization advice. ASSAM, Agenzia Servizi SettorenAgroalimentare delle Marche, Regione Marche, Jesi, Italy. (In Italian)
24- Triantafilis J., Odeh I.O.A., and McBratney A.B. 2001. Five geostatistical models to predict soil salinity from electromagnetic induction data across irrigated cotton. Soil Science Society of America Journal, 65:869–878.
25- Vauclin M., Vieira S.R., Vachaud G., and Nielsen D.R. 1983. The use of cokriging with limited field observations. Soil Science Society of America Journal, 47:175–184.
26- Walkley A., and Black I.A. 1934. An examination of the Degtjareffmethod for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Science, 37:29–38.
27- Wu C., Wu J., Luo Y., Zhang L., and DeGloria S.D. 2009. Spatial prediction of soil organic matter content using cokriging with remotely sensed data. Soil Science Society of America Journal, 73:1202–1208.
28- Yadav V., and Malanson G. 2007. Progress in soil organic matter research: litter decomposition, modeling, monitoring and sequestration. Progress in Physical Geography, 31:131–154.
29- Yates S.R., and Warrick A.W. 1987. Estimating soil water content using cokriging. Soil Science Society of America Journal, 51:23–30.
30- Zhao Z., Yang Q., Benoy G., Chow T.L., Xing Z., Rees H.W., and Meng F.R. 2010. Using artificial neural network models to produce soil organic carbon content distribution maps across landscapes. Canadian Journal of Soil Science, 90:75-87.
CAPTCHA Image