Last updated: 2020-09-01

Checks: 7 0

Knit directory: Cant_eMLR/

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 8f3ce45. 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:    .Rproj.user/
    Ignored:    data/GLODAPv1_1/
    Ignored:    data/GLODAPv2_2016b_MappedClimatologies/
    Ignored:    data/GLODAPv2_2020/
    Ignored:    data/Gruber_2019/
    Ignored:    data/WOCE/
    Ignored:    data/World_Ocean_Atlas_2013_Clement/
    Ignored:    data/World_Ocean_Atlas_2018/
    Ignored:    data/eMLR/
    Ignored:    data/mapping/
    Ignored:    data/pCO2_atmosphere/
    Ignored:    dump/

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/eMLR_model_fitting.Rmd) and HTML (docs/eMLR_model_fitting.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
Rmd 8f3ce45 jens-daniel-mueller 2020-09-01 rebuild without PO4 star selection, oxygen only
html 13a76d5 jens-daniel-mueller 2020-08-28 Build site.
Rmd 2e6a4ca jens-daniel-mueller 2020-08-28 XXX
html 27404de jens-daniel-mueller 2020-08-27 Build site.
html b610bf7 jens-daniel-mueller 2020-08-27 Build site.
Rmd 96ab537 jens-daniel-mueller 2020-08-27 revised
html b6d0e6a jens-daniel-mueller 2020-08-27 Build site.
html f40e48b jens-daniel-mueller 2020-08-26 Build site.
html ec20f40 jens-daniel-mueller 2020-08-24 Build site.
Rmd a804955 jens-daniel-mueller 2020-08-24 split mapping into 2 rmds, po4star selection in parameters, use po4star nitrate
html 5ffe187 jens-daniel-mueller 2020-08-20 Build site.
html 1064ef8 jens-daniel-mueller 2020-08-19 Build site.
Rmd 94f9375 jens-daniel-mueller 2020-08-19 split emlr into data preparation, assumption testing, and model fitting
html 8484d53 jens-daniel-mueller 2020-08-19 Build site.
Rmd eec0022 jens-daniel-mueller 2020-08-19 split emlr into data preparation, assumption testing, and model fitting

library(tidyverse)
library(lubridate)
library(patchwork)
library(GGally)
library(olsrr)
library(knitr)
library(kableExtra)
library(broom)
library(corrr)

1 Required data

Required are:

  • cleaned and prepared GLODAPv2.2020 file
GLODAP <-
  read_csv(
    here::here(
      "data/GLODAPv2_2020/_summarized_data_files",
      "GLODAP_MLR_fitting_ready.csv"
    )
  )

2 Model fitting

Individual linear regression models were fitted for C* as a function of each predictor combination. Fitting was performed separately within each basin, era, and slab. The root mean squared error (RMSE) was calculated for each fitted model.

for (i_basin in unique(GLODAP$basin)) {
  for (i_era in unique(GLODAP$era)) {
    
    # i_basin <- unique(GLODAP$basin)[1]
    # i_era   <- unique(GLODAP$era)[1]
    
    print(i_basin) 
    print(i_era  )
    
    GLODAP_basin_era <- GLODAP %>%
      filter(basin == i_basin, era == i_era)

    for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {

      # i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
      print(i_gamma_slab)
      
      GLODAP_basin_era_slab <- GLODAP_basin_era %>% 
        filter(gamma_slab == i_gamma_slab)
      
      lm_full <- lm(Cstar ~ 
                sal  + 
                tem  + 
                aou  + 
                oxygen  + 
                silicate  + 
                phosphate  + 
                phosphate_star, 
              data = GLODAP_basin_era_slab)
      
      lm_all <- ols_step_all_possible(lm_full)
      
      temp <- as_tibble(lm_all$result)
      
      temp <- temp %>% 
        mutate(basin = i_basin,
               era = i_era,
               gamma_slab = i_gamma_slab)
      
      if (exists("all_lm_stats")) {
        all_lm_stats <- bind_rows(all_lm_stats, temp)
      }
      
      if (!exists("all_lm_stats")) {
        all_lm_stats <- temp
      }
      
      rm(lm_full, lm_all, temp)
      
    }
    
  }
  
}

all_lm_stats %>%
  write_csv(here::here("data/eMLR",
                       "all_lm_stats.csv"))

rm(i_gamma_slab, i_era, i_basin,
   GLODAP_basin_era, GLODAP_basin_era_slab,
   all_lm_stats)

RMSE was plotted against the number of predictors (limited to 2 - 5).

all_lm_stats <- read_csv(here::here("data/eMLR",
                                    "all_lm_stats.csv"))

all_lm_stats <- all_lm_stats %>% 
  filter(n >= parameters$predictors_min,
         n <= parameters$predictors_max)

all_lm_stats %>% 
  ggplot(aes(n, rmse, col = basin)) +
  geom_hline(yintercept = 10) + 
  geom_point(shape = 21) +
  facet_grid(gamma_slab~era) +
  scale_color_brewer(palette = "Set1")

2.1 RMSE alternatives

AIC is an alternative criterion to RMSE to judge model quality, but not (yet) taken into account.

all_lm_stats %>% 
  ggplot(aes(rmse, aic, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(era~basin)

2.2 Predictor selection

Within each basin and slab, the 10 linear regression models with lowest summed RMSE across all eras were selected.

all_lm_stats <- all_lm_stats %>%
  select(basin, era, gamma_slab, predictors, rmse)

all_lm_stats_wide <- all_lm_stats %>%
  pivot_wider(names_from = era,
              values_from = rmse,
              names_prefix = "rmse_")

all_lm_stats_wide <- all_lm_stats_wide %>%
  mutate(
    rmse_sum_JGOFS_GO = rmse_JGOFS_WOCE + rmse_GO_SHIP,
    rmse_sum_GO_new = rmse_GO_SHIP + rmse_new_era
  )

all_lm_stats <- all_lm_stats_wide %>%
  pivot_longer(
    c(rmse_sum_JGOFS_GO, rmse_sum_GO_new),
    names_to = "eras",
    values_to = "rmse_sum",
    names_prefix = "rmse_sum_"
  )

lm_best <- all_lm_stats %>%
  group_by(basin, gamma_slab, eras) %>%
  slice_min(order_by = rmse_sum,
            with_ties = FALSE,
            n = parameters$MLR_number) %>%
  ungroup() %>% 
  arrange(basin, gamma_slab, eras, predictors)

kable(lm_best) %>%
  add_header_above() %>%
  kable_styling() %>%
  scroll_box(width = "100%", height = "400px")
basin gamma_slab predictors rmse_JGOFS_WOCE rmse_GO_SHIP rmse_new_era eras rmse_sum
Atlantic (-Inf,26] aou oxygen phosphate phosphate_star 8.223378 6.151569 2.597069 GO_new 8.748638
Atlantic (-Inf,26] aou oxygen silicate phosphate phosphate_star 7.551415 6.151535 2.585332 GO_new 8.736868
Atlantic (-Inf,26] sal aou oxygen phosphate phosphate_star 7.713991 5.768405 2.539469 GO_new 8.307874
Atlantic (-Inf,26] sal tem aou phosphate phosphate_star 7.585762 5.807751 2.536513 GO_new 8.344264
Atlantic (-Inf,26] sal tem oxygen phosphate phosphate_star 7.547678 5.837633 2.536054 GO_new 8.373687
Atlantic (-Inf,26] tem aou oxygen phosphate phosphate_star 7.150296 6.028173 2.567006 GO_new 8.595179
Atlantic (-Inf,26] tem aou phosphate phosphate_star 8.053366 6.103861 2.593184 GO_new 8.697045
Atlantic (-Inf,26] tem aou silicate phosphate phosphate_star 7.415798 6.103595 2.582598 GO_new 8.686193
Atlantic (-Inf,26] tem oxygen phosphate phosphate_star 7.944696 6.058427 2.586913 GO_new 8.645340
Atlantic (-Inf,26] tem oxygen silicate phosphate phosphate_star 7.358462 6.055636 2.578958 GO_new 8.634594
Atlantic (-Inf,26] aou oxygen silicate phosphate phosphate_star 7.551415 6.151535 2.585332 JGOFS_GO 13.702951
Atlantic (-Inf,26] sal aou oxygen phosphate phosphate_star 7.713991 5.768405 2.539469 JGOFS_GO 13.482396
Atlantic (-Inf,26] sal tem aou oxygen phosphate 7.173663 6.639387 2.699587 JGOFS_GO 13.813050
Atlantic (-Inf,26] sal tem aou phosphate phosphate_star 7.585762 5.807751 2.536513 JGOFS_GO 13.393513
Atlantic (-Inf,26] sal tem oxygen phosphate phosphate_star 7.547678 5.837633 2.536054 JGOFS_GO 13.385311
Atlantic (-Inf,26] sal tem oxygen silicate phosphate 7.248967 6.745194 2.784042 JGOFS_GO 13.994160
Atlantic (-Inf,26] tem aou oxygen phosphate phosphate_star 7.150296 6.028173 2.567006 JGOFS_GO 13.178469
Atlantic (-Inf,26] tem aou silicate phosphate phosphate_star 7.415798 6.103595 2.582598 JGOFS_GO 13.519393
Atlantic (-Inf,26] tem oxygen phosphate phosphate_star 7.944696 6.058427 2.586913 JGOFS_GO 14.003123
Atlantic (-Inf,26] tem oxygen silicate phosphate phosphate_star 7.358462 6.055636 2.578958 JGOFS_GO 13.414099
Atlantic (26,26.5] aou oxygen phosphate phosphate_star 5.479554 7.263998 3.094275 GO_new 10.358274
Atlantic (26,26.5] aou oxygen silicate phosphate phosphate_star 5.391757 7.262779 3.059189 GO_new 10.321968
Atlantic (26,26.5] sal aou oxygen phosphate phosphate_star 5.456763 7.247045 3.078265 GO_new 10.325310
Atlantic (26,26.5] sal tem aou phosphate phosphate_star 5.502307 7.278155 3.085180 GO_new 10.363335
Atlantic (26,26.5] sal tem oxygen phosphate phosphate_star 5.558040 7.317006 3.092409 GO_new 10.409414
Atlantic (26,26.5] tem aou oxygen phosphate phosphate_star 5.452186 7.236607 3.080579 GO_new 10.317186
Atlantic (26,26.5] tem aou oxygen silicate phosphate 5.363554 7.360474 3.065512 GO_new 10.425986
Atlantic (26,26.5] tem aou phosphate phosphate_star 5.509470 7.284165 3.103494 GO_new 10.387658
Atlantic (26,26.5] tem aou silicate phosphate phosphate_star 5.429777 7.282652 3.070020 GO_new 10.352672
Atlantic (26,26.5] tem oxygen silicate phosphate phosphate_star 5.490007 7.317880 3.082355 GO_new 10.400236
Atlantic (26,26.5] aou oxygen phosphate phosphate_star 5.479554 7.263998 3.094275 JGOFS_GO 12.743552
Atlantic (26,26.5] aou oxygen silicate phosphate phosphate_star 5.391757 7.262779 3.059189 JGOFS_GO 12.654537
Atlantic (26,26.5] sal aou oxygen phosphate phosphate_star 5.456763 7.247045 3.078265 JGOFS_GO 12.703808
Atlantic (26,26.5] sal tem aou oxygen phosphate 5.477425 7.346736 3.118023 JGOFS_GO 12.824162
Atlantic (26,26.5] sal tem aou phosphate phosphate_star 5.502307 7.278155 3.085180 JGOFS_GO 12.780462
Atlantic (26,26.5] tem aou oxygen phosphate phosphate_star 5.452186 7.236607 3.080579 JGOFS_GO 12.688793
Atlantic (26,26.5] tem aou oxygen silicate phosphate 5.363554 7.360474 3.065512 JGOFS_GO 12.724027
Atlantic (26,26.5] tem aou phosphate phosphate_star 5.509470 7.284165 3.103494 JGOFS_GO 12.793635
Atlantic (26,26.5] tem aou silicate phosphate phosphate_star 5.429777 7.282652 3.070020 JGOFS_GO 12.712429
Atlantic (26,26.5] tem oxygen silicate phosphate phosphate_star 5.490007 7.317880 3.082355 JGOFS_GO 12.807887
Atlantic (26.5,26.75] sal aou oxygen silicate phosphate 5.121400 6.278369 3.156785 GO_new 9.435155
Atlantic (26.5,26.75] sal tem aou oxygen phosphate 5.318827 6.182233 3.188366 GO_new 9.370599
Atlantic (26.5,26.75] sal tem aou phosphate 5.348852 6.294451 3.190597 GO_new 9.485048
Atlantic (26.5,26.75] sal tem aou phosphate phosphate_star 5.340770 6.241134 3.189793 GO_new 9.430926
Atlantic (26.5,26.75] sal tem aou silicate phosphate 5.080928 6.189493 3.157471 GO_new 9.346963
Atlantic (26.5,26.75] sal tem oxygen phosphate 5.332037 6.268202 3.188583 GO_new 9.456786
Atlantic (26.5,26.75] sal tem oxygen phosphate phosphate_star 5.324850 6.211578 3.188105 GO_new 9.399684
Atlantic (26.5,26.75] sal tem oxygen silicate phosphate 5.073639 6.161326 3.152137 GO_new 9.313462
Atlantic (26.5,26.75] tem aou oxygen silicate phosphate 5.153695 6.171505 3.261512 GO_new 9.433017
Atlantic (26.5,26.75] tem oxygen silicate phosphate phosphate_star 5.141157 6.276629 3.197287 GO_new 9.473917
Atlantic (26.5,26.75] sal aou oxygen silicate phosphate 5.121400 6.278369 3.156785 JGOFS_GO 11.399769
Atlantic (26.5,26.75] sal tem aou oxygen phosphate 5.318827 6.182233 3.188366 JGOFS_GO 11.501060
Atlantic (26.5,26.75] sal tem aou silicate phosphate 5.080928 6.189493 3.157471 JGOFS_GO 11.270421
Atlantic (26.5,26.75] sal tem oxygen phosphate phosphate_star 5.324850 6.211578 3.188105 JGOFS_GO 11.536428
Atlantic (26.5,26.75] sal tem oxygen silicate phosphate 5.073639 6.161326 3.152137 JGOFS_GO 11.234964
Atlantic (26.5,26.75] tem aou oxygen silicate phosphate 5.153695 6.171505 3.261512 JGOFS_GO 11.325200
Atlantic (26.5,26.75] tem aou silicate phosphate 5.173069 6.316120 3.273435 JGOFS_GO 11.489188
Atlantic (26.5,26.75] tem aou silicate phosphate phosphate_star 5.154882 6.309885 3.201271 JGOFS_GO 11.464768
Atlantic (26.5,26.75] tem oxygen silicate phosphate 5.158997 6.281083 3.267692 JGOFS_GO 11.440080
Atlantic (26.5,26.75] tem oxygen silicate phosphate phosphate_star 5.141157 6.276629 3.197287 JGOFS_GO 11.417786
Atlantic (26.75,27] aou oxygen silicate phosphate 4.896401 5.405821 3.681006 GO_new 9.086827
Atlantic (26.75,27] aou oxygen silicate phosphate phosphate_star 4.835819 5.377290 3.675956 GO_new 9.053246
Atlantic (26.75,27] sal aou oxygen silicate phosphate 4.810239 5.344351 3.679790 GO_new 9.024141
Atlantic (26.75,27] sal tem aou silicate phosphate 4.708472 5.279945 3.666486 GO_new 8.946431
Atlantic (26.75,27] sal tem oxygen silicate phosphate 4.689137 5.255316 3.664798 GO_new 8.920113
Atlantic (26.75,27] tem aou oxygen silicate phosphate 4.802227 5.348802 3.665657 GO_new 9.014459
Atlantic (26.75,27] tem aou silicate phosphate 4.826200 5.373558 3.667384 GO_new 9.040943
Atlantic (26.75,27] tem aou silicate phosphate phosphate_star 4.767595 5.343790 3.662341 GO_new 9.006132
Atlantic (26.75,27] tem oxygen silicate phosphate 4.803963 5.356851 3.666006 GO_new 9.022858
Atlantic (26.75,27] tem oxygen silicate phosphate phosphate_star 4.739924 5.321954 3.660579 GO_new 8.982533
Atlantic (26.75,27] sal aou oxygen silicate phosphate 4.810239 5.344351 3.679790 JGOFS_GO 10.154590
Atlantic (26.75,27] sal tem aou oxygen phosphate 4.747838 5.428829 3.705508 JGOFS_GO 10.176667
Atlantic (26.75,27] sal tem aou silicate phosphate 4.708472 5.279945 3.666486 JGOFS_GO 9.988417
Atlantic (26.75,27] sal tem oxygen phosphate phosphate_star 4.737133 5.414265 3.678545 JGOFS_GO 10.151398
Atlantic (26.75,27] sal tem oxygen silicate phosphate 4.689137 5.255316 3.664798 JGOFS_GO 9.944452
Atlantic (26.75,27] tem aou oxygen phosphate phosphate_star 4.761392 5.397500 3.708802 JGOFS_GO 10.158892
Atlantic (26.75,27] tem aou oxygen silicate phosphate 4.802227 5.348802 3.665657 JGOFS_GO 10.151029
Atlantic (26.75,27] tem aou silicate phosphate phosphate_star 4.767595 5.343790 3.662341 JGOFS_GO 10.111385
Atlantic (26.75,27] tem oxygen silicate phosphate 4.803963 5.356851 3.666006 JGOFS_GO 10.160814
Atlantic (26.75,27] tem oxygen silicate phosphate phosphate_star 4.739924 5.321954 3.660579 JGOFS_GO 10.061879
Atlantic (27,27.25] aou oxygen silicate phosphate phosphate_star 5.328621 5.669303 4.933733 GO_new 10.603037
Atlantic (27,27.25] sal tem aou phosphate phosphate_star 5.250498 5.696724 5.034078 GO_new 10.730802
Atlantic (27,27.25] sal tem aou silicate phosphate 5.078648 5.623468 5.051530 GO_new 10.674998
Atlantic (27,27.25] sal tem oxygen phosphate phosphate_star 5.196334 5.658172 5.054511 GO_new 10.712683
Atlantic (27,27.25] sal tem oxygen silicate phosphate 5.018688 5.589122 5.077798 GO_new 10.666920
Atlantic (27,27.25] tem aou oxygen phosphate phosphate_star 5.147498 5.619345 4.971314 GO_new 10.590659
Atlantic (27,27.25] tem aou oxygen silicate phosphate 5.132734 5.610171 5.007417 GO_new 10.617588
Atlantic (27,27.25] tem aou silicate phosphate phosphate_star 5.192439 5.582350 4.937753 GO_new 10.520102
Atlantic (27,27.25] tem oxygen phosphate phosphate_star 5.234406 5.665572 5.060251 GO_new 10.725823
Atlantic (27,27.25] tem oxygen silicate phosphate phosphate_star 5.142108 5.541220 4.952800 GO_new 10.494021
Atlantic (27,27.25] sal tem aou silicate phosphate 5.078648 5.623468 5.051530 JGOFS_GO 10.702117
Atlantic (27,27.25] sal tem oxygen phosphate phosphate_star 5.196334 5.658172 5.054511 JGOFS_GO 10.854506
Atlantic (27,27.25] sal tem oxygen silicate phosphate 5.018688 5.589122 5.077798 JGOFS_GO 10.607810
Atlantic (27,27.25] tem aou oxygen phosphate phosphate_star 5.147498 5.619345 4.971314 JGOFS_GO 10.766844
Atlantic (27,27.25] tem aou oxygen silicate phosphate 5.132734 5.610171 5.007417 JGOFS_GO 10.742905
Atlantic (27,27.25] tem aou silicate phosphate 5.199454 5.645126 5.100418 JGOFS_GO 10.844579
Atlantic (27,27.25] tem aou silicate phosphate phosphate_star 5.192439 5.582350 4.937753 JGOFS_GO 10.774788
Atlantic (27,27.25] tem oxygen phosphate phosphate_star 5.234406 5.665572 5.060251 JGOFS_GO 10.899978
Atlantic (27,27.25] tem oxygen silicate phosphate 5.156878 5.617843 5.133226 JGOFS_GO 10.774722
Atlantic (27,27.25] tem oxygen silicate phosphate phosphate_star 5.142108 5.541220 4.952800 JGOFS_GO 10.683329
Atlantic (27.25,27.5] aou oxygen phosphate phosphate_star 4.928107 4.547030 3.655045 GO_new 8.202075
Atlantic (27.25,27.5] aou oxygen silicate phosphate phosphate_star 4.915736 4.547015 3.648298 GO_new 8.195313
Atlantic (27.25,27.5] sal aou oxygen phosphate phosphate_star 4.902251 4.509664 3.633780 GO_new 8.143443
Atlantic (27.25,27.5] sal tem aou oxygen phosphate 4.685854 4.506963 3.523918 GO_new 8.030880
Atlantic (27.25,27.5] sal tem aou phosphate phosphate_star 4.823760 4.485051 3.663136 GO_new 8.148187
Atlantic (27.25,27.5] sal tem oxygen phosphate phosphate_star 4.777170 4.485453 3.700628 GO_new 8.186081
Atlantic (27.25,27.5] tem aou oxygen phosphate 4.776702 4.554342 3.528315 GO_new 8.082657
Atlantic (27.25,27.5] tem aou oxygen phosphate phosphate_star 4.772737 4.543040 3.421364 GO_new 7.964404
Atlantic (27.25,27.5] tem aou oxygen silicate phosphate 4.708352 4.554250 3.467433 GO_new 8.021683
Atlantic (27.25,27.5] tem aou silicate phosphate phosphate_star 4.856829 4.544790 3.677480 GO_new 8.222270
Atlantic (27.25,27.5] sal tem aou oxygen phosphate 4.685854 4.506963 3.523918 JGOFS_GO 9.192816
Atlantic (27.25,27.5] sal tem aou phosphate 4.825263 4.507063 3.906248 JGOFS_GO 9.332326
Atlantic (27.25,27.5] sal tem aou phosphate phosphate_star 4.823760 4.485051 3.663136 JGOFS_GO 9.308811
Atlantic (27.25,27.5] sal tem aou silicate phosphate 4.824954 4.494151 3.859025 JGOFS_GO 9.319105
Atlantic (27.25,27.5] sal tem oxygen phosphate 4.779201 4.514894 3.962830 JGOFS_GO 9.294095
Atlantic (27.25,27.5] sal tem oxygen phosphate phosphate_star 4.777170 4.485453 3.700628 JGOFS_GO 9.262623
Atlantic (27.25,27.5] sal tem oxygen silicate phosphate 4.779166 4.504624 3.915107 JGOFS_GO 9.283789
Atlantic (27.25,27.5] tem aou oxygen phosphate 4.776702 4.554342 3.528315 JGOFS_GO 9.331044
Atlantic (27.25,27.5] tem aou oxygen phosphate phosphate_star 4.772737 4.543040 3.421364 JGOFS_GO 9.315776
Atlantic (27.25,27.5] tem aou oxygen silicate phosphate 4.708352 4.554250 3.467433 JGOFS_GO 9.262602
Atlantic (27.5,27.75] aou oxygen silicate phosphate 4.181349 4.384925 4.678199 GO_new 9.063125
Atlantic (27.5,27.75] aou oxygen silicate phosphate phosphate_star 4.165555 4.370717 4.544374 GO_new 8.915091
Atlantic (27.5,27.75] sal aou oxygen silicate phosphate 4.170632 4.369707 4.676151 GO_new 9.045858
Atlantic (27.5,27.75] sal aou silicate phosphate phosphate_star 4.297026 4.451836 4.564505 GO_new 9.016341
Atlantic (27.5,27.75] sal tem aou silicate phosphate 4.157749 4.384594 4.678127 GO_new 9.062721
Atlantic (27.5,27.75] tem aou oxygen phosphate phosphate_star 4.221534 4.545312 4.533429 GO_new 9.078741
Atlantic (27.5,27.75] tem aou oxygen silicate phosphate 4.148689 4.355735 4.666207 GO_new 9.021941
Atlantic (27.5,27.75] tem aou silicate phosphate 4.164381 4.398322 4.681511 GO_new 9.079833
Atlantic (27.5,27.75] tem aou silicate phosphate phosphate_star 4.148548 4.383164 4.549982 GO_new 8.933146
Atlantic (27.5,27.75] tem oxygen silicate phosphate phosphate_star 4.139702 4.400058 4.556720 GO_new 8.956778
Atlantic (27.5,27.75] aou oxygen silicate phosphate 4.181349 4.384925 4.678199 JGOFS_GO 8.566275
Atlantic (27.5,27.75] aou oxygen silicate phosphate phosphate_star 4.165555 4.370717 4.544374 JGOFS_GO 8.536272
Atlantic (27.5,27.75] sal aou oxygen silicate phosphate 4.170632 4.369707 4.676151 JGOFS_GO 8.540339
Atlantic (27.5,27.75] sal tem aou silicate phosphate 4.157749 4.384594 4.678127 JGOFS_GO 8.542343
Atlantic (27.5,27.75] sal tem oxygen silicate phosphate 4.148971 4.399401 4.683535 JGOFS_GO 8.548373
Atlantic (27.5,27.75] tem aou oxygen silicate phosphate 4.148689 4.355735 4.666207 JGOFS_GO 8.504424
Atlantic (27.5,27.75] tem aou silicate phosphate 4.164381 4.398322 4.681511 JGOFS_GO 8.562703
Atlantic (27.5,27.75] tem aou silicate phosphate phosphate_star 4.148548 4.383164 4.549982 JGOFS_GO 8.531712
Atlantic (27.5,27.75] tem oxygen silicate phosphate 4.155021 4.417619 4.687587 JGOFS_GO 8.572640
Atlantic (27.5,27.75] tem oxygen silicate phosphate phosphate_star 4.139702 4.400058 4.556720 JGOFS_GO 8.539760
Atlantic (27.75,27.85] aou oxygen silicate phosphate phosphate_star 4.196722 4.307681 7.705766 GO_new 12.013447
Atlantic (27.75,27.85] sal aou oxygen silicate phosphate 4.113685 4.285353 7.686739 GO_new 11.972092
Atlantic (27.75,27.85] sal aou silicate phosphate 4.120126 4.291980 7.691194 GO_new 11.983174
Atlantic (27.75,27.85] sal aou silicate phosphate phosphate_star 4.106682 4.284806 7.682437 GO_new 11.967243
Atlantic (27.75,27.85] sal oxygen silicate phosphate 4.121204 4.285377 7.689071 GO_new 11.974448
Atlantic (27.75,27.85] sal oxygen silicate phosphate phosphate_star 4.109309 4.277623 7.680367 GO_new 11.957991
Atlantic (27.75,27.85] sal tem aou silicate phosphate 4.114571 4.284915 7.687911 GO_new 11.972827
Atlantic (27.75,27.85] sal tem oxygen silicate phosphate 4.112834 4.285301 7.687608 GO_new 11.972910
Atlantic (27.75,27.85] tem aou silicate phosphate phosphate_star 4.191536 4.306259 7.709321 GO_new 12.015580
Atlantic (27.75,27.85] tem oxygen silicate phosphate phosphate_star 4.188746 4.305171 7.708908 GO_new 12.014080
Atlantic (27.75,27.85] aou oxygen silicate phosphate phosphate_star 4.196722 4.307681 7.705766 JGOFS_GO 8.504403
Atlantic (27.75,27.85] sal aou oxygen silicate phosphate 4.113685 4.285353 7.686739 JGOFS_GO 8.399038
Atlantic (27.75,27.85] sal aou silicate phosphate 4.120126 4.291980 7.691194 JGOFS_GO 8.412106
Atlantic (27.75,27.85] sal aou silicate phosphate phosphate_star 4.106682 4.284806 7.682437 JGOFS_GO 8.391489
Atlantic (27.75,27.85] sal oxygen silicate phosphate 4.121204 4.285377 7.689071 JGOFS_GO 8.406581
Atlantic (27.75,27.85] sal oxygen silicate phosphate phosphate_star 4.109309 4.277623 7.680367 JGOFS_GO 8.386932
Atlantic (27.75,27.85] sal tem aou silicate phosphate 4.114571 4.284915 7.687911 JGOFS_GO 8.399486
Atlantic (27.75,27.85] sal tem oxygen silicate phosphate 4.112834 4.285301 7.687608 JGOFS_GO 8.398135
Atlantic (27.75,27.85] tem aou silicate phosphate phosphate_star 4.191536 4.306259 7.709321 JGOFS_GO 8.497795
Atlantic (27.75,27.85] tem oxygen silicate phosphate phosphate_star 4.188746 4.305171 7.708908 JGOFS_GO 8.493917
Atlantic (27.85,27.95] sal aou oxygen silicate phosphate 3.733203 4.162374 5.126559 GO_new 9.288933
Atlantic (27.85,27.95] sal aou silicate phosphate 3.927851 4.640223 5.475140 GO_new 10.115363
Atlantic (27.85,27.95] sal aou silicate phosphate phosphate_star 3.895442 4.625502 5.473952 GO_new 10.099454
Atlantic (27.85,27.95] sal oxygen silicate phosphate 3.793496 4.463692 5.404331 GO_new 9.868023
Atlantic (27.85,27.95] sal oxygen silicate phosphate phosphate_star 3.774106 4.436970 5.397057 GO_new 9.834028
Atlantic (27.85,27.95] sal tem aou silicate phosphate 3.755350 4.186637 5.158031 GO_new 9.344668
Atlantic (27.85,27.95] sal tem oxygen silicate phosphate 3.744962 4.180818 5.155292 GO_new 9.336111
Atlantic (27.85,27.95] sal tem silicate phosphate 3.981432 4.271866 5.181192 GO_new 9.453058
Atlantic (27.85,27.95] sal tem silicate phosphate phosphate_star 3.879289 4.271813 5.168314 GO_new 9.440127
Atlantic (27.85,27.95] tem aou oxygen silicate phosphate 4.212926 4.831915 5.579193 GO_new 10.411108
Atlantic (27.85,27.95] sal aou oxygen silicate phosphate 3.733203 4.162374 5.126559 JGOFS_GO 7.895576
Atlantic (27.85,27.95] sal aou silicate phosphate 3.927851 4.640223 5.475140 JGOFS_GO 8.568074
Atlantic (27.85,27.95] sal aou silicate phosphate phosphate_star 3.895442 4.625502 5.473952 JGOFS_GO 8.520944
Atlantic (27.85,27.95] sal oxygen silicate phosphate 3.793496 4.463692 5.404331 JGOFS_GO 8.257189
Atlantic (27.85,27.95] sal oxygen silicate phosphate phosphate_star 3.774106 4.436970 5.397057 JGOFS_GO 8.211076
Atlantic (27.85,27.95] sal tem aou silicate phosphate 3.755350 4.186637 5.158031 JGOFS_GO 7.941987
Atlantic (27.85,27.95] sal tem oxygen silicate phosphate 3.744962 4.180818 5.155292 JGOFS_GO 7.925780
Atlantic (27.85,27.95] sal tem silicate phosphate 3.981432 4.271866 5.181192 JGOFS_GO 8.253298
Atlantic (27.85,27.95] sal tem silicate phosphate phosphate_star 3.879289 4.271813 5.168314 JGOFS_GO 8.151102
Atlantic (27.85,27.95] tem aou oxygen silicate phosphate 4.212926 4.831915 5.579193 JGOFS_GO 9.044840
Atlantic (27.95,28.05] sal aou oxygen silicate phosphate 3.487770 4.308349 5.366509 GO_new 9.674857
Atlantic (27.95,28.05] sal aou silicate phosphate 3.757788 4.541320 5.618862 GO_new 10.160182
Atlantic (27.95,28.05] sal aou silicate phosphate phosphate_star 3.750648 4.528583 5.554816 GO_new 10.083399
Atlantic (27.95,28.05] sal oxygen silicate phosphate 3.570686 4.354920 5.524972 GO_new 9.879893
Atlantic (27.95,28.05] sal oxygen silicate phosphate phosphate_star 3.570396 4.330161 5.418256 GO_new 9.748417
Atlantic (27.95,28.05] sal tem aou oxygen phosphate 3.755336 4.353971 4.677615 GO_new 9.031586
Atlantic (27.95,28.05] sal tem aou silicate phosphate 3.508118 4.354175 5.504209 GO_new 9.858384
Atlantic (27.95,28.05] sal tem oxygen silicate phosphate 3.498901 4.330126 5.481394 GO_new 9.811519
Atlantic (27.95,28.05] sal tem silicate phosphate phosphate_star 3.755808 4.709496 5.438663 GO_new 10.148158
Atlantic (27.95,28.05] tem aou oxygen silicate phosphate 4.117187 4.860501 5.227002 GO_new 10.087504
Atlantic (27.95,28.05] sal aou oxygen silicate phosphate 3.487770 4.308349 5.366509 JGOFS_GO 7.796119
Atlantic (27.95,28.05] sal aou silicate phosphate 3.757788 4.541320 5.618862 JGOFS_GO 8.299108
Atlantic (27.95,28.05] sal aou silicate phosphate phosphate_star 3.750648 4.528583 5.554816 JGOFS_GO 8.279231
Atlantic (27.95,28.05] sal oxygen silicate phosphate 3.570686 4.354920 5.524972 JGOFS_GO 7.925607
Atlantic (27.95,28.05] sal oxygen silicate phosphate phosphate_star 3.570396 4.330161 5.418256 JGOFS_GO 7.900557
Atlantic (27.95,28.05] sal tem aou oxygen phosphate 3.755336 4.353971 4.677615 JGOFS_GO 8.109307
Atlantic (27.95,28.05] sal tem aou silicate phosphate 3.508118 4.354175 5.504209 JGOFS_GO 7.862294
Atlantic (27.95,28.05] sal tem oxygen phosphate 3.776519 4.617461 5.753398 JGOFS_GO 8.393980
Atlantic (27.95,28.05] sal tem oxygen phosphate phosphate_star 3.771038 4.614698 5.706177 JGOFS_GO 8.385736
Atlantic (27.95,28.05] sal tem oxygen silicate phosphate 3.498901 4.330126 5.481394 JGOFS_GO 7.829027
Atlantic (28.05,28.1] sal tem aou oxygen phosphate 3.539499 4.204534 4.630652 GO_new 8.835185
Atlantic (28.05,28.1] sal tem aou phosphate 3.607164 4.569347 5.205477 GO_new 9.774825
Atlantic (28.05,28.1] sal tem aou phosphate phosphate_star 3.604413 4.562851 5.185582 GO_new 9.748433
Atlantic (28.05,28.1] sal tem aou silicate phosphate 3.586755 4.507157 5.189566 GO_new 9.696724
Atlantic (28.05,28.1] sal tem oxygen phosphate 3.557078 4.445505 5.152437 GO_new 9.597943
Atlantic (28.05,28.1] sal tem oxygen phosphate phosphate_star 3.554312 4.438866 5.121261 GO_new 9.560127
Atlantic (28.05,28.1] sal tem oxygen silicate phosphate 3.531111 4.393917 5.124357 GO_new 9.518275
Atlantic (28.05,28.1] tem aou oxygen phosphate 3.737933 4.326767 4.638593 GO_new 8.965360
Atlantic (28.05,28.1] tem aou oxygen phosphate phosphate_star 3.737328 4.326389 4.601921 GO_new 8.928310
Atlantic (28.05,28.1] tem aou oxygen silicate phosphate 3.643192 4.179970 4.630074 GO_new 8.810044
Atlantic (28.05,28.1] sal aou oxygen silicate phosphate 3.554108 4.406193 5.411507 JGOFS_GO 7.960301
Atlantic (28.05,28.1] sal oxygen silicate phosphate 3.555467 4.420889 5.419153 JGOFS_GO 7.976356
Atlantic (28.05,28.1] sal oxygen silicate phosphate phosphate_star 3.548422 4.401627 5.386608 JGOFS_GO 7.950049
Atlantic (28.05,28.1] sal tem aou oxygen phosphate 3.539499 4.204534 4.630652 JGOFS_GO 7.744033
Atlantic (28.05,28.1] sal tem oxygen phosphate 3.557078 4.445505 5.152437 JGOFS_GO 8.002583
Atlantic (28.05,28.1] sal tem oxygen phosphate phosphate_star 3.554312 4.438866 5.121261 JGOFS_GO 7.993178
Atlantic (28.05,28.1] sal tem oxygen silicate phosphate 3.531111 4.393917 5.124357 JGOFS_GO 7.925028
Atlantic (28.05,28.1] tem aou oxygen phosphate 3.737933 4.326767 4.638593 JGOFS_GO 8.064700
Atlantic (28.05,28.1] tem aou oxygen phosphate phosphate_star 3.737328 4.326389 4.601921 JGOFS_GO 8.063717
Atlantic (28.05,28.1] tem aou oxygen silicate phosphate 3.643192 4.179970 4.630074 JGOFS_GO 7.823162
Atlantic (28.1,28.15] sal tem aou oxygen phosphate 3.453382 3.943365 4.691597 GO_new 8.634961
Atlantic (28.1,28.15] sal tem aou silicate phosphate 3.395913 4.012656 4.782841 GO_new 8.795498
Atlantic (28.1,28.15] sal tem oxygen phosphate 3.459589 3.960398 4.753057 GO_new 8.713455
Atlantic (28.1,28.15] sal tem oxygen phosphate phosphate_star 3.459522 3.960389 4.749213 GO_new 8.709602
Atlantic (28.1,28.15] sal tem oxygen silicate phosphate 3.349985 3.934786 4.738195 GO_new 8.672981
Atlantic (28.1,28.15] tem aou oxygen phosphate 3.547872 3.985209 4.716845 GO_new 8.702054
Atlantic (28.1,28.15] tem aou oxygen phosphate phosphate_star 3.547716 3.984606 4.715557 GO_new 8.700162
Atlantic (28.1,28.15] tem aou oxygen silicate phosphate 3.434468 3.944928 4.715162 GO_new 8.660090
Atlantic (28.1,28.15] tem oxygen silicate phosphate 3.434707 3.965040 4.799940 GO_new 8.764981
Atlantic (28.1,28.15] tem oxygen silicate phosphate phosphate_star 3.427973 3.964974 4.797528 GO_new 8.762502
Atlantic (28.1,28.15] sal tem aou oxygen phosphate 3.453382 3.943365 4.691597 JGOFS_GO 7.396747
Atlantic (28.1,28.15] sal tem aou silicate phosphate 3.395913 4.012656 4.782841 JGOFS_GO 7.408570
Atlantic (28.1,28.15] sal tem oxygen phosphate 3.459589 3.960398 4.753057 JGOFS_GO 7.419988
Atlantic (28.1,28.15] sal tem oxygen phosphate phosphate_star 3.459522 3.960389 4.749213 JGOFS_GO 7.419911
Atlantic (28.1,28.15] sal tem oxygen silicate phosphate 3.349985 3.934786 4.738195 JGOFS_GO 7.284771
Atlantic (28.1,28.15] tem aou oxygen silicate phosphate 3.434468 3.944928 4.715162 JGOFS_GO 7.379396
Atlantic (28.1,28.15] tem aou silicate phosphate 3.459579 4.030995 4.865916 JGOFS_GO 7.490573
Atlantic (28.1,28.15] tem aou silicate phosphate phosphate_star 3.454261 4.030922 4.856988 JGOFS_GO 7.485182
Atlantic (28.1,28.15] tem oxygen silicate phosphate 3.434707 3.965040 4.799940 JGOFS_GO 7.399747
Atlantic (28.1,28.15] tem oxygen silicate phosphate phosphate_star 3.427973 3.964974 4.797528 JGOFS_GO 7.392947
Atlantic (28.15,28.2] sal aou oxygen silicate phosphate 3.726384 3.637891 3.023471 GO_new 6.661362
Atlantic (28.15,28.2] sal tem aou oxygen phosphate 3.726099 3.629154 3.004142 GO_new 6.633295
Atlantic (28.15,28.2] sal tem aou phosphate 3.741007 3.656291 3.004525 GO_new 6.660816
Atlantic (28.15,28.2] sal tem aou phosphate phosphate_star 3.732000 3.655905 2.997692 GO_new 6.653596
Atlantic (28.15,28.2] sal tem aou silicate phosphate 3.739605 3.638041 2.989969 GO_new 6.628010
Atlantic (28.15,28.2] sal tem oxygen phosphate 3.781131 3.701975 3.004907 GO_new 6.706882
Atlantic (28.15,28.2] sal tem oxygen phosphate phosphate_star 3.775184 3.701111 2.996741 GO_new 6.697852
Atlantic (28.15,28.2] sal tem oxygen silicate phosphate 3.770684 3.671003 2.995966 GO_new 6.666968
Atlantic (28.15,28.2] tem aou oxygen phosphate phosphate_star 3.738522 3.641657 3.089357 GO_new 6.731013
Atlantic (28.15,28.2] tem aou oxygen silicate phosphate 3.748534 3.636938 3.094567 GO_new 6.731506
Atlantic (28.15,28.2] sal aou oxygen phosphate 3.727404 3.639684 3.100446 JGOFS_GO 7.367088
Atlantic (28.15,28.2] sal aou oxygen phosphate phosphate_star 3.719971 3.639684 3.099907 JGOFS_GO 7.359655
Atlantic (28.15,28.2] sal aou oxygen silicate phosphate 3.726384 3.637891 3.023471 JGOFS_GO 7.364276
Atlantic (28.15,28.2] sal tem aou oxygen phosphate 3.726099 3.629154 3.004142 JGOFS_GO 7.355253
Atlantic (28.15,28.2] sal tem aou phosphate 3.741007 3.656291 3.004525 JGOFS_GO 7.397297
Atlantic (28.15,28.2] sal tem aou phosphate phosphate_star 3.732000 3.655905 2.997692 JGOFS_GO 7.387904
Atlantic (28.15,28.2] sal tem aou silicate phosphate 3.739605 3.638041 2.989969 JGOFS_GO 7.377646
Atlantic (28.15,28.2] tem aou oxygen phosphate 3.749523 3.641831 3.097179 JGOFS_GO 7.391355
Atlantic (28.15,28.2] tem aou oxygen phosphate phosphate_star 3.738522 3.641657 3.089357 JGOFS_GO 7.380178
Atlantic (28.15,28.2] tem aou oxygen silicate phosphate 3.748534 3.636938 3.094567 JGOFS_GO 7.385473
Atlantic (28.2, Inf] sal aou oxygen phosphate 3.549135 3.174916 3.419848 GO_new 6.594764
Atlantic (28.2, Inf] sal aou oxygen phosphate phosphate_star 3.546339 3.160649 3.411306 GO_new 6.571955
Atlantic (28.2, Inf] sal aou oxygen silicate phosphate 3.487877 3.158126 3.393186 GO_new 6.551312
Atlantic (28.2, Inf] sal oxygen phosphate phosphate_star 3.978629 3.172926 3.414075 GO_new 6.587001
Atlantic (28.2, Inf] sal oxygen silicate phosphate 3.765786 3.165659 3.395547 GO_new 6.561207
Atlantic (28.2, Inf] sal oxygen silicate phosphate phosphate_star 3.756399 3.137716 3.392982 GO_new 6.530697
Atlantic (28.2, Inf] sal tem aou oxygen phosphate 3.522234 3.173511 3.190100 GO_new 6.363611
Atlantic (28.2, Inf] sal tem oxygen phosphate 3.651581 3.173715 3.346153 GO_new 6.519868
Atlantic (28.2, Inf] sal tem oxygen phosphate phosphate_star 3.649876 3.162480 3.345556 GO_new 6.508036
Atlantic (28.2, Inf] sal tem oxygen silicate phosphate 3.544433 3.152200 3.338683 GO_new 6.490883
Atlantic (28.2, Inf] sal aou oxygen phosphate 3.549135 3.174916 3.419848 JGOFS_GO 6.724051
Atlantic (28.2, Inf] sal aou oxygen phosphate phosphate_star 3.546339 3.160649 3.411306 JGOFS_GO 6.706988
Atlantic (28.2, Inf] sal aou oxygen silicate phosphate 3.487877 3.158126 3.393186 JGOFS_GO 6.646004
Atlantic (28.2, Inf] sal aou silicate phosphate 3.530184 3.237823 3.463376 JGOFS_GO 6.768007
Atlantic (28.2, Inf] sal aou silicate phosphate phosphate_star 3.518405 3.220199 3.456579 JGOFS_GO 6.738604
Atlantic (28.2, Inf] sal tem aou oxygen phosphate 3.522234 3.173511 3.190100 JGOFS_GO 6.695745
Atlantic (28.2, Inf] sal tem aou phosphate 3.568388 3.192194 3.467972 JGOFS_GO 6.760581
Atlantic (28.2, Inf] sal tem aou phosphate phosphate_star 3.565166 3.171852 3.465054 JGOFS_GO 6.737018
Atlantic (28.2, Inf] sal tem aou silicate phosphate 3.499887 3.184788 3.462424 JGOFS_GO 6.684675
Atlantic (28.2, Inf] sal tem oxygen silicate phosphate 3.544433 3.152200 3.338683 JGOFS_GO 6.696633
Indo-Pacific (-Inf,26] sal aou oxygen phosphate phosphate_star 11.656608 6.900704 4.748256 GO_new 11.648960
Indo-Pacific (-Inf,26] sal aou oxygen silicate phosphate 11.589172 7.158929 5.649148 GO_new 12.808077
Indo-Pacific (-Inf,26] sal tem aou phosphate phosphate_star 11.727854 6.878761 4.968427 GO_new 11.847189
Indo-Pacific (-Inf,26] sal tem aou silicate phosphate 11.658578 7.235892 5.852129 GO_new 13.088021
Indo-Pacific (-Inf,26] sal tem oxygen phosphate phosphate_star 11.772722 6.912573 5.114563 GO_new 12.027136
Indo-Pacific (-Inf,26] sal tem oxygen silicate phosphate 11.696367 7.284738 5.952533 GO_new 13.237271
Indo-Pacific (-Inf,26] tem aou oxygen phosphate phosphate_star 12.336939 7.286502 5.639494 GO_new 12.925996
Indo-Pacific (-Inf,26] tem aou silicate phosphate phosphate_star 12.281655 7.464928 5.655524 GO_new 13.120452
Indo-Pacific (-Inf,26] tem oxygen phosphate phosphate_star 12.406986 7.414955 5.643696 GO_new 13.058651
Indo-Pacific (-Inf,26] tem oxygen silicate phosphate phosphate_star 12.209892 7.343438 5.628143 GO_new 12.971580
Indo-Pacific (-Inf,26] sal aou oxygen phosphate 11.847558 7.595246 6.038692 JGOFS_GO 19.442804
Indo-Pacific (-Inf,26] sal aou oxygen phosphate phosphate_star 11.656608 6.900704 4.748256 JGOFS_GO 18.557312
Indo-Pacific (-Inf,26] sal aou oxygen silicate phosphate 11.589172 7.158929 5.649148 JGOFS_GO 18.748101
Indo-Pacific (-Inf,26] sal tem aou oxygen phosphate 11.835399 7.592149 6.001991 JGOFS_GO 19.427548
Indo-Pacific (-Inf,26] sal tem aou phosphate 11.936886 7.625869 6.165843 JGOFS_GO 19.562755
Indo-Pacific (-Inf,26] sal tem aou phosphate phosphate_star 11.727854 6.878761 4.968427 JGOFS_GO 18.606615
Indo-Pacific (-Inf,26] sal tem aou silicate phosphate 11.658578 7.235892 5.852129 JGOFS_GO 18.894470
Indo-Pacific (-Inf,26] sal tem oxygen phosphate phosphate_star 11.772722 6.912573 5.114563 JGOFS_GO 18.685295
Indo-Pacific (-Inf,26] sal tem oxygen silicate phosphate 11.696367 7.284738 5.952533 JGOFS_GO 18.981104
Indo-Pacific (-Inf,26] tem oxygen silicate phosphate phosphate_star 12.209892 7.343438 5.628143 JGOFS_GO 19.553330
Indo-Pacific (26,26.5] sal aou oxygen phosphate 7.963659 5.384444 4.264629 GO_new 9.649072
Indo-Pacific (26,26.5] sal aou oxygen phosphate phosphate_star 7.817898 5.209679 3.993262 GO_new 9.202941
Indo-Pacific (26,26.5] sal aou oxygen silicate phosphate 7.855393 5.185799 3.976710 GO_new 9.162509
Indo-Pacific (26,26.5] sal tem aou oxygen phosphate 7.962543 5.283401 4.258474 GO_new 9.541875
Indo-Pacific (26,26.5] sal tem aou phosphate 8.017736 5.295081 4.301171 GO_new 9.596252
Indo-Pacific (26,26.5] sal tem aou phosphate phosphate_star 7.880514 5.108916 4.035301 GO_new 9.144218
Indo-Pacific (26,26.5] sal tem aou silicate phosphate 7.929705 5.143958 4.051154 GO_new 9.195112
Indo-Pacific (26,26.5] sal tem oxygen phosphate 8.048135 5.284498 4.321445 GO_new 9.605943
Indo-Pacific (26,26.5] sal tem oxygen phosphate phosphate_star 7.921449 5.101898 4.066271 GO_new 9.168169
Indo-Pacific (26,26.5] sal tem oxygen silicate phosphate 7.966619 5.153265 4.088387 GO_new 9.241652
Indo-Pacific (26,26.5] sal aou oxygen phosphate 7.963659 5.384444 4.264629 JGOFS_GO 13.348102
Indo-Pacific (26,26.5] sal aou oxygen phosphate phosphate_star 7.817898 5.209679 3.993262 JGOFS_GO 13.027577
Indo-Pacific (26,26.5] sal aou oxygen silicate phosphate 7.855393 5.185799 3.976710 JGOFS_GO 13.041193
Indo-Pacific (26,26.5] sal tem aou oxygen phosphate 7.962543 5.283401 4.258474 JGOFS_GO 13.245944
Indo-Pacific (26,26.5] sal tem aou phosphate 8.017736 5.295081 4.301171 JGOFS_GO 13.312817
Indo-Pacific (26,26.5] sal tem aou phosphate phosphate_star 7.880514 5.108916 4.035301 JGOFS_GO 12.989430
Indo-Pacific (26,26.5] sal tem aou silicate phosphate 7.929705 5.143958 4.051154 JGOFS_GO 13.073663
Indo-Pacific (26,26.5] sal tem oxygen phosphate 8.048135 5.284498 4.321445 JGOFS_GO 13.332633
Indo-Pacific (26,26.5] sal tem oxygen phosphate phosphate_star 7.921449 5.101898 4.066271 JGOFS_GO 13.023347
Indo-Pacific (26,26.5] sal tem oxygen silicate phosphate 7.966619 5.153265 4.088387 JGOFS_GO 13.119885
Indo-Pacific (26.5,26.75] aou oxygen silicate phosphate phosphate_star 6.122193 5.842826 4.471119 GO_new 10.313944
Indo-Pacific (26.5,26.75] sal aou oxygen phosphate 5.972964 5.696797 4.631952 GO_new 10.328749
Indo-Pacific (26.5,26.75] sal aou oxygen phosphate phosphate_star 5.629568 5.431213 4.210965 GO_new 9.642178
Indo-Pacific (26.5,26.75] sal aou oxygen silicate phosphate 5.889764 5.691859 4.563197 GO_new 10.255056
Indo-Pacific (26.5,26.75] sal tem aou oxygen phosphate 5.970854 5.696790 4.615782 GO_new 10.312572
Indo-Pacific (26.5,26.75] sal tem aou phosphate phosphate_star 5.730080 5.471410 4.292707 GO_new 9.764117
Indo-Pacific (26.5,26.75] sal tem aou silicate phosphate 6.032011 5.709349 4.659652 GO_new 10.369001
Indo-Pacific (26.5,26.75] sal tem oxygen phosphate phosphate_star 5.794994 5.512058 4.353305 GO_new 9.865362
Indo-Pacific (26.5,26.75] tem aou oxygen phosphate phosphate_star 6.099615 5.841145 4.577322 GO_new 10.418466
Indo-Pacific (26.5,26.75] tem aou silicate phosphate phosphate_star 6.164107 5.854419 4.582028 GO_new 10.436447
Indo-Pacific (26.5,26.75] sal aou oxygen phosphate 5.972964 5.696797 4.631952 JGOFS_GO 11.669761
Indo-Pacific (26.5,26.75] sal aou oxygen phosphate phosphate_star 5.629568 5.431213 4.210965 JGOFS_GO 11.060781
Indo-Pacific (26.5,26.75] sal aou oxygen silicate phosphate 5.889764 5.691859 4.563197 JGOFS_GO 11.581623
Indo-Pacific (26.5,26.75] sal tem aou oxygen phosphate 5.970854 5.696790 4.615782 JGOFS_GO 11.667644
Indo-Pacific (26.5,26.75] sal tem aou phosphate 6.059163 5.732700 4.712027 JGOFS_GO 11.791863
Indo-Pacific (26.5,26.75] sal tem aou phosphate phosphate_star 5.730080 5.471410 4.292707 JGOFS_GO 11.201489
Indo-Pacific (26.5,26.75] sal tem aou silicate phosphate 6.032011 5.709349 4.659652 JGOFS_GO 11.741360
Indo-Pacific (26.5,26.75] sal tem oxygen phosphate 6.103238 5.761208 4.756842 JGOFS_GO 11.864446
Indo-Pacific (26.5,26.75] sal tem oxygen phosphate phosphate_star 5.794994 5.512058 4.353305 JGOFS_GO 11.307052
Indo-Pacific (26.5,26.75] sal tem oxygen silicate phosphate 6.085625 5.726657 4.715326 JGOFS_GO 11.812282
Indo-Pacific (26.75,27] aou oxygen phosphate phosphate_star 6.567840 5.032576 4.161568 GO_new 9.194144
Indo-Pacific (26.75,27] aou oxygen silicate phosphate phosphate_star 6.488004 4.974323 4.160697 GO_new 9.135019
Indo-Pacific (26.75,27] sal aou oxygen phosphate phosphate_star 6.533269 4.972238 4.063475 GO_new 9.035714
Indo-Pacific (26.75,27] sal aou oxygen silicate phosphate 6.552579 5.062787 4.421716 GO_new 9.484503
Indo-Pacific (26.75,27] sal tem aou oxygen phosphate 6.516743 5.021419 4.250318 GO_new 9.271737
Indo-Pacific (26.75,27] sal tem aou phosphate phosphate_star 6.645887 5.160774 4.264093 GO_new 9.424867
Indo-Pacific (26.75,27] tem aou oxygen phosphate 6.520101 5.023021 4.264277 GO_new 9.287298
Indo-Pacific (26.75,27] tem aou oxygen phosphate phosphate_star 6.482046 4.918858 4.034617 GO_new 8.953476
Indo-Pacific (26.75,27] tem aou oxygen silicate phosphate 6.492756 4.983055 4.235837 GO_new 9.218892
Indo-Pacific (26.75,27] tem aou silicate phosphate phosphate_star 6.519436 5.038575 4.289332 GO_new 9.327907
Indo-Pacific (26.75,27] aou oxygen phosphate phosphate_star 6.567840 5.032576 4.161568 JGOFS_GO 11.600415
Indo-Pacific (26.75,27] aou oxygen silicate phosphate phosphate_star 6.488004 4.974323 4.160697 JGOFS_GO 11.462326
Indo-Pacific (26.75,27] sal aou oxygen phosphate phosphate_star 6.533269 4.972238 4.063475 JGOFS_GO 11.505507
Indo-Pacific (26.75,27] sal aou oxygen silicate phosphate 6.552579 5.062787 4.421716 JGOFS_GO 11.615365
Indo-Pacific (26.75,27] sal tem aou oxygen phosphate 6.516743 5.021419 4.250318 JGOFS_GO 11.538162
Indo-Pacific (26.75,27] tem aou oxygen phosphate 6.520101 5.023021 4.264277 JGOFS_GO 11.543122
Indo-Pacific (26.75,27] tem aou oxygen phosphate phosphate_star 6.482046 4.918858 4.034617 JGOFS_GO 11.400905
Indo-Pacific (26.75,27] tem aou oxygen silicate phosphate 6.492756 4.983055 4.235837 JGOFS_GO 11.475811
Indo-Pacific (26.75,27] tem aou silicate phosphate phosphate_star 6.519436 5.038575 4.289332 JGOFS_GO 11.558011
Indo-Pacific (26.75,27] tem oxygen silicate phosphate phosphate_star 6.555616 5.102646 4.383563 JGOFS_GO 11.658262
Indo-Pacific (27,27.25] sal aou oxygen phosphate phosphate_star 6.376268 4.789831 4.178131 GO_new 8.967962
Indo-Pacific (27,27.25] sal aou oxygen silicate phosphate 6.377139 4.610352 4.154787 GO_new 8.765140
Indo-Pacific (27,27.25] sal aou silicate phosphate 6.502600 4.621347 4.254518 GO_new 8.875865
Indo-Pacific (27,27.25] sal aou silicate phosphate phosphate_star 6.447889 4.618712 4.254365 GO_new 8.873077
Indo-Pacific (27,27.25] sal tem aou oxygen phosphate 5.996869 4.633139 4.167623 GO_new 8.800762
Indo-Pacific (27,27.25] sal tem aou silicate phosphate 6.492264 4.620740 4.179057 GO_new 8.799797
Indo-Pacific (27,27.25] sal tem oxygen silicate phosphate 6.598846 4.662267 4.211183 GO_new 8.873450
Indo-Pacific (27,27.25] tem aou oxygen phosphate 6.024170 4.748844 4.304139 GO_new 9.052983
Indo-Pacific (27,27.25] tem aou oxygen phosphate phosphate_star 6.022258 4.655002 4.175739 GO_new 8.830741
Indo-Pacific (27,27.25] tem aou oxygen silicate phosphate 5.981049 4.537003 4.124776 GO_new 8.661779
Indo-Pacific (27,27.25] aou oxygen silicate phosphate phosphate_star 6.295462 4.772678 4.374179 JGOFS_GO 11.068141
Indo-Pacific (27,27.25] sal aou oxygen phosphate phosphate_star 6.376268 4.789831 4.178131 JGOFS_GO 11.166100
Indo-Pacific (27,27.25] sal aou oxygen silicate phosphate 6.377139 4.610352 4.154787 JGOFS_GO 10.987491
Indo-Pacific (27,27.25] sal aou silicate phosphate 6.502600 4.621347 4.254518 JGOFS_GO 11.123948
Indo-Pacific (27,27.25] sal aou silicate phosphate phosphate_star 6.447889 4.618712 4.254365 JGOFS_GO 11.066601
Indo-Pacific (27,27.25] sal tem aou oxygen phosphate 5.996869 4.633139 4.167623 JGOFS_GO 10.630008
Indo-Pacific (27,27.25] sal tem aou silicate phosphate 6.492264 4.620740 4.179057 JGOFS_GO 11.113004
Indo-Pacific (27,27.25] tem aou oxygen phosphate 6.024170 4.748844 4.304139 JGOFS_GO 10.773014
Indo-Pacific (27,27.25] tem aou oxygen phosphate phosphate_star 6.022258 4.655002 4.175739 JGOFS_GO 10.677260
Indo-Pacific (27,27.25] tem aou oxygen silicate phosphate 5.981049 4.537003 4.124776 JGOFS_GO 10.518052
Indo-Pacific (27.25,27.5] aou oxygen silicate phosphate phosphate_star 5.196364 4.240548 3.681856 GO_new 7.922403
Indo-Pacific (27.25,27.5] sal aou oxygen phosphate phosphate_star 5.168633 4.254234 3.572729 GO_new 7.826964
Indo-Pacific (27.25,27.5] sal aou oxygen silicate phosphate 5.216429 4.136611 3.512906 GO_new 7.649517
Indo-Pacific (27.25,27.5] sal tem aou oxygen phosphate 4.831746 4.162230 3.629966 GO_new 7.792197
Indo-Pacific (27.25,27.5] sal tem aou phosphate phosphate_star 5.509279 4.347958 3.586953 GO_new 7.934910
Indo-Pacific (27.25,27.5] sal tem aou silicate phosphate 5.541110 4.185755 3.535978 GO_new 7.721733
Indo-Pacific (27.25,27.5] sal tem oxygen silicate phosphate 5.764252 4.233532 3.556430 GO_new 7.789962
Indo-Pacific (27.25,27.5] tem aou oxygen phosphate 4.835593 4.192637 3.734596 GO_new 7.927233
Indo-Pacific (27.25,27.5] tem aou oxygen phosphate phosphate_star 4.832492 4.093420 3.684443 GO_new 7.777862
Indo-Pacific (27.25,27.5] tem aou oxygen silicate phosphate 4.804957 3.986913 3.603040 GO_new 7.589953
Indo-Pacific (27.25,27.5] aou oxygen phosphate phosphate_star 5.199698 4.327634 3.762548 JGOFS_GO 9.527332
Indo-Pacific (27.25,27.5] aou oxygen silicate phosphate 5.256328 4.307021 3.757882 JGOFS_GO 9.563349
Indo-Pacific (27.25,27.5] aou oxygen silicate phosphate phosphate_star 5.196364 4.240548 3.681856 JGOFS_GO 9.436912
Indo-Pacific (27.25,27.5] sal aou oxygen phosphate 5.224552 4.399037 3.638993 JGOFS_GO 9.623589
Indo-Pacific (27.25,27.5] sal aou oxygen phosphate phosphate_star 5.168633 4.254234 3.572729 JGOFS_GO 9.422867
Indo-Pacific (27.25,27.5] sal aou oxygen silicate phosphate 5.216429 4.136611 3.512906 JGOFS_GO 9.353040
Indo-Pacific (27.25,27.5] sal tem aou oxygen phosphate 4.831746 4.162230 3.629966 JGOFS_GO 8.993977
Indo-Pacific (27.25,27.5] tem aou oxygen phosphate 4.835593 4.192637 3.734596 JGOFS_GO 9.028231
Indo-Pacific (27.25,27.5] tem aou oxygen phosphate phosphate_star 4.832492 4.093420 3.684443 JGOFS_GO 8.925912
Indo-Pacific (27.25,27.5] tem aou oxygen silicate phosphate 4.804957 3.986913 3.603040 JGOFS_GO 8.791870
Indo-Pacific (27.5,27.75] sal aou oxygen phosphate phosphate_star 4.656269 3.907055 3.757383 GO_new 7.664438
Indo-Pacific (27.5,27.75] sal aou oxygen silicate phosphate 4.655491 3.997240 3.735975 GO_new 7.733215
Indo-Pacific (27.5,27.75] sal tem aou oxygen phosphate 4.409996 4.005957 3.525039 GO_new 7.530996
Indo-Pacific (27.5,27.75] sal tem aou phosphate 4.826836 4.047735 3.639034 GO_new 7.686769
Indo-Pacific (27.5,27.75] sal tem aou phosphate phosphate_star 4.826586 3.894715 3.637182 GO_new 7.531897
Indo-Pacific (27.5,27.75] sal tem aou silicate phosphate 4.822537 4.035117 3.627796 GO_new 7.662913
Indo-Pacific (27.5,27.75] sal tem oxygen phosphate 4.957364 4.095366 3.574533 GO_new 7.669899
Indo-Pacific (27.5,27.75] sal tem oxygen phosphate phosphate_star 4.957076 3.900138 3.561704 GO_new 7.461841
Indo-Pacific (27.5,27.75] sal tem oxygen silicate phosphate 4.949934 4.083693 3.565430 GO_new 7.649123
Indo-Pacific (27.5,27.75] tem aou oxygen phosphate phosphate_star 4.409247 3.996479 3.684163 GO_new 7.680643
Indo-Pacific (27.5,27.75] aou oxygen phosphate phosphate_star 4.663346 3.997608 3.906815 JGOFS_GO 8.660953
Indo-Pacific (27.5,27.75] aou oxygen silicate phosphate phosphate_star 4.660946 3.997539 3.901165 JGOFS_GO 8.658485
Indo-Pacific (27.5,27.75] sal aou oxygen phosphate 4.656383 4.011924 3.758951 JGOFS_GO 8.668307
Indo-Pacific (27.5,27.75] sal aou oxygen phosphate phosphate_star 4.656269 3.907055 3.757383 JGOFS_GO 8.563324
Indo-Pacific (27.5,27.75] sal aou oxygen silicate phosphate 4.655491 3.997240 3.735975 JGOFS_GO 8.652731
Indo-Pacific (27.5,27.75] sal tem aou oxygen phosphate 4.409996 4.005957 3.525039 JGOFS_GO 8.415953
Indo-Pacific (27.5,27.75] sal tem aou phosphate phosphate_star 4.826586 3.894715 3.637182 JGOFS_GO 8.721300
Indo-Pacific (27.5,27.75] tem aou oxygen phosphate 4.410507 4.074666 3.760424 JGOFS_GO 8.485173
Indo-Pacific (27.5,27.75] tem aou oxygen phosphate phosphate_star 4.409247 3.996479 3.684163 JGOFS_GO 8.405726
Indo-Pacific (27.5,27.75] tem aou oxygen silicate phosphate 4.386458 4.074656 3.738199 JGOFS_GO 8.461113
Indo-Pacific (27.75,27.85] aou oxygen silicate phosphate 4.249542 3.356084 2.945707 GO_new 6.301792
Indo-Pacific (27.75,27.85] aou oxygen silicate phosphate phosphate_star 4.245475 3.332871 2.938537 GO_new 6.271408
Indo-Pacific (27.75,27.85] sal aou oxygen silicate phosphate 4.235175 3.347718 2.896917 GO_new 6.244636
Indo-Pacific (27.75,27.85] sal tem aou silicate phosphate 4.258533 3.311684 2.876268 GO_new 6.187952
Indo-Pacific (27.75,27.85] sal tem oxygen silicate phosphate 4.281560 3.291064 2.865015 GO_new 6.156079
Indo-Pacific (27.75,27.85] tem aou oxygen silicate phosphate 4.242882 3.271777 2.903747 GO_new 6.175523
Indo-Pacific (27.75,27.85] tem aou silicate phosphate 4.275388 3.316085 2.924736 GO_new 6.240820
Indo-Pacific (27.75,27.85] tem aou silicate phosphate phosphate_star 4.273516 3.302046 2.918998 GO_new 6.221044
Indo-Pacific (27.75,27.85] tem oxygen silicate phosphate 4.300570 3.296323 2.916579 GO_new 6.212903
Indo-Pacific (27.75,27.85] tem oxygen silicate phosphate phosphate_star 4.299485 3.286256 2.911216 GO_new 6.197472
Indo-Pacific (27.75,27.85] aou oxygen silicate phosphate phosphate_star 4.245475 3.332871 2.938537 JGOFS_GO 7.578346
Indo-Pacific (27.75,27.85] sal aou oxygen silicate phosphate 4.235175 3.347718 2.896917 JGOFS_GO 7.582893
Indo-Pacific (27.75,27.85] sal tem aou oxygen phosphate 4.233667 3.327161 3.001071 JGOFS_GO 7.560827
Indo-Pacific (27.75,27.85] sal tem aou silicate phosphate 4.258533 3.311684 2.876268 JGOFS_GO 7.570217
Indo-Pacific (27.75,27.85] sal tem oxygen silicate phosphate 4.281560 3.291064 2.865015 JGOFS_GO 7.572624
Indo-Pacific (27.75,27.85] tem aou oxygen phosphate 4.242883 3.327650 3.002873 JGOFS_GO 7.570533
Indo-Pacific (27.75,27.85] tem aou oxygen phosphate phosphate_star 4.234736 3.327543 3.000201 JGOFS_GO 7.562279
Indo-Pacific (27.75,27.85] tem aou oxygen silicate phosphate 4.242882 3.271777 2.903747 JGOFS_GO 7.514659
Indo-Pacific (27.75,27.85] tem aou silicate phosphate phosphate_star 4.273516 3.302046 2.918998 JGOFS_GO 7.575562
Indo-Pacific (27.75,27.85] tem oxygen silicate phosphate phosphate_star 4.299485 3.286256 2.911216 JGOFS_GO 7.585741
Indo-Pacific (27.85,27.95] aou oxygen silicate phosphate 4.280926 3.007769 2.912915 GO_new 5.920684
Indo-Pacific (27.85,27.95] aou oxygen silicate phosphate phosphate_star 4.256269 2.993750 2.902527 GO_new 5.896277
Indo-Pacific (27.85,27.95] sal aou oxygen silicate phosphate 4.272593 3.007732 2.896370 GO_new 5.904102
Indo-Pacific (27.85,27.95] sal tem aou silicate phosphate 4.290311 3.000580 2.878296 GO_new 5.878876
Indo-Pacific (27.85,27.95] sal tem oxygen silicate phosphate 4.311405 3.004330 2.869502 GO_new 5.873832
Indo-Pacific (27.85,27.95] tem aou oxygen silicate phosphate 4.280774 3.000763 2.891621 GO_new 5.892384
Indo-Pacific (27.85,27.95] tem aou silicate phosphate 4.295926 3.000801 2.897066 GO_new 5.897867
Indo-Pacific (27.85,27.95] tem aou silicate phosphate phosphate_star 4.274004 2.989254 2.887257 GO_new 5.876511
Indo-Pacific (27.85,27.95] tem oxygen silicate phosphate 4.315597 3.005155 2.892410 GO_new 5.897565
Indo-Pacific (27.85,27.95] tem oxygen silicate phosphate phosphate_star 4.294333 2.994511 2.882504 GO_new 5.877015
Indo-Pacific (27.85,27.95] aou oxygen phosphate phosphate_star 4.256599 3.033388 2.960681 JGOFS_GO 7.289987
Indo-Pacific (27.85,27.95] aou oxygen silicate phosphate 4.280926 3.007769 2.912915 JGOFS_GO 7.288694
Indo-Pacific (27.85,27.95] aou oxygen silicate phosphate phosphate_star 4.256269 2.993750 2.902527 JGOFS_GO 7.250018
Indo-Pacific (27.85,27.95] sal aou oxygen phosphate phosphate_star 4.237357 3.016584 2.959140 JGOFS_GO 7.253941
Indo-Pacific (27.85,27.95] sal aou oxygen silicate phosphate 4.272593 3.007732 2.896370 JGOFS_GO 7.280326
Indo-Pacific (27.85,27.95] sal tem aou phosphate phosphate_star 4.256261 3.030864 2.969139 JGOFS_GO 7.287125
Indo-Pacific (27.85,27.95] tem aou oxygen phosphate phosphate_star 4.255961 3.033021 2.960664 JGOFS_GO 7.288981
Indo-Pacific (27.85,27.95] tem aou oxygen silicate phosphate 4.280774 3.000763 2.891621 JGOFS_GO 7.281537
Indo-Pacific (27.85,27.95] tem aou silicate phosphate phosphate_star 4.274004 2.989254 2.887257 JGOFS_GO 7.263258
Indo-Pacific (27.85,27.95] tem oxygen silicate phosphate phosphate_star 4.294333 2.994511 2.882504 JGOFS_GO 7.288844
Indo-Pacific (27.95,28.05] aou oxygen phosphate phosphate_star 4.378771 2.877804 3.077571 GO_new 5.955375
Indo-Pacific (27.95,28.05] aou oxygen silicate phosphate 4.381257 2.867418 3.113633 GO_new 5.981051
Indo-Pacific (27.95,28.05] aou oxygen silicate phosphate phosphate_star 4.377427 2.861323 3.076725 GO_new 5.938048
Indo-Pacific (27.95,28.05] sal aou oxygen phosphate phosphate_star 4.359852 2.877276 3.073963 GO_new 5.951239
Indo-Pacific (27.95,28.05] sal aou oxygen silicate phosphate 4.359139 2.857154 3.075517 GO_new 5.932671
Indo-Pacific (27.95,28.05] sal tem aou oxygen phosphate 4.359973 2.871354 3.107381 GO_new 5.978735
Indo-Pacific (27.95,28.05] sal tem aou silicate phosphate 4.377696 2.882277 3.064828 GO_new 5.947104
Indo-Pacific (27.95,28.05] tem aou oxygen phosphate phosphate_star 4.378207 2.870660 3.073851 GO_new 5.944512
Indo-Pacific (27.95,28.05] tem aou oxygen silicate phosphate 4.381237 2.866160 3.109301 GO_new 5.975461
Indo-Pacific (27.95,28.05] tem aou silicate phosphate phosphate_star 4.401881 2.886861 3.065987 GO_new 5.952847
Indo-Pacific (27.95,28.05] aou oxygen phosphate phosphate_star 4.378771 2.877804 3.077571 JGOFS_GO 7.256575
Indo-Pacific (27.95,28.05] aou oxygen silicate phosphate 4.381257 2.867418 3.113633 JGOFS_GO 7.248675
Indo-Pacific (27.95,28.05] aou oxygen silicate phosphate phosphate_star 4.377427 2.861323 3.076725 JGOFS_GO 7.238750
Indo-Pacific (27.95,28.05] sal aou oxygen phosphate 4.360181 2.879287 3.107563 JGOFS_GO 7.239469
Indo-Pacific (27.95,28.05] sal aou oxygen phosphate phosphate_star 4.359852 2.877276 3.073963 JGOFS_GO 7.237128
Indo-Pacific (27.95,28.05] sal aou oxygen silicate phosphate 4.359139 2.857154 3.075517 JGOFS_GO 7.216293
Indo-Pacific (27.95,28.05] sal tem aou oxygen phosphate 4.359973 2.871354 3.107381 JGOFS_GO 7.231328
Indo-Pacific (27.95,28.05] tem aou oxygen phosphate 4.383037 2.871954 3.114960 JGOFS_GO 7.254991
Indo-Pacific (27.95,28.05] tem aou oxygen phosphate phosphate_star 4.378207 2.870660 3.073851 JGOFS_GO 7.248867
Indo-Pacific (27.95,28.05] tem aou oxygen silicate phosphate 4.381237 2.866160 3.109301 JGOFS_GO 7.247397
Indo-Pacific (28.05,28.1] aou oxygen phosphate phosphate_star 4.232069 2.875188 2.849672 GO_new 5.724859
Indo-Pacific (28.05,28.1] aou oxygen silicate phosphate phosphate_star 4.227177 2.874648 2.825582 GO_new 5.700230
Indo-Pacific (28.05,28.1] sal aou oxygen phosphate 4.228221 2.863082 2.828633 GO_new 5.691715
Indo-Pacific (28.05,28.1] sal aou oxygen phosphate phosphate_star 4.228217 2.860046 2.799367 GO_new 5.659413
Indo-Pacific (28.05,28.1] sal aou oxygen silicate phosphate 4.225670 2.848200 2.826261 GO_new 5.674461
Indo-Pacific (28.05,28.1] sal aou silicate phosphate 4.362882 2.877440 2.836183 GO_new 5.713623
Indo-Pacific (28.05,28.1] sal aou silicate phosphate phosphate_star 4.362458 2.865816 2.804566 GO_new 5.670382
Indo-Pacific (28.05,28.1] sal tem aou oxygen phosphate 4.225054 2.860776 2.817179 GO_new 5.677955
Indo-Pacific (28.05,28.1] sal tem aou silicate phosphate 4.264970 2.856439 2.836178 GO_new 5.692618
Indo-Pacific (28.05,28.1] tem aou oxygen phosphate phosphate_star 4.227073 2.874723 2.843130 GO_new 5.717853
Indo-Pacific (28.05,28.1] aou oxygen phosphate phosphate_star 4.232069 2.875188 2.849672 JGOFS_GO 7.107257
Indo-Pacific (28.05,28.1] aou oxygen silicate phosphate 4.228213 2.881751 2.870188 JGOFS_GO 7.109964
Indo-Pacific (28.05,28.1] aou oxygen silicate phosphate phosphate_star 4.227177 2.874648 2.825582 JGOFS_GO 7.101825
Indo-Pacific (28.05,28.1] sal aou oxygen phosphate 4.228221 2.863082 2.828633 JGOFS_GO 7.091304
Indo-Pacific (28.05,28.1] sal aou oxygen phosphate phosphate_star 4.228217 2.860046 2.799367 JGOFS_GO 7.088264
Indo-Pacific (28.05,28.1] sal aou oxygen silicate phosphate 4.225670 2.848200 2.826261 JGOFS_GO 7.073870
Indo-Pacific (28.05,28.1] sal tem aou oxygen phosphate 4.225054 2.860776 2.817179 JGOFS_GO 7.085830
Indo-Pacific (28.05,28.1] tem aou oxygen phosphate 4.227670 2.881411 2.893756 JGOFS_GO 7.109081
Indo-Pacific (28.05,28.1] tem aou oxygen phosphate phosphate_star 4.227073 2.874723 2.843130 JGOFS_GO 7.101796
Indo-Pacific (28.05,28.1] tem aou oxygen silicate phosphate 4.226770 2.880685 2.860748 JGOFS_GO 7.107455
Indo-Pacific (28.1, Inf] aou oxygen silicate phosphate 3.640879 3.028695 2.696173 GO_new 5.724868
Indo-Pacific (28.1, Inf] aou oxygen silicate phosphate phosphate_star 3.635023 2.938152 2.688069 GO_new 5.626221
Indo-Pacific (28.1, Inf] sal aou oxygen phosphate 3.599961 3.230532 2.595467 GO_new 5.825999
Indo-Pacific (28.1, Inf] sal aou oxygen phosphate phosphate_star 3.599599 3.136083 2.593341 GO_new 5.729425
Indo-Pacific (28.1, Inf] sal aou oxygen silicate phosphate 3.570328 3.028688 2.582084 GO_new 5.610771
Indo-Pacific (28.1, Inf] sal aou phosphate phosphate_star 3.713417 3.207210 2.648134 GO_new 5.855343
Indo-Pacific (28.1, Inf] sal aou silicate phosphate phosphate_star 3.698036 3.187854 2.643218 GO_new 5.831072
Indo-Pacific (28.1, Inf] sal tem aou oxygen phosphate 3.540683 3.225881 2.586405 GO_new 5.812286
Indo-Pacific (28.1, Inf] sal tem aou phosphate phosphate_star 3.540570 3.173834 2.644631 GO_new 5.818464
Indo-Pacific (28.1, Inf] tem aou oxygen silicate phosphate 3.601769 2.985715 2.648408 GO_new 5.634123
Indo-Pacific (28.1, Inf] aou oxygen silicate phosphate 3.640879 3.028695 2.696173 JGOFS_GO 6.669574
Indo-Pacific (28.1, Inf] aou oxygen silicate phosphate phosphate_star 3.635023 2.938152 2.688069 JGOFS_GO 6.573176
Indo-Pacific (28.1, Inf] sal aou oxygen phosphate phosphate_star 3.599599 3.136083 2.593341 JGOFS_GO 6.735682
Indo-Pacific (28.1, Inf] sal aou oxygen silicate phosphate 3.570328 3.028688 2.582084 JGOFS_GO 6.599015
Indo-Pacific (28.1, Inf] sal tem aou oxygen phosphate 3.540683 3.225881 2.586405 JGOFS_GO 6.766564
Indo-Pacific (28.1, Inf] sal tem aou phosphate phosphate_star 3.540570 3.173834 2.644631 JGOFS_GO 6.714404
Indo-Pacific (28.1, Inf] sal tem aou silicate phosphate 3.542161 3.278460 2.644383 JGOFS_GO 6.820621
Indo-Pacific (28.1, Inf] sal tem oxygen phosphate phosphate_star 3.594199 3.222924 2.710942 JGOFS_GO 6.817124
Indo-Pacific (28.1, Inf] tem aou oxygen silicate phosphate 3.601769 2.985715 2.648408 JGOFS_GO 6.587484
Indo-Pacific (28.1, Inf] tem aou silicate phosphate phosphate_star 3.653818 3.124788 2.839805 JGOFS_GO 6.778606

2.3 Correlations RMSE between eras

2.3.1 All models

all_lm_stats %>%
  ggplot(aes(rmse_GO_SHIP, rmse_sum - rmse_GO_SHIP, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras ~ basin)

# rm(all_lm_stats, all_lm_stats_wide)

2.3.2 10 best models

lm_best %>%
  ggplot(aes(rmse_GO_SHIP, rmse_sum - rmse_GO_SHIP, col = gamma_slab)) +
  geom_point() +
  scale_color_viridis_d() +
  facet_grid(eras ~ basin)

#rm(all_lm_stats, all_lm_stats_wide)
lm_1 <- GLODAP %>% 
  filter(basin == "Atlantic",
         gamma_slab == "(-Inf,26]",
         era == "new_era") %>% 
  lm(Cstar ~ sal + tem + aou + oxygen + phosphate_star, data = .)

lm_2 <- GLODAP %>% 
  filter(basin == "Atlantic",
         gamma_slab == "(-Inf,26]",
         era == "new_era") %>% 
  lm(Cstar ~ sal + tem + aou + oxygen + phosphate, data = .)

lm_1
lm_2

sqrt(mean(lm_1$residuals^2))
sqrt(mean(lm_2$residuals^2))

2.4 Fitting best models

After selecting 10 linear regression models with lowest summed RMSE across all eras, models are fitted again and model coefficients are saved to file.

lm_best <- lm_best %>% 
  mutate(lm_coeff = str_replace_all(predictors, " ", " + "),
         lm_coeff = paste("Cstar ~", lm_coeff))

for (i_basin in unique(GLODAP$basin)) {
  for (i_era in unique(GLODAP$era)) {
    
    # i_basin <- unique(GLODAP$basin)[1]
    # i_era   <- unique(GLODAP$era)[1]
    
    print(i_basin) 
    print(i_era  )
    
    GLODAP_basin_era <- GLODAP %>%
      filter(basin == i_basin, era == i_era)

    lm_best_basin <- lm_best %>%
      filter(basin == i_basin)

    for (i_gamma_slab in unique(GLODAP_basin_era$gamma_slab)) {

      # i_gamma_slab <- unique(GLODAP_basin_era$gamma_slab)[1]
      print(i_gamma_slab)
      
      GLODAP_basin_era_slab <- GLODAP_basin_era %>% 
        filter(gamma_slab == i_gamma_slab)
      
      lm_best_basin_slab <- lm_best_basin %>% 
        filter(gamma_slab == i_gamma_slab)
      
      for (i_eras in unique(lm_best_basin_slab$eras)) {
        
        # i_eras <- unique(lm_best_basin_slab$eras)[1]
        
        lm_best_basin_slab_eras <- lm_best_basin_slab %>%
          filter(eras == i_eras)
      
      for (i_predictors in unique(lm_best_basin_slab_eras$predictors)) {
        
        # i_predictors <- unique(lm_best_basin_slab_eras$predictors)[1]
        
        lm_best_single <- lm_best_basin_slab_eras %>%
          filter(predictors == i_predictors) %>%
          select(lm_coeff) %>%
          pull()
        
        i_rmse_sum <- lm_best_basin_slab_eras %>%
          filter(predictors == i_predictors) %>%
          select(rmse_sum) %>%
          pull()
        
        lm_formula <- as.formula(lm_best_single)
        
        model_single <- lm(lm_formula, 
              data = GLODAP_basin_era_slab)
        
        coefficients <- tidy(model_single)
        coefficients <- coefficients %>% 
          mutate(basin = i_basin,
                 era = i_era,
                 eras = i_eras,
                 gamma_slab = i_gamma_slab,
                 model = lm_best_single,
                 rmse = i_rmse_sum)
        
        print(lm_best_single)
      
      if (exists("all_lm")) {
        all_lm <- bind_rows(all_lm, coefficients)
      }
      
      if (!exists("all_lm")) {
        all_lm <- coefficients
      }
      

      }
      
 
      }
    }
    
  }
  
}


rm(lm_formula, model_single, coefficients, i_rmse_sum,
   GLODAP_basin_era, GLODAP_basin_era_slab,
   lm_best, lm_best_basin, lm_best_basin_slab, lm_best_single,
   i_basin, i_era, i_gamma_slab, i_predictors)

all_lm %>%
  write_csv(here::here("data/eMLR",
                       "all_lm.csv"))


# Prepare model coefficients for mapping

all_lm <- all_lm %>% 
  select(term, estimate, basin, era, eras, gamma_slab, model)

all_lm <- all_lm %>% 
  mutate(estimate = if_else(is.na(estimate), 0, estimate))

all_lm_wide <- all_lm %>% 
  pivot_wider(names_from = era, values_from = estimate,
              names_prefix = "coeff_")

all_lm_wide <- all_lm_wide %>% 
  mutate(JGOFS_GO = coeff_GO_SHIP - coeff_JGOFS_WOCE,
         GO_new = coeff_new_era - coeff_GO_SHIP) %>% 
  select(-c(coeff_JGOFS_WOCE,
            coeff_GO_SHIP,
            coeff_new_era))

all_lm_long <- all_lm_wide %>% 
  pivot_longer(JGOFS_GO:GO_new, names_to = "eras_fit", values_to = "delta_coeff")

all_lm_long <- all_lm_long %>% 
  filter(eras == eras_fit) %>% 
  select(-eras_fit)

all_lm_wide <- all_lm_long %>% 
  pivot_wider(values_from = delta_coeff,
              names_from = term,
              names_prefix = "delta_coeff_",
              values_fill = 0)

all_lm_wide %>%
  write_csv(here::here("data/eMLR",
                       "all_lm_wide.csv"))

rm(all_lm_long, all_lm_wide)

3 Open tasks

  • combine correlation, rmse and model calculation into one routine, ideally using the purrr map approach
  • select exactly best 10 models, even if same rmse sum appears more than once (currently, selection is set two ten, even if doublicated rmse occur)
  • check (phosphate_star) NA’s in coeffcients -> cases of exact colinearity?
  • test model selection based on AIC

4 Open questions


sessionInfo()
R version 4.0.2 (2020-06-22)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 18363)

Matrix products: default

locale:
[1] LC_COLLATE=English_Germany.1252  LC_CTYPE=English_Germany.1252   
[3] LC_MONETARY=English_Germany.1252 LC_NUMERIC=C                    
[5] LC_TIME=English_Germany.1252    

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

other attached packages:
 [1] corrr_0.4.2      broom_0.7.0      kableExtra_1.1.0 knitr_1.29      
 [5] olsrr_0.5.3.9000 GGally_2.0.0     patchwork_1.0.1  lubridate_1.7.9 
 [9] forcats_0.5.0    stringr_1.4.0    dplyr_1.0.0      purrr_0.3.4     
[13] readr_1.3.1      tidyr_1.1.0      tibble_3.0.3     ggplot2_3.3.2   
[17] tidyverse_1.3.0  workflowr_1.6.2 

loaded via a namespace (and not attached):
 [1] fs_1.4.2           webshot_0.5.2      RColorBrewer_1.1-2 httr_1.4.2        
 [5] rprojroot_1.3-2    tools_4.0.2        backports_1.1.8    R6_2.4.1          
 [9] nortest_1.0-4      DBI_1.1.0          colorspace_1.4-1   withr_2.2.0       
[13] tidyselect_1.1.0   gridExtra_2.3      curl_4.3           compiler_4.0.2    
[17] git2r_0.27.1       cli_2.0.2          rvest_0.3.6        xml2_1.3.2        
[21] labeling_0.3       scales_1.1.1       goftest_1.2-2      digest_0.6.25     
[25] foreign_0.8-80     rmarkdown_2.3      rio_0.5.16         pkgconfig_2.0.3   
[29] htmltools_0.5.0    dbplyr_1.4.4       highr_0.8          rlang_0.4.7       
[33] readxl_1.3.1       rstudioapi_0.11    generics_0.0.2     farver_2.0.3      
[37] jsonlite_1.7.0     zip_2.0.4          car_3.0-8          magrittr_1.5      
[41] Rcpp_1.0.5         munsell_0.5.0      fansi_0.4.1        abind_1.4-5       
[45] lifecycle_0.2.0    stringi_1.4.6      whisker_0.4        yaml_2.2.1        
[49] carData_3.0-4      plyr_1.8.6         grid_4.0.2         blob_1.2.1        
[53] promises_1.1.1     crayon_1.3.4       haven_2.3.1        hms_0.5.3         
[57] pillar_1.4.6       reprex_0.3.0       glue_1.4.1         evaluate_0.14     
[61] data.table_1.13.0  modelr_0.1.8       vctrs_0.3.2        httpuv_1.5.4      
[65] cellranger_1.1.0   gtable_0.3.0       reshape_0.8.8      assertthat_0.2.1  
[69] xfun_0.16          openxlsx_4.1.5     later_1.1.0.1      viridisLite_0.3.0 
[73] ellipsis_0.3.1     here_0.1