## Thursday, 17 November 2016

### Effect of the litter layer on runoff and evapotranspiration using the paired watershed method

Published Date
Volume 21, Issue 6pp 306–313

Author
Open AccessShort Communication
DOI: 10.1007/s10310-016-0542-5

Gomyo, M. & Kuraji, K. J For Res (2016) 21: 306. doi:10.1007/s10310-016-0542-5

Abstract

The litter layer on a forest floor can influence both short-term runoff and long-term water balance through modification of various hydrological processes. In this study, we have quantified the watershed-scale effects of the litter layer on runoff and evapotranspiration using a paired watershed method. The removal of the litter layer in a forested watershed with an area of 1.19 ha was conducted annually over the latter half of a 6-year experimental period. An adjacent forested watershed with an area of 1.42 ha was preserved as a control. Our results indicated that litter removal increased the 3-year runoff by 80.3 mm during the post-treatment period. Furthermore, when the peak flow range in the control watershed was 0.4–1.0 mm/h and >1.0 mm/h, peak runoff during flood events was about 1.5 and 1.4 times greater than that observed before litter layer removal, respectively. These data suggest that litter layer removal can decrease litter layer interception and, hence, increase peak flow, particularly during relatively large flood runoff events.

Introduction
There has been a recent decline in harvesting of natural and cultivated forests in developed countries due to the lack of sufficient economic incentive to produce conventional forest products (Onda et al. 2010; Molina and del Campo 2012). Consequently, their woody biomass is increasing and their litter layers are thickening. These forests can exacerbate evapotranspiration (ET) and contribute to a decrease in runoff from forested watersheds (Gallart and Llorens 2003; Komatsu et al. 2015). If this decrease occurs mainly during flood runoff events, flood mitigation could be positively affected; however, if the loss of runoff occurs mainly via the baseflow, this could have negative impacts on the management of water resources.
Three components of ET occur in forests: transpiration by vegetation (35–80 % of ET), rainfall interception by the canopy (10–50 % of precipitation), and forest floor evaporation (10–50 % of throughfall). Of these, forest floor evaporation, which consists of litter layer interception and soil surface evaporation, has received relatively little attention in the literature (Gerrits et al. 2007; Li et al. 2013). Quantifying forest floor evaporation and the effect that the litter layer has on runoff and ET is, therefore, important for managing forests, flood runoff mitigation, and water resources.
A litter layer is typically composed of dead leaves, twigs, small branches, and other fragmented organic material, and has various influences on the hydrological processes that operate in a forested watershed. For example, the litter layer intercepts throughfall and stemflow during periods of rainfall, and causes evaporation during and after them (Hattori 1993; Sato 2007; Gerrits and Savenije 2011; Li et al. 2013). The litter layer covers the ground surface and thus suppresses ground surface evaporation during dry periods (Majima et al. 1990; Murai 1993). In addition, it mitigates the impact energy of raindrops on throughfall, thereby increasing permeability and reducing overland flow (Walsh and Voigt 1977); however, it also creates a rapid-flow component within the litter layer (Okunishi 1963; Walsh and Voigt 1977; Lee and Shibano 1990). The simultaneous operation of these processes causes the litter layer to affect both short-term runoff and long-term water balance within the hydrological cycle.
Previous studies into the effects of the litter layer on runoff and water balance have been limited to individual hydrological processes; for example, by examining the water-holding capacity of the litter layer in order to estimate the factors that control litter interception (Sato et al. 2004; Rao and Zhu 2007; Zhang et al. 2009; Li et al. 2013). Continuous field measurements of litter interception can only be obtained using a sheet-shaped weighing lysimeter (Schaap and Bouten 1997; Gerrits et al. 20072010) or a permeable basin with a tipping bucket rain gauge to continuously monitor the water that drains from the litter layer into the soil (Bulcock and Jewitt 2012). However, these measuring devices are difficult to establish on sloping surfaces because the overlying litter can move downslope under the influence of gravity, which can alter the impact of raindrops on the forest floor. For this reason, all previous studies using these devices have implemented them on flat-lying ground. Moreover, uncertainties arise in extrapolating results obtained from a small flat area to an entire watershed, as the latter includes flat land, slopes, ridges, and the riparian zone. Micro-meteorological conditions, litterfall, and litter decomposition rate differ depending on the topography of each area. In addition, although the flow of streamwater away from a watershed can be measured, it is difficult to measure the groundwater flux from a small flat area to beyond its boundary. As far as the authors are aware, no watershed-scale research has been undertaken on the links between the litter layer, runoff, and ET.

The objective of this study was to quantify the effects of the litter layer on runoff and ET using the paired watershed method (Bosch and Hewlett 1982; Brown et al. 2005; Dung et al. 2012; Komatsu et al. 2012). To achieve this, we analyzed raw data obtained from an experiment in which the litter layer was removed annually from an entire watershed, from 1954 through 1956. The novelty of this research is to quantify the changes in stream runoff and ET before and after mechanical removal of the litter layer on the watershed scale.
Materials and Methods

### Study sites

The paired study sites documented herein comprised two small watersheds: North Creek (NC; watershed area 1.19 ha) and South Creek (SC; 1.42 ha), which are located in the Shirasaka Experimental Watershed (35°13′07′′N, 137°09′54′′E) in the Akazu Research Forest of the Ecohydrology Research Institute (ERI), The University of Tokyo Forests (Fig. 1; Yamaguchi 1963; ERI 2013). Meteorological data collected in an open field located approximately 240 m west–northwest of the outlet of the NC watershed indicate that the annual average temperature between 1935 and 2014 was 12.4 °C, and the annual average rainfall between 1930 and 2014 was 1861 mm/year (Gomyo and Kuraji 2013; Kuraji and Gomyo 2014; data missing in 2005). The bedrock in both watersheds is deeply weathered Cretaceous granite. Their elevations are 320–348 m above sea level and they have slope gradients in the range 3.3–40.8°.
In the 1950s, both watersheds’ vegetation comprised mixed forests of Japanese red pine (Pinus densiflora) and deciduous broadleaved trees, mainly konara oak (Quercus serrata). The highest Japanese red pine and konara oak trees were 28 and 16 m, respectively. In 1954, all trees were measured within 4-m-wide belt transects across both watersheds. Within these areas, the basal area (BA) of Japanese red pine and deciduous broadleaved trees accounted for 54 and 40 % of the BA of all trees in the NC watershed, and 43 and 45 % of those in the SC watershed, respectively (Kataoka et al. 1954).
Litter removal experiments were conducted on three occasions: January 6–15, 1954; January 6–15, 1955; and January 9–13, 1956. All litter from the NC watershed was removed and weighed onsite, whereas the SC watershed was preserved as a control. The weights of litter removed each year from NC are shown in Table 1. The raw data (field notes) produced during these experiments have since been maintained in good condition by the ERI.
Table 1
Date, number of days and total weight of litter removed from the North Creek watershed
Year
Date
Number of days
Litter weight (× 104 kg)
Litter weight (kg m−2)
1954
7–12, 14–15 January*
8
1.800
1.5
1955
10–14 January
5
1.170
1.0
1956
9–13 January
5
1.465
1.2
The litter weight removed in 1954 was the litterfall before December 1953 (long-term) minus litter decomposed during the same period
Litter weights removed in January 1955 and 1956 were the litterfall from January 1954 to December 1955 (1 year) minus litter decomposed during the same period and the litterfall from January 1955 to December 1956 (1 year) minus litter decomposed during the same period, respectively
* No work on 13 January 1954 due to rainy weather

### Hydrological observations

Precipitation in the field was measured at the above-mentioned meteorological observation station using a storage-type rain gauge (20-cm diameter) and a siphon rain gauge (20-cm diameter). The storage-type rain gauge was examined every day at 09:00 local time (LT) where a cylinder was used to measure the volume of water collected in a storage bottle within the gauge. This result was considered the rainfall amount for each day beginning at 09:00 LT. At the siphon rain gauge, recording paper that covered the events of a single day was read at 09:00 LT to determine whether rainfall had occurred during the previous 24 h; if any rainfall had been recorded, the paper was replaced. Hourly rainfall data used in this work were read from the original recording papers, which were maintained in good condition by the ERI, and then corrected such that the 24-h rainfall matched the value obtained from the storage-type rain gauge.
Runoff was observed using weirs installed at the outlets of the watersheds. Each weir consisted of a stilling pool (6.0 × 9.7 m in the NC and 6.0 × 9.0 m in the SC), a 60° V-notch weir, an automatic recording water-level gauge, and a point gauge. The recording paper from the water-level gauge covered a 1-day period and was replaced daily at 09:00 LT. The values from the point gauge were read and recorded on the original recording paper, which were maintained in good condition by the ERI and subsequently examined in this work. Any changes in water level recorded on the recording paper were noted at 1-min intervals. Assuming a linear change in the water level between adjacent points over time, the equation used to calculate flow rate based on the water level was integrated over time to calculate the runoff during this period. These values were then summed each day to obtain the daily runoff. The following formula was used to calculate the runoff based on the water levels at the NC and SC watersheds:

