scholarly journals Efficient Local Error Parameterizations for Square Root or Ensemble Kalman Filters: Application to a Basin-Scale Ocean Turbulent Flow

2011 ◽  
Vol 139 (2) ◽  
pp. 474-493 ◽  
Author(s):  
Jean-Michel Brankart ◽  
Emmanuel Cosme ◽  
Charles-Emmanuel Testut ◽  
Pierre Brasseur ◽  
Jacques Verron
Keyword(s):  
Turbulent Flow ◽  
Covariance Matrix ◽  
Kalman Filters ◽  
Square Root ◽  
Ensemble Kalman Filters ◽  
Error Covariance Matrix ◽  
Square Roots ◽  
Schur Product ◽  
Error Covariance ◽  
Basin Scale

Abstract In large-sized atmospheric or oceanic applications of square root or ensemble Kalman filters, it is often necessary to introduce the prior assumption that long-range correlations are negligible and force them to zero using a local parameterization, supplementing the ensemble or reduced-rank representation of the covariance. One classic algorithm to perform this operation consists of taking the Schur product of the ensemble covariance matrix by a local support correlation matrix. However, with this parameterization, the square root of the forecast error covariance matrix is no more directly available, so that any observational update algorithm requiring this square root must include an additional step to compute local square roots from the Schur product. This computation generates an additional numerical cost or produces high-rank square roots, which may deprive the observational update from its original efficiency. In this paper, it is shown how efficient local square root parameterizations can be obtained, for use with a specific square root formulation (i.e., eigenbasis algorithm) of the observational update. Comparisons with the classic algorithm are provided, mainly in terms of consistency, accuracy, and computational complexity. As an application, the resulting parameterization is used to estimate maps of dynamic topography characterizing a basin-scale ocean turbulent flow. Even with this moderate-sized system (a 2200-km-wide square basin with 100-km-wide mesoscale eddies), it is observed that more than 1000 ensemble members are necessary to faithfully represent the global correlation patterns, and that a local parameterization is needed to produce correct covariances with moderate-sized ensembles. Comparisons with the exact solution show that the use of local square roots is able to improve the accuracy of the updated ensemble mean and the consistency of the updated ensemble variance. With the eigenbasis algorithm, optimal adaptive estimates of scaling factors for the forecast and observation error covariance matrix can also be obtained locally at negligible additional numerical cost. Finally, a comparison of the overall computational cost illustrates the decisive advantage that efficient local square root parameterizations may have to deal simultaneously with a larger number of observations and avoid data thinning as much as possible.

scholarly journals Efficient Parameterization of the Observation Error Covariance Matrix for Square Root or Ensemble Kalman Filters: Application to Ocean Altimetry

2009 ◽  
Vol 137 (6) ◽  
pp. 1908-1927 ◽  
Author(s):  
Jean-Michel Brankart ◽  
Clément Ubelmann ◽  
Charles-Emmanuel Testut ◽  
Emmanuel Cosme ◽  
Pierre Brasseur ◽  
...  
Keyword(s):  
Covariance Matrix ◽  
Kalman Filters ◽  
Signal Amplitude ◽  
Observation Error ◽  
Square Root ◽  
Ensemble Kalman Filters ◽  
Error Covariance Matrix ◽  
Standard Algorithm ◽  
Error Covariance ◽  
North Brazil

Abstract In the Kalman filter standard algorithm, the computational complexity of the observational update is proportional to the cube of the number y of observations (leading behavior for large y). In realistic atmospheric or oceanic applications, involving an increasing quantity of available observations, this often leads to a prohibitive cost and to the necessity of simplifying the problem by aggregating or dropping observations. If the filter error covariance matrices are in square root form, as in square root or ensemble Kalman filters, the standard algorithm can be transformed to be linear in y, providing that the observation error covariance matrix is diagonal. This is a significant drawback of this transformed algorithm and often leads to an assumption of uncorrelated observation errors for the sake of numerical efficiency. In this paper, it is shown that the linearity of the transformed algorithm in y can be preserved for other forms of the observation error covariance matrix. In particular, quite general correlation structures (with analytic asymptotic expressions) can be simulated simply by augmenting the observation vector with differences of the original observations, such as their discrete gradients. Errors in ocean altimetric observations are spatially correlated, as for instance orbit or atmospheric errors along the satellite track. Adequately parameterizing these correlations can directly improve the quality of observational updates and the accuracy of the associated error estimates. In this paper, the example of the North Brazil Current circulation is used to demonstrate the importance of this effect, which is especially significant in that region of moderate ratio between signal amplitude and observation noise, and to show that the efficient parameterization that is proposed for the observation error correlations is appropriate to take it into account. Adding explicit gradient observations also receives a physical justification. This parameterization is thus proved to be useful to ocean data assimilation systems that are based on square root or ensemble Kalman filters, as soon as the number of observations becomes penalizing, and if a sophisticated parameterization of the observation error correlations is required.

scholarly journals Efficient Adaptive Error Parameterizations for Square Root or Ensemble Kalman Filters: Application to the Control of Ocean Mesoscale Signals

2010 ◽  
Vol 138 (3) ◽  
pp. 932-950 ◽  
Author(s):  
Jean-Michel Brankart ◽  
Emmanuel Cosme ◽  
Charles-Emmanuel Testut ◽  
Pierre Brasseur ◽  
Jacques Verron
Keyword(s):  
Kalman Filter ◽  
Covariance Matrix ◽  
Error Estimates ◽  
Forecast Error ◽  
Observation Error ◽  
Square Root ◽  
Error Statistics ◽  
Error Covariance Matrix ◽  
Error Covariance ◽  
Optimal Estimates

Abstract In Kalman filter applications, an adaptive parameterization of the error statistics is often necessary to avoid filter divergence, and prevent error estimates from becoming grossly inconsistent with the real error. With the classic formulation of the Kalman filter observational update, optimal estimates of general adaptive parameters can only be obtained at a numerical cost that is several times larger than the cost of the state observational update. In this paper, it is shown that there exists a few types of important parameters for which optimal estimates can be computed at a negligible numerical cost, as soon as the computation is performed using a transformed algorithm that works in the reduced control space defined by the square root or ensemble representation of the forecast error covariance matrix. The set of parameters that can be efficiently controlled includes scaling factors for the forecast error covariance matrix, scaling factors for the observation error covariance matrix, or even a scaling factor for the observation error correlation length scale. As an application, the resulting adaptive filter is used to estimate the time evolution of ocean mesoscale signals using observations of the ocean dynamic topography. To check the behavior of the adaptive mechanism, this is done in the context of idealized experiments, in which model error and observation error statistics are known. This ideal framework is particularly appropriate to explore the ill-conditioned situations (inadequate prior assumptions or uncontrollability of the parameters) in which adaptivity can be misleading. Overall, the experiments show that, if used correctly, the efficient optimal adaptive algorithm proposed in this paper introduces useful supplementary degrees of freedom in the estimation problem, and that the direct control of these statistical parameters by the observations increases the robustness of the error estimates and thus the optimality of the resulting Kalman filter.

Enhancing the impact of IASI observations through an updated observation-error covariance matrix

2016 ◽  
Vol 142 (697) ◽  
pp. 1767-1780 ◽  
Author(s):  
Niels Bormann ◽  
Massimo Bonavita ◽  
Rossana Dragani ◽  
Reima Eresmaa ◽  
Marco Matricardi ◽  
...  
Keyword(s):  
Covariance Matrix ◽  
Observation Error ◽  
Error Covariance Matrix ◽  
Error Covariance ◽  
The Impact

Impact of AMSU-A Radiances in a Column Ensemble Kalman Filter

2018 ◽  
Vol 146 (12) ◽  
pp. 3949-3976 ◽  
Author(s):  
Herschel L. Mitchell ◽  
P. L. Houtekamer ◽  
Sylvain Heilliette
Keyword(s):  
Covariance Matrix ◽  
Background Error ◽  
Error Covariance Matrix ◽  
Background Error Covariance ◽  
Error Covariance ◽  
Microwave Sounding Unit ◽  
The Matrix ◽  
Microwave Sounding ◽  
Background Temperature ◽  
The Difference

Abstract A column EnKF, based on the Canadian global EnKF and using the RTTOV radiative transfer (RT) model, is employed to investigate issues relating to the EnKF assimilation of Advanced Microwave Sounding Unit-A (AMSU-A) radiance measurements. Experiments are performed with large and small ensembles, with and without localization. Three different descriptions of background temperature error are considered: 1) using analytical vertical modes and hypothetical spectra, 2) using the vertical modes and spectrum of a covariance matrix obtained from the global EnKF after 2 weeks of cycling, and 3) using the vertical modes and spectrum of the static background error covariance matrix employed to initiate a global data assimilation cycle. It is found that the EnKF performs well in some of the experiments with background error description 1, and yields modest error reductions with background error description 3. However, the EnKF is virtually unable to reduce the background error (even when using a large ensemble) with background error description 2. To analyze these results, the different background error descriptions are viewed through the prism of the RT model by comparing the trace of the matrix , where is the RT model and is the background error covariance matrix. Indeed, this comparison is found to explain the difference in the results obtained, which relates to the degree to which deep modes are, or are not, present in the different background error covariances. The results suggest that, after 2 weeks of cycling, the global EnKF has virtually eliminated all background error structures that can be “seen” by the AMSU-A radiances.

RECURSIVE UPDATING OF ERROR COVARIANCE MATRIX IN SUBSPACE METHODS

2006 ◽  
Vol 39 (1) ◽  
pp. 285-290 ◽  
Author(s):  
Yoshinori Takei ◽  
Hidehito Nanto ◽  
Shunshoku Kanae ◽  
Zi-Jiang Yang ◽  
Kiyoshi Wada
Keyword(s):  
Covariance Matrix ◽  
Subspace Methods ◽  
Error Covariance Matrix ◽  
Error Covariance ◽  
Recursive Updating

Intelligent Error Covariance Matrix Resetting for Maneuver Target Tracking

2008 ◽  
Vol 8 (12) ◽  
pp. 2279-2285 ◽  
Author(s):  
M.H. Bahari ◽  
A. Karsaz ◽  
M.B. Naghibi-S
Keyword(s):  
Target Tracking ◽  
Covariance Matrix ◽  
Error Covariance Matrix ◽  
Error Covariance

scholarly journals Best linear unbiased estimation for varying probability with and without replacement sampling

2019 ◽  
Vol 7 (1) ◽  
pp. 78-91
Author(s):  
Stephen Haslett
Keyword(s):  
Parameter Estimation ◽  
Covariance Matrix ◽  
Linear Models ◽  
Unbiased Estimation ◽  
Error Covariance Matrix ◽  
Best Linear Unbiased Estimation ◽  
Error Covariance ◽  
Inclusion Probabilities ◽  
Best Linear Unbiased ◽  
Linear Unbiased Estimation

Abstract When sample survey data with complex design (stratification, clustering, unequal selection or inclusion probabilities, and weighting) are used for linear models, estimation of model parameters and their covariance matrices becomes complicated. Standard fitting techniques for sample surveys either model conditional on survey design variables, or use only design weights based on inclusion probabilities essentially assuming zero error covariance between all pairs of population elements. Design properties that link two units are not used. However, if population error structure is correlated, an unbiased estimate of the linear model error covariance matrix for the sample is needed for efficient parameter estimation. By making simultaneous use of sampling structure and design-unbiased estimates of the population error covariance matrix, the paper develops best linear unbiased estimation (BLUE) type extensions to standard design-based and joint design and model based estimation methods for linear models. The analysis covers both with and without replacement sample designs. It recognises that estimation for with replacement designs requires generalized inverses when any unit is selected more than once. This and the use of Hadamard products to link sampling and population error covariance matrix properties are central topics of the paper. Model-based linear model parameter estimation is also discussed.

scholarly journals Estimating model error covariance matrix parameters in extended Kalman filtering

2014 ◽  
Vol 21 (5) ◽  
pp. 919-927 ◽  
Author(s):  
A. Solonen ◽  
J. Hakkarainen ◽  
A. Ilin ◽  
M. Abbas ◽  
A. Bibov
Keyword(s):  
Covariance Matrix ◽  
Weather Prediction ◽  
Estimation Method ◽  
Model Error ◽  
Tuning Parameter ◽  
Extended Kalman Filtering ◽  
Error Covariance Matrix ◽  
Nonlinear Dynamical ◽  
Error Covariance

Abstract. The extended Kalman filter (EKF) is a popular state estimation method for nonlinear dynamical models. The model error covariance matrix is often seen as a tuning parameter in EKF, which is often simply postulated by the user. In this paper, we study the filter likelihood technique for estimating the parameters of the model error covariance matrix. The approach is based on computing the likelihood of the covariance matrix parameters using the filtering output. We show that (a) the importance of the model error covariance matrix calibration depends on the quality of the observations, and that (b) the estimation approach yields a well-tuned EKF in terms of the accuracy of the state estimates and model predictions. For our numerical experiments, we use the two-layer quasi-geostrophic model that is often used as a benchmark model for numerical weather prediction.

Sign in / Sign up

Export Citation Format

Share Document