1. Introduction
The amount of precipitation intercepted and subsequently evaporated from a forest canopy is considerable, accounting for 29% of the incident precipitation at Plynlimon Forest, United Kingdom (Calder 1990). Canopy interception is a critical component of a forest’s water budget, affecting the amount of water available to the understory and soil. The structural characteristics of the forest canopy reduce and spatially redistribute the overall water input to the underlying area (Lloyd and Marques 1988; Herwitz and Slye 1995; Loescher et al. 2002; Abrahams et al. 2003). Canopy structure has been documented to be a key control on the subcanopy water flux, affecting mineral and trace-gas fluxes, surface runoff, and erosion (Whelan and Anderson 1996; Levia and Frost 2003).
The amount of precipitation lost to canopy interception and evaporation is important in determining the water balance of a forested area (Calder 1976, 1977; Stewart 1977; Dolman 1987). Many forest canopies are highly coupled to the overlying atmosphere, which enhances moisture exchange; a moist canopy will readily evaporate water into the surrounding air (Rutter 1967; Stewart 1978; Harwood et al. 1999). Several recent studies have suggested that certain canopies, particularly those in tropical climates, may be decoupled to the atmospheric boundary layer (e.g., Derbyshire 1999; Kruijt et al. 2000; Delage et al. 2002). The precise mechanism behind this decoupling is still unclear, as is its applicability to canopies other than those in the above studies. Since interception loss varies with species (Lankreijer et al. 1993; Liu 1998; Wei et al. 2005), accurate measurements can aid predictions of the consequences of land-cover changes. Tree growth, streamflow, groundwater recharge, and water stress are closely related to evaporation, which is in turn strongly correlated to the canopy storage capacity (Halldin et al. 1984).
While one can use throughfall and stemflow measurements to calculate interception loss, this requires frequent data collection and offers few opportunities to extrapolate results to other species or locations. To create more applicable formulas, researchers have constructed dynamic, analytic, and numerical models of interception loss. Rutter et al. (1971) was among the first to derive an interception-loss model, which was based on observations of a Corsican pine (Pinus nigra) stand. The Rutter et al. (1971) model was later generalized and provided accurate interception estimates for two coniferus (Douglas fir and Norway spruce) and three hardwood (one hornbeam and two oak) stands (Rutter et al. 1975). Domingo et al. (1998) found that the Rutter model’s estimates of interception agreed with measured values for several semiarid species. While a significant advance, the Rutter model requires a substantial number of input fields, making it difficult to use with routine meteorological data.
Several models developed more recently have overcome the limitations of the Rutter model. This study compares the results from three such canopy interception models against measured interception for a leafless deciduous stand in the mid-Atlantic region of the United States. The three models used are the Mulder (1985) model, the Gash et al. (1995) model, and the WiMo model of Hörmann et al. (1996). The Gash et al. (1995) model is an updated version of the Gash (1979) model that has been reformulated for sparse canopies. By comparing measured and modeled interception, this study will determine the degree to which each model may be generalized to leafless deciduous stands in a temperate winter climate.
Many studies have compared the performance of the Mulder (1985) model and the older Gash (1979) model. Lankreijer et al. (1993) ascertained that both models overestimated interception for a deciduous forest in the Netherlands and a coniferous forest in France. Dolman (1987) compared these models for an oak forest in the Netherlands and found that the Mulder model overestimated interception relative to both the Gash (1979) model and measured values. Hutjes et al. (1990) concluded that both models accurately predicted interception for an evergreen forest in the Ivory Coast. Liu (2001) compared the Liu (1997) model against the Gash (1979) and Rutter models. The Liu model performed favorably in the more than 20 canopies considered, yet required fewer input variables and contained simpler formulas than either the Gash or Rutter model. The WiMo model has received relatively little attention in the literature and has seldom been compared with other interception models. Few studies have compared interception models against measurements taken in a deciduous stand in the leafless season. The canopy intercepts a far lower fraction of the incident precipitation during the leafless season, yet interception in this season is still a critical component of the water balance. Therefore, whether existing interception models can provide accurate predictions during the leafless season is of importance.
2. Study area
The experimental forest plot (800 m2) was located in the Fair Hill Natural Resource Management Area in Fair Hill, Maryland (39°42′N, 75°50′W) at an elevation of 72 m above sea level (Fig. 1). The rectangular study plot was nearly level with a stand density of 225 trees ha−1 and a stand basal area of 36.8 m2 ha−1, indicating that there were relatively few but large trees within the study plot. The canopy trees had a mean diameter at breast height (dbh, 1.37 m above ground level) of 40.8 cm and a mean height of 27.8 m. The mean bark thickness of the canopy trees was 7.0 mm. Fagus grandifolia Ehrh. (American beech) was the most common tree species within the plot, accounting for approximately 50% of the canopy trees. Liriodendron tulipifera L. (yellow poplar) was the second most common species, constituting approximately 25% of the canopy trees. The remaining 25% of canopy trees were composed of Quercus velutina Lam. (black oak), Acer rubrum L. (red maple), and Acer saccharinum L. (silver maple).
The Fair Hill Natural Resource Management Area is close to the Atlantic coast and experiences a temperate climate. The mean annual precipitation is 1221.5 mm and precipitation does not vary greatly with season. Summer is the wettest season (324.1 mm total for June, July, and August) and winter the driest (274.1 mm for December, January, and February). Snowfall occurs almost exclusively in winter, and the annual total (350.5-mm depth on average) is often dominated by only a few storms. The mean January (July) air temperature is −0.6°C (24.6°C). The growing season averages 191 days in April through October.
3. Materials and methods
a. Throughfall and stemflow
Throughfall and stemflow measurements were manually collected approximately two hours after each of 11 rainstorms from 1 December 2004 until 30 April 2005 (Table 1). In the text, we will refer to storms by the date on which they began. The storms cover a wide spectrum—from a warm, driving rain (23 December) to a cold drizzle (7 January) to an intense spring storm (1 April)—and so this dataset may be considered representative of the rainstorms that occur in this region during the leafless season.
Thirty funnel-type throughfall gauges were positioned within the plot and one tree was fitted with a stemflow collar. Eighteen of the gauges had an orifice diameter of 17.8 cm, while the remaining 12 gauges had an orifice diameter of 12.4 cm. The gauges were randomly redistributed throughout the plot after each storm to reduce the probable standard error of the throughfall measurements (Lloyd and Marques 1988). The stemflow collar was placed on a representative F. grandifolia with a similar branching architecture to the others in the plot and a dbh of 39.8 cm, close to the mean stand dbh. The two stemflow collection bins overflowed during two storms: 11 January and 22 March. Stemflow for these storms was estimated with a regression equation comparing gross precipitation (Pg) and stemflow (r 2 = 0.977, p < 0.01) for the other nine events.
b. Meteorological conditions
A clearing about 800 m to the east of the plot housed a meteorological station equipped with a vibrating-wire rain gauge (Geonor A/S; Oslo, Norway) that measured Pg every hour. An alter-type wind shield around the rain gauge reduced collection errors from wind-driven precipitation. Temperature, wind speed and direction, relative humidity, incident radiation, and barometric pressure were measured every 5 min.
c. Wet-canopy evaporation






