Land Suitability for Co ﬀ ee ( Co ﬀ ea arabica ) Growing in Amazonas, Peru: Integrated Use of AHP, GIS and RS

: Peru is one of the world’s main co ﬀ ee exporters, whose production is driven mainly by ﬁve regions and, among these, the Amazonas region. However, a combined negative factor, including, among others, climate crisis, the incidence of diseases and pests, and poor land-use planning, have led to a decline in co ﬀ ee yields, impacting on the family economy. Therefore, this research assesses land suitability for co ﬀ ee production ( Co ﬀ ea arabica ) in Amazonas region, in order to support the development of sustainable agriculture. For this purpose, a hierarchical structure was developed based on six climatological sub-criteria, ﬁve edaphological sub-criteria, three physiographical sub-criteria, four socio-economic sub-criteria, and three restrictions (co ﬀ ee diseases and pests). These were integrated using the Analytical Hierarchy Process (AHP), Geographic Information Systems (GIS) and Remote Sensing (RS). Of the Amazonas region, 11.4% (4803.17 km 2 ), 87.9% (36,952.27 km 2 ) and 0.7% (295.47 km 2 ) are “optimal”, “suboptimal” and “unsuitable” for the co ﬀ ee growing, respectively. It is recommended to orient co ﬀ ee growing in 912.48 km 2 of territory in Amazonas, which presents “optimal” suitability for co ﬀ ee and is “unsuitable” for diseases and pests. This research aims to support co ﬀ ee farmers and local governments in the region of Amazonas to implement new strategies for land management in co ﬀ ee growing. Furthermore, the methodology used can be applied to assess land suitability for other crops of economic interest in Andean Amazonian areas.

ISPRS Int. J. Geo-Inf. 2020, 9, x FOR PEER REVIEW 3 of 21 km 2 ), Bagua (5861.16 km 2 ), Luya (3100.36 km 2 ), Utcubamba (3972.83 km 2 ) and Rodríguez de Mendoza (3714.24 km 2 ). These last four provinces are the main coffee producers, with approximately 34,000-42,744 hectares, producing around 25,000 tons per year, and with average yields ranging from 700 to 1250 kg/ha [31,34].  Figure 2 shows an analysis of coffee land suitability in Amazonas. In summary, (i) it identified and evaluated the climatological, soil, physiographic and socio-economic criteria as well as subcriteria which influence on coffee growth, and (ii) it elaborated the layers with suitability thresholds for each sub-criterion, whilst also looking at (iii) the importance of the criteria and sub-criteria, which was prioritized and weighted using Pairwise Comparison Matrices (PCMs), and (iv) the maps were overlapped based on the importance weights, and finally (v) the suitability map was restricted to the conditions regarding coffee diseases and pests (restrictions).  Figure 2 shows an analysis of coffee land suitability in Amazonas. In summary, (i) it identified and evaluated the climatological, soil, physiographic and socio-economic criteria as well as sub-criteria which influence on coffee growth, and (ii) it elaborated the layers with suitability thresholds for each sub-criterion, whilst also looking at (iii) the importance of the criteria and sub-criteria, which was prioritized and weighted using Pairwise Comparison Matrices (PCMs), and (iv) the maps were overlapped based on the importance weights, and finally (v) the suitability map was restricted to the conditions regarding coffee diseases and pests (restrictions).

Methodological Scheme
ISPRS Int. J. Geo-Inf. 2020, 9,  Restrictions Data collection and processing Figure 2. Flow chart of the methodology used for the analysis of the suitability of the territory for coffee growing in the Amazonas region, Peru.

