# Woody Species Diversity and Structure of Protected Woodlands Adjacent to Free Grazing Land Woodland at Dugda Woreda, Oromia, Ethiopia

Authors: Tsegu Ereso Denbel
Journal: International Journal of Forestry Research
 (2021)

---

## Abstract

The study was conducted in protected woodland and free graze woodland located in Dugda Woreda, Oromia state, Ethiopia. The objective of the study was to explore the floristic composition, structure, and regeneration of woody species. In the park, the vegetation ecology has not been studied up to date, which is necessary for conservation. The systematic sampling technique was used to collect vegetation and human disturbance (presence and absence) data from August to December 2017. The vegetation data were collected from 30 plots from each woodland with a size of 900 m2 (30 m × 300 m) for tree/shrub, while subplots of size 400 m2 (20 m × 20 m) for sapling, respectively, were established in the main plots. Individual tree and shrub diameters at breast height (DBH) ≥2.5 cm and height ≥ 2 m were measured using a tape meter and clinometer, respectively. Diameter at breast height (DBH), frequency, density, basal area, and importance value index (IVI) were used for vegetation structure description, while the densities of mature trees, sapling, and seedling were used for regeneration. A total of 446 individual stems from free grazed woodland and 641 individual stems from protected woodland with a DBH of ≥2.5 cm were encountered from 30 studied sample plots that are protected and free grazed woodlands. Of these, from the total woody species, 68.42% were trees and 31.57% shrubs found in protected woodland; 76.92% were trees and 23.07% shrubs found in free grazed woodland. The total basal area of the woody plant was 3.1 ± 1 m2/ha in free grazed woodland and 4.2 ± 2 m2/ha in protected woodland, calculated for 19 woody species. Fabaceae, Balanitaceae, Capparidaceae, Verbenaceae, and Boraginaceae families were the most abundant families in both woodlands. However, there is a good initiation for the conservation of the park; still, the vegetation of the park was threatened by human-induced fire following intensive farming, gold mining, and overgrazing.


---

## Body

## 1. Introduction

Globally, 40% of tropical forest areas and 14% of total African surface area have been covered by woodland [1]. Woodland is defined as frost-free regions with 500–2000 mm annual precipitation and a dry season of four to seven months [2]; its diverse species composition, leafing phenology, and other functional aspects were studied. Next to forest resources, protected woodland has a huge wealth of biological resources [3]; thus, the woody diversity and vegetation structures are crucial elements to clearly visualize the human activities and environmental factors affecting the vegetation of a given area [4]. However, about 97% of the remaining woodland area is at the risk of being further fragmented and degraded by fire, fuelwood collection [5], and grazing damage [6]. Above half of Ethiopia’s land surface is located in dry areas and associated with woodland [7]. These woodlands are a highly fragile ecosystem adapted to semiarid conditions with erratic rainfall, growing on a complex and vulnerable hydrological system [8]. Most woodlands of the world are concentrated in protected areas [9]. Recently, the depletion of vegetation resources in woodland and environmental degradation has become an issue of national and global concern [10]. Dugda Woreda woodland is found in Girba Korke Adi kebele, East Shewa zone, Ethiopia. The woodland benefited the local community by providing pasture, fuelwood, medicine, and honey and served as a habitat for different biodiversity. However, due to habitat degradation, livestock encroachment, illegal logging, and agricultural expansion, and lacking of scientific and essential baseline investigation on woody species diversity, these woodland areas are being imposed to being destroyed at the highest rate. Therefore, the main objective of this study was to complete the current information gap in woody species diversity and structure of the woodland.


## 2. Materials and Methods

### 2.1. Study Area

#### 2.1.1. Description of the Study Area

The study area, protected woodlands at Dugda Woreda, East Shewa, Ethiopia (Figure1), extends from 8°6′00″ to 8°10′00″°N latitude and 38°42′30″ to 38°51′30″°E longitude. The capital city of the Woreda, Meki, is located 132 km south of Addis Ababa, the capital of Ethiopia. The study area is characterized by varied land features, including mountains, flatland, and agricultural fields. The altitude varies from 1600 to 2020 meters masl. Dugda Woreda has 2 major traditional agroclimatic zones, 45% Woyna Dega (midhighland), and 55% Kola (low land). Temperature and rainfall data were collected from the National Meteorological Agency (NMA). The mean annual temperature of the study area was 12.7°C. The mean maximum and minimum temperatures recorded in March and November were 28°C. The area receives a mean annual rainfall of 750 mm with a monomodal pattern, which falls much between October and November. Dugda Woreda has a total population of 163099 people; of these, 85478 were males and 77621 were females. According to [18], a survey of land in this Woreda shows that 65.25% of its land is cultivated, 8.32% is forest, 3.55% is pasture, 12.54% is water bodies, and 0.31% and 10.03% are swampy and rocky mountain areas [11]. Dugda Woreda protected woodland communities are broadly categorized as Acacia-Commiphora woodland with dominant Acacia tortilis, Vachellia gerrardii, and Balanites aegyptiaca species. Thus, the decline in severity and vegetation cover is higher in these lowland protected areas because they are remote and have a scarcity of resources [12].Figure 1
Map of the study area.


### 2.2. Sampling Design

A reconnaissance survey was undertaken for this study from the first week of November 2017 to have an impression of the woodland sites and vegetation distribution and to decide the reasonable sampling methods. During the surveying period, supportive information was collected from the woreda administrative office and local households living near the woodland. A systematic sampling method was used for collecting vegetation data. For protected and open graze land woodland, 30 plots with 30 m × 30 m (900 m2) quadrants for woody species and 10 m × 10 m (100 m2) at two opposite corner of the plot were laid out for collecting shrub data at every 100 m along 12 parallel transect lines. All transects and plots were located using a compass and Geographical Positioning System (GPS) of a navigation system.


### 2.3. Data Collection

