Acidification
R-bloggers 2025-11-21
[This article was first published on r.iresmi.net, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Day 20 of 30DayMapChallenge: « Water » (previously).
Global ocean acidification mean sea water pH trend map from Multi-Observations Reprocessing from Copernicus.
Data
library(sf)library(ggplot2)library(rnaturalearth)library(glue)library(terra)library(ggspatial)eqearth <- "EPSG:8857"world <- ne_countries() |> st_transform(eqearth)mask <- c(xmin = -179, ymin = -89, xmax = 179, ymax = 89) |> st_bbox() |> st_as_sfc() |> st_set_crs("EPSG:4326") |> st_sf() |> st_segmentize(100) |> st_transform(eqearth) acid_trend <- "global_omi_health_carbon_ph_trend_1985_P20230930.nc" |> rast() |> rotate() |> project(eqearth) |> mask(mask)Map
world |> ggplot() + layer_spatial(data = acid_trend, aes(fill = after_stat(band1))) + geom_sf(data = mask) + geom_sf(color = "grey", fill = "white") + scale_fill_viridis_c(name = bquote(Delta*pH~yr^-1), direction = -1, na.value = "white") + labs(title = "Global ocean acidification", subtitle = "mean sea water pH trend", caption = glue("data: Copernicus / LSCE doi:10.48670/moi-00277 Natural Earth - {st_crs(eqearth)$Name} https://r.iresmi.net - {Sys.Date()}")) + theme_void() + theme(plot.caption = element_text(size = 7, color = "grey40"), plot.margin = unit(c(.2, .2, .2, .2), units = "cm"), legend.position = "bottom", legend.text = element_text(angle = 45, vjust = 1, hjust = 1))
To leave a comment for the author, please follow the link and comment on their blog: r.iresmi.net.
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.