HTML
--> --> -->By applying empirical orthogonal teleconnection (EOT) analysis to observed 1982-2007 pentad rainfall anomalies, (Stephan et al., 2017a) found that 43% of the intraseasonal rainfall variability in China during extended winter is explained by three spatial patterns of temporally coherent rainfall. The leading pattern (EOT-1) is centered along the Yangtze River and was associated with a wave train of Atlantic origin with significant anomalies two pentads in advance. Increased winter precipitation along the coast of Southeast China (EOT-2) was connected to an intensifying low pressure anomaly over Asia and wave propagation from the Mediterranean Sea. Intraseasonal rainfall variability between the Yangtze and Huaihe rivers (EOT-3) was associated with a quasi-stationary wave train that could be traced back to a strong anomalous ridge over northern Europe two pentads in advance. (Stephan et al., 2017a) also performed an analogous analysis for summer. Summer patterns were predominantly associated with tropical convection instead of extratropical wave patterns. Therefore, this study focuses on winter.
Links between subseasonal weather extremes in China and extratropical disturbances have also been reported by others. Cold surges are associated with low air temperatures and strong winds in the midlatitudes and subtropics (Boyle and Chen, 1987). Using empirical orthogonal function analysis of pentad precipitation data, (Yao et al., 2015) found that the two dominant patterns of subseasonal precipitation variability in extended winter over China consist of a monopole in South China and a meridional dipole with opposite-signed precipitation anomalies over the Yangtze River basin and the southern coast. (Yao et al., 2015) linked both patterns to a southward propagating cold surge, triggered by a wave train from the North Atlantic. In their composite analysis, (Takaya and Nakamura, 2005) also found the intraseasonal amplification of the Siberian high to be associated with a blocking ridge that is part of a quasi-stationary Rossby wave train propagating across the Eurasian continent. In a case study of the 2005/06 winter, (Park et al., 2008) examined the causes of two consecutive cold surges. One was related to North Pacific upper-level blocking, and the other to an upper-level wave train across the Eurasian continent. Later, (Park et al., 2011) grouped cold surges into "wave train" and "blocking" types based on circulation features, and discussed their relationship with the Arctic Oscillation (AO; Thompson and Wallace, 1998). They found the blocking type was often associated with the negative AO phase, but no phase preference was found for the wave train type. (Park et al., 2014) further examined the origins of the wave-train cold surge. At a lead time of approximately 12 days they reported negative upper tropospheric height anomalies southeast of Greenland; they suggested that these anomalies may originate in the lower stratosphere over the North Atlantic. (Yao et al., 2015) and (Stephan et al., 2017a) could not find evidence for stratospheric origins of precursor wave trains.
The above indicates that precursor circulations of winter rainfall events become statistically significant at lead times of approximately 10-12 days in regression or composite analyses. The physical nature of the wave trains and their origins remain unclear. Furthermore, previous research has not addressed whether it is possible to clearly identify these precursors for individual extreme events. If such precursors could be identified they could be used for an empirical prediction of rainfall extremes. (Stephan et al., 2017a) and some of the other above-mentioned studies were based entirely on observations and reanalysis. In the present study, we repeat the (Stephan et al., 2017a) EOT analysis on six simulations of the Met Office Unified Model (MetUM) Global Atmosphere 6.0 and Global Coupled 2.0 configurations at resolutions of ～200, 90 and 40 km (in the zonal direction at the equator). This allows us to evaluate the model in terms of its ability to reproduce observed patterns of winter rainfall and to test the sensitivity of precursors to horizontal resolution and air-sea coupling. Secondly, by comparing synthetic time series to observations, we can assess the robustness of precursor circulations.
Section 2 introduces the simulations, data and methods. In section 3 we examine how well the simulations reproduce the observed EOT patterns. The observed and simulated precursors are compared in section 4. Individual events are discussed in section 5. Section 6 examines the dynamics of individual events. Section 7 addresses the wave origins and wave dynamics. A summary is given in section 8.
2.1. MetUM simulations
We analyze two 27-year atmosphere-only simulations of the MetUM Global Atmosphere configuration 6.0 (GA6; Walters et al., 2017) and four 100-year coupled simulations of the Global Coupled configuration 2.0 (GC2; Williams et al., 2015). If the model validates well, long simulations would provide synthetic catalogues of events longer than observations, for analysis of extremes and variability with a larger sample size. We refer to them as A96, A216, C96, C216, C512a and C512b, where "A" and "C" stand for "atmosphere-only" and "coupled", respectively, followed by the nodal number (N96: 1.875°× 1.25°, 208 km ×139 km in longitude and latitude at the equator; N216: 0.83°× 0.55°, 93 km × 62 km; N512: 0.35°× 0.23°, 39 km × 26 km). A96 and A216 use historical forcing; the GC2 simulations are present-day control simulations. Key information about the simulations is summarized in Table 1. A more detailed description of each simulation is given in (Stephan et al., 2017b).2
2.2. Observational data
Precipitation data are obtained from the Asian Precipitation-Highly Resolved Observational Data Integration Toward Evaluation of Water Resources (APHRODITE) dataset (Yatagai et al., 2012). This continental-scale daily product, from 1951-2007, with a resolution of 0.5°× 0.5°, is produced from rain-gauge data after applying an objective quality control procedure (Hamada et al., 2011). We use 500 hPa geopotential height (Z500) for 1979-2007 from the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-Interim; Dee et al., 2011). The resolution of these data is 0.7°× 0.7°. Additional diagnostics are computed from ERA-Interim potential vorticity fields on isothermal surfaces and 200 hPa horizontal wind fields. For simplicity we refer to ERA-Interim as observations.2
2.3. EOTs
The EOT analysis of 1951-2007 gridded APHRODITE precipitation data by (Stephan et al., 2017a) serves as the basis for this study. EOT analysis extracts spatial patterns of temporally coherent precipitation variability and returns time series that are mutually orthogonal (Smith, 2004). Such variability is likely to have the greatest impact on infrastructure and human life due to its spatial organization. The steps of the algorithm are detailed in (Stephan et al., 2017a). Before applying the technique to simulations, MetUM precipitation data are interpolated to the APHRODITE grid. We use linear regression onto EOT time series to link individual patterns to associated atmospheric precursors. To account for the non-Gaussian distribution of rainfall data, we use Spearman's rank correlations to determine the statistical significance (always at the 10% significance level).2
2.4. Rossby wave source function
To diagnose the propagation of Rossby waves we compute the Rossby wave source (RWS) function as defined by (Sardeshmukh and Hoskins, 1988) in their Eq. (3). The 200 hPa absolute vorticity η, horizontal wind divergence D, and the divergent horizontal wind vector vχ are combined to form \begin{equation} RWS=-\eta D-{v}_\chi\cdot\nabla\eta .\ \ (1) \end{equation} The first term describes divergence or convergence, respectively, above ascending or descending air, i.e., the vorticity tendency arising from vortex stretching. The second term describes advection and becomes large where divergent wind encounters sharp horizontal vorticity gradients. The vertical profile of the RWS peaks at 200 hPa because vorticity gradients associated with a vertical peak in the jet stream and divergent horizontal winds associated with convective outflow maximize at that level (Scaife et al., 2017).2
2.5. Rossby wave ray tracing
We perform backward Rossby wave ray tracing to identify regions where waves may have originated. Given a slowly varying zonal flow with velocity U, the dispersion relation of a barotropic Rossby wave is \begin{equation} \omega=Uk-\dfrac{\beta^*k}{K^2} .\ \ (2) \end{equation} Here, ω is the frequency; k and l are the zonal and meridional wavenumbers, respectively; \(K=\sqrt{(k^2+l^2)}\) is the total wavenumber; and β*=β-Uyy is the meridional gradient of the absolute vorticity of the mean flow, which combines the gradient in planetary vorticity β and the curvature of the flow Uyy, i.e., the second meridional derivative of U (Hoskins and Karoly, 1981).For a stationary wave, ω=0, and Eq. (2) becomes \begin{equation} K^{2}=\dfrac{\beta^*}{U} .\ \ (3) \end{equation} Hence, the zonal and meridional group velocities of a stationary wave are given by \begin{equation} c_{\rm gx}=\dfrac{\partial\omega}{\partial k}=\dfrac{2U^2k^2}{\beta^*} \ \ (4)\end{equation} and \begin{equation} c_{\rm gy}=\dfrac{\partial\omega}{\partial l}=\dfrac{2U^2kl}{\beta^*} , \ \ (5)\end{equation} respectively. To trace Rossby waves of a given zonal wavenumber k, we first compute l from Eq. (2). The previous location of the wave front is found from Eqs. (4) and (5), taking into account the spherical geometry of the globe. Points where U=0 form critical lines where the propagation of Rossby waves is not supported; the ray ends. A ray reverses its meridional propagation direction when l approaches zero and thus K2=k2.
The background wind U is the November-April mean climatological 200 hPa zonal wind. The background wind is smoothed using a 60° zonal average. The curvature term Uyy is smoothed using a full 360° zonal average. Our choices of pressure level, a two-hour time step, and smoothing of the background wind field, follow the recommendations in (Scaife et al., 2017). They showed that rays are rather insensitive to the curvature term and time step, but that they are clearly affected by the choice of pressure level and the degree of smoothing of the background wind. We do not perform ray tracing on simulated wind data because 200 hPa wind biases are small and would not significantly alter rays.
2
2.6. Rossby wave initiation identification
Rossby wave initiation (RWI) segments are computed using the (R?thlisberger et al., 2016) algorithm. (R?thlisberger et al., 2016) discussed in detail the choice of the algorithm's tuning parameters. We use their recommended parameters and here only briefly describe the algorithm. The algorithm extracts the geometry of isentropic 2 potential vorticity unit (PVU, 1 PVU = 10-6 m2 s-1 K kg-1) contours, corresponding to the dynamical tropopause. For isentropic levels of 310-320 K (340-350 K), the contour is approximately aligned with the extratropical (subtropical) jet and can be used to measure its waviness (Hoskins et al., 1985; Martius et al., 2010). Every six hours and for 60° longitudinal contour segments, starting every 3° of longitude, the waviness d is measured by integrating the absolute latitudinal variations of the contour position over the length of the segment. Hence, a zonally aligned, straight jet corresponds to small d. The algorithm identifies wave-free segments (d<20°) that show a strong increase (>8°) in d within 30 hours. Candidate segments where increases in waviness may result from downstream development of existing waves are discarded if the waviness in one of the segments starting at 45° longitude upstream of the candidate segment exceeds the waviness in the candidate segment by more than 4°.
Figure1. (a) Observed and (b-g) simulated EOT-1 patterns. Shading shows regressions of November-April precipitation against the normalized EOT time series. Also shown are correlations of the full precipitation-anomaly time series with the EOT base point exceeding 0.8 (magenta), 0.6 (orange), and 0.4 (green). The EOT base point is marked by the orange inverted triangle.Simulated rainfall anomalies match Obs-1, with pattern correlation coefficients exceeding 0.92 (Table 2). The fraction of explained variance in the simulations (17%-25%) is close to the observed value (21%). The standard deviations of the simulated EOT time series are 25%-100% larger than in observations. This is consistent with positive biases in simulated intraseasonal variability and mean precipitation over Southeast China (not shown). Similar results are found for Obs-2 (Fig. 2).
Obs-3 is situated further north (Fig. 3), where biases in winter intraseasonal variability are small (not shown). Here, the simulated amplitude agrees with observations (Table 2). In the atmosphere-only simulations, A96 and A216, areas of covariability extend further northwest; in South China there is a large region with variability in the opposite phase. This explains why the fraction of explained variance for EOT-3 in A96 and A216 is higher (10% and 9%) than observed (7%). Pattern correlations exceed 0.67 for all simulations.
Figure2. (a) Observed and (b-g) simulated EOT-2 patterns. Shading shows regressions of November-April precipitation against the normalized EOT time series. Also shown are correlations of the residual precipitation-anomaly time series with the EOT base point exceeding 0.8 (magenta), 0.6 (orange), and 0.4 (green). The EOT base point is marked by the orange inverted triangle.
Figure3. (a) Observed and (b-g) simulated EOT-3 patterns. Shading shows regressions of November-April precipitation against the normalized EOT time series. Also shown are correlations of the residual precipitation-anomaly time series with the EOT base point exceeding 0.8 (magenta), 0.6 (orange), and 0.4 (green). The EOT base point is marked by the orange inverted triangle.Obs-1 is associated with high Z500 over the eastern United States and central Europe, and low Z500 over southern Greenland, northern Africa and Siberia at Tp=2. At Tp=1, there is a wavenumber-2 pattern of Z500 anomalies between the North Atlantic and eastern Asia. At Tp=0, a zonal dipole is present over Asia and the western North Pacific. The RWS function indicates the genesis of upper tropospheric vorticity anomalies to the east of each of the above-mentioned Z500 anomalies, indicating eastward propagation from the Atlantic sector through Europe to Asia. At Tp=3, observations show only small regions with significant Z500 anomalies: high Z500 over the northeastern US and the North Atlantic, and low Z500 over southern Greenland. All simulations reproduce the Z500 anomalies listed above and the anomalies are statistically significant, even at Tp=3.
Figure4. Observed and simulated Northern Hemisphere precursors of EOT events at lead times of 3-0 pentads (left to right). Contours show Z500 regressed against the normalized observed and C512a EOT time series (blue, negative; orange, positive; intervals of 4 m for EOT-1 and EOT-3; 1 m for EOT-2). The EOT base point is marked by the black inverted triangle. For observations, the regression of the RWS function (shading) is shown. All values are significant at the 90% confidence level. Green areas in the C512a panels show areas where all six simulations agree on the sign of the Z500 regression slope.Obs-2 is associated with low Z500 in the Atlantic, northern Africa and South Asia, and high Z500 over eastern Canada at Tp=3. At Tp=[2,1] there is low Z500 over middle America, the Atlantic, and from northern Africa to South Asia, and high Z500 at high latitudes between eastern Russia and eastern Canada and over western northern Africa. At Tp=0, there are two connected low Z500 areas over the Middle East and East Asia. Prominent RWS anomalies are found in the Atlantic (Tp=3) and across northern Africa (Tp=[2,1]). Simulations produce the same low-latitude anomalies at Tp=[2,1,0]. At Tp=3 they agree with observations on the low Z500 in the Gulf of Mexico and the Atlantic. High Z500 at high latitudes forms later in the simulations (Tp=1 instead of Tp=3), but so do the low Z500 areas over northern Africa and South Asia (Tp=2 instead of Tp=3). This suggests that high-latitude anomalies develop in response to low-latitude anomalies, and that tropical anomalies are the relevant precursors of EOT-2 events in China.
Obs-3 is characterized by high Z500 over northern Europe at Tp=[2,1,0]. At Tp=[2,1], there is low Z500 over northern Canada and Greenland. From Tp=2, low Z500 areas also develop between northern Africa and Siberia, similar to Obs-1. A tripole of high Z500 over northeastern Europe, low Z500 over western Asia, and high Z500 over East Asia, is present at Tp=[1,0]. The strongest RWSs are found to the north of the Arabian Peninsula at Tp=1. All simulations agree on the development of the tripole at Tp=[2,1,0]. Neither observations nor simulations show noteworthy anomalies at Tp=3.
The agreement between observed and simulated anomalies preceding precipitation associated with the leading three EOTs is remarkable and points to the existence of robust physical mechanisms or common statistical artifacts from regression analysis. Below we investigate whether the common anomalies discussed above may be used as potential predictors.
Figure5. Characteristics of observed and simulated top EOT events. (a-c) Monthly occurrence of the top 1/6 EOT events. The color coding is the same as in (d-f). (d-f) Top 10 unmixed EOT events observed in FMA (S10): five-day mean anomaly relative to the pentad climatology averaged over the area where correlations of the full (leading order) or residual (higher order) precipitation-anomaly time series with the EOT base point exceed 0.5.All simulations capture the seasonality of top events and the magnitudes of S10 sufficiently well for the purpose of comparing simulated and observed events and their precursors. Figures 6a-c show maps of precipitation anomalies averaged over S10. The black contour marks the region where correlations of the full or residual precipitation-anomaly time series with the EOT base point exceed 0.5. Consistent with the values in Figs. 5d-f, rainfall anomalies inside this region are ～100%-300%. From the correlations shown in Figs. 1-3, it is expected that the spatial distribution of rainfall differs from event to event. To measure the similarity of S10 rainfall, Figs. 6d-f show the number of events for which anomalies exceeded 100%. Inside the 0.5 correlation region, this number varies between 4 and 10 for EOT-1 and EOT-3; for EOT-2, it exceeds 7 for most points. Thus, Figs. 5d-f and Fig. 6 show that S10 contains a self-consistent sample of EOT events. Figure 6 was also computed for the six simulations to verify that their S10 sets are self-consistent as well (not shown). The linear pattern correlation coefficients of Figs. 6d-f and the corresponding maps for simulations vary between 0.26 and 0.52 for EOT-1, between 0.43 and 0.57 for EOT-2, and between 0.46 and 0.56 for EOT-3. In addition, we examined each event separately to confirm the consistency.
Figure6. Precipitation characteristics of the top 10 unmixed EOT events observed in FMA (S10). (a-c) Average five-day mean precipitation anomaly relative to the pentad climatology. The EOT base point is marked by the white inverted triangle. The black contour marks the region where correlations of the full (leading order) or residual (higher order) precipitation-anomaly time series with the EOT base point exceed 0.5. (d-f) Number of EOT events for which average five-day mean precipitation anomalies exceed 100%.
Figure7. Inter-event variability in the evolution of Z500 for the top 10 unmixed EOT events in FMA (S10) in (a-c) observations and (d-i) two simulations. The upper part of each panel shows the composite of Z500 at the designated lead time; the lower part shows the number of events where Z500 is of the same sign.We next investigate similarities in the precursors associated with S10 for EOT-1. Figures 7a-c are composites of observed Z500 together with maps that show the number of events with Z500 anomalies of the same sign. Z500 anomalies are relative to 1982-2008 pentad climatologies. At Tp=2, we recognize several anomalies that are significant in the regression in Fig. 4: high Z500 over Scandinavia, and low Z500 over northwestern Africa and Siberia and the North Atlantic. Except for the high Z500 over Europe, these anomalies are seen in eight of ten events. At Tp=1, eight events agree on the tripolar Z500 anomaly that extends from central Europe into the Arabian Sea. Mid- and high-latitude anomalies are shifted compared to the regression. In the simultaneous pentad, Z500 anomalies are weak and differ greatly amongst S10.
The same information is given for C512a (Figs. 7d-f) and C512b (Figs. 7g-i). At Tp=2, nine or ten events in C512a exhibit low Z500 to the west of Greenland and high Z500 over northern Europe. The above-mentioned tripolar Z500 low-latitude anomaly is also present in the composite at Tp=[2,1,0]. Unlike observations, all C512a events show high Z500 over the western North Pacific at Tp=0. C512b composites differ substantially from C512a. Similar to observations, they show only weak anomalies over the North Atlantic and Europe. However, C512b agrees on the low-latitude tripole at Tp=2 and all events exhibit a zonal dipole over Asia and the western North Pacific at Tp=0. We chose to show C512a and C512b because their 100-year long records allow us to best isolate strong and unmixed EOT events. Furthermore, their comparison lets us conclude that any agreement in C512a or C512b is not related to biases associated with the C512 configuration; instead, it is random.
If we assumed that in each event and at each location negative and positive Z500 anomalies occurred with equal likelihood, the chances that five, six, seven, eight, nine, and all 10 events agreed on the sign are 0.25, 0.41, 0.23, 0.09, 0.02 and <0.005, respectively. Averaged over observations, C512a and C512b, the fractions of grid points with this agreement are 0.25, 0.42, 0.23, 0.08, 0.02 and <0.005. This confirms that individual EOT-1 events show little agreement with the significant Z500 anomalies of the regression. The dynamical evolution of the atmosphere is very different from event to event and does not resemble the slowly evolving or stationary features of the regression. From the comparison of observations, C512a and C512b, we conclude that there are no robust indicators of EOT-1 events that could serve as useful predictors at any lead time. These conclusions hold for the other four simulations and for EOT-2 and EOT-3. Aside from the Z500 evolution, we also examined 200 hPa geopotential height, surface pressure, 200 hPa wind, the AO (Thompson and Wallace, 1998), North Atlantic Oscillation (Hurrell et al., 2003), Madden-Julian Oscillation (Madden and Julian, 1972) and atmospheric blocking computed as in (Henderson et al., 2016), but could not find significant or systematic indicators for any EOT. Despite this great inter-event variability, our set of MetUM simulations still accurately reproduces the observed EOT rainfall patterns and the observed Z500 evolution for the overall regressions. Next, we consider whether there is a physical basis for this agreement.
At Tp=2, a pattern of high Z500 in the North Pacific, low Z500 over southern Greenland, and high Z500 over central Europe is consistent with k=2 waves with possible origins in a broad region over the Pacific (Fig. 8a). After 5-10 days they reach the Middle East and Asia, such that they may also contribute to the southwest-northeast tilted low Z500 pattern at Tp=1, and to the high Z500 over the Arabian Sea (the target region). The tripole between Europe and the Arabian Sea is consistent with k=3 propagation (Fig. 8b). Note that these waves propagate from Europe to the Arabian Sea within only two to three days and could originate at various longitudes over the Atlantic basin. The EOT-2 target region represents the location of the eastern low Z500 center in the regression at Tp=1. Rossby waves with k=3 may originate over Africa, central Europe or the North Atlantic (Fig. 8c). All waves propagate across the western node of the low Z500 anomaly. Frequent RWIs over Africa are consistent with the strong RWS anomaly in the regression. Recall that EOT-2 had statistically significant Z500 and RWS anomalies over the Atlantic at Tp=3 (Fig. 4). The number of RWIs over the Atlantic is substantially higher than for the other EOTs (Fig. 8d), providing further evidence that processes over the Atlantic may be influencing conditions over Africa in the following pentads. Low Z500 over the target region is already present at Tp=3. Aside from the k=3 pathway, this may be communicated by k=1 waves originating over the western Atlantic (Fig. 8d).
Figure8. Wave dynamics preceding EOT events: regression of Z500 against normalized observed EOT time series (shading; significant at the 90% confidence level); 60°-longitude RWI segments (horizontal lines) on the extratropical (blue) and subtropical (red) jet for the top 100 unmixed EOT events; Rossby wave rays (dots) traced backwards in time from points inside the target region (black box) for 0-5 days (bright brown) and 5-10 days (dark brown). For unmixed EOT events, the considered EOT time series must exceed one standard deviation, while the other two EOT time series must have amplitudes smaller than one standard deviation. RWI segments are drawn at the mean latitude of the RWI segment.The EOT-3 target region extends from North China to Japan (Fig. 8e). At Tp=2, Z500 anomalies are well explained by k=1 waves. These contribute to the northern portion of the high Z500 anomaly over the Barents Sea, the northern portion of the low Z500 anomaly over Siberia, and the anomaly inside the target region itself. At Tp=1, k=2 waves with a large variety of origins are able to propagate into the target region. Possible source regions include the North Atlantic, the central North Pacific, the target region itself, eastern Russia, and North Africa. Given that all endpoints inside the target region are equally spaced in longitude and latitude and otherwise chosen at random, it is remarkable that all rays intersect significant Z500 anomalies in these source regions (Fig. 8f). These rays pass through the southern portion of high Z500 over Europe, the central part of low Z500 over Kazakhstan, and the high Z500 in the target region.
The above shows that waves with different wave numbers and various origins pass through locations of convergent waveguides. These locations are associated with the strongest Z500 anomalies in the regressions. Weaker Z500 anomalies are found in locations that are less frequently crossed by waves associated with EOT events. Rather than coherent waves, the precursors in regressions are statistical amalgamations of more rapidly propagating waves with a variety of origins and properties. The 15-20-day timescale of precursors in regression analysis is consistent with the period of low-frequency modulations and amplifications of alternating hemispheric flow regimes found in previous modeling and observational studies (e.g., Wu, 1993; Petoukhov et al., 2016).
Here, the aim was to better understand the dynamical evolution and the origins of the wave patterns found by (Stephan et al., 2017a), to investigate whether precursors found in regression analysis may be used as predictors of strong extended winter rainfall events in China. We repeated the (Stephan et al., 2017a) EOT analysis on six simulations of the MetUM, using atmosphere-only and coupled configurations at resolutions of ～200, 90 and 40 km (in the zonal direction at the equator). There is excellent agreement in all simulations with the observed EOT patterns and precursor circulation anomalies. Hence, there is no evidence that the physical processes associated with EOT rainfall patterns are sensitive to air-sea coupling or to horizontal resolution. The robust simulation of extratropical wave disturbances would suggest the possibility for empirical prediction of EOT rainfall variability. However, in examining individual severe observed and simulated events, we found that the dynamical evolution of the atmosphere differs severely from event to event and does not resemble the slowly evolving or stationary features of the regression. In individual events, the most prominent anomalies in the regression can be missing, or of opposite sign. There are no robust indicators of EOT events that could serve as useful predictors at any lead time. Rossby wave dynamics are able to explain how the large inter-event variability is consistent with what appear to be slowly evolving or stationary waves. Using diagnostics based on observations and theoretical arguments, we showed that Rossby waves with a large variety of origins, wave numbers and group velocities can contribute to pressure anomalies over East Asia. Along their trajectories these waves produce common Z500 anomalies that appear in composite or regression analysis as slowly evolving features. It is possible that this conclusion also holds for other extratropical teleconnection patterns. While anomalies in composite or regression analysis may not usefully predict rainfall, it is possible that EOT patterns are modulated in a predictable way by the slowly changing atmosphere-ocean coupled background state. Exploring this topic further could be an interesting avenue for future research.
