https://orcid.org/0000-0003-0291-3666
Figures
Abstract
Plant community types are influenced by topographic factors, the physical and chemical properties of soil. Therefore, the study was carried out to investigate the relationships of soil and topographic factors on the distribution of species and plant community formation of the Dega Damot district in Northwestern Ethiopia. Vegetation and environmental data were collected from 86 plots (900 m2). Agglomerative hierarchical cluster analysis and redundancy analysis (RDA) with R software were used to identify plant communities and analyze the relationship between plant community types and environmental variables. Five plant community types were identified: Erica arborea-Osyris quadripartita, Discopodium penninervium-Echinops pappii, Olea europaea -Scolopia theifolia, Euphorbia abyssinica-Prunus africana, Dodonaea anguistifolia-Acokanthera schimperi. The RDA result showed that the variation of species distribution and plant community formation were significantly related to altitude, organic matter, moisture content, slope, sand, pH, EC, total nitrogen and phosphorus. Our results suggest that the variation of plant communities (Community 1, 2, 3, and 4) were closely related to environmental factors, including altitude, moisture content, OM, slope, sand, pH, EC, soil nitrogen, and phosphorus, among which altitude was the most important one. However, all the measured environmental variables are not correlated to Dodonaea anguistifolia-Acokanthera schimperi community type. Therefore, it can be concluded that some other environmental variables may influence the species composition, which is needed to be further investigated.
Citation: Birhanu L, Bekele T, Tesfaw B, Demissew S (2021) Relationships between topographic factors, soil and plant communities in a dry Afromontane forest patches of Northwestern Ethiopia. PLoS ONE 16(3): e0247966. https://doi.org/10.1371/journal.pone.0247966
Editor: Mehdi Heydari, Ilam University, ISLAMIC REPUBLIC OF IRAN
Received: October 6, 2020; Accepted: February 16, 2021; Published: March 12, 2021
Copyright: © 2021 Birhanu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All data are available with in the paper.
Funding: This research was funded by Addis Ababa University. The funders had no role in the study design, data collection and analysis, decision to publish or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The existence of plant communities is due to the interaction between plant species and their environment [1]. Plant community distribution pattern is influenced by many environmental factors such as climate, soil nutrients and topographic features [2–5]. Consequently, the relationship between the distribution of plant communities and environmental factors is one of the most important research problems in plant ecology [6, 7].
At regional and global scales, plant species responses are related to climatic factors [8, 9], while at local scales, topographic and edaphic factors play critical roles in controlling plant community formation [10–12]. According to [13] and [14], soil and topography are known as important factors affecting vegetation distribution. Therefore, environmental variables are important not only in identifying plant community structure and species distribution variations at a spatial scale but also in providing insight into the environmental requirements of the plant species needed for successful ecological restoration and biodiversity protection [15, 16].
In Ethiopia, different plant ecological studies have been conducted on the distribution of plant community types associated with environmental factors [e.g. 17–23]. However, no study is known about the relationships between environmental factors and plant species distribution and community formation in the Dega Damot district of Northwestern Ethiopia. This lack of knowledge hampers targeted conservation and restoration efforts [24]. Therefore, the objective of this study was to study the classification of plant communities and assess the relationship between plant community types and environmental variables in the study area.
Materials and methods
Description of the study area
The study was carried out in a dry Afromontane forest landscape in Dega Damot district, Amhara Regional State, in the Northwestern part of Ethiopia. The major nearby town is Feresbet which is 399 km away Northwestern of Addis Ababa. The district is bordered by Dembecha (in the south), JabiTehnan (in the southwest), Quarite (in the west), HuletIju Enesie (in the north) and Bibugn (in the east) districts. The study area extends between 10°44′41.69″N and 11°00′10.56″N, and 37°26′58.43″E and 37°42′35.42″E (Fig 1). Degadamot district covers approximately 83,124 hectares of land area [25]. The district has forest patches. Of them, five forest patches were selected for the study which are relatively undisturbed in the district. Topographically the study area consists of mountainous areas, undulating, valleys, and plains [25, 26]. The elevation of the study area ranges from 1738 to 3586 m above sea level. On the other hand, the slope classes range from 0–50 degrees [25].
Our study does not need to supply a copy right notice for Fig 1. Because the location map of the study area (Fig 1) shape file data was obtained from Ethiopian Mapping Agency (https://africaopendata.org/dataset/ethiopia-shapefiles) is free and open to researchers.
The average monthly maximum temperatures were 24.7 0C and monthly average minimum temperatures 9.2 0C. Annual rainfall and temperatures were 2113 mm and 15.9 0C, respectively [25]. According to the data obtained from the Geographic Information System (GIS) department of the Ministry of Water, Irrigation and Electricity (MoWIE; [27], the dominant soil type of the study area are Luvisols (28523.31ha) and Cambisols (17552.81 ha) (See Fig 2 and Table 1). The vegetation of the study area belongs to the categories of dry evergreen Afromontane vegetation [28]. The area is covered by agricultural lands, settlements, forestland, shrubland bare lands and grazing lands. Agriculture is common practice in the Dega Damot district, with typical crops of Hordeum vulgare (barley), Triticum sp.(wheat), Vicia faba (faba bean), Eragrostis tef (teff), Zea mays (maize) and Solanum tuberosum (potato) [25, 26]. We obtained written permission from Agriculture and Rural Development Office of the Dega Damot district to collect both vegetation and environmental data and conduct our research.
Because the shape file data for soil map obtained from the Geographic Information System (GIS) department of the Ministry of Water, Irrigation and Electricity (MoWIE), is free and open to researchers.
Reconnaissance survey and sampling technique
A reconnaissance survey was made in January 2017 across the forest patches to get an impression of the site conditions and identify the sampling sites in the study area. The forest patches were located at different altitudes, ranging from 1881–2947 m and also occupy different sizes ranging from 50 ha to 405 ha (Table 2). The actual fieldwork was conducted between February and November 2017. A systematic sampling technique was used for vegetation and environmental data collection. Sampling sites were arranged along transects from the forest patches. The first transect and plot were selected purposely at one side of the forest by avoiding the forest edge.
A total of 22 line transects (6 each in Shangi Derke and Bunsie forest patches, 4 in Aradie forest patch and 3 each in Masikana and Samana forest patches) were laid along elevation gradients. The number of plots per transect differs depending on the length of the transect (Table 2).
The elevation between two consecutive plots and transect were 50 and 200 ma.s.l. apart in the forests, respectively. Accordingly, a total of 86 sampling plots of each measuring 900m2 (30 m × 30 m) were taken from all 5 forest patches. Of which, 29 from Shangi Derke, 32 from Bunsie, 10 from Aradie, 9 from Masikana Medihaniyalem and 6 from Samana forest patches were sampled. When a study mainly focuses on only woody species composition and diversity, research has mainly ignored studies on diversity and not included herbaceous species and rare species which are important to the understanding of a comprehensive list of plant species diversity in the study area. Therefore we have included the herbaceous species. Five smaller subplots of 2 *2 m2, four at the corner, one at the center of the main plot were established for herbaceous plant species collection. In our study the overall cover of herbaceous species estimated by the average of individual cover values of species taken from the subplots (4 m2).
Vegetation data collection
In each plot, all individual trees and shrubs with a diameter at breast height (DBH at 1.30 m above ground) >2.5 cm were measured for DBH classes using a measuring tape. The height of woody individual species with > 3m height was measured using Suunto Clinometers. Within each plot, the cover/abundance values were estimated for all plant species. Three-column data tables were constructed in Microsoft Excel 2007 and saved in the CSV (comma delimited) format. The first, second, and third columns represented by plots, species and abundance of species, respectively. After that, the percent cover of each species was transformed to ordinal scales and assigned to one of the nine cover-abundance classes according to the modified 1–9 Braun-Blanquet scale [29]. The scales for cover-abundance values are 1 ≤ 0.1%, 2 = 0.1 to 1%, 3 = 1 to 2%, 4 = 2 to 5%, 5 = 5 to 10%, 6 = 10 to 25%, 7 = 25 to 50%, 8 = 50 to 75%, and 9 >75%. Finally, the three-column data of the cover-abundance value of each species were imported to R statistical software version 3.5.2 [30] and matrified to carry out the cluster analysis and other parameters such as for ordination, diversity and evenness. The packages labdsv” and vegan” were used to transform the cover values of the vegetation data into the modified Braun Blanquet 1–9 scale and matrified the abundance data.
Moreover, species occurring outside the plots were also recorded to obtain a complete floristic list of the study area. Voucher specimens were collected, coded, pressed and dried for subsequent identification and verification at the National Herbarium (ETH), Addis Ababa University, using Volumes 1–8 of Flora of Ethiopia and Eritrea.
Environmental data collection
Geographical data (altitude, latitude, and longitude) were recorded using GPS for each plot in the forest patches. The slope of each sample plot was measured using Suunto clinometer.
Soil samples were collected from these 86 plots for soil analysis. The samples were collected from topsoil (0–30 cm). From each plot, the soil samples were collected at five locations, one from the center and four from corners and these are mixed to produce a composite sample and only one sample was taken for analysis. The samples were air-dried, grounded, and passed through a 2-mm sieve to remove the stone pieces and large root particles before analysis. The soil analysis was assessed in a soil laboratory of Debremarkos soil research and fertility improvement center.
Soil texture, soil moisture content, bulk density, pH, EC, CEC, total nitrogen, organic matter and available phosphorus were measured for the soil data collected from the selected forest patches. The soil bulk density (BD) was measured using core method prescribed for undisturbed soils. It was was calculated as the ratio of oven-dried soil weight to the volume of the soil:
(1)
where wd is the weight of oven-dry soil (g) and v is the volume of the soil (cm3). The volume was calculated from the volume of the core sampler. The core sample was oven-dried to a constant weight using an oven at 105°C for 24 h. The soil texture was measured using the hydrometer method following [31]. The apparatuses used in measuring soil texture were plunger, hydrometer jar and hydrometer. Soil pH was measured by taking a 1:2.5 soil/water suspension using a glass electrode pH-meter following [32]. The electrical conductivity of the soil sample was measured by preparing a 1:2.5soil/water suspension using a conductivity-meter following [33]. Organic carbon was determined by Walkely and Black method [34]. The organic matter (OM) content was calculated by multiplying the OC content by 1.724 (OM = 1.724 x OC). Total nitrogen was ascertained by the Kjeldhal procedure after digestion with concentrated H2SO4 [35]. Available phosphorus (P) was analyzed by Uv/vis spectrophotometer following the Olson method [36]. Cation exchange capacity (CEC) was determined after extracting the soil samples by ammonium acetate (1N NH4OAc) at pH 7.0 [37].
Data analysis
Plant community classification.
Cluster analysis helps to group a set of observations (plots) based on their floristic similarities [38, 39], i.e. to determine plots that can be classified into the same groups based on the species abundance data. Therefore, for this analysis the data matrix contained 86 plots and 164 plant species were grouped using hierarchical cluster analysis. In this study, an agglomerative hierarchical clustering (AHC) was employed using a similarity ratio (SR) and Ward method (Minimum-variance clustering) to identify plant communities. R statistical software version 3.5.2 [30] and the packages Cluster and Vegan [40, 41] were used during cluster analysis.
An indicator species analysis (ISA) was performed using the indicator value (IndVal) method in R package cluster, Vegan and Labdsv [41, 42] to identify diagnostic species for the naming of plant communities. The indicator value index (IndVal) is based on the abundance of a given species and its occurrence within a given set of samples. The statistical significance of the indicator values for each species was assessed by a Monte Carlo (randomization) test procedure. In this analysis, species with a significant indicator value of P<0.05 is considered to be an indicator species of a community. Finally, the plant community types were named after two of the dominant species that had an indicator value of p< 0.05.
Species diversity.
Shannon-Wiener diversity indices and Shannon’s evenness were computed to describe species diversity, richness and evenness of the study area, sites and plant communities [39].
Where H’ = Shannon diversity index
s = number of species,
pi = proportion of ith species, ln = the natural logarithm.
Shannon’s evenness index (J) was also calculated by using:
(3)
where H’ = Shannon–Wiener Diversity Index; and H’ max = lns where s is the number of species in the plot.
Sørensen’s similarity coefficient was used for comparing the similarity of two communities or forest patches.
(4)
where b and c are the no_of species in communities, forest patches or study areas b and c; a -
the number of species common to both and Ss is the Sørensen’s similarity coefficient.
Ordination.
Ordination is a multivariate method that articulates the relationships between species, plotes and environmental variables in a low-dimensional space using ordination diagrams [38, 43]. A preliminary analysis of the data by Pearson’s correlation coefficient was calculated to find a significant correlation between environmental variables and to remove auto-correlated variables from the ordination analysis using SAS software.
Redundancy Analysis (RDA) ordination method was used for this analysis for testing the vegetation and environmental variables relationships [41, 44]. RDA was selected due to its better visualization of the graph. In this analysis, 86 sample plots, 163 species and 10 environmental variables (altitude, slope, bulk density, moisture content, clay, pH, EC, available Phosphorus, CEC, and OM) were included. Adonis test was performed to determine the significance of the environmental variables on plant community distribution. Consequently, RDA ordination was plotted using cover-abundance values of plant species and data of significant environmental variables.
One-way ANOVA followed by post-hoc Tukey HSD test was used whether there were significant mean differences among plant communities about environmental variables, species richness, diversity and evenness.
Results
Floristic composition
A total of 176 species, belonging to 80 families were recorded in the study area. Asteraceae and Fabaceae were the top most dominant families. Of the total plant species composition of the forest, 18 (10.6%) species are endemic to Ethiopia and Eritrea.
Plant community types
Five plant community types were identified from the hierarchical cluster analysis in the study area (Fig 3). The community is named after one or two dominant indicator tree or shrub species selected by the relative magnitude of their indicator values. In this study, a species is considered as an indicator of a group when its indicator value is significantly higher at p < 0.05 (Table 3). Consequently, the identified communities were Erica arborea-Osyris quadripartita, Discopodium penninervium-Echinops pappii, Olea europaea subsp. cuspidata-Scolopia theifolia, Euphorbia abyssinica-Prunus africana, Dodonaea anguistifolia-Acokanthera schimperi community types.
Erica arborea-Osyris quadripartita community type (C1)
This community is found in the Shangi Derkie forest patch, which has a very steep slope. The altitudinal range of this community type lies between 2516 and 2947 ma.s.l. The community has six indicator species with significant indicator values, namely; Erica arborea, Osyris quadripartita, Bersama abyssinica, Buddleja polystachya, Maesa lanceolata, and Rosa abyssinica. The tree layer was dominated by Ekebergia capensis and Lepidotrichilia volkensii. The shrub layer was dominated by Lippia adoensis. The most dominant species in the herb layer include Kniphofia foliosa, Kalanchoe petitiana, Verbascum sinaiticum and Trifolium decorum. This community type is represented by 108 species, being the richest in a number of species among the five community types.
Discopodium penninervium-Echinops pappii community type (C2)
The altitudinal range of this community type lies between 2375 and 2916 m a.s.l and located in Shangi Derkie, Aradie, Bunise, and Masikana Medhanyalm forest patches. The indicator species are Discopodium penninervium, Echinops pappii and Laggera tomentosa. The most dominant shrub species were Solanecio gigas, Maytenus arbutifolia, Solanum incanum, Rumex nervosus, Gnidia glauca and Acanthus sennii. The tree layer was dominated by Arundinaria alpina. Liana such species such as Asparagus africanus and Phytolacca dodecandra were dominant. The herb layer was dominated by Plectranthus assurgens and Carduus schimperi. This community is comprised of 104 species.
Olea europaea -Scolopia theifolia Community type (C3)
This community lies between altitudes 2475 and 2906 m above sea level. Similar to community 1 and 2, Olea europaea subsp. cuspidata—Scolopia theifolia Community type is also found in Bunise, shangi Derkie, Aradie, and Masikana Medihanyalm forest patches. The indicator tree and shrub species of this community are Olea europaea subsp. cuspidata and Scolopia theifolia, Olinia rochetiana, Myrsine africana, Clutia abyssinica, Pittosporum viridiflorum. The dominant species in the tree layer include Apodytes dimidiata. The shrub species such as Rhus glutinosa, subsp. glutinosa, Dovyalis verrucosa, and Catha edulis were highly dominant. The herb layer comprised of Oplismenus hirtellus, Impatiens hochstetteri, Impatiens rothii and Justicia heterocarpa. This community is comprised of 73 species.
Euphorbia abyssinica-Prunus africana community type (C4)
This community type is located between altitudes 2455 and 2733 m asl, and which is found in 3 forest patches: Bunise, Shangi Derkie, Aradie, and Masikana. The indicator species are Euphorbia abyssinica and Prunus africana. The shrub layer includes Solanecio gigas and liana species Embelia schimperi were also dominant. The herb layer was dominated by Hypoestes forskaolii. It consists of 84 species.
Dodonaea anguistifolia-Acokanthera schimperi community type (C5)
Community 5 found in the Samana forest patch. This community is established at lower elevations (1881–2110 m a.s.l.). Dodonaea anguistifolia and Acokanthera schimperi are the characteristic species of the shrub, and tree layers, respectively. Other dominant tree species include Croton macrostachyus, Acacia seyal, and Combretum molle. The shrub layer was dominated by Euclea schimperi. Dominant species of the herb layer were Bidens macroptera and Satureja abyssinica.
Diversity and evenness of plant communities
Shannon`s diversity indices showed that community 1 had the highest species diversity and richness followed by community 2, while community 5 had the lowest species diversity and richness. On the other hand, community 4 had the highest evenness value followed by communities 3 and 5. Community 1 and 2 had the least evenness value (Table 4).
Similarities between plant communities
Results of plant community similarity indicated that communities 2 and 3, communities 1 and 2, communities 1 and 3 have the same similarity index values. The least similarity was between community 1 and 5, community 3 and 5, community 4 and 5 (Table 5).
Correlation of environmental variables
Pearson’s correlation matrix of the environmental factors is shown in Table 6. Altitude was positively correlated with moisture content, nitrogen and sand are negatively correlated with pH and bulk density. Nitrogen and phosphorus are positively correlated with soil moisture content and organic matter. CEC is positively correlated with organic matter, nitrogen and phosphorus. Soil pH shows a positive correlation with phosphorus while EC shows a positive correlation with organic matter (see Table 6).
Relationship between plant community types and environment variables
RDA ordination.
Out of the 10 environmental variables included in ordination analysis (slope, bulk density, moisture content, elevation, clay, pH, P, EC, OM and CEC), only elevation, moisture content, OM, slope, pH, EC, clay and phosphorus had significant (p< 0.05) correlated on species compositions of communities and their distributions as shown in Table 7 and Fig 4. We described several of these variables are significant but none of them account for much variability except elevation. A high sum of squeres would indicate a lot of variability in the data, while a low sum of squeres (most of the measurement close to mean) would indicate a low amount of variability (Table 7). The RDA diagram indicated that the first axis was primarily correlated with elevation, clay, and moisture content, while the second axis was correlated with P, OM and moisture content (Fig 4, Table 8). The eigenvalue for the first and second axis were 0.54 and 0.19, respectively. The first axis eigenvalue which closer to 1 and show distribution of the plant communities along the axis is good. The cumulative proportion variance explained by the first six RDA axis of the joint biplot was 92%. The variation in patterns of plant species distribution and plant community formation for the first and second axis were 46%, and 16%, respectively. This showed that 62% of the variation in patterns of plant species distribution and plant community formation was explained by axis one and two.
The arrows in the diagram stand for the environmental factors, the length of each arrow indicates the contribution of the factor to ordination axes, The numbers refer to quadrat number, and the angle between the arrows and he axes indicates the correlation between the variable and the ordination axe.
Consequently, the first two axes are sufficient to reflect the relationship between species and environmental factors. In the first RDA axis, elevation was the most important variable and separated communities at high elevation (communities 1, 2, 3, and 4) from communities distributed at lower elevation (communities 5) on the ordination axis (Fig 4). On the other hand, in the second axis, communities 2 and 3 were located on sites with higher organic matter and moisture content. The values of the different environmental variables average for each plant community type are presented in Table 9.
Discussion
We found relatively high plant species richness in Dega Damot district forest patches. However, the study area is fewer species-rich, especially in woody species than similar studies conducted in the Afromontane forest in Ethiopia, e.g. Tara Gedam and Abebaye forests [45], Afromontane forest patches of Awi zone [46] It has higher species richness than other dry Afromontane forests in Ethiopia such as the Zengena forest [47], Kuandisha forest [48], and Amoro forest [26]. Besides, the number of woody species recorded in this study is more or less similar to that of the Zege Peninsula forest [49] whereas [50] recorded low woody species in the vegetation of Kalfou forest reserve, Cameroon. The difference in species composition among the different areas is due to the number of plots sampled and its size can somehow explain this heterogeneity of the species richness. [51] reported that forests with a high degree of human interference and disturbances show relatively lower species richness than others. Thus, the present study suggests that the Dega Damot district had relatively high species composition, compared to other similar vegetation types in Ethiopia.
The dominance of Asteraceae was reported from other studies in dry evergreen Afromontane forests vegetation type [23, 52, 53]. [54] also reported that Asteraceae was the richest family in terms of species numbers for the flora of Ethiopia and Eretria. The dominance of Asteraceae in the present study may indicate that the forest patches might have been under a certain level of disturbances. As reported by [55] and [54], Asteraceae usually have a preference for open and disturbed lands to grow.
According to [56] the proportion of endemic plant species in the montane forests of Ethiopia is high, between 11–15% of the total number of species. Similarly, the results from the study area showed that relatively high endemicity (10.6%). Ethiopia is one of the centers of plant endemism in East Africa [57]. This study is in line with the typical feature of Afromontane forests that house numerous endemic species [58]. However, the Ethiopian Afromontane forests are one of the most degraded forests and continuously shrinking mainly due to anthropogenic disturbances [59, 60].
Plant communities
Based on the hierarchical cluster analysis result, 5 plant communities were identified in the study area. However, the location of the 3 community types (community 2, 3, and 4) had overlapping ranges of elevation. Elevation represents a complex gradient combination of many different environmental factors such as topography, soil, moisture and climate [59] So that it is difficult to separate other environmental factors [61].
The vegetation of the area is classified as dry evergreen Afromontane forest and grassland complex (DAF) [28] and dominated by the species Olea europaea subsp. cuspidata, Juniperus procera, Discopodium penninervium, Pittosporum viridiflorum, Erica arborea, Hypericum revolutum, Dodonaea angustifolia, Acacia abyssinica, Bersama abyssinica, and Prunus africana. The Olea europaea subsp. cuspidata, indicator species was also mentioned to form the upper canopies of different dry Afromontane forests in Ethiopia [19, 21].
The identified plant communities are characterized by their different floristic composition. This could be attributed to variations in environmental factors. According to [62, 63], the variation species composition among plant communities probably associated with the effects of environmental factors. Nevertheless, species diversity and richness were not the same among plant communities. For instance, the highest species richness and diversity were recorded in community 1. On the contrary, community 5 had the least species richness and diversity than the remaining community types. The reason for high species diversity and richness of community 1 may be due to the highest altitudinal ranges in which this community is found (2516–2947 m). Also, community 5 type is the most disturbed community due to settlement, agricultural expansion, overgrazing by livestock and crossing road.
Anthropogenic land-use change is one of the most important factors contributing to global change [64], and changes in landscape structure and agricultural land-use intensity are likely to influence plant community composition and species richness. Changes in land-use practice may result in continuously or more abruptly deteriorating environmental conditions for some plant species, causing their decline in abundance and distribution [65]. For instance, over the last 31 years, the forest cover of Dega Damot district has decreased while areas under farmland, grazing, and settlements have increased [25]. Therefore, these anthropogenic land-use changes might influence the plant community formation of the present study in a different plant species richness and diversity. The result of the pairwise comparison of Sorensen’s similarity coefficient index in species composition between the 3 plant communities (Community 1, 2 and 3) showed high similarity. These due to the location of this community have relatively similar environmental factors (soil and altitudinal range). Community 5 located at the lowest elevation with low content of organic matter, soil moisture content and which may have a less floristic similarity from the remaining plant community types (Community 1, 2 and 3).
Environmental factors and plant community relationship
The distribution of plant communities in the study area reflects the combined influence of altitude and soil factors. Based on the results of RDA, elevation was the major environmental variable in explaining variations in plant species distribution and patterns of plant community formation although there are overlaps among some community types. This might be associated with a continuous change in environmental variables along the elevational gradient [63, 66]. Other studies conducted in Ethiopia [see 18, 23, 67] also noted that altitude is the most important environmental variable for the determinant of vegetation variation. Moreover [68], reported that vegetation distribution influenced by elevation. Besides, elevation is a vital environmental factor that affects the atmospheric pressure, moisture, and temperature which have a strong influence on the growth and development of plants and the distribution of vegetation [69].
Moreover, organic matter was also the most important constraining variable in plant community formation of the study area. Research in the Jibat forest, Ethiopia, by [19] reported that plant distributions were affected by the organic matter at a higher elevation. SOM, its role in soil structure and moisture retention capability are well known [70, 71]; these effects may account for its relevance for determining plant species distribution and community formation in the study area, besides its role as a source of soil nutrients to plants. The mineralization of organic matter is a contributing factor in supplying available mineral nutrients for plant use and through decomposition process the available nitrogen.
Similar to the present findings, nitrogen and P were also attributed to have significant effects on species compositions of plant communities in remnant Afromontane forests on the central plateau of Shewa [20]. Phosphorus is an essential nutrient for the growth of plant. This result is also in agreement with past similar works in in Afromontane and transitional rainforest vegetation of southwestern Ethiopia [67]. Furthermore, a study conducted in Brazil by [72], P was the main factor that distinct plant communities.
Topographic features associated with soil properties are strongly correlated with species distribution and plant community structure on a local scale [33]. For instance, the Dodonaea anguistifolia-Acokanthera schimperi community, with a low content of organic matter, soil moisture content, and N, is differentiated from the other community types in the study area. The anthropogenic disturbances are higher at the lower altitudes (selective cutting of trees, grazing by livestock and expanding of farmlands) and the temperature is also higher at these altitudes. Generally, in the study area, the soil organic matter is higher at higher altitudes. This might be due to the decreased decomposition of organic matter and the long-term accumulation of organic matter [73]. Much of the soil nitrogen is obtained from organic matter [74] and hence is bound in higher quantities at higher elevations.
Conclusion
Five plant communities were identified from this study. Among all of the investigated environmental factors, elevation, moisture content, slope, pH, EC, clay, OM, and phosphorus were found to significantly explain variation in species composition and community formation in the study area. Elevation was the most important environmental factor influencing species distribution and community formation. The lower elevation resulted in a decline in species richness, diversity, organic matter, moisture content, especially in the community 5 in the forest. The anthropogenic disturbances present in this community type may also contribute to the low species richness and diversity recorded. Also, Dodonaea anguistifolia-Acokanthera schimperi community type is not correlated to all the measured environmental variables. Therefore, it can be concluded that other environmental factors may influence the plant community formation, which is required to be additionally studied.
Supporting information
S1 Appendix. Floristic list of Dega Damot district.
https://doi.org/10.1371/journal.pone.0247966.s001
(DOCX)
S2 Appendix. Descriptive statistical analysis of the topography and soil data.
https://doi.org/10.1371/journal.pone.0247966.s002
(DOCX)
Acknowledgments
The authors express their deepest thanks to Agriculture and Rural Development Office of the Dega Damot District as well as to the respective chairpersons for their assistance for field data collection. We extend our sincere thanks to the anonymous reviewers of PLOS ONE for their comments.
References
- 1.
Chapin FS, Matson PA, Mooney HA. Principles of terrestrial ecosystem ecology. Springer, New York. 2002.
- 2. Abella SR, Covington WW. Vegetation–environment relationships and ecological species groups of an Arizona Pinus ponderosa landscape, USA. Plant Ecology. 2006; 185 (2):255–68.
- 3. Eshaghi R. J., Banj shafiei A. The Distribution of Ecological Species Groups in Fagetum. Communities of Caspian Forests: determination of Effective Environmental Factors. Flora. 2010; 205:721–727.
- 4. Liu S L, Ma KM, Fu B J. The relationship between landform, soil characteristics and plant community structure in the Donglingshan Mountain Region, Beijing. Acta Phytoecologica Sinica. 2003; 27(4): 496–502.
- 5. Pugnaire FI, Armas C, Maestre FT. Positive plant interactions in the Iberian Southeast: mechanisms, environmental gradients, and ecosystem function. Journal of Arid Environment. 2011; 75:1310–1320.
- 6. Burke A. Classification and ordination of plant communities of the Naukluft Mountains, Namibia. Journal of Vegetation Science.2001; 12: 53–60.
- 7. Yavitt JB, Harm KE, Garcia MN Spatial heterogeneity of soil chemical properties in a lowland tropical moist forest, Panama. Soil Research. 2009; 47:674–687.
- 8. Liu Y., Zhang Y., He D., Cao M., Zhu H. Climatic control of plant speciesrichness along elevation gradients in the Longitudinal Range-Gorge Region. Chinese Science Bulletin. 2007; 52:50–58.
- 9. Jarema SI, Samson J, McGill BJ. Variation in abundance across a species’ range predictsclimate change responses in the range interior will exceed those at the edge: a case study with North American beaver. Global Change Biology. 2009; 15:508–522.
- 10. Wijesinghe D. K., John E. A., and Hutchings M. J. Does the pattern of soil resource heterogeneity determine plant community structure? An experimental investigation. Journal of Ecology. 2005; 93: 99–112.
- 11. Katabuchi M., Kurokawa H., Davies S.J., Tan S., Nakashizuka T. Soil resource availability shapes community trait structure in species-rich dipterocarp forest. Journal of Ecology. 2012; 100: 643–651.
- 12. Fadrique B, Homeier J. Elevation and topography influence community structure, biomass and host tree interactions of lianas in tropical montane forests of southern Ecuador. Journal of Vegetation Science. 2016; 27 (5):958–68.
- 13. Poulos H. M., Taylor A. H. and Beaty R. M. Environmental controls on dominance and diversity of woody plant species in a Madrean, Sky Island ecosystem, Arizona, USA. Plant Ecology. 2007; 193: 15–30.
- 14. Solon J., Dego´ rski M. and Roo-Zielinska E. Vegetation response to a topographical-soil gradient. Catena. 2007; 71: 309–320.
- 15. Khurana E., Singh JS. Ecology of seed and seedling growth for conservation and restoration of the tropical dry forest: a review. Environ Conserv. 2001; 28:39–52.
- 16. Toledo M., Peña-Claros M., Bongers F., Alarcón A., Balcázar J., Chuviña J. et al. Distribution patterns of tropical woody species in response to climatic and edaphic gradients. Journal of Ecology. 2012; 100:253–263.
- 17. Gebre Egiziabher TB. Vegetation and environment of the mountains of Ethiopia: Implications for utilization and conservation. Mountain Research and Development. 1988; 8: 211–216.
- 18.
Woldu Z, Feoli E, Nigatu L. Partitioning an elevation gradient of vegetation from southeastern Ethiopia by probabilistic methods. In Numerical syntaxonomy (pp. 189–198). Springer, Dordrecht. 1989.
- 19. Bekele T. Vegetation ecology of remnant Afromontane forests on the central plateau of Shewa, Ethiopia. Sv. växtgeografiska sällsk. 1993.
- 20. Bekele T. Phytosociology and ecology of a humid Afromontane forest on the central plateau of Ethiopia. Journal of Vegetation Science. 1994; 5 (1):87–98.
- 21. Teketay D, Bekele T. Floristic composition of Wof‐Washa natural forest, Central Ethiopia: Implications for the conservation of biodiversity. Feddes Repertorium. 1995; 106 (1‐5):127–47.
- 22. Wale HA, Bekele T, Dalle G. Plant community and ecological analysis of woodland vegetation in Metema Area, Amhara National Regional State, Northwestern Ethiopia. Journal of Forestry Research. 2012; 23(4): 599−607.
- 23. Yineger H, Kelbessa E, Bekele T, Lulekal E. Floristic composition and structure of the dry Afromontane forest at Bale Mountains National Park, Ethiopia. SINET: Ethiopian Journal of Science. 2008; 31(2):103–20.
- 24. Sutherland W.J., Pullin A.S., Dolman P.M., Knight T.M. The need for evidence-based conservation. Trends in Ecology and Evolution. 2004; 19: 305–308. pmid:16701275
- 25. Birhanu L, Tesfaw BH, Bekele T, Demissew S. Land use/land cover change along elevation and slope gradient in highlands of Ethiopia. Remote Sensing Applications: Society and Environment. 2019; 16: https://doi.org/10.1016/j.rsase.2019.100260
- 26. Birhanu L, Bekele T, Demissew S. Woody species composition and structure of Amoro Forest in West Gojjam Zone, Northwestern Ethiopia. Journal of Ecology and the Natural Environment. 2018; 10 (4):53–64.
- 27.
FAO. Ethiopian highland reclamation study, Ethiopia. Final report, FAO, Rome. 1986.
- 28.
Friis I., Demissew S. and Breugel P.V. Atlas of the Potential Vegetation of Ethiopia. The Royal Danish Academy of Sciences and Letters, Copenhagen, Denmark. 2010.
- 29. Van der Maarel E. Transformation of cover-abundance values in phytosociology andits effects on community similarity. Vegetation. 1979; 39 (2):97–114.
- 30. R Development Core Team. R: A Language and Environment for Statistical Computing. A Programming Environment for Data Analysis and Graphics. Version 3.5.2. 2018.
- 31. Khurana E., Singh JS. Ecology of seed and seedling growth for conservation and restoration of the tropical dry forest: a review. Environ Conserv. 2001; 28:39–52.
- 32.
Houba V.J.G, Van der Lee J.J, Nouozamssky I and Walinga. I. Soil-plant analysis a serious of syllabi. part 5 soil analysis procedures. Wagengien agricultural university. the Netherlands. 1989. https://doi.org/10.1016/0378-1097(89)90224-3 pmid:2515104
- 33.
Chopra S.L, and KallWar J.S. Analytical agricultural chemistry. Kalyani Publishers, New Delhi. 1982.
- 34.
Nelson D.W., Sommer L.E. Total organic carbon and organic matter. Methods of soil analysis part 2: chemical and microbiological properties. 539-579pp. second edision. A.S.A Modision Wisconsin USA. 1982.
- 35. Dewan H.C., Admasu N. Physical and chemical methods for soil, water, and fertilizer 1. DP/ETH/82/011. Field document NSSP. Adiss Ababa, Ethiopia. 1987.
- 36.
Olsen S.R. and Dean L.A. Phosphorus. In: Methods of Soil Analysis, pp. 1044–1046, (Black C.A., ed.). Agronomy no. 9. American Society of Agronomy, Madison. 1965.
- 37.
Van Reeuwijk L.P. Procedure for Soil Analysis. 3rd Edition. International soil reference and information center Wageningen (ISRIC). The Netherlands. p.o.box 353. 6700 AJ Wageningen. 1992.
- 38. McCune B. and Grace J. B. Analysis of Ecological Communities. MjM Software. 2002.
- 39.
Kent M. Vegetation Description and Data Analysis: A practical approach. JohnWiley and Sons, Chichester. 2012.
- 40. Oksanen J., Blanchet F.G., Kindt R., Legendre P., Minchin P.R., O’Hara R.B., et al. Vegan: Community Ecology Package. R-Package Version 3.2.3. 2013.
- 41.
Woldu Z. Comprehensive analysis of vegetation and ecological data: Basics, concepts, and methods. Addis Ababa University Press, Addis Ababa. 2017.
- 42.
De Cáceres M. How to use the indicspecies package. Centre Tecnol`ogic Forestal de Catalunya, Spain. 2013.
- 43.
Legendre P, Legendre L. Numerical Ecology (2nd English edition). Elsevier ScienceB.V., Netherlands. 1998.
- 44. ter Braak C.J.F., Šmilauer P. Topics in constrained and unconstrained ordination. Plant Ecol. 2015; 216: 683–696.
- 45. Zegeye H, Teketay D, and Kelbessa E. Diversity and regeneration status of woody species in Tara Gedam and Abebaye Forests, north-western Ethiopia. Journal of Forestry Research 2011; 22 (3):315–328.l.
- 46.
Berhanu A. Ecology and Conservation Status of Evergreen Afromontane Forest Patches in Awi Zone of Amhara Region, Northwestern Ethiopia. Ph.D. Dissertation, Addis Ababa University. 2017.
- 47. Tadele D, Lulekal E, Damtie D, Assefa A. Floristic diversity and regeneration status of woody plants in Zengena forest, a remnant montane forest patch in northwestern Ethiopia. Journal of forestry research. 2014; 25 (2):329–336.
- 48. Berhanu A, Demissew S, Woldu Z, Didita D. Woody species composition and structure of Kuandisha Afromontane forest fragment in northwestern Ethiopia. Journal of Forestry Research. 2016; 28: 343–355.
- 49. Alelign A, Teketay D, Yemshaw Y, Sue Edwardss. Diversity and Status of Regeneration of Woody Plants on the Peninsula of Zegie, Northwestern Ethiopia. Tropical Ecology. 2007; 48 (1):37–49.
- 50. Froumsia M, Zapfack L, Mapongmetsem PM, Nkongmeneck BA Woody species composition, structure and diversity of vegetation of Kalfou Forest Reserve, Cameroon. Journal of Ecology and the Natural Environment, 2012; 4(13): 333–343.
- 51. Urban D., Miller C., Halpin P. and Stephenson N. Forest gradient response in Sierran landscapes: the physical template. Landscape Ecology. 2000; 15: 603–620
- 52. Tilahun A, Soromessa T, Kelbessa E, Dibaba A. Floristic composition and community analysis of Menagesha Amba Mariam Forest (Egdu Forest) in the central plateau of Shewa Ethiopia. Ethiopian Journal of Biological Sciences. 2011; 10 (2):111–136.
- 53.
Lulekal E. Plant diversity and ethnobotanical study of medicinal plants in Ankober District, North Shewa Zone of Amhara Region, Ethiopia. Ph.D. Dissertation, Addis Ababa University. 2014.
- 54.
Tadesse M. Asteraceae (Compositae). In: Hedberg I, Friis I, Edwards S, eds. Flora of Ethiopia and Eritrea (vol. 4–2). Addis Ababa and Uppsala: The National Herbarium. 2004.
- 55. Friis I. Forest and Forest Trees of Northeast Tropical Africa: Their natural habitats and distribution pattern in Ethiopia, Djibouti, and Somalia. Kew. Bull. Add. Ser. 15. 1992.
- 56.
Friis, I. Demissew S. Vegetation Maps of Ethiopia and Eritrea. A Review of Existing Map and the Need for the Flora of Ethiopia and Eritrea. In: Biodiversity Research in the Horn of Africa Regions (Eds I. Friis and O. Ryding). Proceedings of theThird InternationalSymposium on the Flora of Ethiopia and Eritrea atCarlsberg Academy. Royal Danish Academy of Sciences and Letters, Copenhagen. 2001.
- 57.
Vivero J.L., Kelbessa K, Demissew S. The Red List of Endemic Trees & Shrubs of Ethiopia and Eritrea. Fauna & Flora International, United Kingdom. 2005.
- 58.
White F. The vegetation of Africa. A Descriptive Memoir to Accompany theUnesco/AETFAT/UNSO Vegetation Map of Africa. UNESCO, Switzerland. 1983.
- 59. Senbeta Feyera, Schmitt Christine, Woldemariam Tadese, Hans Juergen Boehmer and Manfred Denich (2014) Plant diversity, vegetation structure and relationship between plant communities and environmental variables inthe Afromontane forests of Ethiopia. SINET: Ethiopian Journal of Science. 37; (2):113–130.
- 60. Kelbessa E, Soromessa T. Interfaces of regeneration, structure, diversity, and uses of some plant species in Bonga forest: a reservoir for the wild coffee gene pool. Ethiopian Journal of Biological Sciences. 2008; 31 (2):121–134.
- 61. Korner C. Why are there global gradients in species richness? Mountains may hold the answer. Trends Ecology and Evoution. 2000; 15:513–514.
- 62.
Whittaker R.J., Willis K.J. and Field R Climatic–energetic explanations of diversity: a macroscopic perspective. In: Macroecology: concepts and consequences, Blackburn T.M. and Gaston K.J. (eds.), Cambridge University Press; Cambridge. p107–129. 2003.
- 63. Kebede M, Yirdaw E, Luukkanen O, Lemenih M. Plant community analysis and effect of environmental factors on the diversity of woody species in the moist Afromontane forest of Wondo Genet, South Central Ethiopia. Biodiversity: Research and Conservation. 2013; 29:63–80.
- 64. Sala OE, Chapin FS, Armesto JJ, Berlow E, Bloomfeld J, Dirzo R, et al. Global biodiversity scenarios for the year 2100. Science. 2000;287 (5459):1770–177 pmid:10710299
- 65.
Naeem S, Bunker DE, Hector A, Loreau M, Andperrings C. Biodiversity, Ecosystem Functioning, and Human Wellbeing: An Ecological and Economic Perspective. Oxford: Oxford University Press. 2009.
- 66.
Lovett J.C., Wasser S.K. Biogeography and ecology of the rain forests of eastern Africa. Cambridge University Press, Cambridge, UK. 1993.
- 67. Yeshitela K, Bekele T. Plant community analysis and ecology of Afromontane and transitional rainforest vegetation of southwestern Ethiopia. SINET: Ethiopian Journal of Science 2002; 25(2):155–75.
- 68. Zhang JT, Dong Y. Factors affecting species diversity of plant communities and the restoration Process in the loess area of China. Ecological Engineering. 2010 36 (3):345–50.
- 69. Hedberg O. Features of Afro-alpine plant ecology. Acta Phytogeographica Suecica.1964; 49:1–144.
- 70. McLaren KP, McDonald MA. The effects of moisture and shade on seed germination and seedling survival in a tropical dry forest in Jamaica. Forest Ecology and Management. 2003; 183: 61–75.
- 71.
Schaetzl R. J., and Anderson S. Soils: Genesis and Geomorphology. Cambridge: Cambridge University Press. 2005.
- 72. Amorim P. K., and Batalha M.A. Soil-Vegetation Relationships in Hyperseasonal Cerrado, Seasonal Cerrado, and Wet Grassland in Emas National Park (Central Brazil). Acta Oecologica. 2007; 32:319–327.
- 73.
Bolstad PV, Vose JM. The effect of terrain position and elevation on soil C in the southern the Appalachians. In: Lal R, Kimble JM, Follett RF, Stewart BA (eds) Assessment methods for soil carbon, Boca Raton. 2001.
- 74.
Thompson L.M. & Troeh F.R. Soils, and soil fertility. 4thed. TATA McGraw-Hill Publishing company Ltd. New Delhi. 1979.