Bias

Precipitation

French Guiana

Code

cmip6 <- rast("results/evaluation/bias/French-Guiana_CMIP6_chirps.nc")[[1:12]] %>% # nolint
  median()
ggplot() +
  tidyterra::geom_spatraster(data = cmip6) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CMIP6 vs CHIRPS", "Median across month") +
  scale_fill_gradient2("",
    low = "#543005", mid = "white",
    high = "#003c30", midpoint = 0, na.value = "black",
    limits = c(-100, 100)
  )

Code

cordex <- rast("results/evaluation/bias/French-Guiana_CORDEX_chirps.nc")[[1:12]] %>% # nolint
  median()
ggplot() +
  tidyterra::geom_spatraster(data = cordex) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CORDEX vs CHIRPS", "Median across month") +
  scale_fill_gradient2("",
    low = "#543005", mid = "white",
    high = "#003c30", midpoint = 0, na.value = "black",
    limits = c(-100, 100)
  )

Ivory Coast

Code

cmip6 <- rast("results/evaluation/bias/Côte-d'Ivoire_CMIP6_chirps.nc")[[1:12]] %>% # nolint
  median()
ggplot() +
  tidyterra::geom_spatraster(data = cmip6) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CMIP6 vs CHIRPS", "Median across month") +
  scale_fill_gradient2("",
    low = "#543005", mid = "white",
    high = "#003c30", midpoint = 0, na.value = "black",
    limits = c(-100, 100)
  )

Code

cordex <- rast("results/evaluation/bias/Côte-d'Ivoire_CORDEX_chirps.nc")[[1:12]] %>% # nolint
  median()
ggplot() +
  tidyterra::geom_spatraster(data = cordex) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CORDEX vs CHIRPS", "Median across month") +
  scale_fill_gradient2("",
    low = "#543005", mid = "white",
    high = "#003c30", midpoint = 0, na.value = "black",
    limits = c(-100, 100)
  )

New Caledonia

Code

cmip6 <- rast("results/evaluation/bias/New-Caledonia_CMIP6_chirps.nc")[[1:12]] %>% # nolint
  median() %>%
  crop(c(163, 168.5, -23, -19))
ggplot() +
  tidyterra::geom_spatraster(data = cmip6) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CMIP6 vs CHIRPS", "Median across month") +
  scale_fill_gradient2("",
    low = "#543005", mid = "white",
    high = "#003c30", midpoint = 0, na.value = "black",
    limits = c(-200, 200)
  )

Code

cordex <- rast("results/evaluation/bias/New-Caledonia_CORDEX_chirps.nc")[[1:12]] %>% # nolint
  median() %>%
  crop(c(163, 168.5, -23, -19))
ggplot() +
  tidyterra::geom_spatraster(data = cordex) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CORDEX vs CHIRPS", "Median across month") +
  scale_fill_gradient2("",
    low = "#543005", mid = "white",
    high = "#003c30", midpoint = 0, na.value = "black",
    limits = c(-200, 200)
  )

Temperature

French Guiana

Code

cmip6 <- rast("results/evaluation/bias/French-Guiana_CMIP6_gshtd.nc")[[1:12]] %>% # nolint
  median()
ggplot() +
  tidyterra::geom_spatraster(data = cmip6) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CMIP6 vs GSHTD", "Median across month") +
  scale_fill_gradient2("",
    low = "#053061", mid = "white", high = "#67001f",
    midpoint = 0, na.value = "black",
    limits = c(-4, 4)
  )

Code

cordex <- rast("results/evaluation/bias/French-Guiana_CORDEX_gshtd.nc")[[1:12]] %>% # nolint
  median()
ggplot() +
  tidyterra::geom_spatraster(data = cordex) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CORDEX vs GSHTD", "Median across month") +
  scale_fill_gradient2("",
    low = "#053061", mid = "white", high = "#67001f",
    midpoint = 0, na.value = "black",
    limits = c(-4, 4)
  )

Ivory Coast

Code

cmip6 <- rast("results/evaluation/bias/Côte-d'Ivoire_CMIP6_gshtd.nc")[[1:12]] %>% # nolint
  median()
ggplot() +
  tidyterra::geom_spatraster(data = cmip6) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CMIP6 vs GSHTD", "Median across month") +
  scale_fill_gradient2("",
    low = "#053061", mid = "white", high = "#67001f",
    midpoint = 0, na.value = "black",
    limits = c(-4, 4)
  )

Code

cordex <- rast("results/evaluation/bias/Côte-d'Ivoire_CORDEX_gshtd.nc")[[1:12]] %>% # nolint
  median()
ggplot() +
  tidyterra::geom_spatraster(data = cordex) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CORDEX vs GSHTD", "Median across month") +
  scale_fill_gradient2("",
    low = "#053061", mid = "white", high = "#67001f",
    midpoint = 0, na.value = "black",
    limits = c(-4, 4)
  )

New Caledonia

Code

cmip6 <- rast("results/evaluation/bias/New-Caledonia_CMIP6_gshtd.nc")[[1:12]] %>% # nolint
  median() %>%
  crop(c(163, 168.5, -23, -19))
ggplot() +
  tidyterra::geom_spatraster(data = cmip6) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CMIP6 vs GSHTD", "Median across month") +
  scale_fill_gradient2("",
    low = "#053061", mid = "white", high = "#67001f",
    midpoint = 0, na.value = "black",
    limits = c(-2, 2)
  )

Code

cordex <- rast("results/evaluation/bias/New-Caledonia_CORDEX_gshtd.nc")[[1:12]] %>% # nolint
  median() %>%
  crop(c(163, 168.5, -23, -19))
ggplot() +
  tidyterra::geom_spatraster(data = cordex) +
  theme_bw() +
  theme(
    panel.border = element_blank(), panel.grid = element_blank(),
    panel.background = element_rect(fill = "black")
  ) +
  ggtitle("CORDEX vs GSHTD", "Median across month") +
  scale_fill_gradient2("",
    low = "#053061", mid = "white", high = "#67001f",
    midpoint = 0, na.value = "black",
    limits = c(-2, 2)
  )