Journal of Oceanology and Limnology   2020, Vol. 38 issue(4): 1123-1137     PDF       
http://dx.doi.org/10.1007/s00343-020-9329-9
Institute of Oceanology, Chinese Academy of Sciences
0

Article Information

WANG Jing, YUAN Dongliang, XU Tengfei, ZHAO Xia
Wind-forced equatorial wave dynamics of the Pacific Ocean during 2014/2015 and 2015/2016 El Niño events
Journal of Oceanology and Limnology, 38(4): 1123-1137
http://dx.doi.org/10.1007/s00343-020-9329-9

Article History

Received Dec. 20, 2019
accepted in principle Mar. 16, 2020
accepted for publication Apr. 24, 2020
Wind-forced equatorial wave dynamics of the Pacific Ocean during 2014/2015 and 2015/2016 El Niño events
WANG Jing1,2,3, YUAN Dongliang1,2,3, XU Tengfei4, ZHAO Xia1,2,3     
1 CAS Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China;
2 Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266071, China;
3 Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China;
4 First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China
Abstract: The equatorial wave dynamics of interannual sea level variations between 2014/2015 and 2015/2016 El Niño events are compared using the Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics Climate Ocean Model (LICOM) forced by the National Centers for Environmental Prediction (NCEP) reanalysis Ⅰ wind stress and heat flux during 2000-2015. In addition, the LICOM can reproduce the interannual variability of sea surface temperature anomalies (SSTA) and sea level anomalies (SLA) along the equator over the Pacific Ocean in comparison with the Hadley center and altimetric data well. We extracted the equatorial wave coefficients of LICOM simulation to get the contribution to SLA by multiplying the meridional wave structure. During 2014/2015 El Niño event, upwelling equatorial Kelvin waves from the western boundary in April 2014 reach the eastern Pacific Ocean, which weakened SLA in the eastern Pacific Ocean. However, no upwelling equatorial Kelvin waves from the western boundary of the Pacific Ocean could reach the eastern boundary during the 2015/2016 El Niño event. Linear wave model results also demonstrate that upwelling equatorial Kelvin waves in both 2014/2015 and 2015/2016 from the western boundary can reach the eastern boundary. However, the contribution from stronger westerly anomalies forced downwelling equatorial Kelvin waves overwhelmed that from the upwelling equatorial Kelvin waves from the western boundary in 2015. Therefore, the western boundary reflection and weak westerly wind burst inhibited the growth of the 2014/2015 El Niño event. The disclosed equatorial wave dynamics are important to the simulation and prediction of ENSO events in future studies.
Keywords: Pacific Ocean    equatorial waves    western boundary reflection    El Ni?o    
1 INTRODUCTION

El Niño and La Niña events are the strongest anomalous climate events in the tropical Pacific Ocean. El Niño (La Niña) events are warm (cold) phases of the Pacific El Niño-Southern Oscillation (ENSO), which have global repercussions even influences in the Indian Ocean. Figure 1 shows the Niño 3 index (shading) since 2000 using Optimum Interpolation Sea Surface Temperature (OISST) data based on the period of 1982–2015. The anomaly is calculated relative to a monthly climatological seasonal cycle based on the years 1982–2005. Strong El Niño events have been identified by the commonly used definitions in years 1982/1983, 1986/1987, 1997/1998, 2002/2003, 2004/2005, 2006/2007, 2009/2010, 2015/2016 (Santoso et al., 2017). The Scientific community was abuzz that a strong El Niño event did not evolve in 2014 despite westerly wind bursts (WWBs) occurred west of the dateline in early 2014 (Menkes et al, 2014; McPhaden, 2015). In the conventional description of ENSO, WWBs trigger an El Niño event with multiple WWBs initiating the largest El Niño events (McPhaden, 1999; Yu et al., 2013). Wind stress anomalies determine the magnitude and timing of El Niño events (Penland and Sardeshmukh, 1995; Penland, 1996). Westerly anomalies forced downwelling equatorial Kelvin waves play an important role in the coupled process of the equatorial Pacific Ocean and the atmospheric condition associated with ENSO (Lengaigne et al., 2004). The latest studies have demonstrated easterly wind bursts in the boreal summer inhibit the development of an El Niño event (Menkes et al., 2014; Hu and Fedorov, 2016) and also push the 2015– 2016 El Niño to an extreme event (Levine and McPhaden, 2016). The development of the 2014/2015 and 2015/2016 El Niño events forms an important comparison in ENSO dynamics, which has drawn widespread attention recently. Equatorial waves play an important role in the conceptual ENSO oscillators including the delayed oscillator (Schopf and Suarez, 1988; Battisti and Hirst, 1989), advective-reflective oscillator (Picaut et al., 1997) and Western Pacific oscillator (Weisberg and Wang, 1997). However, the role of equatorial wave dynamics during the 2014/2015 and 2015/2016 El Niño events has not been studied before. We will use the 2014/2015 El Niño for convenience in the following text although the 2014 El Niño did not evolve to a strong event.

Fig.1 Comparison of simulated monthly SSTA time series (thick dash) in the Niño 3 region with the OISST data (shading) 2000–2015 The standard deviation is marked in dotted line.

Both Western boundary and eastern boundary reflections play an important role in ENSO cycles in the Pacific Ocean (Gent and Luyten, 1985; Rothstein et al., 1985; Schopf and Suarez, 1988; Kessler and McCreary, 1993) and even in the evolution of the Indian Ocean Dipole in the Indian Ocean (Yuan and Liu, 2009; Wang and Yuan, 2015). There are still uncertainties about the wind-forced and westernboundary-forced equatorial Kelvin waves from the western Pacific Ocean (Boulanger and Menkes, 1999; Delcroix et al., 2000) associated with the approximation made to derive the wave coefficients from altimetry. Some studies ignore the nonlinearities using the linear model of the interannual equatorial waves (McPhaden and Yu, 1999; Delcroix et al., 2000). Yuan et al. (2004) studied the dynamics of ENSO oscillator paradigm in the 1990s using improved wave decomposition method and nonlinear general ocean circulation model (OGCM). However, a tropical Pacific Ocean model with closed western boundary and lateral boundaries is used in their study. Therefore, the process of the western boundary reflection in the Pacific-only OGCM is not realistic. In our study, a nonlinear quasi-global model with real topographic boundary and complete dynamics of the ocean circulation is used to study the equatorial wave dynamics during 2014/2015 and 2015/2016 El Niño events. Therefore, the analyses of equatorial waves based on the nonlinear OGCM are likely reliable in the Pacific Ocean.