$Q\phantom{\rule{1pt}{0ex}}\phantom{\rule{1pt}{0ex}}\alpha \phantom{\rule{1pt}{0ex}}H{\phantom{\rule{1pt}{0ex}}}^{2.5}$
(1)
where Q is runoff [m3/s], H is water level [m], and α is the coefficient of runoff determined by the observation of Q and H (NC watershed: 0.7669, SC watershed: 0.8086).
In this study, the pre- and post-treatment periods were 1951–1953 and 1954–1956, respectively. There were no missing daily rainfall measurements, but daily runoff measurements were missing in the NC watershed for 2 days in 1951, 3 days in 1953, 4 days in 1954, 9 days in 1955, and 2 days in 1956, as well as for 3 days in the SC watershed in 1955. To address these data gaps in each watershed, we multiplied the available daily runoff measurement in the other watershed on the date of the missing data by the ratio of the daily runoff in the NC and SC watersheds on the day before or after that of the missing data. In 1955, there was 1 day on which measurements were missing for both watersheds, although no precipitation occurred; thus, the missing measurements were complemented with the average of the daily runoffs recorded on the previous and following days.

### Paired watershed method

To evaluate the effect of litter removal on the annual runoff in the NC watershed, we estimated the annual runoff during the post-treatment period under the assumption that treatment had not occurred (QENA). A linear regression equation between the annual observed runoff in the NC and SC watersheds during the pre-treatment period was used:

${Q}_{\text{NB}}a{Q}_{\text{SB}}b$
(2)
where QNB and QSB are the annual observed runoffs in the NC and SC watersheds during the pre-treatment period, respectively, and a and b are regression coefficients. We then estimated QENA by the following equation:

${Q}_{\text{ENA}}a{Q}_{\text{SA}}b$
(3)
where QSA is the annual observed runoff in the SC watershed during the post-treatment period, and a and b are the regression coefficients determined by Eq. (2).

### Magnitude of peak flow

We detected how the magnitude of peak flow during storm events changed between the pre- and post-treatment periods by tabulating data from all storm events with a peak-flow magnitude at the SC control watershed (QPS) greater than 0.1 mm/h (=0.028 L/s/ha; Harr and McCorison 1979). If one event with the lag time between the timing of the peak flow in the NC and SC was <60 min, the event was rejected because the rainfall peaks that corresponded to peak flows in the SC and NC may have been different. In total, 41 and 66 peak flow events that met this definition were identified during the pre- and post-treatment periods, respectively. The ratio of the magnitude of peak flow in the NC watershed (QPN) and QPS (hereafter, referred to as “peak flow ratio”) was compared between the pre- and post-treatment periods. To identify the statistical significance of the difference between the pre- and post-treatment periods, we applied an analysis of covariance (ANCOVA). Peak flows were log-transformed to meet the assumption of homoscedasticity and to improve the frequency distribution of the data along the interval of the regression (Wright et al. 1990). The significances of the regression relationship and parallelism were tested before testing for different intercepts.

Results

### Observed hydrograph

The 6-year hyeto-hydrographs observed at daily intervals in the NC and SC watersheds are shown in Fig. 2a, b, respectively. The runoff from both watersheds was similar; however, the observed peak daily runoff in the NC watershed (QON) tended to be greater than that in the SC watershed (QOS) during the pre- and post- treatment periods. Because their rainfall–runoff responses were similar (the correlation coefficients between QON and QOS during the pre- and post-treatment periods were 0.986 and 0.980, respectively), the paired watershed method was considered appropriate. The 6-year cumulative QON and QOS curves are shown in Fig. 2c. There was no clear difference between the cumulative QON and QOS during the pre-treatment period, whereas the cumulative QON was greater than the cumulative QOS during the post-treatment period.

### Annual water balance

Annual precipitation, annual runoff in the NC and SC watersheds, and the difference between the annual runoff in the NC and SC watersheds are shown in Table 2. The pre-treatment 3-year QON was 46.1 mm greater than the 3-year QOS, whereas the post-treatment 3-year QON was 163.2 mm greater than the 3-year QOS.
Table 2
Annual and 3-year precipitation, runoff, and difference between the North Creek (NC) and South Creek (SC) watersheds during the pre- and post-treatment periods, and estimated annual and 3-year runoff, and differences between the estimated and observed annual runoffs in the NC watershed during the post-treatment period
Year
Precipitation (mm)
Discharge (mm)
Difference (mm)
Estimated discharge (mm)
Difference (mm)
Pre-treatment period
QSB
QNB
QNBQSB
1951
1906.4
945.1
972.8
27.7
1952
2176.4
1249.8
1262
12.2
1953
2123.4
1281.9
1288.1
6.2
1951–53
6206.2
3476.7
3522.8
46.1
Post-treatment period
QSA
QNA
QNAQSA
QENA
QENAQNA
1954
2063.8
1044.4
1124.2
79.8
1066.5
57.7
1955
1690.7
734.1
790.7
56.6
774.4
16.3
1956
1941.7
1070.5
1097.3
26.8
1091.0
6.3
1954–56
5696.2
2849.0
3013.1
163.2
2931.9
80.3
From Eqs. (2) and (3), the annual runoff in the NC watershed during the post-treatment period that would be expected in the absence of treatment (QENA) was obtained by the following equation:

${Q}_{\text{ENA}}\phantom{\rule{thinmathspace}{0ex}}0.94\phantom{\rule{thinmathspace}{0ex}}{Q}_{\text{SA}}\phantom{\rule{thinmathspace}{0ex}}83.3\phantom{\rule{thinmathspace}{0ex}}{R}^{2}\phantom{\rule{thinmathspace}{0ex}}0.9999$
(4)
where QSA is the annual runoff in the SC watershed during the post-treatment period (mm/year).
The QENA for 1954, 1955, and 1956 (Table 2) was 57.7, 16.3, and 6.3 mm less than the annual runoff in the NC watershed during the post-treatment period (QNA) for the same years, respectively, indicating that litter removal increased the annual runoff in the NC watershed. The 3-year total runoff increased by 80.3 mm during the post-treatment period.

### Magnitude of peak flow

Figure 3a shows the relationship between the QPN and QPS during the pre- and post-treatment periods. Although there is scattering, the QPN during the post-treatment period was larger than during the pre-treatment period. Figure 3b shows the average peak flow ratio obtained for four QPS ranges during the pre- and post-treatment periods. The results of ANCOVA are shown in Table 3. No significant difference in the peak flow ratio was detectable between the pre- and post-treatment periods when the QPS range was 0.1–0.25 or 0.25–0.4 mm/h; however, the peak flow ratio of the post-treatment period was greater than that in the pre-treatment period when the QPS range was 0.4–1.0 or >1.0 mm/h. For the QPS range of 0.4–1.0 mm/h, the post-treatment peak flow ratio was 1.46, which is about 1.5 times that of the pre-treatment period (0.98). In addition, the post-treatment peak flow ratio was 1.81 during a flood event when the QPS range was >1.0 mm/h, which is about 1.4 times that of the pre-treatment period (1.29). The treatment effect on the QPN was statistically significant when the QPS range was >0.4 mm/h (ANCOVA, Table 3).

Table 3
Result of ANCOVA for the treatment effect
Range of QPS(mm/h)
Pre-treatment, N
Post-treatment, n
Regression p
Parallelism p
Intercept p
0.1–0.25
8
9
0.009*
0.757
0.764
0.25–0.4
7
20
0.005*
0.364
0.643
0.4–1.0
13
19
p < 0.001*
0.185
p < 0.001*
>1.0
13
18
p < 0.001*
0.669
p < 0.001*
Note that the regression p indicates there is a significant regression relationship with the covariance and the parallelism p indicates the regression lines for these individual groups are assumed to be non-parallel; in other words, they have the different slope. The intercept p indicates the intercept of the regression line for these individual groups are assumed to be different
* Significant with α = 0.01

