The Rapid and Participatory Assessment of Land Suitability in Development Cooperation

: Development cooperation in agriculture aims to contribute to the achieving of a large part of the Sustainable Development Goals (SDGs) of the UN 2030 Agenda, especially the ﬁrst three, No Poverty (1), Zero Hunger (2), and Good Health and Well-being (3). Development cooperation in agriculture tries to help local communities to increase their awareness, participation, and skills about the management of land and environmental resources, in order to realize sustainable development. In this context, methods of participatory assessment of land suitability have been widely and successfully applied. The present research took place in the framework of a real development cooperation intervention in Nord Kivu (Democratic Republic of Congo) and aimed to implement a rapid participatory assessment of land suitability. In this context, where ofﬁcial and detailed data are not available, the study fostered the active involvement of local experts and used geographical information systems (GIS) to identify the most suitable crops to be supported in different zones of the study area. Toward this aim, the authors used a procedure based on the following steps: the identiﬁcation of relevant land use types (LUTs), mapping capability factors, describing the responses of each LUT to the different capability factors, mapping potential land suitability for the LUTs, mapping accessibility, mapping land suitability for the LUTs. Resulting maps and tables were used to identify the most suitable areas for the different uses. Globally, forestry was the most suited use (99.6% of the study area is potentially highly suitable), followed by the cropping of manioc, sorghum, banana, oil palm, bean and cattle grazing in decreasing order (62.6% of the study area is potentially highly suitable for grazing). When accessibility is considered, forestry presents the largest decrease in the class of high potential suitability ( − 34.9% equal to a loss of 24,945.5 ha), while less adaptable uses, such as cattle grazing showed lower decreases in highly suitable class ( − 11.2%) and larger increases in scarcely suitable class (+9.5%). At a later stage, the comparison between computed suitability and actual land use helped with identifying the areas where forestry should be the only (or most) supported activity and the areas where to push integrated land uses. Our interpretation of the results allows us to recommend the adoption of agroforestry and intercropping as the main methodologies to integrate multiple aims such as the environmental conservation and the improvement of livelihoods.


Introduction
Land suitability assessment is a method of land evaluation that enables decision makers to determine which type of land use is most appropriate for a particular location [1] and to develop accordingly a crop management system for increasing land productivity [2], in order to enhance the ecosystem services provided [3]. In recent years, land-use suitability has gained the widest diffusion as a decision support system in the evaluation of agricultural land [4], in the determination of land habitats for animal and plant species [5], in landscape evaluation and planning [6], and in regional planning and environmental impact assessment [7][8][9].
If land suitability assessment is used in participatory land use planning, it may be useful for avoiding environmental conflicts by the segregation of competing land uses [10][11][12][13]. When land suitability is used to inform evidence-based decisions about land use and consequent socioecological evolution of a given territorial system, land suitability assessment have to take in consideration both agronomical and logistic factors. In fact, land accessibility analysis allows to understand the difference between the potential, i.e., agronomical, land suitability and the real use potential of land in rural areas [14,15]. This is especially useful in agricultural development cooperation (DC) where initiatives normally take place in rural communities where the management of resources such as land or water is to be improved in order to augment the wellbeing of the population [16].
In such a specific context, participatory land suitability evaluation, with the use of geographical information systems (GIS) technology, has been applied widely and successfully [17][18][19][20][21], mainly for practical and structural reasons. Firstly, agri-environmental and regional planning issues are the typical targets of participative and multicriteria decisionmaking procedures as they involve community resource planning and multidisciplinary knowledge [22]. Consequently, DC practitioners need tools that helps to involve multiple stakeholders, or multiple disciplines, and to integrate the different point of view existing in local communities. Participatory GIS grant the possibility to shift from, and link, local scale perceptions to wider perspectives, facilitating stakeholders' interaction and avoiding misunderstanding and conflict over strategic decisions [23]. Secondly, from a DC structural perspective, existing guidelines demand that DC practitioners adopt planning and evaluation methodologies that support evidence-based decision making at all levels of the results chain. Nevertheless, DC practitioners are usually forced to decide about land uses in scarcely informed environments [24][25][26][27][28][29], where the only available source of knowledge about the land is empirical and is detained by the local communities, may them be the farmers or the local technicians. In addition, from a more technical point of view, most of the land suitability classification systems and mapping are quite demanding in terms of data. In several DC intervention zones, the lack of available data and the time constraints of program cycle management determine an increased need for tools that are able to extract structured information form the perception of the actors of the agricultural system [30].
The present research took place in the framework of a real DC intervention in Nord Kivu (NK), Democratic Republic of Congo (DRC) and aims to implement a rapid and participatory assessment of land suitability that (i) overcome context-related data-shortage, (ii) provide an open-to-participation tool for decision-making and (iii) take in consideration the accessibility of rural areas in land suitability evaluation. Present paper resumes the results of the integration of participatory expert opinion elicitation and GIS for identifying most suitable crops to be supported by the DC project in different zones of the study area.

Case Study
Present study occurred in the frame of an EU-funded DC project called ARDST "Appui au retour de réfugiés et déplacés par le biais de la sécurisation de terres en Diocèse de Goma" ("Support for the return of refugees and displaced persons by securing land in the Diocese of Goma"), led by Caritas Development Goma NGO. The project took place within the European strategy for peace seeking and keeping operations in the African Great Lakes region. The Protocol on the Property Rights of Returnees signed at the International Conference of Great Lakes region (ICGLR) on 30 November 2006. Activities were carried out according to the European Development Fund National Strategy for the DRC, focusing on the in-creased diffusion of the state-of-law and augmented sustainability of the agricultural sector [25]. The research aims to contribute to the evaluation of the agricultural development perspectives and namely to the identification of the most suited land use classes (LUTs) in the Diocese of Goma.

Study Area
The study area corresponds the Collectivity of Katoyi, covering approximately 720 km 2 , and located in the southern part of Masisi territory, NK province, eastern region of the DRC (centroid at LAT 1 • 39 14.03 S, LON 28 • 43 18.60 E, Figure 1). NK is concerned with chronic insecurity due to the presence of several active armed groups [31]. Masisi territory, and Katoyi Collectivity in particular, were the theatre of huge population displacement from 1994 to 2004. During this period, Katoyi Collectivity hosted the headquarters of the main armed group of that period, the FDLR ("Forces démocratiques de libération du Rwanda", "Democratic Forces for the Liberation of Rwanda") [32]. Due to their presence, most of the population fled and left the collectivity scarcely inhabited during the following ten years. Since about 2015, returning internally and externally displaced People (IDP and EDP, respectively) are slowly resettling in the collectivity of Katoyi [33,34].
According to the World Bank [35], NK shows one of the highest population densities in Africa (6,600,000 inhabitants, 110 inhabitants/km 2 ). The primary sector (agriculture, livestock, fishing, forestry, mining) is determinant for the economy of NK, reppresenting about 49.7% of the provincial GDP and employing about 80% of the active population [36]. Poverty rate is high in the area (52.4% of the population in NK live below the poverty line, which grows to 69% in rural areas [37]. According to a previous study [38], the land use/cover in Katoyi shows that 52.2% of the land (37,248 ha) is covered by forest and used for the extraction of timber and other non-timber forest products, 30.5% (21,787 ha) is used for agriculture, and 17.3% (12,418 ha) is dedicated to pastures. Two production systems exist in this setting: (1) small farmers practicing subsistence family agriculture on small plots and (2) large landowners who mainly grow cash crop plantations or staple crops in monocultural cropping systems for business, or breed cattle in extensive systems [38]. Both systems practice rain-fed agriculture and extensive grazing.
In Masisi, about 50% of the rural households have plots of less than 0.3 ha, and only 20-25% have more than 2.5 ha of land [39]. It was evident to the authors during the numerous visits in the rural areas, that land scarcity in the area is exacerbated where IDP or EDP come back to lands they fled from and mostly where large-scale cattle ranches slowly increase their size at the expense of local communities.
The authors' previous surveys and observations allow saying that Katoyi hosts large pastures and forests. Forested areas host some artisanal oil palm (Elaeis guineensis) exploitations. In terms of agriculture, the main staple crops are sorghum (Sorghum vulgare), cassava (Manhiot esculenta), taro (Colocasia esculenta), sweet potatoes (Ipomoea batatas) and potatoes (Solanum tuberosum). Many farmers practice multiple cropping combining staple, fruit and vegetable crops. Small mixed-plantations of banana (mainly Musa acuminata) host other vegetable species such as bean, cabbage, spring onion, carrot and tomato. Previous studies [39] confirm that most traditional crops in Masisi territory are bean, banana, sweet potatoes, green peas and maize.
The orography in Katoyi Collectivity is very complex. The NK province is part of an agroecological zone characterized by a long growing season of nearly 210-365 days per year [40], but the heterogeneity of the relief in Katoyi brings a wide variety of microclimates. Generally, the Goma diocese has a mountain climate, but with a close correlation between altitude and average temperature: the average minimum monthly temperature is 23 • C below 1000 m AMSL (above mean sea level), 19 • C at 1500 m AMSL and 15 • C above 2000 AMSL [41]. Average rainfall varies between 1000 mm and 2000 mm [41]. The lowest monthly precipitation occurs between January and February and between July and August. Four seasons characterize the territory of Masisi: two wet seasons and two dry seasons. The first wet season is between September and December and the second is practically from February to June. As for the two dry seasons, they are very short. The first is observed between July and August and the second in January. The mountain climate and the relief give the soils of Masisi a certain complexity. In Katoyi Collectivity, the soils are very deep and rich in humus. They are mildly clayey and not very compact and contain large reserve of organic matter on the surface, especially where land is covered by forests. In fact, the natural vegetation tends to heterogeneous mountain rain forests.

Research Workflow, Data Sources and Methods
Land suitability mapping needs a first set of information about the study area, which we divided into land capability factors (LCFs) and logistic factors (LFs), and a second set of information related to the effects of each of these factors on the locally relevant land use types (LUTs).
In order to obtain such information in an area characterized by a limited availability of official and reliable data, we defined the six-step procedure presented in Figure 2. The workflow is based on the integration of literature, local expert elicitation methods [42,43], and specific GIS overlay methodologies to collect and merge geo-referred useful information [44]. The experts participating in the land suitability assessment were selected within the ARDST project team and counted five people: three local agronomists, one expert in rural development and one pedologist. Additionally, five farmers were involved in a collective interview during the process.

Identifying Relevant Land Uses
We started the analysis from the information gathered and made available in the framework of the ARDST project, namely by analyzing the field visits reports and the results of a survey aimed at collecting information about existing crops and practices in the Diocese of Goma [45]. Seven most relevant land uses for the whole "petit-nord" region were cross-validated through literature review and expert elicitation and selected to form the list of target crops, hereafter called land use types (LUTs) ( Table 1).

Mapping Land Capability Factors
Because of the lack of a land capability map and information for the study area, we firstly compiled an extensive list of land capability factors found (LCFs) in literature [46][47][48][49][50][51]. The extensive list was then analyzed through experts elicitation devoted to the identification of available sources of data for each LCF. This phase, held in July 2018 in Goma, ended with the validation of a final list of three LCFs for the specific study area (Table 2).

Soil factor
Soil quality (fertility) classes "Dorsale du Kivu" Soil map. Scale 1:500,000 [52] A proxy was defined for each LCF as a variable that approximated the measure of the LCF (e.g., elevation used as a proxy of climatic factor, in order to replace unavailable information such as temperature, humidity, precipitation data).
In the study area, the climate is mostly influenced by altitude, which macroscopically determines usual cropping and other LU choices. Elevation data are obtained from the ASTER (advanced spaceborne thermal emission and reflection radiometer) GDEM (global digital elevation model) version 2, with a pixel resolution of 30 × 30 m. In particular, two altitude zones are usually taken into consideration by local experts, namely, the highlands (indicatively located above 2100 m AMSL) and the lowlands (located below 1600 m AMSL). Building on available relevant studies in the region [53], which confirmed the local experts opinions, we defined three elevation classes ( Figure 3A):
Regarding the topographic factor, we used slope classes as a proxy. Slope data were obtained by computation from the ASTER GDEM version 2 (pixel resolution: 30 × 30 m). The study area is characterized by a continuous series of mountains and high hills. Planes are nearly absent and cropping activities, especially subsistence-oriented ones, are implemented no matter the slope. Based on existing studies focusing on the region [53], three slope classes were defined for land capability assessment ( Figure 3B): • low: slope inferior to 6%, • medium: slope between 6 and 25%, • high: slope over 25%.
Regarding the soil factor, quality classes indicating soil fertility have been used as proxy. Data were obtained from the only available soil map of the region in a 1:500,000 scale [52]. The interpretation of existing soil units descriptions was obtained through local expert elicitation. According to available information, the study area is characterized by 13 different soil types, mostly originating from volcanic parental material. We defined three soil quality classes by evaluating soil depth/age and dominant texture as described in the available sources ( Figure 3C): • low, characterized by yellow or shallow brown soils of variable texture (mainly sandyloamy or sandy-clayey) with weak A1 horizon in hilly or uneven topography. • medium, characterized by light and brown soils, mainly sandy-clayey or clayey, with pronounced A1 horizon in uneven topography. • high, characterized by light and brown clayey or clayey-sandy soils, with pronounced chernozemic A1 horizon in a flat or nearly flat topography.
The thresholds used for classifying the study area according to LCFs and build capability classes were determined with reference to existing literature and local expert elicitation. The land capability assessment allowed for drawing three LCF maps using ESRI ArcGIS 10.8 desktop software ( Figure 3). Regarding elevation, low and medium classes are almost equally present and together occupy 76% of the study area. Regarding the slope, the low class is almost absent (only 309 ha); more than 60% of the study area is characterized by a high slope. Regarding the soil quality, the low class is almost absent (only 850 ha), and the medium class occupies 55% of the study area.

Defining the Suitability Scores for Each LUT
The goal of this phase was to identify the influence of LCFs on the different LUTs defined. Each of the LCFs mapped is present in the study area in three different classes: high, medium, and low. The question to be answered was: How does the presence of a particular LCF class affect the productivity of a specific LUT?
To answer this question, we performed an expert elicitation, organizing an ad hoc collective interview with farmers and experts in Katoyi in August 2018. We asked the experts and the farmers to give a suitability score (SC), on a five-point scale (1 to 5, where 1 is 1 is low suitability and 5 is high suitability), for each combination of LUT and LCF class.
During the collective interview, six farmers collectively filled one survey form, while the experts filled one form each. The scores were synthetized though participatory discussion and final scores were validated for each LUT (Table 3). The potential land suitability for each of the seven LUTs considered was assessed through the topological overlay capability of geographical information systems that allow producing a composite map with new polygons generated from the intersection of the original ones coming from the LCF maps.
The new composite map contains all the possible combinations of the three LCFs and their class. For each LUT, the suitability score indicated in Table 3  For each LUT, the total suitability score of a specific area (TSC xzi ) corresponds to the arithmetic average of the three LCF suitability score assigned to the specific LCF classes of that area: where -TSC xzi is the total suitability score of the i-th area, for the x-th LUT, for the z-th LCF class, -SC xz (SQ) is the suitability score, for the x-th LUT, of the the z-th class of soil quality, -SC xz (SL) is the suitability score, for the x-th LUT, of the the z-th class of slope, -SC xz (EL) is the suitability score, for the x-th LUT, of the the z-th class of elevation.
In order to produce the final potential land suitability (PLS) map for each LUT (Figure 4), the total suitability score of a specific area (TSC xzi ) has been reclassified in the following three classes: -High, greater than 3, -Medium, greater than 2 and up to 3, -Low, up to 2.

Mapping Accessibility
Accessibility is considered one important land-use determinant, especially in agroforestry activities [14,15,54]. The approach used was based on the "access time criterion", suggested for the first time in Italy by Hippoliti [55] and recently updated by other researchers, such as Grigolato et al. [56]. Access time is the time required for a forest/agricultural worker to make a round trip on foot from the nearest road to a given field. This time is based on an assumed average walking speed of 4 km/h on flat terrain, and 400 ma/h on steep terrain, where ma is the difference in elevation in meters [57].
Based on this criterion, we calculated a logistic factor (LF) through a cost-distance method [57], thanks to the GIS specific tools. This tool can determine the weighted distance from each cell to the nearest source location across continuous surfaces using grid data. In this study, the "nearest source locations" have been the roads and the weight (or cost) has been the slope (Figure 5), so that the "cost distance" has been the elevation in meters: where CD i is the cost distance, in meters, of the i-th cell from the nearest road, D i is the dimension, in meters, of the i-th cell (in the study the pixel is 30 m), W i is the weight, equal to the slope of the cell (in %). Based on the average walking speed defined above [57], we calculated five accessibility classes (AC) ( Figure 6): • High, with a CD i less than 200 m, that is 30 min walking, • Medium-high: with a CD i of 200-400 m, that is 1 h walking, • Medium: with a CD i of 400-600 m, that is 1.5 h walking, • Low: with a CD i of 600-800 m, that is 2 h walking, • null: with a CD i greater than 800 m, that is more than 2 h walking.
We used the logistic factor to differently weight the potential suitability according to the distance and elevation gain from existing major roads. Figure 6. Accessibility map.

Mapping Land Suitability
The integration of the potential land suitability with the logistic factor due to the accessibility allowed us to produce the final land suitability (LS) map (Figure 7). To this end, we used the overlay mapping tool with ESRI ArcGIS 10.8 desktop software. Land suitability was calculated as follows: where LS i is the land suitability of the i-th area, PLS j is the potential land suitability of the j-th area, AC k is the accessibility class of the k-th area, calculated as follows: • Accessibility Class High = 1, • Accessibility Class Medium-high = 0.8, • Accessibility Class Medium = 0.6, • Accessibility Class Low = 0.4, • Accessibility Class null = 0.2.
In this way, even if the AC is null (CD i greater than 800 m, that is more than 2 h walking), we have no area with land suitability equal to zero.

Results and Discussion
As shown in Table 4, the study area presents many territories with high potential land suitability for most of the LUTs considered. Among the seven LUTs, forestry is the potentially most suited land use (99.6% of the study area is potentially highly suitable). This was expected because forest represents the climax vegetation of the region; consequently forestry is a potentially ubiquitous activity in the study area. This is a well-known characteristic of the area, also confirmed during the elicitation phase, in which forestry registered the second lowest standard deviation of the SC given by the local experts compared with the other LCFs (see Table 3). Among the crops, manioc is the most suitable (85.9% of the land is highly suitable while only 0.5% in scarcely suitable). The result is congruent with its high diffusion and in line with the opinion expressed by the local experts (Table 3).
Sorghum and banana show a high potential suitability in about 62% of the study area, while oil palm and bean are less suited, with only 36.1% and 17.5% of the land highly suitable, respectively. These results are congruent with locally conducted interviews, which show that there are other more suitable areas than Katoyi for oil palm cropping, while bean is ubiquitous due to its value in the framework of subsistence farming [45].
Finally, cattle grazing was potentially highly suitable in the 62.6% of the study area, even if it received more variable scores when evaluated by the experts (Table 3), probably due to the low attitude of cattle in lowlands where the occurrence of animal pests is higher (i.e., Glossina spp.).
Concerning the standpoint and the operational nature of DC organizations' interest in knowing land suitability, the real suitability, which incorporates accessibility, is much more interesting. In fact, including accessibility in land suitability assessment helps DC organizations to identify, on the one hand, areas where eventual support to specific land use activities could be implemented with much ease and increased sustainability, from a logistic perspective, and, on the other hand, the best land use to be supported in the most accessible areas. Figure 8 shows the differences between potential and real land suitability in the study area. Due to accessibility, potential high land suitability decreases for all the LUTs considered, from 2.9% of bean (that has few areas with high potential suitability) to 34.9% for forestry (that has, on the contrary, many areas with high potential suitability). Potential high suitability is confirmed when accessibility is high and decrease when land become less and less accessible.  Figure 9 shows where the differences between potential and real land suitability for each LUT occur in the study area.
Forestry presents the largest decrease in the class of high potential suitability (−34.9% equal to a loss of 24,945.5 ha), and the largest increase in medium potential suitability class (+27.9% equal to a gain of 19,944.9 ha). This result is congruent with the specificity of the study area, where the products of forestry include both more and less transportable items. For instance, carpentry wood needs steady access to roads for often long-range transportation, while honey or medicinal substances extracted from managed forest are easily transported even from less reachable lands. Consequently, it is difficult to imagine any areas completely losing their suitability to forestry. As a result, forestry confirms itself as the most suited LUT for the study area (with 64.6% or 46,156 ha with high suitability and 28.4% or 20,253.3 ha with medium suitability), with most suitable lands located across the whole study area (except the west part, the worst accessible area).
Conversely, crops and animal grazing are generally less adaptable than forestry, and this could explain the results showing a more relevant increase of scarcely suitable land class.
As far as manioc is concerned, the accessibility factor plays such an important role to pass all the south-west area and a part of the eastern area, previously with high suitability, to low suitability ( Figure 9). Overall, highly suitable land for manioc shows a −29.4% (−20,999 ha) decrease and a relevant increase in both low and medium suitability classes (+19.0% and +10.4%). These modifications in the suitability are probably due to the fact that even though manioc adapts both to high-and lowlands and is quite insensitive to high slopes and poor soils, it is not suited to inaccessible areas because it has high yields that imply heavy loads for transportation and relatively low market price. Nevertheless, manioc remains the second most suited crop in the study area (with 56.5% or 40,327 ha highly suitable and 32.6% or 23,312 ha medium suitable), especially in the northern and eastern parts of the study area. Banana, sorghum and cattle grazing show a similar decrease in highly suitable areas (−17.7%, −11.2% and −11.2%, respectively) Sorghum and cattle grazing show highly suitable areas especially in the north-eastern part while the central part of the study area presents medium suitability. Differently, banana presents a quite large area with medium suitability in the eastern part of the study area. Banana is an important staple crop in the region, both in terms of sweet and plantain bananas. Banana herbaceous plants are adapted to sloppy and poor-soils highlands.
Sorghum is known as a hardy crop for adapting to extreme ecological conditions such as rising temperature, low soil fertility, drought [58,59]. Sorghum is also a very important crop for ensuring food security in the region [60]. Despite this great importance and the region's agricultural potential, sorghum production is still low, probably due to the reduction in the length of fallow and to the practice of monoculture cereals [61], coupled with smallholders poor access to land and the effects of climate change [62]. Consequently, the cultivation of sorghum has lost many traditional farming areas in the region where it is supplanted by other cereals such as maize and rice, which by contrast are particularly vulnerable to drought and soil fertility defects. Sorghum would therefore be the cereal of choice in regions where the rains are becoming more and more limiting for maize such as the remote and scarcely accessible high and lowlands in the northern and eastern parts of Katoyi. Our suitability assessment confirms this conviction.
The diffusion of animal husbandry, namely cattle grazing, is growing in the study area, especially in the northeast, while traditionally the collectivity was known for its plantations of Cinchona spp. and Tanacetum spp. amidst managed forest areas [63]. Moreover, cattle are known to suffer in lowlands where parasite pressure tend to be stronger than in the highlands.
Palm oil shows a decrease of areas with high potential suitability similar to the previous LUTs (−14.9%), but starting from a much lower area of potential high suitability (25,793 ha that become 15,182 ha considering accessibility). This decrease happens especially in the southwest boundary of the study area, while the high suitability is confirmed in the centralsouthern part of the study area ( Figure 9). In the study area, oil palm is a cash crop rarely grown in plantations. Oil harvest is much more common as an artisanal practice, in easily accessible lowlands with fertile soils, where the production of single naturally occurring plants can be valorized. Our results confirm the information collected during filed visits.
Bean shows the lowest decrease of areas with high potential suitability respect to all the other LUTs (−2.9%) but starting from the lowest amount of areas with high potential suitability (12,496 ha that become 10,447 ha). Land suitability decrease significantly when the accessibility factor is taken into consideration: even in the most suitable part in the east, some large areas with medium (and low) suitability appear. Bean is confirmed the least suited LUT in the study area. Nevertheless, during the field visits, widespread cropping of bean was reported, both about climbing and dwarf varieties. Interviewees declared that bean shows good attitude also in poor soils and remote areas, due to the nitrogen fixing ability and to the common local practice of harvesting already-dried beans, which facilitates long-distance transportation. Climbing bean is more adapted to highlands and to sloppy zones while dwarf bean easily adapts to lowlands and are much more indifferent to slopes.
If we consider all the LUTs together, we are able to identify the most versatile areas, i.e., areas suitable for more than one LUT, where farmers, and their supporters, have a wider set of possible crops to choose upon. More than 56% of the study area is simultaneously highly suitable to at least 2 LUTs (Figure 10).
The areas characterized by the highest multi-suitability are located along the eastern border of Katoyi area. The multi-suitability, then, decreases progressively in the north and the central area.
The target LUTs chosen proved to be actually relevant in the collectivity of Katoyi given that almost the 73% of the study area is highly suitable for at least 1 LUT and only the 5% presents low suitability for all the LUTs considered ( Figure 10). A recent study [38] using remote-sensed images for land use analysis in NK province (DRC) allows for a shorthand comparison of the real land suitability results with the ongoing land uses. Figure 11 illustrates the distribution of agriculture, pasture and forest across the study area in 2018. According to the mentioned study [38], 52.2% of the 71,409.8 ha available within Katoyi boundaries are forested. 30.5% of lands is used for agriculture and 17.3% is dedicated to pasture.
The Agricultural Development Cooperation (ADC) organizations active in the study area are interested in identifying the priority zones where to focus interventions that should address the most suitable LUTs by also taking into consideration the present land use. This attitude answers to the need of increasing the efficiency of the intervention by prioritizing the areas where most suited LUTs are already in place, and of supporting land use changes that minimize environmental impact while maximizing socioeconomical returns.
Given the ecological features of the study area, on the outskirts of the second largest forest on planet Earth, we assume the perspective of forest conservation in order to correctly interpret the results and produce useful pieces of advice to ADC organizations.
To this aim, the comparison between actual land use and the obtained results about suitability are shown in the next figures.
Concerning areas suitable for forestry (Figure 12), the observed land use matches highly suitable lands for forestry in 31.8% of the whole study area (and in 16.2% in the case of medium suitable land). As already mentioned, forestry represents the spontaneous climax land cover in the region, and according to local experts, it does not show any constraint in terms of soil quality, slope, and elevation.  Our suitability analysis helps identifying the areas where forestry should be the only (or most) supported activity and the areas where to push mixed land uses. The identification of the areas in which the actual forest cover matches the suitability to forestry helps achieving the first aim. 61% and 31% (22,720.5 and 11,536.2 ha) of the land currently covered by forest is located in areas with high or medium suitability for forestry, respectively. Figure 12 highlights (in dark green) the lands with high suitability for forestry and that are actually, and "consistently", forested. These constitute 61.0% of the land actually covered by forests (which in turns covers 31.8% of the whole study area). On one hand, these areas are easily accessible and highly suitable for forestry. On the other hand, in these areas, the only intervention consistent with international goals of forest conservation could be related to sustainable forest management and valorization. Obviously, it does not mean that forestry or forest management are actually implemented in all the areas covered by forests, but the results show the high potential for the expansion of these activities in these areas, scattered across the study area but generally served by the road network.
In the land actually covered by forest and classified as medium or low suitable for forestry (green and light green), which amount to 20.4% of the study area, inhabitants could profit from advancing agricultural or pastoral practices by eroding the forest cover. In this case, local needs should be accurately weighted against the internationally set environmental goals related to forest and soil conservation. In fact, soils in forested areas are mainly thin and fragile, and their quality is strictly related to the carbon cycle related to forest cover. In such a pedologic context, in order to remain sustainable, the development of these areas critically needs an agroecological approach to land valorization, which could integrate the intrinsic value of the actual forest cover and the increase of human exploitation through agroforestry or sylvo-pastoralism.
The actually forested land featuring low suitability to forestry may be identified as sound areas where to push agricultural or grazing activities, even if it would be important take into account the suitability to the other LUTs. These areas are located mostly in the less accessible part of the collectivity, i.e., the western part.
Areas highly suitable for forestry but actually used for agriculture and pasture (in red in the Figure 12) amount to 32.8% of the study area, 14,912 ha (20.9%) occupied by agriculture and 8523 ha (11.9%) by pasture), that represent the 68.9% and 68.4% of the land currently used for agriculture and for cattle grazing, respectively. These lands could profitably host mixed production systems involving agro-sylvo-pastoral activities on the same lands, instead of specialized land uses. This is particularly important information in the study area where large portion of the land are owned by few landowners [38]. In fact, large land holdings show higher risk of being devoted to only one use with a specialized monocultural cropping system. Nevertheless, in NK plantations of eucalyptus are very common, due to the role of the species in producing construction wood and charcoal. Eucalyptus forestry is suited to highlands fertile soils, and it is indifferent to long distance transportation, even if it must be said that remote plantations are mostly used for charcoal, which have lower market price, while accessible plantation are easily used for carpentry.
In these areas, located mostly within the previously mentioned consistently forested areas, ADC organization may suggest for partial reforestation under a scheme of agroforestry or sylvo-pastoralism. In fact, interventions supporting the inhabitants in the diversification toward mixed production schemes may introduce and valorize the forestry component and increase the overall productivity. In these areas, and especially due to high slopes, it is strongly recommended to adopt soil conservation techniques such as terracing, key-line design, no-tillage, mulching and cover crops.
Areas with medium or low suitability for forestry that are currently used for agriculture or pasture (orange and pink in Figure 12), amount altogether to 15.0% of the whole study area. These areas are the most interesting for agricultural development but appear to be located in the western part of Katoyi, which is the less accessible and less suitable for animal grazing due to the low altitude. ADC organizations could support the ongoing activities or, according to the specific suitability of lands, push for oil palm cropping.
Considering agriculture (Figure 13), which is the second most widespread land use in Katoyi, and focusing on the whole study area, the land highly suitable to at least 1 crop (among the 5 ones in the considered LUTs) matches the actual agricultural (cropped) areas only in 23.9% of the whole study area (orange in Figure 13). It represents 78.3% of the agricultural land. These areas are scattered across the whole collectivity but mainly located in the eastern and central part of Katoyi. These are, from ADC organizations' standpoint, the areas where to support existing activities with the aim of valorizing the right localization of already existing activities, and namely the vocation of land to the target LUTs. From the opposite perspective, it is encouraging to report that only 21.7% of the land currently used for agriculture present no high suitability to crops (light yellow) and represent only the 6.6% of the study area. These areas represent the ones where ADC organizations could profitably advocate for a diversified production and/or a shift away from the actually occurring crops.
Further comparison between current land uses and land suitability to cropping reveals that 68.6% (25,559 ha) of forests (light green) and 75.6% (9359 ha) of pastures (light blue) are located on areas highly suitable for at least one crop, representing the 35.8% and 13.1% of the study area, respectively. The former is mainly located in the central part of the study area, while the latter in the northern and eastern parts.
Considering the high human pressure on natural resources and landscapes, ADC organizations may be particularly interested in the identification of these areas in order to implement specific interventions aiming the integration of agricultural production and cattle grazing on the existing pastures, or the integration of forest conservation and agricultural expansion through agroforestry. Once more, soil management should build on the consciousness of existing slopes and thin fertile soil horizons, which should be managed according to the best available practices for soil conservation (i.e., terracing, key-line design, no-tillage, mulching and cover crops).
Pastures represents the third land use in Katoyi. According to our analysis (Figure 14), 64.6% (7990 ha) and 24.3% (3001 ha) of the land currently used for pasture presents high or medium suitability to cattle grazing, respectively. Unfortunately, the results also show that 54.7% (11,928 ha) of the land currently used for agriculture and 45.1% (16,810 ha) of the land covered by forests present high suitability to grazing (dark red in Figure 14). These areas, located in the northern and eastern parts of the study area, could be the ones where ADC organizations may be interested in supporting the integration of cropping with more rational pastoral activities. Here, given that grazing often takes place in lowlands with relative lower slope, soil management practices could range from the already mentioned conservation techniques to more intensive ones, including subsoiling and minimum tillage.

