(data from Zachos et al 2001, 2008).
data<-read.csv("zachos2008.csv",head=T)
label.O<-expression(paste(delta ,O^18 , sep = ""))
label.C<-expression(paste(delta ,C^13 , sep = ""))
par(tcl=-0.3)
par(mfrow=c(2,1))
par(cex.axis=.8)
par(mgp = c(1.7, .4, 0)) #axis title, axis labels and axis line.default c(3, 1, 0).
par(mar = c(0, 3.5, 1.5, 3)) #c(bottom, left, top, right)
plot(data$Age,data$d18O.5pt,xlim=c(0,65), ylim=c(5,-1),type="l",col="blue",lwd=1.5,
ylab="", axes=F)
axis(4, xlab=label.O,col.axis = "blue",lwd = 2)
mtext(text=label.O, side=4, line=2)
par(mar = c(3, 3.5, 0.5, 3)) #c(bottom, left, top, right)
plot(data$Age,data$d13C.5pt,xlim=c(0,65), ylim=c(-1,2.5),type="l",col="green4",lwd=1.5,
ylab=label.C,xlab="Age (Ma)", axes=F)
axis(1,lwd = 2)
axis(2,col.axis = "green4",lwd = 2)
Sea level reconstructions
(data from Haq et al 1987, Miller et al 2005, Muller et al 2008, De Boer et al 2010)
haq<-read.csv("haq_sealevel_curve.csv",head=T)
plot(haq$age,haq$seaLevel,type="l", xlim=c(0,65),ylim=c(-100,220),
xlab="",ylab="Sea Level (m)", col="darkviolet")
abline(h=0,lty=2,col="grey")
miller<-read.csv("MillerCurve2005.csv",head=T)
points(miller[,2],miller[,1],type="l", xlim=c(0,65),ylim=c(-100,200),
ylab="",col="red")
muller<-read.csv("MullerCurve2008.csv",head=T)
points(muller[,1],muller[,2],type="l", xlim=c(0,65),xlim=c(-100,400),
xlab="",col="black")
points(muller[,1],muller[,2]-muller[,3],type="l", xlim=c(0,65),
ylab="",col="grey")
points(muller[,1],muller[,2]+muller[,4],type="l", xlim=c(0,65),
ylab="",col="grey")
DeBoer<-read.csv("DeBoer_data.csv",head=T)
points(DeBoer[,1]*-0.001,DeBoer[,2],type="l",ylab="",col="green3") #notice *-0.001 to standarize scale
text(31,-18,"Miller",pch=.3,col="red")
text(32,208,"Haq",pch=.3,col="darkviolet")
text(2,101,"Müller",pch=.3)
text(7,-82,"DeBoer",pch=.3,col="green3")
(data from Beerling and Royer 2011).
library(zoo)
Beerling<-read.csv("BeerlingTempCO2.csv",head=T)
plot(Beerling$age,Beerling$CO2,type="p",xlim=c(0,65), ylim=c(0,3500),
ylab="Atmospheric CO2 (ppm)", col="black",pch=19, xlab="Age (Ma)")
arrows(Beerling$age, Beerling$CO2Down,Beerling$age, Beerling$CO2Up,
code=3,angle=90,length = 0.05, col="grey")
points(Beerling$age,Beerling$CO2,pch=20)
ma<-rollmean(Beerling$CO2,3,fill = list(NA, NULL, NA))
lines(Beerling$age,ma,col="orange",lwd=2)
abline(h=390, lty=2)
