Is the weather getting wetter?

February 20, 2019 - 4 minutes
Exploring historical data from 1905 to 2015 from the World Bank
weather lm sf tidyverse

It's gon rain

Intro

Is the weather really getting wetter? It sure does feel like it from where I live. This post is about testing my feelings with linear models.

Data

The first step was to find historical data, which I found on the World Bank’s Climate Change Knowledge Portal. This dataset has country level precipitation records from 1905 to 2015. You are allowed to select up to 30 countries data, so I focused focus on the Western Hemisphere, where I live.

library(tidyverse)

rain <- read_csv("~/Documents/pr_1901_2015.csv") %>% # download available as Excel or CSV
  rename_all(tolower) %>% 
  select(iso_alpha = country, year, month, rainfall = pr)

The data comes labeled with ISO country codes, so I used data from the gapminder pacakges to add in the full country name.

library(magrittr) # viva la %<>%

data("country_codes", package = "gapminder")

rain %<>%
  inner_join(country_codes) %>% 
  mutate(country = gsub("United States", "USA", country)) # wasn't auto-matching

Trend over time

This dataset includes month level records for each country since 1905, but to look at the overall trend for the last century, I aggregated the data into yearly totals by country. The units for rainfall are millimeters (mm).

yearly <- rain %>%
  group_by(year, country) %>%
  summarise(rainfall = sum(rainfall))

Then plot them all out in a big facet panel.

ggplot(yearly, aes(year, rainfall, color = country, group = country)) +
  geom_path(alpha = .5) +
  stat_smooth(method = "lm", aes(fill = country)) +
  scale_x_continuous(breaks = c(1900, 1950, 2000)) +
  facet_wrap(~country, scales = "free", nrow = 6) +
  theme(legend.position = "none") +
  labs("Yearly rainfall by country",
       y = "precipitation (mm)",
       x = NULL)

Annual volumes do vary a lot, but we can still see several clear upward trends, like Argentinia and Canada, where rainfall has increased over the last century.

A map for geographic effect

Finally to see the model data geographically, I made a map with geom_sf().

library(sf)

world <- maps::map('world', plot = FALSE, fill = TRUE) %>%
  st_as_sf() %>%
  rename(country = ID)

inner_join(world, time_table) %>%
  ggplot() +
  geom_sf(aes(fill = estimate, color = p.value < .05)) +
  scale_fill_gradient2(name = "yearly change (mm)") +
  scale_color_manual(values = c("grey", "black")) +
  coord_sf(xlim = c(-180, 0)) +
  theme_void() +
  labs(title = "Yearly change in precipication from 1901-2015")

So most of the countries with ambigous trends are close to the equator, Chile being the exception. It also looks like the countries closest to the south Atlantic see the largest yearly increases.

Conclusion

I’m finishing this post on a day where it is snowing/raining yet again. So next time I have that small-talk conversation about how wet its been lately, I’ll be ready with statistics from the past century, to show that it really is getting wetter.

Geocoded crime reports for Charlottesville Virginia

November 27, 2018 - 5 minutes
Civic Data packages sf tidyverse

Quantifying time in meetings

September 2, 2018 - 5 minutes
Using Google Apps Script and R to analyze Google Calendars.
Business Intelligence googlesheets tidyverse

Maps with the new ggplot2 v3.0.0

August 4, 2018 - 2 minutes
Civic Data ggplot2 tidyverse