An ENSO stability analysis. Part I: results from a hybrid coupled model

In this study, we use the Bjerknes stability (BJ) index as a tool to investigate overall El Niño-Southern Oscillation (ENSO) stability in a hybrid-coupled model (HCM) with various atmosphere and ocean background states. This HCM shows that ENSO growth rates as measured by the BJ index and linear growth rates estimated directly with a time series of the Niño 3.4 indices from the coupled model simulations exhibit similar dependence on background states, coupling strength, and thermodynamic damping intensity. That is, the BJ index and linear growth rates increase with a decrease in the intensity of the background wind, an increase in coupling strength, or a decrease in the intensity of thermodynamic damping, although the BJ index tends to overestimate the growth rate. A detailed analysis of the components of the BJ index formula suggests the importance of model climatological background states and oceanic dynamic parameters in determining ENSO stability. We conclude that the BJ index may serve as a useful tool for qualitatively evaluating the overall ENSO stability in coupled models or in observations without a detailed eigen-analysis that is difficult to perform in models more complex than relatively simple models.     

ENSO in a hybrid coupled model. Part II: prediction with piggyback data assimilation

A hybrid coupled model (HCM) for the tropical Pacific ocean-atmosphere system is employed for ENSO prediction. The HCM consists of the Geophysical Fluid Dynamics Laboratory ocean general circulation model and an empirical atmospheric model. In hindcast experiments, a correlation skill competitive to other prediction models is obtained, so we use this system to examine the effects of several initialization schemes on ENSO prediction. Initialization with wind stress data and initialization with wind stress reconstructed from SST using the atmospheric model give comparable skill levels. In re-estimating the atmospheric model in order to prevent hindcast-period wind information from entering through empirical atmospheric model, we note some sensitivity to the estimation data set, but this is considered to have limited impact for ENSO prediction purposes. Examination of subsurface heat content anomalies in these cases and a case forced only by the difference between observed and reconstructed winds suggests that at the current level of prediction skill, the crucial wind components for initialization are those associated with the slow ENSO mode, rather than with atmospheric internal variability. A “piggyback” suboptimal data assimilation is tested in which the Climate Prediction Center data assimilation product from a related ocean model is used to correct the ocean initial thermal field. This yields improved skill, suggesting that not all ENSO prediction systems need to invest in costly data assimilation efforts, provided the prediction and assimilation models are sufficiently close. Received: 17 April 1998 / Accepted: 22 July 1999     

A Hybrid Coupled Model for the Pacific Ocean–Atmosphere System.Part I: Description and Basic Performance

A hybrid coupled model(HCM) is constructed for El Nino–Southern Oscillation(ENSO)-related modeling studies over almost the entire Pacific basin. An ocean general circulation model is coupled to a statistical atmospheric model for interannual wind stress anomalies to represent their dominant coupling with sea surface temperatures. In addition, various relevant forcing and feedback processes exist in the region and can affect ENSO in a significant way; their effects are simply represented using historical data and are incorporated into the HCM, including stochastic forcing of atmospheric winds, and feedbacks associated with freshwater flux, ocean biology-induced heating(OBH), and tropical instability waves(TIWs). In addition to its computational efficiency, the advantages of making use of such an HCM enable these related forcing and feedback processes to be represented individually or collectively, allowing their modulating effects on ENSO to be examined in a clean and clear way. In this paper, examples are given to illustrate the ability of the HCM to depict the mean ocean state, the circulation pathways connecting the subtropics and tropics in the western Pacific, and interannual variability associated with ENSO. As satellite data are taken to parameterize processes that are not explicitly represented in the HCM, this work also demonstrates an innovative method of using remotely sensed data for climate modeling. Further model applications related with ENSO modulations by extratropical influences and by various forcings and feedbacks will be presented in Part II of this study.     

Separating freshwater flux effects on ENSO in a hybrid coupled model of the tropical Pacific

Climate Dynamics - Freshwater flux (FWF) at the sea surface, defined as precipitation minus evaporation, is a major atmospheric forcing to the ocean that affects sea surface salinity (SSS) and...     

