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.     

Spatial interpolation of snow water equivalency using surface observations and remotely sensed images of snow‐covered area

As demand for water continues to escalate in the western Unites States, so does the need for accurate monitoring of the snowpack in mountainous areas. In this study, we describe a simple methodology for generating gridded‐estimates of snow water equivalency (SWE) using both surface observations of SWE and remotely sensed estimates of snow‐covered area (SCA). Multiple regression was used to quantify the relationship between physiographic variables (elevation, slope, aspect, clear‐sky solar radiation, etc.) and SWE as measured at a number of sites in a mountainous basin in south‐central Idaho (Big Wood River Basin). The elevation of the snowline, obtained from the SCA estimates, was used to constrain the predicted SWE values. The results from the analysis are encouraging and compare well to those found in previous studies, which often utilized more sophisticated spatial interpolation techniques. Cross‐validation results indicate that the spatial interpolation method produces accurate SWE estimates [mean R² = 0·82, mean mean absolute error (MAE) = 4·34 cm, mean root mean squared error (RMSE) = 5·29 cm]. The basin examined in this study is typical of many mid‐elevation mountainous basins throughout the western United States, in terms of the distribution of topographic variables, as well as the number and characteristics of sites at which the necessary ground data are available. Thus, there is high potential for this methodology to be successfully applied to other mountainous basins. Copyright © 2010 John Wiley & Sons, Ltd.     

Improvement of distributed snowmelt energy balance modeling with MODIS‐based NDSI‐derived fractional snow‐covered area data

Describing the spatial variability of heterogeneous snowpacks at a watershed or mountain‐front scale is important for improvements in large‐scale snowmelt modelling. Snowmelt depletion curves, which relate fractional decreases in snow‐covered area (SCA) against normalized decreases in snow water equivalent (SWE), are a common approach to scale‐up snowmelt models. Unfortunately, the kinds of ground‐based observations that are used to develop depletion curves are expensive to gather and impractical for large areas. We describe an approach incorporating remotely sensed fractional SCA (FSCA) data with coinciding daily snowmelt SWE outputs during ablation to quantify the shape of a depletion curve. We joined melt estimates from the Utah Energy Balance Snow Accumulation and Melt Model (UEB) with FSCA data calculated from a normalized difference snow index snow algorithm using NASA's moderate resolution imaging spectroradiometer (MODIS) visible (0·545–0·565 µm) and shortwave infrared (1·628–1·652 µm) reflectance data. We tested the approach at three 500 m² study sites, one in central Idaho and the other two on the North Slope in the Alaskan arctic. The UEB‐MODIS‐derived depletion curves were evaluated against depletion curves derived from ground‐based snow surveys. Comparisons showed strong agreement between the independent estimates. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

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.     

Snow accumulation distribution inferred from time‐lapse photography and simple modelling

The spatial and temporal distribution of snow accumulation is complex and significantly influences the hydrological characteristics of mountain catchments. Many snow redistribution processes, such as avalanching, slushflow or wind drift, are controlled by topography, but their modelling remains challenging. In situ measurements of snow accumulation are laborious and generally have a coarse spatial or temporal resolution. In this respect, time‐lapse photography shows itself as a powerful tool for collecting information at relatively low cost and without the need for direct field access. In this paper, the snow accumulation distribution of an Alpine catchment is inferred by adjusting a simple snow accumulation model combined with a temperature index melt model to match the modelled melt‐out pattern evolution to the pattern monitored during an ablation season through terrestrial oblique photography. The comparison of the resulting end‐of‐winter snow water equivalent distribution with direct measurements shows that the achieved accuracy is comparable with that obtained with an inverse distance interpolation of the point measurements. On average over the ablation season, the observed melt‐out pattern can be reproduced correctly in 93% of the area visible from the fixed camera. The relations between inferred snow accumulation distribution and topographic variables indicate large scatter. However, a significant correlation with local slope is found and terrain curvature is detected as a factor limiting the maximal snow accumulation. Copyright © 2010 John Wiley & Sons, Ltd.     

A comparison of snowmelt‐derived streamflow from temperature‐index and modified‐temperature‐index snow models

Reliable estimation of the volume and timing of snowmelt runoff is vital for water supply and flood forecasting in snow‐dominated regions. Snowmelt is often simulated using temperature‐index (TI) models due to their applicability in data‐sparse environments. Previous research has shown that a modified‐TI model, which uses a radiation‐derived proxy temperature instead of air temperature as its surrogate for available energy, can produce more accurate snow‐covered area (SCA) maps than a traditional TI model. However, it is unclear whether the improved SCA maps are associated with improved snow water equivalent (SWE) estimation across the watershed or improved snowmelt‐derived streamflow simulation. This paper evaluates whether a modified‐TI model produces better streamflow estimates than a TI model when they are used within a fully distributed hydrologic model. It further evaluates the performance of the two models when they are calibrated using either point SWE measurements or SCA maps. The Senator Beck Basin in Colorado is used as the study site because its surface is largely bedrock, which reduces the role of infiltration and emphasizes the role of the SWE pattern on streamflow generation. Streamflow is simulated using both models for 6 years. The modified‐TI model produces more accurate streamflow estimates (including flow volume and peak flow rate) than the TI model, likely because the modified‐TI model better reproduces the SWE pattern across the watershed. Both models also produce better performance when calibrated with SCA maps instead of point SWE data, likely because the SCA maps better constrain the space‐time pattern of SWE.     

A comparison between two statistical and a physically-based model in snow water equivalent mapping

Snow water equivalent (SWE) estimates at the end of the winter season have been compared for the 2002–2006 period in a 200 km² mountainous area in Switzerland, using three different models. The first model, ALPINE3D, is a physically based process-oriented model, which solves the snowpack energy and mass balance equations. The other two models, SWE-SEM and HS-SWE, are statistical algorithms interpolating snow data on a grid. While SWE-SEM interpolates local estimates of SWE, HS-SWE converts interpolated snow depth maps into maps of SWE using a regionally-calibrated conversion model. We discuss similarities and differences among the models’ results, both in terms of total volume, and spatial distribution of SWE. The comparison shows a general good agreement of the results of the three models, with a mean difference in the total volumes between the two statistical models of ∼8%, and between the physical model and the statistical ones of ∼−3% to −10%.     

MODIS Terra Collection 6 fractional snow cover validation in mountainous terrain during spring snowmelt using Landsat TM and ETM+

Daily swath MODIS Terra Collection 6 fractional snow cover (MOD10_L2) estimates were validated with two‐day Landsat TM/ETM + snow‐covered area estimates across central Idaho and southwestern Montana, USA. Snow cover maps during spring snowmelt for 2000, 2001, 2002, 2003, 2005, 2007, and 2009 were compared between MODIS Terra and Landsat TM/ETM + using least‐squared regression. Strong spatial and temporal map agreement was found between MODIS Terra fractional snow cover and Landsat TM/ETM + snow‐covered area, although map disagreement was observed for two validation dates. High‐altitude cirrus cloud contamination during low snow conditions as well as late season transient snowfall resulted in map disagreement. MODIS Terra's spatial resolution limits retrieval of thin‐patchy snow cover, especially during partially cloudy conditions. Landsat's image acquisition frequency can introduce difficulty when discriminating between transient and resident mountain snow cover. Furthermore, transient snowfall later in the snowmelt season, which is a stochastic accumulation event that does not usually persist beyond the daily timescale, will skew decadal snow‐covered area variability if bi‐monthly climate data record development is the objective. As a quality control step, ground‐based daily snow telemetry snow‐water‐equivalent measurements can be used to verify transient snowfall events. Users of daily MODIS Terra fractional snow products should be aware that local solar illumination and sensor viewing geometry might influence fractional snow cover estimation in mountainous terrain. Cross‐sensor interoperability has been confirmed between MODIS Terra and Landsat TM/ETM + when mapping snow from the visible/infrared spectrum. This relationship is strong and supports operational multi‐sensor snow cover mapping, specifically climate data record development to expand cryosphere, climate, and hydrological science applications. Copyright © 2013 John Wiley & Sons, Ltd.     

Melted snow volume control in the snowmelt runoff model using a snow water equivalent statistically based model

Snowmelt is an important component of the river discharge in mountain environments. In the past 40 years, the snowmelt dynamics has been mostly evaluated using degree‐day‐based models like the snowmelt runoff model (SRM). This model has no control on the volume of the melting snow, even if SRM includes as data input the snow‐covered area. This lack explains why the application of SRM may lead to inaccurate snowmelt volume estimations, even if the discharge volumes are accurately reproduced. Here we introduce in SRM the control on the melted snow volume and consider it in the determination of SRM parameters. The total snow volume, accumulated at the end of winter season, is evaluated by a snow water equivalent statistically based model, SWE‐SEM, and used as an estimate of the melting snow during the summer season. The benefit derived from the introduction of the control on the melting snow volume was investigated in the Mallero basin (northern Italy) for the 2003 and 2004 snow melting seasons. The analysis compares the model's results adopting different parameter sets, both considering and ignoring the control on the melting snow volume. Copyright © 2011 John Wiley & Sons, Ltd.     

Snow distribution in a steep mid-latitude alpine catchment

Snow distribution patterns are still poorly understood in steep alpine catchments because of mass redistribution from wind and avalanching. Snow models rarely operate with sufficient resolution, physics or input data to resolve this issue explicitly, and existing sub-grid parameterisations are rarely tested in this type of terrain. To address this issue daily snow cover observations, obtained from a ground-based camera, are combined with a snow melt model to estimate the spatial distribution of snow water equivalent (SWE) in a mountainous alpine catchment. Results show the importance of slope in controlling the spatial distribution of SWE and snow duration. This distribution is linked to the physical process of gravitational transport, where there is removal of snow from steep slopes and preferential deposition on moderate-angle slopes. From a modelling perspective, if sub-grid snow variability is parameterised using a log-normal probability distribution (as is common in hydrological and land-use models) then ignoring steep/shallow slope differences leads to an overestimation of melt at the beginning of the melt season, and a premature end to the snow melt season. When modelling SWE in complex terrain, care should be taken to consider reduced SWE on steep slopes.     

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.     

Role of temporal resolution of meteorological inputs for process‐based snow modelling

Accurate snow accumulation and melt simulations are crucial for understanding and predicting hydrological dynamics in mountainous settings. As snow models require temporally varying meteorological inputs, time resolution of these inputs is likely to play an important role on the model accuracy. Because meteorological data at a fine temporal resolution (~1 hr) are generally not available in many snow‐dominated settings, it is important to evaluate the role of meteorological inputs temporal resolution on the performance of process‐based snow models. The objective of this work is to assess the loss in model accuracy with temporal resolution of meteorological inputs, for a range of climatic conditions and topographic elevations. To this end, a process‐based snow model was run using 1‐, 3‐, and 6‐hourly inputs for wet, average, and dry years over Boise River Basin (6,963 km²), which spans rain dominated (≤1,400 m), rain–snow transition (>1,400 and ≤1,900 m), snow dominated below tree line (>1,900 and ≤2,400 m), and above tree line (>2,400 m) elevations. The results show that sensitivity of the model accuracy to the inputs time step generally decreases with increasing elevation from rain dominated to snow dominated above tree line. Using longer than hourly inputs causes substantial underestimation of snow cover area (SCA) and snow water equivalent (SWE) in rain‐dominated and rain–snow transition elevations, due to the precipitation phase mischaracterization. In snow‐dominated elevations, the melt rate is underestimated due to errors in estimation of net snow cover energy input. In addition, the errors in SCA and SWE estimates generally decrease toward years with low snow mass, that is, dry years. The results indicate significant increases in errors in estimates of SCA and SWE as the temporal resolution of meteorological inputs becomes coarser than an hour. However, use of 3‐hourly inputs can provide accurate estimates at snow‐dominated elevations. The study underscores the need to record meteorological variables at an hourly time step for accurate process‐based snow modelling.     

Estimating the distribution of snow water equivalent using remotely sensed snow cover data and a spatially distributed snowmelt model: A multi-resolution,multi-sensor comparison

Time series of fractional snow covered area (SCA) estimates from Landsat Enhanced Thematic Mapper (ETM+), Moderate Resolution Imaging Spectroradiometer (MODIS), and Advanced Very High Resolution Radiometer (AVHRR) data were combined with a spatially distributed snowmelt model to reconstruct snow water equivalent (SWE) in the Rio Grande headwaters (3419 km²). In this reconstruction approach, modeled snowmelt over each pixel is integrated during the period of satellite-observed snow cover to estimate SWE. Due to underestimates in snow cover detection, maximum basin-wide mean SWE using MODIS and AVHRR were, respectively, 45% and 68% lower than SWE estimates obtained using ETM+ data. The mean absolute error (MAE) of SWE estimated at 100-m resolution using ETM+ data was 23% relative to observed SWE from intensive field campaigns. Model performance deteriorated when MODIS (MAE = 50%) and AVHRR (MAE = 89%) SCA data were used. Relative to differences in the SCA products, model output was less sensitive to spatial resolution (MAE = 39% and 73% for ETM+ and MODIS simulations run at 1 km resolution, respectively), indicating that SWE reconstructions at the scale of MODIS acquisitions may be tractable provided the SCA product is improved. When considering tradeoffs between spatial and temporal resolution of different sensors, our results indicate that higher spatial resolution products such as ETM+ remain more accurate despite the lower frequency of acquisition. This motivates continued efforts to improve MODIS snow cover products.     

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.     

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.     

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.     

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.     

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.     

A comparison of MODIS and NOHRSC snow‐cover products for simulating streamflow using the Snowmelt Runoff Model

Remote sensing is an important source of snow‐cover extent for input into the Snowmelt Runoff Model (SRM) and other snowmelt models. Since February 2000, daily global snow‐cover maps have been produced from data collected by the Moderate Resolution Imaging Spectroradiometer (MODIS). The usefulness of this snow‐cover product for streamflow prediction is assessed by comparing SRM simulated streamflow using the MODIS snow‐cover product with streamflow simulated using snow maps from the National Operational Hydrologic Remote Sensing Center (NOHRSC). Simulations were conducted for two tributary watersheds of the Upper Rio Grande basin during the 2001 snowmelt season using representative SRM parameter values. Snow depletion curves developed from MODIS and NOHRSC snow maps were generally comparable in both watersheds: satisfactory streamflow simulations were obtained using both snow‐cover products in larger watershed (volume difference: MODIS, 2·6%; NOHRSC, 14·0%) and less satisfactory streamflow simulations in smaller watershed (volume difference: MODIS, −33·1%; NOHRSC, −18·6%). The snow water equivalent (SWE) on 1 April in the third zone of each basin was computed using the modified depletion curve produced by the SRM and was compared with in situ SWE measured at Snowpack Telemetry sites located in the third zone of each basin. The SRM‐calculated SWEs using both snow products agree with the measured SWEs in both watersheds. Based on these results, the MODIS snow‐cover product appears to be of sufficient quality for streamflow prediction using the SRM in the snowmelt‐dominated basins. Copyright © 2005 John Wiley & Sons, Ltd.