Édition du 14 août 2010 : à lire aussi mon billet sur le livre Le mythe climatique de Benoît Rittaud, qui nuance quelque peu les données présentées dans billet.

MétéoSuisse propose au téléchargement des données climatiques homogénéisées pour 12 stations de mesures en Suisse. Une bonne occasion pour jeter un œil au réchauffement…

Commençons par télécharger les données pour Genève. Notons que les données démarrent à la ligne 28 de ce fichier : il faudra donc en sauter 27 pour les lire dans R.

clim.ge <- read.table("geneve.txt", header=TRUE, skip=27)

Affichons maintenant les températures :

plot(clim.ge$Temperature, type="l")

On ne voit pas grand chose, si ce n'est des variations sur l'année, qui restent évidemment beaucoup plus importantes qu'un éventuel réchauffement. Nous allons donc moyenner les températures annuelles pour y voir plus clair. Il faudra évidemment supprimer 2010 qui est incomplète :

years <- cut(clim.ge$Year, 1863:2009)

temp.anu <- tapply(clim.ge$Temperature, years, mean)

plot(1864:2009, temp.anu, type="l")

C'est tout de suite beaucoup plus clair. Si les moyennes anuelles oscillaient entre 8 et 10 degrés jusque vers 1940, depuis les années 1990 environ, c'est plutôt entre 10 et 11 !

Regardons maintenant les moyennes sur 10 ans. Pour cela, nous allons devoir couper les années 1864 à 1869 pour regarder de 1870 à 2009, soit 14 décennies :

decades <- cut(clim.ge.hist.dec$Year, seq(1869, 2009, by=10))

temp.dec <- with(clim.ge.hist.dec, tapply(Temperature, decades, mean))

plot(seq(1869, 2009, by=10)[-1], temp.dec, type="l")

Là c'est évident, après un premier pic dans les années 1940, suivies par un tassement jusqu'aux années 1970, les trois dernières décennies ont été les plus chaudes de la période observée.

Bien sûr, cela ne donne aucune indication sur les causes ni les conséquences possibles, et rien ne dit non plus que la tendance va se poursuivre ; mais le réchauffement est bel et bien là, la température a grimpé de 2 degrés en moins d'un siècle.

Pour finir, quelques diagrammes ombrothermique sur plusieurs intervalles de 30 ans, basés sur la fonction diagwl du package CRAN climatol. J'ai modifié la fonction en particulier pour traiter un matrice 2x12 contenant les pluies et températures moyennes plutôt que 4x12 contenant les températures maximales, minimales et minimales absolues, données non communiquées par MétéoSuisse.

source("diagwl.R") # Téléchargez la fonction modifiée avec ce lien
climat <- function(clim, ...) {
  temp.month <- with(clim, tapply(Temperature, Month, mean))
  plui.month <- with(clim, tapply(Precipitation, Month, mean))
  clim.dat <- rbind(plui.month, temp.month, temp.month, rep(NA, 12))
  diagwl(clim.dat, ...)
}
climat(subset(clim.ge, Year < 1894), est="Genève-Cointrin", alt=420, per="1864-1893", mlab="en")
climat(subset(clim.ge, Year >= 1931 & Year <= 1960), est="Genève-Cointrin", alt=420, per="1931-1960", mlab="en")
climat(subset(clim.ge, Year >= 1961 & Year <= 1990), est="Genève-Cointrin", alt=420, per="1961-1990", mlab="en")
climat(subset(clim.ge, Year >= 1980 & Year <= 2009), est="Genève-Cointrin", alt=420, per="1980-2009", mlab="en")