Next Article in Journal
Quantifying Reaeration Rates in Alpine Streams Using Deliberate Gas Tracer Experiments
Previous Article in Journal
A Review of Sustainable Urban Drainage Systems Considering the Climate Change and Urbanization Impacts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Modified Gash Model for Estimating Rainfall Interception Loss of Forest Using Remote Sensing Observations at Regional Scale

1
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Water 2014, 6(4), 993-1012; https://doi.org/10.3390/w6040993
Submission received: 11 February 2014 / Revised: 7 April 2014 / Accepted: 10 April 2014 / Published: 22 April 2014

Abstract

:
Rainfall interception loss of forest is an important component of water balance in a forested ecosystem. The Gash analytical model has been widely used to estimate the forest interception loss at field scale. In this study, we proposed a simple model to estimate rainfall interception loss of heterogeneous forest at regional scale with several reasonable assumptions using remote sensing observations. The model is a modified Gash analytical model using easily measured parameters of forest structure from satellite data and extends the original Gash model from point-scale to the regional scale. Preliminary results, using remote sensing data from Moderate Resolution Imaging Spectroradiometer (MODIS) products, field measured rainfall data, and meteorological data of the Automatic Weather Station (AWS) over a picea crassifolia forest in the upper reaches of the Heihe River Basin in northwestern China, showed reasonable accuracy in estimating rainfall interception loss at both the Dayekou experimental site (R2 = 0.91, RMSE = 0.34 mmd −1) and the Pailugou experimental site (R2 = 0.82, RMSE = 0.6 mmd −1), compared with ground measurements based on per unit area of forest. The interception loss map of the study area was shown to be strongly heterogeneous. The modified model has robust physics and is insensitive to the input parameters, according to the sensitivity analysis using numerical simulations. The modified model appears to be stable and easy to be applied for operational estimation of interception loss over large areas.

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 ( Water 06 00993 i010) 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 Water 06 00993 i010 so that forest canopy is not saturated, and interception loss is equal to cPG, where c is canopy coverage; (ii) the gross rainfall PG is larger than Water 06 00993 i010 but not larger than St/pt so that forest canopy is saturated but the trunk is not, and the interception loss is equal to Water 06 00993 i011 , 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 Water 06 00993 i012, 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 Water 06 00993 i010 is given by:
Water 06 00993 i001
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:
Water 06 00993 i002
The interception loss of trunk (It) can be calculated by:
Water 06 00993 i003
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.

3. RS-Gash Model for Regional Scale Estimation

3.1. The RS-Gash Model Assumptions

As stated above, there are no methods to directly calculate St and pt using other easily measured forest structure parameters, while existing methods to calculate Sc using LAI can only partly represent Sc in the application of the original Gash model. Treating St and pt as constant at regional scale is a possible way with the uniformity forest, however in a heterogeneous forest it is not, and the latter case is more common in our world. Besides, when calculating the interception loss, the heterogeneity within a pixel is still not included in the model. Here, the RS-Gash model was proposed to simplify the Gash (1995) model, so that the interception loss of heterogeneous forest can be obtained using remote sensing observations at regional scale. The major improvement of the RS-Gash model is that the interception loss is calculated based on the whole vegetation (forest), not on canopy and trunk separately, and meanwhile the heterogeneity within a pixel is also considered. The RS-Gash model is based on the following hypotheses:
  • The leaves, branches and trunks of forest are treated as one unit so there is no distinction between canopy and trunk as in the original Gash (1995) model. A term of vegetation storage capacity per unit area of ground (Sveg) is introduced to replace the S and St in the original Gash (1995) model. “Vegetation” here indicates all elements of forest including leaves (green and dry), branches, and trunk. The vegetation storage capacity Sveg is linearly related to Vegetation Area Index (VAI), the latter includes green leaves, dry leaves, branches, and trunk areas;
  • Accordingly, a term of mean evaporation rate per unit vegetation coverage area from saturated vegetation surfaces (EV) is introduced assuming the saturated canopy and saturated trunk have the same evaporation rate. The hypothesis is inherited that EV/R is equal for all storms as used in the original Gash (1995) model;
  • In the RS-Gash model, for one pixel of satellite image, the VAI is divided into several sub-pixels by Poisson distribution. Each sub-pixel is treated as having homogeneous distribution of VAI when calculating the interception loss of it. The interception loss of the entire pixel is the integration of all sub-pixels by the Poisson distribution probability. The interaction effect among different sub-pixels in one pixel is neglected.
