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.

Echinoderms make up a substantial component of Ordovician marine invertebrates, yet their speciation and dispersal history as inferred within a rigorous phylogenetic and statistical framework is lacking. We use biogeographic stochastic mapping (BSM; implemented in the R package BioGeoBEARS) to infer ancestral area relationships and the number and type of dispersal events through the Ordovician for diploporan blastozoans and related species. The BSM analysis was divided into three time slices to analyze how dispersal paths changed before and during the great Ordovician biodiversification event (GOBE) and within the Late Ordovician mass extinction intervals. The best-fit biogeographic model incorporated jump dispersal, indicating this was an important speciation strategy. Reconstructed areas within the phylogeny indicate the first diploporan blastozoans likely originated within Baltica or Gondwana. Dispersal, jump dispersal, and sympatry dominated the BSM inference through the Ordovician, while dispersal paths varied in time. Long-distance dispersal events in the Early Ordovician indicate distance was not a significant predictor of dispersal, whereas increased dispersal events between Baltica and Laurentia are apparent during the GOBE, indicating these areas were important to blastozoan speciation. During the Late Ordovician, there is an increase in dispersal events among all paleocontinents. The drivers of dispersal are attributed to oceanic and epicontinental currents. Speciation events plotted against geochemical data indicate that blastozoans may not have responded to climate cooling events and other geochemical perturbations, but additional data will continue to shed light on the drivers of early Paleozoic blastozoan speciation and dispersal patterns.

The occupation of new environments by evolutionary lineages is frequently associated with morphological changes. This covariation of ecotype and phenotype is expected due to the process of natural selection, whereby environmental pressures lead to the proliferation of morphological variants that are a better fit for the prevailing abiotic conditions. One primary mechanism by which phenotypic variants are known to arise is through changes in the timing or duration of organismal development resulting in alterations to adult morphology, a process known as heterochrony. While numerous studies have demonstrated heterochronic trends in association with environmental gradients, few have done so within a phylogenetic context. Understanding species interrelationships is necessary to determine whether morphological change is due to heterochronic processes; however, research is hampered by the lack of a quantitative metric with which to assess the degree of heterochronic traits expressed within and among species. Here I present a new metric for quantifying heterochronic change, expressed as a heterochronic weighting, and apply it to xiphosuran chelicerates within a phylogenetic context to reveal concerted independent heterochronic trends. These trends correlate with shifts in environmental occupation from marine to nonmarine habitats, resulting in a macroevolutionary ratchet. Critically, the distribution of heterochronic weightings among species shows evidence of being influenced by both historical, phylogenetic processes and external ecological pressures. Heterochronic weighting proves to be an effective method to quantify heterochronic trends within a phylogenetic framework and is readily applicable to any group of organisms that have well-defined morphological characteristics, ontogenetic information, and resolved internal relationships.

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.

Site-selectivity analysis of drilling predation traces may provide useful behavioral information concerning a predator interacting with its prey. However, traditional approaches exclude some spatial information (i.e., oversimplified trace position) and are dependent on the scale of analysis (e.g., arbitrary grid system used to divide the prey skeleton into sectors). Here we introduce the spatial point pattern analysis of traces (SPPAT), an approach for visualizing and quantifying the distribution of traces on shelled invertebrate prey, which includes improved collection of spatial information inherent to drillhole location (morphometric-based estimation), improved visualization of spatial trends (kernel density and hotspot mapping), and distance-based statistics for hypothesis testing (K-, L-, and pair correlation functions). We illustrate the SPPAT approach through case studies of fossil samples, modern beach-collected samples, and laboratory feeding trials of naticid gastropod predation on bivalve prey. Overall results show that kernel density and hotspot maps enable visualization of subtle variations in regions of the shell with higher density of predation traces, which can be combined with the maximum clustering distance metric to generate hypotheses on predatory behavior and anti-predatory responses of prey across time and geographic space. Distance-based statistics also capture the major features in the distribution of traces across the prey skeleton, including aggregated and segregated clusters, likely associated with different combinations of two modes of drilling predation, edge and wall drilling. The SPPAT approach is transferable to other paleoecologic and taphonomic data such as encrustation and bioerosion, allowing for standardized investigation of a wide range of biotic interactions.

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.

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.

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.

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.

Historians have debated whether pathways and events from the past to the present are influenced largely by contingency, the dependence of outcomes on particular prior conditions, or whether there is long-term emergent directional change. Previous arguments for predictability in evolutionary history relied on the high frequency of convergence, but the repeated evolution of widely favored adaptations need not imply long-term directionality. Using evidence from the fossil record and arguments concerning the metabolic evolution of organisms, I show here that power (total energy taken up and expended per unit time) has increased stepwise over time at ecosystem-level and global scales thanks to the ratchet-like, cumulative effects of competition and cooperation and to the disproportionate influence of powerful top competitors and opportunistic species on emergent ecosystem properties and processes. The history of life therefore exhibits emergent directionality at large ecosystem-wide scales toward greater power.

