1. Introduction
Equatorial waves are atmospheric waves whose amplitude is maximized at the equator and decreases exponentially as the distance of the waves from the equator increases. Their dispersion relation and spatial structure were studied theoretically by Matsuno (1966). He derived eigenmodes of equatorial waves by using a plane-wave assumption in the time and zonal directions on the linear shallow-water equation. The equatorial waves that were discovered by Wallace and Kousky (1968) and Yanai and Maruyama (1966) from radiosonde observation were identified as Kelvin waves and Rossby–gravity waves, respectively.
On the other hand, quasi-biennial and semiannual oscillations (QBO and SAO) of the zonal wind exist in the equatorial stratosphere. Previous studies have shown that these oscillations are driven by atmospheric waves. The relation between the waves and the zonal-mean zonal wind can be diagnosed by the transformed Eulerian mean (TEM) equations that were derived by Andrews and McIntyre (1976, 1978). The residual mean flow is expressed as the sum of the Eulerian mean flow and Stokes drift under the small-amplitude assumption and is approximately equal to the zonal-mean Lagrangian mean flow when the wave is linear, steady, and adiabatic and when no dissipation occurs. The Eliassen–Palm (EP) flux is equal to the product of the group velocity and the wave activity density under the Wentzel–Kramers–Brillouin (WKB) approximation and is a useful physical quantity for describing the wave propagation (Edmon et al. 1980). The residual mean flow and the zonal-mean zonal wind acceleration are related to the divergence of EP flux in the zonal momentum equation. When there are no critical levels, the divergence of the EP flux is zero for linear, steady, and conservative waves. Under such conditions, the waves neither drive the residual mean flow nor accelerate the zonal-mean zonal wind. This is known as the nonacceleration theorem (Eliassen and Palm 1961; Charney and Drazin 1961). In studies using the TEM equations and equatorial wave theory, it was shown that the QBO is mainly driven by gravity waves, equatorial Kelvin waves, and Rossby–gravity waves (e.g., Sato and Dunkerton 1997; Haynes 1998; Baldwin et al. 2001), and it is recognized that the SAO is mainly driven by gravity waves, equatorial Kelvin waves, and extratropical Rossby waves (e.g., Hirota 1980; Holton and Wehrbein 1980; Hitchman and Leovy 1988; Sassi and Garcia 1997; Antonita et al. 2007).
Kawatani et al. (2010) recently investigated the zonal variation of wave forcing associated with the equatorial Kelvin waves and inertia–gravity waves using the 3D wave activity flux derived by Miyahara (2006). They showed that this variation in the stratosphere results from zonal variation of the wave sources and from the vertically sheared zonal winds associated with the Walker circulation, depending on the phase of the QBO.
However, it is not clear whether this 3D wave activity flux can describe the equatorial Kelvin waves correctly because this flux is only applicable to inertia–gravity waves, not to equatorial waves. Although Kinoshita and Sato (2013a,b) newly formulated the 3D wave activity flux on the primitive equations, their formulas are not applicable in the equator region. The present study formulates the 3D residual mean flow and wave activity flux applicable to equatorial waves and shows that these formulas are applicable to the equatorial Kelvin waves.
The paper is arranged as follows. In section 2, the 3D Stokes drift for equatorial waves (EQSD) is derived from its definition for the equatorial beta-plane equations. The 3D wave activity flux for equatorial waves (3D-EQW-flux) is formulated by using the EQSD in section 3. It is shown that the 3D-EQW-flux divergence corresponds to the equatorial wave forcing to the mean flow. It is also shown that the meridional integral of the 3D-EQW-flux accords with a product of the group velocity and the meridional integral of the wave activity density for equatorial waves. Moreover, we investigate the 3D wave energy equation for equatorial waves. A summary and concluding remarks are given in section 4.
2. The time-mean 3D Stokes drift applicable to equatorial waves




