Optimal design of snow stake networks to estimate snow depth in an alpine mountain range

Monitoring and estimation of snow depth in alpine catchments is needed for a proper assessment of management alternatives for water supply in these water resources systems. The distribution of snowpack thickness is usually approached by using field data that come from snow samples collected at a given number of locations that constitute the monitoring network. Optimal design of this network is required to obtain the best possible estimates. Assuming that there is an existing monitoring network, its optimization may imply the selection of an optimal network as a subset of the existing one (if there are no funds to maintain them) or enlarging the existing network by one or more stations (optimal augmentation problem). We propose an optimization procedure that minimizes the total variance in the estimate of snowpack thickness. The novelty of this work is to treat, for the first time, the problem of snow observation network optimization for an entire mountain range rather than for small catchments as done in the previous studies. Taking into account the reduced data available, which is a common problem in many mountain ranges, the importance of a proper design of these observation networks is even larger. Snowpack thickness is estimated by combining regression models to approach the effect of the explanatory variables and kriging techniques to consider the influence of the stakes location. We solve the optimization problems under different hypotheses, studying the impacts of augmentation and reduction, both, one by one and in pairs. We also analyse the sensitivity of results to nonsnow measurements deduced from satellite information. Finally, we design a new optimal network by combining the reduction and augmentation methods. The methodology has been applied to the Sierra Nevada mountain range (southern Spain), where very limited resources are employed to monitor snowfall and where an optimal snow network design could prove critical. An optimal snow observation network is defined by relocating some observation points. It would reduce the estimation variance by around 600 cm² (15%).     

Inference and uncertainty of snow depth spatial distribution at the kilometre scale in the Colorado Rocky Mountains: the effects of sample size,random sampling,predictor quality,and validation procedures

Historically, observing snow depth over large areas has been difficult. When snow depth observations are sparse, regression models can be used to infer the snow depth over a given area. Data sparsity has also left many important questions about such inference unexamined. Improved inference, or estimation, of snow depth and its spatial distribution from a given set of observations can benefit a wide range of applications from water resource management, to ecological studies, to validation of satellite estimates of snow pack. The development of Light Detection and Ranging (LiDAR) technology has provided non‐sparse snow depth measurements, which we use in this study, to address fundamental questions about snow depth inference using both sparse and non‐sparse observations. For example, when are more data needed and when are data redundant? Results apply to both traditional and manual snow depth measurements and to LiDAR observations. Through sampling experiments on high‐resolution LiDAR snow depth observations at six separate 1.17‐km² sites in the Colorado Rocky Mountains, we provide novel perspectives on a variety of issues affecting the regression estimation of snow depth from sparse observations. We measure the effects of observation count, random selection of observations, quality of predictor variables, and cross‐validation procedures using three skill metrics: percent error in total snow volume, root mean squared error (RMSE), and R². Extremes of predictor quality are used to understand the range of its effect; how do predictors downloaded from internet perform against more accurate predictors measured by LiDAR? Whereas cross validation remains the only option for validating inference from sparse observations, in our experiments, the full set of LiDAR‐measured snow depths can be considered the ‘true’ spatial distribution and used to understand cross‐validation bias at the spatial scale of inference. We model at the 30‐m resolution of readily available predictors, which is a popular spatial resolution in the literature. Three regression models are also compared, and we briefly examine how sampling design affects model skill. Results quantify the primary dependence of each skill metric on observation count that ranges over three orders of magnitude, doubling at each step from 25 up to 3200. Whereas uncertainty (resulting from random selection of observations) in percent error of true total snow volume is typically well constrained by 100–200 observations, there is considerable uncertainty in the inferred spatial distribution (R²) even at medium observation counts (200–800). We show that percent error in total snow volume is not sensitive to predictor quality, although RMSE and R² (measures of spatial distribution) often depend critically on it. Inaccuracies of downloaded predictors (most often the vegetation predictors) can easily require a quadrupling of observation count to match RMSE and R² scores obtained by LiDAR‐measured predictors. Under cross validation, the RMSE and R² skill measures are consistently biased towards poorer results than their true validations. This is primarily a result of greater variance at the spatial scales of point observations used for cross validation than at the 30‐m resolution of the model. The magnitude of this bias depends on individual site characteristics, observation count (for our experimental design), and sampling design. Sampling designs that maximize independent information maximize cross‐validation bias but also maximize true R². The bagging tree model is found to generally outperform the other regression models in the study on several criteria. Finally, we discuss and recommend use of LiDAR in conjunction with regression modelling to advance understanding of snow depth spatial distribution at spatial scales of thousands of square kilometres. Copyright © 2012 John Wiley & Sons, Ltd.     

The role of meteorological forcing and snow model complexity in winter glacier mass balance estimation,Columbia River basin,Canada

Glaciers are commonly located in mountainous terrain subject to highly variable meteorological conditions. High resolution meteorological (HRM) data simulated by atmospheric models can complement meteorological station observations in order to assess changes in glacier energy fluxes and mass balance. We examine the performance of two snow models, SnowModel and Alpine3D, forced by different meteorological data for winter mass balance simulations at four glaciers in the Canadian portion of the Columbia Basin. The Weather Research and Forecasting model (WRF) with resolution of 1 km and the North American Land Data Assimilation System with ~12 km resolution, provide HRM data for the two snow models. Evaluation is based on the ability of the snow models to simulate snow depth at both point locations (automated snow weather stations) and over the entire glacier surface (airborne LiDAR [Light Detection and Ranging] surveys) during the 2015/2016 winter accumulation. When forced with HRM data, both models can reproduce snow depth to within ±15% of observed values. Both models underestimate winter mass balance when forced by HRM data. When driven with WRF data, SnowModel underestimates winter mass balance integrated over the glacier area by 1 and 10%, whilst Alpine3D underestimates winter mass balance by 12 and 22% compared with LiDAR and stake measurements, respectively. The overall results show that SnowModel forced by WRF simulated winter mass balance the best.     

Updating of snow depletion curve with remote sensing data

A method for using remotely sensed snow cover information in updating a hydrological model is developed, based on Bayes' theorem. A snow cover mass balance model structure adapted to such use of satellite data is specified, using a parametric snow depletion curve in each spatial unit to describe the subunit variability in snow storage. The snow depletion curve relates the accumulated melt depth to snow‐covered area, accumulated snowmelt runoff volume, and remaining snow water equivalent. The parametric formulation enables updating of the complete snow depletion curve, including mass balance, by satellite data on snow coverage. Each spatial unit (i.e. grid cell) in the model maintains a specific depletion curve state that is updated independently. The uncertainty associated with the variables involved is formulated in terms of a joint distribution, from which the joint expectancy (mean value) represents the model state. The Bayesian updating modifies the prior (pre‐update) joint distribution into a posterior, and the posterior joint expectancy replaces the prior as the current model state. Three updating experiments are run in a 2400 km² mountainous region in Jotunheimen, central Norway (61°N, 9°E) using two Landsat 7 ETM+ images separately and together. At 1 km grid scale in this alpine terrain, three parameters are needed in the snow depletion curve. Despite the small amount of measured information compared with the dimensionality of the updated parameter vector, updating reduces uncertainty substantially for some state variables and parameters. Parameter adjustments resulting from using each image separately differ, but are positively correlated. For all variables, uncertainty reduction is larger with two images used in conjunction than with any single image. Where the observation is in strong conflict with the prior estimate, increased uncertainty may occur, indicating that prior uncertainty may have been underestimated. Copyright © 2006 John Wiley & Sons, Ltd.     

Evaluation of snow models in terrestrial biosphere models using ground observation and satellite data: impact on terrestrial ecosystem processes

Snow is important for water management, and an important component of the terrestrial biosphere and climate system. In this study, the snow models included in the Biome‐BGC and Terrestrial Observation and Prediction System (TOPS) terrestrial biosphere models are compared against ground and satellite observations over the Columbia River Basin in the US and Canada and the impacts of differences in snow models on simulated terrestrial ecosystem processes are analysed. First, a point‐based comparison of ground observations against model and satellite estimates of snow dynamics are conducted. Next, model and satellite snow estimates for the entire Columbia River Basin are compared. Then, using two different TOPS simulations, the default TOPS model (TOPS with TOPS snow model) and the TOPS model with the Biome‐BGC snow model, the impacts of snow model selection on runoff and gross primary production (GPP) are investigated. TOPS snow model predictions were consistent with ground and satellite estimates of seasonal and interannual variations in snow cover, snow water equivalent, and snow season length; however, in the Biome‐BGC snow model, the snow pack melted too early, leading to extensive underpredictions of snow season length and snow covered area. These biases led to earlier simulated peak runoff and reductions in summer GPP, underscoring the need for accurate snow models within terrestrial ecosystem models. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluation of snow water equivalent datasets over the Saint‐Maurice river basin region of southern Québec

A 10‐km gridded snow water equivalent (SWE) dataset is developed over the Saint‐Maurice River basin region in southern Québec from kriging of observed snow survey data for evaluation of SWE products. The gridded SWE dataset covers 1980–2014 and is based on manual gravimetric snow surveys carried out on February 1, March 1, March 15, April 1, and April 15 of each snow season, which captures the annual maximum SWE (SWEM) with a mean interpolation error of ±19%. The dataset is used to evaluate SWEM from a range of sources including satellite retrievals, reanalyses, Canadian regional climate models, and the Canadian Meteorological Centre operational snow depth analysis. We also evaluate a number of solid precipitation datasets to determine their contribution to systematic errors in estimated SWEM. None of the evaluated datasets is able to provide estimates of SWEM that are within operational requirements of ±15% error, and insufficient solid precipitation is determined to be one of the main reasons. The Climate System Forecast Reanalysis is the only dataset where snowfall is sufficiently large to generate SWEM values comparable to observations. Inconsistencies in precipitation are also found to have a strong impact on year‐to‐year variability in SWEM dataset performance and spread. Version 3.6.1 of the Canadian Land Surface Scheme land surface scheme driven with ERA‐Interim output downscaled by Version 5.0.1 of the Canadian Regional Climate Model was the best physically based model at explaining the observed spatial and temporal variability in SWEM (root‐mean‐square error [RMSE] = 33%) and has potential for lower error with adjusted precipitation. Operational snow products relying on the real‐time snow depth observing network performed poorly due to a lack of real‐time data and the strong local scale variability of point snow depth observations. The results underscore the need for more effort to be invested in improving solid precipitation estimates for use in snow hydrology applications.     

Small scale spatial variability of snow density and depth over complex alpine terrain: Implications for estimating snow water equivalent

This study analyzes spatial variability of snow depth and density from measurements made in February and April of 2010 and 2011 in three 1–2 km² areas within a valley of the central Spanish Pyrenees. Snow density was correlated with snow depth and different terrain characteristics. Regression models were used to predict the spatial variability of snow density, and to assess how the error in computed densities might influence estimates of snow water equivalent (SWE).The variability in snow depth was much greater than that of snow density. The average snow density was much greater in April than in February. The correlations between snow depth and density were generally statistically significant but typically not very high, and their magnitudes and signs were highly variable among sites and surveys. The correlation with other topographic variables showed the same variability in magnitude and sign, and consequently the resulting regression models were very inconsistent, and in general explained little of the variance. Antecedent climatic and snow conditions prior to each survey help highlight the main causes of the contrasting relation shown between snow depth, density and terrain. As a consequence of the moderate spatial variability of snow density relative to snow depth, the absolute error in the SWE estimated from computed densities using the regression models was generally less than 15%. The error was similar to that obtained by relating snow density measurements directly to adjacent snow depths.     

Dynamics of snow ablation in a small Alpine catchment observed by repeated terrestrial laser scans

The spatio‐temporal distribution of snow in a catchment during ablation reflects changes in the total amount of snow water equivalent and is thus a key parameter for the estimation of melt water run‐off. This study explores possible rules behind the spatial variability of snow depth during the ablation season in a small Alpine catchment with complex topography. The snow depth observations are based on more than 160 000 terrestrial laser scanner data points with a spatial resolution of 1 m, which were obtained from 11 scanning campaigns of two consecutive ablation seasons. The analysis suggests that for estimating cumulative snow melt dynamics from the catchment investigated, assessing the initial snow distribution prior to the melt season is more important than addressing spatial differences in the melt behaviour. Snow volume and snow‐covered area could be predicted well using a conceptual melt model assuming spatially uniform melt rates. However, accurate results were only obtained if the model was initialized with a pre‐melt snow distribution that reflected measured mean and standard deviation. Using stratified melt rates on the other hand did not improve the model results. At least for sites with similar meteorological and topographical conditions, the model approach presented here comprises an efficient way to estimate snow depletion dynamics, especially if persistent snow accumulation pattern between years facilitate the characterization of the initial snow distribution prior to the melt. Copyright © 2011 John Wiley & Sons, Ltd.     

Predicting new snow density in the Italian Alps: A variability analysis based on 10 years of measurements

Despite its strong impact on the time evolution of the snowpack, current estimation of new snow density (ρ_hn) is usually accomplished either by using local empirical techniques or by assuming a constant snow density. Faced with the lack of an estimation model of ρ_hn valid for a wide spatial scale and supported by a suitable number of observations, this study aims to develop simple monthly linear regression models at scale of the entire Italian Alpine chain based on 12,112 snowfall observations at 122 stations, using only air temperature as predictor. Moreover, the remaining variance is investigated in both time and space, also considering some qualitative features of the snowfall events. The daily ρ_hn measurements present a mean value of 115 kg m^?3 (105 and 159 kg m^?3 for dry and wet conditions, respectively). The mean air temperature of the 24 hr preceding the snowfall event has been found to be the best predictor of the ρ_hn, within 31% of uncertainty. The analysis of associated residues allows supporting the idea that the adoption of a more local approach than the one analysed here is not able to substantially increase the predictive capabilities of the model. In fact, the main factor explaining the remaining variance over the air temperature is the wind, but in a complex orography, as mountain regions are, supplying realistic local wind fields is particularly challenging. Therefore, we conclude that using only the daily mean temperature as predictor is a good choice for estimating daily new snow density at scale of Italian Alpine chain, as well as at more regional scale.     

Subgrid variability of snow water equivalent at operational snow stations in the western USA

The spatial distribution of snow water equivalent (SWE) is a key variable in many regional‐scale land surface models. Currently, the assimilation of point‐scale snow sensor data into these models is commonly performed without consideration of the spatial representativeness of the point data with respect to the model grid‐scale SWE. To improve the understanding of the relationship between point‐scale snow measurements and surrounding areas, we characterized the spatial distribution of snow depth and SWE within 1‐, 4‐ and 16‐km² grids surrounding 15 snow stations (snowpack telemetry and California snow sensors) in California, Colorado, Wyoming, Idaho and Oregon during the 2008 and 2009 snow seasons. More than 30 000 field observations of snowpack properties were used with binary regression tree models to relate SWE at the sensor site to the surrounding area SWE to evaluate the sensor representativeness of larger‐scale conditions. Unlike previous research, we did not find consistent high biases in snow sensor depth values as biases over all sites ranged from 74% overestimates to 77% underestimates. Of the 53 assessments, 27 surveys indicated snow station biases of less than 10% of the surrounding mean observed snow depth. Depth biases were largely dictated by the physiographic relationship between the snow sensor locations and the mean characteristics of the surrounding grid, in particular, elevation, solar radiation index and vegetation density. These scaling relationships may improve snow sensor data assimilation; an example application is illustrated for the National Operational Hydrologic Remote Sensing Center National Snow Analysis SWE product. The snow sensor bias information indicated that the assimilation of point data into the National Operational Hydrologic Remote Sensing Center model was often unnecessary and reduced model accuracy. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating the snow water equivalent from snow depth measurements in the Swiss Alps

