Digital soil mapping techniques which incorporate the digital auxiliary environmental data to field observation data using software are more reliable and efficient compared to conventional surveys. Therefore, this study has been conducted to use k- Nearest Neighbors (k-NN) and artificial neural network (ANN) to predict spatial variability of soil salinity in Ardekan district in an area of 700 km², in Yazd province. In this study, 180 soil samples were collected in a grid sampling manner and then soil chemical and physical properties were measured in laboratory. Environmental auxiliary variables were included topographic attributes, remote sensing data (ETM+) and apparent electrical conductivity (ECa). The result of the study showed that the K-mean nearest neighborhood had higher accuracy than ANN models for predicting soil electrical conductivity (ECe). Overall, k-NN models could provide significant relationships between soil salinity data and environmental auxiliary variables. The k-NN model had the root mean square and coefficient of determination of 12.10 and 0.92, respectively, between predicted and observed ECe data. Also, apparent EC, and remotely sensed indices and wetness index were identified as the most important factors for predicating the soil salinity in the studied area.