The cleaned half-hourly psychrometer values have been summarized to daily predawn (mean of 2 hours prior to sunrise, inclusive) and midday (mean of 2 hours following solar noon, inclusive). Similarly, the half-hourly environmental variables have been summarized to daily sums or statistics (mean, min, max), depending on the variable.
load("../data_cleaned/met_daily.Rdata")
load("../data_cleaned/psy_daily.Rdata")
Next, we will exclusively consider the predawn and midday means matched with the daily environmental variables. For plotting purposes, we will also include the maintenance times. After each maintenance, the logger was switched to a new, previously un-instrumented branch. Therefore, the identity of Logger 1 and Logger 2 are arbitrary and considered replicates for the same time period.
met_clip <- met_daily %>%
filter(date >= min(psy_daily$date),
date <= max(psy_daily$date))
psy_all <- psy_daily %>%
select(Tree, Logger, date, type, WP) %>%
pivot_wider(names_from = type,
values_from = WP) %>%
left_join(met_daily, by = "date")
# save(psy_all, file = "../daily_app/psy_all.Rdata")
maint <- read_csv("../data_raw/Maintenance-times.csv") %>%
mutate(st_date = as.Date(Start),
en_date = as.Date(End))
# save(maint, file = "../daily_app/maint.Rdata")
For each tree and logger, we examine how stem water potential varies over time with daily maximum VPD and precipitation. Vertical bars represent maintenance periods.
Next, plot the predawn and midday water potentials as a time series with soil moisture.
Because the loggers were rotated to a fresh branch during each maintenance period, the label of Logger 1 and Logger 2 are better expressed as Branches 1 - x. Let us add the Branch label to each set of water potential observations.
psy_all2 <- psy_all %>%
mutate(period = case_when(date > maint$en_date[1] &
date < maint$st_date[2] ~ 1,
date > maint$en_date[2] &
date < maint$st_date[3] ~ 2,
date > maint$en_date[3] &
date < maint$st_date[4] ~ 3,
date > maint$en_date[4] &
date < maint$st_date[5] ~ 4,
date > maint$en_date[5] ~ 5))
branch_id <- psy_all2 %>%
group_by(Tree, Logger, period) %>%
count() %>%
group_by(Logger, period) %>%
mutate(branch_fixed = cur_group_id()) %>%
select(-n)
psy_all3 <- psy_all2 %>%
left_join(branch_id)
Next, plot relationships between environmental variables and predawn/midday water potential, with individual branches as separate colors.
First, the relationship between daily max VPD and stem water potential differs between trees, but also between branches.
Second, the relationship between soil VWC and stem water potential also shows both between tree and between branch differences.