Mean annual precipitation explains spatiotemporal patterns of Cenozoic mammal beta diversity and latitudinal diversity gradients in North America

The hierarchy of factors predicting the latitudinal diversity gradient

The numerous explanations for why Earth's biodiversity is concentrated at low latitudes fail to explain variation in the strength and even direction of the gradient through deep time. Consequently, we do not know if today's gradient is representative of what might be expected on other planets or is merely an idiosyncrasy of Earth's history. We propose a hierarchy of factors driving the latitudinal distribution of diversity: (i) over geologically long time spans, diversity is largely predicted by climate; (ii) when climatic gradients are shallow, diversity tracks habitat area; and (iii) historical contingencies linked to niche conservatism have geologically short-term, transient influence at most. Thus, latitudinal diversity gradients, although variable in strength and direction, are largely predictable on our planet and possibly others.

Late quaternary biotic homogenization of North American mammalian faunas

Biotic homogenization—increasing similarity of species composition among ecological communities—has been linked to anthropogenic processes operating over the last century. Fossil evidence, however, suggests that humans have had impacts on ecosystems for millennia. We quantify biotic homogenization of North American mammalian assemblages during the late Pleistocene through Holocene (~30,000 ybp to recent), a timespan encompassing increased evidence of humans on the landscape (~20,000–14,000 ybp). From ~10,000 ybp to recent, assemblages became significantly more homogenous (>100% increase in Jaccard similarity), a pattern that cannot be explained by changes in fossil record sampling. Homogenization was most pronounced among mammals larger than 1 kg and occurred in two phases. The first followed the megafaunal extinction at ~10,000 ybp. The second, more rapid phase began during human population growth and early agricultural intensification (~2,000–1,000 ybp). We show that North American ecosystems were homogenizing for millennia, extending human impacts back ~10,000 years.

Biotic homogenization, which is increased similarity in the composition of species among communities, is rising due to human activities. Using North American mammal fossil records from the past 30,000 years, this study shows that this phenomenon is ancient, beginning between 12,000 and 10,000 years ago with the extinction of the mammal megafauna.

Pronghorn (Antilocapra americana) enamel phosphate δ18O values reflect climate seasonality: Implications for paleoclimate reconstruction

Stable oxygen isotope (δ18O) compositions from vertebrate tooth enamel are widely used as biogeochemical proxies for paleoclimate. However, the utility of enamel oxygen isotope values for environmental reconstruction varies among species. Herein, we evaluate the use of stable oxygen isotope compositions from pronghorn (Antilocapra americana Gray, 1866) enamel for reconstructing paleoclimate seasonality, an elusive but important parameter for understanding past ecosystems. We serially sampled the lower third molars of recent adult pronghorn from Wyoming for δ18O in phosphate (δ18OPO4) and compared patterns to interpolated and measured yearly variation in environmental waters as well as from sagebrush leaves, lakes, and rivers (δ18Ow). As expected, the oxygen isotope compositions of phosphate from pronghorn enamel are enriched in 18O relative to environmental waters. For a more direct comparison, we converted δ18Ow values into expected δ18OPO4* values (δ18OW-PO4*). Pronghorn δ18OPO4 values from tooth enamel record nearly the full amplitude of seasonal variation from Wyoming δ18OW-PO4* values. Furthermore, pronghorn enamel δ18OPO4 values are more similar to modeled δ18OW-PO4* values from plant leaf waters than meteoric waters, suggesting that they obtain much of their water from evaporated plant waters. Collectively, our findings establish that seasonality in source water is reliably reflected in pronghorn enamel, providing the basis for exploring changes in the amplitude of seasonality of ancient climates. As a preliminary test, we sampled historical pronghorn specimens (1720 ± 100 AD), which show a mean decrease (a shift to lower values) of 1-2‰ in δ18OPO4 compared to the modern specimens. They also exhibit an increase in the δ18O amplitude, representing an increase in seasonality. We suggest that the cooler mean annual and summer temperatures typical of the 18th century, as well as enhanced periods of drought, drove differences among the modern and historical pronghorn, further establishing pronghorn enamel as excellent sources of paleoclimate proxy data.

The preservation potential of terrestrial biogeographic patterns

Extinction events in the geological past are similar to the present-day biodiversity crisis in that they have a pronounced biogeography, producing dramatic changes in the spatial distributions of species. Reconstructing palaeobiogeographic patterns from fossils therefore allows us to examine the long-term processes governing the formation of regional biotas, and potentially helps build spatially explicit models for future biodiversity loss. However, the extent to which biogeographic patterns can be preserved in the fossil record is not well understood. Here, we perform a suite of simulations based on the present-day distribution of North American mammals, aimed at quantifying the preservation potential of beta diversity and spatial richness patterns over extinction events of varying intensities, and after applying a stepped series of taphonomic filters. We show that taphonomic biases related to body size are the biggest barrier to reconstructing biogeographic patterns over extinction events, but that these may be compensated for by both the small mammal record preserved in bird castings, as well as range expansion in surviving species. Overall, our results suggest that the preservation potential of biogeographic patterns is surprisingly high, and thus that the fossil record represents an invaluable dataset recording the changing spatial distribution of biota over key intervals in Earth History.

Spatial sampling heterogeneity limits the detectability of deep time latitudinal biodiversity gradients

