Benthic stack

The benthic d18O LR4 stack spans 5.3Ma. It is based on an average of 57 benthic d18O records (which measure global ice volume and deep ocean temperature), globally distributed sites along the Earth. Among with the ice cores (from Antartida and Greenland), the benthic stack provides an independent proxy, as a template for global clima and a stratrigraphic tool for paleoceanographers. For example, using the benthic stack we can identify the glacial cycles, climate cycles from Martrat et al. (2004) or marine isotopic stages (so called MIS).

References:
Lisiecki, L. E., and M. E. Raymo (2005), A Pliocene-Pleistocene stack of 57 globally distributed benthic d18O records,
Paleoceanography, 20, PA1003, doi:10.1029/2004PA001071

Here I show a simple code to download it directly from the internet (don't forget to cite), so we can use it in R. No extra packages needed. The plot shown is just a simple example for the first 500 ka (note also the way to write the legend, where d18O is on a correct syntax, I hope :) which I was unable to do in blogger).



# Benthic stack 

 # database created with 1000 corers around the world. 
 # http://lorraine-lisiecki.com/stack.html

 # Please cite: Lisiecki, L. E., and M. E. Raymo (2005), A Pliocene-
 # Pleistocene stack of 57 globally distributed benthic d18O records, 
 # Paleoceanography,20, PA1003, doi:10.1029/2004PA001071.  
 # readLines("http://lorraine-lisiecki.com/LR04stack.txt")
 benthic.stack=read.table("http://lorraine-lisiecki.com/LR04stack.txt", skip=5, col.names=c("age","d18O","se")); 
 benthic.stack$d18O=benthic.stack$d18O*-1
 plot(d18O~c(age), benthic.stack[benthic.stack$age<450,], typ="l", col=4, xlab="Age (ka)", ylab=expression(paste(delta^18~O," Benthic stack")))
 abline(v=seq(round(min(benthic.stack$age)),round(max(benthic.stack$age)), by=20), col=8, lty=3)

Created by Pretty R at inside-R.org