Tree diversity and biomass carbon stock analysis along altitudinal gradients in coffee-based agroforestry system of Western Ethiopia

Abstract Agroforestry systems are thought to reconcile biodiversity protection with food production and as a means of climate change adaptation and mitigation options. The contribution of a coffee-based agroforestry system to tree diversity and carbon stock along altitudinal gradients in Western Ethiopia was assessed. At 500-m intervals, six transect lines were methodically set up throughout the altitudinal gradient. There were made a total of 60 sample plots, each measuring 40 m by 40 m. A total of 34 woody species were identified. Biomass carbon stocks and tree diversity were quantified across altitudinal gradients. In the middle altitude, there were more woody species (28) than in the top altitude, where there were only a few species (16). The tree plants stored around 40.6 t ha−1 of biomass carbon on average. Aboveground biomass had a carbon stock of 32.22 C t ha−1, whereas belowground biomass had a carbon stock of 8.38 C t ha−1. The lower altitude biomass carbon stocks were substantially bigger than the upper altitude, which were 48.4 C t ha−1 and 25.67 C t ha−1, respectively. With increasing altitude, the study found a statistically significant negative link between tree diversity and biomass carbon storage (P < 0.05). The negative link between biomass carbon stock and altitude was that tree parameters that determine the amount of biomass carbon sequestered in a plant, such as basal area, tree diversity, and density, decreased as altitude increased. Despite differences along altitudinal gradients, the systems supported a diverse range of tree species and biomass carbon stocks.


Introduction
Agriculture has become the most vulnerable to the effects of climate change (Mekonnen & Gebreyesus, 2011). As a result, human communities all over the world have devised a wide range of complex tactics for coping with climatic change and the natural environment (Gizachew & Shimelis, 2014). Understanding the true dynamics of climate change consequences at the local level, such as farm households and communities, is crucial (T. Tadesse et al., 2009). Smallholder farmers in a country like Ethiopia, which has a high level of ecological and cultural diversity, such as language, beliefs, and indigenous knowledge of conserving natural resources, are likely to have developed their own indigenous adaption mechanisms Lin, 2020). To combat climate change, conserve biodiversity, and increase food production, subsistence farmers in various regions of the country have developed and adopted the practice of agroforestry (Jamnadass et al., 2013;Reppin et al., 2020).
Agroforestry is an old practice that involves growing trees or shrubs near or amid pastureland or crops in order to benefit from the subsequent ecological and economic interactions Mkonda & He, 2017). They are specific local examples of methods that are relevant to a certain habitat, plant species, arrangement, management, and socioeconomic function, e.g., Coffee-based home-garden. Home-garden is specific due to its distinctive environment, plant species, arrangement and management, social and economic functions .
The likely importance of agroforestry is designated differently. As a result, due to the localityspecific nature of each agroforestry practice, they have evolved differently across the world. In Ethiopia, where arabica coffee originated, coffee-based agroforestry grows naturally as an understory shrub in native woods (G. Tadesse et al., 2014). It is also grown throughout a wide management gradient, but exclusively under shade (Zewdie et al., 2022). The practice is abundantly found in the southwestern and western parts of the country, where a large proportion of the country's forest cover is found. As coffee shrubs are perennial crops, coffee-based agroforestry practices are believed to form a better biodiversity and biomass carbon sequestration system than other agroforestry practices (Buechley et al., 2015;Hylander et al., 2013;Tesfay et al., 2022). Another benefit of coffee-based agroforestry management is that it slows the rate at which ecosystems with trees are being converted to open agricultural landscapes (Hylander et al., 2013).
Due to the composition of trees and shrubs, the system sequesters a high amount of biomass carbon. It also contributes to maintaining food and nutrition security, helps diversify household income, provides fiber and energy to local communities (Jamnadass et al., 2013). Agroforestry systems have a high potential to offset greenhouse gas emissions compared to land resources converted to other more intensive forms of land uses owing to the carbon storage in trees, shrubs, and soils. In the case of traditional coffee farming, retaining a high degree of the canopy increases the capacity of the system to store more carbon in the tree and perennial crop components (Olsson et al., 2019).
In a wide range of ecological settings, the practices maintain biodiversity and provide ecosystem resilience in addition to biomass carbon stocking (Santos et al., 2022). Tropical and subtropical ecological zones are home to the majority of agroforestry systems . According to different researches, many tropical forests are under severe anthropogenic pressure and require management interventions to maintain overall biodiversity, productivity, and sustainability (Kurnar et al., 2006). Despite the fact that species composition often varied widely, examples from around the world demonstrate that tropical agroforestry systems can support high levels of biodiversity, often similar to native forest (Tscharntke et al., 2011). In such regions, subsistence farmers who care for a diverse range of trees to support their families might hearten carbon storage and tree diversity, which will benefit ecosystems (Sistla et al., 2016).
Understanding the places of convergence between carbon storage and biodiversity protection is one of the primary concerns at the intersection of climate change and ecosystem science (Mulatu & Hunde, 2020). Nonetheless, higher carbon storage potential is related to its co-existence with biodiversity conservation, which varies with the type of agroforestry practices (Santos et al., 2022).
Though carbon stocks and biodiversity in Ethiopian agroforestry systems have been researched (Eshetu et al., 2021;Manaye et al., 2021), there has been no research on coffee-based agroforestry systems across altitudinal gradients in Western Ethiopia. The majority of studies on agroforestry systems have mostly focused on food production and soil fertility management Mulatu & Hunde, 2020. However, the impact of coffee-based agroforestry systems to carbon stocks and biodiversity conservation remains unexplored. For that reason, there is a limitation in addressing the contribution of coffee-based agroforestry systems towards tree diversity conservation and carbon sequestration acrossA agro-ecological differences (Tschora & Cherubini, 2020). Therefore, the specific objectives of this study were to analyze tree species diversity in coffeebased agroforestry systems, quantify biomass carbon stock in tree species, and assess the effect of tree diversity on biomass carbon stock along altitudinal gradients. In the study area, the farmers were mostly coffee producers, but they also raise cattle. In addition, population density was varying across the ranges of the altitudinal gradients. The population of the upper land was higher, so they raise many cattle and tend to the production of other grains besides coffee to support the population density. On the other hand, as the altitude increased, the weather variables of the area was varying and decreased at the upper altitude. Therefore, due to these two reasons, we hypothesized that biomass carbon and tree diversity would change amongst the three gradients of coffee-based agroforestry systems. In addition, because of a projected weather variable-related drop with elevation ascent, the high biomass carbon stock would be recorded at lower altitudes and associated to tree diversity; and, as a result, both biomass carbon stock and tree diversity would decrease with elevation.

Study area description
The study was conducted in the Gimbi district in the West Wollega Zone of Oromia National Regional State, Western Ethiopia ( Figure 1). The district is located between 142,000 and 174,000 UTM E and 99,300-1,029,000 UTM N.
Gimbi district has a Weina Dega climate condition. The mean minimum and maximum annual temperature of the area was between 10°C and 30°C. The mean annual rainfall is 1400-1800 ml per year (Abera et al., 2015). In the Gimbi district, June, July, and August are the summer months similar to most of the country's annual seasons. It lies at an altitudinal range of 1200-2222 m.

Figure 1. Map of the study area
Every location in Ethiopia is categorized as lowland, midland, or highland depending on its altitude: below 1500 m asl., between 1500 and 2300 m asl., and beyond 2300 m asl. In this categorization, only altitudinal gradients are taken into consideration. Gimbi district only has lowland and midland altitudinal gradients because it is located between 1200 and 2222 meters above sea level.

Sampling techniques
First, we identified classifications of coffee landscapes based on the farmer's perceptions and experts. They classified the landscapes into three altitudinal gradients namely lower land (1451-1650 m asl.), mid-land (1651-1850 m asl.) and upper land (1851-2050 m asl.). Then, we adopted this classification scheme for our study. We identified three altitudinal gradients of coffee-based agroforestry systems, with an area of 3.2 ha, for each gradient. Across the three-altitudinal gradients, an area of 9.6 ha was identified using Google Earth. In each altitudinal classification, a 1 km long transect was established. Generally, six transects were established across the whole altitudinal gradients over six kebeles 1 in the district.

Sampling design
Data was collected from February 2020 to May 2020. Tree species diversity and biomass carbon stock assessment were carried out using areas of 40 m × 40 m sample plots, along the altitudinal gradients classification. Square plots were selected as they tend to include more of the within-plot heterogeneity and; thus be more representative than circular plots of the same area. A total of 60 sample plots were used for tree parameter measurements. In each plot, tree species diversity, tree height, and tree diameter at breast height (DBH) were recorded following the transect lines.
The coffee farm chosen for the purpose of the study has been producing for a long time, so the trees on the farm were large in height and diameter. During the inventory, the height of most of the trees was >2 m and their thickness was >5 cm. However, a limited number of trees with a height of <2 m and a diameter of <5 cm observed occasionally. Since the number of trees (h < 2 m and d < 5 cm) was small, they have not been included in the estimation of biomass carbon stocks.
Hence, trees with a height ≥2 m and ≥5 cm in diameter at breast height (DBH) were measured. Trees height and diameter were measured using a vertex, diameter tape, and caliper respectively. Branched trees at 1.3 m were measured at the smallest point below 1.3 m, where the stem assumes a nearly cylindrical shape. A woody plant with multiple stems or forked below 1.3 m height was treated as a single individual. Plant identification was done in the field using their local name and useful trees and shrubs for Ethiopia (Bekele-Tesemma, 2007). Nomenclature follows the Flora of Ethiopia and Eritrea (Vol 1-8). To quantify the tree diversity of the coffee-based agroforestry systems, the number and names of tree species (both seedling and sapling) were recorded.

Tree species diversity
For this study, we have used different diversity indices. The weaver indices of the sampled coffeebased agroforestry systems were calculated by using Shannon diversity (Shannon & Weaver, 1963) which is a very widely used index for comparing diversity between various habitats (Okpiliya, 2012). The evenness of a population was estimated following the Shannon-Wiener index. Additionally, Simpson's diversity index was used (Simpson, 1949). Simpson's diversity index is the most sensitive to changes in more abundant species and hence places more weight on the most abundant species in the community. The Simpson's diversity index is derived from probability theory and it is the probability of picking two organisms at random which are of different species (Magurran, 1988).
The Sorensen similarity coefficient was used to estimate the similarity/dissimilarity of species between altitudinal gradients (Sorensen, 1948). The Sorensen similarity coefficient was widely deployed as it gives more weight to the species that are common to the samples rather than to those that only occur in either sample (Gotelli & Chao, 2013).
The IVI of a species was used to express the relative ecological significance of the species in the forest ecosystem. It was calculated by summing up the relative dominance, relative density, and relative frequency of the species. It indicates the significance of species in the system (Mengistu & Asfaw, 2016).

Tree carbon estimation
In this study, the non-destructive allometric equation which is described by Chave et al. (2014) was used to estimate both the aboveground and belowground biomass (ton). This method is applied to the global forest type. It helps calculate the aboveground and belowground biomass contained within each tree that is given as a function of height (m), DBH (cm), and wood density (g/cm −3 ).
The aboveground biomass of trees and shrubs in the coffee-based agroforestry was calculated using the formula in Eq. (1).
AGB ¼ 0:0673XðρD 2 HÞ 0:976 ; (1) Where, AGB = above ground biomass Aboveground live carbon (AGC) was estimated at 50% of AGB, while belowground biomass carbon estimated by multiplying the aboveground biomass by 0.26 as conversion factor (Chave et al., 2014;Manickam et al., 2014). Wood-specific gravity was obtained at the species level from the Global Wood Density Database (Chave et al., 2009). For example, one of the studied trees species (Albizia gummifera) from plot1 in the coffee-based agroforestry system has a height of 28 m and DBH of 65 cm. The wood density of this tree species is 0.580. Therefore, the aboveground biomass of the tree was calculated as follow.
Where, T = ton C = Carbon

AGB = Aboveground biomass
In order to estimate the total carbon sequestered in the agroforestry system, it was necessary to calculate the belowground biomass carbon stock. Therefore, the belowground biomass carbon (BGBC), was calculated by multiplying AGB taking 0.26 as the root-to-shoot ratio used by Sudha et al. (2007). And subsequently, the carbon stock of belowground biomass was calculated by using the equation developed by Manickam et al. (2014).
Furthermore, the formula developed by several scholars is what is used to transform biomass into carbon stock. According to the premise that living biomass (tissues) contains 47% carbon, the carbon stock for both trees and shrubs/bushes was calculated to be 0.47 times the amount of woody biomass (Chave et al., 2005). In other words, the biomass has a carbon content of roughly 47% by dry weight and the carbon stock in the biomass was estimated using the formula in Eq. (5).

