Next Article in Journal
A Tutorial on Model-Assisted Estimation with Application to Forest Inventory
Next Article in Special Issue
Leaf Phenology Drives Spatio-Temporal Patterns of Throughfall under a Single Quercus castaneifolia C.A.Mey.
Previous Article in Journal
Growth-Ring Analysis of Diploknema butyracea Is a Potential Tool for Revealing Indigenous Land Use History in the Lower Himalayan Foothills of Nepal
Previous Article in Special Issue
Influence of Root Distribution on Preferential Flow in Deciduous and Coniferous Forest Soils
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rainfall Partitioning in Chinese Pine (Pinus tabuliformis Carr.) Stands at Three Different Ages

College of Forestry, Beijing Forestry University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Forests 2020, 11(2), 243; https://doi.org/10.3390/f11020243
Submission received: 17 January 2020 / Revised: 12 February 2020 / Accepted: 19 February 2020 / Published: 22 February 2020
(This article belongs to the Special Issue Water Resources Management and Modeling in Forestry)

Abstract

:
Chinese pine (Pinus tabuliformis Carr.) is the main forest species in northern China, with the potential to dramatically affect biotic and abiotic aspects of ecosystems in this region. To discover the rainfall partitioning patterns of different growth periods of Chinese pine forest, we studied the throughfall (Tf), stemflow (Sf) and canopy interception (I) in three stand ages (40-, 50-, 60-year-old) in Liaoheyuan Natural Reserve of Hebei Province during the growing seasons of 2013 and 2014, and analyzed effect of rainfall amount, rainfall intensity, and canopy structure on rainfall partitioning in Chinese pine forest. The results showed that throughfall decreased with the stand age, accounting for 78.8%, 74.1% and 66.7% of gross rainfall in 40-, 50- and 60-year-old Chinese pine forests, respectively. Canopy interception, on the other hand, increased with the stand age (20.4%, 24.8%, and 32.8%, respectively), while the pattern in stemflow was less clear (0.8%, 1.1%, and 0.6%, respectively). As rainfall intensity increased, the Tf and Sf increased and I declined. Additionally, our results showed that leaf area index (LAI) and the diameter at breast height (DBH) increased with age in Chinese pine stands, probably explaining the similar increase in canopy interception (I). On the other hand, the mean leaf angle, openness, gap fraction all decreased with the stand age. Stepwise regression analysis showed that the rainfall amount and LAI were the major determinants influencing the rainfall partition. Our study highlights the importance of stand age in shaping different forest canopy structures, and shows how age-related factors influence canopy rainfall partitioning. This study also significantly adds to our understanding the mechanisms of the hydrological cycle in coniferous forest ecosystems in northern China.

1. Introduction

Forests can influence hydrological processes and alter soil conditions due to their many impacts on water purification, runoff regulation and water and soil conservation [1]. In this capacity, forests can also serve to improve water availability at watershed, regional and worldwide levels [2]. The forest canopy plays an important role in this process, by partitioning rainfall and affecting additional hydrological processes [3], not only changing the spatial distribution of rainfall in the forest [4,5], but also affecting the chemical composition of the falling rain [6,7]. Previous studies of forest stands indicated that forest rainfall interception led to heterogeneous water distribution [8,9,10] and showed a large range of spatial variation relating to different forest structures [11].
Rainfall is redistributed through a forest canopy by throughfall, stemflow and canopy interception. Throughfall is the total volume of raindrops that penetrate through the forest canopy and reach to the ground. Stemflow is the part of rainfall flowing down from stems or trunks, and the amount intercept by the vegetation canopy and evaporate is the interception loss [12]. Throughfall accounts for the highest proportion of rainfall, and it is the main source influencing the spatial distribution of water within the soil [13]. The proportion of stemflow is small, but it can be an important source of ground water [14] and benefit to forest regeneration. Canopy interception returns to the atmosphere through evaporation [15]. Given their importance in ecosystem function and hydrology, many studies about the rainfall partitioning for several vegetation types in different places have conducted [16,17,18].
Rainfall redistribution in forests is a complex process, directly influenced by rainfall characteristics such as amount, intensity and duration [19], as well as environmental factors including wind speed, air humidity and air temperature. Studies have shown that rain events with high rainfall volume are known to induce higher throughfall [20]. In addition, the rainfall amount and the rainfall intensity are positively correlated with the canopy interception [5]. Interannual and intraseasonal variability of rainfall significantly influences the interception loss. Raindrop size and velocity diameter have been shown to influence the throughfall of the maize canopy [21]. Fog in a forest can also increase water input below the canopy [1,15].
Multiple plant characteristics are also important factors in rainfall redistribution. Shrubs present higher stemflow than trees [17]. Branching angle and canopy volume also influence the rainfall distribution in Turkish oak forests, where slender trees with upward thrusting branches producing more stemflow [6]. Stand density is an important factor for all rainfall partitioning elements in coniferous forests [22]. Campo et al. found stemflow ratios decreased with tree density [23]. Tree size, branching angles and bark roughness are the main factors controlling stemflow volumes [24,25]. Spatial variation of canopy interception has a significant relationship with the plant area index, but not with the leaf area index [11]. Other factors including the gap fraction, horizontal and vertical distribution of foliage, present of epiphytes and species composition also influence the interception [26].
Little is known, however, about the relationship between rainfall redistribution patterns and forest stand age in temperate forests especially in a Chinese pine forest, which could represent an important interaction in forest ecosystems. There were some reports on different succession in various places and different forest types. For example, a previous study has shown that stand age, as an indicator of forest successional status, plays a vital role in canopy interception loss in a range of age class of deciduous forest in the southern Appalachian Mountains [27]. In tropical dry forests net rainfall showed an obviously declined trend across successional stages in Costa Rica, but in Brazil no clear trend was found [28]. The annual interception ratio showed dominantly larger in mature stands in Chinese fir stands [29]. Many experiments have been conducted on rainfall partition worldwide, but few consider both stand age and rainfall redistribution patterns.
Liaoheyuan Natural Reserve is the origin of the Liaohe River. Within the reserve, there are large areas of Chinese pine natural forests, which play an important role in water and soil conservation for the Liaohe River Basin. Liaoheyuan Natural Reserve is also an important ecological barrier of wind and sand fixation for Beijing and Tianjin. In this study, we targeted three Chinese pine stand ages in Liaoheyuan Natural Reserve and monitored the precipitation, throughfall and stemflow in two consecutive growing seasons. The objectives of this study were: (1) to study the dynamic rainfall redistribution patterns for different stand ages of Chinese pine forests, and (2) to evaluate the rainfall characteristics and forest structures that affected rainfall partitioning in these pine forests. This study provided basic data for hydrological effects of forest vegetation in Liaoheyuan and a theoretical basis for the construction, management, and protection of water source conservation in Liaoheyuan.

