The roles of tropical and subtropical wind stress anomalies in the El Niño Modoki onset

The El Niño Modoki is separated into two types, El Niño Modoki I and El Niño Modoki II, which are characterized by substantially different mechanism. According to comparisons between the evolution of air-sea coupled processes and heat budget analyses, this study discovered the distinct dynamical mechanism of the onset of El Niño Modoki II compared to that of El Niño Modoki I. For El Niño Modoki I, westerly wind anomalies in the western Pacific appear earlier than warm sea surface temperature (SST) anomalies in the central tropical Pacific; in addition, zonal advection is responsible for changes in mixed-layer temperature in the Niño-4 region during the developing phase in spring and summer. However, warm SST anomalies accompany easterly anomalies in the central tropical Pacific, and there are no westerly anomalies in the western tropical Pacific during the developing phase in spring and summer for El Niño Modoki II. For El Niño Modoki II, Ekman feedback is the major contributor to changes in mixed-layer temperature. According to reanalysis data, the 2.5-layer model results and coupled model simulations, anomalous westerly wind stress in the subtropical North Pacific and easterly wind stress in the central–eastern tropical Pacific can lead to equatorward and westward oceanic current anomalies, thus enhancing meridional and zonal convergences at upper layer. Because anomalous convergence at upper layer can inhibit upwelling motion in the central tropical Pacific, SST anomalies in the central tropical Pacific tend to be warm, even without anomalous westerlies in the western Pacific.