Based on these assumptions, the RS-Gash model can be described as follows. For a single rainfall, according to assumption (3), the interception loss of one pixel is described as:
Water 06 00993 i004
where Ii is the interception loss of the i–th sub-pixel of the pixel, and Pbi is the corresponding probability which will be given in Section 3.4.
For the i–th sub-pixel, the amount of gross rainfall necessary to saturate the vegetation ( Water 06 00993 i013) now can be expressed as:
Water 06 00993 i005
All the variables and parameters in Equation (5) have been defined before.
For low rainfall (PG ˂ Water 06 00993 i013), the interception loss of one sub-pixel (Ii) is given by:
Ii = FVC∙PG
For high rainfall (PG ˃ Water 06 00993 i013), Ii is given by:
Water 06 00993 i006
The ratio of EV over R should be obtained during a single rainfall. However, considering the measurement error especially for a short duration of rainfall, this ratio may not be stable or even not possible to obtain. So, the ratio for a long-term (the whole experimental period in this paper) is a feasible replacement, based on the assumption (2) and the sensitivity analysis of the model in Section 5.4. In the RS-Gash model, all of the parameters, that is mean evaporation rate EV, mean rainfall rate R, vegetation storage capacity Sveg, VAI and FVC, could be obtained from remote sensing observations, using the existing methods or the methods described below.

3.2. Vegetation Storage Capacity

Gash-like models divide the interception loss of forest into two parts: interception loss of canopy and trunk [4,5,8,9,23]. Some researchers thought that the canopy storage capacity was more important because the trunk’s surface area is relatively small and the rain can reach the trunk only when the rainfall is large enough [8]. Gash considered that the canopy storage capacity is linear with canopy coverage (c) [8]. Many researchers have suggested that the canopy storage capacity has a linear relationship with LAI, e.g., Bastiaanssen et al. [13] used one empirical coefficient to represent the canopy storage capacity per unit LAI, Pitman [23] obtained the relationship between canopy storage capacity and LAI by regression using ground measurements, Liu [9] and van Dijk and Bruijnzeel [10] considered linear relationship reasonably exiting between S and LAI. However, this still cannot explain the interception loss of forest in a leafless period, as Llorens and Gallart [24] noted in that the water retained in branches and stems plays a key role in rainfall interception. It is reasonable that all parts of the forest have the ability to hold rainfall and not only the leaves. Assuming a linear relationship existing between Sveg and VAI, the vegetation storage capacity per unit area of ground Sveg that includes canopy storage capacity and trunk storage capacity can be expressed as:
Sveg = SV ∙ VAI
where SV is specific vegetation storage, defined as the depth of water that can be retained by vegetation per unit vegetation area, which can be determined experimentally for a particular species using laboratory methods described for instance by Aston [25] and Liu [9] or the regression method on throughfall data [26]. This implies that the leaf and other parts of forest have the same SV.
At regional scale, the VAI can be calculated from the Moderate Resolution Imaging Spectroradiometer (MODIS) products of LAI by the methods proposed by Zeng et al. [21]:
VAI = LAIg + Is
where LAIg is the green leaf-area index, which can be provided by the MODIS LAI product; Ls is the area index of branches, trunk and dead leaves and can be obtained by:
Water 06 00993 i007
where n denotes the nth component in the Ls time series; Ls,min is the prescribed minimum value of Ls; (1 − α) is the removal rate of dead leaves between two products [21].

3.3. Mean Evaporation Rate from Saturated Vegetation

Evaporation of intercepted rainfall from forest consists of two parts: one is the evaporation from the tree canopy, and the other is from the trunk. Evaporation rate is often calculated by the Penman-Monteith equation by setting the surface resistance to zero [27]:
λE = (ΔRn + ρcpD/ra)(Δ + γ)-1
ra = In((z − d)/z0)2/(κ2u)
where E (kgm−2s−1) is the evaporation rate; λ (Jkg−1) is the latent heat of vaporization; Rn (Wm−2) is the net radiation; cp (Jkg−1K−1) is the specific heat capacity of air at constant pressure; D (Pa) is the vapor pressure deficit; γ (PaK−1) is the psychrometric constant; Δ (PaK−1) is the slope of the saturation vapor pressure curve at air temperature; ρ (kgm−3) is the density of air; ra (sm−1) is the aerodynamic resistance; z (m) is the height of wind speed u (sm−1) measurements; the roughness length parameters are d = 0.75 h and z0 = 0.1 h, h (m) is the vegetation height; and κ (0.41) is the von-Karman constant.
It is reasonable to treat the evaporation rate of the saturated canopy and the saturated trunk equally assuming all parts of forest have similar net radiation (per unit area), vapor pressure deficit and other atmospheric forcing factors. As Gash [8] pointed out, the mean evaporation rate per unit ground is linearly correlated with the canopy cover (c) for spare forest. As mentioned earlier, a new term—the mean evaporation rate from saturated vegetation, EV—is introduced to represent the total evaporation from forest canopy and trunk per unit vegetation cover. Consequently, the FVC is considered more adequate than canopy cover (c) to calculate the mean evaporation per unit ground area (E). Although FVC and c show almost no difference for the leafed period, a large difference for the leafless period may exist between the two. Hence, the mean evaporation rate from saturated canopy per unit ground area E is expressed by:
E = EVFVC
Note that the ratio EV/R is assumed to be constant when FVC is varying.
At regional scale, the FVC can be calculated from the satellite observed parameter Normalized Difference Vegetation Index (NDVI) by the methods used in Bastiaanssen et al. [13]:
Water 06 00993 i008
where NDVIfv and NDVIbs are the threshold values of full vegetation cover and bare soil, respectively.

3.4. Heterogeneous Vegetation at Sub-Pixel Scale

The heterogeneity of vegetation at pixel scale can be described by VAI and FVC. Most of the interception loss models did not take into account the sub-pixel heterogeneity, for they are based on uniformity forest hypothesis that assumes that the forest reaches saturation at almost the same time. Unfortunately, forest often is characterized by very large spatial heterogeneity. Larger VAI implies taking longer time to reach saturation. There are mainly three forms describing the spatial distribution of VAI according to the distribution of canopy elements: regular distribution, random distribution and irregular distribution [28]. The random distribution is the normal form and widely used. The Poisson distribution function is applied to take account of the random distribution to describe effects of heterogeneity of forest density of coverage on interception loss.
The Poisson distribution scheme of a pixel is presented in Figure 1.
If a pixel with an average value of VAI is divided into n sub-pixels according to the Poisson distribution, the pixel value of VAI can be described as the composite of the VAI values of the Poisson-distribution-divided sub-pixels by the equation:
Water 06 00993 i009
where Pbi is the proportion of sub-pixels where VAI = VAIi have the same FVC as the pixel. The interior of sub-pixels in a pixel can be treated as homogeneous when calculating the interception loss. Consequently, the sub-pixels with different VAI values will have a different time and gross rainfall requirement to reach saturation.
Figure 1. Schematic representation of spatial heterogeneity of Vegetation Area Index (VAI) in a pixel.
Figure 1. Schematic representation of spatial heterogeneity of Vegetation Area Index (VAI) in a pixel.
Water 06 00993 g001

4. Study Area and Data

4.1. Study Area

