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

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.

2021 ◽  
Vol 14 (3) ◽  
pp. 1595-1614
Author(s):  
Bertrand Cluzet ◽  
Matthieu Lafaysse ◽  
Emmanuel Cosme ◽  
Clément Albergel ◽  
Louis-François Meunier ◽  
...  

Abstract. Monitoring the evolution of 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 state of the snowpack but usually offer limited spatio-temporal coverage of bulk or surface variables only. In particular, visible–near-infrared (Vis–NIR) reflectance observations can provide information about the snowpack surface properties but are limited by terrain shading and clouds. Snowpack modelling enables the estimation of any physical variable virtually anywhere, but it 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 prone to collapse when assimilating too many observations. Two variants of the PF were specifically implemented to ensure that observational 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. Feasibility testing experiments are carried out in an identical twin experiment setup, with synthetic observations of HS and Vis–NIR reflectances available in only one-sixth of the simulation domain. Results show that compared against runs without assimilation, analyses exhibit an average improvement of the 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 possibility for the assimilation of real observations of reflectance or of any snowpack observations in a spatialised context.


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.


2019 ◽  
Author(s):  
Joel Fiddes ◽  
Kristoffer Aalstad ◽  
Sebastian Westermann

Abstract. Spatial variability in high-relief landscapes is immense, and grid-based models cannot be run at spatial resolutions to explicitly represent important physical processes. This hampers the assessment of the current and future evolution of important issues such as water availability or mass movement hazards. Here, we present a new processing chain that couples an efficient subgrid method with a downscaling tool and data assimilation method with the purpose to improve numerical simulation of surface processes at multiple spatial and temporal scales in ungauged basins. The novelty of the approach is that while we add 1–2 orders of magnitude of computational cost by ensemble simulations, we save 4–5 orders of magnitude over explicitly simulating a high-resolution grid. This approach makes data assimilation at large spatio-temporal scales feasible. In addition, this approach utilises only freely available global datasets and is therefore able to run globally. We demonstrate marked improvements in estimating snow height and snow water equivalent at various experimental scales using this approach. We propose this as a suitable method for a wide variety of operational and research applications where surface models need to be run at large scales with sparse to non-existent ground observations and with the flexibility to assimilate diverse variables retrieved by EO missions.


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.


2012 ◽  
Vol 4 (1) ◽  
pp. 13-21 ◽  
Author(s):  
S. Morin ◽  
Y. Lejeune ◽  
B. Lesaffre ◽  
J.-M. Panel ◽  
D. Poncet ◽  
...  

Abstract. A quality-controlled snow and meteorological dataset spanning the period 1 August 1993–31 July 2011 is presented, originating from the experimental station Col de Porte (1325 m altitude, Chartreuse range, France). Emphasis is placed on meteorological data relevant to the observation and modelling of the seasonal snowpack. In-situ driving data, at the hourly resolution, consist of measurements of air temperature, relative humidity, windspeed, incoming short-wave and long-wave radiation, precipitation rate partitioned between snow- and rainfall, with a focus on the snow-dominated season. Meteorological data for the three summer months (generally from 10 June to 20 September), when the continuity of the field record is not warranted, are taken from a local meteorological reanalysis (SAFRAN), in order to provide a continuous and consistent gap-free record. Data relevant to snowpack properties are provided at the daily (snow depth, snow water equivalent, runoff and albedo) and hourly (snow depth, albedo, runoff, surface temperature, soil temperature) time resolution. Internal snowpack information is provided from weekly manual snowpit observations (mostly consisting in penetration resistance, snow type, snow temperature and density profiles) and from a hourly record of temperature and height of vertically free ''settling'' disks. This dataset has been partially used in the past to assist in developing snowpack models and is presented here comprehensively for the purpose of multi-year model performance assessment. The data is placed on the PANGAEA repository (http://dx.doi.org/10.1594/PANGAEA.774249) as well as on the public ftp server ftp://ftp-cnrm.meteo.fr/pub-cencdp/.


2020 ◽  
Vol 21 (11) ◽  
pp. 2713-2733 ◽  
Author(s):  
Graham A. Sexstone ◽  
Colin A. Penn ◽  
Glen E. Liston ◽  
Kelly E. Gleason ◽  
C. David Moeser ◽  
...  

AbstractThis study evaluated the spatial variability of trends in simulated snowpack properties across the Rio Grande headwaters of Colorado using the SnowModel snow evolution modeling system. SnowModel simulations were performed using a grid resolution of 100 m and 3-hourly time step over a 34-yr period (1984–2017). Atmospheric forcing was provided by phase 2 of the North American Land Data Assimilation System, and the simulations accounted for temporal changes in forest canopy from bark beetle and wildfire disturbances. Annual summary values of simulated snowpack properties [snow metrics; e.g., peak snow water equivalent (SWE), snowmelt rate and timing, and snow sublimation] were used to compute trends across the domain. Trends in simulated snow metrics varied depending on elevation, aspect, and land cover. Statistically significant trends did not occur evenly within the basin, and some areas were more sensitive than others. In addition, there were distinct trend differences between the different snow metrics. Upward trends in mean winter air temperature were 0.3°C decade−1, and downward trends in winter precipitation were −52 mm decade−1. Middle elevation zones, coincident with the greatest volumetric snow water storage, exhibited the greatest sensitivity to changes in peak SWE and snowmelt rate. Across the Rio Grande headwaters, snowmelt rates decreased by 20% decade−1, peak SWE decreased by 14% decade−1, and total snowmelt quantity decreased by 13% decade−1. These snow trends are in general agreement with widespread snow declines that have been reported for this region. This study further quantifies these snow declines and provides trend information for additional snow variables across a greater spatial coverage at finer spatial resolution.


2015 ◽  
Vol 61 (226) ◽  
pp. 267-279 ◽  
Author(s):  
Joshua King ◽  
Richard Kelly ◽  
Andrew Kasurak ◽  
Claude Duguay ◽  
Grant Gunn ◽  
...  

AbstractDuring the 2010/11 boreal winter, a distributed set of backscatter measurements was collected using a ground-based Ku-band (17.2 GHz) scatterometer system at 26 open tundra sites. A standard snow-sampling procedure was completed after each scan to evaluate local variability in snow layering, depth, density and water equivalent (SWE) within the scatterometer field of view. The shallow depths and large basal depth hoar encountered presented an opportunity to evaluate backscatter under a set of previously untested conditions. Strong Ku-band response was found with increasing snow depth and snow water equivalent (SWE). In particular, co-polarized vertical backscatter increased by 0.82 dB for every 1 cm increase in SWE (R2 = 0.62). While the result indicated strong potential for Ku-band retrieval of shallow snow properties, it did not characterize the influence of sub-scan variability. An enhanced snow-sampling procedure was introduced to generate detailed characterizations of stratigraphy within the scatterometer field of view using near-infrared photography along the length of a 5 m trench. Changes in snow properties along the trench were used to discuss variations in the collocated backscatter response. A pair of contrasting observation sites was used to highlight uncertainties in backscatter response related to short length scale spatial variability in the observed tundra environment.


Sign in / Sign up

Export Citation Format

Share Document