Woody vegetation: the detailed vegetation data were collected during the flowering and fruiting season from August 26 to December 30, 2017. Trees and shrubs: in the main plot (600 m2) area, each plant (stem) of all tree and shrub species with a diameter at breast height (DBH) ≥2.5 cm abundance was counted and their circumference was recorded (diameter). Plant species outside plots were also recorded to give a complete list of species. Height of individual trees and shrubs ≥2 m was recorded for every individual woody plant having a diameter at breast height (DBH) ≥2.5 cm [13]. Diameter at breast height (DBH) was measured using a tree caliper and tape meter, while height was measured using a clinometer and visual estimation. Trees with multiple stems arising from the ground level were measured individually and developed a common diameter at breast height (DBH) of all stems by summing up their square roots following the guidelines in [14]. Anthropogenic disturbance: in addition to recording the physiographic variables (altitude, latitude, and longitude), human disturbances such as grazing intensity, illegal logging, and collection firewood signs were recorded in 30 sampled plots. Based on the reconnaissance survey, grazing intensity, illegal logging frequency, and agricultural land explanation were identified as the most pronounced negative impacts in the protected woodland and open graze land; the degree of impact was estimated based on the qualitative presence and absence method./us. If any disturbance sign was present in the plot, it was scaled as presence = one; while if there was no disturbance sign in the plots, it was scaled as absence = zero.Plant species identification was started in the field by recording the local name by asking local elders and referring to the scientific name identification using Flora of Ethiopia and Eritrea Volume 1–Volume 8 [15–21]. Specimens of identified and unidentified species were collected, pressed, and dried properly, following standard Ethiopian Herbarium (ETH) procedures, and taken to Addis Ababa University Herbarium for further confirmation and identification of specimens of those species, which could not be identified in the field.


### 2.4. Data Analysis

#### 2.4.1. Structural Data Analysis

Diameter at breast height, height, basal area, tree density, frequency, and importance value index were used to describe the woody vegetation structure of a given forest. The following formulas were utilized in the Microsoft Excel spreadsheet program and presented in descriptive statistics.Diameter at breast height (DBH): diameter of woody species was arbitrarily arranged in diameter class intervals by referring to the study of Atsbha et al. [21]. The diameter of trees ̸shrubs of plant species of Dugda Woreda protected woodland was classified into nine classes of 7 cm interval. 0 = <2 cm; 1 = 2–5 cm; 2 = 5–10 cm; 3 = 10–15 cm; 4 = 15–20 cm; 5 = 20–25 cm; 6 = 25–30 cm; 7 = >30 cm.Height: height of individuals’ trees ̸shrub plant species was arbitrarily defined by height class intervals [29]. The height of the park was classified into seven classes of 5 m intervals: ≤4, 4.1–9, 9.1–14, 14.1–19, 19.1–24, 24.1–29, and >29.1 m. The densities of individuals falling in the diameter at breast height (DBH) or height classes were summed up. Frequency of species: the probability of finding a species in a given sample area was as follows [22]:(1)frequencyF=number of plots in which a species occurstotal number of plots laid out in the study site∗100,relative frequencyF=frequency of a single speciestotal frequency of all species∗100.Finally, the frequency was summarized by class interval following the study of Lamprecht [23]. The frequency of the park was arranged into seven classes of 15% intervals: ≤5, 5.1–20, 20.1–35, 35.1–50, 50.1–65, 65.1–80, and >80.1%.Density of species: it is a count of the numbers of individuals of each species within the quadrat [30]. The sum of individuals per species is analyzed in terms of species density ha−1 [24]:(2)densityD=number  of  aboveground  stems  of  a  speciesnumber  of  quadrats∗quadrat  area∗100,relative densityRD=density of a single speciestotal density of all species∗100.Density was arranged by class intervals following the dry forest study of Atsbha et al. [21]. The protected woodlands were classified into five density class intervals: ≤2, 2.1–10, 10.1–50, 50.1–100, and >100.1 stems·ha−1.Basal area: it is the area outline of a plant near the ground surface and expressed in m2·ha−1 [22]:(3)basal area=πd24,where π = 3.14 and d is diameter at breast height (m).Dominance: the degree coverage of species occupied a space at ground level was as follows [24]:(4)dominance=the mean basal area per species∗abundancenoof the species,relative dominanceRDO=dominance of a single speciestotal  dominance of all species∗100.Importance value index (IVI): it indicates the relative ecological importance of a given woody species on a particular site [22]:(5)IVI=relative densityRD+relative frequencyRF+relative dominanceRDO.


#### 2.4.2. Statistical Analysis

Then human disturbance variables (Illegal logging, explanation agricultural land, charcoal production, and grazing intensity) were analyzed using regression and correlation statistical method. All analyses were facilitated using the R-statistical package [25]. For qualitative analysis, descriptive statistics were used and these descriptive statistics graphs were performed with the Microsoft Office Excel 2007 software.


## 2.1. Study Area

### 2.1.1. Description of the Study Area

The study area, protected woodlands at Dugda Woreda, East Shewa, Ethiopia (Figure1), extends from 8°6′00″ to 8°10′00″°N latitude and 38°42′30″ to 38°51′30″°E longitude. The capital city of the Woreda, Meki, is located 132 km south of Addis Ababa, the capital of Ethiopia. The study area is characterized by varied land features, including mountains, flatland, and agricultural fields. The altitude varies from 1600 to 2020 meters masl. Dugda Woreda has 2 major traditional agroclimatic zones, 45% Woyna Dega (midhighland), and 55% Kola (low land). Temperature and rainfall data were collected from the National Meteorological Agency (NMA). The mean annual temperature of the study area was 12.7°C. The mean maximum and minimum temperatures recorded in March and November were 28°C. The area receives a mean annual rainfall of 750 mm with a monomodal pattern, which falls much between October and November. Dugda Woreda has a total population of 163099 people; of these, 85478 were males and 77621 were females. According to [18], a survey of land in this Woreda shows that 65.25% of its land is cultivated, 8.32% is forest, 3.55% is pasture, 12.54% is water bodies, and 0.31% and 10.03% are swampy and rocky mountain areas [11]. Dugda Woreda protected woodland communities are broadly categorized as Acacia-Commiphora woodland with dominant Acacia tortilis, Vachellia gerrardii, and Balanites aegyptiaca species. Thus, the decline in severity and vegetation cover is higher in these lowland protected areas because they are remote and have a scarcity of resources [12].Figure 1
Map of the study area.


## 2.1.1. Description of the Study Area

The study area, protected woodlands at Dugda Woreda, East Shewa, Ethiopia (Figure1), extends from 8°6′00″ to 8°10′00″°N latitude and 38°42′30″ to 38°51′30″°E longitude. The capital city of the Woreda, Meki, is located 132 km south of Addis Ababa, the capital of Ethiopia. The study area is characterized by varied land features, including mountains, flatland, and agricultural fields. The altitude varies from 1600 to 2020 meters masl. Dugda Woreda has 2 major traditional agroclimatic zones, 45% Woyna Dega (midhighland), and 55% Kola (low land). Temperature and rainfall data were collected from the National Meteorological Agency (NMA). The mean annual temperature of the study area was 12.7°C. The mean maximum and minimum temperatures recorded in March and November were 28°C. The area receives a mean annual rainfall of 750 mm with a monomodal pattern, which falls much between October and November. Dugda Woreda has a total population of 163099 people; of these, 85478 were males and 77621 were females. According to [18], a survey of land in this Woreda shows that 65.25% of its land is cultivated, 8.32% is forest, 3.55% is pasture, 12.54% is water bodies, and 0.31% and 10.03% are swampy and rocky mountain areas [11]. Dugda Woreda protected woodland communities are broadly categorized as Acacia-Commiphora woodland with dominant Acacia tortilis, Vachellia gerrardii, and Balanites aegyptiaca species. Thus, the decline in severity and vegetation cover is higher in these lowland protected areas because they are remote and have a scarcity of resources [12].Figure 1
Map of the study area.


## 2.2. Sampling Design

A reconnaissance survey was undertaken for this study from the first week of November 2017 to have an impression of the woodland sites and vegetation distribution and to decide the reasonable sampling methods. During the surveying period, supportive information was collected from the woreda administrative office and local households living near the woodland. A systematic sampling method was used for collecting vegetation data. For protected and open graze land woodland, 30 plots with 30 m × 30 m (900 m2) quadrants for woody species and 10 m × 10 m (100 m2) at two opposite corner of the plot were laid out for collecting shrub data at every 100 m along 12 parallel transect lines. All transects and plots were located using a compass and Geographical Positioning System (GPS) of a navigation system.


## 2.3. Data Collection

Woody vegetation: the detailed vegetation data were collected during the flowering and fruiting season from August 26 to December 30, 2017. Trees and shrubs: in the main plot (600 m2) area, each plant (stem) of all tree and shrub species with a diameter at breast height (DBH) ≥2.5 cm abundance was counted and their circumference was recorded (diameter). Plant species outside plots were also recorded to give a complete list of species. Height of individual trees and shrubs ≥2 m was recorded for every individual woody plant having a diameter at breast height (DBH) ≥2.5 cm [13]. Diameter at breast height (DBH) was measured using a tree caliper and tape meter, while height was measured using a clinometer and visual estimation. Trees with multiple stems arising from the ground level were measured individually and developed a common diameter at breast height (DBH) of all stems by summing up their square roots following the guidelines in [14]. Anthropogenic disturbance: in addition to recording the physiographic variables (altitude, latitude, and longitude), human disturbances such as grazing intensity, illegal logging, and collection firewood signs were recorded in 30 sampled plots. Based on the reconnaissance survey, grazing intensity, illegal logging frequency, and agricultural land explanation were identified as the most pronounced negative impacts in the protected woodland and open graze land; the degree of impact was estimated based on the qualitative presence and absence method./us. If any disturbance sign was present in the plot, it was scaled as presence = one; while if there was no disturbance sign in the plots, it was scaled as absence = zero.Plant species identification was started in the field by recording the local name by asking local elders and referring to the scientific name identification using Flora of Ethiopia and Eritrea Volume 1–Volume 8 [15–21]. Specimens of identified and unidentified species were collected, pressed, and dried properly, following standard Ethiopian Herbarium (ETH) procedures, and taken to Addis Ababa University Herbarium for further confirmation and identification of specimens of those species, which could not be identified in the field.


## 2.4. Data Analysis

### 2.4.1. Structural Data Analysis

Diameter at breast height, height, basal area, tree density, frequency, and importance value index were used to describe the woody vegetation structure of a given forest. The following formulas were utilized in the Microsoft Excel spreadsheet program and presented in descriptive statistics.Diameter at breast height (DBH): diameter of woody species was arbitrarily arranged in diameter class intervals by referring to the study of Atsbha et al. [21]. The diameter of trees ̸shrubs of plant species of Dugda Woreda protected woodland was classified into nine classes of 7 cm interval. 0 = <2 cm; 1 = 2–5 cm; 2 = 5–10 cm; 3 = 10–15 cm; 4 = 15–20 cm; 5 = 20–25 cm; 6 = 25–30 cm; 7 = >30 cm.Height: height of individuals’ trees ̸shrub plant species was arbitrarily defined by height class intervals [29]. The height of the park was classified into seven classes of 5 m intervals: ≤4, 4.1–9, 9.1–14, 14.1–19, 19.1–24, 24.1–29, and >29.1 m. The densities of individuals falling in the diameter at breast height (DBH) or height classes were summed up. Frequency of species: the probability of finding a species in a given sample area was as follows [22]:(1)frequencyF=number of plots in which a species occurstotal number of plots laid out in the study site∗100,relative frequencyF=frequency of a single speciestotal frequency of all species∗100.Finally, the frequency was summarized by class interval following the study of Lamprecht [23]. The frequency of the park was arranged into seven classes of 15% intervals: ≤5, 5.1–20, 20.1–35, 35.1–50, 50.1–65, 65.1–80, and >80.1%.Density of species: it is a count of the numbers of individuals of each species within the quadrat [30]. The sum of individuals per species is analyzed in terms of species density ha−1 [24]:(2)densityD=number  of  aboveground  stems  of  a  speciesnumber  of  quadrats∗quadrat  area∗100,relative densityRD=density of a single speciestotal density of all species∗100.Density was arranged by class intervals following the dry forest study of Atsbha et al. [21]. The protected woodlands were classified into five density class intervals: ≤2, 2.1–10, 10.1–50, 50.1–100, and >100.1 stems·ha−1.Basal area: it is the area outline of a plant near the ground surface and expressed in m2·ha−1 [22]:(3)basal area=πd24,where π = 3.14 and d is diameter at breast height (m).Dominance: the degree coverage of species occupied a space at ground level was as follows [24]:(4)dominance=the mean basal area per species∗abundancenoof the species,relative dominanceRDO=dominance of a single speciestotal  dominance of all species∗100.Importance value index (IVI): it indicates the relative ecological importance of a given woody species on a particular site [22]:(5)IVI=relative densityRD+relative frequencyRF+relative dominanceRDO.