In the next section, the EQSD for equatorial Kelvin waves, Rossby–gravity waves, and other types of equatorial waves are formulated from the definition in (2.4).
a. The 3D Stokes drift for equatorial Kelvin waves














b. The 3D Stokes drift for other types of equatorial waves













3. A formulation of the 3D wave activity flux for equatorial waves
a. The 3D residual mean flow and wave activity flux


b. The relation between 3D wave activity flux and group velocity
In the 2D TEM equation system, the meridionally integrated EP flux is equal to a product of the vertical group velocity and the meridionally integrated wave activity density (Andrews et al. 1987). It can be shown that the vertical component of 3D wave activity flux [(3.3c)] satisfies this relation, as in the following.







c. The wave energy equation for equatorial waves
This section examines how the 3D wave activity flux for equatorial waves is related to the wave activity density after the wave energy equation is derived. In this derivation, it is assumed that the time-mean wind shear is negligible.




4. Concluding remarks
In this study, the 3D Stokes drift is formulated from its definition for equatorial beta-plane equations (EQSD) when the slowly varying background field and small-amplitude perturbations are assumed. EQSD is applicable to all equatorial waves. The 3D wave activity flux (3D-EQW-flux) is formulated by substituting the EQSD into the time-mean zonal momentum equation. These expressions are derived using the time mean and are phase independent.
Next, it is shown that the latitudinal integral of 3D-EQW-flux accords with a product of the group velocity and the latitudinally integrated wave activity density in both zonal and vertical directions. This is an extension of the relation for the Eliassen–Palm flux on the 2D TEM equations. The present study also derives the 3D wave energy equation for equatorial waves.
As it is shown that
The EQSD and 3D-EQW-flux are partly different from the 3D Stokes drift and wave activity flux that are applicable to both gravity waves and Rossby waves (Kinoshita and Sato 2013a,b). The difference is due to assumptions of waves. Kinoshita and Sato (2013a,b) assume waves having meridional wavenumbers, and the term
Acknowledgments
We thank Hisashi Nakamura, Toshiyuki Hibiya, Masaaki Takahashi, Keita Iga, and Koutarou Takaya for their helpful comments and important suggestions. We also thank Rolando R. Garcia, Saburo Miyahara, and an anonymous reviewer for providing constructive comments. This study is supported by Grand-Aid for Research Fellow (22-7125) of the JSPS and by Grant-in-Aid for Scientific Research (A) 25247075 of the Ministry of Education, Culture, Sports and Technology, Japan.
APPENDIX A
APPENDIX B
Meridional Integral of 3D Wave Activity Flux Applicable to Inertia–Gravity Waves


APPENDIX C
Another Expression of 

In this section, we introduce another expression of

