Plant Models
NicheMapR can be used to model leaf temperature and water loss explicitly, and root temperature implicitly. The leaf temperature calculations can be linked to photosynthesis with via the plantecophys package. The microclimate model provides many of the environmental inputs required for plant modelling, including soil temperature and moisture, and wavelength-specific irradiance (parameter 'lamb' set to 1 returns the wavelength-specific solar radiation output).
Also, the soil moisture algorithm incorporates a soil-plant-atmosphere-continuum (SPAC) model that is defined by the following plant trait parameters (with default values shown):
- L = c(0, 0, 8.2, 8.0, 7.8, 7.4, 7.1, 6.4, 5.8, 4.8, 4.0, 1.8, 0.9, 0.6, 0.8, 0.4, 0.4, 0, 0) * 10000, root density (m/m3), (19 values descending through soil for specified soil nodes in parameter
- R1 = 0.001, root radius, m
- RW = 2.5e+10, resistance per unit length of root, m3 kg-1 s-1
- RL = 2e+6, resistance per unit length of leaf, m3 kg-1 s-1
- PC = -1500, critical leaf water potential for stomatal closure, J kg-1
- SP = 10, stability parameter for stomatal closure equation, -
- LAI = 0.1, leaf area index, used to partition traspiration/evaporation from potential evapotranspiration (PET)
The microclimate model produces estimates of plant transpiration rate and the leaf and (depth-specific) root water potential from the from the soil moisture calculations (output tables 'plant' and 'shadplant').
The function plantgro takes soil water potential and soil temperature input and computes a crude growing degree day model of potential plant growth and an estimate of plant water content, given threshold values of soil water potential at which the wilting point and permanent wilting point occurs, the depth range of the roots and a threshold root temperature for growth.
A Dynamic Energy Budget plant model is also under development, with a working version that links to the microclimate model outputs available under Julia programming environment (see here).