Title: Deriving optimal weights for combining static and flow-dependent covariance models

Authors: Elizabeth Satterfield (Naval Research Laboratory, Monterey, CA)
Craig H. Bishop (Naval Research Laboratory, Monterey, CA)
David D. Kuhl (Naval Research Laboratory, Washington DC)
Tom Rosmond (Science Application International Corp., Forks, WA)

Hybrid 4D-VAR schemes employ forecast error covariance models that are linear combinations of static covariance models and flow-dependent ensemble based covariance models. Currently, to determine the optimal weights for the flow-dependent and static parts one must perform computationally expensive trials to compare each set of plausible weights. Since it is likely that the accuracy of ensemble based error covariances varies with vertical level and latitude, it also seems reasonable to assume that the optimal weights for a Hybrid covariance model would be a slowly varying function of latitude and height. However, the amount of tuning required to optimize such spatially varying weights would be prohibitively expensive. Here, we introduce and demonstrate a new method for deriving spatially varying weights for Hybrid error covariance models. This method is based on a univariate model of the stochastic relationship between error variance predictions and true error variance (discussed in detail in Bishop et al's presentation). If the assumptions of the univariate model were true, the approach would recover optimal weights for linear combinations of ensemble variances and climatological variances from archives of innovation and ensemble variance pairs. Our first experiment evaluates the assumptions of the univariate model applied to the Lorenz-40-variable model. The second set of experiments creates spatially varying weights for our newly developed NAVDAS-AR Hybrid by dividing the atmosphere into six regions: low-level N.Hemisphere, Tropics and S. Hemisphere and upper-level N. Hemisphere, Tropics and S. Hemisphere. The performance of the NAVDAS-AR Hybrid using these spatially varying weights is then compared against its performance using non-spatially varying weights derived by brute force tuning. At low resolution, we found that the spatially varying weights derived using the univariate theory match the performance the invariant weights obtained from tuning. If available, higher resolution results will also be presented.

GMAO Head: Michele Rienecker Global Modeling and Assimilation Office NASA Goddard Space Flight Center

Curator: Nikki Privé Last Updated: May 27 2011