INSTRUCTIONS (Delete when you compile the submission version of your project.)

Modify your report YAML with the following.

Remove the comments in the dataset_chunk to load the assigned datasets.

If you are including a numerical result in your text, you should use inline code.

All plots should have informative titles and axis labels. You may also change colors.

THIS DOCUMENT WAS EDITED FOR SPELLING / GRAMMAR / CLARITY TUESDAY 02/26

MAUNA LOA CO2 MODELING AND VISUALIZATION

The Question

In 1966, the World Meteorological Organization (WMO) put forth the term “climatic change” to refer to climatic variability on time-scales longer than ten years, regardless of the cause for such change. During the next decade, scientists began to suspect that human activities had the potential to drastically alter the global climate in ways that would have negative impacts on our environment. The term evolved into “climate change” and is now used to describe both the process of change and the perceived problem. Sometimes the term “global warming” is used, though in many ways this fails to adequately describe the variability in impact, since climate change can cause both hot and cold extremes in weather. Anthropogenic climate change is change that is caused by human activity, as opposed to the Earth’s natural processes. However, in the context of environmental policy, the term “climate change” is often used to mean anthropogenic climate change.

Mauna Loa Observatory is a world-renowned atmospheric research facility. It has been continuously monitoring and collecting data since the 1950’s and its remote location makes it very well-suited for monitoring atmospheric components that can contribute to climate change, including the heat-trapping greenhouse gas carbon dioxide (CO2). Carbon overload from burning fossil fuels and deforestation is cited as the primary cause of anthropogenic climate change by proponents of such theories, while opponents assert that natural process (such as photosynthesis) contribute more to atmospheric CO2 than humans and observed changes are simply Earth’s cycle.

Monthly Mean CO2: The Last Five Years

Create your own version of the plot found here. Do not replicate it, but rather design your own. Use one of the themes found in the ggplot2 or ggthemes packages. You are encouraged to make style adjustments to help you informatively display the data.

co2_monthly2015 <- co2_monthly %>% filter(year >= 2015)
#co2_monthly2015



ggplot(co2_monthly2015) +
  geom_line(aes(x = date, y = int_mean_co2), col = "light blue", linetype = 10) +
  geom_point(aes(x = date, y = int_mean_co2), col = "purple", shape = 18) +
  geom_line(aes(x = date, y = trend_mean_co2), col = "black") +
  geom_point(aes(x = date, y = trend_mean_co2), col = "black", shape = 15) +
  scale_x_continuous(breaks = seq(2015, 2020, 1),
                     limits = c(2015, 2020),
                     sec.axis = dup_axis(labels = NULL, name = "")) +
  scale_y_continuous(breaks = seq(395, 415, 5),
                     limits = c(395, 415),
                     sec.axis = dup_axis(labels = NULL, name = "")) +
  ggtitle(expression("RECENT monthly MEAN CO"[2]*" AT MAUNA LOA")) +
  ylab("PARTS PER MILLION") +
  xlab("YEAR") +
  theme_linedraw()

Monthly Mean CO2: A Major Milestone

An atmospheric CO2 level of 400 ppm is considered by many to be a symbolic threshold with regard to climate change. “In the centuries to come, history books will likely look back on September 2016 as a major milestone for the world’s climate. At a time when atmospheric carbon dioxide is usually at its minimum, the monthly value failed to drop below 400 parts per million.” (source)

Adapt your plot above to include a red dashed line at 400 ppm and a large red dot on September 2016, with appropriate annotations to indicate what these additions represent.

ggplot(co2_monthly2015) +
  geom_line(aes(x = date, y = int_mean_co2), col = "light blue", linetype = 10) +
  geom_point(aes(x = date, y = int_mean_co2), col = "purple", shape = 18) +
  geom_line(aes(x = date, y = trend_mean_co2), col = "black") +
  geom_point(aes(x = date, y = trend_mean_co2), col = "black", shape = 15) +
  scale_x_continuous(breaks = seq(2015, 2020, 1),
                     limits = c(2015, 2020),
                     sec.axis = dup_axis(labels = NULL, name = "")) +
  scale_y_continuous(breaks = seq(395, 415, 5),
                     limits = c(395, 415),
                     sec.axis = dup_axis(labels = NULL, name = "")) +
  ggtitle(expression("RECENT monthly MEAN CO"[2]*" AT MAUNA LOA")) +
  ylab("PARTS PER MILLION") +
  xlab("YEAR") +
  theme_linedraw() +
  geom_hline(yintercept = 400, color = "red", linetype = 8) +
  annotate("point", x = 2016.708, y = 401.05, colour = "red", size = 3)

Annual Mean CO2 Since 1959

Replicate as closely as possible the annual mean plot found here. Hint: It uses a ggplot theme for some of the formatting.

### Figure Out How to Do geom_label and geom_text ###

ggplot(co2_annual, aes(x = year, y = mean_co2 )) + 
  labs(title = expression("Atmospheric CO"[2]*" at Mauna Loa Observatory"), 
      subtitle = "with loess smoothed trend curve and estimated historical reference values",
      x = "measurement year", y = expression("CO"[2]*" (ppm)"), 
      caption = "data source: http://www.esrl.noaa.gov/gmd/ccgg/trends/data.html") +   
  geom_bar(stat = "identity", col = "light blue", fill = "light blue", width = .5) +  
  geom_hline(yintercept = 400, color = "red") + 
  annotate("label", label = "crisis threshold", x = 1988, y = 400)+
  geom_smooth() + 
  geom_hline(yintercept = 200) + 
  annotate("label", label = "ice age mean", x = 1988, y = 200)+
  geom_hline(yintercept = 280) + 
  annotate("label", label = "pre-industrial mean", x = 1988, y = 280)+
  scale_x_continuous(breaks = seq(1960, 2020, 5)) + 
  scale_y_continuous(breaks = seq(0, 400, 50)) + 
  theme_classic()