Der Jahrhundertsommer war nicht besonders trocken

Während es in anderen Teilen Deutschlands zu langanhaltender Dürre kam, war die Mange des Niederschlags zumindest in München durchschnittlich. Ein Grund sind sicherlich die regelmäßigen Sommergewitter. In einem Stufendiegramm lässt sich der Verlauf gut visualisieren.
Ein Stufendiagramm des Niederschlags eines bestimmten Ortes lässt sich mit ggplot2 erzeugen. Die größte Schwierigkeit liegt darin, die passenden Daten zu finden. Hier hilft der Deutsche Wetterdienst (www.dwd.de). Unter ftp://ftp-cdc.dwd.de/pub/CDC/observations_germany/climate/daily/kl/historical/ finden sich die Messwerte von 1086 Wetterstationen quer durch Deutschland. München hat die Nummer 3379. Die Daten müssen etwas nachbearbeitet werden, da manche Messwerte doppelt oder mit falscher Jahreszahl vorliegen.

library(ggplot2)
library(dplyr)
library(lubridate)

rain <- read.csv("E:/Source/R-3.3.1/Projekte/daten/precipitation.txt", sep = ";", header = T, row.names = NULL)
rain[rain == -999] = NA
rain$MESS_DATUM <- year(ymd(rain$MESS_DATUM))
rain <- rain %>%
  select(MESS_DATUM, RSK) %>%
  group_by(MESS_DATUM) %>%
  summarise(regen = sum(RSK))
av <- round(mean(rain$regen))
ggplot(rain)+
  geom_step(aes(x = MESS_DATUM, y = regen), colour = "blue")+
  theme_bw()+
  theme(axis.title.x = element_blank(), plot.title = element_text(size=12))+
  labs(title="MÜNCHEN, 48.1N, 11.6E, 530m", y  = "jährlicher Niederschlag (mm)", x = "")+
  scale_x_continuous(breaks=seq(1955, 2020, 5))+
  scale_y_continuous(breaks=seq(600, 1350, 100))+
  geom_hline(yintercept = av, linetype = 2)+
  geom_segment(aes(x = 1954, y = 1200, xend = 1954, yend = av + 10),
               arrow = arrow(length = unit(0.3, "cm")))+
  annotate("text", x = 1954, y = 1210, label = "Durchschnitt = 958mm", size = 3, hjust = 0)

Niederschlag in München 2018

Das Skript: Die Details zu dem Skript finden sich im älteren Beitrag Niederschlag als Stufendiagramm. Die Daten werden mit read.csv() eingelesen. Fehlende Messwerte, in der Quelle als -999 angegeben, wandelt das Skript in NA um, damit R damit umgehen kann. lubridate erleichtert das Extrahieren der Jahre. Der Plot: ggplot2 hält für Stufendiagramme die Funktion geom_step() bereit. Die Parameter sind die X- und Y-Werte sowie die Farbe. Ergebnis: Mit 945 mm Niederschlag lag 2018 ziemlich nahe am Durchschnitt von 958 mm.

Teile diesen Beitrag!