Next Article in Journal
Phytoremediating a Wastewater-Irrigated Soil Contaminated with Toxic Metals: Comparing the Efficacies of Different Crops
Previous Article in Journal
Soil Respiration Is Influenced by Seasonality, Forest Succession and Contrasting Biophysical Controls in a Tropical Dry Forest in Northwestern Mexico
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Aboveground Vegetation on Abundance, Diversity, and Biomass of Earthworms in Selected Land Use Systems as a Model of Synchrony between Aboveground and Belowground Habitats in Mid-Himalaya, India

1
School of Biological and Environmental Sciences, Shoolini University of Biotechnology and Management Sciences, Solan 173229, India
2
Department of Applied Biology, Adama Science and Technology University, Adama 1888, Ethiopia
3
Department of Zoology, Guru Nanak Dev University, Amritsar 143005, India
4
River Basin Research Center, Gifu University, 1-1 Yanagido, Gifu 501-1193, Japan
5
School of Agricultural Sciences, Shoolini University of Biotechnology and Management Sciences, Solan 173229, India
6
Institute of Animal Reproduction and Food Research, Polish Academy of Sciences, 10-748 Olsztyn, Poland
7
Department of Statistics, College of Natural and Computational Sciences, Addis Ababa University, Addis Ababa 1000, Ethiopia
*
Authors to whom correspondence should be addressed.
Soil Syst. 2022, 6(4), 76; https://doi.org/10.3390/soilsystems6040076
Submission received: 7 June 2022 / Revised: 29 July 2022 / Accepted: 29 September 2022 / Published: 4 October 2022

Abstract

:
The population status and biomass of earthworms were studied in three different land use systems of pasture (Pa), silvopasture (SP), and mixed evergreen forest (MEF) from 2019–2020 in the Solan district of Himachal Pradesh, India. The aim of this study was to assess the population status of earthworms and investigate how different land use systems influence their abundance, diversity, and biomass. Earthworms and soil were sampled using the Tropical Soil Biology and Fertility (TSBF) method in all seasons (winter, spring, summer, monsoon, and autumn). The physicochemical properties of the soil were analyzed to evaluate their effects on the diversity, biomass, and density of animals. The diversity status parameters, such as the Shannon diversity index (H′), Margalef richness index (R), evenness (J′), and dominance index (D), were computed. A total of seven earthworm species, belonging to four families, namely, Amynthas corticis, Aporrectodea rosea, Drawida japonica, Eisenia fetida, Metaphire birmanica, Metaphire houlleti, and Lennogaster pusillus, were identified from all three land use systems. The lowest Shannon diversity index (H′), Margalef index (R), and evenness (J′) index values were registered in MEF (H′ = 0.661, R = 0.762, J′ = 0.369) compared to those in Pa (H′ = 1.25, R = 1.165, J′ = 0.696) and SP (H′ = 0.99, R = 0.883, J′ = 0.552), implying that MEF is the least diversified land system. In contrast, the highest dominance index (D) value was registered in MEF (Pa = 0.39, SP = 0.53, MEF = 0.67), which again showed that MEF is the least diversified land system. The highest values of abundance and biomass were recorded in MEF (754.15 individuals m−2 and 156.02 g m−2), followed by SP (306.13 individuals m−2 and 124.84 g m−2) and Pa (77.87 individuals m−2 and 31.82 g m−2). Both the density and biomass of earthworms increased from Pa to MEF (Pa < SP < MEF). This study is novel because it revealed that the diversity and productivity (biomass and abundance) values of earthworms were negatively correlated (as diversity increased, productivity decreased; as diversity decreased, productivity increased). The total values of abundance and biomass of earthworms in the three land use systems indicated perfect synchrony between aboveground and belowground habitats, whereas the diversity values revealed that MEF was dominated by one or two species and the least diversified. Therefore, for sustainable belowground productivity, aboveground conservation is recommended, and vice versa, regardless of diversity.

1. Introduction

Land use change greatly influences diversity, distribution, and abundances of earthworms [1,2,3]. The effect of the plant community on underground soil faunal diversity remains hazy. As a result, understanding the relations among aboveground vegetation and living organisms in soil is of significant interest. The dynamics of earthworm populations in pineapple agroecosystems were evaluated to assess the effect of monocultures on earthworm communities [4]. They found that the community of earthworms in pineapple agroecosystems predominantly consisted of endogeic earthworms. Among all the endogeic species, Drawida assamensis was more dominant with respect to its biomass, density, and relative abundance. In addition, ref. [5] studied the population dynamics of earthworms and found that different land use systems greatly influenced species richness, community organization, abundance, and distribution.
The community structure of both fungi and earthworms was evaluated in four different land use systems, and their community structure was found to change in an increasing manner from conventional arable land to woodland (conventional arable land → no or reduced tillage → grassland → wooded land) [6]. The density and biomass of earthworm populations were examined in a wide range of land use types, such as moderately degraded natural forests, highly degraded natural forests, rehabilitated forest lands, traditional pure crop system, and traditional agroforestry system, and abandoned and rehabilitated farmland in a village landscape of the central Himalayas, India [7]. They reported that the change from the traditional pure crop system to the traditional agroforestry system significantly increased earthworm density and biomass. Both land use intensity and land use type are strong drivers of the abundance and composition of earthworm communities in agricultural ecosystems [8]. They recorded the lowest number of earthworms in coniferous forests and intensively managed agricultural land use, an intermediate number of populations in organic no-till systems and the highest populations in ancient deciduous forest systems.
Studies have reported that changes in land use systems impact soil nutrient dynamics in Amazonia [9] and influence soil physicochemical parameters [10,11]. Earthworms are keystone species of the soil ecosystem, and their distribution, diversity, and abundance assist the total diversity, physicochemical properties, and fertility of the soil [12]. On the other hand, tree plantations affected the abundance of earthworms by altering the physicochemical properties of the soil, such as temperature, humidity, pH, and organic matter content [13,14]. It was also reported that the overall density and biomass of earthworms are higher in older forests due to the accumulation of organic matter in the area occurring through the long-term breakdown of plant and animal material [5]. Moisture has been identified to influence the diversity and abundance of earthworms that are directly affected by temperature. As the temperature increases, the moisture content decreases, which does not allow the survival, growth, and reproduction of earthworms [4]. In studies on earthworm communities in rubber plantations of different ages in West Tripura, India, the population density and total biomass of earthworms increased during the monsoon and post-monsoon seasons, which were characterized by higher humidity, which is suitable for the survival, growth, and reproduction of earthworms [15]; it was noted that both density and biomass increased, while species diversity, species richness, and species evenness decreased with increasing farm age. Compared with forest ecosystems, the diversity of functional guilds was lower in agricultural ecosystems [16].
No such studies have been carried out to reveal how aboveground vegetation influences the physicochemical parameters of soil and in turn directly affects belowground diversity. In addition, there have been no sufficient studies on the belowground diversity of the Himalayan mountain system, which is recognized as a global biodiversity hotspot [17] and requires extensive study. This study aimed to analyze the species diversity, biomass, and density of earthworms in relation to land use systems of pasture (Pa), silvopasture (SP), and mixed evergreen forest (MEF) in the Mid-Himalaya of the Solan district, India. In this study, we used earthworms (keystone species of belowground habitat) as a model to investigate (a) how aboveground vegetation of Pa, SP, and MEF influenced belowground population status (density, biomass, and diversity) of earthworms; (b) alpha diversity and the trend of diversity indices with both biomass and density values in the land use systems. The outcomes of this research could be useful to provide insight on the trend of diversity and productivity (density and biomass) in the soil of the studied land use systems determined by habitat preference of earthworms based on adaptations that may need further investigation.

2. Materials and Methods

