Abstract

Bale Mountains National Park is one of the protected areas in Ethiopia that holds the largest area of Afroalpine habitat in Africa and the second largest stand of moist tropical forest. Nevertheless, human settlements, overgrazing, and recurrent fire are the main problems in the park. This study aimed to determine the effects of human-induced landscape change in floristic composition and structure in the park. The vegetation data were collected systematically from 96 sample plots laid along 24 line transects in the edge and interior habitats of the six land cover types. Vegetation composition and landscape structural analysis were made using R software version 3.5.2 and FRAGSTATS version 4.2.1, respectively. Patch number was strong and positively affected species richness (r = −0.90, ), diversity (r = −0.96, ), and basal area (r = −0.96, ), whereas mean patch size was strong and negatively influenced species richness (r = 0.95, ), diversity (r = 0.87, ), and basal area (r = 0.82, ). The overall species richness, Shannon diversity index, and Margalef index were significantly higher in the edge habitat; however, the mean basal area of woody species was significantly higher in the interior habitat at . This study uncovered that the park is floristically rich and diverse, and it provides a variety of ecological and economic benefits to the surrounding community and to the nation at large. However, these benefits are gradually declining due to the high level of anthropogenic activities in the park. Thus, integrated environmental management strategy that blends with sustainable use of natural resources should be implemented to minimize the threats.

1. Introduction

Landscapes all over the world are alarmingly changed and fragmented due to anthropogenic factors such as urbanization, agricultural expansion, forest fire, and climate change [1, 2]. Most of the global changes responsible for the reduction in population and biodiversity are exacerbated by fragmentation [3, 4]. The primary causes of global biodiversity reduction are the destruction and degradation of natural ecosystems [5]. Predominantly, habitat loss and fragmentation are presently the main threats to terrestrial biodiversity [6]. Moreover, habitat fragmentation can affect the species interactions and community composition, as invasive or pest species, and may substitute the original species pool and increase the transmission and prevalence of the disease in small fragments [7]. Moreover, the species richness and abundance usually decrease with reduced patch size [8]. As landscapes become more fragmented, patch diversity increases with subsequent increase in the edge, exotic, and generalist species and ultimately leads to the reduction in landscape quality as habitat for species [9]. Accordingly, species richness in interior habitat, particularly indigenous and specialist species, tends to decrease [10]. The number of species existing in a patch tends to rise with patch size up to a certain limit, and the types of species found also tend to vary in size [8]. Size and shape interact to influence the amount of interior area remaining in a particular habitat fragment [2].

Tropical montane ecosystem is one of the hot spot ecosystems on Earth that comprises more than 200,000 species of flowering plants [11, 12]. The Ethiopian highland, which is located in the tropical region, encompasses over 50% of the Afromontane vegetation in Africa [13]. A suitable geographical position, a wide range of altitude, a high amount of rainfall, and a wide range of temperature variations equip the area with huge ecological diversity and a wealth of biological resources [14]. However, severe deforestation coupled with the cultivation of steep marginal lands, overgrazing, and sociopolitical uncertainty has resulted in rigorous land degradation over large areas of the country [15]. The overdependence of the country’s economy on agricultural production and the existence of more than 80% of the population in the highlands [16, 17] mainly contribute to the degradation of ecological resources and biodiversity loss.

The mountainous landscape and the mosaic of natural vegetation in the Bale Mountains have considerable economic, recreational, esthetic, and scientific importance [14]. The Bale Mountains National Park (BMNP) is the most significant conservation area situated in this region of Ethiopia and established in 1969 to preserve the endemic and indigenous floras and faunas in the area [18, 19]. It is one of the 34 International Biodiversity Hotspots and meets the requirements for the World Heritage Site and Biosphere Reserve Listing [20]. However, the park is facing a critical challenge from the illegal settlement and overgrazing and that leads to the change in its landscape structure and function. As a result, the habitats in the park are changing and the provision of ecological services from it is substantially reduced. Consequently, no research provides detailed information about the landscape structure and its potential impact on vegetation composition and structure in the park. Therefore, this research was aimed to analyze the potential impact of landscape change in floristic composition, diversity, and structure in the BMNP. Particularly, a comparative analysis was made among the edge and interior habitats of the park.

