The natural history of luck: A synthesis study of structured population models.

Chance pervades life. In turn, life histories are described by probabilities (e.g. survival and breeding) and averages across individuals (e.g. mean growth rate and age at maturity). In this study, we explored patterns of luck in lifetime outcomes by analysing structured population models for a wide array of plant and animal species. We calculated four response variables: variance and skewness in both lifespan and lifetime reproductive output (LRO), and partitioned them into contributions from different forms of luck. We examined relationships among response variables and a variety of life history traits. We found that variance in lifespan and variance in LRO were positively correlated across taxa, but that variance and skewness were negatively correlated for both lifespan and LRO. The most important life history trait was longevity, which shaped variance and skew in LRO through its effects on variance in lifespan. We found that luck in survival, growth, and fecundity all contributed to variance in LRO, but skew in LRO was overwhelmingly due to survival luck. Rapidly growing populations have larger variances in LRO and lifespan than shrinking populations. Our results indicate that luck-induced genetic drift may be most severe in recovering populations of species with long mature lifespan and high iteroparity.

Larval and Adult Traits Coevolve in Response to Asymmetric Coastal Currents to Shape Marine Dispersal Kernels.

AbstractDispersal emerges as an outcome of organismal traits and external forcings. However, it remains unclear how the emergent dispersal kernel evolves as a by-product of selection on the underlying traits. This question is particularly compelling in coastal marine systems, where dispersal is tied to development and reproduction and where directional currents bias larval dispersal downstream, causing selection for retention. We modeled the dynamics of a metapopulation along a finite coastline using an integral projection model and adaptive dynamics to understand how asymmetric coastal currents influence the evolution of larval (pelagic larval duration) and adult (spawning frequency) life history traits, which indirectly shape the evolution of marine dispersal kernels. Selection induced by alongshore currents favors the release of larvae over multiple time periods, allowing long pelagic larval durations and long-distance dispersal to be maintained in marine life cycles in situations where they were previously predicted to be selected against. Two evolutionarily stable strategies emerged: one with a long pelagic larval duration and many spawning events, resulting in a dispersal kernel with a larger mean and variance, and another with a short pelagic larval duration and few spawning events, resulting in a dispersal kernel with a smaller mean and variance. Our theory shows how coastal ocean flows are important agents of selection that can generate multiple, often co-occurring evolutionary outcomes for marine life history traits that affect dispersal.

Severe fire regimes decrease resilience of ectothermic populations.

Understanding populations' responses to environmental change is crucial for mitigating human-induced disturbances. Here, we test hypotheses regarding how three essential components of demographic resilience (resistance, compensation and recovery) co-vary along the distinct life histories of three lizard species exposed to variable, prescribed fire regimes. Using a Bayesian hierarchical framework, we estimate vital rates (survival, growth and reproduction) with 14 years of monthly individual-level data and mark-recapture models to parameterize stochastic integral projection models from five sites in Brazilian savannas, each historically subjected to different fire regimes. With these models, we investigate how weather, microclimate and ecophysiological traits of each species influence their vital rates, emergent life history traits and demographic resilience components in varying fire regimes. Overall, weather and microclimate are better predictors of the species' vital rates, rather than their ecophysiological traits. Our findings reveal that severe fire regimes increase populations' resistance but decrease compensation or recovery abilities. Instead, populations have higher compensatory and recovery abilities at intermediate degrees of fire severity. Additionally, we identify generation time and reproductive output as predictors of resilience trends across fire regimes and climate. Our analyses demonstrate that the probability and quantity of monthly reproduction are the proximal drivers of demographic resilience across the three species. Our findings suggest that populations surpass a tipping point in severe fire regimes and achieve an alternative stable state to persist. Thus, higher heterogeneity in fire regimes can increase the reproductive aspects and resilience of different populations and avoid high-severity regimes that homogenize the environment. Despite being more resistant, species with long generation times and low reproductive output take longer to recover and cannot compensate as much as species with faster paces of life. We emphasize how reproductive constraints, such as viviparity and fixed clutch sizes, impact the ability of ectothermic populations to benefit and recover from disturbances, underscoring their relevance in conservation assessments.

Why we cannot always expect life history strategies to directly inform on sensitivity to environmental change.

