scenarios

library(tidyverse)
library(glue)
library(here)
library(raster)
library(knitr)
#library(bbnj)
devtools::load_all(here("../bbnj")) # devtools::install_local(force=T)
select = dplyr::select

dir_scenarios <- here("../bbnj/inst/app/www/scenarios")
dir_diffs     <- here("../bbnj/inst/app/www/scenarios/diffs")
scenarios <- list(
  s1  = "s01a.bio.now.mol50km",
  s2a = "s02a.bio.alltime.mol50km",
  s2b = "s02b.bio.future.mol50km",
  s3  = "s03a.biofish.now.mol50km",
  s4  = "s04a.biofish.alltime.mol50km")

diff_tbl <- function(x,...){
  x %>% 
    mutate(
      pct = sprintf("%.01f%%", pct*100)) %>% 
    kable(
      format.args = list(big.mark=','), 
      align = c("r","l","r","r","r"), ...)
}

tibble(
  s        = names(scenarios),
  scenario = unlist(scenarios)) %>% 
  kable()
s scenario
s1 s01a.bio.now.mol50km
s2a s02a.bio.alltime.mol50km
s2b s02b.bio.future.mol50km
s3 s03a.biofish.now.mol50km
s4 s04a.biofish.alltime.mol50km

Figure 3: s2a - s4

Biodiversity now & future - Biodiversity now & future with fishing as cost

d <- scenarios_diff(scenarios[c("s2a", "s4")], dir_scenarios, dir_diffs)
diff_tbl(d$tbl)
value label ncells area_km2 pct
-1 loss 3,044 7,642,819 12.7%
0 same 17,504 43,948,719 73.0%
1 gain 3,424 8,596,916 14.3%

Figure 4: s3 - s4

Biodiversity now with fishing as cost - Biodiversity now & future with fishing as cost

d <- scenarios_diff(scenarios[c("s3", "s4")], dir_scenarios, dir_diffs)
diff_tbl(d$tbl)
value label ncells area_km2 pct
-1 loss 583 1,463,786 2.7%
0 same 20,331 51,046,698 94.5%
1 gain 597 1,498,937 2.8%

For discussion: s1 - s2b

Biodiversity now (only) - Biodiversity future (only)

d <- scenarios_diff(scenarios[c("s1", "s2b")], dir_scenarios, dir_diffs)
diff_tbl(d$tbl)
value label ncells area_km2 pct
-1 loss 6,535 16,407,957 24.8%
0 same 14,003 35,158,473 53.2%
1 gain 5,785 14,524,871 22.0%