Finally, we converted E′ to a value of millimeters per hour to match our precipitation measurements. The storm-averaged evaporation rate (
d. Stand and canopy parameters
Values for d, z0,M, and z0,H were calculated as 0.5h (Rutter et al. 1975), 0.36(h–d) (Thom 1971), and 0.14z0,M (Garratt and Francey 1978), respectively, where h is the mean tree height. We estimated that the measurement height z was several meters above h. We assumed that the stand albedo (α) remained constant during the leafless season at the median value found by Moore et al. (1996). That study analyzed 3 yr of observations for a temperate deciduous forest and found that α varied between 0.10 and 0.20 during the leafless season.
Following leaf abscission, a Li-Cor LI-2000 Plant Canopy Analyzer measured a plant-area index (PAI) of 1.19 m2 m−2. The coefficient of free throughfall (ct) was assumed to equal 0.3 for a canopy in full leaf (Mulder 1985), and thereafter to increase linearly as the PAI decreased throughout the autumn; a hypothetical canopy with a PAI of zero should have a ct of unity. We used a linear regression relating ct to PAI to determine ct for our leafless stand.
A value for the canopy storage capacity (Sc) for a leafless deciduous stand was taken from Halldin et al. (1984). The stand in that study contained a substantial fraction of beech trees and had a mean tree height, basal area, and stand density remarkably similar to our stand. Values for the fraction of precipitation diverted to stemflow (cs) and trunk storage capacity (St) were obtained from Hörmann et al. (1996) for a nonfoliated beech forest. These values agree with those given in Kittredge (1948) for beech. Canopy parameter values obtained from previous work in other leafless forests were deemed to be reasonable based on field observations at the study plot. The effects of our estimations and assumptions will be discussed in section 6, where a sensitivity analysis will be conducted for each model.
4. Model descriptions
a. Common terms






b. Mulder (1985) model
Mulder (1985) developed a numerical model (hereafter the “Mulder model”) notable both for its use of standard meteorological variables to derive interception loss, as well as its ability to maintain a moist canopy between storms by calculating a running water balance. For each day, separate averages of hourly Rn and E′ are maintained for rainy and dry periods. Rainy periods are defined as those hours when the Geonor gauge recorded more than 0.1 mm of precipitation. The ability to distinguish between rainy and dry conditions is a key difference between the Mulder model and the Gash and WiMo models.
Precipitation is summed over the day and then evenly divided into small showers of equal duration (tshw). The number of showers is defined as the number of rainy periods bounded by at least 2 h of dry conditions. Showers continuing through the end of the day are considered to have stopped at midnight. Showers beginning in the first (final) two hours of a day are (are not) counted. The model distributes the showers and dry intervals evenly throughout the day. All dry intervals are initially of equal duration (tdry). The length of the last dry interval in each day is averaged with the first dry interval in the following day. Half of the averaged dry interval is placed in each day, so that every day effectively begins and ends with a shortened dry interval.


Each shower begins with the canopy in the wetting-up phase. The canopy can only become saturated if
After the shower ends, the canopy moves into the drying-out phase. The canopy water at the end of the drying-out phase is calculated by Cwd, given Cwf and the mean evaporation rate during dry intervals (
We stipulated that the canopy was completely dry at the beginning of the simulation period (1 December), such that for the very first shower Cwi was equal to zero. Sufficient wetting and drying cycles occurred before the beginning of our first storm (7 December) that the initial condition of the canopy is of little consequence.
c. Gash et al. (1995) model
Gash (1979) simplified the Rutter model into a storm-based analytical model requiring only mean rainfall and evaporation rates, the temporal pattern of rainfall, and certain forest canopy parameters. The Gash (1979) model has been applied to various canopies and climates throughout the world (e.g., Pearce and Rowe 1981; Bruijnzeel and Wiersum 1987; Dolman 1987; Lloyd et al. 1988; Loustau et al. 1992; Návar et al. 1999). While the model performed well under these circumstances, Teklehaimanot et al. (1991) found it to be less accurate in stands where the canopy was sparse. Gash et al. (1995) updated the model to better account for such forests, validating the new formulas with a thin stand of maritime pine in France. Gash et al. (1995) suggested that the updated model (hereafter “the Gash model”) be used in place of the original for all forest types. Carlyle-Moses and Price (1999) confirmed this recommendation by applying the model to hardwood stands in northern North America during the growing season.
Unlike the Mulder model, the Gash model does not maintain a running water balance between storms. Rather, the Gash model assumes that the canopy is completely dry before each storm. The Gash model does not distinguish between rainy and dry periods when calculating average evaporation, but instead takes storm-averaged values (denoted by an overbar). Stemflow is accounted for by the Gash model. The precipitation required to saturate the canopy (Psat) is a function of
d. WiMo model
Herwitz and Slye (1995) recognized the effect of wind-driven precipitation on interception, a factor that neither the Gash nor the Mulder model takes into account. Hörmann et al. (1996) presented a model that does incorporate the influence wind-driven precipitation, named “WiMo.” As described in that study, the WiMo model uses the same underlying equations as the Gash (1979) model, but adds the ability to vary Sc with the maximum sustained wind speed during a rainstorm (umax). Hörmann et al. (1996) found that WiMo proved more accurate than the Gash (1979) model for a beech forest in northern Germany.




A computer program found So for each storm by running the above bucket model for each Sc between 0.1 and 4.0 mm, in increments of 0.01 mm. For each Sc, Tf was compared against the measured value from our stand. The Sc that produced the lowest absolute difference between Tf and the measured value was judged to be the So for that storm. Ordered pairs of So and umax from each storm were used to generate a regression equation for Sc as a function of umax.
5. Results


We obtained interception values from each model for each storm in Table 1. There is little explanation for the models’ pattern of over- and underprediction for individual storms (Table 4). All three models twice overpredicted interception during high-volume events (22 March, 1 April), while underpredicting interception for two others (9 December, 11 January). During these four storms, the high precipitation rate quickly saturated the canopy, which remained saturated for the duration of the event. This reduced the total interception, as a saturated canopy can intercept at most the amount of water that evaporates from it, and during these storms
All models tended to overestimate interception during January, March, and April, while underestimating interception for the February storm and one of the December storms (9 December). The models showed the greatest overestimates when the measured percentage of Pg was low (3 January, 1 April, 6 April) and the greatest underestimates when the measured percentage was high (9 December, 9 February). These overestimates and underestimates indicate a lack of flexibility in each model, which will be discussed further in section 6.
We define PG to be the total precipitation summed over all 11 storms (i.e., ΣPg). To evaluate each model’s overall performance, we calculated 1) the total interception over the study period for models and measurements; 2) the percent error from the total measured interception; 3) the percentage of PG intercepted; and 4) the difference in this percentage between measured and modeled values. We also calculated the percentage of Pg intercepted for each storm and computed the difference between measured and modeled values. We used the difference (in millimeters) between measured and modeled values for each storm to calculate the standard deviation, skewness, kurtosis, and Nash–Sutcliffe model efficiency (Nash and Sutcliffe 1970) (Table 4). A normal Gaussian distribution was assumed to have a kurtosis of zero. The Nash–Sutcliffe efficiency does not include the 9 February storm. The models performed equally poorly for this storm and its inclusion would have resulted in negative Nash–Sutcliffe model efficiencies for all models. Removing this storm adds an equal amount to each model’s efficiency value.
The Mulder model was the most accurate in terms of the total interception, overestimating the measured value by only 3.2%. The WiMo model underestimated the total measured value by 9.5%. The Gash model lagged the other two, underestimating the total measured interception by 17.8%. All models predicted the percentage of PG intercepted to within 1.1% of PG of the measured value (5.7%). Each model’s differences formed a negatively skewed distribution, due to the 9 February outlier. It is interesting to note that while the Mulder model was the most accurate overall, it had the highest standard deviation, the lowest skewness, and a low kurtosis for its differences for individual storms. The Mulder model also had the lowest Nash–Sutcliffe model efficiency. These statistics imply that the Mulder model’s differences for individual storms were typically larger than the other two models’ differences.
When running sums of interceptions were considered, the models tracked the measured values relatively well (Fig. 2). The models oscillated between overestimating the measured running total during late December and January and underestimating it during February and March. It was only in April that the models converged toward to the measured total.
In section 3d, we noted that some of our stand and canopy parameters were adapted from previous studies. While we believe that these parameters are valid for our stand, we conducted a sensitivity analysis to examine whether the assumed parameter values had affected the models’ estimates. We conducted additional runs of each model with modified values of each parameter that was inferred or derived from a previous study. The modified values represent the bounds of a reasonable error range for the parameter in question. The parameters in the Penman–Monteith equation were multiplied and divided by 2, as were the stemflow parameters cs and St; Sc(ct) was modified by adding or subtracting 0.5 (0.05). Only one parameter was changed at a time, as nonlinear interactions between parameters were estimated to be negligible.


We note that changes to individual evaporation parameters had little or no effect on the models’ estimates of total interception. Therefore, the differences between measured and modeled values (Table 4) are not due to our assumed values for the evaporation parameters. The models were more sensitive to changes in Sc and ct. Our values for these two parameters are potentially a significant source of error. We continue to believe that the values specified for Sc and ct are appropriate for our canopy, for the reasons given in section 3d.
6. Discussion


Our stemflow measurements represent an additional source of measurement error. We measured stemflow from only one tree in the plot. This tree had a height and dbh near the stand average and was of the same species (American beech) as the majority of canopy trees. While one tree represents the smallest of samples, stemflow is not a critical component of interception in our stand. Kittredge (1948) and Hörmann et al. (1996) agree that the percentage of Pg diverted to stemflow in a beech forest is at most 10%. Kittredge (1948) found that for a storm with a Pg of 12 mm, only 7.5% became stemflow. Over all 11 storms, measured stemflow in our stand was 5% of PG. While we acknowledge that our limited sample of one tree cannot properly represent the variability in stemflow across the stand, that variability is almost certainly a negligible portion of interception in this stand.
The discrepancies between the models’ predictions of interception must be a direct result of each model’s basic assumptions. The WiMo model uses the same underlying equations as the Gash model, yet it provided a more accurate estimate of the total interception. It is highly probable that this advantage is due to the WiMo’s ability to vary Sc with umax. Whereas Hörmann et al. (1996) found that an increase in umax caused a decrease in Sc, our study found that Sc increased logarithmically with umax [Eq. (10)]. Hörmann et al. (1996) attributed their decrease in Sc to stronger winds shaking more water from the canopy leaves, an effect that is not present in our leafless stand. Instead, we believe that as umax increases, precipitation falls at a lower angle relative to the horizontal, increasing the effective canopy area and allowing the canopy to intercept a greater percentage of the Pg (Herwitz and Slye 1995). Theoretically, the wind effects should be represented as an increased value of c (decreased value of ct). The WiMo model varies Sc, not c, with umax, and so the WiMo model represents the effect of wind speed on interception as an increased value of Sc for a constant c. The logarithmic relationship between Sc and umax is physically reasonable: At low (high) wind speeds, a small increase in the wind speed will have a large (negligible) effect on the angle at which the rain falls relative to the horizontal.
The regression equation relating Sc and umax has a nonsensible boundary condition, as Sc is negative for storms with umax less than 3.1 m s−1. None of the 11 storms considered had a umax below 3.1 m s−1, and so the regression equation in the Hörmann et al. (1996) bucket model fails to properly account for them. Hörmann et al. (1996) experienced similar problems, as their Sc grew improbably large for low values of umax. Evidently, storms with umax less than 3.1 m s−1 are quite rare in this stand during the leafless season, and so this boundary condition is not of great concern.
The Mulder model predicted the total interception more accurately than either the Gash or the WiMo model, and so might be considered the best interception model for this canopy under these conditions. While Mulder (1985) provides a different parameterization than Gash et al. (1995), it is also important to note that the Mulder model maintains a running water balance between storms, while the Gash and WiMo models do not. It is possible that the Gash and WiMo models’ assumption that the canopy is always completely dry prior to each storm is not valid in this stand. An experiment that compared these models’ performance for a set of storms that occurred immediately after one another, such that the canopy would not have the opportunity to dry, would highlight the effects of this assumption.
While all three models predicted the percentage of PG intercepted to within 1.1% of PG of the measured value, it may be that their accuracy can be attributed to the particular set of rainstorms in this study. The models frequently overestimated or underestimated interception for individual events (Table 4). Only 23 December, 27 March, and 1 April were predicted with any real accuracy. The standard deviation of each model’s differences for individual storms is approximately 1 mm. This is quite high considering that the least accurate model overall (the Gash model) underestimated the total precipitation by only 3.4 mm. The WiMo model had the lowest standard deviation and the highest kurtosis value. These statistics indicate that the distribution of the WiMo model’s differences had a higher peak around zero, where zero is the measured value. The Mulder and Gash models’ differences, however, were more uniformly distributed across their space. The WiMo model also had the highest Nash–Sutcliffe model efficiency. The Nash–Sutcliffe efficiency compares the variance of the residuals (i.e., the variance of the difference between measured and modeled values) against the initial variance of the measured values. An efficiency score of unity represents identical measured and modeled values. A high score for the WiMo model indicates that the WiMo model’s estimates for individual storms were, when considered together, reasonably close to the measured values. The Mulder model had the lowest Nash–Sutcliffe efficiency, the highest standard deviation, and a low kurtosis. These metrics imply that the Mulder model’s accurate estimate of the total interception was a result of the offsetting of its inaccuracies for individual storms over the course of the study.
The WiMo model’s ability to optimize Sc improves its predictions of interception, both for individual storms and the seasonal total, but restricts its use in stands where wind measurements are difficult to obtain. The WiMo model thus represents a definite improvement over the Gash model in this stand, but its applicability to other leafless stands is predicated upon the availability of wind data and the determination of a suitable regression equation between Sc and umax.
As mentioned previously, the measured percentage of Pg intercepted varied greatly from storm to storm. The models overpredicted interception when this percentage was much lower than the fraction of canopy covered (c) and underpredicted interception when the percentage was much higher than c. Over all 11 storms, the canopy intercepted 5.8% of PG, which is lower than c (10%) and closer to the models’ predictions. (The canopy should intercept somewhat less than c, since the canopy does not intercept all incident precipitation while saturated.) In other words, the percentage of Pg intercepted converged to a value less than c over the course of the study period. This is evident from Fig. 2. The fixed coefficient of free throughfall (ct)—and consequently c [Eq. (4)]—constrained the models’ ability to accurately determine interception during storms such as 9 February, as the modeled canopy could not possibly intercept more precipitation than the models allowed to strike it [Eq. (5)]. The ct establishes a hard upper boundary on interception in all three models. This limit does not exist in reality, as wind-driven events can increase the effective canopy area and allow the canopy to intercept more than Pυ (Herwitz and Slye 1995).
The models accurately predicted the total interception over the study period, but not the interception for the individual storms. The key difference between the total interception and the individual storms is that the measured percentage of Pg intercepted was less than (greater than) c for the total interception (individual storms). The percentage of PG intercepted is the critical value here, since these models are most often used to estimate interception over an entire season rather than for a single storm. Since all three models predicted the percentage of PG intercepted to within the probable standard error of the measurements, we believe that all three models will perform adequately for a leafless stand of this type in this climate, as long as ct is appropriately specified.
7. Conclusions
This study compared the Mulder (1985), Gash et al. (1995), and WiMo (Hörmann et al. 1996) interception models for a leafless mixed deciduous forest in the eastern United States. Each model predicted interception for 11 rainstorms during the winter and early spring of 2004/05. These values were compared against measurements collected approximately 2 h after each storm.
All three models demonstrated large differences for individual storms (Table 4) that bore no relationship to any applicable meteorological variable (Table 1), but were connected to the measured percentage of Pg intercepted. The models underestimated interception when the percentage of Pg intercepted was much greater than the canopy cover c (10%) and overestimated interception when the percentage of Pg was much less than c. While the percentage of Pg intercepted varied for individual events—leading to the models’ differences—it converged to a value less than c as the sample size of storms grew, which allowed the models to accurately predict the total measured interception. This underscores the importance of selecting an appropriate value of c (or ct) for these models.
The Mulder model overestimated the measured percentage of PG intercepted by 0.1% of PG; the WiMo model underestimated the measured percentage by 0.6% of PG; and the Gash model underestimated by 1.1% of PG. The WiMo model’s advantage over the Gash model is most likely due to the former’s ability to optimize the canopy storage capacity for the maximum sustained wind speed during a storm. In this stand, the storage capacity increased logarithmically with the maximum wind speed. The Mulder model gives the most accurate value of total interception, although its high standard deviation and low kurtosis imply that it was less accurate for individual storms. We hypothesize that the Mulder model’s accuracy is due at least in part to its ability to maintain a moist canopy between storms.
Despite differences in underlying assumptions and complexity, each model successfully predicted the percentage of PG intercepted to within the probable standard error of our measurements. We therefore conclude that all three models are appropriate for use in predicting leafless-season interception in deciduous forests in the eastern United States.
Acknowledgments
The authors thank the Maryland Department of Natural Resources for permission to conduct the study on their property and the Delaware Environmental Observing System (DEOS) for providing the meteorological monitoring equipment. Kevin Brinson assisted us in collecting and formatting the meteorological data. The Maryland State Climatologist’s Office (http://www.atmos.umd.edu/~climate) provided climate data for a location near the Fair Hill Natural Resource Management Area. We would also like to thank our three anonymous reviewers for their comments and suggestions.
REFERENCES
Abrahams, A. D., Parsons A. J. , and Wainwright J. , 2003: Disposition of rainwater under creosote-bush. Hydrol. Processes, 17 , 2555–2566.
Bruijnzeel, L. A., and Wiersum K. F. , 1987: Rainfall interception by a young Acacia auriculiformis (A. Cunn) plantation forest in West Java, Indonesia: Application of Gash’s analytical model. Hydrol. Processes, 1 , 309–319.
Bryant, M. L., Bhat S. , and Jacobs J. M. , 2005: Measurements and modeling of throughfall variability for five forest communities in the southeastern US. J. Hydrol., 312 , 95–108.
Calder, I. R., 1976: The measurement of water losses from a forested area using a “natural” lysismeter. J. Hydrol., 30 , 311–325.
Calder, I. R., 1977: A model of transpiration and interception loss from a spruce forest in Plynlimon, central Wales. J. Hydrol., 33 , 247–265.
Calder, I. R., 1990: Evaporation in the Uplands. Wiley, 148 pp.
Carlyle-Moses, D. E., and Price A. G. , 1999: An evaluation of the Gash interception model in a northern hardwood stand. J. Hydrol., 214 , 103–110.
Carlyle-Moses, D. E., Flores Laureano J. S. , and Price A. G. , 2004: Throughfall and throughfall spatial variability in Madrean oak forest communities of northeastern Mexico. J. Hydrol., 297 , 124–135.
Delage, Y., Bartlett P. A. , and McCaughey J. H. , 2002: Study of “soft” night-time surface-layer decoupling over forest canopies in a land-surface model. Bound.-Layer Meteor., 103 , 253–276.
Derbyshire, S. H., 1999: Boundary-layer decoupling over cold surfaces as a physical boundary-instability. Bound.-Layer Meteor., 90 , 297–325.
Dolman, A. J., 1987: Summer and winter rainfall interception in an oak forest: Predictions with an analytical and a numerical simulation model. J. Hydrol., 90 , 1–9.
Domingo, F., Sanchez G. , Moro M. J. , Brenner A. J. , and Puigdefabregas J. , 1998: Measurement and modelling of rainfall interception by three semi-arid canopies. Agric. For. Meteor., 91 , 275–292.
Garratt, J. R., and Francey R. J. , 1978: Bulk characteristics of heat transfer in the baroclinic atmospheric boundary layer. Bound.-Layer Meteor., 15 , 399–421.
Gash, J. H. C., 1979: An analytical model of rainfall interception by forests. Quart. J. Roy. Meteor. Soc., 105 , 43–55.
Gash, J. H. C., Lloyd C. R. , and Lachaud G. , 1995: Estimating sparse forest rainfall interception with an analytical model. J. Hydrol., 170 , 79–86.
Halldin, S., Saugier B. , and Pontailler J. Y. , 1984: Evapotranspiration of a deciduous forest: Simulation using routine meteorological data. J. Hydrol., 75 , 323–341.
Harwood, K. G., Gillon J. S. , Roberts A. , and Griffiths H. , 1999: Determinants of isotopic coupling CO2 and water vapour within a Quercus petraea forest canopy. Oecologia, 119 , 109–119.
Herwitz, S. R., and Slye R. E. , 1995: Three-dimensional modeling of canopy tree interception of wind-driven rainfall. J. Hydrol., 168 , 205–226.
Hörmann, G., Branding A. , Clemen T. , Herbst M. , Hinrichs A. , and Thamm F. , 1996: Calculation and simulation of wind controlled canopy interception of a beech forest in northern Germany. Agric. For. Meteor., 79 , 131–148.
Hutjes, R. W. A., Wierda A. , and Veen A. W. L. , 1990: Rainfall interception in the Tai Forest, Ivory Coast: Application of two simulation models to a humid tropical system. J. Hydrol., 114 , 259–275.
Jackson, I. J., 1971: Problems of throughfall and interception assessment under tropical forest. J. Hydrol., 12 , 234–254.
Kittredge, J., 1948: Forest Influences. McGraw-Hill, 394 pp.
Kruijt, B., Malhi Y. , Lloyd J. , Nobre A. D. , Miranda A. C. , Pereira M. G. P. , Culf A. , and Grace J. , 2000: Turbulence statistics above and within two Amazon rain forest canopies. Bound.-Layer Meteor., 94 , 297–331.
Lankreijer, H. J. M., Hendricks M. J. , and Klassen W. , 1993: A comparison of models simulating rainfall interception of forests. Agric. For. Meteor., 64 , 187–199.
Levia, D. F., and Frost E. E. , 2003: A review and evaluation of stemflow literature in the hydrologic and biogeochemical cycles of forested and agricultural ecosystems. J. Hydrol., 274 , 1–29.
Liu, S. G., 1997: A new model for the prediction of rainfall interception in forest canopies. Ecol. Modell., 99 , 151–159.
Liu, S. G., 1998: Estimation of rainfall storage capacity in the canopies of cypress wetlands and slash pine uplands in north-central Florida. J. Hydrol., 207 , 32–41.
Liu, S. G., 2001: Evaluation of the Liu model for predicting rainfall interception in forests world-wide. Hydrol. Processes, 15 , 2341–2360.
Lloyd, C. R., and Marques Ade O. , 1988: Spatial variability of throughfall and stemflow measurements in Amazonian rainforest. Agric. For. Meteor., 42 , 63–73.
Lloyd, C. R., Gash J. H. C. , Shuttleworth W. J. , and Marques Ade O. , 1988: The measurement and modelling of rainfall interception by Amazonian rainforest. Agric. For. Meteor., 43 , 277–294.
Loescher, H. W., Powers J. S. , and Oberbauer S. F. , 2002: Spatial variation of throughfall volume in an old-growth tropical wet forest, Costa Rica. J. Trop. Ecol., 18 , 397–407.
Loustau, D., Berbigier P. , and Granier A. , 1992: Interception loss, throughfall, and stemflow in a maritime pine stand. II. An application of Gash’s analytical model of interception. J. Hydrol., 138 , 469–485.
Moore, K., Fitzjarrald D. R. , and Sakai R. K. , 1996: Seasonal variation in radiative and turbulent exchange at a deciduous forest in central Massachussetts. J. Appl. Meteor., 35 , 122–134.
Mulder, J. P. M., 1985: Simulating interception loss using standard meteorological data. The Forest-Atmosphere Interaction, B. Hutchison and B. Hicks, Eds., D. Reidel, 177–196.
Nash, J. E., and Sutcliffe J. V. , 1970: River flow forecasting through conceptual models. Part I: A discussion of principles. J. Hydrol., 10 , 282–290.
Pearce, A. J., and Rowe L. K. , 1981: Rainfall interception in a multi-storied evergreen mixed forest: Estimates using Gash’s analytical model. J. Hydrol., 49 , 341–353.
Rutter, A. J., 1967: Evaporation in forests. Endeavour, 26 , 39–43.
Rutter, A. J., Kershaw K. A. , Robbins P. C. , and Morton A. J. , 1971: A predictive model of rainfall interception in forests. I. Derivation of the model from observations in a plantation of Corsican pine. Agric. For. Meteor., 9 , 367–384.
Rutter, A. J., Morton A. J. , and Robbins P. C. , 1975: A predictive model of rainfall interception in forests. II. Generalization of the model and comparison with observations in some coniferous and hardwood stands. J. Appl. Ecol., 12 , 367–380.
Stewart, J. B., 1977: Evaporation from the wet canopy of a pine forest. Water Resour. Res., 13 , 915–921.
Stewart, J. B., 1978: A micrometeorological investigation into factors controlling the evaporation from a forest. Ph.D. thesis, University of Reading, 211 pp.
Teklehaimanot, Z., Jarvis P. G. , and Ladger D. C. , 1991: Rainfall interception and boundary layer conductance in relation to tree spacing. J. Hydrol., 123 , 261–278.
Thom, A. S., 1971: Momentum absorption by vegetation. Quart. J. Roy. Meteor. Soc., 97 , 414–428.
Wei, X., Liu S. , and Wang C. , 2005: Hydrological processes in major types of Chinese forest. Hydrol. Processes, 19 , 63–75.
Whelan, M. J., and Anderson J. M. , 1996: Modelling spatial patterns of throughfall and interception loss in a Norway spruce (Picea abies) plantation at the plot scale. J. Hydrol., 186 , 335–354.

A map showing the location of Fair Hill Natural Resource Management Area in the context of the mid-Atlantic United States.
Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM564.1

A map showing the location of Fair Hill Natural Resource Management Area in the context of the mid-Atlantic United States.
Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM564.1
A map showing the location of Fair Hill Natural Resource Management Area in the context of the mid-Atlantic United States.
Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM564.1

Cumulative totals of measured and modeled interception (mm) through the study period.
Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM564.1

Cumulative totals of measured and modeled interception (mm) through the study period.
Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM564.1
Cumulative totals of measured and modeled interception (mm) through the study period.
Citation: Journal of Hydrometeorology 8, 4; 10.1175/JHM564.1
A summary of meteorological conditions during each rainstorm: total precipitation (Pg, mm), precipitation rate (


Values held constant in this study.


A comparison of the fundamental equations for each model. The Gash and WiMo models dry out completely after each storm and include a separate term for the “interception” during the drying-out phase. The Mulder model does not account for stemflow.


Measured interception data as compared to the results from each model. Data for individual events are given as the percentage of Pg intercepted. Differences are taken as the measured value minus the modeled value. The statistics for individual differences were calculated using the differences between the measured and modeled values in millimeters. The Nash–Sutcliffe model efficiency does not include the 9 February storm.


A sensitivity test for the three models used in this study. All values are given as ratios of calculated interception. A ratio of 1.000 indicates that the model is perfectly insensitive to the given change in the given parameter. Parameters that do not apply to a model are marked n/a.

