New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
 Agricultural Sciences
 Anthropology
 Applied Biological Sciences
 Biochemistry
 Biophysics and Computational Biology
 Cell Biology
 Developmental Biology
 Ecology
 Environmental Sciences
 Evolution
 Genetics
 Immunology and Inflammation
 Medical Sciences
 Microbiology
 Neuroscience
 Pharmacology
 Physiology
 Plant Biology
 Population Biology
 Psychological and Cognitive Sciences
 Sustainability Science
 Systems Biology
Temporal variation, dormancy, and coexistence: A field test of the storage effect

Edited by Jane Lubchenco, Oregon State University, Corvallis, OR, and approved June 24, 1997 (received for review January 13, 1997)
Abstract
Theoretical models suggest that overlapping generations, in combination with a temporally fluctuating environment, may allow the persistence of competitors that otherwise would not coexist. Despite extensive theoretical development, this “storage effect” hypothesis has received little empirical attention. Herein I present the first explicit mathematical analysis of the contribution of the storage effect to the dynamics of competing natural populations. In Oneida Lake, NY, data collected over the past 30 years show a striking negative correlation between the watercolumn densities of two species of suspensionfeeding zooplankton, Daphnia galeata mendotae and Daphnia pulicaria. I have demonstrated competition between these two species and have shown that both possess longlived eggs that establish overlapping generations. Moreover, recruitment to this longlived stage varies annually, so that both daphnids have years in which they are favored (for recruitment) relative to their competitor. When the longterm population growth rates are calculated both with and without the effects of a variable environment, I show that D. galeata mendotae clearly cannot persist without the environmental variation and prolonged dormancy (i.e., storage effect) whereas D. pulicaria persists through consistently high per capita recruitment to the longlived stage.
The longterm persistence of competitors without competitive exclusion has led to numerous explanations for the maintenance of species diversity (1–3). Included among these hypotheses are those that posit resource partitioning and spatial or temporal separation of potential competitors at equilibrium (4–6). Alternatively, nonequilibrium theory suggests that a fluctuating environment may promote the coexistence of competitors (2, 7–9). Although variation alone cannot prevent the inevitable fixation of a single genotype or species (10, 11), Chesson and colleagues (2, 12–18) have proposed a way in which temporal recruitment fluctuations among species can lead to the stable coexistence of competitors. Their mechanism called the “storage effect,” may occur as long as each competitor can increase from low densities. One way in which this criterion is met is when each taxon has a longlived stage that is immune to the effects of competition and is able to carry the population through periods of poor recruitment. In an environment where recruitment is occasionally low, this longlived stage buffers competitors against exclusion in a way that is not possible in models with nonoverlapping generations.
In its original form, the storage effect results from three components: interspecific competition, a persistent longlived stage, and temporal variation in recruitment to this longlived stage. More recently, the theory has been expanded and can now be expressed as the outcome of three ingredients: speciesspecific responses to the environment, covariance between environment and competition, and subadditive growth rates (16–18). Comprehensive reviews of the general theory are available elsewhere (12–19). Herein, I focus only on the original theory because it provides quantitative methods that are applicable to field data.
The original prerequisites for maintenance of diversity via the storage effect (interspecific competition, a persistent longlived stage, and temporal variation in recruitment to this longlived stage) are found in a wide variety of both terrestrial and aquatic communities (19). In some systems (e.g., perennial plants, fish, and marine invertebrates), the overlapping generations are established by longlived adults and the competition is among juveniles (12, 20, 21). In others (e.g., annual plants, insects, zooplankton, and phytoplankton), competition exists among shortlived active individuals, and the longlived stage is a dormant seed, egg, or cyst (8, 9, 22, 23). In these systems, overlapping generations result from the repeated germination (or hatching) over a number of seasons of the seeds (or eggs) produced in any single year.
To date, nearly all attention to the storage effect hypothesis has been either strictly theoretical or, when applied to natural communities, merely anecdotal. One notable exception is the work of Pake and Venable (8, 9) on Sonoran Desert annuals. These investigators have demonstrated that all the necessary prerequisites of the model are met and suggested that coexistence of these competitors is mediated by environmental fluctuations and a betweenyear seed bank. However, they do not analyze quantitatively the contribution of the storage effect to the persistence of each population.
Herein I report results from, to my knowledge, the first explicit mathematical analysis of the storage effect’s contribution to the maintenance of species diversity. Twentyone years of recruitment data are presented for two species of Daphnia in Oneida Lake, NY. In addition, evidence for the duration of the longlived stage (diapausing egg) in the sediment egg bank and interspecific competition is summarized herein and presented in detail elsewhere (24). Because spatial structure of the eggs in the sediment may be important in the analyses of the longterm population dynamics, the egg bank is analyzed first as a homogeneous and then as a layered system. For Daphnia galeata mendotae, persistence in Oneida Lake, NY, is unambiguously dependent upon the storage effect, whereas its competitor, Daphnia pulicaria, persists through consistently high per capita recruitment.
STUDY SITE AND METHODS
In Oneida Lake, NY, zooplankton data collected over the past 30 years (25) show a striking negative correlation between the watercolumn densities of two zooplankton species, D. galeata mendotae and D. pulicaria [Spearman Rank r_{s} = −0.259; P < 0.001 (24)]. For at least the past three decades, these two species have alternated as the numerically dominant daphnid species. D. pulicaria dominates in some years, D. galeata mendotae dominates in others, and some years involve a seasonal replacement. This negative correlation results at least in part from temporal variation in the intensity of sizeselective fish predation (25) and interspecific competition (24). For Daphnia, like many species of zooplankton, the persistent lifehistory stage is a diapausing egg. In Oneida Lake, both species produce diapausing eggs from late May to early June, and in some years again in the late fall (24). These eggs have accumulated to sediment densities reaching tens of thousands per m^{2} and have the potential to remain viable for more than 125 years (24). Hatching of these eggs back to the water column occurs each spring after iceout but does not continue throughout the summer (24). That these two species are competitors and that both species posses a longlived life history stage satisfies two of the three necessary conditions of the storage effect model. Evidence that the third criterion, temporal fluctuations in recruitment, also exists for the two Daphnia species in Oneida Lake is presented below.
When the longlived life history stage is an egg bank, the number of eggs (X_{i}) in each species’ (i) bank changes through time (t) as: 1 where H_{i} is the hatching fraction, s_{i} is the probability of survival in the egg bank, and R_{i} is the annual per capita recruitment to the egg bank from hatched eggs (which may vary temporally) (19). Eq. 1 reflects only the dynamics of the egg bank and does not explicitly include the dynamics of the Daphnia in the water column. Consequently, “per capita recruitment” is based on the number of dormant eggs that hatch and produce new dormant eggs, irrespective of the number of Daphnia in the plankton. However, the number of Daphnia in the water in any given year has the potential to affect recruitment through competition.
The amount of generational overlap [(1 − H_{i})s_{i}] dictates how long an egg remains in the egg bank. Although hatching from the egg bank can result in the production of new diapausing eggs, it also represents a loss term from the longlived stage. The survivorship term for the dormant eggs reflects three additional loss processes: predation, senescence, and deep burial. Because dormant eggs need to be near the sediment surface to hatch, if an egg is buried below the top few centimeters, it is effectively lost unless sediment mixing returns it to the surface.
Coexistence of competitors occurs when each species is able to increase from low numbers (22, 26), that is, when each population’s theoretical mean boundary growth rate (population growth rate at low density) is positive: 2 In Eq. 2, the amount of generational overlap in the egg bank [(1 − H_{i})s_{i}] has been replaced by the general term (1 − δ_{i}), where δ represents the sum of all four loss processes from the egg bank. In addition, recruitment remains dependent on the hatching fraction but is expressed simply as R. Alternatively, Warner and Chesson (14) express Eq. 2 as: 3 where ρ_{i}(t) = R_{i}(t)/δ_{i}. If ρ_{i}(t) > 1, then recruitment exceeds loss and species i increases and persists, but if ρ_{i}(t) < 1, species i will decrease.
Both the storage of reproductive potential and recruitment fluctuations are likely to occur in many systems; hence, Warner and Chesson (14) provide an explicit analysis for quantifying the contribution of the storage effect mechanism to persistence. The contribution of the nonstorage component to the longterm population growth rate addresses the question of whether or not the population can increase without the extremely good recruitment years (and their subsequent storage) that result from environmental fluctuations. In other words, if a population can increase without recruitment fluctuation, then something other than the storage effect must permit longterm persistence. To quantify this nonstorage component, ρ̃_{i} is defined as the geometric mean of ρ_{i}, which removes the yeartoyear variance in R_{i}. The average over all years (j) gives the population’s mean instantaneous growth rate without environmental variation (Γ): 4 The storage effect’s contribution to the growth rate includes the annual variance in R_{i} and is the overall mean instantaneous growth rate with the nonstorage component (above) removed: 5 According to Warner and Chesson (14), if competition has been demonstrated and Eq. 4 is negative while Eq. 5 is positive for one or both species, then the storage effect is important in maintaining coexisting species in the system. If Eq. 4 is positive, then that species most likely does not depend on storage for its longterm persistence, because even without environmental fluctuations the population still maintains a positive growth rate. In short, the storage effect maintains diversity because the storage protects a rare species from extinction and the recruitment fluctuations lead to periodic large population increases.
Maintenance of diversity via the storage effect implicitly assumes that the species are competitors and that competition affects recruitment. That competition is occurring must be demonstrated experimentally, independent of a mathematical analysis of the longlived stage. Consequently, as a conservative recommendation, Warner and Chesson (14) advised that recruitment fluctuations should be negatively correlated between species. Note however, that this negative correlation is not, in fact, required for coexistence by the storage effect. Rather, interspecific competition must limit recruitment in a way that is affected strongly by the environment.
Parameterizing this model for the Oneida Lake daphnid assemblage requires estimates of both the recruitment and loss rate for each species’ egg bank. Annual estimates of recruitment to the egg bank were calculated from 21 years of plankton samples collected by the Cornell Biological Field Station. For each year, the number of ephippia produced per m^{2} in the water column were estimated from a weekly plankton sample and converted to daily production rates by dividing the number of ephippia per female by the length of the molt cycle (27–29). By integrating daily production rates over each sampling interval, weekly production rates were calculated. I obtained annual production estimates by summing over all sampling dates. Annual estimates of ephippial production were converted to annual diapausing egg production by multiplying by two eggs per case and correcting for the fact that approximately 10% of females shed empty cases (personal observations). In the years 1975 to 1991, one plankton sample from the archived collection of the Cornell Biological Field Station was examined for each week, whereas my own duplicate weekly samples were examined in 1992 to 1995. These numbers were converted to per capita recruitment rates based on the average number of eggs that hatch [10 eggs per m^{2} for D. pulicaria, 5 eggs per m^{2} for D. galeata mendotae (24)].
The second parameter in the model is the loss rate from the egg bank. In the theoretical exploration of this model, Warner and Chesson (14) assumed that loss rates were small and constant across years, because as the per capita death rate approaches 1, the longlived stage is eliminated and there is no storage. The actual loss rate from the Oneida Lake egg bank results from four processes, hatching, deep burial, predation, and senescence, of which only hatching has been quantified (see above). The rate at which eggs are lost due to deep burial depends on the depth to which the hatching cue penetrates. The limited data that exist suggest that eggs may only receive the appropriate light and oxygen cues within the top 0.15–0.33 cm of sediment (N. G. Hairston, Jr., personal communication). However, due to the potential for vertical mixing within the sediment, precisely quantifying the loss due to deep burial is not straightforward. Moreover, although I have evidence of egg mortality due to both predation and senescence (24), quantitative estimates are unavailable.
Lastly, because there are estimates for both the annual input to the egg bank from the water column and the number of eggs in the sediment (24), Eq. 1 must be balanced through the death rate. Because the estimate of eggs produced often greatly exceeded the number found in the top layer of the sediment, the recruitment values may be consistently overestimated (24). If this is the case, this bias can be corrected in Eq. 1 by either reducing all recruitment rates by a constant fraction or by increasing the death rate. Because the storage effect results from the interaction between recruitment variation and generational overlap, growth rates for each species were calculated for 10 loss rates ranging from 0.01 to 0.90. For each estimate of the death rate, growth rates were calculated twice, once using the calculated recruitment values (high value), and once using values that had been reduced by 50% (low value).
The above calculations assume that within the active egg bank the sediment is homogeneous. No allowances are made for the potential of vertical mixing within the sediment [although more recently, Chesson (17) has provided a theoretical analysis of the effects of spatial heterogeneity]. In Oneida lake, the sediment egg bank extends to a depth of 25 cm and exceeds 2.5 × 10^{4} eggs per m^{2} (24). Although some of these deeply buried eggs can be induced to hatch in the laboratory, they cannot contribute to the growth rate of the population unless they are mixed up to the surface. To incorporate the effects of sediment mixing, I constructed a matrix model for a threelayered system and computed the boundary growth rates of each population. The first layer included eggs from 0 to 0.15 cm in the sediment, the second encompassed 0.15 cm to 2 cm, and the third compartment included all eggs below 2 cm. The depth of the first layer was based on the depth to which the hatching cue penetrates (N. G. Hairston, Jr., personal communication) and the second layer extends to the typical depth of bioturbation in freshwater systems (30, 31). The generalized transition probabilities for the matrix are: The terms in the matrix are all annual per capita rates and represent the recruitment rate to the egg bank (R), the rate of upward mixing between the initial (i) and final (f) layer (U_{i,f}), the rate of passive burial downward to a lower layer (B_{i,f}), the rate of active downward mixing (M_{i,f}), and the proportion of eggs in layers two and three that do not move (P). Six separate matrices were analyzed with survivorship (S) ranging between 0.99 and 0.10. Eggs in the top layer either hatch or move downward, and the hatching and mortality are subsumed into the recruitment term.
On the basis of my field data for recruitment, sedimentation rates, and vertical distributions of eggs and literature numbers for mixing (30), the initial values were set at: The rate of passive downward movement (burial) reflects the sedimentation rate in Oneida Lake (24). Because the rate of downward mixing has been shown to be less than upward transport (30, 32), M is always less than U. On the basis of 1.5 years of sediment core data that showed no substantial upward movement into the top 2 cm (24), these rates of active transport by invertebrates were assumed to be quite small.
To calculate the growth rate without environmental variation (analogous to Eq. 4), I used matlab 4.0 to find the logarithm of the dominant eigenvalue of matrix A with R set as the geometric mean of the 21 annual recruitment values (years 1975–1995). In the variable recruitment case, the number of eggs in each compartment was set to low densities at a ratio that reflected the egg distribution at these three depths in the sediment of Oneida Lake (1:19:230 for D. galeata mendotae and 1:19:307 for D. pulicaria). For each time step, a value of R was selected at random from the 21 years of field recruitment data. Simulations of X_{i}(t + 1) = A × X_{i}(t) were run for as many generations as necessary for the longterm growth rate to stabilize (up to 1000). Stabilization of the growth rate was confirmed by calculating two intermediate growth rates within each simulation (e.g., for a 1000generation simulation, growth rates were calculated between the 500th and 750th generations and between the 750th and 1000th generations). For each species, I ran 10 simulations and found the final growth rate by difference from the population over the longest interval (e.g., for a 1000generation run, between population sizes at 500 generations and 1000 generations). These 10 growth rates were averaged to estimate the boundary growth rate. The value obtained from the constant recruitment case (analogous to Eq. 4) was subtracted from this number to give a result analogous to that of Eq. 5.
RESULTS AND DISCUSSION
Annual recruitment to the egg bank varied for both D. pulicaria and D. galeata mendotae (Fig. 1). Although each species experienced years in which it had higher recruitment than its competitor, D. pulicaria was favored more often. Over the 21 years investigated, D. pulicaria produced more diapausing eggs than D. galeata mendotae in 19 years, whereas the reverse was true only in 2 years. Moreover, D. galeata mendotae experienced 8 years of complete recruitment failure, whereas D. pulicaria produced at least a few hundred new diapausing eggs per m^{2} each year. This temporal fluctuation in recruitment satisfies the third necessary condition of the storage effect hypothesis.
This variation in recruitment to the egg bank is explained in part by the density of Daphnia in the water column. For example, for D. pulicaria, which produced the majority of its dormant eggs in the spring, more than 50% of the variation in the springtime production estimates is explained by the average spring D. pulicaria density (regression, logarithm transformed, r^{2} = 0.52, P = 0.0003). The more D. pulicaria that are present in the plankton during dormant egg production, the higher the recruitment to the egg bank. I have shown elsewhere (24) that the production of immediately hatching eggs by D. pulicaria can be reduced by the presence of D. galeata mendotae and vice versa. When fewer immediately hatching eggs are produced by a population, there are fewer individuals in the water column to produce the dormant eggs.
The Warner and Chesson calculations for D. galeata mendotae, using the “low” estimates of recruitment, suggest that without a fluctuating environment (Eq. 4), this species would not persist in Oneida Lake (Fig. 2A). Even with very low loss rates (i.e., high generational overlap), the portion of the growth rate based solely on the mean recruitment is negative. The qualitative results of the analysis did not differ when the “high” estimates of recruitment were used (data not plotted). Without environmental fluctuation, the species cannot persist. The primary reason for this result is that this species often experienced years of complete recruitment failure, eight times for the 21year period studied herein. Moreover, D. galeata mendotae suffered at least four consecutive years of zero recruitment to the egg bank (Fig. 1B). Without a multipleyear dormant stage and occasional years in which it is strongly favored, this species would have been eliminated from the community. In contrast, the D. pulicaria population does not rely on the storage effect for persistence (Fig. 2C). That is, the geometric mean recruitment rate is always greater than the per capita loss rate. This qualitative result is found by using a range (1–10^{3}) of estimates for the number of eggs hatching (and, therefore, values of per capita recruitment).
As pointed out by Warner and Chesson (14), it is essential when using Eqs. 4 and 5 that ρ be measured only at low population densities, when the effects of competition are weak. Because competition occurs among the active individuals in the water column, I found the median springtime daphnid densities and again calculated all values for Eqs. 4 and 5 by using only those years that fell below the median population size (1977–1981, 1983, 1985, 1986, 1993, and 1994). For both D. galeata mendotae and D. pulicaria, the qualitative results do not change.
The estimates of each population’s boundary growth rate obtained from the layered eggbank model support the qualitative results obtained with the unlayered model (Fig. 2 B and D). For D. galeata mendotae, the boundary growth rate is negative without variation. In contrast, D. pulicaria persists through consistently high recruitment.
CONCLUSIONS
The results of these analyses indicate that even though both daphnid species exhibit the necessary prerequisites for persistence via the storage effect, this mechanism may only be necessary in maintaining the weaker competitor in Oneida Lake. However, even though the D. pulicaria population does not rely on the storage effect for persistence, dormant eggs are nevertheless an important stage in its shortterm dynamics. This population is often completely eliminated from the plankton by midsummer and not seen again until the following spring (25). Therefore, without the reliable colonization source provided by dormant eggs, D. pulicaria would not have persisted in Oneida Lake. Moreover, longlived diapausing eggs are likely to be present in many other systems where competing zooplankton are found to coexist (33, 34). However, the role of this dormant stage and the overlapping generations that it creates are almost never considered when the community dynamics are examined. This oversight limits our understanding of community dynamics because diapause may prove to be a common mechanism promoting species coexistence in aquatic systems.
Acknowledgments
I thank S. Ellner for informative conversations regarding both the storage effect and the layered egg bank models and E. L. Mills for allowing me to use his archived plankton samples. Earlier drafts of the manuscript were improved by the comments of A. J. Bohonak, P. L. Chesson, C. R. FreemanGallant, M. A. Geber, N. G. Hairston, Jr., C. D. Harvell, E. L. Mills, and an anonymous reviewer. This work was supported by National Science Foundation Doctoral Dissertation Improvement Grant DEB9410966 and by a National Science Foundation predoctoral fellowship.
Footnotes

