1. Introduction
Rainfall interception loss is the proportion of gross rainfall that is intercepted, stored and subsequently evaporated from all parts of vegetation during or following rainfall. Many studies on rainfall interception loss observations and modeling in forest have been carried out and it has been reported that rainfall interception loss totals were between 10% and 40% of annual gross rainfall [
1,
2,
3]. Therefore, studying the interception loss of forest is very important for improving the modeling of regional and global water balances and understanding the terrestrial water cycle processes [
4].
Until now, many researchers have developed a series of rainfall interception models including: the empirical statistical model [
5] and the physically-based model [
6,
7,
8,
9,
10,
11]. The latter has been paid more attention due to its robustness in physics and less need for calibration of empirical coefficients. Almost all of the physically-based models have been derived from the Rutter (1972) model [
7] and the Gash (1979) and (1995) models [
6,
8]. Applications have covered different types of vegetation including rainforest, conifers, mixed conifer, shrubs, crops,
etc. Most of these applications, however, were at field scale [
12].
Satellite remote sensing observations can provide land surface information with high spatial and temporal resolution. Unfortunately, estimate of interception loss at regional scale using remote sensing information is not straight forward. Empirically based models have been applied to remote sensing observations for regional estimate of interception evaporation. Bastiaanssen
et al. [
13] used the classical von Hoyningen algorithm [
14] to calculate the interception evaporation. Such a method has an implication that maximum interception loss will stay below the vegetation storage capacity for single rainfall and will therefore underestimate interception loss for large rainfalls. Mu
et al. [
15] considered evaporation from a wet canopy surface according to the fractional canopy water cover defined by the relative humidity by setting a threshold of relative humidity in priori. On the contrary, physically-based sophisticated models with solid physical meaning are usually too complex so that parameters needed in such models often cannot be obtained from remote sensing observations.
The Gash (1979 and 1995) models have been widely used and verified in weekly and monthly temporal scales at different sites [
8,
16,
17,
18,
19]. For most researchers only paid attention to the long-term interception loss of forest, the applicability of the Gash model in estimating daily interception loss has not been effectively verified. The Gash (1979) model considered rainfall to occur as a series of discrete events, during which three phases can be distinguished: a wetting phase, a saturation phase, and a drying phase after rainfall has ceased. Based on this, the interception loss of canopy is calculated by the canopy storage capacity, canopy coverage, and the ratio of mean evaporation rate from wet canopy over mean rainfall rate. The interception loss of trunk is calculated by the trunk storage capacity and the proportion of the rainfall diverted to stemflow. The main improvement of the Gash (1995) model is that the wet canopy evaporation rate per unit ground area is linearly dependent on the canopy cover fraction. Since then, there have also been many researchers who have proposed significant improvements at field scale over the Gash (1995) model, e.g., van Dijk [
10] adapted the Gash model for vegetation whose characteristics change markedly during the growing season, by relating the canopy capacity and evaporation rate from a saturated canopy to leaf area index (LAI), but was not given much consideration, as reviewed in [
12].
The results of site-based studies have limited representability for regional scale since forest canopy is often shown to be inhomogeneous with large spatial and temporal variability in its stand distribution, LAI and canopy fractional coverage. There are still many difficulties when applying the Gash model at regional scale. One of the most important reasons is that the model calculates the interception loss of canopy and trunk separately, so that more parameters need to be known as a priori but which are usually difficult to extrapolate from field scale to regional scale. In spite of some researchers suggesting that the canopy storage capacity is linearly linked with LAI [
10,
20], it is still difficult to obtain parameters about trunks at regional scale. In fact, LAI normally refers to the total one-sided area of all green canopy elements per unit ground area [
21], which does not include other parts of canopy, e.g., dead leaves. One case among very rare applications at regional scale is made by Miralles
et al. [
22]. However, they used a single set of parameters derived from statistics from literature for calculation of interception loss over the entire region of their study.
Nowadays, remote sensing can provide many variables or inputs needed by the Gash model, for instance LAI, Fractional Vegetation Cover (FVC). Very few attempts have been made in applying these satellite observed parameters to Gash models, since there are still a few critical parameters which cannot be directly derived from remote sensing observations, e.g., canopy storage capacity, trunk storage capacity and the proportion of the rainfall diverted to stemflow.
The objective of this study is to develop a simple and practical forest interception loss model to satisfy the requirements of interception loss estimate of heterogeneous forest at regional scale. This model is a modified forest interception loss model based on the Gash (1979 and 1995) models, referred to as the RS-Gash model (remote sensing based Gash model), by introducing remote sensing observations of Vegetation Area Index (VAI) and FVC. The RS-Gash model was verified using measurements collected at two forest hydrology experimental sites. The sensitivity analysis of the model was done by varying the inputs one-by-one, to test the robustness of the model. Finally, some conclusions about the RS-Gash model are summarized.
2. Theory of the Gash Model
The Gash model [
6] was widely used initially to estimate interception loss of forest taking into account rainfall rate, forest canopy, trunk and meteorological factors. Improvements were made in a later work applicable to sparse forest canopy [
8].
According to the amount of gross rainfall necessary to saturate the canopy (
) and the amount of gross rainfall necessary to saturate the trunk (
St/
pt), where
St is the trunk storage capacity and
pt is the proportion of the rainfall diverted to stemflow, rainfall can be divided into three categories: (i) the gross rainfall (
PG) is not larger than
so that forest canopy is not saturated, and interception loss is equal to
c∙
PG, where
c is canopy coverage; (ii) the gross rainfall
PG is larger than
but not larger than
St/
pt so that forest canopy is saturated but the trunk is not, and the interception loss is equal to
, where
EC(
mm∙h-1) is the mean evaporation rate per unit canopy cover when the forest canopy is saturated;
R(
mm∙h-1) is the mean rainfall rate for saturated canopy condition, which includes the canopy storage, evaporation from the wetting canopy and the part of rainfall diverted into trunks; (iii) the gross rainfall
PG is larger than
St/
pt so that both forest canopy and trunk are saturated, and the interception loss is equal to
, which includes the canopy storage, trunk storage, and evaporation during rainfall, which implies that the evaporation from wet trunk was neglected. Then, the interception loss of canopy and trunk were calculated separately. So in the Gash analytical model, interception loss of forest consists of interception loss of canopy and interception loss of trunk.
It is generally assumed that there is one rainfall event per day so that the Gash model can be used to calculate the interception loss on a daily scale. For a single rainfall, the amount of gross rainfall necessary to saturate the canopy
is given by:
where,
S is the storage capacity of canopy per unit area of ground and
c is the canopy cover. In the original Gash (1995) model, the mean evaporation rate
EC and mean rainfall rate
R are obtained in the whole experimental period to replace each four-week period as in the Gash (1979) model and it is assumed that the ratio
EC over
R is equal for all storms to still get satisfactory results.
The interception loss of canopy (
Ic) can be calculated by:
The interception loss of trunk (
It) can be calculated by:
Since the Gash model calculates the interception loss of canopy and trunk as two parts separately, it needs a large number of parameters, such as EC, R, c, Sc, St and pt. When applying the Gash model at regional scale, many parameters such as St and pt are difficult to obtain.
5. Results and Analysis
5.1. Parameters for the RS-Gash Model
Firstly, the RS-Gash model was verified using the filed measured rainfall and the parameters based on satellite as described below. After that, an interception loss map was created using the rainfall data from TRMM and the same parameters. Fraction of forest (FF) of the two experimental sites (0.01° × 0.01°, Dayekou and Pailougou experimental sites), was calculated from the Hyperion image used in the field validation.
The VAI and FVC was derived from MODIS LAI and NDVI products after processing using HANTS (Harmonic Analysis of Time Series) [
33] to get a gap-free time series with temporal resolution of eight days. Equations (9) and (10) were used to calculate the VAI with
Ls,min = 1 and
α = 0.85 for 8-day MODIS LAI products of picea crassifolia forest. Equation (14) was used to calculate the FVC with
NDVIfv = 0.8 and
NDVIbs = 0.125. Considering the seasonal variation of VAI and FVC due to seasonally dependent photosynthetic activity with the model running at daily steps, daily VAI and FVC time series were obtained by linear interpolation of the 8-day time series, based on the assumption that daily VAI and FVC would change linearly between the two products.
For specific vegetation storage
SV, we used the value (
SV = 0.35) of picea crassifolia, which is the same species as our forest study just a different location measured by Lankreijer
et al. [
26] with the regression method, since this parameter was not measured during our experiment. The value is maybe a little different from our study, due to differences in canopy structures such as stem density, canopy height, and leaf area index between the two study areas. However, a small difference in
SV has limited influence on our study, according to the sensitivity analysis (see
Section 5.4) which shows that ±20% error of
SV only results in less than 0.2 mm (±0.16
mm) error. The mean evaporation rate
EV was calculated by Equation (11) using meteorological data of AWS (listed in
Table 2) for the periods when rainfall was larger than 0.5
mm∙
h−1 [
6]. Since the study area (1° × 2.5°) is relative small, the difference in meteorological conditions over this study area can also be considered as relatively small. Consequently, it is acceptable to take a unique value of
EV over the study area. In our study,
EV = 0.2
mm∙h−1. .
In the filed validation, we used the measured rainfall at the station, and the mean rainfall rate R is calculated from these data when the rainfall was larger than 0.5 mm∙h−1. After the validation, an interception loss map using TRMM data was created at regional scale, and R was calculated using TRMM 3B42 data during the period when the rainfall was larger than 0.5 mm∙h−1
To compare the estimated interception loss with the ground measurements, the modelled interception loss was extracted from the two pixels where the two experimental sites were located. The fraction of forest (FF) in each one of the site-pixels (0.01° × 0.01°) is calculated by the maximum likelihood method from the fine resolution Hyperion image (30 m × 30 m). The FF in the Dayekou and Pailugou experimental sites was 0.8 and 0.61, respectively.
5.2. Field Validation
To directly verify the RS-Gash model, and avoid the bias due to the difference between TRMM rainfall data and field measured rainfall, the field measured rainfall data together with the other parameters based on satellite data were used in the field validation. In this study the RS-Gash model was evaluated first on a daily scale (single rainfall), and further for the whole period of observations, from 13 June to 7 October 2008.
The interception loss modelled (
Imodelled) by the RS-Gash model using remote sensing observations is presented in mm per unit area of ground, but the measurements are presented in mm per unit area of forest [
22]. So for the validation study, the FF was used to convert the interception loss from per unit area of ground to per unit area of forest, as:
where,
Icorrect is the interception loss after correction by
FF, which can be used to compare with the measured interception loss directly.
The modelled interception loss for single rainfall is in good agreement with the field measurements at the Dayekou experimental site with a
R2 of 0.91 and RMSE of 0.34
mm∙
h−1 (
Figure 3a), at Pailugou experimental site with a
R2 of 0.82 and RMSE of 0.6
mm∙
h−1 (
Figure 3b).
Figure 4 presents the correlation between the modelled interception loss and the ground observed gross rainfall at the two sites. The interception loss, both modelled and measured, increases with gross rainfall. This shows that the gross rainfall is an important factor for the interception loss. For small/light rainfall, the interception loss is a little underestimated, and the main reason might be that the wind in the forest may change the pathway of raindrops during falling thus increasing the interception loss. However, the model does not consider this. Yet, the wind has little effect on a large rainfall. In all, the high
R2 and the low RMSE demonstrate the competence of the modified model at regional scale.
Figure 3.
Validation of the modelled interception loss of a single rainfall using the RS-Gash model against field measurements for the period from 13 June to 7 October 2008. (A) The Dayekou experimental site; (B) The Pailugou experimental site.
Figure 3.
Validation of the modelled interception loss of a single rainfall using the RS-Gash model against field measurements for the period from 13 June to 7 October 2008. (A) The Dayekou experimental site; (B) The Pailugou experimental site.
Figure 4.
Correlation between the interception loss and the ground observed gross rainfall at the two sites for the period from 13 June to 7 October 2008. (A) The Dayekou experimental site; (B) The Pailugou experimental site.
Figure 4.
Correlation between the interception loss and the ground observed gross rainfall at the two sites for the period from 13 June to 7 October 2008. (A) The Dayekou experimental site; (B) The Pailugou experimental site.
The modelled and measured total interception loss are listed in
Table 3 for the period from 13 June to 7 October 2008. From the comparison done at the two experimental sites, modelled interception loss was underestimated by about 14.2% and 13.6% of the measured interception loss, respectively. Considering the uncertainty of the parameters from the satellite due to the different scale, especially the FVC, VAI and FF, the error is acceptable. Although the RS-Gash model can reasonably estimate the interception loss of the most of single rainfall events, there are still some extreme cases with a high relative error. One of the reasons is that
EV/
R is the average of a period, in the case of deviation from the average conditions (some extreme cases in energy and wind speed) large errors in the input might be propagated to the estimated interception loss. Nevertheless, the model is steady enough for a reasonable estimation with acceptable error for the data used in this study as shown in
Section 5.4. On the other hand, the measured errors should not be overlooked. The measured interception loss is very sensitive to the gross rainfall and throughfall. Because interception loss is measured as the relatively small difference between gross and net rainfall (throughfall and steamflow), even small errors in these measurements can result in high relative errors in interception loss and any error of these could be included in the interception loss [
12], the validation data processing can only partly eliminate the error.
Table 3.
The modelled and measured total interception loss at the two sites.
Table 3.
The modelled and measured total interception loss at the two sites.
Site | Gross rainfall (Pg) (mm) | Modelled | Measured | Relative error (% measured) |
---|
mm | % Pg | mm | % Pg |
---|
Dayekou | 162.8 | 33.8 | 20.8 | 39.4 | 24.2 | 14.2 |
Pailugou | 153.5 | 49.7 | 32.4 | 57.5 | 37.5 | 13.6 |
The variation of the interception loss at a regional scale results from the heterogeneity of parameters, such as: VAI, FVC and the mean rainfall rate R. There are two points needed to be stated when using remote sensing to estimate the interception loss at a regional scale. Firstly, there are several unavoidable random errors in results, especially for the mixed pixels. However, with significant classification, the estimated interception loss can also represent the interception loss of the vegetation type. Secondly, the results of modelled interception loss are based on per unit area of ground, but the measurements are based on per unit area of forest. The former has a direct significance for regional water cycle researches and can meet the demand of most of the researches, so there is no need to convert the former to the latter, except for validation.
5.3. Interception Loss of Forest at Regional Scale
Considering the picea crassifolia forest is the main coniferous forest in the study area and information on finer forest classification is not available for the time being, we took the class of coniferous forest in the land cover map to estimate the regional forest interception loss in the study area. The estimated regional distribution of forest rainfall interception loss from 13 June to 7 October 2008 using the RS-Gash model is presented in
Figure 5c,d together with the regional distribution of FVC and VAI (
Figure 5a,b). Both the absolute values of interception loss (
Figure 5c) and the interception loss as a percentage of the gross rainfall (
Figure 5d) showed great spatial variation. The statistical results of FVC, VAI and interception loss are shown in
Table 4.
Table 4.
The variation of Fractional Vegetation Cover (FVC), Vegetation Area Index (VAI) and interception loss (in mm and %) in the spatial from 13 June to 7 October 2008.
Table 4.
The variation of Fractional Vegetation Cover (FVC), Vegetation Area Index (VAI) and interception loss (in mm and %) in the spatial from 13 June to 7 October 2008.
Variables | Mean | Standard deviation |
---|
FVC | 0.47 | 0.16 |
VAI | 1.94 | 0.64 |
Interception loss (mm) | 61.1 | 26.1 |
Interception loss (%) | 20.05 | 6.4 |
The spatial variation of interception loss is clearly attributed to the varying density of the forest as expressed by FVC and VAI in space in addition to rainfall spatial variation. Larger errors might be expected in mixing pixels, in particular at the edge of a large forest area, due to the fact that non-forest pixels might be mis-classed as forest leading to underestimation of the modeled interception loss, and vice versa.
Figure 5.
The spatial distribution of (A) the Fractional Vegetation Cover (FVC); (B) the Vegetation Area Index (VAI); (C) the interception loss (mm); and (D) the interception loss in percentage of gross rainfall (%). Values in the maps are the average over the period from 13 June to 7 October 2008.
Figure 5.
The spatial distribution of (A) the Fractional Vegetation Cover (FVC); (B) the Vegetation Area Index (VAI); (C) the interception loss (mm); and (D) the interception loss in percentage of gross rainfall (%). Values in the maps are the average over the period from 13 June to 7 October 2008.
5.4. Sensitivity Analysis of the Model
There are five parameters used in the RS-Gash model including VAI, FVC, specific vegetation storage SV (used to calculate the vegetation storage capacity Sveg), mean evaporation rate EV, and mean rainfall rate R. The values of these parameters can be determined by remote sensing observations directly or indirectly, by field measurements, or by taking as empirical values when applying the model at regional scale. Errors of these parameters are unavoidable and may have different effects on the estimated interception loss at different gross rainfall levels. It is necessary to evaluate the sensitivity of the model to these parameters. One dataset used in our study on Picea crassifolia forest were taken as the baseline for the sensitivity study. The sensitivity analysis is done by adding ±20% and ±10% of errors to each one of the parameters while keeping the others constant. The difference between the Poisson distribution and the uniform distribution was also analyzed.
The baseline values of the parameters are as follows:
VAI = 2.5;
FVC = 0.6;
SV = 0.35;
EV = 0.2
mm∙h−1;
R = 1.0
mm∙h−1. The sensitivity of the estimated interception loss on each parameter is shown in
Figure 6. The result of interception loss
vs. gross rainfall between the Poisson distribution and the uniform distribution is shown in
Figure 7.
Figure 6.
Sensitivity of the estimated interception loss on: (A) SV; (B) VAI; (C) FVC; (D) R; (E) EV.
Figure 6.
Sensitivity of the estimated interception loss on: (A) SV; (B) VAI; (C) FVC; (D) R; (E) EV.
There are several findings from the results of the sensitivity simulations:
For VAI and
SV, the error of interception loss was a constant and equal to Δ
VAI∙
SV or Δ
sv∙
VAI when
PG ˃
;
For FVC、
EV、
R, the error of interception loss is linear with gross rainfall, and the coefficient is Δ
FVC,Δ
EV and Δ
1/R, when
PG ˃
, respectively;
Interception loss using the Poisson distribution is smaller than using uniform distribution, and there are maximum errors near to SV∙VAI/FVC. However the error can be neglected for relatively larger or smaller rainfall. For forest, SV∙VAI/FVC can range from 1 mm to 10 mm, so the heterogeneity of the pixel cannot be neglected, especially for the arid region where the rainfall is dominated by small rainfall.
In general, estimated interception loss is robust and not overtly sensitive to any of the inputs. The sensitivity analysis has shown that changes in any of the input terms (
x) in the model yield a conservative change (
dI/
dx < 1) in predicted
I (
Figure 7).
Figure 7.
The results of interception loss vs. gross rainfall between Poisson distribution and uniform distribution, where SV∙VAI/FVC is equal to 3 (mm), 6 (mm) and 10 (mm) respectively.
Figure 7.
The results of interception loss vs. gross rainfall between Poisson distribution and uniform distribution, where SV∙VAI/FVC is equal to 3 (mm), 6 (mm) and 10 (mm) respectively.
6. Conclusions
A RS-Gash model was developed based on the Gash (1995) model to estimate the interception loss of heterogeneous forest using remote sensing observations of forest canopy structure parameters and gross rainfall. Compared with the original Gash (1995) model, the improvements in the RS-Gash model are: (1) The forest canopy and trunk are taken as one unity without distinguishing them as in the original Gash (1995) model. In turn, “canopy” related parameters are replaced by “vegetation” related parameters, e.g., “canopy storage capacity” is replaced by “vegetation storage capacity”; (2) Utilizing a simple linear relationship to represent the vegetation storage capacity, the latter can be easily obtained from remote sensing observations; (3) Using FVC to replace c allows the evaporation from saturated trunk to be included in the model; (4) Contrary to the assumption of uniform forest, the heterogeneity of vegetation distribution both at pixel and at sub-pixel scales is considered in the RS-Gash model.
For validation, of a single rainfall, the modelled interception loss is in good agreement with the field measurements at the Dayekou (R2 = 0.91 and RMSE = 0.34 mm∙d−1) and Pailugou (R2 = 0.82 and RMSE = 0.6 mm∙d−1) experimental sites, showing that the model can be used on a daily scale, assuming that there is only one storm per rainday. For the measured period, the modelled interception losses were underestimated by about 14.2% and 13.6% measured at Dayekou and Pailugou experimental sites, respectively, because of the uncertainty of the parameters coming from the satellite at a different scale. For the study area which is in an inland river basin, there are many small rainfall events, so the effect of vegetation heterogeneity should not be neglected.
In this study, we assumed that the EV is a constant. This is reasonable for this study due to the relative small area of the study area. Whereas, when the RS-Gash model is applied on a much larger scale, such as the basin scale, the spatial variation of EV should be considered. Until now, there are several regional near-surface meteorological forcing data products produced, based on model simulation, and this issue may be solved but there remains other challenging work.