Species ethnobotanical values rather than regional species pool determine plant diversity in agroforestry systems.

The conversion of natural systems into farms and agroecosystems is the main cause of biodiversity loss. In human-dominated landscapes, understanding the interactions between agroforestry systems and adjacent natural vegetation is fundamental to developing sustainable agricultural systems. Species can move between these two systems with natural systems providing the regional pool of species that shape the agricultural values and conservation value of the agroforestry systems. We investigated the influence of neighboring natural habitats on traditional agroforestry systems in the buffer zone of Pendjari Biosphere Reserve in Benin to understand the contribution of regional processes on the quality of agroforestry systems. We expected that agroforestry parklands adjacent to natural vegetation with high species diversity will also have higher plant species diversity. We found no similarity in plant species composition between agroforestry systems and adjacent natural habitats. A small proportion of species in adjacent natural habitats were found in agroforestry systems. The proportion of shared species was not significantly influenced by plant diversity in adjacent natural habitats or the distance from the agroforestry systems to the natural adjacent habitat. However, plant diversity in agroforestry systems was strongly associated with site ethnobotanical values indicating that farmers act as a supplemental but severe environmental filter of the regional species pool. Our study suggests that promoting the plantation of plants with high ethnobotanical use-value is a potentially viable strategy for sustainable agriculture and ecological restoration in Biosphere reserves.

Introduction

Understanding how multiple factors govern species diversity in ecosystems and the practical consequences for species diversity is a core question in ecology and conservation biology[1,2]. The spatial patterns of species diversity is driven by a wide range of ecological processes including facilitation[3], competition[4], disturbance[5-7], productivity[6,8] and productivity–disturbance interactions[9]. Communities tend to be highly diversified when ecological processes provide more favorable conditions for species coexistence[9-11]. The environmental filtering hypothesis suggests that local species assemblages are the result of the filtering of regional species pools by constraints posed by abiotic and biotic factors[12-14]. Empirical evidence of environmental filtering is rare[15]. Local species assemblages may result from the single or interactive effects of several drivers including landscape or ecosystem modifications by human activities[16,17].

In human modified landscapes, understanding the interactions between agricultural systems and adjacent natural habitats is fundamental to developing sustainable agricultural systems[18,19]. Environmental filtering can shape diversity in agricultural systems with consequences on the composition of plants that can migrate from nearby natural habitat to act as weed to crops. In addition, plant diversity in adjacent natural habitat can shape arthropod community diversity[20-22] which can affect crop herbivory rate and ultimately crop yield[23]. The strength of environmental filtering in shaping species assemblages in agroforestry systems will directly depends on the dispersal abilities of species from the regional pools[24-27] but also on the distance from these pools. Consistent with the island biogeography theory[28,29], the closer the agroforestry systems are to the regional species pools, the more likely are species going to disperse into agroforestry systems and the higher diversity and shared species richness are expected.

Biodiversity conservation in biosphere reserves is increasingly centered around the effort to minimize the impact of chronic anthropogenic disturbance given the persistent encroachment of croplands into these protected areas[30]. Biosphere reserves are structured into three zones with different functions[31]: (i) a core zone, which is strictly protected and where only research and ecological monitoring activities are allowed, (ii) a buffer zone used for low impact tourism, forestry, and agriculture in line with overall conservation objectives and (iii) the transition zone, often the largest, where a variety of human activities are allowed[32,33]. In the transition zone of most African biosphere reserves, there is an increasing transformation of the natural vegetation into islands of wilderness surrounded by farms and agroforestry systems. This is due to the rapidly increasing world population, land shortage in rural areas and pressure from agricultural intensification[34,35]. Investigating the factors that shape diversity in the traditional agroforestry systems in the transition zone of Biosphere reserves is central to maintaining the double role of protected areas to safeguarding long-term ecological sustainability and satisfying basic human needs[36].

We used the theoretical framework of the environmental filtering and island biogeography theories to examine the effects of various socio-cultural and agricultural (ethnobotanical value, agricultural practices) and proximity to natural habitat (distance from adjacent natural vegetation) on species diversity in traditional agroforestry systems in the Pendjari biosphere reserve in Benin (west Africa, Fig. 1;[37]). We analyzed the changes in species diversity at different levels of agricultural land use intensity. Consistent with the island biogeography theory, we hypothesized that the distance from agroforestry systems to adjacent species pools of natural vegetation will predict the number of species in agroforestry systems and the proportion of shared species between the two plant communities. The closer the agroforestry systems are to their surrounding natural vegetation, the greater number of plant species the two communities will share and the higher the species richness will be in the agroforestry systems. However, we hypothesized that the regional species pool will be filtered by abiotic/biotic factors but also the influence of local farmers who decide what species to clear from or add to their farms based on their ethnobotanical value. We hypothesized that species ethnobotanical values along with the proximity of agroforestry to natural vegetation will predict local species assemblages in agroforestry parklands.

Figure 1

Map of the hunting zone of the Pendjari Biosphere Reserve in Northern Benin showing study sites in agroforestry parklands that are established in the transition zone where sustainable agricultural activities are allowed. The transition zone is composed of two narrow strips of land on each side of the hunting zone (grey). On the Tanguieta-Batia road axis, populations are constrained by the Atacora Chain of Mountains which boarded the eastern side of the hunting zone limiting agricultural expansion. The map was created using ArcGIS 10.5.1 (ESRI, https://www.esri.com).

Results

Similarity in diversity and species composition between agroforestry systems and adjacent natural vegetation

The rank-frequency diagram for the agroforestry systems was different from that of the adjacent natural habitat (Fig. 2). We estimated a total of 8 species in the agroforestry systems against 31 species in adjacent natural habitats indicating a very low diversity in agroforestry systems. In addition, individuals were more fairly distributed among species in the adjacent natural habitat than in the agroforestry systems. However, we found similarities in relative abundance of a few species. For example, Vitellaria paradoxa was the most dominant species in both ecosystems followed by Parkia biglobosa. However, species composition besides these dominant species were different between the two ecosystems. For example, in addition to V. paradoxa and P. biglobosa, the agroforestry systems were represented by only singleton species (5 species) whereas adjacent natural habitat included singleton (12 species) and doubleton (7 species). In addition, regular species found in the adjacent natural habitat included Combretum fragans, Terminalia avicennioides, Burkea africana whereas only Lannea acida was represented in agroforestry systems.

Figure 2

Rank-frequency diagrams of adjacent natural formations and traditional agroforestry systems of Pendjari Biosphere Reserve.

The non-metric multidimensional scaling (NMDS) identified three groups of plant communities (Fig. 3). The first group (A) was composed of two plots that were all established in the natural habitats especially in savannah and rocky outcrop habitat. This first group was highly dissimilar to the others two groups which included plots from moist savannah with unique species such as Pterocarpus erinaceus, Piliostigma thonningii, Terminalia macroptera, Terminalia avicennioides. These two plant communities (B and C) were more similar and included a mixture of agroforestry systems plots (plot #1, 3, 4,5, 7, and 8) and those of adjacent natural vegetation. Typical agroforestry species such as Parkia biglobosa, Vitellaria paradoxa and Anacardium occidentale were shared between these two communities.

Figure 3

Nonmetric dimensional scaling analysis showing the dissimilarity between agroforestry systems and their surrounding natural vegetation. NMDS identified three plant communities (A–C). Numbers represent plots from adjacent natural vegetation (1–8) or from agroforestry systems (8–40). The words in red represent the plant species. Termac = Terminalia macroptera, Pteeri = Pterocarpus erinaceus, Piltho = Piliostigma thonningii, Teravi = Terminalia avicennioides, Comfra = Combretum fragrans, Vitpar = Vitellaria paradoxa, Proafr = Prosopis africana, Anaocc = Anarcadium occidentale, Parbig = Parkia biglobosa, Steset = Sterculia setigera, Comcol = Combretum collinum, Monker = Monotes kerstingii, Acapol = Acacia polyacantha, Lanmic = Lannea microptera.

The influence of regional species pools on agroforestry species assemblages

We found no significant relationship between the agroforestry system and adjacent natural habitat for species richness (β = 0.138 ± 0.259, P = 0.910) and Shannon diversity (β = 0.542 ± 0.348, P = 0.370; Table 1). However, the proportion of species richness shared between agroforestry systems and their adjacent natural habitats was positively affected by the ethnobotanical values (β = 0.362 ± 0.152; P = 0.024, R2 = 0.15, Fig. 4B, Table 2) but was not significantly associated with the regional species pool (β = − 0.426 ± 0.419; P = 0.318, R2 = 0.033, Fig. 4A).

Table 1

Bootstrap regression testing the linear relationships of species richness and Shannon diversity index between agroforestry systems and adjacent natural vegetation.

	Estimate	P value	BootBias	BootSE	BootMed	5%	95%
Intercept	1.576	0.380	− 0.103	0.435	1.530	0.963	2.395
Species Richness	0.138	0.910	0.087	0.259	0.153	− 0.374	0.477
Intercept	1.276	0.000	0.000	0.248	1.293	0.850	1.701
Shannon Diversity	0.542	0.370	0.042	0.348	0.542	− 0.078	1.083

Figure 4

Log–log regression illustrating (A) the non-significant linear decrease in the proportion of shared richness as function of regional species pool (SR) and (B) the linear increase in the proportion of shared richness as function of ethnobotany value (EV). SL is the species richness of agroforestry system.

Table 2

Parameterized models and selection of appropriate models using information theoretic approach.

Models	ΔAICc	df
Distance + ethnobotanical value	0	4
Distance + ethnobotanical value + crop type	1.8	5
Distance + ethnobotanical value + crop type + distance × crop type	1.9	6
Distance + ethnobotanical value + crop type + distance × crop type + ethnobotanical value × crop type	3.5	7
Distance	4.4	3
1	3.6	2

ΔAICc is the small sample corrected difference in Akaike Information Criterion (AIC) between each model and the lowest AIC value. df is the degree of freedom.

The role of human driven environmental filtering on species assemblages in agroforestry systems

The distance between agroforestry systems and neighboring natural habitats did not significantly affect the proportion of the shared species richness (β = − 0.012 ± 0.199, P = 0.950, Table 3). However, plot level species ethnobotanical value positively affected the proportion of the shared species richness (β = 0.701 ± 0.205, P = 0.001) indicating that farmers in this case play an important role in filtering species assemblages from the regional pool of the adjacent natural habitat. Furthermore, we found no significant effect of distance on species richness (β = 0.033 ± 0.093, P = 0.729, Table 4) and Shannon diversity index (β = − 0.012 ± 0.109, P = 0.656, Table 5). Similarly, ethnobotanical value had no significant effect on Shannon diversity index (β = 0.038 ± 0.055, P = 0.250) but positively affected species richness (β = 0.173 ± 0.080, P = 0.030).

Table 3

Beta regression testing the effect of the distance between agroforestry systems and their surrounding vegetation and ethnobotany value on the proportion of shared species.

Variables	Estimate	SE	Z	P
Intercept	− 0.524	0.307	− 1.702	0.088
Distance (scaled)	− 0.012	0.199	− 0.062	0.950
Ethnobotanical value	0.701	0.205	3.417	0.001**

Distance was standardized.

Table 4

Generalized linear model with negative binomial errors testing the effects of the distance between agroforestry systems and their surrounding vegetation and ethnobotanical value on agroforestry systems species richness.

Variables	Estimate	SE	Z	P
Intercept	1.078	0.143	7.545	< 0.0001***
Distance	0.033	0.093	0.359	0.719
Ethnobotanical value	0.173	0.080	2.167	0.030*

Table 5

Bootstrap regression testing the effects of the distance between agroforestry systems and their surrounding vegetation and ethnobotanical value on Shannon diversity index.

	Estimate	P	BootBias	BootSE	BootMed	5%	95%
Intercept	0.663	< 0.001	− 0.019	0.077	0.653	0.487	0.789
Distance	− 0.012	0.656	− 0.027	0.109	− 0.013	− 0.239	0.029
Ethnobotanical value	0.038	0.250	0.013	0.055	0.035	− 0.015	0.206

Discussion

This study examined species diversity patterns in agroforestry systems and their surrounding natural vegetation using the theory of environmental filtering and island biogeography theories as theoretical framework. The agroforestry systems around the Pendjari biosphere reserve have particularly low species diversity as compared to the richness of 60 species identified in semi-natural habitats and the 41 species in agroforestry systems in the nearby Vipalogo region of Burkina Faso[38]. This result may be explained by the increasing charcoal production and agricultural intensification within the transition zone that caused the fast decline in some species over time[37].

Clearly, our rank-frequency analysis showed lower diversity in agroforestry systems than in the natural habitats. Several potential mechanisms can explain such low diversity which may indicate strong environmental filtering of species from the regional species pools. First, the low diversity is perhaps due to the strong selective practice of farmers who prefer to conserve large size trees which provide fodder, fruits as famine food and several other non-timber forest products for local communities. Two main species are spared as a result, including Vitellaria paradoxa for its fruits and seeds used for shea butter and seeds of Parkia biglobosa used to make traditional mustard[37,39,40]. Small diameter trees, considered as non-desirable for this purpose, are instead harvested for fuelwood directly extirpating these species from the assemblages. Previous studies show that disturbance like fire and grazing can induce small-size trees (sapling or juvenile) mortality[41,42]. However, such effect was not evident in our study system since bush fires are not permitted in the biosphere reserve and grazing in the agroforestry systems take place for limited time and in very few places.

We hypothesized that site ethnobotanical values and proximity to natural habitats are positively associated with species diversity in agroforestry systems. We found support for the first hypothesis that ethnobotanical value was positively correlated with agroforestry systems diversity, but we found no significant correlation with distance. The non-significant effect of distance on species diversity in agroforestry systems is in contrast with predictions form the island biogeography theory. The theory predicts that islands/fragments closer to regional pools (mainland), which are more likely to receive more species dispersed from the regional pool, will have higher species diversity than the one that are farther away[28,29]. However, in our study, the benefit of dispersal in closer agroforestry systems was probably cancelled out by the strong human-driven filtering of species in these agroforestry systems. This is further illustrated by our finding that both agroforestry systems and their adjacent natural habitats were largely dominated by two agroforestry species (Parkia biglobosa and Vitellaria paradoxa). Farmers tend to conserve these preferred species because of their high nutritional, medicinal and commercial values[43,44]. The versatility of these species which can be used for multiple purposes make them an extremely valuable species to select for in the human modified landscapes[45,46]. Besides, the species filtering around the Biosphere reserve is also driven by agricultural policies that prioritize cash crops over traditional food crops[47,48] with strong economic and ecological drawbacks[49].

In addition to the lack of influence of farms proximity to natural vegetation on species assemblages in agroforestry systems, we also found no support for the effect of regional species pool. We hypothesized that if regional species pool had strong influence on agroforestry systems, then species diversity in agroforestry systems will increase with species diversity in neighboring natural habitats. Surprisingly, our results revealed a positive but non-significant correlation for species richness and Shannon index between the two systems, indicating a limited influence of regional species pool on agrobiodiversity. However, the proportion of shared richness was significantly and positively correlated with ethnobotanical values. These results, consistent with our previous analysis, suggest that species of higher ethnobotanical values are more likely to be conserved in agroforestry systems than species of lower ethnobotanical values. This finding highlights the key role of ethnobotanical knowledge in biodiversity conservation[50-52] and suggests that effort to increase agrobiodiversity must focus on promoting or adding plant species with high ethnobotanical values.

In this study, we investigated the drivers the species diversity patterns in traditional agroforestry systems around an African Biosphere reserve. We found no significant correlation between agroforestry systems and their adjacent vegetation for species richness and Shannon diversity index. We also found no significant effect of proximity to natural vegetation on plant richness in agroforestry systems. However, we found that ethnobotanical value was the key variable that drives the fraction of species diversity shared meaning that while converting natural vegetation into farm, farmers prefer to conserve species of high medicinal, nutritional, and commercial values. This underlines the central role that ethnobotanical knowledge plays in biodiversity conservation. Documenting and using indigenous knowledge on species of high ethnobotanical values in the ecological restoration of agroforestry system could stimulate farmers to conserve more species in their farms than usual. Such commitment on the part of farmers can sustainably improve crop yields (D. K. N’Woueni and O. G. Gaoue, unpublished data) and contribute to slowing down deforestation within the Biosphere reserve. We recommend ecological restoration efforts to prioritize species of high ethnobotanical values and particularly versatile plants that are used for multiple purposes.

Methods

Study area

The Pendjari Biosphere Reserve covers 4,661.4 km2 and is located in the northwest of Benin between 10° 30′ N–11° 30′ N and 0° 50′–2° 00′ E (Fig. 1). The annual rainfall ranges from 1000 to 1100 mm. The rainy season last nearly five months from mid-May to October followed by a dry season from November to February[53]. The vegetation is dominated by woodlands, savannas and grasslands established on poorly developed and tropical ferruginous soils[54]. The Pendjari Biosphere reserve is dominated by open forests and wooded savannas, wooded, shrub and grassy savannas. These plant communities are increasingly degraded by a rapidly growing population[55].The area is characterized by an annual average temperature of 27 °C. The minima and maxima are respectively of 21 °C and 40 °C[56]. There are rough or little evolved soils, tropical ferruginous soils and hydromorphic soils[54].

The Pendjari Biosphere Reserve is organized into three main zones: a core area or protected zone (The Pendjari National Park), a buffer zone and a transition zone where sustainable agriculture is permitted[57]. The Pendjari hunting zone is bordered by two main roads (Tanguieta-Batia road, and the Tanguieta-Porga road) which separate it from the transition zone[58]. Along the two main roads and inside the Pendjari hunting zone, local populations establish such a traditional agroforestry systems and gather non-timber forest products within the first 5 km perpendicularly from the roads[59]. Two main ethnic groups are established along the two main roads: the Berba ethnic group which accounts for 65% along the northwestern side and Wama and Gourmantche ethnic groups which account for 30% along the northeastern side[60]. Farms extension in the area occupied by the Wama-Gourmantche is constrained by the Atacora Chain of Mountains. Such a constraint makes the Wama-Gourmantche less environment conservators than the Berba[37,55].

Estimating diversity in agroforestry systems and natural vegetation

We conducted a floristic inventory in 32 traditional agroforestry systems and 8 adjacent natural habitats all distributed around 8 villages that border the transition zone of the Biosphere reserve. These villages were selected so that they are closer to the Biosphere reserve and based on the tree density and diversity of agricultural systems. Because the main goal of this study was to understand species assemblages in agroforestry systems, the villages were selected only if they had agroforestry systems that are used by local people. The first set of villages was selected from the Tanguiéta—Batia road and included Tchanwassaga, Pessagou, Tanougou, and Batia (Fig. 1). The second set of villages was selected from the Tanguiéta—Porga road and included Tiélé, Kani, Koupendri, Dassari. The size of the populations in these villages varies from 1169 inhabitants in Tiélé to 7472 inhabitants in Dassari with a mean of 2779 inhabitants per village. Four agroforestry systems were sampled in each village. For a given village, we established one 30 m × 30 m plot in each agroforestry parkland and one plot of the same size was installed in each adjacent natural habitat. The size of our survey plot was chosen based on decades of vegetation surveys in West African savanna which uses plot sizes between 500 and 1000 m2[61,62]. In the agroforestry plot, all crop and adult plant with diameter at breast height (DBH) ≥ 5 cm were identified and recorded. In the adjacent natural habitat, we also identified every plant species and estimated their abundance. The distance between the edge of each agroforestry system and the adjacent natural habitat was measured using a Global Positioning System (GPS).

We estimated trees density by counting individual woody plants within each plot. We used this data to develop a matrix of species abundance, each row representing a unique plot and each column representing a species recorded in our survey. We used this matrix to estimate species richness and Shannon diversity index[63]. Shannon index was estimated at the plot level using packages Hotelling and vegan in R version 3.1.2[64]. We estimated the plot level species richness with Chao abundance-based coverage estimators[65] using the function “estimateR” in package vegan. We also used the matrix of species abundance to construct the rank-abundance curves for agroforestry systems and their surrounding natural systems. The proportion of agroforestry systems species richness shared with their adjacent natural habitats was estimated as SL/(SR-SL) with SL the species richness of the agroforestry system and SR the species richness of the adjacent natural plant community.

Estimating ethnobotanical values

We estimate the total ethnobotanical values of each of the plant species surveyed in the agroforestry systems in the eight villages. We used a snowball sampling method, a sampling approach commonly used in ethnobotany[66-69] to select and individually interviewed a sample of 142 informants using a paper-based questionnaire in eight villages. An average number of 17.75 informants were interviewed per village. Informants were averagely aged of 42.92 years old. The majority (52.82%) were young (20 to 40 years) followed by 28.87% of adults (40 to 65 years old) and 18.31% of old (65 to 95 years). Less than half (40.85%) of informants was female against 59.15% of male. We obtained prior informed consent from informants before each interview and this study received the authorization from the University of Parakou Biomedical Research Ethic Committee (CLERB-UP). All the interviews were conducted in accordance with the relevant guidelines and regulations. For each site, we asked informants to list the different utilizations of the species inventoried in both systems (agroforestry systems and surrounding natural habitats). The ethnobotanical value was estimated for each site as the use value index[70] of species per plot which was computed with the package ethnobotanyR[71] in R. For each species, we asked each informant if the species was used or not for four different types of ethnobotanical uses (food, medicine, religious, construction). Species received 1 if they were used for a specific use category and 0 if not. The use value index for each plant species was calculated following Phillips and Gentry[70] as the sum of the total number of use citations by all informants for a given species, divided by the total number of informants who cited this species. We calculated the use value index at the plot level as the average use value index of all the species recorded in this plot. This represents the ethnobotanical value of the site. When this value is high, the plot has high number of species that are valuable for various purposes including medicinal, food, spiritual and construction uses.

Statistical analysis

We performed the Nonmetric Multidimensional Scaling (NMDS) using the vegan package[72] in R version 3.1.2 to assess the degree of dissimilarity in species composition between agroforestry systems and their adjacent natural community. To further understand how these two plant communities are different, we conduct a rank-frequency analysis, plotting species frequency in our survey plots against their ranks. To test the relationship between species richness or Shannon diversity index in agroforestry systems and in adjacent natural habitats, we used linear bootstrap regressions with 1000 iterations using the package boot[73] in R version 3.1.2. Bootstrap regression was also used to test the effect of the type of crops systems, ethnobotanical value and the distance between agroforestry systems and nearby natural habitat on the Shannon diversity index. To test the effects of the three variables on agroforestry systems species richness (count data), we performed a generalized linear model with a negative binomial error structure given the overdispersion of the species richness. We built several candidate models including all single and interactive effects. For each candidate model we calculated the small sample size corrected Akaike Information Criterion, AICc[74]. We then calculated for each model, their ΔAICc which is calculated as the difference between AICc of each candidate model and the smallest AICc value obtained. We selected the best models as the models with ΔAICc < 2[75].

In addition, we tested the log–log linear relationship between the proportion of shared species richness and the regional species richness to investigate the influence of regional species pools on local species assemblages in agroforestry systems. We used generalized linear models with a beta error structure to test the effect of the distance from agroforestry system to the adjacent natural habitat, ethnobotanical value, and the type of crops systems on the proportion of shared species richness between agroforestry systems and natural vegetation. We parameterized several candidate models including interactions terms. Similarly to previous models, we used small sample size corrected Akaike Information Criterion[74] to select the best models.

Informed consent

This study received the authorization from the University of Parakou Biomedical Research Ethic Committee (CLERB-UP). Prior informed consent was obtained from each participant before the interview.