For sustainable land use planning, evaluating land characteristics and making suitable land use decisions is a priority and critical step. In order to make these evaluations safely, spatial analyzes of many criteria should be made. In this study, the suitability of the land for wheat production was evaluated by Geographical Information Systems (GIS) based Multiple Criteria Decision Analysis (MCDA) in semi-arid conditions. In obtaining the land suitability map; fuzzy set model, Analytical Hierarchy Process (AHP) and GIS are integrated. Ecological criteria weights for agricultural land suitability were determined by AHP. In the suitability analysis, a total of criteria including soil and topographic features were evaluated. Geostatistical analysis approach was applied to determine the spatial variability of soil properties (sand, clay, silt, pH, OM, CEC, ESP, CaCO3, EC). The lowest variation among soil properties was observed in pH (3.8%), while the largest variation was observed in ESP content (107.5%). The nugget/sill ratio is poor for EC and pH, while other soil properties are moderately spatially dependent. According to the results of the analysis, 25.7% (3.226 km2) of the area is highly suitable, while 27.6% (3.457 km2) is moderately suitable and 19.5% (2.440 km2) is marginally suitable for wheat cultivation. In addition, 27.2% (3.415 km2) of the area is not suitable for agricultural production. The use of geostatistical modeling, MCDA and GIS together is very beneficial in making agricultural land management decisions.