Land suitability assessment for major crops by using GIS-based multi-criteria approach in Andit Tid watershed, Ethiopia

Land evaluation is a foundation for sustainable land resource planning and management. In Ethiopia, national economic development is highly dependent on agriculture. And, the productivity of the land is low as a result of land degradation due to mismanagement of the land and intensive cultivation; this situation is worse in Ethiopian Highlands like Andit Tid watershed. Therefore, matching crop requirement with resource available through land suitability analysis has become an urgent need to sustain the productivity of agricultural land in the study area. This study aimed to address land suitability analysis by using geographical information system-based multi-criteria approach in Andit Tid watershed. The analysis of land suitability for major crops (wheat and barley) revealed that, though, there are slight variation in suitable classes for each crop, most part of the watershed was moderately suitable for both crops with 266.55 ha (77.42%) and only very small part of the land 2.43 ha (0.51%) was highly suitable for both crops. On the other hand, 9.3 ha (2%) and 3.6 ha (0.77%) of the land was marginally suitable and not suitable for both wheat and barley crops, respectively. The main limiting factors for land suitability were soil depth, texture, temperature, slope and erosion hazard. The Hamere Yohannes ABOUT THE AUTHOR Hamere Yohannes is a PhD candidate in Center for Environmental Science at Addis Ababa University and lecturer at department of Natural Resource Management in Debre Berehan University. She teaches land degradation and rehabilitation, introduction to environmental science, climatology, environment and development and land-use planning courses and engaged on research and community services. She published more than eight research and review articles in the internationally peer reviewed journals. Teshome Soromessa (PhD) is associate professor of Environmental Science at the College of Natural Science in Addis Ababa University and director for Center for Environmental Science. He specialized in Ecology and Diversity. He teaches vegetation ecology and restoration courses. He actively engaged on society problem solving researches. Dr. Teshome published more than 125 papers in reputable internationally peerreviewed journals and advised more than 90 PhD and MSc. students. PUBLIC INTEREST STATEMENT In Ethiopia, Agriculture is the basis for economy. Wheat, barley, maize, sorghum, oats and millet are major cereals which are essential for this economy. Currently, most of these crops cultivated in less suitable and degraded areas due to mismanagement of the land, shortage of fertile land and pressing demand for food. Consequently, low productivity becomes problematic in most part of Ethiopian Highland areas. To alleviate this problem, land evaluation which identifies the suitable and not suitable land for crop production is significant to plan management. This study highlights the role of GIS-based multi-criteria approach to identify areas which are highly suitable, moderately suitable and not suitable for barley and wheat crops in the study area. The result showed that only small part of the watershed was highly suitable for those crops and most part of the area was moderately suitable. Practically, land has been used against its suitability classes in the study area. Yohannes & Soromessa, Cogent Food & Agriculture (2018), 4: 1470481 https://doi.org/10.1080/23311932.2018.1470481 © 2018 The Author(s). This open access article is distributed under a Creative Commons Attribution (CC-BY) 4.0 license. Received: 27 February 2018 Accepted: 25 April 2018 First Published: 08 May 2018 *Corresponding author: Hamere Yohannes, Center for Environmental Science, College of Natural and Computational Sciences, Addis Ababa University E-mail: hamerenew@gmail.com; hamerey@gmail.com Reviewing editor: Manuel Tejada Moral, University of Seville, Spain Additional information is available at the end of the article

