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
Global biodiversity and the ancient carbon cycle

Communicated by John M. Hayes, Woods Hole Oceanographic Institution, Woods Hole, MA (received for review August 21, 2000)
Related Article
Abstract
Paleontological data for the diversity of marine animals and land plants are shown to correlate significantly with a concurrent measure of stable carbon isotope fractionation for approximately the last 400 million years. The correlations can be deduced from the assumption that increasing plant diversity led to increasing chemical weathering of rocks and therefore an increasing flux of carbon from the atmosphere to rocks, and nutrients from the continents to the oceans. The CO_{2} concentration dependence of photosynthetic carbon isotope fractionation then indicates that the diversification of land plants led to decreasing CO_{2} levels, while the diversification of marine animals derived from increasing nutrient availability. Under the explicit assumption that global biodiversity grows with global biomass, the conservation of carbon shows that the longterm fluctuations of CO_{2} levels were dominated by complementary changes in the biological and fluid reservoirs of carbon, while the much larger geological reservoir remained relatively constant in size. As a consequence, the paleontological record of biodiversity provides an indirect estimate of the fluctuations of ancient CO_{2} levels.
Partly as a result of the episodic nature of sedimentary processes, studies of the Earth's biologic, geologic, and climatic history have typically focused on isolated extreme events (1) such as mass extinctions (2–4) or anomalous excursions of stable chemical isotope records (5). Whereas the importance of unusual events may be nearly selfevident (6), a different, potentially powerful approach is to seek not what makes some events different from others, but that which is common to them all (7–10). In this way, one should expect to learn how fundamental invariants—for example, the conserved mass of certain chemical species—manifest themselves over geologic time.
Here we focus on biological aspects of the carbon cycle during Phanerozoic time [0–543 million years ago (Ma)]. Two types of records are studied: paleontological data for the diversity of marine animal (11) and terrestrial plant (12) life, and the sedimentary record of carbon isotope fractionation between carbon buried in organic form and carbon buried as carbonates (13). We find an unexpected significant correlation between the paleontological and geochemical records for the last 400 million years and provide evidence that it reflects the fluctuations of CO_{2} levels in the atmosphere and oceans. Our analysis shows not only that the diversification of life appears to have depressed CO_{2} levels (14), but also how deep this depression appears to have been.
The paper proceeds as follows. We first show the correlations of interest and subject them to several tests of statistical significance. We then discuss why CO_{2} levels, and specific mechanisms that control them, appear to provide the link between the various records. Simple mathematical models are then formulated to clarify and quantify our arguments. Among the results is an estimate of trends in the longterm evolution of CO_{2} levels, which we proceed to compare with previous efforts (15–19). Finally, we conclude with brief remarks on the relationship of our quantitative conclusions to their mechanistic interpretation.
PaleontologicalGeochemical Correlations
The first of two correlations is shown in Fig. 1. The data for diversity were compiled by Sepkoski† and edited by Bambach‡; they are based on the first and last geologic occurrence of 35,967 marine animal genera, including invertebrates, vertebrates, and animallike protists (foramaniferans and radiolarians). For each of 108 geologic stages and substages that range over Phanerozoic time, the number of genera that originate or become extinct during that stage is tabulated. If a genus passes from one (sub)stage to the next (a condition satisfied by approximately 63% of the genera), it is labeled “continuing.” The time series n(t) in Fig. 1 represents the number of continuing genera at each stage or substage boundary.
The carbon isotope record in Fig. 1 was compiled by Hayes et al. (13). The time series shown is obtained from measurements of the isotopic composition of organic and carbonate carbon found in sedimentary rocks formed at geologic time t. From globally averaged measurements of isotopic abundance ratios R_{x} = (^{13}C/^{12}C)_{x} for carbon in sample x, the isotopic fractionation between sample x and a standard sample, δ_{x} = 1,000[(R_{x} − R_{STD})/R_{STD}], is obtained for carbonate (δ_{a}) and organic (δ_{o}) carbon. The isotopic fractionation ɛ_{toc} between total organic carbon and sedimentary carbonates is then given approximately by ɛ_{toc} = δ_{a} − δ_{o} (13). Fig. 1 compares ɛ_{toc}(t) to n(t); the values of ɛ_{toc}, which were obtained from global averages of thousands of independent measurements (13), are specified at 50 discrete points in time.
The two records displayed in Fig. 1 are strikingly similar after 455 Ma. Not only do both records change sharply (increasing n, decreasing ɛ_{toc}) at the end of the Cretaceous (65 Ma), but visually one can see that many smaller transient fluctuations are exhibited jointly. To quantify this correlation, the more finely sampled paleontological data are resampled by linear interpolation so that they are specified at the same times t_{j} as the geochemical data. The Spearman rank correlation coefficient R_{s} (25, 26) is then computed from the resulting N pairs of n and ɛ_{toc}. We find R_{s} = −0.74 after 455 Ma. Similar results, here and elsewhere in this paper, are obtained if “standing” diversity, the number of genera extant during a given interval, is substituted for continuing diversity.
The second correlation, shown in Fig. 2, is between the diversity of land plant families (12) and the same geochemical data. The plant data§ represent standing diversity at the family level in 67 geologic stages since the rise of land plants in the Silurian (425 Ma). From the late Devonian (370 Ma) onwards, the correlation between the two records is essentially the same as that found in Fig. 1: R_{s} = −0.77.
We use three different methods to assess the statistical significance P of these correlations. The results are summarized in Table 1. For completeness, we report results for the entire available interval in addition to the slightly shorter periods cited above. We argue later that only the latter periods are relevant to our conclusions.
In the first method, a Monte Carlo technique, we construct random realizations of each time series that satisfy the autocorrelations of the original data. Denoting either the diversity or ɛ_{toc} data as x_{j} = x(t_{j}), the method proceeds by Fourier transforming x [after multiplying by a Hanning window (26) and removing the mean] to obtain N complex Fourier coefficients x̂_{k}. Random time series x′_{j} are then generated by computing the inverse Fourier transform such that x′_{j} = ∑_{k}x̂_{k}e^{i(φk}^{+2πjk/N)}, where φ_{k} is uniformly distributed between −π and π while satisfying the symmetries required for x′_{j} to be real. P values are then estimated by computing the fraction of 10^{6} simulated timeseries pairs that yield rank correlation coefficients R_{s} less than those obtained for the actual data.¶
The second method applies a firstorder predictionerror filter (27) to each time series. In principle, this converts the time series to white noise, thereby eliminating the additional variance of the estimated correlation coefficient caused by significant autocorrelation of x_{j} at nonzero lags (28). Denoting the whitenoise sequence by η_{j} = x_{j} + ax_{j−1}, we choose a by minimizing the sum of the squared forward and backward prediction error, ∑(x̃_{j} + ax̃_{j−1})^{2} + (x̃_{j−1} + ax̃_{j})^{2}, where x̃_{j} = x_{j} − N^{−1}∑x_{j}. We obtain a = −0.98 for both diversity series and a = −0.73 for ɛ_{toc}. The rank correlation coefficient between η_{j}(n) and η_{j}(ɛ_{toc}) is then computed, along with the associated onesided statistical significance P obtained from Student's tdistribution with N − 3 degrees of freedom. The Insets in Figs. 1 and 2 show plots of η_{j}(n) vs. η_{j}(ɛ_{toc}).
The third method, included as a reference, constructs η_{j} by computing serial differences, i.e., by setting a = −1 above. This is equivalent to assuming that all three time series are realizations of a random walk.
Table 1 shows that all three methods find that the negative correlations quoted above are statistically significant, usually with P < 0.01. The Monte Carlo method typically yields the lowest P values, whereas serial differencing yields the highest. Because the former method requires the fewest assumptions and the latter the most, it seems reasonable to associate the greatest credence with the Monte Carlo results.
Decomposition of the Geochemical Signal
To understand why the paleontological and geochemical data correlate negatively, consider the following decomposition of ɛ_{toc} (13): 1 Here ɛ_{p} represents the isotopic effects of primary production (i.e., photosynthesis), Δ_{carb} is the isotopic depletion of dissolved CO_{2} in surface waters relative to sedimentary carbonate, and Δ_{2} represents isotopic shifts caused by the reworking of organic carbon before burial. The quantity Δ_{carb} depends mainly on surface water temperature and ranges from about 7 per mil (‰) at 30°C to 10‰ at 3°C (13); this 3‰ range is much less than the 12‰ range of ɛ_{toc}. The secondary processes symbolized by Δ_{2} are defined such that Δ_{2} is positive when the process enriches the ^{13}C content of organic matter. One of the principal components of Δ_{2} derives from the structure of food webs: the ^{13}C content of many animals is enriched by approximately 1‰ relative to its food (29–31). Thus as diversity grows, one might naively expect the size of food webs and therefore the length of food chains to grow, and thus also Δ_{2}, thereby explaining the negative correlation of ɛ_{toc} with n. However, contemporary food chains nearly always contain about four trophic levels, and this length depends so weakly on the size of the food webs within which the chains are embedded (32) that it appears unlikely that food chains are responsible for the correlation. Moreover, the net effect of all secondary biological processes is such that even the sign of Δ_{2} is debated (31), and recent estimates put Δ_{2} at the relatively small and constant value of 1.5‰ throughout the Phanerozoic (13). In summary, the combined fluctuations of Δ_{carb} and Δ_{2} should amount to no more than about 3‰, which is less than onefourth of the variation of ɛ_{toc} during the Phanerozoic.
The photosynthetic processes represented by ɛ_{p} (13, 31, 33, 34) thus encode the only plausible explanations for the bulk of the ɛ_{toc}'s variations. Recent research indicates that plankton photosynthesis produces an isotopic shift in ^{13}C content that may be approximated by the empirical relation (34) 2 where ɛ_{f} ≃ 25–30‰ (35) is the isotopic shift associated with carbon fixation, κ is proportional to the permeability of the algal cell wall, υ_{s} is the volumetosurfacearea ratio of the cell, μ is the specific growth rate, and [CO_{2}] is the concentration of dissolved carbon dioxide. The isotopic enrichment ɛ_{f} derives primarily from the kinetic discrimination against ^{13}C by the enzyme Rubisco in the CalvinBenson or photosynthetic carbon reduction cycle. It may be expected to have remained approximately constant in marine environments since the first appearance of eukaryotes at least two billion years ago (35). However, the physiological parameters κ, υ_{s}, and μ can vary significantly across phytoplankton species (33, 34). Along with [CO_{2}], they control the relative rate at which CO_{2} diffuses into the algal cell compared to the rate at which CO_{2} is consumed by photosynthetic reactions within the cell. Whereas the longterm evolution of [CO_{2}] should be independent of these parameters, they may have evolved in response to changes in [CO_{2}]. For example, decreases in [CO_{2}] could have led to increases in cell permeabilities κ such that the diffusive flux of CO_{2} into the algal cell was relatively unchanged. Provided that such evolutionary innovations are not fully compensating—i.e., so long as κ[CO_{2}] increases with [CO_{2}]—then the correlations between diversity and ɛ_{toc} should reflect correlations between diversity and [CO_{2}]. More generally, this conclusion remains valid for any dependence of the physiological parameters on [CO_{2}], provided that μυ_{s}/(κ[CO_{2}]) decreases with increasing [CO_{2}] and no direct relationship between the physiological parameters and diversity exists.
Mechanisms and Dynamics
Unlike the shortterm carbon cycle, which is dominated by exchanges of carbon between the biosphere, atmosphere, oceans, and soils on time scales ranging from about 10^{0}–10^{4} years, the longterm carbon cycle is dominated by exchanges between rocks and the atmosphere and oceans on time scales of roughly 10^{5}–10^{9} years (14–16, 36, 37). Three slow processes dominate the longterm evolution of CO_{2} levels: chemical weathering of silicate minerals, degassing caused by metamorphic and volcanic processes, and the burial and weathering of organic carbon. Because the latter acts as an independent subcycle (14), the former—weathering and degassing—are considered the major controls on atmospheric and oceanic CO_{2} levels. They are thought to be in rough balance via reactions such as (14, 36, 38) 3 A similar reaction may be written by substituting Mg for Ca. Left to right, such reactions schematically represent the uptake of CO_{2} from the atmosphere, its transformation to dissolved HCO during weathering of silicate rocks, and its eventual precipitation and burial in the oceans as carbonate minerals. Right to left, the reaction represents metamorphism and magmatism and the subsequent transfer of CO_{2} back to the atmosphere and oceans by volcanism and related processes.
The mechanistic basis of our observed correlations is simple: increased plant diversity should result in increased chemical weathering of silicate minerals. We assume, following others (14–18, 39), that the increased chemical and physical alterations of soils—e.g., the mechanical breakdown of rocks and the production of organic acids—induced by the advent of vascular land plants should have led to substantial increases in weathering rates. This notion is augmented with two assumptions: (i) increased plant diversity increases the area occupied by plants (by, for example, migration of plants to higher altitudes), therefore enhancing weathering rates and lowering CO_{2} levels via Eq. 3; and (ii) increased weathering increases the seaward flux of limiting nutrients such as phosphorus, which in turn leads to increased marine animal diversity. The former assumption is sufficient to explain the correlation between land plant diversity and ɛ_{toc}, whereas the latter is also required to explain the correlation with the marine ecosystem. Both of these assumptions may be subsumed in part by the simpler supposition that global diversity increases with global productivity or biomass, a point we return to below.
To complete our argument, we must show not only that increased plant diversity can lead to lower CO_{2} levels, but also that the lower CO_{2} levels are dynamically stable. The rate of uptake of CO_{2} by the weathering reactions should increase with the partial pressure of atmospheric CO_{2}, which at long time scales should be an increasing function of the total mass m_{f} of carbon in the Earth's fluid envelope (i.e., atmospheric and oceanic carbon) (14). We represent the bare uptake rate, in the absence of plants, by the positive monotonically increasing function k(m_{f}) while representing the return flux of CO_{2}, from volcanic and other sources, by υ [which should not depend on m_{f} (14)]. All other processes aside, the rate of change of the mass of atmospheric and oceanic carbon is then 4 where g(n) ≥ 1 is a monotonically increasing function that represents the ratio of the CO_{2} uptake rate when plant diversity equals n to the bare uptake rate k.
The steady state m^{*}_{f} of this system occurs when ṁ_{f} = 0; thus 5 where k^{−1} is the inverse of k such that k^{−1}[k(x)] = x. Because k^{−1} must also be a monotonically increasing function, Eq. 5 shows that increasing plant diversity n causes the steadystate mass of atmospheric and oceanic carbon, m^{*}_{f}, to decrease. The steady state is stable if dṁ_{f}/dm_{f} < 0 when m_{f} = m^{*}_{f}. From Eq. 4, one finds that m^{*}_{f} is indeed stable, because 6 for any m_{f}. This result follows immediately by noting that dk/dm_{f} > 0 and g(n) ≥ 1.
A complete representation of the longterm carbon cycle is of course vastly more complex (14–16, 36). The model expressed by Eq. 4 does, however, serve to formalize our main point: if increasing plant diversity increases weathering rates, then CO_{2} levels fall. We proceed to quantify by how much.
Relative Fluctuations of Carbon Reservoirs
We consider the mass of three reservoirs of carbon: the fluid reservoir, already denoted by m_{f}, the biological reservoir m_{b} (i.e., biomass plus biological debris), and the geological reservoir m_{g} (i.e., rocks). Presentday estimates give m_{b} ≃ 10^{−1}m_{f} ≃ 10^{−5}m_{g} (35, 37).
In globally averaged steadystate conditions, [CO_{2}] increases with m_{f} (40). The inverse relationship is likewise increasing, i.e., 7 where h_{f} is a monotonically increasing function.
A monotonic relationship probably also exists between diversity and m_{b}. Waide et al. (41) have recently compiled a comprehensive review of over 200 published relationships between presentday diversity and primary productivity, where primary productivity is broadly defined to include various surrogates, such as biomass, that vary positively with it. One of the few unambiguous patterns to emerge from their survey is the dependence of the diversity–productivity relationship on scale. Whereas at local scales (<20 km) simple relationships (positive, negative, or unimodal) are usually not found [or when they are, they are often unimodal (42)], at continental to global scales (>4,000 km), a positive correlation dominates the literature, particularly for plants. Additional evidence of such scale dependence is given in ref. 43, and examples of unambiguously positive largescale relationships may be found in ref. 44. The positive largescale relation may derive from the tendency of ecosystems to diversify when more energy is available for consumption (45). Similar relationships appear to hold for the past: Bambach (46) cites a wealth of evidence indicating that the biomass of marine animals has trended upward during the Phanerozoic. Thus it appears reasonable to assume that the mass m_{b} of the biological reservoir of carbon increases with diversity n, i.e., 8 where h_{b} is monotonically increasing. The assumption that increased diversity is associated with increased weathering rates may be viewed as a natural consequence of this relationship.
We may now relate the correlations between n and ɛ_{toc} to m_{b} and m_{f}. First, rewrite Eq. 2 as 9 where γ = μυ_{s}/κ represents the physiological parameters appropriately averaged over all phytoplankton species, and 10 Now assume that ɛ_{0} is constant and γ/[CO_{2}] decreases with increasing [CO_{2}]. In the former case, the only time dependence should come from Δ_{carb}, and as we have already remarked, its meansquare fluctuations must be much smaller than those of γ/[CO_{2}]. Thus this assumption appears reasonable. The latter assumption is satisfied automatically if γ is constant. The less stringent condition used here derives from the aforementioned evolutionary considerations, i.e., γ can change, but never so much that it fully compensates (or overcompensates) changes in [CO_{2}].
Let R_{l}(x, y) and R_{s}(x, y) be the functions that compute the linear [productmoment (25)] and rank correlation coefficients, respectively, between discretely sampled time series x(t) and y(t). Because R_{s} is a correlation of ranks, it is invariant with respect to any monotonically increasing transformation of its arguments (25). From Table 1 and Eqs. 7, 8, and 9, we have 11 12 13 14 These results derive from the invariance property of R_{s}, the condition that [CO_{2}] > 0, and the assumptions that ɛ_{0} is constant, γ/[CO_{2}] decreases with increasing [CO_{2}], and γ is uncorrelated with n. The final equality shows that the mass fluctuations of the biological and fluid reservoirs of carbon are negatively correlated, with the same magnitude as found for diversity and ɛ_{toc}.
This conclusion may be sharpened. In the following, we assume that any measurable coupling between m_{b} and m_{f} requires m_{b} to be sufficiently large, as a consequence of sufficient diversification on land. We therefore explicitly consider only the period after 370 Ma, when both marine animals and land plants exhibit their strongest correlations to ɛ_{toc}.
Because the Earth's total carbon is conserved, 15 The fluctuations m̃_{i} = m_{i} − 〈m_{i}〉, where 〈⋅〉 signifies the mean over the entire interval, satisfy 16 Squaring both sides, rearranging terms, and averaging over time, one obtains 17 We seek the sign of 〈m̃_{b}m̃_{f}〉. Mathematically, it must have the same sign as the linear correlation coefficient R_{l}(m_{b}, m_{f}). R_{l} and its variance may be estimated from R_{s} under the assumption that m̃_{b} and m̃_{f} are Gaussian. We use the relations R_{l} = 2 sin(πR_{s}/6) and var R_{l} = (π^{2}/9)(1 − R/4) var R_{s} (25). Substituting R_{s} = −0.40 from Table 1 (i.e., the estimate obtained by prediction–error filtering) and an estimate of var R_{s} obtained by the bootstrap method (47) (averaged for both marine animals and land plants), we find R_{l} = −0.42 ± 0.13 (one standard deviation). Consequently, 〈m̃_{b}m̃_{f}〉 < 0 with high probability, and therefore, on substitution into Eq. 17, 18 This inequality, which we have obtained here for the longterm millionyear time scale, is an important counterintuitive result. It shows that the meansquare fluctuations of the geologic reservoir are less than those of the biological and fluid reservoirs combined—despite the fact that rocks contain nearly all the Earth's carbon! The longterm depletion of m_{f} is therefore not only due to the growth of m_{b}, but most of m_{f}'s fluctuations are complementary (i.e., of opposite sign) to m_{b}'s. Consequently, an estimate of the changes of either quantity provides an estimate of the changes of the other.
Estimate of Past CO_{2} Levels
We define the dimensionless CO_{2} concentration φ = [CO_{2}]ɛ_{0}/γ and substitute into Eq. 9 to obtain 19 In this formulation, the physiological parameters subsumed into γ no longer appear. We are obliged, however, to choose a value for ɛ_{0}. In the following, ɛ_{0} ranges from 36‰ to 38‰, which is consistent with the known variations of its constituents.
Rather than obtaining φ directly from Eq. 19, we estimate it indirectly from a combination of the geochemical and diversity data of Fig. 1. We fit the estimate of φ obtained from Eq. 19 to a linear model for φ(n): 20 In this model, φ_{0} is the dimensionless CO_{2} concentration before the diversification of animals and land plants, whereas the second term represents the firstorder effect of increasing diversity. For ɛ_{0} = 38‰, we find φ_{0} = 6.6 ± 0.4 and b = (8.0 ± 1.7) × 10^{−4} by linear regression.
Given estimates of φ_{0} and b, we then estimate φ(t) itself by setting n = n(t) in Eq. 20. The result, normalized by φ(0), is shown in Fig. 3 for the cases ɛ_{0} = 36‰ and ɛ_{0} = 38‰. Fig. 3 also shows a comparison with analogous predictions from the numerical model of Berner (15–18). Berner's model is a detailed integration of estimated carbon fluxes caused by chemical weathering, volcanic degassing, and carbon burial; the uncertainty of his estimates is qualitatively reflected by the gray area in Fig. 3.
Comparison of the results from our model with Berner's shows that the pattern of decreasing CO_{2} levels since ∼225 Ma is roughly equivalent in both. This trend is also qualitatively consistent with estimates based on the carbon isotopic composition of paleosols (48–51), plant stomatal ratios (52), carbon isotopic measurements from marine sediments (53), and an alternative reconstruction based on estimates of the Earth's inventory of buried carbon (19). However, the period before 225 Ma shows distinct differences from Berner's predictions. Although our estimate falls within Berner's region of uncertainty, it shows relatively no change compared to Berner's steep decrease at ∼300 Ma. Moreover, recent estimates from paleosols (54, 55) and plant stomatal ratios (52) are in accord with this dip. If the dip is indeed real, the mismatch with our model would presumably be caused by missing ingredients, such as changes in the burial rate of organic carbon accompanied by compensatory changes in algal cell permeabilities.
Conclusion
Surprising correlations exist between paleontological records of biodiversity and the carbon isotope fractionation evident in the sedimentary record for the last 370 million years. Proceeding from the assumptions that (i) chemical weathering rates grow with plant diversity; (ii) diversity grows with biomass at the global scale; and (iii) changes in phytoplankton physiology neither correlate with diversity nor fully compensate changes in CO_{2} levels, we find that the correlations express complementary fluctuations in the size of the organic and inorganic carbon reservoirs within the biosphere, atmosphere, and hydrosphere. Consequently, CO_{2} levels decreased as biodiversity increased. These conclusions imply that fluctuations of CO_{2} levels have been driven primarily by changes within the biosphere and only secondarily by purely geologic and geophysical processes. However, such a causal relationship is not required for our principal quantitative conclusions—the Inequality 18 and the downward trend in Fig. 3—to hold. Indeed, our analysis leaves open the possibility, for example, that tectonically induced reductions of CO_{2} levels led to increased diversity in both continental and marine ecosystems. A definitive statement of causation will require further work that includes not only the study of other geochemical signals (10) but also improved paleontological records.
Acknowledgments
I thank A. Knoll for introducing me to the study of Earth history, A. Knoll and R. Bambach for access to the late Jack Sepkoski's unpublished genuslevel database, and R. Bambach for his generous technical assistance with the paleontological data. I also thank J. Hayes, not only for access to his data before publication, but also for clarifying many of the finer points of carbonisotope geochemistry. Finally, I thank S. Bowring, E. Boyle, P. Dodds, J. Edmond, B. Enquist, D. Erwin, H. Hartman, C. Marshall, D. Raup, N. Schorghofer, and J. Weitz for helpful discussions. This work was supported in part by National Science Foundation Grant DEB0083983.
Footnotes