In this paper, we offer a possible explanation for why the 2014/2015 El Niño event did not evolve a strong ENSO event by using a nonlinear OGCM and by focusing on equatorial wave dynamics in comparison with the 2015/2016 El Niño event. In Section 2, we give a description of the data and model, and then present the model results and wave analyses for the evolution of two El Niño events in Section 3. At last, discussion and conclusions are presented in Section 4 and Section 5, respectively.

2 DATA AND METHOD

The Niño 3 index data are from the OISST weekly data based on the period 1982-2019. The data can be downloaded from the following website: http://stateoftheocean.osmc.noaa.gov/sur/pac/. World Ocean Atlas 2013 version 2 (WOA13) data are used to validate the model simulated vertical profiles of density and buoyancy frequency. The WOA13 data consist of a set of objectively analyzed climatological fields of in situ temperature and salinity etc. (Boyer et al., 2013). The annual mean density and buoyancy frequency squared are calculated based on WOA13 data. The sea surface temperature (SST) data from the Hadley Center data (HadISST, Rayner et al., 2003) are used to validate the simulated SST. The SST climatology was computed based on the monthly average SST data during 2000-2015. Basic data sets and procedures are described in Kessler (1990) and Kessler and McCreary (1993). The monthly sea level anomalies (SLA) from the Archiving, Validation, and Interpretation of Satellite Oceanographic (AVISO) data with a horizontal resolution of 0.25° longitude by 0.25° latitude from 2000 to 2015 are used in this study. The domain covers the region of (0.125°– 359.875°, 89.875°S–89.875°N). The daily forcing of wind stress and heat flux from National Centers for Environmental Prediction (NCEP) reanalysis Ⅰ data is used as the forcing of the OGCM with horizontal resolution of 2.5°×2.5° from 1948 to 2015.

The OGCM used here is the LICOM (LASG/IAP Climate Ocean Model, version 1.1) developed from Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics (LASG/IAP). The model has a free sea surface and nearly leveled vertical coordinates in the deep ocean. There are total 30 levels in the vertical with 12 layers in the upper 300 m. And their depths locate at around -12.5, -37.5, -62.5, -87.5, -112.5, -137.5, -162.5, -187.5, -212.5, -237.5, -262.5, -287.5, -313.5, -344.6, -387.8, -450.7, -540.3, -662.8, -823.6, -1 026.6, -1 274.6, -1 569.0, -1 909.9, -2 295.8, -2 724.0, -3 190.3, -3 689.6, -4 215.6, -4 761.3, and -5 319.1 m. The model domain covers quasi-global ocean (0°–360°, 78.5°S–90°N) with a horizontal resolution of 0.5° longitude by 0.5° latitude. The enhancing mixing in the surface layer follows the mixing scheme used in the paper of Pacanowski and Philander (1981). The model has been spun up for 900 years using climatological forcing of Levitus data from the rest (Levitus and Boyer, 1994). More details of LICOM configuration have been described in Liu et al.(2004a, b). A control run is conducted using the daily forcing of wind stress and heat flux from NCEP reanalysis Ⅰ data from 1948 to 2015. Sea surface temperature and salinity are relaxed to the Levitus climatology data at a time scale of a half-year (Levitus et al., 1998). The analysis of this paper is based on monthly averaged output of the LICOM simulation from 2000–2015. In this study, we also presented the 1997/1998 El Niño event in comparison with the 2014/2015 and 2015/2016 El Niño events in the discussion section. Here, the climatology is calculated based on the monthly averaged output of the LICOM simulation from 1990–2001.

Based on the orthogonal relations of the long equatorial Kelvin and Rossby mode functions, the pressure and zonal currents are projected onto the vertical modes to extract the vertical mode amplitudes, which can be decomposed into the meridional modes. The projection includes the nonlinear effects of the LICOM simulation in the calculation of equatorial wave coefficients by treating the nonlinear terms of the LICOM momentum equations as forcing terms. The pressure and zonal current anomalies from LICOM simulation are decomposed into equatorial waves based on the theory of linear equatorial waves (McCreary, 1981). The Kelvin wave speed was 2.71 m/s for the first baroclinic mode based on the WOA13 data, which is consistent with the calculation of Chelton et al. (1998). The wave decomposition method used here is the same as the method in appendix A in the paper of Yuan et al. (2004).

In order to disentangle the contribution of the wind-forced and reflected wave contribution, we used a simple linear Kelvin wave characteristic-line model to quantify their relative roles. The zonal wind stress anomalies from NCEP reanalysis Ⅰ data is projected onto equatorial Kelvin and Rossby wave modes then is integrated along the characteristic lines with zero condition at the western boundary in the Pacific Ocean for wave solutions. Linear wave model has been introduced in previous studies (Yuan and Han, 2006; Wang and Yuan, 2015). The linear wave model results could estimate the relative contribution of the reflections and wind forcing to the equatorial wave coefficients. In the OGCM output, the contributions from the western boundary reflections and wind forcing cannot be separated quantitatively. The monthly mean zonal wind stress from the NCEP reanalysis Ⅰ daily wind stress data is projected onto the equatorial Kelvin wave and Rossby wave modes, then is integrated along the characteristic lines of equatorial Kelvin waves to obtain the wave coefficients. The eigenvalues and eigen-functions are calculated from the mean vertical profiles of density over the equatorial Pacific Ocean between 5°S and 5°N using WOA13 data. In this paper, we artificially choose the 151°E as the western boundary in the Pacific Ocean since the western boundary reflection is very complicated with one or more narrow gaps at the western boundary (Spall and Pedlosky, 2005).

