Chemical state estimation for the middle atmosphere by four‐dimensional variational data assimilation: A posteriori validation of error statistics in observation space

[1] Chemical state analyses of the atmosphere based on data assimilation may be degraded by inconsistent covariances of background and observation errors. An efficient method to calculate consistency diagnostics for background and observation errors in observation space is applied to analyses of the four-dimensional variational stratospheric chemistry data assimilation system SACADA (Synoptic Analysis of Chemical Constituents by Advanced Data Assimilation). A background error covariance model for the assimilation of Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) ozone retrievals is set up and optimized. It is shown that a significant improvement of the assimilation system performance is attained through the use of this covariance model compared to a simple covariance formulation, which assumes background errors to be a fixed fraction of the field value. The forecast skill, measured by the distance between the model forecast and MIPAS observations, is shown to improve. Further, an evaluation of analyses with independent data from the Halogen Observation Experiment (HALOE), the Stratospheric Aerosol and Gas Experiment II (SAGE II), and ozone sondes reveals that the standard deviation of ozone analyses with respect to these instruments is reduced throughout the middle stratosphere. Compared to the impact of background error variances on analysis quality, it is found that the precise specification of spatial background error correlations appears to be less critical if observations are spatially and temporally dense. Results indicate that ozone forecast errors of a state of the art stratospheric chemistry assimilation system are of the same order of magnitude as MIPAS observation errors.

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

[2]  W. Lahoz,et al.  Dynamical evolution of the 2003 southern hemisphere stratospheric winter using Envisat trace‐gas observations , 2006 .

[3]  John C. Gille,et al.  Assimilation of satellite observations of long-lived chemical species in global chemistry transport models , 2000 .

[4]  Florence Rabier,et al.  Properties and first application of an error‐statistics tuning method in variational assimilation , 2004 .

[5]  M. eld.  Sequential Assimilation of Stratospheric Chemical Observations in a Three-Dimensional Model , 2002 .

[6]  Hermann Oelhaf,et al.  Remote sensing of vertical profiles of atmospheric trace constituents with MIPAS limb emission spectrometers , 1998, Asia-Pacific Environmental Remote Sensing.

[7]  Olivier Talagrand,et al.  Assimilation of Observations, an Introduction (gtSpecial IssueltData Assimilation in Meteology and Oceanography: Theory and Practice) , 1997 .

[8]  James M. Russell,et al.  The Halogen Occultation Experiment , 1993 .

[9]  R. Daley.  Atmospheric Data Analysis , 1991 .

[10]  D. P. DEE,et al.  Bias and data assimilation , 2005 .

[11]  O. Talagrand,et al.  Diagnosis and tuning of observational error in a quasi‐operational data assimilation setting , 2006 .

[12]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

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

[14]  J. R. Eyre,et al.  The optimal density of atmospheric sounder observations in the Met Office NWP system , 2007 .

[15]  J. Nocedal.  Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[16]  S. Cohn.  Dynamics of Short-Term Univariate Forecast Error Covariances , 1993 .

[17]  Y. Trémolet.  Accounting for an imperfect model in 4D‐Var , 2006 .

[18]  A. Segers,et al.  Assimilation of GOME ozone profiles and a global chemistry–transport model using a Kalman filter with anisotropic covariance , 2005 .

[19]  M. McCormick,et al.  New studies of SAGE II and HALOE ozone profile and long-term change comparisons , 2005 .

[20]  A. Hollingsworth,et al.  The statistical structure of short-range forecast errors as determined from radiosonde data Part II: The covariance of height and wind errors , 1986 .

[21]  Clive D Rodgers,et al.  Microwindow selection for high-spectral-resolution sounders. , 2002, Applied optics.

[22]  S. Pawson,et al.  Monitoring of Observation Errors in the Assimilation of Satellite Ozone Data , 2004 .

[23]  J. Lelieveld,et al.  In-situ measurement of reactive hydrocarbons at Hohenpeissenberg with comprehensive two-dimensional gas chromatography (GC×GC-FID): use in estimating HO and NO 3 , 2007 .

[24]  Costas A. Varotsos,et al.  Geophysical validation of MIPAS-ENVISAT operational ozone data , 2007 .

[25]  T. Clarmann,et al.  MIPAS: an instrument for atmospheric and climate research , 2007 .

[26]  Li Wei,et al.  Estimation of Three-Dimensional Error Covariances. Part II: Analysis of Wind Innovation Vectors , 2001 .

[27]  Clive D Rodgers,et al.  Inverse Methods for Atmospheric Sounding: Theory and Practice , 2000 .

[28]  D. Murtagh,et al.  A study of polar ozone depletion based on sequential assimilation of satellite data from the ENVISAT/MIPAS and Odin/SMR instruments , 2006 .

[29]  Manuel López-Puertas,et al.  MIPAS level 2 operational analysis , 2006 .

[30]  Hendrik Elbern,et al.  Emission rate and chemical state estimation by 4-dimensional variational inversion , 2007 .

[31]  H. Oelhaf,et al.  Remote sensing of vertical profiles of atmospheric trace constituents with MlPAS limb-emission spectrometers. , 1996, Applied Optics.

[32]  G. Brasseur,et al.  Assimilation of MIPAS observations using a three‐dimensional global chemistry‐transport model , 2005 .

[33]  Anthony Hollingsworth,et al.  The statistical structure of short-range forecast errors as determined from radiosonde data , 1986 .

[34]  Michael Buchhold,et al.  The Operational Global Icosahedral-Hexagonal Gridpoint Model GME: Description and High-Resolution Tests , 2002 .

[35]  Massimo Bonavita,et al.  The ensemble Kalman filter in an operational regional NWP system: preliminary results with real observations , 2008 .

[36]  E. Kalnay,et al.  Simultaneous estimation of covariance inflation and observation errors within an ensemble Kalman filter , 2009 .

[37]  P. Courtier,et al.  Correlation modelling on the sphere using a generalized diffusion equation , 2001 .

[38]  S. Bekki,et al.  Data assimilation of stratospheric ozone using a high‐resolution transport model , 2002 .

[39]  Derek M. Cunnold,et al.  Assessment of SAGE version 6.1 ozone data quality , 2002 .

[40]  J. Schwinger,et al.  Chemical state estimation for the middle atmosphere by four‐dimensional variational data assimilation: System configuration , 2010 .

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

[42]  Peter M. Lyster,et al.  Assimilation of Stratospheric Chemical Tracer Observations Using a Kalman Filter. Part I: Formulation , 2000 .

[43]  Richard Ménard,et al.  Assimilation of Stratospheric Chemical Tracer Observations Using a Kalman Filter. Part II: χ2-Validated Results and Analysis of Variance and Correlation Dynamics , 2000 .