Spectral Analysis of Multidimensional Stochastic Geophysical Models with an Application to Decadal ENSO Variability

Richard Kleeman Courant Institute of Mathematical Sciences, New York, New York

Search for other papers by Richard Kleeman in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

Simple linear models with additive stochastic forcing have been rather successful in explaining the observed spectrum of important climate variables. Motivated by this, the authors analyze the spectral character of such a general stochastic system of finite dimension. The spectral matrix is derived in the case that the spectrum is a linear combination of dynamical variables and their stochastic forcings. It is found that the most convenient basis for analysis is provided by the normal modes. In general the spectrum consists of two pieces. The first “diagonal” piece is a symmetric Lorentzian curve centered on the normal mode frequencies with breadth and strength determined by the normal mode dissipation. The second cross-spectrum piece derives usually from the coherency of the stochastic forcing of two different normal modes. The cross-spectrum is smaller in magnitude than the corresponding two diagonal pieces. This relative magnitude is controlled by the Wiener coherency, which is equal to the magnitude of the correlation of the stochastic forcings of different normal modes. This new analysis framework is studied in detail for the ENSO case for which a two-dimensional stochastically forced oscillator has been previously suggested as a minimal model. It is found that the observed spectrum is rather easily reproduced given appropriate dissipation. Further, it is found that the cross-spectrum results in a phase-dependent enhancement or suppression of frequencies smaller than the dominant ENSO frequency. This therefore provides a new mechanism for decadal ENSO variability. Since the cross-spectrum is phase dependent, the decadal variability generated has a distinctive spatial character. The significance of the cross-spectrum depends on the Wiener coherency, which in turn depends on the statistics of the stochastic forcing.

Corresponding author address: Richard Kleeman, Courant Institute of Mathematical Sciences, 251 Mercer St., New York, NY 10012. Email: kleeman@cims.nyu.edu

Abstract

Simple linear models with additive stochastic forcing have been rather successful in explaining the observed spectrum of important climate variables. Motivated by this, the authors analyze the spectral character of such a general stochastic system of finite dimension. The spectral matrix is derived in the case that the spectrum is a linear combination of dynamical variables and their stochastic forcings. It is found that the most convenient basis for analysis is provided by the normal modes. In general the spectrum consists of two pieces. The first “diagonal” piece is a symmetric Lorentzian curve centered on the normal mode frequencies with breadth and strength determined by the normal mode dissipation. The second cross-spectrum piece derives usually from the coherency of the stochastic forcing of two different normal modes. The cross-spectrum is smaller in magnitude than the corresponding two diagonal pieces. This relative magnitude is controlled by the Wiener coherency, which is equal to the magnitude of the correlation of the stochastic forcings of different normal modes. This new analysis framework is studied in detail for the ENSO case for which a two-dimensional stochastically forced oscillator has been previously suggested as a minimal model. It is found that the observed spectrum is rather easily reproduced given appropriate dissipation. Further, it is found that the cross-spectrum results in a phase-dependent enhancement or suppression of frequencies smaller than the dominant ENSO frequency. This therefore provides a new mechanism for decadal ENSO variability. Since the cross-spectrum is phase dependent, the decadal variability generated has a distinctive spatial character. The significance of the cross-spectrum depends on the Wiener coherency, which in turn depends on the statistics of the stochastic forcing.

Corresponding author address: Richard Kleeman, Courant Institute of Mathematical Sciences, 251 Mercer St., New York, NY 10012. Email: kleeman@cims.nyu.edu

