http://orcid.org/0000-0003-2020-2171
Figures
Abstract
Human activities cause changes to occur in the environment that affect resource availability for wildlife. The increase in the human population of cities has led to a rise in the amount of waste deposited in landfills, installations that have become a new food resource for both pest and threatened species such as vultures. In this study we used stable isotope analysis (SIA) and conventional identification of food remains from Egyptian Vultures (Neophron percnopterus) to assess the applicability of SIA as a new tool for determining the composition of the diets of vultures, a group of avian scavengers that is threatened worldwide. We focused on an expanding Egyptian Vulture population in NE Iberian Peninsula to determine the part played by landfills and livestock in the diet of these species, and aimed to reduce the biases associated with conventional ways of identifying food remains. We compared proportions of diet composition obtained with isotope mixing models and conventional analysis for five main prey. The greatest agreement between the two methods was in the categories ‘landfills’ and ‘birds’ and the greatest differences between the results from the two methods were in the categories ‘livestock’, ‘carnivores’ and ‘wild herbivores’. Despite uncertainty associated to SIA, our results showed that stable isotope analysis can help to distinguish between animals that rely on waste and so present enriched levels of δ 13C than those that feed on the countryside. Indeed, a high proportion of food derived from landfills (nearly 50%) was detected in some breeding pairs. Furthermore we performed GLMM analyses that showed that high values of δ 13C in Egyptian Vulture feathers (a proxy of feeding in landfills) are related with high levels of humanization of territories. This method has the potential to be applied to other threatened vulture species for which there is a lack of information regarding resources they are consuming, being especially important as the main causes of vultures decline worldwide are related to the consumption and availability of food resources.
Citation: Tauler-Ametller H, Hernández-Matías A, Parés F, Pretus JL, Real J (2018) Assessing the applicability of stable isotope analysis to determine the contribution of landfills to vultures’ diet. PLoS ONE 13(5): e0196044. https://doi.org/10.1371/journal.pone.0196044
Editor: Antoni Margalida, University of Lleida, SPAIN
Received: September 8, 2017; Accepted: April 5, 2018; Published: May 2, 2018
Copyright: © 2018 Tauler-Ametller 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.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: Red Electrica the España—SA (http://www.ree.es/en) financed the project, and Universitat de Barcelona (http://www.ub.edu/web/ub/ca/) supported HTA with a predoctoral grant (APIF 2014).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Human activities have greatly transformed the Earth’s land surface. Natural landscapes have been altered irremediably by human use and management practices in many human-dominated lands are constantly changing [1]. During this transformation of the landscape, humans have modified the natural resources available to wild species with serious implications for their population dynamics and distributions [2–4].
The increase in the human population in cities has led to an accumulation of large quantities of waste, which is deposited in open-air landfills scattered throughout the countryside [5]. These sites provide plentiful and predictable food resources for many generalist or scavenger species such as gulls, raven, fox and rats that are able to feed there. In some cases, these species have undergone a rapid population increase as a result [6–10]. However, resources derived from landfills are also used by other species that are not regarded as pests, some of which are classified as threatened [11,12]. Although it is known that certain species of conservation concern use landfills as feeding sites [13], few studies have quantified this use [14].
Over the last 20 years in the Iberian Peninsula the sight of large groups of vultures feeding in landfills has become commonplace [15] and it has been suggested that these installations are today a key source of food for certain vulture populations. Moreover, vultures also feed in other resources that could act as a threat, as poisoned animals. The consumption of these resources has caused in recent decades the crash of vulture populations in Asia and, more recently, in many parts of Africa [16,17]. The generalized use of veterinary drugs on livestock, pesticides, or intentional poisoning by poachers have pushed some raptor species that just two decades ago were common worldwide to the brink of extinction [18,19]. Indeed, vultures are considered today to be one of the most threatened groups of birds worldwide [20].
Despite these human-driven changes in resource distribution are carrying serious behavioural and demographic responses of vultures populations, our knowledge of the effect of these novel scenarios on vulture diets is still limited (but see [21–24]). This is partly due to complexity to quantify the net contribution in terms of ingested or assimilated biomass in this group of birds. In raptors, most conventional diet studies are based in the analysis of food remains or pellets sampled at nests or resting sites [25–27]. These method allows the identification of prey at a high taxonomic accuracy [24,28]. However, many vulture species such as those of Gyps genus may ingest large amounts of meat from corpses from big sized animals that may contribute scarcely to food sampled remains. In addition, species that use to feed on smaller species, such as Egyptian vultures, commonly carry to nests pieces from large prey while only a small fraction of it was feed; all these making challenging to stablish a correspondence between sampled food remains and ingested biomass. Consequently, the potential biases that have been recognized for conventional diet analysis, such as those linked to prey size or digestibility [27,29], might be present in available quantitative assessments of diet composition in vultures.
As an alternative, stable isotope analysis (SIA) can provide a representation of the food digested and absorbed by an animal [30]. In this field, recent Bayesian isotopic mixing models have been developed to generate potential dietary solutions for multiple dietary sources and to account for uncertainty and variation in model estimates [31], and has been widely used to infer the diet composition of a large range of organisms [32–35]. This method has been proven useful for determining the part played in animals’ diets of important food resources such as wild herbivores or food from landfills [33,36]. In spite of this, SIA has been very little used to study the diet of vultures (but see [37,38]). Thus, it would appear to be a promising tool for investigating vulture diets under the changing scenarios in which vultures live and may be useful as an indicator of on-going environmental changes.
Nevertheless, SIA is not exempt of limitations, as it requires accurate prior information regarding the trophic ecology of the studied species and results usually present a wide range of uncertainty [39,40]. Incorporating prior information obtained from the analysis of stomach in SIA has been shown to provide very reliable results [41]. However, this method requires the regurgitation of ingested food or the death of the animal, so its applicability is very limited when working with species whose individuals are difficult to capture and cannot be sampled repeatedly in time, an issue particularly concerning in endangered species. In addition, despite the great popularity of SIA in ecological studies, formal comparisons between SIA and conventional methods of diet analysis are still limited (but see [33,36]). In that context where more available and feasible methodologies present some limitations, providing quantitative assessments of diet composition of vultures based on the most easily applicable methods is crucial to better understand their potentialities and the possible differences between these methods.
This study focused on the diet of the Egyptian Vulture, a species located at the top of the food web and therefore a good indicator of mature ecosystems. This particularity, together with the fact that it is sensitive to human activities, makes it a good indicator of environmental changes [42–44]. We used both conventional diet analysis and SIA to infer the diet of this endangered vulture, we assumed that if both methodologies agree, then both methods would be providing a good approximation of resources consumed. Unlike the negative global trend in the species [45], the study population in Catalonia (NE Iberian Peninsula) has grown over the past 20 years, and has expanded and colonized new highly humanized areas from where it was not known before [46]. This process has happened in parallel with the increase in open-air landfills that has contributed to a higher probability of recruitment of breeding pairs in the area as well as to a lesser probability of territory disappearance [43,47]. The main objective of our study was to determine the part played by landfills in the diet of Egyptian Vulture breeding pairs, using conventional identification of food remains and stable isotope analysis. We thus hypothesized that breeding pairs found in humanized areas close to landfills are feeding on these facilities and therefore that the contribution of resources from waste tips in their diet will be high.
Methods
Study area and data collection
The study area lies in central eastern Catalonia (NE Iberian Peninsula, Spain) at altitudes of 200–1900 m a.s.l., and encompasses an area of cliffs running from the Prelitoral Mountains in the south to the pre-Pyrenean Mountains in the north (Fig 1). It includes humanized zones with high human population density (more than 150 hab/km2), as well as less populated rural areas (13 hab/km2). Due to this dense human population there are nine landfills within the study area that receive municipal waste, which, before being processed, is exposed and available for different opportunistic species. In 2012–2015, all occupied Egyptian Vulture territories were visited during the breeding season (March–August) by observers equipped with a spotting scope (20–60x). Once nestlings were approximately 45–55 days old, they were caught with the assistance of experienced climbers and three/four mantle feathers from each individual were taken for sampling by SIA (2012–2015; n = 60 chicks from 19 territories). Food remains for conventional diet analysis were collected in two moments during reproduction: (i) from the nests when the nestlings’ feathers were sampled and (ii) after the breeding season at the same nests (2012–2014; n = 32 samples from 18 territories). Both two samples from each nest and year were considered as one single statistical observation for diet analysis. In order to the results of conventional and isotope analysis were comparable, only data of successful breeding pairs was considered for analysis.
Grey gradient represent altitudes from 0 m (white) to 3000 m (black). Two important cities are also represented.
Conventional analysis
Food components were identified to the lowest possible taxonomic level by consulting collections in the Barcelona Museum of Natural Sciences, the Centre de Recursos de Biodiversitat Animal (Biology Faculty, Barcelona University), the collection of one of the authors (J. Real) and a reference guide [48]. Counts of items were carried out following [24]. Numbers of invertebrates and small- and medium-sized vertebrates were recorded as a minimum number of individuals. Carnivores and ungulates were counted as the number of bones or skeletal fragments as it was impossible to determine with precision the number of whole individual specimens for most species. The three distal parts of the limbs, i.e. the tarsus/carpus, metatarsus/metacarpus and digits associated with other smaller bones, were recorded as one skeletal fragment. The same applied to skull fragments, which usually consisted of several connected bones. Food components were categorized into six groups: landfill, livestock, wild herbivores, carnivores, birds and others (mainly micro-mammals, reptiles, amphibians and fish).
To differentiate between extensive livestock remains and landfill remains, cut bones or bones with traces of butchery or cooking were assumed to have been obtained from landfills [24], while the rest of the livestock remains were assumed to have been found as extensive livestock carcasses in the field. We believe this assumption is plausible, first, because due to the strict legislation and controls of intensive farms existing in Central Catalonia, it is very unlikely to find carcasses of intensive livestock abandoned in the field or near farms. Second, because despite there are some supplementary feeding points in our study area, most of them were no longer active during our study period. Also, the few that were active were mostly intended to the Bearded Vulture, so the remains that were provided were goat and sheep limbs. This type of remains are difficult to carry to the nests by Egyptian vultures as they use the beak and, consequently, the size of remains is not the optimal for this species. In addition, the period when remains were provided in these feeding points for the Bearded Vulture was from November to April or beginning of May (Catalan Government data, personal communication). In our study population chicks are born at the end of May or June, so it is very unlikely that these remains could be found in sampled nests. In the case of remains of Oryctolagus cuniculus and Sus scrofa we classified them as landfill origin when they had signs of butchery, while the rest were considered wild individuals. Regarding hair remains, white rabbits were considered domestic as there are not white rabbits in the wild and hair samples of brown rabbits were classified using a collection that allowed us to identify their origin according to their length (wild rabbits present longer hair than domestic). In the case of pigs, in our study area there is not black Iberian pig in farms so we assumed that black hair was from wild boar. When we had doubts about its origin, we classified it as undetermined.
Stable isotope analysis
Nestling feathers were first cleaned in a solution of NaOH (0.25 M) and oven-dried at 40°C for 24h. Feathers were then ground into a fine powder and subsamples of 0.35 mg were loaded into tin receptacles for combustion. Isotopic measurements were performed at the Scientific and Technological Centres of the University of Barcelona using the methods described in [49]. Stable isotope ratios are reported as δ values and expressed in ‰, according to the following equation: δ X = [(Rsample/Rstandard) − 1] X 1000, where X is 13C or 15N and R is the corresponding ratio 13C/12C or 15N/14N. Rstandard is the ratio of the international standards: Pee Dee Belemnite (PDB) for 13C and atmospheric nitrogen (AIR) for 15N. Measurement precisions for δ 13C and δ 15N were ≤ 0.15 ‰ and ≤ 0.25‰.
Isotope analysis can differentiate sources of food given that organisms reliant on plants that use a C4 photosynthesis pathway have higher (more enriched) stable carbon isotope ratios than those from food webs reliant on plants that use a C3 pathway [50]. This relationship is useful in studies of livestock species since it can differentiate between cattle raised on intensive farms that are fed with C4 plants (mainly corn) and cattle raised extensively feeding on pastures [51]. In our study area (Mediterranean and pre-Pyrenean mountains), the species that live in extensive are those aimed to meat production and, as the climate allows it, they graze in pastures during all the year without any supplement based in corn (personal communication from farmers of our study area). More, most plant species that cover pasture lands (Brachypodium, Festuca, Lolium, Agropyron, Koeleria, Oryzopsis, Molinia, Anthoxanthum, Holcus, Dactylis) use a C3 pathway [52,53] so we assumed that any C4-type isotopic signature (reflected in more enriched stable carbon isotope ratios) detected in Egyptian Vulture tissues will be due to food obtained from the landfills where humans dispose of waste from cattle bred intensively.
The isotopic signatures of prey items were obtained from the bibliography, considering main prey determined by conventional analysis and collecting values of our study area or from other regions with similar climatic conditions when possible. Landfill isotopic values were obtained from isotopic analysis of regurgitates composed of meat (mainly chicken and pork) of Larus michaellis feeding on a landfill inside our study area where Egyptian Vultures have been observed foraging [32,54]. These values were coherent with isotopic values of livestock animals raised in intensive and fed with corn in other areas [55–57]. Extensive livestock isotopic values were obtained from beef and lamb fed with C3 plants [51,58]; the wild herbivores analysed were Oryctolagus cuniculus [33] and Sus scrofa [59]. To account for the diet fraction containing carnivores, isotopic signatures of Vulpes vulpes were used [59]. Finally, the isotopic signatures of Columba palumbus and other passerines (Corvidae, Sturnidae and Turdidae) were taken from [33]. To obtain the isotopic value of categories that included more than one prey item (livestock, wild herbivores and birds), the mean and variances of 13δC and 15δN were estimated by assuming stratified random sampling and giving equal weight to each prey type in a given category [60] (Table 1). The mixing model considered that all food sources were taken into account, which means that there was no ‘others’ category.
The Trophic Enrichment Factor (TEF) is necessary to perform stable isotope analysis, as the isotope ratios of a consumer tissue are usually different from their diet due to biochemical processes when proteins are digested, absorbed, and then re-arranged for new tissue synthesis. This is the first-ever analysis of the diet of Egyptian Vultures using isotope analysis and also the first of any vulture from the Accipitridae family. Thus, there are no reference values available in the literature since Egyptian Vultures consume different resources or are taxonomically distant from other species. To solve this shortcoming, we used the SIDER package from R software, which estimates the TEF of species with unknown enrichment factors by taking into account phylogeny, tissue type sampled, general type diet (in this case we treated our study species as a carnivore), isotopic signature of the food source and measurement error [61]. To perform the SIDER model, the values of δ 13C and 15 δ N for the food sources obtained from the literature (see above) were used. TEF values obtained with SIDER model and used in mixing models were 1.11 ‰ ±1.12 for δ 13C and 3.33 ‰ ±1.18 for δ 15N.
Once we had calculated the isotopic mean values and the standard deviations of the nestlings’ feathers, prey sources and TEF, we performed Bayesian isotopic mixing models using the SIMMR package from R software [62] to estimate the relative contribution of each prey type to the diet of the Egyptian Vulture. Models were built over four Markov chains with 10000 steps per chain with a burn in of 1000 iterations. Each nest and year was considered as a single statistical observation by estimating the mean isotopic values of sampled siblings.
Statistical analyses
To assess the agreement between the conventional and stable isotope analyses we performed three different test. First, the weighted Kappa statistic (Kw) was used to compare methods on an ordinal scale by ranking prey categories from higher to lower levels of consumption. Secondly, an intra-class correlation coefficient (ICC) was used to test for the agreement of the two methods in the quantitative prey consumption estimates. We performed a third analysis that took into account uncertainty of prey estimations obtained with mixing models. We used simulated values obtained with mixing models (selecting 1000 random vales for each nest) and calculated the difference between each simulated value and the proportion of each prey obtained with conventional methods for each territory and year. After, we calculated the proportion of these difference values that were below 0 in order to study the distribution of the probability that the estimates of prey consumption by mixing models would be greater than the estimates by conventional methods. We assumed that if there were no differences between methods, the distribution of this probability for each prey category would be symmetrical and would have median values close to 0.5. As much as the distribution would be skewed towards 1 or 0 the methods would be providing different estimates of consumption for the food resource in question. Also, we reported mean differences between methods by subtracting the consumption estimates obtained with conventional analyses from the predicted consumption obtained with the SIA [36]. Only data from 2012–2014 was used to perform the comparisons as there was no information for the conventional analysis from 2015.
Finally, as δ 13C on animal tissues reflects the consumption of C3 or C4 plants and thus if they had fed on intensive (ie landfills) or extensive livestock [51], we assessed whether isotopic values of δ 13C of Egyptian vulture feathers were correlated with the availability of food and the level of humanization in our study area. We performed Generalized Linear Mixed Models (GLMM) to account for the potential non-independence of clustered observations from the same territories and years. We performed four different models, all of them with δ 13C as a response variable. The first two models considered food availability as explanatory variables, using distance to landfill as a proxy of the access of vultures to landfills resources (Model 1) and the number of cattle in extensive farming present in the territory (taken from the municipal census) and potentially available to vultures (Model 2). The other two models used urban surface (Model 3) and forest surface (Model 4) with an 8-km-radius buffer area around the nest also as explanatory variables and as a proxy of the level of humanization of territories (Table 2). We chose these variables as they can simply describe if territories are more urban or rural and so, they could provide information of the typology and distribution of food in the environment. For more details of the explanatory variables, see [43]. Territory and year were included in all four models as random factors.
All models were compared with the null model that did not consider the explanatory variable using a Likelihood Ratio Test and we report the level of significance but also AICc of each model. GLMM were performed with the lmer function from the lme4 package of R software.
Results
Diet estimates
Using conventional analyses we identified 1483 prey items corresponding to at least 62 different species. According to this methodology, the diet of the study population was composed of (mean ± SD) 30.46% (± 18.96) livestock, 19.85% (± 11.09) wild herbivores, 18.12% (± 11.64) resources from landfills, 17.45% (± 10.70) birds, 9.33% (± 13.71) carnivores and 4.79% (± 3.66) others. According to stable isotope mixing models the diet was composed of 15.52% (± 5.86) livestock, 15.82% (± 11.87) wild herbivores, 26.04% (± 5.03) resources from landfills, 19.79% (± 3.52) birds and 22.84% (± 11.88) carnivores (Fig 2). Isotope mixing models for diet inference converged.
Medians, quartiles and outliers of each prey category are represented in the boxplot.
The arithmetic mean isotopic values (±SD) for 60 Egyptian Vulture nestlings were -21.97 (± 1.088) for δ13C and 10.08 (± 1.653) for δ15N. Isotope biplots (δ13C, δ15N) showed that Egyptian Vulture nestlings lay within the space delineated by the main prey categories previously corrected by TEFs (Fig 3).
Same symbols correspond to the isotopic values of nestlings from the same territory but different years.
Individual estimates for diet according to isotope mixing models predicted a high degree of individual variability in regard to resources derived from landfills (Fig 4). Although there was a large degree of overlap in confidence intervals, the results indicate that a number of breeding pairs obtain an important proportion of their diet from landfills and that these pairs are likely to consume fewer livestock resources. Livestock in general is represented in lower proportions than landfills but still plays a significant role, especially in territories where the consumption of prey from landfills is lower (S1 Fig).
Data of each territorial pair but different year is represented with the same number. Boxes represent the credible interval of 50% and error bars the credible interval of 95% obtained with SIMMR.
Comparison of conventional methods and isotope mixing models
When ordering prey categories from higher to lower levels of consumption, the Kappa test found good agreement for ‘landfills’ (Kw = 0.564, P<0.001) and ‘birds’ (Kw = 0.255, P = 0.0129) but disagreement for the rest of prey sources. The ICC test only found agreement between prey remains counts and SIAR estimates for landfills (ICC = 0.565, P<0.001) and no significant agreement between methods for all other categories. According to the test that incorporated uncertainties of mixing models, we obtained symmetrical distribution of probabilities of differences below 0 and median values close to 0.5 for landfills (0.65) and birds (0.48) categories, but asymmetrical distribution of probabilities and median values far from 0.5 for livestock (0.07), wild herbivores (0.21) and carnivores (0.89) (Fig 5a). This results supported findings obtained with Kw and ICC test, showing that landfills and birds were the less biased categories and livestock, wild herbivores and carnivores presented more disagreement between methods. Finally, regarding the mean differences obtained for each prey using the different methods (Fig 5b) in accordance with all performed tests, there were fewest differences for the categories ‘landfill’ and ‘bird’, followed by the ‘wild herbivores’. Most disagreement was observed in the categories ‘livestock’ and ‘carnivores’.
(a) Probability of difference below 0 between stable isotope analysis and conventional analysis for main prey of Egyptian Vulture. (b) Mean differences between stable isotope analysis and conventional methods in estimates of the main prey of Egyptian Vulture. Standard errors of the differences are represented.
δ 13C, food availability and level of humanization
According to the GLMM analyses, there was no significant relationship between availability of food expressed as distance to landfill and extensive livestock availability and δ 13C of feathers of Egyptian Vulture chicks (Model 1 P = 0.0928 and Model 2 P = 0.3308). However, our results revealed a positive relationship between δ 13C and the level of humanization reported with urban surfaces (Model 3, P = 0.0009) being most urbanized areas those with higher levels of δ 13C in Egyptian vulture feathers. Consequently, δ 13C was negative correlated with rural areas reported with forest surface (Model 4, P = 0.0002) (Table 2).
Discussion
In this study we assessed the role played by landfills in the diet of Egyptian Vultures, as well as the usefulness of SIA as a tool for determining vulture diets. Both conventional analysis of food remains and SIA are the most easily applicable methods to quantify diet composition of vultures, though both of them are known to show some limitations [27,39,40,63]. So far, no comparisons have ever been made between isotope and conventional analyses of a scavenger species at the top of the food web. Here, we assumed that the agreement between methods for certain food resources consumption would support the idea that the methods are reliably approximating to the real consumption of those resources and, in this sense, we found that both types of analyses provide similar contributions of food from landfills in this species’ overall diet. Importantly, our results also suggest that landfills are a major food resource for the Egyptian Vulture, thus highlighting that human waste, besides its strong impact on populations of generalist abundant species, may be also of great importance for populations of endangered species able to exploit this type of resources.
Comparison of conventional methods and SIA
By comparing the most easily applicable methods to describe the diet of scavenger bird species we provide valuable information to understand its potential biases and constraints, an issue particularly concerning in our study group for which most available diet assessments where done by conventional diet analysis [28,64]. We grouped food resources into five categories and we considered three approaches to evaluate the agreement between methods. Our results revealed that the two methods showed good agreement for two of the considered categories: ‘landfills’ and ‘birds’. Consequently, we suggest that that isotope analysis can be used as an approximation of the consumption of food from landfills in vultures, a question of major interest in our study. In contrast, the categories that were least in agreement were ‘livestock’, ‘carnivores’ and ‘wild herbivores’. It is likely that conventional methods overestimate ‘livestock’ importance since cattle bones are usually large and easy-to-identify and remain in the nest for long periods of time [29]. In the case of wild herbivores, despite there is less disagreement between methods, it is possible that there were also overestimated by conventional methods due to the easy identification and long permanence of bones in the nest. On the other hand, it is possible that carnivores contribute substantially to the diet of Egyptian vulture as they are frequent roadkill [28,65]. However, these high levels of N detected in some chicks could be due to the use of fertilizers or animal waste products that can increase the levels of N in the environment and so arise N levels in Egyptian Vulture’s prey [66]. In spite of the observed mismatches between methods, uncertainty in parameters estimates provided by SIA were quite wide and, consequently, the percentages of diet provided by conventional methods fell within the credible intervals provided by SIA for all categories (see Fig 5).
Thus, our results show that both methods present limitations to describe de diet of the Egyptian Vulture mainly due to its wide range of prey consumed from different sources [40]. Nevertheless, the combination of both methods seems to provide complementary results. On the one hand conventional methods allow the identification of prey at a high taxonomic accuracy and this information of prey consumed is necessary to be incorporated in SIA. However, the relationship between the type and amount of a given prey that is consumed and the type and amount of a given prey that is quantified from nest remains has not been yet performed, meaning that the magnitude and direction of potential biases in conventional analysis for vulture species is unknown yet.
On the other hand, and as is usual in stable isotope studies of animal diets, our mixing models were based on several assumptions. First, we used bibliographical isotopic values of prey sources. A major constraint in isotopic analysis is to obtain representative samples of consumed prey at the time that the food resource is used, since spatial and temporal variation in isotopic values have been detected for some species [59,67]. Our study species shows a very wide dietary range and we studied it for a fairly large area and time period, so it would have been impossible to sample all potential food sources. Therefore, and according to our aims, we considered the five main sources represented by dominant species according to conventional analysis and selected those values in the literature with the highest geographic, ecological and climatic similarities. In addition, the isotopic values of food resources obtained from other areas corresponded to species or prey types that show low variation in isotopic values across large geographic ranges [68–71], which suggest that the effect of having chosen those values was small in our main results. Second, we estimated TEF values using SIDER models as we had no values of a species phylogenetically close and ecologically similar to the one we are interested in. As a result, uncertainty of TEF values was high. However, it is more correct to use TEF values with large uncertainty than use values from other species with less uncertainty but probably wrong [61]. Overall, any potential biases in our stable isotope models were incorporated as much as possible and, furthermore the models provide us the level of uncertainty in parameter estimates. In this sense, mixing models provided large credible intervals around dietary proportions that should be taken into account when interpreting results.
Despite the difficulty of describing the diet of vultures, it is showed that SIA can help to distinguish between animals that rely on waste and so present enriched levels of C than those that feed on the countryside. However it is worth to mention that in our study area although there are few supplementary feeding points, most of them are intended to Bearded Vulture so its influence to Egyptian Vulture population is predicted to be low [43, Catalan Government data, personal communication]. In this sense, in other areas where supplementary feeding points are intended to Egyptian Vulture where the remains provided could come from intensive farms, it would be difficult to distinguish between landfill and feeding points origin. Therefore, our study provides an initial step towards further research on vultures’ diet that should be refined in the future to disentangle meaningful ecological and conservation questions.
Contribution of landfills to Egyptian Vulture diet
Diet studies have been described as powerful tools for monitoring food resources obtained by raptors on large spatio-temporal scales since they are able to detect changes in resource availability caused by new environmental scenarios [27,34]. Thus, our results are consistent with these studies of other species as our diet analysis shows that an important part of the Egyptian Vulture’s diet is obtained from landfills, a novel food resource that has recently become available in the environment. The diet of Egyptian Vulture has been described to be very variable between different areas. In our case, pairs that feed most on landfills could obtain nearly a 50% of their diet from these sources, but those pairs that consume most livestock this resources represents less than 25% of their whole diet (S1 Fig). Despite the contribution of landfills to diet of this species has been little quantified because most of the studies have focused on taxonomic groups and not on the origin of these prey [24,26], it is known that waste can be important to the diet of Egyptian Vulture in some areas [72]. In addition, the proportion of livestock in our population is less than the livestock consumption described in other areas where cattle has been defined as one of the main sources of food in the Egyptian Vulture [24,64], however it is not the only exception described [25,65]. This result is relevant to conservation as vultures are among the most endangered groups of birds at a global scale and that one of their major threats comes from the type of food resources that they consume [19,73,74]. The main cause of their decline is poisoning due to the consumption of dead animals that contain high levels of pesticides or veterinary drugs such as diclofenac, which is causing a sharp reduction in vulture populations throughout the world, especially in Gyps genus [16,19,75]. Despite it is a first approximation, our results are important as they show that isotope analysis can provide information about the use of livestock resources and the potential change to consuming novel resources derived from landfills, which therefore allows us to assess the new threats to which vultures are currently exposed.
Feeding on landfills can have several conservation implications at both individual and population levels. At individual level, abundant and predictable food supplies should improve the body condition and breeding performance of the individuals that feed at these sites ([5] and references therein). In this sense, it is worth mentioning that we only sampled chicks of successful breeding pairs, so it could be possible that pairs that feed on landfills might be those with higher breeding success. Further analysis of diet of pairs that fail would be necessary to study if feeding of landfills could favour reproduction in our area. However, landfills can also act as a threat as they represent a source of food that potentially contains pollutants and poisons [76–79]. Moreover, the presence of the species in highly humanized areas where landfills are found could also entail an increase in the risk of fatal casualties due to collisions with power lines and other infrastructures [80]. Many of these threats are emerging as major causes of mortality in large avian scavengers and may turn out to be catastrophic if they operate in tandem with other threats such as poisoning [15,81]. At a population level, landfills can promote the settlement of individuals attracted by the availability of food [43] and improve survival rates and recruitment [82], but may also reduce fecundity [83]. Nevertheless, the great predictability and availability of food associated with landfills has also led to rapid population increase of some species and several studies have focused on the ecological and social consequences of this overabundance [84,85].
Moreover, results of GLMM analysis pointed out that the level of δ 13C is linked to the level of humanization of territories. High levels of δ 13C are associated with intensive livestock and so landfill consumption [36,51], therefore territories that present more δ 13C are more likely to be located in humanized areas. However, we didn’t find any relation with δ 13C and our measurement of food availability, estimated as distance to landfill and number of cattle present in the territory, possibly due to the limitations of having an accurate metric for food availability. Obtaining this metrics is particularly challenging in our study species, since breeding pairs of Egyptian Vultures make long daily movements to feed at predictable feeding sources located over 100 km away [86] and so, despite having territories located far from landfills, these pairs could be feeding on these installations.
In our study case, feeding on landfills in addition to the reduced mortality of individuals from our study area [46] could have prompted the Egyptian Vulture population increase that has taken place over the last 30 years in Catalonia. This is a paradoxical situation in which an endangered species has been favoured by human waste when, typically, modern human activities are in fact detrimental to biodiversity. Similar processes could be occurring in other threatened species but could hitherto have been overlooked. However, a new scenario is expected to occur when landfills are closed in compliance with European legislation [87], together with a resolution by the European Parliament (European Parliament resolution of 9 July 2015 on resource efficiency: moving towards a circular economy), urges the European Commission to reduce levels of residual waste to close to zero by 2020. In this context it is thus expected that the amount of waste available to vultures in coming years will fall and the effect of a landfill closure to vulture populations has been little studied [88]. Thus, this study provides a first approximation of a new tool that should be refined in order to monitor the potential diet change caused by this new resource scenario that will help design conservation measures for endangered vultures in the future.
Supporting information
S1 Fig. Percentage of the contribution of livestock in the diet of Egyptian Vulture nestlings estimated with Bayesian mixing models (SIMMR).
Data of each territorial pair but different year is represented with the same number.
https://doi.org/10.1371/journal.pone.0196044.s001
(TIF)
S1 Table. Mean ± SD ‰ values of δ13C and δ15N obtained for the Egyptian Vulture nestlings included in Bayesian mixing models.
It is specified the number of territory and the number of years that each territory was sampled (n).
https://doi.org/10.1371/journal.pone.0196044.s002
(DOCX)
S2 Table. Contribution of different food components to the diet of Egyptian Vulture nestlings obtained with conventional diet analysis.
https://doi.org/10.1371/journal.pone.0196044.s003
(DOCX)
Acknowledgments
We would like to thank the Grup de Suport de Muntanya belonging to the Cos d’Agents Rurals (Departament d’Agricultura, Generalitat de Catalunya) for help with the fieldwork and for climbing the nest-cliffs. We are also grateful for the support received from the Servei de Biodiversitat de la Generalitat de Catalunya, the natural parks of Sant Llorenç del Munt and Muntanya de Montserrat (Diputació de Barcelona) and La Garrotxa Volcanic Zone Natural Park (Generalitat de Catalunya). We would also like to thank F. Fontelles for helping with the conventional diet analysis and the Barcelona Museum of Natural Sciences for allowing us to use their reference collections. Thanks are also due to I. Peragón and the Scientific and Technological Centres of the University of Barcelona for the Stable Isotope Analysis and to J. Resano for his help with the statistical analysis. We are also grateful to Mercedes Gil from Red Eléctrica de España SA and the support received from Grup de Naturalistes d’Osona-ICHN (J. Baucells, C. Martorell). Finally, we are grateful to J. A. Donázar, S. Oppel, two anonymous reviewers and editor A. Margalida for their help in improving the manuscript. All applicable institutional and/or national guidelines for the care and use of animals were followed. We are also grateful to Mercedes Gil from Red Eléctrica de España SA.
References
- 1. Foley JA, DeFries R, Asner GP, Barford C, Bonan G, Carpenter SR, et al. Global Consequences of Land Use. Science (80-). 2005;309: 570–574. pmid:16040698
- 2. Donazar JA, Negro JJ, Hiraldo F, Hiraldo F. Foraging habitat selection, land-use changes and population decline in the lesser kestrel Falco naumanni. J Appl Ecol. 1993;30: 515–522.
- 3. Blanco G, Tella JL, Torre I. for habitats and key foraging farming Traditional in a conservation pyrrhocorax chough Pyrrhocorax landscape Spanish pseudosteppe. J Appl Ecol. 1998;35: 232–239.
- 4. Laiolo P, Dondero F, Ciliento E, Rolando A. Consequences of Pastoral Abandonment for the Structure and Diversity of the Alpine Avifauna Published by: British Ecological Society Consequences of pastoral abandonment for the structure and diversity of the alpine avifauna. J Appl Ecol. 2004;41: 294–304.
- 5. Oro D, Genovart M, Tavecchia G, Fowler MS, Martínez-Abraín A. Ecological and evolutionary implications of food subsidies from humans. Ecol Lett. 2013;16: 1501–1514. pmid:24134225
- 6. Bosch M, Oro D, Ruiz X. Dependence of yellow-legged gulls (Larus cachinnans) on food from human activity in two Western Mediterranean colonies. Avocetta. 1994;18: 135–139.
- 7. Boarman WI, Patten M a., Camp RJ, Collis SJ. Ecology of a population of subsidized predators: Common ravens in the central Mojave Desert, California. J Arid Environ. 2006;67: 248–261.
- 8. Duhem C, Roche P, Vidal E, Tatoni T. Effects of anthropogenic food resources on yellow-legged gull colony size on Mediterranean islands. Popul Ecol. 2008;50: 91–100.
- 9. Bino G, Dolev A, Yosha D, Guter A, King R, Saltz D, et al. Abrupt spatial and numerical responses of overabundant foxes to a reduction in anthropogenic resources. J Appl Ecol. 2010;47: 1262–1271.
- 10. Navarro J, Grémillet D, Afán I, Ramírez F, Bouten W, Forero MG. Feathered Detectives: Real-Time GPS Tracking of Scavenging Gulls Pinpoints Illegal Waste Dumping. PLoS One. 2016;11: 1–9. pmid:27448048
- 11. Martina A, Gallarati M. Use of a garbage dump by some mammal species in the Majella massif (Abruzzo, Italy). Hystrix, Ital J Zool. 1997;9: 23–29.
- 12. Viada C. La Milana reial. Butlletí Soc Història Nat Balear. 1994;37: 101–108.
- 13. Ciucci P, Boitani L, Francisci F, Andreoli G. Home range, activity and movements of a wolf pack in central Italy. J Zool. 1997;243: 803–819.
- 14. Ceballos O, Donázar JA. Roost-tree characteristics, food habits and seasonal abundance of roosting Egyptian Vultures in northern Spain. JRaptor Res. 1990;24: 19–25.
- 15. Martínez-Abraín A, Tavecchia G, Regan HM, Jiménez J, Surroca M, Oro D. Effects of wind farms and food scarcity on a large scavenging bird species following an epidemic of bovine spongiform encephalopathy. J Appl Ecol. 2012;49: 109–117.
- 16. Oaks JL, Gilbert M, Virani MZ, Watson RT, Meteyer CU, Rideout B a, et al. Diclofenac residues as the cause of vulture population decline in Pakistan. Nature. 2004;427: 630–633. pmid:14745453
- 17. Koenig R. Vulture Research Soars as the Scavengers ‘ Numbers Decline. Science (80-). 2006;312: 1591–1592.
- 18. Margalida A, Bogliani G, Bowden CGR, Donázar JA, Genero F, Gilbert M, et al. One Health approach to use of veterinary pharmaceuticals. Science. 2014;346: 1296–1298. pmid:25477214
- 19. Ogada D, Shaw P, Beyers RL, Buij R, Murn C, Thiollay JM, et al. Another Continental Vulture Crisis: Africa’s Vultures Collapsing toward Extinction. Conserv Lett. 2016;9: 89–97.
- 20.
IUCN Red List. The IUCN Red List of Threatened Species [Internet]. [cited 17 Jan 2018]. www.iucnredlist.org
- 21.
Mundy P, Butchart D, Ledger J, Piper S. The vultures of Africa. San Diego: Academic Press; 1992.
- 22.
Donázar JA. Los buitres ibéricos. Reyero JM, editor. Madrid; 1993.
- 23. Margalida A, Bertran J, Heredia R. Diet and food preferences of the endangered Bearded Vulture Gypaetus barbatus: Ibis (Lond 1859). 2009;151: 235–243.
- 24. Milchev B, Spassov N, Popov V. Diet of the Egyptian vulture (Neophron percnopterus) after livestock reduction in Eastern Bulgaria. North West J Zool. 2012;8: 315–323.
- 25. Margalida A, Benítez JR, Sánchez-Zapata J. A., Enrique A, Arenas R, Donázar JA. Short communication Long-term relationship between diet breadth and breeding success in a declining population of Egyptian Vultures Neophron percnopterus. Ibis (Lond 1859). 2012;154: 184–188.
- 26. Donázar JA, Cortés-Avizanda A, Carrete M. Dietary shifts in two vultures after the demise of supplementary feeding stations: consequences of the EU sanitary legislation. Eur J Wildl Res. 2010;56: 613–621.
- 27. Real J. Biases in diet study methods in the Bonelli’s eagle. J Wildl Manage. 1996;60: 632–638.
- 28. Hidalgo S, Zabala J, Zuberogoitia I, Azkona A, Castillo I. Food of the Egyptian Vulture (Neophron percnopterus) in Biscay. Buteo. 2005;14: 23–29.
- 29. Margalida A, Mañosa S, Bertran J, García D. Biases in Studying the Diet of the Bearded Vulture. J Wildl Manage. 2007;71: 1621–1625.
- 30. Hobson K a, Clark RG. Assessing avian diets using Stable Isotopes .1. Turnover of C-13 in tissues. Condor. 1992;94: 189–197.
- 31. Parnell AC, Inger R, Bearhop S, Jackson AL. Source partitioning using stable isotopes: Coping with too much variation. PLoS One. 2010;5: 1–5. pmid:20300637
- 32. Ramos R, Ramírez F, Sanpera C, Jover L, Ruiz X. Feeding ecology of yellow-legged gulls Larus michahellis in the western Mediterranean: A comparative assessment using conventional and isotopic methods. Mar Ecol Prog Ser. 2009;377: 289–297.
- 33. Resano-Mayor J, Hernández-Matías AH, Real J, Parés F, Inger R, Bearhop S. Comparing pellet and stable isotope analyses of nestling Bonelli ‘ s Eagle Aquila fasciata diet. Ibis (Lond 1859). 2014;156: 176–188.
- 34. Resano-Mayor J, Hernández-Matías A, Real J, Moleón M, Parés F, Inger R, et al. Multi-scale effects of nestling diet on breeding performance in a terrestrial top predator inferred from stable isotope analysis. PLoS One. 2014;9: e95320. pmid:24743233
- 35. Payo-Payo A, Oro D, Igual J, Jover L, Sanpera C, Tavecchia G. Population control of an overabundant species achieved through consecutive anthropogenic perturbations. Ecol Appl. 2015;25: 2228–2239. pmid:26910951
- 36. Weiser EL, Powell AN. Evaluating gull diets: A comparison of conventional methods and stable isotope analysis. J F Ornithol. 2011;82: 297–310.
- 37. Chamberlain CP, Waldbauer JR, Fox-Dobbs K, Newsome SD, Koch PL, Smith DR, et al. Pleistocene to recent dietary shifts in California condors. Proc Natl Acad Sci. 2005;102: 16707–16711. pmid:16275902
- 38. Blázquez MC, Delibes-Mateos M, Vargas JM, Granados A, Delgado A, Delibes M. Stable isotope evidence for Turkey Vulture reliance on food subsidies from the sea. Ecol Indic. 2016;63: 332–336.
- 39. Phillips DL, Inger R, Bearhop S, Jackson AL, Moore JW, Parnell AC, et al. Best practices for use of stable isotope mixing models in. Can J Zool. 2014;835: 823–835.
- 40. Robinson BG, Franke A, Derocher AE. Stable isotope mixing models fail to estimate the diet of an avian predator. Auk. 2018;135: 60–70.
- 41. Polito MJ, Trivelpiece WZ, Karnovsky NJ, Ng E, Patterson WP, Emslie SD. Integrating stomach content and stable isotope analyses to quantify the diets of pygoscelid penguins. PLoS One. 2011;6. pmid:22053199
- 42. Margalida A, García D, Avizanda AC. Factors influencing the breeding density of Bearded Vultures, Egyptian Vultures and Eurasian Griffon Vultures in Catalonia (NE Spain): management implications. Anim Biodivers Conserv. 2007;2: 189–200.
- 43. Tauler-Ametller H, Hernández-Matías A, Pretus JL, Real J. Landfills determine the distribution of an expanding breeding population of the endangered Egyptian Vulture Neophron percnopterus. Ibis (Lond 1859). 2017;159: 757–768.
- 44. Zuberogoitia I, Zabala J, Martínez J a., Martínez JE, Azkona a. Effect of human activities on Egyptian vulture breeding success. Anim Conserv. 2008;11: 313–320.
- 45.
BirdLife International. The IUCN Red List of Threatened Species. In: Neophron Percnopterus [Internet]. 2016 [cited 3 Feb 2017]. http://dx.doi.org/10.2305/IUCN.UK.2016-3.RLTS.T22695180A90559652.en.
- 46. Tauler H, Real J, Hernández-Matías A, Aymerich P, Baucells J, Martorell C, et al. Identifying key demographic parameters for the viability of a growing population of the endangered Egyptian Vulture Neophron percnopterus. Bird Conserv Int. © BirdLife International; 2015;25: 426–439.
- 47. Carrete M, Grande JM, Tella JL, Sánchez-Zapata J a., Donázar J a., Díaz-Delgado R, et al. Habitat, human pressure, and social behavior: Partialling out factors affecting large-scale territory extinction in an endangered vulture. Biol Conserv. 2007;136: 143–154.
- 48.
Brown R, Ferguson J, Lawrence M, Lees D. Huellas y señales de las aves de España y Europa. Omega E, editor. Barcelona, Spain; 2003.
- 49. Resano J, Hernández-Matías A, Real J, Parés F. Using Stable Isotopes To Determine Dietary Patterns In Bonelli’s Eagle (Aquila fasciata) Nestlings. J Raptor Res. 2011;45: 342–352.
- 50. Wooller M, Zazula G, Edwards M, Froese D, Boone R, Parker C, et al. Stable carbon isotope compositions of eastern Beringian grasses and sedges: investigating their potential as paleoenvironmental indicators. Arctic, Antarct Alp Res. 2007;39: 318–331.
- 51. Schmidt O, Quilter JM, Bahar B, Moloney AP, Scrimgeour CM, Begley IS, et al. Inferring the origin and dietary history of beef from C, N and S stable isotope ratio analysis. Food Chem. 2005;91: 545–549.
- 52.
De Bolòs O, Vigo J. Flora dels països catalans, vol I-IV. 2015.
- 53. Osborne CP, Salomaa A, Kluyver TA, Visser V, Kellogg EA, Morrone O, et al. A global database of C4 photosynthesis in grasses. New Phytol. 2014;204: 441–446. pmid:25046685
- 54.
Institut Català Ornitologia. Institut Català d’Ornitologia [Internet]. 2016. www.ornitho.cat
- 55. González-Martin I, González-Pérez C, Hernández Méndez J, Marqués-Macias E, Sanz Poveda F. Use of isotope analysis to characterize meat fromIberian-breed swine. Meat Sci. 1999;52: 437–441. pmid:22062708
- 56. Osorio MT, Moloney AP, Schmidt O, Monahan FJ. Multielement isotope analysis of bovine muscle for determination of international geographical origin of meat. J Agric Food Chem. 2011;59: 3285–3294. pmid:21391591
- 57. Heaton K, Kelly SD, Hoogewerff J, Woolfe M. Verifying the geographical origin of beef: The application of multi-element isotope and trace element analysis. Food Chem. 2008;107: 506–515.
- 58. Perini M, Camin F, Bontempo L, Rossmann A, Piasentier E. Multielement (H, C, N, O, S) stable isotope characteristics of lamb meat from different Italian regions. Rapid Commun Mass Spectrom. 2010;24: 3567–3577.
- 59. Nadjafzadeh M, Voigt CC, Krone O. Spatial, seasonal and individual variation in the diet of White-tailed Eagles Haliaeetus albicilla assessed using stable isotope ratios. Ibis (Lond 1859). 2015;158: 1–15.
- 60.
Williams BK, Nichols JD, Conroy MJ. Analysis and management of animal populations. USA: Academic Press; 2002.
- 61. Healy K, Guillerme T, Kelly S, Inger R, Bearhop S, Jackson A. SIDER: An R package for predicting trophic discrimination factors of consumers based on their ecology and phylogenetic relatedness. Ecography (Cop). 2017;
- 62.
Parnell AC. SIMMR: A Stable Isotope Mixing Model. R package version 0.3 [Internet]. 2016. https://cran.r-project.org/package=simmr
- 63. Votier SC, Bearhop S, MacCormick A, Ratcliffe N, Furness RW. Assessing the diet of great skuas, Catharacta skua, using five different techniques. Polar Biol. 2003;26: 20–26.
- 64. Donazar JA, Ceballos O. Alimentacion y tasas reproductoras del alimoche (Neophron percnopterus) en Navarra. Ardeola. 1988;35: 3–14.
- 65. Dobrev V, Boev Z, Arkumarev V, Dobrev D, Kret E, Saravia V, et al. Diet is not related to productivity but to territory occupancy in a declining population of Egyptian Vultures Neophron percnopterus. Bird Conserv Int. 2015; 1–13.
- 66. Hebert CE, Wassenaar LI. Stable nitrogen isotopes in waterfowl feathers reflect agricultural land use in western Canada. Environ Sci Technol. 2001;35: 3482–3487. pmid:11563650
- 67. Dalerum F, Angerbjörn A. Resolving temporal variation in vertebrate diets using naturally occurring stable isotopes. Oecologia. 2005;144: 647–658. pmid:16041545
- 68. Boner M, Förstel H. Stable isotope variation as a tool to trace the authenticity of beef. Anal Bioanal Chem. 2004;378: 301–310. pmid:14634708
- 69. Piasentier E, Valusso R, Camin F, Versini G. Stable isotope ratio analysis for authentication of lamb meat. Meat Sci. 2003;64: 239–247. pmid:22063009
- 70. Jose M. Moreno-Rojas VV, Lanza A, Luciano G, Ladroue V, Guillou C, Priolo A. Stable isotopes to discriminate lambs fed herbage or concentrate both obtained from C3 plants. Rapid Commun Mass Spectrom. 2010;24: 3567–3577.
- 71. Russo G, Danieli PP, Primi R, Amici A, Lauteri M. Stable isotopes in tissues discriminate the diet of free-living wild boar from different areas of central Italy. PLoS One. 2017;12: 1–11. pmid:28817721
- 72. Gangoso L, Agudo R, Anadón JD, de la Riva M, Suleyman AS, Porter R, et al. Reinventing mutualism between humans and wild fauna: insights from vultures as ecosystem services providers. Conserv Lett. 2013;6: 172–179.
- 73. Donázar J a, Margalida A, Carrete M, Sánchez-Zapata J a. Too sanitary for vultures. Science. 2009;326: 664. pmid:19900914
- 74. Hernández M, Margalida A. Poison-related mortality effects in the endangered Egyptian vulture (Neophron percnopterus) population in Spain. Eur J Wildl Res. 2009;55: 415–423.
- 75. Swan GE, Cuthbert R, Quevedo M, Green RE, Pain DJ, Bartels P, et al. Toxicity of diclofenac to Gyps vultures. Biol Lett. 2006;2: 279–282. pmid:17148382
- 76. Quarles W. Protecting Raptors from Rodenticides. Common sense pest Control Q. 2013;XXVII: 1–16.
- 77. de la Casa-Resino I, Hernández-Moreno D, Castellano A, Pérez-López M, Soler F. Breeding near a landfill may influence blood metals (Cd, Pb, Hg, Fe, Zn) and metalloids (Se, As) in white stork (Ciconia ciconia) nestlings. Ecotoxicology. 2014;23: 1377–1386. pmid:25011922
- 78. Carneiro M, Colaço B, Brandão R, Azorín B, Nicolas O, Colaço J, et al. Assessment of the exposure to heavy metals in Griffon vultures (Gyps fulvus) from the Iberian Peninsula. Ecotoxicol Environ Saf. Elsevier; 2015;113: 295–301. pmid:25528380
- 79. Casas-Díaz E, Cristòfol C, Cuenca R, Agustí S, Carneiro M, Marco I, et al. Determination of fluoroquinolone antibiotic residues in the plasma of Eurasian griffon vultures (Gyps fulvus) in Spain. Sci Total Environ. 2016;557–558: 620–626. pmid:27037883
- 80. Carrete M, Sánchez-Zapata J a., Benítez JR, Lobón M, Donázar J a. Large scale risk-assessment of wind-farms on population viability of a globally endangered long-lived raptor. Biol Conserv. 2009;142: 2954–2961.
- 81. Sanz-Aguilar A, Sánchez-Zapata JA, Carrete M, Benítez JR, Ávila E, Arenas R, et al. Action on multiple fronts, illegal poisoning and wind farm planning, is required to reverse the decline of the Egyptian vulture in southern Spain. Biol Conserv. Elsevier Ltd; 2015;187: 10–18.
- 82. Oro D, Margalida A, Carrete M, Heredia R, Donázar JA. Testing the goodness of supplementary feeding to enhance population viability in an endangered vulture. PLoS One. 2008;3. pmid:19115009
- 83. Lieury N, Gallardo M, Ponchon C, Besnard A, Millon A. Relative contribution of local demography and immigration in the recovery of a geographically-isolated population of the endangered Egyptian vulture. Biol Conserv. Elsevier B.V.; 2015;191: 349–356.
- 84. Vidal E, Medail F, Tatoni T. Is the yellow-legged gull a superabundantbird species in the Mediterranean? Impacton fauna and flora, conservation measuresand research priorities. Biodivers Conserv. 1998;1026: 1013–1026.
- 85. Bosch M, Oro D, Cantos F, Zabala M. Short term effects of culling on the ecology and population dynamics of the yellow legged gull. J Appl Ecol. 2000;37: 369–385.
- 86. López-López P, García-Ripollés C, Urios V. Food predictability determines space use of endangered vultures: implications for management of supplementary feeding. Ecol Appl. 2014;24: 938–949. pmid:25154088
- 87.
Comission E. Directive 2008/98/EC of the European Parliament and of the Council of 19 November 2008 on waste and repealing certain Directives. 2008.
- 88. Katzenberger J, Tabur E, Sen B, Isfendiyaroğlu S, Erkol IL, Oppel S. No short-term effect of closing a rubbish dump on reproductive parameters of an Egyptian Vulture population in Turkey. Bird Conserv Int. 2017; 1–12.