Habitat characterization and modeling of the potential distribution of the Military Macaw (Ara militaris) in Mexico

Caracterización del hábitat y modelación de la distribución potencial de la guacamaya verde (Ara militaris) en México

Francisco A. Rivera-Ortíz¹*, Ken Oyama¹, César A. Ríos-Muñoz², Sofía Solórzano³, Adolfo G. Navarro-Sigüenza² and María Del Coro Arizmendi³

¹ Centro de Investigaciones en Ecosistemas, Universidad Nacional Autónoma de México. Antigua carretera a Pátzcuaro Núm. 8701, Colonia Ex-Hacienda de San José de La Huerta, 58190 Morelia, Michoacán, México. *frivera@cieco.unam.mx

³ Facultad de Estudios Superiores-Iztacala, Universidad Nacional Autónoma de México. Avenida de los Barrios Núm. 1, Colonia Los Reyes Iztacala, 54090 Tlalnepantla, Estado de México, México.

Recibido: 17 diciembre 2012
Aceptado: 22 julio 2013

Abstract

Forest structure and composition have been used to assess the habitat characteristics that determine bird distributions. The patterns of distribution have been shaped by historical and ecological factors that play different roles at both temporal and spatial scales. The objectives of this research were to characterize the habitat of the endangered Military Macaw (Ara militaris) and evaluate the potential distribution of this species based on trends of land use changes in Mexico. We characterized the community structure and floristic composition of 8 forests that are currently used by the Military Macaw for breeding and feeding and compared the results with 6 similar forests characterized in other studies but without historical records of the presence of the Military Macaw. The Military Macaw preferred sites with high diversity of plant species dominated by trees from 4 to 15 m in height and from 5 to 90 cm in diameter at breast height. We identified 236 plant species in the 8 forests with 20 species (8.4%) used for nesting and feeding by the Military Macaw. The floristic composition is important for the presence of the Military Macaw because there were significant differences between forests with and without its presence. The potential area of distribution of the Military Macaw had decreased by 32% and the remnant areas are included in only 8 National Protected Areas. The protected areas of natural forests should be increased to preserve the sites of potential distribution and consequently the habitat of the Military Macaw in Mexico.

Key words: Ara militaris, bird conservation, ecological niche modeling, forest community structure, habitat characterization, habitat loss.

La estructura y composición del bosque se ha utilizado para evaluar las características del hábitat que determinan la distribución de las aves. Los patrones de distribución han sido moldeadas por factores históricos y ecológicos que desempeñan diferentes papeles en ambas escalas temporales y espaciales. Los objetivos de esta investigación fueron caracterizar el hábitat de la guacamaya verde en peligro de extinción (Ara militaris) y evaluar su distribución potencial sobre las tendencias de los cambios de uso del suelo en México. Se caracterizó la estructura de la comunidad y la composición florística de 8 fragmentos remanentes de bosque que actualmente utiliza la guacamaya verde para reproducirse y alimentarse, y se compararon estos resultados con los obtenidos en bosques similares en otros estudios pero sin registros de la presencia de guacamaya verde. La guacamaya verde prefiere zonas con una alta diversidad de plantas, dominadas por árboles de 4 a 15 m de altura y de 5 a 90 cm de diámetro a la altura del pecho en las 8 localidades muestreadas. Ser identificaron 236 especies de plantas en los 8 sitios de bosque de las cuales 20 (8.4%) son utilizadas para la anidación y la alimentación de la guacamaya verde. La composición florística es importante para la presencia de estas aves, ya que hubo diferencias significativas en esta composición entre los bosques con y sin su presencia. El área de distribución potencial de esta guacamaya ha disminuido en un 32% y las áreas remanentes están incluidas únicamente en 8 Áreas Naturales Protegidas. Las áreas protegidas de bosques naturales deben de incrementarse para conservar los sitios de distribución potencial y en consecuencia el hábitat de la guacamaya verde en México.

Palabras clave: Ara militaris, conservación de aves, modelación de nicho ecológico, caracterización de hábitat, estructura de la comunidad del bosque, pérdida de hábitat.

Introduction

Forest structure and composition have been used to assess the habitat characteristics that determine bird distributions (Gillespie and Walter, 2001). Forest structure, such as leaf structural diversity, canopy coverage, volume and density of plants, and species composition, has been significantly correlated with bird distribution patterns (Gillespie and Walter, 2001; Warketin et al., 2003). These patterns of distribution have been shaped by historical and ecological factors that play different roles at both temporal and spatial scales (Vuilleumier and Simberloff, 1980; Hutto et al., 1986; Cherril and McClean, 1997; Gaston and Fuller, 2009). Food availability and habitat type determine the geographic distribution of bird species (Hutto, 1985; Orians and Wittenberger, 1991; Luke and Zack, 2001) because habitat selection by birds must ensure the availability of resources for food, nesting areas, and refuge against natural predators (Márquez-Olivas et al., 2002; Canales-del-Castillo et al., 2010; Emrick et al., 2010).

Characterization of preferred habitat of bird species could facilitate the prediction of the species' ability to respond to changes over time and space (Rotenberry, 1978), and eventually, this information may serve to support conservation policies if populations become threatened (Brower et al., 1990).

Moreover, recording the distribution and changes in forest coverage at different scales and times and mapping the physical and environmental characteristics of different ecosystems, allow conservationists to anticipate changes in land use and to suggest plans for appropriate natural resource management (Wadsworth and Treweek, 1999; Ríos-Muñoz and Navarro-Sigüenza, 2009; Contreras-Medina et al., 2010), that enable animal preservation in species threatened with extinction.

The Military Macaw (Ara militaris) is a species of conservation concern in Mexico, with a current estimated population size of less than 10,000 individuals and a clearly declining trend in that number (Collar, 1997; Snyder et al., 2000; Bird Life International, 2011). This species is included in the Appendix I of the Convention on International Trade of Endangered Species of Fauna and Flora (CITES, 1998) and is considered to be globally vulnerable (Bird Life International, 2011). In Mexico, the species is considered endangered according to federal regulations (Norma Oficial Mexicana, Nom-059-Semarnat-2010). Habitat destruction and illegal trade have been recognized as the main threats for survivorship of the Military Macaw (Íñigo-Elías, 2005).

The distribution of the Military Macaw ranges from Mexico (Sonora to Oaxaca states) to South America, with a major distribution gap in Central America (Collar, 1997; Bird Life International, 2011) that is occupied by the great Green Macaw (Ara ambigua). The Military Macaw lives in dry to semi-humid, warm sub-humid (Carreón, 1997; Rivera-Ortiz et al., 2008), and temperate (Juniper and Parr, 1998; Cruz-Nieto et al., 2006; Nocedal et al., 2006) climates with summer rains. During the reproductive season, the species distribution ranges from 0 to 2 500 m in altitude (Carreón, 1997; Rivera-Ortíz et al., 2008; Contreras-González et al., 2009).

