Impact of the equatorial Atlantic sea surface temperature on the tropical Pacific in a CGCM

Many coupled general circulation models (CGCMs) suffer from serious model bias in the zonal gradient of sea surface temperature (SST) in the equatorial Atlantic. The bias of the equatorial Atlantic SST (EASST) may affect the interannual variability of the equatorial Atlantic, which in turn may influence the state of the tropical Pacific. In this paper we investigate the impact of the bias and the interannual variability of the EASST on the tropical Pacific in a CGCM. To determine the impact of the interannual variability of the EASST on the tropical Pacific, we compare a run in a fully coupled mode (CTL run) and a run in which the EASST is nudged toward the climatological monthly mean of the SST in the CTL run, but full air-sea coupling is allowed elsewhere (AT_m run). We find that, when the interannual variability of the EASST is excluded, the thermocline depth in the eastern equatorial Pacific is deepened, and the amplitude of the El Niño/Southern Oscillation is reduced by 30 % compared to the CTL run. The impact of the bias of the EASST on the tropical Pacific is investigated by comparing the AT_m run and a run in which the EASST is nudged toward the observed climatological monthly mean SST (AT_o run). It is found that, when the bias of the EASST is removed (i.e. AT_o run), the Gill–Matsuno type response to the warm SST anomalies in the western equatorial Atlantic induces low-level cyclonic anomalies in the eastern South Pacific, which leads to a deeper thermocline and colder SST in the South Pacific as compared to AT_m. The colder SST in the South Pacific reduces the precipitation along the South Pacific convergence zone. Our results of the model experiments demonstrate the importance of the EASST to the tropical Pacific climate.     

Interdecadal variations of ENSO around 1999/2000

This paper discusses the interdecadal changes of the climate in the tropical Pacific with a focus on the corresponding changes in the characteristics of the El Niño–Southern Oscillation (ENSO). Compared with 1979–1999, the whole tropical Pacific climate system, including both the ocean and atmosphere, shifted to a lower variability regime after 1999/2000. Meanwhile, the frequency of ENSO became less regular and was closer to a white noise process. The lead time of the equatorial Pacific's subsurface ocean heat content in preceding ENSO decreased remarkably, in addition to a reduction in the maximum correlation between them. The weakening of the correlation and the shortening of the lead time pose more challenges for ENSO prediction, and is the likely reason behind the decrease in skill with respect to ENSO prediction after 2000. Coincident with the changes in tropical Pacific climate variability, the mean states of the atmospheric and oceanic components also experienced physically coherent changes. The warm anomaly of SST in the western Pacific and cold anomaly in the eastern Pacific resulted in an increased zonal SST gradient, linked to an enhancement in surface wind stress and strengthening of the Walker circulation, as well as an increase in the slope of the thermocline. These changes were consistent with an increase (a decrease) in precipitation and an enhancement (a suppression) of the deep convection in the western (eastern) equatorial Pacific. Possible connections between the mean state and ENSO variability and frequency changes in the tropical Pacific are also discussed.     

Changes in tropical Atlantic interannual variability from a substantial weakening of the meridional overturning circulation

In response to a substantial weakening of the Atlantic Meridional Overturning Circulation (AMOC)—from a coupled ocean–atmosphere general circulation model experiment—significant changes in the interannual variability are found over the tropical Atlantic, characterized by an increase of variance (by ~150 %) in boreal late spring-early summer and a decrease of variance (by ~60 %) in boreal autumn. This study focuses on understanding physical mechanisms responsible for these changes in interannual variability in the tropical Atlantic. It demonstrates that the increase of variability in spring is a consequence of an increase in the variance of the El Niño-Southern Oscillation, which has a large impact on the tropical Atlantic via anomalous surface heat fluxes. Winter El Niño (La Niña) affects the eastern equatorial Atlantic by decreasing (increasing) cloud cover and surface wind speed which is associated with anomalous downward (upward) short wave radiation and reduced (enhanced) upward latent heat fluxes, creating anomalous positive (negative) sea surface temperature (SST) anomalies over the region from winter to spring. On the other hand, the decrease of SST variance in autumn is due to a deeper mean thermocline which weakens the impact of the thermocline movement on SST variation. The comparison between the model results and observations is not straightforward owing to the influence of model biases and the lack of a major MOC weakening event in the instrumental record. However, it is argued that the basic physical mechanisms found in the model simulations are likely to be robust and therefore have relevance to understanding tropical Atlantic variability in the real world, perhaps with modified seasonality.     