Statistical analysis
The data which was collected from the field inventory was organized and recorded in Microsoft excel 2010 datasheet. The altitudinal gradients were the independent variables while density, basal area, species richness, species evenness, and species diversities were considered as dependent variables. Vegetation data were analyzed using univariate analysis. Variables were compared using one-way analysis of variance (ANOVA) following the linear model (GLM) procedure at P < 0.05 with the help of the R software for Window versions 3.6.1 with statistical significance difference were observed (P < 0.05), Tukey's HSD test were used to separate the means. To investigate the effect of tree diversity on carbon stock across altitudinal gradients, repeated analyses of Variance was performed and the relationship between tree diversity and tree and coffee shrubs biomass were analyzed.

Floristic composition
A total of 36 tree species were recorded from coffee-based agroforestry systems (CBAS). Out of this, 27 species were recorded from the 1451-1650 m altitudinal range and 29 were from the 1651-1850 m altitude, while 17 species were from the 1851-2050 m altitude, respectively (Table  A1, A2, A3). The entire recorded species were belonging to 34 genus and 23 families (  Figure 2). The remaining 15 families altogether accounted for 41.7% of the total species composition each family represented by 1 species (Table 2). Overall, 3,231 tree individuals and 36,424 coffee stems were counted (Table 3).
Species richness, Evenness (E), Shannon-Wiener Diversity Index (H'), and Simpson Diversity (1-D) were significantly different along an altitudinal gradient ( Table 2). The mean species richness,   The standard error (±SE) mean is shown in parenthesis. Note; one-way ANOVA was conducted to see the mean difference between groups and followed by a post hoc test for multiple comparisons. A similar letter indicates have no significant difference but a different letter shows the statistically significant difference between altitudinal gradient at p < 0.05 gradient ( Table 2) Likewise tree species density, the basal area has also shown a significant difference (F = 5.349, P < 0.05) along the altitudinal gradient. Lower altitude showed a higher mean basal area (3.67 m 2 per sample plot) of tree species and the lowest mean basal area was found at 1851-2050 m altitude (1.82 m 2 per sample plot; Table 3).

Carbon stock
The pattern of aboveground and belowground carbon stocks distribution varies across altitudinal gradients. The differences in tree carbon stock, both aboveground and belowground, showed a slight increase from 1451-1650 m altitude to 1651-1850 m while it was decreasing from 1651-1850 m to higher 1851-2050 m altitude in the coffee-based agroforestry system ( Table 4).
The pattern of aboveground and belowground biomass carbon stock in coffee shrub across altitudinal gradients was significantly different along the altitudinal gradients. It was monotonically decreasing as altitude increased. The maximum aboveground and belowground biomass carbon recorded in coffee shrubs were 4.73 t ha −1 and the minimum was 1.11 t ha −1 respectively ( Table 4).
The highest BGC for tree species was recorded in the 1651-1850 m altitudinal class (t ha −1 ) followed by 1451-1650 m altitudinal range (9.01 t ha −1 ). Similarly, the variation of carbon stock in coffee shrubs along the altitudinal gradient was statistically significant (F = 6.0408, P < 0.001). Overall, the highest carbon stock was observed in the 1451-1650 m altitudinal range while the lowest was found in the 1851-2050 m altitudinal scale (Figure 2).

Total carbon stocks and tree diversity relationships
The carbon stock in the study area was influenced by the variation in tree species diversity, abundance, and dominance of trees across the altitudinal gradients in the coffee-based agroforestry system (Table 4; Figure 3). The largest carbon stock was contributed by the tree species, which accounted for 93.3% % of the system for the three altitudinal gradients while the remaining was contributed by coffee shrubs ( Table 4). The mean tree carbon stock was calculated to be 38.68 t ha −1 , and the average coffee shrub carbon stock was estimated to be 2.73 t ha −1 . The overall mean carbon stock stored in the system was estimated to be 41.4 t ha −1 (Table 4).

Floristic composition
Floristic composition and variety were found to be highest in the 1651-1850 m height, followed by the 1451-1650 m altitude, and lower in the 1851-2050 m altitude. Desalegn and Beierkuhnlein (2010) found a similar trend of species richness and pattern. They discovered that woody species richness follows a trend that begins with an increase and peaks about 2000 m, then begins to fall and level off at elevations over 2500 m. In addition, Gebrewahid and Abrehe (2019) stated that species richness had a hump-shaped distribution, with the highest species richness in the intermediate elevation and the lowest species richness at the upper elevation. Altitudes and slopes,   according to Eilu and Obua (2005), influence tree species richness and dispersion behavior. As reported by Sharma et al. (2017), tree species richness was found to be highest at a lower altitude, followed by middle altitude, and lowest at a higher altitude. They also reported that tree species richness (including genera and family richness) shows a similar pattern, reaching a maximum at intermediate elevations in the subtropical Andes of Argentina.

Table 4. Mean (± standard error; n = 20) aboveground and belowground biomass (t ha−1 ) for each of the three altitudinal gradients and results of one-way ANOVA (at α = 0.05, significant differences between the altitudinal gradient of the biomass carbon stocks)
Several studies have shown that species richness patterns can show a monotonic reduction (Gebrehiwot et al., 2019;Sharma et al., 2017). Ali and Yan (2017) and Gebrehiwot et al. (2019) have found a monotonic increase along altitudinal gradient. In contrast, some researchers have discovered that species richness is lowest at mid-elevations, resulting in an inverted hump-shape (Nanda et al., 2018), while others have discovered no apparent association between species richness and elevation. According to Toledo-Garibaldi and Williams-Linera (2014), species, genus, and family richness decreased in a unimodal fashion along the entire elevation gradient, while diversity decreased monotonically.
Several researchers conclude that there is no general pattern along elevation gradients; however, particular patterns, such as the unimodal hump-shaped and monotonic reduction, do exist. Due to a multitude of rules driving the general ecology of forest ecosystems (Lawton, 1999) and the existence of a wide variety of strategies and ways developed by species to adapt to environmental stresses (Grime, 1977). However, there is no general conclusion on the correlation between species distribution and altitude, the fact that tree species show variation along altitudinal gradients. Pausas and Austin (2001) affirmed that the distribution of species richness over any wide area is likely to be influenced by two or more environmental factors.
In the 1851-2050 m altitudinal range, Shannon-index Weiner's values ranged from 0.88 to 1.75 and dropped. Our findings are consistent with those found in other temperate woods (Sharma et al., 2009). Similarly, the important value indexes differed throughout elevation gradients, and different tree species predominated at various altitudes. In this study, Albizia gummifera was the major tree species in the coffee-based agroforestry systems at altitudes of 1451-1650 m and 1651-1850 m, respectively, whereas Croton macrostachyus was the dominating tree species at elevations of 1851-2050 m. According to Curtis and Mcintosh (1950), the highest important value index obtained for a given woody species indicates that it is the dominating species among other plant species in the same area.
According to Nemera et al. (2018), a species' ecological dominance may be linked to pathogen resistance or competition, adaptability to growth in the shade, least preferred by browsing animals, pollinator attraction, and diversity of seed predators that facilitated seed dispersal within the existing environmental conditions. Environmental factors and selective human use of the available species, on the other hand, may be due to the woody species with the least important value.
Tree dominance varies with variation in the coffee-based agroforestry systems across elevation. Simpson's index value was 0.14 in the 1451-1650 m and 1651-1850 m gradients, which gradually increased and reached finally up to 0.41 for the tree species in the 1851-2050 m gradients. This was due to the tolerance-based dominance of only a few tree species under severe environmental conditions. The high importance value index (IVI) of species between 1651 and 1850 m altitudes indicated their dominance and ecological success, good power of regeneration, and ecological amplitude. Similar observations were also recorded by Sharma et al. (2017) and .

Tree diversity and carbon stocks along altitudinal gradients
The mean above and below ground biomass carbon stocks of both trees and coffee shrubs decreased with increasing altitude, which could be attributable to the absence of large trees with maximum DBH, a decrease in tree species diversity, and a decrease in the number of coffee shrubs in the upper altitude. Because tree species richness decreased as altitude increased from 1851 to 2050, relatively small trees (low DBH) with low diversity became dominant. This could be the result of anthropogenic influences in the study area. The farmers did raise animals and grow other crops, but they mostly grew coffee. The population density fluctuated along the altitudinal gradients' ranges as well. For instance, the population of the upper land increased along with the number of cattle and crops produced which displaced coffee production and transformed the land into open agriculture, according to agricultural office report.
On the other hand, in the lower altitudes, farmers cut trees for charcoal making and wood selling. Furthermore, as altitude increased, the weather variables of the area varied with a pattern of decreasing at the upper altitude. Because of these factors, lower and higher tree diversity were observed at lower and upper altitudes, respectively. The study's findings were similar to those of Esubalew et al. (2019) and Tesfaye and Seifu (2016). According to them, the variations in the parameters could be due to a variety of factors that change with altitude. The elements that caused the variation could have been soil humidity or temperature, as well as a study site's dissimilarity. Similar research has found that a number of factors influence the variability of tree variety, which affects both above and belowground carbon reserves in a forest ecosystem. Because of good species diversity, environmental conditions, and soil characteristics in the mid-altitudinal gradient, Eshetu and Hailu (2020) reported that the pattern of total carbon stock in their study had a humped shape with altitudinal gradient and formed a peak carbon stock at the mid-altitude. Yohannes et al. (2015) discovered a declining trend in soil organic carbon stock as elevation increased, which they attribute to canopy cover, litter biomass accumulation, and species diversity. Another reason, as Gedefaw et al. (2014) point out, is that climatic conditions can influence forest carbon stock along an elevation gradient. Farm size, management, socioeconomic needs, species diversity, tree age, local climate, and tree spacing among agroforestry systems could all have a role in carbon stock variation (Kumar, 2011). A continual removal of falling litter, deadwood, and twigs, as well as unlawful cutting for charcoal production, construction, agriculture, and livestock grazing, were also discovered to have an impact on the forest carbon stock balance (Gedefaw et al., 2014). The highest amount of carbon was stored in coffee shrubs at 1451-1650 m, followed by 1651-1850 m, and the lowest at 1851-2050 m. That might be because the upper land forest, which located in the 1851-2050 m altitudinal gradient, is dominated by natural forest, while the 1651-1850 m altitude is mixed, and the 1451-1650 m altitude is dominated by planting because of farmers' shade tree preference intervention.
This research provides valuable information on tree diversity, coffee shrub density, biomass, and carbon stock in coffee-based agroforestry ecosystems over an altitudinal gradient in the study area. The impact of human intervention on tree species variety and coffee shrub density throughout an altitudinal gradient is also highlighted in our research. In comparison to the 1451-1650 m and 1651-1850 m systems, the agroforestry system in the 1851-2050 m altitudinal range is frequently subjected to more intense human disturbance and has a lower density of huge trees, according to local experts. The agroforestry systems in the altitudinal ranges of 1451-1650 m and 1651-1850 m sustained more large trees due to long-term management of larger trees by coffee producers. The 1851-2050 m altitude, which mainly comprises natural forest, was found to contribute low accumulations of carbon stocks as compared to the 1451-1650 m and 1651-1850 m altitudes. The more coffee shrubs, large tree diversity, and high carbon storage potential of 1451-1650 m and 1651-1850 m altitudes show that the two sites are more productive and have more sustainable ecosystems compared to the 1851-2050 m altitude, which has lower coffee shrubs density and small size and less diversity tree with less carbon stock potential.
Thus, comparing these species diversity and carbon storages phenomena, coffee-based agroforestry systems can be seen as an agroecosystem function that is directly interrelated to biodiversity conservation and carbon storage. However, the effect of species diversity on biomass carbon storage depends on management practices, species, and site factors.

Conclusions
This study highlights coffee-based agroforestry systems are important for conserving Ethiopian tree species that are either indigenous or highly endangered. In comparison, the system has a high proportion of native woody species equal to some unmanaged natural forests. Farmers typically establish or keep nearly multipurpose indigenous trees in the coffee-based agroforestry system. Our finding revealed that the system considerably improves native tree conservation.
Compared to the coffee shrubs, shade trees significantly stored large carbon stocks across the altitudinal gradients. Higher carbon stock and large trees with diverse species recorded in the lower and mid-altitudinal gradients indicate that plant density and diversity can actively influence total carbon stock in agroforestry systems than elevation gradients. This shows that the system is vital for carbon sequestration and help in alleviating climate change issues.
The system contributes to the augmentation of productivity of the area by providing food and assisting the community's livelihoods, protecting soil and watersheds. Hence, coffee-based agroforestry systems were found to be vital for biodiversity conservation, carbon sequestration, and ecosystem services that can lead to sustainable ecosystem management.