Variation in life history traits in animals and plants can often be structured along major axes of life history strategies. The position of a species along these axes can inform on their sensitivity to environmental change. For example, species with slow life histories are found to be less sensitive in their long-term population responses to environmental change than species with fast life histories. This provides a tantalizing link between sets of traits and population responses to change, contained in a highly generalizable theoretical framework. Life history strategies are assumed to reflect the outcome of life history tradeoffs that, by their very nature, act at the individual level. Examples include the tradeoff between current and future reproductive success, and allocating energy into growth versus reproduction. But the importance of such tradeoffs in structuring population-level responses to environmental change remains understudied. We aim to increase our understanding of the link between individual-level life history tradeoffs and the structuring of life history strategies across species, as well as the underlying links to population responses to environmental change. We find that the classical association between lifehistory strategies and population responses to environmental change breaks down when accounting for individual-level tradeoffs and energy allocation. Therefore, projecting population responses to environmental change should not be inferred based only on a limited set of species traits. We summarize our perspective and a way forward in a conceptual framework.

Compensatory responses of vital rates attenuate impacts of competition on population growth and promote coexistence.

Competition is among the most important factors regulating plant population and community dynamics, but we know little about how different vital rates respond to competition and jointly determine population growth and species coexistence. We conducted a field experiment and parameterised integral projection models to model the population growth of 14 herbaceous plant species in the absence and presence of neighbours across an elevation gradient (284 interspecific pairs). We found that suppressed individual growth and seedling establishment contributed the most to competition-induced declines in population growth, although vital rate contributions varied greatly between species and with elevation. In contrast, size-specific survival and flowering probability and seed production were frequently enhanced under competition. These compensatory vital rate responses were nearly ubiquitous (occurred in 92% of species pairs) and significantly reduced niche overlap and stabilised coexistence. Our study highlights the importance of demographic processes for regulating population and community dynamics, which has often been neglected by classic coexistence theories.

Evaluating the importance of individual heterogeneity in reproduction to Weddell seal population dynamics using integral projection models.

Identifying and accounting for unobserved individual heterogeneity in vital rates in demographic models is important for estimating population-level vital rates and identifying diverse life-history strategies, but much less is known about how this individual heterogeneity influences population dynamics. We aimed to understand how the distribution of individual heterogeneity in reproductive and survival rates influenced population dynamics using vital rates from a Weddell seal population by altering the distribution of individual heterogeneity in reproduction, which also altered the distribution of individual survival rates through the incorporation of our estimate of the correlation between the two rates and assessing resulting changes in population growth. We constructed an integral projection model (IPM) structured by age and reproductive state using estimates of vital rates for a long-lived mammal that has recently been shown to exhibit large individual heterogeneity in reproduction. Using output from the IPM, we evaluated how population dynamics changed with different underlying distributions of unobserved individual heterogeneity in reproduction. Results indicate that the changes to the underlying distribution of individual heterogeneity in reproduction cause very small changes in the population growth rate and other population metrics. The largest difference in the estimated population growth rate resulting from changes to the underlying distribution of individual heterogeneity was less than 1%. Our work highlights the differing importance of individual heterogeneity at the population level compared to the individual level. Although individual heterogeneity in reproduction may result in large differences in the lifetime fitness of individuals, changing the proportion of above- or below-average breeders in the population results in much smaller differences in annual population growth rate. For a long-lived mammal with stable and high adult-survival that gives birth to a single offspring, individual heterogeneity in reproduction has a limited effect on population dynamics. We posit that the limited effect of individual heterogeneity on population dynamics may be due to canalization of life-history traits.

Modeling the effects of ration on individual growth of Oncorhynchus tshawytscha under controlled conditions.

Fed aquaculture is one of the fastest-growing and most valuable food production industries in the world. The efficiency with which farmed fish convert feed into biomass influences both environmental impact and economic revenue. Salmonid species, such as king salmon (Oncorhynchus tshawytscha), exhibit high levels of plasticity in vital rates such as feed intake and growth rates. Accurate estimations of individual variability in vital rates are important for production management. The use of mean trait values to evaluate feeding and growth performance can mask individual-level differences that potentially contribute to inefficiencies. Here, the authors apply a cohort integral projection model (IPM) framework to investigate individual variation in growth performance of 1625 individually tagged king salmon fed one of three distinct rations of 60%, 80%, and 100% satiation and tracked over a duration of 276 days. To capture the observed sigmoidal growth of individuals, they compared a nonlinear mixed-effects (logistic) model to a linear model used within the IPM framework. Ration significantly influenced several aspects of growth, both at the individual and at the cohort level. Mean final body mass and mean growth rate increased with ration; however, variance in body mass and feed intake also increased significantly over time. Trends in mean body mass and individual body mass variation were captured by both logistic and linear models, suggesting the linear model to be suitable for use in the IPM. The authors also observed that higher rations resulted in a decreasing proportion of individuals reaching the cohort's mean body mass or larger by the end of the experiment. This suggests that, in the present experiment, feeding to satiation did not produce the desired effects of efficient, fast, and uniform growth in juvenile king salmon. Although monitoring individuals through time is challenging in commercial aquaculture settings, recent technological advances combined with an IPM approach could provide new scope for tracking growth performance in experimental and farmed populations. Using the IPM framework might allow the exploration of other size-dependent processes affecting vital rate functions, such as competition and mortality.

The role of demographic compensation in stabilising marginal tree populations in North America.

Demographic compensation-the opposing responses of vital rates along environmental gradients-potentially delays anticipated species' range contraction under climate change, but no consensus exists on its actual contribution. We calculated population growth rate (λ) and demographic compensation across the distributional ranges of 81 North American tree species and examined their responses to simulated warming and tree competition. We found that 43% of species showed stable population size at both northern and southern edges. Demographic compensation was detected in 25 species, yet 15 of them still showed a potential retraction from southern edges, indicating that compensation alone cannot maintain range stability. Simulated climatic warming caused larger decreases in λ for most species and weakened the effectiveness of demographic compensation in stabilising ranges. These findings suggest that climate stress may surpass the limited capacity of demographic compensation and pose a threat to the viability of North American tree populations.

Climate-driven, but dynamic and complex? A reconciliation of competing hypotheses for species' distributions.

Estimates of the percentage of species "committed to extinction" by climate change range from 15% to 37%. The question is whether factors other than climate need to be included in models predicting species' range change. We created demographic range models that include climate vs. climate-plus-competition, evaluating their influence on the geographic distribution of Pinus edulis, a pine endemic to the semiarid southwestern U.S. Analyses of data on 23,426 trees in 1941 forest inventory plots support the inclusion of competition in range models. However, climate and competition together only partially explain this species' distribution. Instead, the evidence suggests that climate affects other range-limiting processes, including landscape-scale, spatial processes such as disturbances and antagonistic biotic interactions. Complex effects of climate on species distributions-through indirect effects, interactions, and feedbacks-are likely to cause sudden changes in abundance and distribution that are not predictable from a climate-only perspective.

Monitoring Demography of Resurrected Populations of Locally Extinct and Extant Species to Investigate Drivers of Species Loss.

AbstractExtinctions are predicted to rise by an order of magnitude over the next century. Although contemporary documented extinctions are uncommon, local extirpations likely provide hints about global extinction risks. Comparing responses to global change of locally extinct versus extant species pairs in a phylogenetic framework could highlight why certain species are more vulnerable to extinction than others and which anthropogenic changes are most relevant to their decline. As anthropogenic changes likely interact to affect population declines, demographic studies partitioning the effects of multifactorial stressors are needed but remain rare. I examine demographic responses to nitrogen addition and deer herbivory, two major drivers of species losses in grasslands, in experimental reintroductions of 14 locally extinct and extant confamilial native plants from Michigan prairies. Nitrogen consistently reduces survival, especially in locally extinct species, and growth of locally extinct species benefits less from nitrogen than growth of extant species. Nitrogen reduces population growth rates, largely via reductions in survival. Deer herbivory, meanwhile, had inconsistent effects on vital rates among species and did not affect population growth. Nitrogen and herbivory rarely interacted to affect vital rates. These results link community-level patterns of species loss under nitrogen addition to the population-level processes underlying those losses.

Pest management in future climates: Warming reduces physical weed management effectiveness.

Climate change alters many aspects of weed performance and may also alter the effectiveness of management practices to control pests. Despite this concern, entire categories of widely used management practices, such as physical control, remain understudied in this context. We conducted a field experiment growing the invasive pest musk thistle (Carduus nutans) at ambient and experimentally elevated temperatures. We tested mowing management strategies that varied in the timing of a single mowing event relative to thistles' stem elongation phenology and compared these with an unmowed control. Results from this experiment informed demographic models to project population growth rates for different warming/mowing scenarios. Compared to plants grown under ambient conditions, warmed thistles were more likely to survive the same mowing treatment, flowered earlier in the season, grew to taller heights, and produced more flowering capitula. Proportional reductions in plant height and capitulum production caused by mowing were smaller under warming. Warming did not change the relative ranking of mowing treatments; mowing late in the growing season (2 weeks after individuals first reached a height of 40 cm) was most effective at ambient temperatures and under warming. Warming caused significant increases in projected local population growth rate for all mowing treatments. For invasive musk thistle, warmed individuals outperformed individuals grown at ambient temperatures across all the mowing treatments we considered. Our results suggest that to achieve outcomes comparable to those attainable at today's temperatures, farmers will need to apply supplemental management, possibly including additional mowing effort or alternative practices such as chemical control. We recommend that scientists test management practices under experimental warming, where possible, and that managers monitor ongoing management to identify changes in effectiveness. Information about changes in managed weeds' mortality, fecundity, and phenology can then be used to make informed decisions in future climates.

DEVELOPING INTEGRAL PROJECTION MODELS FOR ECOTOXICOLOGY.

In many ecosystems, especially aquatic ecosystems, size plays a critical role in the factors that determine an individual's ability to survive and reproduce. In aquatic ecotoxicology, size informs both realized and potential acute and chronic effects of chemical exposure. This paper demonstrates how chemical and nonchemical effects on growth, survival, and reproduction can be linked to population-level dynamics using size-structured integral projection models (IPM). The modeling approach was developed with the goals and constraints of ecological risk assessors in mind, who are tasked with estimating the effects of chemical exposures to wildlife populations in a data-limited environment. The included case study is a collection of daily time-step IPMs parameterized for the life history and annual cycle of fathead minnows (Pimephales promelas), which motivated the development of modeling techniques for seasonal, iteroparous reproduction, density dependent growth effects, and size-dependent over-winter survival. The effects of a time-variable annual chemical exposure were interpreted using a toxicokinetic-toxicodynamic model for acute survival and sub-lethal growth effects model for chronic effects and incorporated into the IPMs. This paper presents a first application of integral projection models to ecotoxicology. Our research demonstrates that size-structured IPMs provide a promising, flexible, framework for synthesizing ecotoxicologically relevant data and theory to explore the effects of chemical and nonchemical stressors and the resulting impacts on exposed populations.

Quantifying the statistical power of monitoring programs for marine protected areas.

Marine Protected Areas (MPAs) are increasingly established globally as a spatial management tool to aid in conservation and fisheries management objectives. Assessing whether MPAs are having the desired effects on populations requires effective monitoring programs. A cornerstone of an effective monitoring program is an assessment of the statistical power of sampling designs to detect changes when they occur. We present a novel approach to power assessment that combines spatial point process models, integral projection models (IPMs) and sampling simulations to assess the power of different sample designs across a network of MPAs. We focus on the use of remotely operated vehicle (ROV) video cameras as the sampling method, though the results could be extended to other sampling methods. We use empirical data from baseline surveys of an example indicator fish species across three MPAs in California, USA as a case study. Spatial models simulated time series of spatial distributions across sites that accounted for the effects of environmental covariates, while IPMs simulated expected trends over time in abundances and sizes of fish. We tested the power of different levels of sampling effort (i.e., the number of 500-m ROV transects) and temporal replication (every 1-3 yr) to detect expected post-MPA changes in fish abundance and biomass. We found that changes in biomass are detectable earlier than changes in abundance. We also found that detectability of MPA effects was higher in sites with higher initial densities. Increasing the sampling effort had a greater effect than increasing sampling frequency on the time taken to achieve high power. High power was best achieved by combining data from multiple sites. Our approach provides a powerful tool to explore the interaction between sampling effort, spatial distributions, population dynamics, and metrics for detecting change in previously fished populations.

Latitudinal gradients in population growth do not reflect demographic responses to climate.

Spatial gradients in population growth, such as across latitudinal or elevational gradients, are often assumed to primarily be driven by variation in climate, and are frequently used to infer species' responses to climate change. Here, we use a novel demographic, mixed-model approach to dissect the contributions of climate variables vs. other latitudinal or local site effects on spatiotemporal variation in population performance in three perennial bunchgrasses. For all three species, we find that performance of local populations decreases with warmer and drier conditions, despite latitudinal trends of decreasing population growth toward the cooler and wetter northern portion of each species' range. Thus, latitudinal gradients in performance are not predictive of either local or species-wide responses to climate. This pattern could be common, as many environmental drivers, such as habitat quality or species' interactions, are likely to vary with latitude or elevation, and thus influence or oppose climate responses.

The projected degradation of subtropical coral assemblages by recurrent thermal stress.

Subtropical coral assemblages are threatened by similar extreme thermal stress events to their tropical counterparts. Yet, the mid- and long-term thermal stress responses of corals in subtropical environments remain largely unquantified, limiting our capacity to predict their future viability. The annual survival, growth and recruitment of 311 individual corals within the Solitary Islands Marine Park (Australia) was recorded over a 3-year period (2016-2018), including the 2015/2016 thermal stress event. These data were used to parameterise integral projection models quantifying the effect of thermal stress within a subtropical coral assemblage. Stochastic simulations were also applied to evaluate the implications of recurrent thermal stress scenarios predicted by four different Representative Concentration Pathways. We report differential shifts in population growth rates (λ) among coral populations during both stress and non-stress periods, confirming contrasting bleaching responses among taxa. However, even during non-stress periods, the observed dynamics for all taxa were unable to maintain current community composition, highlighting the need for external recruitment sources to support the community structure. Across all coral taxa, projected stochastic growth rates (λs ) were found to be lowest under higher emissions scenarios. Correspondingly, predicted increases in recurrent thermal stress regimes may accelerate the loss of coral coverage, species diversity and structural complexity within subtropical regions. We suggest that these trends are primarily due to the susceptibility of subtropical specialists and endemic species, such as Pocillopora aliciae, to thermal stress. Similarly, the viability of many tropical coral populations at higher latitudes is highly dependent on the persistence of up-current tropical systems. As such, the inherent dynamics of subtropical coral populations appear unable to support their future persistence under unprecedented thermal disturbance scenarios.

Different population trajectories of two reef-building corals with similar life-history traits.

Increases in the frequency and intensity of acute and chronic disturbances are causing declines of coral reefs world-wide. Although quantifying the responses of corals to acute disturbances is well documented, detecting subtle responses of coral populations to chronic disturbances is less common, but can also result in altered population and community structures. We investigated the population dynamics of two key reef-building Merulinid coral species, Dipsastraea favus and Platygyra lamellina, with similar life-history traits, in the Gulf of Eilat and Aqaba, Red Sea from 2015 to 2018, to assess potential differences in their population trajectories. Demographic processes, which included rates of survival, growth, reproduction and recruitment were used to parametrize integral projection models and estimate population growth rates and the likely population trajectories of both coral species. The survival and reproduction rates of both D. favus and P. lamellina were positively related to coral colony size, and elasticity analyses showed that large colonies most influenced population dynamics. Although both species have similar life-history traits and growth morphologies and are generally regarded as 'stress-tolerant', the populations showed contrasting trajectories-D. favus appears to be increasing whereas P. lamellina appears to be decreasing. As many corals have long-life expectancies, the process of local and regional decline might be subtle and slow. Ecological assessments based on total living coral coverage, morphological groups or functional traits might overlook subtle, species-specific trends. However, demographic approaches capable of detecting subtle species-specific population changes can augment ecological studies and provide valuable early warning signs of decline before major coral loss becomes evident.

Demographic compensation does not rescue populations at a trailing range edge.

Species' geographic ranges and climatic niches are likely to be increasingly mismatched due to rapid climate change. If a species' range and niche are out of equilibrium, then population performance should decrease from high-latitude "leading" range edges, where populations are expanding into recently ameliorated habitats, to low-latitude "trailing" range edges, where populations are contracting from newly unsuitable areas. Demographic compensation is a phenomenon whereby declines in some vital rates are offset by increases in others across time or space. In theory, demographic compensation could increase the range of environments over which populations can succeed and forestall range contraction at trailing edges. An outstanding question is whether range limits and range contractions reflect inadequate demographic compensation across environmental gradients, causing population declines at range edges. We collected demographic data from 32 populations of the scarlet monkeyflower (Erythranthe cardinalis) spanning 11° of latitude in western North America and used integral projection models to evaluate population dynamics and assess demographic compensation across the species' range. During the 5-y study period, which included multiple years of severe drought and warming, population growth rates decreased from north to south, consistent with leading-trailing dynamics. Southern populations at the trailing range edge declined due to reduced survival, growth, and recruitment, despite compensatory increases in reproduction and faster life-history characteristics. These results suggest that demographic compensation may only delay population collapse without the return of more favorable conditions or the contribution of other buffering mechanisms such as evolutionary rescue.