The study area is a part of the water conservation forest in the Qilian Mountain, in the upper reaches of the Heihe River Basin in northwestern China, which lies between latitude 38°–39° N and longitude 99°–101.5° E, with elevation ranging between 1450 m and 4900 m (Figure 2). The field experiment is a part of Watershed Allied Telemetry Experimental Research (WATER) [29]. This part of the basin exhibits complex hydrological processes due to variability in topography, rainfall and land cover. The annual precipitation varies approximately between 400 mm and 600 mm, and the mean annual air temperature is about 0.7 °C. Picea crassifolia is the major species covering about 76.8% of the forest in this region [29].
Figure 2. Map of the Heihe River Basin (A); study area (B); and Dayekou and Pailugou experimental sites (the two small false color images) (C). The map is from MODIS NDVI on 27 June 2008. The false color image is from Hyperion on 15 June 2008 with resolution of 30 m, within a 0.01° × 0.01° grid.
Figure 2. Map of the Heihe River Basin (A); study area (B); and Dayekou and Pailugou experimental sites (the two small false color images) (C). The map is from MODIS NDVI on 27 June 2008. The false color image is from Hyperion on 15 June 2008 with resolution of 30 m, within a 0.01° × 0.01° grid.
Water 06 00993 g002
The interception loss was measured at two forest experimental sites, shown in Figure 2, from 13 June to 7 October in 2008. One site is named Dayekou forest hydrology experimental site with one experimental plot, and the other experimental site is named Pailugou forest hydrology experimental site with three experimental plots. The two sites had significantly different structures, representing the interception loss of different forest, though they are not very far from each other, and the two experiments were carried out separately. The main features of the experimental plots are shown in Table 1.
Table 1. Main features of the Dayekou and Pailugou experimental sites in the upper reaches of the Heihe River Basin in 2008.
Table 1. Main features of the Dayekou and Pailugou experimental sites in the upper reaches of the Heihe River Basin in 2008.
FeaturesDayekouPailugou
Location(100.258° E, 38.507° N)(100.295° E, 38.543° N)
Elevation (m)31283250
Tree speciespicea crassifoliapicea crassifolia
Tree average height (m)11.58.4
Density (number of trees ha-1)12321880
Average DBH (diameter at breast height) (cm)16.313.2
canopy cover*0.5650.75
Canopy storage*0.551.51
UnderstoreyWith some BryophytaWith some Bryophyta
Plot dimensions (m)25 × 2520 × 15
10 × 15
20 × 25
Note: *: determined using regression method.

4.2. Remote Sensing Data

The VAI and FVC can be retrieved from remote sensing observations. In this paper, we used two MODIS products: MOD15A2 LAI products and MOD13A2 NDVI products in 2008, provided by the MODIS land team of the National Aeronautics and Space Administration (NASA) Earth Observing System (EOS) project. Spatial resolution of all the MODIS data is 0.01° and temporal resolution is eight days.
Land cover data is provided by the Heihe experimental data center that combined the land use data covering many years with the MOD12Q1 product using the International Geosphere Biosphere Programme (IGBP) classification [30].
The fraction of forest (FF) in pixel where the experimental sites were located was estimated from the classification of the Hyperion image, with a resolution of 30 m, on 15 June 2008. The modelled interception loss is based on per unit area of ground, but the measurements are based on per unit area of forest, so FF can be used to convert the former to the latter when validating modelled interception loss.
There are two types of rainfall data used in this study. One is the field measured rainfall data, which was used to validate the RS-Gash model. Another is the product of Tropical Rainfall Measuring Mission (TRMM) TRMM 3B42, with temporal resolution of 3 h and spatial resolution of 0.25°.The linear interpolation method was used to obtain a finer spatial resolution at 0.01°, which was then used to map the interception loss of the study area. TRMM is a joint mission between NASA and the Japan Aerospace Exploration Agency (JAXA) designed to monitor and study tropical rainfall (more information see [31]).

4.3. Meteorological Forcing Data

In the Dayekou forest hydrology experimental site, one Automatic Weather Station (AWS) was built which can measure net radiation, wind speed, air temperature, air pressure and air humidity. The detailed information of the measured variables used in this study is listed in Table 2.
Table 2. The detailed information of the measured variables used in this study from the Automatic Weather Station (AWS) at the Dayakou forest hydrology experimental site.
Table 2. The detailed information of the measured variables used in this study from the Automatic Weather Station (AWS) at the Dayakou forest hydrology experimental site.
VariablesSensorManufacturerObservation accuracyObserved altitude
Tair 24 mHMP45CVaisala±0.2 °C23.75 m
(Helsinki, Finland)
Humidity 24 mHMP45CVaisala±2%23.75 m
(Helsinki, Finland)
Wind speed 24 m034BMetOne±0.11 m/s24.00 m
(Grants Pass, OR, USA)
PressureCS105Campbell±0.5 mb0.50 m
(Logan, UT, USA)
Short-wave radiation (upward and downward)CM3Campbell±10%19.75 m
(Logan, UT, USA)
Long-wave Radiation (upward and downward)CG3Campbell±10%19.75 m
(Logan, UT, USA)

4.4. Validation Data

4.4.1. Dayekou Experimental Site

Gross rainfall was measured using two manual collecting rain-gauges with a 20 cm diameter set up on the open ground beside the tree canopy.
Throughfall was measured with 10 troughs with dimensions 100 cm long, 20 cm wide and 25 cm high. Nine of them were placed in a 3 × 3 grid, 12.5 m apart from each other. The remaining one was placed on the ground in between the 3 × 3 array. The rain-gauges and troughs were cleaned up after the rainfall (sometimes on the next day for night rainfall). For single rainfall, measurements of throughfall larger than gross rainfall were discarded. The average throughfall was calculated by the method of Herbst et al. [32], which adjusted the remainder according to their relative contribution to the average when all troughs were working. It should be noted that this method has a disadvantage: there are some errors when the missing data are the smallest or largest so that the throughfall might be underestimated if the amount of abandoned data is larger. Consequently, the measured interception loss will be slightly overestimated.
Stemflow was collected from five trees using coiled pipes randomly. The stand-scale stemflow volume (in mm) was calculated by taking into account the tree density of the plot.
There were a total of 25 rainfall data in the Dayekou site from 13 June to 7 October 2008, assuming that there was only one storm per rainday as Gash did [6]. The total gross rainfall was 162.8 mm with a total rainfall time of 220.6 h.

4.4.2. Pailugou Experimental Site

Gross rainfall was measured using three siphon rainfall recorders set up on the open ground beside the tree canopy. Throughfall in each plot was measured with six troughs with dimensions 100 cm long, 20 cm wide and 25 cm high. All of them were placed in a fixed grid. For single rainfall, outliers were removed using the same method as for the Dayekou site. Since three siphon rainfall recorders were used at the Pailugou site, additional data screening policy is applied: for one rainfall event, if the difference of gross rainfall measurements between any two recorders among the three is larger than 1mm, the measurements are considered ineffective and screened out.
Stemflow was collected from nine trees using coiled pipes randomly. The stand-scale stemflow volume (in mm) was calculated by taking into account the average tree density of the three plots.
There were a total of 23 rainfall data in the period from 13 June to 7 October 2008, assuming that there was only one storm per rainday as Gash did [6]. The total gross rainfall amount was 153.5 mm and the total rainfall time of 147 h.

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 mmh−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.2mm∙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 mmh−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 mmh−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:
Icorrect = Imodelled/FF
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 mmh−1 (Figure 3a), at Pailugou experimental site with a R2 of 0.82 and RMSE of 0.6 mmh−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.
Water 06 00993 g003
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.
Water 06 00993 g004
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.
SiteGross rainfall (Pg) (mm)ModelledMeasuredRelative error (% measured)
mm% Pgmm% Pg
Dayekou162.833.820.839.424.214.2
Pailugou153.549.732.457.537.513.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.
VariablesMeanStandard deviation
FVC0.470.16
VAI1.940.64
Interception loss (mm)61.126.1
Interception loss (%)20.056.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.
Water 06 00993 g005

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.2mm∙h−1; R = 1.0mm∙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.
Water 06 00993 g006
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 ΔVAISV or ΔsvVAI when PG ˃ Water 06 00993 i013;
  • For FVC、EVR, the error of interception loss is linear with gross rainfall, and the coefficient is ΔFVCEV and Δ1/R, when PG ˃ Water 06 00993 i013, respectively;
  • Interception loss using the Poisson distribution is smaller than using uniform distribution, and there are maximum errors near to SVVAI/FVC. However the error can be neglected for relatively larger or smaller rainfall. For forest, SVVAI/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 SVVAI/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 SVVAI/FVC is equal to 3 (mm), 6 (mm) and 10 (mm) respectively.
Water 06 00993 g007

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.

Acknowledgements

This work was jointly supported by the “China Academy of Sciences, State Bureau of foreign experts and international innovation Partnership Program” project (Grant No. KZZD-EW-TZ-09), the National High Technology Research and Development Program of China (Grant No. 2012AA12A304). The authors thank the anonymous reviewers for their valuable comments. The authors would like to sincerely thank the WATER program for providing the in situ measurements.

Author Contributions

