-
Notifications
You must be signed in to change notification settings - Fork 4
/
README.Rmd
180 lines (132 loc) · 9.34 KB
/
README.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
---
output: github_document
---
<!-- README.md is generated from README.Rmd. Please edit that file -->
```{r, include = FALSE}
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = "man/figures/README-",
out.width = "100%"
)
```
# hockeystick <img src='man/figures/logo.png' align="right" height="139" />
<!-- badges: start -->
[](https://lifecycle.r-lib.org/articles/stages.html#stable)
[](https://CRAN.R-project.org/package=hockeystick)
[](https://cran.r-project.org/package=hockeystick)
<!-- badges: end -->
The goal of `hockeystick` is to make essential Global Heating datasets easily available in R to non-climate experts. `hockeystick` users can download the latest raw data from authoritative sources as well as view it via pre-defined ggplot2 charts. Datasets include atmospheric CO<sub>2</sub> and CH<sub>4</sub>, carbon emissions, instrumental, reconstructed, and paleo ice-core temperature records, sea levels, hurricanes, and Arctic/Antarctic sea-ice.
The choice of data was originally based on [Professor Stefan Rahmstorf's ](http://www.pik-potsdam.de/~stefan/) presentation on [The 5 Most Important Data Sets of Climate Science](http://www.pik-potsdam.de/~stefan/5datasets_rahmstorf.pdf). I came across this on a post on the [Open Mind blog](https://tamino.wordpress.com/2018/11/01/the-5-most-important-data-sets-of-climate-science/). I wrote my own [post](https://rethinking.rbind.io/2018/11/16/the-top-five-climate-charts-using-ggplot2/) on obtaining and visualizing this data (now out of date), which is the basis for this package. Additional datasets and visualizations have been added over The name of the package stems from the well known [hockeystick](https://en.wikipedia.org/wiki/Hockey_stick_graph) temperature chart.
`hockeystick` was highlighted on the RStudio [RViews](https://rviews.rstudio.com/2021/03/19/february-2021-top-40-new-cran-packages/) blog by Joseph Rickert as one of the "Top 40" new packages on CRAN in February 2021.
**New in version 0.8.0:** Daily temperature data from ClimateReanalyzer.com (see below).
New in version 0.7.0: Globally averaged methane (CH~4~) concentration from NOAA.
New in version 0.7.0: Cumulative emissions by country visualization.
New in version 0.6.0: Global CO~2~ emissions by region and country from GCP.
## Installation
To install the latest `hockeystick` release from CRAN type:
``` r
install.packages("hockeystick")
```
You may alternatively install the development version from [https://github.com/cortinah/hockeystick](https://github.com/cortinah/hockeystick) with:
``` r
remotes::install_github("cortinah/hockeystick")
```
## Downloading and viewing global heating data
New: Plot daily global temperatures since 1940 and current anomaly:
```{r dailytemp, fig.retina=2, message=FALSE, out.width='60%', fig.asp=1.03}
library(hockeystick)
globaldaily <- get_dailytemp()
plot_dailytemp(globaldaily)
```
Retrieve NOAA/ESRL Mauna Loa CO<sub>2</sub> Observatory concentration data and plot:
```{r carbon, fig.retina=2, message=FALSE, out.width='60%'}
ml_co2 <- get_carbon()
plot_carbon(ml_co2)
```
Retrieve GCP global CO~2~ emissions and plot:
```{r emissions, out.width='60%', fig.retina=2, message=FALSE}
emissions <- get_emissions()
plot_emissions(emissions)
plot_emissions_with_land(emissions)
```
Visualize cumulative emissions by country:
```{r emissionsmap, out.width='60%', fig.retina=2}
emissions_map()
```
Retrieve NASA/GISS global surface temperature anomaly data and plot:
```{r temp, out.width='60%', fig.retina=2}
anomaly <- get_temp()
plot_temp(anomaly)
```
Visualize warming using Ed Hawkins styled "warming stripes":
```{r stripes, out.width='60%', fig.retina=2}
warming_stripes()
```
```{r stripes2, out.width='53%', fig.retina=2}
warming_stripes(stripe_only = TRUE, col_strip = viridisLite::viridis(11))
```
Retrieve tide gauge and satellite sea level data and plot:
```{r sl, out.width='60%', fig.retina=2}
gmsl <- get_sealevel()
plot_sealevel(gmsl)
```
Retrieve July annual Arctic Sea Ice Index and plot:
```{r si, fig.retina=2, message=FALSE, out.width='60%', warning=FALSE}
seaice <- get_seaice()
plot_seaice(seaice)
```
`get_seaice()` arguments can be modified to download Antarctic sea ice, and allow any month.
You can also visualize sea ice by month and year:
```{r icecurves, fig.retina=2, message=FALSE, out.width='60%'}
arcticice <- get_icecurves()
plot_icecurves(arcticice)
```
Retrieve Common Era temperature reconstruction and plot it with instrumental record:
```{r 2ktemp, out.width='60%', fig.retina=2}
anomaly2k <- get_temp2k()
plot_temp2k(anomaly2k)
```
Retrieve NOAA HURDAT2 hurricane data and plot:
```{r hurricanes, out.width='60%', fig.retina=2, message=FALSE}
hurricanes <- get_hurricanes()
plot_hurricanes(hurricanes)
plot_hurricane_nrg(hurricanes)
```
Retrieve NOAA/ESRL CH<sub>4</sub> Globally averaged mean data and plot:
```{r methane, fig.retina=2, message=FALSE, out.width='60%'}
ch4 <- get_methane()
plot_methane(ch4)
```
Retrieve Vostok paleo ice core data and plot:
```{r paleo, out.width='60%', fig.retina=2, message=FALSE}
vostok <- get_paleo()
plot_paleo(vostok)
```
### Managing the cache
By default, no climate data is cached, and all data is downloaded every time any of the `get_` functions is called. To cache data for future use, use the `write_cache = TRUE` option, available in all of the `get_` functions. To download and cache all data use `hockeystick_update_all()`. To view the files, date, and size of cached data use `hockeystick_cache_details()`. To re-download data from the source use the `use_cache = FALSE` argument in any of the `get_` functions, for example: `get_carbon(use_cache = FALSE, write_cache = TRUE)`. To delete all cached data use `hockeystick_cache_delete_all()`.
Users may also cache data by default by adding `options(hs_write_cache = TRUE)`to their script or `.Rprofile` file.
### All together now: climate data grid
```{r grid, out.width='80%', fig.retina=2, warning=FALSE}
climate_grid()
```
## Acknowledgments
- Carbon Dioxide concentrations: Dr. Pieter Tans, NOAA/GML (https://gml.noaa.gov/ccgg/trends/) and Dr. Ralph Keeling, [Scripps Institution of Oceanography](https://scrippsco2.ucsd.edu/).
- Global temperature anomaly: *GISS Surface Temperature Analysis (GISTEMP), version 4.* GISTEMP Team, 2020: NASA Goddard Institute for Space Studies. https://data.giss.nasa.gov/gistemp/
- Warming Stripes design from Ed Hawkins' Climate Lab. https://www.climate-lab-book.ac.uk/2018/warming-stripes/. In addition grateful to Dr. Dominic Royé for posting his approach to plotting them using ggplot2, which `warming_stripes()` is based on. https://dominicroye.github.io/en/2018/how-to-create-warming-stripes-in-r/
- Sea level data: [NOAA Laboratory for Satellite Altimetry (sat)](https://www.star.nesdis.noaa.gov/socd/lsa/SeaLevelRise/) and [Commonwealth Scientific and Industrial Research Organisation (tide gauges)](https://research.csiro.au/slrwavescoast/sea-level/measurements-and-data/sea-level-data/)
- Sea Ice Index: [National Snow & Ice Data Center](https://nsidc.org/home).
Data Archive: https://nsidc.org/data/explore-data
- Vostok carbon dioxide and temperature data: https://data.ess-dive.lbl.gov/datasets/doi:10.3334/CDIAC/ATG.009
- Common Era reconstructed temperature data: [PAGES2k Consortium and NOAA](https://www.ncei.noaa.gov/access/paleo-search/study/26872)).
- Hurricanes: National Oceanic and Atmospheric Administration HURDAT Atlantic Hurricane Database [Re-analysis Project](https://www.aoml.noaa.gov/hrd/hurdat/Data_Storm.html), particularly Dr. Chris Landsea.
- Carbon Dioxide emissions: [Global Carbon Project](https://www.globalcarbonproject.org/carbonbudget/) and [Our World In Data](https://github.com/owid/co2-data)
- Methane: Lan, X., K.W. Thoning, and E.J. Dlugokencky, NOAA Global Monitoring Laboratory (https://gml.noaa.gov/ccgg/trends_ch4/).
- Thank you to Dirk Eddelbuettel for providing the .isConnected function from his [tint](https://CRAN.R-project.org/package=tint) package to test for internet connectivity.
- Thank you to Chris Campbell and Steven Bernard of the FT Climate Graphics Explained newsletter for sharing their code for climatereanalyzer.org daily temperature chart (Temperature Records Broken...Again, 9 July 2023). This code is the foundation of the get_dailytemp() and plot_dailytemp() functions.
## Notes and resources
- All data is compiled to the best of my ability from reliable and peer-reviewed sources. Please open an issue if you are aware of enhanced or additional data that may be added to the package. Building `hockeystick` is driven by my interest in tracking climate data and making it easily available to the community.
- Here are some online resources I have found very helpful to learn mo re about climate science:
- MIT edX Global Warming Science. https://www.edx.org/learn/global-warming/massachusetts-institute-of-technology-global-warming-science
- SDG Academy: Climate Change: The Science and Global Impact. https://sdgacademy.org/course/climate-change-the-science-and-global-impact/