Comparing Archives

Set Up

library(ggplot2)
library(dplyr)
library(tidyr)
DATA_DIR <- here::here("data")

new_archive <- read.csv(here::here(DATA_DIR, "new_new.csv"), stringsAsFactors = FALSE)

ACCESS-ESM1-5

access_old <- read.csv(here::here(DATA_DIR, "ACCESS-ESM1-5_archive_data.csv"), stringsAsFactors = FALSE)
access_old$source <- "old"
new_archive$source <- "new"
new_archive_access <- dplyr::filter(new_archive, model == "ACCESS-ESM1-5")

access_old %>% 
  bind_rows(new_archive_access) %>%
  group_by(model, experiment, source) %>% 
  summarise(count = n_distinct(ensemble)) %>% 
  ungroup %>% 
  tidyr::spread(source, count) %>% 
  knitr::kable(format = "markdown")

`summarise()` regrouping output by 'model', 'experiment' (override with `.groups` argument)

model	experiment	new	old
ACCESS-ESM1-5	1pctCO2	1	NA
ACCESS-ESM1-5	abrupt-4xCO2	2	NA
ACCESS-ESM1-5	ssp119	30	20
ACCESS-ESM1-5	ssp126	30	20
ACCESS-ESM1-5	ssp245	30	20
ACCESS-ESM1-5	ssp370	30	20
ACCESS-ESM1-5	ssp434	30	20
ACCESS-ESM1-5	ssp460	30	20
ACCESS-ESM1-5	ssp534-over	30	20
ACCESS-ESM1-5	ssp585	30	20

ggplot() + 
  geom_point(data = access_old, aes(fx, dx, color = "old")) + 
  geom_point(data = new_archive_access, aes(fx, dx, color = "new"), alpha = 0.5) + 
  theme_bw() + 
  scale_color_manual(values = c("old" = "blue", "new" = "pink")) +
  facet_wrap("experiment", scales = "free") + 
  NULL + 
  labs(title = "ACCESS-ESM1-5")

CanESM5

canesm_old <- read.csv(here::here(DATA_DIR, "CanESM5_archive_data.csv"), stringsAsFactors = FALSE)
canesm_old$source <- "old"
new_archive_canesm <- dplyr::filter(new_archive, model == "CanESM5")

canesm_old %>% 
  bind_rows(new_archive_canesm) %>%
  group_by(model, experiment, source) %>% 
  summarise(count = n_distinct(ensemble)) %>% 
  ungroup %>% 
  tidyr::spread(source, count) %>% 
  knitr::kable(format = "markdown")

`summarise()` regrouping output by 'model', 'experiment' (override with `.groups` argument)

model	experiment	new	old
CanESM5	1pctCO2	6	NA
CanESM5	abrupt-2xCO2	1	NA
CanESM5	abrupt-4xCO2	2	NA
CanESM5	ssp119	65	25
CanESM5	ssp126	65	25
CanESM5	ssp245	65	25
CanESM5	ssp370	65	25
CanESM5	ssp434	65	25
CanESM5	ssp460	65	25
CanESM5	ssp534-over	65	25
CanESM5	ssp585	65	25

ggplot() + 
  geom_point(data = canesm_old, aes(fx, dx, color = "old")) + 
  geom_point(data = new_archive_canesm, aes(fx, dx, color = "new"), alpha = 0.5) + 
  theme_bw() + 
  scale_color_manual(values = c("old" = "blue", "new" = "pink")) +
  facet_wrap("experiment", scales = "free") + 
  NULL + 
  labs(title = "CanESM5")

MIROC6

miroc_old <- read.csv(here::here(DATA_DIR, "MIROC6_archive_data.csv"), stringsAsFactors = FALSE)
miroc_old$source <- "old"
new_archive_miroc <- dplyr::filter(new_archive, model == "MIROC6")

miroc_old %>% 
  bind_rows(new_archive_miroc) %>%
  group_by(model, experiment, source) %>% 
  summarise(count = n_distinct(ensemble)) %>% 
  ungroup %>% 
  tidyr::spread(source, count) %>% 
  knitr::kable(format = "markdown")

`summarise()` regrouping output by 'model', 'experiment' (override with `.groups` argument)

model	experiment	new	old
MIROC6	1pctCO2	1	NA
MIROC6	abrupt-2xCO2	1	NA
MIROC6	abrupt-4xCO2	1	NA
MIROC6	ssp119	50	50
MIROC6	ssp126	50	50
MIROC6	ssp245	50	50
MIROC6	ssp370	50	50
MIROC6	ssp434	50	50
MIROC6	ssp460	50	50
MIROC6	ssp534-over	50	50
MIROC6	ssp585	50	50