### 2.4.2. Statistical Analysis

Then human disturbance variables (Illegal logging, explanation agricultural land, charcoal production, and grazing intensity) were analyzed using regression and correlation statistical method. All analyses were facilitated using the R-statistical package [25]. For qualitative analysis, descriptive statistics were used and these descriptive statistics graphs were performed with the Microsoft Office Excel 2007 software.


## 2.4.1. Structural Data Analysis

Diameter at breast height, height, basal area, tree density, frequency, and importance value index were used to describe the woody vegetation structure of a given forest. The following formulas were utilized in the Microsoft Excel spreadsheet program and presented in descriptive statistics.Diameter at breast height (DBH): diameter of woody species was arbitrarily arranged in diameter class intervals by referring to the study of Atsbha et al. [21]. The diameter of trees ̸shrubs of plant species of Dugda Woreda protected woodland was classified into nine classes of 7 cm interval. 0 = <2 cm; 1 = 2–5 cm; 2 = 5–10 cm; 3 = 10–15 cm; 4 = 15–20 cm; 5 = 20–25 cm; 6 = 25–30 cm; 7 = >30 cm.Height: height of individuals’ trees ̸shrub plant species was arbitrarily defined by height class intervals [29]. The height of the park was classified into seven classes of 5 m intervals: ≤4, 4.1–9, 9.1–14, 14.1–19, 19.1–24, 24.1–29, and >29.1 m. The densities of individuals falling in the diameter at breast height (DBH) or height classes were summed up. Frequency of species: the probability of finding a species in a given sample area was as follows [22]:(1)frequencyF=number of plots in which a species occurstotal number of plots laid out in the study site∗100,relative frequencyF=frequency of a single speciestotal frequency of all species∗100.Finally, the frequency was summarized by class interval following the study of Lamprecht [23]. The frequency of the park was arranged into seven classes of 15% intervals: ≤5, 5.1–20, 20.1–35, 35.1–50, 50.1–65, 65.1–80, and >80.1%.Density of species: it is a count of the numbers of individuals of each species within the quadrat [30]. The sum of individuals per species is analyzed in terms of species density ha−1 [24]:(2)densityD=number  of  aboveground  stems  of  a  speciesnumber  of  quadrats∗quadrat  area∗100,relative densityRD=density of a single speciestotal density of all species∗100.Density was arranged by class intervals following the dry forest study of Atsbha et al. [21]. The protected woodlands were classified into five density class intervals: ≤2, 2.1–10, 10.1–50, 50.1–100, and >100.1 stems·ha−1.Basal area: it is the area outline of a plant near the ground surface and expressed in m2·ha−1 [22]:(3)basal area=πd24,where π = 3.14 and d is diameter at breast height (m).Dominance: the degree coverage of species occupied a space at ground level was as follows [24]:(4)dominance=the mean basal area per species∗abundancenoof the species,relative dominanceRDO=dominance of a single speciestotal  dominance of all species∗100.Importance value index (IVI): it indicates the relative ecological importance of a given woody species on a particular site [22]:(5)IVI=relative densityRD+relative frequencyRF+relative dominanceRDO.


## 2.4.2. Statistical Analysis

Then human disturbance variables (Illegal logging, explanation agricultural land, charcoal production, and grazing intensity) were analyzed using regression and correlation statistical method. All analyses were facilitated using the R-statistical package [25]. For qualitative analysis, descriptive statistics were used and these descriptive statistics graphs were performed with the Microsoft Office Excel 2007 software.


## 3. Results

### 3.1. Species Composition and Diversity

A total of 446 individual stems from free grazed woodland and 641 individual stems from protected woodland with a diameter at breast height (DBH) of ≥2.5 cm were collected. A total of 19 woody species were identified from 30 studied sample plots from protected and free grazed woodlands. Out of these, 18 woody species belonging to 17 genera and 7 families were identified from protected woodlands and 13 woody species belonging to 14 genera and 4 families were recorded from free grazed woodland (Table1).Table 1
Species composition and diversity.
Woodland typeNumber of speciesNumber of generaNumber of familiesFree grazed woodland446144Protected woodland641177Twelve species are commonly found in both sites. Fabaceae, Balanitaceae, Capparidaceae, Verbenaceae, and Boraginaceae families were the most abundant families in both woodlands. From the total woody species, 68.42% were trees and 31.57% shrubs found in protected woodland; 76.92% were trees and 23.07% shrubs found in free grazed woodland (Figure2).Figure 2
Floristic composition of protected and free grazed woodland.


### 3.2. Diversity of Woody Species

The mean of the Shannon–Wiener diversity index of woody species is 1.69 for free grazed woodland and 2.26 for protected woodland (Table2). The Shannon–Wiener diversity index is higher at protected woodland than free grazed woodland. The Simpson diversity index and evenness of woody species in protected woodland are 0.82 and 0.51, respectively. On the other hand, the Simpson diversity index and evenness of woody species in free protected woodland are 0.86 and 1.69, respectively. The species richness between the two study sites is also slightly different, with protected areas having higher (8.6) species richness than free gazed woodland (7.28) (Table 2).Table 2
Mean diversity indices and richness of woody species in protected and free grazed woodland.
Types of woodlandShannonEvennessSimpsonSpecies richnessFree grazed woodland1.69 ± 0.77b0.74a ± 0.250.86a ± 0.257.28a ± 3.23Protected woodland2.26 ± 0.67a0.51b ± 0.210.82a ± 0.218.63b ± 4.23Grand mean1.97 ± 0.720.62 ± 0.230.73 ± 0.237.96 ± 5.52Different superscript small letters indicate there is significant variation between the two woodland categories (F test, P < 0.05) with Shannon, evenness, and richness of woodland.


