Journal Information
Vol. 19. Issue 1.
Pages 69-80 (January - March 2021)
Download PDF
More article options
Vol. 19. Issue 1.
Pages 69-80 (January - March 2021)
Research Letters
Open Access
Plant diversity conservation in highly deforested landscapes of the Brazilian Atlantic Forest
Julia R.S.A. Mangueiraa,
Corresponding author

Corresponding author at: Rodovia Arquiteto Hélder Cândia, 1684 – casa 38. CEP 78.048-150, Cuiabá, MT, Brazil.
, Leandro T.A. Vieirab, Thaís N. Azevedoc, Ana Paula S. Sabinod, Kátia M.P.M.B. Ferraze, Sílvio F.B. Ferraze, Débora C. Rothere,f, Ricardo R. Rodriguesd
a Plant Biology Department, Biology Institute, University of Campinas, Campinas, SP, 13083-862, Brazil
b MACK Herbarium, Mackenzie Presbyterian University, São Paulo, SP, 01302-907, Brazil
c Paulista University - UNIP, São Paulo, SP, 01533-000, Brazil
d Biological Sciences Department, “Luiz de Queiroz” College of Agriculture, University of São Paulo, Piracicaba, SP 13418-900, Brazil
e Forest Sciences Department, “Luiz de Queiroz” College of Agriculture, University of São Paulo, Piracicaba, SP 13418-900, Brazil
f Department of Ecology, Institute of Biosciences, University of São Paulo, São Paulo, SP, 05508-090, Brazil

  • Matrices influence plant diversity of inferior stratum, being higher in pastureland.

  • Time since abandonment positively influences plant diversity of the superior stratum.

  • Beta diversity is high at the studied highly deforested landscapes in Atlantic Forest.

  • In sugarcane matrix, richness and composition were explained by landscape structure.

  • In pastureland matrix, plant diversity was more related to landscape dynamic metrics.

Article information
Full Text
Download PDF
Figures (6)
Show moreShow less
Tables (1)
Table 1. Variables of forest remnants and landscape structure and dynamic parameters selected in the Corumbataí River Basin, Southeastern, Brazil.

Different types of matrices can perform different disturbance regimes over remnant forests, which can ultimately affect plant diversity. To address these issues, we studied 60 plots in forest remnants embedded within sugarcane plantations and pasturelands in the Atlantic Forest of Southeastern Brazil. Our objectives were to evaluate general patterns of plant diversity and to assess the influence of landscape structure and dynamics on plant richness and diversity. We recorded 297 species in two strata. (DBH superior stratum ≥ 4.8 cm > inferior stratum). Overall, diversity of the inferior stratum was higher in pastureland forest remnants than in sugarcane matrices. In the superior stratum, time since abandonment represented a greater influence. Considering the partition of diversity, more than 40% of total diversity was due to the alpha component alone in all cases, but it was significantly lower than expected. Beta components among forest remnants, landscapes and matrices were higher than expected for both strata, whether together or separately. Habitat heterogeneity showed a higher contribution to the total floristic diversity, but among remnants this contribution was even more significant. Forest remnants in pasturelands had a more dynamic land use history and parameters of landscape dynamics were more related to plant diversity. For remnants in the sugarcane matrices, species richness and composition were explained by landscape structure. Our results highlight that there is an important plant community diffused in the forest remnants of human-modified landscapes. Therefore, conserving the majority of remnants and restoring degraded ones must be a key conservation strategy.

Secondary forests
Multiscale diversity
Human-modified landscapes
Corumbataí River basin
Full Text

Most of the remaining forests in the tropical region is restricted to steep areas, hindered by human occupation, agriculture, and represented by small isolated remnants (Metzger et al., 2009, Ribeiro et al., 2009), which are highly disturbed (Ribeiro et al., 2009) and embedded in agricultural landscapes (Melo et al., 2013). Typically, such landscapes are an agricultural mosaic, containing forest fragments with different degradation histories and drivers, that usually include (i) old-growth forest remnants that have never experienced clear-cutting, but are still degraded by anthropogenic disturbances, (ii) secondary forests that are regenerating after a clear-cut, fire or the abandonment of croplands and pasturelands (Gardner et al., 2009; Molin et al., 2017), and (iii) secondary forests under recurrent events of disturbance (Molin et al., 2018; Reid et al., 2018).

Recent studies have reported the importance of such secondary forests to conserve biodiversity, supply ecosystem services, provide stepping stones, corridors, fauna and flora refuges, and other key conservation elements (Arroyo-Rodríguez et al., 2009; Brancalion et al., 2012; Chazdon et al., 2009; Viani et al., 2015). Nonetheless, different recurrent disturbances play important roles in forest succession when defining the plant community (Arroyo-Rodríguez et al., 2015). Disturbances such as wind, fire, temperature increases and exotic species invasion are important components of many ecosystems, and variations in the disturbance regime can affect an ecosystem's structure and functioning (Hobbs and Huenneke, 1992). The effects of disturbances can be expressed at different levels. Within fragments, environmental (soil conditions and light availability) and ecological filtering (demographic stochasticity, seed dispersal limitations, competition, recruitment, etc.) can alter the species composition and the spatial patterns of plant regeneration (Wright, 2010), by selecting species and defining future biodiversity (Arroyo-Rodríguez et al., 2015; Wills et al., 2006; Wright, 2010). For established woody plants, time since disturbance can play an important role on defining species composition, considering time-lag response that can lead to an extinction debt (Tilman et al., 1994). Such effects on both strata can lead to decreasing local diversity or to biotic homogenization at the local level (Püttker et al., 2015). From a regional perspective, disturbances have been observed to increase the compositional and functional differentiation among forest remnants, especially in highly deforested landscapes with long land use histories (Arroyo-Rodríguez et al., 2015). On the other hand, tropical forests under severe landscape disturbances can also be susceptible to floristic convergence (i.e., species homogenization) (Lôbo et al., 2011) driven by the proliferation and dominance of some fast-growing plant species, such as pioneers and vines (Arroyo-Rodríguez et al., 2015; Liebsch et al., 2008; Lôbo et al., 2011; Melo et al., 2007). Both floristic differentiation and homogenization are dependent on the landscape configuration and spatial scale of analysis (Arroyo-Rodríguez et al., 2013).

Landscape characteristics (e.g., amount of vegetation and matrix context) are known to heavily influence population and metapopulation dynamics (Püttker et al., 2011; Tscharntke and Brandl, 2004). Fragmentation, habitat loss or the integration of both mechanisms has led to several changes in landscape structure and configuration (Fahrig, 2003), driving plant and animal populations to reductions or even to extinctions at the landscape scale in tropical forests (Benchimol et al., 2017; Lima and Mariano-Neto, 2014; Morante-Filho et al., 2015). Species diversity in several tropical plant communities has been strongly and positively related to forest cover and connectivity.

The degree of change in species diversity along with a range of different scales (e.g., samples, habitats or communities) is measured by beta-diversity (β), since it provides a direct link between biodiversity at local (α diversity) and regional scales (γ diversity) (Magurran, 1988). In the Atlantic Forest, some studies have addressed patterns of beta diversity in native vegetation (Carneiro et al., 2016; Machado et al., 2016; Sfair et al., 2016; Bergamin et al., 2017; Benchimol et al., 2017), which remains a crucial field in understanding the patterns related to plant community assembly and its maintenance in highly degraded landscapes. The landscape-divergence hypothesis, postulated by Laurance et al. (2007), states that differential human disturbances among forest remnants, such as the different degrees of reduction of local richness (patches) after habitat loss, can lead to differential trajectories of vegetation succession. Those disturbances may increase the compositional and functional differentiation of communities in fragmented landscapes, leading to high beta diversity between patches in a landscape (Benchimol et al., 2017; Sfair et al., 2016). Along with understanding communities' structures along landscapes, it is also important to understand the historical transformations that habitat patches and biodiversity have been submitted to (Molin et al., 2017). This is particularly important in dynamic landscapes, where deforestation, regeneration, harvesting, grazing, urbanization, and other actions have shaped the distribution of these patches throughout the landscape, playing an important role in selecting species (Arroyo-Rodríguez et al., 2015).

Even with the increase in the number of studies about the Brazilian Atlantic Forest, little is still known about its biodiversity. Less than 1% of the remnant area has been sampled, and most of the current knowledge about this domain comes from private lands and areas outside of forest reserves (Lima et al., 2015). In addition, our understanding of the biological dynamics of such human-modified ecosystems (Arroyo-Rodríguez et al., 2013) has been inadequately studied. Therefore, the combination of beta-diversity analyses with different landscape contexts can inform us about the fundamental mechanisms shaping the remaining diversity. Given such importance, this approach must be investigated in the most varied landscapes possible.