The intraseasonal oscillation in ECHAM4 Part I: coupled to a comprehensive ocean model

This study discusses the representation of the intraseasonal oscillation (ISO) in three simulations with the ECHAM4 atmosphere general circulation model (GCM). First, the model is forced by AMIP sea surface temperatures (SST), then coupled to the OPYC3 global ocean GCM and third forced by OPYC3 SSTs to clarify possible air-sea interactions and connections of the ISO and the ENSO cycle. The simulations are compared to ECMWF reanalysis data and NOAA outgoing longwave radiation (OLR) observations. Although previous studies have shown that the ECHAM4 GCM simulates an ISO-like oscillation, the main deficits are an overly fast eastward propagation and an eastward displacement of the main ISO activity, which is shown with a composite analysis of daily data between 1984 to 1988 for the reanalysis and the AMIP simulation, 25 years of the coupled integration, and a five year subset of the coupled SST output used for the OPYC3 forced atmosphere GCM experiment. These deficits are common to many atmospheric GCMs. The composites are obtained by principal oscillation pattern (POP). The POPs are also used to investigate the propagation speed and the interannual variability of the main ISO activity. The present coupled model version reveals no clear improvements in the ISO simulation compared to the uncoupled version forced with OPYC3 SSTs, although it is shown that the modeled ISO influences the simulated high-frequency SST variability in the coupled GCM. Within the current analysis, ECHAM4 forced by AMIP SSTs provides the most reasonable ISO simulation. However, it is shown that the maximum amplitudes of the annual cycle of the ISO variability in all analyzed model versions are reached too late in the year (spring and summer) compared to the observations (winter and spring). Additionally, the ENSO cycle influences the interannual variability of the ISO, which is revealed by 20 years of daily reanalysis data and 100 years of the coupled integration. The ENSO cycle is simulated by the coupled model, although there is a roughly 1 K cold bias in the East Pacific in the coupled model. This leads to a diminished influence of the ENSO cycle on the spatial variability of the modeled ISO activity compared to observations. This points out the strong sensitivity of the SST on the ISO activity. Small biases in the SST appear to cause large deterioration in the modeled ISO.     

Climate simulations with the global coupled atmosphere-ocean model ECHAM2/OPYC Part I: present-day climate and ENSO events

In this study the global coupled atmosphere-ocean general circulation model ECHAM2/OPYC and its performance in simulating the present-day climate is presented. The model consists of the T21-spectral atmosphere general circulation model ECHAM2 and the ocean general circulation model OPYC with a resolution corresponding to a T42 Gaussian grid, with increasing resolution towards the equator. The sea-ice is represented by a dynamic thermodynamic sea-ice model with rheology. Both models are coupled using the flux correction technique. With the coupled model ECHAM2/OPYC a 210-year integration under present-day greenhouse gas conditions has been performed. The coupled model simulates a realistic mean climate state, which is close to the observations. The model generates several ENSO events without external forcing. Using traditional and advanced (POP-technique) methods these ENSO events have been analyzed. The results are consistent with the delayed action oscillator theory. The model simulates both a tropical and an extra-tropical response to ENSO, which are in good agreement with observations.     

The retrospective prediction of ENSO from 1881 to 2000 by a hybrid coupled model: (II) Interdecadal and decadal variations in predictability

In this study, the retrospective predictions of ENSO (El Niño and Southern Oscillation) were performed for the period from 1881 to 2000 using a hybrid coupled model, which is an ocean general circulation model coupled to a linear statistical atmospheric model, and using a newly developed initialization scheme of SST assimilation by Ensemble Kalman Filter. With the retrospective predictions of the past 120 years, some important issues of ENSO predictability (measured by correlation and RMSE skills of NINO3 sea surface temperature anomaly index) were studied including decadal/interdecadal variations in ENSO predictability and the mechanisms responsible for these variations. Emphasis was placed on investigating the relationship between ENSO predictability and various characteristics of ENSO system such as the signal strength, the irregularity of periodicity, the noise and the nonlinearity. It is found that there are significant decadal/interdecadal variations in the prediction skills of ENSO during the past 120 years. The ENSO events were more predictable during the late nineteenth and the late twentieth centuries. The decadal/interdecadal variations of prediction skills are strongly related to the strength of sea-surface temperature anomaly (SSTA) signals, especially to the strength of SSTA signals at the frequencies of 2–4 year periods. The SSTA persistence, dominated by SSTA signals at frequencies over 4-year periods, also has a positive relationship to prediction skills. The high-frequency noise, on the other hand, has a strong inverse relationship to prediction skills, suggesting that it also probably plays an important role in ENSO predictability.     

Suppression of ENSO in a coupled model without water vapor feedback

 We examine 800-year time series of internally generated variability in both a coupled ocean-atmosphere model where water vapor anomalies are not allowed to interact with longwave radiation and one where they are. The ENSO-like phenomenon in the experiment without water vapor feedback is drastically suppressed both in amplitude and geographic extent relative to the experiment with water vapor feedback. Surprisingly, the reduced amplitude of ENSO-related sea surface temperature anomalies in the model without water vapor feedback cannot be attributed to greater longwave damping of sea surface temperature. (Differences between the two experiments in radiative feedback due to clouds counterbalance almost perfectly the differences in radiative feedback due to water vapor.) Rather, the interaction between water vapor anomalies and longwave radiation affects the ENSO-like phenomenon through its influence on the vertical structure of radiative heating: Because of the changes in water vapor associated with it, a given warm equatorial Pacific sea surface temperature anomaly is associated with a radiative heating profile that is much more gravitationally unstable when water vapor feedback is present. The warm sea surface temperature anomaly therefore results in more convection in the experiment with water vapor feedback. The increased convection, in turn, is related to a larger westerly wind-stress anomaly, which creates a larger decrease in upwelling of cold water, thereby enhancing the magnitude of the original warm sea surface temperature anomaly. In this manner, the interaction between water vapor anomalies and longwave radiation magnifies the air-sea interactions at the heart of the ENSO phenomenon; without this interaction, the coupling between sea surface temperature and wind stress is effectively reduced, resulting in smaller amplitude ENSO episodes with a more limited geographical extent. Received: 26 March 1999 / Accepted: 25 October 1999     

A standard test for AGCMs including their physical parametrizations: I: The proposal

To assist model intercomparison and development of a set of eight numerical experiments is proposed as a test-bed for the interaction of dynamics and physical parameterizations in atmospheric GCMs. The framework for the experiments is that of an aqua-planet and the prescribed sea-surface temperatures (SSTs) are highly idealized.     

Annual ENSO simulated in a coupled ocean–atmosphere model

Using an output from 200-year integration of the Scale Interaction Experiment of EU project-F1 model (SINTEX-F1), the annual ENSO reproduced in the coupled general circulation model is investigated, suggesting the importance of reproducing an annual cycle in realistically simulating ENSO events. Although many features of the annual ENSO are reproduced, the northward expansion of sea surface temperature anomaly (SSTA) in the eastern tropical Pacific stays south of the equator. It is suggested that this model bias is due to the excitation of the too strong Rossby waves in the southeastern tropical Pacific, which reflect at the western boundary and intrude into the eastern equatorial Pacific. The zonal wind stress anomaly along the equator also plays an important role in generating the equatorial Kelvin waves. The amplitude of SSTA for the annual ENSO mode is reproduced, but its variance is only 20% of the observation; this is again due to the lack of northward migration of seasonal SSTA in the equatorial region and weaker coastal Kelvin waves along South America. Remedies for the model bias are discussed.     

ENSO phase-locking in an ocean-atmosphere coupled model FGCM-1.0

