6.2 Loading and visualizing the \(\delta\)D series

We’ll also start by formatting Age in kyr.

dDdf <- read.csv("https://github.com/LinkedEarth/Pyleoclim_util/raw/master/pyleoclim/data/edc3deuttemp2007.csv")

dDdf <- dDdf %>%
  mutate(AgeKy = Age/1000) %>%
  select(AgeKy, Deuterium, Temperature)

head(dDdf)

##        AgeKy Deuterium Temperature
## 1 0.03837379    -390.9        0.88
## 2 0.04681203    -385.1        1.84
## 3 0.05505624    -377.8        3.04
## 4 0.06441511    -394.1        0.35
## 5 0.07315077    -398.7       -0.42
## 6 0.08193244    -395.9        0.05

Let’s see a simple plot

ggplot(data=dDdf, mapping = aes(x=AgeKy, y=Deuterium)) + 
  geom_line() +
  ylab("\u03B4D (\u2030)") +
  xlab("Age_(ky_BP)") +
  ggtitle("EPICA Dome C Deuterium") +
  theme_light()