Chapter 14 Day-level results

In the preceeding chapters we detailed how the model performs canopy and soil energy balances for subdaily time steps and how transpiration and photosynthesis values are determined for any given substep. This chapter indicates how these are aggregated at the daily scale and how other day-level model outputs are calculated.

14.1 Photosynthesis and transpiration

Cohort’s transpiration Trt (eq. (12.1)) are added across subdaily steps to yield daily transpiration (Tr, in mmH2O): Tr=ntt=1Trt An the same for water extraction Exs,t for each soil layer s (eq.(12.2)): Exs=ntt=1Exs,t Exs are substracted from the water content of the corresponding soil layer, closing the soil water balance of the day (eq. (7.1)).

Daily values of gross and net carbon assimilation for plant cohorts are obtained similarly. Ag,t and An,t (eq. (12.3)) are added across subdaily steps to obtain the daily gross and net assimilation values at the cohort level (in gC·m2).

14.2 Plant water potentials and relative water contents

Because the model determines optimum transpiration for every subdaily time step, this leads to a daily sequence of leaf water potential (Ψleaf,t), stem water potential (Ψstem,t), root crown water potential (Ψrootcrown,t) and root surface water potential (Ψrootcrown,s,t) values (for soil layer s in the last case). The model defines the following daily water potentials:

  • Pre-dawn leaf water potential (Ψpd): the maximum of Ψleaf,t values.
  • Pre-dawn shade leaf water potential (Ψshadepd): the maximum of Ψshadeleaf,t values.
  • Pre-dawn sunlit leaf water potential (Ψsunlitpd): the maximum of Ψsunlitleaf,t values.
  • Mid-day leaf water potential (Ψmd): the minimum of Ψleaf,t values.
  • Mid-day shade leaf water potential (Ψshademd): the minimum of Ψshadeleaf,t values.
  • Mid-day sunlit leaf water potential (Ψsunlitmd): the minimum of Ψsunlitleaf,t values.
  • Stem water potential (Ψstem): the minimum of Ψstem,t values.
  • Root-crown water potential (Ψrootcrown): the minimum of Ψrootcrown,t values.
  • Root surface water potentials (Ψrootsurf,s): the minimum of Ψrootsurf,s,t values for each soil layer s.

Analogously, relative water content of stems and leaves is known for every subdaily time step, which results in a daily sequence of leaf relative water content (RWCleaf,t) and stem relative water content (RWCstem,t). These are summarized at the daily level:

  • Leaf relative water content (RWCleaf): the mean RWCleaf,t values.
  • Stem relative water content (RWCstem): the mean of RWCstem,t values.

Finally, the daily sequence of slopes of the supply function (dE/dΨt) is also averaged at the daily level:

  • Slope of the supply function (dE/dΨ): the mean dE/dΨt values.

14.3 Plant drought stress

In order to have an estimate of daily drought stress for the plant cohort, the model uses the complement of relative whole-plant hydraulic conductance. In turn, the relative conductance is determined by dividing the current whole-plant hydraulic conductance by its maximum value:

DDS=ϕ(1.0kplantkmax where \phi is the leaf phenological status, k_{plant} is the current whole-plant hydraulic conductance and k_{\max plant} is the maximum whole-plant conductance. Since k_{plant} varies through the day, average daily values are used for the determination of drought stress. Furthermore, remember than in the Sperry transpiration model we have that the whole-plant conductance is the slope of the vulnerability curve: \begin{equation} k_{plant} = dE/d\Psi \end{equation} In the Sureau sub-model, k_{plant} is estimated from the resistances across the hydraulic network (see 10.5.2).

14.4 Cavitation and conduit refilling

Like with the basic water balance model, the advanced water balance model is normally run assuming that although soil drought may reduce transpiration. Embolized xylem conduits may be automatically refilled when soil moisture recovers if stemCavitationRecovery = "total". Automatic refilling is always assumed for root segments, but not for leaves and stems. There are three other options to deal with leaf and stem xylem cavitation besides automatic refilling. Any of them cause the model to keep track of the proportion of lost conductance in the leaf and stem xylem of the plant cohort (PLC_{leaf} and PLC_{stem}) at successive time steps: \begin{eqnarray} PLC_{leaf, t} &=& \max \{PLC_{leaf, t-1}, 1 - k_{leaf}(\Psi_{leaf})/k_{max,leaf} \} \\ PLC_{stem, t} &=& \max \{PLC_{stem, t-1}, 1 - k_{stem}(\Psi_{stem})/k_{max,stem} \} \end{eqnarray}

In the Sperry sub-model, PLC values are evaluated every subdaily (hourly) time step, whereas in the SurEau sub-model it PLC values are updated more frequently.

For simulations of less than one year one can use stemCavitationRecovery = "none" to keep track of the maximum drought stress. However, for simulations of several years, it is normally advisable to allow recovery. If stemCavitationRecovery = "annual", PLC_{leaf} and PLC_{stem} values are set to zero at the beginning of each year, assuming that embolized plants overcome the conductance loss by creating new xylem tissue. Finally, if stemCavitationRecovery = "rate" the model simulates stem refilling at daily steps as a function of symplasmic water potential. First, a daily recovery rate (r_{refill}; in cm^2 \cdot m^{-2} \cdot day^{-1}) is estimated as a function of \Psi_{symp,stem}: \begin{equation} r_{refill}(\Psi_{symp,stem}) = r_{\max,refill} \cdot \max \{0.0, (\Psi_{symp,stem} + 1.5)/1.5\} \end{equation} Where r_{\max,refill} is the control parameter cavitationRecoveryMaximumRate indicating a maximum refill rate. The right part of the equation normalizes the water potential, so that r_{refill} = r_{refill,\max} if \Psi_{symp,stem} = 0 and r_{refill} = 0 if \Psi_{symp,stem} <= -1.5MPa. The proportion of conductance lost in leaves and stem are then updated using: \begin{eqnarray} PLC_{leaf} &=& \max \{0.0, PLC_{leaf} - (r_{refill}/H_v) \} \\ PLC_{stem} &=& \max \{0.0, PLC_{stem} - (r_{refill}/H_v) \} \end{eqnarray} where H_v is the Huber value in units of cm^2 \cdot m^{-2}.