2. Materials and Methods

2.1. Study Site

This study site was located in Liaoheyuan Nature Reserve, Pingquan County, Hebei Province, China (Figure 1) (latitude 41°01′–41°21′ N, longitude 118°22′–118°37′ E, elevation 625–1738 m a.s.l.). This area is a transition zone from warm to cold temperate zones with a typical continental monsoon climate, and has a humid and warm summer and snowy and cold winter. The average temperature of a year is 7.7 °C and the mean precipitation of a year is 516.9 mm (1980–2018), according to data from Pingquan County Meteorological Bureau, with most rainfalls occurring from May to August. The soil is a typical brown forest soil, 0–100 cm thick [30]. The dominant tree species is the evergreen conifer P. tabuliformis, or Chinese pine, which is widely distributed in north China. Spiraea pubescens and Lespedeza bicolor are the main shrub species. Carex rigescens, Saussurea nivea, and Dianthus chinensis are major herbaceous plant species.
In July 2012, we established three 20 m × 30 m study plots, each of which consisted of Chinese pine forest of a different age, namely the 40-year-old stand (latitude 41°18.522 N, longitude 118°31.169 E, elevation 1004 m a.s.l. (Above sea level)), the 50-year-old stand (latitude 41°18.572 N, longitude 118°31.188 E, elevation 998 m a.s.l.), and the 60-year-old stand (latitude 41°18.585 N, longitude 118°31.301 E, elevation 997 m a.s.l.). We got the plot basic information on vegetation and site characteristic data by field investigation. Species name, height, crown breadth, diameter and density of trees with DBH ≥ 5 cm were recorded for each plot. Nine digital hemispherical photographs (WinScanopy 2010a) were taken for each plot to get LAI and other structure information on overcast days or near sunrise or sunset. The camera was adjusted to horizontal position at 1.3 m above ground [31]. Photographs were analyzed using WinScanopy software for calculating LAI, the gap fraction, the openness, and the mean leaf angle [26]. Our study was conducted based on the National Standards of the People’s Republic of China or “Methodology for field long-term observation of forest ecosystem research” (GB/T 33027-2016). Data for each site are shown in Table 1.

2.2. Collection of P, Tf, and Sf

Precipitation (P), throughfall (Tf) and stemflow (Sf) were collected from May to September of 2013 and 2014. Precipitation was recorded by an automated tipping-bucket gauge (resolution 0.2 mm) at a meteorological station (CR1000, Compbell Scientific Inc., Logan, UT, USA), installed in an open hillside 1000 m away from the plots (Figure 2a). A single rainfall event was classified as rain exceeding 0.5 mm during a period with more than 6 h from the last rainfall event [32]. Rainfall events were divided into 4 intensity classes, including light rain (0.2 mm/d ≤ P < 10 mm), middle rain (10 mm/d ≤ P < 25 mm), heavy rain (25 mm/d ≤ P < 50 mm), rain storm (50 mm/d ≤ P < 100 mm), according to Criteria for Classification of Rainfall Intensity (Inland) produced by China’s Meteorological Administration.
Throughfall (Tf) was monitored manually by using 3 homemade plastic rain gauges (200 × 10 cm) randomly placed in each plot, which were set on a bracket 50 cm above the ground (Figure 2b). Each gauge kept an angle of 5° above ground to make drainage easily, the end of the gauge connected to a plastic tank (16 L), which was in a hole below the ground covered with plastic sheeting to avoid evaporation and splash. The volume of water in the tank was measured soon after each rainfall event by using a graduated cylinder with the resolution of 1 mL. The throughfall depth was calculated by the rain volume divided by the area of the gauge. The average of the three throughfall depths was adopted as the throughfall amount for each plot.
The stemflow (Sf) was collected by using a 5 mm thick white halved rubber hose, which was spirally wrapped around tree trunk and used Vaseline to create a water-tight environment between the stem and the hose. The lower side of the hose was 30 cm above the ground and connected to a plastic bucket (16 L), after each rainfall event the rain volume in the bucket was measured (Figure 2c). Trees for stemflow experiments were selected according to the diameter at breast height (DBH) within the plot. In each plot all the trees were divided into five groups by DBH. DBH groups for three stands were divided as follows: the 40-year-old stand (11 ≤ d < 15, 15 ≤ d < 19, 19 ≤ d < 23, 23 ≤ d < 27, 27 ≤ d < 31), the 50-year-old stand (5 ≤ d < 10, 10 ≤ d < 16, 16 ≤ d < 22, 22 ≤ d < 28, 28 ≤ d < 34), the 60-year-old stand (20 ≤ d < 25, 25 ≤ d < 30, 30 ≤ d < 35, 35 ≤ d < 40, 40 ≤ d < 45). Then one tree was chosen in each group as the standard tree to monitor stemflow for each plot in the studied forest. A total of five trees in each plot were selected as standard trees, information of which was shown in Table 2.
The depth of the stemflow is calculated by the following equation [14].
S f = i = 1 n C i × M i S × 1000
where Sf is the stemflow amount (mm), n is the number of stem groups in this study; Ci is the stemflow amount of monitored tree in ith group (mL); Mi is the number of trees in the ith group (n); S is the plot area (600 m2).
If total precipitation (P), throughfall (Tf), and stemflow (Sf) are known, canopy interception amount (I) can be measured according to the following equation [33]:
I = P T f S f
where I is interception amount (mm), P is the precipitation (mm), Tf is the throughfall amount (mm), Sf is the stemflow amount (mm).
Canopy interception ratio (R) is then calculated as:
R = I P 100 %
where R is interception ratio (%), P is the precipitation (mm), I is the interception amount (mm).

2.3. Data Analysis

