Variable effects of forest canopies on snow processes in a valley of the central Spanish Pyrenees

Snowpacks and forests have complex interactions throughout the large range of altitudes where they co-occur. However, there are no reliable data on the spatial and temporal interactions of forests with snowpacks, such as those that occur in nearby areas that have different environmental conditions and those that occur during different snow seasons. This study monitored the interactions of forests with snowpacks in four forest stands in a single valley of the central Spanish Pyrenees during three consecutive snow seasons (2015/2016, 2016/2017 and 2017/2018). Daily snow depth data from time-lapse cameras were compared with snow data from field surveys that were performed every 10–15 days. These data thus provided information on the spatial and temporal changes of snow–water equivalent (SWE). The results indicated that forest had the same general effects on snowpack in each forest stand and during each snow season. On average, forest cover reduced the duration of snowpack by 17 days, reduced the cumulative SWE of the snowpack by about 60% and increased the spatial heterogeneity of snowpack by 190%. Overall, forest cover reduced SWE total accumulation by 40% and the rate of SWE accumulation by 25%. The forest-mediated reduction of the accumulation rate, in combination with the occasional forest-mediated enhancement of melting rate, explained the reduced duration of snowpacks beneath forest canopies. However, the magnitude and timing of certain forest effects on snowpack had significant spatial and temporal variations. This variability must be considered when selecting the location of an experimental site in a mountainous area, because the study site should be representative of surrounding areas. The same considerations apply when selecting a time period for study.     

Investigating snowpack across scale in the northern Great Lakes–St. Lawrence forest region of Central Ontario,Canada

This study investigates scaling issues by evaluating snow processes and quantifying bias in snowpack properties across scale in a northern Great Lakes–St. Lawrence forest. Snow depth and density were measured along transects stratified by land cover over the 2015/2016 and 2016/2017 winters. Daily snow depth was measured using a time‐lapse (TL) camera at each transect. Semivariogram analysis of the transect data was conducted, and no autocorrelation was found, indicating little spatial structure along the transects. Pairwise differences in snow depth and snow water equivalent (SWE) between land covers were calculated and compared across scales. Differences in snowpack between forested sites at the TL points corresponded to differences in canopy cover, but this relationship was not evident at the transect scale, indicating a difference in observed process across scale. TL and transect estimates had substantial bias, but consistency in error was observed, which indicates that scaling coefficients may be derived to improve point scale estimates. TL and transect measurements were upscaled to estimate grid scale means. Upscaled estimates were compared and found to be consistent, indicating that appropriately stratified point scale measurements can be used to approximate a grid scale mean when transect data are not available. These findings are important in remote regions such as the study area, where frequent transect data may be difficult to obtain. TL, transect, and upscaled means were compared with modelled depth and SWE. Model comparisons with TL and transect data indicated that bias was dependent on land cover, measurement scale, and seasonality. Modelled means compared well with upscaled estimates, but model SWE was underestimated during spring melt. These findings highlight the importance of understanding the spatial representativeness of in situ measurements and the processes those measurements represent when validating gridded snow products or assimilating data into models.     

Observations of late winter Canadian tundra snow cover properties

Tundra snow cover is important to monitor as it influences local, regional, and global‐scale surface water balance, energy fluxes, as well as ecosystem and permafrost dynamics. Observations are already showing a decrease in spring snow cover duration at high latitudes, but the impact of changing winter season temperature and precipitation on variables such as snow water equivalent (SWE) is less clear. A multi‐year project was initiated in 2004 with the objective to quantify tundra snow cover properties over multiple years at a scale appropriate for comparison with satellite passive microwave remote sensing data and regional climate and hydrological models. Data collected over seven late winter field campaigns (2004 to 2010) show the patterns of snow depth and SWE are strongly influenced by terrain characteristics. Despite the spatial heterogeneity of snow cover, several inter‐annual consistencies were identified. A regional average density of 0.293 g/cm³ was derived and shown to have little difference with individual site densities when deriving SWE from snow depth measurements. The inter‐annual patterns of SWE show that despite variability in meteorological forcing, there were many consistent ratios between the SWE on flat tundra and the SWE on lakes, plateaus, and slopes. A summary of representative inter‐annual snow stratigraphy from different terrain categories is also presented. © 2013 Her Majesty the Queen in Right of Canada. Hydrological Processes. © 2013 John Wiley & Sons, Ltd.     

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 accumulation and ablation response to changes in forest structure and snow surface albedo after attack by mountain pine beetle

