Spatial anatomy of species survival: effects of predation and climate-driven environmental variability

The majority of survival analyses focus on temporal scales. Consequently, there is a limited understanding of how species survival varies over space and, ultimately, how spatial variability in the environment affects the temporal dynamics of species abundance. Using data from the Barents Sea, we study the spatiotemporal variability of the juvenile Atlantic cod (Gadus morhua) survival. We develop an index of spatial survival based on changes of juvenile cod distribution through their first winter of life (from age-0 to age-1) and study its variability in relation to biotic and abiotic factors. Over the 25 years analyzed (1980-2004), we found that, once the effect of passive drift due to dominant currents is accounted for, the area where age-0 cod survival was lowest coincided with the area of highest abundance of older cod. Within this critical region, the survival of age-0 cod was negatively affected by its own abundance, by that of older cod, and by bottom depth. Furthermore, during cold years, age-0 cod survival increased in the eastern and coldest portion of the examined area, which was typically avoided by older conspecifics. Based on these results we propose that within the examined area top-down mechanisms and predation-driven density dependence can strongly affect the spatial pattern of age-0 cod survival. Climate-related variables can also influence the spatial survival of age-0 cod by affecting their distribution and that of their predators. Results from these and similar studies, focusing on the spatial variability of survival rates, can be used to characterize species habitat quality of marine renewable resources.     

Assessing tiger population dynamics using photographic capture-recapture sampling

Although wide-ranging, elusive, large carnivore species, such as the tiger, are of scientific and conservation interest, rigorous inferences about their population dynamics are scarce because of methodological problems of sampling populations at the required spatial and temporal scales. We report the application of a rigorous, noninvasive method for assessing tiger population dynamics to test model-based predictions about population viability. We obtained photographic capture histories for 74 individual tigers during a nine-year study involving 5725 trap-nights of effort. These data were modeled under a likelihood-based, "robust design" capture-recapture analytic framework. We explicitly modeled and estimated ecological parameters such as time-specific abundance, density, survival, recruitment, temporary emigration, and transience, using models that incorporated effects of factors such as individual heterogeneity, trap-response, and time on probabilities of photo-capturing tigers. The model estimated a random temporary emigration parameter of gamma" = gamma' = 0.10 +/- 0.069 (values are estimated mean +/- SE). When scaled to an annual basis, tiger survival rates were estimated at S = 0.77 +/- 0.051, and the estimated probability that a newly caught animal was a transient was tau = 0.18 +/- 0.11. During the period when the sampled area was of constant size, the estimated population size N(t) varied from 17 +/- 1.7 to 31 +/- 2.1 tigers, with a geometric mean rate of annual population change estimated as lambda = 1.03 +/- 0.020, representing a 3% annual increase. The estimated recruitment of new animals, B(t), varied from 0 +/- 3.0 to 14 +/- 2.9 tigers. Population density estimates, D, ranged from 7.33 +/- 0.8 tigers/100 km2 to 21.73 +/- 1.7 tigers/100 km2 during the study. Thus, despite substantial annual losses and temporal variation in recruitment, the tiger density remained at relatively high levels in Nagarahole. Our results are consistent with the hypothesis that protected wild tiger populations can remain healthy despite heavy mortalities because of their inherently high reproductive potential. The ability to model the entire photographic capture history data set and incorporate reduced-parameter models led to estimates of mean annual population change that were sufficiently precise to be useful. This efficient, noninvasive sampling approach can be used to rigorously investigate the population dynamics of tigers and other elusive, rare, wide-ranging animal species in which individuals can be identified from photographs or other means.     

Testing sustainability by prospective and retrospective demographic analyses: evaluation for palm leaf harvest.

