Publications by Beni Stocker
Document
Read data This loads data from FluxDataKit. df <- read_rds("~/data/FluxDataKit/rsofun_driver_data_clean.rds") df_fr_pue <- df |> filter(sitename == "FR-Pue") |> select(forcing) |> unnest(forcing) Plot GPP in years 2013 and 2014 is ominously low. ggplot(data = df_fr_pue) + geom_line( aes( date, gpp ), colour = "b...
387 sym 2 img
Publish Document
Code for this notebook is here. Model setup The P-hydro model (Joshi et al., 2022) is used here to illustrate the effect of acclimating P50 (leaf water potential at 50% loss of conductivity) and K (whole plant conductivity) on the response of stomatal conductance and transpiration to VPD and soil moisture. # other environmental factors held consta...
738 sym 1 img
Publish Document
Model setup The P-hydro model (Joshi et al., 2022) is used here to illustrate the effect of acclimating P50 (leaf water potential at 50% loss of conductivity) and K (whole plant conductivity) on the response of stomatal conductance and transpiration to VPD and soil moisture. # other environmental factors held constant kphio = 0.087; # quantum...
649 sym 2 img
Document
Model setup The P-hydro model (Joshi et al., 2022) is used here to illustrate the effect of acclimating P50 (leaf water potential at 50% loss of conductivity) and K (whole plant conductivity) on the response of stomatal conductance and transpiration to VPD and soil moisture. # other environmental factors held constant kphio = 0.087; # quantum...
646 sym 1 img
Document
Land C cycle is a 1-box model. Steady-state in 1958. The sensitivity of the C-input flux to CO2 follows a logarithmic relationship to the CO2-change relative to 1958. A linear multiplication factor is applied to the sensitivity. Then, diagnose the sink flux in years 1980-2012. 1-box ## ////////////////////////////////////////////////// ## 1-BOX MO...
296 sym R (3113 sym/6 pcs) 1 img
Document
This is based on the script shared by Bernhard Schmid and contained also this repo (analysis/TianDiMay23_BS.R). Read and pre-process data Read full species-level data df <- read_csv("~/data/LeafNP_tiandi/Global_total_leaf_N_P_Di/leafnp_data_covariates_20210702.csv") %>% mutate(grass = tree_shrub_Herb == "H") ## Rows: 36414 Columns: 66 ## ─�...
960 sym 1 tbl
Publish Document
It is sometimes stated that “global carbon and water cycles have shifted from a CO2-dominated era into a VPD-dominated one”. The VPD and CO2 effects on photosynthesis are interactive. Therefore, is this statement a good description of the nature of this interaction? Let’s look at this by comparing the CO2 effect on GPP under ambient and eleva...
1686 sym R (3800 sym/12 pcs) 3 img
Document
It is sometimes stated that “global carbon and water cycles have shifted from a CO2-dominated era into a VPD-dominated one”. The VPD and CO2 effects on photosynthesis are interactive. Therefore, is this statement a good description of the nature of this interaction? Let’s look at this by comparing the CO2 effect on GPP under ambient and eleva...
1364 sym R (2066 sym/9 pcs) 1 img
Document
ftemp_kphio <- function(temp, kphio, kphio_par_a, kphio_par_b){ out <- kphio * (1.0 + kphio_par_a * (temp - kphio_par_b)^2) out <- ifelse(out < 0, 0, ifelse(out > 1, 1, out)) } library(ggplot2) ggplot() + geom_function(fun = function(x) ftemp_kphio(x, kphio = ...
6 sym R (1192 sym/2 pcs) 1 img
Document
library(dplyr) ## ## Attaching package: 'dplyr' ## The following objects are masked from 'package:stats': ## ## filter, lag ## The following objects are masked from 'package:base': ## ## intersect, setdiff, setequal, union library(ggplot2) library(readr) library(tidyr) Read data df_temp <- read_csv("~/data/n2o_Yunke/final_obs_dataset/ob...
1951 sym R (3482 sym/25 pcs) 2 img 1 tbl