A Bayesian, global planktic foraminifera core top calibration to sea-surface temperature (SST), for R.
NOTE This repository and package is no longer actively maintained.
bayfoxr is a suite of linear Bayesian calibration models for planktic core top foraminiferal δ18O (δ18Oc) and SSTs. These calibrations are especially useful because they capture the uncertainty in the relationship between modern SSTs and core top δ18Oc. This package is a companion to a paper currently under preparation for the journal "Paleoceanography and Paleoclimatology".
library(bayfoxr)
data(bassriver)The bassriver is example data that comes with the package. It is marine core samples from John et al. (2008). The data.frame has two columns: "depth", giving down-core depth in meters, and "d18o", foraminifera (Morozovella spp.) calcite d18O samples (‰ VPDB). The core samples cover the Paleocene-Eocene thermal maximum (PETM).
Let's run this data through our annual pooled calibration model to make inferences about past SST. Morozovella spp. is a nonexant species so, we're using modern planktic foraminifera as an analog with this pooled calibration.
sst <- predict_seatemp(bassriver$d18o, d18osw = 0.0,
prior_mean = 30.0, prior_std = 20.0)The predict function then spits out a prediction object. Note that we need to specify d18O for seawater (d18osw), and a prior mean and standard deviation for our SST inference. See help(predict_seatemp) for more details, or help(predict_d18oc) for the reversed, "forward" model.
The sst variable contains an ensemble rather than single prediction points because the calibration is a Bayesian regression model. This ensemble is in sst[['ensemble']]. Here we get median and 90% interval for the prediction:
quantile(sst, probs = c(0.05, 0.50, 0.95))We can also make a quick and dirty plot to visualize the inference:
predictplot(x = bassriver$depth, y = sst, ylim = c(20, 40),
ylab = "SST (°C)", xlab = "Depth (m)")Please cite our work if you use bayfoxr in your research. We have a paper currently in preparation and I'll be sure to update this section with the citation as soon as the paper is out.
To cite the code repository directly use:
Malevich, Steven B., 2019. bayfoxr. <https://github.com/brews/bayfoxr >.
Alternatively, you can cite the package in R's CRAN repository. You can see this information by running citation("bayfoxr") in an R session.
bayfoxr is hosted on the offical CRAN repository. Install it from an R session with:
install.packages("bayfoxr")Bleeding edge and development versions of the package can be installed with devtools. Assuming you have devtools installed in R, you can install bayfoxr with:
devtools::install_github("brews/bayfoxr")Documentation is included in the code and can be viewed in R. Please file issues and requests in the bug tracker.
Want to contribute? We're following Hadley's packaging workflow and code style. Please fork away and get in touch if you have a feature or bug fix.