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.     

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.     

Postdispersal seed predation limits the abundance of a long-lived perennial forb (Lithospermum ruderale)

Loss of seeds to consumers is common in plant communities, but the degree to which these losses influence plant abundance or population growth is often unclear. This is particularly the case for postdispersal seed predation by rodents, as most studies of rodent seed predation have focused on the sources of spatiotemporal variation in seed loss but not quantified the population consequences of this loss. In previous work we showed that seed predation by deer mice (Peromyscus maniculatus) substantially reduced seedling recruitment and establishment of Lithospermum ruderale (Boraginaceae), a long-lived perennial forb. To shed light on how rodent seed predation and the near-term effects on plant recruitment might influence longer-term patterns of L. ruderale population growth, we combined experimental results with demographic data in stage-based population models. Model outputs revealed that rodent seed predation had a significant impact on L. ruderale population growth rate (lambda). With the removal of postdispersal seed predation, the projected population growth rates increased between 0.06 and 0.12, depending on site (mean deltalambda across sites = 0.08). Seed predation shifted the projected stable stage distribution of populations from one with a high proportion of young plants to one in which larger adult size classes dominate. Elasticities of vital rates also changed, with germination and growth of seedlings and young plants becoming more important with the removal of seed predation. Simulations varying the magnitude of seed predation pressure while holding other vital rates constant showed that seed predation could lower lambda even if only 40% of available seeds were consumed. These results demonstrate that rodent granivory can be a potent force limiting the abundance of a long-lived perennial forb.     

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.     

Decoupling of component and ensemble density feedbacks in birds and mammals

A component density feedback represents the effect of change in population size on single demographic rates, whereas an ensemble density feedback captures that effect on the overall growth rate of a population. Given that a population's growth rate is a synthesis of the interplay of all demographic rates operating in a population, we test the hypothesis that the strength of ensemble density feedback must augment with increasing strength of component density feedback, using long-term censuses of population size, fertility, and survival rates of 109 bird and mammal populations (97 species). We found that compensatory and depensatory component feedbacks were common (each detected in approximately 50% of the demographic rates). However, component feedback strength only explained <10% of the variation in ensemble feedback strength. To explain why, we illustrate the different sources of decoupling between component and ensemble feedbacks. We argue that the management of anthropogenic impacts on populations using component feedbacks alone is ill-advised, just as managing on the basis of ensemble feedbacks without a mechanistic understanding of the contributions made by its components and environmental variability can lead to suboptimal decisions.     

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.     

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.     

Effects of Harvest of Nontimber Forest Products and Ecological Differences between Sites on the Demography of African Mahogany

Abstract: The demographic impacts of harvesting nontimber forest products (NTFP) have been increasingly studied because of reports of potentially unsustainable harvest. Nevertheless, our understanding of how plant demographic response to harvest is altered by variation in ecological conditions, which is critical for developing realistic sustainable‐use plans, is limited. We built matrix population models to test whether and how variation in ecological conditions affects population responses to harvest. In particular, we examined the effect of bark and foliage harvest on the demography of populations of African mahogany (Khaya senegalensis) in two contrasting ecological regions of Benin, West Africa. K. senegalensis bark and foliage harvest significantly reduced its stochastic population growth rates, but ecological differences between regions had a greater effect on population growth rates than did harvest. The effect of harvest on population growth rates (Δλ) was slightly stronger in the moist than in the drier region. Life‐table response experiments revealed that the mechanism by which harvesting reduced λ differed between ecological regions. Lowered stasis (persistence) of larger life stages lead to a reduction in λ in the drier region, whereas lowered growth of all life stages lowered λ in moist region. Potential strategies to increase population growth rates should include decreasing the proportion of individuals harvested, promoting harvester‐owned plantations of African mahogany, and increasing survival and growth by promoting no‐fire zones in gallery forests. Our results show how population responses to harvest of NTFP may be altered by ecological differences across sites and emphasize the importance of monitoring populations over the climatic range in which they occur to develop more realistic recommendations for conservation.     