This study quantified changes in snow accumulation and ablation with forest defoliation in a young pine stand attacked by mountain pine beetle, a mature mixed species stand, and a clearcut in south‐central British Columbia. From 2006 to 2012, as trees in the pine stand turned from green to grey, average canopy transmittance increased from 27% to 49%. In the mixed stand, transmittance remained constant at 19%. In 2009, the year of greatest needle loss, average snow surface litter cover in the pine stand was 29% (range 4 – 61%), compared to ≤9% in other years and over double that in the mixed stand. By 2012, litter accumulation in the now‐grey pine stand was only a sixth of that in the mixed stand. Inter‐annual variability in the weather had the greatest effect on snow accumulation and ablation, with the greatest differences between both forested stands and the clearcut occurring in 2010, the year of lowest SWE. Differences in snow accumulation between the pine and mixed stand increased in 2010 as a result of decreased snow interception in the young stand after needlefall. Average ablation rates in the attacked stand were most different from the mixed stand in 2009 and 2012, the years with the largest and smallest over‐winter needle loss, respectively. This study shows that grey, non‐pine, and understory trees moderate snow response to changes in the main canopy. It also highlights the complex interrelationships between ecohydrological processes key to assessing watershed response to forest cover loss in snow dominated hydrologic regimes. Copyright © 2012 John Wiley & Sons, Ltd.     

Improvement of snow depth retrieval for FY3B-MWRI in China

The primary objective of this work is to develop an operational snow depth retrieval algorithm for the FengYun3B Microwave Radiation Imager(FY3B-MWRI)in China.Based on 7-year(2002–2009)observations of brightness temperature by the Advanced Microwave Scanning Radiometer-EOS(AMSR-E)and snow depth from Chinese meteorological stations,we develop a semi-empirical snow depth retrieval algorithm.When its land cover fraction is larger than 85%,we regard a pixel as pure at the satellite passive microwave remote-sensing scale.A 1-km resolution land use/land cover(LULC)map from the Data Center for Resources and Environmental Sciences,Chinese Academy of Sciences,is used to determine fractions of four main land cover types(grass,farmland,bare soil,and forest).Land cover sensitivity snow depth retrieval algorithms are initially developed using AMSR-E brightness temperature data.Each grid-cell snow depth was estimated as the sum of snow depths from each land cover algorithm weighted by percentages of land cover types within each grid cell.Through evaluation of this algorithm using station measurements from 2006,the root mean square error(RMSE)of snow depth retrieval is about 5.6 cm.In forest regions,snow depth is underestimated relative to ground observation,because stem volume and canopy closure are ignored in current algorithms.In addition,comparison between snow cover derived from AMSR-E and FY3B-MWRI with Moderate-resolution Imaging Spectroradiometer(MODIS)snow cover products(MYD10C1)in January 2010 showed that algorithm accuracy in snow cover monitoring can reach 84%.Finally,we compared snow water equivalence(SWE)derived using FY3B-MWRI with AMSR-E SWE products in the Northern Hemisphere.The results show that AMSR-E overestimated SWE in China,which agrees with other validations.     

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.     

Estimating the distribution of snow water equivalent and snow extent beneath cloud cover in the Salt–Verde River basin,Arizona

The temporal and spatial continuity of spatially distributed estimates of snow‐covered area (SCA) are limited by the availability of cloud‐free satellite imagery; this also affects spatial estimates of snow water equivalent (SWE), as SCA can be used to define the extent of snow telemetry (SNOTEL) point SWE interpolation. In order to extend the continuity of these estimates in time and space to areas beneath the cloud cover, gridded temperature data were used to define the spatial domain of SWE interpolation in the Salt–Verde watershed of Arizona. Gridded positive accumulated degree‐days (ADD) and binary SCA (derived from the Advanced Very High Resolution Radiometer (AVHRR)) were used to define a threshold ADD to define the area of snow cover. The optimized threshold ADD increased during snow accumulation periods, reaching a peak at maximum snow extent. The threshold then decreased dramatically during the first time period after peak snow extent owing to the low amount of energy required to melt the thin snow cover at lower elevations. The area having snow cover at this later time was then used to define the area for which SWE interpolation was done. The area simulated to have snow was compared with observed SCA from AVHRR to assess the simulated snow map accuracy. During periods without precipitation, the average commission and omission errors of the optimal technique were 7% and 11% respectively, with a map accuracy of 82%. Average map accuracy decreased to 75% during storm periods, with commission and omission errors equal to 11% and 12% respectively. The analysis shows that temperature data can be used to help estimate the snow extent beneath clouds and therefore improve the spatial and temporal continuity of SCA and SWE products. Copyright © 2004 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.     

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.     

Understanding subgrid variability of snow depth at 1‐km scale using Lidar measurements

