Biogeographic and bathymetric determinants of brachiopod extinction and survival during the Late Ordovician mass extinction

Reduced strength and increased variability of extinction selectivity during mass extinctions

Two of the traits most often observed to correlate with extinction risk in marine animals are geographical range and body size. However, the relative effects of these two traits on extinction risk have not been investigated systematically for either background times or during mass extinctions. To close this knowledge gap, we measure and compare extinction selectivity of geographical range and body size of genera within five classes of benthic marine animals across the Phanerozoic using capture-mark-recapture models. During background intervals, narrow geographical range is strongly associated with greater extinction probability, whereas smaller body size is more weakly associated with greater extinction probability. During mass extinctions, the association between geographical range and extinction probability is reduced in every class and fully eliminated in some, whereas the association between body size and extinction probability varies in strength and direction across classes. While geographical range is universally the stronger predictor of survival during background intervals, variation among classes during mass extinction suggests a fundamental shift in extinction processes during these global catastrophes.

Reduction in animal abundance and oxygen availability during and after the end-Triassic mass extinction

The end-Triassic biodiversity crisis was one of the most severe mass extinctions in the history of animal life. However, the extent to which the loss of taxonomic diversity was coupled with a reduction in organismal abundance remains to be quantified. Further, the temporal relationship between organismal abundance and local marine redox conditions is lacking in carbonate sections. To address these questions, we measured skeletal grain abundance in shallow-marine limestones by point counting 293 thin sections from four stratigraphic sections across the Triassic/Jurassic boundary in the Lombardy Basin and Apennine Platform of western Tethys. Skeletal abundance decreased abruptly across the Triassic/Jurassic boundary in all stratigraphic sections. The abundance of skeletal organisms remained low throughout the lower-middle Hettangian strata and began to rebound during the late Hettangian and early Sinemurian. A two-way ANOVA indicates that sample age (p < .01, η²  = 0.30) explains more of the variation in skeletal abundance than the depositional environment or paleobathymetry (p < .01, η²  = 0.15). Measured I/Ca ratios, a proxy for local shallow-marine redox conditions, show this same pattern with the lowest I/Ca ratios occurring in the early Hettangian. The close correspondence between oceanic water column oxygen levels and skeletal abundance indicates a connection between redox conditions and benthic organismal abundance across the Triassic/Jurassic boundary. These findings indicate that the end-Triassic mass extinction reduced not only the biodiversity but also the carrying capacity for skeletal organisms in early Hettangian ecosystems, adding to evidence that mass extinction of species generally leads to mass rarity among survivors.

The role of continental shelf bathymetry in shaping marine range shifts in the face of climate change

As a consequence of anthropogenic climate change, marine species on continental shelves around the world are rapidly shifting deeper and poleward. However, whether these shifts deeper and poleward will allow species to access more, less, or equivalent amounts of continental shelf area and associated critical habitats remains unclear. By examining the proportion of seabed area at a range of depths for each large marine ecosystem (LME), we found that shelf area declined monotonically for 19% of LMEs examined. However, the majority exhibited a greater proportion of shelf area in mid-depths or across several depth ranges. By comparing continental shelf area across 2° latitudinal bands, we found that all coastlines exhibit multiple instances of shelf area expansion and contraction, which have the potential to promote or restrict poleward movement of marine species. Along most coastlines, overall shelf habitat increases or exhibits no significant change moving towards the poles. The exception is the Southern West Pacific, which experiences an overall loss of area with increasing latitude. Changes in continental shelf area availability across latitudes and depths are likely to affect the number of species local ecosystems can support. These geometric analyses help identify regions of conservation priority and ecological communities most likely to face attrition or expansion due to variations in available area.

Accelerated mass extinction in an isolated biota during Late Devonian climate changes

The fossil record can illuminate factors that contribute to extinction risk during times of global environmental disturbance; for example, inferred thermal tolerance was an important predictor of extinction during several mass extinctions that corresponded with climate change. Additionally, members of geographically isolated biotas may face higher risk because they have less opportunity to migrate to suitable climate refugia during environmental disturbances. Here, we investigate how different types of risk intersect in the well-preserved brachiopod fauna of the Appalachian Foreland Basin during the two pulses of the Frasnian-Famennian mass extinction (Late Devonian, ~ 372 Ma). The selectivity of extinction is consistent with climate change (cooling) as a primary kill mechanism in this fauna. Overall, the extinction was mild relative to other regions, despite the many endemic species. However, vulnerable taxa went extinct more rapidly, during the first extinction pulse, such that the second pulse was insignificant. These results suggest that vulnerable taxa in geographically isolated biotas face heightened extinction risk at the initiation of environmental stress, but that taxa in other regions may eventually see elevated extinction risk if environmental stress repeats or intensifies.

Selectivity and the effect of mass extinctions on disparity and functional ecology

Selectivity of mass extinctions is thought to play a major role in coupling or decoupling of taxonomic, morphological, and ecological diversity, yet these measures have never been jointly evaluated within a single clade over multiple mass extinctions. We investigate extinction selectivity and changes in taxonomic diversity, morphological disparity, and functional ecology over the ~160-million-year evolutionary history of diplobathrid crinoids (Echinodermata), which spans two mass extinctions. Whereas previous studies documented extinction selectivity for crinoids during background extinction, we find no evidence for selectivity during mass extinctions. Despite no evidence for extinction selectivity, disparity remains strongly correlated with richness over extinction events, contradicting expected patterns of disparity given nonselective extinction. Results indicate that (i) disparity and richness can remain coupled across extinctions even when selective extinction does not occur, (ii) simultaneous decreases in taxonomic diversity and disparity are insufficient evidence for extinction selectivity, and (iii) selectivity differs between background and mass extinction regimes.

Persistent global marine euxinia in the early Silurian

The second pulse of the Late Ordovician mass extinction occurred around the Hirnantian-Rhuddanian boundary (~444 Ma) and has been correlated with expanded marine anoxia lasting into the earliest Silurian. Characterization of the Hirnantian ocean anoxic event has focused on the onset of anoxia, with global reconstructions based on carbonate δ²³⁸U modeling. However, there have been limited attempts to quantify uncertainty in metal isotope mass balance approaches. Here, we probabilistically evaluate coupled metal isotopes and sedimentary archives to increase constraint. We present iron speciation, metal concentration, δ⁹⁸Mo and δ²³⁸U measurements of Rhuddanian black shales from the Murzuq Basin, Libya. We evaluate these data (and published carbonate δ²³⁸U data) with a coupled stochastic mass balance model. Combined statistical analysis of metal isotopes and sedimentary sinks provides uncertainty-bounded constraints on the intensity of Hirnantian-Rhuddanian euxinia. This work extends the duration of anoxia to >3 Myrs – notably longer than well-studied Mesozoic ocean anoxic events.

The Late Ordovician mass extinction has been attributed to extended marine anoxia. Here, the authors use a metal isotope mass balance model and find the marine anoxic event lasted over 3 million years, notably longer than the anoxic event associated with the Permian-Triassic extinction and Cretaceous ocean anoxic events.

The evolutionary dynamics of the early Palaeozoic marine biodiversity accumulation

The early Palaeozoic Era records the initial biodiversification of the Phanerozoic. The increase in biodiversity involved drastic changes in taxon longevity, and in rates of origination and extinction. Here, we calculate these variables in unprecedented temporal resolution. We find that highly volatile origination and extinction rates are associated with short genus longevities during the Cambrian Period. During the Ordovician and Silurian periods, evolutionary rates were less volatile and genera persisted for increasingly longer intervals. The 90%-genus life expectancy doubled from 5 Myr in the late Cambrian to more than 10 Myr in the Ordovician–Silurian periods. Intervals with widespread ecosystem disruption are associated with short genus longevities during the Cambrian and with exceptionally high longevities during the Ordovician and Silurian periods. The post-Cambrian increase in persistence of genera, therefore, indicates an elevated ability of the changing early Palaeozoic marine ecosystems to sustainably maintain existing genera. This is evidence of a new level of ecosystem resilience which evolved during the Ordovician Period.

Cascading trend of Early Paleozoic marine radiations paused by Late Ordovician extinctions