Identification and Selection of Criteria, Sub-Criteria and Restrictions
Based on a bibliographic review, the main criteria and sub-criteria that determine the development of coffee growing in Peru and the world were identified [35][36][37][38]. For a detailed description of Arabica coffee and its cultivation, see Kuit et al. [39]. A hierarchy, therefore, was constructed based on 4 criteria (climatological, edaphological, physiographic and socioeconomic) which in turn are broken down into 18 sub-criteria ( Figure 3). Accordingly, as Climatological subcriteria, we have considered the following sub-criteria: average annual temperature, annual mean minimum temperature, annual mean maximum temperature, mean annual rainfall, relative humidity and dry periods; as Edaphological sub-criteria, pH, Texture, Stoniness, Soil Organic Matter (SOM) and Cation Exchange Capacity (CEC); in the same way, as Physiographic sub-criteria: Elevation, Terrain slope and Terrain aspect; and finally, as Socioeconomic sub-criteria: Land Use/Land Cover (LULC), distance to the road network, distance to the water network and Protected Natural Areas (PNA). Furthermore, to apply a principle of prevention, the map of potential land suitability was restricted by layers of the most common diseases and pests of coffee growing in Amazonas [40]. Such layers were based on optimal temperature for the development of coffee leaf rust (Hemileia vastatrix), coffee berry borer (Hypothenemus hampei) and leaf spot (Cercospora coffeicola).

Identification and Selection of Criteria, Sub-Criteria and Restrictions
Based on a bibliographic review, the main criteria and sub-criteria that determine the development of coffee growing in Peru and the world were identified [35][36][37][38]. For a detailed description of Arabica coffee and its cultivation, see Kuit et al. [39]. A hierarchy, therefore, was constructed based on 4 criteria (climatological, edaphological, physiographic and socioeconomic) which in turn are broken down into 18 sub-criteria ( Figure 3). Accordingly, as Climatological sub-criteria, we have considered the following sub-criteria: average annual temperature, annual mean minimum temperature, annual mean maximum temperature, mean annual rainfall, relative humidity and dry periods; as Edaphological sub-criteria, pH, Texture, Stoniness, Soil Organic Matter (SOM) and Cation Exchange Capacity (CEC); in the same way, as Physiographic sub-criteria: Elevation, Terrain slope and Terrain aspect; and finally, as Socioeconomic sub-criteria: Land Use/Land Cover (LULC), distance to the road network, distance to the water network and Protected Natural Areas (PNA). Furthermore, to apply a principle of prevention, the map of potential land suitability was restricted by layers of the most common diseases and pests of coffee growing in Amazonas [40]. Such layers were based on optimal temperature for the development of coffee leaf rust (Hemileia vastatrix), coffee berry borer (Hypothenemus hampei) and leaf spot (Cercospora coffeicola).  The sub-criteria were reclassified and weighted, according to suitability thresholds of the territory, in three categories [20]: "Optimal" (3), "Suboptimal" (2) and "Unsuitable" (1) ( Table 1). We considered the same categories for the restrictions; the "optimal" category for diseases or pests resulted in being "unsuitable" for coffee.  The sub-criteria were reclassified and weighted, according to suitability thresholds of the territory, in three categories [20]: "Optimal" (3), "Suboptimal" (2) and "Unsuitable" (1) ( Table 1). We considered the same categories for the restrictions; the "optimal" category for diseases or pests resulted in being "unsuitable" for coffee.

Resources and Mapping
Data on geographic information for climatological sub-criteria and for disease and pest restrictions, with the exception of relative humidity and dry periods, were obtained from WorldClim 2.1, with average data from 1970 to 2000 and a resolution of~1 km [70]. The Relative Humidity (RH) layer was generated on the Google Earth Engine (GEE) platform. This platform is freely accessible, works using the Google cloud, and is capable of processing petabytes of geospatial data on a global scale [71]. GEE contains climate data (evaporation, humidity, temperature, etc.) from the Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS) GEE at 0.1 • spatial resolution [72,73]. A script using JavaScript coding language was implemented to develop Equation (1) using the required bands as input data. The RH (Equation (1)) was calculated as the quotient between the current pressure (e) and the air saturation pressure (e s ) at an air temperature T (in Kelvin), e s is the full denominator of Equation (1). This equation was derived from Claysius Clayperon's law. (1) The dry periods layer was based on methodology used to quantify meteorological droughts [74]. Therefore, the total annual average rainfall data was divided into five quintiles. Values in the first two quintiles 0-20% and 20-40% are considered dry and "unsuitable", values between 40-80% are considered "suboptimal" and values in the last quintile (80-100%) are considered wet or "optimal".
The data for the edaphological sub-criteria were obtained from system for global digital soil mapping SoilGrids, with a spatial resolution of 250 m [75]. The layers of the depths 0-5 cm, 5-10 cm, and 10-15 cm were all averaged out, given that the first 20cm has the highest concentration of coffee roots [76]. The SOM layer was obtained from the product of the Soils Organic Carbon in percentage (SOC%) and a conversion factor of 1.72 [77,78]. SOC% was obtained by converting the Organic Carbon layer into Ton/ha (SOC) of SoilGrids using Equation (2), where dap is the bulk density of the soil (g/cm 3 ) and SD is the sampling depth (cm).
The slope and aspect of the terrain were derived from the 30-m spatial resolution digital elevation model (DEM) from NASA's Shuttle Radar Surveying Mission (SRTM) [79], by using the "slope" and "aspect" tools of QGIS 3.10. The same DEM was reclassified to obtain the elevation suitability map. Regarding the socioeconomic criterion, a raster layer downloaded from The Copernicus Global Land Service (CGLS) was used for the LULC sub-criterion, at a spatial resolution of 100 m [69]. This layer was complemented with the agricultural LU polygons of the LULC map of the ZEE [31] for Amazonas and of Rojas et al. [80] for R. de Mendoza. The distance to the road network was generated using the 2018 national, departmental and neighborhood road network layers downloaded from the Ministry of Transport and Communications (MTC) [81]. The distance to the water network was generated using the IGN's river and stream layer [82]. These last two sub-criteria were updated to 2020 based on a Google Satellite server image as a basemap, using the QGIS 3.10 on-screen digitizing method. Additionally, we used the PNA and their buffer zones updated to 2020 by the National Service of Natural Areas Protected by the State (SERNANP) [83].
The climatological, edaphological, physiographic, socioeconomic and restriction sub criteria were extracted according to the geographical limit of the study area. Subsequently, they were reclassified in the three suitability thresholds (see Table 1), using the simple reclassify values tool of SAGA in QGIS 3.10 and were re-scaled to 30 m pixels, so as to homogenize the geospatial data, using the resampling tool.

Multiple Criteria Evaluation (MCE) and Analytical Hierarchy Process (AHP)
A PCM was constructed to compare and weight the criteria to each other for each hierarchical group (see Figure 3), and, in addition, a PCM was constructed to compare and weight diseases and pests. The construction of the aforementioned PCMs was based on the weighting scale proposed by Saaty [84] (Table 2). For each set of criteria, a group of 10 coffee growing experts assigned a weighted value. This group of experts was made up of national and international researchers in coffee cultivation, as well as professionals from the agricultural coffee cooperatives in the region. The PCM provided by the expert group were processed (for an example of PCM matrix processing, see [85]), which allows to us to obtain the weight of importance of each of the criteria and sub-criteria; the sum of weights should always be 1 per hierarchical group [85]. Because the PMCs of the experts may have inconsistencies, the inconsistency values that may be acceptable are calculated using the Consistency Ratio (CR) [84]. Therefore, it was considered that in a PCM, the highest eigenvalue (λmax) is always greater than or equal to the number of rows and columns (n) of the PCM, therefore we can calculate the Consistency Index (CI) using Equation (4) [86]. Then, CR is calculated as the quotient between the CI and the Random Consistency Index (RI) [84]. The RI values are conditioned by the "n" value of the matrix, shown in Table 3. Table 3. RI values based on the "n" of the matrix according to Saaty [86]. In order for the consistency of the matrix to be considered as acceptable, the CR must be less than 0.1; otherwise, the AHP may give erroneous results unless the assigned values in Table 1 are re-examined by the coffee-growing experts [87][88][89][90]. All of the above calculations were performed using the Easy AHP plugin in QGIS 2.8. Thereafter, each weighting (importance) of the sub-criteria was multiplied by its respective weight, obtaining a local score, which was multiplied by the weight of each criteria (climatological, edaphological, physiographic and socioeconomic) to obtain the global score and identify the most important criteria. We proceeded using a similar method for diseases and pests of coffee growing.

Obtaining Sub-Models and Suitability Model for Coffee Growing
With the weighted weights calculated and the reclassified maps, the ArcMap 10.5 Weighted Overlay tool was used to obtain the sub-models of climatological, edaphological, physiographic, socioeconomic and disease and pest suitability. At the same time, the first four sub-models were integrated with the same tool to generate the final suitability model for coffee cultivation. This final model was then overlaid with the disease and pest sub-model to identify the territory with these phytosanitary restrictions. Furthermore, the "optimal" plots for coffee, which at the same time are "unsuitable" and "suboptimal" for diseases and pests, were counted according to their area: 0.5-1 ha, 1-2 ha, 2-3.5 ha, 3.5-5 ha and >5 ha. These ranges were established considering that coffee cultivation in the Amazonas region is practiced mainly in plots < 1 ha (family and subsistence agriculture) [91]. Plots < 0.5 ha were considered insufficient to be economically viable. Finally, to validate the final suitability model, it was compared with the spatial distribution of the coffee areas of the ZEE (Figure 1) [31] and with farms georeferenced by the coffee cooperatives and by the Research Institute for the Sustainable Development of Ceja de Selva (INDES-CES) of the Toribio Rodríguez National University of Mendoza (UNTRM). For this, it was assumed that the distribution of current coffee areas is the result of a historical and technical selection for the best possible available area [20]. Table 4 shows the weightings obtained for each sub-criterion, based on the estimates of 10 coffee growing experts. The most important sub-criteria for establishing coffee growth were elevation (13.27%), followed by LULC (8.18%). On the other hand, the least important were relative humidity (3.16%) and annual mean minimum temperature (2.69%). When taking the criteria into account, we found that the following sub-criteria obtained the highest value for each criterion group: average annual temperature (22.37%), pH (29.22%), elevation (53.06%) and LULC (32.72%), whereas, the lowest weightings were given to the annual mean minimum temperature (10.76%), SOM (16.32%), aspect of the land (22.75%) and distance to water network (20.22%). Regarding restrictions, H. Vastatrix was the most important (46.7%), followed by H. hampei (32.0%) and C. coffeicola (21.2%).  Figure 4 shows the reclassified maps for each of the climatological, edaphological, physiographic and socioeconomic sub-criteria, according to suitability thresholds (Table 1). Meanwhile, Figure 5 presents the land's suitability for diseases and pests. Furthermore, Table 5 provides the areas for each sub-criteria suitability threshold and diseases and pests. Accordingly, the sub-criteria with the greatest "high suitability" area in relation to their criteria group included: relative humidity (81.8%), SOM (92.8%), aspect of the land (43.9%) and PNA (70.9%). Conversely, the ones with the largest "not suitable" area were: rainfall (59.0%), CEC (25.2%), elevation (55.6%) and LULC (76.4%). Hence, the SOM is the sub-criterion that most favours coffee growth, while the LULC is in fact the sub-criterion that most restricts the territory. Regarding restrictions, H. vastatrix has the largest area (56.6%) that is suitable for its development in Amazonas; on the contrary, H. hampei has the smallest area (37.6%).

Suitability Sub-Model (Criteria) Maps
With the weighted overlap of sub-criteria, the suitability sub-models were generated for each hierarchical group ( Figure 6). The results demonstrate that the edaphological sub-model (60.75%) has the largest "optimal" area for coffee cultivation in Amazonas (Table 6). Nevertheless, the climatological sub-model (3.43%) has the smallest "optimal" area. Most importantly, 37.64% (15,859.40 km 2 ) of Amazonas is not suitable for coffee cultivation, due to its "optimal" conditions for the development of diseases and pests.

