Demographic models inform selection of biocontrol agents for garlic mustard (Alliaria petiolata).

Nonindigenous invasive plants pose a major threat to natural communities worldwide. Biological control of weeds via selected introduction of their natural enemies can affect control over large spatial areas but also risk nontarget effects. To maximize effectiveness while minimizing risk, weed biocontrol programs should introduce the minimum number of host-specific natural enemies necessary to control an invasive nonindigenous plant. We used elasticity analysis of a matrix model to help inform biocontrol agent selection for garlic mustard (Alliaria petiolata (M. Bieb.) Cavara and Grande). The Eurasian biennial A. petiolata is considered one of the most problematic invaders of temperate forests in North America. Four weevil species in the genus Ceutorhynchus (Coleoptera: Curculionidae) are currently considered potential biocontrol agents. These species attack rosettes (C. scrobicollis), stems (C. roberti, C. alliariae), and seeds (C. constrictus) of A. petiolata. Elasticity analyses using A. petiolata demographic parameters from North America indicated that changes in the rosette-to-flowering-plant transition and changes in fecundity consistently had the greatest impact on population growth rate. These results suggest that attack by the rosette-feeder C. scrobicollis, which reduces overwintering survival, and seed or stem feeders that reduce seed output should be particularly effective. Model outcomes differed greatly as A. petiolata demographic parameters were varied within ranges observed in North America, indicating that successful control of A. petiolata populations may occur under some, but not all, conditions. Using these a priori analyses we predict: (1) rosette mortality and reduction of seed output will be the most important factors determining A. petiolata demography; (2) the root-crown feeder C. scrobicollis will have the most significant impact on A. petiolata demography; (3) releases of single control agents are unlikely to control A. petiolata across its full range of demographic variability; (4) combinations of agents that simultaneously reduce rosette survival and seed production will be required to suppress the most vigorous A. petiolata populations. These predictions can be tested using established long-term monitoring sites coupled with a designed release program. If demographic models can successfully predict biocontrol agent impact on invasive plant populations, a continued dialogue and collaboration between empirical and theoretical approaches may be the key to the development of successful biocontrol tactics for plant invaders in the future.     

Sink populations in carnivore management: cougar demography and immigration in a hunted population.

Carnivores are widely hunted for both sport and population control, especially where they conflict with human interests. It is widely believed that sport hunting is effective in reducing carnivore populations and related human-carnivore conflicts, while maintaining viable populations. However, the way in which carnivore populations respond to harvest can vary greatly depending on their social structure, reproductive strategies, and dispersal patterns. For example, hunted cougar (Puma concolor) populations have shown a great degree of resiliency. Although hunting cougars on a broad geographic scale (> 2000 km2) has reduced densities, hunting of smaller areas (i.e., game management units, < 1000 km2), could conceivably fail because of increased immigration from adjacent source areas. We monitored a heavily hunted population from 2001 to 2006 to test for the effects of hunting at a small scale (< 1000 km2) and to gauge whether population control was achieved (lambda < or = 1.0) or if hunting losses were negated by increased immigration allowing the population to remain stable or increase (lambda > or = 1.0). The observed growth rate of 1.00 was significantly higher than our predicted survival/fecundity growth rates (using a Leslie matrix) of 0.89 (deterministic) and 0.84 (stochastic), with the difference representing an 11-16% annual immigration rate. We observed no decline in density of the total population or the adult population, but a significant decrease in the average age of independent males. We found that the male component of the population was increasing (observed male population growth rate, lambda(OM) = 1.09), masking a decrease in the female component (lambda(OF) = 0.91). Our data support the compensatory immigration sink hypothesis; cougar removal in small game management areas (< 1000 km2) increased immigration and recruitment of younger animals from adjacent areas, resulting in little or no reduction in local cougar densities and a shift in population structure toward younger animals. Hunting in high-quality habitats may create an attractive sink, leading to misinterpretation of population trends and masking population declines in the sink and surrounding source areas.     

Population Viability Analysis for an Endangered Plant

Abstract: Demographic modeling is used to understand the population viability of Furbish's lousewort, Pedicularis furbishiae , a perennial plant species endemic to the St. John River Valley in northern Maine. Environment-specific summaries of demographic parameters (survivorship, growth, and fecundity) over four years, organized into stage-based projection matrices, provide predictions of future population dynamics given a deterministic extension of past conditions. Stochastic modeling, using (I) empirically observed variation in demographic parameters, and (2) estimated rates of natural catastrophes, leads to predictions of extinction probability.
P. furbishiae viability has varied widely over the study period Viable populations with finite rates of increase > 1 are found where cover is low, woody plants do not dominate, and disturbance does not occur. Rates of increase vary over time, suggesting that stochastic analyses would be realistic. Stochastic measures of population viability incorporating environmental variation suggest that early successional environments, especially wetter sites, can support viable populations in the absence of disturbance. However; observed rates of natural catastrophe dominate viability estimates of individual populations. Metapopulation dynamics feature extinction rates that are greater than recolonization rates, and may be affected by land use in the watershed Species management needs to consider a large-scale view of the riverine corridor.     

Stochastic population dynamics in populations of western terrestrial garter snakes with divergent life histories

Comparative evaluations of population dynamics in species with temporal and spatial variation in life-history traits are rare because they require long-term demographic time series from multiple populations. We present such an analysis using demographic data collected during the interval 1978-1996 for six populations of western terrestrial garter snakes (Thamnophis elegans) from two evolutionarily divergent ecotypes. Three replicate populations from a slow-living ecotype, found in mountain meadows of northeastern California, were characterized by individuals that develop slowly, mature late, reproduce infrequently with small reproductive effort, and live longer than individuals of three populations of a fast-living ecotype found at lakeshore locales. We constructed matrix population models for each of the populations based on 8-13 years of data per population and analyzed both deterministic dynamics based on mean annual vital rates and stochastic dynamics incorporating annual variation in vital rates. (1) Contributions of highly variable vital rates to fitness (lambda(s)) were buffered against the negative effects of stochastic variation, and this relationship was consistent with differences between the meadow (M-slow) and lakeshore (L-fast) ecotypes. (2) Annual variation in the proportion of gravid females had the greatest negative effect among all vital rates on lambda(s). The magnitude of variation in the proportion of gravid females and its effect on lambda(s) was greater in M-slow than L-fast populations. (3) Variation in the proportion of gravid females, in turn, depended on annual variation in prey availability, and its effect on lambda(s) was 4 23 times greater in M-slow than L-fast populations. In addition to differences in stochastic dynamics between ecotypes, we also found higher mean mortality rates across all age classes in the L-fast populations. Our results suggest that both deterministic and stochastic selective forces have affected the evolution of divergent life-history traits in the two ecotypes, which, in turn, affect population dynamics. M-slow populations have evolved life-history traits that buffer fitness against direct effects of variation in reproduction and that spread lifetime reproduction across a greater number of reproductive bouts. These results highlight the importance of long-term demographic and environmental monitoring and of incorporating temporal dynamics into empirical studies of life-history evolution.     

Elasticity analyses of size-based red and white abalone matrix models: management and conservation.

Prospective elasticity analyses have been used to aid in the management of fished species and the conservation of endangered species. Elasticities were examined for deterministic size-based matrix models of red abalone, Haliotis rufescens, and white abalone, H. sorenseni, to evaluate which size classes influenced population growth (lambda) the most. In the red abalone matrix, growth transitions were determined from a tag recapture study and grouped into nine size classes. In the white abalone matrix, abalone growth was determined from a laboratory study and grouped into five size classes. Survivorship was estimated from tag recapture data for red abalone using a Jolly-Seber model with size as a covariate and used for both red and white abalone. Reproduction estimates for both models used averages of the number of mature eggs produced by female red and white abalone in each size class from four-year reproduction studies. Population growth rate (lambda) was set to 1.0, and the first-year survival (larval survival through to the first size class) was estimated by iteration. Survival elasticities were higher than fecundity elasticities in both the red and white matrix models. The sizes classes with the greatest survival elasticities, and therefore the most influence on population growth in the model, were the sublegal red abalone (150-178 mm) and the largest white abalone size class (140-175 mm). For red abalone, the existing minimum legal size (178 mm) protects the size class the model suggests is critical to population growth. Implementation of education programs for novice divers coupled with renewed enforcement may serve to minimize incidental mortality of the critical size class. For white abalone, conservation efforts directed at restoring adults may have more of an impact on population growth than efforts focusing on juveniles. Our work is an example of how prospective elasticity analyses of size-structured matrix models can be used to quantitatively evaluate research priorities, fishery management strategies, and conservation options.     

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.     

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

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

Environment-specific elasticity and sensitivity analysis of the stochastic growth rate

For matrix population models, analyses of how sensitive the population growth rate is to changes in vital rates (i.e. perturbations) are important for studies of life history evolution as well as for management and conservation of threatened species. There are two types of sensitivity analyses corresponding to absolute (sensitivity) or relative (elasticity) changes in the vital rates and both types can be applied to both deterministic and stochastic matrix population models. To date, most empirical studies of elasticity and sensitivity of the stochastic growth rate have examined the response to perturbations in the vital rates in a complete set of possible environments. However, it is often of interest to examine the response to perturbations occurring in only a subset of the possible environments. This has been done for periodic time-varying models elsewhere, but here we describe a recently published method for calculating the environment-specific sensitivity and elasticity of the stochastic growth rate and apply this method to data. These environment-specific perturbation analyses provide a logical way of dividing the sensitivity and elasticity among the environments. They give important insight into the selection regime in different environments and can provide valuable information for making management decisions and management evaluations in stochastic environments.     

Invasion dynamics of the varnish clam (Nuttallia obscurata): a matrix demographic modeling approach

Marine invaders have become a significant threat to native biodiversity and ecosystem function. In this study, the invasion of the varnish clam (Nuttallia obscurata) in British Columbia, Canada, is investigated using a matrix modeling approach to identify the life history characteristics most crucial for population growth and to investigate population differences. Mark-recapture analyses and field collections from 2003 to 2004 were used to determine individual growth, survival rates, and fecundity for two sites. A multi-state matrix model was used to determine population growth rates and to conduct sensitivity and elasticity analyses. A life table response experiment was also used to determine what life history stage contributed most to observed differences in population growth rates. Population survey data were used in conjunction with the matrix model to determine plausible recruitment levels and to investigate recruitment scenarios. Both populations are currently declining but are likely sustainable because of the pulsed nature of large recruitment events. Survival of larger clams (>40 mm) is the most important for population growth based on elasticity and sensitivity analyses. Adult survival also had the largest influence on observed differences between site-specific population growth rates. The two populations studied differed in recruitment dynamics; one experiencing annual recruitment with higher post-settlement mortality and the other, episodic recruitment and lower post-settlement mortality. The most influential factor for the successful invasion of the varnish clam appears to be survival of the larger size classes. Therefore, any process that decreases adult survival (e.g., predation, commercial harvest) will have the greatest impact on population growth.     

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.     

Conservation of Species in Dynamic Landscapes: Divergent Fates of Silene tatarica Populations in Riparian Habitats

Abstract:  In transient environments, where local extinctions occur as a result of destruction or deterioration of the local habitat, the long-term persistence of a species requires successful colonizations at new, suitable sites. This kind of habitat tracking should be associated with the asynchronous dynamics of local populations, and it can be especially important for the conservation of rare plant species in riparian habitats. We determined spatiotemporal variation in the demography of the perennial Silene tatarica (L.) Pers. in 15 populations (1998–2003) located in periodically disturbed riparian habitats. The habitats differed according to their morphology (flat shores, slopes) and the amount of bare ground (open, intermediate, closed) along a successional gradient. We used elasticity and life-table response analyses and stochastic simulations to study the variation in population demography. Finite population growth rate was higher in intermediate habitats than in open and closed habitats. In stochastic simulations population size increased in most cases, but four populations were projected to become extinct within 12–70 years. The viability of local populations depended most on the survival and growth of juvenile individuals and on the fecundity of large fertile individuals. On a regional scale, the persistence of this species will require a viable network of local populations as protection against local extinctions caused by natural disturbances and succession. Accordingly, the long-term persistence of riparian species may depend on habitat changes; thus, their conservation requires maintenance of natural disturbance dynamics. Along regulated rivers, management activities such as the creation of open habitats for new colonization should be implemented. Similarly, these activities can be rather general requirements for the conservation of endangered species dependent on transient habitats along successional gradients.     

The triple helix of Plantago lanceolata: genetics and the environment interact to determine population dynamics

The theory of evolution via natural selection predicts that the genetic composition of wild populations changes over time in response to the environment. Different genotypes should exhibit different demographic patterns, but genetic variation in demography is often impossible to separate from environmental variation. Here, we asked if genetic variation is important in determining demographic patterns. We answer this question using a long-term field experiment combined with general linear modeling of deterministic population growth rates (lambda), deterministic life table response experiment (LTRE) analysis, and stochastic simulation of demography by paternal lineage in a short-lived perennial plant, Plantago lanceolata, in which we replicated genotypes across four cohorts using a standard breeding design. General linear modeling showed that growth rate varied significantly with year, spatial block, and sire. In LTRE analysis of all cohorts, the strongest influences on growth rate were from year x spatial block, and cohort x year x spatial block interactions. In analysis of genetics vs. temporal environmental variation, the strongest impacts on growth rate were from year and year x sire. Finally, stochastic simulation suggested different genetic composition among cohorts after 100 years, and different population growth rates when genetic differences were accounted for than when they were not. We argue that genetic variation, genotype x environment interactions, natural selection, and cohort effects should be better integrated into population ecological studies, as these processes should result in deviations from projected deterministic and stochastic population parameters.     

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.     

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.     

Conserving Slow-Growing, Long-Lived Tree Species: Input from the Demography of a Rare Understory Conifer, Taxus floridana

Abstract:  Although land preservation and promotion of successful regeneration are important conservation actions, their ability to increase population growth rates of slow-growing, long-lived trees is limited. We investigated the demography of Taxus floridana Nutt., a rare understory conifer, in three populations in different ravine forests spanning its entire geographic range along the Apalachicola River Bluffs in northern Florida (U.S.A.). We examined spatial and temporal patterns in demographic parameters and projected population growth rates by using four years of data on the recruitment and survival of seedlings and established stems, and on diameter growth from cross-sections of dead stems. All populations experienced a roughly 10-fold increase in seedling recruitment in 1996 compared with other years. The fates of seedlings and stems between 8 and 16 mm differed among populations. The fates of stems in two other size classes (the 2- to 4-mm class and the 4- to 8-mm class) differed among both populations and years. Individual stems in all populations exhibited similarly slow growth rates. Stochastic matrix models projected declines in all populations. Stochastic matrix analysis revealed the high elasticity of a measure of stochastic population growth rate to perturbations in the stasis of large reproductive stems for all populations. Additional analyses also indicated that occasional episodes of high recruitment do not greatly affect population growth rates. Conservation efforts directed at long-lived, slow-growing rare plants like Taxus floridana should both protect established reproductive individuals and further enhance survival of individuals in other life-history stages, such as juveniles, that often do not appear to contribute greatly to population growth rates.     

Interactive effects of prey and p,p'-DDE on burrowing owl population dynamics.

We used population models to explore the effects of the organochlorine contaminant p,p'-DDE and fluctuations in vole availability on the population dynamics of Burrowing Owls (Athene cunicularia). Previous work indicated an interaction between low biomass of voles in the diet and moderate levels of p,p'-DDE in Burrowing Owl eggs that led to reproductive impairment. We constructed periodic and stochastic matrix models that incorporated three vole population states observed in the field: average, peak, and crash years. We modeled varying frequencies of vole crash years and a range of impairment of owl demographic rates in vole crash years. Vole availability had a greater impact on owl population growth rate than did reproductive impairment if vole populations peaked and crashed frequently. However, this difference disappeared as the frequency of vole crash years declined to once per decade. Fecundity, the demographic rate most affected by p,p'-DDE, had less impact on population growth rate than adult or juvenile survival. A life table response experiment of time-invariant matrices for average, peak, and crash vole conditions showed that low population growth under vole crash conditions was due to low adult and juvenile survival rates, whereas the extremely high population growth under vole peak conditions was due to increased fecundity. Our results suggest that even simple models can provide useful insights into complex ecological interactions. This is particularly valuable when temporal or spatial scales preclude manipulative experimental work in the field or laboratory.     

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.     

Potential effects of environmental contamination on Yuma Myotis demography and population growth.

Unplanned natural and anthropogenic disasters provide unique opportunities for investigating the influence of perturbations on population vital rates and species recovery times. We investigated the potential effects of a major pesticide spill by comparing annual survival rates using mark-recapture techniques on a riparian bat species, Yuma Myotis (Myotis yumanensis). Demography and population dynamics for most bat species remain poorly understood despite advances in mark-recapture estimation and modeling techniques. We compared survival and population growth rates of two roost populations exposed to a large chemical (metam sodium) spill in the upper Sacramento River in Northern California with two roost populations outside the contaminated area from 1992 to 1996. Hypotheses about long-term effects of the spill on female juvenile and adult survival were tested using an information-theoretic approach (AIC). Working hypotheses included effects of age, chemical spill, and time trend on survival. Female adult survival was higher than female juvenile survival across all sites, suggesting stage-specific mortality risks. Model-averaged estimates of female juvenile survival in the contaminated area (0.50-0.74) were lower than in control roosts (0.60-0.78) for each year in the study, suggesting that the spill may have reduced juvenile survival for several years. Female adult survival (0.73-0.89) did not appear to be strongly affected by the spill during the years of the study. There was an increase in survival for both stage-classes across all populations during the study period, which may have been caused by the end of an extended drought in California in the winter of 1993. The spill-affected population was in decline for the first year of the study as indicated by an estimated growth rate (lambda) < 1, but population growth rates increased during the four-year period.     

Evaluating the importance of demographic connectivity in a marine metapopulation

Recently researchers have gone to great lengths to measure marine metapopulation connectivity via tagging, genetic, and trace-elemental fingerprinting studies. These empirical estimates of larval dispersal are key to assessing the significance of metapopulation connectivity within a demographic context, but the life-history data required to do this are rarely available. To evaluate the demographic consequences of connectivity we constructed seasonal, size-structured metapopulation matrix models for two species of mytilid mussel in San Diego County, California, USA. The self-recruitment and larval exchange terms were produced from a time series of realized connectivities derived from trace-elemental fingerprinting of larval shells during spring and fall from 2003 to 2008. Both species exhibited a strong seasonal pattern of southward movement of recruits in spring and northward movement in fall. Growth and mortality terms were estimated using mark-recapture data from representative sites for each species and subpopulation, and literature estimates of juvenile mortality. Fecundity terms were estimated using county-wide settlement data from 2006-2008; these data reveal peak reproduction and recruitment in fall for Mytilus californianus, and spring for M. galloprovincialis. Elasticity and life-stage simulation analyses were employed to identify the season- and subpopulation-specific vital rates and connectivity terms to which the metapopulation growth rate (lambda) was most sensitive. For both species, metapopulation growth was most sensitive to proportional changes in adult fecundity, survival and growth of juvenile stages, and population connectivity, in order of importance, but relatively insensitive to adult growth or survival. The metapopulation concept was deemed appropriate for both Mytilus species as exchange between the subpopulations was necessary for subpopulation persistence. However, highest metapopulation growth occurred in years when a greater proportion of recruits was retained within the predominant source subpopulation. Despite differences in habitat and planktonic duration, both species exhibited similar overall metapopulation dynamics with respect to key life stages and processes. However, different peak reproductive periods in an environment of seasonal current reversals led to different regional (subpopulation) contributions to metapopulation maintenance; this result emphasizes the importance of connectivity analysis for spatial management of coastal resources.     

Bridging the generation gap in plants: pollination, parental fecundity, and offspring demography

Despite extensive study of pollination and plant reproduction on the one hand, and of plant demography on the other, we know remarkably little about links between seed production in successive generations, and hence about long-term population consequences of variation in pollination success. We bridged this "generation gap" in Ipomopsis aggregata, a long-lived semelparous wildflower that is pollinator limited, by adding varying densities of seeds to natural populations and following resulting plants through their entire life histories. To determine whether pollen limitation of seed production constrains rate of population growth in this species, we sowed seeds into replicated plots at a density that mimics typical pollination success and spacing of flowering plants in nature, and at twice that density to mimic full pollination. Per capita offspring survival, flower production, and contribution to population increase (lambda) did not decline with sowing density in this experiment, suggesting that typical I. aggregata populations freed from pollen limitation will grow over the short term. In a second experiment we addressed whether density dependence would eventually erase the growth benefits of full pollination, by sowing a 10-fold range of seed densities that falls within extremes estimated for the natural "seed rain" that reaches the soil surface. Per capita survival to flowering and age at flowering were again unaffected by sowing density, but offspring size, per capita flower production, and lambda declined with density. Such density dependence complicates efforts to predict population dynamics over the longer term, because it changes components of the life history (in this case fecundity) as a population grows. A complete understanding of how constraints on seed production affect long-term population growth will hinge on following offspring fates at least through flowering of the first offspring generation, and doing so for a realistic range of population densities.