Introduction
Large frugivorous birds that inhabit forests are often among the first species to disappear in fragmented landscapes (e.g. Willis 1979, Aleixo & Vielliard 1995, Galetti & Aleixo 1998, Anjos 2006, Lees & Peres 2008). Ecological traits and requirements such as low density, low fecundity, spatio-temporal variability of food resources, matrix avoidance and area sensitivity increase their vulnerability in fragmented landscapes (Franklin 1980, Strahl & Grajal 1991, Galetti & Pizo 1996, Lynch & Lande 1998, Ragusa-Netto 2010). However, a recent review indicated that when forest fragments differ little the impact of forest fragmentation on frugivorous birds is lower than expected (Bregman et al. 2014). Declines in populations of frugivorous birds may have implications for the biotic integrity of forest fragments due to their ecological function as seed dispersers (Schupp et al. 2002, Terborgh et al. 2002, Galetti et al. 2013). Indeed, if the persistence of frugivorous populations is higher than expected, their role as seed dispersers is crucial for certain plant species (those which depend on birds for dispersal) occupying isolated forest fragments.
Local extinction often occur after a time lag, which makes the evaluation of population persistence more difficult. In several studies, populations may persist on the brink of extinction for long periods before finally becoming extinct (Brooks et al. 1999, Helm et al. 2006, Vellend et al. 2006). This time lag in extinction is called relaxation time (Diamond 1972). Isolated patches of habitat in fragmented or degraded habitats are said to owe an extinction debt for the species for which future extinctions are inevitable (Tilman et al. 1994, Kuussaari et al. 2009). Thus, the mere presence of frugivorous bird species in a forest fragment is no guarantee that the populations will remain in the long-term (Metzger et al. 2009). In this study, we made a preliminary evaluation of the population persistence of large frugivorous birds in forest fragments based on a temporal data series.
Herein we evaluated if large frugivorous birds persisted in six forest fragments in southern Brazil over a period of 10 years. In this case we evaluated only the persistence of the species, not variations in population size. In the largest of those forest fragments, we effectively evaluated variations in the population sizes of the frugivorous birds over the period of 10 years.
Material and Methods
1. Study area
The study was carried out in the Parque Estadual Mata dos Godoy (PEMG) and five other forest fragments (FA, FB, FC, FH and FI) in the north of Paraná State, southern Brazil (Fig. 1). The region is characterized by a humid subtropical climate (Köppen Cfa) with rainfall concentrated in the summer months and annual average temperature and precipitation of 20.8 - 21.6 °C and of 1100 - 1200 mm, respectively (Maack 1981). The soils have the same parent material (basalt), resulting in fertile reddish clay soils such as eutropherric red latosols (Maack 1981).