Harvesting nontimber forest products (NTFPs) is a major economic activity in tropical forests. As many NTFPs are overexploited, sustainability analyses are required to set harvest guidelines. Here we introduce and apply a new approach to evaluating sustainability, which combines prospective (elasticity) and retrospective (Life Table Response Experiments [LTRE]) demographic analyses of matrix population models. We relate the elasticity of vital rates (representing their importance for population growth rate, lamda) to their contribution to harvest-induced change in lamda ("LTRE contribution"). When high-elasticity vital rates have a low LTRE contribution, exploitation is potentially sustainable as negative effects for population growth are buffered. If the reverse is found, there is little scope for sustainability because crucial vital rates are affected. Our approach is less sensitive to chance fluctuations than the commonly used sustainability criterion of lamda = 1.0, as it does not depend on the absolute value of lamda. We applied this analysis to Geonoma deversa, a clustered forest understory palm. We studied three experimentally defoliated and control populations in a Bolivian rainforest during two years. Cutting all leaves of large ramets did not change mortality but strongly affected growth and reproduction. In spite of severe changes in some vital rates, population growth rate was not significantly reduced after defoliation. A literature review revealed that six other understory palms species responded very similarly to defoliation. The combination of LTRE contributions and elasticity analyses showed that low-elasticity vital rates were mainly responsible for the defoliation-induced change in lamda for Geonoma deversa. For two other understory palms (Astrocaryum mexicanum and Chamaedorea radicalis) new demographic analyses yielded very similar results. For Geonoma, the LTRE contribution-elasticity relation strongly changed when we mimicked harvest damage. Adding 5% mortality to defoliated palms caused stronger change in lamda, mainly due to changes in a high-elasticity vital rate (survival). Therefore, harvest practices that involve stem killing are clearly unsustainable. Our results show that commercial leaf cutting in Geonoma deversa is potentially sustainable, and that this is likely the case for understory palms in general. Our approach to evaluating harvest sustainability can be applied to other NTFPs.     

Spatial and temporal variability in somatic growth of green sea turtles (<Emphasis Type="Italic">Chelonia mydas</Emphasis>) resident in the Hawaiian Archipelago

The somatic growth dynamics of green turtles (Chelonia mydas) resident in five separate foraging grounds within the Hawaiian Archipelago were assessed using a robust non-parametric regression modelling approach. The foraging grounds range from coral reef habitats at the north-western end of the archipelago, to coastal habitats around the main islands at the south-eastern end of the archipelago. Pelagic juveniles recruit to these neritic foraging grounds from ca. 35 cm SCL or 5 kg (~6 years of age), but grow at foraging-ground-specific rates, which results in quite different size- and age-specific growth rate functions. Growth rates were estimated for the five populations as change in straight carapace length (cm SCL year^–1) and, for two of the populations, also as change in body mass (kg year^–1). Expected growth rates varied from ca. 0–2.5 cm SCL year^–1, depending on the foraging-ground population, which is indicative of slow growth and decades to sexual maturity, since expected size of first-time nesters is 80 cm SCL. The expected size-specific growth rate functions for four populations sampled in the south-eastern archipelago displayed a non-monotonic function, with an immature growth spurt at ca. 50–53 cm SCL (~18–23 kg) or ca. 13–19 years of age. The growth spurt for the Midway atoll population in the north-western archipelago occurs at a much larger size (ca. 65 cm SCL or 36 kg), because of slower immature growth rates that might be due to a limited food stock and cooler sea surface temperature. Expected age-at-maturity was estimated to be ca. 35–40 years for the four populations sampled at the south-eastern end of the archipelago, but it might well be >50 years for the Midway population. The Hawaiian stock comprises mainly the same mtDNA haplotype, with no differences in mtDNA stock composition between foraging-ground populations, so that the geographic variability in somatic growth rates within the archipelago is more likely due to local environmental factors rather than genetic factors. Significant temporal variability was also evident, with expected growth rates declining over the last 10–20 years, while green turtle abundance within the archipelago has increased significantly since the mid-1970s. This inverse relationship between somatic growth rates and population abundance suggests a density-dependent effect on somatic growth dynamics that has also been reported recently for a Caribbean green turtle stock. The Hawaiian green turtle stock is characterised by slow growth rates displaying significant spatial and temporal variation and an immature growth spurt. This is consistent with similar findings for a Great Barrier Reef green turtle stock that also comprises many foraging-ground populations spanning a wide geographic range.Communicated by P.W. Sammarco, Chauvin     

Synergistic influences of phase, density, and climatic variation on the dynamics of fluctuating populations

Although ecologists have long recognized that certain mammalian species exhibit high-amplitude, often multiannual, fluctuations in abundance, their causes have remained poorly understood and the subject of intense debate. A key contention has been the relative role of density-dependent and density-independent processes in governing population dynamics. We applied capture-mark-recapture analysis to 25 years of monthly trapping data from a fluctuating prairie vole Microtus ochrogaster population in Illinois, USA, to estimate realized population growth rates and associated vital rates (survival and recruitment) and modeled them as a function of vole density and density-independent climatic variation. We also tested for phase dependence and seasonality in the effects of the above processes. Variation in the realized population growth rate was best explained by phase-specific changes in vole density lagged by one month and mean monthly temperatures with no time lags. The underlying vital rates, survival and recruitment, were influenced by the additive and interactive effects of phase, vole density, and mean monthly temperatures. Our results are consistent with the observation that large-scale population fluctuations are characterized by phase-specific changes in demographic and physiological characteristics. Our findings also support the growing realization that the interaction between climatic variables and density-dependent factors may be a widespread phenomenon, and they suggest that the direction and magnitude of such interactive effects may be phase specific. We conclude that density-dependent and density-independent climatic variables work in tandem during each phase of density fluctuations to drive the dynamics of fluctuating populations.     

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.     

Factors contributing to process variance in annual survival of female greater sage-grouse in Montana.

Populations of Greater Sage-Grouse (Centrocercus urophasianus) have declined by 69-99% from historic levels, and information on population dynamics of these birds at a landscape scale is essential to informed management. We examined the relationships between hen survival and a suite of landscape-scale habitat and environmental conditions. We radio-marked 237 female Sage-Grouse and measured 426 vegetation plots during 2001-2004 at four sites in a 3200-km2 landscape in north-central Montana, USA. We used program MARK to model monthly survival rates for 11 seasonal intervals. There was strong support for the best-approximating model (AICc weight = 0.810), which indicated that (1) hen survival varied by season within years and by year within seasons, (2) nesting hens had higher nesting-season survival than non-nesting hens, and (3) individuals at one site had lower hunting-season survival than at other sites. We observed considerable variation in hen survival. Process variation was 0.255, with an expected range of annual survival of 0.12 to 1.0. The ratio of process to total variation was 0.999, indicating that observed variation was real and not attributable to sampling variation. We observed a nearly fourfold difference in maximum and minimum annual survival, ranging from 0.962 +/- 0.024 (mean +/- SE) for nesting hens in 2001-2002 to 0.247 +/- 0.050) for non-nesters in 2003-2004. Low annual survival in 2003 resulted from the compounded effects of a West Nile virus outbreak in August and a severe winter in 2003-2004. Increased hen mortality associated with severe winter weather contrasts with prior beliefs that Sage-Grouse populations are typically unaffected by winter weather conditions and underscores the importance of protecting winter sagebrush (Artemisia spp.) habitats.     

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

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

Density dependence vs. independence, and irregular population dynamics of a swallow-wort fruit fly

Although most long-term studies of consumer-resource (e.g., predator-prey) interactions select species showing cyclic population dynamics, strong consumer-resource interactions can also produce irregular, noncyclic dynamics. Here, we present a case in which a seed predator, the tephritid fruit fly Euphranta connexa, shows fluctuations in density of more than two orders of magnitude over a 22-year period. To explain these fluctuations, we analyzed a stage-specific data set to quantify the density-dependent and density-independent components of larval survivorship and realized fecundity. Both larval survivorship and realized fecundity were strongly density dependent. Larval survivorship dropped from 0.62 at low larval density to 0.081 at high larval density, whereas fecundity dropped from 84.3 to 0.32 eggs per individual, more than a 100-fold decrease. We divided density-independent variation in E. connexa population dynamics into components for variability in (1) larval survivorship, (2) realized fecundity, and (3) annual fruit abundance. Of these components, 96% of the density-independent variance in per capita population growth rates was caused by fluctuations in fruit abundance. This highlights the importance of the strong consumer-resource interactions in driving fluctuations in E. connexa abundance. It also demonstrates that E. connexa dynamics are remarkably simple, and aside from the 4% of unexplained variance in per capita population growth rates, our understanding of E. connexa dynamics is remarkably complete.     

Estimating survival rates with time series of standing age-structure data

It has long been recognized that age-structure data contain useful information for assessing the status and dynamics of wildlife populations. For example, age-specific survival rates can be estimated with just a single sample from the age distribution of a stable, stationary population. For a population that is not stable, age-specific survival rates can be estimated using techniques such as inverse methods that combine time series of age-structure data with other demographic data. However, estimation of survival rates using these methods typically requires numerical optimization, a relatively long time series of data, and smoothing or other constraints to provide useful estimates. We developed general models for possibly unstable populations that combine time series of age-structure data with other demographic data to provide explicit maximum likelihood estimators of age-specific survival rates with as few as two years of data. As an example, we applied these methods to estimate survival rates for female bison (Bison bison) in Yellowstone National Park, USA. This approach provides a simple tool for monitoring survival rates based on age-structure data.     

Spectral analysis unmasks synchronous and compensatory dynamics in plankton communities

Community biomass is often less variable than the biomasses of populations within the community, yet attempts to implicate compensatory dynamics between populations as a cause of this relationship often fail. In part, this may be due to the lack of appropriate metrics for variability, but there is also great potential for large-scale processes such as seasonality or longer-term environmental change to obscure important dynamics at other temporal scales. In this study, we apply a scale-resolving method to long-term plankton data, to identify the specific temporal scales at which community-level variability is influenced by synchrony or compensatory dynamics at the population level. We show that variability at both the population and community level is influenced strongly by a few distinct temporal scales: in phytoplankton, ciliate, rotifer, and crustacean communities, synchronous dynamics are predominant at most temporal scales. However, in phytoplankton and crustacean communities, compensatory dynamics occur at a sub-annual scale (and at the annual scale in crustaceans) leading to substantial reductions in community-level variability. Aggregate measures of population and community variability do not detect compensatory dynamics in these communities; thus, resolving their scale dependence unmasks dynamics that are important for community stability in this system. The methods and results presented herein will ultimately lead to a better understanding of how stability is achieved in communities.     

The interaction between reproductive cost and individual quality is mediated by oceanic conditions in a long-lived bird

Environmental variability, costs of reproduction, and heterogeneity in individual quality are three important sources of the temporal and interindividual variations in vital rates of wild populations. Based on an 18-year monitoring of an endangered, recently described, long-lived seabird, Monteiro's Storm-Petrel (Oceanodroma monteiroi), we designed multistate survival models to separate the effects of the reproductive cost (breeders vs. nonbreeders) and individual quality (successful vs. unsuccessful breeders) in relation to temporally variable demographic and oceanographic properties. The analysis revealed a gradient of individual quality from nonbreeders, to unsuccessful breeders, to successful breeders. The survival rates of unsuccessful breeders (0.90 +/- 0.023, mean +/- SE) tended to decrease in years of high average breeding success and were more sensitive to oceanographic variation than those of both (high-quality) successful breeders (0.97 +/- 0.015) and (low-quality) nonbreeders (0.83 +/- 0.028). Overall, our results indicate that reproductive costs act on individuals of intermediate quality and are mediated by environmental harshness.     

Effective Population Size in Red-Cockaded Woodpeckers: Population and Model Differences

Loss of genetic variability in isolated populations is an important issue for conservation biology. Most studies involve only a single population of a given species and a single method of estimating rate of loss. Here we present analyses for three different Red-cockaded Woodpecker ( Picoides borealis ) populations from different geographic regions. We compare two different models for estimating the expected rate of loss of genetic variability, and test their sensitivity to model parameters. We found that the simpler model (Reed et al. 1988) consistently estimated a greater rate of loss of genetic variability from a population than did the Emigh and Pollak (1979) model. The ratio of effective population size (which describes the expected rate of loss of genetic variability) to breeder population size varied widely among Red-cockaded Woodpecker populations due to geographic variation in demography. For this species, estimates of effective size were extremely sensitive to survival parameters, but not to the probability of breeding or reproductive success. Sensitivity was sufficient that error in estimating survival rates in the field could easily mask true population differences in effective size. Our results indicate that accurate and precise demographic data are prerequisites to determining effective population size for this species using genetic models, and that a single estimate of rate of loss of genetic variability is not valid across populations.     

Estimates of sex- and age-class-specific survival probabilities for a southern Great Barrier Reef green sea turtle population

Sex- and age-class-specific survival probabilities of a southern Great Barrier Reef green sea turtle population were estimated using a capture–mark–recapture (CMR) study and a Cormack–Jolly–Seber (CJS) modelling approach. The CMR history profiles for 954 individual turtles tagged over a 9-year period (1984–1992) were classified into three age classes (adult, subadult, juvenile) based on somatic growth and reproductive traits. Reduced-parameter CJS models, accounting for constant survival and time-specific recapture, fitted best for all age classes. There were no significant sex-specific differences in either survival or recapture probabilities for any age class. Mean annual adult survival was estimated at 0.9482 (95% CI: 0.92–0.98) and was significantly higher than survival for either subadults or juveniles. Mean annual subadult survival was 0.8474 (95% CI: 0.79–0.91), which was not significantly different from mean annual juvenile survival estimated at 0.8804 (95% CI: 0.84–0.93). The time-specific adult recapture probabilities were a function of sampling effort but this was not the case for either juveniles or subadults. The sampling effort effect was accounted for explicitly in the estimation of adult survival and recapture probabilities. These are the first comprehensive sex- and age-class-specific survival and recapture probability estimates for a green sea turtle population derived from a long-term CMR program.Communicated by M.S. Johnson, Crawley     

From metamorphosis to maturity in complex life cycles: equal performance of different juvenile life history pathways

Performance in one stage of a complex life cycle may affect performance in the subsequent stage. Animals that start a new stage at a smaller size than conspecifics may either always remain smaller or they may be able to "catch up" through plasticity, usually elevated growth rates. We study how size at and date of metamorphosis affected subsequent performance in the terrestrial juvenile stage and lifetime fitness of spadefoot toads (Pelobates fuscus). We analyzed capture-recapture data of > 3000 individuals sampled during nine years with mark-recapture models to estimate first-year juvenile survival probabilities and age-specific first-time breeding probabilities of toads, followed by model selection to assess whether these probabilities were correlated with size at and date of metamorphosis. Males attained maturity after two years, whereas females reached maturity 2-4 years after metamorphosis. Age at maturity was weakly correlated with metamorphic traits. In both sexes, first-year juvenile survival depended positively on date of metamorphosis and, in males, also negatively on size at metamorphosis. In males, toads that metamorphosed early at a small size had the highest probability to reach maturity. However, because very few toadlets metamorphosed early, the vast majority of male metamorphs had a very similar probability to reach maturity. A matrix projection model constructed for females showed that different juvenile life history pathways resulted in similar lifetime fitness. We found that the effects of date of and size at metamorphosis on different juvenile traits cancelled each other out such that toads that were small or large at metamorphosis had equal performance. Because the costs and benefits of juvenile life history pathways may also depend on population fluctuations, ample phenotypic variation in life history traits may be maintained.     

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.     

Estimating survival rates of uncatchable animals: the myth of high juvenile mortality in reptiles

Survival rates of juvenile reptiles are critical population parameters but are difficult to obtain through mark-recapture programs because these small, secretive animals are rarely caught. This scarcity has encouraged speculation that survival rates of juveniles are very low, and we test this prediction by estimating juvenile survival rates indirectly. A simple mathematical model calculates the annual juvenile survival rate needed to maintain a stable population size, using published data on adult survival rates, reproductive output, and ages at maturity in 109 reptile populations encompassing 57 species. Counter to prediction, estimated juvenile survival rates were relatively high (on average, only about 13% less than those of conspecific adults) and highly correlated with adult survival rates. Overall, survival rates during both juvenile and adult life were higher in turtles than in snakes, and higher in snakes than in lizards. As predicted from life history theory, rates of juvenile survival were higher in species that produce large offspring, and higher in viviparous squamates than in oviparous species. Our analyses challenge the widely held belief that juvenile reptiles have low rates of annual survival and suggest instead that sampling problems and the elusive biology of juvenile reptiles have misled researchers in this respect.     

Seed availability and insect herbivory limit recruitment and adult density of native tall thistle

Understanding spatial and temporal variation in factors influencing plant regeneration is critical to predicting plant population growth. We experimentally evaluated seed limitation, insect herbivory, and their interaction in the regeneration and density of tall thistle (Cirsium altissimum) across a topographic ecosystem productivity gradient in tallgrass prairie over two years. On ridges and in valleys, we used a factorial experiment manipulating seed availability and insect herbivory to quantify effects of: seed input on seedling density, insect herbivory on juvenile density, and cumulative impacts of both seed input and herbivory on reproductive adult density. Seed addition increased seedling densities at three of five sites in 2006 and all five sites in 2007. Insect herbivory reduced seedling survival across all sites in both years, as well as rosette survival from the previous year's seedlings. In both years, insecticide treatment of seed addition plots led to greater adult tall thistle densities in the following year, reflecting the increase in juvenile thistle densities in the experimental year. Seedling survival was not density dependent. Our analytical projection model predicts a significant long-term increase in adult densities from seed input, with a greater increase under experimentally reduced insect herbivory. While plant community biomass and water stress varied significantly between ridges and valleys, the effects of seed addition and insect herbivory did not vary with gradient position. These results support conceptual models that predict seedling and adult densities of short-lived monocarpic perennial plants should be seed limited. Further, the experiment demonstrates that even at high juvenile plant densities, at which density dependence potentially could have overridden herbivore effects on plant survival, insect herbivory strongly affected juvenile thistle performance and adult densities of this native prairie species.