Yamethin Case Study - Climate Change Scenarios

The maps were produced prompt action to requirement by Universities and public sector. This map series include climate change projection map and (%) change from the baseline on two RCPs (i.e RCP 4.5 and 8.5). The maps were produced using R deploying terra, tidyverse, tmap, raster packages.

Climate Projection - RCP 4.5

The whole country

map1 <- tm_shape(raster(Prec_45_2040_PM))+
    tm_raster(palette = "-viridis", title = "Preci - 2040 RCP45") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map2 <- tm_shape(raster(Prec_45_2060_PM))+
    tm_raster(palette = "-viridis", style = "pretty", title = "Preci - 2060 RCP45") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map3 <- tm_shape(raster(Prec_45_2080_PM))+
    tm_raster(palette = "-viridis", style = "pretty", title = "Preci - 2080 RCP45") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map4 <- tm_shape(raster(Prec_45_2100_PM))+
    tm_raster(palette = "-viridis", style = "pretty", title = "Preci - 2100 RCP45") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
tmap_arrange(map1, map2, map3, map4, nrow = 2)

Yamethin Township

plot(crop(Prec_45_2020_PM, ymt), main = "Precipitation Projection Scenario for 2020, RCP 4.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_45_2040_PM, ymt), main = "Precipitation Projection Scenario for 2040, RCP 4.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_45_2060_PM, ymt), main = "Precipitation Projection Scenario for 2060, RCP 4.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_45_2080_PM, ymt), main = "Precipitation Projection Scenario for 2080, RCP 4.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_45_2100_PM, ymt), main = "Precipitation Projection Scenario for 2100, RCP 4.5\n Yamethin Township")
lines(ymt, col = "blue")

(%) Change from Baseline

Whole country

map1 <- tm_shape(raster(Prec_45_2040_CM))+
    tm_raster(palette = "BuGn", title = "Preci - 2040 RCP45") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map2 <- tm_shape(raster(Prec_45_2060_CM))+
    tm_raster(palette = "BuGn", style = "pretty", title = "Preci - 2060 RCP45") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map3 <- tm_shape(raster(Prec_45_2080_CM))+
    tm_raster(palette = "BuGn", style = "pretty", title = "Preci - 2080 RCP45") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map4 <- tm_shape(raster(Prec_45_2100_CM))+
    tm_raster(palette = "BuGn", style = "pretty", title = "Preci - 2100 RCP45") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
tmap_arrange(map1, map2, map3, map4, nrow = 2)

Yamethin Township

plot(crop(Prec_45_2020_CM, ymt), main = "Precipitation Scenario for 2020-  % Change from baseline, RCP 4.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_45_2040_CM, ymt), main = "Precipitation Scenario for 2040 - % Change from baseline, RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_45_2060_CM, ymt), main = "Precipitation Scenario for 2060 - % Change from baseline, RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_45_2080_CM, ymt), main = "Precipitation Scenario for 2080 - % Change from baseline, RCP 4.5 \nYamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(Prec_45_2100_CM, ymt), main = "Precipitation Scenario for 2100 -  % Change from baseline, RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

RCP 8.5

The Whole Country

map1 <- tm_shape(raster(Prec_85_2040_PM))+
    tm_raster(palette = "-viridis", title = "Preci - 2040 RCP85") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map2 <- tm_shape(raster(Prec_85_2060_PM))+
    tm_raster(palette = "-viridis", style = "pretty", title = "Preci - 2060 RCP85") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map3 <- tm_shape(raster(Prec_85_2080_PM))+
    tm_raster(palette = "-viridis", style = "pretty", title = "Preci - 2080 RCP85") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map4 <- tm_shape(raster(Prec_85_2100_PM))+
    tm_raster(palette = "-viridis", style = "pretty", title = "Preci - 2100 RCP85") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
tmap_arrange(map1, map2, map3, map4, nrow = 2)

Yamethin

