Submit manuscript...
eISSN: 2577-8307

Forestry Research and Engineering: International Journal

Research Article Volume 2 Issue 2

Community structure and recolonization by earthworms in rehabilitated ecosystems in garhwal himalayas, India

Bhadauria T,1 Saxena KG2

1Department of Zoology, Feroze Gandhi P.G. Degree College, India
2School of Environmental Sciences, Jawaharlal Nehru University, India

Correspondence: Tuneera Bhadauria, Department of Zoology, Feroze Gandhi College, Raebareli, Uttar Pradesh, India, Tel 91 9453022105

Received: November 22, 2017 | Published: March 7, 2018

Citation: Bhadauria T, Saxena KG. Community structure and recolonization by earthworms in rehabilitated ecosystems in garhwal himalayas, India. Forest Res Eng Int J. 2018;2(2):35-47. DOI: 10.15406/freij.2018.02.00024

Download PDF

Abstract

Consequences of deforestation arise from site degradation leading to modification of soil properties this significantly affect both incidence and abundance of soil macro fauna. Earthworm communities are more directly altered by these changes. Endemic and exotic species co-existed in the study area following deforestation and intensive cultivation. That native species were dominant in the undisturbed sites and disturbance and degradations leads to invasion by the exotic species holds true in our study. The sites under study represented the degraded areas as none of the species reported from the present experimental plots were endemic to the region, all the species are either peregrine exotic or peregrine endemic to the area as many of the endemic species of this region probably exterminated during the last Quaternary Glaciation. The numbers of species present at our sites ranged between 2 to 5. The presence of litter layer and lower perturbation pressure resulted in the numerical dominance of Lennogaster pussilus in the Oak forest (OF) and the higher biomass of A. alexandri may be because of the larger size of the earthworm. B parvus, L. pussilus and P. excavatus are litter-associated taxa which were more directly affected by OF clearance and the resulting decrease in available litter, thus explaining their disappearance in the changed ecosystems, however improved soil moisture and temperature as well as input of organic matter in rehabilitated agricultural land (RAL) could probably be favourable factor for decolonization anddominance of L. pussilus. A. alexandri had wider ecological amplitudeoccurring under all land use types as reported in our studies, conversion of two other land use types resulted in lower species diversity. Our results show that the abandoned agriculture land (AAL) remains closer to the AL(T) than to the forest because in these land use types the overall vegetation diversity remains low corresponding to a low diversification of the organic resources thus explaining the similarity between the earthworm communities in land use types. Seasonal rhythmic pattern was exhibited by all the species identified.

Keywords: deforestation, rehabilitation, earthworm community, functional guild, ecological categories, diversity index

Introduction

Himalaya is a vast and diverse mountain system. Agro forestry land use covering (20%) of geographical area of Indian Himalayas is distributed as patches in the matrix of forest covering (52%) area. Forest based tree- crop livestock integrated farming is the predominant traditional land use in central Himalayas.1 Hill agriculture appears to be key threat to soil biodiversity and ecosystem service due to huge amount of biomass extraction to sustain live stock and produce manure for managing soil fertility and also in terms of direct loss of forest cover due to agriculture land use.2 This problem has been further accentuated by extensive deforestation and unsustainable land use causing extensive degradation of the mountains. A pilot project was carried out to restore and rehabilitate the abandoned agriculture land site and degraded forest land site through agro forestry and forestry plantations.3 Land degradation implies a decline in the useable natural resource base and therefore loss of biodiversity and ecosystem functions. The interacting functions of soil organisms and the effects of human activities in managing land for agriculture and forestry affect soil health and quality.4 Soil faunal biodiversity is an important resource for environmental monitoring and natural resource management, changes in the variety and abundance of organisms in response to ecosystem disturbance, degradation and rehabilitation provide important management information.5 Soil management options can have dramatic effects upon soil invertebrate communities.6,7 Earthworm communities are the result of both interactions between species8 and sensitivity to ecological factors9 presence or absence of ground vegetation and changes in its composition are known to affect the composition of earthworm communities.10 This is the first systematic study where work is focused primarily on the ecological impact, the extent, causes and consequences of varying land degradation and subsequent rehabilitation strategies on earthworm fauna in Central Himalayas. Micro scale variability in earthworm community in traditional agro ecosystem will also be focused by comparing two types of micro sites viz below tree canopy from (1-3m) the bole and outside tree canopy (1-3m) from canopy margin. Changes in the vegetation cover (forest clearance, tree plantation) change the soil microclimate condition therefore present study also aims to analyze how these changes affected the functional guild of the earthworms.

Soil fauna vary through time as they have seasonal rhythms mainly regulated by temperature and moisture and thus constitute one of the important factors of changes in the species assemblage structure.11,12 Consequently, the way seasonal variations of earthworm communities impacted by land-use changes occurs as well as its impact on biomass and density of earthworms have been examined across the sampling sites.

Materials and methods

  1. description of study site
  2. The study was carried out at Banswara village located at 1200m above mean sea level in Chamoli district, Garhwal (latitude 30°27’N and 79°5’E). The climate is typical monsoon, monthly minimum and maximum temperature varying in the range of 6-21°C and 18-35 °C, respectively, with an average annual rainfall of 1700mm. 80% of total rain fall is received during the monsoon period of July to September (Figure 1A) (Figure 1B). Soil moisture and temperature are as presented in Figure 2(A) & Figure 2(B). The soil is sandy loam to loamy sand in texture (Table 1) derived from felspathic quartz schists ,quartz muscovite schists and quartz chlorite schist.13 To study the impact of changed land use practices on earthworms following land uses were identified in the village Banswara.

    Figure 1(A) Monthly maximum and minimum atmospheric temperature for study area in Garhwal Himalayas during the study period from July 2008-2009.

    Figure 1(B) Monthly Rainfall for study area in Garhwal Himalayas during the study period from July 2008-May 2009.

    Figure 2(A) Soil moisture (%) under different land use.

    Figure 2(B)Soil temparature (°C) under different land use.

  3. Old secondary forest (OF)
  4. The reserved forests are unmanaged old-growth sub tropical warm temperate forest dominated by Quercus samicarpifolia, Q.floribunda Cedrus .

  5. Traditional rainfed farming systems
  6. The farming system is characterized by settled rainfed organic agriculture on terraced slope with scattered multipurpose trees

  7. Microhabitat below agro forestry tree patches AL(T)
  8.  The sampling of soil fauna was done below was done below the tree canopy (1-3mdiameter from the bole).

  9. Microhabitat without agro forestry tree patches AL(WT)
  10. The sampling of soil fauna was done below was done below the tree canopy (1-3m from canopy margin).

  11. Abandoned agriculture land (AAL)
  12. Degraded abandoned agriculture land represented the agriculture terraces which had been abandoned due to poor soil fertility status, they were degraded and exposed to uncontrolled grazing.

  13. Rehabilitated agriculture land (RAL)
  14. The rehabilitated land site was the abandoned agriculture land site where the terraces were repaired and agricultural crops were grown providing supplemental irrigation and organic manure along with the planted trees following traditional farming practices. This site was protected from grazing and other human disturbances.

  15. Degraded forest land (DFL)
  16. Degraded Forest Land represented degraded community forest land which had negligible tree cover and was exposed to uncontrolled grazing. Soils on the degraded land were considerably eroded with respect to abandoned agriculture land.

  17. Rehabilitated forest land (RFL)
  18. The Rehabilitated forest site represented the plantation of trees on degraded community forest land which had more stressful environment as compared to the RAL. This site was protected from grazing and other human disturbances.

  19. Experimental design
  1. Soil sampling: In each land use type a plot of 40 × 50 m 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake GabaagyKqzGeaeaaaaaaaaa8qacaaI0aGaaGimaiaabccacqGHxdaT caqGGaGaaGynaiaaicdacaWGTbWcpaWaaWbaaKqaGeqabaqcLbmape GaaGOmaaaaaaa@4341@ was demarcated for earthworms and for soil sampling, three composite soil samples were prepared in each experimental plot. Six 25cm × 25cm × 30cm MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaaGOmaiaaiwdacaWGJbGaamyBaiaabcca cqGHxdaTcaqGGaGaaGOmaiaaiwdacaWGJbGaamyBaiaabccacqGHxd aTcaqGGaGaaG4maiaaicdacaWGJbGaamyBaaaa@4A23@ deep soil monoliths were randomly sampled from each replicate plot at regular bimonthly intervals. Each monolith was subdivided into 0-10, 10-20, and 20-30cm blocks, the soil samples were air dried and sieved through a 2mm sieve. A representative of the sub-sample was stored for subsequent analysis.14 Soil temperatures were recorded weekly at 0-10 and 10-20cm depths; however, the values presented here are mean monthly values for 0-10 cm depth. Soil moisture was recorded every month at 0-10 and 10-20cm depths, but the values presented here are mean monthly values for 0-10cm depth and are expressed as % oven dry weight at 105 °C. Bulk density was estimated following standard methods.14 The analysis of soil texture was done using a hydrometer method,15 soil pH was measured as 1: 2.5 (soil: water) solution, and organic C through the.16 Soil N was analyzed using the semi-micro Kjeldahl method following the procedures described.17
  2. Soil fauna sampling: Earthworms were sampled using the soil biology and fertility methodology (17). On each site in a plot of 40m × 50m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaaGinaiaaicdacaWGTbGaaeiiaiabgEna 0kaabccacaaI1aGaaGimaiaad2gaaaa@41A0@ earthworms were collected at regular bimonthly intervals over a period of 12 months, i.e. between June 2007 and June 2008. Earthworms were collected from 10 sampling points 5 m apart along a transect with a random origin, they were extracted by hand sorting after digging a trench up to 30 cm deep around a 25cm × 25cm MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaaGOmaiaaiwdacaWGJbGaamyBaiaabcca cqGHxdaTcaqGGaGaaGOmaiaaiwdacaWGJbGaamyBaaaa@4375@ area at each sampling point to get a soil monolith. These soil monoliths were divided into three layers (0-10, 10-20, 20-30cm) and earthworms were extracted from each layer, they were then washed and preserved in 5% formalin for further identification.17
  3. Statistical analysis of the data: Statistical analysis was done following the biostatistical methods described in (18).
  1. Significant differences in physicochemical characteristics of the soil across different sampling sites were carried out using one-way ANOVA.
  2. Significant differences in the density and biomass of earthworm species across different sampling sites (interhabitat variations) were tested using non-parametric one-way ANOVA (F-test), and the New Mann Keul’s multiple range test.
  3. Significant differences in functional guild changes of earthworm’s species in the same site and functional guild changes of earthworm species across different sampling sites were were tested using non-parametric one-way ANOVA (F-test), and the New Mann Keul’s multiple range test.
  4. Seasonal variations in total density and biomass of earthworm species in sampling sites were tested using the non-parametric Kruskal-Wallis test of variance and New Mann Keul’s multiple range test.
  5. Diversity index was calculated as Simpsons Index of Diversity
D=1 i=1 n ni( ni1 )/N( N1 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiraiabg2da9iaaigdacqGHsisljuaG daGfWbGcbeqcbaYdaeaajugWa8qacaWGPbGaeyypa0JaaGymaaqcba YdaeaajugWa8qacaWGUbaan8aabaqcLbsapeGaeyyeIuoaaiaad6ga caWGPbqcfa4aaeWaaOWdaeaajugib8qacaWGUbGaamyAaiabgkHiTi aaigdaaOGaayjkaiaawMcaaKqzGeGaai4laiaad6eajuaGdaqadaGc paqaaKqzGeWdbiaad6eacqGHsislcaaIXaaakiaawIcacaGLPaaaaa a@56A7@

Where ni MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamOBaiaadMgaaaa@3B50@ = number of individuals of species i in the sample. N=total number of individuals in the sample.

S= Number of species in the sample. 1-D^ = Simpsons Index of Diversity

2 Xij Xjk  { ( X 2 ij N 2 j )+( X 2 jk/  N 2 k )  NjNk } MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajuaGqaaaaaaaaaWdbmaalaaak8aabaqcLbmapeGaaGOmaKqba+aa daqfGaGcbeWcbeqaaKqzGeGaaGzaVdqdbaqcLbsapeGaeyyeIuoaaK qzadGaaeiwaiaabMgacaqGQbGaaeiOaiaabIfacaqGQbGaae4Aaiaa bckaaOWdaeaajuaGpeWaaiWaaOWdaeaajuaGpeWaaeWaaOWdaeaaju aGpeWaaSaaaOWdaeaajuaGdaqfGaGcbeWcbeqaaKqzGeGaaGzaVdqd baqcLbsapeGaeyyeIuoaaKqzadGaamiwaSWdamaaCaaajeaqbeqaaK qzadWdbiaaikdaaaGaaeyAaiaabQgaaOWdaeaajugib8qacaWGobqc fa4damaaCaaaleqabaqcLbmapeGaaGOmaaaajugibiaabQgaaaaaki aawIcacaGLPaaajugibiabgUcaRKqbaoaabmaak8aabaqcfa4aaubi aOqabSqabeaajugibiaaygW7a0qaaKqzGeWdbiabggHiLdaajugWai aabIfajuaGpaWaaWbaaKqaGeqabaqcLbmapeGaaGOmaaaacaqGQbGa ae4Aaiaac+cacaqGGcGaamOtaKqba+aadaahaaqcbasabeaajugWa8 qacaaIYaaaaiaabUgaaOGaayjkaiaawMcaaKqzGeGaaeiOaiaabcka caqGobGaaeOAaiaab6eacaqGRbaakiaawUhacaGL9baaaaaaaa@7DE4@
  1. The dominance value was calculated as, Dominance=Total number of individuals of a species/ Total number of individuals of all the species
  2. A simplified Morisita’s Index of Similarity was calculated as

Where

Xij is the number of individuals of species i in group j

Xik is the number of individuals of species i in group k

N is the total number of individuals in the sample j and k

  1. The correlation between soil parameters such as temperature, moisture and organic matter and earthworm species was calculated as a simple correlation coefficient (r). Sample standard error was calculated as the standard error of the Mean (±S E).
  2. The SYSTAT 12 software for Windows package (Systat Software Inc., San Jose, USA) was used for all statistical tests

Result

  1. Soil characteristics
  2. Soil characteristics varied under different land use types, except for RFL soil pH was mildly acidic in all other sites, soil bulk density was higher in AAL and DFL when compared to other land use types under study. Soil organic C (F0.05,6,14= 32.9) and Nitrogen (F0.05,6,14= 151) varied significantly between different sites, organic C was significantly higher (F0.05,6,14= 151) in AL (T) as compared to all other sites but did not vary significantly between RFL and OF, between AL (WT) and OF. Total nitrogen was significantly higher (F0.05,6,14= 151) in AL (WT) as compared to all other sites, but it did not vary significantly between AL (WT) and RFL and between AL (T) and OF. Soil organic carbon (q0.05,21,3 = 3.41) and nitrogen (q0.05, 21,3= 4.76) was significantly higher in AAL as compared to DFL (Table1).

  3. Selected attributes of sampled earthworm species across different land use types
  4. A total of eight species belonging to four families were recorded from land use types under study and the selected attributes of sampled earthworm species (Table 2). All the earthworm species recorded across different land use types were both peregrine endemic or peregrine exotic to the region, none of these were endemic to Garhwal Himalayas (Table 2).

  5. Effect of land use patterns on earthworm species richness
  6. Earthworm species richness did not vary significantly between the OF, AL (T) and AL (WT), but in AAL and DFL number of species was lower; however rehabilitation of these ecosystems led to increase in species number in RAL and RFL (Table 3).

  7. Changes in earthworm species composition across different land use types
  8. Three peregrine exotics Amynthas alexandri, Bimastos parvus, Metaphire anamola and two peregrine native species Lennogaster pussilus, Perionyx excavatus were present in the OF .Conversion of the forest to agro ecosystems led to the change in the community structure with the loss of exotics Bimastos parvus and natives Lennogaster pussilus, Perionyx excavatus in AL (T) and AL (WT), but recolonisation by exotic Metaphire birmanica, Octochaetona beatrix and natives Drawida nepalensis occurred here. M. anamola, M. birmanica and D nepalensis were present in AAL, however rehabilitation of degraded ecosystems RAL resulted in loss of M. birmanica and recolonisation by L. pussilus and A. alexandri. Only exotic A. alexandri and M. birmanica were present in DFL, but A. alexandri, M. anamola and M. birmanica .were present in RFL (Figure 3). B. parvus and P. excavatus did not recolonize any rehabilitated ecosystems whereas A. alexandri and M. anamola were absent in RAL and DFL though they were present on all the other sites under study (Table 4).

    Figure 3 Changes in earthworm species composition across different land use types in Garhwal Himalayas.

  9. Total density and biomass of earthworm species across different land use types
  10. Total density (F0.05,6,14= 228.24) as well as biomass (F0.05,6,14= 403.78) of earthworm’s varied significantly between different land use types due to changes in land use practices. Conversion of forest to agro ecosystem resulted in significantly higher density (q,0.05,,14,7= 122.73) and biomass (q,0.05,14,7= 109 values ) in AL(T). Density (q,0.05,14,4= 114) and biomass of earthworms (q,0.05,14,4= 112) also increased significantly in rehabilitated ecosystems RAL and RFL when compared to AAL and DFL (Table 5).

  11. Endemic and exotics earthworm species
  12. The density and biomass of exotic peregrine earthworm species varied significantly between different sites (Figure 4A) (Figure 4B). These values were significantly higher in AL (T) followed by AL (WT) (q0.05,12,6= 102) and did not vary significantly between RAL, DFL and AAL. The biomass values of exotic species also did not vary significantly between OF and RFL. The density (F0.05,6,14= 211) and biomass (F0.05,6,14= 248) of endemic peregrine species also varied significantly between different sites and were significantly higher in AL(T) as compared to other sites. The biomass values of endemic species did not vary significantly between OF and AAL.

    Figure 4(A) Abundance of exotics/native earthworm species under different landuse changes.

    Figure 4(B) Biomass of exotics/native earthworm species under different land use changes.

    Within the same sites the density of exotic and endemic peregrine earthworm species did not vary significantly in OF, but exotic species had higher biomass (t(2)0.05,4)=7.5 as compared to endemics species. Exotic earthworm species were significantly more abundant (t(2)0.05,4= 57.32) with higher biomass (t(2)0.05,4= 25.5) values in AL(T) than endemic species. In AL (WT) and in AAL endemic species were more abundant numerically and had significantly higher biomass value than exotics species. Endemic species had significantly higher abundance in RAL (t(2)0.05,4= 5.75) as compared to exotic species but biomass values of exotic and endemic earthworm species did not vary significantly in here. Endemic species were absent in DFL and RFL.

  13. Functional guild changes in earthworm communities
  14. Functional guild diversity varied under different land use types (Figure 5A) (Figure 5B), all the three functional categories were present in RAL. Anecics were absent in OF and AAL, in AL (T), AL (WT), DFL and RFL epigeic species were absent, where as in AAL only endogeics species were present. Epigeic were significantly more abundant (F0.05,2,6= 12.23 )with higher biomass values (F0.05,2,6= 7.54 ) in OF. Both endogeics and anecics earthworm species were significantly more abundant (F0.05,2,6= 58.54) and had higher biomass values (F0.05,2,6= 42.63 ) at AL(T) compared to all other sites.

    Figure 5(A) Functional categories of earthworms under different land use types.

    Figure 5(B) Functional categories of earthworms under different land use types.

    Within the same site, numerically epigeics formed the dominant (F0.05,2,6= 37.19) functional group under OF and in RAL (F0.05,2,6= 18.54), whereas endogeic were dominant under AL(T) (F0.05,2,6= 11.7), AL(WT) (F0.05,2,6= 28) and RFL(F0.05,2,6= 7.54), in AAL they were the only functional group present. Endogeics had significantly higher (F0.05,2,6= 19.4) biomass in OF and RAL when compared to epigeic and anecics. Anecics had lower abundance and biomass compared to other functional groups within the same sites.

  15. Seasonal variation for earthworm communities
  16. Seasonal variation for different earthworm species abundance and biomass (Figure 6 (A-H)).

    • Figure (A-H) Seasonal variation of earthworms species abundanceunder different land use types.

    Alexandri All the species showed significant seasonal variation. The density and biomass varied significantly seasonally in OF (F0.05,5,12= 124.6), AL(T) (F0.05,5,12= 133.8), AL(WT) (F0.05,5,12= 23.04) and RFL (F0.05,5,12= 8). It did not vary significantly in RAL and DFL. The species had significantly higher population density (q0.05,12,6= 7.54) and biomass (q0.05,12,6= 8) during Sep at all sites except in AL(WT)and where it had higher biomass (q0.05,12,6= 7.34) during July and in RFL where species showed a minor increase in population density during January. B. parvus was present only in the oak forest and did not show any significant variation seasonally in density as well as biomass values.

    1. D. nepalensis: Density and biomass varied significantly seasonally in AL(T) (F0.05,5,12= 83.83), AL(WT) (F0.05,5,12= 437.33), AAL (F0.05,5,12= 16.6) and in RAL (F0.05,5,12= 7.53). The species showed higher population density (q0.05,12,6= 10.56) and biomass (q0.05,12,6= 9.67) in Sep in AL(T)(q0.05,12,6,= 5.87), AL(WT) (q0.05,12,6= 7.22) and AAL (q0.05,12,6= 7.65), however it also showed a minor peak of growth during March in AL(T). In RAL it had higher density (q0.05,12,6= 12.78) and biomass (q0.05,12,6= 6.58) in July and a minor increase in population was observed in November, though the difference in abundance of the species did not vary significantly between Sep and November.
    2. M. anamola: Population abundance varied significantly seasonally in AL(T) (F0.05,5,12= 130.5) and AL(WT) (F0.05,5,12=14.93). At these sites population improved significantly during the (q0.05,12,6= 9.33) September, however the biomass of these species was higher during November in AL(T) (q0.05,12,6= 6.65), AL(WT) (q0.05,12,6= 7.43). Population abundance of M.anamola did not show significant variation seasonally in AAL, RAL and RFL but the biomass of these species was higher during September in AAL (q0.05,12,6= 5.76), RAL (q0.05,12,6= 6.22) and RFL(q0.05,12,6= 5.49) and declined subsequently.
    3. M.birmanica: Population abundance (F0.05,5,12= 62.32)and biomass (F0.05,5,12= 84.45) of M.birmanica varied significantly seasonally and showed increased population density (q0.05,12,6= 9.02) during September in AL(T) with decline subsequently. Population abundance did not show seasonal variation in other sites. It showed higher biomass during Sep in AL(T) (q0.05,12,6= 11.23), AL(WT) (q0.05,12,6= 8.65) and RAL (q0.05,12,6= 10.56), but in RFL it showed higher biomass during November (q0.05,12,6= 12.67).
    4. L.pussilus: Population abundance varied significantly (F0.05,5,12= 206) seasonally in OF and RAL and it had higher density (q0.05,12,6= 112.67) and biomass (q0.05,12,6= 18.43) during September in both these sites.
    5. O.beatrix: Population abundance varied significantly (F0.05,5,12= 346) seasonally in AL(WT) and it had higher density(q0.05,12,6= 13.58) and biomass (q0.05,12,6= 15.43)values during September.
    6. P.excavatus: Population abundance of P. excavatus also varied significantly (F0.05,5,12=58.40) seasonally and had higher density(q0.05,12,6=14.54) and biomass (q0.05,12,6=16.55) during the September in OF.
  17. Diversity index for earthworm communities between different land use types
  18. The Simpsons Diversity Index within the earthworm communities under different land use types was lowest in AAL and maximum in OF. The diversity index was similar between OF and RAL. Within the agro ecosystem the diversity Index was lower in AL (WT) compared to AL (T), but rehabilitation of degraded ecosystem led to improved diversity Index in RAL and RFL (Figure 7).

    Figure 7 Earthworm species diversity index under different land use changes.

  19. Earthworm species dominance (%) (a) abundance (b) biomass (values in parenthesis) within the same land use types
  20. L. pussilus was dominant numerically in OF but A. alexandri had higher biomass values as compared to other species. At AL(T), AL(WT) and AAL D. nepalensis was dominant numerically as well having a higher biomass percentage, L. pussilus was dominant numerically at RAL but D. nepalensis had higher biomass here. A. alexandri was dominant numerically as well as higher biomass values in DFL and RFL but M. anamola had higher biomass percentage at the RFL site (Table 6).

  21. Morisita’s index of similarity for earthworm communities between different land use types
  22. Morrisitas index of similarity indicated that the earthworm community in AL (T) was more similar to AAL than AL (WT). Rehabilitation of degraded agro ecosystem showed that the community structure of RAL was more similar to the OF where as community structure between RFL and DFL were closer to one another (Table 7).

  23. Correlation coefficient for earthworm communities between different land use types

A significant correlation coefficient was observed between earthworm population and soil moisture and temperature. A. alexandri was positively correlated to soil moisture ( P<0.01 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaigdaaaa@3E28@ ) and temperature ( P<0.05 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaiwdaaaa@3E2C@ ). In AL(WT), RAL, DFL. B. parvus, and L. pussilus showed positive correlation ( P<0.05 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaiwdaaaa@3E2C@ ) to soil moisture and temperature, P. excavatus was positively correlated to moisture ( P<0.01 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaigdaaaa@3E28@ ) but did not show any co relation to soil temperature. D. nepalensis was positively related ( P<0.05 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaiwdaaaa@3E2C@ ) to soil moisture in AL (T), AL (WT) and RAL but it did not show any correlation with soil moisture in AAL. M anamola was positively related to soil moisture ( P<0.01 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaigdaaaa@3E28@ ) and temperature ( P<0.05 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaiwdaaaa@3E2C@ ) in AL (WT) AAL and RAL but in AAL it did not show any correlation to soil moisture. In AL (WT) and AAL M. birmanica exibited similar trend as M. anamola however in DFL ( P<0.01 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaigdaaaa@3E28@ ) and RFL ( P<0.05 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaiwdaaaa@3E2C@ ) it showed positive correlation to soil moisture only. O. beatrix showed positive correlation to soil temperature ( P<0.01 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaigdaaaa@3E28@ ) in AL (T), AL (WT), but it did not show any relation to soil moisture (Table 8).

Land use/cover

Soil property

pH

Bulk density

Organic carbon (g kg-1)

Total nitrogen (g kg-1)

(g cm-3)

Oak forest (OF)

6.1a

1.11a

12.6a

1.2a

Microhabitat below agro forestry tree patches AL(T)

6.4a

1.04a

18.1b

0.9a

Microhabitat without agro forestry tree patches (ALWT)

6.2a

1.10a

13.2a

1.9b

Abandoned agricultural land(AAL)

6.4a

1.4b

8.6c

0.8c

Rehabilitated agricultural land (RAL)

6.3a

1.12a

15.0d

1.4d

Degraded forest land(DFL)

6.2a

1.42b

8.3c

0.7c

Rehabilitated forest land (RFL)

5.8b

1.14a

12.4a

1.7b

Table 1 Soil characteristics of different land use/cover types in central Himalayan ( ±SE MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaeyySaeRaam4uaiaadweaaaa@3CFE@ ), India (Average mean values for one year. Values followed by different superscript letters are significantly ( P<0.05 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaiwdaaaa@3E2C@ ) different.

Species

Color

Length*(mm)

Diameter(mm)*

Family

Peregrine Species

Native region

Ecological Catagories

Amynthas alexandri

Dark Brown

180

8

Megascolecidae

Exotic

S.E.Asia

Endogeic

Bimastos parvus

Dark brown

70

2

Lumbricidae

Exotic

N.America

Epigeic

Drawida nepalensis

Light whitish

140

7

Moniligastridae

Native peregrine

Eastern Himalayas

Endogeic/Anecic

Lennogaster pussilus

Light brown

50

2

Octochaetidae

Native peregrine

Middle and Upper Gangetic plains

Epigeic

Metaphire anamola

Dark brown

140

4

Megascolecidae

Exotic

S.E.Asia

Endogeic

Metaphire birmanica

Blackish brown

160

4

Megascolecidae

Exotic

Burma

Endogeic

Octochaetona Beatrix

Light Pink

120

4

Octochaetidae

Native peregrine

Eastern and Western Plateau region

Endogeic

Perionyx excavatus

dark violet /Purple

76

3

Megascolecidae

Native peregrine

Eastern Himalayas

Epigeic

Table 2 Selected attributes of sampled earthworm species across different land use types in Garhwal Himalayas
Measurements of largest individuals indicating the maximum achievable size

Land Use Pattern

Species Number

Old secondary Forest OF

5

Microhabitat below agro forestry tree patches AL(T)

5

Microhabitat without agro forestry tree patches AL(WT)

5

Abandoned Agriculture LandAAL

3

Rehabilitated Agriculture LandRAL

4

Degraded Forest Land DFL

2

Rehabilitated Forest Land RFL

3

Table 3 Effect of land use patterns on earthworm species richness across different land use types in Garhwal Himalayas

Species

OF

ALT

ALWT

AAL

DFL

RAL

RFL

Amynthas alexandri

Y

Y

Y

N

Y

Y

Y

Bimastos parvus

Y

N

N

N

N

N

N

Drawida nepalensis

N

Y

Y

Y

N

Y

N

Lennogaster pussilus

Y

N

N

N

N

Y

N

Metaphire anamola

Y

Y

Y

Y

N

Y

Y

Metaphire birmanica

N

Y

Y

Y

Y

N

Y

Octochaetona Beatrix

N

Y

Y

N

N

N

N

Perionyx excavatus

Y

N

N

N

N

N

N

Table 4 Occurrence of earthworm species across different land use types in Garhwal Himalayas

Density(individuals/m2/year)

Biomass(g/m2/year)

OF

a97±18

a 79±2.7

AL(T)

b358±25

b 663±35.74

AL(WT)

c 263±21

c 403±16±37

AAL

d 19±1.4

d 27±2.4

RAL

e 59±1.8

e 50±3.2

DFL

f 7±3.4 7

f 15.4±1.2

RFL

d 28±1.1

g 62.5±4.1

Table 5 Density (Individuals/m2/year) and biomass (g/m2/year) of earthworm species across different land use types in Garhwal Himalayas. Numbers followed by the same letter are not significantly different ( P<0.05 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaiwdaaaa@3E2C@ ).

OF

AL(T)

AL(WT)

AAL

RAL

DFL

RFL

A. alexandri

16.5

11.84

4.57

0

12

*57.15

*36

(*39)

(12.7)

(6)

(19)

(*57)

(29)

B.parvus

7.22

0

0

0

0

0

0

(3.1)

(0)

(0)

(0)

(0)

(0)

(0)

D. nepalensis

0

*40

*44.87

*61

21

0

0

0

(*38)

(*53)

(*56)

(*41)

0

0

M.anamola

26.8

16.9

5.71

12

8

0

32.14

(32)

(16)

(6)

(24)

(16)

(0)

(*39)

M.birmanica

0

30.71

8.37

27

0

42.86

32.14

(0)

(33)

(11)

(20)

(31)

(43)

(32)

L.pussilus

*32.99

0

0

0

*59

0

0

(31)

(0)

(0)

(12)

(0)

(0)

O.beatrix

0

0.56

36.5

0

0

0

0

(0)

(0.3)

(24)

(0)

(0)

(0)

(0)

P.excavatus

16.5

0

0

0

0

0

0

(10.1)

(0)

(0)

(0)

(0)

(0)

(0)

Table 6 Earthworm species dominance (%) Abundance & Biomass (Values in parenthesis) within the same land use types in Garhwal Himalayas

OF

AL(T)

AL(WT)

AAL

RAL

DFL

RFL

OF

1

0.248

0.08

0 .088

0.72

0.252

0.51

AL(T)

1

0.688

0.91

0.32

0.491

0.625

AL(WT)

1

0.75

0.273

0.181

0.15

AAL

1

0.315

0.24

0.311

RAL

1

0.142

0.175

DFL

1

0.81

RFL

1

Table 7 Morisita’s index of similarity for earthworm communities between different land use types

 

A. alexandri

B.parvus

D. nepalensis

M.anamola

M.birmanica

L.pussilus

O.beatrix

P.excavatus

 

M
(%)

T
(°C)

M
(%)

T
(°C)

M
(%)

T
(°C)

M
(%)

T
(°C)

M
(%)

T
(°C)

M
(%)

T
(°C)

M
(%)

T
(°C)

M
(%)

T
(°C)

OF

-0.03

0.30

**0.96

 

**0.95

A

A

0.34

0.41

A

A

*0.89

*0.87

A

A

**0.92

*0.9

AL(T)

-0.03

0.40

A

A

*0.88

*0.84

**0.97

0.55

**0.98

*0.96

A

A

0.62

**0.92

A

A

AL(WT)

*0.88

*0.87

A

A

*0.85

*0.87

**0.99

*0.87

0.56

0.18

A

A

0.58

**0.95

A

A

AAL

A

A

A

A

0.19

*0.81

**0.96

*0.97

**0.97

*0.95

A

A

A

A

A

A

RAL

**0.98

*0.95

A

A

*0.88

*0.83

**0.98

*0.96

A

A

*0.85

*0.89

A

A

A

A

DFL

**0.97

*0.96

A

A

A

A

A

A

**0.97

A

A

A

A

A

A

A

RFL

0.31

0.27

A

A

A

A

-0.28

-0.35

*.85

-0.38

A

A

A

A

A

A

Table 8 Correlation coefficient(r) for soil moisture (M%); temperature (T (°C) and earthworm species abundance in Garhwal Himalaya
*indicates P<0.05 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaiwdaaaa@3E2C@ , **indicates P<0.01 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbdfwBIj xAHbqedmvETj2BSbqefm0B1jxALjhiov2DaerbuLwBLnhiov2DGi1B TfMBaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaqaafaaake aajugibabaaaaaaaaapeGaamiuaiabgYda8iaaicdacaGGUaGaaGim aiaaigdaaaa@3E28@

Discussion

Consequences of deforestation arise from site degradation leading to strong modification of soil properties this in turn can significantly affect both incidence and abundance of soil macro fauna. Earthworm communities are more directly altered by this changes.19 Endemic and exotic species co existed in the study area following deforestation and intensive cultivation. That native species dominate the undisturbed sites and disturbance and degradations leads to invasion by the exotic species20 holds true in our study. The sites under study represented the degraded areas as none of the species reported from the present experimental plots were endemic to the region, all the species are either peregrine exotic or peregrine endemic to the area as many of the endemic species of this region probably exterminated during the last Quaternary Glaciation.21 The elimination of old secondary forest and its replacement with agroecosystem also led to changed species composition due to altered habitat with peregrine exotics M. birmanica, O. beatrix and peregrine natives D nepalensis replacing peregrine exotics B. parvus, and natives L. pussilus and P. excavates in AL (T) and in AL (WT), similar results have also been shown through the studies.22 Extreme degradation of agriculture ecosystems due to faulty land management practices probably led to the loss of peregrine exotics A. alaxandri, M. anamola and O. beatrix in AAL as has also been shown through the studies19 in land use types in Northern Kwazulu Natal South Africa. The disappearance of peregrine exotics A. alaxandri and peregrine natives L. pussilus in RAL could probably be related to changed vegetation and edaphic conditions in RAL as has also been reported by Fragoso C et al.23 for earthworm communities in disturbed natural systems of tropical east Mexico and rehabilitation of RAL through various trees and crops probably led to invasion of peregrine exotic M. birmanica and M. anamola in these area. Higher level of degradation leading to replacement of OF to degraded forest land led to loss of exotic B. parvus, M. anamola and natives L. pussilus, P. excavatus. Resistance to invasion by endemic species in DFL and RFL could be a function of physical and chemical characteristics of the site24 and the unsuitability of the habitat probably impeded the invasion by the endemic species in the ecosystems under study.25 The absence of original species composition in RAL and RFL as compared to OF even after a period of 15 to 20 years suggest that this is probably still the secondary successional stage over a time scale of 15 years.26

  1. Total density and biomass
  2. Land-use alteration generally results in changes of vegetation and these changes have a significant effect on soil macrofauna.4 Study done by19 in South Africa showed most land uses supported between five and seven species. but the number of species present at our sites were lower ranging between 2 to 5. The statement that earthworm populations in cultivated lands are generally lower than those found in undisturbed habitats27 does not hold true in our studies, the traditional agriculture probably have a positive effect on earthworms through improved food supply as a result of recycling of crop residues, and organic manure added to the soil resulting in loosening of soil to an extent that facilitates burrowing by earthworms.28 This also explains the increased abundance and biomass of earthworm in AL(WT) as compared to OF. With lower bulk density and higher carbon percentage, coupled with reduced disturbance this effect was more prominent in AL(T). In RAL the supplemental irrigation and organic manure along with the planted trees following traditional farming practices favored the increase in total earthworm density as well as biomass. The smaller community under AAL and DFL than OF is attributable to loss of species when the land was cleared and cultivated. More stress full environment in DFL as compared to AAL due to heavy grazing pressure and lower soil moisture2 probably caused significant difference in the density and biomass of earthworm between AAL and DFL. In the present study the numerical abundance of earthworms in OF is lower compared to similar land use type from Kumoan Himalayas. Besides variation in ecological characteristics probably the land use history of the ecosystems contributes to this variation, the sub-temperate climate in the study area and the relatively lower inputs of organic matter apparently provide poor conditions for earthworm populations to flourish. As a result total numbers were very low ranging from 8 to 348 m-2 as compared to high ranging from 250 to 2400 m-2 in Northern Kwazulu Natal South Africa.19 Earthworm biomass showed broadly similar trends with land use to those for abundance.

  3. Interhabitat variation
  4. Interactions between species and sensitivity to ecological factors, presence or absence of and changes in ground vegetation composition are known to affect the composition of earthworm communities through changes in the distribution and the quality of litter, soil climate, and water availability. The presence of litter layer and lower perturbation pressure probably explains the numerical dominance of L. pussilus in the OF and the higher biomass of A. alexandri could be due to the larger size of the earthworm. B. parvus, L. pussilus and P. excavatus are litter-associated taxa which were more directly affected by OF clearance and the resulting decrease in available litter, thus explaining their disappearance in the changed ecosystems, however improved soil moisture and temperature as well as input of organic matter in RAL could probably be favorable factor for decolonization and dominance of L. pussilus. A. alexandri has wider ecological amplitude occurring under all land use types. Major determinants of earthworm communities’ structure in an agro ecosystem are the quantity and quality of organic matter added, soil type and the perturbation pressure.29 With better adaptation and tolerance to various disturbances during agro forestry practices. D. nepalensis was confined only to agroecosystems and was numerically dominant during cropping in AL (WT), in AL (T) and also in AAL. M. anamola was more directly affected due to perturbation pressure and soil degradation caused through conversion of forest to agro ecosystem,. The increased population density of A. alexandri, D. nepalensis and M. anamola in AL (T) as compared to AL(WT) is likely due to amelioration of the surface soil temperature and moisture by litter, and tree leaf biomass incorporated into the soil. The lower density and biomass of D. nepalensis under AL(WT) as compared to AL(T) could be the result of perturbations caused to the soil due to agricultural practices besides here the decline in D. nepalensis abundance in AAL and its subsequent increase in RAL could probably be because of lower soil bulk density and higher carbon percentage here and also because this land use categories maintain a year round canopy and litter layer.2 The lower biomass of A. alexandri in RFL may be due to species specific competition between M. anamola and A. alexandri as both occupy the same niche. Most changes in the distribution of earthworm species are explained by microclimate variations in soils thus low soil water content, high soil temperature, and incident radiation probably resulted in the decline in abundance and biomass of M. birmanica. in AAL and DFL as has also been shown through the studies of30 in the land use systems in eastern Zambia however improvement of soil conditions through the rehabilitation of forested land led to recurrence of this species. The traditional farming practices under rain fed conditions probably favored the population of O. beatrix species and thus explains its presence in agro ecosystem only.

  5. Similarity index
  6. Our results show that the AAL remains closer to the AL (T) than to the forest, because in these land use types the overall vegetation diversity remains low corresponding to a low diversification of the organic resources thus explaining the similarity between the communities in land use types.31 The shift in plant composition in RAL resulted in a shift in organic input from a below ground pattern in AAL to an above and below ground input here. This increase in above ground litter input created habitat in the plantation floor similar to that of OF and this probably was the reason for community structure of RAL to similar to the OF as has also been shown through the studies done by workers32 who suggested that the community structure under tree plantations is typically similar to that under native forest. The rehabilitation of DFL only caused an increase in the abundance of earthworm species in RFL resulting in community structure of RFL and DFL being closer to one another.

  7. Functional guild
  8. The size and species composition of earthworm communities are important because any shifts in earthworm community may result in significant changes in soil properties, epigeic earthworms, which usually inhabit upper soil and litter layers, may be more exposed than endogeic and anecic species that live primarily in subsurface soil. The conversion of OF to agro ecosystems led to a shift in functional categories from epigeic dominated species community to endogiec and anecic dominated composition, as has also been shown through studies of33 where the loss of the surface litter layer when rain forest is converted to agricultural use resulted in a dramatic decrease in the number of epigeic and anecic species and increase in numbers by endogeics species. The epigeics are litter-associated and thus were more directly affected by forest clearance and the resulting decrease in available litter leading to their loss in all the experimental plots except in RAL where with the deposition of litter the epiges recurred. P. excavates and B. parvus species cannot therefore survive in areas with less plant cover and litter availability and thus can be said to be bioindicator of land use that leads to these conditions. Anecic species need plant litter and specific microclimatic conditions in the soil thus the cropping activities with the removal of litter layer in AL (WT) prior to rice crop plantation resulted in decline of the anecics here. The increase in above ground litter input in AL (T) created microhabitat in the below canopy promoting the colonization of both endoges and anecic species. Further anecics are generally sparse in tropical environments,34 the experimental studies in the rehabilitated plots by Butt KR et al.35 have shown anecics to be slow colonizers in the reclaimed sites and this could probably explain their lower density in the RAL and RFL. Endogeics being geophagus, their predominance in the AL (T), AL (WT) AAL and RAL, could be because probably these habitats offered favorable base resource to the species which migrated from the surrounding plots to colonize and establish themselves in the new habitats.36,37 Dominance of endogeics in the disturbed habitats has also been shown through the studies of.38,39 As a consequence, while epigeics does not significantly alter the soil surface, casts produced by endogeics highly modify soil surface properties.

  9. Diversity index
  10. As reported in our studies conversion of OF to other land use types resulted in lower species diversity, similar result were also obtained through the studies of40 where conversion of forests to pastures resulted in decline of earthworm species diversity. Rehabilitation of these ecosystems in RAL and RFL resulted in improved diversity index of the community structure, this increase in earthworm diversity could probably be related to a shift in organic input from a below ground pattern in AAL and DFL to an above and below ground input in RAL and RFL.41 The lower overall vegetation diversity in AAL probably corresponded to a very low diversification of the organic resources and this can explain a lower diversity of earthworms here as low resource diversity leads to impoverished species diversity.38

  11. Seasonal variation

Seasonal variation has been shown for macro faunal assemblages for temperate,42 for tropical12 and for semiarid region,43 similar seasonal rhythm for all the species was also observed in our studies. Earthworm exhibit temporal variability, population density and biomass as well as the average depth of an organisms’ position in the soil profile is greatly affected by soil temperature and humidity.44,45 such seasonal rhythms of soil macro fauna were clearly expressed in the experimental land use type. Large differences in the number of earthworms were found in the different seasons highest catches were in monsoon, with lower numbers in winter and a large reduction in summer. Many soil organisms display strong seasonality in their life cycles,46 Few earthworms were found in summer compared to monsoon months , suggesting that all species hatch during early monsoon grow to adulthood over the monsoon period and decline in winter months, the lower mean catches in summer can be attributed to a response to drought by some species, the increase in the population density and the biomass of all earthworm species during the rainy season could be attributed to better soil moisture and temperature condition which favored the growth and maturity of earthworms.47 Seasonal temperature variations commonly induce vertical movements of earthworms in as soil profile48,49 as was also observed in D. nepalensis in the AL(T) which retreated to deeper soil50 to escape lower temperature and draught by diapausing during winter.

Conclusion

This study has shown that earthworm communities are directly affected by alteration of land use patterns. The numbers of species present at our sites ranged between 2 to 5. All the ecological categories such as endemic and exotics co-existed in the study area following deforestation and intensive cultivation, and the disturbance and degradations led to invasion by the exotic species and the native species dominated the undisturbed sites. The sites under study represented the degraded areas and none of the species reported from the present experimental plots were endemic to the region, all the species are either peregrine exotic or peregrine endemic to the area, as many of the endemic species of this region probably exterminated during the last Quaternary Glaciation. The presence of litter layer and lower perturbation pressure resulted in the numerical dominance of epigeic L. pussilus in the Oak forest (OF) and the higher biomass of A. alexandri may be because of the larger size of the earthworm. Litter-associated taxa (Epigeics) such as B. parvus, L. pussilus and P. excavatus are more directly affected by primary OF clearance and the resulting decrease in available litter, thus explaining their disappearance in the changed ecosystems. However improved soil moisture and temperature as well as input of organic matter in rehabilitated agricultural land (RAL) were probably one of the favorable factor for decolonization and dominance of epigeic L. pussilus. Endogeic species probably are not much affected by deforestation and degradation process and therefore A. alexandri had a wider ecological amplitude occurring under all land use types, however conversion of primary OF to other land use types resulted in lower species diversity. Our results show that the earthworm diversity in abandoned agriculture land (AAL) remains closer to the AL(T) than to the forest because in these land use types the overall vegetation diversity remains low corresponding to a low diversification of the organic resources thus explaining the similarity between the earthworm communities in land use types. All the species identified at study sites showed seasonal rhythmic pattern.

Acknowledgement

The research was financially supported by TSBF-SARNET programme. Thanks are due to Dr JM Julka (Emeritus Scientist) of Zoological Survey of India, Govt. of India for taxonomic identification of the earthworm species. The authors acknowledge the help extended by the co-ordination unit based at JNU New Delhi for providing the necessary literature.

Conflict of Interest

There is no financial or any conflict of interest exists.

References

  1. Semwal RL, Maikhuri RK, Rao KS, et al. Crop Productivity under Differently Lopped Canopies of Multipurpose Trees in Central Himalaya in India. Agroforestry Systems. 2002;56(1):57–63.
  2. Maikhuri RK, Rao KS, Semwal RL. Rehabilitation of Degraded Lands through Participatory Approach and Action Research: A case study from Central Himalaya. Asscociation Annual Journal. 2001;1(2):36–45.
  3. Maikhuri RK, Semwal RL, Rao KS. Rehabilitation of Degraded Community Lands for Sustainable Development in Himalayas: A case study in Garhwal Himalaya, India. International Journal of Sustainable Development & World Ecology. 1997(4):192–203.
  4. Pauli N, Barrios E, Conacher AJ, et al. Soil Macrofauna in Agricultural Landscapes Dominated by the Quesungual Slash–and–Mulch Agroforestry System, Western Honduras. Applied Soil Ecology. 2010;47(2):119–132.
  5. Barrios E. Soil Biota. Ecosystem Services and Land Productivity. Ecology Economics. 2007;64(2):269–285.
  6. Tao HH, Slade EM, Willis KJ, et al. Effects of Soil Management Practices on Soil Fauna Feeding Activity in an Indonesian oil Palm Plantation. Agriculture, Ecosystems & Environment. 2016;218:133–140.
  7. Baretta D, Bartz MLC, Fachini I, et al. Soil Fauna and its Relation with Environmental variables in Soil Management Systems. Revista Ciencia Agronomica. 2014;45(5):1806–6690.
  8. Capowiecz Y. Differences in Burrowing Behaviour and Spatial Interaction between the two Earthworm Species Aporrectodea nocturna and Allolobophora chlorotica.Biology and Fertility of Soils. 2000;30(4):341–346.
  9. Briones MJI, Mascato R, Mato S. Autoecological Study of some Earthworm Species (Oligochaeta) by Means of Ecological Profiles. Pedobiologia. 1995(39):97–106.
  10. Babel U, Ehrmann O, Krebs M. Relationships between earthworms and some plant species in a meadow. Soil Biology and Biochemistry. 1992;24(12):1477–1481.
  11. Jimenez JJ, Moreno AG, Lavelle P, et al. Population Dynamics and Adaptive Strategies of Martiodrilus carimaguensis (Oligochaeta, Glossoscolecidae), A native species from the well–drained savannas of Colombia. Applied Soil Ecology. 1998;9(1–3):153–160.
  12. Rossi JP, Blanchart E. Seasonal and Land–Use Induced Variations of Soil Macrofauna Composition in the Western Ghats, Southern India. Soil Biology & Biochemistry. 2005;37(6):1093–1104.
  13. Rao PN, Pati UC. Geology and Tectonics of Bhilangana Valley and its Adjoining Parts, Garhwal Himalaya with Special Reference to the Main Central Thrust. Himalayan Geology. 1980;10:220–237.
  14. Okalebo JR, Gathua KW, Woomer PL. Laboratory Methods of Soil and Plant Analysis:A Working Manual. Kenya: ROSTA; 1994. 88 p.
  15. Bouyoucos GJ. A Recalibration of the Hydrometer Method for making Mechanical Analysis of the Soil. Agronomy Journal. 1951;43: 434–438.
  16. Walkey A, Black IA. An Examination of Degtjareff Method for Determining Soil Organic Matter and a Proposed Chromic Acid Titration Method. Soil Science. 1934;37(1):29–38.
  17. Anderson JM, Ingram JS. Tropical Soil Biology and Fertility: A Handbook of Methods. UK: CABInternational;1993:221.
  18. Zar JH. Biostatistical Analysis. USA: Prentice–Hall Inc; 1974:621.
  19. Steinwand M, Steiner BCS, Seebe GUH, et al. Effects of Alpine Land–Use Changes: Soil Macrofauna Community Revisited. Ecology and Evolution. 2017;7(14):5389–5399.
  20. Korell L, Lang BR, Hensen I, et al. Interaction count: plant origin, herbivory and disturbance jointly explain seedling recruitment and community structure. Scientific Reports. 2017;7:82–88.
  21. Julka JM, Paliwal R. Distribution of Earthworms in Different Agroclimate Regions of India. In: Ramakrishnan PS, Saxena KG, et al. editors.Soil Biodiversity: Ecological processes and Landscape Management. India: Oxford and IBH Publishing Co. Pvt. Ltd; 2005.
  22. Myers N, Knoll A. The Biotic Crisis and the Future of Evolution. Proceedings of the National Academy of Sciences of United States of America. 2001;98(10):5389–5392.
  23. Fragoso C, Fernandez PR. Earthworms from South Western Mexico, New Acanthodriline genus and species (Megascolicidae–Oligochaeta). Megadrilogica. 1994;6(1):1–12.
  24. Hendrix PF, Baker GH, Callaham MA, et al. Invasion of Exotic Earthworms into Ecosystems Inhabited by Native Earthworms. Biological Invasion. 2006;8(6):1287–1300.
  25. Hendrix PF, Bohlen P. Exotic earthworm invasion in North America: Ecological and Policy Implication. Bioscience. 2002;52(9):801–811.
  26. Folgarail PJ, Thomas F, Desjardins T, et al. Soil Properties and the Macrofauna Community in Abandoned Irrigated Rice Fields of North Eastern Argentina. Biology and Fertility of Soils. 2003;38(6):349–357.
  27. Tunira B, Pradeep K, Rakesh M, et al. Effect of Application of Vermicompost and conventional Compost Derived from Different Residues on Pea Crop Production and Soil Faunal Diversity in Agricultural System in Garhwal Himalayas India. Natural Science. 2014;6(6):433–446.
  28. Lagerlof G, Goffre B, Vincent C. The Importance of Field Boundaries for Earthworms in the Swedish Agriculture Landscape. Agriculture, Ecosystem Environment. 2002;89(1–2):91–103.
  29. Fraser PM. The Impact of Soil and Crop Management Practices on Soil Macrofauna. In: Pankhurst CE, Doube MM, Gupta VSR, et al. editors. Soil Biota Management in Sustainable Farming Systems, Transactions of 15th World Congress of Soil Science. Mexico: Acapulco; 1994(4a):55–56.
  30. Sileshi G, Mafongoya P. Variation in Macrofauna Communities under Contrasting Land Use Systems in Eastern Zambia. Applied Soil Ecology. 2006;33(1):49–60.
  31. Basu P, Blanchart E, Lepage M. Termite (Isoptera) Community in the Western Ghats, South India: Influence of Anthropogenic Disturbance of Natural Vegetation. European Journal of Soil Biology. 1996;32(3):113–121.
  32. Fragoso C, Kanyonyo J, Moreno A, et al. A Survey of Tropical Earthworms: Taxonomy, Biogeography and Environmental Plasticity. In: Lavelle P, Brussaard L, Hendrix P, editors. Earthworm Management in Tropical Agroecosystems. UK: CABI;1999:1–26.
  33. Fragoso C, Lavelle P. Earthworm Communities of Tropical Rain Forests. Soil Biology and Biochemistry. 1992;24(12):1397–1408.
  34. Barrios I, Lavelle P, Brossard M, et al. Ecology of Earthworm Species with Large Environmental Tolerance and/or Extended Distributions. In: Lavelle P, Brussaard L, Hendrix P, editors. Earthworm Management in Tropical Agroecosystems. UK: CABI;1999:57–85.
  35. Butt KR, Lowe CN, Frederickson J, et al. The Development of Sustainable Earthworm Population at Calvert Landfill site, UK. Land Degradation and Development. 2004;15(1): 27–36.
  36. Pizl V. Earthworm Succession in Afforested Colliery Spoils Heaps in the Sokolov Region, Czech Republic. Restoration Ecology.2001;9(4):359–364.
  37. Winsome T, Epstein L, Hendrix PF, et al. Competition between Native and Exotic Earthworms Species in California Grassland.Applied Soil Ecology. 2006;32(1):38–53.
  38. Blanchart E, Julka JM. Influence of Forest Disturbance on Earthworm (Oligochaeta) Communities in the Western Ghats (South–India). Soil Biology and Biochemistry. 1997;29(3–4):303–306.
  39. Bhadauria T, Ramakrishnan PS. Impact of Land Use Land Cover Change on Earthworm Community Structure in India In: Ramakrishnan PS, Saxena KG, Swift MJ, et al. editors. Soil Biodiversity: Ecological processes and Landscape Management. India: Oxford and IBH Publishing Co. Pvt. Ltd; 2005.
  40. Rombke J, Verhaagh M. About Earthworm Communities in a Rain Forest and an Adjacent Pasture in Peru. Amazoniana. 1992;12(1):29–49.
  41. Zou X, Gonzalez G. Changes in Earthworm Density and Community Structure during Secondary Succession in Abandoned Tropical Pastures. Soil Biology & Biochemistry. 1997;29(3–4):627–629.
  42. Berg MP, Kniese JP, Bedaux JJM, et al. Dynamics and Stratification of Functional Groups of Micro and Mesoarthropods in the Organic Layer of a Scots Pine Forest. Biology and Fertility of Soils. 1998;26(4):268–284.
  43. Miranda ED, Pinero FS, Megias GA. Soil Macro Invertebrate Fauna of a Mediterranean Arid System: Composition and temporal changes in the assemblage. Soil Biology and Biochemistry. (2007);39(8):1916–1925.
  44. Fayolle L, Michaud H, Cluzeau D, et al. Influence of Temperature and Food Source on the Life Cycle of the Earthworm Dendrobaena veneta (Oligochaeta). Soil Biology & Biochemistry. 1997;29(3–4):747–750.
  45. Dibog L, Eggleton P, Forzi F. Seasonality of Soil Termites in a Humid Tropical Forest, Mbalmayo, Southern Cameroon. Journal of Tropical Ecology. 1998;14(6):841–849.
  46. Measey GJ, Gower DJ, Oommen OV, et al. A Subterranean Generalist Predator: Diet of the soil–dwelling caecilian Gegeneophis ramaswamii (Amphibia; Gymnophiona; Caeciliidae) in southern India. Comptes Rendus Biologies. 2004;327(1):65–76.
  47. Campana C, Gauvin S, Ponge JF. Influence of Ground Cover on Earthworm Communities in an Unmanaged Beech Forest: linear gradient studies. European Journal of Soil Biology. 2002;38(2):213–224.
  48. Luxton M. Studies on the Orbatid Mites of a Danish Beech wood Soil IV, Developmental Biology. Paedobiologia. 1981;21:312–340.
  49. Didden WAM. Ecology of Terrestrial Enchytraedae. Pedobiologia. 1993;37: 2–29.
  50. Morgan JE, Morgan AJ. Seasonal Changes in the Tissue–Metal (Cd, Zn and Pb) Concentrations in two Eco–physiologically Dissimilar Earthworm Species: Pollution–monitoring implications. Environmental Pollution. 1993;82(1):1–7.
Creative Commons Attribution License

©2018 Bhadauria, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.