Influence of Elevation and Anthropogenic Disturbance on Woody Species Composition, Diversity, and Stand Structure in Harego Mountain Forest, Northeastern Ethiopia

Environmental variables like elevation afect species composition, diversity, distribution, density, and horizontal and upward growth. Ecologists are constantly working to better understand how species diversity varies along elevational gradients, particularly in mountainous ecosystems. Terefore, the purpose of this research was to examine the species’ horizontal and vertical structural diversity along the Harego Mountain Forest’s elevational gradient. Te area was categorized into lower, middle, and upper elevations. A total of 67 (20m × 20m) plots were created along gradients of elevation 2,079–2,516 meters above sea level (m a.s.l.). Information for foristic composition, diversity, stand structure, and environmental variables were measured and recorded for each plot over the three elevational gradients. Data on anthropogenic disturbances were visually evaluated for every plot in every gradient of elevation. For the diversity analysis, Hill’s diversity statistics were employed. To fnd signifcant variations between the three elevational gradients in terms of Hill’s diversity number, stand structure, and environmental variables, the one-way analysis of variance with SPSS version 26 at the 0.05 level of signifcance was carried out. Te fnding revealed that 50 woody plants that belonged to 35 families and 44 genera in the 67 sample plots with an elevation of 2,079 to 2,516 were identifed. Shrubs were dominant in each elevational gradient. Species richness, abundance, and Hill’s diversity number were all signifcantly ( p < 0 . 05) greater in the upper elevational gradient of the forest. On the other hand, all stand structures were signifcantly ( p < 0 . 05) higher in the middle elevational gradient. Te efect of anthropogenic disturbances and environmental variables were clearly observed in the lower and upper elevational gradients than in the middle elevation. As a result, there were fewer seedlings, saplings, trees, and shrubs in the gradients of lower and higher elevations. For the conservation of the forest, it is crucial to pay special attention to biotic elements at lower elevations and abiotic factors at higher elevations. Accordingly, involving the local community in forest management, reducing anthropogenic pressure in and around the Harego Mountain Forest through tree planting in farmlands and woodlots and implementing physical soil and water conservation structures are recommended


Introduction
Te three most signifcant characteristics of the foral ecosystem are structure, composition, and function.Knowing the temporal and spatial fuctuation patterns in species abundance and variety is among the most signifcant difculties facing ecologists and biogeographers [1,2].Te management of forests, particularly in sub-Saharan Africa, is frequently thought to be severely hampered by an absence of comprehension of the characteristics of forest composition and the processes regulating structural diversity [3].Information on forest conditions, particularly vertical and horizontal structural diversity along environmental gradients [3,4], is critical for sustainable forest ecosystem management [5].As a result, environmental gradients like elevation, slope gradient, aspect, soil conditions (e.g., moisture content, soil depth, and nutrient availability), and other factors hampers species distribution, density, and horizontal and upward growth [6][7][8][9].
Ecologists always strive to comprehend how species diversity might vary [10] of these ecological gradients, particularly in mountainous ecosystems [11].Altitudinal gradients, however, are still served as a starting point for modern studies on diferent issues in ecology and evolution [12].It is discovered to be one of the primary environmental factors structuring both the types of plant communities in the forest region [13,14] and the makeup of plant species [15].Tis is due to the fact that elevation afects the spatiotemporal distribution of other environmental factors, such as temperature, precipitation, air pressure, soil, hydrology, and others [16], which either directly or indirectly afect plant growth and development as well as the dispersal confgurations of fora in a given area [17].Elevation therefore has a signifcant potential to enhance our comprehension of species richness, species distribution patterns, and conservation [13].Terefore, an examination of the species confguration [18] and its relationship to environmental conditions is useful [19,20] in the context of the precursor to developing more efcient ways for management and maintainable exploitation of forest resources [21].
As a result of environmental gradients including elevation, slope, and aspect, certain research carried out elsewhere found a substantial association among plant composition and variety [12,17,20,22,23].Te result of topographic infuences on species variety and organization diversity [3,13,19,24] as well as the interaction between environmental variables and plant species diversity and organizational diversity [4] were also demonstrated in a study carried out in Ethiopia.Numerous research have generally examined existence of species variations along gradients of elevation [2], and the majority of them discovered a "hampered" distribution with peak species richness close to the middle gradient [25] even within a little altitudinal variation [20].Tus, it is not admiring that these intense changes over short distances have led ecologists to articulate many of the basic ecological concepts of elevational gradients [12].
However, instead of looking at the connection between elevation and species richness and diversity, the aforementioned research neglected to consider the efects of environmental gradients, specifcally elevation as well as the impact of anthropogenic disturbances on plant diversity and structure.Consequently, less research has been done on how plant species distribution responds to changes in elevation and human caused disturbances in Ethiopia in general and the northern region in particular [26].As a result, the research region lacks information on plant distributions along elevational gradients, and anthropogenic disturbances which is essential for the efective conservation of biodiversity [13].Studying how vegetation reacts to environmental conditions [4] is crucial for improving our understanding of and management of forest ecosystems [3].Tus, research into the distribution and structural arrangements of species in reaction to external conditions [20] contributes to our indulgent of ecological developments and the managing of ecologies [3,4,27].Terefore, the purpose of this research was to examine how elevation and anthropogenic disturbances afected the diversity, composition, and stand structure of woody species in the Harego Mountain Forest in northeastern Ethiopia.We hypothesized that the variability of plant species richness, diversity, and stand structure was signifcantly infuenced by elevation and human caused disturbances.

