Impacts of seasonal transitions of ENSO on 2 atmospheric river activity over East Asia 3

29 Atmospheric rivers (ARs), narrow water vapor transport bands over the mid-latitudes, 30 often cause great socio-economic impacts over East Asia. While it has been shown that 31 summertime AR activity over East Asia is strongly induced by preceding-winter El Niño 32 development, it remains unclear the extent to which seasonal transitions of El Niño 33 Southern Oscillation (ENSO) from winter to summer affect the AR activity. Here we 34 examine the relationship between the seasonal transitions of ENSO and the summertime 35 AR activity over East Asia using an atmospheric reanalysis and high-resolution 36 atmospheric general circulation model (AGCM) ensemble simulations. A rapid transition 37 from preceding-winter El Niño to summertime La Niña results in more AR occurrence over 38 northern East Asia via northward expansion of an anomalous low-level anticyclone over the 39 western North Pacific compared to sustained or decayed El Niño cases. The northward 40 expansion of the anticyclone is consistent with a steady response of the atmosphere to the 41 anomalous condensation heating over the Maritime Continent and equatorial Pacific. 42 Meridional positions of the extratropical AR occurrence and circulation anomalies are 43 different between the reanalysis and AGCM simulations, which is possibly contributed by a 44 limited sample size and/or AGCM biases and suggests that seasonal prediction of 45 AR-related natural disaster risk over East Asia on a regional scale remains a challenge.


Introduction
Atmospheric rivers (ARs), narrow moisture transport bands developing over the mid-latitudes, have been extensively studied (e.g.Guan and Waliser 2015;Gimeno et al. 2016;Shields et al. 2018; American Meteorological Society 2019) because of their great impacts on society through heavy rainfall and associated natural disasters (e.g.Ralph et al. 2011;Gimeno et al. 2014;Waliser and Guan 2017).The regional occurrence of ARs is highly dependent on seasons, especially in the northeastern and southeastern Pacific, East and South Asia, and the Middle East (Guan and Waliser 2015).Over East Asia, frequent AR genesis especially during summer (Guan and Waliser 2019) often causes extreme rainfall events (Kamae et al. 2017b, c).For example, ARs contribute to heavy rainfall events of August 2014 (Hirota et al. 2016;Tsuji and Takayabu 2019) and July 2018 (Shimpo et al. 2019;Takemura et al. 2019;Tsuguti et al. 2019), which killed 75 and 237 people, respectively.
Recently, interannual variations and future changes in AR occurrence frequency over East Asia have been studied using atmospheric reanalyses and model simulations.Kamae et al. (2017b) found that summertime AR frequency over East Asia tends to increase following wintertime El Niño development with a half-year lag.A key factor connecting El Niño-Southern Oscillation (ENSO) and the East Asian AR frequency is the western North Pacific subtropical high (WNPSH).The local effect of concurrent cold anomaly in underlying sea surface temperature (SST; Wang et al. 2000) and the delayed effect of the North Indian Ocean and South China Sea warming after the wintertime El Niño intensify the WNPSH by increasing surface divergence over the western subtropical Pacific (Xie et al. 2009(Xie et al. , 2016)).
Through this "Indian-western Pacific Ocean Capacitor (IPOC)" effect, AR frequency is increased on the northwestern flank of the anomalous WNPSH in post El Niño summers (Kamae et al. 2017b).Kamae et al. (2019) further pointed out that the Indo-Pacific SST is also the key to future projections of East Asia-western North Pacific AR frequency by analyzing atmospheric general circulation model (AGCM) simulations forced by different SST warming patterns.An increase in atmospheric water vapor under global warming results in an increase in AR frequency over the Northern Hemisphere mid-latitudes (see also Espinoza et al. 2018).In addition, warmer Indian Ocean and South China Sea are responsible for a larger increase in East-Asian ARs via favoring an anomalously strong WNPSH and southwesterly anomalies on its northwestern flank (Kamae et al. 2019).
In contrast to the IPOC effect, summertime AR response over the western North Pacific to the concurrent anomaly in the equatorial Pacific SST has not been well understood.Mundhenk et al. (2016) examined summertime ENSO impact by compositing El Niño and La Niña years and concluded that ENSO has no systematic effects on East Asian ARs (Fig. 8a, c in Mundhenk et al. 2016).This is partly due to limited sample size: reanalysis data used in this study only cover 36 years (i.e., 1979-2014).In contrast, Kamae et al. (2019) pointed out that summertime cool SST anomaly over the equatorial Pacific dominates the second leading mode of the interannual variations in summertime East-Asian AR frequency in a set of 10-member ensemble AGCM simulations.They suggested that the summertime La Niña-like SST anomaly may modulate the East Asian AR frequency via changing North Pacific atmospheric circulation (Wang et al. 2013;Paek et al. 2019).However, the physical connection between the summertime La Niña and the East Asian AR activity has not been systematically explored.
In the present study, the importance of summertime Pacific SST forcing to the East Asian AR frequency is studied using an atmospheric reanalysis and numerical simulations.
Atmospheric and oceanic observations of the recent decades indicate that ENSO exhibits asymmetries in strength, duration and transition (Stein et al. 2010 and references therein).
Strong El Niño developed during winter (sometimes called "super El Niño"), for example in 1997/1998 or 2015/2016, rapidly decays after spring and tends to show a rapid transition to La Niña.In contrast, La Niña tends to have moderate strength and longer duration than El Niño.In this study, the relationship between the seasonal development/transition of ENSO and the summertime AR frequency over East Asia and the physical mechanisms responsible for the relationship are examined.We use both an atmospheric reanalysis and an ensemble of high-resolution AGCM simulations to increase the sample size and to assess the robustness of the results.This paper is organized as follows.Section 2 describes data and methods including an atmospheric reanalysis, model simulations, and AR detection method.Section 3 compares summertime AR frequency over the North Pacific among different ENSO transitions from winter to summer.In section 4, we examine the possible importance of the tropical atmospheric heating/cooling anomalies related to ENSO transitions to the mid-latitude AR frequency.Section 5 presents a summary and discussion.

