scholarly journals Combined assimilation of streamflow and snow water equivalent for mid-term ensemble streamflow forecasts in snow-dominated regions

2016 ◽  
Author(s):  
Jean M. Bergeron ◽  
Mélanie Trudel ◽  
Robert Leconte
Keyword(s):  
Data Assimilation ◽  
State Vector ◽  
Snow Water Equivalent ◽  
Synthetic Data ◽  
Hydrologic Model ◽  
Open Loop ◽  
The State ◽  
Wide Range ◽  
Snow Water ◽  
The Impact

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.

scholarly journals Combined assimilation of streamflow and snow water equivalent for mid-term ensemble streamflow forecasts in snow-dominated regions

2016 ◽  
Vol 20 (10) ◽  
pp. 4375-4389 ◽  
Author(s):  
Jean M. Bergeron ◽  
Mélanie Trudel ◽  
Robert Leconte
Keyword(s):  
Data Assimilation ◽  
State Vector ◽  
Snow Water Equivalent ◽  
Synthetic Data ◽  
Hydrologic Model ◽  
Open Loop ◽  
The State ◽  
Wide Range ◽  
Snow Water ◽  
The Impact

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.

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

2021 ◽  
Author(s):  
Ilaria Clemenzi ◽  
David Gustafsson ◽  
Jie Zhang ◽  
Björn Norell ◽  
Wolf Marchand ◽  
...  
Keyword(s):  
Remote Sensing ◽  
Spatial Distribution ◽  
Data Assimilation ◽  
Snow Water Equivalent ◽  
Distribution Model ◽  
Accurate Estimation ◽  
Snowmelt Runoff ◽  
Snow Water ◽  
Snow Distribution ◽  
The Impact

<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>

scholarly journals Examining cross-scale influences of forcing resolutions in a hillslope-resolving, integrated hydrologic model

2020 ◽  
Author(s):  
Miguel A. Aguayo ◽  
Alejandro N. Flores ◽  
James P. McNamara ◽  
Hans-Peter Marshall ◽  
Jodi Mead
Keyword(s):  
Hydrologic Modeling ◽  
Snow Water Equivalent ◽  
Spatial Scales ◽  
Water Storage ◽  
Hydrologic Model ◽  
Atmospheric Forcing ◽  
Hydrologic Models ◽  
Snow Water ◽  
Integrated Hydrologic Model ◽  
The Impact

Abstract. Water management in semiarid regions of the western United States requires accurate and timely knowledge of runoff generated by snowmelt. This information is used to plan reservoir releases for downstream users and hydrologic models play an important role in estimating the volume of snow stored in mountain watersheds that serve as source waters for downstream reservoirs. Physically based, integrated hydrologic models are used to develop spatiotemporally dynamic estimates of hydrologic states and fluxes based on understanding of the underlying biophysics of hydrologic response. Yet this class of models are associated with many issues that give rise to significant uncertainties in key hydrologic variables of interest like snow water storage and streamflow. Underlying sources of uncertainty include difficulties in parameterizing processes associated with nonlinearities of some processes, as well as from the large variability in the characteristic spatial and temporal scale of atmospheric forcing and land-surface water and energy balance and groundwater processes. Scale issues, in particular, can introduce systematic biases in integrated atmospheric and hydrologic modeling. Reconciling these discrepancies while maintaining computational tractability remains a fundamental challenge in integrated hydrologic modeling. Here we investigate the hydrologic impact of discrepancies between distributed meteorological forcing data exhibiting a range of spatial scales consistent with a variety of numerical weather prediction models when used to force an integrated hydrologic model associated with a corresponding range of spatial resolutions characteristic of distributed hydrologic modeling. To achieve this, we design and conduct a total of twelve numerical modeling experiments that seek to quantify the impact of applied resolution of atmospheric forcings on simulated hillslope-scale hydrologic state variables. The experiments are arranged in such way to assess the impact of four different atmospheric forcing resolutions (i.e., interpolated 30 m, 1 km, 3 km and 9 km) on two hydrologic variables, snow water equivalent and soil water storage, arranged in three hydrologic spatial resolution (i.e., 30 m, 90 m and 250 m). Results show spatial patterns in snow water equivalent driven by atmospheric forcing in hillslope-scale simulations and patterns mostly driven by topographical characteristics (i.e., slope and aspect) on coarser simulations. Similar patterns are observed in soil water storage however, in addition to that, large errors are encountered primarily in riparian areas of the watershed on coarser simulations. The Weather Research Forecasting (WRF) model is used to develop the environmental forcing variables required as input to the integrated hydrologic model. WRF is an open source, community supported coupled land-atmosphere model capable of capturing spatial scales that permit convection. The integrated hydrologic modeling framework used in this work coincides with the ParFlow open-source surface-subsurface hydrology model. This work has important implications for the use of atmospheric and integrated hydrologic models in remote and ungauged areas. In particular, this work has potential ramifications for the design and development of observing system simulation experiments (OSSEs) in complex and snow-dominated landscapes. OSSEs are critical in constraining the performance characteristics of Earth-observing satellites.

