Evolutionary age and risk of extinction in the global avifauna

Species at high risk of extinction are not distributed at random among higher taxa. Here we demonstrate that there is a positive relationship between the proportion of species in a taxon which are considered to be threatened and the evolutionary age of that taxon, both for the global avifauna and the avifauna of the New World. The potential mechanisms and consequences of the relationship are examined.     

Phylogenetic correlates of extinction risk in mammals: species in older lineages are not at greater risk

Phylogenetic information is becoming a recognized basis for evaluating conservation priorities, but associations between extinction risk and properties of a phylogeny such as diversification rates and phylogenetic lineage ages remain unclear. Limited taxon-specific analyses suggest that species in older lineages are at greater risk. We calculate quantitative properties of the mammalian phylogeny and model extinction risk as an ordinal index based on International Union for Conservation of Nature Red List categories. We test for associations between lineage age, clade size, evolutionary distinctiveness and extinction risk for 3308 species of terrestrial mammals. We show no significant global or regional associations, and three significant relationships within taxonomic groups. Extinction risk increases for evolutionarily distinctive primates and decreases with lineage age when lemurs are excluded. Lagomorph species (rabbits, hares and pikas) that have more close relatives are less threatened. We examine the relationship between net diversification rates and extinction risk for 173 genera and find no pattern. We conclude that despite being under-represented in the frequency distribution of lineage ages, species in older, slower evolving and distinct lineages are not more threatened or extinction-prone. Their extinction, however, would represent a disproportionate loss of unique evolutionary history.     

Estimation of capture probabilities using generalized estimating equations and mixed effects approaches

Modeling individual heterogeneity in capture probabilities has been one of the most challenging tasks in capture–recapture studies. Heterogeneity in capture probabilities can be modeled as a function of individual covariates, but correlation structure among capture occasions should be taking into account. A proposed generalized estimating equations (GEE) and generalized linear mixed modeling (GLMM) approaches can be used to estimate capture probabilities and population size for capture–recapture closed population models. An example is used for an illustrative application and for comparison with currently used methodology. A simulation study is also conducted to show the performance of the estimation procedures. Our simulation results show that the proposed quasi‐likelihood based on GEE approach provides lower SE than partial likelihood based on either generalized linear models (GLM) or GLMM approaches for estimating population size in a closed capture–recapture experiment. Estimator performance is good if a large proportion of individuals are captured. For cases where only a small proportion of individuals are captured, the estimates become unstable, but the GEE approach outperforms the other methods.     

Linear analysis solves two puzzles in population dynamics: the route to extinction and extinction in coloured environments

In this paper, we give simple explanations to two unsolved puzzles that have emerged in recent theoretical studies in population dynamics. First, the tendency of some model populations to go extinct from high population densities, and second, the positive effect of autocorrelated environments on extinction risks for some model populations. Both phenomena are given general explanations by simple, linear, sto-chastic models. We emphasize the predictive and explanatory power of such models.     

种群生存力分析（PVA）的方法与应用

随着人们对资源的加速利用,生境丧失和破碎化导致物种濒危问题日益严重.以岛屿生物地理学为理论起源的种群生存力分析（PVA）,通过分析和模拟种群动态过程并建立灭绝概率与种群数量之间的关系,为濒危物种保护提供了重要的理论依据和研究途径.在过去的几十年中,种群生存力分析已成为保护生物学中一项重要的研究内容.目前种群生存力分析发展稳定,但对其实际预测能力和准确性尚存质疑,应用方面也有待进一步发展.种群生存力分析的进一步完善还需要在理论和方法上的创新,特别是籍于景观生态学和可持续性科学的理念,将空间分析手段、经济社会因素纳入到物种和种群的预测和管理上,从而使其 具有更完整的理论基础和更高的实用价值.为此,本文对种群生存力分析的历史、基本概念、研究方法、模型应用和准确性进行了综述,并提出了有关的研究展望.     

Measurement error and estimates of population extinction risk

It is common to estimate the extinction probability for a vulnerable population using methods that are based on the mean and variance of the long‐term population growth rate. The numerical values of these two parameters are estimated from time series of population censuses. However, the proportion of a population that is registered at each census is typically not constant but will vary among years because of stochastic factors such as weather conditions at the time of sampling. Here, we analyse how such sampling errors influence estimates of extinction risk and find sampling errors to produce two opposite effects. Measurement errors lead to an exaggerated overall variance, but also introduce negative autocorrelations in the time series (which means that estimates of annual growth rates tend to alternate in size). If time series data are treated properly these two effects exactly counter balance. We advocate routinely incorporating a measure of among year correlations in estimating population extinction risk.     

Evaluating an icon of population persistence: the Devil's Hole pupfish

The Devil''s Hole pupfish Cyprinodon diabolis has iconic status among conservation biologists because it is one of the World''s most vulnerable species. Furthermore, C. diabolis is the most widely cited example of a persistent, small, isolated vertebrate population; a chronic exception to the rule that small populations do not persist long in isolation. It is widely asserted that this species has persisted in small numbers (less than 400 adults) for 10 000–20 000 years, but this assertion has never been evaluated. Here, we analyse the time series of count data for this species, and we estimate time to coalescence from microsatellite data to evaluate this hypothesis. We conclude that mean time to extinction is approximately 360–2900 years (median 410–1800), with less than a 2.1% probability of persisting 10 000 years. Median times to coalescence varied from 217 to 2530 years, but all five approximations had wide credible intervals. Our analyses suggest that Devil''s Hole pupfish colonized this pool well after the Pleistocene Lakes receded, probably within the last few hundred to few thousand years; this could have occurred through human intervention.     

Population size, weight distribution and food in a persistent population of the rare medicinal leech, Hirudo medicinalis

1. It is important for species recovery and conservation management projects to know the minimum viable population size for rare and endangered species, such as the medicinal leech, Hirudo medicinalis. Therefore, using a catch‐removal method, this study estimated every two years (1986, 1988, 1990, 1992) the total number of medicinal leeches in a tarn in the English Lake District, and the number of mature adults in the population. 2. Four samples were taken each year in June and July, when water temperatures exceeded 20 °C. Population size was estimated both by maximum likelihood and regression methods. All leeches were weighed alive and size groups were separated by polymodal frequency analysis. A small sample of the blood meal in each leech gut was taken before the leeches were returned to the tarn, and was used to estimate the proportion of mammalian and non‐mammalian blood in the meals. 3. Both methods of estimation produced similar values, increasing confidence in the population estimates. Values for the total population in June and July varied among years from 248 to 288, the maximum value being only 16% higher than the minimum. Values for the number of mature leeches varied from 48 to 58 (19–20% of the total population), and this was an estimate of the effective population size. 4. There were four size groups. The largest mature leeches (live weight >5 g) in group IV formed only 1% of the population, and the smallest (0.02–0.5 g) in group I 14–17%. Most leeches were in two overlapping groups of immature (64–67% of population) and mature (18%) leeches with size ranges of 0.4–3.4 g and 2.5–5 g respectively. The percentage of leeches in each size group was very consistent among years. Blood meals were found in 38–44% of the leeches in group I, 45–50% in group II, 70–75% in group III, and 100% in group IV, but mammalian blood was present only in larger mature leeches (>3.5 g). 5. Medicinal leeches were first detected in the tarn in 1980 and are still present in 2007, so the population has persisted for at least 27 years. Compared with minimum viable population sizes for other species, including many endangered species, values for this medicinal leech population are extremely low, but may be typical of some rare freshwater invertebrates in isolated habitats.     

The limits to population density in birds and mammals

We address two fundamental ecological questions: what are the limits to animal population density and what determines those limits? We develop simple alternative models to predict population limits in relation to body mass. A model assuming that within‐species area use increases with the square of daily travel distance broadly predicts the scaling of empirical extremes of minimum density across birds and mammals. Consistent with model predictions, the estimated density range for a given mass, ‘population scope’, is greater for birds than for mammals. However, unlike mammals and carnivorous birds, expected broad relationships between body mass and density extremes are not supported by data on herbivorous and omnivorous birds. Our results suggest that simple constraints on mobility and energy use/supply are major determinants of the scaling of density limits, but further understanding of interactions between dietary constraints and density limits are needed to predict future wildlife population responses to anthropogenic threats.     

Integrating genetic data and population viability analyses for the identification of harbour seal (Phoca vitulina) populations and management units

Identification of populations and management units is an essential step in the study of natural systems. Still, there is limited consensus regarding how to define populations and management units, and whether genetic methods allow for inference at the relevant spatial and temporal scale. Here, we present a novel approach, integrating genetic, life history and demographic data to identify populations and management units in southern Scandinavian harbour seals. First, 15 microsatellite markers and model‐ and distance‐based genetic clustering methods were used to determine the population genetic structure in harbour seals. Second, we used harbour seal demographic and life history data to conduct population viability analyses (PVAs) in the vortex simulation model in order to determine whether the inferred genetic units could be classified as management units according to Lowe and Allendorf's (Molecular Ecology, 19, 2010, 3038) ‘population viability criterion’ for demographic independence. The genetic analyses revealed fine‐scale population structuring in southern Scandinavian harbour seals and pointed to the existence of several genetic units. The PVAs indicated that the census population size of each of these genetic units was sufficiently large for long‐term population viability, and hence that the units could be classified as demographically independent management units. Our study suggests that population genetic inference can offer the same degree of temporal and spatial resolution as ‘nongenetic’ methods and that the combined use of genetic data and PVAs constitutes a promising approach for delineating populations and management units.     

Colour polymorphism is associated with lower extinction risk in birds

Colour polymorphisms have played a major role in enhancing current understanding of how selection and demography can impact phenotypes. Because different morphs often display alternative strategies and exploit alternative ecological niches, colour polymorphism can be expected to promote adaptability to environmental changes. However, whether and how it could influence populations' and species' response to global changes remains debated. To address this question, we built an up‐to‐date and complete database on avian colour polymorphism based on the examination of available data from all 10,394 extant bird species. We distinguished between true polymorphism (where different genetically determined morphs co‐occur in sympatry within the same population) and geographic variation (parapatric or allopatric colour variation), because these two patterns of variation are expected to have different consequences on populations' persistence. Using the IUCN red list, we then showed that polymorphic bird species are at lesser risk of extinction than nonpolymorphic ones, after controlling for a range of factors such as geographic range size, habitat breadth, life history, and phylogeny. This appears consistent with the idea that high genetic diversity and/or the existence of alternative strategies in polymorphic species promotes the ability to adaptively respond to changing environmental conditions. In contrast, polymorphic species were not less vulnerable than nonpolymorphic ones to specific drivers of extinction such as habitat alteration, direct exploitation, climate change, and invasive species. Thus, our results suggest that colour polymorphism acts as a buffer against environmental changes, although further studies are now needed to understand the underlying mechanisms. Developing accurate quantitative indices of sensitivity to specific threats is likely a key step towards a better understanding of species response to environmental changes.     

Using population viability analysis,genomics, and habitat suitability to forecast future population patterns of Little Owl Athene noctua across Europe

The agricultural scene has changed over the past decades, resulting in a declining population trend in many species. It is therefore important to determine the factors that the individual species depend on in order to understand their decline. The landscape changes have also resulted in habitat fragmentation, turning once continuous populations into metapopulations. It is thus increasingly important to estimate both the number of individuals it takes to create a genetically viable population and the population trend. Here, population viability analysis and habitat suitability modeling were used to estimate population viability and future prospects across Europe of the Little Owl Athene noctua, a widespread species associated with agricultural landscapes. The results show a high risk of population declines over the coming 100 years, especially toward the north of Europe, whereas populations toward the southeastern part of Europe have a greater probability of persistence. In order to be considered genetically viable, individual populations must count 1,000–30,000 individuals. As Little Owl populations of several countries count <30,000, and many isolated populations in northern Europe count <1,000 individuals, management actions resulting in exchange of individuals between populations or even countries are probably necessary to prevent losing <1% genetic diversity over a 100‐year period. At a continental scale, a habitat suitability analysis suggested Little Owl to be affected positively by increasing temperatures and urban areas, whereas an increased tree cover, an increasing annual rainfall, grassland, and sparsely vegetated areas affect the presence of the owl negatively. However, the low predictive power of the habitat suitability model suggests that habitat suitability might be better explained at a smaller scale.     

Density dependence and risk of extinction in a small population of sea otters

Sea otters (Enhydra lutris (L.)) were hunted to extinction off the coast of Washington State early in the 20th century. A new population was established by translocations from Alaska in 1969 and 1970. The population, currently numbering at least 550 animals, A major threat to the population is the ongoing risk of majour oil spills in sea otter habitat. We apply population models to census and demographic data in order to evaluate the status of the population. We fit several density dependent models to test for density dependence and determine plausible values for the carrying capacity (K) by comparing model goodness of fit to an exponential model. Model fits were compared using Akaike Information Criterion (AIC). A significant negative relationship was found between the population growth rate and population size (r ²=0.27, F=5.57, df=16, p<0.05), suggesting density dependence in Washington state sea otters. Information criterion statistics suggest that the model is the most parsimonious, followed closely by the logistic Beverton–Holt model. Values of K ranged from 612 to 759 with best-fit parameter estimates for the Beverton–Holt model including 0.26 for r and 612 for K. The latest (2001) population index count (555) puts the population at 87–92% of the estimated carrying capacity, above the suggested range for optimum sustainable population (OSP). Elasticity analysis was conducted to examine the effects of proportional changes in vital rates on the population growth rate (). The elasticity values indicate the population is most sensitive to changes in survival rates (particularly adult survival).     

Life-history characters and phylogeny are correlated with extinction risk in the Australian angiosperms

Aim To determine whether life‐history characters that affect population persistence (e.g. habit and life span) and those that influence reproductive success (e.g. sexual system and fruit type) are non‐randomly correlated with extinction risk (i.e. threat category) in the Australian flora (c. 19,000 species, of which c. 14% is threatened). To identify patterns that present useful conservation directions. To understand patterns of extinction risk in the Australian flora at a broad scale. Location Continental Australia. Methods A country‐wide exploration of four life‐history characters in the Australian flora (n = 18,822 species) was undertaken using reference texts, expert opinion, herbarium records and field work. For each character and threat‐category combination, a G‐test (using a log‐linear model) was performed to test the null hypothesis that the two factors were independent in their effects on count. A generalized linear model (GLM) with a logit link and binomial error distribution was constructed with the proportion of taxa in each extinction risk category as the response variable and the habit, sex and fruit‐type characters as explanatory terms. In a separate approach, we investigated patterns across the threat categories of non‐endangered extant, endangered, and extinct using a multinomial model. We examined whether or not species‐poor genera were more likely to contain threatened or extinct species than species‐rich genera. A GLM with a binomial error distribution and logit link function was constructed to obtain a weighted regression on the proportion of species listed as extinct or endangered within a genus versus the log of the size of the genus. We also used a supertree analysis and character tracing to investigate the role of phylogeny on extinction risk. Results We found that the Australian flora is primarily composed of bisexual shrubs with dry‐dehiscent fruits. Dioecious breeding systems (separate female and male flowers on separate plants) in many floras are the predominant unisexual system, but in Australia there are unexpectedly high levels of monoecy (separate female and male flowers on the same plant). Within the extinct data set of 31 species we detected a significant departure from that expected for habit but not for life span, sexual system or fruit type. There are significantly fewer trees on the extinct list than expected. This may reflect the greater resilience of trees than of other growth habits to extinction processes as well as the observation time‐frame. Within the endangered data set of 450 species we found significant differences in the representation of the observed characters from that expected within sex systems and fruit types. We show that, depending on the life form, unisexual breeding systems can be significantly and positively associated with endangered species compared with non‐threatened species. For example, there are more monoecious species than expected by chance among the tree species listed as endangered but fewer among the herbaceous life forms. Threat category was found to be non‐randomly clustered in some clades. Main conclusions Life‐history characters in certain combinations are predictive of extinction risk. Phylogeny is also an important component of extinction risk. We suggest that specific life‐history characters could be used for conservation planning and as an early warning sign for detecting vulnerability in lists of species.     

Environmental correlates of body size influence range size and extinction risk: A global study in rodents

Aim

Species-level traits, such as body and range sizes, are important correlates of extinction risk. However, both are often related and are driven by environmental factors. Here, we elucidated links between environmental factors, body size, range size and susceptibility to extinction, across the whole order of rodents.

Location

Global.

Time period

Current.

Major taxa studied

Rodents (order Rodentia).

Methods

We compiled an unprecedentedly large database of rodent morphology, phylogeny, range size, conservation status, global climate and normalized difference vegetation index (NDVI), comprising >86% of all described species. Using phylogenetic regressions, we initially explored the environmental factors driving body size. Next, we modelled the relationship between body size and range size. From this relationship, we computed and mapped (at the assemblage level) an index of relative range size, corresponding to the deviation from the expected range size of each species, given its body size. Finally, we tested whether relative range was correlated with the risk of extinction of the species derived from an assessment by the International Union for Conservation of Nature.

Results

We found that, contrary to the expectations of Bergmann's rule, the body size of rodents was mostly influenced by variation in NDVI (rather than latitude/temperature). Body size, in turn, imposed a constraint on species range size, as evidenced by a triangular relationship that was segmented at the lower bound. The relative species range size derived from this relationship highlighted four geographical regions where rodents with small relative range were concentrated globally. We demonstrated that lower relative range size was associated with increased risk of extinction.

Main conclusions

Species that, given their body size, are distributed across ranges that are smaller than expected have elevated extinction risk. Therefore, investigating the relationships between environmental drivers, body size and range size might help to detect species that could become threatened in the near future.     

Geographical variation in predictors of mammalian extinction risk: big is bad, but only in the tropics

Whereas previous studies have investigated correlates of extinction risk either at global or regional scales, our study explicitly models regional effects of anthropogenic threats and biological traits across the globe. Using phylogenetic comparative methods with a newly-updated supertree of 5020 extant mammals, we investigate the impact of species traits on extinction risk within each WWF ecoregion. Our analyses reveal strong geographical variation in the influence of traits on risk: notably, larger species are at higher risk only in tropical regions. We then relate these patterns to current and recent-historical human impacts across ecoregions using spatial modelling. The body–mass results apparently reflect historical declines of large species outside the tropics due to large-scale land conversion. Narrow-ranged and rare species tend to be at high risk in areas of high current human impacts. The interactions we describe between biological traits and anthropogenic threats increase understanding of the processes determining extinction risk.     

Environmental and human factors influencing rare plant local occurrence, extinction and persistence: a 115-year study in the Mediterranean region

Aim Assessing whether environmental and human factors influenced the spatial distribution and the dynamics of regionally rare plant species since the late nineteenth century, and whether these spatial and temporal patterns of rare species occurrences differ according to their chorology (level of endemism and biogeographic affinity). Location An area extending over 6250 km² in the French Mediterranean Region. Methods We used two botanical surveys achieved in 1886 and in 2001, and considered species rare if occurring in only one or two sites in the study area. Each rare species was assigned to a group of endemism level (restricted endemic, non‐endemic), and of biogeographic affinity (Mediterranean, South/Central European, Mountain, Eurosiberian). A 1 × 1 km grid was applied to the study zone. Generalized linear models were developed to study the spatial distribution and the fate of rare species occurrences (local extinction vs. local persistence between 1886 and 2001), as a function of environmental and human variables. Multivariate analyses were used to test whether the spatial distribution and the fate of rare species occurrences differed according to their chorology. Results In 2001, rare species as a whole tended to occur at higher altitude, in zones dominated by semi‐natural open habitats, and where cultivated area had decreased in the last 30 years. Between 1886 and 2001, rare species were the most prone to local extinction in zones where human population density, cultivated area and livestock density had increased the most. Between 1886 and 2001, rare species had a higher probability of local persistence in zones of high altitude and steep slope, on basic bedrocks and with low cultivated area. Rare species with Mountain and Eurosiberian affinities occurred in marginal habitats in the study region, i.e. on gneiss‐micaschist bedrocks and at high altitudes, whereas Mediterranean and South/Central European rare species occupied more varied environmental conditions. Between 1886 and 2001, Eurosiberian rare species showed high rates of local extinction whereas Mediterranean rare species had a significantly higher probability of local persistence. Restricted endemic species mostly occurred in zones of high slope, low human population density, and where cultivated area had decreased in the last 30 years. Occurrences of restricted endemics remained significantly stable between 1886 and 2001. Main conclusions Environmental and land‐use changes that occurred over the twentieth century in the Mediterranean Basin had significant impacts on the spatial distribution and on the long‐term dynamics of rare species occurrences. Urbanization and recent agriculture intensification, occurring mainly in coastal plains and littoral zones, caused most local extinctions of rare species from 1886 to 2001. Local populations of Eurosiberian species, which reach their range limits in marginal zones of the Mediterranean, also appear to be highly vulnerable. Conversely, most restricted endemic species occur in habitats with harsh topography and low human disturbance and have a higher potential of local persistence.     

Data cloning: easy maximum likelihood estimation for complex ecological models using Bayesian Markov chain Monte Carlo methods

We introduce a new statistical computing method, called data cloning, to calculate maximum likelihood estimates and their standard errors for complex ecological models. Although the method uses the Bayesian framework and exploits the computational simplicity of the Markov chain Monte Carlo (MCMC) algorithms, it provides valid frequentist inferences such as the maximum likelihood estimates and their standard errors. The inferences are completely invariant to the choice of the prior distributions and therefore avoid the inherent subjectivity of the Bayesian approach. The data cloning method is easily implemented using standard MCMC software. Data cloning is particularly useful for analysing ecological situations in which hierarchical statistical models, such as state-space models and mixed effects models, are appropriate. We illustrate the method by fitting two nonlinear population dynamics models to data in the presence of process and observation noise.