Save
  • Benzi, R., G. Parisi, A. Sutera, and A. Vulpiani, 1982: Stochastic resonance in climatic change. Tellus, 34 , 1016.

  • DelSole, T., and B. F. Farrell, 1995: A stochastically excited linear system as a model for quasigeostrophic turbulence: Analytic results for one- and two-layer fluids. J. Atmos. Sci., 52 , 25312547.

    • Search Google Scholar
    • Export Citation
  • Flügel, M., P. Chang, and C. Penland, 2004: The role of stochastic forcing in modulating ENSO predictability. J. Climate, 17 , 31253140.

    • Search Google Scholar
    • Export Citation
  • Frankignoul, C., and K. Hasselmann, 1977: Stochastic climate models. Part II. Application to sea-surface temperature anomalies and thermocline variability. Tellus, 29 , 289305.

    • Search Google Scholar
    • Export Citation
  • Gardiner, C. W., 2004: Handbook of Stochastic Methods for Physics, Chemistry, and the Natural Sciences. Springer Series in Synergetics, Vol. 13, Springer, 415 pp.

    • Search Google Scholar
    • Export Citation
  • Hasselmann, K., 1976: Stochastic climate models. Part I. Theory. Tellus, 28 , 473485.

  • Horn, R., and C. Johnson, 1990: Matrix Analysis. Cambridge University Press, 575 pp.

  • Kleeman, R., 2002: Measuring dynamical prediction utility using relative entropy. J. Atmos. Sci., 59 , 20572072.

  • Kleeman, R., cited. 2008a: Stochastic climate models. [Available online at http://www.math.nyu.edu/faculty/kleeman/stochastic.ppt].

  • Kleeman, R., 2008b: Stochastic theories for the irregularity of ENSO. Philos. Trans. Roy. Soc., 366A , 25092524.

  • Kleeman, R., and A. M. Moore, 1999: A new method for determining the reliability of dynamical ENSO predictions. Mon. Wea. Rev., 127 , 694705.

    • Search Google Scholar
    • Export Citation
  • Kloeden, P. E., and E. Platen, 1992: Numerical Solution of Stochastic Differential Equations. Applications of Mathematics Series, Vol. 23, Springer, 636 pp.

    • Search Google Scholar
    • Export Citation
  • Landau, L., and E. Lifshitz, 1981: Mechanics. Vol. 1, 3rd ed. Course of Theoretical Physics, Reed Educational, 173 pp.

  • Latif, M., R. Kleeman, and C. Eckert, 1997: Greenhouse warming, decadal variability, or El Niño? An attempt to understand the anomalous 1990s. J. Climate, 10 , 22212239.

    • Search Google Scholar
    • Export Citation
  • Majda, A. J., C. Franzke, and B. Khouider, 2008: An applied mathematics perspective on stochastic modelling for climate. Philos. Trans. Roy. Soc. London, 366A , 24272453.

    • Search Google Scholar
    • Export Citation
  • Moore, A. M., and R. Kleeman, 1999: The nonnormal nature of El Niño and intraseasonal variability. J. Climate, 12 , 29652982.

  • Moore, A. M., J. Vialard, A. Weaver, D. L. T. Anderson, R. Kleeman, and J. R. Johnson, 2003: The role of air–sea interaction in controlling the optimal perturbations of low-frequency tropical coupled ocean–atmosphere modes. J. Climate, 16 , 951968.

    • Search Google Scholar
    • Export Citation
  • Newman, M., P. D. Sardeshmukh, and C. Penland, 1997: Stochastic forcing of the wintertime extratropical flow. J. Atmos. Sci., 54 , 435455.

    • Search Google Scholar
    • Export Citation
  • Penland, C., and M. Ghil, 1993: Forecasting Northern Hemisphere 700-mb geopotential height anomalies using empirical normal modes. Mon. Wea. Rev., 121 , 23552372.

    • Search Google Scholar
    • Export Citation
  • Penland, C., and L. Matrosova, 1994: A balance condition for stochastic numerical models with applications to the El Niño–Southern Oscillation. J. Climate, 7 , 13521372.

    • Search Google Scholar
    • Export Citation
  • Penland, C., and P. Sardeshmukh, 1995: The optimal growth of tropical sea surface temperature anomalies. J. Climate, 8 , 19992024.

  • Penland, C., and L. Matrosova, 2006: Studies of El Niño and interdecadal variability in tropical sea surface temperatures using a nonnormal filter. J. Climate, 19 , 57965815.

    • Search Google Scholar
    • Export Citation
  • Priestley, M., 1983: Spectral Analysis and Time Series. Vols. 1 and 2, Academic Press, 890 pp.

  • Risken, H., 1989: The Fokker–Plank Equation. 2nd ed. Springer-Verlag, 474 pp.

  • Salby, M., and R. Garcia, 1987: Transient response to localized episodic heating in the tropics. Part I: Excitation and short-time near-field behavior. J. Atmos. Sci., 44 , 458498.

    • Search Google Scholar
    • Export Citation
  • Saravanan, R., and J. McWilliams, 1998: Advective ocean–atmosphere interaction: An analytical stochastic model with implications for decadal variability. J. Climate, 11 , 165188.

    • Search Google Scholar
    • Export Citation
  • Thompson, C. J., and D. S. Battisti, 2001: A linear stochastic dynamical model of ENSO. Part II: Analysis. J. Climate, 14 , 445466.

  • Von Storch, H., T. Bruns, I. Fischer-Bruns, and K. Hasselmann, 1988: Principal oscillation pattern analysis of the 30- to 60-day oscillation in general circulation model equatorial troposphere. J. Geophys. Res., 93 , (D9). 1102211036.

    • Search Google Scholar
    • Export Citation
  • Wheeler, M., and G. Kiladis, 1999: Convectively coupled equatorial waves: Analysis of clouds and temperature in the wavenumber–frequency domain. J. Atmos. Sci., 56 , 374399.

    • Search Google Scholar
    • Export Citation
  • Zavala-Garay, J., A. Moore, C. L. Perez, and R. Kleeman, 2003: The response of a coupled model of ENSO to observed estimates of stochastic forcing. J. Climate, 16 , 28272842.

    • Search Google Scholar
    • Export Citation
  • Zhang, Yuan, J. Wallace, and D. Battisti, 1997: ENSO-like interdecadal variability: 1900–93. J. Climate, 10 , 10041020.

  • Zhang, Yunqing, and I. Held, 1999: A linear stochastic model of a GCM’s midlatitude storm tracks. J. Atmos. Sci., 56 , 34163435.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 238 101 9
PDF Downloads 123 47 8