Bayesian meta-analysis of demographic parameters in three small, temperate passerines

Accurate estimates of population parameters are vital for estimating extinction risk. Such parameters, however, are typically not available for threatened populations. We used a recently developed software tool based on Markov Chain Monte Carlo methods for carrying out Bayesian inference (the BUGS package) to estimate four demographic parameters; the intrinsic growth rate, the strength of density dependence, and the demographic and environmental variance, in three species of small temperate passerines from two sets of time series data taken from a dipper and a song sparrow population, and from previously obtained frequentist estimates of the same parameters in the great tit. By simultaneously modeling variation in these demographic parameters across species and using the resulting distributions as priors in the estimation for individual species, we improve the estimates for each individual species. This framework also allows us to make probabilistic statements about plausible parameter values for small passerines temperate birds in general which is often critically needed in management of species for which little or no data are available. We also discuss how our work relates to recently developed theory on dynamic stochastic population models, and finally note some important differences between frequentist and Bayesian methods.     

Demography of the Pacific walrus (Odobenus rosmarus divergens): 1974–2006

Global climate change may fundamentally alter population dynamics of many species for which baseline population parameter estimates are imprecise or lacking. Historically, the Pacific walrus is thought to have been limited by harvest, but it may become limited by global warming‐induced reductions in sea ice. Loss of sea ice, on which walruses rest between foraging bouts, may reduce access to food, thus lowering vital rates. Rigorous walrus survival rate estimates do not exist, and other population parameter estimates are out of date or have well‐documented bias and imprecision. To provide useful population parameter estimates we developed a Bayesian, hidden process demographic model of walrus population dynamics from 1974 through 2006 that combined annual age‐specific harvest estimates with five population size estimates, six standing age structure estimates, and two reproductive rate estimates. Median density independent natural survival was high for juveniles (0.97) and adults (0.99), and annual density dependent vital rates rose from 0.06 to 0.11 for reproduction, 0.31 to 0.59 for survival of neonatal calves, and 0.39 to 0.85 for survival of older calves, concomitant with a population decline. This integrated population model provides a baseline for estimating changing population dynamics resulting from changing harvests or sea ice.     

Accounting for uncertainty in dormant life stages in stochastic demographic models

Dormant life stages are often critical for population viability in stochastic environments, but accurate field data characterizing them are difficult to collect. Such limitations may translate into uncertainties in demographic parameters describing these stages, which then may propagate errors in the examination of population‐level responses to environmental variation. Expanding on current methods, we 1) apply data‐driven approaches to estimate parameter uncertainty in vital rates of dormant life stages and 2) test whether such estimates provide more robust inferences about population dynamics. We built integral projection models (IPMs) for a fire‐adapted, carnivorous plant species using a Bayesian framework to estimate uncertainty in parameters of three vital rates of dormant seeds – seed‐bank ingression, stasis and egression. We used stochastic population projections and elasticity analyses to quantify the relative sensitivity of the stochastic population growth rate (log λs) to changes in these vital rates at different fire return intervals. We then ran stochastic projections of log λs for 1000 posterior samples of the three seed‐bank vital rates and assessed how strongly their parameter uncertainty propagated into uncertainty in estimates of log λs and the probability of quasi‐extinction, P_q(t). Elasticity analyses indicated that changes in seed‐bank stasis and egression had large effects on log λs across fire return intervals. In turn, uncertainty in the estimates of these two vital rates explained > 50% of the variation in log λs estimates at several fire‐return intervals. Inferences about population viability became less certain as the time between fires widened, with estimates of P_q(t) potentially > 20% higher when considering parameter uncertainty. Our results suggest that, for species with dormant stages, where data is often limited, failing to account for parameter uncertainty in population models may result in incorrect interpretations of population viability.     

Population Dynamics of the Endangered, Long-Lived Perennial Species, Ligularia sibirica

