Sign in
PaleontologybioRxiv
Preprint

Global plate model choice impacts reconstructions of the latitudinal biodiversity gradient

Lewis Alan JonesiD, William GeartyiD, Lucas BuffaniD, Bethany J. AlleniD

Abstract

The latitudinal gradient of declining species richness from the Equator towards the poles is one of the most pervasive macroecological patterns on Earth today. However, the ubiquity of this trend over geological timescales remains unclear. One reason for this uncertainty is that palaeobiologists need Global Plate Models (GPMs) to estimate the latitudinal position of organisms remains at time of deposition. However, as GPMs constitute hypotheses for how tectonic plates have moved over Earths history, reconstructions of the latitudinal biodiversity gradient (LBG) might also vary based on the GPM used. Here, using the fossil record of five major marine invertebrate groups, we evaluate the impact of GPM choice on reconstructions of the LBG over the Phanerozoic. Our results show that GPM choice can lead to different conclusions about the shape and strength of LBGs in deep time, even at a coarse spatial scale. These findings suggest additional caution is needed when reconstructing deep-time biogeographic patterns and macroevolutionary events, such as the origin of the present-day LBG. We therefore advocate for future palaeobiogeographic studies to conduct sensitivity analyses investigating the impact of GPM choice on their conclusions, and for greater interdisciplinary collaboration between palaeobiologists and palaeogeographic modellers to avoid common issues in the use of GPMs.

02 /Specification & Classification

Click any tag to browse matching papers

In-depth summary

✦ Multi-agent verified full-text synthesis

Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 1 (14 Jan 2025) · CC BY 4.0

Across 344,258 Phanerozoic marine invertebrate fossil occurrences, four global plate models often disagree on where latitudinal diversity peaks, with all-model agreement in only ~19% of Palaeozoic stages.

Reconstructing deep-time latitudinal biodiversity gradients (LBGs) requires plate models to place fossils at the latitude where they were deposited, yet few palaeobiological studies have tested how the choice of model changes the inferred gradients. This work measures that sensitivity across the Phanerozoic, which bears on claims that gradient shape tracks climate and on dating the origin of the modern LBG.

Key results

  1. ShownThe four Global Plate Models (GPMs) often fail to agree on which palaeolatitudinal bin holds peak genus richness. Agreement is lowest in the Palaeozoic and slightly higher in the Northern than the Southern Hemisphere.

    All four models agreed on the peak bin in ~19% of Palaeozoic, ~47% of Mesozoic and ~26% of Cenozoic stages. Across the Phanerozoic, agreement was ~49% in the Northern vs ~43% in the Southern Hemisphere (coverage-based estimates, quorum 0.4). · Quantifying the LBG; Palaeogeographic reconstructions

  2. ShownAll four models recover the same long-term shift of peak diversity, from the Southern Hemisphere before the Permian to the Northern Hemisphere after it. Even so, model pairs often rank latitudinal bins by richness differently, and they differ least in the Mesozoic.

    Normalised mean pairwise rank-order difference was ~0.25 in the Mesozoic (no pair above ~0.78), ~0.45 in the Palaeozoic and ~0.27 in the Cenozoic; in the latter two, some pairs disagreed completely (score 1). Raw-richness gradients showed this less often than coverage-based estimates. · Quantifying the LBG; Palaeogeographic reconstructions

  3. ShownIn a subset of stage-hemisphere combinations, the choice of model changes the type of gradient: one GPM gives a tropical peak and another a polar or temperate peak. These conflicts are more frequent in the Southern Hemisphere.

    Tropical vs polar peak conflicts occurred in 25 stage-hemisphere combinations (7 Northern, 18 Southern; e.g. Southern Hemisphere Givetian). Tropical vs temperate peak conflicts occurred in 68 combinations (33 Northern, 35 Southern). · Discussion

  4. ShownDifferences in palaeolatitude between GPMs generally grow with the age of the reconstruction, but not steadily: the Permian has low differences, and the spread is wider before the Late Jurassic.

    Average median pairwise palaeolatitudinal difference per stage was 2.51º in the Cenozoic, 4.49º in the Mesozoic and 10.34º in the Palaeozoic. Some collections differ by more than 25º latitude, e.g. Cambrian (n = 847), Ordovician (n = 680) and Jurassic (n = 488). · Quantifying the LBG; Palaeogeographic reconstructions; Conclusion

  5. ShownEach GPM can place a different set of localities, because of the plate ages it assumes, and this changes bin sample sizes. The authors attribute gradient differences in young intervals mainly to this, not to shifts in reconstructed position.

    Lutetian: of 601 localities, 597 (PALEOMAP), 592 (GOLONKA), 546 (MERDITH2021) and 496 (TorsvikCocks2017) were reconstructed. Anisian: 553 of 879 (~63%) with TorsvikCocks2017. None of 52 Tortonian Panama localities could be placed with GOLONKA or MERDITH2021. Usable bins differed by one bin in 47 and two in 17 of 97 stages. · Palaeogeographic reconstructions; Discussion

Methods, data and reproducibility

Paleobiology Database occurrences of five marine invertebrate clades were binned by geological stage, given palaeocoordinates with four GPMs via the GPlates Web Service, and assigned to six equal-area latitudinal bins. Genus richness was estimated by coverage-based rarefaction and extrapolation (coverage 0.4), and models were compared on peak-bin agreement and normalised rank-order difference.

  • PBDB occurrences of Bivalvia, Brachiopoda, Cephalopoda, Gastropoda and Trilobita, Fortunian to Holocene (541–0 Ma), marine settings and regular preservation only · marine invertebrates (five fossil clades) · n = 344,258 occurrences; 81,180 (pseudo-)collections
  • Unique spatiotemporal fossil localities across 97 geological stages, each reconstructed with four GPMs · marine invertebrates (five fossil clades) · n = 47180 localities
  • Data and code archive for this study · 10.5281/zenodo.14616817
Methods & tools
Paleobiology Database (PBDB) API, GPlates Web Service, palaeoverse R package ver. 1.4.0 (bin_time, palaeorotate, lat_bins_area), iNEXT R package ver. 3.0.1 (coverage-based Hill numbers), deeptime R package ver. 2.1.0, Tidyverse, GOLONKA plate model, PALEOMAP plate model, MERDITH2021 plate model, TorsvikCocks2017 plate model, International Geological Time Scale ver. 2023/09, Shareholder Quorum Subsampling (SQS)-equivalent rarefaction
Software & environment
R with palaeoverse ver. 1.4.0, iNEXT ver. 3.0.1, deeptime ver. 2.1.0 and Tidyverse; International Geological Time Scale ver. 2023/09; palaeorotation via the GPlates Web Service.
Reproducibility gaps
Palaeocoordinates come from the external GPlates Web Service and occurrences from a dated PBDB API download, so exact reruns depend on those services. Processed data and code are deposited on GitHub and Zenodo.
Limitations (6)
  • Coverage-based extrapolation may be overly sensitive: richness estimates can shift when a few collections move between bins. Raw-richness gradients were more consistent across GPMs, which may partly reflect sampling bias instead. · stated by the authors
  • Only six equal-area latitudinal bins were used. The authors state that finer latitudinal resolution yields noisier, patchier gradients driven more by sampling bias than by true diversity, and that subdividing into grid cells would add similar noise. · stated by the authors
  • The normalised rank-order difference depends on how many bins are sampled: with fewer bins, fewer rank changes yield a higher score. Alternative metrics (sum of squares, regression fits, null gradients) were not assessed. · stated by the authors
  • PBDB coverage is patchy, with empty palaeolatitudinal bins and a few entirely empty stages. Sparser Southern Hemisphere sampling likely makes that hemisphere more sensitive to small differences in the data. · stated by the authors
  • The study measures how much the models differ but uses no independent palaeomagnetic or geological benchmark, so it cannot say which of the four GPMs gives more accurate palaeolatitudes. · noted in this summary
  • The gradients appear to pool all five clades at a single coverage level (0.4). Clade-specific patterns and how sensitive the conclusions are to the chosen quorum level are not reported. · noted in this summary

Who should read:Read this if you reconstruct deep-time biogeographic or diversity patterns from fossil occurrences and need evidence on how plate model choice can change your conclusions.

Reproducibility ledger

Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.

Sign in to file a reproducibility report.