Effect of stage-specific vital rates on population growth rates and effective population sizes in an endangered iteroparous plant

Effective population size (N(e)) determines the strength of genetic drift and can influence the level of genetic diversity a population can maintain. Assessing how changes in demographic rates associated with environmental variables and management actions affect N(e) thus can be crucial to the conservation of endangered species. Calculation of N(e) through demographic models makes it possible to use elasticity analyses to study this issue. The elasticity of N(e) to a given vital rate is the proportional change in N(e) associated with a proportional increase in that vital rate. In addition, demographic models can be used to study N(e) and population growth rate (λ) simultaneously. Simultaneous examination is important because some vital rates differ diametrically in their associations with λ and N(e). For example, in some cases increasing these vital rates increases λ and decreases N(e). We used elasticity analysis to study the effect of stage-specific survival and flowering rates on N(e), annual effective population size (N(a)), and λ in seven populations of the endangered plant Austrian dragonhead (Dracocephalum austriacum). In populations with λ ≥ 1, the elasticities of N(e) and N(a) were similar to those of λ. Survival rates of adults were associated with greater elasticities than survival rates of juveniles, flowering rates, or fecundity. In populations with λ < 1, N(e) and N(a) exhibited greater elasticities to juvenile than to adult vital rates. These patterns are similar to those observed in other species with similar life histories. We did not observe contrasting effects of any vital rate on λ and N(e); thus, management actions that increase the λ of populations of Austrian dragonhead will not increase genetic drift. Our results show that elasticity analyses of N(e) and N(a) can complement elasticity analysis of λ. Moreover, such analyses do not require more data than standard matrix models of population dynamics.     

Translating effects of inbreeding depression on component vital rates to overall population growth in endangered bighorn sheep

Evidence of inbreeding depression is commonly detected from the fitness traits of animals, yet its effects on population growth rates of endangered species are rarely assessed. We examined whether inbreeding depression was affecting Sierra Nevada bighorn sheep (Ovis canadensis sierrae), a subspecies listed as endangered under the U.S. Endangered Species Act. Our objectives were to characterize genetic variation in this subspecies; test whether inbreeding depression affects bighorn sheep vital rates (adult survival and female fecundity); evaluate whether inbreeding depression may limit subspecies recovery; and examine the potential for genetic management to increase population growth rates. Genetic variation in 4 populations of Sierra Nevada bighorn sheep was among the lowest reported for any wild bighorn sheep population, and our results suggest that inbreeding depression has reduced adult female fecundity. Despite this population sizes and growth rates predicted from matrix-based projection models demonstrated that inbreeding depression would not substantially inhibit the recovery of Sierra Nevada bighorn sheep populations in the next approximately 8 bighorn sheep generations (48 years). Furthermore, simulations of genetic rescue within the subspecies did not suggest that such activities would appreciably increase population sizes or growth rates during the period we modeled (10 bighorn sheep generations, 60 years). Only simulations that augmented the Mono Basin population with genetic variation from other subspecies, which is not currently a management option, predicted significant increases in population size. Although we recommend that recovery activities should minimize future losses of genetic variation, genetic effects within these endangered populations-either negative (inbreeding depression) or positive (within subspecies genetic rescue)-appear unlikely to dramatically compromise or stimulate short-term conservation efforts. The distinction between detecting the effects of inbreeding depression on a component vital rate (e.g., fecundity) and the effects of inbreeding depression on population growth underscores the importance of quantifying inbreeding costs relative to population dynamics to effectively manage endangered populations.     

Effects of age and sex ratios on offspring recruitment rates in translocated black rhinoceros

Success of animal translocations depends on improving postrelease demographic rates toward establishment and subsequent growth of released populations. Short‐term metrics for evaluating translocation success and its drivers, like postrelease survival and fecundity, are unlikely to represent longer‐term outcomes. We used information theory to investigate 25 years of data on black rhinoceros (Diceros bicornis) translocations. We used the offspring recruitment rate (ORR) of translocated females—a metric integrating survival, fecundity, and offspring recruitment at sexual maturity—to detect determinants of success. Our unambiguously best model (AICω = 0.986) predicted that ORR increases with female age at release as a function of lower postrelease adult rhinoceros sex ratio (males:females). Delay of first postrelease reproduction and failure of some females to recruit any calves to sexual maturity most influenced the pattern of ORRs, and the leading causes of recruitment failure were postrelease female death (23% of all females) and failure to calve (24% of surviving females). We recommend translocating older females (≥6 years old) because they do not exhibit the reproductive delay and low ORRs of juveniles (<4 years old) or the higher rates of recruitment failure of juveniles and young adults (4–5.9 years old). Where translocation of juveniles is necessary, they should be released into female‐biased populations, where they have higher ORRs. Our study offers the unique advantage of a long‐term analysis across a large number of replicate populations—a science‐by‐management experiment as a proxy for a manipulative experiment, and a rare opportunity, particularly for a large, critically endangered taxon such as the black rhinoceros. Our findings differ from previous recommendations, reinforce the importance of long‐term data sets and comprehensive metrics of translocation success, and suggest attention be shifted from ecological to social constraints on population growth and species recovery, particularly when translocating species with polygynous breeding systems.     

An Age-Structured Demographic Model for the Endangered Stephens' Kangaroo Rat

Effective conservation of endangered species often is hampered by inadequate knowledge of demography. We extracted information on survival and fecundity from an 18-month, live-trapping study of Dipodomys stephensi , and from this we developed an age-structured demographic model to assess population viability. Adult Stephens' kangaroo rats persisted longer than juveniles, and adult females persisted longer than adult males. Disappearance rates were high in the first months after initial capture. Thereafter, the fraction of animals persisting decreased slowly and in an approximately linear fashion on a semilogarithmic scale, suggesting age-independent mortality factors such as predation. Juvenile persistence did not differ substantially between two years of strikingly different rainfall. Onset of breeding followed the start of winter rains. Length of the breeding season, average number of litters per female, and the fraction of first-year females breeding were much greater in the year of higher rainfall. We propose a birth-pulse demographic model for D. stephensi that distinguishes juvenile and adult age classes. Temporal environmental variation can be modeled adequately with a constant survivorship schedule and variable fecundity determined by yearly precipitation. Several issues should be resolved, however, before conservation decisions are based on the model. Better estimates of juvenile survivorship are critical, the quantitative relationship between precipitation and fecundity must be determined, and the potential for density dependence and source-sink population dynamics must be evaluated.     

A noninvasive demographic assessment of sea lions based on stage-specific abundances

A pressing need exists to develop new approaches for obtaining information on demographic rates without causing further threats to imperiled animal populations. In this paper, we illustrate and apply a data-fitting technique based on quadratic programming that uses stage-specific abundance data to estimate demographic rates and asymptotic population growth rates (lambda). We used data from seven breeding colonies of California sea lions (Zalophus californianus) in the Gulf of California, Mexico. Estimates of lambda were similar to those from previous studies relying on a diffusion approximation using trends in total abundance. On average, predicted abundances were within 24% of the observed value for the inverse estimation method and within 29% of the observed value for the diffusion approximation. Our results suggest that three of the seven populations are declining (lambda < 1), but as many as six may be at risk. Elasticity and sensitivity analyses suggest that population management in most sites should focus on the protection of adults, whose survival generally contributes the most to lambda. The quadratic programming approach is a promising noninvasive technique for estimating demographic rates and assessing the viability of populations of imperiled species.     

Dependence of the Endangered Black‐Capped Vireo on Sustained Cowbird Management

Conservation‐reliant species depend on active management, even after surpassing recovery goals, for protection from persistent threats. Required management may include control of another species, habitat maintenance, or artificial recruitment. Sometimes, it can be difficult to determine whether sustained management is required. We used nonspatial stochastic population projection matrix simulation and a spatially explicit population model to estimate the effects of parasitism by a brood parasite, the Brown‐headed Cowbird (Moluthrus ater), on a population of endangered Black‐capped Vireos (Vireo atricapilla). We simulated parasitism as a percentage of breeding vireo pairs experiencing decreased fecundity due to cowbirds. We estimated maximum sustainable parasitism (i.e., highest percentage of parasitized vireo breeding pairs for which population growth is ≥1) with the nonspatial model under multiple scenarios designed to assess sensitivity to assumptions about population growth rate, demographic effects of parasitism, and spatial distribution of parasitism. We then used the spatially explicit model to estimate cumulative probabilities of the population falling below the population recovery target of 1000 breeding pairs for a range of parasitism rates under multiple scenarios. We constructed our models from data on vireos collected on the Fort Hood Military Reservation, Texas (U.S.A.). Estimates of maximum sustainable parasitism rates ranged from 9–12% in scenarios with a low (6%) vireo population growth rate to 49–60% in scenarios with a high (24%) growth rate. Sustained parasitism above 45–85%, depending on the scenario, would likely result in the Fort Hood Vireo population dropping below its recovery goal within the next 25 years. These estimates suggest that vireos, although tolerant of low parasitism rates, are a conservation‐reliant species dependent on cowbird management. Dependencia de Vireo atricapilla, Especie en Peligro, hacia el Manejo Sostenido de Moluthurs ater     

Double Allee Effects and Extinction in the Island Fox

Abstract:  An Allee effect (AE) occurs in populations when individuals suffer a decrease in fitness at low densities. If a fitness component is reduced (component AE), per capita population growth rates may decline as a consequence (demographic AE) and extinction risk is increased. The island fox ( Urocyon littoralis ) is endemic to six of the eight California Channel Islands. Population crashes have coincided with an increase in predation by Golden Eagles ( Aquila chrysaetos ). We propose that AEs could render fox populations more sensitive and may be a likely explanation for their sharp decline. We analyzed demographic data collected between 1988 and 2000 to test whether fox density (1) influences survival and reproductive rates; (2) interacts with eagle presence and affects fox fitness parameters; and (3) influences per capita fox population trends. A double component AE simultaneously influenced survival (of adults and pups) and proportion of breeding adult females. The adult survival AE was driven by predation by eagles. These component AEs led to a demographic AE. Multiple-component AEs, a predation-driven AE, and the simultaneous occurrence of both component and demographic AEs in a mammal are all previously unreported processes. Populations below 7 foxes/km² could have suboptimal population growth rates due to the demographic AE, and AEs may have contributed to the dramatic declines in three fox populations. Because fox densities in critically endangered populations are well below this level, removing Golden Eagles appears necessary to prevent a predation-driven AE. Conservationists should also be aware of AEs when planning the release of captive foxes. More generally, our findings highlight the danger of overlooking AEs in the conservation of populations of rare or threatened species.     

Management and Recovery Options for Ural River Beluga Sturgeon

Abstract: Management of declining fisheries of anadromous species sometimes relies heavily on supplementation of populations with captive breeding, despite evidence that captive breeding can have negative consequences and may not address the root cause of decline. The beluga sturgeon (Huso huso), a species threatened by the market for black caviar and reductions in habitat quality, is managed through harvest control and hatchery supplementation, with an emphasis on the latter. We used yield per recruit and elasticity analyses to evaluate the population status and current levels of fishing and to identify the life‐history stages that are the best targets for conservation of beluga of the Ural River. Harvest rates in recent years were four to five times higher than rates that would sustain population abundance. Sustainable rates of fishing mortality are similar to those for other long‐lived marine species such as sharks and mammals. Yield per recruit, which is maximized if fish are first harvested at age 31 years, would be greatly enhanced by raising minimum size limits or reducing illegal take of subadults. Improving the survival of subadult and adult females would increase population productivity by 10 times that achieved by improving fecundity and survival from egg to age 1 year (i.e., hatchery supplementation). These results suggest that reducing mortality of subadults and adult wild fish is a more effective conservation strategy than hatchery supplementation. Because genetics is not factored into hatchery management practices, supplementation may even reduce the viability of the beluga sturgeon.     

An integrated approach for predicting fates of reintroductions with demographic data from multiple populations

We devised a novel approach to model reintroduced populations whereby demographic data collected from multiple sites are integrated into a Bayesian hierarchical model. Integrating data from multiple reintroductions allows more precise population-growth projections to be made, especially for populations for which data are sparse, and allows projections that account for random site-to-site variation to be made before new reintroductions are attempted. We used data from reintroductions of the North Island Robin (Petroica longipes), an endemic New Zealand passerine, to 10 sites where non-native mammalian predators are controlled. A comparison of candidate models that we based on deviance information criterion showed that rat-tracking rate (an index of rat density) was a useful predictor of robin fecundity and adult female survival, that landscape connectivity and a binary measure of whether sites were on a peninsula were useful predictors of apparent juvenile survival (probably due to differential dispersal away from reintroduction sites), and that there was unexplained random variation among sites in all demographic rates. We used the two best supported models to estimate the finite rate of increase (λ) for populations at each of the 10 sites, and for a proposed reintroduction site, under different levels of rat control. Only three of the reintroduction sites had λ distributions completely >1 for either model. At two sites, λ was expected to be >1 if rat-tracking rates were <5%. At the other five reintroduction sites, λ was predicted to be close to 1, and it was unclear whether growth was expected. Predictions of λ for the proposed reintroduction site were less precise than for other sites because distributions incorporated the full range of site-to-site random variation in vital rates. Our methods can be applied to any species for which postrelease data on demographic rates are available and potentially can be extended to model multiple species simultaneously.     

Fitting probability models to population dynamics data

Methods for modeling population dynamics in probability using the generalized point process approach are developed. The life history of these populations is such that seasonal reproduction occurs during a short time. Several models are developed and analyzed. Data about two species: colonial spiders (Stegodyphus dumicola) and a migratory bird (wood thrush, Hylocichla mustelina) are used to estimate model parameters with appropriate log maximum likelihood functions. For the spiders, the model is fitted to provide evolutionary feasible colony size based on maximum likelihood estimates of fecundity and survival data. For the migratory bird species, a maximum likelihood estimates are derived for the fecundity and survival rates of young and adult birds and immigration rate. The presented approach allows computation of quantities of interest such as probability of extinction and average time to extinction.     

Demography of central and marginal populations of monkeyflowers (Mimulus cardinalis and M. lewisii)

Every species occupies a limited geographic area, but how spatiotemporal environmental variation affects individual and population fitness to create range limits is not well understood. Because range boundaries arise where, on average, populations are more likely to go extinct than to persist, range limits are an inherently population-level problem for which a demographic framework is useful. In this study, I compare demographic parameters and population dynamics between central and marginal populations of monkeyflowers, Mimulus cardinalis and M. lewisii, along an elevation gradient spanning both species' ranges. Central and marginal populations of both species differed in survival and fecundity. For M. lewisii, these components of fitness were higher in central than in marginal populations, but for M. cardinalis the converse was true. To assess spatiotemporal variation in population dynamics, I used transition matrix models to estimate asymptotic population growth rates (lambda) and found that population growth rates of M. lewisii were highest at the range center and reduced at the range margin. Population growth rates of M. cardinalis were highest at the range margin and greatly reduced at the range center. Life table response analysis decomposed spatiotemporal variation in lambda into contributions from each transition between life stages, finding that transitions from large nonreproductive and reproductive plants to the seed class and stasis in the reproductive class made the largest contributions to spatial differences in lambda. These transitions had only low to moderate sensitivities, indicating that differences in projected population growth rates resulted mainly from observed differences in transition matrix parameters and their underlying vital rates.     

Effects of Multiple Levels of Social Organization on Survival and Abundance

Abstract: Identifying how social organization shapes individual behavior, survival, and fecundity of animals that live in groups can inform conservation efforts and improve forecasts of population abundance, even when the mechanism responsible for group‐level differences is unknown. We constructed a hierarchical Bayesian model to quantify the relative variability in survival rates among different levels of social organization (matrilines and pods) of an endangered population of killer whales (Orcinus orca). Individual killer whales often participate in group activities such as prey sharing and cooperative hunting. The estimated age‐specific survival probabilities and survivorship curves differed considerably among pods and to a lesser extent among matrilines (within pods). Across all pods, males had lower life expectancy than females. Differences in survival between pods may be caused by a combination of factors that vary across the population's range, including reduced prey availability, contaminants in prey, and human activity. Our modeling approach could be applied to demographic rates for other species and for parameters other than survival, including reproduction, prey selection, movement, and detection probabilities.     