2. Materials and Methods

2.1. Study Area Description

BMNP is located within the geographic bounds of 6°53′08″N latitude and 39°44′03″E longitude and 400 km southeast of Addis Ababa, Ethiopia (Figure 1). It comprehends a wide range of habitats between 1450 m and 4377 m altitude. The park holds the largest area of Afroalpine habitat (about 1000 km2) above 3000 m asl in Africa and the second largest stand of moist tropical forest [21]. It is one of the 34 International Biodiversity Hotspots and also qualifies for World Heritage Site and Biosphere Reserve Listing [22]. It received rainfall that ranged from 520 to 2370 mm annually [23], and the distributional pattern is bimodal with heavy rains from July to October (highest peak in August) and small rains from March to June (highest peak in April). The mean monthly minimum and maximum temperatures are 5.6°C and 21.4°C, respectively. Its soil is fertile silty loam of reddish-brown to black clay soils dominated by Vertic Cambisols and Leptosols [24].

2.2. Vegetation Sampling Design

From 13 to 20 November 2018, a reconnaissance survey was conducted to get insights into the vegetation physiognomy and establish sampling sites in the study area. Following, the actual fieldwork was performed in the dry season between November 2019 and January 2020. A total of 96 sample plots (20 × 20 m) were systematically laid along 24 line transects in eight directions along three altitudinal gradients at 100 m elevational differences as it maximizes the distance between plots and minimizes spatial correlation among the observations [25]. To make a comparison between the vegetation data, an equal number of sample plots have been laid on the edge and interior habitats following Daye [26].

2.3. Species Identification

Plant species in the nested plots were identified at the field with the help of local peoples (for vernacular names) and by referring different volumes of Flora of Ethiopia and Eritrea books [27, 28]. For the species that were difficult to identify in the field, representative specimens were cut, numbered, and pressed at the site. The collections were named using folk taxonomy, and identification of formal taxonomy was determined using the voucher specimens at the National Herbarium, Addis Ababa University.

2.4. Floristic Composition and Structure

The most commonly used diversity indices of species richness (S), Simpson index (D), Shannon–Wiener index (H′), Pielou’s evenness index (J′), Whittaker β-diversity (), Margalef index (DM), and Berger–Parker index (d) were computed to analyze the patterns of plant diversity in the edge and interior habitats following Magurran [29] and Økland [30] using equations (1)–(4):where H′ is the Shannon diversity index, Pi is the proportion of individuals, and ln is the natural logarithm.where Hmax is the maximum level of diversity possible within a given population (ln S) and S is species richness.where a is the number of shared species in two sites and b and c are the numbers of species unique to each site.

The Margalef diversity index (DM) was computed using the following formula:where DM is the Margalef diversity index, S is the number of species, and N is the total number of individuals in the sample.

The woody species density, frequency, dominance, and their relative values in the interior and edge habitats were computed to obtain the important value index and describe the woody species structure following Ellenberg and Mueller-Dombois [31] and Martin [32] using equations (5)–(8). Moreover, DBH, tree height, and basal area were analyzed to determine the population structure following Kitessa et al. [33] and Van der Maarel [34]:where BA is the basal area, π = 3.14, and d is the DBH (cm).where Fr is the frequency of a species and Pi is the number of plots in which the ith species occurred.where Rde is the relative density and ni is the number of individuals of the ith species.where IVI is the importance value index, Rde is the relative density, RFr is the relative frequency, and RDo is the relative dominance.

2.5. Measurement of Landscape Structure

Landsat images of the years 1985, 1995, 2005, and 2017 were processed using ArcGIS version 10.3 to produce time-series datasets of land use/land cover. Then, eight landscape indices were analyzed using the processed land use/land cover data following McGarigal et al. [35] and Smiraglia et al. [36]. The indices include patch number (PN), mean patch size (AREA_MN), total core area (TCA), edge density (ED), area-weighted mean shape index (SHAPE_AM), mean Euclidean nearest neighbor distance (ENN_MN), and interspersion and juxtaposition index (IJI). Edge habitat was identified by deducting 50 m from the edge of each vegetation type. FRAGSTATS software version 4.2.1 was used to compute the landscape patterns in each land cover class and the entire landscape [37]. The two-way analysis of variance (two-way ANOVA) and linear regression analysis were made to test significant differences between fragmentation indices and species composition and structure parameters following the post hoc Tukey’s highly significance difference (Tukey’s HSD) test at 5% significance level using PAST software version 4.02 [38].

3. Results and Discussion

3.1. Landscape Structure Change

The analysis of landscape structure in this study revealed that the habitats in the BMNP are progressively transformed. The area has shown an increase in PN by 40.2% and a decrease in AREA_MN by 28.7% from 1985 to 2017. According to Oertli et al. [39], the high number of separated patches in a habitat indicates a high level of fragmentation. Across the entire study period, SHAPE_AM, which indicates the complexity of patch shape, increased by 18.8%. A higher perimeter-area relationship characterizes the rapid rate of fragmentation in the landscape [40]. Moreover, there was inconsistency in the values of ED; however, it was increased by 22.3% over the study period. As it was emphasized by McGarigal [37], the oscillation of ED indicated a major reduction in the spatial heterogeneity of the landscape. Conversely, the study area has shown a declining trend in TCA by 10.6% from 1985 to 2017. This was due to the escalated level of disturbances in the study area. As it was reported by Kidane et al. [41], the most dominant practices in the Bale Mountains, especially after 1995, were the upward expansion of agriculture and enrichment plantation.

The isolation of patches within the landscape of the study area was increased from 105.22 m to 111.94 m overtime (Table 1). This result is in agreement with the result reported by Tolessa et al. [42] in the central highlands of Ethiopia and Daye [26] in Southwest Ethiopia. Conversely, the intermixing of patches in the study area showed an overall declining trend from 95.38 to 86.77 over the study period. This result showed that the BMNP constitutes more scattered patches compared to other similar areas studied by Posada Posada [43] and Tolessa et al. [42].

3.2. Overall Floristic Composition and Structure

A total of 205 plant species belonging to 71 families and 153 genera were recorded (Table 2). Of these, 50 species were trees, 52 were shrubs, 12 were lianas, and 91 were herbs. Asteraceae was the most dominant family with 31 species, followed by Fabaceae with 11 species. Conversely, Helichrysum was the most abundant genus with 9 species, followed by Alchemilla and Trifolium with 5 species each. Twenty endemic species, including Euphorbia dumalis S. Carter, Lobelia rhynchopetalum Hemsl., and Thymus schimperi subsp. Schimperi Ronniger was identified in this study. The overall Shannon diversity and evenness index of the study area were 4.34 and 0.81, respectively. This indicated that the study area was more diverse compared to other similar vegetation areas including Bonga forest [44], Agama forest [45], and Munessa forest [13]. Conversely, the total density of seedlings, saplings, and mature trees in the study area was 8751, 4413, and 1567 individuals ha−1, respectively. This was lower than other comparable areas such as Kuandisha forest [46] and Wof-Washa forest [47]. The ratios of seedling to mature tree, sapling to mature tree, and seedling to sapling were 5.58, 2.82, and 1.98, respectively. This shows the recruitment potential of the forest is relatively higher [48].

Woody species density with DBH°>°2 cm was 1567 individuals ha−1. This was relatively higher compared to other similar vegetation areas such as the Wof-Washa forest [48] and Agama forest [45]. The most frequent woody species was Croton macrostachyus Hochst. ex Del with 81% frequency followed by Juniperus procera L. (79%), Podocarpus falcatus (Thunb) C.N (63%), and Hagenia abyssinica (Bruce) J.F. Gmel (60%). Conversely, the total basal area of woody species was 170.26 m2 ha−1, and it was considerably higher compared to other similar vegetation areas in Ethiopia. About 75% of the basal area was contributed by five tree species such as Juniperus procera (46.71 m2 ha−1), Syzygium guineense (Willd.) DC (24.76 m2 ha−1), Cordia africana Lam (20.95 m2 ha−1), Hagenia abyssinica (18.47 m2 ha−1), and Ehretia cymose Thonn (15.86 m2 ha−1). Consequently, Juniperus procera was the dominant woody species with an IVI of 26.43. The species with higher IVI values in the study area was among the characteristic species in similar vegetation types in Ethiopia [49, 50].

3.3. Floristic Composition and Structure in the Edge and Interior Habitat

A total of 136 species belonging to 111 genera and 59 families were identified in the edge habitats of the sampled patches, whereas 117 species that belonging to 84 genera and 40 families were recorded in the interior habitats. From the identified life forms, 19 species were trees, 22 species were shrubs, 86 species were herbs, and 7 species were lianas in the edge habitats, whereas 28 species were trees, 21 species were shrubs, 57 species were herbs, and 11 species were lianas in the interior habitats. The overall means (±SE) species richness (35 ± 4.2), Shannon diversity index (2.93 ± 0.17), and Margalef index (5.68 ± 0.69) of the edge habitat were significantly higher compared to the interior habitat at (Table 3). These variations could be due to the differences in site productivity, habitat heterogeneity, and disturbance factors [44, 51] or the invasion of exotic plant species [52]. However, the woody species richness in the interior habitat (28) was significantly higher than the edge (17). Moreover, the evenness index in the interior habitats (0.83 ± 0.04) was higher, but not significant, than the edge habitat (0.79 ± 0.05). This result was in agreement with the finding in [53]. Abiotic factors, seed predation, loss of pollinators and seed dispersers, and tree mortality were reported as the common causes for the differences in woody species composition between the edge and interior habitats [53, 54]. The computed Sorensen’s similarity index depicted that the number of species in the edge habitats was 45% similar to the species in the interior habitats. This value indicated that the similarity between the edge and interior habitat was weak [13]. The mean density of seedling (995.42 ± 19.27 individuals ha−1), sapling (509.29 ± 9.06 individuals ha−1), and mature trees (187.60 ± 4.70 individuals ha−1) in the interior habitat was significantly higher compared to the edge. This indicates that the recruitment potential of the interior forest was significantly higher compared to the edge habitat [48]. This could be due to the increased mortality rates of seedling, sapling, and mature trees in the edge habitats [53, 55].

The mean woody species density in the interior habitat (85 ± 22.17 stems ha−1) was significantly higher compared to the edge habitat (70 ± 16.53 stems ha−1) at (Figure 2(a); Tables 4 and 5). This could be due to the selective cutting of trees for timber production, house construction, and firewood in the edge habitats, which ultimately leads to a reduction in the density of large trees and greater canopy openness [56]. Moreover, the seedlings are most affected by edge effect due to their sensitivity to environmental changes and biotic interactions [57, 58]. Conversely, the mean basal area of woody species in the interior habitat (11.16 ± 1.82 m2 ha−1) was significantly higher than the edge habitat (3.99 ± 0.54 m2 ha−1) at (Figure 2(b); Tables 4 and 5). This was due to the significantly higher mean DBH (78.62 ± 4.56 cm, ) and height (33.63 ± 2.71 m, ) of woody species [59] in the interior habitat than the edge. There were 27.32% of larger diameter individual tree species with DBH > 100 cm recorded in the interior habitat, whereas 4.09% of individuals with DBH > 100 cm were identified in the edge habitat. Microenvironmental conditions such as high temperature, low relative humidity, high wind force, low soil nutrient, and litter moisture in the edge habitats may contribute to the changes in tree abundance and distribution in the forest [60, 61].

Juniperus procera was the dominant woody species in the edge habitat with an IVI of 32.49, whereas Croton macrostachyus was dominant in the interior habitat with an IVI of 40.61 (Tables 5 and 6). Accordingly, Juniperus procera, Hagenia abyssinica, and Syzygium guineense were identified as generalists that abundantly occurred in both edge and interior habitats, whereas Hypericum revolutum Vahl. was identified as marginalized species that characteristically dominated the edge habitats [62, 63]. However, no woody species was found as a specialist that typically occurred in the interior habitats.

3.4. Effects of Landscape Change in Floristic Composition and Structure

The computed regression analysis among the landscape indices and species composition and structural parameters in this study revealed that only PN and AREA_MN significantly affected both the species composition and structural properties of the study area. Accordingly, PN was strong and negatively affected the overall species richness (r = −0.90, ) and Shannon diversity index (r = −0.96, ) (Table 7). Conversely, the overall species richness (r = 0.95, ) and Shannon diversity (r = 0.87, ) have shown strong and positive correlation with AREA_MN. This implies that as the number of fragmented habitats increases, species richness and diversity, particularly interior-dependent species, decreases. However, edge-dependent species comfortably flourished. One of the consequences of habitat fragmentation is an increase in the proportional abundance of the edge influenced habitat and its adverse impacts on interior-sensitive species [64]. Undoubtedly, while some species (e.g., habitat specialists) suffer from fragmentation, others benefit from it (e.g., generalists and edge species) [65]. Consequently, PN was strong and negatively correlated with AREA_MN (r = −0.71, ). This implies that as the PN increases, the area of fragments decreases; as a result, small fragments contain a smaller species richness and lower species density than large fragments [66]. Large areas of habitat tend to support more individuals and, hence, more species [67].

Besides, modifying the spatial pattern of the landscape habitat size reduction and increase in isolation cause an alteration in the dispersal rate, affecting survival and mortality of individuals [8]. Many population and community changes in habitat fragments were commonly attributed to edge effects [66]. Interior species may be affected by the size decrease in their habitat, by edge effect, and by competition with generalists [68, 69]. The most threatened endemic species due to edge effect in the BMNP were Helichrysum harennense Mesfin, Kniphofia insignis Rendle, Rubus erlangeri Engl., and Vepris dainellii Pichi Serm Kokwaro. Also, the most common invasive species in the study area favored by edge effect was Achyranthes aspera L, which is also common in the disturbed forests and forest edges of the dry Afromontane forests and moist Afromontane forests in Ethiopia [70]. The gradual decline in the more sensitive species may induce a species turnover in fragments and cascade effects [62, 63].

Among the landscape indices computed, only PN and AREA_MN significantly affected some of the floristic structural properties assessed. Thus, the PN was strong and negatively affected the woody species density (r = −0.84, ) and basal area (r = −0.96, ), as well as AREA_MN was strong and positively affected the density (r = 0.71, ) and basal area (r = 0.82, ) of woody species. Habitat destruction, isolation, and transformation affect the structure and dynamics of populations, communities, and ecosystems, as well as ecological processes [71]. Generally, as AREA_MN and COA of patches increase, species richness, diversity, evenness, woody species density, basal area, DBH, and height also increase. However, as PN, SHAPE_MN, ED, ENN_MN, and IJI of patches increase, floristic composition and structural variables decrease. This implies that the landscape composition and configuration change may potentially affect the vegetation composition and structure of a particular area.

4. Conclusion

This study recognized that the Bale Mountains National Park has a diverse biodiversity and is an ecologically significant area. It contains a variety of life forms with good ecological integration. It also harbors a number of endemic floras and faunas. However, currently, anthropogenic disturbances strongly impaired the plant species composition and structure as well as the overall ecological integrity of the landscape. The progressive settlement and agricultural land expansion at the expense of natural forest and grassland coupled with human-induced recurrent fire and livestock grazing in park became a potential threat to the landscape structure. This was due to the escalated human and livestock population and their corresponding demands and necessities increment in the park. Both the floristic composition and structure were affected by the expansion of edge habitat and shrinkage of interior habitat. Species richness and diversity were higher in the edge habitat, whereas density, frequency, and basal area were higher in the interior habitat. Therefore, maintenance of the habitats heterogeneity in the park is essential for long-term population persistence. Moreover, human activities in the park should be banned and settlements in the park should be relocated to other areas to avoid their potential impacts on the floras and faunas. Finally, studies on microenvironmental factors such as light availability, air and soil temperature, humidity, and soil nutrients along the edge and interior gradient should be conducted to determine their effect on species richness, composition, and structure.

Data Availability

The data used to support the findings of this study are included within this paper.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors would like to express their special thanks to the Ethiopian Wildlife Authority for giving permission to conduct this study on-site.