R

# 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(linnames(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)