New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Direct evidence of poison-driven widespread population decline in a wild vertebrate
Edited by Gerardo Ceballos, Universidad Nacional Autónoma de México, Mexico City, Mexico, and approved May 20, 2020 (received for review December 19, 2019)

Significance
Pesticides and other toxic substances (e.g. lead ammunition, veterinary drugs) have caused the decline of many animal populations, highlighting the serious threat that wildlife poisoning poses for biodiversity conservation. While compelling evidence demonstrates the harmful effects of poisoning on individuals, its impact on populations rests mostly on indirect evidence. This study directly links poison-induced individual mortality to the population decline of a threatened vertebrate, the red kite (Milvus milvus), across hundreds of locations in Spain, offering strong proof of poisoning contributing to countrywide population decline of a once-abundant vertebrate. A greater understanding of the population-level effects of wildlife poisoning may strengthen effective actions against the use of toxic substances harmful for wildlife.
Abstract
Toxicants such as organochlorine insecticides, lead ammunition, and veterinary drugs have caused severe wildlife poisoning, pushing the populations of several apex species to the edge of extinction. These prime cases epitomize the serious threat that wildlife poisoning poses to biodiversity. Much of the evidence on population effects of wildlife poisoning rests on assessments conducted at an individual level, from which population-level effects are inferred. Contrastingly, we demonstrate a straightforward relationship between poison-induced individual mortality and population changes in the threatened red kite (Milvus milvus). By linking field data of 1,075 poisoned red kites to changes in occupancy and abundance across 274 sites (10 × 10-km squares) over a 20-y time frame, we show a clear relationship between red kite poisoning and the decline of its breeding population in Spain, including local extinctions. Our results further support the species listing as endangered, after a breeding population decline of 31% to 43% in two decades of this once-abundant raptor. Given that poisoning threatens the global populations of more than 2,600 animal species worldwide, a greater understanding of its population-level effects may aid biodiversity conservation through increased regulatory control of chemical substances. Our results illustrate the great potential of long-term and large-scale on-ground monitoring to assist in this task.
Even though populations are a major target for both ecological risk assessments of toxicants (1, 2) and conservation management actions (3), evidence demonstrating straightforward relationships between the effects of toxic substances at individual and population levels remains scarce, especially in vertebrates (4). The paramount cases of DDT (dichlorodiphenyltrichloroethane) and other organochlorine pesticides (1, 5), lead ammunition (6), and diclofenac (7, 8) exemplify how most evidence of poisoning impacts at the population level is retrieved (SI Appendix, Appendix S1). In brief, scattered data of toxic effects from molecular to individual levels are gathered, and inferences on the relationship to observed population declines are determined through, for example, computational methods (4). For instance, strong evidence that identified diclofenac as the major cause of massive mortalities of Gyps vultures across the Indian subcontinent (7, 8) was linked with observed widespread population declines directly through deductive reasoning (7, 8) and the use of demographic simulations (9).
Inferring the effects of toxic substances at a population level from individual responses is not straightforward, as individual parameters known to be affected by toxicants (e.g., survival, fecundity) do not always correlate with population changes (2). Indeed, certain numbers of individuals can be removed from a population (e.g., poisoned) without necessarily leading to its decline due to such processes as density-dependent productivity or immigration, which may compensate for toxic effects at an individual level (1). As a result, forecasting the fate of populations affected by toxicants is a challenging task (10). Alternatively, the use of real-world population changes (e.g., through observed population growth rates), which already incorporate the complexity of population dynamics, has emerged as an ecologically sound option for assessing how toxic effects at the individual level translate into population-level impacts (2, 11). However, collecting field data at individual and population levels simultaneously is extremely costly and time-consuming, particularly at large spatiotemporal scales (10, 11). As a result, the available evidence linking individual responses to toxicants and observed population changes is limited to a few studies at local scales (e.g., one geographic location) and/or over short periods (12, 13). These spatiotemporal limitations, which in turn restrict sampling replication, weaken the strength of available evidence, and preclude drawing sound conclusions about the impact of toxicants on animal populations (11, 14).
Counteracting spatiotemporal limitations, we link field data of toxicant-induced individual mortality and population changes in the threatened red kite (Milvus milvus) across hundreds of localities in a 20-y time frame. To do so, we take advantage of two long-term and large-scale monitoring schemes: 1) Spanish national surveys of the presence and abundance of the species (∼3,500 10 × 10-km squares) in 1994 (15) and 2014 (16) and 2) data on wildlife poisoning events across Spain compiled by the ANTÍDOTO program (17) over the same period (n = 18,500 dead animals and 4,175 baits). The Spanish breeding population of red kite has decreased over the last 20 y, and poisoning, either intentional or accidental, has been suggested as one of the main causes behind this decline (16).
To assess the role of poisoning on the observed changes in the distribution and abundance of breeding red kites in 274 10 × 10-km squares between 1994 and 2014 (Fig. 1), we included the observed data on individual poisoning of red kites and other animal species (e.g., dogs, raptors; Materials and Methods and SI Appendix, Table S1) in each square as explanatory variables in generalized linear models (GLMs). Poisoning was defined here as individual intoxication, mostly by baits illegally used to kill wildlife but also by other toxic compounds used for pest and predator control, as well as anthropogenic pollutants such as lead. Land use changes known to determine the habitat suitability for red kites were also taken into consideration (18).
Distribution and abundance of the red kite breeding population per 10 × 10-km square in Spain. Data were obtained from the species censuses in 1994 (15) (A) and 2014 (16) (B). Changes in occupancy and abundance per square were obtained by combining both censuses (C). Red kite populations in the Balearic Islands and Doñana (within dashed circles) were excluded from the analyses (see Materials and Methods). The global distribution of the species according to the IUCN (24) is shown in the map at the top right. (Red kite image credit: Public Domain Pictures/George Hodan.)
Results and Discussion
A total of 1,075 red kites were registered as poisoned in mainland Spain, including 657 confirmed through toxicologic analyses and 418 suspected but without toxicological confirmation, in the 1995 to 2013 between-census period (Materials and Methods and SI Appendix, Table S2). Aldicarb and carbofuran were detected in >82% of the poisoned red kites analyzed (i.e., 316 and 228 individuals, respectively) (SI Appendix, Table S2). Despite being banned in Europe in 2003 and 2007, respectively, these two insecticides are still being used illegally to kill wildlife (17).
Poison-driven red kite mortality (both suspected and confirmed) showed a negative impact on red kite abundance (Fig. 2 A and B). Our results show a decrease in the number of breeding pairs as the number of poisoned red kites increased in a 10 × 10-km square, highlighting a pernicious effect of poisoning on the species population. This result is further supported by the negative relationship between the number of poisoned dogs, considered a good indicator of poison incidence (19) (see below), and changes in the breeding population of red kite (SI Appendix, Table S3). In contrast, while the presence of poisoned dogs (presumably by illegally poisoned baits) also negatively influenced red kite occupancy (Fig. 2C), the opposite effect was observed for poisoned red kites—that is, a greater probability of occupancy in sites with a higher number of poisoned kites (Fig. 2D). This apparently counterintuitive result may be due to the presence of dead red kites, and thus the probability of locating them is greater at locations with larger breeding populations. Red kite population size—the number of breeding pairs in 1994 and the mean number of breeding pairs in 1994 and 2014—also explains changes in occupancy and abundance, respectively (Fig. 2E and SI Appendix, Table S3). Locations with larger breeding populations seem to be more resistant to local extinction despite suffering high kite mortality, including poisoning, perhaps due to density-dependent or immigration processes. Indeed, red kite populations can expand to some extent even when exposed to toxic substances, as has been observed for some reintroduced or isolated populations (20, 21).
Main factors explaining changes in red kite occupancy and abundance in 1994 to 2014. (A) When raw data were considered, the density of dead red kites (both confirmed and suspected as poisoned) negatively correlated with the observed changes in the number of breeding pairs per 10 × 10-km square (n = 274). Although this linear relationship weakened after removing the influential point, x = 0.004 (i.e., P > 0.05), this point had no effect in the best model, further supporting the observed negative trend. (B) Dead red kites explaining the variation in the proportion of breeding pairs surviving in 2014 when included in the final best model, while keeping the remaining variables constant at their mean values; that is, locations with greater numbers of dead kites had lower proportions of surviving pairs. (C–E) In agreement with the occupancy models for breeding kites, raw data showed that locations where red kite breeding pairs disappeared between 1994 and 2014 (n = 107) had more dogs confirmed as poisoned (C), fewer dead red kites (both confirmed and suspected as poisoned) (D), and a smaller breeding population (E) in 1994. Wilcoxon tests show differences for the raw data. Wald tests for each variable once included in the final best model are provided for comparisons. Significant P values (<0.05) are in bold type. Additional details are provided in Materials and Methods and SI Appendix, Table S3. (Red kite image and dog image credit: Public Domain Pictures/George Hodan and Pixabay/AlbanyColley.)
Our results are in accordance with existing demographic models on the impact of poison on the species (13, 22) and support previous works suggesting a marked impact of poisoning on small and isolated populations (21, 23), while concurring that poisoning represents a major global threat for the species (24). Furthermore, we detected local extinctions in 107 10 × 10-km squares, underscoring the important impact of poisoning on red kite distribution and abundance. Therefore, in the absence of effective measures to eradicate or minimize poisoning, further local extinctions may occur. Along with poison, our models showed a positive association of irrigated crops on red kite presence (SI Appendix, Table S3) that could be due to irrigated farmlands harboring a more heterogeneous landscape with trees for breeding and open areas for foraging. This association is especially strong in the main species strongholds found in northwestern Spain (25). The negative impact of urban and forested areas on the species abundance (SI Appendix, Table S3) would be in agreement with the decreased availability of such heterogeneous farming landscapes in highly urbanized or forested 10 × 10-km squares.
Our results provide straightforward evidence linking the effects of poisoning on individuals (i.e., mortality) with the population trends of a wild species using field data at a fine spatiotemporal scale across hundreds of locations. The response of red kite populations to poisons suggests a strong potential for this species to act as a sentinel of toxic environmental risk in natural ecosystems (26). Interestingly, the density of dogs confirmed as poisoned per square kilometer and per year best predicted the observed changes in occupancy of the red kite breeding population (Fig. 2C and SI Appendix, Table S3). Dog owners are expected to actively look for their missing pets, thereby increasing the detectability of the poisoning cases affecting them. Considering the low detectability of wildlife poisoning events (17, 27), poisoned domestic animals emerge as a reliable index of the true incidence of toxic substances on wildlife (19) (SI Appendix, section S2).
Overall, considering the low detectability of wildlife poisoning (17, 27) and the difficulty of gathering population data over large areas (10), our results highlight poisoning as a major driver of countrywide population decline in a threatened vertebrate. Changes in distribution and abundance are major criteria used to evaluate the conservation status of a species (28). Thus, our findings further support the listing of the species as “endangered” in Spain, where the number of red kite breeding pairs has decreased by at least 31% in 20 y (16), with poisoning a seemingly major and persistent cause of this reduction (16, 28). Poisoning in all its different forms (e.g., poaching, pest and predator control, anthropogenic pollutants) threatens the global populations of at least 2,602 animal species, representing ∼3.5% of the 73,488 animal species and ∼3.5% of the 50,816 vertebrate species included on the International Union for Conservation of Nature (IUCN) Red List (29) (SI Appendix, Table S4). Poisoning threatens ∼756 species (8.4%) of the 9,013 vertebrate species listed under threatened categories (i.e., “Critically Endangered,” “Endangered,” and “Vulnerable”). Two hundred forty-five species of raptors and carnivores (28.5%) are also threatened by poisoning, including almost one-half of the species listed as “Critically Endangered,” “Endangered,” or “Vulnerable” (SI Appendix, Table S4).
Surprisingly, although ample evidence highlights the pernicious effects of poisoning on individuals (1, 4⇓⇓⇓–8), defining the impact of toxicants at high levels of biological organization (e.g., populations) remains a challenging task, resting mostly on the use of inferential methods (4). As a result, we lack straightforward evidence linking field data on the impact of toxic compounds at individual and population levels at a detailed spatiotemporal extent large enough to further support a widespread (rather than local) consistent effect (11). Direct evidence such as that provided here on the effect of poisoning on the population of a wild species at fine spatiotemporal scales across hundreds of locations strengthens the bridge between population dynamics and conservation biology (10), providing further support for taking effective actions against the use of toxic compounds pernicious to wildlife (e.g., lead ammunition, rodenticides) (4, 11). Such evidence may help improve the regulatory processes of chemical substances, including postapproval adaptation (11, 30). In this context, our data demonstrate the strong potential of existing monitoring schemes, such as the ANTÍDOTO program in Spain (17) and the SAGIR monitoring in France (27), to provide straightforward evidence of the impact of toxic substances across species over large spatiotemporal scales (11).
Materials and Methods
Study Species and Study Area.
The red kite (M. milvus) is an endemic raptor of the western Palearctic (Fig. 1), globally listed as “Near Threatened” (24). The European population, estimated as 25,200 to 33,400 breeding pairs, currently accounting for >95% of the global species population, has decreased by ∼30% over the last three generations (34.5 y), mainly due to declines within its core breeding areas in Germany, France, and Spain (26). Spain (our study area) holds roughly 2,000 breeding pairs and an estimated 50,000 wintering individuals (16), representing one of the world’s main breeding strongholds and wintering areas for red kites (31, 32). Nonetheless, the Spanish red kite wintering and breeding populations have experienced sharp decreases (7 to 19% and 31 to 43%, respectively) over the last two decades (15, 16, 33). Although the species trend in the last decade could be considered globally stable (16, 33), in Spain the red kite is legally considered “Endangered,” meaning that the species probability of survival is low as long as the primary threats to the population remain in place (34).
Red kites are opportunistic raptors that feed on a wide range of food resources, from small prey such as rodents, passerines, and young rabbits to all types of meat remains, including organic waste at dumps and carcasses of any size from large ungulates to reptiles, birds, and rodents (24). This makes the species highly susceptible to different toxic compounds, including those legally administered for pest control and deliberately poisoned baits aimed at illegal killing of predators (26, 35⇓–37). Illegal poison use is currently considered among the main threats for the species (24) and has been suggested to be a major factor behind red kite population declines in southwestern Europe (35). Mortality caused by various toxic substances is considered a factor in delaying the expansion of the Scottish reintroduced population of red kites (22), in contrast with reintroduced populations in England (20). Despite the reported high mortality caused by toxicants, these two populations are expanding, however (20, 22). In contrast, extensive anticoagulant rodenticide use in agrarian landscapes in northwestern Spain has been identified as a major driver of short-term population declines at the regional level (38). For instance, high mortality rates and demographic modeling support illegal poison use as a major cause of red kite population declines in small isolated populations in Doñana National Park and the Balearic Islands (13, 21, 23). The suggested major role of poisoning as a driver of large-scale and long-term declines of this species remains unproven, however.
Red Kite Breeding Population.
Data on the distribution and abundance of breeding red kites in mainland Spain were obtained from two national censuses performed in 1994 and 2014, the longest monitoring period available to date (15, 16). The census methodology is described in detail elsewhere (15, 39). In brief, breeding kites were located within 10 × 10-km UTM (Universal Transverse Mercator) squares in March to July of each monitoring year by trained volunteers through both car surveys at low speed (≤40 km/h) to obtain the individuals observed per km and an active search of breeding territories from vantage observation points or intensive nest searching. Data obtained by both methods were converted to population estimates using the equations shown in refs. 15 and 39.
Data on the presence and abundance of breeding red kites in 1994 were obtained for a total of 2,990 10 × 10-km squares (i.e., ∼300,000 km2; 60% of the country’s area). The species presence was confirmed in 386 squares (i.e., 38,600 km2), with an estimated total breeding population of 3,333 to 4,044 pairs (Fig. 1A) (15). In 2014, the species census covered 1,400 10 × 10-km squares (i.e., 140,000 km2), and the presence of red kites was reported in 554 squares, with an estimated population of 2,312 to 2,440 breeding pairs (Fig. 1B) (16). After correction for potential overestimation of breeding pairs detected in the last census, the Spanish breeding population was estimated as ∼2,000 breeding pairs, representing nearly a 40% decrease in its distribution over the two decades (15, 16). Data on red kite presence from both censuses were combined to obtain a response variable of the change in occupancy of breeding red kites per 10 × 10-km square in the 20-y period considered (n = 274; Fig. 1C and SI Appendix, Table S1). Those squares with breeding pairs in both censuses (i.e., 1994 and 2014) were codified as presences (value of 1; n = 167). Squares with breeding kites in 1994 but not in 2014 were considered as absences (or losses, codified as 0; n = 107). Changes in the abundance of breeding pairs among censuses were codified as a variable with two vectors: “number of succeeds,” the number of breeding pairs estimated per 10 × 10-km square in 2014, and “number of failures,” the number of breeding pairs that disappeared between 1994 and 2014 (40). Those squares with no change or an increase in the abundance of breeding pairs in the considered period were assigned 0 failures.
We excluded data from Baleares Islands and Doñana from the analyses, because these populations are geographically isolated from the rest of the Spanish red kite population and seem to present different habitat preferences (20). Nonetheless, both populations are known to have been heavily impacted by poisoning (13, 21, 23, 41).
Wildlife Poisoning Data and Other Environmental Variables.
To assess the potential influence of wildlife poisoning on the observed changes in distribution and abundance of the breeding population of red kite, we used a large and exhaustive database of confirmed or suspected cases of wildlife poisoning compiled over more than 3 decades in Spain by WWF Spain and Spanish Ornithological Society (SEO/BirdLife) (17). Most of these data were gathered by survey programs coordinated by governments of autonomous regions in Spain, and toxicologic analyses were performed at public and private ecotoxicology laboratories designated by the competent authorities. We considered 12 variables related to the incidence of poisoning in the fauna in general and in raptors and kites in particular (SI Appendix, Table S1). According to the low detectability of wildlife poisoning (estimated in 5 to 15% of the total cases) (17, 27), we also considered the number of poisoned dogs per square kilometer and year as a potentially more accurate index of the actual incidence of wildlife poisoning (19). Considering the red kite trophic ecology (i.e., feeding on small carcasses and meat remains) (24), we also calculated the density of poisoned baits as a proxy for the incidence of poisoning in natural ecosystems. Poison-related variables were calculated in two ways: considering all the episodes registered in the database in the period between censuses (i.e., 1995 to 2013; n = 18,500 animals and 4,175 baits) and considering only those episodes in which the presence of a toxic compound susceptible of poisoning wildlife was confirmed (n = 9,562 animals and 3,257 baits). All these variables were initially considered in the models.
To account for other factors able to influence the observed changes in the red kite breeding population, we considered several environmental factors previously identified as important to determine habitat suitability for the species in Spain. Seoane et al. (18) highlighted the importance of different land uses (e.g., cropland, forest, pasture), topographic variables (e.g., elevation, slope), and climatic variables (e.g., temperature) on the distribution and abundance of red kites. Among these variables considered important for the red kite in 1994, we incorporated into our analyses those susceptible of meaningful and noticeable changes in 1994 to 2014. We used CORINE Land Cover information from 1990 and 2012 (42) to calculate the percent change in crop, forest, and pasture lands, as well as urban surfaces, per 10 × 10-km square among the considered red kite censuses (SI Appendix, Table S1). We assumed that the remaining variables considered important for explaining habitat suitability for the red kite (e.g., climate, topography) had not changed noticeably in the period analyzed.
We considered the number of breeding pairs of red kite estimated per 10 × 10-km square in 1994 and the geometric mean of the breeding pairs estimated in 1994 and 2014 as a way to better determine red kite abundance during the considered period (SI Appendix, Table S1). Squares with a higher breeding population could be more resilient to species threats such as poisoning and land use change and thereby have a higher probability of remaining occupied. Red kite population sizes were used also to control the potential effect of population size on the number of poisoned kites recorded per square. The species population size was included in the predictive models in two ways, as an independent explanatory variable and as a divisor of the number of poisoned red kites in each square, to obtain an explanatory variable of poisoned kites per breeding pair (SI Appendix, Table S1).
Analysis.
We used GLMs with binomial error distributions and “logit” links to explain the observed changes in the distribution and abundance of breeding red kites between 1994 and 2014. From the significantly correlated explanatory variables (i.e., those with a Spearman correlation P < 0.05), we included in the models the one with the greater ecological significance or, if not easily discriminated, the variable with best results in a univariate GLM (e.g., lower corrected Akaike information criterion [AICc]) (43).
We used multimodel selection to choose the best models (43). We randomly selected 70% of the data to train the models and the remaining 30% to test the results. We used the training dataset to generate a set of models with all possible combinations of the considered explanatory variables. We tested both linear and quadratic responses, as well as interactions between pairs of variables. The resulting models were ranked by lower AICc and higher relative weight (ωm), which indicates the probability that a model was the best one among all the candidate models evaluated. From these ranked models, we selected a set of best models comprising those whose weights added to 0.95 (Σωm = 0.95). From the set of best models, we depleted the redundant ones, that is, those including the same variables of other more explicative models (i.e., with lower AICc and higher ωm) but with some additional variable that complicates the model without improving its predictive ability. Once filtered, model weights were recalculated to sum to 1. The models resulting from this selection process were combined into a single model considered the final best model (43, 44). The coefficients of the explanatory variables included in the final best model were calculated by weighting each variable coefficient by the weight of the model in which the variable was included.
The relative contribution of each variable to the final best model (ωi) was calculated by summing the weights of the models (ωm) within the set of best models (Σωm = 0.95) that included the considered variable. The most important variables were those with the highest relative weight (ωi = Σωm). We expected that our approach of selecting noncorrelated variables to include in the models would minimize multicollineality, thereby minimizing the issues associated with pondered models (44). Nonetheless, to obtain additional information on the performance of the explanatory variables, we also calculated Z values and P values using the Wald test (SI Appendix, Table S3) for each variable included in the final best model. All statistical analyses were performed in R (45).
Global Assessment of Poisoning as a Species Threat.
We searched the IUCN Red List (29) for animal species threatened by poisoning at global scale according to the IUCN criteria for the inclusion of threats. We considered a species to be threatened by poisoning when threat category “5.1.2: Unintentional Effects (species is not the target)”, “5.1.3: Persecution/Control,” or “9.3.3: Herbicides and Pesticides” was recorded in the IUCN’s assessment of the species. We restricted our search to only global assessments of animal species (i.e., kingdom Animalia). We also considered these figures for vertebrate species (i.e., phylum Chordata) and for raptors (i.e., orders Accipitriformes, Cathartiformes, Falconiformes, and Strigiformes) and carnivores (i.e., order Carnivora), as predators are either frequently prosecuted through illegal poisoning or affected by secondary poisoning (e.g., consuming poisoned prey, such as rodents) (1, 4⇓⇓⇓–8, 27, 35⇓⇓–38, 46, 47). The main results of this search are described in detail in SI Appendix, Table S4.
Acknowledgments
Hundreds of volunteers and personnel from Spanish autonomous governments (Comunidades Autonomas [CCAA]) participated in the national census of red kites, organized by SEO/Birdlife, with financial support from the Royal Society for the Protection of Birds, the Spanish Ministry of Environment (SME), and CCAA. The database of poisoned animals was compiled by CCAA and SME and coordinated by WWF Spain and SEO/Birdlife through the ANTÍDOTO program. This work was funded by WWF Spain under project UCLM-UCTR170245, and by the BBVA Foundation (“Ayudas Fundación BBVA a Equipos de Investigación Científica 2018”) through the TÓXICO project.
Footnotes
↵1P.M.-T. and P.P.O. contributed equally to this work.
- ↵2To whom correspondence may be addressed. Email: rktespejos{at}gmail.com.
Author contributions: P.M.-T., P.P.O., and J.V. designed research; P.M.-T. and P.P.O. performed research; P.M.-T., E.M., R.M., and J.V. contributed new reagents/analytic tools; P.M.-T. and P.P.O. analyzed data; and P.M.-T., P.P.O., R.M., and J.V. wrote the paper.
The authors declare no competing interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.1922355117/-/DCSupplemental.
Published under the PNAS license.
References
- 1.↵
- C. H. Walker,
- R. M. Sibly,
- S. P. Hopkin,
- D. B. Peakall
- 2.↵
- 3.↵
- M. E. Soulé
- 4.↵
- H. R. Köhler,
- R. Triebskorn
- 5.↵
- I. Newton
- 6.↵
- M. E. Finkelstein et al.
- 7.↵
- 8.↵
- 9.↵
- 10.↵
- D. Oro
- 11.↵
- F. Macini,
- B. A. Woodcock,
- N. J. B. Isaac
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- J. Viñuela,
- R. Martí,
- A. Ruíz
- 16.↵
- B. Molina
- 17.↵
- C. Cano,
- D. De la Bodega,
- P. Ayerza,
- E. Mínguez
- 18.↵
- J. Seoane,
- J. Viñuela,
- R. Díaz-Delgado,
- J. Bustamante
- 19.↵
- 20.↵
- F. M. Molenaar et al.
- 21.↵
- G. Tavecchia,
- J. Adrover,
- A. M. Navarro,
- R. Pradel
- 22.↵
- J. Smart et al.
- 23.↵
- F. Sergio et al.
- 24.↵
- BirdLife International
- 25.↵
- F. J. González-Estébanez,
- S. García-Tejero,
- P. Mateo-Tomás,
- P. P. Olea
- 26.↵
- C. Márquez,
- J. M. Vargas,
- R. Villafuerte,
- J. E. Fa
- 27.↵
- 28.↵
- International Union for Conservation of Nature Standards and Petitions Subcommittee
- 29.↵
- International Union for Conservation of Nature
- 30.↵
- S. Krimsky
- 31.↵
- J. Knott,
- P. Newbery,
- B. Barov
- 32.↵
- I. Carter
- 33.↵
- I. E. Cardiel
- 34.↵Spanish Official Gazette (BOE). Spanish Royal Decree on Threatened Wild Species and Wild Species Under Special Protection. BOE-A-2011-3582 [in Spanish]. https://www.boe.es/buscar/pdf/2011/BOE-A-2011-3582-consolidado.pdf. Accessed 10 April 2020.
- 35.↵
- R. Villafuerte,
- J. Viñuela,
- J. C. Blanco
- 36.↵
- 37.↵
- 38.↵
- I. Zuberogoitia,
- J. E. Martínez
- F. Mougeot,
- J. T. García,
- J. Viñuela,
- 39.↵
- J. Viñuela
- 40.↵
- M. J. Crawley
- 41.↵
- F. Sergio et al.
- 42.↵Copernicus, CORINE Land Cover European inventory, 1990 and 2012. https://land.copernicus.eu/pan-european/corine-land-cover. Accessed August 2, 2019.
- 43.↵
- K. P. Burnham,
- D. R. Anderson
- 44.↵
- 45.↵
- R Core Development Team
- 46.↵
- P. Mateo-Tomás,
- P. P. Olea,
- I. S. Sánchez-Barbudo,
- R. Mateo
- 47.↵
- D. L. Ogada
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Biological Sciences
- Environmental Sciences