Phenology
Description and equations for Knorr module
Overview
The leaf area index (LAI) phenology model is based on the scheme developed by Knorr et al (2010) with updates by Norton et al (2023). The Knorr et al (2010) model brings together two concepts. First, the transition of LAI between dormant, active, and senescent states is controlled by three potentially limiting environmental constraints: Temperature, photoperiod, and water availability. Second, spatial variability in growth triggers within a population of plants results in smooth, differentiable functions that describe the transition between dormant, active, and senescent states. This provides a more realistic representation of LAI dynamics, as opposed to unrealistic and impractical step functions. Norton et al (2023) described the coupling of this LAI phenology model to the carbon and water balance in DALEC and implemented it in CARDAMOM, thereby incorporating another potentially constraint of carbon supply. This model demonstrated its skill in capturing LAI and carbon cycle dynamics across diverse ecosystem types. We note that in Norton et al (2023) the LAI phenology sub-model used evapotranspiration as an input to compute the water-limited LAI, considering evaporation and transpiration were not modeled separately. Here, evaporation and transpiration are modeled separately, therefore transpiration is used as the input which is consistent with the original Knorr et al (2010) formulation and adds realism to this version of CARDAMOM.
Background
The LAI phenology model was originally developed by describing a generic differential equation in time for the LAI of an individual plant, which is then integrated in space to represent LAI dynamics of a population of plants. The time evolution of LAI for an individual plants is commonly described using triggers for growth and senescence (e.g. a growing degree day threshold), which are represented by conditional (step) functions. The Knorr et al (2010) model makes the useful development by assuming the triggers for a population of plants can be represented by a Gaussian probability density function. As described below, this enables the time evolution of LAI for the population to be smooth and differentiable.
In the original formulation, which is also implemented into DALEC, there are two triggers considered: temperature (
Where
Where
The above equation can be simplified by representing each integral with a cumulative normal distribution,
Where
Growth and Senescence
The LAI growth function,
Where
Where
Water Limitation
To incorporate water availability constraints on LAI dynamics, a water-limited LAI is defined using the balance of water availability (soil moisture) and water demand (loss via transpiration) as follows:
Where
Maximum Potential Leaf Area Index
The
In practice, the minimum is calculated using a quadratically smoothed minimum function. Additionally, in order for water limitation to account for changes in water availability in the recent past and not just instantaneous changes in
Where
Phenology-Determining Temperature
Typically, the phenology-determining temperature follows the growing-degree days concept. However, for reasons outlined in Knorr et al (2010), the phenology-determining temperature,
Where
Coupling to the Carbon Balance
The equations described above govern how environmental conditions directly constrain the LAI dynamics over time, including temperature, day length and water availability. This must be coupled to carbon balance, whereby a third potentially limiting factor can control LAI dynamics: plant carbon supply.
Note: I think this section may actually be better placed into the carbon allocation and autotrophic respiration section. Makes more sense to me.
References
Knorr, W., T. Kaminski, M. Scholze, N. Gobron, B. Pinty, R. Giering, and P.-P. Mathieu, 2010, Carbon cycle data assimilation with a generic phenology model, J. Geophys. Res., 115, G04017, https://doi.org/10.1029/2009JG001119
Norton, A. J., Bloom, A. A., Parazoo, N. C., Levine, P. A., Ma, S., Braghiere, R. K., and Smallman, T. L., 2023, Improved process representation of leaf phenology significantly shifts climate sensitivity of ecosystem carbon balance, Biogeosciences, 20, 2455–2484, https://doi.org/10.5194/bg-20-2455-2023