Data‐driven modelling of shortwave radiation transfer to snow through boreal birch and conifer canopies

Hydrological and land surface models require simple but accurate methods to predict the solar radiation transmitted through vegetation to snow, backed up by direct comparisons to data. Twenty shortwave pyranometers were deployed in forest plots of varying canopy structures and densities in sparse birch forest near Abisko, Sweden, in spring 2011 and mixed conifer forest near Sodankylä, Finland, in spring 2012. Above‐canopy global and diffuse shortwave irradiances were also measured. These data were used to test a model that uses hemispherical photographs to explicitly estimate both diffuse radiation and direct beam transmission, as well as two models that apply bulk canopy parameters and versions of Beers Law. All three models predict canopy shortwave transmission similarly well for leafless birch forest, but for conifers, the bulk methods perform poorly. In addition, an existing model of multiple reflections between canopy and snow was found to be suitable for birch, but not conifers. A new bulk approach based on empirical relationships with hemisphere‐averaged sky view fraction showed improved performance for both sites; this suggests benefits of avoiding the use of plant area index calculated from optical methods, which can introduce errors. Furthermore, tests using common empirical diffuse radiation models were shown to underestimate shortwave transmission by up to 7% relative to using the data, suggesting that new diffuse models are required for high latitudes. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Influence of canopy density on snow distribution in a temperate mountain range

We analyse spatial variability and different evolution patterns of snowpack in a mixed beech–fir stand in the central Pyrenees. Snow depth and density were surveyed weekly along six transects of contrasting forest cover during a complete accumulation and melting season; we also surveyed a sector unaffected by canopy cover. Forest density was measured using the sky view factor (SVF) obtained from digital hemispherical photographs. During periods of snow accumulation and melting, noticeable differences in snow depth and density were found between the open site and those areas covered by forest canopy. Principal component analysis provided valuable information in explaining these observations. The results indicate a high variability in snow accumulation within forest areas related to differences in canopy density. Maximum snow water equivalent (SWE) was reduced by more than 50% beneath dense canopies compared with clearings, and this difference increased during the melting period. We also found significant temporal variations: when melting began in sectors with low SVF, most of the snow had already thawed in areas with high SVF. However, specific conditions occasionally produced a different response of SWE to forest cover, with lower melting rates observed beneath dense canopies. The high values of correlation coefficients for SWE and SVF (r > 0·9) indicate the reliability of predicting the spatial distribution of SWE in forests when only a moderate number of observations are available. Digital hemispherical photographs provide an appropriate tool for this type of analysis, especially for zenith angles in the range 35–55 . Copyright © 2007 John Wiley & Sons, Ltd.     

Energy balance above a boreal coniferous forest: a difference in turbulent fluxes between snow‐covered and snow‐free canopies

To evaluate the interactive effects of snow and forest on turbulent fluxes between the forest surface and the atmosphere, the surface energy balance above a forest was measured by the eddy correlation method during the winter of 1995–1996. The forest was a young coniferous plantation comprised of spruce and fir. The study site, in Sapporo, northern Japan, had heavy and frequent snowfalls and the canopy was frequently covered with snow during the study period. A comparison of the observed energy balance above the forest for periods with and without a snow‐covered canopy and an analysis using a single‐source model gave the following results: during daytime when the canopy was covered with snow, the upward latent heat flux was large, about 80% of the net radiation, and the sensible heat flux was positive but small. On the other hand, during daytime when the canopy was dry and free from snow, the sensible heat flux was dominant and the latent heat flux was minor, about 10% of the net radiation. To explain this difference of energy partition between snow‐covered and snow‐free conditions, not only differences in temperature but also differences in the bulk transfer coefficients for latent heat flux were necessary in the model. Therefore, the high evaporation rate from the snow‐covered canopy can be attributed largely to the high moisture availability of the canopy surface. Evaporation from the forest during a 60‐day period in midwinter was estimated on a daily basis as net radiation minus sensible heat flux. The overall average evaporation during the 60‐day period was 0·6 mm day⁻¹, which is larger than that from open snow fields. Copyright © 1999 John Wiley & Sons, Ltd.     

