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.

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.

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 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 final 10 Myr of the Paleozoic saw two of the biggest biological crises in Earth history: the middlePermian extinction (often termed the Guadalupian–Lopingian extinction [GLE]) that was followed 7–8 Myr later by Earth's most catastrophic loss of diversity, the Permian–Triassic mass extinction (PTME). These crises are not only manifest as sharp decreases in biodiversity and—particularly for the PTME—total ecosystem collapse, but they also drove major changes in biological morphological characteristics such as the Lilliput effect. The evolution of test size among different clades of foraminifera during these two extinction events has been less studied. We analyzed a global database of foraminiferal test size (volume) including 20,226 specimens in 464 genera, 98 families, and 9 suborders from 632 publications. Our analyses reveal significant reductions in foraminiferal mean test size across the Guadalupian/Lopingian boundary (GLB) and the Permian/Triassic boundary (PTB), from 8.89 to 7.60 log10 μm3 (lg μm3) and from 7.25 to 5.82 lg μm3, respectively. The decline in test size across the GLB is a function of preferential extinction of genera exhibiting gigantism such as fusulinoidean fusulinids. Other clades show little change in size across the GLB. In contrast, all Lopingian suborders in our analysis (Fusulinina, Lagenina, Miliolina, and Textulariina) experienced a significant decrease in test size across the PTB, mainly due to size-biased extinction and within-lineage change. The PTME was clearly a major catastrophe that affected many groups simultaneously, and the GLE was more selective, perhaps hinting at a subtler, less extreme driver than the later PTME.

Cosmopolitanism occurred recurrently during the geologic past, especially after mass extinctions, but the underlying mechanisms remain poorly known. Three theoretical models, not mutually exclusive, can lead to cosmopolitanism: (1) selective extinction in endemic taxa, (2) endemic taxa becoming cosmopolitan after the extinction and (3) an increase in the number of newly originated cosmopolitan taxa after extinction. We analyzed an updated occurrence dataset including 831 middle Permian to Middle Triassic ammonoid genera and used two network methods to distinguish major episodes of ammonoid cosmopolitanism during this time interval. Then, we tested the three proposed models in these case studies. Our results confirm that at least two remarkable cosmopolitanism events occurred after the Permian–Triassic and late Smithian (Early Triassic) extinctions, respectively. Partitioned analyses of survivors and newcomers revealed that the immediate cosmopolitanism event (Griesbachian) after the Permian–Triassic event can be attributed to endemic genera becoming cosmopolitan (model 2) and an increase in the number of newly originated cosmopolitan genera after the extinction (model 3). Late Smithian cosmopolitanism is caused by selective extinction in endemic taxa (model 1) and an increase in the number of newly originated cosmopolitan genera (model 3). We found that the survivors of the Permian–Triassic mass extinction did not show a wider geographic range, suggesting that this mass extinction is nonselective among the biogeographic ranges, while late Smithian survivors exhibit a wide geographic range, indicating selective survivorship among cosmopolitan genera. These successive cosmopolitanism events during severe extinctions are associated with marked environmental upheavals such as rapid climate changes and oceanic anoxic events, suggesting that environmental fluctuations play a significant role in cosmopolitanism.

Throughout their 250 Myr history, archosaurian reptiles have exhibited a wide array of body sizes, shapes, and locomotor habits, especially in regard to terrestriality. These features make Archosauria a useful clade with which to study the interplay between body size, shape, and locomotor behavior, and how this interplay may have influenced locomotor evolution. Here, digital volumetric models of 80 taxa are used to explore how mass properties and body proportions relate to each other and locomotor posture in archosaurs. One-way, nonparametric, multivariate analysis of variance, based on the results of principal components analysis, shows that bipedal and quadrupedal archosaurs are largely distinguished from each other on the basis of just four anatomical parameters (p < 0.001): mass, center of mass position, and relative forelimb and hindlimb lengths. This facilitates the development of a quantitative predictive framework that can help assess gross locomotor posture in understudied or controversial taxa, such as the crocodile-line Batrachotomus (predicted quadruped) and Postosuchus (predicted biped). Compared with quadrupedal archosaurs, bipedal species tend to have relatively longer hindlimbs and a more caudally positioned whole-body center of mass, and collectively exhibit greater variance in forelimb lengths. These patterns are interpreted to reflect differing biomechanical constraints acting on the archosaurian Bauplan in bipedal versus quadrupedal groups, which may have shaped the evolutionary histories of their respective members.

Reconstructing patterns of macroevolution has become a central endeavor in paleobiology, because it offers insight into evolutionary models shaping the history of life. As the most diverse and abundant animals since the Cambrian period, arthropods provide copious data to elucidate the emergence of body plans in metazoan lineages. However, information provided by fossils on the tempo and mode of this phenomenon has lacked a recent synthesis. Here, I investigate macroevolutionary patterns of morphological evolution in Euarthropoda using a combined extinct and extant dataset optimized for multivariate analyses. Overall ordination patterns between the main morphogroups are consistent with another, independently coded, extant-only dataset providing molecular and morphological rates of evolution. Based on a “deep split” phylogenetic framework, total-group Mandibulata and Arachnomorpha emerge as directional morphoanatomical lineages, with basal fossil morphogroups showing heterogeneously spread-out occupations of the morphospace. In addition to a more homogeneous morphological variation, new morphogroups arose by successive reductions of translation distances; this pattern was interrupted only by terrestrialization events and the origin of pancrustaceans. A displaced optimum type of model is proposed to explain the fast assembly of canalized body plans during the Cambrian, with basal fossil morphogroups fitting intermediate fitness peaks in a moving adaptive landscape. Given time constraints imposed by the paleontological evidence, and owing to the interplay between canalization and modularity, as well as a decoupling between molecular and morphological rates, the rise of euarthropods would support the view that the swiftness of the Cambrian explosion was mostly associated with the buildup of genetic regulatory networks.