Discussion

### Effect of litter removal on ET

This paired watershed experiment showed that litter removal increased the 3-year total runoff in the NC watershed by 80.3 mm. Thus, the annual loss of water decreased, as calculated by subtracting annual runoff from annual precipitation. The difference between the annual loss of water and annual ET is the difference in watershed water storage between the start and the end dates of the study period, alongside deep percolation that cannot be measured by a weir (Oda et al. 2008). We estimate that litter removal reduced the annual ET because the difference in watershed water storage might be smaller compared with the 3-year precipitation, runoff, and ET. Although about 5 % of the annual precipitation could possibly have been deep percolation in this study site (Terajima et al. 1993), we suggest that the impact of litter removal on deep percolation might be negligible compared with the impact on runoff and ET.
The observed decrease in annual ET caused by litter removal implies that the presence of a litter layer enhances ET from a forested watershed. One possible mechanism to explain this result could be that the increase in litter interception exceeded the reduction in soil surface evaporation.
Previous studies have suggested that differences in litter weight, rainfall intensity and frequency, and evaporative demand of the forest floor could cause large variations in the volume of forest floor interception and the proportion of throughfall (Gerrits and Savenije 2011; Zagyvainé Kiss et al. 2014). Helvey (1967) reported annual litter interception amounts of 30, 46, and 56 mm in 10-, 30-, and 60-year-old stands of eastern white pine (Pinus strobus), respectively, in which the annual rainfall was about 1500 mm and the litter weight was 1.2 kg/m (similar to our study: 1.0–1.5 kg/m). The forest type in Helvey’s study was similar to that of ours, and the estimated decrease in 3-year ET in our study (80.3 mm) is comparable with that reported by Helvey (1967), which suggests that forest floor evaporation could be accounted for by the decrease in litter interception after litter removal from the entire watershed. Bulcock and Jewitt (2012) reported 3-year litter interceptions of 160, 125, and 231 mm in EucalyptusAcacia, and Pinus stands, respectively. These values are also comparable with those of our study. The litter weights reported by Bulcock and Jewitt (2012; 2.3, 2.4, and 3.3 kg/m for EucalyptusAcacia, and Pinus stands, respectively) were about two or three times greater than our study (1.0–1.5 kg/m). The comparable value of litter interception between our study and Bulcock and Jewitt (2012) may be the difference of rainfall condition. The 3-year rainfall for the study period of Bulcock and Jewitt (2012) was 1885–1910 mm, which is about one-third of that in our study period (5696.2 mm/three-year during the post-treatment period). Thus, the decrease in ET identified herein is comparable with that of previous stand-scale studies; however, because of the rarity of such work, we cannot conclude that the observed forest floor evaporation is consistent at all observational scales (i.e., watershed- or stand-scale studies).

### Peak flow

Figure 3 shows that the peak flow ratio was greater in the post-treatment period than in the pre-treatment period, with the degree varying according to the scale of the flood event.
The effects of a litter layer on runoff in a forested watershed are the occurrence of a rapid-flow component within the litter layer, litter interception, mitigation of the impact of raindrops, suppression of surface flow, and an increase in permeability (Walsh and Voigt 1977). The rapid-flow component within the litter layer increases the peak flow, whereas the other factors might act to reduce it.
The QPN during the post-treatment period had relatively little impact for relatively small events (QPS range of <0.4 mm/h), and the presence or absence of litter had almost no effect on the peak runoff quantity during flood events within this range. By contrast, the removal of litter increased the peak flow by 1.4–1.5 times for relatively greater flood events (QPS range of >0.4 mm/h). For flow events with a QPS range of <0.4 mm/h, either the increased and decreased effects of the litter layer on the peak flow were offset, or both effects could have been insignificant. However, for flood events with a QPS range of >0.4 mm/h, the effects on peak flow reduction by the litter layer (litter interception, raindrop impact mitigation, overland flow suppression, permeability increase; Walsh and Voigt 1977) were remarkably stronger than its effect on the QPN increase. Thus, at least in this study, the loss of these effects in relation to litter removal probably induced the peak flow increase.

Conclusions