Ligularia sibirica is a relict wetland perennial plant species of the Czech and Slovak Republic. Explaining variation in population growth rate and identifying the causes of that variation is important for effective protection of such an endangered species. Matrix models based on four years of data of 11 populations were used to identify the pattern of variation in the demographic vital rates of this species, and to examine the causes of the variation such as population size and habitat type. Further, the matrix model was used to determine the population growth rate, longevity and risk of extinction of each population and to identify the specific vital rates that most affect population growth rate. The results showed that population growth rates were significantly different between years and populations. Temporal variation was mostly due to variable survival of adult individuals, while spatial variation was mainly driven by fertility of one small currently expanding population. Further, most studied populations of L. sibirica are performing well and only those growing in nitrogen-rich habitats have a high extinction risk. The results also indicate that all populations have low adult mortality, long-lived individuals (61.3?years on average) and some populations also show features of remnant populations (i.e., the persistence of populations in severe conditions in spite of no reproduction). Our results imply that detailed demographic data are needed to understand the long-term prospects of these populations. These data may serve as an early warning system for this species long before an obvious decline occurs in the populations.     

Estimating Survival for Elusive Juvenile Pond-Breeding Salamanders

Juvenile vital rates have important effects on population dynamics for many species, but this demographic is often difficult to locate and track. As such, we frequently lack reliable estimates of juvenile survival, which are necessary for accurately assessing population stability and potential management approaches to conserve biodiversity. We estimated survival rates for elusive juveniles of 3 species, the ringed salamander (Ambystoma annulatum), spotted salamander (A. maculatum), and small-mouthed salamander (A. texanum), using 2 approaches. First, we conducted an 11-month (2016–2017) mark-recapture study within semi-natural enclosures and used Bayesian Cormack-Jolly-Seber models to estimate survival and recapture probabilities. Second, we inferred the expected annual juvenile survival rate given published vital rates for pre-metamorphic and adult ambystomatids assuming stable population growth. For all 3 species, juvenile survival probabilities were constant across recapture occasions, whereas recapture probability estimates were time-dependent. Further, survival and recapture probabilities among study species were similar. Post-study sampling revealed that the initial study period median estimate of annual survival probability (0.39) underestimated the number of salamanders known alive at 11 months. We therefore appended approximately 1 year of opportunistic data, which produced a median annual survival probability of 0.50, encompassing salamanders that we knew to have been alive. Calculation from literature values suggested a mean annual terrestrial juvenile ambystomatid survival probability of 0.49. Similar results among our approaches indicated that juvenile survival estimates for the study species were robust and likely comparable to rates in nature. These estimates can now be confidently applied to research, monitoring, and management efforts for the study species and ecologically similar taxa. Our findings indicated that similarly robust vital rate estimates for subsets of ecologically and phylogenetically similar species can provide reasonable surrogate demographic information that can be used to reveal key factors influencing population viability for data-deficient species. © 2020 The Wildlife Society.     

Stochastic LTRE analysis of the effects of herbivory on the population dynamics of a perennial grassland herb

Herbivores can have strong deleterious effects on vital rates (growth, reproduction, and survival) and thus negatively impact the population dynamics of plant species. In practice, however, these effects might be strongly correlated, for example as a result of tradeoffs between vital rates. To get better insights into the effects of herbivory on the population dynamics of the long‐lived grassland plant Primula veris population projection matrices were constructed from demographic data collected between 1999 and 2008 (nine annual transitions). Data were collected in two large grassland populations, each of which was subjected to two treatments (grazing by cattle versus a mowing treatment), yielding a total of 36 matrices. We applied a lower‐level vital rate life table response experiment (LTRE) using the small noise approximation (SNA) of the stochastic population growth rate to disentangle the contributions of changes in mean vital rates, variability in vital rates, correlations between vital rates and vital rate elasticities to the difference in the stochastic growth rate. Stochastic growth rates (a= log λ_S) were significantly lower in grazed than in mown plots (a= 0.0185 and 0.1019, respectively). SNA LTRE analysis showed that contributions of mean vital rates by far made the largest contribution to the observed difference in a between grazed and control plots. In particular, changes in sexual reproduction rates made the largest contributions to lower the stochastic growth rate in grazed plots: both adult flowering probabilities and flower and seed production were importantly lower in grazed populations, but these negative effects were largely buffered by increased establishment and seedling survival rates. Among the stochastic terms of the SNA decomposition, contributions of covariance and correlations between vital rates had the largest impact, whereas contributions of elasticities were smaller. The strongest correlation driver was the association between adult survival and seedling establishment, suggesting that environmental conditions favouring adult survival also are beneficial for seedling establishment. Overall, our results show that herbivory had a strong negative effect on the long‐term population growth rate of P. veris that was primarily mediated by differences in fecundity (flower and seed production) and germination.     

