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.     

Quantifying snow water equivalent using terrestrial ground penetrating radar and unmanned aerial vehicle photogrammetry

This study demonstrates the potential value of a combined unmanned aerial vehicle (UAV) Photogrammetry and ground penetrating radar (GPR) approach to map snow water equivalent (SWE) over large scales. SWE estimation requires two different physical parameters (snow depth and density), which are currently difficult to measure with the spatial and temporal resolution desired for basin-wide studies. UAV photogrammetry can provide very high-resolution spatially continuous snow depths (SD) at the basin scale, but does not measure snow densities. GPR allows nondestructive quantitative snow investigation if the radar velocity is known. Using photogrammetric snow depths and GPR two-way travel times (TWT) of reflections at the snow-ground interface, radar velocities in snowpack can be determined. Snow density (RSN) is then estimated from the radar propagation velocity (which is related to electrical permittivity of snow) via empirical formulas. A Phantom-4 Pro UAV and a MALA GX450 HDR model GPR mounted on a ski mobile were used to determine snow parameters. A snow-free digital surface model (DSM) was obtained from the photogrammetric survey conducted in September 2017. Then, another survey in synchronization with a GPR survey was conducted in February 2019 whilst the snowpack was approximately at its maximum thickness. Spatially continuous snow depths were calculated by subtracting the snow-free DSM from the snow-covered DSM. Radar velocities in the snowpack along GPR survey lines were computed by using UAV-based snow depths and GPR reflections to obtain snow densities and SWEs. The root mean square error of the obtained SWEs (384 mm average) is 63 mm, indicating good agreement with independent SWE observations and the error lies within acceptable uncertainty limits.     

Evaluation of a water equivalent of snow cover map in a small catchment area using a geostatistical approach

We describe an objective method for evaluating the spatial distribution of water equivalents of the snow cover within a small catchment. Regression analysis is used to quantify the relationship between elevation, presence or absence of forest, and potential direct solar radiation as independent variables and water equivalent as measured at a number of sites. First, this regression relationship is used to interpolate water equivalent data all over the basin area. Then we interpolate the residuals of the regression using a geostatistical approach. Superimposing the results obtained by interpolating the regression relationship and the interpolated residuals eventually yields the water equivalent distribution over the test area. The advantages of the interpolation method used lie in the optimal (effective, unbiased) estimation of the interpolated values as well as in the possibility to quantify the associated estimation variances.     

Performance characteristics of a new electronic snow water equivalent sensor in different climates

The US Army ERDC CRREL and the US Department of Agriculture Natural Resources Conservation Service developed a square electronic snow water equivalent (e‐SWE) sensor as an alternative to using fluid‐filled snow pillows to measure SWE. The sensors consist of a centre panel to measure SWE and eight outer panels to buffer edge stress concentrations. Seven 3 m square e‐SWE sensors were installed in five different climate zones. During the 2011–2012 winter, 1.8 and 1.2 m square e‐SWE sensors were installed and operated in Oregon. With the exception of New York State and Newfoundland, the e‐SWE sensors accurately measured SWE, with R² values between the sensor and manual SWE measurements of between 0.86 and 0.98. The e‐SWE sensor at Hogg Pass, Oregon, accurately measured SWE during the past 8 years of operations. In the thin, icy snow of New York during midwinter 2008–2009, the e‐SWE sensors overmeasured SWE because of edge stress concentrations associated with strong icy layers and a shallow snow cover. The New York e‐SWE sensors' measurement accuracy improved in spring 2009 and further improved during the 2011–2012 winter with operating experience. At Santiam Junction, measured SWE from the 1.8 and 1.2 m square e‐SWE sensors agreed well with the snow pillow, 3 m square e‐SWE sensor, and manual SWE measurements until February 2013, when dust and gravel blew onto the testing area resulting in anomalous measurements. © 2014 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.     

Assessing the impact of distributed snow water equivalent calibration and assimilation of Copernicus snow water equivalent on modelled snow and streamflow performance

