A critical appraisal of population viability analysis

Short-term extinction predicted by population viability analysis for a Neotropical salt marsh endemic bird

Salt marshes pose challenges for the birds that inhabit them, including high rates of nest flooding, tipping, and predation. The impacts of rising sea levels and invasive species further exacerbate these challenges. To assess the urgency of conservation and adequacy of new actions, researchers and wildlife managers may use population viability analyses (PVAs) to identify population trends and major threats. We conducted PVA for Formicivora acutirostris, which is a threatened neotropical bird species endemic to salt marshes. We studied the species' demography in different sectors of an estuary in southern Brazil from 2006 to 2023 and estimated the sex ratio, longevity, productivity, first-year survival, and mortality rates. For a 133-year period, starting in 1990, we modeled four scenarios: (1) pessimistic and (2) optimistic scenarios, including the worst and best values for the parameters; (3) a baseline scenario, with intermediate values; and (4) scenarios under conservation management, with increased recruitment and/or habitat preservation. Projections indicated population decline for all assessment scenarios, with a 100% probability of extinction by 2054 in the pessimistic scenario and no extinction in the optimistic scenario. The conservation scenarios indicated population stability with 16% improvement in productivity, 10% improvement in first-year survival, and stable carrying capacity. The disjunct distribution of the species, with remnants concentrated in a broad interface with arboreal habitats, may seal the population decline by increasing nest predation. The species should be considered conservation dependent, and we recommend assisted colonization, predator control, habitat recovery, and ex situ conservation.

How admixed captive breeding populations could be rescued using local ancestry information

This paper asks the question: can genomic information be used to recover a species that is already on the pathway to extinction due to genetic swamping from a related and more numerous population? We show that a breeding strategy in a captive breeding program can use whole genome sequencing to identify and remove segments of DNA introgressed through hybridisation. The proposed policy uses a generalized measure of kinship or heterozygosity accounting for local ancestry, that is, whether a specific genetic location was inherited from the target of conservation. We then show that optimizing these measures would minimize undesired ancestry while also controlling kinship and/or heterozygosity, in a simulated breeding population. The process is applied to real data representing the hybridized Scottish wildcat breeding population, with the result that it should be possible to breed out domestic cat ancestry. The ability to reverse introgression is a powerful tool brought about through the combination of sequencing with computational advances in ancestry estimation. Since it works best when applied early in the process, important decisions need to be made about which genetically distinct populations should benefit from it and which should be left to reform into a single population.

Inferring stability and persistence in the vaginal microbiome: A stochastic model of ecological dynamics

The interplay of stochastic and ecological processes that govern the establishment and persistence of host-associated microbial communities is not well understood. Here we illustrate the conceptual and practical advantages of fitting stochastic population dynamics models to multi-species bacterial time series data. We show how the stability properties, fluctuation regimes and persistence probabilities of human vaginal microbial communities can be better understood by explicitly accommodating three sources of variability in ecological stochastic models of multi-species abundances: 1) stochastic biotic and abiotic forces, 2) ecological feedback and 3) sampling error. Rooting our modeling tool in stochastic population dynamics modeling theory was key to apply standardized measures of a community's reaction to environmental variation that ultimately depends on the nature and intensity of the intra-specific and inter-specific interaction strengths. Using estimates of model parameters, we developed a Risk Prediction Monitoring (RPM) tool that estimates temporal changes in persistence probabilities for any bacterial group of interest. This method mirrors approaches that are often used in conservation biology in which a measure of extinction risks is periodically updated with any change in a population or community. Additionally, we show how to use estimates of interaction strengths and persistence probabilities to formulate hypotheses regarding the molecular mechanisms and genetic composition that underpin different types of interactions. Instead of seeking a definition of "dysbiosis" we propose to translate concepts of theoretical ecology and conservation biology methods into practical approaches for the management of human-associated bacterial communities.

Situating defaunation in an operational framework to advance biodiversity conservation

Anthropogenic pressures are causing the widespread loss of wildlife species and populations, with adverse consequences for ecosystem functioning. This phenomenon has been widely but inconsistently referred to as defaunation. A cohesive, quantitative framework for defining and evaluating defaunation is necessary for advancing biodiversity conservation. Likening defaunation to deforestation, we propose an operational framework for defaunation that defines it and related terms, situates defaunation relative to intact communities and faunal degradation, and encourages quantitative, ecologically reasonable, and equitable measurements. We distinguish between defaunation, the conversion of an ecosystem from having wild animals to not having wild animals, and faunal degradation, the process of losing animals or species from an animal community. The quantification of context-relevant defaunation boundaries or baselines is necessary to compare faunal communities over space and time. Situating a faunal community on the degradation curve can promote Global Biodiversity Framework targets, advancing the 2050 Vision for Biodiversity.

Negative density-dependence buffers against mismatch-induced population decline in the Sinai baton blue butterfly

Phenological mismatches caused by climate change pose a major threat to global biodiversity, yet relatively few studies have reported population declines resulting from mismatch. It has been hypothesised that density effects may underlie this lack of observed responses by buffering against mismatch-induced population decline. We developed an individual-based model of the critically endangered Sinai baton blue butterfly (Pseudophilotes sinaicus) and its hostplant Sinai thyme (Thymus decussatus), parameterised using real field data, to test this hypothesis. Our model showed that the baton blue experiences demographic consequences under only 5 days of phenological mismatch, but that this threshold was increased to 14 days with the inclusion of density-dependent juvenile mortality. The inclusion of density effects also led to the replication of population cycles observed in nature, supporting the ability of our model to accurately represent the baton blue's ecology. These results add to a growing body of literature suggesting that density effects may underlie the observed lack of demographic responses to mismatch in wild populations. However, these buffers may be short-lived in extreme mismatch scenarios, providing a false sense of security against a looming threat of population collapse.

Return to 1616: Multispecies Fauna Reconstruction Requires Thinking Outside the Box

Conservation translocations have become increasingly popular for 'rewilding' areas that have lost their native fauna. These multispecies translocations are complex and need to consider the requirements of each individual species as well as the influence of likely interactions among them. The Dirk Hartog Island National Park Ecological Restoration Project, Return to 1616, aspires to restore ecological function to Western Australia's largest island. Since 2012, pest animals have been eradicated, and conservation translocations of seven fauna species have been undertaken, with a further six planned. Here, we present a synthesis of the innovative approaches undertaken in restoring the former faunal assemblage of Dirk Hartog Island and the key learnings gathered as the project has progressed.

Estimating Wolf Population Size and Dynamics by Field Monitoring and Demographic Models: Implications for Management and Conservation

We estimated the current size and dynamics of the wolf population in Tuscany and investigated the trends and demographic drivers of population changes. Estimates were obtained by two different approaches: (i) mixed-technique field monitoring (from 2014 to 2016) that found the minimum observed pack number and estimated population size, and (ii) an individual-based model (run by Vortex software v. 10.3.8.0) with demographic inputs derived from a local intensive study area and historic data on population size. Field monitoring showed a minimum population size of 558 wolves (SE = 12.005) in 2016, with a density of 2.74 individuals/100 km2. The population model described an increasing trend with an average annual rate of increase λ = 1.075 (SE = 0.014), an estimated population size of about 882 individuals (SE = 9.397) in 2016, and a density of 4.29 wolves/100 km2. Previously published estimates of wolf population were as low as 56.2% compared to our field monitoring estimation and 34.6% in comparison to our model estimation. We conducted sensitivity tests to analyze the key parameters driving population changes based on juvenile and adult mortality rates, female breeding success, and litter size. Mortality rates played a major role in determining intrinsic growth rate changes, with adult mortality accounting for 62.5% of the total variance explained by the four parameters. Juvenile mortality was responsible for 35.8% of the variance, while female breeding success and litter size had weak or negligible effects. We concluded that reliable estimates of population abundance and a deeper understanding of the role of different demographic parameters in determining population dynamics are crucial to define and carry out appropriate conservation and management strategies to address human-wildlife conflicts.

A comparison of population viability measures

The viability of populations can be quantified with several measures, such as the probability of extinction, the mean time to extinction, or the population size. While conservation management decisions can be based on these measures, it has not yet been explored systematically if different viability measures rank species and scenarios similarly and if one viability measure can be converted into another to compare studies. To address this challenge, we conducted a quantitative comparison of eight viability measures based on the simulated population dynamics of more than 4500 virtual species. We compared (a) the ranking of scenarios based on different viability measures, (b) assessed direct correlations between the measures, and (c) explored if parameters in the simulation models can alter the relationship between pairs of viability measures. We found that viability measures ranked species similarly. Despite this, direct correlations between the different measures were often weak and could not be generalized. This can be explained by the loss of information due to the aggregation of raw data into a single number, the effect of model parameters on the relationship between viability measures, and because distributions, such as the probability of extinction over time, cannot be ranked objectively. Similar scenario rankings by different viability measures show that the choice of the viability metric does in many cases not alter which population is regarded more viable or which management option is the best. However, the more two scenarios or populations differ, the more likely it becomes that different measures produce different rankings. We thus recommend that PVA studies publish raw simulation data, which not only describes all risks and opportunities to the reader but also facilitates meta-analyses of PVA studies.

The abundance and persistence of Caprinae populations

Stable or growing populations may go extinct when their sizes cannot withstand large swings in temporal variation and stochastic forces. Hence, the minimum abundance threshold defining when populations can persist without human intervention forms a key conservation parameter. We identify this threshold for many populations of Caprinae, typically threatened species lacking demographic data. Doing so helps triage conservation and management actions for threatened or harvested populations. Methodologically, we used population projection matrices and simulations, with starting abundance, recruitment, and adult female survival predicting future abundance, growth rate (λ), and population trend. We incorporated mean demographic rates representative of Caprinae populations and corresponding variances from desert bighorn sheep (Ovis canadensis nelsoni), as a proxy for Caprinae sharing similar life histories. We found a population’s minimum abundance resulting in ≤ 0.01 chance of quasi-extinction (QE; population ≤ 5 adult females) in 10 years and ≤ 0.10 QE in 30 years as 50 adult females, or 70 were translocation (removals) pursued. Discovering the threshold required 3 demographic parameters. We show, however, that monitoring populations’ relationships to this threshold requires only abundance and recruitment data. This applied approach avoids the logistical and cost hurdles in measuring female survival, making assays of population persistence more practical.

Life-history predicts global population responses to the weather in terrestrial mammals

With the looming threat of abrupt ecological disruption due to a changing climate, predicting which species are most vulnerable to environmental change is critical. The life-history of a species is an evolved response to its environmental context, and therefore a promising candidate for explaining differences in climate-change responses. However, we need broad empirical assessments from across the world's ecosystems to explore the link between life history and climate-change responses. Here, we use long-term abundance records from 157 species of terrestrial mammals and a two-step Bayesian meta-regression framework to investigate the link between annual weather anomalies, population growth rates, and species-level life history. Overall, we found no directional effect of temperature or precipitation anomalies or variance on annual population growth rates. Furthermore, population responses to weather anomalies were not predicted by phylogenetic covariance, and instead there was more variability in weather responses for populations within a species. Crucially, however, long-lived mammals with smaller litter sizes had smaller absolute population responses to weather anomalies compared with their shorter living counterparts with larger litters. These results highlight the role of species-level life history in driving responses to the environment.

Female reproductive skew exacerbates the extinction risk from poaching in the eastern black rhino

Variation in individual demographic rates can have large consequences for populations. Female reproductive skew is an example of structured demographic heterogeneity where females have intrinsic qualities that make them more or less likely to breed. The consequences of reproductive skew for population dynamics are poorly understood in non-cooperatively breeding mammals, especially when coupled with other drivers such as poaching. We address this knowledge gap with population viability analyses using an age-specific, female-only, individual-based, stochastic population model built with long-term data for three Kenyan populations of the Critically Endangered eastern black rhino (Diceros bicornis michaeli). There was substantial reproductive skew, with a high proportion of females not breeding or doing so at very low rates. This had a large impact on the projected population growth rate for the smaller population on Ol Jogi. Moreover, including female reproductive skew exacerbates the effects of poaching, increasing the probability of extinction by approximately 70% under a simulated poaching pressure of 5% offtake per year. Tackling the effects of reproductive skew depends on whether it is mediated by habitat or social factors, with potential strategies including habitat and biological management respectively. Investigating and tackling reproductive skew in other species requires long-term, individual-level data collection.