As the climate changes and ecosystems shift toward novel combinations of species, the methods and metrics of conservation science are becoming less species-centric. To meet this growing need, marine conservation paleobiologists stand to benefit from the addition of new, taxon-free benthic indices to the live–dead analysis tool kit. These indices, which were developed to provide actionable, policy-specific data, can be applied to the readily preservable component of benthic communities (e.g., mollusks) to assess the ecological quality status of the entire community. Because these indices are taxon-free, they remain applicable even as the climate changes and novel communities develop—making them a potentially valuable complement to traditionally applied approaches for live–dead analysis, which tend to focus on maintaining specific combinations of species under relatively stable environmental conditions. Integrating geohistorical data with these established indices has potential to increase the salience of the live–dead approach in the eyes of resource managers and other stakeholders.

Simulations are playing an increasingly important role in paleobiology. When designing a simulation study, many decisions have to be made and common challenges will be encountered along the way. Here, we outline seven rules for executing a good simulation study. We cover topics including the choice of study question, the empirical data used as a basis for the study, statistical and methodological concerns, how to validate the study, and how to ensure it can be reproduced and extended by others. We hope that these rules and the accompanying examples will guide paleobiologists when using simulation tools to address fundamental questions about the evolution of life.

Estimating speciation and extinction rates is essential for understanding past and present biodiversity, but is challenging given the incompleteness of the rock and fossil records. Interest in this topic has led to a divergent suite of independent methods—paleontological estimates based on sampled stratigraphic ranges and phylogenetic estimates based on the observed branching times in a given phylogeny of living species. The fossilized birth–death (FBD) process is a model that explicitly recognizes that the branching events in a phylogenetic tree and sampled fossils were generated by the same underlying diversification process. A crucial advantage of this model is that it incorporates the possibility that some species may never be sampled. Here, we present an FBD model that estimates tree-wide diversification rates from stratigraphic range data when the underlying phylogeny of the fossil taxa may be unknown. The model can be applied when only occurrence data for taxonomically identified fossils are available, but still accounts for the incomplete phylogenetic structure of the data. We tested this new model using simulations and focused on how inferences are impacted by incomplete fossil recovery. We compared our approach with a phylogenetic model that does not incorporate incomplete species sampling and to three fossil-based alternatives for estimating diversification rates, including the widely implemented boundary-crosser and three-timer methods. The results of our simulations demonstrate that estimates under the FBD model are robust and more accurate than the alternative methods, particularly when fossil data are sparse, as the FBD model incorporates incomplete species sampling explicitly.

The taxonomic and ecologic composition of Earth's biota has shifted dramatically through geologic time, with some clades going extinct while others diversified. Here, we derive a metric that quantifies the change in biotic composition due to extinction or origination and show that it equals the product of extinction/origination magnitude and selectivity (variation in magnitude among groups). We also define metrics that describe the extent to which a recovery (1) reinforced or reversed the effects of extinction on biotic composition and (2) changed composition in ways uncorrelated with the extinction. To demonstrate the approach, we analyzed an updated compilation of stratigraphic ranges of marine animal genera. We show that mass extinctions were not more selective than background intervals at the phylum level; rather, they tended to drive greater taxonomic change due to their higher magnitudes. Mass extinctions did not represent a separate class of events with respect to either strength of selectivity or effect. Similar observations apply to origination during recoveries from mass extinctions, and on average, extinction and origination were similarly selective and drove similar amounts of biotic change. Elevated origination during recoveries drove bursts of compositional change that varied considerably in effect. In some cases, origination partially reversed the effects of extinction, returning the biota toward the pre-extinction composition; in others, it reinforced the effects of the extinction, magnifying biotic change. Recoveries were as important as extinction events in shaping the marine biota, and their selectivity deserves systematic study alongside that of extinction.

The typical marine animal has increased in biovolume by more than two orders of magnitude since the beginning of the Cambrian, but the causes of this trend remain unknown. We test the hypothesis that the efficiency of intra-organism oxygen delivery is a major constraint on body-size evolution in marine animals. To test this hypothesis, we compiled a dataset comprising 13,723 marine animal genera spanning the Phanerozoic. We coded each genus according to its respiratory medium, circulatory anatomy, and feeding mode. In extant genera, we find that respiratory medium and circulatory anatomy explain more of the difference in size than feeding modes. Likewise, we find that most of the Phanerozoic increase in mean biovolume is accounted for by size increase in taxa that accomplish oxygen delivery through closed circulatory systems. During the Cambrian, water-breathing animals with closed circulatory systems were smaller, on average, than contemporaries with open circulatory systems. However, genera with closed circulatory systems superseded in size genera with open circulatory systems by the Middle Ordovician, as part of their Phanerozoic-long trend of increasing size. In a regression analysis, respiratory and circulatory anatomy explain far more size variation in the living fauna than do feeding modes, even after accounting for taxonomic affinity at the class level. These findings suggest that ecological and environmental drivers of the Phanerozoic increase in the mean size of marine animals operated within strong, anatomically determined constraints.

The choice of measures used to estimate the richness of species, genera, or higher taxa is a crucial matter in paleobiology and ecology. This paper evaluates four methods called shareholder quorum subsampling, true richness estimated using a Poisson sampling model (TRiPS), squares, and the corrected first-order jackknife (cJ1). Quorum subsampling interpolates to produce a relative richness estimate, while the other three extrapolate to the size of the overall species pool. Here I use routine ecological data to show that squares and cJ1 pass several basic validation tests, but TRiPS does not. First, TRiPS estimates are insensitive to the shape of abundance distributions, being entirely predicted by total counts of species and of individuals regardless of the details. Furthermore, TRiPS tends not to extrapolate at all when sampling is moderate or intense. Second, all three extrapolators yield lower values when they work with small uniform subsamples of large raw inventories. The third test is a split-analyze-and-sum analysis: each inventory is divided between the most common and least common halves of the abundance distribution, the methods are applied to the half-inventories, and the estimates are summed. Squares and cJ1 perform well here, but TRiPS does not extrapolate as long as the full inventories are reasonably well-sampled. It is otherwise not particularly accurate. The extrapolators are largely insensitive to the influence of abundance distribution evenness, as quantified using Pielou's J and a new index called the ratio of means. Quorum subsampling generally performs well, but it stumbles on the split-analyze-and-sum test and is confounded somewhat by evenness.

The Cretaceous–Paleogene (K/Pg) extinction appears to have been geographically heterogeneous for some organismal groups. Southern Hemisphere K/Pg palynological records have shown lower extinction and faster recovery than in the Northern Hemisphere, but no comparable, well-constrained Southern Hemisphere macrofloras spanning this interval had been available. Here, macrofloral turnover patterns are addressed for the first time in the Southern Hemisphere, using more than 3500 dicot leaves from the latest Cretaceous (Maastrichtian) and the earliest Paleocene (Danian) of Argentine Patagonia. A maximum ca. 90% macrofloral extinction and ca. 45% drop in rarefied species richness is estimated across the K/Pg, consistent with substantial species-level extinction and previously observed extirpation of host-specialized leaf mines. However, prior palynological and taxonomic studies indicate low turnover of higher taxa and persistence of general floral composition in the same sections. High species extinction, decreased species richness, and homogeneous Danian macrofloras across time and facies resemble patterns often observed in North America, but there are several notable differences. When compared with boundary-spanning macrofloras at similar absolute paleolatitudes (ca. 50°S or 50°N) from the Williston Basin (WB) in the Dakotas, both Maastrichtian and Danian Patagonian species richnesses are higher, extending a history of elevated South American diversity into the Maastrichtian. Despite high species turnover, our analyses also reveal continuity and expansion of leaf morphospace, including an increase in lobed and toothed species unlike the Danian WB. Thus, both Patagonian and WB K/Pg macrofloras support a significant extinction event, but they may also reflect geographically heterogeneous diversity, extinction, and recovery patterns warranting future study.

Larger body size has long been assumed to correlate with greater risk of extinction, helping to shape body-size distributions across the tree of life, but a lack of comprehensive size data for fossil taxa has left this hypothesis untested for most higher taxa across the vast majority of evolutionary time. Here we assess the relationship between body size and extinction using a data set comprising the body sizes, stratigraphic ranges, and occurrence patterns of 9408 genera of fossil marine animals spanning eight Linnaean classes across the past 485 Myr. We find that preferential extinction of smaller-bodied genera within classes is substantially more common than expected due to chance and that there is little evidence for preferential extinction of larger-bodied genera. Using a capture–mark–recapture analysis, we find that this size bias of extinction persists even after accounting for a pervasive bias against the sampling of smaller-bodied genera within classes. The size bias in extinction also persists after including geographic range as an additional predictor of extinction, indicating that correlation between body size and geographic range does not provide a simple explanation for the association between size and extinction. Regardless of the underlying causes, the preferential extinction of smaller-bodied genera across many higher taxa and most of geologic time indicates that the selective loss of large-bodied animals is the exception, rather than the rule, in the evolution of marine animals.

One of the most remarkable differences between Paleogene penguins and their living relatives is the shape and length of their beaks. Many of the Eocene and Oligocene penguins have a thin and elongated spear-like bill, which contrasts with the proportionally shorter and more robust bill of most living species. These differences suggest an important shift in their feeding strategies. This study explores the morphological disparity on the skull of penguins, emphasizing bill morphology and it relationship with feeding habits. For this, the skulls of 118 species of aquatic birds, including 21 fossil and living penguins, were analyzed using two-dimensional geometric morphometric. The results show that, unlike what has been reported for modern birds overall, in penguins and Aequornithes, bill elongation is related to a reduction of the braincase. The discriminant analysis shows that there are significant differences between penguins that feed near or far from the coast and between those that consume nectonic and planktonic prey, identifying Madrynornis as the only extinct form with a possibly planktonic diet. Additionally, it is clear that Paleogene penguins occupy a region of morphospace unexplored by most diving birds, with the western grebe being their closest modern analogue. This is consistent with the hypothesis that giant penguins hunted by harpooning and not by biting as living forms do, signaling a significant change in the habits of those birds leading to the emergence of their crown group.

Principal coordinates analysis (PCoA) is a statistical ordination technique commonly applied to morphology-based cladistic matrices to study macroevolutionary patterns, morphospace occupation, and disparity. However, PCoA-based morphospaces are dissociated from the original data; therefore, whether such morphospaces accurately reflect body-plan disparity or extrinsic factors, such as body size, remains uncertain. We collated nine character–taxon matrices of dinosaurs together with body-mass estimates for all taxa and tested for relationships between body size and both the principal axis of variation (i.e., PCo1) and the entire set of PCo scores. The possible effects of body size on macroevolutionary hypotheses derived from ordinated matrices were tested by reevaluating evidence for the accelerated accumulation of avian-type traits indicated by a strong directional shift in PCo1 scores in hypothetical ancestors of modern birds. Body mass significantly accounted for, on average, approximately 50% and 16% of the phylogenetically corrected variance in PCo1 and all PCo scores, respectively. Along the avian stem lineage, approximately 30% of the morphological variation is attributed to the reconstructed body masses of each ancestor. When the effects of body size are adjusted, the period of accelerated trait accumulation is replaced by a more gradual, additive process. Our results indicate that even at low proportions of variance, body size can noticeably affect macroevolutionary hypotheses generated from ordinated morphospaces. Future studies should thoroughly explore the nature of their character data in association with PCoA-based morphospaces and use a residual/covariate approach to account for potential correlations with body size.

The gomphotheres were a diverse and widespread group of proboscideans occupying Eurasia, North America, and South America throughout the Neogene. Their decline was temporally and spatially heterogeneous, and the gomphotheres ultimately became extinct during the late Pleistocene; however, the genus Cuvieronius is rarely represented in late Pleistocene assemblages in North America. Two alternative hypotheses have been invoked to explain this phenomenon: (1) competitive exclusion by sympatric mammoths and mastodons or (2) ecologic displacement due to an environmental transition from closed forests to open grasslands. To test whether competition for resources contributed to the demise of North American Cuvieronius, we present herein a large collection of stable isotope and dental microwear data from populations occupying their Pleistocene refugium in the Atlantic Coastal Plain. Results suggest that Cuvieronius consumed a wide range of resources with variable textural and photosynthetic properties and was not specialized on either grasses or browse. Further, we document evidence for the consumption of similar foods between contemporaneous gomphotheres, mammoths, and mastodons. The generalist feeding strategy of the gomphotheres likely facilitated their high Miocene abundance and diversity. However, this “jack of all trades and master of none” feeding strategy may have proved challenging following the arrival of mammoths and likely contributed to the extirpation of Cuvieronius in North America.