### 3.3. Vegetation Structure of Woodland

Except for the mean of the height and diameter branch, the other structural parameters in the studied woodland varied significantly (P < 0.05) (Table 3). The height and branch of the tree are not varied because the management of the woodland and environmental factors influence the growth of the species more or less the same.Table 3
Stand structure of protected and free grazed woodland (mean ± SE).
Stand charactersFree grazed woodlandProtected woodlandOverall meanP valueMean dbh (cm)53.2 ± 7a67.3 ± 9b60.25 ± 80.036Mean diameter branch0.15 ± 2a0.28.1 ± 3a0.43 ± 1.50.23Mean height (m)6.7 ± 2a8.2 ± 1a7.45 ± 1.50.41Mean crown ratio33.2 ± 4.1b42.63 ± 4.02b37.9 ± 4.020.003Mean crown diameter5.4 ± 1a6.7 ± 2b6.02 ± 1.50.001Density ha−1317.24 ± 52.6a346.53 ± 57.3b331.9 ± 54.950.002Basal area (m2 ha−1)3.1 ± 1a4.2 ± 2b3.45 ± 1.50.0024Volume (m3 ha−1)235.17 ± 12.3a355.06 ± 23.6b295.12 ± 17.950.01Canopy cover6.4 ± 1a7.3 ± 2b6.85 ± 1.50.0006Different superscript small letters indicate there is significant variation among the two woodland categories (F test, P < 0.05) with stand structure between woodlands.


### 3.4. Woody Species Frequency

The frequency occurrence of woody species in protected woodland and free grazed woodland is different in Figures3 and 4. In protected woodland, Acacia tortilis, Dichrostachys cinerea, Faidherbia albida, Acacia senegal, Acacia gerrardii, and Balanites aegyptiaca were among the most frequent species frequently. On the other hand, Carissa spinarum, Acacia oerfota, and Rhus vulgaris were among the least recorded woody species (Figure 3).Figure 3
Frequency of occurrence of woody species in protected woodland.Figure 4
Frequency of occurrences of woody species in free grazed woodland.Similarly,Acacia tortilis, Balanites aegyptiaca, Acacia seyal, Acacia amythethophylla, Acacia gerrardii, and Dichrostachys cinerea were commonly found in free grazed woodland (Figure 4). The frequency of woody species in woodland is determined by human intervention and ecological factor.


### 3.5. Importance Value Index

The dominant and ecologically most significant tree species in free grazed woodland wereAcacia tortilis, Dichrostachys cinerea, Acacia amythethophylla, and Acacia gerrardii (Table 4) and in protected woodland were Acacia tortilis, Dichrostachys cinerea, Acacia gerrardii, Balanites aegyptiaca, Acacia Senegal, and Faidherbia albida (Table 5). Acacia tortilis was the most dominant and ecologically significant in both woodlands.Table 4
Summary of IVI of woody species in free grazed woodland in descending order.
Free grazed woodlandSpecies nameFamilyRFRDORDIVIAcacia tortilisFabaceae23.7113495.7618643.6377163.1109Acacia gerrardiiFabaceae9.2783510.9008618.06747318.24668Dichrostachys cinereaFabaceae7.2164950.6133339.71763817.54747Acacia amythethophyllaFabaceae8.2474231.0869086.78401216.11834Balanites aegyptiacaFabaceae11.340210.0587524.03373715.43269Acacia seyalFabaceae9.2783510.4281762.38357212.0901Faidherbia albidaFabaceae6.1855670.0096144.95049511.14568Acacia senegalFabaceae4.1237110.6053216.23395710.96299Acacia niloticaFabaceae6.1855670.0890184.58379210.85838Acacia asakFabaceae4.1237110.0115722.3835726.518855Acacia decurrensFabaceae3.0927840.059822.9336276.08623Acacia brevispicaFabaceae4.1237110.0009351.2834625.408108Ehretia cymosaBoraginaceae2.0618560.3738752.5669235.002654Capparis cartilagineaCapparidaceae1.0309281.87E − 060.5500551.580985Total100100100300RF = relative frequency; RDO = relative dominance; RD = relative density; IVI = important value index.Table 5
Summary of IVI value of woody species in protected woodland in descending order.
Protected woodlandScientific nameFamily nameRFRDORDIVIAcacia tortilisFabaceae9.7902185.8402831.7161127.346Dichrostachys cinereaFabaceae13.286717.3123815.073035.6721Acacia gerrardiiFabaceae7.6923084.02902510.519722.2410Acacia senegalFabaceae8.3916080.4534166.1234114.9684Balanites aegyptiacaBalanitaceae7.6923080.7745356.5944415.0612Faidherbia albidaFabaceae9.0909090.2886475.652315.0319Acacia seyalFabaceae6.9930070.0875564.0822711.1628Acacia amythethophyllaFabaceae5.5944060.2155234.3962910.2062Acacia niloticaFabaceae5.5944060.0902023.925269.60987Acacia decurrensFabaceae4.1958040.8819763.454238.53201Duranta erectaVerbenaceae5.5944060.013552.041137.64909Acacia asakFabaceae3.4965030.0036561.884125.38428Ehretia cymosaBoraginaceae3.4965030.0020211.413094.91161Carissa spinarumApocynaceae2.7972030.0006250.942063.73989Acacia oerfotaFabaceae2.7972030.0003370.942063.73960Capparis cartilagineaCapparidaceae2.7972030.0012430.785053.58349Rhus vulgarisAnacardiaceae0.6993010.0009620.628041.32830Total100100100300RF = relative frequency; RDO = relative dominance; RD = relative density; IVI = important value index.


### 3.6. Size Distribution