Bone healing is an important survival mechanism, allowing vertebrates to recover from injury and disease. Here we describe newly recognized paleopathologies in the hindlimbs of the early tetrapods Crassigyrinus scoticus and Eoherpeton watsoni from the early Carboniferous of Cowdenbeath, Scotland. These pathologies are among the oldest known instances of bone healing in tetrapod limb bones in the fossil record (about 325 Ma). X-ray microtomographic imaging of the internal bone structure of these lesions shows that they are characterized by a mass of trabecular bone separated from the shaft's trabeculae by a layer of cortical bone. We frame these paleopathologies in an evolutionary context, including additional data on bone healing and its pathways across extinct and extant sarcopterygians. These data allowed us to synthesize information on cell-mediated repair of bone and other mineralized tissues in all vertebrates, to reconstruct the evolutionary history of skeletal tissue repair mechanisms. We conclude that bone healing is ancestral for sarcopterygians. Furthermore, other mineralized tissues (aspidin and dentine) were also capable of healing and remodeling early in vertebrate evolution, suggesting that these repair mechanisms are synapomorphies of vertebrate mineralized tissues. The evidence for remodeling and healing in all of these tissues appears concurrently, so in addition to healing, these early vertebrates had the capacity to restore structure and strength by remodeling their skeletons. Healing appears to be an inherent property of these mineralized tissues, and its linkage to their remodeling capacity has previously been underappreciated.

Preferential dissolution of the biogenic carbonate polymorph aragonite promotes preservational bias in shelly marine faunas. While field studies have documented the impact of preferential aragonite dissolution on fossil molluscan diversity, its impact on regional and global biodiversity metrics is debated. Epicontinental seas are especially prone to conditions that both promote and inhibit preferential dissolution, which may result in spatially extensive zones with variable preservation. Here we present a multifaceted evaluation of aragonite dissolution within the Late Cretaceous Western Interior Seaway of North America. Occurrence data of mollusks from two time intervals (Cenomanian/Turonian boundary, early Campanian) are plotted on new high-resolution paleogeographies to assess aragonite preservation within the seaway. Fossil occurrences, diversity estimates, and sampling probabilities for calcitic and aragonitic fauna were compared in zones defined by depth and distance from the seaway margins. Apparent range sizes, which could be influenced by differential preservation potential of aragonite between separate localities, were also compared. Our results are consistent with exacerbated aragonite dissolution within specific depth zones for both time slices, with aragonitic bivalves additionally showing a statistically significant decrease in range size compared with calcitic fauna within carbonate-dominated Cenomanian–Turonian strata. However, we are unable to conclusively show that aragonite dissolution impacted diversity estimates. Therefore, while aragonite dissolution is likely to have affected the preservation of fauna in specific localities, time averaging and instantaneous preservation events preserve regional biodiversity. Our results suggest that the spatial expression of taphonomic biases should be an important consideration for paleontologists working on paleobiogeographic problems.

The estimation of origination and extinction rates and their temporal variation is central to understanding diversity patterns and the evolutionary history of clades. The fossil record provides the only direct evidence of extinction and biodiversity changes through time and has long been used to infer the dynamics of diversity changes in deep time. The software PyRate implements a Bayesian framework to analyze fossil occurrence data to estimate the rates of preservation, origination, and extinction while incorporating several sources of uncertainty. Building upon this framework, we present a suite of methodological advances including more complex and realistic models of preservation and the first likelihood-based test to compare the fit across different models. Further, we develop a new reversible jump Markov chain Monte Carlo algorithm to estimate origination and extinction rates and their temporal variation, which provides more reliable results and includes an explicit estimation of the number and temporal placement of statistically significant rate changes. Finally, we implement a new C++ library that speeds up the analyses by orders of magnitude, therefore facilitating the application of the PyRate methods to large data sets. We demonstrate the new functionalities through extensive simulations and with the analysis of a large data set of Cenozoic marine mammals. We compare our analytical framework against two widely used alternative methods to infer origination and extinction rates, revealing that PyRate decisively outperforms them across a range of simulated data sets. Our analyses indicate that explicit statistical model testing, which is often neglected in fossil-based macroevolutionary analyses, is crucial to obtain accurate and robust results.

It was once thought that the endemic carnivorous mammals of South America, the metatherian sparassodonts, were driven extinct by North American carnivorans through competitive exclusion. However, sparassodonts went extinct before most groups of carnivorans entered South America; only the endemic Cyonasua-group procyonids (Cyonasua and Chapalmalania), which immigrated to South America nearly 4 million years earlier than other carnivorans, significantly overlapped with sparassodonts in time. In this study, we examine the functional morphology of the dentition of Cyonasua and Chapalmalania through quantitative analysis to determine the dietary habits of these taxa and the degree to which they may have ecologically overlapped sparassodonts and large predatory Neogene didelphimorphians. We find Cyonasua and Chapalmalania to be more carnivorous than extant procyonids, other than Bassariscus, in agreement with previous studies, but more omnivorous than most other carnivorans and all meat-eating South American metatherians, including sparassodonts. The extreme ecological dissimilarity between Cyonasua-group procyonids and members of the endemic South American predator guild may explain why procyonids were able to successfully establish themselves in South America several million years earlier than most other northern mammals (including all other carnivorans): they moved into a previously unoccupied ecological niche (large omnivore) and avoided direct competition with incumbent native species, a situation similar to that documented in historical cases of biological invasion. The omnivorous diets and climbing/swimming abilities of procyonids may have increased their chances for a successful over-water dispersal relative to other carnivorans, further favoring their successful establishment in South America.

Previous studies of cranial shape have established a consistent interspecific allometric pattern relating the relative lengths of the face and braincase regions of the skull within multiple families of mammals. In this interspecific allometry, the facial region of the skull is proportionally longer than the braincase in larger species. The regularity and broad taxonomic occurrence of this allometric pattern suggests that it may have an origin near the base of crown Mammalia, or even deeper in the synapsid or amniote forerunners of mammals. To investigate the possible origins of this allometric pattern, we used geometric morphometric techniques to analyze cranial shape in 194 species of nonmammalian synapsids, which constitute a set of successive outgroups to Mammalia. We recovered a much greater diversity of allometric patterns within nonmammalian synapsids than has been observed in mammals, including several instances similar to the mammalian pattern. However, we found no evidence of the mammalian pattern within Therocephalia and nonmammalian Cynodontia, the synapsids most closely related to mammals. This suggests that the mammalian allometric pattern arose somewhere within Mammaliaformes, rather than within nonmammalian synapsids. Further investigation using an ontogenetic series of the anomodont Diictodon feliceps shows that the pattern of interspecific allometry within anomodonts parallels the ontogenetic trajectory of Diictodon. This indicates that in at least some synapsids, allometric patterns associated with ontogeny may provide a “path of least resistance” for interspecific variation, a mechanism that we suggest produces the interspecific allometric pattern observed in mammals.

Organismal metabolic rates reflect the interaction of environmental and physiological factors. Thus, calcifying organisms that record growth history can provide insight into both the ancient environments in which they lived and their own physiology and life history. However, interpreting them requires understanding which environmental factors have the greatest influence on growth rate and the extent to which evolutionary history constrains growth rates across lineages. We integrated satellite measurements of sea-surface temperature and chlorophyll-a concentration with a database of growth coefficients, body sizes, and life spans for 692 populations of living marine bivalves in 195 species, set within the context of a new maximum-likelihood phylogeny of bivalves. We find that environmental predictors overall explain only a small proportion of variation in growth coefficient across all species; temperature is a better predictor of growth coefficient than food supply, and growth coefficient is somewhat more variable at higher summer temperatures. Growth coefficients exhibit moderate phylogenetic signal, and taxonomic membership is a stronger predictor of growth coefficient than any environmental predictor, but phylogenetic inertia cannot fully explain the disjunction between our findings and the extensive body of work demonstrating strong environmental control on growth rates within taxa. Accounting for evolutionary history is critical when considering shells as historical archives. The weak relationship between variation in food supply and variation in growth coefficient in our data set is inconsistent with the hypothesis that the increase in mean body size through the Phanerozoic was driven by increasing productivity enabling faster growth rates.

The hadrosaurids were a successful group of herbivorous dinosaurs. During the Late Cretaceous, 100 to 66 million years ago, hadrosaurids had high diversity, rapid speciation rates, and wide geographic distribution. Most hadrosaurids were large bodied and had similar postcranial skeletons. However, they show important innovations in the skull, including disparate crests that functioned as socio-sexual display structures, and a complex feeding apparatus, with specialized jaws bearing dental batteries. Little is known about the macroevolutionary processes that produced these evolutionary novelties. Here we provide novel perspectives using evolutionary rate and disparity analyses. Our results show that hadrosaurid cranial evolution was complex and dynamic, but their postcranial skeleton and body size were conservative. High cranial disparity was achieved through multiple bursts of phenotypic innovation. We highlight contrasting evolutionary trends within hadrosaurids between the disparate facial skeleton and crests, which both showed multiple high-rate shifts, and the feeding apparatus, which had low variance and high rates on a single phylogenetic branch leading to the diverse Saurolophidae. We reveal that rapid evolutionary rates were important for producing the high disparity of exaggerated crests and present novel evidence that the hadrosaurid diversification was linked to both a key adaptive innovation in the feeding apparatus and multiple bursts of innovation in socio-sexual displays.