Quantifying true patterns of biodiversity change over the Cenozoic has major implications for all of biology and paleontology but is still a source of significant debate. The problem centers on the magnitude and nature of several well-known sampling effects and analytical biases in the fossil record, including the Pull of the Recent. We test the effect of the Pull of the Recent at both generic and species levels on the exemplary New Zealand Cenozoic marine mollusk fossil record. We examine several biological traits of species to determine whether particular attributes of taxa control their likely presence or absence in the youngest fossil record (<2.4 Ma). We demonstrate that, for a tectonically active region, the Pull of the Recent does not exert a strong effect on apparent diversity patterns of genera and species over the Cenozoic at temporal scales typically used in global and regional biodiversity analyses. This result agrees with previous studies quantifying the effect of the Pull of the Recent in the marine and terrestrial realms at the genus level. The effect of the Pull of the Recent, although small, is greatest for the youngest fossil record (<2.4 Ma), particularly for species. This increase cannot easily be explained by effects related to shell mineralogical composition, size, habitat, taxonomic class, or lithification. The small effect that the Pull of the Recent exerts on the New Zealand molluscan fossil record implies that the apparent rise in regional marine diversity during the Cenozoic represents a true biological signal and/or reflects other confounding effects not considered here.

The Red Queen's hypothesis portrays evolution as a never-ending competition for expansive energy, where one species’ gain is another species’ loss. The Red Queen is neutral with respect to body size, implying that neither small nor large species have a universal competitive advantage. Here we ask whether, and if so how, the Red Queen's hypothesis really can accommodate differences in body size. The maximum population growth in ecology clearly depends on body size—the smaller the species, the shorter the generation length, and the faster it can expand given sufficient opportunity. On the other hand, large species are more efficient in energy use due to metabolic scaling and can maintain more biomass with the same energy. The advantage of shorter generation makes a wide range of body sizes competitive, yet large species do not take over. We analytically show that individuals consume energy and reproduce in physiological time, but need to compete for energy in real time. The Red Queen, through adaptive evolution of populations, balances the pressures of real and physiological time. Modeling competition for energy as a proportional prize contest from economics, we further show that Red Queen's zero-sum game can generate unimodal hat-like patterns of species rise and decline that can be neutral in relation to body size.

The Cenozoic genus Terebratula seems to be an exception to the post-Permian trend in brachiopod retreat to offshore habitats, because it was species rich and numerically abundant in warm-temperate shallow-water environments in the Mediterranean and the Paratethys realms. This was so despite the general dominance of bivalves and the pervasive bioturbation and predation pressure during the Neogene. Terebratula, however, went extinct in the Calabrian (Pleistocene). The optimal environmental conditions for Terebratula during its prime are poorly known. The Águilas Basin (SE Spain) is an ideal study area to investigate the habitat of Terebratula, because shell beds of this brachiopod occur there cyclically in early Pliocene deposits. We evaluate the paleoecological boundary conditions controlling the distribution of Terebratula by estimating its environmental tolerances using benthic and planktic foraminiferal and nannoplanktic assemblages and oxygen isotopes of the secondary layer brachiopod calcite. Our results suggest that Terebratula in the Águilas Basin favored oligotrophic to mesotrophic, well-oxygenated environments at water depths of 60–90 m. Planktic foraminiferal assemblages and oxygen isotopes point to sea-surface temperatures between ~16°C and 22°C, and bottom-water temperatures between 17°C and 24°C. The analyzed proxies indicate that Terebratula tolerated local variations in water depth, bottom temperature, oxygenation, productivity, and organic enrichment. Terebratula was probably excluded by grazing pressure from well-lit environments and preferentially occupied sediment-starved, current-swept upper offshore habitats where coralline red algae were absent. Narrow temperature ranges of Terebratula species might have been a disadvantage during the high-amplitude seawater temperature fluctuations that started about 1 Ma, when the genus went extinct.

A thorough understanding of how communities respond to extreme changes, such as biotic invasions, is essential to manage ecosystems today. Here we constructed fossil food webs to identify changes in Late Ordovician (Katian) shallow-marine paleocommunity structure and functioning before and after the Richmondian invasion, a well-documented ancient invasion. Food webs were compared using descriptive metrics and cascading extinction on graphs models. Richness at intermediate trophic levels was underrepresented when using only data from the Paleobiology Database relative to museum collections, resulting in a spurious decrease in modeled paleocommunity stability. Therefore, museum collections and field sampling may provide more reliable sources of data for the reconstruction of trophic organization in comparison to online data repositories. The invasion resulted in several changes in ecosystem dynamics. Despite topological similarities between pre- and postinvasion food webs, species loss occurred corresponding to a minor decrease in functional groups. Invaders occupied all of the preinvasion functional guilds, with the exception of four incumbent guilds that were lost and one new guild, corroborating the notion that invaders replace incumbents and fill preexisting niche space. Overall, models exhibited strong resistance to secondary extinction, although the postinvasion community had a lower threshold of collapse and more variable response to perturbation. We interpret these changes in dynamics as a decrease in stability, despite similarities in overall structure. Changes in food web structure and functioning resulting from the invasion suggest that conservation efforts may need to focus on preserving functional diversity if more diverse ecosystems are not inherently more stable.

The middle (Wuliuan Stage) Cambrian Burgess Shale is famous for its exceptional preservation of diverse and abundant soft-bodied animals through the “thick” Stephen Formation. However, with the exception of the Walcott Quarry (Fossil Ridge) and the stratigraphically older Tulip Beds (Mount Stephen), which are both in Yoho National Park (British Columbia), quantitative assessments of the Burgess Shale have remained limited. Here we first provide a detailed quantitative overview of the diversity and structure of the Marble Canyon Burgess Shale locality based on 16,438 specimens. Located 40 km southeast of the Walcott Quarry in Kootenay National Park (British Columbia), Marble Canyon represents the youngest site of the “thick” Stephen Formation. We then combine paleoecological data sets from Marble Canyon, Walcott Quarry, Tulip Beds, and Raymond Quarry, which lies approximately 20 m directly above the Walcott Quarry, to yield a combined species abundance data set of 77,179 specimens encompassing 234 species-level taxa. Marble Canyon shows significant temporal changes in both taxonomic and ecological groups, suggesting periods of stasis followed by rapid turnover patterns at local and short temporal scales. At wider geographic and temporal scales, the different Burgess Shale sites occupy distinct areas in multivariate space. Overall, this suggests that the Burgess Shale paleocommunity is far patchier than previously thought and varies at both local and regional scales through the “thick” Stephen Formation. This underscores that our understanding of Cambrian diversity and ecological networks, particularly in early animal ecosystems, remains limited and highly dependent on new discoveries.

The Order Spiriferinida spanning the latest Ordovician to Early Jurassic is a small group of brachiopods overshadowed by other taxon-rich clades during the Paleozoic. It diversified significantly after the end-Permian extinction and became one of the four major clades of Triassic brachiopods. However, the phylogeny and recovery dynamics of this clade during the Triassic still remain unknown. Here, we present a higher-level parsimony-based phylogenetic analysis of Mesozoic spiriferinids to reveal their evolutionary relationships. Ecologically related characters are analyzed to indicate the variances in ecomorphospace occupation and disparity of spiriferinids through the Permian–Triassic (P-Tr) transition. For comparison with potential competitors of the spiriferinids, the pre-extinction spiriferids are also included in the analysis. Phylogenetic trees demonstrate that about half of the Mesozoic families appeared during the Anisian, indicating the greatest phylogenetic diversification at that time. Triassic spiriferinids reoccupied a large part of the ecomorphospace released by its competitor spiriferids during the end-Permian extinction; they also fully exploited the cyrtiniform region and developed novel lifestyles. Ecomorphologic disparity of the spiriferinids dropped greatly in the Early Triassic, but it rebounded rapidly and reached the level attained by the pre-extinction spiriferids in the Late Triassic. The replacement in ecomorphospace occupation between spiriferids and spiriferinids during the P-Tr transition clearly indicates that the empty ecomorphospace released by the extinction of Permian spiriferids was one of the important drivers for the diversification of the Triassic spiriferinids. The Spiriferinida took over the empty ecomorphospace and had the opportunity to flourish.