Climate change disrupts local adaptation and favours upslope migration.

Contemporary climate change is proceeding at an unprecedented rate. The question remains whether populations adapted to historical conditions can persist under rapid environmental change. We tested whether climate change will disrupt local adaptation and reduce population growth rates using the perennial plant Boechera stricta (Brassicaceae). In a large-scale field experiment conducted over five years, we exposed > 106 000 transplants to historical, current, or future climates and quantified fitness components. Low-elevation populations outperformed local populations under simulated climate change (snow removal) across all five experimental gardens. Local maladaptation also emerged in control treatments, but it was less pronounced than under snow removal. We recovered local adaptation under snow addition treatments, which reflect historical conditions. Our results revealed that low elevation populations risk rapid decline, whereas upslope migration could enable population persistence and expansion at higher elevation locales. Local adaptation to historical conditions could increase vulnerability to climate change, even for geographically widespread species.

Deer Indirectly Alter the Reproductive Strategy and Operational Sex Ratio of an Unpalatable Forest Perennial.

Environmental conditions impose restrictions and costs on reproduction. Multiple reproductive options exist when increased reproductive costs drive plant populations toward alternative reproductive strategies. Using 4 years of demographic data across a deer impact gradient, where deer alter the abiotic environment, we parameterize a size-dependent integral projection model for a sexually labile and unpalatable forest perennial to investigate the demographic processes driving differentiation in the operational sex ratio (OSR) of local populations. In addition to a relative increase in asexual reproduction, our results illustrate that nontrophic indirect effects by overabundant deer on this perennial result in delayed female sex expression to unsustainably large plant sizes and lead to more pronounced plant shrinkage following female sex expression, effectively increasing the cost of reproduction. Among plants of reproductive age, increased deer impact decreases the size-dependent probability of flowering and reduces reproductive consistency over time. This pattern in sex expression skews populations toward female-biased OSRs at low deer impact sites and male-biased OSRs at intermediate and high deer impact sites. While this shift toward a male-biased OSR may ameliorate pollen limitation, it also decreases the effective population size when coupled with increased asexual reproduction. The divergence of reproductive strategies and reduced lifetime fitness in response to indirect deer impacts illustrate the persistent long-term effects of overabundant herbivores on unpalatable understory perennials.

Demography of the giant monocarpic herb Rheum nobile in the Himalayas and the effect of disturbances by grazing.

BACKGROUND: Perennity of giant rosette species in combination with a single 'big bang' reproduction followed by death of the genetic individual is relatively rare among plants. Such long-lived monocarpic plants are usually slow growing and can be found in deserts, bogs or in alpine regions of the tropics or sub-tropics. Due to their longevity, monocarpic perennials risk losing everything before reproduction, which make them particularly susceptible to disturbances. Because of the inherent difficulties in assessing whether long-lived populations are growing or declining, usually neither their demography nor the consequences of increasing grazing pressure are known. METHODS: We used integral projection modelling (IPM) to measure the growth rate and passage time to flowering of Rheum nobile, a monocarpic perennial, and one of the most striking alpine plants from the high Himalayas. Rosettes which were no longer found due to disturbances or grazing by yaks were either treated as missing or as dead in two series of analysis, thereby simulating demography with and without the impact of grazing cattle. Data were collected from plants at 4500 m a.s.l. in Shangri-la County, Yunnan Province, south-west China. In four consecutive years (2011-2014) and in two populations, 372 and 369 individuals were measured, respectively, and size-dependent growth, survival and fecundity parameters were estimated. In addition, germination percentage, seedling survival and establishment probability were assessed. KEY RESULTS: The probability of survival, flowering and fecundity were strongly size dependent. Time to reach flowering size was 33.5 years [95 % confidence interval (CI) 21.9-43.3, stochastic estimate from pooled transitions and populations]. The stochastic population growth rate (λs) of Rheum nobile was 1.013 (95 % CI 1.010-1.017). When disturbance by grazing cattle (yaks) was accounted for in the model, λs dropped to values <1 (0.940, 95 % CI 0.938-0.943). CONCLUSION: We conclude that natural populations of this unique species are viable, but that conservation efforts should be made to minimize disturbances by grazing and to protect this slow-growing flagship plant from the high Himalayas.