Here, we investigated several ecological processes interacting collectively at multiple landscape contexts and scales of analysis. To discuss these issues, we assessed the floristic composition and diversity of forest remnants embedded in human-modified landscapes which are highly fragmented with a dynamic history of occupation. Our specific goals were (1) to evaluate the general patterns of α-plant diversity of the superior and inferior strata among forest remnants in different periods of regeneration, and among sugarcane and pastureland matrices; (2) to analyze the β-plant diversity patterns at different scales (among forest remnants, among landscapes, and among agricultural matrices); and (3) to assess the influence of landscape structure and dynamics (present and past events) on plant species diversity. Considering the studied forest remnants are embedded in different matrices, isolated from each other, and have been experiencing different disturbance histories over time, we would expect a plant community differentiation among forest remnants, among landscapes and among agricultural matrices resulting in higher β-plant diversity. In contrast, we would expect a lower α-plant diversity within forest remnants (among plots) as a result of the limited seed dispersal inter-forests imposed by the landscape structure.

Material and methodsStudy site

The study areas are situated in the countryside of São Paulo state, Brazil, along the Corumbataí River basin. It is located between the latitudes 22°04′46″S and 22°41′28″S, and longitudes 47°26′23″ W and 47°56′15″W, and covers an area of approximately 1700 km2 (Cassiano et al., 2013). Eight municipalities either totally or partially comprise this basin which contains approximately 650,000 inhabitants. Most of the area is flat relief, disrupted only by the basaltic cuestas, which still maintain a considerable portion of native vegetation. The basin is located in the transition of the Atlantic Forest to the Brazilian Savanna (Cerrado) and its vegetation is composed mostly of seasonal deciduous and semi-deciduous forests (Garcia et al., 2006).

The replacement of the forests by agricultural uses goes back to the early 19th century, and nowadays, the major agricultural land uses are sugarcane plantations and pasturelands. However, this region has a very dynamic history with cycles of different crops, mainly sugarcane, coffee, silviculture, and on a smaller scale, row crops and fruticulture (Garcia et al., 2006). Together with increasing urbanization, these dynamic cycles of high-tech agriculture and recurrent impacts on forest remnants (such as regular fires and pesticides) have influenced the distribution of native vegetation along the basin. In 2001, Valente and Vettorazzi (2005) concluded that 90% of forest remnants were smaller than 5 ha and less than 1% were larger than 85 ha. Ferraz et al. (2014) however, observed an increase in forest cover of 8–15% in sugarcane landscapes and 10% in pasturelands. This expansion mostly occurred in proximity to previously existing forest remnants, resulting in the current remnants being highly heterogeneous. Although the authors observed that forest regeneration had overcome forest loss, old-growth forests are being replaced by young secondary forests representing an overall decrease in habitat quality.

Landscape selection, structure, and dynamic metrics

In 2010, six landscapes were selected to be representative of the predominant agricultural matrices of the Corumbataí River basin: three landscapes in the sugarcane matrix and three in the pastureland matrix. The landscapes are 4 km x 4 km in size, and in 2008, presented at least 10% of native vegetation cover and at least 70% of agricultural matrix from the satellite image available when the selection was performed. Photographic interpretation of the images was used and the following land use classes were considered: sugarcane plantations, pastureland, native vegetation cover, orange plantations, eucalyptus plantations, urban areas, and others. One can see further details of the landscape selection in Ferraz et al. (2014). In order to describe the structure of each landscape, we have employed six variables: Mean Forest Age (in years), Local Forest Neighborhood Dominance (in %), Local Forest Cover (%), Elevation (in meters), Forest Proximity (unitless), and Forest Contiguity (unitless). To assess the landscape dynamic, aerial photographs from 1962, 1978 and 1995, and satellite images from 2008 (used as a reference for mapping) were also photo-interpreted and classified using the same land use classes. Land use transitions of each landscape were analyzed to calculate the following metrics: Mean Annual Forest Change Rate and Forest Change Curve Profile. All variables are described in Table 1 and follow the methods proposed by Ferraz et al. (2014).

Table 1.

Variables of forest remnants and landscape structure and dynamic parameters selected in the Corumbataí River Basin, Southeastern, Brazil.

Variables  Unit  Description 
Mean Forest Age (FA)  Years  Calculated as the overlaying of land cover maps, using the difference between the most recent date and the first years of occurrence in the past. It was used here as a surrogate of forest quality in which there is an increase in provide ecosystem services from early regeneration to old-growth forests. 
Elevation  Meters  Obtained from the digital elevation model ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer), which image was acquired in October 2011. 
Local Forest Neighborhood Dominance (FORNEIGH)  The proportion of surrounding cells around each forest focal cell (1 ha) covered by forest. It was used to express the better performance to provide ecosystem services by the interior than the edges of the forests. 
Local Forest Cover (FORCOVER)  Proportion of each focal cell (1 ha) covered by forest. 
Forest Proximity (PROX)  Unit less  Used as a surrogate of local connectivity, the mean proximity index of forest present in a 2 km buffer around forest cells (1 ha). It is considered that more connected patches provide higher levels of some ecosystem services than more isolated ones. 
Forest Contiguity (FORCONT)  Unit less  Calculated as the relative size of forest patches in relation to each landscape. 
Mean Annual Forest Change rate (q)  %/year  It measures the annual forest change rate on focal cells (1 ha) using the annual rate of forest change equation (FAO, 1995). Positive values represent forest increment and negative values represent forest loss over the years. 
Forest Change Curve Profile (FCCP)  Unit less  Calculated to describe when forest changes are concentrated over time: positive values represent forest changes in earlier years, and negative values represent forest changes concentrated in recent years. 

Collinearity describes a situation where two or more predictor variables in a statistical model are linearly related (Dormann et al., 2013). We performed a pairwise Spearman correlation among all variables in a preliminary analysis and no correlation coefficients were greater than 0.70. Based on this, we decided to keep all variables in the analyses.

Vegetation sampling and biotic data

In 2011, we conducted floristic assessments in 12 Seasonal Semi-deciduous forest remnants of the Corumbataí River basin, two located in each landscape; therefore, six remnants embedded in sugarcane matrices and six in pastureland matrices. We considered forest remnants to be all-natural forest habitats in the landscapes, regardless of their stage of succession. In two forest remnants in each of the six selected landscapes, individuals of shrub and tree species were sampled (Fig. 1). In each remnant, five plots of 30 m × 8 m were established to record the individuals of the superior stratum. Individuals of the inferior stratum were sampled in one sub-plot of 30 m × 2 m. In total, we sampled 1.44 ha for the superior stratum and 0.36 ha for the inferior stratum. All plots were established in the edge-interior direction at least 3 m from the edge to avoid large canopy gaps, grasses, and vines. In total, 60 plots were sampled in the Corumbataí River basin, 30 plots in forest remnants embedded in the sugarcane matrix and 30 plots in forest remnants embedded in pasturelands. The distance among sampled remnants varied from 0.2 to 2.5 km in sugarcane landscapes and 0.6–1.2 km in pasturelands. The size and shape of the forest remnants in the landscapes were very restrictive and did not allow us to increase the distances between the sampled plots, but the spatial autocorrelation is an important issue in spatial data at different scales in community ecology (Diniz-Filho et al., 2012). Thus, before statistical tests, for each response variable, we checked the spatial dependence with Moran's I test through 999 permutations at a 5% significant level. We considered the geographical coordinates of the plots using the binary coding scheme in the neighbor's parameters. This analysis is available in "spdep" package (Bivand and Wong, 2018) in R software (R Core Team, 2016).

Fig. 1.

Distribution of the landscapes dominated by pastureland (1–3) and sugarcane (4–6) in the Corumbataí River basin (Southeastern Brazil), with the corresponding land use classes, with a highlight to the two forest remnants and plots sampled in each landscape. Percentage of forest cover in each landscape: 1 = 19,9%; 2 = 14,7%; 3 = 21,4%; 4 = 15,9%; 5 = 12,1%; 6 = 24,0%.


We sampled two forest layers in each plot, all individuals with a diameter at breast height (DBH) ≥ 4.8 cm were considered superior stratum (including arborescent palms) and all individuals with a DBH < 4.8 cm and at least 30 cm high were considered inferior stratum. All sampled individuals were identified at species level by comparison to the ESA Herbarium (Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo) collection and expert consultation. We excluded exotic species and other life forms (vines and herbs, mainly), individuals that we could not classify at least at the family level, and dead tree individuals from the data.

