Spatial and temporal distribution of species diversity in semideciduous seasonal forests with occurrence of fire

The spatial and temporal distributions of vegetation can be influenced by edaphic and environmental factors as well as disturbances. The present study was undertaken to characterize the distribution and spatial dependence of species diversity in a fragment of semideciduous seasonal forest after a disturbance event by fire and to identify changes during natural regeneration. The fire, occurred in 1999, and surveys were undertaken in 2000, 2005 and 2014 in 26 geo-referenced plots (400 m2) distributed along two transects. The Fire Damage Index (FDI) in each plot was based on a scale from 0 to 5, and all of the arboreal individuals with diameters at breast height (DBH) ≥ 5 cm were identified in each plot. Species diversity was calculated using the Shannon index (H’), and species richness (S) was estimated using species accumulation curves; variations between the surveys compared based on the means and standard errors. The S, H’ and FDI data were submitted to non-spatial exploratory and geostatistical analyses. After adjusting the semivariograms, the levels of spatial dependence were classified, and interpolation of the variables values were performed using ordinary kriging to characterize their spatial distributions in the form of maps. Spatial analysis was used to identify and characterize differences between the post-fire surveys in terms of the configurations of the arboreal community. The fire event influenced the spatial and temporal structures of the variables S and H’. These variables showed spatial dependence and aggregated distributions, with reduction in the distance under spatial influences and a uniformity of individuals distribution in the forest fragment at the different surveys. This research characterized the distribution and spatial dependence of the variables S and H’ in a forest fragment after a fire event and the alterations in the arboreal community structure during natural regeneration. Rev. Biol. Trop. 65 (2): 525-534. Epub 2017 June 01.

Studies of the uncertain, ambiguous, and nonlinear natures of forest community disturbances will require the use of statistical tools capable of measuring those interactions to be able to increase the efficiency of remediation strategies.Geostatistics can be used to represent the spatial distributions of environmental factors and examine the ecological aspects of recovery and these tools have been used by numerous researchers to examine the spatial dependence of forest components (Rufino, Thiersch, Ferreira, Kanegae-Junior, & Fais, 2006;Amaral, Ferreira, Watzlawick, & Genú, 2010;Neves, Lemos, Gonzáles, Vieira, & Siqueira, 2010;Amaral, Ferreira, Lisboa, Longhi, & Watzlawick, 2013a).These studies can help generate a better understanding of the establishment and reproduction of individuals (Condit et al., 2000), the maintenance of species diversity, succession, evolution, and adaptations (Legendre & Fortin, 1989), and optimize techniques for sustainable forest management, and sampling and monitoring of plant species (Pereira, Péllico-Netto, & Carvalho, 2006).Spatial distribution studies require long-term monitoring of the ecological aspects of forest communities and must accompany recuperation processes and their rates long after a disturbance event (Taylor, Hamilton, Lewis, & Nantle, 2008).
Thus, this study was developed to characterize the distribution and spatial dependence of species richness and diversity in a forest fragment, after a disturbance event by fire, and to identify alterations in the arboreal community structure during natural regeneration.

Study area:
The present study was conducted in a remnant of seasonal semideciduous montane forest (Instituto Brasileiro de Geografia e Estatística, 2012), in the municipality of Ibituruna, Minas Gerais State, Brazil (21°09'40'' S -44°50'10'' W).The fragment covers an area of 57.06 ha in a mountainous region, with hillside inclinations between 25° and 40°.The regional climate is humid in the Austral summer and dry in winter (Cwb by the Köppen classification), with an annual mean temperature of 20.4 °C and an annual mean rainfall of 1 460 mm (Dantas, Carvalho, & Ferreira, 2007).The soils are predominantly litholic neosols, haplic cambisols, and red argisols (Empresa Brasileira de Pesquisa Agropecuária, 2006).The forest remnant studied was surrounded by pastureland, coffee plantations, fragments of native vegetation, and a lake formed by a hydroelectric dam.A forest fire occurred in the area in August/1999, affecting an irregularly shaped area in a West to East path through the fragment (Silva, Filho, Venturin, Carvalho, & Gomes, 2005).

Delimitation of the plots and evaluations:
In January 2000, twenty-six 400 m² plots (20 × 20 m) were established along two transects parallel to the Western and Eastern borders of the fragment to represent the environmental variability extant in the areas impacted by fire.The plots were arranged approximately 30 m one from the other along each transect.All of the plots were geo-referenced.
Floristic surveys were undertaken in 2000, 2005 and 2014 by marking, recording, and identifying all individuals with diameters at breast height (DBH) ≥ 5 cm in each plot, and then calculating their Shannon diversity indices (H') (Magurran, 1988) and species richness (S).
The species accumulation curves were plotted and the variations between the different surveys compared, based on their means and standard errors.Comparisons between the curves were based on the 95 % confidence interval (Zar, 2010) using EstimateS version 9.1 software (Cowell, 2013).
The Fire Damage Index (FDI) was calculated based on the scale proposed by Nascimento et al. (2000), with evaluations varying from zero (absence of fire in the plot) to five (maximum fire intensity) for each plot.

Geostatistical analyses:
The S, H' and FDI values were used in the geostatistical analyses in association with the central coordinates of each plot.The data were initially submitted to non-spatial exploratory analyses to calculate the mean, maximum and minimum values, the coefficients of variation, and standard deviations.Based on these exploratory analyses, geostatistical analyses were undertaken to characterize the dependence and spatial variations of the variables.To analyze the structure of their spatial dependence, variographic analyses were made of the S, H' and FDI data.The experimental semivariograms were adjusted to the exponential and spherical models (Journel & Huijbregts, 1978) using the ordinary least square method.The adjusted models were used to estimate the nugget effect (C 0 ), sill (C 0 +C), and range (A) parameters, which were used, in turn, to calculate the degree of spatial dependence (SD) using the equation: [C/(C 0 +C)]×100.Spatial dependence was classified following Biondi, Myers and Avery (1994) with SD considered as weak if ≤ 35.0 %, moderate if > 35.0 % and ≤ 75.0 %, and strong if > 75.0 %.
After adjusting the semivariograms, the data were interpolated by ordinary kriging to characterize, in the form of maps, the spatial dependence structures of the variables in each survey.The R computational system was used for the statistical and geostatistical analyses and for the plotting of maps (R Development Core Team, 2012) using the geoR packages (Diggle & Ribeiro-Junior, 2007) and ArcGis, version 10.1, with the Geostatistical Analyst extension (Environmental Systems Research Institute, 2010).

Species richness and diversity:
No differences were found between the different surveys, based on the accumulation curves of S and H' variables.Some increases in those variables, proportional to species abundances, were observed, although the curves of these variables for all surveys (2000, 2005 and 2014), showed the same patterns and did not differ among themselves (Fig. 1).
After the fire, the number of individuals and species detected was similar: 1 046 individuals were catalogued in the year 2000, 979 in 2005, and 1093 in 2014; these represented 159, 154, and 163 species, respectively.The most abundant species in all three forest surveys were Cupania vernalis, Gymnanthes klotzschiana, and Prunus myrtifolia; but only C. vernalis and G. klotzschiana were among the ten species with the highest densities in both transects.The H' values for the three surveys were 4.41, 4.39, and 4.40 nats.ind.-1 respectively, indicating that there were no statistical differences between them, and making it impossible to characterize the influence of fire on plant diversity, based on the species accumulation curves for the evaluated years after the fire event.
Spatial continuity analyses: Based on the descriptive analyses of S and H' data, there were alterations of those variables after the fire event and heterogeneity of the data, with higher values in the year 2000 variance (s 2 ), standard error (s), and coefficient of variance (CV).The median values of the variables did not oscillate during the surveys, although there were variations in data amplitudes, with the largest values observed in 2000 and the smallest in 2014 (Table 1), characterizing a de-structuring of the variables as a possible consequence of fire, with posterior stability.The reductions of the variance (s 2 ), and coefficient of variance (CV) values of the S and H' values between the 2000 and 2014 surveys (Table 1), confirm this observation.
Spatial analyses demonstrated alterations in the structural variables following the fire.There was spatial dependence of S and H' variables and the spherical and exponential models were adjusted accordingly (Fig. 2; Table 2).The S variable spatial dependence resulted moderate in all three surveys; while H'  2).Both variables demonstrated reductions in their value ranger, nugget effect, and sill during the series of three surveys (Table 2), that is, a reduction in the distance under spatial influence, and the uniformity of tree distributions in the forest fragment.
The variables FDI, S, and H' showed spatial dependence and aggregated patterns of distribution (Fig. 2 and Fig. 3).The fire spread from East to West and some fragment areas were more severely damaged (FDI) (Fig. 3); these focal points likewise demonstrated smaller S and H' values and a gradient was established from them (Fig. 3).In these cases, the spatial distributions of the arboreal species were dependent on their distances from the focal points of the fire.
There was an increase in areas showing greater species richness and a reduction in areas with high species diversity following the first survey (Fig. 3).

DISCUSSION
It was not possible to characterize the influences of the fire on the structure of the arboreal community studied here, based on numerical analyses and the accumulation curves of S and H' variables.The S values in all three surveys were intermediate in relation to other seasonal  semideciduous forest with known fire occurrences; Peixoto et al. (2012) reported from 148 to 168 species, while Abreu, Pinto and Mews (2014) reported from 102 to 100 species under similar conditions, and the latter group was unable to identify alterations in the forest community after a fire event (as was seen in the present study).
The H' values reported here were intermediate in terms of published results from other forest fragments for the same region (Alto Rio Grande) without reported fires (Botrel, Filho, Rodrigues, & Curi, 2002;Nunes, Mendonça, Filho, Botezelli, & Machado, 2003;Dalanesi, Filho, & Fontes, 2004;Guilherme, Filho, Appolinário, & Bearzoti, 2004;Carvalho, Filho, Fontes, & Curi, 2007).As such, the numerical analyses of the H' variable was not useful to characterize the consequences of the event of fire or alterations resulting from regeneration processes, and could only be characterized based on the spatial analyses of the data.
The differential effects of fire in the sampled plots contributed to the observed variability of the data, although the reductions of the values of the coefficient of variation and the measures of dispersion (s 2 and s), indicated the reestablishment of homogeneous distributions of S and H' after the disturbance.
Aggregated patterns of spatial distribution are characteristic results of fires in small forest fragments (Pereira, Barros, Acerbi-Junior, Pereira, & Reis, 2013).After a fire event, areas (or foci) are created with smaller S and H' values due to high fire damage intensity in those sites, resulting in spatial distribution gradients of the variables analyzed, which are dependent on their distance from those foci.As disturbances of equivalent proportions are not repeated, this interference gradually becomes reduced, characterizing natural regeneration.This result could also be seen in the attenuation of spatial dependence in sequential surveys.In terms of the variable H', the spatial dependence was classified as strong in 2000, but alternated to moderate in 2014; this demonstrated a lower degree of dependence in relation to the initial foci.
The occurrence of clearings in forest fragments can initially result in independence of the different regions of the forest community (Lima, 2005;Lima & Moura, 2006), with consequent spatial dependence on the pre-existing population.These canopy opening sites are initially colonized by species that had been excluded by interspecific competition, but returned and became established in those recently disturbed sites (Fox, 2013).This was observed in the present study, with a reduction in the areas that demonstrated small S values in the 2000 survey, with distancing occurring between new individuals and their mother plants during the process of natural regeneration, with consequent reductions in spatial dependence.After forest reconstruction, the renewed dominance of species adapted to lowlight environments resulted in reduced diversity (Bello et al., 2013).
Natural phenomena, in the case of natural regeneration, introduce spatial conditioning (Amaral et al., 2013a).As such, variations in plant distributions are not random and will demonstrate some degree of spatial dependence, and this dependence cannot be identified through classical analyses (Miller et al., 2007;Batista et al., 2016) as observed in the present study.In that sense, as described by Amaral et al. (2013a), analyses that assume independence among samples, have been substituted by geostatistical analyses that provide more reliable and precise data for decisionmaking by forest administrators.
As such, we were able to identify alterations in the configuration of the arboreal community after a fire disturbance by its influence on the spatial and temporal structures of the variables S and H'.These variables demonstrated spatial dependence and aggregated distributions, with reduction in the distance under spatial influences during the evaluation period, characterizing the stability of the arboreal community.During the process of natural regeneration, the distributions of these variables remained aggregated, although the more advanced successional stage demonstrated lower spatial dependence.

ACKNOWLEDGMENTS
The authors would like to thank the team at the Laboratório de Silvicultura e Ecologia of the Departamento de Engenharia Florestal of the Universidade Federal de Lavras (UFLA) for their assistance with the fieldwork.K. G. Silva thanks the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for the doctoral scholarship grant.

Fig. 1 .
Fig. 1.Species accumulation curves as a function of abundance in a semideciduous seasonal forest in Ibituruna, MG, Brazil, in 2000, 2005 and 2014. A. Mean species richness and B. Shannon Diversity Index.

Fig. 2 .
Fig. 2. Spherical and exponential semivariograms for the variables of species richness (S), Shannon diversity index (H') and the Fire Damage Index (FDI) in a fragment of semideciduous seasonal forest in Ibituruna, MG, Brazil, during plant surveys undertaken in 2000, 2005 and 2014.

Fig. 3 .
Fig. 3. Kriging maps for the species richness (S), Shannon diversity index (H'), and Fire Damage Index (FDI) variables in a semideciduous seasonal forest fragment in Ibituruna, MG, Brazil, in three floristic surveys undertaken in 2000, 2005 and 2014.

TABLE 1
Statistical summary of the species richness (S) and Shannon diversity index (H') variables in three floristic surveys undertaken in a fragment of semideciduous seasonal forest in Ibituruna, MG, Brazil was stronger for the 2000 and 2005 surveys, and only moderate in 2014 (Table