The first pattern was represented byAcacia tortilis (Figure 5(a)). This pattern indicates the presence of the highest density in the lower DBH classes, with a gradual decrease in density towards the bigger classes. It represents an inverted J-shaped curve except for a slight decline in the first DBH class, and it suggests good reproduction and recruitment.Figure 5
DBH class distribution of some important woody species in free grazed woodland. 0 = <2 cm; 1 = 2–5 cm; 2 = 5–10 cm; 3 = 10–15 cm; 4 = 15–20 cm; 5 = 20–25 cm; 6 = 25–30 cm; 7 = >30 cm. (a)Acacia tortilis in free gazed woodland. (b) Dichrostachys cinerea in free grazed woodland. (c) Acacia gerrardii in free grazed woodland.
(a)(b)(c)The second pattern is represented byDichrostachys cinerea (Figure 5(b)). This pattern indicates the presence of the highest density in the middle DBH classes, with a gradual decrease in density towards the bigger classes. It represents an inverted J-shaped curve except for a slight decline in the first DBH class, and it suggests good reproduction and recruitment.The third pattern was represented by a single species ofAcacia gerrardii (Figure 5(c)). This pattern indicates the presence of the lowest density in the lower DBH classes with a gradual increase in density towards the bigger classes. It represents a J-shaped curve and it suggests there is not good reproduction and less sustainability (Figure 6).Figure 6
DBH class distribution of some important woody species in protected woodlands. 0 = <2 cm; 1 = 2–5 cm; 2 = 5–10 cm; 3 = 10–15 cm; 4 = 15–20 cm; 5 = 20–25 cm; 6 = 25–30 cm; 7 = >30 cm. (a)Acacia tortilis in PW. (b) Dichrostachys cinerea. (c) Acacia gerrardii.
(a)(b)(c)


## 3.1. Species Composition and Diversity

A total of 446 individual stems from free grazed woodland and 641 individual stems from protected woodland with a diameter at breast height (DBH) of ≥2.5 cm were collected. A total of 19 woody species were identified from 30 studied sample plots from protected and free grazed woodlands. Out of these, 18 woody species belonging to 17 genera and 7 families were identified from protected woodlands and 13 woody species belonging to 14 genera and 4 families were recorded from free grazed woodland (Table1).Table 1
Species composition and diversity.
Woodland typeNumber of speciesNumber of generaNumber of familiesFree grazed woodland446144Protected woodland641177Twelve species are commonly found in both sites. Fabaceae, Balanitaceae, Capparidaceae, Verbenaceae, and Boraginaceae families were the most abundant families in both woodlands. From the total woody species, 68.42% were trees and 31.57% shrubs found in protected woodland; 76.92% were trees and 23.07% shrubs found in free grazed woodland (Figure2).Figure 2
Floristic composition of protected and free grazed woodland.


## 3.2. Diversity of Woody Species

The mean of the Shannon–Wiener diversity index of woody species is 1.69 for free grazed woodland and 2.26 for protected woodland (Table2). The Shannon–Wiener diversity index is higher at protected woodland than free grazed woodland. The Simpson diversity index and evenness of woody species in protected woodland are 0.82 and 0.51, respectively. On the other hand, the Simpson diversity index and evenness of woody species in free protected woodland are 0.86 and 1.69, respectively. The species richness between the two study sites is also slightly different, with protected areas having higher (8.6) species richness than free gazed woodland (7.28) (Table 2).Table 2
Mean diversity indices and richness of woody species in protected and free grazed woodland.
Types of woodlandShannonEvennessSimpsonSpecies richnessFree grazed woodland1.69 ± 0.77b0.74a ± 0.250.86a ± 0.257.28a ± 3.23Protected woodland2.26 ± 0.67a0.51b ± 0.210.82a ± 0.218.63b ± 4.23Grand mean1.97 ± 0.720.62 ± 0.230.73 ± 0.237.96 ± 5.52Different superscript small letters indicate there is significant variation between the two woodland categories (F test, P < 0.05) with Shannon, evenness, and richness of woodland.


## 3.3. Vegetation Structure of Woodland

Except for the mean of the height and diameter branch, the other structural parameters in the studied woodland varied significantly (P < 0.05) (Table 3). The height and branch of the tree are not varied because the management of the woodland and environmental factors influence the growth of the species more or less the same.Table 3
Stand structure of protected and free grazed woodland (mean ± SE).
Stand charactersFree grazed woodlandProtected woodlandOverall meanP valueMean dbh (cm)53.2 ± 7a67.3 ± 9b60.25 ± 80.036Mean diameter branch0.15 ± 2a0.28.1 ± 3a0.43 ± 1.50.23Mean height (m)6.7 ± 2a8.2 ± 1a7.45 ± 1.50.41Mean crown ratio33.2 ± 4.1b42.63 ± 4.02b37.9 ± 4.020.003Mean crown diameter5.4 ± 1a6.7 ± 2b6.02 ± 1.50.001Density ha−1317.24 ± 52.6a346.53 ± 57.3b331.9 ± 54.950.002Basal area (m2 ha−1)3.1 ± 1a4.2 ± 2b3.45 ± 1.50.0024Volume (m3 ha−1)235.17 ± 12.3a355.06 ± 23.6b295.12 ± 17.950.01Canopy cover6.4 ± 1a7.3 ± 2b6.85 ± 1.50.0006Different superscript small letters indicate there is significant variation among the two woodland categories (F test, P < 0.05) with stand structure between woodlands.


## 3.4. Woody Species Frequency

The frequency occurrence of woody species in protected woodland and free grazed woodland is different in Figures3 and 4. In protected woodland, Acacia tortilis, Dichrostachys cinerea, Faidherbia albida, Acacia senegal, Acacia gerrardii, and Balanites aegyptiaca were among the most frequent species frequently. On the other hand, Carissa spinarum, Acacia oerfota, and Rhus vulgaris were among the least recorded woody species (Figure 3).Figure 3
Frequency of occurrence of woody species in protected woodland.Figure 4
Frequency of occurrences of woody species in free grazed woodland.Similarly,Acacia tortilis, Balanites aegyptiaca, Acacia seyal, Acacia amythethophylla, Acacia gerrardii, and Dichrostachys cinerea were commonly found in free grazed woodland (Figure 4). The frequency of woody species in woodland is determined by human intervention and ecological factor.


## 3.5. Importance Value Index

