The purpose of this study was to evaluate the predictive accuracy of habitat suitability models, identifying the potential distribution range of Dorema ammoniacum, and its habitat requirements in the rangelands of Yazd province, central Iran. Bafgh, Mehriz and Nadoushan, were three habitats that were identified, and sampling was conducted in each habitat using a random-systematic method. A set of 10 plots were established (at equal distances) along 350 m long 18 transects. Soil samples (two depths: 0–30 and 30–60 cm from 36 profiles) were collected and measured in the laboratory. Elevation, slope, and aspect maps were derived, and climate information was collected from nearby meteorological stations. The habitat prediction of the species was modeled using Logistic Regression (LR), Maximum Entropy (MaxEnt), and Artificial Neural Network (ANN). The Kappa coefficient and the area under the curve (AUC) were calculated to assess the accuracy of the forecasted maps. The LR model for habitat prediction of the studied species in Mehriz (K = 0.67) and Nadoushan (K = 0.56) habitats were identified as good. The MaxEnt model predicted the habitat distribution for the selected species in Bafgh and Mehriz habitats as excellent (K = 0.89, AUC = 0.76, K = 0.89, AUC = 0.98), and in the Nadoushan habitat as very good (K = 0.78, AUC = 0.85). However, the ANN model predicted Bafgh and Nadoushan habitats as excellent and Mehriz habitat as very good (K = 0.87, K = 0.90, and K = 0.63, respectively). In general, in order to protect species D. ammoniacum, the development of its habitats in other areas of Yazd province and the habitats under study in conservation programs should be given priority.