To predict the soil texture fractions, 115 profiles were identified based on the Latin hypercube sampling technique, the horizons were sampled, and the clay, sand, and silt contents (in percentages) of soil samples were measured. Then equal-area quadratic spline depth functions were used to derive clay, sand, and silt contents at five standard soil depths (0–5, 5–15, 15–30, 30–60, and 60–100 cm). Auxiliary variables used in this study include the terrain attributes (derived from a digital elevation model), Landsat 8 image data (acquired in 2015), geomorphological map, and spectrometric data (laboratory data). Artificial neural network (ANN), regression tree (RT), and neuro-fuzzy (ANFIS) models were used to make a correlation between soil data (clay, sand, and silt) and auxiliary variables. The results of this study showed that the ANFIS model was more accurate in the prediction of the three parameters of clay, silt, and sand than ANN and RT. Moreover, the ability of ANFIS model to estimate the soil texture fractions in the surface layers was higher than the lower layers. The mean coefficient of determination (R²) values calculated by 10-fold cross validation suggested the higher prediction performance in the upper depth intervals and higher prediction error in the lower depth intervals (e.g., R² = 0.91, concordance correlation coefficient (CCC) = 0.90, RMSE = 4.00 g kg⁻¹ for sand of 0–5 cm depth, and R² = 0.68, CCC = 0.60, RMSE = 8.03 g kg⁻¹ for 60–100 cm depth). The results also showed that the most important auxiliary variables are spectrometric data, multi-resolution, valley-bottom flatness index and wetness index. Overall, it is recommended to use ANFIS models for the digital mapping of soil texture fractions in other arid regions of Iran.