Next Article in Journal
Deep Learning-Based Digital Surface Model Reconstruction of ZY-3 Satellite Imagery
Next Article in Special Issue
Impacts of the Middle Route of the South-to-North Water Diversion Project on Land Surface Temperature and Fractional Vegetation Coverage in the Danjiang River Basin
Previous Article in Journal
Histogram-Based Edge Detection for River Coastline Mapping Using UAV-Acquired RGB Imagery
Previous Article in Special Issue
Effects of Soil Moisture and Atmospheric Vapor Pressure Deficit on the Temporal Variability of Productivity in Eurasian Grasslands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Land Suitability Potential Using Ensemble Approaches of Advanced Multi-Criteria Decision Models and Machine Learning for Wheat Cultivation

by
Kamal Nabiollahi
1,2,*,
Ndiye M. Kebonye
2,3,
Fereshteh Molani
1,
Mohammad Hossein Tahari-Mehrjardi
4,
Ruhollah Taghizadeh-Mehrjardi
5,
Hadi Shokati
2 and
Thomas Scholten
2
1
Department of Soil Science and Engineering, College of Agriculture, University of Kurdistan, Sanandaj 6617715175, Iran
2
Chair of Soil Science and Geomorphology, University of Tübingen, 72074 Tübingen, Germany
3
DFG Cluster of Excellence “Machine Learning: New Perspectives for Science”, University of Tübingen, AI Research Building, Maria-von-Linden-Str. 6, 72076 Tübingen, Germany
4
Faculty of Management, University of Tehran, Tehran 141556311, Iran
5
Faculty of Agriculture and Natural Resources, Ardakan University, Ardakan 77871-31587, Iran
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(14), 2566; https://doi.org/10.3390/rs16142566
Submission received: 6 June 2024 / Revised: 8 July 2024 / Accepted: 10 July 2024 / Published: 12 July 2024
(This article belongs to the Special Issue Mapping Essential Elements of Agricultural Land Using Remote Sensing)

Abstract

:
Land suitability assessment, as an important process in modern agriculture, involves the evaluation of numerous aspects such as soil properties, climate, relief, hydrology and socio-economic aspects. The aim of this study was to evaluate the suitability of soils for wheat cultivation in the Gavshan region, Iran, as the country is facing the task of becoming self-sufficient in wheat. Various methods were used to evaluate the land, such as multi-criteria decision-making (MCDM), which is proving to be important for land use planning. MCDM and machine learning (ML) are useful for decision-making processes because they use complicated spatial data and methods that are widely available. Using a geomorphological map, seventy soil profiles were selected and described, and ten soil properties and wheat yields were determined. Three MCDM approaches, including the technique of preference ordering by similarity to the ideal solution (TOPSIS), gray relational analysis (GRA), and simple additive weighting (SAW), were used and evaluated. The criteria weights were extracted using Shannon’s entropy method. Random forest (RF) model and auxiliary variables (remote sensing data, terrain data, and geomorphological maps) were used to represent the land suitability values. Spatial autocorrelation analysis as a statistical method was applied to analyze the spatial variability of the spatial data. Slope, CEC (cation exchange capacity), and OC (organic carbon) were the most important factors for wheat cultivation. The spatial autocorrelation between the key criteria (slope, CEC, and OC) and wheat yield confirmed these results. These results also showed a significant correlation between the land suitability values of TOPSIS, GRA, and SAW and wheat yield (0.74, 0.72, and 0.57, respectively). The spatial distribution of land suitability values showed that the areas classified as good according to TOPSIS and GRA were larger than those classified as moderate and weak according to the SAW approach. These results were also confirmed by the autocorrelation of the MCDM techniques with wheat yield. In addition, the RF model showed its effectiveness in processing complex spatial data and improved the accuracy of land suitability assessment. In this study, by integrating advanced MCDM techniques and ML, an applicable land evaluation approach for wheat cultivation was proposed, which can improve the accuracy of land suitability and be useful for considering sustainability principles in land management.

1. Introduction

One of the most important issues facing humanity is food and health security [1]. The growing world population requires more food, while land resources are decreasing due to population growth and some destructive environmental problems such as climate change, land degradation, and urbanization. Therefore, concern about land capacity is increasing [2].
In many developing countries, land managers use the land for different crops without paying attention to the possibilities and principles of sustainability [3,4]. The unplanned changes in land use have led to destructive phenomena such as deforestation, desertification, global warming, and soil erosion. Therefore, soils in developing countries should be utilized according to their potential in order to increase yield production and conserve soil resources.
Land suitability is determined by the land use planning process [5], and land evaluation is an important step in this process. Land suitability assessment requires a potential evaluation of environmental characteristics such as climate, soil, relief, water, social, and economic factors [6]. Through a detailed land suitability assessment, the most suitable crops or land uses for specific regions can be identified [7,8,9,10].
Until 1973, the systems introduced to classify land suitability or capability were generally used to evaluate land for different uses [11]. Some of these methods have been applied in industrialized countries, but the accessible information on these land uses is sometimes unrelated to regional knowledge and conditions [12,13]. In recent decades, research has focused on the use of multi-attribute models for complex decision-making related to land suitability. Existing multi-criteria decision-making (MCDM) and geographic information system (GIS) techniques allow the combination of different information to manage and design applications.
MCDM, which was developed in the 1960s, is widely used in various fields [14]. MCDM is a framework to generate the best choice for complicated circumstances and helps developers and policymakers select the best solutions based on their constraints [15,16,17]. Determining the priority of effective factors in evaluation and decision-making is always considered one of the basic principles of planning and management. In these decisions, multiple measurement attributes can be used instead of one attribute. This method involves a series of techniques that allow a number of independent attributes to be evaluated, weighted, and then ranked.
MCDM methods have a great potential to reduce costs and increase accuracy in spatial decisions, as well as to provide an applicable solution for problem-solving [16,18]. In this context, techniques such as the technique of preference ordering by similarity to the ideal solution (TOPSIS), gray relational analysis (GRA), and simple additive weighting (SAW) are commonly used [19,20,21]. These methods are a good tool for solving challenging problems with complicated decision conditions in a variety of sectors.
On the other hand, machine learning (ML) methods have been increasingly developed as part of the knowledge of artificial intelligence in the age of technology and information in various scientific fields, especially in digital mapping [22,23,24]. One of the main applications of ML methods is the identification and prediction of existing patterns in several large datasets, such as satellite data, parameters derived from digital elevation models, spectroscopy data, or other sources [25]. Specifically, in the area of land suitability, several studies have successfully applied ML techniques to assess and predict land suitability for various purposes. For example, Ahmed and Hussain [26] applied ML to predict wheat production in northern Pakistan. Their study shows the potential of machine learning in improving agricultural forecasting for developing countries like Pakistan. Similarly, Naghdyzadegan Jahromi [27] used ML models to predict wheat yields based on environmental variables. Their models showed promising results.
The integration of ML algorithms and MCDM methods as an applicable and useful approach has recently been used by researchers in various fields [28,29,30], although there is little research in this regard in the field of land evaluation [28,29]. In determining site suitability, a combination of ML and MCDM can be used to evaluate various parameters, make appropriate land use decisions, and increase the accuracy of the assessment. Saha and Mondal [28] evaluated the suitability of agricultural land in India using MCDM and machine learning techniques. They used the Analytic Hierarchy Process for weighting the criteria and compared the effectiveness of methods such as Random Forest (RF) and FCOPRAS for analyzing land suitability. RF demonstrates high model accuracy through ROC curves and accuracy measurements.
Wheat has the largest acreage of any cereal in the world and, as one of the most important agricultural crops, has an impact on global food security and human nutrition. It is an important source of food for a large part of the world’s population. It provides food for more than 33 percent of the world’s population and is considered a strategic commodity. Iran is one of the largest wheat producers in the Middle East region [30]. However, wheat cultivation in Iran faces difficulties as the lands have not been tested for their suitability, and the plan for “food independence” in Iran has been postponed [31]. Therefore, proper land management is crucial for crop planning and requires a close examination of actual characteristics such as soil quality, relief, water, and climate. It can help to increase wheat yield and minimize undue impacts on the environment. Few studies have been conducted to compare advanced MCDA techniques for land suitability assessment. MCDA methods are particularly considered due to their robust and effective algorithms [32]. Furthermore, ML approaches were combined with the MCDA approach for land suitability mapping of wheat to increase the accuracy of the assessment. The aim of this study was to (1) determine the suitability of land for wheat cultivation in western Iran by combining MCDA and ML, (2) compare the advanced MCDA techniques (TOPSIS, SAW, and GRA), and (3) determine the most important factors for land suitability for wheat cultivation in the studied area.

