The relationship between forest fragmentation and the potential geographical distribution of psittacids (Psittaciformes: Psittacidae) in the Yucatán Peninsula, Mexico

Parrots are a seriously threatened group because of the illegal pet trade, their habitat fragmentation and the destruction of tropical forests. This study aimed to determine the relationship between forest fragmentation and the geographic potential distribution of parrot species in the Yucatan Peninsula, Mexico. The research was carried out in the Mexican portion of the Yucatan Peninsula Biotic region, which includes the states of Campeche, Quintana Roo and Yucatan. We worked with the eight species of psittacids that inhabit the selected areas: Eupsittula nana, Amazona albifrons, Amazona xantholora, Amazona autumnalis, Pionus senilis, Pyrilia haematotis, Amazona oratrix and Amazona farinosa. For these species, we used the potential geographical distribution maps that were obtained through the Maximum Entropy algorithm, published in 2014. To measure the levels of forest fragmentation in the Peninsula, we employed the land use and vegetation maps from series IV (2007-2010) of the National Institute of Statistics and Geography. We calculated different landscape metrics at 100 sites randomly distributed within the predicted presence and absence areas for the species. Multivariate methods were used to shorten dimensionality in the analysis, as well as to explore the observed patterns. We did not find any pattern regarding the number of variables that contribute to the logistic regression models (LR) and the size of the parrots’ potential geographical distribution area within the Peninsula. For Pyrilia haematotis and Amazona oratrix, the fragmentation variables used did not seem to have any relationship with their potential geographical distribution in the Peninsula (intersection, Pr<0.0001). Only for Eupsittula nana, Amazona xantholora and Pionus senilis, we found that the mean patch edge, the mean shape index and the mean perimeterarea ratio were important to determine their potential distribution patterns, respectively; and, in general, the variables indicated the preferences of these parrots to occupy larger habitat patches. Parrots must be given a special care and protection within the states of Campeche and Quintana Roo, since there are areas with high parrot diversity and are importantly affected by human activities. The best recommendation for managers and users of protected areas within the Yucatan Peninsula is that extensive areas of forest should be maintained in order to guarantee the suvival of psittacid populations. Rev. Biol. Trop. 65 (4): 1470-1482. Epub 2017 December 01.

Forest fragmentation is one of the main anthropogenic factors that threaten bird conservation (Bregman, Sekercioglu, & Tobias, 2014).In addition to habitat loss, a reduction in the mean size of the remaining forest patches can also act as a limiting factor on bird populations (Davis, 2004).Although numerous studies have been carried out on fragmentation, these tend to be difficult to interpret since they have been conducted at the patch, but not at the landscape scale (Fahrig, 2003).The landscape, from the ecological point of view, is a heterogeneous surface that contains a mosaic of patches or land use types (Arroyo & Fahrig, 2014).Fragmentation is a process that occurs on the landscape scale (McGarigal & Cushman, 2002), thus, when studies are conducted at the smaller patch level, it is difficult to make interpretations at the larger scale, at which fragmentation is taking place (Brennan, Bender, Contreras, & Fahrig, 2002).
Responses to fragmentation vary between different bird groups or species (Visco et al., 2015;Uezu & Metzger, 2016), and forest bird responses to fragmentation may be highly species/guild-specific (Lees & Peres, 2009).For example, most frugivore species disappear from forest fragments after their isolation and few persist in small fragments surrounded by secondary growth forest (Cintra, Magnusson, & Albernaz, 2013).Some landscape attributes (e.g.composition and spatial arrangement) and behavioral attributes of the species (e.g.mobility or dispersal ability) may interact to determine the species response to fragmentation (Major, Christie, & Gowing, 2001).Moreover, birds with different sensitivity to fragmentation display contrasting responses to landscape dynamics (Metzger et al., 2009;Uezu & Metzger, 2016).It has also been found that, some effects in bird species sensitive to fragmentation, have a 25 year time-lag to landscape changes, reinforcing the idea that it is not enough to conserve the present landscape structure to maintain those species in the long term (Uezu & Metzger, 2016).On the other hand, relative effects of habitat amount and fragmentation are scale dependent, and the effects of habitat fragmentation depend on the size of the landscape considered (Smith, Fahrig, & Francis, 2011).
Among forest birds, psittacids are seriously threatened due to the illegal pet trade (Marín et al., 2012) as well as the destruction and fragmentation of their habitats (Snyder, McGowan, Gilardi, & Grajal, 2000).Nonetheless, parrots are highly mobile species able to move among forest patches (Renton, 2001), therefore studies of fragmentation effects on birds at the patch level may not be applicable to parrots.However, parrots can be affected by changes that occur in the characteristics of their habitats and their extension, even if they can move among distant forest patches (Tella et al., 2016).Different species of parrots and parakeets select areas with a high proportion of forest, and avoid areas with low forest cover (Legault, Chartendrault, Theuerkauf, Rouys, & Barre, 2011).Nevertheless, other species are capable of exploiting some of the resources present in the landscape matrix, surrounding the forest fragments, and even degraded environments, such as agro-ecosystems, which can support viable populations of certain psittacine species (Bucher, 1992;Galindo, 1999;Romero, Naranjo, Morales, & Nigh, 2006).
Identifying how forest fragmentation alters the distribution of different psittacid species is essential for the implementation of measures that are effective in mitigating its consequences.Research needs to be centered on ascertaining how these species respond to different levels of fragmentation, and which variables affect them the most.Parrots must be associated with large tropical forest areas and for the majority, the presence and distribution probability are associated with the most preserved areas characterized by little anthropogenic alteration.However, for some parrot species, habitat fragmentation can have positive effects at one scale and negative effects in another.
Eight species of psittacids inhabit the Yucatan Peninsula, Mexico, (MacKinnon, 2005), and several of these maintain abundant populations (e.g., Galindo, 1999;Macías & Iñigo, 2003;Plasencia & Escalona, 2014a).However, in this region of Southeast Mexico, the intensification of forest degradation and fragmentation (Turner II, Geoghegan, & Foster, 2004) is threatening the viability of parrot populations.In the Peninsula, it is predicted that different levels of fragmentation cannot affect the distribution of parrot species with wide environmental tolerance (e.g.Amazona albifrons or Eupsittula nana), while the most affected will be those species with less environmental tolerance (e.g.Amazona farinosa or Pyrilia haematotis).Within this context, this study objective focused on determining the relationship between forest fragmentation and patterns of potential geographic distribution of psittacids inhabiting the Yucatan Peninsula.