The latitudinal biodiversity gradient (LBG), in which species richness decreases from tropical to polar regions, is a pervasive pattern of the modern biosphere. Although the distribution of fossil occurrences suggests this pattern has varied through deep time, the recognition of palaeobiogeographic patterns is hampered by geological and anthropogenic biases. In particular, spatial sampling heterogeneity has the capacity to impact upon the reconstruction of deep time LBGs. Here we use a simulation framework to test the detectability of three different types of LBG (flat, unimodal and bimodal) over the last 300 Myr. We show that heterogeneity in spatial sampling significantly impacts upon the detectability of genuine LBGs, with known biodiversity patterns regularly obscured after applying the spatial sampling window of fossil collections. Sampling-standardization aids the reconstruction of relative biodiversity gradients, but cannot account for artefactual absences introduced by geological and anthropogenic biases. Therefore, we argue that some previous studies might have failed to recover the ‘true’ LBG type owing to incomplete and heterogeneous sampling, particularly between 200 and 20 Ma. Furthermore, these issues also have the potential to bias global estimates of past biodiversity, as well as inhibit the recognition of extinction and radiation events.

Mammal Community Structure through the Paleocene-Eocene Thermal Maximum

AbstractHuman-mediated species invasion and climate change are leading to global extinctions and are predicted to result in the loss of important axes of phylogenetic and functional diversity. However, the long-term robustness of modern communities to invasion is unknown, given the limited timescales over which they can be studied. Using the fossil record of the Paleocene-Eocene thermal maximum (PETM; ∼56 Ma) in North America, we evaluate mammalian community-level response to a rapid global warming event (5°-8°C) and invasion by three Eurasian mammalian orders and by species undergoing northward range shifts. We assembled a database of 144 species body sizes and created a time-scaled composite phylogeny. We calculated the phylogenetic and functional diversity of all communities before, during, and after the PETM. Despite increases in the phylogenetic diversity of the regional species pool, phylogenetic diversity of mammalian communities remained relatively unchanged, a pattern that is invariant to the tree dating method, uncertainty in tree topology, and resolution. Similarly, body size dispersion and the degree of spatial taxonomic turnover of communities remained similar across the PETM. We suggest that invasion by new taxa had little impact on Paleocene-Eocene mammal communities because niches were not saturated. Our findings are consistent with the numerous studies of modern communities that record little change in community-scale richness despite turnover in taxonomic composition during invasion. What remains unknown is whether long-term robustness to biotic and abiotic perturbation are retained by modern communities given global anthropogenic landscape modification.

The latitudinal diversity gradient of tetrapods across the Permo-Triassic mass extinction and recovery interval

The decline in species richness from the equator to the poles is referred to as the latitudinal diversity gradient (LDG). Higher equatorial diversity has been recognized for over 200 years, but the consistency of this pattern in deep time remains uncertain. Examination of spatial biodiversity patterns in the past across different global climate regimes and continental configurations can reveal how LDGs have varied over Earth history and potentially differentiate between suggested causal mechanisms. The Late Permian-Middle Triassic represents an ideal time interval for study, because it is characterized by large-scale volcanic episodes, extreme greenhouse temperatures and the most severe mass extinction event in Earth history. We examined terrestrial and marine tetrapod spatial biodiversity patterns using a database of global tetrapod occurrences. Terrestrial tetrapods exhibit a bimodal richness distribution throughout the Late Permian-Middle Triassic, with peaks in the northern low latitudes and southern mid-latitudes around 20-40° N and 60° S, respectively. Marine reptile fossils are known almost exclusively from the Northern Hemisphere in the Early and Middle Triassic, with highest diversity around 20° N. Reconstructed terrestrial LDGs contrast strongly with the generally unimodal gradients of today, potentially reflecting high global temperatures and prevailing Pangaean super-monsoonal climate system during the Permo-Triassic.

Structure-function covariation with nonfeeding ecological variables influences evolution of feeding specialization in Carnivora

Skull shape convergence is pervasive among vertebrates. Although this is frequently inferred to indicate similar functional underpinnings, neither the specific structure-function linkages nor the selective environments in which the supposed functional adaptations arose are commonly identified and tested. We demonstrate that nonfeeding factors relating to sexual maturity and precipitation-related arboreality also can generate structure-function relationships in the skulls of carnivorans (dogs, cats, seals, and relatives) through covariation with masticatory performance. We estimated measures of masticatory performance related to ecological variables that covary with cranial shape in the mammalian order Carnivora, integrating geometric morphometrics and finite element analyses. Even after accounting for phylogenetic autocorrelation, cranial shapes are significantly correlated to both feeding and nonfeeding ecological variables, and covariation with both variable types generated significant masticatory performance gradients. This suggests that mechanisms of obligate shape covariation with nonfeeding variables can produce performance changes resembling those arising from feeding adaptations in Carnivora.

Late Cenozoic onset of the latitudinal diversity gradient of North American mammals

The decline of species richness from equator to pole, or latitudinal diversity gradient (LDG), is nearly universal among clades of living organisms, yet whether it was such a pervasive pattern in the geologic past remains uncertain. Here, we calculate the strength of the LDG for terrestrial mammals in North America over the past 65 My, using 27,903 fossil occurrences of Cenozoic terrestrial mammals from western North America downloaded from the Paleobiology Database. Accounting for temporal and spatial variation in sampling, the LDG was substantially weaker than it is today for most of the Cenozoic and the robust modern LDG of North American mammals evolved only over the last 4 My. The strength of the LDG correlates negatively with global temperature, suggesting a role of global climate patterns in the establishment and maintenance of the LDG for North American mammals.