scholarly journals Data Assimilation of Terrestrial Water Storage Observations to Estimate Precipitation Fluxes: A Synthetic Experiment

2021 ◽  
Vol 13 (6) ◽  
pp. 1223
Author(s):  
Manuela Girotto ◽  
Rolf Reichle ◽  
Matthew Rodell ◽  
Viviana Maggioni
Keyword(s):  
Data Assimilation ◽  
Snow Water Equivalent ◽  
Water Storage ◽  
Synthetic Data ◽  
Terrestrial Water Storage ◽  
Grace Data ◽  
Precipitation Estimates ◽  
Terrestrial Water ◽  
Snow Water ◽  
Gravity Recovery

The Gravity Recovery and Climate Experiment (GRACE) mission and its Follow-On (GRACE-FO) mission provide unprecedented observations of terrestrial water storage (TWS) dynamics at basin to continental scales. Established GRACE data assimilation techniques directly adjust the simulated water storage components to improve the estimation of groundwater, streamflow, and snow water equivalent. Such techniques artificially add/subtract water to/from prognostic variables, thus upsetting the simulated water balance. To overcome this limitation, we propose and test an alternative assimilation scheme in which precipitation fluxes are adjusted to achieve the desired changes in simulated TWS. Using a synthetic data assimilation experiment, we show that the scheme improves performance skill in precipitation estimates in general, but that it is more robust for snowfall than for rainfall, and it fails in certain regions with strong horizontal gradients in precipitation. The results demonstrate that assimilation of TWS observations can help correct (adjust) the model’s precipitation forcing and, in turn, enhance model estimates of TWS, snow mass, soil moisture, runoff, and evaporation. A key limitation of the approach is the assumption that all errors in TWS originate from errors in precipitation. Nevertheless, the proposed approach produces more consistent improvements in simulated runoff than the established GRACE data assimilation techniques.

A Simple Data Assimilation System for Complex Snow Distributions (SnowAssim)

2008 ◽  
Vol 9 (5) ◽  
pp. 989-1004 ◽  
Author(s):  
Glen E. Liston ◽  
Christopher A. Hiemstra
Keyword(s):  
Data Assimilation ◽  
Snow Water Equivalent ◽  
Synthetic Data ◽  
Optimal Interpolation ◽  
Data Assimilation Scheme ◽  
Simple Simulation ◽  
Simulation Domain ◽  
Snow Water ◽  
Evolution Modeling ◽  
Assimilation Scheme

Abstract A methodology for assimilating ground-based and remotely sensed snow data within a snow-evolution modeling system (SnowModel) is presented. The data assimilation scheme (SnowAssim) is consistent with optimal interpolation approaches in which the differences between the observed and modeled snow values are used to constrain modeled outputs. The calculated corrections are applied retroactively to create improved fields prior to the assimilated observations. Thus, one of the values of this scheme is the improved simulation of snow-related distributions throughout the entire snow season, even when observations are only available late in the accumulation and/or ablation periods. Because of this, the technique is particularly applicable to reanalysis applications. The methodology includes the ability to stratify the assimilation into regions where either the observations and/or model has unique error properties, such as the differences between forested and nonforested snow environments. The methodologies are introduced using synthetic data and a simple simulation domain. In addition, the model is applied over NASA’s Cold Land Processes Experiment (CLPX), Rabbit Ears Pass, Colorado, observation domain. Simulations using the data assimilation scheme were found to improve the modeled snow water equivalent (SWE) distributions, and simulated SWE displayed considerably more realistic spatial heterogeneity than that provided by the observations alone.