ggplot() + 
  geom_point(data = miroc_old, aes(fx, dx, color = "old")) + 
  geom_point(data = new_archive_miroc, aes(fx, dx, color = "new"), alpha = 0.5) + 
  theme_bw() + 
  scale_color_manual(values = c("old" = "blue", "new" = "pink")) +
  facet_wrap("experiment", scales = "free") + 
  NULL + 
  labs(title = "MIROC6")

MPI-ESM1-2-HR

mpi_old <- read.csv(here::here(DATA_DIR, "MPI-ESM1-2-HR_archive_data.csv"), stringsAsFactors = FALSE)
mpi_old$source <- "old"
new_archive_mpi <- dplyr::filter(new_archive, model == "MPI-ESM1-2-HR")

mpi_old %>% 
  bind_rows(new_archive_mpi) %>%
  group_by(model, experiment, source) %>% 
  summarise(count = n_distinct(ensemble)) %>% 
  ungroup %>% 
  tidyr::spread(source, count) %>% 
  knitr::kable(format = "markdown")

`summarise()` regrouping output by 'model', 'experiment' (override with `.groups` argument)

model	experiment	new	old
MPI-ESM1-2-HR	1pctCO2	1	NA
MPI-ESM1-2-HR	abrupt-4xCO2	1	NA
MPI-ESM1-2-HR	ssp119	10	10
MPI-ESM1-2-HR	ssp126	10	10
MPI-ESM1-2-HR	ssp245	10	10
MPI-ESM1-2-HR	ssp370	10	10
MPI-ESM1-2-HR	ssp434	10	10
MPI-ESM1-2-HR	ssp460	10	10
MPI-ESM1-2-HR	ssp534-over	10	10
MPI-ESM1-2-HR	ssp585	10	10

ggplot() + 
  geom_point(data = mpi_old, aes(fx, dx, color = "old")) + 
  geom_point(data = new_archive_mpi, aes(fx, dx, color = "new"), alpha = 0.5) + 
  theme_bw() + 
  scale_color_manual(values = c("old" = "blue", "new" = "pink")) +
  facet_wrap("experiment", scales = "free") + 
  NULL + 
  labs(title = "MPI-ESM1-2-HR")

MPI-ESM1-2-LR

mpi_old <- read.csv(here::here(DATA_DIR, "MPI-ESM1-2-LR_archive_data.csv"), stringsAsFactors = FALSE)
mpi_old$source <- "old"
new_archive_mpi <- dplyr::filter(new_archive, model == "MPI-ESM1-2-LR")

mpi_old %>% 
  bind_rows(new_archive_mpi) %>%
  group_by(model, experiment, source) %>% 
  summarise(count = n_distinct(ensemble)) %>% 
  ungroup %>% 
  tidyr::spread(source, count) %>% 
  knitr::kable(format = "markdown")

`summarise()` regrouping output by 'model', 'experiment' (override with `.groups` argument)

model	experiment	new	old
MPI-ESM1-2-LR	1pctCO2	1	NA
MPI-ESM1-2-LR	abrupt-4xCO2	1	NA
MPI-ESM1-2-LR	ssp119	10	10
MPI-ESM1-2-LR	ssp126	10	10
MPI-ESM1-2-LR	ssp245	10	10
MPI-ESM1-2-LR	ssp370	10	10
MPI-ESM1-2-LR	ssp434	10	10
MPI-ESM1-2-LR	ssp460	10	10
MPI-ESM1-2-LR	ssp534-over	10	10
MPI-ESM1-2-LR	ssp585	10	10

ggplot() + 
  geom_point(data = mpi_old, aes(fx, dx, color = "old")) + 
  geom_point(data = new_archive_mpi, aes(fx, dx, color = "new"), alpha = 0.5) + 
  theme_bw() + 
  scale_color_manual(values = c("old" = "blue", "new" = "pink")) +
  facet_wrap("experiment", scales = "free") + 
  NULL + 
  labs(title = "MPI-ESM1-2-HR")

NorCPM1

norcpm1_old <- read.csv(here::here(DATA_DIR, "NorCPM1_archive_data.csv"), stringsAsFactors = FALSE)
norcpm1_old$source <- "old"
new_archive_norcpm1 <- dplyr::filter(new_archive, model == "NorCPM1")

norcpm1_old %>% 
  bind_rows(new_archive_norcpm1) %>%
  group_by(model, experiment, source) %>% 
  summarise(count = n_distinct(ensemble)) %>% 
  ungroup %>% 
  tidyr::spread(source, count) %>% 
  knitr::kable(format = "markdown")

`summarise()` regrouping output by 'model', 'experiment' (override with `.groups` argument)

model	experiment	new	old
NorCPM1	1pctCO2	1	NA
NorCPM1	abrupt-4xCO2	1	NA
NorCPM1	ssp119	29	29
NorCPM1	ssp126	29	29
NorCPM1	ssp245	29	29
NorCPM1	ssp370	29	29
NorCPM1	ssp434	29	29
NorCPM1	ssp460	29	29
NorCPM1	ssp534-over	29	29
NorCPM1	ssp585	29	29

ggplot() + 
  geom_point(data = norcpm1_old, aes(fx, dx, color = "old")) + 
  geom_point(data = new_archive_norcpm1, aes(fx, dx, color = "new"), alpha = 0.5) + 
  theme_bw() + 
  scale_color_manual(values = c("old" = "blue", "new" = "pink")) +
  facet_wrap("experiment", scales = "free") + 
  NULL + 
  labs(title = "NorCPM1")

UKESM1-0-LL

ukesm_old <- read.csv(here::here(DATA_DIR, "UKESM1-0-LL_archive_data.csv"), stringsAsFactors = FALSE)
ukesm_old$source <- "old"
new_archive_ukesm <- dplyr::filter(new_archive, model == "UKESM1-0-LL")

ukesm_old %>% 
  bind_rows(new_archive_ukesm) %>%
  group_by(model, experiment, source) %>% 
  summarise(count = n_distinct(ensemble)) %>% 
  ungroup %>% 
  tidyr::spread(source, count) %>% 
  knitr::kable(format = "markdown")

`summarise()` regrouping output by 'model', 'experiment' (override with `.groups` argument)

model	experiment	new	old
UKESM1-0-LL	1pctCO2	4	NA
UKESM1-0-LL	abrupt-4xCO2	1	NA
UKESM1-0-LL	ssp119	19	18
UKESM1-0-LL	ssp126	19	18
UKESM1-0-LL	ssp245	19	18
UKESM1-0-LL	ssp370	19	18
UKESM1-0-LL	ssp434	19	18
UKESM1-0-LL	ssp460	19	18
UKESM1-0-LL	ssp534-over	19	18
UKESM1-0-LL	ssp585	19	18

ggplot() + 
  geom_point(data = ukesm_old, aes(fx, dx, color = "old")) + 
  geom_point(data = new_archive_ukesm, aes(fx, dx, color = "new"), alpha = 0.5) + 
  theme_bw() + 
  scale_color_manual(values = c("old" = "blue", "new" = "pink")) +
  facet_wrap("experiment", scales = "free") + 
  NULL + 
  labs(title = "UKESM1-0-LL")

New Experiments

new_exps <- c("1pctCO2", "abrupt-4xCO2", "abrupt-2xCO2")
old_exps <- setdiff(new_exps, unique(new_archive$experiment))
new_archive %>% 
  mutate(exp = if_else(experiment %in% new_exps, '1%, x4 CO2 & x2 CO2', 'historical + ssps')) -> 
  data
mods <- unique(data$model)
data$exp <- factor(data$exp, c('1%, x4 CO2 & x2 CO2', 'historical + ssps'), ordered = TRUE)
data %>% 
  filter(model %in% mods[1:16]) %>% 
  ggplot(aes(fx, dx, color = exp)) + 
  geom_point() + 
  facet_wrap("model", scales = "free") + 
  theme_bw() + 
  scale_color_manual(values = c('1%, x4 CO2 & x2 CO2' = 'blue', 'historical + ssps' = 'grey'))

data %>% 
  filter(model %in% mods[17:27]) %>% 
  ggplot(aes(fx, dx, color = exp)) + 
  geom_point() + 
  facet_wrap("model", scales = "free") + 
  theme_bw() + 
  scale_color_manual(values = c('1%, x4 CO2 & x2 CO2' = 'blue', 'historical + ssps' = 'grey'))

data %>% 
  filter(model %in% mods[28:33]) %>% 
  ggplot(aes(fx, dx, color = exp)) + 
  geom_point() + 
  facet_wrap("model", scales = "free") + 
  theme_bw() + 
  scale_color_manual(values = c('1%, x4 CO2 & x2 CO2' = 'blue', 'historical + ssps' = 'grey'))