2. Materials and Methods

2.1. Studied Area Description

The research area, whose size is estimated at 5341 hectares, is in the Gavshan region in western Iran (Figure 1). The area is categorized as a cold, semi-arid climate (BWk) according to the Köppen–Geiger climate classification system. Winters are characterized by cold temperatures and heavy rainfall, while summers are quite dry and hot. The average annual precipitation is 370 mm, with more than 90% of it falling between mid-autumn and the end of winter. The average annual temperature is rather cool at 10.8 °C (Figure 2). The soil temperature and humidity regimes are described as Mesic and Xeric [33], respectively. The minimum and maximum altitudes are 1520 and 2366 m above sea level, respectively. The most common soils in the study area belong to the soil reference groups Leptosols, Lithosols, Cambisols, and Calcisols according to the FAO World Reference Base for Soil Resources (WRB 2014) [34] as well as Entisols and Inceptisols according to the USDA soil taxonomy [33]. Cropland (about 80% of the total area) and pastureland are the two main land use categories in this region. Geomorphological features of the area include piedmonts, alluvial fans, hills, and mountains with gentle to steep slopes (Figure 1 and Table 1).

2.2. Soil and Yield Sample Collection

The methodological framework applied to this study is shown in Figure 3. Considering that different parameters have an influence on the spatial distribution of soil, a map of geomorphology was created based on the four-level hierarchical method developed by Tommanian [35]. Then, according to the layers of the geomorphology map at the landform level, the locations of 70 soil profiles were determined using a stratified random sampling method. They were then excavated and described using standard methods and relevant soil properties such as soil structure, soil color, flooding, and drainage characteristics [36]. Soil samples were carefully collected from all soil horizons and then transported to the laboratory. In the next step, the organic carbon content was quantified by a wet combustion method [37]. Soil pH and electrical conductivity (EC) were determined by measurements in a saturated paste using a pH electrode [38] and a conductivity meter [39], respectively. Cation exchange capacity (CEC) was determined using the 1-N-ammonium acetate method at a pH of 7.0 [40]. The calcium carbonate equivalent (CCE) was determined volumetrically [41]. The particle size distribution was determined using the Bouyoucos hydrometer method [42]. The percentage of exchangeable sodium (ESP) was calculated as the ratio of sodium to CEC. In addition, the available phosphorus (P) was determined with the extraction agent Bray-1, a solution of 0.03 M ammonium fluoride (NH4F) and 0.025 M hydrochloric acid (HCl) [43], the available potassium (K) with 1 M ammonium acetate (NH4OAc) as an extraction agent [44], and the gravel content. The ratio of sodium to CEC was used to calculate the percentage of exchangeable sodium in the soil. The gypsum content was considered negligible and, therefore, excluded from further investigation. In addition, slope and climate data from the Gavshan synoptic weather station [45] and a digital elevation model [46] were used.

2.3. Yield Samples Collection

Wheat samples were systematically taken from each soil profile within a 1 m2 to ensure a uniform plant density within the defined area on 23 August 2019. All plant samples were harvested manually and dried at a controlled temperature of 75 °C for 48 h, and the wheat grains were then weighed. The wheat yield was determined for each of the 70 sampling points. The yield (Y) is calculated using the following formula (Equation (1)):
Y = N 1 × 10,000 × N 2 × N 3 1000 1000
Y stands for tons per hectare (t ha−1), N1 for the number of ears per square meter, N2 for the number of grains per ear, which are the components of wheat yield, and N3 for the weight of 1000 grains (in grams).
As a rule, the use of machinery with conventional tillage systems (moldboard plow and disk harrow) [47] on wheat farms ends between 15 October. Wheat was planted on 15 October. In the study area, burning of crop residues is a regular activity. Before sowing, 50 kg ha−1 of triple superphosphate fertilizer was applied to the wheat fields. Before sowing and during the growth of the stalks, 100 kg ha−1 of nitrogen fertilizer (urea) was applied to the wheat fields. The usual herbicides and insecticides are 24D and deltamethrin, respectively. Similar management parameters are used throughout the research region, including machinery, crop rotation, planting density, sowing technique, plant spacing, fertilizer rate, and winter wheat variety.

2.4. Land Suitability Assessment

2.4.1. Diagnostic Attribute

When evaluating land, the first step, which has to do with care, is the selection of the associated attributes. In this study, the suitability of soils for wheat cultivation was assessed on the basis of the following attributes: soil texture, gravel content, soil depth, gypsum, CEC, CCE, ESP (exchange sodium percentage), pH, EC, available phosphorus, available potassium, slope, climatic conditions, susceptibility to flooding, and drainage characteristics. Finally, following the principles for land suitability evaluation formulated by Sys [11], eight attributes were selected for the evaluation of wheat croplands based on a review of the literature and a detailed study of the spatial variability of the parameters along the studied area [12,48,49,50,51,52]. It is worth mentioning that our direct investigations in the studied area showed that this relevance exists. Moreover, these variables proved to be the main constraints for wheat cultivation and are closely related to the specific conditions, soil properties, and soil requirements of the study area.

2.4.2. Diagnostic Attributes Weights

In this phase, the weights of the individual diagnostic attributes were determined using Shannon’s entropy. This method provides a suitable framework for calculating the attribute weights. It was developed by Claude Shannon [53] and is often used to solve the challenge of load determination in MCDM. It assesses the degree of uncertainty or variation within a group of attributes and facilitates the assignment of relative loads to each attribute. In subjective assessment, such as the entropy weighting technique, the subjective weighting is taken directly from the actual data, while the subjective weighting is usually determined by an expert’s judgment based on experience. The advantage of the entropy weighting method is that it increases objectivity while reducing the subjective influence of the decision-makers. Originally, the entropy concept was used to calculate the disorder of a system in the field of thermodynamics. Shannon used it as a strategy to deal with uncertainty by applying it to problems in information theory. Entropy theory states that information can be supplied in proportion to the entropy value. As a result, more weight can be given to the criterion [54]. The idea of entropy weighting has been widely applied in a variety of fields [55,56]. In the following, Shannon’s entropy weight calculation is introduced.
Suppose the decision problem has m options (A1, A2, …, Am) and n criteria (C1, C2, …, Cn).
The first decision matrix is as follows:
Step 1: Standardize the decision matrix.
A = a 11 a 12 a 1 n a 21 a 22 a 2 n a m 1 a m 2 a m n = a i j m × n
where aij donate ith options of jth criterion.
The standardized value of the jth attribute for the ith answer, Ai is given by rij.
r i j = a i j i = 1 m a i j 2 i = 1,2 , , m ;
Step 2: Calculate the entropy.
e j = K i = 1 m r i j l n r i j ,   j = 1,2 , , n ;
where K = 1/ln m.
Step 3: The weights of each attribute are calculated.
r i j = 1 e j i = 1 n ( 1 e j ) i = 1,2 , , n

2.4.3. Rating Values for Selected Attributes

The rating values for the specific parameters of wheat were determined according to the methods described in Table 2. These values were formulated by compiling information from scientific references [11,52]. These scores were then used in the MCDA methods to produce a land suitability rating for wheat cultivation.

2.4.4. An Overview of the Proposed Method