The classical taxonomy of fossil invertebrates is based on subjective judgments of morphology, which can cause confusion, because there are no codified standards for the classification of genera. Here, we explore the validity of the genus taxonomy of 75 species and morphospecies of the Follicucullidae, a late Paleozoic family of radiolarians, using a new method, Hayashi's quantification theory II (HQT-II), a general multivariate statistical method for categorical datasets relevant to discriminant analysis. We identify a scheme of 10 genera rather than the currently accepted 3 genera (Follicucullus, Ishigaconus, and Parafollicucullus). As HQT-II cannot incorporate stratigraphic data, a phylogenetic tree of Follicucullidae was reconstructed for 38 species using maximum parsimony. Six lineages emerged, roughly in concordance with the results of HQT-II. Combined with parsimony ancestral state reconstruction, the ancestral group of this family is Haplodiacanthus. Five other groups were discriminated, the Parafollicucullus, Curvalbaillella, Pseudoalbaillella, Longtanella, and Follicucullus–Cariver lineages. The morphological evolution of these lineages comprises a minimum essential list of eight states of four traits. HQT-II is a novel discriminant analytical multivariate method that may be of value in other taxonomic problems of paleobiology.

Coleoid cephalopods exhibited two distinct reproductive strategies, resulting in small pelagic and large demersal hatchlings, both in the geologic past and recently. In ectocochleate cephalopods, the hatching event is recorded in shell structures (e.g., nepionic constrictions, ultrastructural shifts, or ornamentation differences). In contrast, well-defined hatching markers do not exist on coleoid shells. Changes in septal spacing may be evidence of hatching (e.g., some extant sepiids), but not in all fossil groups. In the present study, we subdivide the early ontogenetic shells of phragmocone-bearing coleoids (belemnoids, spirulids, and sepiids) into key architectural stages and describe their reference to the hatching event. Belemnoids exhibit three key stages, the second of which is here considered to occur shortly before or after hatching. In spirulids and sepiids, there is only one key stage. In Mesozoic belemnoids, spirulids, and sepiids, hatching accordingly occurred with a total shell length of less than 2 mm, which corresponds to mantle lengths of small planktonic hatchlings. Production of small pelagic hatchlings and thus small eggs was therefore the dominant reproductive strategy within the Coleoidea. The first evidence of enlarged hatchlings appeared during the Maastrichtian in Groenlandibelus. During the Eocene, the large-egg strategy apparently became more widespread, particularly in belosaepiids.

Marine bryozoans have been members of benthic skeletal faunas since the Ordovician. These small suspension feeders collect particles in the range of 10 to 100 μm. Specific details of their feeding depend on the morphology of the feeding apparatus, which may be reflected in skeletal characters. While several studies have described the link between the skeletal and soft-body traits of gymnolaemate bryozoans, stenolaemates have received less attention. To fill this gap, we conducted a detailed analysis of morphometry within and across species and attempted to develop robust predictive models that can be used to infer the soft-body morphology from skeletal data. This, in turn, will help with extracting data on ecology of Paleozoic communities of suspension feeders from the extensive bryozoan fossil record. Characters of polypide morphology among New Zealand cyclostomates (single Recent order in Stenolaemata) displayed staggering variability and almost without exception were not connected to skeletal characters at the species level. When this variability is reduced to its central tendency, interspecific trends are more apparent. The relationship is positive, linear, and moderately strong, but the resulting models have wide predictive intervals (plus/minus hundreds of micrometers). A precise estimate of the characters of the feeding apparatus of modern, and especially fossil, stenolaemates may be difficult to attain, at least on the basis of the skeletal traits used here.