This study investigates the ability of different digital soil mapping (DSM) approaches to predict some of physical and chemical topsoil properties in the Shahrekord plain of Chaharmahal-Va-Bakhtiari province, Iran. According to a semi-detailed soil survey, 120 soil samples were collected from 0 to 30 cm depth with approximate distance of 750 m. Particle size distribution, coarse fragments (CFs), electrical conductivity (EC), pH, organic carbon (OC), and calcium carbonate equivalent (CCE) were determined. Four machine learning techniques, namely, artificial neural networks (ANNs), boosted regression tree (BRT), generalized linear model (GLM), and multiple linear regression (MLR), were used to identify the relationship between soil properties and auxiliary information (terrain attributes, remote sensing indices, geology map, existing soil map, and geomorphology map). Root-mean-square error (RMSE) and mean error (ME) were considered to determine the performance of the models. Among the studied models, GLM showed the highest performance to predict pH, EC, clay, silt, sand, and CCE, whereas the best model is not necessarily able to make accurate estimation. According to RMSE%, DSM has a good efficiency to predict soil properties with low and moderate variabilities. Terrain attributes were the main predictors among different studied auxiliary information. The accuracy of the estimations with more observations is recommended to give a better understanding about the performance of DSM approach over low-relief areas.