Land suitability assessment for wheat-barley cultivation in a semi-arid region of Eastern Anatolia in Turkey

The efficient use and sustainability of agricultural lands depend heavily on the characteristics of soil resources in a given area, as different soil properties can significantly impact crop growth and yield. Therefore, land suitability studies play a crucial role in determining the appropriate crops for a given area and ensuring sustainable agricultural practices. This study, conducted in Tusba District-Van, Turkey, represents a significant advancement in land suitability studies for wheat-barley cultivation. Using geographic information systems, the analytical hierarchical process method, and the standard scoring function, lands were determined based on the examined criteria for the suitability of wheat-barley cultivation. One of this study’s main findings is identifying critical factors that influence the suitability of land for wheat-barley cultivation. These factors include slope, organic matter content, available water capacity, soil depth, cation exchange capacity, pH level, and clay content. It is important to note that slope is the most influential factor, followed by organic matter content and available water capacity. A Soil Quality Index map was produced, and the suitability of wheat-barley production in the studied area was demonstrated. More than 28% of the study area was very suitable for wheat-barley production (S2), and more than was 39% moderately suitable (S3). A positive regression (R2 = 0.67) was found between soil quality index values and crop yield. The relationship between soil quality index values and crop yield is above acceptable limits. Land suitability assessment can minimize labor and cost losses in the planning and implementation of sustainable ecological and economic agriculture. Furthermore, land suitability classes play an active role in the selection of the product pattern of the area by presenting a spatial decision support system.


INTRODUCTION
In order to meet the increasing food demand with the rising world population, agricultural production should be brought to a level that meets these needs.Effective and appropriate use of agricultural lands is essential in protecting soil and water resources.Excessive use of soil and water resources causes significant damage to the ecological system.In particular, the decrease in the soil quality of agricultural lands leads to reduced productivity and plant diversity (Jin et al., 2019;Dereumeaux et al., 2020;Jiang et al., 2020;Viana et al., 2022).Soil degradation, the main element of the ecological system, threatens sustainable land use.Sustainable agriculture is defined as the rational and effective use of land resources.For this reason, it is necessary to prepare regional plans and product patterns by considering agricultural policies and environmental concerns to ensure the sustainability of agricultural lands.Land suitability assessment is a scientific planning and management model to establish an ecological-economic balance (Musakwa, 2018;Mihoub et al., 2022;Mohammed et al., 2022).This model guides decision makers and producers by evaluating the suitability of the plant pattern to be grown in the field.Land suitability assessment plays an important role in determining the most appropriate use of lands and creating a sustainable environment for agricultural production.Determining the most appropriate use of the land is aimed for the farmers to produce efficiently and protect natural resources (Talukdar et al., 2022).Scientific analysis, classification, and use plans of soils provide the necessary data for the protection of agricultural areas.It plays an important role in the planning of investments such as soil conservation, irrigation, land use planning, and consolidation by determining the use potential of the soils.Determining the qualities of soil and water resources by evaluating them with participatory methods is used in the development of socio-economic, environmental efficiency, and protection methods (Ennaji et al., 2018;Shakya et al., 2021).Qualitative analysis of environmental factors such as soil properties, topography, and climate and quantitative evaluations based on yield estimates are important criteria in land suitability assessment studies.With these studies, the physical, chemical, and biological properties of the soil are evaluated together, and the soil quality index and suitability of the soil are determined.When parameters are evaluated alone, it is impossible to obtain sufficient information about the extent to which soil quality is affected (Askari & Holden, 2015;Malone et al., 2015;Dengiz & Özyazıcı, 2016;Budak et al., 2017;Shakya et al., 2021).The analytical hierarchy process (AHP) method has gained importance in determining soil quality in recent years.In addition, GIS-based multi-criteria decision-making methods (MCDM) using the AHP approach are successfully applied in determining land suitability in different soil and land conditions (Dengiz & Sarıoğlu, 2013;Dağlı& Çağlayan, 2016;Mihoub et al., 2022).In agricultural suitability analysis, it is possible to analyze and visualize environmental and spatial geographical data accurately and up-to-date by integrating geographic information systems (GIS) with the MCDM approach.With these applications, it is possible to produce objective and consistent information through multi-factor criteria.Considering the current situation, it is essential to study the efficient and effective use of limited agricultural lands to ensure sustainable economic growth and development of countries.The study was carried out in the Alaköy and Atmaca areas of Tuşba-Van district in Turkey.The reason for choosing the region is that there are no similar studies, and organic-sustainable agricultural practices are followed.The wheat-barley growing system is the dominant crop pattern; therefore, it was chosen in the scope of the current research.Suitable areas for ecological-economic cultivation were determined for wheat-barley production using the integration of the AHP approach with the GIS technique.The relationship between productivity and economy was evaluated by creating a model.It is thought that labor and economic losses will be prevented when the producers use the lands according to the determined suitability classes.In addition, similar studies to be carried out in different locations and with crop patterns will contribute to sustainable agriculture without damaging ecological resources.

