Extended assimilation and forecast experiments with a four‐dimensional variational assimilation system

Results of four‐dimensional variational assimilations, 4D‐Var, in cycling mode, over a few two‐week assimilation periods are presented. 4D‐Var is implemented in its incremental formulation, with a high‐resolution model with the full physical parametrization package to compare the atmospheric states with the observations, and a low‐resolution model with simplified physics to minimize the cost‐function. The comparison of 4D‐Var using several assimilation windows (6, 12 and 24 hours) with 3D‐Var (the equivalent of 4D‐Var with no time‐dimension) over a two‐week period shows a clear benefit from using 4D‐Var over a 6 or 12—hour window compared to the static 3D‐Var scheme. It also exhibits some problems with the forecasts started using 4D‐Var over a 24‐hour window. The poorer performance of 4D‐Var over a relatively long assimilation window can be partly explained by the fact that, in these experiments, the tangent‐linear and adjoint models used in the minimization are only approximations of the assimilating model (having lower resolution and crude physics). The error these approximations introduce in the time evolution of a perturbation affects the convergence of the incremental 4D‐Var, with larger discontinuities in the values of the cost‐function when going from low to high resolution for longer assimilation windows. Additional experiments are performed comparing 4D‐Var using a 6‐hour window with the 3D‐Var system. Two additional 2‐week periods show a consistent improvement in extratropical forecast scores with the 4D‐Var system. The main 4D‐Var improvements occur in areas where the 3D‐Var errors were the largest. Local improvement can be as large as 35% for the root‐mean‐square of the 5‐day‐forecast error, averaged over a two‐week period. A comparison of key analysis errors shows that, indeed, 4D‐Var using a 6‐hour window is able to reduce substantially the amplitude of its fast‐growing error components. The overall fit to observations of analyses and short‐range forecasts from 3D‐Var and 4D‐Var is comparable. In active baroclinic areas, the fit of the background to the data is considerably better for the 4D‐Var system, resulting in smaller increments. It appears that in these areas (and in particular over the west Atlantic), 4D‐Var is able to better use the information contained in the observations. The ability of 4D‐Var to extrapolate some aircraft data in the vertical with a baroclinic tilt is illustrated. Problems exist in the tropics and mountainous areas due partly to a lack of physics in the tangent‐linear model. Possible improvements to the system (the introduction of more physics; better behaviour of the incremental approach owing to a line search at high resolution) are also discussed.

[1]  J. M. Lewis,et al.  The use of adjoint equations to solve a variational adjustment problem with advective constraints , 1985 .

[2]  F. L. Dimet,et al.  Variational algorithms for analysis and assimilation of meteorological observations: theoretical aspects , 1986 .

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

[4]  Claude Lemaréchal,et al.  Some numerical experiments with variable-storage quasi-Newton algorithms , 1989, Math. Program..

[5]  J. Derber A Variational Continuous Assimilation Technique , 1989 .

[6]  M. Ghil,et al.  Data assimilation in meteorology and oceanography , 1991 .

[7]  P. Courtier,et al.  Four-dimensional variational data assimilation using the adjoint of a multilevel primitive-equation model , 1991 .

[8]  Philippe Courtier,et al.  Four‐Dimensional Assimilation In the Presence of Baroclinic Instability , 1992 .

[9]  R. Daley Estimating Model-Error Covariances for Application to Atmospheric Data Assimilation , 1992 .

[10]  J. Derber,et al.  Variational Data Assimilation with an Adiabatic Version of the NMC Spectral Model , 1992 .

[11]  Milija Zupanski,et al.  Regional Four-Dimensional Variational Data Assimilation in a Quasi-Operational Forecasting Environment , 1993 .

[12]  Philippe Courtier,et al.  Interactions of Dynamics and Observations in a Four-Dimensional Variational Assimilation , 1993 .

[13]  P. Courtier,et al.  Assimilation of Simulated Wind Lidar Data with a Kalman Filter , 1993 .

[14]  Philippe Courtier,et al.  Use of cloud‐cleared radiances in three/four‐dimensional variational data assimilation , 1994 .

[15]  A. Bennett,et al.  TOPEX/POSEIDON tides estimated using a global inverse model , 1994 .

[16]  Roberto Buizza,et al.  Sensitivity of optimal unstable structures , 1994 .

[17]  P. Courtier,et al.  A strategy for operational implementation of 4D‐Var, using an incremental approach , 1994 .

[18]  Ricardo Todling,et al.  A fixed-lag Kalman smoother for retrospective data assimilation , 1994 .

[19]  Michael Ghil,et al.  Advanced data assimilation in strongly nonlinear dynamical systems , 1994 .

[20]  R. Todling,et al.  Suboptimal Schemes for Atmospheric Data Assimilation Based on the Kalman Filter , 1994 .

[21]  Ying-Hwa Kuo,et al.  Assimilation of Atmospheric Radio Refractivity Using a Nonhydrostatic Adjoint Model , 1995 .

[22]  D. Dee On-line Estimation of Error Covariance Parameters for Atmospheric Data Assimilation , 1995 .

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

[24]  David J. Lary,et al.  Lagrangian four‐dimensional variational data assimilation of chemical species , 1995 .

[25]  L. Riishojgaard On four‐dimensional variational assimilation of ozone data in weather‐prediction models , 1996 .

[26]  Ying-Hwa Kuo,et al.  Variational Assimilation of Precipitable Water Using a Nonhydrostatic Mesoscale Adjoint Model. Part I: Moisture Retrieval and Sensitivity Experiments , 1996 .

[27]  P. Courtier,et al.  Quasi‐continuous variational data assimilation , 1996 .

[28]  Philippe Courtier,et al.  Dynamical structure functions in a four‐dimensional variational assimilation: A case study , 1996 .

[29]  Philippe Courtier,et al.  Sensitivity of forecast errors to initial conditions , 1996 .

[30]  Michael Ghil,et al.  Tracking Atmospheric Instabilities with the Kalman Filter. Part II: Two-Layer Results , 1996 .

[31]  D. Zupanski A General Weak Constraint Applicable to Operational 4DVAR Data Assimilation Systems , 1997 .

[32]  Philippe Courtier,et al.  Unified Notation for Data Assimilation : Operational, Sequential and Variational , 1997 .

[33]  P. Courtier,et al.  The ECMWF implementation of three‐dimensional variational assimilation (3D‐Var). II: Structure functions , 1998 .

[34]  P. Courtier,et al.  The ECMWF implementation of three‐dimensional variational assimilation (3D‐Var). I: Formulation , 1998 .

[35]  Philippe Courtier,et al.  The ECMWF implementation of three-dimensional variational assimilation ( 3 D-Var ) . 111 : Experimental results , 2006 .

[36]  P. Courtier,et al.  Four‐dimensional variational assimilation of SSM/I precipitable water content data , 1998 .

[37]  R. Gelaro,et al.  Estimation of key analysis errors using the adjoint technique , 1998 .