Tracking Down the ENSO Delayed Oscillator with an Adjoint OGCM

Geert Jan van Oldenborgh KNMI, De Bilt, the Netherlands

Search for other papers by Geert Jan van Oldenborgh in
Current site
Google Scholar
PubMed
Close
,
Gerrit Burgers KNMI, De Bilt, the Netherlands

Search for other papers by Gerrit Burgers in
Current site
Google Scholar
PubMed
Close
,
Stephan Venzke Max-Planck Institut für Meteorologie, Hamburg, Germany

Search for other papers by Stephan Venzke in
Current site
Google Scholar
PubMed
Close
,
Christian Eckert Max-Planck Institut für Meteorologie, Hamburg, Germany

Search for other papers by Christian Eckert in
Current site
Google Scholar
PubMed
Close
, and
Ralf Giering Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts

Search for other papers by Ralf Giering in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

According to the delayed-oscillator picture of ENSO, a positive SST anomaly in the eastern tropical Pacific will cause westerly wind anomalies closer to the date line to first give a positive feedback, and later, via planetary wave reflection at the western boundary, a negative feedback. The aim of this study is to follow a chain of sensitivities that lead to a delayed-oscillator mechanism in a general circulation model. To this end, the adjoint of such an ocean model is used for studying sensitivities of ENSO indices.

The ocean model used in this study is the Hamburg Ocean Primitive Equation (HOPE) ocean general circulation model. Its adjoint has been constructed using the Adjoint Model Compiler. Applied to a scalar function computed with a forward model run, an adjoint run goes back in time and calculates sensitivities as the derivatives of this function to forcing fields or ocean state variables at earlier times.

Results from six adjoint runs are reported, tracing the sensitivities of the NINO3 and NINO3.4 indices in October 1987, December 1987, and December 1988, as simulated by a Pacfic-only version of HOPE forced by ECHAM-3 fluxes.

The sensitivities to sea level can be followed back in time for more than a year. They are nonlocal: patterns propagate back in time that are identified as adjoint Kelvin and n = 1, 2, and 3 Rossby waves, with speeds compatible with those obtained from model density profiles. Both the first and the second baroclinic modes seem to play a role. In contrast, the model sensitivities to heat flux, zonal surface currents, and SST are local and decay in about a month.

The sensitivities to the wind stress agree with the wave interpretation of the sea-level sensitivities, but only the n = 1 Rossby wave is visible. Going back in time, the sensitivity to westerly anomalies along the equator changes sign, in agreement with the delayed-oscillator picture.

Finally, a statistical atmosphere model is used to convert sensitivities to wind stress at a given time to sensitivities to SST through the atmosphere at that time. Focusing on the sensitivities to the ENSO index region itself at an earlier time then closes the circle. These sensitivities have a natural interpretation as delayed-oscillator coefficients and show the expected behavior of a positive sensitivity in the recent past changing to a negative sensitivity at longer lags. However, the strength of these feedbacks, and hence the relevance of this mechanism in ENSO simulated in HOPE, cannot be determined accurately.

Corresponding author address: Dr. Geert Jan van Oldenborgh, KNMI, KS/OO, P.O. Box 201, NL-3730 AE De Bilt, the Netherlands.

Email: oldenbor@knmi.nl

Abstract

According to the delayed-oscillator picture of ENSO, a positive SST anomaly in the eastern tropical Pacific will cause westerly wind anomalies closer to the date line to first give a positive feedback, and later, via planetary wave reflection at the western boundary, a negative feedback. The aim of this study is to follow a chain of sensitivities that lead to a delayed-oscillator mechanism in a general circulation model. To this end, the adjoint of such an ocean model is used for studying sensitivities of ENSO indices.

The ocean model used in this study is the Hamburg Ocean Primitive Equation (HOPE) ocean general circulation model. Its adjoint has been constructed using the Adjoint Model Compiler. Applied to a scalar function computed with a forward model run, an adjoint run goes back in time and calculates sensitivities as the derivatives of this function to forcing fields or ocean state variables at earlier times.

Results from six adjoint runs are reported, tracing the sensitivities of the NINO3 and NINO3.4 indices in October 1987, December 1987, and December 1988, as simulated by a Pacfic-only version of HOPE forced by ECHAM-3 fluxes.

The sensitivities to sea level can be followed back in time for more than a year. They are nonlocal: patterns propagate back in time that are identified as adjoint Kelvin and n = 1, 2, and 3 Rossby waves, with speeds compatible with those obtained from model density profiles. Both the first and the second baroclinic modes seem to play a role. In contrast, the model sensitivities to heat flux, zonal surface currents, and SST are local and decay in about a month.

The sensitivities to the wind stress agree with the wave interpretation of the sea-level sensitivities, but only the n = 1 Rossby wave is visible. Going back in time, the sensitivity to westerly anomalies along the equator changes sign, in agreement with the delayed-oscillator picture.

Finally, a statistical atmosphere model is used to convert sensitivities to wind stress at a given time to sensitivities to SST through the atmosphere at that time. Focusing on the sensitivities to the ENSO index region itself at an earlier time then closes the circle. These sensitivities have a natural interpretation as delayed-oscillator coefficients and show the expected behavior of a positive sensitivity in the recent past changing to a negative sensitivity at longer lags. However, the strength of these feedbacks, and hence the relevance of this mechanism in ENSO simulated in HOPE, cannot be determined accurately.

Corresponding author address: Dr. Geert Jan van Oldenborgh, KNMI, KS/OO, P.O. Box 201, NL-3730 AE De Bilt, the Netherlands.

Email: oldenbor@knmi.nl

Save
  • Arakawa, A., and V. R. Lamb, 1977: Computational design of the basic dynamical process of the UCLA general circulation model. Methods Comput. Phys.,17, 173–265.

  • Barnett, T. P., M. Latif, N. 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., and C. F. Ropelewski, 1992: Prediction of ENSO episodes using canonical correlation analysis. J. Climate,5, 1316–1345.

  • Battisti, D. S., and A. C. Hirst, 1989: Interannual variability in a tropical atmosphere–ocean model: Influence of the basic state, ocean geometry and nonlinearity. J. Atmos. Sci.,46, 1687–1712.

  • Bjerknes, J., 1966: A possible response of the atmospheric Hadley circulation to equatorial anomalies of ocean temperature. Tellus,18, 820–829.

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

  • ——, and L.-L. Fu, 1996: Evidence of boundary reflection of Kelvin and first-mode Rossby waves from TOPEX/Poseidon sea level data. J. Geophys. Res.,101, 16 361–16 371.

  • Burgers, G., 1999: The El Niño stochastic oscillator. Climate Dyn., in press.

  • Chao, Y., and S. G. H. Philander, 1993: On the structure of the Southern Oscillation. J. Climate,6, 450–469.

  • Chelton, D. B., and M. G. Schlax, 1996: Global observations of oceanic Rossby waves. Science,272, 234–238.

  • DKRZ, 1992: The ECHAM 3 atmospheric general circulation model. Tech. Rep. 6, Deutsches Klimarechenzentrum, 189 pp. [Available from Deutsches Klimarechenzentrum, Bundesstr. 55, 20146 Hamburg, Germany.].

  • Fischer, M., and M. Latif, 1995: Assimilation of temperature and sea level observations into a primitive equation model of the tropical Pacific. J. Mar. Syst.,6, 31–46.

  • Frey, H., M. Latif, and T. Stockdale, 1997: The coupled GCM ECHO-2. Part I: The tropical Pacific. Mon. Wea. Rev.,125, 703–720.

  • Giering, R., 1994: Adjoint Model Compiler. Users manual, 34 pp. [Available online at http://klima47.dkrz.de/giering/amc.].

  • ——, 1997: Tangent linear and Adjoint Model Compiler, Users manual 1.2. 53 pp. [Available online at http://klima47.dkrz.de/giering/amc.].

  • ——, and T. Kaminski, 1998: Recipes for adjoint code construction. Trans. Math. Software,4, 437–474. [Available from R. Giering, Dept. of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, 77 Massachusetts Ave., Cambridge, MA 02139.].

  • Giese, B. S., and D. E. Harrison, 1990: Aspects of the Kelvin wave response to episodic wind forcing. J. Geophys. Res.,95, 7289–7312.

  • Gill, A. E., 1982: Atmosphere–Ocean Dynamics. Academic Press, 662 pp.

  • Hildebrandsson, H. H., 1897: Quelques recherches sur les centres d’action de l’atmosphere I. Sven. Vet. Akad. Handlingar,29(3), 3–8.

  • Ji, M., A. Leetmaa, and V. E. Kousky, 1996: Coupled model predictions of ENSO during the 1980s and 1990s at the National Centers for Environmental Prediction. J. Climate,9, 3105–3120.

  • Jin, F.-F., 1997: An equatorial recharge paradigm for ENSO. Part I: Conceptual model. J. Atmos. Sci.,54, 811–829.

  • Kessler, W. S., and M. J. McPhaden, 1995: Oceanic equatorial waves and the 1991–93 El Niño. J. Climate,8, 1757–1774.

  • ——, ——, and K. M. Weickman, 1995: Forcing of intraseasonal Kelvin waves in the equatorial Pacific. J. Geophys. Res.,100, 10 613–10 631.

  • Kleeman, R., A. Moore, and N. R. Smith, 1995: Assimilation of subsurface thermal data into a simple ocean model for the initialization of an intermediate tropical coupled ocean–atmosphere forecast model. Mon. Wea. Rev.,123, 3103–3114.

  • Latif, M., T. Stockdale, J. O. Wolff, G. Burgers, E. Maier-Reimer, M. M. Junge, K. Arpe, and L. Bengtsson, 1994: Climatology and variability in the ECHO coupled GCM. Tellus,46A, 351–366.

  • Levitus, S., 1982: Climatological Atlas of the World Ocean. NOAA Prof. Paper 13, 173 pp.

  • McCreary, J. P., and D. L. T. Anderson, 1991: An overview of coupled ocean–atmosphere models of El Niño and the Southern Oscillation. J. Geophys. Res.,96, 3125–3150.

  • Moore, A. M., and R. Kleeman, 1997a: The singular vectors of a coupled ocean–atmosphere model of ENSO. Part I: Thermodynamics, energetics and error growth. Quart. J. Roy. Meteor. Soc.,123, 953–981.

  • ——, and ——, 1997b: The singular vectors of a coupled ocean–atmosphere model of ENSO. Part II: Sensitivity studies and dynamical interpretation. Quart. J. Roy. Meteor. Soc.,123, 983–1006.

  • Neelin, J. D., 1991: The slow sea-surface temperature mode and the fast-wave limit: Analytic theory for tropical interannual oscillations and experiments in a hybrid coupled model. J. Atmos. Sci.,48, 584–606.

  • Parker, D. E., C. K. Folland, A. Bevan, M. N. Ward, M. Jackson, and K. Maskell, 1994: Marine surface data for analysis of climatic fluctuations of interannual to century timescales. Natural Climate Variability on Decade-to-Century Time Scales, D. G. Martinson et al., Eds., National Academy Press, 241–250.

  • Paulson, C. A., and J. J. Simpson, 1977: Irradiance measurements in the upper ocean. J. Phys. Oceanogr.,7, 952–956.

  • Philander, S. G., 1990: El Niño, La Niña, and the Southern Oscillation. Academic Press, 293 pp.

  • Picaut, C., and L. Sombarddier, 1993: Influence of density stratification and bottom depth on vertical mode structure functions in the tropical Pacific. J. Geophys. Res.,98, 14 727–14 737.

  • ——, and T. Delcroix, 1995: Equatorial wave sequence associated with warm pool displacements during the 1986–1989 El Niño–La Niña. J. Geophys. Res.,100, 18 393–18 408.

  • ——, M. Ioulalen, C. Menkes, and T. Delcroix, 1996: Mechanism of the zonal displacement of the Pacific warm pool: Implications for ENSO. Science,274, 1486–1489.

  • Roeckner, E., and Coauthors, 1992: Simulation of the present day climate with the ECHAM model: Impact of model physics and resolution. Rep. 93, 172 pp. [Available from Max-Planck Institut für Meteorologie, Bundesstr. 55, 20146 Hamburg, Germany.].

  • Schneider, N., T. P. Barnett, M. Latif, and T. Stockdale, 1996: Warm pool physics in a coupled GCM. J. Climate,9, 219–239.

  • Schopf, P. S., and M. J. Suarez, 1988: Vacillations in a coupled ocean–atmosphere model. J. Atmos. Sci.,45, 549–566.

  • Sirkes, Z., and E. Tziperman, 1997: Finite difference of adjoint or adjoint of finite difference? Mon. Wea. Rev.,125, 3373–3378.

  • Suarez, M. J., and P. S. Schopf, 1988: A delayed action oscillator for ENSO. J. Atmos. Sci.,45, 3283–3287.

  • Walker, G. T., 1923: Correlation in seasonal variations of weather VIII. Mem. Indian Meteor. Dep.,24, 75–131.

  • ——, 1924: Correlation in seasonal variations of weather IX. Mem. Indian Meteor. Dep.,24, 275–332.

  • Wolff, J.-O., E. Maier-Reimer, and S. Legutke, 1997: The Hamburg Ocean Primitive Equation model HOPE. Tech. Rep. 13, 97 pp. [Available from Deutsches Klimarechenzentrum, Bundesstr. 55, D-20146 Hamburg, Germany.].

  • Wyrtki, K., 1975: El Niño—The dynamic response of the equatorial Pacific Ocean to atmospheric forcing. J. Phys. Oceanogr.,5, 572–584.

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

  • ——, ——, ——, and T. N. Palmer, 1997b: Predictability of a coupled model of ENSO using singular vector analysis. Part II: Optimal growth and forecast skill. Mon. Wea. Rev.,125, 2057–2073.

  • Zebiak, S. E., and M. A. Cane, 1987: A model of 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 122 62 6
PDF Downloads 59 24 1