REFERENCES
Andrews, D. G., and M. E. McIntyre, 1976: Planetary waves in horizontal and vertical shear: The generalized Eliassen–Palm relation and mean zonal acceleration. J. Atmos. Sci., 33, 2031–2048, doi:10.1175/1520-0469(1976)033<2031:PWIHAV>2.0.CO;2.
Andrews, D. G., and M. E. McIntyre, 1978: Generalized Eliassen-Palm and Charney-Drazin theorems for waves on axisymmetric mean flows in compressible atmospheres. J. Atmos. Sci., 35, 175–185.
Andrews, D. G., J. R. Holton, and C. B. Leovy, 1987: Middle Atmosphere Dynamics.Academic Press, 489 pp.
Antonita, T. M., G. Ramkumar, K. K. Kumar, K. S. Appu, and K. V. S. Nambhoodiri, 2007: A quantitative study on the role of gravity waves in driving the tropical stratospheric semiannual oscillation. J. Geophys. Res.,112, D12115, doi:10.1029/2006JD008250.
Baldwin, M. P., and Coauthors, 2001: The quasi-biennial oscillation. Rev. Geophys., 39, 179–229, doi:10.1029/1999RG000073.
Charney, J. G., and P. G. Drazin, 1961: Propagation of planetary-scale disturbances from lower into upper atmosphere. J. Geophys. Res., 66, 83–109, doi:10.1029/JZ066i001p00083.
Edmon, H. J., B. J. Hoskins, and M. E. McIntyre, 1980: Eliassen-Palm cross sections for the troposphere. J. Atmos. Sci., 37, 2600–2616, doi:10.1175/1520-0469(1980)037<2600:EPCSFT>2.0.CO;2.
Eliassen, A., and E. Palm, 1961: On the transfer of energy in stationary mountain waves. Geofys. Publ., 22 (3), 1–23.
Haynes, P. H., 1998: The latitudinal structure of the quasi-biennial oscillation. Quart. J. Roy. Meteor. Soc., 124, 2645–2670, doi:10.1002/qj.49712455206.
Hirota, I., 1980: Observational evidence of the semiannual oscillation in the tropical middle atmosphere—A review. Pure Appl. Geophys.,118, 217–238, doi:10.1007/BF01586452.
Hitchman, M. H., and C. B. Leovy, 1988: Estimation of the Kelvin wave contribution to the semiannual oscillation. J. Atmos. Sci., 45, 1462–1475, doi:10.1175/1520-0469(1988)045<1462:EOTKWC>2.0.CO;2.
Holton, J. R., and W. M. Wehrbein, 1980: A numerical model of the zonal mean circulation of the middle atmosphere. Pure Appl. Geophys., 118, 284–306, doi:10.1007/BF01586455.
Kawatani, Y., K. Sato, T. J. Dunkerton, S. Watanabe, S. Miyahara, and M. Takahashi, 2010: The roles of equatorial trapped waves and internal inertia–gravity waves in driving the quasi-biennial oscillation. Part II: Three-dimensional distribution of wave forcing. J. Atmos. Sci.,67, 981–997, doi:10.1175/2009JAS3223.1.
Kinoshita, T., and K. Sato, 2013a: A formulation of three-dimensional residual mean flow applicable to both inertia–gravity waves and to Rossby waves. J. Atmos. Sci., 70, 1577–1602, doi:10.1175/JAS-D-12-0137.1.
Kinoshita, T., and K. Sato, 2013b: A formulation of unified three-dimensional wave activity flux of inertia–gravity waves and Rossby waves. J. Atmos. Sci., 70, 1603–1615, doi:10.1175/JAS-D-12-0138.1.
Kinoshita, T., Y. Tomikawa, and K. Sato, 2010: On the three-dimensional residual mean circulation and wave activity flux of the primitive equations. J. Meteor. Soc. Japan, 88, 373–394, doi:10.2151/jmsj.2010-307.
Matsuno, T., 1966: Quasi-geostrophic motions in the equatorial area. J. Meteor. Soc. Japan, 44, 25–43.
Miyahara, S., 2006: A three dimensional wave activity flux applicable to inertio-gravity waves. SOLA, 2, 108–111, doi:10.2151/sola.2006-028.
Sassi, F., and R. R. Garcia, 1997: The role of equatorial waves forced by convection in the tropical semiannual oscillation. J. Atmos. Sci., 54, 1925–1942, doi:10.1175/1520-0469(1997)054<1925:TROEWF>2.0.CO;2.
Sato, K., and T. J. Dunkerton, 1997: Estimates of momentum flux associated with equatorial Kelvin and gravity waves. J. Geophys. Res., 102, 26 247–26 261, doi:10.1029/96JD02514.
Wallace, J. M., and V. E. Kousky, 1968: Observation evidence of Kelvin waves in the tropical stratosphere. J. Atmos. Sci., 25, 900–907, doi:10.1175/1520-0469(1968)025<0900:OEOKWI>2.0.CO;2.
Yanai, M., and T. Maruyama, 1966: Stratospheric wave disturbances propagating over the equatorial Pacific. J. Meteor. Soc. Japan, 44, 291–294.