The canopy conductance of a boreal aspen forest,Prince Albert National Park,Canada

Annual fluxes of canopy‐level heat, water vapour and carbon dioxide were measured using eddy covariance both above the aspen overstory (Populus tremuloides Michx.) and hazelnut understory (Corylus cornuta Marsh.) of a boreal aspen forest (53·629 °N 106·200 °W). Partitioning of the fluxes between overstory and understory components allowed the calculation of canopy conductance to water vapour for both species. On a seasonal basis, the canopy conductance of the aspen accounted for 70% of the surface conductance, with the latter a strong function of the forest's leaf area index. On a half‐hour basis, the canopy conductance of both species decreased non‐linearly as the leaf‐surface saturation deficits increased, and was best parameterized and showed similar sensitivities to a modified form of the Ball–Berry–Woodrow index, where relative humidity was replaced with the reciprocal of the saturation deficit. The negative feedback between the forest evaporation and the saturation deficit in the convective boundary layer varied from weak when the forest was at full leaf to strong when the forest was developing or loosing leaves. The coupling between the air at the leaf surface and the convective boundary layer also varied seasonally, with coupling decreasing with increasing leaf area. Compared with coniferous boreal forests, the seasonal changes in leaf area had a unique impact on vegetation–atmosphere interactions. Copyright © 2004 John Wiley & Sons, Ltd.     

Influence of soil moisture hysteresis on the functioning of capillary barriers

Previous experimental studies of capillary barriers have identified highly hysteretic soil moisture retention characteristics in the materials used. In this study, numerical modelling is used to analyse the role of soil moisture hysteresis in capillary barrier functioning. Comparisons between published experimental results and model simulations indicate that soil moisture hysteresis was a necessary inclusion in the modelling approach to adequately reproduce pore water pressure distributions and the timing of breakthrough occurrences. Under hypothetical intermittent infiltration and evaporation conditions, the predicted volumetric water content in the moisture retention layer was significantly different for hysteretic and non‐hysteretic models. The hysteresis effect was found to be dependent on the nature of infiltration–evaporation cycling, although the predicted volume of flow through the hysteretic barrier was lower than that of the non‐hysteretic case, regardless of the nature of the cyclic upper boundary conditions. For practical engineering designs, where the water leakage through the barrier is the primary concern, the inclusion of soil moisture hysteresis in numerical modelling is needed. Copyright © 2009 John Wiley & Sons, Ltd.     

Jet drop enrichment of bacteria,virus, and dissolved organic material

Four parameters that control the enrichment of bacteria in jet drops are bubble scavenging, drop size, drop position in the jet set, and the type of bacteria. Without the scavenging of bacteria as a bubble rises through the water it is doubtful that the observed enrichment factors, EF, greater than 1000 could be obtained. There is a maximum in EF as a function of top jet drop size, and the EF decreases from the top to the bottom drop of the jet set. The efficiency by which bubbles scavenge bacteria varies with species. Presumably these parameters apply in some degree to the EF of virus in jet drops. Dissolved organic material in natural waters can adsorb to bubbles and contribute to a large EF in jet drops, but there is a feedback mechanism whereby changes in bubble surface free energy modify the jet drop-size distribution. However, there is reason to believe this will not significantly influence the jet drop-size distribution produced by bubbles in the sea.     

Temporal fluctuations of water repellency in forest soils of Galicia,NW Spain. Do soil samples dried at laboratory reflect the potential soil water repellency?

Water repellency is a widespread property of Pinus pinaster and Eucalyptus globulus forest soils in NW Spain and is particularly severe during the summer dry conditions. The aim of this work was to compare actual water repellency at field‐moist samples with potential water repellency after drying at 25 and 105 °C in samples collected at different times of year under four forest soils. Also, we investigated whether drying at 25 or 105 °C led to repellency values comparable to the highest levels reached under field conditions in the summer with a view to developing an appropriate sampling protocol towards estimating the maximum possible water repellency of a given soil as a key to establishing its environmental effects. The actual and potential water repellency was determined by using the water drop penetration time (WDPT) and molarity of an ethanol drop (MED) tests. Clear seasonal patterns of water repellency were observed from the results for the four forest soils, peaking in the dry period and disappearing after prolonged wet periods. Water repellency lasts longer in sandy loam soils than in more finely textured soils, and also under eucalyptus than under pine forests. Drying soil samples at 25 or 105 °C increased water repellency, as measured with the WDPT method, in the four soils, but especially in the non‐repellent samples collected during the wet period. The increase was more marked in the sandy loam soils than in the more finely textured soils, and also after drying at 105 °C than at 25 °C. MED measurements exposed a common trait in the four soils; thus, the water repellency values obtained under field conditions in summer invariably exceeded those obtained after drying at 25 or 105 °C. In addition, the repellency values for dried samples collected in the wet period were never comparable to the maximum levels observed under field conditions in the summer. Copyright © 2011 John Wiley & Sons, Ltd.     

The Density Distributions and the Correlations of an Earthquake's Source Parameters

The data analysis of the source parameters of five sets of earthquake sequences, aftershocks and earthquakes scattered in a region shows that the scalar seismic moment is correlated with the linear size of the fault and the static stress drop. We tentatively imply that a correlation also exists between the radius of the faults and the static stress drops and it is suggested that the static stress drop may be decreasing with increasing radius of the source. It is shown that the density distribution of the source radius, calculated through the source rupture duration obtained from the body wave pulse (Boatwright, 1980) using the time T between the P-wave onset and the first zero crossing on the seismogram, may be represented by a power law as the density distribution of the stress drops and of the moment which are also computed. It is also suggested, and tentatively verified, that the density distribution of the areas of the broken barriers on the faults is similar to that of the density distribution of the static stress drops. It is finally suggested that the seismicity of a region may be studied two-dimensionally as a function of the stress drop and the radius of the source instead of the classic b and b₀ values. Concerning the discussion on the range of the values of the static stress drop, whether it is almost constant in a seismic region and varies only from one region to another, it is seen that in the aftershocks of the 1994 Northridge earthquake it covers a range of almost 5 orders of magnitudes. Finally it is ascertained that the density distribution of the source parameters does not give equipartition of seismic moment release.     

Pore-scale contact angle measurements at reservoir conditions using X-ray microtomography

Contact angle is a principal control of the flow of multiple fluid phases through porous media; however its measurement on other than flat surfaces remains a challenge. A new method is presented for the measurement of the contact angle between immiscible fluids at the pore scale at reservoir conditions (10 MPa and 50 °C) inside a quarry limestone through the use of X-ray microtomography. It is applied to a super-critical CO₂–brine–carbonate system by resampling the micro-CT data onto planes orthogonal to the contact lines, allowing for vectors to be traced along the grain surface and the CO₂–brine interface. A distribution of contact angles ranging from 35° to 55° is observed, indicating that the CO₂–brine–carbonate system is weakly water-wet. This range of contact angles can be understood as the result of contact angle hysteresis and surface heterogeneity on a range of length scales. Surface heterogeneity is examined by comparison of micro-CT results with optical thin sections and SEM images.     

The impact of organic material on cloud and fog processes

Airborne surface-active organic substances can form adsorbed films at air-water and air-solid interfaces in the atmosphere. On the basis of considerable laboratory research, it has been suggested that these films can modify the microphysical processes of clouds and fog. The possible alterations of physical processes include retardation of the rate of growth and evaporation of water drops, passivation of cloud condensation nuclei, modification of drop coalescence efficiency, and the promotion of drop instability. The chemical structure of the film-forming compounds and the physical properties of the adsorbed films determine the degree to which the physical processes are modified. Upon consideration of the actual organic composition of the atmosphere, it was concluded that natural surface films will not greatly influence drop growth and evaporation. In addition, there is no definitive evidence from field experiments that other processes are influenced by natural organic films. The potential for useful weather modification through the intentional introduction of film-forming organic substances into the atmosphere was critically evaluated. Although numerous approaches have been suggested by laboratory studies, only the stabilization of airborne drops against evaporation has been unequivocally demonstrated in field experiments.     

地震应力降与岩石破裂应力降

从应力降的静态理论,讨论了地震应力降的精度,并利用野外观测的最大位移和由地震波得到的地震矩,分别计算了十个研究较充分的地震应力降.对不同作者的结果作了比较,可以看出,应力降有因子为2的变化范围.但从统计观点看,地震应力降是分布在几巴到几百巴的数量级,大多数位于20-60巴;总结了地壳中代表性岩石在地壳温度、压力条件下,粘滑及脆性剪切破裂时的应力降,它们与温度、压力、岩性、含水矿物的存在有关.取值范围是从零到几千巴数量级;讨论了地震断层的物理状态,指出地震应力降是一种平均结果,并由于静态理论忽略了断层面破裂所消耗的能量,所以地震应力降比实际的应力降要低.      

三峡水库水面蒸发年内滞后效应分析

水面蒸发在季节变化上相比净辐射等气象要素的相位延迟反映了水体储热对水面蒸发的影响,量化这一对气象要素波动的滞后响应对理解和估算深水水库(湖泊)蒸发非常重要。三峡水库等河道型深水水库的水位和面积具有显著周期性变动,使得其水面蒸发的响应模式更为复杂,而目前对其认识非常薄弱。本文利用2013年8月—2020年7月三峡水库巴东站水面蒸发场和陆面蒸发场的蒸发和气象观测数据,分析了水面蒸发的季节变化及其年内滞后效应。结果表明:水面蒸发场蒸发量在年内的8和12月呈现双峰变化,与只有8月单峰值的陆面蒸发场蒸发量显著不同。水面蒸发场蒸发量相对净辐射、平均气温和水面温度分别存在4、3和2个月的滞后,而陆面蒸发场蒸发量相对滞后时间均在1个月以内;水面与陆面蒸发场相比,水温、蒸发量和水面与大气饱和差之间的滞后时间分别为1、3和4个月,而平均气温和净辐射之间不存在滞后。本文揭示出三峡水库巴东段水面蒸发年内滞后效应主要受到水库水温引起的水面与大气饱和差在季节上滞后的影响,需通过深入分析水温的时空变化来明确整个三峡水库的水面蒸发特征。     

Seasonal variation in the energy and water exchanges above and below a larch forest in eastern Siberia

The water and energy exchanges in forests form one of the most important hydro‐meteorological systems. There have been far fewer investigations of the water and heat exchange in high latitude forests than of those in warm, humid regions. There have been few observations of this system in Siberia for an entire growing season, including the snowmelt and leaf‐fall seasons. In this study, the characteristics of the energy and water budgets in an eastern Siberian larch forest were investigated from the snowmelt season to the leaf‐fall season. The latent heat flux was strongly affected by the transpiration activity of the larch trees and increased quickly as the larch stand began to foliate. The sensible heat dropped at that time, although the net all‐wave radiation increased. Consequently, the seasonal variation in the Bowen ratio was clearly ‘U’‐shaped, and the minimum value (1·0) occurred in June and July. The Bowen ratio was very high (10–25) in early spring, just before leaf opening. The canopy resistance for a big leaf model far exceeded the aerodynamic resistance and fluctuated over a much wider range. The canopy resistance was strongly restricted by the saturation deficit, and its minimum value was 100 s m^?1 (10 mm s^?1 in conductance). This minimum canopy resistance is higher than values obtained for forests in warm, humid regions, but is similar to those measured in other boreal conifer forests. It has been suggested that the senescence of leaves also affects the canopy resistance, which was higher in the leaf‐fall season than in the foliated season. The mean evapotranspiration rate from 21 April 1998 to 7 September 1998 was 1·16 mm day^?1, and the maximum rate, 2·9 mm day^?1, occurred at the beginning of July. For the growing season from 1 June to 31 August, this rate was 1·5 mm day^?1. The total evapotranspiration from the forest (151 mm) exceeded the amount of precipitation (106 mm) and was equal to 73% of the total water input (211 mm), including the snow water equivalent. The understory evapotranspiration reached 35% of the total evapotranspiration, and the interception evaporation was 15% of the gross precipitation. The understory evapotranspiration was high and the interception evaporation was low because the canopy was sparse and the leaf area index was low. Copyright © 2001 John Wiley & Sons, Ltd.     

Long‐term prediction of the water level and salinity in the Dead Sea

The long‐term water level variations of the Dead Sea (DS) were assessed using a previously developed simulation model. The model establishes the condition of the DS by evaluating a series of ordinary differential equations describing mass balances on the water and major chemical species. The DS was modelled as a two‐layer system. The model was modified using up‐to‐date inflow data and recent hypsometric graphs to derive the volume–area–level relationships. Three scenarios were studied: continuation of current conditions; a cessation in industrial activity when the DS water level drops to a certain level; and a simplified weather change scenario. The model predicted that the DS will not dry up, but its level will continue to drop with a decelerating rate with no equilibrium level in 500 years. Changing climate would accelerate the level drop. In the 500 year period, after an initial increase, the DS salinity drops. The opposite behaviour is noted in the evaporation rate, which increases after an initial decrease. Ceasing industrial pumping would eventually restore the DS to its normal level, but with changed conditions. Copyright © 2002 John Wiley & Sons, Ltd.     

Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.     

Source process of the 14 November 2001 western Kunlun Mountain <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 earthquake

Based on digital teleseismic P-wave seismograms recorded by 28 long-period seismograph stations of the global seismic network, source process of the November 14, 2001 western Kunlun Mountain M _S=8.1 (M _W=7.8) earthquake is estimated by a new inversion method. The result shows that the earthquake is a very complex rupture event. The source rupture initiated at the hypocenter (35.95°N, 90.54°E, focal depth 10 km, by USGS NEIC), and propagated to the west at first. Then, in several minutes to a hundred minutes and over a large spatial range, several rupture growth points emerged in succession at the eastern end and in the central part of the finite fault. And then the source rupture propagated from these rupture growth points successively and, finally, stopped in the area within 50 km to the east of the centroid position (35.80°N, 92.91°E, focal depth 15 km, by Harvard CMT). The entire rupture lasted for 142 s, and the source process could be roughly separated into three stages: The first stage started at the 0 s and ended at the 52 s, lasting for 52 s and releasing approximately 24.4% of the total moment; The second stage started at the 55 s and ended at the 113 s, lasting for 58 s and releasing approximately 56.5% of the total moment; The third stage started at the 122 s and ended at the 142 s, lasting for 20 s and releasing approximately 19.1% of the total moment. The length of the ruptured fault plane is about 490 km. The maximum width of the ruptured fault plane is about 45 km. The rupture mainly occurred within 30 km in depth under the surface of the Earth. The average static slip in the underground rocky crust is about 1.2 m with the maximum static slip 3.6 m. The average static stress drop is about 5 MPa with the maximum static stress drop 18 MPa. The maximum static slip and the maximum stress drop occurred in an area within 50 km to the east of the centroid position.     

Contact line friction and surface tension effects on seismic attenuation and effective bulk moduli in rock with a partially saturated crack

The effect of surface phenomena occurring at the interfaces between immiscible fluids and a solid on the seismic attributes of partially saturated rocks has not yet been fully studied. Meanwhile, over the past two decades considerable progress has been made in the physics of wetting to understand effects such as contact line friction, contact line pinning, contact angle hysteresis, and equilibrium contact angle. In this paper, we developed a new rock physics model considering the aforementioned effects on seismic properties of the rock with a partially saturated plane-strain crack. We demonstrated that for small wave-induced stress perturbations, the contact line of the interface meniscus will remain pinned, while the meniscus will bulge and change its shape through the change of the contact angles. When the stress perturbation is larger than a critical value, the contact line will move with advancing or receding contact angle depending on the direction of contact line motion. A critical stress perturbation predicted by our model can be in the range of ∼10²−10⁴ Pa, that is typical for linear seismic waves. Our model predicts strong seismic attenuation in the case when the contact line is moving. When the contact line is pinned, the attenuation is negligibly small. Seismic attenuation is associated with the hysteresis of loading and unloading bulk moduli, predicted by our model. The hysteresis is large when the contact line is moving and negligibly small when the contact line is pinned. Furthermore, we demonstrate that the bulk modulus of the rock with a partially saturated crack depends also on the surface tension and on the contact angle hysteresis. These parameters are typically neglected during calculation of the effecting fluid moduli by applying different averaging techniques. We demonstrate that contact line friction may be a dominant seismic attenuation mechanism in the low frequency limit (<∼10 Hz) when capillary forces dominate over viscous forces during wave-induced two-phase fluid flow.     

Source process of the 14 November 2001 western Kunlun Mountain M S=8.1 earthquake

Based on digital teleseismic P-wave seismograms recorded by 28 long-period seismograph stations of the global seismic network, source process of the November 14, 2001 western Kunlun Mountain M _S=8.1 (M _W=7.8) earthquake is estimated by a new inversion method. The result shows that the earthquake is a very complex rupture event. The source rupture initiated at the hypocenter (35.95°N, 90.54°E, focal depth 10 km, by USGS NEIC), and propagated to the west at first. Then, in several minutes to a hundred minutes and over a large spatial range, several rupture growth points emerged in succession at the eastern end and in the central part of the finite fault. And then the source rupture propagated from these rupture growth points successively and, finally, stopped in the area within 50 km to the east of the centroid position (35.80°N, 92.91°E, focal depth 15 km, by Harvard CMT). The entire rupture lasted for 142 s, and the source process could be roughly separated into three stages: The first stage started at the 0 s and ended at the 52 s, lasting for 52 s and releasing approximately 24.4% of the total moment; The second stage started at the 55 s and ended at the 113 s, lasting for 58 s and releasing approximately 56.5% of the total moment; The third stage started at the 122 s and ended at the 142 s, lasting for 20 s and releasing approximately 19.1% of the total moment. The length of the ruptured fault plane is about 490 km. The maximum width of the ruptured fault plane is about 45 km. The rupture mainly occurred within 30 km in depth under the surface of the Earth. The average static slip in the underground rocky crust is about 1.2 m with the maximum static slip 3.6 m. The average static stress drop is about 5 MPa with the maximum static stress drop 18 MPa. The maximum static slip and the maximum stress drop occurred in an area within 50 km to the east of the centroid position. Foundation item: Joint Seismological Science Foundation of China (103066) and Foundation of the Seismic Pattern and Digital Seismic Data Application Research Office of Institute of Earthquake Science of the China Earthquake Administration.     

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.     

Evaporation from the understorey in the Jarrah (Eucalyptus marginata Don ex Sm.) forest, southwestern Australia

Annual evaporation from groundflora, litter and soil of the jarrah forest was estimated from measurements of daily evaporation by ventilated chambers on several days over two separate 12-month periods. In the first year, when sampling ranged over 0.1 ha of forest, annual evaporation during daylight hours was estimated as 410 mm (0.32 rainfall). In the second year, sampling was more frequent, on a larger scale, and included the night hours. Annual evaporation was estimated at 360 mm (0.36 rainfall).

Similarly, in the second year, annual evaporation from two trees of the dominant middle storey species, Banksia grandis, was estimated at 7500 and 18,9001 respectively. The leaf area of these two trees was 9.6 and 22.4 m², respectively, so that annual evaporation, when expressed as mm³ per mm² leaf area, was similar for both trees (mean = 820 ± 30 mm). Applying that value to all Banksia trees in a hectare of forest, and using a measured estimate of leaf area index of 0.19, the estimated annual evaporation from the Banksia component was 155 mm (0.16 rainfall). For the upland part of the forest sampled, the combined annual evaporation from the lower and middle storeys accounted for about half (0.51) of the annual rainfall.

We conclude that reduced evaporation from the upper storey following clearing or thinning may be strongly counteracted by increased evaporation from the understorey due to increased availability of energy and water.     

Assessment and analysis of microwave emissivity and transmissivity of a deciduous forest towards the estimate of vegetation biomass

Forests play an important role in the global carbon cycle and have a potential impact on global climatic change.Monitoring forest biomass is of considerable importance in understanding the hydrological cycle.Because of the problem of dense forest cover,no reliable method with which to retrieve soil moisture in forest areas from the microwave emission signature has been established.All of these issues relate to the microwave emissivity and transmissivity characteristics of a forest.The microwave emission contribution received by a sensor above a forest canopy comes from both the soil surface and the vegetation layer.To analyze the relationship of forest biomass and forest emission and transmissivity,a high-order emission model,the matrix-doubling model,which consists of both soil and vegetation models,was developed and then validated for a young deciduous forest stand in a field experiment.To simulate the emissivity and transmissivity of a deciduous forest in the L and X bands using the matrix-doubling model,the parameters of components of deciduous trees when the leaf area index varies from 1 to10 were generated by an L-system and a forest growth model.The emissivity and transmissivity of a forest and the relationships of these parameters to forest biomass are presented and analyzed in this paper.Emissivity in the L band when the leaf area index is less than 6 and at viewing angles less than 40°,and transmissivity in the L band are the most sensitive parameters in deciduous forest biomass estimation.     

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.