Rainfall in 2021 summer

今年（2021）夏天的降雨比平时多很多，正好手头上有收集的降雨数据。于是稍作分析，展示出来。

记得往年的汛期都是7月中下旬（有历史数据为证），然而今年6月底到7月初的降雨已经接近150毫米啦！

下面情况今年北京的降雨情况分析。

ggplot(pre, aes(x = Datetime, y = PRE_sum)) + 
  geom_line() +
  scale_x_datetime(date_labels = "%b %d", date_breaks = "3 days",) +
  theme_cowplot() + mytheme +
  theme(axis.text.x = element_text(angle=45, hjust = 1)) 

如图所见，6月初累计降雨已达近150毫米，而且最近两周降雨强度较大（如下图）。降水最多的一天发生在7月2号，单次降水量达到40毫米

ggplot(PRE.daily, aes(x = Date, y = PRE_daily)) + geom_col() +
  scale_x_date(date_breaks = "2 days", labels = date_format("%b %d")) +
  scale_y_continuous(expand = expansion()) +
  # Two weeks
  geom_point(aes(y = PRE_sum14d),size = 3) +
  geom_line(aes(y = PRE_sum14d)) +
  geom_label_repel(aes(y = PRE_sum14d,label = label14d),
                   nudge_x = 1,
                   #nudge_y = 2,
                   box.padding = unit(1.2, "lines"),
                   na.rm = TRUE) +
  theme_cowplot() +theme(axis.text.x = element_text(angle=45, hjust = 1)) +
  mytheme