3 RESULT

In this section, the interannual variability of SLA and SST anomalies simulated by LICOM is validated with the AVISO altimeter SLA data and Hadley Center data, respectively. In addition, the equatorial wave coefficients are extracted from LICOM simulation to analyze the influence of the propagation and reflections of equatorial waves on the evolution of the El Niño events during 2014–2015 and compared to the solution of linear wave model.

3.1 Model validation

El Niño events occurred during the years 2002/2003, 2004/2005, 2006/2007, 2009/2010, and 2015/2016, and La Niña events occur in 2005/2006, 2007/2008, and 2010/2011. These events are simulated well by the LICOM based on the Niño3 index from the OISST data (Fig. 1). The correlation and root-mean-square error (RMSE) between the LICOM simulation and OISST are 0.87 and 0.4. Therefore, it is reasonable to use the LICOM simulation for analyses of ENSO events since this model can basically reproduce the individual ENSO events. Figure 2 shows interannual variability of SST anomalies simulated by the LICOM in comparison with the Hadley Center data averaged between 5°S and 5°N along the equator during 2000–2015. LICOM can reproduce the interannual variability of SST anomalies during ENSO events fairly well. Figure 3 shows the interannual variability of SLA from the LICOM simulation compared with AVISO altimeter data. The interannual anomalies from LICOM are based on monthly climatological averages from 2000 to 2015. Model can reproduce the interannual variability of SLA well in comparison with AVISO altimeter data except that slightly weaker SLA is simulated. In addition, LICOM reproduces fairly well the downwelling anomalies (positive SLA) in the eastern Pacific Ocean and upwelling anomalies (negative SLA) in the western Pacific Ocean during El Niño events. The correlation and RMSE between the simulated SLA and AVISO altimeter data in the Niño3 region are 0.88 and 3.6 cm, respectively. Therefore, the LICOM can be used to study the dynamics of the equatorial Pacific Ocean circulation during 2014–2016 El Niño events, which is also used for equatorial waves analysis in the study of Zhao et al. (2019).

Fig.2 Time-longitude plots of monthly SST anomalies from Hadley centre data (a) and the LICOM simulation (b); averaged between 5°S–5°N along the equator during 2000–2015 (unit: ℃)
Fig.3 Time-longitude plots of monthly SLA (unit: cm) averaged between 5°S–5°N along the equator during 2000-2015 from the AVISO altimeter data (a); the LICOM simulation (b); and time-longitude plot of monthly zonal wind stress anomalies (unit: Pa) (c)

Both model and observations show that the SLA and SST anomalies in the eastern Pacific Ocean in 2015 are much stronger than those in 2014 before mature phase of El Niño event, suggesting that a stronger El Niño event occurred in 2015. Stronger downwelling anomalies occurred in early 2014 (January–April) in comparison with SLA in early 2015 since upwelling anomalies dominated in the eastern equatorial Pacific Ocean in early 2015, which is consistent with the structure of zonal wind anomalies as shown in Fig. 3c. The wind stress anomalies show that Westerly anomalies occurred in early 2014 and easterlies occurred in early 2015 in the western Pacific Ocean. However, stronger westerly anomalies did not result in a stronger El Niño event because the stronger WWBs activity lasted since March 2015, suggesting that the stronger El Niño event in 2015 may be related to the continuous WWBs in comparison with the wind stress anomalies in 2014.

3.2 Decomposed equatorial waves

Figure 4 shows the decomposed equatorial Kelvin wave and the first meridional mode Rossby wave coefficients of the first baroclinic mode during 2000– 2015 using the wave decomposition method (Yuan et al., 2004). Here, we ignore the higher order baroclinic mode results since the second baroclinic mode contribution is energetic in the central and eastern equatorial Pacific Ocean in the mature warming of El Niño events (Dewitte et al., 1999, 2003). In addition, we are focusing the role of the western boundary reflection in the western Pacific Ocean. The decomposed wave coefficients are extracted from the LICOM output to study the wave dynamics of the SLA evolution during the 2014/2015 and 2015/2016 El Niño events. The evolution of the equatorial wave coefficients corresponds to that of SLA in Fig. 3, suggesting the eastward propagation of downwelling equatorial Kelvin waves during El Niño events and of the upwelling equatorial Kelvin waves during La Niña events. The correlation between SLA and wave coefficients (Kelvin wave, Rossby wave) is over 0.8 in the western and eastern Pacific Ocean and slightly low in the central Pacific Ocean. Besides, the correlation between SLA and Rossby wave coefficients is higher than that between SLA and Kelvin wave coefficients.These suggest that the Kelvin wave from the reflection of Rossby wave at the western boundary may play an important role in the evolution of SLA.

Fig.4 Decomposed coefficients of the equatorial Kelvin wave (a) and the Rossby wave (b) of the first baroclinic mode from the LICOM simulation The decomposed coefficients have been smoothed with a 1-2-1 filter and the contour unit corresponds to 0.33 m for Kelvin waves and 0.14 m for the first meridional mode Rossby wave sea level on the equator. The wave coefficients have been detrended.

At the end of the El Niño years such as 2002, 2004, 2006, 2009, 2014, and 2015, upwelling equatorial Kelvin waves from the western boundary propagate eastward, following the propagation of downwelling equatorial Kelvin waves. The arrival of these upwelling equatorial Kelvin waves is coincident with the reversal of the SLA at the eastern boundary in the Pacific Ocean, suggesting that the downwelling anomalies in the eastern Pacific Ocean at the end of the El Niño years are terminated by these upwelling equatorial Kelvin waves from the western boundary. Upwelling equatorial Rossby waves consistently reach the western boundary at the end of El Niño years including 2002, 2004, 2006, 2009, 2014, and 2015. Based on the equatorial wave theory, upwelling equatorial Rossby waves can reflect into the upwelling equatorial Kelvin waves at the western boundary (Clarke, 1983), suggesting that the western boundary reflection may provide weakening impact during the evolution of the downwelling anomalies in the eastern Pacific Ocean and may inhibit the evolution of El Niño event. This phenomenon can be interpreted as the delayed oscillator theory during ENSO events (Schopf and Suarez, 1988). Figure 5 shows the time series of the equatorial wave coefficients at the western and eastern boundary from the LICOM simulation. The decomposed and reflected wave coefficients are plotted in red and green curves, respectively. Here, the reflected wave coefficients are calculated using the linear equatorial wave theory (Clarke, 1991) from the first and the third meridional mode equatorial Rossby wave coefficients decomposed from the model output. The difference between the wave coefficients is plotted in blue curve. The first meridional mode Rossby wave multiplied by its reflection coefficients (0.41) is plotted in orange curve. Here, the reflection coefficients are calculated using the linear theoretical value estimated for a straight coastline and low frequency waves (Clarke, 1983). The use of a 0.41 theoretical reflection coefficient appears to be an overestimation of how much energy can be expected to reflect off the western boundary as suggested by Spall and Pedlosky (2005).Therefore, the role of western boundary reflection may be overestimated. However, large differences can be noticed in Fig. 5a; therefore, we ignored the discussion of the sensitivity of conclusions when using less efficient reflection coefficients. One reason for the differences may be the estimation of the reflected waves from the low-order baroclinic mode at the western boundary in the Pacific Ocean, while another may be from the strong nonlinearity that can be influenced by the complex structure of the nonlinear western boundary currents (Wang and Yuan, 2012, 2014). The reflection process of the Rossby waves at the western boundary is more complicated than the reflection of equatorial Kelvin wave at the eastern boundary because of the strong nonlinearity near the western boundary. Upwelling equatorial Kelvin waves, which propagate eastward and weaken the positive anomalies of SLA in the eastern Pacific Ocean, are detected in early 2014 and 2015 (Fig. 5a, in red). At the eastern boundary, the Rossby waves are in good agreement with the linear reflection of the equatorial Kelvin waves (Fig. 5b). Upwelling equatorial Rossby waves from the reflection of the Kelvin waves can be noticed in June– July 2014 (Fig. 5b). However, only the downwelling equatorial Rossby waves can be observed after May 2015. The difference between 2014 and 2015 at the eastern and western boundaries may induce different downwelling anomalous evolution process during the 2014/2015 and 2015/2016 El Niño events.

Fig.5 Decomposed wave coefficients at the western and eastern boundary a. time series of decomposed (red) and linear reflected (green) Kelvin waves of the first baroclinic mode at the western boundary. The first meridional mode Rossby waves multiply by its reflection ratio is plotted in orange. The difference between decomposed and linear reflected Kelvin wave coefficients is plotted in blue; b. the red and green curves indicate decomposed and linearly reflected Rossby wave coefficients from equatorial Kelvin waves at the eastern boundary. The difference between them is plotted in blue curve.
3.3 Interannual variations of SLA and equatorial wave dynamics during the 2014/2015 and 2015/2016 El Niño events

Figure 6 shows the time-longitude plots of decomposed Kelvin wave coefficients and the first meridional mode Rossby wave coefficients of the first baroclinic mode during 2013–2015. Downwelling Kelvin wave coefficients are seen to propagate eastward, with stronger coefficients seen in 2015 than those in 2014, suggesting that stronger westerly anomalies occur in 2015. This finding is consistent with the Hovmӧller plot of zonal wind stress anomalies along the equator (Fig. 3c). Upwelling equatorial Kelvin waves from the western boundary propagate to the eastern boundary in April 2014, terminating the dowelling Kelvin waves and inhibiting the downwelling anomalies in the eastern Pacific Ocean. However, no such equatorial Kelvin waves could propagate to the eastern boundary in the Pacific Ocean in summer 2015, suggesting that the upwelling equatorial Kelvin waves from western boundary reflection in early 2014 may have played an important role in weakening the downwelling anomalies in the eastern Pacific Ocean.

Fig.6 Time-longitude plots of decomposed wave coefficients, windstress anomalies and Niño 3 index a. d ecomposed Kelvin wave coefficients and b. the first meridional mode Rossby wave coefficients of the first baroclinic mode from LICOM simulation during 2013-2015. The wave coefficient anomalies have been detrended with monthly climatology of 2000-2015; c. zonal wind stress anomalies during 2013-2015; d. Niño 3 index during 2013-2015. A 1-2-1 filter is applied in the plot of wave coefficient.

The evolution of downwelling anomalies instead of the sea surface temperature in the eastern Pacific Ocean is discussed for comparing the 2014/2015 with 2015/2016 El Niño events. For the 2015/2016 El Niño event, downwelling anomalies in the eastern Pacific Ocean occurred in boreal summer, developed in fall and peak in boreal winter. For 2014/2015 El Niño event, downwelling anomalies develop in early 2014 and decrease in July 2014. Although the downwelling anomalies peak in winter, SLA are largely weaker than those in 2015. To analyze the influence of equatorial waves on the evolution of SLA, SLA and wind stress anomalies in the tropical Pacific Ocean are compared from May to August in 2014 and 2015 (Fig. 7).

Fig.7 SLA (unit: cm) and wind stress anomalies (unit: Pa) from May to August in 2014 and 2015

Early 2014 experienced less WWBs activity than early 2015 from the view of longer period, which can also be seen in Fig. 3c. In general, it will take two and a half months for the upwelling equatorial Kelvin waves to propagate to the eastern boundary of the Pacific Ocean based on the equatorial Kelvin wave speed of the first baroclinic mode calculated by Chelton et al. (1998). Upwelling equatorial Kelvin waves generated at the western boundary in April 2014 reached the eastern boundary before July 2014 (Fig. 6). SLA are reduced by at least 5 cm in July 2014 in the eastern Pacific Ocean in comparison with that in June 2014. It is evident that the decreasing of SLA magnitudes is closely related to the arrival of upwelling equatorial Kelvin waves. In contrast, SLA in July 2015 did not decrease in comparison with those in June since no weakening from upwelling equatorial Kelvin waves from the western Pacific Ocean in 2015 (Fig. 6a). In addition, the evolution of SLA as indicated by the AVISO altimeter data is also consistent with the LICOM simulation (Fig. 3a). A clear decrease in SLA occurred in July 2014. The Niño 3 index from the OISST data also shows a decrease in summer 2014, as shown in Fig. 1, suggesting that the arrival of upwelling equatorial Kelvin wave indeed influenced the downwelling anomalies in the eastern Pacific Ocean. In addition, westerlies are stronger in May–August 2015 than those in 2014 in the western and central equatorial Pacific Ocean (Fig. 7). WWBs-forced downwelling equatorial Kelvin waves propagate to the east and increase the downwelling anomalies (positive SLA) in the eastern Pacific Ocean in 2015. The SLA and westerly anomalies in August 2015 is evidently stronger than that in July 2015 (Fig. 7, right panel). However, no evident increase of SLA in August 2014 compared with that in July 2014 is found in the eastern Pacific Ocean especially near the eastern boundary associated with weak WWBs. Therefore, there are two reasons of inducing a stronger El Niño event in 2015: one is the enhancing effects from the westerly anomalies forced downwelling equatorial Kelvin waves; and the other is the absence of weakening effects of upwelling equatorial Kelvin waves from the western boundary reflection in the western Pacific Ocean.

3.4 Processes of western boundary reflection in the Pacific Ocean

To study the generation of the upwelling equatorial Kelvin waves from the western boundary, SLA along 5°N, the equator and zonal wind stress anomalies at the equator are shown in Fig. 8. The upwelling Rossby waves at 5°N are well established at approximately 130°W in the eastern Pacific Ocean and are associated with the westerlies in Fig. 8c, which propagate westward in the Pacific Ocean. In addition, the upwelling equatorial Kelvin waves are generated at the equator in phase with the arrival of the Rossby waves, which can be seen based on the structure of SLA at the equator (Fig. 8b). This finding suggests that upwelling equatorial Kelvin waves are partially generated by the reflection of the upwelling equatorial Rossby waves. The large difference between the decomposed and linearly reflected equatorial Kelvin wave coefficients at the western boundary (Fig. 5a) suggests that not all of the upwelling Kelvin waves are from the reflection of equatorial Rossby waves. A large negative SLA band is perceived in the latitudinal band of 20°S south of the equatorial Pacific Ocean (Fig. 7, left panel). Therefore, the off-equatorial Rossby waves and strong nonlinearity also play a role in the generation of the upwelling Kelvin waves at the western boundary in early 2014. The nonlinear reflection process at the western Pacific Ocean is complicated, and is beyond the scope of the paper. We leave the details of nonlinear processes for future study.

Fig.8 Time-longitude plots of SLA at 5°N (a), the equator (b), and the wind stress anomalies (c) during 2014–2015

Figure 8 shows a qualitative description of the generation of upwelling Kelvin waves at the western boundary. To study the equatorial upwelling Kelvin waves quantitatively, the wind-forced and westernboundary-forced equatorial Kelvin wave coefficients are compared using a linear wave model. The decomposed coefficients at 151°E are integrated along the characteristic lines for obtaining the net forcing from the western boundary without considering the influence of wind forcing in the Pacific Ocean. Figure 9 shows the time-longitude plots of the wind-forced equatorial Kelvin wave coefficients with zero conditions at the western boundary and western-boundary-forced equatorial Kelvin wave coefficients ignoring the influence of wind forcing and the sum of them. Stronger WWBsforced downwelling equatorial Kelvin waves can be observed in 2015 than in 2014, which is consistent with the amplitudes of WWBs (Fig. 3c). In addition, the wind-forced upwelling equatorial Kelvin waves in 2014 and 2015 could not reach the eastern boundary in the Pacific Ocean. In contrast, western-boundaryforced upwelling equatorial Kelvin waves could arrive at the eastern boundary in 2014 and 2015 (Fig. 9b). In addition, a stronger upwelling Kelvin wave could be seen in 2015 (Fig. 6a), which was related to stronger WWBs-forced upwelling equatorial Rossby waves in the western Pacific Ocean (Fig. 6b). However, the upwelling equatorial Kelvin waves from western boundary could not reach the eastern boundary in 2015 (Fig. 6a), a phenomenon which was associated with the overwhelming WWBs-forced downwelling equatorial Kelvin waves seen in 2014. The sum of western-boundary-forced and windforced equatorial Kelvin wave (Fig. 9c) is largely consistent with the decomposed equatorial wave (Fig. 6a), suggesting that linear wave theory is applicable to reproduce the decomposed wave coefficients. The propagation of upwelling equatorial Kelvin waves from the western boundary weakened the downwelling anomalies (positive SLA) in the eastern Pacific Ocean in 2014. Weak WWBs limit the generation of downwelling equatorial Kelvin waves, suggesting no enhancing effects are provided. Therefore, no strong El Niño event can occur in 2014 in comparison with the 2015/2016 El Niño event. In contrast, continuous WWBs-forced downwelling Kelvin waves constantly generated in 2015 accelerated the downwelling anomalies in the eastern Pacific Ocean, which yielded a strong El Niño event (Levine and McPhaden, 2016).

Fig.9 Time-longitude plots of wind-forced equatorial Kevin wave coefficients setting western boundary wave condition (151°E) to zero values (a); western-boundary-forced equatorial Kelvin wave coefficients (b) The wave coefficients have been smoothed with a 1-2-1 filter. c. the sum of (a) and (b).
4 DISCUSSION

A strong El Niño event occurred in 1997. The spring conditions of the zonal wind anomalies in 1997 are similar to those in 2014. However, the strong WWBs occur in the western and central equatorial Pacific Ocean in May 1997 while no such winds occurred in 2014 (Yuan et al., 2004; Chen et al., 2015). Figure 10 shows the time-longitude plots of decomposed Kelvin wave coefficients (a) and the first meridional mode Rossby wave coefficients (b) of the first baroclinic mode from LICOM simulation during 1997–1999. Stronger downwelling equatorial Kelvin waves associated with WWBs propagate eastward during the 1997/1998 El Niño event. This finding is consistent with the zonal wind anomalies in 1997 (Fig. 10c). Upwelling equatorial Kelvin waves from the western boundary could not reach the eastern Pacific Ocean, which is similar to 2015/2016 El Niño event. WWBs-forced equatorial Kelvin wave coefficients are also calculated as shown in Fig. 11. Upwelling equatorial Kelvin waves from WWBs forcing could not arrive at the eastern boundary in 1997. However, upwelling equatorial Kelvin waves from the western boundary could indeed reach the eastern boundary. This finding suggests that the contribution of the strong WWBs-forced downwelling equatorial Kelvin waves overwhelm that of the upwelling equatorial Kelvin waves from the western boundary. The western boundary reflection process in 1997/1998 El Niño is similar to 2015/2016 El Niño. In contrast, western boundary reflections in summer 2014 play an important role in inhibiting the developing of the 2014/2015 El Niño event. An artificial western boundary is set at 151°E in this paper due to consideration that the process of western boundary reflection is complex in the perforated western boundary in the Pacific Ocean. However, nonlinear processes of the western boundary reflection occur and the linear theory cannot be applied to the west of 151°E in the western Pacific Ocean. Therefore, the role of the easterly wind anomalies in the west of 151°E is easily mistaken as the western boundary when calculating the western-boundary-forced Kelvin wave coefficients. Figure 8c shows the time-longitude plot of zonal wind stress anomalies during 2014– 2015. Weak easterly wind anomalies appeared from March 2014 and began to increase in May 2014 to the west of 150°E. In the meantime, easterly wind anomalies also occur near 165°W and were stronger than those anomalies to the west of 150°E. Therefore, stronger upwelling Kelvin waves are purportedly generated near 165°W than west of 151°E in a linear model. Wind-forced upwelling Kelvin waves near 165°W, which can be seen in Fig. 9a using a linear wave model, could only propagate to 140°W and cannot reach the eastern boundary in 2014. Therefore, upwelling equatorial Kelvin waves forced by the easterly wind anomalies to the west of 150°E could not reach the eastern boundary in both a linear and qualitative view. Energies are assumed to dissipate more during the eastward propagation of upwelling equatorial Kelvin waves far to west of 150°E than those near 165°W. From the view of linear equatorial wave dynamics, our study gives a possible explanation the inhibited growth of the 2014/2015 El Niño events, which differs from explanations provided by previous studies (Menkes et al., 2014; Hu and Fedorov, 2016).

Fig.10 Time-longitude plots of decomposed Kelvin wave coefficients (a) and the first meridional mode Rossby wave coefficients (b) of the first baroclinic mode from LICOM simulation during 1997–1998; zonal wind stress anomalies during 1996–1998 (c); Niño 3 index during 1996–1998 (d) The wave coefficient anomalies have been detrended with monthly climatology of 1990-2001. A 1-2-1 filter is applied in the plot of wave coefficient.
Fig.11 Time-longitude plots of wind-forced equatorial Kevin wave coefficients setting western boundary wave condition (151°E) to zero values (a); western-boundary-forced equatorial Kelvin wave coefficients during 1996–1998 (b); sum of the wind-forced and western-boundary-forced equatorial Kelvin wave coefficients (c) The wave coefficients have been smoothed with a 1-2-1 filter.
5 CONCLUSION

Long equatorial wave dynamics of interannual sea level variations during 2014/2015 and 2015/2016 El Niño events are compared using a LICOM simulation forced by wind stress and heat flux data from the daily NCEP reanalysis Ⅰ data during 1948–2015. The LICOM can reproduce the SST anomalies and SLA along the equator over the equatorial Pacific Ocean in comparison with TAO data and AVISO altimeter SLA data well, and can simulate the interannual variability of the SSTA and SLA well including individual El Niño and La Niña events. The equatorial wave coefficients are extracted from LICOM simulation to study the role of the propagation and reflection of the equatorial waves in the evolution of SLA during El Niño events. Significant differences in equatorial waves are analyzed. The responses of the Pacific Ocean to winds are further explored.

Upwelling equatorial Kelvin waves from the western boundary in April 2014 propagated to the eastern Pacific Ocean, which terminated the downwelling equatorial Kelvin waves in the east and weakened the downwelling anomalies (positive SLA) in the eastern Pacific Ocean. However, no such upwelling equatorial Kelvin waves from the western boundary in the Pacific Ocean can reach the eastern boundary during 1997/1998 and 2015/2016 El Niño events. The role of zonal wind forcing and western boundary reflection are compared during 2014/2015 and 2015/2016 El Niño events by calculating the equatorial Kelvin wave coefficients using a linear characteristic wave model. The results demonstrate that the upwelling equatorial Kelvin waves in 2014 are originated partially from the reflection of equatorial Rossby wave at the western boundary in the Pacific Ocean. In addition, off-equatorial Rossby wave and nonlinearity also play a role. Here, we only provide a qualitative explanation. Although upwelling Kelvin wave associated with western boundary reflection in 2015 was stronger than that in 2014, stronger westerly anomalies-forced downwelling equatorial Kelvin waves overwhelmed the role of the western boundary reflection in 2015 since stronger continuous westerly anomalies occurred in 2015 and only weak westerly anomalies occurred after April 2014. Therefore, both the western boundary reflection and weak WWBs inhibit the growth of the 2014/2015 El Niño event.

The analyses presented in this paper suggest that the propagation and reflection of equatorial waves in the Pacific Ocean play an important role in the interannual variability of the downwelling anomalies in the eastern Pacific Ocean during El Niño events. Western boundary reflections provide a weakening effect to the growth of El Niño event. The comparison between linear reflection and decomposed wave coefficients suggests that upwelling equatorial Kelvin waves generated in April 2014 were not only due to the reflection of equatorial Rossby waves but also due to off-equatorial Rossby waves and nonlinearities. This paper provides a possible mechanism for the inhibition of the growth of an El Niño event, and further studies are needed to study the prediction of interannual variability of the tropical climate in the Pacific Ocean and how WWBs activities affect the maintenance or growth of induced warming anomaly or positive SLA. These studies may have useful forecasting implications. In addition, the processes of western boundary reflection to the west of 151°E in the Pacific Ocean need to be analyzed using the model simulation and observational data since we ignore the complicated western boundary reflection process in this paper for simplicity. The disclosed equatorial wave dynamics are important for the evolution and predictability of ENSO events.

6 DATA AVAILABILITY STATEMENT

All data generated and/or analyzed during this study are available from the corresponding author on request.

7 ACKNOWLEDGMENT

We are grateful to the two anonymous reviewers for their insightful comments and suggestions.

References
Battisti D S, Hirst A C. 1989. Interannual variability in a tropical atmosphere-ocean model: influence of the basic state, ocean geometry and nonlinearity. Journal of the Atmospheric Sciences, 46(12): 1 687-1 712. DOI:10.1175/1520-0469(1989)046<1687:IVIATA>2.0.CO;2
Boulanger J P, Menkes C. 1999. Long equatorial wave reflection in the Pacific Ocean from TOPEX/POSEIDON data during the 1992-1998 period. Climate Dynamics, 15(3): 205-225.
Boyer T P, Antonov J I, Baranova O K, Garcia H E, Johnson D R, Mishonov A V, O'Brien T D, Seidov D, Smolyar I, Zweng M M, Paver C R, Locarnini R A, Reagan J R, Forgy C, Grodsky A, Levitus S. World Ocean Database 2013.U.S. Government Printing Office, Washington, USA. 209p.
Chelton D B, DeSzoeke R A, Schlax M G, El Naggar K, Siwertz N. 1998. Geographical variability of the first baroclinic Rossby radius of deformation. Journal of Physical Oceanography, 28(3): 433-460. DOI:10.1175/1520-0485(1998)028<0433:GVOTFB>2.0.CO;2
Chen D K, Lian T, Fu C B, Cane M A, Tang Y M, Murtugudde R, Song X S, Wu Q Y, Zhou L. 2015. Strong influence of westerly wind bursts on El Niño diversity. Nature Geoscience, 8(5): 339-345. DOI:10.1038/ngeo2399
Clarke A J. 1983. The reflection of equatorial waves from oceanic boundaries. Journal of Physical Oceanography, 13(7): 1 193-1 207. DOI:10.1175/1520-0485(1983)013<1193:TROEWF>2.0.CO;2
Clarke A J. 1991. On the reflection and transmission of lowfrequency energy at the irregular western Pacific Ocean boundary. Journal of Geophysical Research, 96(S01): 3 289-3 305. DOI:10.1029/90JC00985
Delcroix T, Dewitte B, duPenhoat Y, Masia F, Picaut J. 2000. Equatorial waves and warm pool displacements during the 1992-1998 El Niño Southern Oscillation events:Observation and modeling. Journal of Geophysical Research-Oceans, 105(C11): 26 045-26 062. DOI:10.1029/2000JC900113
Dewitte B, Reverdin G, Maes C. 1999. Vertical structure of an OGCM simulation of the Equatorial Pacific in 1985-94. Journal of Physical Oceanography, 29(7): 1 542-1 570. DOI:10.1175/1520-0485(1999)029<1542:VSOAOS>2.0.CO;2
Dewitte B, Illig S, Parent L, DuPenhoat Y, Gourdeau L, Verron J. 2003. Tropical Pacific baroclinic mode contribution and associated long waves for the 1994-1999 period from an assimilation experiment with altimetric data. Journal of Geophysical Research-Oceans, 108(C4).
Gent P R, Luyten J R. 1985. How much energy propagates vertically in the equatorial oceans?. Journal of Physical Oceanography, 15(7): 997-1007. DOI:10.1175/1520-0485(1985)015<0997:HMEPVI>2.0.CO;2
Hu S N, Fedorov A V. 2016. Exceptionally strong easterly wind burst stalling El Niño of 2014. Proceedings of the National Academy of Sciences of the United States of America, 113(8): 2 005-2 010. DOI:10.1073/pnas.1514182113
Kessler W S, McCreary J P. 1993. The annual wind-driven Rossby wave in the subthermocline equatorial Pacific. Journal of Physical Oceanography, 23(6): 1 192-1 207. DOI:10.1175/1520-0485(1993)023<1192:TAWDRW>2.0.CO;2
Kessler W S. 1990. Observations of long Rossby waves in the northern tropical Pacific. Journal of Geophysical Research, 95(C4): 5 183-5 217. DOI:10.1029/JC095iC04p05183
Lengaigne M, Guilyardi E, Boulanger J P, Menkes C, Delecluse P, Inness P, Cole J, Slingo J. 2004. Triggering of El Niño by westerly wind events in a coupled general circulation model. Climate Dynamics, 23(6): 601-620. DOI:10.1007/s00382-004-0457-2
Levine A F Z, McPhaden M J. 2016. How the July 2014 easterly wind burst gave the 2015-2016 El Niño a head start. Geophysical Research Letters, 43(12): 6 503-6 510. DOI:10.1002/2016GL069204
Levitus S, Boyer T P, O'Brien T., Antonov J, Stephens C, Stathopolos L, Johnson D, and Gelfeld R. 1998: World Ocean Database 1998. Vol. 1, Introduction, NOAA Atlas NESDIS 18, National Oceanic and Atmopheric Administration. Silver Spring, MD., 346pp.
Levitus S, Boyer T P. 1994. World Ocean Atlas 1994. Volumn 4: Temperature. NOAA Atlas NESDIS 3. U.S.. Department of Commerce, Washington. 117p.
Liu H L, Yu Y Q, Li W, Zhang X H. 2004a. LASG/IAP Climate System Ocean Model (LICOM1, 0): User Manual.Science Publication. Beijing, China. 107p. (in Chinese)
Liu H L, Zhang X H, Li W, Yu Y Q, Yu R C. 2004b. An eddypermitting oceanic general circulation model and its preliminary evaluation. Advances in Atmospheric Sciences, 21(5): 675-690. DOI:10.1007/BF02916365
McCreary J P. 1981. A linear stratified ocean model of the equatorial undercurrent. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 298(1444): 603-635.
McPhaden M J, Yu X R. 1999. Equatorial waves and the 1997-98 El Niño. Geophysical Research Letters, 26(19): 2 961-2 964. DOI:10.1029/1999GL004901
McPhaden M J. 1999. Genesis and evolution of the 1997-98 El Niño. Science, 283(5404): 950-954. DOI:10.1126/science.283.5404.950
McPhaden M J. 2015. Playing hide and seek with El Niño. Nature Climate Change, 5(9): 791-795. DOI:10.1038/nclimate2775
Menkes C E, Lengaigne M, Vialard J, Puy M, Marchesiello P, Cravatte S, Cambon G. 2014. About the role of Westerly Wind events in the possible development of an El Niño in 2014. Geophysical Research Letters, 41(18): 6 476-6 483. DOI:10.1002/2014GL061186
Pacanowski R C, Philander S G H. 1981. Parameterization of vertical mixing in numerical models of tropical oceans. Journal of Physical Oceanography, 11(11): 1 443-1 451. DOI:10.1175/1520-0485(1981)011<1443:POVMIN>2.0.CO;2
Penland C, Sardeshmukh P D. 1995. The optimal growth of tropical sea surface temperature anomalies. Journal of Climate, 8(8): 1 999-2 024. DOI:10.1175/1520-0442(1995)008<1999:TOGOTS>2.0.CO;2
Penland C. 1996. A stochastic model of IndoPacific sea surface temperature anomalies. Physica D: Nonlinear Phenomena, 98(2-4): 534-558. DOI:10.1016/0167-2789(96)00124-8
Picaut J, Masia F, du Penhoat Y. 1997. An advective-reflective conceptual model for the oscillatory nature of the ENSO. Science, 277(5326): 663-666. DOI:10.1126/science.277.5326.663
Rayner N A, Parker D E, Horton E B, Folland C K, Alexander L V, Rowell D P, Kent E C, Kaplan A. 2003. Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century. Journal of Geophysical Research, 108(D14): 4 407. DOI:10.1029/2002JD002670
Rothstein L M, Moore D W, McCreary J P. 1985. Interior reflections of a periodically forced equatorial Kelvin wave. Journal of Physical Oceanography, 15(7): 985-996. DOI:10.1175/1520-0485(1985)015<0985:IROAPF>2.0.CO;2
Santoso A, McPhaden M J, Cai W J. 2017. The defining characteristics of ENSO extremes and the strong 2015/2016 El Niño. Reviews of Geophysics, 55(4): 1 079-1 129.
Schopf P S, Suarez M J. 1988. Vacillations in a coupled oceanatmosphere model. Journal of the Atmospheric Sciences, 45(3): 549-566. DOI:10.1175/1520-0469(1988)045<0549:VIACOM>2.0.CO;2
Spall M A, Pedlosky J. 2005. Reflection and transmission of equatorial Rossby waves. Journal of Physical Oceanography, 35(3): 363-373. DOI:10.1175/JPO-2691.1
Wang J, Yuan D L. 2015. Roles of western and eastern boundary reflections in the interannual sea level variations during negative Indian Ocean dipole events. Journal of Physical Oceanography, 45(7): 1 804-1 821. DOI:10.1175/JPO-D-14-0124.1
Wang Z, Yuan D L. 2012. Nonlinear dynamics of two western boundary currents colliding at a gap. Journal of Physical Oceanography, 42(11): 2 030-2 040. DOI:10.1175/JPO-D-12-05.1
Wang Z, Yuan D L. 2014. Multiple equilibria and hysteresis of two unequal-transport western boundary currents colliding at a gap. Journal of Physical Oceanography, 44(7): 1 873-1 885. DOI:10.1175/JPO-D-13-0234.1
Weisberg R H, Wang C Z. 1997. A western Pacific oscillator paradigm for the El Niño-Southern oscillation. Geophysical Research Letters, 24(7): 779-782. DOI:10.1029/97GL00689
Yu L S, Weller R A, Liu W T. 2013. Case analysis of a role of ENSO in regulating the generation of westerly wind bursts in the western equatorial Pacific. Journal of Geophysical Research, 108(C4): 3 128.
Yuan D L, Han W Q. 2006. Roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean. Journal of Physical Oceanography, 36(5): 930-944. DOI:10.1175/JPO2905.1
Yuan D L, Liu H L. 2009. Long-wave dynamics of sea level variations during Indian Ocean dipole events. Journal of Physical Oceanography, 39(5): 1 115-1 132.
Yuan D L, Rienecker M M, Schopf P S. 2004. Long wave dynamics of the interannual variability in a numerical hindcast of the equatorial Pacific Ocean circulation during the 1990s. Journal of Geophysical Research, 109(C5): C05019.
Zhao X, Yuan D L, Yang G, Wang J, Liu H L, Zhang R F, Han W Q. 2019. Interannual variability and dynamics of intraseasonal wind rectification in the equatorial Pacific Ocean. Climate Dynamics, 52(7): 4 351-4 369.