scholarly journals Generating Observation-Based Snow Depletion Curves for Use in Snow Cover Data Assimilation

Geosciences ◽  
2018 ◽  
Vol 8 (12) ◽  
pp. 484 ◽  
Author(s):  
Kristi Arsenault ◽  
Paul Houser
Keyword(s):  
Data Assimilation ◽  
Snow Cover ◽  
Land Surface ◽  
Snow Depth ◽  
Snow Water Equivalent ◽  
Land Surface Model ◽  
Open Loop ◽  
Surface Model ◽  
Snow Water ◽  
Logarithmic Functions

Snow depletion curves (SDC) are functions that are used to show the relationship between snow covered area and snow depth or water equivalent. Previous snow cover data assimilation (DA) studies have used theoretical SDC models as observation operators to map snow depth to snow cover fraction (SCF). In this study, a new approach is introduced that uses snow water equivalent (SWE) observations and satellite-based SCF retrievals to derive SDC relationships for use in an Ensemble Kalman filter (EnKF) to assimilate snow cover estimates. A histogram analysis is used to bin the SWE observations, which the corresponding SCF observations are then averaged within, helping to constrain the amount of data dispersion across different temporal and regional conditions. Logarithmic functions are linearly regressed with the binned average values, for two U.S. mountainous states: Colorado and Washington. The SDC-based logarithmic functions are used as EnKF observation operators, and the satellite-based SCF estimates are assimilated into a land surface model. Assimilating satellite-based SCF estimates with the observation-based SDC shows a reduction in SWE-related RMSE values compared to the model-based SDC functions. In addition, observation-based SDC functions were derived for different intra-annual and physiographic conditions, and landcover and elevation bands. Lower SWE-based RMSE values are also found with many of these categorical observation-based SDC EnKF experiments. All assimilation experiments perform better than the open-loop runs, except for the Washington region’s 2004–2005 snow season, which was a major drought year that was difficult to capture with the ensembles and observations.

scholarly journals A High-Resolution Data Assimilation Framework for Snow Water Equivalent Estimation across the Western United States and Validation with the Airborne Snow Observatory

2019 ◽  
Vol 20 (3) ◽  
pp. 357-378 ◽  
Author(s):  
Catalina M. Oaida ◽  
John T. Reager ◽  
Konstantinos M. Andreadis ◽  
Cédric H. David ◽  
Steve R. Levoe ◽  
...  
Keyword(s):  
United States ◽  
High Resolution ◽  
Data Assimilation ◽  
Snow Water Equivalent ◽  
High Elevation ◽  
Open Loop ◽  
Western United States ◽  
Olympic Peninsula ◽  
Depletion Curve ◽  
Snow Water

Abstract Numerical simulations of snow water equivalent (SWE) in mountain systems can be biased, and few SWE observations have existed over large domains. New approaches for measuring SWE, like NASA’s ultra-high-resolution Airborne Snow Observatory (ASO), offer an opportunity to improve model estimates by providing a high-quality validation target. In this study, a computationally efficient snow data assimilation (DA) approach over the western United States at 1.75-km spatial resolution for water years (WYs) 2001–17 is presented. A local ensemble transform Kalman filter implemented as a batch smoother is used with the VIC hydrology model to assimilate the remotely sensed daily MODIS fractional snow-covered area (SCA). Validation of the high-resolution SWE estimates is done against ASO SWE data in the Tuolumne basin (California), Uncompahgre basin (Colorado), and Olympic Peninsula (Washington). Results indicate good performance in dry years and during melt, with DA reducing Tuolumne basin-average SWE percent differences from −68%, −92%, and −84% in open loop to 0.6%, 25%, and 3% after DA for WYs 2013–15, respectively, for ASO dates and spatial extent. DA also improved SWE percent difference over the Uncompahgre basin (−84% open loop, −65% DA) and Olympic Peninsula (26% open loop, −0.2% DA). However, in anomalously wet years DA underestimates SWE, likely due to an inadequate snow depletion curve parameterization. Despite potential shortcomings due to VIC model setup (e.g., water balance mode) or parameterization (snow depletion curve), the DA framework implemented in this study shows promise in overcoming some of these limitations and improving estimated SWE, in particular during drier years or at higher elevations, when most in situ observations cannot capture high-elevation snowpack due to lack of stations there.

scholarly journals Effect of snow microstructure variability on Ku-band radar snow water equivalent retrievals

2019 ◽  
Vol 13 (11) ◽  
pp. 3045-3059 ◽  
Author(s):  
Nick Rutter ◽  
Melody J. Sandells ◽  
Chris Derksen ◽  
Joshua King ◽  
Peter Toose ◽  
...  
Keyword(s):  
Spatial Variability ◽  
Layer Thickness ◽  
Snow Water Equivalent ◽  
Estimation Accuracy ◽  
Snow Layer ◽  
Data Set ◽  
Snow Water ◽  
Retrieval Errors ◽  
Snow Microstructure ◽  
The Impact

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.

scholarly journals Assimilation of snow cover and snow depth into a snow model to estimate snow water equivalent and snowmelt runoff in a Himalayan catchment

2017 ◽  
Vol 11 (4) ◽  
pp. 1647-1664 ◽  
Author(s):  
Emmy E. Stigter ◽  
Niko Wanders ◽  
Tuomo M. Saloranta ◽  
Joseph M. Shea ◽  
Marc F. P. Bierkens ◽  
...  
Keyword(s):  
Kalman Filter ◽  
Data Assimilation ◽  
Snow Cover ◽  
Climate Sensitivity ◽  
Snow Depth ◽  
Snow Water Equivalent ◽  
Snowmelt Runoff ◽  
Snow Water ◽  
Snow Model

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.

scholarly journals CrocO_v1.0 : a Particle Filter to assimilate snowpack observations in a spatialised framework

2020 ◽  
Author(s):  
Bertrand Cluzet ◽  
Matthieu Lafaysse ◽  
Emmanuel Cosme ◽  
Clément Albergel ◽  
Louis-François Meunier ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Particle Filter ◽  
Near Infrared ◽  
Snow Water Equivalent ◽  
Sources Of Information ◽  
Average Improvement ◽  
Modelling Uncertainties ◽  
Snow Water ◽  
Spatio Temporal ◽  
Snowpack Properties

Abstract. Monitoring the evolution of the snowpack properties in mountainous areas is crucial for avalanche hazard forecasting and water resources management. In-situ and remotely sensed observations provide precious information on the snowpack but usually offer a limited spatio-temporal coverage of bulk or surface variables only. In particular, visible-near infrared (VIS-NIR) reflectance observations can inform on the snowpack surface properties but are limited by terrain shading and clouds. Snowpack modelling enables to estimate any physical variable, virtually anywhere, but is affected by large errors and uncertainties. Data assimilation offers a way to combine both sources of information, and to propagate information from observed areas to non observed areas. Here, we present CrocO, (Crocus-Observations) an ensemble data assimilation system able to ingest any snowpack observation (applied as a first step to the height of snow (HS) and VIS-NIR reflectances) in a spatialised geometry. CrocO uses an ensemble of snowpack simulations to represent modelling uncertainties, and a Particle Filter (PF) to reduce them. The PF is known to collapse when assimilating a too large number of observations. Two variants of the PF were specifically implemented to ensure that observations information is propagated in space while tackling this issue. The global algorithm ingests all available observations with an iterative inflation of observation errors, while the klocal algorithm is a localised approach performing a selection of the observations to assimilate based on background correlation patterns. Experiments are carried out in a twin experiment setup, with synthetic observations of HS and VIS-NIR reflectances available in only a 1/6th of the simulation domain. Results show that compared against runs without assimilation, analyses exhibit an average improvement of snow water equivalent Continuous Rank Probability Score (CRPS) of 60 % when assimilating HS with a 40-member ensemble, and an average 20 % CRPS improvement when assimilating reflectance with a 160-member ensemble. Significant improvements are also obtained outside the observation domain. These promising results open a way for the assimilation of real observations of reflectance, or of any snowpack observations in a spatialised context.

Sign in / Sign up

Export Citation Format

Share Document