Three MCDM methods, namely SAW, TOPSIS, and GRA, were used to determine the land potential for wheat cultivation. Using these methods, each soil profile was ranked by determining the land suitability based on the weights and scores of the diagnostic attributes.
SAW: In this approach, decision-makers directly weigh each attribute to express its relative importance. The assigned importance for each attribute is then multiplied by the standardized value for wheat cultivation for that attribute to obtain an overall score. The resulting products are then aggregated across all attributes for each soil profile [57,58]. The following methods (3 and 4, respectively) are used to determine the adjusted value for both positive and negative attributes.
n i j = g i j g m a x       i = 1 , . , m       j = 1 , , n
n i j = g m i n g i j       i = 1 , , m       j = 1 , , n
F i n a l   s c o r e = w g i j × n i j   w g i j = 1
where gij is the value of the criterion, gmax is the highest possible score for each positive requirement, gmin is the lowest possible score for each negative requirement, wg is the weighting, and nij is the normalized value.
TOPSIS: TOPSIS is an MCDM method for evaluating and prioritizing choices regarding attributes depending on their distance to the best positive and negative answer. This method was proposed by Hwang and Yoon [59] and soon found its place in the field of MCDM. The term TOPSIS refers to preference methods based on equality for the best answer. The logic of the method is the definition of the best answer. The best answer is the one that increases the advantage of the attribute and decreases the disadvantage of the attribute. The best choice is the one that has the smallest distance to the best answer and, at the same time, the furthest distance to the negative best answer [57,58]. The TOPSIS process comprises the following phases.
Step 1: Create a decision matrix.
r i j = x 11 x 12 x 1 n x 21 x 22 x 2 n x m 1 x m 2 x m n
Step 2: Create a standardized decision matrix.
The standardized value of the jth attribute for the ith answer, Ai is given by rij.
r i j = x i j i = 1 m x i j 2 i = 1 , 2 , , m ;   j = 1 , 2 , n
Step 3: The weighted and standardized decision matrix is calculated.
v i j = w j r i j ,   i = 1 , , m ;   j = 1 , , n
where wj is the weight of the jth attribute.
Step 4: Determine the negative and the best positive answer.
A = v 1 , , v n
A + = v 1 + , , v n +
where A denotes the best negative answer and A+ is the best positive answer. If the jth attribute is an advantageous attribute, then
v j + = m a x v i j , i = 1 , , m , o r v j +   a n d   v j = m i n v i j , i = 1 , , m ,
If, on the other hand, the jth attribute is a cost attribute, then
v j + = m i n v i j , i = 1 , , m ,   o r   v j +   a n d   v j = m a x v i j , i = 1 , , m ,
Step 5: The distances between each answer and the best positive and negative answer are calculated.
d i + = j = 1 n ( v i j v j + ) 2   i = 1 , , m
d i = j = 1 n ( v i j v j ) 2   j = 1 , , n
di+ and di donate the distance between the ith answer and the best positive and negative answer, respectively.
Step 6: The relative estimate of the best answer is calculated.
C i d i d i + + d i ;   0 C i 1 ;   i = 1 , 2 , , m
Step 7: Ranking preference orders are done based on the Ci values; the choices are ranked in descending order.
GRA: This method has shown that it has effectiveness in various challenges in MCDM [60,61]. GRA is founded on the fundamental concepts of the grey systems theory. GRA principles are founded on the grey systems theory and theory that is applied to solve difficult problems [58]. GRA is used to assess the correlation between two series of alternatives, and references use grey relational grades. The procedure for GRA includes the following computational steps.
Steps 1–4 are like Steps 1–4 of the TOPSIS section.
Step 5: The grey relationship coefficients are calculated using the weighted standardized decision matrix. The grey relationship coefficient represents the relationship between the ith answer and the best positive answer for the jth attribute:
γ i j + = γ r i j , v j + = m + ξ M + Δ i j + ξ M + , ξ 0 , 1
Δ i j + ξ M + = r i j v j + , m + = m i n i m i n j Δ i j + , m + = m a x i m a x j Δ i j + , ξ is the coefficient of discrimination.
Using the matrix of grey relationship coefficients of the best positive answer, the grey relationship coefficients between different choices are created.
γ + = γ i j + m × n
By comparing with the negative best answer, the grey relationship coefficient for the ith answer of the jth attribute can be calculated as follows:
γ i j = γ r i j , v j = m ξ M Δ i j ξ M , ξ 0 , 1
where Δ i j ξ M = r i j v j , m ± = m i n i m i n j Δ i j ,   a n d   m ± = m a x i m a x j Δ i j .
If you consider the negative best answer, the grey coefficient matrix is shown here:
γ = γ i j m × n
Step 6: the grey relational grade is determined.
To create an overall evaluation of the solutions based on all attributes, the grey relational score is used. It is defined as the mean of the relational coefficients for the distinctive attributes. For the ith difference, the grey relational score is expressed in terms of the best positive difference as follows:
G i + = 1 n j = 1 n γ i j +
Similarly, the degree of the grey relation for the ith difference in relation to the negative difference is as follows:
G i = 1 n j = 1 n γ i j
Step 7: The relative grey value of the relationship is calculated.
C i D i D i + + D i
Step 8: the alternatives are arranged in descending order based on their Ci.

2.4.5. Comparison and Validation of Methods for Determining Land Suitability Values

In this step, the results of the previous step (i.e., the determination of land suitability values for 70 soil profiles by three MCDM methods, including TOPSIS, SAW, and GRA) were evaluated. In evaluating the effectiveness of the MCDM methods, the analysis focused primarily on identifying the notable correlations between the values calculated by each method and the observed wheat yield. We used the statistical metric of coefficient of determination (R2) to evaluate how well the methods predicted the value of land suitability.

2.4.6. Mapping and Modeling of Land Suitability Values

To assess the spatial variations of land suitability values in our study area, we used a random forest (RF) ML, digital soil mapping (DSM) and a set of environmental covariates. Within the DSM, environmental covariates are used to predict soil attributes and the relationships between soil attributes and environmental covariates are established by ML models. RF is a great machine learning model whose predictive ability has been demonstrated in various fields [62,63,64]. A literature review on RF can be found in Antoniadis [65] and Belgiu and Dragut [66].
The variables were a combination of categorical data (e.g., a geomorphological map), terrain features generated from a model of digital elevation with a resolution of 10 m, and remote sensing data obtained from Landsat imagery with a resolution of 30 m. Together, these data sources contributed to the predictive modeling of suitability values for each geographic unit.
SAGA software 9.5.1 [67] was used to calculate terrain features such as valley floor flatness, slope, LS factor, topographic wetness, aspect, plane curvature, and elevation. The spectral bands and some indices, such as the normalized difference vegetation index, were created from a Landsat 8 Operational Land Imager (OLI) image (23 August 2019). All environmental covariates were resampled to achieve a uniform grid size (30 m) and were then used to predict the determined land suitability values by three MCDM methods.
All data were divided into categories, including train data (70%) and test data (30%). To validate the results of the random forest ML models, a rigorous evaluation strategy was applied. A k-fold cross-validation procedure with 100-fold replication was applied. The statistical indices, including coefficient of determination (R2), mean absolute error (MAE), and root mean square error (RMSE), were applied.

2.4.7. Classification of the Land Suitability Value

To allow a better visual comparison of the map produced by the MCDM methods, the land suitability values calculated by each method were categorized into three classes: high, medium and low. This categorization process included the following steps:
-
Determine the range: The range of land suitability values for each method has been determined. This range is the difference between the maximum and minimum suitability values obtained with each method.
-
Divide into intervals: The total range of values was divided into three equal intervals. These intervals correspond to the categories of land suitability: high, medium, and low.
-
Define interval boundaries: The first interval (low suitability) starts at the minimum score and extends to the value that marks the end of the first third of the range. The second interval (medium suitability) starts at the end of the first interval and extends to the end of the second third of the range. The third interval (high suitability) starts at the end of the second interval and extends to the maximum score.
-
Assign categories: each land suitability value was then assigned to one of the three categories based on the interval in which it fell [68,69].

2.4.8. Spatial Autocorrelation Analysis

We used spatial autocorrelation analysis to further analyze the spatial variability of key attributes and to create a map of land suitability values obtained by the three proposed MCDM approaches (i.e., TOPSIS, SAW, and GRA). A statistical technique used to measure and evaluate the degree of dependency between spatial data is the spatial autocorrelation analysis. In other words, it facilitates the delineation of whether the values are randomly distributed or whether they are close to each other in space [70]. R 4.3.2 software was used to calculate the focal correlation between two raster maps of land suitability assessment produced by the three proposed MCDM approaches and wheat yield.

3. Results

3.1. The Statistical Description of Soil Properties