↵* Email: dan{at}segovia.mit.edu.

See commentary on page 4290.

↵† J. J. Sepkoski, Spring 1996 version of the genuslevel database that appeared in an earlier form in ref. 11.

↵‡ R. K. Bambach, personal communication. The original Sepkoski dataset was timed according to Harland et al. (20), with changes to the Cambrian following Bowring et al. (21). Bambach's updates include retiming of the Cenozoic (22), placing the Permian–Triassic boundary at 251 Ma (23), and the Permian–Carboniferous boundary at 301 Ma (24).

↵§ The compilation is available on the Internet at http://palaeo.gly.bris.ac.uk/frwhole/fr2.html.

↵¶ Because the theory constructed later requires the assumption that diversity and ɛ_{toc} are negatively correlated, we test specifically for the probability P that a negative correlation with magnitude at least as great as that observed could be obtained under the null hypothesis that diversity and ɛ_{toc} are uncorrelated.
Abbreviation
 Ma,
 million years ago
 Received August 21, 2000.
 Accepted January 30, 2001.
 Copyright © 2001, The National Academy of Sciences
References
 ↵
 Walliser O H
 ↵
 Erwin D H

 Kemp T S
 ↵
 Courtillot V
 ↵
 ↵
 Alvarez L W,
 Alvarez W,
 Asaro F,
 Michel W V
 ↵
 Berggren W A,
 Van Couvering J A
 Fischer A G

 Raup D M,
 Sepkoski J J

 Worsley T R,
 Nance D,
 Moody J B
 ↵
 ↵
 Walliser O H
 Sepkoski J J
 ↵
 Benton M J
 ↵
 ↵
 Walker J C G
 ↵
 Berner R A
 ↵
 Berner R A
 ↵
 Berner R A
 ↵
 Berner R A
 ↵
 Budyko M I,
 Ronov A B,
 Yanshin A L
 ↵
 Harland W B,
 Armstrong R,
 Cox A,
 Craig L,
 Smith A G,
 Smith D G
 ↵
 ↵
 Berggren W A,
 Kent D V,
 Aubry MP,
 Hardenbol J
 Berggren W A,
 Kent D V,
 Swisher C C III,
 Aubry MP
 ↵
 ↵
 Rasbury E T,
 Hanson G N,
 Meyers W J,
 Holt W E,
 Goldstein R H,
 Saller A H
 ↵
 Kendall M G,
 Gibbons J D
 ↵
 Press W H,
 Flannery B P,
 Teukolsky S A,
 Vetterling W T
 ↵
 Claerbout J F
 ↵
 Box G,
 Jenkins G M,
 Reinsel G C
 ↵
 ↵
 ↵
 Cohen J E,
 Briand F,
 Newman C M
 ↵
 ↵
 ↵
 Falkowski P G,
 Raven J A
 ↵
 Holland H D
 ↵
 Banfield J F,
 Nealson K H
 Des Marais D J
 ↵
 Urey H C
 ↵
 Algeo T J,
 Scheckler S E
 ↵
 Broecker W S,
 Peng TH
 ↵
 ↵
 Rosenzweig M L
 ↵
 Ricklefs R E,
 Schluter D
 Wright D H,
 Currie D J,
 Maurer B A
 ↵
 ↵
 ↵
 ↵
 Efron B,
 Tibshirani R J
 ↵
 Cerling T E

 Sinha A,
 Stott L D

 Andrews J E,
 Tandon S K,
 Dennis P F
 ↵
 ↵
 McElwain J C
 ↵
 ↵
 ↵
 Yapp C J,
 Poths H