It is well known that snow plays an important role in land surface energy balance; however, modelling the subgrid variability of snow is still a challenge in large‐scale hydrological and land surface models. High‐resolution snow depth data and statistical methods can reveal some characteristics of the subgrid variability of snow depth, which can be useful in developing models for representing such subgrid variability. In this study, snow depth was measured by airborne Lidar at 0.5‐m resolution over two mountainous areas in south‐western Wyoming, Snowy Range and Laramie Range. To characterize subgrid snow depth spatial distribution, measured snow depth data of these two areas were meshed into 284 grids of 1‐km × 1‐km. Also, nine representative grids of 1‐km × 1‐km were selected for detailed analyses on the geostatistical structure and probability density function of snow depth. It was verified that land cover is one of the important factors controlling spatial variability of snow depth at the 1‐km scale. Probability density functions of snow depth tend to be Gaussian distributions in the forest areas. However, they are eventually skewed as non‐Gaussian distribution, largely due to the no‐snow areas effect, mainly caused by snow redistribution and snow melt. Our findings show the characteristics of subgrid variability of snow depth and clarify the potential factors that need to be considered in modelling subgrid variability of snow depth.     

Intercomparison of measurements of bulk snow density and water equivalent of snow cover with snow core samplers: Instrumental bias and variability induced by observers

Manually collected snow data are often considered as ground truth for many applications such as climatological or hydrological studies. However, there are many sources of uncertainty that are not quantified in detail. For the determination of water equivalent of snow cover (SWE), different snow core samplers and scales are used, but they are all based on the same measurement principle. We conducted two field campaigns with 9 samplers commonly used in observational measurements and research in Europe and northern America to better quantify uncertainties when measuring depth, density and SWE with core samplers. During the first campaign, as a first approach to distinguish snow variability measured at the plot and at the point scale, repeated measurements were taken along two 20 m long snow pits. The results revealed a much higher variability of SWE at the plot scale (resulting from both natural variability and instrumental bias) compared to repeated measurements at the same spot (resulting mostly from error induced by observers or very small scale variability of snow depth). The exceptionally homogeneous snowpack found in the second campaign permitted to almost neglect the natural variability of the snowpack properties and focus on the separation between instrumental bias and error induced by observers. Reported uncertainties refer to a shallow, homogeneous tundra-taiga snowpack less than 1 m deep (loose, mostly recrystallised snow and no wind impact). Under such measurement conditions, the uncertainty in bulk snow density estimation is about 5% for an individual instrument and is close to 10% among different instruments. Results confirmed that instrumental bias exceeded both the natural variability and the error induced by observers, even in the case when observers were not familiar with a given snow core sampler.     

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.     

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.     

Snow density variations: consequences for ground‐penetrating radar

Reliable hydrological forecasts of snowmelt runoff are of major importance for many areas. Ground‐penetrating radar (GPR) measurements are used to assess snowpack water equivalent for planning of hydropower production in northern Sweden. The travel time of the radar pulse through the snow cover is recorded and converted to snow water equivalent (SWE) using a constant snowpack mean density from the drainage basin studied. In this paper we improve the method to estimate SWE by introducing a depth‐dependent snowpack density. We used 6 years measurements of peak snow depth and snowpack mean density at 11 locations in the Swedish mountains. The original method systematically overestimates the SWE at shallow depths (+25% for 0·5 m) and underestimates the SWE at large depths (?35% for 2·0 m). A large improvement was obtained by introducing a depth–density relation based on average conditions for several years, whereas refining this by using separate relations for individual years yielded a smaller improvement. The SWE estimates were substantially improved for thick snow covers, reducing the average error from 162 ± 23 mm to 53 ± 10 mm for depth range 1·2–2·0 m. Consequently, the introduction of a depth‐dependent snow density yields substantial improvements of the accuracy in SWE values calculated from GPR data. Copyright © 2005 John Wiley & Sons, Ltd.     

Snow-induced thermal variations around a single conifer tree

The influence of trees on the ground thermal regime is important to the overall winter energy exchange in a snow-covered, forested watershed. In this work, spatial zones around a single conifer tree were defined and examined for their controls on the snow cover, snow-ground interface temperatures and frozen ground extent. A large white spruce (Picea glauca), approximately 18 m tall with a crown diameter of 7.5 m and located in northern Vermont, was the subject of this study. The tree was instrumented with thermistors to measure the snow-ground interface temperature between the tree trunk and 6 m from the tree into undisturbed snow. Four distinct zones around the conifer are defined that affect the snow distribution characteristics: adjacent to the trunk; the tree well; the tree crown perimeter; and the unaffected area away from the tree. At the time of peak snow accumulation and during the ablation season, snow depth and density profiles were measured. The area beneath the canopy accumulated 34% of the snow accumulated in the undisturbed zone. By the end of the ablation season, the depth of snow under the canopy had decreased to 18% of the undisturbed snow depth. The tree and branch characteristics of spruce in this temperate climate resulted in a different snow depth profile compared with previous empirical relationships around a single conifer. A new relationship is presented for snow distribution around conifer trees that has the ability to better fit data from a variety of conifer types than previously published relationships. Less snow beneath the canopy led to colder snow-ground interface temperatures than measured in undisturbed snow. The depth of frozen ground in the different zones was modelled using a simple analytical solution that showed deeper frost penetration in the tree well than beneath the undisturbed snow.     

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.     

Evaluating the potential of a blended passive microwave‐interactive multi‐sensor product for improved mapping of snow cover and estimations of snow water equivalent

The retrieval of Snow Water Equivalent (SWE) from remote sensing satellites continues to be a very challenging problem. In this paper, we evaluate the accuracy of a new SWE product derived from the blending of a passive microwave SWE product based on the Advanced Microwave Sounding Unit (AMSU) with a multi‐sensor snow cover extent product based on the Interactive Multi‐sensor Snow and Ice Mapping System (IMS). The microwave measurements have the ability to penetrate the snow pack, and thus, the retrieval of SWE is best accomplished using the AMSU. On the other hand, the IMS maps snow cover more reliably due to the use of multiple satellite and ground observations. The evolution of global snow cover from the blended, the AMSU and the IMS products was examined during the 2006 snow season. Despite the overall good inter‐product agreement, it was shown that the retrievals of snow cover extent in the blended product are improved when using IMS, with implications for improved microwave retrievals of SWE. In a separate investigation, the skill of the microwave SWE product was also examined for its ability to correctly estimate SWE globally and regionally. Qualitative evaluation of global SWE retrievals suggested dependence on land surface temperature: the lower the temperature, the higher the SWE retrieved. This temperature bias was attributed in part to temperature effects on those snow properties that impact microwave response. Therefore, algorithm modifications are needed with more dynamical adjustments to account for changing snow cover. Quantitative evaluation over Slovakia in central Europe, for a limited period in 2006, showed reasonably good performance for SWE less than 100 mm. Sensitivity to deeper snow decreased significantly. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluation of spatial variability in snow water equivalent for a high mountain catchment

Multivariate statistical analysis was used to explore relationships between catchment topography and spatial variability in snow accumulation and melt processes in a small headwater catchment in the Spanish Pyrenees. Manual surveys of snow depth and density provided information on the spatial distribution of snow water equivalent (SWE) and its depletion over the course of the 1997 and 1998 melt seasons. A number of indices expressing the topographic control on snow processes were extracted from a detailed digital elevation model of the catchment. Bivariate screening was used to assess the relative importance of these topographic indices in controlling snow accumulation at the start of the melt season, average melt rates and the timing of snow disappearance. This suggested that topographic controls on the redistribution of snow by wind are the most important influence on snow distribution at the start of the melt season. Furthermore, it appeared that spatial patterns of snow disappearance were largely determined by the distribution of snow water equivalent (SWE) at the start of the melt season, rather than by spatial variability in melt rates during the melt season. Binary regression tree models relating snow depth and disappearance date to terrain indices were then constructed. These explained 70–80% of the variance in the observed data. As well as providing insights into the influence of topography on snow processes, it is suggested that the techniques presented herein could be used in the parameterization of distributed snowmelt models, or in the design of efficient stratified snow surveys. Copyright © 2003 John Wiley & Sons, Ltd.     

Deriving snow‐cover depletion curves for different spatial scales from remote sensing and snow telemetry data

During the melting of a snowpack, snow water equivalent (SWE) can be correlated to snow‐covered area (SCA) once snow‐free areas appear, which is when SCA begins to decrease below 100%. This amount of SWE is called the threshold SWE. Daily SWE data from snow telemetry stations were related to SCA derived from moderate‐resolution imaging spectroradiometer images to produce snow‐cover depletion curves. The snow depletion curves were created for an 80 000 km² domain across southern Wyoming and northern Colorado encompassing 54 snow telemetry stations. Eight yearly snow depletion curves were compared, and it is shown that the slope of each is a function of the amount of snow received. Snow‐cover depletion curves were also derived for all the individual stations, for which the threshold SWE could be estimated from peak SWE and the topography around each station. A station's peak SWE was much more important than the main topographic variables that included location, elevation, slope, and modelled clear sky solar radiation. The threshold SWE mostly illustrated inter‐annual consistency. Copyright © 2015 John Wiley & Sons, Ltd.