Observations, reanalysis, and ENSO classification
We use the Japanese 55-year Reanalysis (JRA-55; Kobayashi et al. 2015) at 1.25° × 1.25° spatial resolution to examine the historical variations in AR frequency and related atmospheric fields.The SST data used in this study are COBE-SST2 (Hirahara et al. 2014) which were used in the ensemble simulations described in section 2.2 as a boundary condition.The spatial resolution of this SST data set is 1.0° × 1.0°.To compare with the results of the ensemble simulations detailed in section 2.2, we examine the AR frequency during 1958-2010.
We divided 53 years into nine groups based on seasonal ENSO transition: combinations of preceding boreal-winter El Niño, neutral, and La Niña, and concurrent boreal-summer El Niño, neutral, and La Niña.Table 1 shows a summary of the clustering.
Based on SST anomaly over Niño3 domain (150°-90°W; 5°S-5°N), years with December-to-February (DJF)-mean anomaly larger than 1 K compared to the climatology  are defined as El Niño (> 1 K) or La Niña (< -1 K) winters.Totally there are 7 La Niña, 8 El Niño, and 38 neutral winters during 1958-2010 (Table 1).Similarly, June-to-August (JJA) SST anomaly is used for summertime clustering.During summer, the SST threshold used for classification is 0.5 K (> 0.5 K and < -0.5 K for El Niño and La Niña summers, respectively) because in this season Niño3 SST exhibits smaller variability than in boreal winter (Stein et al. 2010).Numbers of obtained La Niña, neutral, and El Niño summers are 14, 26, and 13, respectively (Table 1).Figure 1 shows the composite mean SST anomaly compared to the climatology.In addition to the summertime ENSO signals, lagged effects of wintertime ENSO are also found: for example, SSTs over the North Indian Ocean and South China Sea tend to be higher (lower) after wintertime El Niño (La Niña), consistent with the IPOC effect (Xie et al. 2009(Xie et al. , 2016)).The sustained ENSO years and ENSO transition years have limited samples (1, 2, 1 and 2 for sustained La Niña, sustained El Niño, La Niña-to-El Niño transition, and El Niño-to-La Niña transition years, respectively), suggesting difficulties in examining the effect of ENSO on AR activity by solely using the atmospheric reanalysis.Note that other factors including Atlantic SST may also affect results of the composite analyses especially for the sustained ENSO years or ENSO transition years because of the limited sample size (see sections 3.2 and 5).

Ensemble simulation in an AGCM
To increase the signal-to-noise ratio in ENSO-related variability in mid-latitude AR activity, we use an ensemble model simulation output, similar to Kamae et al. (2017bKamae et al. ( , 2019)).
We use outputs of 10-member simulation (for 1958-2010) with an AGCM.The model in use is the Meteorological Research Institute Atmospheric General Circulation Model (MRI-AGCM) version 3.2 (Mizuta et al. 2012) at a horizontal resolution of TL319 (equivalent to 60-km mesh) with 64 vertical layers.Horizontal resolution of moisture transport data used for AR detection (section 2.3) is identical to JRA-55 (1.25° × 1.25°).The AGCM was driven by historical  variations in SST and sea ice (Hirahara et al. 2014) and radiative forcing agents (anthropogenic greenhouse gases, aerosols, ozone, and natural aerosols).More details of the experiments can be found in Mizuta et al. (2017) and Kamae et al. (2017b).This ensemble simulation output, called the Database for Probabilistic Description of Future Climate Change (d4PDF), is used to evaluate AR response to global SST variability and assess the effect of atmospheric internal variability (Kamae et al. 2017a, b).In this study, the 10-member ensemble mean and inter-member spread are considered as forced atmospheric response to SST and radiative forcing and internal atmospheric variability, respectively (Kamae et al. 2017a, b, Ueda et al. 2018; see section 3.2).Although d4PDF ensemble simulations facilitate the comparison of the two components (i.e., forced atmospheric response and internal variability), they only cover limited ENSO cases during 1958-2010, resulting in a difficulty in obtaining robust atmospheric responses to the sustained or transitioning ENSO forcing (see sections 3.2 and 5).

Detection of atmospheric rivers
We detected ARs from 6-hourly atmospheric moisture transports in JRA-55 and d4PDF (sections 2.1, 2.2).Detection algorithm used in this study is based on Mundhenk et al. (2016) and slightly modified by Kamae et al. (2017b).In this algorithm, we use 6-hourly vertically integrated water vapor transport (IVT) between 300 and 1000 hPa levels calculated using specific humidity and horizontal wind at 1.25° × 1.25° spatial resolution.
IVT is determined as: where g is the acceleration due to gravity, q is specific humidity, and u and v are zonal and meridional components of horizontal winds, respectively.Anomaly in IVT from its daily climatology is used to detect ARs so that the effect of climatological seasonal cycle of water vapor transport is excluded.Anomalous IVT larger than 140 kg m -1 s -1 is detected every 6 hours.Next, the identified water vapor transport bands with small area (< 7.8 × 10 5 km 2 ) or short length (< 1,500 km) or small length/width ratio (< 1.325) are removed.Note that long-term trends of obtained summertime AR frequency over the northwestern Pacific are small compared to the interannual variability (Kamae et al. 2017b).More details of the detection methods can be found in Mundhenk et al. (2016) and Kamae et al. (2017b).

Linear baroclinic model
We use a linear baroclinic model (LBM; Watanabe and Kimoto 2000) to diagnose linear atmospheric response to tropical diabatic heating related to seasonal transition of ENSO.This model is based on primitive equations linearized around the climatology of the observed atmospheric state based on NCEP/NCAR reanalysis (Kalnay et al. 1996).The model resolution is T42 (~2.8°) in the horizontal and 20 levels in the vertical.Details of experimental setups including background state and prescribed forcing are explained in section 4.2.

AR frequency and atmospheric circulation over the North Pacific
Seasonal-mean AR frequency in the nine groups are compared with the climatology for 1958-2010.Figure 2 shows composite anomalies of AR frequency and atmospheric circulation in association with wintertime and summertime ENSO phases in the simulations.-c, 3a-c).In JRA-55, the circulation and AR anomalies are less systematic than d4PDF because of limited samples (Table 1).
In addition to the lagged influence, concurrent ENSO also modulates the North Pacific  , 3c) with small sample sizes (in 1971 and 1976, respectively).

East-Asian AR activity
The analyses in the previous section indicate that the seasonal transition of ENSO from winter to summer is a key for the prediction of summertime AR frequency over the North Pacific.Next, we focus on East Asia because this region often suffers from AR-caused heavy rainfall (Kamae et al. 2017c).Figure 4 shows East-Asian AR and circulation anomalies during post-El Niño summers in JRA-55 and d4PDF.To emphasize the effect of summertime ENSO, we compare the El Niño-to-La Niña transition years (Fig. 4a, c) and sustained El Niño years (Fig. 4b, d).During these years (i.e., both transition and sustained years), the subtropical western North Pacific experiences positive SLP anomalies and correspondingly increased AR frequency, as discussed in the previous section.Figure 5 shows AR frequency anomalies averaged over East Asia (110°-155°E; 25°-55°N, gray rectangle in Fig. 4).To assess year-to-year variability, we show AR frequency anomalies for individual transition and sustained years (see also Table 1).The East-Asian AR frequency (gray bars in Fig. 5) generally exhibit positive anomalies, consistent with Kamae et al. (2017b).No systematic differences are found between the transition years and sustained years (Fig. 5a, b).
In contrast, the regional patterns of the anomalies in SLP and AR frequency over East Asia are slightly different between the two (i.e., transition years and sustained years).The WNPSH is stronger and is more confined to the southern part during the sustained years (Fig. 4b, d) than the transition years (Fig. 4a, c).The WNPSH anomaly in d4PDF simulation is only found over 10°-32°N during the sustained years (Fig. 4d), contrasting to that with a larger meridional extension (10°-40°N) in the transition years (Fig. 4c).Similar difference in meridional extent is also found in JRA-55 (10°-30°N in Fig. 4b and 10°-48°N in Fig. 4a).
Note that the positions of the WNPSH anomalies during these years exhibit non-negligible difference between the JRA-55 and d4PDF, which will be discussed in section 5. Generally, NEA AR frequency in JRA-55 exhibits larger year-to-year variations than that in the 10-member ensemble mean of d4PDF simulations (Fig. 5).This can be partly attributed to atmospheric internal variability (Kamae et al. 2017a, b).Error bars in Fig. 5b are generally large, indicating the great importance of atmospheric internal variability under the fixed radiative forcing and SST boundary conditions in the d4PDF simulations.

ENSO transition and tropical rainfall
The results presented in the previous section indicate that: (1) the regional patterns of  ) is opposite to that in the transition years, consistent with the difference in zonal gradient of underlying SST (Fig. 1i).We assume that the summertime SST forcing and associated condensation heating over the Maritime Continent and the equatorial Pacific are also important for the circulation and AR frequency over the western North Pacific and East Asia in addition to the preceding-winter ENSO forcing (Kamae et al. 2017b).To confirm this hypothesis, we conducted numerical simulations using the LBM (section 2.4).

Atmospheric response to tropical forcing
In this subsection, we explore the possible importance of tropical diabatic heating to the circulation and AR frequency anomaly examined in the previous sections.We use the LBM to examine a linear atmospheric response to the anomalous atmospheric heating in the El Niño transition years relative to the climatology.As shown in the previous section, tropical rainfall exhibits a zonal pair of positive and negative anomalies over the Maritime Continent and western-to-central equatorial Pacific, distinct from what occurs in the sustained years (Fig. 6g, i, Supplement 1).The anomalous convective heating/cooling are prescribed in the LBM.Background atmospheric state is a JJA climatology obtained from the NCEP/NCAR reanalysis (section 2.4).The response at day 20 when the model reaches a quasi-steady state is analyzed.
Figure 7 shows the steady atmospheric response to the tropical diabatic heating.We conducted two experiments: first, the model was forced by a pair of diabatic cooling and heating, which have a deep vertical structure with its rate peaking at 454 hPa (-1 K day -1 and 1 K day -1 ), centered at 180°E, 0°N and 120°E, 0°N, respectively; second, the model was forced solely by the equatorial Pacific cooling.The vertical structure is consistent with the climatological diabatic heating over these regions (Yanai and Tomita 1998).We can examine effects of the two (i.e., the Maritime Continent heating and the Pacific cooling) by comparing the two experiments.The equatorial heating/cooling pair favor a positive WNPSH anomaly (Fig. 7a), largely consistent with the observed and simulated WNPSH anomaly in the El Niño-to-La Niña transition years (Fig. 4a, c).The equatorial heating over the Maritime Continent favors eastward propagating Kelvin waves over the western-to-central equatorial Pacific, similar to the atmospheric response to the Indian Ocean warming (Xie et al. 2009; see section 1).The intensified northeasterly trades over the central tropical Pacific contribute to the enhancement of the WNPSH.
When only the equatorial Pacific cooling is prescribed (Fig. 7b), a pair of zonally-elongated low-level anticyclonic circulation anomalies is favored over the western tropical Pacific to the north and south of the equator.This pair of anticyclones can be simply understood as a Rossby response to the equatorial cooling (Gill 1980).The WNPSH anomaly found in this experiment is consistent with that observed in years with cool SST anomaly over the central Pacific (Wang et al. 2013;Paek et al. 2019).However, the positive response in SLP found here is confined to a narrower meridional extent (0°-30°N) compared to the heating/cooling pair experiment (Fig. 7a).The northern edge of the WNPSH anomaly (Fig. 7b) is not consistent with those found in Figs.4a (48°N) and 4c (40°N), suggesting that the Pacific cooling is not sufficient to explain the circulation and AR frequency anomalies over the NEA found in the El Niño transition years (Fig. 4a, c).Instead, both the heating over the Maritime Continent and the cooling over the Pacific Ocean are important, as demonstrated earlier (Fig. 7a).

Summary and discussion
We examined interannual variability in East Asian AR frequency and associated The latitudinal positions of the WNPSH anomaly in JRA-55 and d4PDF are significantly different, resulting in difficulties in an accurate prediction of regional AR-related disaster risk.One possible factor for the observation-model inconsistency is the difference in sample size.The composite analyses in this study are based on limited samples, for example cases of sustained La Niña year (1971), sustained El Niño years (1983, 1987), La Niña-to-El Niño transition year (1976) and El Niño-to-La Niña transition years (1973, 2010).
Because of the impact of atmospheric internal variability, it is difficult to detect SST-forced signals over the mid-latitudes, especially from the atmospheric reanalysis (Kamae et al. 2017a).Even in the ensemble simulations, the sample size still prevents from obtaining robust atmospheric responses to the ENSO forcing (Fig. 5).To reduce the effects other than the ENSO-related SST (section 3.2), larger samples obtained from longer integrations using coupled atmosphere-ocean general circulation models (e.g.Kay et al. 2016) may be more favorable.
Another important factor for the observation-model inconsistency is the inherent deficiencies in fixed-SST simulations.In the AGCM simulations, prescribed SST boundary condition may result in unrealistic atmospheric responses.Wang et al. (2005) showed that the atmosphere-ocean coupling over the western North Pacific found in the observations is not reproduced in fixed-SST simulations by an AGCM.Zhou et al. (2009)       and solid contours indicate prescribed cooling (-0.6 and -0.2 K day -1 ) and heating (0.2 and 0.6 K day -1 ) at 454 hPa, respectively.