The first 120 million years of Phanerozoic life witnessed significant changes in biodiversity levels. Attempts to correlate these changes to potential short-term environmental drivers have been hampered by the crude temporal resolution of current biodiversity estimates. We present a biodiversity curve for the Early Paleozoic with high temporal precision. It shows that once equatorial sea-surface temperatures fell to present-day levels during the early Mid Ordovician, marine biodiversity accumulation accelerated dramatically. However, this acceleration ceased as increased volcanism commenced during the mid-Late Ordovician. Since biodiversity levels were not restored for at least ∼35 million years, this finding redefines the nature of the end Ordovician mass extinctions and further reframes the Silurian as a prolonged recovery interval.

The greatest relative changes in marine biodiversity accumulation occurred during the Early Paleozoic. The precision of temporal constraints on these changes is crude, hampering our understanding of their timing, duration, and links to causal mechanisms. We match fossil occurrence data to their lithostratigraphical ranges in the Paleobiology Database and correlate this inferred taxon range to a constructed set of biostratigraphically defined high-resolution time slices. In addition, we apply capture–recapture modeling approaches to calculate a biodiversity curve that also considers taphonomy and sampling biases with four times better resolution of previous estimates. Our method reveals a stepwise biodiversity increase with distinct Cambrian and Ordovician radiation events that are clearly separated by a 50-million-year-long period of slow biodiversity accumulation. The Ordovician Radiation is confined to a 15-million-year phase after which the Late Ordovician extinctions lowered generic richness and further delayed a biodiversity rebound by at least 35 million years. Based on a first-differences approach on potential abiotic drivers controlling richness, we find an overall correlation with oxygen levels, with temperature also exhibiting a coordinated trend once equatorial sea surface temperatures fell to present-day levels during the Middle Ordovician Darriwilian Age. Contrary to the traditional view of the Late Ordovician extinctions, our study suggests a protracted crisis interval linked to intense volcanism during the middle Late Ordovician Katian Age. As richness levels did not return to prior levels during the Silurian—a time of continental amalgamation—we further argue that plate tectonics exerted an overarching control on biodiversity accumulation.

Identifying the most surprising victims of mass extinction events: an example using Late Ordovician brachiopods

Mass extinction events are recognized by increases in extinction rate and magnitude and, often, by changes in the selectivity of extinction. When considering the selective fingerprint of a particular event, not all taxon extinctions are equally informative: some would be expected even under a ‘background’ selectivity regime, whereas others would not and thus require special explanation. When evaluating possible drivers for the extinction event, the latter group is of particular interest. Here, we introduce a simple method for identifying these most surprising victims of extinction events by training models on background extinction intervals and using these models to make per-taxon assessments of ‘expected’ risk during the extinction interval. As an example, we examine brachiopod genus extinctions during the Late Ordovician Mass Extinction and show that extinction of genera in the deep-water ‘Foliomena fauna’ was particularly unexpected given preceding Late Ordovician extinction patterns.

Graptolite community responses to global climate change and the Late Ordovician mass extinction

Mass extinctions disrupt ecological communities. Although climate changes produce stress in ecological communities, few paleobiological studies have systematically addressed the impact of global climate changes on the fine details of community structure with a view to understanding how changes in community structure presage, or even cause, biodiversity decline during mass extinctions. Based on a novel Bayesian approach to biotope assessment, we present a study of changes in species abundance distribution patterns of macroplanktonic graptolite faunas (∼447-444 Ma) leading into the Late Ordovician mass extinction. Communities at two contrasting sites exhibit significant decreases in complexity and evenness as a consequence of the preferential decline in abundance of dysaerobic zone specialist species. The observed changes in community complexity and evenness commenced well before the dramatic population depletions that mark the tipping point of the extinction event. Initially, community changes tracked changes in the oceanic water masses, but these relations broke down during the onset of mass extinction. Environmental isotope and biomarker data suggest that sea surface temperature and nutrient cycling in the paleotropical oceans changed sharply during the latest Katian time, with consequent changes in the extent of the oxygen minimum zone and phytoplankton community composition. Although many impacted species persisted in ephemeral populations, increased extinction risk selectively depleted the diversity of paleotropical graptolite species during the latest Katian and early Hirnantian. The effects of long-term climate change on habitats can thus degrade populations in ways that cascade through communities, with effects that culminate in mass extinction.