Study area description
The study was conducted in indigenous agroforestry systems of Gedeo zone, Dilla Zuria District, Southern Ethiopia, lies approximately 396 km south of the Addis Ababa (Fig. 1). The study district, Dilla Zuria, is situated between 6º15’05’’ N- 6º26’35’’N latitude and 38º 15’55’’E and 38º 24’02’’E longitudes (Fig. 1), covering an area of 120 km2, with different land-use types such as agricultural land and AF accounting for 95%, while grassland, wetland, plantations, and others covered the remaining 5% (Mebrate 2007). The altitudinal range of the district ranges between 1350 m to 2550 m with a slope of 39.4–51.5% (Yirefu and Wendawek 2016). The district is characterized by three agro-ecological zones (AEZ), namely highland (Dega) constitutes 23%, Midland (Woyna Dega) cover 70%, and Lowland (Qolla) embraced 7% of the total area. The total population of Dilla Zuria was 137,715 (M = 71262, F = 66453).
Study design
Sampling method and sample size determination
A multi-stage sampling approach was used to gather data on vegetation diversity and management practices. The Dilla Zuria District in Gedeo, Southern Ethiopia, known for its agroforestry systems (Kassa et al. 2015; Tadesse 2002), was chosen first. The area was divided into three elevation categories: lower (1500-1750m), middle (1650-1990m), and upper (2300-2550m). Kebeles within each elevation were stratified by agroforestry types. Three Kebeles—Michile Girisa (higher elevation), Michile Sisota (middle elevation), and Chichu (lower elevation)—were randomly selected. Finally, 286 households were randomly chosen from a master list of 1006 households using proportional allocation (Yamane T 1967) across three agroforestry practices (Table-1).
Table 1
Name of Study Sites, agroforestry practices, populations from which sample were drawn and Sample Households of the study area in the Southern Ethiopia, Gedeo
Site name (Kebeles)
|
Agroforestry practices
|
Total Household
|
sample
|
Elevation gradients
|
Michile-Sitota
|
Coffee-Enset-tree based
|
266
|
76
|
1650–1990
|
Michile-Girisa
|
Enset-based
|
308
|
87
|
2300–2550
|
Chichu
|
Fruit-coffee-tree-based
|
432
|
123
|
1500–1750
|
Total
|
|
1006
|
286
|
|
Methods of data collection
Sample Layout and plant inventory
A vegetation inventory was conducted on 286 randomly selected farms, with one plot per farm. Farms were divided into 10 by 10 plots using ocular estimation, and one plot was randomly chosen (Negash et al., 2011). If a household had multiple farms of the same type, only one was selected1. Quadrat sizes and sampling methods followed recommended practices (Hafte et al., 2024). Due to the small size of farms and considerations of cost and time, the quadrat size was limited to 100 m². Ethiopia’s agro-ecology is divided into three zones: Dega (highland), Woina-Dega (mid-highland), and Kolla (lowland) (Bekele, 2007). The Dega Zone includes highlands over 2,300 meters, while the Woina-Dega Zone covers areas between 1,500 and 2,300 meters4. Locally, this classification is divided into lower midland (1,500-1,750 meters) and upper midland (1,650-1,900 meters). Thus, the study area was stratified into upper, middle, and lower midland agroforestry farms.
All woody species in the plots, whether single-stemmed (Negash et al. 2011) or multi-stemmed (Snowdon et al. 2002), with a DBH of at least 2.5 cm and a height of at least 1.5 m, were measured for their height and DBH. Their uses and growth habits were identified and recorded to the species level with the assistance of local informants.
To measure the coffee shrubs in each plot, the DBH was measured at a height of 40 cm (Negash et al. 2011). For Enset, the basal diameter of the pseudostem of all Enset found in the quadrates was measured at a height of 10 cm (Tesfay et al. 2024). The scientific and vernacular names of plant species, local names, their family, genus, and other relevant information were recorded on the data sheet. For identifying the selectively removed trees number and species in the plots involves a combination of direct observation, counting, identifying species from the stamps left behind, and discussions with farmers.
Socio-economic data collection
A mixed-methods approach was used to collect socioeconomic data through face-to-face and key informant interviews, and focus group discussions (FGDs). Structured questionnaires were administered to 286 rural households, gathering information on socioeconomic and demographic factors, management practices, species use, wealth status, and agroforestry product utilization. Key informant interviews with community leaders, development agents, and knowledgeable farmers provided insights into community wealth status (Crabtree & Miller, 1992). Snowball sampling broadened the sample size. FGDs, involving six to twelve participants each, were conducted in three agroecological settings, focusing on plant diversity, management practices, and socioeconomic factors influencing agroforestry changes. Researchers also reviewed existing research to understand the historical and future dynamics of Gedeo agroforestry.
Data analysis
Floristic composition and diversity
All recorded data of each agroforestry practices were pooled and the total number of species and individuals were tallied. Using the pooled data, number of individual plants, number of species and their scientific names, Name of families, life form per species (trees, shrubs and herbs), identity (Native/Non-Native) and major uses of each species across agroforestry systems were recorded in the plots using supplementary field guide (Bekele-Tesemma 2007), and help of knowledgably key informants (Appendix-I).
The mean diameter, density, and total basal area, including number of individuals per species, ecologically and economically preferred species by farming households were calculated for each species found in sample plots across three agroforestry systems (Appendix-).
Perennial plant species diversity (species richness and Shannon–Wiener diversity, and Simpson Indices) and abundances of individuals.
Shannon Diversity Index (H’)
Shannon diversity index, Simpson's index, and evenness were calculated for all sampled plots to evaluate the richness and diversity of shade trees in the two systems. The Shannon index H’ (Shannon and Weaver 1949) is a measure of the number of species S and their even distribution according to the proportion of species pi and can be calculated as:
$$\:{\text{H}}^{{\prime\:}}=-\sum\:_{\text{i}}^{\text{S}}\text{p}\text{i}\text{*}\left(\text{L}\text{N}\text{p}\text{i}\right)$$
The Simpson index (D)
The Simpson index D is measure of diversity, which takes into account the number of species S and the relative abundance of each species with values starting from one. Higher values obtained with the Simpson index indicate greater diversity. It is calculated as (Simpson 1949)
$$\:\text{D}=\:\:\frac{1}{{\sum\:}_{\text{i}}^{\text{S}}{\text{P}\text{i}}^{2}}$$
The number of species in each plot, known as species richness, was determined by counting the different species present in the area.
Floristic structure
The floristic structure of all perennial plant species was examined in terms of density, frequency, basal area, important value index, and distribution by diameter class. Microsoft Excel was used to conduct the analysis.
Basal area
In order to compare the basal area and density between different DBH classes, we grouped the life form of trees, shrubs and herbs based on the mean DBH classes (0–15, 16–30, 31–45, 46–60, 61–70) of stems for each of agroforestry system.
We then calculated the mean of the basal area and density for each of the DBH classes for trees, shrubs and herbs.
$$\:\text{B}\text{a}\text{s}\text{a}\text{l}\:\text{a}\text{r}\text{e}\text{a}\:\left(\text{B}\text{A}\right)\:\left(\frac{\text{m}2}{\text{h}\text{a}}\right)=\frac{\text{s}\text{u}\text{m}\:\text{o}\text{f}\:\text{c}\text{r}\text{o}\text{s}\text{s}\:\text{s}\text{e}\text{c}\text{t}\text{i}\text{o}\text{n}\text{a}\text{l}\:\text{a}\text{r}\text{e}\text{a}\:\text{f}\text{o}\text{r}\:\text{a}\text{l}\text{l}\:\text{t}\text{r}\text{e}\text{e}\text{s}\:\text{i}\text{n}\:\text{p}\text{l}\text{o}\text{t}\text{s}}{\text{s}\text{t}\text{u}\text{d}\text{y}\:\text{p}\text{l}\text{o}\text{t}\:\text{a}\text{r}\text{e}\text{a}\:\text{i}\text{n}\:\left(\text{h}\text{a}\right)}$$
\(\:\text{C}\text{r}\text{o}\text{s}\text{s}\:\text{s}\text{e}\text{c}\text{t}\text{i}\text{o}\text{n}\text{a}\text{l}\:\text{a}\text{r}\text{e}\text{a}\:\text{o}\text{f}\:\text{a}\:\text{t}\text{r}\text{e}\text{e}={\pi\:}\:\left(\frac{\text{D}\text{B}\text{H}\:\left(\text{c}\text{m}\right)}{200}\right)2\) = 0.00007854* (DBH cm)2
Where
-
Ha stands for hectare; and
-
DBH for diameter at breast height in cm
-
Area in ha = area of plots in ha = 0.01ha (100m2)
Density
Density was determined by adding up all the stems found in each sample plot and then converting this total into a measurement per hectare (Mengistu and Asfaw 2016)
$$\:\text{D}\text{e}\text{n}\text{s}\text{i}\text{t}\text{y}=\frac{\text{t}\text{o}\text{t}\text{a}\text{l}\:\text{n}\text{u}\text{m}\text{b}\text{e}\text{r}\:\text{o}\text{f}\:\text{a}\text{l}\text{l}\:\text{t}\text{r}\text{e}\text{e}\text{s}}{\text{s}\text{a}\text{m}\text{p}\text{l}\text{e}\:\text{s}\text{i}\text{z}\text{e}\:\text{i}\text{n}\:\text{h}\text{e}\text{c}\text{t}\text{a}\text{r}\text{e}\:}$$
Important Value Index (IVI)
It was calculated by summing up the relative dominance, relative Frequency and relative abundance of the species (Mueller-Dombois and Ellenberg 1974). It shows the degree of a certain plant species' dominance, occurrence, and abundance in comparison to other nearby specie (Whittaker 1993). It is calculated as follows:
$$\:\text{R}\text{D}\text{o}=\frac{\text{B}\text{a}\text{s}\text{a}\text{l}\:\text{a}\text{r}\text{e}\text{a}\:\text{o}\text{f}\:\:\text{s}\text{p}\text{e}\text{c}\text{i}\text{e}\text{s}\:\text{A}}{\text{A}\text{r}\text{e}\text{a}\:\text{o}\text{f}\:\text{s}\text{a}\text{m}\text{p}\text{l}\text{e}\text{d}\:}\text{*}100\text{%}$$
RF= \(\:\frac{\text{T}\text{o}\text{t}\text{a}\text{l}\:\text{n}\text{u}\text{m}\text{b}\text{e}\text{r}\:\text{o}\text{f}\:\text{i}\text{n}\text{d}\text{i}\text{v}\text{i}\text{d}\text{u}\text{a}\text{l}\text{s}\:\text{i}\text{n}\:\text{a}\text{l}\text{l}\:\text{q}\text{u}\text{a}\text{d}\text{r}\text{a}\text{t}\text{e}\text{s}}{\text{T}\text{o}\text{t}\text{a}\text{l}\:\text{n}\text{u}\text{m}\text{b}\text{e}\text{r}\:\text{o}\text{f}\:\text{q}\text{u}\text{a}\text{d}\text{r}\text{a}\text{t}\text{s}\:\text{i}\text{n}\:\text{w}\text{h}\text{i}\text{c}\text{h}\:\text{t}\text{h}\text{e}\:\text{s}\text{p}\text{e}\text{c}\text{i}\text{e}\text{s}\:\text{o}\text{c}\text{c}\text{u}\text{r}\text{r}\text{e}\text{d}}\:\)
RA=\(\:\frac{\text{a}\text{b}\text{u}\text{n}\text{d}\text{a}\text{n}\text{c}\text{e}\:\text{f}\text{o}\:\text{o}\text{n}\text{e}\:\text{s}\text{p}\text{e}\text{c}\text{i}\text{e}\text{s}\:}{\text{t}\text{o}\text{t}\text{a}\text{l}\:\text{a}\text{b}\text{u}\text{n}\text{d}\text{a}\text{n}\text{c}\text{e}\:\text{o}\text{f}\:\text{a}\text{l}\text{l}\:\text{s}\text{p}\text{c}\text{e}\text{i}\text{s}\:}\text{x}100\)
$$\:\text{I}\text{V}\text{I}=\text{R}\text{D}0+\text{R}\text{F}+\text{R}\text{A}$$
Where, IVI is importance value index, RD is relative Dominance, RF is relative frequency and RA is relative abundance of a species
Data Analysis of Effects of agroforestry systems, management practices and socio-environmental factors on species Diversity and richness
The research analyzed of the impacts of AFSs, management practices, abundance of plants (specifically maintained for economic purpose: dividing the number of economically important perennial plants by the total number of perennial plants in each site (Rodrigues et al. 2018), removal of canopy tree species (recorded the number of canopy trees removed per plot across AFSs with visual observation at farm level and interview with household heads), Altitude (m.a.s.l), Age, Land size, Educational status, and Wealth status of farm households on species diversity and abundance of woody and non-woody plant species in the study area. The ratio of economically oriented perennial plants in AFSs by categorizing associated plants into Economic (plants maintained for income and domestic needs by AF farming households) and for ecological purpose (Asigbaase et al. 2019) were analyzed .
Statistical analysis
Descriptive statistics are employed to detail the demographic and socioeconomic characteristics of the household sample, including measures like the average, percentage, standard deviation, and frequency. One-way ANOVA tests were conducted to compare the statistical differences between different types of AFPs, as well as between different life forms (Trees, shrubs, and herbs) for normally distributed data. Prior to each test all data was tested for normality by running and plotting a Shapiro-Wilk normality test with a significance level of p < 0.05, and a Levene’s test with a significance level of p < 0.05 was performed to test for homogeneity of variance. A post hoc LSD test was carried out when significant differences were found among the AFPs. Where the data did not follow a normal distribution, a Kruskal-Wallis one-way ANOVA test was used.
Further, a separate multiple linear regression model was built to analysis the impacts of agroforestry systems, management practices and social and environments factors (independent variables). This model is appropriate for investigating the impact of multiple independent variables on a single dependent variable (Montgomery et al., 2012). The model described the effects of the types AFSs, management activities carried out in the study area, abundance of economically managed species, and other socioeconomic factors (such as land size, wealth status, educational status, and age) and environmental factors (altitude: Different altitudes offer varying climatic conditions, impacting species composition and growth on species richness, diversity indices (Shannon and Simpson) and abundance of woody and non-woody species.