Lithofacies and temporal variation predict composition of Siluro-Devonian vertebrate, invertebrate, and plant communities
lead: Professor Corey Bradshaw
collaborators: Richard Cloutier, Clive Burrett, Mongkol Udchachon
While Siluro-Devonian vertebrate communities have been well-described, co-variation in invertebrate and plant communities from the same period has not been adequately assessed. Using methods from community ecology, we tested how the combined vertebrate, invertebrate, and plant communities from 105 fossiliferous sites around the world varied relative to lithofacies and time. We built logistic linear models to examine variation in the three main communities (vertebrates, invertebrates, plants). We found that the type of lithofacies (6 categories) explained the most variation (% deviance explained [%DE] = 8.0), followed by geological period (4 categories); however, community type had little explanatory power (%DE = 0.05). Variation among lithofacies in Silurian communities was driven predominately by low presence probability in limestone, whereas variation in Devonian communities was driven mainly by higher presence probabilities in sandstone and shale. There was little support for a period×lithofacies interaction, nor was much variation explained by any factor when we divided communities by major paleo-environmental category (marine, alluvial/deltaic, freshwater/estuarine). We then applied a stochastic variant of a permutation analysis of variance to the entire community (10 vertebrate, 20 invertebrate, 7 plant taxa) to examine the relative explanatory power of lithofacies and temporal variation on composition. There was an effect of the lithofacies (R2perm = 0.166; pperm; driven mainly by limestone) and categorical period (R2perm = 0.271; pperm = 0.0179; driven mainly by the Middle Devonian) (but not their interaction: pperm = 0.152) on community composition. Treating time as a continuous variable (age) also demonstrated an effect on community composition (R2perm = 0.093; pperm = 0.0052), with incidence tending to rise in most taxa from 435 to 360 million years ago. Model predictions can now be used to show spatial variation in the relative presence probability for any of the taxa we modelled provided detailed distribution maps of lithofacies are available.
silurodev.R
: main code to generate resultsr.squared.R
: functions to estimate goodness of fit for linear modelsnew_lmer_AIC_tables3.r
: functions to compare linear modelspalaeocoords.ipynb
: Jupyter Notebook with Python code to project current-era lat/lon coordinates to palaeo-coordinates, with corresponding plate tectonics, for anytime in the last 1 billion years (modified code supplied by @amer7632; detailed instructions here)
- sildevcomp2.txt: fossiliferous site database
- plate-model-repo: data for hindcasted tectonic reconstruction and generation of palaeo coordinates (maintain sub-folder hierarchy & unzip files in 'Topologies' sub-folder)
- sitecoords: .csv files with site coordinates — all: sitelocs.csv; and divided by major period: Silurian (SIL), Lochkovian (LOCH), Emsian (EMS), Eifelian (EIF), Givetian (GIV), Frasnian (FRAS), Famennian (FAM)
performance
, sjPlot
, lme4
, ggplot2
, stringr
, vegan
, parallel
, vcdExtra
, plyr
, dplyr
Install GPlately
, a Python interface to accelerate spatio-temporal data analysis leveraging pyGPlates
and PlateTectonicTools
(references: Merdith et al. 2021 Earth-Science Reviews 214:103477; Mather et al. 2024 Geosciences Data Journal 11:3-10)
pyGPlates
, plate-model-manager
, Shapely
, NumPy
. SciPy
, Matplotlib
, Cartopy
, Pooch
, GeoPandas
, netCDF4