Location
The study was carried out in an area of approximately 4207 ha in the Alaköy and Atmaca districts of Tusba-Van province, which is included in the 1:25,000 scaled K50c1 and K50d2 sheeted topographic maps.The study area is located between the coordinates 4278,000-4288,000 K and 344,000-352,000 D (Zone 38-UTM, m) (Fig. 1).
According to the long-term  climate data of the Van Regional meteorology station, the annual average temperature is 9.8 • C, and the annual average precipitation is 406.8 mm.The highest precipitation is in April with, 54.5 mm, and the lowest precipitation is in July, with 6.8 mm.The study area is ''Mesic'', and the soil moisture regime is ''Xeric' ', according to Soil Survey Staff (1999).The research area includes sandstone, claystone, marble, pebble, and limestone, formed in different geological periods (Cretaceous-Jurassic, Lower Cretaceous, Upper Cretaceous, Pliocene, Upper Paleozoic, Quaternary and Oligocene -Lower Miocene).

Methodology
The topographic structure in the study area does not contain similar features, showing differences even in a short distance.Sampling was done within the framework of  Lindsay & Norvell (1978) geostatistical methods to determine the point values and spatial distributions of the parameters of soil properties in a total area of 4,207 ha.Surface soil (0-30 cm) samples of 260 points were collected in the study area by making 400 × 400 m gridding in order to determine the soil properties.Degraded and undisturbed soil samples were taken from the selected points.GPS was used to identify the coordinates of each sampling point.This study determined the importance levels of soil parameters for wheat-barley production.Physical (slope, depth, sand, clay, silt, hydraulic conductivity, bulk density, and available water capacity), chemical (pH, electrical conductivity, lime, cation exchange capacity, exchangeable sodium percentage, and organic matter), and productivity properties (nitrogen, phosphorus, available potassium, calcium, magnesium, sodium, manganese, iron, copper, and zinc) of 24 soil parameters were examined.These features were divided into subgroups according to their importance for wheat-barley production with expert opinions and literature support (Dedeoglu & Dengiz, 2019;Günal et al., 2022;Kılıc et al., 2022).The methods used to determine the physical, chemical, and fertility properties of soil samples are given in Table 1.
The soil data (physical, chemical, and productivity characteristics) and data processing flow chart used to determine the land suitability and obtain the maps are given in Fig. 2.
Spatial distribution maps of soil properties were prepared using the most common interpolation methods.The methods used are inverse distance weighting (IDW), the radial based function (RBF), and Kriging methods.While methods 1, 2, and 3 are used in the IDW method, the thin plate spline (TPS), completely regularized spline (CRS), and spline with tension (ST) methods are preferred in the RBF method.In the Kriging method, Natural (Ordinary), Simple (Simple), and Universal (Universal) methods are used.The RMSE approach was chosen, and 15 methods were compared while generating interpolation distribution maps.Among these methods, the method that gives the lowest square root (RMSE) of the mean square error was chosen as the most appropriate method.Equation (1) was used to calculate the RMSE value.

RMSE
Z i = estimated value, Z i = measured value and n = number of samples.

Land suitability parameters and soil requirements for wheat-barley cultivation
Wheat can be grown in many soils ranging from sandy loam to clayey texture.The water and nutrient holding capacity of sandy-loamy bodies is less than that of clayey-loamy soils.
Researchers reported that the average nutrient content of wheat soils would be sufficient for the requirements (Sys Ir et al., 1993;El Baroudy, 2016;Karimi et al., 2018).It is desired that the effective soil depth accepted in wheat cultivation should be between 120 and 150 cm.Under optimum conditions, a soil depth of less than 90 cm limits plant root growth (Fischer et al., 2021;Navidi et al., 2022).The fact that the exchangeable sodium percentage of the soil is less than 15% is very important in terms of plant growth and soil quality index.If the exchangeable sodium percentage is high, the soil pH will increase, and the nutrients in the soil will turn into useless forms (Boyer et al., 2008;Pilevar et al., 2020).In terms of other soil components, the threshold values of the reference values determined for wheat cultivation have been reported by different researchers.The desired limit range values of the 22 reported soil parameters are given in Table 2. Barley cultivation is not sensitive compared to wheat.It is cultivated in more barren soils.Barley can be grown in any soil conditions where wheat is grown.Therefore, our study's methodology was created based on the wheat plant, which has effective breeding criteria.Effective selection criteria were determined over the soil factors demanded by wheat.

Scoring of soil parameters and weighting with AHP
Properties related to land and soil quality indices in terms of wheat-barley farming were scored with values between 0 and 1.Standard scoring functions (SSF) given in Table 3, which are widely used in the literature, were used for the scoring process.If a parameter achieves a high score value, SSF indicates a positive relationship (more is better) between soil quality and this parameter, so positive SSF is used.In other cases, negative SSF is used to obtain the desired parameters in good soil quality (less is better) (Wymore, 1993;Karlen & Stott, 1994;Masto et al., 2008;Liu et al., 2018).
Positive functions for depth, AWC, HC, clay, CEC, OM, N, P, Ca, Mg, Fe, Cu, Zn, and Mn were used as soil properties.In addition, a negative scoring function was preferred for Zn (mg kg −1 ) >1 1-0.5 0.5-0.25 <0.25 FAO, 1976Sys Ir et al., 1993Sillanpaa, 1990El Baroudy., 2016Fischer et al., 2021 slope, BD, sand, silt, pH, EC, CaCO 3 , Na, and ESP.Field suitability analysis of wheat-barley was carried out by evaluating multiple criteria.Land suitability assessment was made using Linear Combination Technique by weighting the criteria that affect the agricultural land use of plants according to the order of importance.Each criterion was divided into sub-criteria and scored by giving a numerical value.The sub-criterion scores were multiplied by the weight value of the criteria they belong to, and the criteria were combined on the same scale.Equation ( 2) was used to calculate the most suitable lands for wheat-barley cultivation.
Here, SQI: total land suitability score for wheat-barley cultivation, Wi: weight value of parameter i (weights with AHP), Xi: sub-criteria score of parameter i (from scoring), n: total number of parameters considered.The values calculated using the Linear Combination Technique for each mapping unit in the soil database were classified according to Table 4.A land suitability map was created for wheat-barley crops by considering this classification, soil, and land quality parameters (Gugino et al., 2009).Multi-Criteria Decision Making
In this way, quantitative weights reflecting the importance of each variable are determined (Srdjevic et al., 2010).Therefore, AHP is a method used in studies where many qualitative and quantitative variables are evaluated together (Panchal & Shrivastava, 2022).The AHP method is based on a pairwise comparison matrix where all criteria are compared with each other (Saaty, 2008).The importance values of the variables that make up the matrix used in the study are given in Table 5.
After the sub-parameters are scored and created, the values of the basic parameters emerge.Revealing the fitness map depends on the main parameter values.Each parameter has different effect levels.In order to achieve this effect, geographical parameters should be compared with each other.The comparison of the parameters was carried out according to Saaty (1980) (Table 5).In binary comparison, two parameters are compared with each other, and the comparison depends on the judgment of the decision maker (Öztürk & Batuk, 2010).Pairwise comparison matrix (n) criteria is made with [n(n-1)/2] equation  (Akıncı, Özalp & Turgut, 2012).A normalized comparison matrix is obtained by dividing the column values of the matrix by the sum of each column.The row values obtained from the binary comparisons are added together, and the weight values are obtained by dividing the total value by the value in the row.While making pairwise comparisons of parameters in the AHP method, a certain amount of inconsistency may occur because it is made with individual or simple decisions.Therefore, the obtained priority vector should be checked for consistency (Mezughi et al., 2012;Öztürk & Batuk, 2010).The consistency index (CI) is calculated according to Equation (3).
The CR (Consistency Ratio) formula is used to calculate the consistency ratio of the binary comparison provisions (Saaty, 1994).In order to obtain the Consistency Ratio (CR), it is necessary to know the Random Consistency Index (RI).This index is given as the mean of the constancy index as a result, according to the order of the matrix presented by Saaty (1980) (Table 6).
The consistency ratio (CR) is calculated by dividing the obtained CI by the RI value (Equation ( 4)): The calculated consistency ratio must be 0.10 (10%) or less to be valid.However, if the consistency ratio is greater than 0.10, the pairwise comparison matrices should be revised (Saaty, 1980).These calculations increase the accuracy of the decisions and ensure the reliability of the analysis results (Arslan & Khisty, 2005;Saaty, 2008).

