Critically evaluating best management practices for preventing freshwater turtle extinctions

Ex situ conservation tools, such as captive breeding for reintroduction, are considered a last resort to recover threatened or endangered species, but they may also help reduce anthropogenic threats where it is difficult or impossible to address them directly. Headstarting, or captive rearing of eggs or neonate animals for subsequent release into the wild, is controversial because it treats only a symptom of a larger conservation problem; however, it may provide a mechanism to address multiple threats, particularly near population centers. We conducted a population viability analysis of Australia's most widespread freshwater turtle, Chelodina longicollis, to determine the effect of adult roadkill (death by collision with motor vehicles), which is increasing, and reduced recruitment through nest predation from introduced European red foxes (Vulpes vulpes). We also modeled management scenarios to test the effectiveness of headstarting, fox management, and measures to reduce mortality on roads. Only scenarios with headstarting from source populations eliminated all risks of extinction and allowed population growth. Small increases in adult mortality (2%) had the greatest effect on population growth and extinction risk. Where threats simultaneously affected other life‐history stages (e.g., recruitment), eliminating harvest pressures on adult females alone did not eliminate the risk of population extinction. In our models, one source population could supply enough hatchlings annually to supplement 25 other similar‐sized populations such that extinction was avoided. Based on our results, we believe headstarting should be a primary tool for managing freshwater turtles for which threats affect multiple life‐history stages. We advocate the creation of source populations for managing freshwater turtles that are greatly threatened at multiple life‐history stages, such as depredation of eggs by invasive species and adult mortality via roadkill.     

Actual and Potential Use of Population Viability Analyses in Recovery of Plant Species Listed under the U.S. Endangered Species Act

Use of population viability analyses (PVAs) in endangered species recovery planning has been met with both support and criticism. Previous reviews promote use of PVA for setting scientifically based, measurable, and objective recovery criteria and recommend improvements to increase the framework's utility. However, others have questioned the value of PVA models for setting recovery criteria and assert that PVAs are more appropriate for understanding relative trade‐offs between alternative management actions. We reviewed 258 final recovery plans for 642 plants listed under the U.S. Endangered Species Act to determine the number of plans that used or recommended PVA in recovery planning. We also reviewed 223 publications that describe plant PVAs to assess how these models were designed and whether those designs reflected previous recommendations for improvement of PVAs. Twenty‐four percent of listed species had recovery plans that used or recommended PVA. In publications, the typical model was a matrix population model parameterized with ≤5 years of demographic data that did not consider stochasticity, genetics, density dependence, seed banks, vegetative reproduction, dormancy, threats, or management strategies. Population growth rates for different populations of the same species or for the same population at different points in time were often statistically different or varied by >10%. Therefore, PVAs parameterized with underlying vital rates that vary to this degree may not accurately predict recovery objectives across a species’ entire distribution or over longer time scales. We assert that PVA, although an important tool as part of an adaptive‐management program, can help to determine quantitative recovery criteria only if more long‐term data sets that capture spatiotemporal variability in vital rates become available. Lacking this, there is a strong need for viable and comprehensive methods for determining quantitative, science‐based recovery criteria for endangered species with minimal data availability. Uso Actual y Potencial del Análisis de Viabilidad Poblacional para la Recuperación de Especies de Plantas Enlistadas en el Acta de Especies En Peligro de E.U.A     

Population Viability Analysis with Species Occurrence Data from Museum Collections

Abstract: The most comprehensive data on many species come from scientific collections. Thus, we developed a method of population viability analysis (PVA) in which this type of occurrence data can be used. In contrast to classical PVA, our approach accounts for the inherent observation error in occurrence data and allows the estimation of the population parameters needed for viability analysis. We tested the sensitivity of the approach to spatial resolution of the data, length of the time series, sampling effort, and detection probability with simulated data and conducted PVAs for common, rare, and threatened species. We compared the results of these PVAs with results of standard method PVAs in which observation error is ignored. Our method provided realistic estimates of population growth terms and quasi‐extinction risk in cases in which the standard method without observation error could not. For low values of any of the sampling variables we tested, precision decreased, and in some cases biased estimates resulted. The results of our PVAs with the example species were consistent with information in the literature on these species. Our approach may facilitate PVA for a wide range of species of conservation concern for which demographic data are lacking but occurrence data are readily available.     

Demographic Viability of a Relict Population of the Critically Endangered Plant Borderea chouardii

Abstract:   In addition to human-caused changes in the environment, natural stochasticity may threaten species persistence, and its impact must be taken into account when priorities are established and management plans are designed. Borderea chouardii is a Tertiary relict at risk of extinction that occurs in only one location in the world, where the probability of human disturbance is low. Its persistence, therefore, is mainly linked to its response to natural threats such as stochasticity. Over 8 years I monitored up to 25% of this rupicolous small geophyte. The population had an unbalanced size structure and 90% failure in seed arrival at appropriate microhabitats, which suggests a problem with recruitment. I used matrix models to describe its population dynamics, conducted hand sowings, and performed stochastic simulations to investigate the effect of environmental stochasticity on population trend and viability. I modeled several scenarios to represent a variety of ecological situations, such as population reduction, episodic or persistent disease, and enhancement or decrease of recruitment. Population growth rate (λ) was never significantly different from unity over the study period. The risk of extinction was null over the next five centuries under current conditions. Increase of mortality and decrease of recruitment reduced stochastic population growth rate, but no factor except a persistent increase of 10% mortality resulted in extinction. These results are the consequence of the plant's extremely long life span (over 300 years) and low temporal variability of key vital rates. Even though hand sowing significantly increased the stochastic population growth rate, other approaches may be more important for the persistence of this species. The extremely slow capacity for recovery following disturbances renders habitat preservation essential. In addition, the founding of new populations would reduce the risk associated with habitat destruction.     

A critical appraisal of population viability analysis

Population viability analysis (PVA) is useful in management of imperiled species. Applications range from research design, threat assessment, and development of management frameworks. Given the importance of PVAs, it is essential that they be rigorous and adhere to widely accepted guidelines; however, the quality of published PVAs is rarely assessed. We evaluated the quality of 160 PVAs of 144 species of birds and mammals published in peer-reviewed journals from 1990 to 2017. We hypothesized that PVA quality would be lower with generic programs than with custom-built programs; be higher for those developed for imperiled species; change over time; and be higher for those published in journals with high impact factors (IFs). Each included study was evaluated based on answers to an evaluation framework containing 32 questions reflecting whether and to what extent the PVA study adhered to published PVA guidelines or contained important PVA components. All measures of PVA quality were generally lower for studies based on generic programs. Conservation status of the species did not affect any measure of PVA quality, but PVAs published in high IF journals were of higher quality. Quality generally declined over time, suggesting the quantitative literacy of PVA practitioners has not increased over time or that PVAs developed by unskilled users are being published in peer-reviewed journals. Only 18.1% of studies were of high quality (score >75%), which is troubling because poor-quality PVAs could misinform conservation decisions. We call for increased scrutiny of PVAs by journal editors and reviewers. Our evaluation framework can be used for this purpose. Because poor-quality PVAs continue to be published, we recommend caution while using PVA results in conservation decision making without thoroughly assessing the PVA quality.     

Factors Leading to Different Viability Predictions for a Grizzly Bear Data Set

Population viability analysis programs are being used increasingly in research and management applications, but there has not been a systematic study of the congruence of different program predictions based on a single data set. We performed such an analysis using four population viability analysis computer programs: GAPPS, INMAT, RAMAS/AGE, and VORTEX. The standardized demographic rates used in all programs were generalized from hypothetical increasing and decreasing grizzly bear ( Ursus arctos horribilis ) populations. Idiosyncracies of input format for each program led to minor differences in intrinsic growth rates that translated into striking differences in estimates of extinction rates and expected population size. In contrast, the addition of demographic stochasticity, environmental stochasticity, and inbreeding costs caused only a small divergence in viability predictions. But, the addition of density dependence caused large deviations between the programs despite our best attempts to use the same density-dependent functions. Population viability programs differ in how density dependence is incorporated, and the necessary functions are difficult to parameterize accurately. Thus, we recommend that unless data clearly suggest a particular density-dependent model, predictions based on population viability analysis should include at least one scenario without density dependence. Further, we describe output metrics that may differ between programs; development of future software could benefit from standardized input and output formats across different programs.     

Density dependence and density independence during the early life stages of four marine fish stocks

Recruitment variability caused by density-dependent and density-independent processes is an important area within the study of fish dynamics. These processes can exhibit nonlinearities and nonadditive properties that may have profound dynamic effects. We investigate the importance of population density (i.e., density dependence) and environmental forcing (i.e., density independence) on the age-0 and age-1 abundance of capelin (Mallotus villosus), northeast Arctic cod (Gadus morhua), northeast Arctic haddock (Melanogrammus aeglefinus), and Norwegian spring spawning herring (Clupea harengus) in the Barents Sea. We use statistical methods that explicitly account for nonlinearities and nonadditive interactions between internal and external variables in the abundance of these two pre-recruitment stages. Our results indicate that, during their first five months of life, cod, haddock, and herring experience higher density-dependent survival than capelin. The abundance of age-0 cod depends on the mean age and biomass of the spawning stock, a result which has implications for the management of the entire cod stock. Temperature is another important factor influencing the abundance at age-0 and age-1 of all four species, except herring at age-1. Between age-0 and age-1, there is an attenuation of density-dependent survival for cod and herring, while haddock and capelin experience density dependence at high and low temperatures, respectively. Predation by subadult cod is important for both capelin and cod at age-1. We found strong indications for interactions among the studied species, pointing to the importance of viewing the problem of species recruitment variability as a community, rather than as a population phenomenon.     

Effects of Human Activity on Global Extinction Risk

Both natural and anthropogenic factors are important in determining a species' risk of extinction. Little work has been done, however, to quantify the magnitude of current anthropogenic influences on the extinction process. The purpose of this study is to determine the extent to which measures of the intensity of human activity are related to the global variability of two measures of species' susceptibility to extinction. We observed six indices of human activities in 90 countries, and we tested their relationships to the proportion of threatened bird and mammal species in each country, as well as to mammalian population density. After correcting for area effects, latitudinal diversity gradients, and body size (for population density), 28 to 50% of the remaining variation was statistically attributable to anthropogenic variables. Different measures of anthropogenic influence were most closely related to extinction risk in birds and mammals. Human population density was the variable most closely related to the proportion of threatened bird species per country, whereas per capita GNP was more important for mammal species. Mammalian population density strongly correlates with the extent of protected area per country. Contrary to suggestions in earlier literature, our work does not support the hypothesis that habitat loss is a prime contributor to species loss because frequencies of threatened birds and mammals are not closely related to patterns of land use.     

The Reliability of Using Population Viability Analysis for Risk Classification of Species

I examine whether or not it is appropriate to use extinction probabilities generated by population viability analyses, based on best estimates for model parameters, as criteria for listing species in Red Data Book categories as recently proposed by the World Conservation Union. Such extinction probabilities are influenced by how accurately model parameters are estimated and by how accurately the models depict actual population dynamics. I evaluate the effect of uncertainty in parameter estimation through simulations. Simulations based on Steller sea lions were used to evaluate bias and precision in estimates of probability of extinction and to consider the performance of two proposed classification schemes. Extinction time estimates were biased (because of violation of the assumption of stable age distribution) and underestimated the variability of probability of extinction for a given time (primarily because of uncertainty in parameter estimation). Bias and precision in extinction probabilities are important when these probabilities are used to compare the risk of extinction between species. Suggestions are given for population viability analysis techniques that incorporate parameter uncertainty. I conclude that testing classification schemes with simulations using quantitative performance objectives should precede adoption of quantitative listing criteria.     

Incorporating Climate and Ocean Change into Extinction Risk Assessments for 82 Coral Species

Many marine invertebrate species facing potential extinction have uncertain taxonomies and poorly known demographic and ecological traits. Uncertainties are compounded when potential extinction drivers are climate and ocean changes whose effects on even widespread and abundant species are only partially understood. The U.S. Endangered Species Act mandates conservation management decisions founded on the extinction risk to species based on the best available science at the time of consideration—requiring prompt action rather than awaiting better information. We developed an expert‐opinion threat‐based approach that entails a structured voting system to assess extinction risk from climate and ocean changes and other threats to 82 coral species for which population status and threat response information was limited. Such methods are urgently needed because constrained budgets and manpower will continue to hinder the availability of desired data for many potentially vulnerable marine species. Significant species‐specific information gaps and uncertainties precluded quantitative assessments of habitat loss or population declines and necessitated increased reliance on demographic characteristics and threat vulnerabilities at genus or family levels. Adapting some methods (e.g., a structured voting system) used during other assessments and developing some new approaches (e.g., integrated assessment of threats and demographic characteristics), we rated the importance of threats contributing to coral extinction risk and assessed those threats against population status and trend information to evaluate each species’ extinction risk over the 21st century. This qualitative assessment resulted in a ranking with an uncertainty range for each species according to their estimated likelihood of extinction. We offer guidance on approaches for future biological extinction risk assessments, especially in cases of data‐limited species likely to be affected by global‐scale threats. Incorporación del Cambio Climático y Oceánico en Estudios de Riesgo de Extinción para 82 Especies de Coral     

Effects of Social Structure and Prey Dynamics on Extinction Risk in Gray Wolves

Extinction models based on diffusion theory generally fail to incorporate two important aspects of population biology—social structure and prey dynamics. We include these aspects in an individual-based extinction model for small, isolated populations of the gray wolf (Canis lupus). Our model predicts mean times to extinction significantly longer than those predicted by more general (diffusion) models. According to our model, an isolated population of 50 wolves has a 95% chance of surviving just 9 years and only a 30% chance of surviving beyond 100 years. Reflecting the influence of social structure, a wolf population initially comprising 50 individuals is expected to persist only a few years longer, on average (71 years), than is a population initially comprising just a single reproductive pair (62 years). In contrast, substantially greater average prey abundance leads to dramatically longer expected persistence times. Autocorrelated prey dynamics result in a more complex distribution of extinction times than predicted by many extinction models. We contend that demographic stochasticity may pose the greatest threat to small, isolated wolf populations, although environmental stochasticity and genetic effects may compound this threat. Our work highlights the importance of considering social structure and resource dynamics in the development of population viability analyses.     

Population ecology of a freshwater turtle Kachuga tentoria near Panchnada (Etawah: U.P.) and its role as water purifier

The present study deals with biology, ecology and population dynamics of freshwater turtle Kachuga tentoria and its role as water purifier. The study area Panchnada is the site, where five important national rivers meet together and is preserving an appreciable population of nine species of fresh water turtles. Kachuga tentoria was located at all the sampling stations surveyed by the authors, and hence selected for the present study. Different activities (nesting, incubation, predation and other reproductive aspects), climatic conditions, habitat, population density and morphometric features were worked out in detail. A time bound conservation strategy is needed to save this species from extinction. In situ conservation will be more helpful for the recruitment of the population of this species. It will help in the "hatch and release programme" to clean different polluted national rivers.     

Sensitivity of population viability to spatial and nonspatial parameters using GRIP.

Metapopulation dynamics are influenced by spatial parameters including the amount and arrangement of suitable habitat, yet these parameters may be uncertain when deciding how to manage species or their habitats. Sensitivity analyses of population viability analysis (PVA) models can help measure relative parameter influences on predictions, identify research priorities for reducing uncertainty, and evaluate management strategies. Few spatial PVAs, however, include sensitivity analyses of both spatial and nonspatial parameters, perhaps because computationally efficient tools for such analyses are lacking or inaccessible. We developed GRIP, a program to facilitate sensitivity analysis of spatial and nonspatial input parameters for PVAs created in RAMAS Metapop, a widely applied software program. GRIP creates random sets of input files by varying parameters specified in the PVA model including vital rates and their correlations among populations, the number and configuration of populations, dispersal rates, dispersal survival, initial population abundances, carrying capacities, and the probability, intensity, and spatial extent of catastrophes, while drawing on specified parameter distributions. We evaluated GRIP's performance as a tool for sensitivity analysis of spatial PVAs and explored the consequences of varying spatial input parameters for predictions of a published PVA model of the sand lizard (Lacerta agilis). We used GRIP output to generate standardized regression coefficients (SRCs) and nonparametric correlation coefficients as indices of the relative sensitivity of predicted conservation status to input parameters. GRIP performed well; with a single analysis we were able to rank the relative influence of input parameters identified as influential by the PVA's original author, S. A. Berglind, who used three separate forms of sensitivity analysis. Our analysis, however, also underscored the value of exploring the relative influence of spatial parameters on PVA predictions; both SRCs and correlation coefficients indicated that the most influential parameters in the sand lizard model were spatial in nature. We provide annotated code so that GRIP may be modified to reflect particular species biology, customized for more complex spatial PVA models, upgraded to incorporate features added in newer versions of RAMAS Metapop, used as a template to develop similar programs, or used as it is for computationally efficient sensitivity analyses in support of conservation planning.     

Habitat Change and Plant Demography: Assessing the Extinction Risk of a Formerly Common Grassland Perennial

Abstract:  An important aim of conservation biology is to understand how habitat change affects the dynamics and extinction risk of populations. We used matrix models to analyze the effect of habitat degradation on the demography of the declining perennial plant Trifolium montanum in 9 calcareous grasslands in Germany over 4 years and experimentally tested the effect of grassland management. Finite population growth rates (λ) decreased with light competition, measured as leaf-area index above T. montanum plants. At unmanaged sites λ was <1 due to lower recruitment and lower survival and flowering probability of large plants. Nevertheless, in stochastic simulations, extinction of unmanaged populations of 100 flowering plants was delayed for several decades. Clipping as a management technique rapidly increased population growth because of higher survival and flowering probability of large plants in managed than in unmanaged plots. Transition-matrix simulations from these plots indicated grazing or mowing every second year would be sufficient to ensure a growth rate ≥1 if conditions stayed the same. At frequently grazed sites, the finite growth rate was approximately 1 in most populations of T. montanum . In stochastic simulations, the extinction risk of even relatively small grazed populations was low, but about half the extant populations of T. montanum in central Germany are smaller than would be sufficient for a probability of survival of >95% over 100 years. We conclude that habitat change after cessation of management strongly reduces recruitment and survival of established individuals of this perennial plant. Nevertheless, our results suggest extinction processes may take a long time in perennial plants, resulting in an extinction debt. Even if management is frequent, many remnant populations of T. montanum may be at risk because of their small size, but even a slight increase in size could considerably reduce their extinction risk.     

Demography and Population Viability of an Endangered Plant Species before and after Protection from Trampling

Land managers often suggest fencing to protect rare plant species from being trampled in heavily used recreation areas, but there are few documented examples of the efficacy of this strategy. In a 7-year demographic study we examined the reproduction, survival, and long-term viability of the endangered sentry milk-vetch (Astragalus cremnophylax var. cremnophylax) before and after protection from trampling. Demographic monitoring and population viability analyses indicated that the population has fluctuated during the 7 years. Before protection the population declined: 26% of individuals died, mortality surpassed natality, and age of first reproduction was significantly older than post-protection. Fifty-eight percent of the population was severely damaged. Population viability analyses of pre-protection years predicted that the population would go extinct within 100 years. Since protection, the population stabilized, grew, and declined again. Seedlings reached reproductive maturity more quickly. Recruitment increased and peaked in 1993 coincident with abundant precipitation, but again declined in 1994. The total numbers of undamaged plants surpassed the numbers of damaged plants. Models of the post-protection population predict stability. Multiple-linear regression analysis indicated that winter and spring precipitation were significantly correlated with lambda. Both "good" and "bad" climatic conditions occurred during the pre- and post-protection periods. Because of small population size and depauperate genetic diversity, climate will continue to influence population growth. Nevertheless, models indicate that where trampling and bad climatic conditions were coupled, extinction was accelerated. Recovery of sentry milk-vetch will depend on continued protection, augmentation, and environmental factors, although risk of extinction remains very high.     

Influence of variability in larval development on recruitment success in the euphausiid Euphausia pacifica : elasticity and sensitivity analyses

Simple life-history models were employed to identify those larval stages of Euphausia pacifica (Hansen) with the greatest potential contribution to variability in recruitment success. Using stage-specific vital parameters (i.e. stage duration and specific mortality) the models explored the influence of developmental variability without ascribing the variability to any particular cause or forcing. Through such a non-mechanistic modelling approach, the critical periods in larval development were determined. The results led to inferences about larval population responses to temperature/food heterogeneity in the natural environment. Peaks in recruitment-success elasticity were evident for variability imposed during the egg and Furcilia I and II stages. The high elasticity of the egg stage suggested the importance of female nutrition and the timing and location of egg release in determining their viability and subsequent recruitment success. Peak elasticity during the Furcilia I and II stages supported the notion that indirect developmental pathways of pleopod development during these stages may be associated with diminished recruitment success. Sensitivity analyses further highlighted the influence of these early furcilia stages in determining recruitment success, and also suggested the importance of the Furcilia III and VII stages. The robustness of elasticity and sensitivity results with respect to the life-history structure utilized is discussed. A criterion is presented for the identification of life-history structures which might obscure the true elasticity relationships of a developmental sequence. Received: 27 April 1998 / Accepted: 14 September 1998     

The scale transition: scaling up population dynamics with field data

Applying the recent developments of scale transition theory, we demonstrate a systematic approach to the problem of scaling up local scale interactions to regional scale dynamics with field data. Dynamics on larger spatial scales differ from the predictions of local dynamics alone because of an interaction between nonlinearity in population dynamics at the local scale and spatial variation in density and environmental factors over the regional population. Our systematic approach to scaling up involves the following five steps. First, define a model for dynamics on the local spatial scale. Second, apply scale transition theory to identify key interactions between nonlinearity and spatial variation that translate local dynamics to the regional scale. Third, measure local-scale model parameters to determine nonlinearities at local scales. Fourth, measure spatial variation. Finally, combine nonlinearity and variation measures to obtain the scale transition. Using field data for the dynamics of grazers and periphyton in a freshwater stream, we show that scale transition terms greatly reduce the growth and equilibrium density of the periphyton population at the stream scale compared to rock scale populations, confirming the importance of spatial mechanisms to stream-scale dynamics.     

Reconstructibility of Density Dependence and the Conservative Assessment of Extinction Risks

Abstract: The probability of extinction is sensitive to the presence and character of density dependence controlling the dynamics of a population. This means that our capacity to estimate a population's risks of extinction under varying environmental conditions or competing management regimes is linked to our ability to reconstruct from data the density-dependence relationships governing the natural dynamics, especially when data do not reveal a trend of population growth or decline. In an example using Gadus morhua , we show that even 10- or 20-year data sets are too short to make precise estimates of these risks. We also observe, however, that under moderate or weak density dependence, the computed risks are lower than when density dependence is not included in the model. We propose, therefore, that when available data sets are insufficient for reconstructing reliable measurements of density dependence, conservative estimates of extinction probabilities can be made from models that simply omit density dependence.     

Hidden demographic impacts of fishing and environmental drivers of fecundity in a sea turtle population

Fisheries bycatch is a critical threat to sea turtle populations worldwide, particularly because turtles are vulnerable to multiple gear types. The Canary Current is an intensely fished region, yet there has been no demographic assessment integrating bycatch and population management information of the globally significant Cabo Verde loggerhead turtle (Caretta caretta) population. Using Boa Vista island (Eastern Cabo Verde) subpopulation data from capture–recapture and nest monitoring (2013–2019), we evaluated population viability and estimated regional bycatch rates (2016–2020) in longline, trawl, purse-seine, and artisanal fisheries. We further evaluated current nesting trends in the context of bycatch estimates, existing hatchery conservation measures, and environmental (net primary productivity) variability in turtle foraging grounds. We projected that current bycatch mortality rates would lead to the near extinction of the Boa Vista subpopulation. Bycatch reduction in longline fisheries and all fisheries combined would increase finite population growth rate by 1.76% and 1.95%, respectively. Hatchery conservation increased hatchling production and reduced extinction risk, but alone it could not achieve population growth. Short-term increases in nest counts (2013–2021), putatively driven by temporary increases in net primary productivity, may be masking ongoing long-term population declines. When fecundity was linked to net primary productivity, our hindcast models simultaneously predicted these opposing long-term and short-term trends. Consequently, our results showed conservation management must diversify from land-based management. The masking effect we found has broad-reaching implications for monitoring sea turtle populations worldwide, demonstrating the importance of directly estimating adult survival and that nest counts might inadequately reflect underlying population trends.