plot(crop(Prec_85_2020_PM, ymt), main = "Precipitation Projection Scenario for 2020, RCP 8.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(Prec_85_2040_PM, ymt), main = "Precipitation Projection Scenario for 2040, RCP 8.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(Prec_85_2060_PM, ymt), main = "Precipitation Projection Scenario for 2060, RCP 8.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(Prec_85_2080_PM, ymt), main = "Precipitation Projection Scenario for 2080, RCP 8.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(Prec_85_2100_PM, ymt), main = "Precipitation Projection Scenario for 2100, RCP 8.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

(%) Change from Baseline - RCP 8.5

The Whole Country

map1 <- tm_shape(raster(Prec_85_2040_CM))+
    tm_raster(palette = "BuGn", style = "pretty", title = "Preci - 2040 RCP85") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map2 <- tm_shape(raster(Prec_85_2060_CM))+
    tm_raster(palette = "BuGn", style = "pretty", title = "Preci - 2060 RCP85") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map3 <- tm_shape(raster(Prec_85_2080_CM))+
    tm_raster(palette = "BuGn", style = "pretty", title = "Preci - 2080 RCP85") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map4 <- tm_shape(raster(Prec_85_2100_CM))+
    tm_raster(palette = "BuGn", style = "pretty", title = "Preci - 2100 RCP85") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
tmap_arrange(map1, map2, map3, map4, nrow = 2)

Yamethin Township

plot(crop(Prec_85_2020_CM, ymt), main = "Precipitation Scenario for 2020-  % Change from baseline, RCP 8.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_85_2040_CM, ymt), main = "Precipitation Scenario for 2040-  % Change from baseline, RCP 8.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_85_2060_CM, ymt), main = "Precipitation Scenario for 2060-  % Change from baseline, RCP 8.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_85_2080_CM, ymt), main = "Precipitation Scenario for 2080-  % Change from baseline, RCP 8.5\n Yamethin Township")
lines(ymt, col = "blue")

plot(crop(Prec_85_2100_CM, ymt), main = "Precipitation Scenario for 2100-  % Change from baseline, RCP 4.5\n Yamethin Township")
lines(ymt, col = "blue")

Mean Temperature Projection RCP 4.5

Whole Country

map1 <- tm_shape(raster(TMean_45_2040_PM))+
    tm_raster(style = "quantile", palette = "Oranges", title = "Tmean - 2040 RCP4.5") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map2 <- tm_shape(raster(TMean_45_2060_PM))+
    tm_raster(style = "quantile", palette = "Oranges", title = "Tmean - 2060 RCP4.5") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map3 <- tm_shape(raster(TMean_45_2080_PM))+
    tm_raster(style = "quantile", palette = "Oranges", title = "Tmean - 2080 RCP4.5") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map4 <- tm_shape(raster(TMean_45_2100_PM))+
    tm_raster(style = "quantile", palette = "Oranges",  title = "Tmean - 2100 RCP4.5") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
tmap_arrange(map1, map2, map3, map4, nrow = 2)

Yamethin Township

par(mfrow = c(1,1))
brks <- seq(0, 1, by=0.15) 
nb <- length(brks)-1 
plot(crop(TMean_45_2020_PM, ymt), col = rev(heat.colors(nb)),lab.breaks = brks, alpha = 1, main = "Mean Temperature Projection Scenario for 2020, RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(TMean_45_2040_PM, ymt), col = rev(heat.colors(nb)),lab.breaks = brks, main = "Mean Temperature Projection Scenario for 2040, RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(TMean_45_2060_PM, ymt), col = rev(heat.colors(nb)),lab.breaks = brks, main = "Mean Temperature Projection Scenario for 2060, RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(TMean_45_2080_PM, ymt), col = rev(heat.colors(nb)),lab.breaks = brks, main = "Mean Temperature Projection Scenario for 2080, RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(TMean_45_2100_PM, ymt), col = rev(heat.colors(nb)),lab.breaks = brks, main = "Mean Temperature Projection Scenario for 2100, RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

#plot(TMean_Ensem_Change_45, col = rev(heat.colors(20)))

ªC change from baseline RCP 4.5

Whole Country

map1 <- tm_shape(raster(TMean_45_2040_CM))+
    tm_raster(style = "quantile", palette = "Oranges", title = "Tmean - 2040 RCP4.5") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map2 <- tm_shape(raster(TMean_45_2060_CM))+
    tm_raster(style = "quantile", palette = "Oranges", title = "Tmean - 2060 RCP4.5") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map3 <- tm_shape(raster(TMean_45_2080_CM))+
    tm_raster(style = "quantile", palette = "Oranges", title = "Tmean - 2080 RCP4.5") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
map4 <- tm_shape(raster(TMean_45_2100_CM))+
    tm_raster(style = "quantile", palette = "Oranges",  title = "Tmean - 2100 RCP4.5") +
    tm_layout(scale=.8, legend.position = c("left","bottom"))
tmap_arrange(map1, map2, map3, map4, nrow = 2)

Yamethin Township

par(mfrow = c(1,1))
brks <- seq(0, 1, by=0.15) 
nb <- length(brks)-1 
plot(crop(TMean_45_2020_CM, ymt), col = rev(heat.colors(20)), alpha = 0.1, main = "Mean Temperature Projection Scenario for 2020, ºC change from baseline RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(TMean_45_2040_CM, ymt), col = rev(heat.colors(20)), main = "Mean Temperature Projection Scenario for 2040, ºC change from baseline RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(TMean_45_2060_CM, ymt), col = rev(heat.colors(20)), main = "Mean Temperature Projection Scenario for 2060, ºC change from baseline RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(TMean_45_2080_CM, ymt), col = rev(heat.colors(20)), main = "Mean Temperature Projection Scenario for 2080, ºC change from baseline RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")

plot(crop(TMean_45_2100_CM, ymt), col = rev(heat.colors(20)), main = "Mean Temperature Scenario for 2100, ºC change from baseline RCP 4.5 \n Yamethin Township")
lines(ymt, col = "blue", type = "l")