The relationship between soil quality index classes and yield values
The accuracy of the land suitability model, which reveals the suitability of the lands for wheat production, was carried out using the wheat yield of 2020 obtained from the Farmer Registration System and the results of the face-to-face surveys with the breeders.In determining the accuracy of the land suitability map, the relationship between the 64 yield values and the land suitability index score obtained by AHP and GIS methods is determined by regression.The (R 2 ) obtained using Equation ( 5) provides information about the success of the model.R 2 value of 1 or close to 1 gives information about the accuracy of the regression (Zhang et al., 2015).
Here, yi and ŷi denote the observed and predicted soil values for wheat-barley yield, respectively.y represents the average of the wheat-barley yield observation values.

Soil properties and maps
Descriptive statistical properties of the 260 surface soil samples taken from the study area are given in Table 7.The clay content of the surface soils in the study area varied between 10.53 and 60.83%, with an average of 28.62%.Considering the silt content, the lowest value was 10.54%, while the highest value was 44.95%.The sand content of the soils in the study area varies between 10.54 and 76.80%, and the average sand content was 46.73%.The hydraulic conductivity values were 1.57 cm h −1 on average, the average bulk density value was 1.49 g cm −3 , and the available water content was 11.10% on average.The slope values varied between 1 and 30%.The soil depth was found to be between 20 and 120 cm.OM content was determined as 0.54% in wheat-barley cultivation areas and 4.04% in pastures.The lime content of the soils varied between 1.23% and 43.67%, and the average lime content was 11.24%.The pH value ranged from neutral (7.09) to strongly alkaline (9.82), and the average pH value (8.27) was determined to be slightly alkaline.The total nitrogen content varied between 0.036 and 0.540%, with an average value of 0.140%.The available phosphorus content ranged from 0.017 to 24.16 mg kg −1 , with an average phosphorus content of 5.62 mg kg −1 .The exchangeable calcium content of the soils varied between 7.17 and 40.07 cmol kg −1 , exchangeable magnesium content between 0.16 and 10.39 cmol kg −1 , exchangeable potassium content between 0.29 and 4.42 cmol kg −1 , and exchangeable sodium content between 0.15 and 5.31 cmol kg −1 .The Fe, Zn, Cu, and Mn contents of the soils vary between 2.57 and 77.83 mg kg −1 , 0.16 and 5.21 mg kg −1 , 0.04 and 5.65 mg kg −1 , and 4.78 and 172.54 mg kg −1 , respectively.The average available Fe, Zn, Cu, and Mn contents were determined as 11.85, 0.61, 1.34, and 33.95 mg kg −1 , respectively.The purpose of giving descriptive statistics is to emphasize the determined value ranges of the criteria examined, to express the changes within the variation with the coefficient of variation, and to make the criteria examined understandable.The investigated soil properties are important in defining and classifying the soil.It constitutes a selection criterion in the selection of cultivable plant patterns on the identified and classified soil structures and in the suitability of their cultivation.For this reason, the criteria we examine in land suitability studies are accepted as basic research criteria (Bayraklı, Dengiz & Kars, 2023).In determining the suitability of wheat-barley agriculture, the limit values required by the plant regarding the soil are known.Considering the score values of the examined criteria, land suitability for wheat-barley cultivation is determined.
The coefficients of variation of the surface soils of the study area are given in Table 7. Hydraulic conductivity, bulk density, and pH had coefficients of variation less than 15%.Sand, clay, silt, available water capacity, CEC, and exchangeable Ca were observed to be between 15 and 35%.It has been observed that the slope, depth, EC, OM, lime, ESP, N, P, Mg, K, Na, Fe, Zn, Cu, and Mn values of the soils have coefficients of variation over 35%.Sand, bulk density, and available water showed negative skewness to the left, while hydraulic conductivity and slope showed positive skewness to the right.Normal distribution was observed in pH and CEC.Chemical parameters, EC, OM, CaCO 3 , and ESP, showed positive skewness to the right.The soil fertility properties N, P, Mg, K, Na, Fe, Zn, Cu, and Mn showed positive skewness to the right.Exchangeable Ca content showed a normal distribution.Statistically, 70.83% of all soil parameters examined in the study area show a positive distribution skewed to the right.In a left-skewed, negative distribution, the features have a higher distribution frequency than the mean, while the opposite happens in a right-skewed, positive distribution (Pacci, Saflı& Dengiz, 2023).Different interpolation methods were examined in the study.The performances of the models, which have different approaches to spatial distribution, were evaluated.Spatial distribution maps were produced by selecting the highest-performing models.The distribution of SQI was examined using different interpolation methods, and the RMSE values are given in Tables 8, 9 and 10.RMSE values are an important criterion in the selection of models.The smaller the RMSE value, the more accurate the model is (Şenol et al., 2020;Alaboz, Odabas & Dengiz, 2023).Fifteen models were used in our study.For each soil criterion, the model representing the smallest RMSE value was selected, and a map was made for that parameter.For each soil criterion examined, six maps were produced using the completely regularized spline (CRS) method whose RMSE values were at most radial based function (RBF).The RMSE values of each model were determined for the physical, chemical, and productivity parameters of the study area surface soils.The RMSE values of the models are given in Tables 8, 9 and 10.Spatial distribution maps of the parameters used in determining the land suitability of the study area for wheat-barley cultivation are given in Fig. 3.

