WF sub-chronic exposure

Import data with (WF) PM2.5 for California ZIP codes (2006-2021) Subset for Southern California

library(dplyr);

load("~/UCSD-Scripps/R Files 2022/WF sub-chronic/wfpm25_06t021.RData")

length(unique(wfpm25$zip))

wfpm25$wf_pm25_gt0 <- ifelse(wfpm25$wf_pm25_imp > 0, wfpm25$wf_pm25_imp, NA)
wfpm25$smoke_gt0 <- ifelse(wfpm25$smoke_intersect > 0, wfpm25$smoke_intersect, NA)

Calculate moving averages (WF PM2.5) for 3 and 5 years, and 1-year rolling sum of smoke days

library(zoo);

wf_subchronic <- wfpm25 %>%
    group_by(zip) %>% 
    arrange(zip, date) %>%
    mutate(smokeday_sum_1yr = zoo::rollapply(smoke_gt0, 365, sum, fill = NA, na.rm = TRUE, 
                                             align = "right", partial = TRUE),
           wfpm25gt0_mean_3yr = zoo::rollapply(wf_pm25_gt0, 1095, mean, fill = NA, na.rm = TRUE, 
                                               align = "right", partial = TRUE),
           wfpm25gt0_mean_5yr = zoo::rollapply(wf_pm25_gt0, 1825, mean, fill = NA, na.rm = TRUE, 
                                               align = "right", partial = TRUE)) #,

wf_subchronic$wfpm25gt0_mean_5yr[is.na(wf_subchronic$wfpm25gt0_mean_5yr)] = 0
wf_subchronic$wfpm25gt0_mean_5yr[is.na(wf_subchronic$wfpm25gt0_mean_5yr)] = 0

Pick one date as an example: March 19 2020

library(sf);
library(ggplot2);

setwd("C:/Users/Rosana/Documents/UCSD-Scripps/R Files 2022/WF sub-chronic/CAzip")
zipCA <- st_read(dsn = ".", layer = "USPS_zipCA_byESRI")

ca.data <- ggplot2::fortify(zipCA, region='ZIPCODE') # shp to df

region <- st_union(ca.data) # shp for state boundary

Merge ZIP code polygons and Wildfire Data Mean Wildfire PM2.5 for the last 5 years on March 19, 2020

one.date <- wfpm25_5yr[wfpm25_5yr$date == "2020-03-19", ]

ca.map <- merge(ca.data, one.date, by.x = "ZIPCODE", by.y = "zip", all.x = TRUE)

my.col <- c("#FFFE9E", "#F69422", "#C53270", "#611163", "#040404")

zip <- ggplot(ca.map) +
  geom_sf(size = 0.1, color = NA, aes(fill= wfpm25gt0_mean_5yr)) +  
  ggtitle(expression(bold(paste("5-year Mean Wildfire ", PM[2.5], " (",mu, "g ", m^bold("-3"),") for March 19, 2020", )))) + 
  theme(plot.title = element_text(hjust = 0.5)) +  # center title
  scale_fill_gradientn(colours = my.col, na.value = "light grey", limits = c(0,45)) +
  
  theme(legend.key.size = unit(0.5, "cm"), legend.key.height= unit(1,"cm") ) +
  theme(plot.title = element_text(size=14, face="bold")) + 
  theme(legend.text=element_text(size=12, face="bold")) +
  theme(axis.title.x=element_blank(), axis.title.y=element_blank()) +
  
  theme(legend.position="right") + theme(legend.title = element_blank()) +
  labs(subtitle = (expression(bold(paste("(averaged for days when WF " , PM[2.5], "> 0)")))), 
       caption = "Unpopulated areas shown in gray") +
  theme(plot.subtitle=element_text(size=12, hjust=0.5, face="italic", color="black")) +
  theme(
    # Hide panel borders and remove grid lines
    panel.border = element_rect(colour = "black", fill= "NA"),
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(), 
    panel.background = element_rect(fill = "white", colour = NA)) +
  theme(axis.title.x=element_blank(),
        axis.text.x=element_blank(),
        axis.ticks.x=element_blank()) +
  theme(axis.text.y=element_blank(), axis.ticks.y=element_blank()) +
  theme(plot.caption = element_text(hjust = 0.5, size = 12, face= "italic"))

zip

Sum of Smoke Days for the last 1 year on October 1, 2008

another.date <- wfpm25_5yr[wfpm25_5yr$date == "2008-10-01", ]

ca.map <- merge(ca.data, another.date, by.x = "ZIPCODE", by.y = "zip", all.x = TRUE)

my.col <- c("#FFFE9E", "#F69422", "#C53270", "#611163", "#040404")

smoke <- ggplot(ca.map) +
  geom_sf(size = 0.1, color = NA, aes(fill= smokeday_sum_1yr)) +  
  ggtitle("1-year Sum of Smoke Days for October 10, 2008") + 
  theme(plot.title = element_text(hjust = 0.5)) +  # center title
  scale_fill_gradientn(colours = my.col, na.value = "light grey", limits = c(0,105)) +
  
  theme(legend.key.size = unit(0.5, "cm"), legend.key.height= unit(1,"cm") ) +
  theme(plot.title = element_text(size=14, face="bold")) + 
  theme(legend.text=element_text(size=12, face="bold")) +
  theme(axis.title.x=element_blank(), axis.title.y=element_blank()) +
  
  theme(legend.position="right") + theme(legend.title = element_blank()) +
  labs(caption = "Unpopulated areas shown in gray") +
  theme(plot.subtitle=element_text(size=12, hjust=0.5, face="italic", color="black")) +
  theme(
    # Hide panel borders and remove grid lines
    panel.border = element_rect(colour = "black", fill= "NA"),
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(), 
    panel.background = element_rect(fill = "white", colour = NA)) +
  theme(axis.title.x=element_blank(),
        axis.text.x=element_blank(),
        axis.ticks.x=element_blank()) +
  theme(axis.text.y=element_blank(), axis.ticks.y=element_blank()) +
  theme(plot.caption = element_text(hjust = 0.5, size = 12, face= "italic"))

smoke

Mean Wildfire PM2.5 for the last 3 years on September 1, 2012

yet.another.date <- wf_subchronic[wf_subchronic$date == "2012-09-01", ]

ca.map <- merge(ca.data, yet.another.date, by.x = "ZIPCODE", by.y = "zip", all.x = TRUE)

my.col <- c("#FFFE9E", "#F69422", "#C53270", "#611163", "#040404")

wf <- ggplot(ca.map) +
  geom_sf(size = 0.1, color = NA, aes(fill=  wfpm25gt0_mean_3yr)) +  
  ggtitle(expression(bold(paste("3-year Mean Wildfire ", PM[2.5], " (",mu, "g ", m^bold("-3"),") for September 1, 2012", )))) + 
  theme(plot.title = element_text(hjust = 0.5)) +  # center title
  scale_fill_gradientn(colours = my.col, na.value = "light grey", limits = c(0,10)) +
  
  theme(legend.key.size = unit(0.5, "cm"), legend.key.height= unit(1,"cm") ) +
  theme(plot.title = element_text(size=14, face="bold")) + 
  theme(legend.text=element_text(size=12, face="bold")) +
  theme(axis.title.x=element_blank(), axis.title.y=element_blank()) +
  
  theme(legend.position="right") + theme(legend.title = element_blank()) +
  labs(subtitle = (expression(bold(paste("(averaged for days when WF " , PM[2.5], "> 0)")))), 
       caption = "Unpopulated areas shown in gray") +
  theme(plot.subtitle=element_text(size=12, hjust=0.5, face="italic", color="black")) +
  theme(
    # Hide panel borders and remove grid lines
    panel.border = element_rect(colour = "black", fill= "NA"),
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(), 
    panel.background = element_rect(fill = "white", colour = NA)) +
  theme(axis.title.x=element_blank(),
        axis.text.x=element_blank(),
        axis.ticks.x=element_blank()) +
  theme(axis.text.y=element_blank(), axis.ticks.y=element_blank()) +
  theme(plot.caption = element_text(hjust = 0.5, size = 12, face= "italic"))

wf