Abundance and diversity of soil invertebrate macro-fauna in different land uses at Shenkolla watershed, South Central Ethiopia

Land use can exert a strong influence on the abundance, diversity, and community composition of soil macro-fauna. This study was conducted to evaluate the effects of four land use types on the abundance and diversity of soil invertebrate macro-fauna communities. These land uses include forest land, grazing land, crop cultivated outfields, and homestead garden fields present at Shenkolla watershed, south central Ethiopia. Monolith sampling of soil macro-fauna was done according to the standard of Tropical Soils Biology and Fertility Institute (TSBF) procedure. Five sampling points were chosen in each land use type and small monolith (25 × 25 × 30 cm) was dug out at 5-m interval along a transect with randomly positioned starting point, but perpendiculars to the slope. A total of 20 monoliths (4 treatments × 5 replications) were taken across all the land use types. Sampling of SIMF was carried out in April 2019 where soil macro-faunas are known to be more active. To evaluate the SIMF community eight parameters were measured: Shannon-Wiener index, Simpson diversity index, Pielou’s measure of evenness, Margalef’s diversity index, the Number of Occurrence Index, Relative abundance, Density (individuals per square meter) of each taxon and density of all SIMF and Bray-Curtis similarity index. The data were further analyzed using ANOVA and a general linear model to determine the variation and the influence of land use type, respectively. In general, 332 individuals, 10 orders, 12 families, and 15 species were identified, from the collected samples. There were significant differences (p < 0.05) among the four land use types for SIMF except wireworm, spiders, and millipedes. Overall abundance and diversity were lowest in the crop cultivated outfields and highest in homestead garden fields and forest land. Bray-Curtis’ similarity was highest between the sampled sites s16 and s19 with in the forest land, and lowest (2%) between sites s1 (crop cultivated out fields) and s14 (homestead garden fields). The results revealed that the diversity of SIMF was positively influenced by forest land and homestead garden fields and negatively influenced in grazing land and cultivated outfields. Therefore, maintenance of a continuous litter cover at the surface and application of a wide range of organic fertilizers (farmyard manure, household refuse, and compost) is very much critical to prevent the decrease in diversity of SIMF.


Background
Land use can exert a strong influence on the abundance, diversity, and community composition of soil macrofauna (Barros, Pashanasi, Constantino, & Lavelle, 2002). As land use changes, so do species compositions and distributions. Land use changes and intensification are severe enough to influence plant communities, and it is likely that they trickle down to soil food webs and change linkages between above and below-ground communities (Wardle et al., 2004). In addition, annual cropping systems decrease the diversity and abundance of soil fauna communities due to soil disturbance and the absence of a permanent soil cover required for the staying alive of the soil organisms (Barrios, 2007;Rossi et al., 2010). When macro-fauna communities of important functional groups are affected by land use and land cover change, biodegradation of organic residues and essential ecosystem functions may also suffer (Moreira, Huising, & Bignell, 2008). Soil macro-fauna are also thought to be influenced by unsustainable land management such as overgrazing, fire, deforestation, pollution, soil erosion, and depletion of fertility (Bignell et al., 2005).
The main features of forest land and homestead garden field are that minimal soil disturbance and better soil cover. These features create more favorable conditions for the development of soil organisms. While these features are practically absent in grazing land and cropcultivated out fields, resulting in higher soil degradation, increased compaction, and the lack of food, especially the absence of cover necessary for the survival of the soil organisms (Moreira et al., 2008;Rossi et al., 2010).
The SIMF distribution, abundance, and diversities are known to vary depending on residue inputs and soil management practices (Brown et al., 2004;Manhães, Gama-Rodrigues, Moço, & Gama-Rodrigues, 2013;Mutema, Mafongoya, Nyagumbo, & Chikukura, 2013). It can therefore be hypothesized that different land management practices have an effect on SIMF differentially. Due to the lack of previous research in Shenkolla watershed, there is information gap on local SIMF assemblages in different land use types. So investigation of the potential of SIMF as a bioindicator in different land use systems under this specific watershed is required to come up with conclusion that the land use change affects the abundance and diversity of SIMF. Thus, in view of this, the present study was conducted to evaluate the pattern of abundance and diversity of SIMF in four land use types at Shenkolla watershed, south central Ethiopia.

Description of the study area
The study was conducted within four main land use types: forest land, homestead garden field, grazing land, and cultivated outfield at Shenkolla watershed which is located in the eastern part of Soro district. Geographically, the area falls within the coordinates of 7°24′ 30″-7°27′ 0″ N latitude and 37°43′ 30″-37°46′ 30″ E longitude (Fig. 1). The altitude ranges from 2200 to 2830 m.a.s.l.
The climate of the watershed is characterized generally as tepid sub-moist mid highland with long-term average rainfall of about 1107 mm with bi-modal pattern having light rainy season is usually from March to May and the heavy rainy season is from June to September. The annual average temperature is 17.2°C (Fig. 2).
The dominant soil type is Nitisol in all land use types. According to the rating proposed by (Landon, 1991), the soil pH of the watershed is moderately acidic in the cultivated outfields (5.80) but neutral under forest (7.20), homestead garden fields (6.83), and grazing lands (6.57) (Bufebo & Elias, 2020). This suggesting the intensive land use including the application of mineral fertilizers in cultivated outfields led to acidification.
Mixed crop-livestock system is the major source of livelihood for the community in the study watershed. The system is noted for its high population densities (200-350 persons per kilometer square). Also, there is a severe land shortage (average holdings of 0.5 ha for a family of 8 persons) along with intensive cultivation. Forest and grazing lands are communally owned and managed while the arable lands are individually owned (Bufebo & Elias, 2020).
Livestock husbandry is based on free grazing on communal grazing lands. Free grazing' is an age-old traditional system which allows owners to indiscriminately graze their livestock on communal land. Especially during the cropping season, all livestock are confined to the scarce grazing lands which for 5-6 months of the year are subject to immense grazing pressure. Large herd size on small grazing lands and poor pasture management increased the pressure on grazing land of the study area (Elias, 2016). As a result, the animals cannot get enough fodder to stay healthy and in good condition; similarly, the natural vegetation has no chance to recover at any time of year. There is no reseeding effect, the most palatable grasses and legumes have disappeared, and bare patches have developed, giving room for accelerated soil erosion and severe dissection by rills and gullies (Bufebo & Elias, 2020).
Soil fertility management is clearly differentiated between the cultivated outfields and homestead garden fields. Hoeing and incorporation of farmyard manure in homestead garden fields are distinctly different from the plough-based complex system in the intensively cultivated outfields (Elias, 2016). The homestead garden fields receive the application of a wide range of organic fertilizers (farmyard manure, household refuse, compost, and leaf litter). Soil fertility in the homestead garden fields is maintained through the application of approximately nine tons per hectare annually of farmyard manure on average, while the crop cultivated outfields are treated with a dose of less than prescribed amounts of mineral fertilizer with an average rate of 50 kg urea and 65 kg/ha DAP (Di-ammonium phosphate: 18% N, 46% P2O5) (Elias, Okothc, & Smalingd, 2019). Crop residue removal is another problem that causes soil fertility decline in cultivated outfields. As the   result, cultivated outfields are largely depleted of soil fertility but homestead garden fields are enriched. The differences in land use and management practices indicated that there was a difference in extent of water erosion in the study area. Field observation indicated that the presence of slight water erosion in forest land and homestead garden fields and accelerated water erosion in grazing land and cultivated outfields at the study site. This shows that the susceptibility of the soils of cultivated outfields and grazing land to water erosion (Bufebo & Elias, 2020).

Sampling methods and experimental design
Sampling was conducted in April 2019 at the beginning of slight rainfall locally called belg when soil macrofaunas are known to be more active (Seeber et al., 2005, Tondoh & Lavelle, 2005. Monolith sampling of soil macro-fauna was done according to the standard of Tropical Soils Biology and Fertility Institute (TSBF) procedure . Five sampling points were chosen at each land use type. At each sampling point, one small monolith (25 × 25 × 30 cm) was dug out at 5m interval along a transect with randomly positioned starting points but perpendiculars to the slope (modified TSBF protocol, Anderson & Ingram, 1993). The total of 20 sampling points was dug out by taking 5 representative monoliths from each land use type. Each sample was taken to the sampling base and hand sorted, removing all the macro-fauna > 2 mm in diameter that were visible to the naked eye (Lavelle, Senapati, & Barros, 2003). Soil macro-faunas were counted and preserved in vials filled with 70% alcohol and labeled. The soil was returned to the site of removal after the extraction of soil organisms to minimize destruction and fragmentation of the vulnerable soils of different land use types.

Identification of soil macro-fauna
The sampled organisms were then taken to the Insect Sciences laboratory of Addis Ababa University for identification. Species of soil invertebrate macro-fauna (SIMF) were identified in laboratory using identification key and pictures (Eaton & Kaufman, 2007;Zettler et al., 2016). Macro-fauna abundance was determined in four land use types according to the following main groups: earthworms (Haplotaxida), termites (Isoptera), ants (Hymenoptera), and other macro-fauna. Finally, soil organisms were separated into 15 taxonomic species; density and abundance of each species in four land use classes were computed.

Statistical analysis
Eight community parameters, including Shannon-Wiener index, Simpson index, Pielou's index of evenness, Margalef's diversity index, Bray-Curtis similarity index, the number of occurrence index, relative abundance, and density, were computed for land use types and used to compare SIMF diversity across the land uses.
Shannon and Weiner index (Shannon & Weaver, 1949) Shannon-Weiner Diversity Index (H') was used as a measure of species abundance and richness to quantify the diversity of the SIMF. The H′ was calculated using the Shannon and Weaver (1949) is given by the equation: where H′ is Shannon-Wiener Diversity Index, Pi is the proportion of each species of SIMF in the sample, and Ln (Pi) is the natural logarithm of this proportion.
Simpson index (Simpson, 1949) Simpson Index measures the probability of any two individuals drawn from noticeably large community belonging to different species. It was measured by the following formula: where n is the total number of SIMF of a particular species and N is the total number of SIMF of all species.
Pielou's measure of evenness (Pielou, 1966) Pielou's evenness index (J′) was calculated using the ratio of observed diversity to maximum diversity using the equation.
where J′ is evenness index, H′ is the Shannon-Wiener Diversity Index, and LnS is the natural log of a total number of observed species.
Margalef's diversity index (Margalef, 1968) Margalef's diversity index is calculated using the formula: where d = Margalef's diversity index, S = number of species, and N = number of individuals.
Bray-Curtis similarity index (Bray & Curts, 1957) Bray and Curts (1957) standardized the Manhattan metric so that it has a range from 0 (similar) to 1 (dissimilar). Bray-Curtis coefficient is computed by dividing the shared abundance by the total abundance. For two plots i and k, B i,k was computed as follows (Bray & Curts, 1957): where B i,k is Bray-Curtis coefficient, sij and skj are the number of species of j in land uses i, and k respectively; p = total number of species.

Number of occurrence index (NOI)
This is the total number of individuals of each species in a sample, expressed as a percentage of the total number of individuals of all species in the sample.
where A = number of individual of each species in the sample and B = total number of individual of all species in the sample.

Relative abundance
The relative abundance of SIMF per land use class was determined using the following formula where n is the total number of SIMF of a particular species and N is the total number of SIMF of all species.

Density of SIMF
The density of SIMF was determined using the following formula where c is the total number of SIMF of a particular species in one square meter of soil and m 2 is the area of soil calculated by multiplying 1 m by 1 m.
One-way analysis of variance (ANOVA) was performed to compare the variations in SIMF abundance among different land uses. A generalized linear models (GLMs) analysis was conducted to determine the influence of land use class (fixed factor), with 4 levels: (i) cultivated out fields, (ii) grazing lands, (iii) homestead garden fields, and (iv) forest lands, on the SIMF abundance (response variable). Statistical package for SPSS v.16.0 (SPSS (Statistical Package for Social Science), 2007) for windows was used to carry out ANOVA and GLMs.

Composition and occurrence of SIMF community
A total of 332 individuals and 15 species, belonging to 12 families and 10 orders were then identified, from the four land use types (Table 1). Ants, beetle adults, earthworms, spiders, and millipedes were found in all land use types. Wireworms (Coleoptera: Elateridae) and earwigs were only sampled in the forest land. Beetle larvae, centipede, and snails were not sampled in the crop cultivated outfield and grazing land use types, while majority were found in the forest lands (15 species) and homestead garden fields (12 species) ( Table 1). This study confirmed that earth worms, termites, beetle adults, and ants were dominant species in the study watershed. The NOI of SIMF in the homestead garden field was 44.58% value 3.4 times larger than that in the cultivated outfields (Table 2).

Density and relative abundance of SIMF
The mean density of earthworms, termites, adult beetles, and ants was relatively high (141.6, 37.60, 19.2, 19.2 individuals/m 2 , respectively) showing their dominance over other species (Table 2). The homestead garden fields had high densities of earthworms (230.4 individuals/m 2 ) ants (48.00 individuals/m 2 ), beetle adult (44.80 individuals/m 2 ), and beetle larva (41.60 individuals/m 2 ) (Fig. 3). The forest lands had a high density of termites (83.20 individuals/m 2 ). However, the lowest densities were observed in the cultivated outfields (Fig. 3). Significant variations among land use types were detected for density of most SIMF groups (P < 0.05), except for wireworm, spiders, and millipedes (P = 0.83, 0.70, 0.74, respectively) ( Table 3). The homestead garden fields had a significantly higher mean density of SIMF than the forest, cultivated outfields and grazing lands ( Table 3). The mean density of earthworms and ants were significantly highest in the homestead garden fields (230.4 ± 4.59, 48.00 ± 25.29, respectively). The mean density of termites reached its maximum (83.23 ± 9.51) in the forest land use (Table 3). The earthworms were ranked 1st as this group has the highest relative abundance (0.533), while wireworms were ranked 10th since they had the lowest relative abundance (0.006) ( Table 2).

Diversity and similarity of SIMF across land use types
The SIMF communities within stress conditions, in the crop cultivated outfields and grazing lands, exhibiting lower diversity than stable communities, in the homestead garden fields and forest lands ( Table 2). Values of species richness, Simpson diversity, evenness, and Margalef diversity were highest for forest land use (15, 0.83, 0.428, 2.56, respectively) while they turned to be lowest for cultivated outfields (5, 0.22, 0.075, 1.06, respectively) ( Table 2). Species richness is high in forest land whereas low in cultivated outfields. Shannon diversity of SIMF showed its maximum value in the homestead garden fields (0.279) followed by forest lands (0.262) and the minimum value in the cultivated outfields (0.121) ( Table 2).
From the data processed according to Bray-Curtis similarity index (Table 4), it is evident that the largest dissimilarity of SIMF communities (0.98) is between the sampled sites s1 (crop cultivated outfields) and s14 (homestead garden fields), and the smallest dissimilarity (0.02) is between s16 and s19 within the forest lands. These means that the species common to s1 and s14 is small (only 2% similarity). The two sample sites (s1 and s14) share low number of species in common as compared with s16 and s19 (98% similarity) ( Table 4).

The effect of land use on abundance of SIMF
The most numerous macro-fauna group present across all land use types were earthworms with a total abundance of 177 across all land use types. This constitutes 53.3% of all macro-fauna. Termites were with a total abundance of 47 (14.2%) of all macro-fauna, and all the rest macrofauna constituted an abundance of 108 (32.5%) of total macro-fauna (Table 5). Earthworms were the dominant group of SIMF, followed by termites. High abundance of earthworms (72) was recorded in the homestead garden fields. The abundance of termites reached its maximum (26) in the forest land use ( Table 5). Values of an abundance of SIMF were in the order of homestead garden fields > forest land > grazing land > cultivated out fields (Table 5). Overall, the abundance of SIMF was lowest in crop cultivated outfields and highest in homestead garden fields followed by forest land (Table 5).
The GLM model indicated that there was no significant effect of land use on the SIMF groups, except earthworm (Table 6). Since we had multiple fixed effects, then the significance of the overall model and the

Composition and occurrence of SIMF community
The SIMF community clearly responded to the environmental disturbance induced by unsustainable land use management. These groups of SIMF living in forest and homestead garden fields were favored, possibly via the large production, and/or better quality of litter. This owing to the fact that high organic matter content under forest and homestead garden fields which provided substrates for the soil organisms, likely reduces the negative effects of soil acidity on soil organisms (Ayuke et al., 2011;Beare, Reddy, Tian, & Srivastava, 1997). Also, the crop cultivated outfields had higher soil acidity concentrations and low organic matter content (Bufebo & Elias, 2020) considered to be "not convenient" for many SIMF.

Diversity and similarity of SIMF across land use types
Homestead garden fields and forest lands had more species diversity explaining the fact that unregulated agricultural expansion into forest land had a negative impact on SIMF species diversity. Diversity decreased gradually with an increasing intensification of land use in the system. This might be hypothesized that the low vegetation diversity led to impoverished SIMF species diversity and high vegetation diversity led to rich SIMF diversity.  The effect of land use on density and abundance of SIMF The density of SIMF communities varied significantly in relation to the land use. In our study, an increase in density of SIMF in the forest and homestead garden fields was largely accounted by variations in litter quality and /or abundance, occurrence of vacant niches, good pH, and moisture content at the soil surface (Bufebo & Elias, 2020). This may be attributed to the variation of disturbance level in the habitats within these land use types on the community composition of SIMF (Ayuke et al., 2011;Barrios, 2007;Rossi et al., 2010). Homestead garden fields and forest land had more species abundance explaining that the fact that unregulated agricultural expansion into forest land had a negative impact on the abundance of SIMF species. The abundance of SIMF was dramatically affected by grazing land and crop cultivated outfields. Similarly, (Decaëns, Lavelle, Jimenez, Escobar, & Rippstein, 1994) reported that the abundance of soil macro-fauna tends to decrease to low levels in crop cultivated lands. The higher the vegetation diversity, so the ability to provide energy and food source for soil macro-fauna would be high. The higher the availability of energy and nutrients for soil macro-fauna, so growth and activity of soil macrofauna would be better (Negasa et al., 2017).
As intensification occurs in crop cultivated outfields, losses and stresses imposed by chemical contamination through use of herbicides and pesticides, as well as chemical imbalances through soil acidification resulted in gradual decrease in SIMF diversity (Elias, 2016). Moreover, continuous grazing and poor pasture management in grazing land of the study area led to soil degradation through nutrient depletion and compaction impaired soil biological functioning (Bufebo & Elias, 2020). Similarly, (Muchane et al., 2012) reported that agricultural activities significantly affected the composition of macro-fauna community. In this study, increasing in abundance of earthworm (39) in crop cultivated out fields showed its more adaptability while earthworms have been found to be sensitive to disturbance . Several studies have noted abundant termite and earthworm populations in agricultural land uses, especially those with limited-or no-tillage management (Ayuke et al., 2009;Mutema et al., 2013). The relatively low abundance of SIMF noted within the crop cultivated field may result from the use of certain agricultural chemical inputs suggesting that technological intensification (Elias, 2016;New, 2005). Modifications of macro-fauna communities are known to have potential negative effects on soil functioning and on the sustainability of ecosystems (Decaëns et al., 2004). Although grazing land and cultivated outfields had negative impacts on the diversity and density of SIMF species, still significant numbers of species do exist on grazing land and cultivated outfields calling for the need to consider these land uses for conservation and maintenance of SIMF species.
SIMF are sensitive indicators of the nature of land use and management (Negasa et al., 2017). The homestead garden fields and forest land had higher abundance and diversity than any other type of land use. Biological activity is concentrated in the forest and homestead garden fields. This might be attributed by the availability of potential substrate and quality of the surface litter (Ayuke et al., 2009;Manhães et al., 2013). Additionally, the homestead garden fields receive the application of a wide range of organic fertilizers and partly decomposed leaves and twigs from multipurpose trees (Elias, 2016). This is an indication that these land uses may sustain sufficiently abundant and diverse communities to optimize the effects of these beneficial organisms (Brown et al., 1999). The quality of organic matter produced and the effects of vegetation on soil water and temperature regimes may vary considerably, with significant effects on the SIMF community . Some important factors influencing the abundance, distribution, and diversity of SIMF communities include soil quality, immediate substrates, and food resource availability (Barros et al., 2002). Ecological imbalance arising from any severe alterations of these factors may affect the natural environment in which SIMF communities live (Dash, 2003;Rybalov, 1990).
The results demonstrate that human disturbances are an important factor influencing macro-fauna community. Generally, this study shows that soil macro-faunas are sensitive to land use and management practices. The reason for this might be connected with the anthropogenic activities that result in disturbance of habitat. The finding of this study is in agreement with (Barrios 2007) who reported that land use can exert a strong influence on the abundance, diversity, and community composition of SIMF. On the other hand, the result of this study disagrees with the findings of  (Manetti et al., 2010) who explained that tillage system does not affect soil macro-fauna in the southeastern Buenos Aires Province, Argentina.
Our study seems to indicate that maintenance of a continuous litter cover at the surface, application of a wide range of organic fertilizers (farmyard manure, household refuse, and compost), and minimizing compaction and over all disturbance of habitat may prevent the decrease in diversity of SIMF.

Conclusions
Results obtained from this study demonstrate that quantitative changes in the abundance and diversity of SIMF communities occur when land is intensively cultivated. The resultant changes are associated with management practices that consequently result in habitat destruction and removal of organic substrate hence reduced availability of food sources for associated SIMFs. Forest land is relatively rich in terms of species of SIMF and high abundance was noted in the homestead garden fields. Of all the land use classes, crop cultivated outfields differed most from the others, recording relatively low taxonomic richness, and evenness. Therefore, practicing of sustainable cropping systems that maintain acceptable levels of SIMF abundance and that may minimize the level of human disturbances are desirable in this watershed to optimize their activities and their impacts on soil fertility. This requires the integration of knowledge of biological processes into the design of land management systems. More detailed studies are also needed to identify the  Signif. codes: 0 "***" 0.001, "**" 0.01, "*" 0.05, "." 0.