Impact of criteria on land suitability
The study area consists of different topographic structures and geological units.Therefore, 24 different soil parameters were used to evaluate the suitability of the land.The AHP method was used to determine the relative weights of the parameters used for land suitability assessment in wheat-barley cultivation.The matrix of the pairwise comparisons of physical, chemical, productivity, and main parameters with the AHP method is given in Table 11.The normalized matrix is presented in Table 12.The CR values indicating the consistency of the AHP were calculated as 0.096, 0.090, 0.090, and 0.055, respectively (Table 13).According to the consistency results obtained, it was seen that the AHP matrices were consistent and within acceptable limits.
The weight parameters obtained as a result of the evaluation of soil properties with AHP are given in Table 14.Efficiency parameters (Hierarchy B3) were represented by the lowest core weight value of 0.1958, while the highest core weight value was 0.4934 for physical properties (Hierarchy B1).For each hierarchy, it was seen that OM, slope, and N contributed the most from physical, chemical, and productivity parameters, respectively.In terms of soil fertility and quality, the chemical properties of the soil and the amount of nutrients are at optimum levels.The physical content of the soil is not in suitable conditions.Although the chemical and fertility parameters of the soils are at the desired level, the unsuitability of the physical parameters negatively affects the soil quality index.Therefore, the contribution rate of physical quality parameters has the greatest effect Productivity 0.5 0.5 1 Total 2 3.5 5 (Şenol et al., 2020).Among the determined physical parameters, the most weight value was sloped.The slope factor affects most agricultural activities, irrigation methods, erosion rate, agricultural input use, and soil quality (Mihoub et al., 2022).Therefore, as the percentage slope decreases, it positively affects the productivity of products such as wheat-barley (Seyedmohammadi et al., 2019).While soil organic matter was represented with the highest weight (0.2982) among the parameters, exchangeable sodium (Na) had the lowest weight with a value of 0.0192 (Table 14).Therefore, in terms of productivity, OM is accepted as one of the most important quality parameters in obtaining soil quality index (Andrews, Karlen & Cambardella, 2004).Organic matter is an important indicator of soil fertility and is of great importance for agricultural production and food security   (Obalum et al., 2017;Mihoub et al., 2022).The positive effects of organic matter on the physical and chemical properties of the soil have been reported in many studies (Erdal et al., 2018;Alaboz et al., 2021).The total percent nitrogen content of soils has a very important role in determining soil fertility scores (Dadhich, Patel & Kalubarme, 2017;Günal et al., 2022).The high level of exchangeable Ca and Na contents causes the productivity parameters to have low weight values in the created hierarchy.It has been observed both in our study and in many literatures that the high pH of the soil causes problems in the chemical quality scores of the soil.For the healthy growth of wheat-barley, it is desirable that the soil pH level be neutral or close to neutral.That is, soil pH is important for the healthy development of wheat-barley (Mandal et al., 2020;Pilevar et al., 2020;Mihoub et al., 2022).It is reported that as the exchangeable sodium content of the soil increases, it has negative effects on the physical structure and permeability level of the soil (Boyer et al., 2008;Pilevar et al.,  2020; Sparks, Singh & Siebecker, 2023).Fe affects yield and quality positively (Havlin, 2020).While Cu in the soil is found to be sufficient for plants at low contents, Cu in excessive amounts in the soil has a phytotoxic effect on plants (Michaud et al., 2007).It has been found that sufficient Zn increases the yield and quality of wheat-barley (Hui et al., 2022).Chlorosis (yellowing) was observed between leaf veins in plants with Mn deficiency (McCauley, Jones & Jacobsen, 2009;Kacar & Katkat, 2010).

