Ecological and floristic study of Hirmi woodland vegetation in Tigray Region, Northern Ethiopia

The dryland area in Ethiopia covers a substantial region endowed with diverse plant resources. However, the landmass has received less attention even if it has high ecological, environmental, and economic uses. The present study was conducted in Hirmi woodland vegetation, which is one of the dryland areas in Ethiopia, with the objective of investigating the floristic composition, plant community types, vegetation structure, community-environment relations and its regeneration status. Vegetation and environmental data were collected from 80 sampling plots with a size of 25 m × 25 m designated as the main plots. Diameter at breast height (DBH), height, basal area, density, vertical structure, importance value index (IVI), and frequency were computed. Species diversity and evenness were analyzed using Shannon diversity and evenness indices. The plant community types and vegetation-environment relationships were analyzed using classification and ordination tools, respectively. A total of 171 vascular plant species belonging to 135 genera and 56 families were recorded. About 5.3% of the species were endemic and near-endemic to Ethiopia. The highest number of species was recorded in families Fabaceae (16.4%) and Poaceae (11.7%) followed by Asteraceae (7.0%), Combretaceae, Lamiaceae, and Moraceae (3.5% each). Five plant communities were identified. According to the results from ordination analysis, the floristic composition of these plant communities was significantly affected by altitude, slope, sand, silt, soil organic matter, total nitrogen, and disturbance. The vegetation structure reveals that a large number of individual species was categorized in the lower DBH, frequency, and height classes. The highest Shannon diversity index and evenness values of the study area were 4.21 and 0.95, respectively. Anogeissus leiocarpa, Combretum hartmannianum, Ziziphus mucronata, Terminalia macroptera, and Acacia polyacantha were the species with high IVI. Some endemic plants were in the IUCN red list categories of Ethiopia and Eritrea. The overall regeneration status of the study area was poor because of anthropogenic disturbances and grazing pressures. Although the study area is endowed with high plant species diversity including endemism, it is under poor regeneration status due to various disturbances. To overcome this challenge, integrated management measures including monitoring and application of restoration techniques by taking into consideration the significant environmental factors associated with species diversity as well as observed regeneration status and IUCN threat level of the species are highly recommended.


Background
Terrestrial biodiversity represents the richest diversity in the tropics and near to the equator (Gaston 2000), which seems to be the result of the warm climate, high primary productivity and diverse topography in the region (Field et al. 2009). About 40% of the global tropical forest area (Mayaux et al. 2005) and 14% of the total African land surface (Eshete et al. 2011) is covered by the dry woodland vegetation. According to Eshete et al. (2011), 25% of the natural vegetation in Africa is found in the dry woodland vegetation. Though the dryland area has high ecological, environmental, and economic importance, currently, the vegetation resources in this landmass are facing serious problems due to human interferences and very poor stakeholders' attention compared to the moist vegetation (Muys et al. 2006). Deforestation and habitat degradation in the dryland area resulted several socioeconomic and environmental challenges that bring a strong impact on the capacity of the forests to provide rich ecosystem services (Nyssen et al. 2007). To mitigate these challenges, conservation of biodiversity in the dryland ecosystems has been receiving greater international attention from time to time (CBD 2006;Lemenih and Bongers 2011). Effective conservation activities in dryland forests definitely demand an exhaustive investigation on environmental determinants such as soil properties, elevation, habitat suitability, habitat diversity, and species interaction (Arponen 2012; Yismaw et al. 2014).
Ethiopia shares more than 50% of the Afromontane region's land areas of Africa which are above 1500 m (Bekele 1993) as well as an extensive dryland region (Teketay 2001). According to a study conducted by Teketay (2001), about 65% of the Ethiopia landmass is covered by dryland ecosystems, where 12% of human and 20% of livestock populations reside. It is also the center of biodiversity, endemism, and grazing (Mengistu et al. 2005).
The wide ranges of agro-ecologies in northern Ethiopia are covered with dry forests or woodland vegetation types (Georgis 2010). However, due to the rapid human growth and agricultural expansions, there were habitat fragmentation and resource exploitations in these areas. Following the long history of land degradation, many land rehabilitation and conservation programs have been carried out in northern Ethiopia in general, Tigray region in particular. However, inadequate investigation of ecological baseline information makes the effectiveness of biodiversity conservation policy and strategy to fall in a challenge for its implementation as it intends (EFAP 1994;Teketay 2001).
Hirmi woodland is one of the dryland vegetation types with potential remnant vegetation stretched in Northern Ethiopia, Tigray region, north to the Tekeze Rift valley (Hagos et al. 1999). An assessment conducted by the Ethiopian REDD+ (MEF 2015) secretariat reveals Hirmi woodland is identified as one of the vegetation potential areas that require ecosystem restoration and conservation by the participation of the local community, the government, and other stakeholders. Investigation of plant community types and the driving environmental and anthropogenic factors shaping the existing vegetation of the dryland area could serve as the basis for taking effective conservation measures of the dwindling natural resource (Pimm and Raven 2000;Brooks et al. 2002). There are no earlier investigations done on plant diversity, composition, structure, threats, and regeneration status of the Hirmi vegetation so far. Thus, this study is conducted to fill the existing gap and serve as a tool for running effective conservation and sustainable utilization of the study vegetation.
The objectives of the research were to: i. examine the floristic composition of Hirmi woodland vegetation; ii. identify the plant community types of the study vegetation; iii. assess the vegetation structure of the study area; iv. examine the relationship between selected environmental determinants and the community types of the study vegetation; v. assess the regeneration status of the study vegetation; vi. identify existing threats to the Hirmi woodland vegetation.