Table 3 contains a statistical description of the soil properties. The amounts of CCE range from 2.37% to 38%, with a mean of 14.3%; the pH values range from 7 to 7.9, with a mean of 7.52; the amounts of OC range from 0.15 to 1.6%, with a mean of 0.81%; gravel values range from 7.47 to 47.42%, with a mean of 27.54%; EC values range from 0.2 to 1.5 dS m−1, with a mean of 0.59 dS m−1; and CEC values range from 2.17 to 23.63 Cmol + kg−1, with a mean of 11.49 Cmol + kg−1. The main soil textural classes are SC (sandy clay), C (clay), SCl (sandy clay loam), CL (clay loam), SL (sandy clay), and L (loam). EC, CCE, sand, and clay have a low positive skewness. pH and silt, on the other hand, have a negative skewness. The coefficient of variation (CV) for EC, CEC, OC, gravel, sand, and CCE is greater than 30% (Table 3). In the study area, silt and clay have a moderate coefficient of variation (10–30%), while pH has a low coefficient of variation (10%).

3.2. Calculated Weighting of the Attributes

Because each attribute has a specific and significant effect on plant growth, decision-makers must assign different importance to these requirements when evaluating land suitability. Table 4 shows the results of using Shannon’s entropy to determine the final weights for all attributes. The results show that slope, OC, and CEC are the most significant variables in evaluating the suitability of land for wheat cultivation, with loadings of 0.439, 0.135, and 0.108, respectively (Table 4). The loadings for CCE, EC, gravel, soil thickness, and soil texture are 0.090, 0.086, 0.077, 0.054, and 0.010, respectively. The spatial variability of key attributes and wheat yield in the studied area is shown in Figure 4.

3.3. Determination and Validation of the Land Suitability Value

Table 5 shows the results of the MCDA methods used to determine the land suitability value of the study region for wheat cultivation. This table also shows the results of the categorized land suitability value. With a mean of 0.31, TOPSIS results range from 0.06 to 0.51, GRA range from 0.14 to 0.44, with a mean of 0.27, and SAW range from 0.06 to 0.21, with a mean of 0.15, and were used to determine the land suitability values.
As can be seen in Figure 5, the performance of each model is evaluated by determining the significant relationship with wheat yield. There are differences between the MCDA approaches, as our comparison shows. The TOPSIS approach gives the best results with a value of 0.74, as shown in Figure 6, followed by GRA (0.72) and SAW.

3.4. Land Suitability Value and Key Attribute Mapping

We used the RF algorithm and a DSM method to assign the suitability values determined by the MCDA approaches. Table 6 shows the result of the RF algorithm for assigning the suitability of wheat cropland using three MCDM methods. The SAW approach (0.003, 0.061, and 0.88) gives the best results as measured by the MAE, RMSE, and R2 statistical indices. This is followed by the GRA (0.081, 0.017 and 0.83) and TOPSIS (0.005, 0.072, and 0.80) methods. Figure 6 shows the results of this study, including maps of machine learning-based land suitability values for wheat cultivation in the study area using the RF model.

3.5. Spatial Autocorrelation Analysis

Figure 7a shows the results of a spatial autocorrelation study comparing wheat yield with three important criteria: slope, CEC, and OC. Spatial autocorrelation can be broadly categorized into three types: positive, negative, and absent. The data from the study area show that slope, CEC, and OC correlate strongly with wheat yield. Within the study area, slope inclination and yield show a negative and significant spatial autocorrelation. On the other hand, the spatial autocorrelation between yield and OC and CEC is positive and high in the study area (Table 2, Figure 3 and Figure 7a). In addition, Figure 7b shows the map of spatial variability of land suitability assessment as determined by the three proposed MCDM approaches (i.e., TOPSIS, SAW, and GRA). In the study area, the spatial autocorrelation of the MCDA methods with the yields is also positive, with the TOPSIS method showing the strongest correlation.

4. Discussion

4.1. Key Attributes for Wheat Cultivation

The results show that soil properties, especially CEC and OC, have undeniable effects on soil productivity and the promotion of agricultural sustainability in the study area (Table 4 and Figure 4). The coefficients of variation (CV) of slope, OC, and CEC (106.75, 49.95, and 49.02, respectively) (Table 2) support this assertion. This study emphasizes the critical importance of topographic factors in developing a reliable land evaluation model, with slope having a particularly strong influence on yield outcomes. As an important regulator of water distribution, slope can directly influence drainage conditions, erosion and water retention in the landscape. On gentle slopes, it is more important for agriculture as it controls the infiltration rate and prevents waterlogging. On steep slopes, however, it can lead to soil erosion and a reduction in soil productivity [71,72].
As an important soil property, CEC has a strong influence on soil productivity and, thus, on crop yields. OC was identified as another important soil property in the study area that significantly influences soil productivity via pH and soil fertilization. It is known that the content of OC in the soil has a positive influence on soil productivity [73]. This variation in OC, CEC, and slope was largely due to variability in topography and parent material. This shows that topography and parent material are the most important factors for soil formation in the study area, as Taghizadeh-Mehrjardi [12] also found.
Dedeoğlu and Dengiz [74] calculated the land suitability index for wheat by integrating AHP and GIS in Turkey, and their results also showed that the most effective attributes were soil strength, soil texture, and slope indicators.
These results are also confirmed by evaluating the spatial autocorrelation of MCDA methods in terms of crop yields (Figure 7a). It is evident that slope, CEC, and OC have a high correlation with crop yields. The spatial autocorrelation of slope and yield within the study area is remarkably negative and high and is observed in most of the region. There is a positive spatial autocorrelation in low slope areas, which include Pi111, Pi121, Pi211, Pi212, Pl111, and pl211. In contrast, spatial autocorrelation is negative and high in high slope areas, which mainly include Hi111, Hi121, Hi131, Hi211, Mo111, and Mo121. Regarding OC, we observed different trends. In areas with high OC and CEC content, there is a positive correlation between yield and OC (Figure 7a).

4.2. Evaluation of the MCDM Techniques

Depending on the validation results (Figure 5), the best results were observed for the TOPSIS method, followed by the GRA and SAW methods. The results show that the TOPSIS and GRA methods have a stronger significant correlation with yield than the SAW method, indicating that the TOPSIS and GRA methods efficiently determined the land suitability assessment in most of the studied areas. The methods (TOPSIS and GRA) have also been used in the context of land suitability by other researchers, such as Ekmekcioglu [75] and Zhang and Duan [76]. These studies have shown that these techniques are a useful tool for decision makers in evaluating land suitability.
The methods help decision-makers make long-term land use and management decisions by considering various parameters [77]. The differences between these methods (TOPSIS and GRA) and SAW lie in the assignment of weights, the use of different mathematical functions, and standardization. Although the SAW method is easy to use and understand, it may not capture the differences in land suitability assessment as well as the other two methods. Seyedmohammadi [57] also evaluated the application of three models, including the fuzzy TOPSIS, TOPSIS, and SAW models, to determine the land suitability values for the cultivation of three crops (soybean, canola, and corn) in the northern part of Iran. Their results proved that the fuzzy TOPSIS method was the best approach for determining land suitability. These results are further supported by examining the spatial autocorrelation of MCDA methods with respect to crop yields (Figure 7b). The spatial autocorrelation of MCDA methods with yields within the study area is positive and high and can be found throughout the region. However, the spatial autocorrelation is lower in areas with a lower number of samples, especially in the northern and western parts.

4.3. Mapping and Modeling the Land Suitability Values

The result of the RF algorithm for mapping the land suitability of wheat is acceptable (Table 6). Random forest provides the best estimate by integrating the results of many decision trees. This method can model complicated spatial patterns and is well suited for mapping land suitability levels. The RF algorithm has also been successfully used in other land suitability studies [12,77].
According to the maps of land suitability for wheat cultivation (Figure 6), suitable land units (good class) are mainly concentrated in the center of the study area. These areas are characterized by a slope of less than 5%, a yield of more than 1.4 t ha−1, a CEC of more than 15 cmol + kg−1 and an OC of more than 1% (Figure 4 and Table 5). These areas account for 33%, 20%, and 6% of the suitability values assessed using the TOPSIS, GRA, and SAW methods, respectively. They correspond to the geomorphological land units Pi111, Pi211, and Pl111 (Figure 1). Moderately suitable land (medium class) is characterized by a slope of 5% to 10%, a yield in the range of 1–1.4 t ha−1, a CEC in the range of 10–15 cmol + kg−1, and an OC of more than 0.5 to 1% (Figure 4 and Table 5). These land categories are concentrated in the central region and partly extend to the southern and eastern parts of this study (Figure 6). They account for 16%, 22%, and 29% of the suitability values determined with the TOPSIS, GRA, and SAW methods, respectively. They are assigned to the geomorphological units Pi121, Pi212, and Pl211 (Figure 1). Overall, the TOPSIS and GRA methods classify more areas as good compared with the SAW method, which classifies more areas as moderately or poorly suitable.

