Regional differences in hydrological response to canopy interception schemes in a land surface model

Two different canopy interception schemes are applied to the parameterization of the hydrological processes in the Community Land Model version 3. One scheme treats rainfall and canopy water storage as spatially uniform within each model grid cell, and the other scheme considers sub‐grid variability of rainfall and water storage in the parameterization of canopy hydrological processes. The hydrological responses to differences between these two schemes in different regions are studied. It is found that the impact of the sub‐grid variability in the tropical regions is generally greater than the extra‐tropical regions. However, such impact can't be negligible for the extra‐tropical regions. Soil water in the total 3.4 m soil depth increases by 3% for Central‐South Europe, and vegetation temperature increases by 0.14 °C for Southeastern United States if the regional averages are considered. The magnitude of the impact is greater if the analysis focuses on the specific grid cells in these regions. The impact is tightly correlated with rainfall amount and vegetation density. The correlation coefficient between such impact and rainfall amount and vegetation density varies with regions and hydrological variables, with the largest value of 0.92 for interception loss in Amazonia. Our results indicate that the impact of the sub‐grid variability on hydrological processes in the extra‐tropical areas is also important, although rainfall amount and vegetation density in these areas are not as high as in the tropical areas. Copyright © 2013 John Wiley & Sons, Ltd.     

Rainfall interception in a lower montane forest in Ecuador: effects of canopy properties

Rainfall interception in forests is influenced by properties of the canopy that tend to vary over small distances. Our objectives were: (i) to determine the variables needed to model the interception loss of the canopy of a lower montane forest in south Ecuador, i.e. the storage capacity of the leaves S and of the trunks and branches S_t, and the fractions of direct throughfall p and stemflow p_t; (ii) to assess the influence of canopy density and epiphyte coverage of trees on the interception of rainfall and subsequent evaporation losses. The study site was located on the eastern slope of the eastern cordillera in the south Ecuadorian Andes at 1900–2000 m above sea level. We monitored incident rainfall, throughfall, and stemflow between April 1998 and April 2001. In 2001, the leaf area index (LAI), inferred from light transmission, and epiphyte coverage was determined. The mean annual incident rainfall at three gauging stations ranged between 2319 and 2561 mm. The mean annual interception loss at five study transects in the forest varied between 591 and 1321 mm, i.e. between 25 and 52% of the incident rainfall. Mean S was estimated at 1·91 mm for relatively dry weeks with a regression model and at 2·46 mm for all weeks with the analytical Gash model; the respective estimates of mean S_t were 0·04 mm and 0·09 mm, of mean p were 0·42 and 0·63, and of mean p_t were 0·003 and 0·012. The LAI ranged from 5·19 to 9·32. Epiphytes, mostly bryophytes, covered up to 80% of the trunk and branch surfaces. The fraction of direct throughfall p and the LAI correlated significantly with interception loss (Pearson's correlation coefficient r = −0·77 and 0·35 respectively, n = 40). Bryophyte and lichen coverage tended to decrease S_t and vascular epiphytes tended to increase it, although there was no significant correlation between epiphyte coverage and interception loss. Our results demonstrate that canopy density influences interception loss but only explains part of the total variation in interception loss. Copyright © 2004 John Wiley & Sons, Ltd.     

Effects of changes in canopy interception on stream runoff response and recovery following clear-cutting of a Japanese coniferous forest in Fukuroyamasawa Experimental Watershed in Japan

Understanding changes in evapotranspiration during forest regrowth is essential to predict changes of stream runoff and recovery after forest cutting. Canopy interception (Ic) is an important component of evapotranspiration, however Ic changes and the impact on stream runoff during regrowth after cutting remains unclear due to limited observations. The objective of this study was to examine the effects of Ic changes on long-term stream runoff in a regrowth Japanese cedar and Japanese cypress forest following clear-cutting. This study was conducted in two 1-ha paired headwater catchments at Fukuroyamasawa Experimental Watershed in Japan. The catchments were 100% covered by Japanese coniferous plantation forest, one of which was 100% clear-cut in 1999 when the forest was 70 years old. In the treated catchment, annual runoff increased by 301 mm/year (14% of precipitation) the year following clear-cutting, and remained 185 mm/year (7.9% of precipitation) higher in the young regrowth forest for 12–14 years compared to the estimated runoff assuming no clear-cutting. The Ic change was −358 mm/year (17% of precipitation) after cutting and was −168 mm/year (6.7% of precipitation) in the 12–14 years old regrowth forest compared to the observed Ic during the pre-cutting period. Stream runoff increased in all seasons, and the Ic change was the main fraction of evapotranspiration change in all seasons throughout the observation period. These results suggest that the change in Ic accounted for most of the runoff response following forest cutting and the subsequent runoff recovery in this coniferous forest.     

Hydrochemical characteristics of throughfall and stemflow in a Moso‐bamboo (Phyllostachys pubescens) forest

To investigate the impacts of the invasion by bamboo on fluxes of nutrients and pollutants, the nutrient/pollutant fluxes and canopy interactions, including neutralization of acidity, leaching and uptake of nitrogen (N), were characterized in conjunction with rainfall partitioning in a Moso‐bamboo (Phyllostachys pubescens) forest. Measurements of precipitation volume, pH, major ions, and silicate (SiO₂) in rainfall, throughfall and stemflow were collected weekly in a Moso‐bamboo forest located in Munakata City, Western Japan for 1 year. Results showed that rainfall partitioning into stemflow was larger than that for other types of forest, which may be due to the properties of Moso‐bamboo forest structure, such as a straight and smooth culm. Inorganic N (NO₃⁻ + NH₄⁺) and S (SO₄²⁻) fluxes of throughfall and stemflow were approximately 1·6 and 1·3 times higher than that of rainfall, respectively. Contribution of stemflow flux to inorganic N and S fluxes to the forest floor was high. This could be due to lower uptake of inorganic N through culm and a higher rainfall partitioning into stemflow than that for other types of forest. The Moso‐bamboo canopy neutralized rainfall acidity, reducing the fluxes of potentially acidifying compounds via throughfall and stemflow. Canopy leaching of K⁺ was distinctly higher than that of Mg²⁺ and Ca²⁺ and could be related to the high mobility of K⁺ in plant tissues. Cl⁻ and SiO₂ were readily leached as for K⁺. The impact of the invasion by bamboo on nutrient cycling was discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

Spatial variations in throughfall in a Moso bamboo forest: sampling design for the estimates of stand‐scale throughfall

We investigated the spatial and seasonal variations in throughfall (Tf) in relation to spatial and seasonal variations in canopy structure and gross rainfall (Rf) and assessed the impacts of the variations in Tf on stand‐scale Tf estimates. We observed the canopy structure expressed as the leaf area index (LAI) once a month and Tf once a week in 25 grids placed in a Moso bamboo (Phyllostachys pubescens) forest for 1 year. The mean LAI and spatial variation in LAI did have some seasonal variations. The spatial variations in Tf reduced with increasing Rf, and the relationship between the spatial variation and the Rf held throughout the year. These results indicate that the seasonal change in LAI had little impact on spatial variations in Tf, and that Rf is a critical factor determining the spatial variations in Tf at the study site. We evaluated potential errors in stand‐scale Tf estimates on the basis of measured Tf data using Monte Carlo sampling. The results showed that the error decreases greatly with increasing sample size when the sample size was less than ～8, whereas it was near stable when the sample size was 8 or more, regardless of Rf. A sample size of eight results in less than 10% error for Tf estimates based on Student's t‐value analysis and would be satisfactory for interception loss estimates when considering errors included in Rf data. Copyright © 2009 John Wiley & Sons, Ltd.     

Spatial and temporal variability of canopy and forest floor interception in a beech forest

Depending on season, rainfall characteristics and tree species, interception amounts to 15–50% of total precipitation in a forest under temperate climates. Many studies have investigated the importance of interception of different tree species in all kinds of different climates. Often authors merely determine interception storage capacity of that specific species and the considered event, and only sometimes a distinction is made between foliated and non‐foliated trees. However, interception is highly variable in time and space. First, since potential evaporation is higher in summer, but secondly because the storage capacity has a seasonal pattern. Besides weather characteristics, such as wind and rain intensity, snow causes large variations in the maximum storage capacity. In an experimental beech plot in Luxembourg, we found storage capacity of canopy interception to show a clear seasonal pattern varying from 0·1 mm in winter to 1·2 mm in summer. The capacity of the forest floor appears to be rather constant over time at 1·8 mm. Both have a standard deviation as high as ± 100%. However, the process is not sensitive to this variability resulting only in 11% variation of evaporation estimates. Hence, the number of raindays and the potential evaporation are stronger driving factors on interception. Furthermore, the spatial correlation of the throughfall and infiltration has been investigated with semi‐variograms and time stability plots. Within 6–7 m distance, throughfall and infiltration are correlated and the general persistence is rather weak. Copyright © 2010 John Wiley & Sons, Ltd.     

A model to consider the spatial variability of rainfall partitioning within deciduous canopy. I. Model description

In this first paper of two, a numerical simulation model capable of simulating the spatial variability of rainfall partitioning within a canopy is presented. The second paper details how the model is parameterized, and some testing of its capabilities. The first stage of the model is to derive the mature canopy structure. This is achieved through simplified individual tree structures and a random placement routine based on a modified Poisson distribution. Following this the spatial discretization for throughfall is attained as a series of layers of triangles with a tree trunk at every apex. The number of layers is derived from the leaf area index through a modified Poisson distribution. Seasonal variation in the deciduous canopy is simulated through a time vector routine. The rainfall partitioning section of the model is based upon the Rutter model which has been modified to each individual triangle layer. The main feature of this model is that it offers a method of simulating rainfall partitioning at a scale that is a function of the size of the elementary physical components (tree and leaf scale). This can be used to investigate soil moisture variations under a canopy, or to study the variations within the forest hydrological processes themselves. © 1997 John Wiley & Sons, Ltd.     

A new approach to estimate canopy evaporation and canopy interception capacity from evapotranspiration and sap flow measurements during and following wetting

Canopy interception and its evaporation into the atmosphere during irrigation or a rainfall event are important in irrigation scheduling, but are challenging to estimate using conventional methods. This study introduces a new approach to estimate the canopy interception from measurements of actual total evapotranspiration (ET) using eddy covariance and estimation of the transpiration from measurements of sap flow. The measurements were conducted over a small‐scale sprinkler‐irrigated cotton field before, during and after sprinkler irrigation. Evaporation and sap flow dynamics during irrigation show that the total ET during irrigation increased significantly because of the evaporation of free intercepted water while transpiration was suppressed almost completely. The difference between actual ET and transpiration (sap flow) during and immediately following irrigation (post irrigation) represents the total canopy evaporation while the canopy interception capacity was calculated as the difference between actual ET and transpiration (sap flow) during drying (post irrigation) following cessation of the irrigation. The canopy evaporation of cotton canopy was calculated as 0.8 mm, and the interception capacity was estimated to be 0.31 mm of water. The measurement uncertainty in both the non‐dimensional ET and non‐dimensional sap flow was shown to be very low. Copyright © 2015 John Wiley & Sons, Ltd.     

Rainfall interception by an isolated evergreen oak tree in a Mediterranean savannah

Redistribution of ground‐level rainfall and interception loss by an isolated Quercus ilex tree were measured over 2 years in a Mediterranean oak savannah. Stemflow, meteorological variables and sap flow were also monitored. Rainfall at ground level was measured by a set of rain‐gauges located in a radial layout centred on the tree trunk and extending beyond the crown limits. Interception loss was computed as the difference between the volume of rainwater that would reach the ground in the absence of the tree and the volume of water that actually fell on the ground sampling area (stemflow included). This procedure provided correct interception loss estimates, irrespective of rainfall inclination. Results have shown a clear non‐random spatial distribution of ground‐level rainfall, with rainwater concentrations upwind beneath the crown and rain‐shadows downwind. Interception loss amounted to 22% of gross rainfall, per unit of crown‐projected area. Stand interception loss, per unit of ground area, was only 8% of gross rainfall and 28% of tree evapotranspiration. These values reflect the low crown cover fraction of the stand (0·39) and the specific features of the Mediterranean rainfall regime (predominantly with few large storms). Nevertheless, it still is an important component of the water balance of these Mediterranean ecosystems. Copyright © 2005 John Wiley & Sons, Ltd.     

Observations of canopy storage capacity and wet canopy evaporation in a humid boreal forest

Evaporation of intercepted rain by a canopy is an important component of evapotranspiration, particularly in the humid boreal forest, which is subject to frequent precipitation and where conifers have a large surface water storage capacity. Unfortunately, our knowledge of interception processes for this type of environment is limited by the many challenges associated with experimental monitoring of the canopy water balance. The objective of this study is to observe and estimate canopy storage capacity and wet canopy evaporation at the sub-daily and seasonal time scales in a humid boreal forest. This study relies on field-based estimates of rainfall interception and evapotranspiration partitioning at the Montmorency Forest, Québec, Canada (mean annual precipitation: 1600 mm, mean annual evapotranspiration: 550 mm), in two balsam fir-white birch forest stands. Evapotranspiration was monitored using eddy covariance sensors and sap flow systems, whereas rainfall interception was measured using 12 sets of throughfall and six stemflow collectors randomly placed inside six 400-m² plots. Changes in the amount of water stored on the canopy were also directly monitored using the stem compression method. The amount of water intercepted by the forest canopy was 11 ± 5% of the total rainfall during the snow-free (5 July–18 October) measurement periods of 2017 and 2018. The maximum canopy storage estimated from rainfall interception measurements was on average 1.6 ± 0.7 mm, though a higher value was found using the stem compression method (2.2 ± 1.6 mm). Taking the average of the two forest stands studied, evaporation of intercepted water represented 21 ± 8% of evapotranspiration, while the contribution of transpiration and understory evapotranspiration was 36 ± 9% and 18 ± 8%. The observations of each of the evapotranspiration terms underestimated the total evapotranspiration observed, so that 26 ± 12% of it was not attributed. These results highlight the importance to account for the evaporation of rain intercepted by humid boreal forests in hydrological models.     

Modelling crop canopy and residue rainfall interception effects on soil hydrological components for semi‐arid agriculture

Crop canopies and residues have been shown to intercept a significant amount of rainfall. However, rainfall or irrigation interception by crops and residues has often been overlooked in hydrologic modelling. Crop canopy interception is controlled by canopy density and rainfall intensity and duration. Crop residue interception is a function of crop residue type, residue density and cover, and rainfall intensity and duration. We account for these controlling factors and present a model for both interception components based on Merriam's approach. The modified Merriam model and the current modelling approaches were examined and compared with two field studies and one laboratory study. The Merriam model is shown to agree well with measurements and was implemented within the Agricultural Research Service's Root Zone Water Quality Model (RZWQM). Using this enhanced version of RZWQM, three simulation studies were performed to examine the quantitative effects of rainfall interception by corn and wheat canopies and residues on soil hydrological components. Study I consisted of 10 separate hypothetical growing seasons (1991–2000) for canopy effects and 10 separate non‐growing seasons (1991–2000) for residue effects for eastern Colorado conditions. For actual management practices in a no‐till wheat–corn–fallow cropping sequence at Akron, Colorado (study II), a continuous 10‐year RZWQM simulation was performed to examine the cumulative changes on water balance components and crop growth caused by canopy and residue rainfall interception. Finally, to examine a higher precipitation environment, a hypothetical, no‐till wheat–corn–fallow rotation scenario at Corvallis, Oregon, was simulated (study III). For all studies, interception was shown to decrease infiltration, runoff, evapotranspiration from soil, deep seepage of water and chemical transport, macropore flow, leaf area index, and crop/grain yield. Because interception decreased both infiltration and soil evapotranspiration, no significant change in soil water storage was simulated. Nonetheless, these findings and the new interception models are significant new contributions for hydrologists. Published in 2006 John Wiley & Sons, Ltd.     

A semi‐theoretical model of canopy rainfall interception for a broad‐leaved tree

Through a series of simulation experiments in the laboratory on the broad‐leaved tree Acer mono Maxim, we obtain interception datasets of individual events under different rainfall intensities and leaf area indexes (LAIs). Based on the data, the relationship between rainfall intensity and maximum interception of per unit LAI is quantified. The variation of interception with canopy wetness index is also identified. Hence, an interception model, in which interception is calculated using rainfall intensity and LAI, is constructed with consideration of canopy wetness. Finally, according to the validation experiments, it is concluded that the precision of the model is 92·7%. Copyright © 2007 John Wiley & Sons, Ltd.     

Measuring interception loss and canopy storage in dryland vegetation: a brief review and evaluation of available research strategies

The interception storage capacity has been measured for a range of dryland plants. Interception losses over time, however, arise in rain events that deliver either less or more than the canopy capacity. The fate of water in these cases depends on the efficiency with which the intercepted water is returned to the atmosphere by evaporation from the plant canopies. Two primary methods to estimate interception losses are (i) calibrated process‐based models of interception and evaporative loss and (ii) direct measurement. Models have been applied only rarely to dryland plant communities, and direct measurement techniques are in need of additional testing and refinement. Most published estimates of interception loss in dryland plant communities therefore appear to be based upon inadequate data and methods. Research needs in this area are highlighted. Copyright © 2000 John Wiley & Sons, Ltd.     

Effects of phenology and meteorological disturbance on litter rainfall interception for a Pinus elliottii stand in the Southeastern United States

Litter layers develop across a diverse array of vegetated ecosystems and undergo significant temporal compositional changes due to canopy phenological phases and disturbances. Past research on temporal dynamics of litter interception has focused primarily on litter thickness and leaf fall, yet forest phenophases can change many more litter attributes (e.g., woody debris, bark shedding, and release of reproductive materials). In this study, weekly changes in litter composition over 1 year were used to estimate litter water storage dynamics and model event‐based litter interception. Litter interception substantially reduced throughfall (6–43%), and litter water storage capacity ranged from 1 to 3 mm, peaking when megastrobili release and liana leaf senescence occurred simultaneously during fall 2015. Tropical storm disturbances occurred during the sampling period, allowing evaluation of how meteorological disturbances altered litter interception. High wind speeds and intense rainfall from 2 tropical storms increased litter interception by introducing new woody debris, which, in this study, stored more water than the pre‐existing woody debris. After 2 extreme weather events, a third (Hurricane Hermine) did not increase woody debris (or litter interception), suggesting that the canopy pool of branches susceptible to breakage had been largely depleted. Needle and bark shedding had minor effects on litter interception. Results suggest that the release of reproductive materials and meteorological disturbances appear to be the major compositional drivers of litter interception beyond their obvious contribution to litter thickness.     

Throughfall and stemflow vary seasonally in different land‐use types in a lower montane tropical region of Panama

Catchment hydrology is influenced by land‐use change through alteration of rainfall partitioning processes. We compared rainfall partitioning (throughfall, stemflow and interception) and soil water content in three land‐use types (primary forest, secondary forest and agriculture) in the Santa Fe region of Panama. Seasonal patterns were typified by larger volumes of throughfall and stemflow in the wet season, and the size of precipitation events was the main driver of variation in rainfall redistribution. Land‐use‐related differences in rainfall partitioning were difficult to identify due to the high variability of throughfall. However, annual throughfall in agricultural sites made up a larger proportion of gross precipitation than throughfall in forest sites (94 ± 1, 83 ± 6 and 81 ± 1% for agriculture, primary and secondary forests, respectively). Proportional throughfall (% of gross precipitation becoming throughfall) was consistent throughout the year for primary forest, but for secondary forest, it was larger in the dry season than the wet season. Furthermore, proportional stemflow in the dry season was larger in secondary forest than primary forest. Stemflow, measured only in primary and secondary forests, ranged between 0.9 and 3.2% of gross precipitation. Relative soil moisture content in agricultural plots was generally elevated during the first half of the dry season in comparison to primary and secondary forests. Because throughfall is elevated in agricultural plots, we suggest careful management of the spatial distribution and spread of this land‐use type to mitigate potential negative impacts in the form of floods and high erosion rates in the catchment. Copyright © 2013 John Wiley & Sons, Ltd.     

Factors affecting rainfall interception determined by a forest simulator and numerical model

To elucidate the factors involved in interception loss, we conducted experiments in which we measured environmental variables such as rainfall intensity, forest structure, and weather conditions. An artificial forest consisting of 24 vinyl trees was used to examine the influences of forest structure and rainfall conditions on interception loss. The interception rate was higher at higher plant area index (PAI) values and wind speeds and lower with greater rainfall intensity. We confirmed the factors affecting interception loss by using an interception model based on the tank model. The artificial forest simulations provide new evidence that interception loss is influenced by the PAI, rainfall intensity, saturation deficit, and wind speed. The effect of the saturation deficit on the interception rate was unclear from the experimental results, but the single‐tank model revealed that wind speed strongly influences the effects of the saturation deficit on interception loss. Thus, whereas interception loss was not significantly affected by the saturation deficit at low wind speeds, it increased significantly with the saturation deficit under windy conditions. The model simulation also showed the sensitivity of each factor with regard to the interception rate. The sensitivity of rainfall intensity decreased as the PAI increased, and the sensitivity of the saturation deficit increased as the wind speed increased. The experiments and model calculations clarified the main elements affecting interception loss and their sensitivities. Compared with previous studies on interception loss, this study revealed a positive relationship between the PAI and interception loss, a negative exponential relationship with rainfall intensity, and the effects of the saturation deficit on interception loss. Copyright © 2008 John Wiley & Sons, Ltd.     

Intra‐storm evaporation as a component of canopy interception loss in dryland shrubs: observations from Fowlers Gap,Australia

Interception losses from the canopies of dryland plant taxa remain poorly understood, especially the relative contributions of intra‐storm and post‐storm evaporative losses. Employing a new measuring apparatus, this study uses low‐intensity simulated rain, matched to the properties of local rain, to explore interception processes in bluebush shrubs at an Australian dryland site. Five shrub specimens were exposed to simulated rain for 60–90 min. Experiments were repeated at three rainfall intensities (10, 15, and 20 mm h^?1). Canopy evaporation was found from the difference between the flux of water delivered to the shrub and the flux of throughfall, once equilibrium had been established. The results show that evaporation from the wet foliage during rain proceeds at an average rate of 3·6 mm h^?1. This figure is for relatively cool spring‐season conditions; evaporation rates in hot summer conditions would be larger. Intra‐storm evaporation is shown to exceed post‐rain evaporation from interception storage on the shrubs, and this differentiates dryland shrub interception processes from those of the better‐studied wet forest environment. Implications of the high dryland shrub canopy evaporation rates for aspects of dryland ecology are highlighted. Copyright © 2007 John Wiley & Sons, Ltd.     

Canopy interception by a spruce forest in the upper reach of Heihe River basin,Northwestern China

The aim of this study is to understand the canopy interception of Qinghai spruce forest under conditions of different precipitation characteristics and canopy structures in the upper reach of Heihe River basin, northwestern China. On the basis of a continuous record covering our investigating period by an automatic throughfall‐collecting system, we analysed the relationships between the canopy interception and the precipitation characteristics. Our results support the well‐established exponential decay relationship between the gross precipitation and the interception percentage after the canopy is saturated. But our results sufficiently illustrate a notable point that the variations in the interception percentage are almost independent from the variations in the gross precipitation before the canopy is saturated. Our examination into the relationship between the interception and the 10‐min average intensity of precipitation demonstrates a divergent relationship, and the divergent relationship is bracketed by an upper ‘dry line’ indicating that 100% of gross precipitation was intercepted before saturation and by a lower ‘wet line’ suggesting that the actual canopy storage capacity reached the maximum and evaporation was the only component of the interception. To search for the relationship between canopy structures and interception, we grouped the canopy covers over the 90 throughfall‐collecting tanks into ten categories ranging from 0 (no cover) to 0.9 (nearly completely covered), and the corresponding canopy interception was calculated by subtracting the averaged throughfall of each canopy‐cover category from the gross precipitation. The results show that the interception percentage increases faster with increasing canopy cover under intermediate rainfall conditions than that under heavy rainfall conditions. Unexpectedly, under light rainfall conditions the increasing rate of interception percentage with increasing canopy cover and also with increasing plant area index is not faster than that under the intermediate rainfall conditions simply because the tank‐measured percentage of interception was extremely high at near‐zero canopy cover conditions. Copyright © 2013 John Wiley & Sons, Ltd.     

Canopy interception of apple orchards should not be ignored when assessing evapotranspiration partitioning on the Loess Plateau in China

Canopy interception is one of the most important processes in an ecosystem, but it is still neglected when assessing evapotranspiration (ET) partitioning in apple orchards on the Loess Plateau in China. To explore the importance of canopy interception, we monitored two neighbouring apple orchards on the Loess Plateau in China, one 8‐year‐old and the other 18‐years old at the start of the study, from May to September for four consecutive years (2013–2016). We measured parameters of canopy interception (I) including precipitation, throughfall, stemflow, leaf area index, transpiration (T), and soil evaporation (S) to quantify ET. The importance of canopy interception was then assessed by comparing the relationship between water supply (precipitation) and water demand (ET), calculated with and without considering canopy interception (T + S and T + S + I, respectively). Tree age clearly influenced canopy interception, as estimates of annual canopy interception during the study years in the younger and older orchards amounted to 22.2–29.4 mm and 26.8–39.9 mm, respectively. Daily incident rainfall and rainfall intensity in both orchards were significantly positively correlated with daily canopy interception in each year. The relationship between annual precipitation and annual ET (calculated with and without consideration of canopy interception) in the younger orchard differed during 2015 and 2016. Ignoring canopy interception would result in underestimation of annual ET in both apple orchards and hence incorrect evaluation of the relationship between water supply and water demand, particularly for the younger orchard during 2015 and 2016. Thus, for a complete understanding of water consumption in apple orchards in this and similar regions, canopy interception should not be ignored when assessing ET partitioning.