Eastward propagating MJO during boreal summer and Indian monsoon droughts

Improved understanding of underlying mechanism responsible for Indian summer monsoon (ISM) droughts is important due to their profound socio-economic impact over the region. While some droughts are associated with ‘external forcing’ such as the El-Niño and Southern Oscillation (ENSO), many ISM droughts are not related to any known ‘external forcing’. Here, we unravel a fundamental dynamic process responsible for droughts arising not only from external forcing but also those associated with internal dynamics. We show that most ISM droughts are associated with at least one very long break (VLB; breaks with duration of more than 10 days) and that the processes responsible for VLBs may also be the mechanism responsible for ISM droughts. Our analysis also reveals that all extended monsoon breaks (whether co-occurred with El-Niño or not) are associated with an eastward propagating Madden–Julian Oscillation (MJO) in the equatorial Indian Ocean and western Pacific extending to the dateline and westward propagating Rossby waves between 10° and 25°N. The divergent Rossby wave associated with the dry phase of equatorial convection propagates westward towards Indian land, couple with the northward propagating dry phase and leads to the sustenance of breaks. Thus, the propensity of eastward propagating MJO during boreal summer is largely the cause of monsoon droughts. While short breaks are not accompanied by westerly wind events (WWE) over equatorial western Pacific favorable for initiating air–sea interaction, all VLBs are accompanied by sustained WWE. The WWEs associated with all VLB during 1975–2005 initiate air–sea interaction on intraseasonal time scale, extend the warm pool eastward allowing the convectively coupled MJO to propagate further eastward and thereby sustaining the divergent circulation over India and the monsoon break. The ocean–atmosphere coupling on interannual time scale (such as El-Niño) can also produce VLB, but not necessary.     

On the Atlantic–Pacific Niños connection: a multidecadal modulated mode

Atlantic and Pacific El Niño are the leading tropical oceanic variability phenomena at interannual timescales. Recent studies have demonstrated how the Atlantic Niño is able to influence on the dynamical processes triggering the development of the Pacific La Niña and vice versa. However, the stationarity of this interbasin connection is still controversial. Here we show for the first time that the Atlantic–Pacific Niños connection takes place at particular decades, coinciding with negative phases of the Atlantic Multidecadal Oscillation (AMO). During these decades, the Atlantic–Pacific connection appears as the leading coupled covariability mode between Tropical Atlantic and Pacific interannual variability. The mode is defined by a predictor field, the summer Atlantic Sea Surface Temperature (SST), and a set of predictand fields which represent a chain of atmospheric and oceanic mechanisms to generate the Pacific El Niño phenomenon: alteration of the Walker circulation, surface winds in western Pacific, oceanic Kelvin wave propagating eastward and impacting on the eastern thermocline and changes in the Pacific SST by internal Bjerknes feedback. We suggest that the multidecadal component of the Atlantic acts as a switch for El Niño prediction during certain decades, putting forward the AMO as the modulator, acting through changes in the equatorial Atlantic convection and the equatorial Pacific SST variability. These results could have a major relevance for the decadal prediction systems.     

Change of El Niño and La Niña amplitude asymmetry around 1980

Amplitude of El Niño and La Niña was significantly different during 1980–2016 but almost same during 1958–1979. The cause of this interdecadal change is investigated through an oceanic mixed-layer heat budget analysis. It was found that this interdecadal change was primarily attributed to the distinctive effects of nonlinear zonal temperature advection between the two periods. During 1980–2016 nonlinear zonal advection, working together with nonlinear meridional advection, contributes to the El Niño and La Niña amplitude asymmetry. During 1958–1979 the nonlinear zonal advection had an opposite effect. The difference in the nonlinear zonal advection between the two interdecadal periods was caused by distinctive longitudinal locations of El Niño centers. Maximum SST anomaly (SSTA) centers were confined near the coast of South America (east of 90° W) during the first period but appear near 110° W during the second period. Because of this difference, an anomalous eastward ocean surface current (caused by a positive thermocline depth anomaly during El Niño) would generate a negative (positive) nonlinear zonal advection before (after) 1980. The distinctive longitudinal locations of El Niño centers are possibly caused by the interdecadal changes of mean thermocline and high-frequency wind variability over the equatorial western-central Pacific. A hypothesis was put forth to understand distinctive initiation locations between El Niño and La Niña.     

Linear trends in sea surface temperature of the tropical Pacific Ocean and implications for the El Niño-Southern Oscillation

A principal component decomposition of monthly sea surface temperature (SST) variability in the tropical Pacific Ocean demonstrates that nearly all of the linear trends during 1950–2010 are found in two leading patterns. The first SST pattern is strongly related to the canonical El Niño-Southern Oscillation (ENSO) pattern. The second pattern shares characteristics with the first pattern and its existence solely depends on the presence of linear trends across the tropical Pacific Ocean. The decomposition also uncovers a third pattern, often referred to as ENSO Modoki, but the linear trend is small and dataset dependent over the full 61-year record and is insignificant within each season. ENSO Modoki is also reflected in the equatorial zonal SST gradient between the Niño-4 region, located in the west-central Pacific, and the Niño-3 region in the eastern Pacific. It is only in this zonal SST gradient that a marginally significant trend arises early in the Northern Hemisphere spring (March–May) during El Niño and La Niña and also in the late summer (July–September) during El Niño. Yet these SST trends in the zonal gradient do not unequivocally represent an ENSO Modoki-like dipole because they are exclusively associated with significant positive SST trends in either the eastern or western Pacific, with no corresponding significant negative trends. Insignificant trends in the zonal SST gradient are evident during the boreal wintertime months when ENSO events typically mature. Given the presence of positive SST trends across much of the equatorial Pacific Ocean, using fixed SST anomaly thresholds to define ENSO events likely needs to be reconsidered.     

An idealized study of the impact of extratropical climate change on El Niño–Southern Oscillation

Extratropical impacts on the tropical El Niño–Southern Oscillation (ENSO) are studied in a coupled climate model. Idealized experiments show that the remote impact of the extratropics on the equatorial thermocline through oceanic tunnel can substantially modulate the ENSO in both magnitude and frequency. First of all, an extratropical warming can be conveyed to the equator by the mean subduction current, resulting in a warming of the equatorial thermocline. Second, the extratropical warming can weaken the Hadley cells, which in turn slow down the mean shallow meridional overturning circulations in the upper Pacific, reducing the equatorward cold water supply and the equatorial upwelling. These oceanic dynamic processes would weaken the stratification of the equatorial thermocline and retard a buildup (purge) of excess heat content along the equator, and finally result in a weaker and longer ENSO cycle. This study highlights a nonlocal mechanism in which ENSO behavior is related to the extratropical climate conditions.     

Mean Flow–Storm Track Relationship and Rossby Wave Breaking in Two Types of El-Nino简

The features of large-scale circulation, storm tracks and the dynamical relationship between them were examined by investigating Rossby wave breaking （RWB） processes associated with Eastern Pacific （EP） and Central Pacific （CP） E1-Nifio. During EP E1-Nino, the geopotential height anomaly at 500 hPa （Z500） exhibits a Pacific-North America （PNA） pattern. During CP EI-Nifio, the Z500 anomaly shows a north positive-south negative pattern over the North Pacific. The anomalous distributions of baroclinicity and storm track are consistent with those of upper-level zonal wind for both EP and CP EI-Nino, suggesting impacts of mean flow on storm track variability. Anticyclonic wave breaking （AWB） oczurs less frequently in EP EI-Nino years, while cyclonic wave breaking （CWB） occurs more frequently in CP EI-Nino years over the North Pacific sector. Outside the North Pacific, more CWB events occur over North America during EP Ei-NiNo. When AWB events occur less frequently over the North Pacific during EP EI-Nino, Z500 decreases locally and the zonal wind is strengthened （weakened） to the south （north）. This is because AWB events reflect a monopoie high anomaly at the centroid of breaking events. When CWB events occur more frequently over the North Pacific under CP EI-Nino conditions, and over North America under EP EI-Nino condition, Z500 increases （decreases） to the northeast （southwest）, since CWB events are related to a northeast-southwest dipole Z500 anomaly. The anomalous RWB events act to invigorate and reinforce the circulation anomalies over the North Pacific-North America region linked with the two types of EI-Nino.     

Impact of the Equatorial Atlantic on the El Ni?o Southern Oscillation

Observations indicate that the Atlantic zonal mode influences El Ni?o Southern Oscillation (ENSO) in the Pacific, as already suggested in previous studies. Here we demonstrate for the first time using partial coupled experiments that the Atlantic zonal mode indeed influences ENSO. The partial coupling experiments are performed by forcing the coupled general circulation model (ECHAM5/MPI-OM) with observed sea surface temperature (SST) in the Tropical Atlantic, but with full air-sea coupling allowed in the Pacific and Indian Ocean. The ensemble mean of a five member simulation reproduces the observational results well. Analysis of observations, reanalysis, and coupled model simulations all indicate the following mechanism: SST anomalies associated with the Atlantic zonal mode affect the Walker Circulation, driving westward wind anomalies over the equatorial Pacific during boreal summer. The wind stress anomalies increase the east-west thermocline slope and enhance the SST gradient across the Pacific; the Bjerknes positive feedback acts to amplify these anomalies favouring the development of a La Ni?a-like anomalies. The same mechanisms act for the cold phase of Atlantic zonal mode, but with opposite sign. In contrast to previous studies, the model shows that the influence on ENSO exists before 1970. Furthermore, no significant influence of the Tropical Atlantic on the Indian Monsoon precipitation is found in observation or model.     

Mechanisms of tropical Pacific interannual-to-decadal variability in the ARPEGE/ORCA global coupled model

Spatial and temporal structures of interannual-to-decadal variability in the tropical Pacific Ocean are investigated using results from a global atmosphere–ocean coupled general circulation model. The model produces quite realistic mean state characteristics, despite a sea surface temperature cold bias and a thermocline that is shallower than observations in the western Pacific. The periodicity and spatial patterns of the modelled El Niño Southern Oscillations (ENSO) compare well with those observed over the last 100 years, although the quasi-biennial timescale is dominant. Lag-regression analysis between the mean zonal wind stress and the 20°C isotherm depth suggests that the recently proposed recharge-oscillator paradigm is operating in the model. Decadal thermocline variability is characterized by enhanced variance over the western tropical South Pacific (~7°S). The associated subsurface temperature variability is primarily due to adiabatic displacements of the thermocline as a whole, arising from Ekman pumping anomalies located in the central Pacific, south of the equator. Related wind anomalies appear to be caused by SST anomalies in the eastern equatorial Pacific. This quasi-decadal variability has a timescale between 8 years and 20 years. The relationship between this decadal tropical mode and the low-frequency modulation of ENSO variance is also discussed. Results question the commonly accepted hypothesis that the low-frequency modulation of ENSO is due to decadal changes of the mean state characteristics.     

Effects of interannual salinity variability on the barrier layer in the western-central equatorial Pacific: A diagnostic analysis from Argo

ABSTRACT In this paper, interannual variations in the barrier layer thickness （BLT） are analyzed using Argo three-dimensional temperature and salinity data, with a locus on the effects of interannually varying salinity on the evolution of the El Nifio Southern Oscillation （ENSO）. The interannually varying BLT exhibits a zonal seesaw pattern across the equatorial Pacific during ENSO cycles. This phenomenon has been attributed to two different physical processes. During E1 Nifio （La Nifia）, the barrier layer （BL） is anomalously thin （thick） west of about 160°E, and thick （thin） to the east. In the western equatorial Pacific （the western part： 130°-160°E）, interannual variations of the BLT indicate a lead of one year relative to those of the ENSO onset. The interannual variations of the BLT can be largely attributed to the interannual temperature variability, through its dominant effect on the isothermal layer depth （ILD）. However, in the central equatorial Pacific （the eastern part： 160~E- 170~W）, interannual variations of the BL almost synchronously vary with ENSO, with a lead of about two months relative to those of the local SST. In this region, the interannual variations of the BL are significantly affected by the interannually varying salinity, mainly through its modulation effect on the mixed layer depth （MLD）. As evaluated by a onedimensional boundary layer ocean model, the BL around the dateline induced by interannual salinity anomalies can significantly affect the temperature fields in the upper ocean, indicating a positive feedback that acts to enhance ENSO.     

Indian Ocean dipole interpreted in terms of recharge oscillator theory

In this paper we use sea surface height (SSH) derived from satellite altimetry and an analytical linear equatorial wave model to interpret the evolution of the Indian Ocean Dipole (IOD) in the framework of recharge oscillator theory. The specific question we address is whether heat content in the equatorial band, for which SSH is a proxy, is a predictor of IOD development as it is for El Niño and the Southern Oscillation (ENSO) in the Pacific. We find that, as in the Pacific, there are zonally coherent changes in heat content along the equator prior to the onset of IOD events. These changes in heat content are modulated by wind-forced westward propagating Rossby waves in the latitude band 5°–10°S, which at the western boundary reflect into Kelvin waves trapped to the equator. The biennial character of the IOD is affected by this cycling of wave energy between 5° and 10°S and the equator. Heat content changes are a weaker leading indicator of IOD sea surface temperature anomaly development than is the case for ENSO in the Pacific though because other factors are at work in generating IOD variability, one of which is ENSO forcing itself through changes in the Walker Circulation.     

Effectiveness of the Bjerknes stability index in representing ocean dynamics

The El Niño-Southern Oscillation (ENSO) is a naturally occurring coupled phenomenon originating in the tropical Pacific Ocean that relies on ocean–atmosphere feedbacks. The Bjerknes stability index (BJ index), derived from the mixed-layer heat budget, aims to quantify the ENSO feedback process in order to explore the linear stability properties of ENSO. More recently, the BJ index has been used for model intercomparisons, particularly for the CMIP3 and CMIP5 models. This study investigates the effectiveness of the BJ index in representing the key ENSO ocean feedbacks—namely the thermocline, zonal advective, and Ekman feedbacks—by evaluating the amplitudes and phases of the BJ index terms against the corresponding heat budget terms from which they were derived. The output from Australian Community Climate and Earth System Simulator Ocean Model (a global ocean/sea ice flux-forced model) is used to calculate the heat budget in the equatorial Pacific. Through the model evaluation process, the robustness of the BJ index terms are tested. We find that the BJ index overestimates the relative importance of the thermocline feedback to the zonal advective feedback when compared with the corresponding terms from the heat budget equation. The assumption of linearity between variables in the BJ index formulation is the primary reason for these differences. Our results imply that a model intercomparison relying on the BJ index to explain ENSO behavior is not necessarily an accurate quantification of dynamical differences between models that are inherently nonlinear. For these reasons, the BJ index may not fully explain underpinning changes in ENSO under global warming scenarios.     

A recipe for simulating the interannual variability of the Asian summer monsoon and its relation with ENSO

This study investigates how accurately the interannual variability over the Indian Ocean basin and the relationship between the Indian summer monsoon and the El Niño Southern Oscillation (ENSO) can be simulated by different modelling strategies. With a hierarchy of models, from an atmospherical general circulation model (AGCM) forced by observed SST, to a coupled model with the ocean component limited to the tropical Pacific and Indian Oceans, the role of heat fluxes and of interactive coupling is analyzed. Whenever sea surface temperature anomalies in the Indian basin are created by the coupled model, the inverse relationship between the ENSO index and the Indian summer monsoon rainfall is recovered, and it is preserved if the atmospherical model is forced by the SSTs created by the coupled model. If the ocean model domain is limited to the Indian Ocean, changes in the Walker circulation over the Pacific during El-Niño years induce a decrease of rainfall over the Indian subcontinent. However, the observed correlation between ENSO and the Indian Ocean zonal mode (IOZM) is not properly modelled and the two indices are not significantly correlated, independently on season. Whenever the ocean domain extends to the Pacific, and ENSO can impact both the atmospheric circulation and the ocean subsurface in the equatorial Eastern Indian Ocean, modelled precipitation patterns associated both to ENSO and to the IOZM closely resemble the observations.     

Can oceanic reanalyses be used to assess recent anthropogenic changes and low-frequency internal variability of upper ocean temperature?

A multivariate analysis of the upper ocean thermal structure is used to examine the recent long-term changes and decadal variability in the upper ocean heat content as represented by model-based ocean reanalyses and a model-independent objective analysis. The three variables used are the mean temperature above the 14°C isotherm, its depth and a fixed depth mean temperature (250?m mean temperature). The mean temperature above the 14°C isotherm is a convenient, albeit simple, way to isolate thermodynamical changes by filtering out dynamical changes related to thermocline vertical displacements. The global upper ocean observations and reanalyses exhibit very similar warming trends (0.045°C per decade) over the period 1965–2005, superimposed with marked decadal variability in the 1970s and 1980s. The spatial patterns of the regression between indices (representative of anthropogenic changes and known modes of internal decadal variability), and the three variables associated with the ocean heat content are used as fingerprint to separate out the different contributions. The choice of variables provides information about the local heat absorption, vertical distribution and horizontal redistribution of heat, this latter being suggestive of changes in ocean circulation. The discrepancy between the objective analysis and the reanalyses, as well as the spread among the different reanalyses, are used as a simple estimate of ocean state uncertainties. Two robust findings result from this analysis: (1) the signature of anthropogenic changes is qualitatively different from those of the internal decadal variability associated to the Pacific Interdecadal Oscillation and the Atlantic Meridional Oscillation, and (2) the anthropogenic changes in ocean heat content do not only consist of local heat absorption, but are likely related with changes in the ocean circulation, with a clear shallowing of the tropical thermocline in the Pacific and Indian oceans.     

Impact of vertical mixing induced by small vertical scale structures above and within the equatorial thermocline on the tropical Pacific in a CGCM

Oceanic vertical mixing is known to influence the state of the equatorial ocean which affects the climate system, including the amplitude of El Niño/Southern Oscillation (ENSO). Recent measurements of ocean currents at high vertical resolution capture numerous small vertical scale structures (SVSs) within and above the equatorial thermocline that contribute significantly to vertical mixing but which are not sufficiently resolved by coarse resolution ocean models. We investigate the impact of the vertical mixing induced by the SVSs on the mean state and interannual variability in the tropical Pacific by using a coupled general circulation model. The vertical mixing induced by the SVSs is represented as an elevated vertical diffusivity from the surface down to the 20 °C isotherm depth, a proxy for the depth of the thermocline. We investigate different forms for the elevated mixing. It is found that the SVS-induced mixing strongly affect the mean state of the ocean leading to a warming of sea surface temperature (SST) and associated deepening and sharpening of the thermocline in the eastern equatorial Pacific. We find that the SST warming induced by the elevated mixing is further strengthened through the Bjerknes feedback and SST-shortwave flux feedback. We also find a reduction in the number of large amplitude ENSO events and in certain cases an increase in the skewness of ENSO.     

