climproxyrecords

The climproxyrecords package supplies published proxy climate record data and compilations. So far it includes data from:

For the two compilations the package contains three dataframes: proxies, dating, and metadata (e.g. marcott.proxies). These contain the proxy data, associated carbon dating information, and metadata respectively.

For the individual proxies there are similar dataframes.

alley.accumulation

alley.temperature

sicre.2011.MD99_2275.age.model

sicre.2011.MD99_2275.metadata

sicre.2011.MD99_2275.temperature

naafs.U1313.temperature

naafs.U1313.metadata

herbert.ODP846.temperature

herbert.ODP846.metadata

Installation

climproxyrecords can be installed from github like this:

if (!require("devtools")) {
  install.packages("devtools")
}

devtools::install_github("EarthSystemDiagnostics/climproxyrecords")

Citing the data

For the compilations, please cite both the compilaton papers (Marcott et al 2013, Shakun et al 2012), and the original sources from which the data were extracted. A full list of sources can be found in the metadata dataframes, e.g. marcott.metadata.

For the other records please cite the original papers. See the help files e.g. help(sicre.2011.MD99_2275.temperature) for the correct citations.

Using the climproxyrecords data package

After loading the package with library(climproxyrecords), all included data has been "lazy loaded". Data are not in memory yet, but are available when the name of an object is passed to the R interpreter.

library(climproxyrecords)
library(dplyr)

marcott.proxies

Subsetting

For each proxy compilation, all records are stored in a single dataframe. To access subsets of proxies use either base R's subsetting

subset(marcott.proxies, Proxy.type == "Uk'37")

or dplyr::filter

marcott.proxies %>% 
  filter(Proxy.type == "Uk'37")

Split to a list

or make a list of dataframes, each a separate proxy record

lst <- plyr::dlply(marcott.proxies, "ID")

lapply(lst[1:2], head)

Warnings

Core.location does not identify a unique proxy record (for the Marcott data at least) because some cores have multiple proxies measured on them.

name.type <- marcott.proxies %>% 
  select(Core.location, Proxy.type) %>% 
  distinct()

name.type %>% 
  group_by(Core.location) %>% 
  mutate(n = n_distinct(Proxy.type)) %>% 
  filter(n > 1)

At present, things like the names of proxy types have been taken "as is" from the supplementary data to the publications. There may be multiple spellings of the same type.

Get metadata

To get additional metadata for the proxies, look at corresponding metadata dataframe and reference/join by ID

dplyr::left_join(marcott.proxies, marcott.metadata)

marcott.metadata %>% 
  filter(Core.location == "GeoB5844-2")


EarthSystemDiagnostics/climproxyrecords documentation built on Feb. 12, 2022, 8:48 p.m.