R/transp_maximumTranspirationModel.R
transp_maximumTranspirationModel.Rd
Builds a model of maximum transpiration (Tmax) over potential evapotranspiration (PET) for increasing leaf area index (LAI) values for each plant cohort.
An object of class spwbInput
, built using the 'Sperry' transpiration mode.
A data frame with daily meteorological data series.
Latitude (in degrees).
Elevation above sea level (in m), slope (in degrees) and aspect (in degrees from North).
Sequence of stand LAI values to be tested.
Logical flag to indicate plotting of results.
Returns a list with as many elements as plant cohorts, each element being a glm
model.
This function performs a metamodelling exercise using the Sperry transpiration model, with the aim to estimate coefficients
for the equation used in the Granier transpiration model (Granier et al. 1999). The model to be fitted is: y ~ a*LAI + b*LAI^2
,
where y
is the ratio between maximum transpiration (Tmax) and Penman's potential evapotranspiration (PET) and LAI
is the stand LAI.
Unlike the original equation of Granier et al. (1999), we fit a zero intercept model so that LAI = 0 translates into zero plant transpiration.
The function fits the model for each cohort separately, assuming it represents the whole stand.
For each stand LAI value in the input sequence, the function uses simulations with Sperry transpiration and the input weather to estimate y = Tmax/PET
as a function of stand's LAI (deciduous stands include leaf phenology).
Once simulations have been conducted for each stand LAI value, the function fits a Generalized Linear Model with the above equation,
assuming a Gamma distribution of residuals and an identity link.
The coefficients of the model can be used to parametrize Granier's transpiration,
since coefficients a
and b
in the equation above correspond to parameters Tmax_LAI
and Tmax_LAIsq
,
respectively (see SpParamsMED
).
Granier A, Bréda N, Biron P, Villette S (1999) A lumped water balance model to evaluate duration and intensity of drought constraints in forest stands. Ecol Modell 116:269–283. https://doi.org/10.1016/S03043800(98)002051.
# \donttest{
#Load example daily meteorological data
data(examplemeteo)
# Load example plot plant data
data(exampleforest)
# Load default species parameters
data(SpParamsMED)
# Initialize soil with default soil params
examplesoil < soil(defaultSoilParams(4))
# Initialize control parameters for 'Sperry' transpiration mode
control < defaultControl(transpirationMode="Sperry")
# Initialize input
x2 < forest2spwbInput(exampleforest,examplesoil, SpParamsMED, control)
# Estimate maximum transpiration ratio models for each cohort
# Weather is subset to speedup results
m < transp_maximumTranspirationModel(x2, examplemeteo[1:10,],
41.82592, elevation = 100,
slope = 0, aspect = 0)
#>

  0%

===  5%

======  9%

==========  14%

=============  18%

================  23%

===================  27%

======================  32%

=========================  36%

=============================  41%

================================  45%

===================================  50%

======================================  55%

=========================================  59%

=============================================  64%

================================================  68%

===================================================  73%

======================================================  77%

=========================================================  82%

============================================================  86%

================================================================  91%

===================================================================  95%

====================================================================== 100%
# Inspect the model for first cohort
m[[1]]
#> LAI I(LAI^2)
#> 0.183385509 0.008626493
# }