+1 Recommend
2 collections
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Potential interactions between herbivorous arthropods and of their natural enemies on Caryocar brasiliense (Caryocaraceae) trees

      Read this article at

          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.


          ABSTRACT The study of arthropod communities foraging on monocultures of native plants (e.g. Caryocar brasiliense) in biodiversity hotspots (e.g. Brazilian Cerrado) is important to understand the processes driving their performance. This study evaluated the possible interactions between galling herbivores, free-feeding arthropods, and of their natural enemies on C. brasiliense trees, in an orchard (monoculture). The possible competition observed were: i) between galling insects with defoliators and phytophagous mites, ii) between sucking insects, iii) between phytophagous mites and, iv) between Aphis gossypii (Hemiptera: Aphididae) with Histiostoma sp. (Acari: Histiostomidae) on C. brasiliense trees. Numbers of predators and prey, parasitoids and hosts, and predators and parasitoids were directly related among them. The competitive interactions between herbivorous predator and parasitoid may reduce Eurytoma sp. (Hymenoptera: Eurytomidae) populations, a potential pest of C. brasiliense. The predator Zelus armillatus (Hemiptera: Reduviidae) and the parasitoids Sycophila sp. (Hymenoptera: Eurytomidae) and Ablerus magistretti (Hymenoptera: Aphelinidae) are important to determine the numbers of Eurytoma sp. and of its galls. These natural enemies may be important to control this galling insect on C. brasiliense trees. The arthropod competition affected their guild associated to C. brasiliense trees.

          Related collections

          Most cited references 67

          • Record: found
          • Abstract: not found
          • Article: not found

          Köppen's climate classification map for Brazil

            • Record: found
            • Abstract: found
            • Article: not found

            Land cover mapping of the tropical savanna region in Brazil.

            The Brazilian tropical savanna (Cerrado), encompassing more than 204 million hectares in the central part of the country, is the second richest biome in Brazil in terms of biodiversity and presents high land use pressure. The objective of this study was to map the land cover of the Cerrado biome based on the segmentation and visual interpretation of 170 Landsat Enhanced Thematic Mapper Plus satellite scenes acquired in 2002. The following land cover classes were discriminated: grasslands, shrublands, forestlands, croplands, pasturelands, reforestations, urban areas, and mining areas. The results showed that the remnant natural vegetation is still covering about 61% of the biome, however, on a highly asymmetrical basis. While natural physiognomies comprise 90% of the northern part of the biome, only 15% are left in its southern portions. Shrublands were the dominant natural land cover class, while pasturelands were the dominant land use class in the Cerrado biome. The final Cerrado's land cover map confirmed the intensive land use pressure in this unique biome. This paper also showed that Landsat-like sensors can provide feasible land cover maps of Cerrado, although ancillary data are required to help image interpretation.
              • Record: found
              • Abstract: found
              • Article: found
              Is Open Access

              Consequences of tropical land use for multitrophic biodiversity and ecosystem functioning

              The transformation from natural ecosystems to agricultural land use and its continued intensification has led to extensive losses in biodiversity and ecosystem services1 resulting in the degradation of human well being2. The transformation of lowland tropical rainforest to oil palm (Elaeis guineensis Jacq.) plantations has gained more recent attention as an especially severe threat to tropical biodiversity3 4. In the last 25 years the total plantation area of oil palm has tripled, with current global estimates of over 15 million hectares3, making this crop one of the world’s most rapidly expanding forms of agriculture5. It is now clear that the expansion of oil palm agriculture is one of the greatest causes of deforestation6 7, and this threat appears to be increasing without respite as Indonesia, one of the world’s leaders in oil palm, makes plans to double production by 2020 (ref. 8). The rapid expansion of such large-scale land-use transformation raises questions about the impending implications for biodiversity and ecosystem functioning in the tropics. Despite a broad consensus that biodiversity is positively correlated with ecosystem functioning in controlled experiments9 10, there are few real-world examples of such biodiversity–ecosystem functioning relationships11 12. In fact, until now there have been no studies that explore the relationship between biodiversity and ecosystem functioning in ecosystems undergoing agricultural land-use transformation to oil palm. Thus, our knowledge of this globally important land-use conversion is strongly limited. Furthermore, over the past decade there have been important advances towards multitrophic approaches in research investigating biodiversity–ecosystem functioning relationships9 13 14 15 16. Despite these advances, however, we are still substantially limited by the lack of clear approaches to quantify single measures of ecosystem functioning that can be compared among any combination of trophic levels. This has resulted in our inability to directly look at whole-community relationships between entire species assemblages and the respective functional processes carried out in these communities. Here we use the total energy flux between functional feeding guilds as a measure of multitrophic ecosystem functioning, as many studies have suggested process rates, such as energy fluxes, to be important proxies for ecosystem functioning10 13 17. Depending on the resource pool that the energy flux comes from, these fluxes can be directly related to ecosystem services such as decomposition18 19, plant biomass production20 21 or biocontrol through predation22. These energy flux calculations are based on metabolic scaling theory23 and principles of food-web energy dynamics18. Using individual metabolic rates that are dependent on body mass, environmental temperature and phylogenetic grouping18 24, combined with resource-specific assimilation efficiencies25 and energy loss to predation18, we present this energy flux calculation as a unified measure of multitrophic ecosystem functioning (Fig. 1). Studies that incorporate diversity across trophic levels to test the relationship between biodiversity and ecosystem functioning have predominantly used only biomass as the measure of ecosystem function26. However, the metabolic activity and thus the energy-processing rates of these biomass pools can vary substantially. Integrating over body mass, phylogeny and temperature with their constraints on metabolic rates, and additionally taking into account assimilation efficiencies and loss to predation, our measure of whole-community energy flux inherently incorporates not only biomass but also other important ecosystem attributes enabling the quantification of emergent functional properties of ecosystems that would otherwise remain undetected. As such, our measure of energy flux provides a comprehensive and robust measure of multitrophic ecosystem functioning that can be utilized for modelling biodiversity–ecosystem functioning relationships for any assemblage of taxonomic groups, while incorporating multiple ecological functions. In the tropical lowland rainforests of Sumatra, Indonesia, which have been undergoing vast land-use transformation to oil palm7, we quantify the impacts of this transformation ranging from tropical secondary rainforest, jungle rubber and intensively managed rubber, to oil palm. We utilize data gathered from 32 sites in Sumatra, Indonesia, comprising 2,415 populations of 871 species. First, we investigate the biodiversity value of jungle rubber, conventional rubber and secondary forest compared with oil palm agriculture by comparing observed species richness, density and biomass of litter-associated macroinvertebrate communities across these systems. Second, as a multitrophic measure of the rate of ecosystem processes carried out by these communities, we calculate total solid fresh mass energy flux in a system by incorporating community metabolism27, resource-specific assimilation efficiencies and biomass loss to predation18 into whole-community energy flux equations (Fig. 1). This provides a quantitative measure of multitrophic ecosystem functioning, defined here as the total flux of energy from any resource pool to consumer trophic levels. In addition, this measure can be attributed to specific functional feeding guilds within communities to look for patterns in ecosystem functioning at different trophic levels. Using the energy-mass flow conversion28, we express energy flux as kilograms per hectare, per year and explore the relationship between total species diversity and energy flux, distinguishing among four transformation systems to test for land-use-dependent biodiversity–ecosystem functioning relationships. Our results demonstrate strong losses in species diversity that, in turn, predict reductions in whole-community energy fluxes. However, these reductions are strongest in oil palm systems, suggesting that land-use conversion from forest to oil palm causes disproportionally strong losses in multitrophic ecosystem functioning. Results Transformation to oil palm leads to biodiversity loss Using generalized linear mixed effects models, we show that transformation of tropical rainforest to oil palm plantations leads to severe losses in species richness (45% decline), animal density (48% decline) and biomass (52% decline; Fig. 2a–c and Supplementary Table 1), supporting previous studies suggesting that land-use transformation to oil palm poses one of the greatest threats to global biodiversity3. Beyond mere diversity effects, land-use transformation altered animal densities and biomass, threatening to not only drive species extinctions but also to eliminate vital ecological functions. The effects of land-use transformation on species richness and animal densities were additionally dependent on functional feeding guilds, with predators decreasing in species richness and density most rapidly (Fig. 2a–c and Supplementary Table 1) as could be expected for higher trophic level feeding guilds29. Such alteration of higher trophic levels is likely to have severe indirect functional impacts on other functional guilds within the trophic network30. Community metabolism Summing up individual metabolic rates, we demonstrate that transformation of forest to oil palm yields a 51% decrease in community metabolism, with jungle rubber and rubber only 16% and 10% below forest levels of community metabolism, respectively. However, all systems yielded significantly higher community metabolism than oil palm (Fig. 2d and Supplementary Table 1). As such, we show that ecosystem energy processing is critically reduced in oil palm plantations. Interestingly, biomass responses to land-use transformation among feeding guilds were not clearly comparable to responses in community metabolism (Fig. 2c,d). This suggests that systematic changes in species composition, body-mass distributions (Supplementary Fig. 1) and biomass exhibited a complex interaction in determining the functional consequences of land-use transformation. Whole-community energy fluxes and ecosystem functioning Aiming to visualize the complex interplay between community biomass dynamics and energy flux, we constructed energy networks for the four transformation systems (Fig. 3) based on total energy fluxes as a promising way to quantify multitrophic ecosystem functioning (Fig. 1). In addition to the general decreases in biomass (node sizes in Fig. 3) and energy-processing rates (arrow widths in Fig. 3), we also found a systematic shift from predator to omnivore dominance when comparing forest and oil palm systems. Specifically, we found that predator biomass in oil palm yielded only 25% of their biomass in forest (0.424 and 1.664 kg ha−1, respectively), while the predator-driven energy flux was reduced to 46% of the energy flux driven by predators in forest (30.697 and 66.816 kg ha−1 per year, respectively). In contrast, omnivore biomass in oil palm was 22% higher than in the forest (0.767 compared with 0.629 kg ha−1), while omnivore-driven energy flux in the oil palm was 47% lower than in forest communities (32.531 compared with 61.900 kg ha−1 per year; Supplementary Table 2), suggesting a considerable mismatch of biomass and energy flux, partly dependent on the trophic group in question. In our analyses, this disparity finds its explanation in varying body-mass distributions (Supplementary Fig. 1) and assimilation efficiencies that strongly modify how biomass translates into total resource assimilation rates (Fig. 1). These results suggest that biomass, alone, may be an unsuitable proxy for general ecosystem functioning in animal communities. Multitrophic biodiversity-ecosystem function relationships Until now, most studies investigating biodiversity–ecosystem function relationships have focused on single trophic levels31 32. We present a new approach to easily quantify multitrophic ecosystem functioning, requiring only information on body mass, phylogeny, temperature and assimilation efficiencies to overcome previous limitations in biodiversity–ecosystem functioning research. Utilizing this approach, we also investigated the relationship between species richness and ecosystem functioning, identifying a clear linear positive effect of diversity on total energy flux (Fig. 4a and Supplementary Table 3). The relationship between diversity and energy flux was dependent on land-use transformation system, whereby oil palm and jungle rubber showed the strongest decrease in energy flux per unit loss in species richness (Fig. 4a and Supplementary Table 3). Our results suggest that each loss of species in oil palm and jungle rubber therefore would be followed by proportionately higher losses in energy flux, compared with equal species losses in forest and rubber. We found the same pattern as in the overall trend for the predator group, which showed transformation system-dependent relationships between species richness and energy flux (Fig. 4b). However, for omnivores, detritivores and herbivores there was a linear effect of diversity on energy flux driven by these groups; however, this effect was independent of transformation system (Fig. 4b and Supplementary Table 3). This implies that studies focusing on single trophic levels, or even specific species, may fail to detect the alteration of ecosystem processes resulting from land-use transformation. These results call for a wider application of multitrophic approaches that not only measure one ecosystem property, such as total productivity or decomposition, but that also aim to assess whole-community ecosystem processes such as total energy flux. Discussion Our study reflects previous findings that the transformation of forest systems to oil palm has severe impacts not only on single animal populations but also on communities as a whole. In particular, species richness and animal biomass are most significantly affected. Furthermore, jungle rubber and rubber appear to represent intermediate steps in land-use intensification. Their higher levels of biodiversity and ecosystem functioning indicate that they potentially provide higher ecological value than oil palm. As such, these rubber land-use systems could present economically viable, lower intensity land-use alternatives. By taking a multitrophic ecosystem functioning approach we demonstrate that, at the community level, species loss leads to a direct linear decrease in ecosystem functioning. This means that any species loss will be followed by a proportionate loss in function, and this relationship becomes proportionately stronger in more intensive transformation systems such as oil palm plantations. Thus, every one of the few species in high-intensity land-use systems is functionally more important than species in low-intensity systems where functional redundancy is likely to be higher33. Without explicit consideration of multiple trophic levels, such emergent properties are likely to be overlooked. Our study demonstrates the crucial implications of tropical land-use intensification for biodiversity and ecosystem functioning across multiple trophic levels, suggesting that these globally important impacts will likely resonate beyond previously explored trophic boundaries. Methods Study site and sampling design Sampling took place in the Jambi province of Sumatra, Indonesia, a region known as a hotspot for biodiversity, but that has also already undergone extensive deforestation6 34. In the second half of the last century, Sumatra’s forests have experienced vast transformation to rubber and oil palm monocultures35 36. This large-scale land-use conversion has left Sumatra with a very limited area of natural forest mainly restricted to national parks and even here, where logging has been reduced, it has not come to a complete halt37. This severe and extensive land-use transformation, that has progressed already further than in most other tropical landscapes, makes Sumatra a unique and ideal example system for studying the impacts of land-use conversion on biodiversity and ecosystem functioning. We sampled secondary rainforest, jungle rubber, rubber and oil palm systems, replicated eight times across two landscapes (n=32; Supplementary Fig. 2). Sites were selected by first looking for landscapes in the Jambi province that still contained secondary rainforest. Second, we identified all lowland areas with little or no slope and then randomly selected two landscapes with 16 sites each. Among all of the 32 sampling sites, we maintained a minimum distance of 120 m to insure independence of the epigaeic invertebrate communities sampled. The secondary-forest regions lie within two protected areas, Bukit Duabelas National Park and Harapan Rainforest, and represent the least influenced land-use system. Jungle rubber—forest stands with a high percentage of rubber trees that are still regularly harvested—represents a low-impact agroforestry system38. Rubber and oil palm plantations serve as locally common36 high-impact monocultures. The 32 sites were carefully selected so that they were all of a similar age and from equal elevations close to the sea level. All agricultural systems (jungle rubber, rubber and oil palm) were treated and harvested by their owners with intensities typical for the respective transformation system. Animal sampling and calculation of response variables Animal sampling took place between early October and early November 2012. All organisms were collected based on Permit No. 51/KKH-5/TRP/2014 issued by the Indonesian Institute of Sciences and the Ministry of Forestry. In all 32 of the 50 × 50-m sites, we sampled once in each of three 5 × 5-m subplots by sieving the leaf litter from 1 m2 through a coarse sieve of 2-cm-width mesh. In all, 7,472 macroinvertebrates were hand-collected from the sieving samples and stored in 65% ethanol. Specimens were identified to morphospecies and assigned to one of four feeding guilds: omnivores, detritivores, predators and herbivores, based on morphology and literature. As biodiversity studies always suffer from undersampling and correlation of sample size with species richness, we compared observed species richness to both extrapolated and rarefied species richness, calculated in the ‘vegan’ package in R39, to assess the accuracy of our species-sampling effort. To extrapolate sampled species richness, we used the nonparametric second-order jacknife estimator40 to calculate extrapolated species richness from the three 1-m2 subsamples at each of the 32 sites, revealing an estimated mean sampling coverage of 56% (s.d. of±2.393%) making the second-order jacknife estimator the most accurate extrapolation method40. In addition, we calculated sample-based rarefaction, whereby rarefaction curves were calculated for each of the 32 sampled sites and then cut off at the sample size of the smallest sample (40 individuals). Because of the very high attrition of data during the rarefaction procedure (a total of 6,192 out of 7,472 individuals, or 83%, were removed), the rarefied species richness yielded very little resemblance to observed species richness when comparing across transformation systems, resulting in almost no pattern of rarefied richness among transformation systems (Supplementary Fig. 3). The jacknife2-extrapolated species richness, however, was extremely closely correlated with observed species richness (Pearson’s ρ=0.993) patterns among transformation systems (Supplementary Fig. 3), suggesting that our observed species richness did in fact accurately capture realistic patterns in total species diversity across the land-use transformation systems. For each of the 7,472 animals collected, we measured individual body length to an accuracy of 0.1 mm using stage micrometres. We then converted all measured individual body lengths to fresh body mass using length-mass regressions and, where necessary, dry mass-fresh mass relationships from the literature (Supplementary Table 4), yielding an estimated fresh mass in milligrams for every collected individual. Where family-specific relationships were not available or animal body lengths in our collection fell outside of the size ranges of published regressions, we then used regressions from higher-order taxonomic groupings. For heavily damaged individuals that could not be measured for body length, we assigned these individuals a fresh body mass from the median body mass of all animals from the same species or order where only one individual of that species was collected. We then calculated community biomass (mg fresh mass m−2) for each of the 32 communities by summing together all individual body masses calculated from length-mass regressions as derived from the individually measured body lengths. We calculated individual metabolic rates for all 7,472 animals using body masses, temperature and phylogeny24 (Supplementary Table 5). Temperature was measured over a period of at least 2.5 months at 30-cm depth below the soil surface in each site and averaged for each transformation system in each of the two landscapes. From this, community metabolism was calculated by summing together all individual metabolic rates within each of the 32 sites, providing the total metabolic demand for each of the 32 communities. Using diet-specific assimilation efficiencies25, energy loss to predation and community metabolism, we analytically calculated energy fluxes for each of these communities18 using the formula where F is the total energy flux into the network node of a feeding guild, e a is the diet-specific assimilation efficiency, X is the metabolic demand of the feeding guild and L is the loss to predation that the feeding guild is subjected to (Fig. 1 and Supplementary Methods). In order to calculate the fluxes between the functional feeding guilds, we constructed a general network of feeding relationships (link structure in Fig. 3) that represents a null model for an energy network structure where no active preferences are assumed. We assumed that, of our four functional feeding guilds, energy fluxes to predators were split up equally into the three animal guilds below them. Energy fluxes to detritivores and herbivores were assumed to come from only detritus and plant material, respectively. Omnivores were assumed to receive energy in equal 25% proportions from the other three functional feeding groups (predators, detritivores and herbivores, making 75%) and the remaining 25% from both plant and detritus material combined (Supplementary Methods). To assess how these assumptions of feeding preferences might affect the calculations of total energy fluxes, we reconstructed the energy networks so that omnivores were assumed to only consume plant and detritus material (50% derived from each) but with no energy derived from animal material. We then recalculated total energy fluxes and found an overall decrease of up to 54%, which appeared to be highly consistent among the different land-use transformation systems. This consistency between models was especially evident after calculating the loss of energy flux in the three agriculturally used systems compared with the forest system, demonstrating a maximum of only 3% disparity between the two models (Supplementary Fig. 4). This sensitivity analysis indicated that our presented method is highly robust in calculating differences in energy fluxes among different systems. Accordingly, the null model was accepted as the simplest model with the least diet preferences assumed. However, we still suggest that studies adopting this method of energy flux calculation should assign feeding preferences with caution, or employ other techniques such as stable isotope analysis to estimate feeding preferences. Statistical analyses Using mixed effects models (GLMM’s), we tested the effects of ‘transformation system’ and its interaction with functional feeding guild on community responses, with ‘landscape’ as a random effect. ‘Density’, ‘biomass’ and ‘community metabolism’ were log10-transformed to meet assumptions of normality and ‘species richness’ (overdispersed poisson-distributed data) was modelled on a negative binomial distribution. We additionally explored biodiversity–ecosystem functioning relationships by first testing for linearity of relationships using untransformed data. Once linearity was established, we then tested for the effects of log10-transformed ‘species richness’ and its interaction with ‘transformation system’ on ‘energy flux’ for overall data and repeated again for data from separate feeding guilds. In addition, because we suspected that our analyses could be affected by spatial autocorrelation, we calculated Moran’s I values for each model’s residuals and tested for spatial autocorrelation using the Moran’s I standard deviate41 in the ‘spdep’ package in R 3.0.2 (ref. 39). Results from these tests provided no support for the spatial autocorrelation of variation in any of the response variables tested (all Moran’s I test results yielded P>0.4). For all GLMM’s, we applied a backward stepwise selection procedure to obtain the model of best fit, based on the Akaike Information Criterion (AIC). In this procedure, we constructed full models that contained all possible predictors and their interactions (‘transformation system’ and ‘feeding guild’ for general community response models; ‘species richness’ and ‘transformation system’ for biodiversity–ecosystem functioning models) and compared these full models and the model of the backward selection procedure to a null, intercept-only model. The model that yielded the lowest AIC score, with a minimum ΔAIC of 2 units, was selected as the model of best fit. All analyses were conducted with the ‘nlme’ and ‘lme4’ packages in R 3.0.2 (ref. 39). Author contributions A.D.B., M.J. and U.B. designed the study; A.D.B., M.J. and S.M. carried out the field and laboratory work; A.D.B. and M.J. prepared and analysed the data; all authors interpreted the results and wrote the paper. Additional information How to cite this article: Barnes, A. D. et al. Consequences of tropical land use for multitrophic biodiversity and ecosystem functioning. Nat. Commun. 5:5351 doi: 10.1038/ncomms6351 (2014). Supplementary Material Supplementary Information Supplementary Figures 1-4, Supplementary Tables 1-5, Supplementary Methods and Supplementary References

                Author and article information

                Revista Brasileira de Entomologia
                Rev. Bras. entomol.
                Sociedade Brasileira De Entomologia (São Paulo, SP, Brazil )
                : 65
                : 2
                Diamantina Minas Gerais orgnameUniversidade Federal dos Vales do Jequitinhonha e Mucuri orgdiv1Departamento de Agronomia Brazil
                Belo Horizonte Minas Gerais orgnameUniversidade Federal de Minas Gerais orgdiv1Departamento de Biologia Geral Brazil
                Montes Claros Minas Gerais orgnameUniversidade Federal de Minas Gerais orgdiv1Instituto de Ciências Agrárias Brazil
                Viçosa Minas Gerais orgnameUniversidade Federal de Viçosa orgdiv1Departamento de Entomologia/BIOAGRO Brazil
                S0085-56262021000200210 S0085-5626(21)06500200210

                This work is licensed under a Creative Commons Attribution 4.0 International License.

                Page count
                Figures: 0, Tables: 0, Equations: 0, References: 67, Pages: 0
                Product Information: website

                Zelus armillatus, Aphids, Sycophila sp., Pequi, Eurytoma sp.


                Comment on this article