Modelling forest snow processes with a new version of WaSiM

ABSTRACT

We present a new model extension for the Water balance Simulation Model, WaSiM, which features (i) snow interception and (ii) modified meteorological conditions under coniferous forest canopies, complementing recently developed model extensions for particular mountain hydrological processes. Two study areas in Austria and Germany are considered in this study. To supplement and constrain the modelling experiments with on-site observations, a network of terrestrial time-lapse cameras was set up in one of these catchments. The spatiotemporal patterns of snow depth inside the forest and at the adjacent open field sites were recorded along with snow interception dynamics. Comparison of observed and modelled snow cover and canopy interception indicates that the new version of WaSiM reliably reconstructs the variability of snow accumulation for both the forest and the open field. The Nash-Sutcliffe efficiency computed for selected runoff events in spring increases from ?0.68 to 0.71 and 0.21 to 0.87, respectively.     

The simulation of heat and water exchange in the boreal spruce forest by the land-surface model SWAP

All previous versions of a physically based land-surface model SWAP have assumed for simplicity that vegetation is fully covered by snow during the cold season. Such assumption is reasonable only for the regions dominated by short vegetation or for warm climates where snow processes are absent. The major goals of this paper are (i) modification of the latest version of SWAP by incorporation of tall vegetation into the cold-season parameterizations to make the model applicable for simulating heat and water transfer within a boreal forest biome and (ii) validation of the modified version using the data from a forested catchment located in the boreal zone. Modification of SWAP required to parameterize radiative and turbulent exchange between the forest crown and forest floor, partitioning snowfall between interception by the canopy (in doing so, snow interception differs from rain interception) and falling to the ground, formation of snow cover on the forest crown and forest floor including snow accumulation (both in solid and liquid fractions), snow evaporation, and snowmelt. The advanced model was validated using a set of hydrometeorological data measured during 18 years (1966–1983) at the Tayozhniy catchment (covered by boreal spruce forest), Valdai, Russia. Simulations of annual and monthly snow/rain interception, daily runoff at the catchment outlet, snow density, snow depth, snow water equivalent, soil water storage in three layers (0–20, 0–50 and 0–100 cm), and monthly evapotranspiration from the catchment were compared with observations. Analysis of the results of validation shows that the new version of the model SWAP reproduces the heat and water exchange processes occurring in mid-latitude boreal forest quite reasonable.     

WINTER RADIATION EXTINCTION AND REFLECTION IN A BOREAL PINE CANOPY: MEASUREMENTS AND MODELLING

Predicting the rate of snowmelt and intercepted snow sublimation in boreal forests requires an understanding of the effects of snow-covered conifers on the exchange of radiant energy. This study examined the amount of intercepted snow on a jack pine canopy in the boreal forest of central Saskatchewan and the shortwave and net radiation exchange with this canopy, to determine the effect of intercepted snow and canopy structure on shortwave radiation reflection and extinction and net radiation attenuation in a boreal forest. The study focused on clear sky conditions, which are common during winter in the continental boreal forest. Intercepted snow was found to have no influence on the clear-sky albedo of the canopy, the extinction of short wave radiation by the canopy or ratio of net radiation at the canopy top to that at the surface snow cover. Because of the low albedo of the snow-covered canopy, net radiation at the canopy top remains positive and a large potential source of energy for sublimation. The canopy albedo declines somewhat as the extinction efficiency of the underlying canopy increases. The extinction efficiency of short wave radiation in the canopy depends on solar angle because of the approximately horizontal orientation of pine branches. For low solar angles above the horizon, the extinction efficiency is quite low and short wave transmissivity through the canopy is relatively high. As the solar angle increases, extinction increases up to angles of about 50°, and then declines. Extinction of short wave radiation in the canopy strongly influences the attenuation of net radiation by the canopy. Short wave radiation that is extinguished by branches is radiated as long wave, partly downwards to the snow cover. The ratio of net radiation at the canopy top to that at the snow cover surface increases with the extinction of short wave radiation and is negative for low extinction efficiencies. For the pine canopy examined, the daily mean net radiation at the snow cover surface became positive when daily mean solar angles exceeded 22° in late March. Hence, canopy structure and solar angle control the net radiation at the snow cover surface during clear sky conditions and will govern the timing and rate of snowmelt. Models of intercepted snow sublimation and forest snowmelt could beneficially incorporate the canopy radiation balance, which can be extrapolated to stands of various canopy densities, coverage and heights in a physically based manner. Such models could hence avoid ‘empirical’ temperature index measures that cannot be extrapolated with confidence.     

Assessing the effects of post‐pine beetle forest litter on snow albedo

The effect of forest litter on snow surface albedo has been subject to limited study, mainly in the hardwood‐dominated forests of the northeastern United States. Given the recent pine beetle infestation in Western North America and associated increases in litter production, this study examines the effects of forest litter on snow surface albedo in the coniferous forests of south‐central British Columbia. Measured changes in canopy transmittance provide an indication of canopy loss or total litterfall over the winter of 2007–2008. Relationships between percent litter cover, an index of albedo, snow depth, and snow ablation during the 2008 melt season are compared between a mature, young, and clearcut coniferous stand. Results indicate a strong feedback effect between canopy loss and subsequent enhanced shortwave transmittance, and litter accumulation on the snow surface from that canopy loss. However, this relationship is confounded by other variables concurrently affecting albedo. While results suggest that a relatively small percent litter cover can have a significant effect on albedo and ablation, further research is underway to extract the litter signal from that of other factors affecting albedo, particularly snow depth. Copyright © 2010 John Wiley & Sons, Ltd.     

Representation of canopy snow interception,unloading and melt in a parsimonious snowmelt model

Recent improvements in the Utah Energy Balance (UEB) snowmelt model are focused on snow–vegetation–atmosphere interactions to understand how different types of vegetation affect snow processes in the mountains of Western USA. This work presents field work carried out in the Rocky Mountains of Northern Utah to evaluate new UEB model algorithms that represent the processes of canopy snow interception, sublimation, mass unloading and melt. Four years' continuous field observations showed generally smaller accumulations of snow beneath the forest canopies in comparison with open (sage and grass) areas, a difference that is attributed to interception and subsequent sublimation and redistribution of intercepted snow by wind, much of it into surrounding open areas. Accumulations beneath the denser forest (conifer) canopies were found to be less than the accumulation beneath the less dense forest (deciduous) canopies. The model was able to represent the accumulation of snow water equivalent in the open and beneath the deciduous forest quite well but without accounting for redistribution tended to overestimate the snow water equivalent beneath the conifer forest. Evidence of redistribution of the intercepted snow from the dense forest (i.e. conifer forest) to the adjacent area was inferred from observations. Including a simple representation of redistribution in the model gave satisfactory prediction of snow water equivalent beneath the coniferous forest. The simulated values of interception, sublimation and unloading were also compared with previous studies and found in agreement. Copyright © 2013 John Wiley & Sons, Ltd.     

Modified treatment of intercepted snow improves the simulated forest albedo in the Canadian Land Surface Scheme

The Canadian Land Surface Scheme (CLASS) was modified to correct an underestimation of the winter albedo in evergreen needleleaf forests. Default values for the visible and near‐infrared albedo of a canopy with intercepted snow, α_VIS,cs and α_NIR,cs, respectively, were too small, and the fraction of the canopy covered with snow, f_snow, increased too slowly with interception, producing a damped albedo response. A new model for f_snow is based on z_I*, the effective depth of newly intercepted snow required to increase the canopy albedo to its maximum, which corresponds in the model with f_snow = 1. Snow unloading rates were extracted from visual assessments of photographs and modelled based on relationships with meteorological variables, replacing the time‐based method employed in CLASS. These parameterizations were tested in CLASS version 3.6 at boreal black spruce and jack pine forests in Saskatchewan, Canada, a subalpine Norway spruce and silver fir forest at Alptal, Switzerland, and a boreal maritime forest at Hitsujigaoka, Japan. Model configurations were assessed based on the index of agreement, d, relating simulated and observed daily albedo. The new model employs α_VIS,cs = 0.27, α_NIR,cs = 0.38 and z_I* = 3 cm. The best single‐variable snow unloading algorithm, determined by the average cross‐site d, was based on wind speed. Two model configurations employing ensemble averages of the unloading rate as a function of total incoming radiation and wind speed, and air temperature and wind speed, respectively, produced larger minimum cross‐site d values but a smaller average. The default configuration of CLASS 3.6 produced a cross‐site average d from October to April of 0.58. The best model employing a single parameter (wind speed at the canopy top) for modelling the unloading rate produced an average d of 0.86, while the two‐parameter ensemble‐average unloading models produced a minimum d of 0.81 and an average d of 0.84. © 2015 Her Majesty the Queen in Right of Canada. Hydrological Processes published by John Wiley & Sons, Ltd.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

A conceptual model of water yield effects from beetle‐induced tree death in snow‐dominated lodgepole pine forests

In regions of western North America with snow‐dominated hydrology, the presence of forested watersheds can significantly influence streamflow compared to areas with other vegetation cover types. Widespread tree death in these watersheds can thus dramatically alter many ecohydrologic processes including transpiration, canopy solar transmission and snow interception, subcanopy wind regimes, soil infiltration, forest energy storage and snow surface albedo. One of the more important causes of conifer tree death is bark beetle infestation, which in some instances will kill nearly all of the canopy trees within forest stands. Since 1996, an ongoing outbreak of bark beetles (Coleoptera: Scolytidae) has caused widespread mortality across more than 600,000 km² of coniferous forests in western North America, including numerous Rocky Mountain headwaters catchments with high rates of lodgepole pine (Pinus contorta) mortality from mountain pin beetle (Dendroctonous ponderosae) infestations. Few empirical studies have documented the effects of MPB infestations on hydrologic processes, and little is known about the direction and magnitude of changes in water yield and timing of runoff due to insect‐induced tree death. Here, we review and synthesize existing research and provide new results quantifying the effects of beetle infestations on canopy structure, snow interception and transmission to create a conceptual model of the hydrologic effects of MPB‐induced lodgepole pine death during different stages of mortality. We identify the primary hydrologic processes operating in living forest stands, stands in multiple stages of death and long‐dead stands undergoing regeneration and estimate the direction of change in new water yield. This conceptual model is intended to identify avenues for future research efforts. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimation of winter leaf area index and sky view fraction for snow modelling in boreal coniferous forests: consequences on snow mass and energy balance

Leaf area index (LAI) and canopy coverage are important parameters when modelling snow process in coniferous forests, controlling interception and transmitting radiation. Estimates of LAI and sky view factor show large variability depending on the estimation method used, and it is not clear how this is reflected in the calculated snow processes beneath the canopy. In this study, the winter LAI and sky view fraction were estimated using different optical and biomass‐based approximations in several boreal coniferous forest stands in Fennoscandia with different stand density, age and site latitude. The biomass‐based estimate of LAI derived from forest inventory data was close to the values derived from the optical measurements at most sites, suggesting that forest inventory data can be used as input to snow hydrological modelling. Heterogeneity of tree species and site fertility, as well as edge effects between different forest compartments, caused differences in the LAI estimates at some sites. A snow energy and mass balance model (SNOWPACK) was applied to detect how the differences in the estimated values of the winter LAI and sky view fraction were reflected in simulated snow processes. In the simulations, an increase in LAI and a decrease in sky view fraction changed the snow surface energy balance by decreasing shortwave radiation input and increasing longwave radiation input. Changes in LAI and sky view fraction affected directly snow accumulation through altered throughfall fraction and indirectly snowmelt through the changed surface energy balance. Changes in LAI and sky view fraction had a greater impact on mean incoming radiation beneath the canopy than on other energy fluxes. Snowmelt was affected more than snow accumulation. The effect of canopy parameters on evaporation loss from intercepted snow was comparable with the effect of variation in governing meteorological variables such as precipitation intensity and air temperature. Copyright © 2012 John Wiley & Sons, Ltd.     

Seasonal snowpack dynamics and runoff in a cool temperate forest: lysimeter experiment in Niigata,Japan

Seasonal snowpack dynamics are described through field measurements under contrasting canopy conditions for a mountainous catchment in the Japan Sea region. Microclimatic data, snow accumulation, albedo and lysimeter runoff are given through the complete winter season 2002–03 in (1) a mature cedar stand, (2) a larch stand, and (3) a regenerating cedar stand or opening. The accumulation and melt of seasonal snowpack strongly influences streamflow runoff during December to May, including winter baseflow, mid‐winter melt, rain on snow, and diurnal peaks driven by radiation melt in spring. Lysimeter runoff at all sites is characterized by constant ground melt of 0·8–1·0 mm day⁻¹. Rapid response to mid‐winter melt or rainfall shows that the snowpack remains in a ripe or near‐ripe condition throughout the snow‐cover season. Hourly and daily lysimeter discharge was greatest during rain on snow (e.g. 7 mm h⁻¹ and 53 mm day⁻¹ on 17 December) with the majority of runoff due to rainfall passing through the snowpack as opposed to snowmelt. For both rain‐on‐snow and radiation melt events lysimeter discharge was generally greatest at the open site, although there were exceptions such as during interception melt events. During radiation melt instantaneous discharge was up to 4·0 times greater in the opening compared with the mature cedar, and 48 h discharge was up to 2·5 times greater. Perhaps characteristic of maritime climates, forest interception melt is shown to be important in addition to sublimation in reducing snow accumulation beneath dense canopies. While sublimation represents a loss from the catchment water balance, interception melt percolates through the snowpack and contributes to soil moisture during the winter season. Strong differences in microclimate and snowpack albedo persisted between cedar, larch and open sites, and it is suggested further work is needed to account for this in hydrological simulation models. Copyright © 2005 John Wiley & Sons, Ltd.     

Insights on stream temperature processes through development of a coupled hydrologic and stream temperature model for forested coastal headwater catchments

Stream temperature controls a number of biological, chemical, and physical processes occurring in aquatic environments. Transient snow cover and advection associated with lateral throughflow inputs can have a dominant influence on stream thermal regimes for headwater catchments in the rain‐on‐snow zone. Most existing stream temperature models lack the ability to properly simulate these processes. We developed and evaluated a conceptual‐parametric catchment‐scale stream temperature model that includes the role of transient snow cover and lateral advection associated with throughflow. The model consists of routines for simulating canopy interception, snow accumulation and melt, hillslope throughflow runoff and temperature, and stream channel energy exchange processes. The model was used to predict discharge and stream temperature for a small forested headwater catchment near Vancouver, Canada, using long‐term (1963–2013) weather data to compute model forcing variables. The model was evaluated against 4 years of observed stream temperature. The model generally predicted daily mean stream temperature accurately (annual RMSE between 0.57 and 1.24 °C) although it overpredicted daily summer stream temperatures by up to 3 °C during extended low streamflow conditions. Model development and testing provided insights on the roles of advection associated with lateral throughflow, channel interception of snow, and surface–subsurface water interactions on stream thermal regimes. This study shows that a relatively simple but process‐based model can provide reasonable stream temperature predictions for forested headwater catchments located in the rain‐on‐snow zone.     

Increasing water losses from snow captured in the canopy of boreal forests: A case study using a 30 year data set

Water losses from snow intercepted by forest canopy can significantly influence the hydrological cycle in seasonally snow‐covered regions, yet how snow interception losses (SIL) are influenced by a changing climate are poorly understood. In this study, we used a unique 30 year record (1986–2015) of snow accumulation and snow water equivalent measurements in a mature mixed coniferous (Picea abies and Pinus sylvestris ) forest stand and an adjacent open area to assess how changes in weather conditions influence SIL. Given little change in canopy cover during this study, the 20% increase in SIL was likely the result of changes in winter weather conditions. However, there was no significant change in average wintertime precipitation and temperature during the study period. Instead, mean monthly temperature values increased during the early winter months (i.e., November and December), whereas there was a significant decrease in precipitation in March. We also assessed how daily variation in meteorological variables influenced SIL and found that about 50% of the variation in SIL was correlated to the amount of precipitation that occurred when temperatures were lower than ?3 °C and to the proportion of days with mean daily temperatures higher than +0.4 °C. Taken together, this study highlights the importance of understanding the appropriate time scale and thresholds in which weather conditions influence SIL in order to better predict how projected climate change will influence snow accumulation and hydrology in boreal forests in the future.     

Snowfall interception in a deciduous Nothofagus forest and implications for spatial snowpack distribution

Native Nothofagus forests in the midlatitude region of the Andes Cordillera are notorious biodiversity hot spots, uniquely situated in the Southern Hemisphere such that they develop in snow‐dominated reaches of this mountain range. Spanning a smaller surface area than similar ecosystems, where forests and snow coexist in the Northern Hemisphere, the interaction between vegetation and snow processes in this ecotone has received lesser attention. We present the first systematic study of snow–vegetation interactions in the Nothofagus forests of the Southern Andes, focusing on how the interplay between interception and climate determines patterns of snow water equivalent (SWE) variability. The Valle Hermoso experimental catchment, located in the Nevados de Chillán vicinity, was fitted with eight snow depth sensors that provided continuous measurements at varying elevations, aspect, and forest cover. Also, manual measurements of snow properties were obtained during snow surveys conducted during end of winter and spring seasons for 3 years, between 2015 and 2017. Each year was characterized by distinct climatological conditions, with 2016 representing one of the driest winters on record in this region. Distance to canopy, leaf area index, and total gap area were measured at each observational site. A regression model was built on the basis of statistical analysis of local parameters to model snow interception in this kind of forest. We find that interception implied a 23.2% reduction in snow accumulation in forested sites compared with clearings. The interception in these deciduous trees represents, on average, 23.6% of total annual snowfall, reaching a maximum measured interception value of 13.8‐mm SWE for all snowfall events analysed in this research.     

Snow process monitoring in montane forests with time‐lapse photography

A camera network with hourly resolution was used to monitor the complex snow processes in montane forest environments. We developed a semi‐automatic procedure to interpret snow depths from the digital images, which exhibited high consistency with manual measurements and station‐based recordings. To extract snow interception dynamics, six binary classification methods were compared. The MaxEntropy classifier demonstrated better performance than the other methods under conditions of varying illumination and was therefore selected as the method used for quantifying snow in tree canopies. Snow accumulation and ablation on the ground, as well as snow loading and unloading in the forest canopies, were investigated using snow parameters derived from the time‐lapse photography monitoring. The influences of meteorologic conditions, forest cover, and elevation on the snow processes were also considered. Time‐lapse photography proved to be an effective and low‐cost approach for collecting useful information on snow processes and facilitating the set‐up of hydrological models.     

Estimating winter evaporation in boreal forests with operational snow course data

Snow course measurements from 11 sites located in eastern and northern Finland were used to estimate the total interception evaporation of a winter season for different forest categories. We categorized the sites based on forest density and tree species. Results showed that interception loss from gross precipitation increased with forest density and approached 30% for a forest with the highest density class. Interception loss for the most common forest density class was 11%. Interception losses were slightly larger in spruce forests than in pine, deciduous, or mixed forests. We provide suggestions as to how to design snow surveys to estimate wintertime interception evaporation better. Rough terrain and transition zones between forest and open areas should be avoided. Since evaporation fraction was strongly dependent on tree crown characteristics, snow course data should include direct estimates of canopy closure. Qualitative observations made by different observers should be given a reference frame to ensure comparability of records from different sites. Copyright © 2003 John Wiley & Sons, Ltd.     

Potential of time‐lapse photography of snow for hydrological purposes at the small catchment scale

Time‐lapse photography provides an attractive source of information about snow cover characteristics, especially at the small catchment scale. The objective of this study was to design and test a monitoring system, which allows multi‐resolution observations of snow cover characteristics. The main aim was to simultaneously investigate the spatio‐temporal patterns of snow cover, snow depth and snowfall interception in the area very close to the camera, and the spatio‐temporal patterns of snow cover in the far range. The multi‐resolution design was tested at three sites in the eastern part of the Austrian Alps (Hochschwab‐Rax region). Digital photographs were taken at hourly time steps between 6:00 and 18:00 in the period November, 2004 to December, 2006. The results showed that the time‐lapse photography allows effective mapping of the snow depths at high temporal resolution in the region close to the digital camera at many snow stake locations. It is possible to process a large number of photos by using an automatic procedure for accurate snow depth readings. The digital photographs can also be used to infer the settling characteristics of the snow pack and snow interception during the day. Although it is not possible to directly estimate the snow interception mass, the photos may indeed give very useful information on the snow processes on and beneath the forest canopy. The main advantage of using time‐lapse photography in the far range of the digital camera is to observe the spatio‐temporal patterns of snow cover over different landscape configurations. The results illustrate that digital photographs can be very useful for parameterising processes such as sloughing on steep slopes, snow deposition in gullies and snow erosion on mountain ridges in a distributed snow model. Copyright © 2011 John Wiley & Sons, Ltd.     

Charred forests accelerate snow albedo decay: parameterizing the post‐fire radiative forcing on snow for three years following fire

As large, high‐severity forest fires increase and snowpacks become more vulnerable to climate change across the western USA, it is important to understand post‐fire disturbance impacts on snow hydrology. Here, we examine, quantify, parameterize, model, and assess the post‐fire radiative forcing effects on snow to improve hydrologic modelling of snow‐dominated watersheds having experienced severe forest fires. Following a 2011 high‐severity forest fire in the Oregon Cascades, we measured snow albedo, monitored snow, and micrometeorological conditions, sampled snow surface debris, and modelled snowpack energy and mass balance in adjacent burned forest (BF) and unburned forest sites. For three winters following the fire, charred debris in the BF reduced snow albedo, accelerated snow albedo decay, and increased snowmelt rates thereby advancing the date of snow disappearance compared with the unburned forest. We demonstrate a new parameterization of post‐fire snow albedo as a function of days‐since‐snowfall and net snowpack energy balance using an empirically based exponential decay function. Incorporating our new post‐fire snow albedo decay parameterization in a spatially distributed energy and mass balance snow model, we show significantly improved predictions of snow cover duration and spatial variability of snow water equivalent across the BF, particularly during the late snowmelt period. Field measurements, snow model results, and remote sensing data demonstrate that charred forests increase the radiative forcing to snow and advance the timing of snow disappearance for several years following fire. Copyright © 2016 John Wiley & Sons, Ltd.     

MID-WINTER STEMFLOW DRAINAGE FROM BIGTOOTH ASPEN (POPULUS GRANDIDENTATA MICHX.) IN CENTRAL MASSACHUSETTS

Under winter conditions, stemflow drainage in forested ecosystems is often assumed to be a negligible component of the hydrological cycle. This paper reports on mid-winter stemflow drainage from the broadleaved deciduous tree species Populus grandidentata. Stemflow volumes from this species at air temperatures of < 0°C were found to be comparable to rainfall-generated stemflow during summer. Over the three-month period January–March 1993, stemflow ranged from 5.4 to 9.9% of the incident gross precipitation. Expressed as depth equivalents per unit trunk basal area, these stemflow inputs ranged from 1.8 to 4.9 m. These concentrated mid-winter inputs of liquid water to the bases of canopy trees were attributable to: (1) snow interception by the leafless woody frame of each tree; (2) snow retention by glazed ice precipitation associated with the snowfall event; (3) increased temperature at the bark/snow interface caused by the low albedo of the bark tissue; and (4) convergence of snowmelt drainage from steeply inclined upthrust primary branches. The hydrological and ecological significance of liquid water inputs to the forest floor under sub-zero conditions are discussed. © 1997 by John Wiley & Sons Ltd.     

Snowfall into a forest and clearing

Experiments in February–March 1987 at Fraser Experimental Forest in Colorado, U.S.A., compared the number flux of snow particles above the center of an 80-m wide clearing with simultaneous samples in and above the forest 85 m upwind. Numbers of snow particles counted with electronic sensors on towers decreased with decreasing height below the canopy top in the forest and also nearer the surface in the clearing. Average fluxes were greater in the clearing than in the forest, and these differences increased with increasing average wind speed, in the 0–6 m s⁻¹ range, measured during the experiments. The greater flux at the tower in the clearing could not be completely explained by interception of snow in the forest canopy, indicating that aerodynamics may have played a part in creating the difference, even in those light winds. A plumer or concentration of snow particles appeared to have developed near the top of the canopy and spread into the clearing.     

Forest ecosystem processes at the watershed scale: dynamic coupling of distributed hydrology and canopy growth

The hydrological recovery of watersheds from disturbances such as fire and harvest can change the magnitude and distribution of flow paths as the canopy regenerates. The spatial distribution of net water input to the soil–topography system is mediated by vegetation patterns through the processes of interception, evapotranspiration and snowmelt. We have previously described RHESSys, a distributed model of water and carbon flux with a prescribed canopy cover. Although the canopy structure varied spatially it did not change through time. We present an expanded model in which carbon and nitrogen are dynamically coupled with distributed hydrology. The model fixes and allocates canopy carbon annually to reflect changes in climate forcing. We test the interactions of the forest ecosystem to distributed hydrology through controlled experiments. In the first experiment, we prescribe canopy cover and examine the sensitivity of the hydrological outputs to the distribution of vegetation. The canopy distribution is found to have significant effects on simulated hydrological outputs where evaporative demand exceeds available water. In a second experiment we simulate the canopy leaf area index (LAI) across the topography and through time. The model is executed over 100 years using repeated 10-year meteorological records to investigate spatial and temporal patterns of LAI. Annual precipitation and temperature differences result in temporally fluctuating LAI about a reasonably stable long-term mean. The topographical position has a strong effect on local forest canopy characteristics. As expected, simulated ecosystem processes are found to be sensitive to rooting depth in more water limited environments. © 1997 John Wiley & Sons, Ltd.