Land suitability assessment for wheat-barley production
The suitability of the study area for wheat-barley cultivation is presented in Fig. 4. The resulting map was created with the AHP-GIS hybrid approach.Areas suitable for wheatbarley cultivation are presented in Table 15.A percentage of 28.29% (1189.94ha) of the studied area was found to be well suited (S2) for wheat-barley production.Moderately suitable (S3) areas cover 1679.95ha with 39.93%.Weak (N1) suitable areas were found to be 1185.20 ha with a rate of 28.17%.2.29% (96.49ha) of the study area was determined as a settlement, 0.94% (39.62 ha) as a stream bed, 0.21% (8.79 ha) as an irrigation pond and 0.16% (6.83 ha) as an irrigation canal.When the distribution map representing the soil quality index values was examined, it was seen that the regions with the lowest quality (N1) in the study area were located in the northwest and south directions (Table 15).The main reason for this situation is due to the high level of slope, sand, and silt content of the soil in the northwest and south directions.Moderately suitable areas (S3) were spread over the western and eastern edges of the study area.This is due to the fact that the physical parameters of the soils representing the well-suited (S2) category in the study area are represented with the highest weight values.Soil texture significantly affects soil aeration, water holding capacity, nutrient content, soil aggregation, structure development, and tillage (Dharumarajan & Hegde, 2022).The value of 68.22% (S2+S3) of the surveyed land represents the acceptable group in wheat-barley farming.When the land structure was examined, the soils clayey loam, silty-clay loam, sandy-clay-loam, silty, silty-loam, silty-clay, and sandy-clay texture classes were also observed (Table 15).This structure is suitable for optimum development of wheat-barley plants.Soil texture is important in wheat-barley cultivation.For example, loamy soils protect the plant against drought because of their high water holding capacity, while sandy soils provide good aeration and regulate the temperature of the soil.Clay soils are important in terms of soil structure and aggregation.It is also effective in the retention and utilization of plant nutrients (Pilevar et al., 2020;Mandal et al., 2020;Mihoub et al., 2022).
The increase in bulk density due to the increase in the amount of sand in the soil is a negative situation regarding root development (Pacci et al., 2022).The increase in the clay content of the soil and the amount of OM has a positive effect on the water retention in the soil.Furthermore, it has a positive effect on the physical quality index.Both are very important in water conservation.The increase in the available water content of the soil is a significant variable in plant development.The increased amount of available water proves that the plant can benefit from more water (Alaboz, Demir & Dengiz, 2020).
The land suitability assessment has been identified as one of the scientific planning and managerial approaches to promote the sustainability of agricultural production.It is known that many researchers use this model similarly in land suitability assessment studies (Zhang et al., 2015;Ostovari et al., 2019;Orhan, 2021;Talukdar et al., 2022).