Study Site Descriptions. Northeastern Ethiopia's Harego
Mountain forest served as the study's location.Harego Mountain Forest is located between latitudes of 11 °4′ 12″ and 11 °6′ and longitudes of 39 °38′ 53″ and 39 °41′ 24″ (Figure 1).Its rough terrain, which spans 341.67 hectares and ranges at elevation from 2,079 to 2,516 meters above sea level (m a.s.l.), is what defnes it.It is made up of high mountains, valleys, and plateaus [28].Cambisols, arenosolsl lithosols, and vertisols are the main soil types in Harego Mountain Forest [29].Te Harego Mountain Forest's average annual temperature was 17.7 °C, with mean minimum and maximum temperatures of 7.7 °C and 27.6 °C, respectively.Te climatic information of Harego Mountain forest was gathered from adjacent stations [30].Te Harego Mountain Forest received 1,040 mm of rain on average year, with signifcant seasonal variations.Te research area generally experiences bimodal rainfall (Figure 2), with the major rainy period existence since July to August and low rainfall occurring from November to February [30].Te populations living near and around the study area (Harego Mountain Forest) are estimated to be 85,367, of whom 41,968 are men and 43,399 women.About 58,667 or 68.72% of the total population are inhabits in Kombolcha city, while the remaining people reside in remote kebeles [31].Diferent types of land use systems existed around the study area.Most parts of the land is covered by crops which spans 6,065 ha (48.17% of the total land) (Table 1) followed by forest land 3,763 ha (29.88%), grazing land 2,144 ha (17.03%), and others 620 ha (4.92%) [32].
Te vegetation of Harego Mountain Forest hosts various species of wild animals.Among the wild animals in the forest, Jib (Crocuta crocuta), Nebr (Panthera pardus), Tera Kebero (Canis aureus), Midakua (Silvicapra grimmia), Ses (Oreotragus oreotragus), Tera Zingero (Papio hamadryas), Tera Tota (Cercopethicus aethiops), Shikoko (Procavia capensis), Tinchel (Lepus starckii), Jart (Hystrix cristata), and Kerkero (Pacochocerus aethiopicus) are commonly found.It also provides wild edible fruits for humans and one of a beautiful area for cultural and recreational services.Hiking and celebrating diferent cultural and spiritual festivities are common in Harego Mountain Forest.However, direct poaching and habitat fragmentation brought on by human development pose a serious threat to the forest's wildlife populations.Around the forest, there are various small communities where people generally cut wood illegally.Crop cultivation is frequently practiced in the vicinity of the forest, as is livestock pressure, particularly in and around the forest's margins.Due to the forest's proximity to the cities of Kombolcha and Dessie, trees were illegally harvested for use in building, as lumber, as fuelwood, and for other uses.Tis last remaining area of forest is on the verge of extinction due to the pressure from the population's growth on other land use systems and the demand for forest products [28].a diameter of ≥2 cm and a height of ≥2 m along the gradient.Plots and transects were spaced an average of 100 and 500 meters apart in each elevational gradient, respectively, in accordance with the procedures that Worku et al. [28] and Yineger et al. [33] advised for mountainous areas.Te scientifc names of all the woody species in each plot were then noted, and references to the published volumes 1-7 of Flora of Ethiopia and Eritrea [34][35][36][37][38][39][40][41] and useful trees and shrubs for Ethiopia [42] were made.A caliper and a Suunto Clinometer have been used to quantity the diameter and height, respectively.Te DBH of each stem was measured and averaged for plants with more than one stem that were less than 1.3 meters in height [43].Te diameter of trees that were aberrant at 1.3 m and buttressed was measured directly above the buttress, at the point where the stem almost takes on a cylindrical shape [43].A diameter tape was used to measure trees that were too big to measure with the calipers.With a graded stick, the heights of small trees and plants were determined.Tree and shrub heights were assessed visually in areas where topography made measurement challenging.Additionally, fve 1 m × 1 m subplots (four at the corners and one in the middle) were used to count the seedlings and saplings in each plot.According to Sighal [44], Temesgen and Werkineh [45], and Worku et al. [28], individual stems with DBH <2 cm and height <1.5 m were counted by species as seedlings, while individuals with DBH <2 and height between 1.5 and 2 m were counted as Measurements of environmental variables were made in addition to vegetation data.Elevation, slope, aspect, soil depth, and litter concentration (depth) are among the environmental variables that are measured in each plot.We measured elevation with a portable GPS unit (GARMIN 72).A GPS compass was used to measure the aspect, and a clinometer was used to measure the slope.Te aspect was classifed using the codes found in [46] as a potential measure of solar energy overall.As a result, the following values are obtained: North

Vegetation Data
and West (W) � 2.5.Plastic metric tape was used to measure the thickness (depth) of the litter.Using a metallic stick, the soil depth was measured in the center and at each of the four corners.Te stand soil depth was then determined by averaging these measurements.In accordance with Woldemichael et al. [21], the measured depth was then categorized as follows: ≤20 cm � very shallow, >20-50 cm � shallow, >50-100 cm � moderately deep, and >100 cm � deep.