Determining Optimal Population Monitoring for Rare Butterflies

Abstract: Determining population viability of rare insects depends on precise, unbiased estimates of population size and other demographic parameters. We used data on the endangered St. Francis' satyr butterfly (Neonympha mitchellii francisci) to evaluate 2 approaches (mark–recapture and transect counts) for population analysis of rare butterflies. Mark–recapture analysis provided by far the greatest amount of demographic information, including estimates (and standard errors) of population size, detection, survival, and recruitment probabilities. Mark–recapture analysis can also be used to estimate dispersal and temporal variation in rates, although we did not do this here. Models of seasonal flight phenologies derived from transect counts (Insect Count Analyzer) provided an index of population size and estimates of survival and statistical uncertainty. Pollard–Yates population indices derived from transect counts did not provide estimates of demographic parameters. This index may be highly biased if detection and survival probabilities vary spatially and temporally. In terms of statistical performance, mark–recapture and Pollard–Yates indices were least variable. Mark–recapture estimates were less likely to fail than Insect Count Analyzer, but mark–recapture estimates became less precise as sampling intensity decreased. In general, count‐based approaches are less costly and less likely to cause harm to rare insects than mark–recapture. The optimal monitoring approach must reconcile these trade‐offs. Thus, mark–recapture should be favored when demographic estimates are needed, when financial resources enable frequent sampling, and when marking does not harm the insect populations. The optimal sampling strategy may use 2 sampling methods together in 1 overall sampling plan: limited mark–recapture sampling to estimate survival and detection probabilities and frequent but less expensive transect counts.     

Simulated Responses of a Forest-Interior Bird Population to Forest Management Options in Central Hardwood Forests of the United States

I used estimates of carrying capacity, survival, fecundity, and edge effects to simulate the responses of a forest-interior bird population to selection cutting clearcutting, and no timber harvest. I also modeled population sensitivity to changes in fecundity, survival, K , and edge relationships. Because model parameters are based on scant data, results should he regarded as hypotheses to be further investigated or measures of the relative impact or sensitivity (given model assumptions). Simulated population size was greater with no timber harvest than with clearcutting and greater with clearcutting than with group selection when edge effects were included in the model. Without edge effects, population levels were only slightly lower under group selection than under no timber harvest, and greater than clearcutting. Edge effects had only a small impact on population levels under clearcutting. Clearcut size did not have much effect on population levels, but longer and shorter rotation ages resulted in higher and lower population levels, respectively. The model was very sensitive to declines in mean fecundity and survival, suggesting that factors affecting mean demographic rates could be more important than local edge effects. Some methods of timber harvest may be compatible with the conservation of forest-interior birds, but better demographic data and information on habitat suitability of selectively cut forests and young even-aged stands is needed to adequately evaluate management options.     

Using Age Structure to Detect Impacts on Threatened Populations: a Case Study with Steller Sea Lions

Abstract:   A delayed response to change is often a characteristic of long-lived species and presents a major challenge to monitoring their status. However, rapid shifts in age structure can occur even while population size remains relatively static. We used time-varying matrix models to study age-structure information as a tool for improving detection of survivorship and fecundity change and status. We applied the methods to Steller sea lions (  Eumetopias jubatus ), a long-lived endangered marine mammal found throughout the North Pacific Rim. Population and newborn counts were supplemented with information on the fraction of the population that was juvenile, obtained by measuring animals in aerial photographs taken during range-wide censuses. By fitting the model to 1976–1998 data, we obtained maximum-likelihood estimates and 95% confidence intervals for juvenile survivorship, adult survivorship, and adult fecundity in the mid-1980s, late 1980s, and 1990s. We used a series of nested models to test whether the data were best fit by a model with one, two, or three temporal changes in demographic rates, and we fit the models to different lengths of data to test the number of years of data needed to detect a demographic change. The declines in the early 1980s were associated with severely low juvenile survivorship, whereas declines in the 1990s were associated with disproportionately low fecundity. We repeated these analyses, fitting only to the count data without the juvenile-fraction information, to determine whether the age-structure information changed the conclusions and/or changed the certainty and speed with which demographic-rate changes could be detected. The juvenile-fraction data substantially improved the degree to which estimates from the model were consistent with field data and significantly improved the speed and certainty with which changes in demographic rates were detected.     

Captive Breeding and Reintroduction Evaluation Criteria: a Case Study of Peninsular Bighorn Sheep

Abstract: Captive breeding and reintroduction programs are rarely evaluated, and assessment criteria vary widely. We used the following criteria to evaluate a bighorn sheep ( Ovis canadensis ) augmentation program: (1) survival and recruitment rates in the captive population, (2) survival of released animals, (3) recruitment of released animals, (4) growth rate of the reintroduced or augmented population, and (5) establishment of a viable wild population. Captive bighorn survival and recruitment was high, averaging 0.98 (SD = 0.05) and 71.0% (SD = 19.4), respectively. Annual survival of free-ranging captive-reared bighorn ( n = 73, x = 0.80, SD = 0.11) did not differ (   Z = −0.85, p = 0.40; n = 14) from survival of wild-reared bighorn ( n = 43, x = 0.81, SD = 0.12). Recruitment was unusually low for both captive-reared (  x = 13.7%, SD = 0.24) and wild-reared ewes (  x = 13.7%, SD = 0.20). Although reintroduction did not result in population growth or establishment of a viable population, it helped prevent extirpation of the reinforced deme, preserved metapopulation linkage, and aided habitat preservation. Chronic low recruitment and low adult survivorship precluded achievement of criteria 3–5. Environmental conditions in the release area also appeared to hinder program success. Standard evaluation criteria for ongoing reintroductions allow for informative assessments and facilitate comparisons needed to refine reintroduction science as a recovery tool for threatened or endangered populations.     

Estimating survival rates with time series of standing age-structure data

It has long been recognized that age-structure data contain useful information for assessing the status and dynamics of wildlife populations. For example, age-specific survival rates can be estimated with just a single sample from the age distribution of a stable, stationary population. For a population that is not stable, age-specific survival rates can be estimated using techniques such as inverse methods that combine time series of age-structure data with other demographic data. However, estimation of survival rates using these methods typically requires numerical optimization, a relatively long time series of data, and smoothing or other constraints to provide useful estimates. We developed general models for possibly unstable populations that combine time series of age-structure data with other demographic data to provide explicit maximum likelihood estimators of age-specific survival rates with as few as two years of data. As an example, we applied these methods to estimate survival rates for female bison (Bison bison) in Yellowstone National Park, USA. This approach provides a simple tool for monitoring survival rates based on age-structure data.     

Bridging gaps in demographic analysis with phylogenetic imputation

Phylogenetically informed imputation methods have rarely been applied to estimate missing values in demographic data but may be a powerful tool for reconstructing vital rates of survival, maturation, and fecundity for species of conservation concern. Imputed vital rates could be used to parameterize demographic models to explore how populations respond when vital rates are perturbed. We used standardized vital rate estimates for 50 bird species to assess the use of phylogenetic imputation to fill gaps in demographic data. We calculated imputation accuracy for vital rates of focal species excluded from the data set either singly or in combination and with and without phylogeny, body mass, and life-history trait data. We used imputed vital rates to calculate demographic metrics, including generation time, to validate the use of imputation in demographic analyses. Covariance among vital rates and other trait data provided a strong basis to guide imputation of missing vital rates in birds, even in the absence of phylogenetic information. Mean NRMSE for null and phylogenetic models differed by <0.01 except when no vital rates were available or for vital rates with high phylogenetic signal (Pagel's λ > 0.8). In these cases, including body mass and life-history trait data compensated for lack of phylogenetic information: mean normalized root mean square error (NRMSE) for null and phylogenetic models differed by <0.01 for adult survival and <0.04 for maturation rate. Estimates of demographic metrics were sensitive to the accuracy of imputed vital rates. For example, mean error in generation time doubled in response to inaccurate estimates of maturation time. Accurate demographic data and metrics, such as generation time, are needed to inform conservation planning processes, for example through International Union for Conservation of Nature Red List assessments and population viability analysis. Imputed vital rates could be useful in this context but, as for any estimated model parameters, awareness of the sensitivities of demographic model outputs to the imputed vital rates is essential.     

The importance of space, time, and stochasticity to the demography and management of Alliaria petiolata

As population modeling is increasingly called upon to guide policy and management, it is important that we understand not only the central tendencies of our study systems, but the consequences of their variation in space and time as well. The invasive plant Alliaria petiolata (garlic mustard) is actively managed in the United States and is the focus of a developing biological control program. Two weevils (Coleoptera: Curculionidae: Ceutorhynchus) that reduce fecundity (C. alliariae) and rosette survival plus fecundity (C. scrobicollis) are under consideration for release pending host specificity testing. We used a demographic modeling approach to (1) quantify variability in A. petiolata growth and vital rates and (2) assess the potential for single- or multiple-agent biocontrol to suppress growth of 12 A. petiolata populations in Illinois and Michigan studied over three plant generations. We used perturbation analyses and simulation models with stochastic environments to estimate stochastic growth rates (lambda(S)) and predict the probability of successful management using either a single biocontrol agent or two agent species together. Not all populations exhibited invasive dynamics. Estimates of lambda(S) ranged from 0.78 to 2.21 across sites, while annual, deterministic growth (lambda) varied up to sevenfold within individual sites. Given our knowledge of the biocontrol agents, this analysis suggests that C. scrobicollis alone may control A. petiolata at up to 63% of our study sites where lambda >1, with the combination of both agents predicted to succeed at 88% of sites. Across sites and years, the elasticity rankings were dependent on lambda. Reductions of rosette survival, fecundity, or germination of new seeds are predicted to cause the greatest reduction of lambda in growing populations. In declining populations, transitions affecting seed bank survival have the greatest effect on lambda. This contrasts with past analyses that varied parameters individually in an otherwise constant matrix, which may yield unrealistic predictions by decoupling natural parameter covariances. Overall, comparisons of stochastic and deterministic growth rates illustrate how analyses of individual populations or years could misguide management or fail to characterize complex traits such as invasiveness that emerge as attributes of populations rather than species.     

Cost‐Effectiveness of Translocation Options for a Threatened Waterbird

Abstract: Reintroduction of captive‐reared animals has become increasingly popular in recent decades as a conservation technique, but little is known of how demographic factors affect the success of reintroductions. We believe whether the increase in population persistence associated with reintroduction is sufficient to warrant the cost of rearing and relocating individuals should be considered as well. We examined the trade‐off between population persistence and financial cost of a reintroduction program for Crested Coots (Fulica cristata). This species was nearly extirpated from southern Europe due to unsustainable levels of hunting and reduction in amount and quality of habitat. We used a stochastic, stage‐based, single‐sex, metapopulation model with site‐specific parameters to examine the demographic effects of releasing juveniles or adults in each population for a range of durations. We parameterized the model with data from an unsuccessful reintroduction program in which juvenile captive‐bred Crested Coots were released between 2000 and 2009. Using economic data from the captive‐breeding program, we also determined whether the strategy that maximized abundance coincided with the least expensive strategy. Releasing adults resulted in slightly larger final abundance than the release of nonreproductive juveniles. Both strategies were equally poor in achieving a viable metapopulation, but releasing adults was 2–4 times more expensive than releasing juveniles. To obtain a metapopulation that would be viable for 30 years, fecundity in the wild would need to increase to the values observed in captivity and juvenile survival would need to increase to almost unity. We suggest that the most likely way to increase these vital rates is by increasing habitat quality at release sites.