Construction and application of covariance functions with variable length‐fields

This article focuses on construction, directly in physical space, of three-dimensional covariance functions parametrized by a length-field, and on an application of these functions to improve the representation of the Quasi-Biennial Oscillation (QBO) in the Goddard Earth Observing System, Version 4 (GEOS-4) data assimilation system. The covariance functions are obtained by fusing collections of auto-covariance functions having different constant length-scales with their associated cross-covariance functions. This construction yields covariance functions with length-scales that can vary arbitrarily over any finite partition of the spatial domain. A simple, and also motivating application of these functions is to the case where the length-scale varies in the vertical direction only. The class of covariance functions with variable length-fields constructed in this article will be called multi-level to associate them with this application. The multi-level covariance functions extend well-known single-level covariance functions depending only on a constant length-scale. Generalizations of the familiar first-and third-order autoregressive covariances in three dimensions are given, providing multi-level covariances with zero and four continuous derivatives at zero separation, respectively. Multi-level piecewise rational covariances with two continuous derivatives at zero separation are also provided. Multi-level power-law covariances are constructed with continuous derivatives of all orders. Additional multi-level covariance functions are constructed using the Schur product of single- and multi-level covariance functions. A multi-variate, multi-level power-law covariance with a large troposphere-to-stratosphere length-field gradient is employed to reproduce the QBO from sparse radiosonde wind observations in the tropical lower stratosphere. This covariance model is described along with details of the assimilation experiments. The new covariance model is shown to represent the vertical wind shear associated with the QBO much more effectively than the multi-variate, multi-level covariance model in the baseline GEOS-4 system. Copyright © 2006 Royal Meteorological Society

[1]  E. J. Hannan,et al.  Multiple time series , 1970 .

[2]  A. Segers,et al.  Assimilation of GOME ozone profiles and a global chemistry–transport model using a Kalman filter with anisotropic covariance , 2005 .

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

[4]  S. Cohn,et al.  Assessing the Effects of Data Selection with the DAO Physical-Space Statistical Analysis System* , 1998 .

[5]  H. Mitchell,et al.  A continuous three‐dimensional model of short‐range forecast error covariances , 1992 .

[6]  J. Holton.  An introduction to dynamic meteorology , 2004 .

[7]  Michael G. Bosilovich,et al.  Documentation and Validation of the Goddard Earth Observing System (GEOS) Data Assimilation System, Version 4 , 2005 .

[8]  Norman A. Phillips,et al.  The spatial statistics of random geostrophic modes and first-guess errors , 1986 .

[9]  A. Hollingsworth,et al.  The statistical structure of short-range forecast errors as determined from radiosonde data Part II: The covariance of height and wind errors , 1986 .

[10]  Philippe Keckhut,et al.  The SPARC Intercomparison of Middle-Atmosphere Climatologies , 2004, Journal of Climate.

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

[12]  Roger Daley,et al.  NAVDAS: Formulation and Diagnostics , 2001 .

[13]  G. Stenchikov,et al.  Longitudinal Variation of the Stratospheric Quasi-Biennial Oscillation , 2004 .

[14]  Harald Cramer,et al.  On the Theory of Stationary Random Processes , 1940 .

[15]  Kevin Hamilton,et al.  The quasi‐biennial oscillation , 2001 .

[16]  Peter M. Lyster,et al.  Assimilation of Stratospheric Chemical Tracer Observations Using a Kalman Filter. Part I: Formulation , 2000 .

[17]  S. Cohn,et al.  Ooce Note Series on Global Modeling and Data Assimilation Construction of Correlation Functions in Two and Three Dimensions and Convolution Covariance Functions , 2022 .

[18]  M. Fiorino,et al.  A comparison of reanalyses in the tropical stratosphere. Part 2: the quasi-biennial oscillation , 1998 .

[19]  John Derber,et al.  The National Meteorological Center's spectral-statistical interpolation analysis system , 1992 .

[20]  A. D. Belmont,et al.  VARIATION WITH LONGITUDE OF THE QUASI-BIENNIAL OSCILLATION1,,2 , 1968 .

[21]  G. Folland.  Fourier analysis and its applications , 1992 .

[22]  R. Daley.  Atmospheric Data Analysis , 1991 .

[23]  Anthony Hollingsworth,et al.  The statistical structure of short-range forecast errors as determined from radiosonde data , 1986 .

[24]  Nedjeljka Žagar,et al.  Balanced tropical data assimilation based on a study of equatorial waves in ECMWF short‐range forecast errors , 2005 .

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