CRAMPON: A Particle Filter to assimilate sparse snowpack observations into a semi-distributed geometry

Author(s):  
Bertrand Cluzet ◽  
Matthieu Lafaysse ◽  
Marie Dumont ◽  
Emmanuel Cosme ◽  
Clément Albergel

<p>In mountainous areas, detailed snowpack models are essential to capture the high spatio-temporal variability of the snowpack. This task is highly challenging, and models suffer from large simulation errors. In these regions, in-situ observations are scarce, while remote sensing observations are generally patchy owing to complex physiographic features (steep slopes, forests, shadows,...) and weather conditions (clouds). This point is stressing the need for a spatially coherent data assimilation system able to propagate the informations into unobserved locations.</p><p>In this study, we present CRAMPON (CRocus with AssiMilation of snowPack ObservatioNs), an ensemble data assimilation system ingesting snowpack observations in a spatialized context. CRAMPON quantifies snowpack modelling uncertainties with an ensemble and reduces them using a Particle Filter. Stochastic perturbations of meteorological forcings and the multi-physical version of Crocus snowpack model (ESCROC) are used to build the ensemble. Two variants of the Sequential Importance Resampling Particle Filter (PF) were implemented to tackle the common PF degeneracy issue that arises when assimilating a large number of observations. In a first approach (so-called global approach), the observations information is spread across topographic conditions by looking for a global analysis. Degeneracy is mitigated by inflating the observation error covariance matrix, with the side effect of reducing the impact of the assimilation. In a second approach (klocal), we propagate the information and mitigate degeneracy by a localisation of the PF based on background correlation patterns between topographic conditions.</p><p>Here, we investigate the ability of CRAMPON to globally benefit from partial observations in a conceptual semi-distributed domain which accounts for the main features of topographic-induced snowpack variability. We compare simulations without assimilation with experiments assimilating synthetic observations of the Height of Snow and VIS/NIR reflectance. This setup demonstrates the ability of CRAMPON to spread the information of various snow observations into unobserved locations.</p>

1990 ◽  
Vol 118 (12) ◽  
pp. 2513-2542 ◽  
Author(s):  
Ross N. Hoffman ◽  
Christopher Grassotti ◽  
Ronald G. Isaacs ◽  
Jean-Francois Louis ◽  
Thomas Nehrkorn ◽  
...  

2020 ◽  
Vol 21 (9) ◽  
pp. 2023-2039
Author(s):  
Dikra Khedhaouiria ◽  
Stéphane Bélair ◽  
Vincent Fortin ◽  
Guy Roy ◽  
Franck Lespinas

AbstractConsistent and continuous fields provided by precipitation analyses are valuable for hydrometeorological applications and land data assimilation modeling, among others. Providing uncertainty estimates is a logical step in the analysis development, and a consistent approach to reach this objective is the production of an ensemble analysis. In the present study, a 6-h High-Resolution Ensemble Precipitation Analysis (HREPA) was developed for the domain covering Canada and the northern part of the contiguous United States. The data assimilation system is the same as the Canadian Precipitation Analysis (CaPA) and is based on optimal interpolation (OI). Precipitation from the Canadian national 2.5-km atmospheric prediction system constitutes the background field of the analysis, while at-site records and radar quantitative precipitation estimates (QPE) compose the observation datasets. By using stochastic perturbations, multiple observations and background field random realizations were generated to subsequently feed the data assimilation system and provide 24 HREPA members plus one control run. Based on one summer and one winter experiment, HREPA capabilities in terms of bias and skill were verified against at-site observations for different climatic regions. The results indicated HREPA’s reliability and skill for almost all types of precipitation events in winter, and for precipitation of medium intensity in summer. For both seasons, HREPA displayed resolution and sharpness. The overall good performance of HREPA and the lack of ensemble precipitation analysis (PA) at such spatiotemporal resolution in the literature motivate further investigations on transitional seasons and more advanced perturbation approaches.


2010 ◽  
Vol 27 (3) ◽  
pp. 528-546 ◽  
Author(s):  
Robert W. Helber ◽  
Jay F. Shriver ◽  
Charlie N. Barron ◽  
Ole Martin Smedstad

Abstract The impact of the number of satellite altimeters providing sea surface height anomaly (SSHA) information for a data assimilation system is evaluated using two comparison frameworks and two statistical methodologies. The Naval Research Laboratory (NRL) Layered Ocean Model (NLOM) dynamically interpolates satellite SSHA track data measured from space to produce high-resolution (eddy resolving) fields. The Modular Ocean Data Assimilation System (MODAS) uses the NLOM SSHA to produce synthetic three-dimensional fields of temperature and salinity over the global ocean. A series of case studies is defined where NLOM assimilates different combinations of data streams from zero to three altimeters. The resulting NLOM SSHA fields and the MODAS synthetic profiles are evaluated relative to independently observed ocean temperature and salinity profiles for the years 2001–03. The NLOM SSHA values are compared with the difference of the observed dynamic height from the climatological dynamic height. The synthetics are compared with observations using a measure of thermocline depth. Comparisons are done point for point and for 1° radius regions that are linearly fit over 2-month periods. To evaluate the impact of data outliers, statistical evaluations are done with traditional Gaussian statistics and also with robust nonparametric statistics. Significant error reduction is obtained, particularly in high SSHA variability regions, by including at least one altimeter. Given the limitation of these methods, the overall differences between one and three altimeters are significant only in bias. Data outliers increase Gaussian statistical error and error uncertainty compared to the same computations using nonparametric statistical methods.


2005 ◽  
Vol 133 (4) ◽  
pp. 829-843 ◽  
Author(s):  
Milija Zupanski ◽  
Dusanka Zupanski ◽  
Tomislava Vukicevic ◽  
Kenneth Eis ◽  
Thomas Vonder Haar

A new four-dimensional variational data assimilation (4DVAR) system is developed at the Cooperative Institute for Research in the Atmosphere (CIRA)/Colorado State University (CSU). The system is also called the Regional Atmospheric Modeling Data Assimilation System (RAMDAS). In its present form, the 4DVAR system is employing the CSU/Regional Atmospheric Modeling System (RAMS) nonhydrostatic primitive equation model. The Weather Research and Forecasting (WRF) observation operator is used to access the observations, adopted from the WRF three-dimensional variational data assimilation (3DVAR) algorithm. In addition to the initial conditions adjustment, the RAMDAS includes the adjustment of model error (bias) and lateral boundary conditions through an augmented control variable definition. Also, the control variable is defined in terms of the velocity potential and streamfunction instead of the horizontal winds. The RAMDAS is developed after the National Centers for Environmental Prediction (NCEP) Eta 4DVAR system, however with added improvements addressing its use in a research environment. Preliminary results with RAMDAS are presented, focusing on the minimization performance and the impact of vertical correlations in error covariance modeling. A three-dimensional formulation of the background error correlation is introduced and evaluated. The Hessian preconditioning is revisited, and an alternate algebraic formulation is presented. The results indicate a robust minimization performance.


2015 ◽  
Vol 2015 ◽  
pp. 1-17 ◽  
Author(s):  
Govindan Kutty ◽  
Xuguang Wang

The impact of observations can be dependent on many factors in a data assimilation (DA) system including data quality control, preprocessing, skill of the model, and the DA algorithm. The present study focuses on comparing the impacts of observations assimilated by two different DA algorithms. A three-dimensional ensemble-variational (3DEnsVar) hybrid data assimilation system was recently developed based on the Gridpoint Statistical Interpolation (GSI) data assimilation system and was implemented operationally for the National Center for Environmental Prediction (NCEP) Global Forecast System (GFS). One question to address is, how the impacts of observations on GFS forecasts differ when assimilated by the traditional GSI-three dimensional variational (3DVar) and the new 3DEnsVar. Experiments were conducted over a 6-week period during Northern Hemisphere winter season at a reduced resolution. For both the control and data denial experiments, the forecasts produced by 3DEnsVar were more accurate than GSI3DVar experiments. The results suggested that the observations were better and more effectively exploited to increment the background forecast in 3DEnsVar. On the other hand, in GSI3DVar, where the observation will be making mostly local, isotropic increments without proper flow dependent extrapolation is more sensitive to the number and types observations assimilated.


Sign in / Sign up

Export Citation Format

Share Document