Land Suitability Models for Coffee Growing
With the weighted superposition of climatological, edaphological, physiographic and socioeconomic sub-models, we generated the model of land suitability for coffee growing (Figure 7a). In fact, in Amazonas, 11.4% (4803.17 km 2 ), 87.9% (36,952.27 km 2 ) and 0.7% (295.47 km 2 ) are "optimal", "suboptimal" and "unsuitable", respectively, for coffee growing (Table 7). At the provincial level, Rodríguez de Mendoza (1130.93 km 2 ), Bagua (913.99 km 2 ) and Condorcanqui (913.90 km 2 ) present the largest areas with "optimal" suitability for coffee growing, while Luya (295.73 km 2 ) and Chachapoyas (282.70 km 2 ) present the least "optimal" areas. Then, with the weighted superposition of the territory suitability model for coffee growing (Figure 7a) and the suitability sub-model for diseases and pests (Figure 6e), the territory suitability model for coffee cultivation with restrictions was generated (Figure 7b). It was found that, 7.6% (1208.18 km 2 ), 18.6% (2628.00 km 2 ) and 8.0% (966.92 km 2 ) of the "optimal", "suboptimal" and "unsuitable" territory, respectively, for the development of diseases and pests, is "optimal" for growing coffee. At the provincial level, Condorcanqui (10,217.15 km 2 ) and Bagua (4059.37 km 2 ) present the largest areas with "optimal" suitability for the development of diseases and pests, while Luya (71.38 km 2 ) and Chachapoyas (15.83 km 2 ) present the least "optimal" areas. Of the territory of the Amazonas region, 6.25% (2628.00 km 2 ) is considered "very optimal", since it has "optimal" suitability for coffee and "suboptimal" suitability for diseases and pests, while 2.30% (966.92 km 2 ) is considered "excellent", since it has "optimal" suitability for coffee and is "unsuitable" for diseases and pests.     Of the "very optimal" and "excellent" territory for coffee cultivation, 0.3% (3.17 km 2 ) and 0.1% (2.52 km 2 ), respectively, have plots with areas < 0.5 ha, so they are not large enough to be economically viable (Table 8). In addition, 75.05-98.90% of both "very optimal" and "excellent" plots have surfaces > 10 ha, in all provinces. Table 8. Area of "very optimal" (VO) and "excellent" (Ex) plots for coffee, according to provinces in the Amazonas region, Peru. The results of the validation show that there is a close correspondence between the patterns of the areas modeled as "optimal" for the cultivation of coffee and the coffee areas of the ZEE and the farms georeferenced by the coffee cooperatives and by the INDES-CES of the UNTRM (Figure 8).

Discussions
Differentiating from previous studies of the suitability of the territory for the growing coffee [15,[20][21][22][23]29,30,55,92], this study included a greater number of sub-criteria. Furthermore, only three [21,29,92] of the previous studies integrated AHP to rank and weigh the importance of the subcriteria. In studies of the suitability of the territory for coffee growing, other crops [24][25][26] or agriculture in general [13,62,63], it is common to consider groups of climatological, edaphological, physiographic and socioeconomic criteria. However, this study included one group of criteria, that is, the most incidental disease and pest restrictions for coffee. These restrictions allow greater

Discussions
Differentiating from previous studies of the suitability of the territory for the growing coffee [15,[20][21][22][23]29,30,55,92], this study included a greater number of sub-criteria. Furthermore, only three [21,29,92] of the previous studies integrated AHP to rank and weigh the importance of the sub-criteria. In studies of the suitability of the territory for coffee growing, other crops [24][25][26] or agriculture in general [13,62,63], it is common to consider groups of climatological, edaphological, physiographic and socioeconomic criteria. However, this study included one group of criteria, that is, the most incidental disease and pest restrictions for coffee. These restrictions allow greater precision in determining the territory highly suitable for growing coffee in Amazonas region. It should be noted that the criteria used in the analysis of the suitability of the territory depend on the focus of the research and the availability of spatial data. For example, future studies may include cost-benefit, productivity, crop rate of return, costs of the land use changes, population benefited by the crop, or other economic sub-criteria [93]. However, the main problem lies in the lack of spatial data for the sub-criteria. In Peru, specifically in the Amazonas region, computer resources are scarce, even more so when it comes to specific studies of biological, environmental and social criteria at detailed local scales [94].
In terms of the development of coffee growing in Amazonas, our study shows that socio-economic (28.35%) and climatological (28.31%) criteria play a major role in coffee growing, followed by edaphological (25.03%) and physiographic (18.31%) criteria. However, for Mighty [21] and y Rono and Mundia [29], climatological criteria are the most important, followed by edaphological, socioeconomic and physiographic criteria. In contrast, although Ochoa et al. [92] they did not analyse socioeconomic criteria, in their study, the edaphological criteria are the most important, followed by the physiographic and climatological criteria. In fact, the number and the different sub-criteria used by each group of criteria, the local reality and the experience of the experts, contributed to differentiating the importance of the criteria. Of course, in our study, the slightly greater importance of the socio-economic criteria over the climatological and edaphological criteria concerning Amazonas may be due to the conservationist approach of the experts. Indeed, of the socio-economic sub-criteria, foremost were LULC (32.72%) and PNA (25.12%). Therefore, both are important here, as agriculture (mainly coffee and cocoa) and livestock farming have driven the greatest changes in LULC, making the region one of the most affected by deforestation [80,95]. In this context, public and private PNAs have been the main strategy for in situ conservation, highlighting that Amazonas is also one of the regions to record the highest number of PNAs [83].
In relation to climatological sub-criteria, the average annual temperature (22.37%) is the highest, followed by rainfall (21.4%); in a similar way, for Mighty [21] and Rono and Mundia [29], the temperature is the most important climatological sub-criterion; for their part, Ochoa et al. [92] state that precipitation is the most important. In fact, coffee crops are highly sensitive to climate, especially to the temperature at which they are grown [96,97]. Additionally, in this type of study, the suitability models are generated under current environmental conditions and must be constantly updated, since climate change is a factor that influences the patterns of climatological sub-criteria [15]. This is important for the agricultural sector in Latin America and the Caribbean because it is one of the most affected by climate change [7].
With regard to the edaphological sub-criteria, Akinci et al. [13] identified soil type as the most decisive factor for agriculture in general, and stoniness and drainage as the least important. Similarly, in our study (coffee growing), for edaphological sub-criteria, soil pH (29.22%) is indeed the most important, being followed by texture (19.06%); this also agrees with Ochoa et al. [92]. Conversely, for Rono and Mundia [29], soil pH is the least important. In fact, soil pH is key for the assimilation of nutrients by the plant, due to the fact that the availability of nutrients varies according to the levels of acidity; more importantly, in very high acidity, some minerals turn out to be toxic for the plants [51].
According to Rono and Mundia [29], elevation is the most important of the physiographic sub-criteria; conversely, Ochoa et al. [92] consider slope to be the most important; however, Mighty [21] consider that there is no significant difference (0.10%) between elevation and slope. In our study, of course, elevation was the most important (53.06%), followed by slope (24.18%) and aspect (22.75%). In addition, like Lara et al. [20], we integrated the aspect of the terrain as a physiographic sub-criterion, because to the fact that the aspect influences sunlight reception and temperature variations at the microclimate level [61]. In two locations in Costa Rica, Avelino et al. [98] found that crops on the easterly-facing slopes produced better quality coffee. The experts, therefore, highlighted the elevation from the point of view of coffee quality, because Duicela et al. [99], determined that the coffee varieties grown at higher elevation present higher organoleptic quality. This altitudinal factor has also contributed to the migration of coffee plantations to higher lands, and consequently more deforestation in the main coffee provinces of the region [80].
Regarding the socio-economic criteria, the distance to roads (21.94%) and distance to the water network (20.22%) have no significant difference (0.43%) for this study, closely related to the findings presented by Mighty [21] (1.40%). Definitely, in our study, the coffee experts prioritize the proximity to a road rather than a water body, in view of the fact that Amazonas has a privileged location in terms of water resources, while it lacks good road access and transportation for the arrival of inputs and the output of the product to the markets.
In the Amazonas region, the area modeled as "excellent" (912.48 km 2 ) for coffee cultivation is 113.5% greater than the area reported in the last national agricultural census of 2012 (427.44 km 2 [34]). Although the areas modeled as "excellent" for coffee cultivation are in line with the reference spatial distribution for validation (coffee areas in the ZEE [31] and geo-referenced farms), there are large "excellent" areas not considered in the ZEE and there are also some georeferenced farms that are not located in "excellent" areas. This is due to the fact that the ZEE is an old reference map (2010) and to the lack of a national (or regional) map of the current distribution of coffee. Additionally, Lara et al. [20] indicate that some validation plantations may be located in less suitable areas as agricultural practices such as planting shade trees improve the suitability of the land when temperatures drop, but coffee plantations should not be located in completely unsuitable areas.
Suitability maps for coffee diseases and pests showed that some areas of Condorcanqui, Bagua, Utcubamba, Luya and Rodriguez de Mendoza have ideal development conditions for coffee leaf rust, coffee berry borer and berry blotch. These provinces (except Condorcanqui) match those reported by the National Agrarian Health Service (SENASA), where they are constantly evaluating and applying control measures for these diseases and pests [40]. Due to this coincidence with SENASA, our temperature-based analysis is reliable. Unfortunately, diseases and pests depend largely on other factors such as relative humidity, photoperiod, elevation, percentage of shade, etc. [65][66][67][68]. Therefore, in future research it is recommended to consider all factors in order to improve the precision in the evaluation of land suitability. Additionally, we suggest using thematic layers with better spatial resolution, given that diseases and pests depend on microclimates [100,101], not well represented in layers with~1 km resolution such as those in WorldClim.
The methodology used in this research is based on the integration of AHP, GIS and RS in order to analyze land suitability for coffee growing; the same can be applied not only to other coffee-growing areas, but also to other crops that have nutritional, economic and environmental importance, adjusting them to their local reality. Furthermore, these approaches have been influential in the field because of land suitability analyses identifying areas in which a crop is suitably developed, thus contributing to not over-exploiting the soil resource and consequently to a sustainable agriculture practice.

Conclusions
For the development of coffee growing in Amazonas region, the most important criteria are socioeconomic and climatological, followed by edaphological and physiographic. The edaphological sub-model has the largest "optimal" area for coffee cultivation in Amazonas region, while the sub-model with the smallest "optimal" area is the climatological one. In Amazonas, 11.4% (4803.17 km 2 ), 87.9% (36,952.27 km 2 ) and 0.7% (295.47 km 2 ) of the territory presented suitability "optimal", "suboptimal" and "unsuitable" characteristics, respectively, for coffee growing. It was found that, 8.0% (966.92 km 2 ), 18.6% (2628.00 km 2 ) and 7.6% (1208.18) present "unsuitable", "suboptimal" and "optimal" territory, respectively, for the development of diseases and pests, is "optimal" for growing coffee. It is recommended to orient coffee growing in 912.48 km 2 of territory in Amazonas, which has "excellent" suitability for coffee.
The integrated use of AHP, GIS and RS allowed the proper identification of areas with potential for coffee cultivation in Amazonas region and aims to guide the proper management of the territory. In addition, a methodological framework of potential application is provided for crops of economic interest in the Andean Amazon areas.