A case study of bats and white‐nose syndrome demonstrating how to model population viability with evolutionary effects

Ecological factors generally affect population viability on rapid time scales. Traditional population viability analyses (PVA) therefore focus on alleviating ecological pressures, discounting potential evolutionary impacts on individual phenotypes. Recent studies of evolutionary rescue (ER) focus on cases in which severe, environmentally induced population bottlenecks trigger a rapid evolutionary response that can potentially reverse demographic threats. ER models have focused on shifting genetics and resulting population recovery, but no one has explored how to incorporate those findings into PVA. We integrated ER into PVA to identify the critical decision interval for evolutionary rescue (DIER) under which targeted conservation action should be applied to buffer populations undergoing ER against extinction from stochastic events and to determine the most appropriate vital rate to target to promote population recovery. We applied this model to little brown bats (Myotis lucifugus) affected by white‐nose syndrome (WNS), a fungal disease causing massive declines in several North American bat populations. Under the ER scenario, the model predicted that the DIER period for little brown bats was within 11 years of initial WNS emergence, after which they stabilized at a positive growth rate (λ = 1.05). By comparing our model results with population trajectories of multiple infected hibernacula across the WNS range, we concluded that ER is a potential explanation of observed little brown bat population trajectories across multiple hibernacula within the affected range. Our approach provides a tool that can be used by all managers to provide testable hypotheses regarding the occurrence of ER in declining populations, suggest empirical studies to better parameterize the population genetics and conservation‐relevant vital rates, and identify the DIER period during which management strategies will be most effective for species conservation.     

Using global sensitivity analysis of demographic models for ecological impact assessment

Population viability analysis (PVA) is widely used to assess population‐level impacts of environmental changes on species. When combined with sensitivity analysis, PVA yields insights into the effects of parameter and model structure uncertainty. This helps researchers prioritize efforts for further data collection so that model improvements are efficient and helps managers prioritize conservation and management actions. Usually, sensitivity is analyzed by varying one input parameter at a time and observing the influence that variation has over model outcomes. This approach does not account for interactions among parameters. Global sensitivity analysis (GSA) overcomes this limitation by varying several model inputs simultaneously. Then, regression techniques allow measuring the importance of input‐parameter uncertainties. In many conservation applications, the goal of demographic modeling is to assess how different scenarios of impact or management cause changes in a population. This is challenging because the uncertainty of input‐parameter values can be confounded with the effect of impacts and management actions. We developed a GSA method that separates model outcome uncertainty resulting from parameter uncertainty from that resulting from projected ecological impacts or simulated management actions, effectively separating the 2 main questions that sensitivity analysis asks. We applied this method to assess the effects of predicted sea‐level rise on Snowy Plover (Charadrius nivosus). A relatively small number of replicate models (approximately 100) resulted in consistent measures of variable importance when not trying to separate the effects of ecological impacts from parameter uncertainty. However, many more replicate models (approximately 500) were required to separate these effects. These differences are important to consider when using demographic models to estimate ecological impacts of management actions.     

Determining the drivers of population structure in a highly urbanized landscape to inform conservation planning

Understanding the environmental contributors to population structure is of paramount importance for conservation in urbanized environments. We used spatially explicit models to determine genetic population structure under current and future environmental conditions across a highly fragmented, human‐dominated environment in Southern California to assess the effects of natural ecological variation and urbanization. We focused on 7 common species with diverse habitat requirements, home‐range sizes, and dispersal abilities. We quantified the relative roles of potential barriers, including natural environmental characteristics and an anthropogenic barrier created by a major highway, in shaping genetic variation. The ability to predict genetic variation in our models differed among species: 11–81% of intraspecific genetic variation was explained by environmental variables. Although an anthropogenically induced barrier (a major highway) severely restricted gene flow and movement at broad scales for some species, genetic variation seemed to be primarily driven by natural environmental heterogeneity at a local level. Our results show how assessing environmentally associated variation for multiple species under current and future climate conditions can help identify priority regions for maximizing population persistence under environmental change in urbanized regions.     

Ability of Matrix Models to Explain the Past and Predict the Future of Plant Populations

Uncertainty associated with ecological forecasts has long been recognized, but forecast accuracy is rarely quantified. We evaluated how well data on 82 populations of 20 species of plants spanning 3 continents explained and predicted plant population dynamics. We parameterized stage‐based matrix models with demographic data from individually marked plants and determined how well these models forecast population sizes observed at least 5 years into the future. Simple demographic models forecasted population dynamics poorly; only 40% of observed population sizes fell within our forecasts’ 95% confidence limits. However, these models explained population dynamics during the years in which data were collected; observed changes in population size during the data‐collection period were strongly positively correlated with population growth rate. Thus, these models are at least a sound way to quantify population status. Poor forecasts were not associated with the number of individual plants or years of data. We tested whether vital rates were density dependent and found both positive and negative density dependence. However, density dependence was not associated with forecast error. Forecast error was significantly associated with environmental differences between the data collection and forecast periods. To forecast population fates, more detailed models, such as those that project how environments are likely to change and how these changes will affect population dynamics, may be needed. Such detailed models are not always feasible. Thus, it may be wiser to make risk‐averse decisions than to expect precise forecasts from models. Habilidad de los Modelos Matriciales para Explicar el Pasado y Predecir el Futuro de las Poblaciones de Plantas     

Interacting Effects of Phenotypic Plasticity and Evolution on Population Persistence in a Changing Climate

Abstract: Climate change affects individual organisms by altering development, physiology, behavior, and fitness, and populations by altering genetic and phenotypic composition, vital rates, and dynamics. We sought to clarify how selection, phenotypic plasticity, and demography are linked in the context of climate change. On the basis of theory and results of recent empirical studies of plants and animals, we believe the ecological and evolutionary issues relevant to population persistence as climate changes are the rate, type, magnitude, and spatial pattern of climate‐induced abiotic and biotic change; generation time and life history of the organism; extent and type of phenotypic plasticity; amount and distribution of adaptive genetic variation across space and time; dispersal potential; and size and connectivity of subpopulations. An understanding of limits to plasticity and evolutionary potential across traits, populations, and species and feedbacks between adaptive and demographic responses is lacking. Integrated knowledge of coupled ecological and evolutionary mechanisms will increase understanding of the resilience and probabilities of persistence of populations and species.     

Inferring the nature of anthropogenic threats from long‐term abundance records

Diagnosing the processes that threaten species persistence is critical for recovery planning and risk forecasting. Dominant threats are typically inferred by experts on the basis of a patchwork of informal methods. Transparent, quantitative diagnostic tools would contribute much‐needed consistency, objectivity, and rigor to the process of diagnosing anthropogenic threats. Long‐term census records, available for an increasingly large and diverse set of taxa, may exhibit characteristic signatures of specific threatening processes and thereby provide information for threat diagnosis. We developed a flexible Bayesian framework for diagnosing threats on the basis of long‐term census records and diverse ancillary sources of information. We tested this framework with simulated data from artificial populations subjected to varying degrees of exploitation and habitat loss and several real‐world abundance time series for which threatening processes are relatively well understood: bluefin tuna (Thunnus maccoyii) and Atlantic cod (Gadus morhua) (exploitation) and Red Grouse (Lagopus lagopus scotica) and Eurasian Skylark (Alauda arvensis) (habitat loss). Our method correctly identified the process driving population decline for over 90% of time series simulated under moderate to severe threat scenarios. Successful identification of threats approached 100% for severe exploitation and habitat loss scenarios. Our method identified threats less successfully when threatening processes were weak and when populations were simultaneously affected by multiple threats. Our method selected the presumed true threat model for all real‐world case studies, although results were somewhat ambiguous in the case of the Eurasian Skylark. In the latter case, incorporation of an ancillary source of information (records of land‐use change) increased the weight assigned to the presumed true model from 70% to 92%, illustrating the value of the proposed framework in bringing diverse sources of information into a common rigorous framework. Ultimately, our framework may greatly assist conservation organizations in documenting threatening processes and planning species recovery. Inferencia la Naturaleza de las Amenazas Antropogénicas para los Registros de Abundancia a Largo Plazo     

Use of Inverse Spatial Conservation Prioritization to Avoid Biological Diversity Loss Outside Protected Areas

Globally expanding human land use sets constantly increasing pressure for maintenance of biological diversity and functioning ecosystems. To fight the decline of biological diversity, conservation science has broken ground with methods such as the operational model of systematic conservation planning (SCP), which focuses on design and on‐the‐ground implementation of conservation areas. The most commonly used method in SCP is reserve selection that focuses on the spatial design of reserve networks and their expansion. We expanded these methods by introducing another form of spatial allocation of conservation effort relevant for land‐use zoning at the landscape scale that avoids negative ecological effects of human land use outside protected areas. We call our method inverse spatial conservation prioritization. It can be used to identify areas suitable for economic development while simultaneously limiting total ecological and environmental effects of that development at the landscape level by identifying areas with highest economic but lowest ecological value. Our method is not based on a priori targets, and as such it is applicable to cases where the effects of land use on, for example, individual species or ecosystem types are relatively small and would not lead to violation of regional or national conservation targets. We applied our method to land‐use allocation to peat mining. Our method identified a combination of profitable production areas that provides the needed area for peat production while retaining most of the landscape‐level ecological value of the ecosystem. The results of this inverse spatial conservation prioritization are being used in land‐use zoning in the province of Central Finland.     

Setting population targets for mammals using body mass as a predictor of population persistence

Conservation planning and biodiversity assessments need quantitative targets to optimize planning options and assess the adequacy of current species protection. However, targets aiming at persistence require population‐specific data, which limit their use in favor of fixed and nonspecific targets, likely leading to unequal distribution of conservation efforts among species. We devised a method to derive equitable population targets; that is, quantitative targets of population size that ensure equal probabilities of persistence across a set of species and that can be easily inferred from species‐specific traits. In our method, we used models of population dynamics across a range of life‐history traits related to species’ body mass to estimate minimum viable population targets. We applied our method to a range of body masses of mammals, from 2 g to 3825 kg. The minimum viable population targets decreased asymptotically with increasing body mass and were on the same order of magnitude as minimum viable population estimates from species‐ and context‐specific studies. Our approach provides a compromise between pragmatic, nonspecific population targets and detailed context‐specific estimates of population viability for which only limited data are available. It enables a first estimation of species‐specific population targets based on a readily available trait and thus allows setting equitable targets for population persistence in large‐scale and multispecies conservation assessments and planning.     

Aligning science and policy to achieve evolutionarily enlightened conservation

There is increasing recognition among conservation scientists that long‐term conservation outcomes could be improved through better integration of evolutionary theory into management practices. Despite concerns that the importance of key concepts emerging from evolutionary theory (i.e., evolutionary principles and processes) are not being recognized by managers, there has been little effort to determine the level of integration of evolutionary theory into conservation policy and practice. We assessed conservation policy at 3 scales (international, national, and provincial) on 3 continents to quantify the degree to which key evolutionary concepts, such as genetic diversity and gene flow, are being incorporated into conservation practice. We also evaluated the availability of clear guidance within the applied evolutionary biology literature as to how managers can change their management practices to achieve better conservation outcomes. Despite widespread recognition of the importance of maintaining genetic diversity, conservation policies provide little guidance about how this can be achieved in practice and other relevant evolutionary concepts, such as inbreeding depression, are mentioned rarely. In some cases the poor integration of evolutionary concepts into management reflects a lack of decision‐support tools in the literature. Where these tools are available, such as risk‐assessment frameworks, they are not being adopted by conservation policy makers, suggesting that the availability of a strong evidence base is not the only barrier to evolutionarily enlightened management. We believe there is a clear need for more engagement by evolutionary biologists with policy makers to develop practical guidelines that will help managers make changes to conservation practice. There is also an urgent need for more research to better understand the barriers to and opportunities for incorporating evolutionary theory into conservation practice.     

Mapping Human and Social Dimensions of Conservation Opportunity for the Scheduling of Conservation Action on Private Land

Abstract Spatial prioritization techniques are applied in conservation‐planning initiatives to allocate conservation resources. Although typically they are based on ecological data (e.g., species, habitats, ecological processes), increasingly they also include nonecological data, mostly on the vulnerability of valued features and economic costs of implementation. Nevertheless, the effectiveness of conservation actions implemented through conservation‐planning initiatives is a function of the human and social dimensions of social‐ecological systems, such as stakeholders’ willingness and capacity to participate. We assessed human and social factors hypothesized to define opportunities for implementing effective conservation action by individual land managers (those responsible for making day‐to‐day decisions on land use) and mapped these to schedule implementation of a private land conservation program. We surveyed 48 land managers who owned 301 land parcels in the Makana Municipality of the Eastern Cape province in South Africa. Psychometric statistical and cluster analyses were applied to the interview data so as to map human and social factors of conservation opportunity across a landscape of regional conservation importance. Four groups of landowners were identified, in rank order, for a phased implementation process. Furthermore, using psychometric statistical techniques, we reduced the number of interview questions from 165 to 45, which is a preliminary step toward developing surrogates for human and social factors that can be developed rapidly and complemented with measures of conservation value, vulnerability, and economic cost to more‐effectively schedule conservation actions. This work provides conservation and land management professionals direction on where and how implementation of local‐scale conservation should be undertaken to ensure it is feasible.     

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     

Estimating Extinction Risk with Metapopulation Models of Large‐Scale Fragmentation

Habitat loss is the principal threat to species. How much habitat remains—and how quickly it is shrinking—are implicitly included in the way the International Union for Conservation of Nature determines a species’ risk of extinction. Many endangered species have habitats that are also fragmented to different extents. Thus, ideally, fragmentation should be quantified in a standard way in risk assessments. Although mapping fragmentation from satellite imagery is easy, efficient techniques for relating maps of remaining habitat to extinction risk are few. Purely spatial metrics from landscape ecology are hard to interpret and do not address extinction directly. Spatially explicit metapopulation models link fragmentation to extinction risk, but standard models work only at small scales. Counterintuitively, these models predict that a species in a large, contiguous habitat will fare worse than one in 2 tiny patches. This occurs because although the species in the large, contiguous habitat has a low probability of extinction, recolonization cannot occur if there are no other patches to provide colonists for a rescue effect. For 4 ecologically comparable bird species of the North Central American highland forests, we devised metapopulation models with area‐weighted self‐colonization terms; this reflected repopulation of a patch from a remnant of individuals that survived an adverse event. Use of this term gives extra weight to a patch in its own rescue effect. Species assigned least risk status were comparable in long‐term extinction risk with those ranked as threatened. This finding suggests that fragmentation has had a substantial negative effect on them that is not accounted for in their Red List category. Estimación del Riesgo de Extinción Mediante Modelos Metapoblacionales de Fragmentación a Gran Escala     

Implications of different population model structures for management of threatened plants

Population viability analysis (PVA) is a reliable tool for ranking management options for a range of species despite parameter uncertainty. No one has yet investigated whether this holds true for model uncertainty for species with complex life histories and for responses to multiple threats. We tested whether a range of model structures yielded similar rankings of management and threat scenarios for 2 plant species with complex postfire responses. We examined 2 contrasting species from different plant functional types: an obligate seeding shrub and a facultative resprouting shrub. We exposed each to altered fire regimes and an additional, species‐specific threat. Long‐term demographic data sets were used to construct an individual‐based model (IBM), a complex stage‐based model, and a simple matrix model that subsumes all life stages into 2 or 3 stages. Agreement across models was good under some scenarios and poor under others. Results from the simple and complex matrix models were more similar to each other than to the IBM. Results were robust across models when dominant threats are considered but were less so for smaller effects. Robustness also broke down as the scenarios deviated from baseline conditions, likely the result of a number of factors related to the complexity of the species’ life history and how it was represented in a model. Although PVA can be an invaluable tool for integrating data and understanding species’ responses to threats and management strategies, this is best achieved in the context of decision support for adaptive management alongside multiple lines of evidence and expert critique of model construction and output.     

Using empirical models of species colonization under multiple threatening processes to identify complementary threat‐mitigation strategies

Approaches to prioritize conservation actions are gaining popularity. However, limited empirical evidence exists on which species might benefit most from threat mitigation and on what combination of threats, if mitigated simultaneously, would result in the best outcomes for biodiversity. We devised a way to prioritize threat mitigation at a regional scale with empirical evidence based on predicted changes to population dynamics—information that is lacking in most threat‐management prioritization frameworks that rely on expert elicitation. We used dynamic occupancy models to investigate the effects of multiple threats (tree cover, grazing, and presence of an hyperaggressive competitor, the Noisy Miner (Manorina melanocephala) on bird‐population dynamics in an endangered woodland community in southeastern Australia. The 3 threatening processes had different effects on different species. We used predicted patch‐colonization probabilities to estimate the benefit to each species of removing one or more threats. We then determined the complementary set of threat‐mitigation strategies that maximized colonization of all species while ensuring that redundant actions with little benefit were avoided. The single action that resulted in the highest colonization was increasing tree cover, which increased patch colonization by 5% and 11% on average across all species and for declining species, respectively. Combining Noisy Miner control with increasing tree cover increased species colonization by 10% and 19% on average for all species and for declining species respectively, and was a higher priority than changing grazing regimes. Guidance for prioritizing threat mitigation is critical in the face of cumulative threatening processes. By incorporating population dynamics in prioritization of threat management, our approach helps ensure funding is not wasted on ineffective management programs that target the wrong threats or species.     

The value of flexibility in conservation financing

Land‐acquisition strategies employed by conservation organizations vary in their flexibility. Conservation‐planning theory largely fails to reflect this by presenting models that are either extremely inflexible—parcel acquisitions are irreversible and budgets are fixed—or extremely flexible—previously acquired parcels can readily be sold. This latter approach, the selling of protected areas, is infeasible or problematic in many situations. We considered the value to conservation organizations of increasing the flexibility of their land‐acquisition strategies through their approach to financing deals. Specifically, we modeled 2 acquisition‐financing methods commonly used by conservation organizations: borrowing and budget carry‐over. Using simulated data, we compared results from these models with those from an inflexible fixed‐budget model and an extremely flexible selling model in which previous acquisitions could be sold to fund new acquisitions. We then examined 3 case studies of how conservation organizations use borrowing and budget carry‐over in practice. Model comparisons showed that borrowing and budget carry‐over always returned considerably higher rewards than the fixed‐budget model. How they performed relative to the selling model depended on the relative conservation value of past acquisitions. Both the models and case studies showed that incorporating flexibility through borrowing or budget carry‐over gives conservation organizations the ability to purchase parcels of higher conservation value than when budgets are fixed without the problems associated with the selling of protected areas.     

Toward Best Practices for Developing Regional Connectivity Maps

Abstract: To conserve ecological connectivity (the ability to support animal movement, gene flow, range shifts, and other ecological and evolutionary processes that require large areas), conservation professionals need coarse‐grained maps to serve as decision‐support tools or vision statements and fine‐grained maps to prescribe site‐specific interventions. To date, research has focused primarily on fine‐grained maps (linkage designs) covering small areas. In contrast, we devised 7 steps to coarsely map dozens to hundreds of linkages over a large area, such as a nation, province, or ecoregion. We provide recommendations on how to perform each step on the basis of our experiences with 6 projects: California Missing Linkages (2001), Arizona Wildlife Linkage Assessment (2006), California Essential Habitat Connectivity (2010), Two Countries, One Forest (northeastern United States and southeastern Canada) (2010), Washington State Connected Landscapes (2010), and the Bhutan Biological Corridor Complex (2010). The 2 most difficult steps are mapping natural landscape blocks (areas whose conservation value derives from the species and ecological processes within them) and determining which pairs of blocks can feasibly be connected in a way that promotes conservation. Decision rules for mapping natural landscape blocks and determining which pairs of blocks to connect must reflect not only technical criteria, but also the values and priorities of stakeholders. We recommend blocks be mapped on the basis of a combination of naturalness, protection status, linear barriers, and habitat quality for selected species. We describe manual and automated procedures to identify currently functioning or restorable linkages. Once pairs of blocks have been identified, linkage polygons can be mapped by least‐cost modeling, other approaches from graph theory, or individual‐based movement models. The approaches we outline make assumptions explicit, have outputs that can be improved as underlying data are improved, and help implementers focus strictly on ecological connectivity.     

The role of life histories and trophic interactions in population recovery

Factors affecting population recovery from depletion are at the focus of wildlife management. Particularly, it has been debated how life‐history characteristics might affect population recovery ability and productivity. Many exploited fish stocks have shown temporal changes towards earlier maturation and reduced adult body size, potentially owing to evolutionary responses to fishing. Whereas such life‐history changes have been widely documented, their potential role on stock's ability to recover from exploitation often remains ignored by traditional fisheries management. We used a marine ecosystem model parameterized for Southeastern Australian ecosystem to explore how changes towards “faster” life histories might affect population per capita growth rate r. We show that for most species changes towards earlier maturation during fishing have a negative effect (3–40% decrease) on r during the recovery phase. Faster juvenile growth and earlier maturation were beneficial early in life, but smaller adult body sizes reduced the lifetime reproductive output and increased adult natural mortality. However, both at intra‐ and inter‐specific level natural mortality and trophic position of the species were as important in determining r as species longevity and age of maturation, suggesting that r cannot be predicted from life‐history traits alone. Our study highlights that factors affecting population recovery ability and productivity should be explored in a multi‐species context, where both age‐specific fecundity and survival schedules are addressed simultaneously. It also suggests that contemporary life‐history changes in harvested species are unlikely to increase their resilience and recovery ability.