Accuracy of the Copernicus snow water equivalent (SWE) product and the impact of SWE calibration and assimilation on modelled SWE and streamflow was evaluated. Daily snowpack measurements were made at 12 locations from 2016 to 2019 across a 4104 km² mixed-forest basin in the Great Lakes region of central Ontario, Canada. Sub-basin daily SWE calculated from these sites, observed discharge, and lake levels were used to calibrate a hydrologic model developed using the Raven modelling framework. Copernicus SWE was bias corrected during the melt period using mean bias subtraction and was compared to daily basin average SWE calculated from the measured data. Bias corrected Copernicus SWE was assimilated into the models using a range of parameters and the parameterizations from the model calibration. The bias corrected Copernicus product agreed well with measured data and provided a good estimate of mean basin SWE demonstrating that the product shows promise for hydrology applications within the study region. Calibration to spatially distributed SWE substantially improved the basin scale SWE estimate while only slightly degrading the flow simulation demonstrating the value of including SWE in a multi-objective calibration formulation. The particle filter experiments yielded the best SWE estimation but moderately degraded the flow simulation. The particle filter experiments constrained by the calibrated snow parameters produced similar results to the experiments using the upper and lower bounds indicating that, in this study, model calibration prior to assimilation was not valuable. The calibrated models exhibited varying levels of skill in estimating SWE but demonstrated similar streamflow performance. This indicates that basin outlet streamflow can be accurately estimated using a model with a poor representation of distributed SWE. This may be sufficient for applications where estimating flow is the primary water management objective. However, in applications where understanding the physical processes of snow accumulation, melt and streamflow generation are important, such as assessing the impact of climate change on water resources, accurate representations of SWE are required and can be improved via multi-objective calibration or data assimilation, as demonstrated in this study.     

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

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

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.     

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.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Hydrological dynamics of snowmelt induced streamflow in a high mountain catchment of the Pyrenees under contrasting snow accumulation and duration years

Snowmelt drives a large portion of streamflow in many mountain areas of the world. However, the water paths from snowmelt to the arrival of the water in the streams are still largely unknown. This work analyzes for first time the influence of snowmelt on spring streamflow with different snow accumulation and duration, in an alpine catchment of the central Spanish Pyrenees. This study presents the water balance of the main melting months (May and June). Piezometric values, water temperature, electrical conductivity and isotope data (δ¹⁸O) allow a better understanding of the hydrological functioning of the basin during these months. Results of the water balance calculations showed that snow represented on average 73% of the water available for streamflow in May and June while precipitation during these months accounted for only 27%. However, rainfall during the melting period was important to determine the shape of the spring hydrographs. On average, 78% of the sum of both the snow water equivalent (SWE) accumulated at the beginning of May and the precipitation in May and June converted into runoff during the May–June melting period. The average evaporation-sublimation during the 2 months corresponded to 8.4% of the accumulated SWE and rainfall, so that only a small part of the water input was ultimately available for soil and groundwater storage. When snow cover disappeared from the catchment, soil water storage and streamflow showed a sharp decline. Consequently, streamflow electrical conductivity, temperature and δ¹⁸O showed a marked tipping point towards higher values. The fast hydrological response of the catchment to snow and meteorological fluctuations, as well as the marked diel fluctuations of streamflow δ¹⁸O during the melting period, strongly suggests short meltwater transit times. As a consequence of this hydrological behaviour, independently of the amount of snow accumulated and of melting date, summer streamflow remained always low, with only small runoff peaks driven by rainfall events.     

Synoptic sampling reveals spatial stability in flow and chemistry patterns despite large seasonal variability in a subarctic mountain catchment

Seasonality plays a critical role in cold mountain regions as variation in air temperature, ground thermal status, and precipitation phase alter the rate, timing and magnitude of hydrological and chemical transport. Additionally, cold mountain catchments can have highly variable topography, geology, permafrost, and landcover, which intrinsically add to this irregularity. Understanding how external and internal variability act to control mass fluxes requires sampling at a high spatial resolution over time, which rarely occurs in cold remote regions. In this work, we conduct five snapshot sampling surveys across 34 subcatchments during the ice-free period in Wolf Creek Research Basin (a mesoscale montane subarctic catchment) and two additional winter surveys across a subset of sites to assess the drivers of variability in stream chemistry and discharge. We sampled for specific conductance (SpC), major ions, and dissolved organic carbon (DOC) and used statistical metrics and Bayesian mixing analysis to quantify patterns of flow and chemistry across space and time. Our results indicate patterns in both flow and chemistry remain largely consistent across seasons for all solutes. However, there was weaker correlation of chemistry between sites, suggesting asynchronous behaviour within the catchment. There was evidence of increasing production of ions and DOC along the stream network during high spring flows but not during low flows. Although concentrations and flows exhibit high seasonality in subarctic mountains, this seasonal variability does not alter spatial patterns that arise from highly variable catchment characteristics.     

Snow accumulation and ablation response to changes in forest structure and snow surface albedo after attack by mountain pine beetle

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

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.     

Comparison of five snow water equivalent estimation methods across categories

Snow water equivalent (SWE) is an important indicator used in hydrology, water resources, and climate change impact. There are various methods of estimating SWE (falling in 3 categories: indirect sensors, empirical models, and process‐based models), but few studies that provide comparison across these different categories to help users make decisions on monitoring site design or method selection. Five SWE estimation methods were compared against manual snow course data collected over 2 years (2015–2016) from the Dorset Environmental Science Centre, including the gamma‐radiation‐based CS725 sensor, 3 empirical estimation models (Sexstone snow density model, McCreight & Small snow density model, and a meteorology‐based model), and the University of British Columbia Watershed Model snow energy‐balance model. Snow depth, density, and SWE were measured at the Dorset Environmental Science Centre weather station in south‐central Ontario, on a daily basis over 6 winters from 2011 to 2016. The 2 snow density‐based models, requiring daily snow depth as input, gave the best performance (R² of .92 and .92 for McCreight & Small and Sexstone models, respectively). The CS725 sensor that receives radiation coming from soil penetrating the snowpack provided the same performance (R² = .92), proving that the sensor is an applicable method, although it is expensive. The meteorology‐based empirical model, requiring daily climate data including temperature, precipitation and solar radiation, gave the poorest performance (R² = .77). The energy‐balance‐based University of British Columbia Watershed Model snow module, only requiring climate data, worked better than the empirical meteorology‐based model (R² = .9) but performed worse than the density models or CS725 sensor. Given differences in application objectives, site conditions, and budget, this comparison across SWE estimation methods may help users choose a suitable method. For ongoing and new monitoring sites, installation of a CS725 sensor coupled with intermittent manual snow course measurements (e.g., weekly) is recommended for further SWE method estimation testing and development of a snow density model.     

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.     

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.     

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.     

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.     

Analysing the contribution of snow water equivalent to the terrestrial water storage over Canada

In this study, the spatial and temporal variabilities of terrestrial water storage anomaly (TWSA) and snow water equivalent anomaly (SWEA) information obtained from the Gravity Recovery and Climate Experiment (GRACE) twin satellites data were analysed in conjunction with multisource snow products over several basins in the Canadian landmass. Snow water equivalent (SWE) data were extracted from three different sources: Global Snow Monitoring for Climate Research version 2 (GlobSnow2), Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E), and Canadian Meteorological Centre (CMC). The objective of the study was to understand whether SWE variations have a significant contribution to terrestrial water storage anomalies in the Canadian landmass. The period was considered from December 2002 to March 2011. Significant relationships were observed between TWSA and SWEA for most of the 15 basins considered (53% to 80% of the basins, depending on the SWE products considered). The best results were obtained with the CMC SWE products compared with satellite-based SWE data. Stronger relationships were found in snow-dominated basins (R_s > = 0.7), such as the Liard [root mean square error (RMSE) = 21.4 mm] and Peace Basins (RMSE = 26.76 mm). However, despite high snow accumulation in the north of Quebec, GRACE showed weak or insignificant correlations with SWEA, regardless of the data sources. The same behaviour was observed in the Western Hudson Bay basin. In both regions, it was found that the contribution of non-SWE compartments including wetland, surface water, as well as soil water storages has a significant impact on the variations of total storage. These components were estimated using the Water-Global Assessment and Prognosis Global Hydrology Model (WGHM) simulations and then subtracted from GRACE observations. The GRACE-derived SWEA correlation results showed improved relationships with three SWEA products. The improvement is particularly important in the sub-basins of the Hudson Bay, where very weak and insignificant results were previously found with GRACE TWSA data. GRACE-derived SWEA showed a significant relationship with CMC data in 93% of the basins (13% more than GRACE TWSA). Overall, the results indicated the important role of SWE on terrestrial water storage variations.