2 min read

Methods

Species information

read_species() %>%
  select(scientific_name, acronym, common_name, biome, leaf_habit) %>%
  knitr::kable()
scientific_name acronym common_name biome leaf_habit
Acer negundo acene box elder temperate deciduous
Acer rubrum aceru red maple temperate deciduous
Acer saccharum acesa sugar maple temperate deciduous
Betula alleghaniensis betal yellow birch temperate deciduous
Betula papyrifera betpa paper birch boreal deciduous
Frangula alnus fraal glossy buckthorn exotic deciduous
Lonicera morrowii lonmo Morrow’s honeysuckle exotic deciduous
Lonicera tatarica lonta tatarian honeysuckle exotic deciduous
Populus tremouloides poptr trembling aspen boreal deciduous
Quercus macrocarpa quema bur oak temperate deciduous
Quercus rubra queru red oak temperate deciduous
Rhamnus cathartica rhaca common buckthorn exotic deciduous
Tilia americana tilam american basswood temperate deciduous
Abies balsamea abiba balsam fir boreal evergreen
Picea glauca picgl white spruce boreal evergreen
Picea mariana picma black spruce boreal evergreen
Pinus banksiana pinba jack pine boreal evergreen
Pinus resinosa pinre red pine temperate evergreen
Pinus strobus pinst white pine temperate evergreen
Thuja occidentalis thuoc white cedar temperate evergreen
Tsuga canadensis tsuca canadian hemlock temperate evergreen

Phenophase data

dat_phenophase <- tidy_phenophase()

setwd("phenologyb4warmed/")
usethis::use_data(dat_phenophase, overwrite = T)
setwd("..")
dat_phenophase_time_clean <- calc_phenophase_time(dat_phenophase, remove_outlier = T)
dat_phenophase_clean <- tidy_phenophase(remove_outlier = T, dat_phenophase_time = dat_phenophase_time_clean)

setwd("phenologyb4warmed/")
usethis::use_data(dat_phenophase_time_clean, overwrite = T)
usethis::use_data(dat_phenophase_clean, overwrite = T)
setwd("..")

Shoot length data

dat_shoot <- tidy_shoot()

usethis::proj_set(str_c(getwd(), "/phenologyb4warmed/"), force = T)
setwd("phenologyb4warmed/")
usethis::use_data(dat_shoot, overwrite = T)
setwd("..")
dat_shoot_extend <- dat_shoot %>% tidy_shoot_extend()
plot_shoot(dat_shoot_extend %>% filter(year == 2016, site == "hwrc") %>% filter(species == "queru"), quantile = F)

Climate data

dat_climate <- read_climate(option = "raw")
dat_climate_daily <- summ_climate_daily(dat_climate)

setwd("phenologyb4warmed/")
usethis::use_data(dat_climate_daily, overwrite = T)
setwd("..")
plot_climate_avail(dat_climate_daily, var = "temp_ag")
plot_climate_avail(dat_climate_daily, var = "mois")
dat_climate_spring <- summ_climate_season(dat_climate_daily, date_start = "Mar 15", date_end = "May 15", rainfall = 1)
dat_climate_spring_ambient <- calc_climate_ambient(dat_climate_spring, dat_phenophase)

plot_climate_change(dat_climate_spring_ambient)
dat_climate_fall <- summ_climate_season(dat_climate_daily, date_start = "Jul 15", date_end = "Sep 15", rainfall = 1)
dat_climate_fall_ambient <- calc_climate_ambient(dat_climate_fall, dat_phenophase)

plot_climate_change(dat_climate_fall_ambient)
dat_multiple <- tidy_multiple_data(dat_shoot, dat_phenophase_clean)

setwd("phenologyb4warmed/")
usethis::use_data(dat_multiple, overwrite = T)
setwd("..")
dat_multiple %>%
  filter(barcode == 6705, year == 2010) %>%
  filter(doy <= 230) %>%
  plot_multiple_data()