The sampling was conducted during the winter, spring, summer, monsoon, and autumn seasons at three different sites, viz., Pa, SP, and MEF land use systems, selected for this study in the same locality of a village called Chemalti, Solan, Himachal Pradesh, India (Figure 1). Solan has an average elevation of 1502 metres (5249.34 feet) at GPS coordinates of 30.9084° N, 77.0999° E and its average annual rainfall is 1262 mm. The GPS coordinates of the sampling site are 30.8628° N, 77.097° E. These different land use systems were characterized by having the following vegetation. The pastoral land use system has Chrysopogon fulvus, Heteropogon contortus, and Apluda spp. Silvopastoral land use systems possess Quercus leucotricophora (Ban Oak), shrubs, and grasses. The mixed evergreen forestland use system comprises Quercus leucotricophora (Ban Oak), Pinus roxburghii (Chir Pine), Myrica nagi (Kafal), Prunus paddus (Paja), Celtis australis (Khirak), Toona ciliata (Toon), Ficus spp. (Anjeer), Indigofera pulchella (Kathi), Beberis asiatica (Rasaunt/Kasgmal), Rubus ellipticus (Yellow Himalayan Raspberry), and Carrisa coronda (Coronda). The soil of the current study area is categorized as the loam soil type.

2.1. Earthworm Sampling

Earthworms were sampled using the Tropical Soil Biology and Fertility (TSBF) method [18,19]. Soil monoliths (25 × 25 × 30 cm) were taken by digging the soil. Three spatially different plots (10 × 10 m) were marked in each land use system (Table 1). Five soil monoliths were randomly taken from each plot, resulting in a total of 15 samples from each land use type. The samples were collected for five seasons (Table 2). All samplings were carried out during different seasons over a period of two years (February 2019–November 2020). For each land use system (Pa, SP, and MEF), triplicates were used (Pa1, Pa2, Pa3; SP1, SP2, SP3; MEF1, MEF2, and MEF3) to ensure statistically significant data.

2.2. Identification of Earthworms

Earthworms were hand-sorted, and the collected earthworms from the sites were labeled with the place of collection (name of site, sampling date) as shown in Figure 2. The collected earthworm samples were thoroughly washed to remove soil debris and preserved in 5% formalin to avoid color loss, as color is an important criterion for earthworm species identification. The earthworms were identified to the species level with the help of keys provided by [20,21,22,23,24,25,26,27,28].
Reasonable amounts of soil were taken from the place where earthworms were hand-sorted in plastic bags and were labeled with the site name and the date of sampling. The physicochemical parameters of the soil samples, such as pH, temperature, moisture, and available carbon, nitrogen, phosphorous, and potassium, were analyzed using different analysis methods. The soil temperature of the different land use types in each plot was immediately measured by a soil thermometer. Soil texture was determined by the textural triangle of the soil classification system defined by the USDA [29], and moisture was detected by the gravimetric method [30]. Soil pH was measured using a digital pH meter (Deluxe pH meter 101). Available organic carbon was analyzed as by [31]. Nitrogen was analyzed as by [32]. Phosphorus was analyzed as by [33], and potassium was analyzed as by [34].

2.3. Statistical Analysis

Direct enumerated numbers of earthworms were input for calculation of density per m−2. The biomass of earthworms per m−2 was calculated from weight measurement of preserved wet earthworms by a digital measuring balance. The mean ± SD of density, mean ± SD biomass, and relative density of earthworms were computed. The density and biomass of earthworms per m−2 were calculated by converting the area of one monolith into m−2. The area of a monolith (25 × 25 cm = 625 cm−2) was converted into 1 m−2 (10,000 cm−2). The conversion factor (16) was determined by dividing 10,000 cm−2 to 625. In one sampling site, there were three plots and from each plot five monoliths were dug up to extract the samples, making a total of 15 monoliths. To calculate density or biomass of earthworms per m−2, the actual counts or measurements of earthworms were multiplied by 16 (conversion factor) and divided to 15 (number of monoliths in one site). The relative density was determined by dividing density (ind·m−2) of each species to the total density of species in each site and multiplied by one hundred. Stata version 14 statistical software was applied to analyze variations in soil nutrient data (C, N, P, and K) between land use systems and seasons by two-way ANOVA.
Index of species richness was estimated by the Margalef species richness index [35].
R = S 1 I n N ,
where R is species richness, S is total number of species, In is natural logarithm of number of species, and N is total number of individuals in the sample.
The index of evenness was quantified by the Pielou evenness index formula [36].
J = H I n s ,
where J′ is evenness index, H′ is Shannon diversity index, s is total number of species in the sample, and In is natural logarithm of number of species.
Species diversity was calculated by the Shannon–Wiener index [37]
H = i = 1 S pi   I n pi ,
where H′ is Shannon diversity index, pi is relative proportion of species, and In is natural logarithm of the number of species.
The index of dominance was calculated by using the Simpson diversity index formula [38]
D = N ( N 1 ) n ( n 1 ) ,
where D is Simpson’s index, n is the total number of organisms of a species, and N is the total number of organisms of all species.

3. Results

3.1. Physicochemical Analysis of Soil

A total of 90 soil samples were analyzed to detect the physicochemical properties of the three land use systems. The data in Table 3 revealed that the soil texture of the three land use systems fell under the same category of loam soil based on the determination of the textural triangle of the soil classification system defined by [29]. The soil temperature (°C) was greater in pasture (Pa), where there was less vegetation canopy that prevents direct sunlight from reaching the soil. The temperature values become lower in SP and much lower in MEF. The soil moisture in Pa was lower, and inversely proportional to the soil temperatures that become higher in SP and MEF. As the decomposition process proceeded, the pH of the soil slightly increased from Pa to SP and MEF. Soil organic carbon (%), available nitrogen (kg/ha), and available phosphorus (kg/ha) also increased from Pa to SP and MEF. The available potassium (kg/ha) was much lower in Pa than in SP and MEF, but its concentration in SP was slightly greater than that in MEF (Table 3, Figure 3 and Figure 4).
As indicated in Table 4, soil temperature and moisture were negatively correlated in both winter and spring (as the soil temperature value increased, the moisture values decreased). During three seasons, viz., summer, monsoon, and autumn, the values of both parameters were positively correlated as presented (as the temperature values decreased or increased, moisture values also decreased or increased). However, land use system-wise, as the temperature values increased the values of moisture decreased (Figure 3 and Figure 4).
The two-way ANOVA output showed that the variation in organic carbon (C) concentration, potassium (K), phosphorus (P), and nitrogen (N) between land use systems were statistically significant at p < 0.0001. The variation in soil nutrients (C, K, and P) between seasons was also statistically significant (p < 0.005) except for N (p < 0.5) (Table 5).
The data of the concentration of nutrients (C, K, P, and N) presented in Figure 5a–d between land use systems revealed less in Pa, an increase in SP, and a great increase in MEF, respectively.
Pair-wise post hoc analyses of soil nutrients were performed to indicate the comparison of soil nutrient (C, K, P, and N) data variation between land use systems (Table 6).
The concentration of nutrients (C, K, P, and N) in different seasons indicated that the concentrations of all nutrients were greater in monsoon, except for nitrogen which was slightly greater in spring than other nutrients (Figure 6a–d).

3.2. Key Used for the Identification of Earthworms

The preserved specimens of earthworms were identified following the key provided by [21,22,23,24,25,26,27,28].
  • Genital markings present, spermathecal atrium pear or finger—Drawida japonica.
  • Pale- or light pink-colored worms, calciferous sacs present in segment 10—Aporrectodea rosea.
  • Worms that are red-colored with yellow transverse stripes, calciferous sacs absent in segment 10—Esenia fetida.
  • Ventral most chaetae on segment 8 not enlarged—Lennogaster pusillus.
  • Spermathecal pores in intersegmental furrows 5/6/7/8; genital markings absent—Metaphire birmanica.
  • Spermathecal pores in furrows 6/7/8/9; genital markings small within copulatory pouches and spermathecal pore invaginations, recognizable internally by the presence of stalked glands—Metaphire houlleti.
  • Spermathecal pores, four pairs, at intersegmental furrows 5/6/7/8/9—Amynthas corticis.

3.3. Density and Biomass of Earthworms

Generally, the mean density of earthworms during the study period of 2019–2020 in Pa, with a value of 77.87 ind·m−2, was lower than mean earthworm density in the other two land use systems. In SP, the mean density, with a value of 306.13 ind·m−2, was much higher than the mean density of Pa, and the highest density was observed in MEF, with a value of 754.15 ind·m−2. Season-wise, the highest mean density (ind·m−2) of earthworms was generally observed in MEF, except in winter. The lowest mean density (ind·m−2) was observed in Pa in all seasons. During winter, the mean density of earthworms was highest in SP, with a value of 87.47 ind·m−2 (Table 7 and Table 8).
The total mean biomass of earthworms was similarly lowest in Pa, with a value of 31.82 g m−2, but in SP, the total mean biomass was much higher than that in Pa, with a value of 124.84 g m−2, and the total mean biomass of MEF was higher than the total mean biomass of SP, with a value of 156.02 g m−2 (Table 7 and Table 8).

3.4. Species-Wise Biomass and Density Status in Different Land Use Systems

Both mean biomass (g m−2) and mean density (ind·m−2) values in MEF were higher than those in Pa and SP. In Pa and SP, A. corticis was denser, and its biomass was much greater than that of the other species. However, in MEF, A. rosea was the densest and most abundantly occurring earthworm species. In the present study, E. fetida was absent in Pa, and L. pusillus occurred only in Pa and was absent in both SP and MEF (Table 8).

3.5. Species-Wise Density and Relative Density Status in Different Land Use Systems

The most dominant earthworm species in Pa was A. corticis, with a value of 45.87 ind·m−2 and RD of 58.91%. D. japonica was the next most abundant in Pa, with a value of 12.8 ind·m−2 and RD of 16.44%, followed by M. birmanica and L. pusillus, with abundance values and RD of 6.4 ind·m−2 and 8.22% and 3.2 ind·m−2 and 4.11%, respectively. M. houlleti occurred with 1.07 ind·m−2 and RD of 1.37%. E. fetida was totally absent in Pa (Table 9).
In SP, A. corticis was again the most abundant, with a density of 217.6 ind·m−2 and RD of 71.08%. The second most abundantly registered earthworm species was M. houlleti, with a density of 43.73 ind·m−2, followed by D. japonica with a density of 20.27 ind·m−2 and RD of 6.62% and E. fetida with a density of 9.6 ind·m−2 and RD of 3.14%. M. birmanica occurred in SP, with an abundance value of 8.53 ind·m−2 and RD of 2.79%, followed by A. rosea, with an abundance value of 6.4 ind·m−2 and RD of 2.09%, but L. pusillus did not occur at all in SP (Table 9).
In MEF, the most abundantly occurring earthworm species was A. rosea, with an abundance value of 605.87 ind·m−2 and RD of 80.34%. This earthworm species possesses a higher degree of adaptation to survive and reproduce in places where high humidity is available in the deep soil of MEF, as it is an endogeic species. The next most abundant earthworm species in MEF was A. corticis, with a density of 109.87 ind·m−2 and RD of 14.57%, followed by D. japonica, with an abundance value of 20.27 ind·m−2 and RD of 2.69%. The next most abundantly observed earthworm species was M. houlleti, with a density value of 12.8 ind·m−2 and RD of 1.69%, followed by E. fetida, with a value of 4.27 ind·m−2 and RD of 0.57%. M. birmanica was the least abundantly recorded species, with a value of 1.07 ind·m−2 and RD of 0.14% (Table 9).

3.6. Diversity Index and Abundance of Earthworm Species in Land Use Systems

The diversity indices were calculated for Pa, SP, and MEF to evaluate the diversity status of the three land use systems (Table 10). According to the output obtained, the lowest value of the Shannon diversity index was registered in MEF at 0.661 and in SP at 0.99, and the highest value was registered in Pa at 1.25. The lowest index of evenness was registered in MEF at 0.369, then higher in SP at 0.552, and the value was the highest in Pa at 0.696. The highest value of the Margalef richness index was registered in Pa (1.165), a lower value in SP (0.883), and the lowest value in MEF (0.762). Similarly, the dominance index was also computed for the three land use systems, and the highest value was registered in MEF at 0.67, then lower in SP with a value of 0.53, and the lowest value (0.39) was registered in Pa.

3.7. Synchrony of Density and Biomass across Land Use Systems

The densities of earthworms were related to the vegetation of the land use systems. As its name indicates, the Pa land use system possesses less dense aboveground vegetation of grass species such as Chrysopogon fulvus, Heteropogon contortus, and Apluda spp. with a little litter accumulation. These vegetation types do not have enough canopy to protect the soil from sunlight. Both temperature and moisture of soil play a vital role in the survival and reproduction of soil organisms. Direct sunlight allows evaporation of water from the soil and results in less moisture content, making land use systems less suitable for the survival and reproduction of earthworms. Hence, less moisture and high temperature in Pa resulted in lower density and biomass of earthworms (77.87 ind·m−2, 31.82 g m−2) (Table 8).
In SP, tree species such as Quercus leucotricophora (Ban Oak), shrubs, and grasses were available which provide more canopy to prevent direct sunlight from reaching the soil than in Pa. The temperature and moisture contents in SP were more suitable for survival and reproduction and hence the density and biomass of earthworms were greater (306.13 ind·m−2, 124.84 g m−2) (Table 8). MEF possesses more vegetation such as Quercus leucotricophora (Ban Oak), Pinus roxburghii (Chir Pine), Myrica nagi (Kafal), Prunus paddus (Paja), Celtis australis (Khirak), Toona ciliata (Toon), Ficus spp. (Anjeer), Indigofera pulchella (Kathi), Beberis asiatica (Rasaunt/Kasgmal), Rubus ellipticus (Yellow Himalayan Raspberry), and Carrisa coronda (Coronda) with much accumulation of litter. The vegetation provides more survival and reproductive advantages for earthworms and hence their density and biomass increased to 754.15 in.m−2, 156.02 g m−2 (Table 8).

3.8. Relationship of Density and Biomass with Diversity

The data of density and biomass of earthworms were negatively correlated with diversity values (Table 11). According to the recorded data in Pa (77.87 ind·m−2 and 31.82 g m−2), in SP (306.13 ind·m−2 and 124.84 g m−2), and in MEF (754.15 ind·m−2 and 156.02 g m−2), both density and biomass (productivities) of earthworms increased from Pa to MEF with the increasing aboveground vegetation (MEF > SP > Pa). However, the diversity index values of Shannon (H′), evenness (J′), and Margalef richness indices (R) increased in Pa and decreased in SP and were the lowest in MEF. The values decreased from Pa to MEF as indicated in Table 9, from Pa (H′ = 1.25, R = 1.165, J′ = 0.696), to SP (H′ = 0.99, R = 0.883, J′ = 0.552), to MEF (H′ = 0.661, R = 0.762, J′ = 0.369) (Pa > SP > MEF). On the other hand, the highest dominance index (D) value was registered in MEF (Pa = 0.39, SP = 0.53, MEF = 0.67). The values of alpha diversity revealed that in the places where there was less vegetation, in, e.g., Pa, high belowground diversity occurred and in the places where more vegetation was observed, e.g., in MEF, a lower diversity of earthworms was recorded.

4. Discussion

4.1. Physicochemical Analysis Results of Three Different Land Use Systems

The physicochemical parameters of soil are independent variables that affect the diversity, density, and biomass (dependent variables) of earthworms. The soil physicochemical parameters of different land use systems (Pa, SP, and MEF) are influenced by the absence and availability of aboveground vegetation [39]. There was an investigation that evaluated how changes in land use systems such as grassland, plantation forest, and natural forest influenced the decrement and increment of soil nutrients in Ethiopia and reported that in natural forest, the concentration of soil nutrients was greatest as compared to other land use systems [40]. The change from natural forest to farmland and to other types decreased soil nutrients such as organic carbon, total nitrogen, and available phosphorus [11]. Another study also reported that conversion of land use systems influenced soil nutrient dynamics in Amazonia [9] and other research revealed changes in land use systems able to affect soil physicochemical parameters [11]. According to the results presented in Table 4, the recorded mean ± SD of temperature in Pa (19.91 ± 2.27 °C) was the highest, lower in SP (18.17 ± 1.72 °C), and the lowest in MEF (16.87 ± 1.91 °C). As the temperature increases, the moisture decreases, and temperature influenced the percentage of moisture in the soil (Figure 3 and Figure 4). The Pa land use system possessed a lower moisture content (21.08 ± 3.59%), as a higher temperature allowed loss of water from the soil through evaporation. The moisture content of the soil in the SP land use system increased (25.47 ± 4.18%) and it was the highest in MEF (29.47 ± 5.69%) (Table 3, Figure 3 and Figure 4).
It was reported that the functional composition and species diversity of plants (above ground) affected the concentration of organic carbon, available nitrogen, and phosphorus in the soil [41]. Accordingly, the concentrations of organic carbon (%), available nitrogen (kg/ha), available phosphorus (kg/ha), and available potassium (kg/ha) in Pa were lower than those in SP and MEF. In Pa, there was less dense aboveground vegetation, resulting in higher temperature (°C) and lower moisture contents (%), which in turn affected the decomposition of organic matter and the availability of these nutrients. In this study, except for available potassium (K) (kg/ha), all soil nutrients (C, N, P) analyzed were highest in MEF. The available potassium (kg/ha) of SP was much greater than that of Pa and slightly greater than that of MEF (Table 3). The data in Figure 5a–d show that the concentrations of nutrients (C, K, P, and N) were lower in Pa, greater in SP, and the greatest in MEF, respectively. The variations in nutrients between land use systems were statistically significant as presented in Table 5 and with a pair-wise post hoc test (Table 6).

4.2. Density and Biomass of Earthworms

The density (ind·m−2), relative density (RD%), and biomass (g m−2) of all earthworms are presented in Table 7, Table 8 and Table 9. The MEF comprised tree species such as Quercus leucotricophora (Ban Oak), Pinus roxburghii (Chir Pine) Myrica nagi (Kafal), Prunus paddus (Paja), Celtis australis (Khirak), Toona ciliata (Toon), and Ficus spp. (Anjeer). The shrub species were Indigofera pulchella (Kathi), Beberisa sciatica (Rasaunt/Kashmal), Rubus ellipticus (Hinr), and Carrisa coronda (Coronda) with much accumulation of leaf litter. Higher densities of these tree and shrub species allowed a higher density of earthworms to occur in MEF. It was found that changes in land use systems resulted in changes in the earthworm community [42]. Earthworm density was mainly influenced by aboveground vegetation [43], which confirmed that the aboveground vegetation supported belowground organisms and vice versa.
The higher density of aboveground vegetation directly influenced the higher density and biomass of belowground organisms by supplying food resources and stabilizing the physicochemical properties of the soil. Canopies of MEF serve to provide shade and conserve moisture in the soil, stabilize temperature, and facilitate the breakdown of organic matter, which supports the survival of belowground organisms [44]. The canopies protect the soil from direct sunlight, creating conditions that are conducive to the survival and reproduction of belowground organisms [13]. Accordingly, the highest density of earthworms was observed in MEF, where a higher density of vegetation was observed as compared to other land use systems (Pa and SP). In addition, earthworm species in the three land use systems (A. rosea, M. houlleti, M. birmanica, D. japonica, and L. pusillus) were endogeic. E. fetida is epigeic and A. corticis is epi-endogeic [45]. Therefore, MEF was the preferred habitat for most species of the study sites. The effect of land use systems having large amounts of aboveground vegetation was confirmed by the data presented in Table 8. The highest mean density (754.13 ind·m−2) of earthworms was recorded in MEF, followed by SP with a value of 306.13 ind·m−2, and the lowest density in Pa with a value of 77.87 ind·m−2.
The variation in the values of the three land use systems was due to the variation in physicochemical parameters caused by the shade provided from vegetation of these different land systems and availability of resources. In the Pa land use system, the vegetation was limited to grass species such as Chrysopogon fulvus, Heteropogon contortus, and Apluda spp., with less litter accumulation, and therefore, the land use system was limited in resources, which resulted in limited density.
As presented in Table 8, the mean biomass of earthworms was lowest, with a value of 31.82 g m−2 in Pa compared to the other land use systems (SP with a value of 124.84 g m−2 and MEF with a value of 156.02 g m−2), as in the case of the mean density. The Pa land use type was devoid of resources due to the lower density of aboveground vegetation, resulting in physicochemical conditions being less suitable for the survival and reproduction of earthworms. SP was more suitable than Pa and subsequently resulted in more biomass. Most earthworm species, such as A. rosea, D. japonica, M. houlletii, and M. birmanica (endogeic), were more adapted to survive in a soil habitat where large amounts of leaf litter were deposited in MEF. The result revealed that MEF was the most suitable land use system for the survival and reproduction of the worms and this caused the biomass to increase. The values of density and biomass of earthworms were high in monsoon (Table 7) and were directly related to the highest values of soil nutrients (Figure 6a–d). Monsoon is the rainy season of north India characterized by highest moisture contents of soil that allow decomposition and availability of soil nutrients to support the survival and reproduction of earthworms. The lowest value of soil nutrients was generally recorded in Pa, while it was greater in SP and the greatest in MEF (Table 3, Figure 5a–d). These values directly influenced the survival and reproduction of earthworms and subsequently affected density and biomass.
According to one study, overall density and biomass of earthworms in different land use systems were determined by plant diversity and soil properties. They divided hilly land use types into annual dominated (AD), mixed (MX), and perennial dominated (PD) and investigated how aboveground and soil properties affected density and biomass of earthworms in the USA and found that both biomass and density of earthworms in mixed (MX) land use were significantly greater than the two land use systems in which plant diversity was the greatest [46]. Greater diversity of vegetation facilitates soil physicochemical properties, conducive for survival and reproduction of earthworms. In addition, [8] investigated the abundance of earthworms in different land use systems such as coniferous forests, intensively managed agricultural land use, organic no-till systems, and ancient deciduous forest systems. The result revealed that the highest population of earthworms was recorded in ancient deciduous forest systems.

4.3. Diversity of Earthworm Species

A total of seven species of earthworms were identified in the three land use systems. In both Pa and SP, A. corticis (epi-endogeic) was dominant, but in MEF, A. rosea (endogeic) was the dominant species (Table 10). Adaptation (pigmentation) helps A. corticis to survive in sunlight-exposed habitat, and it occurred more abundantly in Pa than other species [47]; again, due to its adaptation to dwell in the soil epi-endogeically, A. corticis occurred more abundantly in MEF alongside A. rosea. A. rosea is exclusively endogeic and best adapted to survive in forests where leaf litter accumulates, and rarely occurred in Pa.
Diversity indices were computed for the three land use types separately to compare the diversity status (Table 8). Lower values of the Shannon diversity index (H′), Margalef species richness index (R), and evenness index (J′) imply a lower diversity in the given area. Among the three land use systems, the Shannon diversity index value was lowest in MEF, higher in SP, and the highest in Pa, with values of 0.661, 0.99, and 1.25, respectively, indicating that MEF was the least diversified land use system. Evenness (J′) was also applied, and the lowest evenness value was registered in MEF, a higher value in SP, and the highest value in Pa of 0.369, 0.552, and 0.696, respectively. The lowest value of evenness also indicated the lowest diversity. Another diversity parameter value was also computed to determine species richness by the Margalef richness index (R) with its highest value in Pa (1.165), then SP (0.883) and MEF (0.762). Based on these diversity index values, MEF was the least diversified land use system. In contrast, the greatest Simpson dominance index was registered in MEF with a value of 0.67, while it was lower in SP with 0.53, and lowest in Pa with a value of 0.39 (MEF > SP > Pa), implying again that MEF was less diversified. The diversity values implied that MEF was dominated by one or two species of earthworms that developed adaptations to dwell in the soil epi-endogeically (A. corticis) and endogeically (A. rosea). The study revealed that the dominance of certain species of earthworms in a land use system decreased the diversity index value and made the area less diversified. The area which was dominated by Metaphire tschiliensis greatly reduced the diversity status of the land use system [48].
This study confirmed that the diversity of earthworms was not positively correlated with density of vegetation of land use systems (Pa, SP, and MEF). The density and biomass of earthworms were much lower in the land use system with less dense vegetation, while they increased as the density of vegetation increased [49]. In the three land use systems, the species richness (n) was the same, which matched with the work of [7] which reported that changes in land use systems do not necessarily cause a change in species richness in the same landscape. A study investigated the population status of earthworms in different land use systems and found that the change from the traditional pure crop system to the traditional agroforestry system significantly increased earthworm density and biomass but not species richness in a landscape of the central Himalayas, India [7]. Another study [15] investigated the effect of land use system conversion and reported that density and biomass of earthworms increased in aged farm land, while species diversity, richness, and evenness decreased. According to this study, the farmland was protected from human activities and became more suitable for the survival of certain species of earthworms only, causing the land to be dominated by a few species and minimizing its diversity.

4.4. Trends of Productivity (Density and Biomass) and Diversity across Land Use Systems

Alpha diversity and productivity (biomass and density) values of earthworms in different land use systems revealed that factors that influence earthworm diversity differ from factors that influence their density and biomass. According to the result of this study, MEF was dominated by two species of earthworms (A. corticis and A. rosea) which made the area less diversified, implying dominancy reduces diversity. The alpha diversity index value indicated that the Pa land use system was the most diversified, SP was less diversified, and MEF was the least diversified. Both earthworm species (A. corticis and A. rosea) developed adaptations to survive in the soil epi-endogeically and endogeically, respectively, and were able to survive and reproduce better than other earthworm species in MEF. They possess very little or no cutaneous pigment to survive in less intense light and are able to feed directly on soil materials. These adaptations enable them to survive in the soil and become more dominant in the place where many litter layers accumulated and made MEF less diversified (Table 11).
The productivity (biomass and density) values of earthworms in the three land use systems were the lowest in the Pa land use system, greater in SP, and greatest in MEF (Table 8 and Table 9). In the Pa land use system, the soil temperature value was the highest which does not allow moderate moisture due to a smaller canopy of vegetation. In SP, more vegetation was available to provide a canopy and prevented direct sunlight from reaching the soil, resulting in a more suitable temperature and moisture than in Pa. MEF was the most suitable habitat for the survival and reproduction of most earthworms and hence had more productivity [47]. According to the finding of this study, the trend of productivity decreased from MEF to Pa (MEF > SP > Pa) (Table 8) and diversity increased from MEF to Pa (Pa > SP > MEF) (Table 11). Therefore, temperature and moisture contents of land use systems are the driving factors behind determination of biomass and density of earthworms.

5. Conclusions

The density, diversity, and biomass of earthworms were studied over a period of two years from 2019–2020 in Solan district, Himachal Pradesh, India, in three different land use systems, Pa, SP, and MEF, to evaluate how the change in land use systems impacted the population status of earthworms. The lowest values of H′, R, and J′ were registered in MEF, indicating that it was the least diversified land system. In contrast, the highest value of the dominance index (D) in MEF also revealed that MEF was the least diversified land use system. Among the seven earthworm species identified from the three land use systems, A. corticis was dominant in both Pa and SP, but A. rosea was the dominant species in MEF. The survival of A. corticis and A. rosea depended on the adaptation of habitat preference as epi-endogeic and endogeic, respectively. The species richness (n) in the three land use systems was the same. The mean density and mean biomass of earthworms were lower in Pa, where there was a lower density of aboveground vegetation, and the density and biomass of earthworms subsequently increased in SP and MEF, where the density of aboveground vegetation was the highest. The four measures of alpha diversity (H′, J′, R, and D) showed that MEF was the least diversified land use system, implying that the productivity (biomass and density) values were negatively correlated with the diversity values of the soils in the land use systems (productivity increased while diversity decreased). The results showed that the productivity (biomass and abundant) of earthworms was greatly influenced by changes in land use systems and showed perfect synchrony between aboveground and belowground habitat. Therefore, it is recommended that for belowground productivity, the conservation of aboveground habitat should not be neglected and vice versa.

Author Contributions

K.G.: Conceptualization, Visualization, Software, Methodology, Validation, Writing—review and editing; A.B.: Visualization, Software, Methodology, Validation; J.M.J.: Visualization, Software, Methodology, Validation; S.A.B.: Visualization, Software, Methodology, Validation, Funding; G.K.S.: Visualization, Software, Methodology, Validation; G.B.: Visualization, Software, Methodology, Validation, Writing—review and editing; R.A.: Writing—review and editing; M.B.: Software, Methodology, Validation, Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

The research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Shoolini University for facilitating transportation during sampling in each season of two years and supplying lab facilities for soil analysis and species identification.

Conflicts of Interest

The authors declared no conflict of interest.

References

  1. Whalen, J.K. Spatial and temporal distribution of earthworm patches in corn field, hayfield and forest systems of southwestern Quebec, Canada. Appl. Soil Ecol. 2004, 27, 143–151. [Google Scholar] [CrossRef]
  2. López-Hernández, D. Earthworm populations in savannas of the Orinoco basin. A review of studies in long-term agricultural-managed and protected ecosystems. Agriculture 2012, 2, 87–108. [Google Scholar] [CrossRef]
  3. Bartz, M.L.; Brown, G.G.; da Rosa, M.G.; Klauberg Filho, O.; James, S.W.; Decaëns, T.; Baretta, D. Earthworm richness in land-use systems in Santa Catarina, Brazil. Appl. Soil Ecol. 2014, 83, 59–70. [Google Scholar] [CrossRef]
  4. Dey, A.; Chaudhuri, P.S. Species richness, community organization, and spatiotemporal distribution of earthworms in the pineapple agroecosystems of Tripura, India. Int. J. Ecol. 2016, 2016, 3190182. [Google Scholar] [CrossRef] [Green Version]
  5. Bisht, R.; Pandey, H.; Bharti, D.; Kaushal, B.R. Population dynamics of earthworms (Oligochaeta) in cultivated soils of Central Himalayan tarai region. Trop. Ecol. 2003, 44, 227–232. [Google Scholar]
  6. Spurgeon, D.J.; Keith, A.M.; Schmidt, O.; Lammertsma, D.R.; Faber, J.H. Land-use and land-management change: Relationships with earthworm and fungi communities and soil structural properties. BMC Ecol. 2013, 13, 46. [Google Scholar] [CrossRef] [Green Version]
  7. Bhadauria, T.; Kumar, P.; Kumar, R.; Maikhuri, R.K.; Rao, K.S.; Saxena, K.G. Earthworm populations in a traditional village landscape in Central Himalaya, India. Appl. Soil Ecol. 2012, 53, 83–93. [Google Scholar] [CrossRef]
  8. Smith, R.G.; McSwiney, C.P.; Grandy, A.S.; Suwanwaree, P.; Snider, R.M.; Robertson, G.P. Diversity and abundance of earthworms across an agricultural land-use intensity gradient. Soil Tillage Res. 2008, 100, 83–88. [Google Scholar] [CrossRef]
  9. McGrath, D.A.; Smith, C.K.; Gholz, H.L.; Oliveira, F.D. Effect of land-use change on soil nutrient dynamics in Amazonia. Ecosystems 2001, 4, 625–645. [Google Scholar] [CrossRef]
  10. Haokip, S.L.; Singh, T.B. Diversity and distribution of earthworms in a natural reserved and disturbed sub-tropical forest ecosystem of Imphal-West, Manipur, India. Int. Multidiscip. Res. J. 2012, 2, 28–34. [Google Scholar]
  11. Tellen, V.A.; Yerima, B.P. Effect of land use change on soil physicochemical properties in selected areas in the Northern West region of Cameroon. Environ. Syst. Res. 2018, 7, 3. [Google Scholar] [CrossRef] [Green Version]
  12. Rombke, J.; Jansch, S.; Didden, W. The use of earthworms in ecological soil classification and assessment concepts. Ecotoxicol. Environ. Saf. 2005, 62, 249–265. [Google Scholar] [CrossRef] [PubMed]
  13. Gonzalez, G.; Zou, X.; Borges, S. Earthworm abundance and species composition in abandoned tropical croplands: Comparisons of tree plantations and secondary forests. Pedobiologia 1996, 40, 385–391. [Google Scholar]
  14. Melnichuk, R.D.; Tesimen, H.B.; Gorres, J.H. Do the Invasive Earthworms Amynthas agrestis (Oligochaeta: Megascolecidae) and Lumbricus rubellus (Oligochaeta: Lumbricidae) Stimulate Oxalate-Based Browser Defenses in Jack-in-the-Pulpit (Arisaema triphyllum) by Their Presence or Their Soil Biogeochemical Activity? Soil Syst. 2022, 6, 11. [Google Scholar]
  15. Chaudhuri, P.S.; Bhattacharjee, S.; Dey, A.; Chattopadhyay, S.; Bhattacharya, D. Impact of age of rubber (Hevea brasiliensis) plantation on earthworm communities of West Tripura (India). J. Environ. Biol. 2013, 34, 59. [Google Scholar] [PubMed]
  16. Sinha, B.; Bhadauria, T.; Ramakrishnan, P.S.; Saxena, K.G.; Maikhuri, R.K. Impact of landscape modification on earthworm diversity and abundance in the Hariyali sacred landscape, Garhwal Himalaya. Pedobiologia 2003, 47, 357–370. [Google Scholar] [CrossRef]
  17. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; da Fonseca, G.A.B.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef] [PubMed]
  18. Anderson, J.M.; Ingram, J.S. A handbook of methods. CAB Int. Walling Oxf. Shire 1993, 221, 62–65. [Google Scholar]
  19. Swift, M.; Bignell, D. Standard Methods for Assessment of Soil Biodiversity and Land Use Practice; International Centre for Research in Agroforestry: Bogor, Indonesia, 2001. [Google Scholar]
  20. Julka, J.M. The Fauna of India and the Adjacent Countries. Megadrile Oligochaeta (Earth Worms). Haplotaxida: Lumbricina: Megascolecoidea: Octochaetidae; Zoological Survey of India: Calcutta, India, 1988; 400p. [Google Scholar]
  21. Julka, J.M.; Senapati, B.K. Earthworms (Oligochaeta: Annelida) of Orissa, India. Rec. Zool. Surv. India 1987, 92, 1–44. [Google Scholar]
  22. Gates, G.E. On a taxonomic puzzle and the classification of the earthworms. Bull. Mus. Comp. Zool. Harv. 1959, 121, 229–261. [Google Scholar]
  23. Gates, G.E. Burmese earthworms, an introduction to the systematic and biology of the megadrile oligochaetes with special reference to Southeast Asia. Trans. Am. Philos. Soc. 1972, 62, 1–326. [Google Scholar] [CrossRef]
  24. Blakemore, R.J. Checklist of Indian Earthworms. A Series of Searchable Texts on Earth-Worm Biodiversity, Ecology and Systematics from Various Regions of the World 3rd Edition. 2007. Available online: http://www.annelida.net/earthworm/Indian.pdf (accessed on 15 February 2022).
  25. Blakemore, R.J. Cosmopolitan Earthworms—An Eco-Taxonomic Guide to the Peregrine Species of the World, 5th ed.; Verm Ecology Solutions: Yokohama, Japan, 2012. [Google Scholar]
  26. Ahmed, S.; Julka, J.M.; Kumar, H. Earthworms (Annelida: Clitellata: Megadrili) of Solan, a constituent of Himalayan Biodiversity Hotspot, India. Trav. Muséum Natl. D’histoire Nat. “Grigore Antipa” 2020, 63, 19–50. [Google Scholar] [CrossRef]
  27. Csuzdi, C. Earthworm species, a searchable database. Opuscula Zool. Budapest. 2012, 43, 97–99. [Google Scholar]
  28. Csuzdi, C.; Szederjesi, T.; Sherlock, E. Annotated checklist of earthworm species de- scribed by András Zicsi (Clitellata: Megadrili). Zootaxa 2018, 4496, 11–42. [Google Scholar] [CrossRef]
  29. Devis, R.O.E.; Bennett, H.H. Grouping of Soil on the Basis of Mechanical Analysis; Dep.Circ.419; U.S. Department of Agriculture: Washington, DC, USA, 1927. [Google Scholar]
  30. Reynolds, S.G. Gravimetric method of soil moisture determination Part IA study of equipment, and methodological problems. J. Hydrol. 1970, 11, 258–273. [Google Scholar] [CrossRef]
  31. Walkley, A.; Black, I.A. An examination of the Degtiareff method for determining soil organic matter and proposed modification of the chromic acid titration method. Soil Sci. 1934, 63, 29–38. [Google Scholar] [CrossRef]
  32. Jackson, M.L. Soil Chemical Analysis; Constable & Co., Ltd.: London, UK, 1962. [Google Scholar]
  33. Bray, R.H.; Kurtz, L.T. Determination of total, organic, and available forms of phosphorus in soils. Soil Sci. 1945, 59, 39–45. [Google Scholar] [CrossRef]
  34. Rich, C.I. Elemental analysis by flame photometery. In Methods of Soil Analysis: Part 2 Chemical and Microbiological Properties, 9.2; American Society of Agronomy: Washington, DC, USA, 1965; Volume 9, pp. 849–865. [Google Scholar]
  35. Margalef, R. Perspectives in Ecological Theory; University of Chicago Press: Chicago, IL, USA, 1968; p. 111. [Google Scholar]
  36. Pielou, E.C. Species-diversity and pattern-diversity in the study of ecological succession. J. Theor. Biol. 1966, 10, 370–383. [Google Scholar] [CrossRef]
  37. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  38. Simpson, E.H. Measurment of diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
  39. Xu, S.; Lui, L.L.; Sayer, E.J. Variablity of aboveground litter inputs alters soil physicochemical and biological process: A meta-analysis of litterfall-manipulation experiments. Biogeosciences 2013, 10, 7423–7433. [Google Scholar] [CrossRef] [Green Version]
  40. Selassie, Y.G.; Ayanna, G. Effect of different land use systems on selected physicochemical properties of soils in Northwestern Ethiopia. J. Agric. Sci. 2013, 5, 112. [Google Scholar]
  41. Zhang, W.; Ren, C.; Deng, J.; Zhao, F.; Yang, G.; Han, X.; Tong, X.; Feng, Y. Plant functional composition and species diversity affect soil C, N and P during secondary succession of abandoned farmland on loess plateau. Ecol. Eng. 2008, 122, 91–99. [Google Scholar] [CrossRef]
  42. Ahmed, S.; Marimuthu, N.; Tripathy, B.; Julka, J.M.; Chandra, K. Earthworm community structure and diversity in different land use systems along an elevation gradient in the Western-Himalaya, India. Appl. Soil Ecol. 2022, 176, 104468. [Google Scholar] [CrossRef]
  43. Lavelle, P.; Blanchart, E.; Martin, A.; Spain, A.V.; Martin, S. Impact of Soil fauna on the properties of soils in the humid tropics. Myth. Sci. Soils Trop. 1992, 29, 157–185. [Google Scholar]
  44. Dachaine, J.; Ruan, H.; Sanchez-de Leon, Y.; Zou, X. Correlation between earthworms and plant litter decompositions in a tropical wet forest of Puerto Rico. Pedobiologia 2005, 49, 601–607. [Google Scholar] [CrossRef]
  45. Bottinelli, N.; Hedde, M.; Jouquet, P.; Capowiez, Y. An explicit definition of earthworm ecological categories-Marcel bouche’s triangle revisited. Geoderma 2020, 372, 114361. [Google Scholar] [CrossRef]
  46. Eggleton, P.; Inward, K.; Smith, J.; Jones, D.T.; Sherlock, E. A six year study of earthworms (Lumbricidae) populations in pasture woodland in southern England shows their responses to soil temperature and soil moisture. Soil Biol. Biochem. 2009, 41, 1857–1865. [Google Scholar] [CrossRef]
  47. Xu, S.; Johnson-Maynard, J.L.; Prather, T.S. Earthworm density and biomass in relation to plan diversity and soil properties in a Palouse prairie remnant. Appl. Soil Ecol. 2013, 72, 119–127. [Google Scholar] [CrossRef]
  48. De Deyn, G.B.; Van der Putten, W.H. Linking aboveground and belowground diversity. Trends Ecol. Evol. 2005, 20, 625–633. [Google Scholar] [CrossRef]
  49. Teng, S.K.; Aziz, N.A.; Anang, N.; Mustafa, M.; Ismael, A.; Yan, Y.W. Earthworms diversity and population density in the Kaki Bukit agroecosystem, Perlis, Peninsular Malasia. Trop. Ecol. 2013, 54, 291–299. [Google Scholar]
Figure 1. Sampling sites of the study in Solan, Himachal Pradesh, India.
Figure 1. Sampling sites of the study in Solan, Himachal Pradesh, India.
Soilsystems 06 00076 g001
Figure 2. Samples of collected earthworms from different land use systems (Pa, SP, and MEF) (Photo: Kasahun Gudeta).
Figure 2. Samples of collected earthworms from different land use systems (Pa, SP, and MEF) (Photo: Kasahun Gudeta).
Soilsystems 06 00076 g002
Figure 3. Declining of soil temperature from Pa to SP to MEF.
Figure 3. Declining of soil temperature from Pa to SP to MEF.
Soilsystems 06 00076 g003
Figure 4. Declining of soil moisture from MEF to SP to Pa.
Figure 4. Declining of soil moisture from MEF to SP to Pa.
Soilsystems 06 00076 g004
Figure 5. (a) Land use−wise carbon concentration (%), (b) land use-wise available potassium concentration (kg/ha), (c) land use-wise available phosphorus concentration (kg/ha), (d) land use-wise available nitrogen concentration (kg/ha); the blue dotes (outliers) in all figures are used show some anomalous results were observed in replications.
Figure 5. (a) Land use−wise carbon concentration (%), (b) land use-wise available potassium concentration (kg/ha), (c) land use-wise available phosphorus concentration (kg/ha), (d) land use-wise available nitrogen concentration (kg/ha); the blue dotes (outliers) in all figures are used show some anomalous results were observed in replications.
Soilsystems 06 00076 g005
Figure 6. (a) Season-wise carbon concentration (%), (b) season-wise available potassium concentration (kg/ha), (c) season-wise available phosphorous concentration (kg/ha), (d) season-wise available nitrogen concentration (kg/ha); the blue dotes (outliers) in all figures are used show some anomalous results were observed in replications.
Figure 6. (a) Season-wise carbon concentration (%), (b) season-wise available potassium concentration (kg/ha), (c) season-wise available phosphorous concentration (kg/ha), (d) season-wise available nitrogen concentration (kg/ha); the blue dotes (outliers) in all figures are used show some anomalous results were observed in replications.
Soilsystems 06 00076 g006aSoilsystems 06 00076 g006b
Table 1. Design of sampling site allocation during each sampling time (2019–2020).
Table 1. Design of sampling site allocation during each sampling time (2019–2020).
Study Site PlotsArea of Each PlotNo. of Monoliths from Each PlotTotal No. of Monoliths (25 × 25 × 30 cm)Containers of Earthworm Samples
Pa310 m−251515
SP310 m−251515
MEF310 m−251515
Total9 154545
Table 2. Sampling time (seasons and months) February 2019 to November 2020.
Table 2. Sampling time (seasons and months) February 2019 to November 2020.
SeasonsMonths
WinterDecember–February
SpringMarch–April
SummerMay–June
MonsoonJuly–September
AutumnOctober–November
Table 3. Overall seasonal variation of soil physicochemical parameters in different land use types.
Table 3. Overall seasonal variation of soil physicochemical parameters in different land use types.
Soil ParameterDifferent Land Use Systems
Pa Mean ± SDSP Mean ± SDMEF Mean ± SD
Sand22.6 ± 2.2723.3 ± 2.4920.3 ± 2.21
Silt43.4 ± 3.2043.2 ± 1.9851.2 ± 2.04
Clay34 ± 3.7133.5 ± 3.3028.5 ± 3.40
Temperature (°C)19.91 ± 2.2718.17 ± 1.7216.87 ± 1.91
Moisture (%)21.08 ± 3.5925.47 ± 4.1829.47 ± 5.69
pH5.82 ± 0.126.08 ± 0.126.14 ± 0.12
Organic Carbon (%) (C)1.378 ± 0.091.964 ± 0.122.132 ± 0.09
Available Nitrogen (N) (kg/ha)252.53 ± 6.80272.47 ± 1.64272.83 ± 1.43
Available Phosphorus (P) (kg/ha)12.36 ± 0.2414.84 ± 0.1815.79 ± 0.14
Available Potassium (K) (kg/ha)89.62 ± 9.75150.21 ± 4.16143.25 ± 7.81
Pa = pasture, SP = silvopasture, MEF = mixed evergreen forest, SD = standard deviation.
Table 4. Correlation value of temperature with moisture.
Table 4. Correlation value of temperature with moisture.
SeasonsR
Winter−0.4820
Spring−0.8156
Summer0.0375
Monsoon0.1164
Autumn0.2946
Table 5. The two-way ANOVA result of soil nutrients (C, K, P, and N).
Table 5. The two-way ANOVA result of soil nutrients (C, K, P, and N).
NutrientsSourcePartial SSDfMSFp-Value
CModel14.167682121.18064027.830.0000
Season-wise3.016540.7541255.000.0012
Land use-wise11.15118281.39389789.250.0000
KModel118,784.11129898.67627.560.0000
Season-wise19,008.1944752.04743.630.0092
Land use-wise99,775.925812,471.9919.520.0000
PModel204.797261217.06643823.940.0000
Season-wise11.57364942.89341224.060.0049
Land use-wise193.22361824.15295133.880.0000
NModel18,030.726121502.56059.590.0000
Season-wise623.008754155.752190.990.4160
Land use-wise17,407.71782175.964713.890.0000
Table 6. Land use post hoc analysis results of soil nutrients (C, K, P, and N).
Table 6. Land use post hoc analysis results of soil nutrients (C, K, P, and N).
Nutrients Pair-Wise Comparison of Land Use SystemsContrastSt. ErrorTukey
tp-Value
CPa1 vs. MEF1 −1.2030.189996−6.330.000
Pa2 vs. MEF1−0.8360.189996−4.400.001
Pa1 vs. MEF2−0.9820.189996−5.170.000
Pa2 vs. MEF2 −0.6150.189996−3.240.044
Pa1 vs. MEF3−0.9380.189996−4.940.000
SP1 vs. Pa1 0.8590.1899964.520.001
SP2 vs. Pa1 0.8180.1899964.310.001
SP3 vs. Pa10.9310.1899964.900.001
KSP3 vs. MEF173.76117.204044.290.002
Pa2 vs. MEF2−56.35617.20404−3.280.039
Pa3 vs. MEF2−80.00117.20404−4.650.000
Pa3 vs. MEF3 −75.33317.20404−4.380.001
SP3 vs. Pa181.34417.204044.730.000
SP3 vs. Pa2 99.91917.204045.810.000
SP3 vs. Pa3 123.56417.204047.180.000
SP3 vs. SP1 74.25517.204044.320.001
SP3 vs. SP2 64.87917.204043.770.009
PPa1 vs. MEF1 −3.8970.4051041−9.620.000
Pa2 vs. MEF1−3.6180.4051041−8.930.000
Pa3 vs. MEF1 −3.3470.4051041−8.260.000
SP1 vs. MEF1 −1.5630.4051041−3.860.007
Pa1 vs. MEF2 −2.9020.4051041−7.160.000
Pa2 vs. MEF2 −2.6230.4051041−6.470.000
Pa3 vs. MEF2 −2.3520.4051041−5.810.000
Pa1 vs. MEF3 −4.0510.4051041−10.000.000
Pa2 vs. MEF3 −3.7720.4051041−9.310.000
Pa3 vs. MEF3−3.5010.4051041−8.640.000
SP1 vs. MEF3 −1.7170.4051041−4.240.002
SP1 vs. Pa1 2.3340.40510415.760.000
SP2 vs. Pa1 2.7390.40510416.760.000
SP3 vs. Pa1 3.1820.40510417.850.000
SP1 vs. Pa2 2.0550.40510415.070.000
SP2 vs. Pa2 2.460.40510416.070.000
SP3 vs. Pa2 2.9030.40510417.170.000
SP1 vs. Pa3 1.7840.40510414.400.001
SP2 vs. Pa3 2.189 0.40510415.400.000
SP3 vs. Pa3 2.6320.40510416.500.000
NPa1 vs. MEF1 −38.795.596977−6.930.000
Pa1 vs. MEF2 −45.7365.596977−8.170.000
Pa2 vs. MEF2 −25.7835.596977−4.610.000
Pa3 vs. MEF2 −21.9375.596977−3.920.006
Pa1 vs. MEF3−37.815.596977−6.760.000
Pa2 vs. Pa119.9535.5969773.560.017
Pa3 vs. Pa1 23.7995.5969774.250.002
SP1 vs. Pa1 32.1265.5969775.740.000
SP2 vs. Pa136.5895.5969776.540.000
SP3 vs. Pa1 46.8095.5969778.360.000
SP3 vs. Pa2 26.8565.5969774.800.000
SP3 vs. Pa3 23.015.5969774.110.003
Table 7. Density (ind·m−2) and biomass (g m−2) of earthworms across different land use types.
Table 7. Density (ind·m−2) and biomass (g m−2) of earthworms across different land use types.
Season and Land Use System 2019–20202019–2020
Mean Density (Ind·m−2) ± SEMean Biomass (g m−2) ± SE
Winter
Pa6.4 ± 0.164.5 ± 0.11
SP87.47 ± 1.0520.63 ± 0.29
MEF64 ± 0.8524.7 ± 0.27
Spring
Pa35.2 ± 0.4712.98 ± 0.24
SP60.8 ± 1.3816.75 ± 0.22
MEF120.54 ± 0.5020.08 ± 0.28
Summer
Pa3.2 ± 0.141.54 ± 0.05
SP37.32 ± 0.4424.85 ± 0.39
MEF157.87 ± 1.7635.75 ± 0.31
Monsoon
Pa22.4 ± 0.463.93 ± 0.11
SP68.27 ± 0.9436.91 ± 0.56
MEF230.4 ± 2.3843.06 ± 0.51
Autumn
Pa10.67 ± 0.418.87 ± 0.21
SP52.27 ± 0.5225.7 ± 0.23
MEF181.34 ± 2.7532.43 ± 0.20
Total1138.15 ind·m−2312.68 g m−2
Table 8. Biomass (g m−2) and density (ind·m−2) of different earthworm species in different land use types.
Table 8. Biomass (g m−2) and density (ind·m−2) of different earthworm species in different land use types.
Earthworm SpeciesPaSPMEF
Ind·m−2
Mean ± SE
g m−2
Mean ± SE
Ind·m−2
Mean ± SE
g m−2
Mean ± SE
Ind·m−2
Mean ± SE
g m−2
Mean ± SE
A. corticis45.87 ± 0.2924.65 ± 0.33217.6 ± 2.59102.67 ± 0.80109.87 ± 1.3737.36 ± 0.90
A. rosea8.53 ± 0.231.66 ± 0.076.4 ± 0.190.92 ± 0.02605.87 ± 7.44101.33 ± 2.18
D. japonica12.8 ± 0.272.03 ± 0.0820.27 ± 0.349.84 ± 0.3920.27 ± 0.5811.69 ± 0.20
M. houlleti1.07 ± 0.040.51 ± 0.0243.73 ± 0.377.83 ± 0.2712.8 ± 0.324.4 ± 0.15
M. birmanica6.4 ± 0.1552.24 ± 0.098.53 ± 0.243.01 ± 0.101.07 ± 0.060.17 ± 0.01
E. fetida--9.6 ± 0.240.57 ± 0.024.27 ± 0.151.07 ± 0.06
L. pusillus3.2 ± 0.090.73 ± 0.02----
Total77.8731.82306.13124.84754.15156.02
Table 9. The density (ind·m−2) and relative density (RD%) of earthworm species in different land use types.
Table 9. The density (ind·m−2) and relative density (RD%) of earthworm species in different land use types.
Earthworm Family/SpeciesEcological CategoryPaSPMEF
Ind·m−2RD%Ind·m−2RD%Ind·m−2RD%
MegascolecidaeA. corticesEpi-endogeic45.8758.91217.671.08109.8714.57
LumbricidaeA. roseaEndogeic8.5310.956.42.09605.8780.34
MoniligasteridaeD. japonicaEndogeic12.816.4420.276.6220.272.69
MegascolecidaeM. houlletiEpi-endogeic1.071.3743.7314.2812.81.69
MegascolecidaeM. birmanicaEndogeic6.48.228.532.791.070.14
LumbricidaeE. fetidaEpigeic--9.63.144.270.57
OctochaetidaeL. pusillusEndogeic3.24.11- -
Total 77.87100306.13100754.15100
Table 10. Diversity index of earthworms in the three land use systems.
Table 10. Diversity index of earthworms in the three land use systems.
Earthworm Species and Parameter of Alpha DiversityDifferent Land Use Systems
PaSPMEF
A. cortices43204103
A. rosea86568
D. japonica121919
M. houlleti14112
M. birmanica681
E. fetida-94
L. pusillus3--
Total number of individuals 73287707
Shannon index of diversity (H′)1.250.990.661
Index of evenness (J′)0.6960.5520.369
Margalef species richness index (R) 1.1650.8830.762
Index of dominance (D)0.390.530.67
Species richness (n)666
Table 11. The relationship between diversity, density, and biomass of earthworms.
Table 11. The relationship between diversity, density, and biomass of earthworms.
Population ParameterLand Use System
PaSPMEF
Density (ind·m−2)77.87306.13754.15
Biomass (g m−2)31.82124.84156.02
Diversity parameter
Shannon index of diversity (H′)1.250.990.661
Index of evenness (J′)0.6960.5520.369
Margalef species richness index (R)1.165 0.8830.762
Index of dominance (D)0.390.530.67
Species richness (n)666
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gudeta, K.; Bhagat, A.; Julka, J.M.; Bhat, S.A.; Sharma, G.K.; Bantihun, G.; Amarowicz, R.; Belina, M. Impact of Aboveground Vegetation on Abundance, Diversity, and Biomass of Earthworms in Selected Land Use Systems as a Model of Synchrony between Aboveground and Belowground Habitats in Mid-Himalaya, India. Soil Syst. 2022, 6, 76. https://doi.org/10.3390/soilsystems6040076

AMA Style

Gudeta K, Bhagat A, Julka JM, Bhat SA, Sharma GK, Bantihun G, Amarowicz R, Belina M. Impact of Aboveground Vegetation on Abundance, Diversity, and Biomass of Earthworms in Selected Land Use Systems as a Model of Synchrony between Aboveground and Belowground Habitats in Mid-Himalaya, India. Soil Systems. 2022; 6(4):76. https://doi.org/10.3390/soilsystems6040076

Chicago/Turabian Style

Gudeta, Kasahun, Ankeet Bhagat, Jatinder Mohan Julka, Sartaj Ahmad Bhat, Gopal Krishan Sharma, Getachew Bantihun, Ryszard Amarowicz, and Merga Belina. 2022. "Impact of Aboveground Vegetation on Abundance, Diversity, and Biomass of Earthworms in Selected Land Use Systems as a Model of Synchrony between Aboveground and Belowground Habitats in Mid-Himalaya, India" Soil Systems 6, no. 4: 76. https://doi.org/10.3390/soilsystems6040076

APA Style

Gudeta, K., Bhagat, A., Julka, J. M., Bhat, S. A., Sharma, G. K., Bantihun, G., Amarowicz, R., & Belina, M. (2022). Impact of Aboveground Vegetation on Abundance, Diversity, and Biomass of Earthworms in Selected Land Use Systems as a Model of Synchrony between Aboveground and Belowground Habitats in Mid-Himalaya, India. Soil Systems, 6(4), 76. https://doi.org/10.3390/soilsystems6040076

Article Metrics

Back to TopTop