Comparative Population Biology of Critically Endangered Dracocephalum austriacum (Lamiaceae) in Two Distant Regions

Our study aims to compare the population dynamics of critically endangered species, Dracocephalum austriacum, in two distant regions (Czech and Slovak Karst) with similar habitat conditions, and to evaluate if conclusions concerning factors affecting population performance in one region could help in understanding species dynamics in the other region. Transition matrix models were used to examine population dynamics of the species. Results show that populations in the Slovak Karst are performing better than populations in the Czech Karst. The transitions contributing most to the population growth rates were those of stasis in the small and large adult plant stage. These transitions are, however, stable between populations and years, suggesting that they are not very sensitive to the ongoing changes. Transitions that contribute most to variations in population growth rate include seed and seedling production by large adult plants, seed germination and stasis and growth of small adult plants. These transitions seem to be affected by habitat conditions (soil depth and bare rock cover), genetic parameters of the populations, climate fluctuations (especially severe droughts), and shrub expansion. The transitions contributing most to the population growth rates are very similar between the regions except for contributions of seedling growth and stasis of large adults. These differences need to be considered before using information about the most important life cycle transitions in other regions. We suggest that similar results would be found for other species with narrow habitat requirements occupying very similar habitats, even in very distant regions.     

Linking vital rates to invasiveness of a perennial herb

Invaders generally show better individual performance than non-invaders and, therefore, vital rates (survival, growth, fecundity) could potentially be used to predict species invasiveness outside their native range. Comparative studies have usually correlated vital rates with the invasiveness status of species, while few studies have investigated them in relation to population growth rate. Here, I examined the influence of five vital rates (plant establishment, survival, growth, flowering probability, seed production) and their variability (across geographic regions, habitat types, population sizes and population densities) on population growth rate (λ) using data from 37 populations of an invasive, iteroparous herb (Lupinus polyphyllus) in a part of its invaded range in Finland. Variation in vital rates was often related to habitat type and population density. The performance of the populations varied from declining to rapidly increasing independently of habitat type, population size or population density, but differed between regions. The population growth rate increased linearly with plant establishment, and with the survival and growth of vegetative individuals, while the survival of flowering individuals and annual seed production were not related to λ. The vital rates responsible for rapid population growth varied among populations. These findings highlight the importance of both regional and local conditions to plant population dynamics, demonstrating that individual vital rates do not necessarily correlate with λ. Therefore, to understand the role of individual vital rates in a species ability to invade, it is necessary to quantify their effect on population growth rate.     

Variation in life history and demography of the American black bear

Variation in life history and demography across a species' range informs researchers about regional adaptations and affects whether managers can borrow information from other populations in decision-making. The American black bear (Ursus americanus) is a long-lived game species whose continued persistence depends on management of harvest and removal of habituated bears that come into conflict with humans. Understanding the demography of black bears guides efforts at management and conservation, yet detailed knowledge of many populations is typically lacking. I performed a hierarchical Bayesian meta-analysis of black bear demographic studies across the geographic range of the species to explore how vital rates vary across the range, what information they give us about population growth, and whether managers can justify borrowing information from other studies to inform management decisions. Cub, yearling, and adult survival and fecundity varied between eastern and western North America, whereas subadult survival did not show geographic structuring. Adult survival and fecundity appeared to trade off, with higher survival in the western portions of bears' range and higher fecundity in the east. Although adult survival had the highest elasticity, differences in reproduction drove differences in population growth rate. Mean population growth rate was higher in the east (0.99; 95% credible interval [CrI]: 0.96, 1.03) than the west (0.97; 95%CrI: 0.93, 1.01). Despite declining trends in the west, 34% of the distribution of population growth rate was >1, compared to 55% in the east. Further work needs to be done to address the cause of the apparent trade-off between adult survival and fecundity and explore how the estimated growth rates are likely to affect population status of black bears. Because population growth rates are close to 1 and small deviations could impact whether a population is considered increasing or decreasing, managers need to employ caution in borrowing vital rates from other populations. © 2011 The Wildlife Society.     

The generality of management recommendations across populations of an invasive perennial herb

Demographic models are widely used to produce management recommendations for different species. For invasive plants, current management recommendations to control local population growth are often based on data from a limited number of populations per species, and the assumption of stable population structure (asymptotic dynamics). However, spatial variation in population dynamics and deviation from a stable structure may affect these recommendations, calling into question their generality across populations of an invasive species. Here, I focused on intraspecific variation in population dynamics and investigated management recommendations generated by demographic models across 37 populations of a short-lived, invasive perennial herb (Lupinus polyphyllus). Models that relied on the proportional perturbations of vital rates (asymptotic elasticities) indicated an essential role for plant survival in long-term population dynamics. The rank order of elasticities for different vital rates (survival, growth, retrogression, fecundity) varied little among the 37 study populations regardless of population status (increasing or declining asymptotically). Summed elasticities for fecundity increased, while summed elasticities for survival decreased with increasing long-term population growth rate. Transient dynamics differed from asymptotic dynamics, but were qualitatively similar among populations, that is, depending on the initial size structure, populations tended to either increase or decline in density more rapidly than predicted by asymptotic growth rate. These findings indicate that although populations are likely to exhibit transient dynamics, management recommendations based on asymptotic elasticities for vital rates might be to some extent generalised across established populations of a given short-lived invasive plant species.     

Variability of important vital rates challenges the demographic buffering hypothesis

Selection is assumed to eliminate life-histories showing high variability in vital rates that have the greatest influence on population performance. Therefore, an inverse variability-importance relationship of vital rates is believed to be a universal pattern for diverse life-histories. We tested for such a relationship using multi-year demographic data on a large number of populations of two perennial plant species. Applying different approaches, we first examined the overall variability-importance relationship for the average main vital rates (survival, growth, retrogression, fecundity) per species, and then separately for each population. We found an overall inverse relationship between temporal variation and importance of the average main vital rates for both study species, but these negative species-level correlations were mainly caused by different scales of the examined vital rates. When variability-importance relationships were examined across individual demographic transitions within populations, the abundance of positive and negative correlations depended largely on the method used, and positive correlations were more common after correcting vital rates for sampling variation than when using uncorrected vital rates. Our results cast doubt on the generality of the demographic buffering hypothesis, suggesting that the inverse variability-importance relationship may not be a universal pattern when vital rates are examined for multiple populations of the same plant species.     

Population dynamics and comparative demographics in sympatric populations of the round-leaved orchids Platanthera macrophylla and P. orbiculata

Orchids (Orchidaceae) are a family of flowering plants with a high proportion of threatened taxa making them an important focus of plant conservation. Orchid conservation efforts are most effective when informed by reliable demographic research. We utilized transition matrix models to examine the population dynamics and demography within sympatric populations of a species pair of terrestrial round-leaved orchids, Platanthera macrophylla and P. orbiculata. The models were parameterized from a large data set spanning 9 years from field observations of over 1,000 orchids. Life table response experiments (LTRE) were used to identify which life history transitions, and which vital rates within those transitions, most contributed to observed differences between the two species and most contributed to interannual variation within each species. Results from mean transition matrices projected finite rates of population growth that were not significantly different between the two species, with P. macrophylla near the replacement rate and P. orbiculata below it. LTRE revealed that the difference in population growth rates between the two species was mostly due to differences in fecundity (flowering adult to protocorm transition) driven by differences in fruit set and seed germination into protocorm, which were much greater for P. macrophylla. The two primary contributors to interannual variation in population growth rates for both orchids were adult survival and fruit set, respectively. These findings indicate that any environmental disturbances harming adult survival or fecundity will have a disproportionately negative effect on the orchid populations.     

Accounting for stochasticity in demographic compensation along the elevational range of an alpine plant

Demographic compensation arises when vital rates change in opposite directions across populations, buffering the variation in population growth rates, and is a mechanism often invoked to explain the stability of species geographic ranges. However, studies on demographic compensation have disregarded the effects of temporal variation in vital rates and their temporal correlations, despite theoretical evidence that stochastic dynamics can affect population persistence in temporally varying environments. We carried out a seven‐year‐long demographic study on the perennial plant Arabis alpina (L.) across six populations encompassing most of its elevational range. We discovered demographic compensation in the form of negative correlations between the means of plant vital rates, but also between their temporal coefficients of variation, correlations and elasticities. Even if their contribution to demographic compensation was small, this highlights a previously overlooked, but potentially important, role of stochastic processes in stabilising population dynamics at range margins.     

Pattern of variation in avian population growth rates

A central question in population ecology is to understand why population growth rates differ over time. Here, we describe how the long-term growth of populations is not only influenced by parameters affecting the expected dynamics, for example form of density dependence and specific population growth rate, but is also affected by environmental and demographic stochasticity. Using long-term studies of fluctuations of bird populations, we show an interaction between the stochastic and the deterministic components of the population dynamics: high specific growth rates at small densities r(1) are typically positively correlated with the environmental variance sigma(e)(2). Furthermore, theta, a single parameter describing the form of the density regulation in the theta-logistic density-regulation model, is negatively correlated with r(1). These patterns are in turn correlated with interspecific differences in life-history characteristics. Higher specific growth rates, larger stochastic effects on the population dynamics and stronger density regulation at small densities are found in species with large clutch sizes or high adult mortality rates than in long-lived species. Unfortunately, large uncertainties in parameter estimates, as well as strong stochastic effects on the population dynamics, will often make even short-term population projections unreliable. We illustrate that the concept of population prediction interval can be useful in evaluating the consequences of these uncertainties in the population projections for the choice of management actions.     

Estimates of annual survival,growth, and recruitment of a white-tailed ptarmigan population in Colorado over 43 years

Long-term datasets for high-elevation species are rare, and considerable uncertainty exists in understanding how high-elevation populations have responded to recent climate warming. We present estimates of demographic vital rates from a 43-year population study of white-tailed ptarmigan (Lagopus leucura), a species endemic to alpine habitats in western North America. We used capture-recapture models to estimate annual rates of apparent survival, population growth, and recruitment for breeding-age ptarmigan, and we fit winter weather covariates to models in an attempt to explain annual variation. There were no trends in survival over the study period but there was strong support for age and sex effects. The average rate of annual growth suggests a relatively stable breeding-age population ( \( \bar{\lambda } \)  = 1.036), but there was considerable variation between years for both population growth and recruitment rates. Winter weather covariates only explained a small amount of variation in female survival and were not an important predictor of male survival. Cumulative winter precipitation was found to have a quadratic effect on female survival, with survival being highest during years of average precipitation. Cumulative winter precipitation was positively correlated with population growth and recruitment rates, although this covariate only explained a small amount of annual variation in these rates and there was considerable uncertainty among the models tested. Our results provide evidence for an alpine-endemic population that has not experienced extirpation or drastic declines. However, more information is needed to understand risks and vulnerabilities of warming effects on juveniles as our analysis was confined to determination of vital rates for breeding-age birds.     

Integrating broad‐scale data to assess demographic and climatic contributions to population change in a declining songbird

Climate variation and trends affect species distribution and abundance across large spatial extents. However, most studies that predict species response to climate are implemented at small spatial scales or are based on occurrence‐environment relationships that lack mechanistic detail. Here, we develop an integrated population model (IPM) for multi‐site count and capture‐recapture data for a declining migratory songbird, Wilson's warbler (Cardellina pusilla), in three genetically distinct breeding populations in western North America. We include climate covariates of vital rates, including spring temperatures on the breeding grounds, drought on the wintering range in northwest Mexico, and wind conditions during spring migration. Spring temperatures were positively related to productivity in Sierra Nevada and Pacific Northwest genetic groups, and annual changes in productivity were important predictors of changes in growth rate in these populations. Drought condition on the wintering grounds was a strong predictor of adult survival for coastal California and Sierra Nevada populations; however, adult survival played a relatively minor role in explaining annual variation in population change. A latent parameter representing a mixture of first‐year survival and immigration was the largest contributor to variation in population change; however, this parameter was estimated imprecisely, and its importance likely reflects, in part, differences in spatio‐temporal distribution of samples between count and capture‐recapture data sets. Our modeling approach represents a novel and flexible framework for linking broad‐scale multi‐site monitoring data sets. Our results highlight both the potential of the approach for extension to additional species and systems, as well as needs for additional data and/or model development.     

Effects of sample size on estimates of population growth rates calculated with matrix models

Background

Matrix models are widely used to study the dynamics and demography of populations. An important but overlooked issue is how the number of individuals sampled influences estimates of the population growth rate (λ) calculated with matrix models. Even unbiased estimates of vital rates do not ensure unbiased estimates of λ–Jensen''s Inequality implies that even when the estimates of the vital rates are accurate, small sample sizes lead to biased estimates of λ due to increased sampling variance. We investigated if sampling variability and the distribution of sampling effort among size classes lead to biases in estimates of λ.

Methodology/Principal Findings

Using data from a long-term field study of plant demography, we simulated the effects of sampling variance by drawing vital rates and calculating λ for increasingly larger populations drawn from a total population of 3842 plants. We then compared these estimates of λ with those based on the entire population and calculated the resulting bias. Finally, we conducted a review of the literature to determine the sample sizes typically used when parameterizing matrix models used to study plant demography.

Conclusions/Significance

We found significant bias at small sample sizes when survival was low (survival = 0.5), and that sampling with a more-realistic inverse J-shaped population structure exacerbated this bias. However our simulations also demonstrate that these biases rapidly become negligible with increasing sample sizes or as survival increases. For many of the sample sizes used in demographic studies, matrix models are probably robust to the biases resulting from sampling variance of vital rates. However, this conclusion may depend on the structure of populations or the distribution of sampling effort in ways that are unexplored. We suggest more intensive sampling of populations when individual survival is low and greater sampling of stages with high elasticities.     

Population Viability and Vital Rate Sensitivity of an Endangered Avian Cooperative Breeder,the White-Breasted Thrasher (Ramphocinclus brachyurus)

Social behaviors can significantly affect population viability, and some behaviors might reduce extinction risk. We used population viability analysis to evaluate effects of past and proposed habitat loss on the White-breasted Thrasher (Ramphocinclus brachyurus), a cooperatively breeding songbird with a global population size of <2000 individuals. We used an individual-based approach to build the first demographic population projection model for this endangered species, parameterizing the model with data from eight years of field study before and after habitat loss within the stronghold of the species’ distribution. The recent habitat loss resulted in an approximately 18% predicted decline in population size; this estimate was mirrored by a separate assessment using occupancy data. When mortality rates remained close to the pre-habitat loss estimate, quasi-extinction probability was low under extant habitat area, but increased with habitat loss expected after current plans for resort construction are completed. Post-habitat loss mortality rate estimates were too high for projected populations to persist. Vital rate sensitivity analyses indicated that population growth rate and population persistence were most sensitive to juvenile mortality. However, observed values for adult mortality were closest to the threshold value above which populations would crash. Adult mortality, already relatively low, may have the least capacity to change compared to other vital rates, whereas juvenile mortality may have the most capacity for improvement. Results suggest that improving mortality estimates and determining the cause(s) of juvenile mortality should be research priorities. Despite predictions that aspects of cooperative systems may result in variation in reproduction or juvenile mortality being the most sensitive vital rates, adult mortality was the most sensitive in half of the demographic models of other avian cooperative breeders. Interestingly, vital rate sensitivity differed by model type. However, studies that explicitly modeled the species’ cooperative breeding system found reproduction to be the most sensitive rate.     

Population Dynamics of the Andean Lizard Anolis heterodermus: Fast‐slow Demographic Strategies in Fragmented Scrubland Landscapes

Habitat fragmentation and loss affect population stability and demographic processes, increasing the extinction risk of species. We studied Anolis heterodermus populations inhabiting large and small Andean scrubland patches in three fragmented landscapes in the Sabana de Bogotá (Colombia) to determine the effect of habitat fragmentation and loss on population dynamics. We used the capture‐mark‐recapture method and multistate models to estimate vital rates for each population. We estimated growth population rate and the most important processes that affect λ by elasticity analysis of vital rates. We tested the effects of habitat fragmentation and loss on vital rates of lizard populations. All six isolated populations showed a positive or an equilibrium growth rate (λ = 1), and the most important demographic process affecting λ was the growth to first reproduction. Populations from landscapes with less scrubland natural cover showed higher stasis of young adults. Populations in highly fragmented landscapes showed highest juvenile survival and growth population rates. Independent of the landscape's habitat configuration and connectivity, populations from larger scrubland patches showed low adult survivorship, but high transition rates. Populations varied from a slow strategy with low growth and delayed maturation in smaller patches to a fast strategy with high growth and early maturation in large patches. This variation was congruent with the fast‐slow continuum hypothesis and has serious implications for Andean lizard conservation and management strategies. We suggest that more stable lizard populations will be maintained if different management strategies are adopted according to patch area and habitat structure.     

Buffering and plasticity in vital rates of oldfield rodents

1. Under the hypothesis of environmental buffering, populations are expected to minimize the variance of the most influential vital rates; however, this may not be a universal principle. Species with a life span <1 year may be less likely to exhibit buffering because of temporal or seasonal variability in vital rate sensitivities. Further, plasticity in vital rates may be adaptive for species in a variable environment with reliable cues. 2. We tested for environmental buffering and plasticity in vital rates using stage-structured matrix models from long-term data sets in four species of grassland rodents. We used periodic matrices to estimate stochastic elasticity for each vital rate and then tested for correlations with a standardized coefficient of variation for each rate. 3. We calculated stochastic elasticities for individual months to test for an association between increased reproduction and the influence of reproduction, relative to survival, on the population growth rate. 4. All species showed some evidence of buffering. The elasticity of vital rates of Peromyscus leucopus (Rafinesque, 1818), Sigmodon hispidus Say & Ord, 1825 and Microtus ochrogaster (Wagner, 1842) was negatively related to vital rate CV. Elasticity and vital rate CV were negatively related in Peromyscus maniculatus (Wagner, 1845), but the relationship was not statistically significant. Peromyscus leucopus and M. ochrogaster showed plasticity in vital rates; reproduction was higher following months where elasticity for reproduction exceeded that of survival. 5. Our results suggest that buffering is common in species with fast life histories; however, some populations that exhibit buffering are capable of responding to short-term variability in environmental conditions through reproductive plasticity.