All species were classified into three functional groups (Supplementary material - Table 1), based on literature data and expert consultation. For seed dispersal syndrome, species were classified as Autochoric, Zoochoric, Anemochoric, and Unclassified (sensu Van der Pijl, 1982). For successional groups, we classified species as pioneer, non-pioneer of canopy (which includes species that reaches the canopy and sub-canopy), understory, and unclassified (which includes species lacking appropriate information for classification and also typical species from other physiognomies) (adapted from Gandolfi, 2000; Whitmore, 1989). For dispersal syndrome, some species were classified into two classes, mainly due to primary and secondary syndromes. Finally, for life-form classifications, we used undershrub, shrub, treelet, tree, palm, and unclassified, and some species were also classified into two or three groups due to phenotypic variation in the study region. These mixed categories were included to avoid missing information, to allow the inclusion of contrasting classification of species in the literature, and to consider the phenotypic variation of species. We also classified species according to their threatened status, following the Brazilian National Red List (Portaria MMA 443, 2014).

Data analysesOverall plant diversity (floristic α-diversity among matrices and ages of restoration)

We chose a diversity index insensitive to the sample size to evaluate the general patterns of α-plant diversity in forest remnants in human-modified landscapes (Objective 1). Since the abundance of individuals of the superior and inferior strata varied considerably among plots, we used the Fisher alpha diversity index (Fisher et al., 1943) to describe biodiversity among plots, forest remnants, landscapes, and also to compare diversity with landscape structure and dynamic metrics (see next section). This index is insensitive to sample size (Magurran, 1988) and it has been previously used for ecological studies in fragmented landscapes (Benchimol and Peres, 2015; Laurance et al., 2001). Low values of Fisher's alpha indicate that few species are added as the sampling size increases, which means there is a higher dominance by common species. When the alpha value is high, rare species are more abundant, and increasing the sampling size adds more species (Magurran, 1988). Thus, to freely compare and describe plant diversity, we performed the Fisher alpha diversity index for all 60 plots considering superior and inferior strata together and separately.

To assess whether superior and inferior strata were correlated to each other, we performed the Mantel test considering the Euclidean distance after a Hellinger transformation, which is appropriated to reduce the importance of large abundances (Borcard et al., 2011). We tested the significance of the correlations through 999 random permutations (Legendre and Legendre, 2012). If the two distance matrices are correlated, this means that the results found for the superior stratum would be expected to be the same as those found for the inferior stratum. On the other hand, if the strata are not correlated to each other, this means the ecological processes may act differently on each separate stratum.

Considering that most forest remnants in the landscapes experienced forest suppression and/or regeneration, we assessed differences in floristic diversity among groups of regeneration. We assembled the plots into four groups according to their time of regeneration (time since abandonment from previous anthropic uses), based on aerial photographs and satellite images used to calculate metrics of landscape dynamic: (1) 20 plots with at least 30 years of regeneration (i.e., plots that were already forests in 1978, and for analyses purposes, we included nine plots in this group that were already classified as forests in 1962); (2) 17 plots with at least 13 years of regeneration (i.e., plots that were first classified as forests in 1995); (3) 17 plots with at least eight years of regeneration (i.e., plots that were first classified as forests in 2000); and finally, (4) five plots that are Eucalyptus plantations that have been abandoned and now present an understory of native species regeneration, but still cannot be considered forests. Subsequently, we only used 59 plots in the analysis, since one plot was excluded due to missing landscape data from Ferraz et al. (2014). To test the difference in diversity and total richness between the type of matrix and years, we used Fisher’s alpha diversity index and the total number of spec. All these response variables were not spatially autocorrelated (Supplementary Material - Table 2), providing the independence assumption for statistical analyses per plot, analyzing superior and inferior strata together and separately. We constructed boxplots to visualize the differences between diversities of superior and inferior strata among different ages of regeneration and among agricultural matrices (pastureland and sugarcane). We performed a two-way ANOVA to test the differences in Fisher’s diversity and total richness between the type of matrix and years of regeneration, considering a significance level of 5%.

Partition of diversity

To analyze hierarchical patterns of diversity at different scales (Objective 2), we performed a partition of diversity to determine which scales of our sampling contributes more to the total diversity. In other words, we investigated whether the total diversity is an outcome of differentiation among the matrices, the landscapes, the forest remnants or the plots themselves. To accomplish this, we analyzed our data according to sampling design in four distinct categories: (A) 59 plots, (B) 59 plots arranged in 12 forest remnants, (C) 12 forest remnants arranged in six landscapes, and (D) six landscapes arranged in two matrices. Thus, we performed an Additive Partitioning of diversity using two metrics: (I) the Shannon Wiener index, a diversity metric that accounts for both abundance and evenness of the species present in the community (Fisher’s alpha diversity index was not available in the R package that was used to perform this analysis), and (II) the species richness, that corresponds to the species recorded in the community. Available in the “vegan” package (Oksanen et al., 2016) in R software (R Core Team, 2016), this analysis decomposes gamma (γ) diversity into alpha (α) and beta (β) components that are expressed in the same units (Crist et al., 2003). Thus, the total species diversity and the total richness (γ) found for the entire region is partitioned into the average diversity and richness that occur for each plot (α) and among plots (β), so the contributions of α and β to total diversity can be compared across the hierarchical sampling design (Crist et al., 2003). Thus, we tested whether the gamma diversity and richness (γ) are randomly decomposed into alpha (each plot) and among three levels of beta diversity: β1 – among forest remnants; β2 – among landscapes; and β3 – among matrices. We randomized 999 times to generate a distribution of expected values for each diversity component and richness according to the null model. Null models are typically indicated for analyzing different components of diversity by comparing the observed similarity in community composition to that which is expected by a random assembly (Chase et al., 2011). This approach is especially recommended when changes in beta diversity are simultaneous with differences in α and/or γ diversity, as in the case of habitat loss (Püttker et al., 2015). Specifically, p < 0.05 indicates that alpha and beta diversity or richness is higher or lower than expected according to the null hypothesis, which means that there are differences in diversity or richness between plots (αplot), among forest remnants (βrem), among landscapes (βlan) or among matrices (βmat). Bar charts were plotted with the proportion of contribution to gamma diversity and richness of each component observed against the expected. This analysis was performed on all species together and separately in both the superior and inferior strata.

Influence of landscape structure and dynamic metrics on species distribution

We applied a Canonical Redundancy Analysis (RDA; Rao, 1964) to assess the influence of landscape structure and dynamics on plant richness and diversity in order to understand the role of present and past events on current plant assembly (Objective 3). An RDA is a constrained ordination method that preserves the Euclidean distances among sites in the full-dimensional space. This analysis was performed at the plot level, and we used 58 plots as two were excluded due to a lack of landscape information. We considered only species recorded in more than one plot in the analysis to diminish noise created by rare species, which could mask patterns (Ter Braak, 1995). We analyzed a matrix of 197 species, since 100 of the 297 original species were recorded at one site only, against the 58 plots considering species of the two strata together. For the superior stratum species, we analyzed 99 of the original 161 species, and for inferior stratum, 160 of the original 238 species were analyzed. For longer gradients or landscapes, as in our case study, many species are replaced by others along the gradient, which generates many zeros in the species data matrix (Legendre and Gallagher, 2001) even after rare species are excluded. To minimize this issue, we performed a resemblance matrix among sites using Hellinger’s transformation of species data (Legendre and Gallagher, 2001), which resulted in a matrix of distance that is appropriate for the analysis of community composition (Legendre and Legendre, 1998). This noise reduction in the analysis also reduces the effects of differences in sample size among superior and inferior strata when analyzed together. In the environmental matrix used as a constraining matrix for the RDA, we used the landscape structure, dynamic metrics, and elevation for each plot obtained from the digital elevation model ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer), whose image was acquired in October 2011. Since the environmental matrix has different units, we also performed a standardization on the environmental matrix by scaling each datum to zero mean and unit variance. Both transformations were realized using the deconstand function in R. There were no significant autocorrelations among environmental variables after the Spearman test, considering α = 0.05. We then applied a Monte Carlo test with 999 random permutations to estimate the significance of the relationship between the floristic data set and the environmental variables, assuming H0 as a non-linear relationship between the matrices. To perform these analyses, we used the “vegan” package (Oksanen et al., 2016) in R software (R Core Team, 2016).

ResultsOverall plant diversity (floristic alpha diversity among matrices and ages of restoration)

In total, we sampled 7686 individuals classified into 297 species; of which 161 species and 1803 individuals were sampled in the superior stratum and 5883 individuals from 238 species were in the inferior stratum. Among the individuals from the superior stratum, six species were not confirmed, and two were only classified until the genus level. Among species from the inferior stratum, two were not confirmed, eight were only classified until the genus level, and nine were only classified to the family level.

Fabaceae and Meliaceae were the most abundant families, each with more than 1000 individuals sampled. The most abundant species in the inferior stratum were the shrubs Hybanthus atropurpureus (A.St.-Hil.) Taub. and Acnistus arborescens (L.) Schltdl., and in the superior stratum, the trees Trichilia clausseni C.DC., Eugenia florida DC. and Cupania vernalis Cambess. According to the Brazilian National Red List (Portaria MMA 443, 2014), six of the species sampled are classified as “nearly threatened” (Aspidosperma polyneuron Müll.Arg.; Balfourodendron riedelianum (Engl.) Engl.; Handroanthus impetiginosus (Mart. ex DC.) Mattos; Ocotea puberula (Rich.) Nees; Xylopia aromatica (Lam.) Mart.), two are considered vulnerable (Cedrela fissilis Vell.; Cedrela odorata L.), and one is threatened (Cariniana legalis (Mart.) Kuntze) and all were sampled with less than five individuals.

From more than half the species (171), ten or less individuals were sampled in all plots, and from 51 species, only one individual was sampled, most classified as non-pioneer of canopy trees. Most species were classified from shrub to tree (84%) and mainly classified as non-pioneer of canopy (39%), understory (24%), and pioneer (17%). Almost 65% of species were classified as zoochoric (including a few species also classified as autochoric and anemochoric) and 21% were classified as anemochoric (including a few species also classified as autochoric).

According to the Mantel test, the similarity among plots was correlated between the superior and inferior stratum (r = 0.38, p = 0.001). In other words, the similarity in species composition found among plots for superior stratum species was related in 38% to that found for plants recorded in the inferior stratum. This correlation may indicate that the plant composition of the inferior stratum is partially constituted of regenerating species from the superior stratum. Alternatively, the drivers that act to assemble plant species communities may act equally in both strata.

Regarding the agricultural matrix, pastureland was richer and more diverse (247 species; Fisher alpha = 480.7) than sugarcane (192 species; Fisher alpha = 365.9). In the landscape, the inferior stratum was more diverse than the superior (Fisher alpha = 488.1 and 407.4, respectively), and species density was 111 and 661 species/ha, respectively. At the landscape level, the total density found was 1252 superior stratum ind.ha−1 and 16 341 inferior stratum ind. ha-1. At the plot level, the diversity for superior stratum (Fisher alpha) ranged from 2.23 to 23.15 (8.02 on average). For the inferior stratum, diversity ranged from 2.83 to 22.42 (9.10 on average; Fig. 2). Considering only the superior stratum, we found a difference in diversity between years of regeneration (F = 2.855, p = 0.046), but not for the matrix (F = 3.636, p = 0.062) or their interaction (F = 0.464, p = 0.631). For the inferior stratum, diversity varied among agricultural matrices (F = 6.116, p = 0.0168), but not for the years of regeneration (F = 1.137, p = 0.343) or their interaction (F = 2.037, p = 0.141). Finally, for the superior and inferior strata together, the same results were found: diversity varies between pastureland and sugarcane (F = 4.251, p = 0.044), but no difference was found considering the years of regeneration (F = 1.451, p = 0.238) or interaction between the two variables (F = 0.882, p = 0.42). For species richness, there was no difference in any analyses. Overall, the diversity is higher for forests in the pastureland matrix than in the sugarcane matrix (except for superior stratum), higher in eight years of regeneration than in 30 years, and lower in 13 years of regeneration.

Fig. 2.

Plant diversity among agricultural matrices (pasturelands and sugarcane), years of regeneration (time since abandonment), and abandoned Eucalyptus with native species’ understory, on 60 sampling plots on the Corumbataí River basin, Brazil. A: superior stratum; B: inferior stratum; and C: both strata together.

Partition of diversity

The total diversity and richness differed significantly from the null model for all alpha and beta components, considering the two strata together or separately (Fig. 3 and Supplementary Material - Table 3). More than 40% of total diversity was due to the alpha component alone in all cases (for both strata together, it was close to 60%), but it was significantly lower than expected. Beta components among forest remnants (βrem), among landscapes (βlan) and matrices (βmat) were higher than expected for superior and inferior strata, together or separately. However, βrem was particularly more important for total diversity, comprising 18% for both strata together, 23% for the superior stratum, and 21% for the inferior stratum. Although the portion of total diversity at the local component (plot level) was high, this component was lower than expected according to the null models in all cases. Thus, the differences among forest remnants, among landscapes, and among matrices are more important (significantly greater than expected) for total diversity than the α component. For total species richness, less than 40% was due to α and βrem together in all cases, as both were significantly lower than expected, which contrasted with total diversity. On the other hand, βlan and βmat were higher than expected, and together they were responsible for more than 60% of the total species richness of the superior and inferior strata (together or separately). According to the null models, the differences among the matrices contributed more to total richness than the other components.

Fig. 3.

Additive partition of the total plant species diversity (Shannon index) and the total species richness among 59 plots on the Corumbataí River basin, Brazil. Values are expressed as the percent of the total diversity and total richness explained by each hierarchical level of sampling: αplot — Alpha diversity (plot level); βrem — Beta diversity among forest remnants; βlan — Beta diversity among landscapes; and βmat — Beta diversity among sugarcane and pastureland matrices. The observed partitions are compared to the expected values obtained through randomization.

Influence of landscape structure and dynamic metrics on species distribution

Overall, forest increment under the natural regeneration process has been superior to forest loss in plots for both matrices, indicated by the forest change rate (q%). Although this increment was generally higher for pastureland (Fig. 4), sugarcane plots were predominantly less dynamic over time since some plots were kept without forest suppression or natural regeneration. Processes of deforestation and natural regeneration of the forests were slightly more frequently observed in recent years for the sugarcane matrix indicated by the negative values of FCCP (Fig. 4).

Fig. 4.

Distribution of sampling plots according to Forest Change Profile Curve (FCCP) and Annual Forest Change rate (q). Circles represent plots embedded in sugarcane landscapes, and triangles represent plots embedded in pasturelands, and their size is related to the diversity index. Adapted from Ferraz et al. (2014).


Considering all species together in the superior and inferior strata (Fig. 5), the total inertia was 0.785, the inertia of the constrained axes was 0.148, and the R2 was 18.84. In other words, this value is the proportion of total variance in plant composition explained by landscape variables. Also, the relationship between the floristic data set and the landscape variables was significant (F = 1.422, p = 0.001). Plots were grouped mostly among landscapes in a pattern that plots embedded in pastureland matrices, center to right in ordination, and were positively influenced by Elevation, Mean Forest Age (FA), and Mean Annual Forest Change rate (q). Adversely, plots embedded in sugarcane matrices were positively influenced by Local Forest Neighborhood Dominance (FORNEIGH), Forest Change Curve Profile (FCCP), and Forest Contiguity (FORCONT). Despite this general pattern, few plots of S1 and P3 landscapes were mixed among all of the plots.

Fig. 5.

Relationship between landscape structure and dynamic metrics on plant species composition, in 58 sampling plots on the Corumbataí River basin, Brazil. FA: Mean Forest Age; FCCP: Forest Change Curve Profile; FORCOVER: Local Forest Cover; FORCONT: Forest Contiguity; FORNEIGH: Local Forest Neighborhood Dominance; PROX: Forest Proximity; q: Mean Annual Forest Change rate. S1, S2, S3: Studied landscapes in sugarcane matrix; P1, P2, P3: Studied landscapes in pasturelands.


We found similar results for the superior and inferior stratum analyzed separately with those found for both strata together (Fig. 5). For superior stratum species, the total inertia was 0.838 and for constrained axes it was 0.159 with R² = 18.57. These values for the inferior stratum were 0.805, 0.147, and 18.30, respectively. Also, the relationship between the floristic data set and the landscape variables was significant for both strata (F = 1.398 and 1.373, respectively, p = 0.001).


Previous studies have shown, and our results support, that there is an important plant diversity and richness that remains in the agricultural and deforested landscapes with recurrent impacts on biodiversity (Magnago et al., 2014; Machado et al., 2016; Farah et al., 2017). Plant diversity differed among agricultural matrices, probably as a response to the different disturbance regimes of pasturelands and sugarcane fields. Forest remnants in pasturelands were overall more dynamic, where deforestation and regeneration occurred more frequently over the past few years (Ferraz et al., 2014), and therefore, the highest habitat heterogeneity could also explain the highest diversity (Ricklefs, 1977). The inferior stratum was richer and more diverse than the superior, due to the presence of plants from different growth forms on the understory (shrubs and small trees). Furthermore, in the inferior stratum, diversity varied among agricultural matrices, being higher in the pastureland matrix. Thus, the matrix is likely to be important to forest dynamics. For example, the matrix will often act as a selective filter for species’ movement across the landscape. The type of vegetation, height, or management (e.g., use of herbicides, pesticides, fire and fencing to restrict livestock transit) in the matrix will determine the permeability level for movements of individuals (da Silva et al., 2015). Characteristics of the surrounding matrix can strongly influence the plant regeneration patterns in forest remnants (Grau, 2004). The presence of isolated tress and the scattered distribution of forest remnants, common in pasturelands, are proven to facilitate seed dispersal (Chazdon and Guariguata, 2016; Holl et al., 2016; Toh et al., 1999). Our results support that pastureland matrices seem to provide better conditions for plant regeneration in similar forest edges when compared to sugarcane matrices.

Although individuals from superior stratum are present in lower densities when compared to other fragments of Tropical Semi-deciduous Forests (Ivanauskas et al., 1999), this species pool and high individual density per hectare suggests that these forest remnants still present resilience and therefore may be able to persist in the landscape. The high number of animal-dispersed plant species, as adults and regenerating, suggest that seed dispersal processes are still happening in such anthropic landscapes. However, the local population persistence of such tree species is firmly related to the effective conservation of their frugivorous bird communities responsible for complementary functional processes (Rother et al., 2016). The process of defaunation of large-bodied dispersers observed in human-modified landscapes remains one of the main issues to be addressed and limits the potential persistence of large-seeded plant species in the study sites and other tropical regions (Cramer et al., 2007; Dirzo et al., 2014; Melo et al., 2006).

Previous studies have also shown a great influence of agricultural uses on plant diversity, ecological processes such as seed dispersal and recruitment, and ultimately, forest resilience (Herrera et al., 2010; Jakovac et al., 2015). It is important to recognize that our plots were located near forest edges, where these impacts are stronger, and diversity is expected to be low (Murcia, 1995). Our results showed that diversity is low at the plot level, which might be explained by a myriad of factors that influence plant diversity at different scales. At the plot level, there are natural processes such as environmental filtering (soil conditions and light availability) and ecological filtering (demographic stochasticity, seed dispersal limitations, competition, recruitment, etc.) that have a considerable influence on selecting species and defining biodiversity (Arroyo-Rodríguez et al., 2015; Wills et al., 2006). In degraded environments, these processes are strengthened by anthropogenic disturbances, such as edge effects, competition with hyper abundant native species, mainly vines in our case (Schnitzer and Carson, 2010; Tabarelli et al., 2012), and regenerating species selection due to cattle pressure (Griscom et al., 2009). These disturbances have been shown to decrease diversity, thus leading to biotic homogenization at the local level (Püttker et al., 2015). Disturbances at plot level are even more important as our results showed that plots with eight years of regeneration are more diverse than plots with 13 or 30 years of abandonment. Although previous studies show an increase in forest cover in recent years for other tropical regions along the Corumbataí River basin (Ferraz et al., 2014; Molin et al., 2017), this coincides with a decrease in forest quality, because old-growth forests are being replaced by secondary forests in many stages of succession. This result can be explained by two observations. Firstly, there is an urgent need for the conservation and restoration of old-growth forest patches in human-modified landscapes. These patches are irreplaceable in biodiversity conservation (Gibson et al., 2011), particularly in the protection of groups of species sensitive to fragmentation and anthropogenic disturbances, such as large-seeded species (Melo et al., 2007), endemic species (Liebsch et al., 2008) and rare species (Mouillot et al., 2013). The removal of degradation drivers and restoration of these patches could reverse the trajectory of degradation and increase their potential for housing such sensitive groups (Vidal et al., 2016). Secondly, this result reinforces the importance of forest regeneration for conserving biodiversity at the landscape scale (Latawiec et al., 2016). In our study sites, young forest regeneration accounted for the greatest diversity within communities at the plot level. Natural regeneration has also been included as one of the main methods to restore degraded and open areas in large-scale projects, which would increase forest cover in highly fragmented landscapes (Chazdon and Guariguata, 2016). Together with old-growth forests (and not by replacing them), young forest patches and disturbed patches increase habitat cover for native species, increase habitat heterogeneity at the landscape level, and enhance landscape connectivity, which enhances the potential for biodiversity conservation of human-modified landscapes (Tabarelli et al., 2010).

It is interesting to point out that different matrix types may influence the forest fragment edge(Laurance and Yensen, 1991; Mesquita et al., 1999). The higher diversity observed in eight year plots in the pastureland matrices shows that matrix quality clearly affects plant dynamics and dispersal across the landscape. Pasturelands are widely considered to be a hostile mode of land use (Peres et al., 2010). However, sugarcane is often considered even more detrimental due to the high level of degradation caused by management practices, such as burning the adjacent forests, since sugarcane leaf biomass was typically burned before harvesting in the dry season up until recently (Martinelli and Filoso, 2008). This could explain the differences we observed for species diversity in eight year plots compared to 13 and 30 year plots occurring mainly in sugarcane matrices. Betadiversity was higher than expected by chance for total diversity among plots, among time of regeneration and among matrices, and in most cases for species richness (except among plots). This is consistent with previous studies in other tropical regions (Arroyo-Rodríguez et al., 2013; Machado et al., 2016). In our study area, the impermanence of secondary forests together with recurrent human pressure, implies different successional pathways for such plant communities. Each trajectory implies different forest assemblages, which allows the accumulation of a higher number of species at the landscape level (Arroyo-Rodríguez et al., 2015; Püttker et al., 2015). Essentially, as proposed by Tscharntke et al. (2012), the landscape-moderated dissimilarity of local communities (alpha) determines landscape-wide biodiversity (beta) and overrides any negative local effects of habitat fragmentation on biodiversity.

This understanding of biodiversity in fragmented habitats has a substantial influence on management strategies and the conservation role of forests under the distinct age of regeneration in human-modified landscapes. Previous studies have shown, and our results also support, that these disturbed forests have low plant diversity and richness, especially when compared to large tracts of pristine forests (Santo-Silva et al., 2015). On the other hand, all remaining and regenerating forests together still preserve high diversity at the landscape level (i.e., high beta diversity) even when compared to large conservation units (Farah et al., 2017), and not only for plants but also for other biological groups, such as large-bodied mammals (Beca et al., 2017). In other words, in human-modified landscapes, each forest regardless of its size has an immense value to the conservation of overall biodiversity, and the discussion regarding the conservation role of secondary forests should be migrated from the local to the landscape level (Arroyo-Rodríguez et al., 2015, 2009; Bergamin et al., 2017). Our results reinforce the need for secondary forests to be prioritized in future conservation and restoration efforts at the landscape scale to improve their persistence (Reid et al., 2018).

Notwithstanding, it is not often clear what the best strategy is to conserve remnant biodiversity in highly modified landscapes. Previous studies have reported that in landscapes with a low amount of forest patches, increasing connectivity is ineffective, whereas at low connectivity increasing high patch numbers is useless (Johst et al., 2011). Others have found that 30% of habitat cover is the turning point for species occurrence and the persistence of functional groups sensitive to fragmentation (Estavillo et al., 2013; Pardini et al., 2010). Therefore, landscapes below this threshold of 30% should not be targeted as priority areas for restoration, due to the low cost-effectiveness of these projects (Banks-Leite et al., 2014; Tambosi and Metzger, 2013). Nonetheless, each biological group has a different response to fragmentation thresholds and landscape structure, and for plants, regenerating individuals tend to reflect the current quality and configuration of the landscape, while adults are able to remain in the landscape for a very long time following habitat loss (Rigueira et al., 2013).

In the Atlantic Forest, 88% of forest cover is present in landscapes below the theoretical biodiversity thresholds (Banks-Leite et al., 2014), and therefore would be in danger of not being targeted for restoration programs. In addition, landscapes with different anthropic occupations present different distributions of patches and different historical dynamics. In our study sites, the influence of landscape structure and dynamic metrics on plant diversity differed among agricultural matrices, suggesting that, in highly fragmented landscapes, only the percentage of forest cover is not enough to define restoration and conservation strategies. In our study region, sugarcane landscapes present larger and less connected forest remnants, whereas, in pastureland landscapes, forest remnants are mostly connected but usually smaller and linear (Ferraz et al., 2014). Therefore, different strategies should be applied to each landscape to conserve biodiversity. In landscapes with low forest cover and highly dynamic histories, it might be more effective to increase biodiversity persistence by decreasing patch destruction than by increasing patch creation (Johst et al., 2011), so it is crucial that the remaining biodiversity should be maintained through conservation and restoration of current forests. This strategy must be considered to implement new protected areas, therefore increasing the number of conservation units, along with the assistance of financial incentive programs and specific legislation since biodiversity is spread in forest remnants in the landscape (Bergamin et al., 2017). Although nowadays there are legal instruments that can increase landscape connectivity through the restoration of riparian areas, very few require an increase in quality of the actual remaining forests and prevention of within-forest disturbances (Brancalion et al., 2012; Viani et al., 2015; Barlow et al., 2016; Rother et al., 2018).

Final considerations

Our study showed the high potential for biodiversity conservation of secondary forests within private properties in highly modified landscapes. We found plant diversity to be higher than expected in both strata, but local plant communities differed from forests and agricultural matrices. Therefore, to conserve biodiversity in highly deforested landscapes, conservation plans should be enacted at the regional level and not at the forest patch level. Since every forest remnant presents a distinctive group of species that must be conserved and they play a unique role in biodiversity conservation.

Our results underscore the importance of assessing the structure and function of secondary forests, but we still have a limited understanding of how to provide more options for the persistence of populations of threatened and endangered species in agricultural landscapes. It is also important to continue research on how we can restore these forest remnants in order to improve their role in biodiversity conservation.


We thank the students and partners of the Laboratory of Ecological Restoration (LERF/ESALQ/USP) for helping us in fieldwork and discussion of the data. We are indebted to the Brazilian National Research Council (CNPq) for the productivity fellowship (grant #308503/2014-7), to KMPMBF and the research grant for JRSAM (#140825/2013-4). The authors also thank CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior do Governo do Brasil, #88882.305844/2018-01) for the fellowship grant to DCR. This research was supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo (Biota Project #2013/50718-5 and #1999/09635-0), FAPESP (#2011/06782-5) and CNPq (#561897/2010-7).

[Arroyo-Rodríguez et al., 2009]
V. Arroyo-Rodríguez, E. Pineda, F. Escobar, J. Benítez-Malvido.
Value of small patches in the conservation of plant-species diversity in highly fragmented rainforest.
Conserv. Biol., 23 (2009), pp. 729-739
[Arroyo-Rodríguez et al., 2013]
V. Arroyo-Rodríguez, M. Rös, F. Escobar, et al.
Plant β-diversity in fragmented rain forests: testing floristic homogenization and differentiation hypotheses.
J. Ecol., 101 (2013), pp. 1449-1458
[Arroyo-Rodríguez et al., 2015]
V. Arroyo-Rodríguez, F.P.L. Melo, M. Martínez-Ramos, et al.
Multiple successional pathways in human-modified tropical landscapes: new insights from forest succession, forest fragmentation and landscape ecology research.
[Banks-Leite et al., 2014]
C. Banks-Leite, R. Pardini, L.R. Tambosi, et al.
Using ecological thresholds to evaluate the costs and benefits of set-asides in a biodiversity hotspot.
Science, (2014), pp. 1041
[Barlow et al., 2016]
J. Barlow, G.D. Lennox, J. Ferreira, et al.
Anthropogenic disturbance in tropical forests can double biodiversity loss from deforestation.
Nature, 535 (2016), pp. 144-147
[Beca et al., 2017]
G. Beca, M.H. Vancine, C.S. Carvalho, et al.
High mammal species turnover in forest patches immersed in biofuel plantations.
Biol. Conserv., 210 (2017), pp. 352-359
[Benchimol and Peres, 2015]
M. Benchimol, C.A. Peres.
Edge-mediated compositional and functional decay of tree assemblages in Amazonian forest islands after 26 years of isolation.
J. Ecol., 103 (2015), pp. 408-420
[Benchimol et al., 2017]
M. Benchimol, D.C. Talora, E. Mariano-Neto, T.L.S. Oliveira, A. Leal, M.S. Mielke, D. Faria.
Losing our palms: the influence of landscape-scale deforestation on Arecaceae diversity in the Atlantic forest.
For. Ecol. Manage., 384 (2017), pp. 314-322
[Bergamin et al., 2017]
R.S. Bergamin, V.A.G. Bastazini, E. Vélez-Martin, et al.
Linking beta diversity patterns to protected areas: lessons from the Brazilian Atlantic Rainforest.
Biodivers. Conserv., 26 (2017), pp. 1557-1568
[Bivand and Wong, 2018]
R.S. Bivand, D.W.S. Wong.
Comparing implementations of global and local indicators of spatial association.
Test, 27 (2018), pp. 716-748
[Borcard et al., 2011]
D. Borcard, F. Gillet, P. Legendre.
Numerical Ecology with R.
Springer, (2011),
[Brancalion et al., 2012]
P.H.S. Brancalion, R.A.G. Viani, R.R. Rodrigues, R.G. César.
Estratégias para auxiliar na conservação de florestas tropicais secundárias inseridas em paisagens alteradas.
Boletim do Museu Paraense Emílio Goeldi Ciências Naturais, 7 (2012), pp. 219-234
[Carneiro et al., 2016]
M.S. Carneiro, C.C.F. Campos, F.N. Ramos, F.A.M. Dos Santos.
Spatial species turnover maintains high diversities in a tree assemblage of a fragmented tropical landscape.
Ecosphere, 7 (2016), pp. 1-12
[Cassiano et al., 2013]
C.C. Cassiano, S. Frosini, D.B. Ferraz, et al.
Spatial assessment of water-related ecosystem services to prioritize restoration of forest patches.
Nat. Conserv., 11 (2013), pp. 176-180
[Chase et al., 2011]
J.M. Chase, N.J.B. Kraft, K.G. Smith, M. Vellend, B.D. Inouye.
Using null models to disentangle variation in community dissimilarity from variation in α-diversity.
Ecosphere, 2 (2011), pp. 1-11
[Chazdon and Guariguata, 2016]
R.L. Chazdon, M.R. Guariguata.
Natural regeneration as a tool for large-scale forest restoration in the tropics: prospects and challenges.
Biotropica, 48 (2016), pp. 716-730
[Chazdon et al., 2009]
R.L. Chazdon, C.A. Peres, D. Dent, et al.
The potential for species conservation in tropical secondary forests.
Conserv. Biol., 23 (2009), pp. 1406-1417
[Cramer et al., 2007]
J.M. Cramer, R.C.G. Mesquita, G. Bruce Williamson.
Forest fragmentation differentially affects seed dispersal of large and small-seeded tropical trees.
Biol. Conserv., 137 (2007), pp. 415-423
[Crist et al., 2003]
T.O. Crist, J.A. Veech, J.C. Gering, K.S. Summerville.
Partitioning species diversity across landscapes and regions: a hierarchical analysis of alpha, beta, and gamma diversity.
Am. Nat., 162 (2003), pp. 734-743
[da Silva et al., 2015]
L.G. da Silva, M.C. Ribeiro, É. Hasui, C.A. da Costa, R.G. da Cunha.
Patch size, functional isolation, visibility and matrix permeability influences neotropical primate occurrence within highly fragmented landscapes.
[Diniz-Filho et al., 2012]
J.A.F. Diniz-Filho, T. Siqueira, A.A. Padial, T.F. Rangel, V.L. Landeiro, L.M. Bini.
Spatial autocorrelation analysis allows disentangling the balance between neutral and niche processes in metacommunities.
Oikos, 121 (2012), pp. 201-210
[Dirzo et al., 2014]
R. Dirzo, H.S. Young, M. Galetti, et al.
Defaunation in the anthropocene.
Science, 345 (2014), pp. 401-406
[Dormann et al., 2013]
C.F. Dormann, J. Elith, S. Bacher, C. Buchmann, G. Carl, G. Carré, J.R.G. Marquéz, B. Gruber, B. Lafourcade, P.J. Leitão, T. Münkemüller, C. McClean, P.E. Osborne, B. Reineking, B. Schröder, A.K. Skidmore, D. Zurell, S. Lautenbach.
Collinearity: a review of methods to deal with it and a simulation study evaluating their performance.
Ecography, 36 (2013), pp. 27-46
[Estavillo et al., 2013]
C. Estavillo, R. Pardini, P.L.B. Rocha.
Forest loss and the biodiversity threshold: an evaluation considering species habitat requirements and the use of matrix habitats.
[Fahrig, 2003]
L. Fahrig.
Effects of Habitat Fragmentation on Biodiversity.
Annu. Rev. Ecol. Evol. Syst., 34 (2003), pp. 487-515
[FAO, 1995]
Forest Resources Assessment 1990: Global Synthesis.
Food and Agriculture Organization of the United Nations, (1995),
[Farah et al., 2017]
F.T. Farah, R. Muylaert, L. De, M.C. Ribeiro, J.W. Ribeiro, J.R. Mangueira, S.A. De, V.C. Souza, R.R. Rodrigues.
Integrating plant richness in forest patches can rescue overall biodiversity in human-modified landscapes.
For. Ecol. Manage., 397 (2017), pp. 78-88
[Ferraz et al., 2014]
S.F.B. Ferraz, K.M.P.M.B. Ferraz, C.C. Cassiano, et al.
How good are tropical forest patches for ecosystem services provisioning?.
Landsc. Ecol., 29 (2014), pp. 187-200
[Fisher et al., 1943]
R.A. Fisher, A.S. Corbet, C.B. Williams.
The relation between the number of species and the number of individuals in a random sample of an animal population.
J. Anim. Ecol., 12 (1943), pp. 42-58
[Gandolfi, 2000]
S. Gandolfi.
História natural de uma Floresta Estacional Semidecidual no município de Campinas.
Universidade Estadual de Campinas — UNICAMP, (2000),
[Garcia et al., 2006]
G.J. Garcia, S.L. Antonello, M.G.M. Magalhães.
The environmental atlas of the corumbatai watershed – SP, Brazil.
Revista Brasileira de Cartografia, 58 (2006), pp. 73-79
[Gardner et al., 2009]
T.A. Gardner, J. Barlow, R. Chazdon, et al.
Prospects for tropical forest biodiversity in a human-modified world.
Ecol. Lett., 12 (2009), pp. 561-582
[Gibson et al., 2011]
L. Gibson, T.M. Lee, L.P. Koh, et al.
Primary forests are irreplaceable for sustaining tropical biodiversity.
Nature, 478 (2011), pp. 378-381
[Grau, 2004]
H.R. Grau.
Landscape context and local-scale environment influence regeneration of pioneer tree species in treefall gaps.
Biotropica, 36 (2004), pp. 52-59
[Griscom et al., 2009]
H.P. Griscom, B.W. Griscom, M.S. Ashton.
Forest regeneration from pasture in the dry tropics of Panama: effects of cattle, exotic grass, and forested Riparia.
Restor. Ecol., 17 (2009), pp. 117-126
[Herrera et al., 2010]
J.M. Herrera, D. García, J.M. Morales.
Matrix effects on plant-frugivore and plant-predator interactions in forest fragments.
Landsc. Ecol., 26 (2010), pp. 125-135
[Hobbs and Huenneke, 1992]
R.J. Hobbs, L.F. Huenneke.
Disturbance, diversity, and invasion: implications for conservation.
Conserv. Biol., 6 (1992), pp. 324-337
[Holl et al., 2016]
K.D. Holl, J.L. Reid, J.M. Chaves-Fallas, F. Oviedo-Brenes, R.A. Zahawi.
Local tropical forest restoration strategies affect tree recruitment more strongly than does landscape forest cover.
J. Appl. Ecol., 54 (2016), pp. 1091-1099
[Ivanauskas et al., 1999]
N.M. Ivanauskas, R.R. Rodrigues, A.G. Nave.
Fitossociologia de um trecho de Floresta Estacional Semidecidual em Itatinga, São Paulo, Brasil.
Scientia Florestalis, 56 (1999), pp. 83-99
[Jakovac et al., 2015]
C.C. Jakovac, M. Peña-Claros, T.W. Kuyper, F. Bongers.
Loss of secondary-forest resilience by land-use intensification in the Amazon.
J. Ecol., 103 (2015), pp. 67-77
[Johst et al., 2011]
K. Johst, M. Drechsler, A.J.A.Van Teeffelen, et al.
Biodiversity conservation in dynamic landscapes: trade-offs between number, connectivity and turnover of habitat patches.
J. Appl. Ecol., 48 (2011), pp. 1227-1235
[Latawiec et al., 2016]
A.E. Latawiec, R. Crouzeilles, P.H.S. Brancalion, et al.
Natural regeneration and biodiversity: a global meta-analysis and implications for spatial planning.
Biotropica, 48 (2016), pp. 844-855
[Laurance and Yensen, 1991]
W.F. Laurance, E. Yensen.
Predicting the impacts of edge effects in fragmented habitats.
Biol. Conserv., 55 (1991), pp. 77-92
[Laurance et al., 2001]
W.F. Laurance, S. D’Angelo, D. Pérez-Salicrup, et al.
Rain forest fragmentation and the structure of amazonian liana communities.
Ecology, 82 (2001), pp. 105-116
[Laurance et al., 2007]
W.F. Laurance, H.E. Nascimento, S.G. Laurance, A. Andrade, R.M. Ewers, K.E. Harms, R.C. Luizão, J.E. Ribeiro.
Habitat fragmentation, variable edge effects, and the landscape-divergence hypothesis.
[Legendre and Gallagher, 2001]
P. Legendre, E. Gallagher.
Ecologically meaningful transformations for ordination of species data.
Oecologia, 129 (2001), pp. 271-280
[Legendre and Legendre, 1998]
P. Legendre, L. Legendre.
Numerical Ecology.
2nd ed., Elsevier B.V., (1998),
[Legendre and Legendre, 2012]
P. Legendre, L. Legendre.
Numerical Ecology.
3rd ed., Elsevier B.V., (2012),
[Liebsch et al., 2008]
D. Liebsch, M. Marques, R. Goldenberg.
How long does the Atlantic Rain Forest take to recover after a disturbance? Changes in species composition and ecological features during secondary succession.
Biol. Conserv., 141 (2008), pp. 1717-1725
[Lima and Mariano-Neto, 2014]
M.M. Lima, E. Mariano-Neto.
Extinction thresholds for Sapotaceae due to forest cover in Atlantic Forest landscapes.
For. Ecol. Manage., 312 (2014), pp. 260-270
[Lima et al., 2015]
R.A.F. Lima, D.P. Mori, G. Pitta, et al.
How much do we know about the endangered Atlantic Forest? Reviewing nearly 70 years of information on tree community surveys.
Biodivers. Conserv., 24 (2015), pp. 2135-2148
[Lôbo et al., 2011]
D. Lôbo, T. Leão, F.P.L. Melo, A.M.M. Santos, M. Tabarelli.
Forest fragmentation drives Atlantic forest of northeastern Brazil to biotic homogenization.
Divers. Distrib., 17 (2011), pp. 287-296
[Machado et al., 2016]
F.S. Machado, M.A.L. Fontes, R.M. Santos, P.O. Garcia, C.L. Farrapo.
Tree diversity of small forest fragments in ecotonal regions: why must these fragments be preserved?.
Biodivers. Conserv., 25 (2016), pp. 525-537
[Magnago et al., 2014]
L.F.S. Magnago, D.P. Edwards, F.A. Edwards, A. Magrach, S.V. Martins, W.F. Laurance.
Functional attributes change but functional richness is unchanged after fragmentation of Brazilian Atlantic forests.
J. Ecol., 102 (2014), pp. 475-485
[Magurran, 1988]
A.E. Magurran.
Ecological Diversity and Its Measurement.
1° ed., University Press, (1988),
[Martinelli and Filoso, 2008]
L.A. Martinelli, S. Filoso.
Expansion of sugarcane ethanol production in Brazil: environmental and social challenges.
Ecol. Appl., 18 (2008), pp. 885-898
[Melo et al., 2006]
F.P.L. Melo, R. Dirzo, M. Tabarelli.
Biased seed rain in forest edges: evidence from the Brazilian Atlantic forest.
Biol. Conserv., 132 (2006), pp. 50-60
[Melo et al., 2007]
F.P.L. Melo, D. Lemire, M. Tabarelli.
Extirpation of large-seeded seedlings from the edge of a large Brazilian Atlantic forest fragment.
Ecoscience, 14 (2007), pp. 124-129
[Melo et al., 2013]
F.P.L. Melo, V. Arroyo- Rodríguez, L. Fahrig, M. Martınez-Ramos, M. Tabarelli.
On the hope for biodiversity-friendly tropical landscapes.
Trends Ecol. Evol., 28 (2013), pp. 462-468
[Mesquita et al., 1999]
R.C. Mesquita, P. Delamônica, W.F. Laurance.
Effect of surrounding vegetation on edge-related tree mortality in Amazonian forest fragments.
Biol. Conserv., 91 (1999), pp. 129-134
[Molin et al., 2017]
P.G. Molin, S.E. Gergel, B.S. Soares-Filho, S.F.B. Ferraz.
Spatial determinants of Atlantic Forest loss and recovery in Brazil.
Landsc. Ecol., 32 (2017), pp. 857-870
[Metzger et al., 2009]
J.P. Metzger, A.C. Martensen, M. Dixo, et al.
Time-lag in biological responses to landscape changes in a highly dynamic Atlantic forest region.
Biol. Conserv., 142 (2009), pp. 1166-1177
[Molin et al., 2018]
P.G. Molin, R. Chazdon, S.F.B. Ferraz, P.H. Brancalion.
A landscape approach for cost-effective large-scale forest restoration.
J. Appl. Ecol., 55 (2018), pp. 2767-2778
[Morante-Filho et al., 2015]
J. Morante-Filho, V. Arroyo-Rodríguez, D. Faria.
Patterns and predictors of beta diversity in the fragmented Brazilian Atlantic forest: a multiscale analysis of forest specialist and generalist.
[Mouillot et al., 2013]
D. Mouillot, D.R. Bellwood, C. Baraloto, et al.
Rare species support vulnerable functions in high-diversity ecosystems.
[Murcia, 1995]
C. Murcia.
Edge effects in fragmented forests: implications for conservation.
Trends Ecol. Evol., 10 (1995), pp. 58-62
[Oksanen et al., 2016]
J.F. Oksanen, G. Blanchet, R. Kindt, P. Legendre, P.R. Minchin, R.B. O’hara, G.L. Simpson, P. Solymos, M.H.H. Stevens, H. Wagner.
Vegan: Community Ecology Package. R Package Version 2.3-5.
[Pardini et al., 2010]
R. Pardini, A.D.A. Bueno, T.A. Gardner, P.I. Prado, J.P. Metzger.
Beyond the fragmentation threshold hypothesis: regime shifts in biodiversity across fragmented landscapes.
[Peres et al., 2010]
C.A. Peres, T.A. Gardner, J. Barlow, J. Jansen, F. Michalski, A.C. Lees, I.C.G. Vieira, F.M.D. Moreira, K. Feeley.
Biodiversity conservation in human-modified Amazonian forest landscapes.
Biol. Conserv., 143 (2010), pp. 2314-2327
[Püttker et al., 2011]
T. Püttker, A.A. Bueno, C.D.S. De Barros, S. Sommer, R. Pardini.
Immigration rates in fragmented landscapes–empirical evidence for the importance of habitat amount for species persistence.
[Püttker et al., 2015]
T. Püttker, Ade Arruda Bueno, P.I. Prado, R. Pardini.
Ecological filtering or random extinction? Beta-diversity patterns and the importance of niche-based and neutral processes following habitat loss.
Oikos, 124 (2015), pp. 206-215
[R Core Team, 2016]
R Core Team.
R: A Language and Environment for Statistical Computing.
R Foundation for Statistical Computing, (2016),
[Rao, 1964]
C.R. Rao.
The use and interpretation of principal component analysis in applied research.
Sankhya: Indian J. Stat. Ser. A, 26 (1964), pp. 329-358
[Reid et al., 2018]
J.L. Reid, R.A. Zahawi, M.E. Fagan, J. Lucas, J. Slaughter.
The ephemerality of secondary forests in southern Costa Rica.
Conserv. Lett., e12607 (2018), pp. 1-7
[Ribeiro et al., 2009]
M.C. Ribeiro, J.P. Metzger, A.C. Martensen, F.J. Ponzoni, M.M. Hirota.
The Brazilian Atlantic Forest: how much is left, and how is the remaining forest distributed? Implications for conservation.
Biol. Conserv., 142 (2009), pp. 1141-1153
[Ricklefs, 1977]
R.E. Ricklefs.
Environmental heterogeneity and plant species diversity: a hypothesis.
Am. Nat., 111 (1977), pp. 376-381
[Rigueira et al., 2013]
D.M.G. Rigueira, P.L.B. Rocha, E. Mariano-Neto.
Forest cover, extinction thresholds and time lags in woody plants (Myrtaceae) in the Brazilian Atlantic Forest: resources for conservation.
Biodivers. Conserv., 22 (2013), pp. 3141-3163
[Rother et al., 2016]
D.C. Rother, M.A. Pizo, P. Jordano.
Variation in seed dispersal effectiveness: the redundancy of consequences in diversified tropical frugivore assemblages.
Oikos, 125 (2016), pp. 336-342
[Rother et al., 2018]
D.C. Rother, C.Y. Vidal, I.C. Fagundes, M. Silva, S. Gandolfi, R.R. Rodrigues, G. Nave, R.A.G. Viani, P.H.S. Brancalion.
How legal-oriented restoration programs enhance landscape connectivity? Insights from the brazilian Atlantic Forest.
Trop. Conserv. Sci., 11 (2018), pp. 1-9
[Santo-Silva et al., 2015]
E.E. Santo-Silva, K.D. Withey, W.R. Almeida, et al.
Seedling assemblages and the alternative successional pathways experienced by Atlantic forest fragments.
Plant Ecol. Diversi., (2015), pp. 1-10
[Schnitzer and Carson, 2010]
S.A. Schnitzer, W.P. Carson.
Trepadeiras suppress tree regeneration and diversity in treefall gaps.
Ecol. Lett., 13 (2010), pp. 849-857
[Sfair et al., 2016]
J.C. Sfair, V. Arroyo-Rodríguez, B.A. Santos, M. Tabarelli.
Taxonomic and functional divergence of tree assemblages in a fragmented tropical forest.
Ecol. Appl., 26 (2016), pp. 1816-1826
[Tabarelli et al., 2010]
M. Tabarelli, A.V. Aguiar, M.C. Ribeiro, J.P. Metzger, C.A. Peres.
Prospects for biodiversity conservation in the Atlantic Forest: lessons from aging human-modified landscapes.
Biol. Conserv., 143 (2010), pp. 2328-2340
[Tabarelli et al., 2012]
M. Tabarelli, B.A. Santos, V. Arroyo-Rodríguez, F.P.L. Melo.
Secondary forests as biodiversity repositories in human-modified landscapes: insights from the Neotropics.
Boletim do Museu Paraense Emílio Goeldi Ciências Naturais, 7 (2012), pp. 319-328
[Tambosi and Metzger, 2013]
L.R. Tambosi, J.P. Metzger.
A framework for setting local restoration priorities based on landscape context.
Nat. Conserv., 11 (2013), pp. 152-157
[Ter Braak, 1995]
C.J.F. Ter Braak.
Data Analysis in Community and Landscape Ecology, pp. 91-173
[Tilman et al., 1994]
D. Tilman, R.M. May, C.L. Lehman, M.A. Nowak.
Habitat destruction and the extinction debt.
Nature, 371 (1994), pp. 1
[Toh et al., 1999]
I. Toh, M. Gillespie, L. Lamb.
The role of isolated trees in facilitating tree seedling recruitment at a degraded sub-tropical rainforest site.
Restor. Ecol., 7 (1999), pp. 288-297
[Tscharntke and Brandl, 2004]
T. Tscharntke, R. Brandl.
Plant-insect interactions in fragmented landscapes.
Annu. Rev. Entomol., 49 (2004), pp. 405-430
[Tscharntke et al., 2012]
T. Tscharntke, J.M. Tylianakis, T.A. Rand, et al.
Landscape moderation of biodiversity patterns and processes — eight hypotheses.
Biol. Rev., 87 (2012), pp. 661-685
[Valente and Vettorazzi, 2005]
R.De O.A. Valente, C.A. Vettorazzi.
Avaliação da estrutura florestal na bacia hidrográfica do Rio Corumbataí, SP.
Scientia Florestalis, 68 (2005), pp. 45-57
[Van der Pijl, 1982]
L. Van der Pijl.
Principles of Dispersal in Higher Plants.
2nd ed., Springer-Verlag, (1982),
[Viani et al., 2015]
R.A.G. Viani, F.N.A. Mello, I.E. Chi, P.H.S. Brancalion.
A new focus for ecological restoration: management of degraded forest remnants in fragmented landscapes.
GPL news, 12 (2015), pp. 5-9
[Vidal et al., 2016]
C.Y. Vidal, J.R. Mangueira, F.T. Farah, D.C. Rother, R.R. Rodrigues.
Biodiversity conservation of forests and their ecological restoration in highly-modified landscapes.
Biodiversity in Agricultural Landscapes of Southeastern Brazil,
[Whitmore, 1989]
T.C. Whitmore.
Canopy gaps and the two major groups of forest trees.
Ecology, 70 (1989), pp. 536-538
[Wills et al., 2006]
C. Wills, K.E. Harms, R. Condit, et al.
Nonrandom processes maintain diversity in tropical forests.
Science, 311 (2006), pp. 527-531
[Wright, 2010]
S.J. Wright.
The future of tropical forests.
Ann. N. Y. Acad. Sci., 1195 (2010), pp. 1-27
Copyright © 2021. Associação Brasileira de Ciência Ecológica e Conservação
Perspectives in Ecology and Conservation

Subscribe to our newsletter

Article options