Hierarchical demographic approaches for assessing invasion dynamics of non-indigenous species: An example using northern snakehead (Channa argus)

Models of species’ demographic features are commonly used to understand population dynamics and inform management tactics. Hierarchical demographic models are ideal for the assessment of non-indigenous species because our knowledge of non-indigenous populations is usually limited, data on demographic traits often come from a species’ native range, these traits vary among populations, and traits are likely to vary considerably over time as species adapt to new environments. Hierarchical models readily incorporate this spatiotemporal variation in species’ demographic traits by representing demographic parameters as multi-level hierarchies. As is done for traditional non-hierarchical matrix models, sensitivity and elasticity analyses are used to evaluate the contributions of different life stages and parameters to estimates of population growth rate. We applied a hierarchical model to northern snakehead (Channa argus), a fish currently invading the eastern United States. We used a Monte Carlo approach to simulate uncertainties in the sensitivity and elasticity analyses and to project future population persistence under selected management tactics. We gathered key biological information on northern snakehead natural mortality, maturity and recruitment in its native Asian environment. We compared the model performance with and without hierarchy of parameters. Our results suggest that ignoring the hierarchy of parameters in demographic models may result in poor estimates of population size and growth and may lead to erroneous management advice. In our case, the hierarchy used multi-level distributions to simulate the heterogeneity of demographic parameters across different locations or situations. The probability that the northern snakehead population will increase and harm the native fauna is considerable. Our elasticity and prognostic analyses showed that intensive control efforts immediately prior to spawning and/or juvenile-dispersal periods would be more effective (and probably require less effort) than year-round control efforts. Our study demonstrates the importance of considering the hierarchy of parameters in estimating population growth rate and evaluating different management strategies for non-indigenous invasive species.     

Importance of individual and environmental variation for invasive species spread: a spatial integral projection model

Plant survival, growth, and flowering are size dependent in many plant populations but also vary among individuals of the same size. This individual variation, along with variation in dispersal caused by differences in, e.g., seed release height, seed characteristics, and wind speed, is a key determinant of the spread rate of species through homogeneous landscapes. Here we develop spatial integral projection models (SIPMs) that include both demography and dispersal with continuous state variables. The advantage of this novel approach over discrete-stage spread models is that the effect of variation in plant size and size-dependent vital rates can be studied at much higher resolution. Comparing Neubert-Caswell matrix models to SIPMs allowed us to assess the importance of including individual variation in the models. As a test case we parameterized a SIPM with previously published data on the invasive monocarpic thistle Carduus nutans in New Zealand. Spread rate (c*) estimates were 34% lower than for standard spatial matrix models and stabilized with as few as seven evenly distributed size classes. The SIPM allowed us to calculate spread rate elasticities over the range of plant sizes, showing the size range of seedlings that contributed most to c* through their survival, growth and reproduction. The annual transitions of these seedlings were also the most important ones for local population growth (lambda). However, seedlings that reproduced within a year contributed relatively more to c* than to lambda. In contrast, plants that grow over several years to reach a large size and produce many more seeds, contributed relatively more to lambda than to c*. We show that matrix models pick up some of these details, while other details disappear within wide size classes. Our results show that SIPMs integrate various sources of variation much better than discrete-stage matrix models. Simpler, heuristic models, however, remain very valuable in studies where the main goal is to investigate the general impact of a life history stage on population dynamics. We conclude with a discussion of future extensions of SIPMs, including incorporation of continuous time and environmental drivers.     

Reduced pollinator service and elevated pollen limitation at the geographic range limit of an annual plant

