Assessing the Influence of Soil Quality on Rainfed Wheat Yield

Soil quality assessment based on crop yields and identification of key indicators of it can be used for better management of agricultural production. In the current research, the weighted additive soil quality index (SQIw), factor analysis (FA), and multiple linear regression (MLR) are used to assess the soil quality of rainfed winter wheat fields with two soil orders on 53.20 km2 of agricultural land in western Iran. A total of 18 soil quality indicators were determined for 100 soil samples (0–20 cm depth) from two soil orders (Inceptisols and Entisols). The soil properties measured were: pH, soil texture, organic carbon (OC), cation exchange capacity (CEC), electrical conductivity (EC), soil microbial respiration (SMR), carbonate calcium equivalent (CCE), soil porosity (SP), bulk density (BD), exchangeable sodium percentage (ESP), mean weight diameter (MWD), available potassium (AK), total nitrogen (TN), available phosphorus (AP), available Fe (AFe), available Zn (AZn), available Mn (AMn), and available Cu (ACu). Wheat grain yield for all of the 100 sampling sites was also gathered. The SQIw was calculated using two weighting methods (FA and MLR) and maps were created using a digital soil mapping framework. The soil indicators determined for the minimum data set (MDS) were AK, clay, CEC, AP, SMR, and sand. The correlation between the MLR weighting technique (SQIw-M) and the rainfed wheat yield (r = 0.62) was slightly larger than that the correlation of yield with the FA weighted technique (SQIw-F) (r = 0.58). Results showed that the means of both SQIw-M and SQIw-F and rainfed wheat yield for Inceptisols were higher than for Entisols, although these differences were not statistically significant. Both SQIw-M and SQIw-F showed that areas with Entisols had lower proportions of good soil quality grades (Grades I and II), and higher proportions of poor soil quality grades (Grades IV and V) compared to Inceptisols. Based on these results, soil type must be considered for soil quality assessment in future studies to maintain and enhance soil quality and sustainable production. The overall soil quality of the study region was of poor and moderate grades. To improve soil quality, it is therefore recommended that effective practices such as the implementation of scientifically integrated nutrient management involving the combined use of organic and inorganic fertilizers in rainfed wheat fields should be promoted. Agriculture 2020, 10, 469; doi:10.3390/agriculture10100469 www.mdpi.com/journal/agriculture Agriculture 2020, 10, 469 2 of 18


Introduction
Winter wheat, the second main cereal in the world [1][2][3], is annually grown on over 200 million ha [4]. Rainfed winter wheat is the main crop of semi-arid regions where mean precipitation is 350-500 mm [4]. The production of rainfed winter wheat per unit area is low. One of the main causes for this low production is that the most suitable soil types for planting have not been identified. Evaluation of soil quality under diverse soil types is vital for identifying problematic areas and assessing sustainable soil management [5][6][7][8]. To overcome this problem and also to increase wheat yield, information on soil quality is needed [9].
The quality of soil is measured through several soil indicators, including chemical, physical, and biological features [10] rather than individual soil indicators. The interpretation of the latter's response to ecological change is difficult [11]. Quantifying soil quality can, therefore, improve understanding of soil ecosystems and inform decision-making [12]. The soil management assessment framework (SMAF) as one of the conceptual frameworks to evaluate soil quality has been designed to make a quantitative evaluation of soil quality [9,13]. A soil quality index (SQI) is calculated by applying SMAF scoring algorithms to soil indicator properties. This is the most frequently used process for assessing soil quality. Soil quality indicators generate SQIs in an easy and flexible format based on an integrated valuation of SQIs and their weights [14][15][16][17]. SQIs have been successfully applied at various scales for different purposes such as assessing the impacts of agricultural practices, soil management, salinization, and alkalization on soil quality and crop production [18][19][20][21][22][23].
Minimum data set (MDS) and factor analysis (FA) approaches have been extensively applied for indicator selection in soil quality assessment [24][25][26][27] as they can decrease the time-consuming, labor-intensive and costly processes involved in measuring extra soil properties [13]. Multiple linear regression analysis (MLR) has been applied to assign weights to soil indicators [6]. If the soil indicator weights are assigned based on MLR, the soil indicators that affect the crop yield will have the largest weight. However, MLR has only rarely been used for weighting quality measures [6,28]. Therefore, one aim of this work is to confirm whether MLR is a useful weighting method. Crop yield is commonly used to assess soil quality [29] and could be beneficial for determining key individual soil quality indicators [30]. Although crop yield has been neglected in many previous studies, in sustainable agricultural systems, high soil quality is characterized by high productivity [31]. The development of crops and their production depend on the soil type. On shallow soils (e.g., Entisols) with near-surface lithic contact, which is a root limiting layer with a high degree of erosion, and low fertility, crop growth is limited and crop expansion will decrease. By comparison, cultivated crops can grow vigorously on deeper soils (e.g., Inceptisols), which have greater soil depth, permeability, water retention, fertility, and a lower degree of erosion [5,15,32,33]. Therefore, a systematic soil quality assessment for rainfed winter wheat, as one of the most important crops throughout the world, based on yield under various soil types is essential to develop appropriate soil management strategies and to enhance the rainfed winter wheat yield. This paper applies a soil quality assessment framework to perform soil quality evaluations for rainfed winter wheat for two different soil orders. In so doing it will address several aims: i-to evaluate the quality of soils based on the two weighting methods: multi linear regression and factor analysis; ii-to establish a minimum data set for assessing soil quality for rainfed winter wheat; iii-to evaluate the impacts of soil type on soil quality and rainfed winter wheat yield; iv-to assess soil quality in rainfed winter wheat systems based on rainfed winter wheat yields.

Description of the Study Area
The study area had an area of 53.20 km 2 and was located in the city of Ghorveh, western Iran ( Figure 1). The climate is semi-arid with a mean annual rainfall of 369.8 mm, mainly concentrated in the winter (October to March). The yearly mean temperature (10.8 • C) is relatively cool. The area is categorized by a Xeric soil moisture regime and a Mesic soil temperature regime. Elevation changes from 1833 to 2627 m and slopes range from 1 to 30% (Figure 1). The primary land use is cropland, where farm machinery operations with conventional tillage systems are commonly completed by 15 October. Winter wheat (Triticum aestivum L.) was planted on 15 October at a depth of 4 cm, 17 cm row spacing and a seed rate of 160 kg ha −1 . Burning crop residue is a common practice in the study area. Wheat and fallow is a common 2-year crop rotation. Triple superphosphate fertilizer was applied at a rate of 50 kg ha −1 , on wheat fields before sowing. Nitrogen fertilizer (urea) was applied at a rate of 100 kg ha −1 , on wheat fields before sowing, and at the stem elongation stage. Commonly used herbicides and pesticides are 24D and Deltamethrin, respectively. Management factors such as application rate of fertilizer, crop rotation, variety of winter wheat, seeding method, planting density, setting of planting, machinery etc. are similar throughout the study area.

Description of the Study Area
The study area had an area of 53.20 km 2 and was located in the city of Ghorveh, western Iran ( Figure 1). The climate is semi-arid with a mean annual rainfall of 369.8 mm, mainly concentrated in the winter (October to March). The yearly mean temperature (10.8 °C) is relatively cool. The area is categorized by a Xeric soil moisture regime and a Mesic soil temperature regime. Elevation changes from 1833 to 2627 m and slopes range from 1 to 30% (Figure 1). The primary land use is cropland, where farm machinery operations with conventional tillage systems are commonly completed by 15th October. Winter wheat (Triticum aestivum L.) was planted on 15th October at a depth of 4 cm, 17 cm row spacing and a seed rate of 160 kg ha −1 . Burning crop residue is a common practice in the study area. Wheat and fallow is a common 2-year crop rotation. Triple superphosphate fertilizer was applied at a rate of 50 kg ha −1 , on wheat fields before sowing. Nitrogen fertilizer (urea) was applied at a rate of 100 kg ha −1 , on wheat fields before sowing, and at the stem elongation stage. Commonly used herbicides and pesticides are 24D and Deltamethrin, respectively. Management factors such as application rate of fertilizer, crop rotation, variety of winter wheat, seeding method, planting density, setting of planting, machinery etc. are similar throughout the study area. Based on the USDA soil taxonomy classification method [34] and the World Reference Base for Soil Resources (WRB) of the FAO [35], the majority of the soils in the study area were Entisols (Leptosols) and Inceptisols (Cambisols). The main physiographic units were piedmont plain, fan, plateau, and hill ( Figure 1). The soils of the piedmont plain and the majority of the plateau were Inceptisols and were deep with a surface horizon (Ap) and subsurface horizon (Bk and Bw). Conversely, the soils of the hill and the majority of fans, called Entisols, were thin with surface Based on the USDA soil taxonomy classification method [34] and the World Reference Base for Soil Resources (WRB) of the FAO [35], the majority of the soils in the study area were Entisols (Leptosols) and Inceptisols (Cambisols). The main physiographic units were piedmont plain, fan, plateau, and hill ( Figure 1). The soils of the piedmont plain and the majority of the plateau were Inceptisols and were deep with a surface horizon (Ap) and subsurface horizon (Bk and Bw). Conversely, the soils of the hill and the majority of fans, called Entisols, were thin with surface horizon (Ap) and lithic contact near to the soil surface. Each stage of this work was carried out and Figure 2 summarizes these procedures as a flowchart. The stages of this work were: i-Selection of 100 soil samples based on two soil orders (Entisols and Inceptisols); ii-Measuring 18 soil physical, chemical, and biological characteristics as a total data set. Actual yield of rainfed winter wheat was also measured at 100 sampling points using a 1 m 2 quadrat; iii-Removing the soil indicators that did not demonstrate a consistent correlation with rainfed winter wheat yields from the total data set; iv-Determining the minimum data set (MDS) using principal component analysis (PCA); v-Scoring of each soil indicator in the MDS in the range of 0-1 using two standard scoring functions (SSF) including "More is better" (MO) and "Less is better" (LE) to standardize soil indicators measured on different scales; vi-Weighting each soil indicator in the MDS using two weighting approaches including multi linear regression (MLR) and factor analysis (FA); vii-Calculating the two weighted additive soil quality indices (SQIw): SQIw-F (FA method) and SQIw-M (MLR method); viii-Validating each SQI (SQIw-F and SQIw-M) using the establishment of correlations between them and the actual yields of rainfed winter wheat; ix-Grading each SQI (SQIw-F and SQIw-M) into five grades (very high, high, moderate, low, and very low); x-Mapping each SQI (SQIw-F and SQIw-M) using digital soil mapping (DSM) approaches.
Agriculture 2020, 10, x FOR PEER REVIEW 4 of 18 horizon (Ap) and lithic contact near to the soil surface. Each stage of this work was carried out and Figure 2 summarizes these procedures as a flowchart. The stages of this work were: i-Selection of 100 soil samples based on two soil orders (Entisols and Inceptisols); ii-Measuring 18 soil physical, chemical, and biological characteristics as a total data set. Actual yield of rainfed winter wheat was also measured at 100 sampling points using a 1 m 2 quadrat; iii-Removing the soil indicators that did not demonstrate a consistent correlation with rainfed winter wheat yields from the total data set; iv-Determining the minimum data set (MDS) using principal component analysis (PCA); v-Scoring of each soil indicator in the MDS in the range of 0-1 using two standard scoring functions (SSF) including "More is better" (MO) and "Less is better" (LE) to standardize soil indicators measured on different scales; vi-Weighting each soil indicator in the MDS using two weighting approaches including multi linear regression (MLR) and factor analysis (FA); vii-Calculating the two weighted additive soil quality indices (SQIw): SQIw-F (FA method) and SQIw-M (MLR method); viii-Validating each SQI (SQIw-F and SQIw-M) using the establishment of correlations between them and the actual yields of rainfed winter wheat; ix-Grading each SQI (SQIw-F and SQIw-M) into five grades (very high, high, moderate, low, and very low); x-Mapping each SQI (SQIw-F and SQIw-M) using digital soil mapping (DSM) approaches.

Sampling and Soil Analysis
Soil and winter wheat samples were collected at 100 locations, from the center of 1 m 2 . Soil samples were collected soon after the winter wheat harvest. Figure 1d shows that at the soil order level, 35 and 65 composite samples were collected at a depth of 0-20 cm from the Entisol and Inceptisol orders, respectively. Each composite soil sample (0-20 cm depth) was a combination of five subsamples (one sample from each corner and one from the center of a 1 m 2 quadrat).

Measurement of Actual Yield
From each sample site, wheat samples and plants (planting density within 1 m 2 was uniform) were harvested by hand within a 1 m 2 quadrat (in August 2017). All the plant samples were oven-dried for 48 h at 75°C and then wheat grains were weighed to obtain their dry grain weights. Wheat yield (Equation (1)) for all of the 100 sampling sites were gathered from the farms.
where Y is yield (t ha −1 ), N 1 is ear number per m 2 , N 2 is grain number per ear (yield components in wheat), and N 3 is thousand grain weight.