Interannual linkage between Arctic/North Atlantic Oscillation and tropical Indian Ocean precipitation during boreal winter

In the study authors analyzed the interannual relationship between the Arctic Oscillation (AO)/North Atlantic Oscillation (NAO) and the tropical Indian Ocean (TIO) precipitation in boreal winter for the period 1979–2009. A significant simultaneous teleconnection between them is found. After removing the El Niño/Southern Oscillation and Indian Ocean dipole signals, the AO/NAO and the TIO precipitation (0°–10°S, 60°–80°E) yield a correlation of +0.56, which is also consistent with the AO/NAO-outgoing longwave radiation correlation of ?0.61. The atmospheric and oceanic features in association with the AO/NAO-precipitation links are investigated. During positive AO/NAO winter, the Rossby wave guided by westerlies tends to trigger persistent positive geopotential heights in upper troposphere over about 20°–30°N and 55°–70°E, which is accompanied by a stronger Middle East jet stream. Meanwhile, there are anomalous downward air motions, strengthening the air pressure in mid-lower troposphere. The enhanced Arabian High brings anomalous northern winds over the northern Indian Ocean. As a result the anomalous crossing-equator air-flow enhances the intertropical convergence zone (ITCZ). On the other hand, the anomalous Ekman transport convergence by the wind stress curl over the central TIO deepens the thermocline. Both the enhanced ITCZ and the anomalous upper ocean heat content favor in situ precipitation in the central TIO. The AO/NAO-TIO precipitation co-variations in the IPCC AR4 historical climate simulation (1850–1999) of Bergen Climate Model version 2 were investigated. The Indian Ocean precipitation anomalies (particularly the convective precipitation along the ITCZ), in conjunction with the corresponding surface winds and 200 hPa anticyclonic atmospheric circulation and upper ocean heat contents were well reproduced in simulation. The similarity between the observation and simulation support the physical robustness of the AO/NAO-TIO precipitation links.     

Simulations of two types of El Niño events by an optimal forcing vector approach

In this paper, an optimal forcing vector (OFV) approach is proposed. The OFV offsets tendency errors and optimizes the agreement of the model simulation with observation. We apply the OFV approach to the well-known Zebiak–Cane model and simulate several observed eastern Pacific (EP) El Niño and central Pacific (CP) El Niño events during 1980–2004. It is found that the Zebiak–Cane model with a proper initial condition often reproduces the EP-El Niño events; however, the Zebiak–Cane model fails to reproduce the CP-El Niño events. The model may be much more influenced by model errors when simulating the CP-El Nino events. As expected, when we use the OFV to correct the Zebiak–Cane model, the model reproduces the three CP-El Niño events well. Furthermore, the simulations of the corresponding winds and thermocline depths are also acceptable. In particular, the thermocline depth simulations for the three CP-El Niño events lead us to believe that the discharge process of the equatorial heat content associated with the CP-El Niño is not efficient and emphasizes the role of the zonal advection in the development of the CP-El Nino events. The OFVs associated with the three CP-El Niño events often exhibit a sea surface temperature anomaly (SSTA) tendency with positive anomalies in the equatorial eastern Pacific; therefore, the SST tendency errors occurring in the equatorial eastern Pacific may dominate the uncertainties of the Zebiak–Cane model while simulating CP-El Nino events. A further investigation demonstrates that one of the model errors offset by the OFVs is of a pattern similar to the SST cold-tongue cooling mode, which may then provide one of the climatological conditions for the frequent occurrence of CP-El Nino events. The OFV may therefore be a useful tool for correcting forecast models and then for helping improve the forecast skill of the models.