Results 1  10
of
94
Ensemble Data Assimilation without Perturbed Observations
 MON. WEA. REV
, 2002
"... The ensemble Kalman filter (EnKF) is a data assimilation scheme based on the traditional Kalman filter update equation. An ensemble of forecasts are used to estimate the backgrounderror covariances needed to compute the Kalman gain. It is known that if the same observations and the same gain are ..."
Abstract

Cited by 278 (21 self)
 Add to MetaCart
The ensemble Kalman filter (EnKF) is a data assimilation scheme based on the traditional Kalman filter update equation. An ensemble of forecasts are used to estimate the backgrounderror covariances needed to compute the Kalman gain. It is known that if the same observations and the same gain are used to update each member of the ensemble, the ensemble will systematically underestimate analysiserror covariances. This will cause a degradation of subsequent analyses and may lead to filter divergence. For large ensembles, it is known that this problem can be alleviated by treating the observations as random variables, adding random perturbations to them with the correct statistics. Two important
Ensemble Square Root Filters
, 2003
"... Ensemble data assimilation methods assimilate observations using statespace estimation methods and lowrank representations of forecast and analysis error covariances. A key element of such methods is the transformation of the forecast ensemble into an analysis ensemble with appropriate statistics ..."
Abstract

Cited by 116 (7 self)
 Add to MetaCart
Ensemble data assimilation methods assimilate observations using statespace estimation methods and lowrank representations of forecast and analysis error covariances. A key element of such methods is the transformation of the forecast ensemble into an analysis ensemble with appropriate statistics. This transformation may be performed stochastically by treating observations as random variables, or deterministically by requiring that the updated analysis perturbations satisfy the Kalman filter analysis error covariance equation. Deterministic analysis ensemble updates are implementations of Kalman square root filters. The nonuniqueness of the deterministic transformation used in square root Kalman filters provides a framework to compare three recently proposed ensemble data assimilation methods.
OBSTACLES TO HIGHDIMENSIONAL PARTICLE FILTERING
"... Particle filters are ensemblebased assimilation schemes that, unlike the ensemble Kalman filter, employ a fully nonlinear and nonGaussian analysis step to compute the probability distribution function (pdf) of a system’s state conditioned on a set of observations. Evidence is provided that the ens ..."
Abstract

Cited by 94 (4 self)
 Add to MetaCart
Particle filters are ensemblebased assimilation schemes that, unlike the ensemble Kalman filter, employ a fully nonlinear and nonGaussian analysis step to compute the probability distribution function (pdf) of a system’s state conditioned on a set of observations. Evidence is provided that the ensemble size required for a successful particle filter scales exponentially with the problem size. For the simple example in which each component of the state vector is independent, Gaussian and of unit variance, and the observations are of each state component separately with independent, Gaussian errors, simulations indicate that the required ensemble size scales exponentially with the state dimension. In this example, the particle filter requires at least 1011 members when applied to a 200dimensional state. Asymptotic results, following the work of Bengtsson, Bickel and collaborators, are provided for two cases: one in which each prior state component is independent and identically distributed, and one in which both the prior pdf and the observation errors are Gaussian. The asymptotic theory reveals that, in both cases, the required ensemble size scales exponentially with the variance of the observation loglikelihood, rather than with the state dimension per se. 2
Sampling strategies and square root analysis schemes for the EnKF
"... this paper is to examine how different sampling strategies and implementations of the analysis scheme influence the quality of the results in the EnKF. It is shown that by selecting the initial ensemble, the model noise and the measurement perturbations wisely, it is possible to achieve a signific ..."
Abstract

Cited by 89 (1 self)
 Add to MetaCart
this paper is to examine how different sampling strategies and implementations of the analysis scheme influence the quality of the results in the EnKF. It is shown that by selecting the initial ensemble, the model noise and the measurement perturbations wisely, it is possible to achieve a significant improvement in the EnKF results, using the same number of members in the ensemble
A Local Least Squares Framework for Ensemble Filtering
, 2003
"... Many methods using ensemble integrations of prediction models as integral parts of data assimilation have appeared in the atmospheric and oceanic literature. In general, these methods have been derived from the Kalman filter and have been known as ensemble Kalman filters. A more general class of m ..."
Abstract

Cited by 88 (9 self)
 Add to MetaCart
Many methods using ensemble integrations of prediction models as integral parts of data assimilation have appeared in the atmospheric and oceanic literature. In general, these methods have been derived from the Kalman filter and have been known as ensemble Kalman filters. A more general class of methods including these ensemble Kalman filter methods is derived starting from the nonlinear filtering problem. When working in a joint state observation space, many features of ensemble filtering algorithms are easier to derive and compare. The ensemble filter methods derived here make a (local) least squares assumption about the relation between prior distributions of an observation variable and model state variables. In this context, the update procedure applied when a new observation becomes available can be described in two parts. First, an update increment is computed for each prior ensemble estimate of the observation variable by applying a scalar ensemble filter. Second, a linear regression of the prior ensemble sample of each state variable on the observation variable is performed to compute update increments for each state variable ensemble member from corresponding observation variable increments. The regression can be applied globally or locally using Gaussian kernel methods.
Uncertainty assessment of hydrologic model states and parameters: Sequential data assimilation using the particle filter
 Water Resour. Res. 2005
"... Let us know how access to this document benefits you. ..."
Abstract

Cited by 46 (5 self)
 Add to MetaCart
(Show Context)
Let us know how access to this document benefits you.
R.: Extended versus ensemble Kalman filtering for land data assimilation
 J. Hydrometeor
"... The performance of the extended Kalman filter (EKF) and the ensemble Kalman filter (EnKF) are assessed for soil moisture estimation. In a twin experiment for the southeastern United States synthetic observations of nearsurface soil moisture are assimilated once every 3 days, neglecting horizontal e ..."
Abstract

Cited by 46 (0 self)
 Add to MetaCart
The performance of the extended Kalman filter (EKF) and the ensemble Kalman filter (EnKF) are assessed for soil moisture estimation. In a twin experiment for the southeastern United States synthetic observations of nearsurface soil moisture are assimilated once every 3 days, neglecting horizontal error correlations and treating catchments independently. Both filters provide satisfactory estimates of soil moisture. The average actual estimation error in volumetric moisture content of the soil profile is 2.2 % for the EKF and 2.2 % (or 2.1%; or 2.0%) for the EnKF with 4 (or 10; or 500) ensemble members. Expected error covariances of both filters generally differ from actual estimation errors. Nevertheless, nonlinearities in soil processes are treated adequately by both filters. In the application presented herein the EKF and the EnKF with four ensemble members are equally accurate at comparable computational cost. Because of its flexibility and its performance in this study, the EnKF is a promising approach for soil moisture initialization problems. 1.
A VarianceMinimizing Filter for LargeScale Applications
, 2003
"... A truly varianceminimizing filter is introduced and its performance is demonstrated with the KortewegDeVries (KdV) equation and with a multilayer quasigeostrophic model of the ocean area around South Africa. It is recalled ..."
Abstract

Cited by 42 (1 self)
 Add to MetaCart
A truly varianceminimizing filter is introduced and its performance is demonstrated with the KortewegDeVries (KdV) equation and with a multilayer quasigeostrophic model of the ocean area around South Africa. It is recalled
Ensemblebased atmospheric data assimilation
, 2004
"... Ensemblebased data assimilation techniques are being explored as possible alternatives to current operational analysis techniques such as 3 or 4dimensional variational assimilation. Ensemblebased assimilation techniques utilize an ensemble of parallel data assimilation and forecast cycles. The ..."
Abstract

Cited by 42 (2 self)
 Add to MetaCart
Ensemblebased data assimilation techniques are being explored as possible alternatives to current operational analysis techniques such as 3 or 4dimensional variational assimilation. Ensemblebased assimilation techniques utilize an ensemble of parallel data assimilation and forecast cycles. The backgrounderror covariances are estimated using the forecast ensemble and are used to produce an ensemble of analyses. The backgrounderror covariances are flow dependent and often have very complicated structure, providing a very different adjustment to the observations than are seen from methods such as 3 dimensional variational assimilation. Though computationally expensive, ensemblebased techniques are relatively easy to code, since no adjoint nor tangentlinear models are required, and previous tests in simple models suggest that dramatic improvements over existing operational methods may be possible. A review of the ensemblebased assimilation is provided here, starting from the basic concepts of Bayesian assimilation. Without some simplification, full Bayesian assimilation is computationally impossible for model states of large dimension. Assuming normality of error statistics and linearity of error growth, the state and its error covariance may be predicted optimally using Kalman filter (KF) techniques. The ensemble Kalman filter (EnKF) is then described. The EnKF is an approximation to the KF in that backgrounderror covariances are estimated from a finite ensemble of forecasts. However, no assumptions about linearity of error growth are made. Recent algorithmic variants on the standard EnKF are also described, as well as methods for simplifying the computations and increasing the accuracy. Examples of ensemblebased assimilations are provided in simple and more realistic dynamical systems.
T.: A method for assimilation of Lagrangian data
 Weather Rev
"... Difficulties in the assimilation of Lagrangian data arise because the state of the prognostic model is generally described in terms of Eulerian variables computed on a fixed grid in space, as a result there is no direct connection between the model variables and Lagrangian observations that carry ti ..."
Abstract

Cited by 28 (2 self)
 Add to MetaCart
Difficulties in the assimilation of Lagrangian data arise because the state of the prognostic model is generally described in terms of Eulerian variables computed on a fixed grid in space, as a result there is no direct connection between the model variables and Lagrangian observations that carry timeintegrated information. A method is presented for assimilating Lagrangian tracer positions, observed at discrete times, directly into the model. The idea is to augment the model with tracer advection equations and to track the correlations between the flow and the tracers via the extended Kalman filter. The augmented model state vector includes tracer coordinates and is updated through the correlations to the observed tracers. The technique is tested for point vortex flows: an NF point vortex system with a Gaussian noise term is modeled by its deterministic counterpart. Positions of ND tracer particles are observed at regular time intervals and assimilated into the model. Numerical experiments demonstrate successful system tracking for (NF, ND) � (2, 1), (4, 2), provided the observations are reasonably frequent and accurate and the system noise level is not too high. The performance of the filter strongly depends on initial tracer positions (drifter launch locations). Analysis of this dependence shows that the good launch locations are separated from the bad ones by Lagrangian flow structures (separatrices or invariant manifolds of the velocity field). The method is compared to an alternative indirect approach, where the flow velocity, estimated from two (or more) consecutive drifter observations, is assimilated directly into the model. 1.