Accelerating the spin‐up of Ensemble Kalman Filtering

Ensemble Kalman Filter (EnKF) may have a longer spin-up time to reach its asymptotic level of accuracy than the corresponding spin-up time in variational methods (3D-Var or 4D-Var). During the spin-up EnKF has to fulfill two independent requirements, namely that the ensemble mean be close to the true state, and that the ensemble perturbations represent the ‘errors of the day’. As a result, there are cases, such as radar observations of a severe storm, or regional forecast of a hurricane, where EnKF may spin-up too slowly to be useful. A heuristic scheme is proposed to accelerate the spin-up of EnKF by applying a no-cost Ensemble Kalman Smoother, and using the observations more than once in each assimilation window during spin-up in order to maximize the initial extraction of information. The performance of this scheme is tested with the Local Ensemble Transform Kalman Filter (LETKF) implemented in a quasi-geostrophic model, which requires a very long spin-up time when initialized from random initial perturbations from a uniform distribution. Results show that with the new ‘running in place’ (RIP) scheme the LETKF spins up and converges to the optimal level of error faster than 3D-Var or 4D-Var, even in the absence of any prior information. Additional computations (2 to 12 iterations for each assimilation window) are only required during the initial spin-up, since the scheme naturally returns to the original LETKF after spin-up is achieved. RIP also accelerates spin-up when the initial perturbations are drawn from a well-tuned 3D-Var background-error covariance, rather than being uniform noise, and fewer iterations and RIP cycles are required than in the case without such prior information. Copyright © 2010 Royal Meteorological Society

[1]  Joaquim Ballabrera-Poy,et al.  Response to the discussion on “4-D-Var or EnKF?” by Nils Gustafsson , 2007 .

[2]  Nils Gustafsson,et al.  Discussion on ‘4D-Var or EnKF?’ , 2007 .

[3]  J. Yorke,et al.  Four-dimensional ensemble Kalman filtering , 2004 .

[4]  J. Whitaker,et al.  Ensemble Square Root Filters , 2003, Statistical Methods for Climate Scientists.

[5]  Istvan Szunyogh,et al.  Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter , 2005, physics/0511236.

[6]  J. Whitaker,et al.  Ensemble Data Assimilation without Perturbed Observations , 2002 .

[7]  Istvan Szunyogh,et al.  A Local Ensemble Kalman Filter for Atmospheric Data Assimilation , 2002 .

[8]  Chris Snyder,et al.  A Comparison between the 4DVAR and the Ensemble Kalman Filter Techniques for Radar Data Assimilation , 2005 .

[9]  E. Kalnay,et al.  Comparison of Local Ensemble Transform Kalman Filter, 3DVAR, and 4DVAR in a Quasigeostrophic Model , 2009 .

[10]  Jian-Wen Bao,et al.  A Case Study of Cyclogenesis Using a Model Hierarchy , 1996 .

[11]  John Derber,et al.  The National Meteorological Center's spectral-statistical interpolation analysis system , 1992 .

[12]  Andrew C. Lorenc,et al.  The potential of the ensemble Kalman filter for NWP—a comparison with 4D‐Var , 2003 .

[13]  R Hakkert,et al.  Buenos Aires , 2011 .

[14]  E. Kalnay,et al.  Four-dimensional ensemble Kalman filtering , 2004 .

[15]  Eugenia Kalnay,et al.  Weight interpolation for efficient data assimilation with the Local Ensemble Transform Kalman Filter , 2009 .

[16]  Junjie Liu,et al.  Applications of the LETKF to adaptive observations, analysis sensitivity, observation impact and the assimilation of moisture , 2007 .

[17]  Todd D. Ringler,et al.  Initiation of ensemble data assimilation , 2006 .

[18]  Craig H. Bishop,et al.  Adaptive sampling with the ensemble transform Kalman filter , 2001 .

[19]  E. Kalnay,et al.  C ○ 2007 The Authors , 2006 .

[20]  S. Ravela,et al.  Fast ensemble smoothing , 2007 .