Figure 1. Location of studied areas in Londrina, Parana, Southern Brazil. Studied fragments are presented in dark grey, fragments in light grey and agricultural areas in white. Geographical coordinates are presented in Table 01.
The fragmentation process occurred in the studied region about 50 years ago and reduced the forest to small and isolated fragments immersed in a matrix mainly composed by agriculture. In such a scenario of habitat loss, PEMG remains relatively undisturbed and is considered the largest and best preserved area of semideciduous forest in northern Paraná State with a high biological integrity (Anjos et al. 2009). The other fragments vary in size and isolation levels (Table 1). In a recent study assessing the ecological integrity of Atlantic forest remnants, Medeiros & Torezan (2013) evaluated two of the six areas sampled in this study. According to the authors, the vegetation of the PEMG is of “excellent” integrity category and the vegetation of the fragment FI was classified as “regular”. Furthermore, hunting and selective tree logging are forbidden in the region, in particular in the PEMG and in the five forest fragments located in private areas.
Table 1. Characterization of forest fragments (FF) analyzed by their geographic coordinates, area (ha), connectivity, distance to the PEMG, and proximity index, a proxy for functional connectivity.
FF | Geographic coordinates | Area (ha) | Connectivity | Distance | Proximity index |
---|---|---|---|---|---|
PEMG | 23°27′9,39-S 51°15′17,77-W | 656 | Referential | 0 | 125.96 |
FA | 23°28′8,14-S 51°14′19,85-W | 56 | FC | 500 m | 166.12 |
FB | 23°28′6,31-S 51°15′17,10-W | 25 | PEMG | 400 m | 133.89 |
FC | 23°28′31,27-S 51°15′22,17-W | 11 | FB | 1100 m | 199.17 |
FH | 23°29′17,95-S 51°11′38,35-W | 72 | - | 2400 m | 8.39 |
FI | 23°30′11,35-S 51°4′34,74-W | 184 | - | 14500 m | 499.1 |
Proximity index used a search radius of 1000 meters.
We calculated the proximity index (Uezu et al. 2008) based on a binary forest vs. matrix map obtained from SOS Mata Atlântica (2008) and used by Ribeiro et al. (2009). This index was calculated for the search radius of 1000 meters around a centroid defined within the patches by using ArcGIS (ERSI 2005) and V-LATE extension (LARG 2005). The six forest fragments vary greatly in size and functional connectivity represented by the proximity index (Table 1).
2. Field work
2.1 Previous data
Data from two different studies (Anjos 2001, Anjos et al. 2004) conducted in the six forest fragments were used to evaluate the extinction debt of medium and large frugivorous bird species. Data from both studies were obtained using point counts of unlimited distance (Blondel et al. 1970, Vielliard & Silva 1990). In these studies, the points were 100 m apart and at least 50 m from the edge of forest; the sampling time at each point was 20 minutes.
Data presented in Anjos (2001) were collected monthly from January to December 1996 in the PEMG and fragments FA, FB and FC, totaling 12 samples for each fragment. Data presented in Anjos et al. (2004) were collected monthly between September and December 1997 in the PEMG, FA, FB, and FC, and over the same months in 2000 in FH and FI, totaling four samples for each fragment.
2.2 Recent data
The field data were obtained using the methods described in Anjos et al. (2007). Point counts of fixed radius and distance sampling-line transects were sampled along four 1000 m long trails (TA, TB, TC and TD) in PEMG. For data collected using point counts of limited distance, nine points were established every 100 m on each trail in the PEMG. Thus, 36 sampling points were sampled over 10 days beginning at sunrise, from September 2010 to March 2011 and from August to November 2011. The detection radius at each point was 50 m and the starting point for each transect sample alternated each day, offering the same odds of detection for all species. The sampling period at each point was 15 minutes, thus the total sampling on each day lasted for about four hours. Point counts do not consider differences in detectability or make it possible to calculate the current densities of each species; however, this method can be used as an index to detect significant differences in the abundance of populations of a particular species (Bibby et al. 1992). At each sampling point, each pair or flock of each species was counted once, and precautions were taken not to count the same individual or group more than once (Anjos et al. 2011), making it possible to calculate the Index of Point Abundance (IPA) of each species.
For the distance sampling-line transects, the trails in PEMG were traveled at a constant speed of 2 km/h, allowing the observer to record species through auditory or visual contact at a distance of 50 m. The distance of the bird relative to the observer and the detection angle in relation to the track were measured for subsequent calculation of the relative density for each species. Ten samples per transect from September 2010 to March 2011 were collected, and from August to November 2011, beginning at sunrise to coincide with the beginning of vocal activities of diurnal birds (Gilardi & Munn 1998). In the case of the Psittacidae family, whenever a flock was sighted, the number of individuals was counted visually. Based on the monitored flocks, an average number of individuals per flock of each Psittacidae species was measured giving an estimation of the number of individuals in the PEMG.
In the five forest fragments (FA, FB, FC, FH and FI) we exclusively used the distance sampling-technique with line transects along 1000 m long trails. Each fragment was sampled over a period of 4 days between October and December 2011, to verify the presence or absence of frugivorous birds. Thus, the sampling effort was similar to the lower sampling effort of previous studies carried out 10 years ago, allowing the comparison of the occurrence data in each fragment.
3. Data analysis
Only medium (40-80 g) and large (greater than 80 g) non-passerines frugivorous birds were considered for this study (Pizo 2001). According to Anjos (2001), 24 medium and large non-passerines frugivorous bird species occur in the PEMG. In this study, these species has been treated as large frugivorous birds.
The population density per hectare in the PEMG was estimated by entering data collected in the sampling-line transect into the software DISTANCE 6.0 (Thomas et al. 2010). We also calculated IPA of each species in the PEMG (according to Anjos et al. 2010) and this value was used to analyze the persistence of these populations after 10 years. The IPA was estimated by dividing the total number of species observed by the total number of points sampled in the PEMG. Wilcoxon Signed-Rank analysis was used to compare the IPA values of the species obtained in the present study with those presented for the same species a decade earlier by Anjos et al. (2004). This analysis was carried out using the function wilcox with the R software version 3.0.2 (R Development Core Team 2013).
Results
Twenty-four species of medium and large frugivorous birds were recorded in the PEMG. Of the 20 species of frugivorous birds recorded by Anjos (2001) and Anjos et al. (2004) in the remaining fragments, 11 species disappeared in some of the forest fragments. However, the rate of species that persisted was different in the forest fragments: 87% in FA; 75% in FB; 81% in FC; 72% in FH and 59% in FI (Table 2). There were also some species that had not been recorded in the survey 10 years ago: one species in FA and FC (Aratinga auricapillus) and three species in FB (Aratinga auricapillus, Psittacara leucophthalmus and Amazona aestiva; see Table 2).
Table 2. Presence/absence of medium and large frugivorous birds species in fragments analyzed (PEMG, FA, FC, FB, FH and FI) in the present study and 10 years ago (Anjos 2001, Anjos et al. 2004) (X: present species; -: absent species). Species are in systematic and phylogenetic position following nomenclature of the South American Classification Committee of the American Ornithologists' Union (Remsen et al. 2015).
Present study | Previous study | |||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Species | PEMG | FA | FB | FC | FH | FI | PEMG | FA | FB | FC | FH | FI | ||
Cracidae | ||||||||||||||
Penelope superciliaris Temminck, 1815 | X | - | X | X | X | - | X | X | X | X | X | X | ||
Columbidae | ||||||||||||||
Claravis pretiosa (Ferrari-Pérez, 1886) | X | - | - | - | - | - | X | - | - | - | X | X | ||
Patagioenas picazuro Temminck, 1813 | X | X | X | X | X | X | X | X | X | X | X | X | ||
Patagioenas cayennensis Bonnaterre, 1792 | X | X | X | X | X | X | X | X | X | X | X | X | ||
Patagioenas plumbea Vieillot, 1818 | X | - | - | - | - | - | X | - | - | - | - | - | ||
Leptotila verreauxi Bonaparte, 1855 | X | X | X | X | X | X | X | X | X | X | X | X | ||
Leptotila rufaxilla (Richard & Bernard, 1792) | X | X | X | X | X | - | X | X | X | X | X | X | ||
Geotrygon montana (Linnaeus, 1758) | X | X | X | X | X | - | X | X | X | X | X | X | ||
Trogonidae | ||||||||||||||
Trogon surrucura Vieillot, 1817 | X | X | X | X | X | X | X | X | X | X | X | X | ||
Trogon rufus Gmelin, 1788 | X | - | - | - | - | - | X | - | - | - | - | X | ||
Momotidae | ||||||||||||||
Baryphthengus ruficapillus (Vieillot, 1818) | X | X | X | X | X | X | X | X | X | X | X | X | ||
Ramphastidae | ||||||||||||||
Ramphastos dicolorus Linnaeus, 1766 | X | - | - | - | - | - | X | X | X | X | X | - | ||
Selenidera maculirostris (Lichtenstein, 1823) | X | X | X | X | X | X | X | X | X | X | X | X | ||
Pteroglossus bailloni (Vieillot, 1819) | X | - | - | - | - | - | X | - | - | - | - | X | ||
Pteroglossus aracari (Linnaeus, 1758) | X | - | - | - | - | - | X | - | - | - | - | - | ||
Psittacidae | ||||||||||||||
Primolius maracana (Vieillot, 1816) | X | - | - | - | - | - | X | - | - | X | - | - | ||
Psittacara leucophthalmus (Müller, 1776) | X | X | X | - | X | X | X | - | X | - | X | X | ||
Aratinga auricapillus (Kuhl, 1820) | X | X | X | - | X | - | X | - | - | - | - | - | ||
Pyrrhura frontalis (Vieillot, 1818) | X | - | X | - | - | - | X | X | X | X | X | X | ||
Brotogeris tirica (Gmelin, 1788) | X | - | - | - | - | - | X | - | - | X | - | - | ||
Pionopsitta pileata (Scopoli, 1769) | X | - | - | - | - | X | X | - | - | - | - | X | ||
Pionus maximiliani (Kuhl, 1820) | X | X | X | X | X | X | X | X | X | X | X | X | ||
Amazona aestiva (Linnaeus, 1758) | X | X | - | - | X | - | X | - | X | - | X | - | ||
Triclaria malachitacea (Spix, 1824) | X | - | - | - | - | - | X | - | - | - | - | - |
The Wilcoxon Signed-Rank analysis revealed that the IPA values of the analyzed large frugivorous birds changed over 10 years (V = 285; p = 0.0001). This difference was also found for the 11 species which disappeared from the fragments, but persisted in the PEMG (V = 56; p = 0.0453), and for the remaining 13 species that persisted in the fragments, beside PEMG (V = 90; p = 0.0004). In fact, we found that the mean IPA values in 2004 were higher for all those three groups of species than in the present study. When considering all species, the mean IPA in 2004 was 0.3833 while the present mean IPA was 0.0689. When considering the 11 species that disappeared from the fragments, but persisted in the PEMG, the mean IPA in 2004 was 0.1409 while the present mean IPA was only 0.0116. Finally, when considering the remaining 13 species that persisted in the fragments, the mean IPA in 2004 was 0.5885 while the present mean IPA was only 0.1173.
It was possible to calculate the population density for nine of the 24 sampled species (Table 3; 4). The species with the highest density was Aratinga auricapillus, with 1.04 ind/ha, while the lowest was Patagioenas cayennensis, with 0.231 ind/ha.
Table 3. Index of Point Abundance (IPA) of medium and large frugivorous bird species from PEMG in the present study and 10 years ago (Anjos et al. 2004). Species are in systematic and phylogenetic position following nomenclature of the South American Classification Committee of the American Ornithologists' Union (Remsen et al. 2015).
IPA | |||
---|---|---|---|
Present study | Previous study | ||
Cracidae | |||
Penelope superciliaris Temminck, 1815 | 0.006 | 0.200 | |
Columbidae | |||
Claravis pretiosa (Ferrari-Pérez, 1886) | 0.003 | 0.000 | |
Patagioenas picazuro Temminck, 1813 | 0.264 | 1.850 | |
Patagioenas cayennensis Bonnaterre, 1792 | 0.086 | 0.450 | |
Patagioenas plumbea Vieillot, 1818 | 0.003 | 0.000 | |
Leptotila verreauxi Bonaparte, 1855 | 0.247 | 0.950 | |
Leptotila rufaxilla (Richard & Bernard, 1792) | 0.008 | 0.050 | |
Geotrygon montana (Linnaeus, 1758) | 0.031 | 0.200 | |
Trogonidae | |||
Trogon surrucura Vieillot, 1817 | 0.194 | 0.900 | |
Trogon rufus Gmelin, 1788 | 0.014 | 0.000 | |
Momotidae | |||
Baryphthengus ruficapillus (Vieillot, 1818) | 0.142 | 0.550 | |
Ramphastidae | |||
Ramphastos dicolorus Linnaeus, 1766 | 0.008 | 0.350 | |
Selenidera maculirostris (Lichtenstein, 1823) | 0.236 | 0.900 | |
Pteroglossus bailloni (Vieillot, 1819) | 0.017 | 0.000 | |
Pteroglossus aracari (Linnaeus, 1758) | 0.006 | 0.050 | |
Psittacidae | |||
Primolius maracana (Vieillot, 1816) | 0.008 | 0.200 | |
Psittacara leucophthalmus (Müller, 1776) | 0.114 | 0.200 | |
Aratinga auricapillus (Kuhl, 1820) | 0.131 | 0.750 | |
Pyrrhura frontalis (Vieillot, 1818) | 0.003 | 0.450 | |
Brotogeris tirica (Gmelin, 1788) | 0.014 | 0.000 | |
Pionopsitta pileata (Scopoli, 1769) | 0.014 | 0.150 | |
Pionus maximiliani (Kuhl, 1820) | 0.086 | 0.750 | |
Amazona aestiva (Linnaeus, 1758) | 0.017 | 0.100 | |
Triclaria malachitacea (Spix, 1824) | 0.003 | 0.150 |
Table 4. Population density per hectare (D) of nine species recorded in PEMG. For the families Columbidae, Trogonidae, Momotidae and Ramphastidae, values are number of pairs per ha, while for Psittacidae values are number of individuals per ha. Species are in systematic and phylogenetic position following nomenclature of the South American Classification Committee of the American Ornithologists' Union (Remsen et al. 2015).
Species | D | Individual/band | |
---|---|---|---|
Columbidae | |||
Patagioenas picazuro (Temminck, 1813) | 0.928 | ||
Patagioenas cayennensis (Bonnaterre, 1792) | 0.231 | ||
Leptotila verreauxi Bonaparte, 1855 | 0.439 | ||
Trogonidae | |||
Trogon surrucura Vieillot, 1817 | 0.341 | ||
Momotidae | |||
Baryphthengus ruficapillus (Vieillot, 1818) | 0.415 | ||
Ramphastidae | |||
Selenidera maculirostris (Lichtenstein, 1823) | 0.781 | ||
Psittacidae | |||
Psittacara leucophthalmus (Statius Müller, 1776) | 0.976 | 4,67 ± 3,18 (n = 33) | |
Aratinga auricapillus (Kuhl, 1820) | 1.040 | 4,28 ± 2,01 (n = 40) | |
Pionus maximiliani (Kuhl, 1820) | 0.528 | 3,00 ± 1,24 (n = 27) |
Discussion
The differences found between the values of IPA comparing our recent data with those from ten years ago in the PEMG demonstrates that the effects of forest fragmentation could still have an impact on large frugivorous birds, even in a relatively large forest fragment. Our data suggest that populations of large frugivorous birds could be declining over time in the PEMG. Therefore, the large frugivorous birds in the PEMG may have been paying an extinction debt over the last ten years (Diamond 1972, Tilman et al. 1994, Kuussaari et al. 2009).
Some studies have demonstrated that forest species with large population size in fragmented landscapes tend to persist longer in forest fragments (e.g.Purvis et al. 2000, Manne & Pimm 2001, Henle et al. 2004), probably due to the rescue effect, which increase the chances to colonize empty fragments (Brown & Kodric-Brown 1977, Anjos et al. 2011). This argument is consistent with our data on the 11 species that were not recorded in at least one of the closer fragments; those species reached the lowest IPA values in the PEMG. Therefore, low population size seems to reduce the potential of the rescue effect of large frugivorous birds from surrounding forest fragments. If it is true, those 11 species could be considered as members of the "undead" sensuSodhi et al. (2011), even in the PEMG. Interestingly, the FI fragment that presented the highest proximity index (499.1), also lost the highest number of species. This unexpected result could be explained by the fact that FI is the farthest studied fragment from the PEMG which probably acts as an important source for rescue effects in that landscape. Furthermore, all the forest fragments located within 2400 meters radius from PEMG lost fewer species (two to four species), while FI, which is located 14.500 meters far from PEMG, lost seven species.
Some studies have demonstrated the negative effect of the extinction of large frugivorous birds in forest dynamics (e.g. Silva & Tabarelli 2000, Lord et al. 2002, Schupp et al. 2002, Terborgh et al. 2002, Galetti et al. 2013). Galetti et al. (2013), in a recent study conducted in the Brazilian Atlantic forest showed that functional extinction of large frugivorous birds is associated with the consistent reduction of Euterpe edulis’ seed size resulting in several negative consequences for plant recruitment and population dynamics. Penelope superciliaris, the frugivorous bird species with the largest body mass recorded in PEMG and a major disperser of E. edulis was not found either in the smallest (FB) or in the most isolated (FI) fragments. The level of isolation and the size of the fragments, which among other factors are related to the population decline or local extinction of large frugivorous birds, may be associated with functional loss of these species and its consequences to the ecosystem.
It is important to mention, however, that time difference between the two studies are relatively short (a decade) and such fluctuations could be expected in the seasonal semidecidous forest. Indeed, we do not know about the natural fluctuations of those frugivorous populations over time. For example, if resources were scarce in those particular fragments due to climatic variations in recent years, the occurrence and abundance of the studied birds could be affected. Another aspect that could affect time lag studies is the potential alterations of the matrix habitat due to anthropogenic action. Although the landscape appeared similar during the studied period, it is unclear whether some disturbance, in particular in the matrix habitat, could have affected our data. Therefore, our results should be viewed as a preliminary investigation of the effects of fragmentation over time. Only a long term monitoring campaign over a higher number of forest fragments could uncover bird population trends in the studied landscape.
We also recorded three species of Psittacidae, Aratinga auricapillus, Psittacara leucophthalmus and Amazona aestiva in three fragments of the five analyzed fragments (FA, FB and FC) in the latest survey. Unlike the remaining frugivores, Psittacidae have a high capacity for dispersal across open spaces, using a wide area throughout the day (Gilardi & Munn 1998, Lees & Peres 2009). It is important to mention that A. auricapillus is Near Threatened according to BirdLife International (2013), which could highlight the importance of the PEMG for conservation of this particular species.
In conclusion, our study indicates that forest fragmentation could be negatively influencing the persistence of the remaining populations of large frugivorous birds in the studied fragmented landscapes. Therefore, protecting the PEMG with buffer zones for restoration, as well as increasing connectivity with surrounding forest fragments with forest corridors is crucial for maintenance of frugivorous populations in the studied landscape.