*Article* **Earthworm Diversity, Forest Conversion and Agroforestry in Quang Nam Province, Vietnam**

**Rachmat Mulia <sup>1</sup> , Sam Van Hoang 2,\*, Van Mai Dinh <sup>3</sup> , Ngoc Bich Thi Duong <sup>4</sup> , Anh Duc Nguyen 5,6 , Dang Hai Lam <sup>7</sup> , Duyen Thu Thi Hoang <sup>8</sup> and Meine van Noordwijk <sup>9</sup>**


**Abstract:** The conversion of natural forests to different land uses still occurs in various parts of Southeast Asia with poor records of impact on ecosystem services and biodiversity. We quantified such impacts on earthworm diversity in two communes of Quang Nam province, Vietnam. Both communes are situated within buffer zones of a nature reserve where remaining natural forests are under threat of continued conversion. We identified 25 different earthworm species, out of which 21 were found in natural forests, 15 in agroforestry, 14 in planted forests, and seven each in annual croplands and home gardens. Out of the six species that were omnipresent inhabitants of all observed habitats, *Pontoscolex corethrurus* largely dominated habitats with intensive anthropogenic activities but was rare in natural forests. Natural and regenerated forests had a much denser earthworm population in the top 10 cm of soil rather than in deeper soil layers. We conclude that the conversion of natural forests into different land uses has reduced earthworm diversity which can substantially affect soil health and ecosystem functions in the two communes. Protection of the remaining natural forests is urgent, while the promotion of a tree-based farming system such as agroforestry can reconcile earthworm conservation and local livelihoods.

**Keywords:** land-use change; belowground biodiversity; soil engineers; *Pontoscolex corethrurus*; natural habitats; planted forest

### **1. Introduction**

The inclusion of 'planted forest' in forest statistics has been a major factor in the 'forest transition' concept: the reversal after long periods of decline towards a net gain in forest area [1,2]. Forest statistics can thus mask the ongoing loss of natural forest by larger gains in planted forest area, even though in terms of ecosystem services old-growth natural forests are very different from young planted forests [3]. Beyond the direct benefits that land-use change brings (or is expected to bring) to land users, a wide range of 'externalities' are generated, changes that do not play a role in the land users' decision making [4]. Where such externalities are negative for others in the landscape, corrective action is needed

**Citation:** Mulia, R.; Hoang, S.V.; Dinh, V.M.; Duong, N.B.T.; Nguyen, A.D.; Lam, D.H.; Thi Hoang, D.T.; van Noordwijk, M. Earthworm Diversity, Forest Conversion and Agroforestry in Quang Nam Province, Vietnam. *Land* **2021**, *10*, 36. https://doi.org/ 10.3390/land10010036

Received: 30 October 2020 Accepted: 30 December 2020 Published: 4 January 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

in the form of regulations or incentives; where they are positive, forms of incentives are appropriate [5].

Many of the ecosystem services that, as externalities of decision-making, are influenced by land use are controlled by soil conditions [6]. Beyond vegetation with its aboveground litter and plant root systems and turnover, microbes, fungi, and soil fauna are important for modifying the soil conditions that in turn shape soil structure (pores, aggregates) influencing water flows and storage and net greenhouse gas emissions that determine ecosystem services [7]. As part of the global loss of biodiversity, loss of belowground biodiversity and its impact on ecosystem services is still poorly understood [8,9]. Changes in the soil can be relatively slow to be noticed by soil biota, as the system is buffered in many aspects, but once undesirable changes take place, they are also hard to correct, especially if they reach soil layers below the surface [10], as thresholds are past where both roots and earthworms have difficulties penetrating [11]. Better knowledge of changes in belowground biodiversity in the context of 'forest transition' patterns of land-use change is urgent for many parts of the tropics, especially.

Over the past decades, Vietnam's forest cover has continued to expand from 28% in 1993 to 42% in 2018 [12]. Policy reforms and nation-wide forest protection and replanting programs have been the main drivers [13,14]. Forest expansion, including planted forests, has substantially contributed to economic growth, job creation, and poverty alleviation in the country [12]. However, despite the increasing trend, the conversion and degradation of natural forests still occur in different parts of the country. The existing natural forests in Vietnam are largely in poor condition or regenerating [12]. The forest cover rates overlook forest quality. Competing land uses, exploitation of forest resources, and weak forest governance are serious challenges for the country to protect the remaining and regenerating poor natural forests [12]. Massive conversion of natural forests still occurs for example in the Central Highlands and Central Annamites region of the country. The expansion of commercial crops and infrastructure development have been the main drivers of forest conversion in Central Highlands [15]. In the Central Annamites, considered as one of the largest continuous natural forest areas and key biodiversity landscapes in Asia, forest resources and wildlife are continuously under threat of intensive exploitation by surrounding forest-dependent populations [16,17].

Concerns over the impacts of forest conversions on ecosystem services and biodiversity have recently become more prominent in Vietnam [18,19]. Several sub-national and national policies to include the 2011–2020 National Biodiversity Strategy with vision to 2030 underline the need of taking up existing scientific- and evidence-based information and undertaking further research on such impacts while strengthening forest protection efforts. Several studies from other parts of the world have reported impacts of forest conversions on biodiversity including belowground biodiversity (e.g., [20,21]). Notwithstanding, accounts of such impacts in Vietnam are limited.

The diversity and activities of earthworms as one of the most important soil biotas influence the physical, chemical, and biological properties of soils [22,23]. Their physical movement creates soil pores which ease nutrient and water dynamic in the soils. As recyclers of organic materials, earthworms facilitate microorganisms such as fungi and bacteria to undertake further decomposing process [24]. Through these activities, earthworms influence soil health and the provision of ecosystem services [25]. The population and diversity of earthworm species vary across land habitats due to variation in soil moisture, soil temperature, soil properties, the abundance of surface litter, vegetation types, land use management, and human interventions [20,26–30]. Earthworms are thus sensitive to change in land use. In general, earthworm diversity is reduced in habitats with more intensive anthropogenic interventions [23,31].

Here, we describe the results of a study on earthworm diversity in natural forests and different land uses in two communes of Quang Nam province, South Central Coast region of Vietnam. The two communes are within the buffer zones of the Song Thanh Nature Reserve, one of the main reserves within the Central Annamites mountains, and parts of

main watersheds in the region. In both communes, the conversion of natural forests into shifting cultivations and agricultural systems, such as agroforestry and annual croplands, has been taking place for decades. The purpose of the study was to compare the diversity of earthworm species among natural forest and different land uses to investigate the impact of forest conversion and to recommend necessary measures related to earthworm conservation to maintain soil health and ecosystem functions in the two buffer-zone communes. tigate the impact of forest conversion and to recommend necessary measures related to earthworm conservation to maintain soil health and ecosystem functions in the two buffer-zone communes. **2. Material and Methods**  *2.1. Study Sites* 

parts of main watersheds in the region. In both communes, the conversion of natural forests into shifting cultivations and agricultural systems, such as agroforestry and annual croplands, has been taking place for decades. The purpose of the study was to compare the diversity of earthworm species among natural forest and different land uses to inves-

#### **2. Material and Methods** The study was conducted in Phuoc My commune of Phuoc Son district and Ta Bhing

#### *2.1. Study Sites* commune of Nam Giang district, Quang Nam province (15°35′0″ N–107°55′0″ E), South

The study was conducted in Phuoc My commune of Phuoc Son district and Ta Bhing commune of Nam Giang district, Quang Nam province (15◦3500 00 N–107◦5500 00 E), South Central Coast of Vietnam (Figure 1a). The two communes are within the buffer zones of Song Thanh Nature Reserve and parts of main watersheds in South Central Coast region which supply water for several cities including Da Nang as the biggest city in Central Vietnam. For example, Phuoc My is part of Dak Mi watershed (Figure 1b), while Ta Bhing is part of Vu Gia-Thu Bon as a bigger watershed [32]. Phuoc My is situated at 223–446 m above sea level while Ta Bhing is approximately 100 m above sea level. Both communes are mountainous with flat areas concentrated at the feet of the mountains, along riverbanks. Central Coast of Vietnam (**Error! Reference source not found.**Figure 1a). The two communes are within the buffer zones of Song Thanh Nature Reserve and parts of main watersheds in South Central Coast region which supply water for several cities including Da Nang as the biggest city in Central Vietnam. For example, Phuoc My is part of Dak Mi watershed (Figure 1b), while Ta Bhing is part of Vu Gia-Thu Bon as a bigger watershed [32]. Phuoc My is situated at 223–446 m above sea level while Ta Bhing is approximately 100 m above sea level. Both communes are mountainous with flat areas concentrated at the feet of the mountains, along riverbanks.

*Land* **2021**, *10*, x FOR PEER REVIEW 3 of 19

**Figure 1.** (**a**) Location of Phuoc My and Ta Bhing communes as study sites in Quang Nam province of Vietnam, (**b**) the boundary of Song Thanh Nature Reserve and Dak Mi watershed. **Figure 1.** (**a**) Location of Phuoc My and Ta Bhing communes as study sites in Quang Nam province of Vietnam, (**b**) the boundary of Song Thanh Nature Reserve and Dak Mi watershed.

The two communes have a tropical monsoon climate. The annual rainfall is about 2650 mm in Phuoc My and 2300 mm in Ta Bhing [33]. In both communes, the rainy season usually peaks between September and November, and the dry season lasts from January until March. The highest temperature is usually recorded between May and July. The two communes have similar soil physical and chemical properties, dominated by sandy-loam and sandy-clay-loam textures [33]. The two communes have a tropical monsoon climate. The annual rainfall is about 2650 mm in Phuoc My and 2300 mm in Ta Bhing [33]. In both communes, the rainy season usually peaks between September and November, and the dry season lasts from January until March. The highest temperature is usually recorded between May and July. The two communes have similar soil physical and chemical properties, dominated by sandy-loam and sandy-clay-loam textures [33].

Phuoc My has a smaller land size, but it has a larger area of rich natural forest than Ta Bhing in 2016. The total area of medium, poor, and regenerated natural forests accounts for 21% and 41.5% of the total land size in Phuoc My and Ta Bhing respectively. Land cover in the two communes has been classified into eleven types [16]. The area of each land cover in the two communes are described in Table A1. Phuoc My has a smaller land size, but it has a larger area of rich natural forest than Ta Bhing in 2016. The total area of medium, poor, and regenerated natural forests accounts for 21% and 41.5% of the total land size in Phuoc My and Ta Bhing respectively. Land cover in the two communes has been classified into eleven types [16]. The area of each land cover in the two communes are described in Table A1.

#### *2.2. Land-Use History*

In both communes, massive conversion of natural forests took place more than 20 to 30 years ago [16]. Over the past two decades, natural forests have regenerated in some

areas of the communes, thanks to enhanced protection efforts by the local authorities. Planted forests have become popular in the past decade, especially the short-rotation (3–4 year) of Acacia (mainly the hybrid *Acacia mangium x auriculiformis*) plantation for pulp and paper. Two types of agroforestry have developed in the two communes: a *taungya* or temporary agroforestry in which annual crops such as cassava or maize become intercrops in young planted forests, usually in the first 1–2 years after tree planting before tree canopy closure; and more permanent agroforestry in which timber tree species, such as *Melia azedarach*, *Vernicia montana*, or *Ficus racemosa,* or fruit trees such as *Dimocarpus longan,* were intercropped with annual crops such as peanut, beans, or cassava. In sloping lands, agroforestry practices with rows of trees and rows of pineapple to limit soil erosion also existed. areas of the communes, thanks to enhanced protection efforts by the local authorities. Planted forests have become popular in the past decade, especially the short-rotation (3– 4 year) of Acacia (mainly the hybrid *Acacia mangium x auriculiformis*) plantation for pulp and paper. Two types of agroforestry have developed in the two communes: a *taungya* or temporary agroforestry in which annual crops such as cassava or maize become intercrops in young planted forests, usually in the first 1–2 years after tree planting before tree canopy closure; and more permanent agroforestry in which timber tree species, such as *Melia azedarach*, *Vernicia montana*, or *Ficus racemosa,* or fruit trees such as *Dimocarpus longan,* were intercropped with annual crops such as peanut, beans, or cassava. In sloping lands, agroforestry practices with rows of trees and rows of pineapple to limit soil erosion also existed.

In both communes, massive conversion of natural forests took place more than 20 to 30 years ago [16]. Over the past two decades, natural forests have regenerated in some

*Land* **2021**, *10*, x FOR PEER REVIEW 4 of 19

#### *2.3. Sample Plots by Habitat 2.3. Sample Plots by Habitat*

*2.2. Land-Use History* 

The earthworm's observations were carried out in seven land cover types (hereafter called habitats). Three represent habitats with less anthropogenic interventions namely natural forests, regenerated forests, and grasslands (hereafter referred to as 'natural habitats'); and four with more intensive anthropogenic interventions namely planted forests, agroforestry, mixed annual crops, and home gardens (hereafter referred to as 'human-disturbed habitats'). Among the four, the most frequent and intensive anthropogenic interventions occurred in mixed annual crops and home gardens. The earthworm's observations were carried out in seven land cover types (hereafter called habitats). Three represent habitats with less anthropogenic interventions namely natural forests, regenerated forests, and grasslands (hereafter referred to as 'natural habitats'); and four with more intensive anthropogenic interventions namely planted forests, agroforestry, mixed annual crops, and home gardens (hereafter referred to as 'human-disturbed habitats'). Among the four, the most frequent and intensive

In total, we took 28 and 53 sample plots from the different habitats in Phuoc My and Ta Bhing commune respectively. More sample plots were taken in natural and regenerated forests to anticipate a possibility of higher species diversity in those habitats. Although both communes have similarities in many of the factors that can influence earthworm diversity such as climate condition, soil properties, and landcover types, we took sample plots in both communes, but more samples plots in Ta Bhing than Phuoc My merely for time and resource efficiency. Most of the local assistants involved in the sampling activities in both communes were from Ta Bhing. The aim was to have replications from both communes, not to compare earthworm diversity between the two communes. The number of sample plot by habitat in the two communes are given in Table A2. Figure 2 illustrates the vegetation covers and types in the selected habitats. anthropogenic interventions occurred in mixed annual crops and home gardens. In total, we took 28 and 53 sample plots from the different habitats in Phuoc My and Ta Bhing commune respectively. More sample plots were taken in natural and regenerated forests to anticipate a possibility of higher species diversity in those habitats. Although both communes have similarities in many of the factors that can influence earthworm diversity such as climate condition, soil properties, and landcover types, we took sample plots in both communes, but more samples plots in Ta Bhing than Phuoc My merely for time and resource efficiency. Most of the local assistants involved in the sampling activities in both communes were from Ta Bhing. The aim was to have replications from both communes, not to compare earthworm diversity between the two communes. The number of sample plot by habitat in the two communes are given in Table

A2. Figure 2 illustrates the vegetation covers and types in the selected habitats.

**Figure 2.** *Cont*.

**Figure 2.** Vegetation cover and types in (**a**) natural forest, (**b**) regenerated forest, (**c**) agroforestry, (**d**) plantation of *Acacia mangium x auriculiformis*, (**e**) plantation of *Melia azedarach*, (**f**) plantation for *Machilus odoratissima* Nees, (**g**) home garden, (**h**) grasslands, and (**i**) upland annual crops. **Figure 2.** Vegetation cover and types in (**a**) natural forest, (**b**) regenerated forest, (**c**) agroforestry, (**d**) plantation of *Acacia mangium x auriculiformis*, (**e**) plantation of *Melia azedarach*, (**f**) plantation for *Machilus odoratissima* Nees, (**g**) home garden, (**h**) grasslands, and (**i**) upland annual crops.

#### *2.4. Identification of Earthworm Species 2.4. Identification of Earthworm Species*

Following [34], at each sample plot, earthworm specimens were collected from a 50 cm × 50 cm soil block of 30 cm depth (please see Figure 2e,h,i). A 50 cm × 50 cm supporting frame and a woody measure of 30 cm long were used to ensure identical surface areas and thicknesses of all soil blocks. The specimens were collected by carefully segregating the soils per 10 cm soil depth. All specimens were stored in cloth bags to keep them alive for further treatments in the laboratory. To facilitate species identification, all earthworms were separated from soils by Following [34], at each sample plot, earthworm specimens were collected from a 50 cm × 50 cm soil block of 30 cm depth (please see Figure 2e,h,i). A 50 cm × 50 cm supporting frame and a woody measure of 30 cm long were used to ensure identical surface areas and thicknesses of all soil blocks. The specimens were collected by carefully segregating the soils per 10 cm soil depth. All specimens were stored in cloth bags to keep them alive for further treatments in the laboratory.

cleansing with water. They were subsequently immersed in liquid containing 4% formalin for 6–12 h. Thereafter, all earthworms were transferred to fresh liquid of 4% formalin for long-term storage and morphological studies. We investigated morphological features of collected earthworms, such as intestinal caeca (presence and shape), male pores (presence of copulatory pouches), genital markings (type, location, and number), spermathecal pores (location, number, and type), septa (thin, thick or absent), and additional characteristics such as prostomium (probilobous, epilobous or tanylobous), first dorsal pore (location), intestinal origin, last heart, ovaries, and testes. Earthworm samples were scrutinized using stereo microscope Motic DM143-FBGG-C (Motic Company, Hongkong) with drawing tube and camera attached to the monitor. All specimens were compared and identified using original papers as described in [35–43]. All specimens with different morphological characters compared to all known species were named using their genera and sp. 1, 2, and so on. For biomass measurement, we used paper tissues to remove all liquids (formalin, water) and cast from the preserved earthworms before weighting using the electronic weight Sartorius GM612. To facilitate species identification, all earthworms were separated from soils by cleansing with water. They were subsequently immersed in liquid containing 4% formalin for 6–12 h. Thereafter, all earthworms were transferred to fresh liquid of 4% formalin for long-term storage and morphological studies. We investigated morphological features of collected earthworms, such as intestinal caeca (presence and shape), male pores (presence of copulatory pouches), genital markings (type, location, and number), spermathecal pores (location, number, and type), septa (thin, thick or absent), and additional characteristics such as prostomium (probilobous, epilobous or tanylobous), first dorsal pore (location), intestinal origin, last heart, ovaries, and testes. Earthworm samples were scrutinized using stereo microscope Motic DM143-FBGG-C (Motic Company, Hongkong) with drawing tube and camera attached to the monitor. All specimens were compared and identified using original papers as described in [35–43]. All specimens with different morphological characters compared to all known species were named using their genera and sp. 1, 2, and so on. For biomass measurement, we used paper tissues to remove all liquids (formalin, water) and cast from the preserved earthworms before weighting using the electronic weight Sartorius GM612.

#### *2.5. Indicators for Earthworm Dominance and Diversity 2.5. Indicators for Earthworm Dominance and Diversity*

We used seven indicators to assess the dominance, density, and diversity of earthworm species. Following [34], five indicators represent dominance and density, namely density and biomass per soil layer, quantity and biomass dominance, and occurrence frequency. Two indicators, the Shannon–Wiener index [44] and similarity index [45], measure the level of diversity of earthworm species within habitat and similarity in species composition among habitats, respectively. For the similarity index, we omitted the sub-species component in the calculation. The seven indicators are We used seven indicators to assess the dominance, density, and diversity of earthworm species. Following [34], five indicators represent dominance and density, namely density and biomass per soil layer, quantity and biomass dominance, and occurrence frequency. Two indicators, the Shannon–Wiener index [44] and similarity index [45], measure the level of diversity of earthworm species within habitat and similarity in species composition among habitats, respectively. For the similarity index, we omitted the sub-species component in the calculation. The seven indicators are described in Table A3.

#### described in Table A3. *2.6. Aboveground Litter Biomass as A Covariate*

*2.6. Aboveground Litter Biomass as A Covariate*  In this study, we focused on the impact of conversion from natural to humandisturbed habitats on earthworm diversity. Notwithstanding, to consider the effect of aboveground plant litter biomass on earthworm density (e.g., [21,46,47]), we used the annual production rate of aboveground plant litter biomass (hereafter referred to as "annual litter production") (Table 1) as a covariate when comparing indicators among In this study, we focused on the impact of conversion from natural to human-disturbed habitats on earthworm diversity. Notwithstanding, to consider the effect of aboveground plant litter biomass on earthworm density (e.g., [21,46,47]), we used the annual production rate of aboveground plant litter biomass (hereafter referred to as "annual litter production") (Table 1) as a covariate when comparing indicators among habitats. Due to a lack of data, we estimated the annual litter production from the annual production rate of aboveground

plant biomass (hereafter referred to as "annual biomass production") of the habitats. We used a factor of 0.49 calculated from [48] for the biomass-litter conversion and assumed the same factor applies for all habitats. The biomass production data were obtained from [16] that estimated the standing aboveground plant biomass of all habitats in the two study communes, except NF and AF, under different ages using the Rapid Carbon Stock Appraisal approach [49]. The annual biomass production was calculated from the standing biomass and plot ages, assuming a constant biomass production across the year. We used an annual litter production of 9.3 ton ha−<sup>1</sup> year−<sup>1</sup> from [50] for NF and 2.4 <sup>±</sup> 0.4 ton ha−<sup>1</sup> year−<sup>1</sup> calculated from [51] for AF. The earthworm and plant biomass assessment were not conducted using the same sample plots. Therefore, to create a variation of annual litter biomass among earthworm sample plots, we used random numbers from a uniform distribution to generate higher or lower litter biomass by a maximum of one standard error from the mean. We applied this for all habitats except for NF that lacks input data on the standard error of the mean.

**Table 1.** Annual production of aboveground plant litter biomass by habitat (source: calculated from [16]).


#### *2.7. Statistical Analysis*

We used the IBM SPSS Statistics version 23 to test the difference in earthworm indicators among habitats using one-way analysis of covariance (ANCOVA). For the statistical test, we assumed an independency among sample plots. For all data and statistical analysis, we combined the earthworm data from the two study communes and the difference between the two communes was not investigated. All graphs were produced using the Microsoft Excel.

#### **3. Results**

#### *3.1. Species Occupancy by Habitat*

A total of 25 earthworm species from 3 families and 6 genera were found in the two study communes. Of these, 23 species belong to the *Megascolecidae* family and 17 species to the *Amynthas* genus. Only one genus and one species were found in each of the other two families, *Moniligastridae* and *Rhinodrilidae*. Seven out of the 25 species, namely *Pontoscolex corethrurus*, *Amynthas aspergillum*, *Amynthas divitopapillatus*, *Amynthas modiglianii*, *Amynthas* sp.*1*, *Metaphire houlleti*, and *Polypheretima taprobanae* were omnipresent. However, only *P. corethrurus* dominated the habitats, except in NF and RF. The average frequency of occurrence of this species in human-disturbed habitats was 87%. Based on the occurrence frequency, *P*. *corethrurus* is classified as rare species in NF and uncommon species in RF. In NF, only *Drawida beddardi* had occurrence frequency above 25%. The occurrence frequency of earthworm species by habitat is given in Table A4.

#### *3.2. Earthworm Density and Biomass by Habitat*

The NF tends to have a lower earthworm density compared to other habitats (Figure 3a). However, no significant difference (*p*-value > 0.05) was found among habitats due to large standard errors. In terms of biomass, the natural habitats had higher earthworm biomass compared to human-disturbed habitats (*p*-value < 0.001), likely because of the larger earthworm's physical size. The ratio between earthworm biomass and density was much lower in the human-disturbed habitats (*p*-value < 0.001) with an average of 0.73 ± 0.14 g per

individual compared to 1.62 ± 0.23 g per individual in natural habitats. Related to quantity and biomass dominance, *P. corethrurus* was much more dominant in human-disturbed than natural habitats (*p*-value < 0.01) particularly due to extremely low dominance in NF (Figure 3b). The quantity and biomass dominance of all species by habitat is provided in Table A5. much lower in the human-disturbed habitats (*p*-value < 0.001) with an average of 0.73 ± 0.14 g per individual compared to 1.62 ± 0.23 g per individual in natural habitats. Related to quantity and biomass dominance, *P. corethrurus* was much more dominant in humandisturbed than natural habitats (*p*-value < 0.01) particularly due to extremely low dominance in NF (Figure 3b). The quantity and biomass dominance of all species by habitat is provided in Table A5. larger earthworm's physical size. The ratio between earthworm biomass and density was much lower in the human-disturbed habitats (*p*-value < 0.001) with an average of 0.73 ± 0.14 g per individual compared to 1.62 ± 0.23 g per individual in natural habitats. Related to quantity and biomass dominance, *P. corethrurus* was much more dominant in humandisturbed than natural habitats (*p*-value < 0.01) particularly due to extremely low

larger earthworm's physical size. The ratio between earthworm biomass and density was

*Land* **2021**, *10*, x FOR PEER REVIEW 7 of 19

**Figure 3.** (**a**) Earthworm density and biomass and (**b**) quantity and biomass dominance of *P. corethrurus*. **Figure 3.** (**a**) Earthworm density and biomass and (**b**) quantity and biomass dominance of *P. corethrurus*. (**a**) (**b**) **Figure 3.** (**a**) Earthworm density and biomass and (**b**) quantity and biomass dominance of *P. corethrurus*.

#### *3.3. Earthworm Density and Biomass by Soil Depth 3.3. Earthworm Density and Biomass by Soil Depth*

Both the natural and human-disturbed habitats generally had higher density of earthworm in the first 10 cm than deeper soil layers (*p*-value = 0.001) (Figure 4a). However, the decrease in the density by soil layer was more pronounced in natural habitats, especially in NF and RF, than human-disturbed habitats. In terms of biomass, both natural and human-disturbed habitats generally had comparable earthworm biomass in the first two soil layers, and lower biomass (*p*-value < 0.002) in the 20–30 cm soil depth (Figure 4b). Both the natural and human-disturbed habitats generally had higher density of earthworm in the first 10 cm than deeper soil layers (*p*-value = 0.001) (Figure 4a). However, the decrease in the density by soil layer was more pronounced in natural habitats, especially in NF and RF, than human-disturbed habitats. In terms of biomass, both natural and human-disturbed habitats generally had comparable earthworm biomass in the first two soil layers, and lower biomass (*p*-value < 0.002) in the 20–30 cm soil depth (Figure 4b). *3.3. Earthworm Density and Biomass by Soil Depth*  Both the natural and human-disturbed habitats generally had higher density of earthworm in the first 10 cm than deeper soil layers (*p*-value = 0.001) (Figure 4a). However, the decrease in the density by soil layer was more pronounced in natural habitats, especially in NF and RF, than human-disturbed habitats. In terms of biomass, both natural and human-disturbed habitats generally had comparable earthworm biomass in the first

two soil layers, and lower biomass (*p*-value < 0.002) in the 20–30 cm soil depth (Figure 4b).

The natural habitats had a higher species diversity than the human-disturbed habitats (*p*-value < 0.001) (Figure 5a). The weak dominance of *P*. *corethrurus* especially in

The natural habitats had a higher species diversity than the human-disturbed habitats (*p*-value < 0.001) (Figure 5a). The weak dominance of *P*. *corethrurus* especially in

*3.4. Species Diversity and Similarity Among Habitats*  **Figure 4.** Earthworm (**a**) density and (**b**) biomass by soil depth. **Figure 4.** Earthworm (**a**) density and (**b**) biomass by soil depth.

*3.4. Species Diversity and Similarity Among Habitats* 

#### *3.4. Species Diversity and Similarity Among Habitats*

The natural habitats had a higher species diversity than the human-disturbed habitats (*p*-value < 0.001) (Figure 5a). The weak dominance of *P*. *corethrurus* especially in NF and RF was associated with higher species diversity. Among the human-disturbed habitats, the number of earthworm species found in UC and HG was equal, however, the Shannon– Wiener index of UC was much lower due to the strong dominance of few species. *Land* **2021**, *10*, x FOR PEER REVIEW 8 of 19 NF and RF was associated with higher species diversity. Among the human-disturbed habitats, the number of earthworm species found in UC and HG was equal, however, the Shannon–Wiener index of UC was much lower due to the strong dominance of few

The NF and RF had a similar earthworm species composition (Figure 5b). Among the four human-disturbed habitats, only UC and HG had a level of similarity. In general, the species composition in natural and human-disturbed habitats was different, except between RF and AF. species. The NF and RF had a similar earthworm species composition (Figure 5b). Among the four human-disturbed habitats, only UC and HG had a level of similarity. In general, the species composition in natural and human-disturbed habitats was different, except between RF and AF.

**Figure 5.** (**a**) Shannon-Wiener and (**b**) similarity index among habitats. **Figure 5.** (**a**) Shannon-Wiener and (**b**) similarity index among habitats.

#### *3.5. Variation among Types of Planted Forests*  The Acacia plantation had more diverse earthworm species than Melia or Machilus *3.5. Variation among Types of Planted Forests*

plantation (Table 2). However, *P. corethrurus* was less dominant in Melia plantation because two other species, *Amynthas aspergillum,* and *Metapheretima tiencanhensis*, had an occurrence frequency above 25%. Also, earthworms in Melia plantation had a higher biomass and density ratio that reached 1.58 g per individual, compared to 0.3 and 0.5 g per individual in Acacia and Machilus plantations, respectively. **Table 2.** Earthworm indicators comparing three types of planted forest. The Acacia plantation had more diverse earthworm species than Melia or Machilus plantation (Table 2). However, *P. corethrurus* was less dominant in Melia plantation because two other species, *Amynthas aspergillum,* and *Metapheretima tiencanhensis*, had an occurrence frequency above 25%. Also, earthworms in Melia plantation had a higher biomass and density ratio that reached 1.58 g per individual, compared to 0.3 and 0.5 g per individual in Acacia and Machilus plantations, respectively.

**Earthworm Indicators Acacia Melia Machilus**  No. of sample plots 15 5 4 **Table 2.** Earthworm indicators comparing three types of planted forest.


#### found in tree-based farming systems like AF and PF. Other studies [52–57] reported **4. Discussion**

#### similar evidence. Improved micro-climate, e.g., lower soil temperature and higher soil humidity, minimum plot management practices such as tillage, and a higher supply of *4.1. Earthworm Diversity in Agroforestry*

organic matter from above- and below-ground litter all generate higher earthworm diversity in tree-based systems [58,59]. These factors also create contrast in earthworm diversity between areas nearby and far from the trees within AF, for example between In our study, among the human-disturbed habitats, higher earthworm diversity was found in tree-based farming systems like AF and PF. Other studies [52–57] reported similar evidence. Improved micro-climate, e.g., lower soil temperature and higher soil humidity, minimum plot management practices such as tillage, and a higher supply of organic matter from above- and below-ground litter all generate higher earthworm diversity in tree-based systems [58,59]. These factors also create contrast in earthworm diversity between areas nearby and far from the trees within AF, for example between tree row and crop alley in

the alley-cropping system. In addition, minimal application of chemical inputs such as fertilizer and pesticide contribute to generating higher earthworm diversity within tree rows [58,60]. The variation in micro-climate and plot management practices also influence earthworm density. For example, in poplar AF in Canada, the average earthworm density within tree row was 182 individuals m−<sup>2</sup> , compared to 117 and 95 individuals m−<sup>2</sup> within the crop alley at two and six meters away from tree row, respectively [60]. In our study, some of the observed AFs were alley cropping, and others had scattered trees over the plots. In both cases, soil blocks for earthworm observation were taken relatively close to the trees, to avoid causing damage to annual crops. This likely explains the comparable earthworm diversity between AF and PF. AF can become a solution if larger areas of annual crops are needed without further converting NF or RF which otherwise will further decline the ecosystem health and functions in the two communes.

#### *4.2. Potential Impact of P. corethrurus Dominance*

The dominance of *P. corethrurus* in human-disturbed habitats has also been reported in other studies especially in tropical areas, e.g., see review in [61]. The species was exotic to Vietnam and initially inhabited hilly areas of the country [62,63]. Nowadays, it has a large distribution area including in the coastal areas of the country [64]. Compared to other exotic species which deliberately or inadvertently colonize new habitats, *P. corethrurus* has a high survival rate thanks to its tolerance to a range of biotic and abiotic environment [24,65]. Many other species either fail to survive or survive but not invasive [66].

*P. corethrurus* could bring favorable or unfavorable impacts to new habitats, although generally considered harmful to the environment and native species. As summarized in [61], as favorable impacts, the species could for example stimulate nutrient release in soil and enhance plant's resistance to phyto-parasitic nematodes. However, as unfavorable impacts, it can create soil compaction and alter biogeochemical processes to affect plants, native earthworms, and microbial communities. *P. corethrurus* can compact the soils through its feeding activity which accumulates small soil aggregates. In the absence of intervention from other soil biota, the small aggregates will progressively transform into larges aggregates and become compact [67,68]. In addition, layers of cast produced by this species on the soil surface, can turn into thick crust preventing water and air penetration to the soils [69]. The role of *P. corethrurus* in increasing soil bulk density and reducing soil porosity were reported in (e.g., [52,67,70,71]). Due to the unfavorable impacts, the presence and dominance of *P. corethrurus* have been considered as an indicator for soil health and level of disturbance to ecosystem and ecosystem services [72].

Our study shows the dominance of *P. corethrurus* in human-disturbed habitats namely PF, AF, UC, and HG which occupy about 17% and 20% of total land area in Phuoc My and Ta Bhing commune, respectively. These land uses potentially expand in the two communes driven by augmenting population and livelihood pressure. The two communes are therefore under a serious threat of further loss of earthworm diversity and stronger dominance of *P. corethrurus*, which can further reduce their ecosystem health and functions. Under the increasing dominance of the cosmopolitan species, the soil porosity and ecosystem functions in the two communes will mainly rely on tree roots. Therefore, tree-based farming systems such as AF are preferable to reconcile livelihood and provision of ecosystem services. Moreover, this land-use system generally has higher earthworm diversity than sole annual crop systems. Also, there is a need to select more suitable tree species for PF because of its popularity as land use in the two communes. Although Acacia is currently the most popular type of PF supported by local pulp and paper industries, Melia is worth to get further attention and research due to the lower dominance of *P. corethrurus* in the PF using this tree species. Based on the unpublished data as part of a study by [16], Melia in the communes likely had a superficial rooting system (Figure 6) which can create microand macropores in superficial soil layers through its coarse and fine roots. It has a root shallowness index of 0.83, higher than Acacia (0.69), and Machilus (0.63). The index measures the ratio of cross-sectional areas of horizontal and all proximal roots. Horizontal proximal

roots are defined as those with an angle of less than 45 degrees relative to the soil surface. The number of replications in the root study was however limited, only involved three different trees, and worth for further investigation. Therefore, promoting this species in the communes can contribute to maintaining soil porosity with less dominance of *P. corethrurus*. Also, in terms of topographical, soil, and climate condition, Melia is suitable for the two study communes [16]. horizontal and all proximal roots. Horizontal proximal roots are defined as those with an angle of less than 45 degrees relative to the soil surface. The number of replications in the root study was however limited, only involved three different trees, and worth for further investigation. Therefore, promoting this species in the communes can contribute to maintaining soil porosity with less dominance of *P. corethrurus*. Also, in terms of topographical, soil, and climate condition, Melia is suitable for the two study communes [16].

through its coarse and fine roots. It has a root shallowness index of 0.83, higher than Acacia (0.69), and Machilus (0.63). The index measures the ratio of cross-sectional areas of

(**a**) (**b**) **Figure 6.** Distribution of proximal roots of *Melia azedarach* in (**a**) home garden and (**b**) forest plantation. **Figure 6.** Distribution of proximal roots of *Melia azedarach* in (**a**) home garden and (**b**) forest plantation.

*Land* **2021**, *10*, x FOR PEER REVIEW 10 of 19

#### *4.3. Caveats in Assessing Earthworm-Habitat Linkage 4.3. Caveats in Assessing Earthworm-Habitat Linkage*

In this study, we could demonstrate that natural- and human-disturbed habitats have different levels of earthworm density and diversity. However, we could not do in-depth analysis on the effect of environmental factors to the differences among habitats due to lack of data such as on soil physical and chemical properties by sample plot. Several studies utilized the results of correlation analysis between earthworm and soil indicators for further purpose, for example considering earthworm density or diversity as proxy for levels of soil quality or pollution (e.g., [73–75]). Some of these studies also assessed the level of metal pollution in the soils and in the earthworm tissue or casts [74,76]. Further researches on earthworm in the two study communes are therefore necessary and should focus on assessing the impacts of soil and terrestrial conditions on the earthworm density and diversity. The terrestrial condition includes type and abundance of plant litter biomass. The researches can also associate the results of soil-earthworm analysis with the potential role of earthworm in the two study communes on watershed service. *4.4. Further Studies for the Unidentified Species*  In this study, we could demonstrate that natural- and human-disturbed habitats have different levels of earthworm density and diversity. However, we could not do in-depth analysis on the effect of environmental factors to the differences among habitats due to lack of data such as on soil physical and chemical properties by sample plot. Several studies utilized the results of correlation analysis between earthworm and soil indicators for further purpose, for example considering earthworm density or diversity as proxy for levels of soil quality or pollution (e.g., [73–75]). Some of these studies also assessed the level of metal pollution in the soils and in the earthworm tissue or casts [74,76]. Further researches on earthworm in the two study communes are therefore necessary and should focus on assessing the impacts of soil and terrestrial conditions on the earthworm density and diversity. The terrestrial condition includes type and abundance of plant litter biomass. The researches can also associate the results of soil-earthworm analysis with the potential role of earthworm in the two study communes on watershed service.

#### In our study, seven out of nine unidentified species belong to *Amynthas* genus. The seven species were mostly found in NF and RF with only *Amynthas* sp. 1 was an *4.4. Further Studies for the Unidentified Species*

omnipresent inhabitant of all observed habitats. Also, *Amynthas* sp. 1 had a relatively high occurrence frequency in human-disturbed habitats such as AF and HG. Compared to the In our study, seven out of nine unidentified species belong to *Amynthas* genus. The seven species were mostly found in NF and RF with only *Amynthas* sp. 1 was an omnipresent inhabitant of all observed habitats. Also, *Amynthas* sp. 1 had a relatively high occurrence frequency in human-disturbed habitats such as AF and HG. Compared to the eleven identified species of *Amynthas* genus, three out of the eleven were omnipresent inhabitants, namely *Amynthas aspergillum*, *Amynthas divitopapillatus*, and *Amynthas modiglianii*, while the other eight were inhabitants of some habitats only, or one habitat such as the case for *Amynthas wui*. Therefore, earthworm species of *Amynthas* genus likely have different land colonization patterns or perhaps levels of resistance to the new habitat's condition.

The seven unidentified species have different external and internal morphological features compared to 115 earthworm species of *Amynthas* genus that have been identified from past studies across regions in Vietnam. To ascertain if all the seven species are new to the country, a comparison with four species of the genus *Amynthas* Kinberg, 1867 recently identified in the southeastern part of the country [77] is however necessary. Further studies can also compare morphological features of the unidentified species with those of *Amynthas* species identified in other countries. For example, earthworm species of the genus *Amynthas* Kinberg, 1867 of family Megascolecidae, the same family with *Amynthas* earthworms found in our study, had been known as common inhabitants of natural forests in the northern part of Laos [78]. Recently, three new earthworm species belonging to this genus were found in that part of the country [78]. Apart from morphological features, further studies should compare living habitat and land colonization pattern between the unidentified and known *Amynthas* species. Similar efforts are necessary for the unidentified species of other genus found in the current study, namely one each of *Metaphire* and *Polypheretima* genus of the Megascolecidae family.

#### **5. Conclusions**

The conversion of natural forests to different land uses in the two buffer-zone communes of Song Thanh Nature Reserve had reduced earthworm diversity which can substantially affect soil health and ecosystem functions in the two communes. Also, among the identified species, *P. corethrurus* was omnipresent and largely dominated human-disturbed habitats. The dominance of this species can bring unfavorable impacts such as soil compaction and continuous threat to native earthworm species. This will further affect soil porosity and related ecosystem functions and reduce the role of the two communes to the nature reserve and watershed of which the two communes are part of. The area of human-disturbed habitats in the two communes potentially expand driven by augmenting population and livelihood pressures. This will lead to further loss of earthworm diversity and stronger dominance of *P. corethrurus*. To avoid the further decline of earthworm diversity and rampant expansion of *P. corethrurus* in the two study communes, protection of the remaining natural and regenerated forests is urgent, and tree-based farming systems such as agroforestry should be promoted to reconcile earthworm conservation and local livelihoods. The results of this study enrich current limited knowledge on impacts of forest and land use conversion on earthworm density and diversity in Vietnam. Furthermore, the results can generate stronger concerns with respect to the degradation of belowground diversity in the country and call attention to the urgent need for strengthened forest protection efforts.

**Author Contributions:** Conceptualization, R.M. and M.v.N.; data curation, R.M., S.V.H., V.M.D., N.B.T.D., A.D.N., D.H.L., D.T.T.H., and M.v.N.; formal analysis, S.V.H., V.M.D., N.B.T.D., A.D.N., D.H.L., and D.T.T.H.; investigation, S.V.H., V.M.D., N.B.T.D., A.D.N., D.H.L., and D.T.T.H.; methodology, S.V.H., V.M.D., N.B.T.D., A.D.N., D.H.L., D.T.T.H., and M.v.N.; supervision, M.v.N.; writing original draft, R.M., S.V.H., V.M.D., N.B.T.D., A.D.N., D.H.L., and D.T.T.H.; writing—review & editing, R.M., S.V.H., and M.v.N. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the United States Agency for International Development (USAID) to assess the livelihoods and ecological conditions of buffer-zone communes in Song Thanh Nature Reserve, Quang Nam province, and Phong Dien Nature Reserve, Thua Thien Hue province, grant number TF069018. The APC was funded by the same project.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data available in a publicly accessible repository.

**Acknowledgments:** We thank Le Tho Son, Nguyen Minh Quang, and Dao Van Chuong from the Vietnam National University of Forestry for their great contribution in the field inventory; Delia C. Catacutan, Nguyen Quang Tan, Do Trong Hoan, Nguyen Mai Phuong, Kurniatun Hairiah, and Pham

Thanh Van from the World Agroforestry (ICRAF) and Do Hoang Chung from Thai Nguyen University of Agriculture and Forestry, for their great help in field inventory and valuable inputs during data analysis and manuscript writing. We also thank valuable comments from three anonymous reviewers.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

### **Appendix A**

**Table A1.** Land cover types and areas in the study communes in 2016.


Data source: [16].



\* PM: Phuoc My commune, TB: Ta Bhing commune.



**Table A4.** Occurrence frequency of earthworm species by habitat (codes as in Table A2).



**Table A4.** *Cont.*

\* Number of habitats where the species was found.

**Table A5.** Quantity (n') and biomass (b') dominance (%) of earthworm species by habitat and soil depth (NF: natural forest, RF: regenerated forest, GL: grasslands, PF: planted forest, AF: agroforestry, UC: upland annual crops, HG: home gardens).