A paired watershed method was applied in this study to identify how a litter layer controls runoff and evapotranspiration. Approximately 4.4 × 104 kg of litter was removed from the forest floor of an entire watershed for three consecutive years, which increased the 3-year runoff by a total of 80.3 mm. Results showed that the peak runoff during a large flood event was ~1.4–1.5 times greater due to removal of the litter layer; indeed, increased runoff due to litter layer removal occurred mainly during such flood events. If a litter layer exists in a plantation forest, its removal may have a negative impact on flood mitigation, but a positive impact on water resource management. Conversely, in a plantation forest from which a litter layer has been lost, its recovery may have a positive impact on flood mitigation, but a negative impact on water resource management. Further observational and modeling studies are necessary to clarify the role of a litter layer on the processes and mechanisms of runoff during flood events, and to evaluate the impact of the increase or decrease of litter volume on flood mitigation measures and water resources management.

Acknowledgements

We thank Yoshie Gomyo (project specialist) and Junko Kato (project academic support staff) for their support of this study. The authors express their appreciation to the Ecohydrology Research Institute, The University of Tokyo Forests, and the Graduate School of Agricultural and Life Sciences at The University of Tokyo for their continuous support of our activities. We acknowledge the two anonymous reviewers for providing helpful comments on an initial version of this manuscript. We would like to thank Editage (www.editage.jp) for English language editing.

### Funding

This study was funded by the Japan Society for the Promotion of Science KAKENHI (Grant Number 25·6722).

### Conflict of interest

The authors declare that they have no conflict of interest.

References

1. Bosch JM, Hewlett JD (1982) A review of catchment experiments to determine the effect of vegetation changes on water yield and evapotranspiration. J Hydrol 55:3–23
2. Brown AE, Zhang L, McMahon TA, Western AW, Vertessy RA (2005) A review of paired catchment studies for determining changes in water yield resulting from alterations in vegetation. J Hydrol 310:28–61
3. Bulcock HH, Jewitt GPW (2012) Field data collection and analysis of canopy and litter interception in commercial forest plantations in the KwaZulu-Natal Midlands, South Africa. Hydrol Earth Syst Sci 16:3717–3728
4. Dung BX, Gomi T, Miyata S, Sidle RC, Kosugi K, Onda Y (2012) Runoff responses to forest thinning at plot and catchment scales in a headwater catchment draining Japanese cypress forest. J Hydrol 444–445:51–62
5. ERI (Ecohydrology Research Institute) (2013) Several reports in 1955 on water-level measurements at four gauging stations of Ananomiya, Higashiyama, Kazunari, and Shirasaka. Bull Univ Tokyo For 53:129–152 (in Japanese with English abstract)Google Scholar
6. Gallart F, Llorens P (2003) Observations on land cover changes and water resources in the headwaters of the Ebro catchment, Iberian Peninsula. Phys Chem Earth 29:769–773
7. Gerrits AMJ, Savenije HHG (2011) Forest floor interception: forest hydrology and biogeochemistry: synthesis of past research and future directions. Ecol Stud 216:445–454
8. Gerrits AMJ, Savenije HHG, Hoffmann L, Pfister L (2007) New technique to measure forest floor interception––an application in a beech forest in Luxembourg. Hydrol Earth Syst Sci 11:695–701
9. Gerrits AMJ, Pfister L, Savenije HHG (2010) Spatial and temporal variability of canopy and forest floor interception in a beech forest. Hydrol Process 24:3011–3025
10. Gomyo M, Kuraji K (2013) Long-term variation of annual loss and annual evapotranspiration in small watershed with the forest restoration and succession on denuded hills. J Jpn For Soc 95:109–116 (in Japanese with English abstract)
11. Harr RD, McCorison FM (1979) Initial effects of clearcut logging on size and timing of peak flows in a small watershed in Western Oregon. Water Resour Res 15:90–94
12. Hattori S (1993) Forest and evapotranspiration. In: Japan Society of Erosion Control Engineering (ed) Soil erosion, water discharge and forest influence. Sankaido Co. Ltd, Tokyo, pp 171–198 (in Japanese)Google Scholar
13. Helvey JD (1967) Interception by eastern white pine. Water Resour Res 3:723–729
14. Kataoka J, Takasugi K, Akiya K (1954) unpublished manuscript (in Japanese)
15. Komatsu H, Jaeil C, Matsumoto K, Otsuki K (2012) Simple modelling of the global variation in annual forest evapotranspiration. J Hydrol 420–421:380–390
16. Komatsu H, Shinohara Y, Otsuki K (2015) Models to predict changes in annual runoff with thinning and clearcutting of Japanese cedar and cypress plantations in Japan. Hydrol Process 29:5120–5134. doi:10.1002/hyp.10520
17. Kuraji K, Gomyo M (2014) Estimation of the variation of 78 years’ mean annual temperature in the Ecohydrology Research Institute of The University of Tokyo Forests, and its relation to the temperature increase in Nagoya caused by urbanization with a correction for instrument error. Bull Univ Tokyo For 131:29–40 (in Japanese with English abstract)Google Scholar
18. Lee HH, Shibano H (1990) Hydraulic characteristics of storm water flowing through litter layer. J Jpn For Soc 72:223–229 (in Japanese)Google Scholar
19. Li X, Niu J, Xie B (2013) Study on hydrological functions of litter layers in north China. PLoS One 8:1–11Google Scholar
20. Majima M, Kitahara H, Shimizu A (1990) Measurement of water balance at different ground cover with plain lysimeter. Proc Boreal For Soc 38:233–235 (in Japanese)Google Scholar
21. Molina AJ, del Campo AD (2012) The effects of experimental thinning on throughfall and stemflow: a contribution towards hydrology-oriented silviculture in Aleppo pine plantations. For Ecol Manag 269:206–213
22. Murai H (1993) Forest soil––function of water and soil conservation. In: Japan Society of Erosion Control Engineering (ed) Soil erosion, water discharge and forest influence. Sankaido Co. Ltd, Tokyo, pp 199–218 (in Japanese)Google Scholar
23. Oda T, Asano Y, Suzuki M (2008) Estimating deep percolation in a small catchment in a tertiary formation using the chloride mass balance method. J Jpn Soc Hydrol Water Resour 21:195–204
24. Okunishi K (1963) Rainfall and runoff on mountainside––Part 1. Runoff from fallen-leaf layer. Disaster Prev Res Inst Annu 6:156–165 (in Japanese with English abstract)Google Scholar
25. Onda Y, Gomi T, Mizugaki S, Nonoda T, Sidle RC (2010) An overview of the field and modelling studies on the effects of forest devastation on flooding and environmental issues. Hydrol Process 24:527–534
26. Rao L, Zhu J (2007) Hydrological effects of forest litter and soil in the Simianshan Mountains in Chongqing, China. Front For China 2:157–162
27. Sato Y (2007) Rainfall interception––How to estimate the amount of water reaching to the forest soil. Forest Hydrology––Exploring the fate of water in forest ecosystems. The editing committee of Forest Hydrology, Morikita Publishing Co. Ltd, Tokyo, pp 14–39 (in Japanese)Google Scholar
28. Sato Y, Kumagai T, Kume A, Otuki K, Ogawa S (2004) Experimental analysis of moisture dynamics of litter layers: the effects of rainfall conditions and leaf shapes. Hydrol Process 18:3007–3018
29. Schaap MG, Bouten W (1997) Forest floor evaporation in a dense Douglas fir stand. J Hydrol 193:97–113
30. Terajima T, Mori A, Ishii H (1993) Comparative study of deep percolation amount on two small catchments in granitic mountain. Hydrology (J Jap Assoc Hydrol Sci) 23:105–118 (in Japanese with English abstract)Google Scholar
31. Walsh RPD, Voigt PJ (1977) Vegetation litter: an underestimated variable in hydrology and geomorphology. J Biogeogr 4:253–274
32. Wright KA, Sendek KH, Rice RM, Thomas RB (1990) Logging effects on streamflow: storm runoff at Casper Creek in northeastern California. Water Resour Res 26:1657–1667Google Scholar
33. Yamaguchi I (1963) Fundamental studies on the groundwater flow phenomena at the water source zone. Bull Univ Tokyo For 58:133–285 (in Japanese)Google Scholar
34. Zagyvainé Kiss KA, Kalicz P, Csáfordi P, Gribovszki Z (2014) Forest litter interception model for a sessile oak forest. Acta Silv Lig Hung 10:91–101Google Scholar
35. Zhang Z, Lei Y, Su K, Wang G, Wang D, Ma H (2009) Hydrological characteristics of litter in different forest succession stages at Liuxihe Watershed, South China. Front For China 4:317–322

For further details log on website :