The dominant and ecologically most significant tree species in free grazed woodland wereAcacia tortilis, Dichrostachys cinerea, Acacia amythethophylla, and Acacia gerrardii (Table 4) and in protected woodland were Acacia tortilis, Dichrostachys cinerea, Acacia gerrardii, Balanites aegyptiaca, Acacia Senegal, and Faidherbia albida (Table 5). Acacia tortilis was the most dominant and ecologically significant in both woodlands.Table 4
Summary of IVI of woody species in free grazed woodland in descending order.
Free grazed woodlandSpecies nameFamilyRFRDORDIVIAcacia tortilisFabaceae23.7113495.7618643.6377163.1109Acacia gerrardiiFabaceae9.2783510.9008618.06747318.24668Dichrostachys cinereaFabaceae7.2164950.6133339.71763817.54747Acacia amythethophyllaFabaceae8.2474231.0869086.78401216.11834Balanites aegyptiacaFabaceae11.340210.0587524.03373715.43269Acacia seyalFabaceae9.2783510.4281762.38357212.0901Faidherbia albidaFabaceae6.1855670.0096144.95049511.14568Acacia senegalFabaceae4.1237110.6053216.23395710.96299Acacia niloticaFabaceae6.1855670.0890184.58379210.85838Acacia asakFabaceae4.1237110.0115722.3835726.518855Acacia decurrensFabaceae3.0927840.059822.9336276.08623Acacia brevispicaFabaceae4.1237110.0009351.2834625.408108Ehretia cymosaBoraginaceae2.0618560.3738752.5669235.002654Capparis cartilagineaCapparidaceae1.0309281.87E − 060.5500551.580985Total100100100300RF = relative frequency; RDO = relative dominance; RD = relative density; IVI = important value index.Table 5
Summary of IVI value of woody species in protected woodland in descending order.
Protected woodlandScientific nameFamily nameRFRDORDIVIAcacia tortilisFabaceae9.7902185.8402831.7161127.346Dichrostachys cinereaFabaceae13.286717.3123815.073035.6721Acacia gerrardiiFabaceae7.6923084.02902510.519722.2410Acacia senegalFabaceae8.3916080.4534166.1234114.9684Balanites aegyptiacaBalanitaceae7.6923080.7745356.5944415.0612Faidherbia albidaFabaceae9.0909090.2886475.652315.0319Acacia seyalFabaceae6.9930070.0875564.0822711.1628Acacia amythethophyllaFabaceae5.5944060.2155234.3962910.2062Acacia niloticaFabaceae5.5944060.0902023.925269.60987Acacia decurrensFabaceae4.1958040.8819763.454238.53201Duranta erectaVerbenaceae5.5944060.013552.041137.64909Acacia asakFabaceae3.4965030.0036561.884125.38428Ehretia cymosaBoraginaceae3.4965030.0020211.413094.91161Carissa spinarumApocynaceae2.7972030.0006250.942063.73989Acacia oerfotaFabaceae2.7972030.0003370.942063.73960Capparis cartilagineaCapparidaceae2.7972030.0012430.785053.58349Rhus vulgarisAnacardiaceae0.6993010.0009620.628041.32830Total100100100300RF = relative frequency; RDO = relative dominance; RD = relative density; IVI = important value index.


## 3.6. Size Distribution

The first pattern was represented byAcacia tortilis (Figure 5(a)). This pattern indicates the presence of the highest density in the lower DBH classes, with a gradual decrease in density towards the bigger classes. It represents an inverted J-shaped curve except for a slight decline in the first DBH class, and it suggests good reproduction and recruitment.Figure 5
DBH class distribution of some important woody species in free grazed woodland. 0 = <2 cm; 1 = 2–5 cm; 2 = 5–10 cm; 3 = 10–15 cm; 4 = 15–20 cm; 5 = 20–25 cm; 6 = 25–30 cm; 7 = >30 cm. (a)Acacia tortilis in free gazed woodland. (b) Dichrostachys cinerea in free grazed woodland. (c) Acacia gerrardii in free grazed woodland.
(a)(b)(c)The second pattern is represented byDichrostachys cinerea (Figure 5(b)). This pattern indicates the presence of the highest density in the middle DBH classes, with a gradual decrease in density towards the bigger classes. It represents an inverted J-shaped curve except for a slight decline in the first DBH class, and it suggests good reproduction and recruitment.The third pattern was represented by a single species ofAcacia gerrardii (Figure 5(c)). This pattern indicates the presence of the lowest density in the lower DBH classes with a gradual increase in density towards the bigger classes. It represents a J-shaped curve and it suggests there is not good reproduction and less sustainability (Figure 6).Figure 6
DBH class distribution of some important woody species in protected woodlands. 0 = <2 cm; 1 = 2–5 cm; 2 = 5–10 cm; 3 = 10–15 cm; 4 = 15–20 cm; 5 = 20–25 cm; 6 = 25–30 cm; 7 = >30 cm. (a)Acacia tortilis in PW. (b) Dichrostachys cinerea. (c) Acacia gerrardii.
(a)(b)(c)


## 4. Discussion

### 4.1. Species Composition and Diversity

In the present study, we found that the species composition differed between the protected woodland site and the free grazed woodland site. This difference in species composition might be associated with differences in grazing time and regeneration characteristics of individual tree and soil characteristics between the two sites [26]. The results from analyzed data showed that Fabaceae and Balanitaceae are the most dominant species in Giraba Korke Adii kebele woodland. Vachellia tortilis in free grazed woodland declined as a result of illegal cutting and collecting for fuelwood and Vachellia gerrardii found were highly declined in both woodlands as a result of less regeneration, illegal cutting, and harvest for fuelwood during mourning time in kebele. Those dominant species were gradually replaced with other species year to year, especially with Dichrostachys cinerea. Such changes in species composition have also been observed elsewhere, particularly in relation to uncontrolled grazing [26].


### 4.2. Species Diversity, Evenness, and Richness