Study area
Hirmi woodland vegetation is stretched in three districts (Tahtay Koraro, Medebay Zana and Asgede Tsimbla) of Northwest Tigray National Regional Administration, Ethiopia. The area of Hirmi woodland vegetation is found at 1100 km distance from Addis Ababa (Ethiopian capital) between 13°49'-14°04' latitude and of 38°14'-38°25' longitude, with elevation range between 900 and 2002 m.a.s.l. (Fig. 1) years (1998-2017) climatic data was taken from the nearest Endabaguna metrological station so as to describe the climate of the study area. Analysis of the meteorological data showed that the study area has a maximum and minimum monthly mean temperature of  31.9°C and 13.1°C, respectively. The rainfall of the study site was unimodal with a mean annual of 888 mm. The climate diagram (Fig. 2) showed that the study area receives unimodal rainfall of about 888 mm/year with high precipitation in August and low precipitation in February.

Method of data collection Vegetation data
To evaluate and estimate whether the number of species (floristic composition) in the study area has been sufficiently sampled or not, species accumulation curves were plotted using free statistical software R version 3.6.1. Accordingly, vegetation data were collected from a total of 80 plots along 10 transects laid systematically following the sampling approach as described by Braun-Blanquet (1964) and Mueller-Dombois and Ellenberg (1974). Each of the plots was laid at every 50 m altitudinal drop along the transect lines. The distance between each transect ranged from 500 to 700 m apart. For the collection of woody (shrubs and trees) species, plots of 25 m × 25 m (625 m 2 ) were laid. Herbaceous species, seedlings, and saplings of all woody plants were recorded in a 2 m × 2 m (4 m 2 ) subplot inside each main plot, four from each corner and one at the center. Plants with a height of ≥2 m and DBH ≥2 cm were considered shrubs or trees whereas woody species with height 1-2 m and height ≤1 m were considered saplings and seedlings, respectively (Senbeta and Teketay 2001;Tuxill and Nabhan 2001;Tadesse 2015). Vegetation coverabundance was estimated based on the Van der Maarel scale (1979). Plant specimens were collected, pressed, dried properly during the data collection period at the field site, and brought to the National Herbarium, Addis Ababa University for identification and verification. An attempt was also made to record the vernacular names of species encountered in the study during collection whenever possible.

Environmental data collection
Environmental variables such as altitude, aspect, slope, and geographical coordinates were taken and recorded for each plot. Soil samples were taken from five points (four from each corner and one from the center) measuring 15 × 15 cm, at a depth of 0-30 cm (Tenkir 2006;Tadesse et al. 2008) of the main plot laid for vegetation data collection. These samples were mixed to form a composite sample, and about a kilogram of the sample was taken for analysis. The soil samples were air-dried by spreading on plastic trays crushed and sieved with a mesh size of 2 mm to remove root particles and other organic debris before chemical analyses. Finally, measurements and analysis for soil texture (sand, silt, and clay), pH, electrical conductivity (EC), cation exchange capacity (CEC), soil organic matter (SOM), total nitrogen (N), and available phosphorus (P) were done following the standard procedures outlined in Landon (1991) and Sertsu and Bekele (2000).
Disturbances were recorded as present or absent in the sampled plots. The magnitude of disturbance was rated from 0 to 4 based on visible signs of vegetation disturbance parameters such as tree cutting, debarking, grazing, fire and charcoal production signs following Hadera (2000) and Senbeta et al. (2007). Values were coded as follows: 0: no disturbance; 1: if anyone of the above-mentioned factors was present (slightly disturbed); 2: if any two of the above-mentioned factors were present (moderately disturbed); 3: if any three of the above-mentioned factors were present (highly disturbed); 4: if all of the above-mentioned and other associated natural factors such as tree stumps and landslide/degradation were present (extremely disturbed).