Land suitability assessment for major crops by using GIS-based multi-criteria approach in Andit Tid watershed, Ethiopia Hamere Yohannes 1,2 * and Teshome Soromessa 1 Abstract: Land evaluation is a foundation for sustainable land resource planning and management. In Ethiopia, national economic development is highly dependent on agriculture. And, the productivity of the land is low as a result of land degradation due to mismanagement of the land and intensive cultivation; this situation is worse in Ethiopian Highlands like Andit Tid watershed. Therefore, matching crop requirement with resource available through land suitability analysis has become an urgent need to sustain the productivity of agricultural land in the study area. This study aimed to address land suitability analysis by using geographical information system-based multi-criteria approach in Andit Tid watershed. The analysis of land suitability for major crops (wheat and barley) revealed that, though, there are slight variation in suitable classes for each crop, most part of the watershed was moderately suitable for both crops with 266.55 ha (77.42%) and only very small part of the land 2.43 ha (0.51%) was highly suitable for both crops. On the other hand, 9.3 ha (2%) and 3.6 ha (0.77%) of the land was marginally suitable and not suitable for both wheat and barley crops, respectively. The main limiting factors for land suitability were soil depth, texture, temperature, slope and erosion hazard. The Hamere Yohannes ABOUT THE AUTHOR Hamere Yohannes is a PhD candidate in Center for Environmental Science at Addis Ababa University and lecturer at department of Natural Resource Management in Debre Berehan University. She teaches land degradation and rehabilitation, introduction to environmental science, climatology, environment and development and land-use planning courses and engaged on research and community services. She published more than eight research and review articles in the internationally peer reviewed journals.
Teshome Soromessa (PhD) is associate professor of Environmental Science at the College of Natural Science in Addis Ababa University and director for Center for Environmental Science. He specialized in Ecology and Diversity. He teaches vegetation ecology and restoration courses. He actively engaged on society problem solving researches. Dr. Teshome published more than 125 papers in reputable internationally peerreviewed journals and advised more than 90 PhD and MSc. students.

PUBLIC INTEREST STATEMENT
In Ethiopia, Agriculture is the basis for economy. Wheat, barley, maize, sorghum, oats and millet are major cereals which are essential for this economy. Currently, most of these crops cultivated in less suitable and degraded areas due to mismanagement of the land, shortage of fertile land and pressing demand for food.
Consequently, low productivity becomes problematic in most part of Ethiopian Highland areas. To alleviate this problem, land evaluation which identifies the suitable and not suitable land for crop production is significant to plan management. This study highlights the role of GIS-based multi-criteria approach to identify areas which are highly suitable, moderately suitable and not suitable for barley and wheat crops in the study area. The result showed that only small part of the watershed was highly suitable for those crops and most part of the area was moderately suitable. Practically, land has been used against its suitability classes in the study area. spatial distribution showed that most of the cultivation currently practised in marginally suitable land in the lower part of the watershed. Most of the land parcel being used against its suitability potential in study area. Thus, land-use pattern needs to be modified based on its suitability potential. Subjects: Research article; Agriculture; Agricultural Development; Crop Science; Environmental Sciences; Agriculture and Food Keywords: Andit Tid watershed; Barely; land suitability; multi-criteria decision analysis; Wheat

Background
Worldwide, the amount of cultivated land is limited, and most land has been degraded as irreversibility and become unsuitable for agricultural production (Verheye, 2008). Land evaluation is a foundation for sustainable land resource planning and management since it help us to know whether the resources are degraded or enhanced in quality (Dumanski et al., 2002;Mohana, Mariappan, & Manoharan, 2009). FAO (1976) defines land suitability as "the fitness of a given parcel of land for specific uses". Consequently, land suitability assessment is land evaluation which usually conducted to determine specific land use for a particular location and identify limiting factors for a particular crop production (AbdelRahman, Natarajan, & Hegde, 2016;Mu, 2006). The assessment of land suitability depends on land capability as well as other factors such as land quality, proximity to different accesses, land ownership, customers demand and economic values (Counsel, 1999).
Ethiopia's economic development is highly dependent on agricultural production since agriculture constitutes 46.6% of the national gross domestic product (CSA, 2008). Wheat, barley, maize, sorghum, oats, millet are major cereals which are essential for Ethiopia's agriculture and food economy (Taffesse, Dorosh, & Asrat, 2012). Most of agricultural practice has been done in Ethiopian highlands (Teshome, 2014). In this area, the competition for land is high due to demand for food and space as a result of population pressure. This results in degradation of resource in Ethiopian highlands. Marginal land has been used for cultivation and low productivity remains the major constraint for agriculture sector where yields are less than 1 t ha −1 (Engda, 2009;Pender & Gebremedhin, 2007;Verheye, 2008). Therefore, land suitability analysis becomes a critical step for land management planning and implementation in Ethiopia.
Land suitability evaluation according to FAO standards has been applied in many parts of the world, particularly in the developing countries (FAO, 2007b). In Ethiopia, a number of investigations were conducted to evaluate land suitability for crops and irrigation (Agidew, 2015;Alemu et al., 2013;Ayalew, 2014;Girma, Moges, & Quraishi, 2015;Hailu, Kibret, & Gebrekidan, 2015;Liambila & Kibret, 2016;Mekonen & Kebede, 2011;Melak, 2007;Mengesha, 2014;Motuma, Suryabhagavan, & Balakrishnan, 2016;Nahusenay & Kibebew, 2015;Rabia, 2012;Selassie, Ayalew, Elias, & Getahun, 2014;Teka & Haftu, 2012;Worqlul et al., 2017;Yitbarek, Kibret, Gebrekidan, & Beyene, 2013). These studies are inadequate to made decision on proper utilization of resources in Andit Tid watershed, because locally adapted evidence from the combination of topographic, soil, climatic and management factors is needed for proper management decision where the productivity of the land in the study area is limited; this situation also aggravated through short fallow period as a result of uncontrolled grazing (Engda, 2009). Consequently, agricultural production in the area has become unsustainable and marginal lands are being converted into cultivated lands. Moreover, according to Rabia and Terribile (2013), land suitability is a multi-criteria problem which involve with the number of parameters; therefore, the evaluation needs multi-criteria approach. The availability of GIS and multi-criteria decision analysis (MCDM) methods allow combining knowledge derived from different sources to support land-use planning and management (Malczewski, 2000). The most commonly used MCDM are analytic hierarchy process (AHP), weighted linear combination, ordered weighted averaging, ELECTRE, PROMETHEE, VIKOR and multiple-objective land allocation (Rikalovic, Cosic, & Lazarevic, 2014). Saaty (1996) proposed analytical network process (ANP) to overcome the problem of interdependence and feedback effects observed in AHP method. ANP method has become a new popular planning and decision-making approach that allows to determine the relative weights of the evaluative criteria (Alemu, 2011;Saaty, 2016). It enables to create a network unlike other MCDM methods. It is hardly addressed land evaluation by using by integrating GIS and ANP in Ethiopia in general; only a recent study by Worqlul et al. (2017) assessed land suitability for irrigation for Ethiopia. Therefore, this study aimed at land suitability evaluation by using GIS-based multi-criteria decision tool.

The setting
Andit Tid watershed situated 180 km northeast of Addis Ababa and approximately 60 km east of Debre Berehan adjacent to the Debre Berhan-Mekelle highway at 39°43′E and 9°48′N geographic location ( Figure 1). It is one of the six Soil Conservation Research Program watershed in Ethiopia since 1982 (Alemu et al., 2013). Water Land Resource Centre (WLRC) is responsible for gathering long-term data such as climatic, soil, erosion and sedimentation etc. in the study area. The catchment drains westward in to the Blue Nile basin, and to the east of the catchment, an escarpment separates the Shewa Plateau from the Awash Plains in the Rift Valley. The size of the watershed topographically estimated to be 477.5 ha (Engda, 2009).
The watershed characterized generally by humid climate which have mean annual temperature that is 12.54°C and average rainfall of 1,489 mm per year. The rainfall pattern of the area is bimodal, the smaller (belg) rainy season is usually from March to May and the main (kiremt) rainy season is from July to October (Figure 2).
The geological formation of the watershed belongs to volcanic rocks (rhyolites, trachites, tuffs and basalts) which is the result of tertiary volcanic eruption as similar to most central highland of Ethiopia (Hurni, 2000). Andosols, Cambisols, Fluvisols, Leptosols and Rogosols are different soil types available in the watershed (Figure 3). Among these soils, Cambisols (37%), Andosols (32%) and Rogosols (24%) cover the large portion of the study area. Andosols dominate in the wurch zone while Regosols dominate in the dega zone ( Figure 3).
Field observation in the study catchment revealed that smallholder cereal plough cultivation is the dominant cultivation system in the area and most of the cultivated fields are usually planted with barely. According to Engda (2009), wheat, linseed, peas and beans mix with barely are the common crops growing system in the area. Growing period of most crops in study area is about 175 days and the productivity of the land is limited.

Data sources and general methodology
Data were collected from primary field assessment and observation, questioner and discussion and secondary sources. Eleven expert's opinions from Debre Berehan Research Centre and field  assistance were used to identify the influential factors for suitability of the land for major crops through systematic questioners. The experts were selected based on their qualification and experience on natural resource management and agronomy.
A land evaluation procedure was set based on Food and Agriculture Organization (FAO) 1976 guideline in terms of land suitability. Topographic, soil, climate and management factors were used to evaluate the performance of land in Andit Tid watershed. These were selected based on agronomic knowledge of local experts and reviews of existing literatures such as FAO framework for land evaluation (FAO, 1976(FAO, , 1983(FAO, , 2007b. Eighteen criteria were used for land suitability classification for major crops namely barley and wheat in the watershed. These criteria were soil depth, soil texture, soil drainage (permeability), soil erosion, slope, aspect, altitude fertility and soil chemical characteristics (pH, CEC, electrical conductivity [EC], OM, available phosphorus, total nitrogen and calcium carbonate), climatic parameters (temperature and rainfall) and accessibly (distance from road and water point). The summary of methods used presented at Figure 4.

Current land-use mapping
According to FAO (1976), for land evaluation, delineating the current land-use boundary is the first step. Land-use type is a kind of land use described in terms of its products and management practices (FAO, 1996). The types of use considered are limited to those which appear to be relevant under the general physical, economic and social conditions prevailing in an area. These kinds of land use serve as the subject of land evaluation. They may consist of major kinds of land use such as rainfed agriculture, irrigated agriculture, grassland, forestry or recreation (FAO, 1976).
The current land use/cover was determined from Landsat satellite image downloaded from USGS website (https://earthexplorer.usgs.gov/). Landsat 8 OLI-TIRS of 30 m resolution acquired on 23 September 2017. The image was preprocessed and classified by using supervised classification method with the support of unsupervised classification in ERDAS Imagine 2014 software. In preprocessing stage, radiometric and geometric correction (AbdelRahman et al., 2016;Demissie, Yeshitila, Mengistie, & Thomas, 2017;Yohannes, Mohammed, & Elias, 2017) were performed to correct and enhance the image spectral quality. During image preprocessing, image was also layer stacked and create a subset later image classification was performed by collecting signatures from each class. And, about 125 ground controlling points were used. The study area was classified into five major classes (Table 1).
Finally, the classified image was verified with collected GPS points in the field and Google Earth imagery for some inaccessible control points and about a total of 250 random points were generated. During accuracy assessment, user accuracy, producer accuracy, overall accuracy and kappa coefficient were computed from error matrix as suggested by (Congalton, 2001) (Annex 1). For further analysis, the classified raster image was exported to ArcGIS software version 10.3.

Determination of land-use management conditions
Depending on the nature of the topographic, the combination of free survey and transect walk was used to observe the whole watershed. Observation and recording were done at each land-use type. The management condition depends on different factors such as land goods and service delivery, labour and capital intensity, infrastructure and income level (FAO, 1996). In this study, in addition of the status of management practice recording, ease of management data such as distance from road and river were prepared in GIS environment. Bare land Areas that have little or no vegetation cover mainly include degraded land, gorges and exposed big rocks

Land characteristics and quality assessment
Depending on the objectives of the evaluation, relevant important data sets may include soils, landform, climate, vegetation and surface and/or groundwater reserves which need to be prepared (FAO, 2007b). In this study, soil, climate and topographic data sets were included.

Soil sampling and analysis
To determine the amount and spatial distribution of soil structure and nutrient availability, field soil sampling was used. Stratified random sampling was applied based on soil survey manual (FAO, 2007a). The stratification was based on altitudinal variation. Five stratums were delineated and all the auger observations were georeferenced with the help of hand held GPS ( Figure 5). The core sample was taken in zig-zag manner such that all points in the area have an equal chance of being selected as a sampling point (FAO, 2000). A total of 57 soil samples were collected, by using Edelman auger at the surface layer (0-30 cm), and bulked into nine composite soil samples for surface soil characterization and crop suitability evaluation purpose. At laboratory called Horticoop Ethiopia (Horticultural) PLC Soil and water analysis Laboratory Bishoftu, Ethiopia, the soil samples were air dried, crushed with mortar and pestle, mixed well and passed through a 2-mm sieve for analysis of chemical characteristics and major nutrient contents (pH, organic matter, electric conductivity, cation exchange capacity, available phosphors, calcium carbonate and total nitrogen) of the soil identified by using scientific procedures and standard methods (Table 2).
Moreover, soil depth, drainage and texture data were obtained from WLRC. The maps of these parameters were prepared in GIS environment.

Soil loss estimation
Soil erosion was estimated by using Revised Universal Soil Loss Equation for Complex Terrain (RUSLE3D) model.
The RUSLE3D model calculates potential average soil loss (A) as follows (1) Figure 5. Soil sample points at each stratum.
A is the annual soil loss by water (t ha −1 year −1 ), R is the rainfall-runoff erosivity factor (MJ mm ha −1 year −1 ), K is the soil erodibility factor (t ha −1 per unit R), LS is the slope length and slope steepness factor [-], C is the cover and management factor [-], P is the conservation practices factor [-].

Topographic, climatic data and accessibility data
The 30-m spatial resolution digital elevation model ASTER GDEM was obtained from Aster Global Digital Elevation Map (http://gdex.cr.usgs.gov/gdex/) and used to generate slope, altitude and aspect by using Spatial Analyst Tool in ArcGIS environment. Rainfall and temperature data were obtained from WLRC and imported to GIS environment for map preparation. River and road map obtained from Ministry of Water and Energy.

Matching crops with requirements
Matching crops with their specified requirement which is suitability evaluation, was performed by using the maximum limitation method (FAO, 1983;Sys, Van Ranst, & Debaveye, 1991). This method takes the most limiting factors to define the land suitability classes. Selected parameters were placed in suitability classes by matching with the requirements of selected crops. In this case, barley and wheat cereal crops were selected for crop suitability evaluation (Annex 2). The selections of crop varieties were made based on their dominantly grown, their importance and low productivity as explained by experts and farmers in the studied area. The crop requirements are expressed by defining optimal, marginal and unsuitable conditions for each land attribute that directly or indirectly influences plant growth, performance and biomass production (Sarkar, Ghosh, & Banik, 2014).

Suitability analysis
For land suitability analysis, GIS and ANP software were used in combination, since suitable map was produced from the aggregated value of standardized rated and weighted criteria. Criteria standardization and rating were performed mainly by using GIS and criteria weight was processed using ANP model from the collected expert's opinion.

Criteria standardization and rating
Data are measured on different units as well as on different scales of measurements. Thus, it is necessary to standardize the factors before combination and assure they are transformed, if necessary, so that all factor maps are positively correlated with suitability (Voodg, 1983). Linear scale transformation is the most frequently used GIS-based approach for criteria standardization (Malczewski, 2004). Thus, criteria have been standardized for this study, by using reclassify spatial analyst tool, to make sure that each criterion had an equivalent measurement basis. Simultaneously during reclassification, factor ratings were also assigned for suitability analysis (Table 4). Factor ratings are sets of values that indicate how well each factor is satisfied by particular conditions of the corresponding land quality.

Criteria weight estimation
Weight assignment from pairwise comparison for the criteria was based on expert opinion and processed in superDecisions software by designing ANP model. The pairwise comparison technique developed (Saaty, 1980) in the context of a decision-making process is a ratio (reciprocal matrix) where each factor is compared with the other criteria, relative to its importance on a scale from 1 to 9 (Saaty, 2004). The first step in ANP development was to construct a network; here, we group the criteria into four clusters and each cluster has criteria then the interdependence of the criteria designed. Experts choice obtained from the designed questionnaires to obtain their judgement about the importance of each criterion as suggested by Saaty's (2004), scale (1-9), for instance 1 is equally important and 9 means extremely important. Consequently, 153 possible paired comparisons were made. The questionnaire example of the paired comparison presented in Annex 3. Then, formation of unweighted supermatrix and weighted supermatrix were performed. Finally, as shown in Table 3, priorities weight was calculated from pairwise comparison matrix undertaken and eigenvector values based on the formula below.
where w i is the sum of row for pairwise comparison and n is the size of matrix.
The consistency ratio (CR) needs to be computed to check the consistency of comparisons. The sum of criteria weights should be equal to 1. CR ranges from 0 to 1. A CR close to 1 indicates the probability that the matrix's rating was randomly generated. A CR of 0.1 or less is a reasonable level of consistency (Malczewski, 2000;Saaty, 1980). In this case, the CR of the matrix of paired comparisons between 18 factors for land suitability was 0.048, which indicates logical judgement.
Consistency index (CI) computed from the formula below: λ max is the largest eigenvalue of the pairwise comparison matrix and n is the number of classes.
Then, CR given by the following formula (Saaty, 1980): where RI is ratio index/average value of CI for random matrices using Saaty scale.

Aggregation of the weight and standardized rated criterion map
Weighted overlay method was used to aggregate standardized rated criteria and weighted criteria to map the suitable land based on equation below. These maps were reclassified based on parametric model of land index to generate FAO land classes which converts suitability values into classes to produce the final map (Table 4).
where S is the suitability, W i is the weight of factor i, and X i is the criterion score of factor i.

Land suitability classification for major crops
The spatial land suitability for wheat and barley crops was obtained after the analysis of selected evaluation criteria such as soil physical characteristics (texture, drainage and soil depth), erosion hazard, soil chemical characteristics and fertility status (pH, EC, CEC, Carbonate, OM, TN and AP), topographic factors (slope, altitude and aspect), climatic factors (temperature and rainfall) and ease of management (river/water and road access).

Topographic factors
The topographic feature of the study area stated by attitude, slope and slope direction as presented in Table 5, the elevation of the area ranges from 2,915 to 3,463 m a.s.l distributed in dega (2,915-3,202 m) with 30.32% and 69.70% of the area was wurch (above 3,203 m) altitude zone. The watershed generally characterized by moderately sloppy and only small part (10.51%) of the area was gentle slope. Most of the land slope faced to south and north direction (35.37%) and (33.64%), respectively.

Soil physical characteristics and soil loss
As shown in Table 6, the soil depth of the study area ranges from very shallow <30 cm to very deep (100-150 cm). Most of the soil was moderately deep (49.90%). Sandy loam (361.68 ha) and loamy sand (87.35 ha) were dominant soil textures. Well-drained soil dominated the study area (73.27%). Excessively drained soil covers 25.81% and the rest less than 1% was moderately drained soil. The soil loss ranges from 10 to 864 t ha −1 year −1 . The loss of the soil is generally very high and exceeds very much as compared to tolerable rate (1-12 t ha −1 year −1 ) by FAO (1996).

Soil chemical characteristics and fertility status
The soil chemical properties and fertility status of the soil of Andit Tid watershed are presented in Table 7. The spatial distribution of pH, EC, CEC, calcium carbonate, available phosphorus, organic matter and total nitrogen also presented.

pH and EC
The pH values varied from 5.52 to 6.19; this indicates that the soil of the studies area is moderately acidic. The salinity of the soil was indicated by EC of the soil and the values range between 0.03 and 0.05. The spatial distribution of pH and EC presented in Figure 6. Most of the soil pH ranges  between 5.8 and 5.92 in the middle of the watershed. Similarly, most of the soil EC value was 0.04 mS cm −1 , which distributed in all watershed. If the soil has EC < 2, it considered as non-saline soil and the salinity effect is negligible for most crops (Hazelton & Murphy, 2007).

Cation exchange capacity and calcium carbonate
The CEC of the soil varied between 33.09 and 57.03 meq/100 gm soil. The smallest CEC value found in the lower altitude dega zone and the higher values found at the top of the watershed at wurch zone. According to FAO standard, CEC values ranged from high (25-40) to very high (>40) indicated good OM in the soil. The calcium carbonate content of the soil ranges from 2.52% to 4.04% and most of the soil falls under the range of 2.78-3.02% as indicated in Figure 7 below. Calcium carbonate less than 7 generally favourable for agriculture (Hazelton & Murphy, 2007).

Available phosphorus, total nitrogen and organic matter
As shown in Table 8 total nitrogen content of the soil in Andit Tid watershed varied from 0.1% to 0.58%, this value ranges from medium to high according to Baruah and Barthakur (1997), soil chemical value interpretation. The spatial distribution showed in the upper part of the watershed; highest percentage of total nitrogen was available. On the other hand, in lower part of the watershed, there is limited per cent of total nitrogen as shown in Figure 8. Soil available phosphorus in the watershed ranges from 6.61 to 13.88 mg kg −1 ; therefore, most of the soil available  phosphorus was under medium range, since <5 (very low), 5-10 (low), 10-25 (medium) and 25-50 (high) based on (Barber, 1984) soil interpretation guideline.
Organic matter content of the soil varied from 1.77% to 14.24% (Table 8). If the soil has organic matter rating >4, it is considered as high rate as reported by Tan (1996). The highest organic matter was available at the top of the watershed where disturbance is less and dominated by   grass land. On the other hand, the lowest value distributed mostly in the lower part of the watershed where intensive cultivation was dominated (Figure 9).
In general, for the purpose of agriculture, most of the soil chemical characteristics and fertility indicator parameters showed good rating values as compared to scientific standards (Table 3).

Climate and accessibility factors
Mean monthly rainfall for two growing periods of the watershed ranges from 76.32 to 186.31 mm and annual temperature varied from 11.32 to 12.80°C. The spatial distribution of these climatic variables presented in Figure 10.
Land suitability also depends on the accessibility factor as reported by FAO (2007b), for the purpose of ease of management. In this case, by using the available river and road, the river distance varied from 0 to 776 m and the distance of the major road ranges from 0 to 3,152 m ( Figure 11). 12.9. Land suitable for wheat Figure 12 shows the map of suitability class of the land for wheat crop cultivation in the studied area. The land comprises S1 (highly suitable), S2 (moderately suitable), S3 (marginally suitable) and N2 (permanently not suitable) as described in Table 8.
The largest part of the watershed occupied by S2 class with 428.60 ha (90.56%), which indicated that most of the area was moderately suitable for wheat production. Only 4.39 ha (0.93%) of the area is highly suitable for wheat production (Table 9). Alemu et al. (2013) reported that 46% of the area in Anjeni watershed highly suitable for wheat crop. On the other hand, only 3% of the land estimated as highly suitable for wheat cultivation in Awulaelo District in Tigray (Rabia, 2012). In other study by Motuma et al. (2016), 20.5% of the area were classified as moderately suitable and 79.5% marginally Figure 11. map of distance from river and road. suitable for wheat in South Wello. In this case, 36.68 ha (7.75%) of the land estimated as marginally suitable and less than 1% was not suitable for wheat production. Girma et al. (2015) reported that 61% of the land was not suitable for wheat production and 6% moderately suitable and 33% marginally suitable for wheat crop in Jello watershed, Eastern Ethiopian highland.

Land suitable for barley
Suitable land for barley in Andit Tid watershed comprises S1, S2, S3 and N2 as shown in Figure 13. Most of the suitable lands found in upper part of the watershed. This is due to deep and fertile Andosols. This result consistent with Yohannes (1989) reported that good crop yield was found at the dega part of the watershed due to Andosols and longer fallow period.
Highly suitable class occupied only 3.33 ha (0.70%). Most part of the watershed was moderately suitable 430.98 ha (91.07%) for barley crop production. The rest area occupied by marginally suitable class 35.35 ha (7.47%) and not suitable rock (0.76%) (Table 10). A study by Rabia (2012)   reported that only 1.3.% of the total land occupied by highly suitable class and more than 80.5% of the land was not suitable for barely crop. Hailu et al. (2015) estimated that 9.93% of the study area is marginally suitable, 51.39% of the study area was permanently unsuitable, 38.68% was currently unsuitable for both wheat and barley production in Kabe watershed, South Wello zone. Similarly, 65.13%, 23.62%, and 11.25% of the land is moderately, marginally and not suitable, respectively, for wheat and barley crops (Nahusenay & Kibebew, 2015).
The main limitation factors for wheat and barley crops identified in the study area were soil depth, texture, temperature and slope and erosion hazard (Table 11). Teka and Haftu (2012) reported that the most important limiting factors in wheat and barley productions in the different sites under review included physical properties of the soil such as soil texture and soil depth. Hailu et al. (2015) also reported that OM content, total nitrogen, soil drainage, slope and flooding hazard were found to be the limiting factor for wheat and barley crop production. In other similar study also reported that altitude, soil depth, erosion and surface stoniness were the main limiting factor for wheat and barley production in Wadla Delanta district (Nahusenay & Kibebew, 2015). Specifically, for wheat production, the main limiting factors for wheat production were surface coarse fragments percentage topography and organic carbon (Rabia, 2012) and fertility status (exceedingly available P) (Girma et al., 2015).

Suitability classification for management planning
For detail management planning, the suitability of a land for particular purpose is also important (Kang & Tripathi, 2014;Rossiter, 1994). A single land unit can provide multiple purpose or crops shared similar land unit (FAO, n.d.). As shown in Table 12, 281.88 ha (80.7%) of the land share the same land class for both crops for different land classes (S1 or S2 or S3 or N2); only very small part of the land 2.43 ha (0.51%) was highly suitable for both crops. Furthermore, 266.55 (77.42%) of the land was moderately suitable for both crops. Similarly, 9.3 ha (2%) and 3.6 ha (0.77%) of the land were marginally suitable and not suitable for both wheat and barley crops, respectively. The rest 191.38 ha (19.3%) share different suitability classes for wheat and barley crops; 2.8 ha (0.60%) of the land was highly suitable for barely but moderately suitable for wheat crop. On the other hand, 4.05 ha (0.86%) of the land was highly suitable for wheat and moderately suitable for barely. On the other hand, about 24.70 ha (5.22%) of the land was moderately suitable for barely and marginally suitable for wheat crop. There was 0.4 ha (0.09%) of land which is marginally suitable for barely and highly suitable for wheat crop and 59.53 ha (12.57%) of the area was marginally suitable for barely and moderately suitable for wheat (Table 12).  Figure 14a showed land unit suitability map for wheat and barley crops and as compared to the current land use (Figure 14b), most land unit was not used based on its suitability classes. Most farming practised in marginally and moderately suitable areas for cultivation of these crops. This finding supported by Yohannes (1989) reported that most of cultivation practised on shallow soil depth (<40 cm), which contrast with FAO general recommendation for favourable cultivation (50-60 cm) in the study area. On the other hand, the present small land parcel which claim to be highly suitable for cultivation occupied by eucalyptus-dominated forest land, which might bring a risk of water shortage (Gindaba, Rozanov, & Negash, 2004) and yield reduction (Chanie, Collick, Adgo, Lehmann, & Steenhuis, 2013) due to resource competition and allopathic impact of Eucalyptus trees.
Management measure needed to improve the suitability of land units particularly for marginally suitable lands. As indicated by FAO (n.d.), a particular land mapping unit might be classified as marginally suitable for arable farming and highly suitable for forestry. Therefore, depending on the target of management, decision could be made. Alemu et al. (2013) reported that soil conservation measures are playing an important role in improving the ecosystem services in the form of improving suitability of the soils for most crops. Consequently, different soil and water conservation measures are important to avoid further degradation. Since early 1980s, there were efforts to conserve soil and water through Fanyaa juu, soil bund, grass strip and terracing but not successful (Engda, 2009;Holden, Shiferaw, & Pender, 2005). This situation also confirmed during field observation. There were stone bund terraces in some part of the watershed and grass stripes only established in the upper part of the watershed. In addition, currently uncontrolled grazing practised dominantly in the study area; this was also reported by Engda (2009) in this study area. Teshome (2014) indicated that overgrazing could be one of the major cause for land conversion and fragmentation.

Conclusion and recommendations
Agricultural land suitability of the watershed comprises areas which are highly suitable (S1), moderately suitable (S2), marginally suitable (S3) and permanently not permanently suitable (N2) for wheat and barley crops. Most part of the watershed land is moderately suitable and only small portion are highly suitable, which is without significant limitation to cultivate wheat and barley. Most of marginally suitable land with severe limiting factors such as soil depth, texture, slope, temperature and erosion hazards located in the lower part of the watershed and dominated by the current farm land. This indicates that there is a mismatch between land use and suitability classes (requirement). This situation has become worse as the current resource management practice continues without any improvement. It is recommended that the current land-use pattern need to be modified according to its suitability classes and plan land management measures accordingly. To achieve successful land management, introducing alternative income generation mechanism could help the poor farmers who are totally depend on fragile land for many years. This also possibly will reduce the pressure on the land and provide an opportunity for rehabilitation and restoration of the area.    C < 60s: Clay less than 60% blocky structure; SiCL: silty clay loam; CL: clay loam; Si: silt; L: loam; SC: sandy clay, SCL: sandy clay loam; C > 60v: clay greater than 60% vertisol structure; SL: sandy loam; LS: loamy sand; S: sand; Cm: massive clay; SiCm: massive silty clay; W: well drained; MW: moderately well drained; I: imperfectly drained; P: poorly drained; VP: very poorly drained; E: excessively drained; v: varied.

Citation information
Annex 3. Questionnaire example to show pairwise comparison for the relative importance of the criteria with respect to topography, climate, soil and socio-economic for suitability of a land for major crops (circle the number that mostly represents your viewpoint) Column 1 Column 2