Throughfall depth under each rainfall event for each stand was average amount of the three gauges. We used cumulative data during the study period to reflect the overall rainfall partitioning pattern, and used mean values to distinguish the difference between stand ages and months.
All data preparations including the total precipitation, amount and ratio data were done by Microsoft Excel (2010 version, Microsoft Corporation, Redmond, Washington, USA). A one-way analysis of variance (ANOVA) by SPSS (IBM SPSS Statistics 19 version, and ratio data were done by Microsoft Excel (2010 version) was used to compare the statistical difference between different Chinese pine stand ages. Stepwise multiple linear regression was conducted with SPSS to find out significant factors for Chinese pine rainfall partitioning.

3. Results

3.1. Rainfall Pattern During the Study Period

During May-September of 2013 and 2014, a total of 64 rainfall events produced 857.4 mm of precipitation (P), each event ranging from 0.6 to 58.8 mm of rain with an average of 13.4 mm. The number of rainfall events and amount of precipitation varied across months and years, with 30 total rainfall events and 380.2 mm of rain in 2013, and 34 events and 477.2 mm of rain in 2014. The maximum rainfall amount for a single rainfall event was 58.8 mm, which occurred on 28 June 2013. In 2014, the maximum amount of rain from a single event was 42.6 mm on 5 June.
According to the classification of the “precipitation intensity classification standard (inland part)” of China’s Meteorological Administration, the precipitation events in the study area during the observation period were most often light rain, with 40 times total and rainfall amounting to 194.4 mm, accounting for 22.7% of total rainfall in growing seasons of 2013 and 2014 (Table 3). Increasingly intense rainfall events were less frequent. Rain storms occurred just twice, both in 2013. The gross rainfall amount of heavy rain was the largest, followed by middle rain, light rain, and rain storm.

3.2. Rainfall Partitioning Pattern Across Chinese Pine Stand Age

During the two years of the observation period, the accumulated rainfall division in (mm) and (%) showed different results among the three forest stands (Table 4). In each plot, Tf corresponded to the highest proportion of the total rainfall, while Sf was always the lowest proportion. Tf showed an expected decrease from 40 to 60 years. Stemflow (Sf), on the other hand, did not show an obvious trend, since the 50-year-old stand resulted in the highest value among the three stands. Total interception (I), on the other hand, behaved similarly to Tf. An additional analysis, based on 64 sampled rainfall events (mm), showed I was the only variable that proved to have a significant statistical difference among the three forest stands, being the 60 years plot significantly greater than that of 40 and 50 years.
Monthly rainfall distribution patterns for two years (Figure 3) showed that throughfall, stemflow and interception amount of three stand ages in June and August were larger than the other three months. The Tf, Sf, and I of three stands presented the same trend across five months and the trend among three stands varied. Throughfall in each month decreased as the stand age increases, whereas interception increased with the stand age. The stemflow was largest in the 50-year-old stand in each month.

3.3. Dependence of Rainfall Partitioning on Rainfall Amount

The relationships between rainfall amount and throughfall, stemflow and canopy interception were examined (Figure 4), with the best fitting model and parameters of the model listed in Table 5. Throughfall and stemflow were positively related to the rainfall amount. According to the simulated model, throughfall would generate in the 40-year-old stand when rainfall amount reached more than 1.5 mm, and that for the 50-year-old and the 60-year-old stand was 1.9 and 1.86 mm, respectively. The stemflow occurred when rainfall was more than 4.9 mm for the 40-year-old stand, and the corresponding value for the 50-year-old, and the 60-year-old stand was 4.0 and 5.9 mm, respectively. That meant for 40-year-old stand, in general, the canopy was wetting under rainfall less than 1.5 mm. Forest canopy was saturated and could not hold more water when the rainfall reached 1.5 mm, after that throughfall occurred, and when the rainfall more than 4.9 mm, the water flowed down from the stem namely stemflow. This progress was basically consistent with monitored data. The interception and rainfall had binominal function relationship, with coefficient of determination > 0.67. Interception amount increased with the rainfall amount at first, after that reached to a maximum value then declined as shown in figures.

3.4. Effect of Rainfall Intensity on Rainfall Partitioning

Relative throughfall, stemflow and interception to gross rainfall varied with the rainfall intensity in three Chinese pine stand ages (Figure 5). On the whole, as the rainfall intensity increased, the average Tf% and Sf% for all three stands increased, but I% decreased. The largest interception ratio was under the light rain level for all three stands. Except for rainfall storm class I% of the 60-year-old stand was dominantly larger than the 40-year-old stand, and Tf% was opposite.

3.5. Effect of Canopy Features on Rainfall Partitioning

The rainfall redistribution of the three stands was also affected by the canopy structure (Table 1). As the age of Chinese pine forests increased, the DBH, the tree dimension and the LAI increased, whereas the forest gap fraction, openness and the mean leaf angle decreased. The canopy interception ability of the 60-year-old stand is higher than that of 40- and 50-year-old stand, probably caused by the discrepancy of above mentioned structural characteristics. This result demonstrated that the stand structure in different growing stage of the Chinese pine forest was an important factor for rainfall partitioning.

3.6. Comprehensive Analysis of Factors

A stepwise regression analysis was conducted to understand key factors influencing the rainfall redistribution pattern. We used the measured factors including rainfall amount, rainfall intensity, LAI, density, average tree height, DBH, dimension, and mean leaf angle as independent variables, and used Tf, Sf and I as dependent variables, respectively. Results (Table 6) showed rainfall amount and structural properties (LAI) were dominant factors for Tf and I, while rainfall amount, density and mean leaf angle were dominant factors for Sf.

4. Discussion

4.1. Age Dependence of Rainfall Partitioning

Our results showed a rainfall partitioning pattern in Chinese pine stand ages (40, 50 and 60-year-old). We found that throughfall decreased while the canopy interception increased with increasing age, this pattern was in accordance with conifer species rainfall partitioning trend from young to adult in temperate boreal forests [13] and with deciduous forest A. altissima [31]. The study on a series of succession stands of Liriodendron tulipifera forest in the southern Appalachian Mountains also found that canopy interception increased rapidly with the forest age until reaching a maximum of 21% [27]. Another study on different Chinese fir stand ages got the similar trend that interception ratio in mature stands significantly higher [29]. There seems to be a clear pattern, therefore, that in mature forests of many species the canopy interception ability is higher than that in young forests.
Throughfall, stemflow and interception ratio varied among the three stand ages of Chinese pine forests, all of which were within the range of results of Sun et al. reported [14]. The throughfall in our study was higher than that of A. altissima [31], although lower than redcedar in the USA with an average annual throughfall of 57.3% [8]. We suspect that the discrepancy among different plants is due to differences in climate, tree species characteristics and stand structure.
Compared with previous studies of Chinese pine forest on rainfall partitioning (Table 7), including reports at Miyun [34], at Daqingshan [35], at Hebei [36], at Taiyue Mountain [37] and Yeheshan [38]. I% in this study were higher than Taiyue Mountain reflected the young stand age had less canopy interception ability. Throughfall, stemflow and interception ratio of 33 years old Chinese pine forest in Miyun was similar to corresponding values for the 60-year-old stand in this paper, and Yeheshan result was similar to the 40-year-old forest in our study.
Compared with other coniferous forests worldwide (Table 7), the interception ability in our result was higher than that in mixed forests [39] and P. cembroides forest [9], lower than P. eldarica [40], Pinus nigra in Slovenia [41] and Larix gmelinii forest [42], and similar with other reports listed in the table [43,44]. Throughfall in this paper was higher than P. nigra and Chamaecyparis obtuse forest, but smaller than P. cembroides and mixed (upland) forest, within results range of other reports. The stemflow in the 50-year-old stand was larger than in the other two stand ages, though stemflow accounted for a small part of the rainfall (<2%) in all stands. This low percentage is consistent with the results of Chinese fir plantations that stemflow was the highest in the middle age category (12 years old forest). Stemflow of Chinese pine forest in our study was lower than most of the reports listed in the table, much lower than C. obtuse (22.6%), and higher than L. gmelinii and P. nigra forest. The low stemflow in this study may due to climate and also the accuracy of the experimental device of the stemflow in this study. Previous research confirmed that stemflow was the most variable part of the canopy partitioning (coefficient of variation = 107.8%) [17]. Furthermore, it is known that shrubs generate more stemflow than trees, which has been demonstrated in our results comparing with the results of stemflow in shrub species [45,46].

4.2. Effect of Rainfall Characteristic on Rainfall Partitioning

Rainfall directly influences the rainfall redistribution pattern, in redcedar forest rainfall amount accounts for 60–70% of the variation of the percentage of throughfall and stemflow [8]. In our study rainfall amount is also a dominant factor. According to Table 5, linear regression equations between the rainfall amount (x) and throughfall amount (y) showed positive coefficients of determination (R2) > 0.97, whereas for stemflow coefficients of determination were R2 > 0.79. This result was consistent with previous studies [47,48]. In a single rainfall event canopy gradually becomes moist, until rainfall amount reaches a threshold above which throughfall and stemflow are generated. Binominal function model fitted in this study accurately capture the relationship between the canopy interception and rainfall amount, which was different from previous study result that either the power [19] or the linear [49] relationship. The interception amount increased with rainfall until it reached plateaus. Some previous studies showed different best fit models probably due to the experiment designs and the vegetation features.
Rainfall intensity is also an important factor for the rainfall partitioning pattern. On the whole our results showed that throughfall and stemflow increased with rainfall intensity, while the canopy interception decreased, which meant that rainfall intensity could make more net rainfall under the canopy. The larger interception under light rainfall events in this study demonstrated that forest canopy had a strong interception effect under light rainfall events. Some reports reveal that rainfall partitioning was most influenced by rainfall intensity. For example, Nytch et al. [16] reported that rainfall intensity had a significant effect on delay in Tf reaching the ground and Li et al. found that the maximum and minimum interception storage were only correlated with rainfall intensity in Chinese pine forest [50]. However, Liang et al. found that stemflow generation was mainly affected by canopy structure not by rainfall intensity [51], and Li et al. [46] showed the rainfall intensity increased stemflow with which less than 2 mmh−1 and decreased with which greater than 2 mmh−1. In light of this variation among studies, more research should be done to understand the impact of rainfall intensity on rainfall redistribution.

4.3. Forest Structure Dependence of Rainfall Partitioning

Our results showed that the rainfall partitioning pattern in stands with three different stand ages was correlated with the canopy structure. Under the same rainfall events, the variation of throughfall, stemflow and interception were caused by the difference in the crown structure at different growth stage. The LAI of the 60-year-old stand were larger than the 40-year-old stand, meaning that older trees had larger crowns than younger trees, and are therefore able to intercept more water. Previous studies also analyzed forest structural parameters including basal area, bark sickness, leaf area index across the forest, which can contribute to the rainfall redistribution [52,53,54]. For example, larger canopy surface increases interception ratio, and a denser mid-canopy and taller trees can increase cloud interception [1]. Some studies have found that forest structure and meteorological events importantly influence interception ratios, especially in drier sites [23]. Structural properties including LAI and aggregation index were correlated with canopy interception [14]. Stand density is an informative structural factor for all the rainfall partitioning components in coniferous plantations [22]. Smooth barks in some species perform higher stemflow than with rough bark [13]. Significant changes in rainfall partitioning across an invasive chronosequence of A. altissima were underpinned by canopy structure and trunk parameters [31]. Fang et al. deemed canopy structure as a key element for spatial variation of throughfall in small rainfall events [55]. Our results also revealed that forest structure was an important factor influencing rainfall redistribution as shown in previous studies. Through this study, we show that stand age is a determining factor, causing the divergence of canopy structures that lead to differences in rainfall distribution patterns. Further work should be done to discover the rainfall partitioning laws in a large range of forest ages. Above all, both the rainfall amount and stand structure significantly influenced rainfall partitioning.
The last point, we simulated the relationship between the rainfall amount with Tf, Sf and I, and calculated the tipping points, this is very useful and convenient. But advanced methods should be used like rainfall interception modeling, which includes more intercept progresses and can give deep insight to interception mechanisms [56]. Liang et al. concluded that the revised gash model was able to accurately simulate the weekly canopy interception of Chinese pine forest [36]. The results of canopy interception of Chinese pine from model [36,38] were compared with our results, a similar pattern was found, but more useful information from the models that should be our future work.

5. Conclusions

Our study monitored rainfall partitioning pattern of Chinese pine forest at three stand ages and provided deep insights to understand the rainfall interception dynamics. Specifically, we analyzed the dynamics of three rainfall partitioning aspects (throughfall, stemflow, and interception) reflecting the relationship with rainfall amount and intensity, which helps to better understand the progress and relevant facets. Across three stand ages, the majority of rainfall became throughfall, though this decreased with increasing stand age. We attribute this decrease in throughfall to an increase in interception as canopy structure changes with forest age. Our results also highlighted that the rainfall amount is also an important factor affecting the rainfall partitioning across different stand age of Chinese pine. This study provides insight into how rainfall amount, stand age and elements of forest structure affect rainfall partitioning, which has important implications for water conservation and forest management.

Author Contributions

H.H. and F.K. designed the experiments; J.Z., X.S. investigated the field data; L.D. analyzed and wrote the paper; X.C. revised the paper and modified the language. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Special Fund for Forestry Scientific Research in the Public Interest (No. 201204101) and the CFERN & Bejing Techno Solution Award Funds on excellent academic achievements.

Acknowledgments

We would like to thank Stephen Gaughran at Yale University for his assistance with the English language and grammatical editing. We thank Yitao Li at Beijing Forestry University for grammar checking. We appreciate the help of the editor and three anonymous reviewers for their constructive comments on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brauman, K.A.; Freyberg, D.L.; Daily, G.C. Forest structure influences on rainfall partitioning and cloud interception: A comparison of native forest sites in Kona, Hawai’i. Agric. For. Meteorol. 2010, 150, 265–275. [Google Scholar] [CrossRef]
  2. Llorens, P.; Domingo, F. Rainfall partitioning by vegetation under Mediterranean conditions. A review of studies in Europe. J. Hydrol. 2007, 335, 37–54. [Google Scholar] [CrossRef]
  3. Ellison, D.; Futter, M.N.; Bishop, K. On the forest cover–water yield debate: From demand-to supply-side thinking. Glob. Chang. Biol. 2012, 18, 806–820. [Google Scholar] [CrossRef] [Green Version]
  4. Carlylemoses, D.E.; Gash, J.H.C. Rainfall Interception Loss by Forest Canopies. For. Hydrol. Biogeochem. 2011, 216, 407–423. [Google Scholar]
  5. Wang, D.; Wang, L. Canopy interception of apple orchards should not be ignored when assessing evapotranspiration partitioning on the loess plateau in China. Hydrol. Process. 2019, 33, 372–382. [Google Scholar] [CrossRef]
  6. Corti, G.; Agnelli, A.; Cocco, S.; Cardelli, V.; Masse, J.; Courchesne, F. Soil affects throughfall and stemflow under Turkey oak (Quercus cerris L.). Geoderma 2019, 333, 43–56. [Google Scholar] [CrossRef]
  7. Pryor, S.C.; Barthelmie, R.J. Liquid and chemical fluxes in precipitation, throughfall and stemflow: Observations from a deciduous forest and a red pine plantation in the midwestern USA. Water Air Soil Pollut. 2005, 163, 203–227. [Google Scholar] [CrossRef]
  8. Zou, C.B.; Caterina, G.L.; Will, R.E.; Stebler, E.; Turton, D. Canopy interception for a tallgrass prairie under juniper encroachment. PLoS ONE 2015, 10, e0141422. [Google Scholar] [CrossRef]
  9. Pérez-Suárez, M.; Arredondo-Moreno, J.T.; Huber-Sannwald, E.; Serna-Pérez, A. Forest structure, species traits and rain characteristics influences on horizontal and vertical rainfall partitioning in a semiarid pine-oak forest from central Mexico. Ecohydrology 2014, 7, 532–543. [Google Scholar] [CrossRef]
  10. Yang, X.G.; Chen, L.; Wang, L.; Wang, X.; Gu, J.L.; Qu, W.J.; Song, N.P. Dynamic rainfall-partitioning relationships among throughfall, stemflow, and interception loss by Caragana intermedia. J. Hydrol. 2019, 574, 980–989. [Google Scholar] [CrossRef]
  11. He, Z.B.; Yang, J.J.; Du, J.; Zhao, W.Z.; Liu, H.; Chang, X.X. Spatial variability of canopy interception in a spruce forest of the semiarid mountain regions of China. Agric. Meteorol. 2014, 188, 58–63. [Google Scholar] [CrossRef]
  12. Fathizadeh, O.; Attarod, P.; Keim, R.F.; Stein, A.; Amiri, G.Z.; Darvishsefat, A.A. Spatial heterogeneity and temporal stability of throughfall under individual Quercus brantii trees. Hydrol. Process. 2014, 28, 1124–1136. [Google Scholar] [CrossRef]
  13. Barbier, S.; Balandier, P.; Gosselin, F. Influence of several tree traits on rainfall partitioning in temperate and boreal forests: A review. Ann. For. Sci. 2009, 66, 602. [Google Scholar] [CrossRef] [Green Version]
  14. Sun, J.M.; Yu, X.X.; Wang, H.N.; Jia, G.D.; Zhao, Y.; Tu, Z.H.; Deng, W.P.; Jia, J.B.; Chen, J.G. Effects of forest structure on hydrological processes in China. J. Hydrol. 2018, 561, 187–199. [Google Scholar] [CrossRef]
  15. Holder, C.D. Rainfall interception and fog precipitation in a tropical montane cloud forest of Guatemala. For. Ecol. Manag. 2004, 190, 373–384. [Google Scholar] [CrossRef]
  16. Nytch, C.J.; Meléndez-Ackerman, E.J.; Pérez, M.E.; Ortiz-Zayas, J.R. Rainfall interception by six urban trees in San Juan, Puerto Rico. Urban Ecosyst. 2019, 22, 103–115. [Google Scholar] [CrossRef] [Green Version]
  17. Magliano, P.N.; Whitworth-Hulse, J.I.; Baldi, G. Interception, throughfall and stemflow partition in drylands: Global synthesis and meta-analysisa. J. Hydrol. 2019, 568, 638–645. [Google Scholar] [CrossRef]
  18. Aydm, M.; Sen, S.G.; Celik, S. Throughfall, stemflow, and interception characteristics of coniferous forest ecosystems in the western black sea region of turkey (daday example). Environ. Monit. Assess. 2018, 190, 316. [Google Scholar]
  19. Sheng, H.C.; Cai, T.J. Influence of Rainfall on Canopy Interception in Mixed Broad-Leaved—Korean Pine Forest in Xiaoxing’an Mountains, Northeastern China. Forests 2019, 10, 248. [Google Scholar] [CrossRef] [Green Version]
  20. Cuartas, L.A.; Tomasella, J.; Nobre, A.D.; Hodnett, M.G.; Waterloo, M.J.; Múnera, J.C. Interception water-partitioning dynamics for a pristine rainforest in Central Amazonia: Marked differences between normal and dry years. Agric. For. Meteorol. 2007, 145, 69–83. [Google Scholar] [CrossRef]
  21. Frasson, R.P.D.; Krajewski, M.W.F. Characterization of the drop-size distribution and velocity–diameter relation of the throughfall under the maize canopy. Agric. For. Meteorol. 2011, 151, 1244–1251. [Google Scholar] [CrossRef]
  22. Jeong, S.; Otsuki, K.; Farahnak, M. Stand structures and rainfall partitioning in dense unmanaged Japanese cypress plantations. J. Agric. Meteorol. 2019, 5, 92–102. [Google Scholar] [CrossRef] [Green Version]
  23. Del Campo, A.D.; González-Sanchis, M.; Lidón, A.; Ceacero, C.J.; García-Prats, A. Rainfall partitioning after thinning in two low-biomass semiarid forests: Impact of meteorological variables and forest structure on the effectiveness of water-oriented treatments. J. Hydrol. 2018, 565, 74–86. [Google Scholar] [CrossRef]
  24. Honda, E.A.; Mendonça, A.H.; Durigan, G. Factors affecting the stemflow of trees in the Brazilian Cerrado. Ecohydrology 2015, 8, 1351–1362. [Google Scholar] [CrossRef]
  25. Zhang, Y.F.; Wang, X.P.; Pan, Y.X.; Hu, R. Relative contribution of biotic and abiotic factors to stemflow production and funneling efficiency: A long-term field study on a xerophytic shrub species in Tengger Desert of northern China. Agric. For. Meteorol. 2020, 280, 107781. [Google Scholar] [CrossRef]
  26. Pypker, T.G.; Bond, B.J.; Link, T.E.; Marks, D.; Unsworth, M.H. The importance of canopy structure in controlling the interception loss of rainfall: Examples from a young and an old-growth Douglas-fir forest. Agric. Meteorol. 2005, 130, 113–129. [Google Scholar] [CrossRef]
  27. Brantley, S.T.; Miniat, C.F.; Bolstad, P.V. Rainfall partitioning varies across a forest age chronosequence in the southern appalachian mountains. Ecohydrology 2019, 12, e2081. [Google Scholar] [CrossRef]
  28. Calvo-Alvarado, J.C.; Jiménez-Rodríguez, C.D.; Calvo-Obando, A.J.; do Espírito-Santo, M.M.; Gonçalves-Silva, T. Interception of rainfall in successional tropical dry forests in Brazil and Costa Rica. Geosciences 2018, 8, 486. [Google Scholar] [CrossRef] [Green Version]
  29. Yan, W.; Deng, X.; Chen, X.; Tian, D.; Xiang, W.; Peng, Y. Long-term variations of rainfall interception in different growth stages of Chinese fir plantations. Hydrol. Sci. J. 2015, 60, 2178–2188. [Google Scholar] [CrossRef]
  30. Zhao, J.; Kang, F.; Wang, L.; Yu, X.; Zhao, W. Patterns of biomass and carbon distribution across a chronosequence of Chinese pine (Pinus tabulaeformis) forests. PLoS ONE 2014, 9, e94966. [Google Scholar] [CrossRef]
  31. Sadeghi, S.M.M.; Van Stan, J.T.; Pypker, T.G.; Friesen, J. Canopy hydrometeorological dynamics across a chronosequence of a globally invasive species, Ailanthus altissima (Mill., tree of heaven). Agric. For. Meteorol. 2017, 240, 10–17. [Google Scholar] [CrossRef]
  32. Sun, X.; Onda, Y.; Kato, H.; Gomi, T.; Liu, X. Estimation of throughfall with changing stand structures for Japanese cypress and cedar plantations. For. Ecol. Manag. 2017, 402, 145–156. [Google Scholar] [CrossRef]
  33. Carlyle-Moses, D.E.; Price, A.G. An evaluation of the Gash interception model in a northern hardwood stand. J. Hydrol. 1999, 214, 103–110. [Google Scholar] [CrossRef]
  34. Xiao, Y.; Chen, L.; Yu, X.; Yang, X.; Sun, Q. Influence on precipitation distribution of Pinus tabuleaefomis forest in Miyun Reservoir. J. Soil Water Conserv. 2007, 21, 154–157. [Google Scholar]
  35. Sharentu, Y.; Tian, Y.; Guo, L. Study on characteristics of distribution of rainfall in Pinus tabulaeformis plantation in Daqing-Mountains. J. Arid Land Resour. Environ. 2009, 23, 157–160. [Google Scholar]
  36. Liang, W. Simulation of Gash Model to Rainfall Interception of Pinus tabulaeformis. For. Syst. 2014, 23, 300–303. [Google Scholar] [CrossRef] [Green Version]
  37. Zhou, B.; Han, H.; Kang, F.; Cheng, X.; Song, Y.; Liu, K.; Li, Y. Characteristics of precipitation distribution in Pinus tabulaeformis plantations under different canopy coverage in Taiyue Mountain. Acta Ecol. Sin. 2013, 33, 1645–1653. [Google Scholar] [CrossRef]
  38. Ma, C.; Li, X.; Luo, Y.; Shao, M.; Jia, X. The modelling of rainfall interception in growing and dormant seasons for a pine plantation and a black locust plantation in semi-arid Northwest China. J. Hydrol. 2019, 577, 123849. [Google Scholar] [CrossRef]
  39. Kermavnar, J.; Vilhar, U. Canopy precipitation interception in urban forests in relation to stand structure. Urban Ecosyst. 2017, 20, 1373–1387. [Google Scholar] [CrossRef]
  40. Sadeghi, S.M.M.; Attarod, P.; Grant Pypker, T.; Dunkerley, D. Is canopy interception increased in semiarid tree plantations? Evidence from a field investigation in Tehran, Iran. Turk. J. Agric. For. 2014, 38, 792–806. [Google Scholar] [CrossRef]
  41. Zabret, K.; RakovecM, J.; Šraj, M. Influence of meteorological variables on rainfall partitioning for deciduous and coniferous tree species in urban area. J. Hydrol. 2018, 558, 29–41. [Google Scholar] [CrossRef]
  42. Jiang, H.; Zhao, Y.; Xin, X.; Ma, W.; Li, X.; Sun, C.; Guo, X. Precipitation Distribution Laws of Different Forest Stand in Daxing’anling. J. Soil Water Conser. 2008, 22, 197–201. [Google Scholar]
  43. Saito, T.; Matsuda, H.; Komatsu, M.; Xiang, Y.; Takahashi, A.; ShinoharaK, Y.; Otsuki, K. Forest canopy interception loss exceeds wet canopy evaporation in Japanese cypress (Hinoki) and Japanese cedar (Sugi) plantations. J. Hydrol. 2013, 507, 287–299. [Google Scholar] [CrossRef]
  44. Iida, S.i.; Levia, D.F.; Shimizu, A.; Shimizu, T.; Tamai, K.; Nobuhiro, T.; Kabeya, N.; Noguchi, S.; Sawano, S.; Araki, M. Intrastorm scale rainfall interception dynamics in a mature coniferous forest stand. J. Hydrol. 2017, 548, 770–783. [Google Scholar] [CrossRef]
  45. Zhang, Y.F.; Wang, X.P.; Hu, R.; Pan, Y.X.; Zhang, H. Stemflow in two xerophytic shrubs and its significance to soil water and nutrient enrichment. Ecol. Res. 2013, 28, 567–579. [Google Scholar] [CrossRef]
  46. Li, X.Y.; Liu, L.Y.; Gao, S.Y.; Ma, Y.J.; Yang, Z.P. Stemflow in three shrubs and its effect on soil water enhancement in semiarid loess region of China. Agric. For. Meteorol. 2008, 148, 1501–1507. [Google Scholar] [CrossRef]
  47. Rahmani, R.; Sadoddin, A.; Ghorbani, S. Measuring and modelling precipitation components in an Oriental beech stand of the Hyrcanian region, Iran. J. Hydrol. 2011, 404, 294–303. [Google Scholar] [CrossRef]
  48. Zimmerman, A.; Zimmermann, B.; Elsenbeer, H. Rainfall redistribution in a tropical forest: Spatial and temporal patterns. Water Resour. Res. 2009, 45, W11413. [Google Scholar] [CrossRef] [Green Version]
  49. Gardon, F.R.; Toledo, R.M.D.; Brentan, B.M.; Santos, R.F.D. Rainfall interception and plant community in young forest restorations. Ecol. Indic. 2020, 109, 105779. [Google Scholar] [CrossRef]
  50. Li, X.; Xiao, Q.; Niu, J.; Dymond, S.; van Doorn, N.S.; Yu, X.; Xie, B.; Lv, X.; Zhang, K.; Li, J. Process-based rainfall interception by small trees in Northern China: The effect of rainfall traits and crown structure characteristics. Agric. For. Meteorol. 2016, 218, 65–73. [Google Scholar] [CrossRef]
  51. Liang, W.; Kosugi, K.; Takahisa, M. Characteristics of stemflow for tall stewartia (Stewartia monadelpha) growing on a hillslope. J. Hydrol. 2009, 378, 168–178. [Google Scholar] [CrossRef] [Green Version]
  52. Sadeghi, S.M.M.; Stan, J.T.V.; Pypker, T.G.; Tamjidi, J.; Farahnaklangroudi, M. Importance of transitional leaf states in canopy rainfall partitioning dynamics. Eur. J. For. Res. 2018, 1, 1–10. [Google Scholar] [CrossRef]
  53. César Cisneros Vaca, C.; van der Tol, C.; Ghimire, C.P. The influence of long-term changes in canopy structure on rainfall interception loss: A case study in speulderbos, the netherlands. Hydrol. Earth Syst. Sci. 2018, 22, 1–42. [Google Scholar] [CrossRef] [Green Version]
  54. Liu, C.; Wang, Y.; Ma, C.; Wang, Y.; Zhang, H.; Hui, B. Quantifying the effect of non-spatial and spatial forest stand structure on rainfall partitioning in mountain forests, southern China. For. Chron. 2018, 94, 162–172. [Google Scholar] [CrossRef]
  55. Fang, S.; Zhao, C.; Jian, S. Spatial variability of throughfall in a Pinus tabulaeformis plantation forest in Loess Plateau, China. Scand. J. Forest Res. 2016, 31, 467–476. [Google Scholar] [CrossRef]
  56. Muzylo, A.; Llorens, P.; Valente, F.; Keizer, J.J.; Domingo, F.; Gash, J.H.C. A review of rainfall interception modelling. J. Hydrol. 2009, 370, 191–206. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of Liaoheyuan Natural Reserve, Hebei Province, China.
Figure 1. Location of Liaoheyuan Natural Reserve, Hebei Province, China.
Forests 11 00243 g001
Figure 2. The photos for the experiment plot, (a): the auto precipitation gauge in the study, (b): throughfall gauge used in this study, (c): stemflow gauge in this study.
Figure 2. The photos for the experiment plot, (a): the auto precipitation gauge in the study, (b): throughfall gauge used in this study, (c): stemflow gauge in this study.
Forests 11 00243 g002
Figure 3. The amount of throughfall, stemflow and interception in five different months during the study period.
Figure 3. The amount of throughfall, stemflow and interception in five different months during the study period.
Forests 11 00243 g003
Figure 4. Linear relationship between throughfall, stemflow, interception with gross rainfall for three Chinese pine forest stands. Total of 64 rainfall events. △: 40-year-old stand, ○: 50-year-old stand, ×: 60-year-old stand. —: 40-year-old, ----: 50-year-old, —··—·: 60-year-old.
Figure 4. Linear relationship between throughfall, stemflow, interception with gross rainfall for three Chinese pine forest stands. Total of 64 rainfall events. △: 40-year-old stand, ○: 50-year-old stand, ×: 60-year-old stand. —: 40-year-old, ----: 50-year-old, —··—·: 60-year-old.
Forests 11 00243 g004
Figure 5. Relative throughfall, stemflow, and interception categorized by rainfall classes for three Chinese pine stand ages. Different letters meant the significant difference (p < 0.05) among the three stand ages, and the alphabetical order showed the values from the largest to the smallest.
Figure 5. Relative throughfall, stemflow, and interception categorized by rainfall classes for three Chinese pine stand ages. Different letters meant the significant difference (p < 0.05) among the three stand ages, and the alphabetical order showed the values from the largest to the smallest.
Forests 11 00243 g005
Table 1. Basic structural characteristics of the three sampling plots of the Chinese pine forest of 40, 50 and 60 years old of age.
Table 1. Basic structural characteristics of the three sampling plots of the Chinese pine forest of 40, 50 and 60 years old of age.
Age (Years)Slope (°)Density (Trees hm−2)Mean DBH (cm)Mean Height (m)LAI (m2/m2)GapFraction (%)Openness (%)Mean Leaf Angle (Degree)
4027101716.111.32.514.115.130.6
501763420.615.62.613.314.227.2
601843430.420.72.811.011.624.1
Table 2. Basic information of 15 monitored trees for stemflow in 40-, 50-, 60-year-old Chinese pine stands.
Table 2. Basic information of 15 monitored trees for stemflow in 40-, 50-, 60-year-old Chinese pine stands.
Stand Age (Years)Stemflow Tree Number (n)DBH (cm)Crown (m × m)
40113.83.5 ×3.9
217.84.5 × 4.8
320.64.8 × 5.6
424.64.9 × 5.5
529.84.9 × 6.0
5017.52.0 × 1.5
214.23.4 × 2.8
320.15.1 × 5.8
426.25.6 × 9.2
532.47.5 × 8.9
60124.85.8 × 5.4
227.74.8 × 6.3
333.56.3 × 6.8
436.27.9 × 7.2
540.810.3 × 5.8
Table 3. The total rainfall from May to September of 2013–2014 according to four rainfall classes.
Table 3. The total rainfall from May to September of 2013–2014 according to four rainfall classes.
Rainfall Intensity Class
(mm d−1)
201320142013–2014
Amount
(mm)
Frequency
(N)
Amount
(mm)
Frequency
(N)
Amount
(mm)
Frequency
(N)
0.2 ≤ P < 10113.4238117194.440
10 ≤ P < 25117.86120.88238.614
25 ≤ P < 50129.43178.45307.88
50 ≤ P < 100116.6200116.62
Total rainfall477.234380.230857.464
Table 4. Cumulative of throughfall, stemflow and interception in mm and percentages with concerning gross precipitation of three stand ages of Chinese pine forest from May to September of 2013 and 2014.
Table 4. Cumulative of throughfall, stemflow and interception in mm and percentages with concerning gross precipitation of three stand ages of Chinese pine forest from May to September of 2013 and 2014.
Stand Age (Years)Tf (mm)Tf%Sf (mm)Sf%I (mm)I%
40675.678.86.50.8175.320.4
50635.174.19.61.1212.724.8
60571.866.74.70.6280.832.8
Table 5. Linear relationship between throughfall, stemflow, interception with gross rainfall in three Chinese pine forest stands, with a total of 64 rainfall events.
Table 5. Linear relationship between throughfall, stemflow, interception with gross rainfall in three Chinese pine forest stands, with a total of 64 rainfall events.
ItemModelRegression Parameters
FPSE
Tfy40 = 0.8880x − 1.34139778.1220.0001.0270.997
y50 = 0.8639x − 1.64006188.5250.0001.2560.995
y60 = 0.7744x − 1.44002168.4310.0001.9030.972
Sfy40= 0.0120x − 0.0589360.9830.0000.0720.851
y50 = 0.0159x − 0.0642361.0280.0000.0960.851
y60 = 0.0099x − 0.0582232.5730.0000.0740.786
Iy40 = −0.0008x2 + 0.1406x + 1.167665.7930.0001.0040.673
y50 = −0.0025x2 + 0.2468x + 0.989175.5550.0001.1740.703
y60 = −0.0032x2 + 0.3765x + 0.578495.6500.0001.8300.750
Table 6. Stepwise analysis of rainfall partitioning factors.
Table 6. Stepwise analysis of rainfall partitioning factors.
ItemInfluencing FactorModelN
TfRainfall amount(x1), LAI(x2)y = 12.574 + 0.842x1 − 5.335x2192
SfRainfall amount(x1), density(x2), mean leaf angle(x3)y = 0.013x1 − 0.001x2 + 0.075x3 − 1.567192
IRainfall amount(x1), LAI(x2)y = 0.145x1 + 5.472x2 − 12.871192
Table 7. Rainfall partitioning characteristics of conifer forests in different locations.
Table 7. Rainfall partitioning characteristics of conifer forests in different locations.
Location Forest TypeAge (Years)Study PeriodTf%Sf%I%Literature
Boreal and temperate biogeographical zonesconifers--64.5 04.1028.7 0Barbier et al. 2009
Taiyue Mountain, ChinaPinus tabulaeformis20May to September 201177.892.3019.92Zhou et al. 2013
Miyun Reservoir, ChinaPinus tabulaeformis332004–200667.650.6831.67Xiao et al. 2007
Yeheshan watershed, ChinaPinus tabulaeformis17January–December 201675.400.7023.90Ma et al. 2019
Chitgar Forest Park, IranPinus eldarica4230 January 2011 to 30 January 2012 44.20Sadeghi et al. 2014
urban area of Ljubljana, SloveniaPinus nigra-January 2014 to June 201728.000.0272.00Zabret, Rakovec and Šraj 2018
Bezirgan Basin, TurkeyPinus nigra-2012–201469.802.6027.70Aydm, Şen and Celik 2018
Bezirgan Basin, TurkeyPinus sylvestris-2012–201473.95.9020.20Aydm, Şen and Celik 2018
Bezirgan Basin, Turkeymixed Pinus nigra–Pinus sylvestris-2012–201477.703.1019.20Aydm, Şen and Celik 2018
Daxing’anling, ChinaPinus sylvestris var. mongolica3218 May–September 11 200768.510.5431.00Jiang et al.2008
Sierra San Miguelito, MexicoPinus cembroides81June 2006 to July 200992.233.7618.95Pérez-Suárez et al.2014
Daxing’anling, ChinaPicea koraiensis2018 May–11 September 200766.001.04 33.50Jiang et al. 2008
Daxing’anling, ChinaLarix gmelinii32May 18–September 11 200759.760.3540.20Jiang et al. 2008
Fukuoka, JapanChamaecyparis obtusa33April to October 201749.4022.6028.00Jeong, Kyoichi and Farahnak, 2019
Fukuoka, JapanChamaecyparis obtusa415 June 2010 to 31 December 201165.309.1025.50Saito et al. 2013
Fukuoka, JapanCryptomeria japonica415 June 2010 to 31 December 201167.906.6025.50Saito et al. 2013
Tsukuba Experimental Watershed, JapanCryptomeria japonica502008 and 201076.903.2019.90Iida et al. 2017
Changbai Mountains, Chinamixed forest of broad-leaved and Pinus koaiensis-growing seasons in 2010 and 2011 21.95Sheng and Cai 2019
Ljubljana, Slovenia mixed (upland) forest including Picea abies-1 January 2008 to 31 December 201394.12.003.90Kermavnar et al. 2017
Chitgar Forest Park, IranCupressus arizonica4230 January 2011 to 30 January 2012 34.40Sadeghi et al. 2014

Share and Cite

MDPI and ACS Style

Dong, L.; Han, H.; Kang, F.; Cheng, X.; Zhao, J.; Song, X. Rainfall Partitioning in Chinese Pine (Pinus tabuliformis Carr.) Stands at Three Different Ages. Forests 2020, 11, 243. https://doi.org/10.3390/f11020243

AMA Style

Dong L, Han H, Kang F, Cheng X, Zhao J, Song X. Rainfall Partitioning in Chinese Pine (Pinus tabuliformis Carr.) Stands at Three Different Ages. Forests. 2020; 11(2):243. https://doi.org/10.3390/f11020243

Chicago/Turabian Style

Dong, Lingling, Hairong Han, Fengfeng Kang, Xiaoqin Cheng, Jinlong Zhao, and Xiaoshuai Song. 2020. "Rainfall Partitioning in Chinese Pine (Pinus tabuliformis Carr.) Stands at Three Different Ages" Forests 11, no. 2: 243. https://doi.org/10.3390/f11020243

APA Style

Dong, L., Han, H., Kang, F., Cheng, X., Zhao, J., & Song, X. (2020). Rainfall Partitioning in Chinese Pine (Pinus tabuliformis Carr.) Stands at Three Different Ages. Forests, 11(2), 243. https://doi.org/10.3390/f11020243

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop