Research, based on both observations and ocean models, has indicated that strong East Pacific El Niños are preceded by the development of unusually low sea levels at the western end of the North Equatorial Trough early in the year. This results in an increase in the strength of the North Equatorial Counter Current which, aided by low sea levels due to the annual Rossby wave, then transports West Pacific Warm Pool water to the central and eastern Pacific.
Here, an ocean model is used to investigate the factors affecting sea level prior to the 1982–1983 East Pacific El Niño. The results indicate that during this period the reduction in sea level was caused by Ekman pumping, due to the local winds, acting over a period of many months. The north–south distances involved mean that such upwelling can be connected to the westerly wind phase
of Madden–Julian Oscillations.
Since the 1980s, the quality and quantity of global wind measurement has steadily improved. So, if the Ekman pumping hypothesis is correct, the better-quality wind data available prior to the 1997–1998 and 2015–2016 East Pacific El Niños should also show unusually large Ekman pumping in
the same region, over a similar long period. This is tested and found to be correct.
However, a study of the integrated Ekman pumping for the period 1976 to 2015 indicates that in some years similar events did not develop despite a
comparable amount of pumping during the first part of the year. The results imply that significant Ekman pumping early in the year is a necessary
but not sufficient requirement for the development of a strong East Pacific El Niño.
The North Equatorial Countercurrent (NECC) flows from west to east across the Pacific at latitudes between 5 and 10∘ N. The total transport
is in the range 10–30 Sv, making it comparable to the major Northern Hemisphere currents .
Physically, it is a geostrophic current lying between the sea level minimum of the North Equatorial Trough and the sea level maximum of the North
Equatorial Ridge . However, it is also a shallow current with a depth of only 200–300 m,
the north–south pressure differences near the surface being balanced by large north–south gradients in the near-surface density surfaces
.
The current is important because it carries warm water eastwards out of the West Pacific Warm Pool (also referred to as “Warm Pool” in the following; ). It is also unusual in that it flows in the opposite direction to the steady trade winds of the North Pacific. originally suggested that the current flowed downhill from west to east across the Pacific but, using geostrophic theory, showed that it was the result of the shear in the wind field at latitudes near the Intertropical Convergence Zone (ITCZ) See also. The current also lies in a region where the inverse of the Coriolis term is a rapidly changing function of latitude. As a result, small changes in the position of the trough can have a significant effect on the velocity of the current.
Early studies of the current were limited by lack of data, but, using data from bathythermographs, showed that the depth of the
North Equatorial Trough had a significant annual signal, which he identified as being due to an annual Rossby wave.
Further developments in ocean instrumentation have provided more information on the NECC, satellite observations of sea surface temperature (SST)
showing tropical instability eddies often called waves; that both warm the
Equatorial Current and cool the NECC , the effect possibly being least in El Niño years
.
Satellite altimeter measurements of sea level have errors of only a few centimetres and, as the North Equatorial Trough has a depth of the order of 1 m, this allows the annual and interannual changes of both the trough and the NECC to be studied in detail .
Connection with the El Niño
Originally, El Niño was the name of a southward-flowing current off the coast of Peru which affected fisheries . However, the studies reported in and connected it to the very large changes in surface temperature and thermocline depth observed in the nearby Equatorial Pacific Cold Pool (also referred to as “Cold Pool” in the following).
was the first to make the further connection between this oceanographic El Niño and the strength of the NECC in the western
Pacific. He also suggested that the NECC might trigger an El Niño by transporting West Pacific Warm Pool water to the east.
then showed that the NECC transport increased by between 25 % and 50 % while the strong 1982–1983 El Niño was
developing. They also found that the total transport of warm water (> 28 ∘C) by the NECC was consistent with their estimate of the
loss of warm water from the West Pacific Warm Pool.
Although Wyrtki's hypothesis was not developed further, his work complemented the study by which showed a correlation
between water temperatures in the equatorial Pacific and the Southern Oscillation . This led to further uses of the term
“El Niño” to describe both its worldwide influence and the air–sea interaction (or mechanism) behind the event .
also argued that the east–west temperature gradient along the Equator was involved, and this is usually taken to imply that the temperature of the Equatorial Pacific Cold Pool
is a key part of the El Niño mechanism. However, doubts this, because the east Pacific winds are only weakly affected by an El Niño. If he is correct, temperature changes in the Cold Pool are a consequence of and not a cause of El Niños.
Further studies showed that Bjerknes' correlation was strongest when using ocean temperatures from the central Pacific. As a result, most modern
studies of the Southern Oscillation use the Nino 3.4 index to measure the strength of an El Niño . This is based on the mean sea surface temperature anomaly in the central Pacific between 170 and 120∘ W, and between 5∘ S and 5∘ N.
As well as affecting surface winds in the western and central Pacific, the atmospheric El Niño requires the centre of deep atmospheric convection to
move from the maritime continent and the western Pacific into the central Pacific. , following
and , confirmed that this requires sea surface temperatures of 28 ∘C and
above during wet seasons and even higher temperatures during dry periods.
In the central Pacific, such temperatures are found on the Equator during an El Niño, but eastern Pacific equatorial temperatures never reach this
value Fig. 2. Instead, the eastward limit of warm water on the Equator is closely correlated with the Southern Oscillation
Index .
In the eastern Pacific, temperatures sufficient to trigger deep atmospheric convection are often found further north in the ocean off Central America
and along the path of the NECC close to the ITCZ.
Meteorological studies also showed that the start of an El Niño was often connected with one or more Madden–Julian Oscillations (MJOs)
. In his review, describes MJOs as individual events with lifetimes of between
30 and 100 d, and with the property that there is never more than one major event in existence at a time. Each one consists of an area,
spanning the Equator, in which there is an increased density of short-lived convective cloud groups.
The events progress eastwards at speeds around 5 ms-1, starting in the western Indian Ocean and dying out after passing the western
Pacific. At sea level, they are associated with inflowing winds from the east and west, but in the Pacific the westerly winds may extend past the
centre of the cloud groups.
El Niño models and theories
Following Wyrtki and Bjerknes, most theoretical studies of the ocean's contribution to El Niño have focused on the Cold Pool. This is a narrow band
of low sea surface temperatures, centred on the Equator in the eastern Pacific. During a normal year, this results from the wind-generated upwelling of
cold water from depths of around 200 m, but during an El Niño the surface thermocline becomes thicker and less cold water is upwelled.
Early theoretical studies, using analytic models or numerical models with a few layers , showed that the changes in the Cold Pool could be generated by wind variations in the western and central
Pacific. This occurred via the propagation of baroclinic Kelvin waves in the equatorial wave guide. The studies also showed that only winds close to
the Equator were responsible and that off-equatorial currents, such as the NECC, were not connected with the El Niño changes
to the Cold Pool.
With the development of computer power, it was eventually possible to run ocean models with much more detailed physics and with enough horizontal and
vertical resolution to resolve the key currents and structure of the ocean . The versions developed by , and have been widely used in
stand-alone mode and, coupled to an atmospheric model, for studies of climate.
Unfortunately, the NECC can be poorly represented in such models . For ocean-only models,
found that the poor results were largely due to errors in the atmospheric forcing datasets.
concentrated on the forcing of the NECC and showed that a weak NECC resulted from the atmospheric reanalysis datasets
underestimating the curl of the wind field near the latitudes of the ITCZ and overestimating the strength of easterly winds on
the Equator.
The study also found that, of the datasets studied, the ECMWF reanalysis was best at generating a realistic NECC. This may be because the ECMWF model
was using a spherical harmonic expansion which requires less smoothing than normal finite difference schemes to overcome non-linear instabilities.
There is a related problem with coupled ocean–atmosphere models which fail to generate realistic El Niños . This again may be due to the resolution of the atmospheric component being insufficient to
accurately represent features such as the ITCZ.
Further studies using observational data have indicated that there are two types of El Niño . The most common are the central Pacific (CP) or Modoki events, in which the centre of atmospheric convection moves to the central
Pacific. More rarely, the centre of convection moves even further east to give the strong eastern Pacific (EP) El Niños . These
behave like the classic El Niños , which proposed were triggered by the NECC.
Wyrtki's hypothesis
Although neglected for many years, Wyrtki's hypothesis is supported by more recent studies using altimeter data .
These confirm the increased transport by the NECC during EP El Niños, showing that increases in the transport of the NECC
in the western Pacific occur around 3 months before the corresponding change in the Nino 3.4 index. Also in the western Pacific,
show that there is a delay of between 40 and 90 d between changes in the depth of the North Equatorial Trough, north
of the NECC, and changes in sea level near the Equator.
The results thus leave open the possibility that the increased transport has a role in triggering subsequent EP El Niños, rather than being part of
the ocean's response to the event. The studies also indicate similar changes in the NECC during CP El Niños, indicating that both types of El Niño
may be affected by the NECC in the western Pacific.
Further support for the hypothesis comes from , who used results from a high-resolution global version of the NEMO ocean model to study the strong EP El Niños of 1982–1983 and 1997–1998.
For this study, NEMO was forced by the DRAKKAR dataset , based on the ECMWF reanalysis fields, and produced sea level and sea surface temperature fields in the equatorial Pacific in close agreement with satellite observations . The results indicate that the modified ECMWF reanalysis fields accurately represented the wind stress curl and so generated a realistic NECC.
Unlike previous studies, concentrated on water with temperatures above 28 ∘C, which are sufficient to trigger deep
atmospheric convection. The study showed that, while the strong El Niños were developing, the NECC carried warm water into the eastern Pacific far
in advance of water with similar temperatures on the Equator.
The study also helped explain why the NECC does not carry Warm Pool water into the eastern Pacific every year, by showing that in normal years tropical instability eddies mix in cold upwelled water from near the Equator and that this rapidly reduces temperatures in the NECC See:. In El Niño years, the eddies are weaker so their effect on the NECC is reduced.
It also found that the start time of strong El Niño events resulted in warm water reaching the central and eastern Pacific at a time of year when
the annual Rossby wave meant that the NECC was strongest there. This resulted in the warm water arriving off South America around the end of the year,
thus providing a physical mechanism to explain the name originally given to classic EP El Niños.
Finally, the study showed that both the strong 1982–1983 and 1997–1998 El Niño events started when the depth of the western section of the North Equatorial Trough was unusually low. As a result, the NECC transport, and the flux of Warm Pool water eastwards, were both much stronger than usual.
Low sea levels in the western Pacific
suggested that the low sea level in the trough was due to an unusually strong annual Rossby wave prior to each El Niño. However,
further investigation showed that this is unlikely, as in each case the wave is not unusually large as it crosses the central Pacific.
Another possible explanation is that the stratification in the Pacific, prior to an El Niño, is sufficient to focus the Rossby wave and so cause
lower-than-normal sea levels. Alternatively, the low sea level may not be due to the Rossby wave but instead due to some other process associated
with the wind field. This could be a local feature occurring at or just before the period when sea level drops. It could also be a response to forcing
elsewhere in the ocean, which later propagates into the western Pacific.
To help clarify the cause of the drop in sea level, this paper starts by reporting on series of short ocean model runs which focus on changes in the
western Pacific prior to the 1982–1983 El Niño. As in , this allows the development of hypotheses, which are then tested by
treating them as predictions for the strong El Niños of 1997–1998 and 2015–2016 for which better observations are available.
Occam, the global ocean model used for the tests, is an updated version of the OCCAM project ocean model . In each run, the model is initialised from one of the archive datasets from the original run of the NEMO model and forced with ocean surface stresses calculated during the same run.
The results imply that the low sea levels, at the western end of the North Equatorial Trough, are due to local Ekman pumping acting over a period of
many months.
This hypothesis is tested against the data from the original run of the NEMO 1/12∘ global ocean model, initially during the growth of the
1982–1983 El Niño and then for the strong EP El Niños of 1997–1998 and 2015–2016. A similar comparison covering the whole period from 1976 to 2015
is also carried out.
Structure of the report
In the first part of the paper, Sect. 2 describes the model being used and how the model fields and forcing were converted from the
1/12∘ NEMO grid to the 1/4∘ Occam grid.
Section 3 reports on tests carried out to validate the lower-resolution model. In these, the model was started from the NEMO archive datasets from
early January 1981 and 1982 and then run for a year using the surface wind stresses from the same year.
The following sections report tests with different wind forcing and initial conditions. This includes runs where the model was again started from
early January 1981 and 1982 but with the wind forcing from the opposite year. The results show that the main changes in sea level depend primarily on
the wind field and, to first order, are independent of the stratification at the start of the run.
The tests are repeated but starting later in the year when the annual Rossby wave has developed and is starting to cross the central Pacific. This is
to check whether it is responsible for the sea level drop in the western Pacific. The results are similar to the earlier test, indicating that the
wind-generated Rossby wave is not responsible.
A final test is designed to investigate whether the winds causing the drop in sea level are local to the western equatorial Pacific or propagate in
from other parts of the ocean. The results indicate that it is the local winds that are responsible.
The most likely way in that the local winds could affect sea level in the North Equatorial Trough is through Ekman pumping and the resulting rise in
density surfaces within the ocean. In the second part of the paper, this hypothesis is tested.
Section 7 compares the expected uplift due to Ekman pumping during 1982 with the rise in density surface in the NEMO 1/12∘ model. The results
are promising so in Sect. 8, similar tests are also carried out for the 1997–1998 and 2015–2016 EP El Niños.
Finally, a similar survey of the whole period 1976 to 2015 is carried out to see if strong Ekman pumping in the western section of the North Equatorial Trough during the early part of the year is common and whether it might also be a sufficient condition for the development strong EP El Niños.
The Occam 1/4∘ global ocean model
Occam is a primitive equation model, based on the Bryan–Cox–Semtner series of models . It uses a regular
latitude–longitude grid for all the oceans except the North Atlantic and Arctic. For the latter basins, it uses a second rotated latitude–longitude
grid which is matched to the first grid at the Equator. The Bering Strait between the Arctic and the North Pacific is represented by a simple channel model.
Occam uses 1/4∘ resolution in both longitude and latitude. In the vertical, it has 66 levels, instead of the 75 levels of the NEMO run, and
makes use of an existing global topography that was checked against a database of critical oceanographic sills . Both the Occam
and NEMO models have increased vertical resolution in the surface layers, Occam using 24 layers in the top 300 m and NEMO 34 layers, the
difference being primarily due to NEMO's extra resolution in the top 100 m.
Occam uses harmonic mixing in the horizontal and the scheme of for vertical mixing. For horizontal advection, it uses
the second-order split-quick scheme for both momentum and tracers. For the vertical advection of momentum, it uses the scheme of
.
Occam was chosen for these tests primarily because of its computational efficiency and because the amount of computer time available was limited. The
efficiency arises partly because, unlike NEMO, it uses a regular grid and partly because the code includes fewer of the complex physical options included in NEMO.
However, most ocean models are limited not by the speed of the processor but by the length of time needed to transfer data from main memory. Occam's
main advantage is that it overcomes this by vectorising the code in the vertical
This means that all the variables needed to time step the
cells in a vertical column are held in high-speed cache at the same time and no extra references to main memory are required. Also, when moving from
one column of ocean cells to the next, the cache already contains most of the required data.
.
Initialisation and forcing
The model runs reported here were initialised by averaging the archived data from the high-resolution NEMO model onto the Occam grid. Variable values
within each Occam ocean cell were calculated by averaging over the intersection of the Occam and NEMO cells. In the case of vector quantities, vectors
were rotated to the Occam grid before averaging.
The runs were also carried out with zero flux of heat and salt across the ocean surface. Each run is for only a few months and, although the surface
temperatures and salinities are affected, this approximation allows the analysis to concentrate on the effects of wind stress, and of advection and
diffusion within the ocean.
Each NEMO archive dataset contains the ocean and forcing fields averaged over the previous 5 d. Thus, when initialising the Occam model and
when specifying the wind forcing, the time of each NEMO dataset is set to the central time of the averaging period. However, when specifying particular
NEMO archive files, the original archive date is used.
The use of 5 d averages filters out higher-frequency components of the wind field and also means that the ocean state used to initialise the
Occam model may be unusually smooth. The resulting lack of high-frequency oscillations may affect vertical mixing in the ocean but should otherwise
have little effect on lower-frequency variations in sea level and current velocity, which are the focus of this study.
Validation
Previously, both the NEMO and Occam models have been widely used for successful oceanographic research, but here the use of a lower-resolution Occam
model to investigate features seen in the higher-resolution NEMO run is validated by comparing model results in two runs. Each run was over a full
year and each used the NEMO-archived surface wind stresses from that year.
In the first run, designed to compare performance in a normal year, Occam was initialised from the NEMO archive datasets from 5 January 1981. In
the second, a test of a strong El Niño year, the model was started from 5 January 1982.
The use of 1981 as a representative normal year has the advantage that, there being no other known special considerations that might affect the
results, to a first approximation, developments in 1981 should reflect what would have occurred in 1982 had there been no El Niño.
Sea level (m) plus surface currents on 6 June and 9 September for (a) NEMO in 1981, (b) Occam in 1981, starting in January 1981, (c) Occam in 1982, starting in January 1982 but with 1981 winds.
Comparison of 1981 results
In Fig. , the top two pairs of figures show sea level from the NEMO and Occam models on 6 June and 9 September 1981. 6 June is
chosen because the corresponding date in 1982 can be used to check the early deepening of the North Equatorial Trough in the western Pacific that
occurs prior to the 1982–1983 El Niño. Similarly, 9 September is chosen as a check on how well Occam reproduces the increased depth of the
trough, and corresponding increased strength of the NECC, across the entire Pacific, during the autumn months.
At the largest scales, the Occam model is in close agreement with NEMO, the main differences arising from slightly more extreme maxima and maxima. Thus,
in September, maximum sea level in the South Pacific Gyre is slightly higher in Occam. Similarly, the minimum sea levels in the Cold Pool on the Equator and within the North Equatorial Trough are slightly lower.
This paper is primarily concerned with the behaviour of the NECC and the associated sea levels in the North Equatorial Trough. Here, Occam shows a
similar westward propagation of the trough during the year and a similar growth of the ridge on the southern side of the NECC. The September figures
also show a similar development in Occam of the short wavelength features on the ridge, associated with the growth of tropical instability eddies.
In June, both NEMO and Occam show that the trough is particularly weak in the central Pacific, resulting in an almost non-existent NECC. However, in
September, with a greater contrast between ridge and trough, both models show a stronger NECC.
At the western ends of the trough and ridge, Occam shows a similar development of the high- and low-sea-level regions associated with both the
Halmahera Eddy and Mindanao Eddy and the initial meanders of the NECC. In September, these decay in amplitude, in
contrast to the growth of the trough and ridge features seen in the central and eastern Pacific.
The development of sea level during the year at 6∘ N in both NEMO and Occam is shown in Fig. a and c.
The average slope in sea level in Occam matches that of NEMO, and in the central and eastern Pacific, Occam shows a similar development of the ridges
due to the growth and propagation of tropical instability eddies. The figure also shows the low sea levels due to the propagation on the annual Rossby
wave at 6∘ N. This starts in the eastern Pacific early in the year and reaches the dateline in July and August.
The corresponding sea levels on the Equator are shown in Fig. a and c. Again, the overall east–west slope is similar and Occam reproduces the eastward-propagating equatorial waves that can be generated by westerly wind bursts in the western equatorial Pacific. There are differences between the two model runs but in both models sea level is a maximum either close to the western boundary or near 150∘ E.
Surface temperature (∘C) on 6 June and 9 September for (a) NEMO in 1981, (b) Occam in 1981, starting in January 1981, (c) Occam in 1982, starting in January 1982 but with 1981 winds.
Sea surface temperatures, in June and September 1981, are shown in the top pairs of panels in Fig. . As Occam includes no surface warming
or cooling, it is expected to show no seasonal changes. This is seen to be the case, with Occam showing warmer temperatures in the Southern Hemisphere in
both June and September and cooler temperatures in the north.
The differences are largest off Central America, where Occam fails to reproduce the East Pacific Warm Pool, and in the western section of the South Pacific, where there is little reduction of Occam temperatures between June and September. The regions of cold water, with temperatures at or below 20 ∘C, also behave very differently in Occam and NEMO.
Along the line of the NECC at 6∘ N, Occam fails to show the warming during the first half of the year, but by September it does show a narrow
band of warm water being advected into the eastern Pacific by the NECC. It also shows it being eroded by tropical instability eddies.
Sea level (m) plus surface currents on 6 June and 9 September for (a) NEMO in 1982, (b) Occam in 1982, starting in January 1982 (c) Occam in 1981, starting in January 1981 but with 1982 winds.
Surface temperature (∘C) on 6 June and 9 September for (a) NEMO in 1982, (b) Occam in 1982, starting in January 1982 (c) Occam in 1981, starting in January 1981 but with 1982 winds.
Hovmöller diagrams of sea level (m) at 6∘ N in the Pacific for (a) NEMO model in 1981, (b) NEMO in 1982, (c) Occam starting in January 1981, (d) Occam starting in January 1982, (e) Occam starting in January 1982 but with 1981 winds, (f) Occam starting in January 1981 but with 1982 winds. The figures are based on 1∘ averages of the model data.
Hovmöller diagrams of sea level (m) on the Equator for (a) NEMO model in 1981, (b) NEMO in 1982, (c) Occam starting in January 1981, (d) Occam starting in January 1982, (e) Occam starting in January 1982 but with 1981 winds, (f) Occam starting in January 1981 but with 1982 winds. The figures are based on 1∘ averages of the model data.
Comparison of 1982 results
Corresponding sets of panels for 1982, during the development of the 1982–1983 El Niño, are shown in Figs. –.
In Fig. , Occam shows, in agreement with NEMO, a much deeper North Equatorial Trough in early June together with a well-developed NECC. This is emphasised further in September with a deep trough in the west which extends almost as far as the Mindanao Eddy.
The Hovmöller diagrams of Fig. show that in 1982 the annual Rossby wave, in both Occam and NEMO, links with the region of low sea level
that develops at the western end of the trough. argued that it was this drop which generated the increase the NECC strength which
triggered the development of a classic oceanographic El Niño.
At the Equator (Fig. ), Occam shows, in agreement with NEMO, the movement of the maximum sea level away from the western boundary to a
region near 200∘ E. In Figs. and , this is seen to correspond to the region of high sea level on the Equator where
the warmest equatorial temperatures are also found.
Summary of comparisons
The results show that the Occam model is capturing the key features in the sea level field which affect the strength of the NECC, the strength of
tropical instability eddies and processes occurring on the Equator. It is less effective at capturing changes to sea surface temperature field. This
may partly be due to the lack of surface heat fluxes in this version of the model. It may also be partly due to the much reduced horizontal resolution
of the model.
However, bearing these strengths and weaknesses into account, there is no reason why it cannot be used to study the effect of the winds and the
initial state of the ocean on the development of the 1982–1983 El Niño, as is done here.
Tests using wind stresses from different years
Although identified a number of processes contributing to the development of a strong El Niño in 1982 and 1997, the study was not
able to explain why the events started in those years and not, say, a year earlier or a year later.
One possibility, suggested by , is that there was something different about the stratification of the ocean in early 1982 and 1997,
which focused the annual Rossby waves, so that they had a larger amplitude than usual once they reached the western Pacific. Another possibility is
that the difference was due to the wind. This may have generated a stronger-than-normal annual Rossby wave that year. Alternatively, some independent
wind event may have occurred which lowered sea level in the western Pacific around the time that the leading edge of the annual Rossby wave arrived.
In an attempt to distinguish between these possibilities, the two Occam runs starting in January 1981 and 1982 were repeated, but this time they were forced by
the winds from the opposite year. If focusing is important, then the January 1982 ocean, forced by the 1981 winds, might generate a similar enhanced
annual Rossby wave in the western Pacific. Alternatively, if the winds are important, then the January 1981 run, forced by the 1982 winds, might
generate the low sea level in the western Pacific.
January 1982 ocean forced by 1981 winds
In this test, the model was initialised from the NEMO archive dataset dated 5 January 1982 but then forced with winds from
1981. Figure e shows how sea level developed at 6∘ N. During the first few months of the year, the annual Rossby wave develops and
propagates westwards as normal but in the middle of the ocean the amplitude declines, so by the time the wave reaches 180∘ E the signal is weak.
On the Equator (Fig. e), sea level starts high in the western Pacific, as in the earlier run with 1982 winds forcing the 1982 ocean
(Fig. d), but as the year develops it stays in the west and there is no movement into the central Pacific.
Figures c and c show the sea level and temperature fields for June and September. Previous runs with 1982 winds show that by
June the depth of the trough is starting to increase more than in a normal year but here this is not happening and instead the response is closer to
that of NEMO in 1981. Similarly, in September, the trough fails to develop further and the response again lies closer to that of NEMO in 1981.
The surface temperature field is also closer to the NEMO results from 1981, with the bulk of the Warm Pool water remaining in the west. However, there
is a thin core of warmer NECC water extending into the eastern Pacific, which is not present in the run starting in January 1981 with 1981 winds. This
indicates that there was some change in the structure of the ocean between January 1981 and 1982 which aided the transport of warm water eastwards by
the NECC.
Despite this, the main conclusion from this run is that the ocean state at the start of 1982 was not sufficient to trigger an El Niño in the
middle of the year. Although a reasonable annual Rossby wave was generated early in the year, this was not focused and did not generate or contribute to the
lowering of sea level in the western Pacific.
January 1981 ocean forced by 1982 winds
In the second test, the ocean is initialised from the NEMO-archived dataset dated 5 January 1981 but then forced with 1982
winds. Figure f shows the sea level during the year at 6∘ N. The annual Rossby starts as before, but this time it continues past
180∘ E and links up with a region of low sea level in the western Pacific – as might be expected at the start of a strong El Niño.
On the Equator (Fig. f), the ocean again starts with high sea levels in the west. Mid-year, these move into the central Pacific – again
as might be expected at the start of a strong El Niño.
Figures c and c show that in the middle of the year the distribution of sea surface temperature is similar to the other Occam runs, but sea
level shows a well-developed North Equatorial Trough with a strong NECC developing on its southern slope.
By September, the trough has developed further and the temperature plot shows warm surface water being advected rapidly into the eastern Pacific. The
area of warm water involved is not as large as in the original NEMO run in 1982, but as with Fig. b, this is probably a result of setting
the surface flux of heat to zero in these test runs.
One conclusion from this test is that the ocean will produce an El Niño-like response when forced by the winds from an El Niño year. This has been reported before but based primarily on the large sea level and temperature anomalies that develop in the Equatorial Pacific Cold Pool towards the end of the first year of a classic oceanographic El Niño.
The present study is different in that it shows the winds having a significant effect during the first half of the year, well before any changes in
the Cold Pool region are noticeable. The results also show that by the middle of the year the 1982 winds have produced significant changes in the western section of
the North Equatorial Trough. What is not clear is whether the increased depth of the trough is due to the mean winds during the first part of the year
or due to one or more isolated events.
It is also not clear whether it is due to a local change in the winds along the line of the trough, or whether the key wind forcing occurs elsewhere
and the signal propagates into the region lowering the level of the trough. The remaining model runs reported here attempt to obtain a clearer answer
to these questions.
Tests starting in late April
One possibility that has not been discounted is that, early in 1982, the winds generated a Rossby wave, or similar, which was later responsible for
the sea level drop in the western Pacific. To test this hypothesis, the model was started from the NEMO archive dataset dated 30 April 1982 and
then forced with the 1981 winds. The start date is just before the start of the drop in sea level in the western Pacific. If the Rossby or other waves
are responsible, then by that date they should be established enough to reproduce the drop in sea level despite the change in the wind field.
Hovmöller diagram of sea level (m) at 6∘ N during 1982 for (a) when started from the NEMO model archive from 30 April 1982 but forced by 1981 winds and (b) when started from the NEMO model archive from 30 April 1981 but forced by 1982 winds.
Figure a shows how sea level develops along 6∘ N. The annual Rossby wave starts propagating across the central Pacific as usual,
but near the dateline its amplitude is greatly reduced and there is no connection with the lower sea levels of the far western Pacific. Thus, although
the test does not exclude events occurring before 30 April having some impact, they cannot be the prime cause of the sea level drop that increases
the strength of the NECC near 160∘ E.
The complementary test was also carried out, in which the ocean was started from 30 April 1981 but forced with 1982 winds. The result, shown in
Fig. b, shows that sea level does drop. This implies that it is the winds after 30 April 1982 that are responsible.
Mean wind stress vector anomaly (Nm-2), for the period 16 March to 8 August 1982 relative to the same period in 1981. The blue rectangle shows the region between 140 and 170∘ E and between 5 and 7∘ N.
The role of local winds in the western Pacific
Figure shows the difference between the 1982 and 1981 wind stress vectors for the Pacific when averaged between 16 March and
8 August. It shows that in the western Pacific near the Equator, there is a significant westerly component to the wind stress anomaly. North of New
Guinea, this drops to zero near 10∘ N. This distance is typical of the atmospheric equatorial Rossby radius, a scale which also determines the
northward extent of Madden–Julian Oscillations .
MJOs are reported to be stronger over the western Pacific prior to eastern Pacific El Niños . Thus the anomaly and
the large north–south gradient in the wind stress may be connected with the MJOs.
Further east, the ITCZ often lies close to 10∘ N. It is not clear whether this is related to the equatorial Rossby radius, but the figure
shows that this is also a region where both the wind stress anomaly and its north–south gradient can be large.
To see if the local wind fields in the western Pacific were responsible for the drop in sea level in 1982, a modified wind field was constructed which
combined both the 1981 and 1982 winds. The weighting for the 1982 winds was defined so that it equalled 1 within the region 140 to 180∘ E
and the Equator to 15∘ N and was zero outside the region 130∘ E to 190∘ W and 10∘ S to 25∘ N. Linear
interpolation was used between the two boundaries. Weighting for the 1981 field was set to 1 minus the 1982 weighting.
This results in the 1982 wind field being used for the area where the main sea level drop is observed and the 1981 wind fields being used elsewhere
but with a smooth transition zone. The combined wind field was then used for two runs. As with the previous tests using 1982 winds, the first started
from the ocean state for 3 January 1981 and the second from 30 April 1981. If the drop in sea level is due to local winds, then the change in
sea level should be similar to those obtained using only the 1982 wind field.
Hovmöller diagram of sea level (m) at 6∘ N during 1981, forced by the combination of 1981 and 1982 winds described in the main text when started from the NEMO archive for (a) 5 January 1981 and (b) 30 April 1981.
The resulting sea levels along 6∘ N are shown in Fig. . In both cases, sea level starts dropping in the western Pacific in early
May, confirming the importance of the local winds after 30 April.
Sea level (m) in the different runs averaged between 140 and 170∘ E and between 5 and 7∘ N. Solid red lines correspond to runs forced by 1982 winds; blue correspond to 1981 winds. Thick lines are from the original NEMO run; thin lines are from the Occam runs. The dashed red lines correspond to the runs started in 1981 but with 1982 winds only in the western Pacific.
Intercomparison
Figure shows the changes in sea level from the different runs, averaged over the region 140 to 170∘ E and 5 to 7∘ N, as
outlined in Fig. . The Occam runs start with a slight offset from the corresponding runs of the NEMO model, due to the initial adjustment
to the new coast and topography.
In this figure, colour is used to distinguish the different winds forcing the ocean: blue for 1981 and red for 1982. Even with the different starting
conditions, the runs soon split into two distinct groups: sea levels with the 1981 winds staying roughly constant, whereas the 1982 winds generate a
20 cm drop between March and August. As pressures in the deep ocean remain roughly constant, this drop in sea level must correspond to a
significant rise in the density surfaces within the ocean.
Mean sea levels (m) between 140 and 170∘ E and between 5 and 7∘ N, at the end of April and August for the different runs, identified by the model and start month. The two groups correspond to runs forced by the 1981 winds and by either the 1982 winds or the combined winds.
To quantify the differences between the groups, sea levels from the end of April and August are given in Table , split into runs forced
by 1981 winds and those forced either by 1982 winds or by the combined wind fields. A t test shows that the probability of both sets belonging to the
same group at the end of April is around 0.5, as might be expected. By the end of August the probability drops to around 0.00002. A test using the sea
level change over the period gives a similar value.
Winds and Ekman divergence
Given that local winds appear to be responsible, the most likely cause of the rise in the density surfaces is that it is due to Ekman pumping, the
result of a divergence in the wind-generated Ekman transport in the surface layers of the ocean.
Here, the hypothesis is investigated by calculating the Ekman pumping in the region 140 to 170∘ E and 5 to 7∘ N.
If the Ekman transport vector is E(θ,ϕ), where θ and ϕ are latitude and longitude, then
E(θ,ϕ)=τ(θ,ϕ)∧n^/(ρf(θ)),
where τ(θ,ϕ) is the wind stress vector, n^ the unit vertical vector, ρ is the density and f(θ) is the Coriolis term, equal to
f(θ)=2Ωsin(θ).Ω is the angular rotation rate of the Earth.
The averaged vertical velocity in a region, due the divergence of the Ekman flux, is then given by integrating the outward-flowing Ekman transport
around the region of interest:
P=-(1/A)∮E(θ,ϕ)∧n^⋅ds^,=(1/A)∮τ(θ,ϕ)⋅ds^/(ρf(θ)),
where s^ is the unit vector tangential to the boundary and A the area enclosed.
Integral over time of the vertical velocity due to Ekman divergence in the region bounded by 5–7∘ N and 140–170∘ E, during (black) 1982, (red) 1981 and (green) climatology 1958–2015. In the background are the corresponding plots for the years prior to the 1997–1998 and 2015–2016 El Niños.
Figure shows the results obtained by integrating the vertical velocity over time. It shows that in 1982 the Ekman divergence (given by
the slope of the curve) was positive for most of the period between April and late November, and that it had the potential for raising density
surfaces within the ocean by 80 m. During the period mid-March to mid-August, the potential rise is approximately 50 m.
A noticeable feature of Figs. and is that, during 1982, the drop in sea level and the vertical displacement due to the Ekman
divergence are both relatively steady processes occurring over many months. There are short periods when the rate of change is reduced or reverses,
but overall the results imply that the changes are a result of a long-term systematic change in the wind field.
NEMO 1/12∘ model potential density averaged between 140 and 170∘ E and between 5 and 7∘ N for (a) 1981 and (b) 1982. Contours are at integer values of density (kgm-3).
To check that the Ekman pumping is sufficient to cause the observed change in sea level, Fig. shows density profiles from the original
NEMO 1/12∘ model run averaged over the same area as before.
Figure a shows that during 1981, the changes in the depth of the density surfaces were small and at most depths had balanced out by the
end of the year.
In contrast, during 1982 (Fig. b), there was a significant shallowing of the density surfaces. The figure shows that in mid-March, water
with densities of 1025 kgm-3 (approximately 21.5 ∘C) lay near 150 m and that by mid-August it had risen to near
100 m. This rise of approximately 50 m is comparable to the Ekman pumping estimate. However, after this period, although the negative
Ekman pumping continued, there was no further shallowing in the density surfaces near the surface and there is evidence of a rebound at depth.
The 1997–1998 and 2015–2016 El Niños
The results of the previous section support the hypothesis that Ekman pumping over many months was responsible for the drop of sea level at the start
of the 1992–1993 El Niño. However, if the mechanism is a significant part of all East Pacific El Niños, then similar behaviour should occur prior
to the 1997–1998 and 2015–2016 events. These are also periods when better wind field datasets were available, so they should be better tests of the hypothesis.
NEMO 1/12∘ model potential density averaged between 140 and 170∘ E and between 5 and 7∘ N for (a) 1997 and (b) 2015. Contours are at integer values of density (kgm-3).
The results, shown in Figs. and , give support to the hypothesis, as in each case the vertical displacement is similar to the
1982 value. In both cases, pumping continues almost continuously over many months, the only key difference being that in both cases Ekman pumping
starts earlier in the year.
Figure also contains the results for 1996 and 2014. The results for 1996, like those for 1981, are similar to the climatology
values. However, in 2014, the displacement curve starts by following the curves for 1997 and 2015, and it is only after April that the pumping is
reduced. The displacement then continues to increases until the end of the year.
A number of papers, including , and , suggest that an eastern
Pacific El Niño started to develop in 2014 but was then prevented by a change in the winds. It is possible that the 2014 curve in Fig.
reflects this event.
The main conclusion to be taken from the results is that, between mid-March and mid-August 1982, there was a large amount of Ekman pumping in the
western Pacific around 6∘ N due to a wind anomaly. Similar amounts of Ekman pumping also occurred prior to the strong 1997–1998 and
2015–2016 El Niños. In each case, the results indicate that the pumping was due to a systematic change in the local wind field which lasted for many
months.
Developments between 1976 and 2016
By themselves, the above results imply that strong Ekman pumping is always associated with strong El Niños and that in other years the integrated
Ekman pumping remains small. However, this is not correct.
(a) Integrated Ekman pumping, based on wind observations, and (b) NEMO 1/12∘ model potential density, both averaged between 140 and 170∘ E and between 5 and 7∘ N, between 1976 and 2016. In panel (a) the grey line is the 3-month running average of the Nino 3.4 index multiplied by 100. The shading highlights 1982, 1997 and 2015 when strong eastern Pacific El Niños developed. In panel (b) the contours, interpolated from the density at model levels, are at integer values of density. The densest contour is at 1027 kgm-3.
Figure plots the integrated Ekman pumping, the uplift, for the 40-year period between 1976 and 2016. As before, this is calculated from
the DRAKKAR dataset for the region between 140 and 170∘ E and between 5 and 7∘ N.
Also shown is the average potential density in the region from the NEMO 1/12∘ model and a scaled value of the Nino 3.4 index
covering the same period.
The figure shows the significant amount of uplift that occurred prior to the strong El Niños of 1982–1983, 1997–1998 and 2015–2016 and the
corresponding rise in the ocean density surfaces. It also shows that similar amounts of uplift in 1997 and 2002 may be responsible for the shallowing
of the density surfaces and the increased El Niño index.
Between 1990 and 1996 the figure shows a steady increase in uplift, with a mean slope only slightly lower than that seen prior to the strong EP
El Niño years. However the El Niño index remains low, except for a period near the end of 1992 and a smaller event in 1995. Ekman pumping was also
active early in the year in 1987 and 2002 but no strong EP El Niño developed.
The figure also shows that in late 2009, the El Niño index reached 1.6 and there was a corresponding shallowing of the density surfaces, but the
estimated uplift is small.
Another feature of the figure is that as well as the period of positive Ekman uplift between 1990 and 1996, it shows significant drops in the
integrated pumping between 1997 and 2002 and between 2007 and 2014. It is possible that the change in sign of the decadal average slope around 1996 is
due to changes in the data available or the methods used by the ECMWF reanalysis project. However, there is a small correlation between the sign of the
El Niño index and the slope of the uplift curve, so the change before and after 1996 may be real.
Overall, the results imply that, although Ekman pumping may be required to lower sea level in the North Equatorial Trough before a strong EP El Niño
can develop, other factors are also likely to be involved. If this is correct it implies that Ekman pumping early in the year may be a necessary but
not sufficient condition for the development of EP El Niños.
One obvious additional factor is that the volume of warm water available in the West Pacific Warm Pool needs to be sufficient. Although the
differences in depths are small, this may help explain why the strong El Niño developed from years when the 1023 kgm-3 density surface
lies deeper in the ocean, as opposed to the 1990 to 1995 period when the average depth was shallower.
Conclusions
The first aim of this study was to understand the cause of the low sea levels that developed in the western Pacific, along the line of the North Equatorial Trough, during the growth of the strong 1982–1983 EP El Niño.
A comparison of results from the NEMO 1/12∘ and Occam 1/4∘ global ocean models indicated that the latter was suitable for studying
the development of the 1982–1983 event. The Occam model was then used in a series of short tests, using starting dates in the same 2 years, to determine whether the state of the ocean or the wind field was responsible for the low sea levels.
The results showed that the key feature was the local wind field in the western equatorial Pacific and that, for these years, differences in the
initial state of the ocean were not significant.
The most likely way that the local wind field can affect sea level is through Ekman pumping, which affects the ocean's stratification. The hypothesis
that this was the cause of the sea level drop in the North Equatorial Trough was checked by using data from the original run of the
NEMO 1/12∘ model to compare Ekman pumping in the western section of the trough with the changes in depth of the model's density surfaces.
The comparison showed that during the first few months of the year, the rise in density surfaces was close to that expected from Ekman pumping but
that after mid-year, as the warm surface layer became thinner, the connection was broken, although Ekman pumping continued.
The comparison was then repeated for the period prior to the strong EP El Niños of 1997–1998 and 2015–2016, and similar results were obtained.
As a further check on the hypothesis, the integrated Ekman pumping in the western region of the trough was calculated for the period 1976 to 2016 and
compared with both the movement of the average density surfaces and the Nino 3.4 index during this period.
This showed that the presence of strong Ekman pumping during the first few months of the year was not sufficient by itself to trigger a strong EP
El Niño. Ekman pumping may be a necessary condition but other factors, possibly including the volume of water in the West Pacific Warm Pool, must also
be involved.
El Niños and westerly winds
Periods of positive El Niño index are often thought to be triggered by periods of westerly winds in the western Pacific associated with the
Madden–Julian Oscillation.
One possible mechanism, connecting the two, requires westerly wind bursts to generate equatorial Kelvin waves which eventually result in a warmer Equatorial Pacific Cold Pool. However, there are doubts about this mechanism because the Equatorial Pacific Cold Pool region is never warm enough to trigger deep atmospheric convection.
However, the westerly winds can also reverse the direction of the Equatorial Current in the western Pacific, and so they transport Warm Pool water into the
central Pacific. The temperature of the advected water is sufficient to trigger deep atmospheric convection and its effect on the position of the warm
water front in the central Pacific may explain much of the normal year-to-year variation in the El Niño index.
These mechanisms focus on waves and currents near the Equator, but, as discussed earlier, observations show that water warm enough to cause deep
atmospheric convection is carried even further eastwards by the North Equatorial Counter Current.
investigated the cause of the strong 1982–1983 and 1997–1998 El Niños and argued they were due to the NECC advecting a larger than normal volume of Warm Pool water into the central and eastern Pacific as a result of low sea levels in the North Equatorial Trough early in the year.
The study reported in this paper shows that these low sea levels are the result of local winds upwelling water into the surface layer. Although it is
affected by the topography of the adjacent ocean continent, the northward extent of MJOs and other similar wind features in the western Pacific is
close to the scale of the equatorial Rossby radius in the atmosphere and so includes the latitudes of the NECC.
As a result, it is very possible that the mismatch between the equatorial Rossby radii in the ocean and atmosphere means that MJO westerlies may
trigger strong EP Pacific El Niños as a result of Ekman-driven upwelling, at the latitude of the North Equatorial Trough, in addition to their
impact on currents near the Equator.
Code and data availability
At the time of publication, the NEMO 1/12∘ global ocean model archive datasets are freely available at . The NEMO ocean model code and its recent documentation are available online. The Occam ocean model code used in this study is a two-grid development of the “moma.mpi” branch of the Moma ocean model code .
Competing interests
The author was a founding editor of Ocean Science. The peer-review process was guided by an independent editor, the author no longer has any editorial role with Ocean Science and he has no other competing interests to declare.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
This work contributed to and was aided by the research programme of the Marine Systems Modelling group at the UK National Oceanography Centre (NOC), part of the Natural Environment Research Council (NERC). NERC helped fund the investigation through the NOC National Capability funding. Part of the analysis was carried out using the JASMIN Service at the UK Centre for Environmental Data Analysis, also funded by NERC.
The author also wishes to acknowledge the helpful comments and suggestions from Harry Bryden and four anonymous reviewers.
Financial support
This research has been supported by the Natural Environment Research Council (Gold Standard Open Source Publication Fund).
Review statement
This paper was edited by Anne Marie Tréguier and reviewed by four anonymous referees.
ReferencesAshok, K., Behera, S. K., Rao, S. A., Weng, H., and Yamagata, T.: El Niño Modoki and its possible teleconnection, J. Geophys. Res., 112, C11007, 10.1029/2006JC003798, 2007.
Baturin, N. and Niiler, P.:
Effects of instability waves in the mixed layer of the equatorial Pacific,
J. Geophys. Res.,
102, 27771–27793, 1997.
Bjerknes, J.:
Atmospheric Teleconnections from the Equatorial Pacific,
Mon. Weather Rev.,
97, 163–172, 1969.Bourdallé-Badie, R., Bell, M., Chanut, J., Clementi, E., Coward, A., Drudi, M., Éthé, C., Iovino, D., Lea, D., Lévy, C., Madec, G., Martin, N., Masson, S., Mathiot, P., Mocavero, S., Müller, S., Nurser, G., Samson, G., and Storkey, D.: “NEMO ocean engine”, Scientific Notes of Climate Modelling Center, 27 – ISSN 1288–1619, Institut Pierre-Simon Laplace (IPSL), Zenodo [code], 10.5281/zenodo.1464816, 2019.
Bryan, K.:
A numerical method for the study of the circulation of the world ocean,
J. Comput. Phys.,
4, 1255–1273, 1969.
Bryden, H. L. and Brady, E. C.:
Diagnostic Model of the Three-Dimensional Circulation in the Upper Equatorial Pacific Ocean,
J. Phys. Oceanogr.,
15, 1255–1273, 1985.Byrne, M. P., Pendergrass, A. G., Rapp, A. D., and Wodzicki, K. R.:
Response of the Intertropical Convergence Zone to Climate Change: Location, Width, and Strength,
Current Climate Change Reports,
4, 355–370, 10.1007/s40641-018-0110-5, 2018.Cai, W., Borlace, S., Lengaigne, M., van Rensch, P., Collins, M., Vecchi, G., Timmermann, A., Santoso, A., McPhaden, M. J., Wu, L., England, M. H., Wang, G., Guilyardi, E., and Jin, F.-F.:
Increasing frequency of extreme El Niño events due to greenhouse warming,
Nature Clim. Change,
4, 111–116, 10.1038/nclimate2100, 2014.
Cane, M. A.:
Oceanographic Events during El Niño,
Science,
222, 1189–1195, 2011.
Chelton, D. B., Wentz, F. J., Genremann, C. L., de Szoeke, R. A., and Schlax, M. G.:
Satellite Microwave SST Observations of Transequatorial Tropical Instability Waves,
Geophys. Res. Lett.,
27, 1239–1242, 2000.
Chen, X., Ling, J., and Li, C.:
Evolution of the Madden–Julian Oscillation in Two Types of El Niño,
J. Climate,
29, 1919–1934, 2016.
Chiodi, A. M. and Harrison, D. E.:
Observed El Niño SSTA Development and the Effects of Easterly and Westerly Wind Events in 2014/15,
J. Climate,
30, 1505–1519, 2017.Clarke, A. J.: El Niño Physics and El Niño Predictibility, Annu. Rev. Mar. Sci., 6, 79–99, 10.1146/annurev-marine-010213-135026, 2014.Coward, A. C.: NEMO 1/12∘ global ocean model archive, JASMIN/CEDA [data set], http://gws-access.ceda.ac.uk/public/nemo/runs/ORCA0083-N06/means (last access: 27 October 2021), 2015.
Cox, M.:
A primitive equation 3-dimensional model of the ocean, Ocean Group Technical Report 1,
Geophysical Fluid Dynamics Laboratory, Princeton University, Princeton, NJ 08542, 1989.
de Cuevas, B., Webb, D., Coward, A., Richmond, C., and Rourke, E.:
The UK Ocean Circulation and Climate Advanced Modelling Project (OCCAM),
in: High Performance Computing,
edited by: Allan, R., Guest, M., Simpson, A., Henty, D., and Nicole, D.,
Kluwer Academic, New York, 681 pp., 325–336, 1999.Dussin, R., Barnier, B., Brodeau, L., and Molines, J. M.:
The making of Drakkar forcing set DFS5, Report DRAKKAR/MyOcean Report 01-04-16, Laboratoire de Glaciologie et Géophysique de l'Environnement CNRS, Université de Grenoble, available at: http://archimer.ifremer.fr/doc/00000/11154/14032.pdf
(last access: 25 October 2021), 2016.Evans, J. L. and Webster, C. C.:
A variable sea surface temperature threshold for tropical convection,
Aust. Meteorol. Ocean.,
64, S1–S8, 10.22499/2.6401.007, 2014.
Flato, G., Marotzke, J., Abiodun, B., Braconnot, P., Chou, S., Collins, W., Cox, P., Driouech, F., Emori, S., Eyring, V., Forest, C., Gleckler, P., Guilyardi, E., Jakob, C., Kattsov, V., Reason, C., and Rummukainen, M.:
Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, chap. Evaluation of Climate Models,
Cambridge University Press, Cambridge, UK and New York, NY, USA, 2013.
Gadgil, S., Joseph, P. V., and Joshi, N. V.:
Ocean-atmosphere coupling over monsoon regions,
Nature,
312, 141–143, 1984.Griffies, S. M., Gnanadesikan, A., Dixon, K. W., Dunne, J. P., Gerdes, R., Harrison, M. J., Rosati, A., Russell, J. L., Samuels, B. L., Spelman, M. J., Winton, M., and Zhang, R.: Formulation of an ocean model for global climate simulations, Ocean Sci., 1, 45–79, 10.5194/os-1-45-2005, 2005.
Guilyardi, E., Wittenberg, A., Fedorov, A., Collins, M., Wang, C., Capotondi, A., and van Oldenborg, G. J.:
Understanding El Niño in Ocean–Atmosphere General Circulation Models: Progress and Challenges,
B. Am. Meteorol. Soc.,
90, 325–340, 2009.
Ham, Y.-G. and Kug, J.-S.:
How well do current climate models simulate two types of El Nino?,
Clim. Dynam.,
39, 383–398, 2012.
Hansen, D. V. and Paul, C. A.:
Genesis and Effects of Long waves in the Equatorial Pacific,
J. Geophys. Res.,
89, 10431–10440, 1984.
Harrison, D., Giese, B., and Sarachij, E.:
Mechanisms of SST Change in the Equatorial waveguide during the 1982–83 ENSO,
J. Climate,
3, 173–188, 1990.Hsin, Y.-C. and Qiu, B.:
Seasonal fluctuations of the surface North Equatorial Countercurrent (NECC) across the Pacific basin,
J. Geophys. Res.,
117, 1–17, 10.1029/2011JC007794, 2012.Hsu, C.-W., Yin, J., Griffies, S. M., and Dussin, R.: A mechanistic analysis of tropical Pacific dynamic sea level in GFDL-OM4 under OMIP-I and OMIP-II forcings, Geosci. Model Dev., 14, 2471–2502, 10.5194/gmd-14-2471-2021, 2021.
Hurlburt, H., Kindle, J. C., and O'Brien, J. J.:
A Numerical Simulaiton of the Onset of El Niño,
J. Phys. Oceanogr.,
6, 621–631, 1976.
Jochum, M., Cronin, M. F., Kessler, W. S., and Shea, D.:
Observed horizontal temperature advection by tropical instability waves,
Geophys. Res. Lett.,
34, 1–4, 2007.
Johnson, G. C., Sloyan, B. M., Kessler, W. S., and McTaggart, K. E.:
Direct measurements of upper ocean currents and water properties across the tropical Pacific during the 1990s,
Prog. Oceanogr.,
52, 31–61, 2002.Kashino, Y., Atmadipoera, A., Kuroda, Y., and Lukijanto:
Observed features of the Halmahera and Mindanao Eddies,
J. Geophys. Res.-Oceans,
118, 6543–6560, 10.1002/2013JC009207, 2003.
Kennan, S. C. and Lament, P. J. F.:
Observations of a Tropical Instability Vortex,
J. Phys. Oceanogr.,
30, 2277–2301, 2000.
Knauss, J. A.:
The Structure of the Pacific Equatorial Countercurrent,
J. Geophys. Res.,
66, 143–155, 1961.Larkin, N. K. and Harrison, D. E.: On the definition of El Niño and associated seasonal average U.S. weather anomalies, Geophys. Res. Lett., 32, L13705, 10.1029/2005GL022738, 2005.Lengaigne, M., Boulanger, J.-P. H., Menkes, C., Masson, S., Madec, G., and Delecluse, P.: Ocean response to the March 1997 Westerly Wind Event, J. Geophys. Res., 107, 8015, 10.1029/2001JC000841, 2002.
Love, C.:
EASTROPIC Atlas Vol 1,
Government Printing Office, Washington, DC, 1972.
Love, C.:
EASTROPIC Atlas Vol 9,
Government Printing Office, Washington, DC, 1975.
Madden, R. A. and Julian, P. R.:
Description of a 40–50 day oscillation in the zonal wind in the tropical Pacific,
J. Atmos. Sci.,
28, 702–708, 1971.
Madden, R. A. and Julian, P. R.:
Description of global-scale circulation cells in the tropics with a 40–50 day period,
J. Atmos. Sci.,
29, 1109–1123, 1972.
Madden, R. A. and Julian, P. R.: Historical perspective, in: Intraseasonal Variability in the Atmosphere–Ocean Climate System, 2nd Edn., edited by: Lau, W. K. M. and Waliser, D. E., Springer, Heidelberg, Dordrecht, London, New York Praxis Publishing, Chichester, UK, 1–16, 2012.
Madec, G., Delecluse, P., Imbard, M., and Levy, C.: OPA 8.1 Ocean General Circulation Model reference manual, Note du Pole de modelisation XX, Institut Pierre-Simon Laplace (IPSL), France, 1998.
McCreary, J.:
Eastern Tropical response to Changing Wind Systems: with applicaiton to El Niño,
J. Phys. Oceanogr.,
6, 632–645, 1976.
McCreary, J. P.:
Modelling Equatorial Ocean Circulation,
Annu. Rev. Fluid Mech.,
17, 359–409, 1985.
McPhaden, M. J.:
Continuously Stratified Models of the Steady-State Equatorial Ocean,
J. Phys. Oceanogr.,
11, 337–354, 1981.
McPhaden, M. J.:
Trade Wind Fetch Related Variations in Equatorial Undercurrent Depth, Speed, and Transport,
J. Geophys. Res.,
98, 2555–2559, 1993.
Menkes, C., Vialard, J., Kennan, S., Boulanger, J.-P., and Madec, G.:
A Modeling Study of the Impact of Tropical Instability Waves on the Heat Budget of the Eastern Equatorial Pacific,
J. Phys. Oceanogr.,
36, 847–865, 2006.
Meyers, G.:
On the Annual Rossby Wave in the Tropical North Pacific Ocean,
J. Phys. Oceanogr.,
9, 663–674, 1979.
Meyers, G. and Donguy, J.-R.:
The North Equatorial Countercurrent and heat storage in the western Pacific Ocean during 1982–83,
Nature,
312, 258–268, 1984.
Min, Q., Su, J., Zhang, R., and Rong, X.:
What hindered the El Niño pattern in 2014,
Geophys. Res. Lett.,
42, 6762–6770, 2015.
Montgomery, R. and Palmen, E.: Contribution to the question of the Equatorial Counter Current, J. Mar. Res., 4, 112–133, 1940.
Munk, W. H.:
On the wind driven ocean circulation,
J. Meteorol.,
7, 79–93, 1950.NEMO Consortium: NEMO/XIOS 2.5, IPSL [code], https://forge.ipsl.jussieu.fr/nemo/svn/NEMO/releases/r4.0/r4.0.6, last access: 26 October 2021.
Neumann, G. and Pierson, W. J.: Principals of Physical Oceanography, Prentice-Hall, Inc., Englewood Cliffs, New Jersey, USA, 1966.NOAA: Cold & Warm Episodes by Season, NOAA Climate Prediction Center, available at: http://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.php (last access: 26 October 2021), 2021.
Pacanowski, R. and Philander, S.:
Parameterization of vertical mixing in numerical models of tropical oceans,
J. Phys. Oceanogr.,
11, 1443–1451, 1981.
Philander, S.:
El Niño and La Niña,
J. Atmos. Sci.,
42, 2652–2662, 1985.
Philander, S.:
El Niño and La Niña,
Am. Sci.,
77, 451–459, 1989.Picaut, J., Ioualalen, M., Menkes, C., Delcroix, T., and McPhaden, M. J.:
Mechanism of the Zonal Displacements of the Pacific Warm Pool: Implications for ENSO,
Science,
274, 1486–1489, 10.1126/science.274.5292.1486, 1996.
Semtner, A.:
A general circulation model for the World Ocean, Technical Report No. 9,
Department of Meteorology, University of California, Los Angeles, 1974.
Smith, R., Dukowicz, J., and Malone, R.:
Parallel ocean general circulation modelling,
Physica D,
60, 39–61, 1992.Sun, Z., Liu, H., Lin, P., Tseng, Y., Small, J., and Bryan, F.:
The Modeling of the North Equatorial Countercurrent in the Community Earth System Model and its Oceanic Component,
J. Adv. Model. Earth Sy.,
1, 531–544, 10.1029/2018MS001521, 2019.
Sverdrup, H. U.:
Wind-driven currents in a baroclinic ocean; with application to the equatorial currents of the eastern Pacific,
P- Natl. Acad. Sci. USA,
33, 318–326, 1947.
Taft, B. A. and Kovala, B.: Vertical sections of temperature, salinity, thermosteric anomaly and zonal geostrophic velocity from NORPAX Shuttle Experimant – Part 1, NOAA Data Report ERL PMEL-3, NTIS: PB82-163106, PMEL –
Pacific Marine Environmental Laboratory
Seattle, Washington, USA, 98 pp., 1981.
Taft, B. A., Cantos-Figuerola, A., and Kovala, P. Vertical sections of temperature, salinity, thermosteric anomaly and zonal geostrophic velocity from NORPAX Shuttle Experimant – Part 3, NOAA Data Report ERL PMEL-7, NTIS: PB83-169961, PMEL – Pacific Marine Environmental Laboratory Seattle, Washington, USA, 92 pp., 1982.Tan, S. and Zhou, H.:
The observed impacts of the two types of El Niño on the North Equatorial Countercurrent in the Pacific Ocean,
Geophys. Res. Lett.,
45, 10439–10500, 10.1029/2018GL079273, 2018.
Thompson, S. R.: Sills of the Global Ocean: a compilation, Ocean Model., 109, 7–9, 1995.Wang, G. and Hendon, H. H.:
Why 2015 was a strong El Niño and 2014 was not,
Geophys. Res. Lett.,
44, 8567–8575, 10.1002/2017GL074244, 2017.
Webb, D.:
The vertical advection of momentum in Bryan-Cox-Semtner ocean general circulation models,
J. Phys. Oceanogr.,
25, 3186–3195, 1995.Webb, D., Coward, A., de Cuevas, B., and Gwilliam, C.:
A Multiprocessor Ocean General Circulation Model Using Message Passing,
J. Atmos. Ocean. Tech.,
14, 175–183, 10.1175/1520-0426(1997)014<0175:AMOGCM>2.0.CO;2, 1997.
Webb, D., de Cuevas, B., and Richmond, C.:
Improved advection schemes for ocean models,
J. Atmos. Ocean. Tech.,
15, 1171–1187, 1998.Webb, D. J.: On the role of the North Equatorial Counter Current during a strong El Niño, Ocean Sci., 14, 633–660, 10.5194/os-14-633-2018, 2018.Webb, D. J.: Moma, GitHub [code], https://github.com/djwebb/moma (last access: 26 October 2021), 2019.Webb, D. J., Coward, A. C., and Snaith, H. M.: A comparison of ocean model data and satellite observations of features affecting the growth of the North Equatorial Counter Current during the strong 1997–1998 El Niño, Ocean Sci., 16, 565–574, 10.5194/os-16-565-2020, 2020.Wijaya, Y. J. and Hisaki, Y.:
Differences in the Reaction of North Equatorial Countercurrent to the Developing and Mature Phase of ENSO Events in the Western Pacific Ocean,
Climate,
9, 175–183, 10.3390/cli9040057, 2021.
Wyrtki, K.:
Teleconnections in the Equatorial Pacific Ocean,
Science,
180, 66–68, 1973.
Wyrtki, K.:
Equatorial Currents in the Pacific 1950 to 1970 and Their Relation to the Trade Winds,
J. Phys. Oceanogr.,
4, 372–380, 1974a.Wyrtki, K.:
Sea Level and the Seasonal Fluctuations of the Equatorial Currents in the Western Pacific Ocean,
J. Phys. Oceanogr.,
4, 91–103, 1974b.
Wyrtki, K. and Kilonsly, B.:
Mean Water and Current Structure during the Hawaii-to-Tahiti Shuttle Experiment,
J. Phys. Oceanogr.,
14, 242–254, 1984.
Wyrtki, K., Stroup, E., Patzert, W., Williams, R., and Quinn, W.:
Predicting and Observing El Niño,
Science,
191, 343–346, 1976.Yu, J.-Y. and Liu, W. T.:
A linear relationship between ENSO intensity and tropical instability wave activity in the eastern Pacific Ocean,
Geophys. Res. Lett.,
30, CLM 5-1–5-5, 10.1029/2003GL017176, 2003.
Yu, Z., Mccreary, J. P., Kessler, W. S., and Kelly, K. A.:
Influence of Equatorial Dynamics on the Pacific North Equatorial Countercurrent,
J. Phys. Oceanogr.,
30, 3179–3190, 2000.
Zhang, C.:
Large-Scale Variability of Atmospheric Deep Convection in Relation to Sea Surface Temperature in the Tropics,
J. Climate,
6, 1898–1913, 1993.Zhang, C.:
Madden–Julian Oscillation,
Rev. Geophys.,
43, 1–36, 10.1029/2004RG000158, 2005.
Zhang, R.-H. and Busalacchi, A. J.:
A Possible Link between Off-equatorial Warm Anomalies Propagating along the NECC Path and the Onset of the 1997–98 El Niño,
Geophys. Res. Lett.,
26, 2873–2876, 1999.Zhao, J., Li, Y., and Wang, F.:
Dynamical responses of the west Pacific North Equatorial Countercurrent (NECC) system to El Niño events,
J. Geophys. Res.-Oceans,
118, 2828–2844, 10.1002/jgrc.20196, 2013.