Comparison of equatorial wave activity in the tropical tropopause layer and stratosphere represented in reanalyses

. Equatorial Kelvin and mixed Rossby-gravity (MRG) waves in the tropical tropopause layer and stratosphere represented in recent reanalyses for the period of 1981–2010 are compared in terms of spectral characteristics, spatial structures, long-term variations and their forcing of the quasi-biennial oscillation (QBO). For both wave types, the spectral distributions are broadly similar among most of the reanalyses, while the peak amplitudes exhibit considerable spread. The longitudinal distributions and spatial patterns of wave perturbations show reasonable agreement between the reanalyses. A few exceptions 5 to the similarity of the spectral shapes and spatial structures among them are also noted. While the interannual variations of wave activity appear to be coherent for both the Kelvin and MRG waves, there is substantial variability in long-term trends among

sources of momentum needed to drive the easterly-to-westerly and westerly-to-easterly phase transitions of the quasi-biennial oscillation (QBO), respectively (Holton and Lindzen, 1972;Dunkerton, 1997;Baldwin et al., 2001), although the momentum transported by the MRG waves needs to be further quantified (Randel et al., 1990;Kim and Chun, 2015a).
Investigation of the distribution and variability of large-scale equatorial waves requires datasets with global coverage, and global (re)analyses are extremely useful datasets for this purpose. A reanalysis is a product of a data assimilation system which 15 reconciles the observed atmospheric states from many kinds of measurements with the atmospheric governing equations resolved by a numerical prediction model. Reanalysis products depend on the assimilation method and prediction ("first guess") model, especially in the stratosphere where the assimilated fields are less constrained by the observations due to less density of observations when compared to the troposphere. Around the TTL, there exist abrupt vertical changes in the temperature and stability, which greatly modify characteristics of the equatorial waves, particularly their vertical wavelengths and ampli-20 tudes (e.g., Randel and Wu, 2005;Ryu et al., 2011). Representation of the TTL in reanalyses might be sensitive to the vertical resolution of the prediction model and assimilation techniques used (see Birner et al., 2006, for a case of the extratropical tropopause). Therefore, it is important to identify the difference/spread among various reanalysis products in their representation of the equatorial waves in the TTL and stratosphere. For example, equatorial wave activity in the TTL during 1990TTL during -1999 was compared using seven reanalyses by Fujiwara et al. (2012), and significant differences between the reanalyses were found. 25 In contrast to analyses, which are made using operationally changing, state-of-the-art versions of the prediction model and assimilation system, reanalyses are derived using fixed versions of the data assimilation system and first guess models for the whole period of the product. This helps to generate a temporally homogeneous product, which benefits studies of long-term changes in meteorological variables. However, reanalyses also have a potential for inhomogeneities or discontinuities arising from introduction of new observational data into the assimilation. Examples include the introduction of radi- 30 ance/temperature profiles derived from the Television Infrared Observation Satellite (TIROS) Operational Vertical Sounder (TOVS) suite (Smith et al., 1979) around 1978 and those from the Advanced TOVS (ATOVS) suite around 1998. While these measurements have improved the quality of reanalyses, it has been reported that the transition from the TOVS to ATOVS suites induced temporal discontinuities in assimilated variables such as upper stratospheric global-mean temperature (Onogi et al., 2007;Simmons et al., 2014) and equatorial stratospheric mean wind and temperature (Kawatani et al., 2016).
To the best of our knowledge, impacts of such satellite transitions in reanalyses have so far not been studied in the context of equatorial waves. It is likely more difficult to identify discontinuities or inhomogeneities in wave fields which, by definition, contain smaller-scale variations in space and time than mean fields. Recently, a pair of reanalyses have been identically pro-5 duced using a single assimilation system, with the exception that the satellite data are assimilated in one (JRA-55) but not in the other (JRA-55C; see Table 1 for the abbreviations and references). These datasets can allow us to identify the impact of satellite data in that assimilation system and, in particular, to help further distinguish the impact of the TOVS-ATOVS transition.
In this study, we investigate the characteristics of equatorial Kelvin and MRG waves in the TTL and stratosphere and how they differ between recent reanalyses for the period 1981-2010 (Section 3.1). Spatial distributions and patterns of the waves 10 in the reanalyses are presented in Section 3.2. In Section 3.3, long-term changes in equatorial wave amplitudes are compared among the reanalyses and, based on the comparison between JRA-55 and JRA-55C, the effects of the satellite data upon the assimilated waves are discussed. In addition, spectra of the Eliassen-Palm (EP) flux and its divergence, a measure of wavemean flow interaction, are presented to compare the equatorial wave forcing of the QBO estimated from different reanalyses (Section 3.4). Among the results presented in Section 3.1, we identify a wave signal that has not been dealt with in the literature 15 before, found at 20 hPa along with the well-documented equatorial waves in all the reanalyses studied. This wave spectrum is further discussed in Section 4. A summary of the results is included in Section 5.

Data and method
We examine six reanalyses: ERA-I, MERRA, MERRA-2, CFSR, JRA-55, and JRA-55C (see Table 1 for their full names, vertical resolutions in the lower stratosphere, and key citations). The horizontal and vertical winds, temperature, and geopotential 20 in the TTL and stratosphere for the period of 1981-2010 are used. The data we use are stored at 3-hour intervals for MERRA and MERRA-2 and 6-hour intervals for the others. The results do not significantly change if we use 6-hourly subsampled time series for MERRA and MERRA-2 (not shown) as we analyze equatorial waves with periods longer than 2 days. Both pressure-level datasets, which are interpolated to standard levels (SL, e.g., 100 hPa), and model-level (ML) datasets are used for each reanalysis, except for MERRA of which ML data are not available. At 100 hPa and above in the tropics, pressure 25 variations on a model level are negligibly small or absent in all of the reanalyses used in the present study (see Supplement to Fujiwara et al., 2017). Therefore, the model levels can be regarded as being at nearly constant pressure levels, which allows us to simply perform spectral calculations on horizontal planes without the need to introduce vertical interpolation. For spectral shapes and spatial distributions of Kelvin and MRG waves, the results were very similar between SL and ML datasets (not shown), whereas wave amplitudes differ significantly as will be seen in Fig. 8. Therefore, only the SL results are presented in Zonal wavenumber-frequency (k-ω) spectra are calculated monthly at each latitude and height for each of the symmetric and antisymmetric components of variables with respect to the equator. To obtain the monthly spectra, we use a 90-day time window centered on the target month. The window function is defined as C for the central 30 days and C sin (πt/60) and C cos (πt/60) for the first and last 30 days, respectively, where t is the time (in days) relative to the first day of each 30-day segment and the normalization constant C = 3/2. The window function is determined such that for a long-term mean, the 5 integral of the power spectrum equals the unwindowed variance of the variable. The spectra are calculated using the Fourier transform after removal of zonal mean and application of the time window, and they are averaged over the latitude range 15 • N-15 • S. The spectra are plotted in a variance-preserving form using base-10 logarithm axes in frequency and wavenumber.
The use of a 90-day time window retains intraseasonal variations such as the MJO partly in the spectra in the TTL at the lower-frequency range, as will be seen in Fig. 1. However, we will exclude these signals from our analysis and focus on signals 10 with periods shorter than 20 days.
Following previous studies (Wheeler and Kiladis, 1999;Hendon and Wheeler, 2008;Fujiwara et al., 2012), the background spectra are obtained using the 1-2-1 filter repeatedly in wavenumber and frequency for each reanalysis. In this study, common background spectra for symmetric and antisymmetric components are obtained by averaging spectra of the two components before applying the filter. The filter is applied to the logarithm of the power spectrum; the number of passes is 23 for zonal 15 wavenumber and 7 for frequency (Fujiwara et al., 2012).
k-ω spectra of the EP flux are also calculated monthly for symmetric and antisymmetric modes (Section 3.4). The EP flux formulation defined with the transformed Eulerian mean of the primitive equations is used (Andrews et al., 1987): 20 where A 1 = ∂u ∂z ∂θ ∂z and A 2 = ζ a ∂θ ∂z .
Here, u, v, w, θ, and ζ a are the zonal, meridional, and vertical winds, potential temperature, and absolute vorticity, respectively; ρ 0 = p/(gH) where p and g are the pressure and gravitational acceleration, respectively, and H = 6.6 km. The overbar and prime denote the zonal average and the perturbation from the average, respectively. The EP flux spectra for symmetric 25 (antisymmetric) modes are obtained using u S , θ S , w S , and v A (u A , θ A , w A , and v S ), where the subscripts S and A denote the symmetric and antisymmetric components of each variable, respectively. The flux terms are calculated via co-spectra (e.g., ℜ{ŵ Sû * S } for the term w ′ u ′ for symmetric modes, where the hat indicates the complex Fourier coefficients of the variables, and the asterisk the complex conjugate). The same window function as the above is applied to all perturbation variables before calculating co-spectra. The terms A 1 and A 2 are averaged for the target month to be multiplied by the co-spectra. The EP 30 flux spectra are averaged over 5 • N-5 • S, considering the meridional structures of the Kelvin and MRG wave EP fluxes and their divergence. The EP flux divergence would be largely underestimated if a broader latitude band is used for averaging, in particular for MRG waves for which the EP flux divergence changes its sign around 5 • (e.g., Kim and Chun, 2015b, their Fig. 9c;Garcia and Richter, 2019, their Figs. 9b and 9d).
The buoyancy frequency used to compute equivalent depths of the equatorial waves is set to 0.024 s −1 based on the climatological temperature profile of the tropical lower stratosphere (e.g., Grise et al., 2010).

5
3.1 Spectral characteristics Figure 1 shows k-ω power spectra for the eastward-propagating, symmetric component of temperature (T S ) at 100, 70, 50, and 20 hPa from the SL datasets of the six reanalyses, averaged over the 30-year period (1981-2010). The major portion of the spectral power appears along the Kelvin wave dispersion curves (black diagonal lines) at all altitudes. The spectral characteristics are broadly similar among the reanalyses with peaks at k = 2-3. The reanalyses commonly exhibit a gradual 10 shift of the Kelvin wave spectrum to larger equivalent depths (h) with respect to the altitude: the majority of the spectral power appears at h < 60 m at 100 hPa, at h ∼ 60 m at 70 hPa, and at h = 60-240 m above 70 hPa. There exist variations in the spectral width in h: for example, compared to the other reanalyses, CFSR shows broader spectra at small equivalent depths (h < 60 m) consistently at 70-20 hPa, with peaks at slightly smaller equivalent depths (70-50 hPa). In addition to the Kelvin wave signal, at 100 hPa, another peak exists at k = 2 in the low-frequency range (ω < 0.04 cyc day −1 ). As mentioned in 15 Section 2, this is likely related to the intraseasonal MJO (Hendon and Wheeler, 2008).
The relative magnitude of Kelvin wave spectral power between ERA-I, MERRA, MERRA-2, and CFSR varies with height ( Fig. 1), while JRA-55 and JRA-55C show generally less power below 20 hPa compared to the other four datasets. Note that relatively small variances in JRA-55 and JRA-55C are found in the temperature field but not in the wind fields (see Fig. S1). In Fig. 1, the ratio of the spectral power to the background spectrum is also shown (thin purple contours), which 20 indicates statistical significance of the spectral signals. A ratio of 1.2 was deemed to be statistically significant at the 95% level by Wheeler and Kiladis (1999), and for the majority portions of the Kelvin wave spectra in Fig. 1, the ratios are generally larger than 1.5 (i.e., 50% larger than the background spectral power), conservatively implying that the spectral peaks are all statistically significant.  Note that the lower bound of y-axis at 100 hPa is different from that at the other levels, as the v S spectrum is much broader in frequency at 100 hPa. The 100-hPa spectrum has periods from around 2.5 days (0.4 cyc day −1 ) to longer than 30 days (0.033 cyc day −1 ), with a peak at k = −5, ω = 0.1-0.2 cyc day −1 common to all of the reanalyses. In the upper troposphere, the background zonal wind near the equator is westerly (easterly) in the western (eastern) hemisphere. MRG waves that are generated in the region of westerly (easterly) background flow can have relatively low (high) ground-based frequencies, which 30 might result in the broad spectrum in frequency of v S at 100 hPa. The low-frequency portion of the wave spectrum at 100 hPa seems to roughly follow the MRG wave dispersion curves for a background wind (U ) of +10 m s −1 (Fig. 2, long-dashed).
At 70 hPa, the low-frequency portion (< 0.1 cyc day −1 ) of the spectrum is mostly filtered out, and commonly in all reanalyses, the peaks appear at k = −6, ω ∼ 0.20 cyc day −1 and at k = −5, ω ∼ 0.22-0.25 cyc day −1 (h ∼ 60 m for both, assuming U = 0) (Fig. 2). However, the lower bound in frequency of the spectral power is different between the reanalyses: JRA-55C, JRA-55, and MERRA-2 exhibit broader frequency spectra than the others. Above 70 hPa, another two local peaks appear in a higher-frequency region (k = −4, ω ∼ 0.3 cyc day −1 ; k = −3, ω ∼ 0.4 cyc day −1 ), in addition to the aforementioned peaks, 5 and these peaks become primary at 20 hPa. At 50 hPa, the spectral shapes are generally similar among the reanalyses, although the spectral region where the majority of the powers exist in CFSR is slightly shifted toward lower zonal wavenumbers when compared with the other reanalyses. Figure 2 also exhibits a distinct feature in the v S spectra at 20 hPa, compared to the spectra at the lower altitudes: statistically significant power appears along the narrow spectral region that includes k = −5, ω ∼ 0.5 cyc day −1 and extends into higher 10 wavenumbers and frequencies in all of the reanalyses. Toward lower wavenumbers, it seems to merge into the aforementioned peak at k = −3, ω ∼ 0.4 cyc day −1 . The spectral power along this region is larger in JRA-55C, JRA-55, and CFSR than in the others. This portion of the spectrum is further examined in Section 4. As will be seen therein, the waves with this spectrum do not originate from below, and they appear with different timing and different characteristics from the lower-frequency upward-propagating MRG waves which dominate the v S spectra in the lower altitudes. In Sections 3.2-3.4, we focus on the 15 lower-frequency MRG waves filtered with a cut-off frequency of 0.33 cyc day −1 (period of 3 days).
The spectral shapes of Kelvin and MRG waves obtained from the ML datasets (not shown) are very similar to those from the SL datasets (Figs. 1 and 2) for each reanalysis, whereas their spectral power is larger by up to about 35%, depending on the altitude and reanalysis. This will be further discussed in Section 3.3.
In the following sections, we define Kelvin waves as the symmetric mode with h = 8-240 m, k = 1-10, and periods of 2-20 20 days following Fujiwara et al. (2012), unless otherwise stated. These spectral components include a major portion of the Kelvin wave variances (Fig. 1) while excluding contributions of the other disturbances at low phase-speed ranges at 100 hPa (see also Fig. 9). The MRG waves are defined as the antisymmetric mode with h > 8 m, −10 ≤ k < 0, and periods of longer than 3 days, as previously mentioned, where h is for U = 0. The perturbations filtered for these spectral components are denoted as, for example, T Kelvin for the Kelvin wave temperature or v MRG for the MRG wave meridional wind.

Spatial structure
To investigate the spatial distributions, the Kelvin and MRG signals are reconstructed in physical space by filtering their spectral components. Figure 3a shows the distributions of T Kelvin and v MRG variances at 100, 70, 50, and 20 hPa during the period of 1981-2010, averaged for ERA-I, MERRA-2, CFSR, and JRA-55 (the reanalyses that assimilate satellite data are included as the ensemble members, but MERRA is excluded since it exhibits very similar spatial distributions to MERRA-2, as will be seen in with the equatorial wave theory. Note that, while the phase of MRG waves propagates westward, their wave packet travels eastward when the background flow is westerly or weak easterly which is a preferred condition for MRG wave propagation in the stratosphere. The maxima of variances for the Kelvin (MRG) waves are located in the eastern (western) hemisphere in the lower stratosphere: 70 and 50 hPa, consistent with previous observational studies (Alexander et al., 2008;Yang et al., 2012;Kiladis et al., 2016). At 20 hPa, the Kelvin and MRG waves show rather broad distributions in longitude. The distributions of the wave variances at 100 hPa are closely related to those of background zonal wind in the upper troposphere (Yang et al., 2012;Flannaghan and Fueglistaler, 2013). The easterly (westerly) upper tropospheric wind in the eastern (western) hemisphere 5 allows Kelvin (MRG) waves to more readily propagate vertically, resulting in the hemispheric difference in the wave variances.
In addition, the MRG wave variances at 100 hPa are distributed quite broadly in longitude (Fig. 3a). Yang et al. (2012) showed that during the northern winter the MRG wave variances in the upper troposphere are much larger in the western hemisphere with peaks at about 140 • W and 30 • W than in the eastern hemisphere, whereas during the northern summer, the variance distribution becomes broad, stretching to the western Pacific (see their Fig. 6). The locations of the wave variances vary with respect to the QBO winds in the stratosphere (not shown). For example, MRG wave variances averaged for the 50-hPa easterly QBO are located to the west of those for the westerly QBO. However, the MRG wave variances are much smaller in the easterly QBO phases, and their contribution to the climatological-mean variances 25 shown in Fig. 3 is relatively small. Therefore, the results in Fig. 3 are more representative of the westerly QBO phases for MRG waves and, similarly, easterly phases for Kelvin waves.
The currently used filters for the MRG waves exclude the low-frequency perturbations (< 0.1 cyc day −1 ) at 100 hPa revealed in Fig. 2, of which the spectrum follows the dispersion curves for U ∼ 10 m s −1 . An additional calculation, the same as with Fig. 3, is performed but for the 100-hPa low-frequency components of v S using the filters for 0.033 ≤ ω < 0.1 cyc day −1 and 30 −10 ≤ k < 0 (Fig. S2). It is observed in all the reanalyses that the low-frequency perturbations are located mostly in the western hemisphere where westerlies exist (over the eastern Pacific and Atlantic), consistent with the Doppler shifted dispersion curves discussed in Fig. 2.
To further investigate the spatial structures of the Kelvin and MRG waves, including the circulation patterns and horizontal scales of representative wave modes in each reanalysis, an EOF analysis is used following the technique outlined in 35 in the equatorial meridional wind centered at around the 4.5-day period (Fig. 2). Dynamical fields associated with each EOF were obtained by projecting unfiltered ERA-I data at each grid point onto the associated principal component (PC) time series.
Here we use data interpolated to 2.5 • resolution for the six reanalyses considered. A similar technique as in Kiladis et al. 5 (2016) is used to obtain the statistical structures of the equatorial Kelvin waves, except that a 2-25 day eastward-only filter band is applied to the equatorial zonal wind based on the spectral peaks in Fig. 1. The results are very robust to changes in the filtering, as long as the filter band contains the spectral peaks shown in Figs. 1 and 2, and to changes in latitudinal extent of the EOF basis.
In all cases, EOF "pairs" are obtained with respective PC time series that correlate at better than 0.96, which together 10 represent the propagating pattern of Kelvin or MRG waves. Thus, each mode can be represented by either EOF pattern and its associated PC. Figure (Table 2).
Lag regressions based on the PC time series show the eastward propagation of the Kelvin wave signals at a mean phase speed of around 35 m s −1 (not shown), which is rather faster than the 20-30 m s −1 found in previous studies, although these phase speeds are highly dependent on the state of the QBO (e.g., Randel and Wu, 2005;Ern et al., 2008;Lott et al., 2009Lott et al., , 2014. 25 While these features will be explored in further detail in a future study, the main point here is that the various reanalyses appear to be quite suitable for studying the mean statistical structure of these waves, as well as their variability, as will be demonstrated below. The situation for the MRG wave signals is rather different (Fig. 5). The leading MRG EOFs represent localized wave packets that are confined over the eastern Pacific to Atlantic sector, as expected from the location of the variance peaks in Fig. 3. 30 The EOF patterns have classical structures of equatorial gyres and off-equatorial antisymmetric geopotential perturbations as predicted by Matsuno (1966). The reanalyses each represent structurally very similar MRG wave patterns, except for CFSR which has a zonally broader signal (i.e., longer wavelengths), consistent with the spectra shown in Fig. 2. In addition, the MRG wave patterns in JRA-55, JRA-55C, and CFSR are displaced somewhat to the west of the other reanalyses, again reflecting their variance distributions shown in Fig. 3b. The higher mode EOF pairs of MRG waves (not shown) have wave structures over the Indian and Pacific Oceans, with good agreement among the reanalyses for up to the fourth mode pairs.
Kelvin and MRG signals at other levels from 70 to 10 hPa are broadly similar to those shown here, although the structures change considerably at 100 hPa (not shown). In summary, there is reasonable agreement between the reanalyses in their statistical representation of Kelvin and MRG wave structures, at least for the leading modes.

Long-term change and satellite effects
In this section, we analyze long-term variations in the Kelvin and MRG wave activity and discuss the impact of satellite data upon the assimilated wave activity based on the comparison of JRA-55 and JRA-55C. Figure 6 shows annual-mean time series of 100-hPa T Kelvin and v MRG variances, as defined in Section 3.1. The results from the ML datasets are presented (solid), except for MERRA for which the SL datasets are used (dashed) due to the lack of ML data. The variance of T Kelvin in 10 ERA-I fluctuates between about 0.5 and 0.85 K 2 with an increasing trend, whereas that in JRA-55C fluctuates between ∼0.35 and 0.58 K 2 without such an obvious trend. Despite this difference, much of the interannual variability in the Kelvin wave activity is reflected in both time series, with a Pearson linear correlation between the two of 0.85 (0.88 after trend removal).
Similar correlations are obtained between the other series, with the variances ranging between those in ERA-I and JRA-55C.
By comparing JRA-55 and JRA-55C, it appears that assimilation of satellite data leads to an increase in the Kelvin wave 15 temperature variance, as also found in the previous section (Fig. 3b).
The difference in the variance between JRA-55 and JRA-55C is further investigated in Fig. 7, which presents the time series of the annually-averaged difference at altitudes from 100 to 5 hPa. Because the wave variance in the stratosphere is strongly dependent on the QBO (e.g., Yang et al., 2012), 1-2-1 smoothing is applied to the annual time series to filter out the quasibiennial fluctuations and focus on the longer-term variations. It is evident that the difference (JRA-55 minus JRA-55C) in the 20 T Kelvin variance is generally positive from 100 to 20 hPa during the 30-year period, indicating the enhancement of the Kelvin wave amplitude by assimilation of satellite data in JRA-55. In addition, it is noteworthy that at 100 hPa, the difference is systematically larger from around 2000, compared to that before late 1990s: it is roughly 6% (up to ∼8%) in the years before 1998, whereas it increases to ∼10% in 1998-2000 and becomes ∼20% afterward. A similar systematic change of the difference in the T Kelvin variance is even more evident in the upper stratosphere: at 10 and 5 hPa ( functions of the two suggest that the AMSU-A instruments could be advantageous over the SSU, particularly at an altitude of ∼100 hPa which the latter does not cover. This may explain the systematic change found at 100 hPa around 1998 (Fig. 7). The higher vertical resolution of the AMSU-A instruments is expected to mostly benefit the upper stratosphere where radiosonde sounding observations do not reach (i.e., above ∼10 hPa). This is consistent with the observed systematic change being most prominent in the upper stratosphere. In addition to the TOVS-ATOVS transition, there is also a possibility that the Global Navigation Satellite System Radio Occultation (GNSS-RO) has influenced the assimilated waves in JRA-55 since 2006 (Fig. 10 in Fujiwara et al., 2017), although this cannot be identified by the current analysis. 5 Recalling that JRA-55C does not exhibit a long-term trend in Kelvin wave activity at 100 hPa (Fig. 6), the trend shown in JRA-55 is probably not an actual change in the true atmosphere, but instead an artifact arising due to the satellite transition.
The change in the 100-hPa T Kelvin variance around 1998 is 17% in JRA-55 when it is measured by the difference between the variances averaged for the two periods before and after 1998, i.e., 1981-1997 (P1 hereafter) and 1999-2010 (P2), relative to the 30-year mean variance. The changes in the other reanalyses are comparable to or smaller than that in JRA-55 (e.g., 19% 10 in ERA-I, which is the largest value). This suggests that the long-term changes in the 100-hPa Kelvin wave activity shown in those reanalyses also could result largely from the satellite transition.
The variance of v MRG exhibits generally similar interannual variations at 100 hPa between all of the reanalyses (Fig. 6), while the magnitudes of the variance are different. The v MRG variance has an increasing long-term trend from the early 1990s, even in JRA-55C. The v MRG variance at 100 hPa in JRA-55 is always larger than that in JRA-55C by roughly 6-12% (Figs. 6 15 and 7), reflecting the impact of satellite data on the assimilated MRG waves. In addition, the assimilation of satellite data increases the analyzed MRG wave activity at 5 hPa in JRA-55 (Fig. 7). This increase is up to ∼55%, and it becomes even larger at higher altitudes (see Fig. 8). Notable differences between P1 and P2 in the satellite impact on the assimilated MRG wave activity are not found below 5 hPa (Fig. 7). At 5 hPa, the impact is significantly larger in P2 than in P1 (note that two   The T Kelvin variance decreases with height from these levels to ∼30 hPa, and gradually increases from ∼20 hPa in both periods in all the reanalyses, except in MERRA-2 during P1. Further examination indicates that sometimes the temperature seems not to be properly assimilated in the upper stratosphere during P1 in MERRA-2 ( Fig. S3; note the lack of Kelvin wave peaks in 1984, 1989, and 1991), which may result in the decrease above ∼20 hPa in P1 in MERRA-2. It has previously been reported that the monthly mean zonal wind at 10 hPa in MERRA-2 shows significant differences from observations, in 30 particular with larger annual/semiannual oscillations, until the mid-1990s (Coy et al., 2016;Kawatani et al., 2016), which is perhaps partly related to the under-representation of the Kelvin wave temperature shown in Fig. S3. Errors in a mean state of assimilated fields can cause degradation of any variables describing waves during the first-guess model integration, which in turn degrade the assimilated fields 3 or 6 hours later. In P2, the T Kelvin variance at 70-30 hPa has similar magnitudes between ERA-I, MERRA-2, and CFSR, while JRA-55 exhibits a relatively smaller variance than these reanalyses, as is also evident in Fig. 1. In the upper stratosphere, the variance is notably larger in ERA-I in P2 than in the others (also see Fig. S3).
The difference in the T Kelvin variance between P1 and P2 is positive in the whole stratosphere in all of the reanalyses (Fig. 8).
In JRA-55C, the difference is small at 100 hPa (∼6%), as already discussed (Fig. 6), and above 10 hPa (4-8%), although it reaches ∼20% at 50-20 hPa. Given the small difference in Kelvin wave activity at ∼100 hPa, the increase in the activity at 5 50-20 hPa could result from interaction with the QBO which itself has variability in its morphology. However, since only five (seven) QBO cycles are included in P2 (P1), further investigation of the long-term change associated with the QBO might be less meaningful in a statistical sense. Another possible reason could be the increased number of radiosonde observations around 30 hPa with time (see Fig. 15 in Kawatani et al., 2016). In JRA-55, the difference is large in the upper stratosphere (∼25% or larger), and it is even larger in ERA-I, MERRA-2, and CFSR. Similarly to the case of JRA-55, the large increase in 10 Kelvin wave activity in the middle and upper stratosphere in the other reanalyses could also be in large part due to the transition between satellite instruments. From Fig. S3, it is evident that the difference in the Kelvin wave variances at 10 hPa between ERA-I and the others increases abruptly around 2000, which may support such an impact of the satellite instrument change in ERA-I.
In all the reanalyses, the v MRG variance decreases with height below ∼20 hPa and increases above ∼10 hPa (Fig. 8). In 15 both periods, MERRA-2 and CFSR exhibit larger variances than the others, particularly in the upper stratosphere. In P1, it is seen in the upper stratosphere that the v MRG variance in the ML result of CFSR varies less smoothly in the vertical and that the variance in MERRA-2 is about 2-3 times larger than that in the others. These features are still present to a lesser extent in P2.
Examination of time series of the v MRG variance (Fig. S3) shows that some unexpected peaks appear in P1 in MERRA-2 with exceptionally large magnitudes even during the easterly QBO phase, which are primarily responsible for the large value of the 20 v MRG variance in P1 shown in Fig. 8. The unexpected peaks disappear suddenly starting in 1999 (Fig. S3), likely emphasizing the impact of the AMSU-A observations on the assimilated wave activity in MERRA-2. Regarding the vertical fluctuations exhibited in CFSR, the reason for this is not obvious. The profiles of the difference between P1 and P2 may imply that above 10 hPa, the large increase in the v MRG variance in JRA-55 (10-25%) and ERA-I (∼40%), compared to the small increase in JRA-55C (∼5%), is probably in large part due to the satellite transition. 25 In each reanalysis, the variances obtained from the SL datasets are smaller than those from the ML datasets in both periods ( Fig. 8), by 3-29% (0-38%) for Kelvin (MRG) waves depending on the altitudes and reanalyses. The differences in the variances between the ML and SL results (ML minus SL) normalized by the variances in the ML results at each altitude are reported in Table 3. The smaller amplitudes of the waves in the SL datasets compared to those in the ML datasets result from vertical interpolation of reanalysis output variables, which inevitably damps wave perturbations in the standard-level 30 products of reanalyses. The damping by the interpolation is more significant for waves with smaller vertical wavelengths (e.g., Kim and Alexander, 2013;Kim and Chun, 2015a). In addition, for a given pressure level, the damping of waves in the SL datasets also depends on the distance between the given level and its adjacent model level. That is, the interpolation effect must be less if a model level is very close to the given pressure level. For example, the distance between 100 hPa and its adjacent model level in MERRA-2 is very small (< 50 m), and thus the difference between the SL and ML results of MERRA-2 is only 35 3-4% at 100 hPa (Table 3). The same (opposite) is true for the 10 hPa (50 hPa) in ERA-I, with a small (large) difference at that level. In JRA55 and JRA-55C, the difference in the v MRG variance between the ML and SL results is less than 5% at all levels, implying that the ML fields of these two reanalyses contain less of the MRG wave perturbations with small vertical wavelengths (∼2∆z-4∆z, where ∆z is the vertical grid spacings of their model), compared to the other reanalyses. 5 In this section, we compare the vertical EP flux (F z ) spectra and the wave forcing of the QBO calculated by the EP flux divergence among the reanalyses. As for the figures discussed in Section 3.1, we obtain k-ω spectra for each of the symmetric and antisymmetric modes but for F z averaged over 5 • N-5 • S (Figs. 9 and 10, respectively; see Section 2 for the narrowed latitude band). The spectral characteristics of F z for the symmetric and antisymmetric modes are qualitatively similar to those of T S and v S , respectively, for each reanalysis. For the symmetric mode, the F z spectra are aligned along the Kelvin wave 10 dispersion curves with h values that are similar to those for T S , along with exhibiting a gradual shift of the spectra to the higher h with increasing altitudes (Figs. 1 and 9). On the other hand, the symmetric F z and T S spectra differ in that the F z (T S ) spectra have peaks at k = 3-5 (2-3) in all of the reanalysis and that the magnitudes of F z in JRA-55 and JRA-55C are not smaller compared to those in the others as was the case for temperature (cf. Fig. 1). As already mentioned, JRA-55 and JRA-55C have relatively weak Kelvin wave temperature amplitudes but not wind amplitudes (Fig. S1), and the major term of F z for Kelvin 15 waves is the vertical momentum flux (w ′ u ′ ). The reason that the peak of F z occurs at larger zonal wavenumbers than in the T S spectra is because the vertical wind perturbations tend to have large amplitudes at higher zonal wavenumbers: the ratio of vertical to horizontal wind amplitudes is proportional to k/m (where m is the vertical wavenumber), as from the continuity equation. In the Kelvin wave dispersion relationship, k is independent of m, so that the ratio is proportional to k for given m.

EP flux
For the antisymmetric mode, F z has similar spectral shapes and roughly the same (k, ω) peaks as those of v S for each 20 reanalysis below 20 hPa (Figs. 2 and 10). Consistent with the v S spectra, the frequency bounds of the F z spectra at 70 hPa are quite a bit lower in JRA-55, JRA-55C, and MERRA-2 than in the other reanalyses. At 20 hPa, the F z spectra differ from the v S spectra in that the relatively low frequency and high wavenumber (and thus low phase speed) components are largely suppressed in the F z spectra. Indeed, these low-frequency components are more emphasized in the meridional EP flux (F φ ) spectra (Fig. S4), implying that the low-frequency MRG waves tend to propagate laterally rather than vertically. In addition, 25 even negative F z values appear in the 20-hPa spectra at k = 6-7 in MERRA-2, while this feature is less clear in the other reanalyses. This is suggestive of coexistence of upward and downward propagating waves at this level which could be another reason of the small F z in the low-frequency range.
The EP flux spectra are reconstructed as a function of phase speed to investigate the vertical propagation and dissipation of Kelvin and MRG waves during different QBO phases and to quantify the QBO forcing by those waves as well as the phase-30 speed ranges responsible for the forcing. For this, we construct bins of phase speeds (c j ) with a width of 2 m s −1 (∆c) and integrate the k-ω spectral densities of F φ and F z across the corresponding bins (i.e., (2πa)ω/k ∈ [c j −∆c/2, c j +∆c/2], where a is the mean radius of the Earth). 1 The EP flux divergence is also calculated in each bin. For each wave type, the same filters as described in Section 3.1 are used: h = 8-240 m, k = 1-10, and 0.05 < ω ≤ 0.5 cyc day −1 for Kelvin waves and h > 8 m, |k| = 1-10, and ω < 0.33 cyc day −1 for MRG waves. Figure 11 shows vertical profiles of 5 • N-5 • S averaged EP flux divergence and F z as a function of phase speed, for Kelvin and MRG waves at c > 0 and c < 0, respectively, composited for four different phases of the QBO. The composite is made by 5 selecting one month for each QBO cycle when the zonal wind tendency is largest at a given altitude and then averaging over the 13 cycles available. The observed monthly-mean near-equator wind profiles compiled by the Freie Universität Berlin (FUB) are used to select the months. The four composites presented in Fig. 11 are for the months of maximum westerly and easterly tendencies at 20 and 50 hPa. The composite averages of the FUB wind are also plotted in Fig. 11 (thick green lines) as well as those of each reanalysis wind (thin red lines). The results are generated using the ML datasets of the four reanalyses available 10 for EP flux calculation (ERA-I, MERRA-2, JRA-55, and JRA-55C), since the EP flux could be largely underestimated when using the SL datasets. Furthermore, the height dependence of the amplitude damping, as reported in Table 3, could also affect the estimation of the vertical divergence of EP flux (see Fig. 3 in Kim and Chun, 2015a, for a large difference in the EP flux divergence between ML and SL datasets of ERA-I).
The EP flux divergence of the Kelvin waves is found to be quite different between reanalyses not only in magnitude but also 15 in its peak altitude. At the maximum westerly tendency phase at 20 hPa (Fig. 11, the first row), the vertical EP flux is notably the largest in ERA-I at ∼20 hPa. However, the EP flux divergence at this level is the largest in MERRA-2, implying that the assimilated wave fields exhibit the most severe dissipation in MERRA-2 around this level. ERA-I and JRA-55 show the peaks of the EP flux divergence at ∼15 and 12 hPa, respectively. It is important to note that the easterly-to-westerly transition (i.e., zero wind) occurs at different altitudes between the reanalyses (Fig. 11, red lines; also refer to Fig. 16 in Kawatani et al., 2016, 20 for different time evolutions of the phase transition) and that the maximum Kelvin wave forcing is found to occur at around the zero-wind altitudes in each reanalysis. These suggest that the mean wind in the assimilated field could affect the representation of wave behaviors in reanalyses in the shear zone. The phase speed at which the forcing occurs increases with altitudes from ∼10-30 m s −1 at 20 hPa to 15-40 m s −1 at 10 hPa but with large spread at 10 hPa among the reanalyses. The Kelvin wave forcing integrated for the phase speed has maxima of 9.2, 10.2, and 10.0 m s −1 month −1 in ERA-I, MERRA-2, and JRA-55, 25 respectively. When the westerly tendency is maximized at 50 hPa (Fig. 11, the second row), the peak of the Kelvin wave forcing appears slightly below (above) 40 hPa in ERA-I and MERRA-2 (JRA-55). Similar to the 20-hPa transition phase, the vertical EP flux is the largest in ERA-I, while the EP flux divergence at ∼40 hPa is larger in MERRA-2 (8.6 m s −1 month −1 ), which may also be related to the different mean-wind profiles in the reanalyses. Note that the mean winds in ERA-I and JRA-55 have weaker shear at 40-50 hPa with weaker westerly magnitudes above 40 hPa compared to the observed winds (Fig. 11,   30 thick green lines), and they might affect the wave dissipation rate in these reanalyses. At the other two phases of the QBO (Fig. 11, the third and last rows), the zonal wind is westerly at 70 hPa, causing less upward propagation of Kelvin waves in the lowermost stratosphere. Figure 12 shows the vertical divergence of the EP flux, (ρ 0 cos φ) −1 ∂F z /∂z, as a function of phase speed, along with F z as in Fig. 11. In the main regions of the Kelvin wave forcing during the westerly-shear phases (Fig. 12, the first and second rows), the vertical divergence of the EP flux for the Kelvin waves is roughly 50-75% of the total divergence shown in Fig. 11. 5 The vertical divergence is the largest in ERA-I among the reanalyses, whereas the total divergence of the EP flux is larger in MERRA-2 (Fig. 11). This indicates that the contribution of meridional flux divergence to the Kelvin wave forcing is relatively large in MERRA-2 (40-50%, not shown) compared to that in the other reanalyses (25-45%). The Kelvin wave flux tends to be slanted toward the equator in westerly-shear winds (e.g., Plumb and Bell, 1982;Kim and Chun, 2015b). The larger wind shear in MERRA-2 is likely the cause for the larger contribution of meridional divergence of the flux. 10 MRG waves dissipate largely in the lower stratosphere when the zonal wind is easterly at 70 hPa (Fig. 11, the first and second rows). The MRG wave forcing is about 2-3 m s −1 month −1 below 70 hPa in these QBO phases in the reanalyses.
When the easterly tendency is maximum at 20 hPa (Fig. 11, the third row), MRG waves propagate through the lower stratospheric westerlies and dissipate at around 30-20 hPa. The forcing by the MRG waves at these altitudes is largest in JRA-55 (∼1.6 m s −1 month −1 ), where it occurs at phase speeds between −10 and −30 m s −1 . It is also found in Fig. 11 (the first 15 and last rows) that the EP flux divergence of MRG waves exists in the westerly-shear zone at 10-20 hPa in MERRA-2 when the easterly jet is located at 20-40 hPa. This signal also exists in the other reanalyses but with much weaker magnitudes. The waves responsible for this signal seem not to originate from below, as the vertical EP flux from the lower altitudes is interrupted by the easterly winds. The phase speeds of the signal (around −20 m s −1 ) correspond to the spectral region where F z is suppressed in Fig. 10 (but F φ is emphasized in Fig. S4) mentioned previously. The EP flux divergence in this region is about 3-20 4 m s −1 month −1 in MERRA-2, although it is uncertain to estimate because of the large discrepancy between MERRA-2 and the others and the larger westerly shear in MERRA-2 compared to the observed (Fig. 11). Previously, Maury and Lott (2014) also found that MRG waves with the phase speed of −19 m s −1 appear at 20 hPa when there exist easterlies (westerlies) at 50 hPa (20 hPa), using ERA-I. Recently, using a climate model, Garcia and Richter (2019) simulated MRG waves existing in a westerly-shear zone (see their Fig. 9b), which were generated locally by barotropic instability of the westerly QBO jet. These 25 waves have a spectral peak at around k = 7-8 and ω ∼ 0.2 cyc day −1 (see their Fig. 8a) which is similar to that detected in MERRA-2 (Fig. S4). The simulated MRG waves induced westward forcing of 7.5-15 m s −1 month −1 .
In Fig. 12, it is seen that for MRG waves, the contribution of vertical divergence of the EP flux to the total divergence is very small everywhere (∼10%). In addition, the EP flux divergence and F z for each of the composited 13 QBO cycles, as well as the individual profiles of the mean winds, are included in the Supplement (Fig. S5). While there are substantial cycle-to-cycle 30 variations, the overall structures of the EP flux and its divergence for all cases are generally similar to those of the composite averages shown in Fig. 11.
As mentioned in the Introduction, the momentum budget of the QBO is currently not fully constrained by observations, and it is believed to be essential for modeling the QBO. The different altitudes of Kelvin wave forcing around 15 hPa among the reanalyses, as well as the different magnitudes, shown in Fig. 11 may imply a limitation in quantifying wave forcing 35 of the QBO using a reanalysis dataset. Advances in observation and assimilation of stratospheric waves might be needed to reduce the spread of the assimilated waves and to further improve global models. In addition, given that the representation of equatorial waves in reanalyses is dependent on the assimilated mean winds (Fig. 11), reality of small-scale gravity-wave parameterization in forecast models, which can affect the details of the QBO winds, might also be important for the assessment of the large-scale wave forcing. The rather weak magnitudes of the MRG wave forcing in the easterly-shear zone in all of the 5 reanalyses imply that other waves (especially gravity waves) might play a more important role in driving the easterly phase of the QBO than MRG waves (Kawatani et al., 2010;Evan et al., 2012;Ern et al., 2014) and/or that current observational measurements, as well as models and assimilation methods, cannot fully capture the MRG waves which have relatively small vertical wavelengths. Modeling studies have demonstrated that adequate representation of the stratospheric equatorial waves requires vertical resolutions of ∼700 m or finer in the lower to middle stratosphere (Giorgetta et al., 2006;Richter et al., 2014), 10 while the vertical spacings of the prediction models used for the current reanalyses are about a factor of two larger than that (Table 1). It is noteworthy that the aforementioned simulation of locally generated MRG waves by Garcia and Richter (2019) used a 500-m vertical resolution, in which the MRG waves induced large momentum forcing (∼10 m s −1 month −1 ), although they did not act to drive the QBO but to neutralize the westerly jet.

15
In the k-ω spectrum for the symmetric component of meridional wind (Fig. 2), a distinct spectral peak is seen at 20 hPa, which stretches toward high frequencies and high wavenumbers from around ω ∼ 0.4 cyc day −1 and k ∼ −3. To investigate this further, we filter a part of the spectrum that does not coincide with the spectrum of the lower-frequency MRG waves that are investigated in Section 3. The filtered spectral region is 0.5 < ω < 0.7 cyc day −1 for k = −7 and 0.6 < ω < 0.75 cyc day −1 for k = −8. Figure 13 presents time series of variances of the filtered meridional-wind spectrum at 20 and 10 hPa, along with 20 the observed near-equator zonal wind at each altitude compiled by FUB. For all reanalyses, it is clearly seen that this spectral peak appears only when the background wind is easterly with substantial speeds at each altitude. Furthermore, there exist time lags of the variance peaks between the two altitudes, as well as a lag in the zonal wind due to the downward progression of the QBO. Due to these time lags, there exist periods during which the variance is large at 10 hPa while it is much smaller at 20 hPa at the same time (e.g., late 1999/early 2000 and late 2008/early 2009). This indicates that the spectral signal of these 25 waves does not originate from below. On the other hand, there also exist periods during which the variance is much larger at 20 hPa than at 10 hPa (e.g., late 2007), suggesting that the waves at 20 hPa are, at least in part, generated in situ. In fact, the vertical EP flux of this spectrum is not well identified (see Fig. 10), which could indicate a minimal preference for upward or downward propagation of the waves. This also could be suggestive of in-situ wave generation. This spectral peak is also seen in the upper stratosphere above 10 hPa, although it is less evident here because the peak is not clearly separated from that of 30 the lower-frequency waves (not shown).
Given the dependence of that spectrum filtered for k = −7 and −8 upon the background wind U (Fig. 13), the entire spectrum of those waves can be more clearly identified by compositing for the periods with strong background easterlies. Figure 14 presents two composites of the symmetric meridional-wind spectrum for the months with U < −25 m s −1 and U > −20 m s −1 at 20 hPa, based on the FUB zonal wind. The composite result is shown only for JRA-55C (in which the amplitude of this spectrum is the largest in Fig. 2) but is qualitatively similar for the other reanalyses. For U < −25 m s −1 , it is shown that the spectrum has peaks at k = −4 and −5 with periods of about 2 days. The ground-based zonal phase speeds corresponding to the peaks are between −45 and −60 m s −1 . These are much higher compared to those of the waves detected 5 in the westerly QBO wind in Maury and Lott (2014) and Garcia and Richter (2019) discussed in Section 3.4 (about −20 m s −1 ; see also Fig. 11, the first and last rows). The spectral region is best identified by the dispersion curves of MRG waves (Fig. 14) rather than the other types of waves, provided that U ∼ −30 m s −1 . The spectrum composited for U > −20 m s −1 is similar to that of the lower-frequency MRG waves observed in the lower stratosphere. The variances of the symmetric meridional wind by the waves observed at 20-10 hPa during the easterly QBO phase typically have magnitudes of roughly half of those by the 10 lower-frequency MRG waves during the westerly phase (not shown).
The high-frequency waves detected here have similar characteristics to the quasi-two-day wave (QTDW) observed in the summer mesosphere and uppermost stratosphere: the latter also occurs in strong easterly winds with a ∼2-day period. It has been proposed that the QTDW is a MRG mode triggered by the baroclinic and/or barotropic instabilities associated with the easterly jet in the summer mesosphere (e.g., Salby, 1981;Pfister, 1985;McCormack et al., 2009). However, given that the 15 occurrence of the 20-hPa disturbances is strongly tied to certain phases of the QBO, as seen in Figs. 13 and 14, they might not be necessarily connected from the QTDW above. More details on the structure of these waves and their possible impacts in the tropical stratosphere will be pursued in a future study.

Summary
The equatorial Kelvin and MRG waves in the TTL and stratosphere represented in six reanalyses (ERA-I, MERRA, MERRA- The spatial distributions and patterns of the equatorial waves in the lower stratosphere are investigated . It is shown in all of the reanalyses that Kelvin and MRG wave variances tend to be large in the eastern and western hemispheres, respectively, and that the locations of their maxima tilt eastward in the vertical. However, the longitudinal variations of the Kelvin wave variances are much smaller in CFSR compared to those in the others. The leading mode EOFs of the Kelvin- 30 and MRG-wave filtered perturbations are used to obtain statistically representative patterns of the waves in each reanalysis, following Kiladis et al. (2016). For both waves, the horizontal winds and geopotential perturbations projected onto the leading EOFs show spatial patterns and phase relationships between variables that are consistent with the classical equatorial wave theory. All the reanalyses studied here exhibit remarkably similar patterns for the Kelvin wave leading modes, which have zonal wavenumber 1 structures with larger amplitudes in the eastern hemisphere. The MRG wave leading modes are confined over the eastern Pacific to Atlantic sector and also show reasonable agreement between the reanalyses, although their zonal scales are somewhat larger in CFSR.
From analysis of the time series of Kelvin and MRG wave variances, systematic differences are found between the periods 5 before and after the late 1990s in several aspects (Figs. 6-8).
-The difference in the Kelvin wave variances between JRA-55 and JRA-55C (which stems from the exclusion of satellite data in the assimilation in JRA-55C) shows significant changes after the late 1990s at 100 hPa (from ∼7% to 10-24%) and at 10 and 5 hPa (from between 0 and −10% to 10-20%) but not at ∼50 hPa.
-The Kelvin wave variances in the middle and upper stratosphere in ERA-I exhibit a large increase around 1998 (∼50%), 10 becoming significantly larger than those in the other reanalyses (Figs. 8 and S3).
-In the middle and upper stratosphere in MERRA-2, the peaks of the Kelvin wave variances during westerly-shear phases of the QBO are not represented in some years before the mid-1990s (Fig. S3). The MRG wave variances at 10 hPa show exceptionally large values during the westerly-shear phases until 1998, but not afterward.
-The MRG wave variances in CFSR fluctuate largely in the vertical from ∼30 hPa before 1998, but not afterward.

15
The results listed here demonstrate significant impacts of the TOVS-ATOVS transition starting in 1998 on the assimilated wave amplitudes in the middle and upper stratosphere in all the reanalyses. Below 10 hPa, the satellite impacts on the waves are identifiable in JRA-55 only by direct comparison with JRA-55C at 100 hPa where the SSU instruments in the TOVS suite do not cover.
The time series of the Kelvin wave variances at 100 hPa (Fig. 6) exhibit increasing trends in all the reanalyses except JRA-20 55C. However, the veracity of the trends in Kelvin wave activity is uncertain due to potential satellite effects in the reanalyses.
On the other hand, the MRG wave variances show a long-term increase at 100 hPa, common to all of the reanalyses including JRA-55C which is not affected by satellite transitions.
It is also noted that variances of equatorial waves can be underestimated by up to around 30% when standard-pressure-level datasets are used. The wave perturbations in these datasets are damped out by vertical interpolation, and the damping effect is 25 large for waves with small vertical wavelengths.
The EP flux and its divergence are presented as a function of phase speed and height to compare the QBO forcing measured in the reanalyses (Figs. 11 and 12). In general, the phase speeds at which the Kelvin wave forcing occurs are comparable among ERA-I, MERRA-2, JRA-55, and JRA-55C. The Kelvin wave EP flux is the largest in ERA-I, while its forcing is larger in MERRA-2 with more severe dissipation of the waves. For the QBO phases with westerly acceleration at ∼15 hPa, the height 30 of the peak forcing differs between the reanalyses by up to ∼3 km. The differences in the dissipation rate and peak forcing altitude are found to be related to the different mean-wind profiles in the reanalyses. The MRG wave forcing is in general small in all the reanalyses (up to ∼3 m s −1 month −1 for the 5 • N-5 • S average).
In addition, relatively high-frequency waves detected in the spectra of symmetric meridional winds at 20 hPa (Fig. 2), which are distinguishable from the MRG waves observed below, are further investigated (Figs. 13 and 14). It is found in all the reanalyses that they appear only when the background wind is easterly, unlike the MRG waves propagating from below, with ground-based periods of about 1-3 days. Their spectral shape is best identified by the dispersion curve of MRG waves with the background wind of about −30 m s −1 . More details regarding the spatial pattern, origin, and implication of these waves 5 remain to be studied in the future.
Author contributions. The EOF analysis in Section 3.2 was carried out by GNK, JRA, and JD, and the other calculations were by YHK.

Abbreviation Full name
Vertical resolution         and (right) their differences (1999-2010 minus 1981-1997) relative to the 30-year average (1981-2010). The dashed and solid indicate the results from the standard-level and model-level datasets, respectively.  those composites are also indicated using each reanalysis (thin red) and radiosonde observations (thick green; see the text for the composite method). The contour intervals of the EP flux are 0.008 and 0.0003 mPa/(m s −1 ) for the Kelvin and MRG waves, respectively, and every third contour is distinguished by thicker lines. Note that the contour interval of the EP flux divergence for the Kelvin waves is 4 times larger than that for the MRG waves.  0.5 < ω < 0.7 cyc day −1 with k = −7 and 0.6 < ω < 0.75 cyc day −1 with k = −8 (see Fig. 2). The monthly zonal wind from radiosonde observations is also presented at each level (dashed black).