Exploring the use of gpr and satellite-based snow data for snowmelt runoff predictions

Author(s):  
Ilaria Clemenzi ◽  
David Gustafsson ◽  
Jie Zhang ◽  
Björn Norell ◽  
Wolf Marchand ◽  
...  

<p>Snow in the mountains is the result of the interplay between meteorological conditions, e.g., precipitation, wind and solar radiation, and landscape features, e.g., vegetation and topography. For this reason, it is highly variable in time and space. It represents an important water storage for several sectors of the society including tourism, ecology and hydropower. The estimation of the amount of snow stored in winter and available in the form of snowmelt runoff can be strategic for their sustainability. In the hydropower sector, for example, the occurrence of higher snow and snowmelt runoff volumes at the end of the spring and in the early summer compared to the estimated one can substantially impact reservoir regulation with energy and economical losses. An accurate estimation of the snow volumes and their spatial and temporal distribution is thus essential for spring flood runoff prediction. Despite the increasing effort in the development of new acquisition techniques, the availability of extensive and representative snow and density measurements for snow water equivalent estimations is still limited. Hydrological models in combination with data assimilation of ground or remote sensing observations is a way to overcome these limitations. However, the impact of using different types of snow observations on snowmelt runoff predictions is, little understood. In this study we investigated the potential of assimilating in situ and remote sensing snow observations to improve snow water equivalent estimates and snowmelt runoff predictions. We modelled the seasonal snow water equivalent distribution in the Lake Överuman catchment, Northern Sweden, which is used for hydropower production. Simulations were performed using the semi-distributed hydrological model HYPE for the snow seasons 2017-2020. For this purpose, a snowfall distribution model based on wind-shelter factors was included to represent snow spatial distribution within model units. The units consist of 2.5x2.5 km<sup>2</sup> grid cells, which were further divided into hydrological response units based on elevation, vegetation and aspect. The impact on the estimation of the total catchment mean snow water equivalent and snowmelt runoff volume were evaluated using for data assimilation, gpr-based snow water equivalent data acquired along survey lines in the catchment in the early spring of the four years, snow water equivalent data obtained by a machine learning algorithm and satellite-based fractional snow cover data. Results show that the wind-shelter based snow distribution model was able to represent a similar spatial distribution as the gpr survey lines, when assessed on the catchment level. Deviations in the model performance within and between specific gpr survey lines indicate issues with the spatial distribution of input precipitation, and/or need to include explicit representation of snow drift between model units. The explicit snow distribution model also improved runoff simulations, and the ability of the model to improve forecast through data assimilation.</p>

2016 ◽  
Author(s):  
Thomas Skaugen ◽  
Ingunn H. Weltzien

Abstract. Snow is an important and complicated element in hydrological modelling. The traditional catchment hydrological model with its many free calibration parameters, also in snow sub-models, is not a well-suited tool for predicting conditions for which it has not been calibrated. Such conditions include prediction in ungauged basins and assessing hydrological effects of climate change. In this study, a new model for the spatial distribution of snow water equivalent (SWE), parameterized solely from observed spatial variability of precipitation (SD_G), is compared with the current snow distribution model used in the operational flood forecasting models in Norway. The latter model (SD_LN) has a fixed, calibrated coefficient of variation, which parameterizes a log-normal model for snow distribution. The two models are implemented in the already parameter parsimonious rainfall runoff model Distance Distribution Dynamics (DDD) and their capability for predicting runoff, SWE and snow covered area (SCA) are tested and compared for 71 Norwegian catchments. Results show that SD_G better simulates SCA when compared with MODIS satellite derived snow cover. In addition, SWE is simulated more realistically in that seasonal snow is melted out and the building up of "snow towers" and giving spurious positive trends in SWE, typical for SD_LN, is prevented. The precision of runoff simulations using SD_G is slightly inferior, with a reduction in Nash-Sutcliffe and Kling Gupta Criterion of 0.01, but it is shown that high precision in runoff prediction using SD_LN is accompanied with erroneous simulations of SWE.


2016 ◽  
Vol 10 (5) ◽  
pp. 1947-1963 ◽  
Author(s):  
Thomas Skaugen ◽  
Ingunn H. Weltzien

Abstract. Snow is an important and complicated element in hydrological modelling. The traditional catchment hydrological model with its many free calibration parameters, also in snow sub-models, is not a well-suited tool for predicting conditions for which it has not been calibrated. Such conditions include prediction in ungauged basins and assessing hydrological effects of climate change. In this study, a new model for the spatial distribution of snow water equivalent (SWE), parameterized solely from observed spatial variability of precipitation, is compared with the current snow distribution model used in the operational flood forecasting models in Norway. The former model uses a dynamic gamma distribution and is called Snow Distribution_Gamma, (SD_G), whereas the latter model has a fixed, calibrated coefficient of variation, which parameterizes a log-normal model for snow distribution and is called Snow Distribution_Log-Normal (SD_LN). The two models are implemented in the parameter parsimonious rainfall–runoff model Distance Distribution Dynamics (DDD), and their capability for predicting runoff, SWE and snow-covered area (SCA) is tested and compared for 71 Norwegian catchments. The calibration period is 1985–2000 and validation period is 2000–2014. Results show that SD_G better simulates SCA when compared with MODIS satellite-derived snow cover. In addition, SWE is simulated more realistically in that seasonal snow is melted out and the building up of "snow towers" and giving spurious positive trends in SWE, typical for SD_LN, is prevented. The precision of runoff simulations using SD_G is slightly inferior, with a reduction in Nash–Sutcliffe and Kling–Gupta efficiency criterion of 0.01, but it is shown that the high precision in runoff prediction using SD_LN is accompanied with erroneous simulations of SWE.


2017 ◽  
Vol 11 (4) ◽  
pp. 1647-1664 ◽  
Author(s):  
Emmy E. Stigter ◽  
Niko Wanders ◽  
Tuomo M. Saloranta ◽  
Joseph M. Shea ◽  
Marc F. P. Bierkens ◽  
...  

Abstract. Snow is an important component of water storage in the Himalayas. Previous snowmelt studies in the Himalayas have predominantly relied on remotely sensed snow cover. However, snow cover data provide no direct information on the actual amount of water stored in a snowpack, i.e., the snow water equivalent (SWE). Therefore, in this study remotely sensed snow cover was combined with in situ observations and a modified version of the seNorge snow model to estimate (climate sensitivity of) SWE and snowmelt runoff in the Langtang catchment in Nepal. Snow cover data from Landsat 8 and the MOD10A2 snow cover product were validated with in situ snow cover observations provided by surface temperature and snow depth measurements resulting in classification accuracies of 85.7 and 83.1 % respectively. Optimal model parameter values were obtained through data assimilation of MOD10A2 snow maps and snow depth measurements using an ensemble Kalman filter (EnKF). Independent validations of simulated snow depth and snow cover with observations show improvement after data assimilation compared to simulations without data assimilation. The approach of modeling snow depth in a Kalman filter framework allows for data-constrained estimation of snow depth rather than snow cover alone, and this has great potential for future studies in complex terrain, especially in the Himalayas. Climate sensitivity tests with the optimized snow model revealed that snowmelt runoff increases in winter and the early melt season (December to May) and decreases during the late melt season (June to September) as a result of the earlier onset of snowmelt due to increasing temperature. At high elevation a decrease in SWE due to higher air temperature is (partly) compensated by an increase in precipitation, which emphasizes the need for accurate predictions on the changes in the spatial distribution of precipitation along with changes in temperature.


2016 ◽  
Author(s):  
Jean M. Bergeron ◽  
Mélanie Trudel ◽  
Robert Leconte

Abstract. The potential of data assimilation for hydrologic predictions has been demonstrated in many research studies. Watersheds over which multiple observation types are available can potentially further benefit from data assimilation by having multiple updated states from which hydrologic predictions can be generated. However, the magnitude and time span of the impact of the assimilation of an observation varies according not only to its type, but also to the variables included in the state vector. This study examines the impact of multivariate synthetic data assimilation using the Ensemble Kalman Filter (EnKF) into the spatially distributed hydrologic model CEQUEAU for the mountainous Nechako River located in British-Columbia, Canada. Synthetic data includes daily snow cover area (SCA), daily measurements of snow water equivalent (SWE) at three different locations and daily streamflow data at the watershed outlet. Results show a large variability of the continuous rank probability skill score over a wide range of prediction horizons (days to weeks) depending on the state vector configuration and the type of observations assimilated. Overall, the variables most closely linearly linked to the observations are the ones worth considering adding to the state vector. The performance of the assimilation of basin-wide SCA, which does not have a decent proxy among potential state variables, does not surpass the open loop for any of the simulated variables. However, the assimilation of streamflow offers major improvements steadily throughout the year, but mainly over the short-term (up to 5 days) forecast horizons, while the impact of the assimilation of SWE gains more importance during the snowmelt period over the mid-term (up to 50 days) forecast horizon compared with open loop. The combined assimilation of streamflow and SWE performs better than its individual counterparts, offering improvements over all forecast horizons considered and throughout the whole year, including the critical period of snowmelt. This highlights the potential benefit of using multivariate data assimilation for streamflow predictions in snow-dominated regions.


2016 ◽  
Vol 20 (10) ◽  
pp. 4375-4389 ◽  
Author(s):  
Jean M. Bergeron ◽  
Mélanie Trudel ◽  
Robert Leconte

Abstract. The potential of data assimilation for hydrologic predictions has been demonstrated in many research studies. Watersheds over which multiple observation types are available can potentially further benefit from data assimilation by having multiple updated states from which hydrologic predictions can be generated. However, the magnitude and time span of the impact of the assimilation of an observation varies according not only to its type, but also to the variables included in the state vector. This study examines the impact of multivariate synthetic data assimilation using the ensemble Kalman filter (EnKF) into the spatially distributed hydrologic model CEQUEAU for the mountainous Nechako River located in British Columbia, Canada. Synthetic data include daily snow cover area (SCA), daily measurements of snow water equivalent (SWE) at three different locations and daily streamflow data at the watershed outlet. Results show a large variability of the continuous rank probability skill score over a wide range of prediction horizons (days to weeks) depending on the state vector configuration and the type of observations assimilated. Overall, the variables most closely linearly linked to the observations are the ones worth considering adding to the state vector due to the limitations imposed by the EnKF. The performance of the assimilation of basin-wide SCA, which does not have a decent proxy among potential state variables, does not surpass the open loop for any of the simulated variables. However, the assimilation of streamflow offers major improvements steadily throughout the year, but mainly over the short-term (up to 5 days) forecast horizons, while the impact of the assimilation of SWE gains more importance during the snowmelt period over the mid-term (up to 50 days) forecast horizon compared with open loop. The combined assimilation of streamflow and SWE performs better than their individual counterparts, offering improvements over all forecast horizons considered and throughout the whole year, including the critical period of snowmelt. This highlights the potential benefit of using multivariate data assimilation for streamflow predictions in snow-dominated regions.


Geografie ◽  
2014 ◽  
Vol 119 (3) ◽  
pp. 199-217 ◽  
Author(s):  
Dana Kučerová ◽  
Michal Jeníček

The knowledge of the water volume stored in the snowpack, including its spatial distribution, is vital for many hydrological applications. Such information is useful for hydrological forecasts and it is often used for the calibration of snowmelt runoff models. Data from four field measurements of the snow water equivalent (SWE) carried out in two winter seasons were assessed by ten interpolation methods. Measurements from both snow accumulation and snowmelt periods were evaluated. The ability of methods to predict SWE at unmeasured locations was assessed by the means of cross validation. The best prediction accuracy of SWE was achieved by means of multiple a simple linear regressions, residual kriging and cokriging methods. The accuracy was enhanced by the use of elevation, aspect, slope and vegetation as variables in the calculation of the SWE. Elevation and vegetation show a significant correlation with the SWE in the study area. The multiple regression gave best results for snow accumulation period. However, the spatial variability of SWE was not successfully explained for snowmelt periods.


2006 ◽  
Vol 7 (3) ◽  
pp. 443-457 ◽  
Author(s):  
Michael Durand ◽  
Steven A. Margulis

Abstract A season-long, point-scale radiometric data assimilation experiment is performed in order to test the feasibility of snow water equivalent (SWE) estimation. Synthetic passive microwave observations at Special Sensor Microwave Imager (SSM/I) and Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) frequencies and synthetic broadband albedo observations are assimilated simultaneously in order to update snowpack states in a land surface model using the ensemble Kalman filter (EnKF). The effects of vegetation and atmosphere are included in the radiative transfer model (RTM). The land surface model (LSM) was given biased precipitation to represent typical errors introduced in modeling, yet was still able to recover the true value of SWE with a seasonally integrated rmse of only 2.95 cm, despite a snow depth of around 3 m and the presence of liquid water in the snowpack. This ensemble approach is ideal for investigating the complex theoretical relationships between the snowpack properties and the observations, and exploring the implications of these relationships for the inversion of remote sensing measurements for estimating snowpack properties. The contributions of each channel to recovering the true SWE are computed, and it was found that the low-frequency 10.67-GHz AMSR-E channels contain information even for very deep snow. The effect of vegetation thickness on assimilation results is explored. Results from the assimilation are compared to those from a pure modeling approach and from a remote sensing inversion approach, and the effects of measurement error and ensemble size are investigated.


2004 ◽  
Vol 38 ◽  
pp. 180-186 ◽  
Author(s):  
Thomas Skaugen ◽  
Eli Alfnes ◽  
Elin G. Langsholt ◽  
Hans-Christian Udnæs

AbstractA model is put forward which focuses on the dynamical evolution of the spatial distribution of snow water equivalent (SWE). We make use of the fact that when the accumulation and ablation process of the snow reservoir is modelled as a summation of a gamma-distributed variable, both skewed distributions, typical of alpine areas, and more normal distributions, typical of forested areas, can be accounted for. A particular problem is to represent fractional snow-covered area (SCA) within the distribution framework. The change in SCA as a response to a melting event is explicitly linked to the shape of the distribution of SWE and is estimated as the probability of non-exceedance of the melted amount from a scaled version of the spatial distribution of SWE. An extensive snow-measuring programme, where several snow courses have been measured repeatedly throughout the melting season, justifies the dynamical aspects of the snow distribution in the modelling approach. The modelling approach has been tested with the Swedish rainfall–runoff model, HBV, and estimated values of SWE and SCA are compared with results obtained using the statistical distribution (log-normal) traditionally used in the model.


2019 ◽  
Vol 13 (11) ◽  
pp. 3045-3059 ◽  
Author(s):  
Nick Rutter ◽  
Melody J. Sandells ◽  
Chris Derksen ◽  
Joshua King ◽  
Peter Toose ◽  
...  

Abstract. Spatial variability in snowpack properties negatively impacts our capacity to make direct measurements of snow water equivalent (SWE) using satellites. A comprehensive data set of snow microstructure (94 profiles at 36 sites) and snow layer thickness (9000 vertical profiles across nine trenches) collected over two winters at Trail Valley Creek, NWT, Canada, was applied in synthetic radiative transfer experiments. This allowed for robust assessment of the impact of estimation accuracy of unknown snow microstructural characteristics on the viability of SWE retrievals. Depth hoar layer thickness varied over the shortest horizontal distances, controlled by subnivean vegetation and topography, while variability in total snowpack thickness approximated that of wind slab layers. Mean horizontal correlation lengths of layer thickness were less than a metre for all layers. Depth hoar was consistently ∼30 % of total depth, and with increasing total depth the proportion of wind slab increased at the expense of the decreasing surface snow layer. Distinct differences were evident between distributions of layer properties; a single median value represented density and specific surface area (SSA) of each layer well. Spatial variability in microstructure of depth hoar layers dominated SWE retrieval errors. A depth hoar SSA estimate of around 7 % under the median value was needed to accurately retrieve SWE. In shallow snowpacks <0.6 m, depth hoar SSA estimates of ±5 %–10 % around the optimal retrieval SSA allowed SWE retrievals within a tolerance of ±30 mm. Where snowpacks were deeper than ∼30 cm, accurate values of representative SSA for depth hoar became critical as retrieval errors were exceeded if the median depth hoar SSA was applied.


Sign in / Sign up

Export Citation Format

Share Document