Conclusions
This study assessed the land suitability for different agro-sylvo-pastoral land uses in NK (DRC) in the framework of the EU-funded ARDST Development Cooperation Project led by Caritas Development Goma NGO.
In this context, characterized by a general lack of available and updated data and limited time availability due to the DC programs deadlines, DC organizations need decision support tools that are easily applicable, reliable, and scientifically sound. Local stakeholders, i.e., farmers and technicians, who very often retain in-depth empirical knowledge about local land and agricultural systems, must have access to effective tools for decision-making processes.
In this context, a specific methodology for rapid participatory assessment of land suitability has been developed in order to identify the most suitable crops to be supported in different zones of the study area.
Local farmers and technicians have been involved from the earliest phases of the project in the identification of the relevant land uses (LUTs) of the study area and in the assessment of the productive response of each LUT to the different capability factors and their proxies.
They have also been involved in the quantification of the thresholds for the assessment of accessibility (a very important land-use determinant for agricultural development cooperation in rural communities), based on the "access time criterion".
The use of GIS technology allowed to take into account simultaneously all the parameters involved in the assessment methodology and to integrate all the geographical information produced (thematic maps) in a final overall assessment of the most suitable areas for the different crops.
Finally, actual land uses have been compared with the land suitability in order to propose specific local improvement interventions for the different rural production activities (agriculture, forestry, cattle grazing).
The methodology developed proved to be effective in identifying the areas where it is better to concentrate development cooperation interventions as a priority in order to improve the possibilities of local communities.
In fact, the study shows that about 67% of the study area (47,774 ha) presents a land use type congruent with the land suitability (22,721 ha of forestry in areas with high suitability to forestry, 17,063 ha of crops in areas with high suitability for at least 1 crop, and 7990 ha of pasture in areas with high suitability to cattle grazing).
On the other hand, the remaining 33% of the study area presents a land use type that can be eventually modified for a more rational and sustainable use of environmental resources.
Despite the limitations of the study (related to the different choices we had to make to overcome the lack of data) the results show how development cooperation practitioners can profit from using scientifically sound but still-applicable methodologies to inform strategic evaluation in development planning. Namely, the results obtained have been successfully used in NK by Caritas Development Goma for the strategic evaluation and the increase of sustainability of the agricultural sector. The adoption of the proposed methodology, if supported by technical regulations in DC, may constitute a new standard for improving the sustainability of interventions in line with the European Development Fund National Strategy for the Democratic Republic of Congo, and beyond.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to credit restrictions.

Acknowledgments:
The framework for this study is provided by a three-year EU-funded project called ARDST "Appui au retour de réfugiés et déplacés par le biais de la sécurisation de terres en Diocèse de Goma", led by Caritas Développement Goma NGO. The project started in February 2016 and ended in June 2019. We would like to acknowledge the fundamental role of the whole ARDST project team. Moreover, we acknowledge the role of Giacinto Manfron in the initial conceptualisation of the methodology. The authors acknowledge the support in publishing received from the Department of Agrienvironmental Science, Production, Landscape, Agrienergy of the University of Milan through the APC initiative.