Last updated: 2021-11-19

Checks: 7 0

Knit directory: emlr_obs_v_XXX/

This reproducible R Markdown analysis was created with workflowr (version 1.6.2). The Checks tab describes the reproducibility checks that were applied when the results were created. The Past versions tab lists the development history.


Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

The command set.seed(20200707) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

Great job! Recording the operating system, R version, and package versions is critical for reproducibility.

Nice! There were no cached chunks for this analysis, so you can be confident that you successfully produced the results during this run.

Great job! Using relative paths to the files within your workflowr project makes it easier to run your code on other machines.

Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility.

The results in this page were generated with repository version 5c01583. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:


Ignored files:
    Ignored:    .Rhistory
    Ignored:    .Rproj.user/

Unstaged changes:
    Modified:   code/Workflowr_project_managment.R

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.


These are the previous versions of the repository in which changes were made to the R Markdown (analysis/analysis_anomalous_changes.Rmd) and HTML (docs/analysis_anomalous_changes.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

File Version Author Date Message
html 98d9e33 jens-daniel-mueller 2021-11-11 Build site.
html d3cb92d jens-daniel-mueller 2021-11-08 Build site.
html 3879a6d jens-daniel-mueller 2021-11-08 Build site.
html 0e032dc jens-daniel-mueller 2021-11-04 Build site.
html abcd28f jens-daniel-mueller 2021-11-02 Build site.
html 290c8fc jens-daniel-mueller 2021-11-02 Build site.
html e02acc9 jens-daniel-mueller 2021-11-01 Build site.
html 58da811 jens-daniel-mueller 2021-11-01 Build site.
html 2781a97 jens-daniel-mueller 2021-10-29 Build site.
html 973192c jens-daniel-mueller 2021-10-28 Build site.
html f792cd8 jens-daniel-mueller 2021-10-27 Build site.
html 22750db jens-daniel-mueller 2021-10-27 Build site.
html 3992e41 jens-daniel-mueller 2021-10-27 Build site.
html f9280c7 jens-daniel-mueller 2021-10-21 Build site.
html bfe90e6 jens-daniel-mueller 2021-10-21 Build site.
html 78a956e jens-daniel-mueller 2021-10-21 Build site.
html 1ba14cf jens-daniel-mueller 2021-10-20 Build site.
html dc0723b jens-daniel-mueller 2021-10-19 Build site.
html 722483f jens-daniel-mueller 2021-10-18 Build site.
html cd0ebd8 jens-daniel-mueller 2021-10-18 Build site.
html aedafcb jens-daniel-mueller 2021-10-17 Build site.
html 046c556 jens-daniel-mueller 2021-10-07 Build site.
html ac73c23 jens-daniel-mueller 2021-10-07 Build site.
html 0bd4ea2 jens-daniel-mueller 2021-10-07 Build site.
html 6969c18 jens-daniel-mueller 2021-10-07 check data distribution
html 06afaca jens-daniel-mueller 2021-09-23 Build site.
html 410dd46 jens-daniel-mueller 2021-09-22 Build site.
html 5a36231 jens-daniel-mueller 2021-09-22 Build site.
html c0bbbe4 jens-daniel-mueller 2021-09-22 Build site.
html 9eac969 jens-daniel-mueller 2021-08-30 Build site.
html 3a584b7 jens-daniel-mueller 2021-08-30 Build site.
html 6e17588 jens-daniel-mueller 2021-08-24 Build site.
html deb425c jens-daniel-mueller 2021-08-24 Build site.
html 534854b jens-daniel-mueller 2021-08-19 Build site.
html 77f4ba7 jens-daniel-mueller 2021-08-19 Build site.
html ece96df jens-daniel-mueller 2021-08-19 Build site.
html 27c99b8 jens-daniel-mueller 2021-08-19 Build site.
html a03f2f0 jens-daniel-mueller 2021-08-18 Build site.
html 9335b31 jens-daniel-mueller 2021-08-10 Build site.
html 9943b45 jens-daniel-mueller 2021-08-10 Build site.
html 0b00a2b jens-daniel-mueller 2021-08-09 Build site.
html 755c6b1 jens-daniel-mueller 2021-08-09 Build site.
html cd8e0d5 jens-daniel-mueller 2021-08-06 Build site.
html 15773a0 jens-daniel-mueller 2021-08-06 included calculation of revelle factor
html da61d1a jens-daniel-mueller 2021-08-06 Build site.
html 340d731 jens-daniel-mueller 2021-08-06 Build site.
html 71546e4 jens-daniel-mueller 2021-08-06 test with stricter CANYON-B filtering
html 29444a1 jens-daniel-mueller 2021-08-05 Build site.
html 42e80c0 jens-daniel-mueller 2021-08-04 Build site.
html 48f6eed jens-daniel-mueller 2021-08-04 Build site.
html 7bc0dd5 jens-daniel-mueller 2021-08-04 Build site.
Rmd 537e268 jens-daniel-mueller 2021-08-04 test with intermediate CANYON-B offset threshold
html 81a46a4 jens-daniel-mueller 2021-08-03 Build site.
html b88c61b jens-daniel-mueller 2021-08-03 Build site.
html 0f0d5e5 jens-daniel-mueller 2021-08-03 Build site.
html a53656d jens-daniel-mueller 2021-08-03 Build site.
html 88f7356 jens-daniel-mueller 2021-08-02 Build site.
html d759279 jens-daniel-mueller 2021-08-02 Build site.
html 127b801 jens-daniel-mueller 2021-07-24 Build site.
html 1f65ef1 jens-daniel-mueller 2021-07-23 Build site.
html 912d90e jens-daniel-mueller 2021-07-23 Build site.
html 2477316 jens-daniel-mueller 2021-07-23 rebuild: surface dcant mapping seperate
html c9ccc00 jens-daniel-mueller 2021-07-22 Build site.
html 426b2df jens-daniel-mueller 2021-07-21 Build site.
html 971ce87 jens-daniel-mueller 2021-07-13 Build site.
html c68ae3a jens-daniel-mueller 2021-07-13 Build site.
Rmd 620afd3 jens-daniel-mueller 2021-07-13 complete revision
html 910d64a jens-daniel-mueller 2021-07-02 Build site.
html 1cbf907 jens-daniel-mueller 2021-07-02 Build site.
html 955f5dc jens-daniel-mueller 2021-06-29 Build site.
html 9480f94 jens-daniel-mueller 2021-06-29 Build site.
html d054c39 jens-daniel-mueller 2021-06-24 Build site.
html 4fcd066 jens-daniel-mueller 2021-06-24 Build site.
html aedd245 jens-daniel-mueller 2021-06-23 Build site.
Rmd 31bffa7 jens-daniel-mueller 2021-06-23 revised slab bias and write files
html 00dcc97 jens-daniel-mueller 2021-06-23 Build site.
Rmd 891392a jens-daniel-mueller 2021-06-23 revised bias labels
html b1ddc9b jens-daniel-mueller 2021-06-18 Build site.
Rmd 0fb43db jens-daniel-mueller 2021-06-18 write anomaly files
html c0a1abb jens-daniel-mueller 2021-06-18 Build site.
Rmd e262f99 jens-daniel-mueller 2021-06-18 revise observation anomalies
html 9a8d446 jens-daniel-mueller 2021-06-18 Build site.
Rmd 0d50381 jens-daniel-mueller 2021-06-18 revised bias label
html 3b6387f jens-daniel-mueller 2021-06-18 Build site.
html e03e1ac jens-daniel-mueller 2021-06-16 Build site.
Rmd cdc534f jens-daniel-mueller 2021-06-16 reverted nss sign
html bf487c8 jens-daniel-mueller 2021-06-16 Build site.
Rmd 978e96e jens-daniel-mueller 2021-06-16 non steady state changes added
html b261191 jens-daniel-mueller 2021-06-15 Build site.
Rmd c882921 jens-daniel-mueller 2021-06-15 vertical bias profile
html 7a8d5f8 jens-daniel-mueller 2021-06-15 Build site.
Rmd 3695731 jens-daniel-mueller 2021-06-15 included budget assesment
html 7155166 jens-daniel-mueller 2021-06-15 Build site.
Rmd 03a41f4 jens-daniel-mueller 2021-06-15 included budget assesment
html 71c4934 jens-daniel-mueller 2021-06-15 Build site.
Rmd 0e41426 jens-daniel-mueller 2021-06-15 anomaly inventory maps added
html 3052a6c jens-daniel-mueller 2021-06-15 Build site.
html 694d661 jens-daniel-mueller 2021-06-15 Build site.
Rmd 37f3208 jens-daniel-mueller 2021-06-15 test with pre 2000 and post 2010 data
html 3cb64bd jens-daniel-mueller 2021-06-15 Build site.
html fc03c25 jens-daniel-mueller 2021-06-14 Build site.
Rmd 7f89272 jens-daniel-mueller 2021-06-14 revised delta cant patterns for mod data
html 59a04e1 jens-daniel-mueller 2021-06-14 Build site.
Rmd 807c560 jens-daniel-mueller 2021-06-14 revised delta cant patterns for mod data
html 6d074b8 jens-daniel-mueller 2021-06-14 Build site.
Rmd b9b9a1a jens-daniel-mueller 2021-06-14 revised zonal mean section plots for mod data
html c6b3da6 jens-daniel-mueller 2021-06-14 Build site.
html 48c73fc jens-daniel-mueller 2021-06-14 Build site.
html 7a0e665 jens-daniel-mueller 2021-06-14 Build site.
Rmd 4a3ec7d jens-daniel-mueller 2021-06-14 revised zonal mean section plots for mod data
html b5744ae jens-daniel-mueller 2021-06-14 Build site.
Rmd 12df5f5 jens-daniel-mueller 2021-06-14 revised zonal mean section plots for mod data
html 439ee80 jens-daniel-mueller 2021-06-11 Build site.
Rmd c22eef5 jens-daniel-mueller 2021-06-11 revised scaling of disequilibrium ratio
html 33ffcab jens-daniel-mueller 2021-06-10 Build site.
html 7e1f407 jens-daniel-mueller 2021-06-10 Build site.
html 2cbe18c jens-daniel-mueller 2021-06-10 added zonal mean section control plots
html 69c79d0 jens-daniel-mueller 2021-06-08 Build site.
html 1772903 jens-daniel-mueller 2021-06-07 Build site.
html b9c3ed7 jens-daniel-mueller 2021-06-04 Build site.
Rmd 0f99f02 jens-daniel-mueller 2021-06-04 update anomalous changes
html 594ed9a jens-daniel-mueller 2021-06-04 Build site.
html db7df0e jens-daniel-mueller 2021-06-04 rebuild without overlapping eras
Rmd 71f97a3 jens-daniel-mueller 2021-06-04 rebuild without overlapping eras
html 0ad4865 jens-daniel-mueller 2021-06-04 Build site.
Rmd e331dd3 jens-daniel-mueller 2021-06-04 update anomalous changes
html ebd925f jens-daniel-mueller 2021-06-04 Build site.
Rmd 3459fd2 jens-daniel-mueller 2021-06-04 update anomalous changes
html 2edc791 jens-daniel-mueller 2021-06-04 Build site.
html 1c87a03 jens-daniel-mueller 2021-06-04 Build site.
Rmd b48d553 jens-daniel-mueller 2021-06-04 clean anomalous changes
html a656287 jens-daniel-mueller 2021-06-04 Build site.
Rmd 0e89bd2 jens-daniel-mueller 2021-06-04 clean anomalous changes
html 7e8d839 jens-daniel-mueller 2021-06-03 Build site.
Rmd a81f578 jens-daniel-mueller 2021-06-03 include anomalous changes
html 207339d jens-daniel-mueller 2021-06-03 Build site.
html 315710b jens-daniel-mueller 2021-06-03 include anomalous changes
html d612df8 jens-daniel-mueller 2021-06-03 Build site.
Rmd b8fb70c jens-daniel-mueller 2021-06-03 include anomalous changes
html 406d740 jens-daniel-mueller 2021-06-03 Build site.
Rmd 98e5d73 jens-daniel-mueller 2021-06-03 include anomalous changes
html 7802dd6 jens-daniel-mueller 2021-06-03 Build site.
Rmd dced077 jens-daniel-mueller 2021-06-03 include anomalous changes
html a496b5f jens-daniel-mueller 2021-06-03 Build site.
Rmd 90863c0 jens-daniel-mueller 2021-06-03 include anomalous changes

1 Version ID

The results displayed on this site correspond to the Version_ID: v_XXX

2 Data sources

Required are:

  • tcant from Sabine 2004 (S04)

  • tcant from Gruber 2019 (G19)

  • annual mean atmospheric pCO2

  • Mean eMLR-dcant per grid cell (lat, lon, depth)

dcant_3d_ss <- read_csv(paste(path_version_data,
                             "projected_tcant_3d_tref.csv",
                             sep = ""))

dcant_3d_emlr <-
  read_csv(paste(path_version_data,
                 "dcant_3d.csv",
                 sep = ""))

co2_atm <-
  read_csv(paste(path_preprocessing,
                 "co2_atm.csv",
                 sep = ""))

tref <-
  read_csv(paste(path_version_data,
                 "tref.csv",
                 sep = ""))


dcant_3d_mod_truth <-
  read_csv(paste(path_version_data,
                 "dcant_3d_mod_truth.csv", sep = ""))

dcant_3d_mod_truth_cc <-
  read_csv(paste(path_version_data,
                 "dcant_3d_mod_truth_cc.csv", sep = ""))

3 Prepare delta dcant

3.1 Steady state

# join with basinmask
dcant_3d_ss <- inner_join(dcant_3d_ss, basinmask)

# calculate delta dcant as the difference between total dcant at tref
# this includes the G19 anomaly
dcant_3d_ss <- dcant_3d_ss %>%
  pivot_wider(names_from = era,
              values_from = "tcant_pos") %>% 
  mutate(dcant_pos = !!sym(tref$era[2]) - !!sym(tref$era[1])) %>% 
  select(-c(!!sym(tref$era[2]), !!sym(tref$era[1]))) %>%
  drop_na()

3.2 Model truth

dcant_3d_mod_truth <- dcant_3d_mod_truth %>% 
  select(lon, lat, depth, basin_AIP, dcant_pos, data_source)

dcant_3d_mod_truth_cc <- dcant_3d_mod_truth_cc %>% 
  select(lon, lat, depth, basin_AIP, dcant_pos, data_source)

3.3 eMLR(C*)

dcant_3d_emlr <- dcant_3d_emlr %>% 
  select(data_source, lon, lat, depth, basin_AIP, dcant_pos)

4 Join data

4.1 Model

dcant_mod_3d <- bind_rows(
  dcant_3d_emlr %>%
    filter(data_source == "mod") %>%
    select(-data_source) %>%
    mutate(estimate = "emlr"),
  dcant_3d_mod_truth %>%
    filter(data_source == "mod_truth") %>%
    select(-data_source) %>%
    mutate(estimate = "mod_truth_vc"),
  dcant_3d_mod_truth_cc %>%
    filter(data_source == "mod_truth_cc") %>%
    select(-data_source) %>%
    mutate(estimate = "mod_truth_cc"),
  dcant_3d_ss %>%
    filter(data_source == "mod") %>%
    select(-data_source) %>%
    mutate(estimate = "ss_tss")
)

4.2 Observations

dcant_obs_3d <- bind_rows(
  dcant_3d_emlr %>% 
    filter(data_source == "obs") %>% 
    select(-data_source) %>% 
    mutate(estimate = "emlr"),
  dcant_3d_ss %>% 
    filter(data_source == "obs") %>% 
    select(-data_source) %>% 
    mutate(estimate = "ss_tss")
)

5 Model

5.1 Bias correlation

dcant_mod_3d_bias <- dcant_mod_3d %>%
  pivot_wider(names_from = estimate,
              values_from = dcant_pos) %>%
  drop_na()


dcant_mod_3d_bias <- dcant_mod_3d_bias %>% 
  mutate(emlr_anom = emlr - ss_tss,
         emlr_bias = emlr - mod_truth_vc,
         nss_mod_truth = mod_truth_vc - mod_truth_cc,
         nss_projection = mod_truth_vc - ss_tss) %>% 
  select(-c(emlr:ss_tss))

axis_limit <- dcant_mod_3d_bias %>% 
  pivot_longer(emlr_anom:nss_projection,
               names_to = "estimate",
               values_to = "bias") %>% 
  pull(bias) %>% 
  abs() %>% 
  max()

dcant_mod_3d_bias %>%
  ggplot(aes(emlr_bias, emlr_anom)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap( ~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
00dcc97 jens-daniel-mueller 2021-06-23
9a8d446 jens-daniel-mueller 2021-06-18
3b6387f jens-daniel-mueller 2021-06-18
3052a6c jens-daniel-mueller 2021-06-15
694d661 jens-daniel-mueller 2021-06-15
59a04e1 jens-daniel-mueller 2021-06-14
dcant_mod_3d_bias %>% 
  ggplot(aes(emlr_anom, nss_mod_truth)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap(~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
00dcc97 jens-daniel-mueller 2021-06-23
9a8d446 jens-daniel-mueller 2021-06-18
3b6387f jens-daniel-mueller 2021-06-18
3052a6c jens-daniel-mueller 2021-06-15
694d661 jens-daniel-mueller 2021-06-15
59a04e1 jens-daniel-mueller 2021-06-14
dcant_mod_3d_bias %>% 
  ggplot(aes(nss_mod_truth, nss_projection)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap(~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
c68ae3a jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
00dcc97 jens-daniel-mueller 2021-06-23
9a8d446 jens-daniel-mueller 2021-06-18
3b6387f jens-daniel-mueller 2021-06-18
3052a6c jens-daniel-mueller 2021-06-15
694d661 jens-daniel-mueller 2021-06-15
59a04e1 jens-daniel-mueller 2021-06-14

5.1.1 Vertical patterns

dcant_mod_3d_bias_profile <- dcant_mod_3d_bias %>% 
  pivot_longer(emlr_anom:nss_projection,
               names_to = "estimate",
               values_to = "bias") %>% 
  group_by(depth, estimate, basin_AIP) %>% 
  summarise(bias = mean(bias)) %>% 
  ungroup()

dcant_mod_3d_bias_profile %>% 
  ggplot(aes(bias, depth, col=estimate)) +
  geom_path() +
  scale_y_reverse() +
  facet_grid(~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
00dcc97 jens-daniel-mueller 2021-06-23
9a8d446 jens-daniel-mueller 2021-06-18
3b6387f jens-daniel-mueller 2021-06-18
e03e1ac jens-daniel-mueller 2021-06-16
bf487c8 jens-daniel-mueller 2021-06-16
b261191 jens-daniel-mueller 2021-06-15
dcant_mod_3d_bias_profile <- dcant_mod_3d_bias %>% 
  pivot_longer(emlr_anom:nss_projection,
               names_to = "estimate",
               values_to = "bias") %>% 
  group_by(depth, estimate) %>% 
  summarise(bias = mean(bias)) %>% 
  ungroup()

dcant_mod_3d_bias_profile %>% 
  ggplot(aes(bias, depth, col=estimate)) +
  geom_path() +
  scale_y_reverse()

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
910d64a jens-daniel-mueller 2021-07-02
1cbf907 jens-daniel-mueller 2021-07-02
955f5dc jens-daniel-mueller 2021-06-29
9480f94 jens-daniel-mueller 2021-06-29
d054c39 jens-daniel-mueller 2021-06-24
4fcd066 jens-daniel-mueller 2021-06-24
00dcc97 jens-daniel-mueller 2021-06-23
9a8d446 jens-daniel-mueller 2021-06-18
3b6387f jens-daniel-mueller 2021-06-18
e03e1ac jens-daniel-mueller 2021-06-16
bf487c8 jens-daniel-mueller 2021-06-16
b261191 jens-daniel-mueller 2021-06-15
rm(dcant_mod_3d_bias_profile)

5.2 Zonal sections

5.2.1 delta dcant

# calculate zonal mean section
dcant_mod_zonal <- dcant_mod_3d %>%
  group_by(estimate) %>%
  nest() %>%
  mutate(section = map(.x = data, ~m_zonal_mean_sd(.x))) %>%
  select(-data) %>%
  unnest(section)

dcant_mod_zonal <- dcant_mod_zonal %>% 
  rename(dcant_pos = dcant_pos_mean)

5.2.2 Bias

dcant_mod_zonal_bias <- dcant_mod_zonal %>%
  select(-dcant_pos_sd) %>% 
  pivot_wider(names_from = estimate,
              values_from = dcant_pos) %>% 
  drop_na()


dcant_mod_zonal_bias <- dcant_mod_zonal_bias %>% 
  mutate(emlr_anom = emlr - ss_tss,
         emlr_bias = emlr - mod_truth_vc,
         nss_mod_truth = mod_truth_vc - mod_truth_cc,
         nss_projection = mod_truth_vc - ss_tss) %>% 
  select(-c(emlr:ss_tss))

dcant_mod_zonal_bias_long <- dcant_mod_zonal_bias %>% 
  pivot_longer(emlr_anom:nss_projection,
               names_to = "estimate",
               values_to = "dcant_pos_bias")
  


dcant_mod_zonal_bias_long %>% 
  group_by(basin_AIP, estimate) %>%
  group_split() %>% 
  # head(1) %>%
  map( ~ p_section_zonal(
    df = .x,
    var = "dcant_pos_bias",
    col = "divergent",
    breaks = params_global$breaks_dcant_offset,
    plot_slabs = "n",
    subtitle_text = paste("Basin:", .x$basin_AIP, "| Estimate:", .x$estimate)
  ))
[[1]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[2]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[3]]

Version Author Date
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
0b00a2b jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
88f7356 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c68ae3a jens-daniel-mueller 2021-07-13

[[4]]

Version Author Date
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
0b00a2b jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
88f7356 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c68ae3a jens-daniel-mueller 2021-07-13

[[5]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[6]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[7]]

Version Author Date
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
0b00a2b jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
88f7356 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c68ae3a jens-daniel-mueller 2021-07-13

[[8]]

Version Author Date
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
0b00a2b jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
88f7356 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c68ae3a jens-daniel-mueller 2021-07-13

[[9]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[10]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[11]]

Version Author Date
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
0b00a2b jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
88f7356 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c68ae3a jens-daniel-mueller 2021-07-13

[[12]]

Version Author Date
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
0b00a2b jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
88f7356 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c68ae3a jens-daniel-mueller 2021-07-13
axis_limit <- dcant_mod_zonal_bias %>% 
  pivot_longer(emlr_anom:nss_projection,
               names_to = "estimate",
               values_to = "bias") %>% 
  pull(bias) %>% 
  abs() %>% 
  max()

dcant_mod_zonal_bias %>%
  ggplot(aes(emlr_bias, emlr_anom)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap( ~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
dcant_mod_zonal_bias %>% 
  ggplot(aes(emlr_anom, nss_mod_truth)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap(~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
dcant_mod_zonal_bias %>% 
  ggplot(aes(nss_mod_truth, nss_projection)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap(~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

5.3 Inventories

5.3.1 dcant

dcant_mod_inv <- dcant_mod_3d %>%
  mutate(dcant = 0) %>% 
  group_by(estimate) %>%
  nest() %>%
  mutate(inv = map(.x = data, ~m_dcant_inv(.x))) %>%
  select(-data) %>%
  unnest(inv) %>% 
  select(-c(dcant))

5.3.2 Bias

dcant_mod_inv_bias <- dcant_mod_inv %>%
  filter(inv_depth == params_global$inventory_depth_standard) %>% 
  select(-inv_depth) %>% 
  pivot_wider(names_from = estimate,
              values_from = dcant_pos) %>% 
  drop_na()


dcant_mod_inv_bias <- dcant_mod_inv_bias %>% 
  mutate(emlr_anom = emlr - ss_tss,
         emlr_bias = emlr - mod_truth_vc,
         nss_mod_truth = mod_truth_vc - mod_truth_cc,
         nss_projection = mod_truth_vc - ss_tss) %>% 
  select(-c(emlr:ss_tss))


dcant_mod_inv_bias_long <- dcant_mod_inv_bias %>% 
  pivot_longer(emlr_anom:nss_projection,
               values_to = "dcant_pos_bias",
               names_to = "estimate")
  


dcant_mod_inv_bias_long %>% 
  group_by(estimate) %>%
  group_split() %>% 
  # head(1) %>%
  map( ~ p_map_cant_inv(
    df = .x,
    var = "dcant_pos_bias",
    col = "bias",
    subtitle_text = paste("Estimate:", .x$estimate)
  ))
[[1]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[2]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[3]]

Version Author Date
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
88f7356 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c68ae3a jens-daniel-mueller 2021-07-13

[[4]]

Version Author Date
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
deb425c jens-daniel-mueller 2021-08-24
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
88f7356 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c68ae3a jens-daniel-mueller 2021-07-13
axis_limit <- dcant_mod_inv_bias %>% 
  pivot_longer(emlr_anom:nss_projection,
               names_to = "estimate",
               values_to = "bias") %>% 
  pull(bias) %>% 
  abs() %>% 
  max()

dcant_mod_inv_bias %>%
  ggplot(aes(emlr_bias, emlr_anom)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap( ~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
dcant_mod_inv_bias %>% 
  ggplot(aes(emlr_anom, nss_mod_truth)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap(~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
dcant_mod_inv_bias %>% 
  ggplot(aes(nss_mod_truth, nss_projection)) +
  geom_hline(yintercept = 0) +
  geom_vline(xintercept = 0) +
  geom_bin2d(binwidth = 1) +
  scale_fill_viridis_c(trans = "log10") +
  coord_equal(xlim = c(-axis_limit, axis_limit),
              ylim = c(-axis_limit, axis_limit)) +
  labs(title = "Bias vs bias plots") +
  facet_wrap(~ basin_AIP)

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

5.4 Budgets

dcant_mod_budget <- dcant_mod_inv %>% 
  mutate(surface_area = earth_surf(lat, lon),
         # dcant_inv_grid = dcant_inv*surface_area,
         dcant_pos_grid = dcant_pos*surface_area) %>% 
  group_by(basin_AIP, estimate, inv_depth) %>% 
  summarise(#dcant = sum(dcant_inv_grid)*12*1e-15,
            #dcant = round(dcant,1),
            dcant_pos = sum(dcant_pos_grid)*12*1e-15,
            dcant_pos = round(dcant_pos,1)) %>% 
  ungroup()

duration <- sort(tref$median_year)[2] - sort(tref$median_year)[1]

dcant_mod_budget %>%  
  group_by(estimate, inv_depth) %>% 
  summarise(#dcant = sum(dcant),
            #dcant = round(dcant,1),
            dcant_pos = sum(dcant_pos),
            dcant_pos = round(dcant_pos,1),
            #dcant_rate = dcant / duration,
            dcant_pos_rate = dcant_pos / duration) %>% 
  ungroup()
# A tibble: 20 × 4
   estimate     inv_depth dcant_pos dcant_pos_rate
   <chr>            <dbl>     <dbl>          <dbl>
 1 emlr               100       3.6           0.36
 2 emlr               500      12             1.2 
 3 emlr              1000      15.5           1.55
 4 emlr              3000      18             1.8 
 5 emlr             10000      18.8           1.88
 6 mod_truth_cc       100       3.4           0.34
 7 mod_truth_cc       500      11.6           1.16
 8 mod_truth_cc      1000      15.7           1.57
 9 mod_truth_cc      3000      17.5           1.75
10 mod_truth_cc     10000      17.5           1.75
11 mod_truth_vc       100       3.3           0.33
12 mod_truth_vc       500      11.8           1.18
13 mod_truth_vc      1000      15.7           1.57
14 mod_truth_vc      3000      17.8           1.78
15 mod_truth_vc     10000      18             1.8 
16 ss_tss             100       3.2           0.32
17 ss_tss             500      11.2           1.12
18 ss_tss            1000      14.6           1.46
19 ss_tss            3000      16.1           1.61
20 ss_tss           10000      16.3           1.63

5.4.1 Standard depth

Results integrated over the upper 3000 m

dcant_mod_budget %>%
  filter(inv_depth == params_global$inventory_depth_standard) %>% 
  ggplot(aes(estimate, dcant_pos, fill = basin_AIP)) +
  scale_fill_brewer(palette = "Dark2") +
  geom_col()

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
dcant_mod_budget %>%
  filter(inv_depth == params_global$inventory_depth_standard) %>% 
  gt(rowname_col = "basin_AIP",
     groupname_col = c("estimate", "inv_depth"),
     row_group.sep = " | Depth: ") %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
dcant_pos
emlr | Depth: 3000
Atlantic 5.3
Indian 3.2
Pacific 9.5
total 18.00
mod_truth_cc | Depth: 3000
Atlantic 4.4
Indian 4.7
Pacific 8.4
total 17.50
mod_truth_vc | Depth: 3000
Atlantic 4.5
Indian 4.8
Pacific 8.5
total 17.80
ss_tss | Depth: 3000
Atlantic 4.1
Indian 4.2
Pacific 7.8
total 16.10

5.4.2 Other depths

Results integrated over the upper 100, 500, 1000, 3000, 10^{4} m

dcant_mod_budget %>%
  filter(inv_depth != params_global$inventory_depth_standard) %>% 
  ggplot(aes(estimate, dcant_pos, fill = basin_AIP)) +
  scale_fill_brewer(palette = "Dark2") +
  geom_col() +
  facet_grid(inv_depth ~ ., scales = "free_y")

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
dcant_mod_budget %>%
  filter(inv_depth != params_global$inventory_depth_standard) %>% 
  gt(rowname_col = "basin_AIP",
     groupname_col = c("estimate", "inv_depth"),
     row_group.sep = " | Depth: ") %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
dcant_pos
emlr | Depth: 100
Atlantic 1.1
Indian 0.7
Pacific 1.8
total 3.60
emlr | Depth: 500
Atlantic 3.5
Indian 2.2
Pacific 6.3
total 12.00
emlr | Depth: 1000
Atlantic 4.6
Indian 3.0
Pacific 7.9
total 15.50
emlr | Depth: 10000
Atlantic 5.5
Indian 3.3
Pacific 10.0
total 18.80
mod_truth_cc | Depth: 100
Atlantic 0.9
Indian 0.7
Pacific 1.8
total 3.40
mod_truth_cc | Depth: 500
Atlantic 2.9
Indian 2.7
Pacific 6.0
total 11.60
mod_truth_cc | Depth: 1000
Atlantic 3.9
Indian 4.1
Pacific 7.7
total 15.70
mod_truth_cc | Depth: 10000
Atlantic 4.4
Indian 4.7
Pacific 8.4
total 17.50
mod_truth_vc | Depth: 100
Atlantic 0.9
Indian 0.7
Pacific 1.7
total 3.30
mod_truth_vc | Depth: 500
Atlantic 3.0
Indian 2.8
Pacific 6.0
total 11.80
mod_truth_vc | Depth: 1000
Atlantic 3.9
Indian 4.1
Pacific 7.7
total 15.70
mod_truth_vc | Depth: 10000
Atlantic 4.5
Indian 4.9
Pacific 8.6
total 18.00
ss_tss | Depth: 100
Atlantic 0.8
Indian 0.7
Pacific 1.7
total 3.20
ss_tss | Depth: 500
Atlantic 2.8
Indian 2.6
Pacific 5.8
total 11.20
ss_tss | Depth: 1000
Atlantic 3.6
Indian 3.7
Pacific 7.3
total 14.60
ss_tss | Depth: 10000
Atlantic 4.1
Indian 4.3
Pacific 7.9
total 16.30

6 Observations

6.1 Zonal mean sections

6.1.1 delta dcant

# calculate zonal mean section
dcant_obs_zonal <- dcant_obs_3d %>%
  group_by(estimate) %>%
  nest() %>%
  mutate(section = map(.x = data, ~m_zonal_mean_sd(.x))) %>%
  select(-data) %>%
  unnest(section)

dcant_obs_zonal <- dcant_obs_zonal %>% 
  rename(dcant_pos = dcant_pos_mean)

6.1.2 Bias

dcant_obs_zonal_bias <- dcant_obs_zonal %>%
  select(-dcant_pos_sd) %>% 
  pivot_wider(names_from = estimate,
              values_from = dcant_pos) %>% 
  drop_na()


dcant_obs_zonal_bias <- dcant_obs_zonal_bias %>% 
  mutate(emlr_anom = emlr - ss_tss) %>% 
  select(-c(emlr:ss_tss))

dcant_obs_zonal_bias_long <- dcant_obs_zonal_bias %>% 
  pivot_longer(emlr_anom,
               names_to = "estimate",
               values_to = "dcant_pos_bias")
  


dcant_obs_zonal_bias_long %>% 
  group_by(basin_AIP, estimate) %>%
  group_split() %>% 
  # head(1) %>%
  map( ~ p_section_zonal(
    df = .x,
    var = "dcant_pos_bias",
    col = "divergent",
    breaks = params_global$breaks_dcant_offset,
    plot_slabs = "n",
    subtitle_text = paste("Basin:", .x$basin_AIP, "| Estimate:", .x$estimate)
  ))
[[1]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

[[2]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
7bc0dd5 jens-daniel-mueller 2021-08-04

[[3]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
7bc0dd5 jens-daniel-mueller 2021-08-04

6.2 Inventories

6.2.1 delta dcant

# calculate zonal mean section
dcant_obs_inv <- dcant_obs_3d %>%
  rename(dcant_pos = dcant_pos) %>% 
  mutate(dcant = 0) %>% 
  group_by(estimate) %>%
  nest() %>%
  mutate(inv = map(.x = data, ~m_dcant_inv(.x))) %>%
  select(-data) %>%
  unnest(inv) %>% 
  select(-c(dcant))

6.2.2 Bias

dcant_obs_inv_bias <- dcant_obs_inv %>%
  filter(inv_depth == params_global$inventory_depth_standard) %>% 
  select(-inv_depth) %>% 
  pivot_wider(names_from = estimate,
              values_from = dcant_pos) %>% 
  drop_na()


dcant_obs_inv_bias <- dcant_obs_inv_bias %>% 
  mutate(emlr_anom = emlr - ss_tss) %>% 
  select(-c(emlr:ss_tss))


dcant_obs_inv_bias_long <- dcant_obs_inv_bias %>% 
  pivot_longer(emlr_anom,
               values_to = "dcant_pos_bias",
               names_to = "estimate")
  


dcant_obs_inv_bias_long %>% 
  group_by(estimate) %>%
  group_split() %>% 
  # head(1) %>%
  map( ~ p_map_cant_inv(
    df = .x,
    var = "dcant_pos_bias",
    col = "bias",
    subtitle_text = paste("Estimate:", .x$estimate)
  ))
[[1]]

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13

6.3 Budgets

dcant_obs_budget <- dcant_obs_inv %>% 
  mutate(surface_area = earth_surf(lat, lon),
         # dcant_inv_grid = dcant_inv*surface_area,
         dcant_pos_grid = dcant_pos*surface_area) %>% 
  group_by(basin_AIP, estimate, inv_depth) %>% 
  summarise(#dcant = sum(dcant_inv_grid)*12*1e-15,
            #dcant = round(dcant,1),
            dcant_pos = sum(dcant_pos_grid)*12*1e-15,
            dcant_pos = round(dcant_pos,1)) %>% 
  ungroup()

duration <- sort(tref$median_year)[2] - sort(tref$median_year)[1]

dcant_obs_budget %>%  
  group_by(estimate, inv_depth) %>% 
  summarise(#dcant = sum(dcant),
            #dcant = round(dcant,1),
            dcant_pos = sum(dcant_pos),
            dcant_pos = round(dcant_pos,1),
            #dcant_rate = dcant / duration,
            dcant_pos_rate = dcant_pos / duration) %>% 
  ungroup()
# A tibble: 10 × 4
   estimate inv_depth dcant_pos dcant_pos_rate
   <chr>        <dbl>     <dbl>          <dbl>
 1 emlr           100       4             0.4 
 2 emlr           500      13.6           1.36
 3 emlr          1000      18.9           1.89
 4 emlr          3000      30.1           3.01
 5 emlr         10000      36.1           3.61
 6 ss_tss         100       3.5           0.35
 7 ss_tss         500      12.9           1.29
 8 ss_tss        1000      18.1           1.81
 9 ss_tss        3000      21.4           2.14
10 ss_tss       10000      21.7           2.17

6.3.1 Standard depth

Results integrated over the upper 3000 m

dcant_obs_budget %>%
  filter(inv_depth == params_global$inventory_depth_standard) %>% 
  ggplot(aes(estimate, dcant_pos, fill = basin_AIP)) +
  scale_fill_brewer(palette = "Dark2") +
  geom_col()

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
dcant_obs_budget %>%
  filter(inv_depth == params_global$inventory_depth_standard) %>% 
  gt(rowname_col = "basin_AIP",
     groupname_col = c("estimate", "inv_depth"),
     row_group.sep = " | Depth: ") %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
dcant_pos
emlr | Depth: 3000
Atlantic 11.7
Indian 7.8
Pacific 10.6
total 30.10
ss_tss | Depth: 3000
Atlantic 8.1
Indian 4.8
Pacific 8.5
total 21.40

6.3.2 Other depths

Results integrated over the upper 100, 500, 1000, 3000, 10^{4} m

dcant_obs_budget %>%
  filter(inv_depth != params_global$inventory_depth_standard) %>% 
  ggplot(aes(estimate, dcant_pos, fill = basin_AIP)) +
  scale_fill_brewer(palette = "Dark2") +
  geom_col() +
  facet_grid(inv_depth ~ ., scales = "free_y")

Version Author Date
98d9e33 jens-daniel-mueller 2021-11-11
d3cb92d jens-daniel-mueller 2021-11-08
3879a6d jens-daniel-mueller 2021-11-08
0e032dc jens-daniel-mueller 2021-11-04
abcd28f jens-daniel-mueller 2021-11-02
290c8fc jens-daniel-mueller 2021-11-02
e02acc9 jens-daniel-mueller 2021-11-01
58da811 jens-daniel-mueller 2021-11-01
2781a97 jens-daniel-mueller 2021-10-29
973192c jens-daniel-mueller 2021-10-28
f792cd8 jens-daniel-mueller 2021-10-27
22750db jens-daniel-mueller 2021-10-27
3992e41 jens-daniel-mueller 2021-10-27
f9280c7 jens-daniel-mueller 2021-10-21
bfe90e6 jens-daniel-mueller 2021-10-21
78a956e jens-daniel-mueller 2021-10-21
1ba14cf jens-daniel-mueller 2021-10-20
dc0723b jens-daniel-mueller 2021-10-19
722483f jens-daniel-mueller 2021-10-18
cd0ebd8 jens-daniel-mueller 2021-10-18
aedafcb jens-daniel-mueller 2021-10-17
046c556 jens-daniel-mueller 2021-10-07
ac73c23 jens-daniel-mueller 2021-10-07
0bd4ea2 jens-daniel-mueller 2021-10-07
6969c18 jens-daniel-mueller 2021-10-07
06afaca jens-daniel-mueller 2021-09-23
410dd46 jens-daniel-mueller 2021-09-22
5a36231 jens-daniel-mueller 2021-09-22
c0bbbe4 jens-daniel-mueller 2021-09-22
9eac969 jens-daniel-mueller 2021-08-30
3a584b7 jens-daniel-mueller 2021-08-30
6e17588 jens-daniel-mueller 2021-08-24
deb425c jens-daniel-mueller 2021-08-24
534854b jens-daniel-mueller 2021-08-19
77f4ba7 jens-daniel-mueller 2021-08-19
ece96df jens-daniel-mueller 2021-08-19
27c99b8 jens-daniel-mueller 2021-08-19
a03f2f0 jens-daniel-mueller 2021-08-18
9335b31 jens-daniel-mueller 2021-08-10
9943b45 jens-daniel-mueller 2021-08-10
0b00a2b jens-daniel-mueller 2021-08-09
755c6b1 jens-daniel-mueller 2021-08-09
cd8e0d5 jens-daniel-mueller 2021-08-06
15773a0 jens-daniel-mueller 2021-08-06
da61d1a jens-daniel-mueller 2021-08-06
340d731 jens-daniel-mueller 2021-08-06
71546e4 jens-daniel-mueller 2021-08-06
29444a1 jens-daniel-mueller 2021-08-05
42e80c0 jens-daniel-mueller 2021-08-04
48f6eed jens-daniel-mueller 2021-08-04
81a46a4 jens-daniel-mueller 2021-08-03
b88c61b jens-daniel-mueller 2021-08-03
a53656d jens-daniel-mueller 2021-08-03
88f7356 jens-daniel-mueller 2021-08-02
d759279 jens-daniel-mueller 2021-08-02
127b801 jens-daniel-mueller 2021-07-24
912d90e jens-daniel-mueller 2021-07-23
2477316 jens-daniel-mueller 2021-07-23
c9ccc00 jens-daniel-mueller 2021-07-22
426b2df jens-daniel-mueller 2021-07-21
c68ae3a jens-daniel-mueller 2021-07-13
dcant_obs_budget %>%
  filter(inv_depth != params_global$inventory_depth_standard) %>% 
  gt(rowname_col = "basin_AIP",
     groupname_col = c("estimate", "inv_depth"),
     row_group.sep = " | Depth: ") %>% 
  summary_rows(
    groups = TRUE,
    fns = list(total = "sum")
  )
dcant_pos
emlr | Depth: 100
Atlantic 1.1
Indian 0.8
Pacific 2.1
total 4.00
emlr | Depth: 500
Atlantic 3.7
Indian 3.0
Pacific 6.9
total 13.60
emlr | Depth: 1000
Atlantic 5.6
Indian 4.6
Pacific 8.7
total 18.90
emlr | Depth: 10000
Atlantic 15.1
Indian 8.6
Pacific 12.4
total 36.10
ss_tss | Depth: 100
Atlantic 1.0
Indian 0.8
Pacific 1.7
total 3.50
ss_tss | Depth: 500
Atlantic 3.8
Indian 2.9
Pacific 6.2
total 12.90
ss_tss | Depth: 1000
Atlantic 5.8
Indian 4.4
Pacific 7.9
total 18.10
ss_tss | Depth: 10000
Atlantic 8.2
Indian 5.0
Pacific 8.5
total 21.70

7 Write files

dcant_mod_3d_bias %>%
  write_csv(paste0(path_version_data,
                   "anom_dcant_mod_3d_bias.csv"))

dcant_mod_zonal_bias %>%
  write_csv(paste0(path_version_data,
                   "anom_dcant_mod_zonal_bias.csv"))

dcant_mod_inv_bias %>%
  write_csv(paste0(path_version_data,
                   "anom_dcant_mod_inv_bias.csv"))

dcant_mod_budget %>%
  write_csv(paste0(path_version_data,
                   "anom_dcant_mod_budget.csv"))


dcant_obs_zonal_bias %>%
  write_csv(paste0(path_version_data,
                   "anom_dcant_obs_zonal_bias.csv"))

dcant_obs_inv_bias %>%
  write_csv(paste0(path_version_data,
                   "anom_dcant_obs_inv_bias.csv"))

dcant_obs_budget %>%
  write_csv(paste0(path_version_data,
                   "anom_dcant_obs_budget.csv"))

sessionInfo()
R version 4.0.3 (2020-10-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: openSUSE Leap 15.2

Matrix products: default
BLAS:   /usr/local/R-4.0.3/lib64/R/lib/libRblas.so
LAPACK: /usr/local/R-4.0.3/lib64/R/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] gt_0.2.2        marelac_2.1.10  shape_1.4.5     ggforce_0.3.3  
 [5] metR_0.9.0      scico_1.2.0     patchwork_1.1.1 collapse_1.5.0 
 [9] forcats_0.5.0   stringr_1.4.0   dplyr_1.0.5     purrr_0.3.4    
[13] readr_1.4.0     tidyr_1.1.3     tibble_3.1.3    ggplot2_3.3.5  
[17] tidyverse_1.3.0 workflowr_1.6.2

loaded via a namespace (and not attached):
 [1] fs_1.5.0                 gsw_1.0-5                lubridate_1.7.9         
 [4] RColorBrewer_1.1-2       httr_1.4.2               rprojroot_2.0.2         
 [7] tools_4.0.3              backports_1.1.10         bslib_0.2.5.1           
[10] utf8_1.1.4               R6_2.5.0                 DBI_1.1.0               
[13] colorspace_2.0-2         withr_2.3.0              tidyselect_1.1.0        
[16] compiler_4.0.3           git2r_0.27.1             cli_3.0.1               
[19] rvest_0.3.6              xml2_1.3.2               isoband_0.2.2           
[22] labeling_0.4.2           sass_0.4.0               scales_1.1.1            
[25] checkmate_2.0.0          digest_0.6.27            rmarkdown_2.10          
[28] oce_1.2-0                pkgconfig_2.0.3          htmltools_0.5.1.1       
[31] highr_0.8                dbplyr_1.4.4             rlang_0.4.11            
[34] readxl_1.3.1             rstudioapi_0.13          jquerylib_0.1.4         
[37] generics_0.1.0           farver_2.0.3             jsonlite_1.7.1          
[40] magrittr_1.5             Matrix_1.2-18            Rcpp_1.0.5              
[43] munsell_0.5.0            fansi_0.4.1              lifecycle_1.0.0         
[46] stringi_1.5.3            whisker_0.4              yaml_2.2.1              
[49] MASS_7.3-53              grid_4.0.3               blob_1.2.1              
[52] parallel_4.0.3           promises_1.1.1           crayon_1.3.4            
[55] lattice_0.20-41          haven_2.3.1              hms_0.5.3               
[58] seacarb_3.2.14           knitr_1.33               pillar_1.6.2            
[61] reprex_0.3.0             glue_1.4.2               evaluate_0.14           
[64] RcppArmadillo_0.10.1.2.0 data.table_1.14.0        modelr_0.1.8            
[67] vctrs_0.3.8              tweenr_1.0.2             httpuv_1.5.4            
[70] testthat_2.3.2           cellranger_1.1.0         gtable_0.3.0            
[73] polyclip_1.10-0          assertthat_0.2.1         xfun_0.25               
[76] broom_0.7.9              RcppEigen_0.3.3.7.0      later_1.2.0             
[79] viridisLite_0.3.0        ellipsis_0.3.2           here_0.1