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 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.

Hadrosaurid dinosaurs, the dominant large-bodied terrestrial herbivores in most Laurasian Late Cretaceous ecosystems, have an exceptional fossil record consisting of many species known from partial ontogenetic series, making them an ideal clade with which to conduct life-history studies. Previous research considered the Dinosaur Park Formation (DPF) of Alberta as an attritional, or time-averaged, sample and interpreted size–frequency distribution of long bones collected from the DPF with three size classes to suggest that hadrosaurids from the DPF attained near-asymptotic body size in under 3 years. This conflicted with previously published osteohistological estimates of 6+ years for penecontemporaneous hadrosaurids from the Two Medicine Formation (TMF) of Montana, suggesting either extreme variation in hadrosaurid growth rates or that size–frequency distributions and/or osteohistology and growth modeling inaccurately estimate ontogenetic age.

We tested the validity of the previously proposed size–age relationship of hadrosaurids from the DPF by significantly increasing sample size and combining data from size–frequency distributions and osteohistology across multiple long-bone elements. The newly constructed size–frequency distributions typically reveal four relatively distinct size–frequency peaks that, when integrated with the osteohistological data, aligned with growth marks. The yearling size class was heavily underrepresented in the size–frequency distribution. If not due to preservation, this suggests that either juvenile (<2 years of age) hadrosaurids from the DPF had increased survivorship following an initially high nestling mortality rate or that yearlings were segregated from adults. A growth-curve analysis revealed asymptotic body size was attained in approximately 7 years, which is consistent with hadrosaurids from the TMF. The data suggest size–frequency distributions of attritional samples underestimate age and overestimate growth rates, but when paired with osteohistology can provide unique life-history insights.

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.

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.