Yaokui Cui wrote the paper; Li Jia supervised the study, reviewed and edited the manuscript. All authors read and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Scatena, F.N. Watershed scale rainfall interception on 2 forested watersheds in the luquillo mountains of puerto-rico. J. Hydrol. 1990, 113, 89–102. [Google Scholar] [CrossRef]
  2. Asdak, C.; Jarvis, P.G.; van Gardingen, P.; Fraser, A. Rainfall interception loss in unlogged and logged forest areas of central kalimantan, Indonesia. J. Hydrol. 1998, 206, 237–244. [Google Scholar] [CrossRef]
  3. Levia, D.F.; Frost, E.E. Variability of throughfall volume and solute inputs in wooded ecosystems. Prog. Phys. Geog. 2006, 30, 605–632. [Google Scholar] [CrossRef]
  4. Calder, I.R. Water use by forests, limits and controls. Tree Physiol. 1998, 18, 625–631. [Google Scholar] [CrossRef]
  5. Jackson, I.J. Relationships between rainfall parameters and interception by tropical forest. J. Hydrol. 1975, 24, 215–238. [Google Scholar] [CrossRef]
  6. Gash, J.H. An analytical model of rainfall interception by forests. Q. J. Roy. Meteor. Soc. 1979, 105, 43–55. [Google Scholar] [CrossRef]
  7. Rutter, A.J.; Robins, P.C.; Morton, A.J.; Kershaw, K.A. Predictive model of rainfall interception in forests.1. Derivation of model from observations in a plantation of corsican pine. Agr. Meteorol. 1972, 9, 367–384. [Google Scholar]
  8. Gash, J.H.C.; Lloyd, C.R.; Lachaud, G. Estimating sparse forest rainfall interception with an analytical model. J. Hydrol. 1995, 170, 79–86. [Google Scholar] [CrossRef]
  9. Liu, S.G. Estimation of rainfall storage capacity in the canopies of cypress wetlands and slash pine uplands in North-Central Florida. J. Hydrol. 1998, 207, 32–41. [Google Scholar] [CrossRef]
  10. Van Dijk, A.I.J.M.; Bruijnzeel, L.A. Modelling rainfall interception by vegetation of variable density using an adapted analytical model. Part 1. Model description. J. Hydrol. 2001, 247, 230–238. [Google Scholar] [CrossRef]
  11. Murakami, S. A proposal for a new forest canopy interception mechanism: Splash droplet evaporation. J. Hydrol. 2006, 319, 72–82. [Google Scholar] [CrossRef]
  12. Muzylo, A.; Llorens, P.; Valente, F.; Keizer, J.J.; Domingo, F.; Gash, J.H.C. A review of rainfall interception modelling. J. Hydrol. 2009, 370, 191–206. [Google Scholar] [CrossRef]
  13. Bastiaanssen, W.G.M.; Cheema, M.J.M.; Immerzeel, W.W.; Miltenburg, I.J.; Pelgrum, H. Surface energy balance and actual evapotranspiration of the transboundary indus basin estimated from satellite measurements and the etlook model. Water Resour. Res. 2012, 48, W11512. [Google Scholar]
  14. Von Hoyningen, H.J. Die interception des niederschlags in landwirtschaftlichen bestanden. Schriftenr. DVWK 1983, 57, 1. [Google Scholar]
  15. Mu, Q.Z.; Zhao, M.S.; Running, S.W. Improvements to a modis global terrestrial evapotranspiration algorithm. Remote Sens. Environ. 2011, 115, 1781–1800. [Google Scholar] [CrossRef]
  16. Hutjes, R.W.A.; Wierda, A.; Veen, A.W.L. Rainfall interception in the tai forest, ivory-coast—Application of 2 simulation-models to a humid tropical system. J. Hydrol. 1990, 114, 259–275. [Google Scholar] [CrossRef]
  17. Valente, F.; David, J.S.; Gash, J.H.C. Modelling interception loss for two sparse eucalypt and pine forests in central portugal using reformulated rutter and gash analytical models. J. Hydrol. 1997, 190, 141–162. [Google Scholar] [CrossRef]
  18. Schellekens, J.; Scatena, F.N.; Bruijnzeel, L.A.; Wickel, A.J. Modelling rainfall interception by a lowland tropical rain forest in northeastern puerto rico. J. Hydrol. 1999, 225, 168–184. [Google Scholar] [CrossRef]
  19. Motahari, M.; Attarod, P.; Pypker, T.G.; Etemad, V.; Shirvany, A. Rainfall interception in a pinus eldarica plantation in a semi-arid climate zone: An application of the gash model. J. Agr. Sci. Tech. Iran. 2013, 15, 981–994. [Google Scholar]
  20. Van Dijk, A.I.J.M.; Bruijnzeel, L.A. Modelling rainfall interception by vegetation of variable density using an adapted analytical model. Part 2. Model validation for a tropical upland mixed cropping system. J. Hydrol. 2001, 247, 239–262. [Google Scholar] [CrossRef]
  21. Zeng, X.B.; Shaikh, M.; Dai, Y.J.; Dickinson, R.E.; Myneni, R. Coupling of the common land model to the ncar community climate model. J. Clim. 2002, 15, 1832–1854. [Google Scholar] [CrossRef]
  22. Miralles, D.G.; Gash, J.H.; Holmes, T.R.H.; de Jeu, R.A.M.; Dolman, A.J. Global canopy interception from satellite observations. J. Geophys. Res. Atmos. 2010, 115. [Google Scholar] [CrossRef]
  23. Pitman, J.I. Rainfall interception by bracken in open habitats—Relations between leaf-area, canopy storage and drainage rate. J. Hydrol. 1989, 105, 317–334. [Google Scholar] [CrossRef]
  24. Llorens, P.; Gallart, F. A simplified method for forest water storage capacity measurement. J. Hydrol. 2000, 240, 131–144. [Google Scholar] [CrossRef]
  25. Aston, A.R. Rainfall interception by 8 small trees. J. Hydrol. 1979, 42, 383–396. [Google Scholar] [CrossRef]
  26. Lankreijer, H.; Lundberg, A.; Grelle, A.; Lindroth, A.; Seibert, J. Evaporation and storage of intercepted rain analysed by comparing two models applied to a boreal forest. Agr. Forest Meteorol. 1999, 98–99, 595–604. [Google Scholar]
  27. Monteith, J.L. Evaporation and environment. Symp. Soc. Exp. Biol. 1965, 19, 205–223. [Google Scholar]
  28. Myneni, R.B.; Ross, J.; Asrar, G. A review on the theory of photon transport in leaf canopies. Agr. Forest Meteorol. 1989, 45, 1–153. [Google Scholar] [CrossRef]
  29. Li, X.; Li, X.W.; Li, Z.Y.; Ma, M.G.; Wang, J.; Xiao, Q.; Liu, Q.; Che, T.; Chen, E.X.; Yan, G.J.; et al. Watershed allied telemetry experimental research. J. Geophys. Res. Atmos. 2009, 114. [Google Scholar] [CrossRef]
  30. Ran, Y.H.; Li, X.; Lu, L.; Li, Z.Y. Large-scale land cover mapping with the integration of multi-source information based on the dempster-shafer theory. Int. J. Geogr. Inf. Sci. 2012, 26, 169–191. [Google Scholar] [CrossRef]
  31. Tropical Rainfall Measuring Mission (TRMM). Available online: http://trmm.gsfc.nasa.gov/ (accessed on 16 April 2014).
  32. Herbst, M.; Rosier, P.T.W.; McNeil, D.D.; Harding, R.J.; Gowing, D.J. Seasonal variability of interception evaporation from the canopy of a mixed deciduous forest. Agr. Forest Meteorol. 2008, 148, 1655–1667. [Google Scholar]
  33. Jia, L.; Shang, H.; Hu, G.; Menenti, M. Phenological response of vegetation to upstream river flow in the heihe rive basin by time series analysis of modis data. Hydrol. Earth Syst. Sci. 2011, 15, 1047–1064. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Cui, Y.; Jia, L. A Modified Gash Model for Estimating Rainfall Interception Loss of Forest Using Remote Sensing Observations at Regional Scale. Water 2014, 6, 993-1012. https://doi.org/10.3390/w6040993

AMA Style

Cui Y, Jia L. A Modified Gash Model for Estimating Rainfall Interception Loss of Forest Using Remote Sensing Observations at Regional Scale. Water. 2014; 6(4):993-1012. https://doi.org/10.3390/w6040993

Chicago/Turabian Style

Cui, Yaokui, and Li Jia. 2014. "A Modified Gash Model for Estimating Rainfall Interception Loss of Forest Using Remote Sensing Observations at Regional Scale" Water 6, no. 4: 993-1012. https://doi.org/10.3390/w6040993

APA Style

Cui, Y., & Jia, L. (2014). A Modified Gash Model for Estimating Rainfall Interception Loss of Forest Using Remote Sensing Observations at Regional Scale. Water, 6(4), 993-1012. https://doi.org/10.3390/w6040993

Article Metrics

Back to TopTop