Study area:
This study was carried out in the Mexican portion of the Yucatan Peninsula Biotic region which includes the states of Campeche, Quintana Roo and Yucatán (17°0'-21°35' N & 86°43'-92°25' W) (Barrera, 1962).The vegetation of the Yucatan Peninsula changes gradually along an environmental gradient associated with rainfall distribution patterns that stretch from the semiarid Northwest to the wetter Southeast.The changes in vegetation largely reflect rainfall patterns, with dry low deciduous scrub and semi-deciduous forest dominating towards the Northwest, and semiwet and tropical semi-evergreen forest predominant towards the Southeast (Pennington & Sarukhán, 1998).By the 1990s, more than 60 % of the Yucatan Peninsula was modified as a consequence of land use alterations, which led to widespread changes in vegetation distribution patterns (Carranza, Barreto, Molina, Bezaury, & Loreto, 1996).The disappearance and fragmentation of forest areas within the Yucatan Peninsula continues, principally in the Southwest and Southeast regions (Turner II et al., 2004).The study was made between late June 2014 and September 2016, using database analyzes, without fieldwork.For analyzes, potential geographical distribution maps, land use and vegetation maps were used.

Study species and potential geographical distribution maps:
The eight species of psittacids that inhabit the Yucatan Peninsula were studied: Olive-throated Parakeet (Eupsittula nana), White-fronted Amazon (Amazona albifrons), Yellow-lored Amazon (Amazona xantholora), Red-lored Amazon (Amazona autumnalis), White-crowned Parrot (Pionus senilis), Brown-hooded Parrot (Pyrilia haematotis), Yellow-headed Amazon (Amazona oratrix), and Southern Mealy Amazon (Amazona farinosa) (MacKinnon, 2005).Within the Yucatan Peninsula, the greatest concentration of species occurred in the Southern area near the borders with Guatemala and Belize (Plasencia & Escalona, 2014a).Amazona albifrons and E. nana are present almost in all Yucatan Peninsula.Amazona oratrix was only found in the Southwest portion of the Yucatan Peninsula.Amazona xantholora was present around much of the Peninsula, including the Cozumel Island, but it was absent in areas where A. oratrix was found.Amazona autumnalis, P. haematotis, P. senilis and A. farinosa were only located in wetter areas of the Southern portion of the Yucatan Peninsula.
There are no real distribution maps of the present psittacid species in the Yucatan Peninsula.Therefore, the models of potential geographical distribution published in 2014 were used (Plasencia, Escalona, & Esparza, 2014a;Plasencia & Escalona, 2014b).For A. xantholora and A. oratrix, the potential geographical distribution maps obtained by Plasencia et al. (2014a) were used.The distribution maps obtained by Plasencia and Escalona (2014b) were used for the remaining six species.
Forest fragmentation in Yucatan Peninsula: A total of 15 digital files containing land use and vegetation maps from series IV of INEGI (National Institute of Statistics and Geography) were used to measure levels of forest fragmentation within potential psittacid distribution areas in the Yucatan Peninsula (INEGI, 2010).These maps cover the Mexican portion of the Yucatan Peninsula and were obtained for the period 2007-2010(INEGI, 2010)).The selected forest map included mangroves, petenes (islands of arboreal vegetation associated with freshwater springs), tropical forests and oak forest as these vegetation types are inhabited most frequently by parrot species.The remaining vegetation formations and land uses were eliminated from the study.All of the aforementioned forest types were grouped together, resulting in a map that represented the forest masses of the Yucatan Peninsula.
All the above procedures were performed using the Arc View 3.2 software (ESRI, 1999).
One hundred hexagonal plots, each with 587 km², were randomly distributed onto the potential geographical distribution map of each species.The hexagonal shape minimized the relationship area/perimeter without leaving any spaces (Mateucci & Silva, 2005).The loss of vegetation cover and the processes of fragmentation had varied according to the particular region or state within the Yucatan Peninsula (Roy & Turner II, 2006;Daniels, Painter, & Southworth, 2008).Therefore the 100 plots covered a sufficiently large area (approximately 43 % of total area of the Peninsula) as to represent the different levels of forest fragmentation within the Yucatan Peninsula.For each species, plots with 50 % or more of the potential geographical distribution area occupied by one of the parrot species were classified as presence and the remaining as absence.The number of plots classified as absence (A) or presence (P) for each species was the following: A. albifrons (five A, 95P), E. nana (seven A, 93P), A. xantholora (19A, 81P), A. autumnalis (64A, 36P), P. haematotis (57A, 43P), P. senilis (66A, 34P), A. farinosa (75A, 25P) and A. oratrix (98A, two P).
Each one of the 100 plots that was classified as presence (P) or absence (A) for each parrot species, was used to superimpose onto the forest mass map; therefore, the same plots were obtained, but this time included information on the forest mass.These procedures were performed with the ArcView 3.2 Geoprocessing Wizard Extension (ESRI, 1999).A series of indexes were calculated for each plot using the ArcView 3.2 Patch Analyst Extension (ESRI, 1999).The indexes were selected to describe landscape fragmentation, characterizing the entire landscape as suggested by Fahrig (2003).We calculated the following variables: total forest area (TFA), number of patches (NP), mean patch size (MPS), patch size coefficient of variance (PSCV), patch size standard deviation (PSSD), total edge (TE), edge density (ED), mean patch edge (MPE), mean shape index (MSI), area-weighted mean shape index (AWMSI), mean perimeter-area ratio (MPAR), mean patch fractal dimension (MPFD), area-weighted mean patch fractal dimension (AWMPFD), Shannon's diversity index (SDI) and Shannon's evenness index (SEI) (see Elkie, Rempel, & Carr, 1999).
Furthermore, the distance between patches (DP) and their spatial configuration were calculated.The shortest distance between the edge of one patch to another, and the distance between the centroids of all the patches were calculated for each plot.The sum of the crossed multiples of the matrixes of distances between the patch edges and their centroids was then calculated to obtain a distance summary variable.In order to acquire a measure of the spatial configuration of patches in the landscape, using the matrix of distances between the edges of the patches in each plot, the Moran's I Index (four and eight neighbors) (Moran, 1950) and the Geary's C index (four and eight neighbors) (Geary, 1954) were calculated by applying the ROOKCASE complement for Excel 97/2000 (Sawada, 1999), which determined the global and local spatial autocorrelation.When calculating the spatial autocorrelation, the dispersion level of the patches in each plot was determined, thus obtaining an approximation of how the patches were distributed throughout the landscape and whether any relationships or patterns existed in this distribution.
A Pearson's correlation was performed in order to reduce the collinearity between the fragmentation and spatial configuration variables used.Those variables with coefficients of |r| < 0.7 were selected (Dormann et al., 2012).Between the pairs of highly correlated variables, based on the natural history, habitat and knowledge obtained through field observations of these species, one of the variables that possibly has a strong relationship with psittacine distribution patterns in the Yucatan Peninsula was selected.
The four variables obtained to define the spatial configuration of patches was highlycorrelated (r > 0.7).Only Moran's I Index with four neighbors was selected from these variables; this index was the most appropriate when carrying out analysis at the landscape level, as it produced the best measure of spatial autocorrelation at the global level.Moran's I Index values were not obtained for 12 % of the plots given that they only contained one patch.A total of 51 % and 29 % of the plots, obtained positive and negative values, respectively, using Moran's I Index, indicating that at the scale of this study, there was a trend towards spatial clustering of patches that were found at similar distances from each other.The summary variable for distances was not correlated with any other variables (r≤0.4;p>0.05).Among the fragmentation indexes, only the mean shape index, mean perimeter-area ratio and the mean patch fractal dimension were not highly correlated.
Once the highly-correlated variables were eliminated in order to reduce dimensionality in the analysis, a Principal Component Analysis (PCA) was applied.Out of the 21 initial variables, only 11 remained for the PCA: total forest area, number of patches, mean patch size, patch size standard deviation, mean patch edge, mean shape index, mean perimeter-area ratio, mean patch fractal dimension, Shannon's evenness index, distances between patches and Moran's I Index with four neighbors.The Varimax rotation method was implemented to minimize the number of variables and simplify their interpretation; subsequently, the components that explained 50 % variability were selected.Variables that presented a factorial coefficient of more than 0.7 were then selected from these components.Subsequent to the Varimax rotation, the four first components explained 50 % of the variability.With six variables, factorial coefficients greater than 0.7 were obtained for the first four components: total forest area, patch size standard deviation, mean patch edge, mean shape index, mean perimeter-area ratio and distance between patches.
A logistic regression (LR) was implemented to determine the possible relationship between the matrix of fragmentation and spatial configuration variables and the potential geographical distribution of the psittacid species (presence/absence maps).When applying the LR, a Logit type model with a binary response and 95 % confidence intervals was included.For best model selection, the Akaike information criteria (AIC) was used (Akaike, 1983) with a minimum and maximum of two and six variables, respectively.One hundred iterations with a convergence level of 0.000001 were carried out and second level interactions between variables were considered.In order to assess how good or bad these models were, the Chi 2 test, associated with the Log ratio, was calculated.The Hosmer-Lemeshow test was used as an additional mean for evaluating the goodness of fit of the logistic regression models.For all the above analysis, the Excel XLSTAT (2013) complement was implemented.

RESULTS
No pattern was observed regarding the number of variables provided by the model and the potential geographical distribution area of the eight studied parrots in the Yucatan Peninsula.For the eight species, the best LR model was the one that presented the lowest value according to the Akaike Information Criteria (AIC).The number of resolved variables fluctuated between four and six (Table 1).In the case of P. haematotis and A. oratrix, the variables did not contribute to the regression model, and the Chi 2 value for the intersection was significant (Table 2).At the scale used in this study, the potential geographical distribution of both of these species (P.haematotis and A. oratrix) did not appear to show any relationship with the fragmentation variable; however, for the remaining six species, the Chi 2 test demonstrated this relationship to be significant and different variables contributed significantly to the model (Table 2).But the "odds ratio" values were close to one for most of these variables or their interaction, thus any changes they experienced would have little effect on determining the presence or absence of some of these remaining six psittacids species.Overall, although the interactions between variables appeared to be more important, the "odds ratio" values indicated the contrary.
The distance between patches was the only variable that did not significantly influence the distribution of any of the species.The mean perimeter-area ratio contributed most frequently to parrot distribution; however, this variable only significantly influenced the distribution of P. senilis (Table 2).
Regarding the quality of the logistic regression models, in order to assess how good or bad these models were, one of the most important results is Chi 2 associated with the Log ratio (Table 3).The fragmentation variables provided significant information for explaining the variability related to the presence/absence of six of the species (Table 3).The Hosmer-Lemeshow test proved that these results reaffirm another means of evaluating the goodness of fit of the logistic regression models.With this test, high Chi² values (p<0.0001) were only obtained for A. oratrix and P. haematotis, which in this case signifies that for both of these species, the model does not present a good fit.
contributed was MSI (Table 2).The MSI evaluates the shape complexity of a fragment.When a negative value of the coefficient for this variable is obtained (Table 2), the probability of distribution of A. xantholora is higher in areas with more regular forms, which often coincided with less fragmented sites.For P. senilis, the most related fragmentation variable to its potential distribution was MPAR (Table 2).A decrease in the fragments size is associated with an increase of the relationship perimeter/ area.When a negative value of the coefficient for MPAR is obtained, it can be interpreted that the potential distribution of P. senilis is favored in areas with larger fragments of forest.

