Land suitability classification is a useful management practice to ensure planned and sustainable use of agricultural lands according to their potentials. The main purposes of this study were to analyze land suitability for bread wheat (Triticum aestivum) cultivation and generate a land suitability map for wheat by integrating the analytical hierarchy (AHP)-fuzzy algorithm with the Geographical Information System (GIS) in the Tozanlı sub-basin located in the upper part of Yeşilırmak Basin, Turkey. Topographic (elevation, slope, aspect) characteristics of the basin and some of physical and chemical properties of soils (texture, pH, electrical conductivity, lime, organic matter, and soil depth) were used as criteria in determining the suitability classes. Ninety-two disturbed soil samples were collected from 0 to 20 cm depth in October 2017 using random sampling method. Weighted overlay spatial analysis in GIS was used to combine different thematic layers to identify areas suitable for wheat production. The fuzzy-AHP suitability assessment model was adapted to determine the weights for topographic and soil properties. The highest specific weights were obtained for soil depth (0.232) and elevation (0.218), while the lowest weight was calculated for aspect (0.042). Highly, moderately, and marginally suitable lands for wheat cultivation cover 2.63, 9.85 and 32.59% of the study area, respectively. In addition, the results indicated that 54.92% of the total area is permanently unsuitable for wheat cultivation. The results revealed that integration of AHP-fuzzy algorithm and GIS techniques is a useful method for accurate evaluation of land suitability in planning for specific crop production and decreasing the negative environmental impacts of agricultural practices.