Observations of urban snow properties in Calgary,Canada

Urban winter hydrology has garnered very little attention owing to the general notion that high‐intensity rainfalls are the major flood‐generating events in urban areas. As a result, few efforts have been made to research urban snow and its melt characteristics. This study investigates the characteristics of urban snow that differentiate it from rural snow, and makes recommendations for incorporating these characteristics into an urban snowmelt model. A field study was conducted from the fall of 2001 to the spring of 2002 in the city of Calgary, Canada. Snow depths and densities, soil moisture, soil temperature, snow albedo, net radiation, snow evaporation, and surface temperature were measured at several locations throughout the winter period. The combination of urban snow removal practices and the physical elements that exist in urban areas were found to influence the energy balance of the snowpack profoundly. Shortwave radiation was found to be the main source of energy for urban snow; as a consequence, the albedo of urban snow is a very important factor in urban snowmelt modelling. General observations lead to the classification of snow as one of four types: snow piles, snow on road shoulders, snow on sidewalk edges, and snow in open areas. This resulted in the development of four separate functions for the changing snow albedo values. A study of the frozen ground conditions revealed that antecedent soil moisture conditions had very little impact on frozen ground, and thus frozen ground very nearly always acts as a near impervious area. Improved flood forecasting for urban catchments in cold regions can only be achieved with accurate modelling of urban winter runoff that involves the energy balance method, incorporating snow redistribution and urban snow‐cover characteristics, and using small time steps. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Modelling snow water equivalent and spring runoff in a boreal watershed,James Bay,Canada

The hydrology of boreal regions is strongly influenced by seasonal snow accumulation and melt. In this study, we compare simulations of snow water equivalent (SWE) and streamflow by using the hydrological model HYDROTEL with two contrasting approaches for snow modelling: a mixed degree‐day/energy balance model (small number of inputs, but several calibration parameters needed) and the thermodynamic model CROCUS (large number of inputs, but no calibration parameter needed). The study site, in Northern Quebec, Canada was equipped with a ground‐based gamma ray sensor measuring the SWE continuously for 5 years in a small forest clearing. The first simulation of CROCUS showed a tendency to underestimate SWE, attributable to bias in the meteorological inputs. We found that it was appropriate to use a threshold of 2 °C to separate rain and snow. We also applied a correction to account for snowfall undercatch by the precipitation gauge. After these modifications to the input dataset, we noticed that CROCUS clearly overestimated the SWE, likely as a result of not including loss in SWE because of blowing snow sublimation and relocation. To correct this, we included into CROCUS a simple parameterisation effective after a certain wind speed threshold, after which the thermodynamic model performed much better than the traditional mixed degree‐day/energy balance model. HYDROTEL was then used to simulate streamflow with both snow models. With CROCUS, the main peak flow could be captured, but the second peak because of delayed snowmelt from forested areas could not be reproduced due to a lack of sub‐canopy radiation data to feed CROCUS. Despite the relative homogeneity of the boreal landscape, data inputs from each land cover type are needed to generate satisfying simulation of the spring runoff. Copyright © 2013 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.     

The effects of forest litter on snow energy budget in the Tianshan Mountains,China

Forest litter exerts an impact on the energy budget of snow surfaces, which lie beneath forest canopies. In this study, we measured shortwave and longwave radiation levels, as well as quantities of Asian spruce (Picea schrenkinan ) forest litter, over 3 snow study plots that representing an open environment, 20% forest canopy openness (20% FCO), and 80% forest canopy openness (80% FCO). The fractional litter coverage (lc ) was obtained through the binarization of digital photographs of forest litter. The effects of forest litter on snow surface albedo (α ), snow surface temperature (T _s ), upward shortwave and longwave radiation (K _↑ and L _↑), and sensible heat flux (H ) were then analyzed. According to our results, the energy budget over snow surface influenced by forest litter principally due to forest litter forcing α decrease and T _s increase. The effects of forest litter on the energy budget increased with time and lc . We found that forest litter exerted the most significant impact on K _↑ and L _↑ at daytime during the latter stages of the snowmelt period. The influence of forest litter on H was more apparent on windy days. The presence of forest litter increased gains in shortwave radiation and losses in longwave radiation and decreased gains in H . Compared to the simulated energy (K _↑ + L _↑ + H ) over a snow surface without litter, the calculated energy decreased by ?13.4 W/m² and increased by 9.0 W/m², respectively, at the 20% FCO and 80% FCO sites during the latter stages of the snowmelt period. Overall, forest litter facilitated snow surface energy gains at the 80% FCO site and impeded them at the 20% FCO site during the latter stages of the snowmelt period.     

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.     

Exploring the spatiotemporal variability of the snow water equivalent in a small boreal forest catchment through observation and modelling

In snow-fed catchments, it is crucial to monitor and model the snow water equivalent (SWE), particularly when simulating the melt water runoff. SWE distribution can, however, be highly heterogeneous, particularly in forested environments. Within these locations, scant studies have explored the spatiotemporal variability in SWE in relation with vegetation characteristics, with only few successful attempts. The aim of this paper is to fill this knowledge gap, through a detailed monitoring at nine locations within a 3.49 km² forested catchment in southern Québec, Canada (47°N, 71°W). The catchment receives an annual average of 633 mm of solid precipitation and is predominantly covered with balsam fir stands. Extracted from intensive field campaign and high-resolution LiDAR data, this study explores the effect of fine scale forest features (tree height, tree diameter, canopy density, leaf area index [LAI], tree density and gap fraction) on the spatiotemporal variability in the SWE distribution. A nested stratified random sampling design was adopted to quantify small-scale variability across the catchment and 1810 manual snow samples were collected throughout the consecutive winters of 2016–17 and 2017–18. This study explored the variability of SWE using coefficients of variation (CV) and relating to the LAI. We also present existing spatiotemporal differences in maximum snow depth across different stands and its relationship with average tree diameter. Furthermore, exploiting key vegetation characteristics, this paper explores different approaches to model SWE, such as multiple linear regression, binary regression tree and neural networks (NN). We were unable to establish any relationship between the CV of SWE and the LAI. However, we observed an increase in maximum snow depth with decreasing tree diameter, suggesting an association between these variables. NN modelling (Nash-Sutcliffe efficiency [NSE] = 0.71) revealed that, snow depth, snowpack age and forest characteristics (tree diameter and tree density) are key controlling variables on SWE. Using only variables that are deemed to be more readily available (snow depth, tree height, snowpack age and elevation), NN performance falls by only 7% (NSE = 0.66).     

Estimating sublimation of intercepted and sub‐canopy snow using eddy covariance systems

Direct measurements of winter water loss due to sublimation were made in a sub‐alpine forest in the Rocky Mountains of Colorado. Above‐and below‐canopy eddy covariance systems indicated substantial losses of winter‐season snow accumulation in the form of snowpack (0·41 mm d^?1) and intercepted snow (0·71 mm d^?1) sublimation. The partitioning between these over and under story components of water loss was highly dependent on atmospheric conditions and near‐surface conditions at and below the snow/atmosphere interface. High above‐canopy sensible heat fluxes lead to strong temperature gradients between vegetation and the snow‐surface, driving substantial specific humidity gradients at the snow surface and high sublimation rates. Intercepted snowfall resulted in rapid response of above‐canopy latent heat fluxes, high within‐canopy sublimation rates (maximum = 3·7 mm d^?1), and diminished sub‐canopy snowpack sublimation. These results indicate that sublimation losses from the sub‐canopy snowpack are strongly dependent on the partitioning of sensible and latent heat fluxes in the canopy. This compels comprehensive studies of snow sublimation in forested regions that integrate sub‐canopy and over‐story processes. Copyright © 2007 John Wiley & Sons, Ltd.     

Diagnosing snow accumulation errors in a rain‐snow transitional environment with snow board observations

Diagnosing the source of errors in snow models requires intensive observations, a flexible model framework to test competing hypotheses, and a methodology to systematically test the dominant snow processes. We present a novel process‐based approach to diagnose model errors through an example that focuses on snow accumulation processes (precipitation partitioning, new snow density, and snow compaction). Twelve years of meteorological and snow board measurements were used to identify the main source of model error on each snow accumulation day. Results show that modeled values of new snow density were outside observational uncertainties in 52% of days available for evaluation, while precipitation partitioning and compaction were in error 45% and 16% of the time, respectively. Precipitation partitioning errors mattered more for total winter accumulation during the anomalously warm winter of 2014–2015, when a higher fraction of precipitation fell within the temperature range where partition methods had the largest error. These results demonstrate how isolating individual model processes can identify the primary source(s) of model error, which helps prioritize future research.     

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.     

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.     

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.     

Automated system for measuring snow surface energy balance components in mountainous terrain

The ability to continually monitor several meteorological parameters is needed to estimate snow surface energy balance components in mountainous terrain. In remote mountainous locations, limited accessibility and extreme weather conditions limit the use of delicate meteorological instrumentation. Robust instrumentation and radio telemetry are often needed to measure snow surface energy exchanges. This study examined the practicality and effectiveness of robust instrumentation in estimating radiative and turbulent exchanges in the forested Bear River Mountains of northern Utah. Measurement of reflected shortwave radiation was problematic due to possible selective absorption in the infra-red range. This resulted in overestimates of reflected shortwave radiation and decreased estimates of now surface albedo. During high snowfall, the pyranometer and net radiometer were occasionally covered with snow, resulting in inaccurate radiation measurements. Snow typically melted from instrument surfaces in less than one day under full sun. A relative humidity measurement accuracy of ± 4% may have resulted in a possible error of 20% in the calculation of vapour pressure. Snow depth measurement with an acoustical sensor was affected by new or blowing snow, which resulted in inaccurate snow depth measurement 16.2% of the time. The longest period without a valid snow depth measurement was 19.5 hours. A new snow temperature thermocouple ladder was designed and constructed and provided accurate within-pack temperature measurements throughout the pre-melt and melt season.     

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.     

Simulated watershed responses to land cover changes using the Regional Hydro‐Ecological Simulation System

In this work, we used the Regional Hydro‐Ecological Simulation System (RHESSys) model to examine runoff sensitivity to land cover changes in a mountain environment. Two independent experiments were evaluated where we conducted simulations with multiple vegetation cover changes that include conversion to grass, no vegetation cover and deciduous/coniferous cover scenarios. The model experiments were performed at two hillslopes within the Weber River near Oakley, Utah watershed (USGS gauge # 10128500). Daily precipitation, air temperature and wind speed data as well as spatial data that include a digital elevation model with 30 m grid resolution, soil texture map and vegetation and land use maps were processed to drive RHESSys simulations. Observed runoff data at the watershed outlet were used for calibration and verification. Our runoff sensitivity results suggest that during winter, reduced leaf area index (LAI) decreases canopy interception resulting in increased snow accumulations and hence snow available for runoff during the early spring melt season. Increased LAI during the spring melt season tends to delay the snow melting process. This delay in snow melting process is due to reduced radiation beneath high LAI surfaces relative to low LAI surfaces. The model results suggest that annual runoff yield after removing deciduous vegetation is on average about 7% higher than with deciduous vegetation cover, while annual runoff yield after removing coniferous vegetation is on average as about 2% higher than that produced with coniferous vegetation cover. These simulations thus help quantify the sensitivity of water yield to vegetation change. 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.