Modeling Bat Species Richness and Spatial Distribution in Burkina Faso ()
1. Introduction
Bats are one of the most ecologically diverse groups of vertebrates. They belong to the order of Chiroptera, which has a wide variety of species. It is the second largest group of mammals after rodents. Currently, there are 1.300 bat species described worldwide; this represents about a quarter of known mammals [1]. Bats are considered keystone species in maintaining ecosystem functions [2] [3] [4] [5]. Some bat species are used by local people in their dietary as protein complement and others are useful in medicine [6] [7] [8]. They are also considered as hosts of some parasites which can cause many diseases for human being. Indeed, there are some negative perceptions associated with their zoonotic potential [9] [10] [11]. Some authors argued that the loss of bats could increase human reliance on chemical pesticides, jeopardize ecosystems and harm human economies [12] [13].
Despite the important ecological role of bats worldwide, they are often neglected and few studies have interested in them especially in Africa. As a result, our knowledge of these species is very limited and therefore, does not allow us to protect them and exploit the ecosystem services they can provide. Till the first mention about bats from Burkina Faso in 1969 [14] little is known about these mammals. Nevertheless, important insight given till 2012 [15], based on a large sampling area in 32 localities described 45 species with 15 new records [16] and [17] respectively described the distribution of species and the population’s structure in Burkina Faso. It has been demonstrated that environmental variables influence bat species’ spatial distribution. An adequate understanding of their ecology and distribution area is important in the development and implementation of policy and management actions that are crucial to the long-term survival of these species, especially in Burkina Faso. In the current context of terrorism in Burkina Faso supported by the Covid-19 pandemic, it is becoming increasingly difficult to explore certain sites for biodiversity inventory. Fortunately, Species Distribution Models (SDMs) demonstrated their utility to predict species’ fundamental niche. So, this study was conducted using SDMs to predict bat species richness throughout Burkina Faso to answer three main questions. 1) Is the hotspot of predicted bat diversity corresponding to the protected areas in Burkina Faso? 2) Are they some areas where additional sampling effort can provide more species? 3) What are the main bioclimatic variables that govern the spatial distribution of bat species in Burkina Faso? Addressing these questions, will be of high importance for a) Supporting bats conservation policies; b) Identifying further sampling areas; c) Monitoring biodiversity specifically in this context emerging zoonotic diseases like Covid-19.
2. Methodology
2.1. Study Area
This study covers the entire area of Burkina Faso, a landlocked Sahelian country located in West Africa with an area of 274,200 km2. Most (about 75%) of its territory is based on a Precambrian crystalline basement which gives it a generally flat relief [18]. Burkina Faso’s hydrographic network is enough relatively dense despite the precarious weather conditions. The country is characterized by a Sudano-Sahelian tropical climate, in general, with two seasons: the dry season from October to April and the rainy season from May to September. The larger part of the country is located in the Sudanian climatic zone, including the Center and the South. The northern part of the country is influenced by Sahelian climate [19]. Guinko [20] and Fontès and Guinko [21] distinguish two large phytogeographical areas on the basis of climate, vegetation and fauna; the Sahelian and Sudanian areas, each being divided into two sectors. The vegetation of the Sahelian phytogeographical area includes shrub and tree steppes, grassy steppes, contracted training (tiger bush), riparian and riverine cords formations [22]. The Sudanian phytogeographical area locates south of the thirteenth parallel. The vegetation is characterized by a set of savannas (from woodland to grassland).
2.2. Data Collection
The occurrences (geographical coordinates of localities where the species were observed) and environmental variables are the two types of data used for the SDM. Occurrences used come from a literature review (covering the period of 1964 to 2008) and a field collection from 2009 to 2012. They cover all the national territory (Figure 1). Environmental variables including climate data, topography, land cover and hydrography were collected from various sources (Table 1). Climate data were downloaded from WORLDCLIM version 1.4 [23], land cover data derived from MODIS Vegetation Continuous Fields [=VCF] collection 3 [24] and SPOT VGT mosaic Africa data 2000 [25]. Topography data come
Figure 1. Map of Burkina Faso showing the vegetation zones and the bats collection points.
Table 1. List of environmental variables.
from Shuttle Radar Topography Mission (SRTM) [26]. Data on the river system from the Geographic Information System for Burkina Faso (BF GIS) developed within the framework of the project “SFB 268—History of cultures and languages in the natural area of the African Savannah” to the Institute of Physical Geography of the Johann Wolfgang Goethe University in Frankfurt am Main.
2.3. Modeling Approach
The maximum entropy theory was used through the Maxent software version 3.4.0 for the SDM. Maxent is a general method based on statistical analysis to predict the best environmental conditions for species, from environmental variables and species occurrence. We applied MaxEnt’s default function with 10 replicate cross-validation. Jackknife test was also applied to evaluate the variable importance. To create the Species Richness Distribution (SRD) we first created individual SDM from all the 44 species which present each at least two occurrences. For each individual SDM, eighteen environmental variables were used because of their ecological importance in bats distribution. These are ten from climate variables, five from vegetation cover, two from topography and one from hydrography (Table 1). The spatial resolution was harmonized to 30 arc seconds in order to have a high resolution which is necessary to give best information per grid cell.
The species occurrence data were organized as follows: the national territory has been divided into a grid of approximately 1 kilometer. In cases where we had several records in the same grid, we chose the best records taking into account the quality of the coordinates (taken with a GPS or not), the reliability of the species identification (qualification of the identifier) and the age of the data (old or recent). After using that, if there are still cells with more than one record, we choose one of them by centering it in the grid. Then, records that have not been selected for the training will be used as test points. For species with less than ten records, all records were used for training (Table 2).
Maxent was set up to perform 10 replicates and their average value was used to make the modeling. Area Under the Curve (AUC) was determined for model evaluation and was interpreted following [27]. In order to determine the predictive power of each environmental variable considered in the distribution of species and to identify those that contribute most to the model, the Jackknife test was also performed.
Individual models (by species) obtained with MaxEnt were then imported into R software version 4.0.2 to generate discrete distribution maps for each species as the Maxent output consists of a continuous probability map (from 0 to 1). The division is made from the smallest of the following three probability thresholds obtained with Maxent: the “10 percentile training presence”, the “Equal training sensitivity and specificity” and the “Maximum training sensitivity plus specificity”. Then, the discrete maps for species are added to give the SRD by family, and the overall SRD for the country by using R software version 4.0.2 with package “raster”, “sp” and “rgdal”.
Table 2. List of bat species and the number of used records.
For each species, Maxent gives the contribution of each environmental variable to the SDM. We, therefore, estimated the weight of each environmental variable in the SRD by averaging the contributions of each variable by family and for the entire bat community. The synoptic figure below (Figure 2) shows the main steps of the modeling process.
3. Results
3.1. Distribution of Bat Species Richness by Families and Model Quality
The MaxEnt model had a moderately high performance by presenting medium value of 0.875 and 0.854 respectively for the training data and testing data (Figure 3). The AUC values obtained in this study varied from 0.638 for Mops major to 1.00 for Rousettus aegyptiacus and Doryrhina cyclops. These values of AUC indicate that the SDMs we generated achieved good levels for prediction.
Based on the available occurrence data, we modeled the distribution of 44 of the 51 bat species encountered. Modeled species belong to 9 families. SRD for those families are shown in Figure 4 and the means of environmental variables contribute to the SDM are presented in Table 3. Species richness of Pteropodidae, Hipposideridae, Molossidae, Nycteridae, Rhinolophidae and Vespertilionidae shows clear increases from north to south.
The distributions of Pteropodidae, Hipposideridae, Molossidae, Nycteridae, Rhinolophidae and Vespertilionidae follow clearly the climatic areas with higher richness in the Sudanian region. Pteropodidae reach the maximum richness in the hilly region in the extreme south-west part of the country. Then, maximum precipitation and topography are the two main variables that contribute most to the distribution of this family. The richness of Pteropodidae varies in the
Figure 2. Illustration describing the scope of species richness modeling.
Table 3. Mean of environmental variables contribute to the SDM by family and for the chiropteran community in Burkina Faso. PT: Pteropodidae; EM: Emballonuridae; HI: Hipposideridae; ME: Megadermatidae; MO: Molossidae; NY: Nycteridae; RHL: Rhinolophidae; RHP: Rhinopomatidae; VE: Vespertilionidae; All: All Chiropterans.
Figure 3. AUC value according to training and test.
Figure 4. Maps of the bat species richness distribution by families.
same direction as the two variables, then these variables have a positive effect on the family richness. For Hipposideridae, the two most important variables, with positive correlation, are the standard deviation of precipitation and total annual precipitation. The highest richness in the south-west opposite to Pteropodidae doesn’t correspond to the hilly region and covers up all of this area. Richness of Molossidae shows maximum in the follow the center of the South Sudanian area, in and around the protected area of Nazinga. Hydrography (rivers) are clearly important for this family in the Sahelian region. Minimum precipitation and bare ground cover are the two variables that contribute the most to their distribution. While bare ground cover shows a negative correlation, the minimum precipitation shows a positive correlation. Nycteridae presents the widest area for maximum richness, covering almost the Sudanian zone. In the Sahel region, the rivers play a positive and critical role lead to the presence of more species. But in the Sudanian Region Rivers are likely to decrease the number of species. Then, maximum precipitation and hydrography are shown to be the most important variables. Maximum precipitation shows a positive correlation with the Nycteridae diversity. The extreme south-west is a high diversity area for Rhinolophidae that are absent in the Sahelian regions. The standard deviation of maximum temperature and standard deviation of precipitation are the two main variables that contribute most to Rhinolophidae distribution. But, as distribution of Rhinolophidae shows opposite trend compared to the standard deviation of maximum temperature, we infer that this variable negatively influences their Rhinolophidae distribution. Vespertilionidae is the most diversified family in Burkina Faso, nearly present in all the country but with richness higher than 2 only in the Sudanian region. Highest richness (9 to 10 species) is found in the extreme south of the country at the border with Ghana and Cote d’Ivoire. Rivers are also an important element in the distribution of these species as they determine their presence in the Sahel and even in some parts of the Sudanian zone. Topography also seems an important factor in the extreme south-west, where part of the maximum wealth is observed. The two main environmental variables that contribute most to the distribution of Vespertilionidae are total annual precipitation and standard deviation of minimum temperature. Total annual precipitation shows a positive correlation and the standard deviation of minimum temperature shows a negative correlation to Vespertilionidae distribution.
In opposite to the previous families, Emballonuridae richness increases from south to north. The minimum of minimum temperature and minimum precipitation are the two main variables that contribute most to their distribution. These variables are negatively correlated to the Emballonuridae richness. Megadermatidae is a family with single species. Bare ground cover and standard deviation of minimum temperature are the two variables that contribute the most to this family distribution. These two variables show a negative correlation to the distribution of Megadermatidae family. Rhinopomatidae only present in the highest land of the country. Then topography and annual average of spectral response values in the Near-Infrared, band 2 (vegetation cover) are the two main variables that contribute most to their distribution. The topography therefore positively influences the distribution of Rhinopomatidae, unlike vegetation cover which negatively influences the distribution.
3.2. Overall Pattern of Species Richness Distribution in Burkina Faso
In Burkina Faso, bat species richness increases from north to south with a maximum richness in the hills of south-west region (Figure 5). The general pattern
Figure 5. Bat species richness distribution in Burkina Faso.
shows that we can expect at least one species of bat anywhere in Burkina Faso. Changes in species richness superimpose perfectly to the climatic zones of Burkina Faso. By and large, in the Sahelian region, we can expect a maximum of 7 species with more species (from 6 to 15) following the river beds. The “bands” of richnes from 1 to 10 species, share the North-Sudanian zone while the “band” of 11 to 25 species lies troughout the South Sudanian region. The richnes of 26 to 44 species appear in the South Sudanian region, precisely at the most rugged sites in the country whose cliffs of Banfora and mountain ranges.
In order of importance, total annual precipitation, topography, bare ground cover and standard deviation of maximum temperature are the most contributing variables to the overall bats SRD in Burkina Faso (Table 3). Total annual precipitation and topography are positively correlected to bat species diversity while the bare ground cover and standard deviation of maximum temperature are negatively correlated to bat species diversity.
4. Discussion
4.1. Model Quality and Variables That Govern Bat Species Richness in Burkina Faso
The SDMs results showed a good performance of predicting models. Indeed, according to Araujo et al. [27], the results of SDMs are considered to be good when the AUC value is equal or more than 0.75. Bat species richness in Burkina Faso varies along the north-south gradient, with a positive correlation to the total annual precipitation and topography, the most important variables explaining the species distribution. As already shown, the diversity of bats increases with environmental heterogeneity and complexity of the habitat [28] [29]. According to Tews et al. [30], the majority of studies show a positive correlation between habitat heterogeneity, as the horizontal development of vegetation [31], and diversity of animal species. Such habitat heterogeneity increase with the precipitation level that determines vegetation and the global weather and reaches its highest level in the south Sudanian regions of Burkina Faso. This also explains why the bare ground cover negatively influences bats’ diversity, as this indicates less habitat heterogeneity. Indeed, from the North to the South of Burkina Faso, there is an increase in rainfall and therefore an improvement of the vegetation covers with the passage of shrub-steppe and wooded, grassy steppes to woodland to grassland [32]. This state is supported by the fact that some species as Eidolon helvum migrates annually into savannah regions along a north-south axis [11]. In addition, the South of Burkina Faso presents a set of specific habitats such as the many protected forests, galleries forest that increase the heterogeneity of environment. Seeing the habitat complexity as the vertical development of the vegetation [31], it is clear that habitats complexity reaches its highest level in the Sudanian regions, mainly in the extreme southwest, much richer, where humid Guinean flora appear along the main rivers [33]. We also note that in the Sahelian region of Burkina Faso, rivers contribute to increasing species richness. This is related to the fact that rivers, with their associate gallery forests, contribute to increasing habitat heterogeneity and complexity, thus offering specifically suitable habitats for some species. The rock formations like the cliffs of Banfora, Sindou peaks and the chain of Gobnangou increase habitat heterogeneity and complexity and offer additional lair for bats through many cracks and caves which are home and justify the weight of the topography in the determination of the bats’ diversity. Large variations in the maximum temperature negatively affect the species richness of bats in Burkina Faso. This indicates the sensibility of bat fauna to the potential increase of temperatures and suggests a great vulnerability of these animals to global warming. In the same sense as this observation, Welbergen et al. [34] found that extreme temperatures are important threats to bats in Australia.
At the family level, three distribution patterns emerge. The diversity of Pteropodidae, Hipposideridae, Molossidae, Nycteridae, Rhinolophidae and Vespertilionidae follow the general trend increasing globally from North to South. Quite the reverse, the diversity of Emballonuridae decreases from North to South. The two species of this family (Taphozous nudiventris and T. perforatus) are all known from the north of the Sahara [35]. Then, we can suppose that the distribution area of these species is limited to the Sahelo-Saharian zone. The Rhinopomatidae family (Rhinopoma cystops and Rhinopoma microphyllum), although also present north of the Sahara [35], are strongly and positively related to the topography in Burkina Faso.
Although these variables show a good contribution to the models, we cannot assume that there are the only factors affecting bat diversity. Indeed, many studies like Hawkins et al. [36], Kalko [37] and Kalko [38] [39] showed that variables as the availability, abundance and distribution of food resources are also important factors that influence the organization and dynamics of bats.
4.2. Hotspot of Bat Diversity and Protected Area in Burkina Faso
The highest diversity of bats (from 36 to 44 species) appear in the most rugged area in the country with only narrow protected areas like the classified forests of Beregadougou, Denderesso and Koulima [40]. Sites like the cliffs of Banfora, the peaks of Sindou and the mountain range of Gobnangou where the highest diversity has been predicted are not officially protected sites. On the other hand, areas of good diversity (36 to 44 species), and that of fair diversity (16 to 20 species) are dotted with protected areas. This bodes well for optimal protection for the species found there. Thus, all the protected areas correspond fairly well to areas of high diversity of bat. Better, specific areas such as the Sahel reserve although located in an area of low diversity guarantee the heterogeneity and complexity of protected habitats. However, this complexity of protected habitats would benefit from being supplemented by the classification of larger areas of the Banfora cliffs, the Sindou peaks and the Gobnangou range.
4.3. Importance of Additional Sampling Effort in Burkina Faso
Kangoyé et al. [17] reported a specific richness of 1 to 12 species per sampling site after having compiled data from 164 sites visited over a period of 50 years in Burkina Faso. Similarly, Sawadogo [41] only encountered 11 species in the whole city of Ouagadougou. Several sites covered by these works have higher predicted specific richness levels. This suggests that it is possible to meet even more species there. Indeed, after having listed 51 species for Burkina Faso, Kangoyé et al. [17] using the rarefaction curves and the indices of Jackknife and Chao 2 and Boot underline that the total specific richness could reach 55 to 58 species. Thus, 4 to 7 species remain statistically likely to be encountered and the results of this work allow us to direct our future efforts towards areas of high heterogeneity and high complexity of habitats as well as towards protected areas. This is confirmed by the recent work of Thiombiano et al. [42] that recorded a new species (Taphozous mauritianus E. Geoffroy St.-Hilaire, 1818) in the city of Bobo-Dioulasso. The models of bats richness species produced in this study constitute a useful tool for their conservation management. Because of their economical and environmental benefits that they provide, conservation strategies should strive to maintain bat richness and their population.
5. Conclusion
This study employed species distribution modeling to determine the distribution of bat species richness in Burkina Faso. It shows that bats are everywhere in the country, with increasing diversity from north to south with the highest diversity in the hilly region of the south-west. Using the Jackknife test, we identified that variables as precipitations, topography, vegetation cover, and rivers bed are the principal determinant of the species richness distribution. Moreover, in the context of Burkina Faso diversity hotspot, mainly comprise the protected area. Hilly and rocky formations encountered in the southwest part of Burkina Faso host the largest bat diversity of the country and highlight the importance of this area in the conservation of biodiversity in Burkina Faso. In view of the ecological and therapeutic roles of bats already mentioned in the literature review and in view of their specific richness in Burkina Faso, it is obvious that the country has a great deal of biological wealth that should be preserved. Future studies could explore their multiple uses by local communities. This could contribute to their better sustainable management. This area of vital importance to bats of Burkina Faso should receive more attention from all stakeholder interested in the conservation and protection of biodiversity areas of important bats biodiversity like Banfora cliffs, the Sindou peaks, and the Gobnangou range need more protection measures. Finally, after 51 species further sampling effort in these areas and the other protected areas from the east and the south of Burkina Faso could mead to new records.
Acknowledgements
This article is dedicated to the memory of Professor Elisabeth K. V. Kalko. We thank Jacob Fahr for his support during this project. We are grateful to “BIOTA West Africa” funded this research.