Analysis of Vegetation Phytosociological Characteristics and Soil Physico-Chemical Conditions in Harishin Rangelands of Eastern Ethiopia

The objective of this study is to analyse the phytodiversity, distribution, herb biomass and physico-chemical conditions of the vegetation system in the context of communal continuous open grazing and enclosed grazing management practices in the Harishin rangelands of Eastern Ethiopia. A total of 58 herbaceous species and 11 woody species were recorded in the study area. Analysis of Importance Value Index for two management practices was represented by different combinations of species with varied dominance. The herbs’ diversity–dominance curve revealed a lognormal distribution in both managements practices. The overview of distribution patterns for most of the species layer showed contiguous growth and a clumped distribution pattern. Species diversity, richness, herb biomass, basal cover and soil physico-chemical attributes showed a distinct separation in relation to grazing management practices. Based on the findings, one can conclude that the establishment of enclosures has a positive impact in restoring rangeland vegetation diversity, distribution, in increasing herb productivity and in boosting soil fertility.


Introduction
Rangelands provide about 70% of the global forage for both domestic and wild ungulates [1,2].In Africa, rangelands constitute about 65% of the total land area.They are the major sources of feed for ruminants [3].In Ethiopia, rangelands cover about 62% of the country's total landmass.They are largely found in the lowland areas of the country, particularly in areas below 1500 m above sea level.They support about 12-15 million pastoral and agro-pastoral communities in seven regional states of the country [4].The rangeland ecosystems in the arid and semi-arid areas of Ethiopia have suffered from huge land and vegetation degradation due mainly to poor management, population growth, traditional open grazing systems, deforestation, continuous heavy grazing, as well as land use and climate changes [5,6].These problems have led to a huge decline in rangeland resources and degradation of biodiversity, thereby having negative impacts on the rangeland ecosystems, livestock production, and the livelihoods of the pastoral communities [7][8][9].The Somali Regional State is a major pastoral habitat in Ethiopia.It has a total land area of about 327,000 km 2 of which, about 90% is classified as rangeland [8].Like other arid and semi-arid rangeland ecosystems in the country, the rangelands in the region have long been exposed to extensive livestock grazing pressures.These have gradually caused a reduction in vegetation cover and preferred grass species, depletion of soil nutrients and acceleration of erosion [8].Though extensive livestock grazing pressure is dominant in the region, it is common also to see the practice of enclosing communal open grazing rangelands.For instance, the pastoralists in the Harshin District in the region have been using enclosures in response to declining rangeland resources [10].Currently, approximately 80% of the pasture lands in the district are occupied by individuals [11].In the district, the practice of enclosing rangelands started during the early periods of the second half of the twentieth century when settlers from highland areas started to cultivate in Tuguchal, northwest of the district [12].The intention was to discourage farmers from converting rangelands into farmlands.Some studies reported that rangeland enclosures had begun in the 1970s and showed a dramatic increase in the 1980s and 1990s following intense competition for access to land resources between pastoralists and returnees from Somali-Land [13,14].In relation to this, Napier and Desta [11] reported that enclosures are established to reserve resources for use during dry seasons and to generate income by selling pasture/hay.This way of managing and reserving pastoral resources is common among pastoral and agro-pastoral communities in Africa.For instance, the Gogo and Maasai of Tanzania, the Himba of Namibia, and the Borana of Ethiopia traditionally set aside some of their grazing land during the rainy season and resort to it during the dry season [15].
Sumeet et al. [16] noted that findings from studies that address patterns of plant distribution and dynamics of biological resources can provide a rational basis for planning and making decisions about the environment.Otherwise, restoration of resources in their natural habitat would be very difficult [17].The principal aim of phytosociology is to describe the vegetation, explain or predict its pattern, and classify it in a meaningful manner [18].In connection to this, Warger and Morrel [19] noted that phytosociological analysis is important for understanding the functioning of any community.They strongly argue that it is the basis of the ecological study of any piece of vegetation [20].Various studies (e.g., Brand and Goetz [21], Descheemaeker et al. [22], Yayneshet et al. [23], Bakker et al. [24] Allington and Valone [25], and Haftay et al. [10] were conducted to assess the direction of vegetation change in the context of rangeland enclosures to estimate dry matter yields and to indirectly measure the species composition.However, no study has been conducted using the phytosociological data analysis method for rangeland vegetation layers of the rangelands in Ethiopia in general and the Harishin rangelands in particular.Therefore, this study was designed to analyse the phytosociological characteristics of the Harishin rangelands vegetation and soil parameters.Such an analysis is important to provide detailed information on vegetation social status and to offer insights on rangeland ecosystem function and its restoration.
The objective of the study is to analyze the phytodiversity, distribution, herbs biomass and the soil physico-chemical conditions of the vegetation system in the context of traditional grazing management practices in Harishin rangelands of Ethiopia.

Study Area
The study was conducted in the Harshin rangelands of the Somali Regional State in Eastern Ethiopia (located between 9 • 1 40 N and 9 • 22 20 N and 43 • 32 40 E and 43 • 53 20 E), which is about 950 km east of Addis Ababa (Figure 1).The rainfall of the district is bimodal with a short rainy season in April to May.The main rainy season is from June to August.The mean minimum and maximum temperatures of the area are 20 • C and 35 • C, respectively.The mean annual rainfall is 560 mm.The natural vegetation of the study area is characterized as Acacia-wooded grasslands [26].

Site Selection
A reconnaissance survey was made to identify the rangeland sites.Six enclosures, constituting about 12 ha of total area, and the adjacent communal open grazing areas were systematically selected.Enclosures fenced for about 30 years were selected purposely to reduce variability among ages of enclosures.The opinions of community elders and agricultural office experts in the district were used to evaluate the similarity between the enclosures and the adjacent open areas in terms of their protection history, vegetation management and homogeneity in land use [10].To avoid enclosure edge effects, open grazing area plots that lie at least 50 m from the boundaries of each enclosure were selected.The stocking densities of private enclosures were highly variable and depend, among other things, on enclosure size, level of grass regeneration and household herd size [10].The communal rangelands were continuously grazed throughout the year and subjected to over grazing due to high stocking rates.The regional stocking density was estimated to be about 1.4, 56.6, and 3.5 per km 2 for cattle, small ruminants and camels, respectively [10].The vegetation survey was conducted in September and October 2014 when the flowering vegetation was at a late stage.

Herbaceous and Woody Vegetation Sampling
Two 20 m × 20 m (400 m 2 ) plots were randomly laid out for sampling the private enclosures and the adjacent open access communal grazing areas.Finally, a total of 24,400 m 2 plots (2 management type × 6 sampling sites × 2 plots) were considered.In the 400 m 2 plots, the woody species, including both single stemmed and multi-stemmed ones, were recorded.Then, the number of individuals of each tree and shrub species was counted to estimate the diversity, frequency, density, abundance and Importance Value Index (IVI) of the woody vegetation.The basal cover of the woody species was estimated by measuring two canopy diameters at maximum and minimum diameters perpendicular to each other [27].
Within each of the 400 m 2 plots, five 1 m × 1 m (1 m 2 ) quadrants were nested (four at the corners and one at the center) yielding 120 quadrants (24 plots × 5 quadrants).The number of individuals of each herb species was counted to estimate the diversity, frequency, density, abundance, IVI, biomass and basal area cover.In each of the sample sites, the herbaceous

Site Selection
A reconnaissance survey was made to identify the rangeland sites.Six enclosures, constituting about 12 ha of total area, and the adjacent communal open grazing areas were systematically selected.Enclosures fenced for about 30 years were selected purposely to reduce variability among ages of enclosures.The opinions of community elders and agricultural office experts in the district were used to evaluate the similarity between the enclosures and the adjacent open areas in terms of their protection history, vegetation management and homogeneity in land use [10].To avoid enclosure edge effects, open grazing area plots that lie at least 50 m from the boundaries of each enclosure were selected.The stocking densities of private enclosures were highly variable and depend, among other things, on enclosure size, level of grass regeneration and household herd size [10].The communal rangelands were continuously grazed throughout the year and subjected to over grazing due to high stocking rates.The regional stocking density was estimated to be about 1.4, 56.6, and 3.5 per km 2 for cattle, small ruminants and camels, respectively [10].The vegetation survey was conducted in September and October 2014 when the flowering vegetation was at a late stage.

Herbaceous and Woody Vegetation Sampling
Two 20 m × 20 m (400 m 2 ) plots were randomly laid out for sampling the private enclosures and the adjacent open access communal grazing areas.Finally, a total of 24,400 m 2 plots (2 management type × 6 sampling sites × 2 plots) were considered.In the 400 m 2 plots, the woody species, including both single stemmed and multi-stemmed ones, were recorded.Then, the number of individuals of each tree and shrub species was counted to estimate the diversity, frequency, density, abundance and Importance Value Index (IVI) of the woody vegetation.The basal cover of the woody species was estimated by measuring two canopy diameters at maximum and minimum diameters perpendicular to each other [27].
Within each of the 400 m 2 plots, five 1 m × 1 m (1 m 2 ) quadrants were nested (four at the corners and one at the center) yielding 120 quadrants (24 plots × 5 quadrants).The number of individuals of each herb species was counted to estimate the diversity, frequency, density, abundance, IVI, biomass and basal area cover.In each of the sample sites, the herbaceous vegetation was harvested at the ground level using hand shears from the five 1 m 2 nested quadrants to assess the dry matter biomass.Then, the cut samples were weighed using a simple balance immediately.About 20%-30% of this biomass was retained for oven drying (105 • C for 24 h) and used to estimate dry matter.Herbaceous basal cover was estimated visually based on the proportion of area (soil part) covered by a grass base compared to bare ground in the 1-m 2 quadrants.Identification of all species was done in the field.When it was difficult to identify species in the field, the species were recorded and herbarium specimens collected, pressed and dried properly using plant presses.That dried product was soon transported to Haramaya University for proper identification.The nomenclature of the plant species was designated based on the Flora of Ethiopia [28].

Soil Sampling and Analysis
From each sample site, 10 soil samples were taken from the 1 m 2 quadrant located within 400-m 2 plot yielding 120 soil samples (2 plots × 5 quadrate × 6 sample sites × 2 management practices).The soil samples were dug out with an Auger at a depth of 0-20 cm.The soil samples at each site were pooled to form one composite soil sample per sampling site yielding a total of 12 soil samples (2 management practices × 6 sample sites).Finally, the soil PH, electrical conductivity (EC), Total Nitrogen (N), Soil Organic Mater (SOM), Soil Organic Carbon (SOC), the availability of phosphorous (P), exchangeable potassium (K+), texture and bulk density were analysed in the soil laboratory at Haramaya University.
The soil PH was determined using PH meter in 1:2.5 soil water suspension ratios [29].Electrical Conductivity (EC) was determined in a 1:2.5 soil water suspension following the steps and procedures suggested by Chopra and Kanwar [30].Total percentage of nitrogen was determined using the Kjeldahl procedures suggested by Jackson [31].Organic carbon was determined following the method recommended by Walkey and Black.Then, the organic matter of the sample soil was calculated by multiplying the percentage of the organic carbon by a factor of 1.724 following the method recommended by Brady [32].The level of the available phosphorous was determined using the methods and procedures specified by Olsen et al. [33].Available potassium was determined using Ammonium Acetate method at pH 7 [29].Texture was analyzed with the help of the hydrometer method [34].Finally, soil's bulk density was determined by weighing the volume of the sample soil after drying it at 105 • C for 48 h [35].

Data Analysis
The distribution of species similarity between open grazing and enclosed areas was computed using the Index of Similarity (IS).Community Coefficient or Index of Similarity (IS) indicates the percentage of species similarity across the two rangeland management practices.Thus, Index of Similarity (IS) was calculated as follows: where: c is the number of species common to both the forest stands; a, the number of species unique to the enclosure area; b is the number of species unique to open grazing area.Species phytosociological characters were evaluated by analysing the frequency, density, abundance and Importance Value Index (IVI) using a formula recommended in Mandal and Josh [20] Woody species basal cover refers to the proportion of the ground area covered by the vertical projection of the species canopy cover and was used to determine woody species' relative dominance.It was calculated using the following formula [27]: where: CDMax = canopy diameter maximum; and CDmin = canopy diameter minimum Relative Dominance of Woody Species = basal cover of a species basa cove of all the species × 100 IVI of woody species = Relative frequency + Relative density + Relative basal cover (11) The distribution of dominance of plant species in each of the open grazing and enclosed areas was determined using IVI.The herbs IVI data were graphically plotted in MS-Excel.Herbs dominance-diversity curve was drawn using natural log-transformed by plotting the co-ordinate points of its relative importance index (IVI) on the y-axis and its position in the sequence of species from the highest to the lowest IVI on the X-axis [36].A linear fit was applied to rank VI value to delineate the resource sharing patterns by the species in the community and characterized as highest, medium and least important herbs.Moreover, the ratio of abundance to frequency (A/F) for different species was determined to elicit the distribution pattern.The ratio of abundance to frequency indicates regular random (<0.050), contagious (0.050-1.00) and clump (>1.00) distribution patterns.The diversity of the species was computed using the Shannon-Weiner index (H) calculated as [37]: where s = number of species; pi = proportion of individuals or abundance of the i-th species; and Ln is the natural logarithms to the base e.Shannon-Weiner Index (H') was converted to effective number of species diversity [38] using the formula: N1 = Effective number of species; H' = Shannon-Weiner function.
Principal component analysis (PCA) ordination was used to explain variance of analysed soil attributes across the two rangeland management practices.PCA axis correlation coefficient was used to explain the location of the vegetation and soil attributes across the open grazing and enclosed areas.Moreover, species rarefaction curves were drawn by plotting the average number of species as a function of the number of plots to see differences in the level of richness of the species in the open grazing area and enclosure areas.PCA and species richness rarefaction curve analysis were conducted using PAleontological Statistics (PAST) software package version 2.17 [39].Vegetation and soil variables measured for the two management practices were statistically compared using a two-tailed t-test with each management practice being a categorical predictor and vegetation and soil variables as dependent variables.Vegetation and soil variables' homogeneity of variance were tested using Levene's equality of variance test.The computation was done after observing non-significant differences for all measured vegetation and soil variables.The statistical-test analyses were conducted using Statistical Package for Social Sciences version 20 [40].Differences were considered significant at p < 0.05.

Species Similarity between Management Practices
A total of 58 herbaceous species and 11 woody species were recorded in the open grazing areas and enclosure areas.The species distribution of Similarity Index between the enclosures and the open grazing areas were about 47.9% and 45.5% for herbaceous and woody species, respectively.The identified herbaceous species were greater in the enclosure areas than in the open grazing areas.About 52% of the herbaceous species were unique to enclosure areas.However, the identified woody species were greater in the open grazing areas than in the enclosure areas.About 36% of the woody species were unique to open grazing areas.The species of Eragorostis aspera was found to have the highest density (39.43 m −2 ) followed by Eragrostis cilianensis (14.49m −2 ) and Cynodon dactylon (5.70 m −2 ).Regarding the rank of the IVI distribution plot for herb layer sharing in the open grazing management practices, three species were recorded with the highest dominance followed by 23 species of medium dominance and 12 species of least importance (Table 1).A/F = abundance to frequency ratio; IVI = Importance Value Index.

Diversity and Dominance (IVI) of Herb Species in Enclosure Areas
A total of 55 herbaceous species with a 6.55 diversity value were recorded in the enclosure areas.Analysis of the distribution of the herb species in the enclosure areas showed that the dominant species were Cynodon dactylon and Eragrostis cilianensi with a 33.96IVI value.The co-dominating species were found to be Cenchurus ciliaris (IVI = 28.78),Chrysopogon aucheri (IVI = 15.81),Aristida adscencionis (IVI = 13.84), and Dactyloctenium scindicum (IVI = 10.89).Cynodon dactylon was recorded to have the highest density (24.24 individual m −2 ) followed by Cenchures ciliaris (19.68 individual m −2 ) and Eragrostis cilianensi (18.18 individual m −2 ).An 83.33% frequency was observed in Cenchures ciliaris followed species Cynodon dactylon and Chrysopogon aucheri (66.67%) and Ocimum lamifolium (58.33%).The dominance distribution patterns of herb species in the enclosure areas were ranked as 5 herb layers of high dominance, 38 species of medium dominance and 12 herbs of least dominance (Table 2).were recorded for Securinega virosa.Regarding the distribution patterns of woody species in the open grazing areas, two species were found to have high distribution patterns followed by five species with medium distribution and two species with the least distribution (Table 3).A/F = abundance to frequency ratio; IVI = Importance Value Index.

Diversity and Dominance (IVI) of Woody Species in Enclosure Areas
A total of 7 woody species with a 2.39 diversity value were recorded in the enclosure areas.The maximum IVI distribution analysis of the woody species in the enclosure areas showed that the dominant species was Acacia nilotica with an IVI value of 94.45 and a 100% frequency.The co-dominant species were Acacia etbica (IVI = 88.07) and Acacia busie (IVI = 62.39).Woody species of Acacia bussei and Grewia ferruginea were recorded to have the highest and least distribution patterns, respectively (Table 4).A/F = abundance to frequency ratio; IVI = Importance Value Index.

Vegetation Attributes Ordination and Analysis
The PCA correlation variance of herbs species diversity, richness, biomass, basal cover, and woody species basal cover, diversity and richness showed distinct separation (Figure 2) with principal component 1 accounting for 45.95% (eigenvalue = 3.21) and principal component 2 accounting for 20.89% (eigenvalue = 1.46) of the total explained variance in relation to the open grazing and enclosures areas.
Diversity of herbaceous species was found to be significantly higher (t = 2.2; p = 0.038) in the enclosure areas than in the open grazing areas.However, diversity of woody species was significantly higher (t = −2.42;p = 0.047) in the open grazing areas than in the enclosure areas (Table 1).The biomass of the herbaceous species above the ground was significantly different (t = 7.83; p = 0.000) between the enclosure areas and open grazing areas.The basal cover of herbaceous species was found to be significantly different (t = 17.00; p = 0.000) for the two management practices.Moreover, woody basal cover was also found to be significantly different (t = 2.49; p = 0.020).Both the biomass and basal area cover of the herbs in the enclosure areas was found to be more than twice that in the adjacent open grazing area (Table 5).Diversity of herbaceous species was found to be significantly higher (t = 2.2; P = 0.038) in the enclosure areas than in the open grazing areas.However, diversity of woody species was significantly higher (t = −2.42;P = 0.047) in the open grazing areas than in the enclosure areas (Table 1).The biomass of the herbaceous species above the ground was significantly different (t = 7.83; P = 0.000) between the enclosure areas and open grazing areas.The basal cover of herbaceous species was found to be significantly different (t = 17.00;P = 0.000) for the two management practices.Moreover, woody basal cover was also found to be significantly different (t = 2.49; P = 0.020).Both the biomass and basal area cover of the herbs in the enclosure areas was found to be more than twice that in the adjacent open grazing area (Table 5).The species rarefaction curve of the species richness was found to be higher in the enclosure areas than in the open grazing areas.However, the trends of woody species richness were found to be greater in the open grazing than in the enclosure areas (Figure 3).The species rarefaction curve of the species richness was found to be higher in the enclosure areas than in the open grazing areas.However, the trends of woody species richness were found to be greater in the open grazing than in the enclosure areas (Figure 3).

Soil Physico-Chemical Attribute Ordination and Analysis
Principal component analysis of soil physico-chemical attributes revealed distinct separation between the open grazing and enclosure areas with principal component 1 accounting for 59.77% (eigenvalue = 18.63), principal component 2 accounting for 24.56% (eigenvalue = 7.65) and principal component 3 accounting for 11.39 (eigenvalue = 3.55).Moreover, PCA axis correlation coefficient of physico-chemical attributes across the management practices showed distinct separation (Figure 2) with principal component 1 accounting for 40.17% (eigenvalue = 4.8) and principal component 2 accounting for 16.77% (eigenvalue = 2.01).
Soil parameters of N, P, EC and pH were ordinated in the open grazing areas while the soil samples of organic matter (SOM) and organic carbon (SOC) were ordinated in the enclosure areas.
Moreover, from the measured soil attributes only soil samples K (t = 3.49; P = 0.006), pH (t = −3.57;P = 0.005), SOM (t = 5.31; P = 0.000) and SOC (t = 2.26; P = 0.047) were found to be significantly different between the open grazing and enclosure areas whereas the soil parameters N, P, EC, bulk density and soil texture (clay, sand, silt) were found to be insignificantly different across managements practices (Table 6).

Soil Physico-Chemical Attribute Ordination and Analysis
Principal component analysis of soil physico-chemical attributes revealed distinct separation between the open grazing and enclosure areas with principal component 1 accounting for 59.77% (eigenvalue = 18.63), principal component 2 accounting for 24.56% (eigenvalue = 7.65) and principal component 3 accounting for 11.39 (eigenvalue = 3.55).Moreover, PCA axis correlation coefficient of physico-chemical attributes across the management practices showed distinct separation (Figure 2) with principal component 1 accounting for 40.17% (eigenvalue = 4.8) and principal component 2 accounting for 16.77% (eigenvalue = 2.01).
Soil parameters of N, P, EC and pH were ordinated in the open grazing areas while the soil samples of organic matter (SOM) and organic carbon (SOC) were ordinated in the enclosure areas.
Moreover, from the measured soil attributes only soil samples K (t = 3.49; p = 0.006), pH (t = −3.57;p = 0.005), SOM (t = 5.31; p = 0.000) and SOC (t = 2.26; p = 0.047) were found to be significantly different between the open grazing and enclosure areas whereas the soil parameters N, P, EC, bulk density and soil texture (clay, sand, silt) were found to be insignificantly different across managements practices (Table 6).

Species Phytosociological Characteristics
Analysis of IVI provides information about the social status of a species and can be recognized as patterns of association of dominant species in a community [41].Analysis of IVI in the two management practices represented different combinations of species with different dominants and co-dominants.Eragorostis aspera, Eragrostis cilianensis, Cynodon dactylon, Tragus racemosus, and Cenchurus ciliaris were found to be the most dominant herb species in the open grazing areas while Cynodon dactylon, Eragrostis cilianensis, and Cenchurus ciliaris were found to be the highest IVI and density in the enclosure areas.Supporting the findings of Mandal and Joshi [20] on open grazing areas, the findings in the current study suggest that vegetation experiencing stresses from biotic pressure are under serious threat.However, some of the herb species in the open grazing areas have managed to survive.This could be attributable to their broad ecological amplitude and greater adaptability against biotic influences.In relation to this, McGranahan et al. [42] noted that grazing might be associated with variation in community composition within vegetation states.Vegetation structure and composition may be influenced by novel factors such as fire, herbivore, topography and past management practices that interact in a complex way to determine plant community architecture [43][44][45].
Variation of herb distribution in the enclosure areas may be attributed to the survival and reproduction of maximum species due to moderate level of species competition during early regeneration which has led to the domination of only a few species.This condition has affected the plant species which are susceptible to species' inter-competition.Validating this finding, Grime [46] noted that with an increase in environmental stress, the species adapted to low levels of environmental stress lose their competitive advantage whereas those that are more resistant to environmental stress can increase in abundance.It is uncertain, however, the extent to which the smallest differences between samples of vegetation are environmentally determined and the extent to which they affect the chances of species establishing themselves first [47].
The herbaceous species' dominant diversity curves followed a lognormal distribution in both management practices (Figure 4).The lognormal dominance-diversity curves indicate the heterogeneity of the species [48].Lognormal hypothesis assumes that the importance of species is governed by the interactions between a large numbers of factors determining success in the niche hyperspace [36].In connection to this, Whittaker [49] noted that the log-normal series describes the partitioning of realized niche space among various species and is the consequence of the evolution of a particular species' diversity along the niche parameters which they exploit.Similar types of herb dominance-diversity curves have been reported in the forest ecosystem [50,51].

Species Phytosociological Characteristics
Analysis of IVI provides information about the social status of a species and can be recognized as patterns of association of dominant species in a community [41].Analysis of IVI in the two management practices represented different combinations of species with different dominants and co-dominants.Eragorostis aspera, Eragrostis cilianensis, Cynodon dactylon, Tragus racemosus, and Cenchurus ciliaris were found to be the most dominant herb species in the open grazing areas while Cynodon dactylon, Eragrostis cilianensis, and Cenchurus ciliaris were found to be the highest IVI and density in the enclosure areas.Supporting the findings of Mandal and Joshi [20] on open grazing areas, the findings in the current study suggest that vegetation experiencing stresses from biotic pressure are under serious threat.However, some of the herb species in the open grazing areas have managed to survive.This could be attributable to their broad ecological amplitude and greater adaptability against biotic influences.In relation to this, McGranahan et al. [42] noted that grazing might be associated with variation in community composition within vegetation states.Vegetation structure and composition may be influenced by novel factors such as fire, herbivore, topography and past management practices that interact in a complex way to determine plant community architecture [43][44][45].
Variation of herb distribution in the enclosure areas may be attributed to the survival and reproduction of maximum species due to moderate level of species competition during early regeneration which has led to the domination of only a few species.This condition has affected the plant species which are susceptible to species' inter-competition.Validating this finding, Grime [46] noted that with an increase in environmental stress, the species adapted to low levels of environmental stress lose their competitive advantage whereas those that are more resistant to environmental stress can increase in abundance.It is uncertain, however, the extent to which the smallest differences between samples of vegetation are environmentally determined and the extent to which they affect the chances of species establishing themselves first [47].
The herbaceous species' dominant diversity curves followed a lognormal distribution in both management practices (Figure 4).The lognormal dominance-diversity curves indicate the heterogeneity of the species [48].Lognormal hypothesis assumes that the importance of species is governed by the interactions between a large numbers of factors determining success in the niche hyperspace [36].In connection to this, Whittaker [49] noted that the log-normal series describes the partitioning of realized niche space among various species and is the consequence of the evolution of a particular species' diversity along the niche parameters which they exploit.Similar types of herb dominance-diversity curves have been reported in the forest ecosystem [50,51].Dominance of a few species in both management areas indicated patterns of species distribution.The distribution of patterns (A/F ratio) for most of the species' layers showed contiguous growth patterns followed by a clumped distribution pattern.According to Odum [18], contiguous distribution is the most common pattern in nature and is formed as a result of small but significant variations in the ambient environmental conditions.The current study also suggests that clumped distribution often occurs due to an uneven distribution of nutrients or other resources in the environment.Variations in the distribution pattern among sites and vegetation composition are associated with micro environmental and biotic factors [52].Patterns of distribution depend both on the physicochemical nature of the environment and the biological peculiarities of the organisms [18].

Interaction of Management and Vegetation Attributes
The PCA correlation coefficient of the measured herbaceous vegetation across the management practices showed distinct variation with a greater value in the enclosures areas than in the open grazing areas.In the study area, continuous open grazing has caused a reduction in species richness, plant abundance, standing biomass and basal cover.This is attributed to the soil parameters and pressure of high and unregulated grazing on freely accessed pastures.The total variation explained by the soil parameters in multivariate tests was about 60%.This shows that there are other factors that influence the structure and composition of vegetation in the area.In line with this finding, Yayneshet et al. [23] reported that grazing throughout the year has consistently reduced the above ground biomass production capability of the grazed area.Similarly, Bilotta et al. [53] reported that heavy and uninterrupted grazing can cause a reduction in herbaceous vegetation and a sharp decline in the concentration of soil nutrients, thereby leading to defoliation and soil compaction.Since they consume aboveground biomass, herbivores reduce the amount of litter that falls on the soil surface.Excessive plant defoliation from over grazing damages plant tissues that, as a result, reduces the diversity of plant species and the percentage cover of the herbaceous vegetation [5].
The findings in the current study are in line with similar studies on the pastoral systems of Ethiopia which reported the existence of higher herbage mass and basal cover in the enclosure area [7,9].The high diversity, richness and biomass measured in the enclosures might be attributable to increased litter accumulation, improved soil organic matter and other nutrients inside the enclosures [23,54].According to Tessema et al. [5], higher soil nutrient contents under light grazing conditions are attributed mainly to larger basal cover and greater amounts of standing biomass.Increase in herbage mass in enclosure areas could be linked partly to a seasonal reduction in grazing pressure and partly to the subsequent accumulation of soil organic matter during the wet season resting period [54].In the study area, a reduction in woody species richness inside the enclosure areas is attributed mainly to individual pastoralists' management trends in enclosure areas.Most of the pastoralists in the study area do not prefer woody species in their enclosure areas.They frequently clear most of the shrubs and trees inside the enclosures by leaving a few species, which are used to shade their animals from direct hot sun.A low density of woody species has been shown to increase soil nutrient status and thus encourage grass-growing conditions [55,56].The overall density threshold of woody species in terms of bush encroachment is estimated to be more than 2500 trees equivalent ha -1 [57].However, measurement of the density and basal cover of woody species revealed that they were not large enough to be recognized as being in a woody-encroached state that substantially suppresses the herbaceous biomass.This indicates that the pastoralists in the study area are conscious of the potential threat of woody species and, because of this, they frequently clear most of the shrubs and trees that are not preferred by their livestock.This finding supports previous findings [10] that differences in the herbaceous biomass between enclosures and open-access areas may occur may be attributed more to grazing impacts than to the removal of woody species [10].

Soil Physico-Chemical Attributes
Soil nutrients K, SOM and SOC were significantly higher in the enclosure areas than in the open grazing areas.The higher soil nutrient contents inside enclosure areas are attributed mainly to the higher vegetation coverage and the role of increased biomass in nutrient cycling in enclosure areas, which result in higher litter input and improved content of organic matter in the soil.In line with this study, Reeder and Schuman [58] reported the existence of higher soil organic carbon in lightly grazed sites compared to heavily grazed ones.This happens due to an increase in above ground plant litter and annual forbs and grasses that enhance the formation and accumulation of SOM [59].Among the soil properties, the total organic C is a sensitive soil quality indicator [60].This means that within a narrow range of soil the organic C may serve as a suitable indicator of soil fertility [61].The SOC fraction may offer further insight into soil fertility changes and the sustainability of land-use and management practices [62].
The soil samples from the open access communal grazing area showed a significantly higher pH value than those taken from the enclosure areas.Moreover, high concentrations of soil parameters like N, P, PH, EC and soil silt content were observed in the open grazing areas.Increase in the concentrations of total N and P under open grazing is probably due to the existence of high animal grazing that contributes to increasing animal dung, urine and defoliation of plants during grazing.Validating this finding, Yates et al. [63] reported that concentrations of soil nitrogen and the availability of phosphorous were significantly higher in heavily grazed rangelands.Complex spatial patterns of soil nutrients have been commonly presumed to develop over time as a result of the interactions of climate, parental material, vegetation type and topography [64].In connection to this, Tefera et al. [60] reported that the rangelands of east Africa generally have a low level of soil fertility.

Conclusions
The findings in the current study generally show that grazing management practices have an effect on the dominance distribution, diversity, herb biomass and soil attributes of vegetation species in the Harishin rangelands of Somali Regional State, Ethiopia.The stress of unregulated grazing on open grazing areas is generally linked to the diminished diversity and abundance of herb species due to the dominance of a few species having greater adaptability against grazing influences.This is attributed to a reduction in standing biomass, basal cover and soil nutrients such as soil organic matter.It can be concluded that the establishment of enclosures has a positive effect in restoring the diversity, productivity and soil fertilities of rangeland vegetation.Since the findings of the current study area were based on small samples, the researcher calls for a more broad-based investigation to provide a stronger basis for characterizing the functioning of rangeland systems of the pastoral ecosystem.
Supplementary Materials: Supplementary data related to this article is available online at: www.mdpi.com/2073-445X/6/1/4/s1. Table S1.Attributes for herbaceous species in enclosures areas averaged for 12 plots used for analysing frequency, density, abundance, and Importance Value Index (IVI) in Harishin Rangelands of Ethiopia; Table S2.Attributes for herbaceous species in open grazing areas averaged for 12 plots used for analysing frequency, density, abundance, and Importance Value Index (IVI) in Harishin Rangelands of Ethiopia; Table S3.Attributes for woody species in enclosure areas in 12 plots used for analysing frequency, density, abundance, and Importance Value Index (IVI) in Harishin Rangelands of Ethiopia; Table S4

3. 2 . 1 .
Diversity and Dominance (IVI) of Herb Species in Open Grazing Area A total of 38 herbaceous species with a 4.51 diversity value were recorded in the open grazing areas.Distribution analysis of the herb species in the open grazing areas showed that the dominant species was Eragorostis aspera with an IVI value of 75.29.The co-dominating species were Eragrostis cilianensis IVI =36.11),Cynodon dactylon (IVI = 15.81),Tragus racemosus (IVI = 12.30), and Cenchurus ciliaris (IVI =10.99).A 58.33% frequency was observed in three species, namely Eragorostis aspera, Cynodon dactylon and Cenchurus ciliaris.

3. 2 . 3 .
Diversity and Dominance (IVI) of Woody Species in the Open Grazing Areas A total of 9 woody species with a 3.17 diversity value were recorded in the open grazing areas.The maximum IVI (58.13) value and the highest density (8.10 individual m −2 ) were recorded for Acacia busie.The co-dominating species were Acacia nilotica (IVI = 51.23),Acacia etbica (IVI = 51.23) and Acacia toritilis (IVI = 39.55).The minimum IVI (9.13) value and the lowest density (0.75 individual m −2 )

Figure 3 .
Figure 3. Rarefaction curves of cumulative increase of vegetation species richness (Mean ± SD) at 95% confidence for both of the enclosure and open grazing management practices in Harishin rangelands, Ethiopia: (a) Herbaceous species richness; (b) Woody species richness.

Figure 3 .
Figure 3. Rarefaction curves of cumulative increase of vegetation species richness (Mean ± SD) at 95% confidence for both of the enclosure and open grazing management practices in Harishin rangelands, Ethiopia: (a) Herbaceous species richness; (b) Woody species richness.

Figure 4 .
Figure 4. Herbaceous vegetation dominance-diversity curve based on IVI in Harishin rangelands of Ethiopia: (a) Communal open grazing areas; (b) Traditional private enclosure areas.

Figure 4 .
Figure 4. Herbaceous vegetation dominance-diversity curve based on IVI in Harishin rangelands of Ethiopia: (a) Communal open grazing areas; (b) Traditional private enclosure areas. .

Table 1 .
Herbaceous species distribution in open grazing areas of Harishin rangelands, eastern Ethiopia.

Table 2 .
Herbaceous species distribution in enclosure areas of Harishin rangelands, eastern Ethiopia.

Table 3 .
Woody species distribution in open grazing areas of Harishin rangelands eastern Ethiopia.

Table 4 .
Woody species distribution in enclosure areas of Harishin rangelands, eastern Ethiopia.

Table 5 .
Vegetation attributes (mean ± standard error) of enclosures and open-access communal grazing areas in Harishin rangelands, Ethiopia.

Table 5 .
Vegetation attributes (mean ± standard error) of enclosures and open-access communal grazing areas in Harishin rangelands, Ethiopia.

Table 6 .
Soil physico-chemical attributes (mean ± standard error) of enclosures and open access communal grazing areas in Harishin rangelands, Ethiopia.

Table 6 .
Soil physico-chemical attributes (mean ± standard error) of enclosures and open access communal grazing areas in Harishin rangelands, Ethiopia.

.
Attributes for woody species in open grazing areas in 12 plots used for analysing frequency, density, abundance, and Importance Value Index (IVI) in Harishin Rangelands of Ethiopia; Table S5.Attributes for Vegetation parameters used principal component analysis ordination and statistically compared from six plots in enclosures and open-access communal grazing areas of Harishin rangelands, Ethiopia; Table S6.Attributes for Soil physico-chemical parameters used principal component analysis ordination and statistically compared from twelve plots in enclosures and open-access communal grazing areas of Harishin rangelands, Ethiopia; Table S7.Herbaceous species relative dominance, relative frequency and relative density and IVI distribution analysis in open grazing areas of Harishin rangelands, Ethiopia;

Table S8 .
Herbaceous species relative dominance, relative frequency and relative density and IVI distribution analysis in enclosure areas of Harishin rangelands, Ethiopia; TableS9.Woody species relative dominance, relative frequency and relative density and IVI distribution analysis species in enclosure areas of Harishin rangelands, Ethiopia; TableS10.Woody species relative dominance, relative frequency and relative density and IVI distribution analysis species in open grazing areas of Harishin rangelands, Ethiopia.