Joanne A. Waller

Diagnosing Observation Error Correlations for Doppler Radar Radial Winds in the Met Office UKV Model Using Observation-Minus-Background and Observation-Minus-Analysis Statistics

AbstractWith the development of convection-permitting numerical weather prediction the efficient use of high-resolution observations in data assimilation is becoming increasingly important. The operational assimilation of these observations, such as Doppler radar radial winds (DRWs), is now common, although to avoid violating the assumption of uncorrelated observation errors the observation density is severely reduced. To improve the quantity of observations used and the impact that they have on the forecast requires the introduction of the full, potentially correlated, error statistics. In this work, observation error statistics are calculated for the DRWs that are assimilated into the Met Office high-resolution U.K. model (UKV) using a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. This is the first in-depth study using the diagnostic to estimate both horizontal and along-beam observation error statistics. The new results obtai...

Diagnosing horizontal and inter-channel observation error correlations for SEVIRI observations using observation-minus-background and observation-minus-analysis statistics

It has been common practice in data assimilation to treat observation errors as uncorrelated; however, meteorological centres are beginning to use correlated inter-channel observation errors in their operational assimilation systems. In this work, we are the first to characterise inter-channel and spatial error correlations for Spinning Enhanced Visible and Infrared Imager (SEVIRI) observations that are assimilated into the Met Office high-resolution model. The errors are calculated using a diagnostic that calculates statistical averages of observation-minus-background and observation-minus-analysis residuals. This diagnostic is sensitive to the background and observation error statistics used in the assimilation, although, with careful interpretation of the results, it can still provide useful information. We find that the diagnosed SEVIRI error variances are as low as one-tenth of those currently used in the operational system. The water vapour channels have significantly correlated inter-channel errors, as do the surface channels. The surface channels have larger observation error variances and inter-channel correlations in coastal areas of the domain; this is the result of assimilating mixed pixel (land-sea) observations. The horizontal observation error correlations range between 30 km and 80 km, which is larger than the operational thinning distance of 24 km. We also find that estimates from the diagnostics are unaffected by biased observations, provided that the observation-minus-background and observation-minus-analysis residual means are subtracted.

The conditioning of least squares problems in variational data assimilation

Summary In variational data assimilation a least-squares objective function is minimised to obtain the most likely state of a dynamical system. This objective function combines observation and prior (or background) data weighted by their respective error statistics. In numerical weather prediction, data assimilation is used to estimate the current atmospheric state, which then serves as an initial condition for a forecast. New developments in the treatment of observation uncertainties have recently been shown to cause convergence problems for this least-squares minimisation. This is important for operational numerical weather prediction centres due to the time constraints of producing regular forecasts. The condition number of the Hessian of the objective function can be used as a proxy to investigate the speed of convergence of the least-squares minimisation. In this paper we develop novel theoretical bounds on the condition number of the Hessian. These new bounds depend on the minimum eigenvalue of the observation error covariance matrix and the ratio of background error variance to observation error variance. Numerical tests in a linear setting show that the location of observation measurements has an important effect on the condition number of the Hessian. We identify that the conditioning of the problem is related to the complex interactions between observation error covariance and background error covariance matrices. Increased understanding of the role of each constituent matrix in the conditioning of the Hessian will prove useful for informing the choice of correlated observation error covariance matrix and observation location, particularly for practical applications.