Figure 3
Figure3shows the results for JRA-55.As shown inKamae et al. (2017b), wintertime ENSO atmospheric circulation.During La Niña summers, easterly, westerly, and high-pressure anomalies are found over 20°-30°N, 40°-60°N, and 30°-50°N in the North Pacific, respectively(Figs.2d, g, 3d, g).These are consistent with previous studies on the North Pacific SLP response to the summertime equatorial Pacific forcing(Lau et al. 2005;Paek et al. 2019;Kamae et al. 2019).The central-to-eastern equatorial Pacific cooling related to La Niña (Fig.1a, d, g) suppresses convective heating aloft, favoring the negative phase of a tropics-to-midlatitudes atmospheric teleconnection called the Pacific/North America (PNA) pattern(Wallace and Gutzler 1981).Through this PNA teleconnection, positive SLP anomaly is consistently found over the North Pacific, indicating an intensification of the North Pacific High (NPH; Figs.2d, g, 3d, g).The enhanced NPH results in an increase and decrease in AR frequency over 40°-60°N and 20°-40°N, respectively.Opposite anomalies are also found during the El Niño summers: weakened NPH, a decrease and an increase in AR frequency over respectively (Figs. 2f, i, 3f, i).Exceptions are the AR frequency and atmospheric circulation in a sustained La Niña year(Figs.2a, 3a)      and a La Niña-to-El Niño transition year(Figs.2c East Asia are consistent with the WNPSH anomalies.The positive anomalies in AR frequency on the northwestern flank of the WNPSH is larger in the sustained years (Fig.4b, d) than the transition years (Fig.4a, c), consistent with the stronger WNPSH anomaly during the sustained years.Compared to the sustained years, the increase in AR frequency during the transition years are broader in space over northern part of East Asia (e.g.northeastern China, northern Japan, and around the Sea of Okhotsk; Fig.4a, c), consistent with the meridional extent of the anomaly in the WNPSH.AR frequency averaged over northern East Asia (NEA; 110°-155°E; 40°-55°N, red dashed rectangle in Fig.4) is shown in Fig.5.The NEA region is defined as the northern boundary of the climatological AR occurrence over summertime East Asia (Fig.3cinKamae et al. 2017b).The seasonal-mean ARs over the NEA in the transition years are more frequent (especially in 2010) than the climatology(Figs.4a, c, 5a, b).In contrast, AR frequency over the NEA is greatly reduced in the sustained years.
Kamae     et al. (2017a)  showed that the relative importance of atmospheric internal variability compared to the forced atmospheric response to global SST perturbation is larger over the mid-and high-latitudes than the tropics (Fig.10a, c ofKamae et al. 2017a).The d4PDF ensemble mean can be considered as a forced response to fixed boundary conditions because the effects of atmospheric internal variability cancel out each other(Kamae et al.     2017a, b, Ueda et al. 2018).However, the ensemble mean still exhibits large year-to-year variations under the similar ENSO transition, especially between 1973 and 2010 (Fig.5b).The large variations suggest important effects of factors other than equatorial Pacific SST, for example, Indian Ocean SST, Atlantic SST, North Pacific SST, anthropogenic aerosol emissions, or volcanic eruptions.As an example, variability in the Atlantic SST possibly affects East Asian atmospheric circulation via changing tropical atmospheric circulation (e.g.Li et al. 2016;Kamae et al. 2017d) or mid-latitude circulation pattern (e.g.Liu et al. 2019).In addition, decadal-to-multidecadal variabilities may also affect AR activity discussed here because of their influence on mid-latitude atmospheric circulation (e.g.Newman et al. 2016;     Kamae et al. 2017a; Tokinaga et al. 2019).Environmental factors responsible for the anomalies in AR activity in individual years will be further examined in future studies.
AR frequency over East Asia are different between El Niño transition years and sustained El Niño years; and (2) low-level atmospheric circulation is the key to the difference in AR frequency.The meridional patterns of AR frequency anomaly (broader positive anomaly in the transition years and south-positive north-negative pattern in the sustained years) found in JRA-55 (Fig.4a, b) are largely reproduced in the model simulations (Fig.4c, d).Note that the meridional positions of the AR frequency anomalies are significantly different between the simulations and observations (see section 5).In this section, possible importance of tropical SST forcing and associated diabatic heating to the differences in ARs and atmospheric circulation between transition and sustained years are examined.

Figure 6
Figure6shows summertime rainfall anomaly in the nine ENSO group years relative to summertime atmospheric circulation with a focus on influence of seasonal transition of ENSO from preceding winter to summer.Using JRA-55 reanalysis and d4PDF ensemble simulations, we showed that rapid transitions from wintertime El Niño to summertime La Niña increase the occurrence of ARs over NEA but do not significantly change total occurrence over East Asia compared to sustained El Niño years.In addition to the wintertime El Niño forcing, summertime cool SST over the central-to-eastern Pacific related to La Niña results in strong zonal SST gradient and associated heating gradient in the tropical atmosphere.Linear steady response of the atmosphere to the zonal heating gradient is characterized by a meridionally-elongated anticyclone over the subtropical western North Pacific, consistent with the atmospheric circulation anomaly found in the composite analysis.The combination of the WNPSH anomaly rooted from the preceding-winter El Niño and the meridionally-elongated anticyclonic anomaly favored by the summertime La Niña is effective in increasing AR frequency over the NEA, suggesting the importance of a reliable prediction of ENSO seasonal transitions for the prediction of AR-related disaster risk over the NEA.

Fig. 1 Fig. 2 2 Fig. 3
Fig. 1 Composited sea surface temperature (SST; K) anomaly during June-to-August (JJA) compared to the climatology for 1958-2010 in (a) years with preceding-winter La Niña and concurrent La Niña (La Niña to La Niña; indicated as "LtoL" in the subplot title; Table 1), (b) decaying La Niña (La Niña to Neutral) years, (c) years with a transition from La Niña to El Niño, (d) developing La Niña (Neutral to La Niña) years, (e) neutral years, (f) developing El Niño (Neutral to El Niño) years, (g) years with a transition from El Niño to La Niña, (h) decaying El Niño (El Niño to Neutral) years, and (i) sustained El Niño (El Niño to El Niño) years.Stipples indicate regions with anomalies statistically significant at 0.05 level.

Fig. 7 .
Fig. 7. Atmospheric response to tropical diabatic heating simulated with Linear Baroclinic Model.(a) Atmospheric response to the Maritime Continent heating and the Pacific cooling found in El Niño-to-La Niña transition years (Fig. 4g).(b) Atmospheric response to the Pacific cooling.Shading indicate stream function at 850 hPa (10 6 m 2 s -1 ).Dashed Haarsma et al. 2018)that AGCMs have limited skills in reproducing interannual variation of atmospheric circulation over the extratropical East Asia.The AGCM-based results shown in the present study are consistent with reanalysis in terms of northward expansion of the WNPSH and NEA AR frequency, but other regional features are inconsistent.We plan to conduct further tests using high-resolution atmosphere-ocean coupled model simulations (e.g.,Haarsma et al. 2018)to better understand of the influence of ENSO transition on East-Asian AR activity.