↵† Present address: Illinois Natural History Survey, 607 East Peabody Drive, Champaign, IL 61820. email: caceres{at}mail.inhs.uiuc.edu.

This paper was submitted directly (Track II) to the Proceedings Office.
 Received January 13, 1997.
 Copyright © 1997, The National Academy of Sciences of the USA
References
 ↵
 Connell J H
 ↵
 Diamond J,
 Case T J
 Chesson P L
 ↵
 Rothhaupt K O
 ↵
 ↵
 Schoener T W
 ↵
 ↵
 ↵
 ↵
 ↵
 Pamilo P
 ↵

 Freeman H I,
 Strobeck C
 Chesson P L
 ↵

 Chesson P L,
 Huntly N
 ↵
 ↵
 Chesson P L
 ↵
 ↵
 Rhodes O E Jr,
 Chesser R K,
 Smith M H
 Hairston N G Jr,
 Ellner S,
 Kearns C M
 ↵
 Leigh E G Jr
 Leigh E G Jr
 ↵
 Buttler A J,
 Chesson P L
 ↵
 ↵
 Hairston N G Jr
 ↵
Cáceres, C. E. (1997) Ph.D. dissertation (Cornell Univ., Ithaca, NY).
 ↵
 Carpenter S R
 Mills E L,
 Forney J L
 ↵
 Turelli M
 ↵

 Paloheimo J E
 ↵
Braner, M. (1988) Ph.D. dissertation (Cornell University, Ithaca, NY).
 ↵
 ↵
 McCall P L,
 Tevesz J S
 McCall P L,
 Tevesz J S
 ↵
 Davis R B
 ↵
 ↵
Citation Manager Formats
More Articles of This Classification
Biological Sciences
Related Content
 No related articles found.
Cited by...
 Spatial storage effect promotes biodiversity during adaptive radiation
 Bounded population sizes, fluctuating selection and the tempo and mode of coexistence
 Understanding past, contemporary, and future dynamics of plants, populations, and communities using Sonoran Desert winter annuals
 Functional tradeoffs determine species coexistence via the storage effect