Data analysis Multivariate and diversity analysis
The study vegetation was classified into community types by the Hierarchical Agglomerative Clustering technique using Ward's method employed in R packages. The plant community types from the cluster analysis were also analyzed in a synoptic table based on cover-abundance values (van der Maarel 1979). The community was named after the dominant tree or shrub species in that specific community with higher synoptic values. In each plant community type, the synoptic values of each species were further tested to find significant indicator value of the species performed in R using package labdsv. CCA was used to correlate the identified plant community types with selected environmental and disturbance factors. The reason for using CCA ordination was due to the first axis in the detrended correspondence analysis (DCA) was >4 (7.13) and the study was used heterogeneous environmental datasets (Lepš and Šmilauer 2001). In this analysis, environmental factors including altitude, slope, soil physical, and chemical properties (pH, electrical conductivity, soil texture, available phosphorus, soil organic matter, total nitrogen, and cation exchange capacity) as well as disturbance factors (tree cutting, debarking, grazing, fire and charcoal production sign) were assessed and their association with plant community types were investigated.
Shannon index was used to compute species diversity, richness, and evenness of Hirmi vegetation using R software version 3.6.1 (R Development Core Team 2019) as indicated below: Shannon diversity index (H') is calculated from the equation: Shannon evenness index (J) was calculated using the equation: where H' = Shannon diversity index; Σ = summation symbol; pi = the proportion of individuals or the abundance of ith species expressed as a proportional of total cover in the sample and ln = log base (natural logarithms); J = evenness; H max = lnS; S = total number of species. H´= 0 if there is only one species in the sample and H' is maximum (ranges from 1.5 up to 3.5 and rarely up to 4.5) when all species are represented by the same number of individuals (Kent 2012;Woldu 2016). The values of J ranges normally between 0 and 1, where 1 represents a situation in which all species are equally abundant (Magurran 1988).

Vegetation structure
Vegetation structure of the Hirmi woodland vegetation was determined by computing the following parameters.
Diameter at breast height (DBH): DBH was measured (at 1.3 m above ground) for trees and shrubs having circumferences ≥6 cm and height ≥2 m. For tree/shrub species which branch at breast height or below, DBH of each branch was measured separately and the average was taken and then treated as a single individual for the basal area calculation. There is a direct relationship between DBH and basal area (Kent 2012).
Vertical structure: The vertical structure of the woody species was analyzed using the International Union for Forestry Research Organization (IUFRO) classification scheme (Lamprecht 1989). The IUFRO classification scheme categorizes storey into upper, when the tree/ shrub height is >2/3 of the top height; middle, for the tree/shrub height between 1/3 and 2/3 of the top height and lower; for the tree/shrub height <1/3 of the top height.
Basal area (BA) and relative dominance: Basal area of a species is used to calculate species dominance and it is computed from diameter at breast height (d) as follows: where; π = 3.14 Frequency: Defined as the probability or chance of finding a species in a given plot. It was computed using the following formula: Importance value index (IVI): is useful to compare the ecological significance of species and helps to measure conservation priority among groups of species (Kent 2012). The IVI is the cumulative value of the three parameters (i.e., relative frequency, relative density, and relative dominance), calculated using the following formula:

Regeneration status
Regeneration status was estimated based on the composition and density of seedlings and saplings of all woody species recorded in each plot (Khan et al. 1987). The density of seedlings and saplings were calculated using: Density of sapling=seedling ¼ the number of sapling=seedling area of sample in a hectare The information on the density of seedlings and saplings is important to take conservation measurements based on their status.

Floristic composition
The species accumulation curve graph for the Hirmi vegetation was leveled off as the sampling effort proceeded. This was showing that the rate of discovery to new species would be getting less and less as the further sampling continues (>80 plots). This would be proved that representative samples have been taken to study the floristic composition of the study area. Thus, a total of 171 plant species belonging to 138 genera and 58 families was recorded and identified in the 80 plots of Hirmi woodland vegetation. Regarding the habit of species, herbs, shrubs, trees, and liana were represented by 79 (46.2%), 43 (25.2%), 39 (22.8%), and 10 (5.8%) species respectively. The highest number of species was recorded for the families Fabaceae (28 species, 16.4%) and Poaceae (20 species, 11.7%) followed by Asteraceae (12 species, 7.0%), Lamiaceae, Combretaceae, and Moraceae (7 species, 3.5% each). The remaining 52 families were represented by ≤4 species. Out of the total plant species identified from the study area, 9 (5.3%) species were found endemic (7 species) and near-endemic (2 species) to Ethiopia according to the flora of Ethiopia and Eritrea and Friis et al. (2010). Of these, Acacia venosa, Aloe elegans, Bidens macroptera, Lippia adoensis, Pennisetum glaucifolium, Phragmanthera macrosolen and Urtica simensis were endemic and Combretum hartmannianum and Combretum rochetianum were near-endemic species. Based on the IUCN criteria of threat level, one species (Acacia venosa) was recorded as endangered, while two species (Combretum hartmannianum and Combretum rochetianum) were designated as vulnerable and the remaining endemic species were identified as least concern and not evaluated (Table 1).

Plant community types
Five plant community types were identified from the agglomerative hierarchical cluster analysis (Ward's method) based on the abundance data of the species on the study plots (Fig. 3). Two species with relatively high synoptic and significant indicator values in the group were used to name the corresponding communities ( Table 2). The descriptions of the community types identified from the study vegetation are given below.

Ziziphus mucronata-Acacia polyacantha community type
The altitudinal ranges of this community were from 1170-1565 m.a.s.l. At the edge and lower altitude of the community, high disturbance signs such as cutting, debarking, grazing were documented. In contrast to other community types, this community was situated at a moderate slope which ranges from 2-35%. The shrub and tree layers in this community type were undifferentiated. Ziziphus mucronata and Acacia polyacantha were characteristics species of the community. The common tree species in this community were Acacia polyacantha, Ficus ingens, Terminalia macroptera, Entada abyssinica and Lonchocarpus laxiflorus whereas Ziziphus mucronata, Ziziphus spina-christi, Flueggea virosa, Dichrostachys cinerea, and Senna obtusifolia were the among dominant shrub species. The common species in the field (herb) layer were Bidens pilosa, Achyranthes aspera, Chamaecrista mimosoides, Trifolium campestre, and Setaria megaphylla.

Combretum hartmannianum-Terminalia macroptera community type
This community was found within the altitudinal range of 1098-1300 m.a.s.l. at the south edge of the vegetation. It was the most disturbed community. The community was dominated by undifferentiated Combretum and Terminalia species. Combretum hartmannianum and Terminalia macroptera were the characteristic tree species in the community. Species including Lonchocarpus laxiflorus, Ficus ingens, Anogeissus leiocarpa, and Combretum rochetianum were also commonly found in the tree layer. Lannea fruticosa was the characteristic species in the shrub layer. Grewia ferruginea and Flueggea virosa were other prominent shrub species in the community. The Oxytenanthera abyssinica community in the study area is found in this community type. It was recorded almost in all plots (in 10/11 plots) of the community (Fig. 4). Setaria megaphylla and Oplismenus burmannii were the characteristic species in the field layer. Species such as Senna obtusifolia, Setaria megaphylla, Oplismenus burmannii, Acanthospermum hispidum, and Bidens pilosa were also commonly found in the field layer.

Anogeissus leiocarpa-Ozoroa insignis community type
This community was found within an altitudinal range of 1400-1840 m.a.s.l. with an average slope of 41%. The majority (87.5%) of the plots were taken from (near to) the center of the study vegetation. It was a relatively species-rich community harboring 100 species. Anogeissus leiocarpa, Ozoroa insignis, and Terminalia laxiflora were the characteristic species in the tree layer. The community was dominated by tree species including Anogeissus leiocarpa, Ozoroa insignis, Acacia abyssinica, Terminalia laxiflora and Combretum hartmannianum. Bidens macroptera, Guizotia scabra, Bidens pilosa, Oplismenus burmannii, and Eleusine africana were among the species in the field layer that were codominant with the tree layers. Although there were no prominent species in the shrub layer, however, species such as Vangueria madagascariensis, Rhus retinorrhoea, and Acokanthera schimperi were among the common shrub species in the community.

Euclea racemosa-Acacia abyssinica community type
This community type was found in the altitudinal range between 1405 and 2002 m.a.s.l. and on gentle to strong slopes (8-75%). In other words, the highest altitudes and strongest slopes of the study area were situated in this community. Balanites aegyptiaca, Combretum molle, Ximenia americana, Ficus salicifolia, and Justicia schimperiana were species at the lower altitudes of the community type. At the upper altitudes (>1800 m.a.s.l) plant species that belong to the DAF vegetation type (Friis et al. 2010) such as, Olea europaea subsp. cuspidata, Acacia lahai, Acacia abyssinica, Croton macrostachyus, Erythrina abyssinica, and Calpurnia aurea were recorded. The community was dominated by shrub species that include Euclea racemosa, Galiniera saxifraga, Acokanthera schimperi, Acacia seyal, Rhus retinorrhoea, and Diospyros abyssinica. Out of those Euclea racemosa, Galiniera saxifrage and Acokanthera schimperi were the characteristic species in the shrub layer. Acacia abyssinica and Combretum molle were the characteristic and dominant species in the tree layer. Although no prominent herb species were recorded for this community, few individual herbaceous species including Bidens macroptera, Bidens pilosa, and Desmodium repandum were recorded.

Dodonaea angustifolia-Flueggea virosa community type
This community occurs at altitudes from 1285 to 1790 m.a.s.l. Various disturbances (cutting, grazing, and fire) were recorded in most of the sampled plots. In the elevated and moderate slopes of the community, two  Arthraxon micans, Chamaecrista mimosoides and Dactyloctenium aegyptium were the characteristics of herb species in the field layer.

Ordination
The CCA ordination result (Fig. 4) revealed that environmental variables including altitudes, slope, disturbance, soil organic matter, total nitrogen, sand, and silt soil type were significantly associated with community types identified in the Hirmi woodland vegetation. Total nitrogen, disturbances, and sand were positively correlated with the community type one (Ziziphus mucronata-Acacia polyacantha) and two (Combretum hartmannianum-Terminalia macroptera). Community type three (Anogeissus leiocarpa-Ozoroa insignis) was found in scattered form and affected by different factors including soil organic matter, slope, silt, and total Species richness, evenness, and diversity of the five community types Shannon diversity index and evenness of the identified five community types are shown in Table 3. In terms of the actual number of species per community, community number three (Anogeissus leiocarpa-Ozoroa insignis) was found the most species-rich followed by community number five (Dodonaea angustifolia-Flueggea virosa) and community number four (Euclea racemosa-Acacia abyssinica). Community number two (Combretum hartmannianum-Terminalia macroptera) was the most species-poor community. Similarly, community number five had the highest species diversity (H' = 4.21) while community number two had the lowest diversity (H' = 3.25). In terms of evenness, community V had the highest evenness index (0.95) and community two and three had the highest similar species or the lowest evenness of species (0.90 each) when compared with the other community types.

Vegetation structure DBH and height distribution
Woody species (height ≥2 m and DBH ≥2 cm) were classified into seven DBH (cm) classes: 2-10; 10.1-20; 20.1-30; 30.1-40; 40.1-60; 60.1-80, and >80 cm. As the DBH class size increased, the number of individuals was found gradually decreasing towards the successive higher DBH classes (Fig. 5). Sterculia setigera, Ficus vasta, Ficus palmata, Adansonia digitata, and Terminalia laxiflora were found to be the dominant large-sized DBH trees in Hirmi woodland vegetation. Anogeissus leiocarpa, Ziziphus mucronata, and Combretum hartmannianum were among the species which contributed a high cumulative proportion of DBH (>14% of the total DBH), whereas Buddleja cordata, Clerodendrum myricoides, and Otostegia fruticosa contributed the least to the total DBH (< 0.06% of the total DBH). Individuals of woody species with height ≥2 m and DBH ≥2 cm were classified into seven height classes and their densities were calculated. The first two height classes were accounted for 73% of the total height class while the remaining five classes were accounted for only 27% of the total height class. The cumulative height class of Hirmi woodland vegetation has an inverted J-shape which means the decline in the density of woody plants with increasing height classes (Fig. 6).

Vertical structure
The identified most top height for trees in the Hirmi woodland vegetation was 25 m. According to the IUFRO classification scheme, 65.2% of the species were in the lower storey (<8.3 m) followed by the middle (28.6%) and upper (6.2%) storeys. The lower storey of the Hirmi vegetation was found to contain high numbers of individual species but the upper storey was represented by very few individuals. Woody species present in the lower storey were Acokanthera schimperi, Calpurnia aurea, Dichrostachys cinerea, Dodonaea angustifolia, Flueggea virosa, Maytenus arbutifolia, and Ximenia americana. The wordy plants in the middle storey (8.3-16.6 m) were Albizia malacophylla, Ziziphus mucronata, Entada abyssinica, Diospyros mespiliformis, and Euclea racemosa. Species such as Acacia albida, Combretum hartmannianum, Stereospermum kunthianum, Acacia abyssinica, Anogeissus leiocarpa, Combretum adenogonium, Ficus vasta, Sterculia setigera, and Terminalia macroptera were among the species categorized in the upper storey (>16.6 m).

Density of woody species
A total of 2642 woody species were recorded from 80 plots of the Hirmi woodland vegetation. As indicated in Table 4, only 10 species accounted for 40.4% of the total density. On the other hand, the remaining woody species altogether accounted for 59.6% of the total density in the Hirmi woodland.
The computed frequency values of species were classified into five frequency classes as: 0-10%; 10.1-20%;  20.1-30%; 30.1-40%, and >40% (Fig. 7). The result showed that 59% of the total woody species of Hirmi woodland vegetation were distributed in the first frequency class. Besides, the grand frequency distribution of woody species in Hirmi woodland showed that higher proportions of species were distributed in the first frequency class followed by a rapid decline in the distribution of species across successive higher frequency classes. Anogeissus leiocarpa was the most frequent species (40.5%), followed by Terminalia macroptera (28.5%), Ziziphus mucronata, and Combretum hartmannianum (26.25% each), whereas Acacia nilotica, Olea europaea subsp. cuspidata, Otostegia fruticosa, Phoenix reclinata, Phytolacca dodecandra, and Syzygium guineense were species with low frequencies (2.5% each).

Importance value index (IVI)
The importance value index of woody species in the study area ranged from 0.35 to 20.55 (Table 5). The first top ten woody species that contribute 37.78% of the IVI in decreasing order were Anogeissus leiocarpa, Combretum hartmannianum, Ziziphus mucronata, Terminalia macroptera, Acacia polyacantha, Flueggea virosa, Dodonaea angustifolia, Diospyros abyssinica, Ozoroa insignis, and Diospyros mespiliformis. The 62.22% of the IVI was contributed by the remaining woody species. Some of the species with the least IVI values were Acacia nilotica, Buddleja cordata, Cordia Africana, Olea europaea subsp. cuspidata, Otostegia fruticosa, and Phytolacca dodecandra.

Regeneration status
The regeneration status of Hirmi woodland was characterized by recording and counting of seedlings, saplings, and adults of trees/shrubs. The total density of seedlings, saplings and mature tree, and shrub species in Hirmi woodland vegetation were 2750 ha −1 , 3025 ha −1 , and 528.4 ha −1 respectively (Fig. 8).
Assessments of seedling and sapling densities were used to prioritize conservation measurement (Table 6). Accordingly, urgent/immediate conservation action is required for those species that fall short in their seedlings and saplings.

Floristic composition and diversity
The result of the species accumulation curve plotted for the species recorded in Hirmi vegetation showed that adequate numbers of plots were laid which was confirmed by the presence of rich species composition (171 species). The majority (>85%) of the study vegetation was composed of Combretum-Terminalia and Acacia-Commiphora woodlands. There is also dry evergreen Afromontane forest (DAF) that component in the massifs and elevated parts of the study area. Besides, herbaceous species that belong to the family Poaceae, Asteraceae, and Fabaceae were predominant in the vegetation of the study area. The existence of high species diversity in the study area could be due to the fact of heterogeneity (existing different vegetation patches) within the vegetation community (Murphy and Lugo 1986), the altitudinal difference (Bekele 1993), anthropogenic impacts (disturbance) (Lyaruu et al. 2000) and soil properties (Kooch et al. 2009). The species richness of Hirmi woodland vegetation is higher than other woodland vegetation in Metema for which 87 species were recorded (Adamu et al. 2012) and more or less equivalent to that of the Dello Menna woodland vegetation (Didita et al. 2010) but the recorded species richness was less than that of the Ilu Gelan District Combretum-Terminalia mixed Acacia vegetation reported for 214 species (Tadesse et al. 2018). The recorded high richness of herbaceous species might be attributed to the open canopy cover of the vegetation that enabled for their free growth. This was in agreement with Murphy and Lugo  1986), who reported that an abundance of herbs has an inverse relation with canopy cover. Anthropogenic impacts such as exploiting plant products and selective cutting of trees Yineger et al. 2008) and dominance of shrubby species in dryland area (Lemenih and Bongers 2011) are also reported as factors responsible for the dominant growth of herbaceous species in an area. More than half (60%) of the total recorded plant family were represented by one species only, which indicates the diversity of species in Hirmi woodland vegetation. The presence of high numbers of species for the families Fabaceae, Poaceae, Asteraceae, and Lamiaceae conforms to the fact that they are species-rich families in the flora Ethiopia and Eritrea (Tadesse 2004). Out of the recorded 9 endemic species, Acacia venosa and Lippia adoensis have been included in the IUCN red data list of Ethiopia and Eritrea (Vivero et al. 2005). This could be associated with the existing significant anthropogenic and livestock disturbances in the study area. The dryland vegetation types in Ethiopia are currently under strong environmental stress (Friis et al. 2010). An ecosystem approach of biodiversity conservation and participatory forest management provides an opportunity to conserve large numbers of species that are rare and endemic (Aynekulu 2011; Lemenih and Bongers 2011).
The higher values of Shannon diversity index (4.21) and evenness (0.95) designate that the study area is equipped with high species diversity with more even distribution species within the study plots. This could be associated with the observed topographic and/or altitudinal variation (ranged from 1098-2002 m) and the species richness existing in the study area. In line with this analysis, Suratman (2012) reported that differences in terrain determinants (elevation and slope) cause differences in the soils, water, and microclimate conditions which in turn affect species diversity and adaptability. Magurran (2004) also determine the species richness in  an area that has a substantial attribution for species diversity Plant community and environmental relationship (ordination) The relationship between plant communities and environmental variables is important in understanding the plant communities in a given ecosystem (McCune and Grace 2002). Patterns of species community are associated with environmental factors including climate, topography, and soil as well as anthropogenic and livestock disturbances (Whittaker et al. 2003). Ecological studies conducted in Ethiopia by Woldu et al. (1989), Friis (1992), Bekele (1994), Aynekulu (2011), Kelbessa and Girma (2011) and Gebrehiwot et al. (2019) reported the importance of environmental variables (altitude, slope, landscape, disturbances, and soil physical and chemical properties) in shaping plant communities. CCA analysis result showed that environmental variables including altitude, slope, silt, sand, soil organic matter, total nitrogen, and disturbances are associated with species compositions of plant communities in the Hirmi woodland. However, all those environmental determinants exhibited different effects on the five identified communities. This was in line with a study conducted by Brinkmann et al. (2009) which stated the plant community has its distribution area with a specific combination of environmental variables.
Community number one (Ziziphus mucronata-Acacia polyacantha) and community number two (Combretum hartmannianum-Terminalia macroptera) were strongly influenced by sand and total nitrogen. Nitrogen content of soil has been reported as a major factor limiting species growth and diversity (Brown 2003). Disturbances were also common and significantly affected these communities. The observed disturbance might be related to the vicinity of the two communities close to human settlements (i.e., exposed to edge effect). Disturbance impacts including selective cutting and grazing at the edge of vegetation are frequently revealed (Vieira et.al. 2006).
Anogeissus leiocarpa-Ozoroa insignis (Community number three) was positively correlated with soil organic matter, silt, nitrogen, and slope. This could be related to the fact that soil having a high content of organic matter, nitrogen and moisture-holding capacity of the silt loam soil type is suitable for species growth (Brown 2003;O'Geen 2006). This community was less disturbed and had high species richness due to the fact that plots of the community were taken in sites far from the edge effect and sloping areas. Species found on strong slopes are less disturbed and would be rich in species composition (Tilahun et al. 2011;Gebrehiwot et al. 2019).
In the ordination, the axis altitude was positively correlated with community type number four (Euclea racemosa-Acacia abyssinica). The length of arrows of the environmental determinants in the ordination diagram indicates the strength of their impacts on species composition of plant communities and their distribution (ter Braak 1987). This community was distributed along a wide elevation gradient (1405-2002 m.a.s.l.) and has a high Shannon evenness index (0.93). The altitudinal difference does not only influence the species richness and diversity but also affects the other determinants that play a significant role in their distribution and composition such as the soil texture, atmospheric pressure, moisture, and temperature (Bekele, 1994;Lovett and Wasser 1993). High Shanon diversity (H' = 4.21) and evenness (J = 0.95) in community number five (Dodonaea angustifolia-Flueggea virosa) was positively associated with silt loam soil and disturbances. Silt loam soil contains the perfect combination of sand, silt, and clay particles and moisture-holding capacity that supports the growth of virtually all forms of plant life (Ruwanza 2018). The disturbed drylands of Ethiopia are covered by shrubs like Dodonaea angustifolia and secondary forests (Lemenih and Bongers 2011;Woldu et al. 2019).

Vegetation structure
Height and DBH: the overall height and DBH class distribution showed an inverted J-shape distribution. The existence of a high number of species at the lower DBH and height classes was due to the dominance of a smallsized but a large number of individuals including Dodonaea angustifolia, Flueggea virosa, Maytenus arbutifolia, Acokanthera schimperi, Bersama abyssinica, Ziziphus mucronata, and Ximenia americana. The vegetation has more species that belong to a secondary vegetation type thus lacks a large statured tree. A similar result was reported in Abergele dry woodlands (Eshete et al. 2011), Dawsura-Tembien woodlands  and Ilu Gelan Combretum-Terminalia and Acacia vegetation (Tadesse et.al. 2018). In other aspects, large-sized DBH are selectively harvested and diminished their quantity since they are preferable for various human services (Hadera 2000).
Vertical structure: based on the Lamprecht (1989) scheme of classification, the highest proportion of species were concentrated in the lower storey followed by the middle and upper storey. The finding might be associated with the dominance of shrubby and small-sized species composition of the study area. A study by White (1983) reported shrubby layer vegetation is mostly dominated by lower undifferentiated woody species. A similar result was also reported in different woodland vegetation of Ethiopia (Lemessa 2009;Tadesse et.al (2018).
Density: the overall density of woody species in Hirmi woodland was 528.4 individuals per hectare. This was greater than the density reported from other woodland vegetation studied in Metema, Ethiopia (18.8 stems/ha) (Adamu et al. 2012). However, it was less than the study in woodland vegetation of the Borana lowland (3149 stems/ha) (Dale 2004), dry land vegetation of Welo (3146 stems/ha) (Tadesse et al. 2008) and Shai Hills dryland in Ghana (959/ha) (Swaine et al. 1990). Such variation could be attributed to variations in topographic gradients, habitat preferences and the degree of anthropogenic disturbances (Whittaker et al. 2003).
Basal area: Basal area of the Hirmi woodland (14 m 2 / ha) was higher than the study conducted in the Taltalle woodland area (0.44 m 2 /ha) (Lemessa 2009). However, it was smaller than other woodland vegetation studied in Sire Beggo woodland vegetation (19.3 m 2 /ha) (Dibaba 2014) and Dirki and Jato woodland (18.95 m 2 /ha) (Tadesse 2018). More than 50% of the total basal area of Hirmi woodland was contributed by a few large-sized tree species. However, the richness of these tree species was far less than the shrub species in Hirmi woodland. An area covered with shrub layers and small DBH individuals would have a small total basal area contribution in contrast to tree species (Murphy and Lugo 1986; Lemenih and Bongers 2011).
Frequency: Based on the analysis of frequency by Lamprecht (1989) of frequency, high value in the lower frequency classes showed a heterogeneous species composition of Hirmi woodland. Anogeissus leiocarpa, Combretum hartmannianum, Terminalia macroptera, Ziziphus mucronata, and Ziziphus spina-christi were among the species with high-frequency values, whereas Phoenix reclinata, Phytolacca dodecandra, Olea europaea subsp. cuspidata, Acacia nilotica, and Cordia africana had small frequency values. The frequency of a species always depends on factors that relate to habitat preferences, adaptation, and degree of exploitation and availability of suitable environmental conditions for regeneration (Rey et al. 2000).

Important value index (IVI)
High density and frequency coupled with high basal area indicate the overall dominance or high importance of the species for the site (Lamprecht 1989). Species with high IVI in the present study were different from other woodland vegetation studied in Dello Menna woodland (Didita et al. 2010), dryland vegetation of Welo (Tadesse et al. 2008); but there was some extent of similarities with a species found in woodland vegetation of Gambela region (Awas et al. 2001), Metema and Abergele dry woodlands (Eshete et al. 2011) and Dawsura-Temben by Woldu et al. (2019). This was associated with the topographic gradients and geographical distribution of the species (Clark et al. 1998;Whittaker et al. 2003). Species with high IVI (> 9) values including Anogeissus leiocarpa, Combretum hartmannianum, Ziziphus spina-christi, Terminalia macroptera, Acacia polyacantha, Flueggea virosa, and Dodonaea angustifolia were the most abundant, dominant, and frequent species in the Hirmi woodland vegetation. On the other side, species such as Cordia Africana, Olea europaea, Otostegia fruticosa, Phytolacca dodecandra, Buddleja cordata, and Acacia nilotica that were scored least IVI (<0.5) also were found rarely and restricted to some sites of the study woodland vegetation. The variability of IVI values of the species within the woodland vegetation attributes the conservation status of some of the highly dominant species. IVI reflects the extent of dominance, occurrence and abundance of a given species in relation to other associated species as well as the ecological significance of a given species (Kent 2012).

Regeneration status of Hirmi woodland
The density of saplings was higher than the density of the seedlings and mature trees/shrubs. Accordingly, the regeneration status of Hirmi woodland is categorized as poor, following the assessment methods set by Khan et al. (1987). Dryland vegetation of Ethiopia is characterized by poor and fluctuating patterns of regeneration status due to periodic disturbances (Yebeyen 2006;Tesfaye 2008;Lemenih and Bongers 2011). Frequent signs of disturbances were recorded in the community types near the settlement and edge of the Hirmi vegetation. The ordination result also proves that disturbance was significantly associated with the most identified communities (community numbers 1, 2, and 5). Land use/ land cover dynamics study conducted by Gebrelibanos and Assen (2015) and land use/land cover mapping done by Girma (2014) indicated that human settlements around the present study area were restricted to mountainous landmasses. However, the rapid human population growth from time to time imposed settlement in lowland areas for the demand of large farmland as well as plant and plant products. Thus, the massive introduction of humans and livestock in those areas led to poor regeneration of species. As a result, many species were represented with few/no seedlings or saplings. Moreover, species such as Acacia nilotica, Clerodendrum myricoides, Cordia africana, Justicia schimperiana, Phoenix reclinata, and Phytolacca dodecandra were with no seedlings and saplings. A study conducted by Kelbessa and Girma (2011) reported that human-caused disturbances including intensive grazing, fire, and removal of trees are common threats for the regeneration of dryland vegetation. Therefore, species either with no or few seedlings and saplings (Table 6) require conservation priority. Research programs focusing on these species should be a priority as well.

Conclusion
The Hirmi woodland is dominated by undifferentiated floristic composition that belongs to Combretum-Terminalia and Acacia-Commiphora woodland vegetation types with a small section of DAF vegetation types in the elevated areas. Topographic variation of the study area leads to the present heterogeneous vegetation type and structure. In the study vegetation, the dominance of families Fabaceae, Poaceae, Asteraceae, and Lamiaceae agrees with the fact that they are also species-rich families in the country. Furthermore, the presence of 171 species and also 9 endemic and near-endemic species indicates the area's potential of plant diversity and conservation priority. The multivariate analysis displays the identified five plant communities in the study area that were influenced by altitude, slope, total nitrogen, soil organic matter, sand, silt, and disturbance. The high density of species at the lower DBH class and height class is due to the dominance of shrubs or/and the presence of species in the secondary vegetation species such as Dodonaea angustifolia, Flueggea virosa, Diospyros mespiliformis, Ziziphus mucronata, and Maytenus arbutifolia. There is a high degree of floristic heterogeneity, as reflected by the presence of a large number of species in the lower frequency classes. The assessment for the regeneration status of species in the study area showed that the population of saplings was higher than seedlings and mature trees/shrubs. This indicates that the Hirmi woodland has a poor regeneration status which is associated with various disturbance events. Due to the observed multiple disturbance factors (cutting, grazing, fire, and charcoal production signs), Hirmi woodland was found with plant species either with no seedlings, saplings, or both seedlings and sapling. Acacia nilotica, Clerodendrum myricoides, Cordia africana, Justicia schimperiana, Phoenix reclinata, and Phytolacca dodecandra were among the species found under alarming extinction status and with no seedlings and saplings. Therefore, we underscore the immediate need for conservation measures in Hirmi woodland vegetation by considering the identified significant environmental factors associated with the species diversity (altitude, slope, disturbance, soil organic matter, total nitrogen, sand and silt soil type), IUCN threat status and regeneration status of the species.