1 Dating Profiles

1.1 Getting Started

Load the following packages

library(ggplot2)
library(dplyr)
library(viridis)
library(ggpubr)

1.2 Load the data

I used data from the following files found in the Github repo FileS3_Dating_CSM-IMP.csv and FileS4_Dating_CSM-REF.csv. Access the files here.

Load your data for impact and reference sites into R

dp_IMP2 <- read.csv("Data/FileS2_Dating_CSM-IMP.csv") 
dp_REF3 <- read.csv("Data/FileS3_Dating_CSM-REF.csv") 

Define columns for core depth/year, activity, sedimentation rate (with error) in both files

#core depth
depth_IMP2 <- dp_IMP2$midpoint
depth_REF3 <- dp_REF3$midpoint

#core year
year_IMP <- dp_IMP2$year
year_REF <- dp_REF3$year

#210-Pb activity
act_IMP2 <- dp_IMP2$activity
act_REF3 <- dp_REF3$activity

# sedimentation rate
sed_IMP2 <- dp_IMP2$sedimetation_rate
sed_REF3 <- dp_REF3$sedimetation_rate

error_IMP <- dp_IMP2$sr_error
error_REF <- dp_REF3$sr_error

1.3 Plot 210-Pb dating profiles

Plot the activity and sedimentation rate

activity <- ggplot() +
  geom_point(data = dp_IMP2, (aes(x= depth_IMP2, y=act_IMP2)), shape=22, fill="paleturquoise4", size=3, color = "paleturquoise4", alpha=0.5) +
  geom_line(data = dp_IMP2, (aes(x = depth_IMP2, y=act_IMP2)), color="paleturquoise4") +
  geom_point(data = dp_REF3, (aes(x = depth_REF3, y=act_REF3)), size =3, color="palegreen3") +
  geom_line(data = dp_REF3, (aes(x = depth_REF3, y=act_REF3)), color="palegreen3") +
  labs(x="Midpoint depth (cm)", y="Pb-210 Activity (Bq/kg)") +
  scale_x_continuous(breaks=seq(0,19, by= 1), limits = c(0,19)) +
  theme_classic()+
  theme(text = element_text(size = 12), panel.grid.major.x=element_line(), axis.text.y = element_text(size=12), axis.title.y = element_text(margin = margin(t=0, r= 20, b=0, l=0)))

activity

sed.rate <- ggplot() +
  geom_point(data = dp_IMP2, (aes(x = depth_IMP2, y=sed_IMP2)),shape=21, color = "paleturquoise4") +
  geom_line(data = dp_IMP2, (aes(x = depth_IMP2, y=sed_IMP2)), color="paleturquoise4") +
  geom_point(data = dp_REF3, (aes(x = depth_REF3, y=sed_REF3)), color="palegreen3") +
  geom_line(data = dp_REF3, (aes(x = depth_REF3, y=sed_REF3)), color="palegreen3") +
  labs(x = "Core Depth (cm/yr)", y = "Sedimentation rate (cm/yr)") +
  scale_x_continuous(breaks=seq(0,19, by= 1), limits = c(0,19)) +
  theme_classic() +
  theme(panel.grid.major.x=element_line(), text = element_text(size = 12), axis.text.y = element_text(size=12), 
        axis.text.x = element_text(size=12), axis.title.y = element_text(margin = margin(t=0, r= 20, b=0, l=0)))

par(mar = c(10, 10,10,10))

Plot year based on core depth

year <- ggplot() +
  geom_point(data = dp_IMP2, (aes(x = year_IMP, y= depth_IMP2)), na.rm=TRUE, shape=22, fill="paleturquoise4", size=3, color = "paleturquoise4", alpha=0.5) +
  geom_line(data = dp_IMP2, (aes(x = year_IMP, y= depth_IMP2)), na.rm=TRUE, color = "paleturquoise4") +
  stat_smooth(dp_IMP2, mapping=aes(year_IMP, depth_IMP2), method="lm") +
  geom_errorbar(aes(x=year_IMP, xmin=year_IMP-error_IMP, xmax=year_IMP+error_IMP, y=depth_IMP2, ymax=NULL, ymin=NULL), width=.2) +
  labs(x = "Year (CRS model)", y = "Midpoint depth (cm)") +
  scale_x_continuous(breaks=seq(1800,2023, by= 20)) +
  scale_y_reverse(breaks=seq(0,15, by=1), expand=c(0,0)) +
  theme_classic() +
  theme(panel.grid.major.x=element_line(), text = element_text(size = 12), axis.text.y = element_text(size=12), axis.text.x = element_text(size=12))

year

1.4 Final Figure

ggarrange(activity, year, ncol = 1, nrow = 2, common.legend=TRUE,legend=c("top"))

Figure 1.1: [Bosch et al. 2024]: 210-Pb activity (Bq/kg) over the depth of the impact core (CSM-IMP-2, blue squares) and the reference core (CSM-REF-3, green circles). The 210-Pb dates assigned to the midpoint depth of the cores using a constant rate of supply (CRS) dating model are shown in the smaller, inset graph, where pink trend lines represent the second order polynomial regression for each core.