BGC-Argo Temperature Data
Pasqualina Vonlanthen & Jens Daniel Müller
29 April, 2022
Last updated: 2022-04-29
Checks: 7 0
Knit directory: bgc_argo_r_argodata/
This reproducible R Markdown analysis was created with workflowr (version 1.7.0). 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(20211008) 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 8b582f0. 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: .RData
Ignored: .Rhistory
Ignored: .Rproj.user/
Ignored: output/
Untracked files:
Untracked: code/OceanSODA_argo_extremes.R
Untracked: code/creating_dataframe.R
Untracked: code/creating_map.R
Untracked: code/merging_oceanSODA_Argo.R
Untracked: code/pH_data_timeseries.R
Unstaged changes:
Modified: analysis/_site.yml
Deleted: analysis/load_OceanSODA_biomes.Rmd
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/temp_data.Rmd) and HTML (docs/temp_data.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 | e61c08e | pasqualina-vonlanthendinenna | 2022-04-27 | Build site. |
| Rmd | 9664e0e | pasqualina-vonlanthendinenna | 2022-04-27 | added temp data page, changed double extremes |
Task
Explore BGC-Argo temperature data through timeseries and monthly climatological maps
path_argo <- '/nfs/kryo/work/updata/bgc_argo_r_argodata'
path_emlr_utilities <- "/nfs/kryo/work/jenmueller/emlr_cant/utilities/files/"
path_basin_mask <- "/nfs/kryo/work/updata/reccap2/"
path_argo_preprocessed <- paste0(path_argo, "/preprocessed_bgc_data")Load BGC-SST data
Using only temperature data from flag A profiles which have an associated flag A pH profile
# flag A temperature data in the top 20 m
sst <- read_rds(file = paste0(path_argo_preprocessed, "/bgc_merge_flag_A.rds")) %>%
filter(between(depth, 0, 20)) %>%
mutate(year = year(date),
month = month(date)) %>%
select(-c(ph_in_situ_total_adjusted, ph_in_situ_total_adjusted_error, ph_in_situ_total_adjusted_qc, profile_ph_in_situ_total_qc))
# load in biome separations
nm_biomes <- read_rds(file = paste0(path_argo_preprocessed, "/nm_biomes.rds"))Southern Ocean SST
sst_SO <- sst %>%
filter(lat <= -30)SST offset with depth
Difference between the in-situ measured sst (20 m) and the profile-mean 20m temperature
# calculate the mean sst for each surface profile
mean_profile_sst <- sst_SO %>%
group_by(platform_number, cycle_number) %>%
mutate(mean_prof_sst = mean(temp_adjusted, na.rm = TRUE),
.before = depth) %>%
ungroup() %>%
mutate(offset = temp_adjusted-mean_prof_sst,
.after = mean_prof_sst) # subtract the mean profile sst from the measured in situ sst
mean_profile_sst %>%
ggplot()+
geom_point(aes(x = offset, y = depth, col = as.character(year)), size = 0.3, pch = 19) +
scale_y_reverse()+
geom_vline(xintercept = 0, col = 'red', size = 0.6)+
labs(x = 'offset (ºC)',
y = 'depth (m)',
col = 'year',
title = 'in situ sst - mean profile sst')
| Version | Author | Date |
|---|---|---|
| e61c08e | pasqualina-vonlanthendinenna | 2022-04-27 |
Bin the sst data into 2m-depth intervals and calculate the offset for each sst observation in each depth interval relative to the profile-mean sst
# bin the sst values into 2m bins and calculate the offset for each 2m bin
mean_profile_sst_binned <- sst_SO %>%
mutate(depth = cut(depth, seq(0, 20, 2), seq(1, 19, 2)),
depth = as.numeric(as.character(depth))) %>%
group_by(platform_number, cycle_number) %>%
mutate(mean_prof_sst = mean(temp_adjusted, na.rm = TRUE),
.before = depth) %>%
ungroup() %>%
mutate(offset = temp_adjusted-mean_prof_sst,
.after = mean_prof_sst)
# plot the offset of the depth-binned values
mean_profile_sst_binned %>%
ggplot()+
geom_point(aes(x = offset, y = depth, col = as.character(year)), size = 0.3, pch = 19) +
scale_y_reverse()+
geom_vline(xintercept = 0, col = 'red', size = 0.6)+
labs(x = 'offset (ºC)',
y = 'depth (m)',
col = 'year',
title = 'in situ sst - mean profile sst (2m depth bins)')
| Version | Author | Date |
|---|---|---|
| e61c08e | pasqualina-vonlanthendinenna | 2022-04-27 |
Mean binned offset
# bin the ph values into 2m bins and calculate the offset for each 2m bin
profile_sst_binned_ave <- sst_SO %>%
mutate(depth = cut(depth, seq(0, 20, 2), seq(1, 19, 2)),
depth = as.numeric(as.character(depth))) %>%
group_by(platform_number, cycle_number) %>%
mutate(mean_prof_sst = mean(temp_adjusted, na.rm = TRUE),
.before = depth) %>%
ungroup() %>%
mutate(offset = temp_adjusted-mean_prof_sst,
.after = mean_prof_sst) %>%
group_by(depth) %>%
summarise(mean_offset = mean(offset))
# plot the offset of the depth-binned values
profile_sst_binned_ave %>%
ggplot()+
geom_point(aes(x = mean_offset, y = depth), size = 1, pch = 19) +
geom_line(aes(x = mean_offset, y = depth))+
scale_y_reverse()+
geom_vline(xintercept = 0, col = 'red', size = 1)+
labs(x = 'mean offset (ºC)',
y = 'depth (m)',
col = 'year',
title = 'in situ sst - mean profile sst (2m depth bins)')
| Version | Author | Date |
|---|---|---|
| e61c08e | pasqualina-vonlanthendinenna | 2022-04-27 |
Monthly climatological sst
Map of monthly climatological Argo temperature (BGC floats, flag A pH profiles only), April 2014-December 2021
# average pH values in the top 20 m for each month in each 2 x 2º longitude/latitude grid
sst_clim_SO <- sst_SO %>%
group_by(lat, lon, month) %>%
summarise(sst_clim_month = mean(temp_adjusted))
# read in the map from updata
map <-
read_rds(paste(path_emlr_utilities,
"map_landmask_WOA18.rds",
sep = ""))
# map a monthly climatology of pH (April 2014 - August 2021)
map +
geom_tile(data = sst_clim_SO,
aes(lon, lat, fill = sst_clim_month)) +
lims(y = c(-85, -25)) +
scale_fill_viridis_c() +
labs(x = 'lon',
y = 'lat',
fill = 'SST',
title = 'Monthly climatological \nArgo SST (Apr 2014 - Dec 2021)') +
theme(legend.position = 'right') +
facet_wrap(~month, ncol = 2)
| Version | Author | Date |
|---|---|---|
| e61c08e | pasqualina-vonlanthendinenna | 2022-04-27 |
basemap(limits = -32, data = sst_clim_SO) + # change to polar projection
geom_spatial_tile(data = sst_clim_SO,
aes(x = lon,
y = lat,
fill = sst_clim_month),
linejoin = 'mitre',
col = 'transparent',
detail = 60)+
scale_fill_viridis_c()+
theme(legend.position = 'bottom')+
labs(x = 'lon',
y = 'lat',
fill = 'SST',
title = 'monthly climatological \nArgo SST (Apr 2014 - Dec 2021)')+
facet_wrap(~month, ncol = 2)Monthly timeseries
Timeseries of monthly SST values, for each Mayot biome
# plot the region separations on a map
map +
geom_raster(data = nm_biomes,
aes(x = lon,
y = lat,
fill = biome_name)) +
labs(title = 'Southern Ocean Mayot biomes',
fill = 'biome')
| Version | Author | Date |
|---|---|---|
| e61c08e | pasqualina-vonlanthendinenna | 2022-04-27 |
# plot a timeseries of monthly values over the whole southern ocean south of 30ºS
sst_SO <- inner_join(sst_SO, nm_biomes)
sst_month_SO <- sst_SO %>%
group_by(year, month, biome_name) %>%
summarise(sst_ave = mean(temp_adjusted, na.rm = TRUE))
# timeseries of monthly pH values over 2014-2021 (separate panels for each month)
sst_month_SO %>%
ggplot(aes(x = year,
y = sst_ave,
group = biome_name,
col = biome_name)) +
facet_wrap(~month) +
geom_line() +
geom_point() +
labs(x = 'year',
y = 'SST (ºC)',
title = 'monthly mean Argo SST (Apr 2014-Dec 2021, Southern Ocean)',
col = 'region')
| Version | Author | Date |
|---|---|---|
| e61c08e | pasqualina-vonlanthendinenna | 2022-04-27 |
Monthly average Southern Ocean SST, for each biome
# timeseries of monthly sst values for each year (separate years on the same plot)
sst_month_SO %>%
# filter(year != 2014) %>% # remove the year that is missing data
ggplot(aes(x = month,
y = sst_ave,
group = year,
col = as.character(year)))+
geom_line()+
geom_point()+
scale_x_continuous(breaks = seq(1, 12, 2))+
facet_wrap(~biome_name)+
labs(x = 'month',
y = 'SST (ºC)',
title = 'monthly mean Argo SST (Apr 2014-Dec 2021, Southern Ocean regions)',
col = 'year')
| Version | Author | Date |
|---|---|---|
| e61c08e | pasqualina-vonlanthendinenna | 2022-04-27 |
# calculate a yearly average SST (one SST value per year, for the whole biome)
sst_year_SO <- sst_SO %>%
group_by(year, biome_name) %>%
summarise(sst_ave = mean(temp_adjusted, na.rm = TRUE))
# plot a timeseries of the yearly average SST value (one value per year)
sst_year_SO %>%
ggplot(aes(x = year, y = sst_ave, group = biome_name, col = biome_name))+
geom_line()+
geom_point()+
labs(x = 'year',
y = 'SST (ºC)',
title = 'yearly mean Argo SST (Apr 2014-Dec 2021, south of 30ºS)',
col = 'region')
sessionInfo()R version 4.1.2 (2021-11-01)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: openSUSE Leap 15.3
Matrix products: default
BLAS: /usr/local/R-4.1.2/lib64/R/lib/libRblas.so
LAPACK: /usr/local/R-4.1.2/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] ggOceanMaps_1.2.6 ggspatial_1.1.5 oce_1.5-0 gsw_1.0-6
[5] lubridate_1.8.0 forcats_0.5.1 stringr_1.4.0 dplyr_1.0.7
[9] purrr_0.3.4 readr_2.1.1 tidyr_1.1.4 tibble_3.1.6
[13] ggplot2_3.3.5 tidyverse_1.3.1 workflowr_1.7.0
loaded via a namespace (and not attached):
[1] fs_1.5.2 sf_1.0-5 httr_1.4.2 rprojroot_2.0.2
[5] tools_4.1.2 backports_1.4.1 bslib_0.3.1 rgdal_1.5-28
[9] utf8_1.2.2 R6_2.5.1 KernSmooth_2.23-20 rgeos_0.5-9
[13] DBI_1.1.2 colorspace_2.0-2 raster_3.5-11 withr_2.4.3
[17] sp_1.4-6 tidyselect_1.1.1 processx_3.5.2 compiler_4.1.2
[21] git2r_0.29.0 cli_3.1.1 rvest_1.0.2 xml2_1.3.3
[25] labeling_0.4.2 sass_0.4.0 scales_1.1.1 classInt_0.4-3
[29] callr_3.7.0 proxy_0.4-26 digest_0.6.29 rmarkdown_2.11
[33] pkgconfig_2.0.3 htmltools_0.5.2 highr_0.9 dbplyr_2.1.1
[37] fastmap_1.1.0 rlang_1.0.2 readxl_1.3.1 rstudioapi_0.13
[41] farver_2.1.0 jquerylib_0.1.4 generics_0.1.1 jsonlite_1.7.3
[45] magrittr_2.0.1 Rcpp_1.0.8 munsell_0.5.0 fansi_1.0.2
[49] lifecycle_1.0.1 terra_1.5-12 stringi_1.7.6 whisker_0.4
[53] yaml_2.2.1 grid_4.1.2 parallel_4.1.2 promises_1.2.0.1
[57] crayon_1.4.2 lattice_0.20-45 haven_2.4.3 hms_1.1.1
[61] knitr_1.37 ps_1.6.0 pillar_1.6.4 codetools_0.2-18
[65] reprex_2.0.1 glue_1.6.0 evaluate_0.14 getPass_0.2-2
[69] modelr_0.1.8 vctrs_0.3.8 tzdb_0.2.0 httpuv_1.6.5
[73] cellranger_1.1.0 gtable_0.3.0 assertthat_0.2.1 xfun_0.29
[77] broom_0.7.11 e1071_1.7-9 later_1.3.0 viridisLite_0.4.0
[81] class_7.3-20 units_0.7-2 ellipsis_0.3.2