# Mauna Loa atmospheric CO2 change # multi language plot for Wikipedia # using different possible sources (Scripps if NOAA is not available) #  # Used for https://commons.wikimedia.org/wiki/File:Mauna_Loa_CO2_monthly_mean_concentration.svg # and related files in other languages. See all files in : # https://commons.wikimedia.org/wiki/Template:Other_versions/Mauna_Loa_CO2_monthly_mean_concentration #  # Update this script on : # https://commons.wikimedia.org/w/index.php?title=Template:Other_versions/Mauna_Loa_CO2_monthly_mean_concentration.R&action=edit  # Required packages.   used for : library(readr)         # reading data library(dplyr)         # data processing library(ggplot2)       # plot library(gridExtra)     # plot positionning library(scales)        # pretty breaks on plot axis library(lubridate)     # date processing library(httr)          # check url library(Hmisc)         # ci from bootstrap  # install all locales beforehand, e.g. : # sudo locale-gen ca_ES.UTF-8  old_locale <- Sys.getlocale("LC_TIME")  # Translations ------------------------------------------------------------  language <- list(      ca_ES = list(     locale_lc_time = "ca_ES.UTF-8",     title = bquote("Mitjana mensual de la concentratió de"~CO[2]),     subtitle = "Mauna Loa",     caption_scripps = paste("dades : R. F. Keeling, S. J. Walker, S. C. Piper and A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Visitada", Sys.Date()),     caption_noaa = paste("dades : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) i\n Dr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Visitada ", Sys.Date()),     x = "any",     y = bquote("fracció de"~CO[2]~"en aire sec ("*mu*"mol/mol)"),     x2 = "mes",     y2 = bquote(atop("fracció de"~CO[2]~"en aire sec ("*mu*"mol/mol)",                      "desviació de la mitjana anual")),     title2 = "Variació mensual"   ),      cs_CZ = list(     locale_lc_time = "cs_CZ.UTF-8",     title = "Průměrné měsíční koncentrace oxidu uhličitého",     subtitle = "Mauna Loa",     caption_scripps = paste("data : R. F. Keeling, S. J. Walker, S. C. Piper a A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Přístupné", Sys.Date()),     caption_noaa = paste("data : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) a\nDr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Přístupné", Sys.Date()),     x = "rok",     y = bquote("koncentrace"~CO[2]~"v suchém vzduchu ("*mu*"mol/mol)"),     x2 = "měsíc",     y2 = bquote(atop("koncentrace"~CO[2]~"v suchém vzduchu ("*mu*"mol/mol)", "odchylka od ročního průměru")),     title2 = "Měsíční změna (průměrná roční odchylka)"   ),      de_DE = list(     locale_lc_time = "de_DE.UTF-8",     title = bquote("Monatliche durchschnittliche"~CO[2]*"-Konzentration"),     subtitle = "Mauna Loa",     caption_scripps = paste("Datei : R. F. Keeling, S. J. Walker, S. C. Piper und A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Zugänglich am", Sys.Date()),     caption_noaa = paste("Datei : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) und\nDr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Zugänglich am", Sys.Date()),     x = "Jahr",     y = bquote(CO[2]*"-Anteil in trockener Luft ("*mu*"mol/mol)"),     x2 = "Monate",     y2 = bquote(atop(CO[2]*"-Anteil in trockener Luft ("*mu*"mol/mol)", "Abweichung vom Jahresmittel")),     title2 = "Monatliche Variation"   ),        en_US = list(     locale_lc_time = "en_US.UTF-8",     title = bquote("Monthly mean"~CO[2]~"concentration"),     subtitle = "Mauna Loa",     caption_scripps = paste("Data : R. F. Keeling, S. J. Walker, S. C. Piper and A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Accessed", Sys.Date()),     caption_noaa = paste("Data : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) and\n Dr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Accessed ", Sys.Date()),     x = "Year",     y = bquote(CO[2]~"fraction in dry air ("*mu*"mol/mol)"),     x2 = "Month",     y2 = bquote(atop(CO[2]~"fraction in dry air ("*mu*"mol/mol)", "Departure from yearly average")),     title2 = "Seasonal variation"   ),      es_ES = list(     locale_lc_time = "es_ES.UTF-8",     title = bquote("Media mensual de la concentración de"~CO[2]),     subtitle = "Mauna Loa",     caption_scripps = paste("dato : R. F. Keeling, S. J. Walker, S. C. Piper y A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Visitada", Sys.Date()),     caption_noaa = paste("dato : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) y\nDr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Visitada", Sys.Date()),     x = "Año",     y = bquote("Fraccion de"~CO[2]~"en aire secco ("*mu*"mol/mol)"),     x2 = "Mes",     y2 = bquote(atop("Fraccion de"~CO[2]~"en aire secco ("*mu*"mol/mol)", "Desviación de la media anual")),     title2 = "Variación mensual"   ),        fr_FR = list(     locale_lc_time = "fr_FR.UTF-8",     title = bquote("Moyenne mensuelle de la concentration de"~CO[2]),     subtitle = "Mauna Loa",     caption_scripps = paste("données : R. F. Keeling, S. J. Walker, S. C. Piper et A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Accédé le", Sys.Date()),     caption_noaa = paste("données : Dr Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) et\n Dr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Accédé le", Sys.Date()),     x = "année",     y = bquote("fraction de"~CO[2]~"dans l'air sec ("*mu*"mol/mol)"),     x2 = "mois",     y2 = bquote(atop("fraction de"~CO[2]~"dans l'air sec ("*mu*"mol/mol)", "en écart à la moyenne annuelle")),     title2 = "Variation saisonnière"   ),      nl_NL = list(     locale_lc_time = "nl_NL.UTF-8",     title = bquote("Maandelijkse gemiddelde"~CO[2]*"-concentratie"),     subtitle = "Mauna Loa",     caption_scripps = glue("Gegevens : R. F. Keeling, S. J. Walker, S. C. Piper en A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Geopend {Sys.Date()}"),     caption_noaa = glue("Gegevens : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) en\n Dr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Accessed  {Sys.Date()}"),     x = "Jaar",     y = bquote(CO[2]*"-fractie in droge lucht ("*mu*"mol/mol)"),     x2 = "Mannd",     y2 = bquote(atop(CO[2]*"-fractie in droge lucht ("*mu*"mol/mol)", "Afwijking van jaargemiddelde")),     title2 = "Seizoensgebonden variatie"   ),    nn_NO = list(     locale_lc_time = "nn_NO.UTF-8",     title = bquote("Gjennomsnittlig månedlig"~CO[2]*"-konsentrasjon"),     subtitle = "Mauna Loa",     caption_scripps = paste("data : R. F. Keeling, S. J. Walker, S. C. Piper og A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Vist", Sys.Date()),     caption_noaa = paste("data : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) og\nDr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Vist", Sys.Date()),     x = "År",     y = bquote(CO[2]*"-andel i tørr luft ("*mu*"mol/mol)"),     x2 = "Måned",     y2 = bquote(atop(CO[2]*"-andel i tørr luft ("*mu*"mol/mol)",                      "Avvik fra årlig gjennomsnitt")),     title2 = "Årlig variasjon"   ),      pl_PL = list(     locale_lc_time = "pl_PL.UTF-8",     title = bquote("Średnie miesięczne stężenie"~CO[2]),     subtitle = "Mauna Loa",     caption_scripps = paste("Dane : R. F. Keeling, S. J. Walker, S. C. Piper i A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Dostęp na", Sys.Date()),     caption_noaa = paste("Dane : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) i\nDr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Dostęp na", Sys.Date()),     x = "Rok",     y = bquote("Frakcja"~CO[2]~"w suchym powietrzu ("*mu*"mol/mol)"),     x2 = "Miesiąc",     y2 = bquote(atop("Frakcja"~CO[2]~"w suchym powietrzu ("*mu*"mol/mol)", "Odejście od średniej rocznej")),     title2 = "Wahania sezonowe"   ),    ro_RO = list(     locale_lc_time = "ro_RO.UTF-8",     title = bquote("Concentrația medie lunară de"~CO[2]),     subtitle = "Mauna Loa, Hawaii",     caption_scripps = paste("Date: R. F. Keeling, S. J. Walker, S. C. Piper și A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). Accesat în", Sys.Date()),     caption_noaa = paste("Date : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) și\n Dr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). Accesat în", Sys.Date()),     x = "An",     y = bquote("Ponderea"~CO[2]~"în aer uscat ("*mu*"mol/mol)"),     x2 = "Lună",     y2 = bquote(atop("Abaterea cotei"~CO[2]~"în aer uscat ("*mu*"mol/mol)", "Abatere de la media anuală")),     title2 = "Variație sezonieră"   ),    ru_RU = list(     locale_lc_time = "ru_RU.UTF-8",     title = bquote("Среднемесячная и среднегодовая концентрации"~CO[2]),     subtitle = "Мауна Лоа, Гаваи",     caption_scripps = paste("данные : R. F. Keeling, S. J. Walker, S. C. Piper и A. F. Bollenbacher\nScripps CO2 Program (https://scrippsco2.ucsd.edu/). доступ на", Sys.Date()),     caption_noaa = paste("данные : Dr. Pieter Tans, NOAA/ESRL (https://gml.noaa.gov/ccgg/trends/) и\nDr. Ralph Keeling, Scripps Institution of Oceanography (https://scrippsco2.ucsd.edu/). доступ на ", Sys.Date()),     x = "Годы",     y = bquote("Доля"~CO[2]~"в сухом воздухе ("*mu*"mol/mol)"),     x2 = "Месяцы",     y2 = bquote(atop("Отклонение доли"~CO[2]~"в сухом воздухе", "от среднегодового значения, мкмоль/моль")),     title2 = "Сезонные колебания"   )    )   # Data --------------------------------------------------------------------  # First try NOAA data, if not availble use Scripps data if (!http_error("ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_mm_mlo.txt")) {   # https://gml.noaa.gov/ccgg/trends/data.html   source <- "noaa"   co2ml <- read_table(     "ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_mm_mlo.txt",     col_names = c("year", "month", "decimal", "average", "deseason", "ndays", "stddays", "unc"),     col_types = "iidddidd",     na = c("-99.99", "-1"),     comment = "#") %>%      group_by(year) %>%      mutate(year_mean = mean(average, na.rm = TRUE),            delta = average - year_mean,            vdate = ymd(paste0("2015-", month, "-01")))%>%      rename(co2 = average)    } else {   # https://scrippsco2.ucsd.edu/data/atmospheric_co2/primary_mlo_co2_record   # used during US gov shutdown   source <- "scripps"   co2ml <- read_csv(     "https://scrippsco2.ucsd.edu/assets/data/atmospheric/stations/in_situ_co2/monthly/monthly_in_situ_co2_mlo.csv",     col_names = c("year", "month", "xls_date", "decimal",                   "co2_raw", "co2_raw_seas_adj", "fit", "fit_seas_adj",                   "co2_filled", "co2_filled_seas_adj"),     col_types = "iiiddddddd",     skip = 57,     na = "-99.99",     comment = "\"") %>%      group_by(year) %>%      mutate(year_mean = mean(co2_filled, na.rm = TRUE),            delta = co2_filled - year_mean,            vdate = ymd(paste0("2015-", month, "-01"))) %>%      rename(co2 = co2_filled) }   # Generate the plot for each language -------------------------------------  for (l in names(language)) {   message(l)   current <- language[[l]]      # format the date in local names   Sys.setlocale("LC_TIME", current$locale_lc_time)      # main plot   p1 <- ggplot(co2ml, aes(decimal, co2)) + 	     geom_line(color = "pink") +     geom_point(color = "red", size = 0.6) +     stat_smooth(span = 0.1) +     scale_x_continuous(breaks = pretty_breaks()) +     scale_y_continuous(breaks = pretty_breaks()) +     labs(x = current$x,          y = current$y,          title = current$title,          subtitle = paste(current$subtitle, min(co2ml$year), "-", max(co2ml$year)),          caption = paste(current[[paste0("caption_", source)]],                          "https://w.wiki/4ZWn", sep = "\n")) +     theme_bw() +     theme(plot.caption = element_text(size = 7))      # inset plot   p2 <- ggplot(co2ml, aes(vdate, delta)) +     geom_hline(yintercept = 0) +     stat_smooth(span = 0.4, se = FALSE) +     stat_summary(fun.data = "mean_cl_boot", colour = "red", size = 0.3) +      scale_x_date(breaks = pretty_breaks(4),                   minor_breaks = pretty_breaks(12),                   labels = date_format("%b")) +     labs(x = current$x2,          y = current$y2,          title = current$title2) +     theme_bw()      # merge the plots and export in SVG   p1 + annotation_custom(grob = ggplotGrob(p2),                          xmin = 1957, xmax = 1991, ymin = 361, ymax = 420)      ggsave(file = paste("co2_mauna_loa", l, Sys.Date(), "wp.svg", sep = "_"),           width = 20,           height = 20,           units = "cm",           device = svg) }  # return to former locale Sys.setlocale("LC_TIME", old_locale)