zohreh mosleh; Mohammad hasan Salehi; azam jafari; Abdolmohammad Mehnatkesh; Isa Esfandiarpoor Borujeni
Abstract
Introduction: There is a concern with assessment of land performance when used for specific purposes. Land evaluation analysis is considered as an interface between land resources and land use planning and management. However, the conventional soil surveys are usually not useful for providing quantitative ...
Read More
Introduction: There is a concern with assessment of land performance when used for specific purposes. Land evaluation analysis is considered as an interface between land resources and land use planning and management. However, the conventional soil surveys are usually not useful for providing quantitative information about the spatial distribution of soil properties that are used in many environmental studies. Development of the computers and technology lead to digital and quantitative approaches have been developed. These new techniques rely on finding the relationships between soil and the auxiliary information that explain the soil forming factors or processes and finally predict soil patterns on the landscape. Different types of the machine learning approaches have been applied for digital soil mapping of soil classes, such as the logistic and multinomial logistic regressions, neural networks and classification trees. To our knowledge, most of the previous studiesapplied land suitability evaluation based on the conventional approach. Therefore, the main objective of this study was to assess the performance of digital mapping approaches for the qualitative land suitability evaluation in the Shahrekord plain of Chaharmahal-Va- Bakhtiari province.
Materials and Methods: An area in the Shahrekord plain of Chaharmahal-Va-Bakhtiari Province, Iran, across 32º13′ and 32º 23′N, and 50º 47′ and 51º 00′E was chosen. The soils in the study area have been formed on Quaternary shale and foliated clayey limestone deposits. Irrigated crops such as wheat, potato, maize and alfalfa are the main land uses in the area. According to the semi-detailed soil survey, 120 pedons with approximate distance of 750 m were excavated and soil samples were taken from different soil horizons. Soil physicochemical properties were determined. The average of soil properties was determined by considering the depth weighted coefficient up to 100 and 150 centimeters for annual and perennial crops, respectively. Qualitative land suitability evaluation for main crops of the area including wheat, maize, alfalfa and potato was determined by matching the site conditions (climatic, hydrology, vegetation and soil properties) with studied crop requirement tables presented by Givi (5). Land suitability classes were determined using parametric method. Land suitability classes reflect degree of suitability as S1 (suitable), S2 (moderately suitable), S3 (marginally suitable) and N (unsuitable). Different machine learning techniques, namely artificial neural networks (ANNs), boosted regression tree (BRT), random forest (RF) and multinomial logistic regression (MLR) were used to test the predictive power for mapping the land suitability evaluation. Terrain attributes, normalized difference vegetation index (NDVI), clay index, carbonate index, perpendicular vegetation index (PVI), geology map, existing soil map (1:50000 scale) and geomorphology map were used as auxiliary information. Finally, all of the environmental covariates were projected onto the same reference system (WGS 84 UTM 39 N) and resampled to 50 * 50 m since the soil samples were collected with approximate distance of 750 m (1:50,000 scale). According to the suggested resolutions for digital soil maps, the pixel size 50 *50 m fits to a 1:50,000 cartographic scale. Training the models was done with 80% of the data (i.e., 96 pedons) and their validation was tested by the remaining 20% of the dataset (i.e., 24 pedons) that were split randomly. The accuracy of the predicted soil classes was determined using error matrices and overall accuracy.
Results and Discussion: The results showed that climatic conditions are suitable (S1) for wheat and potato whereas the most important limiting factors for maize and alfalfa were the average of minimum temperature and average temperature, respectively. Results demonstratedthat among the studied models, random forest showed the highest performance to predict the land suitability classes and subclasses. However, different models had the same ability for prediction. In addition, the overall accuracy decreased from class to subclass for all of the crops. The terrain attributes and remote sensing indices (normalized difference vegetation index and perpendicular vegetation index) were the most important auxiliary information to predict the land suitability classes and subclasses.
Conclusion: Results suggest that the DSM approaches have enough accuracy for prediction of the land suitability classes that affecting land use management. Although digital mapping approaches increase our knowledgeabout the variation of soil properties, integrating the management of the sparse lands with different owners should be considered as the first step for optimum soil and land use management.
A. Mehnatkesh; S. Ayoubi; A. Jalalian
Abstract
Introduction: Soil depth is defined as the depth from the surface to more-or-less consolidated material and can be considered as the most crucial soil indicator, affecting desertification and degradation in disturbed ecosystems. Soil depth varies as a function of many different factors, including slope, ...
Read More
Introduction: Soil depth is defined as the depth from the surface to more-or-less consolidated material and can be considered as the most crucial soil indicator, affecting desertification and degradation in disturbed ecosystems. Soil depth varies as a function of many different factors, including slope, land use, curvature, parent material, weathering rate, climate, vegetation cover, upslope contributing area, and lithology. Topography, one of the major soil forming factors, controls various soil properties. Thus, quantitative information on the topographic attributes has been applied in the form of digital terrain models (DTMs). The prediction of soil depth by topographic attributes depends mainly on: i) the spatial scale of topographic variation in the area, ii) the nature of the processes that are responsible for spatial variation in soil depth, and iii) the degree to which terrain-soil relationships have been disturbed by human activities. This study was conducted to explore the relationships of soil depth with topographic attributes in a hilly region of western Iran.
Materials and Methods: The study area is located at Koohrang district between 32°20′ to 32°30′ N latitudes and 50°14′ to 50°24′ E longitudes, in Charmahal and Bakhtiari province, western Iran. The field sites with an area of 30,000 ha are located on the hillslopes at about 20% transversal slope. The soils at the site are classified as Typic Calcixerepts, Typic Xerorthents and Calcic Haploxerepts for the representative excavated profiles in summit, shoulder and backslope, respectively. The soils located at footslope and toeslope were classified as Chromic Calcixererts. Measurements were made in twenty representative hillslopes of the studied area. At the selected site, one hundred points were selected using randomly stratified methodology, considering all geomorphic surfaces including summit, shoulder, backslope, footslope and toeslope during sampling. Overall, 100 profiles were dug and described; and the solum thickness was measured for each profile. DEM data were created by using a 1:2,5000 topographic map. Topographical indices were generated from the DEM using TAS software. Terrain attributes in two categories, primary and secondary (compound) attributes; primary attributes are included elevation, slope, aspect, catchment area, dispersal area, plan curvature, profile curvature, tangential curvature, shaded relief. Secondary or compound attributes such as soil water content or the potential for sheet erosion, stream power index, wetness index, and sediment transport index. Correlation coefficients to define relationships between soil depth and terrain attributes, and analysis of variance by Duncan test were done using the SPSS software. The statistical software SPSS was used for developing multiple linear regression models. Terrain attributes were selected as the independent variables and soil depth was employed as dependent variable in the model. Thirty sampling sites were used to validate the developed soil-landscape model. In testing soil-landscape model, we calculated two indices from the observed and predicted values included mean error (ME) and root mean square error (RMSE).
Results and Discussion: The soil depth in the studied profiles varied from 30 cm to 150 cm with an average of 108.6 cm. Relatively high variability (CV = 76%) was obtained for soil depth in the study area. The linear correlation analysis of the 12 topographic attributes and one soil property (soil depth), showed that there was a significant correlation among 36 of the 77 attribute pairs. Soil depth showed high positive significant correlations with catchment area, plan curvature, and wetness index, and showed high negative correlation with sediment transport index, sediment power index and slope. Low positive significant correlations of soil depth were identified with tangential curvature, and profile curvature. Moreover, soil depth was negatively correlated with elevation. The rest of the topographic attributes including aspect, shaded relief, and dispersal area were not significantly correlated with soil depth. Many of these relationships are similar to those found in other landscapes. The results of analysis of variance showed that there are significant differences for soil depth among the selected slope positions in the studied area. The highest values of soil depth were observed in the downslope positions including footslope and toeslope. The lowest soil depth was observed in shoulder position with the highest rate of soil erosion.
Conclusions: It seems that the high variability for soil depth depends on topography of the field, and the landscape position, causing differential accumulation of water at different positions on the landscape; and moreover the soil erosion and deposition processes, resulting in high variability in the soil depth. We found relatively high correlation coefficients of soil depth with two groups of topographic attributes (erosional processes and water accumulation). Empirical model (MLR) using selected terrain attributes explains 76% of the variation of soil depth in the studied area. The terrain attributes that best predicted soil depth variability in the selected site were mainly the attributes that had significant relationships with soil depth. The dominant attributes in the MLR model included slope, wetness index, catchment area and sediment transport index.
H. Afshar; M.H. Salehi; J. Mohammadi; A. Mehnatkesh
Abstract
Abstract
The quality of soil maps depends upon their ability to show the soils variability. Thus, the accuracy of the maps used for crop recommendations is due to the accuracy of soil maps. This study was performed to investigate the amount of soil properties and crop yield spatial variability in S ...
Read More
Abstract
The quality of soil maps depends upon their ability to show the soils variability. Thus, the accuracy of the maps used for crop recommendations is due to the accuracy of soil maps. This study was performed to investigate the amount of soil properties and crop yield spatial variability in S 2 and S3 units of a semi-detailed quantitative suitability map (1:50000 scale) for irrigated wheat in Shahr-e-Kian area, Chaharmahal-Va-Bakhtiari province. Eighty soil samples were collected in each land unit at 0-30 cm depth using multi-scale sampling method to determine available P, K, total N, %O.M., %CaCO3 equivalent, soil texture and particle size distribution, EC and pH. A 0.5×0.5 m plot of wheat was harvested at each of 160 sites previously sampled to determine crop biomass, 1000 seeds weight and harvest index. The highest CV was related to available potassium (47.43 for S2 and 46.46 for S3 units, respectively) and the lowest one was related to pH (1.07 for S2 and 0.925 for S3 units, respectively). Variography showed a good spatial structure for all variables in both land units. Ranges for variograms were from 17.75 for N to 61.06 m for EC in S2 unit and from 17.47 for P to 62.93 m for 1000 seeds weight in S3 unit. Kriging maps showed high spatial variability of soil properties as well as biomass, wheat yield and harvest index within two land units. This indicates that suitability maps have not enough credibility for precision agriculture. Using information of all pedons as well as representative pedons in land units and combining the information of suitability maps with geostatistical data can be a choice way to improve the accuracy and quality of land suitability maps.
Keywords: Kriging, Precision agriculture, Soil properties, Spatial variability, Suitability map, Wheat yield