Publications by Beni Stocker
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
Document
Code for this is here. library(tidyverse) ## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.2 ── ## ✔ ggplot2 3.3.6 ✔ purrr 0.3.4 ## ✔ tibble 3.1.8 ✔ dplyr 1.0.10 ## ✔ tidyr 1.2.0 ✔ stringr 1.4.1 ## ...
980 sym R (4854 sym/13 pcs) 4 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(tidyr) library(ggplot2) library(rgdal) # for readOGR ## Loading required package: sp ## Please note that rg...
350 sym R (3302 sym/11 pcs) 2 img
Document
Net radiation as forcing? Check whether sufficient sites have net radiation data. filnam <- "~/rsofun/data/df_fluxnet.rds" if (!file.exists(filnam)){ df_fluxnet <- suppressWarnings( suppressMessages( ingestr::ingest( siteinfo = ingestr::siteinfo_fluxnet2015, source = "fluxnet", getvars = lis...
261 sym Python (1549 sym/6 pcs) 2 img 2 tbl
Document
Net radiation as forcing Read data Check whether sufficient sites have net radiation data. filnam <- "./data/df_fluxnet.rds" if (!file.exists(filnam)){ df_fluxnet <- ingest( siteinfo = ingestr::siteinfo_fluxnet2015, source = "fluxnet", getvars = list( netrad = "NETRAD", ppfd = "SW_IN_F", temp = "TA_F_DAY" ...
1921 sym Python (8554 sym/26 pcs) 10 img 1 tbl