Reduction of body size is a common response of organisms to environmental stress. Studying the early Toarcian succession in the Lusitanian Basin of Portugal, we tested whether the shell size of benthic marine communities of bivalves and brachiopods changed at and before the global, warming–related Toarcian oceanic anoxic event (T-OAE). Statistical analyses of shell size over time show that the mean shell size of communities decreased significantly before the T-OAE. This trend is distinct in brachiopods and is caused by larger-sized species becoming less abundant over time, whereas it is not significant in bivalves, suggesting a decoupled response to environmental stress. Reductions in shell size precede the decline in standardized sample-level species richness associated with the early Toarcian extinction event. Such decreases in the shell size of marine invertebrates, well before the onset of biodiversity change, suggest that reductions in body size more generally may be a precursor of a subsequent loss of species and turnover at the community level caused by climate change. Sedimentological evidence is against hypoxia as a driver of extinction and the preceding size decrease in the brachiopod fauna in the studied succession, although low oxygen levels are widely held responsible for elevated early Toarcian extinction rates globally. Reduction of mean shell size in brachiopods but stasis in bivalves is difficult to explain with ocean acidification, because experimental work shows that brachiopods can be resilient to lowered pH, albeit long-term metabolic costs and potential evolutionary adaptations are unknown. Rising early Toarcian temperatures in the Lusitanian Basin seem to be a plausible factor in both diversity decline associated with the T-OAE and the preceding reductions in mean shell size, because thermal tolerances in modern bivalves are among the highest within marine invertebrates.

Recent excavations of Ediacaran assemblages have revealed striking bed-to-bed variation in diversity–abundance structure, offering potential insight into the ecology and taphonomy of these poorly understood early Metazoan ecosystems. Here we compare faunal variability in Ediacaran assemblages to that of younger benthic assemblages, both fossil and modern. We decompose the diversity of local assemblages into within-collection (α) and among-collection (β) components and show that β diversity in Ediacaran assemblages is unusually high relative to younger assemblages. Average between-bed ecological dissimilarities in the Phanerozoic fossil record are comparable to within-habitat dissimilarities typically observed over meter to kilometer scales in modern benthic marine habitats, but dissimilarities in Ediacaran assemblages are comparable to those typically observed over 10–100 km scales in modern habitats. We suggest that the unusually variable diversity–abundance structure of Ediacaran assemblages is due both to their preservation as near snapshots of benthic communities and to original ecological differences, in particular the paucity of motile taxa and the near lack of predation and infaunalization.

The Ordovician was a time of drastic biological and geological change. Previous work has suggested that there was a dramatic increase in global diversity during this time, but also has indicated that regional dynamics and dynamics in specific environments might have been different. Here, we contrast two paleocontinents that have different geological histories through the Ordovician, namely Laurentia and Baltica. The first was situated close to the equator throughout the whole Ordovician, while the latter has traversed tens of latitudes during the same time. We predict that Baltica, which was under long-term environmental change, would show greater average and interval-to-interval origination and extinction rates than Laurentia. In addition, we are interested in the role of the environment in which taxa originated, specifically, the patterns of onshore–offshore dynamics of diversification, where onshore and offshore areas represent high-energy and low-energy environments, respectively. Here, we predict that high-energy environments might be more conducive for originations.

Our new analyses show that the global Ordovician spike in genus richness from the Dapingian to the Darriwilian Stage resulted from a very high origination rate at the Dapingian/Darriwilian boundary, while the extinction rate remained low. We found substantial interval-to-interval variation in the origination and extinction rates in Baltica and Laurentia, but the probabilities of origination and extinction are somewhat higher in Baltica than Laurentia. Onshore and offshore areas have largely indistinguishable origination and extinction rates, in contradiction to our predictions. The global spike in origination rates at the Dapingian/Darriwilian boundary is apparent in Baltica, Laurentia, and onshore and offshore areas, and abundant variability in diversification rates is apparent over other time intervals for these paleocontinents and paleoenvironments. This observation hints at global mechanisms for the spike in origination rates at the Dapingian/Darriwilian boundary but a domination of more regional and local mechanisms over other time intervals in the Ordovician.