The overall diversity of woody plants was much higher in protected woodland (H’ = 2.28) than free grazed woodland (H’ = 1.692), which may be a consequence of the high species richness in protected woodland. It has been noted that the value of obtained empirical data usually falls between 1.5 and 3.5. This implies that the diversity of woody species at free grazed woodland falls at the lowest value of the diversity range. Cattle intervention and woodcutting during weeding and mourning may be responsible for the considerable floristic variations (richness of species, genera, and families and overall diversity) between the two sites.


### 4.3. Population Structure and Important Value Index

Based on the assessment of diameter class distributions, the population structure patterns of the woody species recorded from protected and free grazed were categorized into six groups.Vachellia tortilis and Dichrostachys cinerea species in the protected woodland show inverted J-shaped distribution, which is widely acknowledged to indicate stable population structure, naturally replacing individuals of all age categories together with seedlings and saplings [27]. The results of the size distribution of protected woodland showed a similar trend, which was also reported from the semiarid areas of northern Ethiopia [28, 29]. This indicates a good status of regeneration and stable population of the vegetation in Dugda Woreda with good opportunity for sustainable management. On the other hand, the result of the size distribution of free grazed woodland does not exhibit inverse J-shaped. Tree density is expressed as the number of trees per unit area and it is a crucial parameter for sustainable woodland management. The total density of all woody species was found to be 317.24 and 346.53 stems ha−1 at free gazed woodland and protected woodland use, respectively (Table 2). The total basal area of the woody plant was 3.1 ± 1 m2/ha in free grazed woodland and 4.2 ± 2 m2/ha in protected woodland (Table 2). Importance value index (IVI) is a good measure for summarizing vegetation characteristics of a given habitat and is useful to compare the ecological significance of species and for conservation practices. Vachellia tortilis, Vachellia gerrardii, Dichrostachys cinerea, and Vachellia amythetophylla at free grazed woodland and Vachellia tortilis, Dichrostachys cinerea, Vachellia gerrardii, and Senegalia senegal at protected woodland can be considered the most ecologically important woody species with IVI values based on the contributed high values of density, frequency, and dominance. It is interesting to note that Vachellia tortilis was recorded among the most ecologically important woody species at both study sites.


## 4.1. Species Composition and Diversity

In the present study, we found that the species composition differed between the protected woodland site and the free grazed woodland site. This difference in species composition might be associated with differences in grazing time and regeneration characteristics of individual tree and soil characteristics between the two sites [26]. The results from analyzed data showed that Fabaceae and Balanitaceae are the most dominant species in Giraba Korke Adii kebele woodland. Vachellia tortilis in free grazed woodland declined as a result of illegal cutting and collecting for fuelwood and Vachellia gerrardii found were highly declined in both woodlands as a result of less regeneration, illegal cutting, and harvest for fuelwood during mourning time in kebele. Those dominant species were gradually replaced with other species year to year, especially with Dichrostachys cinerea. Such changes in species composition have also been observed elsewhere, particularly in relation to uncontrolled grazing [26].


## 4.2. Species Diversity, Evenness, and Richness

The overall diversity of woody plants was much higher in protected woodland (H’ = 2.28) than free grazed woodland (H’ = 1.692), which may be a consequence of the high species richness in protected woodland. It has been noted that the value of obtained empirical data usually falls between 1.5 and 3.5. This implies that the diversity of woody species at free grazed woodland falls at the lowest value of the diversity range. Cattle intervention and woodcutting during weeding and mourning may be responsible for the considerable floristic variations (richness of species, genera, and families and overall diversity) between the two sites.


## 4.3. Population Structure and Important Value Index

Based on the assessment of diameter class distributions, the population structure patterns of the woody species recorded from protected and free grazed were categorized into six groups.Vachellia tortilis and Dichrostachys cinerea species in the protected woodland show inverted J-shaped distribution, which is widely acknowledged to indicate stable population structure, naturally replacing individuals of all age categories together with seedlings and saplings [27]. The results of the size distribution of protected woodland showed a similar trend, which was also reported from the semiarid areas of northern Ethiopia [28, 29]. This indicates a good status of regeneration and stable population of the vegetation in Dugda Woreda with good opportunity for sustainable management. On the other hand, the result of the size distribution of free grazed woodland does not exhibit inverse J-shaped. Tree density is expressed as the number of trees per unit area and it is a crucial parameter for sustainable woodland management. The total density of all woody species was found to be 317.24 and 346.53 stems ha−1 at free gazed woodland and protected woodland use, respectively (Table 2). The total basal area of the woody plant was 3.1 ± 1 m2/ha in free grazed woodland and 4.2 ± 2 m2/ha in protected woodland (Table 2). Importance value index (IVI) is a good measure for summarizing vegetation characteristics of a given habitat and is useful to compare the ecological significance of species and for conservation practices. Vachellia tortilis, Vachellia gerrardii, Dichrostachys cinerea, and Vachellia amythetophylla at free grazed woodland and Vachellia tortilis, Dichrostachys cinerea, Vachellia gerrardii, and Senegalia senegal at protected woodland can be considered the most ecologically important woody species with IVI values based on the contributed high values of density, frequency, and dominance. It is interesting to note that Vachellia tortilis was recorded among the most ecologically important woody species at both study sites.


## 5. Conclusion

The Dugda Woreda woodland has a small number of woody species bound with small diversity. The woodland was dominated by woody species, which was dominated by the Fabaceae family. The renewal of species through regeneration was weak; the vulnerability of young plants slowed down the sustainability of woody species. However, human pressures on the majority of the species are represented by woody cuts during death mourning for fuelwood, overgrazing in free grazed woodland and camel intervention are considered as ecological and environmental problems, contributing to the degradation of the diversity in Dugda Woreda woodland.The present study has provided valuable information on woody species diversity, structure, and sustainability status of woody plants in the protected and free grazed woodland that would be used as an input for effective conservation of the area. The variation in population structure and regeneration status of the park indicates the long-time past disturbance of individual species and the whole resources of the park.The importance value index (IVI) values revealed that there is poor regeneration status in woodland due to the activities of human disturbance variables, particularly livestock browsing and grazing and human-induced fire following extensive cultivation and tradition and expansion of agriculture. These factors lead to a decline in the density of seedling and sapling of common tree species, which are economically and ecologically important in the woodland.



---

*Source: 8834892-2021-06-22.xml*