Improved Vertical Covariance Estimates for Ensemble-Filter Assimilation of Near-Surface Observations

Hacker, Joshua P.; Anderson, Jeffrey L.; Pagowski, Mariusz
March 2007
Monthly Weather Review;Mar2007, Vol. 135 Issue 3, p1021
Academic Journal
Strategies to improve covariance estimates for ensemble-based assimilation of near-surface observations in atmospheric models are explored. It is known that localization of covariance estimates can improve conditioning of covariance matrices in the assimilation process by removing spurious elements and increasing the rank of the matrix. Vertical covariance localization is the focus of this work, and two basic approaches are compared: 1) a recently proposed hierarchical filter approach based on sampling theory and 2) a more commonly used fifth-order piecewise rational function. The hierarchical filter allows for dynamic estimates of localization functions and does not place any restrictions on their form. The rational function is optimized for every analysis time of day and for every possible observation and state variable combination. The methods are tested with a column model containing PBL and land surface parameterization schemes that are available in current mesoscale modeling systems. The results are expected to provide context for assimilation of near-surface observations in mesoscale models, which will benefit short-range mesoscale NWP applications. Results show that both the hierarchical and rational function approaches effectively improve covariance estimates from small ensembles. The hierarchical approach provides localization functions that are irregular and more closely related to PBL structure. Analysis of eigenvalue spectra show that both approaches improve the rank of the covariance matrices, but the amount of improvement is not always directly related to the assimilation performance. Results also show that specifying different localization functions for different observation and state variable combinations is more important than including time dependence.


Related Articles

  • Assessing the pattern of covariance matrices via an augmentation multiple testing procedure. Greselin, Francesca; Ingrassia, Salvatore; Punzo, Antonio // Statistical Methods & Applications;2011, Vol. 20 Issue 2, p141 

    This paper extends the scedasticity comparison among several groups of observations, usually complying with the homoscedastic and the heteroscedastic cases, in order to deal with data sets laying in an intermediate situation. As is well known, homoscedasticity corresponds to equality in...

  • Partial linear single index models with distortion measurement errors. Zhang, Jun; Yu, Yao; Zhu, Li-Xing; Liang, Hua // Annals of the Institute of Statistical Mathematics;Apr2013, Vol. 65 Issue 2, p237 

    We study partial linear single index models when the response and the covariates in the parametric part are measured with errors and distorted by unknown functions of commonly observable confounding variables, and propose a semiparametric covariate-adjusted estimation procedure. We apply the...

  • On the Mátern covariance family: a proposal for modeling temporal correlations based on turbulence theory. Kermarrec, Gaël; Schön, Steffen // Journal of Geodesy;Nov2014, Vol. 88 Issue 11, p1061 

    Current variance models for GPS carrier phases that take correlation due to tropospheric turbulence into account are mathematically difficult to handle due to numerical integrations. In this paper, a new model for temporal correlations of GPS phase measurements based on turbulence theory is...

  • Nonparametric estimation of large covariance matrices of longitudinal data. Wei Biao Wu; Mohsen Pourahmadi // Biometrika;Dec2003, Vol. 90 Issue 4, p831 

    Estimation of an unstructured covariance matrix is difficult because of its positive-definiteness constraint. This obstacle is removed by regressing each variable on its predecessors, so that estimation of a covariance matrix is shown to be equivalent to that of estimating a sequence of...

  • Efficient propagation of error covariance matrices of gravitational models: application to GRACE and GOCE. Georges Balmino // Journal of Geodesy;Oct2009, Vol. 83 Issue 10, p989 

    Abstract  We have applied efficient methods for computing variances and covariances of functions of a global gravity field model expanded in spherical harmonics, using the full variance–covariance matrix of the coefficients. Examples are given with recent models derived from GRACE...

  • PROFILE ANALYSIS WITH RANDOM-EFFECTS COVARIANCE STRUCTURE. Srivastava, Muni S.; Singull, Martin // Journal of the Japan Statistical Society;2012, Vol. 42 Issue 2, p145 

    In this paper, we consider a parallel profile model for several groups. Given the parallel profile model we construct tests based on the likelihood ratio, without any restrictions on the parameter space, for testing the covariance matrix for random-effects structure or sphericity. Furthermore,...

  • Improvement in background error covariances using ensemble forecasts for assimilation of high-resolution satellite data. Seung-Woo LEE; Dong-Kyou LEE // Advances in Atmospheric Sciences;Jul2011, Vol. 28 Issue 4, p758 

    Satellite data obtained over synoptic data-sparse regions such as an ocean contribute toward improving the quality of the initial state of limited-area models. Background error covariances are crucial to the proper distribution of satellite-observed information in variational data assimilation....

  • Efficient estimation for the Cox model with varying coefficients. Chen, Kani; Lin, Huazhen; Zhou, Yong // Biometrika;Jun2012, Vol. 99 Issue 2, p379 

    A proportional hazards model with varying coefficients allows one to examine the extent to which covariates interact nonlinearly with an exposure variable. A global partial likelihood method, in contrast with the local partial likelihood method of Fan et al. (2006), is proposed for estimation of...

  • A conditional model for incomplete covariates in parametric regression models. LIPSITZ, STUART R.; IBRAHIM, JOSEPH G. // Biometrika;1996, Vol. 83 Issue 4, p916 

    Incomplete covariate data arise in many data sets. When the missing covariates are categorical, a useful technique for obtaining parameter estimates is the EM algorithm by the method of weights proposed in Ibrahim (1990). This method requires the estimation of many nuisance parameters for the...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics