scholarly journals Estimating snow water equivalent in a Sierra Nevada watershed via spaceborne radiance data assimilation

2017 ◽  
Vol 53 (1) ◽  
pp. 647-671 ◽  
Author(s):  
Dongyue Li ◽  
Michael Durand ◽  
Steven A. Margulis
2015 ◽  
Vol 16 (4) ◽  
pp. 1752-1772 ◽  
Author(s):  
Steven A. Margulis ◽  
Manuela Girotto ◽  
Gonzalo Cortés ◽  
Michael Durand

Abstract This paper presents a newly proposed data assimilation method for historical snow water equivalent SWE estimation using remotely sensed fractional snow-covered area fSCA. The newly proposed approach consists of a particle batch smoother (PBS), which is compared to a previously applied Kalman-based ensemble batch smoother (EnBS) approach. The methods were applied over the 27-yr Landsat 5 record at snow pillow and snow course in situ verification sites in the American River basin in the Sierra Nevada (United States). This basin is more densely vegetated and thus more challenging for SWE estimation than the previous applications of the EnBS. Both data assimilation methods provided significant improvement over the prior (modeling only) estimates, with both able to significantly reduce prior SWE biases. The prior RMSE values at the snow pillow and snow course sites were reduced by 68%–82% and 60%–68%, respectively, when applying the data assimilation methods. This result is encouraging for a basin like the American where the moderate to high forest cover will necessarily obscure more of the snow-covered ground surface than in previously examined, less-vegetated basins. The PBS generally outperformed the EnBS: for snow pillows the PBS RMSE was ~54% of that seen in the EnBS, while for snow courses the PBS RMSE was ~79% of the EnBS. Sensitivity tests show relative insensitivity for both the PBS and EnBS results to ensemble size and fSCA measurement error, but a higher sensitivity for the EnBS to the mean prior precipitation input, especially in the case where significant prior biases exist.


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


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.


2020 ◽  
Author(s):  
Bertrand Cluzet ◽  
Matthieu Lafaysse ◽  
Emmanuel Cosme ◽  
Clément Albergel ◽  
Louis-François Meunier ◽  
...  

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.


2018 ◽  
Author(s):  
Daniel Abel ◽  
Felix Pollinger ◽  
Heiko Paeth

Abstract. Droughts can result in enormous impacts for environment, societies, and economy. In arid or semiarid regions with bordering high mountains, snow is the major source of water supply due to its role as natural water storage. The goal of this study is to examine the influence of snow water equivalent (SWE) on droughts in the United States and find large-scale climatic predictors for SWE and drought. For this, a Maximum Covariance Analysis (MCA), also known as Singular Value Decomposition, is performed with snow data from the ERA–Interim reanalysis and the self-calibrating Palmer Drought Severity Index (sc–PDSI) as drought index. Furthermore, the relationship of resulting principal components and original data with atmospheric patterns is investigated. The leading mode shows the spatial connection between SWE and drought via downstream water/moisture transport. Especially the Rocky Mountains in Colorado (CR) play a key role for the central and western South, but the Sierra Nevada and even the Appalachian Mountains are relevant, too. The temperature and precipitation based sc–PDSI is able to capture this link because increased soil moisture results in higher evapotranspiration with lower sensible heat and vice versa. A time shifted MCA indicates a prediction skill for drought conditions in spring and early summer for the downstream regions of CR on the basis of SWE in March. Furthermore, the phase of the El Niño–Southern Oscillation is a good predictor for drought in the southern US and SWE around Colorado. The influence of the North Atlantic Oscillation and Pacific North American Pattern is not that clear.


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.


2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Yun Xu ◽  
Andrew Jones ◽  
Alan Rhoades

Abstract The simulation of snow water equivalent (SWE) remains difficult for regional climate models. Accurate SWE simulation depends on complex interacting climate processes such as the intensity and distribution of precipitation, rain-snow partitioning, and radiative fluxes. To identify the driving forces behind SWE difference between model and reanalysis datasets, and guide model improvement, we design a framework to quantitatively decompose the SWE difference contributed from precipitation distribution and magnitude, ablation, temperature and topography biases in regional climate models. We apply this framework within the California Sierra Nevada to four regional climate models from the North American Coordinated Regional Downscaling Experiment (NA-CORDEX) run at three spatial resolutions. Models generally predict less SWE compared to Landsat-Era Sierra Nevada Snow Reanalysis (SNSR) dataset. Unresolved topography associated with model resolution contribute to dry and warm biases in models. Refining resolution from 0.44° to 0.11° improves SWE simulation by 35%. To varying degrees across models, additional difference arises from spatial and elevational distribution of precipitation, cold biases revealed by topographic correction, uncertainties in the rain-snow partitioning threshold, and high ablation biases. This work reveals both positive and negative contributions to snow bias in climate models and provides guidance for future model development to enhance SWE simulation.


2012 ◽  
Vol 13 (6) ◽  
pp. 1970-1976 ◽  
Author(s):  
Jonathan D. D. Meyer ◽  
Jiming Jin ◽  
Shih-Yu Wang

Abstract The authors investigated the accuracy of snow water equivalent (SWE) observations compiled by 748 Snowpack Telemetry stations and attributed the systematic bias introduced to SWE measurements to drifting snow. Often observed, SWE outpaces accumulated precipitation (AP), which can be statistically and physically explained through 1) precipitation undercatchment and/or 2) drifting snow. Forty-four percent of the 748 stations reported at least one year where the maximum SWE was greater than AP, while 16% of the stations showed this inconsistency for at least 20% of the observed years. Regions with a higher likelihood of inconsistency contained drier snow and are exposed to higher winds speeds, both of which are positively correlated to drifting snow potential as well as gauge undercatch. Differentiating between gauge undercatch and potential drifting scenarios, days when SWE increased but AP remained zero were used. These drift days occurred on an average of 13.3 days per year for all stations, with 31% greater wind speeds at 10 m for such days (using reanalysis winds). Findings suggest marked consistency between SWE and AP throughout the Cascade Mountains and lower elevations of the interior west while indicating notable inconsistency between these two variables throughout the higher elevations of the Rocky Mountains, Utah mountain ranges, and the Sierra Nevada.


2012 ◽  
Vol 16 (3) ◽  
pp. 815-831 ◽  
Author(s):  
M. He ◽  
T. S. Hogue ◽  
S. A. Margulis ◽  
K. J. Franz

Abstract. The current study proposes an integrated uncertainty and ensemble-based data assimilation framework (ICEA) and evaluates its viability in providing operational streamflow predictions via assimilating snow water equivalent (SWE) data. This step-wise framework applies a parameter uncertainty analysis algorithm (ISURF) to identify the uncertainty structure of sensitive model parameters, which is subsequently formulated into an Ensemble Kalman Filter (EnKF) to generate updated snow states for streamflow prediction. The framework is coupled to the US National Weather Service (NWS) snow and rainfall-runoff models. Its applicability is demonstrated for an operational basin of a western River Forecast Center (RFC) of the NWS. Performance of the framework is evaluated against existing operational baseline (RFC predictions), the stand-alone ISURF and the stand-alone EnKF. Results indicate that the ensemble-mean prediction of ICEA considerably outperforms predictions from the other three scenarios investigated, particularly in the context of predicting high flows (top 5th percentile). The ICEA streamflow ensemble predictions capture the variability of the observed streamflow well, however the ensemble is not wide enough to consistently contain the range of streamflow observations in the study basin. Our findings indicate that the ICEA has the potential to supplement the current operational (deterministic) forecasting method in terms of providing improved single-valued (e.g., ensemble mean) streamflow predictions as well as meaningful ensemble predictions.


Sign in / Sign up

Export Citation Format

Share Document