Snow water equivalent estimation in the Mallero basin using snow gauge data and MODIS images and fieldwork validation

Snow accumulation and ablation rule the temporal dynamics of water availability in mountain areas and cold regions. In these environments, the evaluation of the snow water amount is a key issue. The spatial distribution of snow water equivalent (SWE) over a mountain basin at the end of the snow accumulation season is estimated using a minimal statistical model (SWE‐SEM). This uses systematic observations such as ground measurements collected at snow gauges and snow‐covered area (SCA) data retrieved by remote sensors, here MODIS. Firstly, SWE‐SEM calculates local SWE estimates at snow gauges, then the spatial distribution of SWE over a certain area using an interpolation method; linear regressions of the first two order moments of SWE with altitude. The interpolation has been made by both confining and unconfining the spatial domain by SCA. SWE‐SEM is applied to the Mallero basin (northern Italy) for calculating the snow water equivalent at the end of the winter season for 6 years (2001–2007). For 2007, SWE‐SEM estimates are validated through fieldwork measurements collected during an ‘ad hoc’ campaign on March 31, 2007. Snow‐surveyed measurements are used to check SCA, snow density and SWE estimates. Copyright © 2009 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.     

Investigation of the snow‐cover dynamics in the Upper Euphrates Basin of Turkey using remotely sensed snow‐cover products and hydrometeorological data

The Euphrates and Tigris rivers serve as the most important water resources in the Middle East. Precipitation in this region falls mostly in the form of snow over the higher elevations of the Euphrates Basin and remains on the ground for nearly half of the year. This snow‐covered area (SCA) is a key element of the hydrological cycle, and monitoring the SCA is crucial for making accurate forecasts of snowmelt discharge, especially for energy production, flood control, irrigation, and reservoir‐operation optimization in the Upper Euphrates (Karasu) Basin. Remote sensing allows the detection of the spatio‐temporal patterns of snow cover across large areas in inaccessible terrain, such as the eastern part of Turkey, which is highly mountainous. In this study, a seasonal evaluation of the snow cover from 2000 to 2009 was performed using 8‐day snow‐cover products (MOD10C2) and the daily snow‐water equivalent (SWE) product. The values of SWE products were obtained using an assimilation process based on the Helsinki University of Technology model using equal area Special Sensor Microwave Imager (SSM/I) Earth‐gridded advanced microwave scanning radiometer—EOS daily brightness‐temperature values. In the Karasu Basin, the SCA percentage for the winter period is 80–90%. The relationship between the SCA and the runoff during the spring period is analysed for the period from 2004 to 2009. An inverse linear relationship between the normalized SCA and the normalized runoff values was obtained (r = 0·74). On the basis of the monthly mean temperature, total precipitation and snow depth observed at meteorological stations in the basin, the decrease in the peak discharges, and early occurrences of the peak discharges in 2008 and 2009 are due to the increase in the mean temperature and the decrease in the precipitation in April. 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.     

Probabilistic SWE reanalysis as a generalization of deterministic SWE reconstruction techniques

Snow accumulation and melt is highly variable in space and time in complex mountainous environments. Therefore, it is necessary to provide high‐resolution spatially and temporally distributed estimates of sub‐basin snow water equivalent (SWE) to accurately predict the timing and magnitude of snowmelt runoff. In this study, we compare two reconstruction techniques (a commonly used deterministic reconstruction vs a probabilistic data assimilation framework). The methods retrospectively estimate SWE from a time series of remotely sensed maps of fractional snow‐covered area (FSCA). In testing both methods over the Tokopah watershed in the Sierra Nevada (California), the probabilistic reconstruction approach is shown to be a more robust generalization of the deterministic reconstruction. Under idealized conditions, both probabilistic and deterministic approaches perform reasonably well and yield similar results when compared with in situ verification data, whereas the probabilistic reconstruction was found to be in slightly better agreement with snow‐pit observations. More importantly, the probabilistic approach was found to be more robust: unaccounted for biases in solar radiation impacted the probabilistic SWE estimates less than the deterministic case (4% vs 7% errors for water year (WY)1997 and 0% vs 3% errors for WY1999); the probabilistic reconstruction was found to be less sensitive to the number of available observations (6% vs 10% errors in WY1997 and 13% vs 44% errors in WY1999 from the nominal cases when four fewer FSCA images were available). Finally, results from the probabilistic reconstruction approach, which requires precipitation inputs (unlike the deterministic approach), were found to be relatively robust to bias in prior precipitation estimates, where the nominal case mean estimates were recovered even when an underestimated prior precipitation was used. The additional robustness of the probabilistic SWE reconstruction technique should prove useful in future applications over larger basins and longer periods in mountainous terrain. Copyright © 2013 John Wiley & Sons, Ltd.     

Determination of snow water equivalent over the eastern part of Turkey using passive microwave data

Information on regional snow water equivalent (SWE) is required for the management of water generated from snowmelt. Modeling of SWE in the mountainous regions of eastern Turkey, one of the major headwaters of Euphrates–Tigris basin, has significant importance in forecasting snowmelt discharge, especially for optimum water usage. An assimilation process to produce daily SWE maps is developed based on Helsinki University of Technology (HUT) model and AMSR‐E passive microwave data. The characteristics of the HUT emission model are analyzed in depth and discussed with respect to the extinction coefficient function. A new extinction coefficient function for the HUT model is proposed to suit models for snow over mountainous areas. Performance of the modified model is checked against the original, other modified cases and ground truth data covering the 2003–2007 winter periods. A new approach to calculate grain size and density is integrated inside the developed data assimilation process. An extensive validation was successfully performed by means of snow data measured at ground stations during the 2008–2010 winter periods. The root mean square error of the data set for snow depth and SWE between January and March of the 2008–2010 periods compared with the respective AMSR‐E footprints indicated that errors for estimated snow depth and predicted SWE values were 16.92 cm and 40.91 mm, respectively, for the 3‐year period. Validation results were less satisfactory for SWE less than 75.0 mm and greater than 150.0 mm. An underestimation for SWE greater than 150 mm could not be resolved owing to the microwave signal saturation that is observed for dense snowpack. Copyright © 2012 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.     

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.     

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.     

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.     

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.     

On the role of spatial resolution on snow estimates using a process‐based snow model across a range of climatology and elevation

Hydrological processes in mountainous settings depend on snow distribution, whose prediction accuracy is a function of model spatial scale. Although model accuracy is expected to improve with finer spatial resolution, an increase in resolution comes with modelling costs related to increased computational time and greater input data and parameter information. This computational and data collection expense is still a limiting factor for many large watersheds. Thus, this work's main objective is to question which physical processes lead to loss in model accuracy with regard to input spatial resolution under different climatic conditions and elevation ranges. To address this objective, a spatially distributed snow model, iSnobal, was run with inputs distributed at 50‐m—our benchmark for comparison—and 100‐m resolutions and with aggregated (averaged from the fine to the large resolution) inputs from the 50‐m model to 100‐, 250‐, 500‐, and 750‐m resolution for wet, average, and dry years over the Upper Boise River Basin (6,963 km²), which spans four elevation bands: rain dominated, rain–snow transition, and snow dominated below treeline and above treeline. Residuals, defined as differences between values quantified with high resolution (>50 m) models minus the benchmark model (50 m), of simulated snow‐covered area (SCA) and snow water equivalent (SWE) were generally slight in the aggregated scenarios. This was due to transferring the effects of topography on meteorological variables from the 50‐m model to the coarser scales through aggregation. Residuals in SCA and SWE in the distributed 100‐m simulation were greater than those of the aggregated 750 m. Topographic features such as slope and aspect were simplified, and their gradient was reduced due to coarsening the topography from the 50‐ to 100‐m resolution. Therefore, solar radiation was overestimated, and snow drifting was modified and caused substantial SCA and SWE underestimation in the distributed 100‐m model relative to the 50‐m model. Large residuals were observed in the wet year and at the highest elevation band when and where snow mass was large. These results support that model accuracy is substantially reduced with model scales coarser than 50 m.     

On the evaluation of snow water equivalent estimates over the terrestrial Arctic drainage basin

Comparisons between snow water equivalent (SWE) and river discharge estimates are important in evaluating the SWE fields and to our understanding of linkages in the freshwater cycle. In this study, we compared SWE drawn from land surface models and remote sensing observations with measured river discharge (Q) across 179 Arctic river basins. Over the period 1988‐2000, basin‐averaged SWE prior to snowmelt explains a relatively small (yet statistically significant) fraction of interannual variability in spring (April–June) Q, as assessed using the coefficient of determination (R²). Averaged across all basins, mean R²s vary from 0·20 to 0·28, with the best agreement noted for SWE drawn from a simulation with the Pan‐Arctic Water Balance Model (PWBM) forced with data from the European Centre for Medium‐Range Weather‐Forecasts (ECMWF) Re‐analysis (ERA‐40). Variability and magnitude in SWE derived from Special Sensor Microwave Imager (SSM/I) data are considerably lower than the variability and magnitude in SWE drawn from the land surface models, and generally poor agreement is noted between SSM/I SWE and spring Q. We find that the SWE versus Q comparisons are no better when alternate temporal integrations–using an estimate of the timing in basin thaw–are used to define pre‐melt SWE and spring Q. Thus, a majority of the variability in spring discharge must arise from factors other than basin snowpack water storage. This study demonstrates how SWE estimated from remote sensing observations, or general circulation models (GCMs), can be evaluated effectively using monthly discharge data or SWE from a hydrological model. Copyright © 2007 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.     

Intercomparison of snow water equivalent observations in the Northern Great Plains

In the Northern Great Plains, melting snow is a primary driver of spring flooding, but limited knowledge of the magnitude and spatial distribution of snow water equivalent (SWE) hampers flood forecasting. Passive microwave remote sensing has the potential to enhance operational river flow forecasting but is not routinely incorporated in operational flood forecasting. We compare satellite passive microwave estimates from the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR‐E) to the National Oceanic and Atmospheric Administration Office of Water Prediction (OWP) airborne gamma radiation snow survey and U.S. Army Corps of Engineers (USACE) ground snow survey SWE estimates in the Northern Great Plains from 2002 to 2011. AMSR‐E SWE estimates compare favourably with USACE SWE measurements in the low relief, low vegetation study area (mean difference = ?3.8 mm, root mean squared difference [RMSD] = 34.7 mm), but less so with OWP airborne gamma SWE estimates (mean difference = ?9.5 mm, RMSD = 42.7 mm). An error simulation suggests that up to half of the error in the former comparison is potentially due to subpixel scale SWE variability, limiting the maximum achievable RMSD between ground and satellite SWE to approximately 26–33 mm in the Northern Great Plains. The OWP gamma versus AMSR‐E SWE comparison yields larger error than the point‐scale USACE versus AMSR‐E comparison, despite a larger measurement footprint (5–7 km² vs. a few square centimetres, respectively), suggesting that there are unshared errors between the USACE and OWP gamma SWE data.     

Towards improved snow water equivalent retrieval algorithms for satellite passive microwave data over the mountainous basins of western USA

Space‐borne passive microwave snow water equivalent (SWE) retrieval algorithms are attractive for continuous SWE monitoring over large mountainous areas. The performance of three SWE retrieval algorithms, which were considered relevant for operational purposes, was examined for each month over the Colorado River Basin. In addition, statistical post‐processing was tested as a means of improving the SWE estimates from each algorithm. The evaluation started with the so‐called Chang equation, which was a pioneer algorithm and is still used in practice. Successive attempts were then made to improve the algorithm's performance through the calibration of the equation's coefficient and through the inclusion of brightness temperature data from various frequency channels. The Chang equation consistently underestimated SWE with average bias between 30 mm in November and more than 300 mm in April and root mean square error (RMSE) exceeding 500 mm at many locations in April. The statistical post‐processing effectively removed the bias and reduced the RMSE by half for all the months. When the Chang equation's coefficients were calibrated at each site, biases were reduced by approximately 85%, and RMSE was reduced by 40%–50%. Finally, the multiple channel equations produced unbiased SWE estimates with RMSEs 50%–60% of those from the Chang equation. However, the statistical post‐processing did not reduce RMSE for both calibrated algorithms. The last algorithm produced the most reliable estimates for at‐site analysis, but its skill deteriorated when analyses were performed over larger areal extents; therefore, it is only recommended for SWE monitoring over smaller areas. For larger areas, the calibrated Chang equation is desirable because it only requires interpolations of a calibrated coefficient, which was spatially coherent. Copyright © 2011 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.     

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.     

Evaluation of the snow‐covered area data product from MODIS

The Moderate Resolution Imaging Spectroradiometer (MODIS), flown on board the Terra Earth Observing System (EOS) platform launched in December 1999, produces a snow‐covered area (SCA) product. This product is expected to be of better quality than SCA products based on operational satellites (notably GOES and AVHRR), due both to improved spectral resolution and higher spatial resolution of the MODIS instrument. The gridded MODIS SCA product was compared with the SCA product produced and distributed by the National Weather Service National Operational Hydrologic Remote Sensing Center (NOHRSC) for 46 selected days over the Columbia River basin and 32 days over the Missouri River basin during winter and spring of 2000–01. Snow presence or absence was inferred from ground observations of snow depth at 1330 stations in the Missouri River basin and 762 stations in the Columbia River basin, and was compared with the presence/absence classification for the corresponding pixels in the MODIS and NOHRSC SCA products. On average, the MODIS SCA images classified fewer pixels as cloud than NOHRSC, the effect of which was that 15% more of the Columbia basin area could be classified as to presence–absence of snow, while overall there was a statistically insignificant difference over the Missouri basin. Of the pixels classified as cloud free, MODIS misclassified 4% and 5% fewer overall (for the Columbia and Missouri basins respectively) than did the NOHRSC product. When segregated by vegetation cover, forested areas had the greatest differences in fraction of cloud cover reported by the two SCA products, with MODIS classifying 13% and 17% less of the images as cloud for the Missouri and Columbia basins respectively. These differences are particularly important in the Columbia River basin, 39% of which is forested. The ability of MODIS to classify significantly greater amounts of snow in the presence of cloud in more topographically complex, forested, and snow‐dominated areas of these two basins provides valuable information for hydrologic prediction. Copyright © 2003 John Wiley & Sons, Ltd.     

Calibration of snow parameters in SWAT: comparison of three approaches in the Upper Adige River basin (Italy)

The Soil and Water Assessment Tool (SWAT) model is generally applied in alpine catchments using a unique set of snow parameters for the entire basin, and calibration is based on discharge records only. This technical note presents three calibration procedures for snow parameters of SWAT considering snow water equivalent (SWE) values computed using a dense network of snow depth measurement stations available in the Upper Adige River basin, Italy. The first two procedures calibrate snow parameters according to the average sub-basin SWE: the first one defines a unique set of parameters for the entire basin, while the second allows for sub-basin variability. The last approach includes the elevation band SWE output in the calibration for each sub-basin and qualitatively compares it to the SWE computed from the available snow depth monitoring stations. This last method provides the best agreement between SWAT model results and SWE data.