The snow water equivalent (SWE) characterizes the hydrological significance of snow cover. However, measuring SWE is time-consuming, thus alternative methods of determining SWE may be useful. SWE can be calculated from snow depth if the bulk snow density is known. Thus, a reliable estimation method of snow densities could (a) potentially save a lot of effort by, at least partly, sampling snow depth instead of SWE, and would (b) allow snow hydrological evaluations, when only snow depth data are available. To generate a useful parameterization of the bulk density a large dataset was analyzed covering snow densities and depths measured biweekly over five decades at 37 sites throughout the Swiss Alps. Four factors were identified to affect the bulk snow density: season, snow depth, site altitude, and site location. These factors constitute a convenient set of input variables for a snow density model developed in this study. The accuracy of estimating SWE using our model is shown to be equivalent to the variability of repeated SWE measurements at one site. The technique may therefore allow a more efficient but indirect sampling of the SWE without necessarily affecting the data quality.     

Model complexity and data requirements in snow hydrology: seeking a balance in practical applications

We investigate the problem of balancing model complexity and input data requirements in snow hydrology. For this purpose, we analyze the performance of two models of different complexity in estimating variables of interest in snow hydrology applications. These are snow depth, bulk snow density, snow water equivalent and snowmelt run‐off. We quantify the differences between data and model prediction using 18 years of measurements from an experimental site in the French Alps (Col de Porte, 1325 m AMSL). The models involved in this comparison are a one‐layer temperature‐index model (HyS) and a multilayer model (Crocus). Results show that the expected loss in performance in the one‐layer temperature‐index model with respect to the multilayer model is low when considering snow depth, snow water equivalent and bulk snow density. As for run‐off, the comparison returns less clear indications for identification of a balance. In particular, differences between the models' prediction and data with an hourly resolution are higher when considering the Crocus model than the HyS model. However, Crocus is better at reproducing sub‐daily cycles in this variable. In terms of daily run‐off, the multilayer physically based model seems to be a better choice, while results in terms of cumulative run‐off are comparable. The better reproduction of daily and sub‐daily variability of run‐off suggests that use of the multilayer model may be preferable for this purpose. Variation in performance is discussed as a function of both the calibration solution chosen and the time of year. Copyright © 2016 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.     

Characteristics of snow cover in the Hindukush,Karakoram and Himalaya region using Landsat satellite data

Estimation of snow cover characteristics (snow grain size, snow contamination, snow depth and liquid water content) from satellite data are important components for many hydrological models used for the water resource management. This research aimed to use Landsat ETM+ (Enhanced Thematic Mapper Plus) data to estimate the snow pack characteristics in northern Pakistan. The Normalized Difference Snow Index (NDSI), Snow Contamination Index (SCI) and Snow Grain Size Index (SGI) are applied to estimate the snow cover characteristics in northern Pakistan for the first time. Qualitative maps are generated to show the snow cover distribution, snow contamination concentration and snow grain size distribution over snow cover area. Our results show that NDSI, SCI and SGI can be effectively used to identify snow area, contaminated snow and ageing snow. Furthermore, the results of the current study indicate that in the HKH region 99.8% of the snow is least contaminated whereas 94.50% of the area has fine snow grain size. As no such attempt in the past has been made in northern Pakistan, it is envisaged that the results of this study will be helpful in planning remote sensing data for the water resource management and in characterizing the snow cover for the climate change applications in the region. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling snowpack bulk density using snow depth,cumulative degree-days,and climatological predictor variables

Snowpack water equivalent (SWE) is a key variable for water resource management in snow-dominated catchments. While it is not feasible to quantify SWE at the catchment scale using either field surveys or remotely sensed data, technologies such as airborne LiDAR (light detection and ranging) support the mapping of snow depth at scales relevant to operational water management. To convert snow depth to water equivalent, models have been developed to predict SWE or snowpack density based on snow depth and additional predictor variables. This study builds upon previous models that relate snowpack density to snow depth by including additional predictor variables to account for (1) long-term climatologies that describe the prevailing conditions influencing regional snowpack properties, and (2) the effect of intra- and inter-year variability in meteorological conditions on densification through a cumulative degree-day index derived from North American Regional Reanalysis products. A non-linear model was fit to 114 506 snow survey measurements spanning 41 years from 1166 snow courses across western North America. Under spatial cross-validation, the predicted densities had a root-mean-square error of 47.1 kg m⁻³, a mean bias of −0.039 kg m⁻³, and a Nash-Sutcliffe Efficiency of 0.70. The model developed in this study had similar overall performance compared to a similar regression-based model reported in the literature, but had reduced seasonal biases. When applied to predict SWE from simulated depths with random errors consistent with those obtained from LiDAR or Structure-from-Motion, 50% of the SWE estimates for April and May fell within −45 to 49 mm of the observed SWE, representing prediction errors of −15% to 20%.     

The distribution of daily snow water equivalent in the central Italian Alps

The statistical distribution of the daily Snow Water Equivalent (SWE) is investigated for a network of gauging stations in the Alpine part of Lombardia region, in the central Italian Alps. An event based data analysis is carried out using a 14 year long data set dating back to 1989. SWE is estimated when the new snow depth is greater than 6 cm. The SWE sample average in time is shown to be related to physiographic attributes of the gauging area, thus not being homogeneous in space. The values of SWE scaled by their average, or index value, instead show well approximated homogeneity of the second order moment, or coefficient of variation, in space. This suggests the use of a regional approach for frequency estimation of SWE. The frequency of occurrence of the normalized values of SWE is evaluated and tentatively accommodated by four probability distributions, often adopted in statistical modeling of hydrological variables. The Lognormal distribution shows the best performance. Single site distribution fitting is then carried out using the regional distribution, providing satisfactory results.     

Watershed‐scale controls on snow accumulation in a small montane watershed,southwestern Alberta,Canada

Snow accumulation in mountain headwater basins is a major water source, particularly in semi‐arid environments such as southern Alberta where water resources are stressed and snowmelt supplies more than 80% of downstream runoff. Relationships between landscape predictor variables and snow water equivalent (SWE) were quantified by combining field and LiDar measurements with classification and regression tree analysis over two winter seasons (2010 and 2011) in a small, montane watershed. 2010 was a below average snow accumulation year, while 2011 was well above normal. In both the field and regression tree data, elevation was the dominant control on snow distribution in both years, although snow distribution was driven by melt processes in 2010 and accumulation processes in 2011. The importance of solar radiation and wind exposure was represented in the regression trees in both years. The regression trees also noted the lower importance of canopy closure, slope, and aspect, which was not observed in the field data. This technique could provide an additional method of forecasting annual water supply from melting snow. However, further research is required to address the lack of data collected above treeline, to provide a full‐basin estimate of SWE. Copyright © 2012 John Wiley & Sons, Ltd.     

Inter‐annual variation in the topographic controls on catchment‐scale snow distribution in a maritime alpine catchment,New Zealand

Seasonal snow is a globally important water resource that contributes substantially to upland and lowland water resources. As such, there is a need to understand the controls on the spatial and temporal variation in snow distribution. This study meets this research need by investigating the topographic controls on snow depth distribution in the upper Jollie catchment in the Southern Alps of New Zealand. Furthermore, inter‐annual variation in the importance of the topographic controls is examined and linked to variation in the dominant synoptic‐scale weather patterns over a 4‐year period (2007–2010). Through the use of regression trees, the relative importance of the topographic controls on snow depth was shown to vary between the four study years. In particular, elevation explained the greatest amount of variance in 2007 and 2008 and east‐exposure explained the greatest variance in 2009 and 2010. The other wind exposure variables also had a large effect on the snow depth distribution in 2009 and 2010. Differences in the frequency and duration of synoptic weather patterns were physically consistent with the changing importance of these variables. In particular, a higher frequency of troughing events in 2009 and 2010 is thought to be associated with a reduced importance of elevation and greater influence of wind exposure on snow depth in these years. These findings demonstrate the importance of using multi‐year data sets, and of considering topographic and climatic influences, when attempting to model alpine snow distribution. Copyright © 2014 John Wiley & Sons, Ltd.     

Assimilating snow observations to snow interception process simulations

Snow interception is a crucial hydrological process in cold regions needleleaf forests, but is rarely measured directly. Indirect estimates of snow interception can be made by measuring the difference in the increase in snow accumulation between the forest floor and a nearby clearing over the course of a storm. Pairs of automatic weather stations with acoustic snow depth sensors provide an opportunity to estimate this, if snow density can be estimated reliably. Three approaches for estimating fresh snow density were investigated: weighted post-storm density increments from the physically based Snobal model, fresh snow density estimated empirically from air temperature (Hedstrom, N. R., et al. [1998]. Hydrological Processes, 12, 1611–1625), and fresh snow density estimated empirically from air temperature and wind speed (Jordan, R. E., et al. [1999]. Journal of Geophysical Research, 104, 7785–7806). Automated snow depth observations from adjacent forest and clearing sites and estimated snow densities were used to determine snowstorm snow interception in a subalpine forest in the Canadian Rockies, Alberta, Canada. Then the estimated snow interception and measured interception information from a weighed, suspended tree and a time-lapse camera were assimilated into a model, which was created using the Cold Regions Hydrological Modelling platform (CRHM), using Ensemble Kalman Filter or a simple rule-based direct insertion method. Interception determined using density estimates from the Hedstrom-Pomeroy fresh snow density equation agreed best with observations. Assimilating snow interception information from automatic snow depth measurements improved modelled snow interception timing by 7% and magnitude by 13%, compared to an open loop simulation driven by a numerical weather model; its accuracy was close to that simulated using locally observed meteorological data. Assimilation of tree-measured snow interception improved the snow interception simulation timing and magnitude by 18 and 19%, respectively. Time-lapse camera snow interception information assimilation improved the snow interception simulation timing by 32% and magnitude by 7%. The benefits of assimilation were greatly influenced by assimilation frequency and quality of the forcing data.     

Evaluation of snow products over the Tibetan Plateau

Four satellite‐based snow products are evaluated over the Tibetan Plateau for the 2007–2010 snow seasons. The Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua snow cover daily L3 Global 500‐m grid products (MOD10A1 and MYD10A1), the National Oceanic and Atmospheric Administration Interactive Multisensor Snow and Ice Mapping System (IMS) daily Northern Hemisphere snow cover product and the Advanced Microwave Scanning Radiometer – Earth Observing System Daily Snow Water Equivalent were validated against Thematic Mapper (TM) snow cover maps of Landsat‐5 and meteorological station snow depth observations. The overall accuracy of MOD10A1, MYD10A1 and IMS is higher than 91% against stations observations and than 79% against Landsat TM images. In general, the daily MODIS snow cover products show better performance than the multisensor IMS product. However, the IMS snow cover product is suitable for larger scale (~4km) analysis and applications, with the advantage over MODIS to allow for mitigation for cloud cover. The accuracy of the three products decreases with decreasing snow depth. Overestimation errors are most common over forested regions; the IMS and Advanced Microwave Scanning Radiometer – Earth Observing System Snow Water Equivalent products also show poorer performance that the MODIS products over grassland. By identifying weaknesses in the satellite products, this study provides a focus for the improvement of snow products over the Tibetan plateau. The quantitative evaluation of the products proposed here can also be used to assess their relative weight in data assimilation, against other data sources, such as modelling and in situ measurement networks. Copyright © 2015 John Wiley & Sons, Ltd.