Optimal Estimation of the Climatological Mean

Balachandrudu Narapusetty George Mason University, Fairfax, Virginia, and Center for Ocean–Land–Atmosphere Studies, Calverton, Maryland

Search for other papers by Balachandrudu Narapusetty in
Current site
Google Scholar
PubMed
Close
,
Timothy DelSole George Mason University, Fairfax, Virginia, and Center for Ocean–Land–Atmosphere Studies, Calverton, Maryland

Search for other papers by Timothy DelSole in
Current site
Google Scholar
PubMed
Close
, and
Michael K. Tippett International Research Institute for Climate and Society, Palisades, New York

Search for other papers by Michael K. Tippett in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

This paper shows theoretically and with examples that climatological means derived from spectral methods predict independent data with less error than climatological means derived from simple averaging. Herein, “spectral methods” indicates a least squares fit to a sum of a small number of sines and cosines that are periodic on annual or diurnal periods, and “simple averaging” refers to mean averages computed while holding the phase of the annual or diurnal cycle constant. The fact that spectral methods are superior to simple averaging can be understood as a straightforward consequence of overfitting, provided that one recognizes that simple averaging is a special case of the spectral method. To illustrate these results, the two methods are compared in the context of estimating the climatological mean of sea surface temperature (SST). Cross-validation experiments indicate that about four harmonics of the annual cycle are adequate, which requires estimation of nine independent parameters. In contrast, simple averaging of daily SST requires estimation of 366 parameters—one for each day of the year, which is a factor of 40 more parameters. Consistent with the greater number of parameters, simple averaging poorly predicts samples that were not included in the estimation of the climatological mean, compared to the spectral method. In addition to being more accurate, the spectral method also accommodates leap years and missing data simply, results in a greater degree of data compression, and automatically produces smooth time series.

Corresponding author address: Balachandrudu Narapusetty, Center for Ocean–Land–Atmosphere Studies, 4041 Powder Mill Rd., Suite 302, Calverton, MD 20705. Email: bala@cola.iges.org

Abstract

This paper shows theoretically and with examples that climatological means derived from spectral methods predict independent data with less error than climatological means derived from simple averaging. Herein, “spectral methods” indicates a least squares fit to a sum of a small number of sines and cosines that are periodic on annual or diurnal periods, and “simple averaging” refers to mean averages computed while holding the phase of the annual or diurnal cycle constant. The fact that spectral methods are superior to simple averaging can be understood as a straightforward consequence of overfitting, provided that one recognizes that simple averaging is a special case of the spectral method. To illustrate these results, the two methods are compared in the context of estimating the climatological mean of sea surface temperature (SST). Cross-validation experiments indicate that about four harmonics of the annual cycle are adequate, which requires estimation of nine independent parameters. In contrast, simple averaging of daily SST requires estimation of 366 parameters—one for each day of the year, which is a factor of 40 more parameters. Consistent with the greater number of parameters, simple averaging poorly predicts samples that were not included in the estimation of the climatological mean, compared to the spectral method. In addition to being more accurate, the spectral method also accommodates leap years and missing data simply, results in a greater degree of data compression, and automatically produces smooth time series.

Corresponding author address: Balachandrudu Narapusetty, Center for Ocean–Land–Atmosphere Studies, 4041 Powder Mill Rd., Suite 302, Calverton, MD 20705. Email: bala@cola.iges.org

Save
  • Akaike, H., 1969: Fitting autoregressive models for prediction. Ann. Inst. Stat. Math., 21 , 243247.

  • Barnett, T. P., 1981: On the nature and causes of the large-scale thermal variability in the central North Pacific Ocean. J. Phys. Oceanogr., 11 , 887904.

    • Search Google Scholar
    • Export Citation
  • Box, G. E. P., G. M. Jenkins, and G. C. Reinsel, 1994: Time Series Analysis: Forecasting and Control. 3rd ed., Prentice Hall, 598 pp.

  • Epstein, E., 1988: A spectral climatology. J. Climate, 1 , 88107.

  • Liu, W. T., W. Tang, and X. Xie, 2008: Wind power distribution over the ocean. Geophys. Res. Lett., 35 , L13808. doi:10.1029/2008GL034172.

    • Search Google Scholar
    • Export Citation
  • Press, W. H., B. P. Fannerty, S. A. Teukolsky, and W. T. Vetterling, 1986: Numerical Recipes: The Art of Scientific Computing. 1st ed. Cambridge University Press, 848 pp.

    • Search Google Scholar
    • Export Citation
  • Reynolds, R. W., T. M. Smith, C. Liu, D. B. Chelton, K. S. Casey, and M. G. Schlax, 2007: Daily high resolution blended analysis for sea surface temperature. J. Climate, 20 , 54735496.

    • Search Google Scholar
    • Export Citation
  • Straus, D. M., 1983: On the role of seasonal cycle. J. Atmos. Sci., 40 , 303313.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 529 178 14
PDF Downloads 346 90 5