DISCUSSION
When considering the scale at which this study was carried out, the selected fragmentation variables do not appear to be related to the potential distribution patterns of most parrot species found in the Yucatan Peninsula.This was not expected, as the fragmentation of the forest ecosystems has increased over the entire Peninsula, predominantly due to the loss of large forest tracts (Turner II et al., 2004).It has also been demonstrated, for some psittacid species, that ecosystem modifications caused by anthropogenic activities can actually be beneficial, as they often provide new sources of more accessible food (Sánchez, Ballari, Bucher, & Masello, 2016).In the Calakmul region of the state of Campeche, large flocks of A. albifrons have been observed in perturbed areas; besides, this species is perceived as a pest by local peasant farmers, due to the damage they are thought to inflict on their cornfields (Galindo, 1999;Escalona & Velasco, 2000).However, we cannot ignore that this only occurs for a little number of species, because parrots are among the most threatened bird orders with 28 % of extant species classified as threatened under IUCN criteria (Olah et al., 2016).
Although different reports mention that fragmentation affects psittacid populations (e.g.Snyder et al., 2000;Forshaw, 2006), there are very few studies that state which variables Only for E. nana, A. xantholora and P. senilis some variables were important in determining the distribution patterns (Table 2) and in general, showed that these parrot species had a preference to occupy sites that conserve large patches of vegetation.For E. nana, the variable that contributed the most was MPE (Table 2).The MPE gives us a measure of the degree of fragmentation and connectivity of the landscape and, the fact of obtaining positive values for the coefficient of this variable, means that the probability of distribution of E. nana increases as its values are higher.These higher values of MPE indicated that the areas in which E. nana species are most likely to be distributed, are less fragmented.In the case of A. xantholora, the variable that most are the most important or at what scale do they most frequently affect parrot species (e.g., Renton & Salinas, 2004;Miller, Bianchi, Mullins, & Haig, 2013;Plasencia, Escalona, & Esparza, 2014b).Within the Yucatan Peninsula, it has been speculated that the scarlet macaw (Ara macao) had become locally extinct due to the accelerated destruction and fragmentation of high evergreen rainforest (Berlanga, Wood, Salgado, & Figueroa, 2000); however, there is no empirical evidence that validates this.Most research has concentrated on habitat use or population dynamics of parrots in fragmented landscapes (e.g., Conti & Galetti, 2007;Monterrubio, Ortega, & Marín, 2009;Carneiro, Jiménez, & White, 2012); therefore, since there is no solid basis to this question regarding psittacids, comparing results is problematic.
In the Yucatan Peninsula, the continual and rapid loss of forest areas (Vester et al., 2007) can affect the quality and quantity of potential habitats for Psittacidae species.The forests that have experienced the biggest vegetation cover loss are located in the Southeast and Southwest areas of the states of Quintana Roo and Campeche, respectively (Turner II et al., 2004); and most importantly, it is in these forests, where the highest diversity of parrot species is located (Plasencia & Escalona, 2014a).Some of the species are only distributed in the Southern Yucatan Peninsula because of its humid climate, and the presence of humid to moist forest types.This vegetation loss is causing an increase in isolation, reduction and degradation of the patches of remnant forests in the Yucatan Peninsula.Parrots in these areas may be affected mainly by the reduction and degradation of these limited remnant forest, since they do not fulfill all their requirements to stay and nest.Nevertheless, since parrots have high mobility, they are able to disperse over large areas, and may endure an increase in forest patch isolation (Munshi & Wilkinson, 2006).On the other side, these landscape modifications have a negative influence on food and nesting sites requirements (Renton & Salinas, 2004), and these have to be considered in future research efforts in the area.
Amazona oratrix is one of the most endangered parrot species in Mexico (Cantú, Sánchez, Grosselet, & Silva, 2007), but none of the selected fragmentation variables were related to its potential distribution.In Mexico, one of the most prized species for the pet trade is A. oratrix; in the Yucatan Peninsula this parrot is only found in the Southwest, one of the most fragmented areas of the region, as a consequence of extensive livestock farming (Villalobos & Mendoza, 2010).However, none of the selected fragmentation variables appear to be related to the potential distribution of this species, possibly due to the very small number of plots where this parrot occurs, considering the reduced area they inhabit and the scale of this study.Furthermore, Macías and Iñigo (2003) and Plasencia and Escalona (2014a) demonstrated that in the Yucatan Peninsula, this species presents one of the highest values of relative abundance in Mexico, indicating its capacity to survive even in perturbed areas.Nevertheless, the majority of parrot species presents high longevity (Munshi & Wilkinson, 2006); therefore, changes in population numbers, attributable to the effects of environmental modifications such as fragmentation, only become evident over the long-term.
Before drawing conclusions, it should be taken into account that, as A. oratrix was only represented in a small number of plots, in which the levels of fragmentation were high, it was not possible to establish a comparison between areas with different characteristics.If the analysis is only based on areas with similar characteristics, it is difficult to determine whether fragmentation may or may not influence the distribution of the species.It would be necessary to expand the study area for this species, in order to obtain a larger number of plots on which to perform the analysis.Although, as far as the distribution of A. oratrix is known in Mexico, the remnant areas where it is still found are mostly fragmented and degraded habitats (Macías & Iñigo, 2003), what keeps making this type of analysis difficult.
Eupsittula nana, P. senilis and A. xantholora presented a relationship between some fragmentation variables and their potential distribution, indicating a preference by these species for larger patches of habitat, a characteristic that has also been determined for other species of birds (Davis, 2004).Even though, psittacids have the ability to adapt to anthropic habitat modifications, many of these species are normally associated with areas dominated by conserved forest that provide them with suitable sites for feeding and reproduction (Morales, 2005;Monterrubio & Escalante, 2006).
In the case of A. xantholora, one of the fragmentation variables demonstrated a marked influence on the potential distribution of this species (mean shape index), indicating that at the scale used, this particular species is more sensitive to the fragmentation process.The results determined that the mean shape index negatively affects the presence of this species in the Yucatan Peninsula, indicating that with an increase in patch irregularity there is an associated decrease in the probability of finding this species.In different areas subject to high levels of anthropic pressure throughout the Yucatan Peninsula, forest patches tend to be smaller in area and more irregular in shape.
The results obtained in this study partially support the theory which responds that fragmentation vary among different species (see Visco et al., 2015;Uezu & Metzger, 2016).When performing this type of research, it is necessary to previously determine which scale is most suitable for each species.However, the results also reinforce the fact that fragmentation is scale dependent, and the effects of habitat fragmentation depend on the size of the landscape considered (see Smith et al., 2011).Future research on the relationship between fragmentation and potential geographical distribution of different parrot species in the Yucatan Peninsula, should consider smaller scales of study and the inclusion of other environmental variables (e.g.land use, vegetation types, spatial configuration) in the analysis (e.g.Plasencia et al., 2014b).Fragmentation is a process that is directly related to the loss of vegetation cover (one of its principal causes) and land use changes (surrounding the patches), and, a particular combination of both of these factors, could provide the key to explain psittacid distribution patterns.
The most important recommendation for managers and users of protected areas within the Yucatan Peninsula is that extensive areas of forest should be maintained in order to guarantee the survival of psittacid populations.One of the most important challenges facing the conservation of parrots in the fragmented landscapes of the Yucatán Peninsula, is the existing imbalance between ecosystem transformation and the speed at which decisions are taken at the private, state and federal level.Forest areas sometimes experience extremely rapid changes, while at the same time management decisions require lengthy processes for the assignation of conservation priorities.However, as psittacids are charismatic birds, they could potentially act as "flagship species", promoting actions to preserve the habitats where they are present in the Peninsula.Moreover, protecting these psittacids would indirectly contribute to protect a multitude of other species that make up the ecological community of their habitats, and so acting as "umbrella species".

ACKNOWLEDGMENTS
We are grateful to D. Denis, K. Renton and L. G. Esparza for their helpful comments.Thanks to INEGI for providing the Land Use and Vegetation Maps of the Yucatan Peninsula.Thanks to CONACYT for the economic support provided through the grants 239499 and 21467.Thanks to the manuscript reviewers for their valuable suggestions and comments.Finally, we would like to thank Idea Wild for the donated equipment.

TABLE 1
Best logistic regression models according to the Akaike Information Criterion (AIC) and its corresponding statistics, obtained through the fragmentation variables and the parrots potential geographical distribution models in the Yucatán Peninsula, Mexico Contributions of the fragmentation variables to the parrot's potential geographical distribution in the Yucatán Peninsula, according to the logistic regression models.It just includes the intersections, interactions and independent variables with significant contributions to the models