Species-specific growth adjustment - Nutrients

The nutrient stress growth reducer $N U T_{t x y s}$ [-] is defined as:

N U T_{t x y s} = max (N U T_{A M C, t x y s}, N U T_{R S A, t x y s})

The growth reducers based on arbuscular mycorrhizal colonisation rate $N U T_{A M C, t x y s}$ and the root surface area $N U T_{R S A, t x y s}$ are described below. The maximum of both response curves is used for the nutrient reduction function. It is assumed that the plants needs either many fine roots per above ground biomass or have a strong symbiosis with mycorrhizal fungi. Both functions use the calculation of the plant available nutrients.

Growth reducers

the nutrient stress growth reducer based on the arbuscular mycorrhizal colonisation rate $N U T_{A M C, t x y s}$ [-] is defined as:

\begin{aligned} N U T_{A M C, t x y s} & = {\begin{cases} 0 & if R = 0 \\ 1 / (1 + \exp (- β_{N U T, a m c} \cdot (N_{p, t x y s} - x_{0, N, A M C}))) & if 0 < R < 1 \\ 1 & if R >= 1 \end{cases} \\ x_{0, N, A M C} & = \frac{1}{β_{N U T, a m c}} \cdot (- δ_{N U T, a m c} \cdot (T A M C_{t x y s} - (\frac{1}{δ_{N U T, a m c}} \cdot \log (\frac{1 - α_{N U T, a m c, 05}}{α_{N U T, a m c, 05}}) + ϕ_{T A M C}))) + 0.5 \\ T A M C_{t x y s} & = \frac{B_{B, t x y s}}{B_{t x y s}} \cdot a m c_{s} \end{aligned}

the nutrient stress growth reducer based on the root surface area $N U T_{R S A, t x y s}$ [-] is defined as:

\begin{aligned} N U T_{R S A, t x y s} & = {\begin{cases} 0 & if R = 0 \\ 1 / (1 + \exp (- β_{N U T, r s a} \cdot (N_{p, t x y s} - x_{0, N, R S A}))) & if 0 < R < 1 \\ 1 & if R >= 1 \end{cases} \\ x_{0, N, R S A} & = \frac{1}{β_{N U T, r s a}} \cdot (- δ_{N U T, r s a} \cdot (T R S A_{t x y s} - (\frac{1}{δ_{N U T, r s a}} \cdot \log (\frac{1 - α_{N U T, r s a, 05}}{α_{N U T, r s a, 05}}) + ϕ_{T R S A}))) + 0.5 \\ T R S A_{t x y s} & = \frac{B_{B, t x y s}}{B_{t x y s}} \cdot r s a_{s} \end{aligned}

$ϕ_{T A M C}$ reference trait value [-]
$β_{N U T, a m c}$ slope of response function [-]
$α_{N U T, a m c, 05}$ response at $N_{p, t x y s} = 0.5$ for species with the reference trait value [-]
$δ_{N U T, a m c}$ scales the difference in the growth reducer between species [-]
$ϕ_{T R S A}$ reference trait value [m² g⁻¹]
$β_{N U T, r s a}$ slope of response function [-]
$α_{N U T, r s a, 05}$ response at $N_{p, t x y s} = 0.5$ for species with the reference trait value [-]
$δ_{N U T, r s a}$ scales the difference in the growth reducer between species [g m⁻²]

Visualization

growth reducer based on root surface area per total biomass:

response at Np = 0.5 for species with the reference trait value α_NUT_rsa_05 (strong to weak growth reduction)	0.9
difference between species δ_NUT_rsa (no to strong difference)	10
slope of response β_NUT_rsa	7
reference trait value ϕ_TRSA	0.15

growth reducer based on arbuscular mycorrhizal colonisation rate per total biomass:

response at Np = 0.5 for species with the reference trait value α_NUT_amc_05 (strong to weak growth reduction)	0.9
difference between species δ_NUT_amc (no to strong difference)	10
slope of response β_NUT_amc	7
reference trait value ϕ_TAMC	0.2

API

GrasslandTraitSim.nutrient_reduction! Function

julia

nutrient_reduction!(; container, nutrients, total_biomass)

Reduction of growth based on plant available nutrients and the traits arbuscular mycorrhizal colonisation and root surface area per belowground biomass.

source

Plant available nutrients

The plant available nutrients $N_{p, t x y s}$ [-] are descriped by:

\begin{aligned} N_{p, t x y s} & = \frac{N_{x y}}{α_{N U T, N m a x}} \cdot N U T_{a d j, t x y s} \\ N U T_{a d j, t x y s} & = α_{N U T, m a x a d j} \cdot \exp (\frac{\log (\frac{1}{α_{N U T, m a x a d j}}) \cdot \sum_{i = 1}^{S} T S_{s, i} \cdot B_{t x y i}}{α_{N U T, T S B}}) \\ TS & = [\begin{array}{c} 1 & T S_{1, 2} & \dots & T S_{1, S} \\ T S_{2, 1} & 1 \\ ⋮ & ⋱ \\ T S_{S, 1} & 1 \end{array}] \\ TS & = 1 - \frac{TD}{max (TD)} \\ T D_{s, i} & = \sqrt{{(R S A_{n o r m, s} - R S A_{n o r m, i})}^{2} + {(A M C_{n o r m, s} - A M C_{n o r m, i})}^{2}} \\ R S A_{norm, s} & = \frac{r s a_{s} - mean (R S A)}{sd (R S A)} \\ A M C_{norm, s} & = \frac{a m c_{s} - mean (A M C)}{sd (A M C)} \end{aligned}

$α_{N U T, N m a x}$ maximum total soil nitrogen [g kg⁻¹]
$α_{N U T, m a x a d j}$ maximum nutrient adjustment factor [-]
$α_{N U T, T S B}$ reference value, at $α_{N U T, T S B} = T S_{s, i} \cdot B_{t x y i}$ is the nutrient adjustment factor $N U T_{a d j, t x y s} = 1$ [kg ha⁻¹]

Visualization

maximum nutrient adjustment factor α_NUT_maxadj
reference value for ∑ TS ⋅ B α_NUT_TSB

API

GrasslandTraitSim.similarity_matrix! Function

julia

similarity_matrix!(; container)

Calculates the similarity between plants concerning their investment in fine roots and collaboration with mycorrhiza.

source

GrasslandTraitSim.nutrient_competition! Function

julia

nutrient_competition!(; container, total_biomass)

Models the density-dependent competiton for nutrients between plants.

source

GrasslandTraitSim.input_nutrients! Function

julia

input_nutrients!(; container)

Scale the total soil nitrogen

source

Species-specific growth adjustment - Nutrients ​

Growth reducers ​

Visualization ​

API ​

Plant available nutrients ​

Visualization ​

API ​

Species-specific growth adjustment - Nutrients

Growth reducers

Visualization

API

Plant available nutrients

Visualization

API