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
High-precision timeline for Earth’s most severe extinction
Edited by Dennis Kent, Rutgers University and Lamont-Doherty Earth Observatory, Palisades, NY, and approved January 2, 2014 (received for review September 18, 2013)
This article has a correction. Please see:
Significance
Mass extinctions are major drivers of macroevolutionary change and mark fundamental transitions in the history of life, yet the feedbacks between environmental perturbation and biological response, which occur on submillennial timescales, are poorly understood. We present a high-precision age model for the end-Permian mass extinction, which was the most severe loss of marine and terrestrial biota in the last 542 My, that allows exploration of the sequence of events at millennial to decamillenial timescales 252 Mya. This record is critical for a better understanding of the punctuated nature and duration of the extinction, the reorganization of the carbon cycle, and a refined evaluation of potential trigger and kill mechanisms.
Abstract
The end-Permian mass extinction was the most severe loss of marine and terrestrial biota in the last 542 My. Understanding its cause and the controls on extinction/recovery dynamics depends on an accurate and precise age model. U-Pb zircon dates for five volcanic ash beds from the Global Stratotype Section and Point for the Permian-Triassic boundary at Meishan, China, define an age model for the extinction and allow exploration of the links between global environmental perturbation, carbon cycle disruption, mass extinction, and recovery at millennial timescales. The extinction occurred between 251.941 ± 0.037 and 251.880 ± 0.031 Mya, an interval of 60 ± 48 ka. Onset of a major reorganization of the carbon cycle immediately precedes the initiation of extinction and is punctuated by a sharp (3‰), short-lived negative spike in the isotopic composition of carbonate carbon. Carbon cycle volatility persists for ∼500 ka before a return to near preextinction values. Decamillenial to millennial level resolution of the mass extinction and its aftermath will permit a refined evaluation of the relative roles of rate-dependent processes contributing to the extinction, allowing insight into postextinction ecosystem expansion, and establish an accurate time point for evaluating the plausibility of trigger and kill mechanisms.
The ability to examine the rock record at millennial to decamillennial time scales in rocks that are hundreds of millions of years old permits critical evaluation of the patterns and rates of climate change, biological response to environmental perturbations, and evolution in deep time. This knowledge can give context to our understanding of the scale and rate of current biologic and climate change. In this article, we show that the largest known extinction in the history of animal life occurred in tens of thousands of years, just after a short-lived and major reorganization of the global carbon cycle.
Mass extinctions have long garnered attention, as they are characterized by fundamental restructuring of marine and terrestrial ecosystems and reflect complex feedbacks between environmental change, extinction, and recovery. However, it is only rarely that environmental perturbation leads to global extinction. Proposed drivers of mass extinctions range from asteroid impact to flood basalt volcanism, which are thought to trigger kill mechanisms ranging from global ocean anoxia to high atmospheric pCO2 to high ocean and atmosphere temperatures, for example (1). Although the geologic record is replete with occurrences of all of these, very few lead to mass extinctions. Studying the temporal details of mass extinctions is crucial for understanding how they are triggered and may allow isolation and identification of processes that are associated with a characteristic timescale. These processes in turn may be relevant to current biological and climate change and the timescales of feedbacks between environmental change and extinction.
The last two decades have seen a great deal of interest in the largest Phanerozoic extinction, the end-Permian biotic crisis, and an increased understanding of the patterns and timing of extinction and recovery, the synchronous and rapid perturbation of the global ocean-atmosphere system, and plausible trigger(s) and kill mechanisms (1⇓⇓⇓⇓–6). Attempts to reconcile the patterns and rates of extinction in marine and terrestrial environments has led to some agreement on the nature of severe environmental changes, including increased atmospheric pCO2 and acidification of the oceans, as well as widespread euxinic/anoxic conditions and a sharp spike in sea surface temperature (4, 7⇓⇓⇓⇓–12), and inferentially, kill mechanism(s).
Since 1998, four major U-Pb geochronological studies have attempted to constrain the timing and duration of the extinction (3, 13⇓–15). To better understand the relationship between environmental perturbation and biotic response, accurate and precise age models that integrate geochronology, paleontology, and geochemistry must be developed (8, 11, 16, 17). Recognition of astronomically forced sedimentary cycles (Milankovitch cycles) in late Permian and Triassic sedimentary rocks tuned with available geochronology have been used to refine existing age models of the biotic crisis (18⇓–20). Published estimates of the extinction interval based on radioisotopic dates range from ∼1.5 Mya to approximately <200 ± 100 ka, whereas astrochronological interpretations range from ∼700 ka to as little as ∼10 ka (3, 14, 15, 18, 19). Most recently, Wu et al. (20) used Milankovitch cyclicity and previously published geochronology to constrain the maximum extinction interval at Meishan to 83 ka. Published geochronology is not sufficiently precise to test this estimate.
The geology, biostratigraphy, and chemostratigrahy of the Global Stratotype Section and Point (GSSP) for the Permian-Triassic boundary at Meishan, China, have been previously described in detail (3, 11, 21, 22), and the section contains volcanic ash beds interlayered with fossil-bearing carbonate rocks (Fig. 1 and Fig. S1). In this article, we use the significant progress made in U-Pb geochronology since the acquisition of the data published in Shen et al. (3) (SI Text) and present more precise and accurate dates on zircon crystals isolated from five volcanic layers (beds 22, 25, 28, 33, and 34) that span the main extinction event, the major negative excursion and oscillation in δ13Ccarb, the Permian-Triassic boundary at Meishan as defined by the first appearance datum (FAD) of the conodont, Hindeodus parvus, and the earliest Triassic period (Figs. 1 and 2 and Figs. S1 and S2).
New Age Model
We rely on weighted mean 206Pb/238U dates as the best estimate of eruption/depositional ages for these volcanic rocks. For each sample, we determined a minimum of nine dates on single grains of zircon, allowing recognition of outliers due to either incorporation of older zircon or open system behavior, such as Pb loss (Fig. 2, Table 1, and Tables S1 and S2). The individual 206Pb/238U dates and the calculated weighed mean 206Pb/238U dates we present for beds 22, 25, and 28 are distinctly younger (up to 0.2%) and more precise than dates on the same ash beds published in Shen et al. (3) (Table S1 and Fig. S2). We also present a weighted mean 206Pb/238U date on bed 33, which did not yield a reliable date in the previous study (3). The differences in age and precision reported here reflect adoption of major improvements in the way our U-Pb isotopic data are acquired and reduced, including use of the precisely calibrated EARTHTIME 202Pb-205Pb-233U-235U tracer solution, changes in the isotopic compositions of standards used to calibrate the tracer, new error propagation algorithms, and improved data acquisition and reduction techniques (23⇓–25). Application of these improvements, discussed in detail in the SI Text, yields significantly improved accuracy and precision on the weighted mean and interpolated dates (Table S1 and Fig. S2). Uncertainties associated with weighted mean 206Pb/238U dates are reported as ±x/y/z, where x is the analytical (internal) uncertainties and y and z include the systematic uncertainties associated with tracer calibration (0.03%) and 238U decay constant (0.05%), respectively. If calculated dates are to be compared with other U-Pb laboratories not using the EARTHTIME tracer, then ±y should be used for each laboratory. If compared with other chronometers such as Ar-Ar or astrochronology, then ±z should be used.
The section at Meishan has long been recognized as being highly condensed, implying that calculated sediment accumulation rates may not accurately account for hiatuses between dated ash horizons (3). In agreement with Shen et al. (3), Jin et al. (26), and Wang et al. (27), we define the onset of extinction at the base of bed 25 and the end of the main extinction interval at bed 28. In other sections with higher accumulation rates such as Penglaitan (China) (3) and Gartnerkofel core (Swiss Alps) (28, 29), the extinction appears more abrupt; thus, our duration estimate between bed 25 and bed 28 at Meishan of 61 ± 48 ka is a maximum (Table 1, Fig. 1, and Fig. S1). This estimate is three times shorter than reported by Shen et al. (3) (Figs. S1 and S2) and is more consistent with the recent estimate derived from astrochronology of 83 ka (20).
The new geochronology permits a detailed examination of the relationships between the extinction, the isotopic composition of carbonate carbon, and its rate of change. The carbon isotope record is characterized by a negative shift in composition beginning just above the base of bed 23 (251.999 ± 0.039 Mya), 60 (−17/+56) ka before the beginning of the mass extinction interval, from +3–4‰ toward the lighter values (∼ −1‰) that characterize the earliest Triassic (Fig. 1 and Fig. S1). δ13C(carb) drops off rapidly in the upper 6 cm of bed 24e, from +2 to ∼ −4 ‰ (Fig. 1 and Fig. S1). In many sections that lack detailed paleontology and geochronology, this negative excursion is used to mark the onset of the extinction interval. The negative shift and subsequent rebound has a duration of between 2.1 and 18.8 ka depending on accumulation rate, slightly predating the beginning of the maximum extinction interval. Immediately following the initial large negative excursion, carbonate carbon isotopic composition oscillates (±1–2‰), until ∼1 m above bed 33. Using an accumulation rate derived from interpolating between dated beds 33 and 34, this period of oscillation lasts until 251.572 ± 0.069 Mya, a duration of 427 ± 79 ka. [Uncertainty on interpolated dates is calculated using a Monte Carlo simulation, which exploits stratigraphic superposition of dated rocks (30).] For the remainder of bed 34, ∼5 m of mudstone and micritic limestone, the δ13C(carb) is constant at ∼ −1‰, in contrast to the +3–4‰ that characterizes the preextinction interval (Fig. 1). The δ13C(carb) composition then rises gradually starting at ∼251.5, just below the dated ash within bed 34 and increases to the top of bed 39 where it is interrupted by a sharp, short-lived decrease calculated to have begun ∼251.4 Mya (Fig. 1). Above this perturbation, δ13C(carb) remains at ∼1‰ for the remainder of the Griesbachian. The negative excursion within bed 39 at Meishan was not recognized by Cao et al. (11) or Xie et al. (17) due to sample spacing, although more closely spaced sampling reported in Song et al. (31) recovers it. The second excursion is possibly correlative with one seen in the GK-1 core from the Carnic Alps (32), although this excursion cannot yet be confirmed as a global signal or to be useful in correlation. The total duration of volatility in the carbonate carbon record from the initial negative inflection within bed 23 to the relatively stable positive values in the top of bed 39 is a minimum of 500 ka. Significant changes in lithology above the dated bed suggest that the sediment accumulation rate calculated between beds 32 and 34 is likely not applicable for this interval. As such, the duration represented by this interval of rock is uncertain, and our estimated duration for the entire interval is a minimum (Fig. 1).
Sea surface paleotemperature increases ∼10 °C (∼23–33 °C) over the extinction interval (9, 33), beginning near the base of bed 25 and continuing into the early Triassic (Fig. S1). Sea surface temperatures are estimated to have reached ∼33 °C by bed 28, coinciding with the end of the mass extinction interval, and continued to rise until at least 251.583 ± 0.086 Mya (bed 33) (9). Calcium isotopic composition (δ44/40Ca‰ bulk earth) also varies over the extinction interval, and when coupled with apparent physiological selectivity of the extinction and an absence of reef builders in the early Triassic, these data have been interpreted to support rapid acidification of the surface ocean coincident with the mass extinction (4, 12).
Discussion
The efficacy of many proposed kill mechanisms, such as synchronous sea surface and atmospheric temperature increase, rapid rise in pCO2, and flooding of shelf areas with anoxic and euxinic waters, depends on rate of change and on precisely when they occur relative to the onset of extinction (9, 34, 35). For example, it is crucial to know whether the ∼10 °C increase in sea surface temperature close to the extinction interval slightly predates or postdates the onset of the mass extinction (9, 33) (Fig. S1). More detailed study of the relationship between temperature increase and extinction is needed from less condensed sections than Meishan to evaluate whether temperature leads or lags the extinction and the relationship between temperature rise and changes in the carbonate carbon isotopic record. Using the maximum extinction duration of ∼60 ka, this suggests an ∼1 °C increase per 6,000 y, comparable to the rate and magnitude of the increase at the Paleocene–Eocene Thermal Maximum (PETM) (36) and Pleistocene/Holocene postglacial warming (∼2 °C/5 ka) (37). Tracking with temperature increase is a negative shift in δ44/40Ca, interpreted as resulting in part from acidification of the ocean over this same interval and fluctuations in δ13C(carb), consistent with continued volatility in the carbon cycle after the initial spike toward lighter composition in the top of bed 24e (4, 12) (Fig. 1 and Fig. S1). Thus, in 80 ± 45 ka (base of bed 24e → base of bed 28), there was a short-lived episode of major light carbon addition to the oceans, a major mass extinction, a rapid, dramatic increase in marine and terrestrial temperature, isotopic and biological evidence for ocean acidification, and a major shift in δ13C(carb) composition from an average of approximately +3.5‰ in the late Permian to approximately −1‰ in the earliest Triassic, until 251.502 ± 0.056 Mya. The observation that the terrestrial and marine extinctions occurred simultaneously (3) and the suggestion that the sequence of extinction can be correlated with metabolic rate (1, 4) support the conclusion that rapidly elevated atmospheric pCO2 and ocean/atmosphere temperatures drove a combination of kill mechanisms. However, whether the temperature increase leads, is synchronous with, or postdates the extinction is not yet known with sufficient precision. Although recovery and diversification in Ammonoids began in the earliest Triassic, the broad effects of this short-lived extinction or ecological restructuring persist for 5–10 My after the main extinction interval, emphasizing the evolutionary irreversibility of the event (38⇓⇓–41).
Many have proposed that the end-Permian extinction was triggered by the eruption/intrusion of the Siberian Traps Large Igneous Province, which is hypothesized to have been of short (∼1–2 Ma) duration, to have occurred at approximately the same time as the extinction, and to have generated the large volume of volatiles via degassing of lavas and sediments required to drive such dramatic atmospheric and biotic response (2, 8, 42⇓⇓⇓–46). The end-Permian extinction event occurred suddenly and rapidly (61 ± 48 ka) in an interval much shorter than current estimates for the total duration of Siberian Traps magmatism, suggesting that, similar to the end-Triassic extinction event, a single pulse of magmatism may be the most critical for triggering dramatic environmental change (43, 47, 48). Current U-Pb and Ar-Ar constraints on the timing and tempo of Siberian Traps magmatism are less precise by an order of magnitude than our new constraints on the extinction. With current estimates, it can only be concluded that magmatism either overlaps with or postdates the extinction (43, 47, 49). Additionally, the potential for bias between chronometers and subtle differences in calculated dates generated by single or multiple laboratories using different U-Pb data acquisition and reduction protocols currently prohibits exploring the full details of a causal relationship.
Payne and Kump (45) and Song et al. (50) hypothesized that the large volatility in the carbon cycle that dominates the interval from the beginning of the Dienarian through the Spathian is distinct from the extinction interval represented at the GSSP and likely represents new injection of light carbon and global warming–driven anoxia related to continued activity of the Siberian Traps Large Igneous Province, ∼1 My after the extinction (Fig. 3). However, neither study had sufficient temporal control on the carbon isotope excursions or the age of Siberian volcanism to further evaluate this hypothesis. Here we demonstrate that upper bed 34 (mid-Griesbachian) is 251.495± 0.064, which predates the second negative excursion in δ13C(carb) at Meishan by ∼100 ka if the sediment accumulation rate between beds 32 and 34 is applied. The large positive oscillation in δ13C(carb) (+6) observed by Payne et al. (51) and Meyer et al. (52) in southern China begins at about the Griesbachian/Dienerian boundary and continues through the Dienerian before swinging to values of ∼ −2‰ near the top of the Smithian (Fig. 3). Existing geochronology on the lower Smithian includes an age of 251.22 ± 0.2 (53, 54) (Fig. 3). If correct, this requires the entire Dienerian and part of the Greisbachian to have been deposited in ∼300 ka. As shown here (see SI Text), use of the recalibrated EARTHTIME tracer and adoption of EARTHTIME protocols may in some cases result in dates that are 200–500 ka younger than those reported before adoption of these methods. Until dates from the Smithian are repeated using latest EARTHTIME protocols and tracer, it is probably unwise to combine/compare data. Nonetheless, we are confident that the large positive-negative oscillations that begin at the Griesbachian/Dienerian boundary are −251 Mya or younger, based on upward extrapolation from our dated ash beds at Meishan (Fig. 3). Thus, it is clear that these anomalies are ∼1 Mya or more younger than the main extinction event and separated from the volatility characterizing this interval by a relatively stable plateau of δ13C(carb) at ∼2‰ from the top of the second excursion in the Griesbachian to the Griesbachian/Dienerian boundary, whereupon the reservoir rises to values of 6‰ by the close of the Dinerian (Fig. 3). Meyer et al. (55) use coupled δ13Corg and δ13Ccarb records from South China to support a causal connection between carbon isotope stabilization and enhanced biotic recovery in Middle Triassic time. Whether postextinction carbon cycle dynamics are being driven by intrinsic (biological) or extrinsic (environmental) forces before final restructuring in the early Anisian (∼247 Mya) remains unclear. Integration of a calibrated carbon isotope record into the middle and late Triassic with an improved age model for the Siberian Traps will be required to evaluate this question.
When the end-Permian extinction is compared with other short-lived events such as the end-Triassic and end-Cretaceous extinctions and the PETM, we see in common, a short-lived perturbation of the carbon cycle followed by a rise in atmospheric pCO2 and temperature, evidence for ocean acidification, anoxia, and rapid extinction (10s of thousands of years) (48, 56, 57). Recovery or restructuring (±additional extinction) occurs over an ∼3-Mya timescale for the end-Triassic event and over as much as 5 Mya for the full recovery of marine fauna after the end Cretaceous, whereas full recovery from the end-Permian extinction may have taken as long as 5–10 My (39, 58, 59). A key question is whether the observed sequence is related to a common trigger such as volcanism or whether it less a function of mechanism and more fundamentally related to cascades of multiple feedbacks between intrinsic and extrinsic drivers. For all extinctions, understanding the timescales of extinction, the filling and restructuring of postextinction ecospace, and full recovery in different bathymetric and geographic settings is crucial and may be the key to understanding the truly singular and irreversible nature of the end-Permian event, as well as providing better context for the next millennium. We now have the geochronological tools to explore these feedbacks at the millennial to decamillenial timescales, which will in turn encourage higher-resolution chemostratigraphy to be obtained and allow detailed evaluation of more general models for mass extinction.
Conclusion
Our age model for the end-Permian extinction provides a precise and accurate timeline for the sequence of events at the end of the Permian, including carbon cycle reorganization, the main extinction event, a dramatic increase in global sea surface and atmospheric temperatures, possible ocean acidification, and a framework for exploring the cause and effects of the environmental changes and feedbacks that led to the greatest Phanerozoic mass extinction. The extinction had a duration of 61 ± 48 ka and was preceded by the onset of a rapid reorganization of the carbon cycle, including a rapid negative spike in δ13C(carb) of 3‰ lasting between 2.1 and 18.8 ka and a global shift in δ13C(carb) from approximately +4‰ το −1.5‰, with a duration of 427 ± 79 ka. This record represents a potentially characteristic δ13C(carb) topology for the end-Permian event, which will stimulate refined comparison with other Permian-Triassic sections, although the highly condensed nature of the Meishan section makes comparison with other sections difficult. The timing of the extinction and associated changes in environmental conditions are consistent with a very rapid biological response to environmental change followed by a complex recovery/restructuring period that took some 10 Ma for many species (38⇓⇓–41) and established the ecosystems that would dominate the Mesozoic. Further integration of the extinction timescale with detailed chemostratigraphic, cyclostratigraphic, and paleobiological data should allow many more insights into the dynamics and timing of extinction and restructuring. In addition, it is clear that more and higher precision geochronology from additional stratigraphic sections is needed. We predict that with further work will come the deconvolution of the end-Permian extinction into a cascade of smaller, shorter-lived extinction and recovery events, driven by differences in paleogeography, biology, and environmental degradation. The short-lived nature of the extinction, protracted nature of the recovery, and comparison with other extinction events suggests that environmental conditions preceding the largest of the Phanerozoic mass extinctions must have crossed a critical threshold or “tipping point” from which the biosphere was unable to recover or adapt quickly enough to survive.
Methods
Zircon crystals were separated from bulk samples using a combination of ultrasonic disaggregation and pulverization using a shatterbox, which was followed by magnetic separation, standard heavy-liquid separation, and careful selection of crystals under a microscope. Following thermal annealing at 900 °C for 60 h, each zircon crystal was placed into a 200-μL Teflon microcapsule and leached in 29 M hydrofluoric (HF) inside high-pressure Parr vessels held at 220 °C for 12 h, a procedure modified after the Chemical Abrasion partial-dissolution procedure of Mattinson (60). Grains were then transferred to 3-mL Savillex PFA beakers and rinsed with 16 M HNO3 and 6 M HCl and fluxed in the acid at 80 °C, followed by a 30-min ultrasonic bath. Between acid washes, grains were rinsed with Milli-Q water. Single zircon crystals were loaded with clean water into teflon microcapsules and spiked with the EARTHTIME 202Pb-205Pb-233U-235U (ET2535) tracer solution and dissolved in 29 M HF at 220 °C for 48 h. On dissolution, aliquots were dried down on a hotplate and redissolved under pressure in 6 M HCl overnight at 180 °C. Sample solutions were then dried and redissolved at 80 °C in 3 N HCl. Lead and U were separated using a miniaturized HCl-based ion-exchange chromatography procedure modified from Krogh (61) with 40-μL columns of AG1 × 8 resin. Eluted U and Pb were dried down with H3PO4 and then redissolved in a silica gel emitter solution (62) and loaded onto a zone-refined, outgassed Re filament. Measurement of isotopic ratios was done on an IsotopX ×62 multiple-collector thermal ionization mass spectrometer. Isotopes of Pb were measured by peak-hopping on a single Daly/photomultiplier detector system. Isotopes of U were measured as UO2 on Faraday detectors in static mode. Isotope ratios of U and Pb were corrected for mass fractionation during analysis using the ET2535 tracer solution. Data acquisition and reduction were done using the Tripoli and U-Pb Redux software packages (24, 25).
Acknowledgments
We thank Dan Rothman for discussions and two anonymous reviewers for comments. This work was supported by National Science Foundation Continental Dynamics Grant EAR-0807475 (to S.B.) and Instrumentation and Facilities EAR-0931839 (to S.B.). Additional support was provided from National Aeronautics and Space Administration Astrobiology Grant NNA08CN84A (to S.B.) and National Natural Science Foundation of China Grant 41290260 (to S.-z.S.). This work benefitted from the efforts of the EARTHTIME U-Pb community.
Footnotes
- ↵1To whom correspondence should be addressed. E-mail: sburgess{at}mit.edu.
Author contributions: S.D.B. performed research; S.B. and S.-z.S. collected samples; S.D.B. did the isotopic analyses; S.B. and S.-z.S. contributed to writing and data interpretation; and S.D.B. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
See Commentary on page 3203.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1317692111/-/DCSupplemental.
References
- ↵
- ↵
- ↵.
- Shen SZ,
- et al.
- ↵
- ↵
- ↵.
- Erwin DH
- ↵.
- Gattuso JP,
- Hansson L
- Knoll AH,
- Fischer WW
- ↵.
- Brand U,
- et al.
- ↵.
- Joachimski MM,
- et al.
- ↵
- ↵
- ↵.
- Hinojosa JL,
- et al.
- ↵.
- Bowring SA,
- Erwin DH,
- Davidek K,
- Wang W,
- Jin MW,
- Martin YG
- ↵
- ↵.
- Mundil R,
- Ludwig KR,
- Metcalfe I,
- Renne PR
- ↵.
- Twitchett RJ,
- Looy CV,
- Morante R,
- Visscher H,
- Wignall PB
- ↵.
- Xie S,
- et al.
- ↵.
- Huang C,
- Tong J,
- Hinnov L,
- Chen ZQ
- ↵.
- Rampino MR,
- Kaiho K
- ↵
- ↵.
- Wang Y,
- et al.
- ↵.
- Jin YG,
- et al.
- ↵
- ↵.
- McLean N,
- Bowring J,
- Bowring S
- ↵.
- Bowring JF,
- McLean NM,
- Bowring SA
- ↵.
- Jin YG,
- et al.
- ↵.
- Wang Y,
- et al.
- ↵.
- Rampino MR,
- Prokoph A,
- Adler A
- ↵.
- Gorjan P,
- Kaiho K,
- Chen ZQ
- ↵.
- Guex J,
- et al.
- ↵
- ↵.
- Korte C,
- et al.
- ↵.
- Sun Y,
- et al.
- ↵.
- Kiehl JT,
- Shields CA
- ↵
- ↵
- ↵.
- Lea DW,
- Pak DK,
- Spero HJ
- ↵.
- Payne JL,
- et al.
- ↵
- ↵.
- Song H,
- et al.
- ↵.
- Brayard A,
- et al.
- ↵
- ↵
- ↵
- ↵
- ↵.
- Renne PR,
- Black MT,
- Zichao Z,
- Richards MA,
- Basu AR
- ↵
- ↵.
- Blackburn TJ,
- et al.
- ↵
- ↵.
- Song H,
- Wignall PB,
- Tong J,
- Yin H
- ↵.
- Payne JL,
- et al.
- ↵
- ↵.
- Galfetti T,
- et al.
- ↵
- ↵
- ↵
- ↵.
- Schulte P,
- et al.
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Gradstein FM,
- Ogg JG,
- Schmitz MD,
- Ogg G
- Ogg JG
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Physical Sciences
- Earth, Atmospheric, and Planetary Sciences