Validation of the land suitability for wheat-barley cultivation
In the study area, wheat-barley yield for 2020 varies between 900 and 3,000 kg ha −1 , and the average yield in 2020 was 1,600 kg ha −1 .The relationship between wheat-barley yield grown in Alaköy and Atmaca and SQI was assessed using regression analysis.The regression coefficient value was found to be R 2 = 0.67.This reveals the suitability of the model and shows that the prediction level is successful (Fig. 5).In addition, Sys Ir et al. (1993) reported that grain production of 2,500-3,500 kg ha −1 in rainy conditions and 4,000-6,000 kg ha −1 in irrigation conditions is a good commercial yield.The average yield value was taken since wheat-barley grown in the studied study area was grown under different conditions (wet and dry).It shows that the validation results and the land suitability analysis using AHP and GIS integration were in harmony.Several researchers reported that the AHP method offers more realistic and rational evaluations in the analysis of hierarchical and complex data.At the same time, it is stated that AHP has many more effective and successful advantages than other classical parametric methods (Yalew, Van Griensven & Van der Zaag, 2016;Dedeoglu & Dengiz, 2019;Orhan, 2021).

CONCLUSION
This study reveals that the methods and modeling used to determine land suitability classes are highly usable in production planning.Determining land suitability classes in agricultural lands is considered an important achievement in the effective and sustainable use of land for wheat-barley cultivation.It is seen that it is critical to make land suitability classifications in land and watershed management, planning, and use.Firstly, it is necessary to determine land suitability classes in agricultural lands.Thus, by ensuring the sustainability of land use, it may be possible to protect plant species and diversity by ensuring ecological balance.For this purpose, the effects of the parameters on land suitability classes were weighted using the AHP model, and the standard scoring function was used.The soil quality map created by this method includes three different soil quality index classes as poor, medium, and good.A positive regression (R 2 = 0.67) was found between the soil quality index and the yield values obtained from farmer data records and survey data.
This level of relationship shows the effectiveness of the modeling.According to the suitability analysis, it was determined that 28.29% of the study area is well (S2) suitable, 39.93% moderately (S3) suitable, and 28.17% poorly (N1) suitable for wheat-barley production.These results will benefit farmers in practice and have a positive socio-economic impact.This study will enable the preparation of agriculturally sustainable management plans and the increase of agricultural production by protecting natural resources.At the same time, it is foreseen that it can increase the efficiency of land consolidation studies, product-specific land suitability assessments, and land integration.
Making land suitability classifications is important in determining which land falls into which suitability classes during the production phase, in choosing an effective product pattern, and in economic production.With this study, it was determined which part of the region was suitable for wheat-barley production.It is inevitable that there will be economic losses as a result of insisting on wheat-barley cultivation in areas with low suitability classes.
Land use planning is possible with such studies.The study can be applied as a spatial decision support system in agricultural planning, and effective and economic agricultural planning can be provided.In addition, the obtained parameters are essential for both field and horticultural agriculture.It will be useful in planning the economic cultivation of all cropping patterns.It can be considered a basic study for different cropping patterns other than wheat-barley cultivation.
• Siyami Karaca conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.