Interception capacity with snow interception

If snow interception is included in the simulation, the interception capacity, Simax, can be calculated in two different ways, either as a function of thermal quality or as a function of air temperature (see switch SnowIntUnload). In the latter case interception capacity is calculated as:

                       (3.56)

where iLAIsnow, iLAI and ibase are parameters, and Al is the leaf area index. In this case, thermal quality QI* is assumed to be equal to the thermal quality of precipitation, QP calculated in eq.(4.36).

When the interception capacity is a function of thermal quality, it is instead calculated as:

                                (3.57)

QI* is the thermal quality (fraction of frozen water) of the intercepted water and can either be calculated as a weighted sum of the thermal quality of the intercepted water from the previous time-step, QI, and the thermal quality of new precipitation. Thermal quality is calculated as:

                                          (3.58)

where QP is the thermal quality of precipitation calculated in eq.(4.36). fnew is the fraction of new intercepted precipitation in relation to total intercepted storage:

                                                             (3.59)

where P is precipitation and Si is the interception storage.

When the interception storage, Si, has been calculated in each time-step, a new value on thermal quality of intercepted water, QI, is calculated:

                                                       (3.60)

where QI* is the thermal quality of intercepted water calculated in the beginning of the time-step and Si is the interception storage. The amount of melted intercepted storage, Smelt, is estimated by:

                                               (3.61)

where iscale is a parameter and M(Ssnowthick) is the function for calculating snow melt, eq.(4.32)-(4.34). Ssnowthick replaces Δzsnow and is calculated as:

                            (3.62)

where Sint is the interception storage, ρwater is the density of water, QI* is the thermal quality of intercepted water calculated in the beginning of the time-step and Al is the leaf area index. The figure 100 in the equation is an approximation of the snow density.