Soil Physical Quality and Relationship to Changes in Termite Community in Northwestern Colombian Amazon
Ervin Humprey Duran-Bautista1,2* 
Yolima Muñoz Chilatra1 
Juan Diego Galindo1 
Tania Alexandra Ortiz1 
María Fernanda Bermúdez3,4*- 1Programa de Ingeniería Agroecológica, Universidad de la Amazonia, Florencia, Colombia
- 2Centro de Investigaciones Amazónicas CIMAZ Macagual Cesar Augusto Estrada González, Florencia, Colombia
- 3Laboratorio de Entomología, Grupo de Investigación en Entomología, Universidad de la Amazonia, Florencia, Colombia
- 4Maestría en Ciencias–Entomología, Universidad Nacional de Colombia–Sede Medellín, Medellín, Colombia
Conversion from Amazon forest to low-management pasture or agriculture causes not only degradation of aboveground vegetation but also negative changes in soil properties and ecosystem services. This study aimed to evaluate the impact of physical soil degradation on termite community changes in three contrasting land uses (natural regeneration, rubber plantations, and silvopastoral systems). Soil physical quality was assessed through a set of physical variables, such as bulk density, porosity, soil macro-aggregation state, Visual Evaluation of Soil Structure (VESS) and penetration resistance, which were summarized in an overall synthetic indicator of physical quality. Besides, transects of 20 × 2 m were established in each land use; each transect was divided into four sections of 5 m to search and collect termites during 1 hour in each section; likewise, termites were collected from blocks of soil 25 × 25 × 10 cm (length, width, and depth, respectively) adapted from the Tropical Soil Biology and Fertility (TSBF) method. In total, 60 transects were evaluated, 20 in each land use. A total of 41 species were collected across the three land uses evaluated: natural regeneration presented 60% of the collected species (25 species), silvopastoral systems 53% (22 species), and rubber plantations 39% (16 species). Additionally, composition species from the silvopastoral, agroforestry systems, and natural regeneration were different, and a close association between these last land uses was observed. Soil physical characteristics showed significant variations between land uses. The rubber plantations presented lowest values of soil physical quality, while the natural regeneration showed high soil physical quality. These changes affected termite community and lead to changes in its composition with disproportionate loss of some species; however, there are some that can acclimate well to the decline in the soil physical quality.
Introduction
The Amazon region plays an essential role to the Earth, harboring about 10–15% of land biodiversity (Hubbell et al., 2008; Nobre et al., 2016). It serves as a sink for greenhouse gasses by absorbing up to 5% of annual CO2 emissions (~2 billion tons) (Saatchi et al., 2011), and it has been well-recognized as one of the key components of the Earth's climate system (Malhi et al., 2008). Despite its relevance, this ecosystem is under increasing human pressures related to agriculture and livestock farming, timber extraction, mining, illicit crops, among others, which resulted in unprecedented rates of land cover changes due to forest clearing, degradation, and fragmentation (Sonter et al., 2017; Roque et al., 2018). In Colombia, the Amazon covers 42% of the country's territory (Etter et al., 2006). Despite many efforts to enhance forest landscape conservation and its biodiversity (Furumo, 2020), this region concentrated 70.1% (1,381.8 km2 year−1) of deforestation from a national total in 2018. Deforestation in the department of Caquetá alone accounted for 23.7% (467.6 km2 year−1), becoming an important hotspot of deforestation in Colombian Amazon (SMByC 2020).
Forest clearing not only generates the loss of biological diversity at all levels but also leaves the soil exposed to processes of erosion and degradation. For example, in deforested areas for crops and degraded pastures by domestic livestock, a reduced soil moisture, increase in the bulk density and a decrease in total porosity at soil surface has been observed due to animal grazing, compaction of agricultural machinery, and the impact of the raindrops (Cherubin et al., 2016; de Queiroz et al., 2020).
In tropical rainforests, termites constitute an important part of the soil fauna biomass (Bignell and Eggleton, 2000). Some studies suggest that termites may represent 40–65% of overall soil macrofaunal biomass in some biotopes (Loveridge and Moe, 2004). They contribute to ground decomposition of litter, plant growth, and overall biodiversity (Jouquet et al., 2011; Bottinelli et al., 2015).
Specifically, in tropical and subtropical soils, the role of termites has been widely documented demonstrating its impact on soil structure and properties, such as soil formation and aeration, bioturbation, water content and hydraulic proprieties, soil organic matter, and nutrient cycling (Bignell, 2006; Jouquet et al., 2011; Harit et al., 2017). These alterations are caused by the development of their biogenic structures (nests, foraging tunnels, and formation of galleries) and feeding strategies (ingesting and egesting soil) (Jouquet et al., 2016a).
However, the effect generated by changes in soil characteristics on termite's community distribution has been poorly studied, and it is limited to some efforts aimed to analyze the influence of soil chemical properties on termites (Bourguignon et al., 2015; Jouquet et al., 2015). In this regard, we evaluate the impact of soil physical quality on changes in termite communities and hypothesized that soil physical quality affects different termite communities.
Materials and Methods
Study Area
The study area was located in the Department of Caquetá southeastern Colombia between 1°30′N and 1°10′N, and between 75°35′W and 76°0′W (Figure 1). It has an average yearly temperature of 25°C, yearly rainfall of 3,600 mm, concentrated between April and November and a light dry season from December to March.
Figure 1. Location of sampling sites in Caquetá department, Colombia.
Sampling in this area was carried out during September 2019 in the following land uses: natural forest (NF), rubber plantations (RP), degraded pasture (DP) according to the methodology Corine Land Cover adapted by the Geographic Institute Agustín Codazzi (IGAC) and the Institute of Hydrology, Meteorology and Environmental Studies (IDEAM) for Colombia (Table 1). A brief characterization of the chemical properties of the soil in the evaluated areas is presented in Table 2.
Table 1. Description of land use types evaluated; the photos are representative of the sampled area.
Table 2. Chemical characteristics of the soil in the evaluated areas.
Termite Sampling
Termites were collected from transects of 20 × 2 m (adapted from Swift and Bignell, 2001) located in the central part of each land use type avoiding the edge effects; a minimum distance of 30 m between transects was guaranteed. In the silvopastoral system, the transects were located under the Gliricidia sepium tree line, while in the agroforestry system, they were established under the Hevea brasiliensis and Theobroma cacao intercalated line; a total of 60 transects were evaluated, 20 in each land use.
The transects were divided into four sections of 5 m. In each section, termites were searched for in any specific microhabitat, such as nests, dead wood, trunks, roots, and under rocks, 1 hour per person. When a colony was encountered, collection time could not exceed 3 min to avoid overestimation of the dominant species. Also, termites were collected from extracted blocks of soil by using a 25 × 25 × 10 cm (length, width, and depth, respectively) metallic frame adapted from the International Organization for Standardization (ISO) and the Tropical Soil Biology and Fertility (TSBF) methods (Anderson and Ingram, 1993; ISO, 2011).
Identification was made at genus and species levels, using taxonomic guides and reviews (Constantino, 2002a; Constantino et al., 2006; Constantino and Carvalho, 2012; Rocha et al., 2012; Constantini and Cancello, 2016). When identification was not possible, individuals were separated into morphospecies according to morphological differences. Individuals of the Apicotermitinae subfamily were separated at the species and morphospecies levels by dissecting the enteric valve and comparing its morphology with that described in the literature (Bourguignon et al., 2010, 2013, 2016). Termite samples are kept in 2-ml vials containing 80% ethanol, in the Colección Entomológica, Laboratorio de Entomología, Universidad de la Amazonia (LEUA).
Soil Physical Parameters Determination
A set of five physical variables of the soil was evaluated, which was considered key to have a general idea of the physical quality of the soil namely:
• Bulk density
• Soil moisture
• Soil resistance to penetration
• Soil macro-aggregation state
• Visual evaluation of soil structure
Undisturbed soil samples were collected in the center of the 0–10, 10–20, and 20–30 cm layers using a metallic ring (5 × 5 cm 98 cm2) and disturbed samples from the same soil depths. In the laboratory, undisturbed samples were weighed, dried in a forced-air oven at 105°C for 48 hours, and weighed again. Bulk density (BD, Mg m−3) was calculated by dividing the soil dry mass by the volume of the cylinder, whereas soil moisture (%) was determined by the equation: soil moisture = [(dry soil mass/wet soil mass) – 1] × 100.
Measurements of Soil Resistance to Penetration (SRP) were performed using a hand penetrometer (Eijkelkamp) around the soil sampling trenches down to 30 cm with angle and surface area of a cone of 60° and 3.3 cm2, respectively. Also, total porosity was calculated by using bulk density and particle density of soil through Equation (1):
where TP is the total porosity, Bd and Pd are the bulk density and particle density, respectively.
The state of soil macro-aggregation was determined using the methodology described by Velasquez et al. (2007). In each transect, we also collected a small block of soil (10 × 10 × 10 cm) and separated macroaggregates based on the visual separation of macroaggregates (>4 mm) according to biogenic aggregates produced by macroinvertebrates, or root aggregates made of soil stuck to the roots, physical aggregates produced by physical processes, and non-macroaggregated soil. Separation was done by gently breaking the soil apart into its natural constituents.
The Visual Evaluation of Soil Structure (VESS) was performed following the methodology proposed by Ball et al. (2017). In each transect, a soil block of 20 × 10 × 25 cm (width, thickness, and depth, respectively) was extracted, and each block was divided into two layers: topsoil (0–10 cm) and subsoil (10–25 cm). This method involves the removal and gentle breakup of a spadeful of topsoil by hand to reveal the main structural units and any layers of contrasting aggregation. Each layer is compared to the photographs with identified dimensions and descriptions in a colored chart and allocated to one of five soil quality scoring categories (Sq) from Sq 1 = best to Sq 5 = worst topsoil quality (VESS) and Ssq 1 = best to Ssq 5 = worst subsoil quality (SubVESS).
An overall weighted Sq score was calculated for each sample based on the individual score and thickness of each contrasting soil layer, according to Equation (2)
where VESS Sqscore is the overall VESS score of the sample, Sqi and Ti are the score and thickness of each identified soil layer, respectively, and TT is the total thickness of the soil sample.
Data Analysis
Termite Community Analysis
The relative abundance of each termite species was expressed as the number of encounters established by recording the presence of each species only once in each plot. Encounters and richness species data per transect were analyzed using a mixed generalized linear model (GLMM) with negative binomial and Poisson distribution for abundance and richness, respectively, stating land use as a fixed effect and farms as a random effect. We use the “lme4” package (Bates et al., 2018). Fisher multiple comparison tests (α = 0.05) were also applied for fixed effect using the “multcomp” package (Hothorn et al., 2008). To compare alpha- and beta-diversity across sites, the Shannon exponential (eH′) (Jost, 2006), and Simpson's inversed (1-D) indices, as well as rarefaction and extrapolation curves of richness were calculated using iNext (Hsieh et al., 2020). To test the hypothesis of differences in Shannon exponential and Simpson's inversed between land uses, we performed a generalized linear model (GLM). Non-metric multidimensional scaling (NMDS) based on the Bray–Curtis index of dissimilarity was used to compare composition termite species between land uses; species with a single encounter in all sampling were not considered to avoid the zero-inflated effect. Besides, an analysis of similarities (ANOSIM) was performed with 999 permutations to test statistically whether there were significant differences between land uses. A similarity percentage analysis (SIMPER) was used to identify which species contributed the most to the differences in termite community structure between land use types. These analyses were performed with “vegan” package (Oksanen et al., 2019).
Analysis Between Physical Soil Variables and Land Use Groups
Physical soil variables were analyzed with a generalized linear model (GLM) and a principal component analysis (PCA) associated with a Monte Carlo test from the “ade4” package to assess differences among land use and soil physical characteristics, and their statistical significance (Dray and Dufour, 2007). Besides, physical soil variable dataset was summarized in overall synthetic indicator of physical quality, following the methodology adapted from Velasquez et al. (2007). A principal component analysis (PCA) was used to determine the parameters that best-captured variance within a dataset. We selected those variables with a significant contribution (>50% of the explained variance) in either of the first two main axes. The selected variables were then combined into a single value indicator and scaled to a number ranging from 0.1 to 1.0 using a homothetic transformation.
Analysis Between Physical Soil Variables and Termite Species
To examine the relationship of termite species abundance to physical soil variables, canonical correspondence analysis (CCA) was performed. This analysis was run with two matrices: a dependent matrix containing 24 termite species whose abundance was >1% of the total number of individuals collected, species with abundance lower than 1% were grouped in the category “others” and considered in the analysis. The independent matrix contained the soil physical variables, which we did not keep in the analysis variables that had variance inflation factor (VIF) values higher than 10. The significance of these associations was verified using a permutation test from the “vegan” package (Oksanen et al., 2019); neither were they considered variables that were not significant using permutation tests. All analyses were performed using the R language software, version 3.5.3 (R Core Team, 2018).
Results
Species Richness and Community Structure
A total of 41 species from seven subfamilies and 24 genera were collected across the three land uses evaluated (Table 3) in 135 encounters. Likewise, Nasutitermitinae and Syntermitinae were the most abundant subfamilies with 43 and 24% of the total species collected, respectively. In contrast, Amitermitinae, Coptotermitinae, and Heterotermitinae recorded only one species each in the sampling. The genera with the largest number of species were Nasutitermes, Anoplotermes, and Syntermes containing 10 and 3 morphospecies, respectively, representing 38% of the total specific richness.
Table 3. Termite species encounters per transect (or morpho-species) at individual sampling points, in Natural Regeneration (RN), silvopastoral systems (SPS), rubber plantations (RP) in Colombian Amazon region.
The most frequent species in natural regeneration were Silvestritermes cf. lanei (Canter, 1968; 10 encounters), Nasutitermes guayanae (Holmgren, 1910; 5 encounters), and Anoplotermes sp. 1 (5 encounters). In the silvopastoril system, the most frequent species were Aparatermes silvestrii (Emerson, 1925; seven encounters), Nasutitermes sp. 7 (six encounters), and Embiratermes neotenicus (Holmgren, 1906; four encounters). Finally, in the rubber plantations, the more frequent species were Nasutitermes guayanae (six encounters), Silvestritermes cf. lanei (six encounters), and Heterotermes tenuis (Hagen, 1858; four encounters).
There were no statistically significant differences according to GLMM for species richness by transect among land uses; neither were there statistically significant differences according to GLM in the Shannon exponential (H) and Simpson's inverse (1-D) diversity indices per transect between land uses. When analyzing the accumulated species richness in each land use, natural regeneration presented 60% of the collected species (25 species), silvopastoral systems 53% (22 species), and rubber plantations 39% (16 species). Furthermore, sampling coverage curve including interpolation and extrapolation values for each land use averaged 72% of completeness and suggests that the sampling appropriately represent the termite diversity in the study area (Table 4).
Table 4. Means and standard error for number of species and encounters per transect, Shannon and Simpson indices, as well as total number of species and genus, sample coverage from estimated species richness for Natural Regeneration (RN), silvopastoral systems (SPS), and rubber plantations (RP) in Colombian Amazon region.
The NMDS diagram based on the Bray–Curtis index of dissimilarity showed termite species composition differences (Figure 2). In general, 12 species were found exclusively from natural regeneration, nine from silvopastoral systems, and five in the rubber plantation. Eight termite species were found in all land uses. The ANOSIM demonstrated a significant statistical difference between the land uses evaluated (Global R = 0.29, P = 0.002). The species from the silvopastoral system was different from the species in agroforestry systems and natural regeneration. Moreover, a close association between these last land uses was suggested by the analysis (Figure 2).
Figure 2. Nonmetric multidimensional scaling (NMDS) composition of termite species collected in natural regeneration (RN), silvopastoral systems (SPS), and rubber plantations (RP) in Colombian Amazon region. This ordination analysis was calculated from the Bray–Curtis dissimilarity index's (Stress = 0.07). Ano1, Anoplotermes sp. 1; Apa_sil, Aparatermes silvestrii; Ara_zet, Araujotermes zeteki; Cyl_cap, Cylindrotermes capixaba; Cyr_ang, Cyrilliotermes angulariceps; Emb_neo, Embiratermes neotenicus; Het_ten, Heterotermes tenuis; Micro_arb, Microcerotermes arboreus; Nasu_gua, Nasutitermes guayanae; Nasu1, Nasutitermes sp. 1; Nasu3, Nasutitermes sp. 3; Nasu4, Nasutitermes sp. 4; Nasu6, Nasutitermes sp. 6; Nasu7, Nasutitermes sp. 7; Nasu8, Nasutitermes sp. 8; Nasu9, Nasutitermes sp. 9; Sil_la, Silvestritermes cf. lanei; Spi_lon, Spinitermes longiceps; Sub1, Subulitermes sp. 1; Ter_hisp, Termes hispaniolae.
The similarity percentage analysis (SIMPER) allowed us to identify which termite species mainly contributed to the dissimilarity in the species composition between land uses. The SIMPER showed that species of the genus Nasutitermes, the xylophages species Heterotermes tenuis, Microcerotermes arboreus, and Nasutitermes guayanae and Aparatermes silvestrii were the main species responsible for the differences between land uses (Table 5). Nasutitermes was the more common species in natural regeneration, with a 12.2% contribution; while Heterotermes tenuis, Microcerotermes arboreus, and Nasutitermes guayanae were the more abundant species in rubber plantation, with a 21.3% contribution, and Aparatermes silvestrii dominant in silvopastoral systems, with a 13.25% contribution.
Table 5. Similarity-percentage analysis (SIMPER) comparing species composition of termites collected in natural regeneration (RN), silvopastoral systems (SPS), and rubber plantations (RP) in Colombian Amazon region.
Physical Soil Variables
Physical soil characteristics showed significant variations between land uses according to GLM. Rubber plantations presented lowest values of the overall synthetic indicator of physical quality (0.41 ± 0.04; Figure 3A) associated with highest values for the parameters of bulk density (1.24 ± 0.06 Mg m−3; Figure 3B), soil resistance to penetration (1.17 ± 0.09 Mpa; Figure 3C), and visual evaluation of soil structure (1.49 ± 0.16 Sq score; Figure 3F). Meanwhile, natural regeneration showed high physical soil quality with highest values for soil porosity (57.16 ± 2.69%; Figure 3D) and moisture (44.66 ± 2.28%; Figure 3E).
Figure 3. Physical soil characteristics: (A) overall synthetic indicator of physical quality, (B) bulk density, (C) soil resistance to penetration, (D) total porosity, (E) soil moisture content, and (F) visual evaluation of soil structure in Natural Regeneration (RN), silvopastoral systems (SPS), and rubber plantations (RP) in Colombian Amazon region. The bars represent standard deviation, and the different letters indicate statistically significant differences according to the Fisher Test (P < 0.05).
The variance of physical soil properties was explained by 21.24% by land use according to the PCA. Axes 1 and 2 of the PCA represent 36.2 and 19.2%, respectively. Axis 1 separated on one side the sites with high porosity, macro-aggregation status, and soil moisture (Figure 4A) to which secondary vegetation was associated, at the other end of Axis 1 are the sites with high values for the parameters of bulk density and soil resistance to penetration like rubber plantations (Figure 4B). The formation of Axis 2 was mainly contributed on one side by sites with high remaining material from aggregate morphology and on the other side by high Sq score visual evaluation of the soil structure.
Figure 4. Projection of physical soil characteristics and the 60 sampling points according to land uses (A), and in the plane formed by Axes 1 and 2 of PCA (B) in Colombian Amazon region. Ellipses show sampling points grouped from a given land use centered on their respective bari-centers (P < 0.01: Monte Carlo permutation test, Explained variance: 21.24%). BA, biogenic aggregates; BD, bulk density; NA, soil not macro-aggregate; PA, physical aggregates; RA, root aggregates; RM, remaining material; SRP, soil resistance to penetration; TP, total porosity; VESS, visual evaluation of soil structure.
Soil macro-aggregation on average, 21.5% of the soil volume was not macro-aggregated, biogenic aggregates contributing 21.8% of the total soil mass and physical aggregates comprising 20.9% (Figure 5). Important differences between land uses were suggested in soil macro-aggregation, with high fraction from remaining material representing 33.4% of soil under rubber plantations, while natural regeneration possessed much larger percentages of biogenic aggregates (28.1%).
Figure 5. Proportions of land aggregates classified according to their origin in Natural Regeneration (RN), silvopastoral systems (SPS), and rubber plantations (RP) in Colombian Amazon region. BA, biogenic aggregates; NA, soil not macro-aggregate; PA, physical aggregates; RA, root aggregates; RM, remaining material.
Effect of Physical Soil Variables on Termite Species Abundance
The CCA explained 16.1% of the variation in termite species (P < 0.001) and included variables related to soil physical variables (bulk density, soil resistance to penetration, total porosity, visual evaluation of soil structure, and soil moisture). The first two axes of the biplot accounted for 48.9 and 45.0%, respectively, of the variance explained by the CCA (Table 6).
Table 6. Canonical correspondence analysis (CCA) permutation test results examining the association between termite species composition with physical soil variables.
Some termite species were associated with the physical soil variables, for example, Embiratermes neotenicus associated with high soil moisture, Curvitermes odontognathus and Cyrilliotermes angulariceps with high total porosity, a species group made up of the species Cylindrotermes capixaba, Microcerotermes arboreus, and Termes hispaniolae was related to high values of visual evaluation of soil structure, finally soil resistance to penetration was correlated positively with some species of the genus Nasutitermes and negatively with total species richness (Figure 6).
Figure 6. Canonical correspondence analyses (CCA) ordination plot for termite species composition and physical soil variables with a significant effect in Colombian Amazon region. The analysis deals with 24 termite species related to five physical soil variables. SRP, soil resistance to penetration; TP, total porosity; VESS, visual evaluation of soil structure; Ano1, Anoplotermes sp. 1; Apa_sil, Aparatermes silvestrii; Ara_zet, Araujotermes zeteki; Atli_ster, Atlantitermes stercophilus; Coa_cle, Coatitermes clevelandi; Cur_odon, Curvitermes odontognathus; Cyl_cap, Cylindrotermes capixaba; Cyr_ang, Cyrilliotermes angulariceps; Emb_neo, Embiratermes neotenicus; Het_ten, Heterotermes tenuis; Micro_arb, Microcerotermes arboreus; Nasu_gua, Nasutitermes guayanae; Nasu1, Nasutitermes sp. 1; Nasu3, Nasutitermes sp. 3; Nasu4, Nasutitermes sp. 4; Nasu6, Nasutitermes sp. 6; Nasu7, Nasutitermes sp. 7; Nasu8, Nasutitermes sp. 8; Nasu9, Nasutitermes sp. 9; Rot_bra, Rotunditermes bragantinus; Sil_la, Silvestritermes cf. lanei; Sub1, Subulitermes sp. 1; Syn_mol, Syntermes molestus; Ter_hisp, Termes hispaniolae; Tr, total richness.
Discussion
Termite Community
Although the number of species per transect did not present significant differences, the total number of species did show important variations, the area of natural regeneration to be on the way of recovery of termite diversity with 25 species collected in total, if it is compared with previous studies carried out in the Amazon that report 52 species for forest (Ackerman et al., 2009). The species richness found in the rubber plantation in our study (16 species) was greater compared to others' studies in Colombia for commercial rubber plantation (10 species; Pinzon et al., 2012).
Despite the present study was not evaluated, in the area of natural forest, the total diversity observed was similar to studies carried out by Pinzon et al. (2017) who reported the presence of 40 termite species in riparian forests of the eastern plains of Colombia. However, these data are not completely comparable because of the differences between land uses and the study areas; it can be used as a reference.
Several authors have reported a negative effect on termite assemblages according to the level of habitat disturbance and trampling, with a reduction in termite diversity and abundance (Barros et al., 2004; Vasconcellos et al., 2010). Nevertheless, the inclusion of trees in pastures and the isolation of these to prevent livestock from entering have benefited the diversity of termites in our study site (22 species) exceeding the number of species reported for degraded pastures (Carrijo et al., 2009; Duran-Bautista et al., 2020).
The changes in species composition observed in this study confirm the results of previous work that suggests that simplified habitats are associated with a decrease in termite richness (Vasconcellos et al., 2008; Luke et al., 2014). Besides, the exclusive presence and high dominance of species, such as Microcerotermes arboreus in the rubber plantations, suggest the creation of favorable environments that could turn them into pests (Constantino, 2002b). Likewise, in Colombia, the Hetereotermes spp. have been reported in high densities in rubber crops of Puerto López-Meta (Pinzon et al., 2012) and as a pest in rubber plantations in the Amazon region (Sterling et al., 2011).
The dominance of Aparatermes silvestrii in silvopastoral systems may be associated with selective low regarding the origin of organic matter they feed on, which may explain their presence and significant abundance in artificial pastures (Bourguignon et al., 2009). The dominance of some species of the genus Nasutitermes in the area of natural regeneration shows that this land use is in the process of recovery and is providing the necessary resources for the nesting and feeding of these species. Previous work has shown that Nasutitermes spp are frequent in areas with regeneration process for 5 years (Viana-Junior et al., 2014) as well as in reforestation sites with native species (de Paula et al., 2016).
Land Use Effect on Physical Soil Variables
In general, all land uses presented a good soil physical quality, and the natural regeneration presented the best values. In this sense, it is proposed that changes in forest land use to more intense land uses, such as agroforestry or silvicultural systems generate changes due to disturbance that alters the chemical and physical composition of the soil (Tellen and Yerima, 2018). The time of abandonment and non-intervention in the areas of natural regeneration seems to have generated a positive effect that lead to soil physical characteristic protection (Lisboa et al., 2009), reflected in high values of total porosity and soil moisture, low resistance to penetration, and consequently better overall quality of physical soil characteristics.
However, important differences were observed in rubber plantations with higher values of apparent density and resistance to penetration even surpassing silvopastoral systems. This could be explained by two mechanisms: (1) the silvopastoral systems evaluated were recently established, and (2) the samples were taken below the tree line where the cattle did not have access. de Souza Braz et al. (2013) propose that in areas where grazing pressure is low, it avoids the negative effects on physical soil characteristics generated by cattle trampling, such as increasing bulk density and soil resistance to penetration, reducing soil aeration, aggregate stability, and water infiltration that have been widely documented (Arevalo et al., 1998; Sharrow, 2007; Drewry et al., 2008).
Concerning the above, a study found that rubber agroforestry systems with low diversification, present higher values of apparent density and resistance to penetration than more diversified agroforestry systems due to high abundance of roots mixed with a soft layer composed primarily of porous and rounded aggregates present in these systems (Cherubin et al., 2018). Likewise, an assessment of the effect of different land uses on soil quality suggests that rubber trees are the plant cover causing the smallest change in terms of soil quality (Moreira et al., 2011).
Physical Soil Variables and Relationship With Termite Species
Soil compaction indicated by resistance to penetration had a negative relationship with total termite species richness. In this regard, Jones et al. (2003) suggested that compaction of soil surface and an increase in soil bulk density might have had a negative impact on soil-dwelling termites. This may be due to creation of tunnel and foraging galleries during termite's feeding and foraging activities by moving soil particles of different sizes (Jouquet et al., 2011, 2016b). These activities are difficult by increased soil compaction. Tucker et al. (2004) found that construction speed of tunnel network varied with soil compaction, being much slower in the more compacted soil (1.35 g/cm3). Although this study was conducted under laboratory conditions using a different substrate and moisture content, and it cannot be directly compared with the present study, however, clearly the same relationship was observed.
Nevertheless, some species of the genus Nasutitermes manage to support well in conditions of high resistance to penetration, and this would be related to the ability of these species to adapt to disturbed environments. For example, they can be pests in sugarcane crops (Constantino, 2002b; Miranda et al., 2004) where the physical quality of the soil is low by compaction due to intensive mechanization (Cherubin et al., 2016). The relationship of E. tenicus with high-humidity soil could be explained by two mechanisms: the first is the nesting habits developed for this specie as it builds an epigeal earthen nest of variable shape, sometimes at the bases of trees (Constantino, 1992a,b), and the second is one of the more frequently found termite species in primary and secondary forests (Bandeira et al., 2003; Davies et al., 2003; Souza et al., 2012) where soil moisture is higher.
Cyrilliotermes angulariceps was positively related to total porosity; the species of this genus are soil feeders (Constantino and Carvalho, 2012), and specifically, C. angulariceps feeds on soil particles richer in organic matter (Bourguignon et al., 2009) where soil porosity is generally higher. Although little information is available on the influence of termite foraging activity on soil porosity (Jouquet et al., 2016a), some authors have demonstrated that termite mediated increases in soil porosity (Lavelle, 1997). The relation of Cylindrotermes capixaba, Microcerotermes arboreus, and Termes hispaniolae with high values of soil structure can be explained because these species were presented in high density in the rubber plantations where this parameter presented the highest values. These results coincide with those from Thoumazeau et al. (2019), which reported a VESS gradient from intensive cash crops, the most disturbed system, through different ages of rubber planting to forest as the least disturbed system.
In general, changes in physical soil quality presented in different land uses affect termite communities and lead to changes in its composition with disproportionate loss of some species; however, there are some who can adapt well to the decline in physical soil quality, especially in rubber plantations and silvopastoral systems. Among the parameters evaluated, the resistance to penetration was the one that generated the greatest negative impact, and soil moisture was positively related to the presence of some species.
Data Availability Statement
The datasets presented in this article are not readily available because they were financed by the University of the Amazon. Requests to access the datasets should be directed to Ervin Humprey Duran-Bautista, ervinduranb@gmail.com.
Author Contributions
ED-B: conceiving and designing analyses, data collection, data analysis, and wrote paper. YM, JG, TO, and MB: data collection, data analysis, and wrote paper. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We were grateful to the Universidad de la Amazonia for funding this work as part of the project Provision of Soil Ecosystem Services by Termite Species in the Colombian Amazon Region. The project was led by members of the soil biology research seed. Collection permits for this study are under Resolution 01140 of September 30, 2016 of National Environmental License Authority by which a framework permit for the collection of specimens of wild species of biological diversity for non-commercial scientific research is granted to the Universidad de la Amazonia. We were also grateful to Yimi Katherine Angel for her help in drawing up the map of the study area and Xyomara Carretero for her help in improving the English writing. Finally, we would like to thank the reviewers for their suggestions to improve this document.
References
Ackerman, I. L., Constantino, R., Gauch, H. G. Jr., Lehmann, J., Riha, S. J., Fernandes, E. C. M., et al. (2009). Termite (insecta: isoptera) species composition in a primary rain forest and agroforests in central amazonia. Biotropica 41, 226–233. doi: 10.1111/j.1744-7429.2008.00479.x
Anderson, J. P. E., and Ingram, J. S. I. (1993). Tropical Soil Biology and Fertility. A Handbook of Methods. Wallingford: CAB International.
Arevalo, L. A., Alegre, J., Bandy, D. E., and Szott, L. T. (1998). The effect of cattle grazing on soil physical and chemical properties in a silvopastoral system in the Peruvian Amazon. Agrofor. Syst. 40, 109–124. doi: 10.1023/A:1006075114659
Ball, B. C., Guimarães, R. M. L., Cloy, J. M., Hargreaves, P. R., Shepherd, T. G., and McKenzie, B. M. (2017). Visual soil evaluation: a summary of some applications and potential developments for agriculture. Soil Tillage Res. 173, 114–124. doi: 10.1016/j.still.2016.07.006
Bandeira, A. G., Vasconcellos, A., Silva, M. P. da, and Constantino, R. (2003). Effects of habitat disturbance on the termite fauna in a highland humid forest in the caatinga domain, Brazil. Socibiology 42, 1–11.
Barros, E., Grimaldi, M., Sarrazin, M., Chauvel, A., Mitja, D., Desjardins, T., et al. (2004). Soil physical degradation and changes in macrofaunal communities in central Amazon. Appl. Soil Ecol. 26, 157–168. doi: 10.1016/j.apsoil.2003.10.012
Bates, D., Maechler, M., Bolker, B., Walker, S., Haubo, R., Christensen, B., et al. (2018). Package “lme4.” Available online at: https://cran.r-project.org/web/packages/lme4/lme4.pdf (accessed September 19, 2018).
Bignell, D. E. (2006). “Termites as soil engineers and soil processors,” in Intestinal Microorganisms of Termites and Other Invertebrates, eds H. König and A. Varma (Berlin; Heidelberg: Springer), 183–220. doi: 10.1007/3-540-28185-1_8
Bignell, D. E., and Eggleton, P. (2000). “Termites in ecosystems,” in Termites: Evolution, Sociality, Symbioses, Ecology, eds T. T Abe, D. E. Bignell, M. Higashi, T. Higashi, and Y. Abe (Dordrecht: Kluwer), 363–387. doi: 10.1007/978-94-017-3223-9_17
Bottinelli, N., Jouquet, P., Capowiez, Y., Podwojewski, P., Grimaldi, M., and Peng, X. (2015). Why is the influence of soil macrofauna on soil structure only considered by soil ecologists? Soil Tillage Res. 146, 118–124. doi: 10.1016/j.still.2014.01.007
Bourguignon, T., Drouet, T., Šobotník, J., Hanus, R., and Roisin, Y. (2015). Influence of soil properties on soldierless termite distribution. PLoS ONE 10:e0135341. doi: 10.1371/journal.pone.0135341
Bourguignon, T., Scheffrahn, R. H., Kreček, J., Nagy, Z. T., Sonet, G., and Roisin, Y. (2010). Towards a revision of the neotropical soldierless termites (isoptera:termitidae): redescription of the genus anoplotermes and description of longustitermes, gen. nov. Invertebr. Syst. 24, 357–370. doi: 10.1071/IS10012
Bourguignon, T., Scheffrahn, R. H., Nagy, Z. T., Sonet, G., Host, B., and Roisin, Y. (2016). Towards a revision of the neotropical soldierless termites (isoptera: termitidae): redescription of the genus grigiotermes mathews and description of five new genera. Zool. J. Linn. Soc. 176, 15–35. doi: 10.1111/zoj.12305
Bourguignon, T., Šobotník, J., Hanus, R., Krasulová, J., Vrkoslav, V., Cvačka, J., et al. (2013). Delineating species boundaries using an iterative taxonomic approach: the case of soldierless termites (isoptera, termitidae, apicotermitinae). Mol. Phylogenet. Evol. 69, 694–703. doi: 10.1016/j.ympev.2013.07.007
Bourguignon, T., Šobotník, J., Lepoint, G., Martin, J.-M., and Roisin, Y. (2009). Niche differentiation among neotropical soldierless soil-feeding termites revealed by stable isotope ratios. Soil Biol. Biochem. 41, 2038–2043. doi: 10.1016/j.soilbio.2009.07.005
Carrijo, T. F., Brandão, D., de Oliveira, D. E., Costa, D. A., and Santos, T. (2009). Effects of pasture implantation on the termite (Isoptera) fauna in the Central Brazilian Savanna (Cerrado). J. Insect Conserv. 13, 575–581. doi: 10.1007/s10841-008-9205-y
Cherubin, M. R., Chavarro-Bermeo, J. P., and Silva-Olaya, A. M. (2018). Agroforestry systems improve soil physical quality in northwestern Colombian Amazon. Agrofor. Syst. 93, 1741–1753. doi: 10.1007/s10457-018-0282-y
Cherubin, M. R., Karlen, D. L., Franco, A. L. C., Tormena, C. A., Cerri, C. E. P., Davies, C. A., et al. (2016). Soil physical quality response to sugarcane expansion in Brazil. Geoderma 267, 156–168. doi: 10.1016/j.geoderma.2016.01.004
Constantini, J. P., and Cancello, E. M. (2016). A taxonomic revision of the neotropical termite genus rhynchotermes (isoptera, termitidae, syntermitinae). Zootaxa 4109, 501–522. doi: 10.11646/zootaxa.4109.5.1
Constantino, R. (1992a). Abundance and diversity of termites (insecta: isoptera) in two sites of primary rain forest in Brazilian amazonia. Biotropica 24, 420–430. doi: 10.2307/2388613
Constantino, R. (1992b). Notes on embiratermes fontes (isoptera, termitidae, nasutitermitinae), with descriptions of two new species from Amapá state, Brazil. Bol. Mus. Para. Emílio Goeldi Zool. 8:2.
Constantino, R. (2002a). An illustrated key to Neotropical termite genera (insecta: isoptera) based primarily on soldiers. Zootaxa 67:1. doi: 10.11646/zootaxa.67.1.1
Constantino, R. (2002b). The pest termites of South America: taxonomy, distribution and status. J. Appl. Entomol. 126, 355–365. doi: 10.1046/j.1439-0418.2002.00670.x
Constantino, R., Acioli, A. N. S., Schmidt, K., Cuezzo, C., Carvalho, S. H. C., and Vasconcellos, A. (2006). A taxonomic revision of the neotropical termite genera labiotermes holmgren and paracornitermes emerson (isoptera: termitidae: nasutitermitinae). Zootaxa 1340, 1–44. doi: 10.11646/zootaxa.1340.1.1
Constantino, R., and Carvalho, S. H. C. (2012). A taxonomic revision of the neotropical termite genus cyrilliotermes fontes (isoptera, termitidae, syntermitinae). Zootaxa 3186, 25–41. doi: 10.11646/zootaxa.3186.1.2
Davies, R. G., Hernández, L. M., Eggleton, P., Didham, R. K., Fagan, L. L., and Winchester, N. N. (2003). Environmental and spatial influences upon species composition of a termite assemblage across neotropical forest islands. J. Trop. Ecol. 19, 509–524. doi: 10.1017/S0266467403003560
de Paula, R. C., de Moraes Lima Silveira, R., da Rocha, M. M., and Izzo, T. J. (2016). The restoration of termite diversity in different reforestated forests. Agrofor. Syst. 90, 395–404. doi: 10.1007/s10457-015-9862-2
de Queiroz, M. G., da Silva, T. G. F., Zolnier, S., da Rosa Ferraz Jardim, A. M., de Souza, C. A. A., de Souza, L. S. B., et al. (2020). Spatial and temporal dynamics of soil moisture for surfaces with a change in land use in the semi-arid region of Brazil. Catena 188:104457. doi: 10.1016/j.catena.2020.104457
de Souza Braz, A. M., Fernandes, A. R., and Alleoni, L. R. F. (2013). Soil attributes after the conversion from forest to pasture in Amazon. L. Degrad. Dev. 24, 33–38. doi: 10.1002/ldr.1100
Dray, S., and Dufour, A.-B. (2007). The ade4 package: implementing the duality diagram for ecologists. J. Stat. Softw. 22, 1–20. doi: 10.18637/jss.v022.i04
Drewry, J. J., Cameron, K. C., and Buchan, G. D. (2008). Pasture yield and soil physical property responses to soil compaction from treading and grazing—a review. Soil Res. 46, 237–256. doi: 10.1071/SR07125
Duran-Bautista, E. H., Armbrecht, I., Serrão Acioli, A. N., Suárez, J. C., Romero, M., Quintero, M., et al. (2020). Termites as indicators of soil ecosystem services in transformed amazon landscapes. Ecol. Indic. 117:106550. doi: 10.1016/j.ecolind.2020.106550
Etter, A., McAlpine, C., Wilson, K., Phinn, S., and Possingham, H. (2006). Regional patterns of agricultural land use and deforestation in Colombia. Agric. Ecosyst. Environ. 114, 369–386. doi: 10.1016/j.agee.2005.11.013
Furumo, P. R. (2020). Scaling up zero-deforestation initiatives through public-private partnerships: a look inside post-conflict Colombia. Glob. Environ. Change 62:102055. doi: 10.1016/j.gloenvcha.2020.102055
Harit, A., Moger, H., Duprey, J. L., Gajalakshmi, S., Abbasi, S. A., Subramanian, S., et al. (2017). Termites can have greater influence on soil properties through the construction of soil sheetings than the production of above-ground mounds. Insectes Soc. 64, 247–253. doi: 10.1007/s00040-017-0541-3
Hothorn, T., Bretz, F., and Westfall, P. (2008). Simultaneous inference in general parametric models. Biometric. J. 50, 346–363. doi: 10.1002/bimj.200810425
Hsieh, T. C., Ma, K. H., and Chao, A. (2020). Package “iNEXT” Title Interpolation and Extrapolation for Species Diversity. Hsinchu.
Hubbell, S. P., He, F., Condit, R., Borda-de-Água, L., Kellner, J., and Ter Steege, H. (2008). How many tree species are there in the Amazon and how many of them will go extinct? Proc. Natl. Acad. Sci. U.S.A. 105, 11498–11504. doi: 10.1073/pnas.0801915105
ISO (2011). ISO 23611-5:2011–Soil Quality–Sampling of Soil Invertebrates– Sampling of Soil Invertebrates–Part 5: Sampling and Extraction of Soil Macro-Invertebrates. Geneva.
Jones, D. T., Susilo, F. X., Bignell, D. E., Hardiwinoto, S., Gillison, A. N., and Eggleton, P. (2003). Termite assemblage collapse along a land-use intensification gradient in lowland central Sumatra, Indonesia. J. Appl. Ecol. 40, 380–391. doi: 10.1046/j.1365-2664.2003.00794.x
Jouquet, P., Bottinelli, N., Shanbhag, R. R., Bourguignon, T., Traoré, S., and Abbasi, S. A. (2016a). Termites: the neglected soil engineers of tropical soils. Soil Sci. 181, 157–165. doi: 10.1097/SS.0000000000000119
Jouquet, P., Chintakunta, S., Bottinelli, N., Subramanian, S., and Caner, L. (2016b). The influence of fungus-growing termites on soil macro and micro-aggregates stability varies with soil type. Appl. Soil Ecol. 101, 117–123. doi: 10.1016/j.apsoil.2016.02.001
Jouquet, P., Guilleux, N., Caner, L., Chintakunta, S., Ameline, M., and Shanbhag, R. R. (2015). Influence of soil pedological properties on termite mound stability. Geoderma 262, 45–51. doi: 10.1016/j.geoderma.2015.08.020
Jouquet, P., Traoré, S., Choosai, C., Hartmann, C., and Bignell, D. (2011). Influence of termites on ecosystem functioning. Ecosystem services provided by termites. Eur. J. Soil Biol. 47, 215–222. doi: 10.1016/j.ejsobi.2011.05.005
Lavelle, P. (1997). Faunal activities and soil processes: adaptive strategies that determine ecosystem function. Adv. Ecol. Res. 27, 93–132. doi: 10.1016/S0065-2504(08)60007-0
Lisboa, C. C., Conant, R. T., Haddix, M. L., Cerri, C. E. P., and Cerri, C. C. (2009). Soil carbon turnover measurement by physical fractionation at a forest-to-pasture chronosequence in the Brazilian Amazon. Ecosystems 12, 1212–1221. doi: 10.1007/s10021-009-9288-7
Loveridge, J. P., and Moe, S. R. (2004). Termitaria as browsing hotspots for African megaherbivores in miombo woodland. J. Trop. Ecol. 20, 337–343. doi: 10.1017/S0266467403001202
Luke, S. H., Fayle, T. M., Eggleton, P., Turner, E. C., Davies, R. G., Luke, S. H., et al. (2014). Functional structure of ant and termite assemblages in old growth forest, logged forest and oil palm plantation in Malaysian Borneo. Biodiv. Conserv. 23, 2817–2832. doi: 10.1007/s10531-014-0750-2
Malhi, Y., Roberts, J. T., Betts, R. A., Killeen, T. J., Li, W., and Nobre, C. A. (2008). Climate change, deforestation, and the fate of the Amazon. Science 319, 169–172. doi: 10.1126/science.1146961
Miranda, C. S., Vasconcellos, A., and Bandeira, A. G. (2004). Termites in sugar cane in Northeast Brazil: ecological aspects and pest status. Neotrop. Entomol. 33, 237–241. doi: 10.1590/S1519-566X2004000200015
Moreira, A., Fageria, N. K., and Garcia Y Garcia, A. (2011). Soil fertility, mineral nitrogen, and microbial biomass in upland soils of the central Amazon under different plant covers. Commun. Soil Sci. Plant Anal. 42, 694–705. doi: 10.1080/00103624.2011.550376
Nobre, C. A., Sampaio, G., Borma, L. S., Castilla-Rubio, J. C., Silva, J. S., and Cardoso, M. (2016). Land-use and climate change risks in the amazon and the need of a novel sustainable development paradigm. Proc. Natl. Acad. Sci. U.S.A. 113, 10759–10768. doi: 10.1073/pnas.1605516113
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., Mcglinn, D., et al. (2019). Package “vegan” Title Community Ecology Package Version 2.5-6. Available online at: https://cran.ism.ac.jp/web/packages/vegan/vegan.pdf (accessed October 17, 2020).
Pinzon, O. P., Baquero, L., and Beltran, M. A. (2017). Termite (isoptera) diversity in a gallery forest relict in the Colombian eastern plains. Sociobiology 64, 92–100. doi: 10.13102/sociobiology.v64i1.1184
Pinzon, O. P., Hernández, A. M., and Malagón, A. L. (2012). Diversidad de termitas (isoptera: termitidae, rhinotermitidae) en plantaciones de caucho en puerto lópez (meta, Colombia). Rev. Colomb. Entomol. 38, 291–298.
R Core Team (2018). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Rocha, M. M., Cancello, E. M., and Carrijo, T. F. (2012). Neotropical termites: revision of armitermes wasmann (isoptera, termitidae, syntermitinae) and phylogeny of the syntermitinae. Syst. Entomol. 37, 793–827. doi: 10.1111/j.1365-3113.2012.00645.x
Roque, F. D. O., Menezes, J. F. S., Northfield, T., Ochoa-Quintero, J. M., Campbell, M. J., and Laurance, W. F. (2018). Warning signals of biodiversity collapse across gradients of tropical forest loss. Sci. Rep. 8:1622. doi: 10.1038/s41598-018-19985-9
Saatchi, S. S., Harris, N. L., Brown, S., Lefsky, M., Mitchard, E. T. A., Salas, W., et al. (2011). Benchmark map of forest carbon stocks in tropical regions across three continents. Proc. Natl. Acad. Sci. U.S.A. 108, 9899–9904. doi: 10.1073/pnas.1019576108
Sharrow, S. H. (2007). Soil compaction by grazing livestock in silvopastures as evidenced by changes in soil physical properties. Agrofor. Syst. 71, 215–223. doi: 10.1007/s10457-007-9083-4
Sonter, L. J., Herrera, D., Barrett, D. J., Galford, G. L., Moran, C. J., and Soares-Filho, B. S. (2017). Mining drives extensive deforestation in the Brazilian Amazon. Nat. Commun. 8:1013. doi: 10.1038/s41467-017-00557-w
Souza, H. B. A., Alves, W. F., and Vasconcellos, A. (2012). Termite assemblages in five semideciduous Atlantic forest fragments in the northern coastland limit of the biome. Rev. Bras. Entomol. 56, 67–72. doi: 10.1590/S0085-56262012005000013
Sterling, A., Gómez, C. A., and Campo, A. A. (2011). Patogenicidad de metarhizium anisopliae (deuteromycota: hyphomycetes) sobre heterotermes tenuis (isoptera: rhinotermitidae) en hevea brasiliensis. Rev. Colomb. Entomol. 37, 36–42.
Swift, M., and Bignell, D. (2001). Standard Methods for Assessment of Soil Biodiversity and Land use Practice. Bogor: International Centre for Research in Agroforestry.
Tellen, V. A., and Yerima, B. P. K. (2018). Effects of land use change on soil physicochemical properties in selected areas in the North West region of Cameroon. Environ. Syst. Res. 7, 1–29. doi: 10.1186/s40068-018-0116-y
Thoumazeau, A., Bessou, C., Renevier, M. S., Panklang, P., Puttaso, P., Peerawat, M., et al. (2019). Biofunctool®: a new framework to assess the impact of land management on soil quality. Part B: investigating the impact of land management of rubber plantations on soil quality with the biofunctool® index. Ecol. Indic. 97, 429–437. doi: 10.1016/j.ecolind.2018.10.028
Tucker, C. L., Koehler, P. G., and Oi, F. M. (2004). Influence of soil compaction on tunnel network construction by the eastern subterranean termite (Isoptera: Rhinotermitidae). J. Econ. Entomol. 97, 89–94. doi: 10.1093/jee/97.1.89
Vasconcellos, A., Bandeira, A. G., Almeida, W. O., and Moura, F. M. S. (2008). Térmitas construtores de ninhos conspícuos em duas áreas de mata Atlântica com diferentes níveis de perturbação antrópica. Neotrop. Entomol. 37, 15–19. doi: 10.1590/S1519-566X2008000100003
Vasconcellos, A., Bandeira, A. G., Via, F., Moura, M. S., Farias, V., Araú Jo, P., et al. (2010). Termite assemblages in three habitats under different disturbance regimes in the semi-arid caatinga of NE Brazil. J. Arid Environ. 74, 298–302. doi: 10.1016/j.jaridenv.2009.07.007
Velasquez, E., Pelosi, C., Brunet, D., Grimaldi, M., Martins, M., Rendeiro, A. C., et al. (2007). This ped is my ped: visual separation and near infrared spectra allow determination of the origins of soil macroaggregates. Pedobiologia 51, 75–87. doi: 10.1016/j.pedobi.2007.01.002
Keywords: land use, soil compaction, soil moisture, soil-feeders termites, deforested Amazonia
Citation: Duran-Bautista EH, Muñoz Chilatra Y, Galindo JD, Ortiz TA and Bermúdez MF (2020) Soil Physical Quality and Relationship to Changes in Termite Community in Northwestern Colombian Amazon. Front. Ecol. Evol. 8:598134. doi: 10.3389/fevo.2020.598134
Received: 23 August 2020; Accepted: 09 November 2020;
Published: 18 December 2020.
Edited by:
Alberto Arab, Federal University of ABC, BrazilReviewed by:
Alexandre Vasconcellos, Federal University of Paraíba, BrazilRobin Casalla, Universidad del Norte, Colombia, Colombia
Copyright © 2020 Duran-Bautista, Muñoz Chilatra, Galindo, Ortiz and Bermúdez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ervin Humprey Duran-Bautista, ervinduranb@gmail.com; María Fernanda Bermúdez, mafezz_12@hotmail.com