Research Article Volume 2 Issue 5
Consejo Nacional de Investigaciones Científicas y Tecnológicas, Argentina
Correspondence: Gönc RL, Consejo Nacional de Investigaciones Científicas y Tecnológicas, Ciudad Autónoma de Buenos Aires, CABA, Argentina
Received: August 21, 2018 | Published: September 17, 2018
Citation: Gönc RL, Casaux R. Effects of disturbances generated by different management strategies at Nothofagus antarctica forest on beetles assemblage of Chubut, Argentina. Biodiversity Int J. 2018;2(5):415-423. DOI: 10.15406/bij.2018.02.00094
Nothofagus antarctica is a species of the Andino-patagónico forest with a wide distribution from the north of Neuquén to Tierra del Fuego. N. antartica forests are intensely affected nowadays by human activities such as grazing, firewood, conversion into pastures and the replacement by conifer plantations. These activities produce changes in N. antarcticaforest, mainly by the alteration of its original structure affecting the performance of the ecosystem as a whole, resulting in the loss of biodiversity. In this sense, little is known about the effects of these disturbances on Coleoptera communities that contribute to faunal biodiversity and play key roles in the ecosystem´s dynamics. The response of beetle assemblages to disturbances and the dependence of that response on habitat type were analyzed in this study. Beetles were collected by fluid interception and sticky traps, and by foliage sections shaking; and plant species richness and vegetation cover were measured.
A total of 354 beetles representing 17 families were captured, with 25 species identified and 8 not identified. Multivariate analysis revealed that disturbances affected the composition of the beetle assemblage similarity. According to this analysis, the composition of beetle was different among sampling sites FWL, B25 and PR; while it was similar among sites B15, FWL and PC sampling sites. We found a significant relationship between abundance and structural vegetation characteristics, and this was a positive relationship with shrub stratum cover in PR and B25 sampling sites. Finally, we found an increase in the abundance of beetles with vegetation richness in B25 and PR samplings sites, and these increases had a positive relationship to shrub stratum richness in both sites. These results suggest that in Patagonia, the disturbances on Nothofagus antarctica forests can have a substantial effect on beetle communities, and the consequences of these perturbations depend on habitat type. The pre-disturbed conditions of each habitat type and the distribution of species among the different habitats and how these species use the habitats determined their ability to deal with pos-disturbance conditions. Disturbances produce simplification and homogenization of habitat types and this could reflected by beetle diversity, so a good knowledge of the beetle usage of habitat and how the habitat utilization may change along a regional gradient is necessary to understand the effect of disturbance on native beetle assemblages.
Keywords: coleoptera, Nothofagus antarctica, disturbance, Patagonia Argentina
Biodiversity continues to decline in the face of increasing anthropogenic pressures, such as habitat loss and fragmentation, biological invasions, exploitation, pollution and introduction of exotic species.1 Understanding the factors behind biodiversity loss is one of the most urgent tasks faced by scientists at present. In Patagonia, the consequences resulting from human activity on the forest, mainly Nothofagus,2 and in particular of N. antarctica forest, are the loss of diversity, degradation, fragmentation of habitat, overexploitation of species, introduction of exotic species,3,4 and the like. These losses could be due to changes in the forest structure, microclimatic conditions or nutrient cycles.5,6
The habitat type (distinguished by complexity level) is one of the most important factors as it usually modulates the effect of other agents on natural communities.7,8 In some ecosystems habitat may become less complex after disturbance because of the reduction of vegetation cover,9 so, in consequence, the disturbance may have different effects on communities depending on the type of habitat, and the response of communities will vary depending on the capacity of their species to deal with the post – disturbance conditions.9 Species adapted to more complex habitats, for example habitats with more vegetation stratums, which can mean more environmental niches and more foraging opportunities for the species, will be more severely affected than species that prefer less complex habitats (e.g. habitats with less vegetation stratums).10
Beetles are tightly linked to ecosystem functioning: they can attain large biomass and represent a great proportion of the total species diversity of an ecosystem.11 These taxa have been shown to be sensitive to habitat disturbance, like burning,12 and have also shown to be useful bioindicators for environment monitoring and assessment, with their high diversity and sensitive to environmental conditions, providing a view of ecological change.6,13,14
Few studies have documented changes in post-disturbance insects and arthropods in the Andean forest, but these studies have addressed principally the effects of fire on beetle assemblages.10,15,16 Improving our knowledge of beetle assemblages in the context of forest disturbances and management is necessary for effective arthropod conservation.17 In this study we analyzed the effect of disturbances on beetle (Coleoptera) assemblages in N. antarctica forests. The specific aims of the study were: (1) to compare beetle abundance and richness across different forest disturbance sites, (2) to describe the relationship between beetle diversity and vegetation structure in these habitats and (3) to determine whether beetle abundance, richness and vegetation variables could explain the response after disturbance.
Sampling sites
For the selection of Nothofagus antarctica forest sectors affected by different disturbances we used satellite images shared by the staff of the INTA EEA San Carlos de Bariloche (Río Negro Province) and EEA Esquel (Chubut Province) (Hansen et al., unpublished information) and field information obtained in this study. The selected study sites belong to the same basin (Futaleufú) and the ranges of altitude (between 600 and 840 meters above sea level) and distance among sites (17km between the farthest points) were as narrow as possible. In addition, the original structure of N. antarctica forest in all sampling sites were similar (patches with N. antarctica individuals between 4 and 8meters in height).
The disturbances on N. antarctica forest considered in this study were: 1) firewood extraction by patches (FWP), 2) firewood extraction by landscaped (FWL), 3) conversion of forest into pasture (PC), 4) fires occurring 15 years ago (B15), and 5) fires occurring 25 years ago (B25). As a control, we also considered pristine forests (PR). Two replicates were considered for each disturbance type (10 sampling sites) and for the pristine forest (2 sampling sites), and the minimum distance between replicates was 1,500meters.
Sampling design
Beetle sampling
The samples were obtained in two consecutive seasons. Mobile adult individuals of the order Coleoptera were collected every three months (every sampling was carried out during the central period of each season) from July 2008 to June 2010. The sampling was carried out in dry and calm days. Fluid interception traps, with an area of 1m2, were used to collect beetles at the leaf litter level. Beetles were collected in five trays of 400cm3 containing a solution of water and antifreeze to prevent the escape of the specimens and to slow their decomposition according to the methodology described in News (1998) and Southwood (1989). At each site we selected one specimen of the three dominant species belonging to the shrub or arboreal stratum which were used to sample beetles associated to the foliage (the same species sampled throughout the study). The foliage was sprayed with insecticide and shaken for thirty seconds, and the beetles that fell were collected in a deep trough.18,19 These samples were integrated in a single sample. Finally, beetles associated with tree trunks were collected with three "Sticky" traps. These traps consisted of clear adhesive plastic of 90cm2, and each one was placed with a different cardinal orientation. This approach was repeated at all the selected trees. The samples were integrated into a single. Sticky and fluid interception traps remained in the sampling sites over a period of 24 hours. Then, beetles were extracted from the traps and preserved in ethanol 70º.
At the laboratory the beetles were determined to species or morpho-species, through the observation of the morphological characters normally used for the determination of each group, with the help of guides and taxonomic keys20‒25 and specialized taxonomists. The material was deposited at the laboratory of forest health of the Asentamiento Universitario San Martin de los Andes, Universidad Nacional del Comahue (AUSMA).
Vegetation sampling
The information was also obtained from July 2008 to June 2010. We randomly selected three sampling points at each sampling site where the following information was obtained in a radius of five meters from such points:
Sampling site |
Stratum |
n |
Mean ± S.D. |
B15 |
Herbaceous |
8 |
28.38±19.04 |
Shrub |
17 |
24.71±19.33 |
|
Arboreal |
2 |
49±15.56 |
|
B25 |
Herbaceous |
20 |
12.95±15.29 |
Shrub |
14 |
15.36±18.39 |
|
Arboreal |
12 |
11.08±13.34 |
|
PR |
Herbaceous |
15 |
6.87±10.05 |
Shrub |
21 |
12.95±17.49 |
|
Arboreal |
4 |
7±2.16 |
|
FWL |
Herbaceous |
18 |
19.5±18.2 |
Shrub |
13 |
9.69±14.07 |
|
Arboreal |
6 |
24.33±12.8 |
|
FWP |
Herbaceous |
16 |
27.88±38.66 |
Shrub |
14 |
20.71±28.79 |
|
Arboreal |
4 |
6.5±4.51 |
|
PC |
Herbaceous |
8 |
7.25±7.5 |
Shrub |
3 |
2.67±1.53 |
Annex 1 Vegetation species composition and stratum cover (mean±S.D.) of Nothofagus antarctica forest at each sampling site (FWP, firewood extraction by patches; FWL, firewood extraction by landscaped; PC, conversion of forest into pasture; B15, fires occurred 15 years ago; B25, fires occurred 25 years ago; PR, Pristine forests)
Data analyses
We measured the effect of forest disturbances on the abundance of Coleoptera and vegetation stratum cover using General Linear Model (GLM) ANOVA, with sampling sites (PR, PC, FWL, FWP, B15 and B25) as factors, and sampling season (8) as repeated measures.16 Tukey tests were used for posterior comparisons.30 Abundances were transformed to meet ANOVA assumptions: x’=log(x+1). Given that the estimation of species richness can be influenced by differences in sample size, expected richness was calculated and compared using sample-based rarefaction curves.31 Rarefaction eliminates variation in species richness due to differences in sample size by re-sampling a pool of n individuals repeatedly at random.31 Rarefaction was used to estimate beetle richness in each sampling site.
Analysis of similarity (ANOSIM) with Bray-Curtis similarity index was used to determine if there were significant differences in beetle assemblages among the sampling sites (B15, B25, FWL, FWP, PC and PR).15,32 ANOSIM is a non-parametric permutation procedure applied to rank similarity matrices underlying sample ordinations. The test is based on comparing distances among groups with distances within groups. ANOSIM produces a R-statistic that is an absolute measure of distance among groups. A large positive R (up to 1) means high dissimilarity among groups. Significance is computed by permutation of group membership. To display the relationship between different sampling sites we used non-metric multidimensional scaling (NMDS) with the Bray–Curtis similarity coefficient.15,33 The NMDS algorithm attempts to place data points in a coordinate system such that ranked distances are preserved. This method is well suited for non-normal data and does not assume linear relationships among variables. ANOVA, rarefaction curves, ANOSIM and NMDS were performed with PAST 1.34.34
To describe the relationship between beetles and vegetation cover, correlations were performed (Pearson’s coefficient) between beetles abundance and vegetation structure parameters (richness and cover stratum). This analysis was made with Infostat Profesional v.2014.35
A total of 2743 beetles belonging to 17 families were represented in the samples; 89.13% of them were assigned to 25 species and the remaining was assigned to eight unidentified species (Table 1). The most abundant families were Curculionidae (26.24% of the total of the individuals sampled) and Staphylinidae (17.68%), followed by Scarabaeideae (12.43%) (Table 1). The most abundant species were Entiminae sp. (Shoenherr, 1823), Anthonomini sp. (Thomson, 1859) and Aleocharinae sp. (Fleming, 1821). The highest abundances were registered at FWP and FWL (95 and 70 individuals respectively), while the lowest abundance was registered at PC and PR (48 and 30 individuals respectively) (Table 1). The highest and the lowest richness were registered at FWL (23 species) and at PR (13 species) respectively (Table 1).
Family |
Species |
Sampling site |
Total Abundance |
Mean±S.D. |
Richness± S.D. |
Anthicidae |
sp. 8 |
PC |
1 |
1 |
0 |
Brupestidae |
sp. 1 |
B15 |
1 |
1 |
0 |
Cantharidae |
Chauliognathus sp. |
FWL |
2 |
1 |
0.24 |
FWP |
3 |
1 |
0.24 |
||
PR |
2 |
2 |
0.19 |
||
Carabidae |
Cnemalobus araucanus |
B15 |
6 |
2±1 |
0±0.14 |
FWL |
1 |
1 |
0.14 |
||
FWP |
2 |
2 |
0 |
||
Mimodromius chilensis |
FWP |
1 |
1 |
0.14 |
|
PC |
1 |
1 |
0.14 |
||
Mimodromius nigrotestaceus |
FWP |
3 |
3 |
0.14 |
|
PC |
1 |
1 |
0.14 |
||
Cerambicidae |
Calydon submetallicum |
B15 |
2 |
1 |
0.14 |
Eryphus laetus |
B15 |
1 |
1 |
0.14 |
|
FWL |
7 |
2 |
0 |
||
FWP |
2 |
1 |
0.14 |
||
PC |
2 |
3.5±3.54 |
0.21±0.51 |
||
Chrysomelidae |
Jansonius aeneus |
B25 |
1 |
1 |
0.19 |
PR |
3 |
1 |
0.19 |
||
Kuschelina decorata |
B25 |
8 |
2.67±2.08 |
0.03±0.2 |
|
PR |
1 |
1 |
0.19 |
||
Mylassa crasicollis |
B25 |
10 |
3.33±2.53 |
0.03±0.24 |
|
FWL |
6 |
1.2±0.45 |
0.17±0.04 |
||
Coccinelidae |
Adalia angulifera |
B15 |
7 |
2.33±0.58 |
0.06±0.06 |
B25 |
1 |
1 |
0.19 |
||
FWL |
2 |
1 |
0.19 |
||
FWP |
1 |
1 |
0.19 |
||
PC |
1 |
1 |
0 |
||
Adalia deficiens |
B15 |
1 |
1 |
0.19 |
|
B25 |
1 |
1 |
0.19 |
||
FWL |
3 |
1.5±0.71 |
0.14±0.07 |
||
Coccinella chilena |
B15 |
1 |
1 |
0.19 |
|
FWL |
2 |
1 |
0.19 |
||
PC |
1 |
1 |
0.19 |
||
Curculionidae |
Alastropolus sp. |
B15 |
2 |
1 |
0.21 |
B25 |
2 |
2 |
0.14 |
||
FWL |
2 |
2 |
0.14 |
||
PC |
2 |
2 |
0.14 |
||
Anthonomini sp. |
B15 |
9 |
2.25±1.5 |
0.13±0.11 |
|
B25 |
1 |
1 |
0.21 |
||
FWL |
3 |
1 |
0.21 |
||
FWP |
13 |
2.6±2.3 |
0.1±0.16 |
||
PC |
6 |
3 |
0.07 |
||
PR |
1 |
1 |
0.21 |
||
Cryptorrhynchini sp. |
B15 |
1 |
1 |
0.21 |
|
B25 |
1 |
1 |
0.21 |
||
FWL |
4 |
2±1.41 |
0.14±0.1 |
||
PC |
6 |
6 |
0.14 |
||
PR |
1 |
1 |
0.21 |
||
Entiminae sp. |
B15 |
5 |
2.5±0.71 |
0.11±0.05 |
|
B25 |
4 |
1 |
0.21 |
||
FWL |
10 |
2±1.22 |
0.14±0.09 |
||
FWP |
9 |
1.29±0.49 |
0.19±0.03 |
||
PC |
1 |
1 |
0.21 |
||
PR |
3 |
3 |
0.07 |
||
Elateridae |
sp. 3 |
FWL |
2 |
2 |
0.29 |
Lucanidae |
sp. 2 |
FWL |
1 |
1 |
0 |
Melyridae |
sp. 4 |
B25 |
14 |
3.5±2.38 |
0.71±0.68 |
FWL |
2 |
1 |
0 |
||
FWP |
4 |
2±1.41 |
0.29±0.4 |
||
PC |
1 |
1 |
0 |
||
PR |
3 |
3 |
0.57 |
||
Nitulidae |
sp. 5 |
FWL |
1 |
1 |
0 |
FWP |
2 |
1 |
0 |
||
Salpingidae |
sp. 6 |
B15 |
4 |
3 |
0.57 |
B25 |
1 |
1 |
0 |
||
Scarabaeidae |
Aphodius sp. |
FWP |
4 |
4 |
0.06 |
Brachysternus sp. |
B15 |
1 |
1 |
0.23 |
|
FWL |
5 |
1.25±0.5 |
0.21±0.03 |
||
PC |
2 |
2 |
0.17 |
||
Oryctomorphus sp. |
FWL |
2 |
1 |
0.23 |
|
Phytolasma sp. |
B15 |
1 |
1 |
0.23 |
|
FWP |
1 |
1 |
0.23 |
||
Sericoides sp. |
B15 |
2 |
2 |
0.17 |
|
B25 |
1 |
1 |
0.23 |
||
FWL |
7 |
1.75±0.96 |
0.19±0.05 |
||
FWP |
1 |
1 |
0.23 |
||
PC |
12 |
4±3 |
0.06±0.17 |
||
PR |
3 |
3 |
0.11 |
||
Scolytidae |
sp. 7 |
PR |
1 |
1 |
0 |
Staphylinidae |
Aleocharinae sp. |
B15 |
2 |
1 |
0.19 |
B25 |
7 |
2.33±2.31 |
0.06±0.22 |
||
FWL |
7 |
3.5±0.71 |
0.05±0.07 |
||
FWP |
31 |
7.75±11.53 |
0.45±1.1 |
||
PR |
4 |
2±1.41 |
0.1±0.13 |
||
Leptoglossula sp. |
B15 |
1 |
1 |
0.19 |
|
B25 |
2 |
1 |
0.19 |
||
FWL |
2 |
2 |
0.1 |
||
FWP |
7 |
3.5±2.12 |
0.05±0.2 |
||
PC |
2 |
2 |
0.1 |
||
PR |
3 |
1.5±0.71 |
0.14±0.07 |
||
Thinodromus sp. |
B25 |
11 |
3.67±3.79 |
0.06±0.36 |
|
FWL |
2 |
2 |
0.1 |
||
FWP |
3 |
1 |
0.19 |
||
PC |
2 |
2 |
0.1 |
||
PR |
4 |
2 |
0.1 |
||
Tenebrionidae |
Nyctelia rotundipennis |
FWL |
1 |
1 |
0 |
FWP |
6 |
1.5±0.58 |
0.14±0.17 |
||
PC |
1 |
1 |
0 |
Table 1 Beetle abundance and richness (mean±S.D.) in each sampling sites (FWP, firewood extraction by patches; FWL, firewood extraction by landscaped; PC, conversion of forest into pasture; B15, fires occurred 15 years ago; B25, fires occurred 25 years ago; PR, Pristine forests)
The abundance of beetles was similar among sampling sites (ANOVA, F5;142=1.01; p=0.4167) and sampling seasons (ANOVA, F7;142=1.33; p=0.2398). The interaction term were not significant (F31;142=1.49; p=0.0628). The expected richness of beetles were similar among sampling sites (F5;142=1.35; p=0.2454) and sampling season (F7; 142=0.78; p=0.6056). The interaction among factors were not significant (F31;142=0.60; p=0.9501).
Coleoptera expected richness at FWL was 23% higher than at the remaining sites (Figure 1c), and was followed in importance by FWP and PC (Figure 1d) (Figure1f). The expected richness at B25, B15, and PR were similar (Figure 1b-1e) and lower than at the remaining sites.
Figure 1 Sample-based rarefaction curves rescaled to individuals used to compare expected species richness of beetles (Coleoptera) in sampling sites. Data are means ±SD after randomization.
ANOSIM and NMDS ordination revealed that forest disturbances affected the composition of beetle assemblage (Table 2) (Figure 2). Beetle assemblage composition was similar among B15, FWL and PC sampling sites (Table 2) (Figure 2). In contrast, beetle assemblage was different among FWL, B25 and PR sampling sites (Table 2) (Figure 2).
Figure 2 Non-metric multidimensional scaling ordination of beetle’s assemblage based on sampling sites.
Sampling sites: FWP, PC, X B15, B25, FWL, *PR. For ANOSIM results, see Table 2.
Comparison sampling sites |
R |
p values |
B15 vs. B25 |
0.01 |
0.2126 |
B15 vs. FWL |
-0.02 |
0.8837 |
B15 vs. FWP |
-0.02 |
0.8097 |
B15 vs. PR |
-0.03 |
0.753 |
B15 vs. PC |
0.12 |
0.0414 |
B25 vs. FWL |
0.001 |
0.39 |
B25 vs. FWP |
-0.01 |
0.5914 |
B25 vs. PR |
0.002 |
0.4133 |
B25 vs. PC |
0.09 |
0.088 |
FWL vs. FWP |
-0.04 |
0.8334 |
FWL vs. PR |
-0.02 |
0.6926 |
FWL vs. PC |
0.11 |
0.0346 |
FWP vs. PR |
-0.02 |
0.6412 |
FWP vs. PC |
0.08 |
0.1377 |
PR vs. PC |
0.06 |
0.8958 |
Table 2 Analysis of similarity (ANOSIM) for beetle assemblages between sampling sites (FWP, firewood extraction by patches; FWL, firewood extraction by landscaped; PC, conversion of forest into pasture; B15, fires occurred 15 years ago; B25, fires occurred 25 years ago; PR, Pristine forests)
Statistically significant relationships were observed between beetle abundance and vegetation richness and cover. Beetle abundance increased with vegetation cover at B25 (r=0.35; p=0.0174) and PR (r=0.50; p=0.0011). These positive relationships were also observed when the beetle abundance was compared with the shrub stratum cover at both sites (r=0.54; p=0.0462 at B25; r=0.67; p=0.0008 at PR) (Figure 3). Beetle abundance also increased with vegetation richness at B25 (r=0.35; p=0.0174) and PR (r=0.45; p=0.0039) and with shrub stratum richness at B25 (r=0.54; p=0.0462) and PR (r=0.61; p=0.0030) (Figure 4).
Beetle assemblages were affected by disturbances on N. antarctica forest by different mechanisms. Although not found significant differences in the abundance of individuals among the sampling sites, may be due to the number of species is relatively low compared with previous studies of bettle species.15 However the expected richness and abundance showed the highest value in FWL and FWP among other sampling sites. Sackmann & Farji-Brener15 proposed that the reduction in abundance of some families was compensated for by an increase in another; thus, overall Coleoptera abundance was similar between disturbance and in not disturbance plots, and some species could have been benefited from recent post-disturbance conditions and immigrated to perturb other area; furthermore, some could survive.16
Differences in beetle species composition were found among sampling sites rather than among sampling season. Beetle species composition was similar in B15, FWL and PC, but it was different among FWL, B25 and PR sampling sites. This pattern of beetle succession could occurred at the same time as the vegetation recovery patterns, although vegetation changed among habitats, so that species beetle composition could be influenced by recolonization and recovery rate of communities of plants after disturbances.36
Vegetation cover represents aspects of the productivity of habitats that may have a significant positive impact on beetle abundance.37 The intrinsic heterogeneity of the habitat and spatial segregation of taxa within the habitat (use of different microhabitats), and the consequent simplification post disturbance affect the distribution and frequently of beeltes families.38,15 In this study, shrub cover is the most important predictor of beetle abundance. This association suggests that shrubs may offer foraging and oviposition sites for beetles, thus promoting a local increase in abundance. Some authors39,40 suggest that the positive association of canopy cover with beetle abundance could explain the high sensitivity of beetles to habitat disturbances in vegetation structure caused by humans. Canopy closure has shown to favor distinct, small-scale microhabitats at the soil surface, generating a diverse mosaic of environmental conditions that affect the abundance and distribution of species.41 An increase in shrub cover probably provides more food resources and shelter against microclimatic variation, which may in turn increase beetle abundance in sampling sites. In this study, the effect of stratum cover on beetle abundance varied among sampling sites. Shrub cover affected beetle abundance positively in B25 and PR, and did not affect the abundance of beetles in the remaining sampling sites. Similar results were found by Ruggiero et al.,42 and Sasal et al.,16 in Nothofagus forests, and they suggest that in some environments the shrubs may act as “keystone” structures, which beetles may actively choose, as they provide important resources of food, thermal refuges, and protection against predators.43,44 Bess et al.,36 showed that certain Coleoptera species preferred open, bare ground habitat and that they increased as well after disturbances like fire. On the other hand, we did not find an association among beetle abundance and arboreal and herbaceous stratums. In contrast, Ruggiero et al.42 found a negative association between these variables. So we cannot explain this interesting pattern. For instance, forest beetles are closely associated with variation in understory plant cover in the Northern Hemisphere.37
Pre-disturbance conditions of each habitat type determined different post-disturbance communities’ response. Beetle assemblages pattern was influenced by habitats structure rather than season conditions. Disturbance produced simplifications and homogenization of habitat types, and this was reflected by beetle diversity. Our results support the idea that the consequences of disturbance on beetle assemblages depend on habitat type (e.g. differences in structural complexity of vegetation), so it is possible that beetle assemblages recover as vegetation does.36 These findings suggest that beetles could provide effective indicators of disturbance, but further studies would be necessary to understand the factors underlying in the relationship among beetle abundance vegetation structures. In particular the distribution of species among the different sampling sites and the way that those species use the habitats (e.g. preference for shrub stratums) will determine their ability to deal with post – disturbance conditions. This implies that a good knowledge of the beetle usage of habitat and how the habitat utilization may change along a regional gradient is necessary to understand the effect of disturbance on native beetle assemblages.
We deeply thank M. Archangelsky and S. Tiranti who helped with species determination, and J. Di Rienzo for assistance with the statistical analyses. We are also grateful to EEA INTA Esquel, Héctor Jones, Ronaldo Jones, Julio Sale Rodríguez and Alejandro Garzonio for permission to access to sampling sites, and L. Sosa Pfatschbacher and A. H. Maranghello for helpful comments on the manuscript and their invaluable help and support during the realization of this study.
The author declares that there is no conflict of interest.
©2018 Gönc, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.