Figure 1

Load R packages

library(rainy)
library(patchwork)

Create the plots

p1 <- plot(bremen)
p2 <- plot(bremen, year = 1998, month = 1)
p3 <- plot(bremen, year = 1998, month = 1, day = 19)
p4 <- plot(bremen, year = 2023, month = 1)
p5 <- plot(bremen, year = 2023, month = 1, day = 30)
p1 / (p2 | p3) / (p4 | p5)

Data from Deutscher Wetterdienst (DWD), station 691.