S. Nazari; M. Rostaminia; shamsollah Ayoubi; A. Rahmani; S.R. Mousavi
Abstract
Abstract Background and objectives: High-accuracy of soil maps is a powerful tool for achieving land sustainability in agricultural and natural resources. The present study was conducted in Vargar lands of Abdanan city related to Ilam province for digital mapping of soil classes at two taxonomic level ...
Read More
Abstract Background and objectives: High-accuracy of soil maps is a powerful tool for achieving land sustainability in agricultural and natural resources. The present study was conducted in Vargar lands of Abdanan city related to Ilam province for digital mapping of soil classes at two taxonomic level from subgroup up to family by random forest (RF) and fuzzy logic models. Materials and methods: Study area with 1027 hectare have 628.6 mm and 22.6 C° mean annual precipitation and temperature respectively. Three major physiographic units included Hilland, Piedmont plain and Alluvial plain were observed. Soil moisture and temperature regimes are ustic and hyperthermic calculated based on Newhall model in JNSM 6.1 version software. A total of 44 soil profile observation with random sampling pattern was determined based on standardized soil surveys then digging, description and after sampling from all genetic horizons then soil samples were transferred to laboratory. Finally, all of soil profiles were classified based on soil taxonomy system (2014) up to family level. Geomorphometric covariates as a representative of soil forming factors were prepared from digital elevation model (ALOS PALSAR Satellite,2011) with 12.5 m resolution in SAGA GIS 7.4 version software. Three feature selection approaches included Boruta, Variance inflation factors (VIF) and Mean decrease accuracy (MDA) with two Random forest (RF) and Fuzzy logic data mining algorithms were applied for relating soil-landscape relationship by using “randomforest”, “caret” packages in R 3.5.1 and SoLIM solution version 2015 software. Sample based project used for predicting soil classes in Fuzzy logic modeling process. In totally observation profile split into two data set included 80 percent (n=36) for calibrating and 20 percent for validating (n=8) based on bootstraps sampling algorithm random forest. Internal validation of random forest algorithm was done based on out of bag error percentage (OOB%). The best model performance was determined based on overall accuracy (OA) and kappa index, also for each individual class user accuracy (UA) and producer accuracy (PA) were applied. Results: The results shown that from number of 40 geomorphometrics covariates, six covariates included Terrain classification index for lowlands, Annual insolation, Topographic position Index, Upslope curvature, Real surface area and Terrain surface convexity were selected by MDA as the best environmental covariates. Also, RF-MDA method with overall accuracy 84% and Kappa index 0.56 had the best performance compared to other methods (RF_VIF, RF-BO, Fuzzy-MDA) in subgroup level with 58, 55, 50 and 0.3, 0.67 and 0.18 respectively. Out of bag error results (%OOB) for RF-MDA, RF-VIF and RF-Boruta were obtained that 72.42%, 67.86% and 82.76% for subgroup level and 93.10%, 93.10% and 86.21% for family level respectively. while there was little difference between the accuracy of the method at the family taxonomic level and performed similar results in modeling of soil classes process. The results of the fuzzy approach showed that the kappa index values and overall accuracy of this method were similar to the other three scenarios and there was a slight difference between the accuracy of the results at the soil family level. In the fuzzy method, it was observed that the kappa and overall accuracy values at the subgroup level were lower than the other scenarios. Fuzzy approaches in contrasted to RF modeling prevented continues spatial variability by generating of fuzzy maps for each of soil class in the landscape. These results indicate that the random forest method is superior to the fuzzy method in family class mapping and soil subgroups. Based on MDA sensitivity analysis index, similarly, three geomorphometrics covariate included Terrain surface convexity (convexity), Terrain classification index for lowlands (TCI_Low) and Real surface area (Surface_Ar) had highest importance for predicting soil classes at two taxonomic level. With regarded to final soil predicted maps area, two classes (Fine-silty, carbonatic, hyperthermic Typic Haplustepts) and Typic Calciustolls with 32.70% and 48.90% and (Fine-silty, carbonatic, hyperthermic Typic Calciustolls) and Typic Haplustepts with 0.18% and 1.85% had the highest and lowest content at family and subgroup maps respectively. Conclusion: In general, using different variable selection approaches in situations where soil classes have a relatively imbalanced abundance can increase the accuracy of digital mapping in soil studies. Increasing the number of field observations and the use of other environmental variables affecting soil formation can also be used for gradating in prediction low-accuracy soil classes.