Anthropogenic Disturbances Data Collection.
In the same way, information on anthropogenic disturbances like footpaths, stem cutting, and grazing pressure was gathered from the study area.Tus, the amount of bare ground or ground vegetation cover, the quantity of faecal matter droppings left by herbivores, and the presence of animal trails and signs of trampling and browsing were used to estimate the intensity of grazing.According to Tekle et al. [47] and Woldemichael et al. [21], the grazing rate was then visually categorized as 0 � nil, 1 � slightly grazed, 2 � moderately grazed, and 3 � highly grazed.By evaluating and documenting the existence or lack of stumps, logs, and indications of frewood gathering at each sample plot, the degree of human intervention was also estimated.As a result, in accordance with Woldemichael et al. [21], the impact's magnitude was scaled as follows: 0 � negligible, 1 � low, 2 � moderate, and 3 � heavy.Additionally, the footpath's presence and absence were ascertained, and the frequency and intensity of the footpath were recorded using 0-3 scales.Ultimately, the degree of disturbance was classifed into four (0-3) categories as stated by Aynekulu [48] by taking into account the intensity of the grazing, human impact, and the footpath: 0 � nil, 1 � slightly disturbed, 2 � moderately disturbed, and 3 � highly disturbed.

Data Analysis.
Te diversity of woody species throughout the three elevation gradients (lower, middle, and upper) was evaluated using the following metrics: Fisher alpha, the Simpson diversity index [49], and the Shannon-Wiener diversity and equitability index [50,51].To ascertain the effectively existing species in each elevational gradient, Hill's diversity, or true diversity (true Shannon, true evenness, and true Simpson) [52,53], was also examined.It was determined that the Shannon diversity index was as follows: where H ′ � Shannon-Wiener diversity index, S � the number of species, Pi � the proportion of individuals in the i th species, and ln � natural logarithm � logarithm of the base e Te Shannon-Wiener index has values ranging from 0 to large numbers, which characterize communities with numerous individuals of a single species and communities with numerous species nevertheless few individuals [51,52,54].
Te value of evenness (Shannon equitability index) was calculated as follows: where E � Shannon equitability, H′ max � S, or N 0 � the number of species.Te value of evenness lies in between 0 and 1, with 1 being comprehensive evenness [55].
Te diversity of Simpson index (λ) was computed as follows: where λ = Simpson's diversity index and Pi is labeled in equation (1).Its value ranges from 0 (low diversity) to a maximum of (1 − 1/S), where S is the number of species [54].
Te Shannon-Wiener and Simpson diversity indices, however, are entropies or indexes rather than diversity in and of itself [56].Since the index is based on probabilities of occurrence and measures the uncertainty in predicting the species of a given randomly selected individual from a community, they should be converted into efective numbers of species (true diversities) [52,53].Te efective number of species produces a consistent and credible similarity measure [53].
Te value of true Shannon was calculated by exponential altering of the results gained from Shannon entropy as follows: where 1 D or N 1 � True Shannon and f � the inverse of LN, the natural logarithm of a number.Te value of true evenness (E ′ ) [57] was gained from the exponential result of evenness: International Journal of Forestry Research Te value of true Simpson was gained by inverting the result of Simpson entropy.
where 1/λ � 2 D � is true Simpson diversity index Te study of species frequency and abundance frequency ratio (A/F) [58], density and basal area [59], and Important Value Index (IVI) [55] was used to explain the structural analysis of the vegetation along the gradients.

Frequency(F) �
The total number of quadrats in which a species occurred Total number of sampled quadrats * 100,

Abundance(A) �
The total number of individuals of a species The total number of Quadrats in which the species occurred ,

Density(D) �
The total number of individuals of a species The total number of quadrats studied , where π � 3.14 and DBH � diameter measured at breast height.
Important Value Index(IVI) where RD � relative density, RBA � relative basal area, and RF � relative frequency.
To discover signifcant variations among the elevational gradients in terms of the mean of the diversity indices, stand structure, and environmental parameters, a one-way analysis of variance (ANOVA) was done using SPSS version 26 at the 0.05 level of signifcance.To identify which gradient had the signifcant diference, a post hoc test was performed using the Tukey Honestly Signifcant Diference (Tukey HSD).

Floristic Composition along Elevational
Gradients.In all, 50 woody species that belongs to thirty-fve families and fortyfour genera were found and documented in the 67 sample plots that ranged in elevation from 2,079 to 2,516 m a.s.l.(Table 1).Tere were 40, 45, and 30 species found in the lower, middle, and upper elevational gradients, respectively.Te medium elevation had the most documented number of species, while the upper level had the fewest.Similarly, the middle elevation took the uppermost number of genus (n = 40) and families (n = 31), while the upper elevation contained the lowest number of genus (n = 24) and families (n = 23) (Figure 3).Fabaceae was the greatest species-rich family in each elevation gradient, and the majority of families were exemplifed by one species each.Lamiaceae was also an equally dominant family in the middle elevation, while Anacardiaceae was the second maximum species-rich family in each elevation gradient.In the lower elevation, Fabaceae was denoted by fve species, followed by Anacardiaceae and Lamiaceae, each with three species.Celastraceae, Cupressaceae, Euphorbiaceae, Oleaceae, and Rosaceae were represented with two species each at the same elevation.Each of the remaining 17 families was represented by a single species in this elevational gradient.In the middle elevation, Fabaceae and Lamiaceae were embodied each by four species, followed by Anacardiaceae (three species), and Celastraceae, Cupressaceae, Euphorbiaceae, Oleaceae, and Rosaceae (each with two species).Each of the remaining 23 families was represented by a single species in this elevational gradient.Four species of Fabaceae were found in the upper elevation, followed by Anacardiaceae, Celastraceae, Oleaceae, and Rosaceae each with two species.Each of the remaining 18 families was represented by a single species in this elevational gradient.
In each elevation gradient, shrubs were the predominant life form, accounting for 21, 23, and 18 species of the total species recorded in the lower, middle, and upper elevations, respectively.In the lower, middle, and upper elevation  4).Five of the total species found in the research region were endemic to Ethiopia (Table 2).Every known endemic species in the study area has a shrubby growth habit.Tree out of the total endemic species, namely, Lippia adoensis, Rhus glutinosa, and Maytenus arbutifolia, were found in the entire elevational gradients.Clematis hirsuta was found in both the lower and middle elevation gradients, whereas Solanum marginatum was only found in the lower elevation.
Structural values of woody species in this research area across the three elevational gradients are presented in Table 3. Te three abundant species in the lower elevation in decreasing direction were Dodonaea angustifolia (12.82%),Myrsine africana (9.91%), and Rhus natalensis (6.77%), while Lippia adoensis (8.84%), Dodonaea angustifolia (7.71%) and Rhus natalensis (7.09%) were in the middle elevation.In contrary, Dodonaea angustifolia (9.55%), Myrsine africana (8.23%), and Maytenus arbutifolia (7.55%) were the three abundant species in the upper elevation.In the lower elevation, Euclea racemosa was occurred in all plots (100% of the plots) followed by Rhus natalensis (95.83% of the plots) and Carissa spinarum (87.50% of the plots), while Ekebergia capensis, Prunus africana, Rosa abyssincia, and Acacia seyal (each with 8.33% of the plots) were among the least frequent species.In the middle elevation, Olea europaea (90.91%),Dodonaea angustifolia, and Rhus natalensis (each with 84.85% of the plots) were among the utmost frequent species, while Rhus retinorrhoea and Olinia rochetiana (each with 3.03% of the plots) were among the minimum frequent species.In the upper elevation, Carissa spinarum, Myrsine africana, Dodonaea angustifolia, Rhus natalensis, and Jasminum abyssinicum (occurred each in 90% of the plots) were the utmost frequent species, while Acacia etbaica, Dombeya torrida, and Grewia ferruginea (occurred each in 10% of the plots) were among the least frequent species.Te furthermost densely populated species in the lower, middle, and upper elevational gradients was Dodonaea angustifolia, accounting for 22.21, 15.63, and 16.46% of the total density, respectively.Species like Euphorbia candelabrum, Olea europaea, and Juniperus procera were accounting for 37.44, 34.37, and 34.55% of the total basal area, in the lower, middle and upper elevation, respectively.Based on the species important value index (IVI), Euphorbia candelabrum (49.67% of the total IVI value), Olea europaea (36.89%), and Juniperus procera (45.95%) were also the furthermost ecologically important species in the lower, middle, and upper elevation, respectively.Clerodendrum myricoides, Prunus africana, Rhus glutinosa, and Rosa abyssincia were the rare species recorded in the three elevational gradients.Species like Acacia seyal, Echinops macrochaetus, Pinus patula, and Solanum marginatum were recorded only in the lower elevation, while Clematis hirsuta, Clutia abyssinica, Erica arborea, Festuca simensis, Otostegia integrifolia, Mimusops kummel and Berberis holstii were documented only in the middle elevation.On contrary, Rumex nervosus was only present in the upper elevation.However, Acacia abyssinica, Acacia etbaica, Asparagus africanus, Buddleja polystachya, Calpurnia aurea, Carissa spinarum, Clerodendrum myricoides, Dodonaea angustifolia, Eucalyptus globulus, Euclea racemosa, Grewia ferruginea, Heteromorpha arborescens, Jasminum abyssinicum, Juniperus procera, Maytenus arbutifolia, Maytenus senegalensis, Myrsine africana, Olea europaea, Osyris quadripartita, Pittosporum viridiforum, Prunus africana, Pterolobium stellatum, Rhus glutinosa, Rhus natalensis, and Rosa abyssincia were existed in the entire elevational gradients.International Journal of Forestry Research abundance of the forest under study were 11.37 ± 0.47 and 176.03 ± 11.19, respectively.Te overall Shannon-Wiener diversity, evenness, and Simpson diversity indices of the woody species in the studied forest were 2.04 ± 0.04, 0.86 ± 0.01, and 0.83 ± 0.01, respectively.Te mean woody species richness (observed) in the higher elevation was substantially higher than the lower and middle elevation (p < 0.05).However, the diference in woody species richness (observed) of the forest between the lower and middle elevation was not substantial.Woody species richness or Chao-1 richness (observed plus undetected) in the middle elevation was signifcantly (p < 0.05) higher than that of the lower elevation but not higher elevation.Te variation in woody species richness (observed plus undetected) between lower and upper elevations was not signifcant.Tis implies that undetected species during data collection were high in number than the species in the lower and upper elevation.Te abundance of woody species in the higher and middle elevations were signifcantly (p < 0.001) higher than that in the lower elevation.Generally, woody species richness (observed) and abundance increased as elevation increased.Woody species diversity (Shannon, Simpson, and Fisher alpha) and Hill's diversity number (Hill's N1 and Hill's N2) in the higher elevation was signifcantly (p < 0.05) higher than the diversity at the lower and middle elevation.A least Shannon and Simpson, and their true diversity, in addition to Fisher alpha, were detected in the middle elevation.However, evenness and its true evenness were signifcantly (p < 0.01) lower in the middle elevation than those in the lower and higher elevational gradient.Te highest number of woody species (Hill's N 0 � 26) was documented in the upper elevation at 2378 m, while the lowest number of woody species (Hill's N 0 � 4) was verifed in the middle elevational gradient at 2276 m (Figure 5(a)).Similarly, the upper elevational gradient had the highest Shannon diversity index (H′ � 2.74) at elevation of 2378 m, while the middle elevation had the lowest Shannon diversity index (H′ � 0.72) at elevation 2231 m (Figure 5(b)).In contrast, the highest (J � 0.95) and the lowest (J � 0.29) species evenness were recorded in the middle elevational gradient at 2164 m and 2231 m, respectively (Figure 5(c)).Te overall trend of species richness and the Shannon diversity index over an elevation gradient demonstrates that these parameters grew with elevation until it reached 2378 m, at which point they progressively decreased as elevation increased.On the contrary, as elevation rose, the species evenness distribution pattern shrank.

Stand Structure. Te stand confgurations of Harego
Mountain Forest over each elevational gradient are described in Table 5. Te mean value of tree density and basal area of the studied forest were 4,400.75± 279.83 stems ha −1 and 9.75 ± 1.29 m 2 ha −1 , respectively.Te entire three elevational gradients showed a substantial diference (p < 0.05) in the stand structures.All stand structures were higher in the middle elevational gradient.Te greatest basal area was documented (13.24 m 2 ha −1 ) in the middle elevation, while the lowest (4.88 m 2 ha −1 ) was at the upper elevation.At the middle elevation, the basal area was twofold and threefold higher than the values of the lower and the upper elevation, respectively.Similarly, the basal area in the lower elevation was twofold greater than the basal area of the upper elevational gradient.Te overall stem number (i.e., a combination of seedling, sapling, and tree and shrub density) was significantly (p < 0.01) varied by elevation.A comparable stem number was recorded at the middle and upper elevation, while the number of individual increased greatly from lower (2,976 stems ha −1 ) to middle (5,211 stems ha −1 ), but decreased to upper (5,148 stems ha −1 ) elevation.Tis means the stem number in the lower elevation was substantially (p < 0.01) 43% and 41% lower than the stem number in the middle and upper elevation, respectively.Te stem number of trees and shrubs in the lower elevation was substantially (p < 0.01) lower by 34% and 21% of the stem number of trees and shrubs in the middle and upper elevation.
Te density and basal area of stems augmented as elevation increased (Figures 6(a) and 6(b)).Te uppermost stem density (10,950 stems ha −1 ) in the Harego Mountain Forest was estimated in the middle elevational gradient at 2,316 m, while the lowermost stem density (400 stems ha −1 ) was estimated in the lower elevational gradient at 2,092 m.Similarly, the greatest basal area (57.19 m 2 ha −1 ) was estimated in the middle elevation at 2,286 m, while the lowest basal area (0.85 m 2 ha −1 ) was estimated in the lower elevation at 2,134 m.Te basal area and stem density over the gradient of elevation show an increasing pattern as elevation increased.

Environmental Variables along Elevational Gradients.
Te environmental variables of the Harego Mountain Forest for each elevational gradient are described in Table 6.Except for aspect; all environmental variables had revealed a substantial diference (p < 0.05) along the three elevational gradients.In contrast to the middle elevation, the slope's steepness was substantially (p < 0.05) higher in the lower and upper elevation.Soil depth in the lower and middle    10 International Journal of Forestry Research elevational gradients was substantially (p < 0.05) higher than that in the higher elevation.Te concentration of litter was substantially (p < 0.001) diferent in the entire three elevational gradients.Litter in the middle and higher elevation was four-times and two-times higher than the concentration at the lower elevation.Te moisture content (humidity) of Harego Mountain Forest was also signifcantly (p < 0.05) diferent among the entire three elevational gradients.Te humidity of the middle elevation was higher by 45% and 61% of the lower and higher elevation, respectively while, the humidity in the higher elevation was higher by 28% of the lower elevation.

Anthropogenic Disturbance along Elevational Gradients.
Te Kruskal-Wallis test of anthropogenic disturbance in the three elevational gradients was also described in Table 7. Te extent of grazing pressure and illegal stem cutting showed signifcant (p < 0.01) diferences in the entire three elevational gradients.Grazing pressure and illegal stem cutting were higher and moderate at the lower and higher elevational gradient, respectively, but slight at the middle elevation.Te extent of the footpath at the lower elevation was moderate, but it was slight at the middle and higher elevation.Generally, the level of anthropogenic disturbance (a combination of grazing pressure, stem cutting, and footpath) was moderate at the lower and higher elevation, while a slight disturbance in the middle elevation.

Floristic Composition along Elevational Gradients.
Obviously, foristic composition can be expressed in terms of species richness, genera, families, and life forms.Within 2,079 to 2,516 m a.s.l.elevational ranges, about 50 woody species were recorded and documented in the Harego Mountain Forest.Te fnding of this study is generally in with the Afromontane forest species richness, such as at Jibat mountain forest, with 52 species [60], and Wof-Washa forest, with 51 species [61] in the elevational range of 1,500-3,000 m.Studies also reported elsewhere at Mount Elgon and Mau forest in Kenya also reported 49 and 50 tree species, respectively, at the elevations of 1,200-2,400 m and 2,100-2,700 m [62,63], which are comparable numbers of species to the study area.
Ten again, the number of species found in the present study was higher than that of the species stated at Bale mountain (Rira) forest, 16 woody species, in the elevation range of 3,074-3,274 m [4], at Menagesha forest, 30 species [60], at Naran Valley, Pakistan, 32 woody species, in the gradient range of 2,450-4,100 m [23], at Dawsura exclosure forests, 34 woody species, in the gradient range of 1,670-2,138 m [64].However, it was less than the number of species reported in Mana Angetu forest (117 species) in the elevation range of 1500-3000 m [65], the Great Rift Valley of Tigray, northern Ethiopia (108 species) in the elevation range of 1,000-2,760 m [13], in the Eastern Escarpment of Wollo, Ethiopia (104 species) in the elevation range of 750-1,750 m [19], in the remnant moist Afromontane forest of Wondo Genet, south central Ethiopia, 72 species in the elevation range of 1,800-2,500 m [3], and Harenna forest, 61 species, in the elevation range of 1,500-3,000 m [66].Species number presented here in the current study is lower than species number documented elsewhere, such as at Southern Norway grassland with 141 species in the elevation range of 530-1,230 m [67], at Baihua mountain reserve, Bejing, with 171 species [14].Te possible reason for this variation might be due to the diference in elevation range, which in turn governs other environmental variables such as temperature, precipitation, air pressure, soil, and hydrology [16], which directly or indirectly control the growth and development of plants and the patterns of vegetation distribution in an area [17].Many authors, including Brinkmann et al. [20], Chawla et al. [68], Yirdaw et al. [4], and Zhang et al. [69], have reported on the number of species variations, even in a small elevational range diferences.
Along the elevational gradient of Harego Mountain Forest, the number of species, genera, and families were highest at the middle elevation.Species genera and families were found the highest in the middle elevational gradient, as also reported in the Eastern Escarpment of Wollo [19], in the Great Rift Valley of Tigray, northern Ethiopia [17], in the remnant moist Afromontane forest of Wondo Genet, south central Ethiopia [3], and in the Bhabha Valley in the western Himalaya [68].According to Berhanu et al. [70], the diversity of species, genera, and families in the dry Afromontane forests of Ethiopia in general, are higher at the middle elevational gradient.Te maximum diversity in species, genera, and families in the middle elevational gradient could be explained by the minimum disturbance, which states that the minimum disturbance maximizes the diversity of species, genera, and families [68].As a result of agricultural expansion and road construction in the lower elevational gradient, species distribution in the disturbed area was declined to a minimum fgure.Besides, the species, genera, and families which are found in the transition stage from the disturbance of higher and lower elevational gradients are confned to the middle elevational gradient [68].
Similar to the present study, the dominance of Fabaceae and Lamiaceae along an elevational gradient was reported in the Great Rift Valley of Tigray, northern Ethiopia [13].Furthermore, the dominance of Fabaceae in the forest along an elevation gradient was reported in the Eastern Escarpment of Wollo, Ethiopia, representing 28 species [19], in the Dawsura exclosure forests, Tigray, 8 species [64], in the Yegof mountain forest, 9 species [71], and in the Wof-Washa forest, 6 species [72].Te highest representation of species from the family Fabaceae across the three elevational gradients could be related to the fact that it is the frst largest family in the woody plants of the Afromontane forests of Ethiopia [70].Besides, Fabaceae has the largest number of woody species in the fora of Ethiopia and Eritrea, next to the family Asteraceae [40,73].Tis could also be attributed to its successful dispersal strategies and adaptation potential to the diverse agroecologies of the country [28].In line with this, Chawla et al. [68] reported families that have the potential to grow in a wide range of environmental conditions and had higher recruitment capability with resource limitations would dominate the vegetation areas.

International Journal of Forestry Research
In the current study, the biggest proportion (54%) of life forms was shrubs.A higher proportion of shrubs was also presented in other similar vegetation studies along elevational gradients, such as the Great Rift Valley of Tigray, northern Ethiopia, representing 55% [13], in the Eastern Escarpment of Wollo, Ethiopia, 43% [19], in Yegof mountain forest, 51.3% [71], and in Gra-Kahsu natural vegetation, 39.53% of the total species recorded [74].Additionally, it was observed by that shrub species predominate in Ethiopia's Afromontane forest [70].Te abundance of shrub species was also presented in the vegetation area of the Bhahala Valley in Western Himalaya, with 52% of the total woody species [68].Presence of high shrub species composition was because heavy duty logging of useful woody species used for fuelwood and timber [75].Consequently, the predominance of pioneers and shrubs changed their species composition dramatically [48,76].Furthermore, shrub species can establish themselves early in damaged areas [28,71].Chawla et al. [68] draw the conclusion that although trees are restricted to specifc height gradients, shrubs can potentially be found across a larger geographic range.In the research area, fve endemic plant species that are unique to Ethiopia, cannot be found anywhere else in the world, and require immediate conservation action.Rhus glutinosa was classifed as vulnerable in the IUCN Red List, whereas Maytenus arbutifolia was classifed as near threatened [77].Clematis hirsuta is not afected, while the other two species: Solanum marginatum and Lippia adoensis were the least.

Species Diversity and Stand Structure of Woody Species along the Elevational Gradients.
Te results of the current study showed that stand structure, species richness, and diversity are all signifcantly infuenced by elevation.Several studies have shown that elevation is a signifcant factor infuencing species diversity in mountains regions [28,67,[78][79][80][81][82][83][84][85][86][87].Te result of the current study also supports the results given by Zhang et al. [87], Cui and Zheng [86], Zhang et al. [14] and Woldu et al. [64], who reported elevation was among the most important factors that infuenced species distribution and diversity in the six subtropical mountain forests, Yunnan Province, in the subtropical Broadleaf Forests in Southern China, in the Baihua Mountain Reserve, China, and in the Dawsura exclosure, Tigray, respectively.Te highest value of woody species richness and diversity was found in the highest elevation in this study.Inversely, Austrheim [67] found that the species diversity of vascular plants on a small scale area peaked at midelevation, but diversity on a broad scale decreased continuously with elevation.Likewise, Gracia et al. [82] in the central Pyrenees, Lleida, Spain, Kebede et al. [85] in Wondo Genet forest, and Zhang et al. [69] in Mount Tai and Mount Lao, China, found that species richness and Shannon diversity index decreased continuously as elevation increased.On the other hand, Vetaas and Grytnes [78] in the Himalayan forest in Nepal, Chawla et al. [68] in the Bhabha Valley of the western Himalaya in India, and Khan et al. [23] in the Narran Valley, Pakistan, found that species diversity was higher at the middle elevation than at either the lower or higher elevation.Te diference in species richness along an elevational gradient among the vegetation areas might be due to the variation in elevation range, climatic condition, and anthropogenic disturbance.Species diversity reduction in the lower elevational gradient was mainly due to the existence of high anthropogenic disturbance [68].However, the higher efective number of species for Shannon and Simpson indices at higher and lower elevation might be due to the even distribution of species; while in the middle elevation, high species richness was encountered because of the dominance of Olea europaea and the occurrence of many rare species provides a less efective number of species.In agreement with the current study, Aynekulu [6] found a lower Shannon index within the gradient where the existence of a better dominance of Juniperus procera, which reduces the even distribution of species.In the lower and upper elevational gradients, on the other hand, shrub species that could have dense stems were abundant.Tis may have to do with past  12 International Journal of Forestry Research disturbances in the upper and lower elevational gradients as well as the secondary succession of plants, which includes many shrub species that are among the frst to emerge and survive in the region.Consequently, there are fewer species in the lower and upper elevational gradients but higher species diversity due to past disruption and current protection of the area; in contrast, this was not the case in the middle elevational gradient.In line to this, Chawla et al. [68] reported in an area where human disturbances like road building, habitation, and agricultural operations are prevalent, species richness decreases while species diversity increases.Furthermore, Austrheim [67] and Wondie et al. [79] strengthened land use processes such as farming, grazing, and fuelwood cutting and have leveled out the efects of other variables along the gradient.In line to this, Chawla et al. [68] reported that the lower species richness at the highest elevational gradient might be due to the loss of habitat diversity.
Te presence of elevational variation had also shown a signifcant diference in stand structure, in which high density of seedlings, saplings, tree, and shrubs were found and high basal area, DBH, and height in the middle elevation, followed by the upper elevation, while low in the lower elevation.Tis portrays that the forest had lower DBH of trees at lower and higher elevations.Tis might be due to the high human-caused impacts like expansion of agriculture, grazing pressure, and selective stem cutting in the lower and upper elevational gradients which infuence the survival of seedlings and saplings and reduction of the frequency of larger diameter trees.A similar result was also obtained in the northeastern escarpment of Ethiopia [6] and in the subalpine zone of west Himalaya [88].Te previous study [89] showed that the stand structure and diversity indices could vary as a result of the variation in species composition and the magnitude of disturbance involved.Terefore, in the study site, the high tree parameters in the middle elevation might be due to the dominance of Olea europaea and Eucalyptus globulus which can provide larger stem diameter and height.By contrast, low in lower elevational gradient could be due to the dominance of shrub species despite the presence of Euphorbia candelabrum tree species.In terms of stem density, the lowest in lower elevation could be mainly due to the disturbance efect.As a result, the stem number of seedlings and sapling in the lower elevation was threefold and twofold lower than that in the middle and upper elevation, respectively.In consistent with our fnding, the lower stem number at the lower elevation was also reported by a previous study [4] from the Rira forest.Livestock pressure was also stated to have an adverse impact on the natural regeneration of native woody species and ought to be managed to reverse the current trend [90].In the study area, grazing was not continuous meanwhile, high in the lower and upper elevation, but it had a slight efect in the middle elevation, which led to lower density in all developmental stages at the lower and upper elevational gradients.Besides, stem cutting for frewood and construction purposes was also higher (common) at the lower elevation of Harego Mountain Forest and reduced in the middle elevation.Te lower stem number due to the higher efect of stem cutting in the lower elevation was also reported from Bale mountain (Rira) forest in [4].Furthermore, Hegazy [83] strengthened that the intricate interface of diverse ecological elements in relative to elevation results in a variety of plant communities, vegetation belts, and habitat types.Environmental gradients like slope, aspect, moisture content, soil depth, and nutrient availability are hindering factors of species distribution, density, and population structure [7,64,91].Deep soil depth and high moisture content, slight grazing efect and illegal stem cutting, and the low frequency of footpaths in the middle elevation of the studied site resulted in better stand structure and species richness compared to the other two elevational gradients.As the range of elevation increases, the climatic condition will be also changed and the species coping mechanism will vary with a reduced climate [4], which in turn leads to a lower density of seedlings and saplings.Certain species' distribution might be restricted at the higher elevation due to comparatively damp climatic conditions [13,92].Terefore, giving high concern for biotic factors in the lower elevation while for abiotic factors in the upper elevation should be important for the conservation of the forest.

Conclusion
Harego Mountain Forest's woody species composition is noticeably similar to most of the forest of Ethiopian dry Afromontane.Te fndings revealed that elevation has a signifcant efect on the species composition, diversity, and stand confguration of the forest.Te efect of anthropogenic disturbance is clearly observed in the lower and upper elevational gradients, which results in a lower stem number of seedlings, saplings, and trees and shrubs and a low basal area than in the middle elevation.Te steepness of the slope, low soil moisture content, and shallowness of soil depth in the lower and upper elevations also contributed to the reduction of diversity and stand structure in these elevational gradients.Engaging the local community in forest management through participatory forest management system would raise the awareness of forest ownership, reduce illegal activities in the forest, and improve the regeneration process of indigenous species.Besides, reducing human pressure on forest areas through tree planting in farmlands and woodlots, as well as implementing physical soil and water conservation structures are recommended.Further investigation on seed banks of soil and the result of other physical and human factors not yet investigated on the forest are required.

Data Availability
Te datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Disclosure
Te funders had no involvement in the study's design, data collection, analysis, or interpretation, article preparation, or decision to publish the results.

Figure 3 :Figure 4 :
Figure 3: Floristic composition of collected woody species in the three elevational gradients of Harego Mountain Forest.

Figure 6 :
Figure 6: Stem density (a) and basal area patterns (b) of woody species over elevational gradients in the Harego Mountain Forest.
Collection.From September to November 2018, reconnaissance surveys were conducted in and around the forest area to obtain an overview of the site into three elevational gradients: lower [2,079-2,225], middle (2,225-2,347), and upper (2,347-2,516] elevations with equal distances.Data from the real feld was gathered between December 2018 and April 2019.Following Worku et al. [28], diameter at breast height (DBH) and height in 67 plots (20 m × 20 m) were measured for all woody plants with

Table 1 :
List of the collected plant species in Harego Mountain Forest, including their vernacular name, family, life form, and origin.

Table 2 :
Endemic woody species recorded in the three elevational gradients of Harego Mountain Forest.LC is the least concern, NA data is not available, NT is near threatened, S is a shrub, VU is vulnerable, and X is present.
6International Journal of Forestry Research gradients, trees made up roughly 16, 18, and 9 species of all the species, respectively.Conversely, lianas represent 3, 4, and 3 species of the total species in the lower, middle, and upper elevation gradients, respectively (Figure

Table 3 :
Te RF, RBA, RD, and IVI values of each species encountered in the three elevational gradients of Harego Mountain Forest.
RF � relative frequency, RD � relative density, RBA � relative basal area, and IVI � important value index.International Journal of Forestry Research

Table 4 :
Te mean (mean ± SE) richness, abundance, and diversity indices of woody species in the three elevational gradients of Harego Mountain Forest.Values represented with diferent letters revealed signifcant diference at 0.05 level of signifcance tests.Te plus-minus (±) sign used to separate standard error from the value of mean.Hill's N0 stands for species richness (observed), Hill's E1 stands for true evenness, Hill's N1 stands for true Shannon, Hill's N2 stands for true Simpson, and Chao-1 stands for species richness (observed plus undetected).

Table 5 :
Mean (±SE) DBH, height, basal area, density of seedlings, saplings, trees and shrubs, and overall woody species for each elevational gradient in Harego Mountain Forest.
Te diference between elevational gradients was investigated by analysis of variance, followed by Tukey's HSD test.Values represented with diferent letters revealed signifcant diference at 0.05 level of signifcance tests.Te plus-minus (±) sign used to separate standard error from the value of mean.

Table 6 :
Environmental variables (mean ± SE) of the three elevational gradients located in Harego Mountain Forest.

Table 7 :
Anthropogenic disturbance along the three elevational gradients in the Harego Mountain Forest.
LGP � level of grazing pressure, LISC � level of illegal stem cut, FP � level of the footpath, and LGD � level of general anthropogenic disturbance.