Dynamic vegetation processes

Maximum canopy storage

SPHY allows the user to use the dynamic vegetation module in order to incorporate a time-variable vegetation cover and corresponding rainfall interception. In order to calculate the rainfall interception, the canopy storage needs to be calculated, using a time series of NDVI (Carlson and Ripley 1997). The first step involves the calculation of the fraction photosynthetically active radiation (FPAR). FPAR can be calculated using a relation between NDVI and FPAR, which was found by Peng et al. (2012) and described by Sellers et al. (1996), according to

Equation 4

FPAR=min((SRāˆ’SRmin)(FPARmaxāˆ’FPARmin)SRmaxāˆ’SRmin+FPARmin,0.95)FPAR = min(\frac{(SR-SR_{min})(FPAR_{max}-FPAR_{min})}{SR_{max}-SR_{min}}+FPAR_{min},0.95)

with

Equation 5

SR=1+NVDI1āˆ’NVDISR = \frac{1+NVDI}{1-NVDI}

and FPARmaxFPAR_{max}(-) and FPARminFPAR_{min}(-) having values of 0.95 and 0.001, respectively. An FPAR of 0.95 is equivalent to the maximum LAI for a particular class, and an FPAR of 0.001 is equivalent to a minimum LAI. In order to calculate FPAR, an NDVI time series is required.

The second step is the calculation of the leaf-area index (LAI), which is eventually required to calculate the maximum canopy storage (Scanmax)(Scan_{max}). According to Monteith (1973), LAI for vegetation that is evenly distributed over a surface can be calculated using a logarithmic relation between LAI and FPAR, according to

Equation 6

LAI=LAImaxāˆ—log(1āˆ’FPAR)log(1āˆ’FPARmax)LAI = LAI_{max} *\frac{log(1-FPAR)}{log(1-FPAR_{max})}

with LAI (ā€“) the leaf-area index, and LAImaxLAI_{max} (-) the maximum leaf-area index (vegetation type dependent). This means that the maximum and minimum LAI values are related to the maximum and minimum of FPAR. Table 2 shows the LAImaxLAI_{max}values for a certain number of vegetation types.

Vegetation type

LAImaxLAI_{max}[-]

Broadleaf evergreen trees

7

Broadleaf deciduous trees

7

Mixed trees

7.5

Needleleaf evergreen trees

8

High latitude deciduous trees

8

Grass with 10 - 40% woody cover

5

Grass with <10% woody cover

5

Shrubs and bare soil

5

Moss and lichens

5

Bare

5

Cultivated

6

Table 2: LAImaxLAI_{max} values for different vegetation types (Sellers et al., 1996).

For vegetation that is concentrated in clusters, the linear relation from Goward and Huemmrich (1992) is often used. However, since SPHY is generally applied using grid-cell resolutions between 250m and 1km, we can assume that the effect of having vegetation concentrated in clusters is negligible. Therefore, the calculation of LAI in SPHY is done using the logarithmic relation of Monteith (1973) (Equation 6).

The next step involves the calculation of the maximum canopy storage (Scanmax(mm))(Scan_{max} (mm)). Many different relations between ScanmaxScan_{max} and LAI can be found in the literature, depending on the vegetation type (Jong and Jetten 2010). The best results for crop canopies are shown by Kozak et al. (2007) and are archived by Von Hoyningen-Huene (1981), who derived the following relation between ScanmaxScan_{max} and LAI:

Equation 7

Scanmax=0.935+0.498LAIāˆ’0.00575LAI2Scan_{max} = 0.935 + 0.498LAI-0.00575LAI^2

Interception

Interception is calculated on a daily basis if the dynamic vegetation module is used, and consists of the daily precipitation plus the intercepted water remaining in the canopy storage from the previous day. First, the canopy storage is updated with the amount of precipitation of the current day:

Equation 8

Scant=Scan(tāˆ’1)+PtScan_t=Scan(t-1)+P_t

with Scant(mm)Scan_{t} (mm) the canopy storage on day tt, Scantāˆ’1(mm)Scan_{t-1} (mm) the canopy storage on day tāˆ’1t-1, and Pt(mm)P_{t} (mm) the amount of precipitation on day tt. The portion of precipitation that cannot be stored in the canopy storage is known as precipitation throughfall, or effective precipitation, according to:

Equation 9

Pet=max(0,Scantāˆ’Scanmax,t)Pe_t=max(0,Scan_t-Scan_{max,t})

with Pet(mm)Pe_{t} (mm) the effective precipitation on day tt, and Scant(mm)Scan_{t} (mm) the canopy storage on day tt. This equation shows that precipitation throughfall only occurs if the water stored in the canopy exceeds the maximum canopy storage. After the effective precipitation has been calculated, the canopy storage is updated as:

Equation 10

Scant=Scant=PetScan_t=Scan_t=Pe_t

The remaining amount of water stored in the canopy is available for interception, and the amount of water that will be intercepted depends on the atmospheric demand for open water evaporation. A commonly used value for the atmospheric demand for open water evaporation is 1.5 (Allen et al. 1998), which is derived from the ratio between 1 and the mean pan evaporation coefficient Kp(~0.65). The interception can now be calculated using:

Equation 11

Intt=min(1.5ETr,t,Scant)Int_t=min(1.5ET_{r,t},Scan_t)

with Intt(mm)Int_{t} (mm) the intercepted water on day tt, and ETr,t(mm)ET_{r,t} (mm) the reference evapotranspiration on day tt. Finally, the canopy storage is updated by subtracting the interception:

Equation 12

Scant=Scantāˆ’InttScan_t=Scan_t-Int_t

Last updated