4.4. Practical Applications

Overall, this study offers practical advice with a variety of implications for land use planning and agricultural management. The results of the study have a major impact on the establishment of land use regulations in the area studied. A solid basis for optimal decisions is the understanding and identification of the different soil classes suitable for wheat cultivation. Finding suitable land with certain topographic and soil characteristics, often found in the center of the study area, can provide a great opportunity for wheat cultivation, and this tactic will likely increase wheat yields.
However, information on areas that are neither well nor moderately adapted helps in decision-making, e.g., when it comes to nature conservation or other land uses. Precision agriculture approaches tailored to local conditions are supported by research by identifying critical variables such as soil properties (CEC and OC) and topographic parameters (slope). Precision agriculture builds on this expertise to help farmers manage their land more prudently and reduce the risk of soil degradation. Fertilization and irrigation techniques are adapted based on soil fertility and the ability to store water. This study’s conclusions benefit sustainable land use planning by identifying techniques to preserve the environment in environmentally fragile areas and optimize crop selection to increase productivity. Policymakers, investors and farmers can make informed decisions about agricultural investments, policies and climate change adaptation strategies with the help of decision support tools developed based on land suitability maps. In the future, the use of cutting-edge technologies such as continuous monitoring and remote sensing will further improve the accuracy and relevance of land suitability assessments.

5. Conclusions

In this study, site suitability was determined by advanced MCDM techniques and machine learning for wheat cultivation in the Gavshan region of western Iran. The effects of various factors (such as soil and environment) on wheat yield were evaluated, and the importance of sustainable land management for optimizing agricultural productivity was considered. In this study, three MCDM techniques, including TOPSIS, SAW, and GRA, were presented for land suitability evaluation and their application for wheat cultivation was discussed. The weighting of diagnostic attributes was determined using Shannon’s entropy. It was found that factors such as slope, CEC, and OC were the most important attributes and had the greatest influence on the evaluation of wheat land suitability. The results proved that the TOPSIS and GRA methods were more efficient than the SAW method in determining land suitability and had a higher correlation with the observed wheat yield data. The spatial distribution of suitability values was mapped with RF and categorized into three levels, including good, moderate, and poor areas for wheat cultivation. In general, this study shows the importance of combining advanced MCDM techniques and ML to assess the suitability of wheat areas. It underlines the need for a comprehensive analysis of soil and environmental attributes to optimize agricultural productivity and sustainability. The results provide important insights for land use policy and agricultural management decisions and contribute to more informed and sustainable land use practices. In general, this study has shown the importance of combining advanced MCDM techniques and ML to assess the suitability of wheat land. It underlines the need for a comprehensive analysis of soil and environmental characteristics to optimize agricultural productivity and sustainability. The results provide important insights for decision-making in land use policy and agricultural management and contribute to more informed and sustainable land use practices. In other words, even though this study provides insightful results for wheat cultivation, there are some limitations that need to be considered. First, the models were calibrated using local data, so their direct transferability to other locations with different conditions is limited. The results are unique to the Gavshan region in western Iran. Furthermore, this study is not applicable to temporal variations in environmental conditions that may have a significant impact on wheat cultivation, such as seasonal variations or long-term climatic changes. Finally, additional validation with independent datasets from other locations or time periods could improve the generalizability of the results, even if the MCDM and machine learning models were validated with local wheat production data.

Author Contributions

Conceptualization, K.N.; methodology, K.N., F.M., N.M.K. and R.T.-M.; software, K.N., M.H.T.-M. and N.M.K.; validation, K.N. and R.T.-M.; formal analysis, K.N.; investigation, K.N., F.M. and N.M.K.; writing—K.N., F.M. and H.S.; visualization, K.N., H.S. and N.M.K.; supervision, K.N. and T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the University of Kurdistan and the University of Tubingen. Kamal Nabiollahi was supported by the Alexander von Humboldt Foundation under grant number Ref 3.4-1230353-IRN-GF-E.

Data Availability Statement

Data are available by reasonable email request to the corresponding author at [email protected]. The data are not publicly available due to ethical reasons.

Conflicts of Interest

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.

References

  1. Pozza, L.E.; Field, D.J. The Science of Soil Security and Food Security. Soil Secur. 2020, 1, 100002. [Google Scholar] [CrossRef]
  2. Keating, B.A.; Herrero, M.; Carberry, P.S.; Gardner, J.; Cole, M.B. Food Wedges: Framing the Global Food Demand and Supply Challenge towards 2050. Glob. Food Secur. 2014, 3, 125–132. [Google Scholar] [CrossRef]
  3. Rahmani, M.; Lotfata, A.; Zebardast, E.; Rastegar, S.; Sanchez, T.W.; Goharrizi, B.A.; Landi, S. Land Use Suitability Assessment for Economic Development at the Provincial Level: The Case Study of Yazd Province, Iran. Sustain. Cities Soc. 2022, 87, 104163. [Google Scholar] [CrossRef]
  4. Salem, M.; Tsurusaki, N.; Divigalpitiya, P. Remote Sensing-Based Detection of Agricultural Land Losses around Greater Cairo since the Egyptian Revolution of 2011. Land Use Policy 2020, 97, 104744. [Google Scholar] [CrossRef]
  5. Rossiter, D.G. A Theoretical Framework for Land Evaluation. Geoderma 1996, 72, 165–190. [Google Scholar] [CrossRef]
  6. Abbasi, M.; Pishvaee, M.S.; Bairamzadeh, S. Land Suitability Assessment for Paulownia Cultivation Using Combined GIS and Z-Number DEA: A Case Study. Comput. Electron. Agric. 2020, 176, 105666. [Google Scholar] [CrossRef]
  7. Akpoti, K.; Kabo-bah, A.T.; Zwart, S.J. Review—Agricultural Land Suitability Analysis: State-of-the-Art and Outlooks for Integration of Climate Change Analysis. Agric. Syst. 2019, 173, 172–208. [Google Scholar] [CrossRef]
  8. Dornik, A.; Cheţan, M.A.; Drăguţ, L.; Iliuţă, A.; Dicu, D.D. Importance of the Mapping Unit on the Land Suitability Assessment for Agriculture. Comput. Electron. Agric. 2022, 201, 107305. [Google Scholar] [CrossRef]
  9. Keson, J.; Silalertruksa, T.; Gheewala, S.H. Land Suitability Class and Implications to Land-Water-Food Nexus: A Case of Rice Cultivation in Thailand. Energy Nexus 2023, 10, 100205. [Google Scholar] [CrossRef]
  10. Xue, L.; Cao, P.; Xu, D.; Guo, Y.; Wang, Q.; Zheng, X.; Han, R.; You, A. Agricultural Land Suitability Analysis for an Integrated Rice–Crayfish Culture Using a Fuzzy AHP and GIS in Central China. Ecol. Indic. 2023, 148, 109837. [Google Scholar] [CrossRef]
  11. Sys, C.; Van Ranst, E.; Debaveye, J. Land Evaluation. Part I: Principles in Land Evaluation and Crop Production Calculations; Agricultural Publications Nr. 7, G.A.D.C: Brussels, Belgium, 1991. [Google Scholar]
  12. Taghizadeh-Mehrjardi, R.; Nabiollahi, K.; Rasoli, L.; Kerry, R.; Scholten, T. Land Suitability Assessment and Agricultural Production Sustainability Using Machine Learning Models. Agronomy 2020, 10, 573. [Google Scholar] [CrossRef]
  13. Bagheri Bodaghabadi, M.; Martínez-Casasnovas, J.A.; Khakili, P.; Masihabadi, M.H.; Gandomkar, A. Assessment of the FAO Traditional Land Evaluation Methods, A Case Study: Iranian Land Classification Method. Soil Use Manag. 2015, 31, 384–396. [Google Scholar] [CrossRef]
  14. Kumar, A.; Sah, B.; Singh, A.R.; Deng, Y.; He, X.; Kumar, P.; Bansal, R.C. A Review of Multi Criteria Decision Making (MCDM) towards Sustainable Renewable Energy Development. Renew. Sustain. Energy Rev. 2017, 69, 596–609. [Google Scholar] [CrossRef]
  15. Moisa, M.B.; Tiye, F.S.; Dejene, I.N.; Gemeda, D.O. Land Suitability Analysis for Maize Production Using Geospatial Technologies in the Didessa Watershed, Ethiopia. Artif. Intell. Agric. 2022, 6, 34–46. [Google Scholar] [CrossRef]
  16. Osei-Gyabaah, A.P.; Antwi, M.; Addo, S.; Osei, P. Land Suitability Analysis for Cocoa (Theobroma cacao) Production in the Sunyani Municipality, Bono Region, Ghana. Smart Agric. Technol. 2023, 5, 100262. [Google Scholar] [CrossRef]
  17. Sadiq, F.K.; Ya’u, S.L.; Aliyu, J.; Maniyunda, L.M. Evaluation of Land Suitability for Soybean Production Using GIS-Based Multi-Criteria Approach in Kudan Local Government Area of Kaduna State Nigeria. Environ. Sustain. Indic. 2023, 20, 100297. [Google Scholar] [CrossRef]
  18. Cano, P.B.; Cabrini, S.M.; Peper, A.M.; Poggio, S.L. Multi-Criteria Assessment of Cropping Systems for the Sustainable Intensification in the Pampas. Agric. Syst. 2023, 210, 103723. [Google Scholar] [CrossRef]
  19. Foroozesh, F.; Monavari, S.M.; Salmanmahiny, A.; Robati, M.; Rahimi, R. Assessment of Sustainable Urban Development Based on a Hybrid Decision-Making Approach: Group Fuzzy BWM, AHP, and TOPSIS–GIS. Sustain. Cities Soc. 2022, 76, 103402. [Google Scholar] [CrossRef]
  20. Hooshangi, N.; Mahdizadeh Gharakhanlou, N.; Ghaffari Razin, S.R. Evaluation of Potential Sites in Iran to Localize Solar Farms Using a GIS-Based Fermatean Fuzzy TOPSIS. J. Clean. Prod. 2023, 384, 135481. [Google Scholar] [CrossRef]
  21. Pophali, G.R.; Chelani, A.B.; Dhodapkar, R.S. Optimal Selection of Full Scale Tannery Effluent Treatment Alternative Using Integrated AHP and GRA Approach. Expert Syst. Appl. 2011, 38, 10889–10895. [Google Scholar] [CrossRef]
  22. Taghizadeh-Mehrjardi, R.; Sheikhpour, R.; Zeraatpisheh, M.; Amirian-Chakan, A.; Toomanian, N.; Kerry, R.; Scholten, T. Semi-Supervised Learning for the Spatial Extrapolation of Soil Information. Geoderma 2022, 426, 116094. [Google Scholar] [CrossRef]
  23. Sujatha, M.; Jaidhar, C.D. Machine Learning-Based Approaches to Enhance the Soil Fertility—A Review. Expert Syst. Appl. 2024, 240, 122557. [Google Scholar] [CrossRef]
  24. Yousheng, D.; Keqin, Z.; Zhongju, F.; Wen, Z.; Xinjun, Z.; Huiling, Z. Machine Learning Based Prediction Model for the Pile Bearing Capacity of Saline Soils in Cold Regions. Structures 2024, 59, 105735. [Google Scholar] [CrossRef]
  25. Shahabi, A.; Nabiollahi, K.; Davari, M.; Zeraatpisheh, M.; Heung, B.; Scholten, T.; Taghizadeh-Mehrjardi, R. Spatial Prediction of Soil Properties through Hybridized Random Forest Model and Combination of Reflectance Spectroscopy and Environmental Covariates. Geocarto Int. 2022, 37, 18172–18195. [Google Scholar] [CrossRef]
  26. Ahmed, M.U.; Hussain, I. Prediction of Wheat Production Using Machine Learning Algorithms in Northern Areas of Pakistan. Telecommun. Policy 2022, 46, 102370. [Google Scholar] [CrossRef]
  27. Naghdyzadegan Jahromi, M.; Zand-Parsa, S.; Razzaghi, F.; Jamshidi, S.; Didari, S.; Doosthosseini, A.; Pourghasemi, H.R. Developing Machine Learning Models for Wheat Yield Prediction Using Ground-Based Data, Satellite-Based Actual Evapotranspiration and Vegetation Indices. Eur. J. Agron. 2023, 146, 126820. [Google Scholar] [CrossRef]
  28. Saha, S.; Mondal, P. Estimation of the Effectiveness of Multi-Criteria Decision Analysis and Machine Learning Approaches for Agricultural Land Capability in Gangarampur Subdivision, Eastern India. Artif. Intell. Geosci. 2022, 3, 179–191. [Google Scholar] [CrossRef]
  29. Sarkar, D.; Saha, S.; Maitra, M.; Mondal, P. Site Suitability for Aromatic Rice Cultivation by Integrating Geo-Spatial and Machine Learning Algorithms in Kaliyaganj C.D. Block, India. Artif. Intell. Geosci. 2021, 2, 179–191. [Google Scholar] [CrossRef]
  30. Pakrooh, P.; Kamal, M.A. Modeling the Potential Impacts of Climate Change on Wheat Yield in Iran: Evidence from National and Provincial Data Analysis. Ecol. Model. 2023, 486, 110513. [Google Scholar] [CrossRef]
  31. Mesgaran, M.B.; Madani, K.; Hashemi, H.; Azadi, P. Iran’s Land Suitability for Agriculture. Sci. Rep. 2017, 7, 7670. [Google Scholar] [CrossRef]
  32. Malczewski, J. GIS-based Multicriteria Decision Analysis: A Survey of the Literature. Int. J. Geogr. Inf. Sci. 2006, 20, 703–726. [Google Scholar] [CrossRef]
  33. Keys to Soil Taxonomy, 12th ed.; United Soil Survey Staff: Washington, DC, USA, 2014.
  34. World Reference Base for Soil Resources 2014. In International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; Food and Agriculture Organization of the United Nations: Rome, Italy, 2014; ISBN 978-92-5-108369-7.
  35. Toomanian, N.; Jalalian, A.; Khademi, H.; Eghbal, M.K.; Papritz, A. Pedodiversity and Pedogenesis in Zayandeh-Rud Valley, Central Iran. Geomorphology 2006, 81, 376–393. [Google Scholar] [CrossRef]
  36. Schoeneberger, P.J.; Wysocki, D.A.; Benham, E.C. Field Book for Describing and Sampling Soils; Government Printing Office: Washington, DC, USA, 2012; ISBN 978-0-16-091542-0. [Google Scholar]
  37. Nelson, D.W.; Sommers, L.E. Total Carbon, Organic Carbon, and Organic Matter. In Methods of Soil Analysis; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 1982; pp. 539–579. ISBN 978-0-89118-977-0. [Google Scholar]
  38. Mclean, E.O. Soil pH and Lime Requirement. In Methods of Soil Analysis; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 1982; pp. 199–224. ISBN 978-0-89118-977-0. [Google Scholar]
  39. Rhoades, J.D. Soluble Salts. In Methods of Soil Analysis; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 1982; pp. 167–179. ISBN 978-0-89118-977-0. [Google Scholar]
  40. Buchan, G.D.; Grewal, K.S.; Claydon, J.J.; Mcpherson, R.J. A Ccomparison of Sedigraph and Pipette Methods for Soil Particle-Size Analysis. Soil Res. 1993, 31, 407–417. [Google Scholar] [CrossRef]
  41. Sparks, D.L.; Page, A.L.; Helmke, P.A.; Loeppert, R.H. Methods of Soil Analysis, Part 3: Chemical Methods; John Wiley & Sons: Hoboken, NJ, USA, 2020; ISBN 978-0-89118-825-4. [Google Scholar]
  42. Gee, G.W.; Or, D. 2.4 Particle-Size Analysis. In Methods of Soil Analysis; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2002; pp. 255–293. ISBN 978-0-89118-893-3. [Google Scholar]
  43. Carter, M.R.; Gregorich, E.G. (Eds.) Soil Sampling and Methods of Analysis, 2nd ed.; Canadian Society of Soil Science: Pinawa, MB, Canada; CRC Press: Boca Raton, FL, USA, 2008; ISBN 978-0-8493-3586-0. [Google Scholar]
  44. Brown, J.G.; Lilleland, O. Rapid Determination of Potassium and Sodium in Plant Materials and Soil Extracts by Flame Photometry. Am. Soc. Hortic. Sci. 1946, 48, 341–345. [Google Scholar]
  45. Gavshan Synoptic Weather Station Data Iran Meteorological Organization. Available online: http://www.kurdistanmet.ir/fa/Reading/%D8%A7%DB%8C%D8%B3%D8%AA%DA%AF%D8%A7%D9%87%D9%87%D8%A7%DB%8C-%D8%A8%D8%A7%D8%B1%D8%A7%D9%86-%D8%B3%D9%86%D8%AC%DB%8C (accessed on 1 July 2024).
  46. niafam.com. National Cartographic Center of Iran. Available online: https://en.ncc.gov.ir/Services (accessed on 1 July 2024).
  47. Javadi, A.; Rahmati, M.H.; Tabatabaeefar, A. Sustainable Tillage Methods for Irrigated Wheat Production in Different Regions of Iran. Soil Tillage Res. 2009, 104, 143–149. [Google Scholar] [CrossRef]
  48. Pilevar, A.R.; Matinfar, H.R.; Sohrabi, A.; Sarmadian, F. Integrated Fuzzy, AHP and GIS Techniques for Land Suitability Assessment in Semi-Arid Regions for Wheat and Maize Farming. Ecol. Indic. 2020, 110, 105887. [Google Scholar] [CrossRef]
  49. Bagherzadeh, A.; Gholizadeh, A. Modeling Land Suitability Evaluation for Wheat Production by Parametric and TOPSIS Approaches Using GIS, Northeast of Iran. Model. Earth Syst. Environ. 2016, 2, 126. [Google Scholar] [CrossRef]
  50. Bagherzadeh, A.; Mansouri Daneshvar, M.R. Qualitative Land Suitability Evaluation for Wheat and Barley Crops in Khorasan-Razavi Province, Northeast of Iran. Agric. Res. 2014, 3, 155–164. [Google Scholar] [CrossRef]
  51. Karimi, F.; Sultana, S.; Shirzadi Babakan, A.; Royall, D. Land Suitability Evaluation for Organic Agriculture of Wheat Using GIS and Multicriteria Analysis. Pap. Appl. Geogr. 2018, 4, 326–342. [Google Scholar] [CrossRef]
  52. Givi, J. Qualitative Evaluation of Land Suitability for Field and Fruit Crops. Available online: http://www.swri.ir/en-US/DouranPortal/5223/page/Recent-Publications (accessed on 1 July 2024).
  53. Shannon, C.E. A Mathematical Theory of Communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  54. Lotfi, F.H.; Fallahnejad, R. Imprecise Shannon’s Entropy and Multi Attribute Decision Making. Entropy 2010, 12, 53–62. [Google Scholar] [CrossRef]
  55. Lee, H.-C.; Chang, C.-T. Comparative Analysis of MCDM Methods for Ranking Renewable Energy Sources in Taiwan. Renew. Sustain. Energy Rev. 2018, 92, 883–896. [Google Scholar] [CrossRef]
  56. Hafezalkotob, A.; Hafezalkotob, A. Extended MULTIMOORA Method Based on Shannon Entropy Weight for Materials Selection. J. Ind. Eng. Int. 2016, 12, 1–13. [Google Scholar] [CrossRef]
  57. Seyedmohammadi, J.; Sarmadian, F.; Jafarzadeh, A.A.; Ghorbani, M.A.; Shahbazi, F. Application of SAW, TOPSIS and Fuzzy TOPSIS Models in Cultivation Priority Planning for Maize, Rapeseed and Soybean Crops. Geoderma 2018, 310, 178–190. [Google Scholar] [CrossRef]
  58. Wang, P.; Zhu, Z.; Wang, Y. A Novel Hybrid MCDM Model Combining the SAW, TOPSIS and GRA Methods Based on Experimental Design. Inf. Sci. 2016, 345, 27–45. [Google Scholar] [CrossRef]
  59. Hwang, C.-L.; Yoon, K. Multiple Attribute Decision Making: Methods and Applications A State-of-the-Art Survey; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; ISBN 978-3-642-48318-9. [Google Scholar]
  60. Cao, W.; Xuan, X.; Wu, W.; Chen, L.; Zhang, Q. A New Method for Determining the Optimal Adaptive Thermal Comfort Model with Grey Relational Analysis. Build. Environ. 2022, 221, 109277. [Google Scholar] [CrossRef]
  61. Yuce, B.E.; Nielsen, P.V.; Wargocki, P. The Use of Taguchi, ANOVA, and GRA Methods to Optimize CFD Analyses of Ventilation Performance in Buildings. Build. Environ. 2022, 225, 109587. [Google Scholar] [CrossRef]
  62. Nabiollahi, K.; Taghizadeh-Mehrjardi, R.; Shahabi, A.; Heung, B.; Amirian-Chakan, A.; Davari, M.; Scholten, T. Assessing Agricultural Salt-Affected Land Using Digital Soil Mapping and Hybridized Random Forests. Geoderma 2021, 385, 114858. [Google Scholar] [CrossRef]
  63. Azizi, K.; Ayoubi, S.; Nabiollahi, K.; Garosi, Y.; Gislum, R. Predicting Heavy Metal Contents by Applying Machine Learning Approaches and Environmental Covariates in West of Iran. J. Geochem. Explor. 2022, 233, 106921. [Google Scholar] [CrossRef]
  64. Poppiel, R.R.; Demattê, J.A.M.; Rosin, N.A.; Campos, L.R.; Tayebi, M.; Bonfatti, B.R.; Ayoubi, S.; Tajik, S.; Afshar, F.A.; Jafari, A.; et al. High Resolution Middle Eastern Soil Attributes Mapping via Open Data and Cloud Computing. Geoderma 2021, 385, 114890. [Google Scholar] [CrossRef]
  65. Antoniadis, A.; Lambert-Lacroix, S.; Poggi, J.-M. Random Forests for Global Sensitivity Analysis: A Selective Review. Reliab. Eng. Syst. Saf. 2021, 206, 107312. [Google Scholar] [CrossRef]
  66. Belgiu, M.; Drăguţ, L. Random Forest in Remote Sensing: A Review of Applications and Future Directions. ISPRS J. Photogramm. Remote Sens. 2016, 114, 24–31. [Google Scholar] [CrossRef]
  67. Olaya, V.; Conrad, O. Chapter 12 Geomorphometry in SAGA. In Developments in Soil Science; Hengl, T., Reuter, H.I., Eds.; Geomorphometry, Elsevier: Amsterdam, The Netherlands, 2009; Volume 33, pp. 293–308. [Google Scholar]
  68. Sánchez-Navarro, A.; Gil-Vázquez, J.M.; Delgado-Iniesta, M.J.; Marín-Sanleandro, P.; Blanco-Bernardeau, A.; Ortiz-Silla, R. Establishing an Index and Identification of Limiting Parameters for Characterizing Soil Quality in Mediterranean Ecosystems. Catena 2015, 131, 35–45. [Google Scholar] [CrossRef]
  69. Nabiollahi, K.; Taghizadeh-Mehrjardi, R.; Kerry, R.; Moradian, S. Assessment of Soil Quality Indices for Salt-Affected Agricultural Land in Kurdistan Province, Iran. Ecol. Indic. 2017, 83, 482–494. [Google Scholar] [CrossRef]
  70. Griffith, D.A. Spatial Autocorrelation: A Primer; Resource publications in geography; Association of American Geographers: Washington, DC, USA, 1987; ISBN 978-0-89291-197-4. [Google Scholar]
  71. Liu, S.; Wang, R.; Yang, Y.; Shi, W.; Jiang, K.; Jia, L.; Zhang, F.; Liu, X.; Ma, L.; Li, C.; et al. Changes in Soil Aggregate Stability and Aggregate-Associated Carbon under Different Slope Positions in a Karst Region of Southwest China. Sci. Total Environ. 2024, 928, 172534. [Google Scholar] [CrossRef] [PubMed]
  72. Nabiollahi, K.; Golmohamadi, F.; Taghizadeh-Mehrjardi, R.; Kerry, R.; Davari, M. Assessing the Effects of Slope Gradient and Land Use Change on Soil Quality Degradation through Digital Mapping of Soil Quality Indices and Soil Loss Rate. Geoderma 2018, 318, 16–28. [Google Scholar] [CrossRef]
  73. Ajami, M.; Heidari, A.; Khormali, F.; Gorji, M.; Ayoubi, S. Environmental Factors Controlling Soil Organic Carbon Storage in Loess Soils of a Subhumid Region, Northern Iran. Geoderma 2016, 281, 1–10. [Google Scholar] [CrossRef]
  74. Dedeoğlu, M.; Dengiz, O. Generating of Land Suitability Index for Wheat with Hybrid System Aproach Using AHP and GIS. Comput. Electron. Agric. 2019, 167, 105062. [Google Scholar] [CrossRef]
  75. Ekmekcioğlu, Ö.; Koc, K.; Özger, M. Stakeholder Perceptions in Flood Risk Assessment: A Hybrid Fuzzy AHP-TOPSIS Approach for Istanbul, Turkey. Int. J. Disaster Risk Reduct. 2021, 60, 102327. [Google Scholar] [CrossRef]
  76. Zhang, X.; Duan, X. Evaluating Water Resource Carrying Capacity in Pearl River-West River Economic Belt Based on Portfolio Weights and GRA-TOPSIS-CCDM. Ecol. Indic. 2024, 161, 111942. [Google Scholar] [CrossRef]
  77. Wang, Y.; Li, Y. Mapping the Ratoon Rice Suitability Region in China Using Random Forest and Recursive Feature Elimination Modeling. Field Crops Res. 2023, 301, 109016. [Google Scholar] [CrossRef]
Figure 1. Location map of the studied area in a province in Iran (Refer to Table 1).
Figure 1. Location map of the studied area in a province in Iran (Refer to Table 1).
Remotesensing 16 02566 g001
Figure 2. Water climate diagrams for the study area.
Figure 2. Water climate diagrams for the study area.
Remotesensing 16 02566 g002
Figure 3. The methodological framework of this study.
Figure 3. The methodological framework of this study.
Remotesensing 16 02566 g003
Figure 4. Maps of the spatial variability of key attributes in the study area.
Figure 4. Maps of the spatial variability of key attributes in the study area.
Remotesensing 16 02566 g004
Figure 5. Linear relationships between wheat yield and land suitability values calculated by three MCDM methods.
Figure 5. Linear relationships between wheat yield and land suitability values calculated by three MCDM methods.
Remotesensing 16 02566 g005
Figure 6. ML-based maps of land suitability values and wheat classes for three MCDM methods (note: the maps in the first row represent the model predictions, while the maps in the bottom row are the reclassified versions to represent suitability).
Figure 6. ML-based maps of land suitability values and wheat classes for three MCDM methods (note: the maps in the first row represent the model predictions, while the maps in the bottom row are the reclassified versions to represent suitability).
Remotesensing 16 02566 g006
Figure 7. Autocorrelation of the map of wheat yields with the maps of the key criteria (a) and the values of land suitability (b).
Figure 7. Autocorrelation of the map of wheat yields with the maps of the key criteria (a) and the values of land suitability (b).
Remotesensing 16 02566 g007
Table 1. Geomorphology map legend.
Table 1. Geomorphology map legend.
CodeLandscapeLandformLithologyGeomorphological Surface
Hi111HillHighly eroded hillShale and sandstoneHill with high topography
Hi211 Moderately eroded hillLimestone and sandstoneHill with high topography
Mo111MountainRock outcropLimestoneEroded rock surface
Mo121 Sandstone and conglomerate Eroded rock surface
Pi111PiedmontAlluvial fanAlluvial fanModerate to high slop
Pi121 Moderate slope
Pi211 PedimentMedium-level alluvial depositsLow slope
Pi212 Moderate slope
Pl111PlainFlat Recent alluviumLow slope
Pl211 Alluvial terraceMedium-level alluvial depositsCultivated (moderate to few slopes)
Table 2. Soil and topography demands for wheat farming [11].
Table 2. Soil and topography demands for wheat farming [11].
PropertyRating Scale
100–9595–8585–6060–4040–2525–0
Slope (%)0–22–55–88–1616–25>25
EC (dS m−1)0–44–88–1212–1616–24>24
Soil textureSiCL, SiC, SiL, CSCL, LSC, CLSL, Cm, C, SiCmLS S
Gravel (%)0–33–1515–3535–55->55
CCE (%)3–20<320–3535–5050–60 >60
OC (%)>0.60.6–0.4<0.4---
CEC (Cmol + kg−1)24>24–1616>---
Soil depth (cm)>100100–6030–6020–30 <20
SiCL: silty clay loam; SiC: silty clay; SiL: silty loam; SCL: sandy loam; m: massive (compact); LS: loamy sand; S: sandy.
Table 3. Descriptive statistics for soil properties.
Table 3. Descriptive statistics for soil properties.
SandClay SiltOCGravelCCEECCECpH
%dS m−1cmol + kg−1
Number707070707070707070
Mean42.5228.7828.670.8127.5414.30.5911.497.52
Minimum19.9212.009.200.157.472.370.202.177.00
Maximum66.1945.0043.001.6047.4238.001.5023.637.90
Skewness0.170.20−0.290.200.490.771.300.38−0.53
Kurtosis−0.29−0.141.11−0.27−1.090.822.35−0.271.24
CV44.1721.0224.7238.5140.4849.9541.0541.022.15
Table 4. Results of attribute weighting using Shannon’s entropy method.
Table 4. Results of attribute weighting using Shannon’s entropy method.
SlopeOCCECTextureECCCEGravelSoil Thickness
Weight0.4390.1350.1080.0100.0860.0900.0770.054
Table 5. Values for land suitability, land suitability class, and key attributes for wheat cultivation.
Table 5. Values for land suitability, land suitability class, and key attributes for wheat cultivation.
TOPSISSAWGRASlope
%
OC
%
CEC
(cmol + kg−1)
Yield
(ton ha−1)
Mean0.310.150.27150.8111.491.14
Maximum0.510.210.44791.6023.631.75
Minimum0.060.0601400.152.170.5
Good class range0.34–0.440.38–0.520.34–0.44<51<151.4<
Medium class range0.24–0.340.24–0.380.24–0.345–100.5–110–151–1.4
Weak class range0.14–0.240.10–0.240.14–0.245<<0.5<100.5–1
Table 6. Results of RF models for mapping land suitability values for wheat using three MCDM methods.
Table 6. Results of RF models for mapping land suitability values for wheat using three MCDM methods.
GRATOPSISSAW
RMSE0.0720.0810.061
MAE0.0050.0170.003
R20.830.800.88
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nabiollahi, K.; M. Kebonye, N.; Molani, F.; Tahari-Mehrjardi, M.H.; Taghizadeh-Mehrjardi, R.; Shokati, H.; Scholten, T. Assessment of Land Suitability Potential Using Ensemble Approaches of Advanced Multi-Criteria Decision Models and Machine Learning for Wheat Cultivation. Remote Sens. 2024, 16, 2566. https://doi.org/10.3390/rs16142566

AMA Style

Nabiollahi K, M. Kebonye N, Molani F, Tahari-Mehrjardi MH, Taghizadeh-Mehrjardi R, Shokati H, Scholten T. Assessment of Land Suitability Potential Using Ensemble Approaches of Advanced Multi-Criteria Decision Models and Machine Learning for Wheat Cultivation. Remote Sensing. 2024; 16(14):2566. https://doi.org/10.3390/rs16142566

Chicago/Turabian Style

Nabiollahi, Kamal, Ndiye M. Kebonye, Fereshteh Molani, Mohammad Hossein Tahari-Mehrjardi, Ruhollah Taghizadeh-Mehrjardi, Hadi Shokati, and Thomas Scholten. 2024. "Assessment of Land Suitability Potential Using Ensemble Approaches of Advanced Multi-Criteria Decision Models and Machine Learning for Wheat Cultivation" Remote Sensing 16, no. 14: 2566. https://doi.org/10.3390/rs16142566

APA Style

Nabiollahi, K., M. Kebonye, N., Molani, F., Tahari-Mehrjardi, M. H., Taghizadeh-Mehrjardi, R., Shokati, H., & Scholten, T. (2024). Assessment of Land Suitability Potential Using Ensemble Approaches of Advanced Multi-Criteria Decision Models and Machine Learning for Wheat Cultivation. Remote Sensing, 16(14), 2566. https://doi.org/10.3390/rs16142566

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop