Incorporating landscape stochasticity into population viability analysis.

The importance of incorporating landscape dynamics into population viability analysis (PVA) has previously been acknowledged, but the need to repeat the landscape generation process to encapsulate landscape stochasticity in model outputs has largely been overlooked. Reasons for this are that (1) there is presently no means for quantifying the relative effects of landscape stochasticity and population stochasticity on model outputs, and therefore no means for determining how to allocate simulation time optimally between the two; and (2) the process of generating multiple landscapes to incorporate landscape stochasticity is tedious and user-intensive with current PVA software. Here we demonstrate that landscape stochasticity can be an important source of variance in model outputs. We solve the technical problems with incorporating landscape stochasticity by deriving a formula that gives the optimal ratio of population simulations to landscape simulations for a given model, and by providing a computer program that incorporates the formula and automates multiple landscape generation in a dynamic landscape metapopulation (DLMP) model. Using a case study of a bird population, we produce estimates of DLMP model output parameters that are up to four times more precise than those estimated from a single landscape in the same amount of total simulation time. We use the DLMP modeling software RAMAS Landscape to run the landscape and metapopulation models, though our method is general and could be applied to any PVA platform. The results of this study should motivate DLMP modelers to consider landscape stochasticity in their analyses.     

Demographic stochasticity reduces the synchronizing effect of dispersal in predator-prey metapopulations

Dispersal may affect predator-prey metapopulations by rescuing local sink populations from extinction or by synchronizing population dynamics across the metapopulation, increasing the risk of regional extinction. Dispersal is likely influenced by demographic stochasticity, however, particularly because dispersal rates are often very low in metapopulations. Yet the effects of demographic stochasticity on predator-prey metapopulations are not well known. To that end, I constructed three models of a two-patch predator-prey system. The models constitute a hierarchy of complexity, allowing direct comparisons. Two models included demographic stochasticity (pure jump process [PJP] and stochastic differential equations [SDE]), and the third was deterministic (ordinary differential equations [ODE]). One stochastic model (PJP) treated population sizes as discrete, while the other (SDE) allowed population sizes to change continuously. Both stochastic models only produced synchronized predator-prey dynamics when dispersal was high for both trophic levels. Frequent dispersal by only predators or prey in the PJP and SDE spatially decoupled the trophic interaction, reducing synchrony of the non-dispersive species. Conversely, the ODE generated synchronized predator-prey dynamics across all dispersal rates, except when initial conditions produced anti-phase transients. These results indicate that demographic stochasticity strongly reduces the synchronizing effect of dispersal, which is ironic because demographic stochasticity is often invoked post hoc as a driver of extinctions in synchronized metapopulations.     

Quantifying the importance of local niche-based and stochastic processes to tropical tree community assembly

Although niche-based and stochastic processes, including dispersal limitation and demographic stochasticity, can each contribute to community assembly, it is difficult to quantify the relative importance of each process in natural vegetation. Here, we extend Shipley's maxent model (Community Assembly by Trait Selection, CATS) for the prediction of relative abundances to incorporate both trait-based filtering and dispersal limitation from the larger landscape and develop a statistical decomposition of the proportions of the total information content of relative abundances in local communities that are attributable to trait-based filtering, dispersal limitation, and demographic stochasticity. We apply the method to tree communities in a mature, species-rich, tropical forest in French Guiana at 1-, 0.25- and 0.04-ha scales. Trait data consisted of species' means of 17 functional traits measured over both the entire meta-community and separately in each of nine 1-ha plots. Trait means calculated separately for each site always gave better predictions. There was clear evidence of trait-based filtering at all spatial scales. Trait-based filtering was the most important process at the 1-ha scale (34%), whereas demographic stochasticity was the most important at smaller scales (37-53%). Dispersal limitation from the meta-community was less important and approximately constant across scales (-9%), and there was also an unresolved association between site-specific traits and meta-community relative abundances. Our method allows one to quantify the relative importance of local niche-based and meta-community processes and demographic stochasticity during community assembly across spatial and temporal scales.     

The Impact of Increased Environmental Stochasticity Due to Climate Change on the Dynamics of Asiatic Wild Ass

Abstract:  Theory proposes that increased environmental stochasticity negatively impacts population viability. Thus, in addition to the directional changes predicted for weather parameters under global climate change (GCC), the increase in variance of these parameters may also have a negative effect on biodiversity. As a case study, we assessed the impact of interannual variance in precipitation on the viability of an Asiatic wild ass ( Equus hemionus ) population reintroduced in Makhtesh Ramon Nature Reserve, Israel. We monitored the population from 1985 to 1999 to determine what environmental factors affect reproductive success. Annual precipitation during the year before conception, drought conditions during gestation, and population size determined reproductive success. We used the parameters derived from this model to assess population performance under various scenarios in a Leslie matrix type model with demographic and environmental stochasticity. Specifically, we used a change in the precipitation regime in our study area to formulate a GCC scenario and compared the simulated dynamics of the population with a no-change scenario. The coefficient of variation in population size under the global change scenario was 30% higher than under the no-change scenario. Minor die-offs (≥15%) following droughts increased extinction probability nearly 10-fold. Our results support the idea that an increase in environmental stochasticity due to GCC may, in itself, pose a significant threat to biodiversity.     

Modelling environmental stochasticity in adult survival for a long-lived species

Stochastic matrix population models are often used to help guide the management of animal populations. For a long-lived species, environmental stochasticity in adult survival will play an important role in determining outcomes from the model. One of the most common methods for modelling such stochasticity is to randomly select the value of adult survival for each year from a distribution with a specified mean and standard deviation. We consider four distributions that can provide realistic models for stochasticity in adult survival. For values of the mean and standard deviation that cover the range we would expect for long-lived species, all four distributions have similar shapes, with small differences in their skewness and kurtosis. This suggests that many of the outcomes from a population model will be insensitive to the choice of distribution, assuming that distribution provides a realistic model for environmental stochasticity in adult survival. For a generic age-structured model, the estimate of the long-run stochastic growth rate is almost identical for the four distributions, across this range of values for the mean and standard deviation. Model outcomes based on short-term projections, such as the probability of a decline over a 20-year period, are more sensitive to the choice of distribution.     

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 stochastic model for estimating sustainable limits to wildlife mortality in a changing world

Human-caused mortality of wildlife is a pervasive threat to biodiversity. Assessing the population-level impact of fisheries bycatch and other human-caused mortality of wildlife has typically relied upon deterministic methods. However, population declines are often accelerated by stochastic factors that are not accounted for in such conventional methods. Building on the widely applied potential biological removal (PBR) equation, we devised a new population modeling approach for estimating sustainable limits to human-caused mortality and applied it in a case study of bottlenose dolphins affected by capture in an Australian demersal otter trawl fishery. Our approach, termed sustainable anthropogenic mortality in stochastic environments (SAMSE), incorporates environmental and demographic stochasticity, including the dependency of offspring on their mothers. The SAMSE limit is the maximum number of individuals that can be removed without causing negative stochastic population growth. We calculated a PBR of 16.2 dolphins per year based on the best abundance estimate available. In contrast, the SAMSE model indicated that only 2.3–8.0 dolphins could be removed annually without causing a population decline in a stochastic environment. These results suggest that reported bycatch rates are unsustainable in the long term, unless reproductive rates are consistently higher than average. The difference between the deterministic PBR calculation and the SAMSE limits showed that deterministic approaches may underestimate the true impact of human-caused mortality of wildlife. This highlights the importance of integrating stochasticity when evaluating the impact of bycatch or other human-caused mortality on wildlife, such as hunting, lethal control measures, and wind turbine collisions. Although population viability analysis (PVA) has been used to evaluate the impact of human-caused mortality, SAMSE represents a novel PVA framework that incorporates stochasticity for estimating acceptable levels of human-caused mortality. It offers a broadly applicable, stochastic addition to the demographic toolbox to evaluate the impact of human-caused mortality on wildlife.     

Identifying key species in ecosystems with stochastic sensitivity analysis

The development of approaches to estimate the vulnerability of biological communities and ecosystems to extirpations and reductions of species is a central challenge of conservation biology. One key aim of this challenge is to develop quantitative approaches to estimate and rank interaction strengths and keystoneness of species and functional groups, i.e. to quantify the relative importance of species. Network analysis can be a powerful tool for this because certain structural aspects of ecological networks are good indicators of the mechanisms that maintain co-evolved, biotic interactions. A static view of ecological networks would lead us to focus research on highly-central species in food webs (topological key players in ecosystems). There are a variety of centrality indices, developed for several types of ecological networks (e.g. for weighted and un-weighted webs). However, truly understanding extinction and its community-wide effects requires the use of dynamic models. Deterministic dynamic models are feasible when population sizes are sufficiently large to minimize noise in the overall system. In models with small population sizes, stochasticity can be modelled explicitly. We present a stochastic simulation-based ecosystem model for identification of “dynamic key species” in situations where stochastic models are appropriate. To demonstrate this approach, we simulated ecosystem dynamics and performed sensitivity analysis using data from the Prince William Sound, Alaska ecosystem model. We then compare these results to those of purely topological analyses and deterministic dynamic (Ecosim) studies. We present the relationships between various topological and dynamic indices and discuss their biological relevance. The trophic group with the largest effect on others is nearshore demersals, the species mostly sensitive to others is halibut, and the group of both considerable effect on and sensitivity to others is juvenile herring. The most important trophic groups in our dynamical simulations appear to have intermediate trophic levels.     

Movement Corridors: Conservation Bargains or Poor Investments?

Corridors for movement of organisms between refuges are confounded with corridors designed for other functions, obscuring an assessment of cost-effectiveness. The rationales for movement corridors are (1) to lower extinction rate in the sense of the equilibrium theory, (2) to lessen demographic stochasticity, (3) to stem inbreeding depression, and (4) to fulfill an inherent need for movement. There is a paucity of data showing how corridors are used and whether this use lessens extinction by solving these problems. Small, isolated populations need not be doomed to quick extinction from endogenous forces such as inbreeding depression or demographic stochasticity, if their habitats are protected from humans. In specific instances, corridors could have biological disadvantages. Corridor proposals cannot be adequately judged generically. In spite of weak theoretical and empirical bases, numerous movement corridor projects are planned. In the State of Florida, multi-million-dollar corridor proposals are unsupported by data on which species might use the corridors and to what effect. Similarly, plans for massive corridor networks to counter extinction caused by global warming are weakly supported. Alternative approaches not mutually exclusive of corridors might be more effective, but such a judgment cannot be made without a cost-benefit analysis.     

Extinction-effective population index: incorporating life-history variations in population viability analysis

Viability status of populations is a commonly used measure for decision-making in the management of populations. One of the challenges faced by managers is the need to consistently allocate management effort among populations. This allocation should in part be based on comparison of extinction risks among populations. Unfortunately, common criteria that use minimum viable population size or count-based population viability analysis (PVA) often do not provide results that are comparable among populations, primarily because they lack consistency in determining population size measures and threshold levels of population size (e.g., minimum viable population size and quasi-extinction threshold). Here I introduce a new index called the "extinction-effective population index," which accounts for differential effects of demographic stochasticity among organisms with different life-history strategies and among individuals in different life stages. This index is expected to become a new way of determining minimum viable population size criteria and also complement the count-based PVA. The index accounts for the difference in life-history strategies of organisms, which are modeled using matrix population models. The extinction-effective population index, sensitivity, and elasticity are demonstrated in three species of Pacific salmonids. The interpretation of the index is also provided by comparing them with existing demographic indices. Finally, a measure of life-history-specific effect of demographic stochasticity is derived.     

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.     

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.     

Erosion of Heterozygosity in Fluctuating Populations

Abstract: Demographic, environmental, and genetic stochasticity threaten the persistence of isolated populations. The relative importance of these intertwining factors remains unresolved, but a common view is that random demographic and environmental events will usually drive small populations to the brink of extinction before genetic deterioration poses a serious threat. To evaluate the potential importance of genetic factors, we analyzed a model linking demographic and environmental conditions to the loss of genetic diversity in isolated populations undergoing natural levels of fluctuation. Nongenetic processes—environmental stochasticity and population demography—were modeled according to a bounded diffusion process. Genetic processes were modeled by quantifying the rate of drift according to the effective population size, which was predicted from the same parameters used to describe the nongenetic processes. We combined these models to predict the heterozygosity remaining at the time of extinction, as predicted by the nongenetic portion of the model. Our model predicts that many populations will lose most or all of their neutral genetic diversity before nongenetic random events lead to extinction. Given the abundant evidence for inbreeding depression and recent evidence for elevated extinction rates of inbred populations, our findings suggest that inbreeding may be a greater general threat to population persistence than is generally recognized. Therefore, conservation biologists should not ignore the genetic component of extinction risk when assessing species endangerment and developing recovery plans.     

Sensitivity to assumptions in models of generalist predation on a cyclic prey

Ecological theory predicts that generalist predators should damp or suppress long-term periodic fluctuations (cycles) in their prey populations and depress their average densities. However, the magnitude of these impacts is likely to vary depending on the availability of alternative prey species and the nature of ecological mechanisms driving the prey cycles. These multispecies effects can be modeled explicitly if parameterized functions relating prey consumption to prey abundance, and realistic population dynamical models for the prey, are available. These requirements are met by the interaction between the Hen Harrier (Circus cyaneus) and three of its prey species in the United Kingdom, the Meadow Pipit (Anthus pratensis), the field vole (Microtus agrestis), and the Red Grouse (Lagopus lagopus scoticus). We used this system to investigate how the availability of alternative prey and the way in which prey dynamics are modeled might affect the behavior of simple trophic networks. We generated cycles in one of the prey species (Red Grouse) in three different ways: through (1) the interaction between grouse density and macroparasites, (2) the interaction between grouse density and male grouse aggressiveness, and (3) a generic, delayed density-dependent mechanism. Our results confirm that generalist predation can damp or suppress grouse cycles, but only when the densities of alternative prey are low. They also demonstrate that diametrically opposite indirect effects between pairs of prey species can occur together in simple systems. In this case, pipits and grouse are apparent competitors, whereas voles and grouse are apparent facilitators. Finally, we found that the quantitative impacts of the predator on prey density differed among the three models of prey dynamics, and these differences were robust to uncertainty in parameter estimation and environmental stochasticity.     

Food web structure affects the extinction risk of species in ecological communities

This paper studies the effect of food web structure on the extinction risk of species. We examine 793 different six-species food web structures with different number, position and strength of trophic links and expose them to stochasticity in a model with Lotka–Volterra predator–prey dynamics. The characteristics of species (intrinsic rates of increase as well as intraspecific density dependence) are held constant, but the interactions with other species and characteristics of the food web are varied.     

Utility of Dynamic-Landscape Metapopulation Models for Sustainable Forest Management

Abstract:  We evaluated the utility of combining metapopulation models with landscape-level forest-dynamics models to assess the sustainability of forest management practices. We used the Brown Creeper ( Certhia americana ) in the boreal forests of northern Ontario as a case study. We selected the Brown Creeper as a potential indicator of sustainability because it is relatively common in the region but is dependent on snags and old trees for nesting and foraging; hence, it may be sensitive to timber harvesting. For the modeling we used RAMAS Landscape, a software package that integrates RAMAS GIS, population-modeling software, and LANDIS, forest-dynamics modeling software. Predictions about the future floristic composition and structure of the landscape under a variety of management and natural disturbance scenarios were derived using LANDIS. We modeled eight alternative forest management scenarios, ranging in intensity from no timber harvesting and a natural fire regime to intensive timber harvesting with salvage logging after fire. We predicted the response of the Brown Creeper metapopulation over a 160-year period and used future population size and expected minimum population size to compare the sustainability of the various management scenarios. The modeling methods were easy to apply and model predictions were sensitive to the differences among management scenarios, indicating that these methods may be useful for assessing and ranking the sustainability of forest management options. Primary concerns about the method are the practical difficulties associated with incorporating fire stochasticity in prediction uncertainty and the number of model assumptions that must be made and tested with sensitivity analysis. We wrote new software to help quantify the contribution of landscape stochasticity to model prediction uncertainty.     

Spatial Structure and Population Extinction: A Study with Drosophila Flies

Abstract: The total amount of habitat and also its distribution and subdivision affect the extinction probability of a resident population Two species of Drosophila are studied in spatial configurations of a single large habitat patch, single small habitat patches, and two small but connected habitat patches in which a low rate of migration, roughly one fly per generation, is possible. The single large habitat patch shows the lowest extinction rate lower than the combined rate of two small patches of the same total size. For one of the species, the "corridor" between the pair of small patches seems to produce a "rescue effect" that lowers extinction rates, probably due to a decrease in the coefficient of variation in fluctuations of the population sire in this coupled system. The systems seem to have been influenced by demographic stochasticity, based on the relationship of population size to extinction probability.     

Demographic heterogeneity, cohort selection, and population growth

Demographic heterogeneity--variation among individuals in survival and reproduction--is ubiquitous in natural populations. Structured population models address heterogeneity due to age, size, or major developmental stages. However, other important sources of demographic heterogeneity, such as genetic variation, spatial heterogeneity in the environment, maternal effects, and differential exposure to stressors, are often not easily measured and hence are modeled as stochasticity. Recent research has elucidated the role of demographic heterogeneity in changing the magnitude of demographic stochasticity in small populations. Here we demonstrate a previously unrecognized effect: heterogeneous survival in long-lived species can increase the long-term growth rate in populations of any size. We illustrate this result using simple models in which each individual's annual survival rate is independent of age but survival may differ among individuals within a cohort. Similar models, but with nonoverlapping generations, have been extensively studied by demographers, who showed that, because the more "frail" individuals are more likely to die at a young age, the average survival rate of the cohort increases with age. Within ecology and evolution, this phenomenon of "cohort selection" is increasingly appreciated as a confounding factor in studies of senescence. We show that, when placed in a population model with overlapping generations, this heterogeneity also causes the asymptotic population growth rate lambda to increase, relative to a homogeneous population with the same mean survival rate at birth. The increase occurs because, even integrating over all the cohorts in the population, the population becomes increasingly dominated by the more robust individuals. The growth rate increases monotonically with the variance in survival rates, and the effect can be substantial, easily doubling the growth rate of slow-growing populations. Correlations between parent and offspring phenotype change the magnitude of the increase in lambda, but the increase occurs even for negative parent-offspring correlations. The effect of heterogeneity in reproductive rate on lambda is quite different: growth rate increases with reproductive heterogeneity for positive parent-offspring correlation but decreases for negative parent-offspring correlation. These effects of demographic heterogeneity on lambda have important implications for population dynamics, population viability analysis, and evolution.     

The role of helpers in feeding chicks in cooperatively breeding green (red-billed) woodhoopoes

Summary Observations were made of ten green (red-billed) woodhoopoe Phoeniculus purpureus flocks during the breeding season in order to quantify the relationship between flock size and the amount of food delivered to chicks. The study period was kept short specifically to minimize the effects of environmental stochasticity. Neither woodhoopoe feeding visit rates nor the total amount of food brought to chicks increased with flock size. Although nonbreeders did not increase the net rate of food provisioning to chicks, they reduced parental input in chick rearing, and hence energy expenditure by the breeding pair. However, over an 8-year study period, which includes data for 144 flock years, this did not result in increased breeding frequency or enhanced survival of breeders. There is thus no evidence that helpers' feeding contributions to young per se influence the indirect fitness of helpers.     

A Population Viability Analysis for African Elephant (Loxodonta africana): How Big Should Reserves Be?

We present an age-structured, density-dependent model of elephant population dynamics in a fluctuating environment, drawing primarily upon the life history parameters obtained from studies in semi-arid land at Tsavo National Park, Kenya. Density regulation occurs by changes in the age of first reproduction and calving interval. We model environmental stochasticity with drought events affecting sex- and age-specific survivorships. Results indicate a maximum population growth rate of 3% per year and an equilibrium elephant density of 3.1/mile². Analysis of the demographic results and their sensitivity to changes in juvenile survivorship and drought frequencies, supported by genetic considerations, suggests that in semi-arid regions a minimum reserve size of 1000 mile² is necessary to attain a 99% probability of population persistence for 1000 years. The effect of age-independent culling on population viability is also analyzed.