Introduction
Recent studies suggest that there exist two types of El Niño events: canonical El Niño and El Niño Modoki (Trenberth and Stepaniak 2001;Larkin and Harrison 2005;Ashok et al. 2007;Yu and Kao 2007;Kao and Yu 2009;Kug et al. 2009). During a canonical El Niño, also referred to as the eastern Pacific El Niño, anomalously warm sea surface temperatures (SSTs) are centered at the cold-tongue region in the eastern tropical Pacific. During an El Niño Modoki, also referred to as the dateline El Niño (Larkin and Harrison 2005), or central Pacific El Niño (Yu and Kao 2007), warm SST anomalies are located in the central tropical Pacific. It is noted that El Niño Modoki proposed by Ashok et al. (2007) exhibits dipole-like SST anomalies in the central-eastern tropical Pacific, which is different from warm pool El Niño (Kug et al. 2009). Lian and Chen (2012) suggested that such dipole-like mode of SST anomalies results from the orthogonality constraint of the EOF analysis, and can not represent the real pattern. Regarding these two types of El 1 3 Niño events, not only are the distributions and evolutions of SST anomalies different, but the impacts on global climate are also different (Kumar et al. 2006;Ashok et al. 2007;Taschetto and England 2009;Weng et al. 2011;Liu et al. 2014;Xu et al. 2018). Studies on the El Niño Modoki have attracted an increasing amount of attention (Capotondi et al. 2015;Yu et al. 2017).
Different mechanisms have been proposed to explain the occurrence of El Niño Modoki (Ashok et al. 2007;Kug et al. 2009;Yu et al. 2010;Chen et al. 2015;Lai et al. 2015). Ashok et al. (2007) suggested that warming in the central Pacific is due to deepening of the thermocline, which is related to eastward-propagating down-welling Kelvin waves induced by westerly anomalies in the western equatorial Pacific and westward-propagating down-welling Rossby waves induced by easterly anomalies in the eastern equatorial Pacific (Bjerknes 1969). Kug et al. (2009) noted that wind-driven thermocline variations may not be adequate for producing SST anomalies seen in an El Niño Modoki because the mean thermocline is relatively deep in the central Pacific. Instead, they argued that warming in the central tropical Pacific develops mainly from oceanic advection, and warming in the eastern tropical Pacific is suppressed due to enhanced upwelling and evaporation caused by equatorial easterly anomalies. Yu et al. (2010) emphasized the effects of oceanic advection on the anomalously warm SST in the central tropical Pacific. They stated that SST anomalies first emerge in the northeastern subtropical Pacific, then spread to the central equatorial Pacific through wind-evaporation-SST (WES) feedback (Xie and Philander 1994). Chen et al. (2015) concluded that warm water volume (WWV) and westerly wind bursts (WWBs) in the equatorial Pacific work together to generate the diverse El Niño events. Lai et al. (2015) reported that the diversity of El Niño events is a joint result of zonal wind anomalies in the central-western Pacific (140°E-160°W) during the developing year and subsurface temperature anomalies in the western Pacific during the previous winter. Their study showed that cumulative zonal wind anomalies are more decisive than WWBs.  discovered that El Niño Modoki events can be classified into two types, El Niño Modoki I and El Niño Modoki II, and postulated an index to identify them (Wang et al. 2018). SST anomalies during El Niño Modoki I first appear in the central equatorial Pacific and develop locally, whereas SST anomalies during El Niño Modoki II first appear in the northeastern subtropical Pacific, which later spread towards and peak in the central equatorial Pacific. These two types of El Niño Modoki events show distinct features in SST warming position and regional climatic effects Wang 2013, 2014;Tan et al. 2016). From the observation and model results, El Niño Modoki I results in anticyclone circulation anomalies in the Philippine Sea in developing autumn via Rossby wave responses, while El Niño Modoki II induces cyclone circulation anomalies in the east of the Philippines Wang et al. 2018). The anomalous anticyclone (cyclone) circulations associated with El Niño Modoki I (El Niño Modoki II) can thus lead to the increase (decrease) of rainfall in southern China.  demonstrated that the weakening of the Walker circulation in the Indo-Pacific during the canonical El Niño and El Niño Modoki I is in favor of a positive Indian Ocean Dipole (IOD), the mode that is argued independent of ENSO (Saji et al. 1999). On the other hand, the Walker circulation is enhanced during El Niño Modoki II, and thus tends to induce a negative IOD. Moreover, a weakly positive IOD in 2015 is suggested to be associated with El Niño Modoki II . Tan et al. (2016) showed that during the developing fall, SSTs in the South China Sea (SCS) exhibit negative (positive) anomalies during El Niño Modoki II (canonical El Niño and El Niño Modoki I). They suggested that it is the different position of the anomalous anticyclone in the western Pacific that results in distinct latent heat flux anomalies in the SCS, which creates opposite SST variations.
As shown in previous studies, climate impacts of El Niño Modoki II are distinctly different from those of canonical El Niño and El Niño Modoki I. However, the air-sea coupled preconditions and onset mechanisms of El Niño Modoki II remain unknown. Thus, the goal of this study is to examine the mechanism of El Niño Modoki II events. This paper is organized as follows. Section 2 introduces the datasets, methods and models used in this study. Air-sea coupled evolutions and a mixed-layer heat budget analysis during the developing phase of El Niño Modoki I and II events are compared in Sect. 3. In Sect. 4, the roles of ocean dynamics associated with wind stress anomalies in the subtropical and tropical Pacific during the El Niño Modoki II onset are addressed based on reanalysis datasets, the 2.5-layer model and coupled model results. Section 5 presents the summary and discussion.

Datasets
Oceanic current and temperature data from the Simple Ocean Data Assimilation (SODA) 2.2.4 (Carton and Giese 2008;Giese and Ray 2011) and surface heat flux data from the Twentieth Century Global Reanalysis Version 2 (20CR v2; provided by the NOAA/OAR/ESRL PSD at https ://www.esrl.noaa.gov/ psd/) (Compo et al. 2011) are utilized. The SODA 2.2.4 dataset has a horizontal resolution of 0.5° × 0.5° and 40 unevenly distributed vertical levels (resolution ranges from 5 m near the surface to 250 m at an approximate 5000 m depth). The 20CR v2 data are on a global T62 Gaussian grid. Wind stress fields are taken from SODA 2.2.4. Sea surface wind data are also obtained from the 20CR v2 dataset. SST data are from the Met 1 3 Office Hadley Centre's sea ice and sea surface temperature (HadISST1) dataset with a resolution of 1° × 1° (Rayner et al. 2003). These monthly atmospheric and oceanic variables cover the period of 1911-2010. Anomalies presented in this paper are the deviation from the monthly mean climatology of the last 30 years . The reanalysis data may have their own uncertainties, but we found that these conclusions are insensitive to the selection of datasets or time ranges. Consistent results can be achieved using other reanalysis datasets with different time ranges (e.g., NCEP Global Ocean Data Assimilation System reanalysis data). In addition, monthly potential temperature and potential density from Levitus94 (Levitus and Boyer 1994) are applied for the calculation of model parameters.

Heat budget analysis of mixed-layer temperature
The equation for the mixed-layer temperature heat budget (Huang et al. 2010) is where the temperature tendency (T t ) is decomposed into five terms, including zonal advection (Q u ) , meridional advection (Q v ) , vertical advection (Q w ) , net surface heat flux (Q q ) and the residuals (R). These terms are defined as: Here, T is the mixed-layer temperature, and T t is based on backward time difference (e.g., T t (May) = T(May)−T(April)). u, v, and w represent the zonal, meridional and vertical velocity, respectively. w e is the entrainment velocity below the mixed layer ( written as w here after), ρ is the water density, c p is the specific heat capacity of sea water under constant pressure, and h is the depth of the mixed layer.
By decomposing Eq. (1) into climatology and its anomaly, we obtain the following Eq. (7) for anomalous temperature where the prime symbols indicate monthly anomalies. Using the overbar to refer to climatology, the terms on the right side of the equation can be written as Monthly data are used to calculate the heat budget terms. Several sub-monthly processes that might affect the heat balance, such as tropical instability waves, cannot be sufficiently resolved. However, heat budget analyses based on monthly data have been widely utilized and are adequate for the diagnosis of ENSO (Zhang et al. 2007;Huang et al. 2010;Kug et al. 2010;Graham et al. 2017).

2.5-layer ocean model
A 2.5-layer ocean model (Huang 1987a, b) is used to explore the dynamic influence of wind stress anomalies on the tropical ocean. The 2.5-layer model simplifies the ocean into a rectangular basin with three layers. The first and second layers represent water in the mixed layer and the thermocline, respectively. The bottom layer is assumed to be infinitely deep and motionless. There is no water mixing between each layer.
The momentum and continuity equations for the first and second layers are where f is the Coriolis parameter (for a symmetric twohemisphere model, f = (y − L∕2) , where L is the meridional extent of the basin), h 1 and h 2 are the layer thickness of the first and second layer, respectively, (u 1 , v 1 ) and (u 2 , v 2 ) are the horizontal velocity of each layer, g ′ 1 and g ′ 2 are the reduced gravities, k 1 and k 2 are the vertical friction coefficients (both set to 0.0015 m/s), A m is the momentum mixing coefficient along the isopycnal surface (set to 1.5 × 10 4 m 2 /s), and τ is the zonal wind forcing. Subscripts t , x and y denote the time, zonal and meridional partial derivatives, respectively.
The domain in the model simulation is set to (40°S-40°N, 120°E-103°W), with a horizontal grid of 138 × 82 (grid size set to 111 km, about 1.0° × 1.0° resolution). The thicknesses of the first and second layers of the model, and the corresponding reduced gravity coefficients, are calculated from long-term mean monthly data from Levitus94. The primary domain for El Niño events is the tropical Pacific, and thus the layer thicknesses and reduced gravities are calculated as the mean value in the tropical Pacific (30°S-30°N, 100°E-80°W). The sum of the first and second thicknesses, which represents the depth of the thermocline, is defined as the depth of the 20 °C isotherm. The depth of the mixed layer is defined as the depth where the water temperature is 0.5 °C below the sea surface temperature (Sprintall and Tomczak 1992). Following the above definitions, the depths for the mixed layer and the thermocline from Levitus94 are 66 m and 151 m, respectively. For simplicity, the initial thickness of the first layer is set as H 1 = 65 m, and the initial thickness of the second layer is set as H 2 = 85 m. The reduced gravities of the two layers are calculated as where g is the gravity acceleration, ρ 1 is the density of the first layer, ρ 2 is the density of the second layer, and ρ 3 is the density of the bottom layer. The results are g � 1 = 0.015 m/s 2 and g � 2 = 0.026 m/s 2 .

Coupled model
The In the study, a fully coupled running with pre-industrial conditions is firstly conducted. The model is running 1200 years to reach the stable.

Comparison of evolutions for two types of El Niño Modoki
According to , El Niño Modoki I (El Niño Modoki II) is identified if there are positive (negative) rainfall anomalies in southern China during boreal autumn.
Seven El Niño Modoki I (1914, 1940, 1941, 1963/1964, 1987/1988, 1990/1991, and 2002/2003) and five El Niño Modoki II events (1968/1969, 1979/1980, 1991/1992, 2004/2005, and 2009/2010) are analyzed in this study. Different from the definition in , the 1992/1993 El Niño Modoki II is excluded because its warming peak is not in winter.  illustrated the evolutions of SST anomalies from developing spring (March-April-May, MAM) to the following winter (December-January-February, DJF) during El Niño Modoki I and El Niño Modoki II, respectively. For El Niño Modoki I, warm SST anomalies appear in the central equatorial Pacific, then develop and peak in boreal winter. However, for El Niño Modoki II, anomalously warm SSTs appear in the northern subtropical Pacific in boreal spring, then extend towards the tropics, and peak in the central equatorial Pacific in winter.
In addition to the differences in warm SST anomaly spatial patterns, this study further illustrates that the evolution of air-sea interactions in the tropical Pacific during El Niño Modoki II is distinct from that during El Niño Modoki I. Figure 1 displays the evolutions of composited sea surface wind and SST anomalies in the equatorial Pacific during the development of El Niño Modoki I and El Niño Modoki II, as well as their differences. For El Niño Modoki I, the westerly anomalies first appear in the western equatorial Pacific between 135°E-165°E during early spring and are accompanied by warm SST anomalies in the central equatorial Pacific (Fig. 1a). Such strong westerly anomalies persist into the developing boreal summer and autumn, and are accompanied by prominently warm SST anomalies (warmer than 0.5 °C) in the central and eastern equatorial Pacific. Previous studies have noted that WWBs can trigger downwelling Kelvin waves at the equator, and advect warm surface water eastward and favor the development of El Niño Modoki events (Hu et al. 2014;Lian et al. 2014;Chen et al. 2015;Fedorov et al. 2015;Lai et al. 2015). It is suggested that the onset of El Niño Modoki I is closely associated with air-sea interactions in the tropical Pacific. However, for El Niño Modoki II, there are warm SST anomalies in the central tropical Pacific in developing spring although the intensity is weak, but the westerly anomalies are not seen in the western Pacific until developing summer (Fig. 1b). The easterly anomalies appear between 180°E and 120°W during the developing earlier spring, with warm SST anomalies in the central tropical Pacific. It seems that the warm SST anomalies in the central tropical Pacific occur earlier than the anomalous westerly in the western tropical Pacific in El Niño Modoki II, which is different from that in El Niño Modoki I. Figure 1c illustrates the different evolutions of air-sea interaction between El Niño Modoki I and II. It is seen clearly that the warm SST anomaly intensities during developing spring of El Niño Modoki I are significantly larger than those of El Niño Modoki II, which is similar to 1 3 previous studies Wang et al. 2018). It is also indicated that El Niño Modoki I appears earlier than El Niño Modoki II. The other distinct differences between El Niño Modoki I and II are that the westerly anomalies in the western and central tropical Pacific during developing spring of El Niño Modoki I are larger than those of El Niño Modoki II. That is, the warm SST anomalies in the central tropical Pacific during the developing year of El Niño Modoki II can be developed even with weak westerly wind anomalies in the western tropical Pacific. Here, a question is raised regarding how warm SST anomalies in the central Pacific occur during the developing spring and summer of El Niño Modoki II.
To investigate how warm SST anomalies evolve during various El Niño Modoki events, the tendency of the mixed-layer temperature in the Niño-4 region (5°S-5°N, 160°E-150°W) and associated physical processes are analyzed. During boreal spring and summer (February-March-April-May-June, FMAMJ) of El Niño Modoki I, zonal advection (−uT x ) is a major contributor to warm anomalies in the Niño-4 region, according to the mixedlayer temperature budget analysis (Fig. 2a). This is consistent with previous studies suggesting that zonal advection is the principal contributor (Kug et al. 2009(Kug et al. , 2010Yu et al. 2010;Capotondi et al. 2015). The contribution of zonal advection is twice as much as those of meridional advection. The meridional and vertical advections in mixed-layer temperature change account for approximately 50% and 25% of the mixed-layer temperature change, respectively. The net surface heat flux is the damping term. Distinct from El Niño Modoki I, vertical advection (−wT z ) dominates the SST anomaly tendency during the developing spring and summer of El Niño Modoki II (Fig. 2b). The contribution of vertical advection is as high as about 90%, and is much larger than that of zonal advection (25%) and meridional advection (30%). The net surface heat flux during El Niño Modoki II acts as a damping term as well.
This study focuses on vertical advection during the developing FMAMJ of El Niño Modoki II since the importance of horizontal advection on the development of the El Niño Modoki (so-called El Niño Modoki I in this study) has been revealed by many researchers (Kug et al. 2009(Kug et al. , 2010Yu et al. 2010;Capotondi et al. 2015). Vertical advection during El Niño Modoki II is further decomposed into four terms: subsurface temperature anomalies (−wT � z ) , upwelling anomalies (−w � T z ) and two nonlinear terms of (−w � T � z ) and (w � T � z ) . Their contributions to temperature tendencies are shown in Fig. 3. The two former terms indicate thermocline feedback and Ekman feedback, respectively. It is shown that the amplitude of the Ekman feedback is the largest, which dominates the changes of vertical advection term, and the nonlinear term takes on a secondary role. The thermocline feedback is negative, indicating that it is unfavorable for El Niño Modoki II development. It is illustrated that oceanic dynamic processes in the vertical direction are much more important than thermal stratification changes for SSTA warming during the onset of El Niño Modoki II.

The possible mechanism for the onset of El Niño Modoki II
The heat budget analysis indicates that the contribution of vertical advection to SST tendency in the central tropical Pacific is larger than that of horizontal advection during El Niño Modoki II. The physical processes leading to the large contribution of vertical advection are examined in this section, and a possible mechanism of El Niño Modoki II onset is proposed. Under climatological mean condition, trade winds and the geostrophic effect lead to the poleward flow of the upper layer water in the central equatorial Pacific, inducing the upwelling motions in the central tropical Pacific (Fig. 4a). In the climatological mean state, divergences of upper-layer currents are seen in central and eastern tropical Pacific from 180°E to 140°W (Fig. 4b). During the developing FMAMJ of El Niño Modoki II, downwelling anomalies weakening climatological upwelling motions at the equator appear (Fig. 4c), which is in agreement with the heat budget analysis shown above (Figs. 2b, 3). Moreover, the downwelling anomalies can extend as deep as 150 m near equator, which is the mean depth of thermocline. Therefore, the wind-induced oceanic currents above the thermocline and divergences are investigated in the followings. From Fig. 4c, it is seen that downwelling motion anomalies are associated with southward current anomalies in the upper layer of the Northern Hemisphere, indicating decreases in the subtropical cell (STC). Figure 4d displays equatorial-averaged (5°S and 5°N) divergence anomalies during the developing FMAMJ of El Niño Modoki II. The anomalous convergence of the current velocities appears in most parts of the central tropical Pacific, with a maximum around 160°W (Fig. 4d). Anomalous convergence at the upper layer, induced by horizontal currents, can lead to weakened upwelling at the equator and less cold water upwelled into the mixed layer. Hence, SST anomalies tend to be warm. Since oceanic surface circulations are significantly influenced by wind stress changes (Luo et al. 2015), wind stress anomalies in the developing spring of El Niño Modoki I and El Niño Modoki II are investigated (Fig. 5). The significant westerly wind stress anomalies are seen in the western and central tropical Pacific during developing FMAMJ of El Niño Modoki I accompanying with significant easterly wind stress anomalies in the eastern tropical Pacific (Fig. 5a). Different from El Niño Modoki I, westerly wind stress anomalies are weak at the equator during FMAMJ of El Niño Modoki II. Instead, anomalously strong westerly wind stress is seen in the northern subtropical Pacific (10°N-20°N, blue box in Fig. 5b), which is associated with cyclonic anomalies in the northern subtropical Pacific (Yu et al. 2010;Yeh et al. 2015). In addition, strong easterly wind stress is found in the central-eastern tropical Pacific between 160°W and 130°W for El Niño Modoki II (green box in Fig. 5b), which is consistent with the evolution of surface winds shown in Fig. 1. The location of the significant easterly wind stress anomalies for El Niño Modoki II is more westward than those of El Niño Modoki I.
To examine whether the wind stress anomalies in the northern subtropical Pacific and eastern tropical Pacific in Fig. 5b can result in the changes in upper-layer circulation and decreases in divergent flow, several numerical experiments based on a 2.5-layer model are conducted. A control run experiment and three sensitivity experiments are designed ( Table 1). The model domain covers the Pacific (40°S-40°N, 120°E-103°W). In the control run, the model is forced by the climatological zonal wind stress field for 50 years to reach a quasi-steady state. All three sensitivity experiments were restarted from the final state of the control run; in addition, the model was forced by different wind stress for 3 months.
The zonal wind stress used in the control run is shown in Fig. 6. It exhibits a saddle-shaped double-peak distribution in the meridional direction. Easterlies between 35°S and 30°N reach a maximum of − 0.06 N/m 2 and − 0.09 N/m 2 at   Zonal mean climatology plus wind stress anomalies in the blue box and green box in Fig. 5 1 3 20°S and 15°N, respectively, and a minimum of − 0.025 N/ m 2 at the equator. Figure 7 shows the quasi-steady state at the end of the control run experiment. Due to the easterly wind, water piles up in the west part of the basin, leading to the west-east upward slope of thermocline. The thermocline depth is more than 160 m at the western boundary and less than 60 m at the eastern boundary (Fig. 7a), which is similar to the long-term mean state from reanalysis data. The subtropical gyres in both hemispheres and the equatorial counter-currents are reproduced in the control run (Fig. 7b). Thus, the model is capable of simulating the relationship between wind stress and upper layer circulation in the Pacific.
Three sensitivity experiments, referred to as the subtropical force (SF) experiment, tropical force (TF) experiment and subtropical and tropical force (STF) experiment (Table 1), are conducted. In these sensitivity experiments, the model is forced by the climatological zonal wind stress, plus wind stress anomalies in the blue box, green box and both boxes shown in Fig. 5b, respectively. The differences in ocean velocity above the thermocline and its divergence between the sensitivity experiments and the control run are shown in Fig. 8. In the SF experiment, the strong equatorward current anomalies from the subtropical regions to the tropics appear around 150°W, which are associated with anomalous divergence centered in the subtropical Pacific and anomalous convergence centered near the tropical Pacific (around 150°W) (Fig. 8a). That is, the anomalous westerly wind stress in the northern subtropical Pacific can induce current anomalies equatorward in the upper layer and weaken the STC, which is similar to those inferred from reanalysis (Fig. 4c). Such equatorward flow anomalies result in anomalous convergence in the upper layer of the tropical Pacific near 150°W (Fig. 8d), and reduce climatological divergences over there. Hence, the results of the SF experiment suggest that wind stress in the northern subtropical Pacific (blue box in Fig. 5b) can induce anomalous convergence in the tropical Pacific through weakening the STC. For the TF experiment, anomalous westward currents associated with wind stress anomalies in the central-eastern tropical Pacific (Fig. 5b) are seen in the tropics from 160°W to 120°W. The maximum of westward current anomalies located near 150°W induces the anomalous convergence and divergence in the west and east of 140°W (Fig. 8b), which is consistent with the divergence pattern in the reanalysis as shown in Fig. 4d. There are some differences between the SF and TF experiment results, although both these two experiments simulate convergences in the tropics, which is quite similar to those inferred from the reanalysis. The intensity of the anomalous convergence in the tropical Pacific from the SF experiment (Fig. 8d) is about four times as high as that from the TF experiment (Fig. 8e). In addition, the anomalous convergences in the tropics in the SF experiment result from meridional current anomalies, while zonal current anomalies are responsible for tropical anomalous convergence in the TF experiment.
The results of the STF experiment are linear combinations of results obtained from the SF and TF experiments. The southward current anomalies are from the subtropics to the tropics, and the westward currents in the tropics are strong around 150°W (Fig. 8c). Two divergence centers and one convergence center extend from the subtropics to the tropics, and the intensity of the convergence in the tropics around 150°W is the strongest (Fig. 8c), which is approximately the sum of the convergences from the SF and TF experiments (Fig. 8f). The 2.5-layer model results confirm that anomalous wind stress in the northern subtropical Pacific and central-eastern tropical Pacific (Fig. 5b) can induce convergence anomalies in the upper layer in the central tropical Pacific.
As the 2.5-layer model cannot directly simulate SST and ocean vertical motions, the coupled model, CESM1, is used to investigate the effect of wind stress anomalies in the subtropical Pacific on the onset of El Niño Modoki II. In the observations, the significant westerly anomalies are two times as high as the significant easterly anomalies in the tropics (Fig. 5b). The anomalous convergence induced by the subtropical westerlies in the 2.5-layer model is twice larger than that induced by the equatorial easterlies (Fig. 8d, e). It is indicated that the subtropical westerlies play more important role in development of El Niño Modoki II. Therefore, only the influences of the westerly anomalies in the northern subtropical Pacific are checked in the coupled model experiments.
After 1200 model years of fully coupled running under pre-industrial conditions, the wind stresses from the Coordinated Ocean-Ice Reference Experiments version 2 (COREv2) normal year forcing (Large and Yeager 2004) are used to force ocean in model. Here, the wind stresses are varied on monthly climatological. The simulation is 100 years, and is referred to as the control run. The sensitive experiment is identical to the control run except that the wind stress is imposed by the subtropical (Fig. 5b, blue box) westerly anomaly composites of El Niño Modoki II from January to June. The sensitive experiment is run for 100 years. The last 30 years simulations of control run and sensitive experiment are used in our study. The differences between control run and sensitive experiment could represent the influences of wind stress in northern subtropical Pacific. Figure 9 displays the annual mean SST in the last 30 years of the control run. The distribution of SST in Pacific is well reproduced although there is the cold tongue bias, which is common limitation in current climate model performances (Li and Xie 2012;Zheng et al. 2012). The differences between sensitive experiment and control run during April and June are shown in Fig. 10. Under the forcing of wind stress in northern subtropical Pacific, the significant warm SST anomalies are seen in the central and eastern tropical Pacific and extending to subtropical Pacific (Fig. 10a), which is similar to the SST anomalies pattern of El Nino Modoki II suggested by  and Wang et al. (2018). These warm SST anomalies are associated with upper-layer oceanic circulation (Fig. 10a, b). There are equatorward currents at upper-layer along 180°E-160°W, and the downwelling motions are obvious near the equator (Fig. 10b), which is similar to the observation in Fig. 4c. It is noted that since the anomalous westerlies are imposed in the North Pacific, the responses of meridional circulations in the southern hemisphere are weak. The westerly anomalies in the subtropical Pacific can induce the upper-layer westward current anomalies in (160°E-160°W, 0-10°N), and the convergence of currents are around (180°E-160°W, 0-10°N). Thus, it is confirmed by the coupled model that the westerly anomalies in the northern subtropical Pacific may enhance the equatorward currents in the central tropical Pacific, and then result in the downwelling anomalies and current convergences at upper-layer, which are favorable for warm SST anomaly appearances.

Summary and discussion
The evolution of coupled SST and surface wind anomalies in the equatorial Pacific during El Niño Modoki I and El Niño Modoki II events defined by  are analyzed and compared in this study. In the El Niño Modoki I, westerly wind anomalies appear earlier than warm SST anomalies in the equatorial Pacific, which is in agreement with previous studies; thus, anomalous westerlies in the western Pacific play a key role in triggering El Niño Modoki by extending the warm pool eastward (Chen et al. 2015;Lai et al. 2015). However, SST anomalies might be positive, even in the absence of strong westerly anomalies in the western tropical Pacific during El Niño Modoki II. In addition, easterly anomalies are accompanied by warm SST anomalies in the central tropical Pacific. These differences in SSTs and surface wind evolutions suggest that the dynamical progression of El Niño Modoki II is different from those of El Niño Modoki I. The heat budget analysis of the mixed layer temperature in the Niño-4 region during the El Niño Modoki II developing phase shows that Ekman feedback is a major contributor to the warming SST tendency in the central tropical Pacific, which is quite different from that of El Niño Modoki I. Yang and Zhang (2008) also emphasized the important role of the anomalous upwelling in modulating the ENSO variability. Our analysis further suggests that the Ekman feedback is substantially related to wind stress anomalies in the northern subtropical Pacific and the central-eastern tropical Pacific. From reanalysis datasets and model experiments, subtropical westerly and equatorial easterly wind stress anomalies can influence mixed-layer flow and result in anomalous convergence in the central equatorial Pacific. A schematic diagram is displayed to summarize the proposed El Niño Modoki II mechanism (Fig. 11). Subtropical wind stress anomalies lead to Ekman transport towards the equator, and equatorial wind This study reveals the influence of surface wind in the northern subtropical Pacific and central-eastern equatorial Pacific on mixed-layer temperature changes in the central tropical Pacific during the onset of El Niño Modoki II. This is different from previous studies regarding subtropical-tropical interplays. Gu and Philander (1997) and Kleeman and Klinger (1999) emphasized temperature changes or water transports associated with the Pacific STC that modulate ENSO on interdecadal and decadal scales. However, in this study, we focus on the dynamic effect of subtropical and equatorial wind stress anomalies on the modulation of diverse El Niño events on an interannual scale. In addition, different from the WES feedback mechanism that involves the propagation of SST anomalies and wind anomalies (Vimont et al. 2009;Zhang et al. 2009;Yu et al. 2010;Lu et al. 2017), the impacts of ocean dynamics are emphasized in this study. A coupled model may be needed to further investigate the detail of the SST evolution during the onset of El Niño Modoki II, as well as the subsequent progression afterwards.
This study illustrates the roles of wind stress anomalies in the northern subtropical Pacific and the central-eastern tropical Pacific on the El Niño Modoki II onset. However, the causes of these wind stress anomalies remain unknown. It is speculated that the wind stress anomalies in the northern subtropical Pacific and the central-eastern tropical Pacific may be related to variations in the North Pacific Oscillation (NPO, Yeh et al. 2015) and the southern Pacific/tropical Atlantic (Ham et al. 2013a, b), respectively. These questions will be analyzed in future studies. Fig. 11 A schematic diagram illustrating the appearance of warm SST anomalies in the central tropical Pacific during El Niño Modoki II induced by anomalous westerly anomalies in the northern subtropical Pacific and easterly anomalies in the eastern tropical Pacific. The green and blue thin arrows (1) depict anomalous westerly anomalies in the northern subtropical Pacific and easterly anomalies in the eastern tropical Pacific, respectively. The green and blue bold arrows (2) depict anomalous equatorward currents and westward currents in the upper layer resulting from (1). Black bold arrows (3) show the anomalous downwelling motion induced by anomalous convergence at the surface. The red area at the sea surface (4) indicates warming SST