Mutualisms are well known to influence individual fitness and the population dynamics of partner species, but little is known about whether they influence species distributions and the location of geographic range limits. Here, we examine the contribution of plant-pollinator interactions to the geographic range limit of the California endemic plant Clarkia xantiana ssp. xantiana. We show that pollinator availability declined from the center to the margin of the geographic range consistently across four years of study. This decline in pollinator availability was caused to a greater extent by variation in the abundance of generalist rather than specialist bee pollinators. Climate data suggest that patterns of precipitation in the current and previous year drove variation in bee abundance because of its effects on cues for bee emergence in the current year and the abundance of floral resources in the previous year. Experimental floral manipulations showed that marginal populations had greater outcross pollen limitation of reproduction, in parallel with the decline in pollinator abundance. Although plants are self-compatible, we found no evidence that autonomous selfing contributes to reproduction, and thus no evidence that it alleviates outcross pollen limitation in marginal populations. Furthermore, we found no association between the distance to the range edge and selfing rate, as estimated from sequence and microsatellite variation, indicating that the mating system has not evolved in response to the pollination environment at the range periphery. Overall, our results suggest that dependence on pollinators for reproduction may be an important constraint limiting range expansion in this system.     

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.     

Unprecedented Low Levels of Genetic Variation and Inbreeding Depression in an Island Population of the Black-Footed Rock-Wallaby

Abstract: It has been argued that demographic and environmental factors will cause small, isolated populations to become extinct before genetic factors have a significant negative impact. Islands provide an ideal opportunity to test this hypothesis because they often support small, isolated populations that are highly vulnerable to extinction. To assess the potential negative impact of isolation and small population size, we compared levels of genetic variation and fitness in island and mainland populations of the black-footed rock-wallaby ( Petrogale lateralis [Marsupialia: Macropodidae]). Our results indicate that the Barrow Island population of P. lateralis has unprecedented low levels of genetic variation (  H _e = 0.053, from 10 microsatellite loci) and suffers from inbreeding depression (reduced female fecundity, skewed sex ratio, increased levels of fluctuating asymmetry). Despite a long period of isolation ( ∼ 1600 generations) and small effective population size (  N _e ∼ 15), demographic and environmental factors have not yet driven this population to extinction. Nevertheless, it has been affected significantly by genetic factors. It has lost most of its genetic variation and become highly inbred (  F _e = 0.91), and it exhibits reduced fitness. Because several other island populations of P. lateralis also exhibit exceptionally low levels of genetic variation, this phenomenon may be widespread. Inbreeding in these populations is at a level associated with high rates of extinction in populations of domestic and laboratory species. Genetic factors cannot then be excluded as contributing to the extinction proneness of small, isolated populations.     

Do source-sink dynamics promote the spread of an invasive grass into a novel habitat?

Models of source-sink and other spatial patch dynamics have generated a number of ideas and predictions about species range expansion, the evolution of local adaptation, and the factors influencing population persistence, but relatively few empirical studies have applied these ideas due to the difficulty of measuring both patch-specific demography and movement rates. In this study, I used a combination of mark-recapture experiments, model fitting, and demographic approaches to ask how habitat-specific differences in population growth and dispersal affect spread of the invasive grass Aegilops triuncialis into serpentine environments. A. triuncialis germinated at lower rates but exhibited equivalent survival and greater growth in edge (extreme serpentine) than in core populations, even accounting for density differences between habitats. Estimated growth rates (lambda) for four of five edge subpopulations were strongly positive, ranging from lambda = 1.32 to 2.09 without propagule input from adjacent habitat. Local dispersal was best described by an exponential kernel, with a mean dispersal distance about twice as long on the edge (0.24-0.40 m) as in the core (0.18 m). Twenty-five percent of marked spikes in the edge were not relocated within the patch, suggesting greater rates of either seed predation or long-distance dispersal that reduced population growth. These results suggest that A. triuncialis can successfully spread into extreme serpentine habitats without sustained propagule input from adjacent populations. Further, asymmetric dispersal that may be both habitat- and density-dependent could slow growth rates on the edge. This pattern may also increase the importance of harsh edge patches as a source of long-distance dispersers.     

Emerging prion disease drives host selection in a wildlife population

Infectious diseases are increasingly recognized as an important force driving population dynamics, conservation biology, and natural selection in wildlife populations. Infectious agents have been implicated in the decline of small or endangered populations and may act to constrain population size, distribution, growth rates, or migration patterns. Further, diseases may provide selective pressures that shape the genetic diversity of populations or species. Thus, understanding disease dynamics and selective pressures from pathogens is crucial to understanding population processes, managing wildlife diseases, and conserving biological diversity. There is ample evidence that variation in the prion protein gene (PRNP) impacts host susceptibility to prion diseases. Still, little is known about how genetic differences might influence natural selection within wildlife populations. Here we link genetic variation with differential susceptibility of white-tailed deer to chronic wasting disease (CWD), with implications for fitness and disease-driven genetic selection. We developed a single nucleotide polymorphism (SNP) assay to efficiently genotype deer at the locus of interest (in the 96th codon of the PRNP gene). Then, using a Bayesian modeling approach, we found that the more susceptible genotype had over four times greater risk of CWD infection; and, once infected, deer with the resistant genotype survived 49% longer (8.25 more months). We used these epidemiological parameters in a multi-stage population matrix model to evaluate relative fitness based on genotype-specific population growth rates. The differences in disease infection and mortality rates allowed genetically resistant deer to achieve higher population growth and obtain a long-term fitness advantage, which translated into a selection coefficient of over 1% favoring the CWD-resistant genotype. This selective pressure suggests that the resistant allele could become dominant in the population within an evolutionarily short time frame. Our work provides a rare example of a quantifiable disease-driven selection process in a wildlife population, demonstrating the potential for infectious diseases to alter host populations. This will have direct bearing on the epidemiology, dynamics, and future trends in CWD transmission and spread. Understanding genotype-specific epidemiology will improve predictive models and inform management strategies for CWD-affected cervid populations.     

Living at the Edge: Local versus Positional Factors in the Long‐Term Population Dynamics of an Endangered Orchid

Abstract: Populations at the margin of geographic ranges of distribution have been considered more vulnerable than central ones, but recent reviews have caste doubt on this generalization. We examined the reproductive and demographic performance of a rare Euroasiatic orchid (Cypripedium calceolus) at its southwesterly range limit and compared our findings with those of previous studies of nine central populations at the center of the orchid's range. We sought to test the central‐marginal model and to evaluate factors involved in long‐term performance of forest Eurosiberian species with peripheral populations in southern European mountains. We characterized (structure, temporal fluctuations, herbivory, reproductive success, and recruitment at different habitats) four Pyrenean populations of C. calceolus of different sizes (5–3500 ramets) and monitored three of them for up to 13 years. Two quantitative stochastic models (count data and matrix models) were used to assess population trends and viability and the effect of herbivory. Contrary to expectations, and despite the negative effect of sporadic events of herbivory, the peripheral populations we studied (except the smallest one) performed similarly or better than populations occurring in central part of the species’ range in terms of reproductive success and population growth rates. Landscape changes over the last 50 years suggest that natural reforestation could be involved in the success of this plant at its southern limit. Forest expansion in the mountain regions of southern Europe may provide new opportunities for plants with geographic distributions centered mainly at higher latitudes and give some hope for their recovery in future scenarios dominated by biodiversity loss.     

Using demography and movement behavior to predict range expansion of the southern sea otter

In addition to forecasting population growth, basic demographic data combined with movement data provide a means for predicting rates of range expansion. Quantitative models of range expansion have rarely been applied to large vertebrates, although such tools could be useful for restoration and management of many threatened but recovering populations. Using the southern sea otter (Enhydra lutris nereis) as a case study, we utilized integro-difference equations in combination with a stage-structured projection matrix that incorporated spatial variation in dispersal and demography to make forecasts of population recovery and range recolonization. In addition to these basic predictions, we emphasize how to make these modeling predictions useful in a management context through the inclusion of parameter uncertainty and sensitivity analysis. Our models resulted in hind-cast (1989-2003) predictions of net population growth and range expansion that closely matched observed patterns. We next made projections of future range expansion and population growth, incorporating uncertainty in all model parameters, and explored the sensitivity of model predictions to variation in spatially explicit survival and dispersal rates. The predicted rate of southward range expansion (median = 5.2 km/yr) was sensitive to both dispersal and survival rates; elasticity analysis indicated that changes in adult survival would have the greatest potential effect on the rate of range expansion, while perturbation analysis showed that variation in subadult dispersal contributed most to variance in model predictions. Variation in survival and dispersal of females at the south end of the range contributed most of the variance in predicted southward range expansion. Our approach provides guidance for the acquisition of further data and a means of forecasting the consequence of specific management actions. Similar methods could aid in the management of other recovering populations.     

Estimating the functional form for the density dependence from life history data

Two contrasting approaches to the analysis of population dynamics are currently popular: demographic approaches where the associations between demographic rates and statistics summarizing the population dynamics are identified; and time series approaches where the associations between population dynamics, population density, and environmental covariates are investigated. In this paper, we develop an approach to combine these methods and apply it to detailed data from Soay sheep (Ovis aries). We examine how density dependence and climate contribute to fluctuations in population size via age- and sex-specific demographic rates, and how fluctuations in demographic structure influence population dynamics. Density dependence contributes most, followed by climatic variation, age structure fluctuations and interactions between density and climate. We then simplify the density-dependent, stochastic, age-structured demographic model and derive a new phenomenological time series which captures the dynamics better than previously selected functions. The simple method we develop has potential to provide substantial insight into the relative contributions of population and individual-level processes to the dynamics of populations in stochastic environments.     

Effects of Population Size on Seed Production and Germinability in an Endangered, Fragmented Grassland Plant

Abstract: Fragmentation and isolation of plant populations are thought to affect demographic processes such as seed production and cause reductions in fitness. I followed seed set over a 3-year period in eight populations of the endangered Rutidosis leptorrhynchoides (Asteraceae) that differed in population size from 13 to over 5000 flowering plants. Germinability of the resultant seed was also examined to determine whether small populations had lower fitness than large populations. Seed set was significantly associated with population size in 2 of the 3 years. Small populations (<30 flowering plants) produced significantly fewer seeds per head in 1994 and 1995 than did large populations (500 to over 5000 flowering plants), which did not differ significantly from one another. There was, however, substantial variation within populations. In 1993 seed production did not follow any simple relationship with population size, possibly because environmental stress from low rainfall had an overriding impact. Differences in seed germinability between populations were largely not evident, suggesting that this aspect of fitness has not declined substantially in small populations relative to large populations. This study suggests that nongenetic, demographic factors are of immediate importance to the persistence of small populations of R. leptorrhynchoides because of their potential impacts on seedling recruitment.     

The effect of summer drought on the predictability of local extinctions in a butterfly metapopulation

The ecological impacts of extreme climatic events on population dynamics and community composition are profound and predominantly negative. Using extensive data of an ecological model system, we tested whether predictions from ecological models remain robust when environmental conditions are outside the bounds of observation. We observed a 10-fold demographic decline of the Glanville fritillary butterfly (Melitaea cinxia) metapopulation on the Åland islands, Finland in the summer of 2018 and used climatic and satellite data to demonstrate that this year was an anomaly with low climatic water balance values and low vegetation productivity indices across Åland. Population growth rates were strongly associated with spatiotemporal variation in climatic water balance. Covariates shown previously to affect the extinction probability of local populations in this metapopulation were less informative when populations were exposed to severe drought during the summer months. Our results highlight the unpredictable responses of natural populations to extreme climatic events.