Skip to content

R code for building stochastic demographic models of extinct megafauna of Sahul

License

Notifications You must be signed in to change notification settings

CABAH/MegafaunaSusceptibility

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

50 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

Megafauna Susceptibility

thylacine

R code accompanies article:

BRADSHAW, CJA, CN JOHNSON, J LLEWELYN, V WEISBECKER, G STRONA, F SALTRÉ. 2021. Relative demographic susceptibility does not explain the extinction chronology of Sahul’s megafauna. eLife 10: e63870. doi:10.7554/eLife.63870

AIM: construct plausible stochastic demographic models for main Sahul megafauna to determine relative demographic susceptibility to environmental change & novel predation (human) sources

Abstract

The causes of Sahul’s megafauna extinctions remain uncertain, although several interacting factors were likely responsible. To examine the relative support for hypotheses regarding plausible ecological mechanisms underlying these extinctions, we constructed the first stochastic, age-structured models for 13 extinct megafauna species from five functional/taxonomic groups, as well as 8 extant species within these groups for comparison. Perturbing specific demographic rates individually, we tested which species were more demographically susceptible to extinction, and then compared these relative sensitivities to the fossil-derived extinction chronology. Our models show that the macropodiformes were the least demographically susceptible to extinction, followed by carnivores, monotremes, vombatiform herbivores, and large birds. Five of the eight extant species were as or more susceptible than the extinct species. There was no clear relationship between extinction susceptibility and the extinction chronology for any perturbation scenario, while body mass and generation length explained much of the variation in relative risk. Our results reveal that the actual mechanisms leading to the observed extinction chronology were unlikely related to variation in demographic susceptibility per se, but were possibly driven instead by finer-scale variation in climate change and/or human prey choice and relative hunting success. Diprotodon


Prof Corey J. A. Bradshaw
Global Ecology, Flinders University, Adelaide, Australia
August 2020
e-mail

Groups/genera modelled

Demographic parameters for each species summarised below (Appendix 2, table 1 in the original article)

Appendix 2 table 1

Repository includes the following files

Simosthenurus

- Sahul megafauna demographic susceptibility-base models.R — constructs bases models for all perturbation scenarios (must be run first) - matrixOperators.R — functions to manipulate matrix models - megsuscept.SCENARIO2.juvsurv.R — runs Scenario 2 (reduction in juvenile survival) - megsuscept.SCENARIO3.fertred.R — runs Scenario 3 (reduction in fertility) - megsuscept.SCENARIO4.survred.R — runs Scenario 4 (reduction in all-ages survival) - megsuscept.SCENARIO5.indrem.R — runs Scenario 5 (increasing individual offtake) - megsuscept.SCENARIO6.catinc.R — runs Scenario 6 (increasing frequency of catastrophic die-offs) - megsuscept.SCENARIO7.catMinc.R — runs Scenario 7 (increasing magnitude of catastrophic die-offs)

** NOTE: For Scenario 7, Scenario 6 must be run first to create input .csv file 'catincQpr.csv' **

Sensitivity analysis

  • GRIWM.jk.sensitivity.R — is the Gaussian-Resampled, Inverse-Weighted McInerny (GRIWM) Signor-Lipps algorithm, including a jack-knife estimator to test senstivity of dates to different assumptions. This R code is applied to the various chronologies for the following taxa (.csv files in the 'chronologies' sub-folder): Diprotodon, Palorchestes, Zygomaturus, Phascolonus, Protemnodon, Sthenurus, Simosthenurus, Procoptodon, Metasthenurus, Genyornis, Thylacoleo, Megalibgwilia

Palorchestes

Thylacoleo

About

R code for building stochastic demographic models of extinct megafauna of Sahul

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • R 100.0%