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.     

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.     

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 life history and reproduction on recruitment time lags in reintroductions of rare plants

Reintroductions are important components of conservation and recovery programs for rare plant species, but their long-term success rates are poorly understood. Previous reviews of plant reintroductions focused on short-term (e.g., ≤3 years) survival and flowering of founder individuals rather than on benchmarks of intergenerational persistence, such as seedling recruitment. However, short-term metrics may obscure outcomes because the unique demographic properties of reintroductions, including small size and unstable stage structure, could create lags in population growth. We used time-to-event analysis on a database of unusually well-monitored and long-term (4–28 years) reintroductions of 27 rare plant species to test whether life-history traits and population characteristics of reintroductions create time-lagged responses in seedling recruitment (i.e., recruitment time lags [RTLs]), an important benchmark of success and indicator of persistence in reintroduced populations. Recruitment time lags were highly variable among reintroductions, ranging from <1 to 17 years after installation. Recruitment patterns matched predictions from life-history theory with short-lived species (fast species) exhibiting consistently shorter and less variable RTLs than long-lived species (slow species). Long RTLs occurred in long-lived herbs, especially in grasslands, whereas short RTLs occurred in short-lived subtropical woody plants and annual herbs. Across plant life histories, as reproductive adult abundance increased, RTLs decreased. Highly variable RTLs were observed in species with multiple reintroduction events, suggesting local processes are just as important as life-history strategy in determining reintroduction outcomes. Time lags in restoration outcomes highlight the need to scale success benchmarks in reintroduction monitoring programs with plant life-history strategies and the unique demographic properties of restored populations. Drawing conclusions on the long-term success of plant reintroduction programs is premature given that demographic processes in species with slow life-histories take decades to unfold.     

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.     

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.     

Predicting [corrected] extinction risk in spite of predator-prey oscillations.

Most population viability analyses (PVA) assume that the effects of species interactions are subsumed by population-level parameters. We examine how robust five commonly used PVA models are to violations of this assumption. We develop a stochastic, stage-structured predator-prey model and simulate prey population vital rates and abundance. We then use simulated data to parameterize and estimate risk for three demographic models (static projection matrix, stochastic projection matrix, stochastic vital rate matrix) and two time series models (diffusion approximation [DA], corrupted diffusion approximation [CDA]). Model bias is measured as the absolute deviation between estimated and observed quasi-extinction risk. Our results highlight three generalities about the application of single-species models to multi-species conservation problems. First, our collective model results suggest that most single-species PVA models overestimate extinction risk when species interactions cause periodic variation in abundance. Second, the DA model produces the most (conservatively) biased risk forecasts. Finally, the CDA model is the most robust PVA to population cycles caused by species interactions. CDA models produce virtually unbiased and relatively precise risk estimates even when populations cycle strongly. High performance of simple time series models like the CDA owes to their ability to effectively partition stochastic and deterministic sources of variation in population abundance.     

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.     

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.     

Climate and predation dominate juvenile and adult recruitment in a turtle with temperature-dependent sex determination

Conditions experienced early in life can influence phenotypes in ecologically important ways, as exemplified by organisms with environmental sex determination. For organisms with temperature-dependent sex determination (TSD), variation in nest temperatures induces phenotypic variation that could impact population growth rates. In environments that vary over space and time, how does this variation influence key demographic parameters (cohort sex ratio and hatchling recruitment) in early life stages of populations exhibiting TSD? We leverage a 17-year data set on a population of painted turtles, Chrysemys picta, to investigate how spatial variation in nest vegetation cover and temporal variation in climate influence early life-history demography. We found that spatial variation in nest cover strongly influenced nest temperature and sex ratio, but was not correlated with clutch size, nest predation, total nest failure, or hatching success. Temporal variation in climate influenced percentage of total nest failure and cohort sex ratio, but not depredation rate, mean clutch size, or mean hatching success. Total hatchling recruitment in a year was influenced primarily by temporal variation in climate-independent factors, number of nests constructed, and depredation rate. Recruitment of female hatchlings was determined by stochastic variation in nest depredation and annual climate and also by the total nest production. Overall population demography depends more strongly on annual variation in climate and predation than it does on the intricacies of nest-specific biology. Finally, we demonstrate that recruitment of female hatchlings translates into recruitment of breeding females into the population, thus linking climate (and other) effects on early life stages to adult demographics.     

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

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

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.     

Joint modelling of breeding and survival in the kittiwake using frailty models

Assessment of population dynamics is central to population dynamics and conservation. In structured populations, matrix population models based on demographic data have been widely used to assess such dynamics. Although highlighted in several studies, the influence of heterogeneity among individuals in demographic parameters and of the possible correlation among these parameters has usually been ignored, mostly because of difficulties in estimating such individual-specific parameters. In the kittiwake (Rissa tridactyla), a long-lived seabird species, differences in survival and breeding probabilities among individual birds are well documented. Several approaches have been used in the animal ecology literature to establish the association between survival and breeding rates. However, most are based on observed heterogeneity between groups of individuals, an approach that seldom accounts for individual heterogeneity. Few attempts have been made to build models permitting estimation of the correlation between vital rates. For example, survival and breeding probability of individual birds were jointly modelled using logistic random effects models by [Cam, E., Link, W.A., Cooch, E.G., Monnat, J., Danchin, E., 2002. Individual covariation in life-history traits: seeing the trees despite the forest. Am. Naturalist, 159, in press]. This is the only example in wildlife animal populations we are aware of. Here we adopt the survival analysis approaches from epidemiology. We model the survival and the breeding probability jointly using a normally distributed random effect (frailty). Conditionally on this random effect, the survival time is modelled assuming a lognormal distribution, and breeding is modelled with a logistic model. Since the deaths are observed in year-intervals, we also take into account that the data are interval censored. The joint model is estimated using classic frequentist methods and also MCMC techniques in Winbugs. The association between survival and breeding attempt is quantified using the standard deviation of the random frailty parameters. We apply our joint model on a large data set of 862 birds, that was followed from 1984 to 1995 in Brittany (France). Survival is positively correlated with breeding indicating that birds with greater inclination to breed also had higher survival.     

Life-history and demographic spatial variation in Mediterranean populations of the opportunistic polychaete <Emphasis Type="Italic">Ophryotrocha labronica</Emphasis> (Polychaeta,Dorvilleidae)

The spatial scale of life-history and demographic variation was investigated in the opportunistic polychaete Ophryotrocha labronica La Greca and Bacci. Individuals were collected along the Italian coasts from three thermally different biogeographical regions of the Mediterranean Sea. For each region, populations from four harbours were considered, and for each harbour, two sites were examined. Life-history and demographic traits were investigated after one generation under a common garden experiment, and their variation at the three spatial scales was assessed. All the traits showed high variability with regard to site. A number of life-history and all demographic traits also varied according to the biogeographical region. Conversely, no differences were found between harbours, suggesting that geographical isolation did not contribute to phenotypic variation. Results confirmed the central role of local conditions for the evolution of life history in species colonizing heterogeneous environments, but they also pointed to the importance of large-scale factors in shaping the phenotypic responses of O. labronica, demonstrating the need for a multi-scale approach for obtaining a good measure of natural variation in widespread opportunistic 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.     

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.     

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.     

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.     

Population inertia and its sensitivity to changes in vital rates and population structure

Because the (st)age structure of a population may rarely be stable, studies of transient population dynamics and population momentum are becoming ever more popular. Yet, studies of "population momentum" are restricted in the sense that they describe the inertia of population size resulting from a demographic transition to the stationary population growth rate. Although rarely mentioned, inertia in population size is a general phenomenon and can be produced by any demographic transition or perturbation. Because population size is of central importance in demography, conservation, and management, formulas relating the sensitivity of population inertia to changes in underlying vital rates and population structure could provide much-needed insight into the dynamics of populations with unstable (st)age structure. Here, we derive such formulas, which are readily computable, and provide examples of their potential use in studies of life history and applied arenas of population study.     

Exploring Population Dynamics Patterns in a Rare Fish, Zingel asper, through Capture-Mark-Recapture Methods

Abstract:  We performed a capture-mark-recapture study on one of the last populations of Zingel asper , an endemic percid species of the Rhône River basin in France. The distribution of Z. asper has decreased dramatically during the last century. We sampled three sites in suitable habitats in the Beaume River. No impact of individual tagging on survival was found. The demography of the population was analyzed using capture-recapture methods that allow the estimation of survival, recruitment, and demographic growth rates. Annual survival rates were low (0.35–0.50). The level of transience was high (5% to 25%), suggesting that a significant number of individuals were highly mobile or shifted to suboptimal habitats. Seniority rates suggested random highly variable recruitment between years. The three sites had similar variation patterns in all demographic parameters, indicating broad spatial covariation in population dynamics. We found some local differences in demographic parameters, which could be linked to local habitat quality. Individual tagging allowed for the estimation of demographic parameters that improved our understanding of Z. asper population dynamics and revealed mechanisms that may affect population persistence, such as stochastic recruitment, low survival, and frequent dispersal. The fragmentation of habitat through river damming inhibits dispersal and represents a threat to the persistence of Z. asper in the Rhône basin. Our results offer evidence of the importance of dispersal in nonmigratory fishes and confirm the usefulness of individual tagging methods in rare fish demography.