Figures
Abstract
The application of scientific-based conservation measures requires that sampling methodologies in studies modelling similar ecological aspects produce comparable results making easier their interpretation. We aimed to show how the choice of different methodological and ecological approaches can affect conclusions in nest-site selection studies along different Palearctic meta-populations of an indicator species. First, a multivariate analysis of the variables affecting nest-site selection in a breeding colony of cinereous vulture (Aegypius monachus) in central Spain was performed. Then, a meta-analysis was applied to establish how methodological and habitat-type factors determine differences and similarities in the results obtained by previous studies that have modelled the forest breeding habitat of the species. Our results revealed patterns in nesting-habitat modelling by the cinereous vulture throughout its whole range: steep and south-facing slopes, great cover of large trees and distance to human activities were generally selected. The ratio and situation of the studied plots (nests/random), the use of plots vs. polygons as sampling units and the number of years of data set determined the variability explained by the model. Moreover, a greater size of the breeding colony implied that ecological and geomorphological variables at landscape level were more influential. Additionally, human activities affected in greater proportion to colonies situated in Mediterranean forests. For the first time, a meta-analysis regarding the factors determining nest-site selection heterogeneity for a single species at broad scale was achieved. It is essential to homogenize and coordinate experimental design in modelling the selection of species' ecological requirements in order to avoid that differences in results among studies would be due to methodological heterogeneity. This would optimize best conservation and management practices for habitats and species in a global context.
Citation: Moreno-Opo R, Fernández-Olalla M, Margalida A, Arredondo Á, Guil F (2012) Effect of Methodological and Ecological Approaches on Heterogeneity of Nest-Site Selection of a Long-Lived Vulture. PLoS ONE 7(3): e33469. https://doi.org/10.1371/journal.pone.0033469
Editor: Brock Fenton, University of Western Ontario, Canada
Received: November 3, 2011; Accepted: February 9, 2012; Published: March 8, 2012
Copyright: © 2012 Moreno-Opo et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The European Commission cofunded the field work through the LIFENature financial program (http://ec.europa.eu/environment/life/project/Projects/index.cfm?fuseaction=search.dspPage&n_proj_id=2444). The Spanish Ministry of Environment also cofunded the field work. Funding was received only for field work but not for the meta-analysis design, statistical analysis and publication procedures. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
To date numerous studies have evaluated the relationships between one or various threatened species and the environmental variables at work in the habitats in which they carry out the distinct phases of their life cycles [1]. Of these, one of the commonest areas of study is research into the factors affecting reproductive processes, which have serious repercussions for population dynamics, and which have become one of the most important lines of work in conservation biology [2], [3]. In general, the extent of our knowledge of the factors that determine reproduction depends on natural processes and/or human activities, but it is also influenced by the methodology employed in research [4], [5]. Thus, general conclusions regarding the ecological aspects affecting the choice of reproduction sites of a single species at different spatial scales and geographical locations has only been possible in a very few cases. A possible solution is the application of a meta-analysis in order to combine results of previous studies and draw general conclusions concerning the ecological and human factors that affect habitats and species under study [6]. In order to do so, this type of analysis has to overcome the difficulties posed by the need to standardize heterogeneous information, the deficiencies in data collection in certain analyses, the lack of unifying criteria in data recording and variations in the ecological requirements of the species being modelled [7], [8].
In light of these considerations and taking as a case study the cinereous vulture (Aegypius monachus) we evaluated the factors that determine breeding habitat selection by a new field study in Spain. This species is a good model for evaluating the conservation status of the ecosystems in which it breeds given its role in trophic chains by completing the processing cycle and assimilation of biomass of dead animals [9], [10] and its sensitivity to alterations affecting the landscapes it inhabits, such as non-compatible forestry practices or human disturbances [11], [12].
From a descriptive perspective, precedent studies showed common patterns of nesting-habitat selection by cinereous vulture [13]–[19]: nests were located in forests situated on mountain slopes with large trees and high vegetation cover, far away from the human presence. Nevertheless, there are also differences across the studies in the variables that were statistically significant as well as divergences in the applied methods that could affect the final results. Therefore, we evaluated as hypothesis how the methodological procedures influenced the variability reflected in the results of different studies and, as consequence, the nesting-habitat selection of the cinereous vulture at a global scale, through meta-analysis [20]. Meta-analysis is a statistical procedure for combining data from multiple studies by applying objective formulas with the purpose of evaluating the identification and reasons of the common findings or the variation among the results of the compared studies [20].
The general objectives of this work were thus:
- to know the environmental factors determining the nesting-habitat selection of the cinereous vulture in a breeding colony of central Spain,
- to study the causes of variation in results regarding nest-site preferences from different published studies, together with the present field study, all integrated in a meta-analysis, and
- to evaluate which of the statistically significant factors highlighted in each study are the most relevant to nest-site selection in a Palearctic context and how they are related to the vulnerability and the ecological characteristics of each studied population, as a way of establishing the most appropriate management and conservation measures.
Methods
The study species
The cinereous vulture is classified as near threatened (7 200–10 000 pairs) [21] and breeds from the Iberian Peninsula as far as Eastern Asia. This vulture can be considered as an habitat indicator species due to its large foraging range [9], the specificity of its food requirements [22], [23] and its nest-site selection in large mature trees [13]. This species' habitat is located in areas with high conservation status that are important to many other species, some of which are also threatened [24]–[26].
Study area
The nest-site selection study was conducted in Alcudia and Sierra Madrona Natural Park, Spain (Figure 1), home to a colony of 129 pairs [27]. The site is part of an upland area (736–1 115 m a.s.l.) in which the dominant vegetation consists of typical Mediterranean trees such holm oak (Quercus rotundifolia), cork oak (Quercus suber), strawberry-tree (Arbutus unedo), prickly juniper (Juniperus oxycedrus) and Lusitanian oak (Quercus faginea), associated with a well-developed shrub layer.
(A). Distribution in peninsular Spain of the Special Protected Areas (SPA) with presence of breeding cinereous vultures, specifying those in which nest-site modelling were studied (dark grey) and where the present field study was developed (black) (B).
Field work and studied variables
In October-December 2005, once the breeding season was over, we visited all the cinereous vulture nests in the area (n = 155 nests). All nests, occupied in 2005 or unoccupied but with evidence of occupation in recent years, were studied [14], [28]. We also selected random points (n = 85) in areas within the perimeter of the breeding colony [15], [29]. We visited each point and recorded, from the nearest tree (since cinereous vulture breed in trees), the same data as for nests. The presence or absence of a nest at each point was established as a response variable [30].
The independent variables studied were chosen on the basis of previously evaluated aspects of this species' nest-site selection [15], [16] or as factors that relate to the land management [15]. Specifically, explanatory variables affecting factors relating to two spatial scales (microhabitat, n = 3, and landscape, n = 18) were selected.
As microhabitat variables we measured the tree characteristics in relation to the tree species (Sp_tree), the height (m) from a visual estimation of the tree where nest is present or the tree randomly selected (H_tree) and the diameter (cm) at breast height (dbh) of tree where nest is present or the tree randomly selected (D_tree).
For the landscape scale we considered 1) geomorpholocial variables as the altitude (m asl, Alt), the presence of nest/random tree in a natural scree (yes/no, Scree), the orientation of the slope where nest/random tree is located since all nests in the study area are situated in slopes greater than 15% (N, S, W, E, Orient), the slope of the hillside in a 100 m radius around the location of a nest/random tree (%, Slope) and the distance (m) from nest/random tree to nearest natural scree (D_scree), all of them calculated through GIS (ArcView 3.1 software) and aerial photographs; 2) vegetation variables, as the number of trees taller than 4 m existing in a 25 m radius around the nest/random tree through setting a survey plot and visual estimation of the height of the trees (Rad25_tree), the average high of the shrub in a 100 m radius around the nest/random tree through the measurement of the shrubs existing in four line-transects (H_shrub), the percentage of coverage in a 100 m radius around nest/random tree of trees (%_tree), shrub (%_shrub), pasturelands (%_past), scree or rock outcrop (%_scree-rock), cork oak tree (%_Qsuber), holm oak tree (%_Qrot) and other tree species (%_othersp) through a direct visual assessment in the field; and 3) human disturbance-related variables as the length (m) of unpaved tracks in a 500 m radius around the nest/random tree (Long_tracks) and the distances (m) from the nest/random tree to the nearest paved road (D_road), to the nearest building (D_const) and to the nearest unpaved track (D_track) through the application of geographic information systems (GIS, ArcView 3.1) to measure distances and length of tracks.
Climatic factors were not taken into account due to the relatively small study area (11 115 ha) and the homogeneity of vegetation structure and altitude intervals, which imply similar values of temperature, rainfall, humidity or solar radiation across the studied landscape.
Data compilation from previous studies: meta-analysis
We carried out a bibliographical search of articles published on habitat selection in the cinereous vulture in peer-reviewed journals and official reports. We were able to collate data from seven articles referring to 15 colonies in three different metapopulations (Figure 1, see Appendix S1).
We generated two response variables as means of comparison of the following null-hypothesis: 1) absence of differences in the explained variability ( = variance or deviance) among the studies (Appendix S1). The variance and deviance were selected in order to compare the variability explained and, thus, the robustness of the results from all the analyzed studies [31], [32], but taking into account that deviance is related to the number of studied variables and the sample size (Appendix S1), and 2) no differences in the frequency of a variable resulting statistically significant in the analyzed studies. We assumed that a variable might better explain the general selection patterns if being significant more times. For this later variable we also considered a) its positive/negative (+/−) relationship with the presence of cinereous vultures, and b) the proportion of times that each variable was significant in relation to the number of times it was studied.
Additionally, explanatory variables related to the data-sampling methods (Sampling methods, n = 4) and to aspects regarding the vulnerability and ecology of the species (Ecological-vulnerability, n = 3) were identified from the different studies. Sampling methods variables were intended for a better detailing of the meta-analysis and thus were included into the assessment of the effect of the methodological procedures used in each study on the explained variability. The variables considered were the type of sample considered in the analysis (point or polygon, Sample), the assignment to data set in only one year (yes or no, Year of sampling), the proportion between the number of nest-samples and number of random samples (Nest/random) and the location of random plots (1 = inside the perimeter of breeding colony; 2 = forest habitat in and around breeding colony; 3 = all habitats in and around the breeding colony, Location of random samples).
Studied breeding colonies were classified into vulnerability and ecological categories for a further analysis aiming at illustrating if any of the more significant kind of variables from those evaluated in the studies were related to the ecological and vulnerability characteristics of the studied populations. So, we identified as variables the vulnerability of the study area (1 = less than 40 pairs in the breeding colony and less than 500 pairs at national level; 2 = less than 40 pairs in the breeding colony and more than 500 pairs at national level; 3 = more than 40 pairs in the breeding colony and more than 500 pairs at national level, Threat level), the number of breeding pairs in the colony (1≤30 pairs; 2>30 pairs and ≤100 pairs; 3>100 pairs, Colony size) and the type of vegetation (pine or oak, Habitat).
Statistical analyses
Present field study.
To study the factors that affect nest-site selection, we selected a total of 240 points, that were analyzed at both microhabitat and landscape scales (Microhabitat and landscape). The analyses were performed with the software R.2.8.0 [33].
First, the variables to be included in each model were examined using Spearman's rank correlation (ρ) index to test the correlation between continuous candidate variables. Only non-correlated variables or those with weak correlation (ρ<0.3) were included in the model carried out at each scale. We did not pose a multiple contrast hypothesis and subsequent selection using information criterion (AIC or similar) [34], [35], in order to integrate the results of this study into the meta-analysis presented later in this work. All previous published papers on nesting-habitat selection of the cinereous vulture used the criteria of statistical significance, so we decided to maintain this criterion [13]–[16].
Response variables were binary (nest/random plot) and so we used generalised linear models with binomial family errors and logit-link functions. We looked for overdispersion using the dispersion parameter, which was calculated for each model by dividing the residual deviance by the residual degrees of freedom. Those models showing overdispersion were refitted by quasi-binomial family error [36].
The models were simplified by removing non-significant terms (α = 5%). Once we had determined the statistically significant factors in each model, we subsequently aggregated the non-significant levels of each factor to obtain the “minimal adequate model”, by a stepwise a posteriori procedure. If two levels of a factor did not differ significantly and did not improve the fit of the model, they were grouped together [36].
Meta-analysis.
First, we compared the variability explained by each of the studies, including the present field study results, for evaluating their ability to effectively model habitat selection from a methodological point of view. Hence, a meta-analysis testing the null-hypothesis of no-differences among the mean values of the deviance or variance was performed by using the Cochran's Q statistics of hetereogeneity. Previously, we checked that variance-deviance values of the studies fitted to a Chi-square distribution (χ2 = 136.12; df = 8; p<0.001) [37]. I-squared test quantifies the degree of heterogeneity of the studied values by analysing the percentage of the whole observed differences in the deviance-variance values between studies which are not due to chance [38]. The effect-size of the studied variable (deviance or variance) and its confidence interval at 95% were also evaluated. ANOVA (for categorical variables) and regression (for continuous variables) analyses were performed to assess the influence of sampling procedures on the deviance-variance results. All analyses were performed with the software Comprehensive Meta-Analysis V.2 [20].
On the other hand, we analysed which variables from those studied predicted to be important in habitat selection. Thus, to evaluate the frequency of appearance of each statistically significant variable regarding nest-selection we chose variables that were significant most often by selecting those resulting significant more than once from the whole studies and provided that the variable was significant more than 1/3 of times it was considered. These variables were grouped according to their relationship with microhabitat, to the ecological characteristics at landscape level (landscape), to climate, to the effect of anthropic activities (anthropic) and others.
Subsequently, a Fisher's exact test was used to analyze the proportion of appearance of each type of significant variables (landscape, microhabitat and anthropic) and the characteristics of the studied population (in terms of the threat level, colony size and habitat). The software Statistica 6.1 [39] was used to perform these analyses.
Results
Habitat modelling from present field study
In relation to microhabitat, the variables D_tree and H_tree were highly correlated (ρ = 0.58) and we decided to include the first variable in the analysis, together with Sp_tree and their first-order interaction (Appendix S2). The levels ‘cork oak’ and ‘juniper’ of the factor Sp_tree were not apparently different from each other but were with respect to the level ‘holm oak’. We checked that they could be combined into a single level without any statistically significant variation in the model by analysing changes in null deviance (χ2 = −1.29; df = 2; p = 0.041) to derive the minimal adequate model (Table 1). Vultures bred less often in holm oaks than in cork oaks and junipers. The interaction between species and tree diameter was also significant, and indicated that the diameter effect was greater in the case of holm oaks (see Table 1 for further details).
At landscape scale, we did not include in the model the variables %_Qsuber and %_ othersp, as they were correlated with %_Qrot (ρ = −0.88 and ρ = −0.75, respectively), nor %_scree-rock (correlated with D_scree ρ = −0.61) nor Long_tracks (correlated with D_track ρ = −0.75 and D_road ρ = 0.36). It was necessary to correct for overdispersion. South-facing slopes were selected for cinereous vulture to locate their nests in advance so northern, eastern and western orientation were joined obtaining a more parsimonious model from the precedent without statistical differences (F = 0.91, df = 2, p = 0.40) in order to simplify the number of levels of this variable. The model, after the simplification, it is shown in Table 1. A greater slope and the closeness of screes had a significant positive effect on nest-site selection; southern facing sites were selected in comparison with other orientations. Greater tree coverage within a radius of 100 m around the studied plot and a higher scrub layer showed a positive effect on nest presence. On the other hand, the presence of trees over 4 m high in a radius of 25 m around the nest as well as a greater cover of holm oak and shrub disfavoured selection by the cinereous vulture. Finally, the nests tended to be far from tracks, roads and human buildings (see Table 1 for further details).
Meta-analysis from published articles
In terms of sampling procedures, significant statistical differences in the deviance/variance existed between colonies (Q = 89.18, df = 12, p<0.001, I2 = 86.54, Figure 2). In addition, other factors also influenced the effect size of the deviance/variance reported in each of the analyzed studies: the location of the random plots (Q = 52.81, df = 2, p<0.001, greater deviance/variance when plots are inside the perimeter of breeding colony), whether all data were sampled in a single year (Q = 29.28, df = 1, p<0.001, greater deviance/variance when one-year sampling), the type of random data considered (Q = 47.87, df = 1, p<0.001, greater deviance/variance for points) and the ratio between the number of nests sampled and the number of random samples (Q = 89.18, df = 1, p = 0.002, greater deviance/variance when lesser nest/random ratio).
The studies references, according to Appendix S1, its location (SP = Spain, GR = Greece, GE = Georgia), the values of deviance/variance and their confidence intervals at 95% are presented.
The variables that were most often selected as statistically significant in nest-site selection for the whole Eurasian studies are presented in Table 2, being a greater slope, bigger diameter of the tree, the presence of scree around the nest-tree, orientation of the slope to south and a longer distance to tracks the five more proportionally and positively related to the presence of cinereous vulture nests. The number and type of significant variables regulating nest-site selection did not relate to the vulnerability of the colony (Table 3). Colony size conditioned the number of significant variables in relation to landscape such that in the largest colonies these variables were of greater importance (Table 3). In terms of habitat type, human activities had a greater negative effect in colonies in Mediterranean forests of Quercus sp. than in colonies located in pines (Table 3).
Discussion
Habitat selection
The present field study and the meta-analysis showed that the cinereous vulture selected nest-sites in large trees, on steep, south-facing slopes, close to screes and away from human infrastructures or anthropogenic factors than can provoke disturbances. Steeper slopes determine nest-site selection in the cinereous vulture since disturbances are less likely to occur in rugged areas [40], [41]. The largest trees cope better with the weight of nests and are easier to land on and take off from [16], [17]. This robustness and height, often found around screes [42] and in old cork oaks, as well as the positively selected situation on the upper third of a hillside, ease the detection of predators and other sources of disturbance [12], [43].
South- and east-facing slopes are selected as nest sites, probably because there are higher trees and better climate at a local scale [44]. Nonetheless, no climate-related variables were relevant in the general nest-site choice in the cinereous vulture, either due to the inter-annual heterogeneity of this feature, the fact that at local scale these variables do not predict adequately variations in ecological processes [45] or because not all the same variables were examined in the studies that we analyzed. Anthropic factors are very important in habitat selection, as occurs in other species that are wary of humans [46], [47].
Nevertheless, the impact of human activities is a complex issue that should be case-by-case evaluated in relation to each breeding colony, since the secular management practices, the degree of habituation to human presence, the social awareness to this potential beneficial species or the ecological and biogeographical conditions may nuance its real influence [48], [49]. In this sense, both vegetation and the availability and characteristics of trophic resources, as well as the impact of human activities, vary between regions approximately along a geographical gradient [50], [51], although this species does have a certain plasticity in its ecological requirements [52]. These differences can be very marked as is shown by the fact that this species breeds on cliffs in Mongolia and Russia [53], thereby demonstrating that the factors affecting life histories in a single species with a wide range can be very heterogeneous [54], [55].
Modelling and methodological conclusions
First, the number of studies analyzed (n = 13) was small to be able to draw definitive conclusions so patterns regarding the most efficient sampling approaches (see below) might be considered in relation to this low sample size. This is a common gap for meta-analysis works whose analytical procedures allow to integrate the onset of low sample sizes [32], [37]. Nevertheless, this type of analysis of a single species has never previously been conducted in such a wide geographical context [56].
The total number of variables under study influences the results of the models and the variability that they reflect [57], [58]. This implies that in habitat-modelling studies both the number and type of variables must be carefully chosen beforehand [59], [60]. In this sense, it is worth highlighting the fact that it is important to choose the most explanatory variables and/or those that are easily repeatable [4]; otherwise, there is a risk that the processes will not provide information regarding the proposed objectives. This deficiency is especially relevant in the study of climate, since many different variables are used to study the same factors [e.g. temperature, humidity, rainfall, frost and wind) and provide only scattered and inconclusive information [30], [61]. On the other hand, the influence of other factors such as the effects of global change, the use of integrating variables in ecological processes, diseases and certain biochemical factors [62], [63] could modify the tendencies of the results obtained and could lead to the application of better planned and more efficient conservation policies [64]. Our study followed the patterns of choice of variables commonly considered in previous works, both in number and type, so it was not possible to include more interesting and complex analysis in the field work, which may reveal other significant factors.
The location of the random sampling points influences the variability detected in the study. Thus, if the random plots are situated within the colony the information obtained will be more detailed in terms of factors operating at a local scale. On the other hand, if random points are chosen at a scale that includes heterogeneous types of habitat, some of the more general variables are more likely to be significant [e.g. altitude], [ slope or vegetation cover]; [ 13,57,65]. In our case we selected a priori the inner perimeter of the studied breeding colony as the framework for analysing differences between random and nests plots. One main objective was to assess differences in the selected habitat characteristics at precise level and thus, to show the sensitive factors for the breeding of cinereous vulture and to recommend the most suitable locations for developing land-use practices to local managers.
Our results reveal that if data are gathered during just one breeding season, the variability explained by the models increases, possibly because uncontrolled aspects such as inter-annual change regulating nest-site choice, weather conditions and individual behaviour are avoided [66]. Nevertheless, other studies have shown that bias could be reduced and variability in the results increased when long-term data incorporating temporary dynamics are analyzed [58], [67]. Bias-variance trade-off determines general model fitting [58], [66] and it is not possible to know exactly the bias integrated in each of the meta-analysed studies. Thus, our results should be interpreted with caution in this regard.
The election of points as random plots allows us to reflect greater deviance/variance in the model than if polygons are employed. Thus, results can throw light more accurately on questions regarding nest selection by species [16], [68].
Lastly, it is important that the proportion between the sampled nests and the random plots is as balanced as possible, although if there is an unbalance it should be in favour of the random plots. In this way, when the relationship sample/random approaches 0, the explained variability increases [57], [66], [69]. Therefore, the election of a lesser number of random plots in relation to the nests in our field study could reduce the explained deviance (Table 1).
Implications for conservation
The results obtained show that habitat type and the size of breeding colony affect the type of variables that most influences nest-site selection. Thus, human activities have greater incidence in colonies situated in Mediterranean oak forests. This may be due to the relative ease of access to the colony, to the habitat quality or to the existence of the additional conservation problems [70], [71]. In addition, it is possible that cinereous vultures may have a closer relationship with human activities in specific habitats such as pine forests, which have been exploited for a long time in a sustainable way with respect to the requirements of species present [72], [73]. In Mediterranean oak forests, on the other hand, few human economic activities are undertaken during the breeding season (except cork harvesting) and so it is possible that in these environments the species is more sensitive to human presence [43].
Despite the existence of various forest management models [73], the results of our meta-analysis suggest the need to implement different generalized management policies in temperate forests of the Palaearctic: 1) Mature forests must be given priority in forest protection as they act as source of resources and diversity [48], [72], [74]. Our results showed that areas with greater trees and tree cover are the most valuable type of forests for the breeding cinereous vultures. 2) Economic activities often determine habitat selection by threatened species [75] and so exclusion areas should be established for the most threatening activities and/or agreements should be reached to make human activities compatible [12]. Cinereous vultures tend to locate their nests as far as possible from human presence so one of the management priorities should be the regulation of such activities [71]. 3) It is advisable to coordinate and to standardize the data sampling procedures in advance when planning habitat modelling studies for the same species at different geographical scales. It is thus interesting to make the effort of developing scientific and technical working groups integrated by experts and researchers dealing with species of conservation concern [6], [74]. 4) The analysis of ecological processes that include variables that have not been taken into account to date in habitat modelling like those related to climate change, parasites-diseases or biochemical properties must be encouraged [64] and scientific evidence-based criteria must be applied on the basis of these specific studies [76].
According to these conclusions, the knowledge of habitat selection in indicator-endangered species is very valuable for optimizing evidence-based conservation actions [56], [76], [77]. Specifically, the modelling of species requirements should be undertaken for both conservation actions ex situ and in situ. Species reintroduction programmes should take into account the analysis of global patterns of habitat selection [78], [79] and so studies evaluating ecological requirements are of great relevance for carrying out population viability analyses [80], [81]. In the case of the cinereous vulture, it could be even more important given that one of the main conservation objectives for this species is the establishment of biological corridors that will connect currently isolated Palearctic populations [21], [82] through reintroduction projects (e.g., in the Pyrenees, France, Balkan Peninsula).
Supporting Information
Appendix S1.
Information contained in the articles consulted for the meta-analysis of the variation in the nest-site selection in the cinereous vulture Aegypius monachus. The variables considered for conducting the different analysis are shown (see methods for further information).
https://doi.org/10.1371/journal.pone.0033469.s001
(DOC)
Appendix S2.
Models evaluated for the analysis of nest-site selection by the cinereous vulture Aegypius monachus in the colony in Alcudia Natural Park, Spain. The scale of analysis, categorical (factors included) and continuous variables, the pairs of variables that were found to be correlated using the Spearman test (ρ), the interactions between pairs of variables included in the models, and the number of generalized linear models resulting that were analyzed are shown. Each model was built including all factors, non-correlated continuous variables and those with interactions with ρ<0.30, and only one of correlated variables each model.
https://doi.org/10.1371/journal.pone.0033469.s002
(DOC)
Acknowledgments
L. Carrascosa, Sir G. Grosvenor, F. Landaluce and J. M. Tercero kindly provided all the facilities to perform the fieldwork. L. M. González, N. El Khadir, J. Oria, R. Higuero and J. Guzmán helped in different phases of this study. Technical managers of the Castilla-La Mancha government V. Díez, I. Mosqueda and A. Aranda facilitated the permits. This work was carried out within the framework of the LIFE project 03/NAT/E/0050. The comments of two anonymous reviewers improved a previous version of the manuscript.
Author Contributions
Conceived and designed the experiments: RMO AA FG. Performed the experiments: RMO MFO AM AA FG. Analyzed the data: RMO MFO AM. Contributed reagents/materials/analysis tools: RMO MFO AM. Wrote the paper: RMO MFO AM AA FG.
References
- 1.
Morrison ML, Marcot BG, Mannan RW (2006) Wildlife-Habitat Relationships. Concepts and Applications. 3rd Edition. Washington, DC: Island Press.
- 2. Margalida A, Donázar JA, Bustamante J, Hernández F, Romero-Pujante M (2008) Application of a predictive-model to detect long-term changes in nest-site selection in the Bearded Vulture Gypaetus barbatus: conservation in relation to territory shrinkage. Ibis 150: 242–249.
- 3. López-Bao JV, Palomares F, Rodríguez A, Delibes M (2010) Effects of food supplementation on home-range size, reproductive success, productivity and recruitment in a small population of Iberian lynx. Animal Conservation 13: 35–42.
- 4. Katzner T, Millner-Guland EJ, Bragin E (2007) Using modelling to improve monitoring of structured populations: Are we collecting the right data? Conservation Biology 21: 241–252.
- 5. Magana M, Alonso JC, Martín CA, Bautista LM, Martín B (2010) Nest-site selection by Great Bustards Otis tarda suggests a trade-off between concealment and visibility. Ibis 152: 77–89.
- 6. Pullin AS, Stewart GB (2006) Guidelines for systematic review in conservation and environmental management. Conservation Biology 20: 1647–1656.
- 7. Hoeksema JD, Forde SE (2008) A Meta-Analysis of Factors Affecting Local Adaptation between Interacting Species. The American Naturalist 171: 275–290.
- 8. Benítez-López A, Alkemade R, Verweij PA (2010) The impact of roads and other infraestructures on mammal and bird populations: a meta-analysis. Biological Conservation 143: 1307–1316.
- 9. Carrete M, Donázar JA (2005) Application of central-place foraging theory shows the importance of Mediterranean dehesas for the conservation of cinereous vulture Aegypius monachus. Biological Conservation 126: 582–590.
- 10. Şekercioğlu ÇH, Daily GC, Ehrlich PR (2004) Ecosystem consequences of bird declines. Proceedings of the National Academy of Sciences USA 101: 18042–18047.
- 11. Seddon PJ, Leech T (2008) Conservation short cut, or long and winding road? A critique of umbrella species criteria. Oryx 42: 240–245.
- 12. Margalida A, Moreno-Opo R, Arroyo BE, Arredondo A (2011) Reconciling the conservation of an endangered species with an economically important anthropogenic activity: interactions between cork exploitation and the Cinereous Vulture Aegypius monachus in Spain. Animal Conservation 14: 167–174.
- 13. Morán-López R, Sánchez JM, Costillo E, Villegas A (2005) Nest-site selection of endangered cinereous vulture Aegypius monachus populations affected by anthropogenic disturbance: present and future conservation implications. Animal Conservation 9: 29–37.
- 14. Gavashelishvili A, McGrady MJ, Javakhishvili Z (2006) Planning the conservation of the breeding population of cinereous vulture Aegypius monachus in the Republic of Georgia. Oryx 40: 1–8.
- 15. Poizaridis K, Goutner V, Skartsi T, Stamou G (2004) Modelling nesting habitat as a conservation tool for the Eurasian black vulture Aegypius monachus in Dadia Nature Reserve, northeastern Greece. Biological Conservation 118: 235–248.
- 16. Donázar JA, Blanco G, Hiraldo F, Soto-Largo E, Oria J (2002) Effects of forestry and other land-use practices on the conservation of cinereous vultures. Ecological Applications 12: 1445–1456.
- 17. Fargallo JA, Blanco G, Soto-Largo E (1998) Forest management effects on nesting habitat selected by european black vultures Aegypius monachus in central Spain. Journal of Raptor Research 32: 202–207.
- 18.
Jiménez J (2002) Elaboración de un modelo de hábitat de reproducción de buitre negro Aegypius monachus en los Montes de Toledo. In: Torralvo C, editor. Anuario Ornitológico de Ciudad Real 1995–2001. Ciudad Real: SEO-Ciudad Real.
- 19.
Cuevas JA, De la Puente J (2005) Hábitat potencial del buitre negro Aegypius monachus en la Sierra de Guadarrama. Serie Documentos n° 45. Madrid: CIA Fernando González Bernáldez.
- 20.
Biostat (2006) Comprehensive Meta-Analysis. A computer program for meta-analysis. USA: Englewood.
- 21.
Birdlife International (2008) Aegypius monachus. IUCN 2010. IUCN Red List of Threatened Species. Version 2010.1. Available from www.iucnredlist.com (accesed March 2011).
- 22. Donázar JA, Margalida A, Carrete M, Sánchez-Zapata JA (2009) Too sanitary for vultures. Science 326: 664.
- 23. Moreno-Opo R, Margalida A, Arredondo A, Guil F, Martín M, et al. (2010b) Factors influencing the presence of the cinereous vulture Aegypius monachus at carcasses: food preferences and implications for the management of supplementary feeding sites. Wildlife Biology 16: 25–34.
- 24. Simberloff D (1998) Flagships, umbrellas, and keystones: Is single-species management passé in the landscape era? Biological Conservation 83: 247–257.
- 25. Palomares F, Delibes M, Ferreras P, Fedriani JM, Calzada J, et al. (2000) Iberian lynx in fragmented landscape: predispersal, dispersal and postdispersal habitats. Conservation Biology 14: 809–818.
- 26. González LM, Oria J, Sánchez R, Margalida A, Aranda A, et al. (2008) Status and habitat changes in the endangered Spanish Imperial Eagle Aquila adalberti population during 1974–2004: implications for its recovery. Bird Conservation International 18: 242–259.
- 27.
De la Puente J, Moreno-Opo R, Del Moral JC (2007) El buitre negro en España. Censo nacional 2006. Madrid: SEO/BirdLife.
- 28. Morán-López R, Sánchez JM, Costillo E, Corbacho C, Villegas A (2006) Spatial variation in anthropic and natural factors regulating the breeding success of the cinereous vulture Aegypius monachus in SW Iberian Peninsula. Biological Conservation 130: 169–182.
- 29. Sergio F, Rizzolli F, Marchesi L, Pedrini P (2004) The importance of interspecific interactions for breeding-site selection: peregrine falcons seek proximity to raven nests. Ecography 27: 818–826.
- 30. Guisan A, Zimmermann NE (2000) Predictive habitat distribution models in ecology. Ecological Modelling 135: 147–186.
- 31. Williams DA (1987) Generalized Linear Model diagnostics using the deviance and single case deletions. Applied Statistics 36: 181–191.
- 32. Clarke GM, Morris AP (2010) A Comparison of Sample Size and Power in Case-Only Association Studies of Gene-Environment Interaction. American Journal of Epidemiology 171: 498–505.
- 33.
R Development Core Team (2008) R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.
- 34.
Burnham KP, Anderson DR (1998) Model Selection and Multimodel Inference. A Practical Information-theoretic Approach. New York: Springer.
- 35.
Akaike H (1973) Information theory as an extension of the maximum likelihood principle. Second International Symposium on Information Theory. In: Petrov BN, Csaki F, editors. Budapest: Akademiai Kiado.
- 36.
Crawley MJ (2007) The R book. New York: John Wiley & Sons Ltd.
- 37. Kulinskaya E, Dollinger MB, Bjorkestol K (2011) Testing for Homogeneity in Meta-Analysis I. The One-Parameter Case: Standardized Mean Difference. Biometrics 67: 203–212.
- 38. Huedo-Medina TB, Sánchez-Meca J, Marín-Martínez F, Botella J (2006) Assessing heterogeneity in meta-analysis: Q statistic or I2 index? Psychological Methods 1: 193–206.
- 39.
StatSoft (2002) Statistica 6.1. Tulsa. Available from www.statsoft.com (accessed March 2011).
- 40. González LM, Arroyo BE, Margalida A, Oria J, Sánchez R (2006) Effect of human activities on behaviour and success of breeding Spanish imperial eagles Aquila adalberti. Animal Conservation 9: 85–93.
- 41. Aubad J, Aragon P, Rodríguez MA (2010) Human access and landscape structure effects on Andean forest bird richness. Acta Oecologica 36: 396–402.
- 42.
Thompson JD (2005) Plant evolution in the Mediterranean. Oxford: Oxford University Press.
- 43. Blumstein DT, Fernández-Juricic E, Zollner PA, Garity SC (2005) Inter-specific variation of avian responses to human disturbance. Journal of Applied Ecology 42: 943–953.
- 44. Pérez MC (2007) Diversity in the cork oak forests of the northern Straits of Gibraltar Region (Southern Spain). Plant Ecology 189: 87–100.
- 45. Hallet TB, Coulson T, Pilkington JG, Clutton-Brock TH, Pemberton JM, et al. (2004) Why large-scale climate indices seem to predict ecological processes better than local weather. Nature 430: 71–75.
- 46. Sastre P, Ponce C, Palacín C, Martín CA, Alonso JC (2009) Disturbances to great bustards Otis tarda in central Spain: human activities, bird responses and management implications. European Journal of Wildlife Research 55: 425–432.
- 47. Moreno-Opo R, Fernández-Olalla M, Guil F, Arredondo A, Higuero R, et al. (2011) The role of ponds as feeding habitat for an umbrella species: best management practices for the black stork Ciconia nigra in Spain. Oryx 45: 448–455.
- 48.
Begon M, Townsend CR, Harper JL (2006) Ecology: from individuals to ecosystems. 4th edition. Oxford: Blackwell Publishing.
- 49. Fisher J, Lindenmayer DB, Montague-Drake R (2008) The role of landscape texture in conservation biogeography: a case study on birds in south-eastern Australia. Diversity and Distributions 14: 38–46.
- 50. Doebeli M, Dieckmann U (2003) Speciation along environmental gradients. Nature 421: 259–264.
- 51. Mittelbach GG, et al. (2007) Evolution and the latitudinal diversity gradient: speciation, extinction and biogeography. Ecology Letters 10: 315–331.
- 52. Costillo E, Corbacho C, Morán R, Villegas A (2007) The diet of the black vulture Aegypius monachus in response to environmental changes in Extremadura (1970–2000). Ardeola 54: 197–204.
- 53.
Nyambayar B, Bold A, Fuller M, Watson R (2006) Nesting ecology and breeding success of cinereous vultures Aegypius monachus in central Mongolia. Proceedings of 3rd International Ornithological Conference on Modern problems of ornithology in Siberia and Central Asia, Ulan-Ude.
- 54. Morrison C, Hero JM (2003) Geographic variation in life-history characteristics of amphibians: a review. Journal of Animal Ecology 72: 270–279.
- 55. Grotan V, Saether BE, Lillegard M, Solberg EJ, Engen S (2009) Geographical variation in the influence of density dependence and climate on the recruitment of Norwegian moose. Oecologia 161: 685–695.
- 56. Branton M, Richardson JS (2011) Assessing the value of the umbrella-species concept for conservation planning with meta-analysis. Conservation Biology 25: 9–20.
- 57. Fielding AH, Bell JF (1997) A review of methods for the assessment of prediction errors in conservation presence/absence models. Environmental Conservation 24: 38–49.
- 58.
Drew CA, Wiersma Y, Huettman F (2010) Predictive species and habitat modelling in landscape ecology: concepts and applications. London: Springer.
- 59. Guisan A, Thuiller W (2005) Predicting species distribution: offering more than simple habitat models. Ecology Letters 8: 993–1009.
- 60. Jiménez-Valverde A, Lobo JM, Hortal J (2008) Not as good as they seem: the importance of concepts in species distribution modelling. Diversity and Distributions 14: 885–890.
- 61. Araujo MB, Guisan A (2006) Five (or so) challenges for species distribution modelling. Journal of Biogeography 33: 1677–1688.
- 62. Allen AP, Brown JH, Gillooly JF (2002) Global Biodiversity, Biochemical Kinetics and the Energetic-Equivalence Rule. Science 297: 1545–1548.
- 63. Hijmans RJ, Graham GH (2006) The ability of climate envelope models to predict the effect of climate change on species distributions. Global Change Biology 12: 2272–2281.
- 64. Austin M (2007) Species distribution models and ecological theory: A critical assessment and some possible new approaches. Ecological Modelling 200: 1–19.
- 65. Hirzel A, Guisan A (2002) Which is the optimal sampling strategy for habitat suitability modelling. Ecological Modelling 157: 331–341.
- 66.
Scott JM, Heglund PJ, Morrison ML (2002) Predicting species occurrences: issues of accuracy and scales. Washington, DC: Island Press.
- 67. Brotons L, Herrando S, Pla M (2007) Updating bird species distribution at large spatial scales: applications of habitat modelling to data from long-term monitoring programs. Diversity and Distributions 13: 276–288.
- 68. Segurado P, Araujo MB (2004) An evaluation of methods for modelling species distributions. Journal of Biogeography 31: 1555–1568.
- 69. Engler R, Guisan A, Rechsteiner L (2004) An improved approach for predicting the distribution of rare and endangered species from occurrence and pseudo-absence data. Journal of Applied Ecology 41: 263–274.
- 70. McIntyre S, Hobbs R (1999) A framework for conceptualizing human effects on landscapes and its relevance to management and research models. Conservation Biology 13: 1282–1292.
- 71. Fernández-Juricic E, Jiménez MD, Lucas E (2001) Alert distance as an alternative measure of bird tolerance to human disturbance: implications for park design. Environmental Conservation 28: 263–269.
- 72. Lindenmayer DB, Margules CR, Botkin DB (2000) Indicators of biodiversity for ecologically sustainable forest management. Conservation Biology 14: 941–950.
- 73. Paillet Y, Bergès L, Hjältén J, Ódor P, Avon C, et al. (2009) Biodiversity differences between managed and unmanaged forests: Meta-analysis of species richness in Europe. Conservation Biology 24: 101–112.
- 74. Harris GM, Jenkins CN, Pimm SL (2005) Refining Biodiversity Conservation Priorities. Conservation Biology 19: 1957–1968.
- 75. Martínez-Abraín A, Oro D, Jiménez J, Stewart G, Pullin A (2010) A systematic review of the effects of recreational activities on nesting birds of prey. Basic and Applied Ecology 11: 312–319.
- 76. Sutherland WJ, Pullin AW, Dolman PM, Knight TM (2004) The need for evidence-based conservation. Trends in Ecology and Evolution 19: 305–308.
- 77. Olea PP, Mateo-Tomás P (2011) Spatially explicit estimation of occupancy, detection probability and survey effort needed to inform conservation planning. Diversity and Distributions 17: 714–724.
- 78. Seddon PJ, Armstrong DP, Maloney RF (2007) Developing the science of reintroduction biology. Conservation Biology 21: 303–312.
- 79. Armstrong DP, Seddon PJ (2008) Directions in reintroduction biology. Trends in Ecology and Evolution 23: 20–25.
- 80.
IUCN (1998) Guidelines for Re-Introduction. IUCN/SSC Reintroductions Specialist Group. Gland-Cambridge: IUCN.
- 81. Sutherland WJ, Armstrong D, Butchart SHM, Earnhardt JM, Ewen J, et al. (2010) Standards for documenting and monitoring bird reintroduction projects. Conservation Letters 3: 229–235.
- 82.
Bennet AF (2003) Linkages in the Landscape: The Role of Corridors and Connectivity in Wildlife Conservation. Gland-Cambridge: IUCN.