Sigma-Point Kalman Filter Data Assimilation Methods for Strongly Nonlinear Systems

2009 ◽  
Vol 66 (2) ◽  
pp. 261-285 ◽  
Author(s):  
Jaison Thomas Ambadan ◽  
Youmin Tang

Abstract Performance of an advanced, derivativeless, sigma-point Kalman filter (SPKF) data assimilation scheme in a strongly nonlinear dynamical model is investigated. The SPKF data assimilation scheme is compared against standard Kalman filters such as the extended Kalman filter (EKF) and ensemble Kalman filter (EnKF) schemes. Three particular cases—namely, the state, parameter, and joint estimation of states and parameters from a set of discontinuous noisy observations—are studied. The problems associated with the use of tangent linear model (TLM) or Jacobian when using standard Kalman filters are eliminated when using SPKF data assimilation algorithms. Further, the constraints and issues of SPKF data assimilation in real ocean or atmospheric models are emphasized. A reduced sigma-point subspace model is proposed and investigated for higher-dimensional systems. A low-dimensional Lorenz 1963 model and a higher-dimensional Lorenz 1995 model are used as the test beds for data assimilation experiments. The results of SPKF data assimilation schemes are compared with those of standard EKF and EnKF, in which a highly nonlinear chaotic case is studied. It is shown that the SPKF is capable of estimating the model state and parameters with better accuracy than EKF and EnKF. Numerical experiments showed that in all cases the SPKF can give consistent results with better assimilation skills than EnKF and EKF and can overcome the drawbacks associated with the use of EKF and EnKF.

2009 ◽  
Vol 66 (11) ◽  
pp. 3498-3500 ◽  
Author(s):  
Thomas M. Hamill ◽  
Jeffrey S. Whitaker ◽  
Jeffrey L. Anderson ◽  
Chris Snyder

2007 ◽  
Vol 64 (4) ◽  
pp. 1116-1140 ◽  
Author(s):  
David Kuhl ◽  
Istvan Szunyogh ◽  
Eric J. Kostelich ◽  
Gyorgyi Gyarmati ◽  
D. J. Patil ◽  
...  

Abstract In this paper, the spatiotemporally changing nature of predictability is studied in a reduced-resolution version of the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS), a state-of-the-art numerical weather prediction model. Atmospheric predictability is assessed in the perfect model scenario for which forecast uncertainties are entirely due to uncertainties in the estimates of the initial states. Uncertain initial conditions (analyses) are obtained by assimilating simulated noisy vertical soundings of the “true” atmospheric states with the local ensemble Kalman filter (LEKF) data assimilation scheme. This data assimilation scheme provides an ensemble of initial conditions. The ensemble mean defines the initial condition of 5-day deterministic model forecasts, while the time-evolved members of the ensemble provide an estimate of the evolving forecast uncertainties. The observations are randomly distributed in space to ensure that the geographical distribution of the analysis and forecast errors reflect predictability limits due to the model dynamics and are not affected by inhomogeneities of the observational coverage. Analysis and forecast error statistics are calculated for the deterministic forecasts. It is found that short-term forecast errors tend to grow exponentially in the extratropics and linearly in the Tropics. The behavior of the ensemble is explained by using the ensemble dimension (E dimension), a spatiotemporally evolving measure of the evenness of the distribution of the variance between the principal components of the ensemble-based forecast error covariance matrix. It is shown that in the extratropics the largest forecast errors occur for the smallest E dimensions. Since a low value of the E dimension guarantees that the ensemble can capture a large portion of the forecast error, the larger the forecast error the more certain that the ensemble can fully capture the forecast error. In particular, in regions of low E dimension, ensemble averaging is an efficient error filter and the ensemble spread provides an accurate prediction of the upper bound of the error in the ensemble-mean forecast.


2008 ◽  
Vol 65 (6) ◽  
pp. 1991-2001 ◽  
Author(s):  
Catherine Heyraud ◽  
Wanda Szyrmer ◽  
Stéphane Laroche ◽  
Isztar Zawadzki

Abstract In this paper a simplified UHF-band backscattering parameterization for individual melting snowflakes is proposed. This parameterization is a function of the density, shape, and melted fraction, and is used here in a brightband bulk modeling study. A 1D bulk model is developed where aggregation and breakup are neglected. Model results are in good agreement with detailed bin-model results and simulate the radar brightband observations well. It is shown the model can be seen as an observation operator that could be introduced into a data assimilation scheme to extract information contained in the radar data measurements.


Sign in / Sign up

Export Citation Format

Share Document