27 Oct 2021

27 Oct 2021

Review status: this preprint is currently under review for the journal TC.

Large-scale snow data assimilation using a spatialized particle filter: recovering the spatial structure of the particles

Jean Odry1, Marie-Amélie Boucher1, Simon Lachance-Cloutier2, Richard Turcotte2, and Pierre-Yves St-Louis2 Jean Odry et al.
  • 1Department of Civil and Building Engineering, Université de Sherbrooke, Sherbrooke, Canada
  • 2Quebec Ministère de l’Environnement et de la Lutte contre les Changements Climatiques, Quebec, Canada

Abstract. The use of particle filters for data assimilation is increasingly popular because of its minimal assumptions. Nevertheless, implementing a particle filter over domains of large spatial dimensions remains challenging, as the number of required particles rises exponentially as domain size increases. A common solution to overcome this issue is to localize the particle filter and consider a collection of local applications rather than a single regional one. Although this solution can solve the dimensionality limit, it can also create some spatial discontinuity inside the particles. This issue can become even more problematic when additional data is assimilated. The purpose of this study is to test the possibility of remedying the spatial discontinuities of the particles by locally reordering the particles.

We implement a spatialized particle filter to estimate the snow water equivalent (SWE) over a large territory in eastern Canada by assimilating local manual snow survey observations. We apply two reordering strategies based on 1) a simple ascending order sorting and 2) the Schaake Shuffle and evaluate their ability to maintain the spatial structure of the particles. To increase the amount of assimilated data, we investigate the inclusion of a second data set, in which SWE is indirectly estimated from snow depth. The two reordering solutions maintain the spatial structure of the individual particles throughout the winter season, which significantly reduces the random noise in the distribution of the particles and decreases the uncertainty associated with the estimation. The Schaake Shuffle proves to be a better tool for maintaining a realistic spatial structure for all particles, although we also found that sorting provides a simpler and satisfactory solution. The assimilation of the secondary data set improved SWE estimates in ungauged sites when compared with the open-loop model, but we noted no significant improvement when both snow courses and the SR50 data were assimilated.

Jean Odry et al.

Status: open (until 02 Jan 2022)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse

Jean Odry et al.

Data sets

Meteorological inputs Boucher, Marie-Amélie

HYDROTEL Snow Model Parameters Boucher, Marie-Amélie

Snow Observation Data Boucher, Marie-Amélie

Historical Snow Simulation (Open Loop) Boucher, Marie-Amélie

Model code and software

TheDroplets/Snow_spatial_particle_filter: First release of the codes for the spatial particle filte Odry, Jean

Jean Odry et al.


Total article views: 207 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
168 37 2 207 1 1
  • HTML: 168
  • PDF: 37
  • XML: 2
  • Total: 207
  • BibTeX: 1
  • EndNote: 1
Views and downloads (calculated since 27 Oct 2021)
Cumulative views and downloads (calculated since 27 Oct 2021)

Viewed (geographical distribution)

Total article views: 194 (including HTML, PDF, and XML) Thereof 194 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 26 Nov 2021
Short summary
The research deals with the assimilation of in-situ local snow observation in a large-scale spatialized snow modeling framework over the province of Quebec (eastern Canada). The methodology is based on proposing multiple spatialized snow scenarios using the snow model and weighting then according to the available observations. Especially, the paper focuses on the spatial coherence of the snow scenario proposed in the framework.