Despite its endangered status, a habitat characterization of sites for breeding, foraging and roosting of the Military Macaw has yet to be conducted. This information and a record of the distribution range of this bird species are critical in Mexico to design conservation actions. In this paper, we characterized the habitat of the Military Macaw based on composition and structural traits of forests used by the Military Macaw for feeding and reproduction. We also modeled the potential distribution of the Military Macaw based on its ecological niche traits and the geographic distribution of associated tree species that are considered suitable habitats for the species. We conducted this study to provide information to support conservation efforts for the Military Macaw in Mexico.

Materials and methods

Study areas. This study was carried out in 8 sites in Mexico. The selected sites contained some of the largest populations reported for the Military Macaw and covered most of the distribution range of the species in Mexico (Gaucín, 2000; Gómez, 2004; Rubio et al., 2007; Rivera-Ortiz et al., 2008; Jiménez-Arcos et al., 2012). Five of the sites were located on the Pacific slope: La Sierrita, Sonora; La Reserva de Nuestra Señora del Mineral, Sinaloa; El Mirador del Águila, Nayarit; El Tuito, Jalisco, and Papalutla, Guerrero. Two other areas were located on the Gulf of Mexico slope: El Cielo, Tamaulipas and Santa María de Cocos, Querétaro, and another site was located in central Mexico: Santa María Tecomavaca, Oaxaca (Table 1; Fig. 1).

Habitat structure: characterization. Sampling was conducted in the 8 sites where the Military Macaw was observed nesting, roosting, or foraging in 2010 and 2011. We recorded the tree coverage and the density of plant species, growth form, total plant density, leaf strata diversity (LSD), plant species richness (S), plant diversity (H'), and importance-values index (IVI) (Krebs, 1985; Brower et al., 1990). The sampling included all woody trees with stems with a diameter at breast height (DBH) of 10 cm or more and shrubs taller than 1.5 m, because the Military Macaw uses different layers of the canopy forests (Forshaw, 1989).

We conducted sampling efforts directed to specific zones of nesting, roosting, or feeding the Military Macaw in each site. In these specific areas, we corroborated the presence of the species and measured vegetation cover using transects covering representative areas of vegetation used by the Military Macaw (Fig. 1). In each site, 16 transects of 50 m² were divided into 4 transects of 25 m² each, and oriented to the 4 cardinal points. The plant density was obtained by placing a rod vertically on the forest floor every 1.5 m to record the total number of plants with which the rod made contact. This procedure was repeated until 16 records were obtained along the 4 transects of each plot. For each of the trees, the name, the number of contacts with the rod as well as the height, coverage, and DBH were recorded. Plant specimens were deposited at the Herbarium of FES Iztacala (IZTA) at the Universidad Nacional Autónoma de México (UNAM).

The density of trees and shrubs (total individuals / area) and species coverage (coverage= (n x major diameter x minor diameter) / 2) according to each growth form were estimated. To estimate the LSD, the heights of the contacted plants were grouped in 16 strata of 2 m in height (0-2 m stratum, 2.1-4 m stratum and so on, up to the 26.128 m stratum) and the LSD based on the Shannon-Wiener diversity index was calculated within each stratification (MacArthur and MacArthur, 1961). Using the abundance, frequency, and vegetation coverage data of each species, we calculated the importance value index (Ar + Fr + Cr= 0 to 3) for plant species and growth forms.

An analysis of similarity of structure and composition of plant species between sites with and without presence of the Military Macaw to date was conducted. Information of sites without record of the Military Macaw was obtained from available reports (Table 2).

Statistical analyses: comparison of structure and floristic composition. Data were tested for statistical normality using Shapiro-Wilk test and Levene's homogeneity of variance using SPSS software (SPSS, 2003). The data were log10 transformed when comparisons were made using parametric tests (Sokal and Rohlf, 1979). The differences in height, coverage, and DBH were compared using an Anova (Siegel and Castellan, 2003).

The leaf structural diversity among localities was compared using the chi-square test (Siegel and Castellan, 2003). We used a t-test to compare the similarities in vegetation structure between the locations with and without records of Military Macaw (SPSS, 2003). In addition, the floristic composition between the localities with and without records of the Military Macaw was compared using the Sorensen's similarity index (PAST 2.12, 2001). Distribution models, vegetation cover changes, and environmental overlap. The models were constructed using the Genetic algorithm for rule-set production (Garp) through the desktop garp interface (Scachetti-Pereira, 2001), which has proved to be a useful tool in understanding the ecological and evolutionary processes that explain the distribution of organisms (Peterson and Navarro-Siguenza, 1999; Anderson et al., 2002; Nakazawa et al., 2004). Garp works in an iterative process where there are formation of rules that are evaluated and then considered to pass, or not, to the next generation (Stockwell and Noble, 1992, Stockwell and Peters, 1999).

We generated models of potential distribution of the Military Macaw, the potential distribution of the most important tree plant species that are used by the Military Macaw for feeding or nesting and 4 scenarios of land cover changes as follows.

Potential distribution models were created using 19 environmental variables derived from weather climatic stations stored in the Worldclim Project 1.4 (Hijmans et al., 2005) and 3 topographic layers derived from the Hydro 1k project (http://edcdaac.usgs.gov/gtopo30/hydro). Although it has been stated that the inclusion of all 19 bioclimatic variables will be prone to overfitting (Peterson and Nakazawa, 2008), the use of all variables represents a conservative and a more reliable approach to estimate the potential distribution of the species (Jakob et al., 2009).

The pixel spatial resolution was set as 0.02° by 0.02° (ca. 4 km²). The biological information used for each explanatory and response species was obtained from bibliographic sources. The geographical coordinates of the historical records for Military Macaw were obtained from the "Atlas of the birds of Mexico" (Navarro-Siguenza et al., 2003), which is the largest collection of specimens contained in Mexican and foreign collections, as well as records taken directly from fieldwork. The records obtained in the field were sightings of nesting sites, roosting, and feeding.

For plants, the geographical coordinates were obtained from the Global Biodiversity Information Facility (GBIF, http://www.gbif.org) and specialized references (Gordon, 1981; Mitchell and Daly, 1991; Gale and Pennington, 2004; Pennington and Sarukhan 2005; Rzedowski et al., 2005; Espinosa et al., 2006).

The biological information of each species was divided into training (50% of the data) and validation points (the remaining 50% of the data). In total, 100 replicates for each species were generated using a limit of convergence such that the model rules could not improve in more than 1% or 1 000 iterations. The validation of the models was made using a X² and using the training and testing data to evaluate the predictive capacities of the replicates (Peterson and Shaw 2003; Rios-Munoz and Navarro-Siguenza, 2009). Then, we selected the 10 best models as suggested by Anderson et al. (2002) based on minor errors of omission and values of commission close to the median. These distribution models were summed to obtain only one final consensus model. The consensus model was used to establish a threshold of presence/absence for all the species, which accounted for at least 90% of the biological data (Ríos-Muñoz and Navarro-Sigüenza, 2009) . Once the threshold was established, we created binary models to represent the presence/absence of the species.

The resulting distribution models were linked to watersheds and biogeographical provinces where the Military Macaw and plant species have been recorded. Hence, based upon the biogeographic history of each taxon, potential distribution models were obtained (Illoldi-Rangel and Escalante, 2008; Ríos-Muñoz and Navarro-Sigüenza, 2009). The resulting distribution scenarios were overlapped to depict coincident areas, and the final model considered as the potential distribution of the Military Macaw was the intersection with the plant species potential distributions.

To evaluate the change in vegetation cover, we used 4 databases of land use change in the country, the first for 1973-1976 was based on serial photographs (Peterson et al., 2006, Ríos-Muñoz and Navarro-Sigüenza, 2009), and the series II (1990's), III (2005), and IV (2010) of land use change in Mexico (Inegi, 2000, 2005, 2010). The use of the database linkage permitted an identification of the Military Macaw habitat loss assessment (Ríos-Muñoz and Navarro-Sigüenza, 2009). In this sense, the 4 databases were reclassified and were considered unsuitable zones for the distribution of the plant species and the macaw. These zones were urban, agricultural, forestry, livestock, grassland, and non-vegetated areas (Sánchez-Cordero et al., 2005; Peterson et al., 2006; Contreras-Medina et al., 2010) . Finally, we calculated the percentage of the area occupied in each temporary stage to obtain the pattern of habitat loss for the area associated with tree species richness. For this calculation, we considered the areas remaining in the year 2011 that were contained in the Natural Protected Areas system (NPAs) and the Important Bird Areas (IBAs) for Mexico.

Data from the same climatic layers used in the creation of the models of potential distribution for each recording site were extracted to construct a matrix to determine the environmental overlap between the Military Macaw and the most important associated tree species. With this matrix, the variation of 19 environmental parameters and altitude were subjected to Principal Component Analysis (Pca) (Novak et al., 2010; Janzekovic and Novak, 2012) and to a discriminant analysis to identify if there was a separation between species (Military Macaw and tree species). Ellipses at 95% confidence for each species were estimated (Novak et al., 2010; Janzekovic and Novak, 2012). The overlap between the areas of the ellipses was calculated using the Jacquard index ranging from 0-1 (Real and Vargas, 1996). The statistics were performed in R 3. 0. 1. (R Development Core Team, 2008).

Results

Habitat characterization. A total of 236 plant species were recorded in the 8 sampled sites. We quantified a total of 1 353 trees and 424 shrubs in the 8 sampling sites. The sites with the highest plant density were Salazares (297 ind/ha) and El Tuito (291 ind/ha). In contrast, the sites with the lowest densities were La Sierrita (177 ind/ha) and Papalutla (121 ind/ha). However, no significant differences among sites were detected (F_7,56= 0.95, p> 0.05). The tree growth form prevailed in all the sites (Table 3).

The vertical forest structure composed of 16 strata showed significant differences among the 8 sampled sites (X²= 36.124, D. F.= 15, p< 0.001); the height strata varied from 0 to 28 m across the sites (Table 3). Trees and shrubs ranging from 2 to 10 m in height dominated the vertical forest stratification in the 8 sites; however, in the localities of El Mirador del Aguila and El Cielo, the tallest trees reached over 26 m (Fig. 2).

The highest species richness was documented in El Tuito (63 species), followed by Papalutla (59 species) and Nuestra Señora del Mineral (46 species); the site with the lowest species richness was Santa María de Cocos (22 species) (Table 3). The sites with the highest plant diversity were Papalutla (H' = 3.8) and El Tuito (H' = 3.5), while Santa María de Cocos had the lowest diversity (H'= 2.2). The analysis of permutational variance indicated that the diversity of plant communities was not significantly different among the sampling sites (F_7,40= 0.83, p> 0.05) (Table 3).

The tree coverage significantly differed among the sites (F _7,38= 0.56, p< 0.001) (Table 3). The sites with the greatest tree coverage were Mirador del Águila (162.85 m²) and El Cielo (118.28 m²). In contrast, the lowest tree coverage was documented in Papalutla (39.26 m²) and Santa María de Cocos (50.73 m²). The tree growth form had the highest coverage values in all of the sampling sites. The areas with plant species with greater height and larger DBH were Mirador del Águila, El Cielo, and El Tuito (Table 3). We found significant differences in height (F_7,38= 20.17, p< 0.001) and DBH (F_7,38= 5.63, p< 0.001) among the sites.

The IVI values showed that plants sampled in all of the sites were highly variable (Appendix 1, supporting information). A total of 14 tree species (Brosimum alicastrum, Bursera simaruba, Ceiba aescutifolia, Ceiba pentandra, Cyrtocarpa procera, Guaiacum coulteri, Guazuma ulmifolia, Hura polyandra, Haematoxylon brassileto, Ipomea arborences, Lysiloma divaricata, Lysiloma microphylla, Plumeria rubra, and Taxodium mucronatum) had an IVI above 0.20 and were used for modeling their distribution in association with the modeling of the Military Macaw (see the corresponding section below).

The plants that showed the highest values of IVI were Lysiloma divaricata, L. microphyla, Brosimum alicastrum, Hura polyandra, and Cyrtocarpa procera. Important plant species that were present in more than one site were: L. divaricata, B. alicastrum, H. polyandra, Taxodium mucronatum, Bursera simaruba, and Guazuma ulmifolia. With the results of structure and composition obtained, it was observed that plant species with highest IVIs are those that the Military Macaw uses for feeding and nesting (Appendix 1).

Comparison of structure and floristic composition. In the vegetation structure, no significant differences were found comparing the cover height and DBH between sites with and without the Military Macaw [coverage (t ₍₁₁₎= 0.987, p> 0.05), height (t ₍₁₂₎= 1780, p> 0.05) and DBH (t ₍₁₂₎= 15, p> 0.05)], indicating that the forests were structurally similar. A comparison of Sorensen's similarity index among the sites confirmed 2 clearly separated groups; one contained those sites with records of the presence of the Military Macaw, and a second group formed by sites without the bird (Fig. 3).

Distribution models, land cover changes, and environmental overlap. All the models obtained presented predictions above the expected by random (X² test, all models: p< 0.01, D. F.= 1), Also, the potential distribution of the Military Macaw showed low levels of omission (i.e., the model was successful in predicting most of the primary source data), indicating a predictive power above 90%. Figure 4 shows the modeled distribution potential map of the Military Macaw and the most important plant species for feeding and nesting.

Analyzing the changes in land use from those observed in the original map (without land use) in the 1976 scenario indicated that areas with 2 to 6 species have been the most affected by the change in land use, with a reduction of 32% to 48% of their original distribution. In the Series III and Series IV, it is shown that the areas with 4 and 6 species have had a decrease of 2% and 3% respectively with respect the Series II, showing a decrease of 50%-51% of the potential distribution in comparison with the original distribution. This finding is in contrast to other areas that had 7 and 14 plant species, which were not significantly affected by land use changes, with only 10% of the original distribution reduced under the 4 scenarios (1976, 2000, 2005, 2010) (Figs. 4, 5).

The potential distribution of Military Macaw in Mexico suggests the existence of 226 000 km² of suitable climatic area without considering any impact caused by changing land use. When changes were considered, the estimated remaining area was 182 000 km², a 21.12% reduction of the original area in the 1976 scenario. For Series II, the estimation was 160 000 km² (28.82% reduction) and for Series III, the estimated remaining habitat was 158 000 km² (30.23% reduction), similar to Series IV with an estimated potential distribution of 154 000 km² (32%) (Fig. 4). This pattern showed a drastic decrease in the percentage of forest cover reaching up to 32% for the species. In 2011, the calculation of protected areas available for the Military Macaw in NPAs and IBAs only accounted for 5% and 15%, respectively, of 100% (154 000 km²) of the area distributed in 26 NPAs and 43 IBAs, along the Sierra Madre Occidental and 5 NPAs and 19 IBAs in the Sierra Madre Oriental. In the western zone, the potential area was in Sinaloa, Durango, and Guerrero, and did not include any NPAs. In the eastern zone, the NPAs and IBAs were well represented through the potential distribution of the Military Macaw.

In the PCA, component 1 explained 37.8% and component 2 explained 25.8% of the total variance of 19 environmental variables and altitude; PCA showed a single group (Fig. 6). In the discriminant analysis, the component LD1 explained 53.0 % and LD2 component explained 18.0 %; there was a clear overlap of environmental requirements of the Military Macaw with the 14 most important tree species associated to its distribution (Fig. 6). The projections of the environmental dimensions of Military Macaw and the 14 tree species are represented by ellipses in Figure 6. According to the Jacquard index, tree species distributions that showed the highest overlap with those of the Military Macaw were: Lysiloma microphylla (0.64), Lysiloma divaricata (0.53), Guaiacum coulteri (0.50), Ipomea arborences (0.50), Hura polyandra (0.46), Plumeria rubra (0.45), Guazuma ulmifolia (0.39), Haematoxylon brassileto (0.37), and Ceiba aescutifolia (.36). Species that showed lower overlap with Military Macaw were: Cyrtocarpa procera (0.27), Taxodium mucronatum (0.26), Ceiba pentandra (0.22), Bursera simaruba (0.16), and Brosimum alicastrum (0.14).

Discussion

Habitat characterization. The structural variables of the Military Macaw habitat indicated that the type of vegetation influenced the habitat selection. The Military Macaw is considered a canopy species (Inigo-Elias, 1996; Loza, 1997; Gómez, 2004) because it requires large canopy trees of deciduous and subdeciduous forests for feeding, breeding, and nesting behavior as well as protection against predators and thermal cover (Forshaw, 1989; Collar and Juniper, 1992; Collar, 1997; Loza, 1997; Íñigo-Elías, 1999; Salazar, 2001; Peterson et al., 2004; Rivera-Ortiz et al., 2008; Contreras-González et al., 2009). This species nests in trees of at least 15 m in height and the nests are 90 cm wide. In the nesting sites of El Mirador del Aguila, El Tuito, and El Cielo, the trees had the required structural characteristics for nesting (Collar, 1997; Loza 1997). The Military Macaw has the ability to shift its nesting sites to inaccessible sites such as steep cliffs in well-preserved areas: in La Sierrita, Papalutla, Santa María Tecomavaca, and Santa María de Cocos (Carreón, 1997; Gómez, 2004; Rivera-Ortiz et al., 2008).

The suitability of habitats for the Military Macaw requires the presence of certain genera of trees, such as Brosimum, Cyrtocarpa, Celtis, Hura, Quercus, Bunchonsia, Lysiloma, and Bursera; plant species of these genera have been reported in the distribution of the Military Macaw in Mexico as important sources either for nesting or as food supply by different authors (Carreón, 1997; Loza, 1997; Gaucín, 2000 and Contreras-González et al., 2009). In populations of Colombia and Peru, species of Hura and Bursera are also reported as important trees for feeding (Flores and Sierra, 2004); these plant species contain a large amount of nutrients, such as lipids, carbohydrates, and proteins, that are important for egg laying and the development of chicks (Contreras-Gonzalez et al., 2009)

Comparing the vegetation structure and floristic composition in sites with and without presence of the Military Macaw, we found significant differences in the floristic composition but structural similarities. These findings indicate the reliance of the Military Macaw on specific floristic composition, commonly found in bird specialists (such as the Military Macaw). This pattern is due to the close relationship between the availability of food resources and reproductive effort (Saunders, 1977; Saunders, 1990; Collar and Juniper, 1992) with significant implications for the conservation of this species (Ruth et al., 2003).

The information available to establish conservation strategies for the Military Macaw has been based mainly on the effects of illegal traffic and other biological and ecological aspects, such as abundance, demography, and reproduction (Carreon, 1997; Loza, 1997; Gaucin, 2000; inigo-Elias, 2000; Rivera-Ortiz et al., 2008; Contreras-González et al., 2009). Feeding and reproductive habitat modification has not taken into account in the analysis of land-use changes (Jetz and Rahbek, 2002).

Distribution models, vegetation cover changes, and environmental overlap. Ecological niche modeling represents a conceptualization of the distribution of favorable environmental conditions in which a species could be found (Peterson, 2001). Our models indicate that Military Macaw and 14 arboreal plant species are found in areas with similar characteristics, at least in a coarse environmental space; this is reinforced by the high overlap environmental found in the discriminant analysis. The reciprocal prediction of environmentally based overlap could indicate few ecological differences between the Military Macaw and tree species.

It is important to identify and preserve the habitats of endangered species with particular requirements such as the Military Macaw. We estimated a reduction of 32% in the potential distribution of the Military Macaw comparing 4 land-use change scenarios since 1976 to 2010. These changes were particularly dramatic when only 6 of the plant species that the Military Macaw relies on were present (Lysiloma microphylla, Lysiloma divaricata, Hura polyandra, Ceiba aescutifolia, Guaiacum coulteri, and Ipomea arborences). These findings indicated the potential negative impacts on the survival of the Military Macaw if reductions of available habitats occur as land-cover changes continue in the future (Peterson et al., 2006; Rios-Munoz and Navarro-Siguenza, 2009; Contreras-Medina et al., 2010). This is supported by previous studies. Rios-Munoz and Navarro-Siguenza (2009) reported a reduction of 28.5% in the available habitat of the Military Macaw by the year 2000. Marin-Togo et al. (2011) and Monterrubio-Rico et al. (2010) declared the Military Macaw locally extinct in the Mexican Pacific Coast (i.e., Michoacan, Guerrero, and Oaxaca states) and in coastal areas of more than 400 m in altitude, with a decrease of 16% of the distribution as of 2000.

The land-cover change in tropical rain forests has caused the highest rates of deforestation in the country (Trejo and Dirzo, 2000), and as a consequence, Mexican parrots have suffered severe habitat declines. Specifically, a drastic decrease has been reported in habitat occupied by Ara macao (Scarlet Macaw) (86% reduction), Aratinga astec (Aztec Parakeet) (48%), and Pionus senilis (White-crowned Parrot) (49%) (Rios-Munoz and Navarro-Siguenza et al., 2009; Marin-Togo et al., 2011). Renton and Salinas-Melgoza (2004) found that fragmentation and climatic variations of habitats in seasonally dry forests could adversely affect the reproductive success ofAmazona finschi (Liliac-crowned Parrot).

According to our results, the habitat of the Military Macaw in tropical dry forests has already been reduced drastically by almost 32%, endangering the viability of its populations. In addition, the illegal international trade of wild species has also seriously affected populations of the Military Macaw and this directly affects the loss of species distribution (Gaucín, 2000; Marín-Togo et al., 2011). Although models based on the intended habitat are very important to detect changes in the potential distribution of the Military Macaw in different scenarios, we must take into account the use of updated cartographic information of land-cover change and factors such as hunting and illegal capture to make better predictions for this species (MarinTogo et al., 2011; Monterrubio-Rico et al., 2011). Conservation implications. The present study provides information regarding the type of vegetation and species composition that is critical for the preservation of the Military Macaw. Our findings suggest the importance of knowing the floristic composition of the habitat of endangered species and the impact of land-use variation over time on the potential distribution of those species as a tool to direct conservation efforts. It is worth noting that the use of ecological niche models and geographic data of land-use change are fundamental tools to be considered in the conservation efforts of the Military Macaw. Therefore, the protection of suitable habitats and the implementation of sustainable activities should be prioritized in conservation strategies for the Military Macaw. Habitat degradation and capture of the Military Macaw for illegal trade must be stopped and the size and number of natural protected areas must be increased.

Acknowledgments

F. A. Rivera-Ortíz and C. A. Ríos-Muñoz acknowledge Conacyt for doctoral scholarships for their studies in the Posgrado en Ciencias Biológicas of the Universidad Nacional Autónoma de México. Financial support was provided by the projects of Conacyt 60270 (S. Solórzano) and DT006 (M. C. Arizmendi). Logistical support was also provided by the project SDEI-PTID-02-UNAM to P. Dávila. To local authorities for the facilities to carry out the fieldwork and plant sampling. E. Vega, L. Vázquez and R. Aguilar helped with the statistical analyses. O. Téllez and I. Calzada assisted with the plant identification. Many people participated in the field trips, but we are particularly grateful for the enthusiastic support of V.García, A. M. Contreras-González, Y. Rubio, B. Jáuregui, H. Verdugo, and E. Berrones.

Literature cited

Anderson, M. J. 2001. A new method for non-parametric multivariate analysis of variance. Austral Ecology 26:32-46.         [ Links ]

Anderson, R. P., M. Gómez-Laverde and A. Peterson. 2002. Geographical distributions of spiny pocket mice in South America: insights from predictive models. Global Ecology and Biogeography 11:131-141.         [ Links ]

Anderson, M. J. 2005. Permanova: A Fortran computer program for permutational multivarate analysis of variance. Department Statistical. University Auckland. New Zealand. 24p.         [ Links ]

Arizmendi, M. C. and L. Márquez-Valdelamar. 2000. Áreas de importancia para la conservación de las aves en México. CIPAMEX, México, D. F. 280 p.         [ Links ]

Basáñez, A. J., J. L. Alanís and E. Badillo. 2008. Composición florística y estructura arbórea de la selva mediana subperennifolia del ejido "El Remolino", Papantla Veracruz. Avances en Investigación Agropecuaria 12:3-21.         [ Links ]

Bird Life International. 2011. BirdLife's online world bird database: the site for bird conservation. Version 2.0. http://www.birdlife.org; last access: 28.III.2011.         [ Links ]

Bonadie, W. A. and P. R. Bacon. 2000. Year-round utilization of fragmented palm swamp forest by Red-bellied Macaws (Ara manilata) and Orange-winged Parrots (Amazona amazonica) in the Nariva Swamp (Trinidad). Biological Conservation 95:1-5.         [ Links ]

Brower, J., J. Zar and E. C. Von. 1990. Field and laboratory methods for general ecology. Brown Publishers. Dubuque. 220 p.         [ Links ]

Brightsmith, D. J. 2005. Competition, predation and nest niche shifts among tropical cavity nesters: phylogeny and natural history evolution of parrots (Psittaciformes) and trogons (Trogoniformes). Journal of Avian Biology 36:64-73.         [ Links ]

Canales-del Castillo, R., J. L. González-Rojas, I. Ruvalcaba-Ortega and A. García-Ramírez. 2010. New breeding localities of Worthen's Sparrows in northeastern Mexico. Journal of Field Ornithology 81:5-12.         [ Links ]

Carreón, A. G. 1997. Estimación poblacional, biología reproductiva y ecología de la nidificación de la guacamaya verde (Ara militaris) en una selva estacional del oeste del estado de Jalisco. Thesis, Facultad de Ciencias, Universidad Nacional Autónoma de México. México, D. F. 67 p.         [ Links ]

Cherril, A. and C. McClean. 1997. The impact of landscape and adjacent land cover upon linear boundary features. Landscape Ecology 12:255-260.         [ Links ]

CITES. 1998. Appendices I, II and III to the Convention on International Trade in Endangered Species of Wild Fauna and Fora. U.S. Fish and Wildlife Service, Department of the Interior, Washington D. C.         [ Links ]

Collar, N. J. and A. Juniper. 1992. Dimensions and causes of the parrot conservation crisis. In New world parrots in crisis. Solution from conservation biology, S. R. Bessingerand and N. F. R. Snyder (eds.). Smithsonian Institution Press, Washington, D. C. 1-25 p.         [ Links ]

Collar, N. J. 1997. Family Psittacidae (parrots). In Handbook of the birds of the world, 4. J. E. del Hoyo and J. Sargatal (eds.). Lynx Editions, Barcelona. 280-477 p.         [ Links ]

Conanp (Comisión Nacional de Áreas Naturales Protegidas). 2007. http://www.conanp.gob.mx; last access: 02.XII.2010.         [ Links ]

Contreras-González, A. N., F. A. Rivera-Ortíz, C. A. Soberanes-González, A. Valiente-Banuet and C. Arizmendi. 2009. Feeding ecology of Military Macaws (Ara militaris) in a semi-arid region of central Mexico. The Wilson Journal of Ornithology 121:384-391.         [ Links ]

Contreras-Medina, R. I., I. Luna-Vega and C. A. Ríos-Muñoz. 2010. Distribución de Taxus globosa (Taxaceae) en México: modelos ecológicos de nicho, efectos del cambio del uso de suelo y conservación. Revista Chilena de Historia Natural 83:421-433.         [ Links ]

Corcuera, P. J. and E. L. Butterfield. 1999. Bird communities of dry forest and oak woodland of western Mexico. Ibis 141:240-255.         [ Links ]

Cruz-Nieto, J., G. Ortiz-Maciel, M. Cruz-Nieto, M. Bujanda-Rico and E. Enkerlin. 2006. Military Macaw nesting cliff in Otachique, Chihuahua, Mexico. PsittaScene 2:18.         [ Links ]

Emrick, V. R., S. Tweddale and M. S. Germain. 2010. Characterization of Golden cheeked warbler Dendroica chrysoparia habitat at Fort Hood, Texas, USA. Endangered Species Research 11:215-220.         [ Links ]

Enkerlin-Hoeflich, E. C. and K. M. Hogan. 1997. Red-crowned Parrot (Amazona viridigenalis). In The Academy of Natural Sciences. Philadelphia, and the American Ornithologists Union, A. Poole and F. Gill (eds.). Washington D. C. 380-392 p.         [ Links ]

Espinosa, D., J. Llorente and J. J. Morrone. 2006. Historical biogeographic patterns of the species off Bursera (Burseraceae) and their taxonomical implications. Journal of Biogeography 33:1945-1958.         [ Links ]

Flores, P. and A. Sierra. 2004. Iniciativa para la conservación de la Guacamaya verde (Ara militaris) y su hábitat en el occidente de Antoquia-Colombia. Informe parcial. Fundación Proaves.         [ Links ]

Forshaw, J. M. 1989. Parrots of the world. Third Edition. Lansdowne Press. Melbourne. 180p.         [ Links ]

Gale, S. W. and T. D. Pennington. 2004. Lysiloma (Leguminosae: Mimosoideae) in Mesoamerica. Kew Bulletin 59:453-467.         [ Links ]

Galicia, L., A. García-Romero, L. Gómez-Mendoza and M. I. Ramírez. 2007. Cambio de uso del suelo y degradación ambiental. Ciencia 58:50-59.         [ Links ]

García, R. E. 1995. Estudio agroecológico de la microrregión que corresponde a la asociación agrícola local de productores de hortalizas, frutas y legumbres del municipio de Nayarit. Thesis, Facultad de Agricultura, Universidad Autónoma de Nayarit. Tepic. 46 p.         [ Links ]

Gaston, K. J. and R. A. Fuller. 2009. The sizes of species geographic ranges. Journal of Applied Ecology 46:1-9.         [ Links ]

Gaucín, R. N. 2000. Biología de la conservación de la Guacamaya verde (Ara militaris) en el Sótano del Barro, Querétaro. Facultad de Ciencias Naturales, Universidad Autónoma de Querétaro. México, Querétaro. 42 p.         [ Links ]

Gillespie, T. W. and H. Walter. 2001. Distribution of bird species richness at a regional scale in tropical dry forest of Central America. Journal of Biogeography 28:651-662.         [ Links ]

Gómez, J. O. 2004. Ecología reproductiva y abundancia relativa de la guacamaya verde en Jocotlán, Jalisco México. Thesis, Facultad de Estudios Superiores Zaragoza, Universidad Nacional Autónoma de México. México, D. F. 60 p.         [ Links ]

Gordon, J. H. 1981. Persimmons under black walnuts. Pomona 4:220-235.         [ Links ]

Herkert, J. R. 1997. Bobolink (Dolichonyx oryzivorus) population decline in agricultural landscape in the Midwestern USA. Biological Conservation 80:107-112.         [ Links ]

Hijmans, R., J. Cameron, S. E. Parra, P. G. Jones and A. Jarvis. 2005. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25:1965-1978.         [ Links ]

Hutto, R. L. 1985. Habitat selection by non- breeding, migratory land birds. In Habitat selection in birds, L. C. Cody (ed.). Academic Press, New York. 455-476 p.         [ Links ]

Hutto, R. L., S. M. Pletschet and P. Hendricks. 1986. A fixed-radius point count method for nonbreeding and breeding season use. The Auk 103:593-602.         [ Links ]

Illoldi-Rangel, P. and T. Escalante. 2008. De los modelos de nicho ecológico a las áreas de distribución geográfica. Biogeografia 3:7-12.         [ Links ]

INEGI (Instituto Nacional de Estadística y Geografia). 2000. Conjunto Nacional de Uso de Suelo y Vegetación a escala 1:250,000. Serie II. DGG-INEGI. México.         [ Links ]

INEGI (Instituto Nacional de Estadistica y Geografia). 2005. Conjunto Nacional de Uso de Suelo y Vegetación a escala 1:250,000. Serie III. DGG-INEGI. México.         [ Links ]

INEGI (Instituto Nacional de Estadistica y Geografia). 2010. Conjunto Nacional de Uso de Suelo y Vegetación a escala 1:250,000. Serie IV. DGG-INEGI. México.         [ Links ]

Íñigo-Elias, E. 1996. Ecology and breeding biology of the Scarlet Macaw (Ara macao) in the Usumacinta drainage basin of Mexico and Guatemala. Ph.D. Dissertation, University of Florida, Gainesville. 117 p.         [ Links ]

Íñigo-Elias, E. 1999. Las guacamayas verde y escarlata en México. Biodiversitas 25:7-11.         [ Links ]

Íñigo-Elias, E. 2000. Estado de conservación de las guacamayas vede (Ara militaris) y escarlata (Ara macao) en México. Audubon Latin Americana 3:1-3 p.         [ Links ]

Íñigo-Elias, E. 2005. Species assessment of resident and migrant birds in Mexico. Project Final Report to National Fish and Wildlife Foundation. Cornell Lab of Ornithology, Ithaca NY. Unpublished Report. 185 p.         [ Links ]

IUCN (International Union for Conservation ofNature). 2011. The World Conservation Union Red List of Threatened Species. http://www.iucnredlist.org; last access: 20.XII.2010.         [ Links ]

Jackson, S. F. and K. J. Gaston. 2008. Land use change and the dependence of national priority species on protected areas. Global Change Biology 14:2132-2138.         [ Links ]

Jakob, S. S., E. Martinez-Meyer and F. R. Blattner. 2009. Phylogeographic analyses and paleodistribution modeling indicate Pleistocene in situ survival of Hordeum species (Poaceae) in southern Patagonia without genetic or spatial restriction. Molecular Biology and Evolution 26:907-923.         [ Links ]

Janzekovic, F. and T. Novak. 2012. PCA: a powerful method to analyze ecological niches. In Principal component analysis-multidisciplinary applications, P. Sanguansat (ed.). Intech editions, Croatia. 127-142 p.         [ Links ]

Jetz, W. and C. Rahbek. 2002. Geographic range size and determinants of avian species richness. Science 297:1548-1551.         [ Links ]

Jiménez-Arcos, V. H., S. A. Santa Cruz-Padilla, A. Escalona-López, M. C. Arizmendi and L. Vázquez. 2012. Ampliación de la distribución y presencia de una colonia reproductiva de la guacamaya verde (Ara militaris) en el alto Balsas de Guerrero, México. Revista Mexicana de Biodiversidad 83:864-867.         [ Links ]

Juniper, T. and M. Parr. 1998. Parrots: a guide to parrots of the world. Yale University Press, New Haven. 350 p.         [ Links ]

Krebs, C. J. 1985. Ecología. Estudio de la distribución y la abundancia. 2^a edición. Harla. México, D. F. 540 p.         [ Links ]

Lanning, D. V. and J. T. Shiflett. 1983. Nesting ecology of thick-billed Parrots. The Condor 85:66-73.         [ Links ]

López-Medellín, X. and E. Íñigo-Elías. 2009. La captura de aves silvestres en México: una tradición milenaria y las estrategias para regularla. Biodiversitas 83:11-15.         [ Links ]

Loza, S. C. 1997. Patrones de abundancia, uso de hábitat y alimentación de la guacamaya verde (Ara militaris) en la Presa Cajón de Peña, Jalisco, México. Thesis, Facultad de Ciencias, Universidad Nacional Autónoma de México. México, D. F. 60 p.         [ Links ]

Luke, G.T. and S. Zack. 2001. Spatial and temporal considerations in restoring habitat for wildlife. Restoration Ecology9:272-279.         [ Links ]

MacArthur, R. H. and J. W. MacArthur. 1961. On bird species diversity. Ecology 42:594-598.         [ Links ]

Marín-Togo, M. C., T. C. Monterrubio-Rico, K. Renton, Y. Rubio-Rocha, C. Macías-Caballero, J. M. Ortega-Rodríguez and R. Cancino-Murillo. 2012. Reduced current distribution of Psittacidae on the Mexican Pacific coast: potential impacts of habitat loss and capture for trade. Biodiversity Conservation 21:451-473.         [ Links ]

Márquez-Olivas, M., L. A. Turango-Arámbula and G. D.Mendoza-Martínez. 2002. Habitat characteristics of Mexican spoteed owl (Strix occidentalis lucida) Sierra Fría. Aguascalientes. Agrociencia 36:541-546.         [ Links ]

Marsden, S. J. and A. H. Fielding. 1999. Habitat associations of parrots on the Wallacean islands of Buru, Seram and Sumba. Journal of Biogeography 26:439-446.         [ Links ]

Mitchell, J. D. and D. C. Daly. 1991. Cyrtocarpa Kunth (Anacardiaceae) in South America. Annals of the Missouri Botanical Garden 78:184-189.         [ Links ]

Monterrubio-Rico, C. T., M. J. Labrada-Hernández, J. M. Ortega-Rodríguez, R. Cancino-Murillo and J. F. Villaseñor. 2010. Distribución potencial y actual de la guacamaya verde en Michoacán, México. Revista Mexicana de Biodiversidad 82:1311-1319.         [ Links ]

Morales-Pérez, L. 2002. Efecto de la modificación del hábitat sobre la avifauna terrestre de la Reserva de la Biosfera Chamela-Cuixmala y sus alrededores. Thesis, Facultad de Ciencias, Universidad Nacional Autónoma de México. México, D. F. 65 p.         [ Links ]

Nakazawa, Y. J., A. T. Peterson, E. Martínez-Meyer and A. G. Navarro-Sigüenza. 2004. Seasonal niches of Nearctic-Neotropical migratory birds: implications for the evolution of migration. The Auk 121:610-618.         [ Links ]

Navarro-Sigüenza, A. G., A. T. Peterson and A. Gordillo-Martínez.2003. Museums working together: the atlas of the birds of Mexico. Bulletin British Ornithologists' Club 123:207-225.         [ Links ]

Nocedal, J., D. Sierra and S. Arroyo. 2006. Es la guacamaya verde realmente un ave tropical? Nidificación y alimentación en bosques templados de pino-encino del sur de Durango, México. Proceedings ofthe IV North American Ornithological Conference, Veracruz, México. October 5th, 2006. 242 p.         [ Links ]

Novak, T., C. Thirion and F. Janzekovic. 2010. Hypogean ecophase of three hymenopteran species in Central European caves. Italian Journal of Zoology 77:469-475.         [ Links ]

Ordóñez, M. J. and O. Flores. 1995. Áreas naturales protegidas. Pronatura A. C. México, D. F. 43 p.         [ Links ]

Orians, G. H. and J. F. Wittenberger. 1991. Spatial and temporal scales in habitat selection. The American Naturalist 137:829-849.         [ Links ]

Palomera-García, C., S. Contreras-Martínez and R. Amparan, R. 2007. Jalisco. In Avifaunas estatales de México. CIPAMEX, R. Ortiz-Pulido, A. Navarro-Sigüenza, H. Gómez de Silva, O. Rojas-Soto and T. A. Peterson (eds.). Pachuca, Hidalgo México. 1-48 p.         [ Links ]

Pennington, T. D. and J. Sarukhán. 2005. Árboles tropicales de México: manual para la identificación de las principales especies. 3rd edición. UNAM, FCE, México. 523 p.         [ Links ]

Peterson, A. T. 2001. Predicting species geographic distributions based on ecological niche modeling. The Condor 103:599-605.         [ Links ]

Peterson, A. T. and A. G. Navarro-Sigüenza. Alternate species concepts as bases for determining priority conservation areas. Biological Conservation 13:427-431.         [ Links ]

Peterson, A. T. and J. Shaw. 2003. Lutzomyia vectors for cutaneous leishmaniasis in Southern Brazil: ecological niche models, predicted geographic distributions, and climate change effects. International Journal of Parasitology 33:919-931.         [ Links ]

Peterson, T. M., C. Jiménez, C. Escalona-Segura, G. Flores-Villela, J. García-López, O. Hernández-Baños, B. Ruíz, A. León-Paniagua, L. Amaro, M. Navarro-Sigüenza, G. Sánchez-Cordero and D. Willard. 2004. A preliminary biological survey of Cerro Piedra Larga, Oaxaca, México: birds, mammals, reptiles, amphibians and plants. Anales del Instituto de Biología, Universidad Nacional Autónoma de México. Serie Zoología 75:439-466.         [ Links ]

Peterson, A. T., V. Sánchez-Cordero, E. Martínez-Meyer and A. G. Navarro-Sigüenza. 2006. Tracking population extirpations via melding ecological niche modeling with land-cover information. Ecological Modelling 95:229-236.         [ Links ]

Peterson, A. T. and Y. Nakazawa. 2008. Environmental data sets matter in ecological niche modelling: an example with Solenopsis invicta and Solenopsis richteri. Global Ecology and Biogeography 17:135-144.         [ Links ]

R Development Core Team. 2008. R: Un lenguaje y un entorno para cálculo estadístico. Fundación R estadística. Viena, Austria. (http://www.r-project.org.         [ Links ]).

Real, R. and J. M. Vargas. 1996. The probabilistic basis of Jaccard's index of similarity. Systematic Biology 45:380-385.         [ Links ]

Renton, K. 2001. Lilac-crowned Parrot diet and food resource availability: resource tracking by a parrot seed predator. The Condor 103:62-69.         [ Links ]

Renton, K. 2006. Diet of adult and nestling Scarlet Macaws in southwest Belize, Central America. Biotropica 38:280-283.         [ Links ]

Renton, K. and A. Salinas-Melgoza. 1999. Nesting behavior of the Lilac-crowned Parrot. Wilson Bulletin 111:488-493.         [ Links ]

Renton, K. and A. Salinas-Melgoza. 2004. Climatic variability, nest predation and reproductive output of Lilac-crowned Parrots (Amazona finschi) in tropical dry forest of western Mexico. The Auk 121:1214-1225.         [ Links ]

Ríos-Muñoz, C. A. and A. G. Navarro-Sigüenza. 2009. Efectos del cambio de uso de suelo en la disponibilidad hipotética de hábitat para los psitácidos de México. Ornitología Neotropical 20:491-509.         [ Links ]

Rivera-Ortíz, F. A., A. M. Contreras-González, C. A. Soberanes-González, A. Valiente-Banuet and M. C. Arizmendi. 2008. Seasonal abundance and breeding chronology of the Military Macaw (Ara militaris) in a semi-arid region of central Mexico. Ornitología Neotropical 19:255-263.         [ Links ]

Rodríguez-Estrella, R., E. Mata and L. Rivera. 1992. Ecological notes on the Green Parakeet of Isla Socorro, Mexico. The Condor 94:523-525.         [ Links ]

Rotenberry, J. T. 1978. Components of avian diversity along a multifactorial climatic gradient. Ecology 59:693-699.         [ Links ]

Rubio, Y., A. Beltrán, F. Aviléz, B. Salomón y M. Ibarra. 2007. Conservación de la guacamaya verde (Ara militaris) y otros psitácidos en una reserva ecológica universitaria, Cosalá, Sinaloa, México. Mesoamericana 11:58-64.         [ Links ]

Ruth, J. M., D. R. Petit, J. R. Sauer, M. D. Samuel, F. A. Johnson, M. D. Fornwall, C. E. Korschgen and J. P. Bennett. 2003. Science for avian conservation: priorities for the new millennium. The Auk 120:204-211.         [ Links ]

Rzedowski, J., L. Medina and G. Calderón-de Rzedowski. 2005. Inventario del conocimiento taxonómico, así como de la diversidad y del endemismo regionales de las especies mexicanas de Bursera (Burseraceae). Acta Botanica Mexicana 70:85-111.         [ Links ]

Salas-Morales, S. H. 2002. Relaciones entre la heterogeneidad ambiental y la variabilidad estructural de las selvas tropicales secas de la costa de Oaxaca, México. Master's thesis (Biología Ambiental) Facultad de Ciencias, Universidad Nacional Autónoma de México. México, D. F. 101 p.         [ Links ]

Salazar, J. M. 2001. Registro de la guacamaya verde (Ara militaris) en los cañones del Río Sabino y Río Seco, Santa María Tecomavaca, Oaxaca, México, Huitzil 2:18-20.         [ Links ]

Sánchez-Cordero, V., P. M. Illoldi-Rangel, A. Linaje, S. Sarkar and A. T. Peterson. 2005. Deforestation and extant distributions of Mexican endemic mammals. Biological Conservation 126:465-473.         [ Links ]

Saunders, D. A. 1977. Food and movements of the short Villeda form of the White-tailed Black Cackatoo. Australian Wildlife Research 7:257-269.         [ Links ]

Saunders, D. A. 1990. Problems of survival in an extensively cultivated landscape: the case of the Carnaby's Cackatoo Calyptorhynchus furereus latirostris. Biological Conservation 54:277-290.         [ Links ]

Scachetti-Pereira, R. 2001. Desktop GARP; http://www.lifemapper.org/desktopgarp/; last access: 23.VII.2010.         [ Links ]

Semarnat, 2010. Norma Oficial Mexicana NOM-059-SEMARNAT-2010, Protección ambiental - Especies nativas de México de flora y fauna silvestres - Categorías de riesgo y especificaciones para su inclusión, exclusión o cambio -Lista de especies en riesgo. Diario Oficial de la Federación. 30 de diciembre de 2010, Segunda Sección, México. http://www.conabio.gob.mx/informacion/catalogo_autoridades/NOM-059-SEMARNAT-2001/NOM-059-SEMARNAT-2001.pdf; last access: 15.X.2010.         [ Links ]

Siegel, S. and A. T. Castellan. 2003. Estadística no paramétrica aplicada a ciencias de la conducta. 4^a edición. Trillas. México, D. F. 519 p.         [ Links ]

Snyder, N., P. Mcgowan, J. Gilardi and A. Grajal. 2000. Parrots, status survey and conservation action plan 2000-2004. IUCN, Gland, Switzerland. 180 p.         [ Links ]

Sokal, R. R. and F. J. Rholf. 1979. Biometría. Principios y métodos estadísticos en la investigación biológica. H. Blume Ediciones. España. 426 p.         [ Links ]

Stockwell, D. and I. R. Noble. 1992. Induction of sets of rules from animal distribution data: A robust and informative method of data analysis. Mathematics and Computers in Simulation 33:385-390.         [ Links ]

Stockwell, D. and D. Peters. 1999. The GARP Modeling system: problems and solutions to automated spatial prediction. International Journal of Geographical Information Science 13:143-158.         [ Links ]

Trejo, I, 1998. Distribución y diversidad de selvas bajas de México: relaciones con el clima y el suelo. Ph.D. Thesis. División de Posgrado. Facultad de Ciencias, Universidad Nacional Autónoma de México. México, D. F. 206p.         [ Links ]

Trejo, I. and R. Dirzo. 2000. Deforestation of seasonally dry tropical forest: a national and local analysis in Mexico. Biological Conservation 94:133-142.         [ Links ]

Vázquez R. L. 2007. Relación entre la estructura de la vegetación y la comunidad de aves en la selva baja caducifolia de Santa María Tecomavaca, Oaxaca. Thesis, Facultad de Estudios Superiores Iztacala, Universidad Nacional Autónoma de México. México, D. F. 55 p.         [ Links ]

Vuilleumier, F. and D. Simberloff. 1980. Ecology versus history as determinants of patchy and insular distributions in high Andean birds. In Evolutionary biology, M. K. Hecht and W. C. Steere (eds.). Plenum Publishing Corporation, New York. 235-379 p.         [ Links ]

Wadsworth, R. A. and K. J. Trewee. 1999. Geographical information systems for ecology: an introduction. Addison Wesley Longman, Harlow. 245 p.         [ Links ]

Warketin, I. G., J. M. Reed and S. M. Dunham. 2003. Nest site characteristics of American Robins breeding in desertriparian habitat. Wilson Bulletin 1:16-23.         [ Links ]

Yamashita, C. and Y. Machado de Barros. 1997. The Blue-throated Macaw Ara glaucogularis: characterization of its distinctive habitats in savannahs of the Beni, Bolivia. Ararajuba 5:141-150.         [ Links ]

Zamora, P., G. García, J. S. Flores and J. J. Ortiz. 2008. Estructura y composición florística de la selva mediana subcaducifolia en el sur del estado de Yucatán, México. Polibotánica 26:39-66.         [ Links ]