Determining the Minimum Data Set (MDS)
The total data set comprised the following soil properties: soil texture, pH, OC, BD, PO, EC, MWD, CEC, CCE, SMR, ESP, TN, AK, AP, AFe, AMn, AZn, and ACu. These soil properties were selected based on previous literature [5,[15][16][17]22,24,26,27,33,[50][51][52][53][54]. Initially, the soil quality indicators that did not demonstrate a significant correlation with wheat yields were removed and were not included in the MDS [6]. Then, PCA was applied to select soil quality indicators for the MDS [9,13]. In the PCA method, those PCs with large Eigenvalues (≥1) were selected as soil quality indicators of the MDS. For each PC, variables received a weight representing their contributions to the PC. Only highly weighted variables were selected from each PC to define the minimum data set [13,31].

Indicator Scoring of Soil Quality Indicators
A linear scoring function for each soil quality indicator in the MDS was used to eliminate the effect of soil quality indicators measured on different scales [13,51]. Two standard scoring functions (SSF), "More is better" (Equation (2)) and "Less is better" (Equation (3)) were applied to allocate a score in the 0-1 range to each soil quality indicator. The "More is better" SSF was used for all of the MDS indicators except sand where the "Less is better" SSF was used. These standard scoring functions were selected according to previous literature (Table 1). For the "More is better" SSF, each soil quality indicator was divided by the maximum observed value such that maximum values received a score of 1. Conversely, for "Less is better", each soil quality indicator was divided by the minimum value such that the lowest value received a score of 1. Table 1 shows the maximum and minimum threshold values for all soil quality indicators.
where Y is the linear score in the range of 0 and 1, X is the variable value, L is the lowest value, and H is the highest value [55].

Weighting of Soil Quality Indicators
Soil quality indicator weights in the MDS were assigned by two weighting approaches (MLR and FA). In the MLR method, modeling was carried out by using the wheat yield and the MDS soil quality indicators as dependent and independent variables, respectively. The weight of each soil quality indicator was defined as the ratio of its standardized regression coefficient divided using the summation of the standardized regression coefficients of all MDS soil quality indicators. In the FA method, soil quality indicators were weighted based on communality. Communality of each soil quality indicator explained a certain percentage of the variation in the MDS data. For all soil quality indicators, weight values resulted from the ratio of communality of all soil quality indicators in the MDS to the summation of each soil quality indicator in the MDS approaches [56].

Calculating of Soil Quality Index
According to the method proposed by [58,59], the chosen soil indicators were used to calculate the Weighted Additive Soil Quality Index (SQIw) (Equation (4)): where Wi is the weighting factor for soil indicators and Ni is the indicator score. Two SQIs, the SQIw-F (FA method) and the SQIw-M (MLR method) were computed using the two weighting methods.

Validation of the Soil Quality Index
The measured wheat yield was used for validation of soil quality indices. Correlations between each SQIw (SQIw-F SQIw-M) and actual yield of rainfed winter wheat were analyzed.

Grading and Mapping of Soil Quality Index
The two SQIs, SQIw-F (FA method) and the SQIw-M (MLR method), were separated into five grades (very high, high, moderate, low, and very low). Then, the two SQIs maps were plotted from predictions of soil quality at unmeasured locations using digital soil mapping (DSM) approaches. DSM approaches are useful tools for mapping soil properties [60]. In DSM, quantitative models have been applied to link between soil observations and auxiliary variables to understand spatial variation in soil properties [60]. In this approach a set of auxiliary variables (e.g., terrain attributes) is applied to map different soil properties [17,25,61]. In the current research, an artificial neural network (ANN) model and a set of auxiliary variables (i.e., soil map, terrain attributes, and remotely sensed data) were used to predict and map the SQIs (SQIw-M and SQIw-F) in the R software environment. Validation of the DSM maps and the ANN model was performed using the coefficient of determination (R 2 ) and root mean square error (RMSE) and through leave-one-out (LOO) cross-validation.

Statistical Analysis
Prior to statistical analysis, the normality of all data was examined using the Kolmogorov-Smirnov test and visual examination of histograms. All statistical analyses including the two weighting methods (multi linear regression and factor analysis) and principal component analysis to determine minimum data set were performed with the statistical software SPSS 17.0 [62]. Significant differences between all of the measured soil properties were tested using an independent samples t-test. The statistical significance of the difference was p < 0.05. Table 2 shows the descriptive statistics of the soil properties measured and rainfed wheat yield at each sampling point. Almost all of the soils were basic (pH 7.60-8.92) with a mean pH of 8.09. EC values were low (0.1 and 0.30 dSm −1 ) with a mean value of 0.17 dSm −1 . The OC and TN contents were also low in the studied soils, with mean values of 0.15% (range = 0.06 to 0.52%) and 0.02% (range = 0.01 to 0.10%), respectively. Based on the general evaluation of the coefficient of variation (CV), the CVs of OC and TN were high (more than 35%), suggesting a high variability across the study area. This large variation in OC and TN can be mostly due to the variability in topography. CEC, ESP, MWD, BD, SMR, and CCE were in the range of 6.25-29.00 cmoL ckg −1 , 1.11-5.72%, 0.55-5.10 mm, 1.1-1.6 g cm −3 , 75.70-205 mg CO2 g day −1 , and 1.25-30.25%, respectively. The average sand, silt, and clay contents were 48.66%, 28.18%, and 23.15%, respectively. The mean AK, AP, ACu, AZn, AFe, and AMn contents were 128. 67, 13.14, 0.43, 0.28, 3.49, and 1.1 mg kg −1 , respectively, and were randomly dispersed across the study area (with a CV of more than 35%, except for AP, AK). The mean rainfed wheat yield was 1.16 ton ha −1 and ranged between 0.5 and 1.85 ton ha −1 . The CV of rainfed wheat yield was moderate (25.71%). Electrical conductivity (EC); soil microbial respiration (SMR); available potassium (AK); available phosphorus (AP); cation exchange capacity (CEC); soil acidity (pH); organic carbon (OC); bulk density (BD); exchangeable sodium percentage (ESP); mean weight diameter (MWD); total Nitrogen (TN); soil porosity (SP); calcium carbonate equivalent (CCE); available Fe (AFe); available Zn (AZn); available Cu (ACu); available Mn (AMn).

The Minimum Data Set (MDS)
Based on the Pearson's correlation coefficients among the rainfed wheat yield and all soil indicators (Table 3), soil indicators including silt, SP, BD, AZn, ACu, AFe, CCE, pH, and EC were deleted because they were not significantly correlated with rainfed wheat yield. The soil indicators retained for PCA comprised OC, sand, clay, MWD, CEC, SMR, ESP, TN, AK, AP, and AMn. The first six PCs had eigenvalues of more than 1 and explained 98% of the variance of the original data ( Table 4). The PCA results showed that SMR had the highest loadings (0.866) with respect to PC1 compared to other soil quality indicators. AK (−0.855), sand (0.882), CEC (−0.889), clay (−0.816), and AP (0.900) also had the highest loadings with respect to PC2, PC3, PC4, PC5, and PC6, respectively, compared to other soil quality indicators (Table 4). OC also had a high loading (−0.603) for PC4, but had a strong correlation (r > 0.60) with CEC and therefore it was not used as a property to represent PC4.

Weighting
MLR and FA were applied to allocate the soil indicator weights. In the MLR method, clay and AK had the highest weights (0.20 and 0.193, respectively), sand and SMR had the lowest weights (0.132 and 0.134, respectively), and AP and CEC had intermediate weights (0.163 and 0.177, respectively), (Table 5). In the FA method, according to communality analysis (Table 5), clay and CEC had maximum weights (0.217 and 0.196, respectively), sand and SMR had minimum weights (0.100 and 0.136, respectively), and AK and AP had intermediate weights (0.179 and 0.172, respectively).

SQI and Weighting Methods
Using the MLR weighting technique (SQIw-M), the values of SQIw-M for the studied soils varied from 0.16 to 0.73, with an average of 0.559 and a CV of 21.16%. Using the FA weighting technique (SQIw-F), the values of SQIw-F varied from 0.17 to 0.77, with a mean of 0.579 and a CV of 22.09%. Table 6 shows the grades of soil quality indices (SQIw-M and SQIw-F). Percentage of grades I (very high quality), II (high quality), III (moderate quality), IV (low quality), and V (very low quality) in the study area for the SQIw-F, were 1%, 26.03%, 48.03%, 18.27%, and 6.67% respectively, and for the SQIw-M were 0.25%, 23.55%, 50.06%, 11.04%, and 15.10%, respectively. For SQIw-F and SQIw-M, grade III was dominant, accounting for 48.03 and 50.06% of all locations, respectively, therefore, the majority of the study area has moderate soil quality (  Figure 3). Conversely, there were lower percentages of grade IV and V of SQIw-M for Inceptisols (3.05% and 5.79% for SQIw-M, and 3.85% and 2.10% for SQIw-F, respectively), compared to Entisols (7.99% and 9.31% for SQIw-M, and 14.42% and 4.57% for SQIw-F, respectively) (Table 6, Figure 3). Table 6. Classification and percentage of soil quality grades for two soil quality indices (SQIw-M and SQIw-F) in two soil orders.

SQI According to Soil Orders
The SQIs were calculated using two weighting methods, SQIw-M and SQIw-F. For Entisols, the SQIw-M ranged from 0.22 to 0.72 and SQIw-F ranged from 0.20 to 0.77, with an average of 0.54 and 0.56, respectively. For Inceptisols, the SQIw-M and SQIw-F values ranged between 0.16 to 0.73 and 0.16 to 0.76, averaging 0.57 and 0.59, respectively ( Table 6). The mean values for both SQIs (SQIw-F and SQIw-M) were lower in Entisols than in Inceptisols although these values were not significantly different (Table 7). For both SQIw-M and SQIw-F, the majority of grades I, II, and III were located in the piedmont plain and plateau geomorphology units. In contrast, the majority of grades IV and V of both SQIw-M and SQIw-F were located in the hill geomorphology unit (Figures 1 and 3). Table 7. Results of independent samples t-test for soil properties and SQIs.

SQI According to Soil Orders
The SQIs were calculated using two weighting methods, SQIw-M and SQIw-F. For Entisols, the SQIw-M ranged from 0.22 to 0.72 and SQIw-F ranged from 0.20 to 0.77, with an average of 0.54 and 0.56, respectively. For Inceptisols, the SQIw-M and SQIw-F values ranged between 0.16 to 0.73 and 0.16 to 0.76, averaging 0.57 and 0.59, respectively ( Table 6). The mean values for both SQIs (SQIw-F and SQIw-M) were lower in Entisols than in Inceptisols although these values were not significantly different (Table 7). For both SQIw-M and SQIw-F, the majority of grades I, II, and III were located in the piedmont plain and plateau geomorphology units. In contrast, the majority of grades IV and V of both SQIw-M and SQIw-F were located in the hill geomorphology unit (Figures 1 and 3).

Different Weighting Methods
For validation of the SQIs, the relations between both rainfed wheat yield and SQIs were examined. The rainfed winter wheat yield showed a significant correlation with SQIw-M and SQIw-F for the total data set. It also showed that the MLR weighting technique was slightly more accurate (r = 0.62) than the FA (r = 0.58) as a weighting method (Figure 4). For validation of the SQIs, the relations between both rainfed wheat yield and SQIs were examined. The rainfed winter wheat yield showed a significant correlation with SQIw-M and SQIw-F for the total data set. It also showed that the MLR weighting technique was slightly more accurate (r = 0.62) than the FA (r = 0.58) as a weighting method (Figure 4).

Different Soil Orders
There is a significant correlation with the rainfed winter wheat yield and SQIw-F and SQIw-M for both Entisols and Inceptisols. The MLR weighting technique in Entisols had slightly higher accuracy (r = 0.47) than the FA weighting technique (r = 0.43) ( Figure 5) and the MLR weighting method was also more accurate for Inceptisols (r = 0.67) than the FA weighting method (r = 0.63) ( Figure 5). The values for two SQIs (SQIw-M and SQIw-F) and rainfed wheat yield were lower for Entisols than for Inceptisols; however, these values were not significantly different (Table 6).

Different Soil Orders
There is a significant correlation with the rainfed winter wheat yield and SQIw-F and SQIw-M for both Entisols and Inceptisols. The MLR weighting technique in Entisols had slightly higher accuracy (r = 0.47) than the FA weighting technique (r = 0.43) ( Figure 5) and the MLR weighting method was also more accurate for Inceptisols (r = 0.67) than the FA weighting method (r = 0.63) ( Figure 5). The values for two SQIs (SQIw-M and SQIw-F) and rainfed wheat yield were lower for Entisols than for Inceptisols; however, these values were not significantly different (Table 6).

The Minimum Data Set (MDS)-
Eighteen physical, chemical and biological soil properties (as total data set) were measured. Then, an MDS was established to assess soil quality for rainfed wheat growth across soil orders at the regional scale. The MDS selection procedure has been commonly suggested by other researchers to assess soil quality [27,57,63]. The temporal and financial resources for soil sampling and analysis can be reduced by using the MDS method. However, in comparison to the total data set, some soil quality information is lost. The selected six soil properties (clay, sand, CEC, SMR, AP, and AK) used in SQIs formed with the MDS method are similar to the findings of the other studies [5,6,26,27,51,52,54].

Evaluation of Approaches to Weighting SQIs
The results of the current study showed that the MDS soil indicators (clay, sand, CEC, SMR, AP, and AK) and rainfed wheat yield had higher significant correlations with the calculated soil quality index using the MLR weighting method (SQIw-M) than with the soil quality index calculated using the FA weighting method (SQIw-F) ( Table 3, Figures 3 and 5). Therefore, the MLR weighting technique is more accurate than the FA weighting technique and should be used for future soil quality evaluations for rainfed wheat in western Iran. Cheng et al. [6] indicated that SQIs calculated using the MLR weighting method (SQIw-M) had a significantly higher correlation with orange yield (R 2 = 0.43) compared to the FA weighting method SQIw-F (R 2 = 0.35). Similar results were also reported in Northeast and Central China by Chen et al. [28] when studying a soybean production system.

The Minimum Data Set (MDS)
Eighteen physical, chemical and biological soil properties (as total data set) were measured. Then, an MDS was established to assess soil quality for rainfed wheat growth across soil orders at the regional scale. The MDS selection procedure has been commonly suggested by other researchers to assess soil quality [27,57,63]. The temporal and financial resources for soil sampling and analysis can be reduced by using the MDS method. However, in comparison to the total data set, some soil quality information is lost. The selected six soil properties (clay, sand, CEC, SMR, AP, and AK) used in SQIs formed with the MDS method are similar to the findings of the other studies [5,6,26,27,51,52,54].

Evaluation of Approaches to Weighting SQIs
The results of the current study showed that the MDS soil indicators (clay, sand, CEC, SMR, AP, and AK) and rainfed wheat yield had higher significant correlations with the calculated soil quality index using the MLR weighting method (SQIw-M) than with the soil quality index calculated using the FA weighting method (SQIw-F) ( Table 3, Figures 3 and 5). Therefore, the MLR weighting technique is more accurate than the FA weighting technique and should be used for future soil quality evaluations for rainfed wheat in western Iran. Cheng et al. [6] indicated that SQIs calculated using the MLR weighting method (SQIw-M) had a significantly higher correlation with orange yield (R 2 = 0.43) compared to the FA weighting method SQIw-F (R 2 = 0.35). Similar results were also reported in Northeast and Central China by Chen et al. [28] when studying a soybean production system.

Soil Quality Indicators and Suggestions for Soil Quality Management
In the study region, the mean AK contents in both Entisols and Inceptisols (122.79 and 131.48 mg kg −1 , respectively) ( Table 7), were smaller than the critical amount of soil AK content (200 mg kg −1 ) needed for rainfed wheat based on the soil nutrient classification criteria of the Soil and Water Research Institute of Iran [64][65][66]. The potassium status of different soil orders (Inceptisols and Entisols) depends on the parent material, texture, mineralogy, and management practices [67]. Traditionally in the study region, application of K fertilizer is not a common practice of farmers. Consequently, it is recommended that K fertilizer is used in the studied area, particularly in fields where AK contents are less than the critical limit. It should be noted that based on the critical amount of soil AK content (200 mg kg −1 ), response probability to potash fertilizer consumption is lower than 50%. In contrast, the mean AP contents in Entisols and Inceptisols were 12.32 and 13.59 mg kg −1 , respectively, which were more than the moderate critical amount of soil AP content (7 mg kg −1 ) needed for rainfed wheat based on the categorization criteria of soil nutrients of the Soil and Water Research Institute of Iran [64][65][66]. In this study, AP content was more than the critical amount due to annual P fertilizer being commonly applied in rainfed wheat fields.
In this study, the soil textures were sandy clay loam, loam, sandy loam, and clay loam; and sand content ranged between 35.20 to 75.68%. The variability of sand percentage affects the capability of soil tillage, subsequently affecting wheat growth and yield. Most of the studied rainfed wheat was found on soils underlain by granite and granodiorite, which have coarse textures; therefore, soils derived from these rocks also have a coarse texture. Soil texture is a permanent soil property, which controls the retention and uptake of soil water, nutrients, and oxygen. Coarse texture means a higher percentage of sand as compared to silt or clay, therefore such soils have greater aeration and poor water retention which is not suitable for rainfed wheat sites.
In general, as SOC serves as food in the soil, instead of pure fertilizers an application of organic fertilizer can lead to an increase in nutrients, accumulation of SOC, and can amend soil texture.

Links between Soil Orders and SQIs
The mean SQI values showed higher quality for Inceptisols than for Entisols; however, this difference was not statistically significant ( Table 7). The mean values of rainfed winter wheat yield were also greater for Inceptisols than Entisols due to slope and erodibility differences between these soil types. Nevertheless, this difference between the soil types was not statistically significant (Table 7). Entisols were characterized as younger, less developed and shallower soils than Inceptisols. Entisols were mostly located on the hill geomorphology units, which are more susceptible to soil erosion and soil quality degradation due to steeper slopes. In comparison, Inceptisols orders were located on the lower slopes (piedmont plain and plateau geomorphology units), which could influence many soil processes that are important for soil formation (Figure 1). In several studies, it has been demonstrated that slope has a marked influence on soil quality [25,68,69]. This result demonstrates the impact of soil type on the soil quality of the study region. Many studies have noted that soil type has a significant impact on soil quality [7,70,71]. Biswas et al. [5] showed that a SQI in cropping systems (rice-rice) in South Asia was greatest for Inceptisols (0.66 to 0.89), then Entisols (0.23 to 0.76), and was lowest for Alfisols (0.37 to 0.60) as the Alfisols were highly weathered, while Inceptisols/Entisols were relatively younger in nature. In contrast, Santos-Frances et al. [15] showed that Alfisols, Inceptisols, and Entisols had a high, medium, and low soil quality index, respectively, in a semi-arid ecosystem of Spain. Alfisols-covered flat surfaces (plateaus), were not very stony, and were thick.
Results also showed that the areas of the soil order Entisol had lower proportions of high soil quality grades (Grades I and II) ( Table 6). Conversely, grades with low soil quality (Grades IV and V) were more abundant in this soil order compared to those of Inceptisols. The soil quality of Entisols was significantly restricted by high slopes and erodibility; it therefore had low overall quality. In the study area, the lower SQI value for Entisols was due to its intrinsic properties, high slopes, and erosion potential. Moreover, this soil type generally has lower development and did not have sufficient quantities of macro-and micro-nutrients (Table 7). This characteristic was a main reason for these soils mostly being classified as low-quality soil [15]. Consequently, for Entisols management practices including application of organic matter, chemical fertilizers and soil amendments could increase wheat yield [72,73].
The main factors affecting soil quality include climate, topography (geomorphological unit), soil type, plant species, and agricultural management practices like irrigation, fertilization, and tillage. Apart from soil type and topography, the factors influencing soil quality are similar amongst all the sampling sites. Therefore, the variation of soil types and topography could explain the majority of soil quality variability in this area and soil type must be considered for maintaining and increasing soil quality and yield sustainability in the wheat cropping systems for semi-arid Iran.

Conclusions
This study indicated that the MLR weighting technique is superior compared to the FA weighting technique with regard to its accuracy. Accordingly, it is suggested that MLR weighting be used for future soil quality evaluations in this region. The soil quality indicators retained in the MDS included AK, clay, CEC, AP, SMR, and sand. The study region's soil quality is mostly poor and moderate; therefore, more attempts should be made to apply organic inputs like farm yard manure, compost or green manure and potassium (K) fertilizer in rainfed wheat fields. Generally, the soils of the Inceptisol order exhibited greater soil quality and rainfed wheat yield than soils of the Entisol order. For the Entisol order, the low SQIs especially depend on chemical, physical, and biological soil properties. Consequently, the improvement of soil properties with sustainable management and agricultural practices is necessary to enhance the SQI in the soils of the study area to increase wheat yield. Therefore, soil type must be taken into account for assessing soil quality and implementing appropriate practices. However, it is recommended for future studies to consider other driving factors such as soil quality, cost, work efficiency, occupational safety, and environmental load. Furthermore, it is suggested that in future studies soil moisture should be measured as one of the important soil quality indicators due to its effect on the ability of crops to effectively absorb nutrients from fertilizers.