Numerical Aspects of the Application of Recursive Filters to Variational Statistical Analysis. Part II: Spatially Inhomogeneous and Anisotropic General Covariances

R. James Purser Science Applications International Corporation, Beltsville, Maryland

Search for other papers by R. James Purser in
Current site
Google Scholar
PubMed
Close
,
Wan-Shu Wu Science Applications International Corporation, Beltsville, Maryland

Search for other papers by Wan-Shu Wu in
Current site
Google Scholar
PubMed
Close
,
David F. Parrish National Centers for Environmental Prediction/Environmental Modeling Center, Washington, D.C

Search for other papers by David F. Parrish in
Current site
Google Scholar
PubMed
Close
, and
Nigel M. Roberts Joint Centre for Mesoscale Meteorology, University of Reading, Reading, and Met Office, Bracknell, United Kingdom

Search for other papers by Nigel M. Roberts in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

In this second part of a two-part study of recursive filter techniques applied to the synthesis of covariances in a variational analysis, methods by which non-Gaussian shapes and spatial inhomogeneities and anisotropies for the covariances may be introduced in a well-controlled way are examined. These methods permit an analysis scheme to possess covariance structures with adaptive variations of amplitude, scale, profile shape, and degrees of local anisotropy, all as functions of geographical location and altitude.

First, it is shown how a wider and more useful variety of covariance shapes than just the Gaussian may be obtained by the positive superposition of Gaussian components of different scales, or by further combinations of these operators with the application of Laplacian operators in order for the products to possess negative sidelobes in their radial profiles.

Then it is shown how the techniques of recursive filters may be generalized to admit the construction of covariances whose characteristic scales relative to the grid become adaptive to geographical location, while preserving the necessary properties of self-adjointness and positivity. Special attention is paid to the problems of amplitude control for these spatially inhomogeneous filters and an estimate for the kernel amplitude is proposed based upon an asymptotic analysis of the problem.

Finally, a further generalization of the filters that enables fully anisotropic and geographically adaptive covariances to be constructed in a computationally efficient way is discussed.

Corresponding author address: R. James Purser, NOAA/NCEP, W/NP2 Rm. 207, WWBG, 5200 Auth Rd., Camp Springs, MD 20746-4304. Email: jim.purser@noaa.gov

Abstract

In this second part of a two-part study of recursive filter techniques applied to the synthesis of covariances in a variational analysis, methods by which non-Gaussian shapes and spatial inhomogeneities and anisotropies for the covariances may be introduced in a well-controlled way are examined. These methods permit an analysis scheme to possess covariance structures with adaptive variations of amplitude, scale, profile shape, and degrees of local anisotropy, all as functions of geographical location and altitude.

First, it is shown how a wider and more useful variety of covariance shapes than just the Gaussian may be obtained by the positive superposition of Gaussian components of different scales, or by further combinations of these operators with the application of Laplacian operators in order for the products to possess negative sidelobes in their radial profiles.

Then it is shown how the techniques of recursive filters may be generalized to admit the construction of covariances whose characteristic scales relative to the grid become adaptive to geographical location, while preserving the necessary properties of self-adjointness and positivity. Special attention is paid to the problems of amplitude control for these spatially inhomogeneous filters and an estimate for the kernel amplitude is proposed based upon an asymptotic analysis of the problem.

Finally, a further generalization of the filters that enables fully anisotropic and geographically adaptive covariances to be constructed in a computationally efficient way is discussed.

Corresponding author address: R. James Purser, NOAA/NCEP, W/NP2 Rm. 207, WWBG, 5200 Auth Rd., Camp Springs, MD 20746-4304. Email: jim.purser@noaa.gov

Save
  • Benjamin, S. G., 1989: An isentropic meso-α-scale analysis system and its sensitivity to aircraft and surface observations. Mon. Wea. Rev., 117 , 15861603.

    • Search Google Scholar
    • Export Citation
  • Coxeter, H. S. M., 1973: Regular Polytopes. Dover, 321 pp.

  • da Silva, A., J. Pfaendtner, J. Guo, M. Sienkiewicz, and S. Cohn, 1995: Assessing the effects of data selection with DAO's physical-space statistical analysis system. Second International Symposium on Assimilation of Observations in Meteorology and Oceanography, WMO Tech. Rep. TD 651, 273–278.

    • Search Google Scholar
    • Export Citation
  • Dee, D. P., and A. M. da Silva, 1999: Maximum-likelihood estimation of forecast and observation error covariance parameters. Part I: Methodology. Mon. Wea. Rev., 127 , 18221834.

    • Search Google Scholar
    • Export Citation
  • Derber, J., and A. Rosati, 1989: A global oceanic data assimilation system. J. Phys. Oceanogr., 19 , 13331347.

  • Desroziers, G., 1997: A coordinate change for data assimilation in spherical geometry of frontal structure. Mon. Wea. Rev., 125 , 30303038.

    • Search Google Scholar
    • Export Citation
  • Fujita, T., and R. J. Purser, 2001: Parallel implementation of compact numerical schemes. NOAA/NCEP Office Note 434, 34 pp. [Available from NCEP, 5200 Auth Rd., Camp Springs, MD 20746.].

    • Search Google Scholar
    • Export Citation
  • Gaspari, G., and S. E. Cohn, 1998: Construction of correlation functions in two and three dimensions. Office Note Series on Global Modeling and Data Assimilation, DAO Office Note 96-03R1, DAO, GSFC, 53 pp.

    • Search Google Scholar
    • Export Citation
  • Gaspari, G., and S. E. Cohn, 1999: Construction of correlation functions in two and three dimensions. Quart. J. Roy. Meteor. Soc., 125 , 723757.

    • Search Google Scholar
    • Export Citation
  • Gneiting, T., 1999: Correlation functions for atmospheric data analysis. Quart. J. Roy. Meteor. Soc., 125 , 24492464.

  • Hayden, C. M., and R. J. Purser, 1995: Recursive filter objective analysis of meteorological fields: Applications to NESDIS operational processing. J. Appl. Meteor., 34 , 315.

    • Search Google Scholar
    • Export Citation
  • Lorenc, A. C., 1981: A global three-dimensional multivariate statistical interpolation scheme. Mon. Wea. Rev., 109 , 701721.

  • Lorenc, A. C., 1997: Development of an operational variational assimilation scheme. J. Meteor. Soc. Japan, 75 , 339346.

  • Otte, T. L., N. L. Seaman, and D. R. Stauffer, 2001: A heuristic study on the importance of anisotropic error distributions in data assimilation. Mon. Wea. Rev., 129 , 766783.

    • Search Google Scholar
    • Export Citation
  • Parrish, D. F., and R. J. Purser, 1998: Anisotropic covariances in 3DVAR: Application to hurricane Doppler radar observations. Proc. HIRLAM Workshop on Variational Analysis, Toulouse, France, Météo-France, 57–65. [Available from Met Eireann, Glasnevin Hill, Dublin 9, Ireland.].

    • Search Google Scholar
    • Export Citation
  • Press, W. H., S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, 1992: Numerical Recipes in Fortran 77. 2d ed. Cambridge University Press, 933 pp.

    • Search Google Scholar
    • Export Citation
  • Purser, R. J., and R. McQuigg, 1982: A successive correction analysis scheme using recursive numerical filters. Met Office Tech. Note 154, British Meteorological Office, 17 pp.

    • Search Google Scholar
    • Export Citation
  • Purser, R. J., and D. F. Parrish, 2003: A Bayesian technique for estimating continuously varying statistical parameters of a variational assimilation. Meteor. Atmos. Phys., 82 , 209226.

    • Search Google Scholar
    • Export Citation
  • Purser, R. J., W-S. Wu, D. F. Parrish, and N. M. Roberts, 2001: Numerical aspects of the application of recursive filters to variational statistical analysis. NOAA/NCEP Office Note 431, 34 pp. [Available from NCEP, 5200 Auth Rd., Camp Springs, MD 20746.].

    • Search Google Scholar
    • Export Citation
  • Purser, R. J., W-S. Wu, D. F. Parrish, and N. M. Roberts, 2003: Numerical aspects of the application of recursive filters to variational statistical analysis. Part I: Spatially homogeneous and isotropic Gaussian covariances. Mon. Wea. Rev., 131 , 15241535.

    • Search Google Scholar
    • Export Citation
  • Riishøjgaard, L-P., 1998: A direct way of specifying flow-dependent background error correlations for meteorological analysis systems. Tellus, 50A , 4257.

    • Search Google Scholar
    • Export Citation
  • Schoenberg, I. J., 1938: Metric spaces and completely monotone functions. Ann. Math., 39 , 811841.

  • Shapiro, M. A., and J. T. Hastings, 1973: Objective cross section analysis by Hermite polynomial interpolation on isentropic surfaces. J. Appl. Meteor., 12 , 753762.

    • Search Google Scholar
    • Export Citation
  • Thépaut, J-N., P. Courtier, G. Belaud, and G. Lemaître, 1996: Dynamical structure functions in a four-dimensional variational assimilation: A case study. Quart. J. Roy. Meteor. Soc., 122 , 535561.

    • Search Google Scholar
    • Export Citation
  • Weaver, A., and P. Courtier, 2001: Correlation modelling on the sphere using a generalized diffusion equation. Quart. J. Roy. Meteor. Soc., 127 , 18151846.

    • Search Google Scholar
    • Export Citation
  • Wu, W-S., R. J. Purser, and D. F. Parrish, 2002: Three-dimensional variational analysis with spatially inhomogeneous covariances. Mon. Wea. Rev., 130 , 29052916.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1057 402 85
PDF Downloads 310 115 14