### Fate of CO2 Emissions ###################################################################### # D Kelly O'Day http://chartsgraphs.wordpress.com Sept 8, 2009 # # Reproduce Hansen & IPCC charts that show time series of mass of # # CO2 emissions & atmosp CO2 in PgC # ################################################################################################ # Set pars par(las=1); par(ps=10) par(oma=c(2,2,1,1)); par(mar=c(2,4,1,1)) par(font.main = 4) # Read Data file link <- "http://spreadsheets.google.com/pub?key=tCVSIFNBYVBySwfqHDtDksg&output=csv" df <- read.table(link, header=T, skip=0,sep=",", na.strings = c("*", "NA")) names(df) <- c("yr", "emis", "change") # Plot setup plot_title <- expression(paste("Fossil Fuel Emissions and Atmospheric ", CO[2], " Changes (PgC/yr)")) ## Create df for polygon df1 <- data.frame(df$yr,df$emis) df2 <- data.frame(df$yr, df$change) df2 <- df2[order(-df$yr),] px<- append(df1$df.yr, df2$df.yr,after=nrow(df1)) py <- append(df1$df.emis, df2$df.change, after=nrow(df2)) p <- data.frame(px, py) # Notes for plot a_note <- "Annual Emissions" b_note <- "Soaked-Up\n Amount" c_note <- "Airborne \nAmount" sum_emis <- sum(df$emis, na.rm = T) sum_change <- sum(df$change, na.rm=T) d_note <- paste("Emissions: ", sum_emis, " PgC") e_note <- paste("CO2 Airborne: ", sum_change, "PgC") f_note <- paste("CO2 Airborne: ", signif(100*sum_change/sum_emis, 3),"% of Emissions") g_note <- " Summary: 1960 - 2006" # Plot plot(df$yr, df$emis,type = "n", xlab="", ylab = "PgC/yr", ylim=c(0,8), xlim=c(1960, 2006), main=plot_title, cex.main=0.85, xaxs="i", yaxs="i") polygon(px,py, border=NA, col="lightgreen") points(df$yr, df$change, type="l", col = "white") points(df$yr, df$emis, type="l", col = "red", lwd=2) text(1980, 6, a_note,cex=0.8) rect(1979.5,4.2, 1987,4.8, col = "white", border=NA) text(1983, 4.5, b_note,cex=0.8) text(1985, 2, c_note,cex=0.8) arrows(1980, 5.85, 1985, 5.5, col = "darkgrey", length = 0.05) text(1962, 7.6, d_note,cex=0.75, adj=0 ) text(1962, 7.4, e_note, cex=0.75, adj = 0) text(1962, 7.2, f_note, cex = 0.75, adj = 0) text(1962, 7.8, g_note, cex=0.85, adj=0) # Chart Annotation my_date <- format(Sys.time(), "%m/%d/%y") mtext("D Kelly O'Day - http://chartsgraphs.wordpress.com", side = 1, line = .5, cex=0.8, outer = T, adj = 0) mtext(my_date, side = 1, line =.5, cex = 0.8, outer = T, adj = 1)## Add script name to each plot # Annual airborne amount calcs ap <- df$change/df$emis max_ap <- max(ap, na.rm=T) mean_ap <- mean(ap, na.rm=T) min_ap <- min(ap, na.rm=T) cat(c("max", max_ap), fill = 30) cat(c("mean", mean_ap), fill = 30) cat(c("min", min_ap),fill = 30)