Next Article in Journal
Ecology of Predation and Scavenging and the Interface: A Special Issue
Next Article in Special Issue
Increasing Survival of Wild Macaw Chicks Using Foster Parents and Supplemental Feeding
Previous Article in Journal
Alternative Development Strategies of Clinostomum chabaudi (Digenea) Metacercariae in Frog Hosts (Hyperolius spp.)
Previous Article in Special Issue
Parrot Ownership and Capture in Coastal Ecuador: Developing a Trapping Pressure Index
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 13
Issue 2
10.3390/d13020094
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting the Future Distribution of Ara rubrogenys, an Endemic Endangered Bird Species of the Andes, Taking into Account Trophic Interactions

by
Alain Hambuckers
1,*,
Simon de Harenne
1,2,
Eberth Rocha Ledezma
3,
Lilian Zúñiga Zeballos
1,4 and
Louis François
2
1
Behavioural Biology Lab, Unit of Research SPHERES, University of Liège, Quai Van Beneden 22, 4020 Liège, Belgium
2
Modelling of Climate and Biogeochemical Cycles, Unit of Research SPHERES, University of Liège, Sart Tilman, 4000 Liège, Belgium
3
Centro de Biodiversidad y Genética, Facultad de Ciencas y Tecnología, Universidad Mayor de San Simón, Cochabamba, Bolivia
4
Museo de Historia Natural Alcide d’Orbigny, Cochabamba, Bolivia
*
Author to whom correspondence should be addressed.
Diversity 2021, 13(2), 94; https://doi.org/10.3390/d13020094
Submission received: 22 January 2021 / Revised: 15 February 2021 / Accepted: 18 February 2021 / Published: 21 February 2021
(This article belongs to the Special Issue Ecology and Conservation of Parrots in Their Native and Non-Native Ranges)
Browse Figures
Versions Notes

Abstract

:
Species distribution models (SDMs) are commonly used with climate only to predict animal distribution changes. This approach however neglects the evolution of other components of the niche, like food resource availability. SDMs are also commonly used with plants. This also suffers limitations, notably an inability to capture the fertilizing effect of the rising CO2 concentration strengthening resilience to water stress. Alternatively, process-based dynamic vegetation models (DVMs) respond to CO2 concentration. To test the impact of the plant modelling method to model plant resources of animals, we studied the distribution of a Bolivian macaw, assuming that, under future climate, DVMs produce more conservative results than SDMs. We modelled the bird with an SDM driven by climate. For the plant, we used SDMs or a DVM. Under future climates, the macaw SDM showed increased probabilities of presence over the area of distribution and connected range extensions. For plants, SDMs did not forecast overall response. By contrast, the DVM produced increases of productivity, occupancy and diversity, also towards higher altitudes. The results offered positive perspectives for the macaw, more optimistic with the DVM than with the SDMs, than initially assumed. Nevertheless, major common threats remain, challenging the short-term survival of the macaw.

1. Introduction

The problems of habitat destruction and climate change are the main threat to tropical mountain birds. Mountain bird species in the tropics are particularly at risk because they are isolated by hotter lowland zones which often makes them sedentary. In addition, when shifting their distribution up to higher altitudes, the new area of occupancy narrows [1]. The structure of the mountains itself also appears to be a constraining factor limiting the distribution shift with possible decline of habitat quality, for instance, the absence of suitable nesting sites or even vertical gaps between actual and potential future areas of distribution [2].
The slopes of the Andes are recognized as supporting the highest avian diversity in the world combined with high endemism rate but also more than 20% of threatened species [3]. In Bolivia, the red-fronted macaw (Ara rubrogenys Lafresnaye, 1847) is one of the 15 endemic species of this country [4]. Less than 30 years ago, A. rubrogenys was a little-known species [5]. It lives on the east Andean slopes of south-central Bolivia from 553 m up to 3094 m a.s.l. (Figure 1) and breeds between 1188 and 2696 m [6]. Its natural habitat is mainly semi-deciduous dry forest but this is most often severely degraded by pastoralism and by timber extraction into thorny scrubs with scattered trees [7]. While it was estimated that the threats were of limited extent in the early nineties [8], the status of the species worsens over the course of time with land conversion for agriculture, with poaching for illegal trade, with killing by farmers who consider them a pest and with poisoning with pesticides when they feed on the crops [6,9,10,11]. The small breeding population (only 67 to 136 pairs) in eight close areas was also pointed out as major risk which increases their extinction risk due to correlated environmental fluctuations [6]. It was found that the birds use agriculture-scrub ecotones more than the forests for foraging, probably because the forests do not offer enough resources. A. rubrogenys is now ranked as “Critically Endangered” in Bolivia [11] and on the International Union for Conservation of Nature Red list [12]. In addition, it could be particularly threatened by climate change. In its area of distribution, A. rubrogenys uses only terrains along river valleys for roosting, feeding, resting and nesting [9]. Most of the nesting sites are located in steep river cliffs but such environments are not necessarily available at higher altitudes given the magnitude of warming predictions. Climate change is supposed to particularly affect the tropical Andes and notably Bolivia [13]; while warming already averaged 0.1 °C/decade between 1939 and 1998, it accelerated to 0.33 °C/decade between 1980 and 2005. Climate change scenarios suggest warming as high as 7.5 °C by 2080 and important modifications of the precipitation regime with respect to pre-industrial times.
Species distribution models (SDMs) are based on the computation of an empirical relationship between the presence of a species (a sample of its distribution) and the actual values of the selected explaining factors. Range projections are obtained by computing a probability of presence over the study area using the relationship. The methodology has been applied to parrots several times with different objectives. For instance, comparing the projection of models fitted to historical data with actual ranges allowed to study the conservation status of Andean Pyrrhura in Colombia [14] or of Amazona in Venezuela [15]. Thanks to climate driven SDMs, the substitution of the actual climate with future climate allows projection of decreasing and shifting ranges in Amazona pretrei [16]. Climate driven SDMs also permitted computing of habitat suitability of 13 parrot species in invaded countries and to test the consequences of two successive trade bans in the US and the EU on the invasion success [17]. In some situations, however, climate factors are not sufficient explaining variables, they only improve the fit of models describing the niches in combination with other factors reflecting species requirements like habitat characteristics. Authors reached this conclusion for modelling Bonelli’s eagle nesting sites with topography, disturbance, land-use or climate variable at several geographic scales [18]. Bad results were obtained with climate variables only, but climate significantly improves the quality of the prediction offered by the other sets of variables. In light of their results, the authors suggest that snow and low winter temperature may cause physiological stress hampering breeding success; however, they underlined that with more complex models, the interpretation of the effect of each explaining climate factor could become hard to find. Another interesting approach to model bird distribution consists of including biotic interactions as well as abiotic factors. For instance, in [19] SDMs driven by climate were used to simulate several shrub species making up the habitat of a bird species and the outputs of these models were set up as input variables to model animal presence with climate and topography. The authors found that this approach outperforms climate only models, which stressed the importance of taking into account, as far as possible, the different types of niche components to produce consistent simulations. This SDM approach was applied to refine the mapping of the suitability area of Amazona tucumana in Argentina and Bolivia. Here, the niche was defined with climate, land-use and the output of another SDM forecasting the distribution of a key plant resource for nestlings, Podocarpus parlatorei, providing niche cavities and food [20].
However, while SDMs driven by climate variables are now considered as a standard method to predict plant species distribution under future climate, this approach fails to consider the effect of the increasing CO2 concentration in air on plant physiology. Indeed, it is well established that increased air CO2 concentration improves the capacity of plants to resist water stress because the plants can minimize transpiration while still satisfying their CO2 requirements [21]. Contrary to SDMs, dynamic vegetation models (DVMs) are commonly able to reproduce this effect. Although, questions remain about the acclimation (organism trait responses occurring in days to weeks) and adaptation processes of plants (response occurring through evolutionary processes) to new climates which could lower future changes [22]. Forecasting under future climates with these types of models, forcing, or not, the increasing air CO2 concentration, gives contrasting results and the projected distributions of plant species or habitats under future climate with increased CO2 concentrations appear better conserved than with SDMs (e.g., [23,24]).
The objectives of this study are to evaluate the potential impact of climate change on the distribution of A. rubrogenys. We compare the results produced by SDMs driven by climate variables for A. rubrogenys and for 17 resource plant species over the area of A. rubrogenys and the results produced by a DVM for the same plant species, under present conditions and future climates (2070–2100) under RCP2.6 and RCP8.5 forcing. We predict that for the future, approaches with SDMs should conserve less of the original distribution area of A. rubrogenys than obtained with DVM.

2. Materials and Methods

2.1. Species Occurrences

Censuses of A. rubrogenys were conducted between 2008 and 2010 by two of us (E. Rocha Ledezma and L. Zúñiga Zeballos). The birds were observed between 5 and 11 a.m., in the course of linear surveys and near known nesting, feeding and roosting sites, while prospections were also realized in fields, forests and scrubs. Observations of A. rubrogenys in its area of distribution indicated that it mainly feeds on the fruits and seeds of the following wild species [25,26]: Anadenanthera colubrina (Vell.) Brenan, Anisocapparis speciosa (Griseb.) Cornejo and Iltis, Aspidosperma quebracho-blanco Schltdl., Celtis ehrenbergiana (Klotzsch) Liebm., Cenchrus (i.e., the species occurring in Bolivia in GBIF: C. alopecuroides J. Presl., C. bambusiformis (E. Fourn.) Morrone, C. brevisetus E. Fourn., C. brownii Roem. and Schult., C. chilensis (É. Desv.) Morrone, C. ciliaris L., C. echinatus L., C. insularis Scribn., C. latifolius (Spreng.) Morrone, C. longisetus M.C. Johnst, C. myosuroides Kunth, C. nervosus (Nees) Kuntze, C. polystachios (L.) Morrone, C. purpureus (Schumach.) Morrone, C. setosus Sw., C. viridis Spreng.), Cnidoscolus tubulosus I.M. Johnst., Jatropha hieronymi Kuntze, Lithraea molleoides Engl. (only occurrences in South America), Loxopterygium grisebachii Hiern ex Griseb., Neoraimondia herzogiana (Backeb.) Buxb. and Krainz, Parasenegalia visco (Lorentz ex Griseb.) Seigler and Ebinger, Parkinsonia praecox (Ruiz and Pav.) Hawkins, Prosopis chilensis Stuntz, Prosopis kuntzei Harms ex Kuntze, Sarcomphalus mistol (Griseb.) Hauenschild, Schinopsis marginata Engl. and Selaginella sellowii Hieron. We obtained the plant species coordinates of occurrences by querying the Global Diversity Information Facility site database (GBIF) in January 2017 and in August 2020, further checking for duplicates (Table 1). A second check for duplicates was conducted after combining with climate factors for coordinates belonging to the same pixels. This had limited consequences on plant sample size but not for A. rubrogenys sample size which dropped to 63. A. rubrogenys’ area of distribution was defined with the alpha-hull polygon method [27] using the ashape function of the R package “alphahull” [28].

2.2. Climate Data

For the current climate (monthly values for temperature, difference between maximum and minimum daily temperatures, precipitation, relative humidity, sunshine hours, wind speed), we used Worldclim version 2 for the time period between 1970 and 2000 at 2.5 arc-minutes [29]. For the future climate, we used the CMIP5 projections of the HadGEM2-AO global circulation model [30] under the representative concentration pathways of greenhouse gases corresponding to end of 21st century radiative forcings of 2.6 and 8.5 W/m2 (RCP2.6 and RCP8.5, [31]). RCP2.6 would require a decline of greenhouse gas to reach no emission after 2072. With respect to the pre-industrial period, this scenario would keep global temperature rise below 2 °C in 2100. RCP8.5 is the worst hypothesis with emission of greenhouse gases continuing to increase throughout the 21st century producing global temperature increase between 2.6 and 4.8 °C for the period 2081–2100 [32].

2.3. Dynamic Vegetation Modelling

The DVM CARAIB (CARbon Assimilation In the Biosphere) was mainly described in [23,33,34,35,36]. This model was initially conceived to simulate vegetation at global or continental scale and its response to climate change in the future or in the past [37,38,39,40,41,42,43]. The model was also applied to agricultural systems [44,45,46] or to tree species [24,47,48,49]. It is a grid point model composed of 5 main interacting modules (hydrology budget, photosynthesis and stomatal regulation, carbon allocation and growth, heterotrophic respiration and carbon dynamics in the soil, competition between ecosystem strata and biogeography) and it is also possible to activate natural fire and migration modules. Input data are spatial monthly climates (minimal and maximal temperature, precipitation, relative humidity, sunshine hours and wind speed), CO2 air concentration, soil texture and color, elevation and a set of information describing the morpho-physiological characteristics of the plant species (traits), like the specific leaf area, leaf and sapwood C:N, plant height, deciduousness nature, et cetera, and climatic thresholds extracted from the distribution samples. Since trait information for the species was lacking, it was replaced by the values of the plant functional type to which the species belongs (Table 4 of Electronic Supplementary Material of [49]). Threshold values controlling germination and mortality under stress conditions are extracted from prescribed percentiles in their actual climate distribution extracted from the occurrence samples [36,40]. For computing the fitness statistics, we also ran the simulation for the present on the coordinates of the sets of pseudo-absences drawn for SDM modelling (see below). As output of the model, we used the net primary productivities of the species fractions (fNPP, gC/m2/y) and we selected thresholds of presence maximizing the true skill statistic (TSS, [50]). The threshold of presence allows us to compute the sensitivity (proportion of presences correctly predicted) and specificity (proportion of absences correctly predicted). We also computed the receiver operating characteristic curve (AUC). Positive TSS values and AUC larger than 0.7 indicate better agreement than random.

2.4. A. rubrogenys and Plant Species SDM

We used multiple logistic regressions, also called logit models. Studies have showed that it is difficult to rank the SDM methods according to their performances because they vary with dataset properties [51,52,53], but logistic regression has the best theoretical background and it is used here as a reference methodology for SDM. The models were computed in R with the glm function of the R package “stat” (The R Core Team) and we interpreted the output as a probability of presence (Px). Since absence data are also needed, we generated pseudo-absences for each species (A. rubrogenys and plant species). The pseudo-absences were sets of points randomly drawn around the occurrences containing twice as many points as occurrences. Our strategy is based on the fact that for rare events, rather than randomly sampling, it is more efficient to collect positive cases (presences), as much as possible, and to complete them with a limited number of random cases as absences [54]. There is no common rule to select the ratio of absences to presences and this varied considerably in ecological studies (see [55]). For each plant species, we built 10 datasets containing the occurrence coordinates and pseudo-absence coordinates drawn in radius buffers comprised of between 200 and 2000 km around presences. Furthermore, drawing pseudo-absences within the presences is justified by the fact that the species do not occupy their entire range. For each dataset, we identified the best combination of climate factor effects, in other words, a maximum of 6 (only 4 for N. herzogiana, due to smaller sample size) linear or quadratic effects, after exhaustive screening based on the Akaike information criterion (AIC), using the glmulti function of the R package “glmulti” [56]. This procedure allowed us to select the shortest distance for the pseudo-absences providing evident AUC increase, which may give more accuracy and meaningful fit of the models [57]. We finally used the datasets giving the highest AUC and the model selected in the above procedure with the lowest AIC. Possibly, we chose a model with slightly higher AIC so that most of the model coefficients may have Z-test p-values lower than 0.05. For validation we tested the selected models against their null models using the likelihood ratio test (with critical p-value = 0.05) and the AUC (auc function of “SDMTools” [58]). The cutting thresholds were also selected thanks to TSS. A linear effect alone reveals a strictly positive or negative ecological response to the considered factor. The combination of a linear effect with the quadratic effect of the same factor produces a bell-shaped ecological response curve. Owing to the logit link function, a quadratic effect alone acts as a threshold. As climate variables, we first selected annual mean temperature, temperature seasonality, maximum temperature of the warmest month, minimum temperature of the coldest month, annual precipitation, precipitation seasonality, precipitation of the wettest quarter and precipitation of the driest quarter. Then, to minimize collinearity, we computed the matrix of Pearson correlation coefficients in the species datasets and pointed the couple of variables with coefficients >0.7 [59]. Thus, we had to drop the maximum temperature of the warmest month, minimum temperature of the coldest month and annual precipitation for the plant species, and additionally, precipitation of the driest quarter for A. rubrogenys. For projections, we set missing value pixels with at least one of the explanatory factors outside the range encountered in the calibration datasets.
It should be noted that we tested the approach consisting of including outputs of the plant models to drive A. rubrogenys SDM using restricted estimation of the linear coefficients of the plant model output to impose positive coefficients. However, this approach failed mainly due to the fact that the productivity (DVM outputs) or the probability of presence (SDM outputs) were not uniformly higher over the A. rubrogenys area than outside. We understood this negative result as that animals are able to rely on resources made by plants growing in sub-optimal conditions, which seemed evident.

3. Results

Climate over the A. rubrogenys area varied a lot, exemplified by mean annual temperatures which fall between 12.31 and 23.00 °C and precipitation of the driest quarter which ranged between 6 and 62 mm (Table S1). Under RCP2.6 forcing, temperature and precipitation over the A. rubrogenys area increased while seasonality decreased with nevertheless increases of the ranges between extreme higher and lower values. Under RCP8.5 forcing, temperature and precipitation followed the same direction but more intensively than under RCP2.6 forcing.

3.1. Dynamic Vegetation Modelling

Fitness statistics (Table 2) suggested limited to acceptable agreements for plant species DVM modelling. It is normal that the mechanistic model had significantly lower AUC than those obtained with the SDM. Indeed, first there is no training step in the DVM computation and, second, the AUC of the SDM has been maximized in the generation of pseudo-absences by varying the distance to the points of presence. Thus, comparing the performances of the two models in terms of AUC, or any other evaluator using pseudo-absences, cannot be made without bias. In this respect, it should be noted that the performance of the DVM in terms of sensitivity (Se) (which do not use pseudo-absences) is generally quite high, since they are generally larger than 0.8 or even 0.9. Otherwise, owing to excessive computing time, it was not possible to compute maps over the entire areas of distribution of the plant species. All the species but C. tubulosus were simulated as occurring over the A. rubrogenys presence area (Figure 2 and Figure 3 and Figure S1). The range of fNPP also varied considerably between the species and the climate conditions.

3.2. SDM Modelling

It was possible to obtain models with acceptable to very good agreement for each of the plant species and for A. rubrogenys, with fitness indicators generally higher than those obtained for the DVM results (Table 3), except for the sensitivity (Se). For plants, projection maps revealed large areas of potential presence without any presence point which could be due to uneven sampling but also to natural barriers to migration, to competition or to speciation (Figure 4 and Figure 5, Figures S2 and S3). Ten of the plant species were largely simulated as occurring over the A. rubrogenys’ present area (Figure 6). The simulated distributions over Bolivia sometimes showed large similarities with those produced by the DVM. For A. rubrogenys, the model defined a close region over the alpha hull defining the present area (Figure 5). The original range was conserved under the future climates, the Px increased (Figure 6), while new suitable areas appeared connected with the former one (Figure 5).

3.3. Comparisons of Plant Model Predictions for the Future

Under future conditions, changes predicted by the DVM over Bolivia were limited but plant presences tended rather to spread (Figure 2 and Figure 3 and Figure S1). The trend was more pronounced over A. rubrogenys’ area, particularly under RCP8.5 forcings while the fNPP of the pixel presence over the A. rubrogenys area behaved more or less similarly but with some exceptions (Cenchrus, P. visco, Figure 3). The mean number of plant species per pixel (plant diversity, Figure S4) computed for present climate was 10.3. It increased to 12.0 under RCP2.6 and 14.2 under RCP8.5 forcings. Furthermore, 80% (RCP2.6) to 96% (RCP8.5) of the pixels lost no species while a majority of them gained up to 11 new species (Figure 7). In addition, the maximal altitude increased for the majority of the species under the future climates (Figure 8).
With the SDMs, changes predicted over Bolivia were more significant than with the DVM (Figure 5 and Figure 6 and Figure S3) but generally the presence also tended to spread. Over A. rubrogenys’ area, no clear trend emerged for occupancy or Px. Diversity over pixels (Figure S4) for the present was higher than with the DVM (mean: 13.2) but decreased to 12.6 under RCP2.6 and to 13 under RCP8.5 forcings. Indeed, more pixels lost species while the gains were more modest (Figure 7). Furthermore, less species than with the SDM had maximal altitude increase under the future climates (Figure 8).

4. Discussion

The SDM allowed accurately modelling of A. rubrogenys presence with climate. Under future climate, occupancy, probability of presence and maximal altitude increased and new areas of occupancy appeared, connected with its present area. For the plants, DVM and SDM modelling allowed us to successfully compute the distribution of 17 plant species which are feeding resources of A. rubrogenys. Under present conditions, only one species, C. tubulosus, was not simulated as present over the area of A. rubrogenys by the DVM. SDMs gave better fits than the DVM, except for sensitivity (Se). The DVM mainly produced fNPP, and occupancy increases over the area of A. rubrogenys, also towards higher altitudes, under the RCP2.6 and RCP8.5 forcings. With DVM, we also obtained pixel diversity increases resulting from few species loss and more species gains. By contrast, the plant SDMs had no overall response under the future climate conditions for occupancy, probability of presence, altitude or diversity.
The lowest agreement of the DVM simulations compared to the plant SDMs were the result of lower specificity of the DVM. Here, we added an optimization step in which we selected a threshold for each species thanks to TSS. However, these thresholds were low compared to commonly fixed values. They probably reflected the fact that we used fNPP and thus the competition between species for water and light. Lowering the thresholds increases sensitivity at the expense of specificity. It was already noted in previous studies that the DVM may tend to simulate the fundamental niche of the species rather than the realized niche, in the absence of some critical biotic interactions. This could ultimately produce wider distributions compared to those produced by the SDMs as observed elsewhere [49,60]. Nevertheless, the results of both approaches, DVM and SDMs, were rather congruent for present conditions. For the future, the results were mostly in accordance with literature, in other words, the SDMs generally predict more widespread distribution shifts than DVM [24,49,61,62,63] while this is not absolute. Therefore, it is not surprising that in limited parts of their present area, the Px of some species computed with SDM increased under the RCP forcings. In the Andes, SDMs predicted with dramatic reductions of the Andean vascular plant diversity for the future without excluding some small areas of better climatic stability where species numbers could increase [64]. Otherwise, the DVM results showing increases in diversity and the migrations towards higher altitudes seemed in accordance with the observations from the recent past. Permanent plots in the Andes showed thermophilization, in other words, the increasing number of tropical and subtropical tree species from lower altitudes showing shifts to higher altitudes [65]. This phenomenon, was also observed on European mountain summits. However, it should finally threaten the mountains species of the former communities because they are composed of slower growing species more adapted to harsher conditions compared to species from lower altitudes [66].
It is established that birds are directly sensitive to climate factors. While one of the main advantages of body temperature regulation in homeotherms is to allow an optimal functioning of the organisms under a wide temperature range, the limit of the mechanism is clearly the dissipation of heat in excess during exercise [67]. Thus, air temperature has to be considered as an effective component of the bird’s niche. Homeotherms could also be directly sensitive to other climate factors than temperature, for instance, evaporation can limit the survival of birds in deserts during heat waves [68]. Furthermore, it was showed that birds have following their realized climate niche during the course of the last century [69]. For these reasons, the use of species distribution models (SDMs) with climate factors only to predict their future distribution has to be considered as a useful tool. The common response of birds in the Andes to climate change would be a decline [64,70]. Species with restricted distribution, like endemic species, are generally characterized by narrow climate niche and they seem to resist extinction by relatively high local abundance and good demographic resilience resulting from the accumulation of local adaptions [71]. The study of species, including parrots, naturalized outside their original area of distribution indicates however that the climate tolerance could be higher than estimated in the natural area, particularly for the species occupying narrow ranges of climate conditions or marginal climates in their native region [72,73,74]. With the climate factor, we obtained a very well-defined geographic range (Figure 5). We supposed that the climate niche of A. rubrogenys would be broad compared to the other endemic species on the basis of the range of the climate factors over its present area provoked by the steepness of the climate gradient in the mountains (Table S1). Therefore, it would not be surprising that the SDM is driven by climate predicted occupancy increase of the present area under the RCP forcings and also on areas directly connected, at higher altitudes. However, this kind of situation would be rare. Another possibility would be that the area of A. rubrogenys is climatically stable enough to guarantee no species loss or even species gain even if those situations are rare [64]. If so, it would also explain the limited changes found for the plant species under the RCP forcings with both modelling approaches. Thus, as long as the non-climate components of the niche of A. rubrogenys remain available, we could venture to forecast limited risk of extinction due to climate change provided that the diversity of the plant resources could increase and that species turnover could satisfy animal needs.
While the DVM approach is more promising than the SDM one for evaluating plant future, DVMs are nevertheless perfectible. DVMs need to elaborate an optimal strategy for parametrization [45,48,75]. It is indeed limited by the knowledge of the species-specific traits but also by the response of those traits to environmental factors and results could be improved by determining those responses and integrating them into the DVM [76,77,78]. The collection of plant traits with biogeographic information is a fast-growing field of knowledge [79]. Since the DVM simulates competition for light and water resources, another interesting point could be to integrate the species of the former communities with the potential newcomers to compute the emerging communities. The newcomers, in addition, could also constitute potential resources.

5. Conclusions

Under the milder RCP2.6 and the harsher RCP8.5 forcings, the SDM for A. rubrogenys and SDMs and the DVM for plant feeding resources of A. rubrogenys showed that the current area of suitability would be mostly preserved or even broadened. The actual feeding resources could be more precisely evaluated by taking into account the land-use and the productivity of the agricultural species which constitute important feeding resources for A. rubrogenys [6,25,26]. Nevertheless, predictions of the evolution of those resources in the course of time would require additional effort since uncertainties increase with the characteristics of the cultivated varieties and the decisions of the landowners or the farmers. The use of the DVM for predicting the future of plant species distribution and productivity is thus improvable. Our result boded well because A. rubrogenys would not be subjected to unprecedented climate conditions in the present area and the feeding resources so far would be preserved.
From topology, it could be possible to examine whether new suitable areas include cliffs appropriate for nesting. However, is it possible that the new suitable areas could be colonized even without appropriate cliffs thanks to behavioral flexibility? For instance, A. rubrogenys have been observed nesting in the palm Parajubaea torallyi [80] but, this tree species has an extremely small range distribution, and there are no other tree species that could offer large cavities for nesting within the range of the red-fronted macaw. Moreover, recent population genetic analyses show a low capacity for the species to colonize distant areas. Despite its restricted range, the population is structured in genetic clusters with low or null gene flow among them, despite that some colonies are separated by few tens of kilometers, that there are no ecological barriers and that macaws make larger daily and seasonal movements for feeding. Therefore, it is highly improbable that the species could colonize very distant areas in its future niche suitability [81].
Thus, while bird decline at global scale seems mainly provoked by increasing temperature [82], the main concerns about A. rubrogenys’ future remained the direct human threats, in other words, habitat conversion, poaching and killing, which need to be urgently solved before the species rapidly goes extinct in nature [6,11,12].

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-2818/13/2/94/s1, Table S1: Bioclimate variables over Ara rubrogenys area and changes under the RCP2.6 and RCP8.5 forcings, Table S2: Coefficients of the climate logistic models for Ara rubrogenys and the 17 plant species, Figure S1: Plant species distribution over Bolivia predicted by the DVM CARAIB for present and future under RCP2.6 andRCP8.5 forcings, Figure S2: Plant species continental distribution predicted by the SDMs for the present, Figure S3: Plant species distribution over predicted by the SDM for present and future under RCP2.6 andRCP8.5 forcings, Figure S4: Number of plant species over the area of Ara rubrogenys computed with the DVM or the SDMs for present and future climate under RCP2.6 and RCP8.5 forcings.

Author Contributions

Conceptualization, methodology and investigation, A.H. and L.F.; investigations, E.R.L., S.d.H. and L.Z.Z.; data curation, L.Z.Z., writing—original draft preparation and editing, A.H.; review, L.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Information on data from publicly archived datasets used in this study is given in the Section 2.

Acknowledgments

We thank the reviewers for their comments and suggestions and their sharing of information on the red-fronted macaw which allowed us to improve our study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Şekercioğlu Çağan, H.; Primack, R.B.; Wormworth, J. The effects of climate change on tropical birds. Biol. Conserv. 2012, 148, 1–18. [Google Scholar] [CrossRef]
  2. La Sorte, F.A.; Jetz, W. Tracking of climatic niche boundaries under recent climate change. J. Anim. Ecol. 2012, 81, 914–925. [Google Scholar] [CrossRef]
  3. Orme, C.D.L.; Davies, R.G.; Burgess, M.; Eigenbrod, F.; Pickup, N.; Olson, V.A.; Webster, A.J.; Ding, T.-S.; Rasmussen, P.C.; Ridgely, R.S.; et al. Global hotspots of species richness are not congruent with endemism or threat. Nat. Cell Biol. 2005, 436, 1016–1019. [Google Scholar] [CrossRef]
  4. Herzog, S.K.; Maillard, O.; Embert, D.; Caballero, P.; Quiroga, D. Range size estimates of Bolivian endemic bird species revisited: The importance of environmental data and national expert knowledge. J. Ornithol. 2012, 153, 1189–1202. [Google Scholar] [CrossRef]
  5. Lanning, D.V. Distribution and breeding biology of the Red-fronted Macaw. Wilson Bull. 1991, 103, 357–365. [Google Scholar]
  6. Tella, J.; Rojas, A.; Carrete, M.; Hiraldo, F. Simple assessments of age and spatial population structure can aid conservation of poorly known species. Biol. Conserv. 2013, 167, 425–434. [Google Scholar] [CrossRef]
  7. Herzog, S.K.; Kessler, M.; Maijer, S.; Hohnwald, S. Distributional notes on birds of Andean dry forests in Bolivia. Bull. Br. Ornithol. Club 1997, 117, 223–235. [Google Scholar]
  8. Boussekey, M.; Saint-Pie, J.; Morvan, O. Observations on a population of Red-fronted Macaws Ara rubrogenys in the Río Caine valley, central Bolivia. Bird Conserv. Int. 1991, 1, 335–350. [Google Scholar] [CrossRef] [Green Version]
  9. Pitter, E.; Christiansen, M.B. Ecology, status and conservation of the Red-fronted Macaw Ara rubrogenys. Bird Conserv. Int. 1995, 5, 61–78. [Google Scholar] [CrossRef] [Green Version]
  10. Herrera, M.; Hennessey, B. Quantifying the illegal parrot trade in Santa Cruz de la Sierra, Bolivia, with emphasis on threatened species. Bird Conserv. Int. 2007, 17, 295–300. [Google Scholar] [CrossRef] [Green Version]
  11. Rojas, A.; Zeballos, A.; Rocha, E.; Balderrama, J.A. Ara rubrogenys Lafrenaye 1847; Ministerio de Medio Ambiente y Agua: La Paz, Bolivia, 2009; pp. 332–334.
  12. BirdLife International. Ara rubrogenys. The IUCN Red List of Threatened Species. 2018. Available online: https://dx.doi.org/10.2305/IUCN.UK.2018-2.RLTS.T22685572A131382876.en (accessed on 21 September 2020).
  13. Rangecroft, S.; Harrison, S.; Anderson, K.; Magrath, J.; Castel, A.P.; Pacheco, P. Climate Change and Water Resources in Arid Mountains: An Example from the Bolivian Andes. Ambio 2013, 42, 852–863. [Google Scholar] [CrossRef] [Green Version]
  14. Botero-Delgadillo, E.; Páez, C.A.; Bayly, N. Biogeography and conservation of Andean and Trans-Andean populations of Pyrrhura parakeets in Colombia: Modelling geographic distributions to identify independent conservation units. Bird Conserv. Int. 2012, 22, 445–461. [Google Scholar] [CrossRef] [Green Version]
  15. Ferrer-Paris, J.R.; Sánchez-Mercado, A.; Rodríguez-Clark, K.M.; Rodríguez, J.P.; Rodríguez, G.A. Using limited data to detect changes in species distributions: Insights from Amazon parrots in Venezuela. Biol. Conserv. 2014, 173, 133–143. [Google Scholar] [CrossRef]
  16. Marini, M.Â.; Barbet-Massin, M.; Martinez, J.; Prestes, N.P.; Jiguet, F. Applying ecological niche modelling to plan conservation actions for the Red-spectacled Amazon (Amazona pretrei). Biol. Conserv. 2010, 143, 102–112. [Google Scholar] [CrossRef]
  17. Cardador, L.; Lattuada, M.; Strubbe, D.; Tella, J.L.; Reino, L.; Figueira, R.; Carrete, M. Regional Bans on Wild-Bird Trade Modify Invasion Risks at a Global Scale. Conserv. Lett. 2017, 10, 717–725. [Google Scholar] [CrossRef] [Green Version]
  18. López-López, P.; García-Ripollés, C.; Aguilar, J.M.; Garcia-López, F.; Verdejo, J. Modelling breeding habitat preferences of Bonelli’s eagle (Hieraaetus fasciatus) in relation to topography, disturbance, climate and land use at different spatial scales. J. Ornithol. 2005, 147, 97–106. [Google Scholar] [CrossRef]
  19. Preston, K.L.; Rotenberry, J.T.; Redak, R.A.; Allen, M.F. Habitat shifts of endangered species under altered climate conditions: Importance of biotic interactions. Glob. Chang. Biol. 2008, 14, 2501–2515. [Google Scholar] [CrossRef]
  20. Pidgeon, A.M.; Rivera, L.; Martinuzzi, S.; Politi, N.; Bateman, B.L. Will representation targets based on area protect critical resources for the conservation of the Tucuman Parrot? Condor 2015, 117, 503–517. [Google Scholar] [CrossRef] [Green Version]
  21. Fitter, A.; Hay, R. Environmental Physiology of Plants; Elsevier: Amsterdam, The Netherlands, 2002; p. 367. [Google Scholar]
  22. Smith, N.G.; Dukes, J.S. Plant respiration and photosynthesis in global-scale models: Incorporating acclimation to temperature and CO. Glob. Chang. Biol. 2012, 19, 45–63. [Google Scholar] [CrossRef]
  23. Dury, M.; Hambuckers, A.; Warnant, P.; Henrot, A.; Favre, E.; Ouberdous, M.; François, L. Responses of European forest ecosystems to 21st century climate: Assessing changes in interannual variability and fire intensity. iForest Biogeosci. For. 2011, 4, 82–99. [Google Scholar] [CrossRef]
  24. Raghunathan, N.; François, L.; Huynen, M.-C.; Oliveira, L.C.; Hambuckers, A. Modelling the distribution of key tree species used by lion tamarins in the Brazilian Atlantic forest under a scenario of future climate change. Reg. Environ. Chang. 2014, 15, 683–693. [Google Scholar] [CrossRef]
  25. Blanco, G.; Hiraldo, F.; Rojas, A.; Dénes, F.V.; Tella, J.L. Parrots as key multilinkers in ecosystem structure and functioning. Ecol. Evol. 2015, 5, 4141–4160. [Google Scholar] [CrossRef] [Green Version]
  26. Montesinos-Navarro, A.; Hiraldo, F.; Tella, J.L.; Blanco, G. Network structure embracing mutualism–antagonism continuums increases community robustness. Nat. Ecol. Evol. 2017, 1, 1661–1669. [Google Scholar] [CrossRef]
  27. Burgman, M.A.; Fox, J.C. Bias in species range estimates from minimum convex polygons: Implications for conservation and options for improved planning. Anim. Conserv. 2003, 6, 19–28. [Google Scholar] [CrossRef] [Green Version]
  28. Pateiro-López, B.; Rodríguez-Casal, A. Generalizing the Convex Hull of a Sample: TheRPackagealphahull. J. Stat. Softw. 2010, 34, 1–28. [Google Scholar] [CrossRef] [Green Version]
  29. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  30. The HadGEM2 Development Team; Martin, G.M.; Bellouin, N.; Collins, W.J.; Culverwell, I.D.; Halloran, P.R.; Hardiman, S.C.; Hinton, T.J.; Jones, C.D.; McDonald, R.E.; et al. The HadGEM2 family of Met Office Unified Model climate configurations. Geosci. Model Dev. 2011, 4, 723–757. [Google Scholar] [CrossRef] [Green Version]
  31. Moss, R.H.; Edmonds, J.A.; Hibbard, K.A.; Manning, M.R.; Rose, S.K.; Van Vuuren, D.P.; Carter, T.R.; Emori, S.; Kainuma, M.; Kram, T.; et al. The next generation of scenarios for climate change research and assessment. Nature 2010, 463, 747–756. [Google Scholar] [CrossRef] [PubMed]
  32. Stocker, T.F.; Qin, D.; Plattner, G.-K.; Tignor, M.; Allen, S.K.; Boschung, J.; Nauels, A.; Xia, Y.; Bex, V.; Midgley, P.M. Climate Change 2013: The Physical Science Basis; Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013. [Google Scholar] [CrossRef] [Green Version]
  33. Warnant, P.; François, L.; Strivay, D.; Gérard, J.-C. CARAIB: A global model of terrestrial biological productivity. Glob. Biogeochem. Cycles 1994, 8, 255–270. [Google Scholar] [CrossRef]
  34. Gérard, J.C.; Nemry, B.; François, L.M.; Warnant, P. The interannual change of atmospheric CO2: Contribution of subtropical ecosystems? Geophys. Res. Lett. 1999, 26, 243–246. [Google Scholar] [CrossRef]
  35. Otto, D.; Rasse, D.; Kaplan, J.; Warnant, P.; François, L. Biospheric carbon stocks reconstructed at the Last Glacial Maximum: Comparison between general circulation models using prescribed and computed sea surface temperatures. Glob. Planet. Chang. 2002, 33, 117–138. [Google Scholar] [CrossRef]
  36. Laurent, J.; Francois, L.; Barhen, A.; Bel, L.; Cheddadi, R. European bioclimatic affinity groups: Data-model comparisons. Glob. Planet. Chang. 2008, 61, 28–40. [Google Scholar] [CrossRef]
  37. François, L.; Goddéris, Y.; Munhoven, G. Modelling the interactions between biospheric and weathering processes: Towards a mechanistic description of the land environment. Mineral. Mag. 1998, 62A, 468–469. [Google Scholar] [CrossRef]
  38. François, L.M.; Goddéris, Y.; Warnant, P.; Ramstein, G.; De Noblet, N.; Lorenz, S. Carbon stocks and isotopic budgets of the terrestrial biosphere at mid-Holocene and last glacial maximum times. Chem. Geol. 1999, 159, 163–189. [Google Scholar] [CrossRef]
  39. François, L.; Ghislain, M.; Otto, D.; Micheels, A. Late Miocene vegetation reconstruction with the CARAIB model. Palaeogeogr. Palaeoclim. Palaeoecol. 2006, 238, 302–320. [Google Scholar] [CrossRef] [Green Version]
  40. François, L.; Utescher, T.; Favre, E.; Henrot, A.-J.; Warnant, P.; Micheels, A.; Erdei, B.; Suc, J.-P.; Cheddadi, R.; Mosbrugger, V. Modelling Late Miocene vegetation in Europe: Results of the CARAIB model and comparison with palaeovegetation data. Palaeogeogr. Palaeoclim. Palaeoecol. 2011, 304, 359–378. [Google Scholar] [CrossRef]
  41. Henrot, A.-J.; François, L.; Favre, E.; Butzin, M.; Ouberdous, M.; Munhoven, G. Effects of CO2, continental distribution, topography and vegetation changes on the climate at the Middle Miocene: A model study. Clim. Past 2010, 6, 675–694. [Google Scholar] [CrossRef] [Green Version]
  42. Henrot, A.-J.; Utescher, T.; Erdei, B.; Dury, M.; Hamon, N.; Ramstein, G.; Krapp, M.; Herold, N.; Goldner, A.; Favre, E.; et al. Middle Miocene climate and vegetation models and their validation with proxy data. Palaeogeogr. Palaeoclim. Palaeoecol. 2017, 467, 95–119. [Google Scholar] [CrossRef]
  43. Chang, J.; Ciais, P.; Wang, X.; Piao, S.; Asrar, G.; Betts, R.; Chevallier, F.; Dury, M.; François, L.; Frieler, K.; et al. Benchmarking carbon fluxes of the ISIMIP2a biome models. Environ. Res. Lett. 2017, 12, 045002. [Google Scholar] [CrossRef] [Green Version]
  44. Fontaine, C.M.; Dendoncker, N.; De Vreese, R.; Jacquemin, I.; Marek, A.; Van Herzeleg, A.; Devillet, G.; Mortelmans, D.; François, L. Towards participatory integrated valuation and modelling of ecosystem services under land-use change. J. Land Use Sci. 2014, 9, 278–303. [Google Scholar] [CrossRef]
  45. Minet, J.; Laloy, E.; Tychon, B.; Francois, L. Bayesian inversions of a dynamic vegetation model at four European grassland sites. Biogeosciences 2015, 12, 2809–2829. [Google Scholar] [CrossRef] [Green Version]
  46. Fronzek, S.; Pirttioja, N.; Carter, T.R.; Bindi, M.; Hoffmann, H.; Palosuo, T.; Ruiz-Ramos, M.; Tao, F.; Trnka, M.; Acutis, M.; et al. Classifying multi-model wheat yield impact response surfaces showing sensitivity to temperature and precipitation change. Agric. Syst. 2018, 159, 209–224. [Google Scholar] [CrossRef]
  47. Cheddadi, R.; Henrot, A.-J.; François, L.; Boyer, F.; Bush, M.; Carré, M.; Coissac, E.; De Oliveira, P.E.; Ficetola, F.; Hambuckers, A.; et al. Microrefugia, Climate Change, and Conservation of Cedrus atlantica in the Rif Mountains, Morocco. Front. Ecol. Evol. 2017, 5. [Google Scholar] [CrossRef] [Green Version]
  48. Dury, M.; Mertens, L.; Fayolle, A.; Verbeeck, H.; Hambuckers, A.; François, L. Refining Species Traits in a Dynamic Vegetation Model to Project the Impacts of Climate Change on Tropical Trees in Central Africa. Forests 2018, 9, 722. [Google Scholar] [CrossRef] [Green Version]
  49. Raghunathan, N.; François, L.; Dury, M.; Hambuckers, A. Contrasting climate risks predicted by dynamic vegetation and ecological niche-based models applied to tree species in the Brazilian Atlantic Forest. Reg. Environ. Chang. 2018, 19, 219–232. [Google Scholar] [CrossRef]
  50. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  51. Segurado, P.; Araújo, M.B. An evaluation of methods for modelling species distributions. J. Biogeogr. 2004, 31, 1555–1568. [Google Scholar] [CrossRef]
  52. Meynard, C.N.; Quinn, J.F. Predicting species distributions: A critical comparison of the most common statistical models using artificial species. J. Biogeogr. 2007, 34, 1455–1469. [Google Scholar] [CrossRef]
  53. Bedia, J.; Busqué, J.; Gutiérrez, J.M. Predicting plant species distribution across an alpine rangeland in northern Spain. A comparison of probabilistic methods. Appl. Veg. Sci. 2011, 14, 415–432. [Google Scholar] [CrossRef] [Green Version]
  54. King, G.; Zeng, L.; Tomz, M. Logistic Regression in Rare Events Data. J. Stat. Softw. 2003, 8, 1–27. [Google Scholar] [CrossRef]
  55. Elith, J.; Graham, C.H.; Anderson, R.P.; Dudík, M.; Ferrier, S.; Guisan, A.; Hijmans, R.J.; Huettmann, F.; Leathwick, J.R.; Lehmann, A.; et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 2006, 29, 129–151. [Google Scholar] [CrossRef] [Green Version]
  56. Calcagno, V.; De Mazancourt, C. glmulti: AnRPackage for Easy Automated Model Selection with (Generalized) Linear Models. J. Stat. Softw. 2010, 34, 1–29. [Google Scholar] [CrossRef] [Green Version]
  57. Vanderwal, J.; Shoo, L.P.; Graham, C.; Williams, S.E. Selecting pseudo-absence data for presence-only distribution modeling: How far should you stray from what you know? Ecol. Model. 2009, 220, 589–594. [Google Scholar] [CrossRef]
  58. VanDerWal, J.; Lorena Falconi, L.; Januchowski, S.; Shoo, L.; Storlie, C. Species Distribution Modelling Tools: Tools for Pro-Cessing Data Associated with Species Distribution Modelling Exercises. 2014. Available online: https://www.rforge.net/SDMTools/ (accessed on 26 November 2020).
  59. Dormann, C.F.; Elith, J.; Bacher, S.; Buchmann, C.; Carl, G.; Carré, G.; Marquéz, J.R.G.; Gruber, B.; Lafourcade, B.; Leitão, P.J.; et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 2012, 36, 27–46. [Google Scholar] [CrossRef]
  60. Bolliger, J.; Kienast, F.; Bugmann, H. Comparing models for tree distributions: Concept, structures, and behavior. Ecol. Model. 2000, 134, 89–102. [Google Scholar] [CrossRef]
  61. Morin, X.; Thuiller, W. Comparing niche- and process-based models to reduce prediction uncertainty in species range shifts under climate change. Ecolology 2009, 90, 1301–1313. [Google Scholar] [CrossRef]
  62. Cheaib, A.; Badeau, V.; Boe, J.; Chuine, I.; Delire, C.; Dufrêne, E.; François, C.; Gritti, E.S.; Legay, M.; Pagé, C.; et al. Climate change impacts on tree ranges: Model intercomparison facilitates understanding and quantification of uncertainty. Ecol. Lett. 2012, 15, 533–544. [Google Scholar] [CrossRef]
  63. Takolander, A.; Hickler, T.; Meller, L.; Cabeza, M. Comparing future shifts in tree species distributions across Europe projected by statistical and dynamic process-based models. Reg. Environ. Chang. 2018, 19, 251–266. [Google Scholar] [CrossRef] [Green Version]
  64. Ramirez-Villegas, J.; Cuesta, F.; Devenish, C.; Peralvo, M.; Jarvis, A.; Arnillas, C.A. Using species distributions models for designing conservation strategies of Tropical Andean biodiversity under climate change. J. Nat. Conserv. 2014, 22, 391–404. [Google Scholar] [CrossRef] [Green Version]
  65. Fadrique, B.; Báez, S.; Duque, Á.; Malizia, A.; Blundo, C.; Carilla, J.; Osinaga-Acosta, O.; Malizia, L.; Silman, M.; Farfán-Ríos, W.; et al. Widespread but heterogeneous responses of Andean forests to climate change. Nat. Cell Biol. 2018, 564, 207–212. [Google Scholar] [CrossRef]
  66. Steinbauer, M.J.; Grytnes, J.-A.; Jurasinski, G.; Kulonen, A.; Lenoir, J.; Pauli, H.; Rixen, C.; Winkler, M.; Bardy-Durchhalter, M.; Barni, E.; et al. Accelerated increase in plant species richness on mountain summits is linked to warming. Nat. Cell Biol. 2018, 556, 231–234. [Google Scholar] [CrossRef]
  67. Heinrich, B. Why Have Some Animals Evolved to Regulate a High Body Temperature? Am. Nat. 1977, 111, 623–640. [Google Scholar] [CrossRef]
  68. McKechnie, A.E.; Wolf, B.O. Climate change increases the likelihood of catastrophic avian mortality events during extreme heat waves. Biol. Lett. 2009, 6, 253–256. [Google Scholar] [CrossRef] [Green Version]
  69. Tingley, M.W.; Monahan, W.B.; Beissinger, S.R.; Moritz, C. Birds track their Grinnellian niche through a century of climate change. Proc. Natl. Acad. Sci. USA 2009, 106, 19637–19643. [Google Scholar] [CrossRef] [Green Version]
  70. Avalos, V.D.R.; Hernández, J. Projected distribution shifts and protected area coverage of range-restricted Andean birds under climate change. Glob. Ecol. Conserv. 2015, 4, 459–469. [Google Scholar] [CrossRef]
  71. Williams, S.; Williams, Y.M.; Vanderwal, J.; Isaac, J.L.; Shoo, L.P.; Johnson, C.N. Ecological specialization and population size in a biodiversity hotspot: How rare species avoid extinction. Proc. Natl. Acad. Sci. USA 2009, 106, 19737–19741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Early, R.; Sax, D.F. Climatic niche shifts between species’ native and naturalized ranges raise concern for ecological forecasts during invasions and climate change. Glob. Ecol. Biogeogr. 2014, 23, 1356–1365. [Google Scholar] [CrossRef]
  73. Bocsi, T.; Allen, J.M.; Bellemare, J.; Kartesz, J.; Nishino, M.; Bradley, B.A. Plants’ native distributions do not reflect climatic tolerance. Divers. Distrib. 2016, 22, 615–624. [Google Scholar] [CrossRef] [Green Version]
  74. Abellán, P.; Tella, J.L.; Carrete, M.; Cardador, L.; Anadón, J.D. Climate matching drives spread rate but not establishment success in recent unintentional bird introductions. Proc. Natl. Acad. Sci. USA 2017, 114, 9385–9390. [Google Scholar] [CrossRef] [Green Version]
  75. Hartig, F.; Dyke, J.G.; Hickler, T.; Higgins, S.I.; O’Hara, R.B.; Scheiter, S.; Huth, A. Connecting dynamic vegetation models to data—An inverse perspective. J. Biogeogr. 2012, 39, 2240–2252. [Google Scholar] [CrossRef]
  76. Bloomfield, K.J.; Cernusak, L.A.; Eamus, D.; Ellsworth, D.S.; Prentice, I.C.; Wright, I.J.; Boer, M.M.; Bradford, M.G.; Cale, P.; Cleverly, J.; et al. A continental-scale assessment of variability in leaf traits: Within species, across sites and between seasons. Funct. Ecol. 2018, 32, 1492–1506. [Google Scholar] [CrossRef]
  77. Gillison, A.N. Plant functional indicators of vegetation response to climate change, past present and future: I. Trends, emerging hypotheses and plant functional modality. Flora Morphol. Distrib. Funct. Ecol. Plants 2019, 254, 12–30. [Google Scholar] [CrossRef]
  78. Modeling past plant species’ distributions in mountainous areas: A way to improve our knowledge of future climate change impacts? Past Glob. Chang. Mag. 2020, 28, 16–17. [CrossRef]
  79. Kattge, J.; Bönisch, G.; Díaz, S.; Lavorel, S.; Prentice, I.C.; Leadley, P.; Tautenhahn, S.; Werner, G.D.A.; Aakala, T.; Abedi, M.; et al. TRY plant trait database – enhanced coverage and open access. Glob. Chang. Biol. 2019, 26, 119–188. [Google Scholar] [CrossRef] [Green Version]
  80. Rojas, A.; Yucra, E.; Vera, I.; Requejo, A.; Tella, J. A new population of the globally Endangered Red-fronted Macaw Ara rubrogenys unusually breeding in palms. Bird Conserv. Int. 2012, 24, 389–392. [Google Scholar] [CrossRef] [Green Version]
  81. Blanco, G.; Morinha, F.; Roques, S.; Hiraldo, F.; Rojas, A.; Tella, J.L. Fine-scale genetic structure in the critically endangered red-fronted macaw in the absence of geographic and ecological barriers. Sci. Rep. 2021, 11, 1–17. [Google Scholar] [CrossRef] [PubMed]
  82. Spooner, F.E.B.; Pearson, R.G.; Freeman, R. Rapid warming is associated with population decline among terrestrial birds and mammals globally. Glob. Chang. Biol. 2018, 24, 4521–4531. [Google Scholar] [CrossRef] [Green Version]
Diversity 13 00094 g001 550
Figure 1. Distribution area (alpha hull) in Bolivia of Ara rubrogenys at altitudes with occurrences. 96% of the pixels over the area have an altitude between 1100 and 3000 m.
Figure 1. Distribution area (alpha hull) in Bolivia of Ara rubrogenys at altitudes with occurrences. 96% of the pixels over the area have an altitude between 1100 and 3000 m.
Diversity 13 00094 g001
Diversity 13 00094 g002 550
Figure 2. Examples of plant species distributions over Bolivia (Ara rubrogenys area delimited by orange lines) predicted by the DVM CARAIB (CARbon Assimilation In the Biosphere) for present and future climates under RCP2.6 andRCP8.5 forcings (+: occurrences of the species, light to dark color variation shows plant net productivity of the species fraction, in blue-gray between 0 and threshold of presence, in green between threshold and maximal value; the other plant species are shown in Figure S1).
Figure 2. Examples of plant species distributions over Bolivia (Ara rubrogenys area delimited by orange lines) predicted by the DVM CARAIB (CARbon Assimilation In the Biosphere) for present and future climates under RCP2.6 andRCP8.5 forcings (+: occurrences of the species, light to dark color variation shows plant net productivity of the species fraction, in blue-gray between 0 and threshold of presence, in green between threshold and maximal value; the other plant species are shown in Figure S1).
Diversity 13 00094 g002
Diversity 13 00094 g003a 550Diversity 13 00094 g003b 550
Figure 3. Area of Ara rubrogenys occupancy and fraction of the net primary productivity (fNPP) of the resource plant species predicted by the DVM CARAIB for present and future climates under the RCP2.6 andRCP8.5 forcings (Cnidoscolus tubulosus was not predicted on this area).
Figure 3. Area of Ara rubrogenys occupancy and fraction of the net primary productivity (fNPP) of the resource plant species predicted by the DVM CARAIB for present and future climates under the RCP2.6 andRCP8.5 forcings (Cnidoscolus tubulosus was not predicted on this area).
Diversity 13 00094 g003aDiversity 13 00094 g003b
Diversity 13 00094 g004 550
Figure 4. Examples of continental distribution of resource plant species predicted by the SDMs for the present (+: present occurrences of the species, light to dark color variation shows the probability of presence, in blue-grey between 0 and threshold of presence, in green between threshold and maximal value, pink color masks the area where at least one of the climate factors was out of the range of the model computation dataset, Ara rubrogenys area in orange; other plant species shown in Figure S2).
Figure 4. Examples of continental distribution of resource plant species predicted by the SDMs for the present (+: present occurrences of the species, light to dark color variation shows the probability of presence, in blue-grey between 0 and threshold of presence, in green between threshold and maximal value, pink color masks the area where at least one of the climate factors was out of the range of the model computation dataset, Ara rubrogenys area in orange; other plant species shown in Figure S2).
Diversity 13 00094 g004
Diversity 13 00094 g005a 550Diversity 13 00094 g005b 550
Figure 5. Examples of distributions of plant species and distribution of Ara rubrogenys over Bolivia predicted by SDMs for present and future climates under 2.6 and RCP8.5 forcings (+: present occurrences of the species, light to dark color variation shows the probability of presence, in blue-grey between 0 and threshold of presence, in green between threshold and maximal value, pink color masks the area where at least one of the climate factors was out of the range of the model computation dataset, Ara rubrogenys area delimited by orange lines; other plant species shown in Figure S3).
Figure 5. Examples of distributions of plant species and distribution of Ara rubrogenys over Bolivia predicted by SDMs for present and future climates under 2.6 and RCP8.5 forcings (+: present occurrences of the species, light to dark color variation shows the probability of presence, in blue-grey between 0 and threshold of presence, in green between threshold and maximal value, pink color masks the area where at least one of the climate factors was out of the range of the model computation dataset, Ara rubrogenys area delimited by orange lines; other plant species shown in Figure S3).
Diversity 13 00094 g005aDiversity 13 00094 g005b
Diversity 13 00094 g006a 550Diversity 13 00094 g006b 550
Figure 6. Area of Ara rubrogenys occupancy and mean probability of presence (Px) of A. rubrogenys and of the 18 resource plant species predicted with SDM for present and future climates under the RCP2.6 andRCP8.5 forcings.
Figure 6. Area of Ara rubrogenys occupancy and mean probability of presence (Px) of A. rubrogenys and of the 18 resource plant species predicted with SDM for present and future climates under the RCP2.6 andRCP8.5 forcings.
Diversity 13 00094 g006aDiversity 13 00094 g006b
Diversity 13 00094 g007 550
Figure 7. Proportion of pixels exhibiting a given number (0, 1, 2, etc.) of species lost (a,c) or gained (c,d) over the area of Ara rubrogenys under the RCP2.6 andRCP8.5 forcings compared to the present, for studies with the dynamic vegetation model (DVM: a,b) and the species distribution model (SDM: c,d).
Figure 7. Proportion of pixels exhibiting a given number (0, 1, 2, etc.) of species lost (a,c) or gained (c,d) over the area of Ara rubrogenys under the RCP2.6 andRCP8.5 forcings compared to the present, for studies with the dynamic vegetation model (DVM: a,b) and the species distribution model (SDM: c,d).
Diversity 13 00094 g007
Diversity 13 00094 g008a 550Diversity 13 00094 g008b 550
Figure 8. Maximal altitudes of the species over the area of Ara rubrogenys for present and future climates under the RCP2.6 andRCP8.5 forcings predicted by the dynamic vegetation model (DVM) and the species distribution model (SDM).
Figure 8. Maximal altitudes of the species over the area of Ara rubrogenys for present and future climates under the RCP2.6 andRCP8.5 forcings predicted by the dynamic vegetation model (DVM) and the species distribution model (SDM).
Diversity 13 00094 g008aDiversity 13 00094 g008b
Table
Table 1. Species occurrences and data sources (Global Diversity Information Facility site database: GBIF).
Table 1. Species occurrences and data sources (Global Diversity Information Facility site database: GBIF).
SpeciesOccurrencesSources
Anadenanthera colubrina3272GBIF.org (12 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.xccvu7
Anisocapparis speciosa218GBIF.org (12 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.xxw4c9
Ara rubrogenys159This study
Aspidosperma quebracho-blanco394GBIF.org (12 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.xccvu7
Celtis ehrenbergiana2528GBIF.org (11 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.ph4g3s
Cenchrus8069GBIF.org (12 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.3s3af9
Cnidoscolus tubulosus.159GBIF.org (6 January 2017) GBIF Occurrence Download http://doi.org/10.15468/dl.wwofy0
Jatropha hieronymi96GBIF.org (6 January 2017) GBIF Occurrence Download http://doi.org/10.15468/dl.2ziu73,
Lithraea molleoides13,271GBIF.org (12 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.m45p2c
Loxopterygium grisebachii86GBIF.org (11 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.3zkvc4
Neoraimondia herzogiana29GBIF.org (12 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.xag67z
Parasenegalia visco34GBIF.org (12 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.mnw4qm
Parkinsonia praecox1123GBIF.org (12 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.fk3as9
Prosopis chilensis81GBIF.org (6 January 2017) GBIF Occurrence Download http://doi.org/10.15468/dl.dh9ski
Prosopis kuntzei Kuntze68GBIF.org (6 January 2017) GBIF Occurrence Download http://doi.org/10.15468/dl.dh9ski
Sarcomphalus mistol217GBIF.org (11 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.mk3r74
Schinopsis marginata Engl. 41GBIF.org (6 January 2017) GBIF Occurrence Download http://doi.org/10.15468/dl.ngmuf4
Selaginella sellowii924GBIF.org (21 August 2020) GBIF Occurrence Download https://doi.org/10.15468/dl.5d5hwt
Table
Table 2. Dynamic vegetation modelling (DVM) information: sample size for statistics (N), area under the receiver operating characteristic curve (AUC), true skill statistic (TSS), threshold for TSS (net primary productivity of the species fraction: fNPP, gC/m2/y), maximal net primary productivity of the species fraction over the area (Max fNPP, gC/m2/y), sensitivity (Se), specificity (Sp).
Table 2. Dynamic vegetation modelling (DVM) information: sample size for statistics (N), area under the receiver operating characteristic curve (AUC), true skill statistic (TSS), threshold for TSS (net primary productivity of the species fraction: fNPP, gC/m2/y), maximal net primary productivity of the species fraction over the area (Max fNPP, gC/m2/y), sensitivity (Se), specificity (Sp).
Plant SpeciesNAUCTSSThresholdMax fNPPSeSp
Anadenanthera colubrina16430.675690.410830.0015318.030.94340.46744
Anisocapparis speciosa1470.604540.408164.57381131.690.959180.44898
Aspidosperma quebracho-blanco2180.65080.389910.07061064.810.944950.44495
Celtis ehrenbergiana14440.634660.477150.19821253.740.966760.51039
Cenchrus44490.668340.350420.4236930.890.825130.52529
Cnidoscolus tubulosus1450.686280.419.58461071.890.765520.63448
Jatropha hieronymi540.731310.611117.02161049.120.925930.68519
Lithraea molleoides18700.757290.611760.30591184.870.947590.66417
Loxopterygium grisebachii490.750520.571431.62231112.630.918370.65306
Neoraimondia herzogiana250.47840.40.00531111.650.96000.4400
Parasenegalia visco290.66290.4482847.26071417.990.79310.65517
Parkinsonia praecox5500.537710.3812.7344851.750.958180.42182
Prosopis chilensis760.471520.171050.02381067.330.894740.27632
Prosopis kuntzei670.623640.447761.47211009.610.910450.53731
Sarcomphalus mistol1490.637990.416110.76051206.760.959730.45638
Schinopsis marginata380.769390.5789510.20591246.940.815790.76316
Selaginella sellowii2770.725570.38628274.332910.910.660650.72563
Table
Table 3. Species distribution models (SDM) information: sample size for model estimation and statistics (N), maximal distance from presence for pseudo-absences (D, km), area under the receiver operating characteristic curve (AUC), true skill statistic (TSS), selected threshold (Thresh.), sensitivity (Se), specificity (Sp). p-value of likelihood ratio test was lower than 10−5 for each model; formula, i.e., polynomial parts of the logistic models are given in Table S2).
Table 3. Species distribution models (SDM) information: sample size for model estimation and statistics (N), maximal distance from presence for pseudo-absences (D, km), area under the receiver operating characteristic curve (AUC), true skill statistic (TSS), selected threshold (Thresh.), sensitivity (Se), specificity (Sp). p-value of likelihood ratio test was lower than 10−5 for each model; formula, i.e., polynomial parts of the logistic models are given in Table S2).
SpeciesNDAUCTSSThresh.SeSp
Ara rubrogenys3673200.97750.9430.27560.96230.9557
Anadenanthera colubrina817318000.85840.60160.16850.78230.7830
Anisocapparis speciosa58318000.94450.79870.21240.87760.8807
Aspidosperma quebracho-blanco101818000.91950.7330.21130.86090.8591
Celtis ehrenbergiana649518000.90080.68890.23630.83220.8322
Cenchrus20,71320000.77940.42570.20400.69820.6974
Cnidoscolus tubulosus46314000.81430.45270.39830.71720.7179
Jatropha hieronymi24612000.92920.78760.26050.88890.8802
Lithraea molleoides475816000.93380.74010.55550.86980.8696
Loxopterygium grisebachii22118000.97620.8720.20440.93880.9302
Neoraimondia herzogiana8314000.99030.94830.32780.96000.9655
Parasenegalia visco978000.91130.68050.29140.82760.8235
Parkinsonia praecox280718000.92030.69870.17880.84190.8418
Prosopis chilensis23820000.85300.6940.46650.84210.8395
Prosopis kuntzei20318000.92310.78570.35980.89550.8897
Sarcomphalus mistol58316000.93670.77630.24890.88000.8799
Schinopsis marginata12014000.95700.87230.37620.92110.9268
Selaginella sellowii107114000.77940.66690.34960.82670.8262
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hambuckers, A.; de Harenne, S.; Rocha Ledezma, E.; Zúñiga Zeballos, L.; François, L. Predicting the Future Distribution of Ara rubrogenys, an Endemic Endangered Bird Species of the Andes, Taking into Account Trophic Interactions. Diversity 2021, 13, 94. https://doi.org/10.3390/d13020094

AMA Style

Hambuckers A, de Harenne S, Rocha Ledezma E, Zúñiga Zeballos L, François L. Predicting the Future Distribution of Ara rubrogenys, an Endemic Endangered Bird Species of the Andes, Taking into Account Trophic Interactions. Diversity. 2021; 13(2):94. https://doi.org/10.3390/d13020094

Chicago/Turabian Style

Hambuckers, Alain, Simon de Harenne, Eberth Rocha Ledezma, Lilian Zúñiga Zeballos, and Louis François. 2021. "Predicting the Future Distribution of Ara rubrogenys, an Endemic Endangered Bird Species of the Andes, Taking into Account Trophic Interactions" Diversity 13, no. 2: 94. https://doi.org/10.3390/d13020094

APA Style

Hambuckers, A., de Harenne, S., Rocha Ledezma, E., Zúñiga Zeballos, L., & François, L. (2021). Predicting the Future Distribution of Ara rubrogenys, an Endemic Endangered Bird Species of the Andes, Taking into Account Trophic Interactions. Diversity, 13(2), 94. https://doi.org/10.3390/d13020094

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop