Representing model error in ensemble data assimilation

Abstract. The paper investigates a method to represent model error in the ensemble data assimilation (EDA) system. The ECMWF operational EDA simulates the effect of both observations and model uncertainties. Observation errors are represented by perturbations with statistics characterized by the observation error covariance matrix whilst the model uncertainties are represented by stochastic perturbations added to the physical tendencies to simulate the effect of random errors in the physical parameterizations (ST-method). In this work an alternative method (XB-method) is proposed to simulate model uncertainties by adding perturbations to the model background field. In this way the error represented is not just restricted to model error in the usual sense but potentially extends to any form of background error. The perturbations have the same correlation as the background error covariance matrix and their magnitude is computed from comparing the high-resolution operational innovation variances with the ensemble variances when the ensemble is obtained by perturbing only the observations (OBS-method). The XB-method has been designed to represent the short range model error relevant for the data assimilation window. Spread diagnostic shows that the XB-method generates a larger spread than the ST-method that is operationally used at ECMWF, in particular in the extratropics. Three-dimensional normal-mode diagnostics indicate that XB-EDA spread projects more than the spread from the other EDAs onto the easterly inertia-gravity modes associated with equatorial Kelvin waves, tropical dynamics and, in general, model error sources. The background error statistics from the above described EDAs have been employed in the assimilation system. The assimilation system performance showed that the XB-method background error statistics increase the observation influence in the analysis process. The other EDA background error statistics, when inflated by a global factor, generate analyses with 30–50% smaller degree of freedom of signal. XB-EDA background error variances have not been inflated. The presented EDAs have been used to generate the initial perturbations of the ECMWF ensemble prediction system (EPS) of which the XB-EDA induces the largest EPS spread, also in the medium range, leading to a more reliable ensemble. Compared to ST-EDA, XB-EDA leads to a small improvement of the EPS ignorance skill score at day 3 and 7.

[1]  Jean-Jacques Morcrette,et al.  Linearized radiation and cloud schemes in the ECMWF model: Development and evaluation , 2002 .

[2]  Jean-Noël Thépaut,et al.  The Spatial Structure of Observation Errors in Atmospheric Motion Vectors from Geostationary Satellite Data , 2003 .

[3]  H.-L. Huang,et al.  Estimating effective data density in a satellite retrieval or an objective analysis , 1993 .

[4]  Gérald Desroziers,et al.  Diagnosis and adaptive tuning of observation‐error parameters in a variational assimilation , 2001 .

[5]  Renate Hagedorn,et al.  Ensemble prediction: a pedagogical perspective , 2008 .

[6]  Fuqing Zhang,et al.  Tests of an Ensemble Kalman Filter for Mesoscale and Regional-Scale Data Assimilation. Part I: Perfect Model Experiments , 2006 .

[7]  Adrian M. Tompkins,et al.  A cloud scheme for data assimilation: Description and initial tests , 2004 .

[8]  J. Whitaker,et al.  Accounting for the Error due to Unresolved Scales in Ensemble Data Assimilation: A Comparison of Different Approaches , 2005 .

[9]  G. Shutts.  A kinetic energy backscatter algorithm for use in ensemble prediction systems , 2005 .

[10]  David J. Stensrud,et al.  Surface Data Assimilation Using an Ensemble Kalman Filter Approach with Initial Condition and Model Physics Uncertainties , 2005 .

[11]  J. Derber,et al.  A reformulation of the background error covariance in the ECMWF global data assimilation system , 1999 .

[12]  T. Palmer,et al.  Stochastic representation of model uncertainties in the ECMWF ensemble prediction system , 2007 .

[13]  J. Tribbia,et al.  Balance properties of the short‐range forecast errors in the ECMWF 4D‐Var ensemble , 2013 .

[14]  Martin Leutbecher,et al.  A Spectral Stochastic Kinetic Energy Backscatter Scheme and Its Impact on Flow-Dependent Predictability in the ECMWF Ensemble Prediction System , 2009 .

[15]  Feng Gao,et al.  Adaptive Tuning of Numerical Weather Prediction Models: Randomized GCV in Three- and Four-Dimensional Data Assimilation , 1995 .

[16]  Roberto Buizza,et al.  The new VarEPS-monthly forecasting system: A first step towards seamless prediction , 2008 .

[17]  W. Briggs.  Statistical Methods in the Atmospheric Sciences , 2007 .

[18]  Riccardo Benedetti,et al.  Scoring Rules for Forecast Verification , 2010 .

[19]  Lars Isaksen,et al.  Potential use of an ensemble of analyses in the ECMWF Ensemble Prediction System , 2008 .

[20]  Fuqing Zhang,et al.  Tests of an Ensemble Kalman Filter for Mesoscale and Regional-Scale Data Assimilation. Part III: Comparison with 3DVAR in a Real-Data Case Study , 2008 .

[21]  mith,et al.  Evaluating Probabilistic Forecasts Using Information Theory , 2002 .

[22]  Jeffrey L. Anderson,et al.  Balance of the Background-Error Variances in the Ensemble Assimilation System DART/CAM , 2011 .

[23]  P. L. Houtekamer,et al.  Ensemble Kalman filtering , 2005 .

[24]  Frederic Vitart,et al.  An Ensemble Generation Method for Seasonal Forecasting with an Ocean-Atmosphere Coupled Model , 2005 .

[25]  D. S. Sivia,et al.  Data Analysis , 1996, Encyclopedia of Evolutionary Psychological Science.

[26]  Paul Poli,et al.  Diagnosis of observation, background and analysis‐error statistics in observation space , 2005 .

[27]  L. Raynaud,et al.  Spatial averaging of ensemble‐based background‐error variances , 2008 .

[28]  R. Buizza,et al.  A Comparison of the ECMWF, MSC, and NCEP Global Ensemble Prediction Systems , 2005 .

[29]  Daniel S. Wilks,et al.  Statistical Methods in the Atmospheric Sciences: An Introduction , 1995 .

[30]  Carla Cardinali,et al.  Observation Influence Diagnostic of a Data Assimilation System , 2013 .

[31]  Mats Hamrud,et al.  The new ECMWF VAREPS (Variable Resolution Ensemble Prediction System) , 2007 .

[32]  O. Talagrand,et al.  A Posteriori Validation of Assimilation Algorithms , 2003 .

[33]  L. Raynaud,et al.  Objective filtering of ensemble‐based background‐error variances , 2009 .

[34]  A. Simmons,et al.  The ECMWF operational implementation of four‐dimensional variational assimilation. I: Experimental results with simplified physics , 2007 .

[35]  Andrew C. Lorenc,et al.  Analysis methods for numerical weather prediction , 1986 .

[36]  F. Molteni,et al.  The ECMWF Ensemble Prediction System: Methodology and validation , 1996 .

[37]  J. Tukey.  Data analysis, computation and mathematics , 1972 .

[38]  M. Buehner,et al.  Intercomparison of Variational Data Assimilation and the Ensemble Kalman Filter for Global Deterministic NWP. Part II: One-Month Experiments with Real Observations , 2010 .

[39]  Gérald Desroziers,et al.  An extended specification of flow‐dependent background error variances in the Météo‐France global 4D‐Var system , 2011 .

[40]  JAN HASELER,et al.  The new Ensemble of Data Assimilations , 2011 .

[41]  G. Brier.  VERIFICATION OF FORECASTS EXPRESSED IN TERMS OF PROBABILITY , 1950 .

[42]  Philippe Lopez,et al.  A convection scheme for data assimilation: Description and initial tests , 2005 .

[43]  T. Palmer,et al.  Stochastic parametrization and model uncertainty , 2009 .

[44]  X. Deng,et al.  Model Error Representation in an Operational Ensemble Kalman Filter , 2009 .

[45]  Erik Andersson,et al.  Influence‐matrix diagnostic of a data assimilation system , 2004 .

[46]  M. Buehner,et al.  Atmospheric Data Assimilation with an Ensemble Kalman Filter: Results with Real Observations , 2005 .

[47]  M. Fisher,et al.  Background Error Covariance Modelling , 2003 .

[48]  Gérald Desroziers,et al.  Accounting for model error in the Météo‐France ensemble data assimilation system , 2012 .

[49]  P. L. Houtekamer,et al.  A System Simulation Approach to Ensemble Prediction , 1996 .

[50]  Roy E. Welsch,et al.  Efficient Computing of Regression Diagnostics , 1981 .