Editorial Type:
Article
Article Type:
Research Article

Application of a Reduced-Order Kalman Filter to Initialize a Coupled Atmosphere–Ocean Model: Impact on the Prediction of El Niño

Joaquim Ballabrera-Poy Universities Space Research Association, NASA Goddard Space Flight Center, Greenbelt, Maryland

Search for other papers by Joaquim Ballabrera-Poy in
Current site
Google Scholar
PubMed Close
,
Antonio J. Busalacchi Laboratory for Hydrospheric Processes, NASA Goddard Space Flight Center, Greenbelt, Maryland

Search for other papers by Antonio J. Busalacchi in
Current site
Google Scholar
PubMed Close
, and
Ragu Murtugudde ESSIC, University of Maryland, College Park, Maryland, and NASA Goddard Space Flight Center, Greenbelt, Maryland

Search for other papers by Ragu Murtugudde in
Current site
Google Scholar
PubMed Close
Print Publication:
15 Apr 2001
DOI:
https://doi.org/10.1175/1520-0442(2001)014<1720:AOAROK>2.0.CO;2
Page(s):
1720–1737
Restricted access

Abstract

A reduced-order Kalman filter is used to assimilate observed fields of the surface wind stress, sea surface temperature, and sea level into the coupled ocean–atmosphere model of Zebiak and Cane. The method projects the Kalman filter equations onto a subspace defined by the eigenvalue decomposition of the error forecast matrix, allowing its application to high-dimensional systems.

The Zebiak and Cane model couples a linear, reduced-gravity ocean model with a single, vertical-mode atmospheric model. The compatibility between the simplified physics of the model and each observed variable is studied separately and together. The results show the ability of the empirical orthogonal functions (EOFs) of the model to represent the simultaneous value of the wind stress, SST, and sea level, when the fields are limited to the latitude band 10°S–10°N, and when the number of EOFs is greater than the number of statistically significant modes.

In this first application of the Kalman filter to a coupled ocean–atmosphere prediction model, the sea level fields are assimilated in terms of the Kelvin and Rossby modes of the thermocline depth anomaly. An estimation of the error of these modes is derived from the projection of an estimation of the sea level error over such modes.

The ability of the method to reconstruct the state of the equatorial Pacific and to predict its time evolution is shown. The method is quite robust for predictions up to 6 months, and able to predict the onset of the 1997 warm event 15 months before its occurrence.

Corresponding author address: Joaquim Ballabrera-Poy, NASA Goddard Space Flight Center, Mailcode 970, Greenbelt, MD 20771.

Email: joaquim@neptune.gsfc.nasa.gov

Abstract

A reduced-order Kalman filter is used to assimilate observed fields of the surface wind stress, sea surface temperature, and sea level into the coupled ocean–atmosphere model of Zebiak and Cane. The method projects the Kalman filter equations onto a subspace defined by the eigenvalue decomposition of the error forecast matrix, allowing its application to high-dimensional systems.

The Zebiak and Cane model couples a linear, reduced-gravity ocean model with a single, vertical-mode atmospheric model. The compatibility between the simplified physics of the model and each observed variable is studied separately and together. The results show the ability of the empirical orthogonal functions (EOFs) of the model to represent the simultaneous value of the wind stress, SST, and sea level, when the fields are limited to the latitude band 10°S–10°N, and when the number of EOFs is greater than the number of statistically significant modes.

In this first application of the Kalman filter to a coupled ocean–atmosphere prediction model, the sea level fields are assimilated in terms of the Kelvin and Rossby modes of the thermocline depth anomaly. An estimation of the error of these modes is derived from the projection of an estimation of the sea level error over such modes.

The ability of the method to reconstruct the state of the equatorial Pacific and to predict its time evolution is shown. The method is quite robust for predictions up to 6 months, and able to predict the onset of the 1997 warm event 15 months before its occurrence.

Corresponding author address: Joaquim Ballabrera-Poy, NASA Goddard Space Flight Center, Mailcode 970, Greenbelt, MD 20771.

Email: joaquim@neptune.gsfc.nasa.gov

Save
Cover Journal of Climate
Journal of Climate

  • Barnett, T. P., N. E. Graham, M. A. Cane, S. E. Zebiak, S. C. Dolan, J. J. O’Brien, and D. M. Legler, 1988: On the prediction of the El Niño of 1986–1987. Science,241, 192–196.

  • ——, M. Latif, N. E. Graham, M. Flügel, S. Pazan, and W. White, 1993: ENSO and ENSO-related predictability. Part I: Prediction of equatorial Pacific sea surface temperature with a hybrid coupled ocean–atmosphere model. J. Climate,6, 1545–1566.

  • Barnston, A. G., M. H. Glantz, and Y. He, 1999: Predictive skill of statistical and dynamical climate models in SST forecasts during the 1997–98 El Niño episode and the 1998 La Niña onset. Bull. Amer. Meteor. Soc.,80, 217–243.

  • Battisti, D. S., 1988: Dynamics and thermodynamics of a warming event in a coupled tropical atmosphere–ocean model. J. Atmos. Sci.,45, 2889–2919.

  • Behringer, D. W., M. Ji, and A. Leetmaa, 1998: An improved coupled model for ENSO prediction and implications for ocean initialization. Part I: The ocean data assimilation system. Mon. Wea. Rev.,126, 1013–1021.

  • Boulanger, J.-P., and C. Menkes, 1995: Propagation and reflection of long equatorial waves in the Pacific Ocean during 1992–1993 El Niño. J. Geophys. Res.,100, 25 041–25 059.

  • Brasseur, P., J. Ballabrera-Poy, and J. Verron, 1999: Assimilation of altimetric data in the mid-latitude oceans using the SEEK filter with an eddy-resolving, primitive equation model. J. Mar. Syst.,22, 269–294.

  • Busalacchi, A. J., M. J. McPhaden, and J. Picaut, 1994: Variability in equatorial Pacific sea surface topography during the verification phase of the TOPEX/POSEIDON Mission. J. Geophys. Res.,99, 25 725–25 738.

  • Cane, M. A., and R. J. Patton, 1984: A numerical model for low-frequency equatorial dynamics. J. Phys. Oceanogr.,14, 1853–1863.

  • ——, S. E. Zebiak, and S. C. Dolan, 1986: Experimental forecasts of El Niño. Nature,321, 827–832.

  • ——, A. Kaplan, R. N. Miller, B. Tang, E. C. Hackert, and A. J. Busalacchi, 1996: Mapping tropical Pacific sea level: Data assimilation via a reduced state space Kalman filter. J. Geophys. Res.,101, 22 599–22 617.

  • Chen, D., S. E. Zebiak, A. J. Busalacchi, and M. A. Cane, 1995: An improved procedure for El Niño forecasting: Implications for predictability. Science,269, 1699–1702.

  • ——, M. A. Cane, S. E. Zebiak, and A. Kaplan, 1998: The impact of sea level data assimilation on the Lamont model prediction of the 1997/98 El Niño Geophys. Res. Lett.,25, 2837–2840.

  • ——, ——, and ——, 1999: The impact of NSCAT winds on predicting the 1997/98 El Niño: A case study with the Lamont-Doherty Earth Observatory model. J. Geophys. Res.,104, 11 321–11 327.

  • Dewitte, B., and C. Perigaud, 1996: El Niño–La Niña events simulated with Cane and Zebiak’s model and observed with satellite or in situ data. Part II: Model forced with observations. J. Climate,9, 1188–1207.

  • Gill, A. E., 1980: Some simple solutions for heat-induced tropical circulations. Quart. J. Roy. Meteor. Soc.,106, 447–462.

  • Goswami, B. N., and J. Shukla, 1991: Predictability of a coupled ocean–atmosphere model. J. Climate,4, 3–22.

  • Graham, N. E., J. Michaelson, and T. P. Barnett, 1987: An investigation of the El Niño–Southern Oscillation cycle with statistical models, 1, Precursor field characteristics. J. Geophys. Res.,92, 14 251–14 270.

  • Ide, K., A. F. Bennett, P. Courtier, M. Ghil, and A. C. Lorenc, 1997:Unified notation for data assimilation: Operational, sequential and variational. J. Meteor. Soc. Japan,75, 181–189.

  • Kalman, R. E., and R. S. Bucy, 1960: A new approach to linear filtering and prediction problems. J. Basic Eng. (ASME),82D, 35–45.

  • Kaplan, A., Y. Kushnir, M. A. Cane, and M. B. Blumenthal, 1997: Reduced space optimal analysis for historical data sets: 136 years of Atlantic sea surface temperatures. J. Geophys. Res.,102, 27 835–27 860.

  • Kirtman, B. P., J. Shukla, B. Huang, Z. Zhu, and E. K. Schneider, 1997: Multiseasonal predictions with a coupled Tropical Ocean–Global Atmosphere system. Mon. Wea. Rev.,125, 789–808.

  • Latif, M., and M. Flügel, 1991: An investigation of short-range climate predictability in the tropical Pacific. J. Geophys. Res.,96, 2661–2673.

  • ——, A. Sterl, E. Maier-Reimer, and M. M. Junge, 1993: Climate variability in a coupled GCM. Part I: The tropical Pacific. J. Climate,6, 5–21.

  • ——, and Coauthors, 1998: A review of the predictability and prediction of ENSO. J. Geophys. Res.,103, 14 375–14 393.

  • Leon, S. J., 1990: Linear Algebra with Applications. Macmillan Publishing Company, 458 pp.

  • McPhaden, M. J., 1999: Genesis and evolution of the 1997–98 El Niño. Science,283, 950–954.

  • ——, and X. Yu, 1999: Equatorial waves and the 1997–98 El Niño. Geophys. Res. Lett.,26, 2961–2964.

  • ——, and Coauthors, 1998: The tropical ocean–global atmosphere observing system: A decade of progress. J. Geophys. Res.,103, 14 169–14 240.

  • Neelin, J. D., M. Latif, and F.-F. Jin, 1994: Dynamics of coupled ocean–atmosphere models: The tropical problem. Annu. Rev. Fluid Mech.,26, 617–659.

  • Overland, J. E., and R. W. Preisendorfer, 1982: A significance test for principal components applied to a cyclone climatology. Mon. Wea. Rev.,110, 1–4.

  • Pham, D. T., J. Verron, and M.-C. Roubaud, 1998: Singular evolutive extended Kalman filter with EOF initialization for data assimilation in oceanography. J. Mar. Syst.,16, 323–340.

  • Press, W. H., S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, 1992: Numerical Recipes in Fortran. Cambridge University Press, 963 pp.

  • Reynolds, R. W., and T. M. Smith, 1994: Improved global sea surface temperature analyses. J. Climate,7, 929–948.

  • Rosati, A., K. Miyakoda, and R. Gudgel, 1997: The impact of ocean initial conditions on ENSO forecasting with a coupled model. Mon. Wea. Rev.,125, 754–772.

  • Shriver, J. F., and J. J. O’Brien, 1995: Low-frequency variability of the equatorial Pacific using a new pseudostress dataset: 1930–1989. J. Climate,8, 2762–2768.

  • Smith, T. M., R. W. Reynolds, R. E. Livezey, and D. C. Stokes, 1996:Reconstruction of historical sea surface temperatures using Empirical Orthogonal Functions. J. Climate,9, 1403–1420.

  • Stricherz, J. N., J. J. O’Brien, and D. M. Legler, 1992: Atlas of the Florida State University tropical pacific winds for TOGA. Mesoscale Air–Sea Interaction Group Tech. Report, The Florida State University, Tallahassee, FL, 275 pp.

  • Trenberth, K. E., and T. J. Hoar, 1997: El Niño and climate change. Geophys. Res. Lett.,24, 3057–3060.

  • Verron, J., L. Gourdeau, D. T. Pham, R. Murtugudde, and A. J. Busalacchi, 1999: An extended Kalman filter to assimilate satellite altimeter data into a non-linear numerical model of the tropical Pacific: Method and validation. J. Geophys. Res.,104, 5441–5458.

  • Xu, J.-S., and H. von Storch, 1990: Principal oscillation pattern: Prediction of the state of ENSO. J. Climate,3, 1316–1329.

  • Xue, Y., M. A. Cane, and S. E. Zebiak, 1997: Predictability of a coupled model of ENSO using singular vector analysis. Part II:Optimal growth in seasonal background and ENSO cycles. Mon. Wea. Rev.,125, 2043–2056.

  • ——, A. Leetmaa, and M. Ji, 2000: ENSO prediction with Markov models: The impact of sea level. J. Climate,13, 849–871.

  • Zebiak, S. E., 1986: Atmospheric convergence feedback in a simple model for El Niño. Mon. Wea. Rev.,114, 1263–1271.

  • ——, 1989: Oceanic heat content variability and El Niño cycles. J. Phys. Oceanogr.,19, 475–486.

  • ——, and M. A. Cane, 1987: A model El Niño/Southern Oscillation. Mon. Wea. Rev.,115, 2262–2278.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 247 105 2
PDF Downloads 70 28 4