The mean climatology and the basic characteristics of the ENSO cycle simulated by a coupled model FGCM-1.0 are investigated in this study. Although with some common model biases as in other directly coupled models, FGCM-1.0 is capable of producing the interannual variability of the tropical Pacific, such as the ENSO phenomenon. The mechanism of the ENSO events in the coupled model can be explained by “delayed oscillator” and “recharge-discharge” hypotheses. Compared to the observations, the simulated ENSO events show larger amplitude with two distinctive types of phase-locking： one with its peak phase-locked to boreal winter and the other to boreal summer. These two types of events have a similar frequency of occurrence, but since the second type of event is seldom observed, it may be related to the biases of the coupled model. Analysis show that the heat content anomalies originate from the central south Pacific in the type of events peaking in boreal summer, which can be attributed to a different background climatology from the normal events. The mechanisms of their evolutions are also discussed.     

Parameter sensitivity of a coupled atmosphere-ocean model

 The sensitivity of a coupled model to the oceanic vertical diffusion coefficient κ _v is examined. This is compared to the sensitivity of an ocean-only model forced by mixed boundary conditions (BC). The atmospheric component of the coupled model is a moist energy balance model. The ocean component is a 12-level geostrophic model, defined on a midlatitude β-plane. Atmosphere and ocean are coupled through the fluxes of heat and moisture at their interface. The coupled model contains a number of feedback processes which are not represented in the ocean-only model. This results in a temperature and salinity response to κ _v which is stronger in the coupled model than in the ocean-only model. On the other hand, there is a weaker response in oceanic processes such as meridional heat transport, deep-water formation at high latitudes, etc. Ocean-only sensitivity experiments were also performed with modified BCs, which parametrise the feedback processes included in the coupled model. These are the modified thermal BC of Rahmstorf and Willebrand and a modified freshwater BC proposed in the present study. Large-scale features of the response in oceanic surface fields are well represented with modified BCs. However, the sensitivity of the deep ocean temperature is only partly captured due to local differences in the surface response. The scaling behavior of the zonal overturning stream function was found to depend on the surface BCs. In contrast to this, the meridional overturning stream function basically scales with κ^0.5 _v in all sensitivity experiments. Differences in the heat transport response among the experiments are thus primarily related to differences in the temperature response. Received: 28 February 1997/Accepted: 12 September 1997     

Impact of intra-daily SST variability on ENSO characteristics in a coupled model

This paper explores the impact of intra-daily Sea Surface Temperature (SST) variability on the tropical large-scale climate variability and differentiates it from the response of the system to the forcing of the solar diurnal cycle. Our methodology is based on a set of numerical experiments based on a fully global coupled ocean–atmosphere general circulation in which we alter (1) the frequency at which the atmosphere sees the SST variations and (2) the amplitude of the SST diurnal cycle. Our results highlight the complexity of the scale interactions existing between the intra-daily and inter-annual variability of the tropical climate system. Neglecting the SST intra-daily variability results, in our CGCM, to a systematic decrease of 15% of El Ni?o—Southern Oscillation (ENSO) amplitude. Furthermore, ENSO frequency and skewness are also significantly modified and are in better agreement with observations when SST intra-daily variability is directly taken into account in the coupling interface of our CGCM. These significant modifications of the SST interannual variability are not associated with any remarkable changes in the mean state or the seasonal variability. They can therefore not be explained by a rectification of the mean state as usually advocated in recent studies focusing on the diurnal cycle and its impact. Furthermore, we demonstrate that SST high frequency coupling is systematically associated with a strengthening of the air-sea feedbacks involved in ENSO physics: SST/sea level pressure (or Bjerknes) feedback, zonal wind/heat content (or Wyrtki) feedback, but also negative surface heat flux feedbacks. In our model, nearly all these results (excepted for SST skewness) are independent of the amplitude of the SST diurnal cycle suggesting that the systematic deterioration of the air-sea coupling by a daily exchange of SST information is cascading toward the major mode of tropical variability, i.e. ENSO.     

ENSO cycle in a coupled ocean-atmosphere model and its negative feedback mechanism

Summary A coupled ocean-atmosphere anomaly model has been developed for simulating ENSO cycle and its mechanism-study in this paper. After a long model run, the coupled model is successful in demonstrating ENSO-like irregular interannual variability and corresponding horizontal spatial structures. Based on the simulated results, the dynamics and the thermodynamics of the model ENSO cycle have been investigated, and in particular the negative feedback mechanisms that act to oppose instability of air-sea interaction, inducing termination of warm and cold events, have been examined. A detailed analysis of the oceanic wave dynamical properties and heat budget of the SST changes in a representative cycle suggest that the negative feedback mechanism to check the unstable growth of a warm event obviously differs from that of a cold event. The mechanism that induces decay and termination of a cold event is closely related to the negative, delayed feedback effect produced by the oceanic dynamical wave reflection at the western boundary. However, independent of the wave reflection effect, the negative feedback mechanism by which the coupled system returns from a warm event is associated with a slowly eastward-propagating coupling mode. Accompanied with the strong unstable development of the equatorial positive SST anomaly, the anomalous upwelling of cold water generated off the equator and the nonlinear anomalous meridional advection generated in the equator west of instability area jointly restrain the instability and finally plunge the system from a mature warm phase into a weak cold phase. A comparison between the results from the present model and the previous works is also discussed in this paper.With 16 Figures     

The impact of new physical parametrizations in the Hadley Centre climate model: HadAM3

 Results are presented from the latest version of the Hadley Centre climate model, HadAM3 (Hadley Centre Atmospheric Model version 3). It represents a significant improvement over the previous version, HadAM2b. This is demonstrated using a series of ten year integrations with AMIP (Atmospheric Model Intercomparison Project) boundary conditions. The work covers three aspects of model performance: (1) it shows the improvements in the mean climate in changing from HadAM2b to HadAM3; (2) it demonstrates that the model now compares well with observations and (3) it isolates the impacts of new physical parametrizations. Received: 17 August 1998 / Accepted: 20 July 1999     

CLIMBER-2: a climate system model of intermediate complexity. Part II: model sensitivity

 A set of sensitivity experiments with the climate system model of intermediate complexity CLIMBER-2 was performed to compare its sensitivity to changes in different types of forcings and boundary conditions with the results of comprehensive models (GCMs). We investigated the climate system response to changes in freshwater flux into the Northern Atlantic, CO₂ concentration, solar insolation, and vegetation cover in the boreal zone and in the tropics. All these experiments were compared with the results of corresponding experiments performed with different GCMs. Qualitative, and in many respects, quantitative agreement between the results of CLIMBER-2 and GCMs demonstrate the ability of our climate system model of intermediate complexity to address diverse aspects of the climate change problem. In addition, we used our model for a series of experiments to assess the impact of some climate feedbacks and uncertainties in model parameters on the model sensitivity to different forcings. We studied the role of freshwater feedback and vertical ocean diffusivity for the stability properties of the thermohaline ocean circulation. We show that freshwater feedback plays a minor role, while changes of vertical diffusivity in the ocean considerably affect the circulation stability. In global warming experiments we analysed the impact of hydrological sensitivity and vertical diffusivity on the long-term evolution of the thermohaline circulation. In the boreal and tropical deforestation experiments we assessed the role of an interactive ocean and showed that for both types of deforestation scenarios, an interactive ocean leads to an additional cooling due to albedo and water vapour feedbacks. Received: 28 May 2000 / Accepted: 9 November 2000     

Climate sensitivity to changes in solar insolation in a simple coupled climate model

A simple coupled ocean, atmosphere and sea-ice model is presented. The idealised model consists of a zonally averaged land and ocean strip of constant angular width extending from pole to pole. The meridional energy transport in the ocean is modelled by contributions from the large scale thermohaline overturning cells and from horizontal diffusive fluxes. The atmospheric meridional energy transports are parametrised as diffusive fluxes in addition to advective transports in the Hadley domain. This parametrisation resolves the equatorward moisture transport as well as the poleward transport of potential energy in the upper branch of the Hadley circulation. The model reproduces the annual averaged meridional energy transports in the climate system with a small number of free model parameters. The basic feedbacks between the three climatic components are studied by investigating the model's sensitivity towards reductions in the solar insolation. It is found that the meridional energy transport in the ocean does not amplify the ice albedo feedback. This has important implications for modelling the climate sensitivity in atmosphere-only models, as these would exaggerate the sensitivity to changes in the solar insolation if their parametrisations of the meridional energy transport are constrained by surface temperatures. The role of the dependence of the atmospheric transports on the meridional temperature gradient is shown to have a significant influence on the sensitivity on the coupled model, and the inclusion of seasonal cycles greatly increase the models sensitivity. The Hadley circulation does significantly alter the strength of the ice-albedo feedback in the coupled model. The idealised configuration of the model makes it a useful tool for studying the feedbacks in the ocean-atmosphere-sea ice system in the context of the "Snowball Earth" hypothesis.     

Impacts of shortwave radiation forcing on ENSO: a study with a coupled tropical ocean-atmosphere model

 We describe a coupled tropical ocean-atmosphere model that represents a new class of models that fill the gap between anomaly coupled models and fully coupled general circulation models. Both the atmosphere and ocean are described by two and half layer primitive equation models, which emphasize the physical processes in the oceanic mixed layer and atmospheric boundary layer. Ocean and atmosphere are coupled through both momentum and heat flux exchanges without explicit flux correction. The coupled model, driven by solar radiation, reproduces a realistic annual cycle and El Nino-Southern Oscillation (ENSO). In the presence of annual mean shortwave radiation forcing, the model exhibits an intrinsic mode of ENSO. The oscillation period depends on the mean forcing that determines the coupled mean state. A perpetual April (October) mean forcing prolongs (shortens) the oscillation period through weakening (enhancing) the mean upwelling and mean vertical temperature gradients. The annual cycle of the solar forcing is shown to have fundamental impacts on the behavior of ENSO cycles through establishing a coupled annual cycle that interacts with the ENSO mode. Due to the annual cycle solar forcing, the single spectral peak of the intrinsic ENSO mode becomes a double peak with a quasi-biennial and a low-frequency (4–5 years) component; the evolution of ENSO becomes phase-locked to the annual cycle; and the amplitude and frequency of ENSO become variable on an interdecadal time scale due to interactions of the mean state and the two ENSO components. The western Pacific monsoon (the annual shortwave radiation forcing in the western Pacific) is primarily responsible for the generation of the two ENSO components. The annual march of the eastern Pacific ITCZ tends to lock ENSO phases to the annual cycle. The model's deficiencies, limitations, and future work are also discussed. Received: 15 June 1999 / Accepted: 11 December 1999     

The variation of ENSO characteristics associated with atmospheric parameter perturbations in a coupled model

We analyse the differences in the properties of the El Niño Southern Oscillation (ENSO) in a set of 17 coupled integrations with the flux-adjusted, 19-level HadCM3 model with perturbed atmospheric parameters. Within this ensemble, the standard deviation of the NINO3.4 deseasonalised SSTs ranges from 0.6 to 1.3 K. The systematic changes in the properties of the ENSO with increasing amplitude confirm that ENSO in HadCM3 is prevalently a surface (or SST) mode. The tropical-Pacific SST variability in the ensemble of coupled integrations correlates positively with the SST variability in the corresponding ensemble of atmosphere models coupled with a static mixed-layer ocean (“slab” models) perturbed with the same changes in atmospheric parameters. Comparison with the respective coupled ENSO-neutral climatologies and with the slab-model climatologies indicates low-cloud cover to be an important controlling factor of the strength of the ENSO within the ensemble. Our analysis suggests that, in the HadCM3 model, increased SST variability localised in the south-east tropical Pacific, not originating from ENSO and associated with increased amounts of tropical stratocumulus cloud, causes increased ENSO variability via an atmospheric bridge mechanism. The relationship with cloud cover also results in a negative correlation between the ENSO activity and the model’s climate sensitivity to doubling CO₂.